版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

yongning88

銅蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 683.8
散金: 18
紅花: 1
帖子: 81
在線: 40.6小時(shí)
蟲號(hào): 2140323
注冊(cè): 2012-11-21
專業(yè): 金屬基復(fù)合材料

[求助] 勢(shì)函數(shù)的構(gòu)件已有1人參與

我按照文獻(xiàn)中的公式和參數(shù)計(jì)算了一下pair-potentiial, 算出的值和文獻(xiàn)中給出的圖一樣（如圖3所示），可為什么當(dāng) r 小于1.5 A 時(shí)，對(duì)勢(shì)就很大（如圖5所示），文獻(xiàn)中沒給出這部分的值，而LAMMPS中給出的勢(shì)文件卻是有個(gè)峰值的（如圖4所示）。跪求大牛們的幫助。

function.PNG

parmeters.PNG

pair_potential_in-paper.PNG

pair_potential_for_LAMMPS.PNG

pair_potenial_in_paper.PNG

回復(fù)此樓

» 猜你喜歡

292求調(diào)劑已經(jīng)有15人回復(fù)
304求調(diào)劑已經(jīng)有5人回復(fù)
材料與化工考研調(diào)劑已經(jīng)有17人回復(fù)
321求調(diào)劑已經(jīng)有5人回復(fù)
329求調(diào)劑已經(jīng)有5人回復(fù)
材料求調(diào)劑一志愿哈工大324 已經(jīng)有8人回復(fù)
286求調(diào)劑已經(jīng)有7人回復(fù)
083000學(xué)碩274求調(diào)劑已經(jīng)有7人回復(fù)
311求調(diào)劑已經(jīng)有8人回復(fù)
一志愿C9材料與化工專業(yè)總分300求調(diào)劑已經(jīng)有9人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

Si原子，O原子，C原子，H原子的L-J勢(shì)函數(shù)的參數(shù) 已經(jīng)有6人回復(fù)
非能量最小點(diǎn)的頻率計(jì)算(任意結(jié)構(gòu)的頻率計(jì)算) 已經(jīng)有14人回復(fù)
想用GULP擬合純鉬的morse勢(shì)函數(shù)參數(shù) 已經(jīng)有13人回復(fù)
Mg金屬結(jié)構(gòu)分子模擬已經(jīng)有5人回復(fù)
關(guān)于勢(shì)函數(shù)選取的問題。。。已經(jīng)有9人回復(fù)
分子模擬的關(guān)于勢(shì)函數(shù)以及初始構(gòu)型的幾個(gè)小問題已經(jīng)有15人回復(fù)
求助：哪里能找到各種物質(zhì)的勒納德瓊斯勢(shì)函數(shù) 已經(jīng)有5人回復(fù)
對(duì)于MD，其中的勢(shì)函數(shù)參數(shù)選取，如何自己計(jì)算已經(jīng)有5人回復(fù)
如何在lammps里加入新的勢(shì)函數(shù)形式？已經(jīng)有6人回復(fù)
Gulp 擬合勢(shì)函數(shù) 初值設(shè)置的問題已經(jīng)有12人回復(fù)
國(guó)內(nèi)有開發(fā)分子動(dòng)力學(xué)勢(shì)函數(shù)的課題組嗎？已經(jīng)有12人回復(fù)
關(guān)于勢(shì)函數(shù)的求助已經(jīng)有23人回復(fù)
日本開發(fā)出在磁場(chǎng)中變硬的橡膠新材料已經(jīng)有4人回復(fù)
Si的S-W勢(shì)函數(shù) 已經(jīng)有3人回復(fù)
【求助】?jī)?yōu)化的結(jié)構(gòu)非常奇怪，Au10 簇被優(yōu)化成了4組小簇，為什么，輸入及結(jié)構(gòu)已貼出已經(jīng)有12人回復(fù)
【求助】ms中如何假如勢(shì)函數(shù)，比如EAM，比如自己修改過的勢(shì)函數(shù)。已經(jīng)有12人回復(fù)
【討論】銅和水之間的LJ勢(shì)函數(shù)的參數(shù)是多少？已經(jīng)有19人回復(fù)

1樓 2014-04-29 09:51:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

