5.26. Guessing unknown Topology information — MDAnalysis.topology.guessers¶
In general guess_atom_X returns the guessed value for a single value, while guess_Xs will work on an array of many atoms.
-
MDAnalysis.topology.guessers.get_atom_mass(element)[source]¶ Return the atomic mass in u for element.
Masses are looked up in
MDAnalysis.topology.tables.masses.Warning
Unknown masses are set to 0.0
Changed in version 0.20.0: Try uppercase atom type name as well
-
MDAnalysis.topology.guessers.guess_angles(bonds)[source]¶ Given a list of Bonds, find all angles that exist between atoms.
Works by assuming that if atoms 1 & 2 are bonded, and 2 & 3 are bonded, then (1,2,3) must be an angle.
Returns: List of tuples defining the angles. Suitable for use in u._topology Return type: list of tuples See also
-
MDAnalysis.topology.guessers.guess_atom_charge(atomname)[source]¶ Guess atom charge from the name.
Warning
Not implemented; simply returns 0.
-
MDAnalysis.topology.guessers.guess_atom_element(atomname)[source]¶ Guess the element of the atom from the name.
Looks in dict to see if element is found, otherwise it uses the first character in the atomname. The table comes from CHARMM and AMBER atom types, where the first character is not sufficient to determine the atom type. Some GROMOS ions have also been added.
See also
-
MDAnalysis.topology.guessers.guess_atom_mass(atomname)[source]¶ Guess a mass based on the atom name.
guess_atom_element()is used to determine the kind of atom.Warning
Anything not recognized is simply set to 0; if you rely on the masses you might want to double check.
-
MDAnalysis.topology.guessers.guess_atom_type(atomname)[source]¶ Guess atom type from the name.
At the moment, this function simply returns the element, as guessed by
guess_atom_element().
-
MDAnalysis.topology.guessers.guess_bonds(atoms, coords, box=None, **kwargs)[source]¶ Guess if bonds exist between two atoms based on their distance.
Bond between two atoms is created, if the two atoms are within
\[d < f \cdot (R_1 + R_2)\]of each other, where \(R_1\) and \(R_2\) are the VdW radii of the atoms and \(f\) is an ad-hoc fudge_factor. This is the same algorithm that VMD uses.
Parameters: - atoms (AtomGroup) – atoms for which bonds should be guessed
- coords (array) – coordinates of the atoms (i.e., AtomGroup.positions))
- fudge_factor (float, optional) – The factor by which atoms must overlap eachother to be considered a bond. Larger values will increase the number of bonds found. [0.55]
- vdwradii (dict, optional) – To supply custom vdwradii for atoms in the algorithm. Must be a dict
of format {type:radii}. The default table of van der Waals radii is
hard-coded as
MDAnalysis.topology.tables.vdwradii. Any user defined vdwradii passed as an argument will supercede the table values. [None] - lower_bound (float, optional) – The minimum bond length. All bonds found shorter than this length will be ignored. This is useful for parsing PDB with altloc records where atoms with altloc A and B maybe very close together and there should be no chemical bond between them. [0.1]
- box (array_like, optional) – Bonds are found using a distance search, if unit cell information is
given, periodic boundary conditions will be considered in the distance
search. [
None]
Returns: List of tuples suitable for use in Universe topology building.
Return type: Warning
No check is done after the bonds are guessed to see if Lewis structure is correct. This is wrong and will burn somebody.
Raises: ValueErrorif inputs are malformed or vdwradii data is missing.New in version 0.7.7.
Changed in version 0.9.0: Updated method internally to use more
numpy, should work faster. Should also use less memory, previously scaled as \(O(n^2)\). vdwradii argument now augments table list rather than replacing entirely.
-
MDAnalysis.topology.guessers.guess_dihedrals(angles)[source]¶ Given a list of Angles, find all dihedrals that exist between atoms.
Works by assuming that if (1,2,3) is an angle, and 3 & 4 are bonded, then (1,2,3,4) must be a dihedral.
Returns: - list of tuples – List of tuples defining the dihedrals. Suitable for use in u._topology
- .. versionadded 0.9.0
-
MDAnalysis.topology.guessers.guess_improper_dihedrals(angles)[source]¶ Given a list of Angles, find all improper dihedrals that exist between atoms.
Works by assuming that if (1,2,3) is an angle, and 2 & 4 are bonded, then (2, 1, 3, 4) must be an improper dihedral. ie the improper dihedral is the angle between the planes formed by (1, 2, 3) and (1, 3, 4)
Returns: - List of tuples defining the improper dihedrals.
- Suitable for use in u._topology
-
MDAnalysis.topology.guessers.guess_masses(atom_types)[source]¶ Guess the mass of many atoms based upon their type
Parameters: atom_types – Type of each atom Returns: atom_masses Return type: np.ndarray dtype float64