如何构建Si和C的过渡层SixC(1-x)

jasondeng

各位大佬好，我想请教一下我想在Si和C之间构建一个过渡层SixC(1-x)（x=0时刚好完全是Si,x=1时刚好完全是C），类似于下面的结构，然后计算他的声子态密度，请问我应该如何构建啊，是利用VESTA还是用python自己设置和替换原子位置啊，不知道vaspkit有没有这种功能。