[gmx-users] sequence of atoms different for rmsd reference structure

2014-09-17 Thread gromacs query
Hi All I have one PDB structure which I want to use as reference for rmsd calculation (with best fitting). But I realize in my gromacs trajectory the same structure (if written in PDB) has different sequence of atoms and some different atom names (especially H, OP1 vs O1P and O2P vs OP2) compared

[gmx-users] make_ndx bug?

2014-09-19 Thread gromacs query
Hi All I have some pdb file and I removed hydrogens. When I try to make ndx file using some residue option then rather printing atom numbers in that residue it prints line number where that residue starts and ends. make_ndx -f no_H.pdb -o test.ndx selecting as: r 300-330 Then in ndx file I get

Re: [gmx-users] make_ndx bug?

2014-09-19 Thread gromacs query
g. if you use r 2 then it complains that residue is missing but I think atom numbers are not read from 'atom number column' itself but as line numbers On Fri, Sep 19, 2014 at 9:37 PM, Justin Lemkul wrote: > > > On 9/19/14 8:43 AM, gromacs query wrote: > >> Hi All >>

Re: [gmx-users] make_ndx bug?

2014-09-19 Thread gromacs query
n such case. best regards On Fri, Sep 19, 2014 at 10:20 PM, Justin Lemkul wrote: > > > On 9/19/14 5:15 PM, gromacs query wrote: > >> Hi Justin >> >> Hope it makes sense now :) Here is the link (small piece of some pdb >> without hydrogens) >> >>

[gmx-users] can not split pdb beyond 1019

2014-10-06 Thread gromacs query
Hi All, I am trying to split my trajectory in single frame pdbs for some analysis using trjconv but it ends with an error after ouput 1019 pdbs. I tried with different version of gromacs (4.6.3 and 5.0) its the same error. trjconv -s ../em.tpr -f ../image_fitted.trr -o test.pdb -split 10 -b 4

[gmx-users] trr file time/frame information

2014-12-08 Thread gromacs query
Dear All, I combined the splitted trajectories (trr) from 50ns to 100ns and got combined trajectoty trr file. I am wondering whether the new combined trajectory will start from 0ns or it will retain time/frame information which 50ns in this case? Thanks JIom -- Gromacs Users mailing list * Plea

Re: [gmx-users] trr file time/frame information

2014-12-08 Thread gromacs query
ems no more supported) Thanks On Mon, Dec 8, 2014 at 3:36 PM, Mark Abraham wrote: > Hi, > > That depends which tool you used, how you used it, etc. You should use gmx > check and perhaps gmx dump to see what is there. > > Mark > > On Mon, Dec 8, 2014 at 4:31 PM, gromacs

[gmx-users] constraining charges

2015-10-17 Thread gromacs query
Hi All, I have some amino-acid residue in which some extra side chain(s) attached by replacing H from amino-acid. I want to simulate this with CHARMM ff. Somehow I can choose best bonded params. But in order to keep integral charge on this modified aminoacid I was thinking if there is any way to

[gmx-users] parameter for repeating units with cgenff

2015-12-11 Thread gromacs query
Hi All, I am using cgenff to derive parameters for some polymer (A-...-C). The repeating block B has some Oxygen O atom which joins the next residue (-BO-). When I derive charges for B residue with cgenff then I get net -1 charge on B residue because of this unsatsified O atom. I want to keep

Re: [gmx-users] parameter for repeating units with cgenff

2015-12-11 Thread gromacs query
al. best regards, Vishal On Fri, Dec 11, 2015 at 4:43 PM, Justin Lemkul wrote: > > > On 12/11/15 10:17 AM, gromacs query wrote: > >> Hi All, >> >> I am using cgenff to derive parameters for some polymer (A-...-C). The >> repeating block B has some Oxyg

Re: [gmx-users] parameter for repeating units with cgenff

2015-12-11 Thread gromacs query
caps in topology will give non-integral charge or in other words this part -CO-CH(X)-NH- will always be non-integral. Sorry am confused. Thanks On Fri, Dec 11, 2015 at 5:50 PM, Justin Lemkul wrote: > > > On 12/11/15 12:17 PM, gromacs query wrote: > >> Hi Justin >> >>

