According to cp2k documentation <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/SUBSYS/KIND.html#POTENTIAL">https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/SUBSYS/KIND.html#POTENTIAL</a> it is possible to set POTENTIAL variable value as ECP, in order to use Effective core potential. I've converted basis set from http://www.tc.uni-koeln.de/PP/clickpse.en.html according to CP2K documentation. The section KIND in my input file is:<div><i> &KIND V<br>BASIS_SET STOLL-PBE<br>POTENTIAL ECP<br>&POTENTIAL<br>V nelec 10<br>V ul<br>2      1.000000000            0.000000000      <br>V S<br>2     14.490000000          178.447971000      <br>2      6.524000000           19.831375000      <br>V P<br>2     14.300000000          109.529763000      <br>2      6.021000000           12.570310000      <br>V D<br>2     17.480000000          -19.219657000      <br>2      5.709000000           -0.642775000      <br>&END<br>&END KIND</i><br></div><div>The basis set in new file is:</div><div><i>V STOLL-PBE<br> 4<br> 2  0  0  7  1  0 # Gauss type functions l=0 set 1 (3s orbital)<br>  12.8432080  0.460448<br>  11.3757530 -0.495309<br>   5.4069740 -0.350821<br>   1.4659270  0.654979<br>   0.5980800  0.582548<br>   0.0887900  0.017081<br>   0.0353180 -0.007455<br> 2  0  0  7  1  0 </i><i># Gauss type functions l=0 set 2 (4s orbital)</i><i><br>  12.8432080  0.126932<br>  11.3757530 -0.139770<br>  5.4069740  -0.089862<br>  1.4659270   0.204669<br>  0.5980800   0.321169<br>  0.0887900  -0.657043<br>  0.0353180  -0.474893<br> 2  1  1  6  0  1 </i><i># Gauss type functions l=1 set 1 (3p orbital)</i><i><br>  31.8898680  0.004658<br>   8.2371780 -0.119249<br>   4.3283730  0.125195<br>   1.5405260  0.620704<br>   0.5280810  0.416780<br>   0.0899620  0.014015<br> 2  2  2  5  1 </i><i># Gauss type functions l=2 set 1 (3d orbital)</i><i><br>  22.6804330  0.031628<br>   6.8613120  0.153807<br>   2.2754450  0.378661<br>   0.7319220  0.507276<br>   0.2007460  0.275921</i><br></div><div><br></div><div>CP2K starts normally but exit with error message:</div><div> <i>CPASSERT failed   qs_mo_types.F:295</i><br></div><div>By checking source code I have error here:</div><div> <i>CPASSERT(nmo >= nomo)</i><br></div><div><br></div><div>Could someone kindly give a working example for force calculations with ECP potential.</div><div><br></div><div>Many thx in advance,<br>Eugene</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/a65b2506-cc62-475b-8514-3ad646aea6a3n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/a65b2506-cc62-475b-8514-3ad646aea6a3n%40googlegroups.com</a>.<br />