Shape distributions in simulations are not due to fluctuations around local minima.

A: The evolution of three simulations with the same parameters (pulling force λF = 20 and γ(c, M) = 15) and three in vitro time lapses in the L2/L1 − L3/L1 space. The simulations are performed for a prolonged time of 500,000 Monte Carlo steps. The thin solid lines show the L2/L1 and L3/L1 values every 5,000 Monte Carlo steps. The thick solid lines indicate values averaged over 25 measurements (each 500 Monte Carlo steps apart). The dots highlight the values at 100,000 Monte Carlo steps, which is the simulation time in all other figures. The dashed lines represent the shapes of gastruloids measured by in vitro time-lapse microscopy from 72 h to 94 h, 91 h 20 m and 88 h 40 m respectively for the orange, golden and brown curves. The gastruloid shapes were obtained manually. The thin lines show the raw data (each 40 minutes apart), and the thick lines show time averages over 200 minutes. B: Distribution of the L2/L1 and L3/L1 values for 100 simulations with the same parameters at different time points, up to 500,000 Monte Carlo steps. The values at 100,000 Monte Carlo steps are indicated with dots. C: Distribution of the L2/L1 and L3/L1 values for 100 simulations with the same parameters at different time points without pulling, up to 50,000 Monte Carlo steps. D: Same as B, but with two cell types and parameters γ(1, 2) = −4 and λF = 20 as in (Fig 4D). E: Same as C, but with two cell types and parameter γ(1, 2) = −4. (EPS)