def __init__(self, probe_symbol, surface_symbol, ): self.__probe_symbol = probe_symbol surface = BodyCenteredCubic111(surface_symbol) l = vectorLength(surface.cell()[1])/3 x, y = self.probePlacement(l) positions = [(x,y,2.210), (x, y + l, 0.92), (x + l*np.cos(np.pi/6), y-l*np.sin(np.pi/6), 0.92), (x - l*np.cos(np.pi/6),y-l*np.sin(np.pi/6), 0.92), ] added_atoms = Atoms(probe_symbol*len(positions), positions) SurfacePointDefect.__init__( self, surface=surface, added_atoms=added_atoms, )