Above the equality scale the galaxy bispectrum will be a key probe for measuring primordial non-Gaussianity which can help differentiate between different inflationary models and other theories of the early universe. On these scales a variety of relativistic effects come into play once the galaxy number-count fluctuation is projected onto our past lightcone. By decomposing the Fourier-space bispectrum into invariant multipoles about the observer's line of sight we examine in detail how the relativistic effects contribute to these. We show how to perform this decomposition analytically, which is significantly faster for subsequent computations. While all multipoles receive a contribution from the relativistic part, odd multipoles arising from the imaginary part of the bispectrum have no Newtonian contribution, making the odd multipoles a smoking gun for a relativistic signature in the bispectrum for single tracers. The dipole and the octopole are significant on equality scales and above where the Newtonian approximation breaks down. This breakdown is further signified by the fact that the even multipoles receive a significant correction on very large scales.

The Fourier-space galaxy bispectrum is complex, with the imaginary part arising from leading-order relativistic corrections, due to Doppler, gravitational redshift and related line-of-sight effects in redshift space. The detection of the imaginary part of the bispectrum is potentially a smoking gun signal of relativistic contributions. We investigate whether next-generation spectroscopic surveys could make such a detection. For a Stage IV spectroscopic $H\alpha$ survey similar to Euclid, we find that the cumulative signal to noise of this relativistic signature is $\mathcal{O}(10)$. Long-mode relativistic effects couple to short-mode Newtonian effects in the galaxy bispectrum, but not in the galaxy power spectrum. This is the basis for detectability of relativistic effects in the bispectrum of a single galaxy survey, whereas the power spectrum requires multiple galaxy surveys to detect the corresponding signal.

The bispectrum will play an important role in future galaxy surveys. On large scales it is a key probe for measuring primordial non-Gaussianity which can help differentiate between different inflationary models and other theories of the early universe. On these scales a variety of relativistic effects come into play once the galaxy number-count fluctuation is projected onto our past lightcone. We show for the first time that the leading relativistic correction from these distortions in the galaxy bispectrum generates a significant dipole, mainly from relativistic redshift space distortions. The amplitude of the dipole can be more than 10% of the monopole even on equality scales. Such a dipole is absent in the Newtonian approximation to the redshift space bispectrum, so it offers a clear signature of relativistic effects on cosmological scales in large scale structure.

The galaxy bispectrum on scales around and above the equality scale receives contributions from relativistic effects. Some of these arise from lightcone deformation effects, which come from local and line-of-sight integrated contributions. Here we calculate the local contributions from the generated vector and tensor background which is formed as scalar modes couple and enter the horizon. We show that these modes are sub-dominant when compared with other relativistic contributions.

The galaxy bispectrum is affected on equality scales and above by relativistic observational effects, at linear and nonlinear order. These lightcone effects include local contributions from Doppler and gravitational potential terms, as well as integrated contributions like lensing, together with all the couplings at nonlinear order. We recently presented the correction to the galaxy bispectrum from all local lightcone effects up to second order in perturbations, using a plane-parallel approximation. Here we update our previous result by including the effects from relativistic nonlinear dynamical evolution. We show that these dynamical effects make a significant contribution to the projection effects.

General relativistic imprints on the galaxy bispectrum arise from observational (or projection) effects. The lightcone projection effects include local contributions from Doppler and gravitational potential terms, as well as lensing and other integrated contributions. We recently presented for the first time, the correction to the galaxy bispectrum from all local lightcone projection effects up to second order in perturbations. Here we provide the details underlying this correction, together with further results and illustrations. For moderately squeezed shapes, the correction to the Newtonian prediction is ~30% on equality scales at z ~ 1. We generalise our recent results to include the contribution, up to second order, of magnification bias (which affects some of the local terms) and evolution bias.

Next-generation galaxy surveys will increasingly rely on the galaxy bispectrum to improve cosmological constraints, especially on primordial non-Gaussianity. A key theoretical requirement that remains to be developed is the analysis of general relativistic effects on the bispectrum, which arise from observing galaxies on the past lightcone, {as well as from relativistic corrections to the dynamics}. {As an initial step towards a fully relativistic analysis of the galaxy bispectrum, we compute for the first time the local relativistic lightcone effects on the bispectrum,} which come from Doppler and gravitational potential contributions. For the galaxy bispectrum, the problem is much more complex than for the power spectrum, since we need the lightcone corrections at second order. Mode-coupling contributions at second order mean that relativistic corrections can be non-negligible at smaller scales than in the case of the power spectrum. In a primordial Gaussian universe, we show that the local lightcone corrections for squeezed shapes at $z\sim1$ mean that the bispectrum can differ from the Newtonian prediction by $\gtrsim 10\%$ when the short modes are $k\lesssim (50\,{\rm Mpc})^{-1}$. These relativistic projection effects, if ignored in the analysis of observations, could be mistaken for primordial non-Gaussianity. For upcoming surveys which probe equality scales and beyond, {all relativistic lightcone effects and relativistic dynamical corrections should be included} for an accurate measurement of primordial non-Gaussianity.