24小时热门版块排行榜    

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

digghost

木虫 (著名写手)

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

如题:,最好带文献
回复此楼
呵呵了
已阅   回复此楼   关注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的回帖
查看全部 6 个回答

dahezhanjian

金虫 (小有名气)

【答案】应助回帖

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

dahezhanjian

金虫 (小有名气)

我只是给你提供一个思路,其实还有更简单的方法,现在的计算量有点大,计算两要大至少10000倍,,你直接用这个程序算需要30000秒。你动动你的死脑筋吧,实在不行我再告诉你简单方法
5楼2013-03-06 17:05:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见