We have performed a density functional theory investigation of hydroformylation of ethylene for monosubstituted rhodium-carbonyl catalysts, HRh(CO)3L, where the modifying ligand, L, is a phosphite (L = P(OMe)3, P(OPh)3, or P(OCH2CF3)3), a phosphine (L = PMe3, PEt3, PiPr3, or PPh3), or a N-heterocyclic carbene (NHC) based on the tetrahydropyrimidine, imidazol, or tetrazol ring, respectively. The study follows the Heck and Breslow mechanism. Excellent correspondence between our calculations and existing experimental information is found, and the present results constitute the first example of a realistic quantum chemical description of the catalytic cycle of hydroformylation using ligand-modified rhodium carbonyl catalysts. This description explains the mechanistic and kinetic basis of the contemporary understanding of this class of reaction and offers unprecedented insight into the electronic and steric factors governing catalytic activity. The insight has been turned into structure-activity relationships and used as guidelines when also subjecting to calculation phosphite and NHC complexes that have yet to be reported experimentally. The latter calculations illustrate that it is possible to increase the electron-withdrawing capacity of both phosphite and NHC ligands compared to contemporary ligands through directed substitution. Rhodium complexes of such very electron-withdrawing ligands are predicted to be more active than contemporary catalysts for hydroformylation.