We first solve the model by performing calculations in an armchair-edged sample with the width \(W=40\sqrt{3} a\) . Two terminals L and R are
assumed to be semi-infinite. A relatively-short scattering region with the length \(d=10a\) 
is fixed, and the main electrically-tunable parameters are the bias voltage \(U\)  and the Fermi energy \(E_{\rm F}\)  (see details in Appendix B) and Ref. [GavLaz]). We keep \(\vert U \vert \)  the same for P and AP configurations and fix the parameters \(t=2.7\)  eV and \(\gamma =0.39\)  eV from bare GBL [1]} and \(M_0=27\)  meV, \(\Delta =13.5\)  meV and \(\lambda _{\rm R}=0.9\)  meV.
