A procedure that combines a one‐point probability density function (pdf) and two‐point correlation functions (TPCF’s) is introduced to simulate the multispecies chemical reactions carried by a homogeneous turbulence field. In this procedure, the pdf equation and the evolution of TPCF’s are treated in parallel using time splitting. The one‐point pdf equation requires closure of only the molecular diffusion terms. At any instant the time scale for the linear mean square estimate closure is obtained from the parallel solutions of the equations for TPCF’s, using a physical‐space EDQNM (eddy damped quasinormal Markovian) closure for the advection terms and a similarity assumption for the reaction terms. With this procedure, a study of chemical reactions with four or five species is possible. The model’s predictions are in good agreement with DNS data and a fully two‐point pdf description of simpler kinetic schemes. Results for a range of Damko¨hler numbers, Schmidt numbers, and Reynolds numbers that are beyond the reach of DNS are also shown.