We have used UV and visible spectrophotometry to study self-association of aromatic riboflavin molecules (RFN, vitamin B 2 ,7,isoalloxazine) in aqueous solution (pH 6.86) at T = 298 K, using a dimer model. We have determined the equilibrium dimerization constant for riboflavin, K dB = 125 ± 40 M −1 . We have studied heteroassociation in the system of molecules of 7,8-dimethyl-10-ribitylisoalloxazine with 1,3,7-trimethylxanthine (caffeine) and sodium salicylate (NAS) in aqueous solution (pH 6.86; T = 298 K). We have determined the heteroassociation constants for RFN-NAS and RFN-caffeine molecules in the absence and in the presence of urea in solutions using a modified Benesi-Hildebrand equation: 25 ± 4, 17 ± 3, and 74 ± 11, 53 ± 7 M −1 respectively. We have determined the dimerization constants for NAS (2.7 ± 0.5 M −1 ) and caffeine (17.0 ± 1.5 M −1 ). We conclude that heteroassociation of the aromatic molecules leads to a lower effective riboflavin concentration in solution, and the presence of urea in mixed solutions leads to an decrease in the complexation constants for the RFN-NAS and RFN-caffeine systems.