A surrogate-based optimization routine is coupled with a Multi-Zone Model (MZM) for optimizing the combustion efficiency of a Reactivity Controlled Compression Ignition engine. In the first step, some model parameters are calibrated against experimental data. In the second step, initial reactive fuel distribution is determined by optimizing the engine’s combustion efficiency. The results suggest that a highly stratified initial fuel distribution should be aimed for to maximize combustion efficiency.