Proton affinity is a major factor in the atmospheric pressure chemical ionization of illicit drugs. The detection of illicit drugs by mass spectrometry and ion mobility spectrometry relies on the analytes having greater proton affinities than background species. Evaluating proton affinities for fentanyl and its analogues is informative for predicting the likelihood of ionization in different environments and for optimizing the compounds' ionization and detection, such as through the addition of dopant chemicals. Herein, density functional theory was used to computationally determine the proton affinity and gas-phase basicity of 15 fentanyl compounds and several relevant molecules as a reference point. The range of proton affinities for the fentanyl compounds was from 1018 to 1078 kJ/mol. Fentanyl compounds with the higher proton affinity values appeared to form a bridge between the oxygen on the amide and the protonated nitrogen on the piperidine ring based on models and calculated bond distances. Experiments with fragmentation of proton-bound clusters using atmospheric flow tube-mass spectrometry (AFT-MS) provided estimates of relative proton affinities and showed proton affinity values of fentanyl compounds >1000 kJ/mol, which were consistent with the computational results. The high proton affinities of fentanyl compounds facilitate their detection by ambient ionization techniques in complex environments. The detection limits of the fentanyl compounds with AFT-MS are in the low femtogram range, which demonstrates the feasibility of trace vapor drug detection.