Re: [gmx-users] parameter for repeating units with cgenff

2015-12-11 Thread gromacs query
c 11, 2015 at 6:33 PM, Justin Lemkul wrote: > > > On 12/11/15 1:28 PM, gromacs query wrote: > >> Hi Justin >> >> I am bit lost I think. For e.g. in the amino acid residues library we have >> -CO-CH(R)-NH- which we can combine these amino acids in any way and in t

Re: [gmx-users] parameter for repeating units with cgenff

2015-12-12 Thread gromacs query
have tried different combinations but this is the minimal formal charge I was able to achieve. Thanks again for help. On Fri, Dec 11, 2015 at 7:23 PM, Justin Lemkul wrote: > > > On 12/11/15 2:20 PM, gromacs query wrote: > >> Hi Justin >> >> Thanks for the explana

Re: [gmx-users] parameter for repeating units with cgenff

2015-12-12 Thread gromacs query
Due to reasons, would it be possible for you to suggest off-list? I will really appreciate. Thanks. On Sat, Dec 12, 2015 at 6:35 PM, Justin Lemkul wrote: > > > On 12/12/15 1:25 PM, gromacs query wrote: > >> Hi Justin >> >> Somehow I can not use available units in

Re: [gmx-users] radial density profile

2015-12-19 Thread gromacs query
Hi, Am wondering how you can have an 'axis' for single stranded DNA simulation? ssDNA should collapse quickly if starting from helical type configuration unless you are considering only one strand from ds-DNA or some restraints or very small ss fragment. If you are using RDF then I assume your si

[gmx-users] identifying terminal residue error

2015-12-20 Thread gromacs query
Hi All I have a modified base in normal DNA chain. Thus I added that residue in rtp file and updated residuetypes.dat as well (CHARMM36). When I use pdb2mx with -ter interactive mode then it wrongly identifies the ending terminal at the residue just before the modified base. e.g. : : TER DC (169

Re: [gmx-users] identifying terminal residue error

2015-12-20 Thread gromacs query
Ah! Thanks Justin, I just need extra coffee. Required file was not in working directory. On Sun, Dec 20, 2015 at 9:17 PM, Justin Lemkul wrote: > > > On 12/20/15 4:13 PM, gromacs query wrote: > >> Hi All >> >> I have a modified base in normal DNA chain. Thus I add

[gmx-users] COM motion removal: diffusion analysis

2015-12-28 Thread gromacs query
Hi All I have protein-membrane system and gromacs tutorial by Justin suggests to do COM separately for protein_membrane and water. I am bit meesed up COM removal understanding. 1) I want to do diffusion analysis of protein in a membrane (2D/3D) how it will influence my results if I perform it o

Re: [gmx-users] COM motion removal: diffusion analysis

2015-12-28 Thread gromacs query
>> It is not a matter of resetting positions; it is about subtracting falsely accumulated kinetic energy. Thanks, this is it :) On Mon, Dec 28, 2015 at 8:58 PM, Justin Lemkul wrote: > > > On 12/28/15 10:28 AM, gromacs query wrote: > >> Hi All >> >> >>

[gmx-users] charmm36 ions fix

2016-06-20 Thread gromacs query
Hi All, I am using charmm36 downloaded from (MacKerell website). Somehow I was enlightened while reading some papers applying NBFIX by Roux et al. I can see there is indeed a nbfix.itp file in CHARMM36 (not 27) and used as default (via forcefield.itp) but the file does not mention any reference wh

Re: [gmx-users] charmm36 ions fix

2016-06-20 Thread gromacs query
PM, Justin Lemkul wrote: > > > On 6/20/16 12:30 PM, gromacs query wrote: > >> Hi All, >> >> I am using charmm36 downloaded from (MacKerell website). Somehow I was >> enlightened while reading some papers applying NBFIX by Roux et al. I can >> see there is i

Re: [gmx-users] charmm36 ions fix

