【LAMMPS学习】八、基础知识(6.5)PyLammps 教程

8. 基础知识

此部分描述了如何使用 LAMMPS 为用户和开发人员执行各种任务。术语表页面还列出了 MD 术语,以及相应 LAMMPS 手册页的链接。 LAMMPS 源代码分发的 examples 目录中包含的示例输入脚本以及示例脚本页面上突出显示的示例输入脚本还展示了如何设置和运行各种模拟。

8.1.通用基础知识

8.2. 设置入门

8.3. 分析入门

8.4. 力场入门

8.5. 软件包入门

8.6.操作教程

8.6.1.将 CMake 与 LAMMPS 结合使用

8.6.2. LAMMPS GitHub 教程

8.6.3.使用 LAMMPS GUI

8.6.4.Moltemplate教程

8.6.5. PyLammps 教程

  • 概述
  • PyLammps 是 LAMMPS 的 Python 包装类,可以单独创建或使用现有的lammps Python 对象。与使用 Python ctypes 编写的 C 风格 LAMMPS 库接口的 lammps 包装器相比,它为常见的 LAMMPS 功能创建了一个更简单、更“Pythonic”的接口。使用 Python 与 LAMMPS 文档页面讨论了 lammps 包装器。

    与扁平的 ctypes 界面不同,PyLammps 公开了一个可发现的 API。它不再需要了解底层 C++ 代码实现的知识。最后, IPyLammps 包装器构建在 PyLammps 之上,并添加了一些用于将 IPython 集成到 Jupyter 笔记本中的附加功能,例如用于转储样式图像的嵌入式可视化输出。

  • lampmps 和 PyLammps 接口的比较
  • lammps.lammps 
  • 使用 ctypes

  • 直接内存访问本机 C++ 数据

  • 提供向 LAMMPS 发送和接收数据的函数

  • 需要了解 LAMMPS 内部如何工作(C 指针等)

  • lammps.PyLammps
  • 建立在原始 ctypes 接口之上的更高级别的抽象

  • Python 对象的操作

  • API 用户无法看到与 LAMMPS 的通信

  • 更短、更简洁的Python

  • 更好的 IPython 集成,专为快速原型设计而设计

  • 快速开始
  • 整个系统的安装
 第 1 步:将 LAMMPS 构建为共享库

要在 Python 中使用 LAMMPS,必须将其编译为共享库。然后,Python 接口会加载该库。在此示例中,我们启用 MOLECULE 包并编译 LAMMPS,并启用 PNG、JPEG 和 FFMPEG 输出支持。

步骤 1a:对于基于 CMake 的构建系统,步骤如下:

mkdir $LAMMPS_DIR/build-shared
cd  $LAMMPS_DIR/build-shared

# MPI, PNG, Jpeg, FFMPEG are auto-detected
cmake ../cmake -DPKG_MOLECULE=yes -DBUILD_LIB=yes -DBUILD_SHARED_LIBS=yes
make

步骤 1b:对于遗留的、基于 make 的构建系统,步骤如下:

cd $LAMMPS_DIR/src

# add packages if necessary
make yes-MOLECULE

# compile shared library using Makefile
make mpi mode=shlib LMP_INC="-DLAMMPS_PNG -DLAMMPS_JPEG -DLAMMPS_FFMPEG" JPG_LIB="-lpng -ljpeg"
步骤2:安装LAMMPS Python包

PyLammps 是 lampmps Python 包的一部分。要安装它,只需将该包安装到当前的 Python 安装中:

make install-python

注:重新编译共享库需要重新安装Python包

在 virtualenv 内安装

您可以使用 virtualenv 创建专门针对您的工作流程进行调整的自定义 Python 环境。

使用 virtualenv 的好处
  • 将系统 Python 安装与开发安装隔离

  • 安装可以在您的用户目录中进行,无需 root 访问权限(对于 HPC 集群很有用)

  • 通过 pip 安装软件包可以让您获得比通过 apt-get 或 yum 软件包管理器(并且没有 root 访问权限)更新版本的软件包

  • 如果需要,您甚至可以安装特定旧版本的软件包

