| 查看: 3925 | 回复: 11 | ||||
yashumimi银虫 (小有名气)
|
[交流]
【求助】二维材料k点的设置 已有1人参与
|
|||
|
今天在算k点不同对transiesta计算结果的影响 k点设置为 1 0 0 0 1 0 0 0 i 其中i为变量 我将其设为2－20 因为开始用1算了 出错 是个二维的结构 x方向有一定的宽度 没有周期性 加了真空层 y方向不变 z方向具有周期性 请问k点应该怎么设置才合理呢？ 另外我更改k点后用tbtrans计算了I－V 同样的操作 两次 电流数量级有小幅度变化 趋势变化规律也不一样 这说明计算的结果不稳定么？ |
回复此楼» 收录本帖的淘帖专辑推荐
 transport |
» 猜你喜欢
93BebMhtakh前后11位开头都是大写
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有86人回复
-
感觉是下周放榜了
已经有11人回复
-
filecode=后面第一个是大写字母
已经有10人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1974.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
★ ★ ★ ★ ★ ★
cenwanglai(金币+3):谢谢广平兄答复！ 2010-11-09 20:29:56
yashumimi(金币+5):谢谢 2010-11-22 18:05:16
zzy870720z(金币+3):谢谢zhang兄指教 2010-11-22 18:39:08
cenwanglai(金币+3):谢谢广平兄答复！ 2010-11-09 20:29:56
yashumimi(金币+5):谢谢 2010-11-22 18:05:16
zzy870720z(金币+3):谢谢zhang兄指教 2010-11-22 18:39:08
|
对于scattering region一般是通过优化得到中间结构的。但是一旦做输运计算的时候，就不在考虑偏压下的几何结构变化了。因为这个问题Transiesta是处理不了的。再说即使处理得了，这个时间也是负担不起的。所以这个地方就有一个近似了。至于transiesta计算中输出的力，已经不再具有参考价值。不过听说atk-2010.8可以在偏压下进行几何结构的优化。具体的执行方法我不是很清楚，但是明确的告诉你，transiesta不能胜任这个工作。 所以在做输运的时候，要将MD.NumCGsteps 设置为0，也就是说所有的关于MD的选项不再使用了。 关于这个你可以到siesta自带的\Tests\transiesta目录下，查看一下所有的散射区输运计算的输入文件都是这样设置的。 关于这个我就知道这么写，如果有疑问，可以到siesta的mail_list上提问。 |
11楼2010-11-09 20:15:29
haowenping
木虫 (著名写手)
- 应助: 6 (幼儿园)
- 金币: 2330.2
- 散金: 384
- 红花: 5
- 帖子: 1104
- 在线: 244.5小时
- 虫号: 1124850
- 注册: 2010-10-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2010-11-05 22:37:18
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1974.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
★ ★ ★
yashumimi(金币+1):谢谢参与
mazuju028(金币+2):很详细。 2010-11-09 17:28:22
yashumimi(金币+1):谢谢参与
mazuju028(金币+2):很详细。 2010-11-09 17:28:22
|
你把你的输入文件贴出来，高手可以用你的输入文件具体查看一下，有时候你只是这么一说很难排查。我也想知道这个问题的答案。 我一直感觉对于碳纳米管这种一维体系，k点设置如下就行： %block kgrid_Monkhorst_Pack 1 0 0 0.0 0 1 0 0.0 0 0 6 0.0 %endblock kgrid_Monkhorst_Pack 其中第一二行都设置为1，因为这两个方向是没有周期性的。（碳管沿着z方向）。第三行是为了计算电极部分的碳管用的，对于碳管这种体系，可能这个地方需要的k点数目比较多。需要做一个收敛测试。但是对于输运计算的时候，第三行就没设么用了。只看第一二行。 |
3楼2010-11-08 09:40:28
yashumimi
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 387.5
- 散金: 40
- 帖子: 68
- 在线: 71.2小时
- 虫号: 461676
- 注册: 2007-11-18
- 专业: 原子和分子物理
4楼2010-11-08 13:56:48
yashumimi
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 387.5
- 散金: 40
- 帖子: 68
- 在线: 71.2小时
- 虫号: 461676
- 注册: 2007-11-18
- 专业: 原子和分子物理
5楼2010-11-08 13:57:13
yashumimi
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 387.5
- 散金: 40
- 帖子: 68
- 在线: 71.2小时
- 虫号: 461676
- 注册: 2007-11-18
- 专业: 原子和分子物理
|
SystemName SystemLabel %block ChemicalSpeciesLabel 1 1 H 2 6 C %endblock ChemicalSpeciesLabel PAO.EnergyShift 50 meV PAO.BasisType split PAO.BasisSize DZP %block kgrid_Monkhorst_Pack 1 0 0 0.0 0 1 0 0.0 0 0 20 0.0 %endblock kgrid_Monkhorst_Pack xc.functional GGA # Exchange-correlation functional xc.authors PBE # Exchange-correlation version SpinPolarized F # Logical parameters are: yes or no FixSpin F MeshCutoff 300. Ry # Equivalent plane wave cutoff for the grid # SCF options DM.Tolerance 1.d-4 # Tolerance in maximum difference between input and output DM MaxSCFIterations 1000 # Maximum number of SCF iter DM.UseSaveDM T # to use continuation files DM.MixingWeight 0.02 # New DM amount for next SCF cycle DM.NumberPulay 6 DM.MixSCF1 F DM.PulayOnFile F # Store in memory ('F') or in files ('T') # NeglNonOverlapInt T # 'F'=do not neglect SolutionMethod Diagon # OrderN or Diagon ElectronicTemperature 300 K # Temp. for Fermi smearing # MD options MD.TypeOfRun CG # Type of dynamics: MD.VariableCell F MD.NumCGsteps 100 # Number of CG steps for coordinate optimization %include POSITION.fdf 输入文件差不多就是这个样子的 有个问题用一个小的架构优化后supercell 力还是很大 只能再优化么 有没有好的办法提高优化速度 另外用siesta优化和transiesta优化的区别在哪里? |
6楼2010-11-08 14:25:10
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1974.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
7楼2010-11-08 15:15:05
yashumimi
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 387.5
- 散金: 40
- 帖子: 68
- 在线: 71.2小时
- 虫号: 461676
- 注册: 2007-11-18
- 专业: 原子和分子物理
8楼2010-11-08 19:01:26
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1974.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
9楼2010-11-08 19:49:18
yashumimi
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 387.5
- 散金: 40
- 帖子: 68
- 在线: 71.2小时
- 虫号: 461676
- 注册: 2007-11-18
- 专业: 原子和分子物理
|
Dear Transiesta user: I have one question about the calculation of Transiesta. I am using Transiesta to calculate the transmission in nanowires. I am not sure about the optimization of the scattering region using siesta or transiesta. Since if I use siesta to optimize the scattering region, the results are not correct and if I use the Transiesta to optimize the scattering region, it can not converge . Could anybody help me? Thanks in advance. Lihong 我觉得我和他遇到的问题是一样的 就是加上偏压后力变得很大 |
10楼2010-11-09 17:09:46
