Search

The molecular orbital (MO) concept is a useful tool, which relates the molecular ground-state energy with the energies (and occupations) of the individual orbitals. However, analysis of the excited states from linear response computations is performed in terms of the initial state MOs or some other forms of orbitals, e.g., natural or natural transition orbitals. Because these orbitals lack the respective energies, they do not allow developing a consistent orbital picture of the excited states. Herein, we argue that Dyson’s orbitals enable description of the response states compatible with the concepts of molecular orbital theory. The Dyson orbitals and their energies obtained by mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) for the response ground state are remarkably similar to the canonical MOs obtained by the usual DFT calculation. For excited states, the Dyson orbitals provide a chemically sensible picture of the electronic transitions, thus bridging the chasm between orbital theory and response computations.

Internal conversion (IC) between the two lowest singlet excited states, 11Bu+ and 21Ag-, of s-trans-butadiene and s-trans-hexatriene is investigated using a series of single- and multi- reference wave function and density functional theory (DFT) methodologies. Three independent types of the equation-of-motion coupled-cluster (EOMCC) theory capable of providing an accurate and balanced description of one- as well as two-electron transitions, abbreviated as δ-CR-EOMCC(2,3), DIP-EOMCC(4h2p){No}, and DEA-EOMCC(4p2h){Nu} or DEA-EOMCC(3p1h,4p2h){Nu}, consistently predict that the 11Bu+/21Ag- crossing in both molecules occurs along the bond length alternation coordinate. However, the analogous 11Bu+ and 21Ag- potentials obtained with some multireference approaches, such as CASSCF and MRCIS(D), as well as with the linear-response formulation of time-dependent DFT (TDDFT), do not cross. Hence, caution needs to be exercised when studying the low-lying singlet excited states of polyenes with conventional multiconfigurational methods and TDDFT. The multistate many-body perturbation theory methods, such as XMCQDPT2, do correctly reproduce the curve crossing. Among the simplest and least expensive computational methodologies, the DFT approaches that incorporate the contributions of doubly excited configurations, abbreviated as MRSF (mixed reference spin-flip) TDDFT and SSR(4,4), accurately reproduce our best EOMCC results. This is highly promising for nonadiabatic molecular dynamics simulations in larger systems.

Computation of the ionization energies and of the respective Dyson orbitals based on the use of the extended Koopmans theorem (EKT) is implemented in connection with an ensemble density functional theory (eDFT) method, the state-interaction state-averaged spin-restricted ensemble-referenced Kohn–Sham (SI-SA-REKS or SSR) method. The new methodology enables fast computation of the ionization energies and evaluation of the respective Dyson orbitals, the square norms of which are related with the ionization probabilities, in the ground and excited electronic states of molecules. As the application of EKT recycles the intermediate quantities from the SSR analytical energy gradient, evaluation of the ionization energies and probabilities can be carried out on-the-fly during the nonadiabatic molecular dynamics simulations. This opens up a perspective for fast theoretical simulation of the time-resolved photoelectron spectroscopy observations. In the present work, the new methodology is tested in the computation of the ionization energies and Dyson orbitals of several molecules in the ground and excited electronic states, including strongly correlated species, such as the ozone molecule, dissociating chemical bonds, and conical intersections.

The analytic derivative algorithm for the effective contact densities obtained by the exact two-component (X2C) relativistic Hamiltonian is extended to the local approximations to X2C to achieve a higher computational efficiency without losing accuracy. The new algorithm has been implemented in a standalone program, which can utilize the molecular orbitals from state-of-the-art ab initio or density functional theory (DFT) calculations by other quantum chemistry programs in connection with various relativistic Hamiltonians. With the help of the utility program, the effective contact densities as well as the related Mössbauer isomer shifts can be studied by various advanced single-reference and multi-reference ab initio methods as long as the canonical or natural orbitals are available. Using the developed algorithm, the effective contact densities and the Mössbauer isomer shifts in a series of iron compounds and in HgFn (n = 1, 2, 4, and 6) molecules were studied. The obtained results show that (1) adequate account of the static electron correlation significantly improves the agreement of the theoretical 57Fe effective contact densities with the experimental isomer shifts, and (2) the non-monotonous changes of the effective contact density in a series of HgFn compounds are caused by the increasing screening effect due to shrinking of the Hg 5d orbitals.

Time-resolved photoelectron spectra during the photochemical ring-opening reaction of 1,3-cyclohexadiene (CHD) are modeled by an ensemble density functional theory (eDFT) method. The computational methodology employed in this work is capable of correctly describing the multi-reference effects arising in the ground and excited electronic states of molecules, which is important for the correct description of the ring-opening reaction of CHD. The geometries of molecular species along the non-adiabatic molecular dynamics (NAMD) trajectories reported in a previous study of the CHD photochemical ring-opening were used in this work to calculate the ionization energies and the respective Dyson orbitals for all possible ionization channels. The obtained theoretical time-resolved spectra display decay characteristics in a reasonable agreement with the experimental observations; i.e., the decay (and rise) of the most mechanistically significant signals occurs on the timescale of 100–150 fs. This is very different from the excited state population decay characteristics (τ _(S1) = 234 ± 8 fs) obtained in the previous NAMD study. The difference between the population decay and the decay of the photoelectron signal intensity is traced back to the geometric transformation that the molecule undergoes during the photoreaction. This demonstrates the importance of including the geometric information in interpretation of the experimental observations.

Internal conversion (IC) between the two lowest singlet excited states, 1

