-- dipole operator for -- light polarized in the -- {ex,ey,ez} direction NF = 16 NB = 0 IndexDn_2p = {0,2,4} IndexUp_2p = {1,3,5} IndexDn_3d = {6,8,10,12,14} IndexUp_3d = {7,9,11,13,15} ex = 1 ey = 2 ez = 3 exn = ex / math.sqrt(ex^2 + ey^2 + ez^2) eyn = ey / math.sqrt(ex^2 + ey^2 + ez^2) ezn = ez / math.sqrt(ex^2 + ey^2 + ez^2) Akm = {{1, -1,( exn + I * eyn)/math.sqrt(2)}, {1, 0, ezn}, {1, 1,(-exn + I * eyn)/math.sqrt(2)}} OppT = NewOperator("CF", NF, IndexUp_3d, IndexDn_3d, IndexUp_2p, IndexDn_2p, Akm)