Analysis of large eddy simulation data reveals that Langmuir circulation (LC) induces a significant enhancement of the mixed layer deepening, only if the mixed layer depth (MLD) h is shallow and the buoyancy jump across it Delta B is small, when simulations are initiated by applying the wind stress to a motionless mixed layer with stratification. The difference in the entrainment rate between the cases with and without LC decreases with h Delta B/upsilon(2)(L), where upsilon(L) is the velocity scale of LC. The ratio of the mixing length scale / between the cases with and without LC is close to 1 for larger Rt [=(Nl(0)/q)(2); Rt > similar to 1], but it increases to above 10 with the decrease of Rt, where N is the Brunt-Vaisala frequency and q and l(0) are the velocity and length scales of turbulence in the homogeneous layer. It is also found that, in the presence of LC, the effect of stratification on vertical mixing should be parameterized in terms of Rt instead of Ri (=(N/S)(2)), because velocity shear S is no longer a dominant source of turbulence. The parameterization is provided by l/l(0) = (1 + alpha Rt)(-1/2) with alpha similar to 50, regardless of the presence of LC. However, LC makes l(0) much larger than conventionally used for the boundary layer.