24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 1336  |  回复: 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):谢谢您的参与,欢迎您常来计算模拟版!
所谓自由边界条件就是取无限大真空层 实际操作起来就是将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的回帖

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的回帖

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的回帖

honghuzlj

木虫 (著名写手)

谢谢你们,
Fortran程序中,原子N被排列在固定晶格上形成一个Box,像lei0736所说的10X10X10这样一个盒子,那么初始的Configuration形成;而且周期性边界条件是在计算力的程序中体现的,       例如x方向:xr=x(i)-x(j)
                                   xr=xr-box*nint(xr/box)     周期性边界条件
可是你说的“一般取50或者100是足够了 你要取1000当然更可以
那么Z方向是不是有很大的真空层?”,所以我不太明白?或者你的意思是说初始的configuration不是10X10X10而是10X10X50(或100,1000)?如果是这样也可以实现初始的configuration(不过感觉那只是盒子的形状问题啊,要不正方体,要不长方体),那么在计算力的时候,Fortran该如何编写呢?再次感谢
6楼2008-03-26 09:48:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

honghuzlj

木虫 (著名写手)

期待回复,谢谢
7楼2008-03-27 09:10:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

honghuzlj

木虫 (著名写手)

即是大牛又是好心人的何在?
8楼2008-03-28 09:56:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lei0736

荣誉版主 (职业作家)

优秀版主

★ ★
zzgyb(金币+2,VIP+0):谢谢你的参与,欢迎光临计算模拟版!
我的意思是模拟盒子可取10X10X50 而其中原子占据空间为 10X10X10 那么Z方向是不是就出来真空层了? 至于力的作用一样的算撒 没区别 这个不是你要考虑的 势的作用在一定范围内会截断
9楼2008-03-28 10:07:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

honghuzlj

木虫 (著名写手)

我明白你(lei0736)的意思了。

[ Last edited by honghuzlj on 2008-3-29 at 16:51 ]
10楼2008-03-28 11:02:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 honghuzlj 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程281还有调剂机会吗 +41 xaw. 2026-04-11 42/2100 2026-04-14 20:05 by honglizhao
[考研] 290调剂生物0860 +33 哇哈哈,。 2026-04-11 36/1800 2026-04-14 19:56 by fenglj492
[考研] 求调剂 +9 小聂爱学习 2026-04-11 13/650 2026-04-14 19:18 by Art1977
[考研] 211本科材料化工求调剂 +18 YHLAH 2026-04-11 22/1100 2026-04-14 17:52 by lhj2009
[考研] 材料相关专业344求调剂双非工科学校或课题组 +21 hualkop 2026-04-12 22/1100 2026-04-14 17:44 by lhj2009
[考研] 22408 312求调剂 +15 门路摸摸 2026-04-14 16/800 2026-04-14 17:15 by zs92450
[考研] 一志愿沪9,326求生物学调剂 +10 刘墨墨 2026-04-13 10/500 2026-04-14 15:16 by zs92450
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +6 zju2000 2026-04-14 17/850 2026-04-14 14:34 by jurkat.1640
[考研] 电气专硕320求调剂 +6 小麻子111 2026-04-10 6/300 2026-04-12 10:54 by lemon6009
[考研] 296求调剂 +14 汪!?! 2026-04-10 16/800 2026-04-12 10:48 by zhouyuwinner
[考研] 326求调剂 +6 Shansyn 2026-04-10 6/300 2026-04-12 09:46 by hammer3
[考研] 296求调剂 +6 汪!?! 2026-04-09 6/300 2026-04-11 11:25 by zhq0425
[考研] 085410-273求调剂 +6 X1999 2026-04-10 6/300 2026-04-11 10:32 by Delta2012
[考研] 311求调剂 +13 xyp想读书 2026-04-10 14/700 2026-04-11 09:41 by 猪会飞
[考研] 一志愿北理工298英一数二已上岸,感谢各位老师 +14 Reframe 2026-04-10 16/800 2026-04-10 23:07 by caotw2020
[考研] 一志愿中南大学物理学,英一66,求调剂 +4 长烟旖旎 2026-04-08 5/250 2026-04-10 10:31 by 颖果儿
[考研] 292求调剂 +9 笑笑袁 2026-04-09 9/450 2026-04-10 10:05 by LHGeng
[考研] 调剂 +12 月@163.com 2026-04-08 12/600 2026-04-09 14:27 by rl1980
[考研] 机械专硕273请求调剂 +6 庚申壬申 2026-04-07 6/300 2026-04-08 22:41 by bljnqdcc
[考研] 求调剂,现在还能填的 +3 上岸小莹加油 2026-04-08 3/150 2026-04-08 14:30 by zhq0425
信息提示
请填处理意见