“…Technically, a first-principles-based effective Hamiltonian 34 is used within Monte-Carlo (MC) simulations to determine the energetics and local electric dipoles in each perovskite five-atom cell of these supercells. The validity of this approach was demonstrated by previous theoretical studies of ultrathin PZT films under compressive strains that (1) yield 180 degrees up and down stripe domains that periodically alternate along [100] (or along [010]) for their ground-state 34,35 , in agreement with experimental observation 36 (note that 'up' (respectively, 'down') domains refer to domains in which the z-component of the dipole is parallel (respectively, antiparallel) to the z-axis, respectively); (2) predict a linear dependency between the width of these periodic stripes and the square root of the film's thickness 37 , as consistent with recent measurements 38 ; and (3) have also led to the prediction of various topological defects such as vortices 39 , dipolar waves 40 , bubbles 41 and merons (or convex disclinations) 42 in ferroelectrics, that have been experimentally confirmed [42][43][44] . Note that the predicted temperature has to be rescaled by a factor of ∼ 1.6 with respect to measurements 45 .…”