“…In dipole-coupled molecular arrays, it is necessary to take into account the dipole–dipole interactions, which can be written as the following form , V i , j = μ i · μ j − 3 ( bold-italicμ i · n ) · ( bold-italicμ boldj · n ) | r i − r j | 3 Here, r i ( j ) denotes the position vector of the i ( j )th molecule, and n represents a unit vector along r i , j with r i , j = r j – r i . In addition, following refs and , the vector r i , j is assumed to be orthogonal to the static electric field in the rest of this paper. Consequently, for N dipole–dipole-coupled ultracold polar molecules in an electric field, the whole systemic Hamiltonian reads H = ∑ i = 1 N H i + ∑ i = 1 N − 1 ∑ j = i + 1 N V i , j …”