While processing large set of .cif files I encountered an interesting issue regarding density calculations. By looking at the source of the Crystal.calculated_density command (look below), I found it interesting that it doesn't proceed with calculations if certain Z' criteria is not met. Am I missing something, as I thought that you only need Z (and not Z' values) for crystallographic density calculations.
Can anyone elaborate on this?
if self.z_value is None or self.z_prime is None or self.z_prime == 0.0: return None
I am currently working on a research project, where I need to process quite large amount of .cif files to remove both free solvent molecules and part of the bound solvent molecules (defined) to further use the "modified" crystallographic data for set of theoretical calculations. I am guessing that such procedure requires modeling of a "grid" and then initially sorting components by their weight (lighter matrix should be related to the non-bound solvent molecules). However, this has turned out to be a very daunting task for me to achieve.
Maybe someone has already done something similar or at least can provide some ideas to deal with the task at hand?
Any help would be much appreciated!