Transient Terahertz Stark Effect: A Dynamic Probe of Electric Interactions in Polar Liquids.

Electric forces acting on molecules in liquids at ambient temperature fluctuate at terahertz (THz) frequencies with a direct impact on their electronic and optical properties. We introduce the transient THz Stark effect to modify the electronic absorption spectra of dye molecules and, thus, elucidate and determine the underlying molecular interactions and dynamics. Picosecond electric fields of megavolts/cm induce a nonequilibrium response of the prototypical Betaine-30 in polar solution that is probed via transient absorption changes. The field-induced broadening of the absorption band follows the THz intensity in time, with a minor impact of solvent dynamics. The ground and excited state dipole energies in the THz field govern this response, allowing for a quantification of electric forces in a structurally frozen molecular environment.

Short-Range Cooperative Slow-down of Water Solvation Dynamics Around SO4 2--Mg2+ Ion Pairs.

The presence of ions affects the structure and dynamics of water on a multitude of length and time scales. In this context, pairs of Mg2+ and SO4 2- ions in water constitute a prototypical system for which conflicting pictures of hydration geometries and dynamics have been reported. Key issues are the molecular pair and solvation shell geometries, the spatial range of electric interactions, and their impact on solvation dynamics. Here, we introduce asymmetric SO4 2- stretching vibrations as new and most specific local probes of solvation dynamics that allow to access ion hydration dynamics at the dilute concentration (0.2 M) of a native electrolyte environment. Highly sensitive heterodyne 2D-IR spectroscopy in the fingerprint region of the SO4 2- ions around 1100 cm-1 reveals a specific slow-down of solvation dynamics for hydrated MgSO4 and for Na2SO4 in the presence of Mg2+ ions, which manifests as a retardation of spectral diffusion compared to aqueous Na2SO4 solutions in the absence of Mg2+ ions. Extensive molecular dynamics and density functional theory QM/MM simulations provide a microscopic view of the observed ultrafast dephasing and hydration dynamics. They suggest a molecular picture where the slow-down of hydration dynamics arises from the structural peculiarities of solvent-shared SO4 2--Mg2+ ion pairs.

Field-Induced Electron Generation in Water: Solvation Dynamics and Many-Body Interactions.

The solvated electron represents an elementary quantum system in a liquid environment. Electrons solvated in water have raised strong interest because of their prototypical properties, their role in radiation chemistry, and their relevance for charge separation and transport. Nonequilibrium dynamics of photogenerated electrons in water occur on ultrafast time scales and include charge transfer, localization, and energy dissipation processes. We present new insight into the role of fluctuating electric fields of the liquid for generating electrons in the presence of an external terahertz field and address polaronic many-body properties of solvated electrons. This Perspective combines a review of recent results from experiment and theory with a discussion of basic electric interactions of electrons in water.

Coherent polaron dynamics of electrons solvated in polar liquids.

An electron solvated in a polar liquid is an elementary quantum system with properties governed by electric interactions with a fluctuating molecular environment. In the prevailing single particle picture, the quantum ground and excited states are determined by a self-consistent potential, as defined by the particular local configuration of the solvation shell. This description neglects collective many-body excitations, which arise from the coupling of electronic degrees of freedom and nuclear motions of the environment. While recent experiments have demonstrated collective nonequilbrium electronic-nuclear motion, i.e. polaron excitations in liquid water, their relevance in the broader context of polar liquids has remained unexplored. Here, we study the nonequilibrium dielectric response of the, compared to water, less polar alcohols isopropanol, and ethylene glycol, that also display a different hydrogen bond pattern. We demonstrate that ultrafast relaxation of photogenerated electrons impulsively induces coherent charge oscillations, which persist for some 10 ps. They emit electric waves in a frequency range from 0.1 to 2 THz, depending on electron concentration. Oscillation frequencies and line shapes are reproduced by a unified polaron picture for alcohols and water, which is based on a Clausius-Mossotti local field approach for the THz dielectric function. The analysis suggests a longitudinal character of many-body polaron excitations and a weak coupling to transverse excitations, supported by the underdamped character of charge oscillations. Polaron dynamics are governed by the long-range Coulomb interaction between an excess electron and several thousands of polar solvent molecules, while local electron solvation geometries play a minor role.

Electron Transfer-Induced Active Site Structural Relaxation in 64-Photolyase of Drosophila melanogaster.

