“…) of an oper-atorÔ for a set of imaginary times τ i (i = 0, 1, · · · , N τ ) with statistical errors can be obtained from DQMC simulations. By SAC method [48][49][50][51], the corresponding real-frequency spectral function A(ω) can be obtained from them according to the relationship of G(τ) = ∫ ∞ −∞ dωA(ω)K(τ, ω), where the kernel K(τ, ω) depends on the type of the spectral function, i.e., fermionic or bosonic, finite or zero temperature.…”