the matrix elements can be found in Ref and are expressed using 3j-symbols. The effective Hamiltonians in ( ) read (7) where is the reduced molecular mass and are the electronic potential energy surfaces. The induced dipole moment, , of each molecule is evaluated and the macroscopic polarization is computed following the mean field approximation in accordance to (8) where is the number density of molecules. It is coupled back to Maxwell's equations ( ) via polarization current. The complete model employed in this work relies on the system of coupled partial differential equations ( ) and ( ). We note that equations ( ) are driven by a local electric field, which in principle depends sensitively on spatial coordinates. We thus need to propagate quantum equations at every grid point and couple their dynamics to Maxwell's equations.