| 查看: 3155 | 回復(fù): 17 | ||||
qyshang1989新蟲 (小有名氣)
|
[交流]
請教DPD計(jì)算界面張力問題
|
|||
| 請大家指教:在MS中用DPD模擬出結(jié)果后,得到的結(jié)果都是簡化后的數(shù)據(jù),例如結(jié)果中的溫度,壓力,表面張力這些數(shù)據(jù)該如何轉(zhuǎn)換成現(xiàn)實(shí)中的物理量呢? |
xuexijisuan |
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗(yàn): +49 |
|
你要搞清楚是用哪個密度泛函理論(DFT)算表面張力。人們通常所說的DFT是指量子化學(xué)里面的DFT,也就是把能量表示為電子密度的泛函,據(jù)我了解量子化學(xué)里的DFT很難用來計(jì)算表面張力。MS和高斯所算的DFT都是指這種DFT。 另外一種DFT是指液體理論里的DFT,基于統(tǒng)計(jì)力學(xué)原理,是把自由能表示為粒子密度函數(shù)的泛函,可以計(jì)算體系很多熱力學(xué)性質(zhì),計(jì)算表面張力也是可行的。據(jù)我所了解很少有計(jì)算化學(xué)軟件包含了這種DFT的功能,很多課題組都有自己的代碼來進(jìn)行這方面的計(jì)算,你可以嘗試聯(lián)系相關(guān)課題組看別人是否愿意給代碼。 |
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗(yàn): +49 |
|
你這人怎么老是做伸手黨?自己搜一下material studio surface tension就會有東西出來。比如這個 http://accelrys.com/events/webin ... 20(NXPowerLite).pdf 參見第29張幻燈片。 或者參考這篇文獻(xiàn) http://pubs.acs.org/doi/abs/10.1021/la702286h Method部分中的公式。 簡而言之如果垂直于界面的方向定為z軸,需要計(jì)算的量就是壓力張量在xx, yy, zz 方向上的分量,然后表面張力 \gamma = (P_zz - 0.5*(P_xx + P_yy ) )*L_z L_z is the box length in z direction. 至于你怎么從MS的輸出中提取壓力張量在xx, yy, zz 方向上的分量請你自己研究吧,我沒用過MS,只能幫你到這里。 |
新蟲 (著名寫手)
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗(yàn): +49 |
|
"自由能密度泛函理論應(yīng)該就是您說的液體理論里的DFT吧?" 是。 我不清楚你到底想算什么,一會兒在MS里跑DPD,一會兒又想著DFT,一會兒又想用gromacs。 流體DFT不屬于分子模擬范疇,而是計(jì)算液體結(jié)構(gòu)的一種數(shù)值方法,所以gromacs這種分子模擬軟件是不包含流體DFT的相關(guān)計(jì)算的。網(wǎng)上有時看到的gromacs和DFT結(jié)合是指gromacs和量化計(jì)算中的DFT結(jié)合進(jìn)行QM/MM計(jì)算。 看你之前的帖子,你似乎是用MS做了個DPD模擬,想求表面張力,那你就應(yīng)該盡量想辦法從你現(xiàn)有的MS的輸出文件中求出表面張力,如果像樓上所說只是單位轉(zhuǎn)換的話,那難度很小。 另外使用別人的DFT代碼或者使用gromacs重新模擬計(jì)算表面張力都需要額外的時間去學(xué)習(xí),你要考慮你是否有這么多時間。如果你花了很長時間去學(xué)習(xí)另外一個軟件最后得到的效果和你處理MS數(shù)據(jù)得到的效果是一樣的,那只能說是得不償失。 |
新蟲 (小有名氣)
新蟲 (著名寫手)
新蟲 (小有名氣)
至尊木蟲 (文壇精英)

新蟲 (著名寫手)
新蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑 +7 | 妮妮ninicgb 2026-03-15 | 11/550 |
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +3 | FY_yy 2026-03-14 | 3/150 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 9/450 |
|
|
[教師之家] 焦慮 +8 | 水冰月月野兔 2026-03-13 | 12/600 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | pupcoco 2026-03-17 | 7/350 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 機(jī)械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 0703一志愿211 285分求調(diào)劑 +5 | ly3471z 2026-03-13 | 5/250 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
[考研] 求調(diào)劑 +5 | 一定有學(xué)上- 2026-03-12 | 5/250 |
|
|
[考研] 土木第一志愿276求調(diào)劑,科研和技能十分豐富,求新興方向的導(dǎo)師收留 +3 | 土木小天才 2026-03-12 | 3/150 |
|