In a previous study we have shown that chemical space of MOST candidates can be searched efficiently by genetic algorithms and we apply a similar approach here. The chemical space shown in Figure is encoded by a gene like that shown in Figure . Each of the four positions is represented by a label ("A", "B", and "C") followed by five integers ranging from 0 to 17. The labels represent the three types of ligands shown in Figure , b, and c, respectively, while the five integers represent the five different positions on the phenyl rings shown in Figure and c (if the label is "A", only the first digit is read). The GA search generates an initial population of 300 genes corresponding to random singly substituted NBD/QC systems (phenyl and phenylacetylene groups are also singly substituted and random mutations are applied if duplication occurs) after which a new generation of same size is created from crossovers between two parent genes. In the crossover, between one and three ligands from the first parent are randomly selected and combined with the remaining ligands from the second parent to produce a child. Furthermore, random mutations of the children occurs with a mutation rate of 25 %. These mutations allows one of the ligands to be modified by either changing the ligand type ("A", "B", or "C" from Figure ) or one of the attached "A" ligands. The different parents are selected with a probability that is proportional to their score (roulette selection). This score is computed as sum of three thresholded-linear functions ranging from 0 to 1 with thresholds of 0.6 MJ/kg, {375, 400, 450, 525} nm, and 250 kJ/mol for the energy density, absorption, and TBR barrier, respectively, predicted by the GCN models. To favour a redshift of the absorption spectrum, we multiplied the absorption score by two resulting in a maximum GA score of four. Hereafter, 300 of the worst scoring genes are discarded, which results in a new population of 300 genes. Each GA search runs for 300 generations and the final output is the 50 highest scoring molecules of the final population. We select four different absorption thresholds resulting in four runs of a thousand GA searches to also promote higher energy densities and TBR barriers.