[gmx-users] interaction energy

2013-04-18 Thread fatemeh ramezani


 HI
I'm simulating gold-protein interaction by gromacs. after MD simulation I want 
to calculate interaction energy of each aminoacids with 

AU surface. 
Enon-bond = E(vdw) + E(elec) 


gold atoms  charge in simulation were considered 0, then E non-bonded=E(vdw)
Is this true?
The interaction energy is equal to the van der Waals energy?


Fatemeh Ramezani
--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] Plotting data from g_spatial

2013-04-18 Thread Venkat Reddy
Dear users

I have a question regarding the cube file output we get from g_spatial
analysis. My system comprises of hundreds of ionic liquids (composed of
cation and anion). If I would like to see the spatial distribution of
anions around cation, how do I get an average distribution. E.g, the way I
am using it now, spatial distribution of all the anions in simulation box
is shown, which is very inconclusive. I cannot specify a particular set of
cations and anions, as they separate over time.

Please help


-- 
With Best Wishes
Venkat Reddy Chirasani
PhD student
Laboratory of Computational Biophysics
Department of Biotechnology
IIT Madras
Chennai
INDIA-600036
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] installing gromacs-4.6.1 with CMake

2013-04-18 Thread 中野佑香
 Hi. I'm trying to build a build with cmake-2.8.10.2-linux for
'gromacs-4.6.1' as the build target.

First,
I extracted gromacs-4.6.1.tar.gz  and, create a new folder; build inside
the extracted folder.
Then, inside build folder, I configured using CMake.

$ tar -xvzf gromacs-4.6.1.tar.gz
$ cd gromacs-4.6.1
$ mkdir build
$ cd build
$ cmake ..

However, I got an error.
--
CMake Error: The source directory /home/okayuka/Gromacs does not appear
to contain CMakeLists.txt.
Specify --help for usage, or press the help button on the CMake GUI.
-
Inside gromacs-4.6.1 directory, there is certainly CMakeLists.txt.

For this, I think cmake mistaking  /home/okayuka/Gromacs for
/home/okayuka/Gromacs/gromacs-4.6.1 as ' the source directory'.

Would you please tell me how correct this mismatch of the source directory?

-Yuka

-- 
*
中野佑香
立教大学 理学研究科化学専攻
博士前期課程 1年
理論創薬・分子設計 常盤研究室
〒171-8501 東京都豊島区西池袋3-34-1
4号館4230教室
Email: yuk...@rikkyo.ac.jp 09cc0...@rikkyo.ac.jp
*
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] Regarding top2psf script

2013-04-18 Thread Venkat Reddy
Dear all,
I want to reverse transform my CG lipid system to AA. I have been using
g_fg2cg but I am getting some errors (like Atoms in the .top are not
numbered consecutively from 1). So, I intend to use the script
top2psf.plwritten by Dr.Lemkul, so that I can avail the help of VMD in
cg2aa
conversion. Can we use it for pure lipid systems (say
POPC+CHOL+triglyceride+CE)?

When I tried to execute it, I am getting the following error
 Cannot open angles for reading: No such file or directory

Thanks for your concern

-- 
With Best Wishes
Venkat Reddy Chirasani
PhD student
Laboratory of Computational Biophysics
Department of Biotechnology
IIT Madras
Chennai
INDIA-600036
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] interaction energy

2013-04-18 Thread Justin Lemkul



On 4/18/13 2:00 AM, fatemeh ramezani wrote:



  HI
I'm simulating gold-protein interaction by gromacs. after MD simulation I want 
to calculate interaction energy of each aminoacids with

AU surface.
Enon-bond = E(vdw) + E(elec)


gold atoms  charge in simulation were considered 0, then E non-bonded=E(vdw)
Is this true?
The interaction energy is equal to the van der Waals energy?



If charges are zero, then yes.

-Justin

--


Justin A. Lemkul, Ph.D.
Research Scientist
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin


--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] installing gromacs-4.6.1 with CMake

2013-04-18 Thread Justin Lemkul



On 4/18/13 3:09 AM, 中野佑香 wrote:

  Hi. I'm trying to build a build with cmake-2.8.10.2-linux for
'gromacs-4.6.1' as the build target.

First,
I extracted gromacs-4.6.1.tar.gz  and, create a new folder; build inside
the extracted folder.
Then, inside build folder, I configured using CMake.

$ tar -xvzf gromacs-4.6.1.tar.gz
$ cd gromacs-4.6.1
$ mkdir build
$ cd build
$ cmake ..

However, I got an error.
--
CMake Error: The source directory /home/okayuka/Gromacs does not appear
to contain CMakeLists.txt.
Specify --help for usage, or press the help button on the CMake GUI.
-
Inside gromacs-4.6.1 directory, there is certainly CMakeLists.txt.

