24小时热门版块排行榜    

查看: 1138  |  回复: 3
【悬赏金币】回答本帖问题,作者@七安将赠送您 10 个金币

@七安

新虫 (小有名气)

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

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

卡开发发

专家顾问 (著名写手)

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的回帖

卡开发发

专家顾问 (著名写手)

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的回帖
相关版块跳转 我要订阅楼主 @七安 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085602高分子方向求调剂 +7 tlgudy 2026-03-04 7/350 2026-03-06 14:17 by guosr9609
[考研] 一志愿哈尔滨工业大学0856材料与化工,前三科206,总分283,求调剂 +4 26考研求调剂 2026-03-06 4/200 2026-03-06 12:05 by leocomeon
[考研] 材料工程调剂 +5 zbcm_zbcm 2026-03-05 5/250 2026-03-06 10:16 by 好好好1233
[考研] 材料277分求调剂 +12 饭饭星球 2026-03-04 12/600 2026-03-06 08:13 by 杨杨杨紫
[考研] 271求调剂 +7 月色c 2026-03-05 8/400 2026-03-05 19:51 by wangjihu
[考研] 085602 293分求调剂 +3 SivanNano. 2026-03-05 3/150 2026-03-05 19:41 by caszguilin
[考研] 347求调剂 +6 啊欧欧欧 2026-03-03 8/400 2026-03-05 11:40 by 0202liuyan
[考研] 化学 0703求调剂 总分293 一志愿211 +9 土土小虫 2026-03-03 11/550 2026-03-05 09:26 by vfdbghn
[考研] 0703化学 学硕 理工科均可 不区分研究方向 总分279求调剂 +7 1一11 2026-03-03 7/350 2026-03-05 09:23 by xwxstudy
[考研] 22408-273求调剂-担任3个项目的负责人-1篇国际期刊论文(一作)1篇核心期刊论文在投。 +3 没想好取什么名 2026-03-03 3/150 2026-03-05 08:38 by Destiny688
[考研] 求调剂 +7 博斯特525 2026-03-04 7/350 2026-03-04 18:09 by houyaoxu
[考研] 085602化学工程350,调剂,有没有211的 +5 利好利好. 2026-03-02 9/450 2026-03-03 17:06 by 利好利好.
[考研] 材料工程专硕283求调剂 5+8 ,!? 2026-03-02 10/500 2026-03-03 12:01 by EBSD
[考研] 材料工程求调剂 +3 1431251 2026-03-03 3/150 2026-03-03 11:58 by EBSD
[考研] 0854复试调剂 276 +5 wmm9 2026-03-01 7/350 2026-03-03 02:49 by xiadaiyang
[考研] 306分材料调剂 +5 chuanzhu川烛 2026-03-01 6/300 2026-03-02 19:51 by 张晓芳0105
[考研] 275求调剂 +7 明远求学 2026-03-01 7/350 2026-03-02 19:22 by zhukairuo
[考博] 诚招农业博士 +3 心欣向荣 2026-02-28 3/150 2026-03-02 13:33 by 时间不狗
[考研] 材料调剂 +6 爱擦汗的可乐冰 2026-02-28 7/350 2026-03-02 10:42 by Jy?
[考研] 272求调剂 +6 田智友 2026-02-28 6/300 2026-03-01 21:40 by 公瑾逍遥
信息提示
请填处理意见