Hirota bilinear formalism and Supersymmetry

Extending the gauge-invariance principle for $\tau$ functions of the standard bilinear formalism to the supersymmetric case, we define $\mathcal{N}=1$ supersymmetric Hirota bilinear operators. Using them we bilinearize supersymmetric nonlinear evolution equations. The super-soliton solutions are discussed. As a quite strange paradox it is shown that the Lax integrable supersymmetric KdV of Manin-Radul-Mathieu equation does not possesses N super-soliton solution for $N\ge 3$ for arbitrary parameters. Only for a particular choice of them the N super-soliton solution exists.