Simulations of time-dependent fluorescence in nano-confined solvents

Deprotonation of formic, acetic acids and bicarbonate ion in slit silica nanopores at infinite dilution and in the presence of electrolytes

Dielectric effects and the coupled electrostatics between the nanoconfined and the internal/external aqueous media contribute to the observed deviations of chemistry within the nanoconfined environment when compared with unconfined systems. A systematic understanding has remained elusive, especially with respect to background salt concentration and boundary condition effects like the nanopore surface chemistry and the reference state used to calculate free energies. We utilize molecular dynamics simulations along with thermodynamic integration to determine the free energy difference associated with acid-base chemistry in 2 nm and 4 nm slit pores open to a bulk-like reservoir. pKa increases are predicted when confining acetic acid, formic acid, and bicarbonate in the slits at infinite dilution conditions. We find that confinement weakens the acids, and the modulation of outer pore surface dipole magnitudes can tune the pKa shift values, suggesting that purely "intrinsic" electrostatic effect on confinement may not exist. At sufficiently high salt concentrations, the dielectric/electrostatic effects on pKa values diminish due to charge screening effects. These discoveries enable future modifications of nanopore chemistries to achieve desirable properties for industrial applications.

A Model Electron Transfer Reaction in Confined Aqueous Solution

Liquid water confined within nanometer-sized channels exhibits a strongly reduced local dielectric constant perpendicular to the wall, especially at the interface, and this has been suggested to induce faster electron transfer kinetics at the interface than in the bulk. We study a model electron transfer reaction in aqueous solution confined between graphene sheets with classical molecular dynamics. We show that the solvent reorganization energy is reduced at the interface compared to the bulk, which explains the larger rate constant. However, this facilitated solvent reorganization is due to the partial desolvation by the graphene sheet of the ions involved in the electron transfer and not to a local dielectric constant reduction effect.

Signatures of nanoconfinement on the linear and nonlinear vibrational spectroscopy of a model hydrogen-bonded complex dissolved in a polar solvent

The one-dimensional IR (1D-IR) absorption and IR pump-probe spectra of a hydrogen stretch in a model hydrogen-bonded complex dissolved in a polar solvent confined in spherical hydrophobic cavities of different sizes were simulated using ground-state mixed quantum-classical dynamics. Due to a thorough analysis of key properties of the complex and solvent from equilibrium trajectory data, we were able to gain insight into the microscopic details underlying the spectra. Both the 1D-IR and IR pump-probe spectra manifested the effects of confinement on the relative stabilities of the covalent and ionic forms of the complex through pronounced changes in their peak intensities and numbers. However, in contrast to the 1D-IR spectra, the time-resolved pump-probe spectra were found to be uniquely sensitive to the changes in the molecular dynamics as the cavity size is varied. In particular, it was found that the variations in the time evolutions of the peak intensities in the pump-probe spectra reflect the differences in the solvation dynamics associated with the various forms of the complex in different locations within the cavities. The ability to detect these differences underscores the advantage of using pump-probe spectroscopy for studying nanoconfined systems.

Solvation dynamics and proton transfer in nanoconfined liquids

Nanoconfined liquids are of interest because of both their fundamental properties and their potential utility in an array of applications. The structure and dynamics of the liquid can be dramatically impacted by the geometrical constraints and the interactions with the interface. Understanding the molecular-level origins of these changes and how they are determined by the characteristics of the confining framework is the subject of ongoing experimental and theoretical studies. The progress and remaining challenges in these efforts are reviewed in the context of solvation dynamics and proton transfer reactions, processes that are strongly affected by nanoscale confinement.

Molecular simulation study of water mobility in aerosol-OT reverse micelles

