如何使用 PyIron 选择用于声子计算的网格
How to choose the mesh for phonon calculation with PyIron
我想用 phononpy 包计算 pyiron 的声子态密度和能带结构。
我按照教程创建了一份工作:
phono = pr.create_job(pr.job_type.PhonopyJob,"pDOS")
我可以运行这个工作,但是因为网格太密,所以需要很多时间。有没有一种方法可以选择我想要使用的网格?
另外,我想计算给定路径的声子能带结构,用pyiron可以吗?
您可以指定输入:
phono.input
这里可以设置网格为:
phono.input["dos_mesh"]
最佳,
一月
要解决有关能带结构的评论 - 您可以直接使用 phonopy API:
bands = []
q_start = np.array([0.5, 0.5, 0.0])
q_end = np.array([0.0, 0.0, 0.0])
band = []
for i in range(51):
band.append(q_start + (q_end - q_start) / 50 * i)
bands.append(band)
q_start = np.array([0.0, 0.0, 0.0])
q_end = np.array([0.5, 0.0, 0.0])
band = []
for i in range(51):
band.append(q_start + (q_end - q_start) / 50 * i)
bands.append(band)
phon.phonopy.set_band_structure(bands)
phon.phonopy.plot_band_structure().show()
我想用 phononpy 包计算 pyiron 的声子态密度和能带结构。 我按照教程创建了一份工作:
phono = pr.create_job(pr.job_type.PhonopyJob,"pDOS")
我可以运行这个工作,但是因为网格太密,所以需要很多时间。有没有一种方法可以选择我想要使用的网格?
另外,我想计算给定路径的声子能带结构,用pyiron可以吗?
您可以指定输入:
phono.input
这里可以设置网格为:
phono.input["dos_mesh"]
最佳,
一月
要解决有关能带结构的评论 - 您可以直接使用 phonopy API:
bands = []
q_start = np.array([0.5, 0.5, 0.0])
q_end = np.array([0.0, 0.0, 0.0])
band = []
for i in range(51):
band.append(q_start + (q_end - q_start) / 50 * i)
bands.append(band)
q_start = np.array([0.0, 0.0, 0.0])
q_end = np.array([0.5, 0.0, 0.0])
band = []
for i in range(51):
band.append(q_start + (q_end - q_start) / 50 * i)
bands.append(band)
phon.phonopy.set_band_structure(bands)
phon.phonopy.plot_band_structure().show()