  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.
