Creating fullerene molecules using BuckyballGraph

– you have earned Featured Contributor Badge
Your exceptional post has been selected for our editorial column Staff Picks http://wolfr.am/StaffPicks and Your Profile is now distinguished by a Featured Contributor Badge and is displayed on the Featured Contributor Board. Thank you!

@J.M. - this is a great application of FindIndependentEdgeSet.

It’s interesting to note that you can convert a BuckyballGraph to a molecule using the XYZ importer,

In[109]:= bucky = BuckyballGraph\[\];
xyzString = StringJoin[IntegerString[VertexCount[bucky]], "\n\n",
   ExportString[Prepend[#, "C"] & /@ GraphEmbedding[bucky], "Table"]
   ];
ImportString[xyzString, "XYZ"] // OutputForm

Out[111]//OutputForm= Molecule[<C60>]

In[112]:= Counts[%["BondType"]]

Out[112]= <|"Aromatic" -> 90|>