In [19]:
import sys
from IPython.display import Image
from siman import header
from siman.calc_manage import smart_structure_read
from siman.geo import create_supercell, create_surface2, supercell
%matplotlib inline
In [4]:
st = smart_structure_read('Cu/POSCARCU.vasp') # read required structure
In [7]:
Image(filename='figs/Thompson-tetrahedron-notation-for-FCC-slip-systems.png')
Out[7]:
In [10]:
# create supercell using chosen directions, the *mul* allows to choose one half of the third vector
sc = create_supercell(st, [ [-1,-1,-1], [0,1,-1], [2,-1,-1]], mul = (1,1,0.5))
In [18]:
# here we choose [100] normal in supercell, which is equivalent to [111]cub
# combinations of *min_slab_size* and *cut_thickness* (small cut of slab from one side) allows create symmetrical slab
st_suf = create_surface2(sc, [1, 0, 0], min_vacuum_size = 10,
min_slab_size = 16, cut_thickness = 3, oxidation = {'Cu':'Cu0+' }, return_one = 1, surface_i = 0)
In [20]:
st_sufsc112 = supercell(st_suf, [10,10,32]) # make 2x2 slab
In [21]:
st_sufsc112.write_poscar() # save file as POSCAR for VASP
Out[21]: