Oksana Tsendra, Andrea Michalkova Scott, Leonid Gorb, Adrian Daniel Boese, Frances C. Hill, Mykola M. Ilchenko, Danuta Leszczynska, Jerzy Leszczynski

• A cluster approach extended to the ONIOM methodology has been applied using several density functionals and Moller-Plesset perturbation theory (MP2) to simulate the adsorption of selected nitrogen-containing compounds [NCCs, 2,4,6-trinitrotoluene (TNT), 2,4-dinitrotoluene (DNT), 2,4-dinitroanisole (DNAN), and 3-nitro-1,2,4-triazole-5-one (NTO)] on the hydroxyated (100) surface of a-quartz. The structural properties were calculated using the M06-2X functional and 6-31G(d,p) basis set. The M06-2X-D3, PBE-D3, and MP2 methods were used to calculate the adsorption energies. Results have been compared with the data from other studies of adsorption of compounds of similar nature on silica. Effect of deformation of the silica surface and adsorbates on the binding energy values was also studied. The atoms in molecules (AIM) analysis was employed to characterize the adsorbate-adsorbent binding and to calculate the bond energies. The silica surface shows different sorption affinity toward the chemicals considered depending on their electronicA cluster approach extended to the ONIOM methodology has been applied using several density functionals and Moller-Plesset perturbation theory (MP2) to simulate the adsorption of selected nitrogen-containing compounds [NCCs, 2,4,6-trinitrotoluene (TNT), 2,4-dinitrotoluene (DNT), 2,4-dinitroanisole (DNAN), and 3-nitro-1,2,4-triazole-5-one (NTO)] on the hydroxyated (100) surface of a-quartz. The structural properties were calculated using the M06-2X functional and 6-31G(d,p) basis set. The M06-2X-D3, PBE-D3, and MP2 methods were used to calculate the adsorption energies. Results have been compared with the data from other studies of adsorption of compounds of similar nature on silica. Effect of deformation of the silica surface and adsorbates on the binding energy values was also studied. The atoms in molecules (AIM) analysis was employed to characterize the adsorbate-adsorbent binding and to calculate the bond energies. The silica surface shows different sorption affinity toward the chemicals considered depending on their electronic structure. All target NCCs are physisorbed on the modeled silica surface. Adsorption occurs due to the formation of multiple hydrogen bonds between the functional groups of NCCs and surface silanol groups. Parallel orientation of NCCs interacting with the silica surface was found to be favorable when compared with perpendicularly oriented NCCs. NTO was found to be the most strongly adsorbed on the silica surface among all of the considered compounds. Dispersion correction was shown to play an important role in the DFT calculations of the adsorption energies of silica-NCC systems.…