A quantum mechanics-based scoring function for halogen bonding interaction, namely XBScore(QM), is developed based on 18,135 sets of geometrical and energetical parameters optimized at M06-2X/aug-cc-pVDZ level. Applying the function on typical halogen bonding systems from Protein Data Bank demonstrates its strong ability of predicting halogen bonding as attractive interaction with strength up to -4 kcal mol(-1). With a diverse set of proteins complexed with halogenated ligands, a systematic evaluation demonstrates the integrative advantage of XBScore(QM) over 12 other scoring functions on halogen bonding in four aspects, viz. pseudo docking power, ranking power, scoring power, and genuine docking power. Thus, this study not only provides a practicable scoring function of halogen bonding for high throughput virtual screening, but also serves as a benchmark for evaluating the performance of current scoring functions on characterizing halogen bonding.