alexa Molecular Dynamics Simulations of Atomistic Detailed Hydration Structures of Poly(vinyl methyl ether) | Open Access Journals
ISSN: 2161-0398
Journal of Physical Chemistry & Biophysics
Like us on:
Make the best use of Scientific Research and information from our 700+ peer reviewed, Open Access Journals that operates with the help of 50,000+ Editorial Board Members and esteemed reviewers and 1000+ Scientific associations in Medical, Clinical, Pharmaceutical, Engineering, Technology and Management Fields.
Meet Inspiring Speakers and Experts at our 3000+ Global Conferenceseries Events with over 600+ Conferences, 1200+ Symposiums and 1200+ Workshops on
Medical, Pharma, Engineering, Science, Technology and Business

Molecular Dynamics Simulations of Atomistic Detailed Hydration Structures of Poly(vinyl methyl ether)

Rongliang Wu1,2*, Xinlong Qiu1, Xiguo Zeng2, Bin Kong2 and Xiaozhen Yang2
1College of Material Science and Engineering, State Key Laboratory for Modification of Chemical Fibers and Polymer Materials, Donghua University, Shanghai, China
2Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Polymer Science and Engineering, Joint Laboratory of Polymer Science and Materials, Beijing, China
Corresponding Author : Rongliang Wu
College of Material Science and Engineering
State Key Laboratory for Modification of Chemical Fibers
and Polymer Materials, Donghua University
Shanghai 201620, China
Tel: +862167792914
Received: October 27, 2015 Accepted: November 07, 2015 Published: November 10, 2015
Citation: Wu R, Qiu X, Zeng X, Kong B, Yang X (2015) Molecular Dynamics Simulations of Atomistic Detailed Hydration Structures of Poly(vinyl methyl ether). J Phys Chem Biophys 5:192. doi:10.4172/2161-0398.1000192
Copyright: © 2015 Wu R, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Related article at Pubmed, Scholar Google

Visit for more related articles at Journal of Physical Chemistry & Biophysics


Molecular dynamics simulations have been performed on the aqueous solutions of Poly(vinyl methyl ether) (PVME) at various concentrations. Both radial and spatial distribution functions are used to investigate the detailed hydration structures. The structures of water are found to get increasingly concentrated when polymers are introduced and the water motions are severely hindered by the polymer matrix. At low concentrations, larger populations of tt conformers in meso dyads than those at higher concentrations are found and this phenomenon is believed to be due to the increase of water molecule bonding to two ether oxygens in meso dyad. At higher concentrations the sizes and conformations of polymers are quite similar with those in bulk. A transition of hydrogen bond fractions between PVME and water at around the concentration of 0.3 is observed and this value coincides perfectly well with the results of conformational analysis and Raman spectra. Second neighbor hydrogen bond statistics showed that complicated hydrogen bond networks take the lead at low concentrations, but single hydrogen bonds as well as isolated clusters composed of 2-4 water molecules are found common around each polymer repeat unit.

Molecular dynamics; PVME; Hydration structure; Solvation
Aqueous solutions of poly(vinyl methyl ether) (PVME) exhibit lower critical solution temperature (LCST) miscibility behavior with moderate critical solution temperatures around 35 ~ 37℃ [1-3]. The polymer dissolves below this temperature but phase separates when heated. In common cases, the LCST behaviors of aqueous polymer solutions can experimentally be sorted into three types [2,4]: type I represents the classical Flory-Huggins (FH) phase behavior, whose liquid-liquid phase separation critical composition shifts towards lower polymer concentrations for longer polymers; while the critical composition of type II hardly changes even for polymers with infinite molar masses; the LCST curves of type III are bimodal and have two stable critical compositions, one at lower polymer concentration (type I) and the other at higher polymer concentration (type II). The phase behavior of the PVME/H2O system belongs to type III at physiological temperatures [2].
This bimodal feature of the atypical phase behavior should, from the atomistic point of view, originates from the special hydrophilic and hydrophobic interactions between water and PVME [3,5]. The hydrophilic sites (ether oxygens) stabilize the aqueous solution by forming hydrogen bonds (HB) with neighboring water molecules while the hydrophobic groups (particularly the methyl groups due to the polarization of ether oxygens) destabilize the solution by altering the normal hydrogen bond network structure of water close to the polymer chains. The resulting structures of water around the hydrophobic groups were referred to as the “icebergs” by Nemethy and Scheraga [6], and also by Frank and Evans [7]. These two interactions compete with each other with increasing temperature until phase transition occurs [3]. Wu et al. also utilized FTIR spectroscopy to investigate the microdynamics mechanism of the PVME aqueous solution, and they found the hydrogen bonds between hydrophilic groups and water and the hydration interactions between hydrophobic groups and water are completely reversible [5].
There has been extensive investigation [1,3-5,8-16] on the nonclassical phase behavior and the above mentioned two sorts of interactions in the PVME/H2O system. From changes in the nearinfrared and calorimetric parameters with the variation of concentration, Maeda [3,10] proposed that an average of 2.7 water molecules per chain repeat unit in PVME solutions formed stable molecular complex through hydrogen bonding to the ether oxygens and additional 2.3 water molecules formed “icebergs” around the methyl groups, over 5 water molecules were “free”, namely they had nearly no interaction with the polymer. Nies et al. [1,4,14] concluded that experimental data from small angle neutron scattering did not support the existence of a stable molecular complex as the mixtures were characterized by Ornstein- Zernike-like composition fluctuations even at the temperature as low as 10℃. The bimodal LCST phase behavior of the PVME/H2O system was also predicted, by way of the Wertheim Lattice Thermodynamics Perturbation Theory (Wertheim-LTPT), together with two adjacent narrow upper critical solution temperatures (UCST) miscibility gaps at low temperature, and these results were in agreement with experimental observations. The predicted solid-liquid equilibrium of the solvent showed a pronounced drop in melting temperature at intermediate weight compositions of 40 ~ 70 wt %. This peculiar concentration dependence of the melting line of water provided a new explanation for the double melting endotherm observed in (MT)DSC rather than being linked to the melting of bound and free water, namely the formation of a stable intermolecular complex.
The phenomenological results from experiments and theoretical predictions could not provide detailed molecular insights into the observed anomalous phase behavior. Previous work in our group [12] utilized quantum mechanics (QM) calculations to investigate the interactions between the PVME model compound and clusters of 2 ~ 5 water molecules, and the detailed molecular packing structures with relatively lower energies were obtained. The most stable configurations obtained were those with single HBs for clusters with no more than 3 water molecules and those with double HBs for clusters with more than 3 water molecules. The HBs between water and the ether oxygens were proved to be stronger than those between waters, namely these HBs are shorter in length and lower in energy. The hydrophobic interactions, called “QHBs” (quasi-hydrogen bonds), were found to coexist and associate with the HBs in many of those configurations. Different hydration features and mechanisms were analyzed and an uneven distribution of water adsorption on the polymer was found to be preferred. An average of 1.5 water molecules per repeat unit was suggested if all units of the polymer were bound with water.
Despite that the detailed configurations of possible packing structures can be successfully derived, QM calculations are still unable to describe the influence of polymer chain connectivity and conformation on the hydration structures. With the fast development of computational capabilities, atomistic molecular dynamics (MD) simulations have been extensively used to investigate the aqueous solutions of organic compounds, such as alcohols [17], Urea [18], amides [19] and dimethylsulfoxide [20]. Tamai et al. [8] utilized the AMBER/OPLS force field [21] and SPC/E water model [22] to simulate the hydrogel model structures of PVA (poly vinyl alcohol), PVME and PNiPAM (poly N-isopropylacrylamide). It was found in the hydrophilic regions the average number of hydrogen bonds between water molecules was smaller by 1.0 ~ 1.9 than those in pure water, and the defects of water-water hydrogen bonds were completely compensated by the water-polymer HBs for PVA but not for PVME and PNiPAM; In the hydrophobic regions, the average number of HB between water molecules was only slightly smaller than those in pure water, the water-water HBs were stabilized in the hydrophobic regions especially for PVME and PNiPAM, and this was deemed to be owing to the hydrophobic hydration effects. It was also found that in the regions out of the first hydration shell, the HB numbers were not at all influenced by the polymer matrix. In their simulations, the stability of HBs was much more enhanced in the systems with lower water contents than those with more water, and this phenomenon was attributed to the constraints on the mutual orientation of water and polar groups imposed by the hydrophobic side group at large polymer concentrations. While in their simulations, the concentration dependence of the hydration structures of PVME was not fully described.
In the present work, MD simulations have been performed on the aqueous solutions of PVME at concentrations varying from 0 to 100 wt %. The OPLS-AA force field [23,24] and the TIP4P water model [25] are used, for TIP4P is the only one with a realistic dimmer structure and its quadrupole is the most realistic among the most frequently used water models [26-28]. In addition, the most recent OPLS-AA force field parameters have been optimized with the TIP4P water model for better description of the interactions between water and organic substances [27].
Molecular models and force field
Classical molecular dynamics simulations of aqueous solutions of PVME were carried out using the GROMACS simulation package [29] with the OPLS-AA force field and the TIP4P water model. The potential functions adopt the form in which the former two terms represent the harmonic vibrations of bond lengths and bond angles, respectively; The third term stands for the dihedral potentials, whose OPLS parameters are translated into the Ryckaert-Bellemans [30] type in GROMACS; The latter two are non-bonded interactions including the Lennard-Jones (L-J 12-6) and the electrostatic interactions. Table 1 lists the diagonal non-bonded parameters for PVME and the TIP4P water, in which ε and σ are the well depth and collision diameter for the L-J 12-6 interaction, q is the polar charge, and m is the atomic mass. The standard combination rules [24] are used for the off diagonal non-bonded atomic pairs
image (2)
where image. Intra-molecular interactions between all pairs of atoms separated by more than three bonds are regarded as ordinary non-bonded interactions, while the 1-4 interactions are monitored by both the dihedral potentials and the ordinary non-bonded interactions scaled by 0.5 [24].
Simulation details
All PVME chains used in the current study contain 44 backbone carbons and all repeat units are connected head to tail. The chains are atactic and were generated according to the Monte Carlo method described in the appendix with steric ratio meso:racemic (m:r)=1:1. Table 2 lists the number of PVME chains np, number of water molecules nw, the number of water molecules each polymer repeat unit can ideally be assigned nw/nmon and the weight concentration cp of each system. The initial amorphous structures were built using the Amorphous Builder in Accelrys Materials Studio 3.1, 20 structures were built for each system and the ones with the lowest potential energies were chosen for subsequent procedures. These chosen structures were then “solvated” in water using the “genbox” command in GROMACS and each system underwent hundreds of steps of energy minimization with the conjugate gradient (CG) method.
The constant particle, constant pressure and constant temperature (NPT) MD runs were then carried out for 4-10 ns depending on the concentrations of the solution systems, and the equilibrium was monitored by the potential energies, densities, hydrogen bond numbers and polymer conformations. The canonical ensemble (NVT) sampling MD simulations followed and extended for another 20 ns. The temperatures of all solution systems were maintained at 300 K with the Nose-Hoover [31,32] temperature coupling method, and the pressures were isotropically coupled to 1 bar with the Parrinello- Rahman [33] pressure coupling method. The bulk PVME system was firstly equilibrated at 600 K then annealed to 300 K and another 10 ns further equilibration was performed before NVT runs. A cutoff of 1.2 nm was used for the Van der Waals interactions. The Particle Mesh Ewald (PME) method [34,35] was used for long range electrostatic corrections and the long range dispersion corrections were also implemented for both energy and pressure. All bonds in polymers were constraint with the Lincs method [36] and the rigid structure of water molecules were maintained with the Settle method [37] during all MD simulations. A time step of 1 fs was used and the trajectories were recorded every 1 ps.
Radial and spatial distribution functions
The structure of liquid can often be well characterized by the radial distribution function (RDF) of atomic density,
image (3)
which represents the probability of finding particle B within the range r+dr around particle A. The spatial distribution function (SDF) takes into account both radial and angular correlations between particles, and allows three dimensional (3D) visualization of the local atomic densities, thus it provides much more detailed information on the local structural environment than RDF does. The visualization of SDF is often described by iso-surfaces, which show the distribution of particles with densities several times that of average around a reference atomic group. To calculate SDF, a local coordinate system has to be defined with a reference structure. For water, the local coordinate is defined so that oxygen stays at the origin and the whole molecule lies in the XY plane. For PVME repeat unit, the ether oxygen is set as the origin and the plane formed by the oxygen and the two connected carbons is regarded as the reference, so that the backbone configurations and conformations do not influence this local coordinate system.
Definitions for HB and QHB
The formation of a hydrogen bond needs an electronegative proton donor (D) and an acceptor (A) with at least one lone pairs of electrons. The atomic pairs in our simulations are considered to be hydrogen bonded if the distance between D and A is within 0.35 nm, and simultaneously the angle between the DA axis and one of the DH bonds is smaller than 30°. This geometrical criteria is the same as that used by Luzar and Chandler [38,39]. For QHBs, there exist various definitions in literatures and our definition refers to that of Li et al. [19]. Based on the Van der Waals diameter [40] and quantum mechanics calculations [12], the QHBs are defined with the distance r(C-O)<0.372 nm and the angle C-H…O>150°.
Diffusion coefficient of water
The diffusion of water can be estimated from the Mean Square Displacement (MSD) within a time period of t
image (5)
The experimental [41] diffusion coefficient of pure water is (2.11- 2.66) × 10-5 cm2/s at 298 K. The calculated value at 300 K is 3.78 × 10-5 cm2/s, which is slightly larger than the experimental result but still corresponds well with the simulation values (3.73 × 10-5 cm2/s at 298 K) using the same water model.
Measurements of the PVME Raman spectra
Poly (vinyl methyl ether) with Mw=60.5 kg/mol and Mw/Mn=3 was synthesized in the group of Prof. Nies in Universiteit of Leuven. Aqueous solutions with concentrations cp=0.1, 0.3, 0.6 and 0.8 as well as the bulk were prepared for use. The Raman spectra were recorded at room temperature using a JY-HR800 (Jobin Yvon) Raman micro spectrometer. A Nd:YAG laser (532 nm and 6 mW) was used as the excitation source. The Isotropic Raman scattering intensities were calculated from the parallel and perpendicular components of the experimentally measured Raman scattered light: image
image (5)
The 314 cm-1 band of carbon tetrachloride was used to verify the accuracy of the scattering geometry and light collection optics used, and gave the expected value of 0.75 within 2%.
Results and Discussion
Densities and excess volumes of mixing
From NPT runs, the obtained equilibrium density of pure water (0.993 ± 0.006 g/cm3) corresponds well with experiments (~ 0.997 g/ cm3). While the density of bulk PVME at 300 K is 1.019 ± 0.004 g/cm3, slightly smaller than the experimental value of 1.061 g/cm3 [42]. This is probably caused by the relatively shorter chains simulated than used in experiments [43].
The binary mixtures exhibit strong non-ideal behavior as reflected by the excess specific volumes in Figure 1. The values are calculated from the equilibrium densities of each solution system and those of pure water and bulk PVME. All excess specific volumes are negative and their standard deviations are smaller than the symbols shown in the Figure. The largest excess specific volume is anticipated to be reached at cp=0.5 ~ 0.7, which corresponds well with the experimental values inferred by the Redlich-Kister [44] polynomial [13].
Comparing with the experimental data in Figure 1, all calculated excess specific volumes are smaller due to the conformations of the relatively short chains simulated. In aqueous polymer solutions, not only hard core volume packing and free volume effects are operational but also volume effects related to the directionality of the hydrogen bonds between water molecules and those between water molecules and PVME repeat units must be considered [13]. The conformations of the shorter chains simulated should in principle have led to the smaller excess volumes of mixing. Despite these differences, the simulated systems have in general well captured the essential features of PVME aqueous solutions.
Structure and diffusion of water in solutions
Local water structure variation with concentration: The results of RDF for atomic pairs Ow-Ow at all concentrations investigated are shown in Figure 2. All RDFs have prominent first peaks at around 0.276 nm, and nearly no position changes are observed. The RDF of Ow-Ow for pure water (not plotted for clarity) is nearly identical with that for cp=0.03. All curves show clear peaks of HB, second neighbor HB and even third peaks of longer correlations, longer correlations are diminished because of the diffusion of water molecules. The near neighbor correlations show almost similar behavior for systems below cp=0.6, while in the water-poor mixtures namely cp=0.74 and 0.86 the g(r) curves show sliding down beyond the second minimum. Together with the enhanced first and second peaks, these phenomena are deemed as the results from the association of water molecules as dimers or clusters [45,46], and such effects get stronger with the increase of concentration in these systems.
The SDFs in Figure 3 have provided more detailed and clearer pictures of the local water HB structures. The first neighbor HBs are characterized by the three density caps: two caps over the H atoms are due to proton accepting neighbors, and the single large density cap over the O atoms is due to proton donating neighbors. The density caps are plotted with the same iso-surface thresholds and the enhanced waterwater correlations with concentration are clearly observed when more polymer chains are introduced. It is obviously shown in the figure that the low concentration one has relatively thinner density caps and the high concentration one has thicker density caps. This indicates that in the high concentration, with small amount of water molecules, most are confined within the first neighbor shell of polymers, taking nearly all interaction configurations even with higher energies. And in the low concentration, large amount of water molecules has huge interaction configurations to be selected according to lower energies, avoiding the small population configurations in the first neighbor shell.
Water diffusion in solutions: The mean square displacement of water molecules are plotted in Figure 4 for all systems including pure water. At low concentrations, these water molecules show diffusion behavior, for their MSDs show linear relation with time. The diffusion coefficient characterized by the slope of these curves show monotonic decrease as the polymer concentration is increased. At high concentrations, the MSD curves even approach plateau values at longer time period, which shows that the presence of polymer matrix significantly decreases the mobility of water molecules. The cages and HBs in polymers, discussed in the previous section, might have trapped the motion of these water molecules.
The solvation structures of the Op and CH3 groups: For the investigation on the solvation structures of PVME repeat unit, the RDFs of atomic pairs Op-Ow and C3-Ow at various concentrations are plotted in Figure 5. With the increase of concentration, HBs between PVME and water also get concentrated after a concentration of 0.38. Only slight migration of the first peak is observed to increase from 0.269 nm to 0.273 nm (given in Figure 6), which are smaller than those of Ow-Ow (see Figure 2) which corresponds well with our previous quantum calculations [12]. At high concentrations, the HB having a relatively longer bond length is manifested by SDF for concentration 0.86 (Figure 3) under confinement, taking the high energy configurations. While at low concentrations with abundant water molecules, the association between hydrophilic and hydrophobic interactions was found to shorten the HBs [12]. But a small part of the ether oxygens have chances to form double hydrogen bonds, which makes the HBs longer, later statistical results of HB will demonstrate that such double HBs are scarcely seen at all concentrations.
In addition, the first minima shift toward longer distances and the first peaks get broader at high concentrations, suggesting that there are changes of HB structures. The second peaks are much broader than the first ones and at low concentrations, the heights of second peaks of Op-Ow RDFs are even larger than those of the first peaks, which might indicate that water molecules from the second hydration shell form hydrogen bonds to the water molecules already donating their protons to the ether oxygens. But this is not necessarily the case, specific structures around adjacent atoms may also cause such welldefined peaks. Compared with the slight increase of peak positions of the Op-Ow RDFs, The peaks of the C3-Ow RDFs decrease more significantly from 0.367 nm to 0.355 nm. These QHBs are much weaker than HBs, and the hydration effects will have much influence on their structures.
For a clearer view of how water molecules are arranged around the PVME repeat units, SDF corresponding to the atomic pairs referred to above are plotted in Figure 7. Only the SDF at the concentration of 0.54 is shown, for clear distinctions are not found at different concentrations since the peak height differences of g(r) at other concentrations are much smaller than those of water as seen from comparisons between Figures 2 and 5. With the polymer repeat unit as reference (C1 is plotted just for reference), the density map of water around the repeat unit is shown. The iso-surface covering the reference ether oxygen represents the strong HBs between water molecules and the ether oxygens, such iso-surface corresponds to the first peak in Figure 5 a. As for the second peaks in the Op-Ow RDFs, there exit two sorts of interactions: the second neighbor HBs and the QHBs around the methyl group. It is also evident from Figure 7 that the second neighbor HB iso-surface mixes with the QHBs, suggesting the cooperative association between HB and QHB, which results in complex water clusters around PVME. This result corresponds well with our previous quantum caculations [12].
Polymer backbone conformations: The polymers in aqueous solutions often have significant conformational changes due to the specific hydrophilic and hydrophobic hydrations. As an overall estimate of the polymer conformations, the radii of gyrations of PVME in solutions are shown in Figure 8. At low concentrations, the polymers adopt relatively extended conformations due to the specific hydrophilic and hydrophobic hydrations. With the increase of concentration, polymers aggregate and interpenetrate each other, and their sizes approach those in bulk.
We now have a look at the experimental Raman spectra at various concentrations in Figure 9. The peaks with most significant variations are denoted as a, b, c, d. Similar shapes of the peaks are observed for concentrations above 30 wt %, but for the spectra at low concentrations drastic differences are observed, the peaks of c and d level off and form plateaus. According to our previous experience [12,47,48], the peaks of a and b are related to the stretching modes of C-C and C-O; while c and d are mainly caused by the bending modes. The bending modes of C-C-C are directly related to the backbone conformations and the spectra around 400~500 cm-1 have been regarded as the disordered longitudinal acoustic mode (D-LAM) regions which well represent the variations of the long sequence of PVME conformations [47,48]. As far as the radius of gyration is concerned, the Raman spectra are in good correspondence with our MD simulations.
In order to discover the detailed conformations of PVME in solutions, Figure 10 shows the fractions of the major isomeric dihedral pairs found for both meso and racemic dyads along the chains obtained in the simulation, and the ones with very small populations are not shown for clarity. For the meso dyads, the gt conformers are preferred and the drastic drop of tt conformers around the concentration of 0.3 are the most apparent. In the bulk state, the tt conformers in meso dyads are not at all preferred, which takes only ~5% in the conformer pairs in meso dyad. But in the hydration, the significant change of the tt conformer around 0.3 demonstrates that the hydration of PVME helps formation of the tt conformer through HB interaction of at least a water molecule bonding to two ether oxygens in meso dyad as shown in Figure 11. Compared with the meso dyads, the tt conformers are preferred in the racemic dyads, and the tt conformer pair show some increase when the concentration increase above 0.3. In general, the polymer conformations are similar with that in bulk for solutions of high concentrations; but the conformations of dihedral pairs at low concentrations have more significant variations than those at high concentrations, and the conformations change at a critical concentration of 0.3. The variation of this conformer change around the polymer concentration of 0.3 shows similar behavior of the Raman spectra, in which the emergence of peaks also occurs around this concentration.
During the hydration of PVME, it has changes 1) in the radius of gyration, 2) in the conformational distribution for the tt conformers, and 3) in the Raman spectral profiles of D-LAM. Obviously, the three changes are correlated. The global size increase is due to the local conformational variation, which is predominated by the tt conformation. Geometrically, it is reasonable, because among the end-to-end distances (within four backbone bonds) of various conformation dyads, the tt conformation dyad has the longest end-to-end distance. It leads to local structure expansion and the global size increase eventually. D-LAM shows the conformational structure transition occurs accordingly.
HB and QHB statistics: For aqueous solutions of small compounds, water should be meanly distributed around the solute and the structure of the solute should not have much difference, but for polymers, the local chain conformation changes upon solvation. For PVME, each repeat unit is only possible to form at most two hydrogen bonds with water, and the hydrogen bonding in polymers should not be uniform along the chains. In order to get more quantitative picture of the solvation structures of PVME solutions, We have counted the factions of units forming i (i=0,1,2) HBs or QHBs according to the geometrical definitions defined in section 2.4.
Figure 12 a gives the fractions of Op with i (i=0,1,2) HBs with water. At low concentrations, most of the ether oxygens (over 90%) are able to form HB with water, but even at the lowest concentration investigated, there still exist ~10% ether oxygens which are not able to form efficient HBs with water due to the close contacts between local chain segments caused by the flexibility of the polymer chains. With the increase of concentration, the fraction of such HBs decrease drastically, which is, in our opinion, due greatly to the conformational change as mentioned in the previous section. The conformational changes will cause significant variations of the solvation structures as HBs or the association between HBs and QHBs. The transition at around 0.3 matches perfectly well with the results of conformational analysis and Raman spectra as previously discussed. The further decrease of HB fractions at higher concentration is conceptually attributed to the deficient number of water surrounding the polymer chains. Due to the existence of negative cooperativity effect [49] for HBs, the fractions of repeat units with double HBs are very small, but still possess a transition around the concentration of 0.3. And moreover, recent work of Nies [4] measured the demixing enthalpy of the aqueous solutions of a mono-dispersed low molecular weight (about 140 backbone carbones) PVME utilizing modulated temperature differential scanning calorimetry (MTDSC), and found there also exist a transition around the concentration of 0.3. But their results level off at low concentrations and decrease linearly over 0.3.
Figure 12 b and c depicts the fractions of Ow with i (i=0,1,2) HBs with PVME and water respectively. With the increase of concentration, the fractions of HBs with PVME increase and at high concentrations significant amount of water have been observed to be hydrogen bonded to two polymer repeat units, which might be adjacent or from two different chains. Here, the numbers of double HBs have higher fractions at high concentrations. Water can serve both as acceptors and donors, but in Figure 12 c we just count those serve as donors or acceptors and thus the water molecules can also just form at most two HBs with water, so that they can be directly compared with those HBs between polymer and water. With the increase of concentration, the fractions of water molecules that form only one HB with other water molecules increase while the fractions of those forming two HBs decrease, which indicates the severe interruption of normal water HB networks and the tendency to form isolated water clusters.
Figure 12 d describes the fractions of the CH3 group with i (i=0,1,2) QHBs with water, the fractions of QHBs are quite minor compared with those of HBs. The increased fraction of QHBs at low concentrations is strong evidence that the associations between HBs and QHBs have lead to the increase of tt conformers in the meso dyads.
Cluster analysis: The possible cluster structures with 1-5 water molecules around PVME repeat units have been surveyed with QM in our previous investigations [12]. In this study we have taken into account the polymer chain connectivity and the dynamic features of real systems, the hydration structures might not well be completely the same with the cluster structures obtained in vacuum. Based on the statistical results in Figure 12, we could make assumptions that all water molecules are within the clusters around the polymer, then at the concentration of 0.86 each repeat unit that has HB bonding with water can be assigned 1.56 water molecules, namely the clusters have on average 1.56 water molecules, this number is consistent with predictions (1.5) in our previous study [12]. With the same assumption, the clusters will have 2.77 and 5.18 water molecules at concentrations of 0.74 and 0.54 respectively. But it is apparent that not all water molecules in these solutions are in the clusters.
Through the RDFs and SDFs in previous sections, we have observed that around the polymer repeat units water molecules form at most two hydration shells. A representative snapshot of the HB bonded PVME-water cluster is shown in Figure 13 a. Most ether oxygen form first shell HB with just one water molecule, namely very few form hydrogen bond with two water molecules as is the case for the negative cooperativity effect [49]. As far as the second neighbor HB is concerned, some PVME repeat units form clusters with 2, 3 and 4 water molecules as noted by the circles, the oval-shaped circle designates the chain effect in HB between water molecules. The clusters inter-correlate with each other, forming relatively stable configurations.
We have also counted the fractions of repeat units with i (i=0-5) water molecules which are hydrogen bonded to the first shell water molecules already donating their protons to the ether oxygen. These results are added to those shown in Figure 12 a, and the resulting cluster size estimations are plotted against concentration in Figure 13 b. Here, only clusters with no more than 5 water molecules are shown, because the fractions of clusters with more water molecules are negligibly small. At low concentrations, clusters with 4 water molecules take the lead, in such structures all hydrogens and lone pairs of the first shell water molecules are hydrogen bonded and the entire structures form complicated 3D HB networks. The number of such structures decreases as the concentration increase due to the change of polymer conformations as demonstrated in section 3.4, and they do not take overwhelming proportions at high concentrations. The fractions of clusters with 5 water molecules are barely perceivable at low concentrations and look very much similar with those of double HBs in Figure 12 a, since only double HBs can form more than 3 second neighbor HBs. With increase of concentration, fractions of 3 water clusters show slight decrease while those of 2 water clusters obtained slight increase, at high concentrations the regional clusters composed of 2-4 water molecules are commonly seen, but in these systems most of the repeat units are wrapped in polymers due to interchain packing. What is worth mentioning that isolated HBs between water and ether oxygen are seen to take ~ 3% at the concentration of 0.86.
Through atomistic molecular dynamics simulations of the PVME solutions with various concentrations at 300K, structure variations of water, polymer and their correlations were analyzed. The water structures get increasingly concentrated when polymers are introduced, and the water motions were severely hindered by the polymer matrix. The Op-Ow RDFs show two well defined peaks, the first peaks represent the HBs between the polymer and water and they have shorter lengths than those between water molecules; the second peaks are much broader and are found to be contributed from both the second neighbor HBs and the QHBs around the CH3 groups.
At low concentrations, the polymers adopt relatively extended conformations and the tt conformers in meso dyads undergo drastic decrease at the concentration of around 0.3, which is believed to be due to the increase of water molecule bonding to two ether oxygens in meso dyad at low concentrations; while other conformers have less apparent variations. At higher concentrations the sizes and conformations are quite similar with those in bulk. And the structural variations of polymers are in qualitative agreement with the observations from the Raman spectra.
Solvation mechanisms different from small molecules were proposed for PVME solutions through the counting of HBs and QHBs. The transition of fractions of Op with HBs with water at around 0.3 coincides perfectly well with the results of conformational analysis and Raman spectra, which indicates the hydration structures are greatly influenced by the conformational changes in polymer. Fractions of water molecules hydrogen bonded to PVME increase with concentration and at high concentrations significant amount of water have been observed to be hydrogen bonded to two polymer repeat units which might be within the same chain or from two different chains.
The fractions of clusters around PVME repeat units with i (i=0-5) water molecules are well designated and it is found that the 3D HB networks take the majority in solutions at low concentrations, and regional clusters composed of 2-4 water molecules are common and single isolated HBs between water and ether oxygen are also seen at high concentrations.
The detailed method for the generation of atactic polymers with the most probable steric sequences is as follows. First, 1000 configurations of chains, each with 10000 repeat units, are generated as prototypes (we think so many polymer chains with such length are sufficient for the generation of the most probable steric sequences) with the specified steric ratio using the Monte Carlo method. The fractions of the continuous sequences of stereo dyads (such as mr, mmrr, …) are counted, for efficiency only sequences with utmost 4 dyads are counted. The same method is used for the generation of 1000 chains with 44 repeat units, their fractions of the continuous sequences of stereo dyads are compared with those of the prototypes by the root mean squared differences (RMSD), and the chain with the smallest RMSD will be regarded as the most probable configuration with the limited chain length.
This work was supported by the National Science Foundation of China (No. 21304017), the Fundamental Research Funds for the Central Universities and Innovation Program of Shanghai Municipal Education Commission. The calculations in this work were performed in the supercomputing center of Chinese Academy of Sciences. We are fully grateful to the provision of the PVME sample from Prof. Nies in KULeuven.
  1. Van Durme K,Loozen E,Nies E, Van Mele B (2005)Phase Behavior of Poly(vinyl methyl ether) in Deuterium Oxide. Macromolecules 38: 10234-10243.

  2. Schafer-Soenen H,Moerkerke R,Berghmans H,Koningsveld R,Dusek K,et al. (1997) Zero and Off-Zero Critical Concentrations in Systems Containing Polydisperse Polymers with Very High Molar Masses. 2. The System Water−Poly(vinyl methyl ether). Macromolecules 30: 410-416.

  3. Maeda HJ (1994)Interaction of water with poly(vinyl methyl ether) in aqueous solution. PolymSci Part B: PolymPhys 32: 91-97.

  4. Van Durme K, Van Assche G, Nies E, Van Mele B (2007) Phase transformations in aqueous low molar mass poly(vinyl methyl ether) solutions: theoretical prediction and experimental validation of the peculiar solvent melting line, bimodal LCST, and (adjacent) UCST miscibility gaps. J PhysChem B 111: 1288-1295.

  5. Sun B, Lai H, Wu P (2011) Integrated microdynamics mechanism of the thermal-induced phase separation behavior of poly(vinyl methyl ether) aqueous solution. J PhysChem B 115: 1335-1346.

  6. Némethy G,Scheraga HA (1962) Structure of Water and Hydrophobic Bonding in Proteins. II. Model for the Thermodynamic Properties of Aqueous Solutions of Hydrocarbons. J ChemPhys 36: 3401.

  7. Frank HS, Evans MW (1945) Free Volume and Entropy in Condensed Systems III. Entropy in Binary Liquid Mixtures; Partial Molal Entropy in Dilute Solutions; Structure and Thermodynamics in Aqueous Electrolytes. J ChemPhys 13: 507.

  8. Tamai YTH, Nakanishi K (1996)Molecular Dynamics Study of Polymer−Water Interaction in Hydrogels. 1. Hydrogen-Bond Structure. Macromolecules 29: 6750-6760.

  9. Meeussen F,Bauwens Y,Moerkerke R,Nies E,Berghmans H (2000)Molecular complex formation in the system poly(vinyl methyl ether)/water. Polymer 41: 3737-3743.

  10. Maeda Y (2001) IR Spectroscopic Study on the Hydration and the Phase Transition of Poly(vinyl methyl ether) in Water. Langmuir 17: 1737-1742.

  11. Zhang J, Berge B,Meeussen F,Nies E,Berghmans H,et al. (2003) Influence of the Interactions in Aqueous Mixtures of Poly(vinyl methyl ether) on the Crystallization Behavior of Water. Macromolecules 36: 9145-9153.

  12. Zeng X, Yang X (2004) A Study of Interaction of Water and Model Compound of Poly(vinyl methyl ether). J PhysChem B 108: 17384-17392.

  13. Loozen E,Nies E, Heremans K, Berghmans H (2006) The influence of pressure on the lower critical solution temperature miscibility behavior of aqueous solutions of poly(vinyl methyl ether) and the relation to the compositional curvature of the volume of mixing. J PhysChem B 110: 7793-7802.

  14. Nies E, Li T, Berghmans H, Heenan RK, King SM (2006) Upper critical solution temperature phase behavior, composition fluctuations, and complex formation in poly (vinyl methyl ether)/D2O solutions: small-angle neutron-scattering experiments and wertheim lattice thermodynamic perturbation theory predictions. J PhysChem B 110: 5321-5329.

  15. Xavier P, Bose S (2013) Multiwalled-carbon-nanotube-induced miscibility in near-critical PS/PVME blends: assessment through concentration fluctuations and segmental relaxation. J PhysChem B 117: 8633-8646.

  16. BharatiA, Xavier P,Kar GP, Madras G, Bose S (2014)Nanoparticle-Driven Intermolecular Cooperativity and Miscibility in Polystyrene/Poly(vinyl methyl ether) Blends. The Journal of Physical Chemistry B 118: 2214-2225.

  17. van der Spoel D, van Maaren PJ, Larsson P, Tîmneanu N (2006) Thermodynamics of hydrogen bonding in hydrophilic and hydrophobic media. J PhysChem B 110: 4393-4398.

  18. Kokubo H,Pettitt BM (2007) Preferential solvation in urea solutions at different concentrations: properties from simulation studies. J PhysChem B 111: 5233-5242.

  19. Zhang R, Li H, Lei Y, Han S (2005) All-atom molecular dynamic simulations and relative NMR spectra study of weak C-H...O contacts in amide-water systems. J PhysChem B 109: 7482-7487.

  20. Vishnyakov A,Lyubartsev AP,Laaksonen AJ (2001) Molecular Dynamics Simulations of Dimethyl Sulfoxide and Dimethyl Sulfoxide−Water Mixture. PhysChemA 105: 1702-1710.

  21. Jorgensen WL, Tirado-Rives J(1988) The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. J Am ChemSoc 110: 1657- 1666.

  22. Berendsen HJC,Grigera JR,Straatsma TP (1987)The missing term in effective pair potentials. J PhysChem 91: 6269-6271.

  23. Kaminski GA,Friesner RA,Tirado-Rives J, Jorgensen WL (2001) Evaluation and Reparametrization of the OPLS-AA Force Field for Proteins via Comparison with Accurate Quantum Chemical Calculations on Peptides. J PhysChem B 105: 6474-6487.

  24. Jorgensen WL, Maxwell DS,Tirado-Rives J (1996)Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids. J Am ChemSoc 118: 11225-11236.

  25. Jorgensen WL, Chandrasekhar J, Madura JD,Impey RW, Klein ML (1983)Comparison of simple potential functions for simulating liquid water. J ChemPhys 79: 926.

  26. Maaren PJ,SpoelD (2001) Molecular Dynamics Simulations of Water with Novel Shell-Model Potentials. J PhysChem B 105: 2618-2626.

  27. SpoelD,Maaren PJ (2006) The Origin of Layer Structure Artifacts in Simulations of Liquid Water. Journal of Chemical Theory and Computation 2: 1-11.

  28. Guillot B(2002) A reappraisal of what we have learnt during three decades of computer simulations on water. J MolLiq 101: 219-260.

  29. Van Der Spoel D,Lindahl E, Hess B, Groenhof G, Mark AE, et al. (2005) GROMACS: fast, flexible, and free. J ComputChem 26: 1701-1718.

  30. Ryckaert JP,Bellemans A (1978) Molecular dynamics of liquid alkanes. Far Disc ChemSoc 66: 95-106.

  31. Hoover WG (1985) Canonical dynamics: Equilibrium phase-space distributions. Phys Rev A 31: 1695-1697.

  32. Nosé S (1984) A molecular dynamics method for simulations in the canonical ensemble.MolPhys 52: 255-268.

  33. Parrinello M,Rahman A (1981) Polymorphic transitions in single crystals: A new molecular dynamics method. J ApplPhys52: 7182-7190.

  34. Linse B,Linse P (2014) Tuning the smooth particle mesh Ewald sum: application on ionic solutions and dipolar fluids. J ChemPhys 141: 184114.

  35. Darden T, York D, Pedersen L (1993) Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems. J ChemPhys 98: 10089.

  36. Hess B,Bekker H,Berendsen HJC,Fraaije JGEM (1997)LINCS: A linear constraint solver for molecular simulations. J ComputChem 18: 1463.

  37. Miyamoto S,Kollman PA (1992) Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models. JComputChem 13: 952-962.

  38. Luzar AM Chandler D (1996)Hydrogen-bond kinetics in liquid water. Nature 379: 55-57.

  39. Luzar A, Chandler D (1996) Effect of environment on hydrogen bond dynamics in liquid water. Phys Rev Lett 76: 928-931.

  40. Wu R,Ji Q,Kong B, Yang X (2008)Molecular dynamics simulations of the hydration of poly(vinyl methyl ether): Hydrogen bonds and quasi-hydrogen bonds. Sci China Ser B-Chem 51: 736-742.

  41. Eisenberg D,Kauzmann W (1969)The Structure and Properties of Water; Oxford University Press: London, UK.

  42. Shiomi T, Hamada F,Nasako T,Yoneda K, Imai K,et al. (1990) Thermodynamics of polymer blends of poly(vinyl methyl ether) and polystyrene. Macromolecules 23: 229-233.

  43. Byutner OG, Smith GD (2000) Quantum Chemistry Based Force Field for Simulations of Poly(vinylidene fluoride). Macromolecules 33: 4264-4270.

  44. Redlich O,Kister AT (1948) Algebraic Representation of Thermodynamic Properties and the Classification of Solutions.IndEngChem 40: 345-348.

  45. Gubskaya AV,Kusalik PG (2004) Molecular Dynamics Simulation Study of Ethylene Glycol, Ethylenediamine, and 2-Aminoethanol. 2. Structure in Aqueous Solutions. J PhysChemA 108: 7165-7178.

  46. Laaksonen A,Kusalik PG,Svishchev IM (1997) Three-Dimensional Structure in Water−Methanol Mixtures. J PhysChemA101: 5910-5918.

  47. Yang X, Kang S, Hsu SL,Stidham HD, Smith PB,et al. (2001) A Spectroscopic Analysis of Chain Flexibility of Poly(lactic acid). Macromolecules 34: 5037-5041.

  48. Yang X, Su Z, Wu D, Hsu SL,Stidham HD (1997)Raman Analysis of a Conformational Distribution of Poly(ethylene oxide) and Its Model Compound in the Liquid State. Macromolecules 30: 3796-3802.

  49. Scheiner S (1997) Hydrogen Bonding; Oxford University Press: New York.

Tables and Figures at a glance

Table icon Table icon
Table 1 Table 2


Figures at a glance

Figure Figure Figure Figure Figure
Figure 1 Figure 2 Figure 3 Figure 4 Figure 5


Figure Figure Figure Figure Figure
Figure 6 Figure 7 Figure 8 Figure 9 Figure 10


Figure Figure Figure
Figure 11 Figure 12 Figure 13
Select your language of interest to view the total content in your interested language
Post your comment

Share This Article

Article Usage

  • Total views: 8047
  • [From(publication date):
    December-2015 - Oct 24, 2017]
  • Breakdown by view type
  • HTML page views : 7978
  • PDF downloads :69

Post your comment

captcha   Reload  Can't read the image? click here to refresh

Peer Reviewed Journals
Make the best use of Scientific Research and information from our 700 + peer reviewed, Open Access Journals
International Conferences 2017-18
Meet Inspiring Speakers and Experts at our 3000+ Global Annual Meetings

Contact Us

Agri, Food, Aqua and Veterinary Science Journals

Dr. Krish

1-702-714-7001 Extn: 9040

Clinical and Biochemistry Journals

Datta A

1-702-714-7001Extn: 9037

Business & Management Journals


1-702-714-7001Extn: 9042

Chemical Engineering and Chemistry Journals

Gabriel Shaw

1-702-714-7001 Extn: 9040

Earth & Environmental Sciences

Katie Wilson

1-702-714-7001Extn: 9042

Engineering Journals

James Franklin

1-702-714-7001Extn: 9042

General Science and Health care Journals

Andrea Jason

1-702-714-7001Extn: 9043

Genetics and Molecular Biology Journals

Anna Melissa

1-702-714-7001 Extn: 9006

Immunology & Microbiology Journals

David Gorantl

1-702-714-7001Extn: 9014

Informatics Journals

Stephanie Skinner

1-702-714-7001Extn: 9039

Material Sciences Journals

Rachle Green

1-702-714-7001Extn: 9039

Mathematics and Physics Journals

Jim Willison

1-702-714-7001 Extn: 9042

Medical Journals

Nimmi Anna

1-702-714-7001 Extn: 9038

Neuroscience & Psychology Journals

Nathan T

1-702-714-7001Extn: 9041

Pharmaceutical Sciences Journals

John Behannon

1-702-714-7001Extn: 9007

Social & Political Science Journals

Steve Harry

1-702-714-7001 Extn: 9042

© 2008-2017 OMICS International - Open Access Publisher. Best viewed in Mozilla Firefox | Google Chrome | Above IE 7.0 version