In this work, we present results from molecular dynamics simulations on the single-molecule relaxation of water within reverse micelles (RMs) of different sizes formed by the surfactant aerosol-OT (AOT, sodium bis(2-ethylhexyl)sulfosuccinate) in isooctane. Results are presented for RM water content w(0) = [H(2)O]/[AOT] in the range from 2.0 to 7.5. We show that translational diffusion of water within the RM can, to a good approximation, be decoupled from the translation of the RM through the isooctane solvent. Water translational mobility within the RM is restricted by the water pool dimensions, and thus, the water mean-squared displacements (MSDs) level off in time. Comparison with models of diffusion in confined geometries shows that a version of the Gaussian confinement model with a biexponential decay of correlations provides a good fit to the MSDs, while a model of free diffusion within a sphere agrees less well with simulation results. We find that the local diffusivity is considerably reduced in the interfacial region, especially as w(0) decreases. Molecular orientational relaxation is monitored by examining the behavior of OH and dipole vectors. For both vectors, orientational relaxation slows down close to the interface and as w(0) decreases. For the OH vector, reorientation is strongly affected by the presence of charged species at the RM interface and these effects are especially pronounced for water molecules hydrogen-bonded to surfactant sites that serve as hydrogen-bond acceptors. For the dipole vector, orientational relaxation near the interface slows down more than that for the OH vector due mainly to the influence of ion-dipole interactions with the sodium counterions. We investigate water OH and dipole reorientation mechanisms by studying the w(0) and interfacial shell dependence of orientational time correlations for different Legendre polynomial orders.

Model silica pores with controllable surface chemistry for molecular dynamics simulatinos

Model amorphous silica pores have been developed for use in molecular dynamics simulations. Specifically, roughly cylindrical pores have been constructed with hydrophilic, hydroxyl-terminated surfaces. The approach is designed to allow systematic variation of the pore radius and surface functionality. Thus, these pores are suitable for studying the variability in solvent structure, energy transfer and reaction dynamics occurring inside the pore due to surface modification. The method is described and the properties of the generated pores are discussed.

Water-mediated ordering of nanoparticles in an electric field

Interfacial polar molecules feature a strongly anisotropic response to applied electric field, favoring dipole orientations parallel to the interface. In water, in particular, this effect combines with generic orientational preferences induced by spatial asymmetry of water hydrogen bonding under confined geometry, which may give rise to a Janus interface. The two effects manifest themselves in considerable dependence of water polarization on both the field direction relative to the interface and the polarity (sign) of the field. Using molecular simulations, we demonstrate strong field-induced orientational forces acting on apolar surfaces through water mediation. At a field strength comparable to electric fields around a DNA polyion, the torques we predict to act on an adjacent nanoparticle are sufficient to overcome thermal fluctuations. These torques can align a particle with surface as small as 1 nm2. The mechanism can support electrically controlled ordering of suspended nanoparticles as a means of tuning their properties and can find application in electro-nanomechanical devices.

Interaction of ionic liquid with water with variation of water content in 1-butyl-3-methyl-imidazolium hexafluorophosphate ([bmim][PF6])/TX-100/water ternary microemulsions monitored by solvent and rotational relaxation of coumarin 153 and coumarin 490

The interaction of water with room temperature ionic liquid (RTIL) [bmim][PF6] has been studied in [bmim][PF6]/TX-100/water ternary microemulsions by solvent and rotational relaxation of coumarin 153 (C-153) and coumarin 490 (C-490). The rotational relaxation and average solvation time of C-153 and C-490 gradually decrease with increase in water content of the microemulsions. The gradual increase in the size of the microemulsion with increase in w0 (w0=[water]/[surfactant]) is evident from dynamic light scattering measurements. Consequently the mobility of the water molecules also increases. In comparison to pure water the retardation of solvation time in the RTIL containing ternary microemulsions is very less. The authors have also reported the solvation time of C-490 in neat [bmim][PF6]. The solvation time of C-490 in neat [bmim][PF6] is bimodal with time constants of 400 ps and 1.10 ns.

On the connection between Gaussian statistics and excited-state linear response for time-dependent fluorescence