Sudden polarization (SP) is one of the manifestations of electron transfer in the electronically excited states of molecules. Proposed initially to explain the unusual reactivity of photoexcited olefins, SP often occurs in the excited states of molecules possessing strongly correlated diradical ground state. Theoretical description of SP involves mixing between the singly excited and the doubly excited zwitterionic states, which makes it inaccessible with the use of the popular linear-response time-dependent density functional theory methods. In this work, an extended variant of the state-interaction state-averaged spin-restricted ensemble-referenced Kohn–Sham (SI-SA-REKS, or SSR) method is applied to study SP in a number of organic diradical systems. To this end, the analytical derivative formalism is derived and implemented for the SSR(3,2) method (see the main text for explanation of the acronym), which enables the automatic geometry optimization and obtains the relaxed density matrices as well as the electron binding energies and respective Dyson’s orbitals. Application of the new method to SP in the lowest singlet excited state of ethylene agrees with the results obtained previously with the use of multireference methods of wavefunction theory. A number of interesting manifestations of SP are observed, such as the charge transfer in photoexcited tetramethyleneethene (TME) diradical mediated by the vibrational motion and conductivity switching in the excited state of a donor–acceptor dyad placed in an external electric field.

The mixed-reference spin-flip (MRSF) time-dependent density functional theory (TDDFT) method eliminates the notorious spin contamination of SF-TDDFT, thus enabling identification of states of proper spin-symmetry for automatic geometry optimization and molecular dynamics simulations. Here, we analyze and optimize the MRSF-TDDFT in the calculations of the vertical excitation energies (VEEs) and the singlet-triplet (ST) gaps. The dependence of the obtained VEEs and ST gaps on the intrinsic parameters of the MRSF-TDDFT method is investigated, and prescriptions for the proper use of the method are formulated. For VEEs, MRSF-TDDFT displays similar or better accuracy than SF-TDDFT (ca. 0.5 eV), while considerably outperforming the LR-TDDFT for the ST gaps. As a result, a new functional of STG1X (dubbed here), especially for ST gaps is suggested on the basis of splitting between the components of the atomic multiplets.

The mixed-reference spin-flip time-dependent density functional theory (MRSF-TD-DFT) method eliminates the erroneous spin contamination of the SF-TD-DFT methodology, while retaining the conceptual and practical simplicity of the latter. The availability of the analytic gradient of the energy of the MRSF-TD-DFT response states enables automatic geometry optimization of the targeted states. Here, we apply the new method to optimize the geometry of several S1/S0 conical intersections occurring in typical organic molecules. We demonstrate that MRSF-TD-DFT is capable of producing the correct double-cone topology of the intersections and describing the geometry of the lowest-energy conical intersections and their relative energies with accuracy matching that of the best multireference wavefunction ab initio methods. In this regard, MRSF-TD-DFT differs from many popular single-reference methods, such as, e.g., the linear response TD-DFT method, which fail to produce the correct topology of the intersections. As the new methodology completely eliminates the ambiguity with the identification of the response states as proper singlets or triplets, which is plaguing the SF-TD-DFT calculations, it can be used for automatic geometry optimization and molecular dynamic simulations not requiring constant human intervention.

Photodecomposition of cyclopropanone is investigated by static quantum chemical calculations and non-adiabatic molecular dynamics (NAMD) simulations. The quantum chemical calculations are carried out by an ensemble density functional theory (eDFT) method capable of delivering high quality results for the ground and excited electronic states of molecules with dissociating bonds. In the NAMD simulations, this method is combined with a novel trajectory surface hopping (TSH) methodology derived from the exact factorization of the electronic-nuclear wavefunction. An ultrafast biexponential decay of the S1 state of cyclopropanone is predicted, where the short (ca. 30 fs) decay time is due to the trajectories reaching the conical intersection (CI) seam on the first approach and the long (ca. 120 fs) decay time is due to recrossing of the CI seam. The experimentally observed dependence of the dissociation (C3H4O* → C2H4 + CO) quantum yield on the excitation wavelength is correctly reproduced by the NAMD simulations. The dependence is explained by the necessity to excite certain vibrational normal modes (e.g., a ring stretching mode with the frequency of 769 cm-1) to break a lateral C-C bond remaining intact at the geometries of the CI seam.

We present a fast and accurate numerical algorithm for computing the first-order nonadiabatic coupling matrix element (NACME). The algorithm employs the truncated Leibniz formula (TLF) approximation within the finite-difference method, which makes it easily applicable in connection with any wave function-based methodology. In this work, we used the algorithm in connection with the recently developed mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT, MRSF for brevity). The accuracy is assessed for NACME between the singlet electronic states of a dissociating hydrogen molecule. It is demonstrated that an intermediate approximation, TLF(1), affords a negligible numeric error on the order of ∼10⁻¹⁰ a.u. while enabling a fast computation of NACME. As the MRSF method yields the correct description of the dissociation curves of H₂ for all the electronic states involved, the numeric TLF(1)/MRSF NACME values are in excellent agreement with the reference analytical values obtained by the full configuration interaction. For polyatomic molecules, the MRSF NAC vectors agree very closely with the MRCISD NAC vectors. Hence, the proposed protocol is a promising tool for the evaluation of NACMEs.

Non-radiative relaxation of the photoexcited thymine in the gas phase shows an unusually long excited-state lifetime, and, over the years, a number of models, i.e., S1-trapping, S2-trapping, and S1&S2-trapping, have been put forward to explain its mechanism. Here, we investigate this mechanism using non-adiabatic molecular dynamics (NAMD) simulations in connection with the recently developed mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) method. We show that the previously predicted S2-trapping model was due to an artifact caused by an insufficient account of the dynamic electron correlation. The current work supports the S1-trapping mechanism with two lifetimes, τ1 = 30 ± 1 fs and τ2 = 6.1 ± 0.035 ps, quantitatively consistent with the recent time-resolved experiments. Upon excitation to the S2 (ππ*) state, thymine undergoes an ultrafast (ca. 30 fs) S2→S1 internal conversion and resides around the minimum on the S1 (nOπ*) surface, slowly decaying to the ground state (ca. 6.1 ps). While the S2→S1 internal conversion is mediated by fast bond length alternation distortion, the subsequent S1→S0 occurs through several conical intersections, involving a slow puckering motion.

