| 5 | 1/1 | 返回列表 |
| 查看: 1767 | 回復(fù): 7 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
ShawnXu1986新蟲(chóng) (初入文壇)
|
[交流]
水分子模擬,如何解決由于正負(fù)電荷點(diǎn)過(guò)于靠近導(dǎo)致的平均動(dòng)能過(guò)大問(wèn)題?
|
||
|
各位前輩,大家好。 小弟自己寫(xiě)了個(gè)水分子模擬程序,NVE ensemble,使用了TIP3P, TIP5P水分子模型。在計(jì)算庫(kù)侖力時(shí)使用的是Ewald summation。 在這兩個(gè)模型中,正電荷位于H原子上,負(fù)電荷位于O原子上或者M site上。當(dāng)O 原子和H原子特別靠近時(shí),由于正負(fù)電荷吸引,庫(kù)侖力會(huì)使得這兩個(gè)原子靠近,但越靠近庫(kù)侖力越大,則兩個(gè)原子更加靠近! 由于LJ勢(shì)能只應(yīng)用在O-O 原子之間,所以沒(méi)有排斥力使得這一對(duì)OH原子分開(kāi),這兩種原子只能由極高的速度逃脫彼此,避免碰撞在一起。 但這就導(dǎo)致了整個(gè)系統(tǒng)的平均動(dòng)能變得非常大,使得模擬的溫度不合常理的高(大于10000K)! 我想到的一個(gè)辦法是在每一步都對(duì)速度進(jìn)行等比縮小,使得系統(tǒng)平均動(dòng)能給出的溫度能夠靠近我想要的溫度(比如300K)。但這個(gè)辦法并不可行。因?yàn)閮蓚(gè)電荷點(diǎn)避免碰撞在一起需要兩個(gè)電荷點(diǎn)有較高的速度,如果沒(méi)有速度,那么這兩個(gè)電荷就最終重疊到一起了,使得計(jì)算崩潰。 請(qǐng)各位能夠幫幫我,支出我的理解是否有誤,并請(qǐng)給我指出正確的方向,非常感謝 |
專(zhuān)家顧問(wèn) (正式寫(xiě)手)
![]() |
專(zhuān)家經(jīng)驗(yàn): +14 |
至尊木蟲(chóng) (文壇精英)
銀蟲(chóng) (初入文壇)
新蟲(chóng) (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料工程322分 +7 | 哈哈哈吼吼吼哈 2026-04-01 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 一志愿安徽大學(xué)計(jì)算機(jī)科學(xué)與技術(shù)學(xué)碩,331分求調(diào)劑 +3 | 蔣昌鵬qtj 2026-04-01 | 3/150 |
|
|
[考研] 295材料工程專(zhuān)碩求調(diào)劑 +19 | 1428151015 2026-03-27 | 19/950 |
|
|
[考研] 085602化學(xué)工程268分蹲調(diào)劑 +8 | 月照花林。 2026-04-01 | 8/400 |
|
|
[考研] 286求調(diào)劑 +5 | Sa67890. 2026-04-01 | 7/350 |
|
|
[考研] 290求調(diào)劑 +5 | dfffsar 2026-03-29 | 5/250 |
|
|
[考研] 26考研調(diào)劑 +3 | Wnz.20030617 2026-04-01 | 3/150 |
|
|
[考研] 309求調(diào)劑 +19 | 誰(shuí)不是少年 2026-03-29 | 19/950 |
|
|
[考研] 考研調(diào)劑 +11 | Amber00 2026-03-31 | 11/550 |
|
|
[考研] 生物學(xué)296求調(diào)劑 +10 | 湯圓包 2026-03-29 | 14/700 |
|
|
[考研] 318求調(diào)劑 +8 | 七憶77 2026-04-01 | 8/400 |
|
|
[考研] 考研材料工程351分調(diào)劑 +5 | 整個(gè)好的 2026-03-31 | 5/250 |
|
|
[考研] 085602 307分 求調(diào)劑 +10 | 不知道叫什么! 2026-03-26 | 10/500 |
|
|
[考研] 求化學(xué)調(diào)劑 +12 | wulanna 2026-03-28 | 12/600 |
|
|
[考研]
|
小羊36 2026-03-30 | 4/200 |
|
|
[考研] 367求調(diào)劑 +7 | 芋泥啵啵… 2026-03-28 | 7/350 |
|
|
[考研] 一志愿西電085401數(shù)一英一299求調(diào)劑 六級(jí)521 +4 | 愛(ài)吃大鴨梨 2026-03-31 | 4/200 |
|
|
[有機(jī)交流] 考研調(diào)劑 +8 | watb 2026-03-26 | 8/400 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動(dòng)285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
[考研]
|
平樂(lè)樂(lè)樂(lè) 2026-03-26 | 4/200 |
|