VASPスレッド１

VASPを卒論のテーマにしてしまったんだけど
使い方がさっぱりわからん
教えてエロイ人！！

何それ？２げと

VASPはWien大学Hafner教授のグループが開発した第一原理計算ソフトウェアである。ポテンシャルは超ソフト擬ポテンシャルであり、最新バージョンでは全電子計算に近い結果が得られるPAW法もサポートしている。
だってさ・・・
半導体のバンドギャップやら 最適な格子定数がもとめられるらすい
ちなみにマニュアル
ttp://cms.mpi.univie.ac.at/vasp/vasp/vasp.html
HPはここ
ttp://cms.mpi.univie.ac.at/vasp/

そもそも卒論で入って
VASPが研究室にすでにあった時点で、
周りの誰かが使い方を知っているのではないか。
ちなみに私はMD社が出しているMedeA経由で使っているので、
使い方を覚えるのは簡単でした。

>4
MedeAってはじめて聞いたなぁ〜
それは入力ファイル作成のサポートソフト？
ただなら つかってみたい・・・が
VASPの出力ファイルのマニュアルの日本語訳とか
ってどっかにないかなぁ・・

MedeAはたしかVASPとLASWにＧＵＩが付いて
企業が買う場合一千万するとか聞いたが...
うちの会社では買ってもらえなかったよ。
仕方ないから自作するはめに。

VASPはpseudo potentialが自作できないらしい
(pseudo potentialが無い元素は、Vienna大に特注するようだ)
から、ご注意。

CdSe（Wulzite)をマニュアルと見よう見まね
でやってみたらバンドが２５本でてきたｙｐ
8本よねぇ 普通・・・・
今日は帰ろう
それからうちの研究室にはMedeA
はなかったよ。
＞７
thanks
でも何のことか理解できない漏れにorz

Abinitっていうフリーソフトつかって
みたら？　けっこういいよ。

abinit ユーザーがいた！
わーい。

>>11
あなたもAbinitユーザ？　よろしこ。
どんな計算していますか？

売り物のVASPとABINITの比較って
誰かしてくれませんか？　
僕はABINITしかしらないです。

Abinitおんなじ物質でVASPと比較しながら　やってみるとおもしろいかもなぁ　多分入力ファイルとか一緒の使えそうだし、ちょっとダウンロードしてみるよ〜

VASP : GGA-PAWを使える
ABINIT : 基本的には LDA-norm conserving PP
　　　　　(GGAなどは一部の原子について公開されている。
　　　　　他は、FHIコードなどで自分で作らないといけません。
　　　　　GGA-PAWも徐々に使えるようになっているみたいです。)

さすがに有料ソフトだけあって、
シングルポイント計算で比べると、VASPの方が速い。
構造最適化はVASPのほうが圧倒的に上手。
速いのはPAWで平面波のカットオフエネルギーを
小さくできるのが効いてると思う

ちなみに入力ファイルは全然違うよ。
ABINITのほうが、分かりやすい。

構造最適化がVASPが上手いってのは
どういう点でしょう。
僕はVASPを使ったことがないので
ぜひ知りたいです。
とにかく収束しやすいのですか？

差がはっきりと分かるのは、
格子サイズ(a b c軸長など)の最適化です。
収束しやすいです。

VASP
FHI98md
ABINIT
CASTEP
PHASE・STATE(和製)
OSAKA2000(和製)
などなと
どれが優秀なのか議論しよう。

すまんそ　よく見たら
シミュレート板があったのね
内容的にそっちっぽいので
そっちに移動するね

ttp://science2.2ch.net/test/read.cgi/sim/1070718487/

>>8
input fileでバンドの本数とか指定してない？
まぁoutput fileでfermi準位が出てるから多く出てても無問題だと思われ・・・

oh,yeah!

>>　16　あ、そういやABINITでやってみたら（体積は一定）
収束しなかったんだな。あと、VASPコードについている
ELASTIC　BAND　MEDHOD　っていったいどういうものですか？
なんか拡散とかを計算するときに便利なようですが。

