[CP2K-user] Ionic crystals with GFN-xTB
mpol
polynsk... at gmail.com
Thu Jul 4 08:57:12 UTC 2019
Dear CP2k developers and users,
Have anyone encountered a problem with GFN-xTB when modeling ionic
crystals? If I set
CHECK_ATOMIC_CHARGES .TRUE.
When performing CELL_OPT of a CaF_2 supercell, the program stops with the
following warning:
*** WARNING in xtb_matrices.F:1088 :: Atomic charges outside chemical ***
*** range were detected. Switch-off CHECK_ATOMIC_CHARGES keyword in the ***
*** &xTB section if you want to force to continue the calculation. ***
Setting CHECK_ATOMIC_CHARGES .FALSE. makes the calculation divergent. As I
get it, the charges on Ca and F become more and more unphysical in every
SCF procedure, and, concomitantly, the total energy becomes lower and
lower. I have attached the input and output of the divergent calculation to
this message.
Could anyone give a hint on what to do in this situation? I am interested
in MD of solid-liquid interfaces with ionic crystals. GFN-xTB seems the
best choice.
With best regards,
Mikhail
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20190704/4e6a2af2/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test1.inp
Type: chemical/x-gamess-input
Size: 2604 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20190704/4e6a2af2/attachment.inp>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: opt.cif
Type: chemical/x-cif
Size: 346 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20190704/4e6a2af2/attachment.cif>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: opt.xyz
Type: chemical/x-xyz
Size: 4715 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20190704/4e6a2af2/attachment.xyz>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test1.out
Type: application/octet-stream
Size: 391815 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20190704/4e6a2af2/attachment.obj>
More information about the CP2K-user
mailing list