A study on orientation and absorption spectrum of interfacial molecules by using continuum model

In this work, a numerical procedure based on the continuum model is developed and applied to the solvation energy for ground state and the spectral shift against the position and the orientation of the interfacial molecule. The interface is described as a sharp boundary separating two bulk media. The polarizable continuum model (PCM) allows us to account for both electrostatic and nonelectrostatic solute-solvent interactions when we calculate the solvation energy. In this work we extend PCM to the interfacial system and the information about the position and orientation of the interfacial molecule can be obtained. Based on the developed expression of the electrostatic free energy of a nonequilibrium state, the numerical procedure has been implemented and used to deal with a series of test molecules. The time-dependent density functional theory (TDDFT) associated with PCM is used for the electron structure and the spectroscopy calculations of the test molecules in homogeneous solvents. With the charge distribution of the ground and excited states, the position- and orientation-dependencies of the solvation energy and the spectrum have been investigated for the interfacial systems, taking the electrostatic interaction, the cavitation energy, and the dispersion-repulsion interaction into account. The cavitation energy is paid particular attention, since the interface portion cut off by the occupation of the interfacial molecule contributes an extra part to the stabilization for the interfacial system. The embedding depth, the favorable orientational angle, and the spectral shift for the interfacial molecule have been investigated in detail. From the solvation energy calculations, an explanation has been given on why the interfacial molecule, even if symmetrical in structure, tends to take a tilting manner, rather than perpendicular to the interface.     

A molecular-dynamics study of structural and physical properties of nitromethane nanoparticles

The structural and physical properties of nanoparticles of nitromethane are studied by using molecular dynamics methods with a previously developed force field. [Agrawal et al., J. Chem. Phys. 119, 9617 (2003).] This force field accurately predicts solid- and liquid-state properties as well as melting of bulk nitromethane. Molecular dynamics simulations of nanoparticles with 480, 240, 144, 96, 48, and 32 nitromethane molecules have been carried out at various temperatures. The carbon-carbon radial distribution function, dipole-dipole correlation function, core density, internal enthalpy, and atomic diffusion coefficients of the nanoparticles were calculated at each temperature. These properties were used to characterize the physical phases and thus determine the melting transitions of the nanoparticles. The melting temperatures predicted by the various properties are consistent with one another and show that the melting temperature increases with particle size, approaching the bulk limit for the largest particle. A size dependence of melting points has been observed in experimental and theoretical studies of atomic nanoparticles, and this is a further demonstration of the effect for large nanoparticles of complex molecular materials.     

X-ray absorption spectroscopy study of the hydrogen bond network in the bulk water of aqueous solutions

We utilized X-ray absorption spectroscopy (XAS) and X-ray Raman scattering (XRS) in order to study the ion solvation effect on the bulk hydrogen bonding structure of water. The fine structures in the X-ray absorption spectra are sensitive to the local environment of the probed water molecule related to the hydrogen bond length and angles. By varying the concentration of ions, we can distinguish between contributions from water in the bulk and in the first solvation sphere. We show that the hydrogen bond network in bulk water, in terms of forming and breaking hydrogen bonds as detected by XAS/XRS, remains unchanged, and only the water molecules in the close vicinity to the ions are affected.     

Structure and dynamics at the surface of a concentrated aqueous solution of CsF

A molecular dynamics simulation of a liquid layer of a concentrated CsF solution in water has been performed in order to compare the results with those obtained in an experimental study of our group. The main result of the experiment was the existence of a monolayer of nearly pure water constituting the surface and a homogeneous mixture constituting the bulk of the system. The simulation reveals the same phenomena which can be explained by the circumstance that the ions near the surface mostly keep their first solvation shell intact. The water molecules belonging to these shells and being placed on the vapor side constitute this monolayer. The density profiles of the ions indicate that the Cs ions penetrate further into the surface than the F ions. The orientational structure of the first shell of water molecules around an ion is the same for ions in the surface and ions in the bulk in contrast to the dynamics which is altered. The spectra of the librational motion are shifted to lower frequencies. In addition to that the spectra belonging to libration which involves motion of the dipole moment develop a peak in the low frequency range irrespective of whether the water molecules are bonded to Cs or to F ions. This can be correlated with an overall preferred orientation of the water molecules in the surface which is most pronounced for the dipole moment. The calculation of the diffusion coefficients shows that the top surface layer of nearly pure water is a region of enhanced and extremely anisotropic mobility. The mean residence time of water molecules in the surface in the first shell of an ion is reduced according to the enhanced mobility.     

