# M94 as a Unique Testbed for Black Hole Mass Estimates and AGN Activity at Low Luminosities

Authors: Anca Constantin; Anil C. Seth
Journal: Advances in Astronomy
 (2012)

---

## Abstract

We discuss the peculiar nature of the nucleus of M94 (NGC 4736) in the context of new measurements of the broadHα emission from HST-STIS observations. We show that this component is unambiguously associated with the high-resolution X-ray, radio, and variable UV sources detected at the optical nucleus of this galaxy. These multiwavelength observations suggest that NGC 4736 is one of the least luminous broad-line (type 1) LINERs, with Lbol=2.5×1040 erg s−1. This LINER galaxy has also possibly the least luminous broad-line region known (LHα=2.2×1037 erg s−1). We compare black hole mass estimates of this system to the recently measured ~7×106M⨀ dynamical black hole mass measurement. The fundamental plane and M-σ* relationship roughly agree with the measured black hole mass, while other accretion-based estimates (the M-FWHM(Hα) relation, empirical correlation of BH mass with high-ionization mid-IR emission lines, and the X-ray excess variance) provide much lower estimates (~105M⨀). An energy budget test shows that the AGN in this system may be deficient in ionizing radiation relative to the observed emission-line activity. This deficiency may result from source variability or the superposition of multiple sources including supernovae.


---

## Body

## 1. Introduction: LINERs and M94

Most high-mass galaxies are known to host massive black holes (BH), some passively lurking in their centers while others are actively accreting surrounding material (e.g., [1] and references therein). The mechanism that causes a BH’s activity to turn on and off is still largely unknown. Understanding the structure of the active galactic nuclei (AGN) at their lowest luminosities is crucial to determining the physical and possibly evolutionary links between the most luminous galaxy centers and the passive ones. However, at low luminosities it is difficult to disentangle the various emission mechanisms that could be concurrently present in galaxy centers. As a consequence, the dominant power source of a large majority of actively line-emitting galaxies remains ambiguous (e.g., [2] for a review).Diagnostic diagrams (e.g., [3–6]) are relatively successful in separating out bona fide accretion sources (Seyferts) from nuclei whose emission-line activity is mainly powered by young, hot stars (H ii galaxies), based on emission line ratios. At least 50% of the strong line emitters fall easily onto the H ii  locus; however, only less than 10% are of the Seyfert type [5, 7, 8]. A large fraction of the objects “in between” these two categories, that exhibit relatively low levels of ionization (i.e., low values of [O III]/Hβ), maintain reasonably strong forbidden line activity (i.e., high values of [N II]/Hα) and are classified as low-ionization nuclear emission regions (LINERs). The reminder are usually called transition objects (Ts). Whether Ts and Ls are powered, at least partly, by accreting BHs, and thus could be called AGN, is a matter of continuous debate [2].There are some typical emission characteristics that are considered to be particularly good indications that accretion onto a massive BH is an important, if not the dominant source of ionization in some Ls, and possibly Ts. The detection of broadHα emission, regardless of its strength or luminosity, is generally considered to be the clue to AGN emission. Some LINERs (and maybe Ts as well) exhibit these features; however, the majority of these systems show only narrow emission, which could be generated by shocks, poststarbursts, or other processes unrelated with accretion. Observations outside of the optical wavelengths often reveal AGN signatures in ambiguous and even starburst galaxy nuclei. X-rays are particularly good tracers of accretion; however, they are not efficient in distinguishing AGN at LX≲1042 erg s-1, where contamination by X-ray binaries can be significant. X-rays are also unlikely to detect heavily absorbed AGN (i.e., Compton thick; NH>1.5×1024cm-2). On the other hand, mid-IR high-ionization emission lines like [Ne V] λ14.32 μm, 24.32 μm (97.1 eV) appear to be a trustworthy indicator of AGN activity [9–11] due to the extreme conditions (i.e., very hard ionizing radiation) required to produce them; because these features have considerably lower optical depth, their detection can also reveal Compton-thick AGN ([12], e.g., NGC 1068). This technique has now been applied to reveal new and large numbers of optically unidentified AGN [13, 14], providing thus sensitive improvements on previous AGN censuses. These studies, along with X-ray and radio studies of nearby galaxies (e.g., [15–17]). suggest that a majority of LINERs and a large fraction of transition galaxies might in fact host accreting black holes. The presence of an accreting black hole does not however guarantee that the accretion power is the dominant source of ionization of those galaxy nuclei. A more recent assessment of the energy budget of LINERs by Eracleous et al. [18] argues that in 85% of LINERs the AGN ionizing photons are not sufficient for producing the observed nebular emission, and thus other power sources are likely to dominate.Some new potentially powerful insights into the excitation mechanism of the low-luminosity AGN (LLAGN), and in particular the ambiguous sources, come from recent studies of large statistical samples of nearby galaxy nuclei, which reveal a potential Hii→S→T→LINER→Passive Galaxies evolutionary sequence in the process of BH growth within galaxies [19, 20]. This sequence traces trends in (1) increasing host halo mass, (2) increasing environmental density, (3) increasing central BH mass and host stellar mass, (4) decreasing BH accretion rate, (5) aging of the stellar population associated with their nuclei, and (6) decreasing in the amount of dust obscuration, which might translate into a decrease in the amount of material available for star forming or accretion. In this picture, Seyferts and Ts are transition phases between the initial onset of accretion, usually swamped by the star-forming gas and associated dust, which is seen optically as an H ii system, and the final phase of accretion observed as LINERs of already massive BHs. While this idea is supported by various other independent observational studies of low-luminosity AGN and starburst galaxies [21, 22], along with state-of-the-art hydrodynamical models for the life cycles of the most luminous AGN (e.g., [23, 24]) it is very probable that not all sources fit into this scenario. It is also important that this evolutionary sequence idea is tested on samples that span a narrow distribution in Hubble types; because most of the H iis and Ts are relatively late-type disk galaxies that, likely, never experienced a (recent) major merger, it is possible that the trigger of such a sequence is different from the merger that initiates a similar life cycle at high luminosities. Nevertheless, investigating the challenges certain objects bring to this idea is useful for identifying and quantifying the caveats associated with this sequence. These objects may also be at interesting stages in their galaxy evolution.NGC 4736 (or M94, UGC 7996) is a captivating example of an ambiguous galaxy nucleus, which poses challenges to the general understanding of AGN phenomena, including the above-mentioned sequence. This object is one of the closest (d=4.3 Mpc; [25]) nearly face-on spiral, with a SAab Hubble Type. Its proximity enables study of details that would be unobservable in more distant systems. Its nucleus has a low-luminosity LINER spectrum, but has been also included in catalogs of transition objects [26] or Seyfert 2s (e.g., [27]). The AGN nature of this object has been constantly debated, an aging starburst being a compelling alternative [28–31]. The galaxy presents a ring of H ii regions at a radius of ~50′′, red arcs at ~15′′, a high-surface-brightness nuclear region, and high far-infrared bulge emission [32, 33]. Its intricate structure of off-nuclear compact source detections in X-ray [30], radio [31], and UV [29], that do not necessarily match with each other, certainly increases the ambiguity associated with the nature of the main nuclear ionization mechanism. A common implied scenario in all of these studies is that this system is probably in the final stages of a merger.We reexamine here NGC 4736 in the context of additional evidence for its AGN nature, which is the detection of a broadHα component in its nuclear spectrum, as observed by the Hubble Space Telescope (HST) with the Space Telescope Imaging Spectrograph (STIS). We also gather multiwavelength data and show that source of the broad H-alpha emission line is coincident with a compact X-ray and radio source. These observations suggest that NGC 4736 hosts a broad-line region of significantly low luminosity, which makes this object one of the least luminous LINERs with strong evidence for BH accretion. Interestingly, BH mass indicators calibrated on rapidly accreting Seyfert galaxies give highly discrepant mass estimates in this more quiescent system. The ~107M⨀ value given by the M-σ* is two orders of magnitude higher than the values obtained via estimators based on the observed emission (X-ray variability, scaling relations, mid-IR emission) which, although physically independent of each other, give a consistent result of ~105M⨀. We are thus facing the following conundrum: either (a) the standard AGN BH mass indicators do not necessarily apply to sources emitting in this low-luminosity regime as, probably, the emission mechanism is fundamentally different from that associated with higher-luminosity AGN, or (b) the emission signatures do not trace accretion onto the central BH. We discuss possible resolutions of this discrepancy in Section 5, where we propose some rather exotic scenarios.


## 2. Data Compilation and Analysis

NGC 4736 has been quite extensively observed across the whole electromagnetic spectrum. We present in this section the multiwavelength observations of this galaxy nucleus, in connection to new measurements of the broadHα emission detected with HST-STIS.


### 2.1. The BroadHα Emission Observed with HST-STIS