For this, I think cmake mistaking  /home/okayuka/Gromacs for
/home/okayuka/Gromacs/gromacs-4.6.1 as ' the source directory'.

Would you please tell me how correct this mismatch of the source directory?



The sequence of commands you have entered above is correct, but that's 
apparently not what's actually happening in your terminal.  The error indicates 
you are trying to use CMake to build within /home/okayuka/Gromacs, not 
/home/okayuka/Gromacs/gromacs-4.6.1.  That would happen if perhaps you had 
forgotten to cd into gromacs-4.6.1.


-Justin

--


Justin A. Lemkul, Ph.D.
Research Scientist
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin


--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] Regarding top2psf script

2013-04-18 Thread Justin Lemkul



On 4/18/13 5:15 AM, Venkat Reddy wrote:

Dear all,
I want to reverse transform my CG lipid system to AA. I have been using
g_fg2cg but I am getting some errors (like Atoms in the .top are not
numbered consecutively from 1). So, I intend to use the script
top2psf.plwritten by Dr.Lemkul, so that I can avail the help of VMD in
cg2aa
conversion. Can we use it for pure lipid systems (say
POPC+CHOL+triglyceride+CE)?

When I tried to execute it, I am getting the following error
  Cannot open angles for reading: No such file or directory



My script only works on a single molecule topology.  It will not work for a 
complex system like yours.  Perhaps you can process each .itp file separately, 
which may be of some use, but that will only potentially give you an AA 
representation of a single molecule.


-Justin

--


Justin A. Lemkul, Ph.D.
Research Scientist
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin


--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] Applying a Lennard-Jones Gauss potential function for TIP5P water

2013-04-18 Thread Andrish Reddy
Greetings, 

I am simulating a box of TIP5P waters and I would like to know if it is 
possible to change the standard LJ potential function to a more complex one 
(LJ-Gauss), while keeping the coulombic interaction standard. The custom 
potential will act between the oxygen atoms only and has the form: 
eps*[(sig/r)^12 - 2*(sig/r)^6] - epsG*[exp-(((r-rG)^2) / (2*sigG^2))] where 
epsG, sigG and rG are new parameters. 

Can this be accomplished using tables, a custom force field / topology, or 
other? 
I would appreciate any suggestions relating to how this can be done.

Thanks,
Andrish


-- 
This message is subject to the CSIR's copyright terms and conditions, e-mail 
legal notice, and implemented Open Document Format (ODF) standard. 
The full disclaimer details can be found at 
http://www.csir.co.za/disclaimer.html.

This message has been scanned for viruses and dangerous content by MailScanner, 
and is believed to be clean.

Please consider the environment before printing this email.

-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Re: [gmx-users] Applying a Lennard-Jones Gauss potential function for TIP5P water

2013-04-18 Thread Justin Lemkul



On 4/18/13 7:31 AM, Andrish Reddy wrote:

Greetings,

I am simulating a box of TIP5P waters and I would like to know if it is 
possible to change the standard LJ potential function to a more complex one 
(LJ-Gauss), while keeping the coulombic interaction standard. The custom 
potential will act between the oxygen atoms only and has the form:
eps*[(sig/r)^12 - 2*(sig/r)^6] - epsG*[exp-(((r-rG)^2) / (2*sigG^2))] where 
epsG, sigG and rG are new parameters.

Can this be accomplished using tables, a custom force field / topology, or 
other?
I would appreciate any suggestions relating to how this can be done.



Tabulated potentials are the best bet.

-Justin

--


Justin A. Lemkul, Ph.D.
Research Scientist
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin


--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] installing gromacs-4.6.1 with CMake

2013-04-18 Thread 中野佑香
Dear Justin,

Thank you for your fast reply.

cmake was successfully worked after I removed CMakeCache.txt from
gromacs-**4.6.1
directory and
 corrected CMake path in .tcshrc .

Thank you!

-Yuka


