Computational Calorimetry of PNIPAM Cononsolvency in Water/Methanol Mixtures

Mechanism of poly(N-isopropylacrylamide) cononsolvency in aqueous methanol solutions explored via oxygen K-edge X-ray absorption spectroscopy

The cononsolvency mechanism of poly(N-isopropylacrylamide) (PNIPAM), dissolving in pure methanol (MeOH) and water (H2O) but being insoluble in MeOH-H2O mixtures, was investigated by O K-edge X-ray absorption spectroscopy (XAS). The cononsolvency emerges from the aggregation of PNIPAM with MeOH clusters, leading to the collapse of the hydrophobic hydration of PNIPAM.

Distinctly different solvation behaviors of poly(N,N-diethylacrylamide) gels in water/acetone and water/DMSO mixtures

The solvation behaviors and intermolecular interactions of a poly(N,N-diethylacrylamide) (PDEA) gel network in water/DMSO and water/acetone mixtures have been investigated by variable-temperature high-resolution ¹H MAS NMR. Unlike decreasing volume phase transition temperature (VPTT) of the typical thermosensitive poly(N-isopropylacrylamide) (PNIPAM) gel induced by both acetone and DMSO in a water-rich region, distinct phase transition behaviors are revealed for the PDEA gel. That is, acetone is found to increase the VPTT of PDEA directly in the water-rich region while DMSO is also found to increase the VPTT of PDEA at a very low concentration but then decrease the VPTT as the concentration further increases. The above distinctly different VPTT shifts of PDEA are attributed to the different polymer-cosolvent interactions in water/acetone and water/DMSO systems. DMSO molecules with a strong water affinity are preferentially excluded by the PDEA gel network, and can promote the volume phase transition by favoring the collapse of the PDEA gel network, while acetone molecules preferentially adsorbed on the polymer interact with PDEA via non-specific van der Waals interaction, which favors the swollen state of the PDEA gel.

Poly(sulfobetaine)-Based Diblock Copolymer Thin Films in Water/Acetone Atmosphere: Modulation of Water Hydration and Co-nonsolvency-Triggered Film Contraction

The water swelling and subsequent solvent exchange including co-nonsolvency behavior of thin films of a doubly thermo-responsive diblock copolymer (DBC) are studied via spectral reflectance, time-of-flight neutron reflectometry, and Fourier transform infrared spectroscopy. The DBC consists of a thermo-responsive zwitterionic (poly(4-((3-methacrylamidopropyl) dimethylammonio) butane-1-sulfonate)) (PSBP) block, featuring an upper critical solution temperature transition in aqueous media but being insoluble in acetone, and a nonionic poly(N-isopropylmethacrylamide) (PNIPMAM) block, featuring a lower critical solution temperature transition in water, while being soluble in acetone. Homogeneous DBC films of 50-100 nm thickness are first swollen in saturated water vapor (H₂O or D₂O), before they are subjected to a contraction process by exposure to mixed saturated water/acetone vapor (H₂O or D₂O/acetone-d6 = 9:1 v/v). The affinity of the DBC film toward H₂O is stronger than for D₂O, as inferred from the higher film thickness in the swollen state and the higher absorbed water content, thus revealing a pronounced isotope sensitivity. During the co-solvent-induced switching by mixed water/acetone vapor, a two-step film contraction is observed, which is attributed to the delayed expulsion of water molecules and uptake of acetone molecules. The swelling kinetics are compared for both mixed vapors (H₂O/acetone-d6 and D₂O/acetone-d6) and with those of the related homopolymer films. Moreover, the concomitant variations of the local environment around the hydrophilic groups located in the PSBP and PNIPMAM blocks are followed. The first contraction step turns out to be dominated by the behavior of the PSBP block, whereas the second one is dominated by the PNIPMAM block. The unusual swelling and contraction behavior of the latter block is attributed to its co-nonsolvency behavior. Furthermore, we observe cooperative hydration effects in the DBC films, that is, both polymer blocks influence each other's solvation behavior.

Thermoreversible Gelation of Associating Polymers in Hydrogen-Bonding Mixed Solvents

