def slab(self, name):
from pymatgen.analysis.adsorption import AdsorbateSiteFinder, plot_slab, reorient_z
from pymatgen.core.surface import Slab, SlabGenerator, generate_all_slabs, Structure, Lattice, ReconstructionGenerator
from pymatgen.symmetry.analyzer import SpacegroupAnalyzer
from pymatgen.core.structure import Structure
from pymatgen.ext.matproj import MPRester
from matplotlib import pyplot as plt
from pymatgen.io.vasp.inputs import Poscar
from pymatgen.io.vasp.sets import MVLSlabSet
import shutil
import os
os.chdir(self.dire)
print(os.getcwd())
poscar = Poscar.from_file("POSCAR")
relax = poscar.structure
custom_settings = {
"NPAR": 4,
'ICHARG': 2,
} # 用户的INCAR 设置
relax = MVLSlabSet(relax, user_incar_settings=custom_settings)
#Vasp输入文件生成器
dire = str(name)
#设置一个用作存储输入文件的名称
relax.write_input(dire) #将生成的VASP输入文件写入存储
diree = os.chdir("./" + dire)
#定义一个更改当前目录的变量
dire2 = './vaspstd_sub'
#确立脚本名称
shutil.copy(r"C:\Users\41958\.spyder-py3\vaspstd_sub", dire2)
#将脚本写入VASP输入文件所在文件夹
# os.chdir("../")
#将当前目录改为默认目录
os.chdir(r"D:\Desktop\VASP practical\workdir")
print(os.getcwd())
print('finished')
def get_wf_surface(slabs,
molecules=[],
bulk_structure=None,
slab_gen_params=None,
vasp_cmd="vasp",
db_file=None,
ads_structures_params={},
add_molecules_in_box=False):
"""
Args:
slabs (list of Slabs or Structures): slabs to calculate
molecules (list of Molecules): molecules to place as adsorbates
bulk_structure (Structure): bulk structure from which generate slabs
after reoptimization. If supplied, workflow will begin with
bulk structure optimization.
slab_gen_params (dict): dictionary of slab generation parameters
used to generate the slab, necessary to get the slab
that corresponds to the bulk structure if in that mode
ads_structures_params (dict): parameters to be supplied as
kwargs to AdsorbateSiteFinder.generate_adsorption_structures
add_molecules_in_box (boolean): flag to add calculation of molecule
energies to the workflow
db_file (string): path to database file
vasp_cmd (string): vasp command
Returns:
Workflow
"""
fws, parents = [], []
if bulk_structure:
vis = MVLSlabSet(bulk_structure, bulk=True)
fws.append(
OptimizeFW(bulk_structure,
vasp_input_set=vis,
vasp_cmd="vasp",
db_file=db_file))
parents = fws[0]
for slab in slabs:
name = slab.composition.reduced_formula
if getattr(slab, "miller_index", None):
name += "_{}".format(slab.miller_index)
fws.append(
get_slab_fw(slab,
bulk_structure,
slab_gen_params,
db_file=db_file,
vasp_cmd=vasp_cmd,
parents=parents,
name=name + " slab optimization"))
for molecule in molecules:
ads_slabs = AdsorbateSiteFinder(
slab).generate_adsorption_structures(molecule,
**ads_structures_params)
for n, ads_slab in enumerate(ads_slabs):
ads_name = "{}-{} adsorbate optimization {}".format(
molecule.composition.formula, name, n)
fws.append(
get_slab_fw(ads_slab,
bulk_structure,
slab_gen_params,
db_file=db_file,
vasp_cmd=vasp_cmd,
parents=parents,
name=ads_name))
if add_molecules_in_box:
for molecule in molecules:
# molecule in box
m_struct = Structure(Lattice.cubic(10),
molecule.species_and_occu,
molecule.cart_coords,
coords_are_cartesian=True)
m_struct.translate_sites(
list(range(len(m_struct))),
np.array([0.5] * 3) - np.average(m_struct.frac_coords, axis=0))
vis = MVLSlabSet(m_struct)
fws.append(
OptimizeFW(molecule,
job_type="normal",
vasp_input_set=vis,
db_file=db_file,
vasp_cmd=vasp_cmd))
# TODO: add analysis framework
return Workflow(fws, name="")
def get_slab_fw(slab,
bulk_structure=None,
slab_gen_params={},
db_file=None,
vasp_input_set=None,
parents=None,
vasp_cmd="vasp",
name=""):
"""
Function to generate a a slab firework. Returns a TransmuterFW if bulk_structure is specified,
constructing the necessary transformations from the slab and slab generator parameters,
or an OptimizeFW if only a slab is specified.
Args:
slab (Slab or Structure): structure or slab corresponding
to the slab to be calculated
bulk_structure (Structure): bulk structure corresponding to slab, if
provided, slab firework is constructed as a TransmuterFW using
the necessary transformations to get the slab from the bulk
slab_gen_params (dict): dictionary of slab generation parameters
used to generate the slab, necessary to get the slab
that corresponds to the bulk structure
vasp_input_set (VaspInputSet): vasp_input_set corresponding to
the slab calculation
parents (Fireworks or list of ints): parent FWs
db_file (string): path to database file
vasp_cmd (string): vasp command
Returns:
Firework
"""
vasp_input_set = vasp_input_set or MVLSlabSet(slab)
# If a bulk_structure is specified, generate the set of transformations, else
# just create an optimize FW with the slab
if bulk_structure:
if not isinstance(slab, Slab):
raise ValueError(
"structure input to get_slab_fw requires slab to be a slab object!"
)
slab_trans_params = {
"miller_index": slab.miller_index,
"shift": slab.shift
}
slab_trans_params.update(slab_gen_params)
# Get supercell parameters
trans_struct = SlabTransformation(**slab_trans_params)
slab_from_bulk = trans_struct.apply_transformation(bulk_structure)
supercell_trans = SupercellTransformation.from_scaling_factors(
round(slab.lattice.a / slab_from_bulk.lattice.a),
round(slab.lattice.b / slab_from_bulk.lattice.b))
# Get adsorbates for InsertSitesTransformation
if "adsorbate" in slab.site_properties.get("surface_properties",
[None]):
ads_sites = [
site for site in slab
if site.properties["surface_properties"] == "adsorbate"
]
else:
ads_sites = []
transformations = [
"SlabTransformation", "SupercellTransformation",
"InsertSitesTransformation", "AddSitePropertyTransformation"
]
trans_params = [
slab_trans_params, {
"scaling_matrix": supercell_trans.scaling_matrix
}, {
"species": [site.species_string for site in ads_sites],
"coords": [site.frac_coords for site in ads_sites]
}, {
"site_properties": slab.site_properties
}
]
return TransmuterFW(name=name,
structure=bulk_structure,
transformations=transformations,
transformation_params=trans_params,
copy_vasp_outputs=True,
db_file=db_file,
vasp_cmd=vasp_cmd,
parents=parents,
vasp_input_set=vasp_input_set)
else:
return OptimizeFW(name=name,
structure=slab,
vasp_input_set=vasp_input_set,
vasp_cmd=vasp_cmd,
db_file=db_file,
parents=parents,
job_type="normal")
def sol(self,
EB_K_2,
miller_index_2,
min_slab_size_2=8.0,
min_vacuum_size_2=15):
from pymatgen.analysis.adsorption import AdsorbateSiteFinder, plot_slab, reorient_z
from pymatgen.core.surface import Slab, SlabGenerator, generate_all_slabs, Structure, Lattice, ReconstructionGenerator
from pymatgen.symmetry.analyzer import SpacegroupAnalyzer
from pymatgen.core.structure import Structure
from pymatgen.ext.matproj import MPRester
from matplotlib import pyplot as plt
from pymatgen.io.vasp.inputs import Poscar
from pymatgen.io.vasp.sets import MVLSlabSet
import shutil
import os
mpr = MPRester() #密钥
mp_id = self.mp_id #通过mp_id来索引结构
struct = mpr.get_structure_by_material_id(mp_id)
#获取结构信息
struct = SpacegroupAnalyzer(
struct).get_conventional_standard_structure()
#空间群分析
need_miller_index = miller_index_2 #通过米勒指数，确定要切的晶面
slab = SlabGenerator(struct, miller_index=need_miller_index, min_slab_size=min_slab_size_2,\
min_vacuum_size=min_vacuum_size_2, center_slab=True)
#晶面生成器参数
for n, slabs in enumerate(slab.get_slabs()):
slabs_bak = slabs.copy() #可能的晶面
slabs.make_supercell(self.supercell)
#晶胞扩充
A = Poscar(slabs) #将切面转换为Poscar
relax = A.structure #将Poscar 转换为结构信息
custom_settings = {"NPAR": 4} # 用户的INCAR 设置
relax = MVLSlabSet(relax, user_incar_settings=custom_settings)
#Vasp输入文件生成器
fig = plt.figure() #绘图--确立画布
ax = fig.add_subplot(111) #绘图--确立位置
plot_slab(slabs, ax, adsorption_sites=False) #绘图
dire = str(mp_id) + "---" + "sol" + str(EB_K_2) + str(
need_miller_index) + '----' + str(n)
#设置一个用作存储输入文件的名称
plt.savefig(dire) #将该名称用于保存图片
relax.write_input(dire) #将生成的VASP输入文件写入存储
os.chdir("./" + dire)
#定义一个更改当前目录的变量
dire2 = './vaspstd_sub'
#确立脚本名称
shutil.copy(r"C:\Users\41958\.spyder-py3\vaspstd_sub", dire2)
#将脚本写入VASP输入文件所在文件夹
eb = str(EB_K_2)
ls = str('TURE')
with open('INCAR', 'a') as file_object:
file_object.write('LSOL = ' + ls + '\n' + 'EB_K = ' + eb)
# os.chdir("../")
#将当前目录改为默认目录
os.chdir(r"D:\Desktop\VASP practical\workdir")
print(os.getcwd())
print('finished')
def fix_absorbed(self,
need_miller_index,
mole,
num,
selective_dynamic,
min_slab_size_1=8.0,
min_vacuum_size_1=15,
judge='fuchdi',
appendage=""):
from pymatgen import Structure, Lattice, MPRester, Molecule
import pymatgen.core.structure
import pymatgen.core.sites
from pymatgen.analysis.adsorption import AdsorbateSiteFinder, reorient_z, plot_slab
from pymatgen.core.surface import generate_all_slabs
from pymatgen.symmetry.analyzer import SpacegroupAnalyzer
from matplotlib import pyplot as plt
from pymatgen.ext.matproj import MPRester
from pymatgen.io.vasp.inputs import Poscar
from pymatgen.io.vasp.sets import MVLSlabSet
from pymatgen.io.cif import CifWriter
import os
import shutil
from openbabel import openbabel
from pymatgen.core.surface import Slab, SlabGenerator, generate_all_slabs, Structure, Lattice, ReconstructionGenerator
mp_id = self.mp_id
os.chdir(r"F:\VASP practical\Input")
print(os.getcwd())
# Note that you must provide your own API Key, which can
# be accessed via the Dashboard at materialsproject.org
mpr = MPRester()
struct = mpr.get_structure_by_material_id(mp_id)
struct = SpacegroupAnalyzer(
struct).get_conventional_standard_structure()
# fcc_ni = Structure.from_spacegroup("Fm-3m", Lattice.cubic(3.5), ["Ni", "Ni"],
# [[0, 0, 0], [0.5, 0.5, 0.5]])
slab = SlabGenerator(struct,
miller_index=need_miller_index,
min_slab_size=min_slab_size_1,
min_vacuum_size=min_vacuum_size_1,
center_slab=True)
for n, slabs in enumerate(slab.get_slabs()):
if str(n) in str(num):
slabs_bak = slabs.copy() #可能的晶面
slabs.make_supercell(self.supercell)
print(n)
#晶胞扩充
asf_ni_111 = AdsorbateSiteFinder(
slabs, selective_dynamics=selective_dynamic)
ads_sites = asf_ni_111.find_adsorption_sites()
# print(ads_sites)
assert len(ads_sites) == 4
fig0 = plt.figure()
ax = fig0.add_subplot(111)
plot_slab(slabs, ax, adsorption_sites=False)
fig1 = plt.figure()
ax = fig1.add_subplot(111)
os.chdir(r"D:\Desktop\VASP practical\Cif library")
print(os.getcwd())
obConversion = openbabel.OBConversion()
obConversion.SetInAndOutFormats("pdb", "gjf")
mol = openbabel.OBMol()
print(mol)
c = obConversion.ReadFile(mol, "CH3OH.pdb")
obConversion.WriteFile(mol, "CH3OH.pdb" + '1.gjf')
adsorbate = Molecule.from_file("CH3OH.pdb" + '.gjf')
os.chdir(r"F:\VASP practical\Input")
print(os.getcwd())
print(adsorbate.sites)
ads_structs = asf_ni_111.add_adsorbate(
adsorbate,
(20, 20, 20),
translate=False,
)
# ads_structs = asf_ni_111.generate_adsorption_structures(adsorbate,
# repeat=[1, 1, 1])
# A = Poscar(ads_structs[0])
A = Poscar(reorient_z(ads_structs)) #将切面转换为Poscar
open('POSCAR', 'w').write(str(A))
p = Poscar.from_file('POSCAR')
# w = CifWriter(A.struct)
# w.write_file('mystructure.cif')
path = r'F:\VASP practical\Input\POSCAR' # 文件路径
if os.path.exists(path): # 如果文件存在
# 删除文件，可使用以下两种方法。
os.remove(path)
#os.unlink(path)
else:
print('no such file:%s' % my_file) # 则返回文件不存在
# w = CifWriter(A.struct)
# w.write_file('mystructure.cif')
relax = p.structure #将Poscar 转换为结构信息
custom_settings = {"NPAR": 4} # 用户的INCAR 设置
relaxs = MVLSlabSet(relax, user_incar_settings=custom_settings)
# Vasp输入文件生成器
dire = str(mp_id) + str(selective_dynamic) + str(mole) + str(
need_miller_index).replace(" ", "") + str(n)
# print (relax)
relaxs.write_input(dire)
os.chdir("./" + dire)
print(os.getcwd())
fig0.savefig('slab.png',
bbox_inches='tight',
transparent=True,
dpi=600,
format='png')
plot_slab(ads_structs, ax, adsorption_sites=False, decay=0.09)
fig1.savefig('slab_adsobate.png',
bbox_inches='tight',
transparent=True,
dpi=600,
format='png')
#定义一个更改当前目录的变量
dire2 = './vaspstd_sub'
#确立脚本名称
shutil.copy(r"C:\Users\41958\.spyder-py3\vaspstd_sub", dire2)
eb = appendage #添加其他INCAR参数
with open('INCAR', 'r') as f1:
lines = f1.readlines()
with open('INCAR', 'w') as f2:
for line in lines:
if judge in line:
continue
f2.write(line)
with open('INCAR', 'a') as f3:
f3.write(eb)
# open('POSCAR001', 'w').write(str(Poscar(reorient_z(ads_structs[0]))))
os.chdir(r"D:\Desktop\VASP practical\workdir")
print(os.getcwd())
print('finished')
# my_lattace = Lattace('mp-698074')#半水石膏
# # my_lattace.phase_out()#生成晶胞优化的输入文件
# go = my_lattace.phase_sol(66,judge='LWAVE', appendage= '\nLWAVE = Ture')
# print('yoo')
def slab(self, miller_index_1, min_slab_size_1=8.0, min_vacuum_size_1=15):
from pymatgen.analysis.adsorption import AdsorbateSiteFinder, plot_slab, reorient_z
from pymatgen.core.surface import Slab, SlabGenerator, generate_all_slabs, Structure, Lattice, ReconstructionGenerator
from pymatgen.symmetry.analyzer import SpacegroupAnalyzer
from pymatgen.core.structure import Structure
from pymatgen.ext.matproj import MPRester
from matplotlib import pyplot as plt
from pymatgen.io.vasp.inputs import Poscar
from pymatgen.io.vasp.sets import MVLSlabSet
import shutil
import os
mpr = MPRester() #密钥
mp_id = self.mp_id #通过mp_id来索引结构
struct = mpr.get_structure_by_material_id(mp_id)
#获取结构信息
struct = SpacegroupAnalyzer(
struct).get_conventional_standard_structure()
#空间群分析
need_miller_index = miller_index_1 #通过米勒指数，确定要切的晶面
slab = SlabGenerator(struct, miller_index=need_miller_index, min_slab_size=min_slab_size_1,\
min_vacuum_size=min_vacuum_size_1, center_slab=True)
#晶面生成器参数
gh = str(miller_index_1).replace(" ", "")
for n, slabs in enumerate(slab.get_slabs()):
slabs_bak = slabs.copy() #可能的晶面
slabs.make_supercell(self.supercell)
#晶胞扩充
cc = slabs.surface_area
os.chdir(r"F:\VASP practical\Input")
print(os.getcwd())
print(n)
A = Poscar(slabs) #将切面转换为Poscar
relax = A.structure #将Poscar 转换为结构信息
custom_settings = {"NPAR": 4} # 用户的INCAR 设置
relax = MVLSlabSet(relax, user_incar_settings=custom_settings)
#Vasp输入文件生成器
fig = plt.figure() #绘图--确立画布
ax = fig.add_subplot(111) #绘图--确立位置
plot_slab(slabs, ax, adsorption_sites=False) #绘图
dire = str(mp_id) + "---" + str(gh) + '----' + str(n)
#设置一个用作存储输入文件的名称
plt.show()
# plt.savefig(dire)#将该名称用于保存图片
relax.write_input(str(gh) + '--' + str(n)) #将生成的VASP输入文件写入存储
dire_1 = str(gh) + '--' + str(n)
diree = os.chdir("./" + dire_1)
fig.savefig('slab.png',
bbox_inches='tight',
transparent=True,
dpi=600,
format='png')
#定义一个更改当前目录的变量
dire2 = './vaspstd_sub'
#确立脚本名称
shutil.copy(r"C:\Users\41958\.spyder-py3\vaspstd_sub", dire2)
#将脚本写入VASP输入文件所在文件夹
with open('surface_area', 'w') as f:
f.write(str(cc))
print(cc)
# os.chdir("../")
#将当前目录改为默认目录
os.chdir(r"D:\Desktop\VASP practical\workdir")
print(os.getcwd())
print('finished')
def absorbed(self,
millerindex_1,
absorbate_1,
absorba,
judge='',
appendage=""):
from pymatgen import Structure, Lattice, MPRester, Molecule
import pymatgen.core.structure
import pymatgen.core.sites
from pymatgen.analysis.adsorption import AdsorbateSiteFinder, reorient_z, plot_slab
from pymatgen.core.surface import generate_all_slabs
from pymatgen.symmetry.analyzer import SpacegroupAnalyzer
from matplotlib import pyplot as plt
from pymatgen.io.cif import CifParser
from pymatgen.io.vasp.inputs import Poscar
from pymatgen.io.vasp.sets import MVLSlabSet
from pymatgen.io.cif import CifWriter
import os
import shutil
ass = self.cif_route
print(ass)
# os.chdir(r"E:\VASP practical\Input")
# print (os.getcwd())
# Note that you must provide your own API Key, which can
# be accessed via the Dashboard at materialsproject.org
#mpr = MPRester()#密钥
struct = CifParser(ass)
structure = struct.get_structures()[0]
print(structure)
os.chdir(r"E:\VASP practical\Input")
print(os.getcwd())
# fcc_ni = Structure.from_spacegroup("Fm-3m", Lattice.cubic(3.5), ["Ni", "Ni"],
# [[0, 0, 0], [0.5, 0.5, 0.5]])
slabs = generate_all_slabs(structure,
max_index=1,
min_slab_size=8.0,
min_vacuum_size=10.0)
millerindex = millerindex_1
struct_111 = [
slab for slab in slabs if slab.miller_index == millerindex_1
][0]
asf_ni_111 = AdsorbateSiteFinder(struct_111)
ads_sites = asf_ni_111.find_adsorption_sites()
# print(ads_sites)
assert len(ads_sites) == 4
fig = plt.figure()
ax = fig.add_subplot(111)
plot_slab(struct_111, ax, adsorption_sites=True)
fig = plt.figure()
ax = fig.add_subplot(111)
adsorbate = Molecule(absorbate_1, absorba)
ads_structs = asf_ni_111.generate_adsorption_structures(
adsorbate, repeat=[1, 1, 1])
A = Poscar(reorient_z(ads_structs[0])) #将切面转换为Poscar
open('POSCAR', 'w').write(str(A))
p = Poscar.from_file('POSCAR')
# w = CifWriter(A.struct)
# w.write_file('mystructure.cif')
path = r'E:\VASP practical\Input\POSCAR' # 文件路径
if os.path.exists(path): # 如果文件存在
# 删除文件，可使用以下两种方法。
os.remove(path)
#os.unlink(path)
else:
print('no such file:%s' % my_file) # 则返回文件不存在
# w = CifWriter(A.struct)
# w.write_file('mystructure.cif')
relax = p.structure #将Poscar 转换为结构信息
custom_settings = {"NPAR": 4} # 用户的INCAR 设置
relaxs = MVLSlabSet(relax, user_incar_settings=custom_settings)
# Vasp输入文件生成器
dire = str(ass) + "---" + str(absorbate_1) + str(millerindex_1)
# print (relax)
relaxs.write_input(dire)
os.chdir("./" + dire)
print(os.getcwd())
#定义一个更改当前目录的变量
dire2 = './vaspstd_sub'
#确立脚本名称
shutil.copy(r"C:\Users\41958\.spyder-py3\vaspstd_sub", dire2)
eb = appendage #添加其他INCAR参数
with open('INCAR', 'r') as f1:
lines = f1.readlines()
with open('INCAR', 'w') as f2:
for line in lines:
if judge in line:
continue
f2.write(line)
with open('INCAR', 'a') as f3:
f3.write(eb)
plot_slab(ads_structs[0], ax, adsorption_sites=False, decay=0.09)
# open('POSCAR001', 'w').write(str(Poscar(reorient_z(ads_structs[0]))))
os.chdir(r"D:\Desktop\VASP practical\workdir")
print(os.getcwd())
print('finished')
# my_lattace = Lattace('mp-698074')#半水石膏
# # my_lattace.phase_out()#生成晶胞优化的输入文件
# go = my_lattace.phase_sol(66,judge='LWAVE', appendage= '\nLWAVE = Ture')
# print('yoo')
def slab(self, miller_index_1, min_slab_size_1=8.0, min_vacuum_size_1=15):
from pymatgen.analysis.adsorption import AdsorbateSiteFinder, plot_slab, reorient_z
from pymatgen.core.surface import Slab, SlabGenerator, generate_all_slabs, Structure, Lattice, ReconstructionGenerator
from pymatgen.symmetry.analyzer import SpacegroupAnalyzer
from pymatgen.core.structure import Structure
from pymatgen.io.cif import CifParser
from matplotlib import pyplot as plt
from pymatgen.io.vasp.inputs import Poscar
from pymatgen.io.vasp.sets import MVLSlabSet
import shutil
import os
##mpr = MPRester()#密钥#密钥
ass = self.cif_route
print(ass)
# os.chdir(r"E:\VASP practical\Input")
# print (os.getcwd())
# Note that you must provide your own API Key, which can
# be accessed via the Dashboard at materialsproject.org
#mpr = MPRester()#密钥
struct = CifParser(ass)
structure = struct.get_structures()[0]
print(structure)
os.chdir(r"E:\VASP practical\Input")
print(os.getcwd())
structure = SpacegroupAnalyzer(
structure).get_conventional_standard_structure()
#空间群分析
need_miller_index = miller_index_1 #通过米勒指数，确定要切的晶面
slab = SlabGenerator(structure, miller_index=need_miller_index, min_slab_size=min_slab_size_1,\
min_vacuum_size=min_vacuum_size_1, center_slab=True)
#晶面生成器参数
for n, slabs in enumerate(slab.get_slabs()):
slabs_bak = slabs.copy() #可能的晶面
slabs.make_supercell(self.supercell)
#晶胞扩充
A = Poscar(slabs) #将切面转换为Poscar
relax = A.structure #将Poscar 转换为结构信息
custom_settings = {"NPAR": 4} # 用户的INCAR 设置
relax = MVLSlabSet(relax, user_incar_settings=custom_settings)
#Vasp输入文件生成器
fig = plt.figure() #绘图--确立画布
ax = fig.add_subplot(111) #绘图--确立位置
plot_slab(slabs, ax, adsorption_sites=False) #绘图
dire = str(ass) + "---" + str(need_miller_index) + '----' + str(n)
#设置一个用作存储输入文件的名称
relax.write_input(dire) #将生成的VASP输入文件写入存储
os.chdir(r"E:\VASP practical\Input")
print(os.getcwd())
plt.savefig('002', format='png') #将该名称用于保存图片
os.chdir("./" + dire)
#定义一个更改当前目录的变量
dire2 = './vaspstd_sub'
#确立脚本名称
shutil.copy(r"C:\Users\41958\.spyder-py3\vaspstd_sub", dire2)
#将脚本写入VASP输入文件所在文件夹
# os.chdir("../")
#将当前目录改为默认目录
os.chdir(r"D:\Desktop\VASP practical\workdir")
print(os.getcwd())
print('finished')