<div dir="ltr">Dear Developers<o:p></o:p> I would like to draw your attention to a possible bug of CP2K. I am trying to make periodic RTP by using hybrid functionals also.<o:p></o:p> However, when I set the keyword PERIODIC = TRUE in the REAL_TIME_PROPAGATION block, the perturbed state that is to be propagated is wrong. <o:p></o:p> Here below I report the electronic energy for an H2 molecule after a Kick (input at the end of this mail). As you can see, the propagated state has energies which differ a lot from the ground state. Furthermore, also the calculated electronic spectrum is very strange.<o:p></o:p> On the other side, when using PERIODIC = FALSE the energy of the propagated state slightly differs from the ground state and the calculated electronic spectrum is the right one.<o:p></o:p> I have to underline that the bug appears just when one uses hybrid functionals. With GGA everything goes fine. Is this a bug or something already known?.<o:p></o:p> <o:p> </o:p> PERIODIC TRUE :<o:p></o:p> Total energy: -1.16466508784488<o:p></o:p> Total energy: -0.76496701712800<o:p></o:p> Total energy: -0.76311027473197<o:p></o:p> Total energy: -0.76310900639003<o:p></o:p> Total energy: -0.76310899786364<o:p></o:p> Total energy: -0.76310899785477<o:p></o:p> Total energy: -0.76325873454457<o:p></o:p> Total energy: -0.76316886201251<o:p></o:p> Total energy: -0.76316802986251<o:p></o:p> Total energy: -0.76316802392354<o:p></o:p> Total energy: -0.76316802386850<o:p></o:p> <o:p> </o:p> <o:p> </o:p> PERIODIC FALSE:<o:p></o:p> GS Total energy: -1.16466508784488<o:p></o:p> Total energy: -1.16466508783571<o:p></o:p> Total energy: -1.16466508783580<o:p></o:p> Total energy: -1.16466508783580<o:p></o:p> Total energy: -1.16466508783579<o:p></o:p> Total energy: -1.16466508783579<o:p></o:p> <o:p> </o:p> ---------------------------------------------------------------------------------<o:p></o:p> INPUT<o:p></o:p> <o:p> </o:p> &FORCE_EVAL<o:p></o:p> METHOD QUICKSTEP<o:p></o:p> &DFT<o:p></o:p> &REAL_TIME_PROPAGATION<o:p></o:p> APPLY_DELTA_PULSE<o:p></o:p> DELTA_PULSE_DIRECTION 1 0 0<o:p></o:p> DELTA_PULSE_SCALE 1.0E-5<o:p></o:p> PERIODIC .TRUE.<o:p></o:p> MAX_ITER 10<o:p></o:p> MAT_EXP TAYLOR<o:p></o:p> PROPAGATOR ETRS<o:p></o:p> EXP_ACCURACY 1.0E-10<o:p></o:p> EPS_ITER 1.0E-8<o:p></o:p> INITIAL_WFN SCF_WFN<o:p></o:p> &END<o:p></o:p> BASIS_SET_FILE_NAME …./data/BASIS_SET<o:p></o:p> POTENTIAL_FILE_NAME …./data/POTENTIAL<o:p></o:p> &MGRID<o:p></o:p> CUTOFF 240<o:p></o:p> &END MGRID<o:p></o:p> &QS<o:p></o:p> EPS_DEFAULT 1.0E-10<o:p></o:p> MAP_CONSISTENT<o:p></o:p> &END QS<o:p></o:p> &SCF<o:p></o:p> MAX_SCF 100<o:p></o:p> EPS_SCF 1.0E-9<o:p></o:p> SCF_GUESS ATOMIC<o:p></o:p> &OUTER_SCF<o:p></o:p> EPS_SCF 1.0E-6<o:p></o:p> MAX_SCF 10<o:p></o:p> &END &END SCF<o:p></o:p> &XC<o:p></o:p> &XC_FUNCTIONAL B3LYP<o:p></o:p> &END XC_FUNCTIONAL<o:p></o:p> <o:p> </o:p> &HF<o:p></o:p> &INTERACTION_POTENTIAL<o:p></o:p> POTENTIAL_TYPE TRUNCATED<o:p></o:p> CUTOFF_RADIUS 5.00<o:p></o:p> T_C_G_DATA …/data/t_c_g.dat<o:p></o:p> &END<o:p></o:p> &END &END XC<o:p></o:p> &END DFT<o:p></o:p> &SUBSYS<o:p></o:p> &CELL<o:p></o:p> ABC 11.00 11.00 11.00<o:p></o:p> &END CELL<o:p></o:p> &COORD<o:p></o:p> H -0.375 0.0 0.0<o:p></o:p> H +0.375 0.0 0.0<o:p></o:p> &END COORD<o:p></o:p> &KIND H<o:p></o:p> BASIS_SET DZV-GTH-PBE<o:p></o:p> POTENTIAL GTH-PBE-q1<o:p></o:p> &END KIND<o:p></o:p> &END SUBSYS<o:p></o:p> &END FORCE_EVAL<o:p></o:p> &GLOBAL<o:p></o:p> PROJECT H2-RTP<o:p></o:p> RUN_TYPE RT_PROPAGATION<o:p></o:p> PRINT_LEVEL MEDIUM<o:p></o:p> &END GLOBAL<o:p></o:p> &MOTION<o:p></o:p> &MD<o:p></o:p> ENSEMBLE NVE<o:p></o:p> STEPS 2<o:p></o:p> TIMESTEP [au_t] 0.25<o:p></o:p> TEMPERATURE 300.0<o:p></o:p> &END MD<o:p></o:p> &END MOTION Luca <o:p></o:p></div>