[gmx-users] Temperature or Pressure Coupling

2016-01-25 Thread sun
Hello I want to simulate a protein in water and observe its behavior alone and then in the presence of ligand. I have read somewhere that NPT MD is best for pro-lig complexes as it resembles the in vitro and in vivo conditions (conformational changes,biomolecular reactions, binding etc.). So My

Re: [gmx-users] g_correlation

2016-01-25 Thread Felix W.-H. Weng
Hi Dr. Bowerman: Thanks for replying. Happy holidays! Best, Wei-Hsiang -- Wei-Hsiang Weng (翁偉翔) Department of Life Sciences Tzu-Chi University No.701, Sec. 3, Zhongyang Rd. Hualien 97004, Taiwan Tel.: +886-975-232-245 2016-01-26 11:18 GMT+08:00 Samuel Bowerman : > Hello Wei-Hsiang, > > G_c

Re: [gmx-users] g_correlation

2016-01-25 Thread Samuel Bowerman
Hello Wei-Hsiang, G_correlation will print a 0-1 value automatically. If you provide the "-mi" flag, then the output is raw mutual information (range 0 to infinity). You can not "set a range" for the values; they are what they naturally are. You simply choose what formalism to work in (correla

[gmx-users] g_correlation

2016-01-25 Thread Felix W.-H. Weng
Dear all: I'm using the generalized correlation analysis tool developed by Lange and Grubmüller. Does anyone know the option for setting the range of correlation values? Say, 0 to 1? Any advice would be appreciated, thanks in advance. Best, Wei-Hsiang -- Wei-Hsiang Weng (翁偉翔) Department of Life

Re: [gmx-users] GPU efficiency

2016-01-25 Thread Szilárd Páll
Hi Albert, I have to concur with Mark's point, it's not hard at all to pick (or happen to benchmark) circumstances under which one code is significantly faster than the other. In order to understand the differences, please share: - GROMACS logs [ - full hardware & software specs, compilers, CUDA,

[gmx-users] Brefeldin A and GDP parameters

2016-01-25 Thread Jeremy T First
Hello all, Has anyone parameterized Brefeldin A or GDP for Amber force fields? I wish to simulate a protein with BFA and GDP as ligands, but have so far been unable to find parameters for these molecules. I have a procedure to do so using RED and Gamess, but understand this is a difficult process

Re: [gmx-users] Temperature Annealing

2016-01-25 Thread Justin Lemkul
On 1/25/16 4:33 PM, abdus sabuj wrote: Hi Justin, I have used these following commands in npt and final .tpr production. grompp -f npt.mdp -p topol.top -c nvt.gro -o npt.tpr grompp -f final.mdp -p topol.top -c npt.gro -o final.tpr ...which explains everything. You're not preserving the

Re: [gmx-users] Temperature Annealing

2016-01-25 Thread abdus sabuj
Hi Justin, I have used these following commands in npt and final .tpr production. grompp -f npt.mdp -p topol.top -c nvt.gro -o npt.tpr grompp -f final.mdp -p topol.top -c npt.gro -o final.tpr Thanks, Sabuj On Mon, Jan 25, 2016 at 1:55 PM, Justin Lemkul wrote: > > > On 1/25/16 2:41 PM, abdu

Re: [gmx-users] Changing a quantity at each time step using a formula

2016-01-25 Thread Mark Abraham
Hi, It's not clear what you intend to do to preserve neutrality, but you can get charges to change incrementally over time with the slow-growth free-energy perturbation code. Whether that's a useful model remains to be seen... Mark On Mon, Jan 25, 2016 at 9:06 PM AshutoshAkshay Shah wrote: > H

[gmx-users] Changing a quantity at each time step using a formula

2016-01-25 Thread AshutoshAkshay Shah
Hello GROMACS users, I am simulating a graphene based supercapacitor. I want to change the charge of each particle at the electrode surface such that the voltage of the electrochemical cell remains constant. I know the formula to be used for this, but I have no idea how to implement it. I want to

Re: [gmx-users] Temperature Annealing

2016-01-25 Thread Justin Lemkul
On 1/25/16 2:41 PM, abdus sabuj wrote: Hi Justin, Here is the mdp file I used for NPT equilibration... integrator = md nsteps = 1000 dt = 0.001 nstenergy= 1000 nstxout = 1000 nstvout

Re: [gmx-users] Error - syntax

2016-01-25 Thread Justin Lemkul
On 1/25/16 2:50 PM, Chang Woon Jang wrote: Dear Mark, Thank you for the information. I have changed the topol.top file now like I put another three lined for posres right below each itp, DGA.itp and J230.itp. ; Include forcefield parameters ;#include "forcefield.itp" ; Include chain top

Re: [gmx-users] Error - syntax

2016-01-25 Thread Chang Woon Jang
Dear Mark, Thank you for the information. I have changed the topol.top file now like I put another three lined for posres right below each itp, DGA.itp and J230.itp. ; Include forcefield parameters ;#include "forcefield.itp" ; Include chain topologies #include "DGA.itp" #ifdef POSRES #includ

Re: [gmx-users] Temperature Annealing

2016-01-25 Thread abdus sabuj
Hi Justin, Here is the mdp file I used for NPT equilibration... integrator = md nsteps = 1000 dt = 0.001 nstenergy= 1000 nstxout = 1000 nstvout = 1000 vdw-type

Re: [gmx-users] gromacs.org_gmx-users Digest, Vol 141, Issue 106

2016-01-25 Thread lloyd riggs
Dear Sabuj,   It would be difficult to help without more info (size of system, contents of system, .mdp file, etc...).  I have not used simulated annealing in years, and when I used the software in the past, it was mostly structure determination from unknown phases (cooling an unknown partial st

Re: [gmx-users] Build fails with AVX2_256 on Intel Xeon CPU E5-2680 v3

2016-01-25 Thread Koji Tanaka
Hi Nicolas, have you installed the latest binutils before going through the installation? I used gcc/5.3.0, but the versions you've tried should work fine, too. Koji On Mon, Jan 25, 2016 at 12:56 PM, Nicolas Cheron < nicolas.cheron.bou...@gmail.com> wrote: > Hi, > > All my installations are made

Re: [gmx-users] Temperature Annealing

2016-01-25 Thread Justin Lemkul
On 1/25/16 1:51 PM, abdus sabuj wrote: Hi Mark, I have equilibrated my system very well before the production run at 300K, but at time zero, the temperature is at 321K. Here is the fluctuations I am getting. time (ns) temperature (K) 0 321 2 315 4

Re: [gmx-users] Temperature Annealing

2016-01-25 Thread abdus sabuj
Hi Mark, I have equilibrated my system very well before the production run at 300K, but at time zero, the temperature is at 321K. Here is the fluctuations I am getting. time (ns) temperature (K) 0 321 2 315 4 333 6 358

Re: [gmx-users] Temperature Annealing

2016-01-25 Thread Mark Abraham
Hi, What fluctuations were you seeing over the time segments? Mark On Mon, 25 Jan 2016 19:25 abdus sabuj wrote: > Dear gmx users, > > I am trying to get viscosity as a function of temperature. So I have put > the following annealing inputs to the mdp file. > > ;simulated anneling > annealing

Re: [gmx-users] Error - syntax

2016-01-25 Thread Mark Abraham
Hi, See http://www.gromacs.org/Documentation/Errors#Invalid_order_for_directive_xxx Mark On Mon, 25 Jan 2016 19:19 Chang Woon Jang wrote: > Dear Gromacs Users, > > I have two types of molecules. The topology file was separated by one > (topol.top) and two (DGA.itp, and J230.itp). > > The form

Re: [gmx-users] Error in comparison of Coulomb-14 potential?

2016-01-25 Thread Mark Abraham
Hi, That could be evidence of a bug. If you can upload your grompp inputs, tpr and log files to a new issue on redmine.gromacs.org we can consider further. Cheers, Mark On Mon, 25 Jan 2016 15:27 Andreas Mecklenfeld < a.mecklenf...@tu-braunschweig.de> wrote: > Dear all, > > I want to perform an

[gmx-users] Temperature Annealing

2016-01-25 Thread abdus sabuj
Dear gmx users, I am trying to get viscosity as a function of temperature. So I have put the following annealing inputs to the mdp file. ;simulated anneling annealing = single annealing-npoints = 6 annealing-time= 0 2000 4000 6000 8000 1 annealing-temp= 300 310 320 330 340 350

[gmx-users] Error - syntax

2016-01-25 Thread Chang Woon Jang
Dear Gromacs Users, I have two types of molecules. The topology file was separated by one (topol.top) and two (DGA.itp, and J230.itp). The format I made in topol.top is #include "DGA.itp" #include "J230.itp" [ system ] ; Name DGA and J230 [ molecules ] ; Compound#mols Other_chain_A

Re: [gmx-users] Build fails with AVX2_256 on Intel Xeon CPU E5-2680 v3

2016-01-25 Thread Mark Abraham
Hi, Thanks. Indications are that you also need an updated binaries package. Specifically, there's a test that the compiler accepts the avx2 flag, which is passing. But some other part of the infrastructure isn't working (probably the linker) because the next test fails. Mark On Mon, 25 Jan 2016

Re: [gmx-users] Build fails with AVX2_256 on Intel Xeon CPU E5-2680 v3

2016-01-25 Thread Nicolas Cheron
Hi, All my installations are made from a clean build directory. I am using CentOS release 6.6 and the CPU is Intel Xeon CPU E5-2650 v3 @ 2.30GHz. Here is the exact list of command that I have used (as I said, I have also tried with gcc4.9.3 and gcc5.2): module purge module load gcc/4.8.5 cd /home

Re: [gmx-users] Gromos96 53a6 and Berger Lipids

2016-01-25 Thread Piggot T.
Hi, Neither of the two papers you mentioned use the Berger lipid parameters, rather they both use purely GROMOS lipid force fields. This is plain GROMOS 45A3 for the Chandrasekhar et al. paper (using older, original, GROMOS DPPC lipid charges). This force field doesn't reproduce the properties

[gmx-users] Error in comparison of Coulomb-14 potential?

2016-01-25 Thread Andreas Mecklenfeld
Dear all, I want to perform an alchemical change of state of a molecule in TIP3P, where I do only change the partial charges of the atomtypes. To control my results, I perform two MD runs with the different charge sets respectively. I'm interested in the Coulomb-14 interactions of the solute

[gmx-users] Gromos96 53a6 and Berger Lipids

2016-01-25 Thread Kindlein, Moritz
Hi all, I would appreciate your opinion on a quite often discussed topic. I am doing simulations of triglycerides and phosphoglycerides, and I am currently using the combination of gromos96 53a6 and Berger Lipid parameters as suggested in Justin Lemkuls Tutorial. If I am right, this combination

Re: [gmx-users] How to get an index file containing non polar atoms in gromacs?

2016-01-25 Thread Catarina A. Carvalheda dos Santos
I think you can figure this out after playing 5 min with the make_ndx tool, but here it goes: "r ILE | r LEU | r VAL" will select all the isoleucines, leucines and valines in your system. Cheers, On 25 January 2016 at 12:23, Shima ebrahimi wrote: > > > Dear Users, > > I want to determine hydr

[gmx-users] How to get an index file containing non polar atoms in gromacs?

2016-01-25 Thread Shima ebrahimi
Dear Users, I want to determine hydrophobic interaction between ligand and protein using g_mindist. However, how can I identify non-polar atoms and specify them in the index file.Could anyone help me? Regards, Shima -- Gromacs Users mailing list * Please search the archive at http://www.