Time-dependent fluorescence (TDF) of a chromophore in a polar or nonpolar solvent is frequently simulated using linear-response approximations. It is shown that one such linear-response-type approximation for the TDF Stokes shift derived by Carter and Hynes [J. Chem. Phys. 94, 5961 (1991)] that is based on excited-state dynamics gives the same result as that obtained by assuming Gaussian statistics for the energy gap. The derivation provides insight into the much discussed relationship between linear response and Gaussian statistics. In particular, subtle but important differences between the two approximations are illuminated that suggest that the result is likely more generally applicable than suggested by the usual linearization procedure. In addition, the assumption of Gaussian statistics directly points to straightforward checks of the validity of the approximation with essentially no additional computational effort.

Excitation Wavelength Dependence of Solvation Dynamics in a Supramolecular Assembly: PEO−PPO−PEO Triblock Copolymer and SDS

The triblock copolymer (PEO)20-(PPO)70-(PEO)20 (P123) forms a supramolecular aggregate with sodium dodecyl sulfate (SDS). The solvation dynamics and anisotropy decay of coumarin 480 (C480) in different regions of a P123-SDS aggregate are studied through variation of the excitation wavelength (lambdaex) using femtosecond upconversion. In a P123 micelle, because of the drastic differences in polarity between the hydrophilic corona region (PEO block) and the hydrophobic PPO core, C480 exhibits a pronounced red edge excitation shift (REES) of emission maximum by 24 nm. In the P123-SDS aggregate, SDS penetrates the core of the P123 micelle. This increases the polarity of the core and reduces the difference in the polarity between the core and the corona region. In a P123-SDS aggregate, the REES is much smaller (5 nm) which suggests a reduced difference between the core and the corona. Solvation dynamics in a P123 micelle displays a bulklike ultrafast component (<0.3 and 1 ps) in the PEO corona region, a 200 ps component arising from dynamics of polymer segments, and a very long component (5000 or 3000 ps) due to the highly restricted PPO core. In a P123-SDS aggregate, at lambdaex = 375 and 405 nm, the solvation dynamics is found to be faster than that in P123 micelle. In this case, the component (3000 ps) arising from the core region is faster than that (5000 ps) in P123 micelle. In both P123 micelle and P123-SDS aggregate, the relative contribution of the core region decreases and that of the corona region increases with an increase in lambdaex. At lambdaex = 435 nm, which probes the hydrophilic corona, the solvation dynamics for both P123 micelle and P123-SDS aggregate are almost similar.

Ultrafast Electron Transfer in a Nanocavity. Dimethylaniline to Coumarin Dyes in Hydroxypropyl γ-Cyclodextrin

Photoinduced electron transfer (PET) from N,N-dimethylaniline (DMA) to four coumarin dyes (C151, C481, C153, and C480) inside the cavity of hydroxypropyl gamma-cyclodextrin (hpCD) is studied using femtosecond upconversion. The rate of PET is found to be nonexponential and to vary significantly with the coumarin dyes. The rate for C481 is 100 times faster than that for C480. The PET rate displays a bell-shaped dependence on the free energy change and thus reveals a Marcus-type inverted region. The anisotropy decay of the four dyes in hpCD are found to be very similar, and hence the observed variation in the rate of PET is not due to variation in rotational diffusion of the acceptors (coumarin dyes).

Aqueous Solvation Dynamics at Metal Oxide Surfaces

Broadband transient absorption (TA) spectroscopy, three-pulse photon echo peak shift (3PEPS), and anisotropy decay measurements were used to study the solvation dynamics in bulk water and interfacial water at ZrO(2) surfaces, using Eosin Y as a probe. The 3PEPS results show a multiexponential behavior with two subpicosecond components that are similar in bulk and interfacial water, while a third component of several picoseconds is significantly lengthened at the interface. The bandwidth correlation function from TA spectra exhibits the same behavior, and the TA spectra are well reproduced using the doorway-window picture with the time constants from PEPS. Our results suggest that interfacial water is restricted to a thickness of less than 5 A. Also the high-frequency collective dynamics of water does not seem to be affected by the interface. On the other hand, the increase of the third component may point to a slowing down of diffusional motion at the interface, although other effects, may play a role, which are discussed.

Testing the Core/Shell Model of Nanoconfined Water in Reverse Micelles Using Linear and Nonlinear IR Spectroscopy

