| 3 | 1/1 | 返回列表 |
| 查看: 1134 | 回復(fù): 2 | |||
smartcheng鐵蟲 (小有名氣)
|
[求助]
castep中限制原子只沿某個(gè)方向運(yùn)動(dòng)
|
|
我看前人的帖子說,只有修改.cell 文件,才能做到按方向性限制。如下所示,如果要限制第二個(gè)Al原子不能沿著X方向移動(dòng) 4 Al 2 1.0000000000 0.0000000000 0.0000000000 5 Al 2 0.0000000000 0.0000000000 1.0000000000 這樣怎么會(huì)達(dá)到限制x方向的目的,這里面的1和0都代表什么啊。 原帖: Q3:如何在CASTEP計(jì)算中限制某個(gè)原子的移動(dòng)方向? A3:CASTEP計(jì)算中只支持Fraction坐標(biāo),不支持Coordinate坐標(biāo)。采用常規(guī)方式無法限制原子的移動(dòng)方向,需要手動(dòng)修改CASTEP的Cell文件,具體方式如下,在 %BLOCK POSITIONS_FRAC Al 0.0000000000 0.0000000000 0.0000000000 Al 0.2270171283 0.2094837802 0.5009850469 Al 0.5000000000 0.0000000000 0.5000000000 Al 0.5000000000 0.5000000000 0.0000000000 %ENDBLOCK POSITIONS_FRAC 中包括每個(gè)原子的分?jǐn)?shù)坐標(biāo),請(qǐng)留意每個(gè)原子的編號(hào),在下面的計(jì)算中會(huì)有所幫助。 在Cell文件中,關(guān)于原子坐標(biāo)限制的部分如下: FIX_COM : false %BLOCK IONIC_CONSTRAINTS 1 Al 1 1.0000000000 0.0000000000 0.0000000000 2 Al 1 0.0000000000 1.0000000000 0.0000000000 3 Al 1 0.0000000000 0.0000000000 1.0000000000 4 Al 3 1.0000000000 0.0000000000 0.0000000000 5 Al 3 0.0000000000 1.0000000000 0.0000000000 6 Al 3 0.0000000000 0.0000000000 1.0000000000 7 Al 4 1.0000000000 0.0000000000 0.0000000000 8 Al 4 0.0000000000 1.0000000000 0.0000000000 9 Al 4 0.0000000000 0.0000000000 1.0000000000 %ENDBLOCK IONIC_CONSTRAINTS 其中包括每行編號(hào)(紫色),原子類型(紅色),原子編號(hào)(綠色)和需要限制的分?jǐn)?shù)坐標(biāo)(黑色)。如果要限制第二個(gè)Al原子不能沿著X方向移動(dòng),那么在該部分的設(shè)置如下: 1 Al 1 1.0000000000 0.0000000000 0.0000000000 2 Al 1 0.0000000000 1.0000000000 0.0000000000 3 Al 1 0.0000000000 0.0000000000 1.0000000000 4 Al 2 1.0000000000 0.0000000000 0.0000000000 5 Al 2 0.0000000000 0.0000000000 1.0000000000 6 Al 3 1.0000000000 0.0000000000 0.0000000000 7 Al 3 0.0000000000 1.0000000000 0.0000000000 8 Al 3 0.0000000000 0.0000000000 1.0000000000 9 Al 4 1.0000000000 0.0000000000 0.0000000000 10 Al 4 0.0000000000 1.0000000000 0.0000000000 11 Al 4 0.0000000000 0.0000000000 1.0000000000 注意,編號(hào)要順次后移,否則會(huì)出錯(cuò),黑色部分分別對(duì)應(yīng)于x、y、z三個(gè)方向,完成后,使用Run Files方式即可進(jìn)行計(jì)算。 此外,在 %BLOCK CELL_CONSTRAINTS 1 2 3 4 5 6 %ENDBLOCK CELL_CONSTRAINTS 中可以限制晶胞向量的變化,第一行分別表示三個(gè)晶胞參數(shù)a、b、c, 第二行表示晶胞角度α、β、γ。上述結(jié)構(gòu)表明參數(shù)可以被優(yōu)化,如果需要加以限制,則需要改為0。 |
計(jì)算MS等 |
金蟲 (正式寫手)

新蟲 (小有名氣)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑 +4 | 兔兔兔111223314 2026-03-29 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 各位老師好,我的一志愿為北京科技大學(xué)085601材料專碩 +9 | Koxui 2026-03-28 | 9/450 |
|
|
[考研] 307求調(diào)劑 +9 | 超級(jí)伊昂大王 2026-03-24 | 10/500 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5 | hanamiko 2026-03-29 | 5/250 |
|
|
[考研] 一志愿雙一流機(jī)械285分求調(diào)劑 +4 | 幸運(yùn)的三木 2026-03-29 | 5/250 |
|
|
[考研] 297求調(diào)劑 +11 | 田洪有 2026-03-26 | 11/550 |
|
|
[考研] 375求調(diào)劑 +4 | 雨夏整夜 2026-03-29 | 4/200 |
|
|
[考研] 343求調(diào)劑 +6 | 愛羈絆 2026-03-29 | 6/300 |
|
|
[考研] 學(xué)碩274求調(diào)劑 +9 | Li李魚 2026-03-26 | 9/450 |
|
|
[考研] 346求調(diào)劑 一志愿070303有機(jī)化學(xué) +3 | 蘿卜燉青菜 2026-03-28 | 3/150 |
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
[考研] 286求調(diào)劑 +4 | 丟掉懶惰 2026-03-27 | 7/350 |
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +9 | 丹青奶蓋 2026-03-26 | 10/500 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
[考研] 0703化學(xué)338求調(diào)劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-23 | 4/200 |
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|