Please ｔeach me run AkaiKKR with supercell 48 atoms (in Japanese)

Posted on : July 04, 2017 (Tue) 07:25:21

by 岩堀よしひろ

48原子のSupercellを作りDOSを計算したいのですが，上手くspecxが走ってくれません。条件を以下に示しますのでどなたかアドバイスを頂けないでしょうか。

specxのコンパイル条件(54原子個で作成しています)
(gfortranでコンパイル成功)
parameter
& (natmmx=54, ncmpmx=54, msizmx=486, mxlmx=3, nk1x=40, nk3x=21,
& msex=35, ngmx=15, nrpmx=250, ngpmx=250, npmx=100, msr=400)

inファイルの中身
bzqlty=1で小さく取っています。計算実行するとエラーは出ず，すぐにプロンプトに戻ってしまいます。outファイル等は作成されますが，0 kBで中身は空白です。

c---(V0.9375Nb0.0625)O2
c---Tetra(P 42/m n m)
dos data/VNbO2_T_dos
c---
c brvtyp a c/a b/a alpha beta gamma
st 17.127 , 0.6333 , , , , ,
c---
c---edelt ewidth reltyp sdftyp magtyp record
0.001 1.0 sra mjw mag 2nd
c---
c outtyp bzqlty maxitr pmix
update 1 100 0.025
c---
c ntyp
3
c---
c type ncmp rmt field mxl anclr conc
V 1 1 0.0 2 23 100
O 1 1 0.0 2 8 100
Nb 1 1 0.0 2 41 100
c---
c natm
48
c---
c atmicx atmtyp
0.0000 0.0000 0.0000 Nb
0.2500 0.2500 0.2500 V
0.0000 0.0000 0.5000 V
0.2500 0.2500 0.7500 V
0.0000 0.5000 0.0000 V
0.7500 0.2500 0.2500 V
0.2500 0.7500 0.2500 V
0.5000 0.0000 0.0000 V
0.0000 0.5000 0.5000 V
0.7500 0.2500 0.7500 V
0.2500 0.7500 0.7500 V
0.5000 0.0000 0.5000 V
0.5000 0.5000 0.0000 V
0.7500 0.7500 0.2500 V
0.5000 0.5000 0.5000 V
0.7500 0.7500 0.7500 V
0.1525 0.1525 0.0000 O
0.8475 0.8475 0.0000 O
0.0975 0.4025 0.2500 O
0.4025 0.0975 0.2500 O
0.1525 0.1525 0.5000 O
0.8475 0.8475 0.5000 O
0.0975 0.4025 0.7500 O
0.4025 0.0975 0.7500 O
0.1525 0.6525 0.0000 O
0.8475 0.3475 0.0000 O
0.5975 0.4025 0.2500 O
0.9025 0.0975 0.2500 O
0.0975 0.9025 0.2500 O
0.4025 0.5975 0.2500 O
0.6525 0.1525 0.0000 O
0.3475 0.8475 0.0000 O
0.1525 0.6525 0.5000 O
0.8475 0.3475 0.5000 O
0.5975 0.4025 0.7500 O
0.9025 0.0975 0.7500 O
0.0975 0.9025 0.7500 O
0.4025 0.5975 0.7500 O
0.6525 0.1525 0.5000 O
0.3475 0.8475 0.5000 O
0.6525 0.6525 0.0000 O
0.3475 0.3475 0.0000 O
0.5975 0.9025 0.2500 O
0.9025 0.5975 0.2500 O
0.6525 0.6525 0.5000 O
0.3475 0.3475 0.5000 O
0.5975 0.9025 0.7500 O
0.9025 0.5975 0.7500 O
c---

よろしくお願いします。

[Re:01] Please ｔeach me run AkaiKKR with supercell 48 atoms (in Japanese)

Posted on : July 04, 2017 (Tue) 14:45:42

by Hitoshi GOMI

岩堀よしひろさま

コンパイルに成功しているように見えて実行すると落ちるのは、実装メモリの量の問題のときが多い気がしています。
コンパイル時に -mcmodel=medium のオプションをつけておられますでしょうか？
下記のスレッドを参照ください。
Increase natmmx without reducing other parameters
http://kkr.issp.u-tokyo.ac.jp/bbs/thread.php?id=252


またdos計算はgo計算で得られたポテンシャルファイルを必要としますが、go計算は正しく収束していますか？


それから msex=35 としているようですが、下記スレッドの報告もあるので msex=201 としておいたほうがよいかもしれません。
http://kkr.issp.u-tokyo.ac.jp/bbs/message.php?id=190




以降、落ちる問題とは関係ないと思いますが、入力ファイルを見て気づいた点をいくつか。


正方晶なので原子位置はfractional coordinateで与えるべきと思います。
0.0000a 0.0000b 0.0000c Nb
元ファイルはおそらく意図していない位置に原子があるはずです。


スーパーセル法をあまり使わないので自信ないのですが type は全ての原子位置に対して別々に定義しておくほうが安全と思います。
c ntyp
48
c---
c type ncmp rmt field mxl anclr conc
V1 1 1 0.0 2 23 100
V2 1 1 0.0 2 23 100
...
Nb 1 1 0.0 2 41 100
c---
c natm
48
c---
c atmicx atmtyp
0.0000a 0.0000b 0.0000c Nb
0.2500a 0.2500b 0.2500c V1
0.0000a 0.0000b 0.5000c V2
...
といった感じで。


ただしこれをやると色々なサブルーチンから怒られるかもしれません。
chklat.f とか gsdatp.f とか totalw.f とか。
それらのソースの ntypmx や ncmpmx を大きくすればいいだけだとは思いますが、あまりやってみたこと無いので分かりません。

五味