High-resolution optical spectra of the NGC 4736 nucleus were obtained withHST-STIS on July 2002. Data are publicly available (Prop ID 8591) but have not been published. The observations were obtained with the 52′′×0.′′1 aperture oriented at PA=49.65°, with the slit centered along the major axis of the starlight distribution; two cosmic-ray split exposures were obtained, one being slightly shifted in the slit direction. The total combined exposure time is ~4000 s. The G750M grating was set at 6581 Å, with a scale of 0.′′05/pixel, with no binning. We reduced the spectra using IRAF (IRAF is distributed by NOAO, which is operated by AURA Inc., under contract with the National Science Foundation) and the STIS reduction pipeline maintained by the Space Telescope Science Institute [34]. This reduction included image combination and cosmic-ray rejection, flux calibration, and correction of the wavelength to the heliocentric frame. To measure the nuclear nebular line-emission properties we extracted the 1-dimensional aperture spectrum five pixels wide (0.′′25) centered on the continuum peak. The extracted spectrum thus consists of the central emission convolved with the STIS spatial point-spread function (PSF) and sampled over a rectangular aperture of 0.′′25×0.′′1. The measurements span 6295–6867 Å  with a resolution of 0.87 Å (σinst=17km s−1).Figure1 shows the resulting spectrum together with the best matching continuum model, spectral fits of the pure emission-line component, and the associated residuals. The continuum model of the underlying stellar population is obtained via a χ2 minimization of a nonnegative least-squares fit between the observed spectrum and a sum of discrete star bursts of different ages, adopted from the Bruzual and Charlot [35] stellar population synthesis templates, together with dust attenuation modeled as an additional free parameter. The continuum fitting is performed using an adaptation to our data of Christy Tremonti’s code [36]. The pure emission-line spectrum, obtained by subtracting the modeled continuum from the observed spectrum, is fit by a combination of linear continuum and Gaussian components. In addition to narrow emission, the Hα+[NII] feature shows clear evidence for a broader feature. A flux ratio of 1 : 3 was assumed for the [N II] doublet, as dictated by the branching ratio [37]; the [O I] feature and the [N II] and [S II] doublets were assumed to share common velocity (red) shifts and widths. The best-fitting Gaussian parameters were derived via an interactive χ2 minimization, using SPECFIT [38]. Because of the generally low signal to noise of the 2-d spectrum we are not able to test whether our measurements of the broad Hα feature could be corrupted by a possible rotating disk as in the case of M84 (e.g., [39]); the match in the widths of the [S II] and [N II] narrow features (as long as that of the narrow Hα) argue however against such a significant effect.Figure 1
Spectral fits to theHST-STIS spectrum of NGC 4736, based on a 5-pixel extraction aperture, which corresponds to a total emission area of 0.025 arcsec2. Upper panel: The thin continuous line shows the observed data, while the green line shows the best-fit stellar population model obtained by clipping out the hashed emission-line region. Middle panel: The continuum subtracted spectrum, that is, the emission-line component, together with the corresponding spectral fit, that have been shifted by an additive constant to match the initial continuum level, for illustration purposes. The thin continuous lines show the individual Gaussian components (1 per narrow line), the thick grey continuous line represents the broad Hα feature, and the red line is the final fit to the observed spectrum. Lower panel: The residuals after the subtraction of the model fit.Measurements related to the broadHα feature from the STIS measurements are listed in Table 1. The fractional contribution of this broad component to the total flux of the  Hα+[NII]blend is 80%. Compared to other galaxies with broad-line emission in the Palomar survey [40], the width of the broad Hα is typical (1570 km s−1). However, at just 9.1×1037 ergs s−1, the Hα luminosity is lower than any known broad-line sources in the Palomar sample or the Sloan Digital Sky Survey [40, 41]. This includes the famous low luminosity Seyfert 1 galaxy, NGC 4395, which has a luminosity of 1.2×1038 ergs s−1 [42]  and is at a similar distance to NGC 4736. Note also that the broad Hα feature is redshifted by ~140±20 km s−1 relative to the narrow emission, that is locked at the systemic velocity of the object; velocity offsets between the broad Hα and narrow lines are not unusual [43].Table 1
BroadHα emission measurements.
QuantityValueNotesfblend0.80Fraction of broadHα toHα+[NII] blend.fHα0.93Fraction of broadHα to total HαFWHM(Hαbroad)1570±110; inkms-1Δv140±20; inkms-1; broad relative to narrow HαlogF(Hαbroad)−13.38Observed flux (ergs-1cm-2)logL(Hαbroad)37.96Observed luminosity(ergs-1)aaassuming distanced=4.3Mpc.Evidence for broadHα emission has also been presented in a recent PCA tomography study applied to this nucleus [44]. The observations reported in this case come from the Gemini Multi Object Spectrograph (GMOS)-IFU data cube and have been obtained 4 years after the HST-STIS spectrum had been acquired. Thus, the broad Hα line associated with the nucleus of M94 appears to be persistent for at least 4 years. They derive a broad Hα luminosity of 6×1038 ergs s−1, ≳6 times brighter in the GMOS data than the one we detect in the STIS spectrum. The difference could result from true variability in the broad-line luminosity or could be due to differences in aperture size or placement. We note that Steiner et al. [44] suggest the broad-line region is offset from the photocenter of the galaxy by 0.′′15, while our spectral 0.′′25×0.′′1 aperture coincides with the photocenter. Even with their higher Hα luminosity, this source is among the least luminous known broad-line AGN.


### 2.2. Ground-Based Optical Spectroscopy and the Spectral Classification

Ground-based optical spectra of the NGC 4736 nucleus are available from the Palomar spectroscopic survey by Ho et al. [43] and the integrated spectrophotometric survey of Moustakas and Kennicutt [45] conducted with the 2.3 m Bok telescope, as well as from new data we acquired on February 2008 at the MMT Observatory. The MMT spectrum is obtained with the Blue Chanel Spectrograph, the 500 grooves/mm grating used in first order with the 1′′ slit, and covers λλ3800–7000 with 3.6 Å resolution. None of these data show indications of broad Hα emission. This outcome is not surprising given the feature’s flux in the HST spectrum, which would be very difficult to discern in the 1′′ or 2′′ apertures employed in these ground-based observations, which are at least one order of magnitude larger than that used in the HST-STIS observations. Through a 2′′ aperture (~40 pc), NGC 4736’s emission complex Hα+[N II] is generally heavily swamped by the host stellar light.Probably as expected, measurements of NGC 4736’s narrow-line emission are consistent with this system being powered at least partially by a nonstellar source, regardless of the resolution of its observation. Figure2 shows the location of this nucleus in the 3-dimensional diagnostic diagram usually employed in classifying emission-line galaxies (e.g., [3, 4, 6]), for different sets of data mapping ~10 to ~100 pc. As with the HST-STIS observations, the ground-based measurements are performed after the host stellar contribution has been subtracted from the observed spectrum by means of absorption galaxy template fits. Our MMT spectrum offers only an upper limit for the [O III ]/Hβ while the Palomar line flux measurements are at least 50% uncertain. We list the measurements of the nebular emission of this nucleus as observed by Palomar, Bok, MMTO, and HST-STIS in Table 2. Because measurements of the [O III]/Hβ ratio are not available in the HST spectrum, we show the high-resolution measurements in Figure 2 using the [O III]/Hβ from the Palomar catalog; it is readily apparent that the high-resolution data are consistent with a LINER or a Seyfert classification for this object and reveal its type 1 (broad-lined) character.Table 2
Emission linemeasurementsa.
line nameHST-STISPalomarbBokcMMTOd[O I] λ63005.03±0.316.032.8±0.97.0±1.0Hα (narrow)3.00±0.32> 25.1213.5±1.922.0±2.4[N II] λ65835.25±0.2554.0121.5±1.539.3±2.9[S II] λ67162.48±0.1818.598.3±1.115.5±0.3[S II] λ67312.16±0.1816.337.2±1.19.8±0.3aAll fluxes are in units of  10-15ergs-1cm-2 and represent the observed values, not corrected for reddening.bnonphotometric conditions; line ratios are at least 50% uncertain.c2.′′5 slit, 2.3 m telescope; photometric conditions [45].d1′′slit; fluxes are in units of  10-15ergs-1cm-2.Figure 2
NGC 4736 within emission-line diagnostic diagrams. Filled circles reflect measurements based on the Palomar observations. Diamonds represent measurements from Moustakas and Kennicutt Jr. [45], while crosses show measurements from our MMT spectrum. The vertical bands indicate the HST-STIS measurements and the associated errors (with no correction for reddening; [O III] and Hβ are not available in the HST spectra). The solid (black) curves indicate the Kewley et al. [6] classification, while the dotted (blue) lines indicate criteria used by Ho et al. [43]. The background grey points correspond to measurements of SDSS nearby galaxies from Constantin et al. [19].


### 2.3. The Radio, UV, and X-Ray Observations and the Astrometric Coincidence with the BroadHα Emission

Across the electromagnetic spectrum, the center of this galaxy exhibits a complex morphology. There is a plethora of bright radio, UV, and X-ray sources detected in the center of NGC 4736. There is a nuclear compact (15 GHz) radio source measured by Nagar et al. [46] that appears to be associated with that of the brightest of the two close (8.49 GHz) radio sources revealed by Körding et al. [31], and with that of the brightest of the two ultraviolet point sources, which varies on a 10-year time scale [29]. Chandra observations [30] reveal numerous discrete X-ray sources in the inner galaxy; the second brightest X-ray source, X2 (LX,2–10keV=5.9×1038 erg s-1), coincides within the errors with the nucleus position. The off-nuclear radio, UV, and X-ray sources are apparently unrelated to each other.The relation to possible optical counterparts of these observations is not well constrained in the literature. We thus investigated this issue, with a particular interest in the degree to which our newly detected broadHα coincides in position with the multi-λ detections. We used for this purpose a variety of archival images from the HST: WFPC2/PC data in the F555W filter data (PID: 5741 and 10402), HRC data in the F250W and F330W filters (PID: 9454), and NICMOS/NIC3 data in the F160W filter (PID: 9360). All data were downloaded from the HST archive, and images were drizzled together when required. Absolute astrometry was performed on these data as well as the STIS observations (taken with the slit out) using a ground-based V-band image obtained from the Spitzer Infrared Nearby Galaxies Survey (SINGS) ancillary data by Kennicutt et al. [48]. These SINGS observations were aligned to the USNO-B system using ~80 stars, and then astrometry of the HST data was obtained by degrading the resolution of the F555W image to match the SINGS image. All other HST images were then matched with the astrometry-corrected F555W frame to an accuracy of <0.′′05. The absolute error on this astrometry is about 0.′′2 and is dominated by scatter of stars in the SINGS image relative to the USNO-B positions.Figure3 illustrates the result of this data compilation and the corresponding radio, UV, and X-ray source matches. It is clear that, within<0.′′2(4 pc), there is an obvious astrometric match in the nuclear X-ray, UV, optical, radio compact sources, and the newly detected broadHαemission line. The nucleus position as observed by STIS has RA: 12 h 50 m 53 s.20, DEC: 41°07′13.′′40. The off-nuclear X1, that is the brightest compact X-ray source detected in this galaxy nucleus, has no counterpart at other wavelengths. Same is true for the off-nuclear 8.49 GHz detection, that is only 1′′ (~20 pc) away from the nucleus. The off-nuclear UV source is the only one detected in optical light.Figure 3
The mosaic-ed F555W (V-band; combination of 5 WFPC2 images taken at two different epochs) image of the NGC4736 nucleus, with intensity contours overlaid (green). The position of the STIS-nucleus is indicated (by the yellow slit) along with the locations of the 2 most luminous and closest to the nucleus hard X-ray sources, X1 and X2, as detected with Chandra (red squares), along with the position of the 15 GHz VLA radio core ([46], orange circle), and those of the two radio compact sources detected by Körding et al. [31] (white circles). The blue and cyan contours correspond to the F330W and F250W observations from Maoz et al. [29]; note the presence of NGC 4736b, ≈2.′′5 to the north of the nucleus. The absolute astrometry is good to within 0.′′2 (pixel size is 0.′′05 for F555W, and 0.′′025 for F250W and F330W).


