24小时热门版块排行榜    

查看: 1607  |  回复: 28
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

honghuzlj

木虫 (著名写手)

[交流] 【求助】自由边界条件

请问自由边界条件怎么讲,比方说如果把z方向用自由边界条件,x.y方向用周期性边界条件,是否就是说,不把z方向定为周期性边界,那z自然就为自由边界?
否则,自由边界条件的Fortran代码如何编写?谢谢

[ Last edited by lei0736 on 2008-3-25 at 22:38 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lei0736

荣誉版主 (职业作家)

优秀版主

★ ★ ★ ★
zzgyb(金币+2,VIP+0):谢谢您的参与,欢迎您常来计算模拟版!
zzgyb(金币+2,VIP+0):谢谢您的参与,欢迎您常来计算模拟版!
晕 看来你的基本概念还需要加强
比如说原子占据的是10X10X10的立方盒子大小
那么对于模拟盒子 XY方向长度就取10 保证周期边界
Z方向取大于10 具体大多少看你的作用势范围
一般取50或者100是足够了 你要取1000当然更可以
那么Z方向是不是有很大的真空层?
这样就实现了Z方向的自由边界条件
够清楚不?
4楼2008-03-26 09:12:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 29 个回答

lei0736

荣誉版主 (职业作家)

优秀版主

★ ★
zzgyb(金币+2,VIP+0):谢谢您的参与,欢迎您常来计算模拟版!
所谓自由边界条件就是取无限大真空层 实际操作起来就是将Z方向的盒子取得比你的原子位置的最大最小值大很多或小很多 保证Z方向的表面不跟上下盒子相互作用即可
取原子位置最大值的100 ,1000倍都可以 如果作用力不是长程的话 取10倍足够了
2楼2008-03-25 22:38:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

honghuzlj

木虫 (著名写手)

zzgyb(金币+0,VIP+0):不明白,应该说出你不明白的地方,这样别人就更有针对性的回答你
谢谢你,可是我还是不明白啊
3楼2008-03-25 23:08:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cuihang

铁杆木虫 (正式写手)

VASPer

★ ★
lei0736(金币+2,VIP+0):谢谢 呵呵 我可不能这么提炼啊
用数学语言说,就是U(x)=U(x+L),U(y)=U(y+L),U(z)则没有限制
聪明人知道自己是傻子,傻子以为自己是聪明人
5楼2008-03-26 09:17:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 4/200 2026-11-02 00:50 by xzcv0Tehd8lQ
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-02 00:24 by xzcv0Tehd8lQ
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-02 00:08 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 19/950 2026-11-02 00:02 by xzcv0Tehd8lQ
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +7 X5bHxOULxFqb 2026-10-30 23/1150 2026-11-01 23:49 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:27 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[基金申请] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 22:01 by xzcv0Tehd8lQ
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 21:45 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 8/400 2026-11-01 21:35 by xzcv0Tehd8lQ
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 21:35 by xzcv0Tehd8lQ
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 21:19 by xzcv0Tehd8lQ
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 16:52 by BCdvBAEpt3zm
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:45 by BCdvBAEpt3zm
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 12:40 by BCdvBAEpt3zm
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 X5bHxOULxFqb 2026-10-30 26/1300 2026-11-01 10:22 by BCdvBAEpt3zm
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
信息提示
请填处理意见