An increasing number of numerical simulations and experiments describing the turbulent spectrum of Rayleigh-Taylor (RT) mixing layers came to light over the past few years. Results reported in recent studies allow to rule out a turbulenceà la Kolmogorov as a mechanism acting on a self similar RT turbulent mixing layer. A different mechanism is presented, which complies with both numerical and experimental results and relates RT flow to other buoyant flows.PACS numbers: 52.35. Py, 47.20.Bp, 47.27.eb, 47.27.te A Rayleigh-Taylor (RT) instability [1, 2] occurs whenever a light fluid, ρ 1 , pushes a heavy fluid, ρ 2 , or similarly, when a heavy fluid on top of a lighter fluid is subject to a gravitational acceleration. The understanding of such instability in the developed turbulent regime is of primary interest to many fields of physics and technology since it is an important cause of mixing between two fluids of different density. In astrophysics for instance, it is responsible for the outward acceleration of a thermonuclear flame in type Ia supernovae [3], but it also plays an important role in shaping the interstellar medium so that new stars can be born [4]. The technology of confinement fusion also relies on a good understanding of RT mixing [5] and ways to reduce it [6]. The RT flow of two incompressible fluids in the low Atwood limit, A = (ρ 2 − ρ 1 )/(ρ 2 + ρ 1 ) ≪ 1 (Boussinesq approximation), is governed by a concentration equation (1), the Navier Stokes equation supplemented with a buoyant source term (2) and the incompressibility constraint (3)where g is a stationary and uniform gravitational acceleration vector field (i.e. planar symmetry is assumed). The coefficient κ is the molecular diffusion coefficient and ν is the kinematic viscosity of the mixture. They are both supposed constant. Without loss of generality, g is parallel to the z-axis. This is why, for any generic physical value Φ, the average (so defined numerically and experimentally) will be Φ (z) = 1 S S dx dy Φ throughout this paper. The fluctuating part will be denoted with a prime and defined as Φ ′ = Φ − Φ . With the increasing capacity of super computers, many simulations of RT flows in the developed turbulent regime have been performed, which describe the velocity spectrum E(k) [7,8,9], defined in such a way that u ′ 2 = dk E(k), or the concentration spectrum E c (k) [10,11,12,13], defined as c ′ 2 = dk E c (k), or both [14,15]. In the same time, although fewer in number, experimental investigations of E(k) [16] and E c (k) [10,17,18] have been carried out. A quick inspection of these results shows that no consensus arises concerning the shape of these spectrum. From a theoretical point of view, the situation is not more satisfactory. In [19] it is claimed that the Kolmogorov-Obukhov scheme, E(k) ∼ k −5/3 , holds in 3D RT mixing given that the effect of buoyancy on turbulence, although fundamental at the largest scale, becomes irrelevant at smaller scales. In [20] the particular RT time scale 1/ √ kgA at wave number k has been postu...