Femtosecond solvation dynamics in acetonitrile: Observation of the inertial contribution to the solvent response

On the nature of initial solvation in bulk polar liquids: Gaussian or exponential?

Measurement of time evolution of fluorescence of a probe solute has been a quintessential technique to quantify how dipolar solvent molecules dynamically minimize the free energy of an electronically excited probe. During such solvation dynamics in bulk liquids, a substantial part of relaxation was shown to complete within sub-100 fs from time-gated fluorescence measurements, as also predicted by molecular dynamics simulation studies. However, equivalent quantification of solvation timescales by femtosecond pump-probe and broadband fluorescence measurements revealed an exponential nature of this initial relaxation having quite different timescales. Here, we set out to unveil the reason behind these puzzling contradictions. We introduce a method for estimating probe wavelength-dependent instrument response and demonstrate that the observation of the Gaussian vs exponential nature of initial relaxation is indeed dependent on the method of data analysis. These findings call for further experimental investigation and parallel development of theoretical models to elucidate the molecular-level mechanism accounting for different types of early time solvation.

Excited-State Symmetry Breaking in Solvent Mixtures

A large number of multipolar dyes undergo excited-state symmetry breaking (ESSB) in polar media. During this process, electronic excitation, initially distributed evenly over the molecule, localizes, at least partially, on one donor-acceptor branch. To resolve its initial stage, ESSB is investigated with a donor-acceptor-donor dye in binary mixtures of nonpolar and polar solvents using time-resolved infrared absorption spectroscopy. The presence of a few polar molecules around the dye is sufficient to initiate ESSB. Although the extent of asymmetry in a mixture is close to that in a pure solvent of similar polarity, the dynamics are slower and involve translational diffusion. However, preferential solvation in the mixtures leads to a larger local polarity. Furthermore, inhomogeneous broadening of the S1 ← S0 absorption band of the dye is observed in the mixtures, allowing for a photoselection of solutes with different local environments and ESSB dynamics.

Photodynamics of azaindoles in polar media: the influence of the environment

We have studied the relaxation dynamics of a family of azaindole (AI) structural isomers, 4-, 5-, 6- and 7-AI, by steady-state and time-resolved methods (fs-transient absorption and fluorescence up-conversion), in solvents of different polarity. The measurements in aprotic solvents show distinctive fluorescence yields and excited state lifetimes among the isomers, which are tuned by the polarity of the medium. Guided by simple TD-DFT calculations and based on the behavior observed in the isolated species, it has been possible to address the influence of the environment polarity on the relaxation route. According to the obtained picture, the energy of the nπ* state, which is strongly dependent on the position of the pyridinic nitrogen, controls the rate of the internal conversion channel that accounts for the distinctive photophysical behavior of the isomers. On the other hand, preliminary measurements in protic media (methanol) show a very different photodynamical behavior, in which the anomalous measured fluorescent patterns are very likely the result of reactive channels (proton transfer) triggered by the electronic excitation.

Observation of parallel intersystem crossing and charge transfer-state dynamics in [Fe(bpy)3]2+ from ultrafast 2D electronic spectroscopy

Transition metal-based charge-transfer complexes represent a broad class of inorganic compounds with diverse photochemical applications. Charge-transfer complexes based on earth-abundant elements have been of increasing interest, particularly the canonical [Fe(bpy)3]2+. Photoexcitation into the singlet metal-ligand charge transfer (1MLCT) state is followed by relaxation first to the ligand-field manifold and then to the ground state. While these dynamics have been well-studied, processes within the MLCT manifold that facilitate and/or compete with relaxation have been more elusive. We applied ultrafast two-dimensional electronic spectroscopy (2DES) to disentangle the dynamics immediately following MLCT excitation of this compound. First, dynamics ascribed to relaxation out of the initially formed 1MLCT state was found to correlate with the inertial response time of the solvent. Second, the additional dimension of the 2D spectra revealed a peak consistent with a ∼20 fs 1MLCT → 3MLCT intersystem crossing process. These two observations indicate that the complex simultaneously undergoes intersystem crossing and direct conversion to ligand-field state(s). Resolution of these parallel pathways in this prototypical earth-abundant complex highlights the ability of 2DES to deconvolve the otherwise obscured excited-state dynamics of charge-transfer complexes.

