| 查看: 2412 | 回复: 10 | ||||
[交流]
【求助】ATK中的电极温度能表征系统的真实温度吗？ 已有5人参与
|
||||
|
ATK中的电极温度决定了费米分布函数 它能表征真实的温度吗？ 如果不能得话，这个参数的意义在哪？ |
回复此楼» 收录本帖的淘帖专辑推荐
 SIESTA |
» 猜你喜欢
315求调剂
已经有15人回复
-
070300化学学硕311分求调剂
已经有9人回复
-
一志愿南昌大学，085600，344分求调剂
已经有7人回复
-
一志愿 江南大学 085602 化工专硕 338分求调剂
已经有13人回复
-
化学0703-一志愿211-338分求调剂
已经有8人回复
-
0703调剂，一志愿天津大学319分
已经有8人回复
-
085602调剂 初试总分335
已经有9人回复
-
336材料与化工085600求调剂
已经有8人回复
-
085600材料与化工301分求调剂院校
已经有6人回复
-
071000生物学调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 变温粉末表征中峰在不同温度下出现峰的强弱是什么原因
已经有3人回复
- 电极表面修饰物如何进行SEM、IR表征？
已经有19人回复
- 参比电极室的温度要和工作电极室温度一致么
已经有3人回复
- 【交流】谁有Au电极在硫酸中表征的标准CV图吗？
已经有8人回复
- 【求助】DSC和DMA表征玻璃化转变温度
已经有29人回复
- 【求助】修饰金电极红外表征
已经有5人回复
- 【求助】求教：表征裸金电极
已经有6人回复
★ ★ ★
zzy870720z(金币+3):感谢专家指导 2010-05-20 12:26:16
zzy870720z(金币+3):感谢专家指导 2010-05-20 12:26:16
|
T当然不是真实体系温度。 T只是 fermi-dirac smear参数，你可以把它看做展宽时的能量(kB*T)。 这么做，是为了费米面上下电子的占据分数平滑过渡，特别对于金属，费米面处，是光滑的。 这样在计算时候，也能加速收敛。 再次说明， 第一原理是进行基态(绝对零度下)的计算。不可能算高温激发态。你可能说看到好多paper算啊。你接着往下看。 有限温度下的计算，常用的近似方法如准粒子GW近似等。 在很多软件中都有。 但是，注意，这也仅限于一种近似方法。 在有限温度下，各种元激发，如声子、电子、磁子相互作用。 是无法确切计算出来的。 有人号称是计算的高温如 100K 300K 或更高温度的，那也只是一种在有限温度下进行用分子动力学的方法进行晶体结构优化，然后还得回到零温下进行第一原理的自洽计算。 [ Last edited by wenjlu on 2010-5-20 at 11:39 ] |
5楼2010-05-20 11:35:05
2楼2009-08-05 11:22:26
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
★ ★
zzy870720z(金币+2):谢谢交流 2010-05-20 12:24:47
zzy870720z(金币+2):谢谢交流 2010-05-20 12:24:47
|
我觉得费米分布函数里面的温度指的就是电子所处的周围环境的温度！分布函数中的T的含义指的是温度引起的电子分布在费米能级处的展宽。但是设置了这个温度，计算的结果却不能说是体系在该温度下的性质。这个要看研究的是什么性质了。要是研究振动性质。那么可以说结果就是相应温度的结果，但是要是研究电子的性质，300K的和0K的，对体系分子能级的影响是可以忽略的，但是设置了一个较大的温度之后，对计算中收敛是很有利的。这就是我们在计算的时候为什么设置一个不是0K的温度，来预测低温的性质。这个是有前提的，要看你研究的是什么性质，温度对这个性质的影响有多大。在计算电子结构的时候，我们可以把温度看成是收敛的一个技巧。但是和直接说这个温度不是体系的温度还是有差异的。因为分布函数中温度的意义已经摆在那个地方了。那样说不是一种矛盾吗？ 另外，在siesta中还有一种电子占据的分布方式，不是费米分布，而是别的一种分布方式，那里面的温度绝对是一种收敛手法。 这是我个人的理解，不对的地方敬请指教。以后多交流讨论。谢谢。 |
3楼2010-05-20 09:57:39
夕阳西下
金虫 (著名写手)
- 应助: 138 (高中生)
- 贵宾: 0.041
- 金币: 9247.1
- 散金: 73
- 红花: 28
- 沙发: 1
- 帖子: 1268
- 在线: 1164.4小时
- 虫号: 220750
- 注册: 2006-03-20
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 物理
4楼2010-05-20 11:23:44
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
6楼2010-05-23 16:26:39
acridine
木虫 (正式写手)
- 1ST强帖: 11
- 应助: 13 (小学生)
- 金币: 4901.6
- 红花: 13
- 帖子: 706
- 在线: 389.5小时
- 虫号: 305478
- 注册: 2006-12-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★
gavinliu7390(金币+2):谢谢交流！ 2010-06-01 10:53:27
gavinliu7390(金币+2):谢谢交流！ 2010-06-01 10:53:27
|
这个问题仔细琢磨琢磨还是很有意思的。 其实本质上看，这个参数所的设置东西，无非是在DFT的scf自洽计算的同时，人为的约束电子的排布情况，让电子按照某个温度下的分布函数来占据能级。但这种约束DFT计算的结果是否就是真实中某个温度下的情况？应该没有那么简单，如果这么简单的话，那么激发态的计算也不用研究了，在DFT计算时人为的约束Ef之下的某个能级是空的，之上的某个能级是占据的就可以了。 按照这个方式计算激发态的一个重要问题在于，没有考虑激发上去的电子与留下的空位之间的作用；这里的问题也应该类似吧，按照有限温度的Fermi分布重新排布电子之后，在Fermi能级之上有分数占据的电子，在Fermi能级之下有分数非占据的空位，这两者之间的作用也是使用这种约束DFT计算所没有考虑的吧。所以与实际情况还是有一定的差别的。这种差别具体是怎样的，对分析问题有多大的影响，俺就不懂了。 |
7楼2010-05-23 21:16:01
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
8楼2010-05-23 22:00:01
acridine
木虫 (正式写手)
- 1ST强帖: 11
- 应助: 13 (小学生)
- 金币: 4901.6
- 红花: 13
- 帖子: 706
- 在线: 389.5小时
- 虫号: 305478
- 注册: 2006-12-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★
gavinliu7390(金币+2):谢谢交流！ 2010-06-01 11:06:17
gavinliu7390(金币+2):谢谢交流！ 2010-06-01 11:06:17
|
我上面的回复是说这个虚拟出来的温度与真实有限温度的情况是有差别的 对于分子、半导体体系，这个参数当然没有什么影响，所以在计算这类材料 smearing设置成0也没有问题。但在金属时候，Ef附近电子态很密集，这个参数将影响电子排布。如你所说，使用这个参数是可以使分布函数平滑，进而加速收敛的。在这一点上大家应该有共识的。 但这并不与我的回复矛盾吧。对于分子体系，如你所说，这个参数对计算电子结构基本没有影响，那么它能否反映真实温度也就没有意义了吧。对于金属，如同我上面的理解，由于方法中并没有考虑某一些相互作用，不能反映有限温度下的真实情况。但作为一个加速收敛的工具，它还是很有价值的。 另外，在DFT的scf过程中当然使用这个参数了，否则何来这个参数能加速收敛一说呢？  |
9楼2010-05-23 22:52:57
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
★ ★
gavinliu7390(金币+2):谢谢交流！ 2010-06-01 11:08:03
gavinliu7390(金币+2):谢谢交流！ 2010-06-01 11:08:03
|
我只知道这个费米分布函数在算输运的电流时在朗道公式中出现。但是单纯的siesta计算中也需要设置这个参数，所以在DFT计算中也使用了这个参数。使用地方是这个位置：  （点击一下图片，会在新窗口打开，那样图片就会变清楚！） 这里面说的“DFT计算的结果原则上给出体系的基态”，这里的基态应该是电子的基态，因为在计算的过程中没有考虑振动效应。所以也就是体系的基态了！ 对于半导体和分子体系，这样做给出的结果确实和不做费米分布的时候相差不大，但是，对于金属体系，能带结构，能级连续，所以计算结果还是有一定的差别的。差别在于费米能级附近的电子布局。但是这种差别在0.03eV之内（300K）。 [ Last edited by zhangguangping on 2010-6-1 at 03:48 ] |
10楼2010-06-01 09:09:21
