The chemical complexity and diversity of an Athabasca asphaltene sample was described using a series of molecular representations. The molecular representations were created with a Monte Carlo construction method that represented molecules with a series of aromatic and aliphatic groups. After the groups were randomly sampled for a molecule, a connection algorithm linked them together to form molecules consisting of aromatic groups connected by aliphatic chains and thioethers. A sequential nonlinear optimization algorithm was used to select a small subset of molecules that were consistent with elemental, molecular weight, and NMR spectroscopy (both C-13 and H-1) data. To accurately represent the analytical data for the asphaltene sample, a minimum of five molecules was needed. On the basis of the results of the sequential optimization, at least 50 molecules in the starting population were required to produce an analytically consistent molecular representation.