我试图做一个二维带结构计算,其中涉及多个方向的平面波。对于一个简单的正方形系统,我有两个可能的方向,以及{n,m}到其中一个的耦合项:{n,m},{n+1,m+1},{n+1,m-1},{n-1,m-1},{n-1,m+1},{n-2,m},{n+2,m},{n,m+2},{n,m+2},{n,m-2}。
我在想怎么做一个数组来做这个。我不知道如何定义某种基础。在我的例子中,我可以在某个任意的max n,m处切断,矩阵的维数将是(2n(max)+1)(2m(max)+1),因为我想从-max运行到+max。
所以本质上,我试图编写一个矩阵m嫒i j,这个矩阵m嫒ij给了我ij和i和j之间的耦合,其中i和j是{n,m}组合的一部分。
希望这是有意义的。这是我第一次体验python的有趣任务。