老虎大王

木蟲 (著名寫手)

模擬EPI: 17
應(yīng)助: 26 (小學(xué)生)
貴賓: 0.17
金幣: 4774.1
散金: 8
紅花: 42
帖子: 1361
在線: 215.2小時(shí)
蟲號(hào): 659094
注冊(cè): 2008-11-21
專業(yè): 金屬結(jié)構(gòu)材料

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
月只藍(lán): 金幣+10, 感謝指導(dǎo)！ 2014-05-16 20:11:55

引用回帖:

17樓: Originally posted by yongning88 at 2014-05-14 21:28:31
謝謝您的回復(fù)，其他的都明白了，可現(xiàn)在還是不明白的是，如果勢(shì)函數(shù)不能取很大，也就是 r 的值不能太小，但構(gòu)件勢(shì)函數(shù)文件時(shí)，r得以零為起點(diǎn)，那對(duì)于這個(gè)部分的勢(shì)函數(shù)值怎么得到？...

這個(gè)，我想，我已經(jīng)說得比較清楚了�？墒悄愫孟襁€是不大清楚�，F(xiàn)在我給你一個(gè)例子，比如說，我要對(duì)某種立方晶體擬合L-J勢(shì)函數(shù)的參數(shù)： E=(A/r)^6-(B/r)^12，要擬合的就是A和B兩個(gè)參數(shù)。我選擇晶體的結(jié)合能Ec，晶胞體積V0，體彈性模量B，以及彈性系數(shù)C11和C12等來做為擬合目標(biāo)。擬合之后：
1）我獲得了A和B的值。現(xiàn)在我的勢(shì)函數(shù)就完整了。
2）現(xiàn)在我做E-r圖。做為一個(gè)數(shù)學(xué)函數(shù)，r 可以從零取到無窮大，都有對(duì)應(yīng)的E值。
3）但是對(duì)于我的模擬工作來說，r很小的時(shí)候，雖然有對(duì)應(yīng)的E值，但這些值在模擬上是用不到的。r很大的時(shí)候，同樣如此。
4）我發(fā)表這個(gè)勢(shì)函數(shù)的曲線圖的時(shí)候，我不需要從r=0開始畫圖，并且我要避免從r=0開始畫圖。因?yàn)閞很小的時(shí)候，勢(shì)能非常大，比如在30000以上，這樣我就需要把縱坐標(biāo)范圍定得很大，這就把后面實(shí)際上有用的部分掩蓋掉了，后面的E值的波動(dòng)在圖上就反映不出來了。
5）所以我畫圖的時(shí)候，要從某個(gè)r值開始畫，這個(gè)時(shí)候，E值還是比較大，但是能夠保證我后面的E值的波動(dòng)可以清楚地顯示。
6）這就是你看到的文獻(xiàn)的圖上，r不從零畫的原因。

總結(jié)一下，
我沒有說勢(shì)函數(shù)不能取很大，我是說兩點(diǎn)：
A) 作圖的時(shí)候，如果把r 很小的時(shí)候的勢(shì)函數(shù)曲線作出來，縱坐標(biāo)范圍就要定得很大，后面的E值的波動(dòng)在圖上就反映不出來了，而只有后面的E值才對(duì)于模擬是真正有用的。
B） r很小的時(shí)候，模擬上用不到，所以你根本不要去糾結(jié)r很小時(shí)候的勢(shì)能的情況，你計(jì)算勢(shì)函數(shù)值的時(shí)候，要把精力放在后面的部分。后面的部分和文獻(xiàn)一致了，你就算對(duì)了，后面的部分和文獻(xiàn)不一致，那么說明你和文獻(xiàn)作者之間至少有一方出錯(cuò)了。但是前面的r很小的部分的一致性是不需要考慮的。

贊一下(1人)

回復(fù)此樓

18樓2014-05-15 08:25:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

lsloneil

專家顧問 (正式寫手)

專家經(jīng)驗(yàn): +49
模擬EPI: 2
應(yīng)助: 148 (高中生)
金幣: 8640.9
散金: 507
紅花: 30
帖子: 505
在線: 598.4小時(shí)
蟲號(hào): 556471
注冊(cè): 2008-05-11
性別: GG
專業(yè): 理論和計(jì)算化學(xué)
管轄: 分子模擬

