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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 【求助】声子在siesta中的计算问题
北京石油化工学院2026年研究生招生接收调剂公告
81/1返回列表
查看: 762  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

dalongmao

木虫 (正式写手)

[交流] 【求助】声子在siesta中的计算问题

我使用siesta计算声子，MD.TypeOfRun=FC，
这是我采用的特征点
1  0.000  0.000  0.000  G                                                
15  1.000  0.000  0.000       X                                                
25  1.000  0.000  1.000       L                                                
20  0.000  0.000  0.000       G      
得到vectors文件在附件中，其中Eigenvector，Eigenmode (real part)和Eigenmode (imaginary part)分别代表什么意思？siesta中虚频指的是什么？怎么区分？我的结构稳定吗？如果出现虚频，怎样消除使其稳定？
然后我使用vib2xsf使其可视化，得到了.XSF file 和.AXSF文件，但是其中有NaN，怎么消除NaN？
下面是.XSF文件
# ---- XSF block for ---- iev =     6  freq =    -171.135485 cm-1
ATOMS
   0   0.0000000   0.0000000   0.0000000   NaN         NaN         NaN        
   0   8.2053923   0.0000002   9.9912885   NaN         NaN         NaN        
   0   8.2053922   8.2053927   9.9912887   NaN         NaN         NaN        
   1   8.2053925   2.0169873   2.7592412     0.0000000   0.0000000   0.0000000
   1   4.1026965   6.1196834   7.7548862     0.0000000   0.0000000   0.0000000
   5   5.9871737   4.1026963   3.4703564     0.0000000   0.0000000   0.0000000
下面是.AXSF文件
# --- AXSF block for ---- iev =     6  freq =    -171.135485 cm-1
ANIMSTEPS   10
ATOMS    1
   0 NaN         NaN         NaN        
   0 NaN         NaN         NaN        
   0 NaN         NaN         NaN        
   1   8.2053925   2.0169873   2.7592412
   1   4.1026965   6.1196834   7.7548862
   5   5.9871737   4.1026963   3.4703564
ATOMS    2
   0 NaN         NaN         NaN        
   0 NaN         NaN         NaN        
   0 NaN         NaN         NaN        
   1   8.2053925   2.0169873   2.7592412
   1   4.1026965   6.1196834   7.7548862
   5   5.9871737   4.1026963   3.4703564
ATOMS    3
   0 NaN         NaN         NaN        
   0 NaN         NaN         NaN        
   0 NaN         NaN         NaN        
   1   8.2053925   2.0169873   2.7592412
   1   4.1026965   6.1196834   7.7548862
   5   5.9871737   4.1026963   3.4703564
ATOMS    4
   0 NaN         NaN         NaN        
   0 NaN         NaN         NaN        
   0 NaN         NaN         NaN        
   1   8.2053925   2.0169873   2.7592412
   1   4.1026965   6.1196834   7.7548862
   5   5.9871737   4.1026963   3.4703564
ATOMS    5
   0 NaN         NaN         NaN        
   0 NaN         NaN         NaN        
   0 NaN         NaN         NaN        
   1   8.2053925   2.0169873   2.7592412
   1   4.1026965   6.1196834   7.7548862
   5   5.9871737   4.1026963   3.4703564
ATOMS    6
   0 NaN         NaN         NaN        
   0 NaN         NaN         NaN        
   0 NaN         NaN         NaN        
   1   8.2053925   2.0169873   2.7592412
   1   4.1026965   6.1196834   7.7548862
   5   5.9871737   4.1026963   3.4703564
ATOMS    7
   0 NaN         NaN         NaN        
   0 NaN         NaN         NaN        
   0 NaN         NaN         NaN        
   1   8.2053925   2.0169873   2.7592412
   1   4.1026965   6.1196834   7.7548862
   5   5.9871737   4.1026963   3.4703564
