Successes and failures of time-dependent density functional theory for the low-lying excited states of chlorophylls

Combined Multireference-Multiscale Approach to the Description of Photosynthetic Reaction Centers

A first-principles description of the primary photochemical processes that drive photosynthesis and sustain life on our planet remains one of the grand challenges of modern science. Recent research established that explicit incorporation of protein electrostatics in excited-state calculations of photosynthetic pigments, achieved for example with quantum-mechanics/molecular-mechanics (QM/MM) approaches, is essential for a meaningful description of the properties and function of pigment-protein complexes. Although time-dependent density functional theory has been used productively so far in QM/MM approaches for the study of such systems, this methodology has limitations. Here we pursue for the first time a QM/MM description of the reaction center in the principal enzyme of oxygenic photosynthesis, Photosystem II, using multireference wave function theory for the high-level QM region. We identify best practices and establish guidelines regarding the rational choice of active space and appropriate state-averaging for the efficient and reliable use of complete active space self-consistent field (CASSCF) and the N-electron valence state perturbation theory (NEVPT2) in the prediction of low-lying excited states of chlorophyll and pheophytin pigments. Given that the Gouterman orbitals are inadequate as a minimal active space, we define specific minimal and extended active spaces for the NEVPT2 description of electronic states that fall within the Q and B bands. Subsequently, we apply our multireference-QM/MM protocol to the description of all pigments in the reaction center of Photosystem II. The calculations reproduce the electrochromic shifts induced by the protein matrix and the ordering of site energies consistent with the identity of the primary donor (ChlD1) and the experimentally known asymmetric and directional electron transfer. The optimized protocol sets the stage for future multireference treatments of multiple pigments, and hence for multireference studies of charge separation, while it is transferable to the study of any photoactive embedded tetrapyrrole system.

Exciton interactions of chlorophyll tetramer in water-soluble chlorophyll-binding protein BoWSCP

The exciton interaction of four chlorophyll a (Chl a) molecules in a symmetrical tetrameric complex of the water-soluble chlorophyll-binding protein BoWSCP was analyzed in the pH range of 3-11. Exciton splitting ΔE = 232 ± 2 cm-1 of the Qy band of Chl a into two subcomponents with relative intensities of 78.1 ± 0.7 % and 21.9 ± 0.7 % was determined by a joint decomposition of the absorption and circular dichroism spectra into Gaussian functions. The exciton coupling parameters were calculated based on the BoWSCP atomic structure in three approximations: the point dipole model, the distributed atomic monopoles, and direct ab initio calculations in the TDDFT/PCM approximation. The Coulomb interactions of monomers were calculated within the continuum model using three values of optical permittivity. The models based on the properties of free Chl a in solution suffer from significant errors both in estimating the absolute value of the exciton interaction and in the relative intensity of exciton transitions. Calculations within the TDDFT/PCM approximation reproduce the experimentally determined parameters of the exciton splitting and the relative intensities of the exciton bands. The following factors of pigment-protein and pigment-pigment interactions were examined: deviation of the macrocycle geometry from the planar conformation of free Chl; the formation of hydrogen bonds between the macrocycle and water molecules; the overlap of wave functions of monomers at close distances. The most significant factor is the geometrical deformation of the porphyrin macrocycle, which leads to an increase in the dipole moment of Chl monomer from 5.5 to 6.9 D and to a rotation of the dipole moment by 15° towards the cyclopentane ring. The contributions of resonant charge-transfer states to the wave functions of the Chl dimer were determined and the transition dipole moments of the symmetric and antisymmetric charge-transfer states were estimated.

Physicochemical, antioxidant properties of carotenoids and its optoelectronic and interaction studies with chlorophyll pigments

The physicochemical and antioxidant properties of seven carotenoids: antheraxanthin, β-carotene, neoxanthin, peridinin, violaxanthin, xanthrophyll and zeaxanthin were studied by theoretical means. Then the Optoelectronic properties and interaction of chlorophyll-carotenoid complexes are analysed by TDDFT and IGMPLOT. Global reactivity descriptors for carotenoids and chlorophyll (Chla, Chlb) are calculated via conceptual density functional theory (CDFT). The higher HOMO-LUMO (HL) gap indicated structural stability of carotenoid, chlorophyll and chlorophyll-carotenoid complexes. The chemical hardness for carotenoids and Chlorophyll is found to be lower in the solvent medium than in the gas phase. Results showed that carotenoids can be used as good reactive nucleophile due to lower µ and ω. As proton affinities (PAs) are much lower than the bond dissociation enthalpies (BDEs), it is anticipated that direct antioxidant activity in these carotenoids is mainly due to the sequential proton loss electron transfer (SPLET) mechanism with dominant solvent effects. Also lower PAs of carotenoid suggest that antioxidant activity by the SPLET mechanism should be a result of a balance between proclivities to transfer protons. Reaction rate constant with Transition-State Theory (TST) were estimated for carotenoid-Chlorophyll complexes in gas phase. Time dependent Density Functional Theory (TDDFT) showed that all the chlorophyll (Chla, Chlb)-carotenoid complexes show absorption wavelength in the visible region. The lower S1-T1 adiabatic energy gap indicated ISC transition from S1 to T1 state.

