| 查看: 1138 | 回复: 3 | ||||
| 【悬赏金币】回答本帖问题，作者@七安将赠送您 10 个金币 | ||||
[求助]
如何提取siesta输出的hamiltonian和overlap
|
||||
|
在输入文件中已设置Save HS True并且计算已输出label.HSX文件，现需要查看其中计算的哈密顿量，如何才能从该文件中提取哈密顿量？ 发自小木虫手机客户端 |
回复此楼» 猜你喜欢
欢迎211本科同学，过A区国家线，A区非偏远一本，交叉学科课题组
已经有26人回复
-
070300调剂一志愿南京大学化学专业，初试273分，有论文有专利有竞赛，求调剂B区211
已经有10人回复
-
新疆大学地质与矿业工程学院招生
已经有12人回复
-
材料学硕318求调剂
已经有12人回复
-
085602高分子方向求调剂
已经有7人回复
-
进入个人成果库好难，一下午都没进去
已经有3人回复
-
材料调剂
已经有11人回复
-
A区一本交叉课题组，低分调剂，招收机械电子信息通信等交叉方向
已经有37人回复
-
招调剂
已经有12人回复
-
一志愿山东大学105500药学专硕，总分302求调剂
已经有8人回复
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3503.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
|
虽然我没具体研究和折腾过，可行的方式有如下几种： 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
|
siesta编译已具有netcdf4模块，只是在siesta生成的TSHS或HSX文件中读取哈密顿量时有些困难，即使用Python语言提取的哈密顿量也是错误提取 发自小木虫手机客户端 |
3楼2025-08-18 17:55:52
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3503.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
4楼2025-08-20 10:15:33
5