Determination of viscoelastic properties by analysis of probe-particle motion in molecular simulations

Rheological Considerations of Pharmaceutical Formulations: Focus on Viscoelasticity

Controlling rheological properties offers the opportunity to gain insight into the physical characteristics, structure, stability and drug release rate of formulations. To better understand the physical properties of hydrogels, not only rotational but also oscillatory experiments should be performed. Viscoelastic properties, including elastic and viscous properties, are measured using oscillatory rheology. The gel strength and elasticity of hydrogels are of great importance for pharmaceutical development as the application of viscoelastic preparations has considerably expanded in recent decades. Viscosupplementation, ophthalmic surgery and tissue engineering are just a few examples from the wide range of possible applications of viscoelastic hydrogels. Hyaluronic acid, alginate, gellan gum, pectin and chitosan are remarkable representatives of gelling agents that attract great attention applied in biomedical fields. This review provides a brief summary of rheological properties, highlighting the viscoelasticity of hydrogels with great potential in biomedicine.

Nanorheology of Entangled Polymer Melts

We use molecular simulations to probe the local viscoelasticity of an entangled polymer melt by tracking the motion of embedded nonsticky nanoparticles (NPs). As in conventional microrheology, the generalized Stokes-Einstein relation is employed to extract an effective stress relaxation function G_{GSE}(t) from the mean square displacement of NPs. G_{GSE}(t) for different NP diameters d are compared with the stress relaxation function G(t) of a pure polymer melt. The deviation of G_{GSE}(t) from G(t) reflects the incomplete coupling between NPs and the dynamic modes of the melt. For linear polymers, a plateau in G_{GSE}(t) emerges as d exceeds the entanglement mesh size a and approaches the entanglement plateau in G(t) for a pure melt with increasing d. For ring polymers, as d increases towards the spanning size R of ring polymers, G_{GSE}(t) approaches G(t) of the ring melt with no entanglement plateau.

Determination of linear viscoelastic properties of an entangled polymer melt by probe rheology simulations

Particle rheology is used to extract the linear viscoelastic properties of an entangled polymer melt from molecular dynamics simulations. The motion of a stiff, approximately spherical particle is tracked in both passive and active modes. We demonstrate that the dynamic modulus of the melt can be extracted under certain limitations using this technique. As shown before for unentangled chains [Karim et al., Phys. Rev. E 86, 051501 (2012)PLEEE81539-375510.1103/PhysRevE.86.051501], the frequency range of applicability is substantially expanded when both particle and medium inertia are properly accounted for by using our inertial version of the generalized Stokes-Einstein relation (IGSER). The system used here introduces an entanglement length d_{T}, in addition to those length scales already relevant: monomer bead size d, probe size R, polymer radius of gyration R_{g}, simulation box size L, shear wave penetration length Δ, and wave period Λ. Previously, we demonstrated a number of restrictions necessary to obtain the relevant fluid properties: continuum approximation breaks down when d≳Λ; medium inertia is important and IGSER is required when R≳Λ; and the probe should not experience hydrodynamic interaction with its periodic images, L≳Δ. These restrictions are also observed here. A simple scaling argument for entangled polymers shows that the simulation box size must scale with polymer molecular weight as M_{w}^{3}. Continuum analysis requires the existence of an added mass to the probe particle from the entrained medium but was not observed in the earlier work for unentangled chains. We confirm here that this added mass is necessary only when the thickness L_{S} of the shell around the particle that contains the added mass, L_{S}>d. We also demonstrate that the IGSER can be used to predict particle displacement over a given timescale from knowledge of medium viscoelasticity; such ability will be of interest for designing nanoparticle-based drug delivery.

Active one-particle microrheology of an unentangled polymer melt studied by molecular dynamics simulation

We present molecular dynamics simulations for active one-particle microrheology of an unentangled polymer melt. The tracer particle is forced to oscillate by an oscillating harmonic potential, which models an experiment using optical tweezers. The amplitude and phase shift of this oscillation are related to the complex shear modulus of the polymer melt. In the linear response regime at low frequencies, the active microrheology gives the same result as the passive microrheology, where the thermal motion of a tracer particle is related to the complex modulus. We expand the analysis to include full hydrodynamic effects instead of stationary Stokes friction only, and show that different approaches suggested in the literature lead to completely different results, and that none of them improves on the description using the stationary Stokes friction.

Trajectories of probe spheres in generalized linear viscoelastic complex fluids

We have developed a fast simulation that generates a random walk of an isolated probe sphere in a generalized linear viscoelastic complex fluid over a highly extended dynamic range. We introduce a coupled harmonically bound Brownian particle (c-HBBP) model, in which the relaxation modes of the viscoelastic medium are treated as harmonic wells. These wells are coupled to the probe sphere and perform Brownian motion in bound harmonic potentials corresponding to the next-longer relaxation mode, according to the relaxation spectrum of the viscoelastic material. We implement this c-HBBP model by generating variable temporal step sizes that have a uniform distribution in logarithmic time. We create and analyze trajectories for several different viscoelastic complex fluids: a polymer system at its gel point, a dense emulsion system, a blend of two monodisperse polystyrene polymers for which the relaxation spectrum has been measured, and a model anisotropic soft system that shows dense emulsion-like and gel-point behaviors along two orthogonal directions. Except for unusual viscoelastic materials, such as the polymer system at its gel point, the generated trajectories are neither self-similar nor self-affine. The resulting mean square displacements predicted by the c-HBBP model are consistent with the single-particle generalized Stokes-Einstein relation of linear passive microrheology.

Passive one-particle microrheology of an unentangled polymer melt studied by molecular dynamics simulation

We present a molecular dynamics simulation study of the possibility of performing a microrheological analysis of a polymer melt by following the Brownian motion of a dispersed nanoparticle. We study the influence of the size of the nanoparticle, taken to be comparable to the radius of gyration of the chains, and of the strength of the interaction between the nanoparticle and the repeat units of the polymer chains. The influence of the presence of the nanoparticle on the melt mechanical behavior is analyzed, and the importance of effects of different levels of hydrodynamic analysis on the frequency-dependent dynamic shear modulus derived from the particle motion is worked out.