Communication: On the locality of hydrogen bond networks at hydrophobic interfaces

The formation of structured hydrogen bond networks in the solvation shells immediate to hydrophobic solutes is crucial for a large number of water mediated processes. A long lasting debate in this context regards the mutual influence of the hydrophobic solute into the bulk water and the role of the hydrogen bond network of the bulk in supporting the solvation structure around a hydrophobic molecule. In this context we present a molecular dynamics study of the solvation of various hydrophobic molecules where the effect of different regions around the solvent can be analyzed by employing an adaptive resolution method, which can systematically separate local and nonlocal factors in the structure of water around a hydrophobic molecule. A number of hydrophobic solutes of different sizes and two different model potential interactions between the water and the solute are investigated.     

Hybrid supermolecule-polarizable continuum approach to solvation: Application to the mechanism of the Stevens rearrangement

Semiempirical molecular orbital theory has been used to study the effects of solvation by acetonitrile on the Stevens rearrangement of methylammonium formylmethylide to 2-aminopropanal. Three methods of solvation have been used to investigate both the electrostatic and specific solvent–solute effects of solvation: a supermolecule calculation involving the complete geometry optimization of up to six solvent molecules about the solute, the conductor-like screening model (COSMO) polarizable continuum method which allows for geometry optimization of the solute in a solvent defined by its dielectric constant, and a hybrid method in which up to five solvent molecules are incorporated inside the solute cavity and complete geometry optimization of the complex is carried out within the polarizable continuum. A comparison of the calculated geometries, rearrangement activation energies, and enthalpies of solvation from these approaches is presented, and the explicit versus bulk solvation effects are discussed. The overall effect of all methods for incorporating solvation effects is that the radical pair pathway is perferred over the concerted mechanism. © 1996 by John Wiley & Sons, Inc.     

Molecular dynamics simulations of melting and the glass transition of nitromethane

Molecular dynamics simulations have been used to investigate the thermodynamic melting point of the crystalline nitromethane, the melting mechanism of superheated crystalline nitromethane, and the physical properties of crystalline and glassy nitromethane. The maximum superheating and glass transition temperatures of nitromethane are calculated to be 316 and 160 K, respectively, for heating and cooling rates of 8.9 x 10(9) Ks. Using the hysteresis method [Luo et al., J. Chem. Phys. 120, 11640 (2004)] and by taking the glass transition temperature as the supercooling temperature, we calculate a value of 251.1 K for the thermodynamic melting point, which is in excellent agreement with the two-phase result [Agrawal et al., J. Chem. Phys. 119, 9617 (2003)] of 255.5 K and measured value of 244.73 K. In the melting process, the nitromethane molecules begin to rotate about their lattice positions in the crystal, followed by translational freedom of the molecules. A nucleation mechanism for the melting is illustrated by the distribution of the local translational order parameter. The critical values of the Lindemann index for the C and N atoms immediately prior to melting (the Lindemann criterion) are found to be around 0.155 at 1 atm. The intramolecular motions and molecular structure of nitromethane undergo no abrupt changes upon melting, indicating that the intramolecular degrees of freedom have little effect on the melting. The thermal expansion coefficient and bulk modulus are predicted to be about two or three times larger in crystalline nitromethane than in glassy nitromethane. The vibrational density of states is almost identical in both phases.     

Crystal growth of para-sexiphenyl on clean and oxygen reconstructed Cu(110) surfaces

The formation of crystalline para-sexiphenyl (6P) films on Cu(110) and Cu(110)-(2 × 1)O (Cu-O) has been studied by low energy electron diffraction, X-ray absorption spectroscopy and both in situ and ex situ X-ray diffraction methods to elucidate the transition from the initial monolayers to crystalline thin films. It is found that, for Cu-O, a single and, for Cu(110), a double wetting layer is formed which then acts as a template for the subsequent 3D crystal growth. For both substrates the orientation of the long molecular axes of the 6P molecules in the first layers is conserved for the molecules in the bulk crystals growing on them. The main difference between both systems is that on Cu-O the first monolayer assembles in a form close to that of a 6P bulk plane which can be easily continued by crystallites grown upon them, while on the Cu(110) surface the 6P mono- and bi-layers differ substantially from the bulk structure. The bi-layer forms a complex periodically striped phase. Thin 6P films grow with the 6P(203) crystal plane parallel to the Cu-O substrate surface. For this orientation, the 6P molecules are stacked in layers and the molecules demonstrate only one tilt of the mean molecular plane with respect to the sample surface. On clean Cu(110), a more complex 6P(629) plane is parallel to the substrate surface and this orientation is likely a consequence of the super-molecular long-range periodicity of the second molecular layer striped phase.     

