This work discusses applications of computational simulations to a wide variety of chemical systems, to investigate intermolecular interactions to develop force field parameters and gain new insights into chemical reactivity and structure stability. First, we cover the characterization of hydrogen-bonding interactions in pyrazine tetracarboxamide complexes employing quantum topological analyses. Second we describe the use of quantum mechanical energy decomposition analysis (EDA) and non-covalent interactions (NCIs) analysis to investigate hydrogen-bonding and intermolecular interactions in a series of representative 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([bmim][Tf2N]) ion pairs extracted from classical equilibrium and non-equilibrium molecular dynamics simulations. Thirdly, we describe the use of multipolar/polarizable AMOEBA force field to study the extraction of benzene from a gasoline model employing 1,3-dimethylimidazolium tetrafluorobrorate, [DMIM][BF4], and ethylmethylimidazolium tetrafluorobrorate, [EMIM][BF4]. Fourthly, we cover the recent improvements and new capabilities of the QM/MM code "LICHEM". Finally, we describe the use of polarizable ab initio QM/MM calculations and study the reaction mechanism of N-tert-butyloxycarbonylation of aniline in [EMIm][BF4], and ground state destabilization in uracil DNA glycosylase (UDG).