Vibronic coupling and double excitations in linear response time-dependent density functional calculations: Dipole-allowed states of N2

Time dependent vibrational electronic coupled cluster (VECC) theory for non-adiabatic nuclear dynamics

A time-dependent vibrational electronic coupled-cluster (VECC) approach is proposed to simulate photo-electron/UV-VIS absorption spectra as well as time-dependent properties for non-adiabatic vibronic models, going beyond the Born-Oppenheimer approximation. A detailed derivation of the equations of motion and a motivation for the ansatz are presented. The VECC method employs second-quantized bosonic construction operators and a mixed linear and exponential ansatz to form a compact representation of the time-dependent wave-function. Importantly, the method does not require a basis set, has only a few user-defined inputs, and has a classical (polynomial) scaling with respect to the number of degrees of freedom (of the vibronic model), resulting in a favorable computational cost. In benchmark applications to small models and molecules, the VECC method provides accurate results compared to multi-configurational time-dependent Hartree calculations when predicting short-time dynamical properties (i.e., photo-electron/UV-VIS absorption spectra) for non-adiabatic vibronic models. To illustrate the capabilities, the VECC method is also successfully applied to a large vibronic model for hexahelicene with 14 electronic states and 63 normal modes, developed in the group by Aranda and Santoro [J. Chem. Theory Comput. 17, 1691, (2021)].

Doubly Tuned Exchange-Correlation Functionals for Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory

It is demonstrated that significant accuracy improvements in MRSF-TDDFT can be achieved by introducing two different exchange-correlation (XC) functionals for the reference Kohn-Sham DFT and the response part of the calculations, respectively. Accordingly, two new XC functionals of doubly tuned Coulomb attenuated method-vertical excitation energy (DTCAM-VEE) and DTCAM-AEE were developed on the basis of the "adaptive exact exchange (AEE)" concept in the framework of the Coulomb-attenuating XC functionals. The values by DTCAM-VEE are in excellent agreement with those of Thiel's set [mean absolute errors (MAEs) and the interquartile range (IQR) values of 0.218 and 0.327 eV, respectively]. On the other hand, DTCAM-AEE faithfully reproduced the qualitative aspects of conical intersections (CIs) of trans-butadiene and thymine and the nonadiabatic molecular dynamics (NAMD) simulations on thymine. The latter functional also remarkably exhibited the exact 1/R asymptotic behavior of the charge-transfer state of an ethylene-tetrafluoroethylene dimer and the accurate potential energy surfaces (PESs) along the two torsional angles of retinal protonated Schiff base model with six double bonds (rPSB6). Overall, DTCAM-AEE generally performs well, as its MAE (0.237) and IQR (0.41 eV) are much improved as compared to BH&HLYP. The current idea can also be applied to other XC functionals as well as other variants of linear response theories, opening a new way of developing XC functionals.

Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory: Multireference Advantages with the Practicality of Linear Response Theory

The density functional theory (DFT) and linear response (LR) time-dependent (TD)-DFT are of the utmost importance for routine computations. However, the single reference formulation of DFT suffers in the description of open-shell singlet systems such as diradicals and bond-breaking. LR-TDDFT, on the other hand, finds difficulties in the modeling of conical intersections, doubly excited states, and core-level excitations. In this Perspective, we demonstrate that many of these limitations can be overcome by recently developed mixed-reference (MR) spin-flip (SF)-TDDFT, providing an alternative yet accurate route for such challenging situations. Empowered by the practicality of the LR formalism, it is anticipated that MRSF-TDDFT can become one of the major workhorses for general routine tasks.

High-performance strategies for the recent MRSF-TDDFT in GAMESS

Multiple ERI (Electron Repulsion Integral) tensor contractions (METC) with several matrices are ubiquitous in quantum chemistry. In response theories, the contraction operation, rather than ERI computations, can be the major bottleneck, as its computational demands are proportional to the multiplicatively combined contributions of the number of excited states and the kernel pre-factors. This paper presents several high-performance strategies for METC. Optimal approaches involve either the data layout reformations of interim density and Fock matrices, the introduction of intermediate ERI quartet buffer, and loop-reordering optimization for a higher cache hit rate. The combined strategies remarkably improve the performance of the MRSF (mixed reference spin flip)-TDDFT (time-dependent density functional theory) by nearly 300%. The results of this study are not limited to the MRSF-TDDFT method and can be applied to other METC scenarios.