The non-adiabatic dynamics of furan excited in the ππ* state (S2 in the Franck–Condon geometry) was studied using non-adiabatic molecular dynamics simulations in connection with an ensemble density functional method. The time-resolved photoelectron spectra were theoretically simulated in a wide range of electron binding energies that covered the valence as well as the core electrons. The dynamics of the decay (rise) of the photoelectron signal were compared with the excited-state population dynamics. It was observed that the photoelectron signal decay parameters at certain electron binding energies displayed a good correlation with the events occurring during the excited-state dynamics. Thus, the time profile of the photoelectron intensity of the K-shell electrons of oxygen (decay constant of 34 ± 3 fs) showed a reasonable correlation with the time of passage through conical intersections with the ground state (47 ± 2 fs). The ground-state recovery constant of the photoelectron signal (121 ± 30 fs) was in good agreement with the theoretically obtained excited-state lifetime (93 ± 9 fs), as well as with the experimentally estimated recovery time constant (ca. 110 fs). Hence, it is proposed to complement the traditional TRPES observations with the trXPS (or trNEXAFS) measurements to obtain more reliable estimates of the most mechanistically important events during the excited-state dynamics.

Conical intersections control excited state reactivity, and thus, elucidating and predicting their geometric and energetic characteristics are crucial for understanding photochemistry. Locating these intersections requires accurate and efficient electronic structure methods. Unfortunately, the most accurate methods (e.g., multireference perturbation theories such as XMS-CASPT2) are computationally challenging for large molecules. The state-interaction state-averaged restricted ensemble referenced Kohn-Sham (SI-SA-REKS) method is a computationally efficient alternative. The application of SI-SA-REKS to photochemistry was previously hampered by a lack of analytical nuclear gradients and nonadiabatic coupling matrix elements. We have recently derived analytical energy derivatives for the SI-SA-REKS method and implemented the method effectively on graphical processing units. We demonstrate that our implementation gives the correct conical intersection topography and energetics for several examples. Furthermore, our implementation of SI-SA-REKS is computationally efficient, with observed sub-quadratic scaling as a function of molecular size. This demonstrates the promise of SI-SA-REKS for excited state dynamics of large molecular systems.

In different branches of time-dependent density functional theory (TDDFT), the static and dynamic electron correlation enters in different ways. The standard spin-conserving linear response (LR-TDDFT) methodology includes explicitly the contributions of the singly-excited configurations; however, it relies on an implicit account of the electron correlation through an (approximate) exchange-correlation (XC) functional. In the mixed-reference spin-flip TDDFT (MRSF-TDDFT), a number of doubly-excited (DE) configurations are explicitly included in the description of their response states. Here, the importance of the explicit account of DE is investigated for the lowest four excited singlet states of all-trans-polyenes up to C₂₄H₂₆. For the optically bright 1B_u⁺ state, the DE contribution in MRSF-TDDFT approaches 10% with the increasing system size. For the optically dark 2Ag⁻ state, the DE contribution increases from ca. 13% (C₄H₆) to nearly 30% (C₂₄H₂₆). An even more considerable DE contribution (∼50%) is observed in the higher 1B_u⁻ states. As LR-TDDFT misses these contributions entirely, its ability to accurately describe the excited states is limited by the XC functional. The hybrid XC functionals with a small fraction of the exact exchange, e.g., B3LYP, may mimic certain effects of DE through the self-interaction error (SIE). However, the description of the 1B_u⁺ state by LR-TDDFT remains poor. On the other hand, MRSF-TDDFT can flexibly take an implicit (through the XC functional) and an explicit (through DE) account of the electron correlation, which enables a more balanced description of various types of the excited states regardless of their character, thus reducing the chances of failure.

Due to their multiconfigurational nature featuring strong electron correlation, accurate description of diradicals and diradicaloids is a challenge for quantum chemical methods. The recently developed mixed-reference spin-flip (MRSF)-TDDFT method is capable of describing the multiconfigurational electronic states of these systems while avoiding the spin-contamination pitfalls of SF-TDDFT. Here, we apply MRSF-TDDFT to study the adiabatic singlet-triplet (ST) gaps in a series of well-known diradicals and diradicaloids. On average, MRSF displays a very high prediction accuracy of the adiabatic ST gaps with the mean absolute error (MAE) amounting to 0.14 eV. In addition, MRSF is capable of accurately describing the effect of the Jahn-Teller distortion occurring in the trimethylenemethane diradical, the violation of the Hund rule in a series of the didehydrotoluene diradicals, and the potential energy surfaces of the didehydrobenzene (benzyne) diradicals. A convenient criterion for distinguishing diradicals and diradicaloids is suggested on the basis of the easily obtainable quantities. In all of these cases, which are difficult for the conventional methods of density functional theory (DFT), MRSF shows results consistent with the experiment and the high-level ab initio computations. Hence, the present study documents the reliability and accuracy of MRSF and lays out the guidelines for its application to strongly correlated molecular systems.

