Consider the long-range models on Z d of random walk, self-avoiding walk, percolation and the Ising model, whose translation-invariant 1-step distribution/coupling coefficient decays as |x| −d−α for some α > 0. In the previous work [15], we have shown in a unified fashion for all α = 2 that, assuming a bound on the "derivative" of the n-step distribution (the compound-zeta distribution satisfies this assumed bound), the critical two-point function G pc (x) decays as |x| α∧2−d above the uppercritical dimension d c ≡ (α ∧ 2)m, where m = 2 for self-avoiding walk and the Ising model and m = 3 for percolation.In this paper, we show in a much simpler way, without assuming a bound on the derivative of the n-step distribution, that G pc (x) for the marginal case α = 2 decays as |x| 2−d / log |x| whenever d ≥ d c (with a large spread-out parameter L). This solves the conjecture in [15], extended all the way down to d = d c , and confirms a part of predictions in physics [10]. The proof is based on the lace expansion and new convolution bounds on power functions with log corrections.