\section{Introduction} \label{sec:intro}
Galaxies with active star formation are dominated by strong emission features at 3.3, 6.2, 7.7, 8.6, 11.3, 12.7, 16.4, and 17\,\micron\ complex generally attributed to Polycyclic Aromatic Hydrocarbons \citep[PAHs;][]{Leger1984, Allamandola1985, Smith2004, Werner2004}. These molecules typically consist of 20--1000 carbon atoms \citep{Leger1984, Tielens2008}, and have been found to be an abundant and energetically important component of interstellar dust, contributing $\sim5$--20\% of the total infrared emission from a galaxy \citep{Smith2007}. In addition, they are extremely effective in reprocessing radiative energy and coupling photon energy to the gas through photoelectric heating \citep{DHendecourt1987}. PAH grains of different sizes and ionization states vary widely in the efficiency with which they emit at different bands \citep{Draine2007}. However, the underlying mechanism that governs such variation had not yet been well understood.

Spanning the wavelength range of 2.4--45\,\micron\, the Short Wavelength Spectrometer (SWS) on board the $\textit{Infrared Space Observatory}$ \citep[ISO: ][]{Kessler1996} offered for the first time the sensitivity to exploit the diagnostic power of many of the most luminous PAH bands \citep{Spoon2000, Peeters2002}. The advent of the sensitive Spitzer/IRS instrument \citep{Houck2004} permitted detailed explorations of PAH emission at wavelengths of 5–-20\,\micron\, and provided crucial insights regarding the nature of PAH populations (i.e., size, charge state) within differing interstellar environments and on galaxy-wide scales \citep{Smith2007, Galliano2008, Gordon2008, Hunt2010, Sandstrom2012}. The metal content and radiation field of the ISM have been suggested as the primary drivers of variations in PAH emission. PAH abundance was found to decline rapidly as metallicity drops, vanishing entirely below 12+log(O/H)$\sim$8.1 \citep{Engelbracht2005, Draine2007, Smith2007, Sandstrom2012}. Also, \citet{Smith2007} found a very low PAH 7.7/11.3 ratio in galaxies that host an AGN \citep[see also][]{ODowd2009, Diamond-Stanic2010}. These deficits have been interpreted as preferential destruction of small PAHs by strong radiation fields, leading to an enhanced decrease of the shorter-wavelength band relative to the longer-wavelength band. In those studies, however, 3.3\,\micron\ PAH emission was left unexplored, due to the limited spectral coverage of \textit{Spitzer}.

Putting the 3.3\,\micron\ PAH in context with all other PAH bands at longer wavelengths is crucial. The 3.3\,\micron\ C-H stretching feature is known to originate from the smallest members of the emitting PAH population \citep{Desert1990, Schutte1993, Draine2007}, making it the most sensitive probe to the size distribution at the grain/molecule interface given their larger surface-to-volume ratios.  Very recently, \citet{Maragkoudakis2020} utilized the NASA Ames PAH IR Spectroscopic Database \citep[PAHdb;][]{Bauschlicher2010} to compare 3.3\,\micron\ PAH to other PAH bands and concluded that the 3.3\,\micron\ band is most sensitive to the sizes of the smallest PAH molecules. Many theoretical and laboratory studies have shown that the ratio of 3.3/11.3 should be a powerful diagnostic to the size distribution as both bands originated from neutral PAHs \citep{Allamandola1999, Croiset2016, Maragkoudakis2018, Maragkoudakis2020}.

The 3.3\,\micron\ PAH emission has been observed by ground-based facilities in various galaxies, including starbursting and infrared luminous galaxies, as well as obscured AGNs associated with star-forming activity \citep{Moorwood1986, Imanishi2000, Imanishi2002, Rodriguez_Ardila2003, Spoon2003, Imanishi2004, Imanishi2006, Imanishi2006b, Risaliti2010}. Taking advantage of the redshifted spectra, \citet{Sajina2009} even observed the 3.3\,\micron\ feature in galaxies at z\,$\sim\,$2 with \textit{Spitzer}/IRS. While \textit{ISO}\ also provided access to the 3.3\,\micron\ PAH emission \citep{Spoon2000}, the interpretation of these spectral features was limited by sensitivity to sampling only bright, nearby sources \citep{Sturm2000, Lu2003, VanDiedenhoven2004}, and the 3.3\,\micron\ PAH was typically only marginally detected.

The 3.3\,\micron\ PAH observed by the Infrared Camera \citep[IRC:][2.5--5\,\micron]{Onaka2007} on board the \textit{AKARI}\ Space Telescope \citep{Murakami2007} has been applied to a suite of diagnostics to identify obscured AGNs in (U)LIRGs \citep{Imanishi2008, Imanishi2010, Lee2012, Ichikawa2014}. These authors found AGN-dominated sources have relatively low 3.3\,\micron\ equivalent width (EQW), attributing it to either the destruction of small PAHs by radiation in the AGN environment or the enhanced integrated infrared luminosity (\mbox{8--1000\,\micron}; \,L\ensuremath{_\mathrm{IR}}) from the AGN activity. Recently, a starburst/AGN diagnostic technique using the 3.3\,\micron\ PAH and the slope of the IRC continuum has also been reported \citep{Inami2018}. The fact that 3.3\,\micron\ PAH is bright in star-forming galaxies while suppressed in star-forming sources hosting AGNs draws into question whether 3.3\,\micron\ PAH can serve as a reliable star formation rate indicator \citep{Yano2016, Inami2018}. Several studies have shown a correlation between L\ensuremath{_\mathrm{PAH\,3.3}}\ and \,L\ensuremath{_\mathrm{IR}}\ in local LIRGs, but also find a 3.3\,\micron\ PAH deficit at the higher luminosity of ULIRGs \citep{Imanishi2010, Kim2012, Y2013,  Murata2017}. The measurement of the 3.3\,\micron\ PAH in these studies, however, dealt only with the attenuation from a nearby 3.05\,\micron\ water ice absorption band, but not the steeply rising continuum opacity at short mid-infrared (MIR) wavelengths near 3\,\micron, which can be substantially larger.

To investigate the role of attenuation and environment on the full MIR spectrum, and model all the prominent vibrational emission bands of PAH molecules in a widely selected sample of galaxies, we present the \textit{AKARI}-\textit{Spitzer}\ Extragalactic Spectral Survey (ASESS). ASESS is an exhaustive cross-archival match between the \textit{AKARI}/IRC dataset and the \textit{Spitzer}/IRS Infrared Database of Extragalactic Observables from \textit{Spitzer}\ (IDEOS) catalog \citep{Hernan-Caballero2016}. With a comprehensive spectral range of 2.5--38 $\mu m$, this survey gains access to $\textit{all}$ PAH bands in the MIR and to model them consistently using a systematic approach. Using a subset of the full \textit{AKARI}-\textit{Spitzer}\ cross-archival sample, we focus on a collection of star-forming galaxies that exhibit bright PAH emission over a broad range of total infrared luminosities (10$^{9.5-12.5}$ L$_\odot$). We present our sample selection in \S~\ref{sec:sample_selection}, followed by a detailed discussion of the spectral decomposition model in \S~\ref{sec:pahfit_decomposition}. We present our results in \S~\ref{sec:result} and discuss their interpretation in \S~\ref{sec:discussion}, summarizing these findings in \S~\ref{sec:summary}. Throughout this paper, a concordance cosmology of H$_\mathrm{0}$=75 km s$^{-1}$ Mpc$^{-1}$, $\Omega_\mathrm{M}$=0.3, and $\Omega_\mathrm{\Lambda}$=0.7 is adopted.

