Nuclear Magnetic Resonance Analysis of Hydrothermal Reactions of Ethyl- and Octylamine in Sub- and Supercritical Water.

Hydrothermal reactions of aliphatic amines have recently gained importance in relation to the application of amines as film-forming corrosion inhibitors for steam-water cycles. The kinetics and mechanism of the hydrothermal reactions of ethylammomiun cation (EtAH+) and n-octylammonium cation (OctAH+) were studied for comparison with the corresponding neutral amines to elucidate their reaction products and pathways at sub- and supercritical temperatures of 300-400 °C as model reactions of aliphatic ammonium cations. We analyzed the reaction of 13C-15N-labeled EtAH+ using NMR spectroscopy and revealed that the initial hydrolysis to ethanol, known as the main path, is followed by the elimination reaction producing ethene and the disproportionation reaction giving diethylammonium cation. The OctAH+ yields octene and octanol, each of which isomerizes to thermodynamically more stable species as the major products. Comparisons were made between the reactions of the neutral amines and ammonium cations to highlight their different reactivity. The hydrolysis, alkene formation, and dehydration of alcohols to alkenes were all found to be accelerated at low pH. The formation of low-molecular-weight organic acids such as acetic acid and formic acid was not observed. These results indicate that the corrosion protection effect of film-forming amines will be maintained under practical conditions with pH values as high as around 9 to 10, and hence side reactions involving byproducts will be suppressed.

Self-diffusion of water-cyclohexane mixtures in supercritical conditions as studied by NMR and molecular dynamics simulation.

The self-diffusion coefficients of water (Dw) and cyclohexane (Dch) in their binary mixtures were determined using the proton pulsed field gradient spin-echo method from medium to low densities in subcritical and supercritical conditions. The density (ρ), temperature (T), and water mole fraction (xw) are studied in the ranges 0.62-6.35 M (M = mol dm-3), 250-400 °C, and 0.109-0.994, respectively. A polynomial fitting function was developed for a scaled value of Ξ = ρDT-1/2 with ρ, T, and xw as variables in combination with a comprehensive molecular dynamics (MD) simulation. The NMR and MD results agree within 5% for water and 6% for cyclohexane, on average. The differences between Dw and Dch in the dependence on ρ, T, and xw are characterized by the activation energy Ea and the activation volume ΔVΞ ‡ expressed by the scaled fitting function. The decrease in the ratio Dw/Dch and the increase in the Ea of water with increasing xw are related to the increase in the number of hydrogen bonds (HBs). The Dw value for a solitary water molecule at a low xw is controlled by the solvation shell, most of which is occupied by nonpolar cyclohexane molecules that provide less friction as a result of weaker interactions with water. A microscopic diffusion mechanism is discussed based on an analysis of the HB number as well as the first-peak height of the radial distribution functions that are taken as measures of the potential of the mean field controlling self-diffusion.

Effect of water on hydrolytic cleavage of non-terminal α-glycosidic bonds in cyclodextrins to generate monosaccharides and their derivatives in a dimethyl sulfoxide-water mixture.

Hydrolytic cleavage of the non-terminal α-1,4-glycosidic bonds in α-, ß-, and γ-cyclodextrins and the anomeric-terminal one in d-maltose was investigated to examine how the cleavage rate for α-, ß-, and γ-cyclodextrins is slower than that for d-maltose. Effects of water and temperature were studied by applying in situ (13)C NMR spectroscopy and using a dimethyl sulfoxide (DMSO)-water mixture over a wide range of water mole fraction, xw = 0.004-1, at temperatures of 120-180 °C. The cleavage rate constant for the non-anomeric glycosidic bond was smaller by a factor of 6-10 than that of the anomeric-terminal one. The glycosidic-bond cleavage is significantly accelerated through the keto-enol tautomerization of the anomeric-terminal d-glucose unit into the d-fructose one. The smaller the size of the cyclodextrin, the easier the bond cleavage due to the ring strain. The remarkable enhancement in the cleavage rate with decreasing water content was observed for the cyclodextrins and d-maltose as well as d-cellobiose. This shows the important effect of the solitary water whose hydrogen bonding to other water molecules is prohibited by the presence of the organic dipolar aprotic solvent, DMSO, and which has more naked partial charges and higher reactivity. A high 5-hydroxymethyl-2-furaldehyde (5-HMF) yield of 64% was attained in a non-catalytic conversion by tuning the water content to xw = 0.30, at which the undesired polymerization by-paths can be most effectively suppressed. This study provides a step toward designing a new optimal, earth-benign generation process of 5-HMF starting from biomass.

