Molecular properties of excited electronic state: Formalism, implementation, and applications of analytical second energy derivatives within the framework of the timedependent density functional theory/molecular mechanics
Abstract
This work extends our previous works [J. Liu and W. Z. Liang, J. Chem. Phys. 135, 014113 (2011); J. Liu and W. Z. Liang, J. Chem. Phys. 135, 184111 (2011)] on analytical excitedstate energy Hessian within the framework of timedependent density functional theory (TDDFT) to couple with molecular mechanics (MM). The formalism, implementation, and applications of analytical first and second energy derivatives of TDDFT/MM excited state with respect to the nuclear and electric perturbations are presented. Their performances are demonstrated by the calculations of adiabatic excitation energies, and excitedstate geometries, harmonic vibrational frequencies, and infrared intensities for a number of benchmark systems. The consistent results with the full quantum mechanical method and other hybrid theoretical methods indicate the reliability of the current numerical implementation of developed algorithms. The computational accuracy and efficiency of the current analytical approach are also checked and the computational efficient strategies are suggested to speed up the calculations of complex systems with many MM degrees of freedom. Finally, we apply the current analytical approach in TDDFT/MM to a realistic system, a red fluorescent protein chromophore together with part of its nearby protein matrix. The calculated results indicate that the rearrangement of the hydrogen bond interactions betweenmore »
 Authors:
 State Key Laboratory of Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Province Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005 (China)
 Hefei National Laboratory for Physical Sciences at Microscale, University of Science and Technology of China, Hefei 230026 (China)
 Publication Date:
 OSTI Identifier:
 22253548
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Chemical Physics; Journal Volume: 140; Journal Issue: 18; Other Information: (c) 2014 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 37 INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL CHEMISTRY; ACCURACY; DEGREES OF FREEDOM; DENSITY FUNCTIONAL METHOD; EFFICIENCY; EXCITATION; EXCITED STATES; FLUORESCENCE; HYBRIDIZATION; INTERACTIONS; PERTURBATION THEORY; PROTEINS; TIME DEPENDENCE
Citation Formats
Zeng, Qiao, Liang, WanZhen, Email: liangwz@xmu.edu.cn, and Liu, Jie. Molecular properties of excited electronic state: Formalism, implementation, and applications of analytical second energy derivatives within the framework of the timedependent density functional theory/molecular mechanics. United States: N. p., 2014.
Web. doi:10.1063/1.4863563.
Zeng, Qiao, Liang, WanZhen, Email: liangwz@xmu.edu.cn, & Liu, Jie. Molecular properties of excited electronic state: Formalism, implementation, and applications of analytical second energy derivatives within the framework of the timedependent density functional theory/molecular mechanics. United States. doi:10.1063/1.4863563.
Zeng, Qiao, Liang, WanZhen, Email: liangwz@xmu.edu.cn, and Liu, Jie. 2014.
"Molecular properties of excited electronic state: Formalism, implementation, and applications of analytical second energy derivatives within the framework of the timedependent density functional theory/molecular mechanics". United States.
doi:10.1063/1.4863563.
@article{osti_22253548,
title = {Molecular properties of excited electronic state: Formalism, implementation, and applications of analytical second energy derivatives within the framework of the timedependent density functional theory/molecular mechanics},
author = {Zeng, Qiao and Liang, WanZhen, Email: liangwz@xmu.edu.cn and Liu, Jie},
abstractNote = {This work extends our previous works [J. Liu and W. Z. Liang, J. Chem. Phys. 135, 014113 (2011); J. Liu and W. Z. Liang, J. Chem. Phys. 135, 184111 (2011)] on analytical excitedstate energy Hessian within the framework of timedependent density functional theory (TDDFT) to couple with molecular mechanics (MM). The formalism, implementation, and applications of analytical first and second energy derivatives of TDDFT/MM excited state with respect to the nuclear and electric perturbations are presented. Their performances are demonstrated by the calculations of adiabatic excitation energies, and excitedstate geometries, harmonic vibrational frequencies, and infrared intensities for a number of benchmark systems. The consistent results with the full quantum mechanical method and other hybrid theoretical methods indicate the reliability of the current numerical implementation of developed algorithms. The computational accuracy and efficiency of the current analytical approach are also checked and the computational efficient strategies are suggested to speed up the calculations of complex systems with many MM degrees of freedom. Finally, we apply the current analytical approach in TDDFT/MM to a realistic system, a red fluorescent protein chromophore together with part of its nearby protein matrix. The calculated results indicate that the rearrangement of the hydrogen bond interactions between the chromophore and the protein matrix is responsible for the large Stokes shift.},
doi = {10.1063/1.4863563},
journal = {Journal of Chemical Physics},
number = 18,
volume = 140,
place = {United States},
year = 2014,
month = 5
}

