In Eq. 3, i is the imaginary unit number (i 2 = -1), ω is the angular frequency, and t is time. ϵ(0) and ϵ ∞ are the low-frequency and high-frequency limits, respectively. Since ϕ(t) is subject to thermal fluctuations, its derivative is easily corrupted. Earlier works usually employ discrete relaxation models or perform the simulation repeatedly to reduce the noise. The former method can distort the overall dielectric behavior, while the later one is time-consuming.