Accretion disk versus jet orientation in H$_{2}$O megamaser galaxies

An essential part of the paradigm describing active galactic nuclei is the alignment between the radio jet and the associated rotation axis of the sub-pc sized accretion disks. Because of the small linear and angular scales involved, this alignment has not yet been checked in a sufficient number of Low Luminosity Active Galactic Nuclei (LLAGNs). The project intends to examine the validity of this paradigm by measuring the radio continuum on the same physical scale as the accretion disks, to investigate any possible connection between these disks and the radio continuum. We observed a sample of 18 LLAGNs in the 4.8 GHz (6 cm) radio continuum using the Very Long Baseline Array (VLBA) with 3.3 to 6.5 milliarcseconds resolution. The sources were selected to show both an edge-on accretion disk revealed by 22 GHz H$_{2}$O megamaser emission and signatures of a radio jet. Furthermore, the sources were previously detected in 33GHz radio continuum observations made with the Very Large Array. Five out of 18 galaxies observed were detected at 8 sigma or higher levels (Mrk0001, Mrk1210, Mrk1419, NGC2273 and UGC3193). While all these sources are known to have maser disks, four of these five sources exhibit a maser disk with known orientation. For all four sources, the radio continuum is misaligned relative to the rotation axis of the maser disk, but with a 99.1% confidence level, the orientations are not random and are confined to a cone within 32 degree of the maser disk's normal. Among the four sources the misalignment of the radio continuum with respect to the normal vector to the maser disk is smaller when the inner radius of the maser disk is larger. Furthermore, a correlation is observed between the 5 GHz VLBA radio continuum and the [OIII] luminosity and also with the H$_{2}$O maser disk's inner radius.


Introduction
Accretion of material onto the central engine of an Active Galactic Nucleus (AGN) is often accompanied by ejection of material, either in the form of outflows or (collimated) jets. Conservation of angular momentum implies that the distribution of any mate-always perpendicular to these disks (Schmitt et al. 2002). This misalignment indicates that gas on different scales has different angular momentum orientation. Whether this is a result of past mergers or warping of the accretion disk, or a general physical condition required for accretion, is still a matter of debate (see Greene et al. 2010, and references therein). With the discovery of H 2 O maser emission, emitted from the water vapor J K − K + =6 16 -5 23 rotational transition in the sub−pc disk surrounding the SMBH in NGC 4258 (e.g., Nakai et al. 1993;Greenhill et al. 1995;Miyoshi et al. 1995;Herrnstein et al. 1999), an opportunity arose to observe the accretion disk of AGNs with the extremely high angular resolution only offered by Very Long Baseline Interferometry (VLBI) observations. The other unique feature of galaxies with H 2 O masers in their nuclear accretion disk is that these accretion disks are viewed within a few degrees edge-on, otherwise linear scales for coherent amplification would be too short, and thus we expect the putative jets to be on the plane of the sky and free of viewing angle biases. Therefore, the nuclear region of galaxies hosting maser disks are good laboratories for studying the alignment between the jet and the normal to the pc-size disks. Greene et al. (2013) and Pjanka et al. (2017) investigated this alignment by observing the radio continuum jets on >50 pc scales. They also observed the structures such as galaxy-scale disks, bars and spiral arms on <5 kpc scales and compared the alignment of these structures with the smaller scale sub-pc maser disks. They reported that the orientation of the 100 pc scale disks and galaxy scale disks are distributed randomly with respect to the orientation of the maser disks. However, the radio continuum jets (observed on tens of pc to kpc scale) are found to align within < 15 • of the maser disk's rotation axis . With the purpose of investigating the jet-disk alignment on smaller scales (∼ 2 pc) and also providing a larger sample compared to that already existing in the literature, in this work we attempt to study the alignment of the jet and the rotation axis of the accretion disk in low luminosity active galactic nuclei (LLAGN) with H 2 O megamaser disks seen approximately edge-on. Recently we have observed a sample of 24 such galaxies with a resolution of 0. 2 -0. 5 using the Karl G. Jansky Very Large Array (VLA) at 33 GHz in a pilot project (Kamali et al. 2017). Radio synchrotron emission was detected on kpc scales in 21 of these galaxies. In the follow-up observations presented here, we look for radio emission on scales of <50 pc (∼5 mas or ∼2 pc resolution at a fiducial distance of 90 Mpc) using the Very Long Baseline Array (VLBA). We chose C-band because the synchrotron emission is supposedly brighter at this lower frequency compared to the previous VLA Ka-band observations. We aim at measuring the radio continuum and investigate the relation between the jet axis and the maser disk orientation. This paper is organized in the following way: In Sect. 2 we present our sample. In Sect. 3 we present the observations and data reduction process. The description of the observational results is presented in Sect. 4, followed by analysis and subsequent discussions in Sect. 5. A summary is given in Sect. 6.

