FES

600 views
Skip to first unread message

T. Salavati fard

unread,
Mar 23, 2014, 4:52:18 PM3/23/14
to cp...@googlegroups.com
Dear all

I am CPMD user. I am doing metadynamics calculations and I need to reconstruct FES. I know there is a code by Teo to reconstruct free energy surface using colvar and parvar files.
My problem is finding FES program source!! I couldn't find that code to compile and try..
I appreciate your helps--Thanks

Ari Paavo Seitsonen

unread,
Mar 23, 2014, 5:39:08 PM3/23/14
to cp...@googlegroups.com
Dear T. Salavati fard,

  In <cp2k-root>/src/metadyn_tools/, if you are using the latest version? It used to be in <cp2k-root>/tools/metadyn/ earlier. One compiles by typing "make fes" in '<cp2k-root>/makefiles'.

    Greetings from Zurich,

       apsi


--
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...@googlegroups.com.
To post to this group, send email to cp...@googlegroups.com.
Visit this group at http://groups.google.com/group/cp2k.
For more options, visit https://groups.google.com/d/optout.



--
-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-
  Ari Paavo Seitsonen / Ari.P.S...@iki.fi / http://www.iki.fi/~apsi/
  Institut für Chemie der Universität Zürich
  Tel: +41 44 63 55 44 97  /  Mobile: +41 79 71 90 935

Ole Schütt

unread,
Mar 23, 2014, 7:00:28 PM3/23/14
to cp...@googlegroups.com
Hi T. Salavati fard,

> One compiles by typing "make fes" in '<cp2k-root>/makefiles'.

Actually, if you are using the latest svn version you now have to type
"make graph", because the names of the make-targets now match the names
of the executables. Alternatively, you can just call make without any
targets, which will build all executables including the fes tool.


-Ole



On 2014-03-23 22:39, Ari Paavo Seitsonen wrote:
> Dear T. Salavati fard,
>
>   In <cp2k-root>/src/metadyn_tools/, if you are using the latest
> version? It used to be in <cp2k-root>/tools/metadyn/ earlier. One
> compiles by typing "make fes" in '<cp2k-root>/makefiles'.
>
>     Greetings from Zurich,
>
>        apsi
>
> 2014-03-23 21:52 GMT+01:00 T. Salavati fard
> <t.salav...@gmail.com>:
>
>> Dear all
>>
>> I am CPMD user. I am doing metadynamics calculations and I need to
>> reconstruct FES. I know there is a code by Teo to reconstruct free
>> energy surface using colvar and parvar files.
>> My problem is finding FES program source!! I couldn't find that code
>> to compile and try..
>> I appreciate your helps--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+uns...@googlegroups.com.
>> To post to this group, send email to cp...@googlegroups.com.
>> Visit this group at http://groups.google.com/group/cp2k [1].
>> For more options, visit https://groups.google.com/d/optout [2].
>
> --
>
> -=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-=*=-
>   Ari Paavo Seitsonen / Ari.P.S...@iki.fi /
> http://www.iki.fi/~apsi/ [3]
>   Institut für Chemie der Universität Zürich
>   Tel: +41 44 63 55 44 97  /  Mobile: +41 79 71 90 935
>
> --
> 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...@googlegroups.com.
> To post to this group, send email to cp...@googlegroups.com.
> Visit this group at http://groups.google.com/group/cp2k [1].
> For more options, visit https://groups.google.com/d/optout [2].
>
>
> Links:
> ------
> [1] http://groups.google.com/group/cp2k
> [2] https://groups.google.com/d/optout
> [3] http://www.iki.fi/~apsi/

Teodoro Laino

unread,
Mar 24, 2014, 2:36:48 AM3/24/14
to cp...@googlegroups.com
On 24 Mar 2014, at 00:00, Ole Schütt <o...@schuett.name> wrote:


Actually, if you are using the latest svn version you now have to type "make graph", because the names of the make-targets now match the names of the executables. Alternatively, you can just call make without any targets, which will build all executables including the fes tool.


Not entirely true - using the SVN version, one can still compile FES, by typing 

> make fes

alternatively also 

> make graph

works..

Teo


Visit this group at http://groups.google.com/group/cp2k.
For more options, visit https://groups.google.com/d/optout.

T. Salavati fard

unread,
Mar 24, 2014, 7:50:37 PM3/24/14
to cp...@googlegroups.com
Thanks all for your reply.
I've compiled fes and it seems working. But the problem is when I use it, in fes.dat all values are zero.
I've used following synatax:
fes.sopt [-cpmd] in a folder which contains colvar and parvar files.
any ideas?
Thanks in advance--

Teodoro Laino

unread,
Mar 25, 2014, 3:13:39 AM3/25/14
to cp...@googlegroups.com
fes.sopt -help

will give you an help.. Try to go though it.

You’ll have to specify -ndim and also -ndw (if I remember correctly).

Teo

Taha Salavati-fard

unread,
Mar 25, 2014, 11:10:20 AM3/25/14
to cp...@googlegroups.com
Dear Teo,

I've tried that but still zeros!

FES|  Parsing file:   <colvar_mtd>

**********************************************************************
FES|  Parameters for FES:

