バンド計算について語るスットレ
■ このスレッドは過去ログ倉庫に格納されています
0001名無し三等兵
NGNGｵﾌﾞｲｴｰｸﾄ
0002断固自営業を阻止するぞ！有限会社
NGNG0003名も無きマテリアルさん
NGNG0004ガイシュツ ◆atMKiSyUTU
NGNGttp://www.geocities.co.jp/Technopolis/4765/index.html
0005名も無きマテリアルさん
NGNG一箇所だけさっぱり分からんかった．
0006名も無きマテリアルさん
NGNG0007名も無きマテリアルさん
NGNG有効質量も計算したかったからk≠0も必要で，しかも(110)面...
まだ，スプリットオフバンド考慮してなかったとはいえ
そんなもん地方大学の学部生に分かるわけないやん
0008名も無きマテリアルさん
NGNG0009名も無きマテリアルさん
NGNG>>7
材料系は？
0010名も無きマテリアルさん
NGNG圧縮ひずみGaAsとInGaAs/InP
0011名も無きマテリアルさん
NGNG0012名も無きマテリアルさん
NGNG3次元のしか書いてないYO｡･ﾟ･(ﾉД`)･ﾟ･｡
0013名も無きマテリアルさん
NGNG0014名も無きマテリアルさん
NGNGそれぞれの電子のエネルギーを真空準位でそろえた上で比較しないといけないように思うのですが、
間違っていますでしょうか？
同じ結晶の計算でも、マフィンティン半径を変えるとエネルギーの基準が変わりますよね？
0015名も無きマテリアルさん
NGNG漏れもそうかと思ったんだが・・・これだけの情報じゃわからんよな
0016名も無きマテリアルさん
NGNG違う。一電子波動関数の波数がk。
バンド構造は、その一電子波動関数のkと固有値の関係。
0017名も無きマテリアルさん
NGNGThx.
でもその波数ｋってのがいまいち物理的に何なのかピンとこないんだ。
局在電子系ならその電子って、d性とかp性とか言ってその状態を指定できるでしょ？
平面波基底（ブロッホ関数てこれの足し合わせみたいなものって理解でOK?）における波数ｋはどういうものに対応するかが良く判らないんだ。
0018名も無きマテリアルさん
NGNGブロッホ状態の波動関数phiは、
ブロッホ関数uと、波数kの平面波の掛け算になっている。
phi(r) = exp(ikr)×u(r)
（u(r)は結晶周期の関数＝逆格子ベクトルの平面波の足し合わせ。
もちろん、u(r)はkごとにちがう）
この平面波exp(ikr)の波数kが、波動関数phiの波数になる。
∵kの演算子は-i∂/∂rだから。
要は波動関数の波数kとは、波動関数からブロッホ関数uを除いた、
ベースの平面波に相当している。
（添削希望）
001914
NGNG0020名も無きマテリアルさん
NGNG先輩や先生に聞いてみては？
0021名も無きマテリアルさん
NGNGまた、そのY_lmを変更するにはどうすればいいですか？
マニュアル読んでも分からなかったよ…
0022名も無きマテリアルさん
NGNG0023名も無きマテリアルさん
NGNG次のワークショップは日本ですが。
0024名も無きマテリアルさん
NGNGまじでつか！
あと、WIEN2ｋって確か将来的には輸送についても議論できるよう拡張するとかBlahaさんが言ってたよね。どうなるんだろう？
0025名も無きマテリアルさん
NGNGこれですね。
ttp://cms.mtl.kyoto-u.ac.jp/wien2k-ws.htm
0026vaspman
NGNG0027名も無きマテリアルさん
NGNGE-kのバンド構造から状態密度を計算したいのですが。
自分では精度良くプログラムできないみたいです。
0028名も無きマテリアルさん
NGNG0029名も無きマテリアルさん
05/02/17 16:05:470030名も無きマテリアルさん
05/02/24 19:22:340031名も無きマテリアルさん
05/03/06 13:16:10マジお勧め
0032名も無きマテリアルさん
2005/04/02(土) 08:24:16量子井戸の歪みってどんなふうに計算にとりいれればいんですかね？？
だれか計算した人いませんか～～～？
0033名も無きマテリアルさん
2005/04/03(日) 00:21:35あとは普通に量子準位を計算すればいいんじゃね
0034名も無きマテリアルさん
2005/04/03(日) 02:25:50だれかした人はいますか？
003532
2005/04/03(日) 02:56:26ありがとうございます。
ほんとそのとおりだと思います。
ただ有効質量なんかも変わってこないのかな～～？
って思ったんですけど、そんこたぁないですかね？？
003632
2005/04/03(日) 02:57:12連続書き込み、失礼。
003732
2005/04/03(日) 03:37:38わかりました。
バンドオフセットの変化量も歪みのシフトから求められますね。
お恥ずかしいかぎりです。
ギモンは解決できました。
どうもありがとうございました。>33
またなにかありましたら書き込ませていただきます。
0038名も無きマテリアルさん
2005/04/19(火) 23:32:39これとかいじってみた人いる？
0039名も無きマテリアルさん
2005/04/20(水) 01:19:09化学屋はGaussianって聞くけど、物性屋にも使える？
誰か他に情報くれ～。
0040名も無きマテリアルさん
2005/04/20(水) 12:19:51俺も気になる。MSとGaussian、他には？
どっちがいいの？
0041名も無きマテリアルさん
2005/04/20(水) 18:34:470042名も無きマテリアルさん
2005/04/21(木) 00:35:560043名も無きマテリアルさん
2005/04/23(土) 16:15:21他の大学は何使ってんの？プログラムじゃなくソフトは何がいいんだろう。
0044名も無きマテリアルさん
2005/04/23(土) 22:37:59ＯｐｅｎＭＸ
つか、京大の全研究室はMaterial Studioしか使っていけないと決められている訳ではないだろうに。
0045名も無きマテリアルさん
2005/04/25(月) 01:11:110046名も無きマテリアルさん
2005/04/29(金) 01:49:29ところで材料系の人はPCシミュレーションではどんな事やってんの？
バンド計算以外でなんか使えるのある？
0047名も無きマテリアルさん
2005/04/29(金) 12:34:35結局ポテンシャルの数値自体の評価方法を洗練させないと、どう信じてよいのか判らないので、漏れの研究テーマの一貫としては没った。
最近考えているのだけれど、超伝導の電子構造をバンド計算みたいな方法で評価する方法って誰か知らない？
0048名も無きマテリアルさん
2005/04/30(土) 01:55:21MSのCASTEP並に普及度が高いと思うのだが．
オープンソースならABINITが比較的有名かな．
Ultrasoftな擬ポテンシャル(PAW)も使えるみたいだが
実際総合的に使える代物なのかどうかは、わからん．
Ｇａｕｓｓｉａｎは周期系の計算がしょぼいので
クラスター（しかも数十原子クラス）の計算にしか実用上適していない．
でも高次電子相関が扱える点では魅力的だ．
価格は比較的良心的だけど個人で買える額ではないだろう．
化学系では他にＧＡＭＥＳＳってのがあって，
こちらはフリーで使えるが周期系の計算には対応していない．
でも相対論計算や励起状態計算にはＧａｕｓｓｉａｎよりも
優れた機能を搭載しているらしい．
化学系第一原理計算ソフトはここで聞いたら詳しく教えてくれると思うよ．
http://science3.2ch.net/test/read.cgi/bake/1104533051/l50
>>46
みんなで盛り上げていこう． というわけでage．
0049名も無きマテリアルさん
2005/04/30(土) 15:28:45あんさん、いい人だね。涙でた…。
005048
2005/05/04(水) 19:45:52GAMESSの相対論計算は国産(東大？)の相対論用基底が
使えるというだけで，特にGaussianとはあまり変わらないみたいだ．
>>46
自分は第一原理計算で得られた電子密度分布やPDOS，
安定化エネルギー等から，無機材料が関係する化学反応を
解析したりしている．
なので最近化学系のソフトも使い始めているんだけど
なかなかこれが難しい・・・
0051名も無きマテリアルさん
2005/05/14(土) 20:14:11005251
2005/06/07(火) 23:44:380053名も無きマテリアルさん
2005/06/08(水) 21:16:36行きたかったが、一週間休むと言ったら上司に不許可下された。orz
0054名も無きマテリアルさん
2005/06/26(日) 01:35:13扱ってるのは銅の合金など。
次のパラメータの意味がよく判りません。
・相対論の有無
・交換相関ポテンシャル（mjw,vbh,vwn）
判る方います？
0055名も無きマテリアルさん
2005/06/26(日) 08:00:10漏れはAKAI-KKRは使ったこと無いが…
相対論は、内殻軌道だとか重元素の計算する時に「有り」にするもの。f系も入れた方が良いだろう。
交換相関ポテンシャルは、LDAかGGAの選択の事だろうか？
一般にGGAのほうが計算精度が優れると思えば良いが、GGAと呼ばれる方法は色々提案されててこれ使っておけというものはあまり無い。
つーか、>>4のサイトを小一時間眺める事を勧める。
0056名も無きマテリアルさん
2005/07/04(月) 21:04:30nrlで非相対論計算
sraでスカラー相対論計算
sralsでz方向を量子化軸とした相対論計算（スピン軌道相互作用）
mjw,vbh,vwn
これらはすべてLDAでタイプが異なるだけ。
ところでAkaiKKR以外でCPA計算ができるバンド計算パッケージがあれば
紹介してほしいのですが。
0057名も無きマテリアルさん
2005/07/05(火) 16:16:42gga91とするとGGAも使える。
しかし赤井先生がコーディングしたものでないらしく、結果に責任はもてないとおっしゃってた。
0058名も無きマテリアルさん
2005/07/05(火) 16:52:20高齢化ビジネスにも最適！
初心者でも簡単運営
お客さんにも便利屋さんにも喜ばれるサイトです
http://money4.2ch.net/test/read.cgi/manage/1120100670/120
0059名も無きマテリアルさん
2005/08/08(月) 17:20:21出てきた結果に一気一優する。それが大半のソフトユーザー。
第一原理計算の大きな問題点だな。
密度汎関数理論とKohn-sham方程式から勉強してください。
それが、どのような位置づけなのか？多体理論から見せかけの一電子方程式に
焼き直して。。。。。。そんな事も知らずに、結果を議論しちゃだめですよ。
0060名も無きマテリアルさん
2005/09/30(金) 23:10:05一から勉強したいんだけど、時間がないよ。
職場では、計算結果を出す事だけを期待されている。
適用範囲外のパラメータで計算してしまって、
嘘の結果を出してしまったとしても、自分も上司も判断できない状況。
かなりヤバいと思う。
006160
2005/09/30(金) 23:22:49今の会社に入って、VASPで色々計算しろって言われてるんだけど、
僕がやる理由は職場で量子力学を知ってるから。
とは言ってもさ、僕も大学では流体をやってたんだよね。
量子論も教科書レベルの知識しか無いしさ。
とにかく、バンド理論がまだ理解できないんだよね。
VASPの計算はガンマ点だけで計算すれば良いんだよ、
なんて聞いたんだけど、何で良いのかさっぱり分かんないよ。
会社では、論文も教科書もいっさい読むなと言われてるしね。
勉強は仕事じゃないので家でやれってさ。
なんかヤバいよね。こんなんでいいのかしら。
0062名も無きマテリアルさん
2005/10/01(土) 00:48:15ちがうじゃねーか。バカ。
0063名も無きマテリアルさん
2005/10/02(日) 09:53:42VASPのマニュアルを隅々読んだほうがいいよ。あとNIMSのK林先生のサイトを
ヒマな時にでも眺めることを勧める。最近は格子振動を考慮して熱力学的物性値
も出せるから、結果については現実の再現性を重視した方が良いかもだよ。
勿論、どんな系を計算するかにもよるけど。あと、セルのサイズによって取るべき
K点数は変わりますよ。
0064名も無きマテリアルさん
2005/10/07(金) 00:12:05CMDワークショップとかCAMMフォーラムに参加してみるのは？
まぁ>>61の様子からすると上司から許可が下りないかもしれないけど・・・
孤立分子とか流体などを平面派基底で計算する場合はGamma点
だけで良いので，多分そういう話を耳にしたのでは。
無機結晶等の周期系ならk点を増やしていって計算結果
(エネルギー・格子定数etc)があまり変動しなくなる点を選んで
計算するのが妥当な方法だと，とりあえずはそう覚えておいて，
あとはやりながらその理由を勉強すればいいんじゃないでしょうか。
0065名も無きマテリアルさん
2005/10/08(土) 22:41:24ドープしたり欠陥作ったりした系では必要だとは思いますが、
対称性を落として構造緩和すべきでしょうか？。
0066名も無きマテリアルさん
2005/10/10(月) 23:40:55どんな物理量を評価するかによるのでは？
0067名も無きマテリアルさん
2005/10/11(火) 00:16:26基本的にはTotal energyとpartial DOSです。
うまく進めば実験系との比較をいろいろとやりたいとは思ってますが。
006866
2005/10/12(水) 23:58:38一点計算でも定性的な比較はできると思う．
でももしあとでドープした系や結果のある系も計算するならば
ドープしてない系も構造緩和したほうが同じ基準で比較できるよ．
006966
2005/10/13(木) 00:01:10結果のある系 ⇒ 欠陥のある系 の誤りです。
0070名も無きマテリアルさん
2005/10/16(日) 14:05:44実際どの程度効果ある？
0071名も無きマテリアルさん
2005/10/19(水) 16:16:58なんか神様の道具みたいだけど、そんなに信頼性ない部分もあるし。
特に、計算が速いソフトはそれだけ精度を犠牲にしてる訳で、
その精度が微妙な所が、実はみんなが議論したい所だからね～。
実験と合わなくても胸をはって
『第一原理の結果です、間違ってるかもしれません！』
と言えるかどうかですよね。
構造緩和って言ったって、勝手にやってくれるしな～。
昔、某有名ソフトを使って同じものを計算して
結果が違ってたときはかなり笑ったが。
使い手が一番重要っすよ！
0072名も無きマテリアルさん
2005/10/20(木) 00:44:51で？
0073名も無きマテリアルさん
2005/10/21(金) 02:37:40VASP : 超お勧め
ABINIT: いろんな人がいろんな機能をばんばん取り込むので多機能
USPP/PAWは最近やっと取り込まれたがまだまだ使えない。
Pwscf: 構造緩和しないならUSPPも使えるのでおすすめ。MD得意そう。
PHASE: かなりよさげ マニュアル日本語だし USPPポテンシャルが充実。
OSAKA2000 マニュアルがすばらしい 一回読んだらわかった気になれる（かも）
ま、ユーザーとしてはざっくり理論を勉強した後で、
高機能な測定装置 くらいのスタンスで使うのがよいかなと。
0074名も無きマテリアルさん
2005/10/21(金) 08:47:34今、Wien2kつかっててそこそこいいんだが、
構造緩和も不器用。
対称性落としたりsupercell組んだりするとしゃれにならないほど時間かかる。
当たり前だが…
自分がやる系はフルポテンシャルにこだわる必要はおそらく無いから
psodopotentialで計算したい。
やっぱりCASTEPかVASPあたりなんだろうな。
PhaseはHubberd potentialとか考慮できるのかな？
0075名も無きマテリアルさん
2005/10/21(金) 14:24:36通常、擬ポテンシャルと対になるのは全電子だよなー。
ちなみにスペルも違ってるよ。
0076名も無きマテリアルさん
2005/10/21(金) 14:27:25適応範囲についてはむしろLDAの方が広いぞ。
電荷分布の位置依存が強い場合にはGGA、ってのが一般的な思想だよ。
007773
2005/10/21(金) 23:01:43Pwscf : LDA+U
ABINIT: GW近似
VASP : LDA+U / GGA+U
0078名も無きマテリアルさん
2005/10/21(金) 23:44:53今までこんなミスしてたのかorz
ついでに全電子=fullpotentialだと思ってた。
なら全電子に対応する英語は?
0079名も無きマテリアルさん
2005/10/22(土) 02:25:58そのまんま。
full potential と対を成すのは球対称近似。
マフィンティン(muffin-tin approx.)とかASA(atomic-sphere approx.)とか。
0080名も無きマテリアルさん
2005/10/22(土) 14:58:50PWscfはUSPPで構造緩和できるよ
FPMDがノルム保存型じゃないとできないだけ。
008173
2005/10/22(土) 21:44:27やってますけど、
セルの緩和（a,b,c,alpha,beta,gamma)が
"へたくそ"なんです。CG法とかないから。
0082名も無きマテリアルさん
2005/10/25(火) 13:25:57遅いって事かな？
0083名も無きマテリアルさん
2005/10/28(金) 13:29:020084名も無きマテリアルさん
2005/10/28(金) 18:47:170085名も無きマテリアルさん
2005/10/28(金) 22:53:26crystal とかはどうでしょう。
周期境界条件使えるガウス関数基底ってのは
意外と珍しいのでは。
0086名も無きマテリアルさん
2005/10/29(土) 17:55:25それで何ができてどんな嬉しい事があるのかを
語った方が、いいかも…
0087名も無きマテリアルさん
2005/11/01(火) 00:52:21結晶構造、逆格子、バンド計算、電子密度
などに関しての自分に対する教育効果が高かったこと。
それに自分で気になった物質を買ってきたり、
合成したりすることと比べると遥かに手軽に計算できるしね。
「おお 確かにここの電子密度が高くなってる」とかさ。
0088名も無きマテリアルさん
2005/11/01(火) 22:27:50うんこ、うんこへ。
このうんこはうんこだ。まずうんこを飲んでうんこ。
うんこ、「うんこ」なんだ。うんこ。
うんこの顔もうんこまでと言ううんこ、うんこしてもらおうとも思っていない。
うんこ、このうんこを見たとき、うんこは、きっとうんこでは言い表せない
「うんこ」のようなうんこをうんこしてくれたはずだ。
うんこな世の中で、そういううんこをうんこしないでうんこ、そう思って
このうんこをうんこしたんだ。
さあ、うんこを食おうか。
このレスを見た人は酷い下痢になる
0089名も無きマテリアルさん
2005/12/14(水) 00:21:29プログラミング言語何使ってます？
0090名も無きマテリアルさん
2005/12/18(日) 17:05:18Fortran90 > C > C++ の順じゃない？
俺はCで書いてるけど．
0091名も無きマテリアルさん
2005/12/19(月) 19:02:40つーかFortran90，sma4，Wordって感じだわな
俺的には C++，gnuplot，LaTeX をよく使うんでちょっと嫌なんだけど仕方ないと思ってる．
Compaq Visual Fortranってのがあるんだが，これでGUIなプログラム書いたことある人いるのかね．
0092名も無きマテリアルさん
2005/12/20(火) 14:11:03一から自分で書くときは90/95で書くけど、それよりも
古より伝来のソースをいじくる方が多いから・・・
ちょっと出力を解析したりするときなんかはPythonを使う。
0093名も無きマテリアルさん
2005/12/21(水) 01:33:500094名も無きマテリアルさん
2005/12/21(水) 03:00:56目的次第。目的を十分達成できるなら手段なんか
なんでもいい。ツールを使いこなすのが目的ではなかんべ？
0095名も無きマテリアルさん
2005/12/21(水) 03:19:25osaka2kは使える？やってみようと思っているのだが。
VASPはちょっとかじったが。
自力プログラムは77。出力制御は気になるが。
Visual Fortranは気になってた（このスレッドとぜんぜん関係ないが）。
どんなの？
0096名も無きマテリアルさん
2005/12/21(水) 20:03:15自分で作れるし、いまどきねぇ。ノルム保存ならabinitにすりゃいいんじゃない？
あとは、ウルトラソフトなSTATEとかPWSCFがいいんじゃない？
STATEはオープンなんかどうか知らんけど。
PWSCFなら色々計算出来るし、速いしね。
0097名も無きマテリアルさん
2006/01/04(水) 21:26:330098名も無きマテリアルさん
2006/03/18(土) 21:53:170099名も無きマテリアルさん
2006/03/19(日) 07:09:54もうちょっと考えて名前付けろよな
0100名も無きマテリアルさん
2006/04/05(水) 02:11:22実際にはマッフィンティン半径を設定してるんだけど
どういうこと？
0101名も無きマテリアルさん
2006/04/05(水) 19:14:10MT近似ってのは球対称近似のことだから。
MT球を使ってても、球対称近似をしなければフルポテンシャル。
0102名も無きマテリアルさん
2006/04/06(木) 00:30:13それならどうやって、球内と外のポテンシャルを連続にしてるんだろ？
0103sage
2006/04/07(金) 01:26:28LAPWでは、
球内と球外の波動関数の２次微分までが球面上で連続となるように展開係数が決められる。
なので、電子密度（波動関数の２乗）は球面上で１次微分まで連続。
なので、ポテンシャルも１次微分まで連続。
0104名も無きマテリアルさん
2006/04/07(金) 20:07:12御説明ありがとうございます。
球面上でちゃんと連続になるようにしてるんですね
当然といえば当然ですけど…
あと、細かいことになりますが、
Wien2kはAPWとLAPWを使い分けられます。
この使い分けもどうしたらよいのかまだ良くわかっていないorz
0105名も無きマテリアルさん
2006/04/08(土) 01:35:00球面での、いわゆる接続問題に落とし込んでるわけ。
APWとLAPWの違いはL=線形化（動径関数のエネルギー依存性が）されてるかどうか。
線形化しないと行列の対角化法が使えないので計算が大変になる。
逆に線形化すると、ゴースト・バンドが出たりすることがある。
使い分けは・・・Wienを使ったことないから自信ないけど、通常はLAPWでやって、
どうも結果がおかしいってときだけAPWを試してみる、みたいな感じでいいのでは。
0106名も無きマテリアルさん
2006/04/09(日) 18:45:18解説ありがとうございます。
まだまだ使いこなすには程遠いですががんばってみます。
0107名も無きマテリアルさん
2006/04/11(火) 17:36:010108名も無きマテリアルさん
2006/04/11(火) 19:49:20お前アホか。LAPWでは、波動関数は一次微分まで連続だぞ。
二階微分まで連続なのは普通はSLAPW-3だよ。
D.Singhの論文の話してんの？
これは、マルチウィンドゥだからゴーストバンドの恐れはひとまず無いよ。
やばそうなバンドは全部Valenceにすりゃいいんだから。
方法論的にはね。
フルポテンシャルは
球内のポテンシャルの球対称成分を用いてLAPW動径関数をつくって
その後に、球内のポテンシャルの非球対称成分の評価と
Weinertの方法で球外のポテンシャルをフルポテンシャルにしてるんだよ。
>>105
APW+loは、LAPWみたいにエネルギーパラメータの所で
動径関数と、エネルギー微分関数を作って、そこから
球内に局在したloを作る。あとは、これに加えてAPWをあるエネルギーパラメーターで
解いた動径関数を合わせて基底関数を作るんでしょ。。
本質的には線形法だよ。APWをもろにやってる訳じゃないよ。
APWでフルポテンシャルにするには恐ろしい困難があるんだよ。
Wien2kはLAPW+LOとAPW+loを混ぜてるんだよ。
各lに対してどちらを選択するかの自由度があるだけ。
てか、論文読めよ。KoellingとArbmanのやつとか。
Andersenのとか、D.SinghのLAPW+LOとか、
E.WimmerとM.Weinertのフルポテンシャルとか、
最近のAPW+loとかさ。こんなとこで、嘘ばっかり教えられてもしょうがないよ。
103は本気であほ。勘弁してくれ。
0109名も無きマテリアルさん
2006/04/11(火) 23:59:52105だが、APWとLAPWの違いに対する回答としては間違っちゃいまいｗ
つか、104の言うAPWが、実はAPW+loであって、それが"An alternative way of linearizing"だなんて
んな細かいこと知らないよorz
線形化してんなら、ちゃんとLつけろよヽ(`Д´)ﾉ
0110名も無きマテリアルさん
2006/04/12(水) 02:31:400111名も無きマテリアルさん
2006/04/12(水) 02:39:23煽りも程ほどに。
なんでも論文まで漁って調べろってのは乱暴では？
甘えと切り捨てればそれまでだけど、
限られた時間の中であらゆる知識取り入れて
成果出していかなきゃならんのだから
ずべてに対してそんなに丁寧にやってられんことぐらい
あんたもわかるだろう。
まぁ、他人から聞いた情報を鵜呑みにして確認を怠るのは以ての外だが。
0112名も無きマテリアルさん
2006/04/12(水) 09:50:24普通は出ないだろうけどね。
99.9%は。
>>111
すまんです、言い過ぎな部分ありました。
まあただ、ここや、ウィキペディアしかりインターネットの
間違った情報をあさるよりは、専門家（決してパッケージユーザーでない！）
や、ソフトの説明書を読んだ方がいいですよ。
0113名も無きマテリアルさん
2006/04/12(水) 15:56:43そんな事言い出したらこのスレの存在意味がなくなってしまうではないか。
間違った情報がここに記されたならそれを訂正し切磋琢磨しようくらいの気持ちで居てください。
0114名も無きマテリアルさん
2006/04/12(水) 17:03:39問題は切磋琢磨しあうだけの有識者の数
108ぐらいの人が頻繁に覗いて書き込んでくれるといいんだけど…
0115103
2006/04/15(土) 01:24:12皆様、大変失礼いたしました。
0116名も無きマテリアルさん
2006/04/15(土) 15:07:24ところで我々人類がタイムワープの技術を開発し、過去の世界に
干渉できるようになり、そして過去に干渉した時に現在の世界に
影響があるのではなく、平行時間軸ができてそこに新たな歴史が
刻まれるとしたら、それは４次元世界ではなく５次元世界になる。
人類がタイムワープの技術を手に入れれば、人類が支配できるのは
４次元世界ではなく５次元世界になるわけだ。
そしてこの５次元世界が神の作った仮想世界であり、神は独自の
３次元空間の中で我々の５次元世界を観察しているとすれば、神の
世界は我々から見て８次元の世界にいることになる。
神は８次元世界に存在する！
0117名も無きマテリアルさん
2006/04/16(日) 01:13:32やっぱりCAPTEPかVASP買ってもらうしか…
0118名も無きマテリアルさん
2006/04/22(土) 00:11:43意外となんとかなるよ。 コンパイルには苦戦するかもしれんけど。
VASPに比べると遅いけどVASPではできないNPTでの
第一原理ＭＤができるのがｲｲ
CASTEPは・・・高すぎて並列計算ムリ。 ポスト処理楽だけど。
0119名も無きマテリアルさん
2006/04/22(土) 00:24:47値段が高いってこと？まさか1ノードごとにライセンス必要なの？
0120名も無きマテリアルさん
2006/04/22(土) 01:47:13結構、悪名高いよね。
ノードライセンスって聞いた事あるけど。。
とりあえず、VASPもあった方がいいんじゃない？
0121名も無きマテリアルさん
2006/04/22(土) 07:17:11Wien2kでスラブ入れてとか…
0122名も無きマテリアルさん
2006/04/22(土) 07:21:01そうですか…なんとかなりますか…orz
なんとかコンパイルは出来てます。一応exampleはreferenceと同じ
出力を得られるものが出来てます。ifor9+mkl721でチョット苦労しましたが
とりあえずNEBが出来るとこまでがんばってみます。
とりあえずMLのアーカイブ漁るかな…
0123名も無きマテリアルさん
2006/04/22(土) 07:23:49げ…なんてこった…
となるとVASP一択だな
どうりでみんなVASP使うわけだ…
0124名も無きマテリアルさん
2006/04/22(土) 10:48:20CASTEPのコード自体はUKではアカデミックに無料で提供してるみたいだし
そんなこと無いんじゃないの？
実際どうかはしらんが…
0125名も無きマテリアルさん
2006/04/22(土) 17:18:37計算するものによるんじゃないですか。
金属表面に酸素の吸着とかならFLAPWでもクラスタなしで
できそうな感じがしますが。。（遅いですけど）
いや、やっぱりPAWでも十分ですよ。
全てFLAPWで計算する必要はないでしょう。
USPPだと、まだまだ磁性の系は怪しいですからね。
FLAPWだと、計算重すぎるし。
するとやっぱ、PAWですかね。てか、VASPか。
別に、Kresseと仲いい訳じゃないが、いいソフトだよね。
0126名も無きマテリアルさん
2006/04/22(土) 18:50:35ウーン
やっぱりVASP買ってもらうか…
3000EUROか
やっぱり年度始めだし先延ばしにされそうorz
0127名も無きマテリアルさん
2006/04/23(日) 09:40:13常温での物性を計算したい場合、
平均構造がちゃんとわかってるなら
構造緩和はしないほうがいいよね？
常温での構造とは違ってて当たり前だし…
0128名も無きマテリアルさん
2006/04/23(日) 14:49:22常温での構造がどれくらい違うかも問題ですよね。
構造相転移みたいな事があって大きく構造が違えば
常温の構造で色々計算する事は意味ある。
（当然、計算する物理量によりますが。。）
大きな構造の変化があると普通、構造緩和すれば
低温の構造に行き着くから常温とはかけ離れた
結果に多分、なるでしょう。
常温と構造が殆ど変わらないなら構造緩和しても
いいと思います。実験的に見積もられた平均構造と
大きく違わないと思うし、実験のある意味チェックにもなるし。
計算は低温の状態を考えている訳だけど、
普通は相転移がなければ常温との
構造はあんまり変わりませんし、
内部座標もあんまり変わらないと思いますが。。。
分かりません、場合によるかも。
0129名も無きマテリアルさん
2006/04/26(水) 00:02:54CASTEPアカデミック版は「イギリスのacademic機関」
のみ利用可能だから実質利用可能なのはMSのCASTEPだけ
>>119,123
MS CASTEPは1CPUで1ライセンス消費する
そのくせVASPよりも遅いしPAWが使えないが
全てのランタノイド元素の計算ができる
余談だがもともとはVASPとCASTEPは
同じPayne先生のコードから出発している
0130名も無きマテリアルさん
2006/04/26(水) 14:31:48Payneのが元っていってもね、
もう中身は別もんだよね。
昔はCASTEPほど洗練されたソフトでなかったわけだし。
VASPは、元はどうあれKresseのソフトだよね。
全てのランタノイドの計算が出来るけど
いい結果は出るんだろうか？
まぁ、擬ポテンシャルを作っときました程度では？
もう少し良いのかもしれないけど、
ポテンシャル作っとけば後はどうにでもなる訳で、
あんまり利点では無い気がしますが。
さすがにf電子系は擬ポテンシャル、PAWでは無理な気がします。
0131名も無きマテリアルさん
2006/05/04(木) 18:09:29あると聞いたんだけど
何を基準にして決めるの?
total energyだとUに影響受けちゃうし…
0132名も無きマテリアルさん
2006/05/08(月) 23:19:14LDA+Uで、 Uを決めて何が欲しいの？
あなたはkohn-shamの全エネルギーが欲しいのか？
励起状態、バンド構造が欲しいのか？
前者なら変分原理の範囲内でLDA+Uの範囲内で既に意味あり。
後者はGWのセルフコンシステントやらない限り、
Hartree-Fockレベルの精度は求められない。
ちなみに、abinitやら擬ポテンシャルレベルのLDA+Uをいくら
線形応答やなんかで+Uを決めても意味は薄いよ。
方法論を知ってる人ならね。
LAPWでベストとは言えないけど、 かなりmore betterよ。
まあ、気にせずLDA+Uを適当なパラメータでやっと来なさい。
最近のcond-matでも今後後悔しそうな論文は一杯発表されてるから。。。
VASP,wien買いました。LDA+Uで最適化しました。
こうなりました。第一原理の結果ですから。。。。正しいです。
あんたら、そこらの安いソフトがどんな近似や数値計算の横着やってるか知っての事？
そもそも、あんた等の精度でその議論出来るのって？
パソコンって楽でしょ。待ってれば結果でるから。
みんな、よく考えて。そんな簡単な研究があり得るか？
0133名も無きマテリアルさん
2006/05/08(月) 23:20:15頑張ろう！
0134131
2006/05/09(火) 00:49:17私は固体を扱う化学系の実験屋の学生です。
LDA+UでUを決めたいのは、Uを実験結果のみを
よりどころとした試行錯誤で決めたくないなと思っただけでして、
Uを計算の身から見積もる方法があると小耳に挟み、
どんな方法なのかと興味を持ったしだいです。
LDA+Uにこだわっているわけではないんですが、
今対象にしている系では電子相関の近似がＧＧＡでは
不十分であることが指摘されており、
Uを計算に取り入れることで著しい改善が見られるといった
報告がなされております。
それが、その今後後悔しそうな論文なのかもしれませんが…
まだ駆け出しなので既報のデータと比較検討できるといった事情もありました。
今使用しているコード(Wien2k,PWSCF,CPMD)でも
GW近似等を利用して、計算コストに見合う
より優位性のある結果が得られるのであればむしろ喜ばしいです。
結局やりたいのは、今時分が合成して特性を評価したり、
その特性の本質を実験的に探る過程で、
一つの測定手法のように、これでしか引き出せない情報を得る。
あるいは実験から予想する結晶構造や電子密度、スピン状態が
第一原理計算の結果と対比しても妥当なものであるのかといった
検討をしてみたいと考えています。
0135名も無きマテリアルさん
2006/05/13(土) 22:36:46どのくらいのオーダーで妥協してる？
論文見てると、R_MT=1.9ぐらいでRK_max=7で計算してる。
マニュアルとかチュートリアル見てもR_MTの4倍位が妥当だとか書いてあるけど
自分でやるとRK_max=10でもまだTotal Energyは減少傾向にある。
結局やりたいのは実験データとの比較だし、
絶対的に正しい値が得られるとは思ってないから、
ある程度は適当でいいんだけどどこまでやるのが
流儀なのか知りたい。
0136名も無きマテリアルさん
2006/05/14(日) 13:14:34そこらへんは、David.J.Singhの
planewaves , pseudopotential and the LAPW methodをお読みになると良いかもしれません。
計算方法に詳しいです。
RK_max=7~9くらいでよいのでは無いでしょうか？
減少とはどういう事でしょうか、収束に向かっていないのですか？
MT半径と基底関数の数から算出されるRK_maxでも、
問題によってはこれでは不十分な場合もありますし十分な場合もあります。
比較対象と、収束度合いの兼ね合いになります。
しかし、だいたいRK_maxを7~9程度にすれば十分と思いますし、
その上でk-pointに対する物理量、エネルギーの収束性を見れば良いのではないでしょうか？
何meVまでの精度まで計算で収束させたいかという事だと思います。
ただ、RK_max=10でも20でもk-pointが足りなければ比較すべき量でない場合があります。
例えばCuの計算でもk-pointが1pointではRK_max=100でも
wien2KではL/Apw+loを採用しているのでかなりRK_maxの収束性は速いと思います。
なかなか難しいですが先にRK_maxを固定してk-pointでの収束性を見るのが普通だと思います。
ただ、wienコードはよく知りませんがK_maxは
基底関数の波動関数に対する展開を決めているので
電子密度の展開を直接決めている訳ではありません。
ですから、厳密には電子密度の展開数に対する収束性を見る必要があると思います。
この値がK_maxと同時に変化すれば、それぞれのK_maxでの全エネルギーは
比較すべき量ではなくなります。電子密度の展開数を多くとっておいてから
それを固定してK_maxを変化させる必要があると思います。
細かい話ですが。。
0137136
2006/05/14(日) 13:15:55>>Uを計算に取り入れることで著しい改善が見られるといった 報告がなされております。
なかなかデリケートな問題ですし、何を議論しているか分かりませんので何とも言えません。
上司の方等に相談されるのが良いと思います。研究題材を公開して頂けるなら
なにかアドバイスは出来るかもしれませんが。。
+Uを決めることや、GW近似をやる事よりも
（いや～GWも終了しましたね。。）
その系で扱われている適当なUの値を用いてLDA+Uの計算、（GGA+Uもほぼ同義です）
GGA,LDAの計算で実験結果を比較されてはいかがでしょうか？
あまりに実験結果とかけ離れていれば、常識の範囲内で+Uの値を変化させ
実験値に近づける事も様々な論文で行われています。
現時点ではそのレベルで十分だと思います。
GW近似はレベルの方法が色々ありますし、どれが正しいか、精度は？
など、私は手を付けません。５年後くらいには普通に使える様になってるかも知れませんが。
0138名も無きマテリアルさん
2006/05/21(日) 00:28:34RH8.0 cluster machine(20台） で CASTEP3.0 を使っていたのですが、4.0にver. up
したらEnterprise AS3.0/4.0 でないと駄目とのこと。
全部OSを入れ替えるとどのくらいの手間と金額がかかるでしょうか?
（それにしてもVASPの方が安上がりだな。）
0139名も無きマテリアルさん
2006/05/23(火) 12:02:500140名も無きマテリアルさん
2006/05/25(木) 13:22:58普通にCENTOSとかでいいんでないでしょうか？
0141名も無きマテリアルさん
2006/05/28(日) 00:08:440142名も無きマテリアルさん
2006/05/28(日) 14:02:21有難うございます。
「RH8.0のままでよいか，CENT OSにしても良いか」と
accelrys と cluster vendorに聞いたら
「駄目だ」
という返事が来たものですから...。
駄目もとでやってみませうか。
0143名も無きマテリアルさん
2006/06/05(月) 12:06:400144名も無きマテリアルさん
2006/07/14(金) 19:59:52立場で、できるだけ簡単に（というかＧＵＩ充実してるほうがいい）バンド計算
したいのですが、お勧めありますか。
GAUSSIAN03あたりを考えてはいるのですが
0145名も無きマテリアルさん
2006/07/15(土) 08:49:46ある程度使える人間雇った方がいいよ
それに関わるあなたの人件費が無駄になる
0146名も無きマテリアルさん
2006/07/15(土) 20:41:11レスサンクス
で、お勧めは？
0147名も無きマテリアルさん
2006/07/16(日) 01:34:00使える人間が習熟してる奴
0148名も無きマテリアルさん
2006/07/16(日) 20:38:37そりゃ過疎化するわな
0149名も無きマテリアルさん
2006/07/16(日) 21:00:49どんな物理量に注目してるか書かんとな…
あといちいちあげんなよｶｽ
0150名も無きマテリアルさん
2006/07/17(月) 02:29:53VASP
FSIS
Materials Studio
WIEN2k
PWSCF
CAMP
0151名も無きマテリアルさん
2006/07/17(月) 19:05:21プログラム自分で組んでやってる人って
どのくらいいるの？？
0152名も無きマテリアルさん
2006/07/24(月) 18:44:22Windows上でFORTRAN90を使って、プログラムを組んでみたいという
意味です・・・
0153名も無きマテリアルさん
2006/07/24(月) 23:23:36Intel
Absoft
あたりじゃない？
0154名も無きマテリアルさん
2006/07/24(月) 23:39:44MS-Fortranってのは、Fortran90とは違うの
0155名も無きマテリアルさん
2006/08/08(火) 13:11:37どっちがいいんでしょう？
値段は高いが、GUI があって便利なのが CASTEPで、
使いにくいけど、安い(3000 EURO？)のが VASP、
って認識なんですけど、ほかに VASP と CASTEP の機能の
違いはありますか？
どなたかおっしゃってましたが、元は同じものだから機能
はそんなに変わらないんですかね？
0156名も無きマテリアルさん
2006/08/09(水) 04:04:10が必要なので高いです。
最終的な精度、将来性を見てもVASPが有利です。
CASTEPは全ての擬ポテンシャルが装備されているそうですが
かなり怪しいです。
0157名も無きマテリアルさん
2006/08/10(木) 10:20:40MedeAは、ぶっちゃけあまりちゃんとしてない･･･orz
0158名も無きマテリアルさん
2006/08/14(月) 19:06:550159名も無きマテリアルさん
2006/08/15(火) 14:56:01VASPは多いだろうなぁ。
とにかく、今後は全電子計算の意義も問われるだろうね。
VASP-PAWの様にfrozen-core All-electronは、
All-electronと言うはったりがきくし、
（結局の意味は擬ポテンシャルだからね）
次のVASPで予定されているrelaxed core
（SCFサイクルの中で擬ポテンシャルを作り直す事に対応）
がうまく行けばFLAPWと同じ程度の精度になって、
k点や平面波のカットオフを考えると最終的な精度は
VASP-PAWが勝かもしれない。
これで、DFPTや動力学やCPで構造最適化ができれば文句無いしね。
今時、FLAPWだって最強じゃないよ。
lo入れなければ線形法のエラーが心配だし、
コアまで入れてるから数値上のエラーも大きい。
価電子のみ考えればいい擬ポテンシャルと違って
何かと問題も多いしね。コード開発だって計算だって時間がかかる。
実際の応用でFLAPWでないといけない精度を追う状況ってそうは無いと思うしな。
いやいや、VASPの回し者じゃないですよ。。
ま～、日本でPAWの本質を理解している人が
どれだけいるのか疑問だけどね。
0160名も無きマテリアルさん
2006/08/15(火) 18:06:49沢山いるけど、その多くは昔から作ってる人達。
今からひとりで作るのはちょっと無謀じゃないかなー。
オリジナリティーが無いと研究にならないし、
これからの研究者は任期制とかで時間もないし。
0161名も無きマテリアルさん
2006/08/15(火) 20:41:14作れるとしてもノルム保存でしょ。
あれなら学生レベルでも可能。
これでカーパリネロぐらいまでなら趣味だよねぇ。
そこから先のいろんな応用とか物理量とか、その中でバグ取りとか。。
abinitとかのプロジェクトに参加していくのが普通だと思う。
ただ、使ってる人はある意味幸せよ。
バグがあったって、ソフトのせいに出来るし
新しいプログラムは人が作ってくれるんだから。
リターンキーとデータ整理の達人こそ今生き残ってるんだよね。
うらやましいなぁ、多分一番簡単な研究方法だしね。
0162名も無きマテリアルさん
2006/08/25(金) 17:01:160163名も無きマテリアルさん
2006/08/29(火) 04:06:550164名も無きマテリアルさん
2006/08/31(木) 22:17:04ただ既製ソフトを使ってるだけの
計算屋っているの?
0165名も無きマテリアルさん
2006/09/01(金) 00:01:570166名も無きマテリアルさん
2006/09/02(土) 00:23:33他人が作ったコードを大学や研究機関所有の
スパコンやクラスタで動かすだけなら
実験屋で十分じゃん
むしろ実験屋が使った方が有益
0167名も無きマテリアルさん
2006/09/03(日) 10:58:48実験装置を作る人も入れば、どっかから買ってきて
実験して結果を出す人もいるし。
実験屋ほど、理論に弱いと計算は出来ないし、
理論ばっかりに強くてもだめ。
バランスよくないといけないから、
実験屋が使えばいいとは言えない。
それこそ、間違った使い方をするよ。
第一、全部自分で計算コードを一から作ってる人は
こんな所に来ないと思う。
なんかのコードを参考にしたり、一部使ったり
もしくはabinitなんかに参加してるだけでしょ。
全部作るって言う事はかなり大変だよ。
どうでもいいことでも、実装するのは難しい場合が多いしね。
0168名も無きマテリアルさん
2006/09/03(日) 17:25:55手法の改良やコードの開発に寄与せずに
コードを使用するだけの計算屋てのはちょっとどうかと思う…
0169名も無きマテリアルさん
2006/09/03(日) 22:10:20それはそれで大変だと思います。
手法の改良はかなり大変です。
やっぱりしっかり基礎理論や計算手法を
理解していないといけないですし、
その後に実際にコードを開発して計算しないといけませんから
企業なんかでは厳しいでしょう。
コード開発といっても色々あるとは思いますが
単なるデータ整理レベルのものなら可能だと思います。
自分で１から開発は無理でしょうし、abinit系のプロジェクトに入るとしても
簡単な事はやられていますから。。
アイデアをもって計算する事も重要だと思いますよ。
基礎でなく応用と考えればそれはそれで十分価値があると思います。
個人的にはfortranとかいう最悪な言語のソースを読みたくないし、
書き加えたくないと思うこともあるので、
コード開発にはあんまり興味なしです。
0170名も無きマテリアルさん
2006/09/04(月) 08:14:33そうか。
それでもやっぱりあなたの書き込みを見ても
167みたいな計算屋は怠慢だと思う。
企業なら事情はわからんでも無いが…
0171名も無きマテリアルさん
2006/09/04(月) 08:24:03アイディアをちゃんとアルゴリズムとして書き下して
コード書ける人に任せりゃ良いじゃん
コード書けるけどアイディアが無い人は結構いるよ
一から書かなきゃ行けないほどのコードでも
アイディアを掲げて人を集めるか
開発できるところに持ち込めば、
価値があるなら何とかなる。
何にしたって一人で開発なんて無理なんだから
0172170
2006/09/04(月) 08:27:02167ではなく168です
0173名も無きマテリアルさん
2006/09/04(月) 10:47:04「素人は手を出すな」ではなく、「どんどん使ってください」という
オープンな態度でいるほうが、分野の発展のためにも良いと思う。
ただし、コードをいじってる人はそれだけ大変な仕事をしてるわけで、
（しかも、コードをいじったこと自体は物理の仕事としてあまり認められないし）
コードを使うだけの研究の場合は、質・量的に優れたものであってほしいとは思う。
0174169
2006/09/04(月) 12:01:39少し前なら重宝されてたし、今でも知ってる人の中では
非常に評価は高いと思います。
ただ、VASPみたいに安くて高性能なソフトが増えると
いままで地道に大変な作業をしてソフトを作ってきた人の
ソフトよりいいものが、簡単に手に入り、
それを使って計算した方が学会発表、論文発表、結果という点で
いつも有利に立つ事は事実ですからね。
実際にこの分野に流入してきている理論、実験家は
そういうところがおいしいと思ってきてる訳ですから、
その人たちが大多数になってきて、
実際に作ることへの評価はどんどん下がってきていると思います。
一杯いますよ、コードの中身や理論も中途半端に理解してるだけなのに
計算結果出して素人だましてるお偉い先生も。
>>171
fortranをいじりたくはないですが、
アイデアや式を書いたらコードを作ってくれる様な
偉い身分にはありません。。
自分で作れって感じですね。
やりますよ、いざとなったら。
そんなに難しい事じゃないので。大変は大変ですけどね。
0175名も無きマテリアルさん
2006/09/04(月) 13:45:28コード開発の苦労は誰も否定してない。
そこばっかり強調してもしょうがないだろ
それに、既にある優秀なコードが簡単に手に入るといっても
どっちにしたって、現存のモノに比べて新奇性、優位性がなければ
つくったって意味は無いんだから。
スペシャリストならそれなりの仕事をしろって言ってるんだよ
0176名も無きマテリアルさん
2006/09/04(月) 14:02:13どうでも良い屁理屈ばっかりで言い訳にもなってないな
0177169
2006/09/04(月) 16:44:12評価されるべきではあると思います。例えばコードを作る方なら、中身が分かってるので
ちょっといじるだけで、いろんな事を引き出せる自由度があります。
計算のプロにしても、計算のアイデアでいいものを出していけば良いと思います。
ただ、表向きに見てどうなんでしょうか？もしよろしければ、その点についてご意見いただきたいです。
コード開発のスペシャリストの仕事というのは手法開発や上記の様な点でしょうか？それは、マニアックな評価は受けますが
バンド計算をやりたいと思っている方でそんなことを求めている方はごく少数でしょう。
私は結局、物理量などを色々出していく方が一般的評価は得られますし、
いやらしい話お金も呼び込めると思います。現在の優秀なソフトはコード開発の先生方がつくってきたものより
優秀で様々な機能がついていたり、基本部分もしっかりしています。
それだけ、人の労力がかかっているからです。その優秀なソフトを売って、お金を稼いで、
また人を雇うという良いフィードバックがかかります。
すると、そのソフトを買って計算した方が新奇で優位な計算が出来る事になるのでは無いでしょうか？
日本以外ではソフト開発なんぞのんびりできる環境ではないとも聞きます。それほど、結果が求められているのだと思います。
日本でも、コードでスペシャリストとして食える人はどんどん減ってしかるべきだと考えるのですが、いかがでしょうか？
結果だけみれば、肩身は狭いですよ。裾野を広げてみたら、肩身が狭くなったって訳ですね。
やはり、コードは買って（もしくは少しは自前でいじれるものを持つ）計算のプロとしてもやっていける時代ではないですか？
どうやって計算するか？も重要ですが、何を計算するか？も重要だと思います。
色々、御託を並べてすみません。。。
0178名も無きマテリアルさん
2006/09/04(月) 17:17:00予算なんてつくの？
お金なんて要らないじゃん
0179名も無きマテリアルさん
2006/09/04(月) 19:46:24なぜコードを使うだけという怠慢をそんなに肯定しようとする！
よくわかってない人がいるのを良いことに甘えてやしないか？
そもそも、バンド計算しただけで価値ある仕事って何よ
0180169
2006/09/04(月) 20:42:16新しい超伝導物質は計算しただけでも論文になりますねぇ。
しかも、みんな出来るから競争ですよね。
誰でもある程度計算が出来るから、何を計算するかを
みんなが狙っているんではないでしょうか？
昔はバンド計算する事もすごかった訳ですから
それ自身重要な仕事だとは思います。
当然、バンド構造に対する議論は必要ですが。
いや、私もコードだけを使う事は反対です。もちろん。
みんながコードを作って、それで競い合う事が最もうれしい事だと思います。
ただ、現実、作ってる人はごく少数派。
コード作ってなくても大きな予算もついているようですよ。
計算機を買わないといけないですし、電気代もありますからね。
（そこらへんの議論は危ないのでやめましょう。）
計算を出来て、内容の議論が出来るっていうのは
必ずしもコードを全て把握している必要は無いと思います。
私も嫌ですよ、そりゃ。
コード書く事は楽しい事ですから。
（fortranは嫌いですけどね。）
バンド計算はコードを書いて計算するというのは
今となっては非主流派で、それを否定せず、
我々も移行していかなければならないと私は感じますが。。
実は書き込みしてる人も、知り合いだったりして。。。
まあ、多分そうでしょうけど。
0181169
2006/09/04(月) 20:50:35数的な意味で、決して立場的意味ではありません。
分野ではコードを作る人の評価はすばらしく高いと思います。
0182名も無きマテリアルさん
2006/09/04(月) 22:15:38書き込みから察するに、あなたもわかってるようだが、
あなたが挙げた例はバンド計算屋がやらなくたって良い仕事ですよね。
優秀なコードが手軽に使えるようになった今、
計算を齧った実験屋で十分にできてしまう。
あなたの書き込みから、
計算屋自らコードを書くような風潮では
なくなってきていることはよくわかった。
でもそれに迎合しなくても良いだろう。
物質に関する計算科学については要求が大きく
具体的になってきているから
進化は続きはするだろうが、
それに寄与せず、極端に言えば誰でもできる
ような成果を垂れ流しながら
この業界にぶら下がられるのは迷惑だよ。
失礼だが
あんたが泥臭いことしたくないだけだろ
0183名も無きマテリアルさん
2006/09/07(木) 23:20:410184名も無きマテリアルさん
2006/09/08(金) 14:03:33無駄に予算食いつぶしてるのかー
告発とかされたらどうなるんだろ
0185名も無きマテリアルさん
2006/09/08(金) 18:32:59優遇されてるところも、場所によってはある。
0186名も無きマテリアルさん
2006/09/08(金) 18:58:38と思っていたら、巨大な計算機の場所、空調、騒音、電気代が
かかることが問題になってきてるんじゃないですか？
0187名も無きマテリアルさん
2006/09/08(金) 21:16:370188優太７１５
2006/09/11(月) 18:57:43ベースとか、安いのありませんか？
0189名も無きマテリアルさん
2006/09/11(月) 19:11:47ガウシアンでも平面波でも何でもあるよ。
0190名も無きマテリアルさん
2006/09/15(金) 09:52:01ええーっと、ツッコミが欲しいのか？
0191名も無きマテリアルさん
2006/09/15(金) 13:10:28多分もういないよ
0192名も無きマテリアルさん
2006/09/16(土) 13:54:19Wien2k俺には使いこなせんorz
lapw1の設定難しすぎるよBlaha先生…
0193優太７１５
2006/09/17(日) 15:38:310194名も無きマテリアルさん
2006/09/23(土) 11:39:28欲しい物理量が実験値に合うけど
格子定数とか原子位置(結合距離)
なんかが実験値と対応しない場合
どう考えれば良いんでしょう。
少し縮む位なら温度の影響とかかと思うんだけど
格子が膨らんでたりすると…
なんだかいかにも無理矢理突っ込んだポテンシャルのせいで
膨らんでる気がする…。
0195名も無きマテリアルさん
2006/09/24(日) 16:40:030196194
2006/09/24(日) 20:28:05格子定数が1%位大きくなります。
Ue=0だと3桁は一致する程度。
しかし欲しい物理量は明らかに違う
やっぱりGGAのみだと合わないけど
だからってオンサイトのポテンシャルでは
現実のモデルにあってないってことですかね…
LDA+U使ってる論文とか読んでも
これでこれで計算しました。
いろいろポストプロセッシングして合わないところもいっぱいあるけど
この位のUが良いっぽいです。とか
Uを入れて計算したらこんなに挙動がかわって
変になったところもあるけど注目してるところは改善されてます。
とかそんなんばっかりorz
ぶっちゃけ欲しい物理量さえ再現できてれば
もうpublishしちゃっていいって感じなの？
曖昧にしか言及できない結果ばかりじゃないですか？
0197名も無きマテリアルさん
2006/09/25(月) 12:51:15を参照していただくと、結構経験的な事が書いてあります。
おそらく、U=0で格子定数がそのオーダーで合う事は
何か計算がまずいと考えてください。
擬ポテンシャルならばポテンシャルの質、収束方法の問題、
平面波などの量。
おそらく、ベースになる構造最適化の時点でまずい事が起こっている気がします。
というわけで、LDA+Uの問題ではないのでは？
結構、計算結果はそれっぽく出てきても、根本的な間違いも結構ありますよ。
0198194
2006/09/26(火) 14:35:13なるほど
参考になります。
GGAは経験的に少し格子定数を過大評価してしまう
(つまりストレスを過大評価？)
から、現実の格子定数と近すぎるのは逆におかしい訳ですね。
全電子のコード使ってるんで計算のミスがあったってのは結構辛いです。
カットオフやk点数はそれなりに追い込んだはずなんですが…
トータルエネルギーの格子体積依存性もきれいにでてたのに…
やっぱりデフォルトではゴーストバンドに捕まるからと
いろいろいじった当たりに問題があるのか…
もうこうなったら素人の学生がどうこうできる問題ではなくなってきてる
ような気がする…むしろ逃げたい…
やっぱり先生も含めてまわりに素人しかいない状況で
たった一人、独学で進めるのは無理がありますね
web上には紹介いただいた情報含めいろいろありますけど
専門家レベルの相談役の人がいないと辛すぎます…orz
そもそもこんな構造も複雑で3d遷移金属も入ってる局在系で
何をやりたかったんだろう…
調べれば調べるほどネガティブな情報しか出てこない
先生の思いつきにまともにつきあいすぎたか…
0199名も無きマテリアルさん
2006/09/26(火) 22:00:30全電子計算法だとすると、ストレスは簡単に計算出来ませんよね？
全エネルギーを生に計算して有限差分から出すしかありません。
そこに、原子内座標も同時に最適化するとなると、
計算は結構大変だと思いますが、その辺りで精度はどうだったのでしょうか？
おそらく使われているコードは wien2kですね？
LDA+Uが出来て、forceが計算出来て、全電子計算となると
ばれてしまうかもしれません。それか、FP-LOなどですかね？
wien2kであれば、LOをきっちりすれば
ゴーストの心配は少なくなるでしょう。
k点や、カットオフもよっぽどケチってない限り、
まずい事は無いと思います。
FLAPWなら、当然ながら、平面波のカットオフは、
平面波＋擬ポテンシャルより少なく出来ます。（USPPとかは知らん。）
まず、状態密度などは計算されましたか？
部分状態密度、バンド構造をみて、おかしな所はないかを見てください。
その上で、物理量がおかしければ、見ている物理量が
LDA,GGAではあまりうまく行かないものなのかもしれません。
第一原理って言ったって、何でも魔法みたいにうまくいくとは限りません。
もし良ければ、見ている物理量は何ですか？
教授も困り者です。
そんな、ソフトだからって簡単に出来ないですよねぇ。
0200194
2006/09/27(水) 11:26:50ご明察です。
構造最適化は格子ごとにジオメトリの
最適化も行うことになるので
膨大な時間がかかります。
k点一個ごとにCPU割り振っても通常のSCFサイクルが
収束するのに
半日かかる系ですから…(コードはBenchmarkの性能程度には最適化してあるにもかかわらず)
はじめはこれにさらに不純物とか欠陥入れたりも考えてました…
今考えるととてもやる気にはなれません
計算してる系は
3d遷移金属を含むリン酸塩で
LOは基本的にsemi coreのみくわえて
酸素だけ伝導帯のあたりにLO追加してます(defaultの設定)
Pの2sがMT球からリークしやすいのでその辺りをいじったりしてます。
scfのログをチェックしてみたところ
計算自体は順調おこなわれているようです。
状態密度は計算してみて部分状態密度をまとめて
教授を喜ばせていたりしてたんですが、
おかしい所というのが良くわかっていません。
これは私自身があるべき状態密度を
想像できないのが問題ですね。
バンド構造の図に至っては解釈のしかたすら勉強不足で
理解できていません…orz
0201194
2006/09/27(水) 11:29:03疑似生成エンタルピーとK-edge,L-edgeのスペクトルです。
核にかかるforceをつかってphonon計算
(PHONONは持っていないのでfrophoをなんとか使って)
もやりたいとかいろいろと考えていましたが
結局欲しいのは実験のスペクトルの再現(ケミカルシフトやL2,3の相対強度等)、
構造の安定性とそのエネルギーです。
0202名も無きマテリアルさん
2006/09/27(水) 13:09:09私が、要求しすぎましたね。すいません。
＞Pの2sがMT球からリークしやすいのでその辺りをいじったりしてます。
？？？MT球の大きさを変えたという事ですか？
Coreにせず、価電子にしたということですか？
また、MT球が重なっている事は無いでしょうか？
＞見ている物理量は基本的にはtotal nergyから見積もれる
＞疑似生成エンタルピーとK-edge,L-edgeのスペクトルです。
なるほど。書きすぎると狙いがばてれしまいますよ。
エンタルピーは、結構エラーが含まれてるのかもしれませんね。
細かい手法は色々あるとは思いますが、原子位置の最適化と
格子の最適化から出てくる全エネルギ－をフィッティングして、
出すものなので、例えば、フィッティングした線は
図にした時にキレイになっていますか？
遷移金属元素のedgeですね。
スペクトルは、ご存知の通り、スペクトル形状は大まかには比較出来ますが、
細かい形状までは比較は難しいです。所詮、一電子状態密度の差みたいなもんなので。
ただ、積分量としては比較しやすいと思います。（形状は関係ないので）
wienでどうやってるか知りませんが、L2,L3の相対強度を出す場合に、
どの部分まで、L2として、L3として扱っているんですかね？
この様な量だと、数%オーダーの実験との精度一致要求は厳しいかもしれませんね。
まあ、あまり力にはなれないかもしれませんが。。
0203名も無きマテリアルさん
2006/09/29(金) 11:42:100204194
2006/09/29(金) 14:15:01やはり大まかな比較しかできないということで
満足するしか無いみたいですね
MT半径の設定はWien2kでは重ならない
最大の半径を自動で決定できるので、
動かしたい量に合わせてMT半径を縮めているので
大丈夫です。
Pにかんしてはちょっとしたノウハウがあるのでご勘弁を。
Wien2kではL2,L3の強度比とスプリットはマニュアルで設定します。
実験のスペクトルを再現するために試行錯誤で
決めてたんですが…
まぁXASにかんしては元よりかなり情報がつぶれてしまっているので
スペクトル形状については大まかな再現ができれば問題なかったんですが、
問題はケミカルシフトの大きさでして…
やはり基底状態の平衡格子定数が気になるので
もう一度基礎から計算を見直してみます。
丁寧な対応ありがとうございました。
0205名も無きマテリアルさん
2006/09/30(土) 02:18:49当方は、FLAPWの専門家ではないのですが、
L2,L3の強度比をマニュアルで設定する？というのはどういう事ですか？
双極子遷移を仮定すれば、一電子近似の範囲内で、
L2,L3のスペクトルは一意に決まるはずです。
基本的に遷移金属？のcoreの2p1/2と2p3/2の
スピン軌道相互作用などによる分裂によって、
L2,L3のスプリット幅は決まるはずですが？
違いますか？
wienはCOREもFull relativisticですよね？
強度比も同様にして求まると思いますが。。。
これが、マニュアルで設定とは？？？？
つまり、この量はLDA（その他諸々の交換相関）による
コーンシャム固有関数と固有値から、第一原理的に決まるはずです。
あと、スペクトルの計算で、細かい形状まであわせようとすると、
k点が結構必要です。セルフコンシステントな計算の後に、
一度、k点を増やして。。。という過程を入れた方がいいと思います。
いずれにせよ、マニュアルというのは？
基底状態の格子定数の話は、ちょっと擬ポテンシャルを仮定した話でした。
おそらく、全電子法のエラーの範囲内なら、大きな差異は出ないと思います。
擬ポテンシャルの質が悪すぎて、変なバンド構造が出てきているのかと心配しただけです。
0206名も無きマテリアルさん
2006/10/15(日) 20:47:27こないだESCAの束縛エネルギー評価しようとマグネの単体、酸化物、金属間化合物で比較したのだが、
Mg1sの固有値は微動だにせず悩んだ気がする。
とりあえずフェルミ準位からの差分見積もったけど、実験結果とは定性的一致もしなかったんで
DV-Xa使ってしのいだ。
0208名も無きマテリアルさん
2006/10/18(水) 14:15:56技術的に簡単だし。
Mg1sは殆ど変わらないのでは？コアですし。
バレンスがスピン分極したりすると、変化があるかも。
DV-Xa。。。。
まだまだ、現役なんですね。
はぁ～、計算終わらねーよ。
0209名も無きマテリアルさん
2006/10/18(水) 15:33:13マジすかって感じの説明もあるし．
1119ﾍﾟｰｼﾞの(A2.20)と(A2.21)の間に書いてる文とか…
う～む．下手すっと専攻科修了できんかもしれん．泣きてぇorz
0210194
2006/10/18(水) 16:21:52Feでも1sシフトした経験有り
0211名も無きマテリアルさん
2006/10/18(水) 20:19:50やめとき、KKRはマニアックだから情報も無いし。
簡単な系なら別だけど。
それに特化したプログラムで済むしね。
ASWとかLMTOの方がいいと思うよ。線形化しただけでも全然違う方法だから。
>>210
確かに。可能性大ですね。失礼しました。
Feは、Fe自身が磁性を持ちやすいですからね。
周りの環境に引っ張られるかとオモいまし。
0212名も無きマテリアルさん
2006/10/18(水) 20:47:21それ、Kohn&Rostokerの奴でしょ。
KKRを勉強する上で読まなきゃならない原典ではあるけれど、
実は、それを読んでもプログラムを書くどころか、既存のプログラムの
内容を理解するのにも、あまり役にたたないと思われ。
だって1954年って、Hohenberg&KohnとかKohn&Shamの論文より前なわけだし。。。
0213209
2006/10/18(水) 21:51:47なるほど．Wikipediaの説明も無いに等しいのはそれが理由ですか．
この論文の次にLMTOの論文読む予定でした．
ﾊﾞﾝﾄﾞ計算についてはAPWしか知らない素人なもんで…
>>212
それですそれ．
論文の中であっちこっち飛べって書いてて錯綜してます．
しかし最近の人はあんまり自分でﾌﾟﾛｸﾞﾗﾑ書いたりしないのかなぁ．
いやぁそれにしてもKohn先生はﾀﾀﾞﾓﾝじゃないな．
0214206
2006/10/18(水) 21:52:59ただESCAでMgOの実験値は怪しいのしかなくて、実際測定したら炭酸Mg(と思われる)成分ばかりぐちゃぐちゃ出てきてねぇ…（涙
DV-Xa…すまぬいまだたまぁに使ってる（笑）
外に出す機会はあまり無いけれど、ちょっとした計算機実験して推定したりするのに使うのだ、お手軽だし定性的説明するなら役に立つし。
0215名も無きマテリアルさん
2006/10/18(水) 22:30:21いや、まだかなり初級レベルね。
O.K.Andersenの線形法の論文は嫌らしいよ。
最近の解説とかを呼んだ方が良いとおもう。
LMTOだったら、Andersenらの解説があるし、
藤原先生や小口先生の教科書、Richard Martinの教科書も良いと思います。
KKRは。。どうなんかなあ。
KKRとLMTOはもう全然別物だからね。
LMTOはスタンダードな固有値問題を解くけど、
KKRはグリーン関数を求める訳だし。
赤井先生の論文とか、FP-KKRの論文が参考になるかも。
ただ、重要なのは以下の二点。
１、論文を読んでも、式は途中までしかない。
論文に書いてある式を直接計算すれば良いってことは殆どないよ。
そこから、色々と進めないといけない。
ここは、自分の力だし、初めてやると添字が色々出てきたり
結構ごちゃごちゃする。これも時間がかかる事。
ここがしっかりしてないとプログラムも作れないよ。
２、実際のプログラムを作るにはノウハウがかなり重要。
色々細かなテクニックとか、常識とかあるので
それを一から作り上げるのは至難の業だと思う。
ただ、原子一つの、bcc、fccとかならまだ楽かも。
複数の原子でとか、構造が複雑で、、とか言い出すと
難しいと思うよ。頑張ってね。
ただ、今時KKRを作らないと卒業出来ないなんて、
目的が分からんなあ。
0216名も無きマテリアルさん
2006/10/18(水) 22:30:26Arスパッタしてみた？
0217名も無きマテリアルさん
2006/10/18(水) 22:33:33材料物性板のスレなのにリアルタイムにすすんでる…
0218名も無きマテリアルさん
2006/10/18(水) 22:34:30結構違いますね。参考になります。
難しい物質で、1eVの精度っていうと、、。比較も難しいですね。
DV-Xaでも、いいですよね。
結構目安とか、いきなり高尚な事やる前にって考えると。
しかも、結局解いてる式は同じですからね。
slater-exchangeかぁ。
まあ、考えてみると交換相関項の所だけ、いじってると思えば、
最近のハイブリッド系の交換相関の話にも似てますね。
原理は全然違いますが。。
0219名も無きマテリアルさん
2006/10/18(水) 22:35:52ごめん、KKRなら、赤井先生のノートがあるじゃん。
あれ読めば、殆ど分かるよ。
プログラムはそのだいぶ先だけど。
0220206
2006/10/19(木) 04:00:39実は試薬が粉末だった罠orz...
>>218
枯れている技術って安心できるから楽ですな（笑）
しかし、Fe1sのピークシフトか…XASといい(XAFSかな？)ESCAといい、SR施設に出入りしている不良が何人かいますねきっと？ぢつは私も今ビームラインでくぁｗせｄｆｒｔｇｙふじこｌｐ
0221名も無きマテリアルさん
2006/10/30(月) 18:28:110222名も無きマテリアルさん
2006/10/30(月) 18:36:54つgoogle
0223名も無きマテリアルさん
2006/10/30(月) 21:19:55ありがとうございます。
こんなサイト知りませんでした。
沢山出て来たので参考になりそうです。
0224名も無きマテリアルさん
2006/11/01(水) 20:04:10たとえばLegendre陪多項式や，Neumann関数の定義が(-1)^mだけ異なっているとか
はよく見るんですが，他にこのような定義の違いって沢山あるんですか？
今の場合は
Legendre陪多項式：Messiah，岩波公式集等とﾆｭｰﾒﾘｶﾙﾚｼﾋﾟ
Neumann関数：MessiahとSchiff
などなど．
Legendre陪多項式の定義が書いてないのに球面調和関数とLegendre陪多項式の関係で
(-1)^mだけずれてると不安になります・・・・
0225名も無きマテリアルさん
2006/11/02(木) 01:09:300226名も無きマテリアルさん
2006/11/02(木) 09:10:39なんか５０原子を超えると勝手に対角項しか吐き出さないので。。。非対角項を見たいのに。。。
0227名も無きマテリアルさん
2006/11/04(土) 14:51:27皆さん変数の精度はどうしてますか？
一応，倍精度で統一しているのですが，
実行速度に影響がなければ4倍精度で書いたほうが
いいのかなと思ってるんですが．．
0228名も無きマテリアルさん
2006/11/04(土) 15:45:430229名も無きマテリアルさん
2007/01/05(金) 01:33:48ＮａＣｌの構造を入力して、周期的境界条件も設定して計算したはずなのに、
計算結果に表示されるＨＯＭＯとＬＵＭＯの形がどうも・・・
その周期的境界条件を満たしてないのですが、左端と右端がつながらないつか。
こんなもんなんでしょうか。それとも
わたしゃ的外れなことをやってるのでしょーか。
そもそもＮａＣｌのＨＯＭＯとかＬＵＭＯとか言ってＯＫなのか？
0230名も無きマテリアルさん
2007/01/06(土) 02:57:410231名も無きマテリアルさん
2007/01/06(土) 14:42:03伝導体の下端をＬＵＭＯと呼んでいいのか
ということが一つだと思われ
0232名も無きマテリアルさん
2007/01/14(日) 01:30:520233名も無きマテリアルさん
2007/01/22(月) 00:18:240234名も無きマテリアルさん
2007/01/22(月) 21:04:550235バンド計算初心者
2007/01/22(月) 22:31:53実験値使う以外になんか方法ありますか
0236バンド計算初心者
2007/01/22(月) 22:33:190237名も無きマテリアルさん
2007/01/25(木) 21:57:41カットオフからすぐ出る。
>>235
マフィンティン半径は数値計算パラメータ。
原子も格子も変化させないなら、
限界まで大きくすりゃ良い。
>>236
長距離に置ける振る舞いを正しく
かつ、表式を簡単にしたのが
PBE(GGA96) PRLの論文読め。
0238名も無きマテリアルさん
2007/01/26(金) 01:11:56スットレ の検索結果 約 29 件
スレッド の検索結果 約 47,200,000 件
0239バンド計算初心者
2007/01/26(金) 19:27:100240バンド計算初心者
2007/01/26(金) 19:54:32カットオフからすぐ出るのはGmaxではなく
ｋmaxではないんですか？
またマフィンティン半径ですが
フルポテンシャル化した場合最大値ではなく
フルポテンシャルしていないものに比べ
少し小さい適当な値を取るとあるので困っています
ご鞭撻をお願いします。
0241名も無きマテリアルさん
2007/01/27(土) 15:27:25wienの事はよく知りません。
Gmaxってのは、chargeのカットオフの事ですね。
勘違いしておりました。
普通、バンド計算ではGmaxはkmaxの４倍で十分とされています。
平面波基底関数に対して電子密度は平面波＊平面波で勘定するからです。
ただ、ウルトラソフトの様にMT球的な物の中に
固いchargeをしこんでおく場合は１０倍くらいとります。
さらに、交換相関ポテンシャルが非線形なので４倍では不十分な事、
FLAPWでは、ワイナートの擬電子密度の部分にGmaxが関わるので
一概に４倍で十分とは言えません。
Gmaxはkmaxと別にとるべき収束パラメータです。
MT半径ですが、格子間領域をフルポテンシャルで扱わなくとも、
MT球内のみフルポテンシャルにする事は可能です。
おそらく、これをフルポテンシャル化してない場合と言っていますね？
MT球内は高々、LMAXと呼ばれる上限値をもうけたLAPW基底関数から
その範囲内で行列要素を評価しています。
非球対称成分は式を追っていくと、基底関数のLMAXまでの制限から
2LMAXまでの非球対称ポテンシャルしか評価しません。
基底ですので、まあLMAXに対して収束性をみればよいと思いますが、
格子間領域では非常に高いLに関するポテンシャルが存在する為に
MT球を大きくとってその成分を無視するよりは、
ワイナートの方法で平面波表現したポテンシャルから評価すべきだと
思ってるんじゃないですかねぇ？よくシランです。
そう考えると、少し小さめにするのがベストでしょう。
まあ、1~2Aの中で、自分の計算コストにあった値にすればいいんでは？
僕わ、LMTOなんで。
0242名も無きマテリアルさん
2007/02/01(木) 14:33:57でもバンド計算はポテンシャルをセルフコンシステントに決めるんですよね
ということは最初に用意する擬ポテンシャルってのは、初期状態だけを与えるものなんですか？
0243名も無きマテリアルさん
2007/02/01(木) 14:34:410244名も無きマテリアルさん
2007/02/01(木) 14:59:17コアのポテンシャルを置き換えるの。
そしたら、価電子のポテンシャルはもともと滑らかだから
滑らかな全体のポテンシャルで解ける。
普通の擬ポテンシャルはそのままコアポテンシャルは同じ。
価電子から来るポテンシャルは電子密度のセルフコンシステントな過程で変化する。
コアコレクションとかは自分で勉強しろ！
PAWやUSPPでは、もうちょっと価電子の変化による
ポテンシャルの変化を正確に扱える。
0245名も無きマテリアルさん
2007/02/01(木) 17:48:30ありがとうございます
通常はコアは同じものを使い続けるんですね
コア補正はそうではないという事ですかね…
イメージできましたので、ちょっと自分で調べてみます
0246名も無きマテリアルさん
2007/02/02(金) 16:03:00全電子密度が必要になるけど、どうしますかってこと。
擬ポテンシャル自身を改善するもんじゃないよ。
0247名も無きマテリアルさん
2007/02/02(金) 16:39:24内核の補正を行うのは擬ポテンシャルの生成時で、
内核と価電子帯との相互作用が大きいとき
内核補正をすると擬ポテンシャルの
トランスフェラビリティーがよくなる
…と理解していましたが、これは間違いですか？
0248名も無きマテリアルさん
2007/02/02(金) 17:30:580249名も無きマテリアルさん
2007/02/05(月) 12:40:520250名も無きマテリアルさん
2007/02/07(水) 07:51:53波数ｋの定義って
波動関数をL2空間の元と見て
exp（-ik・x）で平面波展開したときに出てきた
ｋのことなのか
Blochの定理から存在が示されるｋなのか
どっちなんですか？
というか一般的な定義ってなんなんでしょうか？
よくわかりません
0251名も無きマテリアルさん
2007/02/07(水) 17:35:07LDAにスピンの自由度を考慮したLSDAでは必ずスピン分極した解がでるのですか？
あと、LDA+UとLSDA+Uは同じですか？
0252名も無きマテリアルさん
2007/02/07(水) 17:48:26普通の定義ならブロッホの定理から来るもの。
>>251
必ず出る必要は無い。
LDA+UとLSDA+Uは厳密には違う。
非磁性で+Uするか磁性ありで+Uやるか。
式も少々違う。
0253251
2007/02/07(水) 18:21:31ありがとうございます。ついでにもう一つ質問があります。
実際にはスピン分極してない物質でもLSDA(+U)でスピン分極した解が出るのは、
LSDAの近似が悪いということでしょうか？
0254名も無きマテリアルさん
2007/02/07(水) 19:13:08使ってるソフトは？
0255251
2007/02/08(木) 11:16:05>>251にも書いたように、自分はバンド計算屋ではなく、バンド計算したことはありません。
あるスピン分極してない物質のLSDA、LSDA+Uのバンド計算の論文を見たら
スピン分極した解しか論文に載ってなかったので質問してみました。
論文を一通り見てみましたが、使ってるソフトが何かはわかりませんでした。
0256名も無きマテリアルさん
2007/02/08(木) 18:47:27小林先生のサイトにコメントが載ってるけど、
追加もしておきます。
まず、LSDAなどでスピン分極した解があるというのは
先生も言っている様に単に準安定状態の結果が出ただけで
本当に安定な解は非磁性の解という事もあります。
ただ、論文になっているのであれば、
（最近はソフトユーザーが多いので間違いもあるかもしれませんが）
スピン分極した結果と、非磁性解の全エネルギーを比較して、
分極した方がエネルギーが低いとなっているから
スピン分極のみの結果を出しているはずです。
この意味ではLSDA（もしくはLSDA+U）で安定な結果は
スピン分極した結果だという事です。
実際の実験などと合わない事は
１、LSDA,GGAなどの近似が悪い。
２、LSDA+UでUの値などがめちゃくちゃ。
３、コアコレクションなどの話もあるし、
擬ポテンシャルで磁性の議論は危険です。
VASPでは旨く行くと言う論文もありますが、
FLAPWなどの結果にあうように各物質で入念に作った
擬ポテンシャルを用いての結果なので
実は微妙な所は非常に危険です。
USPPやNCPPでは全くだめな場合も多いです。
（NCPPで平面波を増やしたからといって
PAWと同等にはならない）
４、計算した人が間違っている。
私の経験上、スピン分極しているものがスピン分極していない結果を第一原理計算が出して
それをGGAやLSDA+Uやるとスピン分極した正しい結果を生むと言う例は沢山あると思いますが、
スピン分局してないものを、スピン分極していると結果を出す事は、なかなかないと思います。
LSDA+Uの使い方の間違い、計算方法の問題、
もしくは実験結果が間違っている。だと思います。
0257名も無きマテリアルさん
2007/02/09(金) 11:51:41くわしい解説ありがとうございます。
私が言ってた物質の場合、実験の方は複数のグループで一致した結果なので間違ってないと思います。
となると、近似が悪いのか、計算が間違ってるのかどちらかということになりそうです。
0258名も無きマテリアルさん
2007/02/10(土) 16:49:250259名も無きマテリアルさん
2007/03/26(月) 00:31:24もってて計算できるね。最近つかってます。
0260名も無きマテリアルさん
2007/04/05(木) 00:46:36対称化するプログラムを書いてるんですが，困ったことがおきました．
球面調和関数を
(方位量子数，磁気量子数)
で表わしたとき，元々の基底は
(3,-3) (3,-2) (3,-1) (3,0) (3,1) (3,2) (3,3)
と取っています．つまり L = 3 の基底で基底関数の個数は7個です．
これをN点で対称化して立方調和関数を作ったところ，基底関数が11個できてしまいました．
基底関数が増えるのはおかしいですよね？
0261名も無きマテリアルさん
2007/04/05(木) 01:22:030262名も無きマテリアルさん
2007/05/05(土) 12:52:330263名も無きマテリアルさん
2007/05/05(土) 16:21:470264名も無きマテリアルさん
2007/05/05(土) 16:47:44ってことはDFPTか何かフォノンの計算も出来るの？
平面波なら簡単だけど、ああいう基底関数でもやってるのかね？
しかも結構まじめに計算できそうだし。
0265名も無きマテリアルさん
2007/05/05(土) 20:02:280266名も無きマテリアルさん
2007/05/08(火) 11:37:43擬ポテンシャルの間違いで誤った結果を、論文にしたみたいね。
PRBにあったよ。
Kresse＋Blahaがそれに噛み付いて、思いっきり恥かいてるね。
replyでは、間違いをみとめて、でも擬ポテンシャルを作り直したら
うまくいったって書いてあるけど、計算に合わせるように作っただけだろ～。
しかも、恥ずかしいからってLDA+Uの議論にすり替えてるし。。
ノルム保存やウルトラソフトで磁性とか、複雑な構造をやるとエラーが出るね。
ああ、怖いわ。
0267名も無きマテリアルさん
2007/05/08(火) 12:34:59競合ソフトを蹴落とそうとしてるみたいで…
0268名も無きマテリアルさん
2007/05/08(火) 18:34:06Kresseとしては、擬ポテンシャルで適当に計算されて
（PAWも似たようなもんだから）
ちょっとイライラしたんだろうなぁ。
競合ソフトではあるけど、開発力の面からしてVASPの勝ちかな。
pwscfもオープンだけど、開発協力者は居ないだろうし、
きちんと開発できるやつを集めて、クローズでやってるVASPの方が
まっとうな気がする。
方法論やるだけなら、一番簡単なノルム保存のabinitでやってしまえば良い訳だし。
てか、pwscfとか擬ポテンシャルで、磁性とかありえないし。
なんでもOKになってるよね。
0269名も無きマテリアルさん
2007/05/16(水) 18:54:11もちろんBSW記号付きで．
0270名も無きマテリアルさん
2007/08/24(金) 11:21:010271名も無きマテリアルさん
2007/08/25(土) 18:04:000272名も無きマテリアルさん
2007/08/27(月) 13:31:440273名も無きマテリアルさん
2007/09/24(月) 09:47:180274名も無きマテリアルさん
2007/09/27(木) 11:39:35orderN系だから、速いけど固体の細かいところは全然だめ。
プログラムソースもだめだね。
いろいろ計算できる物理量はあるけど、
基本的な性能に問題ありと思いますが。
ところで、wien2kがかなり高速化したみたいですね。
charge mixingと、iterative diagonalization.
RMM-DIISでは無い感じだけど、誰かリファレンス教えていただけませんか？
0275名も無きマテリアルさん
2007/09/29(土) 12:40:35mixerへのMulti-Secant法の採用による
サイクル数の大幅な低減が主。
構造最適化もサイクル数が減ったみたいだけど
使ってないから良くわからん。
Wien2kでの構造最適化は元が不器用で使いづらいから…
0276名も無きマテリアルさん
2007/10/01(月) 09:24:17構造最適化の方は、きっちり論文になってて
どうもいい感じみたいね。ロバストなようだし。
対角化のほうも誰かご存知でしたら教えてください。
0277名も無きマテリアルさん
2007/10/09(火) 11:21:310278名も無きマテリアルさん
2007/10/17(水) 20:43:30辛口コメント、本人が見たら泣くぞ。
0279名も無きマテリアルさん
2007/10/18(木) 23:54:320280名も無きマテリアルさん
2007/10/19(金) 00:07:26理由: 実用的な存在意義が無い
0281名も無きマテリアルさん
2007/10/19(金) 12:07:18OpenMXは、ユーザが増えてきちんといろんな系で
チェックをしていくうちに、だめな点が見つかって改善はされていくと思うが、
現在の感想は、まだ発展途上という感じかな。
OZAさんは、もうちょっと応用面でなく基礎面の強化に力を入れてほしいかな。
projectorとか、速さじゃなくて、基本的な精度。
>>279
こんな所にかける訳無いでしょ。
>>280
同意。ただ、日本にはその類いのプログラムは沢山あるけどね。
いまさら、ノルム保存で開発を始めてるとかいうグループもあるらしいし。
0282名も無きマテリアルさん
2007/10/25(木) 23:32:01とはいえ、すでにいろんなコードが手に入るからねぇ・・・
>>281
精度ってのは何の事？
0283名も無きマテリアルさん
2007/10/26(金) 01:41:23まぁ、O(N)だから仕方ないけど。。
この業界、0.1Aくらいずれても、平気で精度良いっていう業界だし。
色々あるにせよ。OpenMXのターゲットはLarge scaleだけじゃなくって
普通のバルクみたいな系でも使われてるじゃないですか。
そこで、チェックしてみたら。ってことですよ。
いろんなコードと比較してみたら、ヤバいかな～と感じました。
あと、NCPS2Kの話ですけど、例えば、自分がパソコンを作りたいと思ったときに、
真空管パソコンから作るかってことですよね。
ボードかって、そこから、自分のオリジナリティをつけるみたいな。
確かに、自分のコードって魅力的だけど、研究や開発には使えませんよね。
0284名も無きマテリアルさん
2007/10/28(日) 11:43:00てのがあるみたいだが
これってWien2kよりすごくないか？
0285名も無きマテリアルさん
2007/11/01(木) 22:03:14基本、全く同じなので好き勝手にdraxolがやってる感じか。
wienほど、ユーザーフレンドリーでないし、安定もしてないけど、
どんどんいろんな物理量を計算できるようにしてるってとこでしょう。
しかし、あの開発スピードはすごいと思いますね、FLAPWですからね、
式やコード開発は偽ポテンシャルとは比較にならないほどむずかしいですよ。
0286名も無きマテリアルさん
2007/11/01(木) 23:41:59いろんな研究で活用され始めて、
うまく結果を引き出して議論に組み込んでる研究者も増えてきた。
使えなさもそれと同時に大体、周知されてきたように思う。
となったら、EXCITING FP-LAPWみたいに、
出せる物性はとにかくいろいろユーザーが扱いやすい形で
出してやるコードが求められるだろう。
波動関数と固有値は与えてやるからあとは御勝手にってんじゃ
活用されんわな
0287名も無きマテリアルさん
2007/11/01(木) 23:49:49どこらへんが開発協力なんだ？
そうするつもりならGPLで開発なんかしないだろう
LAPW + local orbitalのフルポテンシャルコードだからってだけで
適当に言ってるんじゃないだろうな？
0288名も無きマテリアルさん
2007/11/02(金) 08:01:39ホームページ読んで。
あと、彼女はblaha一派。
0289名も無きマテリアルさん
2007/11/02(金) 08:07:37そうでしょうねぇ。今時、バンド構造計算できます！なんて
なんのウリにもならないし。
abinit,pwscf,VASPが人気なのはとにかく、
速くていろんな物理量が計算できるからだし。
EXCITINGは、案外みんなが欲しい物理量は計算できないかもしれないけどね～。
しかし、中国、アジアの研究者はひどいね。
新参者がabinitとかで、適当に計算して論文にしてるし。。
パソコンかってフリーのソフトでちょろちょろ計算かよ。。
やっぱり、コード開発しないんなら、計算屋とはいえないし、
もしくは、実験屋が色々評価するために実験＋計算ツールとして
使うんなら納得もいくけどねぇ。
0290名も無きマテリアルさん
2007/11/05(月) 00:18:46多くの場合、安価でそこそこ正しいUSPPとかPAWとかでことが足りる？
EXCITING FP-LAPWのメリットは何だろう。
0291名も無きマテリアルさん
2007/11/05(月) 19:34:12NMRとかXAS、EELS、Mossbauer、磁性とかじゃないか？
NMRは疑ポテンシャルで計算できるようにしてあったりするけど、
実際、ちゃんとした値が出てくるもんなのかね？
磁性がクリティカルなケースだとやっぱり駄目だと思うんだが…
0292名も無きマテリアルさん
2007/11/06(火) 23:11:24磁性の計算は一電子近時を使った計算手法では皆同程度ではなかろうか。
全電子計算が磁性に優れている理由は？
0293名も無きマテリアルさん
2007/11/07(水) 07:44:47それ以外にもやっぱり、精度の問題もあるよね。
PAWも精度はいいけど、あくまで偽ポテンシャル。
（frozen core all electronは真の意味でall electronじゃないから）
磁性もPAWでは結構、いいもんだと思うけど。
ポテンシャルの作り方でいくらでも結果はかわるし、
そこで心配になるのが普通でしょう。
もともと、いい偽ポテンシャルはFLAPWの結果と同じだという
センスがほとんど。ただし、FLAPWとLDAでのバンド構造とか、格子定数が一致したというだけなので、
そこから、GWとかをやって、うまくいくポテンシャルかどうかは不明。
その辺り、beyond LDA系にEXCITINGは強いですよ。
kresseもPAWポテンシャル作るときには、徹底的にwien2kとチェックをしてるみたい。
0294名も無きマテリアルさん
2007/12/17(月) 03:01:290295名も無きマテリアルさん
2007/12/26(水) 17:55:21バンド線図が書けるんだから使えるんじゃないの
http://www.openmx-square.org/openmx_man/node41.html
0296名も無きマテリアルさん
2007/12/28(金) 12:27:330297名も無きマテリアルさん
2008/02/03(日) 22:12:050298名も無きマテリアルさん
2008/02/04(月) 00:04:20バンド計算をするソフトじゃないよ
CASTEPのことを言いたいのか？
0299名も無きマテリアルさん
2008/02/06(水) 01:20:440300名も無きマテリアルさん
2008/02/08(金) 00:20:36All electron じゃないときちんと出ない理由っておしえて。
0301名も無きマテリアルさん
2008/02/16(土) 12:35:26基底関数が、スピン比非依存（基底関数を増やせばいい）
コアコレクションの問題（結局一番大きな問題）
0302名も無きマテリアルさん
2008/03/03(月) 02:42:58【上戸彩】アテンションプリーズ スペシャル [テレビドラマ]
0303名も無きマテリアルさん
2008/03/15(土) 23:56:51磁性の計算で疑似ポテンシャルでうまくいってない例教えて下さい。
0304名も無きマテリアルさん
2008/03/16(日) 01:54:57コアコレクションの効果が大きいというよりも、
コアに現実を十分に反映してないモデルを
使用することによって、精度が落ちる。
その結果として、磁性をはじめとする
比較的微妙なエネルギー差を議論する問題において
信頼できる結果を得るのが難しい。
結局のところ、要所要所でFLAPWなどの方法や詳細な実験データと
十分に比較検討する必要がある。
疑ポテンシャルだけで何でも言おうとしちゃ駄目だよってことだよ。
ソフトが吐き出す結果を鵜呑みにするのは論外だし、
どこまでの精度があり、どこまでの結論が導けるかを
十分にわかった上でないと、議論に使ってはいけません。
まぁ、結果を出すのは馬鹿みたいに簡単なんだし、それくらいは頭使え
0305名も無きマテリアルさん
2008/03/17(月) 10:27:58>>304
まぁ、コアを十分に反映していないけれど
そこそこ、いいってのがコアコレクションですね。
PAWとか、USPPは、コアはある程度無視しても、
SCFで偽ポテンシャルモドキが決まっていくので精度も良くなってる訳ですねーーーー。
うまくギポテンシャルをつくれば、なかなかよい結果は出てくるし。
ちなみに、FLAPWデ計算出来ない大きな系での
偽ポテンシャルとの比較はどうするの？
すでに偽ポテンシャルの構築はかなり精度が高くなってるし、
おおよそフローズンコアでまずい状況は減ってきてると思いますが？
FLAPWの方が数値的計算精度を計算負荷のために犠牲にしているから、
精度悪い場合が多いと思いますが～。
0306名も無きマテリアルさん
2008/03/17(月) 16:30:48そのうまく疑ポテンシャルを作ればってのが曲者
うまく出来てるかどうかはどうやって判断するの？
大きな系での疑ポテンシャルは
小さな系を切り出して考察すればよろし
0307名も無きマテリアルさん
2008/03/17(月) 17:35:06>精度悪い場合が多いと思いますが～。
具体例希望。
FLAPWの利点は高い精度だけ（言いすぎ？）なんだから、
だからこそFLAPWを使う人って精度評価に凄くシビアな印象があるんだが。
0308名も無きマテリアルさん
2008/03/17(月) 17:36:170309名も無きマテリアルさん
2008/03/18(火) 01:21:13そしたら、十分に小さな系でFLAPWと整合性をとった擬ポテンシャルは
十分に信頼出来る結果が出るでしょ？
VASPとかは、ポテンシャル構築段階でしっかりチェックされてるから
ユーザーとしては、細かい精度のチェックは基底関数の数くらいで十分なのでは？
そこまで、FLAPWジャナイトダメっていう時代じゃないと思う。
>>307
例えば、ブリルアンゾーンのkpointの数とか収束精度ね。対称性が下がると基底関数をケチったり、
少なめのkpointで計算していたりしますね。
擬ポテンシャルだと結構、しっかり平面波を持ってきて、kpointも多いし、
力とかの収束も十分にしてあるし～。どうでしょ？オレの気のせい？
0310名も無きマテリアルさん
2008/03/19(水) 19:34:19FLAPWじゃないとだめなんていってないが…
なにかコンプレックスでもあるのかね？
構築段階のチェックなんてどこまで信用できるかわからんし、
万能になってしまう程度に小さい系だと
不十分なケースがでてくるでしょう。
興味をもたれる物質ってのは特殊なものが多いしね。
後半で挙げてるのはユーザーの使い方の問題では？
コードや手法の本質的な問題とは程遠いｗ
0311名も無きマテリアルさん
2008/03/20(木) 05:46:06> 少なめのkpointで計算していたりしますね。
> 擬ポテンシャルだと結構、しっかり平面波を持ってきて、kpointも多いし、
> 力とかの収束も十分にしてあるし～。
むしろとりあえずこんだけやっとけば大丈夫だろって感じで
何も考えてなさそうなのが多いと思うけどな
0312名も無きマテリアルさん
2008/03/22(土) 20:52:12通常の一般的なLAPWなんかよりは
平面波の数をかなり少なく出来るし
何でもたくさん突っ込めば良いってもんでもないだろうに…
0313名も無きマテリアルさん
2008/03/29(土) 00:20:41表面上しってても、なんで平面波が少なくできて、
どんなプログラムを書かなきゃいけなくってとか。
では、質問です。私は分かりませんが、
LAPW+LOは平面波減らせないのですか？
0314名も無きマテリアルさん
2008/04/01(火) 10:23:060315名も無きマテリアルさん
2008/04/03(木) 08:48:17減らせません。論文読みなさい。
なんで、APWにしてるか考えなさい。
0316名も無きマテリアルさん
2008/04/13(日) 16:53:28何で減らせないの？
0317名も無きマテリアルさん
2008/04/18(金) 00:57:33apw+loのlocal orbitalとLAPW+LOのLOCAL ORBITALは、
本質的に別ものということを理解してると仮定して話すと。
apwをfixed energyで解いて+local orbitalを付け加えることで
エネルギー的に線形化しつつ、LAPWのようにfixed energy近傍の
解を正しく表現出来る。
これに加えて、マッフィンティン球表面で微分が連続と言う
拘束条件をapwでは課さないので、球内の動径関数の
変分自由度が同じ基底関数の量でも改良される。
0318名も無きマテリアルさん
2008/04/19(土) 15:57:12それはAPW+loの話では？
0319名も無きマテリアルさん
2008/04/20(日) 15:54:24> 変分自由度が同じ基底関数の量でも改良される
それはそうだが、
fixed energy近傍の線形化だけでなく、
単に、local orbitalが加わることによる効果ってのがあるわけなんだけど…
そこのところは考えてないのか？
0320名も無きマテリアルさん
2008/04/21(月) 12:30:03おれはFLAPWそこまで詳しくないが、
local orbitalが加わると確かに混合基底になるから、
平面波は減るけど、基底関数自身が減るかどうかは不明。
apw+loのlocal orbitalは、混合基底にして基底関数を減らす目的ではないはず。
なぜなら、LAPW法のように、マッフィンティンの中は
原子的な動径関数を持ってくれば、わざわざ混合基底にする必要なく
平面波を減らせるし、apw+loのlocal orbitalは、
apwで使った動径関数と同じようなもんだし。
混合基底の話にすれば、apwで使った動径関数以外の
変分自由度をもった関数を持ってきた方が基底関数は減るでしょう。
ところで、local orbitalが加わる効果ってなんでしょうか？
0321名も無きマテリアルさん
2008/04/21(月) 23:58:55自分が馬鹿だといわれない範囲で
腹の探り合いをしてる感じだな…
0322名も無きマテリアルさん
2008/04/22(火) 00:00:18基底関数は減らないってどういうこと？？？
0323名も無きマテリアルさん
2008/04/23(水) 13:45:34local orbitalが増えるというセンスに持っていくと、
平面波は減っても、基底関数が減るとは限らないってことでしょ。
local orbitalは原子ごとに作るから、原子が増えると一気に増えるからね。
0324名も無きマテリアルさん
2008/04/26(土) 00:28:150325名も無きマテリアルさん
2008/04/26(土) 16:45:350326名も無きマテリアルさん
2008/04/29(火) 09:36:19理由は？
0327名も無きマテリアルさん
2008/04/29(火) 22:14:590328名も無きマテリアルさん
2008/05/02(金) 01:37:44バンド計算することと、その材料の紫外可視吸収スペクトルを
シミュレーションすることは、同じことでしょうか？
0329名も無きマテリアルさん
2008/05/04(日) 16:52:45バンド計算の結果を使って
シミュレーションすることは出来る
0330名も無きマテリアルさん
2008/05/05(月) 18:27:18う～ん。減るかどうかは、かなり微妙だとは思うけど。
元々、LAPW+LOは基底関数を減らすために作られた訳じゃないから。
0331名も無きマテリアルさん
2008/05/06(火) 18:45:10波数ごとの複素誘電率が求まれば、紫外可視吸収スペクトルの
シミュレーションになると思うのですが、１対１で対応しないのは
膜厚の影響など、干渉などの効果のためでしょうか？
0332名も無きマテリアルさん
2008/05/06(火) 19:06:47のプロセスって、簡単ですか？？
0333名も無きマテリアルさん
2008/05/06(火) 22:09:43バンド計算が何をしていて、
紫外可視吸収スペクトルがどういった
物理量を反映しているのかはわかってる？
0334名も無きマテリアルさん
2008/05/06(火) 23:17:18ところまではわかってます。バンド計算で吸収端？というのかが
求まるところまではわかってる気がするのですが、そこから先が・・・
0335名も無きマテリアルさん
2008/05/07(水) 04:15:22そこから、電磁波の電場成分との相互作用を考えて
複素誘電率を計算しようと思えばできるが、
自分でかなりがんばる必要があるかもね。
光学特性を計算してくれるコードもあるみたいだけど
使ったことが無くて
関心も無いのでよくわからん
0336名も無きマテリアルさん
2008/05/07(水) 12:13:18結構、計算できる。
簡単かすると、状態密度にしか、依存しなくなる。
でも、あんまり細かい事をやろうとすると、
バンド計算の精度（励起状態）とか、吸収過程の細かい話
なんかでてきて、一致しずらくなるかもね。
wien2kとかでは、簡単にして計算してるよ。
0337328
2008/05/07(水) 23:27:04前からの疑問なのですが、バンドギャップが110nmのUVのエネルギーに
対応している場合って、110nm以上の波長では吸収なしで一応いいと
思うのですが･･･110nm以下の波長、つまりバンドギャップ以上のエネルギーの
吸収って、どう考えるべきなのでしょう・・・。
これって実は簡単じゃない話？
0338名も無きマテリアルさん
2008/05/08(木) 11:57:29110nmしか吸収しない。
量子力学的な離散エネルギー準位だから、
遷移もそのエネルギーでしか起こらない。
ただ、
いろんな過程とか、考えれば110nm以上でも以下でも吸収する。
温度とかもあるし、スペクトル自身は、急峻にはもちろんならないし、
ダブルピークとかも出てくる。
0339名も無きマテリアルさん
2008/05/10(土) 08:48:42このレスをみたあなたは・・・３日から７日に
ラッキーなことが起きるでしょう。片思いの人と両思いになったり
成績や順位が上ったりetc...でもこのレスをコピペして別々のスレに
５個貼り付けてください。貼り付けなかったら今あなたが１番起きて
ほしくないことが起きてしまうでしょう。
コピペするかしないかはあなた次第...
○。・。○。・。○。・。○。・。○。・。○。・。○。・。○。・。○
0340永遠の幻想 ◆UNIVERSEok
2008/07/06(日) 06:12:53【ﾉﾝｼﾞｬﾝﾙ】UNIVERSE PROJECT【ｵﾙﾀﾅﾃｨｳﾞ】
http://bubble6.2ch.net/test/read.cgi/band/1214494216/
0341名も無きマテリアルさん
2008/08/16(土) 20:36:04バンド計算を活用しようかと考えてるんだけど、
coreからの励起だとcore holeが重要になるケースも多いよね？
でも、Valenceでそういった検討をあまり見ない。
技術的に難しいのだろうか？
系にもよるんだろうけど…
0342名も無きマテリアルさん
2008/09/10(水) 23:27:290343名も無きマテリアルさん
2008/09/10(水) 23:45:330344名も無きマテリアルさん
2008/09/11(木) 18:11:38と言うことは一部の元素でしかできませんか？
0345名も無きマテリアルさん
2008/09/20(土) 16:58:27いろいろ調べたのですが・・・
誰か教えてください。
0346名も無きマテリアルさん
2008/09/20(土) 19:59:10使わないほうが良いよ
お金の無駄
0347345
2008/09/21(日) 07:00:37今度はもっと具体的に聞くので・・・
ttp://cms.mpi.univie.ac.at/vasp/
の
＞if you are interested in this package please contact Prof. Jurgen Hafner Juergen.Hafner_at_univie.ac.at.
に、使いたいよとのメールを送ればいいということでしょうか？
0348名も無きマテリアルさん
2008/10/11(土) 22:49:22バンド図描けるソフト
ハンド計算とは少し違うけど…もしよければどうぞ。
0349名も無きマテリアルさん
2008/10/31(金) 23:58:26金あるんなら業者にまかせりゃいいじゃん。
金無いんなら国産のソフト使っとけ。今じゃVASPに似たフリーソフトなんか山ほどある。
0350名も無きマテリアルさん
2008/11/07(金) 07:19:57VASPと似たフリーｿﾌﾄというのは、VASPと同じくらい精度があって、
色々計算出来てって言ってるの？
無いよ。しかも、国産ではまずあり得ない。
OpenMXとか？精度のレベルが違うし、計算出来ることのレベルも違う。
0351名も無きマテリアルさん
2008/11/11(火) 17:05:30小細工が次から次へと必要になるような仮説はたいてい根本から間違っているから、
御破算にして出発し直すべきである。
これに反して正しい仮説をうまく見つければ、あらゆる謎が次から次へと解けてくる。
南部陽一郎
0352名も無きマテリアルさん
2008/11/11(火) 17:07:45とりあえず
ABINITかPWSCFで良いんじゃね？
真似事ならできるし、金を払ってまで導入する価値が
有るかどうかもわかるようになるでしょう。
研究費は無駄に使うべきではない。
この業界の人間は大して研究費が必要でもないのに
とってきたりするから根本的に考え方が間違ってる奴が多い。
0353名も無きマテリアルさん
2009/01/31(土) 19:33:23他人の作ったコードで自動的にでてくる結果見てよろこんでるの？
科学に寄与してるの？
0354345
2009/02/07(土) 20:20:10abinitやpwscfってvaspよりも劣るの？
0355名も無きマテリアルさん
2009/02/11(水) 15:53:21いかがわしい、いや、コードを作らない計算だけする研究者はすぐに淘汰されるよ。
化学のGaussianよりマシでしょ。
というか、コードを一から作ることがいつも正しいとは言えないしね。
自作して、新しい結果を出せるようなコードを作るには、それだけ時間も必要だし。
科学への寄与を考えれば、abinitプロジェクトなんかに参加して、
今ある基盤を発展させることが重要だと思う。
>>354
vaspは、金がかかってるだけあって、
ポテンシャルの精度とかがよい。
pwscfは超ｿﾌﾄだし、結構危うい。
abinitは無料だけど、pawで何でも出来るって訳じゃない。
『第一原理計算』を実行するならどれも同じだけど、
結果の精度はかなり違うと思うな。
0356名も無きマテリアルさん
2009/02/11(水) 20:42:48決め手の精度ってのが主観なんですねｗ
VASPにしろ現行の擬ポテンシャル系コード(全電子でもそれほどでもないが)の
信頼性を考えて有意な差なんですかね？
0357名も無きマテリアルさん
2009/02/11(水) 20:45:58> 今ある基盤を発展させることが重要だと思う。
何故、今ある国内プロジェクトに寄与せず海外のプロジェクト？
理屈ばっかりこねて実践しない奴でも
国内のしょぼいプロジェクトなら活躍(笑)できるでしょｗ
0358名も無きマテリアルさん
2009/02/11(水) 21:01:22東北大の井上明久総長らの研究論文に再現性がないと指摘された問題で、
ナノ学会（会長・川添良幸東北大教授）は２２日、
同大２教授が出した質問書に回答するよう、井上総長らに電子メールで要請した。
問題は、井上総長らが９０年代に材料学の学術誌（英文）で発表した
金属ガラスの複数の論文に不正があるとする匿名の告発で浮上。
大学の対応委員会が０７年末、「事実無根」とする報告書を公表している。
これを不服とした東北大の大村泉・経済学研究科教授と
高橋禮二郎・国際文化研究科客員教授が昨年９月から今年１月２０日にかけて、
この学術誌を共同刊行している日本金属学会やナノ学会など１０団体に「質問書」を提出。
井上総長自身の説明を学会誌に掲載するよう求めていた。
ナノ学会・編集委員長の尾上順・東京工業大准教授は「言いがかりや誹謗中傷でない限り、
論文著者は質問に答える責任があると学会として判断した」としている。
朝日新聞宮城版３１頁（２００９年１月２３日）
井上総長の研究不正疑惑の解消を要望する会ホームページ
http://sites.google.com/site/httpwwwforumtohoku/Home
東北大総長が研究不正③
http://science6.2ch.net/test/read.cgi/rikei/1221116260/
0359名も無きマテリアルさん
2009/03/11(水) 13:52:450360345
2009/04/22(水) 11:37:250361名も無きマテリアルさん
2009/07/03(金) 00:40:15量子化学分野では、GaussianやGAMESSなど無料の入力デッキが普及して
大きく進歩しましたがバンド計算ソフトは
使いやすい入力デッキやバンド計算用空間群の無料ソフトが普及してないので
なかなか発展していないのが実情です。
0362名も無きマテリアルさん
2009/07/03(金) 21:42:16なあw
マフィンチン半径の実験値教えてや(笑)
判らんかったら原子半径で良いよ、ゼロ点振動程度は許すわ、でも原子外殻と原子内は厳密に定義付けして区別してくれよ。
もちろん実験屋だからこその絶大な信頼がおける、計算屋が裸足で逃げ出す精度でいっちょ頼むわ。
0363名も無きマテリアルさん
2009/07/08(水) 00:36:000364名も無きマテリアルさん
2009/07/26(日) 01:30:06わかってるNaCl型のNaClのバンド計算をやろうと思ったら、
ソフトのインストールから始めて何日くらいかかるでしょう？
キッテルレベルは理解してるつもりです。
0365名も無きマテリアルさん
2009/07/31(金) 20:05:000366名も無きマテリアルさん
2009/08/05(水) 01:33:46初心者向けの。
0367名も無きマテリアルさん
2009/08/15(土) 03:34:05キッテルみたいな迷著を使っているのか。
金があれば菱化ってところで買えばいいよ。
使い方も教えてくれる。数百万あれば。
0368名も無きマテリアルさん
2009/08/18(火) 14:48:46格子からの寄与の誘電テンソルと電子からの寄与の誘電率を求めたいのですが、
周りにシュミレーションを行っている人間がおらず、計算の仕方もさっぱりの状態です。
どのようにしたら計算できますか？
参考になる本等ありましたら、教えていただければ幸いです。
よろしくお願いします。
0369名も無きマテリアルさん
2009/08/18(火) 18:09:120370名も無きマテリアルさん
2009/08/18(火) 21:06:55講習会に出たいのは山々なのですが、参加は難しい状況で困っております。
非常に欲張っているというのは重々承知のつもりです。
何とぞよろしくお願いします。
0371名も無きマテリアルさん
2009/08/18(火) 22:51:11電子のヤツは、線形応答で出す。
CASTEPなら、こんな感じでは？
0372名も無きマテリアルさん
2009/08/20(木) 10:43:36phonon計算を行い、そこから固有振動数を出してもらう。
それを使ってLST式により導くという解釈でよろしいのでしょうか？
あと、線形応答で出すというのがいまいち分からないのですが
ご教授いただければ幸いです。
0373名も無きマテリアルさん
2009/08/21(金) 22:09:48Core i7 950 (3.06GHz)。FSB忘れた。
Intel Fortran, MKLつかって、
bench.HgがTotal CPU 138.521, Elapsed time 40.403だった。
みんなどのぐらい？
0374名も無きマテリアルさん
2009/08/22(土) 01:21:27第一原理計算を考えれば、それは早いオーダーだお。
0375名も無きマテリアルさん
2009/08/22(土) 01:22:32LSTって何？
linear responseは調べようぜ。
0376名も無きマテリアルさん
2009/08/23(日) 02:28:55Lyddane-Sacks-Teller式のことです。
線形応答は計算方法のことで、直接値として計算してもらうものではないと認識しております。
MS CASTEPで電子の寄与する誘電率を求める場合に直接値として出してもらうPropertyが線形応答！？
といった感じで混乱しております。
いろいろ試行錯誤した結果Optical Propertyで電子の寄与する誘電率が求められそうなのですが・・・
未熟者故ミスがあると思います。是非添削お願いいたします。
0377名も無きマテリアルさん
2009/08/23(日) 15:11:11あれを横軸、状態密度の図に変える計算て、
どんな計算でしょう？
ふと疑問をもったもんで。
0378名も無きマテリアルさん
2009/08/29(土) 13:59:54でもバンド図で見てるのは逆空間の一部分だけだから
本当はテトラヘドロン法とかを使って全空間を積分しなきゃいけないんだけどね
0379名も無きマテリアルさん
2009/08/29(土) 14:01:57フローズンフォノンは原子を直接動かしてモデル作るので大変
フォノンのバンド図をみるなら線形応答
フォノンの振動スペクトルを見るならフローズンフォノンでしょうか
0380名も無きマテリアルさん
2009/08/29(土) 14:03:17アブイニットとかオープンエムエックスとかだとフリーで
しかもマニュアルもインストール方法も充実してるよ
0381名も無きマテリアルさん
2009/08/29(土) 14:05:16あれは初期値を自動的に選択してくれるので初心者にはありがたいけど
最先端の研究に使用するならよほど気をつけないと駄目でしょうね
0382名も無きマテリアルさん
2009/08/29(土) 21:01:23それはCASTEPにこそあてはまる
最先端の研究に使われてるのは事実上
ほとんどがVASP
0383名も無きマテリアルさん
2009/08/30(日) 22:00:39バンドを波数で足し算って、線が横切ったところの波数の値を
足し算するってことそれとも、横切った回数をってこと？
0384名も無きマテリアルさん
2009/09/03(木) 15:02:14クレッセ一人 vs. 数十人の研究者って感じだね。
0385名も無きマテリアルさん
2009/09/03(木) 16:46:210386名も無きマテリアルさん
2009/10/03(土) 00:14:25シミュレーション種類は色々とあるけど、
どう考えて使い分ければいいの？
0387名も無きマテリアルさん
2009/10/05(月) 14:48:55横切った回数
0388名も無きマテリアルさん
2009/10/05(月) 14:49:56話題提供して
0389名も無きマテリアルさん
2009/10/06(火) 00:45:02>どう考えて使い分ければいいの？
実験にあうものは何かしら本質をついてると信じてそれを使う。
いや、マジでこの業界そんな感じよ
0390名も無きマテリアルさん
2009/10/21(水) 22:44:46elnesの計算てwien2k以外でもできるの？
0391名も無きマテリアルさん
2009/10/22(木) 09:28:08なんで出来ないと思うの？
0392名も無きマテリアルさん
2009/12/25(金) 20:59:240.49eV <- M.Shishkin and G.Kresse, Phys.Rev.B75, 2007, 235102
約0.3eV <- J.Hafne, J. Computational Chemistry, vol.29, 2008, 2044
0393名も無きマテリアルさん
2009/12/26(土) 10:39:190394名も無きマテリアルさん
2010/02/24(水) 01:50:120395名も無きマテリアルさん
2010/02/24(水) 04:24:39self-consistent GWとかの問題じゃないか？
バンドギャップはLDAとかGGAのDFTじゃまともにでないはず。
そこでGWでってことなんだろうが
GWのformalismで結果は変わるし
Zをかけたりしたらまた変わるし…
0396名も無きマテリアルさん
2010/04/27(火) 20:03:05第一原理やる意味ってある？
これで何かを開発できた例はあるのか？
0397名も無きマテリアルさん
2010/05/05(水) 13:58:58■ このスレッドは過去ログ倉庫に格納されています