Computational Model for Electrochemical Surface-Enhanced Raman Scattering: Key Role of the Surface Charges and Synergy between Electromagnetic and Charge-Transfer Enhancement Mechanisms

We present a computational model for electrochemical surface-enhanced Raman scattering (EC-SERS). The surface excess of charge induced by the electrode potential (V_el) was introduced by applying an external electric field to a set of clusters [Ag_n]^q with (n, q) of (19, ±1) or (20, 0) on which a molecule adsorbs. Using DFT/TD-DFT calculations, these metal-molecule complexes were classified by the adsorbate partial charge, and the main V_el-dependent properties were simultaneously studied with the aid of vibronic resonance Raman computations, namely, changes on the vibrational wavenumbers, relative intensities, and enhancement factors (EFs) for all SERS mechanisms: chemical or nonresonant, and resonance Raman with bright states of the adsorbate, charge-transfer (CT) states, and plasmon-like excitations on the metal cluster. We selected two molecules to test our model, pyridine, for which V_el has a remarkable effect, and 9,10-bis((E)-2-(pyridin-4-yl)vinyl)anthracene, which is almost insensitive to the applied bias. The results nicely reproduced most of the experimental observations, while the limitations of our approach were critically evaluated. We detected that accounting explicitly for the surface charges is key for EC-SERS models and that the highest calculated EFs, up to 10⁷ to 10⁸, are obtained by interstate coupling of bright local excitations of the metal cluster and CT states. These results highlight the importance of nonadiabatic effects in SERS and the capabilities of EC-SERS as a technique with potential to study excited-state coupling by tuning the CT and plasmon-like states by manipulating V_el.

Surface Hopping Dynamics on Vibronic Coupling Models

The simulation of photoinduced non-adiabatic dynamics is of great relevance in many scientific disciplines, ranging from physics and materials science to chemistry and biology. Upon light irradiation, different relaxation processes take place in which electronic and nuclear motion are intimately coupled. These are best described by the time-dependent molecular Schrödinger equation, but its solution poses fundamental practical challenges to contemporary theoretical chemistry. Two widely used and complementary approaches to this problem are multiconfigurational time-dependent Hartree (MCTDH) and trajectory surface hopping (SH). MCTDH is an accurate fully quantum-mechanical technique but often is feasible only in reduced dimensionality, in combination with approximate vibronic coupling (VC) Hamiltonians, or both (i.e., reduced-dimensional VC potentials). In contrast, SH is a quantum–classical technique that neglects most nuclear quantum effects but allows nuclear dynamics in full dimensionality by calculating potential energy surfaces on the fly. If nuclear quantum effects do not play a central role and a linear VC (LVC) Hamiltonian is appropriate—e.g., for stiff molecules that generally keep their conformation in the excited state—then it seems advantageous to combine the efficient LVC and SH techniques. In this Account, we describe how surface hopping based on an LVC Hamiltonian (SH/LVC)—as recently implemented in the SHARC surface hopping package—can provide an economical and automated approach to simulate non-adiabatic dynamics. First, we illustrate the potential of SH/LVC in a number of showcases, including intersystem crossing in SO₂, intra-Rydberg dynamics in acetone, and several photophysical studies on large transition-metal complexes, which would be much more demanding or impossible to perform with other methods. While all of the applications provide very useful insights into light-induced phenomena, they also hint at difficulties faced by the SH/LVC methodology that need to be addressed in the future. Second, we contend that the SH/LVC approach can be useful to benchmark SH itself. By the use of the same (LVC) potentials as MCTDH calculations have employed for decades and by relying on the efficiency of SH/LVC, it is possible to directly compare multiple SH test calculations with a MCTDH reference and ponder the accuracy of various correction algorithms behind the SH methodology, such as decoherence corrections or momentum rescaling schemes. Third, we demonstrate how the efficiency of SH/LVC can also be exploited to identify essential nuclear and electronic degrees of freedom to be employed in more accurate MCTDH calculations. Lastly, we show that SH/LVC is able to advance the development of SH protocols that can describe nuclear dynamics including explicit laser fields—a very challenging endeavor for trajectory-based schemes. To end, this Account compiles the typical costs of contemporary SH simulations, evidencing the great advantages of using parametrized potentials. The LVC model is a sleeping beauty that, kissed by SH, is fueling the field of excited-state molecular dynamics. We hope that this Account will stimulate future research in this direction, leveraging the advantages of the SH/VC schemes to larger extents and extending their applicability to uncharted territories.

