Forces and torques on rigid inclusions in an elastic environment: Resulting matrix-mediated interactions, displacements, and rotations

Elasticity-controlled jamming criticality in soft composite solids

Soft composite solids are made of inclusions dispersed within soft matrices. They are ubiquitous in nature and form the basis of many biological tissues. In the field of materials science, synthetic soft composites are promising candidates for building various engineering devices due to their highly programmable features. However, when the volume fraction of the inclusions increases, predicting the mechanical properties of these materials poses a significant challenge for the classical theories of composite mechanics. The difficulty arises from the inherently disordered, multi-scale interactions between the inclusions and the matrix. To address this challenge, we systematically investigated the mechanics of densely filled soft elastomers containing stiff microspheres. We experimentally demonstrate how the strain-stiffening response of the soft composites is governed by the critical scalings in the vicinity of a shear-jamming transition of the included particles. The proposed criticality framework quantitatively connects the overall mechanics of a soft composite with the elasticity of the matrix and the particles, and captures the diverse mechanical responses observed across a wide range of material parameters. The findings uncover a novel design paradigm of composite mechanics that relies on engineering the jamming properties of the embedded inclusions.

Density functional approach to elastic properties of three-dimensional dipole-spring models for magnetic gels

Magnetic gels are composite materials consisting of a polymer matrix and embedded magnetic particles. Those are mechanically coupled to each other, giving rise to the magnetostrictive effects as well as to a controllable overall elasticity responsive to external magnetic fields. Due to their inherent composite and thereby multiscale nature, a theoretical framework bridging different levels of description is indispensable for understanding the magnetomechanical properties of magnetic gels. In this study, we extend a recently developed density functional approach from two spatial dimensions to more realistic three-dimensional systems. Along these lines, we connect a mesoscopic characterization resolving the discrete structure of the magnetic particles to macroscopic continuum parameters of magnetic gels. In particular, we incorporate the long-range nature of the magnetic dipole-dipole interaction and consider the approximate incompressibility of the embedding media and relative rotations with respect to an external magnetic field breaking rotational symmetry. We then probe the shape of the model system in its reference state, confirming the dependence of magnetostrictive effects on the configuration of the magnetic particles and on the shape of the considered sample. Moreover, calculating the elastic and rotational coefficients on the basis of our mesoscopic approach, we examine how the macroscopic types of behavior are related to the mesoscopic properties. Implications for real systems of random particle configurations are also discussed.

Effect of boundaries on displacements and motion in two-dimensional fluid or elastic films and membranes

Thin fluid or elastic films and membranes are found in nature and technology, for instance, as confinements of living cells or in loudspeakers. When applying a net force, the resulting flows in an unbounded two-dimensional incompressible low-Reynolds-number fluid or displacements in a two-dimensional linearly elastic solid seem to diverge logarithmically with the distance from the force center, which has led to some debate. Recently, we have demonstrated that such divergences cancel when the total (net) force vanishes. Here, we illustrate that if a net force is present, the boundaries play a prominent role. Already a single no-slip boundary regulates the flow and displacement fields and leads to their decay to leading order inversely in distance from a force center and the boundary. In other words, it is the boundary that stabilizes the system in this situation, unlike the three-dimensional case, where an unbounded medium by itself is able to absorb a net force. We quantify the mobility and displaceability of an inclusion as a function of the distance from the boundary, as well as interactions between different inclusions. In the case of free-slip boundary conditions, a kinked boundary is necessary to achieve stabilization.

A simple model for elastic wave propagation in hard sphere-filled random composites

A simple model is proposed to study wave propagation in hard sphere-reinforced elastic random composites. Compared to existing related models, the proposed model is featured by a modified form of classical elastodynamic equations in which the inertia term is substituted by the acceleration field of the mass centre of a representative unit cell, supplied with a derived simple differential relation between the displacement field of the composite and the displacement field of the mass centre of a representative unit cell. The present model enjoys conceptual and mathematical simplicity although it is restricted to hard sphere-filled elastic composites in which the elastic moduli of embedded spheres are much (at least 4-5 times) stiffer than those of a softer matrix. Explicit formulas are derived for the attenuation coefficient and the effective phase velocity of plane longitudinal P-waves and transverse S-waves. The efficiency and reasonable accuracy of the present model are demonstrated by reasonably good agreement between the predicted results and some established known data. The proposed model could offer a potential general method to study various three-dimensional dynamic problems of hard sphere-filled elastic random composites.

