The Bayesian optimization algorithm and the PI stopping criterion was developed and discussed in previous work . In brief, it is composed of Gaussian process (GP) model and acquisition functions (AF). GP was a mixed kernel (Eq. S3), combining the Maté rn52 kernel (Eq. S1) with the categorical kernel (Eq. S2), to handle the reaction's design space, which includes both continuous and categorical