Create a PDB file representing the cone models from the centre of mass (CoM) analysis.
cone_type: The type of cone model to represent.
scale: Value for scaling the pivot-CoM distance which the size of the cone defaults to.
file: The name of the PDB file.
dir: The directory where the file is located.
force: A flag which, if set to True, will overwrite the any pre-existing file.
WARNING: This analysis is now defunct!
This creates a PDB file containing an artificial geometric structure to represent the various cone models. These models include:
|`diff in cone'|
|`diff on cone'|
The model can be selected by setting the cone type to one of these values. The cone is represented as an isotropic cone with its axis parallel to the average pivot-CoM vector, the vertex placed at the pivot point of the domain motions, and the length of the edge of the cone equal to the pivot-CoM distance multiplied by the scaling factor. The resultant PDB file can subsequently read into any molecular viewer.
There are four different types of residue within the PDB. The pivot point is represented as as a single carbon atom of the residue `PIV'. The cone consists of numerous H atoms of the residue `CON'. The average pivot-CoM vector is presented as the residue `AVE' with one carbon atom positioned at the pivot and the other at the head of the vector (after scaling by the scaling factor). Finally, if Monte Carlo have been performed, there will be multiple `MCC' residues representing the cone for each simulation, and multiple `MCA' residues representing the varying average pivot-CoM vector for each simulation.
To create the diffusion in a cone PDB representation, a uniform distribution of vectors on a sphere is generated using spherical coordinates with the polar angle defined from the average pivot-CoM vector. By incrementing the polar angle using an arccos distribution, a radial array of vectors representing latitude are created while incrementing the azimuthal angle evenly creates the longitudinal vectors. These are all placed into the PDB file as H atoms and are all connected using PDB CONECT records. Each H atom is connected to its two neighbours on the both the longitude and latitude. This creates a geometric PDB object with longitudinal and latitudinal lines representing the filled cone.