Understanding Charge Dynamics in Dense Electronic Manifolds in Complex Environments

Photoinduced charge transfer (CT) excited states and their relaxation mechanisms can be highly interdependent on the environment effects and the consequent changes in the electronic density. Providing a molecular interpretation of the ultrafast (subpicosecond) interplay between initial photoexcited states in such dense electronic manifolds in condensed phase is crucial for improving and understanding such phenomena. Real-time time-dependent density functional theory is here the method of choice to observe the charge density, explicitly propagated in an ultrafast time domain, along with all time-dependent properties that can be easily extracted from it. A designed protocol of analysis for real-time electronic dynamics to be applied to time evolving electronic density related properties to characterize both in time and in space CT dynamics in complex systems is here introduced and validated, proposing easy to be read cross-correlation maps. As case studies to test such tools, we present the photoinduced charge-transfer electronic dynamics of 5-benzyluracil, a mimic of nucleic acid/protein interactions, and the metal-to-ligand charge-transfer electronic dynamics in water solution of [Ru(dcbpy)2(NCS)2]4-, dcbpy = (4,4'-dicarboxy-2,2'-bipyridine), or "N34-", a dye sensitizer for solar cells. Electrostatic and explicit ab initio treatment of solvent molecules have been compared in the latter case, revealing the importance of the accurate modeling of mutual solute-solvent polarization on CT kinetics. We observed that explicit quantum mechanical treatment of solvent slowed down the charge carriers mobilities with respect to the gas-phase. When all water molecules were modeled instead as simpler embedded point charges, the electronic dynamics appeared enhanced, with a reduced hole-electron distance and higher mean velocities due to the close fixed charges and an artificially increased polarization effect. Such analysis tools and the presented case studies can help to unveil the influence of the electronic manifold, as well as of the finite temperature-induced structural distortions and the environment on the ultrafast charge motions.

Ultrafast Excited State Dynamics of Forward and Reverse trans-cis Photoisomerization of Red-Light-Absorbing Indigo Derivatives

Femtosecond time-resolved transient absorption (TRTA) spectroscopy was carried out to investigate the ultrafast excited state dynamics of both trans → cis and trans ← cis photoisomerization of red-light-absorbing indigo derivatives. For N,N'-bis(tert-butyloxycarbonylmethyl)indigo (tBOMI), the excited state lifetime of the trans-form was measured to be 41 ps while that of the cis-form was as short as 730 fs in acetonitrile (Acn). The excited state lifetime of trans-N,N'-dimethylindigo (DMI) in Acn was also measured to be as short as 10 ps. These values are much shorter than those of the blue-light-absorbing trans-forms of indigo derivatives such as N,N'-diacetylindigo (DAI) and thioindigo (ThI). The chromophore of indigo is composed of two pairs of electron donor and acceptor substituents conjugated in the shape of a letter "H" (so-called "H-chromophore"), although DFT and TDDFT calculations suggest that the charge transfer (CT) character is not very significant. Nevertheless, when a weak CT within the H-chromophore is promoted, the absorption band shifts to longer wavelengths and the excited state lifetime shortens. For the photoisomerization of DAI and ThI, a relatively long excited state lifetime is required for the photoisomerization, while for tBOMI and DMI, a vibrationally hot ground state that overcomes the energy barrier in the ground state is produced by rapid nonradiative decay through conical intersection. In the case of cis-tBOMI, the repulsion between the two adjacent negatively charged carbonyl groups and the weakening of the central C═C double bond in the S₁ state twist the molecule, shorten the excited state lifetime, and increase the quantum yield of the trans ← cis photoisomerization.

In Silico Ultrafast Nonlinear Spectroscopy Meets Experiments: The Case of Perylene Bisimide Dye