先决条件(例如在 Ubuntu 上)

apt-get install python-virtualenv
创建安装了lammps的virtualenv
# create virtualenv named 'testing'
virtualenv $HOME/python/testing

# activate 'testing' environment
source $HOME/python/testing/bin/activate

现在,如上所述配置并编译 LAMMPS 共享库。当使用 CMake 并且共享库已经构建时,您需要重新运行 CMake 将 python 可执行文件的位置更新为虚拟环境中的位置:

cmake . -DPython_EXECUTABLE=$(which python)

# install LAMMPS package in virtualenv
(testing) make install-python

# install other useful packages
(testing) pip install matplotlib jupyter mpi4py

...

# return to original shell
(testing) deactivate
创建 PyLammps 的新实例

要创建 PyLammps 对象,您需要首先从 lammps 模块导入该类。通过使用默认构造函数,创建一个新的 lampmps 实例。

from lammps import PyLammps
L = PyLammps()

命令您还可以在此现有的 lampmps 对象之上初始化 PyLammps:​​​​​​​

from lammps import lammps, PyLammps
lmp = lammps()
L = PyLammps(ptr=lmp)
命令

使用lammps对象实例的command方法可以通过现有的库接口发送LAMMPS命令。

例如,我们采用以下 LAMMPS 命令:

region box block 0 10 0 5 -0.5 0.5

在原始界面中,如果 L 是 lammps 实例,则可以使用以下 Python 代码执行此命令:

L.command("region box block 0 10 0 5 -0.5 0.5")

使用 PyLammps 接口,任何命令都可以分为由空格分隔的任意部分,作为单独的参数传递给区域方法。

L.region("box block", 0, 10, 0, 5, -0.5, 0.5)

请注意,每个参数都设置为 Python 文字浮点数。在 PyLammps 界面中,每个命令都采用任意参数列表,并将其透明地合并为单个命令字符串,并用空格分隔各个参数。

这种方法的好处是避免冗余命令调用并且更容易参数化。在原始接口中,参数化需要通过创建格式化字符串来手动完成。

L.command("region box block %f %f %f %f %f %f" % (xlo, xhi, ylo, yhi, zlo, zhi))
相比之下,PyLammps 的方法直接接受参数,并将它们自动转换为最终的命令字符串。
L.region("box block", xlo, xhi, ylo, yhi, zlo, zhi)
系统状态

除了直接通过 PyLammps 对象调度命令之外,它还提供了几个允许您查询系统状态的属性。

L.system

是描述系统的字典,例如边界框或原子数

L.system.xlo, L.system.xhi

沿 x 轴的边界框限制

L.system.ylo、L.system.yhi

沿 y 轴的边界框限制

L.system.zlo、L.system.zhi

沿 z 轴的边界框限制

L.communication

通信子系统的配置,例如线程或处理器的数量

L.communication.nthreads

每个 LAMMPS 进程使用的线程数

L.communication.nprocs

LAMMPS 使用的 MPI 进程数

L.fixes 
当前系统中的fixes列表

L.computes

当前系统中活动计算的列表

L.dump

当前系统中活动转储的列表

L.groups

当前系统中存在的组列表

使用 LAMMPS 变量

LAMMPS 变量可以通过 PyLammps 接口定义和访问。

要定义变量,您可以使用变量命令:

L.variable("a index 2")

L.variables 返回所有变量的字典

您可以通过按名称从 L.variables 字典中检索变量对象来访问单个变量

a = L.variables['a']

然后可以通过访问该对象的 value 属性轻松读取和写入变量值。

print(a.value)
a.value = 4
检索任意 LAMMPS 表达式的值

LAMMPS 表达式可以使用 eval 方法立即求值。传递的字符串参数可以是包含全局热值、变量、计算或固定数据的任何表达式。

result = L.eval("ke") # kinetic energy
result = L.eval("pe") # potential energy