Protein Matrix Control of Reaction Center Excitation in Photosystem II

Photosystem II (PSII) is a multisubunit pigment-protein complex that uses light-induced charge separation to power oxygenic photosynthesis. Its reaction center chromophores, where the charge transfer cascade is initiated, are arranged symmetrically along the D1 and D2 core polypeptides and comprise four chlorophyll (PD1, PD2, ChlD1, ChlD2) and two pheophytin molecules (PheoD1 and PheoD2). Evolution favored productive electron transfer only via the D1 branch, with the precise nature of primary excitation and the factors that control asymmetric charge transfer remaining under investigation. Here we present a detailed atomistic description for both. We combine large-scale simulations of membrane-embedded PSII with high-level quantum-mechanics/molecular-mechanics (QM/MM) calculations of individual and coupled reaction center chromophores to describe reaction center excited states. We employ both range-separated time-dependent density functional theory and the recently developed domain based local pair natural orbital (DLPNO) implementation of the similarity transformed equation of motion coupled cluster theory with single and double excitations (STEOM-CCSD), the first coupled cluster QM/MM calculations of the reaction center. We find that the protein matrix is exclusively responsible for both transverse (chlorophylls versus pheophytins) and lateral (D1 versus D2 branch) excitation asymmetry, making ChlD1 the chromophore with the lowest site energy. Multipigment calculations show that the protein matrix renders the ChlD1 → PheoD1 charge-transfer the lowest energy excitation globally within the reaction center, lower than any pigment-centered local excitation. Remarkably, no low-energy charge transfer states are located within the "special pair" PD1-PD2, which is therefore excluded as the site of initial charge separation in PSII. Finally, molecular dynamics simulations suggest that modulation of the electrostatic environment due to protein conformational flexibility enables direct excitation of low-lying charge transfer states by far-red light.

Accurate Computation of the Absorption Spectrum of Chlorophyll a with Pair Natural Orbital Coupled Cluster Methods

The ability to accurately compute low-energy excited states of chlorophylls is critically important for understanding the vital roles they play in light harvesting, energy transfer, and photosynthetic charge separation. The challenge for quantum chemical methods arises both from the intrinsic complexity of the electronic structure problem and, in the case of biological models, from the need to account for protein–pigment interactions. In this work, we report electronic structure calculations of unprecedented accuracy for the low-energy excited states in the Q and B bands of chlorophyll a. This is achieved by using the newly developed domain-based local pair natural orbital (DLPNO) implementation of the similarity transformed equation of motion coupled cluster theory with single and double excitations (STEOM-CCSD) in combination with sufficiently large and flexible basis sets. The results of our DLPNO–STEOM-CCSD calculations are compared with more approximate approaches. The results demonstrate that, in contrast to time-dependent density functional theory, the DLPNO–STEOM-CCSD method provides a balanced performance for both absorption bands. In addition to vertical excitation energies, we have calculated the vibronic spectrum for the Q and B bands through a combination of DLPNO–STEOM-CCSD and ground-state density functional theory frequency calculations. These results serve as a basis for comparison with gas-phase experiments.

Interactions of copper(ii) and zinc(ii) with chlorophyll: insights from density functional theory studies

The most favored reaction of chlorophyll is computed to be substitution for Cu²⁺ and peripheral chelation for Zn²⁺.

Putting David Craig’s Legacy to Work in Nanotechnology and Biotechnology

David Craig (1919–2015) left us with a lasting legacy concerning basic understanding of chemical spectroscopy and bonding. This is expressed in terms of some of the recent achievements of my own research career, with a focus on integration of Craig’s theories with those of Noel Hush to solve fundamental problems in photosynthesis, molecular electronics (particularly in regard to the molecules synthesized by Maxwell Crossley), and self-assembled monolayer structure and function. Reviewed in particular is the relation of Craig’s legacy to: the 50-year struggle to assign the visible absorption spectrum of arguably the world’s most significant chromophore, chlorophyll; general theories for chemical bonding and structure extending Hush’s adiabatic theory of electron-transfer processes; inelastic electron-tunnelling spectroscopy (IETS); chemical quantum entanglement and the Penrose–Hameroff model for quantum consciousness; synthetic design strategies for NMR quantum computing; Gibbs free-energy measurements and calculations for formation and polymorphism of organic self-assembled monolayers on graphite surfaces from organic solution; and understanding the basic chemical processes involved in the formation of gold surfaces and nanoparticles protected by sulfur-bound ligands, ligands whose form is that of Au0-thiyl rather than its commonly believed AuI-thiolate tautomer.