Spectroscopy simulations are of paramount importance for the interpretation of experimental electronic spectra, the disentangling of overlapping spectral features, and the tracing of the microscopic origin of the observed signals. Linear and nonlinear simulations are based on the results drawn from electronic structure calculations that provide the necessary parameterization of the molecular systems probed by light. Here, we investigate the applicability of excited-state properties obtained from linear-response time-dependent density functional theory (TDDFT) in the description of nonlinear spectra by employing the pseudowavefunction approach and compare them with benchmarks from highly accurate RASSCF/RASPT2 calculations and with high temporal resolution experimental results. As a test case, we consider the prediction of femtosecond transient absorption and two-dimensional electronic spectroscopy of a perylene bisimide dye in solution. We find that experimental signals are well reproduced by both theoretical approaches, showing that the computationally cheaper TDDFT can be a suitable option for the simulation of nonlinear spectroscopy of molecular systems that are too large to be treated with higher-level RASSCF/RASPT2 methods.

Ultrafast Fluctuations in PM6 Domains of Binary and Ternary Organic Photovoltaic Thin Films Probed with Two-Dimensional White-Light Spectroscopy

We present two-dimensional white-light spectroscopy (2DWL) measurements of binary and ternary bulk heterojunctions of the polymer donor PM6 mixed with state-of-the-art nonfullerene acceptors Y6 or IT4F. The ternary film has a shorter lifetime and faster spectral diffusion than either of the binary films. 2D line shape analysis of the PM6 ground state bleach with a Kubo model determines that all three films have similar amplitudes of fluctuations (Δ = 0.29 fs^-1) in their transition frequencies, but different relaxation times (ranging from 102 to 24 fs). The ternary film exhibits faster dynamics than either of the binary films. The short lifetime of the ternary blend is consistent with increased photoexcitation transfer and the fast frequency fluctuations are consistent with structural dynamics of aliphatic side chains. These results suggest that the femtosecond fluctuations of PM6 are impacted by the choice of the acceptor molecules. We hypothesize that those dynamics are either indicative, or perhaps the initial source, of structural dynamics that ultimately contribute to solar cell operation.

Electron-Nuclear Dynamics Accompanying Proton-Coupled Electron Transfer

Although photoinduced proton-coupled electron transfer (PCET) plays an essential role in photosynthesis, a full understanding of the mechanism is still lacking due to the complex nonequilibrium dynamics arising from the strongly coupled electronic and nuclear degrees of freedom. Here we report the photoinduced PCET dynamics of a biomimetic model system investigated by means of transient IR and two-dimensional electronic-vibrational (2DEV) spectroscopies, IR spectroelectrochemistry (IRSEC), and calculations utilizing long-range-corrected hybrid density functionals. This collective experimental and theoretical effort provides a nuanced picture of the complicated dynamics and synergistic motions involved in photoinduced PCET. In particular, the evolution of the 2DEV line shape, which is highly sensitive to the mixing of vibronic states, is interpreted by accurate computational modeling of the charge separated state and is shown to represent a gradual change in electron density distribution associated with a dihedral twist that occurs on a 120 fs time scale.

Solvation dynamics in electronically polarizable solvents: Theoretical treatment using solvent-polarizable three-dimensional reference interaction-site model theory combined with time-dependent density functional theory

The theory of solvation structure in an electronically polarizable solvent recently proposed by us, referred to as the "solvent-polarizable three-dimensional reference interaction-site model theory," is extended to dynamics in this study through the combination with time-dependent density functional theory. Test calculations are performed on model charge-transfer systems in water, and the effects of electronic polarizability on solvation dynamics are examined. The electronic polarizability slightly retards the solvation dynamics. This is ascribed to the decrease in the curvature of the nonequilibrium free energy profile along the solvation coordinate. The solvent relaxation is bimodal, and the faster and the slower modes are assigned to the reorientational and the translational modes, respectively, as was already reported by the surrogate theory combined with the site-site Smoluchowski-Vlasov equation. The relaxation path along the solvation coordinate is a little higher than the minimum free energy path because the translational mode is fixed in the time scale of the reorientational relaxation.

Lifetime Broadening and Impulsive Generation of Vibrational Coherence Triggered by Ultrafast Electron Transfer

