Hi,
I have been using freud for the past few months to compute the radial distribution function. So far I have had an equal number of particles for every data frame. I have a 3d array of shape 2001, 450, 3 (2001 data frames, 450 particles, and x,y,z coordinates) which I fed to freud and computed the rdf and it worked.
Now, I have variable particle numbers for every data frame. So I tried to make an empty 3d zero array where for every data frame, I replaced the zeros with the particle coordinates. But the rdf from this method gives a misleading picture in the sense that I have values along the y-axis for x value between 0 and 1 which indicates that my particles are overlapping which is not the case.
Now, if I find the frame that has the lowest number of particles and then truncate the whole data (eliminating the zeros along with some particle coordinates), then the rdf gives the correct picture. But with this method, I lose many particle coordinates.
So, does freud have any other way to compute the rdf of particles where the number of particles changes with every data frame?
Thank you for your time and consideration.
Regards,
Jenis Thongam.