Harmonic Models in Cartesian and Internal Coordinates to Simulate the Absorption Spectra of Carotenoids at Finite Temperatures

• Published in: Journal of chemical theory and computation Vol. 9; no. 11; pp. 4947 - 4958
• Main Authors: Cerezo, Javier, Zúñiga, José, Requena, Alberto, Ávila Ferrer, Francisco J, Santoro, Fabrizio
• Format: Journal Article
• Language: English
• Published: United States American Chemical Society 12-11-2013
• ISSN: 1549-9618; 1549-9626
• DOI: 10.1021/ct4005849

When large structural displacements take place between the ground state (GS) and excited state (ES) minima of polyatomic molecules, the choice of a proper set of coordinates can be crucial for a reliable simulation of the vibrationally resolved absorption spectrum. In this work, we study two carotenoids that undergo structural displacements from GS to ES minima of different magnitude, from small displacements for violaxanthin to rather large ones for β-carotene isomers. Their finite-temperature (77 and 300 K) spectra are simulated at the harmonic level, including Duschinsky effect, by time-dependent (TD) and time-independent (TI) approaches, using (TD)DFT computed potential energy surfaces (PES). We adopted two approaches to construct the harmonic PES, the Adiabatic (AH) and Vertical Hessian (VH) models and, for AH, two reference coordinate frames: Cartesian and valence internal coordinates. Our results show that when large displacements take place, Cartesian coordinates dramatically fail to describe curvilinear displacements and to account for the Duschinsky matrix, preventing a realistic simulation of the spectra within the AH model, where the GS and ES PESs are quadratically expanded around their own equilibrium geometry. In contrast, internal coordinates largely amend such deficiencies and deliver reasonable spectral widths. As expected, both coordinate frames give similar results when small displacements occur. The good agreement between VH and experimental line shapes indicates that VH model, in which GS and ES normal modes are both evaluated at the GS equilibrium geometry, is a good alternative to deal with systems exhibiting large displacements. The use of this model can be, however, problematic when imaginary frequencies arise. The extent of the nonorthogonality of the Dushinsky matrix in internal coordinates and its correlation with the magnitude of the displacement of the GS and ES geometries is analyzed in detail.