Bayesian Optimization of Hubbard U’s for Investigating InGaN Superlattices