Ab initio modeling of boron clustering in silicon