Nevertheless, the sensitivity of the results to the correlation level in calculations of Dyson orbitals does not affect our main conclusion, that the electron-dipole interactions must be taken into account for accurate modeling of PADs. We emphasize that for this effect to be prominent, the dipole moment D of the neutral residue does not have to exceed the critical value Dc for binding an electron in the dipole field. Nonetheless, as D is increased towards Dc, the parametric dependence of PADs on the former increases sharply. In the future it would be interesting to incorporate the (non-Born-Oppenheimer) effects associated with molecular rotation and examine these effects in the regime of D  Dc, which is not accessible in the current model implementation.