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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » Siesta 胞形可变 怎么设置？？
81/1返回列表
查看: 1484  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

wyding

木虫 (正式写手)

[求助] Siesta 胞形可变 怎么设置？？

fdf文件有这一行，没有注释掉
Variable cell         T
在优化之后胞行没变，是不是还有别的设置
输出有应力的存在：
siesta: Stress tensor (static) (eV/Ang**3):  
siesta:    -0.071167    0.000803   -0.000176
siesta:     0.000803   -0.071355    0.001454
siesta:    -0.000175    0.001462   -0.114067

本来打算不管，可Z方向应力有点大，也就是说Z方向的晶格常数要优化，可优化起来比较麻烦。
敢问各位同行怎么解决这个问题的，小弟不胜感激！！！
回复此楼

» 收录本帖的淘帖专辑推荐

siesta问题

» 猜你喜欢
向强者看齐，即使高度不够也是昂首挺胸！
1楼 2012-07-12 15:10:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xhzha

铁杆木虫 (正式写手)
★ ★
fzx2008: 金币+2, 谢谢指教！ 2012-07-12 18:52:50
MD.MaxForceTol        0.0001 eV/Ang  
MD.MaxStressTol       0.01 GPa      
MD.VariableCell       T
加上这几条算出来的应力应该比较小，可能要注意下结构有没有变形，有些结构优化后对称性会有变化
赞一下 回复此楼
2楼2012-07-12 15:52:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wyding

木虫 (正式写手)
引用回帖:
2楼: Originally posted by xhzha at 2012-07-12 15:52:14
MD.MaxForceTol        0.0001 eV/Ang  
MD.MaxStressTol       0.01 GPa      
MD.VariableCell       T
加上这几条算出来的应力应该比较小，可能要注意下结构有没有变形，有些结构优化后对称性会有变化

茶壶Siesta不错嘛，呵呵！谢谢了！没想到是你回答的！我试试看吧！
赞一下 回复此楼
向强者看齐，即使高度不够也是昂首挺胸！
3楼2012-07-12 16:35:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

redskywei

木虫 (著名写手)

【答案】应助回帖

★
感谢参与，应助指数 +1
WDD880227: 金币+1, 多谢指教 2012-08-18 08:56:54
有一个简单的笨办法：z方向应力大（正值），那就手动增大z方向晶格参数，反之（应力为负值）减小相应晶格参数。三个方向同时优化太慢了点。
赞一下 回复此楼
坠入无边的网~~
4楼2012-07-15 10:54:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wyding

木虫 (正式写手)
引用回帖:
4楼: Originally posted by redskywei at 2012-07-15 10:54:36
有一个简单的笨办法：z方向应力大（正值），那就手动增大z方向晶格参数，反之（应力为负值）减小相应晶格参数。三个方向同时优化太慢了点。

你说反了，是应力正值减小晶格参数，应力负值增大晶格参数。
赞一下(1人) 回复此楼
向强者看齐，即使高度不够也是昂首挺胸！
5楼2012-08-17 20:06:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

redskywei

木虫 (著名写手)
引用回帖:
5楼: Originally posted by wyding at 2012-08-17 20:06:33
你说反了，是应力正值减小晶格参数，应力负值增大晶格参数。...

我很久以前做过和这样的测试，又找出原来的结果，确实记错了，但是你以计算结果为准就是，呵呵。那个调整使得趋近于零就是对的。
赞一下 回复此楼
坠入无边的网~~
6楼2012-08-17 22:23:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wyding

木虫 (正式写手)
引用回帖:
6楼: Originally posted by redskywei at 2012-08-17 22:23:30
我很久以前做过和这样的测试，又找出原来的结果，确实记错了，但是你以计算结果为准就是，呵呵。那个调整使得趋近于零就是对的。...

只是觉得这样是笨方法，以为可以智能的调胞的大小
赞一下 回复此楼
向强者看齐，即使高度不够也是昂首挺胸！
7楼2012-08-18 08:53:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

redskywei

木虫 (著名写手)
引用回帖:
7楼: Originally posted by wyding at 2012-08-18 08:53:41
只是觉得这样是笨方法，以为可以智能的调胞的大小...

确实这方法笨，但是其实做起来比智能的还快些，人的智力还是高于机器滴，
赞一下 回复此楼
坠入无边的网~~
8楼2012-08-19 10:31:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wyding 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 02:52 by l48hFd8Nwae2
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 02:52 by l48hFd8Nwae2
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 5sJFDKp2kwuI 2026-11-06 11/550 2026-11-09 02:48 by l48hFd8Nwae2
[公派出国] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +10 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 17:18 by 5sJFDKp2kwuI
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 16:43 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 16:23 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +7 ncTVBM6AOJHF 2026-11-04 17/850 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +7 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:47 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +6 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-06 13:47 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +8 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:36 by BkA9Bi5cMhkD
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 09:01 by 2lRNluKrJpwL
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 04:47 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +6 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
信息提示
关闭
请填处理意见
关闭