No Arabic abstract
The computational investigation of photochemical processes often entails the calculation of excited state geometries, energies, and energy gradients. The nuclear-electronic orbital (NEO) approach treats specified nuclei, typically protons, quantum mechanically on the same level as the electrons, thereby including the associated nuclear quantum effects and non-Born-Oppenheimer behavior into quantum chemistry calculations. The multicomponent density functional theory (NEO-DFT) and time-dependent DFT (NEO-TDDFT) methods allow efficient calculations of ground and excited states, respectively. Herein, the analytical gradients are derived and implemented for the NEO-TDDFT method and the associated Tamm-Dancoff approximation (NEO-TDA). The programmable equations for these analytical gradients, as well as the NEO-DFT analytical Hessian, are provided. The NEO approach includes the anharmonic zero-point energy and density delocalization associated with the quantum protons, as well as vibronic mixing, in geometry optimizations and energy calculations of ground and excited states. The harmonic zero-point energy associated with the other nuclei can be computed via the NEO Hessian. This approach is used to compute the 0-0 adiabatic excitation energies for a set of nine small molecules with all protons quantized, exhibiting slight improvement over the conventional electronic approach. Geometry optimizations of two excited state intramolecular proton transfer systems are performed with one and two quantized protons, respectively. The NEO calculations for these systems produce electronically excited state geometries with stronger intramolecular hydrogen bonds and similar relative stabilities compared to conventional electronic methods. This work provides the foundation for nonadiabatic dynamics simulations of fundamental processes such as photoinduced proton transfer and proton-coupled electron transfer.
Time-dependent orbital-free density functional theory (TD-OFDFT) is an efficient ab-initio method for calculating the electronic dynamics of large systems. In comparison to standard TD-DFT, it computes only a single electronic state regardless of system size, but it requires an additional time-dependent Pauli potential term. We propose a nonadiabatic and nonlocal Pauli potential whose main ingredients are the time-dependent particle and current densities. Our calculations of the optical spectra of metallic and semiconductor clusters indicate that nonlocal and nonadiabatic TD-OFDFT performs accurately for metallic systems and semiquantitatively for semiconductors. This work opens the door to wide applicability of TD-OFDFT for nonequilibrium electron and electron-nuclear dynamics of materials.
A set of density functionals coming from different rungs on Jacobs ladder are employed to evaluate the electronic excited states of three Ru(II) complexes. While most studies on the performance of density functionals compare the vertical excitation energies, in this work we focus on the energy gaps between the electronic excited states, of the same and different multiplicity. Excited state energy gaps are important for example to determine radiationless transition probabilities. Besides energies, a functional should deliver the correct state character and state ordering. Therefore, wavefunction overlaps are introduced to systematically evaluate the effect of different functionals on the character of the excited states. As a reference, the energies and state characters from multi-state second-order perturbation theory complete active space (MS-CASPT2) are used. In comparison to MS-CASPT2, it is found that while hybrid functionals provide better vertical excitation energies, pure functionals typically give more accurate excited state energy gaps. Pure functionals are also found to reproduce the state character and ordering in closer agreement to MS-CASPT2 than the hybrid functionals.
Time-dependent orbital-free DFT is an efficient method for calculating the dynamic properties of large scale quantum systems due to the low computational cost compared to standard time-dependent DFT. We formalize this method by mapping the real system of interacting fermions onto a fictitious system of non-interacting bosons. The dynamic Pauli potential and associated kernel emerge as key ingredients of time-tependent orbital-free DFT. Using the uniform electron gas as a model system, we derive an approximate frequency-dependent Pauli kernel. Pilot calculations suggest that space nonlocality is a key feature for this kernel. Nonlocal terms arise already in the second order expansion with respect to unitless frequency and reciprocal space variable ($frac{omega}{q, k_F}$ and $frac{q}{2, k_F}$, respectively). Given the encouraging performance of the proposed kernel, we expect it will lead to more accurate orbital-free DFT simulations of nanoscale systems out of equilibrium. Additionally, the proposed path to formulate nonadiabatic Pauli kernels presents several avenues for further improvements which can be exploited in future works to improve the results.
A very specific ensemble of ground and excited states is shown to yield an exact formula for any excitation energy as a simple correction to the energy difference between orbitals of the Kohn-Sham ground state. This alternative scheme avoids either the need to calculate many unoccupied levels as in time-dependent density functional theory (TDDFT) or the need for many self-consistent ensemble calculations. The symmetry-eigenstate Hartree-exchange (SEHX) approximation yields results comparable to standard TDDFT for atoms. With this formalism, SEHX yields approximate double-excitations, which are missed by adiabatic TDDFT.
The time-dependent density functional theory (TDDFT) has been broadly used to investigate the excited-state properties of various molecular systems. However, the current TDDFT heavily relies on outcomes from the corresponding ground-state density functional theory (DFT) calculations which may be prone to errors due to the lack of proper treatment in the non-dynamical correlation effects. Recently, thermally-assisted-occupation density functional theory (TAO-DFT) [J.-D. Chai, textit{J. Chem. Phys.} textbf{136}, 154104 (2012)], a DFT with fractional orbital occupations, was proposed, explicitly incorporating the non-dynamical correlation effects in the ground-state calculations with low computational complexity. In this work, we develop time-dependent (TD) TAO-DFT, which is a time-dependent, linear-response theory for excited states within the framework of TAO-DFT. With tests on the excited states of H$_{2}$, the first triplet excited state ($1^3Sigma_u^+$) was described well, with non-imaginary excitation energies. TDTAO-DFT also yields zero singlet-triplet gap in the dissociation limit, for the ground singlet ($1^1Sigma_g^+$) and the first triplet state ($1^3Sigma_u^+$). In addition, as compared to traditional TDDFT, the overall excited-state potential energy surfaces obtained from TDTAO-DFT are generally improved and better agree with results from the equation-of-motion coupled-cluster singles and doubles (EOM-CCSD).