24小时热门版块排行榜    

查看: 1094  |  回复: 7

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问题

» 猜你喜欢

向强者看齐,即使高度不够也是昂首挺胸!
已阅   回复此楼   关注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方向晶格参数,反之(应力为负值)减小相应晶格参数。三个方向同时优化太慢了点。

你说反了,是应力正值减小晶格参数,应力负值增大晶格参数。
向强者看齐,即使高度不够也是昂首挺胸!
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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +3 大火山小火山 2026-03-16 5/250 2026-03-16 16:54 by barlinike
[考研] 070303一志愿西北大学学硕310找调剂 +5 d如愿上岸 2026-03-12 8/400 2026-03-16 15:19 by peike
[考研] 一志愿华中师范071000,325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +7 Ncdx123456 2026-03-13 8/400 2026-03-16 12:15 by karry wen
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 331求调剂(0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 0805,333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 332求调剂 +3 zjy101327 2026-03-11 6/300 2026-03-13 22:48 by JourneyLucky
[考研] 0703化学调剂 +4 快乐的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 310求调剂 +3 【上上签】 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
信息提示
请填处理意见