FES|  NDIM         ::           2
FES|  NWD          ::           2
FES|  HILLS        ::        2156
FES|  COLVAR #   1 ::     (  2.294 , 10.688)               (   MAPPED)
FES|  COLVAR #   2 ::     (  2.278 , 10.693)               (   MAPPED)
FES|  NGRID        ::         100    100
FES|  DX           ::       0.085  0.085
FES|  CUTOFF       ::     .10000E-06
**********************************************************************


FES|  Computing Free Energy Surface
FES|  Mapping Gaussians ::  10 %
FES|  Mapping Gaussians ::  20 %
FES|  Mapping Gaussians ::  30 %
FES|  Mapping Gaussians ::  40 %
FES|  Mapping Gaussians ::  50 %
FES|  Mapping Gaussians ::  60 %
FES|  Mapping Gaussians ::  70 %
FES|  Mapping Gaussians ::  80 %
FES|  Mapping Gaussians ::  90 %
FES|  Mapping Gaussians :: 100 %
FES|  Dumping FES structure in file: < fes.dat >

FES| NORMAL FES TERMINATION.



--
You received this message because you are subscribed to a topic in the Google Groups "cp2k" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/cp2k/NgPWFuuJ3ak/unsubscribe.
To unsubscribe from this group and all its topics, send an email to cp2k+uns...@googlegroups.com.

Teodoro Laino

unread,
Mar 25, 2014, 11:16:50 AM3/25/14
to cp...@googlegroups.com
can you post your files?

Taha Salavati-fard

unread,
Mar 25, 2014, 11:25:35 AM3/25/14
to cp...@googlegroups.com
Yes of course,
Here they are!
colvar_mtd
parvar_mtd

Teodoro Laino

unread,
Mar 26, 2014, 1:09:39 AM3/26/14
to cp...@googlegroups.com
On my laptop, when I run:

fes.sopt -cpmd -ndim 2

I get the fes.dat file attached. Which is everything but not zero.
Can you please attach your fes.dat output file generated by fes.sopt ?

Teo

fes.dat

Taha Salavati-fard

unread,
Mar 26, 2014, 9:21:46 AM3/26/14
to cp...@googlegroups.com
Thanks Teo.
Problem might be due to old linux on the cluster I am using. I am going to try it again on my laptop and will let you know.
Thanks for great helps.



--
You received this message because you are subscribed to a topic in the Google Groups "cp2k" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/cp2k/NgPWFuuJ3ak/unsubscribe.
To unsubscribe from this group and all its topics, send an email to cp2k+uns...@googlegroups.com.
To post to this group, send email to cp...@googlegroups.com.
Visit this group at http://groups.google.com/group/cp2k.
For more options, visit https://groups.google.com/d/optout.

<colvar_mtd><parvar_mtd>



Teodoro Laino

unread,
Mar 26, 2014, 10:53:57 AM3/26/14
to cp...@googlegroups.com
please, post that file.

To unsubscribe from this group and all its topics, send an email tocp2k+uns...@googlegroups.com.

To post to this group, send email to cp...@googlegroups.com.
Visit this group at http://groups.google.com/group/cp2k.
For more options, visit https://groups.google.com/d/optout.

-- 
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 tocp2k+uns...@googlegroups.com.
To post to this group, send email to cp...@googlegroups.com.
Visit this group at http://groups.google.com/group/cp2k.
For more options, visit https://groups.google.com/d/optout.


-- 

You received this message because you are subscribed to a topic in the Google Groups "cp2k" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/cp2k/NgPWFuuJ3ak/unsubscribe.
To unsubscribe from this group and all its topics, send an email to cp2k+uns...@googlegroups.com.
To post to this group, send email to cp...@googlegroups.com.
Visit this group at http://groups.google.com/group/cp2k.
For more options, visit https://groups.google.com/d/optout.

-- 
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 tocp2k+uns...@googlegroups.com.
To post to this group, send email to cp...@googlegroups.com.
Visit this group at http://groups.google.com/group/cp2k.
For more options, visit https://groups.google.com/d/optout.
<colvar_mtd><parvar_mtd>



Taha Salavati-fard

unread,
Mar 26, 2014, 11:15:14 AM3/26/14
to cp...@googlegroups.com
here it is
fes.dat

Teodoro Laino

unread,
Mar 26, 2014, 11:32:48 AM3/26/14
to cp...@googlegroups.com
mmm… strange..  it could really be some issue with the compiler..

Teo


For more options, visit https://groups.google.com/d/optout.
<fes.dat>

Taha Salavati-fard

unread,
Mar 26, 2014, 11:41:13 AM3/26/14
to cp...@googlegroups.com
I'll try with different compiler and let you know.
thanks--

Taha Salavati-fard

unread,
Mar 26, 2014, 1:26:58 PM3/26/14
to cp...@googlegroups.com
Dear Teo,

I've compiled it again on different machine with latest ubuntu.
fes.dat this time contains zeros and NaNs !!
I used exactly same files that I have sent to you!

Teodoro Laino

unread,
Mar 26, 2014, 3:42:48 PM3/26/14
to cp...@googlegroups.com
which compiler ?

Taha Salavati-fard

unread,
Mar 26, 2014, 3:45:25 PM3/26/14
to cp...@googlegroups.com
I guess it is gfortran
Reply all
Reply to author
Forward
0 new messages