【ASE+python学习】批量实现将含有空位的结构进行氢饱和(实际就是在某个位置添加氢原子)

本文主要是介绍【ASE+python学习】批量实现将含有空位的结构进行氢饱和(实际就是在某个位置添加氢原子),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

批量实现将含有空位的结构进行氢饱和

  • 任务大纲
  • 任务思路
  • 氢原子相对位置计算
  • 代码实现-单个结构加氢
    • 导入所需库/包
    • 计算氢原子坐标
    • 给结构加氢
    • 运行结果
  • 代码实现-批量结构加氢
    • 导入相关库/包
    • 遍历文件夹内所有结构,并逐次给结构加氢
    • 批量实现的逻辑与结果

在搬砖过程中,我经常要对大批量的含有空位的结构进行氢饱和处理。如果一个个打开结构,然后进行加氢,我估计要累死,对颈椎及其不友好。而且这种简单且没有技术含量的工作,我想可以用更科学的方式进行实现。

任务大纲

本文以含有一个空位的石墨烯团簇为例,通过ASE与python程序对该结构进行加氢。单个结构加氢的任务大纲见下,如果是多个任务,结合for循环执行即可。
在这里插入图片描述

任务思路

实现给结构加氢的任务思路为:

  1. 首先需要基于已有的已经进行加氢的结构,计算在空位处氢原子与三个碳原子的相对位置信息,包括氢原子偏离结构平面的距离、与碳相连的键长、以及氢原子与结构平面的夹角、氢原子与碳相连的键在结构平面的投影在x方向与y方向的夹角(大概计算就行,后期会进行结构优化,优化氢原子的位置)
  2. 接着使用ASE读取结构信息,获取空位处三个碳原子的坐标信息
  3. 基于三个碳原子的坐标信息,计算三个氢原子的坐标信息
  4. 然后通过ASE将三个氢原子逐次添加到结构
  5. 将新结构写入新的文件中

氢原子相对位置计算

这里计算的参数大概就行,只要能让氢与相应的碳原子连接在一起,后续进行结构优化会优化到实际位置
在这里插入图片描述

代码实现-单个结构加氢

导入所需库/包

import os
import math
import shutil
from ase.io import read,write
from ase.visualize import view
from ase.build import molecule
from ase.build import add_adsorbate

计算氢原子坐标

def addHsite(C95Vac1_path=r'D:\software output files\quntum_experiment\POSCAR\non_periodic_CN\C96_molecule\4-C95Vac1_str_1e-2',C95Vac1_file='POSCAR9'):'''该函数用于获取空位附近三个碳原子的位置,并计算出三个氢的位置思路设计:读入含有空位的结构获取空位处三个碳原子的位置在MS打开含有空位且已进行氢饱和的结构,量碳氢键长,计算该键与xy平面的夹角,以及在xy平面的投影接着根据空位的排布,计算碳碳氢二面角,从而计算碳氢键长在xy平面的投影在x方向与y方向的分量基于空位三个碳原子的坐标,计算三个氢原子的坐标'''C95Vac1=read(os.path.join(C95Vac1_path,C95Vac1_file),format='vasp')#读入结构信息,转为atoms objectview(C95Vac1)#获取空位碳原子的位置,仅提取x,y方向position_62=C95Vac1[62].position[:2]position_61=C95Vac1[61].position[:2]position_114=C95Vac1[114].position[:2]#计算x,z方向偏移量z=1.06*math.sin(49*math.pi/180)#1.06为氢键的长度,这里计算的是氢原子偏离结构平面的高度x=1.06*math.cos(49*math.pi/180)#这里计算的是氢原子与碳原子所成的键,在xy平面的投影长度#分别计算与C62,C61,C114相连氢原子的位置#C62-H位置计算position_62_newx=[i+x for i in position_62 if i ==position_62[0]]#计算在x方向氢原子与C62位置的偏移量position_62_new=[]#建立空列表用于存储氢原子在平面的坐标position_62_new.append(position_62_newx[0])position_62_new.append(position_62[1])#计算C61与C114的氢原子在x,y方向的偏移量xx=x*math.cos(71*math.pi/180)#71为二面角xy=x*math.sin(71*math.pi/180)#C61-H位置计算position_61_newx=[i-xx for i in position_61 if i ==position_61[0]]position_61_newy=[i+xy for i in position_61 if i ==position_61[1]]position_61_new=[]position_61_new.append(position_61_newx[0])position_61_new.append(position_61_newy[0])#C114-H位置计算position_114_newx=[i-xx for i in position_114 if i ==position_114[0]]position_114_newy=[i-xy for i in position_114 if i ==position_114[1]]position_114_new=[]position_114_new.append(position_114_newx[0])position_114_new.append(position_114_newy[0])return x,z,position_62_new,position_61_new,position_114_new

