Linear and nonlinear mechanical responses can be quite different in models for biological tissues.

The fluidity of biological tissues - whether cells can change neighbors and rearrange - is important for their function. In traditional materials, researchers have used linear response functions, such as the shear modulus, to accurately predict whether a material will behave as a fluid. Similarly, in disordered 2D vertex models for confluent biological tissues, the shear modulus becomes zero precisely when the cells can change neighbors and the tissue fluidizes, at a critical value of control parameter s0* = 3.81. However, the ordered ground states of 2D vertex models become linearly unstable at a lower value of control parameter (3.72), suggesting that there may be a decoupling between linear and nonlinear response. We demonstrate that the linear response does not correctly predict the nonlinear behavior in these systems: when the control parameter is between 3.72 and 3.81, cells cannot freely change neighbors even though the shear modulus is zero. These results highlight that the linear response of vertex models should not be expected to generically predict their rheology. We develop a simple geometric ansatz that correctly predicts the nonlinear response, which may serve as a framework for making nonlinear predictions in other vertex-like models.