Solvation Correlation Function Research Articles

Solute dynamics of a hydrophobic molecule in a menthol-thymol based type-V deep eutectic solvent: effect of composition of the components.

Type-V deep eutectic solvents (DESs) are a newly emerging unique class of solvents obtained by physical mixing and heating of non-ionic components. These solvents show deviation from the thermodynamic ideality. Compared to type-I to IV DESs, type-V DESs are less explored and their physical chemistry is in its nascent stage. In this work, we have chosen a type-V DES based on menthol-thymol (MT) for our working media. Solvent and rotational dynamics were studied with varying temperature using a well-known solvatochromic probe, Coumarin 153 (C153). We prepared the MT-based DES using a reported procedure at three molar ratios: 1 : 1 (M1T1), 1 : 1.5 (M1T1.5), and 2 : 1 (M2T1) of menthol (M) and thymol (T). Time-resolved emission spectra (TRES) were constructed with varying temperature. Utilizing TRES, the decay of the solvent correlation function (C(t)) was plotted. We have correlated the solvent relaxation time in these DESs as a function of viscosity. The time-resolved anisotropy decays were also collected to perceive the rotational relaxation dynamics of C153 as a function of temperature. The decay of solvent relaxation was found to be bi-exponential, and the average solvation time (〈τs〉) in M2T1 was found to be longer than those of M1T1.5 and M1T1. The rotational reorientation times (〈τrot〉) also follow the same trend. We have analysed the rotational dynamics of C153 in type-V DESs employing the Stokes-Einstein-Debye (SED) hydrodynamic model. The rotational dynamics in DESs demonstrate a good correlation with the SED model with a little deviation. In MT-based DESs, the solute's rotational relaxation times approach hydrodynamic stick boundary condition at low viscosity (or at high temperatures) for all molar compositions. Using the Arrhenius-type equations, we have correlated the activation energies for the rotational motion of C153, along with the viscous flow and non-radiative pathways for all the DESs.

Sigma enlarging bridge correction of three dimensional Ornstein–Zernike theory for solvation free energy of polyatomic solutes immersed in Lennard-Jones monatomic solvent

We have studied the applicability of the sigma enlarging bridge (SEB) correction method to the solvation of polyatomic solute molecules in a Lennard-Jones monatomic solvent using the three dimensional Ornstein–Zernike (3D-OZ) theory. It is found that the SEB correction improves the solvation free energy (SFE) significantly. It has been concluded from the analysis of the radial distribution function (RDF) that the parameter included in the SEB function can be transferred from the values of the monatomic solute to those of the polyatomic one for the 3D-OZ theory. The one dimensional reference interaction site model (1D-RISM) theory has also been examined. The SEB correction is found to be applicable to improve the SFE of the 1D-RISM theory. Except for buried or sterically hindered atoms of the solute, the transferability of the SEB parameter has been confirmed in terms of the RDF obtained by the 1D-RISM theory. This paper also examines the applicability of the hybrid closure between the molecular dynamics simulation and the one of the following closure equations—the hyper-netted chain, Kovalenko–Hirata, or Kobryn–Gusarov–Kovalenko equation—for preparing the solvent–solvent correlation function. Using the results of the hybrid closure in addition to the SEB-corrected closure, we discuss the effect of the quality of the correlation functions for the bulk solvent on the accuracy of the SFE.

Extending correlation functions of molecular dynamics simulation by Kovalenko–Hirata and Kobryn–Gusarov–Kovalenko closures for monatomic Lennard-Jones solvent and its application to a calculation of solvation

We have tried to extend the short-ranged radial distribution function of the molecular dynamics (MD) simulation for a monatomic Lennard-Jones solvent by applying the Ornstein–Zernike theory with two hybrid closures. One was the hybrid with the Kovalenko–Hirata or the KH closure (hybrid MD–KH closure), whereas the other was the hybrid with the Kobryn–Gusarov–Kovalenko or the KGK closure (hybrid MD–KGK closure). As long as the switching distance between the MD and the KH (or the KGK) is chosen appropriately, the direct correlation functions obtained by these hybrid closures were almost identical to each other, which also agree with those of the hybrid closure with the hypernetted chain or the HNC (hybrid MD–HNC closure). The calculations on the solute–solvent correlation function and the solvation free energy also supported the practical equivalence among the hybrid MD–KH, hybrid MD–KGK, and hybrid MD–HNC closures.