Fragment Diabatization Linear Vibronic Coupling Model for Quantum Dynamics of Multichromophoric Systems: Population of the Charge-Transfer State in the Photoexcited Guanine-Cytosine Pair

We introduce a method (FrD-LVC) based on a fragment diabatization (FrD) for the parametrization of a linear vibronic coupling (LVC) model suitable for studying the photophysics of multichromophore systems. In combination with effective quantum dynamics (QD) propagations with multilayer multiconfigurational time-dependent Hartree (ML-MCTDH), the FrD-LVC approach gives access to the study of the competition between intrachromophore decays, like those at conical intersections, and interchromophore processes, like exciton localization/delocalization and the involvement of charge-transfer (CT) states. We used FrD-LVC parametrized with time-dependent density functional theory (TD-DFT) calculations, adopting either CAM-B3LYP or ωB97X-D functionals, to study the ultrafast photoexcited QD of a guanine-cytosine (GC) hydrogen-bonded pair, within a Watson-Crick arrangement, considering up to 12 coupled diabatic electronic states and the effect of all of the 99 vibrational coordinates. The bright excited states localized on C and, especially, on G are predicted to be strongly coupled to the G → C CT state, which is efficiently and quickly populated after an excitation to any of the four lowest energy bright local excited states. Our QD simulations show that more than 80% of the excited population on G and ∼50% of that on C decay to this CT state in less than 50 fs. We investigate the role of vibronic effects in the population of the CT state and show that it depends mainly on its large reorganization energy so that it can occur even when it is significantly less stable than the bright states in the Franck-Condon region. At the same time, we document that the formation of the GC pair almost suppresses the involvement of dark nπ* excited states in the photoactivated dynamics.

Parameterization of a linear vibronic coupling model with multiconfigurational electronic structure methods to study the quantum dynamics of photoexcited pyrene

With this work, we present a protocol for the parameterization of a Linear Vibronic Coupling (LVC) Hamiltonian for quantum dynamics using highly accurate multiconfigurational electronic structure methods such as RASPT2/RASSCF, combined with a maximum-overlap diabatization technique. Our approach is fully portable and can be applied to many medium-size rigid molecules whose excited state dynamics requires a quantum description. We present our model and discuss the details of the electronic structure calculations needed for the parameterization, analyzing critical situations that could arise in the case of strongly interacting excited states. The protocol was applied to the simulation of the excited state dynamics of the pyrene molecule, starting from either the first or the second bright state (S₂ or S₅). The LVC model was benchmarked against state-of-the-art quantum mechanical calculations with optimizations and energy scans and turned out to be very accurate. The dynamics simulations, performed including all active normal coordinates with the multilayer multiconfigurational time-dependent Hartree method, show good agreement with the available experimental data, endorsing prediction of the excited state mechanism, especially for S₅, whose ultrafast deactivation mechanism was not yet clearly understood.

Ultrafast Dynamics of the Two Lowest Bright Excited States of Cytosine and 1-Methylcytosine: A Quantum Dynamical Study