The computation of electron attachment energies (electron affinities) was implemented in connection with an ensemble density functional theory method, the state-interaction state-averaged spin-restricted ensemble-referenced Kohn-Sham (SI-SA-REKS or SSR) method. With the use of the extended Koopmans' theorem, the electron affinities and the respective Dyson orbitals are obtained directly for the neutral molecule, thus avoiding the necessity to compute the ionized system. Together with the EKT-SSR (extended Koopmans' theorem-SSR) method for ionization potentials, which was developed earlier, EKT-SSR for electron affinities completes the implementation of the EKT-SSR formalism, which can now be used for obtaining electron detachment as well as the electron attachment energies of molecules in the ground and excited electronic states. The extended EKT-SSR method was tested in the calculation of several closed-shell molecules. For the molecules in the ground states, the EKT-SSR energies of Dyson's orbitals are virtually identical to the energies of the unoccupied orbitals in the usual single-reference spin-restricted Kohn-Sham calculations. For the molecules in the excited states, EKT-SSR predicts an increase of the most positive electron affinity by approximately the amount of the vertical excitation energy. The electron affinities of a number of diradicals were calculated with EKT-SSR and compared with the available experimental data. With the use of a standard density functional (BH&HLYP), the EKT-SSR electron affinities deviate on average by ca. 0.2 eV from the experimental data. It is expected that the agreement with the experiment can be improved by designing density functionals parametrized for ionization energies.

In this work, a direct mixed quantum-classical dynamics approach is presented, which combines two new computational methodologies. The nuclear dynamics is solved by the decoherence-induced surface hopping based on the exact factorization (DISH-XF) method, which is derived from the exact factorization of the electronic-nuclear wave function and correctly describes quantum decoherence phenomena. The state-interaction state-averaged spin-restricted ensemble-referenced Kohn-Sham (SI-SA-REKS, or SSR, for brevity) electronic structure method is based on ensemble density functional theory (eDFT) and provides correct description of real crossings between the ground and excited Born–Oppenheimer electronic states. The new combined approach has been applied to the excited-state nonadiabatic dynamics of the trans-penta-2,4-dieniminium cation (PSB3). The predicted S1 lifetime of trans-PSB3, τ = 99 ± 51 fs, and the quantum yield of the cis conformation, ϕ = 0.63, agree with the results obtained previously in nonadiabat...

Since carbon tetrahalides CX4 (X = Cl/Br) are adsorbed on carbon nanotubes and graphene sheets, we have studied the structures, adsorption energies, and electronic properties of CX4 adsorbed on benzene, coronene, and graphene using dispersion corrected density functional theory (DFT) with hybrid functionals. As compared with the benzene–CX4 complexes (with binding energy of ∼14/15 kJ/mol) where electrostatic energy is significant due to the halogen bonding effect, the graphene–CX4 complexes show about three times the benzene–CX4 binding energy (∼40/45 kJ/mol) where the dispersion interaction is overwhelming with insignificant electrostatic energy. Since the X atoms in CX4 are slightly positively charged and the X atom’s ends are particularly more positively charged due to the σ-hole effect, CX4 behaves as an electron acceptor. This results in electron transfer from locally negatively charged C sites of benzene/coronene to CX4. In contrast, no electron transfer occurs from graphene to CX4 because of the la...

Various types of interactions between halogen (X) and π moiety (X−π interaction) including halogen bonding play important roles in forming the structures of biological, supramolecular, and nanomaterial systems containing halogens and aromatic rings. Furthermore, halogen molecules such as X2 and CX4 (X = Cl/Br) can be intercalated in graphite and bilayer graphene for doping and graphene functionalization/modification. Due to the X−π interactions, though recently highly studied, their structures are still hardly predictable. Here, using the coupled-cluster with single, double, and noniterative triple excitations (CCSD(T)), the Moller–Plesset second-order perturbation theory (MP2), and various flavors of density functional theory (DFT) methods, we study complexes of benzene (Bz) with halogen-containing molecules X2 and CX4 (X = Cl/Br) and analyze various components of the interaction energy using symmetry adapted perturbation theory (SAPT). As for the lowest energy conformers (S1), X2–Bz is found to have the...

Ensemble density functional theory (DFT) furnishes a rigorous theoretical framework for describing the non-dynamic electron correlation arising from (near) degeneracy of several electronic configurations. Ensemble DFT naturally leads to fractional occupation numbers (FONs) for several Kohn-Sham (KS) orbitals, which thereby become variational parameters of the methodology. The currently available implementation of ensemble DFT in the form of the spin-restricted ensemble-referenced KS (REKS) method was originally designed for systems with only two fractionally occupied KS orbitals, which was sufficient to accurately describe dissociation of a single chemical bond or the singlet ground state of biradicaloid species. To extend applicability of the method to systems with several dissociating bonds or to polyradical species, more fractionally occupied orbitals must be included in the ensemble description. Here we investigate a possibility of developing the extended REKS methodology with the help of the generalized valence bond (GVB) wavefunction theory. The use of GVB enables one to derive a simple and physically transparent energy expression depending explicitly on the FONs of several KS orbitals. In this way, a version of the REKS method with four electrons in four fractionally occupied orbitals is derived and its accuracy in the calculation of various types of strongly correlated molecules is investigated. We propose a possible scheme to ameliorate the partial size-inconsistency that results from perfect spin-pairing. We conjecture that perfect pairing natural orbital (NO) functionals of reduced density matrix functional theory (RDMFT) should also display partial size-inconsistency.

