Parameter Scans

Generate a directory tree of simulation files (one directory per parameter point), then run with the queue system (see dqmc-run), then analyze (see dqmc-analyze).

Temperature Scan

Vary L while adjusting dt to maintain Trotter error bound:

from dqmc_util import gen_1band_hub
import numpy as np

U = 4.0
step = 5  # L must be divisible by n_matmul and period_eqlt (defaults: 5)
for T in [0.1, 0.2, 0.5, 1.0]:
    beta = 1.0 / T
    dt = min((0.05/U)**0.5, beta / 10)
    L = int(np.ceil(beta / dt / step) * step)
    dt = beta / L