2013/4/18 Justin Lemkul jalem...@vt.edu



 On 4/18/13 3:09 AM, 中野佑香 wrote:

   Hi. I'm trying to build a build with cmake-2.8.10.2-linux for
 'gromacs-4.6.1' as the build target.

 First,
 I extracted gromacs-4.6.1.tar.gz  and, create a new folder; build inside
 the extracted folder.
 Then, inside build folder, I configured using CMake.

 $ tar -xvzf gromacs-4.6.1.tar.gz
 $ cd gromacs-4.6.1
 $ mkdir build
 $ cd build
 $ cmake ..

 However, I got an error.
 --**--**
 --**
 CMake Error: The source directory /home/okayuka/Gromacs does not appear
 to contain CMakeLists.txt.
 Specify --help for usage, or press the help button on the CMake GUI.
 --**--**
 --**---
 Inside gromacs-4.6.1 directory, there is certainly CMakeLists.txt.

 For this, I think cmake mistaking  /home/okayuka/Gromacs for
 /home/okayuka/Gromacs/gromacs-**4.6.1 as ' the source directory'.

 Would you please tell me how correct this mismatch of the source
 directory?


 The sequence of commands you have entered above is correct, but that's
 apparently not what's actually happening in your terminal.  The error
 indicates you are trying to use CMake to build within
 /home/okayuka/Gromacs, not /home/okayuka/Gromacs/gromacs-**4.6.1.  That
 would happen if perhaps you had forgotten to cd into gromacs-4.6.1.

 -Justin

 --
 ==**==

 Justin A. Lemkul, Ph.D.
 Research Scientist
 Department of Biochemistry
 Virginia Tech
 Blacksburg, VA
 jalemkul[at]vt.edu | (540) 231-9080
 http://www.bevanlab.biochem.**vt.edu/Pages/Personal/justinhttp://www.bevanlab.biochem.vt.edu/Pages/Personal/justin

 ==**==
 --
 gmx-users mailing listgmx-users@gromacs.org
 http://lists.gromacs.org/**mailman/listinfo/gmx-usershttp://lists.gromacs.org/mailman/listinfo/gmx-users
 * Please search the archive at http://www.gromacs.org/**
 Support/Mailing_Lists/Searchhttp://www.gromacs.org/Support/Mailing_Lists/Searchbefore
  posting!
 * Please don't post (un)subscribe requests to the list. Use the www
 interface or send it to gmx-users-requ...@gromacs.org.
 * Can't post? Read 
 http://www.gromacs.org/**Support/Mailing_Listshttp://www.gromacs.org/Support/Mailing_Lists




-- 
*
中野佑香
立教大学 理学研究科化学専攻
博士前期課程 1年
理論創薬・分子設計 常盤研究室
〒171-8501 東京都豊島区西池袋3-34-1
4号館4230教室
Email: yuk...@rikkyo.ac.jp 09cc0...@rikkyo.ac.jp
*
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] Re: Applying a Lennard-Jones Gauss potential function for TIP5P water

2013-04-18 Thread Andrish Reddy
Thanks Justin, 

I had a look at the manual regarding tabulated potentials... for clarity, 
The table values must contain: x, f(x), −f(x), g(x), −g(x), h(x), −h(x). 

For my function then would f(x) = 1/r   g(x) = eps*[(sig/r)^12 -
2*(sig/r)^6] 
and h(x) = -epsG*[exp-(((r-rG)^2) / (2*sigG^2))] ?

Also from the manual:
The parameters A, B, and C are read from the topology.  The type of
interaction, Lennard-Jones or Buckingham, is determined by the nbfunc
parameter in the [ defaults ] directive of the topology.

How would the A,B  C terms be applied to my function with 5 constants? What
would 'nbfunc' have to be set to?

Regards,
Andrish






--
View this message in context: 
http://gromacs.5086.x6.nabble.com/Applying-a-Lennard-Jones-Gauss-potential-function-for-TIP5P-water-tp5007415p5007418.html
Sent from the GROMACS Users Forum mailing list archive at Nabble.com.
--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] Running simulation over openmpi in rsh

2013-04-18 Thread 라지브간디
Dear Gmx,


Is there anyway I can run simulation using mpirun command via rsh? I am using 
openmpi 1.6 version. 
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

[gmx-users] How to fix interfacial tension in NPT

2013-04-18 Thread Souilem Safa
Dear Gromacs users,
I have runned a semiisotropic NPT simulation for a membrane, I have fixed
the pressure in x and y axis and I varied it only in z axis to keep a
stable interface. But after simulation the obtained interfacial tension is
far bigger than the experimental value. Experimental is 27 mN.m-1 and the
simulated one was around 38 mN.m-1.
I will be very grateful if you can suggest me which parameters should I
check in order to regulate the interfacial tension value.
Thank you very much
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] Membrane Simulation

2013-04-18 Thread Giuseppe


Dear all
 
I'm trying to start a simulation with a protein embedded in a bilayer of
POPC molecules, but I have some problems when I pack the lipids around the
protein. Infact I place the protein in the proper position using the
editconf tool with the “-rotate” and “-center” options on, but when I scale
the lipid positions with the perl script (before the energy minimization),
the protein moves up, de facto off from the membrane! Obviously I've set
very strong position-restraining force on protein, but it moves the same.
Probably the problem is trivial, but I can't see it, so any help is much
appreciated.
 
Thanks a lot,
 
GIuseppe Cimicata



--
View this message in context: 
http://gromacs.5086.x6.nabble.com/Membrane-Simulation-tp5007421.html
Sent from the GROMACS Users Forum mailing list archive at Nabble.com.
--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] ffnonbonded value of charmm to charmm27 ff

2013-04-18 Thread 라지브간디
Dear gmx,


I m having a confusion of using the literature value of non-bonded ( charmm 
format) to use in gromacs over charmm27.ff ? 


