The infrared (IR) spectroscopy and ligand migration of photodissociated nitric oxide (NO) in and around the active sites in myoglobin (Mb) are investigated. A distributed multipolar model for open-shell systems is developed and used, which allows one to realistically describe the charge distribution around the diatomic probe molecule. The IR spectra were computed from the trajectories for two conformational substates at various temperatures. The lines are narrow (width of 3–7 cm^–1 at 20–100 K), in agreement with the experimental observations where they have widths of 4–5 cm^–1 at 4 K. It is found that within one conformational substate (B or C) the splitting of the spectrum can be correctly described compared with recent experiments. Similar to photodissociated CO in Mb, additional substates exist for NO in Mb, which are separated by barriers below 1 kcal/mol. Contrary to full quantum mechanical calculations, however, the force field and mixed QM/MM simulations do not correctly describe the relative shifts between the B- and C-states relative to gas-phase NO. Free energy simulations establish that NO preferably localizes in the distal site and the barrier for migration to the neighboring Xe4 pocket is ΔG_B→C = 1.7–2.0 kcal/mol. The reverse barrier is ΔG_B←C = 0.7 kcal/mol, which agrees well with the experimental value of 0.7 kcal/mol, estimated from kinetic data.