A model to quantitatively predict ion abundances from atmospheric pressure chemical ionization (APCI) between hydrated protons and a volatile organic compound (VOC) was extended to binary mixtures of VOCs. The model includes differences in vapor concentrations, rate coefficients, and reaction times and is enhanced with cross reactions between neutral vapors and protonated monomers. In this model, two specific VOCs were considered, a ketone, 6-methyl-5-hepten-2-one (M, and an amine, 2,6-di-tert-butyl-pyridine (N), with measured "conditional rate coefficients" (in cm3·s-1) of kM = 1.11 × 10-9 and kN = 9.17 × 10-10, respectively. The cross reaction of MH+(H2O)x to NH+(H2O)y was measured as kcr = 1.31 × 10-12 at 60 °C. Cross reactions showed an impact on ion abundances at t > 30 ms for equal vapor concentrations of 100 ppb for M and N. In contrast, this impact was negligible for vapor concentrations of 1 ppb and did not exceed 5% change in product ion abundance up to 1000 ms reaction times. The model was validated with laboratory measurements to within ∼10% using an ion mobility spectrometer and effective reaction time obtained from computational fitting of experimental findings. This was necessitated by complex flow patterns in the ion source volume and was determined as ∼10.5 ms. The model has interpretative and predictive value for quantitative analysis of responses with ambient pressure ion sources for mass spectrometry and ion mobility spectrometry.