shows $w_{2}\left(T\right)$ with fitting parameters $\gamma=1.66$ and $T_{O}=256.7$ K. Similar to the PZN-0.13PT case, the results here also shows $\chi_{2}\ll\chi_{1}$.

Discussion\[subsec:Discussion\]
===============================

In the statistical model we divide the dipoles inside ferroelectrics relaxors into two groups, one group being confined in potential wells, while the other having can overcome the potential confinement and exhibiting a more vibrant dynamics. It has been demonstrated that the Maxwell-Boltzmann distribution plays a significant role in determining the profile of $\chi\left(T\right)$. To address a particular type of relaxor, one may only need to adjust the dielectric response function associated with each group of dipoles, while keeping the rest unchanged. In this section, we discuss a few issues of this approach and its limitations.

Characteristic temperature $T_{m}$
----------------------------------

The present analysis helps us to understand why the susceptibility of a relaxor reaches its peak value at some temperature, $T_{m}$. For BZT, the function $\chi_{1}\mathcal{L}\left(\frac{\theta}{T}\right)$ [\[]{}see Eq. (\[eq:BZT-susceptibility\]) and Fig. \[fig:components\][\]]{} or $\chi_{1}/\left[1+b\exp\left(-\theta/T\right)\right]$ [\[]{}see Eq. (\[eq:BZT-frequency\])[\]]{} describes the contribution to susceptibility from dipoles with kinetic energy higher than $E_{b}$. These two functions are both monotonically decreasing with $T$, reflecting the fact that thermal motions prevent dipoles from aligning with the applied electric field. On the other hand, 