The absorption band shape of chromophores in liquid solution at room temperature is usually dominated by pure electronic dephasing dynamics, which occurs on the sub-100 fs time scale. Herein, we report on a series of dyads consisting of a naphthalenediimide (NDI) electron acceptor with one or two phenyl-based donors for which photoinduced intramolecular electron transfer is fast enough to be competitive with pure electronic dephasing. As a consequence, the absorption band of the π-π* transition of these dyads is broader than that of the NDI alone to an extent that scales with the electron transfer rate. Additionally, this reaction is so fast that it leads to the impulsive excitation of a low-frequency vibrational mode of the charge-separated product. Quantum-chemical calculations suggest that this vibration involves the C-N donor-acceptor bond, which shortens considerably upon electron transfer.

Proton Transfer Dynamics in the Aprotic Proton Accepting Solvent 1-Methylimidazole

The dynamics of proton transfer to the aprotic solvent 1-methylimidazole (MeIm, proton acceptor) from the photoacid 8-hydroxypyrene-1,3,6-trisulfonic acid trisodium salt (HPTS) was investigated using fast fluorescence measurements. The closely related molecule, 8-methoxypyrene-1,3,6-trisulfonic acid trisodium salt (MPTS), which is not a photoacid, was also studied for comparison. Following optical excitation, the wavelength-dependent population dynamics of HPTS in MeIm resulting from the deprotonation process were collected over the entire fluorescence emission window. Analysis of the time-dependent fluorescence spectra revealed four distinct fluorescence bands that appear and decay on different time scales. We label these four states as protonated (P), associated I (A_I), associated II (A_II), and deprotonated (D). We find that the simple kinetic scheme of P → A_I → A_II → D is not consistent with the data. Instead, the kinetic scheme that describes the data has P decaying into A_I, which mainly goes on to deprotonation (D), but A_I can also feed into A_II. A_II can return to A_I or decay to the ground state, but does not deprotonate within experimental error. Quantum chemistry and excited state QM/MM Born-Oppenheimer molecular dynamics simulations indicate that A_I and A_II are two H-bonding conformations of MeIm to the HPTS hydroxyl, axial, and equatorial, respectively.

Solvation Dynamics in Water. 4. On the Initial Regime of Solvation Relaxation

It is shown, by means of numerical and analytic work, that initial molecular momenta play little significant role in the initial fast solvation relaxation that follows electronic excitation of, and charge creation for, a standard model system of a solute in water. Instead, the nonequilibrium dynamics are predominantly described by noninertial "steering" by the torques directly generated by the newly created charge distribution. It is this process that largely overcomes inertia and drives the relaxation dynamics on a time scale of a few tens of femtoseconds in the key initial regime of the dynamics. These results are discussed in the context of commonly employed descriptions such as inertial, Gaussian, and underdamped dynamical behavior.

Quantum yield and energy efficiency of photoinduced intramolecular charge separation

Kinetics of photoinduced intramolecular charge separation (CS) and the ensuing ultrafast charge recombination (CR) in electron-donor-acceptor dyads are studied numerically, taking into account the excitation of charge-transfer active intramolecular vibrations and multiple relaxation time scales of the surrounding polar solvent. Both energetic and dynamic properties of intramolecular and solvent reorganization are considered, and their influence on the CS/CR kinetics and quantum yield of ultrafast CS is explored. Particular attention is paid to the energy efficiency of CS, as one of the most important parameters indicating the promise of using a molecular compound as a basis for emerging optoelectronic devices. The CS quantum yield and the energy efficiency of CS are shown to depend differently on the key model parameters. Necessary conditions for the highly efficient CS are evaluated using analytic formulae for the electron transfer rates and derived from numerical simulation data. The reasons why low-exergonic CS taking place in the Marcus normal region can be much slower than CR in the deep inverted region are discussed.

Time-Resolved Probing of the Nonequilibrium Structural Solvation Dynamics by the Time-Dependent Stokes Shift

The time-dependent fluorescence Stokes shift monitors the relaxation of the polarization of a polar solvent in the surroundings of a photoexcited solute molecule but also the structural variation of the solute following photoexcitation and the subsequent molecular charge redistribution. Here, we formulate a simple nonequilibrium quantum theory of solvation for an explicitly time-dependent continuous solvent. The time-dependent solvent induces nonequilibrium fluctuations on the solvent dynamics which are directly reflected in different time components in the time-dependent Stokes shift. We illustrate the structural dynamics in the presence of an explicitly time-dependent solvent by the example of a dynamically shrinking solute which leads to a bimodal Stokes shift. Interestingly, both contributions are mutually coupled. Furthermore, we can explain a prominent long-tail decay of the Stokes shift associated with slow structural dynamical variations.