A core/shell model has often been used to describe water confined to the interior of reverse micelles. The validity of this model for water encapsulated in AOT/isooctane reverse micelles ranging in diameter from 1.7 to 28 nm (w0 = 2-60) and bulk water is investigated using four experimental observables: the hydroxyl stretch absorption spectra, vibrational population relaxation times, orientational relaxation rates, and spectral diffusion dynamics. The time dependent observables are measured with ultrafast infrared spectrally resolved pump-probe and vibrational echo spectroscopies. Major progressive changes appear in all observables as the system moves from bulk water to the smallest water nanopool, w0 = 2. The dynamics are readily distinguishable for reverse micelle sizes smaller than 7 nm in diameter (w0 = 20) compared to the response of bulk water. The results also demonstrate that the size dependent absorption spectra and population relaxation times can be quantitatively predicted using a core-shell model in which the properties of the core (interior of the nanopool) are taken to be those of bulk water and the properties of the shell (water associated with the headgroups) are taken to be those of w0 = 2. A weighted sum of the core and shell components reproduces the size dependent spectra and the nonexponential population relaxation dynamics. However, the same model does not reproduce the spectral diffusion and the orientational relaxation experiments. It is proposed that, when hydrogen bond structural rearrangement is involved (orientational relaxation and spectral diffusion), dynamical coupling between the shell and the core cause the water nanopool to display more homogeneous dynamics. Therefore, the absorption spectra and vibrational lifetime decays can discern different hydrogen bonding environments whereas orientational and spectral diffusion correlation functions predict that the dynamics are size dependent but not as strongly spatially dependent within a reverse micelle.

Temperature dependence of solvation dynamics and anisotropy decay in a protein: ANS in bovine serum albumin

Temperature dependence of solvation dynamics and fluorescence anisotropy decay of 8-anilino-1-naphthalenesulfonate (ANS) bound to a protein, bovine serum albumin (BSA), are studied. Solvation dynamics of ANS bound to BSA displays a component (300 ps) which is independent of temperature in the range of 278-318 K and a long component which decreases from 5800 ps at 278 K to 3600 ps at 318 K. The temperature independent part is ascribed to a dynamic exchange of bound to free water with a low barrier. The temperature variation of the long component of solvation dynamics corresponds to an activation energy of 2.1 kcal mol(-1). The activation energy is ascribed to local segmental motion of the protein along with the associated water molecules and polar residues. The time scale of solvation dynamics is found to be very different from the time scale of anisotropy decay. The anisotropy decays are analyzed in terms of the wobbling motion of the probe (ANS) and the overall tumbling of the protein.

Dynamics of Solvent and Rotational Relaxation of Glycerol in the Nanocavity of Reverse Micelles

The dynamics of solvent and rotational relaxation of Coumarin 480 and Coumarin 490 in glycerol containing bis-2-ethyl hexyl sulfosuccinate sodium salt (AOT) reverse micelles have been investigated with steady-state and time-resolved fluorescence spectroscopy. We observed slower solvent relaxation of glycerol confined in the nanocavity of AOT reverse micelles compared to that in pure glycerol. However, the slowing down in the solvation time on going from neat glycerol to glycerol confined reverse micelles is not comparable to that on going from pure water or acetonitrile to water or acetonitrile confined AOT reverse micellar aggregates. While solvent relaxation times were found to decrease with increasing glycerol content in the reverse micellar pool, rotational relaxation times were found to increase with increase in glycerol content.

Femtosecond Study of Partially Folded States of Cytochrome C by Solvation Dynamics

Using femtosecond time-resolved fluorescence spectroscopy, it is shown that the solvation dynamics in the two partially folded states (IS' and IS' ') of a protein, cytochrome C, are very different. In the case of IS' (formed by the addition of 2 mM sodium dodecyl sulfate, SDS) almost the entire dynamic solvent shift of coumarin 153 (C153) is captured in a picosecond setup and the contribution of the ultrafast component (0.5 ps) is very small (5%). Solvation dynamics of IS' ' (formed by 2 mM SDS and 5 M urea) displays a major component (47%) of 1.3 ps. This indicates that the structure of IS' ' is much more open and exposed compared to that of IS'. The difference in the dynamics of IS' and IS' ' is attributed to differences in their structure, particularly near the heme region, and the presence of urea in IS' '.

Fluorescence Anisotropy Decay and Solvation Dynamics in a Nanocavity: Coumarin 153 in Methyl β-Cyclodextrins

Fluorescence anisotropy decay and solvation dynamics of coumarin 153 (C153) are studied in dimethyl beta-cyclodextrin (DIMEB) and trimethyl beta-cyclodextrin (TRIMEB) nanocavity in water. C153 binds to DIMEB and TRIMEB to form both 1:1 and 1:2 (C153:cyclodextrin) complexes. The anisotropy decays of C153 in DIMEB and TRIMEB are found to be biexponential. The fast component of anisotropy decay (approximately 1000 ps) is attributed to the 1:1 complex and the slower one (approximately 2500 ps) to the 1:2 complex. From the components of the anisotropy decay, the length of the 1:1 and 1:2 complexes are estimated. Solvation dynamics of C153 in DIMEB exhibits a very fast (2.4 ps) component (41%) and two slower components of 50 ps (29%) and 1450 ps (30%). Solvation dynamics in TRIMEB is described by three slow components of 10.3 ps (24%), 240 ps (45%), and 2450 ps (31%). Possible origins of the ultraslow components are discussed.

Temperature Dependence of Anisotropy Decay and Solvation Dynamics of Coumarin 153 in γ-Cyclodextrin Aggregates

Effect of temperature on the fluorescence anisotropy decay and the ultraslow component of solvation dynamics of coumarin 153 (C153) in a gamma-cyclodextrin (gamma-CD) nanocavity are studied using a picosecond set up. The steady-state anisotropy (0.13 +/- 0.01) and residual anisotropy (0.14 +/- 0.01) in fluorescence anisotropy decay in an aqueous solution containing 7 microM C153 and 40 mM gamma-CD are found to be quite large. This indicates formation of large linear nanotube aggregates of gamma-CD linked by C153. It is estimated that >53 gamma-CD units are present in each aggregate. In these aggregates with rise in temperature, the average solvation time (<tau(s)>(obs)) decreases markedly from 680 ps at 278 K to 160 ps at 318 K. The dynamic Stokes shift is found to decrease from 800 cm(-1) at 278 K to 250 cm(-1) at 318 K. The fraction of dynamic Stokes shift (f(d)) detected in a picosecond set up is calculated using the Fee-Maroncelli procedure. The corrected solvation time (<tau(s)>(corr) = f(d)<(tau(s)>(obs)) displays an Arrhenius type temperature dependence. From the temperature variation, the activation energy and entropy of the solvation process are determined to be 12.5 kcal M(-1) and 28 cal M(-1) K(-1), respectively. The ultraslow component and its temperature dependence are ascribed to a dynamic exchange between bound and free water molecules.

Solvation Dynamics in Reverse Micelles: The Role of Headgroup−Solute Interactions

We present molecular dynamics simulation results for solvation dynamics in the water pool of anionic-surfactant reverse micelles (RMs) of varying water content, w(0). The model RMs are designed to represent water/aerosol-OT/oil systems, where aerosol-OT is the common name for sodium bis(2-ethylhexyl)sulfosuccinate. To determine the effects of chromophore-headgroup interactions on solvation dynamics, we compare the results for charge localization in model ionic diatomic chromophores that differ only in charge sign. Electronic excitation in both cases is modeled as charge localization on one of the solute sites. We find dramatic differences in the solvation responses for anionic and cationic chromophores. Solvation dynamics for the cationic chromophore are considerably slower and more strongly w(0)-dependent than those for the anionic chromophore. Further analysis indicates that the difference in the responses can be ascribed in part to the different initial locations of the two chromophores relative to the surfactant interface. In addition, slow motion of the cationic chromophore relative to the interface is the main contributor to the longer-time decay of the solvation response to charge localization in this case.