| 查看: 711 | 回复: 5 | |||
[求助]
如何根据晶格常数求lennard jones 势函数的两个参量
|
如题： ，最好带文献 |
回复此楼» 猜你喜欢
中南大学易小艺课题组诚招2026申请-考核制博士生
已经有0人回复
-
海南师范大学招收化学博士（光电功能材料课题组招收博士研究生）
已经有10人回复
-
无机化学论文润色/翻译怎么收费?
已经有239人回复
dahezhanjian
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 120.9
- 红花: 1
- 帖子: 84
- 在线: 122小时
- 虫号: 2280526
- 注册: 2013-02-13
- 专业: 催化化学
2楼2013-03-06 07:08:20
3楼2013-03-06 13:36:33
dahezhanjian
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 120.9
- 红花: 1
- 帖子: 84
- 在线: 122小时
- 虫号: 2280526
- 注册: 2013-02-13
- 专业: 催化化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
digghost: 金币+300, ★有帮助 2013-03-06 17:14:14
digghost: 金币+300, ★有帮助 2013-03-06 17:14:14
|
Rn[mx_, my_, mz_] := If[mx*my*mz == 0, 10000, (((mx)*ax)^2 + ((my)*ay)^2 + ((mz)*az)^2)^0.5; U[oab_, R_] := ((oab/R)^12 - (oab/R)^6)]; radio = 10 F[oab_] := Sum[U[oab, Rnm[mx, my, mz]], {mx, 0, radio}, {my, 0, radio}, {mz, 1, radio}] Solve[(F[oab])' == 0, oab] |
4楼2013-03-06 17:01:50
dahezhanjian
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 120.9
- 红花: 1
- 帖子: 84
- 在线: 122小时
- 虫号: 2280526
- 注册: 2013-02-13
- 专业: 催化化学
，实在不行我再告诉你简单方法
5楼2013-03-06 17:05:44
|
我自己弄出来了 Clear["Global`*"] radius = 15; a = 4.38/2 Rnm[mx_, my_, mz_] := ((mx*a)^2 + (my*a)^2 + (mz*a)^2)^0.5; Aa = (Sum[ If[Mod[(mx + my + mz), 2] == 0 && (mx != 0 || my != 0 || mz != 0), 1/(Rnm[mx, my, mz])^6, 0], {mx, 0, radius}, {my, 0, radius}, {mz, 0, radius}]/ Sum[If[Mod[(mx + my + mz), 2] == 0 && (mx != 0 || my != 0 || mz != 0), 1/(Rnm[mx, my, mz])^12, 0], {mx, 0, radius}, {my, 0, radius}, {mz, 0, radius}]/2)^(1/ 6) |
6楼2013-03-07 05:57:13