Ultrafast Excited-State Dynamics of Tricarbocyanine Dyes Probed by Two-Dimensional Electronic Spectroscopy: Polar Solvation vs Photoisomerization.

Photophysical properties of tricarbocyanine dyes in various solvents have been widely investigated using a variety of spectroscopic tools. However, the presence of several ground-state isomers and interconversion between these isomers on an ultrafast timescale upon photoexcitation render unambiguous assignment of spectral features quite difficult. In this work, ultrafast excited-state dynamics of two tricarbocyanine dyes in two solvents, DNTTCI and IR140, in ethanol and ethylene glycol, are studied by two-dimensional electronic spectroscopy (2DES). We present a detailed discussion on design and calibration of the 2DES apparatus and on the method for data processing by phase-cycling. For DNTTCI we report a method to obtain solvation correlation function, the nature of which is found to be strongly dependent on the excitation frequencies; a blue-shifted spectrum at early time is observed and explained based on preferential emission from a subset among various isomers having overlapping spectral features. For IR140 in ethanol, four isomers with distinct spectral features are identified, and most importantly, three of these isomers were found to interconvert upon photoexcitation which completes within 100 fs and is explained based on a kinetic model of consecutive chemical reaction. Density functional theory calculations show the presence of several ground-state isomers for both these dyes. Through this work we demonstrate how 2DES can help us to decipher distinct excited-state photophysics in two carbocyanine dyes, polar solvation and photoisomerization, by resolving spectral congestion without sacrificing time resolution.

Spectrum Correction in Study of Solvation Dynamics by Fluorescence Non-collinear Optical Parametric Amplification Spectroscopy

Femtosecond time-resolved fluorescence non-collinear optical parametric amplification spectroscopy can extract the curve of spectral gain from its parametric superfluorescence. This unique spectrum correction method enables fluorescence non-collinear optical parametric amplification spectroscopy acquiring the genuine transient fluorescence spectrum of the studied system. In this work we employ fluorescence non-collinear optical parametric amplification spectroscopy technique to study the solvation dynamics of DCM dye in ethanol solution, and confirm that genuine solvation correlation function and shift of peak frequency can be derived from transient fluorescence spectra after the spectral gain correction. It demonstrates that fluorescence non-collinear optical parametric amplification spectroscopy can benefit the research fields, which focuses on both fluorescence intensity dynamics and fluorescence spectral shape evolution.

Measuring correlated electronic and vibrational spectral dynamics using line shapes in two-dimensional electronic-vibrational spectroscopy.

Two-dimensional electronic-vibrational (2DEV) spectroscopy is an experimental technique that shows great promise in its ability to provide detailed information concerning the interactions between the electronic and vibrational degrees of freedom in molecular systems. The physical quantities 2DEV is particularly suited for measuring have not yet been fully determined, nor how these effects manifest in the spectra. In this work, we investigate the use of the center line slope of a peak in a 2DEV spectrum as a measure of both the dynamic and static correlations between the electronic and vibrational states of a dye molecule in solution. We show how this center line slope is directly related to the solvation correlation function for the vibrational degrees of freedom. We also demonstrate how the strength with which the vibration on the electronic excited state couples to its bath can be extracted from a set of 2DEV spectra. These analytical techniques are then applied to experimental data from the laser dye 3,3'-diethylthiatricarbocyanine iodide in deuterated chloroform, where we determine the lifetime of the correlation between the electronic transition frequency and the transition frequency for the backbone C = C stretch mode to be ∼1.7 ps. Furthermore, we find that on the electronic excited state, this mode couples to the bath ∼1.5 times more strongly than on the electronic ground state.

How fast can a proton-transfer reaction be beyond the solvent-control limit?

