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

相关推荐
egoist202341 分钟前
【Linux仓库】命令行参数与环境变量【进程·伍】
linux·运维·服务器·环境变量·命令行参数·内建命令
ζั͡山 ั͡有扶苏 ั͡✾43 分钟前
CentOS/RHEL LVM 磁盘扩展完整教程
linux·运维·centos
Fireworkitte1 小时前
Linux 中替换sed
linux·运维·服务器
gooxi_hui1 小时前
性能狂飙 Gooxi 8卡5090服务器重新定义高密度算力
运维·服务器
weixin_420571871 小时前
Windos服务器升级MySQL版本
运维·服务器·mysql
dessler1 小时前
Kafka-日常运维命令
linux·运维·kafka
誰能久伴不乏2 小时前
深入解析 TCP 连接状态与进程挂起、恢复与关闭
服务器·网络·tcp/ip
Dreamboat¿2 小时前
小皮面板搭建pikachu
运维·服务器·网络
weixin_443290692 小时前
【云服务器安全相关】堡垒机、WAF、防火墙、IDS 有什么区别?
服务器·安全·php
北北~Simple2 小时前
第一次搭建数据库
服务器·前端·javascript·数据库