Effect of water content on conversion of D-cellobiose into 5-hydroxymethyl-2-furaldehyde in a dimethyl sulfoxide-water mixture.

Noncatalytic conversion of D-cellobiose (at 0.5 M) into 5-hydroxymethyl-2-furaldehyde (5-HMF), a platform chemical for fuels and synthetic materials, was analyzed at 120-200 °C over a wide range of water mole fraction, xw = 0.007-1 in a binary dimethyl sulfoxide (DMSO)-water mixture by means of the in situ (13)C NMR spectroscopy. Effects of the water content were revealed as follows: (i) The tautomerization of the anomeric residue of D-cellobiose from D-glucose to D-fructose type was not initially observed at a lower water content, in contrast to the significant tautomerization into the D-fructose type in a higher water content and pure water. (ii) The lower the water content, the faster the glycosidic-bond cleavage by hydrolysis, because of the high reactivity of solitary water molecules with the large partial charges more naked as in supercritical water clusters due to the isolation by the organic solvent DMSO. (iii) The amount of D-fructose as the intermediate product was larger at the higher xw; despite the increase of D-fructose, the production of 5-HMF from D-fructose was slowed down. (iv) A high 5-HMF yield of 71% was reached at xw = 0.20-0.30 that was 6-10 times the initial D-cellobiose concentration. The best yield of 5-HMF was attained in the low xw region when the polymerization paths into NMR-undetectable species via 5-HMF and anhydromonosaccharides were effectively suppressed. This study provides a new framework to design optimal, noncatalytic reaction process to produce 5-HMF from cellulosic biomass by tuning the water content as well as the temperature and the reaction time.

Effect of heavy hydrogen isotopes on the vibrational line shape for supercritical water through rotational couplings.

The rotational couplings, which determine the infrared spectral line shape in the low-density supercritical water, were analyzed as functions of the density and the temperature by applying molecular dynamics simulation to a flexible water model, SPC∕Fw and by varying the moment of inertia of the water through substitution for the H atom in H2O by heavy hydrogen isotopes. The differences in the frequency and the relative intensity between the sharp center peak and the rotational broad side-bands were analyzed for the O-H, O-D, and O-T stretch spectra. The frequency differences between the sharp center peak and the rotational broad side-bands are linearly correlated with the inverse of the moment of inertia of the isotope-substituted water species. The intensity of the sharp peak is associated with the long-time component of the reorientational time correlation function for the stretching bond vector. At 400 °C, where a substantial amount of hydrogen bonds are dynamically persisting, an intensity decrease in the rotational broad side-bands was observed with increasing density from 0.01 to 0.40 g cm(-3), respectively, corresponding to 0.56 and 22.2 M (=mol dm(-3)), orders of magnitude higher than the ideal gas densities. This arises from the decrease in the correlation time of the angular velocity and the rotational couplings due to an increase in the hydrogen-bonding perturbation. The intensity decrease of the rotational side-bands with increasing density is more significant for the water isotopes with heavier hydrogens. At a high temperature of 1200 °C, the rotational side-bands at 0.01 to 0.05 g cm(-3) were more distinct than those at 400 °C, and even at a medium density of 0.40 g cm(-3) a significant signal broadening due to the rotational couplings was clearly observed because of the accelerated rotational momentum. The vibrational spectrum cannot be decomposed into definite chemical clusters for the thermodynamic and kinetic analysis because of the dynamic origin.

