On the molecular origin of the cooperative coil-to-globule transition of poly(N-isopropylacrylamide) in water

Entropy Rules: Molecular Dynamics Simulations of Model Oligomers for Thermoresponsive Polymers

We attempted to attain atomic-scale insights into the mechanism of the heat-induced phase transition of two thermoresponsive polymers containing amide groups, poly(N-isopropylacrylamide) (PNIPAM) and poly(2-isopropyl-2-oxazoline) (PIPOZ), and we succeeded in reproducing the existence of lower critical solution temperature (LCST). The simulation data are in accord with experimental findings. We found out that the entropy has an important contribution to the thermodynamics of the phase separation transition. Moreover, after decomposing further the entropy change to contributions from the solutes and from the solvent, it appeared out that the entropy of the solvent has the decisive share for the lowering of the free energy of the system when increasing the temperature above the LCST. Our conclusion is that the thermoresponsive behavior is driven by the entropy of the solvent. The water molecules structured around the functional groups of the polymer that are exposed to contact with the solvent in the extended conformation lower the enthalpy of the system, but at certain temperature the extended conformation of the polymer collapses as a result of dominating entropy gain from “released” water molecules. We stress also on the importance of using more than one reference molecule in the simulation box at the setup of the simulation.

Thermosensitive Hydration of Four Acrylamide-Based Polymers in Coil and Globule Conformations

To characterize the thermosensitive coil–globule transition in atomistic detail, the conformational dynamics of linear polymer chains of acrylamide-based polymers have been investigated at multiple temperatures. Therefore, molecular dynamic simulations of 30mers of polyacrylamide (AAm), poly-N-methylacrylamide (NMAAm), poly-N-ethylacrylamide (NEAAm), and poly-N-isopropylacrylamide (NIPAAm) have been performed at temperatures ranging from 250 to 360 K for 2 μs. While two of the polymers are known to exhibit thermosensitivity (NEAAm, NIPAAm), no thermosensitivity is observed for AAm and NMAAm in aqueous solution. Our computer simulations consistently reproduce these properties. To understand the thermosensitivity of the respective polymers, the conformational ensembles at different temperatures have been separated according to the coil–globule transition. The coil and globule conformational ensembles were exhaustively analyzed in terms of hydrogen bonding with the solvent, the change of the solvent accessible surface, and enthalpic contributions. Surprisingly, independent of different thermosensitive properties of the four polymers, the surface affinity to water of coil conformations is higher than for globule conformations. Therefore, polymer–solvent interactions stabilize coil conformations at all temperatures. Nevertheless, the enthalpic contributions alone cannot explain the differences in thermosensitivity. This clearly implies that entropy is the distinctive factor for thermosensitivity. With increasing side chain length, the lifetime of the hydrogen bonds between the polymer surface and water is extended. Thus, we surmise that a longer side chain induces a larger entropic penalty due to immobilization of water molecules.

A molecular brush with thermoresponsive poly(2-ethyl-2-oxazoline) side chains: a structural investigation

The thermoresponsive behavior of a poly(2-oxazoline)-based molecular brush is investigated in aqueous solution. The molecular brush under study, PiPOx100-g-PEtOx17, has a poly(2-isopropenyl-2-oxazoline) (PiPOx) backbone grafted with thermoresponsive poly(2-ethyl-2-oxazoline) (PEtOx) side chains. Since the backbone degree of polymerization is only a factor of ~ 6 higher than the ones of the side chains, it features an architecture between a star-like polymer and a comb-like polymer. Its aqueous solution exhibits lower critical solution temperature (LCST) behavior with a cloud point temperature Tcp = 40.5 °C at 30 g L−1. The temperature-dependent structural evolution is disclosed using dynamic light scattering (DLS) and small-angle neutron scattering (SANS). An increase of the molecular brush size is found upon heating from room temperature to Tcp, which is attributed to the extension of the backbone resulting from the dehydration and collapse of the side chains. Above Tcp, the size decreases again, which indicates the collapse of the whole molecular brush. Large aggregates are found to be present in the solution in the temperature range 25–50 °C. These become more compact, as the temperature is increased across Tcp.

