Search papers, labs, and topics across Lattice.
This paper explores the use of orbital-optimized density functional theory (DFT) to calculate the valence and Rydberg excited states of CO2, focusing on bond dissociation curves. The study compares the performance of several density functionals with varying degrees of exchange, using both real and complex-valued orbitals optimized via saddle point searches. The results show that PBE functional with complex orbitals achieves excitation energies within 0.3 eV of multireference configuration interaction, and dissociation curves closely match high-level calculations, demonstrating the method's accuracy and efficiency.
Orbital-optimized DFT can now accurately model CO2 excited states and dissociation curves at a fraction of the computational cost of traditional multireference methods.
Calculations of the lowest valence {\pi}* as well as the 3s and higher energy 3p{\sigma} Rydberg excited states of the CO2 molecule are carried out using density functionals with variational optimization of the orbitals, an approach involving relatively little computational effort. Five functionals with varying degree of exchange are used in combination with real or complex-valued orbitals that are optimized by finding saddle points on the electronic energy surface corresponding to the excited states. When the PBE functional is used in combination with complex orbitals, the calculated excitation energy is found to be within 0.3 eV of multireference configuration interaction reference values, and the results are further improved with hybrid functionals. In contrast, linear-response time-dependent density functional theory calculations give errors up to 1.9 eV for the most diffuse 3p{\sigma} excitation and exhibit stronger dependence on both the excitation character and the functional used. Calculated C-O dissociation curves using the PBE functional and the orbital-optimized approach compare remarkably well with the reported multireference configuration interaction and equation-of-motion coupled-cluster singles and doubles calculations. Thanks to the low computational cost, these results demonstrate that orbital-optimized density functional calculations can be a promising route for modelling photorelaxation in condensed-phase CO2, for example in the context of interstellar cosmic-ray radiation driven process involving high-energy Rydberg states.
