Search

Hydrogels are particularly versatile materials that are widely found in both Nature and industry. One key reason for this versatility is their high water content, which lets them dramatically change their volume and many of their mechanical properties -- often by orders of magnitude -- as they swell and dry out. Currently, we lack techniques that can precisely characterize how these properties change with water content. To overcome this challenge, here we develop Gel-Freezing Osmometry (GelFrO): an extension of freezing-point osmometry. We show how GelFrO can measure a hydrogel's mechanical response to compression and osmotic pressure, while only using small, $O(100\mu$L$)$ samples. Because the technique allows measurement of properties over an unusually wide range of water contents, it allows us to accurately test theoretical predictions. We find simple, power-law behavior for both mechanical response to compression, and osmotic pressure, while these are not well-captured by classical Flory-Huggins theory. We interpret this power-law behavior as a hallmark of a microscopic fractal structure of the gel's polymer network, and propose a simple way to connect the gel's fractal dimension to its mechanical and osmotic properties. This connection is supported by observations of hydrogel microstructures using small-angle x-ray scattering. Finally, our results motivate us to propose an updated constitutive model describing hydrogel swelling, and mechanical response.

We propose the Self Returning Excluded Volume (SR-EV) model for the structure of chromatin based on stochastic rules and physical interactions. The SR-EV rules of return generate conformationally-defined domains observed by single cell imaging techniques. From nucleosome to chromosome scales, the model captures the overall chromatin organization as a corrugated system, with dense and dilute regions alternating in a manner that resembles the mixing of two disordered bi-continuous phases. This particular organizational topology is a consequence of the multiplicity of interactions and processes occurring in the nuclei, and mimicked by the proposed return rules. Single configuration properties and ensemble averages show a robust agreement between theoretical and experimental results including chromatin volume concentration, contact probability, packing domain identification and size characterization, and packing scaling behavior. Model and experimental results suggest that there is an inherent chromatin organization regardless of the cell character and resistant to an external forcing such as Rad21 degradation.

Fracture phenomena in soft materials span multiple length- and timescales. This poses a major challenge in computational modeling and predictive materials design. To pass quantitatively from molecular- to continuum scales, a precise representation of the material response at the molecular level is vital. Here, we derive the nonlinear elastic response and fracture characteristics of individual siloxane molecules using molecular dynamics (MD) studies. For short chains, we find deviations from classical scalings for both the effective stiffness and mean chain rupture times. A simple model of a non-uniform chain of Kuhn segments captures the observed effect and agrees well with MD data. We find that the dominating fracture mechanism depends on the applied force scale in a non-monotonic fashion. This analysis suggests that common polydimethylsiloxane (PDMS) networks fail at crosslinking points. Our results can be readily lumped into coarse-grained models. Although focusing on PDMS as a model system, our study presents a general procedure to pass beyond the window of accessible rupture times in MD studies employing mean first passage time theory, which can be exploited for arbitrary molecular systems.

Rod-shaped particles embedded in certain matrices have been reported to exhibit an increase in their center of mass diffusivity upon increasing the matrix density. This increase has been considered to be caused by a kinetic constraint in analogy with tube models. We investigate a mobile rod-like particle in a sea of immobile point obstacles using a kinetic Monte Carlo scheme equipped with a Markovian process, that generates gas-like collision statistics, so that such kinetic constraints do essentially not exist. Even in such a system, provided the particle's aspect ratio exceeds a threshold value of about 24, the unusual increase in the rod diffusivity emerges. This result implies that the kinetic constraint is not a necessary condition for the increase in the diffusivity., Comment: 9 pages, 10 figures

By combining molecular dynamics simulations and topological analyses with scaling arguments, we obtain analytic expressions that quantitatively predict the entanglement length $N_e$, the plateau modulus $G$, and the tube diameter $a$ in melts that span the entire range of chain stiffnesses for which systems remain isotropic. Our expressions resolve conflicts between previous scaling predictions for the loosely entangled [Lin-Noolandi: $G\ell_K^3/k_\textrm{B}T \sim (\ell_K/p)^3$], semiflexible [Edwards/de Gennes: $G\ell_K^3/k_\textrm{B}T \sim (\ell_K/p)^2$], and tightly-entangled [Morse: $G\ell_K^3/k_\textrm{B}T \sim (\ell_K/p)^{1+\epsilon}$] regimes, where $\ell_K$ and $p$ are respectively the Kuhn and packing lengths. We also find that maximal entanglement (minimal $N_e$) coincides with the onset of local nematic order., Comment: Accepted for publication in Physical Review Letters

