Modeling Intermolecular Coulombic Decay with Non-Hermitian Real-Time Time-Dependent Density Functional Theory

In this work, we investigate the capability of using real-time time-dependent density functional theory (RT-TDDFT) in conjunction with a complex absorbing potential (CAP) to simulate the intermolecular Coulombic decay (ICD) processes following the ionization of an inner-valence electron. We examine the ICD dynamics in a series of noncovalent bonded dimer systems, including hydrogen-bonded and purely van der Waals (VdW)-bonded systems. In comparison to previous work, we show that RT-TDDFT simulations with a CAP correctly capture the ICD phenomenon in systems exhibiting a stronger binding energy. The calculated time scales for ICD of the studied systems are in the range of 5–50 fs, in agreement with previous studies. However, there is a breakdown in the accuracy of the methodology for the pure VdW-bonded systems. Overall, the presented RT-TDDFT/CAP methodology provides a powerful tool for differentiating between competing electronic relaxation pathways following inner-valence or core ionization without necessitating any a priori assumptions.