国产主播欧美精品,在线视频cao,中文字幕免费一区二区,日韩三级影视

Document

當前位置：文庫百科 ? 文章詳情

基于Quantum Espresso的第一性原理計算建模教程（三）：晶胞和原胞的相互轉(zhuǎn)換

來源：本站時間：2020-11-11 14:52:26 瀏覽：12764次

原胞是保持平移對稱性的最小單元，所以，在計算能帶、聲子色散時，研究對象是原胞（聲子的有限位移方法需要超胞，但是，這時的超胞是一個輔助系統(tǒng)，色散仍然是針對原胞畫的）。用比原胞更大的單元計算能帶會造成band folding現(xiàn)象，即能帶折疊，相當于布里淵區(qū)變小了，布里淵區(qū)的形狀也可能發(fā)生了變化，從超胞到原胞的能帶可以通過(band unfolding)[https://www.ifm.liu.se/theomod/compphys/band-unfolding]還原到原胞的能帶。

文獻通常按照以下約定：晶胞是定義晶面、晶向、超胞的參照，而不是原胞或其他單元。

首先，定義一般的周期性單元CELL的變換。按照文獻[1]的約定，將（分數(shù)）坐標寫為列矢量，基矢量

也各為列矢量。點X在基矢

（O為原點）下的坐標

考慮晶格靜止不動，選擇不同的基矢，即選取不同的CELL，同一個點X對新的基矢有坐標

。下面給出有撇號和無撇號的基矢選擇下基矢和坐標的變換關(guān)系。

周期性單元CELL的變換是一種特殊的仿射變換（數(shù)學(xué)上一般的仿射變換形成的新單元并不保證具有晶格周期性），可以分解為線性部分和平移兩個部分。

線性部分包括基矢方向和長度的改變，由一個矩陣P表示。

基于Quantum Espresso的第一性原理計算建模教程（三）：晶胞和原胞的相互轉(zhuǎn)換-5

基于Quantum Espresso的第一性原理計算建模教程（三）：晶胞和原胞的相互轉(zhuǎn)換-5

知道了矩陣P和逆矩陣

，就可以進行CELL之間的變換，常見的晶胞轉(zhuǎn)換為原胞的變換矩陣由表2給出，反之交換P和Q得到。注意這里的矩陣選取并不是唯一的，這里選擇的原胞基矢的原點和晶胞基矢的原點是重合的，即沒有平移。

表2 常見單元變換矩陣

基于Quantum Espresso的第一性原理計算建模教程（三）：晶胞和原胞的相互轉(zhuǎn)換-7

基于Quantum Espresso的第一性原理計算建模教程（三）：晶胞和原胞的相互轉(zhuǎn)換-7

對于菱方布拉伐格子，見注2的7種空間群，cif生成的是六方的晶胞，體積是菱方的3倍，原胞計算用（ibrav=5，同時定義celldm(1)和celldm(4)），用ibrav=0時也存著晶胞轉(zhuǎn)原胞的問題，轉(zhuǎn)換矩陣如下，參見注1 ：

基于Quantum Espresso的第一性原理計算建模教程（三）：晶胞和原胞的相互轉(zhuǎn)換-8

基于Quantum Espresso的第一性原理計算建模教程（三）：晶胞和原胞的相互轉(zhuǎn)換-8

注1：注意晶胞和原胞的區(qū)別，對于非簡單格子ibrav適合于設(shè)置原胞（對于簡單格子ibrav≠0當然也是設(shè)置了原胞），布拉伐格子中的7個簡單格子本身就是原胞，而且，除了菱方外的6個簡單格子，不僅是原胞，同時也是晶胞，菱方的布拉伐格子是原胞但不是晶胞，菱方的晶胞是六方的簡單格子，體積是原胞的3倍，而底心、面心、體心的7個布拉伐格子本身是晶胞，存在體積更小的原胞。

注2：trigonal三方晶系有兩種布拉伐格子，一種是ibrav=5，菱方（rhombohedral）布拉伐格子，另一種是ibrav=4，六方（hexgonal）布拉伐格子，晶體屬于菱方還是六方要看具體的空間群，在hexgonal和trigonal晶系中，7個空間群：具有菱方布拉伐格子的原胞，其余的45個空間群具有六方布拉伐格子的原胞。這里的菱方和六方是指晶體的格點系統(tǒng)lattice system，而非晶系，格點系統(tǒng)是按照布拉伐格子分類的，晶系是按照晶體點群分類的。

周期性單元CELL的變換，還可以包含平移，平移用變換前的基矢定義為：

。平移的逆變換

分數(shù)坐標的變換公式為：基于Quantum Espresso的第一性原理計算建模教程（三）：晶胞和原胞的相互轉(zhuǎn)換-11

將晶胞轉(zhuǎn)換為原胞，因為晶胞和原胞中的原子個數(shù)不同，坐標變換后會有重復(fù)或相差一個原胞格子，需要將重復(fù)的刪去。反之，將原胞轉(zhuǎn)換為晶胞，需要將原胞格子重復(fù)足夠大以覆蓋晶胞，重復(fù)后的原子分別變換到新的分數(shù)坐標，最后保留晶胞內(nèi)部的原子（分數(shù)坐標在0到1之間）。

以底心單斜碳的同素異形體為例[4]，用TransCell（https://github.com/yyyu200/SlabMaker）變換CELL和原子坐標，并用VESTA驗證。

先下載cif文件（https://yyyu200.github.io/DFTbook/img/A_mC16_12_4i.cif），用VESTA打開，畫出晶胞如下：