Dears lammps users ,

I simulate a system containing 64 polymers, and I get a dump file containing angles, using compute dihedral/local command, but I don’t know how to calculate the dihedral angle distribution function (distribution probability) , (can I use ovito as we do to calculate the radial distribution function!! ) or there is a lammps command reserved to this!!

I am really stuck and don’t know what to do.

thanks in advance .

this sounds more like you want a histogram not a distribution function.

for a histogram you can use:

https://lammps.sandia.gov/doc/fix_ave_histo.html

thanks, but i want the probability distribution and get a graph like the figure below.

for my case i have the angle and i wanna get f(angle).

that is a mere plotting / postprocessing issue. that graph *is* a histogram, only each column is divided by the sum of all column values.