[CP2K-user] [CP2K:21380] Re: Struggling to converge single atom calculation

Frederick Stein f.stein at hzdr.de
Thu Apr 17 11:42:02 UTC 2025


Dear Harry,
Consider them an inspiration as wB97X-2 is a double-hybrid functional, i.e. 
it has an additional MP2 correlation contribution and requires a rather 
different amount of HF exchange. Read the paper defining wB97X (and 
potentially also wB97X-2) how to set up the parameters.
HTH,
Frederick

Harry Richardson schrieb am Donnerstag, 17. April 2025 um 13:24:50 UTC+2:

> Thank you very much for your response Frederick, 
>
> I see there are some xc_sections that use the omega97X-D functional, 
> Do you think either:
> WB97X_2LP.sec 
> <https://github.com/cp2k/cp2k/blob/master/data/xc_section/WB97X_2LP.sec> 
> or WB97X_2TQZ.sec 
> <https://github.com/cp2k/cp2k/blob/master/data/xc_section/WB97X_2TQZ.sec>
> may be suitable, or would it be better to define it myself. 
>
>
> Thanks again 
> Harry 
> On Thursday, April 17, 2025 at 12:04:05 PM UTC+1 Frederick Stein wrote:
>
>> Dear Harry,
>> You have to setup a suitable HF section (see 
>> https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT/XC/HF.html). 
>> CP2K will not do it by itself except of an available shortcut as stated in 
>> the manual (see 
>> https://manual.cp2k.org/trunk/CP2K_INPUT/FORCE_EVAL/DFT/XC/XC_FUNCTIONAL.html). 
>> Compare the already available xc_sections in the CP2K data directory (
>> https://github.com/cp2k/cp2k/tree/master/data/xc_section) for their 
>> setup.
>> Best,
>> Frederick
>>
>> Harry Richardson schrieb am Donnerstag, 17. April 2025 um 11:15:16 UTC+2:
>>
>>> i have been suggested to use multiplicity 3, I have tried and does not 
>>> affect convergence, please see attached. 
>>>
>>> On Thursday, April 17, 2025 at 9:36:10 AM UTC+1 Harry Richardson wrote:
>>>
>>>> Hi all, I am pretty new to DFT and trying to run some single atom (C) 
>>>> DFT to use as a subtraction when fine tuning a machine learned force field 
>>>> for MD. As such it needs to be run at the same level of theory as my full 
>>>> system calculations. 
>>>>
>>>> I just need single point energies. 
>>>>
>>>> At the moment it is not converging, I have attached my input and output 
>>>> files, i am using the omega97X-D functional. 
>>>>
>>>> Some thing I have tried include: 
>>>>
>>>>    - Changing the box from periodic to non periodic
>>>>    - Reducing mixing
>>>>    - Increasing & decreasing the size of the box 
>>>>    - Changing multiplicity 
>>>>    - Reducing cutoff 
>>>>    
>>>> The energy change per step seems to reduce over steps but does not 
>>>> converge. 
>>>>
>>>> I am using the DZVP-MOLOPT basis set and GTH-PBE-q4 potential as 
>>>> suggested in here: https://www.cp2k.org/exercises:2018_uzh_cmest:faq 
>>>> potentially that is the issue? 
>>>>
>>>> Would appreciate any help you can provide 
>>>>
>>>> Thanks
>>>>
>>>

-- 
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 visit https://groups.google.com/d/msgid/cp2k/a1841aef-5cc3-437d-bd75-a5f898849ac0n%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20250417/34989a07/attachment.htm>


More information about the CP2K-user mailing list