result = L.eval("v_t/2.0")
访问原子数据

可以使用 L.atoms 列表访问当前模拟中的所有原子。该列表的每个元素都是一个公开其属性(id、类型、位置、速度、力等)的对象。

# access first atom
L.atoms[0].id
L.atoms[0].type

# access second atom
L.atoms[1].position
L.atoms[1].velocity
L.atoms[1].force

还可以使用一些属性来设置:

# set position in 2D simulation
L.atoms[0].position = (1.0, 0.0)

# set position in 3D simulation
L.atoms[0].position = (1.0, 0.0, 1.)
评估热数据

每次模拟运行通常都会根据系统状态、计算、修复或变量产生热输出。运行后可以通过 L.runs 列表查询这些值的轨迹。该列表包含不断增长的运行数据列表。第一个元素是第一次运行的输出,第二个元素是第二次运行的输出。

L.run(1000)
L.runs[0] # data of first 1000 time steps

L.run(1000)
L.runs[1] # data of second 1000 time steps

每次运行都包含所有轨迹的字典。每个轨迹都可以通过其热名称访问:

L.runs[0].thermo.Step # list of time steps in first run
L.runs[0].thermo.Ke   # list of kinetic energy values in first run

与 matplotlib 一起,从 LAMMPS 中绘制数据变得简单:

import matplotlib.plot as plt
steps = L.runs[0].thermo.Step
ke    = L.runs[0].thermo.Ke
plt.plot(steps, ke)
使用 PyLammps 进行错误处理

在 LAMMPS 中使用 C++ 异常来处理错误,可以在 C++ 端捕获它们并在 Python 端重新抛出它们。这样就可以通过Python异常处理机制来处理LAMMPS错误。

注:在Python中捕获LAMMPS异常仍然可能意味着当前LAMMPS进程处于非法状态,必须终止。建议尽快保存数据并终止Python实例。

在 IPython 笔记本和 Jupyter 中使用 PyLammps

如果为与 IPython 相同的 Python 解释器安装 LAMMPS Python 包,则可以直接在 Jupyter 内的 IPython 笔记本中使用 PyLammps。 Jupyter 是一个强大的集成开发环境 (IDE),适用于 Python、Julia 等许多动态语言,可在任何 Web 浏览器中运行。除了自动完成和语法突出显示之外,它还允许您使用与可执行 Python 代码混合的标记、数学公式、图形和动画来创建格式化文档。这是一种很棒的教程和展示您最新研究的格式。

要启动 Jupyter 实例,只需在 Python 环境中运行以下命令(假设您遵循快速入门说明):

jupyter notebook
IPyLamps 示例

IPython 笔记本的示例可以在 python/examples/pylammps 子目录中找到。要打开这些笔记本,请在此目录中启动 jupyter 笔记本并导航到其中之一。如果您编译并安装了 LAMMPS 共享库(但有例外)、PNG、JPEG 和 FFMPEG 支持,您应该能够重新运行所有这些笔记本。

验证二面势

此示例展示了如何使用 IPython Notebook 将调和二面势的简单 LAMMPS 模拟与其解析解进行比较。模拟中放置了四个原子,并使用数据文件对它们应用二面势。然后,通过从 Python 设置其位置,使其中一个原子沿中心轴旋转,从而改变二面角。

phi = [d \* math.pi / 180 for d in range(360)]

pos = [(1.0, math.cos(p), math.sin(p)) for p in phi]

pe = []
for p in pos:
    L.atoms[3].position = p
    L.run(0)
    pe.append(L.eval("pe"))

通过评估每个位置的势能,我们可以用解析公式验证该轨迹。为了比较这两个解决方案,我们使用 matplotlib 绘制了两个轨迹,它将生成的绘图嵌入到 IPython 笔记本中。

运行蒙特卡罗弛豫

第二个示例展示了如何使用 PyLammps 创建 2D 蒙特卡罗松弛模拟、计算和绘制能量项,甚至嵌入视频输出。

