Search

A comprehensive computational study of the gas-phase photodissociation dynamics of methanol is presented. Using a multiconfigurational active space based method (RASSCF) to obtain multidimensional potential energy surfaces (PESs) on-the-fly, direct quantum dynamics simulations were run using the variational multi-configurational Gaussian method (DD-vMCG). Different initial excitation energies were simulated to investigate the dependence of the branching ratios on the electronic state being populated. A detailed mechanistic explanation is provided for the observed differences with respect to the excitation energy. Population of the lowest lying excited state of methanol leads to rapid hydroxyl hydrogen loss as the main dissociation channel. This is rationalized by the strongly dissociative nature of the PES cut along the O-H stretching coordinate, confirmed by the broad feature in the absorption spectrum. In contrast, more energetic excitations lead mainly to C-O bond breaking. Again, analysis of the diabatic surfaces offers a clear explanation in terms of the nature of the electronic states involved and the coupling between them. The type of calculations presented, as well as the subsequent analysis of the results, should be seen as a general workflow for the modeling of photochemical reactions.

Quantum dynamics simulations are becoming a powerful tool for understanding photo-excited molecules. Their poor scaling, however, means that it is hard to study molecules with more than a few atoms accurately, and a major challenge at the moment is the inclusion of the molecular environment. Here, we present a proof of principle for a way to break the two bottlenecks preventing large but accurate simulations. First, the problem of providing the potential energy surfaces for a general system is addressed by parameterizing a standard force field to reproduce the potential surfaces of the molecule's excited-states, including the all-important vibronic coupling. While not shown here, this would trivially enable the use of an explicit solvent. Second, to help the scaling of the nuclear dynamics propagation, a hierarchy of approximations is introduced to the variational multi-configurational Gaussian method that retains the variational quantum wavepacket description of the key quantum degrees of freedom and uses classical trajectories for the remaining in a quantum mechanics/molecular mechanics like approach. The method is referred to as force field quantum dynamics (FF-QD), and a two-state ππ*/nπ* model of uracil, excited to its lowest bright ππ* state, is used as a test case., (© 2024 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).)