Optical and radiationless intramolecular electron transitions in nonpolar fluids: Relative effects of induction and dispersion interactions

Reorganization energy of electron transfer

The theory of electron transfer reactions establishes the conceptual foundation for redox solution chemistry, electrochemistry, and bioenergetics. Electron and proton transfer across the cellular membrane provide all energy of life gained through natural photosynthesis and mitochondrial respiration. Rates of biological charge transfer set kinetic bottlenecks for biological energy storage. The main system-specific parameter determining the activation barrier for a single electron-transfer hop is the reorganization energy of the medium. Both harvesting of light energy in natural and artificial photosynthesis and efficient electron transport in biological energy chains require reduction of the reorganization energy to allow fast transitions. This review article discusses mechanisms by which small values of the reorganization energy are achieved in protein electron transfer and how similar mechanisms can operate in other media, such as nonpolar and ionic liquids. One of the major mechanisms of reorganization energy reduction is through non-Gibbsian (nonergodic) sampling of the medium configurations on the reaction time. A number of alternative mechanisms, such as electrowetting of active sites of proteins, give rise to non-parabolic free energy surfaces of electron transfer. These mechanisms, and nonequilibrium population of donor-acceptor vibrations, lead to a universal phenomenology of separation between the Stokes shift and variance reorganization energies of electron transfer.

Electron transfer in nonpolar media

Electron transfer in nonpolar media violates the temperature scaling predicted by the fluctuation–dissipation theorem.

Polarizability of the active site of cytochrome c reduces the activation barrier for electron transfer

Enzymes in biology’s energy chains operate with low energy input distributed through multiple electron transfer steps between protein active sites. The general challenge of biological design is how to lower the activation barrier without sacrificing a large negative reaction free energy. We show that this goal is achieved through a large polarizability of the active site. It is polarized by allowing a large number of excited states, which are populated quantum mechanically by electrostatic fluctuations of the protein and hydration water shells. This perspective is achieved by extensive mixed quantum mechanical/molecular dynamics simulations of the half reaction of reduction of cytochrome c. The barrier for electron transfer is consistently lowered by increasing the number of excited states included in the Hamiltonian of the active site diagonalized along the classical trajectory. We suggest that molecular polarizability, in addition to much studied electrostatics of permanent charges, is a key parameter to consider in order to understand how enzymes work.

Protein electron transfer: is biology (thermo)dynamic?

Simple physical mechanisms are behind the flow of energy in all forms of life. Energy comes to living systems through electrons occupying high-energy states, either from food (respiratory chains) or from light (photosynthesis). This energy is transformed into the cross-membrane proton-motive force that eventually drives all biochemistry of the cell. Life's ability to transfer electrons over large distances with nearly zero loss of free energy is puzzling and has not been accomplished in synthetic systems. The focus of this review is on how this energetic efficiency is realized. General physical mechanisms and interactions that allow proteins to fold into compact water-soluble structures are also responsible for a rugged landscape of energy states and a broad distribution of relaxation times. Specific to a protein as a fluctuating thermal bath is the protein-water interface, which is heterogeneous both dynamically and structurally. The spectrum of interfacial fluctuations is a consequence of protein's elastic flexibility combined with a high density of surface charges polarizing water dipoles into surface nanodomains. Electrostatics is critical to the protein function and the relevant questions are: (i) What is the spectrum of interfacial electrostatic fluctuations? (ii) Does the interfacial biological water produce electrostatic signatures specific to proteins? (iii) How is protein-mediated chemistry affected by electrostatics? These questions connect the fluctuation spectrum to the dynamical control of chemical reactivity, i.e. the dependence of the activation free energy of the reaction on the dynamics of the bath. Ergodicity is often broken in protein-driven reactions and thermodynamic free energies become irrelevant. Continuous ergodicity breaking in a dense spectrum of relaxation times requires using dynamically restricted ensembles to calculate statistical averages. When applied to the calculation of the rates, this formalism leads to the nonergodic activated kinetics, which extends the transition-state theory to dynamically dispersive media. Releasing the grip of thermodynamics in kinetic calculations through nonergodicity provides the mechanism for an efficient optimization between reaction rates and the spectrum of relaxation times of the protein-water thermal bath. Bath dynamics, it appears, play as important role as the free energy in optimizing biology's performance.

