Flows in low Prandtl number liquid pools are relevant for various technical applications, and have so far only been investigated for the case of pure fluids, i.e. with a constant, negative surface tension temperature coefficient ∂γ/∂T . Real-world fluids containing surfactants have a temperature dependent ∂γ/∂T > 0, which may change sign to ∂γ/∂T < 0 at a critical temperature T c . Where thermocapillary forces are the main driving force, this can have a tremendous effect on the resulting flow patterns and the associated heat transfer.Here we investigate the stability of such flows for five Marangoni numbers in the range of 2.1 × 10 6 ≤ M a ≤ 3.4 × 10 7 using dynamic large eddy simulations (LES), which we validate against a high resolution direct numerical simulation (DNS). We find that the five cases span all flow regimes, i.e. stable laminar flow at M a ≤ 2.1 × 10 6 , transitional flow with rotational instabilities at M a = 2.8 × 10 6 and M a = 4.6 × 10 6 and turbulent flow at M a = 1.8 × 10 7 and M a = 3.4 × 10 7 .