The bipolar polarization hysteresis loops reveal the expected ferroelectric nature of the samples, without any notable leakage-current contributions. The subtle differences between the loops can be quantified by using the squareness parameter, 𝑅 𝑆𝑞 (Eq. 6 ), and the full-width at half maximum of the current peak around the coercive field, IFWHM.