Selective Molecular Transport in Thermoresponsive Polymer Membranes: Role of Nanoscale Hydration and Fluctuations

ROS-responsive hydrogels: from design and additive manufacturing to biomedical applications

Hydrogels with intricate 3D networks and high hydrophilicity have qualities resembling those of biological tissues, making them ideal candidates for use as smart biomedical materials. Reactive oxygen species (ROS) responsive hydrogels are an innovative class of smart hydrogels, and are cross-linked by ROS-responsive modules through covalent interactions, coordination interactions, or supramolecular interactions. Due to the introduction of ROS response modules, this class of hydrogels exhibits a sensitive response to the oxidative stress microenvironment existing in organisms. Simultaneously, due to the modularity of the ROS-responsive structure, ROS-responsive hydrogels can be manufactured on a large scale through additive manufacturing. This review will delve into the design, fabrication, and applications of ROS-responsive hydrogels. The main goal is to clarify the chemical principles that govern the response mechanism of these hydrogels, further providing new perspectives and methods for designing responsive hydrogel materials.

Feedback-controlled solute transport through chemo-responsive polymer membranes

Polymer membranes are typically assumed to be inert and nonresponsive to the flux and density of the permeating particles in transport processes. Here, we theoretically study the consequences of membrane responsiveness and feedback on the steady-state force-flux relations and membrane permeability using a nonlinear-feedback solution-diffusion model of transport through a slab-like membrane. Therein, the solute concentration inside the membrane depends on the bulk concentration, c₀, the driving force, f, and the polymer volume fraction, ϕ. In our model, the solute accumulation in the membrane causes a sigmoidal volume phase transition of the polymer, changing its permeability, which, in return, affects the membrane's solute uptake. This feedback leads to nonlinear force-flux relations, j(f), which we quantify in terms of the system's differential permeability, P_sys ^Δ∝dj/df. We find that the membrane feedback can increase or decrease the solute flux by orders of magnitude, triggered by a small change in the driving force and largely tunable by attractive vs repulsive solute-membrane interactions. Moreover, controlling the inputs, c₀ and f, can lead to the steady-state bistability of ϕ and hysteresis in the force-flux relations. This work advocates that the fine-tuning of the membrane's chemo-responsiveness will enhance the nonlinear transport control features, providing great potential for future (self-)regulating membrane devices.

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.

Nonequilibrium Processes in Polymer Membrane Formation: Theory and Experiment

Porous polymer and copolymer membranes are useful for ultrafiltration of functional macromolecules, colloids, and water purification. In particular, block copolymer membranes offer a bottom-up approach to form isoporous membranes. To optimize permeability, selectivity, longevity, and cost, and to rationally design fabrication processes, direct insights into the spatiotemporal structure evolution are necessary. Because of a multitude of nonequilibrium processes in polymer membrane formation, theoretical predictions via continuum models and particle simulations remain a challenge. We compiled experimental observations and theoretical approaches for homo- and block copolymer membranes prepared by nonsolvent-induced phase separation and highlight the interplay of multiple nonequilibrium processes─evaporation, solvent-nonsolvent exchange, diffusion, hydrodynamic flow, viscoelasticity, macro- and microphase separation, and dynamic arrest─that dictates the complex structure of the membrane on different scales.

Theory of the effect of external stress on the activated dynamics and transport of dilute penetrants in supercooled liquids and glasses

We generalize the self-consistent cooperative hopping theory for a dilute spherical penetrant or tracer activated dynamics in dense metastable hard sphere fluids and glasses to address the effect of external stress, the consequences of which are systematically established as a function of matrix packing fraction and penetrant-to-matrix size ratio. All relaxation processes speed up under stress, but the difference between the penetrant and matrix hopping (alpha relaxation) times decreases significantly with stress corresponding to less time scale decoupling. A dynamic crossover occurs at a critical "slaving onset" stress beyond which the matrix activated hopping relaxation time controls the penetrant hopping time. This characteristic stress increases (decreases) exponentially with packing fraction (size ratio) and can be well below the absolute yield stress of the matrix. Below the slaving onset, the penetrant hopping time is predicted to vary exponentially with stress, differing from the power law dependence of the pure matrix alpha time due to system-specificity of the stress-induced changes in the penetrant local cage and elastic barriers. An exponential growth of the penetrant alpha relaxation time with size ratio under stress is predicted, and at a fixed matrix packing fraction, the exponential relation between penetrant hopping time and stress for different size ratios can be collapsed onto a master curve. Direct connections between the short- and long-time activated penetrant dynamics and between the penetrant (or matrix) alpha relaxation time and matrix thermodynamic dimensionless compressibility are also predicted. The presented results should be testable in future experiments and simulations.

Tuning the permeability of regular polymeric networks by the cross-link ratio