For example, the literature reported the Carbon nonbonded as ( sigma =3.83, 
epsilon -0.0262 ), whereas the charmm27.ff in gromacs- ffnonbonded value listed 
for Carbon listed as (sigma 0.356359487256, epsilon  0.46024). 


so whats is the procedure to convert this literature value favor to charmm27 in 
gromacs? 


Thanks in advance.




Regards


Rajiv




-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

RE: [gmx-users] Building Single and Double Precision in 4.6.1?

2013-04-18 Thread Mike Hanby
Thanks for the reply, so the next question, after I finish building single 
precision non parallel, is there an efficient way to kick off the double 
precision build, then the single precision mpi and so on?

Or do I need to delete everything from my build directory before running cmake 
again?

One concern is that each time I run cmake, it will download fftw and the 
regression tests.


-Original Message-
From: gmx-users-boun...@gromacs.org [mailto:gmx-users-boun...@gromacs.org] On 
Behalf Of Mark Abraham
Sent: Monday, April 15, 2013 11:05 AM
To: Discussion list for GROMACS users
Subject: Re: [gmx-users] Building Single and Double Precision in 4.6.1?

That's best. GROMACS will take care of a _d suffix for double precision for you.

Mark


On Mon, Apr 15, 2013 at 5:39 PM, Mike Hanby mha...@uab.edu wrote:

 Howdy,

 What's the proper procedure for building both single and double 
 precision and installing them to the same directory?

 For example, is this the proper set of steps?

 #Build and install single precision
 cmake .. \
-DGMX_BUILD_OWN_FFTW=ON \
-DCMAKE_INSTALL_PREFIX=/share/apps/gromacs/4.6.1 \
-DREGRESSIONTEST_DOWNLOAD=ON
 make
 make check
 make install

 cmake .. \
-DGMX_BUILD_OWN_FFTW=ON \
-DGMX_DOUBLE=ON \
-DCMAKE_INSTALL_PREFIX=/share/apps/gromacs/4.6.1 \
-DREGRESSIONTEST_DOWNLOAD=ON
 make
 make check
 make install

 Or is there a more efficient way? Thanks, Mike
 --
 gmx-users mailing listgmx-users@gromacs.org
 http://lists.gromacs.org/mailman/listinfo/gmx-users
 * Please search the archive at
 http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
 * Please don't post (un)subscribe requests to the list. Use the www 
 interface or send it to gmx-users-requ...@gromacs.org.
 * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the www interface 
or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] Building Single and Double Precision in 4.6.1?

2013-04-18 Thread Szilárd Páll
On Thu, Apr 18, 2013 at 6:17 PM, Mike Hanby mha...@uab.edu wrote:
 Thanks for the reply, so the next question, after I finish building single 
 precision non parallel, is there an efficient way to kick off the double 
 precision build, then the single precision mpi and so on?

 Or do I need to delete everything from my build directory before running 
 cmake again?

Well, most setting can just be changed in the build tree and rerunning
cmake build trigger a rebuild, e.g.:
$ cmake path/to/source -DGMX_DOUBLE=OFF
-DCMAKE_INSTALL_PREFIX=/my/install/path; make install -j; cmake
.-DGMX_DOUBLE=ON; make install -j;

However, I recommend that you do your builds outside the source tree
in separate build directories - or at least in separate directories in
the source tree - and just use the same CMAKE_INSTALL_PREFIX for every
build.


 One concern is that each time I run cmake, it will download fftw and the 
 regression tests.

If you want to do regression tests of multiple builds, to avoid the
downloading, you should download the regressiontests tarball manually
and just point cmake to it instead of having it download every time.

While you'll need a different FFTW build for double precision, you are
right that single precision MPI vs thread-MPI builds would not need a
different FFTW. If you want to avoid this (IMHO minor) overhead,
compile FFTW manually and just point CMake to it (via
CMAKE_PREFIX_PATH).

Cheers,
--
Szilárd




 -Original Message-
 From: gmx-users-boun...@gromacs.org [mailto:gmx-users-boun...@gromacs.org] On 
 Behalf Of Mark Abraham
 Sent: Monday, April 15, 2013 11:05 AM
 To: Discussion list for GROMACS users
 Subject: Re: [gmx-users] Building Single and Double Precision in 4.6.1?

 That's best. GROMACS will take care of a _d suffix for double precision for 
 you.

 Mark


 On Mon, Apr 15, 2013 at 5:39 PM, Mike Hanby mha...@uab.edu wrote:

 Howdy,

 What's the proper procedure for building both single and double
 precision and installing them to the same directory?

 For example, is this the proper set of steps?

 #Build and install single precision
 cmake .. \
-DGMX_BUILD_OWN_FFTW=ON \
-DCMAKE_INSTALL_PREFIX=/share/apps/gromacs/4.6.1 \
-DREGRESSIONTEST_DOWNLOAD=ON
 make
 make check
 make install

 cmake .. \
