simulated annealing using CP2K

jing... at outlook.com jing... at outlook.com
Sat Aug 20 20:23:18 UTC 2016


Hi, Janos. Recently, I am learning simulated annealing to search the global 
minimum of small molecule. I saw your nice suggestion, but I have some 
questions,
For annealing, how to set the temperature? I tested my input file, when I 
use the keywords temperature and annealing, the error implied *thermostat *is 
needed. 
But thermostat and annealing are against, so I wonder how to set the 
temperature? And  in every input files(heating and cooling), should I do 
NVE simulations or
just do the optimization? 
      Thanks.
       
       Jing

在 2015年7月27日星期一 UTC-7下午5:53:32,Geng Sun写道:
>
> Janos,
>
>     Thanks ! that's very helpful. I will give it a try.
>
> Geng Sun
>
> 2015-07-23 15:57 GMT+08:00 Jano... at googlemail.com <javascript:> <
> jano... at gmail.com <javascript:>>:
>
>> Dear Geng,
>>
>> It can be done by using the EXT_RESTART option and some simple input 
>> files.
>>
>> template: simple input file with EXT_RESTART from projectname-1.restart, 
>> restarting everything, and annealing parameter:
>> ANNEALING ${PAR}
>>
>> Then, you can also use the SET and INCLUDE commands in three different 
>> inputs:
>>
>> heating.inp:
>> @SET PAR 1.000001 #check the number for your needs
>> @INCLUDE ./template
>>
>> const.inp
>> @SET PAR 1.000000
>> @INCLUDE ./template
>>
>> cooling.inp:
>> @SET PAR 0.999999 #check the number for your needs
>> @INCLUDE ./template
>>
>> then you simply write a script:
>> #! /usr/bin/bash
>>
>> for i in {1..99999}; do  #again adjust to your needs
>>   cp2k.popt const.inp >>out
>>   cp2k.popt heating.inp >>out
>>   cp2k.popt const.inp >>out
>>   cp2k.popt cooling.inp >>out
>> done
>>
>> You can set different simulation times by the SET command as well...
>> If you want, you can set different project names for different 
>> temperatures as well, so you can have the different temperature structures 
>> in different xyz files. And it is also a good idea to set the temperature 
>> by @SET for the input files.
>>
>> HTH,
>>
>> Janos
>>
>>
>> On Tuesday, July 14, 2015 at 10:28:37 AM UTC+2, Geng Sun wrote:
>>>
>>> Hello CP2K users 
>>>
>>>       Could anyone give me some hints that how to do fix cell simulated 
>>> annealing to find the stable structures in CP2K ?
>>>
>>>
>>> Thanks in advance !
>>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20160820/c627e919/attachment.htm>


More information about the CP2K-user mailing list