给结构加氢

def str_addH(C95Vac1_path=r'D:\software output files\quntum_experiment\POSCAR\non_periodic_CN\C96_molecule\4-C95Vac1_str_1e-2',C95Vac1_file='POSCAR9'):'''该函数用于给结构的空位位置附件添加氢原子思路:基于函数addHsite()所得到的三个氢原子的坐标,依次添加到结构表面,接着将修改后的结构写入新的POSCAR文件中'''x,z,position_62_new,position_61_new,position_114_new=addHsite()H=molecule('H')C95Vac1=read(os.path.join(C95Vac1_path,C95Vac1_file),format='vasp')#读入结构信息,转为atoms objectadd_adsorbate(C95Vac1,H,position=position_62_new,height = z)#在结构上依次添加氢原子add_adsorbate(C95Vac1,H,position=position_61_new,height = z)add_adsorbate(C95Vac1,H,position=position_114_new,height = z)write(os.path.join(C95Vac1_path,C95Vac1_file+'_addH'),C95Vac1,format='vasp')#将添加了氢原子的结构写入文件view(C95Vac1)

运行结果

调用函数运行

str_addH()

运行结果
在这里插入图片描述

代码实现-批量结构加氢

批量给结构加氢,实际上只需要在读入结构那部分代码改为使用for循环遍历文件夹中的结构即可。

导入相关库/包

import os
import math
import shutil
from ase.io import read,write
from ase.visualize import view
from ase.build import molecule
from ase.build import add_adsorbate

遍历文件夹内所有结构,并逐次给结构加氢

def batchstr_addH(file_path,save_str,save_folder):'''该函数实现批量读入含有空位的结构,并将三个氢原子添加到空位附近,实现氢饱和,并将新结构写入文件,存入新文件夹file_path:结构文件所处的文件夹路径save_str:添加氢原子后的结构存储文件夹所在路径save_folder:添加氢原子后的结构存储文件夹的名字'''save_path=os.path.join(save_str,save_folder)str_files=os.listdir(file_path)x,z,position_62_new,position_61_new,position_114_new=addHsite()H=molecule('H')if os.path.exists(save_path):shutil.rmtree(save_path)if not os.path.exists(save_path):os.makedirs(save_path)for str_file in str_files:#遍历结构文件structure=read(os.path.join(file_path,str_file),format='vasp')add_adsorbate(structure,H,position=position_62_new,height = z)#在结构上依次添加氢原子add_adsorbate(structure,H,position=position_61_new,height = z)add_adsorbate(structure,H,position=position_114_new,height = z)write(os.path.join(save_path,str_file),structure,format='vasp')
file_path=r'D:\software output files\quntum_experiment\POSCAR\non_periodic_CN\C96_molecule\5-C94Vac1N1_str_1e-2(based_on_poscar9)_ASE\Npe'
save_str=r'D:\software output files\quntum_experiment\POSCAR\non_periodic_CN\C96_molecule\5-C94Vac1N1_str_1e-2(based_on_poscar9)_ASE'
save_folder='Npe_addH'
batchstr_addH(file_path,save_str,save_folder)

批量实现的逻辑与结果

在这里插入图片描述
如此,便实现了自动批量加氢的操作。

这篇关于【ASE+python学习】批量实现将含有空位的结构进行氢饱和(实际就是在某个位置添加氢原子)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/306099

相关文章

HarmonyOS学习(七)——UI(五)常用布局总结

自适应布局 1.1、线性布局(LinearLayout) 通过线性容器Row和Column实现线性布局。Column容器内的子组件按照垂直方向排列,Row组件中的子组件按照水平方向排列。 属性说明space通过space参数设置主轴上子组件的间距,达到各子组件在排列上的等间距效果alignItems设置子组件在交叉轴上的对齐方式,且在各类尺寸屏幕上表现一致,其中交叉轴为垂直时,取值为Vert