Solvent effect on pathways and mechanisms for D-fructose conversion to 5-hydroxymethyl-2-furaldehyde: in situ 13C NMR study.

Noncatalytic reactions of D-fructose were kinetically investigated in dimethylsulfoxide (DMSO), water, and methanol as a function of time at temperatures of 30-150 °C by applying in situ (13)C NMR spectroscopy. The products were quantitatively analyzed with distinction of isomeric species by taking advantage of site-selective (13)C labeling technique. In DMSO, D-fructose was converted first into 3,4-dihydroxy-2-dihydroxymethyl-5-hydroxymethyltetrahydrofuran having no double bond in the ring, subsequently into 4-hydroxy-5-hydroxymethyl-4,5-dihydrofuran-2-carbaldehyde having one double bond through dehydration, and finally into 5-hydroxymethyl-2-furaldehyde (5-HMF) having two double bonds. No other reaction pathways were involved, as shown from the carbon mass balance. In water, 5-HMF, the final product in DMSO, was generated with the precursors undetected and furthermore transformed predominantly into formic and levulinic acids and slightly into 1,2,4-benzenetriol accompanied by polymerization. D-glucose was also produced through the reversible transformation of the reactant D-fructose. In methanol, some kinds of anhydro-D-fructoses were generated instead of 5-HMF. The reaction pathways can thus be controlled by taking advantage of the solvent effect. The D-fructose conversion reactions are of the first order with respect to the concentration of D-fructose and proceed on the order of minutes in DMSO but on the order of hours in water and methanol. The rate constant was three orders of magnitude larger in DMSO than in water or methanol.

Nuclear magnetic resonance study on rotational dynamics of water and benzene in a series of ionic liquids: anion and cation effects.

The rotational correlation times (τ(2R)) for polar water (D(2)O) molecule and apolar benzene (C(6)D(6)) molecule were determined in ionic liquids (ILs) by means of the (2)H (D) NMR spin-lattice relaxation time (T(1)) measurements. The solvent IL was systematically varied to elucidate the anion and cation effects separately. Five species, bis(trifluoromethylsulfonyl)imide (TFSI(-)), trifluoromethylsulfonate (TfO(-)), hexafluorophosphate (PF(6)(-)), chloride (Cl(-)), and formate (HCOO(-)), were examined for the anion effect against a fixed cation species of 1-butyl-3-methyl-imidazolium (bmim(+)). Four species, bmim(+), N-methyl-N-butylpyrrolidinium (bmpy(+)), N,N,N-trimethyl-N-propylammonium (N(1,1,1,3)(+)), and P,P,P-trihexyl-P-tetradecylphosphonium (P(6,6,6,14)(+)), were employed for the cation effect against a fixed anion species of TFSI(-). The τ(2R) ratio of water to benzene, expressed as τ(W/B), was used as a probe to characterize the strength of Coulombic solute-solvent interaction in ILs beyond the hydrodynamic limit based on the excluded-volume effect. The τ(W/B) value was found to strongly depend on the anion species, and the solute dynamics are sensitive not only to the size but also to the chemical structure of the component anion. The cation effect was rather weak, in contrast. The largest and most hydrophobic P(6,6,6,14)(+) cation was exceptional and a large τ(W/B) was observed, indicating a unique solvation structure in [P(6,6,6,14)(+)]-based ILs.

Density effect on infrared spectrum for supercritical water in the low- and medium-density region studied by molecular dynamics simulation.