ATOMS    8
   0 NaN         NaN         NaN        
   0 NaN         NaN         NaN        
   0 NaN         NaN         NaN        
   1   8.2053925   2.0169873   2.7592412
   1   4.1026965   6.1196834   7.7548862
   5   5.9871737   4.1026963   3.4703564
ATOMS    9
   0 NaN         NaN         NaN        
   0 NaN         NaN         NaN        
   0 NaN         NaN         NaN        
   1   8.2053925   2.0169873   2.7592412
   1   4.1026965   6.1196834   7.7548862
   5   5.9871737   4.1026963   3.4703564
ATOMS   10
   0 NaN         NaN         NaN        
   0 NaN         NaN         NaN        
   0 NaN         NaN         NaN        
   1   8.2053925   2.0169873   2.7592412
   1   4.1026965   6.1196834   7.7548862
   5   5.9871737   4.1026963   3.4703564
还有我想做出声子态密度，怎样得到？请高手不惜赐教！

[ Last edited by dalongmao on 2009-6-16 at 15:17 ]
回复此楼

» 猜你喜欢
1楼 2009-06-15 10:37:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dummy1

木虫 (正式写手)
★ ★ ★
yjr(金币+1,VIP+0):多谢！！ 6-26 23:11
dalongmao(金币+2,VIP+0):3ks 6-27 09:16
Eigenvector: 在这里表示第几个模。随后是这个振动模式的频率值 。
Eigenmode (real part，imaginary part)：应该是表示本征向量的值，实虚我也不确定。
虚频一般说明结构不稳定（至少是在VASP和gaussian）。需要重新优化一下
NaN：一般在fortran中表示出错或者是数值无穷，不能正常输出。

另外：你的  vector文件是怎么得到？能否发个输入文件的例子？我刚刚学siesta，看到网上说需要一个vibra辅助，是不是呢？是的话能不能发给我一个，我也想做一下频率。谢谢！
赞一下(10人) 回复此楼
2楼2009-06-26 21:30:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hchfox

木虫 (正式写手)
★
dalongmao(金币+1,VIP+0):3ks 6-27 09:16
建议你别用SIESTA算声子，换其他的吧。。。
引用回帖:
Originally posted by dalongmao at 2009-6-15 10:37:
我使用siesta计算声子，MD.TypeOfRun=FC，
这是我采用的特征点
1  0.000  0.000  0.000  G                                                
15  1.000  0.000  0.000       X                               ...

赞一下(4人) 回复此楼
3楼2009-06-26 23:31:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dummy1

木虫 (正式写手)
如果体系比较大的话你觉得什么程序可行呢？
谢谢！
引用回帖:
Originally posted by hchfox at 2009-6-26 23:31:
建议你别用SIESTA算声子，换其他的吧。。。
赞一下 回复此楼
4楼2009-06-27 09:23:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hchfox

木虫 (正式写手)
★ ★ ★ ★
zeoliters(金币+3,VIP+0):感谢回贴交流 6-27 13:28
dalongmao(金币+1,VIP+0):3ks 6-29 10:50
体系有多大？
pwscf，abinit实现了DFPT；VASP可采用冻结声子（直接方法）考虑（好像vasp5.2也实现DFPT了吧，不清楚）。。。
我用直接方法，超胞原子数都达到了240个原子，得到的结果还是很好（当然我的体系电子数是比较少的）。
不管怎样，声子计算是很耗时的，有时特别是体系大了，虚频可能很难消除。。。
引用回帖:
Originally posted by dummy1 at 2009-6-27 09:23:
如果体系比较大的话你觉得什么程序可行呢？
谢谢！
赞一下(13人) 回复此楼
5楼2009-06-27 13:18:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dummy1

木虫 (正式写手)
你是用多少个CPU算的，我实在是有点忍受不了了，慢的要是。
引用回帖:
Originally posted by hchfox at 2009-6-27 13:18:
体系有多大？
pwscf，abinit实现了DFPT；VASP可采用冻结声子（直接方法）考虑（好像vasp5.2也实现DFPT了吧，不清楚）。。。
我用直接方法，超胞原子数都达到了240个原子，得到的结果还是很好（当然我的体系电子 ...

