“…Since most of the 2D Dirac system is unlike the QED, which, in non-relativistic case, the momentum integral is ultraviolet convergent and doesn't need the cutoff, the most of the 2D Dirac systems as we discussed need the ultraviolet cutoff Λ (i.e., the Fourier transform of the Λ r ), and we estimate the ultraviolet curoff as the bandwidth of the silicene, which is about 4.8 eV, and such estimation is valid and enough for our computation in this article [21,22,20]. While for the effective mass m which is related to both the interlayer and intralayer hopping, we use the typical value of bilayer silicene which is m = 0.298m 0 [23,24,25,26,27], while for the bilayer graphene, our results are applicable after replace the effective mass as m = 0.033m 0 [28] or m = 0.029m 0 [29]. For the 3D Dirac or Weyl system, the longitudinal susceptibility as well as the current-current correlation function is needed to taking the chiral anomaly and the monopole formed by the Weyl nodes into account. )…”