Statistical field theory of polarizable polymer chains with nonlocal dipolar interactions

The electromechanical response of polymeric soft matter to applied electric fields is of fundamental scientific interest as well as relevant to technologies for sensing and actuation. Several existing theoretical and numerical approaches for polarizable polymers subject to a combined applied electric field and stretch are based on discrete monomer models. In these models, accounting for the interactions between the induced dipoles on monomers is challenging due to the nonlocality of these interactions. On the other hand, the framework of statistical field theory provides a continuous description of polymer chains that potentially enables a tractable way to account for these interactions. However, prior formulations using this framework have been restricted to the case of weak anisotropy of the monomer polarizability. This paper formulates a general approach based in the framework of statistical field theory to account for the nonlocal nature of the dipolar interactions without any restrictions on the anisotropy or nonlinearity of the polarizability of the monomer. The approach is based on three key elements: (1) the statistical field theory framework, in which the discrete monomers are regularized to a continuous dipole distribution, (2) a replacement of the nonlocal dipole-dipole interactions by the local electrostatics partial differential equation with the continuous dipole distribution as the forcing, and (3) the use of a completely general relation between the polarization and the local electric field. Rather than treat the dipole-dipole interactions directly, the continuous description in the field theory enables the computationally tractable nonlocal-to-local transformation. Further, it enables the use of a realistic statistical-mechanical ensemble wherein the average far-field applied electric field is prescribed, rather than prescribing the applied field at every point in the polymer domain. The model is applied, using the finite element method, to study the electromechanical response of a polymer chain in the ensemble with fixed far-field applied electric field and fixed chain stretch. The nonlocal dipolar interactions are found to increase, over the case where dipole-dipole interactions are neglected, the magnitudes of the polarization and electric field by orders of magnitude as well as significantly change their spatial distributions. Next, the effect of the relative orientation between the applied field and the chain on the local electric field and polarization is studied. The model predicts that the elastic response of the polymer chain is linear, consistent with the Gaussian approximation, and largely unchanged by the orientation of the applied electric field, though the polarization and local electric field distributions are significantly impacted.

A dimensionally-reduced nonlinear elasticity model for liquid crystal elastomer strips with transverse curvature

Liquid crystalline elastomers (LCEs) are active materials that are of interest due to their programmable response to various external stimuli such as light and heat. When exposed to these stimuli, the anisotropy in the response of the material is governed by the nematic director, which is a continuum parameter that is defined as the average local orientation of the mesogens in the liquid crystal phase. This nematic director can be programmed to be heterogeneous in space, creating a vast design space that is useful for applications ranging from artificial ligaments to deployable structures to self-assembling mechanisms. Even when specialized to long and thin strips of LCEs - the focus of this work - the vast design space has required the use of numerical simulations to aid in experimental discovery. To mitigate the computational expense of full 3-d numerical simulations, several dimensionally-reduced rod and ribbon models have been developed for LCE strips, but these have not accounted for the possibility of initial transverse curvature, like carpenter's tape spring. Motivated by recent experiments showing that transversely-curved LCE strips display a rich variety of configurations, this work derives a dimensionally-reduced 1-d model for pre-curved LCE strips. The 1-d model is validated against full 3-d finite element calculations, and it is also shown to capture experimental observations, including tape-spring-like localizations, in activated LCE strips.

Statistical field theory for nonlinear elasticity of polymer networks with excluded volume interactions

Polymer networks formed by cross linking flexible polymer chains are ubiquitous in many natural and synthetic soft-matter systems. Current micromechanics models generally do not account for excluded volume interactions except, for instance, through imposing a phenomenological incompressibility constraint at the continuum scale. This work aims to examine the role of excluded volume interactions on the mechanical response. The approach is based on the framework of the self-consistent statistical field theory of polymers, which provides an efficient mesoscale approach that enables the accounting of excluded volume effects without the expense of large-scale molecular modeling. A mesoscale representative volume element is populated with multiple interacting chains, and the macroscale nonlinear elastic deformation is imposed by mapping the end-to-end vectors of the chains by this deformation. In the absence of excluded volume interactions, it recovers the closed-form results of the classical theory of rubber elasticity. With excluded volume interactions, the model is solved numerically in three dimensions using a finite element method to obtain the energy, stresses, and linearized moduli under imposed macroscale deformation. Highlights of the numerical study include: (i) the linearized Poisson's ratio is very close to the incompressible limit without a phenomenological imposition of incompressibility; (ii) despite the harmonic Gaussian chain as a starting point, there is an emergent strain-softening and strain-stiffening response that is characteristic of real polymer networks, driven by the interplay between the entropy and the excluded volume interactions; and (iii) the emergence of a deformation-sensitive localization instability at large excluded volumes.