Energetics and kinetics of primary charge separation in bacterial photosynthesis

We report the results of molecular dynamics (MD) simulations and formal modeling of the free-energy surfaces and reaction rates of primary charge separation in the reaction center of Rhodobacter sphaeroides. Two simulation protocols were used to produce MD trajectories. Standard force-field potentials were employed in the first protocol. In the second protocol, the special pair was made polarizable to reproduce a high polarizability of its photoexcited state observed by Stark spectroscopy. The charge distribution between covalent and charge-transfer states of the special pair was dynamically adjusted during the simulation run. We found from both protocols that the breadth of electrostatic fluctuations of the protein/water environment far exceeds previous estimates, resulting in about 1.6 eV reorganization energy of electron transfer in the first protocol and 2.5 eV in the second protocol. Most of these electrostatic fluctuations become dynamically frozen on the time scale of primary charge separation, resulting in much smaller solvation contributions to the activation barrier. While water dominates solvation thermodynamics on long observation times, protein emerges as the major thermal bath coupled to electron transfer on the picosecond time of the reaction. Marcus parabolas were obtained for the free-energy surfaces of electron transfer by using the first protocol, while a highly asymmetric surface was obtained in the second protocol. A nonergodic formulation of the diffusion-reaction electron-transfer kinetics has allowed us to reproduce the experimental results for both the temperature dependence of the rate and the nonexponential decay of the population of the photoexcited special pair.

Energetics of electron-transfer reactions in soft condensed media

The coupling of electronic transitions within molecules to condensed-phase media involves a complex hierarchy of spatial and dynamical scales. Thermodynamics of activation is related to the length scale of microscopic interactions reflected in the non-Arrhenius reaction kinetics. Solvent dynamics make a particularly strong impact on the activation barrier when the time scale of the reaction is comparable to the relaxation time of the solvent, and the reaction barrier becomes nonergodic. Finally, molecular polarizability is responsible for complex nonparabolic free energy surfaces for electron transfer. We discuss the application of these ideas to soft condensed solvents such as supercooled liquids, liquid crystals, and photosynthetic reaction centers.

Dynamical Arrest of Electron Transfer in Liquid Crystalline Solvents

We argue that electron transfer reactions in slowly relaxing solvents proceed in the nonergodic regime, making the reaction activation barrier strongly dependent on the solvent dynamics. For typical dielectric relaxation times of polar nematics, electron transfer reactions in the subnanosecond time scale fall into nonergodic regime in which nuclear solvation energies entering the activation barrier are significantly lower than their thermodynamic values. The transition from isotropic to nematic phase results in weak discontinuities of the solvation energies at the transition point and the appearance of solvation anisotropy weakening with increasing solute size. The theory is applied to analyze experimental kinetic data for the electron transfer kinetics in the isotropic phase of 5CB liquid crystalline solvent. We predict that the energy gap law of electron transfer reactions in slowly relaxing solvents is characterized by regions of fast change of the rate at points where the reaction switches between the ergodic and nonergodic regimes. The dependence of the rate on the donor-acceptor separation may also be affected in a way of producing low values for the exponential falloff parameter.

Quadrupolar solvatochromism: 4-amino-phthalimide in toluene

We present calculations of the temperature dependence of the solvent reorganization energy of 4-amino-phthalimide chromophore in quadrupolar toluene. The reorganization energy is a sum of the contributions from quadrupolar and induction solvation. We employ several calculation formalisms in order to evaluate their performance against the experiment. The point-dipole and full atomic distributions of solute charge are compared to show that the point-dipole approximation works well for this chromophore. We also show that most of the reorganization entropy comes from the quadrupolar response. Induction solvation amounts to about 10% of the entropy. Both the reorganization energy and the reorganization entropy are greatly affected by the local solute-solvent density profile (density reorganization) which contributes about half of their values. The induction reorganization energy is strongly affected by the microscopic, nonlocal nature of the density fluctuations of the solvent around the solute.

