Traction force microscopy for linear and nonlinear elastic materials as a parameter identification inverse problem
Gesa Sarnighausen, Tram Thi Ngoc Nguyen, Thorsten Hohage, Mangalika Sinha, Sarah Koester, Timo Betz, Ulrich Sebastian Schwarz, Anne Wald
TL;DR
This work formulates the traction force microscopy inverse problem for linear 2.5D and nonlinear pure 2D substrate models within a rigorous elasticity framework, defining forward maps $\hat{A}$ and $\hat{S}$ that relate traction stresses and force densities to substrate displacements. It proves well-posedness, derives the Fréchet derivative and its adjoint via energy methods and Gelfand triples, and introduces a polyconvex, isotropic material law to guarantee existence and stability of solutions in the nonlinear case. Numerical results on simulated and experimental data demonstrate stable, regularized reconstructions using CGNE for linear problems and truncated Newton-CG for nonlinear problems, with artifacts reduced by appropriate $H^1$-type penalties and homotopy strategies. The approach provides a robust mathematical basis for extending TFM to nonlinear substrate materials and for systematic comparison with standard physics-based methods like FTTC, with ready-to-use software tools for practitioners.
Abstract
Traction force microscopy is a method widely used in biophysics and cell biology to determine forces that biological cells apply to their environment. In the experiment, the cells adhere to a soft elastic substrate, which is then deformed in response to cellular traction forces. The inverse problem consists in computing the traction stress applied by the cell from microscopy measurements of the substrate deformations. In this work, we consider a linear model, in which 3D forces are applied at a 2D interface, called 2.5D traction force microscopy, and a nonlinear pure 2D model, from which we directly obtain a linear pure 2D model. All models lead to a linear resp. nonlinear parameter identification problem for a boundary value problem of elasticity. We analyze the respective forward operators and conclude with some numerical experiments for simulated and experimental data.
