Probing the nucleation mechanism for the binary n-nonane/1-alcohol series with atomistic simulations

The AVUS-HR approach, which combines histogram reweighting with aggregation-volume-bias Monte Carlo nucleation simulations using self-adaptive umbrella sampling, was extended to multicomponent nucleation systems. It was applied to investigate the homogeneous vapor-liquid nucleation for the binary n-nonane/1-alcohol series, including the n-nonane/methanol, n-nonane/ethanol, n-nonane/1-propanol, n-nonane/1-butanol, n-nonane/1-hexanol, and n-nonane/1-decanol systems. The simple transferable potentials for phase equilibria-united atom force field was used in this investigation. It was found that the nucleation free energy (NFE) contour plots obtained for these binary n-nonane/1-alcohol nucleation systems exhibit rather interesting mechanistic features, some of which are distinct from other binary systems previously studied (such as water/ethanol and water/n-nonane). In addition, the NFE profiles show a subtle evolution with the increase in alcohol chain length, from a somewhat two-pathway type of shape as observed for the n-nonane/methanol system to a more normal single-pathway one for systems involving longer alcohols (1-hexanol and 1-decanol). In contrast, the NFE maps obtained for the other three binary systems involving those medium-length alcohols display the most striking feature with the saddle point stretched almost all the way from the n-nonane-enriched to the alcohol-enriched domain, implying that multiple pathways coexist on the nucleation map. These free energy profiles were shown to be consistent with the non-ideal nucleation behavior observed experimentally for this binary series, namely, a rather reluctant conucleation of the alcohols with n-nonane. In particular, this non-ideal behavior becomes more severe with a decrease in the alcohol chain length. Also, analysis of the compositions of the critical nuclei indicates a reluctant mixing behavior between these two species, i.e., depletion of the alcohol at low alcohol activity or depletion of n-nonane at low n-nonane activity, in agreement with the experimental interpretations. Furthermore, a microscopic inhomogeneity is present inside these critical nuclei, that is, alcohols aggregate via hydrogen bonds forming alcohol-enriched domains.     

Molecular content and structure of aqueous organic nanodroplets from the vapor-liquid nucleation study of the water/n-nonane/1-alcohol series

Prompted by a previous finding of unusual mixing behavior for the critical clusters involved in the vapor-liquid nucleation of the ternary water/n-nonane/1-butanol mixture, atomistic simulations employing the AVUS-HR technique were carried out to extend such investigations to include both shorter and longer alcohols, namely, the water/n-nonane/CiH2i+1OH mixture with i = 2, 4, 6, and 8. It is clear from this extensive investigation that the miscibility between water and n-nonane can be further improved by increasing the chain length of the alcohol (surfactant). In fact, for the water/n-nonane/1-octanol mixture at an intermediate gas-phase activity composition, the nucleation can proceed via fully mixed critical nuclei containing a roughly equal amount of all three components, which is in contrast to the dominantly binary-like nucleation channels observed for such mixtures involving shorter alcohols. Structural analysis revealed that these mixed nuclei take on a multilayered structural motif of the core-shell (water-alcohol) type with n-nonane distributed outside, forming an additional layer, more or less uniformly, compared to the one-sided deposition found for systems involving shorter alcohols. This structure provides a microscopic origin for the enhanced miscibility of water with n-nonane observed in the presence of 1-alcohol. These results may also have important implications for atmospheric organic aerosols.     

Towards understanding the nucleation mechanism for multi-component systems: an atomistic simulation of the ternary nucleation of water/n-nonane/1-butanol

Inspired by the previous finding of some unusual vapour/liquid nucleation results on the ternary water/n-nonane/1-butanol system, atomistic simulations were carried out for a detailed investigation of this mixture. These simulations reproduced the experimentally-reported non-ideal nucleation behaviour for this system, including both onset activities and the average compositions of the critical nuclei. Close examination of the nucleation free energy data and the structure of the critical nuclei reveals two types of phase separation. One occurs internally inside the cluster via formation of a multi-layered structure. The other takes place externally, leading to the coexistence of multiple nucleation channels, characterized by critical clusters of different compositions. Such mechanistic and structural heterogeneity is the microscopic origin of the complex nucleation behaviour observed for this ternary mixture.     

Dumbbells and onions in ternary nucleation

