Assessment of the Ab Initio Bethe-Salpeter Equation Approach for the Low-Lying Excitation Energies of Bacteriochlorophylls and Chlorophylls

Femtosecond optical studies of the primary charge separation reactions in far-red photosystem II from Synechococcus sp. PCC 7335

Primary processes of light energy conversion by Photosystem II (PSII) were studied using femtosecond broadband pump-probe absorption difference spectroscopy. Transient absorption changes of core complexes isolated from the cyanobacterium Synechococcus sp. PCC 7335 grown under far-red light (FRL-PSII) were compared with the canonical Chl a containing spinach PSII core complexes upon excitation into the red edge of the Qy band. Absorption changes of FRL-PSII were monitored at 278 K in the 400-800 nm spectral range on a timescale of 0.1-500 ps upon selective excitation at 740 nm of four chlorophyll (Chl) f molecules in the light harvesting antenna, or of one Chl d molecule at the ChlD1 position in the reaction center (RC) upon pumping at 710 nm. Numerical analysis of absorption changes and assessment of the energy levels of the presumed ion-radical states made it possible to identify PD1+ChlD1- as the predominant primary charge-separated radical pair, the formation of which upon selective excitation of Chl d has an apparent time of ∼1.6 ps. Electron transfer to the secondary acceptor pheophytin PheoD1 has an apparent time of ∼7 ps with a variety of excitation wavelengths. The energy redistribution between Chl a and Chl f in the antenna occurs within 1 ps, whereas the energy migration from Chl f to the RC occurs mostly with lifetimes of 60 and 400 ps. Potentiometric analysis suggests that in canonical PSII, PD1+ChlD1- can be partially formed from the excited (PD1ChlD1)* state.

Quantum-Quantum and Quantum-Quantum-Classical Schemes for Near-Gap Excitations with Projection-Based-Embedded GW-Bethe-Salpeter Equation

We present quantum-quantum and quantum-quantum-classical schemes based on many-body Green's functions theory in the GW approximation with the Bethe-Salpeter equation (GW-BSE) employing projection-based-embedding (PbE). Such approaches allow defining active and inactive subsystems of larger, complex molecular systems, with only the smaller active subsystem being explicitly treated by GW-BSE offering significant computational advantages. However, as PbE can modify the single-particle states in the Kohn-Sham (KS) ground state calculation and screening effects from the inactive region are not automatically included in GW-BSE, results from such PbE-GW-BSE calculations can deviate from a full-system reference. Here, we scrutinize in detail, e.g., the individual and combined effects of different choices of active regions, the influence of omitting the screening from the inactive region, and strategies for basis set truncation on frontier orbital and near-gap electron-hole excitation energies. As prototypical systems, we consider a diketopyrrolopyrrole bicyclic ring including side-chains, a polarity-sensitive dye (prodan) in aqueous environment, and a π-stacked dimer of benzene and tetracyanoethylene in water, respectively, covering a variety of excitation characters in molecular systems with complex chemical environments and photoinduced processes. Our results suggest that to obtain agreement of approximately 0.1 eV between near-gap excitation energies from embedded and full calculations, the active region should be chosen based on the Mulliken population of the full highest-occupied molecular orbital and that careful benchmarking should be done on the KS level before the actual GW-BSE steps when basis set truncation is used. We find that PbE-GW-BSE offers significant reductions in computation times and, more importantly, memory requirements, making calculations for considerably larger systems tractable.

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.

Investigation on absorption cross-section of photosynthetic pigment molecules based on a mechanistic model of the photosynthetic electron flow-light response in C3, C4 species and cyanobacteria grown under various conditions