Robust Multiscale-Oriented Thermoresponsive Fibrous Hydrogels with Rapid Self-Recovery and Ultrafast Response Underwater

Hydrogels with ultrafast response to environmental stimuli, possessing robust structural integrity and rapid self-recovery, have been considered as promising platforms for numerous applications, for example, in biomimetic materials and nanomedicine. Inspired by the bundled fibrous structure of actin, we developed a robust and ultrafast thermoresponsive fibrous hydrogel (TFH) by fully utilizing the weak noncovalent bonds and strong covalently cross-linked semiflexible electrospun fibrous nets. The TFH exhibits an ultrafast response (within 10 s), rapid self-recovery rate (74% within 10 s), tunable tensile strength (3-380 kPa), and high toughness (∼1560 J/m²) toward temperature. A multiscale orientation is considered to play a key role in the excellent mechanical properties at the fibrous mesh, fiber, and molecular scales. Furthermore, to take advantage of this TFH adequately, a novel kind of noodle-like hydrogel for thermo-controlled protein sorption based on the TFH is prepared, which exhibits high stability and ultrafast sorption properties. The bioinspired platforms hold promise as artificial skins and "smart" sorption membrane carriers, which provide a unique bioactive environment for tissue engineering and nanomedicine.

Atomic scale investigation of the volume phase transition in concentrated PNIPAM microgels

Combining elastic incoherent neutron scattering and differential scanning calorimetry, we investigate the occurrence of the volume phase transition (VPT) in very concentrated poly-(N-isopropyl-acrylamide) (PNIPAM) microgel suspensions, from a polymer weight fraction of 30 wt. % up to dry conditions. Although samples are arrested at the macroscopic scale, atomic degrees of freedom are equilibrated and can be probed in a reproducible way. A clear signature of the VPT is present as a sharp drop in the mean square displacement of PNIPAM hydrogen atoms obtained by neutron scattering. As a function of concentration, the VPT gets smoother as dry conditions are approached, whereas the VPT temperature shows a minimum at about 43 wt. %. This behavior is qualitatively confirmed by calorimetry measurements. Molecular dynamics simulations are employed to complement experimental results and gain further insights into the nature of the VPT, confirming that it involves the formation of an attractive gel state between the microgels. Overall, these results provide evidence that the VPT in PNIPAM-based systems can be detected at different time- and length-scales as well as under overcrowded conditions.

Influence of crowding on hydrophobic hydration-shell structure

The influence of molecular crowding on water structure, and the associated crossover behavior, is quantified using Raman multivariate curve resolution (Raman-MCR) hydration-shell vibrational spectroscopy of aqueous tert-butyl alcohol, 2-butyl alcohol and 2-butoxyethanol solutions of variable concentration and temperature. Changes in the hydration-shell OH stretch band shape and mean frequency are used to identify the temperature at which the hydration-shell crosses over from a more ordered to less ordered structure, relative to pure water. The influence of crowding on the crossover is found to depend on solute size and shape in a way that is correlated with the corresponding infinitely dilute hydration-shell structure (and the corresponding first hydration-shell spectra are invariably very similar to pure water). Analysis of the results using a Muller-like two-state equilibrium between more ordered and less ordered hydration-shell structures implies that crossover temperature changes are dictated primarily by enthalpic stabilization of the more ordered hydration-shell structures.

Accurate and rapid calculation of hydration free energy and its physical implication for biomolecular functions

Here we review a new method for calculating a hydration free energy (HFE) of a solute and discuss its physical implication for biomolecular functions in aqueous environments. The solute hydration is decomposed into processes 1 and 2. A cavity matching the geometric characteristics of the solute at the atomic level is created in process 1. Solute-water van der Waals and electrostatic interaction potentials are incorporated in process 2. The angle-dependent integral equation theory combined with our morphometric approach is applied to process 1, and the three-dimensional reference interaction site model theory is employed for process 2. Molecular models are adopted for water. The new method is characterized by the following. Solutes with various sizes including proteins can be treated in the same manner. It is almost as accurate as the molecular dynamics simulation despite its far smaller computational burden. It enables us to handle a solute possessing a significantly large total charge without difficulty. The HFE can be decomposed into a variety of physically insightful, energetic, and entropic components. It is best suited to the elucidation of mechanisms of protein folding, pressure and cold denaturation of a protein, and different types of molecular recognition.