Adaptive self-assembly: environment-induced formation and reversible switching of polynuclear metallocyclophanes

Ligand 3 has been shown to self-assemble under coordination of copper(II) cations in a 1:1 ratio in acetonitrile to give equilibrating mixtures of a [2 x 2] grid-type tetranuclear structure 1 and a hexanuclear achitecture of hexagonal shape 2. The latter was confirmed by determination of the crystal structure which further indicated that 2 contained acetonitrile molecules and hydroxo groups bound to the copper(II) centers, which are therefore five-coordinate. The structures assigned to 1 and 2 were further supported by the spectral (mass, UV/Vis) data. The self-assembly process is strongly dependent on the conditions of the medium. An increase in concentration in acetronitrile increases the relative amount of hexamer 2, which appears to be the favored entity at the highest concentrations that can be reached before precipitation occurs. On the other hand, in nitromethane only the tetranuclear complex 1 was detected by mass spectrometry. Replacement of nitromethane by acetonitrile and vice versa indicated the reversible switching between a solution containing either 1 alone or an equilibrium mixture of 1 and 2, respectively. In conclusion, the system described presents several remarkable features: 1) self-assembly with substrate binding, 2) dynamic combinatorial structure generation, and 3) environment-induced structural switching amounting in effect to a process of adaptive self-assembly.     

Structural and thermodynamic aspects of ionic solvation in concentrated aqueous poly(ethylene glycol)

Solvation of ions in concentrated aqueous poly(ethylene glycol) (PEG) has been studied from thermodynamic and structural viewpoints using ion-transfer voltammetry at the interface between aqueous and nitrobenzene phases and X-ray absorption fine structure (XAFS). Systematic changes in the ion-transfer potential from water to aqueous PEG have been confirmed for several ions relative to the corresponding potential of tetraethylammonium ion (Et4N+), which is almost independent of PEG concentration. The results obtained for alkali cations strongly suggest the involvement of their complexation with PEG even in relatively diluted PEG solutions. It has been implied that the solvation circumstances of Br- and ClO4- are drastically altered when the PEG concentration becomes higher than particular critical values (e.g., 30-50% PEG200), where free water molecules are diminished because of the hydration of PEG. XAFS measurements have also been performed for K+ and Br- to get direct evidence for these findings. Although the spectra at the K K-edge clearly indicate the presence of a PEG complex of K+ in relatively diluted PEG solutions ( approximately 33% PEG200), an obvious increase in its ion-transfer potential has been detected at lower PEG concentrations, indicating that complexes formed at the interface rather than in bulk solution are transferred into an organic phase. Br- is fully hydrated in 0-50% PEG solutions, whereas some water molecules are replaced by PEG when the PEG concentration increases. Increasing the PEG concentration causes decreases in the coordination number from 6 in water to 2-3 in neat PEG. Thus, the present approach not only has elucidated the structural and thermodynamic aspects of ionic solvation in aqueous PEG but also has provided the information of the hydration of PEG.     

Interfacial water on hydrophobic surfaces recognized by ions and molecules

Recent spectrophotometric and molecular dynamics simulation studies have shown that the physicochemical properties and structures of water in the vicinity of hydrophobic surfaces differ from those of the bulk water. However, the interfacial water acting as a separation medium on hydrophobic surfaces has never been detected and quantified experimentally. In this study, we show that small inorganic ions and organic molecules differentiate the interfacial water formed on the surfaces of octadecyl-bonded (C(18)) silica particles from the bulk water and the chemical separation of these solutes in aqueous media with hydrophobic materials can be interpreted with a consistent mechanism, partition between the bulk water phase and the interfacial water formed on the hydrophobic surface. Thermal transition behaviour of the interfacial water incorporated in the nanopores of the C(18) silica materials and the solubility parameter of the water calculated from the distribution coefficients of organic compounds have indicated that the interfacial water may have a structure of disrupted hydrogen bonding. The thickness of the interfacial water or the limit of distance from the hydrophobic surface at which molecules and ions can sense the surface was estimated to be 1.25 ± 0.13 nm from the volume of the interfacial water obtained by a liquid chromatographic method and the surface area, suggesting that the hydrophobic effect may extend beyond the first solvation shell of water molecules directly surrounding the surfaces.     

Structural and dynamic properties of water within the solvation layer around various conformations of the glycine-based polypeptide

Several conformations of the solvated glycine-based polypeptides were investigated using molecular dynamics simulations. Some properties of water in the neighboring space around these molecules were investigated. It was found that water forms a well-defined layer-the first solvation shell-around the peptide molecule, and thickness of this layer is independent of the peptide structure and is equal to approximately 0.28 nm. Within this layer, water molecules show marked orientations relative to a peptide surface. Using the two-particle contribution to entropy as a measure of structural ordering of water, we found that the first solvation shell contributes 95% or more to the total water ordering around the peptide molecule. In investigating the dynamic properties of water, diffusion coefficients and lifetime of the hydrogen bond, clear differences between solvation layer and the bulk water were observed. It was found that the translational diffusion coefficient, D(T), decreases by 30% or more compared to bulk water; also, the lifetime of the water-water hydrogen bond clearly increases. The rotational diffusion coefficient, however, decreases only slightly, no more than approximately 10%. These differences correspond to the slightly higher energy of the hydrogen bond, and to its slightly distorted geometry. Analyzing the translational dynamics of water in the vicinity of the peptide molecule, it was deduced that the structure of the first solvation shell becomes more rigid than the structure of the bulk water. Investigation of a "pure hydrophobic" form of the polypeptide shows that the structure and the properties of water within the solvation shell are predominantly determined by the hydrophobic effect. The specific interactions between water molecules and various charge groups of the peptide molecule modifies this effect only slightly.     

Probing the transition from hydrophilic to hydrophobic solvation with atomic scale resolution

Picosecond and femtosecond X-ray absorption spectroscopy is used to probe the changes of the solvent shell structure upon electron abstraction of aqueous iodide using an ultrashort laser pulse. The transient L(1,3) edge EXAFS at 50 ps time delay points to the formation of an expanded water cavity around the iodine atom, in good agreement with classical and quantum mechanical/molecular mechanics (QM/MM) molecular dynamics (MD) simulations. These also show that while the hydrogen atoms pointed toward iodide, they predominantly point toward the bulk solvent in the case of iodine, suggesting a hydrophobic behavior. This is further confirmed by quantum chemical (QC) calculations of I(-)/I(0)(H(2)O)(n=1-4) clusters. The L(1) edge sub-picosecond spectra point to the existence of a transient species that is not present at 50 ps. The QC calculations and the QM/MM MD simulations identify this transient species as an I(0)(OH(2)) complex inside the cavity. The simulations show that upon electron abstraction most of the water molecules move away from iodine, while one comes closer to form the complex that lives for 3-4 ps. This time is governed by the reorganization of the main solvation shell, basically the time it takes for the water molecules to reform an H-bond network. Only then is the interaction with the solvation shell strong enough to pull the water molecule of the complex toward the bulk solvent. Overall, much of the behavior at early times is determined by the reorientational dynamics of water molecules and the formation of a complete network of hydrogen bonded molecules in the first solvation shell.     

The importance of the shape of the protein-water interface of a kinesin motor domain for dynamics of the surface atoms of the protein

A single kinesin motor domain immersed in water has been investigated using molecular dynamics. It has been found that local properties of water in the solvation shell change along with the nature of the neighboring protein surface. However, a detailed analysis leads to the conclusion that the geometrical features of hydrogen bonds and overall structure of kinesin hydration water are not very different from bulk water. The local values of diffusion coefficients (translational and rotational) of water adjacent to specific patches on the protein surface seem not to be correlated to the orientational ordering of hydration water, but instead they depend on spatial roughness and degree of exposure of the patch to the solvent. Finally, a relationship between the mobility of various surface atoms of the protein and the mean values of the diffusion coefficient of the adjacent water molecules has been observed. The latter finding suggests a close relationship between the dynamics of the inner kinesin movements and the behavior of solvation water which is in turn determined by the topography of the contact surface between the protein and the surrounding water molecules.     

有机改性剂对氯代酚毛细管区带电泳分离的影响