2016-06-20 Thread gromacs query
ne-protein simulations. Thanks for any suggestions. JIomm On Mon, Jun 20, 2016 at 6:23 PM, Justin Lemkul wrote: > > > On 6/20/16 1:19 PM, gromacs query wrote: > >> Thanks Justin, basically am stuck with an old question of using high salt >> conc. with membrane simulations a

Re: [gmx-users] charmm36 ions fix

2016-06-21 Thread gromacs query
Hi Justin >> Follow the analysis of the Venable paper I think you mean this paper: http://pubs.acs.org/doi/abs/10.1021/jp401512z JIomm -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read

[gmx-users] translate problem

2015-03-02 Thread gromacs query
Hi All, I have membrane thickness along Z and want to merge peptide in membrane. I need to translate peptide (say 10 Angs) along -Z to adjust its position and used command like this: editconf -f pep.pdb -c -princ -rotate 0 90 0 -translate -10 0 0 -o new_pep.pdb It seems that it does not move at

Re: [gmx-users] translate problem

2015-03-02 Thread gromacs query
your > editconf command. > Your command seems to translate along x. > > Chandan > > On Mon, Mar 2, 2015 at 5:16 PM, gromacs query > wrote: > > > Hi All, > > > > I have membrane thickness along Z and want to merge peptide in membrane. > I > > nee

Re: [gmx-users] translate problem

2015-03-02 Thread gromacs query
Hi Justin, >> If that doesn't work, trjconv can translate Thanks a lot, trjconv works fine. For some reason I could not achieve this with editconf. regards, JIom On Mon, Mar 2, 2015 at 1:48 PM, Justin Lemkul wrote: > > > On 3/2/15 6:46 AM, gromacs query wrote: >

[gmx-users] gmx_mtop_atomlookup_settle_init called without settles

2015-03-02 Thread gromacs query
Hi, I am using g_membed (for CG) to insert peptide with following steps: 1) grompp generates a tpr file; with no issues 2) g_membed generates membed.dat; no issues 3) lastly mdrun gives this error: gmx_mtop_atomlookup_settle_init called without settles Justin has suggested this may arise due

[gmx-users] Error: see option -rdd, for pairs and tabulated bonds also see option -ddcheck

2015-03-09 Thread gromacs query
Hi All, I am doing CG membrane with protein. After running long simulation I am getting this error: Fatal error: 2 of the 13319 bonded interactions could not be calculated because some atoms involved moved further apart than the multi-body cut-off distance (1.4 nm) or the two-body cut-off distanc

Re: [gmx-users] Error: see option -rdd, for pairs and tabulated bonds also see option -ddcheck

2015-03-10 Thread gromacs query
5 at 5:23 PM, Mark Abraham wrote: > Hi, > > Very likely your run is blowing up from (long-time) instability - two atoms > have moved some crazy amount. You can probably confirm that from the time > series of some observables. > > Mark > > On Mon, Mar 9, 2015 at 12:54 PM,

[gmx-users] no error on missing cpt file

2015-04-29 Thread gromacs query
Dear All, mdrun is not giving any error about missing cpt file. It runs using tpr from initial time zero. Sometimes my job get killed and I need to use cpt file in some script but if cpt is not found my job starts running from zero time in tpr file which is harmful as it consumes my processors to

Re: [gmx-users] no error on missing cpt file

2015-04-29 Thread gromacs query
e names as such without any compromise. regards, Jiom On Wed, Apr 29, 2015 at 12:36 PM, Mark Abraham wrote: > On Wed, Apr 29, 2015 at 1:20 PM Kutzner, Carsten wrote: > > > Hi JIom, > > > > > On 29 Apr 2015, at 11:50, gromacs query > wrote: > > > >

[gmx-users] all-atom pca

2015-05-11 Thread gromacs query
Hi All I am doing PCA and its working fine if I choose few atoms only (heavy) and I get extreme pdb file having heavy atoms for, say, eigenvector 1. But I want to have extremes of all-atom structures, so I thought of doing PCA on all-atoms but I end up with having Segmentation default error which