## 2.1. The BroadHα Emission Observed with HST-STIS

High-resolution optical spectra of the NGC 4736 nucleus were obtained withHST-STIS on July 2002. Data are publicly available (Prop ID 8591) but have not been published. The observations were obtained with the 52′′×0.′′1 aperture oriented at PA=49.65°, with the slit centered along the major axis of the starlight distribution; two cosmic-ray split exposures were obtained, one being slightly shifted in the slit direction. The total combined exposure time is ~4000 s. The G750M grating was set at 6581 Å, with a scale of 0.′′05/pixel, with no binning. We reduced the spectra using IRAF (IRAF is distributed by NOAO, which is operated by AURA Inc., under contract with the National Science Foundation) and the STIS reduction pipeline maintained by the Space Telescope Science Institute [34]. This reduction included image combination and cosmic-ray rejection, flux calibration, and correction of the wavelength to the heliocentric frame. To measure the nuclear nebular line-emission properties we extracted the 1-dimensional aperture spectrum five pixels wide (0.′′25) centered on the continuum peak. The extracted spectrum thus consists of the central emission convolved with the STIS spatial point-spread function (PSF) and sampled over a rectangular aperture of 0.′′25×0.′′1. The measurements span 6295–6867 Å  with a resolution of 0.87 Å (σinst=17km s−1).Figure1 shows the resulting spectrum together with the best matching continuum model, spectral fits of the pure emission-line component, and the associated residuals. The continuum model of the underlying stellar population is obtained via a χ2 minimization of a nonnegative least-squares fit between the observed spectrum and a sum of discrete star bursts of different ages, adopted from the Bruzual and Charlot [35] stellar population synthesis templates, together with dust attenuation modeled as an additional free parameter. The continuum fitting is performed using an adaptation to our data of Christy Tremonti’s code [36]. The pure emission-line spectrum, obtained by subtracting the modeled continuum from the observed spectrum, is fit by a combination of linear continuum and Gaussian components. In addition to narrow emission, the Hα+[NII] feature shows clear evidence for a broader feature. A flux ratio of 1 : 3 was assumed for the [N II] doublet, as dictated by the branching ratio [37]; the [O I] feature and the [N II] and [S II] doublets were assumed to share common velocity (red) shifts and widths. The best-fitting Gaussian parameters were derived via an interactive χ2 minimization, using SPECFIT [38]. Because of the generally low signal to noise of the 2-d spectrum we are not able to test whether our measurements of the broad Hα feature could be corrupted by a possible rotating disk as in the case of M84 (e.g., [39]); the match in the widths of the [S II] and [N II] narrow features (as long as that of the narrow Hα) argue however against such a significant effect.Figure 1
Spectral fits to theHST-STIS spectrum of NGC 4736, based on a 5-pixel extraction aperture, which corresponds to a total emission area of 0.025 arcsec2. Upper panel: The thin continuous line shows the observed data, while the green line shows the best-fit stellar population model obtained by clipping out the hashed emission-line region. Middle panel: The continuum subtracted spectrum, that is, the emission-line component, together with the corresponding spectral fit, that have been shifted by an additive constant to match the initial continuum level, for illustration purposes. The thin continuous lines show the individual Gaussian components (1 per narrow line), the thick grey continuous line represents the broad Hα feature, and the red line is the final fit to the observed spectrum. Lower panel: The residuals after the subtraction of the model fit.Measurements related to the broadHα feature from the STIS measurements are listed in Table 1. The fractional contribution of this broad component to the total flux of the  Hα+[NII]blend is 80%. Compared to other galaxies with broad-line emission in the Palomar survey [40], the width of the broad Hα is typical (1570 km s−1). However, at just 9.1×1037 ergs s−1, the Hα luminosity is lower than any known broad-line sources in the Palomar sample or the Sloan Digital Sky Survey [40, 41]. This includes the famous low luminosity Seyfert 1 galaxy, NGC 4395, which has a luminosity of 1.2×1038 ergs s−1 [42]  and is at a similar distance to NGC 4736. Note also that the broad Hα feature is redshifted by ~140±20 km s−1 relative to the narrow emission, that is locked at the systemic velocity of the object; velocity offsets between the broad Hα and narrow lines are not unusual [43].Table 1
BroadHα emission measurements.
QuantityValueNotesfblend0.80Fraction of broadHα toHα+[NII] blend.fHα0.93Fraction of broadHα to total HαFWHM(Hαbroad)1570±110; inkms-1Δv140±20; inkms-1; broad relative to narrow HαlogF(Hαbroad)−13.38Observed flux (ergs-1cm-2)logL(Hαbroad)37.96Observed luminosity(ergs-1)aaassuming distanced=4.3Mpc.Evidence for broadHα emission has also been presented in a recent PCA tomography study applied to this nucleus [44]. The observations reported in this case come from the Gemini Multi Object Spectrograph (GMOS)-IFU data cube and have been obtained 4 years after the HST-STIS spectrum had been acquired. Thus, the broad Hα line associated with the nucleus of M94 appears to be persistent for at least 4 years. They derive a broad Hα luminosity of 6×1038 ergs s−1, ≳6 times brighter in the GMOS data than the one we detect in the STIS spectrum. The difference could result from true variability in the broad-line luminosity or could be due to differences in aperture size or placement. We note that Steiner et al. [44] suggest the broad-line region is offset from the photocenter of the galaxy by 0.′′15, while our spectral 0.′′25×0.′′1 aperture coincides with the photocenter. Even with their higher Hα luminosity, this source is among the least luminous known broad-line AGN.


## 2.2. Ground-Based Optical Spectroscopy and the Spectral Classification

Ground-based optical spectra of the NGC 4736 nucleus are available from the Palomar spectroscopic survey by Ho et al. [43] and the integrated spectrophotometric survey of Moustakas and Kennicutt [45] conducted with the 2.3 m Bok telescope, as well as from new data we acquired on February 2008 at the MMT Observatory. The MMT spectrum is obtained with the Blue Chanel Spectrograph, the 500 grooves/mm grating used in first order with the 1′′ slit, and covers λλ3800–7000 with 3.6 Å resolution. None of these data show indications of broad Hα emission. This outcome is not surprising given the feature’s flux in the HST spectrum, which would be very difficult to discern in the 1′′ or 2′′ apertures employed in these ground-based observations, which are at least one order of magnitude larger than that used in the HST-STIS observations. Through a 2′′ aperture (~40 pc), NGC 4736’s emission complex Hα+[N II] is generally heavily swamped by the host stellar light.Probably as expected, measurements of NGC 4736’s narrow-line emission are consistent with this system being powered at least partially by a nonstellar source, regardless of the resolution of its observation. Figure2 shows the location of this nucleus in the 3-dimensional diagnostic diagram usually employed in classifying emission-line galaxies (e.g., [3, 4, 6]), for different sets of data mapping ~10 to ~100 pc. As with the HST-STIS observations, the ground-based measurements are performed after the host stellar contribution has been subtracted from the observed spectrum by means of absorption galaxy template fits. Our MMT spectrum offers only an upper limit for the [O III ]/Hβ while the Palomar line flux measurements are at least 50% uncertain. We list the measurements of the nebular emission of this nucleus as observed by Palomar, Bok, MMTO, and HST-STIS in Table 2. Because measurements of the [O III]/Hβ ratio are not available in the HST spectrum, we show the high-resolution measurements in Figure 2 using the [O III]/Hβ from the Palomar catalog; it is readily apparent that the high-resolution data are consistent with a LINER or a Seyfert classification for this object and reveal its type 1 (broad-lined) character.Table 2
Emission linemeasurementsa.
line nameHST-STISPalomarbBokcMMTOd[O I] λ63005.03±0.316.032.8±0.97.0±1.0Hα (narrow)3.00±0.32> 25.1213.5±1.922.0±2.4[N II] λ65835.25±0.2554.0121.5±1.539.3±2.9[S II] λ67162.48±0.1818.598.3±1.115.5±0.3[S II] λ67312.16±0.1816.337.2±1.19.8±0.3aAll fluxes are in units of  10-15ergs-1cm-2 and represent the observed values, not corrected for reddening.bnonphotometric conditions; line ratios are at least 50% uncertain.c2.′′5 slit, 2.3 m telescope; photometric conditions [45].d1′′slit; fluxes are in units of  10-15ergs-1cm-2.Figure 2
NGC 4736 within emission-line diagnostic diagrams. Filled circles reflect measurements based on the Palomar observations. Diamonds represent measurements from Moustakas and Kennicutt Jr. [45], while crosses show measurements from our MMT spectrum. The vertical bands indicate the HST-STIS measurements and the associated errors (with no correction for reddening; [O III] and Hβ are not available in the HST spectra). The solid (black) curves indicate the Kewley et al. [6] classification, while the dotted (blue) lines indicate criteria used by Ho et al. [43]. The background grey points correspond to measurements of SDSS nearby galaxies from Constantin et al. [19].


## 2.3. The Radio, UV, and X-Ray Observations and the Astrometric Coincidence with the BroadHα Emission