Photophysics of perylene monoimide-labelled organocatalysts

A fluorophore-tagged organocatalyst undergoes electron transfer in polar solvents allowing to sense the presence of its free quinuclidine catalytic site.

Changes in protein hydration dynamics by encapsulation or crowding of ubiquitin: strong correlation between time-dependent Stokes shift and intermolecular nuclear Overhauser effect

The local changes in protein hydration dynamics upon encapsulation of the protein or macromolecular crowding are essential to understand protein function in cellular environments. We were able to obtain a spatially-resolved picture of the influence of confinement and crowding on the hydration dynamics of the protein ubiquitin by analyzing the time-dependent Stokes shift (TDSS), as well as the intermolecular Nuclear Overhauser Effect (NOE) at different sites of the protein by large-scale computer simulation of single and multiple proteins in water and confined in reverse micelles. Besides high advanced space resolved information on hydration dynamics we found a strong correlation of the change in NOE upon crowding or encapsulation and the change in the integral TDSS relaxation times in all investigated systems relative to the signals in a diluted protein solution.

Changes in local protein hydration dynamics caused by encapsulation or crowding are reflected in the TDSS and the intermolecular NOE alike.

Spectroscopic and Calorimetric Studies of Molecular Recognitions in a Dendrimer-Surfactant Complex

Molecular recognitions, causing supramolecular complex formation between a hyperbranched polymer molecule (polyamidoamine (PAMAM) dendrimer generation 3) with oppositely charged surfactant sodium dodecyl sulfate (SDS) in aqueous solution, were studied by using various spectroscopic techniques and calorimetric titration of heat change measurements. Spectroscopic measurements were performed using dynamic Stokes shift (DSS), rotational anisotropy decay, and translational diffusion of a fluorescent probe molecule coumarin 153 (C153) noncovalently attached to the dendrimer-surfactant complex. All these studies unanimously confirm that the critical aggregation concentration (CAC) of SDS falls to ∼0.8 mM (from its critical micelle concentration (CMC) ∼ 8 mM) in the presence of ∼0.2 mM dendrimer. Further studies of isothermal titration calorimetry (ITC) measurement show that the CAC of SDS in the presence of dendrimer remains invariant to the dendrimer concentration. Complexation reaction between SDS and dendrimer is highly exothermic in nature. A maximum heat release (ΔH∼ -6.6 kJ/mol of SDS binding) was observed at a SDS-to-dendrimer mole ratio of ∼3-5; where up to 3 to 5 SDS molecules were encapsulated by one dendrimer molecule to form dendrimer-SDS encapsulation complex. When negatively charged SDS was replaced with a positively charged surfactant dodecyl-trimethylammonium-bromide (DTAB), we found that the DTAB hardly interacted with positively charged dendrimer due to the charge-charge repulsions.

The effect of solvent relaxation time constants on free energy gap law for ultrafast charge recombination following photoinduced charge separation

Ultrafast low exergonic charge recombination following photoinduced charge separation proceeds in a non-equilibrium mode and its rate constant is nearly independent of the free energy gap.

Ultrafast Solvation Dynamics and Vibrational Coherences of Halogenated Boron-Dipyrromethene Derivatives Revealed through Two-Dimensional Electronic Spectroscopy

Boron-dipyrromethene (BODIPY) chromophores have a wide range of applications, spanning areas from biological imaging to solar energy conversion. Understanding the ultrafast dynamics of electronically excited BODIPY chromophores could lead to further advances in these areas. In this work, we characterize and compare the ultrafast dynamics of halogenated BODIPY chromophores through applying two-dimensional electronic spectroscopy (2DES). Through our studies, we demonstrate a new data analysis procedure for extracting the dynamic Stokes shift from 2DES spectra revealing an ultrafast solvent relaxation. In addition, we extract the frequency of the vibrational modes that are strongly coupled to the electronic excitation, and compare the results of structurally different BODIPY chromophores. We interpret our results with the aid of DFT calculations, finding that structural modifications lead to changes in the frequency, identity, and magnitude of Franck-Condon active vibrational modes. We attribute these changes to differences in the electron density of the electronic states of the structurally different BODIPY chromophores.