Re: [gmx-users] all-atom pca

2015-05-11 Thread gromacs query
; Tsjerk > > On Mon, May 11, 2015 at 10:32 PM, gromacs query > wrote: > > > Hi All > > > > I am doing PCA and its working fine if I choose few atoms only (heavy) > and > > I get extreme pdb file having heavy atoms for, say, eigenvector 1. But I > > want t

Re: [gmx-users] all-atom pca

2015-05-11 Thread gromacs query
erk Wassenaar wrote: > Hi Jio, > > You'll have to look at the projections (-proj), and see at which time these > have an extreme value. > > Cheers, > > Tsjerk > > On Mon, May 11, 2015 at 11:00 PM, gromacs query > wrote: > > > Hi Tsjerk > > > &g

Re: [gmx-users] all-atom pca

2015-05-12 Thread gromacs query
re not > real structures. But you wanted the corresponding all-atom ones. > > Cheers, > > Tsjerk > On May 12, 2015 00:45, "gromacs query" wrote: > > > Hi Tsjerk > > > > In -proj plot of v1 vs Time I used extreme values of V1 and saw > structures >

[gmx-users] Fatal error: You selected groups with differen number of atoms.

2015-06-22 Thread gromacs query
Hi All I want to align small residues based on 3 atom names only (ignoring rest of the atoms) and confrms does not require same number of atoms: first.pdb with 7 atoms and second.pdb with 3 atoms. echo "0" "0"| gmx confrms -f1 first.pdb.pdb -f2 second.pdb -no -name -o fit.pdb I am getting error

Re: [gmx-users] Fatal error: You selected groups with differen number of atoms.

2015-06-24 Thread gromacs query
Hi All Am just thinking if my mail was spammed. Thanks for any help. On Mon, Jun 22, 2015 at 4:58 PM, gromacs query wrote: > Hi All > > I want to align small residues based on 3 atom names only (ignoring rest > of the atoms) and confrms does not require same number of atoms: firs

[gmx-users] polymer in octahedron box pbc problem

