EON安装ASE Interface

EON安装

我的eon路径于/eon/。

则环境为

复制代码
export PYTHONPATH=/eon/:$PYTHONPATH
export PATH=/eon/bin:$PATH

source: https://theory.cm.utexas.edu/eon/installation.html

ASE

测试系统ubuntu。如果你python2和python3总是纠缠不清,可以sudo apt install python-is-python3直接解决。

经检查,我PC的

python地址为:

/usr/include/python3.8/

pybind11地址为:

/usr/include/pybind11/

已确认python地址下存有python.h,pybind11地址下存有pybind.h。

即得到

复制代码
CXXFLAGS += -DASE_POT -I/usr/include/pybind11/ -I/usr/include/python3.8/

第二部分,libpython3.8.so地址在:

/usr/lib/x86_64-linux-gnu/libpython3.8.so

得到

复制代码
LDFLAGS += -L/usr/lib/x86_64-linux-gnu/ -lpython3.8

如果什么都找不到可以通过以下代码搜索:

find / -name libpython3.8.so

vim 位于$EON_DIR/client/的Rules.mk,修改即可。

此时在client目录下make ASE_POT=1即可查看结果。

使用

教程表示要eon client可以连到Python dynamic library,就是LDFLAGS上的地址

复制代码
export LD_LIBRARY_PATH=/path/to/python/rootdir/lib:$LD_LIBRARY_PATH

在程序的config.ini里面添加如下内容即可让EON利用ASE。

复制代码
[Potential]
potential = ase
ext_pot_path = /full/path/to/the/ase_script.py

不过我觉得我系统已经有类似export LD_LIBRARY_PATH=/path/to/python/rootdir/lib:$LD_LIBRARY_PATH这样的操作了,便跳过了这一步。

这里最好用absolute path而不是relative path。
然后这个python code我则放到了eon的位置。

复制代码
from ase import Atoms
from ase.calculators.lj import LennardJones

def ase_calc():
    # MODIFY THIS SECTION
    calc = LennardJones(epsilon=0.0103, sigma=3.40, rc=10.0, ro=0.0, smooth=True)
    return calc

#=======================================================================
# DO NOT EDIT
def _calculate(R, atomicNrs, box, calc):
    system = Atoms(symbols=atomicNrs, positions=R, pbc=True, cell=box)
    system.calc = calc
    forces = system.get_forces()
    energy = system.get_potential_energy()
    return energy, forces
#=======================================================================


if __name__ == '__main__':
    # This is just to verify that the script runs properly with `python3 <this_script>`
    from ase.io import read
    atoms = read("...")  # structure file name
    pos = atoms.positions  # array of positions (n x 3)
    atomicNrs = atoms.get_atomic_numbers()  # array of atomic numbers (n)
    cell = atoms.cell  # cell array (3 x 3)
    calc = ase_calc()
    e, f = _calculate(pos, atomicNrs, cell, calc)
    print(f"E = {e}")
    print("F =")
    print(f)

测试

Source

https://theory.cm.utexas.edu/eon/ase_interface.html#ase-interface

相关推荐
YC运维1 小时前
网络配置综合实验全攻略(对之前学习的总结)
linux·服务器·网络
MarkGosling2 小时前
【开源项目】网络诊断告别命令行!NetSonar:开源多协议网络诊断利器
运维·后端·自动化运维
无敌的牛2 小时前
Linux基础开发工具
linux·运维·服务器
Edingbrugh.南空2 小时前
实战指南:用pmap+gdb排查Linux进程内存问题
linux·运维·服务器
朱颜辞镜花辞树‎3 小时前
GitHub 操作指南:项目协作与自动化工作流实践
运维·自动化·github
孙克旭_3 小时前
day051-ansible循环、判断与jinja2模板
linux·运维·服务器·网络·ansible
总有刁民想爱朕ha4 小时前
零基础搭建监控系统:Grafana+InfluxDB 保姆级教程,5分钟可视化服务器性能!
运维·服务器·grafana
Mr_Orangechen4 小时前
Linux 下使用 VS Code 远程 GDB 调试 ARM 程序
linux·运维·arm开发
撰卢5 小时前
【个人笔记】负载均衡
运维·笔记·负载均衡
Shartin5 小时前
Can201-Introduction to Networking: Application Layer应用层
服务器·开发语言·php