| 10 | 1/1 | 返回列表 |
| 查看: 1175 | 回復(fù): 9 | |||
Masterarbeit新蟲 (初入文壇)
|
[交流]
鉆石分子動態(tài)建模 已有2人參與
|
|
研究生畢業(yè)論文教授要求用mathematica進行分子動力學(xué)建模。 看了兩本書有了大概的思路。 疑惑的是 建模第一步計算立場 計算一個分子的力場書上有基本的計算力場的方法也有一些由經(jīng)驗來的比如mm amber 等簡化的立場。 但是鉆石是晶體結(jié)構(gòu) 沒有分子式。我應(yīng)該怎樣處理這個問題,怎樣來計算鉆石的力場。還有公式里所有的參數(shù)在哪里能查到? 希望前輩們能給提供幫助,初來乍到給沒有金幣可以貢獻。本人學(xué)的是機械系統(tǒng)工程教授給了我一個對于我來說全新的題目。希望能和大家多交流。 發(fā)自小木蟲IOS客戶端 |
新蟲 (初入文壇)
新蟲 (初入文壇)
版主 (著名寫手)
|
能否再多給點信息? 鉆石就是金剛石,一種結(jié)構(gòu)的C,標(biāo)準(zhǔn)的金剛石結(jié)構(gòu)是固定的,也很好建,如果單純建模也不需要力場,所以: 1、你的要求就只是建出模型就可以了?還是要有后續(xù)計算? 2、你所說的模型要求多大?比如多少個原子或尺寸? |
新蟲 (初入文壇)
送紅花一朵 |
謝謝你的回答,我是要用分子動力學(xué)來模擬金剛石內(nèi)部原子的運動。 我看了一些關(guān)于分子動力學(xué)建模的書。 第一步就是計算分子勢能也就是立場,然后用verlet來求解。其中用了周期邊界條件。 關(guān)于你補充的兩點: 1. 需要力場來進行后續(xù)模擬。我很困惑因為金剛石是一個晶體結(jié)構(gòu),無限的延伸下去,我要怎么計算才能算是一個分子內(nèi)部勢能,還是有什么簡化的方法? 2.按照周期邊界條件是取一個九宮格,每個里面放置一些分子使每個格子里的分子密度和材料的分子密度保持一致?墒俏疫@個不是一個分子啊,他是一個c一個c無限延伸下去的啊。 以上就是我很困惑的地方。謝謝指導(dǎo)。 發(fā)自小木蟲IOS客戶端 |
版主 (著名寫手)
|
我們一般是先建立一個小模型,比如8*8*8,每個單位是一個標(biāo)準(zhǔn)的金剛石結(jié)構(gòu),原子坐標(biāo)為 0 0 0 0 0.5 0.5 0.5 0 0.5 0.5 0.5 0 0.25 0.25 0.25 0.25 0.75 0.75 0.75 0.25 0.75 0.75 0.75 0.25 對于無限延展,是要用周期性邊界條件,這個概念和算法在所有的md說明中都有介紹。 至于勢函數(shù),一般是使用現(xiàn)成的,對于金剛石來說,常用的是sw和tersoff勢,具體你也自己查吧 |
新蟲 (初入文壇)
版主 (著名寫手)
|
每個單位是這樣的。我所說的8*8*8是指在xyz三個方向上重復(fù)8次,當(dāng)然你也可以多或者少 |
新蟲 (初入文壇)
送紅花一朵 |
明白你的意思了,弄懂了邊界條件的建立,比如我要模擬3.3.3個晶胞就應(yīng)該建立一個5.5.5的模型。我最近看了很多講tersoff的資料大多都是講了計算的原理和參數(shù)。對于怎樣建立坐標(biāo)系模擬的具體步驟并不清楚。因為我是要用mathmatica來建模程序都要自己編寫。所以想知道坐標(biāo)系怎么建立的。比如我要模擬一個5.5.5的金剛石模型大概有1000多個原子,每個原子的坐標(biāo)我都要自己計算出來嗎? 發(fā)自小木蟲IOS客戶端 |
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 石河子大學(xué)(211、雙一流)碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 尋找調(diào)劑 +4 | 倔強芒? 2026-03-21 | 4/200 |
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +3 | niko- 2026-03-22 | 3/150 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研]
|
.6lL 2026-03-18 | 8/400 |
|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +5 | 材化逐夢人 2026-03-18 | 5/250 |
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
[考研] 環(huán)境工程調(diào)劑 +9 | 大可digkids 2026-03-16 | 9/450 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|