| 查看: 1675 | 回复: 6 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
[求助]
md分子模拟 简化单位问题 已有2人参与
|
|||||
|
硕士毕业论文用mathematica自己编程进行md对金刚石分子模拟。使用了tersoff势函数。程序写了大部分了，现在需要化简单位。查了很多资料基本上都是讲了lj势能的单位化简。请问tersoff怎么简化单位以便计算。 还有tersoff对金刚石的势能参数在哪里能查到呢？我找了几篇论文看到了一个简单的参数表但是教授说不对。有没有具体的书籍能够查阅？ 本人机械专业却选到了这么一个题目。希望研究分子模拟的大师们能多多指导。 发自小木虫IOS客户端 |
回复此楼» 猜你喜欢
310求调剂
已经有6人回复
-
面上5B能上会吗？
已经有5人回复
-
279求调剂
已经有3人回复
-
材料与化工272求调剂
已经有17人回复
-
材料工程专硕求调剂
已经有3人回复
-
337求调剂
已经有5人回复
-
材料与化工328分调剂
已经有9人回复
-
318一志愿吉林大学生物与医药 求调剂
已经有5人回复
-
一志愿南开大学0710生物学359求调剂
已经有4人回复
-
求调剂
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
brucefan
专家顾问 (著名写手)
-
专家经验: +218 - 模拟EPI: 4
- 应助: 422 (硕士)
- 金币: 2294.8
- 散金: 1720
- 红花: 88
- 帖子: 1423
- 在线: 670.2小时
- 虫号: 3048965
- 注册: 2014-03-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 分子模拟
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
abinitio: 金币+5, 感谢应助 2016-05-30 10:21:38
感谢参与，应助指数 +1
abinitio: 金币+5, 感谢应助 2016-05-30 10:21:38
|
这是我自己写分子动力学模拟程序时用的单位制，供参考： /*----------------------------------------------------------------------------80 The units used in this program are as follows: (1) The basic units are chosen to be (A) Mechanics -- energy: eV (electron volt) -- length: A (Angstrom) -- mass: amu (atomic mass unit) (B) Thermodynamics -- temperature: K (Kelvin) (C) Electricity -- charge: e (elementary charge) (2) Some derived units are (A) Mechanics -- time: A amu^(1/2) eV^(-1/2) ~ 10 fs -- velocity: eV^(1/2) amu^(-1/2) ~ 0.1 A fs^(-1) -- acceleration: eV A^(-1) amu^(-1) ~ 0.01 A fs^(-2) -- force: eV A^(-1) (B) Thermodynamics Boltzmann constant: K_B ~ 0.8625 * 10^(-4) eV K^(-1) (C) Electricity --electrostatic constant: K_C = 1 / (4 * PI * epsilon_0) ~ 14.4 eV A e^(-2) ------------------------------------------------------------------------------*/ 另外， (1) 关于tersoff势的参数，在Tersoff的1989年的文章中就有。 不知你是否觉得Tersoff原来的参数不好？ (2)半年时间显然不能掌握第一性原理计算，入门都难。 |
6楼2016-05-28 15:14:20
abinitio
版主 (著名写手)
- 模拟EPI: 2
- 应助: 221 (大学生)
- 贵宾: 0.093
- 金币: 7365.2
- 散金: 3917
- 红花: 57
- 帖子: 1317
- 在线: 485.2小时
- 虫号: 368545
- 注册: 2007-05-11
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
2楼2016-05-27 08:38:48
3楼2016-05-27 14:18:20
abinitio
版主 (著名写手)
- 模拟EPI: 2
- 应助: 221 (大学生)
- 贵宾: 0.093
- 金币: 7365.2
- 散金: 3917
- 红花: 57
- 帖子: 1317
- 在线: 485.2小时
- 虫号: 368545
- 注册: 2007-05-11
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 分子模拟
|
我觉得这个势函数不应该仅针对于石墨结构，金刚石结构也同样适用。如果有异议，可以用金刚石结构跑下看看是否稳定，或看看晶格参数。我感觉应该没问题。至于单位简化，我没用过，不了解。 发自小木虫Android客户端 |
» 本帖已获得的红花（最新10朵）
Masterarbeit4楼2016-05-27 15:35:44
