24小时热门版块排行榜    

查看: 1113  |  回复: 3
【悬赏金币】回答本帖问题,作者@七安将赠送您 10 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

@七安

新虫 (小有名气)

[求助] 如何提取siesta输出的hamiltonian和overlap

在输入文件中已设置Save HS  True并且计算已输出label.HSX文件,现需要查看其中计算的哈密顿量,如何才能从该文件中提取哈密顿量?
发自小木虫手机客户端
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
3楼: Originally posted by @七安 at 2025-08-18 17:55:52
siesta编译已具有netcdf4模块,只是在siesta生成的TSHS或HSX文件中读取哈密顿量时有些困难,即使用Python语言提取的哈密顿量也是错误提取
...

那你看看方案1或2呢?不过2的方法可能比较土也比较折腾。
不一定挂在论坛,计算问题问题欢迎留言。
4楼2025-08-20 10:15:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

虽然我没具体研究和折腾过,可行的方式有如下几种:
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的库并写程序读取之。
不一定挂在论坛,计算问题问题欢迎留言。
2楼2025-08-04 06:26:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

@七安

新虫 (小有名气)

引用回帖:
2楼: Originally posted by 卡开发发 at 2025-08-04 06:26:17
虽然我没具体研究和折腾过,可行的方式有如下几种:
1、直接读取HSX。但HSX是Fortran写入的二进制文件,可能需要去了解HSX读取或者写入逻辑。例如SIESTA程序提供的Utils当中有hsx2hs.f90,其中有读取HSX文件的逻辑 ...

siesta编译已具有netcdf4模块,只是在siesta生成的TSHS或HSX文件中读取哈密顿量时有些困难,即使用Python语言提取的哈密顿量也是错误提取
发自小木虫手机客户端
3楼2025-08-18 17:55:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0857调剂 +5 一ll半 2026-02-28 6/300 2026-03-03 00:52 by 柒槿levana
[考研] 求调剂 +11 yunziaaaaa 2026-03-01 13/650 2026-03-02 21:59 by sunny81
[考研] 321求调剂一志愿东北林业大学材料与化工英二数二 +5 虫虫虫虫虫7 2026-03-01 9/450 2026-03-02 21:33 by sunny81
[考研] 085600材料工程一志愿中科大总分312求调剂 +9 吃宵夜1 2026-02-28 11/550 2026-03-02 20:14 by hypershenger
[考研] 302材料工程求调剂 +5 Doleres 2026-03-01 6/300 2026-03-02 19:53 by 张晓芳0105
[考研] 306分材料调剂 +5 chuanzhu川烛 2026-03-01 6/300 2026-03-02 19:51 by 张晓芳0105
[考研] 261求调剂 +3 陆lh 2026-03-01 3/150 2026-03-02 19:32 by zhukairuo
[考研] 275求调剂 +7 明远求学 2026-03-01 7/350 2026-03-02 19:22 by zhukairuo
[考博] 博士自荐 +4 kkluvs 2026-02-28 5/250 2026-03-02 19:19 by 轻松不少随
[考研] 一志愿东北大学化学314分求调剂 +3 lr1212.. 2026-03-02 3/150 2026-03-02 17:36 by yeahyou
[考研] 一志愿东北大学材料专硕328,求调剂 +3 shs1083 2026-03-02 3/150 2026-03-02 17:27 by houyaoxu
[考研] 0856化工专硕求调剂 +15 董boxing 2026-03-01 15/750 2026-03-02 15:06 by 晃晃不许晃
[考研] 材料类求调剂 +11 wana_kiko 2026-02-28 14/700 2026-03-02 08:46 by 聪明的大松鼠
[考研] 322求调剂 +3 熊境喆 2026-03-01 3/150 2026-03-02 08:44 by houyaoxu
[考博] 26申博 +4 想申博! 2026-02-26 6/300 2026-03-01 17:32 by 想申博!
[考研] 313求调剂 +3 水流年lc 2026-02-28 3/150 2026-03-01 16:01 by 新能源达人
[考研] 311求调剂 +6 亭亭亭01 2026-03-01 6/300 2026-03-01 15:41 by 324616
[考研] 寻找调剂 +4 LYidhsjabdj 2026-02-28 4/200 2026-03-01 10:56 by sunny81
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
请填处理意见