HORTENSIA, a program package for the simulation of nonadiabatic autoionization dynamics in molecules

TL;DR

HORTENSIA is a computational package that simulates ultrafast autoionization dynamics in molecular anions using nonadiabatic surface-hopping, enabling detailed analysis of electron ejection processes and molecular configurations.

Contribution

This work introduces HORTENSIA, a novel program integrating nonadiabatic surface-hopping with quantum chemical calculations for simulating autoionization in molecules.

Findings

Successfully simulated time- and angle-resolved autoionization spectra

Provided insights into molecular geometries influencing autoionization

Enabled analysis of electron ejection dynamics in molecular anions

Abstract

We present a program package for the simulation of ultrafast vibration-induced autoionization dynamics in molecular anions in the manifold of the adiabatic anionic states and the discretized ionization continuum. This program, called HORTENSIA ($\underline{Ho}$pping $\underline{r}$eal-time $\underline{t}$rajectories for $\underline{e}$lectron-ejection by $\underline{n}$onadiabatic $\underline{s}$elf-$\underline{i}$onization in $\underline{a}$nions), is based on the nonadiabatic surface-hopping methodology, wherein nuclei are propagated as an ensemble along classical trajectories in the quantum-mechanical potential created by the electronic density of the molecular system. The electronic Schr\"odinger equation is numerically integrated along the trajectory, providing the time evolution of electronic state coefficients, from which switching probabilities into discrete electronic states are determined. In the case of a discretized continuum state, this hopping event is interpreted as the ejection on an electron. The derived diabatic and nonadiabatic couplings in the time-dependent electronic Schr\"odinger equation are calculated from anionic and neutral wavefunctions obtained from quantum chemical calculations with commercially available program packages interfaced with our program. Based on this methodology, we demonstrate the simulation of autoionization electron kinetic energy spectra that are both time- and angle-resolved. In addition, the program yields data that can be interpreted easily with respect to geometric characteristics such as bonding distances and angles, which facilitates the detection of molecular configurations important for the autoionization process. Moreover, useful extensions are included, namely generation tools for initial conditions and input files as well as for the evaluation of output files both through console commands and a graphical user interface.