The nonadiabatic quantum dynamics (QD) of cytosine and 1-methylcytosine in the gas phase is simulated for 250 fs after a photoexcitation to one of the first two bright states. The nuclear wavepacket is propagated on the coupled diabatic potential energy surfaces of the lowest seven excited states, including ππ*, nπ*, and Rydberg states along all the vibrational degrees of freedom. We focus in particular on the interplay between the bright and the dark nπ* states, not considering the decay to the ground electronic state. To run these simulations, we implemented an automatic general procedure to parametrize linear vibronic coupling (LVC) models with time-dependent density functional theory (DFT) computations and interfaced it with Gaussian package. The wavepacket was propagated with the multilayer version of the multiconfigurational time dependent Hartree method. Two different density functionals, PBE0 and CAM-B3LYP, which provide a different description of the relative stability of the lowest energy dark states, were used to parametrize the LVC Hamiltonian. Part of the photoexcited population on lowest HOMO-LUMO transition (π_Hπ_L^*) decays within less than 100 fs to a nπ* state which mainly involves a promotion of an electron from the oxygen lone pair to the LUMO (n_Oπ_L^*). The population of the second ππ* state decays almost completely, in <100 fs, not only to π_Hπ_L^* and to n_Oπ_L^* states but also to another nπ_L^* state involving the nitrogen lone pair. The efficiency of the adopted protocol allowed us to check the accuracy of the predictions by repeating the QD simulations with different LVC Hamiltonians parametrized either at the ground-state minimum or at stationary structures of different relevant excited states.

Mixed Quantum/Classical Method for Nonadiabatic Quantum Dynamics in Explicit Solvent Models: The ππ*/nπ* Decay of Thymine in Water as a Test Case

We present a novel mixed quantum classical dynamical method to include solvent effects on internal conversion (IC) processes. All the solute degrees of freedom are represented by a wavepacket moving according to nonadiabatic quantum dynamics, while the motion of an explicit solvent model is described by an ensemble of classical trajectories. The mutual coupling of the solute and solvent dynamics is included within a mean-field framework and the quantum and classical equations of motions are solved simultaneously. As a test case we apply our method to the ultrafast ππ* → nπ* decay of thymine in water. Solvent dynamical response modifies IC yield already on the 50 fs time scale. This effect is due to water librational motions that stabilize the most populated state. Pure static disorder, that is, the existence of different solvent configurations when photoexcitation takes place, also has a remarkable impact on the dynamics.

The adiabatic strictly-correlated-electrons functional: kernel and exact properties

We investigate a number of formal properties of the adiabatic strictly-correlated electrons (SCE) functional, relevant for time-dependent potentials and for kernels in linear response time-dependent density functional theory. Among the former, we focus on the compliance to constraints of exact many-body theories, such as the generalised translational invariance and the zero-force theorem. Within the latter, we derive an analytical expression for the adiabatic SCE Hartree exchange-correlation kernel in one dimensional systems, and we compute it numerically for a variety of model densities. We analyse the non-local features of this kernel, particularly the ones that are relevant in tackling problems where kernels derived from local or semi-local functionals are known to fail.

Combined spectroscopic and computational analysis of the vibrational properties of vitamin B12 in its Co3+, Co2+, and Co1+ oxidation states

While the geometric and electronic structures of vitamin B12 (cyanocobalamin, CNCbl) and its reduced derivatives Co(2+)cobalamin (Co(2+)Cbl) and Co(1+)cobalamin (Co(1+)Cbl(-)) are now reasonably well established, their vibrational properties, in particular their resonance Raman (rR) spectra, have remained quite poorly understood. The goal of this study was to establish definitive assignments of the corrin-based vibrational modes that dominate the rR spectra of vitamin B12 in its Co(3+), Co(2+), and Co(1+) oxidation states. rR spectra were collected for all three species with laser excitation in resonance with the most intense corrin-based π → π* transitions. These experimental data were used to validate the computed vibrational frequencies, eigenvector compositions, and relative rR intensities of the normal modes of interest as obtained by density functional theory (DFT) calculations. Importantly, the computational methodology employed in this study successfully reproduces the experimental observation that the frequencies and rR excitation profiles of the corrin-based vibrational modes vary significantly as a function of the cobalt oxidation state. Our DFT results suggest that this variation reflects large differences in the degree of mixing between the occupied Co 3d orbitals and empty corrin π* orbitals in CNCbl, Co(2+)Cbl, and Co(1+)Cbl(-). As a result, vibrations mainly involving stretching of conjugated C-C and C-N bonds oriented along one axis of the corrin ring may, in fact, couple to a perpendicularly polarized electronic transition. This unusual coupling between electronic transitions and vibrational motions of corrinoids greatly complicates an assignment of the corrin-based normal modes of vibrations on the basis of their rR excitation profiles.