Molecular simulations for a ternary nucleation system (water/n-nonane/1-butanol) demonstrate a more complex nucleation mechanism than previously thought, where critical nuclei with different compositions are present even for a given vapour-phase composition; the spatial distribution in these critical nuclei is heterogeneous and dumbbell and onion motifs are found; in the former, water and nonane nano-droplets are connected through a butanol handle, while in the latter a water core is surrounded by a nonane corona with an interfacial butanol shell.     

Dynamic structure of methane/n-nonane clusters during nucleation and growth

We report results on nucleation, growth, and structure formation of methane/n-nonane clusters in an expanding system investigated by molecular dynamics simulation. From bulk phase equilibria data, it is expected that the concentration of the less volatile substance n-nonane in the clusters is very high. However, analyses of experimental data in the literature suggest somewhat higher methane content at onset of nucleation. Our simulations show that the methane mole fraction is actually very high and increases even further at the beginning of the cluster growth. On the other hand, in this transient state after nucleation the methane mole fraction in the cluster core decreases, leaving a n-nonane rich core, i.e., we observe the phase separation inside the growing cluster. Methane is squeezed out from the core to the surface and then evaporates from the surface shell during expansion of the system.     

Experimental investigation of the Soret effect in acetone/water and dimethylsulfoxide/water mixtures

The thermal diffusion behavior of acetone/water and dimethylsulfoxide(DMSO)/water mixtures has been experimentally investigated by a transient holographic grating technique named thermal diffusion forced Rayleigh scattering (TDFRS). For both systems a sign change of the Soret coefficient S(T) with varying water content has been predicted by simulations [C. Nieto Draghi et al., J. Chem. Phys. 122, 114503 (2005)]. The sign change of S(T) is confirmed by the experiment. Except for equimolar concentrations of acetone/water the agreement between the experimental and simulation data is reasonable.     

Evaporation of alcohol/water mixtures through hydrophobic porous membranes

The evaporation of methanol/water, ethanol/water and propanol/water mixtures across hydrophobic porous membranes (Gore-Tex polytetrafluoroethylene membranes) with different nominal pore radii was studied at ambient temperatures. The evaporation process leads to a higher concentration of the alcohols in the vapor phase. There exists only a limited composition range (dilute solutions of alcohols in water) in which evaporation experiments are possible. The experimental results are interpreted in terms of the average pore diameter of the membranes, the vaporliquid equilibrium data of the systems and the measured contact angles between methanol/water, ethanol/water and n-propanol/water mixtures and polytetrafluoroethylene, respectively.     

Mixtures of sodium dodecyl sulfate/dodecanol at the air/water interface by computer simulations

Molecular dynamics simulations of monolayers of surfactant mixtures at the air/water interface were performed where the binary mixture was composed of sodium dodecyl sulfate (SDS) and dodecanol molecules. At the same ratio of SDS and dodecanol molecules, two monolayer mixtures were prepared. In the first monolayer, all the dodecanol molecules were placed together in the center of the simulation box, whereas in the second monolayer, those molecules were uniformly distributed in the surface area in such a way that they were far from each other. Simulations of both systems indicate that the dodecanol tails in the first monolayer are straighter and more ordered than those in the second monolayer. From the present results, we observed new insights of how the different molecules should array or distribute at the interface in real systems. Finally, studies of the interfacial water around the different surfactants were also analyzed, showing that they are closer to the polar headgroups of dodecanol than to the SDS headgroups.     

Analyzing the molecular basis of enzyme stability in ethanol/water mixtures using molecular dynamics simulations

One of the drawbacks of nonaqueous enzymology is the fact that enzymes tend to be less stable in organic solvents than in water. There are, however, some enzymes that display very high stabilities in nonaqueous media. In order to take full advantage of the use of nonaqueous solvents in enzyme catalysis, it is essential to elucidate the molecular basis of enzyme stability in these media. Toward this end, we performed μs-long molecular dynamics simulations using two homologous proteases, pseudolysin, and thermolysin, which are known to have considerably different stabilities in solutions containing ethanol. The analysis of the simulations indicates that pseudolysin is more stable than thermolysin in ethanol/water mixtures and that the disulfide bridge between C30 and C58 is important for the stability of the former enzyme, which is consistent with previous experimental observations. Our results indicate that thermolysin has a higher tendency to interact with ethanol molecules (especially through van der Waals contacts) than pseudolysin, which can lead to the disruption of intraprotein hydrophobic interactions and ultimately result in protein unfolding. In the absence of the C30-C58 disulfide bridge, pseudolysin undergoes larger conformational changes, becoming more open and more permeable to ethanol molecules which accumulate in its interior and form hydrophobic interactions with the enzyme, destroying its structure. Our observations are not only in good agreement with several previous experimental findings on the stability of the enzymes studied in ethanol/water mixtures but also give an insight on the molecular determinants of this stability. Our findings may, therefore, be useful in the rational development of enzymes with increased stability in these media.     