-DGMX_BUILD_OWN_FFTW=ON \
-DGMX_DOUBLE=ON \
-DCMAKE_INSTALL_PREFIX=/share/apps/gromacs/4.6.1 \
-DREGRESSIONTEST_DOWNLOAD=ON
 make
 make check
 make install

 Or is there a more efficient way? Thanks, Mike
 --
 gmx-users mailing listgmx-users@gromacs.org
 http://lists.gromacs.org/mailman/listinfo/gmx-users
 * Please search the archive at
 http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
 * Please don't post (un)subscribe requests to the list. Use the www
 interface or send it to gmx-users-requ...@gromacs.org.
 * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

 --
 gmx-users mailing listgmx-users@gromacs.org
 http://lists.gromacs.org/mailman/listinfo/gmx-users
 * Please search the archive at 
 http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
 * Please don't post (un)subscribe requests to the list. Use the www interface 
 or send it to gmx-users-requ...@gromacs.org.
 * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
 --
 gmx-users mailing listgmx-users@gromacs.org
 http://lists.gromacs.org/mailman/listinfo/gmx-users
 * Please search the archive at 
 http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
 * Please don't post (un)subscribe requests to the list. Use the
 www interface or send it to gmx-users-requ...@gromacs.org.
 * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] Re: Applying a Lennard-Jones Gauss potential function for TIP5P water

2013-04-18 Thread Justin Lemkul



On 4/18/13 8:44 AM, Andrish Reddy wrote:

Thanks Justin,

I had a look at the manual regarding tabulated potentials... for clarity,
The table values must contain: x, f(x), −f(x), g(x), −g(x), h(x), −h(x).

For my function then would f(x) = 1/r   g(x) = eps*[(sig/r)^12 -
2*(sig/r)^6]
and h(x) = -epsG*[exp-(((r-rG)^2) / (2*sigG^2))] ?

Also from the manual:
The parameters A, B, and C are read from the topology.  The type of
interaction, Lennard-Jones or Buckingham, is determined by the nbfunc
parameter in the [ defaults ] directive of the topology.

How would the A,B  C terms be applied to my function with 5 constants? What
would 'nbfunc' have to be set to?



Sorry, I have no experience with complex cases like this.  Set up a small test 
case that can be easily verified and try your luck.  Tabulated potentials are 
definitely the way to go though, unless you want to completely overhaul the 
Gromacs code to implement the new potential equation expression :)


-Justin

--


Justin A. Lemkul, Ph.D.
Research Scientist
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin


--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] How to fix interfacial tension in NPT

2013-04-18 Thread Justin Lemkul



On 4/18/13 10:51 AM, Souilem Safa wrote:

Dear Gromacs users,
I have runned a semiisotropic NPT simulation for a membrane, I have fixed
the pressure in x and y axis and I varied it only in z axis to keep a
stable interface. But after simulation the obtained interfacial tension is
far bigger than the experimental value. Experimental is 27 mN.m-1 and the
simulated one was around 38 mN.m-1.
I will be very grateful if you can suggest me which parameters should I
check in order to regulate the interfacial tension value.


You haven't posted your .mdp file or even stated what force field you're using. 
 Is your force field designed to reproduce such surface tension values?  The 
results are only as good as the model itself, and your usage of it.


-Justin

--


Justin A. Lemkul, Ph.D.
Research Scientist
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin


--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] Membrane Simulation

2013-04-18 Thread Justin Lemkul



On 4/18/13 10:59 AM, Giuseppe wrote:



Dear all

I'm trying to start a simulation with a protein embedded in a bilayer of
POPC molecules, but I have some problems when I pack the lipids around the
protein. Infact I place the protein in the proper position using the
editconf tool with the “-rotate” and “-center” options on, but when I scale
the lipid positions with the perl script (before the energy minimization),
the protein moves up, de facto off from the membrane! Obviously I've set
very strong position-restraining force on protein, but it moves the same.
Probably the problem is trivial, but I can't see it, so any help is much
appreciated.



You probably haven't set up the protein's location in a box with the same 
dimensions as the membrane, so it's centered within some different unit cell.


-Justin

--


Justin A. Lemkul, Ph.D.
Research Scientist
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin


--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] Free energy on GPU?

2013-04-18 Thread Олег Титов
Dear all.

I'm trying to estimate relative solvation enegry of brombenzene (relative
to benzene) via thermodynamic integration. Can I perform such calculation
in GROMACS with GPU acceleration?

As I currently understand to make use of GPU one need to use Verlet cutoff
scheme. When I do so with pme electrostatics everything works fine until
free energy calculation stage. Something strange happens there: the
temperature of the system increases with increase of lambda value with no
respect to thermostat (300 at lambda = 0.0 - 500 at lambda = 1.0).
Interesting fact that with reaction-field electrostatics everything seemed
to work fine but now I cannot be sure that these results are reliable.

