“…Also, at high frequency, the complete recovery of the gain or absorption saturation is not allowed, which means that the carrier density modulation is not as deep as the optical modulation. In this case, a frequency domain simulation using the multiple modes as the basis is more attractive than the time domain simulation that can include the spectral dependencies, at best, by gain-spectrum filtering and the group velocity dispersion [7], [9]. The frequency domain method has been used either for qualitative descriptions [2], [12] or in a semi-quantitative way including only the nearest neighbor mode-coupling [13].…”