Cellular Automata (CA) models can represent dynamic systems which are discrete in space and time that reflects the effect of intrinsic parameters where individual events are considered to occur from randomness. A CA model of two agents' chemical kinetics has been optimized earlier using NSGA-II based on Evolutionary Algorithm (EA). But the stochastic nature of the CA model along with its high sensitivity on the model parameters requires extensive investigation using different optimization algorithms. For this purpose, in the current study, four more recently developed and popular optimization algorithms based on EA, called NSGA-IIr, NSGA-IIa, AbYSS and MOEA/D, have been considered for investigation based on various performance measuring parameters. The study also compares the performances of the algorithms for different computational efforts with an objective to minimize the required number of objective function evaluations. Simulation results and Friedman rank statistical test show NSGA-IIa and NSGA-IIr as the best choices to optimize the CA stochastic model across any number of objective function evaluations. Though the choice of optimization algorithm does not change with function evaluations, higher function evaluations improve the pseudo-pareto front for the CA optimization problem. Such results will facilitate the use of stochastic CA models to represent complex (bio)-chemical networks.