Across the electromagnetic spectrum, the center of this galaxy exhibits a complex morphology. There is a plethora of bright radio, UV, and X-ray sources detected in the center of NGC 4736. There is a nuclear compact (15 GHz) radio source measured by Nagar et al. [46] that appears to be associated with that of the brightest of the two close (8.49 GHz) radio sources revealed by Körding et al. [31], and with that of the brightest of the two ultraviolet point sources, which varies on a 10-year time scale [29]. Chandra observations [30] reveal numerous discrete X-ray sources in the inner galaxy; the second brightest X-ray source, X2 (LX,2–10keV=5.9×1038 erg s-1), coincides within the errors with the nucleus position. The off-nuclear radio, UV, and X-ray sources are apparently unrelated to each other.The relation to possible optical counterparts of these observations is not well constrained in the literature. We thus investigated this issue, with a particular interest in the degree to which our newly detected broadHα coincides in position with the multi-λ detections. We used for this purpose a variety of archival images from the HST: WFPC2/PC data in the F555W filter data (PID: 5741 and 10402), HRC data in the F250W and F330W filters (PID: 9454), and NICMOS/NIC3 data in the F160W filter (PID: 9360). All data were downloaded from the HST archive, and images were drizzled together when required. Absolute astrometry was performed on these data as well as the STIS observations (taken with the slit out) using a ground-based V-band image obtained from the Spitzer Infrared Nearby Galaxies Survey (SINGS) ancillary data by Kennicutt et al. [48]. These SINGS observations were aligned to the USNO-B system using ~80 stars, and then astrometry of the HST data was obtained by degrading the resolution of the F555W image to match the SINGS image. All other HST images were then matched with the astrometry-corrected F555W frame to an accuracy of <0.′′05. The absolute error on this astrometry is about 0.′′2 and is dominated by scatter of stars in the SINGS image relative to the USNO-B positions.Figure3 illustrates the result of this data compilation and the corresponding radio, UV, and X-ray source matches. It is clear that, within<0.′′2(4 pc), there is an obvious astrometric match in the nuclear X-ray, UV, optical, radio compact sources, and the newly detected broadHαemission line. The nucleus position as observed by STIS has RA: 12 h 50 m 53 s.20, DEC: 41°07′13.′′40. The off-nuclear X1, that is the brightest compact X-ray source detected in this galaxy nucleus, has no counterpart at other wavelengths. Same is true for the off-nuclear 8.49 GHz detection, that is only 1′′ (~20 pc) away from the nucleus. The off-nuclear UV source is the only one detected in optical light.Figure 3
The mosaic-ed F555W (V-band; combination of 5 WFPC2 images taken at two different epochs) image of the NGC4736 nucleus, with intensity contours overlaid (green). The position of the STIS-nucleus is indicated (by the yellow slit) along with the locations of the 2 most luminous and closest to the nucleus hard X-ray sources, X1 and X2, as detected with Chandra (red squares), along with the position of the 15 GHz VLA radio core ([46], orange circle), and those of the two radio compact sources detected by Körding et al. [31] (white circles). The blue and cyan contours correspond to the F330W and F250W observations from Maoz et al. [29]; note the presence of NGC 4736b, ≈2.′′5 to the north of the nucleus. The absolute astrometry is good to within 0.′′2 (pixel size is 0.′′05 for F555W, and 0.′′025 for F250W and F330W).


## 3. NGC 4736’s Nuclear Emission across the Electromagnetic Spectrum andLBol

With the wealth of data available for this galaxy nucleus, we are able to build its least contaminated X-ray to radio nuclear spectral energy distribution (SED). The multiwavelength observations of the sources detected in the very central regions of NGC 4736, plotted inνLν units, are displayed in Figure 4. The X-ray detection X2 is represented here by a power law, estimated based on its photon indices Γ=1.6, where the absorbing column density is fixed at the Galactic value. Both F250W and F330W UV observations are plotted, only for the nuclear detections. The optical data are represented by the HST-STIS spectrum, featuring the strongHα emission line. Measurements corresponding to all the radio observations discussed above are indicated. For the sake of completeness, we also include in this plot lower spatial-resolution observations from SINGS IRS [48] as well. We show all of these measurements superposed on the SEDs of LLAGN from Ho [47]. No artificial normalization has been performed, and there is no correction for absorption in either the Galactic extinction (except for the X-ray data) or intrinsic to NGC 4736.Figure 4
Radio to X-ray nuclear SED of NGC 4736 superposed on SEDs of LLAGN from Ho [47]. The optical STIS spectrum stands out as an amorphous blob featuring the strong Hα emission feature. The 2–10 keV Chandra X-ray detection X2 is depicted as a power-law corresponding to Γ=1.6. The low spatial resolution Spitzer IRS observations are shown for comparison.With X2 and n4736-b as the X-ray and the 8.5 GHz counterparts of the broadHα detection, respectively and assuming that the continuum could be described as simple power laws between the points we present data for, this object’s nuclear SED corresponds to Lbol≈2.5×1040 erg s-1. Spitzer data is not included in estimating Lbol because they do not reflect the IR emission of the uncontaminated nucleus; the aperture used in these observations is very large (see Table 3), including the entire field of view shown in Figure 3. This luminosity makes NGC 4736 one of the least luminous LINERs with strong evidence of BH accretion. This source is thus a critical signpost of BH accretion at extremely low levels. Note that this object’s weak emission is most likely not caused by obscuration; Ho et al. [43] list a Balmer decrement of Hα/Hβ=3.1 (albeit highly uncertain, with a probable error of±100%), and Eracleous et al. [30] provide an upper limit for the neutral hydrogen column density of NH<3.3×1020 cm-2.Table 3
Nuclear SED data.
log(ν/Hz)log(νLν/erg/s)ResolutionObs. dateInstr.18.3838.690.′′155/2000Chandra-ACIS17.6838.420.′′155/2000Chandra-ACIS15.0839.410.′′56/2003ACS F250W14.9639.720.′′56/2003ACS F330W14.6839.910.′′16/2002HST-STIS14.6439.930.′′16/2002HST-STIS13.4840.584.′′75/2004IRS-SH13.1940.3311.′′15/2004IRS-SH/LH12.9140.8711.′′15/2004IRS-LH10.1835.790.′′151/2001VLA, A config.9.9335.470.′′246–10/2003VLA, A config.


### 3.1. Comparison with Other Observed LLAGN

Figure4 shows that NGC 4736’s SED is very similar to those of other LLAGN [47]. This suggests that, in spite of the significant difference in Lbol, there is no fundamental transition in the accretion mode in this source compared to other LLAGN.The similarity of the SED means that previously proposed bolometric corrections for LLAGN appear to work quite well for NGC 4736. The value estimated based on the correction to the [O III] line luminosity that Heckman et al. [49] proposed to work well for Seyferts is Lbol≈3500×L[O iii]=4.5×1040 erg s-1 (here L[O III] is measured from the ground within the Palomar survey and is not corrected for intrinsic reddening). If the more recent assessment of the Lbol=600×L[O III] bolometric correction of kauffmann and Heckman [50] is used, with L[O III] corrected for reddening based on the Balmer decrement listed above and a τ∝λ-0.7 attenuation law [51], then Lbol=0.93×1040 erg s-1, which compares satisfactorily with the measured Lbol. The average bolometric correction to the observed 2–10 keV X-ray luminosity proposed by Ho [2], with L/LX=16, results in a somewhat lower value, 0.94×1040 erg s-1; with the more recent L/LX=50 average bolometric correction to the 2–10 keV luminosity of Eracleous et al. [52], Lbol=1.8×1040 erg s-1 and, thus, very close to our integrated value. Within a typical uncertainty of ~50%, all of these estimates are consistent with the measured Lbol value. Other, more uncertain bolometric indicators are also consistent with our integration: the correlation between the mid-IR [Ne V] emission-line luminosity and bolometric luminosity derived by Satyapal et al. [53] for a small sample of much brighter nearby AGN, that has a large scatter (~1order of magnitude), gives a Lbol=1.6×1040 erg s-1, which is very close to our measured value.This good match among quite a variety of bolometric estimators derived independently from the multiwavelength properties of AGN supports the AGN interpretation for this system. However, there are two details of its emission spectrum which are unlike typical broad-lined LLAGN.First, despite the fact that the broadHα emission of NGC 4736 is one of the weakest measured among type 1 AGN, the total Hα emission is peculiarly strong relative to the X-ray counterpart. The ratio LX(2–10 keV)/LHα is ~6, and thus lower than the median of ~15 exhibited by the type 1 AGN (and low-z quasars) included in, for example, Ho et al. [54] study that revealed a relatively tight correlation between the two types of emission. With the higher broad Hα flux from Steiner et al. [44], the ratio is even lower, LX(2–10 keV)/LHα≈1. This finding is surprising given that the low-luminosity nearby galaxies that deviate from the LX-LHα linear scaling toward lower values of the LX/LHα ratios are the type 2 sources, mostly transition objects, where the ionization mechanism is not necessarily dominated by an AGN type of source. Following, for example, Ho et al. [54] arguments, the unusually low LX/LHα ratio measured in the nucleus of NGC 4736 suggests that either: (1) the optical line emission is not powered exclusively by a central AGN, or (2) the X-ray emission in this system arises, at least partially, from a non-AGN source, for example, an X-ray binary. A final possibility is that the source is highly variable (as suggested by the UV, X-ray, and Hα observations), and thus, the unusual LX/LHα is the result of measurements made at different times.Second, NGC 4736’s nebular emission shows a number of peculiarities. To start with, the electron density in NGC 4736’s nebular emission is very low, and it does not show the typical gradient exhibited by AGN, that is, increase toward the more nuclear regions (Constantin et al. 2012, in prep.). In both the Palomar andHST-STIS observations, the ratio [S II] λλ 6716/6731 is ~1.26 (Table 2), implying an electron density of ne~102 cm-3, which lies at the lowest end in the distribution of particle densities measured in the Palomar objects [55].Given the lack of a density gradient in the line-emitting region it is then rather peculiar to observe a significant increase (by a factor of~6) in the [O I]/Hα line flux ratio in the HST spectrum relative to the large aperture observations, which would usually be interpreted as an indication for a more pronounced AGN-like ionization in the more central regions. There are two reasons for this: (1) [O I] requires a significantly hard radiation field, that is, that of an AGN, to sustain a sufficiently extensive partially ionized zone in clouds optically thick to Lyman continuum and thus to produce such a strong feature. Since the ionization potential of [O I] matches that of H very well, large differences in the [O I]/Hα ratios are expected between accretion and nonaccretion sources. (2) given the lack of a density stratification in this nucleus, the degree of ionization of the emitting gas is expected to diminish with radius as the ionizing radiation emerging from a nuclear source falls off in density as r-2 and thus produce strong gradients in the [O I]/Hαratio.In summary, NGC 4736 may be a broad-lined AGN, with an SED similar to other LLAGN, but it is atypical in several aspects whose physical origin remains unclear. Given the above-listed possible explanations for the peculiarLX/LHα ratio, we evaluate in the next subsection whether photoionization by the weak AGN in this system is sufficiently powerful to balance the emission cooling in this system.


