[gmx-users] LJ - SR and Coulbom - SR mdp options

2013-10-14 Thread Steven Neumann
Dear Gmx Users, I want to obtain average Protein-Ligand SR electrostatics and LJ energy values using g_energy from NPT equlibration with position restraints. I used energygrps = Protein LIG However, when I run g_energy -f npt.edr -s npt.tpr I see neither Protein-LIG for any of them. Would you pl

Re: [gmx-users] carbon nanotube - ifnite in length ?

2013-10-10 Thread Steven Neumann
Thank you! Steven On Thu, Oct 10, 2013 at 3:02 PM, Justin Lemkul wrote: > > > On 10/10/13 10:01 AM, Steven Neumann wrote: > >> Thanks a lot. I will use position restraints then with a strong force >> constant, no bonds and place the edge atoms within half of distance

Re: [gmx-users] carbon nanotube - ifnite in length ?

2013-10-10 Thread Steven Neumann
Thanks a lot. I will use position restraints then with a strong force constant, no bonds and place the edge atoms within half of distance between them from the box edge. Is that correct? Steven On Thu, Oct 10, 2013 at 2:56 PM, Justin Lemkul wrote: > > > On 10/10/13 9:54 AM, Steve

Re: [gmx-users] carbon nanotube - ifnite in length ?

2013-10-10 Thread Steven Neumann
Thank you. I do not have any explicit solvent in my system. I included the solvent in nonbonded parameters so not even implicit. Steven On Thu, Oct 10, 2013 at 2:47 PM, Justin Lemkul wrote: > > > On 10/10/13 9:44 AM, Steven Neumann wrote: > >> Thanks. I will place th

Re: [gmx-users] carbon nanotube - ifnite in length ?

2013-10-10 Thread Steven Neumann
Thank you. Would both be equal in terms of gaining computational time? Steven On Thu, Oct 10, 2013 at 2:46 PM, Justin Lemkul wrote: > > > On 10/10/13 9:39 AM, Steven Neumann wrote: > >> And also ... my tube has 1200 atoms and I wish to apply [ exclusions ] - >> is

Re: [gmx-users] carbon nanotube - ifnite in length ?

2013-10-10 Thread Steven Neumann
, 2013 at 2:37 PM, Justin Lemkul wrote: > > > On 10/10/13 9:34 AM, Steven Neumann wrote: > >> Thanks a lot. You the bond as a distance between atoms? I wish to avoid >> bonds as they are not necessary...just position restraint. What would be >> the force constant? I tried

Re: [gmx-users] carbon nanotube - ifnite in length ?

2013-10-10 Thread Steven Neumann
And also ... my tube has 1200 atoms and I wish to apply [ exclusions ] - is there any gmx tool to exclude interactions within the given molecule so 1 with all 1200, 2 with all 1200...etc... 1200 with all 1200? Steven On Thu, Oct 10, 2013 at 2:34 PM, Steven Neumann wrote: > Thanks a lot.

Re: [gmx-users] carbon nanotube - ifnite in length ?

2013-10-10 Thread Steven Neumann
: > > > On 10/10/13 9:21 AM, Steven Neumann wrote: > >> How about applying position restarints with a strong force constant? What >> is less computationally expensive: position restrained, freezing the whole >> molecule? The nanotube should be rigid... Shall I pla

Re: [gmx-users] carbon nanotube - ifnite in length ?

2013-10-10 Thread Steven Neumann
onic > interaction approach, both the CNT and the molecule to which it is tethered > need to be in the same [moleculetype], so run pdb2gmx (or whatever you are > using) with care. > > -Justin > > > >> Dr. Vitaly V. Chaban >> >> >> On Thu, Sep 26, 2013

[gmx-users] Steered MD

2013-10-09 Thread Steven Neumann
Dear Gromacs Users, I am trying to look for some references regarding the SMD. I found one which tells about logarithmically dependency between the Rupture force (maximum pulling force) obtained from SMD and the pulling rate. Just wondering whether you are aware (or tested) the dependency between

Re: [gmx-users] carbon nanotube - ifnite in length ?

2013-09-26 Thread Steven Neumann
oms should be (already) seen as neighboring. This looks the same as > the solvent molecule, one atom of which crossed the box boundary. > > No? > > > > Dr. Vitaly V. Chaban > > > On Thu, Sep 26, 2013 at 3:59 PM, Justin Lemkul wrote: > > > > > > On 9/26/1

[gmx-users] carbon nanotube - ifnite in length ?

2013-09-26 Thread Steven Neumann
Dear Gmx Users, I have my carbon nanotube and I wish to make it infinite in lenght. Which mdp options whall be used? pbc = xy and z is the infinite dimension? another issue: Would you apply bonds between carbon atoms within the nanotube or constraints using LINCS? Which of them is less computation

Re: [gmx-users] The optimal PME mesh load for parallel simulations is below 0.5