The sample
We selected a sample of 18 LLAGNs with declinations > -5 • , including 14 sources from our 33 GHz survey (Kamali et al. 2017) plus four more from the Megamaser Cosmology Project (MCP, e.g., Reid et al. 2009Reid et al. , 2013Braatz et al. 2010;Greene et al. 2010;Kuo et al. 2011Kuo et al. , 2013Pesce et al. 2015;Gao et al. 2017). These sources have signatures of megamasers disks, i.e., they show three groups of H 2 O megamaser emissions in their spectrum; one set of maser features have velocities grouped around the recession velocity of the galaxy and those of the two other sets of maser features are offset from this recession velocity, representing the approaching and receding sides of the disk. In our selected sample, 12 sources are considered "clean disks", where the maser emission from the disk dominates over any other maser emission from the jets or outflows (Pesce et al. 2015). These sources are shown in boldface in all tables. The maser emission in six other sources of our sample may not arise from the disk, but rather from the nuclear outflows or star forming regions in these galaxies. Three out of these six non-clean disks are detected in this work. The maser distributions in two of these three galaxies (Mrk 0001 and Mrk 1210) have been mapped using the VLBI and the consequent fitting of disk models does not rule out a disk maser scenario (we refer the reader to Sect. 4.1 for further detail on the maser disk observations for the individual sources). Therefore, since it is unclear whether either of the aforementioned scenarios applies (the maser emission may actually arise from the disks), we consider all the sources as maser disks and our analysis is based on this assumption. The MCP is an ongoing project and the maser disk properties such as disk sizes and orientations are not yet measured for all the candidates.
For instance, the distances measured by the MCP (which are among the most accurate ones), are only available for four sources in our sample. These distances are in agreement within uncertainties with the distances listed in the NASA/IPAC Extragalactic Database (NED) 1 . For consistency we adopted the NED distances since they are available for all sources in our sample. The NED distances were obtained using H 0 =70.0 km s −1 Mpc −1 , Ω matter =0.27 and Ω vacuum =0.73 as the cosmological parameters. H 2 O maser luminosities from the literature were also rescaled to H 0 =70.0 km s −1 Mpc −1 to be consistent with other luminosities in this work.

