24小时热门版块排行榜    

查看: 768  |  回复: 5

digghost

木虫 (著名写手)

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

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

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的回帖

dahezhanjian

金虫 (小有名气)

我只是给你提供一个思路,其实还有更简单的方法,现在的计算量有点大,计算两要大至少10000倍,,你直接用这个程序算需要30000秒。你动动你的死脑筋吧,实在不行我再告诉你简单方法
5楼2013-03-06 17:05:44
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 digghost 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 309求调剂 +4 gajsj 2026-03-25 5/250 2026-03-26 00:27 by Dyhoer
[考研] 调剂310 +3 温柔的晚安 2026-03-25 4/200 2026-03-25 23:16 by peike
[考研] 07化学303求调剂 +5 睿08 2026-03-25 5/250 2026-03-25 22:46 by 418490947
[考研] 335分 | 材料与化工专硕 | GPA 4.07 | 有科研经历 +6 cccchenso 2026-03-23 6/300 2026-03-25 22:25 by 544594351
[考研] 求b区院校调剂 +4 周56 2026-03-24 5/250 2026-03-25 17:12 by yishunmin
[考研] 材料学硕333求调剂 +4 北道巷 2026-03-24 4/200 2026-03-25 14:16 by mapenggao
[考研] 各位老师您好:本人初试372分 +5 jj涌77 2026-03-25 6/300 2026-03-25 14:15 by mapenggao
[考研] 0854电子信息求调剂 324 +4 Promise-jyl 2026-03-23 4/200 2026-03-25 11:36 by Sugarlight
[考研] 材料学求调剂 +6 Stella_Yao 2026-03-20 6/300 2026-03-25 00:37 by baoball
[考研] 300分,材料,求调剂,英一数二 +5 超赞的 2026-03-24 5/250 2026-03-24 21:07 by 星空星月
[考研] 一志愿南航材料专317分求调剂 +5 炸呀炸呀炸薯条 2026-03-23 5/250 2026-03-24 16:52 by 星空星月
[考研] 一志愿山东大学药学学硕求调剂 +3 开开心心没烦恼 2026-03-23 4/200 2026-03-24 00:06 by 开开心心没烦恼
[考研] 求老师收我 +3 zzh16938784 2026-03-23 3/150 2026-03-23 12:56 by ztnimte
[考研] 070300,一志愿北航320求调剂 +3 Jerry0216 2026-03-22 5/250 2026-03-23 09:16 by 。。堂堂
[考研] 材料求调剂 +5 @taotao 2026-03-21 5/250 2026-03-21 20:55 by lbsjt
[考研] 材料 271求调剂 +5 展信悦_ 2026-03-21 5/250 2026-03-21 17:29 by 学员8dgXkO
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 288求调剂 +16 于海海海海 2026-03-19 16/800 2026-03-20 22:28 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
信息提示
请填处理意见