Magnetization of magnetoactive elastomers under the assumption of breakable adhesion at the particle/matrix interface

In this work we study the magnetization of magnetoactive elastomers (MAE) in which the interface between the matrix and magnetic particles is unstable and allows for slipping of the particles against the wall of their elastomer cavities. The estimate of the maximal angle at which each particle can decline its axis from the initial position is made based on cyclic measurement of several consecutive hysteresis loops at different maximal magnetic fields. A model of magnetization of magnetically hard multigrain particles in an elastic environment with allowance for their possible slipping is proposed. Results of modelling is in fair agreement with the experimental data obtained on MAEs whose polymeric matrix is made of polydimethylsiloxane and the magnetic filler is NdFeB spherical particles.

Mediated interactions between rigid inclusions in two-dimensional elastic or fluid films

Interactions between rigid inclusions in continuous three-dimensional linearly elastic solids and low-Reynolds-number viscous fluids have largely been quantified in the past. Prime example systems are given by functionalized elastic composite materials or fluid colloidal suspensions. Here, we address the significantly less frequently studied situation of rigid inclusions in two-dimensional elastic or low-Reynolds-number fluid films. We concentrate on the situation in which disklike inclusions remain well separated from each other and do not get into contact. Specifically, we demonstrate and explain that the logarithmic divergence of the associated Green's function is removed in the absence of net external forces on the inclusions, in line with physical intuition. For instance, this situation applies when only pairwise mutual interactions between the inclusions prevail. Our results will support, for example, investigations on membranes functionalized by appropriate inclusions, both of technical or biological origin, or the dynamics of active microswimmers in appropriately prepared thin films.

Collision and separation of nickel particles embedded in a polydimethylsiloxan matrix under a rotating magnetic field: A strong magneto active function

In order to function as soft actuators, depending on their field of use, magnetorheological elastomers (MREs) must fulfill certain criteria. To name just a few, these can include rapid response to external magnetic fields, mechanical durability, mechanical strength, and/or large deformation. Of particular interest are MREs which produce macroscopic deformation for small external magnetic field variations. This work demonstrates how this can be achieved by just a small change in magnetic field orientation. To achieve this, (super)paramagnetic nickel particles of size ≈ 160 μm were embedded in a non-magnetic polydimethylsiloxan (PDMS) (661–1301 Pa) and their displacement in a stepwise rotated magnetic field (170 mT) recorded using a video microscope. Changes in particle aggregation resulting from very small variations in magnetic field orientation led to the observation of a new strongly magneto-active effect. This configuration is characterized by an interparticle distance in relation to the angle difference between magnetic field and particle axis. This causes a strong matrix deformation which in turn demonstrates hysteresis on relaxation. It is shown that the occurrence strongly depends on the particle size, particle distance, and stiffness of the matrix. Choosing the correct parameter combination, the state can be suppressed and the particle-matrix system demonstrates no displacement or hysteresis. In addition, evidences of non-negligible higher order magnetization effects are experimentally ascertained which is qualitatively in agreement with similar, already theoretically described, particle systems. Even at larger particle geometries, the new strongly magneto-active configuration is preserved and could create macroscopic deformation changes.

Magneto-Mechanical Coupling in Magneto-Active Elastomers

In the present work, the magneto-mechanical coupling in magneto-active elastomers is investigated from two different modeling perspectives: a micro-continuum and a particle–interaction approach. Since both strategies differ significantly in their basic assumptions and the resolution of the problem under investigation, they are introduced in a concise manner and their capabilities are illustrated by means of representative examples. To motivate the application of these strategies within a hybrid multiscale framework for magneto-active elastomers, their interchangeability is then examined in a systematic comparison of the model predictions with regard to the magneto-deformation of chain-like helical structures in an elastomer surrounding. The presented results show a remarkable agreement of both modeling approaches and help to provide an improved understanding of the interactions in magneto-active elastomers with chain-like microstructures.

Modeling and theoretical description of magnetic hybrid materials—bridging from meso- to macro-scales