Assignment of the Q-bands of the chlorophylls: coherence loss via Qx - Qy mixing

We provide a new and definitive spectral assignment for the absorption, emission, high-resolution fluorescence excitation, linear dichroism, and/or magnetic circular dichroism spectra of 32 chlorophyllides in various environments. This encompases all data used to justify previous assignments and provides a simple interpretation of unexplained complex decoherence phenomena associated with Q_x → Q_y relaxation. Whilst most chlorophylls conform to the Gouterman model and display two independent transitions Q_x (S₂) and Q_y (S₁), strong vibronic coupling inseparably mixes these states in chlorophyll-a. This spreads x-polarized absorption intensity over the entireQ-band system to influence all exciton-transport, relaxation and coherence properties of chlorophyll-based photosystems. The fraction of the total absorption intensity attributed to Q_x ranges between 7% and 33%, depending on chlorophyllide and coordination, and is between 10% and 25% for chlorophyll-a. CAM-B3LYP density-functional-theory calculations of the band origins, relative intensities, vibrational Huang-Rhys factors, and vibronic coupling strengths fully support this new assignment.

Density functional theory for charge transfer: the nature of the N-bands of porphyrins and chlorophylls revealed through CAM-B3LYP, CASPT2, and SAC-CI calculations

While density functional theory (DFT) has been proven to be extremely useful for the prediction of thermodynamic and spectroscopic properties of molecules, to date most functionals used in common implementations of DFT display a systematic failure to predict the properties of charge-transfer processes. While this is explicitly manifest in Rydberg transitions of atoms and molecules and in molecular charge-transfer spectroscopy, it also becomes critical for systems containing extended conjugation such as polyenes and other conducting polymers, porphyrins, chlorophylls, etc. A new density functional, a Coulomb-attenuated hybrid exchange-correlation functional (CAM-B3LYP), has recently been developed specifically to overcome these limitations, and it has been shown to properly predict molecular charge-transfer spectra. Here, we demonstrate that it predicts qualitatively reasonable spectra for porphyrin, some oligoporphyrins, and chlorophyll. However, alternate density functionals developed to overcome the same limitations such as current-density functional theory are shown, in their present implementation, to remain inadequate. The CAM-B3LYP results are shown to be in excellent agreement with complete-active-space plus second-order Møller-Plesset perturbation theory and symmetry-adapted cluster configuration interaction calculations: These depict the N and higher bands of porphyrins and chlorophylls as being charge-transfer bands associated with localization of molecular orbitals on individual pyrrole rings. The validity of the basic Gouterman model for the spectra of porphyrins and chlorophylls is confirmed, rejecting modern suggestions that non-Gouterman transitions lie close in energy to the Q-bands of chlorophylls. As porphyrins and chlorophylls provide useful paradigms for problems involving extended conjugation, the results obtained suggest that many significant areas of nanotechnology and biotechnology may now be realistically treated by cost-effective density-functional-based computational methods.

Assignment of the Qy absorption spectrum of photosystem-I from Thermosynechococcus elongatus based on CAM-B3LYP calculations at the PW91-optimized protein structure

The Qy absorption spectrum of Photosystem-I from Thermosynecochoccus elongatus (formerly Synecochoccus elongatus) is calculated using the CAM-B3LYP density functional and INDO schemes based on a quantum-mechanically refined structure for the entire photosystem obtained using the PW91 density functional. These methods present a priori predictions of the absorption and linear dichroism spectra and include protein electrostatic effects, short range inductive effects, long-range and short-range exciton couplings, and superexchange effects involving aromatic residues and carotenes. CAM-B3LYP is used as it is the only known density functional that correctly describes the Q bands of chlorophylls, all other methods contaminating them with erroneous charge-transfer excitations. A critical feature is found to be the use of fully optimized heavy-atom coordinates, with those obtained from just X-ray crystallography providing a poor description of the electronic properties of the chromophores. The result is a realistic first-principles prediction of the observed absorption band that identifies the nature of the red-shifted chlorophylls as well as the energies of the reaction-center chlorophylls and the exciton couplings acting between them. The "special pair" appears more like a dimer of dimers than a self-contained functional unit, with the exciton couplings between its members and the accessory chlorophylls exceeding the internal coupling.

Pressure-induced spectral changes for the special-pair radical cation of the bacterial photosynthetic reaction center