### 3.2. Comparison of Ionizing and Emission-Line Power

The multiwavelength observations of this system allow for a relatively rigorous assessment of whether the photoionization by this system’s AGN can power the measured emission-line luminosities and in particular that of the broad Hα component. Following Eracleous et al. [18], we can run an energy budget test via a direct comparison of the Hα luminosity and count rate with the ionizing luminosity Li=L1Ry-100keV and the ionizing photon rate Qi=Q1Ry-100keV.It is important to treat separately the broad and narrow emission line features as there is strong evidence that they originate from regions of significantly different physical conditions. The broad Hα comes from a much more compact and much denser emitting gas than the narrow Balmer and forbidden lines. The difference in density is at least 3 orders of magnitude; the critical density for collisional excitation of the [O I] (~2×106 cm-3) can be used to estimate the gas density of the broad-line region, as [O I] does not exhibit a broad component in this nucleus’ emission. The difference in size is expected to be ~2 orders of magnitude (e.g., reverberation mapping of nearby Seyfert galaxies measured to be <1 week, [56]). Thus, the mechanisms that can operate in these two regions are expected to be qualitatively different. The photon and energy balance conditions should reflect these differences and thus should differ as well; in particular, potential contributions to the Hαemission via collisional excitation should be minimal for the narrow-line-emitting region, but important for the broad component.Based on photoionization models (i.e., Cloudy, v94.0; see [57]) computed by Lewis et al. [58] for a wide range of ionization parameters, densities, and metallicities, Eracleous et al. [18] find that energy balance in a line-emitting nebula requires that Li>18(±2)LHα/fc, where fc is the covering factor or the fraction of the ionizing luminosity of the AGN that is absorbed by the line-emitting gas. These models are covering electron densities that appear to encompass both the broad- and the narrow-line-emitting regions in NGC 4736; thus, with only a fraction of 10% of ionizing photons being absorbed by the line-emitting gas, the minimum energy balance condition for AGN ionization is given by Li/LHα>180, for both narrow and broad emission features. In the same time, however, a minimum requirement for photon balance can be quite different for the two emission regions: for the narrow-line component Qi>2.2QHα, corresponding to the case B recombination (i.e., one Hα photon is emitted for every 2.2 recombinations); for the denser broad-line-emitting region, the number of H-alpha photons that can be produced for each ionization can be at least 7-8 times higher than the standard case B estimate or Qi≳0.25QHα (e.g., [37]).With theLi and Qi values already calculated by Eracleous et al. [18] by integrating M94’s nuclear SED assuming that pairs of points could be connected by a power law, and with the Hα measurements from all HST and ground-based observations presented above, we can proceed with the comparison. Table 4 lists the LHα, QHα, and associated ratios Li/LHα, Qi/QHα for all of these optical spectroscopic observations, where the narrow and broad Hα measurements are shown separately.Table 4
Ionizing and emission-line power.
logLHαalogQHαaLi=LHαbQi=QHαcPalomar>37.7>49.3<93<2.75Bok37.549.0165.94.9MMTO37.549.2104.73.1HST-STIS, narrow36.848.475922.4HST-STIS, broad37.949.557.51.6PCA tomography38.850.38.70.25athe luminosities are measured in ergs-1 and the photon rates in s-1.bthe minimum energy balance condition for AGN ionization is given byLi/LHα≳180 (when a fraction of 10% of ionizing photons photons is absorbed by the line-emitting gas).cthe minimum photon balance condition for AGN ionization is given by Qi/QHα>2.2 for the narrow line regions, and Qi/QHα≳0.25 for the broad component.It is readily apparent that whileLi/LHα is in general >18 (fc=1), it is almost never >180 (fc=0.1); the only exception is the narrow Hα emission measured in the HST aperture. Thus, a dominant AGN ionization of the narrow-line region in this nucleus is definitely possible in the HST aperture, but only for fc≳18% at larger radial distances. For the HST broad Hα feature, an AGN ionization is possible only if fc≳30%. The PCA tomography measurement of the broad Hα argues, however, against a balanced energy budget originating entirely in an AGN-like power mechanism, even when a maximum covering factor is considered; the AGN-produced energy falls short of the required amount by at least 50%. Interestingly, the measured Qi/QHα ratio is well within the required photon balance corresponding to an AGN excitation for both the narrow and the broad emission features. Thus, with one clear exception, the AGN in NGC 4736 appears to be capable of providing enough photons to explain the observed Hα luminosity but only for relatively high fc values.Simply because the AGN-like SED of NGC 4736 can explain the majority of the ionization energy and photon rate does not imply that the actual mechanism is an AGN, alternative excitation mechanisms must be explored. Possible options are (i) we are missing ionizing photons from accretion onto the central BH, or (ii) there are other power sources that could make up the power deficit in this system, particularly for producing the broadHα feature.The first alternative could be possible if we were observing an “echo” of a previous epoch of more violent accretion, a few hundred years ago. This idea has been explored by Eracleous et al. [59] who showed that the reverberation of an ionizing flare in the nebula can produce LINER-like emission-line ratios. In this scenario, it is expected that the central UV source and the [O III] line would also follow the decay of the ionizing continuum; while the UV observations are not providing clear evidence for such a decay over the course of one decade, the multiple optical spectroscopic observations that we present in Figure 2 are consistent with a possible decrease in time (a few years) in the [O III] flux (in this model, Hβ, [S II], [O I] are expected to decay very slowly, in 60–250 years). Nevertheless, this duty-cycle hypothesis also requires that the broad-line emission fades immediately if the ionizing continuum declines, and the PCA tomography observations show that this is not happening in this source. The echo of such an ionizing continuum flare should also be detectable in a narrow-band [O III ]λ5007 image in the form of a ring, which is not readily observed in M94.The second alternative of power sources other than accretion has been often proposed in explaining the emission-line spectra of LINERs, with mixed success (see Section1). The most probable alternative sources appear to be the mechanical power delivered by compact radio jets, along with photoionization by young or post-AGB stars from old or intermediate-age stellar populations. Shock models are highly unlikely to produce broad emission features and, thus, are not favored in this case. Recent star-formation activity remains, however, a viable option, particularly in light of relatively new discoveries of peculiar supernovae with broad Hα features that do not appear to fade in time. This idea is discussed in more detail in Section 5.


## 3.1. Comparison with Other Observed LLAGN

Figure4 shows that NGC 4736’s SED is very similar to those of other LLAGN [47]. This suggests that, in spite of the significant difference in Lbol, there is no fundamental transition in the accretion mode in this source compared to other LLAGN.The similarity of the SED means that previously proposed bolometric corrections for LLAGN appear to work quite well for NGC 4736. The value estimated based on the correction to the [O III] line luminosity that Heckman et al. [49] proposed to work well for Seyferts is Lbol≈3500×L[O iii]=4.5×1040 erg s-1 (here L[O III] is measured from the ground within the Palomar survey and is not corrected for intrinsic reddening). If the more recent assessment of the Lbol=600×L[O III] bolometric correction of kauffmann and Heckman [50] is used, with L[O III] corrected for reddening based on the Balmer decrement listed above and a τ∝λ-0.7 attenuation law [51], then Lbol=0.93×1040 erg s-1, which compares satisfactorily with the measured Lbol. The average bolometric correction to the observed 2–10 keV X-ray luminosity proposed by Ho [2], with L/LX=16, results in a somewhat lower value, 0.94×1040 erg s-1; with the more recent L/LX=50 average bolometric correction to the 2–10 keV luminosity of Eracleous et al. [52], Lbol=1.8×1040 erg s-1 and, thus, very close to our integrated value. Within a typical uncertainty of ~50%, all of these estimates are consistent with the measured Lbol value. Other, more uncertain bolometric indicators are also consistent with our integration: the correlation between the mid-IR [Ne V] emission-line luminosity and bolometric luminosity derived by Satyapal et al. [53] for a small sample of much brighter nearby AGN, that has a large scatter (~1order of magnitude), gives a Lbol=1.6×1040 erg s-1, which is very close to our measured value.This good match among quite a variety of bolometric estimators derived independently from the multiwavelength properties of AGN supports the AGN interpretation for this system. However, there are two details of its emission spectrum which are unlike typical broad-lined LLAGN.First, despite the fact that the broadHα emission of NGC 4736 is one of the weakest measured among type 1 AGN, the total Hα emission is peculiarly strong relative to the X-ray counterpart. The ratio LX(2–10 keV)/LHα is ~6, and thus lower than the median of ~15 exhibited by the type 1 AGN (and low-z quasars) included in, for example, Ho et al. [54] study that revealed a relatively tight correlation between the two types of emission. With the higher broad Hα flux from Steiner et al. [44], the ratio is even lower, LX(2–10 keV)/LHα≈1. This finding is surprising given that the low-luminosity nearby galaxies that deviate from the LX-LHα linear scaling toward lower values of the LX/LHα ratios are the type 2 sources, mostly transition objects, where the ionization mechanism is not necessarily dominated by an AGN type of source. Following, for example, Ho et al. [54] arguments, the unusually low LX/LHα ratio measured in the nucleus of NGC 4736 suggests that either: (1) the optical line emission is not powered exclusively by a central AGN, or (2) the X-ray emission in this system arises, at least partially, from a non-AGN source, for example, an X-ray binary. A final possibility is that the source is highly variable (as suggested by the UV, X-ray, and Hα observations), and thus, the unusual LX/LHα is the result of measurements made at different times.Second, NGC 4736’s nebular emission shows a number of peculiarities. To start with, the electron density in NGC 4736’s nebular emission is very low, and it does not show the typical gradient exhibited by AGN, that is, increase toward the more nuclear regions (Constantin et al. 2012, in prep.). In both the Palomar andHST-STIS observations, the ratio [S II] λλ 6716/6731 is ~1.26 (Table 2), implying an electron density of ne~102 cm-3, which lies at the lowest end in the distribution of particle densities measured in the Palomar objects [55].Given the lack of a density gradient in the line-emitting region it is then rather peculiar to observe a significant increase (by a factor of~6) in the [O I]/Hα line flux ratio in the HST spectrum relative to the large aperture observations, which would usually be interpreted as an indication for a more pronounced AGN-like ionization in the more central regions. There are two reasons for this: (1) [O I] requires a significantly hard radiation field, that is, that of an AGN, to sustain a sufficiently extensive partially ionized zone in clouds optically thick to Lyman continuum and thus to produce such a strong feature. Since the ionization potential of [O I] matches that of H very well, large differences in the [O I]/Hα ratios are expected between accretion and nonaccretion sources. (2) given the lack of a density stratification in this nucleus, the degree of ionization of the emitting gas is expected to diminish with radius as the ionizing radiation emerging from a nuclear source falls off in density as r-2 and thus produce strong gradients in the [O I]/Hαratio.In summary, NGC 4736 may be a broad-lined AGN, with an SED similar to other LLAGN, but it is atypical in several aspects whose physical origin remains unclear. Given the above-listed possible explanations for the peculiarLX/LHα ratio, we evaluate in the next subsection whether photoionization by the weak AGN in this system is sufficiently powerful to balance the emission cooling in this system.


