六方格子系に垂直に磁場を与えることでHofstadter's butterflyが現れる. 昔かいたFortranコードを発掘したのでアップロード. 自然界に現れるフラクタル構造って良いですよね. ちなみに,Hofstadter's butterflyは六方格子系以外の系に垂直磁場を与えてもみられます. 作ったハミルトニアンをLapackで対角化しています. program zheevtest implicit none integer,parameter :: d=60 integer,parameter :: N = 2*d*d !対角化したい行列の次元 integer :: lda,lw…