最初,二维系统是在能量最小的状态下创建的。

然后通过将每个原子移动随机增量来使其无序。

random.seed(27848)
deltaperturb = 0.2

for i in range(L.system.natoms):
    x, y = L.atoms[i].position
    dx = deltaperturb \* random.uniform(-1, 1)
    dy = deltaperturb \* random.uniform(-1, 1)
    L.atoms[i].position  = (x+dx, y+dy)

L.run(0)

最后用Python实现了蒙特卡罗算法。它以随机增量连续移动随机原子,并且只接受某些移动。

estart = L.eval("pe")
elast = estart

naccept = 0
energies = [estart]

niterations = 3000
deltamove = 0.1
kT = 0.05

natoms = L.system.natoms

for i in range(niterations):
    iatom = random.randrange(0, natoms)
    current_atom = L.atoms[iatom]

    x0, y0 = current_atom.position

    dx = deltamove \* random.uniform(-1, 1)
    dy = deltamove \* random.uniform(-1, 1)

    current_atom.position = (x0+dx, y0+dy)

    L.run(1, "pre no post no")

    e = L.eval("pe")
    energies.append(e)

    if e <= elast:
        naccept += 1
        elast = e
    elif random.random() <= math.exp(natoms\*(elast-e)/kT):
        naccept += 1
        elast = e
    else:
        current_atom.position = (x0, y0)

每次迭代的能量都收集在 Python 列表中,并最终使用 matplotlib 进行绘制。

IPython 笔记本还展示了如何使用转储命令以及如何在 IPython 笔记本中嵌入视频文件。

使用 PyLammps 和 mpi4py(实验)

PyLammps 可以使用 mpi4py 并行运行。这个 python 包可以使用以下命令安装

pip install mpi4py

以下是一个简短的示例,它读取现有的 LAMMPS 输入文件并并行执行。您可以在 example/melt 文件夹中找到 in.melt。

from mpi4py import MPI
from lammps import PyLammps

L = PyLammps()
L.file("in.melt")

if MPI.COMM_WORLD.rank == 0:
    print("Potential energy: ", L.eval("pe"))

MPI.Finalize()

要使用 4 个 MPI 进程并行运行此脚本 (melt.py),我们调用以下 mpirun 命令:

mpirun -np 4 python melt.py

注:任何命令都必须由所有 MPI 进程执行。但是,评估和查询系统状态仅适用于等级 0。

