[CP2K-user] [CP2K:17997] Re: QM/MM

Kamil unixrulez.73 at gmail.com
Tue Nov 8 09:44:43 UTC 2022


Dear Marcella,

thank you for help. However, the main problem I encounter is slightly 
different - as I said: my input file specifies temperature of 300K for 
velocities initialization. Nonetheless, QM/MM is not starting from that 
temperature...

 MD_VEL| Velocities initialization
 MD_VEL| Initial temperature [K]                                     
 207.386439
 MD_VEL| COM velocity             0.0000002385    -0.0000001532   
 -0.0000000931

Abovementioned part of the calculations does not change between the first 
(equilibration) run and the consecutive restarts. So, keyword TEMPERATURE 
in input, as well as velocities gathered in *.restart files, are somehow 
skipped by CP2K. Did I forget about sth in the input? 

Best regards,
K

niedziela, 6 listopada 2022 o 10:48:01 UTC+1 Marcella Iannuzzi napisał(a):

> Hi
>
> According to the input that you posted the calculation is restarted from 
> ${FILERESTART}
> Velocities and thermostat should also be restarted from the value in that 
> file, hence the temperature.
> To avoid that, you can use RESTART_VEL FALSE
> Regards
> Marcella
>
>
> On Friday, November 4, 2022 at 7:03:20 PM UTC+1 Kamil wrote:
>
>> Dear developers, 
>>
>> I have a problem regarding QM/MM simulation of biochemical system - after 
>> preparation and equilibration in AMBER I started QM/MM in cp2k, however I 
>> noticed strange behaviour:
>>
>> i) velocities are initialized for the temperature ca. 207 K - it's 
>> strange, because in the input file temperature is set to ca. 300K,
>>
>> ii) for every restart initialization is performed again at 207K...
>>
>> What can be the cause of that behaviour? 
>>
>> P.S. My input file is in the attachment. 
>>
>> Best regards
>>
>>

-- 
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/8f00c27a-d0b8-4ada-a052-afec0f0892f4n%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20221108/308a6e00/attachment.htm>


More information about the CP2K-user mailing list