Magnetic gels and elastomers consist of magnetic or magnetizable colloidal particles embedded in an elastic polymeric matrix. Outstanding properties of these materials comprise reversible changes in their mechanical stiffness or magnetostrictive distortions under the influence of external magnetic fields. To understand such types of overall material behavior from a theoretical point of view, it is essential to characterize the substances starting from the discrete colloidal particle level. It turns out that the macroscopic material response depends sensitively on the mesoscopic particle arrangement. We have utilized and developed several theoretical approaches to this end, allowing us both to reproduce experimental observations and to make theoretical predictions. Our hope is that both these paths help to further stimulate the interest in these fascinating materials.

Magnetic torque-driven deformation of Ni-nanorod/hydrogel nanocomposites

Nickel (Ni) nanorods were prepared by the anodized aluminum oxide (AAO) template method and dispersed in poly(acrylamide) (PAM) hydrogels. The deformation of the magnetoresponsive composites was studied with particular attention to the consequences of finite magnetic shape anisotropy as compared to rigid dipoles on the field-dependent torque. For comparison with experiments, the composite was described as an elastic continuum with a local magnetic torque density, applied by discrete particles and determined by the local orientation of their magnetic anisotropy axis with respect to the magnetic field. The mean magnetic moment of the single domain particles m and their volume density in the composite φ vol were derived from the static field-dependent optical transmission (SFOT) of linear polarized light. The mechanical coupling between the particles and their viscoelastic environment was retrieved from the rotational dynamics of the nanorods using oscillating field-dependent optical transmission (OFOT) measurements. Field- and orientation-dependent magnetization measurements were analyzed using the Stoner–Wohlfarth (SW) model and a valid parameter range was identified by introducing an effective anisotropy constant K A as a new empirical model parameter. This adapted SW-model for quantitative description of the field- and orientation dependence of the magnetic torque was validated by measuring the local rotation of nanorods in a soft elastic hydrogel. Finally, torsional and bending deformation of thin magnetically textured composite filaments were computed and compared with experiments.

Ordering of adsorbed rigid rods mediated by the Boussinesq interaction on a soft substrate

Orientational ordering driven by mechanical distortion of soft substrates plays a major role in material transformation processes such as elastocapillarity and surface anchoring. We present a theoretical model of the orientational response of anisotropic rods deposited onto a surface of a soft, elastic substrate of finite thickness. We show that anisotropic rods exhibit a continuous isotropic-nematic phase transition, driven by orientational interactions between surface deposited rods. This interaction is mediated by the deformation of the underlying elastic substrate and is quantified by the Boussinesq solution adapted to the case of slender, surface deposited rods. From the microscopic rod-rod interactions, we derive the appropriate Maier-Saupe mean-field description, which includes the Boussinesq elastic free energy contribution due to the substrate elasticity, derive the conditions for the existence of a continuous orientational ordering transition, and discuss the implication of results in the soft (bio)system context.

Classical density functional theory for a two-dimensional isotropic ferrogel model with labeled particles

In this study, we formulate a density functional theory (DFT) for systems of labeled particles, considering a two-dimensional bead-spring lattice with a magnetic dipole on every bead as a model for ferrogels. On the one hand, DFT has been widely studied to investigate fluidlike states of materials, in which constituent particles are not labeled as they can exchange their positions without energy cost. On the other hand, in ferrogels consisting of magnetic particles embedded in elastic polymer matrices, the particles are labeled by their positions as their neighbors do not change over time. We resolve such an issue of particle labeling, introducing a mapping of the elastic interaction mediated by springs onto a pairwise additive interaction (pseudosprings) between unlabeled particles. We further investigate magnetostriction and changes in the elastic constants under altered magnetic interactions employing the pseudospring potential. It is revealed that there are two different response scenarios in the mechanical properties of the dipole-spring systems: While systems at low packing fractions are hardened as the magnetic moments increase in magnitude, at high packing fractions softening due to diminishing effects from the steric force, associated with increases in the volume, is observed. The validity of the theory is also verified by Monte Carlo simulations with both real springs and pseudosprings. We expect that our DFT approach may promote our understanding of materials with particle inclusions.

The loss of mechanical stability and the critical magnetization of a ferromagnetic particle in an elastomer

The rotation of a spherical magnetically soft ferromagnetic particle with uniaxial magnetic anisotropy and located in an elastomer through a certain angle under the action of a magnetic field has been studied theoretically. It was found that if the particle loses its mechanical stability, its rotation angle becomes critically dependent on the magnetic field magnitude. The magnetically induced critical rotation of the particle in the elastomer has a magneto-elastic origin and is accompanied by a critical change in the magnetization of the particle multidomain state, as occurs in the case of second-order magnetic phase transitions. The transition of a particle in a soft-elastic elastomer from the multidomain state into a uniformly magnetized one was shown to be independent of the magnitude of the particle magnetic anisotropy field. The particle rotation was found to result in the appearance of a maximum in the field dependence of the magnetic susceptibility.

Displacement field around a rigid sphere in a compressible elastic environment, corresponding higher-order Faxén relations, as well as higher-order displaceability and rotateability matrices

An efficient route to the displacement field around a rigid spherical inclusion in an infinitely extended homogeneous elastic medium is presented in a slightly alternative way when compared to some common textbook methods. Moreover, two Faxén relations of next-higher order beyond the stresslet are calculated explicitly for compressible media. They quantify higher-order moments involving the force distribution on a rigid spherical particle in a deformed elastic medium. As a consequence, additional contributions to the distortions of the deformed elastic medium are identified that are absent to lower order. Furthermore, the displaceability and rotateability matrices for an ensemble of rigid spheres are calculated up to (including) sixth order in inverse particle separation distance. These matrices describe the interactions mediated between the rigid embedded particles by the elastic environment. In this way, additional coupling effects are identified that are absent to lower order, particularly when rotations and torques are involved. All methods and results can formally be transferred to the corresponding case of incompressible hydrodynamic low-Reynolds-number Stokes flow by considering the limit of an incompressible environment. The roles of compressibility of the embedding medium and of the here additionally derived higher-order contributions are highlighted by some selected example configurations.

Mechanical Model of Nuclei Ordering in Drosophila Embryos Reveals Dilution of Stochastic Forces

During the initial development of syncytial embryos, nuclei go through cycles of nuclear division and spatial rearrangement. The arising spatial pattern of nuclei is important for subsequent cellularization and morphing of the embryo. Although nuclei are contained within a common cytoplasm, cytoskeletal proteins are nonuniformly packaged into regions around every nucleus. In fact, cytoskeletal elements like microtubules and their associated motor proteins exert stochastic forces between nuclei, actively driving their rearrangement. Yet, it is unknown how the stochastic forces are balanced to maintain nuclear order in light of increased nuclear density upon every round of divisions. Here, we investigate the nuclear arrangements in Drosophila melanogaster over the course of several nuclear divisions starting from interphase 11. We develop a theoretical model in which we distinguish long-ranged passive forces due to the nuclei as inclusions in the elastic matrix, namely the cytoplasm, and active, stochastic forces arising from the cytoskeletal dynamics mediated by motor proteins. We perform computer simulations and quantify the observed degree of orientational and spatial order of nuclei. Solely doubling the nuclear density upon nuclear division, the model predicts a decrease in nuclear order. Comparing results to experimental recordings of tracked nuclei, we make contradictory observations, finding an increase in nuclear order upon nuclear divisions. Our analysis of model parameters resulting from this comparison suggests that overall motor protein density as well as relative active-force amplitude has to decrease by a factor of about two upon nuclear division to match experimental observations. We therefore expect a dilution of cytoskeletal motors during the rapid nuclear division to account for the increase in nuclear order during syncytial embryo development. Experimental measurements of kinesin-5 cluster lifetimes support this theoretical finding.

Memory-based mediated interactions between rigid particulate inclusions in viscoelastic environments

Many practically relevant materials combine properties of viscous fluids and elastic solids to viscoelastic behavior. Our focus is on the induced dynamic behavior of damped finite-sized particulate inclusions in such substances. We explicitly describe history-dependent interactions that emerge between the embedded particles. These interactions are mediated by the viscoelastic surroundings. They result from the flows and distortions of the viscoelastic medium when induced by the rigid inclusions. Both viscoelastic environments of terminal fluidlike flow and of completely reversible damped elastic behavior are covered. For illustration and to highlight the role of the formalism in potential applications, we briefly address the relevant examples of dragging a rigid sphere through a viscoelastic environment together with subsequent relaxation dynamics, the switching dynamics of magnetic fillers in elastic gel matrices, and the swimming behavior of active microswimmers in viscoelastic solutions. The approach provides a basis for more quantitative and extended investigations of these and related systems in the future.

Reversible magnetomechanical collapse: virtual touching and detachment of rigid inclusions in a soft elastic matrix

Soft elastic composite materials containing particulate rigid inclusions in a soft elastic matrix are candidates for developing soft actuators or tunable damping devices. The possibility to reversibly drive the rigid inclusions within such a composite together to a close-to-touching state by an external stimulus would offer important benefits. Then, a significant tuning of the mechanical properties could be achieved due to the resulting mechanical hardening. For a long time, it has been argued whether a virtual touching of the embedded magnetic particles with subsequent detachment can actually be observed in real materials, and if so, whether the process is reversible. Here, we present experimental results that demonstrate this phenomenon in reality. Our system consists of two paramagnetic nickel particles embedded at finite initial distance in a soft elastic polymeric gel matrix. Magnetization in an external magnetic field tunes the magnetic attraction between the particles and drives the process. We quantify our experimental results by different theoretical tools, i.e., explicit analytical calculations in the framework of linear elasticity theory, a projection onto simplified dipole-spring models, as well as detailed finite-element simulations. From these different approaches, we conclude that in our case the cycle of virtual touching and detachment shows hysteretic behavior due to the mutual magnetization between the paramagnetic particles. Our results are important for the design and construction of reversibly tunable mechanical damping devices. Moreover, our projection on dipole-spring models allows the formal connection of our description to various related systems, e.g., magnetosome filaments in magnetotactic bacteria.

Tunable dynamic moduli of magnetic elastomers: from characterization by x-ray micro-computed tomography to mesoscopic modeling

Ferrogels and magnetorheological elastomers are composite materials obtained by embedding magnetic particles of mesoscopic size in a crosslinked polymeric matrix. They combine the reversible elastic deformability of polymeric materials with the high responsivity of ferrofluids to external magnetic fields. These materials stand out, for example, for significant magnetostriction as well as a pronounced increase of the elastic moduli in the presence of external magnetic fields. By means of x-ray micro-computed tomography, the position and size of each magnetic particle can be measured with a high degree of accuracy. We here use data extracted from real magnetoelastic samples as input for coarse-grained dipole-spring modeling and calculations to investigate internal restructuring, stiffening, and changes in the normal modes spectrum. More precisely, we assign to each particle a dipole moment proportional to its volume and set a randomized network of springs between them that mimics the behavior of the polymeric elastic matrix. Extending our previously developed methods, we compute the resulting structural changes in the systems as well as the frequency-dependent elastic moduli when magnetic interactions are turned on. Particularly, with increasing magnetization, we observe the formation of chain-like aggregates. Interestingly, the static elastic moduli can first show a slight decrease with growing amplitude of the magnetic interactions, before a pronounced increase appears upon the chain formation. The change of the dynamic moduli with increasing magnetization depends on the frequency and can even feature nonmonotonic behavior. Overall, we demonstrate how theory and experiments can complement each other to learn more about the dynamic behavior of this interesting class of materials.

Dynamics in a one-dimensional ferrogel model: relaxation, pairing, shock-wave propagation

Based on theory and simulations, we elucidate the relaxation dynamics of a one-dimensional ferrogel model and provide classification scenarios.

A density functional approach to ferrogels

Ferrogels consist of magnetic colloidal particles embedded in an elastic polymer matrix. As a consequence, their structural and rheological properties are governed by a competition between magnetic particle-particle interactions and mechanical matrix elasticity. Typically, the particles are permanently fixed within the matrix, which makes them distinguishable by their positions. Over time, particle neighbors do not change due to the fixation by the matrix. Here we present a classical density functional approach for such ferrogels. We map the elastic matrix-induced interactions between neighboring colloidal particles distinguishable by their positions onto effective pairwise interactions between indistinguishable particles similar to a 'pairwise pseudopotential'. Using Monte-Carlo computer simulations, we demonstrate for one-dimensional dipole-spring models of ferrogels that this mapping is justified. We then use the pseudopotential as an input into classical density functional theory of inhomogeneous fluids and predict the bulk elastic modulus of the ferrogel under various conditions. In addition, we propose the use of an 'external pseudopotential' when one switches from the viewpoint of a one-dimensional dipole-spring object to a one-dimensional chain embedded in an infinitely extended bulk matrix. Our mapping approach paves the way to describe various inhomogeneous situations of ferrogels using classical density functional concepts of inhomogeneous fluids.