24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁

1ST强帖: 1
应助: 2 (幼儿园)
贵宾: 0.08
金币: 921.3
散金: 149
帖子: 910
在线: 218.4小时
虫号: 366989
注册: 2007-05-10
性别: MM
专业: 凝聚态物性I:结构、力学和

[交流] 【求助】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

这样的固定方法正确。

回复此楼

1楼2010-01-19 09:36:21

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

yb2006

铁杆木虫 (正式写手)

应助: 6 (幼儿园)
金币: 8799.1
红花: 4
帖子: 754
在线: 335.9小时
虫号: 622347
注册: 2008-10-10
专业: 催化化学

★
minmin_0082003(金币+1):谢谢参与

你这样写是对的，但我通常是把要固定的原子坐标移了靠在一起，只写一个from...to...就行了。。。

赞一下(1人)

回复此楼

6楼2010-05-09 21:30:44

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

查看全部 10 个回答

zhangguangping

木虫 (著名写手)

1ST强帖: 23
应助: 71 (初中生)
贵宾: 0.031
金币: 1973.1
散金: 9059
红花: 79
帖子: 2681
在线: 2128.8小时
虫号: 529624
注册: 2008-03-20
性别: GG
专业: 原子和分子物理

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表示放开相应的自由度！这个你可以从手册找到！我用过这种固定方法，并且，还可以加限制！

赞一下(3人)

回复此楼

弘德明志博学笃行