A universal black-box quantum Monte Carlo approach to quantum phase transitions

We derive exact, universal, closed-form quantum Monte Carlo estimators for finite-temperature energy susceptibility and fidelity susceptibility, applicable to essentially arbitrary Hamiltonians. Combined with recent advancements in Monte Carlo, our approach enables a black-box framework for studying quantum phase transitions–without requiring prior knowledge of an order parameter or the manual design of model-specific ergodic quantum Monte Carlo update rules. We demonstrate the utility of our method by applying a single implementation to the transverse-field Ising model, the XXZ model, and an ensemble of models related by random unitaries.

💡 Research Summary

This paper introduces a groundbreaking “black-box” framework for studying quantum phase transitions (QPTs) using quantum Monte Carlo (QMC) simulations. The core achievement is the derivation of exact, universal, closed-form QMC estimators for two key quantities: the energy susceptibility (ES) and the fidelity susceptibility (FS). These estimators are applicable to essentially arbitrary Hamiltonians, removing a major limitation of previous methods which often required the perturbation term to have a specific form (e.g., proportional only to diagonal or off-diagonal parts of the Hamiltonian).

The power of this approach lies in its integration with the recently developed Permutation Matrix Representation QMC (PMR-QMC) framework. A key advance in PMR-QMC is its ability to automatically generate ergodic Monte Carlo update rules that satisfy detailed balance for a wide variety of Hamiltonians, including arbitrary spin-1/2, bosonic, and fermionic systems. By combining these auto-generated updates with the new universal estimators, the method eliminates the need for manual, model-specific tuning or prior knowledge of an order parameter. This creates a truly automated, black-box tool for probing QPTs.

The authors rigorously demonstrate the framework’s utility and accuracy through a series of numerical experiments on diverse models. First, they validate the method on a simple two-spin model with an avoided level crossing, showing perfect agreement with exact numerical calculations for the FS across different temperatures and its convergence to the ground-state value. Second, they apply it to the well-studied 2D transverse-field Ising model (TFIM). They successfully reproduce known results for ES and FS computed via other QMC methods and demonstrate the ability to compute these quantities within a fixed parity subspace—a methodological convenience highlighting PMR-QMC’s versatility.

Third, the framework is applied to the 2D XXZ model to investigate a finite-temperature phase transition. By computing ES and FS for increasing system sizes, the method correctly identifies the critical temperature, confirming its capability for finite-size scaling studies relevant to critical phenomena. The most compelling demonstration involves an ensemble of 100-spin models generated by applying random unitary rotations to the simple two-spin Hamiltonian. Each resulting Hamiltonian is highly complex, comprising hundreds of random Pauli operator terms, making it intractable for traditional QMC approaches. Remarkably, the single, unmodified implementation accurately estimates the imaginary-time correlator, ES, and FS for this ensemble, matching the expected exact results. This experiment powerfully showcases the method’s universality and black-box nature.

Finally, the paper discusses computational performance. The algorithm is memory-efficient and exhibits near-perfect strong scaling for parallel execution. Empirical tests show that simulation wall-clock time scales approximately as N^2 (with N being the number of spins). While the FS estimator has a higher computational complexity per measurement than ES, its impact on total runtime can be managed by adjusting the measurement frequency. In summary, this work provides a robust, automated, and general-purpose tool that significantly lowers the barrier to studying quantum phase transitions, particularly for novel or complex systems where order parameters are unknown. It represents a significant step towards the automation of quantum many-body simulations.