The efficient fixation and utilization of CO 2 has been consistently pursued by chemists for decades. Although Cu-based catalysts, e.g., Cu/ZnO/Al 2 O 3 , have been widely used in industry for methanol synthesis from CO 2 hydrogenation (CO 2 + 3H 2 !H 3 COH + H 2 O), many issues on the mechanism and the kinetics remain largely uncertain. For example, the surface site for CO 2 activation and the synergetic effect between Cu and oxide have been hotly debated in literature. In the past few years, theoretical modeling on pure Cu surfaces and Cu/oxide interfaces has been utilized to provide insight into these important questions. Here we will review the recent theoretical advances on simulating this complex heterogeneous catalytic process with first principles density functional theory (DFT) calculations and kinetics modeling. The theoretical results on the mechanism and the kinetics are compared and summarized.