Thermodynamically favorable 4B10BV (BV = Sb, Bi) clustering, reducing the internal strain energy, is presented in highly mismatched InAs-rich BxIn1-xBVyAs1-y alloys. The cluster formation decreases the sum of the internal energies of the constituent compounds and internal strain energy. The cohesive energies, enthalpies of formation, bond stretching, and bond bending elastic constants of zinc blende BBi and InBi are calculated. The clustering conditions are estimated up to 800 °C within the compositions 5 × 10−5 ≤ x ≤ 0.01 and x ≤ y ≤ 100x. The internal strain energies of BxIn1-xSbyAs1-y and BxIn1-xBiyAs1-y (y ≥ 2.5x) in which 99.99% boron atoms are in the clusters are more than three and seven times, respectively, less than in alloys without clusters. The formation conditions for such weakly strained highly mismatched alloys have been obtained. Almost all boron atoms should be in 4B10Sb and 4B10Bi clusters at 800 °C in BxIn1-xSbyAs1-y with x ≥ 4.2 × 10−4 and y ≥ 10x and in BxIn1-xBiyAs1-y with x ≥ 6.6 × 10−5 and y ≥ 10x.