[gmx-users] Dihedral transition times for a non-protein structure
giu.tesei at gmail.com
Tue Oct 4 16:31:17 CEST 2011
I'm sorry, actually I forgot to add the flag in the email but I added it
when I started the program.
g_angle -f traj.trr -s topol.tpr -n dihe.ndx -oh histo.xvg -type dihedral
I get the following output:
Group 0 (Phi=0.0_3_5.92) has 232 elements
There is one group in the index
Reading file diheA3.tpr, VERSION 3.3.3 (single precision)
trn version: GMX_trn_file (single precision)
Last frame 10000 time 20000.000
Found points in the range from 0 to 359 (max 360)
< angle > = 0.374572
< angle^2 > = 262.338
Std. Dev. = 16.1925
Order parameter S^2 = 0.344698
but it generates only a angdist.xvg file.
I'm wondering why it's not generating the histo.xvg file.
Thanks for your help!
2011/10/4 Justin A. Lemkul <jalemkul at vt.edu>
> Giulio Tesei wrote:
>> I'm performing the analysis of a 30 residues PVA polymer chain in SPC
>> I'm trying to find the dihedral transition times of each residue (possibly
>> as an histogram 'number of transitions vs. residue number').
>> I've already tried:
>> 1) g_angle -f traj.trr -s topol.tpr -n angle.ndx -oh histo.xvg
>> but it gives no output besides angdist.xvg
> Add "-type dihedral" to this command, otherwise the analysis is for normal
> 2) g_dih -f traj.trr -s topol.tpr -o dihe.out
>> but it gives no output besides an empty dihe.out
>> 3) g_chi (it says there are no dihedral angles in the structure. I guess
>> that this program is for proteins only and it doesn't recognize PVA backbone
>> I've tried these programs on both gromacs 4.5.3 and gromacs 3.3.3.
>> Thanks for any incoming help!
> Justin A. Lemkul
> Ph.D. Candidate
> ICTAS Doctoral Scholar
> MILES-IGERT Trainee
> Department of Biochemistry
> Virginia Tech
> Blacksburg, VA
> jalemkul[at]vt.edu | (540) 231-9080
> gmx-users mailing list gmx-users at gromacs.org
> Please search the archive at http://www.gromacs.org/**
> Support/Mailing_Lists/Search<http://www.gromacs.org/Support/Mailing_Lists/Search>before posting!
> Please don't post (un)subscribe requests to the list. Use the www interface
> or send it to gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/**Support/Mailing_Lists<http://www.gromacs.org/Support/Mailing_Lists>
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the gromacs.org_gmx-users