The amount of cross-linking in the design of polymer materials is a key parameter for the modification of numerous physical properties, importantly, the permeability to molecular solutes. We consider networks with a diamond-like architecture and different cross-link ratios, concurring with a wide range of the polymer volume fraction. We particularly focus on the effect and the competition of two independent component-specific solute-polymer interactions, i.e., we distinguish between chain-monomers and cross-linkers, which individually act on the solutes and are altered to cover attractive and repulsive regimes. For this purpose, we employ coarse-grained, Langevin computer simulations to study how the cross-link ratio of polymer networks controls the solute partitioning, diffusion, and permeability. We observe different qualitative behaviors as a function of the cross-link ratio and interaction strengths. The permeability can be tuned ranging over two orders of magnitude relative to the reference bulk permeability. Finally, we provide scaling theories for the partitioning and diffusion that explicitly account for the component-specific interactions as well as the cross-link ratio and the polymer volume fraction. These are in overall good agreement with the simulation results and grant insight into the underlying physics, rationalizing how the cross-link ratio can be exploited to tune the solute permeability of polymeric networks.

Activated penetrant dynamics in glass forming liquids: size effects, decoupling, slaving, collective elasticity and correlation with matrix compressibility

We employ the microscopic self-consistent cooperative hopping theory of penetrant activated dynamics in glass forming viscous liquids and colloidal suspensions to address new questions over a wide range of high matrix packing fractions and penetrant-to-matrix particle size ratios. The focus is on the mean activated relaxation time of smaller tracers in a hard sphere fluid of larger particle matrices. This quantity also determines the penetrant diffusion constant and connects directly with the structural relaxation time probed in an incoherent dynamic structure factor measurement. The timescale of the non-activated fast dissipative process is also studied and is predicted to follow power laws with the contact value of the penetrant-matrix pair correlation function and the penetrant-matrix size ratio. For long time penetrant relaxation, in the relatively lower packing fraction metastable regime the local cage barriers are dominant and matrix collective elasticity effects unimportant. As packing fraction and/or penetrant size grows, much higher barriers emerge and the collective elasticity associated with the correlated matrix dynamic displacement that facilitates penetrant hopping becomes important. This results in a non-monotonic variation with packing fraction of the degree of decoupling between the matrix and penetrant alpha relaxation times. The conditions required for penetrant hopping to become slaved to the matrix alpha process are determined, which depend mainly on the penetrant to matrix particle size ratio. By analyzing the absolute and relative importance of the cage and elastic barriers we establish a mechanistic understanding of the origin of the predicted exponential growth of the penetrant hopping time with size ratio predicted at very high packing fractions. A dynamics-thermodynamics power law connection between the penetrant activation barrier and the matrix dimensionless compressibility is established as a prediction of theory, with different scaling exponents depending on whether matrix collective elasticity effects are important. Quantitative comparisons with simulations of the penetrant relaxation time, diffusion constant, and transient localization length of tracers in dense colloidal suspensions and cold viscous liquids reveal good agreements. Multiple new predictions are made that are testable via future experiments and simulations. Extension of the theoretical approach to more complex systems of high experimental interest (nonspherical molecules, semiflexible polymers, crosslinked networks) interacting via variable hard or soft repulsions and/or short range attractions is possible, including under external deformation.

How the Shape and Chemistry of Molecular Penetrants Control Responsive Hydrogel Permeability

The permeability of hydrogels for the selective transport of molecular penetrants (drugs, toxins, reactants, etc.) is a central property in the design of soft functional materials, for instance in biomedical, pharmaceutical, and nanocatalysis applications. However, the permeation of dense and hydrated polymer membranes is a complex multifaceted molecular-level phenomenon, and our understanding of the underlying physicochemical principles is still very limited. Here, we uncover the molecular principles of permeability and selectivity in hydrogel permeation. We combine the solution-diffusion model for permeability with comprehensive atomistic simulations of molecules of various shapes and polarities in a responsive hydrogel in different hydration states. We find in particular that dense collapsed states are extremely selective, owing to a delicate balance between the partitioning and diffusivity of the penetrants. These properties are sensitively tuned by the penetrant size, shape, and chemistry, leading to vast cancellation effects, which nontrivially contribute to the permeability. The gained insights enable us to formulate semiempirical rules to quantify and extrapolate the permeability categorized by classes of molecules. They can be used as approximate guiding ("rule-of-thumb") principles to optimize penetrant or membrane physicochemical properties for a desired permeability and membrane functionality.

Tuning the selective permeability of polydisperse polymer networks

