Stabilization of specific crystal polymorphs of an active pharmaceutical ingredient is crucial for preventing uncontrollable interconversion of various crystalline forms, which affects physicochemical properties as well as physiological activity. Co-crystallization with various excipients is an emerging productive way of achieving such stabilization in the solid state. In this work, we identified an opportunity for co-crystallization of antiviral drug nevirapine (NVP) with a classical XB donor, 1,2,4,5-tetrafluoro-3,6-diiodobenzene (1,4-FIB), as well as 1,3-diiodobenzene (1,3-DIB), which has been seldom employed as an XB donor to date. In the X-ray structures of NVP·1,4-FIB and NVP·1,3-DIB co-crystals, different hydrogen and halogen bonding modes were detected and further investigated via DFT calculations as well as topological analysis of the electron density distribution within the framework of the QTAIM method at the M06/DZP-DKH level of theory. Estimated energies of these supramolecular contacts vary from 0.6 to 5.7 kcal/mol.