We present a new method to accurately describe the ionization fraction and the size distribution of polycyclic aromatic hydrocarbons (PAHs) within astrophysical sources. To this purpose, we have computed the mid-infrared emission spectra of 308 PAH molecules of varying sizes, symmetries, and compactness, generated in a range of radiation fields. We show that the intensity ratio of the solo CH out-of-plane bending mode in PAH cations and anions (referred to as the 11.0 $\mu$m band, falling in the 11.0-11.3 $\mu$m region for cations and anions) to their 3.3 $\mu$m emission, scales with PAH size, similarly to the scaling of the 11.2/3.3 ratio with the number of carbon atoms (N$_{\mathrm{C}}$) for neutral molecules. Among the different PAH emission bands, it is the 3.3 $\mu$m band intensity which has the strongest correlation with N$_{\mathrm{C}}$, and drives the reported PAH intensity ratio correlations with N$_{\mathrm{C}}$ for both neutral and ionized PAHs. The 6.2/7.7 intensity ratio, previously adopted to track PAH size, shows no evident scaling with N$_{\mathrm{C}}$ in our large sample. We define a new diagnostic grid space to probe PAH charge and size, using the (11.2+11.0)/7.7 and (11.2+11.0)/3.3 PAH intensity ratios respectively. We demonstrate the application of the (11.2+11.0)/7.7 - (11.2+11.0)/3.3 diagnostic grid for galaxies M82 and NGC 253, for the planetary nebula NGC 7027, and the reflection nebulae NGC 2023 and NGC 7023. Finally, we provide quantitative relations for PAH size determination depending on the ionization fraction of the PAHs and the radiation field they are exposed to.

We present 5-20 micron spectral maps of the reflection nebula NGC2023 obtained with the Infrared Spectrograph SL and SH modes on board the Spitzer Space Telescope which reveal emission from polycyclic aromatic hydrocarbons (PAHs), C60, and H2 superposed on a dust continuum. We show that several PAH emission bands correlate with each other and exhibit distinct spatial distributions revealing a spatial sequence with distance from the illuminating star. We explore the distinct morphology of the 6.2, 7.7 and 8.6 micron PAH bands and find that at least two spatially distinct components contribute to the 7--9 micron PAH emission in NGC2023. We report that the PAH features behave independently of the underlying plateaus. We present spectra of compact oval PAHs ranging in size from C_66 to C_210, determined computationally using density functional theory, and investigate trends in the band positions and relative intensities as a function of PAH size, charge and geometry. Based on the NASA Ames PAH database, we discuss the 7--9 micron components in terms of band assignments and relative intensities. We assign the plateau emission to very small grains with possible contributions from PAH clusters and identify components in the 7--9 micron emission that likely originates in these structures. Based on the assignments and the observed spatial sequence, we discuss the photochemical evolution of the interstellar PAH family as they are more and more exposed to the radiation field of the central star in the evaporative flows associated with the PDRs in NGC2023.