or interacting with the relaxation modes, it becomes necessary to use coupled modes to model the dielectric response [@Wang2014; @Wang2016]. On the other hand, there are also many investigations on how the dielectric response, $\chi$, depends on the temperature, $T$, at given frequencies. In addition to the well known Curie law for $\chi\left(T\right)$ at high temperature, most useful equation for fitting around the dielectric peak appears to be the square law. In particular, the formula $$\begin{aligned}
\frac{1}{\varepsilon\left(T\right)}= & \frac{1}{\varepsilon_{A}}+\frac{\left(T-T_{A}\right)^{^{\eta}}}{B}\label{eq:square-law}\end{aligned}$$ was proposed to describe the permittivity at $T>T_{m}$ [@Smolenskii1970; @Kirillov1973]. Initially, $\eta$ was found to be 2, but later was shown to be between 1 and 2 [@Clarke1974; @Uchino1982; @Santos2001; @Bokov2006a; @Bokov2012]. Here, we further the investigation in this direction and attempt to address some important questions regarding relaxor behaviors. We will explain why the dielectric constant has a peak value at $T_{m}$, and what causes the asymmetry around the peak. Moreover, by constructing a statistical model that properly describes how dipoles behave in relaxors, we propose formula to fit experimental results, which further illuminates the physics behind relaxation behavior.

This paper is organized as the follows. In Sec. \[sec:Statistical-modeling\], we introduce the statistical model. In Sec. \[sec:Results\], we apply this 