“…The order of each term is defined as X r + X v À 2. A similar construction is performed for the dipole moment (l) and polarizability (a) operators that are used for absorption and Raman intensity calculations, respectively [18,19,6,8]. We do not develop here the technical details of these construction any further, since it is described in various papers.…”