In this article, we review the field of photoacids. The rate of excited-state proton transfer (ESPT) to solvent spans a wide range of time scales, from tens of nanoseconds for the weakest photoacids to short time scales of about 100 fs for the strongest photoacids synthesized so far. We divide the photoacid strength into four regimes. Regime I includes the weak photoacids 0 < pKa* < 3. These photoacids can transfer a proton only to water or directly to a mild-base molecule in solution. The ESPT rate to other protic solvents, like methanol or ethanol, is too small in comparison with the radiative rate. The second regime includes stronger photoacids whose pKa*'s range from -4 to 0. They are capable of transferring a proton to other protic solvents and not only to water. The third regime includes even stronger photoacids. Their pKa* is ∼ -6, and the ESPT rate constant, kPT, is limited by the orientational time of the solvent which is characterized by the average solvation correlation function ⟨S(t)⟩. The fourth regime sets a new time limit for the ESPT rate of the strongest photoacids synthesized so far. The kPT value of such photoacids is 10(13) s(-1), and τPT = 100 fs. We attribute this new time limit (beyond the solvent control) to intermolecular vibration between the two heavy atoms of the proton donor and the proton acceptor, which assist the ESPT by lowering the height and width of the potential barrier, thus enhancing the ESPT rate.

Interfacial bond-breaking electron transfer in mixed water-ethylene glycol solutions: reorganization energy and interplay between different solvent modes.

We explore solvent dynamics effects in interfacial bond breaking electron transfer in terms of a multimode approach and make an attempt to interpret challenging recent experimental results (the nonmonotonous behavior of the rate constant of electroreduction of S2O82– from mixed water–EG solutions when increasing the EG fraction; see Zagrebin, P.A. et al. J. Phys. Chem. B2010, 114, 311). The exact expansion of the solvent correlation function (calculated using experimental dielectric spectra) in a series predicts the splitting of solvent coordinate in three independent modes characterized by different relaxation times. This makes it possible to construct a 5D free-energy surface along three solvent coordinates and one intramolecular degree of freedom describing first electron transfer at the reduction of a peroxodisulphate anion. Classical molecular dynamics simulations were performed to study the solvation of a peroxodisulphate anion (S2O82–) in oxidized and reduced states in pure water and ethylene glycol (EG) as well as mixed H2O–EG solutions. The solvent reorganization energy of the first electron-transfer step at the reduction of S2O82– was calculated for several compositions of the mixed solution. This quantity was found to be significantly asymmetric. (The reorganization energies of reduction and oxidation differ from each other.) The averaged reorganization energy slightly increases with increasing the EG content in solution. This finding clearly indicates that for the reaction under study the static solvent effect no longer competes with solvent dynamics. Brownian dynamics simulations were performed to calculate the electron-transfer rate constants as a function of the solvent composition. The results of the simulations explain the experimental data, at least qualitatively.

Nanoconfinement of Water Layers in Lamellar Structures Prepared in the Presence and Absence of Organic Solvent

An attempt is made to draw a line of comparison between the extent of rigidity of the hydration layers bound to the interfacial region of lamellar structures of Aerosol OT (AOT, sodium bis(2-ethylhexyl) sulfosuccinate) in water, in the presence and absence of an organic solvent using POM, SAXS, cryo-TEM, and time-resolved fluorescence spectroscopy. These systems are ternary mixtures of AOT, water, and n-heptane containing lamellar structures in an aqueous layer at higher w(0) values (w(0) = 300 and 150) and a binary solution of 20 and 50% AOT in neat water (w/w). The solvation shells residing at the vicinity of these lamellar structures are monitored using two different coumarin probes (C153 and C500). It is intended to envisage a comparative solvation dynamics study of the restricted aqueous region confined in lamellar structures formed in ternary mixture and binary solution. Though steady state measurements show a similar microenvironment probed by the fluorophores in lamellar structures formed in the two different aqueous phases, temporal evolution of the solvent correlation function C(t) unveils the existence of lamellar structures with different degrees of confinement of water layers in these two systems. A slower relaxation of the restricted aqueous region in lamellar structures of binary solution signifies the presence of more rigid interfacially bound water layers at the lamellar interface than in the ternary mixture having a similar weight percentage of AOT in water. The present investigation concludes that the lamellar structures formed under two different conditions provide a similar hydrophobic environment with different extents of localized water populations at the lamellar interface as manifested by the solvent relaxation time in agreement with SAXS and cryo-TEM images.

Ultrafast Solvation Dynamics of Subtilisin-Polyethylene Glycol Interaction for Protein Crystallization