While catalytic electron flow and photoreactivation of CPD-photolyases are increasingly understood, the microscopic details of the 64-photolyase repair mechanism are perpetually debated. Here, we investigate in long-time (µs) molecular dynamics simulations combined with extensive quantum mechanical/molecular mechanical (QM/MM) simulations the primary electron transfer (ET) reactions in 64-photolyase of Drosophila melanogaster (D. melanogaster). The characterization of the relative energetics of locally excited and charge separated states in the (6-4) photoproduct enzyme repair complex reveals a charge-separated state involving the adenine moiety of the FADH- cofactor that facilitates reduction of the photoproduct. Microscopic details of the collective reaction coordinate of ET reactions are identified that involve the reorganization of the hydrogen bond network and structural relaxation of the active site. The simulations reveal complex active site relaxation dynamics involving distinguished amino acids (Lys246, His365, and His369), conformational reorganization of the hydroxyl group of the (6-4) photoproduct, and a strengthening of hydrogen bonds with immobilized water molecules. In particular, rotation of the Lys246 side chain is found to impose a double-well character along the reaction coordinate of the ET reaction. Our findings suggest that the primary ET reactions in the (6-4) photoproduct enzyme repair complex of D. melanogaster are governed by a complex multi-minima active site relaxation dynamics and potentially precede the equilibration of the protein. ET pathways mediated by the adenine moiety and the 5' side of the photoproduct are proposed to be relevant for triggering the catalytic (6-4) photoproduct reactivation.

Infrared and NMR Spectroscopic Fingerprints of the Asymmetric H7 + O3 Complex in Solution.

The front cover artwork is provided by the groups of Prof. Ehud Pines (BGU, Israel) and Dr. Benjamin Fingerhut (MBI, Berlin). The image shows a scientist integrating experiments with theory for resolving the structural diffusion of the aqueous proton in acetonitrile providing a novel view on the Grotthuss mechanism. Read the full text of the Article at 10.1002/cphc.202001046.

Phosphate Vibrations Probe Electric Fields in Hydrated Biomolecules: Spectroscopy, Dynamics, and Interactions.

Electric interactions have a strong impact on the structure and dynamics of biomolecules in their native water environment. Given the variety of water arrangements in hydration shells and the femto- to subnanosecond time range of structural fluctuations, there is a strong quest for sensitive noninvasive probes of local electric fields. The stretching vibrations of phosphate groups, in particular the asymmetric (PO2)- stretching vibration νAS(PO2)-, allow for a quantitative mapping of dynamic electric fields in aqueous environments via a field-induced redshift of their transition frequencies and concomitant changes of vibrational line shapes. We present a systematic study of νAS(PO2)- excitations in molecular systems of increasing complexity, including dimethyl phosphate (DMP), short DNA and RNA duplex structures, and transfer RNA (tRNA) in water. A combination of linear infrared absorption, two-dimensional infrared (2D-IR) spectroscopy, and molecular dynamics (MD) simulations gives quantitative insight in electric-field tuning rates of vibrational frequencies, electric field and fluctuation amplitudes, and molecular interaction geometries. Beyond neat water environments, the formation of contact ion pairs of phosphate groups with Mg2+ ions is demonstrated via frequency upshifts of the νAS(PO2)- vibration, resulting in a distinct vibrational band. The frequency positions of contact geometries are determined by an interplay of attractive electric and repulsive exchange interactions.

Terahertz Polaron Oscillations of Electrons Solvated in Liquid Water.

The terahertz (THz) response of solvated electrons in liquid water is studied in nonlinear ultrafast pump-probe experiments. Free electrons with concentrations from c_{e}=4 to 140×10^{-6} moles/liter are generated by high-field THz or near-infrared multiphoton excitation. The time-resolved change of the dielectric function as mapped by broadband THz pulses exhibits pronounced oscillations persisting up to 30 ps. Their frequency increases with electron concentration from 0.2 to 1.5 THz. The oscillatory response is assigned to impulsively excited coherent polarons involving coupled electron and water shell motions with a frequency set by the local electric field.

Infrared and NMR Spectroscopic Fingerprints of the Asymmetric H7+ O3 Complex in Solution.

