HELLO
I want simulate gold-protein complex by Golp forcefield. but when I use pdb_gmx
I came across this error "There is a dangling bond at at least one of the
terminal ends and the force field does not provide terminal entries or files.
Edit a .n.tdb and/or .c.tdb
file."
PLEASE HELP ME
r
Dear all,
I am running a 9000 atom system with GBSA (Gromacs 4.5.5) in a
Blue Gene/Q cluster. I got the speed 1.002 ns/day with 8 cores.
However, in my own workstation with 8 cores the same system can reach
nearly 10 ns/day (Intel(R) Xeon(R) CPU E5620 @ 2.40GHz). Can anyone
tell me what's wr
Hi Justin,
I have the following Notes during NPT equilibration.
NOTE 1 [file pr_NPT.mdp]:
nstcomm < nstcalcenergy defeats the purpose of nstcalcenergy, setting
nstcomm to nstcalcenergy
NOTE 2 [file pr_NPT.mdp]:
leapfrog does not yet support Nose-Hoover chains, nhchainlength reset to 1
Wh
On 17/07/2012 4:32 PM, cuong nguyen wrote:
Dear Gmx-users,
I created a box size 4 4 2 and named layer.gro. Then "genconf" was
used to doulble this box:
genconf -f layer.gro -o 2layers.gro -nbox 1 1 2 -dist 0 0 8
However, the copied box has the same direction as the original box.
Could you please
On 17/07/2012 5:00 PM, DeChang Li wrote:
Dear all,
I am running a 9000 atom system with GBSA (Gromacs 4.5.5) in a
Blue Gene/Q cluster. I got the speed 1.002 ns/day with 8 cores.
However, in my own workstation with 8 cores the same system can reach
nearly 10 ns/day (Intel(R) Xeon(R) CPU E56
Dear cuong nguyen..
I think use following commands.
Try editconf -rotate
along these use -center < co-ordinate > if you want to place canter of
box at particular position
Try editconf -translate
On Tue, Jul 17, 2012 at 2:07 PM, Mark Abraham wrote:
> On 17/07/2012 4:32 PM, cuong nguyen wro
>--
>
>Message: 8
>Date: Tue, 17 Jul 2012 18:40:05 +1000
>From: Mark Abraham
>Subject: Re: [gmx-users] why Blue Gene/Q is so slow?
>To: Discussion list for GROMACS users
>Message-ID: <500524e5.9050...@anu.edu.au>
>Content-Type: text/plain; charset=ISO-8859-1; format=fl
On 7/17/12 4:05 AM, J Peterson wrote:
Hi Justin,
I have the following Notes during NPT equilibration.
NOTE 1 [file pr_NPT.mdp]:
nstcomm < nstcalcenergy defeats the purpose of nstcalcenergy, setting
nstcomm to nstcalcenergy
http://manual.gromacs.org/online/mdp_opt.html#out
nstcalcene
Dear All,
I rephrase my question (for the original see below).
Is there any (easy) way to give gromacs (mdrun -rerun) accurate coordinates for
a single point calculation (to get accurate energy and forces)?
The accuracy of .pdb and .gro formats is not sufficient and .g96
format reading seems to
Hello,
I've been trying some time now to run a QMMM calculation with gromacs and
orca on a system of 5000 atoms with the QM system being ~300 atoms.
Setting up the calculation with bOpt = yes and 1 minim step in the .mdp file
produces some problems.
Orca does the calculation without taking into
On 7/17/12 7:16 AM, Markus Kaukonen wrote:
Dear All,
I rephrase my question (for the original see below).
Is there any (easy) way to give gromacs (mdrun -rerun) accurate coordinates for
a single point calculation (to get accurate energy and forces)?
The accuracy of .pdb and .gro formats is n
> I have the following Notes during NPT equilibration.
>
> NOTE 1 [file pr_NPT.mdp]:
> nstcomm < nstcalcenergy defeats the purpose of nstcalcenergy, setting
> nstcomm to nstcalcenergy
>
>
> NOTE 2 [file pr_NPT.mdp]:
> leapfrog does not yet support Nose-Hoover chains, nhchainlength reset to 1
Dear gmxers,
I am trying to generate one polymer melt from one big enough box using
NPT MD in gmx. I find that the density of system varies very slow.
Could you please give me some hints about how to speed up this
process? Thanks a lot for any reply.
Yours sincerely,
Chaofu Wu.
--
gmx-users mai
Well,
the problem was solved at the end, please don't bother answering.
Following the advise of one of the Orca developers, Christoph, I changed
QMMMscheme to normal instead of ONIOM and disabled periodic boundary
conditions.
Now the files are extracted normally.
for further reference, one sh
Dear all,
I have a question about internal bond interactions within a molecule.
I understand that 1-4 interactions can be scaled with "nbfunc". But I wonder:
1. What is the default for 1-5 and larger interactions in Gromacs? Are the LJ
and charge interactions included (i.e., not scaled)? I'm s
A variabvle thermostat that increases from 180 K to 300 K over cycles of
100-1
Original-Nachricht
> Datum: Tue, 17 Jul 2012 20:40:12 +0800
> Von: Wu Chaofu
> An: gmx-users@gromacs.org
> Betreff: [gmx-users] How to speed up equilibrating the density of bulk
> system?
On 17/07/2012 7:06 PM, DeChang Li wrote:
--
Message: 8
Date: Tue, 17 Jul 2012 18:40:05 +1000
From: Mark Abraham
Subject: Re: [gmx-users] why Blue Gene/Q is so slow?
To: Discussion list for GROMACS users
Message-ID: <500524e5.9050...@anu.edu.au>
Content-Type: text/pl
>Convert the .g96 file to .trr format. Double precision may be necessary and
>will certainly give the greatest possible accuracy.
>-Justin
I try
>trjconv -f gromacs.g96 -o test.trr
and get
>Program trjconv, VERSION 4.5.5
>Source code file: /build/buildd/gromacs-4.5.5/src/gmxlib/trxio.c, line: 693
On 7/17/12 12:51 PM, Markus Kaukonen wrote:
Convert the .g96 file to .trr format. Double precision may be necessary and
will certainly give the greatest possible accuracy.
-Justin
I try
trjconv -f gromacs.g96 -o test.trr
and get
Program trjconv, VERSION 4.5.5
Source code file: /build/build
Hi,
You can scale the system with editconf to some density and energy
minimize, and equilibrate 1000 steps or so in NVT. Do this with small
scaling factors, repeating until tje density is correct. You might
want to run at a relatively high temperature.
Cheers,
Tsjerk
On Tue, Jul 17, 2012 at 2:4
Hi all,
I have a system that consists of two distinct reservoirs filled with
water+NaCl connected by a CNT. I want to measure the average number of
hydrogen bonds per water in each reservoir as a function of time, but
my current workflow for this is far too time-consuming to get any
significant re
On 7/17/12 4:37 PM, Alex Marshall wrote:
Hi all,
I have a system that consists of two distinct reservoirs filled with
water+NaCl connected by a CNT. I want to measure the average number of
hydrogen bonds per water in each reservoir as a function of time, but
my current workflow for this is far
Hi Justin,
I am just curious, do you have experience writing such a shell script to
iterate over the index files (one per frame)? I am just curious, because I
have been trying this, and I have found it very difficult to do using bash.
Do you use a different shell language?
Also, I wonder if an
On 7/17/12 4:57 PM, Andrew DeYoung wrote:
Hi Justin,
I am just curious, do you have experience writing such a shell script to
iterate over the index files (one per frame)? I am just curious, because I
have been trying this, and I have found it very difficult to do using bash.
Do you use a dif
Thanks Justin, this is great.
On Tue, Jul 17, 2012 at 5:07 PM, Justin Lemkul wrote:
>
>
> On 7/17/12 4:57 PM, Andrew DeYoung wrote:
>>
>> Hi Justin,
>>
>> I am just curious, do you have experience writing such a shell script to
>> iterate over the index files (one per frame)? I am just curious,
Hi,
I tried again simulating the same peptide with OPLS FF for 2000 ns time
duration. I used VMD for the secondary structure analysis and trajectory
visualization. This time I was able to get folding events but the residues
that are should be beta strands were helical for most of the time during
s
and this did not happen with EM?
Did you add new atoms/FF parameters to the existing C36 set?
On 2012-07-13 02:05:47AM -0700, Shima Arasteh wrote:
>
> Dear gmx users,
>
>
> My system is composed of a protein and water. I am working with CHARMM36 and
> the current version of Gromacs, 4.5.5.
>
I added a new residue by using defined atom types.
I did not get this error with EM. This error happened when I just entered the
command of #deffnm NVT.
Thanks in advance.
Sincerely,
Shima
From: Peter C. Lai
To: gmx-users@gromacs.org
Sent: Wednesday, July 1
Has anyone been recently using the coupling feature of mdrun for
parameterising? I used it many versions ago and have just been trying to use
it today with 4.5.5, however I can't get it to work as it did previously.
Currently you can't couple more than one set of charges (seems a bit weird th
29 matches
Mail list logo