The polymers are simulated using LAMMPS in a 2D plane and a Langevin thermostat. The specific parameters used in this model are varied to explore the total landscape of rupturing behaviour and discussed in more detail in each section. Where not otherwise specified, the simulations use an attractive Lennard-Jones well depth of ε HH = 16.5 × 10 -21 J (1200 K k B ),a repulsive Lennard-Jones energy parameter of ε BB = 64.0 × 10 -21 J and an angular energy parameter of K θ = 200 × 10 -21 J. The Langevin thermostat is set in the range 100 K to 300 K with a damping parameter of 1.00 ns. However, rupture of a covalent bond in a single tropocollagen molecule cannot be repaired and must be modelled differently. To account for the effect of bonds breaking, we use fix bond/break feature of LAMMPS. This checks all bonds every N steps to check if any bond lengths are greater than a cutoff distance r c . If a bond is longer than r c , it is removed and the neighbour lists representing bonds and angular bonds are recalculated. In this work, we choose r c to be equal to the cutoff distance for repulsive interactions between body beads, σ BB = 138 nm.