| 查看: 1583 | 回复: 6 | ||||
[求助]
md分子模拟 简化单位问题 已有2人参与
|
|
硕士毕业论文用mathematica自己编程进行md对金刚石分子模拟。使用了tersoff势函数。程序写了大部分了，现在需要化简单位。查了很多资料基本上都是讲了lj势能的单位化简。请问tersoff怎么简化单位以便计算。 还有tersoff对金刚石的势能参数在哪里能查到呢？我找了几篇论文看到了一个简单的参数表但是教授说不对。有没有具体的书籍能够查阅？ 本人机械专业却选到了这么一个题目。希望研究分子模拟的大师们能多多指导。 发自小木虫IOS客户端 |
回复此楼» 猜你喜欢
国自然申请面上模板最新2026版出了吗？
已经有12人回复
-
售中科院一区文章，我:8 O 5 5 1 O 5 4，含JCR-TOP
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有5人回复
-
溴的反应液脱色
已经有7人回复
-
推荐一本书
已经有12人回复
-
基金申报
已经有4人回复
-
计算机、0854电子信息（085401-058412）调剂
已经有4人回复
-
纳米粒子粒径的测量
已经有7人回复
-
常年博士招收(双一流，工科)
已经有4人回复
-
参与限项
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
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朵）
4楼2016-05-27 15:35:44
送红花一朵 |
ab initio 你好。我碰到一个问题想问问你的想法。 又找过一次教授。对于单位简化的问题估计他也没什么头绪。然后他让我用如你名字一样的ab initio计算。我之前三个月都在软件和精典力学的方法中寻求思路并且程序完成了大部分。现在基本上是要从头开始。我下个月合同到期，可以续签五个月。六个月的时间里你认为对于一个零基础的人来说是否有一定的困难。这两天大概翻了一下量子力学的书感觉要花长时间才能吸收。我的研究生专业是系统工程我在考虑要不要重新选择一个方向一致的题目。因为我感觉我就算是潜心研究第一性原理从头计算的方法，花大量的时间看了量子力学书对于以后的工作并没有什么帮助。 站在你的角度有什么看法呢。我是继续坚持下去，还是中止合同另寻出路呢。 发自小木虫IOS客户端 |
5楼2016-05-28 14:20:50
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:结构、力学和
- 管辖: 分子模拟
7楼2016-05-30 10:42:38