 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3788 | 回復(fù): 8 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
15216105346金蟲 (著名寫手)
|
[求助]
為什么修改Dreiding力場形式用文獻(xiàn)的勢參數(shù)，體系勢能大的離譜！ 已有1人參與
|
||
|
我把Dreiding力場的默認(rèn)L_J-12-6形式改成需要的Buckingham形式，勢參數(shù)選用某篇文獻(xiàn)，最后使用自己改進(jìn)的力場對體系進(jìn)行幾何優(yōu)化， 發(fā)現(xiàn)體系無論在優(yōu)化前還是優(yōu)化后都大的離奇，71707286144820604000000000000000000000000000000000000000.000000 kcal/mol。 請問哪里錯了，是不是修改Dreiding力場過程中有什么問題？還是別的什么 Forcite ------- Task : Geometry Optimization Version : 7.0 Build date : Nov 12 2013 Host : xuzhenming Threads : Parallel 4 Operating system : Windows Task started : Tue Jan 27 23:43:32 2015 ---- Geometry optimization parameters ---- Algorithm : Smart Convergence tolerance: Energy : 0.001 kcal/mol Force : 0.5 kcal/mol/A Dreiding ---- Energy parameters ---- Forcefield : \xuzhenming（Dreiding的改動） Charges : Charge using QEq Electrostatic terms: Summation method : Ewald Accuracy : 0.001 kcal/mol Buffer width : 0.5 A van der Waals terms: Summation method : Atom based Truncation method : Cubic spline Cutoff distance : 12.5 A Spline width : 1 A Long range correction : YES Buffer width : 0.5 A ---- Initial structure ---- Total enthalpy : 71707286144820604000000000000000000000000000000000000000.000000 kcal/mol External pressure term : 0.000000 kcal/mol Total energy : 71707286144820600000000000000000000000000000000000000000.000000 kcal/mol Contributions to total energy (kcal/mol): Non-bond energy : 71707286144820600000000000000000000000000000000000000000.000 van der Waals : 20608469379.584 Long range correction : 71707286144820600000000000000000000000000000000000000000.000 Electrostatic : -1033129.589 rms force : 6.592E+004 kcal/mol/A max force : 2.059E+005 kcal/mol/A Cell parameters: a: 33.799800 A b: 33.763800 A c: 47.771400 A alpha: 90.000 deg beta: 90.000 deg gamma: 90.000 deg WARNING No progress. Optimization with current method stopped. ---- Final structure ---- Total enthalpy : 71707286144820604000000000000000000000000000000000000000.000000 kcal/mol External pressure term : 0.000000 kcal/mol Total energy : 71707286144820600000000000000000000000000000000000000000.000000 kcal/mol Contributions to total energy (kcal/mol): Non-bond energy : 71707286144820600000000000000000000000000000000000000000.000 van der Waals : 20029421699.475 Long range correction : 71707286144820600000000000000000000000000000000000000000.000 Electrostatic : 91604.397 rms force : 7.466E+004 kcal/mol/A max force : 4.661E+005 kcal/mol/A Cell parameters: a: 33.799800 A b: 33.763800 A c: 47.771400 A alpha: 90.000 deg beta: 90.000 deg gamma: 90.000 deg WARNING Convergence criteria are not satisfied. 范德華項(xiàng).jpg 庫倫項(xiàng).jpg 我的目標(biāo)勢函數(shù)形式.png |
回復(fù)此樓» 猜你喜歡
362求調(diào)劑
已經(jīng)有14人回復(fù)
-
求調(diào)劑 302分初試 0854
已經(jīng)有4人回復(fù)
-
299求調(diào)劑
已經(jīng)有4人回復(fù)
-
266分，一志愿電氣工程，本科材料，求材料專業(yè)調(diào)劑
已經(jīng)有3人回復(fù)
-
312求調(diào)劑
已經(jīng)有4人回復(fù)
-
315求調(diào)劑
已經(jīng)有7人回復(fù)
-
土木304求調(diào)劑
已經(jīng)有3人回復(fù)
-
316求調(diào)劑
已經(jīng)有16人回復(fù)
-
一志愿華東理工大學(xué)，080500學(xué)碩，317分，求調(diào)劑
已經(jīng)有13人回復(fù)
-
材料與化工306分找調(diào)劑
已經(jīng)有9人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- MS的力場文件函數(shù)如何調(diào)用
已經(jīng)有5人回復(fù)
- 關(guān)于MS力場編輯
已經(jīng)有11人回復(fù)
- 如何利用GULP擬合分子間勢能曲線得出力場參數(shù)
已經(jīng)有10人回復(fù)
- LJ勢怎么用，有圖片
已經(jīng)有6人回復(fù)
- ds中的分子優(yōu)化
已經(jīng)有5人回復(fù)
15216105346
金蟲 (著名寫手)
- 應(yīng)助: 13 (小學(xué)生)
- 金幣: 2062.4
- 散金: 2693
- 紅花: 7
- 帖子: 1205
- 在線: 298小時
- 蟲號: 3441372
- 注冊: 2014-09-25
- 專業(yè): 固體無機(jī)化學(xué)
|
只考慮Na、Al、F這三種元素。 謝謝��！ 勢參數(shù)，原子單位.gif 文獻(xiàn)參數(shù)形式.jpg |
» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : Structure_and_Raman_Spectra_in_Cryolitic_Melts_Simulations_with_an_ab_Initio_Interaction_Potential.pdf
2015-01-28 11:11:42, 1005.86 K
5樓2015-01-28 11:13:01
pigrass
專家顧問 (著名寫手)
伊卡斯特的小豬
-
專家經(jīng)驗(yàn): +503 - 模擬EPI: 4
- 應(yīng)助: 377 (碩士)
- 金幣: 10371
- 散金: 3286
- 紅花: 166
- 帖子: 2101
- 在線: 1968.4小時
- 蟲號: 1463681
- 注冊: 2011-10-27
- 性別: GG
- 專業(yè): 高分子材料結(jié)構(gòu)與性能
- 管轄: 分子模擬
2樓2015-01-28 08:47:24
15216105346
金蟲 (著名寫手)
- 應(yīng)助: 13 (小學(xué)生)
- 金幣: 2062.4
- 散金: 2693
- 紅花: 7
- 帖子: 1205
- 在線: 298小時
- 蟲號: 3441372
- 注冊: 2014-09-25
- 專業(yè): 固體無機(jī)化學(xué)
|
單位不會存在問題的，文獻(xiàn)是原子單位a.u，我已經(jīng)轉(zhuǎn)換成kcal/mol了。 剛用這個修改的力場跑完動力學(xué)，結(jié)果震驚了，密度未變（為晶體密度），溫度始終為0，壓強(qiáng)超大。。晶體熔融后原子盡然分層，幫我再利用這些信息分析一下錯在哪里？ ---- Dynamics parameters ---- Ensemble : NPT Temperature : 1300.00 K Control method : Nose Q ratio : 0.0100000 Pressure : 0.0001 GPa Control method : Berendsen Decay constant : 0.1000000 ps Timestep : 1.00 fs Number of steps : 100000 Duration : 100 ps Initial velocities : Random Random number seed : 1422376112 ---- Energy parameters ---- Forcefield : \xuzhenming Charges : Charge using QEq Electrostatic terms: Summation method : Ewald Accuracy : 0.001 kcal/mol Buffer width : 0.5 A van der Waals terms: Summation method : Atom based Truncation method : Cubic spline Cutoff distance : 12.5 A Spline width : 1 A Long range correction : YES Buffer width : 0.5 A ---- Thermostat ---- Nose mass parameter, Q mass: 334.727 kcal/mol*ps^2 ---- Dynamics summary ---- Initial Final Average Std. Dev. ------- ----- ------- --------- Tot. energy (kcal/mol) 71707286144820600000000000000000000000000000000000000000.000 71707286144820600000000000000000000000000000000000000000.000 71707286144820600000000000000000000000000000000000000000.000 0.000 Pot. energy (kcal/mol) 71707286144820600000000000000000000000000000000000000000.000 71707286144820600000000000000000000000000000000000000000.000 71707286144820600000000000000000000000000000000000000000.000 0.000 Kin. energy (kcal/mol) 16736.333 0.000 1674.217 523823.583 Tot. enthalpy (kcal/mol) 143414572289641230000000000000000000000000000000000000000.000 143414572289641230000000000000000000000000000000000000000.000 143414572289641230000000000000000000000000000000000000000.000 0.000 Temperature (K) 1300.000 0.000 130.045 40688.165 Pressure (GPa) 9138409243595302800000000000000000000000000000000000.000 9138409243595302800000000000000000000000000000000000.000 9138409243595302800000000000000000000000000000000000.000 0.000 Volume (A^3) 54517.184 54517.184 54517.184 0.000 Density (g/cm^3) 2.762 2.762 2.762 0.000 Task terminated : Wed Jan 28 10:24:16 2015 Total CPU time used by Forcite: 9:13:32 hours (33211.5s) Termination status : Normal 密度始終為晶體密度.jpg 溫度始終為0.jpg 原子分層.jpg |
3樓2015-01-28 10:55:29
pigrass
專家顧問 (著名寫手)
伊卡斯特的小豬
-
專家經(jīng)驗(yàn): +503 - 模擬EPI: 4
- 應(yīng)助: 377 (碩士)
- 金幣: 10371
- 散金: 3286
- 紅花: 166
- 帖子: 2101
- 在線: 1968.4小時
- 蟲號: 1463681
- 注冊: 2011-10-27
- 性別: GG
- 專業(yè): 高分子材料結(jié)構(gòu)與性能
- 管轄: 分子模擬
4樓2015-01-28 11:01:21
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 296材料專碩求調(diào)劑 +10 | 202451007219 2026-04-02 | 11/550 |
|
|---|---|---|---|---|
|
|
[考研] 26考研調(diào)劑 +4 | Wnz.20030617 2026-04-01 | 5/250 |
|
|
|
[考研] 材料專碩322分 +9 | 哈哈哈吼吼吼哈 2026-04-02 | 9/450 |
|
|
|
[考博] 材料工程專業(yè)碩士申博 +3 | 麟正宇 2026-03-30 | 3/150 |
|
|
[考研] 326求調(diào)劑 +3 | 崽崽仔 2026-04-02 | 3/150 |
|
|
|
[考研] 272求調(diào)劑，接受跨專業(yè)調(diào)劑！ +4 | 閑魚盧 2026-03-31 | 4/200 |
|
|
|
[考研] 324分 085600材料與化工 +20 | 呆鵝oor 2026-03-27 | 20/1000 |
|
|
|
[考研] 354求調(diào)劑 +4 | lxb598 2026-03-31 | 5/250 |
|
|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +14 | c297914 2026-03-28 | 14/700 |
|
|
|
[考研] 食品學(xué)碩362求調(diào)劑 +3 | xuanxianxian 2026-04-01 | 3/150 |
|
|
|
[考研] 求調(diào)劑0703 +5 | 周嘉堯 2026-03-31 | 8/400 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-30 | 16/800 |
|
|
|
[考研] 311求調(diào)劑 +10 | 李芷新1 2026-03-31 | 10/500 |
|
|
|
[考研] 一志愿北交材料工程總分358 +5 | cs0106 2026-04-01 | 7/350 |
|
|
|
[考研] 土木304求調(diào)劑 +5 | 頂級擦擦 2026-03-31 | 5/250 |
|
|
|
[考研] 一志愿西交大080500材料學(xué)碩349 +6 | jqx1258 2026-03-31 | 7/350 |
|
|
|
[考研] 375求調(diào)劑 +7 | 雨夏整夜 2026-03-29 | 7/350 |
|
|
|
[考研] 085600，專業(yè)課化工原理，320分求調(diào)劑 +6 | 大饞小子 2026-03-29 | 6/300 |
|
|
|
[考研]
 一志愿大連理工大學(xué)材料求調(diào)劑
+6
|
Gymno 2026-03-30 | 6/300 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|