| 1 | 1/1 | 返回列表 |
| 查看: 1071 | 回復(fù): 0 | |||
A寶祿計(jì)算銀蟲 (初入文壇)
|
[交流]
分子動(dòng)力學(xué)計(jì)算
|
|
要運(yùn)行基本的分子動(dòng)力學(xué)計(jì)算,請(qǐng)執(zhí)行以下步驟: 選擇包含足夠大的超級(jí)電池的POSCAR 。 如果執(zhí)行連續(xù)運(yùn)行,則將CONTCAR復(fù)制到POSCAR或可能在POSCAR文件中提供初始速度。它們是在 Wycoff 立場之后的一個(gè)單獨(dú)段落中寫成的。如果沒有提供初始速度,則在計(jì)算開始時(shí)假定為隨機(jī)速度。這完全沒問題,但用戶應(yīng)該知道,由于初始隨機(jī)速度,從不同計(jì)算獲得的軌跡很難比較。 設(shè)置主要INCAR標(biāo)簽: oIBRION =0:通過將IBRION標(biāo)簽設(shè)置為 0啟用分子動(dòng)力學(xué)計(jì)算。 oPOTIM:此標(biāo)簽設(shè)置 fs 中分子動(dòng)力學(xué)運(yùn)行的時(shí)間步長。 oNSW:此標(biāo)簽設(shè)置分子動(dòng)力學(xué)運(yùn)行中執(zhí)行的步驟數(shù)。 oTEBEG:如果使用恒溫器,請(qǐng)定義分子動(dòng)力學(xué)計(jì)算應(yīng)運(yùn)行的所需溫度。 oISIF(可選)。 oMDALGO:此標(biāo)簽決定使用哪個(gè)恒溫器執(zhí)行分子動(dòng)力學(xué)計(jì)算。對(duì)于常規(guī)分子動(dòng)力學(xué)計(jì)算,恒溫器由一位數(shù)字選擇(例如,1 代表 Andersen,2 代表 Nose-Hoover 等)。對(duì)于偏置分子動(dòng)力學(xué)、元?jiǎng)恿W(xué)等,恒溫器的選擇方式與 VASP 6 或更高版本相同。在 VASP 5.x 中,它由兩位數(shù)字選擇,其中第一個(gè)數(shù)字對(duì)應(yīng)于類似于常規(guī)分子動(dòng)力學(xué)的恒溫器,第二個(gè)數(shù)字對(duì)應(yīng)于分子動(dòng)力學(xué)類型(例如,Andersen 恒溫器的 11 元?jiǎng)恿W(xué),Nose-Hoover 的 21 元?jiǎng)恿W(xué)溫控器等)。NVE系綜是一個(gè)特例。通過選擇 Andersen 恒溫器并設(shè)置不與熱浴發(fā)生碰撞(ANDERSEN_PROB =0)。 oISIF:在分子動(dòng)力學(xué)計(jì)算中,此標(biāo)簽用于選擇NVT 系綜或NpT 系綜(NVE 系綜是一個(gè)特例。。對(duì)于ISIF =2,體積保持不變并使用NVT 系綜。使用此標(biāo)簽可以計(jì)算應(yīng)力張量,因此可以監(jiān)控壓力。對(duì)于ISIF =3,應(yīng)力張量(壓力)保持恒定并使用NpT 系綜。使用此標(biāo)簽可以計(jì)算體積并對(duì)其進(jìn)行監(jiān)控。 決定使用 哪個(gè)集成: oNVT 合奏:設(shè)置ISIF = 2。 oNpT 系綜:設(shè)置ISIF =3。 oNVE 集成:設(shè)置MDALGO =1 和ANDERSEN_PROB =0.0。 決定使用哪個(gè)恒溫器(恒溫器和整體的組合在表中給出): oAndersen 恒溫器:設(shè)置MDALGO =1。還設(shè)置ANDERSEN_PROB >0.0 以控制恒溫器的隨機(jī)更新頻率。 oNose-Hoover 恒溫器:設(shè)置MDALGO =2。還設(shè)置SMASS >0.0 以控制與熱浴的耦合。 oLangevin 恒溫器:設(shè)置MDALGO =3。同時(shí)設(shè)置LANGEVIN_GAMMA >0.0 來控制摩擦參數(shù)。如果使用NpT 系綜(通過設(shè)置ISIF =3),還必須提供晶格LANGEVIN_GAMMA_L的摩擦系數(shù)。 以下恒溫器和恒壓器的組合是可能的: Thermostat Ensemble Andersen Nose-Hoover Langevin Multiple Andersen NVE MDALGO =1, ANDERSEN_PROB =0.0 NVT MDALGO =1 MDALGO =2 MDALGO =3 MDALGO =13 ISIF =2 ISIF =2 ISIF =2 ISIF =2 NpT not available not available MDALGO =3 not available ISIF =3 NpH MDALGO =3,ISIF =3,LANGEVIN_GAMMA_L =0.0 匯編 要運(yùn)行分子動(dòng)力學(xué)計(jì)算,必須使用makefile.include 文件中的-Dtbdyn預(yù)編譯器標(biāo)志編譯 VASP。使用此標(biāo)簽的示例輸入如下所示: CPP = $(CPP_) -DHOST=\"IFC9_fftw\" \ -Dkind8 -DNGXhalf -DCACHE_SIZE=12000 -DPGF90 -Davoidalloc \ -Dtbdyn |
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工304求調(diào)劑 +3 | 邱gl 2026-03-05 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-04 | 3/150 |
|
|
[考研] 320材料一志愿西工大專碩英二數(shù)二 有兩年光伏工作經(jīng)驗(yàn) +5 | 錘某人 2026-03-04 | 8/400 |
|
|
[考研] 化學(xué)工程求調(diào)劑 +10 | 化工人999 2026-03-04 | 10/500 |
|
|
[考研] 材料復(fù)試調(diào)劑 +7 | 學(xué)材料的點(diǎn) 2026-03-01 | 8/400 |
|
|
[考研] 一志愿985材料與化工 326分求調(diào)劑 +3 | Hz795795 2026-03-04 | 3/150 |
|
|
[考研] 一志愿中石油,化學(xué)工程,301求調(diào)劑 +3 | 豚饅枳枳 2026-03-04 | 3/150 |
|
|
[考研] 求調(diào)劑 +7 | 博斯特525 2026-03-04 | 7/350 |
|
|
[考研] 0703化學(xué) 一志愿東北大學(xué) 279分求調(diào)劑 +4 | 聶樂平 2026-03-03 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | G212 2026-03-03 | 5/250 |
|
|
[考研] 0854總分272 +5 | 打小就是老實(shí)人 2026-03-02 | 6/300 |
|
|
[考研] 環(huán)境工程學(xué)碩288求助調(diào)劑 +7 | 多吃億口芝士 2026-03-02 | 7/350 |
|
|
[考研] 江蘇省農(nóng)科院招調(diào)劑1名 +5 | Qwertyuop 2026-03-01 | 5/250 |
|
|
[考研] 278求調(diào)劑 +3 | 滿天星11_22 2026-03-02 | 3/150 |
|
|
[考研] 課題組接收材料類調(diào)劑研究生 +6 | gaoxiaoniuma 2026-02-28 | 9/450 |
|
|
[考研] 0854復(fù)試調(diào)劑 276 +5 | wmm9 2026-03-01 | 7/350 |
|
|
[考研] 284求調(diào)劑 +10 | 天下熯 2026-02-28 | 11/550 |
|
|
[考研] 一志愿鄭大材料學(xué)碩298分,求調(diào)劑 +6 | wsl111 2026-03-01 | 6/300 |
|
|
[考研] 調(diào)劑 +3 | 13853210211 2026-03-02 | 4/200 |
|
|
[考研] 291分工科求調(diào)劑 +9 | science餓餓 2026-03-01 | 10/500 |
|