| 5 | 1/1 | 返回列表 |
| 查看: 1765 | 回復(fù): 7 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
ShawnXu1986新蟲 (初入文壇)
|
[交流]
水分子模擬,如何解決由于正負(fù)電荷點過于靠近導(dǎo)致的平均動能過大問題?
|
||
|
各位前輩,大家好。 小弟自己寫了個水分子模擬程序,NVE ensemble,使用了TIP3P, TIP5P水分子模型。在計算庫侖力時使用的是Ewald summation。 在這兩個模型中,正電荷位于H原子上,負(fù)電荷位于O原子上或者M site上。當(dāng)O 原子和H原子特別靠近時,由于正負(fù)電荷吸引,庫侖力會使得這兩個原子靠近,但越靠近庫侖力越大,則兩個原子更加靠近!! 由于LJ勢能只應(yīng)用在O-O 原子之間,所以沒有排斥力使得這一對OH原子分開,這兩種原子只能由極高的速度逃脫彼此,避免碰撞在一起。 但這就導(dǎo)致了整個系統(tǒng)的平均動能變得非常大,使得模擬的溫度不合常理的高(大于10000K)! 我想到的一個辦法是在每一步都對速度進行等比縮小,使得系統(tǒng)平均動能給出的溫度能夠靠近我想要的溫度(比如300K)。但這個辦法并不可行。因為兩個電荷點避免碰撞在一起需要兩個電荷點有較高的速度,如果沒有速度,那么這兩個電荷就最終重疊到一起了,使得計算崩潰。 請各位能夠幫幫我,支出我的理解是否有誤,并請給我指出正確的方向,非常感謝 |
至尊木蟲 (文壇精英)
銀蟲 (初入文壇)
新蟲 (初入文壇)
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +14 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 337求調(diào)劑 +4 | 《樹》 2026-03-29 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +8 | 柳溪邊 2026-03-26 | 8/400 |
|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +12 | c297914 2026-03-28 | 12/600 |
|
|
[考研] 277跪求調(diào)劑 +6 | 1915668 2026-03-27 | 10/500 |
|
|
[考研] 求收留 +5 | 1943443204 2026-03-28 | 5/250 |
|
|
[考研] 375求調(diào)劑 +4 | 雨夏整夜 2026-03-29 | 4/200 |
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
[考研] 0856求調(diào)劑 +7 | 楒桉 2026-03-28 | 7/350 |
|
|
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4 | IceGuo 2026-03-26 | 5/250 |
|
|
[考研] 332求調(diào)劑 +4 | @MZB382400 2026-03-28 | 4/200 |
|
|
[考研] 一志愿華理,數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 12/600 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
[考研]
|
18419759900 2026-03-25 | 8/400 |
|
|
[考博] 26申博 +3 | 加油沖! 2026-03-26 | 3/150 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-26 | 5/250 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|
|
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +4 | 17501029541 2026-03-23 | 6/300 |
|