The stabilization of lifted flames involves a complex competition between small-scale mixing and propagation of flame kernels at the base of the flame. Here, the effect of species diffusivity on this stabilization process is explored. For this purpose, a large-eddy simulation (LES)/probability density function (PDF) methodology accounting for differential diffusion is developed. A dynamic model for scalar mixing time-scale is formulated to accurately describe the turbulence-driven mixing of scalars at the small-scales. Autoignition-stabilized H 2 /N 2 lifted turbulent flames in a vitiated coflow burner are chosen as test cases. Detailed chemical kinetics is used along with in situ adaptive tabulation (ISAT) to accelerate the computations. Four configuration test cases are considered: neglecting differential diffusion (Cases 1 and 2), considering differential diffusion (Case 3) and considering differential diffusion but neglecting molecular transport of H and H 2 (Case 4). A dynamic model for scalar mixing time-scale appearing in the mixing model is used in Cases 2, 3 and 4, which removes the need to specify the value of model constant. It is found that the lifted flames are very sensitive to differential diffusion. The predictions of major species mass fractions, temperatures and lift-off heights are in good agreement with the experimental data. Further, such agreement could be achieved without changing any boundary conditions. In particular, increased diffusion of species upstream allows ignition kernels to survive turbulent dissipation, causing stabilization to occur upstream. H 2 and H are identified as two critical species that play a dominant role in the differential diffusion effect. The model is also able to capture the change in lift-off height with change in coflow temperature. The results show that even in turbulent flames, stabilization of the flame zone can be sensitive to molecular diffusion. In LES computations where a significant part of this diffusion effect is resolved, the use of the enhanced LES/PDF/ISAT approach is seen as a promising approach for capturing the lift-off physics accurately.2 for in the solution to the one-dimensional flamelet equations. On the other hand, the spatial transport of this flame structure is through a conserved scalar (mixture fraction) and a reactive scalar representation. These scalars are necessarily transported using a single diffusivity. Consequently, the effect of differential diffusion at scales larger than the flame scales are neglected (or at best, inaccurately modeled). Blanquart et al. [31] have sought to address this issue by constructing mixture fraction equations that contain source terms to account for the change in species balance arising from differential diffusion. However, application of this approach to detailed chemistry based models is still being pursued. In the context of PDF methods, differential diffusion manifests in two different ways again. At the resolved scales, the PDF is altered by molecular diffusion through a second-order diffusion...