Follow the procedure in [1]} or [2]} and use the results (REF ), (REF ) of Appendix in our paper, we know that for each \(\xi _j\)
\(m^{mod,j}(z)=I+\frac{\left(2t\epsilon _j\theta ^{\prime \prime }(\xi _j)\right)^{-\frac{1}{2}}}{z-\xi _j}\left(\begin{array}{cc}0 & \epsilon _j i\beta ^{\xi _j}_{12}\\-\epsilon _j i\beta ^{\xi _j}_{21} & 0\end{array}\right)+\mathcal {O}(t^{-1}).\) 
