This paper investigates the multilayer Rayleigh–Taylor instability (RTI) using statistically stationary experiments conducted in a gas tunnel. Employing diagnostics such as particle image velocimetry (PIV) and planar laser induced fluorescence (PLIF), we make simultaneous velocity–density measurements to study how dynamics and mixing are linked in this variable density flow. Experiments are conducted in a newly built, blow-down three-layer gas tunnel facility. Mixing between three gas streams is studied, where the top and bottom streams are comprised of air, and the middle stream is an air–helium mixture. Shear is minimized between these streams by matching their inlet velocities. The four experimental conditions investigated here consist of two different density ratios (Atwood numbers 0.3 and 0.6), each investigated at two instability development times (or equivalently, two streamwise locations), and all experiments are with the same middle stream thickness of 3 cm. The growth of the middle layer is measured using laser-based planar Mie scattering visualization. The mixing width is found to grow linearly with time at late times. Various quantitative measures of molecular mixing indicate a very high degree of molecular mixing at late times in the multilayer RTI flow. The vertical turbulent mass flux
$a_y$
is calculated. In addition to mostly negative values of
$a_y$
, typical of buoyancy-dominated flows due to negative correlation between velocity and density fluctuations, positive regions are also observed in profiles of
$a_y$
due to entrainment and erosion at the lower edge of the mixing region. Global energy budgets are calculated for the multilayer RTI flow at late times and it is found that the majority of potential energy released has been dissipated due to viscous effects, and a large value of mixing efficiency (
$\sim$
60 %) is observed.