Even though tremendous achievement has been made experimentally in the performance of Li-S battery, theoretical studies in this area are lagging behind due to the complexity of the Li-S systems and effects of solvent. For this purpose, we have developed a new methodology for investigating the 2D hexaaminobenzene-based coordination polymers (2D-HAB-CP) as cathode candidate materials for Li-S battery via density functional theory calculation in combination with in-house developed charge polarized solvent model and genetic algorithm structure global search code. With high ratios of transition metal atoms and 2-coordinated nitrogen atoms, excellent electric conductivity, and structural porosity, the 2D-HAB-CP is able to address all of the three main challenges facing Li-S battery: confining the lithium polysulfides from dissolution, facilitating the electron conductivity and buffering the volumetric expansion during the lithiation process. In addition, the theoretical energy density of this system is as high as 1395 Wh/Kg. These results demonstrated that the 2D-HAB-CP is a promising cathode material for Li-S battery. Our proposed computational framework not only opens a new avenue in understanding the key role played by solution and liquid electrolytes in Li-S battery, but also can be generally applied to other processes with liquids involved.