 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3790 | 回復(fù): 8 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
15216105346金蟲 (著名寫手)
|
[求助]
為什么修改Dreiding力場形式用文獻的勢參數(shù)，體系勢能大的離譜！ 已有1人參與
|
||
|
我把Dreiding力場的默認L_J-12-6形式改成需要的Buckingham形式，勢參數(shù)選用某篇文獻，最后使用自己改進的力場對體系進行幾何優(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. 范德華項.jpg 庫倫項.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è): 固體無機化學(xué)
|
單位不會存在問題的，文獻是原子單位a.u，我已經(jīng)轉(zhuǎn)換成kcal/mol了。 剛用這個修改的力場跑完動力學(xué)，結(jié)果震驚了，密度未變（為晶體密度），溫度始終為0，壓強超大。。晶體熔融后原子盡然分層，幫我再利用這些信息分析一下錯在哪里？ ---- 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)驗: +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
pigrass
專家顧問 (著名寫手)
伊卡斯特的小豬
-
專家經(jīng)驗: +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
15216105346
金蟲 (著名寫手)
- 應(yīng)助: 13 (小學(xué)生)
- 金幣: 2062.4
- 散金: 2693
- 紅花: 7
- 帖子: 1205
- 在線: 298小時
- 蟲號: 3441372
- 注冊: 2014-09-25
- 專業(yè): 固體無機化學(xué)
|
只考慮Na、Al、F這三種元素。 謝謝��！ 勢參數(shù)，原子單位.gif 文獻參數(shù)形式.jpg |
» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責(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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
 一志愿北京科技材料科學(xué)與工程288分，求調(diào)劑
+13
|
是辰啊 2026-04-02 | 13/650 |
|
|---|---|---|---|---|
|
|
[考研] 310求調(diào)劑 +17 | 爭取九點睡 2026-03-30 | 17/850 |
|
|
|
[考研] 336求調(diào)劑 +6 | kiyy 2026-04-01 | 6/300 |
|
|
[考博] 材料工程專業(yè)碩士申博 +3 | 麟正宇 2026-03-30 | 3/150 |
|
|
|
[考研] 學(xué)碩機械工程303求調(diào)劑 +4 | 無名所以叫吳明 2026-03-30 | 4/200 |
|
|
|
[考研] 314求調(diào)劑 +11 | 1xiaojun23 2026-03-31 | 12/600 |
|
|
|
[考研] 311求調(diào)劑 +14 | 藍月亮亮 2026-03-30 | 14/700 |
|
|
|
[考研] 292求調(diào)劑 +14 | 是妍子也是研子 2026-03-30 | 14/700 |
|
|
|
[基金申請] esi高被引論文是不是能對中標(biāo)有所加分和幫助呢 +4 | redcom 2026-04-01 | 4/200 |
|
|
|
[考研] 266求調(diào)劑 +10 | 陽陽哇塞 2026-04-01 | 10/500 |
|
|
|
[考研] 303分 0807學(xué)碩求調(diào)劑 +3 | TYC3632 2026-04-01 | 3/150 |
|
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +7 | 阿姨復(fù)古不過 2026-03-27 | 7/350 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +9 | hyzzzzzzz. 2026-04-01 | 9/450 |
|
|
|
[考研] 080500-315分復(fù)試調(diào)劑 +9 | 上岸3821 2026-03-31 | 9/450 |
|
|
|
[考研] 求收留 +8 | 1943443204 2026-03-28 | 8/400 |
|
|
|
[考研] 一志愿西電085401數(shù)一英一299求調(diào)劑 六級521 +4 | 愛吃大鴨梨 2026-03-31 | 4/200 |
|
|
|
[考研]
 一志愿大連理工大學(xué)材料求調(diào)劑
+6
|
Gymno 2026-03-30 | 6/300 |
|
|
|
[考研] 11408總分309，一志愿東南大學(xué)求調(diào)劑，不挑專業(yè) +5 | 天賦帶到THU 2026-03-29 | 6/300 |
|
|
|
[考研] 085602 化學(xué)工程專碩 340分求調(diào)劑 +4 | qianbai11 2026-03-29 | 4/200 |
|
|
|
[考研] 356求調(diào)劑 +4 | gysy?s?a 2026-03-28 | 4/200 |
|