それと、コード自体VASPって読みやすいんでしょうか？
ABINITなんてずいぶんしっかりしたコード書いているな
って思ったんですけど。

>19
なるほどINCARでバンドの本数指定できるのね・・・
まだINCARファイルのマニュアル全部読めてなくて

CdSeのWulziteだと4つの原子で計算することになるけど
合計16個の荷電子が存在するからバンドは16本出て
くるのかな・・・？

VASPのマニュアルが読むの
面倒なひとはここちょっといい
かなと思ったよ。（まぁ自分以外はいないと思うけど）
ttp://www.aki.che.tohoku.ac.jp/~tsukidate/index-list.html
入力ファイルと出力ファイルの説明が分かりやすくでてる
細かいところをマニュアルで補足していけば使いこなせて
いけそう。

>20
基底状態のトータルエネルギー計算とかなら
上はそんなに必要ないのでは？

バンドギャップをもとめてるんだけど・・・
VASPってバンドギャップは正確にはでて
こないらしいね（だいたい１／２になるみたい）

CdS（Wulzite)4×4×4をいまやってるんだけど
（原子256個）無理っぽいわ・・・・
メモリが足りないのかしら・・・
POSCARファイル作るの大変だったのに・・

>>26　他のコードでもバンドギャップを正確に出すのは
難しいよ。そもそもLDAの限界もあるし。

>>8
VASPスレなんてあったのか…。それはともかく、CdSeの計算例なんてマニュアルにはなかったと
思うが。まあVASPの英文マニュアルは結構いいこと書いてあるから読んでおいて損はないね。

スレ立てた卒論生は、無事にVASPを使いこなせているのでしょうか。

無事に使いこなせてるような　こなせてないような・・
INCARファイルの中の
ISTART　＝　1
で設定した時　CONTCARファイルの中の
原子位置をつかって計算をスタートする
見たい事をかいてたので

一回目
ISTART＝0で普通に計算

2回目
ISTART＝1で計算したら
計算時間は短くなってもいいような
気がするんだけど実際は時間
が余計にかかってた・・・
WARNING: random wavefunctions but no delay for mixing, default for NELMDL
ってでるから　WAVECARファイルは出力しとかないと　だめなのかな・・
どなたか　アドバイスをお願いします・・

ヘルプでISTART=1や2に関する記述を見ると、
"Continuation job -- read wave functions from the file WAVECAR"
って書いてあるから、WAVECARファイルを出力するようにしておくことは必須ですよ。

でもどうして続きの計算をやらなくちゃいけないの？もしかしてバックグラウンドで走らせておく
という技をしらない？

csh系を使っている場合には sh Enterとやってsh系(またはbash)にしておいて、
$ nohup vasp & Enter
と、最後に&を入れるのを忘れずにnohupで起動しておけば、nohup.outにログ
を出力するよというメッセージが出て、裏でタスクが走る。終わったことをしりたかったら、
スクリプトを書いておいて、最後に
echo "finished" > mail [email protected]
とでも入れておけば、ジョブが終了したらケータイに連絡が入るよ。

共役方向法でも収束しないときは、
一回切って再計算した方がよいらしい。

>>32
それはIALGO=48でなかなか収束しないときはIALGO=38でやるべきって意味？

CGは一回調べた方向は二度と調べない
ためだそうです。
OSAKA2000 のマニュアルに書いてありました。

こんなスレがあったんだね・・・

俺はExamplesというところにある、
「Surface calculations」の「bulk calculation」というのをしているんだけど
何か「wrap around error」が出ないように計算しないといけない。
vaspで計算を始めるとほとんどといって良いほどこのエラーが出るのかな。
でもこの計算ではこれを消さないといけないのよね。

とは言ってもかなり使い方が難しい・・・・・本当に誰か教えて欲しいですよ（泣）
>>24
何だか便利なサイトだねｗ



