Computational studies of liquid water and diluted water in carbon tetrachloride

Molecular dynamics simulations were carried out to study solvent effects on the energetic and dynamical properties of water molecules in liquid water and in carbon tetrachloride (CCl4). In these studies, the free-energy profiles or potentials of mean force (PMF) for water dimers in both solvents were computed. The computed PMF results showed a stable minimum near 3 A for the O-O separation, with a minimum free energy of about -2.8 kcal/mol in CCl4, as compared to a value of -0.5 kcal/mol in liquid water. The difference in free energy in water as compared to that in CCl4 was expected and is the result of competition from surrounding water molecules that are capable of forming hydrogen bonds in the liquid water. This capability is absent in the diluted water found in CCl4. We found that the rotational motions of H2O/D2O were nonisotropic, with the out-of-plane vector correlation times in H2O/D2O varying from 5.6/5.8 ps at 250 K to 0.57/0.56 ps at 350 K and the corresponding OH/OD bond vectors varying from 6.5/7.7 ps to 0.75/0.75 ps. The results compare reasonably well to the available NMR experimental and computer simulation data on the same system (Farrar; Skinner; et al. J. Am. Chem. Soc. 2001, 123, 8047). For diluted water in CCl4, we found the computed rotational correlation times also were nonisotropic and much longer than the corresponding NMR experimental values at the same concentration (Farrar; et al. J. Phys. Chem. A 2007, 111, 6146). Upon analyzing the water hydrogen-bonding patterns as a function of water concentration, we conclude that the differences in the rotational correlation times mainly result from the formation of water hydrogen-bonding networks as the water concentration is increased in liquid CCl4. In addition, we found the rotational correlation times to be substantially faster in liquid CCl4 than in liquid water.     

Rotational motion in liquid water is anisotropic: a nuclear magnetic resonance and molecular dynamics simulation study

Experimental NMR measurements of the deuterium and (17)O T(1) relaxation times in deuterium-enriched liquid water have been performed from 275 to 350 K. These relaxation times can yield rotational correlation times of appropriate molecule-fixed unit vectors if the quadrupole coupling constants and asymmetry parameters are known. We determine the latter from ab initio studies of water clusters and experimental chemical shift measurements. We find that the rotational correlation time for the OD bond vector in D(2)(16)O varies from 5.8 ps at 275 K to 0.86 ps at 350 K, and that the rotational correlation time for the out-of-plane vector of dilute D(2)(17)O in D(2)(16)O varies from 4.4 ps at 275 K to 0.64 ps at 350 K. These results indicate that the rotational motion of water is anisotropic. Molecular dynamics simulations of liquid water are in good agreement with these experiments at the higher temperatures, but the simulation results are considerably faster than experiment at the lower temperatures.     

17O nuclear quadrupole coupling constants of water bound to a metal ion: a gadolinium(III) case study

Rotational correlation times of metal ion aqua complexes can be determined from 17O NMR relaxation rates if the quadrupole coupling constant of the bound water oxygen-17 nucleus is known. The rotational correlation time is an important parameter for the efficiency of Gd3+ complexes as magnetic resonance imaging contrast agents. Using a combination of density functional theory with classical and Car-Parrinello molecular dynamics simulations we performed a computational study of the 17O quadrupole coupling constants in model aqua ions and the [Gd(DOTA)(H2O)]- complex used in clinical diagnostics. For the inner sphere water molecule in the [Gd(DOTA)(H2O)]- complex the determined quadrupole coupling parameter chi square root of (1 + eta2/3) of 8.7 MHz is very similar to that of the liquid water (9.0 MHz). Very close values were also predicted for the the homoleptic aqua ions of Gd3+ and Ca2+. We conclude that the 17O quadrupole coupling parameters of water molecules coordinated to closed shell and lanthanide metal ions are similar to water molecules in the liquid state.     

Conformational equilibrium isotope effects in glucose by (13)C NMR spectroscopy and computational studies.

