We employ a set of recent, theoretically motivated fits to nonperturbative unquenched gluon propagators to check on how far double gluon exchange can be used to describe the soft sector of ppscattering data (total and differential cross section). In particular, we use the refined Gribov-Zwanziger gluon propagator (as arising from dealing with the Gribov gauge fixing ambiguity) and the massive Cornwall-type gluon propagator (as motivated from Dyson-Schwinger equations) in conjunction with a perturbative quark-gluon vertex, next to a model based on the nonperturbative quark-gluon Maris-Tandy vertex, popular from Bethe-Salpeter descriptions of hadronic bound states. We compare the cross sections arising from these models with older ISR and more recent TOTEM and ATLAS data. The lower the value of total energy √s, the better the results appear to be.