[gmx-users] g_mdmat

2009-02-18 Thread rob yang
Hello list, does g_mdmat ouput the standard deviation that goes along with the mean.xpm file? Thank you. -Rob _ ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org

RE: [gmx-users] g_hbond angle def

2009-02-06 Thread rob yang
om: sp...@xray.bmc.uu.se > To: gmx-users@gromacs.org > Subject: Re: [gmx-users] g_hbond angle def > > rob yang wrote: > > At the risk of exposing my density in triangle geometry here, why does > > 180-DHA = ADH? > > H > >/ \ > > D - A > >

RE: [gmx-users] g_hbond angle def

2009-02-06 Thread rob yang
At the risk of exposing my density in triangle geometry here, why does 180-DHA = ADH? H / \ D - A Thank you. -Rob > Obviously the DHA angle in a linear Hbond is > 180, so it should either be 180-DHA or ADH. The ADH angle is what's > implemented in gromacs > 3.2.1, however in

[gmx-users] g_hbond angle def

2009-02-05 Thread rob yang
Hi list, What's the geometric relation b/w the commonly used Donor-H-Acceptor angle and the gromacs implementation of Acceptor-Donor-H angle? I guess I don't quite follow why it's valid to define the ADH cutoff as 30 degrees, the same angle most commonly used to define DHA, which suggests ADH==

RE: [gmx-users] trjconv behaving differently in 3.3.3 vs previous versions

2008-03-07 Thread rob yang
Okay, my mistake. Here's what each option gave me (all with -ur compact): -pbc whole: broken protein, rectangular box -pbc mol: intact protein, rectangular box, part of protein sticking outside the box -pbc atom: intact protein, dodecahedron box > but -pbc mol -ur compact should give you a dod

RE: [gmx-users] .gro files generation for a certain time-step interval in a simulation

2008-03-07 Thread rob yang
trjconv > Date: Fri, 7 Mar 2008 15:54:11 + > From: [EMAIL PROTECTED] > To: gmx-users@gromacs.org > Subject: [gmx-users] .gro files generation for a certain time-step interval > in a simulation > > Hi All, > > I wonder if anyone knows what is the command to u

RE: [gmx-users] trjconv behaving differently in 3.3.3 vs previous versions

2008-03-07 Thread rob yang
Dave, I am using it for visualization for now. -pbc mol gave me the same output as -pbc whole, which was a cubic box with broken protein. -pbc atom gave the dodecahedron box and intact protein. But since the functionality of -pbc mol has changed, and may not be applicable to visualization purpo

RE: [gmx-users] trjconv behaving differently in 3.3.3 vs previous versions

2008-03-07 Thread rob yang
> Subject: Re: [gmx-users] trjconv behaving differently in 3.3.3 vs previous > versions > > rob yang wrote: >> Hi list, >> I have a system solvated in dodecahedron box. To visualize it, I used: >> trjconv -f system.gro -s system.gro -o visual.pdb -ur compact -pbc who

Re: [gmx-users] TIP5P and cg minimization

2008-03-06 Thread rob yang
TED]: Re: RE: [gmx-users] TIP5P and cg minimization Rob,I wonder whether you have got the flexible water model, TIP5P. I want to run simulations using the flexible TIP5P water model. Well, how does the results from the MD be affected if I use rigid water models in place of a flexible one?regards

[gmx-users] trjconv behaving differently in 3.3.3 vs previous versions

2008-03-06 Thread rob yang
Hi list, I have a system solvated in dodecahedron box. To visualize it, I used: trjconv -f system.gro -s system.gro -o visual.pdb -ur compact -pbc whole Version 3.3.3 gave me a cubic box with broken molecules while the previous versions gave me the right visualization on the same system.gro Does

RE: [gmx-users] TIP5P and cg minimization

2008-03-06 Thread rob yang
Thanks Mark for the pointer. The TIP5P.itp has [settles] only, whereas TIP4P.itp (as well as spc.itp, and TIP3P) has parameters for flexible water [FLEXIBLE] as well as rigid [settles]. This, by definition would restrict any systems using TIP5P to the SHAKE algorithm. I am wondering: 1) are t