The dynamics of the ring opening in the S1 state of cyclohexa-1,3-diene (CHD) is studied by a new direct mixed quantum-classical non-adiabatic dynamics approach which employs the decoherence-induced surface hopping based on the exact factorisation (DISH-XF) molecular dynamics method in connection with the state-interaction state-averaged spin-restricted ensemble-referenced Kohn–Sham (SI-SA-REKS, or SSR) electronic structure method. The critical species on the S0 and S1 PESs of CHD were studied using the SSR method and the minimum energy pathways (MEPs) were optimised. The obtained vertical excitation energies are in good agreement (within ca. 5–6 kcal/mol) with the experimental values. The optimised geometry of the S1/S0 minimum energy conical intersection (MECI) agrees well with the previously obtained MSPT2 geometry. The DISH-XF/SSR non-adiabatic molecular dynamics (NAMD) simulations of ring opening in CHD predict the S1 exponential decay constant τ=234±8 fs in a reasonable agreement with an experimental estimate (230±30 fs). The calculated product branching ratio (CHD:HT = 64:36) is in agreement with the recent experimental measurement (70:30). The NAMD trajectories are analysed in terms of the vibrational normal modes and the obtained branching ratio is explained by persistent stretching of the fissile bond when the trajectories propagate on the S1 PES.

An extended variant of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method, the REKS(4,4) method, designed to describe the ground electronic states of strongly multireference systems is modified to enable calculation of excited states within the time-independent variational formalism. The new method, the state-interaction state-averaged REKS(4,4), i.e., SI-SA-REKS(4,4), is capable of describing several excited states of a molecule involving double bond cleavage, polyradical character, or multiple chromophoric units. We demonstrate that the new method correctly describes the ground and the lowest singlet excited states of a molecule (ethylene) undergoing double bond cleavage. The applicability of the new method for excitonic states is illustrated with π stacked ethylene and tetracene dimers. We conclude that the new method can describe a wide range of multireference phenomena.

Ensemble density functional theory (DFT) is a novel theoretical approach that is capable of exact treatment of non-dynamic electron correlation in the ground and excited states of many-body fermionic systems. In contrast to ordinary DFT, ensemble DFT has not found so far a way to the repertoire of methods of modern computational chemistry, probably owing to the lack of practically affordable implementations of the theory. The spin-restricted ensemble-referenced Kohn–Sham (REKS) method represents perhaps the first computational scheme that makes ensemble DFT calculations feasible. The REKS method is based on the rigorous ensemble representation of the energy and the density of a strongly correlated system and provides for an accurate and consistent description of molecular systems the electronic structure of which is dominated by the non-dynamic correlation. This includes the ground and excited states of molecules undergoing bond breaking/bond formation, the low-spin states of biradicals and polyradicals, symmetry forbidden chemical reactions and avoided crossings of potential energy surfaces, real intersections between the energy surfaces of the ground and excited states (conical intersections), and many more. The REKS method can be employed in connection with any local, semi-local and hybrid (global and range-separated) functional and affords calculations of large and very large molecular systems at a moderate mean-field cost. WIREs Comput Mol Sci 2015, 5:146–167. doi: 10.1002/wcms.1209 Conflict of interest: The author has declared no conflicts of interest for this article. For further resources related to this article, please visit the WIREs website.

WOS:000340351200020; International audience; We report and characterize ground-state and excited-state potential energy profiles using a variety of electronic structure methods along a loop lying on the branching plane associated with a conical intersection (Cl) of a reduced retinal model, the penta-2,4-dieniminium cation (PSB3). Whereas the performance of the equation-of-motion coupled-duster, density functional theory, and multireference methods had been tested along the excited- and ground-state paths of PSB3 in our earlier work, the ability of these methods to correctly describe the potential energy surface shape along a CI branching plane has not yet been investigated. This is the focus of the present contribution. We find, in agreement with earlier studies by others, that standard time-dependent DFT (TDDFT) does not yield the correct two-dimensional (i.e., conical) crossing along the branching plane but rather a one-dimensional (i.e., linear) crossing along the same plane. The same type of behavior is found for SS-CASPT2(IPEA=0), SS-CASPT2(IPEA=0.25), spin-projected SF-TDDFT, EOM-SF-CCSD, and, finally, for the reference MRCISD+Q method. In contrast, we found that MRCISD, CASSCF, MS-CASPT2(IPEA=0), MS-CASPT2(IPEA=0.25), XMCQDPT2, QD-NEVPT2, non-spin-projected SF-TDDFT, and SI-SA-REKS yield the expected conical crossing. To assess the effect of the different crossing topologies (i.e., linear or conical) on the PSB3 photoisomerization efficiency, we discuss the results of 100 semiclassical trajectories computed by CASSCF and SS-CASPT2(IPEA=0.25) for a PSB3 derivative. We show that for the same initial conditions, the two methods yield similar dynamics leading to isomerization quantum yields that differ by only a few percent.

WOS:000361155700013; International audience; Conical intersections are perhaps the most significant mechanistic features of chemical reactions occurring through excited states. By providing funnels for efficient non-adiabatic population transfer, conical intersections govern the branching ratio of products of such reactions, similar to what the transition states do for ground-state reactivity. In this regard, intersections between the ground and the lowest excited states play a special role, and the correct description of the potential energy surfaces in their vicinity is crucial for understanding the mechanism and dynamics of excited-state reactions. The methods of density functional theory, such as time-dependent density functional theory, are widely used to describe the excited states of large molecules. However, are these methods suitable for describing the conical intersections or do they lead to artifacts and, consequently, to erroneous description of reaction dynamics? Here we address the first part of this question and analyze the ability of several density functional approaches, including the linear-response time-dependent approach as well as the spin-flip and ensemble formalisms, to provide the correct description of conical intersections and the potential energy surfaces in their vicinity. It is demonstrated that the commonly used linear-response time-dependent theory does not yield a proper description of these features and that one should instead use alternative computational approaches.

