Conformation change of an isotactic poly (N-isopropylacrylamide) membrane: Molecular dynamics

Effect of Tacticity Sequence of the Poly(N-isopropylacrylamide) Oligomer on Phase Transition Behavior in Aqueous Solution

The tacticity of poly(N-isopropylacrylamide) (PNIPAM) has a strong impact on the lower critical solution temperature (LCST) in aqueous solution. The sequence of meso diads (m) and racemo diads (r) further contributes to such an effect. In this work, the phase transition behaviors of poly(N-isopropylacrylamide) pentamers with four kinds of sequences, i.e., rrmm, rmmr, mrrm, and rmrm, in water were studied applying replica exchange molecular dynamics with a modified OPLS/AA force field. The difference in local component concentration in the system was used as an order parameter to quantitatively describe the phase separation extent. It was found that the phase separation degree of rrmm and rmmr is higher than that of mrrm and rmrm at the same temperature. The LCSTs of rrmm and rmmr are lower than those of mrrm and rmrm. The radial distribution function and hydrogen bond analysis revealed that the average values of hydrogen bonds between pentamers for rrmm and rmmr are greater than those of mrrm and rmrm, whereas the average values of hydrogen bonds between pentamers and water for rrmm and rmmr are less than those of mrrm and rmrm. It was demonstrated that the isotactic triad (mm) plays an important role in the thermosensitive behaviors of the PNIPAM pentamer. The increase of isotactic triad (mm) content in the PNIPAM chain promotes the formation of intermolecular hydrogen bonds between amide and amide and leads to a higher aggregation of the pentamer with the sequence of rrmm or rmmr. Finally, the effect of the isotactic triad was qualitatively explained with the mean-field theory.

Phenol release from pNIPAM hydrogels: scaling molecular dynamics simulations with dynamical density functional theory

We employed molecular dynamic simulations (MD) and the Bennett's acceptance ratio method to compute the free energy of transfer, ΔG_trans, of phenol, methane, and 5-fluorouracil (5-FU), between bulk water and water-pNIPAM mixtures of different polymer volume fractions, ϕ_p. For this purpose, we first calculate the solvation free energies in both media to obtain ΔG_trans. Phenol and 5-FU (a medication used to treat cancer) attach to the pNIPAM surface so that they show negative values of ΔG_trans irrespective of temperature (above or below the lower critical solution temperature of pNIPAM, T_c). Conversely, methane switches the ΔG_trans sign when considering temperatures below (positive) and above (negative) T_c. In all cases, and contrasting with some theoretical predictions, ΔG_trans maintains a linear behavior with the pNIPAM concentration up to large polymer densities. We have also employed MD to compute the diffusion coefficient, D, of phenol in water-pNIPAM mixtures as a function of ϕ_p in the diluted limit. Both ΔG_trans and D as a function of ϕ_p are required inputs to obtain the release halftime of hollow pNIPAM microgels through Dynamic Density Functional Theory (DDFT). Our scaling strategy captures the experimental value of 2200 s for 50 μm radius microgels with no cavity, for ϕ_p ≃ 0.83 at 315 K.

Modeling Solution Behavior of Poly(N-isopropylacrylamide): A Comparison between Water Models

Water is known to play a fundamental role in determining the structure and functionality of macromolecules. The same crucial contribution is also found in the in silico description of polymer aqueous solutions. In this work, we exploit the widely investigated synthetic polymer poly(N-isopropylacrylamide) (PNIPAM) to understand the effect of the adopted water model on its solution behavior and to refine the computational setup. By means of atomistic molecular dynamics simulations, we perform a comparative study of PNIPAM aqueous solution using two advanced water models: TIP4P/2005 and TIP4P/Ice. The conformation and hydration features of an atactic 30-mer at infinite dilution are probed at a range of temperature and pressure suitable to detect the coil-to-globule transition and to map the P–T phase diagram. Although both water models can reproduce the temperature-induced coil-to-globule transition at atmospheric pressure and the polymer hydration enhancement that occurs with increasing pressure, the PNIPAM–TIP4P/Ice solution shows better agreement with experimental findings. This result can be attributed to a stronger interaction of TIP4P/Ice water with both hydrophilic and hydrophobic groups of PNIPAM, as well as to a less favorable contribution of the solvent entropy to the coil-to-globule transition.

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.

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.

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).

Aqueous Nanoclusters Govern Ion Partitioning in Dense Polymer Membranes