你是如何用lammps得到圖4和圖5的？

贊一下

回復(fù)此樓

2樓2014-04-29 10:31:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yongning88

銅蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 683.8
散金: 18
紅花: 1
帖子: 81
在線: 40.6小時(shí)
蟲號(hào): 2140323
注冊(cè): 2012-11-21
專業(yè): 金屬基復(fù)合材料

引用回帖:

2樓: Originally posted by lsloneil at 2014-04-29 10:31:08
你是如何用lammps得到圖4和圖5的？

這個(gè)不是用LAMMPS得到的，是用MATLAB 算的，我是想自己構(gòu)架那個(gè)可以在LAMMPS中用的勢(shì)函數(shù)

贊一下

回復(fù)此樓

3樓2014-04-29 10:37:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fanarsenal

木蟲 (正式寫手)

應(yīng)助: 47 (小學(xué)生)
金幣: 2023.2
紅花: 23
帖子: 515
在線: 659.5小時(shí)
蟲號(hào): 1069215
注冊(cè): 2010-08-04
專業(yè): 金屬結(jié)構(gòu)材料

★
月只藍(lán): 金幣+1, 鼓勵(lì)交流！ 2014-04-30 13:22:50

mark一下。
樓主的那張圖跟文獻(xiàn)一樣？我覺得從這個(gè)尺度看哪有一樣的，或者樓主標(biāo)清楚一下圖幾。。。
最好把文獻(xiàn)貼出來，帖子里面只放自己的圖。

贊一下

回復(fù)此樓

4樓2014-04-30 09:37:56

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

老虎大王

木蟲 (著名寫手)

模擬EPI: 17
應(yīng)助: 26 (小學(xué)生)
貴賓: 0.17
金幣: 4774.1
散金: 8
紅花: 42
帖子: 1361
在線: 215.2小時(shí)
蟲號(hào): 659094
注冊(cè): 2008-11-21
專業(yè): 金屬結(jié)構(gòu)材料

【答案】應(yīng)助回帖

★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
月只藍(lán): 金幣+3, 感謝指導(dǎo)！ 2014-04-30 11:06:35

原子離得近，當(dāng)然勢(shì)能就很大了，但是原子不可能離那么近的。像你這么做圖，把勢(shì)函數(shù)的細(xì)節(jié)全部掩蓋掉了。

贊一下

回復(fù)此樓

5樓2014-04-30 10:20:36

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yongning88

銅蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 683.8
散金: 18
紅花: 1
帖子: 81
在線: 40.6小時(shí)
蟲號(hào): 2140323
注冊(cè): 2012-11-21
專業(yè): 金屬基復(fù)合材料

引用回帖:

5樓: Originally posted by 老虎大王 at 2014-04-30 10:20:36
原子離得近，當(dāng)然勢(shì)能就很大了，但是原子不可能離那么近的。像你這么做圖，把勢(shì)函數(shù)的細(xì)節(jié)全部掩蓋掉了。

不好意思這幾天放假，沒及時(shí)回復(fù)，那對(duì)于原子很近的情況是怎么處理的呢？

贊一下

回復(fù)此樓

6樓2014-05-03 19:02:45

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

老虎大王

木蟲 (著名寫手)

模擬EPI: 17
應(yīng)助: 26 (小學(xué)生)
貴賓: 0.17
金幣: 4774.1
散金: 8
紅花: 42
帖子: 1361
在線: 215.2小時(shí)
蟲號(hào): 659094
注冊(cè): 2008-11-21
專業(yè): 金屬結(jié)構(gòu)材料

【答案】應(yīng)助回帖

★ ★ ★
月只藍(lán): 金幣+3, 感謝指導(dǎo)！ 2014-05-16 20:10:35

引用回帖:

6樓: Originally posted by yongning88 at 2014-05-03 19:02:45
不好意思這幾天放假，沒及時(shí)回復(fù)，那對(duì)于原子很近的情況是怎么處理的呢？...

