A moderately strong vector repulsion between quarks in dense quark matter is needed to explain how a quark core can support neutron stars heavier than two solar masses. We study this repulsion, parametrized by a four-fermion interaction with coupling gV , in terms of non-perturbative gluon exchange in QCD in the Landau gauge. Matching the energy of quark matter, gV n 2 q (where nq is the number density of quarks) with the quark exchange energy calculated in QCD with a gluon propagator parametrized by a finite gluon mass mg and a frozen coupling αs, at moderate quark densities, we find that gluon masses mg in the range 200 -600 MeV and αs = 2 -4 lead to a gV consistent with neutron star phenomenology. Estimating the effects of quark masses and a colorflavor-locked (CFL) pairing gap, we find that gV can be well approximated by a flavor-symmetric, decreasing function of density. We briefly discuss similar matchings for the isovector repulsion and for the pairing attraction.