Ilya-AI分享的他在OpenAI学习到的15个提示工程技巧

Ilya(不是本人,claude AI)在社交媒体上分享了他在OpenAI学习到的15个Prompt撰写技巧。 以下是详细的内容: 提示精确化:在编写提示时,力求表达清晰准确。清楚地阐述任务需求和概念定义至关重要。例:不用"分析文本",而用"判断这段话的情感倾向:积极、消极还是中性"。 快速迭代:善于快速连续调整提示。熟练的提示工程师能够灵活地进行多轮优化。例:从"总结文章"到"用

python: 多模块(.py)中全局变量的导入

文章目录 global关键字可变类型和不可变类型数据的内存地址单模块(单个py文件)的全局变量示例总结 多模块(多个py文件)的全局变量from x import x导入全局变量示例 import x导入全局变量示例 总结 global关键字 global 的作用范围是模块(.py)级别: 当你在一个模块(文件)中使用 global 声明变量时,这个变量只在该模块的全局命名空

【前端学习】AntV G6-08 深入图形与图形分组、自定义节点、节点动画(下)

【课程链接】 AntV G6:深入图形与图形分组、自定义节点、节点动画(下)_哔哩哔哩_bilibili 本章十吾老师讲解了一个复杂的自定义节点中,应该怎样去计算和绘制图形,如何给一个图形制作不间断的动画,以及在鼠标事件之后产生动画。(有点难,需要好好理解) <!DOCTYPE html><html><head><meta charset="UTF-8"><title>06

学习hash总结

2014/1/29/   最近刚开始学hash,名字很陌生,但是hash的思想却很熟悉,以前早就做过此类的题,但是不知道这就是hash思想而已,说白了hash就是一个映射,往往灵活利用数组的下标来实现算法,hash的作用:1、判重;2、统计次数;

hdu1043(八数码问题,广搜 + hash(实现状态压缩) )

利用康拓展开将一个排列映射成一个自然数,然后就变成了普通的广搜题。 #include<iostream>#include<algorithm>#include<string>#include<stack>#include<queue>#include<map>#include<stdio.h>#include<stdlib.h>#include<ctype.h>#inclu

【C++】_list常用方法解析及模拟实现

相信自己的力量,只要对自己始终保持信心,尽自己最大努力去完成任何事,就算事情最终结果是失败了,努力了也不留遗憾。💓💓💓 目录   ✨说在前面 🍋知识点一:什么是list? •🌰1.list的定义 •🌰2.list的基本特性 •🌰3.常用接口介绍 🍋知识点二:list常用接口 •🌰1.默认成员函数 🔥构造函数(⭐) 🔥析构函数 •🌰2.list对象

【Prometheus】PromQL向量匹配实现不同标签的向量数据进行运算

✨✨ 欢迎大家来到景天科技苑✨✨ 🎈🎈 养成好习惯,先赞后看哦~🎈🎈 🏆 作者简介:景天科技苑 🏆《头衔》:大厂架构师,华为云开发者社区专家博主,阿里云开发者社区专家博主,CSDN全栈领域优质创作者,掘金优秀博主,51CTO博客专家等。 🏆《博客》:Python全栈,前后端开发,小程序开发,人工智能,js逆向,App逆向,网络系统安全,数据分析,Django,fastapi

让树莓派智能语音助手实现定时提醒功能

最初的时候是想直接在rasa 的chatbot上实现,因为rasa本身是带有remindschedule模块的。不过经过一番折腾后,忽然发现,chatbot上实现的定时,语音助手不一定会有响应。因为,我目前语音助手的代码设置了长时间无应答会结束对话,这样一来,chatbot定时提醒的触发就不会被语音助手获悉。那怎么让语音助手也具有定时提醒功能呢? 我最后选择的方法是用threading.Time

【Python编程】Linux创建虚拟环境并配置与notebook相连接

1.创建 使用 venv 创建虚拟环境。例如,在当前目录下创建一个名为 myenv 的虚拟环境: python3 -m venv myenv 2.激活 激活虚拟环境使其成为当前终端会话的活动环境。运行: source myenv/bin/activate 3.与notebook连接 在虚拟环境中,使用 pip 安装 Jupyter 和 ipykernel: pip instal