2013-09-04 Thread Steven Neumann
10:20 AM, Steven Neumann wrote: > >> Sorry it is a vacuum but I included implicit solvent in vdw >> parameters...So >> I need pbc as well. >> >> >> > Sorry, this doesn't make much sense to me. If you're using implicit > solvent (GB), then it's

Re: [gmx-users] The optimal PME mesh load for parallel simulations is below 0.5

2013-09-04 Thread Steven Neumann
Thank you. But with rwdv = 0 and vdw_type =User the vdw parameters will be taken into account at infinite cutoff or omitted? On Wed, Sep 4, 2013 at 3:37 PM, Justin Lemkul wrote: > > > On 9/4/13 10:35 AM, Steven Neumann wrote: > >> I am not using any solvent. I mimic the pr

Re: [gmx-users] The optimal PME mesh load for parallel simulations is below 0.5

2013-09-04 Thread Steven Neumann
Thanks a lot! On Wed, Sep 4, 2013 at 3:46 PM, Justin Lemkul wrote: > > > On 9/4/13 10:44 AM, Steven Neumann wrote: > >> Thank you. But with rwdv = 0 and vdw_type =User the vdw parameters will be >> taken into account at infinite cutoff or omitted? >> >> >

Re: [gmx-users] The optimal PME mesh load for parallel simulations is below 0.5

2013-09-04 Thread Steven Neumann
Sorry it is a vacuum but I included implicit solvent in vdw parameters...So I need pbc as well. On Wed, Sep 4, 2013 at 3:18 PM, Steven Neumann wrote: > Thank you. i am using my own vdw tables so need a cut off. > > > > > On Wed, Sep 4, 2013 at 3:13 PM, Justin Lemkul wrote: &g

Re: [gmx-users] The optimal PME mesh load for parallel simulations is below 0.5

2013-09-04 Thread Steven Neumann
Thank you. i am using my own vdw tables so need a cut off. On Wed, Sep 4, 2013 at 3:13 PM, Justin Lemkul wrote: > > > On 9/4/13 10:11 AM, Steven Neumann wrote: > >> Thank you! Would you suggest just a cut-off for coulmb? >> >> > Not a finite one. The best in

Re: [gmx-users] The optimal PME mesh load for parallel simulations is below 0.5

2013-09-04 Thread Steven Neumann
Thank you! Would you suggest just a cut-off for coulmb? Steven On Wed, Sep 4, 2013 at 3:09 PM, Justin Lemkul wrote: > > > On 9/4/13 10:03 AM, Steven Neumann wrote: > >> DEa Users, >> >> My system involves protein in vacuum - 80 atoms in box of 9x9x9 nm3. I &

[gmx-users] The optimal PME mesh load for parallel simulations is below 0.5

2013-09-04 Thread Steven Neumann
DEa Users, My system involves protein in vacuum - 80 atoms in box of 9x9x9 nm3. I want to use PME in my mdp: rcoulomb = 2.0 coulombtype = PME pme_order= 4 fourierspacing = 0.12 The cutoff needs to stay like this, I have my own tables with VDW, bonds, angles and

[gmx-users] Position Restrained of ions

2013-09-03 Thread Steven Neumann
Dear Gmx Users, i wish I could restrain some ions in my system (NA). I tried to include it in ions.itp: [ moleculetype ] ; molnamenrexcl NA1 #ifdef POSRES_NA [ position_restraints ] ; atom type fx fy fz 1 1 1000 1000 1000 #endif But it does not work. Woul

Re: [gmx-users] Umbrella Sampling

2013-07-31 Thread Steven Neumann
I am do not want to choose different pulling conditions as I build a model for specific force constant and pulling rate in given force filed. I think restraining would help much more to then exclude the ions impact. Steven On Wed, Jul 31, 2013 at 12:49 PM, Steven Neumann wrote: > Thank yo

Re: [gmx-users] Umbrella Sampling

2013-07-31 Thread Steven Neumann
Thank you a lot! On Wed, Jul 31, 2013 at 12:46 PM, Justin Lemkul wrote: > > > On 7/31/13 6:52 AM, Steven Neumann wrote: > >> But even though on the other hand that could be more realistic free energy >> which could be compare to experiment which also involves ions. W

Re: [gmx-users] Umbrella Sampling

2013-07-31 Thread Steven Neumann
But even though on the other hand that could be more realistic free energy which could be compare to experiment which also involves ions. Would Justin please comment on this? On Wed, Jul 31, 2013 at 11:46 AM, Steven Neumann wrote: > They do not dissociate...Are you sure? My mdp specifies o

Re: [gmx-users] Umbrella Sampling

2013-07-31 Thread Steven Neumann
binding? > > It would be the ligand and two ions unless the ions also at some point > discossiate from the ligand once in solvent. Could add positional restraint > for them, but dont know how that effects the calculation? > *Gesendet:* Mittwoch, 31. Juli 2013 um 09:29 Uhr &g

[gmx-users] Umbrella Sampling

