No Arabic abstract
We study the partitioning of cosolute particles in a thin film of a semi-flexible polymer network by a combination of coarse-grained (implicit-solvent) stochastic dynamics simulations and mean-field theory. We focus on a wide range of solvent qualities and cosolute-network interactions for selected polymer flexibilities. Our investigated ensemble (isothermal-isobaric) allows the network to undergo a volume transition from extended to collapsed state while the cosolutes can distribute in bulk and network, correspondingly. We find a rich topology of equilibrium states of the network and transitions between them, qualitatively depending on solvent quality, polymer flexibility, and cosolute-network interactions. In particular, we find a novel `cosolute-induced collapsed state, where strongly attractive cosolutes bridge network monomers albeit the latter interact mutually repulsive. Finally, the cosolutes global partitioning `landscape, computed as a function of solvent quality and cosolute-network interactions, exhibits very different topologies depending on polymer flexibility. The simulation results are supported by theoretical predictions obtained with a two-component mean-field approximation for the Helmholtz free energy that considers the chain elasticity and the particle interactions in terms of a virial expansion. Our findings have implications on the interpretation of transport processes and permeability in hydrogel films, as realized in filtration or macromolecular carrier systems.
Deformations of amorphous polymer networks prepared with significant concentrations of liquid crystalline mesogens have been recently reported to undergo mechanotropic phase transitions. Here, we report that these mechanotropic phase transitions are accompanied by an elastocaloric response ($Delta T = 2.9 text{ K}$). Applied uniaxial strain to the elastomeric polymer network transitions the organization of the material from a disordered, amorphous state (order parameter $Q=0$) to the nematic phase ($Q=0.47$). Both the magnitude of the elastocaloric temperature change and mechanically induced order parameter are dependent on the concentration of liquid crystal mesogens in the material. While the observed temperature changes in these materials are smaller than those observed in shape memory alloys, the responsivity, defined as the temperature change divided by the input stress, is larger by an order of magnitude.
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 monoatomic (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 $Ksimeq>$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 surrounding, as is known for aqueous liquid-vapor interfaces. This potential, which we find concealed in cases of symmetric monoatomic 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, $Ksimeq>$10-30 (depending on the cation type) or even $10^3$ 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.
The effect of excluded volume interactions on the structure of a polymer in shear flow is investigated by Brownian Dynamics simulations for chains with size $30leq Nleq 300$. The main results concern the structure factor $S({bf q})$ of chains of N=300 Kuhn segments, observed at a reduced shear rate $beta=dot{gamma}tau=3.2$, where $dot{gamma}$ is the bare shear rate and $tau$ is the longest relaxation time of the chain. At low q, where anisotropic global deformation is probed, the chain form factor is shown to match the form factor of the continuous Rouse model under shear at the same reduced shear rate, computed here for the first time in a wide range of wave vectors. At high q, the chain structure factor evolves towards the isotropic equilibrium power law $q^{-1/ u}$ typical of self-avoiding walk statistics. The matching between excluded volume and ideal chains at small q, and the excluded volume power law behavior at large q are observed for ${bf q}$ orthogonal to the main elongation axis but not yet for ${bf q}$ along the elongation direction itself, as a result of interferences with finite extensibility effects. Our simulations support the existence of anisotropic shear blobs for polymers in good solvent under shear flow for $beta>1$ provided chains are sufficiently long.
We extend the Rouse model of polymer dynamics to situations of non-stationary chain growth. For a dragged polymer chain of length $N(t) = t^alpha$, we find two transitions in conformational dynamics. At $alpha= 1/2$, the propagation of tension and the average shape of the chain change qualitatively, while at $alpha = 1 $ the average center-of-mass motion stops. These transitions are due to a simple physical mechanism: a race duel between tension propagation and polymer growth. Therefore they should also appear for growing semi-flexible or stiff polymers. The generalized Rouse model inherits much of the versatility of the original Rouse model: it can be efficiently simulated and it is amenable to analytical treatment.
The topological effects on the thermal properties of several knot configurations are investigated using Monte Carlo simulations. In order to check if the topology of the knots is preserved during the thermal fluctuations we propose a method that allows very fast calculations and can be easily applied to arbitrarily complex knots. As an application, the specific energy and heat capacity of the trefoil, the figure-eight and the $8_1$ knots are calculated at different temperatures and for different lengths. Short-range repulsive interactions between the monomers are assumed. The knots configurations are generated on a three-dimensional cubic lattice and sampled by means of the Wang-Landau algorithm and of the pivot method. The obtained results show that the topological effects play a key role for short-length polymers. Three temperature regimes of the growth rate of the internal energy of the system are distinguished.