In the spectator models of the nucleon with scalar and axial-vector diquarks, we show the effect of Pauli coupling in the photon-quark vertex to the parton distribution functions (PDFs) of nucleon and azimuthal asymmetry in the unpolarized semi-inclusive deep inelastic scattering (SIDIS). This anomalous coupling gives an obvious contribution to unpolarized and polarized PDFs and also leads to a cos 2ϕ h azimuthal asymmetry proportional to the squared Pauli form factor, due to the helicity flip of the struck quark. After determining the model parameters by fitting PDFs to the global fits, this new distribution for cos 2ϕ h asymmetry is given numerically. In the framework of transverse momentum dependence (TMD), we find that it is positive and of a few percent in the kinematical regime of HERMES and COMPASS Collaborations, in the same order of magnitude as the Cahn effect.