2013-07-31 Thread Steven Neumann
Dear Gmx Users, I run SMD to extract the windows for US calculations. The system involves negatively charged ligand and protein. I generated the protein-ligand complex within self assembly MD simulations. I pulled my molecule away and two ions were also detached from the protein surface being att

[gmx-users] Verlet cut off scheme with tabulated potentials

2013-07-30 Thread Steven Neumann
Dear Gmx Users, I want to run my simulations with tabulated non bonded and bonded parameters on Gromacs 4.6.3. When I ask for 16 cpus it says: that I am not able to use 16 or 8 with cut off scheme Group. So I set up cutoff-scheme = Verlet verlet-buffer-drift = 0.005 ERROR 1 [file grompp.mdp]:

Re: [gmx-users] Mass in GMX topology

2013-07-11 Thread Steven Neumann
I have the same feeling, thank you. But in general mass can influence the equilibrium property so I guess yes. Steven On Thu, Jul 11, 2013 at 1:27 PM, Justin Lemkul wrote: > > > On 7/11/13 8:17 AM, Steven Neumann wrote: > >> Thank you. Another question: Does RDF depends

Re: [gmx-users] Mass in GMX topology

2013-07-11 Thread Steven Neumann
7/11/13 7:54 AM, Steven Neumann wrote: > >> Dear Users, >> >> Can you please write me where gromacs does read the mass of each atom: Is >> that the [atoms] under [moleculetype] or from [atomtypes] ? I wish to >> assign different mass for two different beads of the

[gmx-users] Mass in GMX topology

2013-07-11 Thread Steven Neumann
Dear Users, Can you please write me where gromacs does read the mass of each atom: Is that the [atoms] under [moleculetype] or from [atomtypes] ? I wish to assign different mass for two different beads of the same type in my topology. Steven -- gmx-users mailing listgmx-users@gromacs.org htt

[gmx-users] Re: Angle Distribution

2013-07-01 Thread Steven Neumann
o use it as a function which will describe the angular potential? Can I just refine and smooth it? But how I am going to calculate the derivative? Steven On Mon, Jul 1, 2013 at 8:19 AM, Steven Neumann wrote: > Dear All, > > Do you know how can I calculate angular distribution of al

[gmx-users] Angle Distribution

2013-07-01 Thread Steven Neumann
Dear All, Do you know how can I calculate angular distribution of all angles in my system? Shall I specify in one index group all e.g. 80 atoms so g_angle will calculate all possibile distributions and plot it as a sum? Steven -- gmx-users mailing listgmx-users@gromacs.org http://lists.groma

[gmx-users] Tabulated angles - how to use?

2013-06-26 Thread Steven Neumann
Dear Gmx Users, It is described how to use tabulated bonds/angles/dihedrals in 4.2.13 manual. I wish to use tables with angles table_a1.xvg, table_a2.xvg However it is not described which function to use in [ angles ]. It is said about [ bonds ] function 8 or 9 but no angles... Can anyone te

Re: [gmx-users] Can't read 1048576 bytes of 'pullx500.xvg' to compute checksum.

2013-06-25 Thread Steven Neumann
same space numbers, or > tab spacing, or etc...the only thing I noted was with xmgrace theirs a & or > something at the end, but putting even this in manually to the files did > not work...dont know if that helps... > > Stephan > > *Gesendet:* Dienstag, 25. Juni 2013 um

[gmx-users] Can't read 1048576 bytes of 'pullx500.xvg' to compute checksum.

2013-06-25 Thread Steven Neumann
Dear Users, I know this error has been discussed many times but the outcome from mdrun -pf and -px stopped at the same time which is 39470 ps. Somehow gromacs caanot read pullx500.xvg but no clue why. I tried dos2gmx and still the same error occur. As I do not care about pullx500.xvg I run grompp

Re: [gmx-users] Quartic angle potential - how to use?

2013-06-19 Thread Steven Neumann
lation. Yet another reason to avoid modifying the C > code :-) > > Mark > > On Wed, Jun 19, 2013 at 10:06 AM, Steven Neumann > wrote: > > I will speak to someone who is familiar with the code. However, my > Gromacs > > version is installed on the cluster, can I crea

Re: [gmx-users] Quartic angle potential - how to use?

2013-06-19 Thread Steven Neumann
I will speak to someone who is familiar with the code. However, my Gromacs version is installed on the cluster, can I create a file in my directory which will use this code with a given potential? On Wed, Jun 19, 2013 at 8:41 AM, Steven Neumann wrote: > > > > On Tue, Jun 18, 201

Re: [gmx-users] Quartic angle potential - how to use?

2013-06-19 Thread Steven Neumann
On Tue, Jun 18, 2013 at 6:07 PM, Mark Abraham wrote: > On Tue, Jun 18, 2013 at 5:39 PM, Steven Neumann >wrote: > > > On Mon, Jun 10, 2013 at 1:22 PM, Justin Lemkul wrote: > > > > > > > > > > > On 6/10/13 5:16 AM, Steven Neumann wrote: > >

