The density field of underexpanded supersonic free jets issuing from a choked circular nozzle was measured using a Rayleigh scattering-based technique. This reliable and nonintrusive technique is particularly suitable for high-speed flows and is fundamentally superior to the intrusive probes and particle-based techniques such as laser Doppler velocimetry. A continuous wave laser and photon counting electronics were employed for time and phase-averaged density measurements. The use of dust-free air for the entrained flow allowed measurements in the shear layer region. The free jets were produced in the plenum to ambient pressure ratio range of 1.88–5.75, which corresponded to a fully expanded Mach number range of 0.99⩽Mj⩽1.8. A comparative study of schlieren photographs and time-averaged density data provided insight into the shock-cell structures. The radial profiles obtained at various axial stations covering a downstream distance of 10 jet diameters show the development of the jet shear layer and the decay of the shock–cells. The supersonic free jets produced screech sound. A phase-averaged photon counting technique, using the screech tone as the trigger source, was used to measure the unsteady density variation. The phase-averaged density data show the evolution of the large-scale turbulent vortices that are found to be modulated periodically along the flow direction. A comparison with previously obtained data showing near-field pressure fluctuation and convective speed of the organized vortices reveals many interesting dynamics. All quantities show regular spatial modulation. The locations of local maxima in density fluctuations are found to coincide with the high convective speed and the antinode points in the near-field pressure fluctuation. Interestingly, the periodicity of modulation is found to be somewhat different from the shock spacing. Instead it shows that the standing wave system, known to exist in the near-field pressure fluctuation, extends into the jet shear layer. The standing wave is formed between the downstream moving Kelvin–Helmholtz instability waves and the upstream propagating part of sound waves. A detailed field measurement of the unsteady density fluctuation was conducted for the Mj=1.19 and 1.42 jets for which the near-field pressure fluctuation data were obtained previously. The phase-matched, combined plots of the density fluctuation present inside the jet flow, and the pressure fluctuation present just outside the jet boundary provide a charming insight into the shock–vortex interaction leading to the sound wave generation.