 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 760 | 回復(fù): 0 | |||
hhwj340銀蟲 (初入文壇)
|
[交流]
【求助】Ni的聲子計(jì)算
|
|
我想計(jì)算Ni的聲子譜。 第一步，先做了結(jié)構(gòu)優(yōu)化計(jì)算，輸入文件如下： # Ni fcc structure with two atoms per unit cell #Ni is a ferromagnatic material #spin definition nsppol 2 nspden 2 spinat 0.0 0.0 4.0 prtdos 1 prtden 1 #Structure optimize optcell 1 ionmov 2 ntime 15 dilatmx 1.05 tolmxf 1.0e-6 ecutsm 0.5 ####################################################### #Common data acell 3*6.65184 rprim 0.0 0.5 0.5 0.5 0.0 0.5 0.5 0.5 0.0 ntypat 1 znucl 28 natom 1 typat 1 xred 0.0 0.0 0.0 ixc 11 occopt 3 tsmear 0.01 nband 10 kptopt 1 ngkpt 8 8 8 nshiftk 4 shiftk 0.5 0.5 0.5 0.5 0.0 0.0 0.0 0.5 0.0 0.0 0.0 0.5 ecut 60.0 tolvrs 1.0d-18 iscf 5 nstep 70 diemac 12 計(jì)算得到的結(jié)果是：acell=6.73,查閱文獻(xiàn)的結(jié)果是：6.65. 下面計(jì)算聲子， 第一步，先做了，GS 計(jì)算,輸入文件是： # Ni fcc structure with two atoms per unit cell #Ni is a ferromagnatic material #spin definition nsppol 2 nspden 2 spinat 0.0 0.0 4.0 prtdos 1 prtden 1 #Structure optimize #optcell 1 #ionmov 2 #ntime 15 #dilatmx 1.05 #tolmxf 1.0e-6 ecutsm 0.5 ####################################################### #Common data acell 3*6.65 rprim 0.0 0.5 0.5 0.5 0.0 0.5 0.5 0.5 0.0 ntypat 1 znucl 28 natom 1 typat 1 xred 0.0 0.0 0.0 ixc 11 occopt 3 tsmear 0.01 nband 10 kptopt 1 ngkpt 8 8 8 nshiftk 4 shiftk 0.5 0.5 0.5 0.5 0.0 0.0 0.0 0.5 0.0 0.0 0.0 0.5 ecut 50 tolvrs 1.0d-18 iscf 5 nstep 70 diemac 12 計(jì)算得到WFK文件，以上兩步計(jì)算都沒有問題，no warning and no comment. 第二步計(jì)算Gamma點(diǎn)的聲子。 輸入文件如下， # Crystalline Ni : computation of the second drivative of the total energy # #Ni is a ferromagnatic material #spin definition nsppol 2 nspden 2 spinat 0.0 0.0 4.0 #Response-function calculatio, with q=0 rfphon 1 # Will consider phonon-type perturbation rfatpol 1 1 # all the atoms wiil be displaced rfdir 1 1 1 # Along the first reduced coordinate axis nqpt 1 # One wavevector is to be considered qpt 0 0 0 # This wavevector is q=0 (Gamma) kptopt 2 # Automatic generation of k points, taking # into account the time-reversal symmetry only tolvrs 1.0d-8 # SCF stopping criterion iscf 5 # Self-consistent calculation, using algorithm 5 irdwfk 1 # Read the ground-state wavefunctions ######################################################################### #common input variables #Definition of the unit cell acell 3*6.65 rprim 0.0 0.5 0.5 # In lessons 1 and 2, these primitive vectors 0.5 0.0 0.5 # (to be scaled by acell) were 1 0 0 0 1 0 0 0 1 0.5 0.5 0.0 # that is, the default. #Definition of the atom types ntypat 1 # There is only one type of atom znucl 28 # The keyword "znucl" refers to the atomic number of the # possible type(s) of atom. The pseudopotential(s) # mentioned in the "files" file must correspond # to the type(s) of atom. Here, the only type is Silicon. #Definition of the atoms natom 1 # There are two atoms typat 1 # They both are of type 1, that is, Silicon. xred # This keyword indicate that the location of the atoms # will follow, one triplet of number for each atom 0.0 0.0 0.0 # Triplet giving the REDUCED coordinate of atom 1. # Triplet giving the REDUCED coordinate of atom 2. # Note the use of fractions (remember the limited # interpreter capabilities of ABINIT) #Definition of the planewave basis set ecut 50.0 # Maximal kinetic energy cut-off, in Hartree #Definition of the k-point grid ngkpt 8 8 8 # This is a 8x8x8 grid based on the primitive vectors nshiftk 4 # of the reciprocal space (that form a BCC lattice !), # repeated four times, with different shifts : shiftk 0.5 0.5 0.5 0.5 0.0 0.0 0.0 0.5 0.0 0.0 0.0 0.5 # In cartesian coordinates, this grid is simple cubic, and # actually corresponds to the # so-called 4x4x4 Monkhorst-Pack grid ixc 11 occopt 3 tsmear 0.05 nband 10 #Definition of the SCF procedure nstep 80 # Maximal number of SCF cycles diemac 12.0 # Although this is not mandatory, it is worth to # precondition the SCF cycle. The model dielectric # function used as the standard preconditioner # is described in the "dielng" input variable section. # Here, we follow the prescription for bulk silicon. 問題就來了，得不到合理的結(jié)果， 而且log中有一大堆這樣的WARNING: cgwf3.F90:701:WARNING New trial energy at line 4 = -8.174248E+02 is higher than former: -8.174248E+02 得到Gamma點(diǎn)聲子竟然為： Phonon wavevector (reduced coordinates) : 0.00000 0.00000 0.00000 Phonon energies in Hartree : -3.212615E-03 -3.212615E-03 -3.212615E-03 Phonon frequencies in cm-1 : - -7.050876E+02 -7.050876E+02 -7.050876E+02 有哪位達(dá)人能夠幫我�。扛嬖V我問題出在哪里？ 計(jì)算中所使用的贗勢，可以在這里下載： http://www.sas.upenn.edu/rappegroup/htd ... p_gga.html |
回復(fù)此樓» 猜你喜歡
08工學(xué)調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有25人回復(fù)
-
085600材料與化工
已經(jīng)有6人回復(fù)
-
265求調(diào)劑
已經(jīng)有8人回復(fù)
-
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有21人回復(fù)
-
295材料求調(diào)劑，一志愿武漢理工085601專碩
已經(jīng)有4人回復(fù)
-
0856調(diào)劑，是學(xué)校就去
已經(jīng)有9人回復(fù)
-
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑
已經(jīng)有4人回復(fù)
-
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？
已經(jīng)有7人回復(fù)
找到一些相關(guān)的精華帖子，希望有用哦~
今天用CASTEP計(jì)算聲子譜，出現(xiàn)下列提示，有能看懂的大俠幫忙解釋下
已經(jīng)有11人回復(fù)
- 為何用晶格常數(shù)的實(shí)驗(yàn)值計(jì)算聲子？
已經(jīng)有18人回復(fù)
- 關(guān)于用超胞方法計(jì)算聲子的一些經(jīng)驗(yàn)
已經(jīng)有21人回復(fù)
- Abinit計(jì)算聲子問題
已經(jīng)有10人回復(fù)
- ZnO的聲子計(jì)算
已經(jīng)有18人回復(fù)
- 聲子計(jì)算又負(fù)頻的出現(xiàn)�。�！
已經(jīng)有24人回復(fù)
- MS計(jì)算晶體聲子譜的時候，橫坐標(biāo)有負(fù)能量值，為什么？
已經(jīng)有11人回復(fù)
- castep計(jì)算Fe2O3聲子譜失敗，原因求教~~
已經(jīng)有10人回復(fù)
- 自旋體系的聲子計(jì)算應(yīng)該用什么軟件呢？
已經(jīng)有5人回復(fù)
- 【求助】castep計(jì)算聲子譜，為什么不能進(jìn)行振動分析
已經(jīng)有9人回復(fù)
- 【求助】使用vasp進(jìn)行帶心凍結(jié)聲子計(jì)算的幾個問題~~~
已經(jīng)有9人回復(fù)
- 【求助】請教CASTEP計(jì)算聲子譜后出現(xiàn)虛頻
已經(jīng)有6人回復(fù)
- 【求助】有人用有限差分方法計(jì)算過G點(diǎn)聲子譜嗎？
已經(jīng)有10人回復(fù)
- 【求助】請問大家在做聲子計(jì)算時，這個 qpt是按照什么設(shè)置的？？？？
已經(jīng)有9人回復(fù)
- 【求助】用MS中的castep計(jì)算聲子振動模頻率，具體怎么算
已經(jīng)有10人回復(fù)
點(diǎn)擊這里搜索更多相關(guān)資源科研從小木蟲開始，人人為我，我為人人
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 hhwj340 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
295復(fù)試調(diào)劑
+7
|
簡木ChuFront 2026-03-19 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +4 | 西南交通專材355 2026-03-19 | 4/200 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +6 | z1z2z3879 2026-03-18 | 6/300 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
|
[考研] 一志愿天大材料與化工（085600）總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯！ 2026-03-15 | 18/900 |
|
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
|
[基金申請] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|