2015-10-01 Thread gromacs query
Hi All, I have linear polymer in octahedron box and when I try to visualise a polymer in octahedraon box then water octahedron box looks as expected but my polymer looks broken. trjconv_avx -f input.xtc -o out.xtc -pbc mol -ur compact -s input.tpr -b 10 I have tried to use -center (whole po

Re: [gmx-users] polymer in octahedron box pbc problem

2015-10-01 Thread gromacs query
Hi Victor, With 'pbc whole' polymer remains broken in pbc and the octahedron box looks weird too. JIom On Thu, Oct 1, 2015 at 2:52 PM, Victor Rosas Garcia wrote: > Have you tried the option "-pbc whole" instead of "-pbc mol"? > > Victor > > 2015-

Re: [gmx-users] polymer in octahedron box pbc problem

2015-10-01 Thread gromacs query
Hi, I think its solved, I used -center option but now rather selecting full polymer I selected some of the central residue in polymer chain. Now box looks good and polymer is not broken. JIom On Thu, Oct 1, 2015 at 3:43 PM, gromacs query wrote: > Hi Victor, > > With 'pbc

Re: [gmx-users] charmm36 ions fix

2016-06-21 Thread gromacs query
Thanks Justin, your inputs are really helpful. I did basic analysis on thickness which increased slightly with salt conc. (with NBFix) but I do need to see other things more carefully. JIomm On Tue, Jun 21, 2016 at 12:28 PM, Justin Lemkul wrote: > > > On 6/21/16 6:29 AM, gromacs qu

[gmx-users] rotation while pulling

2017-03-28 Thread gromacs query
Hi All, {apologies if posted multiple times; seems like some trouble with my account} I want to pull linear molecule across the membrane in two ways allowing: - translation along Z and rotation around Z - translation along Z and rotation in all axes My current mdp-setting pulls molecular along Z

[gmx-users] xtc with perl

2017-06-30 Thread gromacs query
Hi all, I would like my perl script to analyse xtc trajectory file directly instead of analysing human readable extracted frame files (pdb/gro etc.) to save disk space. I could not find much on mailing list but there was a brief reference about using Swig which I could not understand how to start

Re: [gmx-users] xtc with perl

2017-06-30 Thread gromacs query
t; files, or using the GROMACS C++ analysis tools template, rather than > wasting time writing your own code to read them. If you really must, then > you should wrap this library > http://www.gromacs.org/Developer_Zone/Programming_Guide/XTC_Library. > > Mark > > On Fri, Jun

[gmx-users] pullx (gromacs 4.6 vs 2016)

2017-08-25 Thread gromacs query
Hi All, I am using 'equivalent' mdp settings for pull run in Gromacs 4.6 and 2016. With gromacs 4.6 I get time, 0Z and 1dZ; which I can understand as explained here: https://www.mail-archive.com/gmx-users@gromacs.org/ msg22453.html With Gromacs 2016 I get only time and some distance and xvg fil

Re: [gmx-users] pullx (gromacs 4.6 vs 2016)

2017-08-29 Thread gromacs query
Hi All, Am just wondering if somebody has seen this in different gmx versions? Thanks. JIom On Fri, Aug 25, 2017 at 11:00 AM, gromacs query wrote: > Hi All, > > I am using 'equivalent' mdp settings for pull run in Gromacs 4.6 and 2016. > > With gromacs 4.6 I get time

[gmx-users] performance

2017-09-13 Thread gromacs query
Hi All, I have some issues with gromacs performance. There are many nodes and each node has number of gpus and the batch process is controlled by slurm. Although I get good performance with some settings of number of gpus and nprocs but when I submit same job twice on the same node then the perfor

Re: [gmx-users] performance

2017-09-13 Thread gromacs query
Sorry forgot to add; we thought the two jobs are using same GPU ids but cuda visible devices show both jobs are using different ids (0,1 and 2,3) - J On Wed, Sep 13, 2017 at 5:33 PM, gromacs query wrote: > Hi All, > > I have some issues with gromacs performance. There are many nodes

Re: [gmx-users] performance

2017-09-13 Thread gromacs query
t; -- > Szilárd > > > On Wed, Sep 13, 2017 at 6:35 PM, gromacs query > wrote: > > Sorry forgot to add; we thought the two jobs are using same GPU ids but > > cuda visible devices show both jobs are using different ids (0,1 and 2,3) > > > > - > > J >

Re: [gmx-users] performance

2017-09-13 Thread gromacs query
your mdrun launch command. > > > I hope one of these will deliver the desired results :) > > Cheers, > -- > Szilárd > > > On Wed, Sep 13, 2017 at 7:47 PM, gromacs query > wrote: > > Hi Szilárd, > > > > Thanks for your reply. This is useful but now am think

Re: [gmx-users] performance

2017-09-13 Thread gromacs query
iably less than 50%) as expected from a single independent job. Now am wondering if its still related to overlap of cores as pin on should lock the cores for the same job. -J On Wed, Sep 13, 2017 at 7:33 PM, gromacs query wrote: > Hi Szilárd, > > Thanks, option 3 was in my mind but I need

Re: [gmx-users] performance

2017-09-13 Thread gromacs query
a node? > > How many threads/cores does you node have? Can you share log files? > > Cheers, > -- > Szilárd > > > On Wed, Sep 13, 2017 at 9:14 PM, gromacs query > wrote: > > Hi Szilárd, > > > > Sorry I was bit quick to say its working with pinoffset. I j

Re: [gmx-users] performance

2017-09-14 Thread gromacs query
M, Szilárd Páll wrote: > On Wed, Sep 13, 2017 at 11:14 PM, gromacs query > wrote: > > Hi Szilárd, > > > > Thanks again. I tried now with -multidir like this: > > > > mpirun -np 16 gmx_mpi mdrun -s test -ntomp 2 -maxh 0.1 -multidir t1 t2 > t3 t4 > > &g

Re: [gmx-users] performance