Observations and data reduction
Our sample was observed between October and November 2015 using the VLBA in C-band. The synthesized beamwidth ranges from 3.3 to 6.5 mas. We used a total bandwidth of 128 MHz (2×64 MHz IFs, 4.644 to 4.708 GHz and 4.836 to 4.900 GHz) and right-hand circular polarization. The sources were grouped into 6 sets of three for observation, with a total time of 7 hours per group and an integration time of ∼1.5 hours per target after (i.e., excluding) calibration. Our observations included two scans of fringe finders, one placed at the beginning and one at the end of the observations. The sources and their associated calibrators were observed alternating in scans of ∼1 minute. Table 1 indicates fringe finders and calibrators, and the J2000 coordinates obtained from previous VLA observations. From the 18 observed sources, five were detected at signal-to-noise ratios (S/Ns) of 8 or higher levels. In the case of four of these detected sources, Mrk 0001, Mrk 1210, Mrk 1419 and UGC 3193, the initial images did not have a sufficiently high quality. Therefore, to reassure our detections and obtain higher S/Ns, we re-observed them with the VLBA in C-band, between October and November 2017. This time they were grouped in 2 sets of 2 sources, with a total time of 8 hours per group (2.5 hours on source per target, after calibration), bandwidth of 256 MHz (8×32 MHz IFs, from 4.836 to 5.092 GHz) and dual polarization.
The data were calibrated using the Astronomical Image Processing System (AIPS) software package developed by the Na- tional Radio Astronomy Observatory (NRAO). At first, the ionospheric dispersive delays and the Earth orientation parameters were corrected using AIPS tasks VLBATECR and VLBAEOPS respectively. Then, using task CLCOR, the parallactic angle cor-rection was applied. The visibilities were normalized by the auto-correlations with the task ACCOR. The instrumental phases and delays were removed using the fringe finders listed in Table 1 and task VLBAMPCL. Amplitude corrections were applied using system temperature and the antenna gain information by the APCAL task in AIPS. We solved phases, delays and rates by performing fringe fitting on the calibrators listed in Table 1 assuming a point source model, and applied the solutions to the targets. For imaging, we used the task IMAGR in AIPS. Different weighting functions (including tapering of the uv-distances) were tried to find the optimum solution for a better S/N.

Results
Out of 18 sources observed, 5 were detected: Mrk 0001, Mrk 1210, Mrk 1419, NGC 2273 and UGC 3193. Figure 1 shows the 5 GHz radio continuum maps. For all images, phase selfcalibration was employed to improve the image quality, using a solution interval of 1 minute. For Mrk 0001 the S/N was improved from 8 to 32, for Mrk 1210 from 36 to 62, for Mrk 1419 from 12 to 18, for NGC 2273 from 56 to 61 and for UGC 3193 from 8 to 36. For detected sources we fitted 2-dimensional Gaussians on the individual components using the task JMFIT in AIPS, in order to obtain source properties such as peak and integrated flux, position and its uncertainty (in case the continuum components are more complicated, multi-Gaussians are fitted over emission > 3 σ). The result is presented in Table 2. For detected sources, in case a source has more than one component, the flux reported in Table 2 is the total flux density which includes all the components 2 . For positions reported in Tables 1 2 The radio luminosities are also obtained from the total flux densities unless mentioned otherwise. and 2, the uncertainties include the statistical uncertainty, the uncertainty in the phase calibrator position or the calibrator's size in case the calibrator is resolved, and the systematic errors which are estimated from (Ψ 2 + beam 2 ) 1/2 × (1/[S/R] + 1/20), where Ψ is the source major axis, beam is the beam size, and S/R is the signal-to-noise ratio (White et al. 1997).

Continuum images and previous radio observations
Here we review our maps and other radio observations available in the literature for the five detected sources.

-Mrk 0001
We observed compact radio emission in our 33 GHz radio maps with a PA of 120 • ± 2 • . In the C-band VLBA data we see extended emission, also in the northwest-southeast direction. The maser disk has a PA of 6.4 • ± 10 • and its rotation axis is by 17 • ± 11 • misaligned from the jet propagation direction. The systemic masers are not observed in this source (Kuo et. al. in. prep.).

