The crystal order is determined by a modified bond-orientation order
analysis [1]}, [2]}, [3]}. A complex order parameter
\(q_{lm}(i,t) \equiv 1/n_b \Sigma _{j=1}^{n_b}Y_{lm}(\textbf {r}_{\alpha \beta }(t))\) 
