 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3822 | 回復: 2 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
xun梟水草新蟲 (初入文壇)
|
[交流]
[VASP] 求助：DFPT法算出來的有LO-TO劈裂的聲子譜在G點斷開了，要如何解決呢 已有2人參與
|
||
|
如題，希望能有聲子計算方面的高人幫我解決一下這個問題，我真的技窮了 壓縮包里的SPOSCAR是3*3*3擴胞的結果，是我計算所使用的。在phonopy進行后處理的時候額外考慮了LO-TO劈裂 INCAR： NELM=100 EDIFF = 1e-9 ! It is a recommended value. When calculating phonon, 1e-9 at least. AMIX=0.7 BMIX=0.01 #NBANDS = 96 ! Default = NELECT/2 + NIONS/2 #NPAR=8 #NCORE=6 KPAR=10 ! When irreducible k-points are large in amount, may increase KPAR, which will increase memory usage in most instances and sharply increase speed. ###Parallel settings are related to convergence, but have no certain relation. If electron step fail, change the above three may solve problem. ISTART = 0 ICHARG = 2 ! For band structures,with ISTART=0 and ICHARG=11. Copying WAVECAR and using ISTRAT=1 may lead to fault in high excited states and lower speed #LORBIT = 11 #LAMXMIX=4 ! system includes d-element, 4 ; system includes f-element, 6 LWAVE = .FALSE. LCHARG = .FALSE. ! Sometimes it should be commented #NWRITE=3 #LREAL=A ! When calculating phonon, do not change it. Default setting is .FALSE. #ISYM = -1 SYMPREC = 1e-6 ! Default 1e-5 #GGA = PE GGA_COMPAT = .FALSE. ! Apply spherical cutoff on gradient field. #ALGO = Fast ! When precision is high in many aspects, usually need to comment it. NEVER use Exact except some kinds of GW. Exact will decrease speed by thousands of times. When electron step goes wrong, see wiki or comment it to solve. ISMEAR = 0 ! Gaussian smearing usually, 1 for metal, -5 for DOS or accurate Etot SIGMA = 0.01 ! Default=0.2, for metals; other should be smaller. For phonon calculation, 0.03 at least PREC = Accurate ! Use it in phonon calculation ENCUT = 450 ! Please set it manually KSPACING=0.15 ADDGRID=.TRUE. #LSORBIT = .TRUE. #LNONCOLLINEAR = .TRUE. #ISPIN = 2 ! Default=1; 2 for spin-polarized calculations #MAGMOM = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 #LDAU = .TRUE. #LDAUU = [real array] ! Coulomb interaction parameters #LDAUJ = [real array] ! Exchange interaction parameters #LHFCALC = .TRUE. ! Hartree-Fock Calculation. Recommended all bands simultaneous algorithm #HFSCREEN = [real] #SCISSOR= #LVEL = .TRUE. #LOPTICS = .TRUE. CSHIFT=0.05 NEDOS=4000 ! Number of grid points for the electronic density of states and dielectric function #OMEGAMAX=70 ! It can not change the frequency range of calculation for ionic dielectricfunction whatever IBRION=6 or 7 or 8, but can change range for electron dielectricfunction and cause little(or nothing) difference to dielectricfunction's values when the change of OMEGAMAX is not large, so in this case, OMEGAMAX can be used. #For ionic dielectricfunction, it is possible for VASP that whole range(highest frequency) can be calculated will be calculated, so to change range, should change accuracy to include higher frequency phonon modes, which means change supercell size. K density, ISYM, NBAND, parallel setting have very little effects on upper limit #OMEGATL ! do not effect electron and ion dielectricfunction in every aspect, should not use it. #DIPOL = [real array] ! center for cell in direct coordinates #LVTOT = .TRUE. ! Total local potential #LVHAR = .TRUE. ! Containing entire local potential #LELF = .TRUE. ! Electron localization function #LBERRY = .TRUE. ! #LCALCPOL = .TRUE. #LCALCEPS = .TRUE. #EFIELD_PEAD = 0 0 0 ! unit eV/Ang #LASPH=.TRUE. ! For Raman's latter half series calculations, use it and LEPSILON. IBRION=8 LEPSILON=.TRUE. ! Ionic dielectric function. DFPT need no more empty bands #use supercell that: 1three lattice constants have little difference, 2lattice constants > 10A, 3atom number > 90. Every result relying on IBRION=5-8 must use supercell. #When ISYM=2, do not support NCORE!=1, but support kpoints parallel. Kpoints parallel do not change dielectricfunction's values. When ISYM=-1, support band parallel. This role also applies to Raman. #LPHON_POLAR=.TRUE. ! Static dielectric tensor and Born charges that are needed for LO-TO splitting are also need supercell calculation. #PHON_DIELECTRIC= #PHON_BORN_CHARGES= |
回復此樓» 猜你喜歡
308求調劑
已經有6人回復
-
材料調劑，307分
已經有3人回復
-
085600材料與化工調劑 280分
已經有15人回復
-
求材料讀博院校
已經有8人回復
-
申博
已經有8人回復
-
復試調劑
已經有3人回復
-
中科大材料299求調劑
已經有16人回復
-
296求調劑
已經有3人回復
-
A區(qū)一本交叉課題組，低分調劑，招收機械電子信息通信等交叉方向
已經有65人回復
-
一志愿中科大材料與化工。353分能調劑到哪？
已經有3人回復
小蟲蟲LG
木蟲 (正式寫手)
- 應助: 114 (高中生)
- 金幣: 1775.9
- 散金: 1871
- 紅花: 25
- 帖子: 975
- 在線: 678.5小時
- 蟲號: 2128430
- 注冊: 2012-11-15
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
3樓2025-05-16 16:53:23
dxcharlary
專家顧問 (職業(yè)作家)
-
專家經驗: +83 - 1ST強帖: 2
- 應助: 251 (大學生)
- 金幣: 10416.6
- 散金: 155
- 紅花: 65
- 帖子: 3137
- 在線: 757.4小時
- 蟲號: 936718
- 注冊: 2010-01-04
- 專業(yè): 凝聚態(tài)物性I:結構、力學和
- 管轄: 第一性原理
2樓2025-05-03 20:34:29
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿東北大學080500，英一數二343分求調劑 +7 | Seele_v5 2026-03-04 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 0703化學求調劑 +3 | 很老實人 2026-03-09 | 3/150 |
|
|
|
[考研] 0703化學求調劑，總分309分，一志愿華南師范 +3 | 花與葉@ 2026-03-08 | 3/150 |
|
|
|
[考研] 求調劑 +3 | 拾柒12。 2026-03-08 | 3/150 |
|
|
|
[考研] 一志愿中國石油大學（華東） 本科齊魯工業(yè)大學 求調劑 +3 | snw石 2026-03-07 | 3/150 |
|
|
|
[考研] 一志愿211 085600 280數二英二求調劑 +3 | 月山斜 2026-03-06 | 3/150 |
|
|
|
[考研] 085600材料與化工 298 調劑 +11 | 小西笑嘻嘻 2026-03-03 | 11/550 |
|
|
|
[考研] 求調劑 +4 | 呼呼？~+123456 2026-03-05 | 5/250 |
|
|
|
[考研] 268求調劑 +4 | 劉合華 2026-03-05 | 4/200 |
|
|
|
[考研] 材料考研339求調劑 +3 | Karry*^_^* 2026-03-04 | 3/150 |
|
|
|
[考研] 328求調劑 +6 | vuzhdkfjkx 2026-03-04 | 7/350 |
|
|
|
[考研] 282求調劑 +7 | 夕～日 2026-03-05 | 8/400 |
|
|
|
[考研] 調劑材料學碩 +6 | 詞凝Y 2026-03-02 | 6/300 |
|
|
|
[考研] 一志愿鄭州大學，學碩，物理化學， 333求調劑 +5 | 李魔女斗篷 2026-03-04 | 5/250 |
|
|
|
[考研] 312求調劑 +8 | 醋精華了一下發(fā)?/a> 2026-03-03 | 10/500 |
|
|
|
[考研] 一志愿西交大材料學碩 346 求調劑 +3 | zju51 2026-03-04 | 3/150 |
|
|
|
[論文投稿]
EST拒稿重投
 5+3
|
15102603076 2026-03-02 | 3/150 |
|
|
|
[考研] 290求調劑 +9 | ErMiao1020 2026-03-02 | 9/450 |
|
|
|
[考研] 278求調劑 +3 | 滿天星11_22 2026-03-02 | 3/150 |
|
|
|
[考研] 化工京區(qū)271求調劑 +7 | 11ing 2026-03-02 | 7/350 |
|