The origin of the line shape of the O-H stretch vibrational spectrum is analyzed for supercritical water in the low- and medium-density region by using classical molecular dynamics simulation for the flexible point-charge model, SPC/Fw. The spectrum calculated for the water model is in good agreement with the experimental one in the low-density region. The spectral origins in the low-density region of 0.01-0.04 g cm(-3) are assigned to a sharp peak due to the bond oscillation along the O-H vector and two broad bands due to the rotational coupling, by taking an isolated single molecule as a reference in the low-density limit. The bands due to the rotational coupling reduce in intensity with increasing density as the rotations are more hindered by the hydrogen-bonding interactions, and their intensities increase with increasing temperature due to the accelerated rotational motion. The O-H stretch oscillation in the time correlation function attenuates in a timescale comparable with the lifetime of the hydrogen bonds, and the spectra conditioned by the number of hydrogen bonds are dominantly controlled by the local solvation structure.

Noncatalytic hydrothermal elimination of the terminal D-glucose unit from malto- and cello-oligosaccharides through transformation to D-fructose.

Noncatalytic hydrothermolyses of malto- and cello-oligosaccharides (di-, tri-, tetraose), linked by α- and ß-1,4-glycosidic bonds, respectively, were investigated at 100-140 °C. In situ (13)C NMR spectroscopy was applied to elucidate the position and pathways of the glycosidic bond breakage and the dependence of the hydrolysis rate on the bond type. Spectral analysis was carried out quantitatively as a function of time with the mass balance confirmed, and it was shown for both the malto- and the cello-oligosaccharides that the terminal D-glucose unit with a free anomeric carbon is selectively eliminated after transformation to D-fructose. Site-selective breakage of the glycosidic bonds proceeded on the order of hours. The initial apparent rates for terminal hydrolysis were found to be independent of the degree of oligomerization but dependent on the type of glycosidic bond. Rate constants were larger for the α-1,4-linked malto-oligosaccharides by a factor of 3-4 than for the ß-1,4-linked cello-oligosaccharides. The pathways and mechanisms for the malto- and cello-oligosaccharide hydrothermolyses are common and can be understood in terms of the elementary reactions of the di- and monosaccharides.

Investigation of in situ oxalate formation from 2,3-pyrazinedicarboxylate under hydrothermal conditions using nuclear magnetic resonance spectroscopy.

We have investigated the assembly of a two-dimensional coordination polymer, Nd(2)(C(6)H(2)N(2)O(4))(2)(C(2)O(4))(H(2)O)(2), that has been prepared from the hydrothermal reaction of Nd(NO(3))(3)·6H(2)O and 2,3-pyrazinedicarboxylic acid (H(2)pzdc). In situ oxalate formation as observed in this system has been been investigated using (1)H and (13)C nuclear magnetic resonance spectroscopy, and a pathway for C(2)O(4)(2-) anion formation under hydrothermal conditions has been elucidated. The oxalate ligands found in Nd(2)(C(6)H(2)N(2)O(4))(2)(C(2)O(4))(H(2)O)(2) result from the oxidation of H(2)pzdc, which proceeds through intermediates, such as 2-pyrazinecarboxylic acid (2-pzca), 2-hydroxyacetamide, 3-amino-2-hydroxy-3-oxopropanoic acid, 2-hydroxymalonic acid, 2-oxoacetic acid (glyoxylic acid), and glycolic acid. The species are generated through a ring-opening that occurs via cleavage of the C-N bond of the pyrazine ring, followed by hydrolysis/oxidation of the resulting species.

Rotational dynamics of benzene and water in an ionic liquid explored via molecular dynamics simulations and NMR T1 measurements.