在分离19种氯代酚的过程中,考察了不同的有机添加剂对其毛细管区带电泳分离的影响,发现除了缓冲溶液的pH值外,缓冲溶液添加剂对氯代酚的电泳分离也有较大影响。这种影响与添加剂和氯代酚形成氢键的能力有关。     

Correlation of the vibrations of the aqueous azide ion with the O-H modes of bound water molecules

Dual frequency two-dimensional infrared spectroscopy (2D-IR) has been used to investigate the dynamics of the azide-water solvation shell. The memory of the azide transition frequencies is detected in the echo emitted by the OH stretching mode of the ion-bound water molecules. There is a significant positive correlation of the two frequency distributions that decays on a 140 fs time scale. The result confirms that the O-H bond of water molecules in the solvent shell have frequency fluctuations that are considerably slowed from those that are known in bulk water. The positive correlation is attributed to cooperative interactions of coordinated water molecules with an azide ion.     

Coupling of complex aromatic ring vibrations to solvent through hydrogen bonds: effect of varied on-ring and off-ring hydrogen-bonding substitutions

In this study, we examine the coupling of a complex ring vibration to solvent through hydrogen-bonding interactions. We compare phenylalanine, tyrosine, l-dopa, dopamine, norepinephrine, epinephrine, and hydroxyl-dl-dopa, a group of physiologically important small molecules that vary by single differences in H-bonding substitution. By examination of the temperature dependence of infrared absorptions of these molecules, we show that complex, many-atom vibrations can be coupled to solvent through hydrogen bonds and that the extent of that coupling is dependent on the degree of both on- and off-ring H-bonding substitution. The coupling is seen as a temperature-dependent frequency shift in infrared spectra, but the determination of the physical origin of that shift is based on additional data from temperature-dependent optical experiments and ab initio calculations. The optical experiments show that these small molecules are most sensitive to their immediate H-bonding environment rather than to bulk solvent properties. Ab initio calculations demonstrate H-bond-mediated vibrational coupling for the system of interest and also show that the overall small molecule solvent dependence is determined by a complex interplay of specific interactions and bulk solvation characteristics. Our findings indicate that a full understanding of biomolecule vibrational properties must include consideration of explicit hydrogen-bonding interactions with the surrounding microenvironment.     

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.     

Hydration and dewetting near graphite-CH(3) and graphite-COOH plates

The dynamics of water near the nanoscale hydrophobic (graphite-CH(3)) and hydrophilic (graphite-COOH) plates has been studied in detail with molecular dynamics simulations in this paper. It is shown that these designed surfaces (by growing a layer of methyl or carboxyl groups on top of graphite) can have a significant impact on the neighboring water dynamics, with the hydrophilic carboxyl surface having even more profound effects. The water hydrogen bond lifetime is much longer near both types of surfaces than that in the bulk, while on the other hand the water diffusion constant is much smaller than that in the bulk. The difference in the diffusion constant can be as large as a factor of 8 and the difference in the hydrogen bond lifetime can be as large as a factor of 2, depending on the distance from the surface. Furthermore, the water molecules in the first solvation shell of surface atoms show a strong bias in hydroxyl group orientation near the surface, confirming some of the previous findings. Finally, the possible water dewetting transition between two graphite-CH(3) plates and the effect of the strength of the solute-solvent attractions on the water drying transition are investigated. The relationship among the dewetting transition critical distance, van der Waals potential well depth, and water contact angle on the graphite-CH(3) surface is also analyzed on the basis of a simple macroscopic theory, which can be used to predict the dewetting transition critical distance.     

Strength of hydrogen bonds of water depends on local environment

In-depth knowledge of water-water potential is important for devising and evaluating simple water models if they are to accurately describe water properties and reflect various solvation phenomena. Water-water potential depends upon inter-molecular distance, relative orientation of water molecules, and also local environment. When placed at a favorable distance in a favorable orientation, water molecules exhibit a particularly strong attractive interaction called hydrogen bond. Although hydrogen bond is very important for its effects on the elements of life, industrial applications, and bulk water properties, there is no scientific consensus on its true nature and origin. Using quantum-mechanical methods, hydrogen bond strength was calculated in different local environments. A simple empirical linear relationship was discovered between maximum hydrogen bond strength and the number of water molecules in the local environment. The local environment effect was shown to be considerable even on the second coordination shell. Additionally, a negative linear correlation was found between maximum hydrogen bond strength and the distance, at which it was observed. These results provide novel insights into the nature of hydrogen bonding.