We compute numerically the second- and third-order nonlinear magnetic susceptibilities of an Ising ladder model in context two-dimensional coherent spectroscopy by using infinite time-evolving block decimation method. The couples a quantum chain to bath polarized spins, thereby effecting decay spinon excitations. show that its susceptibility contains robust echo signal weak-coupling regime, whi...