赞一下 回复此楼
6楼2009-06-27 17:01:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hchfox

木虫 (正式写手)
16个CPU。。。
引用回帖:
Originally posted by dummy1 at 2009-6-27 17:01:
你是用多少个CPU算的，我实在是有点忍受不了了，慢的要是。
回复此楼
7楼2009-06-27 17:28:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

y1ding

铁杆木虫 (著名写手)
★ ★
dalongmao(金币+2,VIP+0):3ks,在哪里可以看见hchfox的意见? 6-30 09:58
嘿嘿,也看到有人用siesta算声子了.
我用siesta算过虚频怎么也消除不了,改用pw计算就不存在.
用siesta算声子需要很好的基组,用一般的DZP不太好,只能算算孤立体系G点.
还是看hchfox的意见比较有实际意义.
赞一下(3人) 回复此楼
8楼2009-06-30 09:55:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dalongmao 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 288求调剂，一志愿华南理工大学071005 +3 ioodiiij 2026-04-08 3/150 2026-04-08 06:21 by belle15
[考研] 308求调剂 +12 墨墨漠 2026-04-06 12/600 2026-04-08 00:29 by JourneyLucky
[考研] 11408 325分 +3 jgtxuxgkx 2026-04-07 3/150 2026-04-07 23:10 by lbsjt
[考研] 材料工程322 +14 哈哈哈吼吼吼哈 2026-04-07 15/750 2026-04-07 22:51 by 来看流星雨10
[考研] 081700，311，求调剂 +17 冬十三 2026-04-04 18/900 2026-04-07 12:50 by Sammy2
[考研] 一志愿华东理工085601材料工程303分求调剂 +8 a1708 2026-04-06 8/400 2026-04-07 11:20 by 诗与自由
[论文投稿] Decision: Revise for Editor还会送审吗 100+3 CccccccccFD 2026-04-04 5/250 2026-04-07 10:58 by 北京莱茵润色
[考研] 071000生物学调剂 +7 拉提桃 2026-04-06 7/350 2026-04-06 18:55 by 52305043001
[考研] 346分的生物与医药08600求调剂 +6 常雨阳上岸 2026-04-05 7/350 2026-04-06 12:36 by lys0704
[考研] 材料调剂 +14 壹贰贰亿 2026-04-04 14/700 2026-04-05 23:31 by 来看流星雨10
[考研] 262求调剂 +7 天下第一文 2026-04-04 8/400 2026-04-05 21:31 by 激流勇渡
[考研] 材料调剂 +9 革微桂 2026-04-04 9/450 2026-04-05 08:27 by 544594351
[考研] 本科211 分数293请求调剂 +4 莲菜就是藕吧 2026-04-01 4/200 2026-04-04 22:32 by hemengdong
[考研] 320求调剂 +3 一样圆 2026-04-04 3/150 2026-04-04 22:29 by 啵啵啵0119
[考研] 可跨专业调剂 +3 周的得地 2026-04-04 6/300 2026-04-04 22:21 by barlinike
[考研] 考研调剂 +5 小sun要好运 2026-04-03 5/250 2026-04-03 21:43 by 啵啵啵0119
[考研] 262求调剂 +6 励志一定发文章 2026-04-02 7/350 2026-04-03 09:54 by linyelide
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +8 双马尾痞老板2 2026-04-02 9/450 2026-04-02 14:45 by 5896
[硕博家园] 考研调剂 +5 骆驼男人 2026-04-01 5/250 2026-04-01 14:28 by syjjj0321
[考研] 一志愿北交材料工程总分358 +5 cs0106 2026-04-01 7/350 2026-04-01 11:45 by wangjy2002
信息提示
关闭
请填处理意见
关闭