24小时热门版块排行榜    

查看: 786  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

赵奇一

新虫 (小有名气)

[求助] 请教一下关于siesta优化层状体系的问题 已有2人参与

我建好了一个由两层原子构成的体系,想进行弛豫。希望在弛豫过程中,只发生两层原子间距的变化,而每层原子的相对位置不变,也就是说所有原子只在Z方向上坐标发生变化。这个如何实现?看了手册感觉还是不太明白,希望大神能够解答一下,谢谢了...
回复此楼
天使可以飞起来是因为他们把自己看得很轻
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

【答案】应助回帖


liliangfang: 金币+1, 谢谢交流 2015-11-22 16:12:55
%block Zmatrix                                                               
cartesian
1     -3.438229    6.506952    8.241026      0   0   0
1      2.329288    6.506952    8.241026      0   0   0
1      8.096804    1.512136    8.241026      0   0   0
1      6.654925    4.009544    8.241026      0   0   0
1      2.329288    1.512136    8.241026      0   0   0
1     -0.554470    6.506952    8.241026      0   0   0
1      5.213046    6.506952    8.241026      0   0   0
1      0.887409    9.004360    8.241026      0   0   0
1     -1.996350    9.004360    8.241026      0   0   0
1     -4.880108    9.004360    8.241026      0   0   0
1     -0.554470    1.512136    8.241026      0   0   0
1      5.213046    1.512136    8.241026      0   0   0
1      3.771167    9.004360    8.241026      0   0   0
1     -1.996350    4.009544    8.241026      0   0   0
1      3.771167    4.009544    8.241026      0   0   0
1      0.887409    4.009544    8.241026      0   0   0
2     -1.351749    6.858092    z1            0   0   1
2     -2.793628    9.355500    z1            0   0   1
2      4.415768    1.863277    z1            0   0   1
2      5.857647    9.355500    z1            0   0   1
2      2.973889    4.360684    z1            0   0   1
2      1.532010    6.858092    z1            0   0   1
2      7.299526    6.858092    z1            0   0   1
2      0.090131    9.355500    z1            0   0   1
2      7.299526    1.863277    z1            0   0   1
2      2.973889    9.355500    z1            0   0   1
2      5.857647    4.360684    z1            0   0   1
2      8.741405    4.360684    z1            0   0   1
2      4.415768    6.858092    z1            0   0   1
2      1.532010    1.863277    z1            0   0   1
2     10.183280    1.863277    z1            0   0   1
2      0.090131    4.360684    z1             0   0   1
variables
z1    37.271882
%endblock Zmatrix
上面这个zmatrix的示例就可以将两层(分别是1号原子和2号原子组成)的x,y坐标固定,在优化中只允许第二层的原子整体沿着z方向移动。第二层原子的z坐标的初值是37.271882
这样做的前提是这个双层结构平面要在xy平面内,两层间的距离在z方向上。
弘德明志博学笃行
4楼2015-11-22 15:50:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

charleslian

木虫 (小有名气)

【答案】应助回帖

★ ★
liliangfang: 金币+2, 谢谢交流 2015-11-22 08:31:37
使用zmatrix可以进行设置,详见其例子

发自小木虫IOS客户端
2楼2015-11-13 22:45:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

赵奇一

新虫 (小有名气)

引用回帖:
2楼: Originally posted by charleslian at 2015-11-13 22:45:16
使用zmatrix可以进行设置,详见其例子

谢谢您的回复,听说stress 1 2 4 5 6 这个关键词也可达到效果,是这样吗?因为我现在就是这么做的
天使可以飞起来是因为他们把自己看得很轻
3楼2015-11-14 19:46:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 澳大利亚 Murdoch University 全奖博士招生(3个名额)地质化工冶金领域 +51 AI8RaGaPaSCR 2026-08-07 54/2700 2026-08-14 03:15 by 51OlFEHbqXvR
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +9 Tide man 2026-08-10 10/500 2026-08-13 20:33 by 温淼淼淼
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +4 majunge000 2026-08-11 6/300 2026-08-13 16:50 by 四季常青藤
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
信息提示
请填处理意见