He-droplet-induced changes of the hyperfine structure constants of alkali-metal atoms are investigated by a combination of relativistically corrected ab initio methods with a simulation of the helium density distribution based on He density functional theory. Starting from an accurate description of the variation of the hyperfine structure constant in the M-He diatomic systems (M=Li, Na, K, Rb) as a function of the interatomic distance we simulate the shifts induced by droplets of up to 10,000 (4)He atoms. All theoretical predictions for the relative shifts in the isotropic hyperfine coupling constants of the alkali-metal atoms attached to helium droplets of different size are then tied to a single, experimentally derived parameter of Rb.

Analytical second derivatives for the normalized elimination of the small component (NESC) method are derived for the first time and implemented for the routine calculation of NESC vibrational frequencies and other second order molecular properties using the scalar relativistic form of NESC. Using response theory, the second derivatives of the transformation matrix U connecting the large and the pseudolarge components of the relativistic wave function are correctly derived. The 24 derivative terms involving the NESC Hamiltonian and the NESC renormalization matrix are individually tested, and their contributions to the energy Hessian are calculated. The influence of a finite nucleus model and that of the picture change is determined. Different ways of speeding up the calculation of the NESC second derivatives are tested. It is shown that second order properties can routinely be calculated in combination with Hartree-Fock, density functional theory, Moller-Plesset perturbation theory, and any other electron correlation corrected quantum chemical method provided analytic second derivatives are available in the nonrelativistic case. The general applicability of the analytic NESC Hessian is demonstrated by benchmark calculations for NESC/DFT calculations involving up to 1500 basis functions.

The ground state electronic structure of the high-temperature (HT) and the low-temperature (LT) phases of (EDO-TTF)(2)PF6 is investigated using the embedded cluster approach in combination with the density functional method designed to describe the strong non-dynamic electron correlation. It is found that, in the HT phase, the unpaired electron spins located on pairs of neighbouring EDO-TTF molecules are antiferromagnetically coupled along the stacking direction with the Heisenberg exchange integral J = -655 cm(-1). In the LT phase, the unpaired spins located on the cationic EDO-TTF molecules are coupled antiferromagnetically with J values strongly alternating along the stacking axis of the crystal thus rendering it diamagnetic. The parameters of the extended Hubbard model are evaluated and the conductance properties of the two phases are estimated using these parameters. It is suggested to investigate the charge and spin excitations in the two phases of (EDO-TTF)(2)PF6 with the use of angle-resolved photoemission spectroscopy.

The use of the mixed reference (MR) reduced density matrix, which combines reduced density matrices of the MS = +1 and −1 triplet-ground states, is proposed in the context of the collinear spin-flip-time-dependent density functional theory (SF-TDDFT) methodology. The time-dependent Kohn-Sham equation with the mixed state is solved by the use of spinor-like open-shell orbitals within the linear response formalism, which enables to generate additional configurations in the realm of TD-DFT. The resulting MR-SF-TDDFT computational scheme has several advantages before the conventional collinear SF-TDDFT. The spin-contamination of the response states of SF-TDDFT is nearly removed. This considerably simplifies the identification of the excited states, especially in the “black-box” type applications, such as the automatic geometry optimization, reaction path following, or molecular dynamics simulations. With the new methodology, the accuracy of the description of the excited states is improved as compared to the collinear SF-TDDFT. Several test examples, which include systems typified by strong non-dynamic correlation, orbital (near) degeneracy, and conical intersections, are given to illustrate the performance of the new method.

Mossbauer spectroscopy is a widely used analytic too] which provides information about local electronic Structure of solid materials on an atomic scale. The isomer shift of resonance nuclear gamma transition is a sensitive parameter which depends on the charge and spin state of the resonating atom as well as on its chemical environment. Theory underlying the isomer shift is reviewed and its connection to the local electronic structure is discussed. A review of advances made in the ob initio calculation of isomer shift is presented. The importance of Careful calibration of the parameters of nuclear gamma transitions on the basis of high-level quantum chemical calculations With the inclusion of both relativistic effects and electron correlation is underlined. With the help of accurate theoretical calculations of the isomer shift over a wide range of chemical environments deeper understanding of a relationship between the observed spectroscopic parameters and the electronic Structure of materials will be gained. (C) 2008 Elsevier B.V. All rights reserved.

A time-independent density functional approach to the calculation of excitation energies from the ground states of molecules typified by the strong nondynamic electron correlation is suggested. The new method is based on the use of the spin-restricted ensemble-referenced Kohn-Sham formalism [Filatov, M.; Shaik, S. Chem. Phys. Lett. 1999, 304, 429] for the calculation of the ground state. In the new method, the average energy of the ground state and a state created by a single excitation thereof is minimized with respect to the Kohn-Sham orbitals and their fractional occupation numbers. The lowest singlet excitation energies obtained with the help of the new formalism for a number of model systems, such as the hydrogen molecule with stretched bond, twisted ethylene, and twisted hexa-1,3,5-triene, are compared with the results of the time-dependent density functional theory, with the results of ab initio CASSCF/CASPT2 calculations, and with the experimental data. Applicability of the new method to the description of photochemical reactions is discussed.

With the help of a recently suggested computational scheme [J. Chem. Phys. 2007, 127, 084101], Mossbauer isomer shifts are calculated within the context of density functional theory, for a series of iron containing compounds. The influence of the choice of a density functional and of the truncation of a basis set on the results of calculations is analyzed. It has been observed that the hybrid density functionals, especially BH&HLYP, provide better correlation with experimental results than pure density functionals. The analysis of basis set truncation reveals that the addition (or removal) of the tightmost primitive functions to a large uncontracted basis set has only a minor influence on the calculated isomer shift values. It is observed that, with the use of a small contracted basis set, a reasonable accuracy for the calculated isomer shifts can be achieved.