The uptake and sorption of charged molecules by responsive polymer membranes and hydrogels in aqueous solutions is of key importance for the development of soft functional materials. Here, we investigate the partitioning of simple monatomic (Na⁺, K⁺, Cs⁺, Cl^-, I^-) and one molecular ion (4-nitrophenolate; NP^-) within a dense, electroneutral poly(N-isopropylacrylamide) membrane using explicit-water molecular dynamics simulations. Inside the predominantly hydrophobic environment, water distributes in a network of polydisperse water nanoclusters. The average cluster size determines the mean electrostatic self-energy of the simple ions, which preferably reside deeply inside them; we therefore find substantially larger partition ratios K ≃10^-1 than expected from a simple Born picture using a uniform dielectric constant. Despite their irregular shapes, we observe that the water clusters possess a universal negative electrostatic potential with respect to their surroundings, as is known for aqueous liquid-vapor interfaces. This potential, which we find concealed in cases of symmetric monatomic salts, can dramatically impact the transfer free energies of larger charged molecules because of their weak hydration and increased affinity to interfaces. Consequently, and in stark contrast to the simple ions, the molecular ion NP^- can have a partition ratio much larger than unity, K ≃10-30 (depending on the cation type) or even 10³ in excess of monovalent salt, which explains recent observations of enhanced reaction kinetics of NP^- reduction catalyzed within dense polymer networks. These results also suggest that ionizing a molecule can even enhance the partitioning in a collapsed, rather hydrophobic gel, which strongly challenges the traditional simplistic reasoning.

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.

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.

Transfer Free Energies and Partitioning of Small Molecules in Collapsed PNIPAM Polymers

A central quantity in the design of functional hydrogels used as nanocarrier systems, for instance, for drug delivery or adaptive nanocatalysis, is the partition ratio, which quantifies the uptake of a molecular substance by the polymer matrix. By employing all-atom molecular dynamics simulations, we study the solvation and partitioning (with respect to bulk water) of small subnanometer-sized solutes in a dense matrix of collapsed poly( N-isopropylacrylamide) (PNIPAM) polymers above the lower critical solution temperature in aqueous solution. We examine the roles of the solute's polarity and its size on the solubility properties in the thermoresponsive polymer. We show that the transfer free energies of nonpolar solutes from bulk water into the polymer are favorable and scale in a good approximation with the solute's surface area. Even for small solute size variation, partitioning can vary over orders of magnitude. A polar nature of the solute, on the other hand, generally opposes the transfer, at least for alkyl solutes. Finally, we find a strong correlation between the transfer free energies in the gel and the adsorption free energies on a single extended polymer chain, which enables us to relate the partition ratios in the swollen and collapsed state of a PNIPAM gel.

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.

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

The cooperativity of PNIPAM coil-to-globule transition in water arises from the structuring of solvent in proximity to hydrophobic groups.

By means of atomistic molecular dynamics simulations we investigate the behaviour of poly(N-isopropylacrylamide), PNIPAM, in water at temperatures below and above the lower critical solution temperature (LCST), including the undercooled regime. The transition between water soluble and insoluble states at the LCST is described as a cooperative process involving an intramolecular coil-to-globule transition preceding the aggregation of chains and the polymer precipitation. In this work we investigate the molecular origin of such cooperativity and the evolution of the hydration pattern in the undercooled polymer solution. The solution behaviour of an atactic 30-mer at high dilution is studied in the temperature interval from 243 to 323 K with a favourable comparison to available experimental data. In the water soluble states of PNIPAM we detect a correlation between polymer segmental dynamics and diffusion motion of bound water, occurring with the same activation energy. Simulation results show that below the coil-to-globule transition temperature PNIPAM is surrounded by a network of hydrogen bonded water molecules and that the cooperativity arises from the structuring of water clusters in proximity to hydrophobic groups. Differently, the perturbation of the hydrogen bond pattern involving water and amide groups occurs above the transition temperature. Altogether these findings reveal that even above the LCST PNIPAM remains largely hydrated and that the coil-to-globule transition is related with a significant rearrangement of the solvent in the proximity of the surface of the polymer. The comparison between the hydrogen bonding of water in the surrounding of PNIPAM isopropyl groups and in the bulk displays a decreased structuring of solvent at the hydrophobic polymer–water interface across the transition temperature, as expected because of the topological extension along the chain of such interface. No evidence of an upper critical solution temperature behaviour, postulated in theoretical and thermodynamics studies of PNIPAM aqueous solution, is observed in the low temperature domain.

Development and Application of a Coarse-Grained Model for PNIPAM by Iterative Boltzmann Inversion and Its Combination with Lattice Boltzmann Hydrodynamics

The polymer poly(N-isopropylacrylamide) (PNIPAM) is studied using a novel combination of multiscale modeling methodologies. We develop an iterative Boltzmann inversion potential of concentrated PNIPAM solutions and combine it with lattice Boltzmann as a Navier-Stokes equation solver for the solvent. We study in detail the influence of the methodology on statics and dynamics of the system. The combination is successful and significantly simpler and faster than other mapping techniques for polymer solution while keeping the correct hydrodynamics. The model can semiquantitatively describe the correct phase behavior and polymer dynamics.