Infrared (IR) absorption in the 1000-3700 cm-1 range and 1 H NMR spectroscopy reveal the existence of an asymmetric protonated water trimer, H7+ O3, in acetonitrile. The core H7+ O3 motif persists in larger protonated water clusters in acetonitrile up to at least 8 water molecules. Quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations reveal irreversible proton transport promoted by propagating the asymmetric H7+ O3 structure in solution. The QM/MM calculations allow for the successful simulation of the measured IR absorption spectra of H7+ O3 in the OH stretch region, which reaffirms the assignment of the H7+ O3 spectra to a hybrid-complex structure: a protonated water dimer strongly hydrogen-bonded to a third water molecule with the proton exchanging between the two possible shared-proton Zundel-like centers. The H7+ O3 structure lends itself to promoting irreversible proton transport in presence of even one additional water molecule. We demonstrate how continuously evolving H7+ O3 structures may support proton transport within larger water solvates.

Magnesium Contact Ions Stabilize the Tertiary Structure of Transfer RNA: Electrostatics Mapped by Two-Dimensional Infrared Spectra and Theoretical Simulations.

Ions interacting with hydrated RNA play a central role in defining its secondary and tertiary structure. While spatial arrangements of ions, water molecules, and phosphate groups have been inferred from X-ray studies, the role of electrostatic and other noncovalent interactions in stabilizing compact folded RNA structures is not fully understood at the molecular level. Here, we demonstrate that contact ion pairs of magnesium (Mg2+) and phosphate groups embedded in local water shells stabilize the tertiary equilibrium structure of transfer RNA (tRNA). Employing dialyzed tRNAPhe from yeast and tRNA from Escherichia coli, we follow the population of Mg2+ sites close to phosphate groups of the ribose-phosphodiester backbone step by step, combining linear and nonlinear infrared spectroscopy of phosphate vibrations with molecular dynamics simulations and ab initio vibrational frequency calculations. The formation of up to six Mg2+/phosphate contact pairs per tRNA and local field-induced reorientations of water molecules balance the phosphate-phosphate repulsion in nonhelical parts of tRNA, thus stabilizing the folded structure electrostatically. Such geometries display limited sub-picosecond fluctuations in the arrangement of water molecules and ion residence times longer than 1 µs. At higher Mg2+ excess, the number of contact ion pairs per tRNA saturates around 6 and weakly interacting ions prevail. Our results suggest a predominance of contact ion pairs over long-range coupling of the ion atmosphere and the biomolecule in defining and stabilizing the tertiary structure of tRNA.

On the role of non-diagonal system-environment interactions in bridge-mediated electron transfer.

Bridge-mediated electron transfer (ET) between a donor and an acceptor is prototypical for the description of numerous most important ET scenarios. While multi-step ET and the interplay of sequential and direct superexchange transfer pathways in the donor-bridge-acceptor (D-B-A) model are increasingly understood, the influence of off-diagonal system-bath interactions on the transfer dynamics is less explored. Off-diagonal interactions account for the dependence of the ET coupling elements on nuclear coordinates (non-Condon effects) and are typically neglected. Here, we numerically investigate with quasi-adiabatic propagator path integral simulations the impact of off-diagonal system-environment interactions on the transfer dynamics for a wide range of scenarios in the D-B-A model. We demonstrate that off-diagonal system-environment interactions can have profound impact on the bridge-mediated ET dynamics. In the considered scenarios, the dynamics itself does not allow for a rigorous assignment of the underlying transfer mechanism. Furthermore, we demonstrate how off-diagonal system-environment interaction mediates anomalous localization by preventing long-time depopulation of the bridge B and how coherent transfer dynamics between donor D and acceptor A can be facilitated. The arising non-exponential short-time dynamics and coherent oscillations are interpreted within an equivalent Hamiltonian representation of a primary reaction coordinate model that reveals how the complex vibronic interplay of vibrational and electronic degrees of freedom underlying the non-Condon effects can impose donor-to-acceptor coherence transfer on short timescales.

Regulatory Impact of the C-Terminal Tail on Charge Transfer Pathways in Drosophila Cryptochrome.

Interconnected transcriptional and translational feedback loops are at the core of the molecular mechanism of the circadian clock. Such feedback loops are synchronized to external light entrainment by the blue light photoreceptor cryptochrome (CRY) that undergoes conformational changes upon light absorption by an unknown photoexcitation mechanism. Light-induced charge transfer (CT) reactions in Drosophila CRY (dCRY) are investigated by state-of-the-art simulations that reveal a complex, multi-redox site nature of CT dynamics on the microscopic level. The simulations consider redox-active chromophores of the tryptophan triad (Trp triad) and further account for pathways mediated by W314 and W422 residues proximate to the C-terminal tail (CTT), thus avoiding a pre-bias to specific W-mediated CT pathways. The conducted dissipative quantum dynamics simulations employ microscopically derived model Hamiltonians and display complex and ultrafast CT dynamics on the picosecond timescale, subtly balanced by the electrostatic environment of dCRY. In silicio point mutations provide a microscopic basis for rationalizing particular CT directionality and demonstrate the degree of electrostatic control realized by a discrete set of charged amino acid residues. The predicted participation of CT states in proximity to the CTT relates the directionality of CT reactions to the spatial vicinity of a linear interaction motif. The results stress the importance of CTT directional charge transfer in addition to charge transfer via the Trp triad and call for the use of full-length CRY models including the interactions of photolyase homology region (PHR) and CTT domains.

