We calculate the next-to-leading order (NLO) quantum chromodynamics (QCD) corrections to double charmonium production processes e + e − → γ * → ηc+hc(1P )/ψ1,2(1D) within the non-relativistic QCD (NRQCD) factorization framework. We find that the corrections to ηc+hc(1P ) production are positive, while those to ηc+ψ1,2(1D) are negative. Unlike the J/ψ + ηc case, all the corrections here are not large. Uncertainties in the renormalization scale, quark mass and running energy of center-of-mass are discussed, and the scale dependence of these processes is found to be greatly reduced with the NLO QCD corrections.