  The pion scalar, vector and tensor form factors are calculated within a
symmetry-preserving contact interaction model (CI) of quantum chromodynamics
(QCD), encompassed within a Dyson-Schwinger and Bethe-Salpeter equations
approach. In addition to the traditional rainbow-ladder truncation, a modified
interaction kernel for the Bethe-Salpeter equation is adopted. The implemented
kernel preserves the vector and axial-vector Ward-Takahashi identities, while
also providing additional freedom. Consequently, new tensor structures are
generated in the corresponding interaction vertices, shifting the location of
the mass poles appearing in the quark-photon and quark tensor vertex and
yielding a notorious improvement in the final results. Despite the simplicity
of the CI, the computed form factors and radii are compatible with recent
lattice QCD simulations.