Field-Induced Tunneling Ionization and Terahertz-Driven Electron Dynamics in Liquid Water.

Liquid water at ambient temperature displays ultrafast molecular motions and concomitant fluctuations of very strong electric fields originating from the dipolar H2O molecules. We show that such random intermolecular fields induce the tunnel ionization of water molecules, which becomes irreversible if an external terahertz (THz) pulse imposes an additional directed electric field on the liquid. Time-resolved nonlinear THz spectroscopy maps charge separation, transport, and localization of the released electrons on a few-picosecond time scale. The highly polarizable localized electrons modify the THz absorption spectrum and refractive index of water, a manifestation of a highly nonlinear response. Our results demonstrate how the interplay of local electric field fluctuations and external electric fields allows for steering charge dynamics and dielectric properties in aqueous systems.

Change of Hydration Patterns upon RNA Melting Probed by Excitations of Phosphate Backbone Vibrations.

The water hydration shell has a decisive impact on the structural and functional properties of RNA. Changes of the RNA structure upon melting and in biochemical processes are accompanied by a change of hydration patterns, a process that is barely characterized. To discern hydration geometries around the backbone phosphate groups of an RNA double helix at the molecular level, we combine two-dimensional infrared spectroscopy of phosphate vibrations with theoretical simulations. There are three distinct coexisting hydration motifs of the RNA A-helix: an ordered chain-like arrangement of water molecules with links between neighboring phosphate groups, separate local hydration shells of up to six water molecules, and hydrated phosphate/counterion contact pairs. RNA disordering upon melting is connected with a transition from predominantly ordered water structures to local hydration shells around phosphate units. Structural fluctuations are dominated by librational water motions occurring on a 300 fs time scale, without exchange between hydration motifs.

Contact Ion Pairs of Phosphate Groups in Water: Two-Dimensional Infrared Spectroscopy of Dimethyl Phosphate and ab Initio Simulations.

The interaction of phosphate groups with ions in an aqueous environment has a strong impact on the structure and folding processes of DNA and RNA. The dynamic variety of ionic arrangements, including both contact pairs and water separated ions, and the molecular coupling mechanisms are far from being understood. In a combined experimental and theoretical approach, we address the properties of contact ion pairs of the prototypical system dimethyl phosphate with Na+, Ca2+, and Mg2+ ions in water. Linear and femtosecond two-dimensional infrared (2D-IR) spectroscopy of the asymmetric (PO2)- stretching vibration separates and characterizes the different species via their blue-shifted vibrational signatures and 2D-IR line shapes. Phosphate-magnesium contact pairs stand out as the most compact geometry while the contact pairs with Ca2+ and Na+ display a wider structural variation. Microscopic density functional theory simulations rationalize the observed frequency shifts and reveal distinct differences between the contact geometries.

Hydrated Excess Protons in Acetonitrile/Water Mixtures: Solvation Species and Ultrafast Proton Motions.

The solvation structure of protons in aqueous media is highly relevant to electric properties and to proton transport in liquids and membranes. At ambient temperature, polar liquids display structural fluctuations on femto- to picosecond time scales with a direct impact on proton solvation. We use two-dimensional infrared (2D-IR) spectroscopy to follow proton dynamics in acetonitrile/water mixtures with the Zundel cation H5O2+ prepared in neat acetonitrile as a benchmark. The 2D-IR spectra of the proton transfer mode of H5O2+ demonstrate stochastic large-amplitude motions in the double-minimum proton potential, driven by fluctuating electric fields. In all cases, the excess proton is embedded in a water dimer, forming an H5O2+ complex as the major solvation species. This observation is rationalized by quantum mechanics/molecular mechanics molecular dynamics simulations including up to four water molecules embedded in acetonitrile. The Zundel motif interacts with its closest water neighbor in an H7O3+ unit without persistent proton localization.

Coupled excitation energy and charge transfer dynamics in reaction centre inspired model systems.