## 3.2. Comparison of Ionizing and Emission-Line Power

The multiwavelength observations of this system allow for a relatively rigorous assessment of whether the photoionization by this system’s AGN can power the measured emission-line luminosities and in particular that of the broad Hα component. Following Eracleous et al. [18], we can run an energy budget test via a direct comparison of the Hα luminosity and count rate with the ionizing luminosity Li=L1Ry-100keV and the ionizing photon rate Qi=Q1Ry-100keV.It is important to treat separately the broad and narrow emission line features as there is strong evidence that they originate from regions of significantly different physical conditions. The broad Hα comes from a much more compact and much denser emitting gas than the narrow Balmer and forbidden lines. The difference in density is at least 3 orders of magnitude; the critical density for collisional excitation of the [O I] (~2×106 cm-3) can be used to estimate the gas density of the broad-line region, as [O I] does not exhibit a broad component in this nucleus’ emission. The difference in size is expected to be ~2 orders of magnitude (e.g., reverberation mapping of nearby Seyfert galaxies measured to be <1 week, [56]). Thus, the mechanisms that can operate in these two regions are expected to be qualitatively different. The photon and energy balance conditions should reflect these differences and thus should differ as well; in particular, potential contributions to the Hαemission via collisional excitation should be minimal for the narrow-line-emitting region, but important for the broad component.Based on photoionization models (i.e., Cloudy, v94.0; see [57]) computed by Lewis et al. [58] for a wide range of ionization parameters, densities, and metallicities, Eracleous et al. [18] find that energy balance in a line-emitting nebula requires that Li>18(±2)LHα/fc, where fc is the covering factor or the fraction of the ionizing luminosity of the AGN that is absorbed by the line-emitting gas. These models are covering electron densities that appear to encompass both the broad- and the narrow-line-emitting regions in NGC 4736; thus, with only a fraction of 10% of ionizing photons being absorbed by the line-emitting gas, the minimum energy balance condition for AGN ionization is given by Li/LHα>180, for both narrow and broad emission features. In the same time, however, a minimum requirement for photon balance can be quite different for the two emission regions: for the narrow-line component Qi>2.2QHα, corresponding to the case B recombination (i.e., one Hα photon is emitted for every 2.2 recombinations); for the denser broad-line-emitting region, the number of H-alpha photons that can be produced for each ionization can be at least 7-8 times higher than the standard case B estimate or Qi≳0.25QHα (e.g., [37]).With theLi and Qi values already calculated by Eracleous et al. [18] by integrating M94’s nuclear SED assuming that pairs of points could be connected by a power law, and with the Hα measurements from all HST and ground-based observations presented above, we can proceed with the comparison. Table 4 lists the LHα, QHα, and associated ratios Li/LHα, Qi/QHα for all of these optical spectroscopic observations, where the narrow and broad Hα measurements are shown separately.Table 4
Ionizing and emission-line power.
logLHαalogQHαaLi=LHαbQi=QHαcPalomar>37.7>49.3<93<2.75Bok37.549.0165.94.9MMTO37.549.2104.73.1HST-STIS, narrow36.848.475922.4HST-STIS, broad37.949.557.51.6PCA tomography38.850.38.70.25athe luminosities are measured in ergs-1 and the photon rates in s-1.bthe minimum energy balance condition for AGN ionization is given byLi/LHα≳180 (when a fraction of 10% of ionizing photons photons is absorbed by the line-emitting gas).cthe minimum photon balance condition for AGN ionization is given by Qi/QHα>2.2 for the narrow line regions, and Qi/QHα≳0.25 for the broad component.It is readily apparent that whileLi/LHα is in general >18 (fc=1), it is almost never >180 (fc=0.1); the only exception is the narrow Hα emission measured in the HST aperture. Thus, a dominant AGN ionization of the narrow-line region in this nucleus is definitely possible in the HST aperture, but only for fc≳18% at larger radial distances. For the HST broad Hα feature, an AGN ionization is possible only if fc≳30%. The PCA tomography measurement of the broad Hα argues, however, against a balanced energy budget originating entirely in an AGN-like power mechanism, even when a maximum covering factor is considered; the AGN-produced energy falls short of the required amount by at least 50%. Interestingly, the measured Qi/QHα ratio is well within the required photon balance corresponding to an AGN excitation for both the narrow and the broad emission features. Thus, with one clear exception, the AGN in NGC 4736 appears to be capable of providing enough photons to explain the observed Hα luminosity but only for relatively high fc values.Simply because the AGN-like SED of NGC 4736 can explain the majority of the ionization energy and photon rate does not imply that the actual mechanism is an AGN, alternative excitation mechanisms must be explored. Possible options are (i) we are missing ionizing photons from accretion onto the central BH, or (ii) there are other power sources that could make up the power deficit in this system, particularly for producing the broadHα feature.The first alternative could be possible if we were observing an “echo” of a previous epoch of more violent accretion, a few hundred years ago. This idea has been explored by Eracleous et al. [59] who showed that the reverberation of an ionizing flare in the nebula can produce LINER-like emission-line ratios. In this scenario, it is expected that the central UV source and the [O III] line would also follow the decay of the ionizing continuum; while the UV observations are not providing clear evidence for such a decay over the course of one decade, the multiple optical spectroscopic observations that we present in Figure 2 are consistent with a possible decrease in time (a few years) in the [O III] flux (in this model, Hβ, [S II], [O I] are expected to decay very slowly, in 60–250 years). Nevertheless, this duty-cycle hypothesis also requires that the broad-line emission fades immediately if the ionizing continuum declines, and the PCA tomography observations show that this is not happening in this source. The echo of such an ionizing continuum flare should also be detectable in a narrow-band [O III ]λ5007 image in the form of a ring, which is not readily observed in M94.The second alternative of power sources other than accretion has been often proposed in explaining the emission-line spectra of LINERs, with mixed success (see Section1). The most probable alternative sources appear to be the mechanical power delivered by compact radio jets, along with photoionization by young or post-AGB stars from old or intermediate-age stellar populations. Shock models are highly unlikely to produce broad emission features and, thus, are not favored in this case. Recent star-formation activity remains, however, a viable option, particularly in light of relatively new discoveries of peculiar supernovae with broad Hα features that do not appear to fade in time. This idea is discussed in more detail in Section 5.


## 4. The Black Hole Mass

Given the unusual energetics of M94’s nucleus, it is important to investigate whether the BH mass estimators derived for rapidly accreting Seyferts, which appear to be widely used for AGN, also work in the low-luminosity regime flagged by this particular system.There are a variety of indirect methods that can be used to estimate the mass of BHs in galaxy centers. The available multiwavelength measurements of the nuclear emission for NGC 4736 allow the calculation ofMBH based on four different techniques, as well as with the M-σ* relation. In this section, we explore and compare the results of these five methods, along with their consequences for this object’s energetics, and a comparison with a recent dynamical measurement of 6.68(5.14–8.22)×106M⨀ for this object, which is listed in Kormendy et al. [60] as obtained from (Gebhardt et al. 2012, in prep.). We first present the methods and associated MBH calculations and then discuss the shortcomings of each measurement.(1) Using theM-σ* relation, established for quiescent nearby galaxies, including both ellipticals and spirals with classical bulges, as quantified by Gültekin et al. [1], for σ*=110±5km s-1 [61], MBH=1.05±0.64×107M⨀; within errors, this value agrees well with Gebhardt et al. dynamical measurement. With this value, Lbol/LEdd≈2×10-5, which is consistent with the range of values within LINERs are expected to (e.g., [2, 62]).(2) Using theHST measurements of the FWHM and the luminosity of the broad Hα component (see Table 1) within the scaling relation based solely on observations of this broad emission feature, as derived by Greene and Ho [63], we obtain MBH≈3×104M⨀; the fractional uncertainty associated with this measurement is ~30%, and includes both the scatter in the scaling relation and the errors in the line measurements. With the broad Hα line luminosity measured in the PCA tomography study [44], the BH mass would increase (by a factor of 6.5) to MBH=1.9×105M⨀. The difference in these two values could be considered the most conservative uncertainty associated with this BH mass estimate. The corresponding Eddington ratios for these BH masses are Lbol/LEdd≈1–7×10-3.(3) Mid IR detection and measurements of the[Ne V] (14.32 μm) and [O IV] (25.89 μm) emission lines [64] give, via the empirical correlations between the MIR line luminosities and reverberation mapping-based MBH values presented by Dasyra et al. [65], a black hole mass of 2.3±0.4×105M⨀ and 1.7±0.5×105M⨀, respectively. The corresponding Eddington ratio is in this case Lbol/LEdd≈1.0×10-3. Note that the scatter adopted for these relations is only a lower limit of the real value; thus, the uncertainty may be larger.(4) The normalized X-ray excess variance method, as described in Papadakis et al. [66], applied to Chandra observations [30] gives MBH=2.5±1.7×105M⊙. To be specific, we used for this calculation the excess variance σ=0.06±0.04, Lbol=2.5×1040 erg s-1 as derived in Section 3, and νlf=1/T=1/14h-1 (where T is the length of the light curve), to estimate the break frequency νbf and then MBH via equations 4 and 6, respectively, of Papadakis et al. [66].(5) TheMBH of this system can also be obtained via the “fundamental plane of black hole activity” that relates black hole mass to the emitted compact radio LR=νLν(5 GHz) and hard X-ray luminosities LX(2–10 keV) and spans nine orders of magnitude in black hole mass. We calculate LR=1.7×1035 erg s-1 using the two 8.5 GHz and 15 GHz radio measurements presented in this paper (Figure 3, Table 3), with a flux modeled by a powerlaw Sν∝ν-α. Using the empirical fits of Merloni et al. [67] for the fundamental plane relation, log(MBH/M⊙)=5.9±1.1 for the black hole in this system; note that this measurement and its associated errors embrace all of the estimates presented above and, thus, do not provide any additional constraint to the MBH. Interestingly, the latest derivation of the fundamental plane relation [68], applied to very low nuclear galactic luminosities, provides a log(MBH/M⊙)=7.2±0.4 for this system, which remains consistent only with the M-σ* value, as it departs considerably from those provided by the relations employing AGN emission.There is a significant inconsistency between the value given by theM-σ* relation, which is supported by Gebhardt’s dynamical measurement, and those based on the AGN emission. The M-σ* relation suggests a BH mass in NGC 4736 of 1×107M⊙ and a correspondingly low Lbol/LEdd of ~2×10-5. Three of the estimates based on the AGN emission converge to black hole masses of ~105M⊙ (and Lbol/LEdd~10-3), showing a surprisingly consistent departure of two orders of magnitude from the M-σ* estimate. The fundamental plane relation provides a MBH value right in between these two different situations, however, with no real additional constraint, due to its associated large uncertainty. These differences are somewhat puzzling given that the BH mass estimates based on nuclear emission properties are all calibrated to follow the M-σ* for high-mass black holes (MBH≳106M⊙). It is however true that the calibrators are biased toward nearby Seyfert galaxies with much higher Eddington ratios than that of NGC 4736. We briefly discuss in the following subsections more specific weaknesses of each measurement.