The rotational dynamics of benzene and water in the ionic liquid (IL) 1-butyl-3-methylimidazolium chloride are studied using molecular dynamics (MD) simulation and NMR T(1) measurements. MD trajectories based on an effective potential are used to calculate the (2)H NMR relaxation time, T(1) via Fourier transform of the relevant rotational time correlation function, C(2R)(t). To compensate for the lack of polarization in the standard fixed-charge modeling of the IL, an effective ionic charge, which is smaller than the elementary charge is employed. The simulation results are in closest agreement with NMR experiments with respect to the temperature and Larmor frequency dependencies of T(1) when an effective charge of ±0.5e is used for the anion and the cation, respectively. The computed C(2R)(t) of both solutes shows a bi-modal nature, comprised of an initial non-diffusive ps relaxation plus a long-time ns tail extending to the diffusive regime. Due to the latter component, the solute dynamics is not under the motional narrowing condition with respect to the prevalent Larmor frequency. It is shown that the diffusive tail of the C(2R)(t) is most important to understand frequency and temperature dependencies of T(1) in ILs. On the other hand, the effect of the initial ps relaxation is an increase of T(1) by a constant factor. This is equivalent to an "effective" reduction of the quadrupolar coupling constant (QCC). Thus, in the NMR T(1) analysis, the rotational time correlation function can be modeled analytically in the form of aexp (-t/τ) (Lipari-Szabo model), where the constant a, the Lipari-Szabo factor, contains the integrated contribution of the short-time relaxation and τ represents the relaxation time of the exponential (diffusive) tail. The Debye model is a special case of the Lipari-Szabo model with a = 1, and turns out to be inappropriate to represent benzene and water dynamics in ILs since a is as small as 0.1. The use of the Debye model would result in an underestimation of the QCC by a factor of 2-3 as a compensation for the neglect of the Lipari-Szabo factor.

In situ kinetic study on hydrothermal transformation of D-glucose into 5-hydroxymethylfurfural through D-fructose with 13C NMR.

Kinetics of hydrothermal reaction of D-glucose was investigated at 0.02 M over a temperature range of 120-160 °C by applying in situ (13)C NMR spectroscopy. D-Glucose was found to be reversibly transformed first into D-fructose (intermediate) and successively into 5-hydroxymethylfurfural (5-HMF) through dehydration. The carbon mass balance has been kept within the detection limit, and no other reaction pathways are present. The hydrothermal reaction of d-glucose is thus understood as that of D-fructose in the sense that the D-glucose reaction proceeds only through D-fructose. All the isomers of D-glucose and D-fructose were detected by the in situ (13)C NMR in D(2)O: they are the open chains and the pyranoses and furanoses of α- and ß-types. The ß-forms are the most stable due to the hydration. For both D-glucose and D-fructose, the isomers are in a rapid equilibrium for each monosaccharide, and they are treated collectively in the kinetic analysis of the slower hydrothermal reactions. The reactions are of the first order with respect to the concentrations of D-glucose and D-fructose, and D-glucose converts to 5-HMF on the order of hours. The kinetic parameters were determined by the in situ method.

NMR-NOE and MD simulation study on phospholipid membranes: dependence on membrane diameter and multiple time scale dynamics.

Motional correlation times between the hydrophilic and hydrophobic terminal groups in lipid membranes are studied over a wide range of curvatures using the solution-state (1)H NMR-nuclear Overhauser effect (NOE) and molecular dynamics (MD) simulation. To enable (1)H NMR-NOE measurements for large vesicles, the transient NOE method is combined with the spin-echo method, and is successfully applied to a micelle of 1-palmitoyl-lysophosphatidylcholine (PaLPC) with diameter of 5 nm and to vesicles of dipalmitoylphosphatidylcholine (DPPC) with diameters ranging from 30 to 800 nm. It is found that the NOE intensity increases with the diameter up to ∼100 nm, and the model membrane is considered planar on the molecular level beyond ∼100 nm. While the NOE between the hydrophilic terminal and hydrophobic terminal methyl groups is absent for the micelle, its intensity is comparable to that for the neighboring group for vesicles with larger diameters. The origin of NOE signals between distant sites is analyzed by MD simulations of PaLPC micelles and DPPC planar bilayers. The slow relaxation is shown to yield an observable NOE signal even for the hydrophilic and hydrophobic terminal sites. Since the information on distance and dynamics cannot be separated in the experimental NOE alone, the correlation time in large vesicles is determined by combining the experimental NOE intensity and MD-based distance distribution. For large vesicles, the correlation time is found to vary by 2 orders of magnitude over the proton sites. This study shows that NOE provides dynamic information on large vesicles when combined with MD, which provides structural information.

