This small python3 script is about postprocessing of CP2K AIMD output with BEC, for polarization changes during ferroelectric switching under external electric fields.
Example: The abnormal continuous ferroelectric switching (homogeneous) of Pca21-HfO2, via UCP-aN (N = 2, 4) and UCP-b paths[1,2]. This figure is plotted in the general scheme of P-E loop, if it was.
Here will come some introductions…