Re: [gmx-users] Quartic angle potential - how to use?

2013-06-18 Thread Steven Neumann
On Mon, Jun 10, 2013 at 1:22 PM, Justin Lemkul wrote: > > > On 6/10/13 5:16 AM, Steven Neumann wrote: > >> Thank you. >> Do you know whether it is possible to use the 5th order polynomial for >> angles in Gromacs? I know I can use tables but wish to fit m

Re: [gmx-users] Quartic Angle

2013-06-12 Thread Steven Neumann
The constant angle has a weird units of kJ /mol rad^n while the angle is in degrees. Does not make sense but at least I know where I made a mistake... On Wed, Jun 12, 2013 at 10:43 PM, Justin Lemkul wrote: > > > On 6/12/13 10:01 AM, Steven Neumann wrote: > >> Dear Gmx Us

[gmx-users] Quartic Angle

2013-06-12 Thread Steven Neumann
Dear Gmx Users, I run a simulation with quartic angle potential. These function has two minima at 100 at 120 degrees. However, seems from my simulation that it is not applied... all angles are around 180 degrees and my chain is a line chain (straight). I have in my topology [ angles ] ; i j k

Re: [gmx-users] Nonbonded parameters

2013-06-12 Thread Steven Neumann
So how come 1830 1-4 interaction? There should be 57 of 1-4 interactions -- 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 pos

Re: [gmx-users] Nonbonded parameters

2013-06-12 Thread Steven Neumann
60 types of atoms and 60 atoms in total belonging to one residue On Wed, Jun 12, 2013 at 12:59 PM, Mark Abraham wrote: > How many atom *types* are there? > > > On Wed, Jun 12, 2013 at 1:57 PM, Steven Neumann >wrote: > > > and 1830 1-4 interactions... there are "57

Re: [gmx-users] Nonbonded parameters

2013-06-12 Thread Steven Neumann
and 1830 1-4 interactions... there are "57 all possible 1-4 interactions" On Wed, Jun 12, 2013 at 12:55 PM, Steven Neumann wrote: > "All possible interactions = 60*60/2 = 1800. Where grompp takes another 30 > from? > > > On Wed, Jun 12, 2013 at 12:52 PM, Mark Abra

Re: [gmx-users] Nonbonded parameters

2013-06-12 Thread Steven Neumann
offs... grompp is not doing > a neighbor search! > > Mark > > > On Wed, Jun 12, 2013 at 1:48 PM, Steven Neumann >wrote: > > > On Wed, Jun 12, 2013 at 12:45 PM, Mark Abraham > >wrote: > > > > > grompp is just enumerating the possible combinations

Re: [gmx-users] Nonbonded parameters

2013-06-12 Thread Steven Neumann
; > Mark I understand...but why for 60 mer polypeptide there are 1830 nonbonded interactions ? > > > On Wed, Jun 12, 2013 at 12:40 PM, Steven Neumann >wrote: > > > I got it. However, > > > > I have 60 atoms in my chain all with given C6 and C12 with a combination &g

Re: [gmx-users] Nonbonded parameters

2013-06-12 Thread Steven Neumann
An why do I have 1830 non bonded? 60*60/2 = 1800. There are 57 of 1-4 combinations so it should be lower than 1800... And why 1830 1-4 interactions as I have only 57... Please, help/ Steven On Wed, Jun 12, 2013 at 11:40 AM, Steven Neumann wrote: > I got it. However, > > I have 60 at

Re: [gmx-users] Nonbonded parameters

2013-06-12 Thread Steven Neumann
Number of degrees of freedom in T-Coupling group rest is 118.00 This run will generate roughly 2 Mb of data So if I exclude 1-3 interactions I have the same number of nonbonded parameters as well as 1-4 interactions. Can someone explain me this please? On Wed, Jun 12, 2013 at 11:12 AM,

Re: [gmx-users] Nonbonded parameters

2013-06-12 Thread Steven Neumann
On Wed, Jun 12, 2013 at 11:05 AM, Justin Lemkul wrote: > > > On 6/12/13 6:00 AM, Steven Neumann wrote: > >> Thank you. That means that e.g. when two LJ atoms of a different type >> approaching each other the non bonded LJ potential energy is a sum of two >> potenti

Re: [gmx-users] Nonbonded parameters

2013-06-12 Thread Steven Neumann
Thank you. That means that e.g. when two LJ atoms of a different type approaching each other the non bonded LJ potential energy is a sum of two potentials of those atoms? On Wed, Jun 12, 2013 at 10:44 AM, Justin Lemkul wrote: > > > On 6/12/13 5:30 AM, Steven Neumann wrote: > >&

[gmx-users] Nonbonded parameters

2013-06-12 Thread Steven Neumann
Dear Gmx Users, Please, correct me if I am wrong. In ffnonbonded.itp [ atomtypes ] - these are nonbonded parameters between atoms of the same type [ nonbond_params ] these are nonbonded parameters between atoms of different type [ pairtypes ] - 1-4 interactions If that is correct why amber force

