Exploring the limit of the Lattice-Bisognano-Wichmann form describing the Entanglement Hamiltonian: A quantum Monte Carlo study

The entanglement Hamiltonian (EH) encapsulates the essential entanglement properties of a quantum many-body system and serves as a powerful theoretical construct. 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 complex quantum many-body systems.