Solvation dynamics of N-substituted acrylamide polymers and the importance for phase transition behavior

Functional groups present in thermo-responsive polymers are known to play an important role in aqueous solutions by manifesting their coil-to-globule conformational transition in a specific temperature range. Understanding the role of these functional groups and their interactions with water is of great interest as it may allow us to control both the nature and temperature of this coil-to-globule transition. In this work, polyacrylamide (PAAm), poly(N-isopropylacrylamide) (PNIPAm), and poly(N-isopropylmethacrylamide) (PNIPMAm) solvated in water are studied with the goal of discovering the structure of the solvent and its interaction with these polymers in determining the polymer conformations. Specifically, all-atom molecular dynamics (MD) simulations were performed on polymer chains with 30 monomer units (30-mers) at 295 K, 310 K and 320 K, which is below and above the lower critical solution temperature (LCST) of PNIPAm (LCST = 305 K) and PNIPMAm (LCST = 315 K), respectively. The MD simulation trajectories suggest that changes in the functional groups in the backbone and side-chains alter the water solvation shell around the polymer. This results in a change in the residence time probability and hydrogen bond characteristics of water at simulated temperatures. Specifically, water molecules reside for longer times near PAAm (no LCST) and PNIPMAm (LCST = 315 K) chains as compared to PNIPAm. This might be one of the possible causes for the higher LCST of PNIPMAm as compared to that of PNIPAm. These results can guide experimentalists and theoreticians to design new polymer structures with tailor-made LCST transitions while controlling the water solvation shell around the functional group.

Reduction of water-mediated repulsion drives poly(N-vinylcaprolactam) collapse upon heating

Thermo-sensitive aqueous polymers undergo a coil-to-globule transition on heating, with drastic chemical and structural changes. We performed molecular dynamics simulations for PVCL in water to study the driving forces for the polymer's collapse.

Fractography of poly(N-isopropylacrylamide) hydrogel networks crosslinked with mechanofluorophores using confocal laser scanning microscopy

Employing mechanofluorophores in polymer fractography to obtain new information on force-induced events when analyzed by confocal laser scanning microscopy.

A coil-to-globule transition capable coarse-grained model for poly(N-isopropylacrylamide)

We present a model for mesoscopic molecular dynamics simulations of poly(N-isopropyl-acrylamide) (pNIPAM).

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.

Functional Surfaces through Controlled Assemblies of Upper Critical Solution Temperature Block and Star Copolymers

Endowing surfaces with multiple advanced functionalities, such as temperature-controlled swelling or the triggered release of functional small molecules, is attractive for a large variety of applications ranging from smart textiles to advanced biomedical applications. This Invited Feature Article summarizes recent advances in the development of upper critical solution temperature (UCST) behavior of copolymers in aqueous solutions and compares the fundamental differences between lower critical solution temperature (LCST) and UCST transitions. The effect of polymer chemistry and architecture on UCST transitions is discussed for block copolymer micelles (BCMs) and star polymers in solution and assembled at surfaces. The inclusion of such nanocontainers (i.e., BCMs and star polymers) in layer-by-layer (LbL) coatings and how to control their responsive behavior through deposition conditions and binding partners is explored. Finally, the inclusion and temperature-triggered release of functional small molecules is explored for nanocontainers in LbL coatings. Taken together, UCST nanocontainers containing LbL films are promising building blocks for the development of new generations of practical, functional surface coatings.

Preparation of thermo-responsive PNIPAAm-MWCNT membranes and evaluation of its antifouling properties in dairy wastewater

