Large eddy simulation (LES) with the approximate deconvolution model is combined with Lagrangian particles simulation (LPS) for simulating turbulent reactive flows at high Schmidt numbers Sc. The LES is used to simulate velocity and nonreactive scalar while reactive scalars are simulated by the LPS using the mixing volume model for molecular diffusion. The LES-LPS is applied to turbulent scalar mixing layers with a second-order isothermal irreversible reaction at Sc5600. The mixing volume model is implemented with the IEM, Curl's, and modified Curl's mixing schemes. The mixing volume model provides a correct decay rate of nonreactive scalar variance at high Sc independently of the number of particles. The statistics in the LES-LPS with the IEM or modified Curl's mixing scheme agree well with the experiments for both moderately-fast and rapid reactions. However, the LPS with the Curl's mixing scheme overpredicts the effects of the rapid reaction. (C) 2016 American Institute of Chemical Engineers