A. Models of antiferroelectricsComplexity of antiferroelectrics resulted in a rather large number of different approaches and perspectives in their modeling. To help the reader navigate that, a brief systematization follows.First model of antiferroelectrics is due to Kittel [1], who reduced an (at that time -abstract) antiferroelectric to two inter-penetrating lattices with anti-parallel polarizations. The most studied real antiferroelectric, PbZrO 3 , is much more complex than that, in part -due to the more complex dipole ordering ↑↑↓↓ [24], due to the possibility of alternative incommensurate orderings [9] and due to the presence of anti-phase octahedral tilts as a second type of distortion [24], usually coexisting with the anti-polar ordering of cations.Modern models fall, roughly, into three categories. First -ab-initio analysis [25][26][27] of energy landscape at zero Kelvin, sometimes extended to modeling of finite temperature behavior via molecular dynamics [22,28,29]. Usually this extension needs a parametrization of energy in terms of local ionic displacements around