Re: [gmx-users] [ defaults ] and nonbonded

2013-06-12 Thread Steven Neumann
Whether it would be that easy i will find it. I specified my own potential (not LJ) - so how come I should choose option 1 or 2 in nbfunct? I changed h(x) in my table 6th column. On Wed, Jun 12, 2013 at 1:18 AM, Justin Lemkul wrote: > > > On 6/11/13 4:12 PM, Steven Neumann wrote:

[gmx-users] [ defaults ] and nonbonded

2013-06-11 Thread Steven Neumann
Dear Gromacs Users, I got really confused: In manual [defualts ]: "nbfunc is the non-bonded function type. Use 1 (Lennard-Jones) or 2 (Buckingham)" I want to use mdrun -table table.xvg with my own potential, which one I should choose? "gen-pairs - is for pair generation. The default is ‘no’, i.

Re: [gmx-users] Constraints of distance

2013-06-11 Thread Steven Neumann
Thank you Justin. So using nrexcl 0 or 1 (should be the same with constraints) and type 1 will allow bonded atoms to interact (vdW) with each other? On Tue, Jun 11, 2013 at 6:10 PM, Justin Lemkul wrote: > > > On 6/11/13 10:46 AM, Steven Neumann wrote: > >> I was not sure

Re: [gmx-users] Constraints of distance

2013-06-11 Thread Steven Neumann
all is correct. > > Why are you asking? People normally report problems. > > Dr. Vitaly Chaban > > > > > > > On Tue, Jun 11, 2013 at 12:30 PM, Steven Neumann >wrote: > > > Dear Gmx Users, > > > > I am running CG simulation and I wish my beads

[gmx-users] Constraints of distance

2013-06-11 Thread Steven Neumann
Dear Gmx Users, I am running CG simulation and I wish my beads to be constraint - one away from each other of 0.4 nm. I wan to use Lincs for this purpose. I do not have any bonds in my topology or rtp entry. I just add: [ constraints ] 1 2 1 0.4 2 3 1 0.4 ... 31 32

Re: [gmx-users] Angles not read

2013-06-11 Thread Steven Neumann
On Mon, Jun 10, 2013 at 1:44 PM, Justin Lemkul wrote: > > > On 6/10/13 8:42 AM, Justin Lemkul wrote: > >> >> >> On 6/10/13 8:40 AM, Steven Neumann wrote: >> >>> On Mon, Jun 10, 2013 at 1:27 PM, Justin Lemkul wrote: >>> >>> >>

Re: [gmx-users] Angles not read

2013-06-10 Thread Steven Neumann
On Mon, Jun 10, 2013 at 2:24 PM, Justin Lemkul wrote: > > > On 6/10/13 9:22 AM, Steven Neumann wrote: > >> On Mon, Jun 10, 2013 at 2:14 PM, Justin Lemkul wrote: >> >> >>> >>> On 6/10/13 9:11 AM, Steven Neumann wrote: >>> >>> Th

Re: [gmx-users] Angles not read

2013-06-10 Thread Steven Neumann
On Mon, Jun 10, 2013 at 2:14 PM, Justin Lemkul wrote: > > > On 6/10/13 9:11 AM, Steven Neumann wrote: > >> Thank you Justing. Last question - each of my amino acid is one bead. How >> can I specify [angles ] between atoms corresponding to different residues? >> S

Re: [gmx-users] Angles not read

2013-06-10 Thread Steven Neumann
3 9:04 AM, Steven Neumann wrote: > >> On Mon, Jun 10, 2013 at 1:55 PM, Justin Lemkul wrote: >> >> >>> >>> On 6/10/13 8:53 AM, Steven Neumann wrote: >>> >>> On Mon, Jun 10, 2013 at 1:44 PM, Justin Lemkul wrote: >>>> >>>>

Re: [gmx-users] Angles not read

2013-06-10 Thread Steven Neumann
On Mon, Jun 10, 2013 at 1:55 PM, Justin Lemkul wrote: > > > On 6/10/13 8:53 AM, Steven Neumann wrote: > >> On Mon, Jun 10, 2013 at 1:44 PM, Justin Lemkul wrote: >> >> >>> >>> On 6/10/13 8:42 AM, Justin Lemkul wrote: >>> >>> >

Re: [gmx-users] Angles not read

2013-06-10 Thread Steven Neumann
On Mon, Jun 10, 2013 at 1:44 PM, Justin Lemkul wrote: > > > On 6/10/13 8:42 AM, Justin Lemkul wrote: > >> >> >> On 6/10/13 8:40 AM, Steven Neumann wrote: >> >>> On Mon, Jun 10, 2013 at 1:27 PM, Justin Lemkul wrote: >>> >>> >>

Re: [gmx-users] Angles not read

