Multivariate Curve Resolution (MCR) was used to determine the phase purity of pharmaceutical co-crystals from mid infrared spectra. An in-silico coformer screening was used to choose one of ten potential coformers. This analysis used quantum chemistry simulation to predict which coformers are thermodynamically inclined to form cocrystals with the model drug, hydrochlorothiazide. The coformer chosen was nicotinamide. An experimental solvent screening by ultrasound assisted slurry co-crystallization was performed to evaluate the capacity of the method to determine phase purity. Afterwards, slurry and slow evaporation co-crystallizations were performed at 10, 25, and 40 °C using 7 solvent systems, and two levels of agitation for the evaporation co-crystallization (on and off). Mid infrared spectroscopy (MIRS) analysis of the products of these co-crystallizations was used to develop an MCR model to determine co-crystal phase purity. The MCR results were compared with a reference co-crystal. Experimental design (DoE) was used to investigate the effect of solvents, temperature, and agitation on the purity of co-crystals produced by slurry and evaporation co-crystallization. DoE revealed that evaporation co-crystallization with agitating at 65 rpm formed co-crystals with greater phase purity. The optimal temperature varied with the solvent used.