As I understood from manual there is some issue with Verlet cutoff scheme
and free energy calculations. Can I perform TI calculation in GROMACS with
GPU acceleration or should I use group cutoffs and forget about GPU?

Thanks a lot,
Oleg Titov
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] Free energy on GPU?

2013-04-18 Thread Justin Lemkul



On 4/18/13 2:30 PM, Олег Титов wrote:

Dear all.

I'm trying to estimate relative solvation enegry of brombenzene (relative
to benzene) via thermodynamic integration. Can I perform such calculation
in GROMACS with GPU acceleration?

As I currently understand to make use of GPU one need to use Verlet cutoff
scheme. When I do so with pme electrostatics everything works fine until
free energy calculation stage. Something strange happens there: the
temperature of the system increases with increase of lambda value with no
respect to thermostat (300 at lambda = 0.0 - 500 at lambda = 1.0).
Interesting fact that with reaction-field electrostatics everything seemed
to work fine but now I cannot be sure that these results are reliable.

As I understood from manual there is some issue with Verlet cutoff scheme
and free energy calculations. Can I perform TI calculation in GROMACS with
GPU acceleration or should I use group cutoffs and forget about GPU?



As far as I am aware, free energy calculations are not supported on GPU.

-Justin

--


Justin A. Lemkul, Ph.D.
Research Scientist
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin


--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] ffnonbonded value of charmm to charmm27 ff

2013-04-18 Thread Mark Abraham
Identify the equations into which they are inserted. Apply scaling
parameters to change the units. Learn to check others' work before
investing months of effort ;-)

Mark


On Thu, Apr 18, 2013 at 5:39 PM, 라지브간디 ra...@kaist.ac.kr wrote:

 Dear gmx,


 I m having a confusion of using the literature value of non-bonded (
 charmm format) to use in gromacs over charmm27.ff ?


 For example, the literature reported the Carbon nonbonded as ( sigma
 =3.83, epsilon -0.0262 ), whereas the charmm27.ff in gromacs- ffnonbonded
 value listed for Carbon listed as (sigma 0.356359487256, epsilon  0.46024).


 so whats is the procedure to convert this literature value favor to
 charmm27 in gromacs?


 Thanks in advance.




 Regards


 Rajiv





 --
 gmx-users mailing listgmx-users@gromacs.org
 http://lists.gromacs.org/mailman/listinfo/gmx-users
 * Please search the archive at
 http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
 * Please don't post (un)subscribe requests to the list. Use the
 www interface or send it to gmx-users-requ...@gromacs.org.
 * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] Autocorrelation of C-N

2013-04-18 Thread Mark Abraham
Chapter 8 is your friend. Find a tool to feed data to g_analyze.

Mark


On Wed, Apr 17, 2013 at 4:23 PM, Steven Neumann s.neuman...@gmail.comwrote:

 Dear Users,

 Could you advise me please how to calculate vector C-N autocorrelation
 function  in my protein along the simulation time?

 Steven
 --
 gmx-users mailing listgmx-users@gromacs.org
 http://lists.gromacs.org/mailman/listinfo/gmx-users
 * Please search the archive at
 http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
 * Please don't post (un)subscribe requests to the list. Use the
 www interface or send it to gmx-users-requ...@gromacs.org.
 * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] Running simulation over openmpi in rsh

2013-04-18 Thread Mark Abraham
No idea. Configuring your MPI system is not really on topic for the GROMACS
mailing list. I suggest you read your OpenMPI documentation.

Mark


On Thu, Apr 18, 2013 at 2:58 PM, 라지브간디 ra...@kaist.ac.kr wrote:

 Dear Gmx,


 Is there anyway I can run simulation using mpirun command via rsh? I am
 using openmpi 1.6 version.

 --
 gmx-users mailing listgmx-users@gromacs.org
 http://lists.gromacs.org/mailman/listinfo/gmx-users
 * Please search the archive at
 http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
 * Please don't post (un)subscribe requests to the list. Use the
 www interface or send it to gmx-users-requ...@gromacs.org.
 * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] pdb2gmx for small molecules?

2013-04-18 Thread Warren Gallin
Hi,

I'm not clear on how to use pdb2gmx to set up a small molecule for 
simulation, or perhaps the limitations on this capability.  My understanding 
from the latest documentation and the 2013 publication is that v4.5 is capable 
of handling non-protein/nucleic acid as input.

I have created a pdb file for my molecule of interest, but when I 
invoke pdb2gmx as follows:

pdb2gmx -f BrMT_mono_nores.pdb -water tip4

I get a Fatal error message because there is no residue name (it's 
neither protein or nucleic acid):

---
Program pdb2gmx, VERSION 4.5.4
Source code file: /Users/wgallin/Desktop/gromacs-4.5.4/src/kernel/resall.c, 
line: 581