2017-09-14 Thread gromacs query
t convenient > sometimes when job crashes in the middle and automatic restart from cpt > file would be difficult. > > -J > > > On Thu, Sep 14, 2017 at 11:26 AM, Szilárd Páll > wrote: > >> On Wed, Sep 13, 2017 at 11:14 PM, gromacs query >> wrote: >> > Hi S

[gmx-users] documentation error?

2017-09-15 Thread gromacs query
Hi All, Though am naive, and much appreciate the efforts required to put things altogether but I think I found few things confusing or wrong by mistake or repeated here: http://manual.gromacs.org/documentation/2016/user-guide/mdrun-performance.html The replaced corrections are in [CAPITAL] (plea

Re: [gmx-users] performance

2017-09-18 Thread gromacs query
Hi Szilárd, {I had to trim the message as my message is put on hold because only 50kb allowed and this message has reached 58 KB! Not due to files attached as they are shared via dropbox}; Sorry seamless reading might be compromised for future readers. Thanks for your replies. I have shared log f

Re: [gmx-users] charmm 36 force field for DNA

2017-09-21 Thread gromacs query
Hi Qinghua, I am not sure about any too but If you have pdb working fine with AMBER (as you said) then its not much of work doing it 'manually' ; takes less than two minutes. You can do this: 1) First remove hydrogens and let gromacs add by itself at the final step (matching from charm36). You c

Re: [gmx-users] performance

2017-09-21 Thread gromacs query
and how to lay out jobs in a node to > get best throughput. Start with run configs testing settings with > -multi to avoid pinning headaches and fill at least half a node (or a > full node) with #concurrent simulations >= #GPUs. > I will see if I get some node free. I need to wait. Th

[gmx-users] AMBER ff10 with Gromacs

2013-11-26 Thread gromacs query
Dear All I never used Gromacs for nucleic acids. In AMBER there is ff10 (ff99SB+parmbasco+ also one needs to add ions according to water used by Thomas Cheatham et al.) does gromacs include this force field (AMBER ff implemented in gromacs)? regards JIo -- Gromacs Users mailing list * Please s

[gmx-users] flags for internal coordinates restraints

2013-11-26 Thread gromacs query
Dear All I want to simulate two polymer chains as rigid bodies, in other words allowing two polymer to interact with each other but keeping their internal coordinates restrained. Also it can be seen as two rods interacting with each other. Which flags/options can be used in Gromacs? regards JIom

Re: [gmx-users] AMBER ff10 with Gromacs

2013-11-28 Thread gromacs query
ilárd Páll wrote: > >> > >> Hi, > >> > >> If you look at share/gromacs/top/ in the GROMACS installation > >> directory you can see which FF-s are included and amberff10 is not > >> there, so the answer is no. > >> > >> Cheers,

Re: [gmx-users] flags for internal coordinates restraints

2013-12-04 Thread gromacs query
LINCS or SHAKE (all atoms) can be used for this? or some better options available? it can be seen as two protein helix interacting with each other as rigid bodies... thanks On Wed, Nov 27, 2013 at 12:35 AM, gromacs query wrote: > Dear All > > I want to simulate two polymer chains

[gmx-users] DNA with charmm27 works fine but why

2014-02-04 Thread gromacs query
Hi All I have built DNA with NAB in AMBER and it provide names as D*, D*3, D*5 (where * = A,T,G,C) and want to use charmm27 in gromacs. The charmm27.ff/dna.rtp does not have D*5 or D*3 and has only D* residues, so I renamed D*5 or D*3 to just corresponding D* residues. The output pdb with pdb2gmx

Re: [gmx-users] DNA with charmm27 works fine but why

2014-02-04 Thread gromacs query
d use some .top file? (and if I understood it correctly; if I am using from rtp file, which is a sort of residue library file, then order does NOT matter) thanks On Tue, Feb 4, 2014 at 5:03 PM, Justin Lemkul wrote: > > > On 2/4/14, 11:50 AM, gromacs query wrote: > >> Hi All >

Re: [gmx-users] DNA with charmm27 works fine but why

