应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 【求助】态密度出现一些很小的负值是怎么回事？
北京石油化工学院2026年研究生招生接收调剂公告
182/2返回列表上一页12
查看: 2047  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

qscao

木虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
请问楼主，你能否传一个ATK得Lic文件给我啊，谢谢！！
邮箱：qscao@qq.com
赞一下(1人) 回复此楼
11楼2009-10-18 18:27:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamddba

铁杆木虫 (正式写手)
★ ★ ★ ★
zdhlover(金币+2,VIP+0):谢谢热情的你 10-18 13:17
wuchenwf(金币+2,VIP+0):谢谢回帖交流 10-22 21:09
回楼上的 qscao ，我用的也是原来宏剑公司发布的试用license，早就过期了，没办法只好改系统时间为2008年，呵呵，无奈之举啊，主要是没钱买啊，太贵了。不过我还是把以前我用的许可文件给你发过去了，请查收
赞一下(20人) 回复此楼
12楼2009-10-19 08:58:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

freshgirl

木虫 (正式写手)
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuchenwf(金币+5,VIP+0):呵呵，列的真详细啊，呵呵 10-22 21:09
引用回帖:
Originally posted by xhsh at 2009-8-17 09:05:


与自旋极化没有关系，不要误导人。

还是用文献来说话吧：
1. Phys. Rev. B 78, 155130 (2008)
"Negative density of states: Screening, Einstein relation, and negative diffusion"

2. Physica B: Condensed Matter, Volume 169, Issues 1-4, 2 February 1991, Pages 559-560
"Negative density of states of a two-dimensional electron gas"

不过个人同意之前说的，ATK中出现负DOS，是积分能量点不够的原因。
赞一下(11人) 回复此楼
13楼2009-10-22 10:34:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhsh

木虫 (小有名气)
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+5,VIP+0):感谢您的积极讨论 10-23 08:38
引用回帖:
Originally posted by freshgirl at 2009-10-22 10:34:

还是用文献来说话吧：
1. Phys. Rev. B 78, 155130 (2008)
"Negative density of states: Screening, Einstein relation, and negative diffusion"

2. Physica B: Condensed Matter, Volume 169 ...

应该强调，我们这里说的态密度指的是能态密度，也就是体系中单位能量间能级的 数目。在某个能量范围内，只可能有”有能级“或”没有能级“两种情况。有就是正 的，没有就是0。 说这个能级数目是负的，本身是没有意义的。能态密度的定义是：

N(E)=lim(Δ Z/ΔE）

其中ΔE是能量间隔，ΔZ是这个能量间隔中的能态数目。

楼上所列的两篇文献中讨论的是所谓的thermodynamic density of states，虽然译出来也有“态密度”的含义，但是并不是我们所指的“能态密度”。其定义是：

A(n)=lim(Δμ/Δn)

其中μ为体系中的化学势，n为电子密度分布。

本人对thermodynamic density of states知之不多，但可以肯定的是，两个态密度不是一回事。
赞一下(7人) 回复此楼
14楼2009-10-23 08:36:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

freshgirl

木虫 (正式写手)
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+5,VIP+0):谢谢解答，但个人认为在态密度图中上下自旋分为正负只是便于区分和讨论 10-23 14:38
引用回帖:
Originally posted by xhsh at 2009-10-23 08:36:


应该强调，我们这里说的态密度指的是能态密度，也就是体系中单位能量间能级的 数目。在某个能量范围内，只可能有”有能级“或”没有能级“两种情况。有就是正 的，没有就是0。 说这个能级数目是负的，本身是没 ...

谢谢纠正。

我同意你用有和没有来解释DOS的正值以及零值。

但是，之前我看过一些文章，APL水平的，会把spin up和spin down分别画成正负。我们之所以一直认为态密度是正的，是由于纵坐标写成"abs. counts"，但是你想啊，unoccupied DOS，也就是虚频部分算出来的结果，波函数平方后的结果就是负值。也就是说，虽然从物理角度说，DOS不应该有负。但是数学角度算出来，确实有负值。所以这个态密度的正负问题，我个人比较困惑。

有兴趣好好讨论一下。
赞一下(7人) 回复此楼
15楼2009-10-23 14:28:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aylayl08

荣誉版主 (文坛精英)

至尊天蝎

优秀版主
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
qasd(金币+2,VIP+0):谢谢~ 10-25 12:50
引用回帖:
Originally posted by freshgirl at 2009-10-23 14:28:

谢谢纠正。

我同意你用有和没有来解释DOS的正值以及零值。

但是，之前我看过一些文章，APL水平的，会把spin up和spin down分别画成正负。我们之所以一直认为态密度是正的，是由于纵坐标写成"abs. co ...

从态密度单位来看，electrons/eV就可以看出电子数不会为负值。个人想法
赞一下(11人) 回复此楼
珍惜拥有，快乐人生!欢迎光临计算模拟区(*^__^*)。。。。。。。。。。。。。。
16楼2009-10-23 14:40:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhsh

木虫 (小有名气)
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
aylayl08(金币+5,VIP+0):感谢解释，欢迎常来 10-24 09:41
引用回帖:
Originally posted by freshgirl at 2009-10-23 14:28:

谢谢纠正。

我同意你用有和没有来解释DOS的正值以及零值。

