We consider the presence of oscillations in the primordial bispectrum, inspired by three different cosmological models; features in the primordial potential, resonant type non-Gaussianities and deviation from the standard Bunch Davies vacuum. In order to put constraints on their bispectra, a logical first step is to put these into factorized form which can be achieved via the recently proposed method of polynomial basis expansion on the tetrahedral domain. We investigate the viability of such an expansion for the oscillatory bispectra and find that one needs an increasing number of orthonormal mode functions to achieve significant correlation between the expansion and the original spectrum as a function of their frequency. To reduce the number of modes required, we propose a basis consisting of Fourier functions orthonormalized on the tetrahedral domain. We show that the use of Fourier mode functions instead of polynomial mode functions can lead to the necessary factorizability with the use of only 1/5 of the total number of modes required to reconstruct the bispectra with polynomial mode functions. Moreover, from an observational perspective, the expansion has unique signatures depending on the orientation of the oscillation due to a resonance effect between the mode functions and the original spectrum. This effect opens the possibility to extract informa- tion about both the frequency of the bispectrum as well as its shape while considering only a limited number of modes. The resonance effect is independent of the phase of the reconstructed bispectrum suggesting Fourier mode extraction could be an efficient way to detect oscillatory bispectra in the data.

17 pages, 12 figures. Matches published version