We formulate a general hybrid quantum-classical technique to describe the interaction of diatomic molecules with XUV pulses. We demonstrate the accuracy of our model in the context of the interaction of the
O2 molecule with an XUV pulse with photon energy ranging from 20–42 eV. We account for the electronic structure and electron ionization quantum mechanically employing accurate molecular continuum wave functions. We account for the motion of the nuclei using classical equations of motion. However, the force of the nuclei is computed by obtaining accurate potential-energy curves of
O2 up to
O22+, relevant to the 20–42 eV photon-energy range, using advanced quantum-chemistry techniques. We find the dissociation limits of these states and the resulting atomic fragments and employ the velocity-Verlet algorithm to compute the velocities of these fragments. We incorporate both electron ionization and nuclear motion in a stochastic Monte Carlo simulation and identify the ionization and dissociation pathways when
O2 interacts with an XUV pulse. Focusing on the
O++O+ dissociation pathway, we obtain the kinetic-energy release (KER) distributions of the atomic fragments and find very good agreement with experimental results. Also, we explain the main features of the KER in terms of ionization sequences consisting of two sequential single-photon absorptions resulting in different
O+ and
O2+ electronic state configurations involved in the two transitions. Published by the American Physical Society 2025