We study the ultrafast solvation dynamics of protein-precipitant complexes. Protein subtilisin carlsberg (SC) was mixed with several polyethylene glycol (PEG) precipitants for protein crystallization. Picosecond-resolved emission spectra from single intrinsic tryptophan residue (Trp-113) are recorded to construct solvation correlation functions. For precipitant concentrations with various crystallization effects, we observe drastically different solvation relaxation processes. These differences in solvation dynamics are correlated with the local protein structural integrity and water-network stability upon interaction with the precipitants. The solvation dynamics at the protein surface is proposed as a new perspective to study precipitant-protein interactions.

Excited-State Intramolecular Hydrogen Atom Transfer of Curcumin in Surfactant Micelles

Femtosecond fluorescence upconversion experiments were performed on the naturally occurring medicinal pigment, curcumin, in anionic, cationic, and neutral micelles. In our studies, the micelles are composed of sodium dodecyl sulfate (SDS), dodecyl trimethyl ammonium bromide (DTAB), and triton X-100 (TX-100). We demonstrate that the excited-state kinetics of curcumin in micelles have a fast (3-8 ps) and slow (50-80 ps) component. While deuteration of curcumin has a negligible effect on the fast component, the slow component exhibits a pronounced isotope effect of approximately 1.6, indicating that micelle-captured curcumin undergoes excited-state intramolecular hydrogen atom transfer. Studies of solvation dynamics of curcumin in a 10 ps time window reveal a fast component (< or = 300 fs) followed by a 8, 6, and 3 ps component in the solvation correlation function for the TX-100, DTAB, and SDS micelles, respectively.

Transport properties of water at functionalized molecular interfaces

Understanding transport properties of solvent such as diffusion and viscosity at interfaces with biomacromolecules and hard materials is of fundamental importance to both biology and biotechnology. Our study utilizes equilibrium molecular dynamics simulations to calculate solvent transport properties at a model peptide and microarray surface. Both diffusion and selected components of viscosity are considered. Solvent diffusion is found to be affected near the peptide and surface. The stress-stress correlation function of solvent near the hard surface exhibits long time memory. Both diffusion and viscosity are shown to be closely correlated with the density distribution function of water along the microarray surface.

Considerations for the Construction of the Solvation Correlation Function and Implications for the Interpretation of Dielectric Relaxation in Proteins

The dielectric response of proteins is conveniently measured by monitoring the time-dependent Stokes shift of an associated chromophore. The interpretation of these experiments depends critically upon the construction of the solvation correlation function, C(t), which describes the time-dependence of the Stokes shift and hence the dielectric response of the medium to a change in charge distribution. We provide an analysis of various methods of constructing this function and review selected examples from the literature. The naturally occurring amino acid, tryptophan, has been frequently used as a probe of solvation dynamics in proteins. Its nonexponential fluorescence decay has stimulated the generation of an alternative method of constructing C(t). In order to evaluate this method, we have studied a system mimicking tryptophan. The system is comprised of two coumarins (C153 and C152) having different fluorescence lifetimes but similar solvation times. The coumarins are combined in different proportions in methanol to make binary probe mixtures. We use fluorescence upconversion spectroscopy to obtain wavelength-resolved kinetics of the individual coumarins in methanol as well as the binary mixtures of 75:25, 50:50, and 25:75 of C153:C152. The solvation correlation functions are constructed for these systems using different methods and are compared.

Absorption and emission lineshapes and ultrafast solvation dynamics of NO in parahydrogen

We perform a theoretical study on the electronic spectroscopy of dilute NO impurity embedded in parahydrogen (p-H(2)). Absorption and emission lineshapes for the A (2)Sigma(+)<--X (2)Pi Rydberg transition of NO in parahydrogen have been previously measured and simulated, which yielded results for the NO/p-H(2) ground and excited state pair potentials [L. Bonacina et al., J. Chem. Phys. 125, 054507 (2006)]. Using these potentials, we performed molecular dynamics simulation, theoretical statistical mechanical calculations of absorption and emission lineshapes, and both equilibrium and nonequilibrium solvation correlation functions for NO chromophore in parahydrogen. Theory was shown to be in good agreement with simulation. Linear response treatment of solvation dynamics was shown to break down due to a dramatic change in the solute-solvent microstructure upon solute excitation to the Rydberg state and the concomitant increase of the solute size.