-Mrk 1210
As mentioned before, four of our VLBA detected sources were re-observed to obtain higher S/Ns. While for three sources morphologies are very similar in both observations, for Mrk 1210 the southeastern component was not detected in the data set from the first observations and the northwestern extended emission is less prominent in this data, likely due to insufficient sensitivities. A former 6 cm VLBA observation of this galaxy has shown a similar structure as our re-observations which is shown in   Kuo et al. (2018). Column 8: types of nuclear activity following NED; U stands for unidentified. Fig. 1, except that the northwestern component is resolved into four sub-components (Middelberg et al. 2004). The flux reported by Middelberg et al. (2004) for the central component in Fig. 1 is the same as in our observations (5 mJy). However the flux reported for the northwestern component is less than ours (9.1±1.1 mJy compared to 13.7±0.4 mJy). The southeastern component is considered a tentative detection in the work of Middelberg et al. (2004) and therefore no flux is reported for this component. Identifying the core of the galaxy from where the jet is emanating is not easy since both northwestern and central components have steep spectra (Middelberg et al. (2004) reported a spectral index of α 18 cm 6 cm =1.26 for the northwestern and α 18 cm 6 cm =0.78 for the central component obtained from VLBI observations with matching beam sizes using the S ∝ ν −α convention). The maser configuration is consistent with a a rotating disk. However, only one single maser component is detected at the systemic velocity of the galaxy. The high velocity masers also deviate slightly from an edge-on, flat, Keplerian disk model. This implies that if the masers arise from a disk, its inclination is deviating ∼10 • from an edge-on disk (inclination=100.68 +1.13 −0.90 • ), and there is potential warping or eccentricity in the disk (Zhao et al. 2018).
Mrk 1210 shows long-term variability in the X-rays and has changed from Compton-thick to Compton-thin over a period of  approximately 17 years (Masini et al. 2017). This variability was interpreted as a combination of the intrinsic variability of the AGN and presence of an obscuring cloud passing through the line of sight.

-Mrk 1419
This galaxy poses a clean maser disk (Kuo et al. 2011). Sun et al. (2013) reported extended emission with a PA of 125 • ± 10 • from 20 cm radio data. In our VLA maps we see a weak extension in the southeast-northwest direction with a PA of 117 • ±84 • . Gaussian deconvolution of our 5 GHz map shows a compact component that is not clearly resolved and may contain a weak extension with a PA of 159 • ±2 • . The outer part of the maser disk shows some degree of warping (Kuo et al. 2011).

-NGC 2273
This galaxy poses a clean maser disk (Kuo et al. 2011). NGC 2273 was previously observed with the VLBA at C-band in 2004 by Lal et al. (2004) and showed a compact radio core. Anderson & Ulvestad (2005) observed this galaxy at 8.4 GHz with the VLBA and detected radio emission extended along the northwest to southeast direction with a beam width of approximately 2 mas × 1 mas. The structure is similar to our map shown in the lower right panel of Fig. 1. The distribution of the masers shows a hint of a warp in the disk (Kuo et al. 2011).
-UGC 3193 In our 33 GHz data, we see extended emission in the north-south direction, the same as what we observe in the C-band VLBA data. The southern component is not resolved in the Gaussian deconvolution. The maser disk orientation and position is not yet known for UGC 3193.

A note on the non-detections
About 72% of the sources in our original sample of 18 sources were not detected. Considering that ∼70% of the undetected sources were previously detected in our VLA observations and assuming a flat spectrum or rising flux densities with decreasing frequency, one might expect to detect these sources in the VLBA observations. Figure 2 shows the expected flux at 5 GHz (assuming a flat spectrum) and the 5×rms noise of our VLBA maps (for a 5σ detection limit of a single spatially unresolved continuum component). The five detected sources are among the brightest ones as shown in Fig. 2. However, there are similarly bright sources among the undetected ones. Either these sources have inverted spectra, or extended emission is resolved out. Time variability of AGNs is another issue that could affect the detec-tion rates. We searched for trends of common properties among the detected sources, such as X-ray, IR, radio or maser luminosity limits, but did not succeed in finding any. We report 3σ upper limits for the undetected sources, rms of our maps, as well as hard X-ray and 33 GHz fluxes in Table 3.

