The epoch of reionization (EoR) 21-cm signal is expected to be highly
non-Gaussian in nature and this non-Gaussianity is also expected to evolve with
the progressing state of reionization. Therefore the signal will be correlated
between different Fourier modes (k). The power spectrum will not be able
capture this correlation in the signal. We use a higher-order estimator -- the
bispectrum -- to quantify this evolving non-Gaussianity. We study the
bispectrum using an ensemble of simulated 21-cm signal and with a large variety
of k triangles. We observe two competing sources driving the non-Gaussianity
in the signal: fluctuations in the neutral fraction (xHI) field and
fluctuations in the matter density field. We find that the non-Gaussian
contribution from these two sources vary, depending on the stage of
reionization and on which k modes are being studied. We show that the sign of
the bispectrum works as a unique marker to identify which among these two
components is driving the non-Gaussianity. We propose that the sign change in
the bispectrum, when plotted as a function of triangle configuration
cosθ and at a certain stage of the EoR can be used as a confirmative
test for the detection of the 21-cm signal. We also propose a new consolidated
way to visualize the signal evolution (with evolving xHI or
redshift), through the trajectories of the signal in a power spectrum and
equilateral bispectrum i.e. P(k)−B(k,k,k) space.Comment: 18 pages, 11 figures. Accepted for publication in MNRAS. Replaced to
match the accepted versio