Surface calculation

Posted on : November 13, 2002 (Wed) 15:47:16

by OKUNO

小谷様。

奥野（日本総合研究所）です。
たいへんごぶたさしております。　

2月近く前になろうとしておりますが、奈良でのCMDの際にはお世話になりました。
またいつのまにかこのような役にたつサイトができていて（日本の
第一原理計算関連では掲示板までついて、特定のプログラムまで
計算を議論できるサイトははじめてではないのでしょうか？）
業務の方がくそ忙しくなかなかAKAI?KKRをまわしたりする時間がないのですが
（LaGeO_{3}の計算もできておりません。ちょっと前までMOのプログラムやら高分子関連の計算やら
していました。）
私信モードはこれまでにして、、、、

燃料電池関連の「マテリアルデザイン」に関する質問なのですが
電池に使われている触媒　（　Pt単独が主。　最近Pt?Fe、Pt?Ni　Pt-Co　Pt-Ruなど
合金系が検討されている。　Phys.Chem.Chem.Phys. vol(3) p306 (2001)
（Pt合金）表面上の耐CO被毒性を調べるための（一酸化炭素）吸着エネルギーや、
酸素還元特性をしらべるためのPt（合金）?O_{2}表面の計算など、AKAI?KKRでは
できるものでしょうか？　擬ポテンシャル＋平面波では類する計算は多いようですが
合金（例えばPt-Fe　50atom%Fe です。）系の表面、もしくは原子層
slabの表面層だけPt-Feなどの合金にするなどの計算ならCPA?KKRの方が少ない原子数
で計算できるのではないかと、、。（計算する系として考えているのは表面＋吸着原子
になりますが。）う?ん、あんまりCPAでするメリットはないですか、、？
それとも、やはり表面系の計算は平面波をつかった計算の方が
よろしいものでしょうか？


あと、先ほどから（やっとですが）LaGeO_{3}の計算の計算をしてみようかと考えている
のですが、以下のファイルをつくってみたのですが、
Segmentation fault　
をはいて、計算がとまってしまいますが、どこかプログラム中の配列を大きくとるか、
なにか間違った書き方をしてしまっているかしているでしょうか？

こちらの勉強不足な質問で申し訳ありません。お時間あるときにお返事もらえると
うれしいです。

奥野幸洋

#---------------------------------- 0 ----------------------------------
# now # can be used as comment lines as well
#--------------------------------------------
go data/GaLa0.9Sr0.1O0.9Vc0.1O2
so 10.3651,1.00693,1.41659,90.0,90.0,90.0,
# edelt,ewidth, reltyp, sdftyp magtype record
0.001 1.2 sra vwn nmag init
update 5 100 0.020
4
#--------------------------------------------
# name of type, number of components, rtr,
# field, lmax, atomic number, concentration
#--------------------------------------------
Ga 1 0 0.0 2 31 100
LaSr 2 0 0.0 3 57 90
38 10
O_aVc 2 0 0.0 2 8 95
0 5
O_b 1 0 0.0 2 8 100
20
0.00000000 0.00000000 0.50346400 Ga
0.00000000 0.70829535 0.50346400 Ga
0.50000000 0.70829535 0.00000000 Ga
0.50000000 0.00000000 0.00000000 Ga
0.99010000 0.35414767 1.00400790 LaSr
0.00990000 1.06244303 0.00292009 LaSr
0.49010000 0.35414767 0.50638409 LaSr
0.50990000 1.06244303 0.50054390 LaSr
0.49800000 0.35414767 0.06041568 O_aVc
0.50200000 1.06244303 0.94651231 O_aVc
0.99800000 0.35414767 0.44304832 O_aVc
0.00200000 1.06244303 0.56387967 O_aVc
0.22900000 0.05241386 0.76727913 O_b
0.77100000 0.76070921 0.23964886 O_b
0.72900000 0.65588149 0.74311286 O_b
0.27100000 1.36417684 0.26381513 O_b
0.77100000 1.36417684 0.23964886 O_b
0.22900000 0.65588149 0.76727913 O_b
0.27100000 0.76070921 0.26381513 O_b
0.72900000 0.05241386 0.74311286 O_b