Density functional theory calculations on the reaction mechanisms of the direct arylation of benzene with chlorobenzene mediated by a series of low-valent iron complexes, in which the Fe(II) center is surrounded by different electron-donor ligands (acetate anion (OAc), baphophenanthroline (baph), 1,10-phenanthroline (phen), and redox active ligand amidophenolate (ap)) using density functional theory. Fe(II) models, 1b Fe II (baph), 1p Fe II (phen), 1d Fe II (diimine), 2o Fe II (OAc) 2 , 2po Fe II (OAc) (phen), 2p Fe II (phen) 2 as well as 2a Fe II (ap) 2 were established. According to our calculations, 1b and 2a are promising candidates for the direct arylation transformation. The complexes under different ligands show their unique mechanism characteristics. Furthermore, a correlation has been established among the activation barriers, the energy gaps of frontier orbitals, the distortion energies, as well as the reaction enthalpies. The knowledge obtained herein not only deepens our mechanistic understanding of iron-mediated direct arylation but may also provide guidance for the rational design of catalysts. K E Y W O R D S arylation, C-C coupling, computational chemistry, iron catalysis, ligand effect 1 | INTRODUCTIONLigands play a crucial role in metal-based catalytic reactions. With tunable steric and electronic properties, they are capable of enhancing/reducing the redox potentials of metal centers, increasing/lowering the associated reaction barriers, and so on. [1][2][3][4][5][6][7][8][9][10][11][12][13][14] Modifying the local environment of metal centers in a controlled way enables the modification of local charge and of spin states of the active site; correspondingly, the reactivity could be changed dramatically. [15,16] Quintessentially, Schwarz and coworkers have demonstrated, by using prototypical models, that how and