Hydrogenolysis is an important approach for depolymerization of lignin, which provides attractive new sustainable platforms of fuels, chemicals, and materials. The theory of lignin hydrogenolysis is, however, still unsound, which limits the development of this approach and causes inconsistencies among experimental studies. In this paper, density functional theory is employed to investigate the initial hydrogenolytic cleavages of recognized five different types of interaromatic unit linkages of lignin, assuming the presence of hydrogen free radicals. The relative free energies of reactant complexes, reaction free energy changes, and rate constants for candidate reactions are calculated comprehensively at 298-538 K On the basis of the results of calculation and a rapid equilibrium hypothesis, the major reaction channel is decided for each linkage, and its kinetics is assessed. It is concluded that the hydrogenolysis occurs at beta-O-4 ether, diphenyl ether 4-O-5', and beta-1' diphenylmethane linkages instantaneously if these are accessible to hydrogen free radicals, while beta-5 phenylcoumaran and beta-beta' pinoresinol linkages are virtually inert to hydrogenolysis.