 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3791 | 回復: 8 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
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 我的目標勢函數(shù)形式.png |
回復此樓» 猜你喜歡
362求調劑
已經有14人回復
-
求調劑 302分初試 0854
已經有4人回復
-
299求調劑
已經有4人回復
-
266分，一志愿電氣工程，本科材料，求材料專業(yè)調劑
已經有3人回復
-
312求調劑
已經有4人回復
-
315求調劑
已經有7人回復
-
土木304求調劑
已經有3人回復
-
316求調劑
已經有16人回復
-
一志愿華東理工大學，080500學碩，317分，求調劑
已經有13人回復
-
材料與化工306分找調劑
已經有9人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- MS的力場文件函數(shù)如何調用
已經有5人回復
- 關于MS力場編輯
已經有11人回復
- 如何利用GULP擬合分子間勢能曲線得出力場參數(shù)
已經有10人回復
- LJ勢怎么用，有圖片
已經有6人回復
- ds中的分子優(yōu)化
已經有5人回復
15216105346
金蟲 (著名寫手)
- 應助: 13 (小學生)
- 金幣: 2062.4
- 散金: 2693
- 紅花: 7
- 帖子: 1205
- 在線: 298小時
- 蟲號: 3441372
- 注冊: 2014-09-25
- 專業(yè): 固體無機化學
7樓2015-01-28 18:49:18
pigrass
專家顧問 (著名寫手)
伊卡斯特的小豬
-
專家經驗: +503 - 模擬EPI: 4
- 應助: 377 (碩士)
- 金幣: 10371
- 散金: 3286
- 紅花: 166
- 帖子: 2101
- 在線: 1968.4小時
- 蟲號: 1463681
- 注冊: 2011-10-27
- 性別: GG
- 專業(yè): 高分子材料結構與性能
- 管轄: 分子模擬
2樓2015-01-28 08:47:24
15216105346
金蟲 (著名寫手)
- 應助: 13 (小學生)
- 金幣: 2062.4
- 散金: 2693
- 紅花: 7
- 帖子: 1205
- 在線: 298小時
- 蟲號: 3441372
- 注冊: 2014-09-25
- 專業(yè): 固體無機化學
|
單位不會存在問題的，文獻是原子單位a.u，我已經轉換成kcal/mol了。 剛用這個修改的力場跑完動力學，結果震驚了，密度未變（為晶體密度），溫度始終為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
專家顧問 (著名寫手)
伊卡斯特的小豬
-
專家經驗: +503 - 模擬EPI: 4
- 應助: 377 (碩士)
- 金幣: 10371
- 散金: 3286
- 紅花: 166
- 帖子: 2101
- 在線: 1968.4小時
- 蟲號: 1463681
- 注冊: 2011-10-27
- 性別: GG
- 專業(yè): 高分子材料結構與性能
- 管轄: 分子模擬
4樓2015-01-28 11:01:21
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 土木304求調劑 +3 | 兔突突突， 2026-04-02 | 3/150 |
|
|---|---|---|---|---|
|
[考研]
 一志愿北京科技材料科學與工程288分，求調劑
+13
|
是辰啊 2026-04-02 | 13/650 |
|
|
|
[考研] 085601一志愿北理325分求調劑 +3 | 找調劑，， 2026-04-02 | 3/150 |
|
|
|
[考研] 310求調劑 +17 | 爭取九點睡 2026-03-30 | 17/850 |
|
|
|
[考研] 288求調劑 一志愿哈工大 材料與化工 +30 | 洛神哥哥 2026-03-31 | 32/1600 |
|
|
|
[考研] 材料專業(yè)求調劑 +10 | 月月鳥木 2026-04-01 | 10/500 |
|
|
|
[考研] 0856初試324分求調劑 +6 | 想上學求調 2026-04-01 | 6/300 |
|
|
|
[考研] 08生物與醫(yī)藥專碩初試346找調劑 +6 | dianeeee 2026-04-01 | 7/350 |
|
|
|
[考研] 江蘇科技大學招材料研究生 +4 | Su032713. 2026-04-01 | 5/250 |
|
|
|
[考研] 材料專碩322分 +9 | 哈哈哈吼吼吼哈 2026-04-01 | 9/450 |
|
|
|
[考研] 086000生物與醫(yī)藥298調劑求助 +4 | 元元青青 2026-03-31 | 6/300 |
|
|
|
[考研] 一志愿武理材料工程302調劑環(huán)化或化工 +15 | Doleres 2026-03-31 | 16/800 |
|
|
|
[考研] 求調劑 生物學 377分 +6 | zzll03 2026-03-31 | 6/300 |
|
|
|
[考研] 江蘇蘇北高校誠邀調劑同學 +3 | zzll406 2026-03-31 | 3/150 |
|
|
|
[考研] 343求調劑 +8 | 愛羈絆 2026-03-28 | 8/400 |
|
|
|
[考研] 本科211安全工程，初試290分，求調劑 +3 | 2719846834 2026-03-28 | 3/150 |
|
|
|
[考研] 085600材料與化工調劑 +16 | kikiki7 2026-03-30 | 16/800 |
|
|
|
[考研] 福建理工大學材料學院先進合金團隊招收考研調劑學生 +3 | 大華金商都 2026-03-30 | 4/200 |
|
|
|
[考研] 348求調劑 +6 | 小懶蟲不懶了 2026-03-28 | 6/300 |
|
|
|
[考研] 272求調劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|