  The Crab nebula is an astrophysical system that exhibits complex
morphological patterns at different observing frequencies. We carry out a
systematic investigation of the structural complexity of the nebula using
publicly available imaging data at radio and infrared frequencies. For the
analysis, we use the well-known multifractal detrended fluctuation analysis
(MFDFA) in two dimensions. We study the Crab nebula at various length scales to
search for local scaling behaviors. We find that both radio and infra-red data
exhibit long-range correlations, as expected from the underlying physics of the
supernova explosion and evolution. The correlations follow a power-law scaling
with length scales. The structural complexity is found to be multifractal in
nature, as evidenced by the dependence of the generalized Hurst exponent on the
order of the moments of the detrended fluctuation function. The nature and
strength of the multifractality vary with the observing frequency. By repeating
the analysis on shuffled data, we further probe the origin of the
multifractality in both the imaging data. For the radio data, we find that the
probability density function (PDF) is close to a Gaussian form, and hence the
multifractal behaviour is due to the differing nature of long-range
correlations of the large and small detrended fluctuation field values. In
contrast, we find that the infra-red data has a broad non-Gaussian PDF.
Consequently, its multifractal properties originate from the broad PDF as well
as the different correlations of large and small fluctuation values. Our
analysis thus provides a fresh perspective on the morphology of the Crab nebula
from a statistical physics viewpoint.
