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

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

| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 264求調劑 +7 | 26調劑 2026-03-03 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 085600 材料與化工 298 +8 | 小西笑嘻嘻 2026-03-03 | 8/400 |
|
|
[考研] 歡迎采礦、地質、巖土、計算機、人工智能等專業(yè)的同學報考 +8 | pin8023 2026-02-28 | 11/550 |
|
|
[考研] 環(huán)境工程學碩288求助調劑 +7 | 多吃億口芝士 2026-03-02 | 7/350 |
|
|
[基金申請] 請問大家,研究風險與應對措施那里, 大家都怎么寫呢 ? +3 | cauasen 2026-03-02 | 3/150 |
|
|
[考研] 085700資環(huán)求調劑,初始279,六級已過,英語能力強 +3 | 085700資環(huán)調劑 2026-03-03 | 4/200 |
|
|
[考研] 085602化學工程350,調劑,有沒有211的 +5 | 利好利好. 2026-03-02 | 9/450 |
|
|
[考研] 一志愿天津大學085600 319分 材料與化工 金屬方向 求調劑 +4 | 青科11 2026-03-02 | 4/200 |
|
|
[考研] 材料工程求調劑 +3 | 1431251 2026-03-03 | 3/150 |
|
|
[考博] 26申博 +4 | north, 2026-02-28 | 4/200 |
|
|
[考研] 課題組接收材料類調劑研究生 +6 | gaoxiaoniuma 2026-02-28 | 9/450 |
|
|
[考研]
|
好好好1233 2026-02-28 | 16/800 |
|
|
[考研] 0856材料求調劑 +12 | hyf hyf hyf 2026-02-28 | 13/650 |
|
|
[考研] 0856求調劑285 +11 | 呂仔龍 2026-02-28 | 11/550 |
|
|
[考研] 302材料工程求調劑 +5 | Doleres 2026-03-01 | 6/300 |
|
|
[考研] 306分材料調劑 +5 | chuanzhu川燭 2026-03-01 | 6/300 |
|
|
[考研] 0856化工專碩求調劑 +15 | 董boxing 2026-03-01 | 15/750 |
|
|
[考研] 材料學調劑 +10 | 提神豆沙包 2026-02-28 | 12/600 |
|
|
[考研] 322求調劑 +3 | 熊境喆 2026-03-01 | 3/150 |
|
|
[考博] 26申博 +4 | 想申博! 2026-02-26 | 6/300 |
|