where Ai, εi, and Γi are the integrated intensity, peak position, and half width at half maximum (HWHM) of the Lorentzian function of peak i, respectively. s represents the full width at half maximum (FWHM) of the Gaussian function, which is the instrumental function. We set M = 1 for NPC and NPG synthesized with θ-Al2O3, and zeolite-templated carbon (ZTC), M = 2 for stacked NPC and NPG, M = 3 for NPG synthesized with γ-Al2O3, and M = 7 for NPC synthesized with γ-Al2O3 to improve the numerical stability of the deconvolution procedure. We fixed s at 13.7 cm -1 , according to the energy resolution of the spectrometer. The other parameters are fully optimized using a leastsquares fitting procedure. Each peak intensity Ii (εi) was evaluated as Ii (εi) = AiVi(εi).