The performance of density functional theory in estimating the magnetic coupling constant in a series of Cu(II) binuclear complexes is investigated by making use of two open shell formalisms: the broken symmetry and the spin-restricted ensemble-referenced Kohn-Sham methods. The strong dependence of the calculated magnetic coupling constants with respect to the exchange-correlation functional is confirmed and found to be independent of whether spin symmetry is imposed or not. The use of a method which guarantees the spin state does not improve the correlation with the experiment and indeed shows some worsening due to an overestimation of the ferromagnetic interactions. However, with the present exchange-correlation functionals, a rather systematic deviation is found. Therefore, it would be possible to develop improved density functionals which will allow for a rigorous treatment of open shell systems in density functional theory.

Ensemble density functional theory (DFT) is a novel time-independent formalism for obtaining excitation energies of many-body fermionic systems. A considerable advantage of ensemble DFT over the more common Kohn–Sham (KS) DFT and time-dependent DFT formalisms is that it enables one to account for strong non-dynamic electron correlation in the ground and excited states of molecular systems in a transparent and accurate fashion. Despite its positive aspects, ensemble DFT has not so far found its way into the repertoire of methods of modern computational chemistry, probably because of the perceived lack of practically affordable implementations of the theory. The spin-restricted ensemble-referenced KS (REKS) method is perhaps the first computationally feasible implementation of the ideas behind ensemble DFT which enables one to describe accurately electronic transitions in a wide class of molecular systems, including strongly correlated molecules (biradicals, molecules undergoing bond breaking/formation), extended π-conjugated systems, donor–acceptor charge transfer adducts, etc.

A new method for calculating the indirect nuclear spin-spin coupling constant within the regular approximation to the exact relativistic Hamiltonian is presented. The method is completely analytic in the sense that it does not employ numeric integration for the evaluation of relativistic corrections to the molecular Hamiltonian. It can be applied at the level of conventional wave function theory or density functional theory. In the latter case, both pure and hybrid density functionals can be used for the calculation of the quasirelativistic spin-spin coupling constants. The new method is used in connection with the infinite-order regular approximation with modified metric (IORAmm) to calculate the spin-spin coupling constants for molecules containing heavy elements. The importance of including exact exchange into the density functional calculations is demonstrated.

95 Mo and 183 W, respectively, which is mainly due to improvements in the paramagnetic contributions. The relationship between electronic structure of a molecule and the relativistic paramagnetic corrections is discussed. Relativistic effects for the diamagnetic part of the magnetic shielding caused by a relativistic contraction of the s, p orbitals in the core region concern only the shielding values, however, have little consequence for the shift values because of the large independence from electronic structure and a cancellation of these effects in the shift values. It is shown that the relativistic corrections can be improved by level shift operators and a B3LYP hybrid functional, for which Hartree‐Fock exchange is reduced to 15%. © 2003 American Institute of Physics. @DOI: 10.1063/1.1580091#

The 2A″ state of ClOO and BrOO is investigated using single- and multi-reference DFT (density functional theory). Unrestricted DFT (UDFT) carried out with the BLYP functional exaggerates the ionic character of the X–OO bond (X=Cl, Br) and by this its bond dissociation energy, while UDFT with the B3LYP functional underestimates ionic character and bond dissociation energies. In previous investigations, this was overlooked and has led to a misleading interpretation of single determinant UDFT results for XOO(2A″). The correct description of the two radicals can only be achieved by an appropriate configurational state mixing between the doublet ground state (leading to ionic character) and the first excited state, a quartet state leading to covalent character of the X–OO bond. REKS (spin-restricted ensemble-referenced Kohn–Sham) theory provides a multireference description and is capable of describing bonding in XOO correctly provided it is carried out with the hybrid functional B3LYP rather than BLYP, which is the result of the self-interaction error of pure exchange functionals and the inclusion of non-specific non-dynamic electron correlation effects suppressing the specific ones introduced by REKS. The relevance of results for the DFT description of related compounds (HOOO radical, MOO with M =Cu, Ag or Au) or multi-reference problems in general is discussed.

Unrestricted density functional theory (UDFT) can be used for the description of open-shell singlet (OSS) biradicals provided a number of precautions are considered. Biradicals that require a two-determinantal wave function (e.g. OSS state of carbenes) cannot be described by UDFT for principal reasons. However, if the overlap between the open-shell orbitals is small (the single electrons are located at different atomic centers) errors become small and, then, the principal failure of UDFT in these cases is not apparent and may even be disguised by the fact that UDFT has the advantage of describing spin polarization better than any restricted open shell DFT method. In the case of OSS biradicals with two- or multiconfigurational character (but a onedeterminantal form of the leading configuration), reasonable results can be obtained by broken-symmetry (BS)-UDFT, however in each case this has to be checked. In no case is it reasonable to lower the symmetry of a molecule to get a suitable UDFT description. Hybrid functionals such as B3LYP perform better than pure DFT functionals in BS-UDFT calculations because the former reduce the self-interaction error of DFT exchange functionals, which mimics unspecified static electron correlation effects, so that the inclusion of specific static electron correlation effects via the form of the wavefunction becomes more effective.

A recently developed variationally stable quasi-relativistic method, which is based on the low-order approximation to the method of normalized elimination of the small component, was incorporated into density functional theory (DFT). The new method was tested for diatomic molecules involving Ag, Cd, Au, and Hg by calculating equilibrium bond lengths, vibrational frequencies, and dissociation energies. The method is easy to implement into standard quantum chemical programs and leads to accurate results for the benchmark systems studied.

