A Simple One-Body Approach to the Calculation of the First Electronic Absorption Band of Water

Analysis of Site Energies and Excitonic Couplings: The Role of Symmetry and Polarization

An excitonic coupling model is developed based on an equation-of-motion coupled cluster combined with the fragment molecular orbital method. The effects of polarization and excitonic coupling on the splitting of quasi-degenerate levels in systems containing multiple chromophores are elucidated on dimers of formaldehyde, water, formic acid, hydrogen fluoride, and carbon monoxide. It is shown that the level structure is mainly determined by the mutual polarization of chromophores and to a lesser extent by the excitonic coupling. The role of symmetry in excitonic coupling in dimers is discussed. The excitonic coupling between all residues in the photoactive yellow protein (PDB: 2PHY) is analyzed.

Detangling Solvatochromic Effects by the Effective Fragment Potential Method

Understanding molecular interactions in complex systems opens avenues for the efficient design of new materials with target properties. Energy decomposition methods provide a means to obtain a detailed picture of intermolecular interactions. This work introduces a molecular modeling approach for decomposing the solvatochromic shifts of the electronic excited states into the contributions of the individual molecular fragments of the environment surrounding the chromophore. The developed approach is implemented for the QM/EFP (quantum mechanics/effective fragment potential) model that provides a rigorous first-principles-based description of the electronic states of the chromophores in complex polarizable environments. On the example of two model systems, water pentamer and hydrated uracil, we show how the decomposition of the solvatochromic shifts into the contributions of individual solvent water molecules provides a detailed picture of the intermolecular interactions in the ground and excited states of these systems. The analysis also demonstrates the nonadditivity of solute-solvent interactions and the significant contribution of solute polarization to the total values of solvatochromic shifts.

Towards a quantitative description of excitonic couplings in photosynthetic pigment-protein complexes: quantum chemistry driven multiscale approaches

A structure-based quantitative calculation of excitonic couplings between photosynthetic pigments has to describe the dynamical polarization of the protein/solvent environment of the pigments, giving rise to reaction field and screening effects. Here, this challenging problem is approached by combining the fragment molecular orbital (FMO) method with the polarizable continuum model (PCM). The method is applied to compute excitonic couplings between chlorophyll a (Chl a) pigments of the water-soluble chlorophyll-binding protein (WSCP). By calibrating the vacuum dipole strength of the 0-0 transition of the Chl a chromophores according to experimental data, an excellent agreement between calculated and experimental linear absorption and circular dichroism spectra of WSCP is obtained. The effect of the mutual polarization of the pigment ground states is calculated to be very small. The simple Poisson-Transition-charge-from-Electrostatic-potential (Poisson-TrEsp) method is found to accurately describe the screening part of the excitonic coupling, obtained with FMO/PCM. Taking into account that the reaction field effects of the latter method can be described by a scalar constant leads to an improvement of Poisson-TrEsp that is expected to provide the basis for simple and realistic calculations of optical spectra and energy transfer in photosynthetic light-harvesting complexes. In addition, we present an expression for the estimation of Huang-Rhys factors of high-frequency pigment vibrations from experimental fluorescence line-narrowing spectra that takes into account the redistribution of oscillator strength by the interpigment excitonic coupling. Application to WSCP results in corrected Huang-Rhys factors that are less than one third of the original values obtained by the standard electronic two-state analysis that neglects the above redistribution. These factors are important for the estimation of the dipole strength of the 0-0 transition of the chromophores and for the development of calculation schemes for the spectral density of the exciton-vibrational coupling.

Cooperation and Environment Characterize the Low-Lying Optical Spectrum of Liquid Water

The optical spectrum of liquid water is analyzed by subsystem time-dependent density functional theory. We provide simple explanations for several important (and so far elusive) features. Due to the disordered environment surrounding each water molecule, the joint density of states of the liquid is much broader than that of the vapor, thus explaining the red-shifted Urbach tail of the liquid compared to the gas phase. Confinement effects provided by the first solvation shell are responsible for the blue shift of the first absorption peak compared to the vapor. In addition, we also characterize many-body excitonic effects. These dramatically affect the spectral weights at low frequencies, contributing to the refractive index by a small but significant amount.

Ab initio calculation of the electronic absorption spectrum of liquid water

The electronic absorption spectrum of liquid water was investigated by coupling a one-body energy decomposition scheme to configurations generated by classical and Born-Oppenheimer Molecular Dynamics (BOMD). A Frenkel exciton Hamiltonian formalism was adopted and the excitation energies in the liquid phase were calculated with the equation of motion coupled cluster with single and double excitations method. Molecular dynamics configurations were generated by different approaches. Classical MD were carried out with the TIP4P-Ew and AMOEBA force fields. The BLYP and BLYP-D3 exchange-correlation functionals were used in BOMD. Theoretical and experimental results for the electronic absorption spectrum of liquid water are in good agreement. Emphasis is placed on the relationship between the structure of liquid water predicted by the different models and the electronic absorption spectrum. The theoretical gas to liquid phase blue-shift of the peak positions of the electronic absorption spectrum is in good agreement with experiment. The overall shift is determined by a competition between the O-H stretching of the water monomer in liquid water that leads to a red-shift and polarization effects that induce a blue-shift. The results illustrate the importance of coupling many-body energy decomposition schemes to molecular dynamics configurations to carry out ab initio calculations of the electronic properties in liquid phase.

Insights into the ultraviolet spectrum of liquid water from model calculations: the different roles of donor and acceptor hydrogen bonds in water pentamers

With a view toward a better understanding of changes in the peak position and shape of the first absorption band of water with condensation or temperature, results from electronic structure calculations using high level wavefunction based and time-dependent density functional methods are reported for water pentamers. Excitation energies, oscillator strengths, and redistributions of electron density are determined for the quasitetrahedral water pentamer in its C(2v) equilibrium geometry and for many pentamer configurations sampled from molecular simulation of liquid water. Excitations associated with surface molecules are removed in order to focus on those states associated with the central molecule, which are the most representative of the liquid environment. The effect of hydrogen bonding on the lowest excited state associated with the central molecule is studied by adding acceptor or donor hydrogen bonds to tetramer and trimer substructures of the C(2v) pentamer, and by sampling liquid-like configurations having increasing number of acceptor or donor hydrogen bonds of the central molecule. Our results provide clear evidence that the blueshift of excitation energies upon condensation is essentially determined by acceptor hydrogen bonds, and the magnitudes of these shifts are determined by the number of such, whereas donor hydrogen bonds do not induce significant shifts in excitation energies. This qualitatively different role of donor and acceptor hydrogen bonds is understood in terms of the different roles of the 1b(1) monomer molecular orbitals, which establishes an intimate connection between the valence hole and excitation energy shifts. Since the valence hole of the lowest excitation associated with the central molecule is found to be well localized in all liquid-like hydrogen bonding environments, with an average radius of gyration of ~1.6 Å that is much lower than the nearest neighbor O-O distance, a clear and unambiguous connection between hydrogen bonding environments and excitation energy shifts can be established. Based on these results, it is concluded that peak position of the first absorption band is mainly determined by the relative distribution of single and double acceptor hydrogen bonding environments, whereas the shape of the first absorption band is mainly determined by the relative distribution of acceptor and broken acceptor hydrogen bonding environments. The temperature dependence of the peak position and shape of the first absorption band can be readily understood in terms of changes to these relative populations.