On the microscopic theory of polar solvation dynamics

A microscopic model of the time-resolved Stokes shift is developed. The model calculates the solvation dynamics by combining the atomic resolution of the solute structure with dipolar dynamics from the polarization structure factors of the homogeneous solvent. Calculations are made for coumarin 153 and quinoxaline optical dyes with atomic geometries and charge distributions taken from quantum calculations. Stokes shift dynamics is calculated and compared to experiment in high-temperature acetonitrile and methanol and in low-temperature 2-methyl-tetrahydrofurane using dielectric relaxation data from experiment.     

Excited state proton transfer and solvent relaxation of a 3-hydroxyflavone probe in lipid bilayers

The photophysics of a ratiometric fluorescent probe, N-[[4'- N, N-diethylamino-3-hydroxy-6-flavonyl]methyl]- N-methyl- N-(3-sulfopropyl)-1-dodecanaminium, inner salt (F2N12S), incorporated into phospholipid unilamellar vesicles is presented. The reconstructed time-resolved emission spectra (TRES) unravels a unique feature in the photophysics of this probe. TRES exhibit signatures of both an excited-state intramolecular proton transfer (ESIPT) and a dynamic Stokes shift associated with solvent relaxation in the lipid bilayer. The ESIPT is fast, being characterized by a risetime of approximately 30-40 ps that provides an equilibrium to be established between the excited normal (N*) and the ESIPT tautomer (T*) on a time scale of 100 ps. On the other hand, the solvent relaxation displays a bimodal decay kinetics with an average relaxation time of approximately 1 ns. The observed slow solvent relaxation dynamics likely embodies a response of nonspecific dipolar solvation coupled with formation of probe-water H-bonds as well as the relocation of the fluorophore in the lipid bilayer. Taking into account that ESIPT and solvent relaxation are governed by different physicochemical properties of the probe microenvironment, the present study provides a physical background for the multiparametric sensing of lipid bilayers using ESIPT based probes.     

Solvation in supercritical water

Solvation in supercritical water under equilibrium and nonequilibrium conditions is studied via molecular dynamics simulations. The influence of solute charge distributions and solvent density on the solvation structures and dynamics is examined with a diatomic probe solute molecule. It is found that the solvation structure varies dramatically with the solute dipole moment, especially in low-density water, in accord with many previous studies on ion solvation. This electrostrictive effect has important consequences for solvation dynamics. In the case of a nonequilibrium solvent relaxation, if there are sufficiently many water molecules close to the solute at the outset of the relaxation, the solvent response measured as a dynamic Stokes shift is almost completely governed by inertial rotations of these water molecules. By contrast, in the opposite case of a low local solvent density near the solute, not only rotations but also translations of water molecules play an important role in solvent relaxation dynamics. The applicability of a linear response is found to be significantly restricted at low water densities.     

Solvent structural relaxation dynamics in dipolar solvation studied by resonant pump polarizability response spectroscopy

Resonant pump polarizability response spectroscopy (RP-PORS) was used to study the isotropic and anisotropic solvent structural relaxation in solvation. RP-PORS is the optical heterodyne detected transient grating (OHD-TG) spectroscopy with an additional resonant pump pulse. A resonant pump excites the solute-solvent system and the subsequent relaxation of the solute-solvent system is monitored by the OHD-TG spectroscopy. This experimental method allows measuring the dispersive and absorptive parts of the signal as well as fully controlling the beam polarizations of incident pulses and signal. The experimental details of RP-PORS were described. By performing RP-PORS with Coumarin 153(C153) in CH(3)CN and CHCl(3), we have successfully measured the isotropic and anisotropic solvation polarizability spectra following electronic excitation of C153. The isotropic solvation polarizability responses result from the isotropic solvent structural relaxation of the solvent around the solute whereas the anisotropic solvation polarizability responses come from the anisotropic translational relaxation and orientational relaxation. The solvation polarizability responses were found to be solvent-specific. The intramolecular vibrations of CHCl(3) were also found to be coupled to the electronic excitation of C153.     

Solvent response and dielectric relaxation in supercooled butyronitrile

We have measured the dynamics of solvation of a triplet state probe, quinoxaline, in the glass-forming dipolar liquid butyronitrile near its glass transition temperature T(g)=95 K. The Stokes shift correlation function displays a relaxation time dispersion of considerable magnitude and the optical linewidth changes along the solvation coordinate in a nonmonotonic fashion. These features are characteristic of solvation in viscous solvents and clearly indicate heterogeneous dynamics, i.e., spatially distinct solvent response times. Using the dielectric relaxation data of viscous butyronitrile as input, a microscopic model of dipolar solvation captures the relaxation time, the relaxation dispersion, and the amplitude of the dynamical Stokes shift remarkably well.     

Bandwidth analysis of solvation dynamics in a simple liquid mixture

The time-dependent energy distribution of solvation dynamics is studied by molecular dynamics simulations of a Lennard-Jones mixture. We calculate the response functions of the average and the variance which correspond to the spectral peak shift and bandwidth. Our calculation shows that the variance relaxation is slower than that of the average. The result agrees qualitatively with the experimental results. Dividing the obtained response functions into subcomponents caused by each solvent, we find that the relaxation is dominated by that solvent which strongly interacts with the solute. Extracting the redistribution component from the response functions, we find that it causes the slower relaxation of the response function. Thus, we conclude that the difference of the slower relaxations between the average and variance is caused by the redistribution process.     

Viscosity dependent radiationless relaxation rate of cyanine dyes. A picosecond laser spectroscopy study

The viscosity dependent radiationless relaxation of several cyanine dyes has been studied by picosecond laser spectroscopy. It was found that the relaxation rate is proportional to η^?α. The value of α, however, is not constant for a certain dye molecule, but is strongly dependent on the kind of solvent used. In n-alcohols for instance α is typically about 1. In glycerol/methanol or glycerol/water mixtures on the other hand α ≈ 0.5. A comparison is made with literature data on orientational relaxation lifetimes of some dyes in similar solvents. It is shown that the radiationless relaxation of cyanine dyes and the orientational relaxation of for instance xanthene dyes changes in roughly the same way as the solvent is changed. This is taken as proof of the proposal that a torsional motion of the heterocyclic quinolyl rings is the main course of the viscosity dependent relaxation of the cyanine dyes studied.     

Experimental and theoretical investigations of solvation dynamics of ionic fluids: appropriateness of dielectric theory and the role of DC conductivity

An analysis is provided of the subnanosecond dynamic solvation of ionic liquids in particular and ionic solutions in general. It is our hypothesis that solvation relaxation in ionic fluids, in the nonglassy and nonsupercooled regimes, can be understood rather simply in terms of the dielectric spectra of the solvent. This idea is suggested by the comparison of imidazolium ionic liquids with their pure organic counterpart, butylimidazole (J. Phys. Chem. B 2004, 108, 10245-10255). It is borne out by a calculation of the solvation correlation time from frequency dependent dielectric data for the ionic liquid, ethylammonium nitrate, and for the electrolyte solution of methanol and sodium perchlorate. Very good agreement is obtained between these theoretically calculated solvation relaxation functions and those obtained from fluorescence upconversion spectroscopy. Our comparisons suggest that translational motion of ions may not be the predominant factor in short-time solvation of ionic fluids and that many tools and ideas about solvation dynamics in polar solvents can be adapted to ionic fluids.     

Two-dimensional measurements of the solvent structural relaxation dynamics in dipolar solvation

Resonant-pump polarizability response spectroscopy (RP-PORS) is based on an optical heterodyne detected transient grating (OHD-TG) method with an additional resonant pump pulse. In RP-PORS, the resonant pump pulse excites the solute-solvent system and the subsequent relaxation of the solute-solvent system is monitored by the OHD-TG spectroscopy. RP-PORS is shown to be an excellent experimental tool to directly measure the solvent responses in solvation. In the present work, we extended our previous RP-PORS (Park et al., Phys. Chem. Chem. Phys., 2011, 13, 214-223) to measure time-dependent transient solvation polarizability (TSP) spectra with Coumarin153 (C153) in acetonitrile. The time-dependent TSP spectra showed how the different solvent intermolecular modes were involved in different stages of the solvation process. Most importantly, the inertial and diffusive components of the solvent intermolecular modes in solvation were found to be spectrally and temporally well-separated. In a dipolar solvation of C153, high-frequency inertial solvent modes were found to be driven instantaneously and decay on a subpicosecond timescale while low-frequency diffusive solvent modes were induced slowly and decayed on a picosecond timescale. Our present result is the first experimental manifestation of frequency-dependent solvent intermolecular response in a dipolar solvation.     

Solvatochromism effects on the dipole moments and photo-physical behavior of some azo sulfonamide dyes

Absorption and emission spectra of three azo sulfonamide dyes with various molecular structures have been studied in different solvents. The solute photo-physical behavior depends strongly on the solvent-solute interactions and solvent microenvironment. In order to understand the effect of intermolecular interactions on spectral behaviors of these dyes in different solvents and to conceive nature and extent of solvent-solute interactions the spectral variations were analyzed by the linear solvation energy relationships concept. In addition, by means of solvatochromic method the dipole moments of these dyes, in ground and excited states, were investigated.     

Dynamics of solvation and rotational relaxation of Coumarin 153 in ionic liquid confined nanometer-sized microemulsions

The effects of confinement of the ionic liquid 1-butyl-3-methylimidazolium tetrafluoroborate on solvation dynamics and rotational relaxation of Coumarin 153 (C-153) in Triton X-100/cyclohexane microemulsions have been explored using steady-state and picosecond time-resolved emission spectroscopy. The steady-state and rotational relaxation data indicate that C-153 molecules are incorporated in the core of the microemulsions. The average rotational relaxation time increases with increase in w ([bmim][BF(4)]/[TX-100]) values. The solvent relaxation in the core of the microemulsion occurs on two different time scales and is almost insensitive to the increase in w values. The solvent relaxation is retarded in the pool of the microemulsions compared to the neat solvent. Though, the retardation is very small compared to several-fold retardation of the solvation time of the conventional solvent inside the pool of the microemulsions.     

Solvation dynamics by coherence period resolved transient grating

We report a third-order nonlinear time-domain method, coherence period (τ) resolved transient grating (TRTG), that gives accurate solvation dynamics free from population relaxation in a short data acquisition time. The validity of TRTG is established by theory and experiment. The TRTG signal is shown to follow the transition frequency correlation function by an analytic expression based on the response function theory for delta function pulses and by model numerical calculations including finite pulse durations. TRTG is demonstrated for two cyanine dyes IR144 and IR125 in methanol by using a diffractive-optics based four wave mixing apparatus. Solvation dynamics in methanol obtained from the TRTG are consistent with those reported previously confirming the validity of TRTG.     

Preferential solvation and solvation shell composition of free base and protonated 5, 10, 15, 20-tetrakis(4-sulfonatophenyl)porphyrin in aqueous organic mixed solvents

The solvatochromic properties of the free base and the protonated 5, 10, 15, 20-tetrakis(4-sulfonatophenyl)porphyrin (TPPS) were studied in pure water, methanol, ethanol (protic solvents), dimethylsulfoxide, DMSO, (non-protic solvent), and their corresponding aqueous-organic binary mixed solvents. The correlation of the empirical solvent polarity scale (E(T)) values of TPPS with composition of the solvents was analyzed by the solvent exchange model of Bosch and Roses to clarify the preferential solvation of the probe dyes in the binary mixed solvents. The solvation shell composition and the synergistic effects in preferential solvation of the solute dyes were investigated in terms of both solvent-solvent and solute-solvent interactions and also, the local mole fraction of each solvent composition was calculated in cybotactic region of the probe. The effective mole fraction variation may provide significant physico-chemical insights in the microscopic and molecular level of interactions between TPPS species and the solvent components and therefore, can be used to interpret the solvent effect on kinetics and thermodynamics of TPPS. The obtained results from the preferential solvation and solvent-solvent interactions have been successfully applied to explain the variation of equilibrium behavior of protonation of TPPS occurring in aqueous organic mixed solvents of methanol, ethanol and DMSO.     

Solvation dynamics of the electron produced by two-photon ionization of liquid polyols. 1. Ethylene glycol

Solvated electrons have been produced in ethylene glycol by two-photon ionization of the solvent with 263 nm femtosecond laser pulses. The two-photon absorption coefficient of ethylene glycol at 263 nm is determined to be beta = (2.1 +/- 0.2) x 10(-11) m W(-1). The dynamics of electron solvation in ethylene glycol has been studied by pump-probe transient absorption spectroscopy. So, time-resolved absorption spectra ranging from 430 to 710 nm have been measured. A blue shift of the spectra is observed for the first tens of picoseconds. Using the Bayesian data analysis method, the observed solvation dynamics are reconstructed with different models: stepwise mechanisms, continuous relaxation models, or combinations of stepwise and continuous relaxation. Comparison between models is in favor of continuous relaxation, which is mainly governed by solvent molecular motions.     

A molecular dynamics computer simulation study of room-temperature ionic liquids. II. Equilibrium and nonequilibrium solvation dynamics

The molecular dynamics (MD) simulation study of solvation structure and free energetics in 1-ethyl-3-methylimidazolium chloride and 1-ethyl-3-methylimidazolium hexafluorophosphate using a probe solute in the preceding article [Y. Shim, M. Y. Choi and H. J. Kim, J. Chem. Phys. 122, 044510 (2005)] is extended to investigate dynamic properties of these liquids. Solvent fluctuation dynamics near equilibrium are studied via MD and associated time-dependent friction is analyzed via the generalized Langevin equation. Nonequilibrium solvent relaxation following an instantaneous change in the solute charge distribution and accompanying solvent structure reorganization are also investigated. Both equilibrium and nonequilibrium solvation dynamics are characterized by at least two vastly different time scales--a subpicosecond inertial regime followed by a slow diffusive regime. Solvent regions contributing to the subpicosecond nonequilibrium relaxation are found to vary significantly with initial solvation configurations, especially near the solute. If the solvent density near the solute is sufficiently high at the outset of the relaxation, subpicosecond dynamics are mainly governed by the motions of a few ions close to the solute. By contrast, in the case of a low local density, solvent ions located not only close to but also relatively far from the solute participate in the subpicosecond relaxation. Despite this difference, linear response holds reasonably well in both ionic liquids.     

Solvatochromism,halochromism, and preferential solvation of new dipolar guaiazulenyl 1,4-benzoquinone methides

The synthesis and the solvatochromic properties of five dyes, obtained by condensation of guaiazulene with 4-hydroxybenzaldehydes, are described. Crystal structures of a quinoid dye and a phenolic dye precursor are presented. The dyes are sensitive to the dipolarity-polarizability of the medium and to the hydrogen-bond donor ability of protic solvents. Their solvatochromism is discussed in terms of Kamlet-Taft's pi* and alpha scales, and their difference in behaviour is interpreted. Alkali and alkaline earth metal salts effect halochromism, with one exception due to extreme steric hindrance. Thus, this dye is capable of measuring solvent polarities without sensing the presence of electrolytes. Preferential solvation of the dyes in a series of binary solvent mixtures is explained quantitatively by solvent-exchange models.     

Dynamics of a supercooled ionic liquid studied by optical and dielectric spectroscopy

We have measured the dynamics of solvation of a triplet state probe, quinoxaline, in the glass-forming ionic liquid dibutylammonium formate near its glass transition temperature Tg=153 K. The Stokes-shift correlation function displays a relaxation time dispersion of considerable magnitude and the optical line width changes systematically along the solvation coordinate. The solvent dynamics in the viscous regime is compared with the rotational behavior of the solute and with the dielectric relaxation of the ionic liquid. Among the different quantities derived from the dielectric experiments, the relaxation of the macroscopic electric field, i.e., the modulus Mt, matches best the solvent response Ct regarding time scale and stretching exponent. Many other properties are reminiscent of the behavior of polar molecular liquids which lack the ionic character.     

Anisotropic reorientation of 9-methylpurine and 7-methylpurine molecules in methanol solution studied by combining 13C and 14N nuclear spin relaxation data and quantum chemical calculations

Reorientation of 9-(trideuteromethyl)purine and 7-(trideuteromethyl)purine molecules in methanol-d4 solutions has been investigated on the basis of the interpretation of the nuclear spin relaxation rates of their 14N (or 1H) and 13C nuclei. The transverse quadrupole relaxation rates of 14N nuclei have been obtained from the line shape analysis of their 14N NMR spectra. Alternatively, the information on the longitudinal 14N relaxation rates has been obtained via the scalar relaxation of the second kind of protons coupled to 14N. The longitudinal dipolar relaxation rates of the protonated 13C nuclei in the investigated molecules have been determined by measuring their overall relaxation rates and NOE enhancement factors. The molecular geometries, scalar coupling constants, and EFG tensors needed for quantitative interpretation of the above data have been calculated theoretically [DFT B3LYP/6-311++G(2d,p) or B3PW91/6-311+G(df,pd)] including the impact of the solvent by using discrete solvation and the polarizable continuum model. The reorientation of the investigated purines has been described as rotational diffusion of an asymmetrical top. It has been found that to get a fully consistent interpretation of the relaxation data, effective C-H bond lengths being 3% longer than the calculated ones had to be used in analysis to compensate for the ground-state vibrations. The obtained rotational diffusion coefficients and orientations of the principal diffusion axes show that the investigated molecules reorient anisotropically and that the mode of their solvation is remarkably different, in spite of their structural similarity.     

Anchor points for the unified Brønsted acidity scale: the rCCC model for the calculation of standard Gibbs energies of proton solvation in eleven representative liquid media

The COSMO cluster-continuum (CCC) solvation model is introduced for the calculation of standard Gibbs solvation energies of protons. The solvation sphere of the proton is divided into an inner proton-solvent cluster with covalent interactions and an outer solvation sphere that interacts electrostatically with the cluster. Thus, the solvation of the proton is divided into two steps that are calculated separately: 1) The interaction of the proton with one or more solvent molecules is calculated in the gas phase with high-level quantum-chemical methods (modified G3 method). 2) The Gibbs solvation energy of the proton-solvent cluster is calculated by using the conductor-like screening model (COSMO). For every solvent, the solvation of the proton in at least two (and up to 11) proton-solvent clusters was calculated. The resulting Gibbs solvation energies of the proton were weighted by using Boltzmann statistics. The model was evaluated for the calculation of Gibbs solvation energies by using experimental data of water, MeCN, and DMSO as a reference. Allowing structural relaxation of the proton-solvent clusters and the use of structurally relaxed Gibbs solvation energies improved the accordance with experimental data especially for larger clusters. This variation is denoted as the relaxed COSMO cluster-continuum (rCCC) model, for which we estimate a 1σ error bar of 10 kJ mol(-1) . Gibbs solvation energies of protons in the following representative solvents were calculated: Water, acetonitrile, sulfur dioxide, dimethyl sulfoxide, benzene, diethyl ether, methylene chloride, 1,2-dichloroethane, sulfuric acid, fluorosulfonic acid, and hydrogen fluoride. The obtained values are absolute chemical standard potentials of the proton (pH=0 in this solvent). They are used to anchor the individual solvent specific acidity (pH) scales to our recently introduced absolute acidity scale.