三角格子の設定 前回と同じく原子の配列から。設定が若干面倒くさい。 la = 0.4 # 格子定数 row = 4 col = 6 as = vec([Atom([x,y],[0.0,0.0],m) for x=0.0:la:la*(col-1),y=0.0:la*√3:la*√3*(row-1)]) bs = vec([Atom([x,y],[0.0,0.0],m) for x=la/2:la:la*(col-1),y=la*√3/2:la*√3:la*√3*(row-1)]) append!(as,bs) x軸方向には格子定数周期で、y軸方向は格子定数の√3倍周期で格子を作成する。 これに…