Equilibrium solvation in quadrupolar solvents

We present a microscopic theory of equilibrium solvation in solvents with zero dipole moment and nonzero quadrupole moment (quadrupolar solvents). The theory is formulated in terms of autocorrelation functions of the quadrupolar polarization (structure factors). It can be therefore applied to an arbitrary dense quadrupolar solvent for which the structure factors are defined. We formulate a simple analytical perturbation treatment for the structure factors. The solute is described by coordinates, radii, and partial charges of constituent atoms. The theory is tested on Monte Carlo simulations of solvation in model quadrupolar solvents. It is also applied to the calculation of the activation barrier of electron transfer reactions in a cleft-shaped donor-bridge-acceptor complex dissolved in benzene with the structure factors of quadrupolar polarization obtained from molecular-dynamics simulations.

Anomalous Microscopic Dielectric Response of Dipolar Solvents and Water

This paper reports measurements of static microscopic dielectric response of several dipolar solvents to charge redistribution in a fluorescent probe. Contrary to recent predictions of dielectric theories and computer simulations of bulk liquids, the observed dielectric response of most solvents conforms to the macroscopic continuum description even at atomic distances, as if these solvents had no spatial intermolecular structure. Such conformance is observed for several probes when the contribution of specific probe-solvent interactions to the response is negligible. However, water, formamide, and glycerol exhibit anomalous responses even though such a probe is used. We discuss a possible reason for the macroscopic-like behavior and a connection between the anomaly and fluctuating structures formed by anomalous solvents near the hydrophobic surface of the probe.

Solvent reorganization energy of electron-transfer reactions in polar solvents

A microscopic theory of solvent reorganization energy in polar molecular solvents is developed. The theory represents the solvent response as a combination of the density and polarization fluctuations of the solvent given in terms of the density and polarization structure factors. A fully analytical formulation of the theory is provided for a solute of arbitrary shape with an arbitrary distribution of charge. A good agreement between the analytical procedure and the results of Monte Carlo simulations of model systems is achieved. The reorganization energy splits into the contributions from density fluctuations and polarization fluctuations. The polarization part is dominated by longitudinal polarization response. The density part is inversely proportional to temperature. The dependence of the solvent reorganization energy on the solvent dipole moment and refractive index is discussed.

Inertial solvent dynamics and the analysis of spectral line shapes: Temperature-dependent absorption spectrum of β-carotene in nonpolar solvent

The influence of solvent dynamics on optical spectra is often described by a stochastic model which assumes exponential relaxation of the time-correlation function for solvent-induced frequency fluctuations. In contrast, theory and experiment suggest that the initial (subpicosecond) phase of solvent relaxation, resulting from inertial motion of the solvent, is a Gaussian function of time. In this work, we employ numerical and analytical calculations to compare the predicted absorption line shapes and the derived solvent reorganization energies obtained from exponential (Brownian oscillator) versus Gaussian (inertial) solvent dynamics. Both models predict motional narrowing as the ratio kappa = Lambda/Delta is increased, where Lambda and Delta are the frequency and variance, respectively, of the solvent-induced frequency fluctuations. However, the motional narrowing limit is achieved at lower values of kappa for the Brownian oscillator model compared to the inertial model. For a given line shape, the derived value of the solvent reorganization energy lambdasolv is only weakly dependent on the solvent relaxation model employed, though different solvent parameters Lambda and Delta are obtained. The two models are applied to the analysis of the temperature-dependent absorption spectrum of beta-carotene in isopentane and CS2. The derived values of lambdasolv using the Gaussian model are found to be in better agreement with the high temperature limit of Delta2/2kBT than are the values obtained using the Brownian oscillator model. In either approach, the solvent reorganization energy is found to increase slightly with temperature as a result of an increase in the variance Delta of the solvent-induced frequency fluctuations.