OK, here is Update. Actually I'm testing ECP with V2O5 crystall 1x3x2 Supercell. Now when I switch to GTH pseudopotential I have the following initial values for number of MO:<div><i>  Total number of            - Atomic kinds:                                   2<br>                             - Atoms:                                         84<br>                             - Shell sets:                                    84<br>                             - Shells:                                       684</i><br>There are 24 V atoms and 60 oxygens. I.e. total number of shells is: (24*13+60*6)=672 (for V there is 13 valence electrons for O it's 6). Which is different then 684 for some reason.</div><div>In case of ECP the numbers are:</div><div>  Total number of            - Atomic kinds:                                   2<br>                             - Atoms:                                         84<br>                             - Shell sets:                                   336<br>                             - Shells:                                       192<br></div><div>Well, Huston we have a problem. How to resolve the issue?</div><div class="gmail_quote"><div dir="auto" class="gmail_attr">понедельник, 7 февраля 2022 г. в 19:30:11 UTC+3, Eugene: <br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">According to cp2k documentation <a href="https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/SUBSYS/KIND.html#POTENTIAL" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=ru&q=https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/SUBSYS/KIND.html%23POTENTIAL&source=gmail&ust=1644344452815000&usg=AFQjCNH-e2exVhsWMXDsM9_zaJaID2mgxQ">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 <a href="http://www.tc.uni-koeln.de/PP/clickpse.en.html" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=ru&q=http://www.tc.uni-koeln.de/PP/clickpse.en.html&source=gmail&ust=1644344452815000&usg=AFQjCNF9hMKqijZHorW-uqGThQZ9t8NERw">http://www.tc.uni-koeln.de/PP/clickpse.en.html</a> 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></blockquote></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/f4fcedd7-4dbe-4403-9d54-f67fee1b9b94n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/f4fcedd7-4dbe-4403-9d54-f67fee1b9b94n%40googlegroups.com</a>.<br />