Functional operation conditions of reaction centre core complexes require the tight coupling of exciton states to concomitant charge separation. Rigorous theoretical treatment of such integrated excitation energy transfer (EET) and charge transfer (CT) dynamics is particularly challenging due to (i) appreciable system sizes, (ii) inter-site and system-bath couplings of similar magnitude that render the Born-Markov approximation invalid, (iii) substantial reorganization energies of CT states, and (iv) the presence of complex structured spectral densities due to vibrational modes of the surroundings. We present numerical simulations on bacterial reaction centre (bRC) inspired model systems that utilize the recently developed MACGIC-iQUAPI method [Richter et al., J. Chem. Phys., 2017, 146, 214101]. The simulations demonstrate that the method provides a rigorous framework for the investigation of such integrated EET-CT dynamics. First, the applicability of the MACGIC-iQUAPI method is explored for a transition from monotonically decaying to oscillatory system-bath influence coefficients, a behavior inherently imposed by structured bath spectral densities. Tightly coupled EET and CT dynamics is further addressed for an excitonic subsystem that resembles strong coupling of special pair states and serves as donor towards a generic bridge-acceptor system. By solving the dissipative quantum dynamics of such bRC inspired model systems, the quenching of excitonic coherence on the hundreds of femtoseconds timescale is explored via a variation of the bridge state energetics, resembling a continuous transition from sequential to superexchange mediated CT regimes. Further, the simulations explore the influence of resonant vibrational modes on the quenching of excitonic coherence via CT. The results reveal a moderate influence of vibrational mode on charge separation dynamics in regimes of biologically relevant EET and CT dynamics.

Phosphate-Magnesium Ion Interactions in Water Probed by Ultrafast Two-Dimensional Infrared Spectroscopy.

Electric interactions between ions and ionic molecular groups in aqueous solution play a fundamental role in chemistry and biology. While Mg2+ ions are known to strongly affect the structure and folding dynamics of biomolecules, the relevance of different solvation geometries and the underlying interactions are mainly unresolved. We study dynamics and couplings between the hydrated Mg2+ and the dimethylphosphate anion, an established model system for the DNA and RNA backbone. The asymmetric (PO2-) stretching vibration serves as a sensitive noninvasive probe of phosphate-ion interactions. Femtosecond two-dimensional infrared (2D-IR) spectroscopy directly maps Mg2+ ions in contact with the phosphate groups via a distinct blue-shifted signature in the 2D spectrum. Data for different Mg2+ concentrations are analyzed by microscopic density functional theory modeling of cluster geometries and associated spectroscopic features, providing spatial assignments of the observed 2D-IR signatures. Phosphate-ion interactions arising from electrostatic Coulomb forces and exchange repulsion are the predominant origin of the observed frequency shifts.

Decay Pathways of Thymine Revisited.

The decay of electronically excited states of thymine (Thy) and thymidine 5'-monophosphate (TMP) was studied by time-resolved UV/vis and IR spectroscopy. In addition to the well-established ultrafast internal conversion to the ground state, a so far unidentified UV-induced species is observed. In D2O, this species decays with a time constant of 300 ps for thymine and of 1 ns for TMP. The species coexists with the lowest triplet state and is formed with a comparably high quantum yield of about 10% independent of the solvent. The experimentally determined spectral signatures are discussed in the light of quantum chemical calculations of the singlet and triplet excited states of thymine.

Vibrational Dynamics and Couplings of the Hydrated RNA Backbone: A Two-Dimensional Infrared Study.

The equilibrium structure of the RNA sugar-phosphate backbone and its hydration shell is distinctly different from hydrated DNA. Applying femtosecond two-dimensional infrared (2D-IR) spectroscopy in a range from 950 to 1300 cm-1, we elucidate the character, dynamics, and couplings of backbone modes of a double-stranded RNA A-helix geometry in its aqueous environment. The 2D-IR spectra display a greater number of backbone modes than for DNA, with distinctly different lineshapes of diagonal peaks. Phosphate-ribose interactions and local hydration structures are reflected in the complex coupling pattern of RNA modes. Interactions with the fluctuating water shell give rise to spectral diffusion on a 300 fs time scale, leading to a quasi-homogeneous line shape of the symmetric (PO2)- stretching mode of the strongly hydrated phosphate groups. The RNA results are benchmarked by 2D-IR spectra of DNA oligomers in water and analyzed by molecular dynamics and quantum mechanical molecular mechanics simulations.