We study the time correlation functions of coupled linear Langevin dynamics without and with inertia effects, both analytically and numerically. The model equation represents the physical behavior of a harmonic oscillator in two or three dimensions in the presence of friction, noise, and an external field with both rotational and deformational components. This simple model plays pivotal roles in understanding more complicated processes. The presented analytical solution serves as a test of numerical integration schemes, its derivation is presented in a fashion that allows to be repeated directly in a classroom. While the results in the absence of fields (equilibrium) or confinement (free particle) are omnipresent in the literature, we write down, apparently for the first time, the full nonequilibrium results that may correspond, e.g., to a Hookean dumbbell embedded in a macroscopically homogeneous shear or mixed flow field. We demonstrate how the inertia results reduce to their noninertia counterparts in the nontrivial limit of vanishing mass. While the results are derived using basic integrations over Dirac delta distributions, we mention its relationship with alternative approaches involving (i) Fourier transforms, that seems advantageous only if the measured quantities also reside in Fourier space, and (ii) a Fokker--Planck equation and the moments of the probability distribution. The results, verified by numerical experiments, provide additional means of measuring the performance of numerical methods for such systems. It should be emphasized that this manuscript provides specific details regarding the derivations of the time correlation functions as well as the implementations of various numerical methods, so that it can serve as a standalone piece as part of education in the framework of stochastic differential equations and calculus.

We investigate the conformation, position, and dynamics of core-shell nanoparticles (CSNPs) composed of a silica core encapsulated in a cross-linked poly-N-isopropylacrylamide shell at a water-oil interface for a systematic range of core sizes and shell thicknesses. We first present a free-energy model that we use to predict the CSNP wetting behavior at the interface as a function of its geometrical and compositional properties in the bulk phases, which gives good agreement with our experimental data. Remarkably, upon knowledge of the polymer shell deformability, the equilibrium particle position relative to the interface plane, an often elusive experimental quantity, can be extracted by measuring its radial dimensions after adsorption. For all the systems studied here, the interfacial dimensions are always larger than in bulk and the particle core resides in a configuration wherein it just touches the interface or is fully immersed in water. Moreover, the stretched shell induces a larger viscous drag at the interface, which appears to depend solely on the interfacial dimensions, irrespective of the portion of the CSNP surface exposed to the two fluids. Our findings indicate that tailoring the architecture of CSNPs can be used to control their properties at the interface, as of interest for applications including emulsion stabilization and nanopatterning., Comment: 38 pages, 12 figures

We discuss the design of state-of-the-art numerical methods for molecular dynamics, focusing on the demands of soft matter simulation, where the purposes include sampling and dynamics calculations both in and out of equilibrium. We discuss the characteristics of different algorithms, including their essential conservation properties, the convergence of averages, and the accuracy of numerical discretizations. Formulations of the equations of motion which are suited to both equilibrium and nonequilibrium simulation include Langevin dynamics, dissipative particle dynamics (DPD), and the more recently proposed "pairwise adaptive Langevin" (PAdL) method, which, like DPD but unlike Langevin dynamics, conserves momentum and better matches the relaxation rate of orientational degrees of freedom. PAdL is easy to code and suitable for a variety of problems in nonequilibrium soft matter modeling, our simulations of polymer melts indicate that this method can also provide dramatic improvements in computational efficiency. Moreover we show that PAdL gives excellent control of the relaxation rate to equilibrium. In the nonequilibrium setting, we further demonstrate that while PAdL allows the recovery of accurate shear viscosities at higher shear rates than are possible using the DPD method at identical timestep, it also outperforms Langevin dynamics in terms of stability and accuracy at higher shear rates.