| 3 | 1/1 | 返回列表 |
| 查看: 721 | 回復(fù): 2 | ||
龍山獅崖鐵桿木蟲 (正式寫手)
|
[求助]
晶體坐標(biāo)變換問(wèn)題
|
| 有一個(gè)晶體,晶胞三個(gè)方向的晶格長(zhǎng)度分別為a,b,c,我想將晶體的前兩個(gè)方向變一下,也就是晶格長(zhǎng)度變?yōu)閎,a,c,請(qǐng)問(wèn)晶體內(nèi)的各個(gè)原子坐標(biāo)該如何變化?用那些軟件能辦到?謝謝! |

鐵桿木蟲 (正式寫手)

木蟲 (著名寫手)
|
用線性代數(shù)里的坐標(biāo)變換方法。 寫出新的基矢在舊基矢下面的轉(zhuǎn)換矩陣,然后用每個(gè)原子對(duì)應(yīng)的位置矢量乘以這個(gè)矩陣。這樣雖然麻煩點(diǎn),但是卻是萬(wàn)能的。這個(gè)腳本可以用matlab來(lái)寫。 設(shè)原來(lái)的基矢是v1,v2,v3,新基矢是u1,u2,u3,某個(gè)原子在舊基矢下的分量表示成[v|a],對(duì)應(yīng)的新基矢下的分量是[u|a],這個(gè)轉(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)換的時(shí)候原子坐標(biāo)也是列矢量。 |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 083000學(xué)碩274求調(diào)劑 +5 | Li李魚 2026-03-26 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 322求調(diào)劑 +4 | 宋明欣 2026-03-27 | 4/200 |
|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +3 | 小張zxy張 2026-03-26 | 5/250 |
|
|
[考研] 321求調(diào)劑 +4 | 璞玉~~ 2026-03-25 | 4/200 |
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[考研] 349求調(diào)劑 +4 | 李木子啊哈哈 2026-03-25 | 4/200 |
|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +4 | 柳溪邊 2026-03-26 | 4/200 |
|
|
[考研] 324求調(diào)劑 +4 | wysyjs25 2026-03-21 | 4/200 |
|
|
[考研] 【雙一流院校新能源、環(huán)境材料,材料加工與模擬招收大量調(diào)劑】 +4 | Higraduate 2026-03-22 | 8/400 |
|
|
[考研] 291 求調(diào)劑 +7 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 8/400 |
|
|
[考研] 打過(guò)很多競(jìng)賽,085406控制工程300分,求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
[考研] 0854AI CV方向招收調(diào)劑 +4 | 章小魚567 2026-03-23 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
[考研]
|
黃粱一夢(mèng)千年 2026-03-24 | 3/150 |
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|
|
[考研] 一志愿國(guó)科過(guò)程所081700,274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 6/300 |
|