2013-06-10 Thread Steven Neumann
On Mon, Jun 10, 2013 at 1:27 PM, Justin Lemkul wrote: > > > On 6/10/13 8:23 AM, Steven Neumann wrote: > >> Dear Gmx Users, >> >> I created my own CG force field and i process my structure to pdb2gmx. I >> process 3 beads I created to check whether

[gmx-users] Angles not read

2013-06-10 Thread Steven Neumann
Dear Gmx Users, I created my own CG force field and i process my structure to pdb2gmx. I process 3 beads I created to check whether the topology is properly created: Using the CG force field in directory ./CG.ff No file 'watermodels.dat' found, will not include a water model Reading 60central_C

Re: [gmx-users] Quartic angle potential - how to use?

2013-06-10 Thread Steven Neumann
Thank you. Do you know whether it is possible to use the 5th order polynomial for angles in Gromacs? I know I can use tables but wish to fit my data into such a function. Steven On Wed, Jun 5, 2013 at 5:55 PM, Justin Lemkul wrote: > > > On 6/5/13 8:40 AM, Steven Neumann wrote: >

Re: [gmx-users] User-define non bonded table

2013-06-07 Thread Steven Neumann
values so the potential will equal zero? On Fri, Jun 7, 2013 at 2:17 PM, Steven Neumann wrote: > > > > On Fri, Jun 7, 2013 at 1:53 PM, Justin Lemkul wrote: > >> >> >> On 6/7/13 8:21 AM, Steven Neumann wrote: >> >>> Thank you. >>> Another ques

Re: [gmx-users] User-define non bonded table

2013-06-07 Thread Steven Neumann
On Fri, Jun 7, 2013 at 1:53 PM, Justin Lemkul wrote: > > > On 6/7/13 8:21 AM, Steven Neumann wrote: > >> Thank you. >> Another question... I specify in my table functions e.g. g(x) and h(x) and >> it is written that I need to setup parameters C6 and C12. But wher

Re: [gmx-users] User-define non bonded table

2013-06-07 Thread Steven Neumann
[ nonbond_params ] ; i j func V(c6) W(c12) O O 1 0.22617E-02 0.74158E-06 O OA 1 0.22617E-02 0.13807E-05 Or maybe both? On Fri, Jun 7, 2013 at 12:46 PM, Justin Lemkul wrote: > > > On 6/7/13 7:36 AM, Steven Neumann wrote: > >> Thank you, just getting into this. I just dont underst

Re: [gmx-users] User-define non bonded table

2013-06-07 Thread Steven Neumann
check the site given below: > http://www.gromacs.org/Documentation/How-tos/Tabulated_Potentials > > Thanks, > Mohan > > > On Fri, Jun 7, 2013 at 3:38 PM, Steven Neumann >wrote: > > > Dear Gmx Users, > > > > I want to specify a table to mdrun for non bond

[gmx-users] User-define non bonded table

2013-06-07 Thread Steven Neumann
Dear Gmx Users, I want to specify a table to mdrun for non bonded parameters. I wish to set it up for all atoms with same potential. Is there any example of the table like this? I want use my specific potential so table should have two columns: X distance [nm] and Y Potential energy {kJ/mol]. Is t

[gmx-users] diherdal psi vs angle teta plot

2013-06-05 Thread Steven Neumann
Dear Gmx Users, I wish to obtain the 2D scatter plot in which I will have my protein Ca-Ca-Ca tetta angle and the psi angle N(i)-CA(i)-C(k)-N(i+1) of the central residue? Would you please help? Steven -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/g

[gmx-users] Quartic angle potential - how to use?

2013-06-05 Thread Steven Neumann
Dear Gmx Users, I wish to use quartic angle potential and specify all constants. I know it is a function 6 of [ angles ] but do not know how to place my polynomial constants? Thank you, Steven -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-user

[gmx-users] CG angle distribution

2013-06-03 Thread Steven Neumann
Dear Gmx Users, I am building coarse-grained model for my protein. I group them int0 5-10 atoms. I wish to reproduce atomistic equilibrium angles as well as the spring constant. Shall I use g_angle -od option and calulate distributions? Then I will use equilibrium th0 as the centre of the gaussian

[gmx-users] FEP of ion binding protein

2013-05-27 Thread Steven Neumann
Dear Gmx Users, I wish to compare binding free energy obtained from US which I have sucsessfully conducted and the one obtained using FEP. Would you suggest any tutorial? Steven -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users * Please search

[gmx-users] RDF - water and protein

2013-05-13 Thread Steven Neumann
Dear Gmx Users, I run long simulation of my protein with 50 small molecules in water. I calculated the RDF (Protein - Water) using -surf mol and -rdf mol_com. Please, take a look at my plot: http://speedy.sh/tmJbD/rdf-P-W.png Could you please, explain me why the second peak is so high? Shall I u

[gmx-users] Re: Illegal instruction (core dumped) - trjconv

