where N is the number of monomers, m is the monomer number within a given simulation, α and β represent the Cartesian x, y and z axes, delta is the Kronecker delta, a is a vector that describes the molecular long axis, which is computed for each monomer as the eigenvector associated with the smallest eigenvalue of the inertia tensor. The director at each frame was defined as the eigenvector associated with the largest eigenvalue of the ordering tensor. The order parameter <P2> corresponds to the largest eigenvalue of 𝑄 𝛼𝛽 . For polar/ferroelectric nematics we introduce the polar order parameter, <P1>, which is calculated according to Eq. 8