Anomeric equilibrium isotope effects for dissolved sugars are required preludes to understanding isotope effects for these molecules bound to enzymes. This paper presents a full molecule study of the alpha- and beta-anomeric forms of D-glucopyranose in water using deuterium conformational equilibrium isotope effects (CEIE). Using 1D (13)C NMR, we have found deuterium isotope effects of 1.043 +/- 0.004, 1.027 +/- 0.005, 1.027 +/- 0.004, 1.001 +/- 0.003, 1.036 +/- 0.004, and 0.998 +/- 0.004 on the equilibrium constant, (H/D)K(beta/alpha), in [1-(2)H]-, [2-(2)H]-, [3-(2)H]-, [4-(2)H]-, [5-(2)H]-, and [6,6'-(2)H(2)]-labeled sugars, respectively. A computational study of the anomeric equilibrium in glucose using semiempirical and ab initio methods yields values that correlate well with experiment. Natural bond orbital (NBO) analysis of glucose and dihedral rotational equilibrium isotope effects in 2-propanol strongly imply a hyperconjugative mechanism for the isotope effects at H1 and H2. We conclude that the isotope effect at H1 is due to n(p) --> sigma* hyperconjugative transfer from O5 to the axial C1--H1 bond in beta-glucose, while this transfer makes no contribution to the isotope effect at H5. The isotope effect at H2 is due to rotational restriction of OH2 at 160 degrees in the alpha form and 60 degrees in the beta-sugar, with concomitant differences in n --> sigma* hyperconjugative transfer from O2 to CH2. The isotope effects on H3 and H5 result primarily from syn-diaxial steric repulsion between these and the axial anomeric hydroxyl oxygen in alpha-glucose. Therefore, intramolecular effects play an important role in isotopic perturbation of the anomeric equilibrium. The possible role of intermolecular effects is discussed in the context of recent molecular dynamics studies on aqueous glucose.     

How does internal motion influence the relaxation of the water protons in Ln(III)DOTA-like complexes?

Taking advantage of the Curie contribution to the relaxation of the protons in the Tb(III) complex, and the quadrupolar relaxation of the 17O and 2H nuclei on the Eu(III) complex, the effect of the internal motion of the water molecule bound to [Ln(DOTAM)(H2O)]3+ complexes was quantified. The determination of the quadrupolar coupling constant of the bound water oxygen chi(Omicron)(1 + eta(Omicron)2/3)1/2 = 5.2 +/- 0.5 MHz allows a new analysis of the 17O and 1H NMR data of the [Gd(DOTA)(H2O)]- complex with different rotational correlation times for the Gd(III)-O(water) and Gd(III)-H(water) vectors. The ratio of the rotational correlation times for the Ln(III)-H(water) vector and the overall rotational correlation time is calculated tau(RH)/tau(RO) = 0.65 +/- 0.2. This could have negative consequences on the water proton relaxivity, which we discuss in particular for macromolecular systems. It appears that the final effect is actually attenuated and should be around 10% for such large systems undergoing local motion of the chelating groups.     

Chirality transfer through hydrogen-bonding: experimental and ab initio analyses of vibrational circular dichroism spectra of methyl lactate in water

The infrared vibrational absorption (VA) and vibrational circular dichroism (VCD) spectra of methyl lactate were measured in the 1000-1800 cm(-1) region in the CCl(4) and H(2)O solvents, respectively. In particular, the chirality transfer effect, i.e. the H-O-H bending bands of the achiral water subunits that are hydrogen-bonded to the methyl lactate molecule exhibit substantial VCD strength, was detected experimentally. A series of density functional theory calculations using B3PW91 and B3LYP functionals with 6-311++G(d,p) and aug-cc-pVTZ basis sets were carried out to simulate the VA and VCD spectra of the methyl lactate monomer and the methyl lactate-(H(2)O)(n) complexes with n = 1, 2, 3. The population weighted VA and VCD spectra of the methyl lactate monomer are in good agreement with the experimental spectra in CCl(4). Implicit polarizable continuum model was found to be inadequate to account for the hydrogen-bonding effect in the observed VA and VCD spectra in H(2)O. The methyl lactate-(H(2)O)(n) complexes with n = 1, 2, 3 were used to model the explicit hydrogen-bonding. The population weighted VA and VCD spectra of the methyl lactate-H(2)O binary complex are shown to capture the main spectral features in the observed spectra in aqueous solution. The theoretical modeling shows that the extent of chirality transfer depends sensitively on the specific binding sites taken by the achiral water molecules. The observation of chirality transfer effect opens a new spectral window to detect and to model the hydrogen-bonding solvent effect on VCD spectra of chiral molecules.     

Dimerization and double proton transfer-induced tautomerism of 4(3H)-pyrimidinone in solution studied by IR spectroscopy and quantum chemical calculations

The tautomerism and dimerization of 4(3H)-pyrimidinone (4(3H)Pyr) in carbon tetrachloride (CCl4) and chloroform (CHCl3) solutions were investigated using IR spectroscopy and quantum chemical calculations. The observed IR spectra in the NH and OH stretching regions clearly revealed the predominance of the keto tautomer in both solvent systems. The enol form only exists in a very small proportion in the CCl4 solution. The tautomeric constant for the two monomers KT[OH/NH] = 0.012 and DeltaE = 2.62 kcal/mol were estimated at 25 degrees C. This result was supported by the self-consistent reaction field/polarizable continuum (SCRF/PCM) calculation at the MP4(full, SDQ)/aug-cc-pVDZ level, which predicted DeltaE = 3.06 kcal/mol in CCl4. In the C=O stretching region, two bands were observed, suggesting the coexistence of two keto structures at equilibrium. The calculated IR spectra indicated that the bands at 1711 and 1675 cm(-1) arise from the keto monomer and keto-keto (KK) ring dimer, respectively. At elevated temperature, the populations of both the keto and enol monomers increased for the CCl4 solution. The present study revealed that the keto <--> enol tautomerization does not occur in the isolated monomer molecule. The double proton transfer (DPT) reaction in the KK ring dimer presumably plays a substantial role in the population increase of the enol monomer. To our knowledge, this may be the first observation of the tautomerization in a model base pair via the temperature-induced ground-state DPT reaction under a nonpolar liquid environment reported so far. This tautomerism can serve as a mimic circumstance for the spontaneous mutations induced by proton transfer in the DNA base pairs.     

On-line determination of the deuterium abundance in breath water vapour by flowing afterglow mass spectrometry with applications to measurements of total body water

We have developed a new method for the on-line quantification of deuterium in water vapour. We call this method flowing afterglow mass spectrometry (FA-MS). A swarm of H3O+ precursor ions is created in flowing helium carrier gas by a microwave discharge. These precursor ions react with the H2O, HDO, H2(17)O and H2(18)O molecules in a water vapour sample that is introduced into the carrier gas/H3O+ ion swarm. The hydrated ions, H3O+.(H2O)3 at m/z 73, and their isotopic variant ions H8DO4(+) and H9(17)OO(3)(+) at m/z 74 and H9(18)OO(3)(+) at m/z 75, are thus formed. By adopting the known fractional abundance of 18O in water vapour, and accounting for the contribution of the isotopic ions H9(17)OO(3)(+) to the ion signal at m/z 74, a measurement of the 74/75 ion signal ratio under equilibrium conditions provides the fractional deuterium abundance in the water vapour sample. Using this technique, the deuterium abundance in the water vapour present in single exhalations of breath can be determined. Thus, from the temporal variations of breath deuterium following the ingestion of a known quantity of D(2)O, we show that total body water can be determined non-invasively and the kinetics of water flow around the body can be tracked.     

Treatment of dilute clusters of methanol and water by ab initio quantum mechanical calculations

Large molecular clusters can be considered as intermediate states between gas and condensed phases, and information about them can help us understand condensed phases. In this paper, ab initio quantum mechanical methods have been used to examine clusters formed of methanol and water molecules. The main goal was to obtain information about the intermolecular interactions and the structure of methanol/water clusters at the molecular level. The large clusters (CH(4)O...(H(2)O)(12) and H(2)O...(CH(4)O)(10)) containing one molecule of one component (methanol or water) and many (12, 10) molecules of the other component were considered. M?ller-Plesset perturbation theory (MP2) was used in the calculations. Several representative cluster geometries were optimized, and nearest-neighbor interaction energies were calculated for the geometries obtained in the first step. The results of the calculations were compared to the available experimental information regarding the liquid methanol/water mixtures and to the molecular dynamics and Monte Carlo simulations, and good agreement was found. For the CH(4)O...(H(2)O)(12) cluster, it was shown that the molecules of water can be subdivided into two classes: (i) H bonded to the central methanol molecule and (ii) not H bonded to the central methanol molecule. As expected, these two classes exhibited striking energy differences. Although they are located almost the same distance from the carbon atom of the central methanol molecule, they possess very different intermolecular interaction energies with the central molecule. The H bonding constitutes a dominant factor in the hydration of methanol in dilute aqueous solutions. For the H(2)O...(CH(4)O)(10) cluster, it was shown that the central molecule of water has almost three H bonds with the methanol molecules; this result differs from those in the literature that concluded that the average number of H bonds between a central water molecule and methanol molecules in dilute solutions of water in methanol is about two, with the water molecules being incorporated into the chains of methanol. In contrast, the present predictions revealed that the central water molecule is not incorporated into a chain of methanol molecules, but it can be the center of several (2-3) chains of methanol molecules. The molecules of methanol, which are not H bonded to the central water molecule, have characteristics similar to those of the methane molecules around a central water molecule in the H(2)O...(CH(4))(10) cluster. The ab initio quantum mechanical methods employed in this paper have provided detailed information about the H bonds in the clusters investigated. In particular, they provided full information about two types of H bonds between water and methanol molecules (in which the water or the methanol molecule is the proton donor), including information about their energies and lengths. The average numbers of the two types of H bonds in the CH(4)O...(H(2)O)(12) and H(2)O...(CH(4)O)(10) clusters have been calculated. Such information could hardly be obtained with the simulation methods.     

Rotational dynamics in supercooled water from nuclear spin relaxation and molecular simulations

Structural dynamics in liquid water slow down dramatically in the supercooled regime. To shed further light on the origin of this super-Arrhenius temperature dependence, we report high-precision (17)O and (2)H NMR relaxation data for H(2)O and D(2)O, respectively, down to 37 K below the equilibrium freezing point. With the aid of molecular dynamics (MD) simulations, we provide a detailed analysis of the rotational motions probed by the NMR experiments. The NMR-derived rotational correlation time τ(R) is the integral of a time correlation function (TCF) that, after a subpicosecond librational decay, can be described as a sum of two exponentials. Using a coarse-graining algorithm to map the MD trajectory on a continuous-time random walk (CTRW) in angular space, we show that the slowest TCF component can be attributed to large-angle molecular jumps. The mean jump angle is ～48° at all temperatures and the waiting time distribution is non-exponential, implying dynamical heterogeneity. We have previously used an analogous CTRW model to analyze quasielastic neutron scattering data from supercooled water. Although the translational and rotational waiting times are of similar magnitude, most translational jumps are not synchronized with a rotational jump of the same molecule. The rotational waiting time has a stronger temperature dependence than the translation one, consistent with the strong increase of the experimentally derived product τ(R)?D(T) at low temperatures. The present CTRW jump model is related to, but differs in essential ways from the extended jump model proposed by Laage and co-workers. Our analysis traces the super-Arrhenius temperature dependence of τ(R) to the rotational waiting time. We present arguments against interpreting this temperature dependence in terms of mode-coupling theory or in terms of mixture models of water structure.     

Synthesis and characterization of the novel monoamide derivatives of Gd-TTDA

For this study, the N'-monoamide derivatives of TTDA (3,6,10-tri(carboxymethyl)-3,6,10-triazadodecanedioic acid), N'-methylamide (TTDA-MA), N'-benzylamide (TTDA-BA), and N'-2-methoxybenzylamide (TTDA-MOBA), were synthesized. Their protonation constants and stability constants (log K(ML)'s) formed with Ca(2+), Zn(2+), Cu(2+), and Gd(3+) were determined by potentiometric titration in 0.10 M Me(4)NCl at 25.0 +/- 0.1 degrees C. The relaxivity values of [Gd(TTDA-MA)](-), [Gd(TTDA-BA)](-), and [Gd(TTDA-MOBA)](-) remained constant with respect to pH changes over the range 4.5-12.0. The (17)O NMR chemical shift of H(2)O induced by [Dy(TTDA-MA)(H(2)O)](-) at pH 6.80 showed 0.9 inner-sphere water molecules. Water proton relaxivity values for [Gd(TTDA-MA)(H(2)O)](-), [Gd(TTDA-BA)(H(2)O)](-), and [Gd(TTDA-MOBA)(H(2)O)](-) at 37.0 +/- 0.1 degrees C and 20 MHz are 3.89, 4.21, and 4.25, respectively. The water-exchange lifetime (tau(M)) and rotational correlation time (tau(R)) of [Gd(TTDA-MA)(H(2)O)](-), [Gd(TTDA-BA)(H(2)O)](-), and [Gd(TTDA-MOBA)(H(2)O)](-) are obtained from reduced the (17)O relaxation rate and chemical shifts of H(2)(17)O. The (2)H NMR longitudinal relaxation rates of the deuterated diamagnetic lanthanum complexes for the rotational correlation time were also thoroughly investigated. The water-exchange rates (K(298)(ex) for [Gd(TTDA-MA)(H(2)O)](-), [Gd(TTDA-BA)(H(2)O)](-), and [Gd(TTDA-MOBA)(H(2)O)](-) are lower than that of [Gd(TTDA)(H(2)O)](2)(-) but significantly higher than those of [Gd(DTPA)(H(2)O)](2)(-) and [Gd(DTPA-BMA)(H(2)O)]. The rotational correlation times for [Gd(TTDA-BA)(H(2)O)](-) and [Gd(TTDA-MOBA)(H(2)O)](-) are significantly longer than those of [Gd(TTDA)(H(2)O)](2)(-) and [Gd(DTPA)(H(2)O)](2)(-) complexes. The marked increase of the relaxivity of [Gd(TTDA-BA)(H(2)O)](-) and [Gd(TTDA-MOBA)(H(2)O)](-) results mainly from their longer rotational correlation time. The noncovalent interaction between human serum albumin (HSA) and [Gd(TTDA-BA)(H(2)O)](-) and [Gd(TTDA-MOBA)(H(2)O)](-) complexes containing a hydrophobic substituent was investigated by measuring the water proton relaxation rate of the aqueous solutions. The binding association constant (K(A)) values are 1.0 +/- 0.2 x 10(3) and 1.3 +/- 0.2 x 10(3) M(-1) for [Gd(TTDA-BA)(H(2)O)](-) and [Gd(TTDA-MOBA)(H(2)O)](-), which indicates a stronger interaction of [Gd(TTDA-BA)(H(2)O)](-) and [Gd(TTDA-MOBA)(H(2)O)](-) with HSA.     

Supramolecular assembly of an amphiphilic Gd(III) chelate: tuning the reorientational correlation time and the water exchange rate

We report the synthesis and characterization of the novel ligand H(5)EPTPA-C(16) ((hydroxymethylhexadecanoyl ester)ethylenepropylenetriaminepentaacetic acid). This ligand was designed to chelate the Gd(III) ion in a kinetically and thermodynamically stable way while ensuring an increased water exchange rate (kappa(ex)) on the Gd(III) complex owing to steric compression around the water-binding site. The attachment of a palmitic ester unit to the pendant hydroxymethyl group on the ethylenediamine bridge yields an amphiphilic conjugate that forms micelles with a long tumbling time (tau(R)) in aqueous solution. The critical micelle concentration (cmc = 0.34 mM) of the amphiphilic [Gd(eptpa-C(16))(H(2)O)](2-) chelate was determined by variable-concentration proton relaxivity measurements. A global analysis of the data obtained in variable-temperature and multiple-field (17)O NMR and (1)H NMRD measurements allowed for the determination of parameters governing relaxivity for [Gd(eptpa-C(16))(H(2)O)](2-); this is the first time that paramagnetic micelles with optimized water exchange have been investigated. The water exchange rate was found to be kappa(298)(ex) = 1.7 x 10(8) s(-1), very similar to that previously reported for the nitrobenzyl derivative [Gd(eptpa-bz-NO(2))(H(2)O)](2-) kappa(298)(ex) = 1.5 x 10(8) s(-1)). The rotational dynamics of the micelles were analysed by using the Lipari-Szabo approach. The micelles formed in aqueous solution show considerable flexibility, with a local rotational correlation time of tau(298)(l0) = 330 ps for the Gd(III) segments, which is much shorter than the global rotational correlation time of the supramolecular aggregates, tau(298)(g0) = 2100 ps. This internal flexibility of the micelles is responsible for the limited increase of the proton relaxivity observed on micelle formation (r(1) = 22.59 mM(-1) s(-1) for the micelles versus 9.11 mM(-1) s(-1) for the monomer chelate (20 MHz; 25 degrees C)).     

Novel synthetic method of phosphonamidate peptides and its application in peptide sequencing via multistage mass spectrometry

Phosphonamidate peptides were prepared in good yields by reaction of ethoxy(phenyl)phosphinate with free peptides in a mixed solution of H2O, C2H5OH, Et3N and CCl4 at room temp.; their in situ electrospray ionization mass spectra exhibited high sensitivity compared with the free peptides, and sequential loss of amino acid residues of the sodiated phosphonamidate peptides from the C-terminus was found in multistage ESI mass spectrometry. The results show that N-phosphonyl derivatization combined with multistage ESI-MS is a powerful method for peptide sequencing.     

Molecular dynamics simulation of [Gd(egta)(H(2)O)](-) in aqueous solution: internal motions of the poly(amino carboxylate) and water ligands, and rotational correlation times

Molecular dynamics simulations of [Gd(egta)(H(2)O)](-) (egta(4-)=3,12-bis(carboxymethyl)-6,9-dioxa-3,12-diazatetradecanedioate(4-)) have been performed without any artificial constraint on the first coordination sphere, such as covalent bonds between the Gd(3+) and the coordination sites. Two new crystallographic structures were determined for this gadolinium chelate and used to start two molecular dynamics simulations. [Gd(egta)(H(2)O)](-) and [Gd(egta)](-) were both observed during the simulations, with a mean volume for the reaction of dissociation [Gd(egta)(H(2)O)](-)-->[Gd(egta)](-)+H(2)O of +7.2 cm(3)mol(-1), which corroborates the previously published experimental value of +10.5 cm(3)mol(-1). Changes in the conformation of the complex with the inversion of several dihedral angles are observed in the simulations independently from the water dissociation. Very fast changes of the third-order rotation axis direction of the Gd(3+) coordination polyhedron (of symmetry D(3h)) are observed during the simulations and are related to the mechanism of electronic relaxation of the complex. Different rotational correlation times (tau(R)) were calculated from the simulations on various observables of the complex. Protons of the inner sphere have different tau(R). The mean tau(R) of the two Gd-HW(HW=hydrogen of water molecule) vectors is 72% lower than tau(R) of the complex, and 75% lower than tau(R) of the vector Gd-OW (OW=oxygen of water molecule). This discrimination of the tumbling rates should be taken into account in future global (17)O NMR, EPR and NMRD (nuclear magnetic relaxation dispersion) data analysis.     

Solvent (1)H NMRD study of hexaaquochromium(III): inferences on hydration and electron relaxation

The water proton nuclear magnetic relaxation dispersions (NMRD) of hexaaquochromium(III) in water and in water-glycerol solutions were obtained at several temperatures and viscosities. The data were analyzed in terms of the available theories by taking into account the contributions from first sphere, second sphere, and outer sphere water molecules. A meaningful analysis was possible by taking advantage of the structural model obtained from (17)O relaxation data, which was recently made available in the literature (Bleuzen, A.; Foglia, F.; Furet, E.; Helm, L.; Merbach, A.; Weber, J. J. Am. Chem. Soc. 1996, 118, 12 777). Dynamic parameters, like the molecular rotational time, the exchange time of the water protons of the first coordination sphere, the correlation time for electron relaxation, and the magnetic field dependence of electron relaxation were obtained. The possible contribution to water proton relaxivity of second sphere water molecules for some other hexaaqua complexes is also discussed.     

Water exchange in fluoroaluminate complexes in aqueous solution: a variable temperature multinuclear NMR study

An 17O, 19F, and 27Al NMR study of fluoroaluminate complexes (AlFn(H2O)6-n((3-n)+), n = 0, 1, and 2) in aqueous solution supports the idea that for each substitution of a bound water molecule by a fluoride anion, the exchange rate of bound water with free water increases by about 2 orders of magnitude. New rate coefficients for exchange of inner-sphere water molecules in AlF(H2O)5(2+) are kex(298) = 230(+/-20) s(-1), DeltaH(dagger) = 65(+/-3) kJ mol(-1), and DeltaS(dagger) = 19(+/-10) J mol(-1) K(-1). The corresponding new values for the AlF2(H2O)4(+) complex are: kex(298) = 17 100(+/-500) s(-1), DeltaH(dagger) = 66(+/-2) kJ mol(-1), and DeltaS(dagger) = 57(+/-8) J mol(-1) K(-1). When these new results are combined with those of our previous study,(4) we find no dependence of the solvent exchange rate, in either AlF(H2O)5(2+) or AlF2(H2O)4(+), on the concentration of fluoride or protons over the range of SigmaF = 0.06-0.50 M and [H(+)] = 0.01-0.44 M. A paramagnetic shift of 27Al resonances results from addition of Mn(II) to the aqueous solution as a relaxation agent for bulk waters. This shift allows resolution of the AlFn(H2O)6-n((3-n)+) species in 27Al NMR spectra and comparison of the speciation determined via thermodynamic calculations with that determined by 27Al, 19F, and 17O NMR.     

Unexpected aggregation of neutral, xylene-cored dinuclear GdIII chelates in aqueous solution

We have synthesized ditopic ligands L(1), L(2), and L(3) that contain two DO3A(3-) metal-chelating units with a xylene core as a noncoordinating linker (DO3A(3-) = 1,4,7,10-tetraazacyclododecane-1,4,7-triacetate; L(1) = 1,4-bis{[4,7,10-tris(carboxymethyl)-1,4,7,10-tetraazacyclododecane-1-yl]methyl}benzene; L(2) = 1,3-bis{[4,7,10-tris(carboxymethyl)-1,4,7,10-tetraazacyclododecane-1-yl]methyl}benzene; L(3) = 3,5-bis{[4,7,10-tris(carboxymethyl)-1,4,7,10-tetraazacyclododecane-1-yl]methyl}benzoic acid). Aqueous solutions of the dinuclear Gd(III) complexes formed with the three ligands have been investigated in a variable-temperature, multiple-field (17)O NMR and (1)H relaxivity study. The (17)O longitudinal relaxation rates measured for the [Gd(2)L(1-3)(H2O)(2)] complexes show strong field dependence (2.35-9.4 T), which unambiguously proves the presence of slowly tumbling entities in solution. The proton relaxivities of the complexes, which are unexpectedly high for their molecular weight, and in particular the relaxivity peaks observed at 40-50 MHz also constitute experimental evidences of slow rotational motion. This was explained in terms of self-aggregation related to hydrophobic interactions, pi stacking between the aromatic linkers, or possible hydrogen bonding between the chelates. The longitudinal (17)O relaxation rates of the [Gd(2)L(1-3)(H2O)(2)] complexes have been analysed with the Lipari-Szabo approach, leading to local rotational correlation times tau(1)(298) of 150-250 ps and global rotational correlation times tau(g)(298) of 1.6-3.4 ns (c(Gd): 20-50 mM), where tau(1)(298) is attributed to local motions of the Gd segments, while tau(g)(298) describes the overall motion of the aggregates. The aggregates can be partially disrupted by phosphate addition; however, at high concentrations phosphate interferes in the first coordination sphere by replacing the coordinated water. In contrast to the parent [Gd(DO3A)(H2O)(1.9)], which presents a hydration equilibrium between mono- and dihydrated species, a hydration number of q = 1 was established for the [Ln(2)L(1-3)(H2O)(2)] chelates by (17)O chemical shift measurements on Ln = Gd and UV/Vis spectrophotometry for Ln = Eu. The exchange rate of the coordinated water is higher for [Gd(2)L(1-3)(H2O)(2)] complexes k(ex)(298) = 7.5-12.0 x 10(6) s(-1)) than for [Gd(DOTA)(H2O)](-). The proton relaxivity of the [Gd(2)L(1-3)(H2O)(2)] complexes strongly decreases with increasing pH. This is related to the deprotonation of the inner-sphere water, which has also been characterized by pH potentiometry. The protonation constants determined for this process are logK(OH) = 9.50 and 10.37 for [Gd(2)L(1)(H2O)(2)] and [Gd(2)L(3)(H2O)(2)], respectively.     

