I am reading this article: Edge superconducting correlation in the attractive-U Kane-Mele-Hubbard model.
Considering just the first part of the article, where a negative-U Hubbard model with the intrinsic Spin-Orbital term is considered on the infinite single layer graphene. In deriving the gap equation and the so-called number equation. The author assumed that the superconductivity order parameter in sub-lattice A and sub-lattice B is equal. which is $$\Delta_A=\frac{U}{N}\sum_k\langle a_{-k\downarrow}a_{k\uparrow}\rangle=\Delta_B=\frac{U}{N}\sum_k\langle b_{-k\downarrow}b_{k\uparrow}\rangle=\Delta$$
However, according to my calculation, using the mean-field Hamiltonian with $\Delta_A=\Delta_B=\Delta$ in the article, I found that $\Delta_A \neq \Delta_B$. The gap-equation in the paper can only obtained when I let $\Delta=\frac{\Delta_A+\Delta_B}{2}$; similar to number equation, I have to average $n_a$ and $n_b$.
Because this is a rather lengthy calculation, I am not 100 percent sure that my calculation is right. I confirm myself that the Spin-Orbital term is not invariant under exchange of A and B, so there is no A,B sub-lattice symmetry.
So my question is: am I right? If so, is it proper in the article to let the order parameter on the two sub-lattice equal?