Quantum simulation in imaginary time for gauge-invariant models

This repository contains Wolfram Language code and Jupyter Notebooks for simulating gauge-invariant models in imaginary time, computing thermal averages.

Ground–State Evaluation

The package allows you to obtain both the approximate ground–state energy (from the circuit-based imaginary–time method) and the exact ground–state energy (from diagonalization of the stored Hamiltonian operator). Check gauge_simulation_run.py for other accessible quantities.

from gauge_simulation import models

H = models.build_hamiltonian(model = "tfim", lattice_size=[1,3])

exact_groundstate = H.exact_ground_state()

approx_groundstate = H.approximate_ground_state()

print(f"Exact:  {exact_groundstate:.4f}")
print(f"Approx: {approx_groundstate:.4f}")

Exact:  -2.7441
Approx: -2.5255

The bond parameters are currently set to π/4.