Absorption and emission lineshapes and solvation dynamics of NO in supercritical Ar

We perform a theoretical study of electronic spectroscopy of dilute NO in supercritical Ar fluid. Absorption and emission lineshapes for the A(2)Sigma(+)<--X(2)Pi Rydberg transition of NO in argon have been previously measured and simulated, which yielded results for the NO/Ar ground- and excited-state pair potentials [Larregaray et al., Chem. Phys. 308, 13 (2005)]. Using these potentials, we have performed molecular dynamics simulations and theoretical statistical mechanical calculations of absorption and emission lineshapes and nonequilibrium solvation correlation functions for a wide range of solvent densities and temperatures. Theory was shown to be in good agreement with simulation. Linear response treatment of solvation dynamics was shown to break down at near-critical temperature due to dramatic change in the solute-solvent microstructure upon solute excitation to the Rydberg state and the concomitant increase of the solute size.

Dynamic Solvation in Phosphonium Ionic Liquids: Comparison of Bulk and Micellar Systems and Considerations for the Construction of the Solvation Correlation Function, C(t)

Dynamic solvation of the dye coumarin 153 is studied in a phosphonium ionic liquid: hexadecyltributylphosphonium bromide, [(C4)3C16P+][Br-]. It forms micelles in water, and the bulk also exists as a liquid under our experimental conditions. This system permits a comparison with an imidazolium ionic liquid studied earlier, which also formed micelles in water (J. Phys. Chem. A 2006, 110, 10725-10730). We conclude that our analysis of the comparable situation in a phosphonium liquid is not as definitive as we had proposed earlier, i.e., that the majority of the early-time solvation arises from the organic cation. Part of the difficulty in performing this analysis is most likely due to the amount of water that is associated with the micelle. In the course of this work, we have focused on the calculation of the solvation correlation function, C(t), and investigated how it depends upon the methods with which the "zero-time" spectrum is constructed.

Ultrafast Proton Transfer to Solvent: Molecularity and Intermediates from Solvation- and Diffusion-Controlled Regimes

Photoinduced proton transfer (PT) from cations 6-hydroxyquinolinium (6HQc) and 6-hydroxy-1-methylquinolinium (6MQc) to water and alcohols, and solvation of the zwitterionic conjugate base 1-methylquinolinium-6-olate (6MQz) were studied with stationary and transient absorption spectroscopy and by quantum chemical calculations. Transient emission spectra from 6MQz in acetonitrile and protic solvents shift dynamically to the red without changing their shape and intensity. The shift matches the solvation correlation function C(t) either measured with known solvatochromic probes coumarin 343 and coumarin 153 or derived from infrared/dielectric-loss data on neat solvents. This indicates that 6MQz monitors the solvation dynamics and that no intramolecular electron transfer occurs on a subpicosecond or longer time scale. The PT dynamics S(t) from 6HQc and 6MQc closely follows C(t), being initially 2-3 times slower. This allows for the conclusion that PT is controlled by solvation, with a barrier of 2 kJ/mol. In water, a pre-condition of this ultrafast reaction seems to be hydrogen-bonding between the negatively charged oxygen and two water molecules, resulting in a complex 6HQc:H2O:H2O. The complex is stable due to a high (47 kJ/mol) bonding energy between 6HQc and a water molecule. In acetonitrile, the reaction equilibrium is strongly shifted to the cation. There an intermediate PT state was detected, which may be ascribed to the cationic form 6HQc:H2O due to residual water impurities. In water-acetonitrile mixtures, the ultrafast solvent-controlled PT is followed by a diffusion-controlled reaction; the measured rate kD approximately 1010 s-1 M-1 is characteristic for simple bimolecular diffusion. The dependence of the short-time PT signal on water concentration can be fitted with a Poisson distribution of water molecules around the cation. Altogether, the short-time and long-time behaviors provide strong evidence that diffusion of only one water molecule is sufficient to detach the proton. Subsequent solvent stabilization of the products completes the PT reaction.

Ultrafast Solvation Dynamics of Human Serum Albumin: Correlations with Conformational Transitions and Site-Selected Recognition