2013-04-29 Thread Steven Neumann
On Mon, Apr 29, 2013 at 10:11 AM, Steven Neumann wrote: > Dear Gmx Users, > > I produced my trajectory using Gromacs 4.6 on GPUs. > When I try: > > trjconv -s md298_gpu.tpr -f md298_gpu.xtc -pbc mol -ur compact -o > md298_gpuURcomp.xtc > > Select a group: 0 > Select

[gmx-users] Illegal instruction (core dumped) - trjconv

2013-04-29 Thread Steven Neumann
Dear Gmx Users, I produced my trajectory using Gromacs 4.6 on GPUs. When I try: trjconv -s md298_gpu.tpr -f md298_gpu.xtc -pbc mol -ur compact -o md298_gpuURcomp.xtc Select a group: 0 Selected 0: 'System' Reading frame 0 time0.000 Precision of md298_gpu.xtc is 0.001 (nm) Using output p

Re: [gmx-users] Re: SMD - reproducibility

2013-04-26 Thread Steven Neumann
Thanks anyway. I will try tu pull it 20-40 times to tell something. Steven On Fri, Apr 26, 2013 at 1:23 PM, Thomas Schlesier wrote: > Don't know. One idea i have: Take a flexible and a relative rigid system > and perform simulations with the same starting conditions (-> using -t > *.cpt). I woul

Re: [gmx-users] Re: SMD - reproducibility

2013-04-26 Thread Steven Neumann
Thanks for this. I think option 2 is more reasonable. However, still do not know why I get sometimes 3 types of profiels and sometimes 10 for 10 SMD simulations... On Fri, Apr 26, 2013 at 12:46 PM, Thomas Schlesier wrote: > Think i now understand your question. Forget what i wrote before. > I cou

[gmx-users] SMD - reproducibility

2013-04-26 Thread Steven Neumann
Dear Users, I am running my puling simulations of ligand with constant velocity. First I minimize and equilibrate my system: grompp -f minim.mdp -c Solvions.gro -p topol.top -o em.tpr mdrun -s em.tpr -deffnm em grompp -f nvt298US.mdp -n index.ndx -c em.gro -p topol.top -o nvt298.tpr mdrun -s nvt2

Re: [gmx-users] Interaction across PBC

2013-04-23 Thread Steven Neumann
Thanks a lot Steven On Tue, Apr 23, 2013 at 2:18 PM, Justin Lemkul wrote: > > > On 4/23/13 9:15 AM, Steven Neumann wrote: > >> Shall I specify one index group for two regions or 2 seprate? g_mindist >> asks just for one group. >> >> > If it only t

[gmx-users] which tool returns vectors?

2013-04-23 Thread Steven Neumann
Dear Users, Does any one know which command is capable to return the vector of a specified group of 2 atoms (e.g. C=O in protein) over the simulation time? Steven -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users * Please search the archive at

Re: [gmx-users] Interaction across PBC

2013-04-23 Thread Steven Neumann
Shall I specify one index group for two regions or 2 seprate? g_mindist asks just for one group. would twice as cutoff would be sufficent to assess they do not interact? On Tue, Apr 23, 2013 at 1:54 PM, Steven Neumann wrote: > Thank you. > > Steven > > > On Tue, Apr 23, 2013

Re: [gmx-users] Interaction across PBC

2013-04-23 Thread Steven Neumann
Thank you. Steven On Tue, Apr 23, 2013 at 1:23 PM, Justin Lemkul wrote: > > > On 4/23/13 6:37 AM, Steven Neumann wrote: > >> Dear Gmx users, >> >> >> My protien has got some strong acidic and strong basic parts. I fold and >> unfold my protein w

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

2013-04-23 Thread Steven Neumann
Thanks. Which tool would provide me vectors over a time? Steven On Thu, Apr 18, 2013 at 8:17 PM, Mark Abraham wrote: > 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 >wrote: > > >

[gmx-users] Interaction across PBC

2013-04-23 Thread Steven Neumann
Dear Gmx users, My protien has got some strong acidic and strong basic parts. I fold and unfold my protein with different temperaturss. I bserved high affinity of those regions towards each other, they are very close to each other over the simulation. How can I possibly check whether my two regi

Re: [gmx-users] trjcat set timestep/amb2gmx?

2013-04-17 Thread Steven Neumann
I read it... did not know styles are equal to points :) On Wed, Apr 17, 2013 at 4:08 PM, Justin Lemkul wrote: > > > On 4/17/13 11:06 AM, Steven Neumann wrote: > >> What do you mean by style 8? >> >> > Did you read trjconv -h? There's an enumerated

Re: [gmx-users] trjcat set timestep/amb2gmx?

2013-04-17 Thread Steven Neumann
d, but it could be more useful using as x the > frame > > number. > > This could, at least, permitting to visualize the plot with xmgrace > without > > postprocessing > > the .xvg > > > > > > > > Francesco > > > > > > 2013

[gmx-users] Autocorrelation of C-N

2013-04-17 Thread Steven Neumann
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.o

