Ab initio molecular dynamics* to study the dissociation of 2-thiouracil and its tautomers for low-energy electron
Zhang Ying,Wang Xing,Xu Zhong-feng,Ren Jie-ru,Zhang Yan-ning,Zhou xian-ming,Liang chang-hui,Zhang xiao-an,,,
DOI: https://doi.org/10.7498/aps.73.20231304
IF: 0.906
2024-01-01
Acta Physica Sinica
Abstract:When biomolecules interact with high-energy particles and rays, they are directly ionized or dissociated, after which a large number of low-energy electrons are formatted as secondary particles. These low-energy electrons will attach to biomolecules, and trigger secondary dissociation, form free radicals and ions with high reactivity, which can damage the structure and function of the biomolecule and cause irreversible radiation damage to the biomolecule. It is important to study the low-energy dissociative electron attachment (DEA) processes of biomolecules for understanding radiation damage in biological organisms. Currently, the theoretical studies of DEA have mainly focused on the bound states of negative ions and the types of resonances during dissociation processes. While the dissociation process is well described by quantum computational methods, the diversity and complexity of dissociation channels present in the dissociation process of 2-thiouracil molecule also poses a great computational challenge for these methods. In addition, the quantum computational methods are not ideal for dealing with the discrete states of chemical bonds and the problem of continuity coupling of electrons. The dissociation dynamics of biomolecules are mainly caused by ionization and electron attachment. Ab initio molecular dynamics simulation can reasonably describe these processes. In the light of these considerations, Ab initio molecular dynamics simulation is used in this work to study dynamic variation processes in DEA. The low-energy electrons dissociative attachment to 2-thiouracil in the gas phase is studied by using the Born-Oppenheimer molecular dynamics model combined with density functional theory. It is found that an important dehydrogenation phenomenon of 2-thiouracil and its tautomers during DEA processes, the N-H and C-H bond are broken at specific locations. Due to loss the hydrogen atoms at the N and C sites, forming the closed-shell dehydrogenated negative ion (TU-H) - , which is the most important negative ion fragments in the dissociation processes. The potential energy curves, the bond dissociation energy and the electron affinity energy of the broken bond show the N-H bond is the most likely to break, indicating the formation of the negative ion (TU-H) - is mainly from the N-H bond broken. The results of the theoretical calculations in this work are in good agreement with the available experimental results, indicating that the chosen calculation method is fully reliable. BOMD simulations can not only dynamically recover the process of dissociative attachment of low-energy electrons to 2-thiouracil, but, more importantly, provide insight into the mechanisms of dehydrogenation and dissociation channels of 2-thiouracil molecules in DEA processes.
physics, multidisciplinary