Fatal error:
Residue '' not found in residue topology database
For more information and tips for troubleshooting, please check the GROMACS
website at http://www.gromacs.org/Documentation/Errors
---



 and just above that message there is also a Warning in the message saying:

Warning: Starting residue 1 in chain not identified as Protein/RNA/DNA.
Problem with chain definition, or missing terminal residues.
This chain does not appear to contain a recognized chain molecule.
If this is incorrect, you can edit residuetypes.dat to modify the behavior.
8 out of 8 lines of specbond.dat converted successfully


Is there some documentation that specifically addresses the methods and issues 
with using pdb2gmx for non-protein/nucleic acid molecules?

Thanks.

Warren Gallin--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] pdb2gmx for small molecules?

2013-04-18 Thread Justin Lemkul



On 4/18/13 3:32 PM, Warren Gallin wrote:

Hi,

I'm not clear on how to use pdb2gmx to set up a small molecule for 
simulation, or perhaps the limitations on this capability.  My understanding 
from the latest documentation and the 2013 publication is that v4.5 is capable 
of handling non-protein/nucleic acid as input.

I have created a pdb file for my molecule of interest, but when I 
invoke pdb2gmx as follows:

pdb2gmx -f BrMT_mono_nores.pdb -water tip4

I get a Fatal error message because there is no residue name (it's 
neither protein or nucleic acid):

---
Program pdb2gmx, VERSION 4.5.4
Source code file: /Users/wgallin/Desktop/gromacs-4.5.4/src/kernel/resall.c, 
line: 581

Fatal error:
Residue '' not found in residue topology database
For more information and tips for troubleshooting, please check the GROMACS
website at http://www.gromacs.org/Documentation/Errors
---



  and just above that message there is also a Warning in the message saying:

Warning: Starting residue 1 in chain not identified as Protein/RNA/DNA.
Problem with chain definition, or missing terminal residues.
This chain does not appear to contain a recognized chain molecule.
If this is incorrect, you can edit residuetypes.dat to modify the behavior.
8 out of 8 lines of specbond.dat converted successfully


Is there some documentation that specifically addresses the methods and issues 
with using pdb2gmx for non-protein/nucleic acid molecules?



pdb2gmx only knows how to handle what it knows about.  There is no generic 
support for arbitrary compounds.


http://www.gromacs.org/Documentation/How-tos/Adding_a_Residue_to_a_Force_Field

-Justin

--


Justin A. Lemkul, Ph.D.
Research Scientist
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin


--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] pdb2gmx for small molecules?

2013-04-18 Thread David van der Spoel

On 2013-04-18 21:35, Justin Lemkul wrote:



On 4/18/13 3:32 PM, Warren Gallin wrote:

Hi,

I'm not clear on how to use pdb2gmx to set up a small molecule for
simulation, or perhaps the limitations on this capability.  My
understanding from the latest documentation and the 2013 publication
is that v4.5 is capable of handling non-protein/nucleic acid as input.

I have created a pdb file for my molecule of interest, but when I
invoke pdb2gmx as follows:

pdb2gmx -f BrMT_mono_nores.pdb -water tip4

I get a Fatal error message because there is no residue name (it's
neither protein or nucleic acid):

---
Program pdb2gmx, VERSION 4.5.4
Source code file:
/Users/wgallin/Desktop/gromacs-4.5.4/src/kernel/resall.c, line: 581

Fatal error:
Residue '' not found in residue topology database
For more information and tips for troubleshooting, please check the
GROMACS
website at http://www.gromacs.org/Documentation/Errors
---



  and just above that message there is also a Warning in the message
saying:

Warning: Starting residue 1 in chain not identified as Protein/RNA/DNA.
Problem with chain definition, or missing terminal residues.
This chain does not appear to contain a recognized chain molecule.
If this is incorrect, you can edit residuetypes.dat to modify the
behavior.
8 out of 8 lines of specbond.dat converted successfully


Is there some documentation that specifically addresses the methods
and issues with using pdb2gmx for non-protein/nucleic acid molecules?



pdb2gmx only knows how to handle what it knows about.  There is no
generic support for arbitrary compounds.

http://www.gromacs.org/Documentation/How-tos/Adding_a_Residue_to_a_Force_Field


-Justin


Try GAFF using the Ambertools and acpype.

--
David van der Spoel, Ph.D., Professor of Biology
Dept. of Cell  Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205.
sp...@xray.bmc.uu.sehttp://folding.bmc.uu.se
--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] Announcement: new version of GridMAT-MD available

2013-04-18 Thread Justin Lemkul


Greetings all,

I wanted to announce that we have released version 2.0 of GridMAT-MD.  Many of 
you use this program for membrane analysis and I am pleased to note that we have 
introduced many new features based on your feedback, including trajectory 
support (multi-frame .pdb and .gro files).  Please visit our webpage for 
details: http://www.bevanlab.biochem.vt.edu/GridMAT-MD/