原子離得很近的情況，比如r<1.0ang，那基本上是沒有什么實(shí)際意義的，原子不會(huì)靠得這么近，如果你模擬中原子真要是靠得這么近，那一定是弄錯(cuò)了。這一部分的勢(shì)函數(shù)不需要擬合，因?yàn)閷?shí)際上用不到。只要在r 比較小的的時(shí)候勢(shì)能足夠大，能夠阻止原子繼續(xù)靠近就可以了。你忽略掉這一部分就行了，不必糾結(jié)。主要是看后面的是否一致。

贊一下

回復(fù)此樓

7樓2014-05-05 09:13:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yongning88

銅蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 683.8
散金: 18
紅花: 1
帖子: 81
在線: 40.6小時(shí)
蟲號(hào): 2140323
注冊(cè): 2012-11-21
專業(yè): 金屬基復(fù)合材料

引用回帖:

7樓: Originally posted by 老虎大王 at 2014-05-05 09:13:09
原子離得很近的情況，比如r<1.0ang，那基本上是沒有什么實(shí)際意義的，原子不會(huì)靠得這么近，如果你模擬中原子真要是靠得這么近，那一定是弄錯(cuò)了。這一部分的勢(shì)函數(shù)不需要擬合，因?yàn)閷?shí)際上用不到。只要在r 比較小的 ...

謝謝您的解釋，既然r很小是可以忽略不用擬合，那為什么在有些已有的勢(shì)函數(shù)文件中，對(duì)勢(shì)值是從零開始增長(zhǎng)，直到出現(xiàn)峰值，他這有事怎么處理的呢？？

贊一下

回復(fù)此樓

8樓2014-05-05 09:42:33

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

老虎大王

木蟲 (著名寫手)

模擬EPI: 17
應(yīng)助: 26 (小學(xué)生)
貴賓: 0.17
金幣: 4774.1
散金: 8
紅花: 42
帖子: 1361
在線: 215.2小時(shí)
蟲號(hào): 659094
注冊(cè): 2008-11-21
專業(yè): 金屬結(jié)構(gòu)材料

★ ★ ★
月只藍(lán): 金幣+3, 感謝指導(dǎo)！ 2014-05-16 20:11:11

引用回帖:

8樓: Originally posted by yongning88 at 2014-05-05 09:42:33
謝謝您的解釋，既然r很小是可以忽略不用擬合，那為什么在有些已有的勢(shì)函數(shù)文件中，對(duì)勢(shì)值是從零開始增長(zhǎng)，直到出現(xiàn)峰值，他這有事怎么處理的呢？？...

他既然提出了一個(gè)函數(shù)，這個(gè)函數(shù)的參數(shù)擬合出來之后，自然會(huì)有一個(gè)函數(shù)值隨r變化的過程，r從零開始增長(zhǎng)，不久出現(xiàn)峰值，這是這個(gè)函數(shù)的性質(zhì)所決定的，是一個(gè)數(shù)學(xué)問題。而這個(gè)函數(shù)在r很小的時(shí)候的性質(zhì)，是不影響它的使用的。一個(gè)勢(shì)函數(shù)的可用性和勢(shì)函數(shù)的好壞，不是看這個(gè)函數(shù)在r很小的時(shí)候的情況，而是看這個(gè)函數(shù)在原子最近鄰距離的附近開始之后的情況，比如r＝1.0ang 或者1.5ang以后的情況。

贊一下

回復(fù)此樓

9樓2014-05-06 09:36:17

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yongning88

銅蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 683.8
散金: 18
紅花: 1
帖子: 81
在線: 40.6小時(shí)
蟲號(hào): 2140323
注冊(cè): 2012-11-21
專業(yè): 金屬基復(fù)合材料

引用回帖:

9樓: Originally posted by 老虎大王 at 2014-05-06 09:36:17
他既然提出了一個(gè)函數(shù)，這個(gè)函數(shù)的參數(shù)擬合出來之后，自然會(huì)有一個(gè)函數(shù)值隨r變化的過程，r從零開始增長(zhǎng)，不久出現(xiàn)峰值，這是這個(gè)函數(shù)的性質(zhì)所決定的，是一個(gè)數(shù)學(xué)問題。而這個(gè)函數(shù)在r很小的時(shí)候的性質(zhì)，是不影響它的 ...