A novel MWCNT-PNIPAAm nanocomposite membrane was developed with an excellent cleaning efficiency of thermo-responsive surface. The thermo-responsive N-isopropyle acryleamide (NIPAAm) monomer was polymerized on the surface of MWCNT via free radical polymerization. The prepared MWCNT-PNIPAAm nanocomposite was characterized by FTIR, SEM and TGA analyses. Various amounts of the prepared nanocomposite were incorporated into the membrane matrix by the physical blending method. The resultant membranes showed better surface wettability and pure water flux compared to pristine Polyethersulfone (PES) membrane. Furthermore, after filtration, the COD value of dairy wastewater was reduced to around 90% for all membranes. The thermo-responsive cleaning method was employed to investigate the cleaning efficiency of MWCNT-PNIPAAm membrane for dairy wastewater. The 99.9% flux recovery ratio was obtained for MWCNT-PNIPAAm-0.05% membranes. All these results confirmed that the presence of MWCNT-PNIPAAm nanocomposite in the membrane matrix improves the membrane hydrophilicity and antifouling properties.

Plant leaves inspired sunlight-driven purifier for high-efficiency clean water production

Natural vascular plants leaves rely on differences in osmotic pressure, transpiration and guttation to produce tons of clean water, powered by sunlight. Inspired by this, we report a sunlight-driven purifier for high-efficiency water purification and production. This sunlight-driven purifier is characterized by a negative temperature response poly(N-isopropylacrylamide) hydrogel (PN) anchored onto a superhydrophilic melamine foam skeleton, and a layer of PNIPAm modified graphene (PG) filter membrane coated outside. Molecular dynamics simulation and experimental results show that the superhydrophilicity of the relatively rigid melamine skeleton significantly accelerates the swelling/deswelling rate of the PNPG-F purifier. Under one sun, this rational engineered structure offers a collection of 4.2 kg m⁻² h⁻¹ and an ionic rejection of > 99% for a single PNPG-F from brine feed via the cooperation of transpiration and guttation. We envision that such a high-efficiency sunlight driven system could have great potential applications in diverse water treatments.

Natural leaves can purify water under sunlight through a combination of osmotic pressure, transpiration, and guttation effects. Here the authors design a composite material mimicking these combined effects, achieving sunlight-driven pure water production from brine with high collection rate.

Understanding the role of water in the aggregation of poly(N,N-dimethylaminoethyl methacrylate) in aqueous solution using temperature-dependent near-infrared spectroscopy

For understanding the role of water in the aggregation of polymers, the variation of water structures with the structural change of polymers in the process of aggregation was studied by temperature-dependent near-infrared (NIR) spectroscopy. The NIR spectra of the aqueous poly(N,N-dimethylaminoethyl methacrylate) (PDMAEMA) solutions of different concentrations were measured at different temperatures. The spectral changes of the polymer and water with temperature were analyzed by N-way principal component analysis (NPCA). It was found that, at low concentration, the chains of the polymer tend to form a loose hydrophobic structure below 36 °C and then aggregate into a micelle at a lower critical solution temperature (LCST) of around 39 °C. In the process of the aggregation, the water species with two hydrogen bonds (S2) increases gradually before 36 °C and then a sudden decrease occurs after that temperature. The results clearly indicate that water species S2 plays an important role in the formation of the intermediate, i.e., the loose hydrophobic structure of the polymer chains linked by the two hydrogen bonds of S2 water. When the temperature increases, the dissociation of the hydrogen bonds enables the intermediate to be destroyed to form a micelle structure. For the high concentration solution, however, the spectral information of S2 was not found in the aggregation, suggesting direct formation of the micelle from the dehydrated chains.

Water-Polymer Coupling Induces a Dynamical Transition in Microgels

