Based on the nonrelativistic QCD factorization approach, O（αsv2） corrections to J/ψ plus ηc production in e＋e- annihilation at √s= 10.6 GeV are calculated in this work. The numerical results show that the correction at αsv2 order is only about a few percent of the total theoretical result. This indicates that the perturbative expansions become convergent and that a higher order correction will be smaller. The uncertainties from the long-distance matrix elements, renormalization scale and the measurement in the experiment are also discussed. Our result is in agreement with the previous result by Jia.
