<div dir="ltr">Hi,<div><br></div><div>Is it possible for CP2K to output the transition dipole moment <phi_A | r | phi_B>  between two KS states, say HOMO to LUMO transition? </div><div>I tried to calculate the integral by printing the MO_CUBE files which are the phi_A and phi_B and then do the integration. This is fine for </div><div>molecules or in other word non-periodic systems. But for periodic systems, the result was not as expected, maybe due to the position </div><div>operator is ill-defined in that case. So how to evaluate the transition dipole integral for periodic systems with the MO_CUBE files known?</div><div><br></div><div>Best,</div><div>Xiaoming   </div></div>