The effect of pressure up to 6 kbars on the near to mid infrared absorption spectrum (7500-14,300 cm(-1) or 1333-700 nm) of the oxidized reaction center of Rhodobacter sphaeroides is measured and interpreted using density-functional B3LYP, INDO, and PM5 calculations. Two weak electronic transition origins at approximately 8010 and approximately 10,210 cm(-1) are unambiguously identified. The first transition is assigned to a Qy tripdoublet band that involves, in the localized description of the excitation, a triplet absorption on one of the bacteriochlorophyll molecules (PM) in the reaction center's special pair intensified by the presence of a radical cation on the other (PL). While most chlorophyll transition energies decrease significantly with increasing pressure, the tripdoublet band is found to be almost pressure insensitive. This difference is attributed to the additional increase in the tripdoublet-band energy accompanying compression of the pi-stacked special pair. The second band could either be the anticipated second Qy tripdoublet state, a Qx tripdoublet state, or a state involving excitation from a low-lying doubly occupied orbital to the half-occupied cationic orbital. A variety of absorption bands that are also resolved in the 8300-9600 cm(-1) region are assigned as vibrational structure associated with the first tripdoublet absorption. These sidebands are composites that are shown by the calculations to comprise many unresolved individual modes; while the calculated pressure sensitivity of each individual mode is small, the calculated pressure dependence of the combined sideband structure is qualitatively similar to the observed pressure dependence, preventing the positive identification of possible additional electronic transitions in this spectral region.

The Green's Function Density Functional Tight-Binding (gDFTB) Method for Molecular Electronic Conduction

A review is presented of the nonequilibrium Green's function (NEGF) method "gDFTB" for evaluating elastic and inelastic conduction through single molecules employing the density functional tight-binding (DFTB) electronic structure method. This focuses on the possible advantages that DFTB implementations of NEGF have over conventional methods based on density functional theory, including not only the ability to treat large irregular metal-molecule junctions with high nonequilibrium thermal distributions but perhaps also the ability to treat dispersive forces, bond breakage, and open-shell systems and to avoid large band lineup errors. New results are presented indicating that DFTB provides a useful depiction of simple gold-thiol interactions. Symmetry is implemented in DFTB, and the advantages it brings in terms of large savings of computational resources with significant increase in numerical stability are described. The power of DFTB is then harnessed to allow the use of gDFTB as a real-time tool to discover the nature of the forces that control inelastic charge transport through molecules and the role of molecular symmetry in determining both elastic and inelastic transport. Future directions for the development of the method are discussed.

Quantum chemical insights in energy dissipation and carotenoid radical cation formation in light harvesting complexes

Light harvesting complexes (LHCs) have been identified in all photosynthetic organisms. To understand their function in light harvesting and energy dissipation, detailed knowledge about possible excitation energy transfer (EET) and electron transfer (ET) processes in these pigment proteins is of prime importance. This again requires the study of electronically excited states of the involved pigment molecules, in LHCs of chlorophylls and carotenoids. This paper represents a critical review of recent quantum chemical calculations on EET and ET processes between pigment pairs relevant for the major LHCs of green plants (LHC-II) and of purple bacteria (LH2). The theoretical methodology for a meaningful investigation of such processes is described in detail, and benefits and limitations of standard methods are discussed. The current status of excited state calculations on chlorophylls and carotenoids is outlined. It is focused on the possibility of EET and ET in the context of chlorophyll fluorescence quenching in LHC-II and carotenoid radical cation formation in LH2. In the context of non-photochemical quenching of green plants, it is shown that replacement of the carotenoid violaxanthin by zeaxanthin in its binding pocket of LHC-II can not result in efficient quenching. In LH2, our computational results give strong evidence that the S(1) states of the carotenoids are involved in carotenoid cation formation. By comparison of theoretical findings with recent experimental data, a general mechanism for carotenoid radical cation formation is suggested.

The role of axial ligands for the structure and function of chlorophylls

We have studied the effect of axial ligation of chlorophyll and bacteriochlorophyll using density functional calculations. Eleven different axial ligands have been considered, including models of histidine, aspartate/glutamate, asparagine/glutamine, serine, tyrosine, methionine, water, the protein backbone, and phosphate. The native chlorophylls, as well as their cation and anion radical states and models of the reaction centres P680 and P700, have been studied and we have compared the geometries, binding energies, reduction potentials, and absorption spectra. Our results clearly show that the chlorophylls strongly prefer to be five-coordinate, in accordance with available crystal structures. The axial ligands decrease the reduction potentials, so they cannot explain the high potential of P680. They also redshift the Q band, but not enough to explain the occurrence of red chlorophylls. However, there is some relation between the axial ligands and their location in the various photosynthetic proteins. In particular, the intrinsic reduction potential of the second molecule in the electron transfer path is always lower than that of the third one, a feature that may prevent back-transfer of the electron.