Using the self-consistent Landau-Ginzburg-Devonshire approach we simulate and analyze the spontaneous formation of the domain structure in thin ferroelectric films covered with the surface screening charge of the specific nature (Bardeen-type surface states). Hence we consider the competition between the screening and the domain formation as alternative ways to reduce the electrostatic energy and reveal unusual peculiarities of distributions of polarization, electric and elastic fields conditioned by the surface screening length and the flexocoupling strength. We have established that the critical thickness of the film and its transition temperature to a paraelectric phase strongly depend on the Bardeen screening length, while the flexocoupling affects the polarization rotation and closure domain structure and induces ribbon-like nano-scale domains in the film depth far from the top open surface. Hence the joint action of the surface screening (originating from e.g. the adsorption of ambient ions or surface states) and flexocoupling may remarkably modify polar and electromechanical properties of thin ferroelectric films. (A.N.M.) 2
I. INTRODUCTIONFerroelectric materials remain the object of endless fascination for applied and fundamental science alike. The fundamental aspect of these materials is the presence of surface and interface bound charges due to the discontinuity of the spontaneous polarization. Since the early days of ferroelectricity, these charges were recognized as the key aspect of the physics of ferroelectric surfaces and interfaces.Indeed, if the polarization charge were uncompensated, it would provide bulk-like contributions to the free energy of materials, i.e. the corresponding energy would diverge with the system size. These considerations stimulated the search for mechanisms for lowering of this depolarization energy. One such mechanism is formation of ferroelectric domains, extensively analyzed in classical textbooks [1,2 ]. The second is charge screening, either internal "bulk" screening by free charges inside a ferroelecric [3,4] or external "surface" screening by the free charges in the case of the open or electroded ferroelectric film surface. Surprisingly, until now the domain formation and surface screening mechanisms were considered separately, and the competition between these effects almost escaped the attention of scientific community, despite the fact that these processes are intertwined in thin films.Surface screening of the bound charges is typically provided by the mobile charges adsorbed from the ambience in the case of high or normal humidity [5,6,7,8,9] or by internal mobile charges of defect nature [10,11]. In a specific case of the very weak screening, or its artificial absence due to the experimental treatment (cleaned surface in dry atmosphere, ultra-high vacuum or thick dielectric layer at the surface) the screening charges can be localized at surface states caused by the strong bandbending by depolarization field [12,13,14,15,16,17]. For both aforementioned cases the s...