### 4.1. Caveats of theM-σ* Relation

TheM-σ* relation is expected to provide a reliable estimate of the BH mass as it is based on the strong correlation between dynamical mass measurements of supermassive BHs and their host properties [1, 69–71]. This relationship is derived primarily from ellipticals and spirals with classical bulges (formed during major mergers). Recent observations suggest it may not be valid for samples of later-type spirals which more commonly host pseudobulges (formed via secular disk processes) (e.g., [72]). Because the distinction between classical and pseudobulges is based on formation, it does not simply correlate with observable properties [73, 74]. In general pseudobulges are less luminous, have lower bulge-to-total ratios, and have on-going star formation and lower Sérsic indices than classical bulges. The possibility that the NGC 4736 is in fact a pseudobulge provides a solution to the apparent conflict in BH mass estimates.Nevertheless, pseudobulges are difficult to identify, and there is not yet a consensus as to what defines them. To complicate things further, classical and pseudobulges can exist within the same galaxy (e.g., NGC 2787; [75]). The classification of the NGC 4736 bulge is ambiguous: Fisher and Drory [74] classifies it as a pseudobulge based on its nuclear spiral and bar and low Sérsic index (n=1.3). However, they also find it has a low star formation rate, more typical of classical bulges, and thus classify it as an “Inactive Pseudobulge;’’ if NGC 4736 hosts a pseudobulge, it is an atypical one. We have created a surface brightness profile from NICMOS and 2MASS H-band data from the Large Galaxy Atlas [76] and found results that conflict with the fits of Fisher and Drory [74]. Specifically, we find Sérsic indices for the bulge of n=2.3–3.0 depending on the radial range and type of fit (single versus double Sérsic), which are consistent with a classical one (as shown in [74]).Our surface brightness profile fits also reveal the presence of a nuclear star cluster within the central~0.6′′ (~12 pc), with an H-band magnitude of ~12.5. Such nuclear star clusters are common in early-type spiral galaxies [77]. The luminosity and mass of nuclear star clusters are known to scale with bulge luminosity and mass [78–80], and the NGC 4736 nuclear star cluster has a luminosity that is 0.1% of its bulge, typical for nuclear star clusters [81]. Nuclear star clusters commonly coexist with black holes, but there are a very limited number of cases where masses for both can be estimated [82]. In these cases, including the Milky Way, the BH mass is similar to the mass of the nuclear star cluster within an order of magnitude. There is also some evidence that the ratio of BH mass to nuclear star cluster mass increases with spheroid mass [83]. For NGC 4736, assuming an old population with an H-band M/L~0.7, the nuclear cluster would have a mass of ~2×107M⊙. This mass is quite similar to the M-σ* BH mass estimate, suggesting thus the presence of a similar-sized BH.


### 4.2. Caveats of the BLR Scaling Relation

The BLR scaling relations were derived using high-luminosity systems (i.e., Seyferts, not LINERs) and have been scaled to match the dynamical black hole detections of BHs with masses >106M⊙. Thus, their applicability to NGC 4736, where they yield an estimate of ~105M⊙, is a (perhaps unwise) extrapolation. These relations assume that the BLR is virialized due to proximity to the BH [84] and thus the 5100 Å continuum luminosity correlates with the emissivity-weighted radius of the BLR and thus with the BH mass [85]. The overall errors associated with BH mass measurements based on these relations do not exceed 0.5 dex [86–88], and there appears to be good consistency with BH masses obtained via the M-σ* relation for relatively bright AGN (L>1042 erg s-1), with MBH≳106M⊙ [89]. More recent calibrations of the radius-luminosity relationship on which these techniques are based infer that BH masses have been overestimated, however, only by up to a factor of ~3 [90]. This latter study also indicates a trend toward larger uncertainties and larger amount of overestimation in the BH with decreasing luminosity; however, no conclusive results are yet available for this regime.


### 4.3. Caveats of the MIR Line Correlation

The Dasyra et al. [65] empirical relation between the MIR line emission properties and the BH mass is derived using reverberation mapping BH masses and, thus, as with the previous method, may not apply to systems with MBH<106M⊙. The relation holds for systems with Lbol/LEdd>0.003, but not necessarily beyond this range. The MBH≈105M⊙ given by this method places this object at the low end of this Lbol/LEdd range; on the other hand, a more massive and thus quiescent BH would correspond to a [Ne V] luminosity of a few orders of magnitude lower than that measured. The scatter associated with this relation remains in average 0.5 dex and thus cannot account for the ~2-3 orders of magnitude difference between this value and that obtained via M-σ*.


### 4.4. Issues with Estimating the BH Mass from X-Rays

The BH mass derivation based on X-ray variability, or more precisely on the relation between the excess variance andMBH, relies on the hypothesis of a universal power spectral density function (PSD) shape and amplitude in AGN, which is based on the idea that the X-ray variability mechanism and the accretion efficiency are the same for all AGN, at all redshifts. These assumptions appear to hold for the objects involved in deriving and testing this relation [66]; that sample is, however, small and rather biased toward luminous (X-ray) sources, with log LX/(ergs-1)>41.5. None of the 2–10 keV sources detected in the nucleus of NGC 4736 are in this luminosity range. Nevertheless, Galactic BHs (Mbh<102M⊙) in their hard states show variability properties that match well those of AGN, both of these types of sources falling on the same projection of the TB-Mbh-Lbol/LEdd plane, where TB=1/νbf, suggesting that an extrapolation to the intermediate mass BH (or lower L) level is practicable.


### 4.5. Problems with the Fundamental Plane

Finally, it appears that the main problem with the fundamental plane is that the relation is not sharpened enough to provide strong constraints on BH masses. Because the relation spans nine orders of magnitude in BH mass, it is expected to equally apply to any value ofMbh in the range we are interested in. It is also the case that a wide variety of BH accretion models (e.g., with efficient and inefficient flows for the X-ray emission or associating X-ray flux with synchrotron emission near the base of a jet) are consistent with this relation (e.g., [91, 92]), suggesting that a large diversity of accretion modes or rates are accommodated.Nevertheless, the fundamental plane remains best constrained only for systems withMBH≳106M⊙ and for very low nuclear luminosities, that is, with negligible or zero AGN contribution that allow a dynamical measurement of their BH mass [68]. In this regime, the latest derivation of the fundamental plane relation provides for M94 a MBH value consistent with the M-σ* estimate and the dynamical measurement, while it would not reliably constrain a MBH≲106M⊙ value. The scatter in the fundamental plane relation increases for lower BH masses, with higher Eddington ratios [68, 93].


## 4.1. Caveats of theM-σ* Relation

TheM-σ* relation is expected to provide a reliable estimate of the BH mass as it is based on the strong correlation between dynamical mass measurements of supermassive BHs and their host properties [1, 69–71]. This relationship is derived primarily from ellipticals and spirals with classical bulges (formed during major mergers). Recent observations suggest it may not be valid for samples of later-type spirals which more commonly host pseudobulges (formed via secular disk processes) (e.g., [72]). Because the distinction between classical and pseudobulges is based on formation, it does not simply correlate with observable properties [73, 74]. In general pseudobulges are less luminous, have lower bulge-to-total ratios, and have on-going star formation and lower Sérsic indices than classical bulges. The possibility that the NGC 4736 is in fact a pseudobulge provides a solution to the apparent conflict in BH mass estimates.Nevertheless, pseudobulges are difficult to identify, and there is not yet a consensus as to what defines them. To complicate things further, classical and pseudobulges can exist within the same galaxy (e.g., NGC 2787; [75]). The classification of the NGC 4736 bulge is ambiguous: Fisher and Drory [74] classifies it as a pseudobulge based on its nuclear spiral and bar and low Sérsic index (n=1.3). However, they also find it has a low star formation rate, more typical of classical bulges, and thus classify it as an “Inactive Pseudobulge;’’ if NGC 4736 hosts a pseudobulge, it is an atypical one. We have created a surface brightness profile from NICMOS and 2MASS H-band data from the Large Galaxy Atlas [76] and found results that conflict with the fits of Fisher and Drory [74]. Specifically, we find Sérsic indices for the bulge of n=2.3–3.0 depending on the radial range and type of fit (single versus double Sérsic), which are consistent with a classical one (as shown in [74]).Our surface brightness profile fits also reveal the presence of a nuclear star cluster within the central~0.6′′ (~12 pc), with an H-band magnitude of ~12.5. Such nuclear star clusters are common in early-type spiral galaxies [77]. The luminosity and mass of nuclear star clusters are known to scale with bulge luminosity and mass [78–80], and the NGC 4736 nuclear star cluster has a luminosity that is 0.1% of its bulge, typical for nuclear star clusters [81]. Nuclear star clusters commonly coexist with black holes, but there are a very limited number of cases where masses for both can be estimated [82]. In these cases, including the Milky Way, the BH mass is similar to the mass of the nuclear star cluster within an order of magnitude. There is also some evidence that the ratio of BH mass to nuclear star cluster mass increases with spheroid mass [83]. For NGC 4736, assuming an old population with an H-band M/L~0.7, the nuclear cluster would have a mass of ~2×107M⊙. This mass is quite similar to the M-σ* BH mass estimate, suggesting thus the presence of a similar-sized BH.


