“…The free space propagator for spherical vector waves (dyadic Green's function) is the solution of the inhomogeneous Helmholtz equation ( g 2 + k2)G(r,r') = -47716(r -r') (Al) in spherical coordinates and can be written with respect to the molecular sites as If we expadd the unit dyadic, 1, with respect to a spherical basis (Rose, 1955) as and assume Ir, -Rap 1 2 rj, the Neumann expansion of G is (Goertzel and Tralli, 1960) Now the addition theorem for scalar Helmholtz harmonics can be applied to h,(k)r, -RnP I)YLM(rP -RaB) to obtain (Nozawa, 1966) where we assumed r, < Rap, and thus define Now we recall that the spherical harmonics and the irreducible spherical tensors are related through the equation (Peterson and Strom, 1973) so that equation (A4) an easy but lengthy calculation, with the help of the formulas of Borghese et al (1980) shows that the H's and the K's are the matrix elements of G with respect to M and N. The other matrix elements of G do not appear in the present work because we deal with solenoidal fields, which require M and N only for their description. Finally, we notice that the above procedure also allows us to define J;~P, , , , and L$~,,, as the matrix elements of G with respect to M and N. It is, in fact, sufficient to assume r, > R a p and consequently substitute in GLjWLM, jA(kRap) to hA(kRaB) and, in M and N, h , , (kr,) to /L,(kra).…”