2014-02-04 Thread gromacs query
some common top file then even I can not use same .top file? On Tue, Feb 4, 2014 at 5:39 PM, Justin Lemkul wrote: > > > On 2/4/14, 12:35 PM, gromacs query wrote: > >> Hi Justin >> >> Presence of the correct atoms is the only absolute requirement >>>> &g

[gmx-users] dna with charmm36

2014-02-05 Thread gromacs query
Hi All I have simple DNA pdb which has just 5 DA residues built from NAB in AMBER. I am trying to use pdb2gmx with charmm36-jan2014.ff downloaded from: http://mackerell.umaryland.edu/CHARMM_ff_params.html pdb2gmx -f a.pdb -o gro.pdb -ff charmm36-jan2014 I have changed all names according to

Re: [gmx-users] dna with charmm36

2014-02-05 Thread gromacs query
Hi Justin Thanks it worked, just to add I have to use -ignh as well. I think it does not detect some H at 5 end built from AMBER NAB programme. Jiom On Wed, Feb 5, 2014 at 3:13 PM, Justin Lemkul wrote: > > > On 2/5/14, 9:49 AM, gromacs query wrote: > >> Hi All >>

[gmx-users] No default U-B types charmm

2014-03-03 Thread gromacs query
Hi All I am trying to use charmm36 (charmm36-jan2014 from charmm website) with popc membrane built using charmm-gui (have water and ions). I used commands as follows: pdb2gmx -f step5_assembly.pdb -o popc.gro -water tip3p -ff charmm36-jan2014 editconf -f popc.gro -o popc_box.gro -c -d 0.0 grompp

Re: [gmx-users] No default U-B types charmm

2014-03-03 Thread gromacs query
error. On Mon, Mar 3, 2014 at 5:46 PM, Justin Lemkul wrote: > > > On 3/3/14, 10:57 AM, gromacs query wrote: > >> Hi All >> >> I am trying to use charmm36 (charmm36-jan2014 from charmm website) with >> popc membrane built using charmm-gui (have water and ion

Re: [gmx-users] No default U-B types charmm

2014-03-03 Thread gromacs query
and trying, will let you know; I am newbie to gromacs) thanks On Mon, Mar 3, 2014 at 6:48 PM, Mark Abraham wrote: > On Mon, Mar 3, 2014 at 7:25 PM, gromacs query >wrote: > > > Hi Justin > > > > Its problem with the waters I think. > > > Can't be, they

Re: [gmx-users] dna with charmm36

2014-03-03 Thread gromacs query
t works fine and all Hs are added without error. thanks JIom On Wed, Feb 5, 2014 at 3:53 PM, gromacs query wrote: > need a coffee! > > > On Wed, Feb 5, 2014 at 3:44 PM, Justin Lemkul wrote: > >> >> >> On 2/5/14, 10:42 AM, gromacs query wrote: >> >&g

Re: [gmx-users] dna with charmm36

2014-03-03 Thread gromacs query
Hi Justin Yes I chose termini 5TER or 3TER correctly but it only works fine when D*5 replaced by D* and D*3 by D* otherwise it ends with errors as said previously. thanks JIom On Tue, Mar 4, 2014 at 1:00 AM, Justin Lemkul wrote: > > > On 3/3/14, 7:57 PM, gromacs query wrote: &g

Re: [gmx-users] No default U-B types charmm

2014-03-05 Thread gromacs query
Mar 3, 2014 at 9:07 PM, Justin Lemkul wrote: > > > On 3/3/14, 3:56 PM, gromacs query wrote: > >> What atom numbers are on line 184419? What atom types are they in the >>>> >>> [atoms] section for that [moleculetype]? >> >> These are waters only

[gmx-users] why not coordinates from cpt file

2014-03-11 Thread gromacs query
Hi all I have very simple query. While continuing simulations why we need to use *.gro (-c) with grompp as *.cpt (-t) has all the information (as checked with gmxcheck)? cpt file should suffice all the purposes. I tried using grompp providing -t *.cpt file but without -c *.gro file, it does not wo

[gmx-users] looping job: time not updating

