[CP2K-user] [CP2K:17228] Best practices for large non periodic molecules
Marco Severi
marcoseveri.c at gmail.com
Fri Jun 24 12:31:42 UTC 2022
Thank you very much for your comment. I saw in the documentation the
constraint about the shape of the cell for the wavelet solver and was
thinking about trying different solvers and cell shapes, now I am sure that
it was not a bad idea. Thanks again!
Marco
Il giorno giovedì 23 giugno 2022 alle 16:24:20 UTC+2 mattwa... at gmail.com ha
scritto:
> Quick note on the POISSON solver:
> The wavelet solver you have chosen is great because it only needs the
> density to drop to zero at the box edges - so size of molecule + about 10A.
> However, it only works for cubic cells so it needs to be largest molecule
> dimension + 10A in all 3 directions. This could be inefficient for very
> long or disc shaped molecules.
> Alternatively the MT or analytic solvers can be used but in each direction
> they should be at least twice the minimum box size (I believe). This could
> still be beneficial in some cases compared to above...
> Matt
>
> On Wednesday, 22 June 2022 at 16:01:28 UTC+1 marcos... at gmail.com wrote:
>
>> Thank you very much, I will try as soon as possible.
>>
>> Best,
>> Marco
>>
>> Il giorno mercoledì 22 giugno 2022 alle 16:45:14 UTC+2 hehn... at gmail.com
>> ha scritto:
>>
>>> Dear Marco,
>>>
>>> regarding convergence I would recommend to use the OUTER_SCF keyword and
>>> SCF_GUESS RESTART, giving a relatively small number of inner SCF cycles,
>>> maybe 20 to 50. Also choosing DIIS as minimizer can help to speed up
>>> convergence.
>>>
>>> Best regards,
>>>
>>> Anna
>>>
>>> Am Mi., 22. Juni 2022 um 14:11 Uhr schrieb Marco Severi <
>>> marcos... at gmail.com>:
>>>
>>>> Good morning,
>>>>
>>>> I am a fairly new CP2K user and I am trying to perform some
>>>> calculations on a non-periodic molecular system of 657 atoms.
>>>>
>>>> I am writing because I do not fully understand which keywords I should
>>>> use and I have the feeling I am doing some rookie mistakes.
>>>>
>>>> At the moment I am trying to do some DFT calculations with the PBE
>>>> functional and the dispersion correction. The molecule is about 35 Angstrom
>>>> long and so I chose a 50x50x50 cell. The system is organometallic, contains
>>>> 8 Fe atoms and the BF4- counterions. I tried to use the DZVP basis set but
>>>> I quickly reached 100GB of occupied memory, is it normal? I am using the
>>>> non "PERIODIC NONE" and "PSOLVER WAVELET" keywords and the FULL_KINETIC
>>>> preconditioner. I am using 5 grids with cutoff equal to 400 and rel_cutoff
>>>> equal to 60. I tried to maintain all the other options as “default” as
>>>> possible. I enclose the input and the output of the calculation. (I am
>>>> aware that the attached example did not converge and that I can change the
>>>> number of SCF iterations)
>>>>
>>>> I would like to know if those options are adequate for the type of
>>>> calculation I am performing and if there is room for any improvement in the
>>>> efficiency or in the convergence, since this calculation seems to converge
>>>> particularly slowly. Any help or reference will be greatly appreciated.
>>>>
>>>> Thanks in advance for your time and patience.
>>>>
>>>> Kind regards,
>>>>
>>>> Marco Severi
>>>>
>>>> --
>>>> You received this message because you are subscribed to the Google
>>>> Groups "cp2k" group.
>>>> To unsubscribe from this group and stop receiving emails from it, send
>>>> an email to cp2k+uns... at googlegroups.com.
>>>> To view this discussion on the web visit
>>>> https://groups.google.com/d/msgid/cp2k/838e8bc8-a443-4cd2-a8fa-b591ddf1d8ffn%40googlegroups.com
>>>> <https://groups.google.com/d/msgid/cp2k/838e8bc8-a443-4cd2-a8fa-b591ddf1d8ffn%40googlegroups.com?utm_medium=email&utm_source=footer>
>>>> .
>>>>
>>>
--
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/72d33217-854a-4338-be05-9b651e51bd96n%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20220624/020ce574/attachment-0001.htm>
More information about the CP2K-user
mailing list