The perturbed-chain statistical associating fluid theory (PC-SAFT) equation of state (EoS) was combined with the electrostatic free energy from the mean spherical approximation (MSA) theory, and applied to represent densities of pure imidazolium ionic liquids (ILs) with anion [BF4]−. The PC-SAFT parameters of cations were linearized with their molar mass and obtained by simultaneously fitting the model predictions to experimental densities of some ILs. The PC-SAFT-MSA model provides accurate correlations and predictions of densities comparing with experimental data. Then a classical density functional theory (DFT) was developed based on PC-SAFT-MSA. The DFT model was applied to explore the structure and differential capacitance of the electrical double layer (EDL) in ILs on graphite. The model predicts similar density profiles for both cation and anion on a neutral surface, and layered structure with alternating layers of cations and anions on a charged surface. The charge inversion phenomena were also studied based on ion distributions. We further studied the effects of the alkyl chain length, temperature and non-electrostatic solid-fluid interactions on the differential capacitance of the EDL. The model provides bell-shaped differential capacitance curves. The peak positions of differential capacitance curves shift toward positive potentials as ions size asymmetry increases in agreement with previous experiments and simulations studies. The maximum capacitance decreases with increasing alkyl chain length as well as increasing temperature.