The thermoreversible gelation of associating polymers in hydogen-bonding mixed solvents is studied theoretically on the basis of the model ternary solution in which polymers carrying functional group A are dissolved in a primary solvent SB (mainly water) and cosolvent SC. Functional group A forms intermolecular cross-links (A-A) leading to gelation. It can also form hydrogen bonds with either solvent molecule SB or SC (A-SB, A-SC). If bound to a solvent molecule, A becomes inert. Solvent molecules form complexes with variable association numbers (SB-SC). Ternary phase diagrams for sol-gel transition lines and their cross sections on the temperature-solvent composition plane are constructed. It is shown that a minimum in the polymer concentration along the sol-gel transition line appears at a special solvent composition which depends upon the structure of the complexes. At such an optimal gel point, the sol-gel transition temperature becomes the highest, thus leading to stable, strong gels. The degree of adsorption, preferential adsorption coefficient, average molecular weight of the complexes, and cross-link probability are calculated as functions of the solvent composition. At the optimal gel point, the preferential adsorption coefficient changes its sign and the average molecular weight of the complexes reaches a maximum. Results are compared with the experimental data on poly(vinyl alcohol) in a mixed solvent of water/dimethyl sulfoxide, dimethyl sulfoxide/urea, or water/methanol.

Solvation shell thermodynamics of extended hydrophobic solutes in mixed solvents

The ability of various cosolutes and cosolvents to enhance or quench solvent density fluctuations at solute–water interfaces has crucial implications on the conformational equilibrium of macromolecules such as polymers and proteins. Herein, we use an extended hydrophobic solute as a model system to study the effect of urea and methanol on the density fluctuations in the solute’s solvation shell and the resulting thermodynamics. On strengthening the solute–water/cosolute repulsive interaction, we observe distinct trends in the mutual affinities between various species in, and the thermodynamic properties of, the solvation shell. These trends strongly follow the respective trends in the preferential adsorption of urea and methanol: solute–water/cosolute repulsion strengthens, urea accumulation decreases, and methanol accumulation increases. Preferential accumulation of urea is found to quench the density fluctuations around the extended solute, leading to a decrease in the compressibility of the solvation shell. In contrast, methanol accumulation enhances the density fluctuations, leading to an increase in the compressibility. The mode of action of urea and methanol seems to be strongly coupled to their hydration behavior. The observations from this simple model is discussed in relation to urea driven swelling and methanol induced collapse of some well-known thermo-responsive polymers.

Cononsolvency of thermoresponsive polymers: where we are now and where we are going

Cononsolvency is an intriguing phenomenon where a polymer collapses in a mixture of good solvents. This cosolvent-induced modulation of the polymer solubility has been observed in solutions of several polymers and biomacromolecules, and finds application in areas such as hydrogel actuators, drug delivery, compound detection and catalysis. In the past decade, there has been a renewed interest in understanding the molecular mechanisms which drive cononsolvency with a predominant emphasis on its connection to the preferential adsorption of the cosolvent. Significant efforts have also been made to understand cononsolvency in complex systems such as micelles, block copolymers and thin films. In this review, we will discuss some of the recent developments from the experimental, simulation and theoretical fronts, and provide an outlook on the problems and challenges which are yet to be addressed.

Linear shrinkage of hydrogel coatings exposed to flow: interplay between dissolution of water and advective transport

We investigate the shrinkage of a surface-grafted water-swollen hydrogel under shear flows of oils by laser scanning confocal microscopy. Interestingly, external shear flows of oil lead to linear dehydration and shrinkage of the hydrogel for all investigated flow conditions irrespective of the chemical nature of the hydrogel. The reason is that the finite solubility of water in oil removes water from the hydrogel continuously by diffusion. The flow advects the water-rich oil, as demonstrated by numerical solutions of the underlying convection-diffusion equation. In line with this hypothesis, shear does not cause gel shrinkage for water-saturated oils or non-solvents. The solubility of water in the oil will tune the dehydration dynamics.

Nonadditive ion effects on the coil-globule equilibrium of PNIPAM: A computer simulation study

The combined effect of a weakly hydrated and a strongly hydrated anion on the lower critical solution temperature (LCST) of poly(N-isopropylacrylamide)(PNIPAM) is nonadditive (Bruce et al. J. Am. Chem. Soc....

Systematic Comparison of the Structural and Dynamic Properties of Commonly Used Water Models for Molecular Dynamics Simulations

Water is a unique solvent that is ubiquitous in biology and present in a variety of solutions, mixtures, and materials settings. It therefore forms the basis for all molecular dynamics simulations of biological phenomena, as well as for many chemical, industrial, and materials investigations. Over the years, many water models have been developed, and it remains a challenge to find a single water model that accurately reproduces all experimental properties of water simultaneously. Here, we report a comprehensive comparison of structural and dynamic properties of 30 commonly used 3-point, 4-point, 5-point, and polarizable water models simulated using consistent settings and analysis methods. For the properties of density, coordination number, surface tension, dielectric constant, self-diffusion coefficient, and solvation free energy of methane, models published within the past two decades consistently show better agreement with experimental values compared to models published earlier, albeit with some notable exceptions. However, no single model reproduced all experimental values exactly, highlighting the need to carefully choose a water model for a particular study, depending on the phenomena of interest. Finally, machine learning algorithms quantified the relationship between the water model force field parameters and the resulting bulk properties, providing insight into the parameter-property relationship and illustrating the challenges of developing a water model that can accurately reproduce all properties of water simultaneously.

Length-Scale Effects in Hydrophobic Polymer Collapse Transitions

The study of molecular mechanisms for cosolvent-driven hydrophobic polymer collapse transitions in water is of pivotal importance in the field of smart responsive materials. Computational studies together with complementary experimental data have led to the discovery and understanding of new phenomena in recent years. However, elementary mechanisms, generally contributing to polymer coil-globule transitions in different classes of cosolvent-water systems, remain elusive due to compensating energy-entropy effects. Herein, I discuss the role of length scales in polymer solubility problems. New ideas on surfactant mechanisms are discussed based on examples in which these mechanisms drive polymer swelling or collapse.

An interplay of excluded-volume and polymer-(co)solvent attractive interactions regulates polymer collapse in mixed solvents

Cosolvent effects on the coil-globule transitions in aqueous polymer solutions are not well understood, especially in the case of amphiphilic cosolvents that preferentially adsorb on the polymer and lead to both polymer swelling and collapse. Although a predominant focus in the literature has been placed on the role of polymer-cosolvent attractive interactions, our recent work has shown that excluded-volume interactions (repulsive interactions) can drive both preferential adsorption of the cosolvent and polymer collapse via a surfactant-like mechanism. Here, we further study the role of polymer-(co)solvent attractive interactions in two kinds of polymer solutions, namely, good solvent (water)-good cosolvent (alcohol) (GSGC) and poor solvent-good cosolvent (PSGC) solutions, both of which exhibit preferential adsorption of the cosolvent and a non-monotonic change in the polymer radius of gyration with the addition of the cosolvent. Interestingly, at low concentrations, the polymer-(co)solvent energetic interactions oppose polymer collapse in the GSGC solutions and contrarily support polymer collapse in the PSGC solutions, indicating the importance of the underlying polymer chemistry. Even though the alcohol molecules are preferentially adsorbed on the polymer, the trends of the energetic interactions at low cosolvent concentrations are dominated by the polymer-water energetic interactions in both the cases. Therefore, polymer-(co)solvent energetic interactions can either reinforce or compensate the surfactant-like mechanism, and it is this interplay that drives coil-to-globule transitions in polymer solutions. These results have implications for rationalizing the cononsolvency transitions in real systems such as polyacrylamides in aqueous alcohol solutions where the understanding of microscopic driving forces is still debatable.

Is the Microgel Collapse a Two-Step Process? Exploiting Cononsolvency to Probe the Collapse Dynamics of Poly-N-isopropylacrylamide (pNIPAM)

Many applications of responsive microgels rely on the fast adaptation of the polymer network. However, the underlying dynamics of the de-/swelling process of the gels have not been fully understood. In the present work, we focus on the collapse kinetics of poly-N-isopropylacrylamide (pNIPAM) microgels due to cononsolvency. Cononsolvency means that either of the pure solvents, e.g., pure water or pure methanol, act as a so-called good solvent, leading to a swollen state of the polymer network. However, in mixtures of water and methanol, the previously swollen network undergoes a drastic volume loss. To further elucidate the cononsolvency transition, pNIPAM microgels with diameters between 20 and 110 μm were synthesized by microfluidics. To follow the dynamics, pure water was suddenly exchanged with an unfavorable mixture of 20 mol% methanol (solvent-jump) within a microfluidic channel. The dynamic response of the microgels was investigated by optical and fluorescence microscopy and Raman microspectroscopy. The experimental data provide unique and detailed insight into the size-dependent kinetics of the volume phase transition due to cononsolvency. The change in the microgel's diameter over time points to a two-step process of the microgel collapse with a biexponential behavior. Furthermore, the dependence between the two time constants from this biexponential behavior and the microgel's diameter in the collapsed state deviates from the square-power law proposed by Tanaka and Fillmore [ J. Chem. Phys. 1979, 70, 1214-1218]. The deviation is discussed considering the adhesion-induced deformation of the gels and the physical processes underlying the collapse.

Conformational Stability of Poly (N-Isopropylacrylamide) Anchored on the Surface of Gold Nanoparticles

To verify the temperature sensitive failure of poly (N-isopropylacrylamide) (PNIPAM) anchored on the surface of gold nanoparticles (AuNPs), the UV-Vis spectra with temperature variations of the following aqueous solutions respectively containing AuNPs-PNIPAM, Au-PNIPAM/PNIPAM, PNIPAM, in different media (including salt, ethanol, HCl and cetyltrimethylammoniumbromide (CTAB)), were systematically determined. The results indicated that the UV-Vis spectrum of AuNPs-PNIPAM suspension hardly changed even above the Lower Critical Solution Temperature (LCST) of PNIPAM, but that of Au-PNIPAM/PNIPAM sharply increased only in absorbance intensity. A possible mechanism of the failed temperature sensitivity of PNIPAM anchored on the surface of AuNPs was proposed. Being different from free PNIPAM molecules, a strong interaction exists among PNIPAM molecules anchored on the surface of AuNPs, restraining the change in conformation of PNIPAM. The temperature sensitivity of Au-PNIPAM/PNIPAM originates from the free PNIPAM molecules rather than the anchored PNIPAM one. The changing electrostatic interaction could effectively regulate the aggregation behavior of AuNPs-PNIPAM and enhance its sensitivity to temperature.

A cosolvent surfactant mechanism affects polymer collapse in miscible good solvents

The coil-globule transition of aqueous polymers is of profound significance in understanding the structure and function of responsive soft matter. In particular, the remarkable effect of amphiphilic cosolvents (e.g., alcohols) that leads to both swelling and collapse of stimuli-responsive polymers has been hotly debated in the literature, often with contradictory mechanisms proposed. Using molecular dynamics simulations, we herein demonstrate that alcohols reduce the free energy cost of creating a repulsive polymer-solvent interface via a surfactant-like mechanism which surprisingly drives polymer collapse at low alcohol concentrations. This hitherto neglected role of interfacial solvation thermodynamics is common to all coil-globule transitions, and rationalizes the experimentally observed effects of higher alcohols and polymer molecular weight on the coil-to-globule transition of thermoresponsive polymers. Polymer-(co)solvent attractive interactions reinforce or compensate this mechanism and it is this interplay which drives polymer swelling or collapse.

On-Off of Co-non-solvency for Poly(N-vinylcaprolactam) in Alcohol-Water Mixtures: A Molecular Dynamics Study

Poly(N-vinylcaprolactam) (PVCL) exhibits co-non-solvency in aqueous solutions of 2-propanol but not in methanol. What distinguishes the impact of these two cosolvents on the polymer conformational stability? We report a molecular dynamics simulation study on PVCL 50-mer and monomers dissolved in methanol-water and 2-propanol-water mixtures. We show that the alcohol-concentration dependence of the effective attraction between a pair of PVCL monomers closely resembles the conformational changes in a single PVCL 50-mer as well as the experimentally observed behavior for PVCL chains. We also found that, at the co-non-solvency maximum, the monomer-monomer attraction works over a long-range beyond the solvent-separated distance. Then, we correlate the long-range attraction to the appearance of a dense alcohol concentration accumulated between the monomers. Furthermore, we distinctly demonstrate that the co-non-solvency of PVCL monomers can be switched on/off by artificially tuning the alcohol size while keeping the energetic parameters. Thus, we conclude that the magnitude of the excluded volume effect in alcohol accompanying the gain in translational entropy of the solvent is crucial to the occurrence of PVCL polymer co-non-solvency.

Small crowder interactions can drive hydrophobic polymer collapse as well as unfolding

Biomolecules evolve and function in the intracellular crowded environment that is densely packed with macromolecules. Yet, a microscopic understanding of the effects of such an environment on the conformational preferences of biomolecules remains elusive. While prior investigations have attributed crowding effects mainly to the excluded volume (size) effects of the crowders, very little is known about the effects exerted due to their chemical interactions. In this study, crowding effects of tri-alanine peptides on the collapse equilibria of generic hydrophobic polymer are investigated using molecular dynamics simulations. The role of weak, non-specific, attractive polymer-crowder interactions in modulating the polymer collapse equilibria is examined. The results highlight that crowding effects can lead to polymer compaction as well as unfolding depending on the strength of polymer-crowder interaction energy. Strongly interacting crowders weaken hydrophobic collapse (or unfold the polymer) at high volume fractions and induce polymer collapse only under dilute conditions. Weakly interacting crowders induce polymer collapse at all crowder concentrations. Interestingly, the thermodynamic driving forces for polymer collapse are remarkably different in the two cases. Strongly and weakly interacting crowders induce collapse by preferential adsorption and preferential depletion respectively. The findings provide new insights into the possible effects of interplay of intermolecular interactions in a crowded environment. The results have implications in understanding the impact of crowding in altering free energy landscapes of proteins.

SJ, van der Vegt NFA. Characterizing Polymer Hydration Shell Compressibilities with the Small-System Method

The small-system method (SSM) exploits the unique feature of finite-sized open systems, whose thermodynamic quantities scale with the inverse system size. This scaling enables the calculation of properties in the thermodynamic limit of macroscopic systems based on computer simulations of finite-sized systems. We herein extend the SSM to characterize the hydration shell compressibility of a generic hydrophobic polymer in water. By systematically increasing the strength of polymer-water repulsion, we find that the excess inverse thermodynamic correction factor (Δ1/Γs∞) and compressibility (Δχs) of the first hydration shell change sign from negative to positive. This occurs with a concurrent decrease in water hydrogen bonding and local tetrahedral order of the hydration shell water. The crossover lengthscale corresponds to an effective polymer bead diameter of 0.7 nm and is consistent with previous works on hydration of small and large hydrophobic solutes. The crossover lengthscale in polymer hydration shell compressibility, herein identified with the SSM approach, relates to hydrophobic interactions and macromolecular conformational equilibria in aqueous solution. The SSM approach may further be applied to study thermodynamic properties of polymer solvation shells in mixed solvents.

Thermodynamics of N-Isopropylacrylamide in Water: Insight from Experiments, Simulations, and Kirkwood-Buff Analysis Teamwork

The behavior of thermoresponsive polymer poly(N-isopropylacrylamide) (PNiPAM), an essential building block in the design of smart soft materials, in aqueous solutions has attracted much interest, which contrasts with our knowledge of N-isopropylacrylamide (NiPAM) monomer. Strikingly, the physicochemical properties of aqueous NiPAM are similarly rich, and their understanding is far from being complete. This stems from the lack of accurate thermodynamic data and quantitative model for atomistic simulations. In this joint study, we have probed the thermodynamic behavior of aqueous NiPAM by experimental methods, molecular dynamics (MD) simulations, and Kirkwood-Buff (KB) analysis at ambient conditions. From the partial molar volumes and simultaneously correlated osmotic coefficients, with excess partial molar enthalpies of NiPAM in water, the concentration and temperature dependence of KB integrals was determined. For the purpose of this work, we have developed and employed a novel NiPAM force field, which not only reproduces KB integrals (G_ij) and adequately captures macroscopic thermodynamic quantities but also provides more accurate structural insight than the original force fields. We revealed in the vicinity of NiPAM the competing effect of amide hydration with interaction between nonpolar regions. This microscopic picture is reflected in the experimentally observed NiPAM-NiPAM association, which is present from highly dilute conditions up to the solubility limit and is evidenced by G₂₂. From intermediate concentrations, it is accompanied by the existence of apparent dense-water regions, as indicated by positive G₁₁ values. The here-employed KB-based framework provided a mutually consistent thermodynamic and microscopic insight into the NiPAM solution and may be further extended for ion-specific effects. Moreover, our findings contribute to the understanding of thermodynamic grounds behind PNiPAM collapse transition.

Effect of sodium thiocyanate and sodium perchlorate on poly(N-isopropylacrylamide) collapse

The T(collapse) of poly(N-isopropylacrylamide), PNIPAM, shows a nonlinear dependence on the concentration of NaSCN or NaClO4; in the case of NaClO4, for example, at very low concentrations of the salt, T(collapse) increases with the concentration, while it has an opposite trend at higher NaClO4 concentrations [J. Am. Chem. Soc., 2005, 127, 14505]. These puzzling experimental data can be rationalized by considering that low charge density and poorly hydrated ions, such as thiocyanate and perchlorate, interact preferentially with the surface of the polymer, and cause an increase of the magnitude of the energetic term that stabilizes swollen conformations at low salt concentrations. However, as both swollen and collapsed PNIPAM conformations are accessible to such ions in view of their large conformational freedom, the difference in the number of ions bound to PNIPAM surface upon collapse changes little on increasing the salt concentration. Thus, the energetic term that favors swollen conformations increases with salt concentration to a lesser extent than the solvent-excluded volume term (linked to the density increase caused by salt addition to water), that favors collapsed conformations, leading to a nonlinear trend of T(collapse).

Distinct and Nonadditive Effects of Urea and Guanidinium Chloride on Peptide Solvation

Using enhanced-sampling replica exchange fully atomistic molecular dynamics simulations, we show that, individually, urea and guanidinium chloride (GdmCl) denature the Trpcage protein, but remarkably, the helical segment ¹NLYIQWL⁷ of the protein is stabilized in mixed denaturant solutions. GdmCl induces protein denaturation via a combination of direct and indirect effects involving dehydration of the protein and destabilization of stabilizing salt bridges. In contrast, urea denatures the protein through favorable protein-urea preferential interactions, with peptide-specific indirect effects of urea on the water structure around the protein. In the case of the helical segment of Trpcage, urea "oversolvates" the peptide backbone by reorganizing water molecules from the peptide side chains to the peptide backbone. An intricate nonadditive thermodynamic balance between GdmCl-induced dehydration of the peptide and the urea-induced changes in solvation structure triggers partial counteraction to urea denaturation and stabilization of the helix.

Many-chain effects on the co-nonsolvency of polymer brushes in a good solvent mixture

Polymer brushes normally swell in a good solvent and collapse in a poor solvent. An abnormal response of polymer brushes, so-called co-nonsolvency, is the phenomenon where the brush counter-intuitively collapses in a good solvent mixture. In this work, we employed molecular dynamics simulations to investigate the structural properties of the grafted polymers in the occurrence of co-nonsolvency. Brushes with various grafting densities were considered to study the effect of topologically excluded volumes on the co-nonsolvency. We found that the brush height follows a novel scaling behavior of the grafting density h ∼ σg0.71 in the co-nonsolvent mixture. Using the scaling exponent and Alexander-de Gennes theory, an analytic function that predicts the monomer density was obtained. The many-chain effects in the co-nonsolvent lead to the formation of both intermolecular and intramolecular bridging structures. Increasing the grafting density entails lower looping events occuring because of the intermolcular bridging, causing diverse structural properties. We report how the average thickness, the polymer orientation, and the looping probability vary as the grafting density increases. Based on these observations, we constructed a phase diagram of the polymer brush system using the average thickness and orientation as order parameters. Our simulations and analytical results reveal the nature of co-nonsolvency in polymer brushes in an explicit way and will help to provide practical guidelines for applications such as drug delivery and sensor devices.

Coil-Globule Transition Thermodynamics of Poly(N-isopropylacrylamide)

Thermosensitive polymers such as poly(N-isopropylacrylamide) (PNIPAM) undergo a phase transition in aqueous solution from a random-coil structural ensemble to a globule structural ensemble at the lower critical solution temperature (LCST). Above this temperature, PNIPAM agglomerates and becomes insoluble, whereas it is soluble below the temperature. Thus, thermosensitive polymers represent essential targets for several applications, e.g., in drug delivery. Although their ability to change structure in response to a temperature alteration is highly relevant for industrial processes, their thermodynamic properties are mostly qualitatively understood, and the quantitative thermodynamic picture is still elusive. In this study, we used a combined atomistic molecular dynamics and well-tempered metadynamics simulation approach to estimate coil-globule transition thermodynamics. An isotactic 30-mer of PNIPAM was investigated over a broad temperature range between 200 and 360 K. The transition from the globule to the random-coil structure was observed with well-tempered metadynamics. For the first time, the free energy surface of PNIPAM was estimated and it is shown that the simulation results are in line with the experimentally observed thermosensitive behavior. Below the LCST, the random-coil ensemble represents the global energy minimum and is thermodynamically favored by 21 ± 9 kJ/mol compared to the globule ensemble; both are separated by a barrier of 49 ± 14 kJ/mol. In contrast, above the LCST, the globule ensemble is thermodynamically favored by 21 ± 8 kJ/mol over the random-coil ensemble. The barrier from random-coil to globule is 17 ± 10 kJ/mol.

Poly(2-propyl-2-oxazoline)s in Aqueous Methanol: To Dissolve or not to Dissolve

At room temperature, poly(N-isopropylacrylamide) (PNIPAM) is soluble in water and methanol, but it is not soluble in certain water/methanol mixtures. This phenomenon, known as cononsolvency, has been explored in great detail experimentally and theoretically in an attempt to understand the complex interactions occurring in the ternary PNIPAM/water/co-nonsolvent system. Yet little is known about the effects of the polymer structure on cononsolvency. To address this point, we investigated the temperature-dependent solution properties in water, methanol, and mixtures of the two solvents of poly(2-cyclopropyl-2-oxazoline) (PcyPOx) and two structural isomers of PNIPAM (M_n ∼ 11 kg/mol): poly(2-isopropyl-2-oxazoline) (PiPOx) and poly(2-n-propyl-2-oxazoline) (PnPOx). The phase diagram of the ternary water/methanol/poly(2-propyl-2-oxazolines) (PPOx) systems, constructed based on cloud point (T_CP) measurements, revealed that PnPOx exhibits cononsolvency in water/methanol mixtures. In contrast, methanol acts as a cosolvent for PiPOx and PcyPOx in water. The enthalpy, ΔH, and temperature, T_max, of the coil-to-globule transition of the three polymers in various water/methanol mixtures were measured by high-sensitivity differential scanning calorimetry. T_max follows the same trends as T_CP, confirming the cononsolvency of PnPOx and the cosolvency of PiPOx and PcyPOx. ΔH decreases linearly as a function of the methanol content for all PPOx systems. Ancillary high-resolution ¹H NMR spectroscopy studies of PPOx solutions in D₂O and methanol-d₄, coupled with DOSY and NOESY experiments revealed that the n-propyl group of PnPOx rotates freely in D₂O, whereas the rotation of the isopropyl and cyclopropyl groups of PiPOx and PcyPOx, respectively, is limited due to steric restriction. This factor appears to play an important role in the case of the PPOxs/water/methanol ternary system.

Improved Temperature Behavior of PNIPAM in Water with a Modified OPLS Model

We test the OPLS/AA force field for a single PNIPAM 40-mer in aqueous solution using replica exchange molecular dynamics simulations and find that the force field fails to reproduce the experimental temperature behavior. To resolve this issue, we apply a modification on the partial charges previously suggested to reproduce the liquid-liquid phase separation of NIPAM aqueous solutions. The modified force field features stronger amide-water electrostatic interactions than the original OPLS model, predicts a weaker water-mediated monomer-monomer attraction, and reproduces the experimental coil-globule collapse enthalpy of PNIPAM in water. We revisit the cononsolvency problem of PNIPAM in methanol/water mixtures with the modified model and show that the dependence of the coil-globule collapse enthalpy on methanol concentration follows the experimental trend of the lower critical solution temperature. The calculations with the modified force field confirm that polymer dehydration is the determining factor for chain collapse in the cononsolvency regime.

Self-Assembly of Short Chain Poly- N-isopropylacrylamid Induced by Superchaotropic Keggin Polyoxometalates: From Globules to Sheets

We show here for the first time that short chain poly( N-isopropylacrylamide) (PNIPAM), one of the most famous thermoresponsive polymers, self-assembles in water to form (i) discrete nanometer-globules and (ii) micrometric sheets with nm-thickness upon addition of the well-known Keggin-type polyoxometalate (POM) H₃PW₁₂O₄₀ (PW). The type of self-assembly is controlled by PW concentration: at low PW concentrations, PW adsorbs on PNIPAM chains to form globules consisting of homogeneously distributed PWs in PNIPAM droplets of several nm in size. Upon further addition of PW, a phase transition from globules to micrometric sheets is observed for PNIPAMs above a polymer critical chain length, between 18 and 44 repeating units. The thickness of the sheets is controlled by the PNIPAM chain length, here from 44 to 88 repeating units. The PNIPAM sheets are electrostatically stabilized PWs accumulated on each side of the sheets. The shortest PNIPAM chain with 18 repeating units produces PNIPAM/PW globules with 5-20 nm size but no sheets. The PW/PNIPAM self-assembly arises from a solvent mediated mechanism associated with the partial dehydration of PW and of the PNIPAM, which is related to the general propensity of POMs to adsorb on neutral hydrated surfaces. This effect, known as superchaotropy, is further highlighted by the significant increase in the lower critical solubilization temperature (LCST) of PNIPAM observed upon the addition of PW in the mM range. The influence of the POM nature on the self-assembly of PNIPAM was also investigated by using H₄SiW₁₂O₄₀ (SiW) and H₃PMo₁₂O₄₀ (PMo), i.e. changing the POM's charge density or polarizability in order to get deeper understanding on the role of electrostatics and polarizability in the PNIPAM self-assembly process. We show here that the superchaotropic behavior of POMs with PNIPAM polymers enables the formation and the shape control of supramolecular organic-inorganic hybrids.

Enrichment of methanol inside pNIPAM gels in the cononsolvency-induced collapse

From Raman, we determined an enrichment of methanol inside the polymer in the cononsolvency-induced collapse and donor-type hydrogen-bonding of methanol with pNIPAM.

Cross-linker effect on solute adsorption in swollen thermoresponsive polymer networks

Molecular dynamics study on the solute adsorption to thermoresponsive polymers estimating the cross-link impact on particle partitioning in swollen hydrogels.

P-NIPAM in water–acetone mixtures: experiments and simulations

The lower critical solution temperature (LCST) of poly-N-isopropylacrylamide (p-NIPAM) diminishes when a small volume of acetone is added to the aqueous polymer solution, and then increases for further additions, producing a minimum at a certain acetone concentration.

Studying polymer solutions with particle-based models linked to classical density functionals: co-non-solvency

We demonstrate the potential of hybrid particle-based models, where interactions are introduced through functionals of local order parameters, in describing multicomponent polymer solutions. The link to a free-energy-like functional is advantageous for controlling the thermodynamics of the model. We focus on co-non-solvency - the collapse of polymer chains in dilute mixtures with two miscible good solvents, having different affinities towards the polymer. We employ a simple model where polymers and solvents are represented, respectively, by worm-like chains and single particles. Non-bonded interactions are captured by a polynomial which is third order in local densities and can, therefore, describe liquid-vapour coexistence. The parameterisation of the functional benefits from an elementary mean-field approximation to the statistical mechanics of the model. The model provides a framework for Monte Carlo simulations using a particle-to-mesh algorithm. Studies with conventional generic bead-spring and all-atom models have demonstrated that co-non-solvency is caused by preferential binding of the better solvent (termed cosolvent) with polymer. Hence, segmental loops bridged by cosolvent molecules are formed, initiating polymer collapse. The mesoscopic hybrid model differs conceptually from the conventional microscopic descriptions. Yet, it reproduces the same co-non-solvency mechanism supporting its universality. Films of adsorbed ternary solutions, showing co-non-solvency in the dilute regime, are considered at high concentrations. In this case, chains do not collapse. The properties of loops and tails of the adsorbed polymer agree with early theoretical predictions obtained for concentrated binary solutions.

Cononsolvency Transition of Polymer Brushes: A Combined Experimental and Theoretical Study

In this study, the cononsolvency transition of poly(N-isopropylacrylamide) (PNiPAAm) brushes in aqueous ethanol mixtures was studied by using Vis-spectroscopic ellipsometry (SE) discussed in conjunction with the adsorption-attraction model. We proved that the cononsolvency transition of PNiPAAm brushes showed features of a volume phase transition, such as a sharp collapse, reaching a maximum decrease in thickness for a very narrow ethanol volume composition range of 15% to 17%. These observations are in agreement with the recently published preferential adsorption model of the cononsolvency effect.

Phase Behavior of Co-Nonsolvent Systems: Poly( N-isopropylacrylamide) in Mixed Solvents of Water and Methanol

Cloud points of poly( N-isopropylacrylamide) in aqueous mixed solvents, with methanol as the cosolvent, are experimentally measured for polymer concentrations varied up to as high as the weight fraction 0.25. They are shown to form closed loops on the ternary phase plane in the temperature region between 5 and 30 °C, and hence co-nonsolvency is complete. Miscibility loops shrink by cooling, or equivalently, they exhibit lower critical solution temperature behavior. For a fixed polymer concentration, there is a composition of the mixed solvent at which the cloud-point temperature takes the lowest value. This minimum cloud-point temperature composition of the mixed solvent turned out to be almost independent of the polymer concentration, at least within the measured dilute region below the weight fraction 0.25. On the basis of the assumption that the phase separation is closely related to the preferential adsorption of the solvents by hydrogen bonding, we employ a model solution of Flory-Huggins type, augmented with direct and cooperative polymer-solvent hydrogen bonds, to construct the ternary phase diagrams. Theoretical calculation of the spinodal curves is performed, and the results are compared with the obtained experimental cloud-point data. The effect of molecular volume of the cosolvent is also studied within the same theoretical framework. Possibility for a upper critical solution temperature co-nonsolvency to appear for cosolvents with larger molecular volume is discussed.

Cosolvent Effects on Polymer Hydration Drive Hydrophobic Collapse

Water-mediated hydrophobic interactions play an important role in self-assembly processes, aqueous polymer solubility, and protein folding, to name a few. Cosolvents affect these interactions; however, the implications for hydrophobic polymer collapse and protein folding equilibria are not well-understood. This study examines cosolvent effects on the hydrophobic collapse equilibrium of a generic 32-mer hydrophobic polymer in urea, trimethylamine- N-oxide (TMAO), and acetone aqueous solutions using molecular dynamics simulations. Our results unveil a remarkable cosolvent-concentration-dependent behavior. Urea, TMAO, and acetone all shift the equilibrium toward collapsed structures below 2 M cosolvent concentration and, in turn, to unfolded structures at higher cosolvent concentrations, irrespective of the differences in cosolvent chemistry and the nature of cosolvent-water interactions. We find that weakly attractive polymer-water van der Waals interactions oppose polymer collapse in pure water, corroborating related observations reviewed by Ben-Amotz ( Annu. Rev. Phys. Chem. 2016, 67, 617-638). The cosolvents studied in the present work adsorb at the polymer/water interface and expel water molecules into the bulk, thereby effectively removing the dehydration energy penalty that opposes polymer collapse in pure water. At low cosolvent concentrations, this leads to cosolvent-induced stabilization of collapsed polymer structures. Only at sufficiently high cosolvent concentrations, polymer-cosolvent interactions favor polymer unfolding.

On the cononsolvency behaviour of hydrophobic clusters in water–methanol solutions

The solvent-excluded volume effect cannot rationalize the cononsolvency phenomenon of polymers, such as PNIPAM, in water–methanol solutions.