[CP2K-user] [CP2K:20237] aimd-1

贾建峰 jjf_sxnu at 163.com
Sun May 26 07:36:38 UTC 2024


Dear CP2K experts,
 I want ot use cp2k to do a metadynamics calculation to get the free energy change for a reaction. I use the coordination number of a  O atom to all H atoms as CV (this O-H bond will be formed in the reaction) to do the metadynamic calculation.  
the CV was defined as:
         &COLVAR
          &COORDINATION
             ATOMS_FROM  61
             KINDS_TO H
             R0 [angstrom] 0.93
             NN 8
             ND 14
          &END COORDINATION
          &END COLVAR


For metadynamics:


    &METADYN
      DO_HILLS
      NT_HILLS 80
      WW 3.0e-3


      &METAVAR
        SCALE 0.2
        COLVAR 1
      &END METAVAR


When the job was ended, I get the fes.dat file by using the command: graph.psmp -cp2k -ndim 1 -ndw 1 -file aimd-1.restart. This file looks like:
       -0.5800822901       -0.0005584043
       -0.5464051991       -0.0009289081
       -0.5127281080       -0.0015031205
       -0.4790510169       -0.0023662732
       -0.4453739259       -0.0036245468
       -0.4116968348       -0.0054030940
       -0.3780197438       -0.0078403220
       -0.3443426527       -0.0110778295
       -0.3106655617       -0.0152462075
........


My question is why the CV vaules in the first column can be negative?  My CV is the coordination number. I think it should be positive anyway.  The whole input file and corresponding xyz file and some output file were attached in this letter. Please help me.


Jianfeng Jia





-- 
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/e541e0e.12cc.18fb3d47c11.Coremail.jjf_sxnu%40163.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20240526/788ca3e5/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: aimd-1.restart
Type: application/octet-stream
Size: 71721 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20240526/788ca3e5/attachment-0006.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: aimd-COLVAR.metadynLog
Type: application/octet-stream
Size: 1640082 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20240526/788ca3e5/attachment-0007.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: aimd-HILLS.metadynLog
Type: application/octet-stream
Size: 13000 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20240526/788ca3e5/attachment-0008.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: fes.dat
Type: application/octet-stream
Size: 4100 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20240526/788ca3e5/attachment-0009.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: md.inp
Type: application/octet-stream
Size: 4578 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20240526/788ca3e5/attachment-0010.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: xyz
Type: application/octet-stream
Size: 14105 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20240526/788ca3e5/attachment-0011.obj>


More information about the CP2K-user mailing list