We introduce a general optimization algorithm based on an interpolation of property values on a hypercube. Each vertex of the hypercube represents a molecule, while the interior of the interpolation represents a virtual superposition (alchemical mutation) of molecules. The resultant algorithm is similar to branch-and-bound/tree-search methods. We apply the algorithm to the optimization of the first electronic hyperpolarizability for several tolane libraries. The search includes structural and conformational information. Geometries ...