应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » 【求助】自由边界条件
291/3返回列表123下一页
查看: 1357  |  回复: 28
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

honghuzlj

木虫 (著名写手)

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

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

[ Last edited by lei0736 on 2008-3-25 at 22:38 ]
回复此楼

» 猜你喜欢
1楼 2008-03-25 20:51:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lei0736

荣誉版主 (职业作家)

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

honghuzlj

木虫 (著名写手)
zzgyb(金币+0,VIP+0):不明白，应该说出你不明白的地方，这样别人就更有针对性的回答你
谢谢你，可是我还是不明白啊
赞一下(10人) 回复此楼
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方向的自由边界条件
够清楚不？
赞一下(20人) 回复此楼
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)则没有限制
赞一下(10人) 回复此楼
聪明人知道自己是傻子，傻子以为自己是聪明人
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方向是不是就出来真空层了？ 至于力的作用一样的算撒 没区别 这个不是你要考虑的 势的作用在一定范围内会截断
赞一下(10人) 回复此楼
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 的主题更新
291/3返回列表123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +3 yexuqing 2026-04-19 3/150 2026-04-20 11:29 by jurkat.1640
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考研] 26药学专硕105500求调剂 +7 喽哈加油 2026-04-13 8/400 2026-04-19 20:21 by Equinoxhua
[考研] 297，工科调剂? +11 河南农业大学-能 2026-04-14 11/550 2026-04-19 20:07 by Equinoxhua
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 291求调剂 +11 关忆北. 2026-04-14 11/550 2026-04-19 17:16 by 中豫男
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂推荐 +9 小聂爱学习 2026-04-14 9/450 2026-04-19 17:03 by 中豫男
[考研] 085404 22408 309分求调剂 +10 lzmk 2026-04-14 11/550 2026-04-19 16:42 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 307中医考研调剂 +9 于以采蘩 2026-04-14 9/450 2026-04-19 08:41 by 烟雨流涯
[考研] 0854求调剂 +23 门路摸摸 2026-04-15 27/1350 2026-04-19 01:59 by 烟雨流涯
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 一志愿中科大材料与化工，353分还有调剂学校吗 +10 否极泰来2026 2026-04-15 12/600 2026-04-17 17:54 by mapenggao
[考研] 295分求调剂 +5 ?要上岸? 2026-04-17 5/250 2026-04-17 16:51 by fenglj492
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
[考研] 一志愿沪9，生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
[考研] 一志愿A区211，22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[考研] 085600材料与化工329分求调剂 +24 叶zilin 2026-04-13 25/1250 2026-04-14 09:20 by 试管破裂
信息提示
关闭
请填处理意见
关闭