Measuring the viscoelastic relaxation function of cells with a time-dependent interpretation of the Hertz-Sneddon indentation model.

The Hertz-Sneddon elastic indentation model is widely adopted in the biomechanical investigation of living cells and other soft materials using atomic force microscopy despite the explicit viscoelastic nature of these materials. In this work, we demonstrate that an exact analytical viscoelastic force model for power-law materials, can be interpreted as a time-dependent Hertz-Sneddon-like model. Characterizing fibroblasts (L929) and osteoblasts (OFCOLII) demonstrates the model's accuracy. Our results show that the difference between Young's modulus EY obtained by fitting force curves with the Hertz-Sneddon model and the effective Young's modulus derived from the viscoelastic force model is less than 3%, even when cells are probed at large forces where nonlinear deformation effects become significant. We also propose a measurement protocol that involves probing samples at different indentation speeds and forces, enabling the construction of the average viscoelastic relaxation function of samples by conveniently fitting the force curves with the Hertz-Sneddon model.

Sublinear drag regime at mesoscopic scales in viscoelastic materials.

Stressed soft materials commonly present viscoelastic signatures in the form of power-law or exponential decay. Although exponential responses are the most common, power-law time dependencies arise peculiarly in complex soft materials such as living cells. Understanding the microscale mechanisms that drive rheologic behaviors at the macroscale shall be transformative in fields such as material design and bioengineering. Using an elastic network model of macromolecules immersed in a viscous fluid, we numerically reproduce those characteristic viscoelastic relaxations and show how the microscopic interactions determine the rheologic response. The macromolecules, represented by particles in the network, interact with neighbors through a spring constant k and with fluid through a non-linear drag regime. The dissipative force is given by γvα, where v is the particle's velocity, and γ and α are mesoscopic parameters. Physically, the sublinear regime of the drag forces is related to micro-deformations of the macromolecules, while α ≥ 1 represents rigid cases. We obtain exponential or power-law relaxations or a transitional behavior between them by changing k, γ, and α. We find that exponential decays are indeed the most common behavior. However, power laws may arise when forces between the macromolecules and the fluid are sublinear. Our findings show that in materials not too soft not too elastic, the rheological responses are entirely controlled by α in the sublinear regime. More specifically, power-law responses arise for 0.3 ⪅ α ⪅ 0.45, while exponential responses for small and large values of α, namely, 0.0 ⪅ α ⪅ 0.2 and 0.55 ⪅ α ⪅ 1.0.

Double power-law viscoelastic relaxation of living cells encodes motility trends.

Living cells are constantly exchanging momentum with their surroundings. So far, there is no consensus regarding how cells respond to such external stimuli, although it reveals much about their internal structures, motility as well as the emergence of disorders. Here, we report that twelve cell lines, ranging from healthy fibroblasts to cancer cells, hold a ubiquitous double power-law viscoelastic relaxation compatible with the fractional Kelvin-Voigt viscoelastic model. Atomic Force Microscopy measurements in time domain were employed to determine the mechanical parameters, namely, the fast and slow relaxation exponents, the crossover timescale between power law regimes, and the cell stiffness. These cell-dependent quantities show strong correlation with their collective migration and invasiveness properties. Beyond that, the crossover timescale sets the fastest timescale for cells to perform their biological functions.

Fracturing highly disordered materials.

We investigate the role of disorder on the fracturing process of heterogeneous materials by means of a two-dimensional fuse network model. Our results in the extreme disorder limit reveal that the backbone of the fracture at collapse, namely, the subset of the largest fracture that effectively halts the global current, has a fractal dimension of 1.22 ± 0.01. This exponent value is compatible with the universality class of several other physical models, including optimal paths under strong disorder, disordered polymers, watersheds and optimal path cracks on uncorrelated substrates, hulls of explosive percolation clusters, and strands of invasion percolation fronts. Moreover, we find that the fractal dimension of the largest fracture under extreme disorder, d(f) = 1.86 ± 0.01, is outside the statistical error bar of standard percolation. This discrepancy is due to the appearance of trapped regions or cavities of all sizes that remain intact till the entire collapse of the fuse network, but are always accessible in the case of standard percolation. Finally, we quantify the role of disorder on the structure of the largest cluster, as well as on the backbone of the fracture, in terms of a distinctive transition from weak to strong disorder characterized by a new crossover exponent.