甘氨酸在DMF/水和乙醇/水混合溶剂中的焓对相互作用

利用LKB2277生物活性检测仪对298.15K时甘;氨酸在纯水,DBF/水和乙醇/水混合溶剂DBF(或乙醇的质量分数为0.05~0.45,以10%递增)中的稀释焓进行了测定,并利用维里展开式法关联得到焓对相互作用系数值.结果表明,焓对相互作用系数h~2与溶剂的结构性质密切相关.在DBF/水混合物溶剂中,h~2随DBF含量的增加而变得更负;而在乙醇/水混合溶剂中,h~2与乙醇含量近似呈抛物线型关系,在乙醇质量分数为0.25时h~2达到一个极小值。     

Solvation phenomena of a tetrapeptide in water/trifluoroethanol and water/ethanol mixtures: a diffusion NMR,intermolecular NOE,and molecular dynamics study

Solvation of a tetrapeptide, NAc-Ser-Phe-Val-Gly-OMe (1), in water and in water/alcohol mixtures with 2,2,2-trifluoroethanol (TFE)/water or ethanol (ETH)/water has been studied by diffusion NMR and intermolecular NOE measurements. The experimental results were compared with those obtained from detailed Molecular Dynamics (MD) calculations. Independently, all three methods revealed preferential solvation on the surface of the peptide by TFE in the water/TFE mixtures, but not by ETH in the water/ETH mixtures. The MD calculations show that the TFE concentration coating the peptide is higher than that in the bulk, while for ethanol, the concentration is nearly equal to that in the bulk. Calculated site-specific preferential solvation data between TFE, ETH, and water with the different peptide groups have been compared with the NMR data and shown to be in general agreement with the experimental facts.     

Systematic study of the thermal diffusion in associated mixtures

We performed systematic temperature and concentration dependent measurements of the Soret coefficient in different associated binary mixtures of water, deuterated water, dimethyl sulfoxide (DMSO), methanol, ethanol, acetone, methanol, 1-propanol, 2-propanol, and propionaldehyde using the so-called thermal diffusion forced Rayleigh scattering method. For some of the associating binary mixtures such as ethanol/water, acetone/water, and DMSO/water, the concentration xw+/- at which the Soret coefficient changes its sign does not depend on temperature and is equal to the concentration xw x where the Soret coefficient isotherms intersect. For others such as 1-propanol/water, 2-propanol/water, and ethanol/DMSO, the sign change concentration is temperature dependent, which is the typical behavior observed for nonassociating mixtures. For systems with xw+/-=xw x, we found that xw+/- depends linearly on the ratio of the vaporization enthalpies of the pure components. Probably due to the similarity of methanol and DMSO, we do not observe a sign change for this mixture. The obtained results are related to structural changes in the fluid observed by nuclear magnetic resonance, mass spectrometric, and x-ray experiments in the literature. Furthermore, we discuss the influence of hydrophilic and hydrophobic interactions and the solubility on thermal diffusion behavior.     

Molecular dynamics simulation of quasi-two-dimensional water clusters on ice nucleation protein

We performed molecular dynamics simulations of systems that consisted of the ice nucleation protein and the quasi-two-dimensional water cluster on it. The angle distributions, percolation probabilities, mean cluster sizes, cluster size distributions, and hydrogen bond relaxation times were analyzed. We concluded that the behavior of the water clusters on the ice nucleation protein was elaborately intertwined by the interaction between the ice nucleation protein and water, the interaction between the water molecules and the effect of temperature. The percolation probability and mean cluster size depended on the interactions and temperatures.     

Selective adsorption from methanol/water mixtures by C60 fullerene nanospheres

