Considering the great potential of deep eutectic solvent (DES) for CO2 capture, this work studies the CO2 solubility in DESs by combining experimental measurement and molecular dynamics (MD) simulation. First, four phosphonium-based DESs have been prepared in laboratory, involving two types of hydrogen bond acceptors (HBAs), namely tetrabutylphosphonium bromide (TBPB) and allyltriphenylphosphonium bromide (ATPPB), and two types of hydrogen bond donors (HBDs), namely phenol (PhOH) and diethylene glycol (DEG). The CO2 solubility in the obtained DESs has been measured under 313.15 K -333.15 K and pressure below 2000 kPa, and compared with that of previously reported DESs and ionic liquids. Second, MD simulations have been performed to study the microscopic behaviors of the involved DESs and mixtures. Through the analyses of radial distribution functions (RDFs), spatial distribution functions (SDFs) and intermolecular interaction energies, the eutectic formation and CO2 absorption mechanisms, as well as the effect of HBA/HBD type and molar ratio are interpreted.