The adsorption isotherms of CO2 in several porous aromatic frameworks (PAFs) have been simulated with Grand Canonical Monte Carlo technique, to support the synthesis of new materials for efficient carbon dioxide capture and storage. The simulations covered the 0-60 bar pressure range and were repeated at 273, 298, and 323 K. The force field employed in the simulations was optimized to fit the correct behavior of the free gas and to reproduce the CO2-phenyl interactions computed at high quantum mechanical level. PAFs are based on the diamond structure, with polyaromatic chains inserted in C-C bonds. We examined four PAF-30n (n being the number of phenyl rings in the aromatic linkers), finding that PAF-302 is overall the best performing, although PAF-301 provides higher adsorbed densities at very low pressure. The CO2 adsorption then was simulated in a number of modified PAF-302, with different functional groups (aminomethane, toluene, pyridine, and imidazole) attached to the phenyl chains; different degrees of substitution (25%, 50%, and 100% derivatized rings) were considered. The effects of functionalization and the dependence on the substitution degree are carefully discussed, to determine the most promising materials at low, intermediate, and high pressures.