sys = read_in {
    int_file = "INTDUMP",
}

targets = {2*10^3, 4*10^3, 8*10^3, 1.6*10^4, 3.2*10^4, 6.4*10^4, 1.28*10^5}

for i,target in ipairs(targets) do
    fciqmc {
        sys = sys,
        qmc = {
            tau = 0.01,
            rng_seed = 8,
            init_pop = target/20,
            mc_cycles = 5,
            nreports = 3*10^3,
            tau_search = true,
            target_population = target,
            excit_gen = "heat_bath",
            initiator = true,
            real_amplitudes = true,
            spawn_cutoff = 0.1,
            state_size = -1000,
            spawned_state_size = -100,
        },
    }
end
