| 24小時(shí)熱門版塊排行榜 |
| 4 | 1/1 | 返回列表 |
| 查看: 1159 | 回復(fù): 3 | ||||
| 【懸賞金幣】回答本帖問題,作者@七安將贈(zèng)送您 10 個(gè)金幣 | ||||
@七安新蟲 (小有名氣)
|
[求助]
如何提取siesta輸出的hamiltonian和overlap
|
|||
|
在輸入文件中已設(shè)置Save HS True并且計(jì)算已輸出label.HSX文件,現(xiàn)需要查看其中計(jì)算的哈密頓量,如何才能從該文件中提取哈密頓量? 發(fā)自小木蟲手機(jī)客戶端 |
專家顧問 (著名寫手)
Ab Initio Amateur
![]() |
專家經(jīng)驗(yàn): +224 |
|
雖然我沒具體研究和折騰過,可行的方式有如下幾種: 1、直接讀取HSX。但HSX是Fortran寫入的二進(jìn)制文件,可能需要去了解HSX讀取或者寫入邏輯。例如SIESTA程序提供的Utils當(dāng)中有hsx2hs.f90,其中有讀取HSX文件的邏輯,或者你可以參考DeepH程序,當(dāng)中preprocess/siesta_get_data.py提供了python語言處理HSX的功能。 2、對(duì)SIESTA源碼部分進(jìn)行修改,將HSX的寫入格式從unformatted改成formatted后重新編譯,這樣你能得到一份“人類可讀”的HSX,然后寫程序讀取之,只不過可能其他模塊并不兼容。 3、在SIESTA編譯時(shí)掛上NetCDF4,然后計(jì)算時(shí)使用Write.DMHS.History.NetCDF,然后HS可以按照NetCDF4進(jìn)行存儲(chǔ),使用相應(yīng)語言的NetCDF4的庫(kù)并寫程序讀取之。 |

新蟲 (小有名氣)
|
siesta編譯已具有netcdf4模塊,只是在siesta生成的TSHS或HSX文件中讀取哈密頓量時(shí)有些困難,即使用Python語言提取的哈密頓量也是錯(cuò)誤提取 發(fā)自小木蟲手機(jī)客戶端 |
專家顧問 (著名寫手)
Ab Initio Amateur
![]() |
專家經(jīng)驗(yàn): +224 |

| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 復(fù)試調(diào)劑 +6 | 呼呼?~+123456 2026-03-08 | 7/350 |
|
|---|---|---|---|---|
|
[碩博家園] 木蟲好像不熱鬧了,是不是? +3 | 偏振片 2026-03-10 | 3/150 |
|
|
[考研] 一志愿南大化學(xué)339分求調(diào)劑,四六級(jí)已過,有比賽,有文章 +7 | Gallantzhou 2026-03-07 | 7/350 |
|
|
[考研] 0703求調(diào)劑 +10 | c297914 2026-03-06 | 11/550 |
|
|
[考研] 0856求調(diào)劑 +3 | squirtle11 2026-03-07 | 3/150 |
|
|
[考研] 0703化學(xué)求調(diào)劑,總分309分,一志愿華南師范 +3 | 花與葉@ 2026-03-08 | 3/150 |
|
|
[考博] 26申博自薦 +9 | 烷基八氮- 2026-03-06 | 11/550 |
|
|
[考研] 求調(diào)劑,一志愿江南大學(xué),食品科學(xué)與工程,總分,320 +3 | yyyyyukino 2026-03-07 | 3/150 |
|
|
[考研] 理學(xué)287分求調(diào)劑 +4 | 看看我. 2026-03-05 | 4/200 |
|
|
[考研] 0701-322 求調(diào)劑 +3 | jiliuxian 2026-03-06 | 8/400 |
|
|
[教師之家] 交大前校長(zhǎng)王樹國(guó):現(xiàn)在最先進(jìn)的科技并不在大學(xué)實(shí)驗(yàn)室,而是在企業(yè)研究院 +4 | zju2000 2026-03-08 | 6/300 |
|
|
[考研] 0817化學(xué)工程與技術(shù)312分求調(diào)劑 +7 | T123 tt 2026-03-04 | 7/350 |
|
|
[考研] 材料專碩290求調(diào)劑 +8 | 杰尼龜aaa 2026-03-04 | 8/400 |
|
|
[考研] 08工科求調(diào)劑 +3 | 隆LLL 2026-03-06 | 4/200 |
|
|
[考研] 304求調(diào)劑 +4 | 曼殊2266 2026-03-05 | 4/200 |
|
|
[考研] 085701環(huán)境工程 求調(diào)劑 +7 | xiiiia 2026-03-04 | 7/350 |
|
|
[考研] 320材料一志愿西工大專碩英二數(shù)二 有兩年光伏工作經(jīng)驗(yàn) +5 | 錘某人 2026-03-04 | 8/400 |
|
|
[考研] 武理材料工程302調(diào)劑 +3 | Doleres 2026-03-03 | 6/300 |
|
|
[考研] 312求調(diào)劑 +8 | 醋精華了一下發(fā)?/a> 2026-03-03 | 10/500 |
|
|
[考研] 一志愿西交大材料學(xué)碩 346 求調(diào)劑 +3 | zju51 2026-03-04 | 3/150 |
|