Unexpected adsorption behavior of nonionic surfactants from glycol solvents

Adsorption and interfacial properties of model methyl-capped nonionic surfactants C8E4OMe [C8H17O(C2H4O)4CH3] and C10E4OMe [C10H21O(C2H4O)4CH3] were studied in water and water/ethylene glycol mixtures as well as pure ethylene glycol. Critical micellar concentrations (cmc's), surface tensions, and surface excess were determined using surface tension (ST) and neutron reflection (NR) as a function of solvent type and surfactant tail length. The ST results show a strong dependence on solvent type in terms of cmc. The NR data were analyzed using a single-layer model for the adsorbed surfactant films. Surprisingly, the adsorption parameters obtained in both water and pure ethylene glycol were very similar, and variations in film thickness or area per molecule are negligible in respect of the uncertainties. Similarly, for C10E4OMe, estimates for the free energies of adsorption and micellization show only a weak solvent dependence. These results suggest that for such model nonionic surfactants dilute solution properties are dictated by solvophobicity, which is quite similar for this class of water, glycol, and water-glycol mixtures. More specifically, the nature of the adsorption layer appears to be hardly affected by the type of solvent subphase. The findings highlight the significance of solvophobicity and show that model nonionic surfactants can behave very similarly in hydrogen-bonding glycol solvents and water.     

Reactions associated with ionization in water: a direct ab initio dynamics study of ionization in (H2O)17

Quasiclassical ab initio simulations of the ionization dynamics in a (H(2)O)(17) cluster, the first water cluster that includes a fourfold coordinated (internally solvated) water molecule, have been carried out to obtain a detailed picture of the elementary processes and energy redistribution induced by ionization in a model of aqueous water. General features observable from the simulations are the following: (i) well within 100 fs following the ionization, one or more proton transfers are seen to take place from the "ionized molecule" to neighboring molecules and beyond, forming a hydronium ion and a hydroxyl radical; (ii) two water molecules close to the ionized water molecule play an important role in the reaction, in what we term a "reactive trimer." The reaction time is gated by the encounter of the ionized water molecule with these two neighboring molecules, and this occurs anytime between 10 and 50 fs after the ionization. The distances of approach between the ionized molecule and the neighboring molecules indeed display best the time characteristics of the transfer of a proton, and thus of the formation of a hydronium ion and a OH radical. These findings are consistent with those for smaller cyclic clusters, albeit the dynamics of the proton transfer displays more varieties in the larger cluster than in the small cyclic clusters. We used a partitioning scheme for the kinetic energy in the (H(2)O)(17) system that distinguishes between the reactive trimer and the surrounding "medium." The analysis of the simulations indicates that the kinetic energy of the surrounding medium increases markedly right after the event of ionization, a manifestation of the local heating of the medium. The increase in kinetic energy is consistent with a reorganization of the surrounding medium, electrostatically forced in a very short time by the water cation and in a longer time by the formation of the hydronium ion.     

Properties of ternary insulating systems: the electronic structure of MgSO4.H2O

Structural and electronic properties of (100)-oriented MgSO(4) and MgSO(4).H(2)O surfaces and the adsorption of water on the latter were investigated theoretically with a combination of ab initio and semiempirical methods. Ab initio electronic structure calculations were based on a density functional theory (DFT)-Hartree-Fock (HF) hybrid approach. The semiempirical method MSINDO was used for the determination of the local adsorption geometry of the water molecule. With the hybrid method good agreement was obtained with the experimental band gap of 7.4 eV determined with electron energy loss spectroscopy of polycrystalline MgSO(4).H(2)O samples under ultrahigh vacuum conditions. The valence bands of the (100) surfaces of both MgSO(4) and MgSO(4).H(2)O are formed mainly by the O2p levels, whereas the S2p states contribute to the lower part of the conduction band. The preferred adsorption site of water at MgSO(4).H(2)O (100) is above a surface Mg atom. The water molecule is stabilized by two additional hydrogen bonds with surface atoms. Only small differences between the electronic structure of MgSO(4).H(2)O and MgSO(4) were observed. Also, the molecular adsorption of water on the MgSO(4).H(2)O surface leads to only small shifts of the electronic energy levels.