Full Time-dependent Density Functional Theory Research Articles

Time-dependent density functional theory based upon the fragment molecular orbital method

Time-dependent density functional theory (TDDFT) was combined with the two-body fragment molecular orbital method (FMO2). In this FMO2-TDDFT scheme, the system is divided into fragments, and the electron density for fragments is determined self-consistently. Consequently, only one main fragment of interest and several fragment pairs including it are calculated by TDDFT. To demonstrate the accuracy of FMO2-TDDFT, we computed several low-lying singlet and triplet excited states of solvated phenol and polyalanine using our method and the standard TDDFT for the full system. The BLYP functional with the long-range correction (LC-BLYP) was employed with the 6-31G(*) basis set (some tests were also performed with 6-311G(*), as well as with B3LYP and time-dependent Hartree-Fock). Typically, FMO2-TDDFT reproduced the full TDDFT excitation energies within 0.1 eV, and for one excited state the error was about 0.2 eV. Beside the accurate reproduction of the TDDFT excitation energies, we also automatically get an excitation energy decomposition analysis, which provides the contributions of individual fragments. Finally, the efficiency of our approach was exemplified on the LC-BLYP6-31G(*) calculation of the lowest singlet excitation of the photoactive yellow protein which consists of 1931 atoms, and the obtained value of 3.1 eV is in agreement with the experimental value of 2.8 eV.

Time-dependent density functional theory with the multilayer fragment molecular orbital method

Time-dependent density functional theory (TDDFT) is combined with the multilayer fragment molecular orbital (FMO) method. The excitation energies of the system are estimated from those of its small active part, whereas the remaining part is described by its electron density distribution and the corresponding Coulomb field. The efficiency of this approach is illustrated on the lowest excitations of several solvated systems and phenylalanine–(alanine) n . The full TDDFT excitation energies calculated are reproduced within 0.01 eV with LC-BLYP and the 6-31G ∗ basis set.

Mechanism of hydrogen bond assisted cis-trans isomerization in the excited state of 2-naphthol

Stationary points on the ground and first two excited state potential energy surfaces of 2-naphthol and the hydrogen bonded complex 2-naphtholNH(3) have been located using an approximate method (time dependent density functional based tight-binding) based on time dependent density functional theory (TDDFT). For the obtained extrema, single point calculations at the full TDDFT level (PW916-311+G(**)) were performed for validation purposes. The energetical ordering of the cis and trans forms of 2-naphthol in the different electronic states is found to be in good agreement with available experimental data. Interestingly, also the optical spectrum is correctly predicted by TDDFT in contrast to calculations on unsubstituted acenes in the literature. Besides the stationary points, also the minimum energy paths connecting them are calculated using an adaptive nudged elastic band approach in order to estimate isomerization barriers. While the increase of barrier height going from the ground to the first excited state is already documented, the authors find as a new result a significant barrier reduction upon hydrogen bond formation in the S(2). The observed effect can contribute to the explanation of recent experiments on infrared laser induced cis-trans isomerization in 2-naphthol.