Human serum albumin, the most abundant protein found in blood plasma, transports a great variety of ligands in the circulatory system and undergoes reversible conformational transitions over a wide range of pH values. We report here our systematic studies of solvation dynamics and local rigidity in these conformations using a single intrinsic tryptophan (W214) residue as a local molecular probe. With femtosecond resolution, we observed a robust bimodal distribution of time scales for all conformational isomers. The initial solvation occurs in several picoseconds, representing the local librational/rotational motions, followed by the dynamics, in the tens to hundreds of picoseconds, which result from the more bonded water in the tryptophan crevice. Under the physiological condition of neutral pH, we measured approximately 100 ps for the decay of the solvation correlation function and observed a large wobbling motion at the binding site that is deeply buried in a crevice, revealing the softness of the binding pocket and the large plasticity of the native structure. At acidic pH, the albumin molecule transforms to an extended conformation with a large charge distribution at the surface, and a similar temporal behavior was observed. However, at the basic pH, the protein opens the crevice and tightens its globular structure, and we observed significantly faster dynamics, 25-45 ps. These changes in the solvation dynamics are correlated with the conformational transitions and related to their structural integrity.

Generalized Langevin theory on the dynamics of simple fluids under external fields

A theory on the time development of the density and current fields of simple fluids under an external field is formulated through the generalized Langevin formalism. The theory is applied to the linear solvation dynamics of a fixed solute regarding the solute as the external field on the solvent. The solute-solvent-solvent three-body correlation function is taken into account through the hypernetted-chain integral equation theory, and the time correlation function of the random force is approximated by that in the absence of the solute. The theoretical results are compared with those of molecular-dynamics (MD) simulation and the surrogate theory. As for the transient response of the density field, our theory is shown to be free from the artifact of the surrogate theory that the solvent can penetrate into the repulsive core of the solute during the relaxation. We have also found a large quantitative improvement of the solvation correlation function compared with the surrogate theory. In particular, the short-time part of the solvation correlation function is in almost perfect agreement with that from the MD simulation, reflecting that the short-time expansion of the theoretical solvation correlation function is exact up to t(2) with the exact three-body correlation function. A quantitative improvement is found in the long-time region, too. Our theory is also applied to the force-force time correlation function of a fixed solute, and similar improvement is obtained, which suggests that our present theory can be a basis to improve the mode-coupling theory on the solute diffusion.

Ultrafast dynamics in a nanocage of enzymes: Solvation and fluorescence resonance energy transfer in reverse micelles

In this contribution we report studies of the nature of solvation and resonance energy transfer processes in a reverse micelle (RM) upon encapsulation of a digestive enzyme, α-chymotrypsin (CHT). We have used one donor, Coumarin 500 (C500), and three acceptors Rhodamine 123 (R123, cationic), ethidium bromide (EtBr, cationic), and Merocyanine 540 (MC540, anionic). By selectively exciting the donor at the surface of the RM with a proper excitation wavelength we have examined solvation dynamics in the microenvironment. The solvation correlation function in the RM without CHT exhibits single-exponential decay with time constant ∼660 ps, which is similar to that of the CHT-included RM. However, in the case of CHT-included RM ( w 0 = 10 ) , the time-resolved anisotropy and spectral linewidth analysis of the surface-bound donor reveal the existence of an annular aqueous channel of thickness ∼2.5 Å between the enzyme surface and the inner surface of the RM. The aqueous channel is a potential host for the water-soluble substrate and also is involved in maintaining the proper functionality of RM encapsulated CHT. The studies use both steady-state and time-resolved fluorescence resonance energy transfer (FRET) techniques to measure donor–acceptor distances in the RM and also emphasize the danger of using steady-state fluorescence quenching as a method in careful estimation of the distances. The local geometrical restriction on the donor and acceptor molecules was estimated from time-resolved polarization (anisotropy) measurements. The time-resolved anisotropy of the donor and acceptor molecules also revealed significant randomization of the relative orientation of transition dipoles of the donor and acceptor, justifying the use of 2 / 3 as the value of the orientation factor κ 2 . These studies attempt to elucidate the excellence of the RM as a nanohost of biological macromolecules.