A protocol for the investigation of the IVR problem: The isomerization of HONO as a case of study

The conformational isomerization of nitrous acid (HONO) promoted by excitation of the $\nu_{OH}$ or $\nu_{N=O}$ stretching normal coordinates is the first observed case of an infrared-induced photochemical reaction. The energy captured by the excited normal modes is redistributed into a highly excited vibrational level of the $\tau_{OH}$ torsion normal coordinate, which is the isomerization reaction coordinate. Herein, we present simple numerical methods to qualitatively investigate the coupling between the normal coordinates and the possible gateways for vibrational energy redistribution leading to the isomerization process. Our methodology involves the generation of the relevant 2D Potential Energy Surface (PES), by spanning the reaction coordinate and one of the $3N-7$ projected normal coordinates along the Intrinsic Reaction Coordinate (IRC). Once the PES has been obtained, the time-independent wavefunctions are calculated using the standard Discrete Variable Representation (DVR) approach. The reaction barrier is investigated using the Interacting Quantum Atoms (IQA) decomposition scheme, evidencing an important contribution from the exchange-correlation energy to the isomerization. Coupling between normal coordinates indicates preferential normal modes to redistribute the vibrational energy. Thermal activated 1D tunneling rates were found to be negligible.