24小时热门版块排行榜    

查看: 1047  |  回复: 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的回帖

redskywei

木虫 (著名写手)

引用回帖:
5楼: Originally posted by wyding at 2012-08-17 20:06:33
你说反了,是应力正值减小晶格参数,应力负值增大晶格参数。...

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

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的回帖
信息提示
请填处理意见