“…By now, many efficient approaches have been available for finding soliton solutions, some of which include inverse scattering transformation [1,2], Darboux transformation [3][4][5][6], Bäcklund transformation [7], Riemann-Hilbert method [8], and Hirota's bilinear method [9][10][11][12]. In recent years, there has been an increasing interest in exploring abundant multiple soliton solutions of NLEEs via the Riemann-Hilbert method, including the coupled derivative Schrödinger equation [13], the Kundu-Eckhaus equation [14], and others [15][16][17][18][19][20][21].…”