As analogies to magnetic skyrmions, polar skyrmions in ferroelectric superlattices and multilayers have garnered widespread attention for their non-trivial topology and novel properties like negative capacitance and nonlinear optical effect. So far, they have only been theoretically predicted to be able to assemble ordered hexagonal skyrmion lattices (SkLs) in ferroelectric thin films. Here, based on phase-field simulations, we report the critical roles of flexoelectricity playing in the stabilization and transformation of polar SkLs. Different polar SkL patterns can emerge in the ferroelectric thin films, including tetragonal-SkL, and hexagonal-SkLs with diverse orientations, as summarized by phase diagrams. These emergent SkL states are attributed to the material anisotropy modified by the flexoelectric effect. Interestingly, we further found that the hexagonal-SkLs can be rotated by applying strain gradient or in-plane electric field to the films. Moreover, a nonreciprocal bending response of tetragonal-SkL is also induced by the flexoelectric effect. Our results provide useful guidelines for the implementation of polar skyrmion lattices in experiments.
