Recent advances in analog and digital quantum-simulation platforms have enabled exploration of the spectrum of entanglement Hamiltonians via variational algorithms. In this work we analyze the convergence properties of the variationally obtained solutions and compare them to numerically exact calculations in quantum critical systems. We demonstrate that interpreting the cost functional as an integral permits the deployment of iterative quadrature schemes, thereby reducing the required number of measurements by more than an order of magnitude even in the presence of noise. We further show that a modified ansatz captures deviations from the Bisognano-Wichmann form in lattice models, improves convergence, improves trainability and provides a cost-function-level diagnostic for quantum phase transitions. Finally, we establish that a low cost value does not by itself guarantee convergence in trace distance. Nevertheless, it faithfully reproduces degeneracies and spectral gaps, which are essential for applications to topological phases.
Recent advances in analog and digital quantum-simulation platforms have enabled exploration of the spectrum of entanglement Hamiltonians via variational algorithms. In this work we analyze the convergence properties of the variationally obtained solutions and compare them to numerically exact calculations in quantum critical systems. We demonstrate that interpreting the cost functional as an integral permits the deployment of iterative quadrature schemes, thereby reducing the required number of measurements by more than an order of magnitude even in the presence of noise. We further show that a modified ansatz captures deviations from the Bisognano-Wichmann form in lattice models, improves convergence, improves trainability and provides a cost-function-level diagnostic for quantum phase transitions. Finally, we establish that a low cost value does not by itself guarantee convergence in trace distance. Nevertheless, it faithfully reproduces degeneracies and spectral gaps, which are essential for applications to topological phases.
As a powerful theoretical construct, the entanglement Hamiltonian (EH) encapsulates the essential entanglement properties of a quantum many-body system. From the EH, one can extract a variety of entanglement quantities, such as entanglement entropies, negativity, and the entanglement spectrum. However, its general analytical form remains largely unknown. While the Bisognano-Wichmann theorem gives an exact EH form for Lorentz-invariant field theories, its validity on lattice systems is limited, especially when Lorentz invariance is absent. In this work, we propose a general scheme based on the lattice-Bisognano-Wichmann (LBW) ansatz and multi-replica-trick quantum Monte Carlo methods to numerically reconstruct the entanglement Hamiltonian in two-dimensional systems and systematically explore its applicability to systems without translational invariance, going beyond the original scope of the primordial Bisognano-Wichmann theorem. Various quantum phases--including gapped and gapless phases, critical points, and phases with either discrete or continuous symmetry breaking--are investigated, demonstrating the versatility of our method in reconstructing entanglement Hamiltonians. Furthermore, we find that when the entanglement boundary of a system is ordinary (i.e., free from surface anomalies), the LBW ansatz provides an accurate approximation well beyond Lorentz-invariant cases. Our work thus establishes a general framework for investigating the analytical structure of entanglement in the complex quantum many-body systems.
Learning the structure of the entanglement Hamiltonian (EH) is central to characterizing quantum many-body states in analog quantum simulation. We describe a protocol where spatial deformations of the many-body Hamiltonian, physically realized on the quantum device, serve as an efficient variational ansatz for a local EH. Optimal variational parameters are determined in a feedback loop, involving quench dynamics with the deformed Hamiltonian as a quantum processing step, and classical optimization. We simulate the protocol for the ground state of Fermi-Hubbard models in quasi-1D geometries, finding excellent agreement of the EH with Bisognano-Wichmann predictions. Subsequent on-device spectroscopy enables a direct measurement of the entanglement spectrum, which we illustrate for a Fermi Hubbard model in a topological phase.
The modular (or entanglement) Hamiltonian correspondent to the half-space bipartition of a quantum state uniquely characterizes its entanglement properties. However, in the context of lattice models, its explicit form is analytically known only for the two spin chains and certain free theories in one dimension. In this work, we provide a thorough investigation of entanglement Hamiltonians in lattice models obtained via the Bisognano-Wichmann theorem, which provides an explicit functional form for the entanglement Hamiltonian itself in quantum field theory. Our study encompasses a variety of one- and two-dimensional models, supporting diverse quantum phases and critical points, and, most importantly, scanning several universality classes, including Ising, Potts, and Luttinger liquids. We carry out extensive numerical simulations based on the density matrix renormalization group method, exact diagonalization, and quantum Monte Carlo. In particular, we compare the exact entanglement properties and correlation functions to those obtained applying the Bisognano-Wichmann theorem on the lattice. We carry out this comparison on both the eigenvalues and eigenvectors of the entanglement Hamiltonian, and expectation values of correlation functions and order parameters. Our results evidence that as long as the low-energy description of the lattice model is well captured by a Lorentz-invariant quantum field theory, the Bisognano-Wichmann theorem provides a qualitatively and quantitatively accurate description of the lattice entanglement Hamiltonian. The resulting framework paves the way to direct studies of entanglement properties utilizing well-established statistical mechanics methods and experiments.
Quantum variational algorithms are one of the most promising applications of near-term quantum computers; however, recent studies have demonstrated that unless the variational quantum circuits are configured in a problem-specific manner, optimization of such circuits will most likely fail.In this paper, we focus on a special family of quantum circuits called the Hamiltonian variational ansatz (HVA), which its takes inspiration from the quantum approximate optimization algorithm and adiabatic quantum computation.Through the study of its entanglement spectrum and energy-gradient statistics, we find that the HVA exhibits favorable structural properties such as mild or entirely absent barren plateaus and a restricted state space that eases their optimization in comparison to the well-studied "hardware-efficient ansatz."We also numerically observe that the optimization landscape of the HVA becomes almost trap free, i.e., there are no suboptimal minima, when the ansatz is overparametrized.We observe a size-dependent "computational phase transition" as the number of layers in the HVA circuit is increased where the optimization crosses over from a hard to an easy region in terms of the quality of the approximations and the speed of convergence to a good solution.In contrast to the analogous transitions observed in the learning of random unitaries, which occur at a number of layers that grows exponentially with the number of qubits, our variational-quantum-eigensolver experiments suggest that the threshold to achieve the overparametrization phenomenon scales at most polynomially in the number of qubits for the transverse-field Ising and XXZ models.Lastly, as a demonstration of its entangling power and effectiveness, we show that the HVA can find accurate approximations to the ground states of a modified Haldane-Shastry Hamiltonian on a ring, which has long-range interactions and has a power-law entanglement scaling.
Abstract Results about entanglement (or modular) Hamiltonians of quantum many‐body systems in field theory and statistical mechanics models, and recent applications in the context of quantum information and quantum simulation, are reviewed. In the first part of the review, what is known about entanglement Hamiltonians of ground states (vacua) in quantum field theory is summarized, based on the Bisognano–Wichmann theorem and its extension to conformal field theory. This is complemented with a more rigorous mathematical discussion of the Bisognano–Wichmann theorem, within the framework of Tomita–Takesaki theorem of modular groups. The second part of the review is devoted to lattice models. There, exactly soluble cases are first considered and then the discussion is extended to non‐integrable models, whose entanglement Hamiltonian is often well captured by the lattice version of the Bisognano–Wichmann theorem. In the last part of the review, recently developed applications in quantum information processing that rely upon the specific properties of entanglement Hamiltonians in many‐body systems are summarized. These include protocols to measure entanglement spectra, and schemes to perform state tomography.
We perform a systematic investigation of variational forms (wave-function Ans\"atze), to determine the ground-state energies and properties of two-dimensional model fermionic systems on triangular lattices (with and without periodic boundary conditions), using the variational quantum eigensolver (VQE) algorithm. In particular, we focus on the nature of the entangler blocks which provide the most efficient convergence to the system ground state inasmuch as they use the minimal number of gate operations, which is key for the implementation of this algorithm in noisy intermediate-scale quantum computers. Using the concurrence measure, the amount of entanglement of the register qubits is monitored during the entire optimization process, illuminating its role in determining the efficiency of the convergence. Finally, we investigate the scaling of the VQE circuit depth as a function of the desired energy accuracy. We show that the number of gates required to reach a solution within an error $\ensuremath{\varepsilon}$ follows the Solovay-Kitaev scaling, $O[{log}_{10}^{c}(1/\ensuremath{\varepsilon})]$, with an exponent $c=1.31\ifmmode\pm\else\textpm\fi{}0.13$.