2014-03-13 Thread gromacs query
Hi All Though its mentioned here how to restart jobs: http://www.gromacs.org/Documentation/How-tos/Extending_Simulations but I tried to run something like this (by scripting) but copying here how script is looping jobs: #stage 1 grompp -f npt.mdp -c conf_1ns.gro -t state_1ns.cpt -p topol.top -o

Re: [gmx-users] looping job: time not updating

2014-03-13 Thread gromacs query
Hi Jutsin Thanks, I was thinking that .cpt has all the information which I assumed would work like AMBER restart file to read time from .cpt (~ amber restart) file. regards Jiom -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_

[gmx-users] Bond length not finite error

2014-03-13 Thread gromacs query
Hi All I am simulating charmm36 popc membrane with ions and got this error (NPT with position restraints) before this I did minimization and NVT with posres. I could not find "Bond length not finite" error in archive. Step Time Lambda 00.0

Re: [gmx-users] Bond length not finite error

2014-03-13 Thread gromacs query
Hi Justin I reduced dt from 0.002 to 0.001, till now its working fine. I will post again if get some trouble. thanks Jiom -- Gromacs Users mailing list * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting! * Can't post? Read http://www.gro

[gmx-users] gen pairs missing from vmd top

2014-05-09 Thread gromacs query
Hi All I am using topo writegmxtop output.top in VMD and it writes fake sort of top file which needs further user based adjustments. The top file so generated has all [bonds] [angles] [dihedral] information but [pairs] are missing. I am using charmm36 which says gen-pair =yes in forcefield.itp.

[gmx-users] g_anaeig -extr option

2014-05-11 Thread gromacs query
Hi All In g_anaeig we have this option -extr: calculate the two extreme projections along a trajectory It provides me two pdb files which are extreme (I don't want to use interpolate -nframes option). I have simple doubt: say if I am doing this for Eigenvector 1 (EV1) then do these two pdbs simpl

[gmx-users] truncated octahedron box vector angles and number of waters

2014-06-13 Thread gromacs query
Dear All I am creating a truncated octahedron box and in vmd it appears as a rectangular box. In mailing list it has been suggested to use trjconv -ur compact option. But I have two issues, 1) I get box vector angle 70 90 70 with editconf and should be around 109 109 109 for a truncated octahedro

[gmx-users] -extend problem

2014-07-03 Thread gromacs query
Hi All I generated my npt.tpr with: tinit=0 dt=0.002 nsteps=125000 I am trying to extend simulation for 250 ps like this: tpbconv_mpi -s npt.tpr -o npt_new.tpr -extend 250 But when I used this npt_new.tpr simulation it ran for 500 ps. Have I missed something? Also one more thing why the

Re: [gmx-users] -extend problem

2014-07-03 Thread gromacs query
:07 PM, Justin Lemkul wrote: > > > On 7/3/14, 6:05 PM, gromacs query wrote: > >> Hi All >> >> I generated my npt.tpr with: >> >> >> tinit=0 >> >> dt=0.002 >> >> nsteps=125000 >> >> I am trying to extend simu

Re: [gmx-users] -extend problem

2014-07-03 Thread gromacs query
Hi forgive my ignorance but I could not find whether my cpt was read successfully or not as I dont see anything in log file as which files were taken as input. Or any keyword I should search for? On Thu, Jul 3, 2014 at 11:36 PM, Justin Lemkul wrote: > > > On 7/3/14, 6:34 PM, grom

Re: [gmx-users] -extend problem

2014-07-03 Thread gromacs query
2014 at 11:58 PM, Justin Lemkul wrote: > > > On 7/3/14, 6:51 PM, gromacs query wrote: > >> Hi >> >> forgive my ignorance but I could not find whether my cpt was read >> successfully or not as I dont see anything in log file as which files were >>

Re: [gmx-users] -extend problem

2014-07-03 Thread gromacs query
Hi Sorry I think I found the error. I was using wrong name of -cpi file. But its strange that mdrun does not complain even if it could not find the file and it still ran the job. On Fri, Jul 4, 2014 at 12:31 AM, gromacs query wrote: > Hi Justin > > Indeed this key word 'Read