| 查看: 637 | 回复: 0 | ||
[求助]
流体力学气固界面作用，在建立固体作用势时遇到了问题，弄过碳材料的各位来帮咱看看
|
|
本人专业是流体力学，因为要做气固相互作用所以才开始看固体的书，现在还处于扫盲阶段，但是一方面工作也必须展开，所以先在这就一些基础问题向各位请教。 因为暂时不做反应、烧蚀等内容，所以涉及的物理过程比较简单，MD程序是自己编的，气体和气固都是查文献取的LJ势或Morse势。但是原先为了跑通程序，固体是假想的FCC结构，自己随便给的LJ势。但现在要算个可以对比验证的算例，且材料是碳，所以取石墨作为对象，由于不反应不烧蚀，且气固相互作用力不大，不会产生太大的变形或破坏，所以按文献取C-C的Morse势，里面仅包含键能和键角能，如下图：  那么问题来了，一般不是说石墨的键长是1.42嘛，但是这个势里r0却给的是1.39，请问是不是要把1.39改成1.42才对？ 第二个问题，上面是石墨的键能，但是远处的碳原子或者石墨层间应该还存在范德华力吧？百度到的碳原子间范德华力是L-J势： U =B/r^12- A/r^6 （对于2 个碳原子间，其参数值为B =11.5 ×10-6 kJnm^12/mol ；A=5.96 × 10-3 kJnm^6/mol） 上式开导算出来的平衡距离是3.96，但是石墨层间距不是3.35吗？个人感觉就算不是完全堆积，错位堆积也不应该差这么多，所以我还是怀疑势参数给的不对。如果不是势参数的问题，那么就是势本身的问题？然后继续查文献提到石墨层间相互作用叫什么PI-PI相互作用，是一种键力，但有些资料又说PI-PI相互作用就是范德华力，我现在彻底晕了，到底石墨层间相互作用是什么？恳切希望弄过碳材料的大伙能给说个明白，或者直接上个表达式也行  C-C Morse.JPG |
回复此楼» 猜你喜欢
咱们一起用铁证分析2026国家社科基金中标与否
已经有27人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有3人回复
-
2027广东省杰青
已经有9人回复
-
朋友圈看到的
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有110人回复
-
明天放榜？
已经有5人回复
-
时间戳变了，能看出什么问题？
已经有11人回复
-
重要消息，中午系统在维护
已经有10人回复
-
投稿咨询
已经有6人回复
-
时间戳又变了8-15
已经有26人回复
» 本主题相关商家推荐: (我也要在这里推广)