あ〜あ・・・
結構疲れますわ・・これ。

というか、誰かNGXを変更する方法を知っておられる方いますか？

>>36
INCARファイルでNGX=31とか指定する。NGXはINCARファイルのパラメータ。

>>36
だいたいそういうのは先生に聞いてみたら？

>>37
アドバイスどうもありがとうございます（ペコリ）
ただ、それをするとエラーが発生するんですよ（泣）
↓以下がそれです
[*******]$ ./vaspa
vasp.4.6.12 27Jun03 complex
POSCAR found : 1 types and 2 ions


Input/Output Error 152: File does not exist

In Procedure: pseudo..rd_pseudo
At Line: 158

Statement: OPEN
Unit: 10
File: POTCAR
[*******]$
ファイルが存在しないって・・・Σ（￣□￣；）


結局、NGXは変化させずK-pointsを増加させることによって
wrap around errorは回避できたッぽいですが・・・教授からNGXの値も変化させるとか聞いていたので不安です・・・


>>38
それが・・そういう状況でもなくて（泣）
他の事はちょくちょく教えていただきましたし。

すみません。
39のエラーメッセージは私の勘違いでした。
POTCARファイルを用意していませんでした・・・・・

再度行ってみたところ
↓以下のようなエラーが出ました
[******]$ ./vaspa
vasp.4.6.12 27Jun03 complex
POSCAR found : 1 types and 2 ions
Error reading item 'NGX' from file INCAR.
Error code was IERR=5 ... . Found N= 0 data.
[******]$

う〜む、マニュアルでも見てみるか・・・

>>40
sage進行でお願いします。

どんなINCARファイルを書いたのかが問題なのに、それを示していなかったらアドバイス
も受けられないよ。

こういうエラーに対しては、例えば1行に複数のパラメータを書いてあったりしたら、
それを1行ずつに書き直すとかして少しずつ原因を追求するしかないです。

>>41
あと、最近はsourceforgeでvaspviewっていうソフトが公開されているので、使ってみたら。
このプログラムはELFCARファイルを使ってvaspの計算結果を表示します。ELFCAR
ファイルはデフォルトでは出力されないので、INCARファイルに"LELF = .TRUE."っていう
行を入れることをお忘れなく。

>>41
>>42
どうもありがとうございます。
失敗したときのINCARファイルの中身は
↓
System = ****

ISMEAR = -5
PREC = High
NGX = **

となっています。


取り合えず今はk-pointsを増加させて何度か計算を行っているところなんです。
エラーは回避しているのでひとまずはこれで試してみようかと・・。
取り合えずtetrahedron methodを用いて計算しています。

>>43
どこがまちがっていたのかさっぱりわからんがまあいいか。
数字を漢字で書いてたりでもしたんだろうな。

>>44
というかNGXはINCARでは指定できないぞ？何言ってんだ
無知だから分からないのは当たり前

>>45
In VASP.4.X all parameters are determined during runtime, either defaults are used - Sec. 7.7 or 6.1 - or NGX etc. are read from the INCAR file, see Sec. 7.2).

>>45
は、はずかしいやつ。

auto meshの種類のなかにある
Monkhorst-packってなんなんでしょ？


age

特殊点法のMonkhost-Packのことでしょう。
PRB 13,5188(1976).

すっごいくだらないこと教えてくれ。
俺は初心者だが、k点のシフトってどういう場合に
指定すればいいんだ？　Monkhost-Packってのが関係するのか？？

>>50
Monkhost-Pack タイプのやりかたが標準みたいですが、
FCCやBCCなんかだと別のk点シフトの方が良いこともあるらしい。
自動でやってくれるし、
あんまり気にしないようにしてます。

おまえ等全員

みなさんどれくらいのスペックのマシンで計算していますか？
計算している系と条件、計算時間も付けて頂けるとありがたいです。

>>53
そんな身元が割れるようなレスつけるわけがない。