24小时热门版块排行榜    

查看: 714  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

digghost

木虫 (著名写手)

[求助] 如何根据晶格常数求lennard jones 势函数的两个参量

如题:,最好带文献
回复此楼
呵呵了
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

digghost

木虫 (著名写手)

我自己弄出来了
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

dahezhanjian

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
digghost: 金币+20, ★★★★★最佳答案, 果然还是你厉害呀 2013-03-06 13:34:59
只能求指数的大小,把所有键能加起来得到一个指数的函数,取最小值对应的指数就行了
2楼2013-03-06 07:08:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

digghost

木虫 (著名写手)

引用回帖:
2楼: Originally posted by dahezhanjian at 2013-03-06 07:08:20
只能求指数的大小,把所有键能加起来得到一个指数的函数,取最小值对应的指数就行了

请问能提供一个程序么?
呵呵了
3楼2013-03-06 13:36:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dahezhanjian

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见