-- 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!

Regarding my mention of "Wolfram Summer School 2020" in the original version of my notebook, I was not a student of this year's session, but I did apply for admission. I began working on my 1st notebook while I was (optimistically) waiting to hear if I would be accepted into WSS 2020. I wanted to get a head start, just in case. But there were so many great applicants, and not everyone could be squeezed in, and I didn't make the cut. But I totally understand. I've enjoyed watching the WSS 2020 student project presentations. Very exciting minds working at the cutting edge!

Nice and interesting! Looks like this notebook was published as private. Can you republish it with public permissions so we an take a look? Thank you.

This is great work! I look forward to reading through the code more thoroughly.

In the long run we definitely need a chemical reaction data structure that can operate on a list of molecules to make this easier on you. Until we get that, another approach you could take for chemical reactions is using the "ReplaceSubstructure" method:

alkyne = Molecule["1-octyne"]
primaryalkylhalide = Molecule["benzyl bromide"]

The trick here is to join the two molecules into one disconnected structure, then search for a single substructure matching the bromly group and terminal alkyne, and replace it with the joined version:

joinMolecules[list : {__Molecule}] := 
 Molecule@StringRiffle[MoleculeValue[list, "SMILES"], "."]

MoleculeModify[
 joinMolecules[{alkyne, primaryalkylhalide}], {"ReplaceSubstructure", 
  MoleculePattern["[CH2][Br].[CH1]#[C]"] -> "CC#C"}]

Hi Leonardo, interesting work! ..do you have SOCIAL (ex: like LinkedIn), so we can connect?