A parallel implementation of analytical timedependent density functional theory gradients is presented for the quantum chemistry program NWChem. The implementation is based on the Lagrangian approach developed by Furche and Ahlrichs. To validate our implementation, we first calculate the Stokes shifts for a range of organic dye molecules using a diverse set of exchangecorrelation functionals (traditional density functionals, global hybrids, and rangeseparated hybrids) followed by simulations of the onephoton absorption and resonance Raman scattering spectrum of the phenoxyl radical, the wellstudied dye molecule rhodamine 6G, and a molecular host–guest complex (TTFcCBPQT ^{4+}). The study of organic dye molecules illustrates thatmore »

Excitedstate forces within timedependent densityfunctional theory: A frequencydomain approach
We introduce a frequencydomain formalism for computing excitedstate forces using linear response theory within the framework of timedependent densityfunctional theory. We present an implementation of the formalism within a realspace computational framework. We demonstrate the validity and usefulness of the approach by comparing its results to those obtained with a Green's function approach using the GW+BetheSalpeterequation method, for the simple case of the CO molecule. We expect the method to be advantageous for computing excitedstate dynamics in large systems. 
Analytical derivatives of the individual state energies in ensemble density functional theory method. I. General formalism
The stateaveraged (SA) spin restricted ensemble referenced KohnSham (REKS) method and its state interaction (SI) extension, SISAREKS, enable one to describe correctly the shape of the ground and excited potential energy surfaces of molecules undergoing bond breaking/bond formation reactions including features such as conical intersections crucial for theoretical modeling of nonadiabatic reactions. Until recently, application of the SAREKS and SISAREKS methods to modeling the dynamics of such reactions was obstructed due to the lack of the analytical energy derivatives. Here, the analytical derivatives of the individual SAREKS and SISAREKS energies are derived. The final analytic gradient expressions are formulated entirelymore » 
Configuration interaction singles, timedependent HartreeFock, and timedependent density functional theory for the electronic excited states of extended systems
A general formalism for timedependent linear response theory is presented within the framework of linearcombinationofatomicorbital crystalline orbital theory for the electronic excited states of infinite onedimensional lattices (polymers). The formalism encompasses those of timedependent HartreeFock theory (TDHF), timedependent density functional theory (TDDFT), and configuration interaction singles theory (CIS) (as the TammDancoff approximation to TDHF) as particular cases. These singleexcitation theories are implemented by using a trialvector algorithm, such that the atomicorbitalbased twoelectron integrals are recomputed as needed and the transformation of these integrals from the atomicorbital basis to the crystallineorbital basis is avoided. Convergence of the calculated excitation energies withmore » 
Analytic energy gradient of excited electronic state within TDDFT/MMpol framework: Benchmark tests and parallel implementation
The timedependent density functional theory (TDDFT) has become the most popular method to calculate the electronic excitation energies, describe the excitedstate properties, and perform the excitedstate geometric optimization of medium and largesize molecules due to the implementation of analytic excitedstate energy gradient and Hessian in many electronic structure software packages. To describe the molecules in condensed phase, one usually adopts the computationally efficient hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) models. Here, we extend our previous work on the energy gradient of TDDFT/MM excited state to account for the mutual polarization effects between QM and MM regions, which is believed to holdmore »