Ionization of pyridine: Interplay of orbital relaxation and electron correlation

The valence shell ionization spectrum of pyridine was studied using the third-order algebraic-diagrammatic construction approximation scheme for the one-particle Green's function and the outer-valence Green's function method. The results were used to interpret angle resolved photoelectron spectra recorded with synchrotron radiation in the photon energy range of 17-120 eV. The lowest four states of the pyridine radical cation, namely, ²A₂(1a₂^-1), ²A₁(7a₁^-1), ²B₁(2b₁^-1), and ²B₂(5b₂^-1), were studied in detail using various high-level electronic structure calculation methods. The vertical ionization energies were established using the equation-of-motion coupled-cluster approach with single, double, and triple excitations (EOM-IP-CCSDT) and the complete basis set extrapolation technique. Further interpretation of the electronic structure results was accomplished using Dyson orbitals, electron density difference plots, and a second-order perturbation theory treatment for the relaxation energy. Strong orbital relaxation and electron correlation effects were shown to accompany ionization of the 7a₁ orbital, which formally represents the nonbonding σ-type nitrogen lone-pair (nσ) orbital. The theoretical work establishes the important roles of the π-system (π-π* excitations) in the screening of the nσ-hole and of the relaxation of the molecular orbitals in the formation of the 7a₁(nσ)^-1 state. Equilibrium geometric parameters were computed using the MP2 (second-order Møller-Plesset perturbation theory) and CCSD methods, and the harmonic vibrational frequencies were obtained at the MP2 level of theory for the lowest three cation states. The results were used to estimate the adiabatic 0-0 ionization energies, which were then compared to the available experimental and theoretical data. Photoelectron anisotropy parameters and photoionization partial cross sections, derived from the experimental spectra, were compared to predictions obtained with the continuum multiple scattering approach.