| 5 | 1/1 | 返回列表 |
| 查看: 1605 | 回復: 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 |
木蟲 (小有名氣)
木蟲 (小有名氣)
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 291求調劑 +7 | Y-cap 2026-03-29 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 277跪求調劑 +6 | 1915668 2026-03-27 | 10/500 |
|
|
[考研] 349求調劑 +6 | 李木子啊哈哈 2026-03-25 | 6/300 |
|
|
[考研] 289求調劑 +13 | 新時代材料 2026-03-27 | 13/650 |
|
|
[考研] 315求調劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
[考研] 調劑 +3 | 好好讀書。 2026-03-28 | 3/150 |
|
|
[有機交流]
高溫高壓反應求助
10+4
|
chibby 2026-03-25 | 4/200 |
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
[考研] 一志愿南師大0703化學 275求調劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
[考研] 一志愿華東理工大學081700,初試分數(shù)271 +6 | kotoko_ik 2026-03-23 | 7/350 |
|
|
[考研] 求調劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
[考研] 0703化學一志愿南京師范大學303求調劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
[考研] 0703化學338求調劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研]
|
鐘llll 2026-03-26 | 4/200 |
|
|
[考研] 材料專碩 335 分求調劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
[考研] 0854AI CV方向招收調劑 +4 | 章小魚567 2026-03-23 | 4/200 |
|
|
[考研] 285求調劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
[考研] 340求調劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
[考研] 344求調劑 +3 | desto 2026-03-24 | 3/150 |
|
|
[考研] 328求調劑 +4 | LHHL66 2026-03-23 | 4/200 |
|