 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 676 | 回復(fù): 0 | ||
[求助]
abinit 計(jì)算聲子譜無法運(yùn)行求助，急急急！
|
|
救急，請(qǐng)教各位大神幫忙，計(jì)算鋁的聲子譜，為什得不到DDB文件，下面是in文件幫忙修改一下，非常感謝！ # Crystalline al : computation of the phonon spectrum ndtset 18 #Set 1 : ground state self-consistency getwfk1 0 # Cancel default kptopt1 1 # Automatic generation of k points, taking # into account the symmetry nqpt1 0 # Cancel default tolvrs1 1.0d-18 # SCF stopping criterion (modify default) rfphon1 0 # Cancel default #Q vectors for all datasets #Complete set of symmetry-inequivalent qpt chosen to be commensurate # with kpt mesh so that only one set of GS wave functions is needed. #Generated automatically by running GS calculation with kptopt=1, # nshift=0, shiftk=0 0 0 (to include gamma) and taking output kpt set # file as qpt set. Set nstep=1 so only one iteration runs. nqpt 1 # One qpt for each dataset (only 0 or 1 allowed) # This is the default for all datasets and must # be explicitly turned off for dataset 1. qpt2 0.00000000E+00 0.00000000E+00 0.00000000E+00 qpt3 0.00000000E+00 0.00000000E+00 0.00000000E+00 qpt4 1.66666667E-01 0.00000000E+00 0.00000000E+00 qpt5 3.33333333E-01 0.00000000E+00 0.00000000E+00 qpt6 5.00000000E-01 0.00000000E+00 0.00000000E+00 qpt7 1.66666667E-01 1.66666667E-01 0.00000000E+00 qpt8 3.33333333E-01 1.66666667E-01 0.00000000E+00 qpt9 5.00000000E-01 1.66666667E-01 0.00000000E+00 qpt10 -3.33333333E-01 1.66666667E-01 0.00000000E+00 qpt11 -1.66666667E-01 1.66666667E-01 0.00000000E+00 qpt12 3.33333333E-01 3.33333333E-01 0.00000000E+00 qpt13 5.00000000E-01 3.33333333E-01 0.00000000E+00 qpt14 -3.33333333E-01 3.33333333E-01 0.00000000E+00 qpt15 5.00000000E-01 5.00000000E-01 0.00000000E+00 qpt16 5.00000000E-01 3.33333333E-01 1.66666667E-01 qpt17 -3.33333333E-01 3.33333333E-01 1.66666667E-01 qpt18 -3.33333333E-01 5.00000000E-01 1.66666667E-01 #Set 2 : Response function calculation of d/dk wave function iscf2 -3 # Need this non-self-consistent option for d/dk kptopt2 2 # Modify default to use time-reversal symmetry rfphon2 0 # Cancel default rfelfd2 2 # Calculate d/dk wave function only tolwfr2 1.0d-22 # Use wave function residual criterion instead #Set 3 : Response function calculation of Q=0 phonons and electric field pert. getddk3 2 # d/dk wave functions from last dataset kptopt3 2 # Modify default to use time-reversal symmetry rfelfd3 3 # Electric-field perturbation response only #Sets 4-10 : Finite-wave-vector phonon calculations (defaults for all datasets) getwfk 1 # Use GS wave functions from dataset1 kptopt 3 # Need full k-point set for finite-Q response rfphon 1 # Do phonon response rfatpol 1 1 # Treat displacements of all atoms rfdir 1 1 1 # Do all directions (symmetry will be used) tolvrs 1.0d-8 # This default is active for sets 3-10 ####################################################################### #Common input variables #Definition of the unit cell acell 3*7.50 # This is equivalent to 7.6 7.6 7.6 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 are two types of atom znucl 13 # 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, type 1 is the Aluminum, # type 2 is the Arsenic. #Definition of the atoms natom 1 # There are one atom typat 1 # The atom is of type 1 (Al) xred 0.0 0.0 0.0 #Gives the number of band, explicitely (do not take the default) nband 4 #Exchange-correlation functional ixc 1 # LDA Teter Pade parametrization #Definition of the planewave basis set ecut 6.0 # Maximal kinetic energy cut-off, in Hartree #Definition of the k-point grid ngkpt 8 8 8 nshiftk 4 # Use one copy of grid only (default) shiftk 0.0 0.0 0.5 # This gives the usual fcc Monkhorst-Pack grid 0.0 0.5 0.0 0.5 0.0 0.0 0.5 0.5 0.5 #Definition of the SCF procedure iscf 5 # Self-consistent calculation, using algorithm 5 nstep 25 # Maximal number of SCF cycles diemac 9.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. # The dielectric constant of AlAs is smaller that the one of Si (=12). # add to conserve old < 6.7.2 behavior for calculating forces at each SCF step optforces 1 #%%<BEGIN TEST_INFO> #%% [setup] #%% executable = abinit #%% test_chain = trf2_1.in, trf2_3.in, trf2_4.in, trf2_5.in, trf2_6.in, trf2_7.in #%% [files] #%% files_to_test = #%% trf2_1.out, tolnlines= 13, tolabs= 4.294e-05, tolrel= 5.000e-04, fld_options=-medium #%% psp_files = 13al.981214.fhi, 33as.pspnc #%% [paral_info] #%% max_nprocs = 1 #%% [extra_info] #%% author = #%% keywords = NC, DFPT #%% description = #%%<END TEST_INFO> |
回復(fù)此樓» 猜你喜歡
085405 考的11408求各位老師帶走
已經(jīng)有3人回復(fù)
-
339求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料與化工272求調(diào)劑
已經(jīng)有4人回復(fù)
-
085701環(huán)境工程，267求調(diào)劑
已經(jīng)有15人回復(fù)
-
材料求調(diào)劑一志愿哈工大324
已經(jīng)有9人回復(fù)
-
292求調(diào)劑
已經(jīng)有15人回復(fù)
-
304求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料與化工考研調(diào)劑
已經(jīng)有17人回復(fù)
-
321求調(diào)劑
已經(jīng)有5人回復(fù)
-
329求調(diào)劑
已經(jīng)有5人回復(fù)
找到一些相關(guān)的精華帖子，希望有用哦~
用Phonon計(jì)算fcc聲子譜不能選擇[110]方向？
已經(jīng)有4人回復(fù)
- ubuntu12.04用gfortran調(diào)用Lapack95函數(shù)庫出錯(cuò)。
已經(jīng)有3人回復(fù)
- openmpi安裝問題
已經(jīng)有17人回復(fù)
- Abinit計(jì)算聲子問題
已經(jīng)有10人回復(fù)
- 【求助】abinit聲子計(jì)算
已經(jīng)有7人回復(fù)
- 【轉(zhuǎn)帖】VASP和ABINIT優(yōu)缺點(diǎn)比較
已經(jīng)有18人回復(fù)
- 【討論】交流討論下用VASP計(jì)算不同溫度下材料的熱力學(xué)性能
已經(jīng)有19人回復(fù)
點(diǎn)擊這里搜索更多相關(guān)資源科研從小木蟲開始，人人為我，我為人人
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zbo007 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|---|---|---|---|---|
|
[考研]
 一志愿廈門大學(xué)化學(xué)學(xué)碩307求調(diào)劑
