Motivated by TeV-scale neutrino mass models, we propose a systematic treatment of heavy neutrino (N ) production at hadron colliders. Our simple and efficient modeling of the vector boson fusion (VBF) W ± γ → N ℓ ± and N ℓ ± + nj signal definitions resolve collinear and soft divergences that have plagued past studies, and is applicable to other color-singlet processes, e.g., associated Higgs (W ± h), sparticle (l ±ν ℓ ), and charged Higgs (h ±± h ∓ ) production. We present, for the first time, a comparison of all leading N production modes, including both gluon fusion (GF)ν ℓ and VBF. We obtain fully differential results up to next-to-leading order (NLO) in QCD accuracy using a Monte Carlo tool chain linking FeynRules, NLOCT, and MadGraph5 aMC@NLO. Associated model files are publicly available. At the 14 TeV LHC, the leading order GF rate is small and comparable to the NLO N ℓ ± + 1j rate; at a future 100 TeV Very Large Hadron Collider, GF dominates for mN = 300 − 1500 GeV, beyond which VBF takes lead.