Let's understand this model briefly. If we set \(\Delta = 0\) , we see that this model is two copies of an untilted Weyl node in a magnetic field with Hamiltonian \(H_0 = \frac{\omega _B}{\sqrt{2}} (\sigma ^+ a + \sigma ^- a^{\dagger }) + \sigma ^z k_z\) , at different energies \(\pm b_0\) . The \(\Delta \)  term is a node hybridization term introduced by [1]}, and it gaps out the WPs at \(\Delta \ne 0\) . The \(\sigma \)  matrices represent the valley degrees of freedom, and the \(\tau \)  matrices denote the node d.o.f. Motivated by this, we first rewrite the Hamiltonian in Eqn.(REF ) in a more compact form:
\(H_{s,B}^n = s\Omega [\sigma ^+ a^n + \sigma ^- (a^{\dagger })^n] + sv\sigma ^z k_z + C_s k_z.\) 
