【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中使用pytz库处理时区信息,包括获取当前UTC时间,转换为特定时区等,有需要的小伙伴可以参考一下... 目录时区基本信息python datetime使用timezonepandas处理时区数据知识延展时区基本信息

Python位移操作和位运算的实现示例

《Python位移操作和位运算的实现示例》本文主要介绍了Python位移操作和位运算的实现示例,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一... 目录1. 位移操作1.1 左移操作 (<<)1.2 右移操作 (>>)注意事项:2. 位运算2.1

如何在 Spring Boot 中实现 FreeMarker 模板

《如何在SpringBoot中实现FreeMarker模板》FreeMarker是一种功能强大、轻量级的模板引擎,用于在Java应用中生成动态文本输出(如HTML、XML、邮件内容等),本文... 目录什么是 FreeMarker 模板?在 Spring Boot 中实现 FreeMarker 模板1. 环

Qt实现网络数据解析的方法总结

《Qt实现网络数据解析的方法总结》在Qt中解析网络数据通常涉及接收原始字节流,并将其转换为有意义的应用层数据,这篇文章为大家介绍了详细步骤和示例,感兴趣的小伙伴可以了解下... 目录1. 网络数据接收2. 缓冲区管理(处理粘包/拆包)3. 常见数据格式解析3.1 jsON解析3.2 XML解析3.3 自定义

使用Python和Pyecharts创建交互式地图

《使用Python和Pyecharts创建交互式地图》在数据可视化领域,创建交互式地图是一种强大的方式,可以使受众能够以引人入胜且信息丰富的方式探索地理数据,下面我们看看如何使用Python和Pyec... 目录简介Pyecharts 简介创建上海地图代码说明运行结果总结简介在数据可视化领域,创建交互式地

SpringMVC 通过ajax 前后端数据交互的实现方法

《SpringMVC通过ajax前后端数据交互的实现方法》:本文主要介绍SpringMVC通过ajax前后端数据交互的实现方法,本文给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价... 在前端的开发过程中,经常在html页面通过AJAX进行前后端数据的交互,SpringMVC的controll

Spring Security自定义身份认证的实现方法

《SpringSecurity自定义身份认证的实现方法》:本文主要介绍SpringSecurity自定义身份认证的实现方法,下面对SpringSecurity的这三种自定义身份认证进行详细讲解,... 目录1.内存身份认证(1)创建配置类(2)验证内存身份认证2.JDBC身份认证(1)数据准备 (2)配置依

利用python实现对excel文件进行加密

《利用python实现对excel文件进行加密》由于文件内容的私密性,需要对Excel文件进行加密,保护文件以免给第三方看到,本文将以Python语言为例,和大家讲讲如何对Excel文件进行加密,感兴... 目录前言方法一:使用pywin32库(仅限Windows)方法二:使用msoffcrypto-too

C#使用StackExchange.Redis实现分布式锁的两种方式介绍

《C#使用StackExchange.Redis实现分布式锁的两种方式介绍》分布式锁在集群的架构中发挥着重要的作用,:本文主要介绍C#使用StackExchange.Redis实现分布式锁的... 目录自定义分布式锁获取锁释放锁自动续期StackExchange.Redis分布式锁获取锁释放锁自动续期分布式

springboot使用Scheduling实现动态增删启停定时任务教程

《springboot使用Scheduling实现动态增删启停定时任务教程》:本文主要介绍springboot使用Scheduling实现动态增删启停定时任务教程,具有很好的参考价值,希望对大家有... 目录1、配置定时任务需要的线程池2、创建ScheduledFuture的包装类3、注册定时任务,增加、删