We present an automatized workflow which, starting from molecular dynamics simulations, identifies reaction events, filters them, and prepares them for accurate quantum chemical calculations using, for example, Density Functional Theory (DFT) or Coupled Cluster methods. The capabilities of the automatized workflow are demonstrated by the example of simulations for the combustion of some polycyclic aromatic hydrocarbons (PAHs). It is shown how key elementary reaction candidates are filtered out of a much larger set of redundant reactions and refined further. The molecular species in question are optimized using DFT and reaction energies, barrier heights, and reaction rates are calculated. The setup is general enough to include at this stage configurational sampling, which can be exploited in the future. Using the introduced machinery, we investigate how the observed reaction types depend on the gas atmosphere used in the molecular dynamics simulation. For the re-optimization on the DFT level, we show how the additional information needed to switch from reactive force-field to electronic structure calculations can be filled in and study how well ReaxFF and DFT agree with each other and shine light on the perspective of using more accurate semi-empirical methods in the MD simulation.

Using diode laser vibration-rotation-tunneling spectroscopy near 15 Thz (500 cm−1), we have measured and assigned 142 transitions to three a-type librational subbands of the water hexamer-d12 prism. These subbands reveal dramatically enhanced (ca. 1000×) tunneling splittings relative to the ground state. This enhancement is in agreement with that observed for the water dimer, trimer, and pentamer in this same frequency region. The water prism tunneling motion has been predicted to potentially describe the motions of water in interfacial and confined environments; hence, the results presented here indicate that excitation of librational vibrations has a significant impact on the hydrogen bond dynamics in these macroscopic environments.