本論文では、3次元トポロジカル絶縁体の2次元表面に現れる質量のないディラックフェルミ粒子の振る舞いをシミュレートするために、2次元正方格子上でディラック方程式を離散化する、フェルミオンダブリング問題を起こさない数値解法を提案している。