## 4.2. Caveats of the BLR Scaling Relation

The BLR scaling relations were derived using high-luminosity systems (i.e., Seyferts, not LINERs) and have been scaled to match the dynamical black hole detections of BHs with masses >106M⊙. Thus, their applicability to NGC 4736, where they yield an estimate of ~105M⊙, is a (perhaps unwise) extrapolation. These relations assume that the BLR is virialized due to proximity to the BH [84] and thus the 5100 Å continuum luminosity correlates with the emissivity-weighted radius of the BLR and thus with the BH mass [85]. The overall errors associated with BH mass measurements based on these relations do not exceed 0.5 dex [86–88], and there appears to be good consistency with BH masses obtained via the M-σ* relation for relatively bright AGN (L>1042 erg s-1), with MBH≳106M⊙ [89]. More recent calibrations of the radius-luminosity relationship on which these techniques are based infer that BH masses have been overestimated, however, only by up to a factor of ~3 [90]. This latter study also indicates a trend toward larger uncertainties and larger amount of overestimation in the BH with decreasing luminosity; however, no conclusive results are yet available for this regime.


## 4.3. Caveats of the MIR Line Correlation

The Dasyra et al. [65] empirical relation between the MIR line emission properties and the BH mass is derived using reverberation mapping BH masses and, thus, as with the previous method, may not apply to systems with MBH<106M⊙. The relation holds for systems with Lbol/LEdd>0.003, but not necessarily beyond this range. The MBH≈105M⊙ given by this method places this object at the low end of this Lbol/LEdd range; on the other hand, a more massive and thus quiescent BH would correspond to a [Ne V] luminosity of a few orders of magnitude lower than that measured. The scatter associated with this relation remains in average 0.5 dex and thus cannot account for the ~2-3 orders of magnitude difference between this value and that obtained via M-σ*.


## 4.4. Issues with Estimating the BH Mass from X-Rays

The BH mass derivation based on X-ray variability, or more precisely on the relation between the excess variance andMBH, relies on the hypothesis of a universal power spectral density function (PSD) shape and amplitude in AGN, which is based on the idea that the X-ray variability mechanism and the accretion efficiency are the same for all AGN, at all redshifts. These assumptions appear to hold for the objects involved in deriving and testing this relation [66]; that sample is, however, small and rather biased toward luminous (X-ray) sources, with log LX/(ergs-1)>41.5. None of the 2–10 keV sources detected in the nucleus of NGC 4736 are in this luminosity range. Nevertheless, Galactic BHs (Mbh<102M⊙) in their hard states show variability properties that match well those of AGN, both of these types of sources falling on the same projection of the TB-Mbh-Lbol/LEdd plane, where TB=1/νbf, suggesting that an extrapolation to the intermediate mass BH (or lower L) level is practicable.


## 4.5. Problems with the Fundamental Plane

Finally, it appears that the main problem with the fundamental plane is that the relation is not sharpened enough to provide strong constraints on BH masses. Because the relation spans nine orders of magnitude in BH mass, it is expected to equally apply to any value ofMbh in the range we are interested in. It is also the case that a wide variety of BH accretion models (e.g., with efficient and inefficient flows for the X-ray emission or associating X-ray flux with synchrotron emission near the base of a jet) are consistent with this relation (e.g., [91, 92]), suggesting that a large diversity of accretion modes or rates are accommodated.Nevertheless, the fundamental plane remains best constrained only for systems withMBH≳106M⊙ and for very low nuclear luminosities, that is, with negligible or zero AGN contribution that allow a dynamical measurement of their BH mass [68]. In this regime, the latest derivation of the fundamental plane relation provides for M94 a MBH value consistent with the M-σ* estimate and the dynamical measurement, while it would not reliably constrain a MBH≲106M⊙ value. The scatter in the fundamental plane relation increases for lower BH masses, with higher Eddington ratios [68, 93].


## 5. Discussion: Alternative Power Generation Mechanisms

The detection of broadHα emission, combined with the spatial coincidence of this emission with the detection of X-ray, UV, and radio compact sources, provides strong evidence that this galaxy hosts an accreting massive black hole. The SED of this nucleus makes this object one of the lowest luminosity LINER with a distinct contribution to the total emission by black hole accretion. In this scenario, the nucleus is an AGN, and the presence of off-nuclear sources, particularly the radio and the UV detections, may result from remnant jet activity emerging from the nucleus. It is definitely exciting to detect AGN activity at energy levels equal to that of several young supernova remnants of the Cas A variety [94], an OB association which hosts a high-mass X-ray binary [30] or simply a group of five late O supergiants [28].Nevertheless, the presence of a number of unusual off-nuclear sources, coupled with the apparent deficit in the photoionizing photons, and with the fact that the BH mass estimates based on AGN emission appear to fail for this object, encourages us to explore alternative scenarios for the NGC 4736 nucleus. In particular, it is very likely that the BH mass estimates that exploit the multiwavelength AGN characteristics do not work for this system because at least some of this emission is not the result of BH accretion.There are certain (peculiar) kinds of core-collapse SNe that present multiwavelength observations and in particular broad Hα emission components with characteristics that are very similar to those we measure in the nucleus of M94. The SN 2005ip, presented by Smith et al. [95], is a very good example. The so-called intermediate Hα component associated with this SN ejecta presents the same FWHM and brightness level as the broad Hα detected in the STIS aperture, and it does not show any sign of diminishing its strength over more than 3-year period (see their Figure  7; a very broad component with FWHM≳10000 km s-1 like the one exhibited by SN 2005ip would not be measurable in the galaxy spectrum as it would be completely swamped in the continuum stellar light, even in the HST-STIS observations).SNe as luminous as X2, with intrinsically hard X-ray spectra (photon indexΓ≲1), have certainly been encountered, for cases observed few years after the explosion (e.g., ATe #1023; [96]); for these cases, LX/Lbol ratios are high relative to more standard SN cases [97]. Another particular example of a SN that matches well the measurements of the nuclear emission in NGC 4736 is SDSSJ09529.56+214313.3 [98], which is believed (but not confirmed) to be a SN type IIn. For this system, the broad Hα component stays strong for at least three years, and its LX(2–10 keV)/LHα≈1, matching thus very well the surprisingly low value measured for NGC 4736.The peculiar nebular, radio, and UV characteristics of the nucleus of NGC 4736 appear to also compare well the SN phenomenology. The electron densities (or [S II] line flux ratios) measured in this galaxy center are matching exactly the ones measured in the environments of the extraordinary type IIn SNe we compare here with, for example, SDSSJ09529.56+214313.3. The two different band radio observations of this nucleus are consistent with emission from extragalactic SNe, in both intensity and decline rate in flux density [99–101]. Moreover, the compact radio off-nuclear detection, only 1′′(20 pc) away from the nuclear one [31], could be interpreted as the result of shock emission associated with a nuclear core-collapse SN. Measurements of the brightness of the off-nuclear UV detection in F250W and F330W place this object into the O star of late (5-ish) type spectral category, fitting thus well into the idea that this nucleus could simply be a star-forming site and that NGC 4736a’s emission includes significant contributions from a SN which exploded close to the weakly active (and massive, MBH~10(6-7)M⊙) central BH.This SN contribution scenario may have its own drawbacks. If the broad Hα emission has actually gotten brighter by a factor of 6, as suggested by the comparison between the HST-STIS and GMOS observations, the SN interpretation becomes problematic; variability of the broad-line region originating in the AGN could, in principle, account for this effect. Also, the observed intraday (hour-scale) X-ray variability measured for X2 remains yet to be detected in a supernova and conflicts with the physical scale over which X-ray emission is expected in SN remnants; thus, this particular behavior may remain strictly associated with the AGN. Additional observations would be necessary to fully confirm or rule out the SN scenario. Specifically, new high-resolution observations in the optical, UV, or X-rays would be able to confirm if there is a fading in the light curve, as expected from a SN, and would also allow accurate localization of the source. High S/N optical spectra would much better resolve the emission-line profile and a possible temporal evolution.


## 6. Conclusion

We have presented here an exhaustive multiwavelength analysis of the nuclear emission properties of NGC 4736, prompted by new measurements of a broad Hα emission component detected in its high-resolution HST-STIS optical spectrum. This broad Hα component, with a luminosity of 9×1037 ergs s-1, is one of the lowest luminosity broad line known. This broad Hα is coincident with a compact bright X-ray and radio source. Our measurements of this object’s spectral energy distribution reveal a bolometric luminosity of Lbol≈2.5×1040 erg s-1, that categorizes NGC 4736 as one of the least luminous LINERs with strong evidence for BH accretion. Our comparison of five independent BH mass estimates reveals a discrepancy of two orders of magnitude between the value ~107M⊙ predicted by the M-σ* relation and the value ~105M⊙ toward which methods based on AGN emission activity in optical, mid-IR, and X-ray, seem to converge; the fifth method is provided by the fundamental plane relation, which, however, due to its large associated uncertainties, does not offer any additional constraint to this comparison.We conclude that this system’s BH mass cannot be reliably estimated via standard AGN BH mass indicators because the nuclear emission in this system is not entirely tracing the accretion onto the central BH. Our assessment of the energy budgets of the ionizing and emission-line power suggests a possible deficit in the AGN ionization and production of a broad Hα emission feature which can be made up by a peculiar kind of Type IIn SN that matches well the nuclear emission of NGC 4736 over the whole electromagnetic spectrum and supports this galaxy nucleus’ general aging starburst-like appearance.



---

*Source: 178060-2012-02-07.xml*