但是，之前我看过一些文章，APL水平的，会把spin up和spin down分别画成正负。我们之所以一直认为态密度是正的，是由于纵坐标写成"abs. co ...

不只是APL水平的，PRL水平的都有，为了表示清楚和对比，通常都把自旋向上和自 旋向下分别画成正负，但请注意，这只是画成正负，并不代表实际的值就是负的。 这只是在作图的时候的一种处理手法，千万不要因为你看到文献的图里画的是负的就认为是负的。

“unoccupied DOS”，指的是没有占据的态密度，也就是费米能级以上的态密度。怎 么会是虚频计算出来的结果呢？这个概念是不是还有别的含义？愿意听你详细解 释。但有一点，算电子密度的时候（你说的是求态密度，这是错的，纠正一下）， 可不是直接对波函数求平方，而是波函数乘以波函数的复共轭，这个永远都是实的 且是非负的。态密度只与能级有关，具体计算的时候采用类似如下的计算公式：

N(E)= Σδ/((E-Ei)^2+δ^2)

式中Σ为对所有能级求和，Ei为能级，δ为展宽系数。
赞一下(7人) 回复此楼
17楼2009-10-24 09:01:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

freshgirl

木虫 (正式写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
qasd(金币+2,VIP+0):谢谢~ 10-25 12:50
引用回帖:
Originally posted by xhsh at 2009-10-24 09:01:


不只是APL水平的，PRL水平的都有，为了表示清楚和对比，通常都把自旋向上和自 旋向下分别画成正负，但请注意，这只是画成正负，并不代表实际的值就是负的。 这只是在作图的时候的一种处理手法，千万不要因为你 ...

ATK计算DOS时会出现负值。我自己算出过这样的结果。虽然有办法避免，但是由于不知道ATK具体是如何处理的，所以很困惑。

你给的DOS计算公式，只是一个大概的算法，具有普遍通用性。但是不同的模型，都有不同的DOS计算公式。怎样用来解释现在碰到这个问题呢？

看Gaussian输出文件的时候，判别HOMO和LUMO，是看虚频和实频分界。我不是指一个复函数的虚频部分。是不是这里的理解有点差异？

你说电子密度是波函数的模的平方。我赞同。我的意思是，有没有这种可能，在计算中引入了乘法计算，导致了负数结果。或者是积分的时候。

请多指教。
赞一下(11人) 回复此楼
18楼2009-10-25 12:35:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 iamddba 的主题更新
182/2返回列表上一页12
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物调剂 +3 橙子橙子橙子啊 2026-04-05 6/300 2026-04-07 14:06 by 18828373951
[考研] 生物医药调剂｜SCI中科院三区一作+多项科研成果 +3 likangxing 2026-04-07 5/250 2026-04-07 13:13 by guanxin1001
[考研] 265求调剂 +17 小木虫085600 2026-04-06 19/950 2026-04-07 12:04 by Sammy2
[考研] 385分 生物学（071000）求调剂 +12 qf626 2026-04-01 12/600 2026-04-06 17:36 by 52305043001
[考研] 化学调剂求助 +8 LULONG1 2026-04-03 8/400 2026-04-06 10:26 by dongzh2009
[考研] 0857大类环境工程B区求调剂 +3 龚禹铭 2026-04-05 3/150 2026-04-06 10:22 by 蓝云思雨
[考研] 283分求调剂 +7 小聂爱学习 2026-04-03 7/350 2026-04-04 21:51 by hemengdong
[考研] 0835学硕299求调剂 08大类可接受 +5 useryy 2026-04-03 5/250 2026-04-04 20:07 by 蓝云思雨
[考研] 。 +5 雾与海 2026-04-02 6/300 2026-04-04 19:53 by 蓝云思雨
[考研] 一志愿华南师范361分，化学求调剂 +7 Nicole88888 2026-04-01 7/350 2026-04-04 18:28 by macy2011
[考研] 321求调剂 +13 认真求上学 2026-04-02 13/650 2026-04-04 18:23 by macy2011
[考研] 一志愿0817化学工程与技术，求调剂 +24 我不是只因 2026-04-02 28/1400 2026-04-04 15:15 by dongzh2009
[考研] 372分材料与化工（085600）一志愿湖南大学求调剂 +3 蓝笺片 2026-04-03 4/200 2026-04-03 17:58 by Jimmyandyou
[硕博家园] 求老师收留 +9 lllq123 2026-04-03 9/450 2026-04-03 13:48 by 呼吸都是减肥
[考研] 071000生物学调剂 +8 知昭蔓 2026-04-02 8/400 2026-04-03 10:36 by macy2011
[考研] 312求调剂 +4 赊月色 2026-04-02 5/250 2026-04-03 08:21 by fangshan711
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-04-02 3/150 2026-04-02 15:06 by cal0306
[考研] 298求B区调剂 +4 zzz，，r 2026-04-02 5/250 2026-04-02 12:17 by 土木硕士招生
[考研] 266求调剂 +4 学员97LZgn 2026-04-02 4/200 2026-04-02 09:52 by yulian1987
[考研] 292求调剂 +17 木虫er12138 2026-04-01 17/850 2026-04-01 21:37 by 七度不信任
信息提示
关闭
请填处理意见
关闭