relaxors remains an open question.

Third, at high temperature, Curie law is observed in many ferroelectrics. For the static susceptibility of BZT [\[]{}Fig. \[fig:Fitting-susceptibility-of\][\]]{}, this law can be recovered from the proposed equation, Eq. (\[eq:BZT-susceptibility\]). On the other hand, for Eqs. (\[eq:BZT-frequency\]) and (\[eq:PZN\]), the Curie law cannot be directly recovered. For Eq. (\[eq:PZN\]), we have the asymptotic relation $\text{\ensuremath{\chi}}\sim A/\left(T-T_{C}\right)^{\gamma}+B/T^{3/2}$ at very large $T$. It is unclear how good this relation can fit the Curie law. Therefore, in fitting experimental data at high temperatures, one needs to bear in mind that Eqs. (\[eq:BZT-frequency\]) and (\[eq:PZN\]) should be used with care.

Finally, with Eqs. (\[eq:BZT-susceptibility\]) and (\[eq:PZN\]), in principle we can obtain the relation between $T_{m}$ and $\nu$, which can then be compared to the well-known Vogel-Fulcher law [@Wang2016]. However, we have failed to obtain analytical expressions for $T_{m}\left(\nu\right)$ and believe numerical calculation seems to be the only feasible way to establish the relation between $T_{m}$ and $\nu$.

Conclusion\[sec:Conclusion\]
============================

Instead of working on the atomic level, the present work employs a macroscopic statistical approach to help understanding dielectric properties of relaxors. The effects of disorder, PNRs, and random electric fields are considered statistically by introducing the average potential well, 