Hello, everyone:<div><br><div>I was trying to do some excited state calculations involving some charge transfer states using the LR-TDDFPT method in cp2k2022.1. </div><div><br></div><div>Actually, I have tried this before using cp2k7.1. </div><div><a href="https://groups.google.com/g/cp2k/c/PcPUy6FouGE/m/CKOUiGxUDAAJ">https://groups.google.com/g/cp2k/c/PcPUy6FouGE/m/CKOUiGxUDAAJ</a></div><div>Even though all global hybrid functionals such as PBE0-TC and B3LYP-TC can give results comparable to the G09 program, the results obtained by some well known range separated functionals such as wB97XD-TC with MIX_CL_TRUNC are very problematic and the excited state energy difference between the results given by G09 and cp2k7.1 can be large than 1 eV. Such problem remain the same even in cp2k9.0.</div><div><br></div><div>Recently, I have updated the cp2k version from cp2k7.1 to cp2k2022.1 and I tried the same calculation using the same input as before. Luckily, the newest version of cp2k2022.1 seems give excitation energies close to those in G09, which excites me a lot for a period of time. However, when I do some tests on a donor-acceptor system composed of ZnPc and C60, I found that the lowest CT states corresponds to S1 state predicted by G09 is not the case in cp2k2022.1, in which the CT states have much higher energies than the LE states.</div><div><br></div><div>Since the MIX_CL_TRUNC interaction potential is used for range separated functional calculation, I guess the small CUTOFF_RADIUS in the HF section of about 5.0 ang might be the reason responsible for such inconsistency. So I increased the CUTOFF_RADIUS from original 5.0 ang to 10.0 ang. Unfortunately, whenever I use a larger CUTOFF_RADIUS for the truncated calculation, the convergence problem arises and I have tried reading a pre-converged wavefunction, tried to increase the CUTOFF in GPW calculation, switched to the OT algorithm, etc. However, none of them works  and the convergence remain fail.</div><div><br></div><div>I wonder whether anyone can give me some advices on the convergence problem or the reason for correcting such inconsitency between G09 and cp2k2022.1 since it is important to correcly predict the charge transfer state in organic solar cells.</div><div><br></div><div>Thanks a lot.</div><div>PS: I have attached my input file and coordinate file for LR-TD-wB97XD caulations.</div><div><br></div><div><br></div><div><br></div></div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/6e6e41e2-dd2c-4f15-969a-021379ba547en%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/6e6e41e2-dd2c-4f15-969a-021379ba547en%40googlegroups.com</a>.<br />