We study the permeability and selectivity ('permselectivity') of model membranes made of polydisperse polymer networks for molecular penetrant transport, using coarse-grained, implicit-solvent computer simulations. In our work, permeability P is determined on the linear-response level using the solution-diffusion model, P = KDin, i.e., by calculating the equilibrium penetrant partition ratio K and penetrant diffusivity Din inside the membrane. We vary two key parameters, namely the network-network interaction, which controls the degree of swelling and collapse of the network, and the network-penetrant interaction, which tunes the selective penetrant uptake and microscopic energy landscape for diffusive transport. We find that the partitioning K covers four orders of magnitude and is a non-monotonic function of the parameters, well interpreted by a second-order virial expansion of the free energy of transferring one penetrant from a reservoir into the membrane. Moreover, we find that the penetrant diffusivity Din in the polydisperse networks, in contrast to highly ordered membrane structures, exhibits relatively simple exponential decays. We propose a semi-empirical scaling law for the penetrant diffusion that describes the simulation data for a wide range of densities and interaction parameters. The resulting permeability P turns out to follow the qualitative behavior (including maximization and minimization) of partitioning. However, partitioning and diffusion are typically anti-correlated, yielding large quantitative cancellations, controlled and fine-tuned by the network density and interactions, as rationalized by our scaling laws. We finally demonstrate that even small changes of network-penetrant interactions, e.g., by half a kBT, modify the permselectivity by almost one order of magnitude.

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.

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.

Atomistic molecular dynamics simulations of the LCST conformational transition in poly(N-vinylcaprolactam) in water

Thermoresponsive poly(N-vinylcaprolactam) (PVCL) has received growing interest due to a temperature-induced phase transition, which switches its solubility in aqueous solutions. However, the lower critical solution temperature (LCST) of PVCL is greatly influenced by the molecular weight, morphology and the environment. Therefore, despite of numerous experimental studies of the thermal response of PVCL, a driving force and a molecular origin of conformation transitions in solution remain far less studied. To get a better understanding of the coil-to-globule conformation transition of PVCL in aqueous solution, we examined the structure and conformation dynamics of a single-chain PVCL₃₀ in a temperature range of 280-360 K by using atomistic molecular dynamics (MD) simulations. The united-atom GROMOS G53a6 force field was re-parameterized and fine-tuned by DFT calculations to reproduce the experimental LCST transition of PVCL. Our MD model reproduces the LCST transition of PVCL₃₀ to occur within a temperature range of 34.6-38.5°. MD simulation results suggest a significant difference between the hydration state of the carbonyl group of PVCL below and above the LCST threshold. The analysis of the number of hydrogen bonds of PVCL with water molecules demonstrates that dehydration of the polymer plays an important role and drives the temperature-induced polymer collapse. Finally, the developed MD model and FF parameters were successfully tested for large-scale systems, such as mixture PVCL₃₀ oligomer and single-chain PVCL₈₁₆ polymer, respectively.

Tuning the Permeability of Dense Membranes by Shaping Nanoscale Potentials

Permeability is one of the most fundamental transport properties in soft matter physics, material engineering, and nanofluidics. Here, we report by means of Langevin simulations of ideal penetrants in a nanoscale membrane made of a fixed lattice of attractive interaction sites, how the permeability can be massively tuned, even minimized or maximized, by tailoring the potential energy landscape for the diffusing penetrants, depending on the membrane attraction, topology, and density. Supported by limiting scaling theories we demonstrate that the observed nonmonotonic behavior and the occurrence of extreme values of the permeability is far from trivial and triggered by a strong anticorrelation and substantial (orders of magnitude) cancellation between penetrant partitioning and diffusivity, especially within dense and highly attractive membranes.

Nonequilibrium Uptake Kinetics of Molecular Cargo into Hollow Hydrogels Tuned by Electrosteric Interactions

Hollow hydrogels represent excellent nano- and microcarriers due to their ability to encapsulate and release large amounts of cargo molecules (cosolutes) such as reactants, drugs, and proteins. In this work, we use a combination of a phenomenological effective cosolute-hydrogel interaction potential and dynamic density functional theory to investigate the full nonequilibrium encapsulation kinetics of charged and dipolar cosolutes by an isolated charged hollow hydrogel immersed in a 1:1 electrolyte aqueous solution. Our analysis covers a broad spectrum of cosolute valences ( z_c) and electric dipole moments (μ_c), as well as hydrogel swelling states and hydrogel charge densities. Our calculations show that, close to the collapsed state, the polar cosolutes are predominantly precluded and the encapsulation process is strongly hindered by the excluded-volume interaction exerted by the polymer network. Different equilibrium and kinetic sorption regimes (interface versus interior) are found depending on the value and sign of z_c and the value of μ_c. For cosolutes of the same sign of charge as the gel, the superposition of steric and electrostatic repulsion leads to an "interaction-controlled" encapsulation process, in which the characteristic time to fill the empty core of the hydrogel grows exponentially with z_c. On the other hand, for cosolutes oppositely charged to the gel, we find a "diffusion-controlled" kinetic regime, where cosolutes tend to rapidly absorb into the hydrogel membrane and the encapsulation rate depends only on the cosolute diffusion time across the membrane. Finally, we find that increasing μ_c promotes the appearance of metastable and stable surface adsorption states. For large enough μ_c, the kinetics enters an "adsorption-hindered diffusion", where the enhanced surface adsorption imposes a barrier and slows down the uptake. Our study represents the first attempt to systematically describe how the swelling state of the hydrogel and other leading physical interaction parameters determine the encapsulation kinetics and the final equilibrium distribution of polar molecular cargo.

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.

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.