RESUMEN

An analytical formulation for the geometrical derivatives of excitation energies within the time-dependent density-functional tight-binding (TD-DFTB) method is presented. The derivation is based on the auxiliary functional approach proposed in [Furche and Ahlrichs, J Chem Phys 2002, 117, 7433]. To validate the quality of the potential energy surfaces provided by the method, adiabatic excitation energies, excited state geometries, and harmonic vibrational frequencies were calculated for a test set of molecules in excited states of different symmetry and multiplicity. According to the results, the TD-DFTB scheme surpasses the performance of configuration interaction singles and the random phase approximation but has a lower quality than ab initio time-dependent density-functional theory. As a consequence of the special form of the approximations made in TD-DFTB, the scaling exponent of the method can be reduced to three, similar to the ground state. The low scaling prefactor and the satisfactory accuracy of the method makes TD-DFTB especially suitable for molecular dynamics simulations of dozens of atoms as well as for the computation of luminescence spectra of systems containing hundreds of atoms.

Asunto(s)

RESUMEN

Rhodopsins can modulate the optical properties of their chromophores over a wide range of wavelengths. The mechanism for this spectral tuning is based on the response of the retinal chromophore to external stress and the interaction with the charged, polar, and polarizable amino acids of the protein environment and is connected to its large change in dipole moment upon excitation, its large electronic polarizability, and its structural flexibility. In this work, we investigate the accuracy of computational approaches for modeling changes in absorption energies with respect to changes in geometry and applied external electric fields. We illustrate the high sensitivity of absorption energies on the ground-state structure of retinal, which varies significantly with the computational method used for geometry optimization. The response to external fields, in particular to point charges which model the protein environment in combined quantum mechanical/molecular mechanical (QM/MM) applications, is a crucial feature, which is not properly represented by previously used methods, such as time-dependent density functional theory (TDDFT), complete active space self-consistent field (CASSCF), and Hartree-Fock (HF) or semiempirical configuration interaction singles (CIS). This is discussed in detail for bacteriorhodopsin (bR), a protein which blue-shifts retinal gas-phase excitation energy by about 0.5 eV. As a result of this study, we propose a procedure which combines structure optimization or molecular dynamics simulation using DFT methods with a semiempirical or ab initio multireference configuration interaction treatment of the excitation energies. Using a conventional QM/MM point charge representation of the protein environment, we obtain an absorption energy for bR of 2.34 eV. This result is already close to the experimental value of 2.18 eV, even without considering the effects of protein polarization, differential dispersion, and conformational sampling.

Asunto(s)

RESUMEN

This work investigates the capability of time-dependent density functional response theory to describe excited state potential energy surfaces of conjugated organic molecules. Applications to linear polyenes, aromatic systems, and the protonated Schiff base of retinal demonstrate the scope of currently used exchange-correlation functionals as local, adiabatic approximations to time-dependent Kohn-Sham theory. The results are compared to experimental and ab initio data of various kinds to attain a critical analysis of common problems concerning charge transfer and long range (nondynamic) correlation effects. This analysis goes beyond a local investigation of electronic properties and incorporates a global view of the excited state potential energy surfaces.

Asunto(s)