A linear analysis of electrostatic waves propagating parallel to the ambient field in a four component homogeneous, collisionless, magnetised plasma comprising fluid protons, fluid He++, electron beam, and suprathermal electrons following kappa distribution is presented. In the absence of electron beam streaming, numerical analysis of the dispersion relation shows six modes: two electron acoustic modes (modes 1 and 6), two fast ion acoustic modes (modes 2 and 5), and two slow ion acoustic modes (modes 3 and 4). The modes 1, 2 and 3 and modes 4, 5, and 6 have positive and negative phase speeds, respectively. With an increase in electron beam speed, the mode 6 gets affected the most and the phase speed turns positive from negative. The mode 6 thus starts to merge with modes 2 and 3 and generates the electron beam driven fast and slow ion acoustic waves unstable with a finite growth. The electron beam driven slow ion-acoustic waves occur at lower wavenumbers, whereas fast ion-acoustic waves occur at a large value of wavenumbers. The effect of various other parameters has also been studied. We have applied this analysis to the electrostatic waves observed in lunar wake during the first flyby of the ARTEMIS mission. The analysis shows that the low (high) frequency waves observed in the lunar wake could be the electron beam driven slow (fast) ion-acoustic modes.