Communication: exploring the reorientation of benzene in an ionic liquid via molecular dynamics: effect of temperature and solvent effective charge on the slow dynamics.

The rotational time correlation function (RTCF) of solute benzene molecules in the ionic liquid (1-butyl-3-methylimidazolium chloride) has been studied using classical molecular dynamics simulation. The effect of solvent charge on the functional form of RTCF was investigated by comparing four force fields for the solvent where the total charge on the anion and the cation was set to ±1e, ±0.7e, ±0.5e, and 0, respectively. For all three charged solvent models, the RTCF exhibits a long-time tail where the relaxation rate exhibits a significant slowdown. This feature is strengthened by higher solvent charges as well as lower temperatures, indicating the influence of the strong Coulombic fields arising from the solvent charges. The long-time tail is caused by the extraordinarily slow solvent structural relaxation of ionic liquids compared to the time scale of their local vibrational and librational dynamics.

Controlling the equilibrium of formic acid with hydrogen and carbon dioxide using ionic liquid.

The equilibrium for the reversible decomposition of formic acid into carbon dioxide and hydrogen is studied in the ionic liquid (IL) 1,3-dipropyl-2-methylimidazolium formate. The equilibrium is strongly favored to the formic acid side because of the strong solvation of formic acid in the IL through the strong Coulombic solute-solvent interactions. The comparison of the equilibrium constants in the IL and water has shown that the pressures required to transform hydrogen and carbon dioxide into formic acid can be reduced by a factor of approximately 100 by using the IL instead of water. The hydrogen transformation in such mild conditions can be a chemical basis for the hydrogen storage and transportation using formic acid.

Water as an in situ NMR indicator for impurity acids in ionic liquids.

A sensitive in situ NMR spectroscopic method for detecting acids contaminating ionic liquids (ILs) has been developed. The chemical shift and the spectral width of water added to ILs were used as indicators to measure the impurity acid level. Owing to the high resolution power of NMR, the detection limit is below the level of 10(-3) mol kg(-1). A new method is applicable to a number of commonly used ILs such as the imidazolium- and ammonium-based ILs except for those composed of acidic cations or anions. The method was utilized to monitor the purification efficiency in the recrystallization of a typical hydrophilic IL, 1-butyl-3-methylimidazolium methanesulfonate from acetone. It was demonstrated that impurity acids can be almost perfectly removed by single or double recrystallization.

Self-diffusion coefficients for water and organic solvents at high temperatures along the coexistence curve.

The self-diffusion coefficients D for water, benzene, and cyclohexane are determined by using the pulsed-field-gradient spin echo method in high-temperature conditions along the liquid branch of the coexistence curve: 30-350 degrees C (1.0-0.58 g cm(-3)), 30-250 degrees C (0.87-0.56 g cm(-3)), and 30-250 degrees C (0.77-0.48 g cm(-3)) for water, benzene, and cyclohexane, respectively. The temperature and density effects are separated and their origins are discussed by examining the diffusion data over a wide range of thermodynamic states. The temperature dependence of the self-diffusion coefficient for water is larger than that for organic solvents due to the large contribution of the attractive hydrogen-bonding interaction in water. The density dependence is larger for organic solvents than for water. The difference is explained in terms of the van der Waals picture that the structure of nonpolar organic solvents is determined by the packing effect due to the repulsion or exclusion volumes. The dynamic solvation shell scheme [K. Yoshida et al., J. Chem. Phys. 127, 174509 (2007)] is applied for the molecular interpretation of the translational dynamics with the aid of molecular dynamics simulation. In water at high temperatures, the velocity relaxation is not completed before the relaxation of the solvation shell (mobile-shell type) as a result of the breakdown of the hydrogen-bonding network. In contrast, the velocity relaxation of benzene is rather confined within the solvation shell (in-shell type).