The state-averaged (SA) spin restricted ensemble referenced Kohn-Sham (REKS) method and its state interaction (SI) extension, SI-SA-REKS, 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 non-adiabatic reactions. Until recently, application of the SA-REKS and SI-SA-REKS methods to modeling the dynamics of such reactions was obstructed due to the lack of the analytical energy derivatives. In this work, the analytical derivatives of the individual SA-REKS and SI-SA-REKS energies are derived. The final analytic gradient expressions are formulated entirely in terms of traces of matrix products and are presented in the form convenient for implementation in the traditional quantum chemical codes employing basis set expansions of the molecular orbitals. The implementation and benchmarking of the derived formalism will be described in a subsequent article of this series.

Electron transfer in the ground and excited states of a model donor–acceptor (D–A) system is investigated using the single-reference and multi-reference density functional theory (DFT) methods. To analyze the results of the calculations, a simple two-site multi-reference model was derived that predicts a stepwise electron transfer in the S 0 state and a wave-like dependence of the S 1 electron transfer on the external stimulus. The standard single-reference Kohn-Sham (KS) DFT approach and the time-dependent DFT (TDDFT) method failed to describe the correct dependence of the S 0 and S 1 electron transfer on the external electric field applied along the donor–acceptor system. The multi-reference DFT approach, the spin-restricted ensemble-referenced KS (REKS) method, was able to successfully reproduce the correct behavior of the S 0 and S 1 electron transfer on the applied field. The REKS method was benchmarked against experimentally measured gas phase charge transfer excitations in a series of organic donor–acceptor complexes and displayed its ability to describe this type of electronic transitions with a very high accuracy, mean absolute error of 0.05 eV with the use of the standard range separated density functionals. On the basis of the calculations undertaken in this work, it is suggested that the non-adiabatic coupling between the S 0 and S 1 states may interfere with the electron transfer in a weakly coupled donor–acceptor system. It is also suggested that the electronic excitation of a D+–A− system may play a dual role by assisting the further electron transfer at certain magnitudes of the applied electric field and causing the backward transfer at lower electric field strengths.

WOS:000342844100027; International audience; Quantum-chemical computational methods are benchmarked for their ability to describe conical intersections in a series of organic molecules and models of biological chromophores. Reference results for the geometries, relative energies, and branching planes of conical intersections are obtained using ab initio multireference configuration interaction with single and double excitations (MRCISD). They are compared with the results from more approximate methods, namely, the state-interaction state-averaged restricted ensemble-referenced Kohn-Sham method, spin-flip time-dependent density functional theory, and a semiempirical MRCISD approach using an orthogonalization-corrected model. It is demonstrated that these approximate methods reproduce the ab initio reference data very well, with root-mean-square deviations in the optimized geometries of the order of 0.1 angstrom or less and with reasonable agreement in the computed relative energies. A detailed analysis of the branching plane vectors shows that all currently applied methods yield similar nuclear displacements for escaping the strong non-adiabatic coupling region near the conical intersections. Our comparisons support the use of the tested quantum-chemical methods for modeling the photochemistry of large organic and biological systems. (C) 2014 AIP Publishing LLC.

Density functional theory (DFT) predicts that bicyclo[3.1.0]hexatriene ( 2 ) is more stable than its isomer m -benzyne ( 1 ). Hess [Eur. J. Org. Chem. (2001) 2185] has argued that experimental findings suggesting 1 can equally or even better be associated with 2 . However, high level ab initio calculations (CCSD(T), CASPT2) show that 2 does not exist and that the previously measured infrared spectrum is correctly assigned to 1 . Bond stretch isomers are possible for p -benzynes but not for m -benzynes. The electrophilic character of m -benzynes is in line with 1 but not with 2 .

REKS-type (Filatov, M.; Shaik, S. Chem. Phys. Lett. 1999, 304, 429) density functional calculations were carried out for the lowest energy singlet and triplet states of tetramethyleneethane (TME) diradical. The calculations indicate that the ground state of TME in the gas phase is the singlet state, whereas the triplet state should be metastable at low temperatures. The triplet state metastability derives from the energetic preference for the triplet state at the optimal triplet molecular geometry and from the extremely small (

A recently proposed spin-restricted open-shell Kohn–Sham (ROKS) method is applied to investigate various atomic and molecular multiplet states. A wide range of multiplets is considered: multiplet terms for which the spin-restricted open-shell theory of Roothaan applies, as well as state situations which cannot be described by Roothaan’s theory (e.g., states of square cyclobutadiene, etc.). Problems associated with the use of approximate density functionals and possible perspectives of the ROKS method are discussed.

Density functional theory (DFT) is used to study model ferryl species of cytochrome P-450 and horseradish peroxidase (HRP), as well as of the product complex due to oxidation of H2 by the P-450 species (1–4 and 7). The ferryl species studied include neutral and cation radical states of the porphyrin, as well as high- and low-spin situations. A few issues are addressed concerning the mechanism of alkane hydroxylation, and theoretical support is provided for: (i) the contention that spin inversion occurs along the reaction path, (ii) that the cation radical state of the porphyrin is an essential feature required to accommodate an excess electron from the ferryl moiety and thereby stabilize the ground state of the hydroxylation product, and (iii) that the donor property of the proximal ligand has a significant influence on the energy of the ferryl-to-ring charge-transfer states which are essential to convert the reactant state to the hydroxylation product state. In this sense, our study sheds some light on the difference between the oxidized and reduced HRP forms, HRP(I) and HRP(II), and suggests that the combination of a cation radical porphyrin state and a good π-donor proximal ligand like thiolate, could be the underlying reason for the potent hydroxylation ability of the P-450 ferryl-complex.