Micro-Raman spectroscopy was used to investigate the selective adsorption of methanol/water mixtures on the surface of [60] fullerene nanospheres. C₆₀ molecules were dispersed in methanol/water mixtures with different methanol molar fractions ranging between 1 and 0.5. The Raman active pentagon pinch mode shifted significantly (±4 cm⁻¹) as the mixture composition was changed. The shift in the Raman mode was sinusoidal in nature indicating that methanol then water is adsorbed preferentially on the fullerene surface at different mixture compositions. The observed behavior is attributed to structure forming effects of alcohol/water mixtures and the shape and size effect of fullerene surface.     

Solvent structures of mixed water/acetonitrile mixtures at chromatographic interfaces from computer simulations

Atomistic simulations are used to characterize the molecular dynamics (MD) of alkyl chains with different functionalizations in different water/acetonitrile mixtures (80/20 and 50/50). Starting from fully equilibrated solvent systems (flat density profile for both components), microheterogeneous structuring of the solvent in the chromatographic system is found for both mixtures. Depending on the functionalization of the alkyl chain (nitrile, amide, nitro, phenyl), differences in the density profiles of the two solvents (water/acetonitrile), the effective width of the stationary phase and the solvent gradients in the overlap region are observed. The solvent mixture (mobile phase) in RPLC is a liquid which is directly involved in the physical process and must be included explicitly. Far from the surface, the solvent displays bulk properties; closer to it the mixed solvent partitions due to the presence of the stationary phase. This creates a gradient in solvent strength perpendicular to the surface which influences the motions of the analyte. The surface is found to define the amount of water that can bind to it and defines its hydrophilic character. Proposals from the literature, such as the existence of persistent water filaments extending from the functionalized silica layer towards the bulk solvent, are discussed. Simulations of acridine orange near a -NH(2)- and -phenol-functionalized surface highlight the different dynamical behaviour (insertion vs. adsorption) of an analyte depending on the functionalization of the surface.     

Preferential solvation of phenol in binary solvent mixtures. A molecular dynamics study

Molecular dynamics computer simulations were carried out to study the preferential solvation of phenol in equimolar acetonitrile-water and ethanol-water binary mixtures. Two water models were used to investigate the model dependence of preferential solvation. The results are compared to recent intermolecular 1H NOESY experiments reported on the same systems. In the case of acetonitrile-water the local mole fraction obtained from simulations agrees quite well with experiments. In the case of ethanol-water there was a qualitative difference, which was observed for both water models. However, when comparing the degree of preferential solvation of the two cosolvents ethanol and acetonitrile with each of the two water models, the trend obtained from the simulations agrees with experimental data.     

Binary nucleation rates for ethanol/water mixtures in supersonic Laval nozzles

Although the conditions corresponding to the onset of condensation of aqueous-alcohol mixtures have been measured in supersonic nozzles [B. E. Wyslouzil et al., J. Chem. Phys. 113, 7317 (2000)], the true nucleation rates have not. Here, we propose a new analytical method to estimate the temperature, the concentrations of condensable species in both the vapor and the liquid phases, and the amount of the condensate using only the measured static pressure profiles in the nozzle. We applied the method to ethanol/water (CH(3)CH(2)OH/D(2)O or CH(3)CH(2)OD/D(2)O) mixtures and confirmed that the aerosol volume fractions derived from pressure measurements and small angle neutron scattering measurements are in very good agreement when this method is used. Combining the results from the pressure measurements with the number densities of the condensed droplets, measured either by small angle neutron or small angle x-ray scattering, we determined the first quantitative ethanol/water binary nucleation rates in the supersonic nozzle at a temperature of 229±1?K.     

Ethanol, acetic acid, and water adsorption from binary and ternary liquid mixtures on high-silica zeolites

Adsorption isotherms were measured for ethanol, acetic acid, and water adsorbed on high-silica ZSM-5 zeolite powder from binary and ternary liquid mixtures at room temperature. Ethanol and water adsorption on two high-silica ZSM-5 zeolites with different aluminum contents and a high-silica beta zeolite were also compared. The amounts adsorbed were measured using a recently developed technique that accurately measures the changes in adsorbent/liquid mixture density and liquid concentration. This technique allows the adsorption of each compound in a liquid mixture to be measured. Adsorption data for binary mixtures were fit with the dual-site extended Langmuir model, and the parameters were used to predict ternary adsorption isotherms for each compound with reasonable accuracy. In ternary mixtures, acetic acid competed with ethanol and water for adsorption sites and reduced ethanol adsorption more than it reduced water adsorption.