| 查看: 792 | 回复: 5 | ||
[求助]
如何根据晶格常数求lennard jones 势函数的两个参量
|
如题： ，最好带文献 |
回复此楼» 猜你喜欢
江西科技师范大学物理化学课题组接收调剂生，名额还有，速来
已经有0人回复
-
求碳酸钠和碳酸氢钠的晶胞结构图或晶体结构图
已经有0人回复
-
无机化学论文润色/翻译怎么收费?
已经有84人回复
-
求碳酸钠和碳酸氢钠的晶胞结构图或晶体结构图
已经有0人回复
-
江西科技师范大学物理化学课题组接收调剂生，12点开启，速来
已经有0人回复
-
【2026 考研调剂】哈尔滨工程大学 招收学硕调剂生（英语一、数学一）
已经有0人回复
-
四川轻化工大学化学与环境工程学院 付昱教授课题组招收调剂硕士研究生
已经有11人回复
-
大连化物所招收有机背景的博士候选人
已经有0人回复
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