Hierarchically structured bioinspired nanocomposites

Next-generation structural materials are expected to be lightweight, high-strength and tough composites with embedded functionalities to sense, adapt, self-repair, morph and restore. This Review highlights recent developments and concepts in bioinspired nanocomposites, emphasizing tailoring of the architecture, interphases and confinement to achieve dynamic and synergetic responses. We highlight cornerstone examples from natural materials with unique mechanical property combinations based on relatively simple building blocks produced in aqueous environments under ambient conditions. A particular focus is on structural hierarchies across multiple length scales to achieve multifunctionality and robustness. We further discuss recent advances, trends and emerging opportunities for combining biological and synthetic components, state-of-the-art characterization and modelling approaches to assess the physical principles underlying nature-inspired design and mechanical responses at multiple length scales. These multidisciplinary approaches promote the synergetic enhancement of individual materials properties and an improved predictive and prescriptive design of the next era of structural materials at multilength scales for a wide range of applications.

Nonlinear statistical mechanics drives intrinsic electrostriction and volumetric torque in polymer networks

Statistical mechanics is an important tool for understanding polymer electroelasticity because the elasticity of polymers is primarily due to entropy. However, a common approach for the statistical mechanics of polymer chains, the Gaussian chain approximation, misses key physics. By considering the nonlinearities of the problem, we show a strong coupling between the deformation of a polymer chain and its dielectric response, that is, its net dipole. When chains with this coupling are cross linked in an elastomer network and an electric field is applied, the field breaks the symmetry of the elastomer's elastic properties and, combined with electrostatic torque and incompressibility, leads to intrinsic electrostriction. Conversely, deformation can break the symmetry of the dielectric response, leading to volumetric torque and asymmetric actuation. Both phenomena have important implications for designing high-efficiency soft actuators and soft electroactive materials, and the presence of mechanisms for volumetric torque, in particular, can be used to develop higher degree of freedom actuators and to achieve bioinspired locomotion.

Torque-dense photomechanical actuation

Contactless actuation powered using light is shown to generate torque densities approaching 10 N m kg-1 at angular velocities ∼102 rad s-1: metrics that compare favorably against tethered electromechanical systems. This is possible even though the extinction of actinic light limits the characteristic thickness of photoresponse in polymers to tens of μm. Confinement of molecularly patterned developable shells fabricated from azobenzene-functionalized liquid crystalline polymers encodes torque-dense photoactuation. Photostrain gradients from unstructured irradiation segment this geometry into two oppositely curved regions connected by a curved crease. A monolithic curved shell spontaneously bifurcates into a jointed, arm-like mechanism that generates flexure over sweep angles exceeding a radian. Strain focusing at the crease is hierarchical: an integral crease nucleates at smaller magnitudes of the prebiased curvature, while a crease decorated with point-like defects emerges at larger curvatures. The phase-space of morphogenesis is traceable to the competition between stretch and bending energies and is parameterizable as a function of the geometry. The framework for generating repetitive torque-dense actuation from slender light-powered actuators holds broader implications for the design of soft, remotely operated machines. Here, it is harnessed in illustrative mechanisms including levers, lifters and grabbers that are powered and regulated exclusively using light.

Network topologies dictate electromechanical coupling in liquid metal-elastomer composites

Elastomers embedded with micro- and nanoscale droplets of liquid metal (LM) alloys like eutectic gallium-indium (EGaIn) can exhibit unique combinations of elastic, thermal, and electrical properties that are difficult to achieve using rigid filler. For composites with sufficient concentrations of liquid metal, the LM droplets can form percolating networks that conduct electricity and deform with the surrounding elastomer as the composite is stretched. Surprisingly, experimental measurements performed on LM-embedded elastomers (LMEEs) show that the total electrical resistance of the composite increases only slightly even as the elastomer is stretched to several times its natural length. In contrast, Pouillet's law would predict an exponential increase in resistance (Ω) with stretch (λ) due to the incompressibility of liquid metal and elastomer. In this manuscript, we perform a computational analysis to examine the unique electromechanical properties of conductive LMEE composites. Our analysis suggests that the gauge factor that quantifies electromechanical coupling (i.e. G = {ΔΩ/Ω0}/λ) decreases with increasing tortuosity of the conductive pathways formed by the connected LM droplets. A dimensionless parameter for path tortuosity can be used to estimate G for statistically homogeneous LMEE composites. These results rationalize experimental observations and provide insight into the influence of liquid metal droplet assembly on the functionality of the composite.

Mechanobiology predicts raft formations triggered by ligand-receptor activity across the cell membrane

Clustering of ligand-binding receptors of different types on thickened isles of the cell membrane, namely lipid rafts, is an experimentally observed phenomenon. Although its influence on cell's response is deeply investigated, the role of the coupling between mechanical processes and multiphysics involving the active receptors and the surrounding lipid membrane during ligand-binding has not yet been understood. Specifically, the focus of this work is on G-protein-coupled receptors (GPCRs), the widest group of transmembrane proteins in animals, which regulate specific cell processes through chemical signalling pathways involving a synergistic balance between the cyclic Adenosine Monophosphate (cAMP) produced by active GPCRs in the intracellular environment and its efflux, mediated by the Multidrug Resistance Proteins (MRPs) transporters. This paper develops a multiphysics approach based on the interplay among energetics, multiscale geometrical changes and mass balance of species, i.e. active GPCRs and MRPs, including diffusion and kinetics of binding and unbinding. Because the obtained energy depends upon both the kinematics and the changes of species densities, balance of mass and of linear momentum are coupled and govern the space-time evolution of the cell membrane. The mechanobiology involving remodelling and change of lipid ordering of the cell membrane allows to predict dynamics of transporters and active receptors -in full agreement with experimentally observed cAMP levels- and how the latter trigger rafts formation and cluster on such sites. Within the current scientific debate on Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2) and on the basis of the ascertained fact that lipid rafts often serve as an entry port for viruses, it is felt that approaches accounting for strong coupling among mechanobiological aspects could even turn helpful in better understanding membrane-mediated phenomena such as COVID-19 virus-cell interaction.

Statistical mechanical analysis of the electromechanical coupling in an electrically-responsive polymer chain

Polymeric materials that couple deformation and electrostatics have the potential for use in soft sensors and actuators with applications ranging from robotic, biomedical, energy, aerospace and automotive technologies. In contrast to the mechanics of polymers that has been studied using statistical mechanics approaches for decades, the coupled response under deformation and electrical field has largely been modeled only phenomenologically at the continuum scale. In this work, we examine the physics of the coupled deformation and electrical response of an electrically-responsive polymer chain using statistical mechanics. We begin with a simple anisotropic model for the electrostatic dipole response to electric field of a single monomer, and use a separation of energy scales between the electrostatic field energy and the induced dipole field energy to reduce the nonlocal and infinite-dimensional statistical averaging to a simpler local finite-dimensional averaging. In this simplified setting, we derive the equations of the most likely monomer orientation density using the maximum term approximation, and a chain free energy is derived using this approximation. These equations are investigated numerically and the results provide insight into the physics of electromechanically coupled elastomer chains. Closed-form approximations are also developed in the limit of small electrical energy with respect to thermal energy; in the limit of small mechanical tension force acting on the chain; and using asymptotic matching for general chain conditions.

Fractional-order nonlinear hereditariness of tendons and ligaments of the human knee

In this paper the authors introduce a nonlinear model of fractional-order hereditariness used to capture experimental data obtained on human tendons of the knee. Creep and relaxation data on fibrous tissues have been obtained and fitted with logarithmic relations that correspond to power-laws with nonlinear dependence of the coefficients. The use of a proper nonlinear transform allows one to use Boltzmann superposition in the transformed variables yielding a fractional-order model for the nonlinear material hereditariness. The fundamental relations among the nonlinear creep and relaxation functions have been established, and the results from the equivalence relations have been contrasted with measures obtained from the experimental data. Numerical experiments introducing polynomial and harmonic stress and strain histories have been reported to assess the provided equivalence relations. This article is part of the theme issue 'Advanced materials modelling via fractional calculus: challenges and perspectives'.

Designing soft pyroelectric and electrocaloric materials using electrets

A temperature variation can electrically polarize a pyroelectric material. In its converse manifestation, the electrocaloric effect entails a change in temperature due to the application of an electric field. These phenomena have wide applications ranging from infrared detection sensors and solid-state refrigeration to energy harvesting. However, the pyroelectric-electrocaloric effect is typically observed in certain classes of hard, brittle crystalline materials that must satisfy a stringent set of lattice symmetry conditions. Some limited experiments have however demonstrated that embedding immobile charges and dipoles in soft foams (thus creating an electret state) may lead to a pyroelectric-like response as well as large deformations desired from soft matter. In this work, we develop a systematic theory for coupled electrical, thermal and mechanical responses of soft electrets. Using simple illustrative examples, we derive closed-form explicit expressions for the pyroelectric and electrocaloric coefficients of electrets. While pyroelectricity in electrets has been noted before, our derived expressions provide a clear quantitative basis to interpret (and eventually design) this effect as well as insights into how the geometrically nonlinear deformation and Maxwell stress give rise to its emergence. We present conditions to obtain a larger pyroelectric and electrocaloric response. In particular, the electrocaloric effect is predicted for the first time in such materials and we show that a proper design and a reasonable choice of materials can lead to a temperature reduction of as much as 1.5 K under the application of electrical fields of 10 MV cm-1.

Computing the Diamagnetic Susceptibility and Diamagnetic Anisotropy of Membrane Proteins from Structural Subunits

The behavior of large, complex molecules in the presence of magnetic fields is experimentally challenging to measure and computationally intensive to predict. This work proposes a novel, mixed-methods approach for efficiently computing the principal magnetic susceptibilities and diamagnetic anisotropy of membrane proteins. The hierarchical primary (amino acid), secondary (α helical and β sheet), and tertiary (α helix and β barrel) structure of transmembrane proteins enables analysis of a complex molecule using discrete subunits of varying size and resolution. The proposed method converts the magnetic susceptibility tensor for all protein subunits to a unit coordinate system and sums them to build the magnetic susceptibility tensor for the membrane protein. Using this approach, we calculate the diamagnetic anisotropy for all transmembrane proteins of known structure and investigate the effect of different subunit resolutions on the resulting predictions of diamagnetic anisotropy. We demonstrate that amino acid residues with aromatic side groups exhibit higher diamagnetic anisotropies. On average, high percentages of aromatic amino acid subunits, a β barrel tertiary structure, and a small volume are correlated with high volumetric diamagnetic anisotropy. Finally, we demonstrate that accounting for the spatial position of the residues with respect to one another is critical to accurately computing the magnetic properties of the complex protein molecule.

Steering with light: indexable photomotility in liquid crystalline polymers

Harnessing light to achieve manipulation and motility in meso and mm-scale systems offers the ability to remotely trigger actuation without requiring on-board power.

Coherent Motion of Monolayer Sheets under Confinement and Its Pathological Implications

Coherent angular rotation of epithelial cells is thought to contribute to many vital physiological processes including tissue morphogenesis and glandular formation. However, factors regulating this motion, and the implications of this motion if perturbed, remain incompletely understood. In the current study, we address these questions using a cell-center based model in which cells are polarized, motile, and interact with the neighboring cells via harmonic forces. We demonstrate that, a simple evolution rule in which the polarization of any cell tends to orient with its velocity vector can induce coherent motion in geometrically confined environments. In addition to recapitulating coherent rotational motion observed in experiments, our results also show the presence of radial movements and tissue behavior that can vary between solid-like and fluid-like. We show that the pattern of coherent motion is dictated by the combination of different physical parameters including number density, cell motility, system size, bulk cell stiffness and stiffness of cell-cell adhesions. We further observe that perturbations in the form of cell division can induce a reversal in the direction of motion when cell division occurs synchronously. Moreover, when the confinement is removed, we see that the existing coherent motion leads to cell scattering, with bulk cell stiffness and stiffness of cell-cell contacts dictating the invasion pattern. In summary, our study provides an in-depth understanding of the origin of coherent rotation in confined tissues, and extracts useful insights into the influence of various physical parameters on the pattern of such movements.

Epithelial and endothelial cells that line various cavities and the vasculature in our bodies, are tightly connected to each other and exist as sheets. Upon confinement in two-dimensional geometries, these cells exhibit rotational motion, which has also been observed in vivo and implicated in physiological processes. However, how this rotational motion is achieved remains unclear. We show that a simple rule wherein preferred direction of motion (i.e., polarization) of cells tends to align with the direction of their velocity is sufficient to induce such coherent movement in confined geometries. We also show that the number of cells within the confinement, the size of the tissue, cell motility and physical properties of the cell and cell-cell connections regulate this coherent motion, and the pattern of invasion when the confinement is relaxed.

Fractional hereditariness of lipid membranes: Instabilities and linearized evolution

In this work lipid ordering phase changes arising in planar membrane bilayers is investigated both accounting for elasticity alone and for effective viscoelastic response of such assemblies. The mechanical response of such membranes is studied by minimizing the Gibbs free energy which penalizes perturbations of the changes of areal stretch and their gradients only (Deseri and Zurlo, 2013). As material instabilities arise whenever areal stretches characterizing homogeneous configurations lie inside the spinoidal zone of the free energy density, bifurcations from such configurations are shown to occur as oscillatory perturbations of the in-plane displacement. Experimental observations (Espinosa et al., 2011) show a power-law in-plane viscous behavior of lipid structures allowing for an effective viscoelastic behavior of lipid membranes, which falls in the framework of Fractional Hereditariness. A suitable generalization of the variational principle invoked for the elasticity is applied in this case, and the corresponding Euler-Lagrange equation is found together with a set of boundary and initial conditions. Separation of variables allows for showing how Fractional Hereditariness owes bifurcated modes with a larger number of spatial oscillations than the corresponding elastic analog. Indeed, the available range of areal stresses for material instabilities is found to increase with respect to the purely elastic case. Nevertheless, the time evolution of the perturbations solving the Euler-Lagrange equation above exhibits time-decay and the large number of spatial oscillation slowly relaxes, thereby keeping the features of a long-tail type time-response.

A frequency-based hypothesis for mechanically targeting and selectively attacking cancer cells

Experimental studies recently performed on single cancer and healthy cells have demonstrated that the former are about 70% softer than the latter, regardless of the cell lines and the measurement technique used for determining the mechanical properties. At least in principle, the difference in cell stiffness might thus be exploited to create mechanical-based targeting strategies for discriminating neoplastic transformations within human cell populations and for designing innovative complementary tools to cell-specific molecular tumour markers, leading to possible applications in the diagnosis and treatment of cancer diseases. With the aim of characterizing and gaining insight into the overall frequency response of single-cell systems to mechanical stimuli (typically low-intensity therapeutic ultrasound), a generalized viscoelastic paradigm, combining classical and spring-pot-based models, is introduced for modelling this problem by neglecting the cascade of mechanobiological events involving the cell nucleus, cytoskeleton, elastic membrane and cytosol. Theoretical results show that differences in stiffness, experimentally observed ex vivo and in vitro, allow healthy and cancer cells to be discriminated, by highlighting frequencies (from tens to hundreds of kilohertz) associated with resonance-like phenomena—prevailing on thermal fluctuations—that could be helpful in targeting and selectively attacking tumour cells.