Femtosecond and attosecond X-ray transient absorption experiments are becoming increasingly sophisticated tools for probing nuclear dynamics. In this work, we explore and develop theoretical tools needed for interpretation of such spectra,in order to characterize the vibrational coherences that result from ionizing a molecule in a strong IR field. Ab initio data for F2 is combined with simulations of nuclear dynamics, in order to simulate time-resolved X-ray absorption spectra for vibrational wavepackets after coherent ionization at 0 K and at finite temperature. Dihalogens pose rather difficult electronic structure problems, and the issues encountered in this work will be reflective of those encountered with any core–valence excitation simulation when a bond is breaking. The simulations reveal a strong dependence of the X-ray absorption maximum on the locations of the vibrational wave packets. A Fourier transform of the simulated signal shows features at the overtone frequencies of both the neutral and the cation, which reflect spatial interferences of the vibrational eigenstates. This provides a direct path for implementing ultrafast X-ray spectroscopic methods to visualize coherent nuclear dynamics.

A summary of the technical advances that are incorporated in the fourth major release of the Q-Chem quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and open-shell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Møller-Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Some specific chemical examples include calculations on the strongly correlated Cr2 dimer, exploring zeolite-catalysed ethane dehydrogenation, energy decomposition analysis of a charged ter-molecular complex arising from glycerol photoionisation, and natural transition orbitals for a Frenkel exciton state in a nine-unit model of a self-assembling nanotube.