Investigation on intrinsic properties of photosynthetic pigment molecules participating in solar energy absorption and excitation, especially their eigen-absorption cross-section (σ ik) and effective absorption cross-section (σ ' ik), is important to understand photosynthesis. Here, we present the development and application of a new method to determine these parameters, based on a mechanistic model of the photosynthetic electron flow-light response. The analysis with our method of a series of previously collected chlorophyll a fluorescence data shows that the absorption cross-section of photosynthetic pigment molecules has different values of approximately 10-21 m2, for several photosynthetic organisms grown under various conditions: (1) the conifer Abies alba Mill., grown under high light or low light; (2) Taxus baccata L., grown under fertilization or non-fertilization conditions; (3) Glycine max L. (Merr.), grown under a CO2 concentration of 400 or 600 μmol CO2 mol-1 in a leaf chamber under shaded conditions; (4) Zea mays L., at temperatures of 30°C or 35°C in a leaf chamber; (5) Osmanthus fragrans Loureiro, with shaded-leaf or sun-leaf; and (6) the cyanobacterium Microcystis aeruginosa FACHB905, grown under two different nitrogen supplies. Our results show that σ ik has the same order of magnitude (approximately 10-21 m2), and σ ' ik for these species decreases with increasing light intensity, demonstrating the operation of a key regulatory mechanism to reduce solar absorption and avoid high light damage. Moreover, compared with other approaches, both σ ik and σ ' ik can be more easily estimated by our method, even under various growth conditions (e.g., different light environment; different CO2, NO2, O2, and O3 concentrations; air temperatures; or water stress), regardless of the type of the sample (e.g., dilute or concentrated cell suspensions or leaves). Our results also show that CO2 concentration and temperature have little effect on σ ik values for G. max and Z. mays. Consequently, our approach provides a powerful tool to investigate light energy absorption of photosynthetic pigment molecules and gives us new information on how plants and cyanobacteria modify their light-harvesting properties under different stress conditions.

Recent progress in atomistic modeling of light-harvesting complexes: a mini review

In this mini review, we focus on recent advances in the atomistic modeling of biological light-harvesting (LH) complexes. Because of their size and sophisticated electronic structures, multiscale methods are required to investigate the dynamical and spectroscopic properties of such complexes. The excitation energies, in this context also known as site energies, excitonic couplings, and spectral densities are key quantities which usually need to be extracted to be able to determine the exciton dynamics and spectroscopic properties. The recently developed multiscale approach based on the numerically efficient density functional tight-binding framework followed by excited state calculations has been shown to be superior to the scheme based on pure classical molecular dynamics simulations. The enhanced approach, which improves the description of the internal vibrational dynamics of the pigment molecules, yields spectral densities in good agreement with the experimental counterparts for various bacterial and plant LH systems. Here, we provide a brief overview of those results and described the theoretical foundation of the multiscale protocol.

Quasiparticle Self-Consistent GW-Bethe-Salpeter Equation Calculations for Large Chromophoric Systems

The GW-Bethe–Salpeter equation (BSE) method is promising for calculating the low-lying excitonic states of molecular systems. However, so far it has only been applied to rather small molecules and in the commonly implemented diagonal approximations to the electronic self-energy, it depends on a mean-field starting point. We describe here an implementation of the self-consistent and starting-point-independent quasiparticle self-consistent (qsGW)-BSE approach, which is suitable for calculations on large molecules. We herein show that eigenvalue-only self-consistency can lead to an unfaithful description of some excitonic states for chlorophyll dimers while the qsGW-BSE vertical excitation energies (VEEs) are in excellent agreement with spectroscopic experiments for chlorophyll monomers and dimers measured in the gas phase. Furthermore, VEEs from time-dependent density functional theory calculations tend to disagree with experimental values and using different range-separated hybrid (RSH) kernels does change the VEEs by up to 0.5 eV. We use the new qsGW-BSE implementation to calculate the lowest excitation energies of the six chromophores of the photosystem II (PSII) reaction center (RC) with nearly 2000 correlated electrons. Using more than 11,000 (6000) basis functions, the calculation could be completed in less than 5 (2) days on a single modern compute node. In agreement with previous TD-DFT calculations using RSH kernels on models that also do not include environmental effects, our qsGW-BSE calculations only yield states with local characters in the low-energy spectrum of the hexameric complex. Earlier works with RSH kernels have demonstrated that the protein environment facilitates the experimentally observed interchromophoric charge transfer. Therefore, future research will need to combine correlation effects beyond TD-DFT with an explicit treatment of environmental electrostatics.

An optimally tuned range-separated hybrid starting point for ab initio GW plus Bethe–Salpeter equation calculations of molecules

The ab initio GW plus Bethe–Salpeter equation (GW-BSE, where G is the one particle Green's function and W is the screened Coulomb interaction) approach has emerged as a leading method for predicting excitations in both solids and molecules with a predictive power contingent upon several factors. Among these factors are the (1) generalized Kohn–Sham eigensystem used to construct the GW self-energy and to solve the BSE and (2) the efficacy and suitability of the Tamm–Dancoff approximation. Here, we present a detailed benchmark study of low-lying singlet excitations from a generalized Kohn–Sham (gKS) starting point based on an optimally tuned range-separated hybrid (OTRSH) functional. We show that the use of this gKS starting point with one-shot G0W0 and G0W0-BSE leads to the lowest mean absolute errors (MAEs) and mean signed errors (MSEs), with respect to high-accuracy reference values, demonstrated in the literature thus far for the ionization potentials of the GW100 benchmark set and for low-lying neutral excitations of Thiel’s set molecules in the gas phase, without the need for self-consistency. The MSEs and MAEs of one-shot G0W0-BSE@OTRSH excitation energies are comparable to or lower than those obtained with other functional starting points after self-consistency. Additionally, we compare these results with linear-response time-dependent density functional theory (TDDFT) calculations and find GW-BSE to be superior to TDDFT when calculations are based on the same exchange-correlation functional. This work demonstrates tuned range-separated hybrids used in combination with GW and GW-BSE can greatly suppress starting point dependence for molecules, leading to accuracy similar to that for higher-order wavefunction-based theories for molecules without the need for costlier iterations to self-consistency.

The atomistic modeling of light-harvesting complexes from the physical models to the computational protocol

The function of light-harvesting complexes is determined by a complex network of dynamic interactions among all the different components: the aggregate of pigments, the protein, and the surrounding environment. Complete and reliable predictions on these types of composite systems can be only achieved with an atomistic description. In the last few decades, there have been important advances in the atomistic modeling of light-harvesting complexes. These advances have involved both the completeness of the physical models and the accuracy and effectiveness of the computational protocols. In this Perspective, we present an overview of the main theoretical and computational breakthroughs attained so far in the field, with particular focus on the important role played by the protein and its dynamics. We then discuss the open problems in their accurate modeling that still need to be addressed. To illustrate an effective computational workflow for the modeling of light harvesting complexes, we take as an example the plant antenna complex CP29 and its H111N mutant.

The Seamless Connection of Local and Collective Excited States in Subsystem Time-Dependent Density Functional Theory

The theoretical understanding of photoinduced processes in multichromophoric systems requires, as an essential ingredient, the possibility of accurately describing their electronically excited states. However, the size of these systems often prohibits the usage of conventional electronic-structure methods, so that often multiscale approaches based on phenomenologically motivated models are employed. In contrast, subsystem time-dependent density functional theory (sTDDFT) allows for a subsystem-based ab initio description of multichromophoric systems and therefore allows for, in principle, an exact description of photoinduced processes. This Perspective aims to outline the theoretical foundations and commonly used practical realizations as well as to illustrate benefits of recent developments and open issues in the field of sTDDFT. Prospective, potential future applications and possible methodological developments are discussed.

Low-Order Scaling Quasiparticle Self-Consistent GW for Molecules

Low-order scaling GW implementations for molecules are usually restricted to approximations with diagonal self-energy. Here, we present an all-electron implementation of quasiparticle self-consistent GW for molecular systems. We use an efficient algorithm for the evaluation of the self-energy in imaginary time, from which a static non-local exchange-correlation potential is calculated via analytical continuation. By using a direct inversion of iterative subspace method, fast and stable convergence is achieved for almost all molecules in the GW100 database. Exceptions are systems which are associated with a breakdown of the single quasiparticle picture in the valence region. The implementation is proven to be starting point independent and good agreement of QP energies with other codes is observed. We demonstrate the computational efficiency of the new implementation by calculating the quasiparticle spectrum of a DNA oligomer with 1,220 electrons using a basis of 6,300 atomic orbitals in less than 4 days on a single compute node with 16 cores. We use then our implementation to study the dependence of quasiparticle energies of DNA oligomers consisting of adenine-thymine pairs on the oligomer size. The first ionization potential in vacuum decreases by nearly 1 electron volt and the electron affinity increases by 0.4 eV going from the smallest to the largest considered oligomer. This shows that the DNA environment stabilizes the hole/electron resulting from photoexcitation/photoattachment. Upon inclusion of the aqueous environment via a polarizable continuum model, the differences between the ionization potentials reduce to 130 meV, demonstrating that the solvent effectively compensates for the stabilizing effect of the DNA environment. The electron affinities of the different oligomers are almost identical in the aqueous environment.