24小时热门版块排行榜    

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

minmin_0082003

金虫 (正式写手)

莫愁莫愁

[交流] 【求助】siesta如何让固定衬底? 已有7人参与

在计算多层的surface的时候,常常会固定除表层外的其他层面,请问在固定的过程中是,直接固定每层的每个原子就可以吗?
比如说,我一共要算4层的,但是要固定底下面的3层。我把每个原子都固定住了,如下:
%block GeometryConstraints
Position 1
position from 3 to 5
position from 5 to 9
position from 11 to 13
position from 15 to 17
position from 19 to 21
position from 23 to 25
position from 27 to 29
position from 31 to 33
position from 35 to 36
%endblock GeometryConstraints

这样的固定方法正确。
回复此楼

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

siesta问题

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

★ ★
qasd(金币+2):xiexie 2010-04-10 20:00
你可以参考一下seista-3.0-b的manual的第33-34页!
弘德明志博学笃行
3楼2010-03-24 13:55:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

zhangguangping

木虫 (著名写手)

siesta如何让固定衬底?

★ ★ ★
minmin_0082003(金币+1):谢谢参与
ice_rain(金币+2):xiexie!~ 2010-03-24 21:42
你的这种方法可行@
还有一种方法,就是使用zmtrix,如下:
%block Zmatrix                                                               
cartesian                                                                     
4     0.832500     1.441900     0.070000    0    0    0                       
4     0.832500     4.325700     0.070000    0    0    0                       
4     0.832500     7.209500     0.070000    0    0    0                       
4     3.329900     0.000000     0.070000    0    0    0                       
4     3.329900     2.883800     0.070000    0    0    0                       
4     3.329900     5.767600     0.070000    0    0    0                       
4     5.827400    -1.441900     0.070000    0    0    0                       
4     5.827400     1.441900     0.070000    0    0    0                       
4     5.827400     4.325700     0.070000    0    0    0                       
1     3.602290     1.731278     6.212320    1    1    1                     
1     3.319832     2.314888     7.483680    1    1    1                       
1     2.095981     1.995008     8.160953    1    1    1                       
1     1.179880     1.055287     7.579063    1    1    1                       
1     1.463109     0.470284     6.307043    1    1    1                       
1     2.686364     0.792157     5.629784    1    1    1                       
2     4.560824     2.006388     5.644557    1    1    1                       
2     4.060912     3.044923     7.962504    1    1    1                       
2     0.218799     0.778868     8.141618    1    1    1                       
2     0.720051    -0.260182     5.831012    1    1    1                       
2     2.264637     3.524194     9.544223    1    1    1                       
2     0.774280     2.665518     9.602215    1    1    1                       
2     4.021096     0.084236     4.221899    1    1    1                       
2     2.491712    -0.710195     4.216641    1    1    1                       
3     1.840356     2.531978     9.498710    0    0    1                       
3     2.958183     0.262423     4.292613    0    0    1                     
4     2.438337     1.576108     2.385995    0    0    1                       
4     2.578104     1.329674    11.405459    0    0    1                       
4     1.665000     0.000000    13.722500    0    0    0                     
4     1.665000     2.883800    13.722500    0    0    0                       
4     1.665000     5.767600    13.722500    0    0    0                     
4     4.162400    -1.441900    13.722500    0    0    0                     
4     4.162400     1.441900    13.722500    0    0    0                     
4     4.162400     4.325700    13.722500    0    0    0                       
4     6.659800    -2.883800    13.722500    0    0    0                     
4     6.659800     0.000000    13.722500    0    0    0                       
4     6.659800     2.883800    13.722500    0    0    0                       
%endblock Zmatrix            
第一列是元素序号,后面三列是相应的原子的cartesian坐标,后面的0表示固定相应的自由度,1表示放开相应的自由度!这个你可以从手册找到!我用过这种固定方法,并且,还可以加限制!
弘德明志博学笃行
2楼2010-03-24 11:33:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lionz

银虫 (初入文坛)

minmin_0082003(金币+1):谢谢参与
aylayl08(金币-1):禁止无意义回帖,扣回所得 2010-05-09 18:33:25
Thanks Zhangguangping
4楼2010-04-09 10:38:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yb2006

铁杆木虫 (正式写手)

★
minmin_0082003(金币+1):谢谢参与
你这样写是对的,但我通常是把要固定的原子坐标移了靠在一起,只写一个from...to...就行了。。。
6楼2010-05-09 21:30:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯,如何抉择 +3 Isak0320 2026-10-12 5/250 2026-10-12 16:42 by Isak0320
[教师之家] 课题组会帮你通过考核么? +3 好运来~9527 2026-10-10 6/300 2026-10-12 14:55 by 简单化xn
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +4 wangyl_123 2026-10-10 7/350 2026-10-12 14:47 by chuans0923
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +5 13108017953 2026-10-08 5/250 2026-10-10 11:15 by wdwd123
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见