A numerical study of boson star binaries