Role of Dispersive Fluorous Interaction in the Solvation Dynamics of the Perfluoro Group Containing Molecules

Perfluoro group containing molecules possess an important self-aggregation property through the fluorous (F···F) interaction which makes them useful for diverse applications such as medicinal chemistry, separation techniques, polymer technology, and biology. In this article, we have investigated the solvation dynamics of coumarin-153 (C153) and coumarin-6H (C6H) in ethanol (ETH), 2-fluoroethanol (MFE), and 2,2,2-trifluoroethanol (TFE) using the femtosecond upconversion technique and molecular dynamics (MD) simulation to understand the role of fluorous interaction between the solute and solvent molecules in the solvation dynamics of perfluoro group containing molecules. The femtosecond upconversion data show that the time scales of solvation dynamics of C6H in ETH, MFE, and TFE are approximately the same whereas the solvation dynamics of C153 in TFE is slow as compared to that of ETH and MFE. It has also been observed that the time scale of solvation dynamics of C6H in ETH and MFE is higher than that of C153 in the same solvents. MD simulation results show a qualitative agreement with the experimental data in terms of the time scale of the slow components of the solvation for all the systems. The experimental and simulation studies combined lead to the conclusion that the solvation dynamics of C6H in all solvents as well as C153 in ETH and MFE is mostly governed by the charge distribution of ester moieties (C═O and O) of dye molecules whereas the solvation of C153 in TFE is predominantly due to the dispersive fluorous interaction (F···F) between the perfluoro groups of the C153 and solvent molecules.

Time-resolved signatures across the intramolecular response in substituted cyanine dyes

The optically populated excited state wave packet propagates along multidimensional intramolecular coordinates soon after photoexcitation. This action occurs alongside an intermolecular response from the surrounding solvent. Disentangling the multidimensional convoluted signal enables the possibility to separate and understand the initial intramolecular relaxation pathways over the excited state potential energy surface. Here we track the initial excited state dynamics by measuring the fluorescence yield from the first excited state as a function of time delay between two color femtosecond pulses for several cyanine dyes having different substituents. We find that when the high frequency pulse precedes the low frequency one and for timescales up to 200 fs, the excited state population can be depleted through stimulated emission with efficiency that is dependent on the molecular electronic structure. A similar observation at even shorter times was made by scanning the chirp (frequencies ordering) of a femtosecond pulse. The changes in depletion reflect the rate at which the nuclear coordinates of the excited state leave the Franck-Condon (FC) region and progress towards achieving equilibrium. Through functional group substitution, we explore these dynamic changes as a function of dipolar change following photoexcitation. Density functional theory calculations were performed to provide greater insight into the experimental spectroscopic observations. Complete active space (CAS) self-consistent field and CAS second order perturbation theory calculated potential energy surfaces tracking twisting and pyramidalization confirm that the steeper potential at the FC region leads to the observation of faster wave packet dynamics.

Probe dependence on polar solvation dynamics from fs broadband fluorescence

Solvation dynamics is remarkably independent of the probe as long as specific interactions remain similar.

Broad-Band Pump-Probe Spectroscopy Quantifies Ultrafast Solvation Dynamics of Proteins and Molecules

In this work, we demonstrate the use of broad-band pump-probe spectroscopy to measure femtosecond solvation dynamics. We report studies of a rhodamine dye in methanol and cryptophyte algae light-harvesting proteins in aqueous suspension. Broad-band impulsive excitation generates a vibrational wavepacket that oscillates on the excited-state potential energy surface, destructively interfering with itself at the minimum of the surface. This destructive interference gives rise to a node at a certain probe wavelength that varies with time. This reveals the Gibbs free-energy changes of the excited-state potential energy surface, which equates to the solvation time correlation function. This method captures the inertial solvent response of water (∼40 fs) and the bimodal inertial response of methanol (∼40 and ∼150 fs) and reveals how protein-buried chromophores are sensitive to the solvent dynamics inside and outside of the protein environment.