如果您发现此 Python 界面有用,请随时向 Richard Berger (richardberger@outlookcom) 提供有关如何改进它的反馈和想法。我们还希望鼓励人们编写教程风格的 IPython 笔记本,展示 LAMMPS 的使用情况,也许还有他们的最新研究成果。

  • 31
    点赞
  • 9
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
### 回答1: 以下是LAMMPS在Windows 10上的安装教程: 1. 下载和安装MinGW-w64 MinGW-w64是一个Windows上的GNU编译器集合,可以用于编译LAMMPS。您可以从MinGW-w64官网下载安装程序,并按照提示进行安装。 2. 下载和解压LAMMPS 您可以从LAMMPS官网下载最新版本的LAMMPS源代码,并将其解压到您选择的目录中。 3. 配置和编译LAMMPS 打开命令提示符窗口,并进入LAMMPS源代码目录。运行以下命令: mingw32-make yes-all mingw32-make -j4 这将配置和编译LAMMPS。请注意,-j4选项指定使用4个线程进行编译,您可以根据需要进行调整。 4. 运行LAMMPS 编译完成后,您可以在命令提示符窗口中运行LAMMPS。进入LAMMPS源代码目录,并运行以下命令: lmp_serial.exe -in input_file.in 其中,input_file.in是您的输入文件。您可以根据需要进行修改。 希望这个教程能够帮助您在Windows 10上安装LAMMPS。如果您有任何问题,请随时联系我们。 ### 回答2: LAMMPS是一种分子动力学软件,在物理化学、材料科学等领域广泛应用。本文将介绍在Windows10平台上安装LAMMPS的步骤。 安装步骤如下: 1. 下载软件 首先需要下载LAMMPS软件,可以从官网(https://lammps.sandia.gov/download.html)或者GitHub(https://github.com/lammps/lammps)下载安装包。下载完成后,可以将其解压到任意目录(比如D盘根目录下),以便安装时使用。 2. 安装依赖库 LAMMPS在Windows系统上需要一些依赖库的支持,比如MSYS2、OpenMPI和CMake等。这些依赖库需要在安装LAMMPS前先行安装。可以分别从官网上下载安装包,安装过程可能因系统版本、硬件配置等因素而有差异,在安装过程中需要耐心调试。 3. 编译安装 安装依赖库后,可以开始编译安装LAMMPS了。首先需要进入LAMMPS的安装目录(假设是D:\lammps-30Aug19)下,在该目录下创建build文件夹,用于存放编译后的LAMMPS程序。进入build目录后,使用以下命令进行编译: cmake -C ../cmake/presets/most.cmake ../src make mpi -j4 以上命令分别进行了两个步骤,第一个命令使用了CMake工具来进行配置,指定了软件包的位置,第二个命令是进行编译,其中-j4参数指定了使用4个核心(若电脑核心数不足,也可以设置为2,或者不设置,自动分配核心)进行编译。 4. 测试运行 编译后的LAMMPS程序位于build目录下,名为lmp_mpi.exe。运行该程序时需要指定输入文件,以及可选的各种参数。如以下命令: lmp_mpi.exe -in input.file 其中input.file是LAMMPS输入文件的名称。若执行不成功,可以检查输入文件和其它参数是否正确,或者重新进行编译。 总结 通过以上步骤进行LAMMPS在Windows10系统下的安装和编译,即可通过CMD或者PowerShell进行LAMMPS程序的运行。由于LAMMPS的编译过程与依赖库以及系统有关,因此可能会有误差,需要根据自身实际情况进行调试。 ### 回答3: LAMMPS是一个免费的分子动力学模拟程序,可以用于研究各种不同材料的物理和化学行为。在Windows 10上安装LAMMPS需要完成以下步骤: 1.安装cygwin Cygwin是一个允许在Windows机器上运行UNIX软件的开源工具。在Windows 10上安装cygwin需要先下载cygwin安装包,然后根据安装向导完成安装过程。在安装过程中,需要选择安装以下必要软件包:mingw64-x86_64-gcc-core、make、wget、curl、git、python、perl和OpenMPI等。并将cygwin添加到系统环境变量中。 2.安装编译器 安装完cygwin后,需要在cygwin上安装编译器,以便能够编译LAMMPS。可以通过运行以下命令安装编译器: apt-cyg install gcc-g++ 3.下载LAMMPS源代码 LAMMPS可以从其官方网站下载。选择最新版本的LAMMPS源代码,并下载到本地。 4.编译和安装LAMMPS 进入LAMMPS源代码目录,依照LAMMPS的安装指南进行编译和安装。在进行编译之前,需要根据所需的功能和功能组件配置MAKEFILE文件,以便进行定制化编译。在编译之前,还需要配置LAMMPS环境变量。此外,在编译过程中,还需确保系统配置正确,并且没有缺失任何必要的文件和库。一旦编译成功,就可以将LAMMPS添加到系统环境变量中,以便于使用。 5.测试LAMMPS 为了确保LAMMPS安装成功,在安装LAMMPS之后,需要对其进行测试。选择LAMMPS源代码目录下的test目录,依照README文件中的说明,进行测试。如果测试无误,则说明LAMMPS已经顺利安装并可以正常使用。 总之,在Windows 10上安装LAMMPS需要下载cygwin和LAMMPS源代码,配置编译器、环境变量和MAKEFILE文件,并进行编译、安装和测试等步骤。如果以上步骤都按照指导完成,则可以在Windows 10上成功安装LAMMPS

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值