Effect of surface ionic screening on the polarization reversal scenario in ferroelectric thin films: Crossover from ferroionic to antiferroionic states

Nonlinear electrostatic interaction between the surface ions of electrochemical nature and ferroelectric dipoles gives rise to the coupled ferroionic states in nanoscale ferroelectrics. Here, we investigate the role of the surface ion formation energy on the polarization states and its reversal mechanisms, domain structure, and corresponding phase diagrams of ferroelectric thin films. Using 3D finite element modeling, we analyze the distribution and hysteresis loops of ferroelectric polarization and ionic charge, and the dynamics of the domain states. These calculations performed over large parameter space delineate the regions of single-and polydomain ferroelectric, ferroionic, antiferroionic, and nonferroelectric states as a function of surface ion formation energy, film thickness, applied voltage, and temperature. We further map the analytical theory for 1D systems onto an effective Landau-Ginzburg free energy and establish the correspondence between the 3D numerical and 1D analytical results. This approach allows us to perform an overview of the ferroionic system phase diagrams and explore the specifics of polarization reversal and domain evolution phenomena.