The long debated protein dynamical transition was recently found also in nonbiological macromolecules, such as poly- N-isopropylacrylamide (PNIPAM) microgels. Here, by using atomistic molecular dynamics simulations, we report a description of the molecular origin of the dynamical transition in these systems. We show that PNIPAM and water dynamics below the dynamical transition temperature T _d are dominated by methyl group rotations and hydrogen bonding, respectively. By comparing with bulk water, we unambiguously identify PNIPAM-water hydrogen bonding as mainly responsible for the occurrence of the transition. The observed phenomenology thus crucially depends on the water-macromolecule coupling, being relevant to a wide class of hydrated systems, independently from the biological function.

Numerical modelling of non-ionic microgels: an overview

Microgels are complex macromolecules. These colloid-sized polymer networks possess internal degrees of freedom and, depending on the polymer(s) they are made of, can acquire a responsiveness to variations of the environment (temperature, pH, salt concentration, etc.). Besides being valuable for many practical applications, microgels are also extremely important to tackle fundamental physics problems. As a result, these last years have seen a rapid development of protocols for the synthesis of microgels, and more and more research has been devoted to the investigation of their bulk properties. However, from a numerical standpoint the picture is more fragmented, as the inherently multi-scale nature of microgels, whose bulk behaviour crucially depends on the microscopic details, cannot be handled at a single level of coarse-graining. Here we present an overview of the methods and models that have been proposed to describe non-ionic microgels at different length-scales, from the atomistic to the single-particle level. We especially focus on monomer-resolved models, as these have the right level of details to capture the most important properties of microgels, responsiveness and softness. We suggest that these microscopic descriptions, if realistic enough, can be employed as starting points to develop the more coarse-grained representations required to investigate the behaviour of bulk suspensions.

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.

Molecular dynamics study of coil-to-globule transition in a thermo-responsive oligomer bound to various surfaces: hydrophilic surfaces stabilize the coil form

The structural and dynamical properties of 40-mer of thermo-responsive polymer PNIPAM covalently bound to different surfaces have been studied, at different temperatures, by means of molecular dynamics simulations. Evolution of the radius of gyration, Rg, of the polymer chain and radial distribution functions (RDFs) calculated for the carbon atoms of the PNIPAM backbone with water oxygens and for the hydrogen atom of the amide groups with water oxygens indicate that functionalized surfaces affect the coil-to-globule transition of PNIPAM, by means of electrostatic interactions, increasing the lower critical solution temperature (LCST) of the polymer. Such interactions, mainly represented by a H-bond, hinder the transition in the globular form while hydrophobic groups on the surface, such as -OCH3, contribute to the globular collapse. A significant alteration in the arrangement of water molecules around the polymer is testified by: (i) the absence of the second peak in the RDF between the C atoms of the PNIPAM backbone and the O atoms of water at the same temperature at which the radius of gyration decreases; (ii) the height of both the first and the second peak of the RDF between the H atom of the amide groups and water O atoms decreases when the temperature increases above the LCST. Finally, the H-bond autocorrelation function indicates that: (i) hydrogen bonds between the bound-to-surface PNIPAM acceptor groups (O[double bond, length as m-dash]C[double bond splayed right]) and the H atoms of water molecules are less persistent than H-bonds formed between the free PNIPAM acceptor groups and water; (ii) H-bonds between the PNIPAM acceptor groups and hydroxyl groups on the quartz surface are longer lived than those formed on graphene oxide.

Modelling realistic microgels in an explicit solvent

Thermoresponsive microgels are polymeric colloidal networks that can change their size in response to a temperature variation. This peculiar feature is driven by the nature of the solvent-polymer interactions, which triggers the so-called volume phase transition from a swollen to a collapsed state above a characteristic temperature. Recently, an advanced modelling protocol to assemble realistic, disordered microgels has been shown to reproduce experimental swelling behavior and form factors. In the original framework, the solvent was taken into account in an implicit way, condensing solvent-polymer interactions in an effective attraction between monomers. To go one step further, in this work we perform simulations of realistic microgels in an explicit solvent. We identify a suitable model which fully captures the main features of the implicit model and further provides information on the solvent uptake by the interior of the microgel network and on its role in the collapse kinetics. These results pave the way for addressing problems where solvent effects are dominant, such as the case of microgels at liquid-liquid interfaces.