| 5 | 1/1 | 返回列表 |
| 查看: 1608 | 回復: 4 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
[求助]
如何計算clamped-ion 和free_stress 狀態(tài)下的介電常數(shù)
|
||
|
這是本人的.in文件,請問告訴要如何修改才能計算上面所說的介電常數(shù) ndtset 3 #set 1 : Initial self-consistent run iscf1 5 kptopt1 1 tolvrs1 1.0d-18 #need excellent convergence of GS quantities for RF runs # Set 2 : Calculate the ddk wf's - needed for piezoelectric tensor and # Born effective charges in dataset 3 getwfk2 -1 iscf2 -3 #this option is needed for ddk kptopt2 2 #use time-reversal symmetry only for k points nqpt2 1 #one wave vector will be specified qpt2 0 0 0 #need to specify gamma point rfelfd2 2 #set for ddk wf's rfdir2 1 1 1 #full set of directions needed tolwfr2 1.0d-20 #only wf convergence can be monitored here # Set 3 : response-function calculations for all needed perturbations getddk3 -1 getwfk3 -2 iscf3 5 kptopt3 2 #use time-reversal symmetry only for k points nqpt3 1 qpt3 0 0 0 rfphon3 1 #do atomic displacement perturbation rfatpol3 1 4 #do for all atoms rfstrs3 3 #do strain perturbation rfelfd3 3 efield3 0.0001 0.0001 0.0001 rfdir3 1 1 1 #the full set of directions is needed tolvrs3 1.0d-10 #need reasonable convergence of 1st-order quantities #Common input data # acell COPY RELAXED RESULT FROM PREVIOUS CALCULATION # Here is a default value, for automatic testing : suppress it and fill the previous line acell 6.0591285222E+00 6.0591285222E+00 9.7876758521E+00 rprim sqrt(0.75) 0.5 0.0 #hexagonal primitive vectors must be -sqrt(0.75) 0.5 0.0 #specified with high accuracy to be 0.0 0.0 1.0 #sure that the symmetry is recognized #and preserved in the optimization #process #Definition of the atom types and atoms ntypat 2 znucl 8 30 natom 4 typat 1 1 2 2 #Starting approximation for atomic positions in REDUCED coordinates xred 3.3333333333E-01 6.6666666667E-01 3.8063479051E-01 -3.3333333333E-01 -6.6666666667E-01 8.8063479051E-01 3.3333333333E-01 6.6666666667E-01 1.9652094936E-03 -3.3333333333E-01 -6.6666666667E-01 5.0196520949E-01 #Gives the number of bands, explicitely (do not take the default) nband 20 # For an insulator (if described correctly as an # insulator by DFT), conduction bands should not # be included in response-function calculations #Definition of the plane wave basis set ecut 50 # Maximum kinetic energy cutoff (Hartree) ecutsm 0.5 # Smoothing energy needed for lattice paramete # optimization. This will be retained for # consistency throughout. #Definition of the k-point grid kptopt 1 # Use symmetry and treat only inequivalent points ngkpt 8 8 8 # 4x4x4 Monkhorst-Pack grid nshiftk 1 # Use one copy of grid only (default) shiftk 0.0 0.0 0.5 # This choice of origin for the k point grid # preserves the hexagonal symmetry of the grid, # which would be broken by the default choice. #Definition of the self-consistency procedure diemac 3.0 # Model dielectric preconditioner iscf 5 # Use conjugate-gradient SCF cycle nstep 100 # Maxiumum number of SCF iterations # enforce calculation of forces at each SCF step optforces 1 |
木蟲 (小有名氣)
木蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 282求調(diào)劑 不挑專業(yè) 求收留 +4 | Yam. 2026-03-30 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 070300一志愿211,312分求調(diào)劑院校 +10 | 小黃鴨寶 2026-03-30 | 10/500 |
|
|
[考研] 070300化學354求調(diào)劑 +13 | 101次希望 2026-03-28 | 13/650 |
|
|
[考研] 學碩274求調(diào)劑 +16 | Li李魚 2026-03-26 | 16/800 |
|
|
[考研] 323分 食品與營養(yǎng)調(diào)劑 +3 | 嘿ooo 2026-03-31 | 3/150 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學 ,080500材料科學與工程學碩,總分289分 +9 | @taotao 2026-03-29 | 9/450 |
|
|
[考研] 一志愿華東師范大學有機化學專業(yè),初試351分,復試被刷求調(diào)劑! +7 | 真名有冰 2026-03-29 | 8/400 |
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
[考研] 抱歉 +3 | 田洪有 2026-03-30 | 3/150 |
|
|
[考研] 0703一志愿9,初試成績:338,四六級已過,有科研經(jīng)歷,求調(diào)劑! +7 | Zuhui0306 2026-03-25 | 7/350 |
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +18 | Mr. Z 2026-03-25 | 19/950 |
|
|
[考研] 化學0703 調(diào)劑 306分 一志愿211 +8 | 26要上岸 2026-03-28 | 8/400 |
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-24 | 6/300 |
|
|
[考研] 291求調(diào)劑 +5 | Y-cap 2026-03-29 | 6/300 |
|
|
[考研] 352分-085602-一志愿985 +5 | 海納百川Ly 2026-03-29 | 5/250 |
|
|
[考研] 330分求調(diào)劑 +5 | qzenlc 2026-03-29 | 5/250 |
|
|
[考研]
|
y7czhao 2026-03-26 | 10/500 |
|
|
[考研] 考研調(diào)劑 +10 | 呼呼?~+123456 2026-03-24 | 10/500 |
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|