[gmx-users] run.mdp on 4.6-gpu

2013-04-17 Thread Steven Neumann
Dear Gmx Users, Could you please assess my mdp file on GPUs: title = Protein-ligand complex MD simualation ; Run parameters integrator = md; leap-frog integrator nsteps = 1; 200 ns dt = 0.002 ; 2 fs ; Output control nstxout = 0 ; suppress .trr

Re: [gmx-users] trjcat set timestep/amb2gmx?

2013-04-16 Thread Steven Neumann
Got this. Thank you On Tue, Apr 16, 2013 at 11:21 AM, Justin Lemkul wrote: > > > On 4/16/13 5:32 AM, Steven Neumann wrote: > >> Thank you all. >> >> Would you suggest any link to a script like amb2gmx.pl ? I cannot fina >> anything like this - basically I w

Re: [gmx-users] trjcat set timestep/amb2gmx?

2013-04-16 Thread Steven Neumann
s could, at least, permitting to visualize the plot with xmgrace without > postprocessing > the .xvg > > > > Francesco > > > 2013/4/16 Mark Abraham > > > On Apr 15, 2013 6:27 PM, "Justin Lemkul" wrote: > > > > > > > > > &g

[gmx-users] trjcat set timestep/amb2gmx?

2013-04-15 Thread Steven Neumann
Dear Gmx Users, I obtained dcd trajectory from simulation in another software. I wish to merge many trajectories using trjcat with a proper timestep. Is that option possible using gromacs or shall use a script to produce tpr file from my prmtop file e.g. amb2gmx ? Any links for such a script? tha

Re: [gmx-users] Peptide distance matrix

2013-04-15 Thread Steven Neumann
output is a symmetrical matrix of smallest distances between all residues. And nothing else... Steven On Mon, Apr 15, 2013 at 4:24 PM, Steven Neumann wrote: > Thank you for this. > > Steven > > On Mon, Apr 15, 2013 at 4:22 PM, Justin Lemkul wrote: > > > > >

Re: [gmx-users] Peptide distance matrix

2013-04-15 Thread Steven Neumann
Thank you for this. Steven On Mon, Apr 15, 2013 at 4:22 PM, Justin Lemkul wrote: > > > On 4/15/13 11:08 AM, Steven Neumann wrote: >> >> On Mon, Apr 15, 2013 at 3:55 PM, Justin Lemkul wrote: >>> >>> >>> >>> On 4/15/13 10:06 AM, Steven N

Re: [gmx-users] Peptide distance matrix

2013-04-15 Thread Steven Neumann
And another question: I want to analyze all 5 nanoseconds every 20 ns of my trajectory. Would you suggest using trjcat to create one trajectory first and then process to g_mdmat or can I specify time periods (frames) I wish to analyze? Steven On Mon, Apr 15, 2013 at 4:08 PM, Steven Neumann

Re: [gmx-users] Peptide distance matrix

2013-04-15 Thread Steven Neumann
On Mon, Apr 15, 2013 at 3:55 PM, Justin Lemkul wrote: > > > On 4/15/13 10:06 AM, Steven Neumann wrote: >> >> On Mon, Apr 15, 2013 at 3:01 PM, Justin Lemkul wrote: >>> >>> >>> >>> On 4/15/13 9:59 AM, Steven Neumann wrote: >>>&g

Re: [gmx-users] Peptide distance matrix

2013-04-15 Thread Steven Neumann
On Mon, Apr 15, 2013 at 3:01 PM, Justin Lemkul wrote: > > > On 4/15/13 9:59 AM, Steven Neumann wrote: >> >> Dear Gmx Users, >> >> I want to calculate a distance matrix of each amino acid (1, 2, ...25) >> averaged over simulation time with all amino acids. S

[gmx-users] Peptide distance matrix

2013-04-15 Thread Steven Neumann
Dear Gmx Users, I want to calculate a distance matrix of each amino acid (1, 2, ...25) averaged over simulation time with all amino acids. So matrix of 25x25: 1) is there a tool which can do this or just the use of g_dist fof 600 (25x25 - 25) times? 2) Would you recommend any nice visualisation t

Re: [gmx-users] Desired density - genbox

2013-03-19 Thread Steven Neumann
Thanks. So in this case no matter what density I will start with e.g. 480 kg/m3 presuming the force filed is correct I should get at given conditions the density of interest? Steven On Mon, Mar 18, 2013 at 10:34 PM, Justin Lemkul wrote: > > > On 3/18/13 6:14 PM, Steven Neumann wrote:

Re: [gmx-users] Desired density - genbox

2013-03-18 Thread Steven Neumann
On Mon, Mar 18, 2013 at 8:40 PM, Justin Lemkul wrote: > > > On 3/18/13 2:56 PM, Steven Neumann wrote: >> >> Dear Gmx Users, >> >> I am trying to obtain given density for my system for a given molecule >> - its a cubic box of 5 nm in dimension. >> >

  1   2   3   4   >