Furthermore, it is evident that strong correlations exist along the z-direction, i.e. even when the overall polarisation is zero, each molecule in a chain along the z-direction will prefer to have similar z-components of their polarisation. It would be of interest to probe the local order in the high temperature phase using diffuse X-ray (or neutron) scattering to further understand the merits of the model Hamiltonian beyond prediction of the phase transition temperature.