Orientation of the jets with respect to the maser disks
In NGC 4258, the prototypical H 2 O megamaser galaxy, the rotation axis of the maser disk is aligned with the inner part of the extended jet (Cecil et al. 1992). Here we focus on identifying the orientation of the jet with respect to that of the maser disk in our sample. Once again we emphasize that two of the detected sources, Mrk 1419 and NGC 2273, are among the clean maser disks (see Sect. 2), while the other three detected galaxies, Mrk 0001, Mrk 1210 and UGC 3193, are among those exhibiting a the non-clean disk, i.e., the masers might arise from either disks or outflows. In the following discussion, we con-sider the possibility that the masers arise from the circumnuclear disk surrounding the SMBHs 3 . Then for four out of the five detected galaxies in our sample, Mrk 0001, Mrk 1210, Mrk 1419 and NGC 2273, the maser disk orientation is known (see Sect. 4 for more information). As mentioned before, in our sample the maser disk is viewed edge-on (±10 • ) and the normal vector to the disk is in the plane of the sky. We do not apply any inclination corrections to the measured position angles (PAs). Instead we assume that the PA (measured east of north) of the radio continuum shows the PA of the elongation of the jet. In case that more than one component is detected, we assume that a line going through the central position of each component (obtained from Gaussian fits) is representing the jet propagation direction. The same convention for measuring the PAs (east of north) is used for all other reported PAs. For all four sources, the jets are misaligned with respect to the normal vector to the maser disk. However they are confined to a cone within ∼0.6 radians (32 • ) of the disks' normal. Table 4 presents PAs obtained in this work and from the literature and Fig. 3 shows a comparison between the PAs of different scale disks and the radio continuum. For the null hypothesis that the jets have random distribution in space (which means they do not care about the maser disk), the probability that we observe these four measured PA offsets is 0.009, indicating that the jet orientations are not random with a 99.1% confidence level (see Appendix A for more details). On the other hand, the PA offset may be related to the outer and inner radii of the disk: the larger the disk size, the smaller the offset (see Figs. 4 and 5). We note, however, that our conclusion is based on a sample of four sources only (where two of them are non-clean disks), so that a study of a larger sample is needed to verify or falsify our finding.

Multi scale position angle
The PAs of disks on different scales are available for three of the VLBA-detected sources Pjanka et al. 2017). Therefore we can compare the alignment of the jet with the angular momentum orientation of the disks on different scales. Our VLBA radio continuum in Mrk 1210 is misaligned by ∼26 • with the maser disk's normal. However, it is misaligned with the nuclear region of the galaxy (with size of 170 pc) by only ∼8 • . The radio continuum orientation is displaced by ∼57 • from the angular momentum orientation of the large scale disk. For NGC 2273, the VLBA radio continuum is not aligned with any of the rotation axes: ∼32 • from the maser disk's normal, ∼36 • from the angular momentum orientation of the nuclear 150 pc scale disk and ∼58 • from the normal vector to the galaxy's large scale disk. In Mrk 1419, which has likely experienced a merger or some interaction (Läsker et al. 2016), the maser disk and the nuclear region with a size of 270 pc are aligned and the jet is ∼20 • from the angular momentum direction of this disk, while the galaxy's large scale disk has a normal vector that is ∼29 • away from the jet direction.

Radio continuum positions relative to the maser positions
Compact radio emission is detected in the nuclei of most galaxies (Herrnstein et al. 1997, and references therein). Since the centers of the maser disks are the gravitational centers of the nuclear regions, we expect the centers of the maser disks to coincide with their radio continuum core representing the nucleus of the galaxy. For three sources among the detected ones, Mrk 0001, Mrk 1419 and NGC2273, the maser positions are available. For our radio continuum data the central position from the Gaussian fit of the component with the highest brightness temperature (if there are more than one) is considered as the position of the nuclear source (see Sect.5.4). In Table 5 we present the spatial offsets of the maser disk with respect to the radio continuum. For Mrk 1419 and NGC 2273, the maser positions are from VLA observations. Within the uncertainties, the maser positions for all three sources are consistent with those of the radio continuum. For Mrk 0001, however, in our second observations the phase calibrator was resolved with a major axis of ∼ 15 mas and therefore we report the position offset with respect to the continuum position in our first observations which has a higher position accuracy. The proper motion of the galaxies can be neglected. In case of a speed of 1000 km s −1 in the plane of the sky and a distance D = 50 Mpc, the expected shift is only 4 µarcsec yr −1 .

