| 查看: 769 | 回复: 5 | |||
[求助]
如何根据晶格常数求lennard jones 势函数的两个参量
|
如题： ，最好带文献 |
回复此楼» 猜你喜欢
求助
已经有1人回复
-
0703化学26考研调剂，一志愿南昌大学
已经有4人回复
-
无机化学论文润色/翻译怎么收费?
已经有150人回复
-
有没有化学、材料专业的同学需要调剂 考虑天津高效的可以邮件或者私信
已经有1人回复
-
26年博士招生
已经有16人回复
-
江西理工大学功能晶态材料方向刘遂军课题组招收2026年秋季入学博士研究生
已经有11人回复
-
宁夏大学团簇新材料团队招收材料/化学/化工专业博士生
已经有0人回复
-
MIL-101有单晶样品吗
已经有1人回复
-
已删除
已经有3人回复
-
收调剂研究生
已经有2人回复
-
【2026 考研调剂】哈尔滨工程大学核科学与技术学院核化工系 招收调剂生
已经有2人回复
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