Double excitations from modified Hartree Fock subsequent minimization scheme

Doubly excited states have nowadays become important in technological applications, e.g., in increasing the efficiency of solar cells and therefore, their description using ab initio methods is a great theoretical challenge as double excitations cannot be described by linear response theories based on a single Slater determinant. In the present work we extend our recently developed Hartree-Fock (HF) approximation for calculating singly excited states [M. Tassi, I. Theophilou, and S. Thanos, Int. J. Quantum Chem. 113, 690 (2013)] in order to allow for the calculation of doubly excited states. We describe the double excitation as two holes in the subspace spanned from the occupied HF orbitals and two particles in the subspace of virtual HF orbitals. A subsequent minimization of the energy results to the determination of the spin orbitals of both the holes and the particles in the occupied and virtual subspaces, respectively. We test our method, for various atoms, H2 and polyene molecules which are known to have excitations presenting a significant double excitation character. Importantly, our approach is computationally inexpensive.

The analog of Koopmans’ theorem for virtual Kohn–Sham orbital energies

An analog of Koopmans’ theorem is formulated for the energies, εa, of virtual Kohn–Sham (KS) molecular orbitals (MOs) from the requirement that the KS theory provides, in principle, not only the exact electron density, but also its exact response. The starting point is the Kohn–Sham analog of Koopmans’ theorem, relating the vertical ionization energies, Ii, to the energies, εi, of the occupied MOs ( Chong, D.P.; Gritsenko, O.V.; Baerends, E.J. J. Chem. Phys. 2002, 116, 1760 ). Combining this with the coupled-perturbed equations of time-dependent density functional theory (TDDFT), exact relations between the energies, εa, of virtual KS MOs and the excitation energies, ωia, and vertical ionization energies (VIPs), Ii, are obtained. In the small matrix approximation for the coupling matrix K of TDDFT, two limiting cases of these relations are considered. In the limit of a negligible matrix element, Kia,ia, the energy, εa, can be interpreted as (minus) the energy of ionization from the ?i → ?a excited state, εa ≈ –Ia, where –Ia is defined from the relation Ii = ωia + Ia. This relation breaks down in special cases, such as charge-transfer transitions and the HOMO–LUMO (highest occupied molecular orbital – lowest unoccupied molecular orbital) transition of a dissociating electron-pair bond (also of charge-transfer character). The present results highlight the important difference between virtual orbital energies in the Kohn–Sham model (εa ≈ –Ia) and in the Hartree–Fock model (εa ≈ –Aa). Kohn–Sham differences εa – εi approximate the excitation energy, ωia, while Hartree–Fock differences [Formula: see text] do not approximate excitation energies but approximate the difference of an ionization energy and an electron affinity, Ii – Aa.

Charge transfer, double and bond-breaking excitations with time-dependent density matrix functional theory

Time-dependent density functional theory (TDDFT) in its current adiabatic implementations exhibits three remarkable failures: (a) completely wrong behavior of the excited state surface along a bond-breaking coordinate; (b) lack of doubly excited configurations; (c) much too low charge transfer excitation energies. These TDDFT failure cases are all strikingly exhibited by prototype two-electron systems such as dissociating H2 and HeH+. We find for these systems with time-dependent density matrix functional theory that: (a) Within previously formulated simple adiabatic approximations, the bonding-to-antibonding excited state surface as well as charge transfer excitations are described without problems, but not the double excitations; (b) An adiabatic approximation is formulated in which also the double excitations are fully accounted for.