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

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

| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 295求調(diào)劑 +4 | 小賽不吃香菜 2026-03-04 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 本科太原理工采礦工程,求調(diào)劑 +3 | onlx 2026-03-01 | 3/150 |
|
|
[考研]
085600 英一數(shù)二272求調(diào)劑
5+6
|
vida_a 2026-03-01 | 47/2350 |
|
|
[考博] 26申博 +5 | north, 2026-02-28 | 5/250 |
|
|
[考博] 26申博-目前4篇SCI一作 +3 | chen_2024 2026-03-02 | 3/150 |
|
|
[考研] 289求調(diào)劑 +8 | yang婷 2026-03-02 | 10/500 |
|
|
[考研] 材料工程269求調(diào)劑 +7 | 白刺玫 2026-03-02 | 7/350 |
|
|
[考研] 化工專碩調(diào)劑 +4 | 利好利好. 2026-03-03 | 7/350 |
|
|
[考研] 284求調(diào)劑 +6 | 天下熯 2026-03-02 | 6/300 |
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
[考研] 考研復(fù)試調(diào)劑,過(guò)國(guó)家線的同學(xué)都可報(bào)名 +7 | 黑!在干嘛 2026-02-28 | 8/400 |
|
|
[考研] 化工京區(qū)271求調(diào)劑 +7 | 11ing 2026-03-02 | 7/350 |
|
|
[考研] 290分材料工程085601求調(diào)劑 數(shù)二英一 +8 | llx0610 2026-03-02 | 9/450 |
|
|
[考研] 求調(diào)劑 +7 | repeatt?t 2026-02-28 | 7/350 |
|
|
[考研] 材料工程274求調(diào)劑 +5 | Lilithan 2026-03-01 | 5/250 |
|
|
[考研] 一志愿東北大學(xué)材料專碩328,求調(diào)劑 +3 | shs1083 2026-03-02 | 3/150 |
|
|
[考研] 材料與化工328求調(diào)劑 +3 | 。,。,。,。i 2026-03-02 | 3/150 |
|
|
[考研] 哈工大計(jì)算機(jī)劉劼團(tuán)隊(duì)招生 +4 | hit_aiot 2026-03-01 | 6/300 |
|
|
[考研] 295求調(diào)劑 +8 | 19171856320 2026-02-28 | 8/400 |
|
|
[考研] 272求調(diào)劑 +6 | 田智友 2026-02-28 | 6/300 |
|