[gmx-users] TIP5P and cg minimization

2008-03-05 Thread rob yang
Hi list, I have solvated my protein in TIP5P box that's provided by 3.3.3. However, when I try to run an energy minimization using cg (at the grompp stage), I get: ERROR: can not do Conjugate Gradients with constraints (218328) I don't have any problem using the same em.mdp file on TIP4P, so I

RE: [gmx-users] ffG43a1.rtp file error

2006-03-21 Thread rob yang
Hi Ken, #define gd_9 0.000 3.14 2 0.000 = phase shift. it takes on values of 0 or 180. 3.14 = force constant in KJ/mol 2 = multiplicity of the torsion term the torsion equation is: Sum of {(Force constant)/2*(1+cos(phase shift)*cos(multiplicity * your angle))} it can the

[gmx-users] oplsaa2005?

2006-03-14 Thread rob yang
Hi list, Just wondering if anybody have gromacs-formatted oplsaa2005 parameters? Thank you. Rob ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please don't post (un)subscribe requests to the lis

Re: [gmx-users] calculation without box

2006-03-10 Thread rob yang
7:26:43.0416 (UTC) FILETIME=[CCC10580:01C64467] It can be possibly done as you said... =) You may have your difficulty in this but not specified in the mail. Use nstep=1 to calculate the energy. Yang Ye rob yang wrote: Hi list, is there a way to calculate the energy of a given molecule wi

[gmx-users] calculation without box

2006-03-10 Thread rob yang
Hi list, is there a way to calculate the energy of a given molecule without having to put it in a periodic box?.. thanks Rob ___ gmx-users mailing listgmx-users@gromacs.org http://www.gromacs.org/mailman/listinfo/gmx-users Please don't post (un)

Re: [gmx-users] grompp single vs double precisions

2006-03-04 Thread rob yang
o/gmx-users>,<mailto:[EMAIL PROTECTED]> Errors-To: [EMAIL PROTECTED] Return-Path: [EMAIL PROTECTED] X-OriginalArrivalTime: 04 Mar 2006 05:16:52.0925 (UTC) FILETIME=[D91BDED0:01C63F4A] rob yang wrote: hi list, i am comparing the relative energy differences between different conformat

Re: [gmx-users] grompp single vs double precisions

2006-03-04 Thread rob yang
o/gmx-users>,<mailto:[EMAIL PROTECTED]> Errors-To: [EMAIL PROTECTED] Return-Path: [EMAIL PROTECTED] X-OriginalArrivalTime: 04 Mar 2006 05:16:52.0925 (UTC) FILETIME=[D91BDED0:01C63F4A] rob yang wrote: hi list, i am comparing the relative energy differences between different conformat

Re: [gmx-users] grompp single vs double precisions

2006-03-04 Thread rob yang
alTime: 04 Mar 2006 01:15:18.0051 (UTC) FILETIME=[197D9F30:01C63F29] rob yang wrote: hi list, i am comparing the relative energy differences between different conformations of a same molecule in vaccuum. i am using both gromacs and another software to compare the results. the forcefield for b

Re: [gmx-users] grompp single vs double precisions

2006-03-04 Thread rob yang
alTime: 04 Mar 2006 01:15:18.0051 (UTC) FILETIME=[197D9F30:01C63F29] rob yang wrote: hi list, i am comparing the relative energy differences between different conformations of a same molecule in vaccuum. i am using both gromacs and another software to compare the results. the forcefield for b

[gmx-users] grompp single vs double precisions

2006-03-03 Thread rob yang
hi list, i am comparing the relative energy differences between different conformations of a same molecule in vaccuum. i am using both gromacs and another software to compare the results. the forcefield for both softwares is oplsaa (2001), and i've manually verified that they are indeed the sam

[gmx-users] grompp single vs double precisions

2006-03-03 Thread rob yang
hi list, i am comparing the relative energy differences between different conformations of a same molecule in vaccuum. i am using both gromacs and another software to compare the results. the forcefield for both softwares is oplsaa (2001), and i've manually verified that they are indeed the sam