Turbulence mixing models of different degree of complexity are investigated for Rayleigh-Taylor mixing flows with reference to high-resolution implicit large eddy simulations. The models considered, in order of increasing complexity, comprise the (i) two-equation K-L, (ii) three-equation K-L-a, (iii) four-equation K-L-a-b, and (iv) Besnard-Harlow-Rauenzahn (BHR-2). The above models are implemented in the same numerical framework to minimize the computational uncertainty. The impact of the various approximations represented by the different models is investigated for canonical one-dimensional (1D) Rayleigh-Taylor mixing and for the more complex (2D on average) case of the tilted-rig experiment, aiming to understand the balance between accuracy and complexity. The results provide guidance on the relative merits of various turbulence models over a variety of conditions.