| 查看: 862 | 回复: 5 | ||
[求助]
如何根据晶格常数求lennard jones 势函数的两个参量
|
如题： ，最好带文献 |
回复此楼» 猜你喜欢
在MOF中，有没有羧基参加配位的情况
已经有2人回复
-
招2026年入学上海交通大学博士生，急！！！
已经有1人回复
-
无机化学论文润色/翻译怎么收费?
已经有84人回复
-
理论计算合作
已经有0人回复
-
求高校合作
已经有1人回复
-
27级无机化学申博自荐
已经有2人回复
-
Chemdraw免费软件
已经有0人回复
-
单晶结构解析，各类疑难晶体数据处理和答疑
已经有38人回复
-
这个水相反应有没有催化剂促进发生
已经有4人回复
-
Mn Me3TACN
已经有5人回复
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
28