Author Topic: outputing a file describing bond connectivity  (Read 3887 times)

0 Members and 1 Guest are viewing this topic.

Offline Ehsan

  • Regular QuantumATK user
  • **
  • Posts: 5
  • Country: ir
  • Reputation: 0
    • View Profile
outputing a file describing bond connectivity
« on: December 1, 2013, 16:00 »
hello all;
I have a xyz format of a structure. i can see the atoms and bonds of this structure but i need a text file that describe the bonds, i mean i need to know which atom is bonded to another. does atk have this ability?

Offline Anders Blom

  • QuantumATK Staff
  • Supreme QuantumATK Wizard
  • *****
  • Posts: 5574
  • Country: dk
  • Reputation: 96
    • View Profile
    • QuantumATK at Synopsys
Re: outputing a file describing bond connectivity
« Reply #1 on: December 1, 2013, 22:06 »
Not immediately, but it would be possible to do it as long as you can come up with a reasonable definition of a "bond". Bonds are not physical, it's something chemists like to use as a way to characterize the interactions between atoms. Sometimes they are mainly used for nice pictures, and sometimes they are crucial for computations since specific models, e.g. force fields, contain bond-specific parameters.

But under any circumstance the bond is a fictitious object, and what you mean by two atoms being bonded or not becomes a matter of definition.

In ATK you can compute the bond overlap population between two atoms; see http://quantumwise.com/documents/manuals/latest/ReferenceManual/index.html/ref.mullikenpopulation.html and search for "bond". If you can figure out a way to use that number to determine the bonding table, and very importantly the bond order, then this may be enough. Writing the output file is the easy part :)

Another way, which is just based on the overlaps of covalent radii, is also possible. It may be useful for converting files to specific formats, but one should note that the definition of bonds in this case is much more uncertain - which not least would be reflected in the fact that one could not assign any bond order to the pairs of atoms a priori. This method is not officially supported, in the sense that the API may change in the future, so it may not work in some newer version of ATK, but I could provide an example that works in 13.8 at least. The bonds printed if you don't change anything will be those that are used in VNL for visualizing the molecules, since the fuzz-factor (the extension of the covalent radii) is set to 1.1 which is what VNL uses.

Now, there is an added complication due to periodicity, since most file formats that I know which contain bonding information are for molecules. So the example attached is for a molecule (tetralin, to be specific).
« Last Edit: December 1, 2013, 22:08 by Anders Blom »

Offline Shinji Usui

  • QuantumATK Staff
  • Heavy QuantumATK user
  • *****
  • Posts: 44
  • Country: jp
  • Reputation: 2
    • View Profile
Re: outputing a file describing bond connectivity
« Reply #2 on: December 2, 2013, 03:08 »
if you just want to get a bond index table, you can obtain it by converting xyz file to mol file.

1. Output your structure in the XYZ format
2. Install Openbabel from http://openbabel.org/wiki/Main_Page
3. Launch your terminal or cmd.exe
4. Go to  the directory the XYZ file located
5. Command will be
> babel -ixyz C60.xyz -omol C60.mol

You will get the bond index table at the below part of the Molfile.
Like
  2  3  2  0  0  0  0
means Number 2 atom and Number 3 atom are bonding by a double bond.

As Anders mentioned these bond descriptions are not physical and just a conventional description.
« Last Edit: December 2, 2013, 03:17 by Shinji Usui »

Offline Ehsan

  • Regular QuantumATK user
  • **
  • Posts: 5
  • Country: ir
  • Reputation: 0
    • View Profile
Re: outputing a file describing bond connectivity
« Reply #3 on: December 2, 2013, 11:35 »
thank you Dr. Blom and Dr. Shinji, i got my answer, i need bond index table as you told me.