Spectral indices and brightness temperatures
We obtained the spectral indices for the detected sources between 33 GHz and 5 GHz. They are presented in Table 2. The indices were obtained assuming a power-law dependence for the fluxes (S ∝ ν −α ). We should note that the VLA observations have linear beam sizes about 100 times those of the VLBA observations. The spectral indices have a mean of α=0.21±0.20.
We also obtained the brightness temperatures (T B ), using where k B is the Boltzmann constant, S tot is the total flux density, λ is the wavelength of observation, and θ ma j and θ min are the deconvolved sizes of the major and minor axes of the source. When the component is unresolved, the beam size is used as an upperlimit for the source size. Therefore the corresponding T B is a lower limit. The obtained brightness temperatures are reported in Table 2. For sources with more than one component, we report the measurement for the component with highest T b , which we identify as the putative core of the AGN. It should be noted that we do not know the spectral indices for the individual components, therefore the brightness temperatures are the only tool by which we can tentatively define the AGN's core. All sources have T B > 10 6 , consistent with non-thermal emission.
5.5. Radio luminosity versus other properties of the galaxies 5.5.1. Radio luminosity vs. X-ray luminosity and the SMBH mass: fundamental plane of black hole activity The fundamental plane (FP) of black hole activity is an empirical relation between the radio continuum luminosity, X-ray luminosity and black hole (BH) mass. This relation holds for 10 orders of magnitude in radio luminosity and ranges from stellar mass BHs (in X-ray binaries) to SMBHs, supporting the idea that jet physics is scale invariant (Merloni et al. 2003). For many years there have been efforts to understand the accretion mechanism in the SMBHs and how the accretion rate is related to the SMBH mass and the radio and X-ray luminosities. Studies have suggested that the radio luminosity depends on both SMBH mass and accretion rate of material onto the black hole (Falcke & Biermann 1995). The SMBH mass can be determined observationally, and the accretion rate can be inferred through measuring   the luminosity at those frequencies at which the accretion process dominates. For example, the accretion flow can be inferred from the X-ray luminosity since the X-ray emission is produced by the accretion of material onto the SMBHs. The scaling of the X-ray luminosity with the accretion rate depends on the physics of the disk, e.g., whether it is radiatively efficient or not (a quan-titative description of these models is beyond the scope of this paper; see Merloni et al. 2003, and references therein). The Xray emission could also be produced by the inverse Compton effect involving synchrotron radiation. Summarizing everything, we expect the core radio luminosity to be correlated with the X-ray luminosity and also with the black hole mass.  Table 4). The disk sizes are to scale with respect to each other. , assuming a constant X-ray luminosity (mean value of our data) and the solid lines present ξ M obtained in this work. Right panel: 5 GHz radio luminosity versus X-ray luminosity. Dashed lines present ξ X available in the literature, assuming a constant black hole mass (mean value of our data) and the solid lines present ξ X obtained in this work.
The following plane equation describes the FP of BH activity: where L R , L X and M are the radio luminosity, X-ray luminosity and the BH mass respectively, b R is the zero intercept, and ξ X and ξ M are the linear regression coefficients. The coefficients in Eq. 2 depend on the accretion efficiency, electron spectral index, viscosity of the disk and initial conditions at the base of the jet such as dependence of magnetic energy density on the SMBH mass and accretion rate. Others have followed the lead of Merloni et al. (2003) and investigated the FP using different subsamples (e.g., Saikia et al. 2015;Körding et al. 2006;Gültekin et al. 2009;Bonchi et al. 2013;Nisbet & Best 2016). All these studies, fit similar coefficients, which is somewhat surprising since the SMBH mass determinations differ from sample to sample. We note that a sample for which accurate BH mass measurements have been obtained from edge-on megamaser disks has not yet been analyzed.
Adopting X-ray luminosities from Litzinger et al. (in. prep), obtained using data from the Swift/BAT instrument in a range of 20-100 keV (see Table 3), the black hole masses from the MCP and radio luminosities from this work, we fitted the same function to our data using the least square method. The fitting parameters are shown in the first row of Table 6. We considered increasing the number of sources for this analysis by extrapolating the 5 GHz luminosities from the 33 GHz radio luminosities, using the mean of the indices between our 33 GHz and 5 GHz observations (see Table. 2). This sample with extrapolated fluxes includes 12 sources that have known SMBH masses. The result  (1) Kuo et al. in prep.; (2) Kuo et al. (2011).
of the fit is presented in the second row of Table 6. We should note that the mentioned studies that probe the FP use soft X-ray (2-10 keV) data and this can cause a discrepancy to our result. We also tried to fit Eq.2 to a sample of 5 sources for which we have soft X-ray data (Masini et al. 2016) and extrapolated 5 GHz fluxes. The fitting parameters are shown in the third row of Table 6. We present the fitting parameters for the FP in Fig. 6. A selection of parameters reported in the literature is also presented with dashed lines. As seen in Fig. 6, in our sample the black hole mass decreases with increasing radio luminosity. This is opposite to what is found in other studies and could be responsible for the significant difference between this work and other studies for the obtained value of ξ M . On the other hand, galaxies with H 2 O megamaser disks may not follow the well-known scaling relation between the BH mass and the stellar dispersion velocity of the galaxies. They may show an offset with respect to the early-type galaxies, while other spiral galaxies with non-maser dynamically determined BH masses do not show such an offset ). This could be because the maser disks are preferentially formed in spirals with particularly low mass nuclear sources, or because the non-maser dynamical measurements are missing the low-mass end of the BH distribution due to their too small spheres of influence. It should be noted that these two above mentioned interpretations are only two possible scenarios among other possible ones.