\begin{figure*}
    \centering
        \includegraphics[width=0.495\textwidth]{sample_slc}
        \includegraphics[width=0.495\textwidth]{stern_wedge_sample_slc}
    \caption{Selection process of our bright-PAH sample. (a) Following the galaxy classification proposed by \citet{Spoon2007}, we present our parent sample based on the PAH 6.2 $\mathrm{\mu} m$ EQW and the 9.7 $\mathrm{\mu} m$ silicate strength. We first place a silicate cut (-1.4\,$<$ S$_\mathrm{Si}<$\,0) to include galaxies only subjected to moderate obscuration. (b) To exclude a number of low-obscuration AGNs, we further leave out galaxies lie within the ``Stern Wedge" (gray area). Consequently, with a spectral quality cut (3.3\,\micron\ PAH $>$\,3$\sigma$ detection), our bright-PAH sample consists of 113 galaxies satisfying both criteria as shown in red in (a) (see \S~\ref{subsec:bright_PAH} for detail). The two galaxies in 1A are II~Zw~40 and IRAS 01004-2237.}

    \label{fig:sample_selection}

\end{figure*}

\section{Sample Selection, Observations, and Data reduction}
\label{sec:sample_selection}
With $\sim$100$\times$ better sensitivity than \textit{ISO}/SWS, the combined power of the \textit{Spitzer}\ with \textit{AKARI}\ allows us to extend deep MIR spectral coverage to more distant and less luminous sources. The ASESS sample is drawn from an exhaustive cross-archival comparison between the Infrared Database of Extragalactic Observables from \textit{Spitzer}\ (IDEOS) catalog \citep{Hernan-Caballero2016} and the public data provided by \textit{AKARI}\ Infrared Camera (IRC; \citealt{Onaka2007}) pointed spectroscopic observations compiled in the DARTS database Provided by ISAS, JAXA\footnote{https://darts.isas.jaxa.jp/astro/akari/pointing/}.

The IRC offers two observing modes in the NIR: NG (high-resolution grism) and NP (low-resolution prism) mode. The grism mode covers a spectral range from 2.5--5$\mu m$ with R$\sim$120 at 3.6\,\micron\ for point sources \citep{Onaka2007}, whereas the prism mode covers a spectral range from 1.5--5.5\,\micron\ with low-resolution of R$\sim$20 at 3.6\,\micron\ \citep{Ohyama2007}. In this cross-archival comparison, only the spectroscopic observations of 2220 galaxies with higher resolution NG/grism mode are considered. Also considered are 3361 IDEOS galaxies with \textit{Spitzer}/IRS low resolution (R\,$\simeq$\,60--125) staring mode observations covering the spectral range from 5.2--38\,\micron\ \citep{Hernan-Caballero2016}.

In this study, we focus primarily on star-forming galaxies exhibiting strong PAH emission with modest obscuration. In this section, we describe the selection of our samples, followed by a detailed discussion of data reduction, \textit{AKARI}+\textit{Spitzer}\ spectral combination, and end with our selected subset that consists of exclusively star-forming galaxies.

\subsection{Parent Sample}
The parent sample is selected by matching sources from the IDEOS catalog to the IRC catalog. For each IRC spectrum, we create spatial profiles by averaging the spectral trace along the wavelength direction and fit up to two Gaussians. For the systems that consist of multiple nuclei within the IRC window, i.e., galaxy mergers, defined by an overlap of the Gaussian profiles within 2$\sigma$, are excluded. We further remove extended sources and spectra suffering from hot/cold pixels that are unable to recover by our profile fitting method (see \S~\ref{subsec:AKARI} for detail). Finally, we only include those spectra that have full coverage of 2.5--38\,\micron\ for spectral decomposition, since a small number of IRS spectra have missing segments. These criteria result in a sample of 378 galaxies in our parent sample.

\subsection{IDEOS}
IDEOS provides the principle mid-infrared diagnostics for more than 3000 galaxies/galactic nuclei from spectroscopic measurements obtained by the IRS on the \textit{Spitzer}\ Space Telescope. The spectra have been drawn and vetted from the Combined Atlas of Spitzer-IRS Sources \citep[CASSIS;][]{Lebouteiller2011} a repository of 13500 homogeneously extracted, publication-quality, low-resolution spectra of all \textit{Spitzer}-IRS observations performed in staring-mode (as opposed to mapping mode). For more detail about the observation and data reduction of IDEOS, we refer the reader to a paper by \citet{Hernan-Caballero2016}.

\subsection{AKARI/IRC}
\label{subsec:AKARI}
Our 2.5--5\,\micron\ spectroscopic sample is compiled from the archival data of the \textit{AKARI}\ \citep{Murakami2007} mission from the Infrared Camera (IRC), obtained by using the Astronomical Observing Template (AOT) IRC04 and IRCZ4 with the configuration b;Np, where b refers to the NG grism mode and Np refers to point source spectroscopy. While this is a slit-less mode, an aperture of $1\arcmin \times 1\arcmin$ is designed to block lights from nearby sources. To optimize the process of cosmic-ray removal, each target was observed by one to ten pointings, depending on the brightness of the source, and each pointing consists of eight or nine sub-frames. For each pointing, the total on-source exposure time was $\sim$ 6 minutes \citep{Ohyama2007}.

Our parent sample includes three phases of observations during the \textit{AKARI}\ mission -- phase 1 \& 2 were observations taken during the cryogenic phase, while phase 3 was post-cryogenic. The number of galaxies in phase 1 \& 2 and 3 are 50 and 328, respectively. Data reduction was carried out using \textit{IDL} packages prepared by the \textit{AKARI}\ team: ``IRC Spectroscopy Toolkit Version 20080528'' for phase 1 \& 2 data and ``IRC Spectroscopy Toolkit for Phase 3 data Version 20170225'' for phase 3 data. Both packages involve dark subtraction, linearity correction, flat-fielding, sub-frame combination, source detection, background/source subtraction, and wavelength and flux calibration. Furthermore, we applied a correction for second-order contamination as detailed in \citet{Baba2016, Baba2019}. We derived the spatial profile of each source from the 2D spectrum, using the 2$\sigma$ range for the spatial extraction width, which varies from 4-12 pixels (5$\farcs$85--17$\farcs$52). Extended galaxies that have a broad spatial profile with 2$\sigma$ range $>$12 pixels (17$\farcs$52) are excluded from our sample.

Several additional issues emerged in phase 3 after cryogen exhaustion: 1) the detector sensitivity was reduced by $\sim 70 \%$ compared to phase 2\footnote[2]{\textit{AKARI}\ observer's manual version 1.2}, and 2) the occurrence of \emph{hot} pixels increased to $\sim 5 \%$ of the total. These effects make 2D spectral images contaminated. To correct these contaminations, we follow the prescription from \citep{Horne1986} to construct a normalized spatial profile at each wavelength, and average these over wavelength to construct a single characteristic profile. The defect pixels, which typically shown as NaNs in the 2D spectral image, are then replaced by values scaled from the characteristic profile at the corresponding wavelength. To reach a higher S/N, we combine available spectra from different pointings (i.e., with similar primary ID but different sub-IDs) by calculating the weighted average.

\subsection{Spitzer-AKARI Spectra Stitching}
The $\sim$0.2\,\micron\ gap in between the long-wavelength end of \textit{AKARI}/IRC and the start of \textit{Spitzer}/IRS complicates the process of combining the two spectral segments. For all 378 \textit{AKARI}-\textit{Spitzer}\ combined spectra in our parent sample, we scaled the flux of each segment to match their synthetic photometry to the \textit{WISE} photometry obtained from the profile-fitted magnitude in the ALLWISE Source Catalog \citep{Cutri2013}. The values were converted into fluxes in Jy based on the WISE calibration provided by \citet{Jarrett2011}. Only the W1 band is used to scale the IRC spectrum because the band transmission profile of W2 falls beyond the IRC spectral range. For scaling the IRS spectrum, we first follow the method of \citet{Wright2010} to correct the W4 fluxes, which appear to be $\sim$8--10\% too bright in the red sources, and further calculate the scaling factor using the average of W3 and W4 bands. For three sources lacking WISE photometry (IRAS 15250+3609, WKK 2031, and IC 4518 A), we used \textit{IRAC} photometry in a similar manner from the Spitzer Enhanced Imaging Products (SEIP) database to scale our spectra. For two additional sources (UGC 02369 and NGC 5990), where no photometry is provided in either database, we performed a linear extrapolation from feature-free regions on each side of the gap, scaling the IRC spectrum such that the two lines meet at gap center (5.1\,\micron). We find the scaling factors of both the \textit{AKARI}\ and \textit{Spitzer}\ spectra lie mostly within 20\%, similar to the values reported by \citet{Baba2018}.

\subsection{Star-Forming Subset: Bright-PAH Sample} \label{subsec:bright_PAH}
Our parent sample consists of sources including (i) star-forming galaxies with prominent PAH emissions, (ii) heavily obscured AGNs, with strong 9.7 and 18\,\micron\ silicate absorptions and PAH-diluting warm continuum, and (iii) low obscured AGN and Quasi-Stellar Objects (QSO) with smooth continuum, mainly dominated by dust re-emission with weaker PAH emissions. For the purpose of this study, we further narrow down our sample to galaxies exhibiting bright PAH emissions and dub them the ``bright-PAH sample". \citet{Spoon2007} provided a galaxy classification scheme based on the equivalent width (EQW) of the 6.2\,\micron\ PAH and the silicate strength, S$_\mathrm{Si}$, defined by S$_\mathrm{Si}$ = ln(f$\mathrm{_{9.7 \mu m}/C_{9.7 \mu m})}$, where $f_\mathrm{9.7 \mu m}$ is the central flux intensity of the silicate absorption feature and C$_\mathrm{9.7 \mu m}$ is the estimated continuum flux intensity adopted for the local MIR spectrum. We follow this methodology and present our sample in Fig.~\ref{fig:sample_selection}(a) using the 6.2\,\micron\ PAH EQW and the silicate strength both obtained from the IDEOS catalog. The average spectral template of each class will be presented in \S~\ref{sec:template}. Galaxies exhibiting strong PAH features are typically subjected to less silicate absorption. We further select only those galaxies with -1.4 $<$ S$_\mathrm{Si}$ $<$ 0 (circles), thus avoiding both deep obscuration (triangles) and silicate \emph{emission} (squares) contributed by inner AGN regions \citep[e.g.,][]{Fritz2006, Nenkova2008}. Galaxies satisfying the above criterion lie in between the two dashed lines in Fig.~\ref{fig:sample_selection}(a), but with this simple silicate cut, a modest number of lower-obscuration AGN sources remain.

To exclude these sources, we use the IRAC photometric AGN classification scheme of \citet{Stern2005}. In general, AGNs are redder in the [3.6]-[4.5] color and restricted to a narrower distribution of [5.8]-[8.0] color since those sources tend to have less prominent PAH emission diluted by the strong AGN continuum. We use synthetic photometry from IRAC response curves, excluding from our bright-PAH sample those IRAC-selected AGNs that lie inside the ``Stern Wedge" (shaded region in Fig.~\ref{fig:sample_selection}(b)), which are added with cross signs. The symbols used here are similar to Fig.~\ref{fig:sample_selection}(a) and are color-coded by strength of the 6.2\,\micron\ PAH EQW. Galaxies exhibiting strong PAH emission tend to cluster at the bottom right corner as opposed to another galaxy cluster that lies inside the AGN region. Finally, we limit our bright-PAH sample to have robust ($>$\,3$\sigma$) detection of the 3.3\,\micron\ PAH band. Our bright-PAH sample of 113 galaxies that meet all the imposed selection criteria are colored in red in Fig.~\ref{fig:sample_selection}(a). The two bright-PAH galaxies that lie within 1A are II~Zw~40 (See \S \ref{subsec:IIZw40} for detail) and IRAS~01004-2237, which hosts a moderately obscured AGN \citep{Imanishi2007} leading to a weak 6.2\,\micron\ PAH emission. Even with this selection method, substantially buried AGNs which do not significantly contribute to the MIR continuum may still contribute to the total infrared luminosity, especially for galaxies with \,L\ensuremath{_\mathrm{IR}}$>$10$^{12}$\,L\ensuremath{_\mathrm{\odot}}\ \citep[][see also \S~\ref{sec:all_the_PAHs}]{Imanishi2008, Nardini2010, Petric2011, Alonso-Herrero2012, Imanishi2018, Imanishi2019}.

To illustrate how infrared luminosity may affect the overall shape of the spectrum, we divide bright-PAH galaxies into three categories according to \,L\ensuremath{_\mathrm{IR}}. Fig. \ref{fig:brightPAH_specs} shows these averaged spectra normalized at 5\,\micron, with primary PAH emissions indicated. The \textit{AKARI}/IRC spectra of galaxies with log\,\,L\ensuremath{_\mathrm{IR}}\,$<12$ show dominant stellar components, while galaxies with log\,\,L\ensuremath{_\mathrm{IR}}\,$\geq 12$ have a steeper continuum slope towards red. Table~\ref{tab:bright_PAH_sample} summarizes basic information of galaxies in the bright-PAH sample, with Table~\ref{tab:bright_PAH_id} summarizes the observational details (See Appendix). The redshifts of our subset span a range 0 $\leq$ z $\leq$ 0.2 with a median of z=0.03.

\input{tbl1_brightPAH_basic_info}

\begin{figure}
 	\includegraphics[width=1.0\columnwidth]{brightPAH_specs}
    \caption{The averaged spectra of our bright-PAH sample binned by log\,\,L\ensuremath{_\mathrm{IR}}. Spectra are normalized at 5\,\micron\ and the width of the stripe indicates the FWHM of the main PAH features centered at 3.3, 6.2, 7.7, 8.6, 11.3, and 17\,\micron. The numbers of galaxies used in the blue, yellow, and red spectra are 16, 73, and 23, respectively.}
    \label{fig:brightPAH_specs}
\end{figure}

\section{Decomposition of galaxy spectra from 2.5$-$38 $\mu m$}
\label{sec:pahfit_decomposition}
The infrared wavelengths from 2.5--5\,\micron\ host a rich inventory of dust and line emission features, together with continuum and ice absorption. To assess how these features vary with the ISM environment, we use \textsc{Pahfit}\ to decompose the spectrum and to extract their integrated strength. \textsc{Pahfit}\ is a robust decomposition model first employed for estimating the strengths of weak, blended features of the low-resolution \textit{Spitzer}/IRS spectra in the MIR regime \citep[5.2--38 $\mu m$;][]{Smith2007}. IRS low-resolution spectra have a 5.2\,\micron\ short wavelength cutoff, and therefore the 3.3\,\micron\ PAH band was not included in \textsc{Pahfit}. Our main objective in extending \textsc{Pahfit}\ model coverage down to 2.5\,\micron\ is to simultaneously model \emph{all} the PAH emission features with their underlying continua and absorption features. In this section, we first discuss the newly added emission components, followed by the new absorption components together with the adopted geometries. Fig.~\ref{fig:pahfit_example} shows the decomposition result of a typical star-forming galaxy.

\begin{figure*}
 	\includegraphics[width=1.0\textwidth]{PAHFIT_example_1C_snr_gt60}
    \caption{Detailed decomposition of the normalized type 1C spectrum ranging from 2.5--6\,\micron. The newly added features in the enhanced \textsc{Pahfit}\ are labelled, including dust emission features (blue) of the 3.3\,\micron\ PAH, 3.4\,\micron\ aliphatic, and 3.47\,\micron\ plateau, line emission features (purple) of the 3.74\,\micron\ Pf$\gamma$, 4.05\,\micron\ Br$\alpha$, and 4.65\,\micron\ Pf$\beta$ (Pf$\beta$ are not obvious in this case). The dashed line indicates the attenuation caused by the fully mixed extinction with absorption features of H$_\mathrm{2}$O ice at 3.05\,\micron, CO$_{2}$ at 4.27\,\micron, and CO at 4.60\,\micron. The full decomposition result is shown in the inset.}
    \label{fig:pahfit_example}
\end{figure*}

\subsection{New Emission Components}
The profile of the 3.3\,\micron\ PAH is set following the suggestions from \citet{Tokunaga1991} and \citet{VanDiedenhoven2004}. The 3.3\,\micron\ PAH feature sometimes shows a broader wing towards longer wavelengths which is often attributed to a 3.4\,\micron\ aliphatic band, together with an even broader underlying plateau centered at 3.47\,\micron\ \citep[][see \S~\ref{sec:aliphatic}]{Hammonds2015}. The modest range of the redshift in our sample fills in the gap between the \textit{AKARI}/IRC and the \textit{Spitzer}/IRS spectra, which allows us to investigate the weaker PAH feature at $\sim$\,5.27\,\micron\ that had been reported in \textit{ISO}\ spectra \citep[e.g.][]{Boersma2009}. We find the 5.27\,\micron\ PAH feature can further be separated into two individual bands centered at 5.24 and 5.33\,\micron. All the newly added/modified dust feature parameters can be found in Table~\ref{Tab:additional_features}(a).

To account for the unresolved recombination lines within the \textit{AKARI}/IRC spectra, we added Br$\beta$ at 2.63\,\micron, Pf$\gamma$ at 3.74\,\micron, Br$\alpha$ at 4.05\,\micron, and Pf$\beta$ at 4.65\,\micron\ to \textsc{Pahfit}. Since the Pf$\beta$ line at 4.65\,\micron\ may be blended with an absorption band arising from CO at $\sim$ 4.67\,\micron\ (see \citet{Spoon2003} for a high resolution spectrum), to avoid such ambiguity, we tied Pf$\beta$ with Br$\alpha$ using the theoretical line ratio of Pf$\beta$/Br$\alpha$=0.20 under the Case B condition, with n$_{e} \sim 10^{2}-10^{7}$ cm$^{-3}$ and T$_{e} \sim (3-30) \times 10^{3}$K \citep{Storey1995}. The newly added line feature parameters can be found in Table~\ref{Tab:additional_features}(b).

The original \textsc{Pahfit}\ model treats starlight continuum as a single temperature blackbody spectrum with T$_{\star}$=5000K. This assumption is sufficient for $\lambda>5$\,\micron, but with wavelengths down to 2.5\,\micron, additional features are needed to account for the non-stellar NIR excess that has an average temperature $\sim$ 1000K from the hot grains \citep{Lu2003} and the divergence of the mean starlight spectrum at $\lambda \sim$2--5\,\micron. To address this and gain flexibility in the 2.5--4\,\micron\ continuum fit, we add 800K and 1500K modified blackbody dust components, with the $\sim$1500K component representative of the dust sublimation temperature in AGNs \citep{Marchall2007, Nenkova2008}. Accurately separating stellar and non-stellar continuum sources is challenging given our spectral coverage and simplified treatment, so no physical significance is ascribed to individual continuum components.

\input{tbl2_new_add_features}

\pagebreak
\subsection{New Absorption Features}
The expanded coverage of the \textit{AKARI}/IRC also includes several absorption features: a broad 3.05\,\micron\ H$_{2}$O ice feature \citep[e.g.][]{Gibb2004}, 4.27\,\micron\ CO$_{2}$ \citep{Spoon2000, Spoon2003}, and a very broad ($\sim$\,0.2\,\micron) feature of CO absorption centered at 4.67\,\micron\ \citep{Baba2018}. To precisely extract the 3.3\,\micron\ PAH, it is necessary to differentiate the band from the nearby 3.05\,\micron\ H$_{2}$O ice. Another H$_{2}$O ice absorption features at 6.02\,\micron\ is also included in the model, with its peak optical depth tied to the 3.05\,\micron\ H$_{2}$O ice. All of the newly added absorption features are approximated as Drude profiles with the widths and relative strengths of \citet{Gao2013} (see Table~\ref{Tab:additional_features}(c)), except for the FWHM of the 3.05\,\micron\ H$_{2}$O ice. We found our sample have a profile on average 30\% wider than reported by Gao et al., with FWHM = 0.41\,\micron. Since the profile of H$_{2}$O absorption is temperature dependent and its wing extends towards longer wavelength as temperature increases \citep{Spoon2002}, the upper limit of the FWHM is set to be 0.53\,\micron, slightly less than the even broader ice profile (FWHM = 0.56\,\micron) reported by \citet{Hammonds2015} when analyzing a large number of Galactic objects using \textit{AKARI}/IRC. We found no correlation between optical depths of H$_{2}$O and silicate, similar to what was reported in \citet{Doi2019}, so \ensuremath{\tau_\mathrm{H_{2}O}}\ and \ensuremath{\tau_\mathrm{Si}}\ were left uncoupled.

The CO ro-vibrational absorption profile (4.67 $\mu m$) has been used to model the column density and gas temperature of obscured AGNs \citep{Spoon2004, Shirahata2013, Baba2018}. As pointed by Baba et al., double-branched features of CO are observed, with the possible width of each branch extends to $\sim$\,0.2\,\micron, depending on the CO gas temperature. Here, for simplicity, a Drude profile with a central wavelength constrained to 4.5\,$<\lambda<$\,4.69 is adopted to account for such deep and wide feature. Neither CO$_2$ nor CO absorption feature affects dust emission component (only Pf $\beta$ can be affected), so these simplifying treatments have limited impact on our results.

To demonstrate the full range (2.5--38\,\micron) extinction profile that includes the newly added absorption features, we show PAHFIT results of the two example spectra (UGC~1385 and IRAS~F06076-2139) that contain the minimum and maximum \ensuremath{\tau_\mathrm{Si}}\ values in our bright-PAH sample in Fig.~\ref{fig:galaxy_ext_curve}. IRAS~F06076-2139 passed our selection criteria (\S~\ref{subsec:bright_PAH}) because its spline-based silicate strength is very close to but just below the cutoff \mbox{(S$_\mathrm{Si}$=-1.38)}.

\begin{figure}
 	\includegraphics[width=1.\columnwidth]{galaxy_ext_curve.pdf}
    \caption{Two PAHFIT results of sample spectra that exhibit the minimum/maximum \ensuremath{\tau_\mathrm{Si}}\ (mixed) in our bright-PAH sample. Note that  \ensuremath{\tau_\mathrm{Si}}, \ensuremath{\tau_\mathrm{H_{2}O}}, and $\tau_\mathrm{CO+CO_2}$ vary independently.}
    \label{fig:galaxy_ext_curve}
\end{figure}

\subsection{Attenuation in Different Geometries} \label{subsec:attenuation}
With a variety of smooth and spectrally resolved extinction components, the relative geometry of emission sources and obscuring material plays important roles in the interpretation of the MIR spectra of galaxies. For normal star-forming galaxies, it often considered appropriate to assume stars and dust are spatially well mixed, with both PAH+line and continuum emission subject to a similar range of attenuation\footnote{Here and throughout we use the term ``attenuation'' to refer to the combined effects of photon absorption and scattering, as well as the impact of geometry on these phenomena, reserving ``extinction'' for the geometry independent and wavelength dependent curves measured along individual lines of sight (typically in the MW).}. For galaxies powered by buried AGN or compact central starburst with only H\,\textsc{ii}-regions, however, the obscuration comes predominately from the central part of the system, while the PAH/line emissions excited by stars in the outer part are mostly not attenuated \citep{Imanishi2007, Imanishi2008}. To account for this latter scenario, we add a new geometry, hereafter the ``obscured continuum'' geometry, upon the fully mixed geometry that was already developed in \textsc{Pahfit}. For the extinction profile, we retain the custom curve developed with \textsc{Pahfit}\ \citep{Smith2007}, which derives from the Galactic center of \citet{Kemper2004}. With both the mixed and obscured continuum geometries, it is our attempt to bracket the range of realistic emission and absorption by these two simple schemes. In other words, the ``real" measurements should then be bracketed by these two extremes, and we interpret our results in both cases throughout this paper. We decide not to use the screen geometry as an alternative to the mixed geometry because the difference of the optical depths derived from these two schemes is simply a constant factor within our attenuation range, with no impact on our decomposition results.

Including the additional obscured continuum geometry in the model is necessary since \textsc{Pahfit}\ was trained primarily on a sample of nearby galaxies from the SINGS survey \citep{Kennicutt2003}, in which the attenuation is moderate \citep[with mixed silicate opacity, $\tau_{\mathrm{Si}} $$\lesssim$ 0.35;][]{Smith2007}. Even our trimmed sample, however, contains galaxies with higher attenuation (see \S~\ref{subsec:bright_PAH}) with the median $\tau_{\mathrm{Si}}\!\sim\!1.1$. This amount of attenuation may affect our measurement of the 3.3\,\micron\ PAH since it suffers the strongest effects of dust and ice attenuation (see \S~\ref{sec:3.3PAH} and \S~\ref{subsec:geometry_disc} for detailed discussion).

In the new obscured continuum geometry, all dust and line emission features sit atop an attenuated continuum as illustrated schematically in Fig.~\ref{fig:geometry}, where both geometries in mixed and obscured continuum are demonstrated. While \textsc{Pahfit}\ by default outputs recovered fluxes corrected from dust attenuation, the measurements obtained from the obscured continuum geometry are closer to the observed fluxes measured directly by integrating the line/band profile above the assumed local continuum. Note that in both geometry schemes, optical depth acts in a similar fashion with a scaling factor of (1-$e^{-\tau_{\lambda}}$)/$\tau_{\lambda}$.

\begin{figure}
 	\includegraphics[width=1.0\columnwidth]{geometry.pdf}
    \caption{The two bracketing geometries considered in this study. In the mixed geometry, emitting sources and dust are fully mixed, whereas in the obscured continuum geometry, PAH and line emission occurs outside an obscured central continuum zone, and is not subjected to attenuation.
    }
    \label{fig:geometry}
\end{figure}

\section{Results} \label{sec:result}
In this section, we first present a set of galaxy template spectra produced from ASESS and the spectral decomposition result of the galaxies in our bright-PAH sample. We then discuss the properties of 3.3\,\micron\ PAH together with star formation rate calibration using 3.3\,\micron\ and how different geometries may come into play. A special case of low-metallicity galaxy, II~Zw~40, is discussed. We then study the relationship between the aromatic and aliphatic molecules using the 3.3 and 3.4\,\micron\ bands and a size diagnostic using 3.3, 11.3, and 17\,\micron\ PAHs, followed by a comprehensive comparison between 3.3\,\micron\ PAH and other prominent PAH features. We end this section with an outlook for probing the 3.3\,\micron\ PAH strength using the \textit{JWST}\ photometric filters.

\subsection{The Composite Template Spectra} \label{sec:template}
From the complete cross-archival sample of galaxies with compatible \textit{Spitzer}/IRS and \textit{AKARI}/IRC coverage (378 objects), we construct average spectral templates covering the complete spectral range of 2.6--28\,\micron\ from the \citet{Spoon2007} galaxy classification
discussed in Fig.~\ref{fig:sample_selection}(a). Galaxy spectra in the eight classes are presented in Fig.~\ref{fig:template}, with each produced by the median of high-quality spectra (S/N $\geq$ 60 at 3.29\,\micron) in the respective class.

\begin{figure}
 	\includegraphics[width=1.0\columnwidth]{template_stripe}
    \caption{Galaxy templates of the eight classes defined in Fig.~\ref{fig:sample_selection}(a). All spectra are plotted in log scale and normalized to 6.2\,\micron. The gray areas indicate the 1$\sigma$ dispersion. Vertical dashed lines denote the central wavelengths of the main PAH emission features at 3.3, 6.2, 7.7, 8.6, 11.3, and 17\,\micron. Most galaxies in our moderate attenuation, bright-PAH sample are found within classes 1B, 2B, 1C, and 2C. Spectra are color-coded correspondingly to classes in Fig.~\ref{fig:sample_selection}(a). Template of class 3C is not provided since only one galaxy (NGC 4945) belongs to this class. The number of sources averaged for the templates range from 10 in 2A to 116 in 1A.}
    \label{fig:template}
\end{figure}

The bulk of galaxies in our bright-PAH sample lies within classes  1B, 1C, 2B, or 2C, with more than half in 1C, as shown in Fig.~\ref{fig:sample_selection}(a). The 1C spectrum in Fig.~\ref{fig:template} is a typical star-forming galaxy, showing prominent PAH emission, with weak to moderate 9.7\,\micron\ silicate absorption. As the spectrum shifts from class 1C to 2B, the 9.7\,\micron\ silicate absorption becomes more prominent and the PAH emission bands become more suppressed. For classes 3A and 3B, the absorption features, including H$_{2}$O, CO$_{2}$, CO, and silicates all reach their maximum depths, while most PAH features are diminished. But as an exception, the 3.3\,\micron\ PAH band remains apparent even when other strong PAH bands at 6.2, 7.7, and 8.6\,\micron\ disappear in the 3A class. The behavior of the 3.3\,\micron\ PAH in the most highly obscured systems will be a topic for a future study, but the impact of attenuation on this shortest wavelength PAH feature is considerable even at moderate levels of obscuration (see \S~\ref{subsec:PAH2Geoms}). The digital format of these templates can be found in CDS format in the online version of this publication.

\subsection{\textsc{Pahfit}\ Result}
\label{sec:result_pahfit}
In Fig.~\ref{fig:pahfit_example}, we demonstrate the capabilities of our enhanced \textsc{Pahfit}\ model on the composite, high-quality 1C spectrum (see Fig.~\ref{fig:template}). The data are shown in squares overlaid by the fitted profile in green. Dust features (blue) are modeled using Drude profile and emission lines (purple) are modeled with Gaussians. Newly added features are labelled. The \textsc{Pahfit}\ results of the bright-PAH sample in both the mixed and obscured continuum geometries can be found in Table~\ref{tab:PAHFIT_result_mx} and \ref{tab:PAHFIT_result_oc}, respectively.

\subsection{The 3.3\,\micron\ PAH} \label{sec:3.3PAH}
Thanks to \textit{Spitzer}, studies of the PAH bands at wavelengths $>$ 5\,\micron\ abound. However, investigations of 3.3\,\micron\ PAH has been limited due to the lack of \textit{Spitzer}/IRS wavelength coverage, the hindrance by the atmospheric background emission for ground observation, and the low sensitivity at 3\,\micron\ in \textit{ISO}. To put the 3.3\,\micron\ PAH in context, we use the \textit{AKARI}-\textit{Spitzer}\ combined spectra obtained from the 113 galaxies in our bright-PAH sample and investigate how the 3.3\,\micron\ PAH measurement varies under the assumption of different geometries, and as infrared luminosity changes.

\subsubsection{3.3\,\micron\ PAH and Attenuation Geometry} \label{subsec:PAH2Geoms}
Among all the PAH emission features, the 3.3\,\micron\ PAH has the shortest wavelength, and it sets itself apart from other prominent PAH emissions at longer wavelengths. Unlike the PAH band complexes at 6--9\,\micron, for which measuring band strengths can be challenging due to the blending of multiple components, measuring the strength of the 3.3\,\micron\ PAH seems to be straightforward. However, the coincidence with broad 3.05\,\micron\ H$_{2}$O ice absorption \citep{Gibb2004}, together with the 3.4\,\micron\ aliphatic feature, if present, affects the 3.3\,\micron\ PAH measurement. In particular, the H$_{2}$O ice feature has an impact since its optical depth is not negligible even for some star-forming galaxies ---  the average H$_{2}$O opacity for our bright-PAH sample is $\tau_\mathrm{H_{2}O}=$0.20. Note that in this study, we only consider the 3.4\,\micron\ aliphatic feature in emission, even though absorption feature can also be found in the ISM \citep{Pendleton2002, Doi2019}. For most standalone \textit{AKARI}/IRC studies, the 3.3\,\micron\ PAH strength is measured by placing a linear continuum and integrating the band atop \citep{Imanishi2008, Imanishi2010, Ichikawa2014, Inami2018}. The result produced by this approach only takes into account the attenuation within the 2.5--5\,\micron\ spectral range, but not the overall MIR attenuation predominately constrained by the 9.7\,\micron\ absorption arising from the amorphous silicates and raised steeply towards the NIR at $\sim$\,3\,\micron.

\input{tbl34_pahfit.tex}

\begin{figure}
 	\includegraphics[width=\columnwidth]{L33_LIR}
    \caption{Comparison of the integrated infrared luminosity (L$\mathrm{_{IR}}$) and the 3.3\,\micron\ PAH luminosity (L$\mathrm{_{PAH \ 3.3}}$) under the two different geometries shown in Fig.~\ref{fig:geometry}. Circles and triangles denote galaxies with \,L\ensuremath{_\mathrm{IR}} $< 10^{12}$ L$_{\odot}$ and \,L\ensuremath{_\mathrm{IR}}\ $\geq 10^{12}$ L$_{\odot}$ (ULIRGs), respectively. The solid lines are the median of the two binned distributions, with the 25--75\% ranges shaded. Both distributions fall between the dashed and dotted lines, representing L\ensuremath{_\mathrm{PAH\,3.3}}/\,L\ensuremath{_\mathrm{IR}}\ = 0.1\% and 0.05\%, a typical range also found in \citet{Imanishi2010} and \citet{Y2013}. Note that II~Zw~40 is not included in this plot due to its exceptionally low 3.3\,\micron\ PAH luminosity (L\ensuremath{_\mathrm{PAH\,3.3}}$\,\sim\,3 \times$10$^{5}$\,L\ensuremath{_\mathrm{\odot}}).}
    \label{fig:L33_LIR}
\end{figure}

To explore how attenuation can alter inferred PAH strengths, we treated two different geometric scenarios as mentioned in \S~\ref{subsec:attenuation} -- fully mixed and obscured continuum. Given the steep rise of the extinction curve towards $\sim$\,3\,\micron, we expect the 3.3\,\micron\ PAH to suffer the greatest impacts of dust obscuration. Even with our bright-PAH sample, we found the 3.3\,\micron\ PAH measured under the mixed geometry assumption can be boosted up to 5 times compared to the obscured continuum geometry. This ratio, r\,=\,L\ensuremath{_\mathrm{PAH\,3.3}}(mixed)/L\ensuremath{_\mathrm{PAH\,3.3}}(obscured continuum), can be approximated by a functional form: r\,=\,exp[(0.320\,$\pm$\,0.004)(\ensuremath{\tau_\mathrm{Si}}+\,\ensuremath{\tau_\mathrm{H_{2}O}})], with optical depths measured in the mixed geometry. In our bright-PAH sample, the median of \ensuremath{\tau_\mathrm{Si}}+\,\ensuremath{\tau_\mathrm{H_{2}O}}\ is 1.31, indicating, in principle, a factor of 1.5 in between the measured strengths in the mixed and obscured continuum schemes, in which the limit of the ``true'' value is set.

\begin{figure*}
 	\includegraphics[width=1.0\textwidth]{pah33_to_totpah_S24}
    \caption{The strength of the fractional 3.3\,\micron\ PAH against the infrared luminosity in two different adopted geometries. The fractional 3.3\,\micron\ PAH strength increases with $\mathrm{L_{IR}}$ in the mixed geometry (\textit{left}) but remains constant in the obscured continuum geometry (\textit{right}). The data points are color-coded by the f$_{\nu}$(24\,\micron)/f$_{\nu}$(2.7\,\micron) ratio. The horizontal dashed lines indicate the nominal value adopted in the dust model of \citet{Draine2007}.}
    \label{fig:pah33_to_totpah}
\end{figure*}

\subsubsection{3.3\,\micron\ PAH and \,L\ensuremath{_\mathrm{IR}}}
The upper bound of the observed ratio of 3.3\,\micron\ PAH to the integrated infrared luminosity (L\ensuremath{_\mathrm{PAH\,3.3}}/\,L\ensuremath{_\mathrm{IR}}) for star-forming LIRGs has been found to be $\sim$ 10$^{-3}$ \citep{Mouri1990, Imanishi2010}. For ULIRGs (L$_\mathrm{IR} \geqslant 10^{12}$ \,L\ensuremath{_\mathrm{\odot}}), however, a declining trend is seen, with  bolometric fractional luminosity of the uncorrected 3.3\,\micron\ PAH (L\ensuremath{_\mathrm{PAH\,3.3}}/\,L\ensuremath{_\mathrm{IR}}) dropping due to dust attenuation or the overall destruction of PAH molecules under intensive UV radiation \citep{Imanishi2010, Y2013}. Our results broadly agree with this trend. Fig.~\ref{fig:L33_LIR} shows binned 3.3\,\micron\ PAH luminosity in two different geometries. In the obscured continuum geometry, the distribution lies in between the 0.1 and 0.05\% lines and gradually enters below the 0.05\% line at \,L\ensuremath{_\mathrm{IR}}$\,\sim\,5 \times 10^{11}$ \,L\ensuremath{_\mathrm{\odot}}. However, in the mixed geometry, rather than showing a 3.3\,\micron\ PAH downturn, L\ensuremath{_\mathrm{PAH\,3.3}}/\,L\ensuremath{_\mathrm{IR}}\ remains relatively constant over the whole range of \,L\ensuremath{_\mathrm{IR}}. This is because the corrected strengths of the 3.3\,\micron\ PAH in ULIRGs typically increase $\sim\,2 \times$ in the mixed geometry compared to the obscured continuum geometry. The difference between the two distributions suggests the high optical depth in ULIRGs likely causes the PAH deficit shown in the obscured continuum geometry.

\subsubsection{Fractional 3.3\,\micron\ PAH}
Changes in PAH luminosity as a fraction of bolometric power can depend on both PAH power as well as variation in non-stochastic dust heating (AGN, etc.), but fractional PAH power relative to total PAH luminosity is more closely tied to properties of the PAH grains and their heating conditions. Fig.~\ref{fig:pah33_to_totpah} shows the relationship between the fractional 3.3\,\micron\ PAH (L\ensuremath{_\mathrm{PAH\,3.3}}/L\ensuremath{_\mathrm{\Sigma\,PAH}}) and \,L\ensuremath{_\mathrm{IR}}. The fractional 3.3\,\micron\ PAH strength increases with \,L\ensuremath{_\mathrm{IR}}\ in the fully mixed geometry. For galaxies with \,L\ensuremath{_\mathrm{IR}}\,$<10^{12}$~\,L\ensuremath{_\mathrm{\odot}}, the averaged fractional 3.3\,\micron\ PAH is 2.06$\pm$0.59\%, whereas for \,L\ensuremath{_\mathrm{IR}}\,$\geq\,10^{12}$~\,L\ensuremath{_\mathrm{\odot}}, this value increases to 2.74$\pm$0.97\%. On the contrary, for the obscured continuum geometry, the fractional 3.3\,\micron\ PAH power show little change with luminosity, with a similar average above and below $10^{12}$\,L\ensuremath{_\mathrm{\odot}}\ of 1.50$\pm$0.38\%. We found galaxies with a steeper continuum slope ($f_\nu(24\,\mu m)/f_\nu(2.7\,\mu m)\,>\,60$) are associated with larger discrepancies among the two geometries because they are more likely dusty, obscured sources. The horizontal line in the figure indicates the nominal ratio of L\ensuremath{_\mathrm{PAH\,3.3}}\ to L\ensuremath{_\mathrm{\Sigma\,PAH}}\ in the \citet{Draine2007} dust model, in which the integrated strength of the 3.3\,\micron\ PAH was increased by a factor of 2 to match the theoretical calculation. The fractional 3.3\,\micron\ PAH obtained from the dust model is relatively high compared to our result, despite only some highly obscured ULIRGs are likely to have such a strong 3.3\,\micron\ PAH emission measured in the mixed geometry.

We note that low metallicity galaxy II~Zw~40 is a special case, which exhibits a steep 24/2.7 continuum slope, but the least \,L\ensuremath{_\mathrm{IR}}\ among our sample. The L\ensuremath{_\mathrm{PAH\,3.3}}/L\ensuremath{_\mathrm{\Sigma\,PAH}}\ in II~Zw~40 in the mixed and obscured continuum geometries are 2.8\% and 2.5\%, respectively. It lies at the upper end of the distribution in the obscured continuum geometry, despite being being a blue compact, low-metallicity galaxy, where small grains are expected to be under-abundant (see \S~\ref{subsec:IIZw40}).

\subsection{3.3\,\micron\ PAH as a Star Formation Rate Indicator}
\label{Sec:SFR}
One key to assessing the evolutionary history of galaxies is the link between their star formation rate and the physical condition in the ISM. Investigations into PAH emission as an indicator of star formation have been extensive \citep{Forster_Schreiber2004, Peeters2004a, Shipley2016, Xie2019}, but caution is needed when applying to low-metallicity galaxies, as PAH abundance is mostly driven by metallicity \citep{Calzetti2007}. By carefully treating galaxies with low-metallicity, \citet{Shipley2016} used a sample of star-forming galaxies and calibrated the extinction-corrected PAH bands at 6.2, 7.7, and 11.3\,\micron\ to the SFR measured by H$\alpha$+24\,\micron. More recently, \citet{Xie2019} provided a calibration using integrated 5--15\,\micron\ PAH emission against the SFR measured by [Ne\,\textsc{ii}]\ 12.8\,\micron\ and [Ne\,\textsc{iii}]\ 15.6\,\micron. These studies, however, used stronger PAH bands accessible from \textit{Spitzer} but not the 3.3\,\micron\ PAH. Studies of the reliability of 3.3\,\micron\ PAH as a SFR indicator are relatively scarce \citep{Kim2012, Y2013, Murata2017}. All these studies found 3.3\,\micron\ PAH correlates with \,L\ensuremath{_\mathrm{IR}}, but the correlation deviates at high luminosities (\,L\ensuremath{_\mathrm{IR}}\ $\gtrsim 10^{12}$\,L\ensuremath{_\mathrm{\odot}}), where \,L\ensuremath{_\mathrm{IR}}\ may no longer be a direct SFR indicator but be contaminated by the contribution of AGN to the infrared luminosity. Alternatively, small PAHs traced by the 3.3\,\micron\ PAH may be destroyed due to the energetic radiation field originating from the AGN in the host galaxy \citep{Genzel1998, Smith2007, ODowd2009, Diamond-Stanic2010, Wu2010}. In this study, we calibrate the SFR using the 3.3\,\micron\ PAH and study how the results may differ for different assumptions of emission/attenuation geometry. It is important to note that this analysis only pertains to our pre-selected star-forming sample, and therefore does not pertain to galaxies that are metal-poor and/or AGN-dominated. For example, low-metallicity galaxy, II~Zw~40, shows a very different star-forming property as discussed in the next section.

\subsubsection{A New Calibration of SFR using 3.3\,\micron\ PAH} \label{subsubsec:SFR_calibration}
Given the selection criteria described in \S~\ref{sec:sample_selection}, galaxies in our sample are primarily star-forming galaxies, with at most minor AGN contribution. To assess whether the 3.3\,\micron\ PAH is a reliable SFR indicator is crucial and will benefit the study offered by the impending \textit{JWST}/NIRSpec with wavelength coverage ranging from 0.6--5\,\micron\ with an unprecedented spectral resolution up to R=2700.
Other accessible SFR indicators in our \mbox{2.5--38\,\micron} range are the neon lines ([Ne\,\textsc{ii}]\ 12.81\,\micron\ and [Ne\,\textsc{iii}]\ 15.56 \,\micron). To investigate the behavior of 3.3\,\micron\ PAH emission as with SFR, we utilize the 102 galaxies from our bright-PAH subset with $>$3$\sigma$ detection in both [Ne\,\textsc{ii}]\ and [Ne\,\textsc{iii}]. To estimate SFR inferred from [Ne\,\textsc{ii}]\ and [Ne\,\textsc{iii}]\ (SFR$_\mathrm{Ne}$), we adopt the neon-based SFR relation of \citet{Zhuang2019}, an updated version of \citet{Ho2007} which explicitly includes the effect of metallicity. Here we adopt solar metallicity ($Z=Z_{\odot}$) for simplicity and normalize the SFR indicator to a Kroupa IMF \citep{Kroupa2001} with mass cutoffs at 0.1 and 100\,M$_{\odot}$.

\begin{figure}
 	\includegraphics[width=1.0\columnwidth]{SFR_pah33_Ne}
    \caption{\textit{Top}: The relationship between the line and PAH luminosities L\ensuremath{_\mathrm{[Ne\,II] + [Ne\,III]}}\ and L\ensuremath{_\mathrm{PAH\,3.3}}, in both the mixed and obscured continuum geometries. The two color-coded lines are best fits for the corresponding geometry. \textit{Bottom}: The calibration of the SFR using 3.3\,\micron\ PAH. Each data points is the weighted average of the values obtained from the mixed and obscured continuum geometries. Galaxies subject to only moderate optical depth (\ensuremath{\tau_\mathrm{Si}}+\ensuremath{\tau_\mathrm{H_{2}O}} $<$ 0.9) are shown as filled circles, and the dashed line and shading show the best weighted fit and 1$\sigma$ scatter (see Eq.~\ref{eqn:SFR}). The dotted line (at the upper end of the shading) indicates the average value from galaxies with (\ensuremath{\tau_\mathrm{Si}}+\ensuremath{\tau_\mathrm{H_{2}O}} $>$ 0.9). Two galaxies (II~Zw~40 and IRAS 15206+3342) exhibiting extremely low L\ensuremath{_\mathrm{PAH\,3.3}}/SNR ratios (below 3$\sigma$ from the calibrate correlation) are omitted from our calibration.}
    \label{fig:SFR_pah33_Ne}
\end{figure}

The top panel of Fig.~\ref{fig:SFR_pah33_Ne} shows the correlation between L\ensuremath{_\mathrm{[Ne\,II] + [Ne\,III]}}\ and L\ensuremath{_\mathrm{PAH\,3.3}}\ in the mixed and obscured continuum geometries, with scatters of 0.20 and 0.17 dex, respectively about the best fits. The bottom panel shows the relation of L\ensuremath{_\mathrm{[Ne\,II] + [Ne\,III]}}\ and the ratio of L\ensuremath{_\mathrm{PAH\,3.3}}\ to SFR$_\mathrm{Ne}$, calculated from the weighted average of the values obtained from both geometries. Aside from several outliers, the ratio of L\ensuremath{_\mathrm{PAH\,3.3}}\ to SFR$_\mathrm{Ne}$ shows no variation as L\ensuremath{_\mathrm{[Ne\,II] + [Ne\,III]}}\ spans $\sim$3 orders of magnitude, suggesting L\ensuremath{_\mathrm{PAH\,3.3}}\ and SFR$_\mathrm{Ne}$ are well correlated. The horizontal dashed line shows the best fit to the SFR calibration, which can be translated to a new calibration of SFR using the 3.3\,\micron\ PAH:

\begin{equation}
\mathrm{log} \bigg(\mathrm{\frac{SFR}{M_{\odot}yr^{-1}} \bigg) = -(6.80 \pm 0.18) + log \bigg(\frac{L_{PAH \ 3.3}}{L_{\odot}} \bigg)}
\label{eqn:SFR}
\end{equation}

To avoid the effects of attenuation on this calibration, only galaxies with modest obscuration \ensuremath{\tau_\mathrm{Si}}+\ensuremath{\tau_\mathrm{H_{2}O}}\,$<$\,0.9 are used (see \S \ref{subsubsec:SFR_geometry} for detail). Galaxies with \ensuremath{\tau_\mathrm{Si}}+\ensuremath{\tau_\mathrm{H_{2}O}}\,$\geq$\,0.9 are not considered for SFR calibration since collectively they have the ratio of L\ensuremath{_\mathrm{PAH\,3.3}}/SFR$_\mathrm{Ne}$ that is $\sim$20\% higher (dotted line) than the value adopted. Also omitted are galaxies with peculiar low L\ensuremath{_\mathrm{PAH\,3.3}}/SFR$_\mathrm{Ne}$ ratios, below the 3$\sigma$ limit, namely II~Zw~40 and IRAS 15206+3342. We stress that our calibration applies only to star-forming, metal-rich galaxies. Metal-poor galaxies such as II~Zw~40 appears to have a ratio of L\ensuremath{_\mathrm{PAH\,3.3}}/SFR$_\mathrm{Ne}$ that is a factor of $\sim$\,8 lower than the average of our bright-PAH sample (see \S~\ref{subsec:IIZw40}). IRAS~15206+3342 lies very close to the AGN selection area in Fig.~\ref{fig:sample_selection}(b), and its nucleus is dominated by H\,\textsc{ii}\ regions and surrounded by shock-induced LINER ionization \citep{Arribas2002} that may suppress PAH emission. Clearly, PAHs lose their diagnostic power as a SFR indicator in environments such as low-metallicity or AGN/LINER-dominated regions.

\subsubsection{SFR Calibration in Different Geometries}
\label{subsubsec:SFR_geometry}
To explore how geometry affects 3.3\,\micron\ PAH emission as a tracer for ongoing star formation, we examine how the ratio of L\ensuremath{_\mathrm{PAH\,3.3}}\ to SFR$_\mathrm{Ne}$ varies with respect to $\tau_\mathrm{Si}$+$\tau_\mathrm{H_{2}O}$ in both geometries. Fig.~\ref{fig:SFR_calibration} shows the two median trends in the binned data. Both trends exhibit no difference at low $\tau$ and progressively diverge as $\tau$ increases. The 3.3\,\micron\ PAH luminosity per unit star formation rate stays rather consistent in the obscured continuum geometry within $0\,<\,\tau\,<\,4$, but monotonically increases together with $\tau$ in the mixed geometry. To test up to which $\tau$ value can the two populations be regarded as statistically similar, we perform a two-sample t-test, choosing a p-value significance level equal to 0.05, and find $\tau$=0.9 is the limit where the two populations can still count as identical. we therefore conclude the two trends exhibit no significant difference at $\tau <$ 0.9, indicated by the dashed line. Note that the ratio of L\ensuremath{_\mathrm{PAH\,3.3}}/\,L\ensuremath{_\mathrm{IR}}\ for those galaxies with \ensuremath{\tau_\mathrm{Si}}+\ensuremath{\tau_\mathrm{H_{2}O}} $<$ 0.9 in the mixed and obscured continuum geometries are 0.11$\pm$0.05\% and 0.09$\pm$0.04\%, respectively. The offset between the two trends becomes as large as $\sim3\times$ in the most obscured cases in our bright PAH sample (\ensuremath{\tau_\mathrm{Si}}+\ensuremath{\tau_\mathrm{H_{2}O}}\,$\sim$\,4). Whether the two trends will continue diverging at even greater optical depths requires further study on more deeply embedded galaxies.

\begin{figure}
 	\includegraphics[width=1.0\columnwidth]{SFR_calibration.pdf}
    \caption{The dependence of the L\ensuremath{_\mathrm{PAH\,3.3}}/SFR ratio to the sum of silicate and water optical depths, in the mixed (red) and obscured continuum (blue) geometries. The trends are the median of the binned data with inner two quartiles shaded. Our SFR calibration is measured using galaxies with moderate obscuration ($\tau_\mathrm{Si}$+$\tau_\mathrm{H_{2}O}< 0.9$, left of the dashed line), where the two trends converge. The low metallicity outlier galaxy II~Zw~40 is indicated. }
    \label{fig:SFR_calibration}
\end{figure}

\begin{figure*}
 	\includegraphics[width=1.0\textwidth]{IIZw40_ExtraEmissionLine}
    \caption{The extended \textsc{Pahfit}\ model of low metallicity galaxy II~Zw~40. The prominent 3.3\,\micron\ PAH together with other emission lines and their corresponding wavelengths are labelled (with the dashed line indicating He I line at 4.30\,\micron). Our 3.3\,\micron\ PAH flux in this galaxy has been corrected by removing the contribution of the indistinguishable recombination line Pf~$\delta$ at 3.297\,\micron\ that blends with the 3.3\,\micron\ feature, contributing $\sim$\,21\% of the uncorrected 3.3\,\micron\ PAH power.}
    \label{fig:IIZw40}
\end{figure*}

\subsection{Low Metallicity Case Study: II~Zw~40}
\label{subsec:IIZw40}
Large PAHs are expected to have a higher survival rate than small PAHs in a low-metallicity ISM \citep{Hunt2010}, since PAHs as a whole are subjected to enhanced photon interaction in the hard, intense interstellar radiation field which characterize as dust-poor environments \citep{Engelbracht2005, Galliano2005, Madden2006, Wu2006, Gordon2008}. However, a counter example has been found in SINGS galaxies, where L\ensuremath{_\mathrm{PAH\,17}}/L\ensuremath{_\mathrm{\Sigma\,PAH}}\ decreases with decreasing metallicity \citep{Smith2007}. Also, supporting such scenario, study by \citet{Sandstrom2012} found the size distribution of PAHs are apparently shifted to smaller sizes in the Small Magellanic Cloud relative to higher metallicity systems. These authors argued that the size shift towards smaller PAHs cannot be the result of ISM processing because small PAHs are more prone to photo-destruction under such conditions. Rather, PAHs at low-metallicity although experience more overall destruction, they have inherently smaller size distribution, making the PAH index $q_\mathrm{PAH}$, the fraction of PAHs relative to the total amount of dust, to have a positive correlation with the metallicity. Galaxies with 12+log$_{10}$(O/H)\,$<$\,8.1 typically have $q_\mathrm{PAH}$=1.0\%, a value that is at least three times less than that measured in galaxies with metallicity 12+log$_{10}$(O/H)\,$>$\,8.1 \citep{Drain2007b}. Given the many impact of low metal content on the conditions of the ISM, ascribing the PAH emission deficiency at low-Z to an exact mechanism may be challenging, but our sample galaxy II~Zw~40 provides a case for understanding the effect of metallicity on the formation/destruction of PAH molecules incorporating all the PAH emission bands.

II~Zw~40 is a nearby (D=10 Mpc), prototypical blue compact dwarf galaxy (BCD), containing a super star cluster and dominated by a giant H\,\textsc{ii}\ region with diameter $\sim$0.5 kpc \citep{Galliano2005}. Its metallicity is log(O/H)+12\,=\,8.09 \citep{Guseva2000}, or Z\,$\sim$\,1/4 Z$_{\odot}$. Fig.~\ref{fig:IIZw40} shows the extended \textsc{Pahfit}\ result of II~Zw~40. The overall PAH emission is highly suppressed in II~Zw~40, such that its ratio of L\ensuremath{_\mathrm{\Sigma\,PAH}}\ to \,L\ensuremath{_\mathrm{IR}}\ is less than that measured in other sampled galaxies by a factor of $\sim$\,10. However, we find relatively strong 3.3\,\micron\ PAH emission (L\ensuremath{_\mathrm{PAH\,3.3}}/L\ensuremath{_\mathrm{\Sigma\,PAH}}\,$\sim$\,2.8\%) with the absence of 17\,\micron\ PAH, suggesting a size distribution shifted towards the smaller end.

Because the 5--9\,\micron\ PAH features are weak compared to the continuum, and due to the limited sensitivity of PAHFIT to opacities \ensuremath{\tau_\mathrm{Si}}\,$\lesssim$\,0.2 \citep{Smith2007}, the total PAH power in II~Zw~40 is somewhat uncertain.  We estimate \ensuremath{\tau_\mathrm{Si}}\,$\simeq$\,0.2, consistent with the CO$_{2}$ and CO absorption bands near 4.5\,\micron. As a simpler alternative with no dependence on attenuation or continuum modeling, we compare our result to a spline-based decomposition similar to the method employed in \citet{Peeters2002}. The estimated fractional 3.3\,\micron\ PAH power with a spline-fitting technique is  L\ensuremath{_\mathrm{PAH\,3.3}}/L\ensuremath{_\mathrm{\Sigma\,PAH}}\,=\,6.6\%, suggesting that L\ensuremath{_\mathrm{PAH\,3.3}}/L\ensuremath{_\mathrm{\Sigma\,PAH}}\ from the full decomposition may be a lower limit. Additional evidence is found in the band ratios of isolated PAH features. Unlike ratios involving the integrated power of broad PAH features, whose strengths are often underestimated by spline-fitting methods  \citep[e.g.][see Table 6]{Smith2007}, the ratio between the 3.3 and 11.3\,\micron\ PAH bands is robust against decomposition method, since these features are narrower and more isolated from other prominent PAH bands. The 3.3/11.3 ratio of 0.36 in II~Zw~40 is among the highest in the sample, and this ratio shows $<$\,20\% difference between the two estimation methods.

II~Zw~40 has an unusually strong series of recombination lines. In addition to lines that have incorporated in \textsc{Pahfit}\ (see Table~\ref{Tab:additional_features}), five line features, including Pf 11, Pf$\epsilon$, Hu 13, He I, and Hu 12, are identified in the spectrum as shown in Fig.~\ref{fig:IIZw40}. The 3.3\,\micron\ PAH strength may be overestimated since it is blended with the Pf$\delta$ line at 3.297\,\micron. Because differentiating Pf$\delta$ line from the 3.3\,\micron\ PAH is difficult with our spectral resolution, we correct the 3.3\,\micron\ PAH strength in II~Zw~40 assuming the strength of Pf$\delta$ to be 69.71\% that of Pf$\gamma$ \citep{Hummer1987}. In such case, this results in a reduction of 3.3\,\micron\ flux by 20.6\%. This is an unusually large correction compared to the typical value inferred from our 1C template spectrum of just 0.9\%. We note that this Pf$\delta$ correction of the 3.3\,\micron\ PAH power has \emph{only} been applied to II~Zw~40 due to its exceptionally strong recombination lines exhibited in the spectrum. Also, this correction in II~Zw~40 has been applied throughout the results of this paper, and the corrected 3.3\,\micron\ luminosity is in Table~\ref{tab:PAHFIT_result_mx} and Table~\ref{tab:PAHFIT_result_oc}.

\subsection{Aromatic Versus Aliphatic} \label{sec:aliphatic}
Understanding the abundance of aromatic and aliphatic components of the ISM can shed light on the evolution of dust grains and provide important constraints on dust models. The 3.3\,\micron\ PAH emission is sometimes accompanied by a weak 3.4\,\micron\ feature often attributed to the C-H vibrational modes in aliphatic hydrocarbons \citep{Joblin1996, Pendleton2002, Yang2016}. Anharmonicity of the aromatic C-H stretch has also been suggested \citep{Barker1987}, and more recently, based on high-resolution IR spectra of jet-cooled ``decorated'' PAHs, \citet{Maltseva2018} suggested that these hydrogenated and methylated PAHs are likely the carriers of the 3.4\,\micron\ feature. In addition, the plateau feature centered at 3.47\,\micron\ that extends up to $\sim$\,3.6\,\micron\ (see Fig.~\ref{fig:pahfit_example}) is tightly correlated with 3.3\,\micron\ PAH emission, suggesting the 3.47\,\micron\ feature also arises from aromatic carriers \citep{Hammonds2015, Pilleri2015}.

\begin{figure*}
 	\includegraphics[width=1.0\textwidth]{aromatic_and_aliphatic.pdf}
    \caption{Correlations between the measured 3.3\,\micron\ PAH, 3.4\,\micron\ aliphatic emission, and 3.47\,\micron\ plateau. Dispersions of the three correlations are 0.18, 0.25, and 0.19 dex, respectively. The dashed line denotes the fit of the data, with the red line indicating the proportional relationship. A tighter correlation between 3.3\,\micron\ PAH and 3.47\,\micron\ plateau is found, suggesting the plateau is dominated by the aromatic emission. The contours represent 90\%, 70\%, 50\%, 30\%, and 10\% enclosed probability.}
    \label{fig:aromatic_and_aliphatic}
\end{figure*}

We investigate the relationships among the 3.3\,\micron\ aromatic emission, 3.4\,\micron\ aliphatic emission, and the broad 3.47\,\micron\ plateau in our bright-PAH sample. Figure~\ref{fig:aromatic_and_aliphatic} shows the correlations between all three combinations of the 3.3, 3.4, and 3.47\,\micron\ band fluxes. In each panel, the correlation is calculated for galaxies with $>$\,3$\sigma$ detections in both bands. The sample size and the Pearson correlation coefficient are listed. We find the 3.3 and 3.47\,\micron\ bands are best correlated (with r=0.95), supporting the claim that the plateau feature originates from aromatic bands \citet{Pilleri2015}. The 3.4\,\micron\ aliphatic feature is relatively weak in our sample --- L$\mathrm{_{3.4}}$/L$\mathrm{_{3.3}}=0.10\,\pm\,0.04$ (either geometry).  This limits the sample size with robust 3.4\,\micron\ detection given the sensitivity of \textit{AKARI}/IRC to $\sim$\,33\% of the galaxies in our bright-PAH sample. The 3.4\,\micron\ aliphatic is less correlated with the 3.47\,\micron\ plateau and the 3.3\,\micron\ PAH band, with r=0.87 and r=0.91, respectively.

\begin{figure}
 	\includegraphics[width=1.0\columnwidth]{34over33_30to15.pdf}
    \caption{Correlation between L$_{3.4}$/L$_{3.3}$ (aliphatic to aromatic) and the monochromatic continuum at 15 and 30\,\micron, which correlates with radiation field intensity. We divide our data into three equal-size sub-samples of the continuum slope, with the crosses indicating the weighted average. The error bars indicate the range of the continuum slope bins (horizontal) and the 1$\sigma$ scatter (vertical). We find a weak correlation with decreasing L$_{3.4}$/L$_{3.3}$ as  radiation fields increases. Galaxies with $<3\sigma$ detections of the L$_{3.4}$/L$_{3.3}$ ratio are displayed as limits.}
    \label{fig:34over33_30to15}
\end{figure}

\citet{Yamagishi2012} found that the 3.4/3.3 ratio decreases towards the center of starburst galaxy M82, suggesting that the relative abundance of the aliphatic band carriers drops towards regions with more intense star formation. The decline of the 3.4/3.3 ratio in PDRs is consistent with this picture, suggesting an efficient photochemical process leads to the destruction of the aliphatic subgroups present on the periphery of PAHs \citep{Joblin1996, Pilleri2015}. In Fig.~\ref{fig:34over33_30to15}, we compare the 3.4/3.3 ratio with the intensity of the radiation field inferred from the continuum slope between 15 and 30\,\micron, where wavelengths are free from PAH and line emissions. The slope of the continuum is an often used surrogate for radiation field intensity, with steeper slopes implying stronger radiation in starburst galaxies \citep{Groves2008}. We find a mild correlation between the two. The 3.4/3.3 ratio ranges between $\sim$\,0.03--0.2 and drops by 40\% between the subsets with the weakest and strongest radiation fields. This agrees with the findings of \citet{Hammonds2015} and \citet{Pilleri2015}, who pointed out that aliphatics, with their fragile C-H bonds, are prone to efficient photo-destruction, resulting in higher relative abundance of the aromatics in intense radiation environments.

\begin{figure*}
 	\includegraphics[width=\textwidth]{pah33andpah17_marginal.pdf}
    \caption{Size-sensitive PAH band ratios of 3.3/11.3 versus 17/11.3 in the mixed (left) and obscured continuum (right) geometries, color-coded with the integrated infrared luminosity. The marginalized distributions are also shown, with red and blue bars indicating the corresponding medians in the mixed and obscured continuum geometries, respectively, together with the yellow and white asterisks indicating the intersections of the bars. The dashed lines indicate the nominal values for individual PAH molecules from the model of \citet{Draine2007}, with tracks indicating neutral (purple) and ionized (red) PAHs with varying number of carbon atoms, N$_\mathrm{c}$. Detail discussion pertaining to the offset is in \S~\ref{sec:PAH_band}. Also shown are the attenuation tracks for these band ratios, parameterized by the silicate and water optical depths adopting the fully mixed dust geometry.
    }
    \label{fig:pah33andpah17}
\end{figure*}

\subsection{Size diagnostics with the 3.3\,\micron\ PAH} \label{sec:PAH_band}
The smallest PAH grains have the least heat capacities, and can therefore be stochastically heated to higher temperatures by single photon. As a result, small PAHs tend to emit efficiently at shorter wavelengths (3--9\,\micron), while their larger counterparts tend to dominate at longer wavelengths \citep[15--20 $\mu m$;][]{Draine2007}. A widely used method for probing the grain size is to use the band intensity ratio 6.2/7.7 \citep{Draine2001}. However, since these two bands are so close in wavelength to each other, the range of grain sizes they probe is similar. With the full range of PAHs accessible in our study, we consider the relative strengths of the 3.3 and 17\,\micron\ PAH bands as a probe of the size distribution, as these bands trace the smallest and largest ends of the PAH size distribution \citep{Schutte1993}. We note that In \textsc{Pahfit}, the 17\,\micron\ complex is a combination of a group of PAH bands at 16.4, 17.0, 17.4 and 17.9\,\micron, which is a broad PAH band feature thought to originate from PAH C-C-C banding modes \citep{Smith2004, Werner2004, Peeters2004b} and arise in PAHs with a few 1000 carbon atoms \citep{Draine2007}.

Fig.~\ref{fig:pah33andpah17} shows the ratios of both 3.3 and 17\,\micron\ PAH normalized to the 11.3\,\micron\ band for our two bracketing geometries, color-coded by \,L\ensuremath{_\mathrm{IR}}. Overlaid are tracks of the model PAH band ratios from \citet{Draine2007} for individual neutral and ionized PAH grains, of varying carbon atom count (N$_\mathrm{c}$). The bulk of our sample lies at $\sim$\,80\,$<$\,N$_\mathrm{c}$\,$<$\,200. Galaxies with higher \,L\ensuremath{_\mathrm{IR}}\ are likely to have a higher 3.3/11.3 ratio. We divide our sample into two groups and find galaxies with \,L\ensuremath{_\mathrm{IR}} $<$ 10$^{12}$\,L\ensuremath{_\mathrm{\odot}}\ and \,L\ensuremath{_\mathrm{IR}}$ \geq 10^{12}$\,L\ensuremath{_\mathrm{\odot}}\ have median of 3.3/11.3 as 0.20$^{+0.05}_{-0.04}$ and 0.29$^{+0.10}_{-0.07}$ in the mixed geometry, and 0.16$^{+0.03}_{-0.02}$ and 0.18$^{+0.04}_{-0.03}$ in the obscured continuum geometry, respectively (with the uncertainties indicating the inner two quartiles). Also shown are the dashed lines, indicating the nominal value derived from an averaged spectrum of PAH molecules with different N$_\mathrm{c}$'s, for a radiative heating intensity $U=1$, and assuming the size distribution from \citet{Weingartner2001} as adapted by \citet{Draine2007}\footnote{The average band ratios do not fall along the single-grain tracks due to large grains which produce weak PAH bands with moderate 17\,\micron\ power (and hence high 17/11.3 ratios).}. Only a few number of galaxies with boosted 3.3\,\micron\ PAH power in the mixed geometry can reach the nominal value in the x-axis, while no galaxy appear to have such high 3.3/11.3 in the obscured continuum geometry. In general, our 3.3\,\micron\ PAH flux measurements are at least a factor of 2 less than the nominal value, as we have demonstrated in Fig.~\ref{fig:pah33_to_totpah}.

\begin{figure*}
 	\includegraphics[width=1.0\textwidth]{PAH_trend.pdf}
    \caption{Changes of the main PAH band luminosity with respect to different \,L\ensuremath{_\mathrm{IR}}\ under the two geometries. The trends show the median of the binned data in our bright-PAH sample, with the error bars indicating the inner two quartiles. While all L$_\mathrm{PAH}$/L$_\mathrm{IR}$ decline with increasing \,L\ensuremath{_\mathrm{IR}}, diverse trends are shown in L$_\mathrm{PAH}$/L\ensuremath{_\mathrm{\Sigma\,PAH}}. }
    \label{fig:pah_trend}
\end{figure*}

More galaxies have higher 3.3/11.3 ratio in the mixed geometry than in the obscured continuum geometry suggest the 3.3\,\micron\ PAH is boosted by attenuation correction, leading to a spread towards the right side in the fully mixed geometry. To study impact of the attenuation on migrating the points between the two adopted geometries, we use the PAH band profiles decomposed from the 1C spectrum and impose attenuation on those features to produce the attenuation tracks that depict the change of the band ratios with varying \ensuremath{\tau_\mathrm{Si}}\ and \ensuremath{\tau_\mathrm{H_{2}O}}. We show the black and blue lines indicate the attenuation trends assuming \ensuremath{\tau_\mathrm{H_{2}O}}=0 and \ensuremath{\tau_\mathrm{H_{2}O}}=1.5, respectively, with corresponding \ensuremath{\tau_\mathrm{Si}}\ shown in each tick. These $\tau$ values are all parameterized by assuming a mixed geometry. The optical depth of H$\mathrm{_{2}O}$ has impact solely on 3.3/11.3 since the 3.3\,\micron\ PAH is the only band here that affected by \ensuremath{\tau_\mathrm{H_{2}O}}. The tracks indicate galaxies in the obscured continuum geometry are likely migrate towards the bottom right in the mixed geometry, depending on the levels the attenuation. We note that these tracks depict the normal behaviors of the galaxies, certain individuals may not follow such tracks, e.g., ULIRGs that are heavily obscured.

\subsection{PAH 3.3\,\micron\ Emission in Context} \label{sec:all_the_PAHs}
A deficiency of PAH emission with increasing IR luminosity has been found in a variety of environments, including starburst galaxies and galaxies with moderate to strong AGN contribution; in particular, L$_\mathrm{PAH}$/\,L\ensuremath{_\mathrm{IR}}\ drops significantly at ULIRG luminosities \citep{Wu2010, Pope2013, Stierwalt2014}. To assess whether a similar degree of decline is found between the 3.3\,\micron\ PAH power, and how the emission geometry impacts the result, in Fig.~\ref{fig:pah_trend} (a) and (b), we show the relationship of \,L\ensuremath{_\mathrm{IR}}\ and L$_\mathrm{PAH}$/\,L\ensuremath{_\mathrm{IR}}\ in the two geometries, in three luminosity bins -- infrared galaxies (IRGs: \,L\ensuremath{_\mathrm{IR}}\,$<$\,10$^{11}$\,\,L\ensuremath{_\mathrm{\odot}}), LIRGs, and ULIRGs. Comparing \,L\ensuremath{_\mathrm{IR}}\,$<10^{11}$\,\,L\ensuremath{_\mathrm{\odot}}\ and \,L\ensuremath{_\mathrm{IR}}\,$>10^{12}$\,\,L\ensuremath{_\mathrm{\odot}}, we find the drop of the PAH power as a whole (black line) shows no difference in the two geometries, both decrease by 61\%. In general, the PAH band luminosities of individual bands also behave similarly in both geometries except for the 3.3\,\micron\ PAH. While the drop of the L\ensuremath{_\mathrm{PAH\,3.3}}/\,L\ensuremath{_\mathrm{IR}}\ is greater in the obscured continuum geometry (60\%) than in the mixed geometry (42\%), both drop less than L\ensuremath{_\mathrm{\Sigma\,PAH}}/\,L\ensuremath{_\mathrm{IR}}\ (61\%). This suggests that the smallest PAHs are certainly not subjected to enhanced suppression, even though they are often considered the most vulnerable PAH population \citep{Allain1996}.

While all the PAH bands exhibit deficits as L$_\mathrm{PAH}$/\,L\ensuremath{_\mathrm{IR}}\ declines with increasing \,L\ensuremath{_\mathrm{IR}}, they show diverse trends in L$_\mathrm{PAH}$/L\ensuremath{_\mathrm{\Sigma\,PAH}}, as shown in Fig.~\ref{fig:pah_trend}(c) and (d). Both the 3.3 and 8.6\,\micron\ PAH exhibit noticeable increases in fractional luminosity (by 40\% and 20\%, respectively), whereas the 17\,\micron\ PAH decreases significantly by over 30\%. The increase in the fractional 8.6\,\micron\ PAH is due to the fact that galaxies with higher infrared luminosity are likely exhibiting stronger silicate absorption, and the depth of such feature may affect the measurement of nearby 8.6\,\micron\ PAH. The increase in fractional 3.3\,\micron\ PAH power is rather interesting since it implies the size distribution shifts towards smaller end as \,L\ensuremath{_\mathrm{IR}}\ increases. Even in the obscured continuum geometry 3.3\,\micron\ PAH fractional power decreases only insignificantly while the 17\,\micron\ PAH drops by 34\%. Note that the true, attenuation independent behavior is likely bracketed by the results in (c) and (d). Hence, we conclude there is no statistically significant downward trend in L\ensuremath{_\mathrm{PAH\,3.3}}/L\ensuremath{_\mathrm{\Sigma\,PAH}}\ as \,L\ensuremath{_\mathrm{IR}}\ increases, suggesting rather non-intuitively that small PAH populations have an advantage in surviving in environments with strong global radiation fields, or are otherwise able to emit more effectively in exciting radiation (see \S~\ref{subsec:size_distribution} for details).

With our sample selection criteria using the AGN classification scheme of \citet[][see \S~\ref{subsec:bright_PAH}]{Stern2005}, a modest number of galaxies harboring deeply-embedded AGNs might remain in the sample \citep[e.g.][]{Petric2011}, especially at \,L\ensuremath{_\mathrm{IR}}\,$>5 \times10^{12}$\,L\ensuremath{_\mathrm{\odot}} \citep{Nardini2010, Alonso-Herrero2012}. Even though 20\% of our sample are ULIRGs (23 galaxies), none have \,L\ensuremath{_\mathrm{IR}}\,$>5 \times10^{12}$\,L\ensuremath{_\mathrm{\odot}}, with most below $3\times10^{12}$\,L\ensuremath{_\mathrm{\odot}}, and only 2 galaxies in the range $=3$--$5\times10^{12}$\,L\ensuremath{_\mathrm{\odot}}. We note that the possible inclusion of AGN in our sample could contribute to the drop of L$_\mathrm{PAH}$/\,L\ensuremath{_\mathrm{IR}}\ in Fig.~\ref{fig:pah_trend}, but does not alter conclusions drawn from the absence of a decreasing trend of L\ensuremath{_\mathrm{PAH\,3.3}}/L\ensuremath{_\mathrm{\Sigma\,PAH}}\ with increasing \,L\ensuremath{_\mathrm{IR}}\ (Fig.~\ref{fig:pah_trend} bottom panels).

\begin{deluxetable}{ccccc}
\tabletypesize{\footnotesize}
\tablecolumns{4}
\tablewidth{0pt}
\tablecaption{Best-fit Photometric Calibration Parameters for 3.3\,\micron\ PAH using \textit{JWST}/NIRCam Medium Width Filters}
\label{tab:JWST_params}
\tablehead{
\colhead{z} &  \colhead{S} &  \colhead{$\mathrm{\sigma_S}$} & \colhead{$W_{300}$} & \colhead{$W_{360}$} \\
}
\startdata
0.00 & 10.78 & 0.44 & 0.35 & 0.65 \\
0.01 & 11.56 & 0.42 & 0.35 & 0.65 \\
0.02 & 12.88 & 0.52 & 0.34 & 0.66 \\
0.03 & 15.55 & 0.73 & 0.30 & 0.70 \\
\enddata
\vspace{0.cm}
\end{deluxetable}

\subsection{\textit{JWST}\ Photometric 3.3\,\micron\ PAH Strengths} \label{sec:JWST_filter}
\textit{JWST}\ will provide routine photometric and spectroscopic access to the 3.3\,\micron\ PAH feature, which has not been as extensively studied as the other brighter PAH emission features found at longer wavelengths. The three \textit{JWST}/NIRCam medium width filters F300M, F335M, and F360M\footnote[3]{https://jwst-docs.stsci.edu/near-infrared-camera/nircam-instrumentation/nircam-filters} are customized for detecting the 3.3\,\micron\ PAH. To investigate the reliability of inferring 3.3\,\micron\ PAH power from \textit{JWST}/NIRCam photometry, in Fig.~\ref{fig:JWST_filter} we compare synthetic photometry derived from our full spectra with the 3.3\,\micron\ PAH strength recovered by \textsc{Pahfit}\ in the obscured continuum geometry (i.e. without attenuation correction). The synthetic photometry is calculated using the response curves of the respective \textit{JWST}/NIRCam medium filters with the rest-frame 1C spectrum. The x-axis represents the result of the synthetic photometry of the 3.3\,\micron\ PAH band, derived from a subtraction of the continuum contribution to F335M. The continuum is represented by a weighted combination of F300M and F360M, with weighting coefficients optimized to provide the minimum scatter against the 3.3\,\micron\ PAH band power measured in \textsc{Pahfit}. We only provide calibration up to z=0.03 since this calibration does not hold for higher redshifts as the 3.3\,\micron\ band falls out of the F335M filter. The photometric calibration of the 3.3\,\micron\ PAH feature using the photometric bands with respect to different redshifts can be parameterized by:

\begin{equation}
    \bigg(\frac{f_\mathrm{PAH \ 3.3}}{10^{-14} \ \mathrm{erg \ s^{-1} \ cm^{-2}}}\bigg) = (\mathrm{S} \ \pm \ \sigma_\mathrm{S}) \cdot x,
\end{equation}
where
\begin{equation}
    x = \mathrm{F335M} - (\mathrm{W_{300}}\times\mathrm{F300M} + \mathrm{W_{360}}\times\mathrm{F360M}) \ \mathrm{[mJy]},
\end{equation}
$\mathrm{S}$ is the slope of the relation with an uncertainty of $\mathrm{\sigma_{S}}$, and $\mathrm{W_{300}}$ together with $\mathrm{W_{360}}$ are the weightings for the F300M and F360M filters. The corresponding best-fit parameters in different redshifts are listed in Table~\ref{tab:JWST_params}. Note that this calibration is derived from star-forming galaxies with moderate obscuration; heavily embedded galaxies may not be applicable. The reason for the local continuum at around 3.3\,\micron\ to have a slightly lower weighting toward F300M may be due to the fact that F300M is subjected to the water ice absorption at 3.05\,\micron\ that can induce a greater dispersion.

\begin{figure}
 	\includegraphics[width=\columnwidth]{JWST_filter_infer_pah33}
    \caption{The comparison between full spectroscopic and \textit{JWST}/NIRCam synthetic photometry of the 3.3\,\micron\ PAH feature. The response curves of the filters F300M, F335M, and F360M are inset. The x-axis is the optimized formula to narrow the dispersion at z=0 (see Table~\ref{tab:JWST_params}).
    \label{fig:JWST_filter}}
\end{figure}

Because the 3.3\,\micron\ PAH has the shortest wavelength among all the PAH bands in the MIR, it is the final band to fall out of \textit{JWST}\ NIRCam or MIRI filters as redshift increases into the distant Universe, providing an advantage for the detection of resolved ISM features in galaxies up to z\,$\sim$\,7. Here we investigate the contribution of the 3.3\,\micron\ PAH in different NIRCam filters with increasing redshift. In Fig.~\ref{fig:JWST_NIRCam_1C}, we show the contribution of the 3.3\,\micron\ PAH in the respective NIRCam filters as a function of redshift. At z=0, 3.3\,\micron\ PAH contributes $\sim$15\% to the F335M filter, and such value evolves with redshift. Each curve shows the median calculated from the star-forming galaxies, which were classified as 1C in Fig.~\ref{fig:sample_selection}(a), with shading region indicating the 25--75\% range. The three filters with longer central wavelengths, i.e. F430M, F460M and F480M, reach higher peaks since they are narrower compared to filters at shorter wavelengths.

\begin{figure}
 	\includegraphics[width=\columnwidth]{JWST_NIRCam_1C.pdf}
    \caption{The fraction of the 3.3\,\micron\ PAH that contributes to the corresponding \textit{JWST}/NIRCam medium filter with respect to the redshift. Each curve is the median of  spectra obtained from galaxies in class 1C. The shaded regions indicate the inner two quartiles of our 1C sample.}
    \label{fig:JWST_NIRCam_1C}
\end{figure}

\section{Discussion} \label{sec:discussion}
With extended wavelength coverage down to 2.5\,\micron, our cross-archival \textit{Spitzer}+\textit{AKARI}\ spectroscopic sample has opened a comprehensive view of the mid-infrared resolved dust emission features in star-forming galaxies. The features at 3.3--3.4\,\micron, although not major energetic contributors, serve as an important diagnostic of the smallest dust grains and large aromatic and aliphatic molecules with 10's of carbon atoms. And yet as the shortest wavelength resolved dust features in the MIR, they suffer most significantly from the effects of dust attenuation, which can impact their interpretation. Here we discuss the effects of dust emission geometry on the intrepretation of 3.3\,\micron\ PAH emission, and then consider the 3.3 and 17\,\micron\ PAH features together as wide dynamic range tracers of opposite ends of the PAH size distribution. The relatively bright 3.3\,\micron\ PAH emission together with the absence of 17\,\micron\ PAH in II~Zw~40 leads to a discussion of possible mechanisms that conduce the stability of small dust grains under intense radiation fields.

\subsection{Impact of the Geometry: Mixed and Obscured Continuum} \label{subsec:geometry_disc}
Among all PAH bands, even considering the localized silicate absorption band, the 3.3\,\micron\ feature suffers the most from attenuation. This is due both to the steep rise of continuum opacity in the extinction curve towards $\lambda$\,$\sim$\,3\,\micron\ as well as 3.05\,\micron\ H$_{2}$O ice absorption. Even at the limited attenuations of $\tau_\mathrm{mixed} \lesssim$ 3.0 of our PAH-bright sample, the interpretation of 3.3\,\micron\ PAH behavior can depend on assumptions of emission/attenuation geometry, and accordingly we have explored the bracketing cases of \emph{fully mixed} and \emph{obscured continuum} geometries.

For our adopted silicate attenuation profile, the optical depth at 3.3\,\micron\ is $\sim$\,70\% of the peak of the primary silicate band at 9.7\,\micron. While H$_{2}$O ice also has an impact at 3.3\,\micron, the dominant attenuation component is the smooth continuum opacity associated with and constrained by silicate absorption. For the sample average $\tau_\mathrm{H_{2}O}$/$\tau_\mathrm{Si}$=0.25 in the mixed geometry case, H$_{2}$O ice only contributes an average $\sim$\,18\% of the total attenuation at 3.3\,\micron. Among all PAH bands, 3.3\,\micron\ PAH has the largest attenuation correction factor, $\sim$\,1.5, for an extinction curve characterized by typical parameters \ensuremath{\tau_\mathrm{Si}}=1 and \ensuremath{\tau_\mathrm{H_{2}O}}=0.2, followed by a factor of 1.2 for both the 8.6 and 11.3\,\micron\ bands that flank the silicate feature.

Without prior knowledge of the geometry, it is challenging to constrain the \emph{intrinsic} power of 3.3\,\micron\ PAH emission, especially for highly obscured galaxies. Even for high metallicity systems powered entirely by star formation, SFRs inferred from the 3.3\,\micron\ PAH can differ by a factor of $\sim$\,3 compared to SFR$_\mathrm{Ne}$ in galaxies with $\tau_\mathrm{Si}$+$\tau_\mathrm{H_{2}O}\!\sim\!4$ (see Fig.~\ref{fig:SFR_calibration}). Statistically, only at $\tau _\mathrm{Si}$+$\tau_\mathrm{H_{2}O}$ $<$ 0.9 is there no significant difference in the ratios of 3.3\,\micron\ PAH to neon-based SFR between the two geometries (see \S~\ref{subsubsec:SFR_geometry}). Another extreme beyond the mixed geometry is the screen geometry, which assumes attenuation arise from a monolithic foreground layer of dust. A similar optical depth in screen geometry results in more pronounced attenuation than in a mixed geometry. The ``convergent" attenuation $\tau$=0.90 in the mixed geometry is equivalent to $\tau$=0.42 in a screen geometry, and $\tau_\mathrm{screen}/\tau_\mathrm{mixed} \simeq 0.47$ is approximately constant within our attenuation range. It is important to note that \emph{the effective correction for this attenuation is nearly identical} between mixed and screen geometries, since both apply to all emission components (lines, continuum, and resolved dust emission bands).

To conclude, applying both mixed and obscured continuum geometries permits constraining the impact of attenuation on our results. L\ensuremath{_\mathrm{PAH\,3.3}}/\,L\ensuremath{_\mathrm{IR}}\ is most sensitive to the obscuration geometry, increasing by up to $\sim$\,5$\times$ for the most embedded sources in the mixed geometry. For L\ensuremath{_\mathrm{PAH\,3.3}}/L\ensuremath{_\mathrm{\Sigma\,PAH}}, since attenuation impacts all PAH bands to differing degrees, the difference is not as pronounced ---  at most a factor of $\sim$3$\times$. Similarly, as we compare L\ensuremath{_\mathrm{PAH\,3.3}}\ to SFR$_\mathrm{Ne}$, the strength of the neon lines are also subjected to the applied geometry, which results in a difference up to 3$\times$ at high obscuration in our sample. The aliphatic-to-aromatic ratio (3.4/3.3) is not significantly affected by the applied geometry due to wavelength proximity, with variation $\lesssim 20$\%.

\subsection{Impact of the extinction curve}
The impact of the selected extinction curve $A_\lambda$ on our study is minor compared to the geometry. In particular, the \emph{obscured continuum} geometry yield feature strengths which are largely independent of $A_\lambda$ (which only affects the underlying continuum). To investigate the impact on our \emph{mixed} geometry results, we compared PAH band ratios derived using the PAHFIT-default Galactic center extinction curve of \citet{Kemper2004}, to the newly published extinction curve of \citet{Hensley2020} measured from Galactic blue hypergiant Cyg OB2-12. The main difference between the two curves normalized at 9.7\,\micron\ is that \citeauthor{Hensley2020} has a flatter (grayer) profile between 5-8\,\micron\ with broader silicate features at 9.7 and 18\,\micron. Despite these differences, A$_{3.3\,\micron}$/A$_{9.7\,\micron}$ is nearly identical between the two profiles.

For the median galaxy in our sample with \ensuremath{\tau_\mathrm{Si}}=1.10 and \ensuremath{\tau_\mathrm{H_{2}O}}=0.20, applying the extinction curve of \citet{Hensley2020} reduces the fractional 3.3\,\micron\ PAH power by $\sim$10\% due to the fact that the extinction-corrected PAH strengths between 5-8\,\micron\ can be enhanced due to the flatter profile in that range. Our provided 3.3\,\micron\ PAH star formation rate (Eq.~\ref{eqn:SFR}) is not affected by the choice of $A_\lambda$, since it is calibrated against galaxies with limited obscuration (\ensuremath{\tau_\mathrm{Si}}\ + \ensuremath{\tau_\mathrm{H_{2}O}}\ $<$ 0.9). Using the extinction curve of \citeauthor{Hensley2020} does not mitigate the offset between measured and modeled band ratios, as the ratio of 17/11.3 remains unchanged ($<$6\%) and the ratio of 3.3/11.3 drops only by $\sim$9\%. In the worst case of IRAS~F06076-2139 --- the most ``obscured" galaxy in our bright-PAH sample, with \ensuremath{\tau_\mathrm{Si}}=4.21 and \ensuremath{\tau_\mathrm{H_{2}O}}=0.23 in mixed geometry --- the drop of 3.3/11.3 is at most $\sim$30\%. Such a drop actually further \emph{increases} the discrepancies in Fig.~\ref{fig:pah33andpah17} between measurement and model. To gain a comprehensive understanding of how different extinction curves impact recovered PAH bands, further investigation on more embedded sources is needed.

\begin{figure}
 	\includegraphics[width=\columnwidth]{PAH33and17}
    \caption{The change of the 3.3 (circle) and 17 (triangle)\,\micron\ PAH in the L$\mathrm{_{PAH}}$/\,L\ensuremath{_\mathrm{IR}}\ and L$\mathrm{_{PAH}}$/L\ensuremath{_\mathrm{\Sigma\,PAH}}\ space as infrared luminosity increases (directions of the arrows). Solid (empty) point indicates measurements in the mixed (obscured continuum) geometry. Apparently 17\,\micron\ PAH band decreases in both geometries. While the behavior of 3.3\,\micron\ PAH track depends on the adopted geometry, there is \emph{no} sign of decrease of the fractional 3.3\,\micron\ PAH power as \,L\ensuremath{_\mathrm{IR}}\ increases (see \S~\ref{subsec:size_distribution}).
    }
    \label{fig:PAH33and17}
\end{figure}

\subsection{The Size Distribution of PAHs} \label{subsec:size_distribution}
Among all the PAH emission bands, the 3.3\,\micron\ band and 17\,\micron\ complex trace PAH molecules at opposite ends of the PAH size distribution, probing the smallest and largest grains, respectively \citep{Schutte1993, Draine2007}. The variation between the two bands therefore allows us to probe the overall size shift in different interstellar environments across a wide range of sizes, from 10's to 1000's of carbon atoms. In Fig.~\ref{fig:pah_trend}, we demonstrated the behaviors of the dominant PAH emission features with increasing \,L\ensuremath{_\mathrm{IR}}. Here in Fig.~\ref{fig:PAH33and17}, we specifically compare the changing contributions of the 3.3\,\micron\ (circle) and 17\,\micron\ (triangle) PAH bands for the entire sample. Both L\ensuremath{_\mathrm{PAH\,3.3}}/\,L\ensuremath{_\mathrm{IR}}\ and L\ensuremath{_\mathrm{PAH\,17}}/\,L\ensuremath{_\mathrm{IR}}\ decrease as \,L\ensuremath{_\mathrm{IR}}\ increases (indicating by arrows all pointing left), regardless of the geometry. While the ratio of L\ensuremath{_\mathrm{PAH\,17}}\ to \,L\ensuremath{_\mathrm{IR}}\ and L\ensuremath{_\mathrm{\Sigma\,PAH}}\ decreases in a similar fashion in both geometries, the behavior of 3.3\,\micron\ band differs due to the larger impact of attenuation. In the fully mixed scheme, as \,L\ensuremath{_\mathrm{IR}}\ increases, the trend of 3.3\,\micron\ PAH emission tracks towards the upper left, suggesting an increase of the fractional 3.3\,\micron\ PAH abundance. In contrast, the trend moves broadly towards the left in the obscured continuum scheme, with no sign of significant change in the fractional 3.3\,\micron\ PAH power. A precise projection of the trend for an individual galaxy requires a better handle of the geometry --- a challenging task requiring a sophisticated radiative transfer model that is beyond the scope of this study. Nevertheless, these two bracketing trends allow us to narrow the possibilities ---  L\ensuremath{_\mathrm{PAH\,3.3}}/L\ensuremath{_\mathrm{\Sigma\,PAH}}\ either increases or is relatively unchanged as \,L\ensuremath{_\mathrm{IR}}\ increases.

The fact that the fractional 17\,\micron\ PAH decreases together with the increase of the fractional 3.3\,\micron\ PAH supports theories of top-down fragmentation \citep[e.g.][]{Jones2016}. Enhanced photo-destruction of PAHs may change the size distribution \citep{Allain1996}, leading to the breakdown of larger PAHs into smaller PAH segments. A scenario with the photochemical erosion process is also supported by the decline of the ratio of 3.4/3.3 with the increasing \,L\ensuremath{_\mathrm{IR}}\ (\S~\ref{sec:aliphatic}).

Another potential cause for shifting PAH power from longer to shorter wavelengths is the enhanced radiation field that makes larger PAHs ``act'' like smaller PAHs. In a high UV radiation density environment, the mean time between photon absorptions become shorter than the radiative cooling time for PAHs, allowing molecules to reach higher energy level by encountering multiple photons during the relaxation process, which in turn followed by short-wavelength emission that resembles emission of small PAHs \citep{Draine2007}. A detailed model of the PAH size distribution with respect to changing starlight fields will shed light on whether radiation could be the main driver of the enhancement of the fractional 3.3\,\micron\ PAH in environments with high infrared luminosity or at low metallicity.

\subsection{The Stability of Small Grains} \label{subsec:grain_stability}
With a strong 3.3\,\micron\ emission, together with non-detection of the 17\,\micron\ PAH, II~Zw~40 casts considerable doubt on the suggestion that only larger PAHs are able to survive in harsh low-metallicity environments \citep{Madden2006, Gordon2008, Hunt2010}. \citet{Sandstrom2012} proposed an alternative scenario suggesting a size distribution shifted towards smaller PAHs in low-metallicity environment is due to the intrinsically enhanced formation of small grains rather than dust processing. That said, bottom-up formation process seemed be favored in low-metallicity environments as opposed to the top-down fragmentation process. The fact that strong 3.3\,\micron\ PAH is detected in II~Zw~40 demonstrates that small grains are able to survive even in the harsh radiation environments of metal-poor starburst galaxies.

Despite the fact that in general small PAHs are readily destroyed by UV photons, if the absorbed energy can be channeled through a relaxation process other than dissociation, photo-destruction can be avoided. PAH molecules have three relaxation channels: dissociation, IR vibrational emission, and recurrent fluorescence (RF). Dissociation becomes efficient only when excitation energy is higher than a certain level (depending on the size of the molecule). At lower energies, IR emission and recurrent fluorescence are the two main relaxation processes. The relaxation timescale for IR emission is on the order of $\sim$s, whereas recurrent fluorescence is a relatively fast radiative cooling process with relaxation time on the order of $\sim$ms \citep{Bernard2017}. According to \citet{Leger1988}, molecule consists of fewer carbon atoms produces RF photons more efficiently. Hence, it is worth noting that RF can be a potentially powerful mechanism to significantly increases the ability for small particles to survive under high UV radiation densities. Experimentally, such process has been confirmed in the laboratory that RF can be a main radiative cooling mechanism for PAH molecules by observing it through anthracene cations trapped in a compact electrostatic storage ring, the Mini-Ring \citep{Martin2013, Martin2015, Bernard2017}. Observationally, it has been proposed to conduce to the survival rate of the carriers of the extended red emission and diffuse interstellar bands \citep{Lai2017, Lai2020}. Therefore, RF may also play an important role in maintaining the abundance of small PAHs in the ISM, especially in environments with strong and high energy prevailing UV radiation fields such as II~Zw~40.

\section{Summary}
\label{sec:summary}
We present a sample of 378 galaxies drawn from the \textit{AKARI}-\textit{Spitzer}\ Extragalactic Spectral Survey (ASESS), an exhaustive cross-archival comparison between the \textit{AKARI}/IRC dataset and the \textit{Spitzer}/IRS IDEOS extragalactic catalog, with combined wavelength coverage from 2.5--38\,\micron. A set of composite template spectra adopted from the classification scheme of \citet{Spoon2007} is released. With extended wavelength coverage down to 2.5\,\micron, we place the 3.3\,\micron\ PAH feature in context of the complete MIR emission spectrum, and explore it as a powerful diagnostic of the smallest PAH molecules in a diverse range of galaxies. A subset of 113 PAH-bright galaxies with low-to-moderate dust attenuation were investigated using an extended version of the spectral decomposition tool \textsc{Pahfit}. This has enabled, for the first time, modeling \textit{all} PAH emission bands in the MIR simultaneously. Our major findings are summarized as follows:

\smallskip

\begin{itemize}
\item An extended version of \textsc{Pahfit}\ with new features, including the aromatic and aliphatic bands, recombination lines, and several absorption components, found to be highly successful in modelling the complete MIR spectrum in a diverse sample of galaxies. We find the 5.27\,\micron\ PAH feature can further be divided into two individual bands centered at 5.24 and 5.33\,\micron.

\item A mild correlation exists between radiation heating intensity and the aliphatic-to-aromatic ratio (3.4/3.3), which suggests that aliphatics are more prone to photo-destruction under intense radiation environments. We also find that a broad ``plateau'' feature centered at 3.47\,\micron\ is required in the model. This plateau typically carries 20\% of the 3.3\,\micron\ PAH power, and is likely aromatic in nature, as it correlates better with 3.3\,\micron\ PAH than with 3.4\,\micron\ aliphatic emission.

\item As the shortest of the PAH vibrational emission features, we find the 3.3\,\micron\ PAH is more susceptible to the effects of attenuation than any other PAH band. To bracket the impact of attenuation, two geometrical schemes are adopted: \emph{fully mixed geometry}, in which dust and emission components are mixed and \emph{obscured continuum geometry}, in which line and PAH emission lie atop an attenuated underlying continuum. We find that diagnostics including 3.3\,\micron\ PAH emission diverge between these two geometries at moderate optical depths \ensuremath{\tau_\mathrm{Si}}+\ensuremath{\tau_\mathrm{H_{2}O}}\,$\gtrsim$\,1 (equivalent to 0.5 in screen geometry).

\item The bolometric fraction L\ensuremath{_\mathrm{PAH\,3.3}}/\,L\ensuremath{_\mathrm{IR}}\ is $\sim$\,0.1\%, although this varies between the two geometries by a median $\sim$\,1.5$\times$, and this discrepancy increases to $\sim$\,5$\times$ in the most obscured galaxies in our sample. The fractional 3.3\,\micron\ PAH power depends less sensitively on attenuation geometry, and is typically L\ensuremath{_\mathrm{PAH\,3.3}}/L\ensuremath{_\mathrm{\Sigma\,PAH}}\,$\simeq$\,1.5--3\%.

\item Ratios among the main non-ionized PAH bands at 3.3, 11.3, and 17\,\micron\ follow well the expected tracks for individual neutral PAH molecules, but the ratio 3.3/11.3 is typically $2\times$ smaller than the nominal value adopted in the dust model of \citet{Draine2007}, indendent of the adopted extinction curve.  The ratio of 3.3\,\micron\ to 11.3\,\micron\ power rises mildly in galaxies with higher infrared luminosity, independent of the adopted geometry.

\item We provide a new star formation rate calibration based on 3.3\,\micron\ PAH luminosity, constructed from a carefully selected sample of low-attenuation star forming galaxies. However, we find that obscuration plays an important role in this use, with a factor of $\sim$\,3 divergence between different geometries at high obscuration, even in this pre-selected sample. Metallicity plays an even larger role, with SFR biased low by $\gtrsim8\times$ at $1/4\,Z_\sun$.

\item Although we confirm the significant drop in power of all PAH features as a fraction of bolometric infrared luminosity as \,L\ensuremath{_\mathrm{IR}}\ increases from 10$^{9}$ to 10$^{12}$\,L\ensuremath{_\mathrm{\odot}}, surprisingly, the 3.3\,\micron\ feature, which traces the smallest and presumably most easily destroyed molecules, constitutes a flat or even \emph{rising} portion of total PAH emission, depending on the geometry adopted. In contrast, the 17\,\micron\ feature, which traces the largest PAH grains with thousands of carbon atoms, drops most acutely at the highest infrared luminosities, independent of the assumed geometry.

\item We observed highly suppressed PAH emission in low metallicity galaxy II~Zw~40, with L\ensuremath{_\mathrm{\Sigma\,PAH}}/\,L\ensuremath{_\mathrm{IR}}\ less than other sampled galaxies by a factor of $\sim$\,10. Surprisingly, however, 3.3\,\micron\ power \emph{increases} relative to other PAH bands, with a complete absence of the 17\,\micron\ PAH. This strongly suggests that caution is needed when using PAH bands as SFR indicators in moderate to low-metallicity environments, and that the smallest grains are more robust against photo-destruction in harsh radiation environments than previously believed.

\item  We suggest that the likely mechanisms leading to strong fractional 3.3\,\micron\ PAH power are either enhanced stability of small PAH molecules from efficient relaxation processes (e.g. recurrent fluorescence) or the migration of PAH power to shorter wavelengths under environments with intense and/or hard radiation fields.

\item Although 3.3\,\micron\ PAH emission contributes just 15--25\% of the power in JWST/NIRCam medium band observations out to $z\,\sim\,0.5$, using the flanking bands for continuum subtraction permits recovery of the 3.3\,\micron\ PAH power up to z=0.03 within 5\%.

\end{itemize}

\textit{JWST}\ will soon re-open the MIR window with more than an order of magnitude increase in sensitivity and spatial resolution, and the all-sky IR spectroscopic survey \textit{SPHEREx} will offer unprecedented access to more than 100 million galaxies in the 3\,\micron\ wavelength regime. Together with comprehensive modeling of the complete suite of PAH emission, the diagnostic potential of these powerful bands will be fully realized, opening new windows into the study of star formation, radiation environment, metallicity, AGN processing, dust physics, and the physical properties of the ISM on sub-galactic scales and across much of cosmic time.
\\
\acknowledgments
We thank the anonymous referee for the insightful comments that helped to improve this manuscript. This research is based on observations with \textit{AKARI}, a JAXA project with the participation of ESA. This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration. TSYL gratefully acknowledges support from the University of Toledo Graduate Fellowship and the NASA ERS Grant and the host of ISAS, JAXA. HWWS would like to acknowledge support for the IDEOS project from NASA grant NNX16AF26G.
\\
\vspace{5mm}
\facilities{AKARI, Spitzer}

\software{Astropy \citep{Astropy2013, Astropy2018},
          Bokeh \citep{Bokeh2018},
          Matplotlib \citep{Hunter2007},
          Numpy \citep{VanderWalt2011},
          \textsc{Pahfit}\ \citep{Smith2007},
          SciPy \citep{Virtanen2020}
          }