Then, the ϕ is fitted with two exponential decays and two damped harmonic oscillators as basis functions to include not only the Debye relaxations but also the high-frequency resonances, which have been identified in earlier experimental works . To perform the fitting, we employ a Markov chain Monte Carlo algorithm implemented in PyOECP to obtain the fitting parameters for the discrete relaxation model, which is given as