Skip to main content

Gnuplot to plot histogram

GNUPLOT can be used to generate histograms.
The script below allows one to plot variations in angle from 0° to 180°.
Each histogram box width correspond to 3° (covered here under bw option).
To execute this script; type this in the terminal:


gnuplot histogram.p

You should be in the folder where below script is saved in a file name "histogram.p" and a dat file (here, it is N---H-O.dat) should be present.

histogram.p
------------------
set term gif
set output "angle_N---H-O.gif"
set xrange [0.0:180.0]
set boxwidth 2.0 absolute
set style fill solid 1.0 border -1
set xtics 10
set title "molecule 3.0A Angle N---H-O"
set ylabel "occurance"
set xlabel "angle"
bw = 3   # substitute what you want
bin(x,width)=width*floor(x/width)
plot 'N---H-O.dat' using (bin($2,bw)):(1.0) smooth freq with boxes

------------------

The dat file is an input file with data in the format:
0    23.680344
1    73.766251
2    67.482254
3    52.086418
4    72.886803
5    54.009365
6    36.297497
7    35.184067
8    162.915512
9    52.081470
10    49.975559
11    46.514202
12    45.990093 ..............

Two columns where the second column refers to bond angles. First column here gives molecular dynamics steps.

On execution a gif file (here,angle_N---H-O.gif) will be created which looks something like this:


Labels:

Comments

Post a Comment

Popular posts from this blog

Material Studio - selective molecule movement

So, working with Material Studio is awesome. You get to handle bigger, complicated and even unimaginable system. But I am facing problems in very minute functioning. I recently faced two problems, but managed to resolve on my own. Somehow, tutorials doesn't help in resolving them and you end up doing some exercise on your own to get the solution. Among these, the problem of bringing two molecular structures in one window and then adjusting one molecule by moving it around, while other structure remains untouched, as it is. Actually, this should not be a problem to do. There is a simple way of doing this. Select the molecule to be moved and then there is a symbol (button) to select, which allows you to selectively move that molecule only. This is the same way it can be done in Material studio too.  But among all the dozens of symbols, the new person handling MS, kind of gets lost. So, what he does is check the tutorials, for any help. But the tutorials don't talk about suc...
Post a Comment
Read more

Using Packmol

Packmol is a software that allows you to generate a single combined PDB file for all the molecules you are considering. It is a free software and can be downloaded from this link: http://www.ime.unicamp.br/~martinez/packmol/ The installation guide is also provided on the page and its really easy to follow. http://www.ime.unicamp.br/~martinez/packmol/ You can also find some examples but I tried on my own. Anyway what you need is an input file which looks like: ------------------ tolerance 2.0 output trial.pdb add_amber_ter filetype pdb structure 2a.pdb   number 1   inside cube 0. 0. 0. 10. end structure structure 2d.pdb   number 1   inside cube 0. 0. 0. 10. end structure --------------------------------- IT includes some information like: how close the molecules considered can be.(tolerance) name of output pdb file "add_amber_ter" adds term "TER" after every molecule thus defining the start and end of every molecule. After mentioning the...
Post a Comment
Read more

COVID-19 Projections -till 22nd July

How would you call the COVID-19 situation in India for the last week. In contrast to what I hear around, I would say, the situation in India is improving. Lets see the reasons why I feel so: - Recovery rate has gone up to 63%. So, there is almost 2/3rd chance of survival for every person getting infected by COVID - The mortality rate has now reduced down to ~2.5%, way down from 3.2% in a 10 days duration.This is happening because of more and more testing. Today's report said, we did record 3.2 lakh testing in a day. This implies that the possibility of death because of COVID-19 is getting low. It may mean,  most of the people around can do fine with a little low panic level. Of course, precautions are necessary, if you don't want to fall in that very low 2.5% mortality set. Coming down to numbers. I am happy to see that India deviated significantly from my mid-July projection of 10.25 lakhs. Good to see the lowered down number. This shows the kind of ...
Post a Comment
Read more