| 5 | 1/1 | 返回列表 |
| 查看: 1607 | 回復(fù): 4 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
[求助]
如何計(jì)算clamped-ion 和free_stress 狀態(tài)下的介電常數(shù)
|
||
|
這是本人的.in文件,請(qǐng)問(wèn)告訴要如何修改才能計(jì)算上面所說(shuō)的介電常數(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 |
木蟲(chóng) (小有名氣)
木蟲(chóng) (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +7 | 研研,接電話 2026-03-24 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 287求調(diào)劑 +11 | land xuxu 2026-03-26 | 11/550 |
|
|
[考研] 0703化學(xué) +14 | 妮妮ninicgb 2026-03-27 | 14/700 |
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +14 | jhybd 2026-03-23 | 19/950 |
|
|
[考研] 0703 化學(xué) 求調(diào)劑,一志愿山東大學(xué) 342 分 +6 | Shern—- 2026-03-28 | 6/300 |
|
|
[考研] 0856求調(diào)劑 +8 | 楒桉 2026-03-28 | 8/400 |
|
|
[考研] 327求調(diào)劑 +6 | 汲亦昊 2026-03-29 | 6/300 |
|
|
[考研] 一志愿鄭州大學(xué),080500學(xué)碩,總分317分求調(diào)劑 +8 | 舉個(gè)栗子oi 2026-03-24 | 9/450 |
|
|
[考研] 332求92調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
[考研] 295材料工程專碩求調(diào)劑 +7 | 1428151015 2026-03-27 | 7/350 |
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
[考研] 266求調(diào)劑 +11 | 陽(yáng)陽(yáng)哇塞 2026-03-27 | 12/600 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
[考研] 求調(diào)劑323材料與化工 +7 | 1124361 2026-03-24 | 7/350 |
|
|
[考研] 機(jī)械學(xué)碩310分,數(shù)一英一,一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
[考研] 打過(guò)很多競(jìng)賽,085406控制工程300分,求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|