We present the calculation of the dominant next to leading order QCD corrections to Higgs boson production in association with three jets via vector boson fusion in the form of a NLO parton-level Monte Carlo program. QCD corrections to integrated cro
ss sections are modest, while the shapes of some kinematical distributions change appreciably at NLO. Scale uncertainties are shown to be reduced at NLO for the total cross section and for distributions. We consider a central jet veto at the LHC and analyze the veto probability for typical vector boson fusion cuts. Scale uncertainties of the veto probability are sufficiently small at NLO for precise Higgs coupling measurements at the LHC.