Radio luminosity vs. [OIII] luminosity
Previous studies have shown a strong correlation between the optical narrow emission-line luminosity and the radio luminosity of radio galaxies (e.g., Baum & Heckman 1989;Rawlings & Saunders 1991). This correlation is valid over four to five orders of magnitude for both radio and line luminosities and could be the result of a common energy source for optical line and radio emissions, suggesting that the luminosity depends on the properties of the central engine. It is also possible that both optical line and radio luminosities depend on a third parameter, e.g., the amount of cold gas present on kpc scales (Baum & Heckman 1989). A correlation was reported by Baum & Heckman (1989) between the Hα+[N II] line luminosity, L line , and total radio luminosity at the VLA scales for a sample of radio galaxies where L line is proportional to L 0.73 radio . Other studies have compared the radio core luminosity and the luminosities from the narrow line region, where the core is defined as the unresolved component with highest flux density in the radio observations. Buttiglione et al. (2010) reported L [OIII] ∝ L 0.75 core for a sample of high excitation galaxies and Baldi et al. (2018) reported L [OIII] ∝ L 0.35±0.20 core for a sample of LINER galaxies. Therefore, we check for a possible linear correlation between the optical emission-line luminosities and radio luminosities in our sample. We obtained the 5007 Å[OIII] luminosities from NED taken by the SDSS. Four of our five detected sources, Mrk 0001, Mrk 1210, Mrk 1419 and NGC 2273, have [OIII] luminosities. Using a Spearman rank correlation test, we find a correlation between 5 GHz (derived from VLBI) and the [OIII] luminosities, L [OIII] ∝ L 0.87±0.09 5 GHz . The correlation coefficient is 0.80, and the P-value, indicating that the likelihood that parameters are unrelated, is 0.2 (see Fig. 7). While this suggests reasonable agreement with previously mentioned studies, we note the result is based on very few sources only. We also note that the correlation of radio luminosity with emission-line luminosity extends over five orders of magnitude in L radio and four orders of magnitude in L line , i.e., there is ∼ one order of magnitude scatter in the radio luminosity for a given line luminosity, suggesting that other criteria such as environment or different accretion modes play a role in determining the radio and optical emission-line luminosities of these galaxies (Baum & Heckman 1989

Radio luminosity vs. maser disk radius
In our pilot project (Kamali et al. 2017) we demonstrated a correlation between the disk's inner and outer radii and the 33 GHz radio luminosity of galaxies with H 2 O masers in their nuclear disks (i.e., both inner and outer radii of the disk increase as the radio luminosity increases). In this work, we find a similar correlation for inner and outer disk radius versus the 5 GHz luminosities (Fig. 8). Correlations between the maser disk size and luminosity at other frequencies have been investigated formerly. For example, Gao et al. (2017) reported a correlation between the infrared (IR) luminosity and the maser disk's outer radius. Notes. Column 1: coefficient of X-ray luminosity in Eq.2. Column 2: coefficient of SMBH masses in Eq.2. Column 3: intercept of Eq.2. Column 4: notes. Extrapolated 5 GHz fluxes are obtained using spectral index 0.21 between 33 GHz VLA observation and 5 GHz VLBA observations assuming a power law dependence of S ∝ ν −α .
This correlation exists because the gas temperature should be more than the minimum kinetic temperature of ∼ 400 K for masing to occur. The H 2 O molecules are mixed with dust and the dust temperature depends on the bolometric luminosity. Therefore the outer radius of the maser disk also depends on the bolometric luminosity of the AGN. In Sect. 5.5.2 we showed a correlation between the 5 GHz luminosity and the [OIII] luminosity (which is a proxy for the AGN's bolometric luminosity, e.g., Bassani et al. 1999;Heckman et al. 2005). Thus a correlation between the outer radius of the maser disk and 5 GHz luminosities is expected. Furthermore, the pressure outside the outer disk is less than the critical pressure which permits the existence of the molecules (Neufeld et al. 1994) and this critical pressure depends on the X-ray luminosity (for detailed analysis see Neufeld & Maloney 1995). Although a correlation between the X-ray luminosity and the maser disk size is expected, former studies did not find such a correlation (Kamali et al. 2017). High energy electrons are responsible for both synchrotron radio emission and the synchrotron self Compton X-ray radiation (Falcke & Biermann 1995). However, our sources are Compton thick and a misestimate of the column densities can affect the estimation of the intrinsic unobscured X-ray luminosities. Radio emission does not suffer from extinction and might be a better tool to estimate the impact of these high energy electrons that are confining the disk.

Summary
We observed a sample of 18 LLAGNs with the VLBA at C-band (6 cm) and detected 5 sources: Mrk 0001, Mrk 1210, Mrk 1419, NGC 2273 and UGC 3193. We aimed to measure the orientations of the putative jets with respect to those of the maser disks. Excecpt for UGC 3193, for the other four detected sources, the maser disk orientation is known. All four sources have jet directions that are misaligned with the disk's normal by < 0.6 radian.
For the null hypothesis that the jets have random distribution in space, the probability that all four sources have a misalignment of < 0.6 radian is 0.009, indicating that the jet orientations are not random. The PA differences seem to be related to the disk size: the larger the disk size, the more aligned is the radio continuum to the disk's normal. A tight correlation with correlation coefficient of 0.87 is observed between the 5007 Å[OIII] line luminosity and the 5 GHz radio continuum luminosity. The disk's inner and outer radii show a linear correlation with the 5 GHz continuum luminosity: the larger the inner radius of the maser disk the brighter is the central region of the galaxy at 5 GHz. Note, however, that all these results are based on a very small sample of sources. Further observations and maps of maser disks and of the nuclear regions of maser host galaxies will provide a larger sample for a more detailed investigation on the alignment of the accretion disk and jets in the galaxies where both phenomena are directly observable.