Hydrothermal C-C bond formation and disproportionation of acetaldehyde with formic acid.

Reaction pathways and kinetics of C2 (carbon-two) aldehyde, acetaldehyde (CH3CHO), and formic acid HCOOH or HOCHO, are studied in neutral and acidic subcritical water at 200-250 degrees C. Acetaldehyde is found to exhibit (i) the acid-catalyzed C-C bond formation between acetaldehyde and formic acid, which generates lactic acid (CH3CH(OH)COOH), (ii) the cross-disproportionation, where formic acid reduces acetaldehyde into ethanol, and (iii) the aldol condensation. The lactic acid formation is a green C-C bond formation, proceeding without any organic solvents or metal catalysts. The new C-C bond formation takes place between formic acid and aldehydes irrespective of the presence of alpha-hydrogens. The hydrothermal cross-disproportionation produces ethanol without base catalysts and proceeds even in acidic condition, in sharp contrast to the classical base-catalyzed Cannizzaro reaction. The rate constants of the reactions (i)-(iii) and the equilibrium constant of the lactic acid formation are determined in the temperature range of 200-250 degrees C and at HCl concentrations of 0.2-0.6 M (mol/dm(3)). The reaction pathways are controlled so that the lactic acid or ethanol yield may be maximized by tuning the reactant concentrations and the temperature. A high lactic acid yield of 68% is achieved when acetaldehyde and formic acid are mixed in hot water, respectively, at 0.01 and 2.0 M in the presence of 0.6 M HCl at 225 degrees C. The ethanol yield attained 75% by the disproportionation of acetaldehyde (0.3 M) and formic acid (2.0 M) at 225 degrees C in the absence of added HCl.

Free-energy analysis of the molecular binding into lipid membrane with the method of energy representation.

A statistical-mechanical treatment of the molecular binding into lipid membrane is presented in combination with molecular simulation. The membrane solution is viewed as an inhomogeneous, mixed solvent system, and the free energy of solvation of a solute in membrane is computed with a realistic set of potential functions by the method of energy representation. Carbon monoxide, carbon dioxide, benzene, and ethylbenzene are adopted as model solutes to analyze the binding into 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC) membrane. It is shown that the membrane inside is more favorable than bulk water and that the solute distribution is diffuse throughout the membrane inside. The membrane-water partition coefficient is then constructed with the help of the Kirkwood-Buff theory from the solvation free energy obtained separately in the hydrophobic, glycerol, headgroup, and aqueous regions. To discuss the role of repulsive and attractive interactions, the solvation free energy is partitioned into the DMPC and water contributions and the effect of water to stabilize the benzene and ethylbenzene solutes within the membrane is pointed out.

Partial pair correlation functions of low-density supercritical water determined by neutron diffraction with the H/D isotopic substitution method.

Neutron diffraction measurements were carried out on H/ D isotopically substituted water in the low-density supercritical condition (T = 673 K, P = 26.3 MPa, and rho = 0.0068 molecules.A-3) in order to obtain direct information on the intermolecular partial structure functions, gHH inter(r), gOH inter(r), and gOO inter(r). In correspondence to the high-density supercritical water previously reported, the intermolecular nearest neighbor peaks in gHH inter(r), gOH inter(r), and gOO inter(r) are diffuse compared with the ambient ones. The nearest neighbor O...O distance (3.3 A) and its coordination number (2.6) were determined from the observed gOO inter(r). These results indicate that the orientational correlation between neighboring water molecules is considerably lost in low-density supercritical water. Small clusters involving a few water molecules are preferentially formed in low-density supercritical water, which is consistent with results obtained by previous IR and NMR studies.