對(duì)于我問題中給出的對(duì)勢(shì)公式和擬合出來的參數(shù)值，我用MATLAB從r=0開始算了下，如圖5所示，這個(gè)關(guān)于r的函數(shù)就沒有峰值，也就是從無窮大趨向于零。如果我要用這個(gè)公式來建立一個(gè)勢(shì)函數(shù)文件，需要從r=0，dr=0.0134A，截?cái)喟霃饺?.68A，也就是500個(gè)對(duì)勢(shì)值。選取r=0.67A時(shí)，對(duì)勢(shì)值為910eV，這個(gè)為最小的原子間距，也就是最大的對(duì)勢(shì)值，那對(duì)于r<0.67A的對(duì)勢(shì)值怎么選��？

贊一下

回復(fù)此樓

10樓2014-05-06 11:14:56

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 yongning88 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 一志愿廈門大學(xué)化學(xué)學(xué)碩307求調(diào)劑 +9	y7czhao 2026-03-26	9/450	2026-03-28 07:53 by Iveryant
[考研] 266分，求材料冶金能源化工等調(diào)劑 +6	哇呼哼呼哼 2026-03-27	8/400	2026-03-28 07:15 by mmm just
[考研] 330一志愿中國(guó)海洋大學(xué) 化學(xué)工程 085602 有讀博意愿求調(diào)劑 +3	wywy.. 2026-03-27	4/200	2026-03-28 03:32 by fmesaito
[考研] 277跪求調(diào)劑 +5	1915668 2026-03-27	8/400	2026-03-27 19:53 by WYUMater
[考研] 一志愿上海理工能源動(dòng)力（085800）310分求調(diào)劑 +3	zhangmingc 2026-03-27	4/200	2026-03-27 19:01 by 給你你注意休息
[考研] 07化學(xué)280分求調(diào)劑 +10	722865 2026-03-23	10/500	2026-03-27 15:51 by Plutoqq
[考研] 322求調(diào)劑 +4	舊吢 2026-03-24	4/200	2026-03-27 15:38 by 不吃魚的貓
[考研] 298調(diào)劑 +3	jiyingjie123 2026-03-27	3/150	2026-03-27 11:57 by wxiongid
[考研] 334求調(diào)劑 +3	雨清天晴 2026-03-21	3/150	2026-03-27 11:30 by 不吃魚的貓
[考研] 317求調(diào)劑 +7	蛋黃咸肉粽 2026-03-26	7/350	2026-03-27 02:29 by fmesaito
[考研] 349求調(diào)劑 +4	李木子啊哈哈 2026-03-25	4/200	2026-03-26 22:49 by fmesaito
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +4	柳溪邊 2026-03-26	4/200	2026-03-26 21:14 by XPU李慶
[考研] 352求調(diào)劑 +4	大米飯！ 2026-03-22	4/200	2026-03-26 16:40 by 不吃魚的貓
[考研] 資源與環(huán)境調(diào)劑申請(qǐng)(333分) +9	holy J 2026-03-21	9/450	2026-03-26 15:47 by 161765490
[考研] 07化學(xué)303求調(diào)劑 +5	睿08 2026-03-25	5/250	2026-03-25 22:46 by 418490947
[有機(jī)交流] 有機(jī)合成求助 20+3	FENGSHUJEI 2026-03-23	5/250	2026-03-24 19:31 by 88817753
[考研] 336化工調(diào)劑 +4	王大坦1 2026-03-23	5/250	2026-03-23 18:32 by allen-yin
[考研] 一志愿070300浙大化學(xué)358分，求調(diào)劑！ +4	酥酥魚.. 2026-03-21	4/200	2026-03-23 08:12 by Iveryant
[考研] 求調(diào)劑 +5	Zhangbod 2026-03-21	7/350	2026-03-22 13:13 by Zhangbod
[考研] 285求調(diào)劑 +6	ytter 2026-03-22	6/300	2026-03-22 12:09 by 星空星月

24小時(shí)熱門版塊排行榜

[求助] 勢(shì)函數(shù)的構(gòu)件 已有1人參與

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

[求助] 勢(shì)函數(shù)的構(gòu)件已有1人參與

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助: