| 5 | 1/1 | 返回列表 |
| 查看: 722 | 回復(fù): 2 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
龍山獅崖鐵桿木蟲 (正式寫手)
|
[求助]
晶體坐標變換問題
|
||
| 有一個晶體,晶胞三個方向的晶格長度分別為a,b,c,我想將晶體的前兩個方向變一下,也就是晶格長度變?yōu)閎,a,c,請問晶體內(nèi)的各個原子坐標該如何變化?用那些軟件能辦到?謝謝! |

木蟲 (著名寫手)
|
用線性代數(shù)里的坐標變換方法。 寫出新的基矢在舊基矢下面的轉(zhuǎn)換矩陣,然后用每個原子對應(yīng)的位置矢量乘以這個矩陣。這樣雖然麻煩點,但是卻是萬能的。這個腳本可以用matlab來寫。 設(shè)原來的基矢是v1,v2,v3,新基矢是u1,u2,u3,某個原子在舊基矢下的分量表示成[v|a],對應(yīng)的新基矢下的分量是[u|a],這個轉(zhuǎn)換關(guān)系寫成矩陣式就是[u|a] = [u|v]*[v|a]。[u|v]就是基組v在基組u下的表示矩陣。 例如v1 = 1 * u1 + 2 * u2 + 3 * u3,轉(zhuǎn)換矩陣的第一列就是transpose([1,2,3])。注意這里是列矢量而不是行矢量,轉(zhuǎn)換的時候原子坐標也是列矢量。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 江蘇科技大學(xué)招材料研究生 +4 | Su032713. 2026-04-01 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)338求調(diào)劑! +8 | Zuhui0306 2026-03-26 | 9/450 |
|
|
[考研] 275求調(diào)劑 +12 | waltzh 2026-04-01 | 12/600 |
|
|
[考研] 286求調(diào)劑 +5 | lim0922 2026-03-26 | 5/250 |
|
|
[考研]
|
廈大化工 2026-04-01 | 4/200 |
|
|
[考研] 314求調(diào)劑 +4 | 溪云珂 2026-03-26 | 4/200 |
|
|
[考研] 285求調(diào)劑 +5 | FZAC123 2026-03-30 | 5/250 |
|
|
[考研] 求調(diào)劑 +4 | DADA怪 2026-03-31 | 4/200 |
|
|
[考研] 352分-085602-一志愿985 +6 | 海納百川Ly 2026-03-29 | 6/300 |
|
|
[考研] 070300化學(xué)354求調(diào)劑 +15 | 101次希望 2026-03-28 | 15/750 |
|
|
[考研] 生醫(yī)工0831調(diào)劑求推薦 +4 | 小熊睿睿_s 2026-03-27 | 6/300 |
|
|
[考研]
|
小羊36 2026-03-30 | 4/200 |
|
|
[考研] 085601 329分調(diào)劑 +6 | yzsa12 2026-03-31 | 6/300 |
|
|
[考研] 085601一志愿西北工業(yè)大學(xué)初試346 +4 | 085601初試346 2026-03-30 | 4/200 |
|
|
[考研] 福建理工大學(xué)材料學(xué)院先進合金團隊招收考研調(diào)劑學(xué)生 +3 | 大華金商都 2026-03-30 | 4/200 |
|
|
[考研] 抱歉 +4 | 田洪有 2026-03-30 | 4/200 |
|
|
[考研] 11408總分309,一志愿東南大學(xué)求調(diào)劑,不挑專業(yè) +5 | 天賦帶到THU 2026-03-29 | 6/300 |
|
|
[考研] 一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑 +5 | 兔兔兔111223314 2026-03-29 | 7/350 |
|
|
[考研] 環(huán)境科學(xué)與工程334分求調(diào)劑 +6 | 王一一依依 2026-03-30 | 8/400 |
|
|
[考研] 296求調(diào)劑 +10 | 彼岸t 2026-03-29 | 10/500 |
|