+9
|
y7czhao 2026-03-26 | 9/450 |
|
|
|
[考研] 286求調(diào)劑 +10 | PolarBear11 2026-03-26 | 10/500 |
|
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +6 | hanamiko 2026-03-27 | 6/300 |
|
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
|
[考研] 291求調(diào)劑 +7 | 孅華 2026-03-22 | 7/350 |
|
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
|
[考研] 330一志愿中國海洋大學(xué) 化學(xué)工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
|
[考研] 學(xué)碩274求調(diào)劑 +5 | Li李魚 2026-03-26 | 5/250 |
|
|
|
[考研] 一志愿上海理工能源動(dòng)力（085800）310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
|
[考研] 307求調(diào)劑 +8 | 超級(jí)伊昂大王 2026-03-24 | 9/450 |
|
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
|
[考研] 340求調(diào)劑 +3 | Amber00 2026-03-26 | 3/150 |
|
|
|
[考研] 一志愿211 初試270分 求調(diào)劑 +6 | 谷雨上岸 2026-03-23 | 7/350 |
|
|
|
[考研] 271求調(diào)劑 +6 | 生如夏花… 2026-03-22 | 6/300 |
|
|
|
[考研] 一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑 +9 | 稻妻小編 2026-03-26 | 9/450 |
|
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|
|
|
[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +5 | leaflight 2026-03-24 | 5/250 |
|
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 6/300 |
|