Virtual quaternary alloy compositions were created of the form NbaXxYyZz where X, Y, Z represents any of the 14 alloying elements B, C, Co, Cr, Hf, Mo, N, Re, Si, Ta, Ti, V, W, Zr, and x, y, and z represent their corresponding weight percent value and 𝑎 = [100 -(𝑥 + 𝑦 + 𝑧)] represents the weight percent of Nb. These 14 alloying elements were selected based on domain knowledge of the promising candidates to provide the solid solution, precipitation, and dispersion strengthening. The concentration granularity of the alloying elements was 1 weight percent. The concentration of B, C, and N varied from 1 weight percent to 5 weight percent, and for the rest of the elements, concentration varied from 1 weight percent to 20 weight percent. This combination yielded 1,666,625 (approximately 1.6 million) unique quaternary alloy compositions. With a similar approach, quinary alloy compositions of the form NbaWwXxYyZz have been generated where W, X, Y, and Z represent any of the 14 alloying elements and w, x, y, and z represent their corresponding weight percent value and 𝑎 = [100 -(𝑤 + 𝑥 + 𝑦 + 𝑧)] represents the weight percent of Nb. A total of 74,277,500 (approximately 74 million) quinary alloy compositions were created. The next step was to screen out optimal alloy composition from this vast virtual candidate search space.