【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

相关文章

Python基础文件操作方法超详细讲解(详解版)

《Python基础文件操作方法超详细讲解(详解版)》文件就是操作系统为用户或应用程序提供的一个读写硬盘的虚拟单位,文件的核心操作就是读和写,:本文主要介绍Python基础文件操作方法超详细讲解的相... 目录一、文件操作1. 文件打开与关闭1.1 打开文件1.2 关闭文件2. 访问模式及说明二、文件读写1.

基于SpringBoot+Mybatis实现Mysql分表

《基于SpringBoot+Mybatis实现Mysql分表》这篇文章主要为大家详细介绍了基于SpringBoot+Mybatis实现Mysql分表的相关知识,文中的示例代码讲解详细,感兴趣的小伙伴可... 目录基本思路定义注解创建ThreadLocal创建拦截器业务处理基本思路1.根据创建时间字段按年进

Python将博客内容html导出为Markdown格式

《Python将博客内容html导出为Markdown格式》Python将博客内容html导出为Markdown格式,通过博客url地址抓取文章,分析并提取出文章标题和内容,将内容构建成html,再转... 目录一、为什么要搞?二、准备如何搞?三、说搞咱就搞!抓取文章提取内容构建html转存markdown

Python获取中国节假日数据记录入JSON文件

《Python获取中国节假日数据记录入JSON文件》项目系统内置的日历应用为了提升用户体验,特别设置了在调休日期显示“休”的UI图标功能,那么问题是这些调休数据从哪里来呢?我尝试一种更为智能的方法:P... 目录节假日数据获取存入jsON文件节假日数据读取封装完整代码项目系统内置的日历应用为了提升用户体验,

SpringBoot3实现Gzip压缩优化的技术指南

《SpringBoot3实现Gzip压缩优化的技术指南》随着Web应用的用户量和数据量增加,网络带宽和页面加载速度逐渐成为瓶颈,为了减少数据传输量,提高用户体验,我们可以使用Gzip压缩HTTP响应,... 目录1、简述2、配置2.1 添加依赖2.2 配置 Gzip 压缩3、服务端应用4、前端应用4.1 N

SpringBoot实现数据库读写分离的3种方法小结

《SpringBoot实现数据库读写分离的3种方法小结》为了提高系统的读写性能和可用性,读写分离是一种经典的数据库架构模式,在SpringBoot应用中,有多种方式可以实现数据库读写分离,本文将介绍三... 目录一、数据库读写分离概述二、方案一:基于AbstractRoutingDataSource实现动态

Python FastAPI+Celery+RabbitMQ实现分布式图片水印处理系统

《PythonFastAPI+Celery+RabbitMQ实现分布式图片水印处理系统》这篇文章主要为大家详细介绍了PythonFastAPI如何结合Celery以及RabbitMQ实现简单的分布式... 实现思路FastAPI 服务器Celery 任务队列RabbitMQ 作为消息代理定时任务处理完整

使用Jackson进行JSON生成与解析的新手指南

《使用Jackson进行JSON生成与解析的新手指南》这篇文章主要为大家详细介绍了如何使用Jackson进行JSON生成与解析处理,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 目录1. 核心依赖2. 基础用法2.1 对象转 jsON(序列化)2.2 JSON 转对象(反序列化)3.

Python Websockets库的使用指南

《PythonWebsockets库的使用指南》pythonwebsockets库是一个用于创建WebSocket服务器和客户端的Python库,它提供了一种简单的方式来实现实时通信,支持异步和同步... 目录一、WebSocket 简介二、python 的 websockets 库安装三、完整代码示例1.

揭秘Python Socket网络编程的7种硬核用法

《揭秘PythonSocket网络编程的7种硬核用法》Socket不仅能做聊天室,还能干一大堆硬核操作,这篇文章就带大家看看Python网络编程的7种超实用玩法,感兴趣的小伙伴可以跟随小编一起... 目录1.端口扫描器:探测开放端口2.简易 HTTP 服务器:10 秒搭个网页3.局域网游戏:多人联机对战4.