Automated quantum chemistry based molecular dynamics simulations of electron ionization induced fragmentations of the nucleobases uracil, thymine, cytosine, and guanine

The gas-phase decomposition pathways of electron ionization (EI)-induced radical cations of the nucleobases uracil, thymine, cytosine, and guanine are investigated by means of mixed quantum-classical molecular dynamics. No preconceived fragmentation channels are used in the calculations. The results compare well to a plethora of experimental and theoretical data for these important biomolecules. With our combined stochastic and dynamic approach, one can access in an unbiased way the energetically available decomposition mechanisms. Additionally, we are able to separate the EI mass spectra of different tautomers of cytosine and guanine. Our method (previously termed quantum chemistry electron ionization mass spectra) reproduces free nucleobase experimental mass spectra well and provides detailed mechanistic in-sight into high-energy unimolecular decomposition processes.