We hope you continue to find our software useful.  Please contact me if any 
issues arise.


-Justin

--


Justin A. Lemkul, Ph.D.
Research Scientist
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin


--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] A beginner's question ...

2013-04-18 Thread Venkatesh Ramakrishnan
Dear All

I am new to GROMACS, I have the following question:

What is the difference between Pres. DC and Pressure in g_energy analysis ?

Thanks in advance

Venkatesh
--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] A beginner's question ...

2013-04-18 Thread Justin Lemkul



On 4/18/13 7:22 PM, Venkatesh Ramakrishnan wrote:

Dear All

I am new to GROMACS, I have the following question:

What is the difference between Pres. DC and Pressure in g_energy analysis ?



Pres. DC is the correction applied to the pressure term in dispersion 
correction when using DispCorr = EnerPres.


-Justin

--


Justin A. Lemkul, Ph.D.
Research Scientist
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin


--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


[gmx-users] Different energy using the identical tpr file

2013-04-18 Thread Fahimeh Baftizadeh
Hello,

I have a 2 tpr files which are identical. Then I am doing a single step MD
(I put n_step=0) simulation just to compute the total energy of a certain
configuration.

These tpr files give me different energy values ... when I run it on my
computer or on another computer. Even the gromacs version is the same in
both computers. I compiled gromacs identical in both computers.

Can you help me to figure it out?
Thanks
Fahimeh
-- 
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists


Re: [gmx-users] Announcement: new version of GridMAT-MD available

2013-04-18 Thread Albert
It seems to be good. However, I always have problem to compile matpack 
which is needed by xmatrix... I sent the author an email, but never 
replied..


here is the log file:


creating libpng12.la
(cd .libs  rm -f libpng12.la  ln -s ../libpng12.la libpng12.la)
make[4]: Leaving directory 
`/home/albert/Desktop/matpack/source/3rdparty/libpng'
make[3]: Leaving directory 
`/home/albert/Desktop/matpack/source/3rdparty/libpng'

(cd Matutil ; make )
make[3]: Entering directory `/home/albert/Desktop/matpack/source/Matutil'
g++ -ansi -std=c++0x -m64 -DXPM_INCLUDE=X11/xpm.h -c -s -Wall -O3 
-fforce-addr -funroll-loops -felide-constructors -I ../../include 
mpromannumber.cpp
g++ -ansi -std=c++0x -m64 -DXPM_INCLUDE=X11/xpm.h -c -s -Wall -O3 
-fforce-addr -funroll-loops -felide-constructors -I ../../include 
mpcontextsave.cpp
g++ -ansi -std=c++0x -m64 -DXPM_INCLUDE=X11/xpm.h -c -s -Wall -O3 
-fforce-addr -funroll-loops -felide-constructors -I ../../include 
mpparsetool.cpp
g++ -ansi -std=c++0x -m64 -DXPM_INCLUDE=X11/xpm.h -c -s -Wall -O3 
-fforce-addr -funroll-loops -felide-constructors -I ../../include 
mpgetopt.cpp
mpgetopt.cpp: In member function ‘bool 
MATPACK::MpGetopt::assign(MATPACK::MpGetopt::OptNode*, const char**, 
std::string)’:

mpgetopt.cpp:314: error: ‘strtol’ was not declared in this scope
mpgetopt.cpp:330: error: ‘strtoul’ was not declared in this scope
mpgetopt.cpp:345: error: ‘strtod’ was not declared in this scope
g++ -ansi -std=c++0x -m64 -DXPM_INCLUDE=X11/xpm.h -c -s -Wall -O3 
-fforce-addr -funroll-loops -felide-constructors -I ../../include 
mptimerqueue.cpp

make[3]: *** [mpgetopt.o] Error 1
make[3]: *** Waiting for unfinished jobs
make[3]: Leaving directory `/home/albert/Desktop/matpack/source/Matutil'
make[2]: *** [static_lib] Error 2
make[2]: Leaving directory `/home/albert/Desktop/matpack/source'
make[1]: *** [lib] Error 2
make[1]: Leaving directory `/home/albert/Desktop/matpack'
make: *** [allstatic] Error 2


thank you very much

Albert


On 04/18/2013 11:19 PM, Justin Lemkul wrote:


Greetings all,

I wanted to announce that we have released version 2.0 of GridMAT-MD.  
Many of you use this program for membrane analysis and I am pleased to 
note that we have introduced many new features based on your feedback, 
including trajectory support (multi-frame .pdb and .gro files).  
Please visit our webpage for details: 
http://www.bevanlab.biochem.vt.edu/GridMAT-MD/


We hope you continue to find our software useful.  Please contact me 
if any issues arise.


-Justin



--
gmx-users mailing listgmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-requ...@gromacs.org.

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists