Aqueous solvation of amphiphilic solutes: concentration and temperature dependent study of the ultrafast polarisability relaxation dynamics

An understanding of the influence of hydrophilic and hydrophobic interactions on the dynamics of solvating water molecules is important in a diverse range of phenomena. The polarisability anisotropy relaxation dynamics of aqueous solutions of the amphiphiles TBA (t-butyl alcohol) and TMAO (trimethylamine N-oxide) have been measured as a function of concentration and temperature. TMAO is shown to have a greater effect on the picosecond relaxation dynamics of water than TBA. This result is consistent with hydrophilic interactions being mainly responsible for the slowing down the polarisability relaxation in aqueous solutions. The room temperature Raman spectral densities of the two solutions are remarkably similar to that of bulk water, an effect which is tentatively ascribed to the formation of nanoscale structure in the solutions, allowing the formation of bulk-like water pools. The temperature dependent spectral density of TMAO remains similar to that of bulk water at all temperatures, while that for TBA shows a marked decrease in the amplitude of the response usually ascribed to a water-water stretch with increasing temperature. This is discussed in terms of the temperature dependent structure of TBA aggregates in solution.     

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.     

A surface science approach to ultrafast electron transfer and solvation dynamics at interfaces

Excess electrons in polar media, such as water or ice, are screened by reorientation of the surrounding molecular dipoles. This process of electron solvation is of vital importance for various fields of physical chemistry and biology as, for instance, in electrochemistry or photosynthesis. Generation of such excess electrons in bulk water involves either photoionization of solvent molecules or doping with e.g. alkali atoms, involving possibly perturbing interactions of the system with the parent-cation. Such effects are avoided when using a surface science approach to electron solvation: in the case of polar adsorbate layers on metal surfaces, the substrate acts as an electron source from where photoexcited carriers are injected into the adlayer. Besides the investigation of electron solvation at such interfaces, this approach allows for the investigation of heterogeneous electron transfer, as the excited solvated electron population continuously decays back to the metal substrate. In this manner, electron transfer and solvation processes are intimately connected at any polar adsorbate-metal interface. In this tutorial review, we discuss recent experiments on the ultrafast dynamics of photoinduced electron transfer and solvation processes at amorphous ice-metal interfaces. Femtosecond time-resolved two-photon photoelectron spectroscopy is employed as a direct probe of the electron dynamics, which enables the analysis of all elementary processes: the charge injection across the interface, the subsequent electron localization and solvation, and the dynamics of electron transfer back to the substrate. Using surface science techniques to grow and characterize various well-defined ice structures, we gain detailed insight into the correlation between adsorbate structure and electron solvation dynamics, the location (bulk versus surface) of the solvation site, and the role of the electronic structure of the underlying metal substrate on the electron transfer rate.     

Molecular dynamics study of the solvation of an alpha-helical transmembrane peptide by DMSO

A 10-ns molecular dynamics study of the solvation of a hydrophobic transmembrane helical peptide in dimethyl sulfoxide (DMSO) is presented. The objective is to analyze how this aprotic polar solvent is able to solvate three groups of amino acid residues (i.e., polar, apolar, and charged) that are located in a stable helical region of a transmembrane peptide. The 25-residue peptide (sMTM7) used mimics the cytoplasmic proton hemichannel domain of the seventh transmembrane segment (TM7) from subunit a of H(+)-V-ATPase from Saccharomyces cerevisiae. The three-dimensional structure of peptide sMTM7 in DMSO has been previously solved by NMR spectroscopy. The radial and spatial distributions of the DMSO molecules surrounding the peptide as well as the number of hydrogen bonds between DMSO and the side chains of the amino acid residues involved are extracted from the molecular dynamics simulations. Analysis of the molecular dynamics trajectories shows that the amino acid side chains are fully embedded in DMSO. Polar and positively charged amino acid side chains have dipole-dipole interactions with the oxygen atom of DMSO and form hydrogen bonds. Apolar residues become solvated by DMSO through the formation of a hydrophobic pocket in which the methyl groups of DMSO are pointing toward the hydrophobic side chains of the residues involved. The dual solvation properties of DMSO cause it to be a good membrane-mimicking solvent for transmembrane peptides that do not unfold due to the presence of DMSO.     

Quantitative distinction between competing intramolecular bond twisting and solvent relaxation dynamics: an ultrafast study

Often an intramolecular relaxation process takes place in a time scale similar to that of the solvent relaxation process. Under these circumstances the dynamic Stokes' shift of the probe can be modulated by the combined effect of these two relaxation processes. In the present article we have studied ultrafast solvent relaxation using three different coumarin dyes and proposed a methodology for the quantitative separation of the dynamics of two competing processes, namely, solvent relaxation and bond twisting, that take place simultaneously in the present systems.     

2-Naphthyl(carbomethoxy)carbene revisited: combination of ultrafast UV-vis and infrared spectroscopic study

Ultrafast laser flash photolysis (310 nm) of methyl 2-napthyldiazoacetate (2-NpCN2CO2CH3) in acetonitrile or cyclohexane produces a diazo excited state which absorbs broadly in the visible region (tau = 300 fs). The decay of the excited diazo compound is accompanied by growth of the vibrationally excited singlet 2-naphthyl(carbomethoxy)carbene ((1)NpCCO2CH3). The singlet carbene absorbs at 360 and 470 nm. In acetonitrile these bands do not decay over 3 ns, but they do decay by approximately 50% of their original intensity in cyclohexane in 3 ns. It is concluded that (1)NpCCO2CH3 has a singlet ground state in acetonitrile but a triplet ground state in cyclohexane. Related experiments reveal a singlet ground state in Freon-113 and chloroform. This interpretation is supported by ultrafast IR spectroscopy, which confirms that only (1)NpCCO2CH3 is formed within 50 ps of the laser pulse rather than a singlet-triplet equilibrium mixture of carbene. The planar singlet relaxes to the preferred perpendicular singlet over a few tens of picoseconds, as evidenced by a red shift of the carbonyl stretching vibration. Although our data agrees with previous studies, its interpretation is somewhat altered.     

Water dynamics at protein interfaces: ultrafast optical Kerr effect study

The behavior of water molecules surrounding a protein can have an important bearing on its structure and function. Consequently, a great deal of attention has been focused on changes in the relaxation dynamics of water when it is located at the protein surface. Here we use the ultrafast optical Kerr effect to study the H-bond structure and dynamics of aqueous solutions of proteins. Measurements are made for three proteins as a function of concentration. We find that the water dynamics in the first solvation layer of the proteins are slowed by up to a factor of 8 in comparison to those in bulk water. The most marked slowdown was observed for the most hydrophilic protein studied, bovine serum albumin, whereas the most hydrophobic protein, trypsin, had a slightly smaller effect. The terahertz Raman spectra of these protein solutions resemble those of pure water up to 5 wt % of protein, above which a new feature appears at ~80 cm(-1), which is assigned to a bending of the protein amide chain.     

Characterization of the solvation dynamics of an ionic liquid via molecular dynamics simulation

The solvation dynamics of ionic liquids have been the subject of intense experimental study but remain poorly understood. We present the results of molecular dynamics simulations of the solvation dynamics of the ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate in response to photoexcitation of the fluorescent dye coumarin-153. We reproduce the time-resolved fluorescence Stokes shift using linear response theory, then use novel statistical techniques to analyze cation and anion contributions to the signal. We find that the solvation dynamics are dominated by collective ionic motion and characterize the time scale for various features of the collective response. Further, we use the Steele analysis [Mol. Phys. 61, 1031 (1987)] to characterize the contributions to the observed Stokes shift made by translational and rovibrational degrees of freedom. Our results indicate that in contrast to molecular liquids, the rovibrational response is trivial and the observed fluorescence response arises almost entirely from ionic translation. Our results resolve previously open questions in the literature about the nature of the rapid dynamics in room-temperature ionic liquids and offer insight into the physical principles governing ionic liquid behavior on longer time scales.     

A comparative study of solvation dynamics in room-temperature ionic liquids

The solvation dynamics of ionic liquids have been the subject of many experimental and theoretical studies but remain poorly understood. We analyze these dynamics by modeling the time-resolved fluorescence response of coumarin 153 in two room-temperature ionic liquids: 1-butyl-1-methylpyrrolidinium bromide and 1-butyl-1-methylpyrrolidinium bis(trifluoromethylsulfonyl)imide. Our results demonstrate that phenomena such as electrostatic screening operate significantly differently in the two liquids, and the relative importance of translational and rovibrational components of the ionic response depends significantly on the character of the ions involved. However, collective motion dominates the response of both ionic liquids, and the qualitative features of this collective behavior are strikingly similar in both cases.     

Molecular dynamics study of 2rotaxanes: influence of solvation and cation on co-conformation

The conformational preference of a [2]rotaxane system has been examined by molecular dynamics simulations. The rotaxane wheel consists of two bridged binding components: a cis-dibenzo-18-crown-6 ether and a 1,3-phenyldicarboxamide, and the penetrating axle consists of a central isophthaloyl unit with phenyltrityl capping groups. The influence of solvation on the co-conformation of the [2]rotaxane was evaluated by comparing the conformational flexibility in two solvents: chloroform and dimethyl sulfoxide. Attention was also paid to the effect of cation binding on the dynamical properties of the [2]rotaxane. The conformational stability of the [2]rotaxane was calculated using a MM/PB-SA strategy, and the occurrence of specific motions was examined by essential dynamics analysis. The changes in the co-conformational properties in the two solvents and upon cation binding are discussed in light of the available NMR data. The results indicate that in chloroform solution the [2]rotaxane system exists as a mixture of co-conformational states including some that have hydrogen bonds between axle C=O and wheel NH groups. Analysis of the simulations allow us to hypothesize that the [2]rotaxane's circumrotation motion can occur as the result of a dynamic process that combines a preliminary axle sliding step that breaks these hydrogen bonds and a conformational change in the ester group more distant from the wheel. In contrast, no hydrogen-bonded co-conformation was found in dimethyl sulfoxide, which appears to be due to the preferential formation of hydrogen bonds between the wheel NH groups with solvent molecules. Moreover, the axle experiences notable changes in anisotropic shielding, which would explain why the NMR signals are broadened in this solvent. Insertion of a sodium cation into the crown ether reduces co-conformational flexibility due to an interaction of the axle with the cation. Overall, the results reveal how both solvent and ionic atmosphere can influence the co-conformational preferences of rotaxanes.     

Ab initio study of ultrafast photochemical reaction dynamics of phenol blue

Phenol blue (PB) is a primary skeletal structure part of indoaniline dyes and well-known as a solvatochromic dye. It has been recently observed by pump-probe (PP) transient absorption measurements that PB shows ultrafast ground state recovery within a few hundred femtoseconds after photoexcitation. In this work, the ultrafast photochemical reaction mechanism of PB has been investigated using direct ab initio (CASSCF) nonadiabatic molecular dynamics with the trajectory surface hopping (TSH) method. The swarm of trajectories starting from the S1 Franck-Condon (FC) point has mostly shown surface hops (nonadiabatic transitions) from the S1 state to the S0 state at 110-120 fs in the vicinity of an S1/S0 conical intersection and after decay to the S0 state bifurcated into two (Reverse and Forward) directions with almost the same branching ratio and reached the vicinity of the S0 minimum energy point at 200-300 fs, which is in good agreement with the fast time component of the ground state recovery in the PP measurements. After reaching the vicinity of the S0 minimum energy point, the trajectories showed a coherent vibration of bending motion between quinoneimine and aniline rings with a low frequency of 43 cm-1, which presumably corresponds to a coherently photoexcitation-induced vibrational mode with a low frequency recently observed by the PP measurements.     

Molecular dynamics study of the solvation of calcium carbonate in water

We performed molecular dynamics simulations of diluted solutions of calcium carbonate in water. To this end, we combined and tested previous polarizable models. The carbonate anion forms long-living hydrogen bonds with water and shows an amphiphilic character, in which the water molecules are expelled in a region close to its C(3) symmetry axis. The calcium cation forms a strongly bound ion pair with the carbonate. The first hydration shell around the CaCO(3) pair is found to be very similar to the location of the water molecules surrounding CaCO(3) in ikaite, the hydrated mineral.     

Direct dynamics study of ultrafast vibrational energy relaxation in ice Ih

Car-Parrinello molecular dynamics (CPMD) and a previously developed wave packet model are used to study ultrafast relaxation in water clusters. Water clusters of 15 water molecules are used to represent ice Ih. The relaxation is studied by exciting a symmetric or an asymmetric stretch mode of the central water molecule. The CPMD results suggest that relaxation occurs within 100 fs. This is in agreement with experimental work by Woutersen and Bakker and the earlier wave packet calculations. The CPMD results further indicate that the excitation energy is transferred both intramolecularly and intermolecularly on roughly the same time scale. The intramolecular energy transfer occurs predominantly between the symmetric and asymmetric modes while the bend mode is largely left unexcited on the short time scale studied here.     

Structure and ultrafast dynamics of liquid water: a quantum mechanics/molecular mechanics molecular dynamics simulations study

A quantum mechanics/molecular mechanics molecular dynamics simulation was performed for liquid water to investigate structural and dynamical properties of this peculiar liquid. The most important region containing a central reference molecule and all nearest surrounding molecules (first coordination shell) was treated by Hartree-Fock (HF), post-Hartree-Fock [second-order Moller-Plesset perturbation theory (MP2)], and hybrid density functional B3LYP [Becke's three parameter functional (B3) with the correlation functional of Lee, Yang, and Parr (LYP)] methods. In addition, another HF-level simulation (2HF) included the full second coordination shell. Site to site interactions between oxygen-oxygen, oxygen-hydrogen, and hydrogen-hydrogen atoms of all ab initio methods were compared to experimental data. The absence of a second peak and the appearance of a shoulder instead in the gO-O graph obtained from the 2HF simulation is notable, as this feature has been observed so far only for pressurized or heated water. Dynamical data show that the 2HF procedure compensates some of the deficiency of the HF one-shell simulation, reducing the difference between correlated (MP2) and HF results. B3LYP apparently leads to too rigid structures and thus to an artificial slow down of the dynamics.     

Ar solvation shells in K(+)-HFBz: from cluster rearrangement to solvation dynamics

The effect of some leading intermolecular interaction components on specific features of weakly bound clusters involving an aromatic molecule, a closed shell ion, and Ar atoms is analyzed by performing molecular dynamics simulations on potential energy surfaces properly formulated in a consistent way. In particular, our investigation focuses on the three-dimensional Ar distributions around the K(+)-hexafluorobenzene (K(+)-HFBz) dimer, in K(+)-HFBz-Ar(n) aggregates (n ≤ 15), and on the gradual evolution from cluster rearrangement to solvation dynamics when ensembles of 50, 100, 200, and 500 Ar atoms are taken into account. Results indicate that the Ar atoms compete to be placed in such a way to favor an attractive interaction with both K(+) and HFBz, occupying positions above and below the aromatic plane but close to the cation. When these positions are already occupied, the Ar atoms tend to be placed behind the cation, at larger distances from the center of mass of HFBz. Accordingly, three different groups of Ar atoms are observed when increasing n, with two of them surrounding K(+), thus, disrupting the K(+)-HFBz equilibrium geometry and favoring the dissociation of the solvated cation when the temperature increases. The selective role of the leading intermolecular interaction components directly depending on the ion size repulsion is discussed in detail by analyzing similarities and differences on the behavior of the Ar-solvated K(+)-HFBz and Cl(-)-Bz aggregates.     

Polar solvation dynamics of lysozyme from molecular dynamics studies

The solvation dynamics of a protein are believed to be sensitive to its secondary structures. We have explored such sensitivity in this article by performing room temperature molecular dynamics simulation of an aqueous solution of lysozyme. Nonuniform long-time relaxation patterns of the solvation time correlation function for different segments of the protein have been observed. It is found that relatively slower long-time solvation components of the α-helices and β-sheets of the protein are correlated with lower exposure of their polar probe residues to bulk solvent and hence stronger interactions with the dynamically restricted surface water molecules. These findings can be verified by appropriate experimental studies.     

Femtosecond infrared study of the dynamics of solvation and solvent caging.

The ultrafast reaction dynamics following 295-nm photodissociation of Re2CO10 were studied experimentally with 300-fs time resolution in the reactive, strongly coordinating CCl4 solution and in the inert, weakly coordinating hexane solution. Density-functional theoretical (DFT) and ab initio calculations were used to further characterize the transient intermediates seen in the experiments. It was found that the quantum yield of the Re-Re bond dissociation is governed by geminate recombination on two time scales in CCl4, approximately 50 and approximately 500 ps. The recombination dynamics are discussed in terms of solvent caging in which the geminate Re(CO)5 pair has a low probability to escape the first solvent shell in the first few picoseconds after femtosecond photolysis. The other photofragmentation channel resulted in the equatorially solvated dirhenium nonacarbonyl eq-Re2(CO)9(solvent). Theoretical calculations indicated that a structural reorganization energy cost on the order of 6-7 kcal/mol might be required for the unsolvated nonacarbonyl to coordinate to a solvent molecule. These results suggest that for Re(CO)5 the solvent can be treated as a viscous continuum, whereas for the Re2(CO)9 the solvent is best described in molecular terms.     

Probing polar solvation dynamics in proteins: a molecular dynamics simulation analysis

Measurements of time-resolved Stokes shifts on picosecond to nanosecond time scales have been used to probe the polar solvation dynamics of biological systems. Since it is difficult to decompose the measurements into protein and solvent contributions, computer simulations are useful to aid in understanding the details of the molecular behavior. Here we report the analysis of simulations of the electrostatic interactions of the rest of the protein and the solvent with 11 residues of the immunoglobulin binding domain B1 of protein G. It is shown that the polar solvation dynamics are position-dependent and highly heterogeneous. The contributions due to interactions with the protein and with the solvent are determined. The solvent contributions are found to vary from negligible after a few picoseconds to dominant on a scale of hundreds of picoseconds. The origin for the latter is found to involve coupled hydration and protein conformational dynamics. The resulting microscopic picture demonstrates that a wide range of possibilities have to be considered in the interpretation of time-resolved Stokes shift measurements.