Skip to main content

Automorphism groups and fragment graphs

Structure generation involves not just graph theory, but group theory. Or, I should say, it does in some of the papers I have read. For example, in this paper by J.L.Faulon, there is the sentence:
"The two main steps are to compute the orbits of the automorphism group of G and to saturate all the atoms of a chosen orbit
which may well be incomprehensible to many readers, except if the reader is a mathematician.

I am no mathematician, but thanks to some books on groups, I now understand both what an automorphism group is and what an orbit is. On the other hand, I also believe that this definition of how the algorithm works is overly complex. A more simple term might just be "fragment sets" - as it is fairly clear, if not mathematically exact. So, for the fragment graph [CH3, CH3, CH2, CH2, CH, CH] the fragment set is [CH3, CH2, CH].

Anyway, here is a short analysis of the automorphism group of the fragment graph [CH2, CH2]. This first image shows the tiny group of permutations that swaps the two fragments:


The notation is taken from an excellent book called "Visual Group Theory" that is also associated with some software called group explorer on sourceforge. It might be quite general, I suppose (and I hope I'm using it right), but it shows the permutation that swaps the fragments as a circled s. This is an automorphism with respect to the edges - in other words, after the swap, there are still bonds between [1, 2], [2, 3], [4, 5], and [5-6].

Another part of the automorphism group is a 'flip' like:

which is a little more complex, but shows how 'flipping' each fragment separately combines to form four possible permutations. If this does not seem particularly tricky, consider what happens if you take the direct product of these two groups:

Assuming I have done it right, this should show most (all?) of the automorphisms of the fragment graph. It does look pretty cool, but I don't think that it gets me any closer to implementing the cursed algorithm :)

Comments

Anonymous said…
I really like your blog and i really appreciate the excellent quality content you are posting here for free for your online readers. thanks peace claudia.
Anonymous said…
Give never the wolf the wether to keep.

Popular posts from this blog

Adamantane, Diamantane, Twistane

After cubane, the thought occurred to look at other regular hydrocarbons. If only there was some sort of classification of chemicals that I could use look up similar structures. Oh wate, there is . Anyway, adamantane is not as regular as cubane, but it is highly symmetrical, looking like three cyclohexanes fused together. The vertices fall into two different types when colored by signature: The carbons with three carbon neighbours (degree-3, in the simple graph) have signature (a) and the degree-2 carbons have signature (b). Atoms of one type are only connected to atoms of another - the graph is bipartite . Adamantane connects together to form diamondoids (or, rather, this class have adamantane as a repeating subunit). One such is diamantane , which is no longer bipartite when colored by signature: It has three classes of vertex in the simple graph (a and b), as the set with degree-3 has been split in two. The tree for signature (c) is not shown. The graph is still bipartite accordin

1,2-dichlorocyclopropane and a spiran

As I am reading a book called "Symmetry in Chemistry" (H. H. Jaffé and M. Orchin) I thought I would try out a couple of examples that they use. One is 1,2-dichlorocylopropane : which is, apparently, dissymmetric because it has a symmetry element (a C2 axis) but is optically active. Incidentally, wedges can look horrible in small structures - this is why: The box around the hydrogen is shaded in grey, to show the effect of overlap. A possible fix might be to shorten the wedge, but sadly this would require working out the bounds of the text when calculating the wedge, which has to be done at render time. Oh well. Another interesting example is this 'spiran', which I can't find on ChEBI or ChemSpider: Image again courtesy of JChempaint . I guess the problem marker (the red line) on the N suggests that it is not a real compound? In any case, some simple code to determine potential chiral centres (using signatures) finds 2 in the cyclopropane structure, and 4 in the

General Graph Layout : Putting the Parts Together

An essential tool for graph generation is surely the ability to draw graphs. There are, of course, many methods for doing so along with many implementations of them. This post describes one more (or perhaps an existing method - I haven't checked). Firstly, lets divide a graph up into two parts; a) the blocks, also known as ' biconnected components ', and b) trees connecting those blocks. This is illustrated in the following set of examples on 6 vertices: Trees are circled in green, and blocks in red; the vertices in the overlap between two circles are articulation points. Since all trees are planar, a graph need only have planar blocks to be planar overall. The layout then just needs to do a tree layout  on the tree bits and some other layout on the embedding of the blocks. One slight wrinkle is shown by the last example in the image above. There are three parts - two blocks and a tree - just like the one to its left, but sharing a single articulation point. I had