The aqueous solutions of bisulfite (SO3H-) and sulfonate (HSO3-) were simulated by the ab initio quantum mechanical charge field molecular dynamics (QMCF MD) formalism. All superimposed trajectories for the atomic coordinates of solutes with three-dimensional alignment here illustrated the reactivities of the ions. Power spectra were evaluated on the basis of the velocity autocorrelation functions (VACFs) with the normal-mode analysis, presenting a higher frequency of the symmetric SO3 deformation (delta(s)(SO3)) than the asymmetric SO3 deformation (delta(as)(SO3)) modes for the sulfonate ion. The different influence of solvent on the frequency of the 0 H and S H stretching suggests a higher stability of hydrated sulfonate ion. The bisulfite shows a slightly stronger molecular hydration shell than the sulfonate ion with the average number of ion solvent hydrogen bonds (H-bonds) of 5.3 and 5.0, respectively. Extra water molecules within the molecular hydration shell are found for bisulfite (1.2) and for sulfonate (1.6). The mean residence times for the water ligands classify each ion as a structure maker, while the S-H bond within the sulfonate ion displays a hydrophobic behavior. No tautomerization was observed within the simulation period.