We have theoretically studied the possibility to control the equilibrium solubility of dopants in semiconductor alloys, by strategic tuning of the alloy concentration. From the modeled cases of C 0 in Si x Ge 1ÿx , Zn ÿ and Cd ÿ in Ga x In 1ÿx P it is seen that under certain conditions the dopant solubility can be orders of magnitude higher in an alloy or multilayer than in either of the elements of the alloy. This is found to be due to the solubility's strong dependence on the lattice constant for size mismatched dopants. The equilibrium doping concentration in alloys or multilayers could therefore be increased significantly. More specifically, Zn ÿ in a Ga x In 1ÿx P multilayer is found to have a maximum solubility for x 0:9, which is 5 orders of magnitude larger than that of pure InP. [8][9][10]. From a theoretical point of view the raised nonequilibrium concentrations have been successfully explained by the increased range of the chemical potential of the dopant during epitaxial growth [11][12][13].In this Letter we consider ways to increase the solid solubility to achieve higher doping concentrations under equilibrium conditions. We study the dependence of the equilibrium solubility of substitutional defects, including dopants, in compound semiconductors A x B 1ÿx . We demonstrate that the equilibrium solubility of dopants need not be monotonic in x and lie between those of the pure materials A and B. Instead, we will show that careful tuning of x can result in equilibrium solubilities many orders of magnitude higher in the alloy than in either of the pure materials.It is well known that high doping can alter the lattice constant of a material [14]. Here we investigate the possibility of using changes in the lattice constant to control the doping concentration. Initial indications of this have been noted in previous experimental work [15], in theoretical studies of dopants in silicon under applied biaxial strain [16] and for vacancies and self-interstitials using a phenomenological method [17]. However, the full implications of this possibility have not so far been realized.Our hypothesis builds on the way in which the solubility depends on the difference between the dopant's and the substituted host atom's equilibrium bond lengths. The solubility of a dopant atom of the same size and equilibrium bonding distance as the host atoms will have little dependence on the strain of the material since the difference in bond energy will be constant. If the equilibrium bonding distances differ, however, the solubility of the dopant may have a strong dependence on the lattice constant, such that the equilibrium solubility of smaller (larger) dopants increases under compressive (expansive) strain. Our results show for the first time that the equilibrium concentration of a defect can sometimes be orders of magnitude higher in an alloy or multilayer than in either of the constituents. This leads to the possibility of increasing the equilibrium solubilities of dopants through strategically composing multilayers or alloys...