\section{Introduction}

The discovery of {\it{Hanny's Voorwerp}} (HV hereafter) by a citizen scientist taking part in the Galaxy Zoo project (\citealt{Lintott2008,Lintott2009}) opened the field for studying Active Galactic Nuclei (AGN) variability on $10^{4}-10^{5}$ yr time-scales.
HV is an extended emission line region ($11 \times 16$ kpc in projected extent) located at $\sim 20$ kpc in projection from the core of IC\,2497, a massive ($M_* \sim 10^{11} M_{\odot}$), nearby ($z = 0.0502$) postmerger galaxy. The optical spectrum and the Sloan Digital Sky Survey (SDSS) $g$-band luminosity of HV are dominated by [$\mathrm{O}$\textsc{iii}]$\lambda 5007$ emission. The presence of emission lines of high-ionization species such as [$\mathrm{Ne}$\textsc{v}]$\lambda \lambda 3346, 3426$ and [\heii]$\lambda4616$, together with the relatively quiet kinematics (line widths $<100$ km s$^{-1}$) and the low electron temperature ($T_{e} \sim 13,500$K) suggests that the gas is photoionized by the hard continuum of an AGN in IC\,2497 rather than ionization from star formation or shocks \citep{Lintott2009}.
In order to produce sufficient ionizing photons to power the observed emission, the AGN output should be at least $L_{\rm ion} \sim 10^{45}$\,erg\,s$^{-1}$ between 1 and 4 Ryd (13.6-54.4 eV), corresponding to a bolometric luminosity $L_{\rm bol} \sim 10^{46}$\,erg\,s$^{-1}$.

Although the photoionization state of HV requires the presence of a strong AGN, the optical spectrum of the core of IC\,2497 only shows weak optical emission lines (\citealt{Lintott2009}; \citealt{Keel2012b}) and a weak radio source ($\sim10^{38}$\,erg\,s$^{-1}$; \citealt{Jozsa2009}). The two possible explanations to reconcile the observed and expected AGN output are: (1) the dust in the central part of the galaxy is arranged in such a way that the quasar appears obscured only along our line of sight but not in the direction of HV, or (2) the AGN dropped in luminosity in the last $\sim 100$ kyr, the travel time needed for the ionizing photons to reach the cloud \textcolor{black}{{\rm{(we note that this time is a lower limit since the measured distance is a projected distance)}}}. The first scenario was favoured by \cite{Jozsa2009}: using multifrequency radio observations they found that HV is embedded in a massive $\sim 10^9 M_{\odot}$ $\mathrm{H}$\textsc{i} reservoir crossed by an outflow connecting the centre of the galaxy to HV (see also \citealt{Rampadarath2010}), and derived a lower limit for the nuclear $\mathrm{H}$\textsc{i} column density of $2.8 \pm 0.4 \times 10^{21}$ cm$^{-2}$. Based on these data they suggested a scenario where the AGN is strongly obscured and the outflow cleared a path for the ionizing radiation to escape in the direction of HV.
On the other hand, the fading scenario was supported by \cite{Lintott2009}, \cite{Schawinski2010b}, \cite{Keel2012b}, \cite{Sartori2016} and \cite{Keel2017}. First, the IR data from {\it{IRAS}} show no evidence of the reprocessed luminosity expected from an obscured quasar of the required luminosity. Secondly, the available soft X-ray {\it{Chandra}} \, and {\it{XMM-Newton}} \, spectra are well fitted by an unabsorbed power law plus a component for the hot gas in the galaxy, and the quality of the fit is not improved by the addition of an obscured AGN component. Finally, {\it{Hubble Space Telescope}} ({\it{HST}}) imaging shows a complex dust structure in the centre of IC\,2497, but the view of the nucleus is not hindered by the detected dust lanes.

Despite intensive efforts to measure it, the intrinsic current luminosity of the AGN in IC\,2497 is still not known. In this paper we present new hard X-ray observations obtained with the {\it{Nuclear Spectroscopic Telescope Array}} ({\it{NuSTAR}}; \citealt{Harrison2013}). Due to the high sensitivity and the hard X-ray coverage of {\it{NuSTAR}} \, ($3-79$ keV) we are able, for the first time, to put strong constraints on the current intrinsic AGN luminosity and to test directly the fading and obscuration scenarios. \textcolor{black}{{\rm{The paper is organized as follows. In section 2 we describe the data and data reduction. In section 3 we describe the performed analysis and how we estimate the past and current AGN bolometric luminosity. The results are discussed in Section 4.}}}

\section{Observations}

\begin{table*}
 \begin{tabular}{lcccc}
  \hline
  \hline
  \noalign{\smallskip}
   Obs. \# & \textcolor{black}{\rm{Observatory}} & Observation date & Observation ID & Net exposure \\
   \noalign{\smallskip}
  \hline
  \noalign{\smallskip}
  1 & {\it{Chandra}} \, & 2012 Jan 8 & 13966 & 59.35 ks \\
  \noalign{\smallskip}
  2 & {\it{Chandra}} \, & 2012 Jan 11 & 14381 & 53.4 ks \\
  \noalign{\smallskip}
  3 & {\it{NuSTAR}} \, & 2016 Sep 29 & 60201017002 & 31.4 ks \\
  \noalign{\smallskip}
  \hline
 \end{tabular}
 \caption{X-ray observations log for the data analysed in this paper. The archival {\it{Chandra}} \, observations cover the soft part of the X-ray spectrum ($0.5-8.0$ keV) while the newly obtained {\it{NuSTAR}} \, observations probe the hard part of the spectrum ($3-24$ keV).}
 \label{tab:log}
\end{table*}

We present the analysis of new {\it{NuSTAR}} \, observations of IC\,2497 ($3-24$ keV, \textcolor{black}{\rm{see Section \ref{sec:NuSTAR}}}), together with previously obtained {\it{Chandra}} \, X-ray observations of the galaxy ($0.5-8.0$ keV), and optical longslit spectroscopy \textcolor{black}{{\rm{and narrow-band imaging}}} of both the galaxy and HV. Details about the {\it{Chandra}} \, and {\it{NuSTAR}} \, observations are listed in Table\,\ref{tab:log}.

\subsection{Chandra}
{\it{Chandra}} \, \citep{Weisskopf2000} observations of IC\,2497 were obtained with the Advanced CCD Imaging Spectrometer S-array (ACIS-S; \citealt{Garmire2003}) in very faint (VFAINT) time-exposure (TE) mode on 2012 January 8 (ObsID 13966, 60 ks) and 2012 January 11 (ObsID 14381, 55 ks; PI K. Schawinski, Cycle 13). We performed standard data reduction starting from the level 1 event files using the {\textsc{CIAO}} \, 4.7 software \citep{Fruscione2006} provided by the {\it{Chandra}} \, X-ray Center (CXC). For both exposures we extracted the spectra and generated the response curves from a circular aperture with a 4 arcsec radius around the centroid of the emission using the {\textsc{CIAO}} \, task \textsc{specextract}. For the background we considered three source-free apertures $\sim 20$\,arcsec away from the target. We grouped the spectra with a minimum of 3 photons per bin using the {\textsc{HEAsoft}} \, task \textsc{grppha}. Details about the {\it{Chandra}} \, data, data reduction and spectral extraction are provided in \cite{Sartori2016}.

\subsection{NuSTAR} \label{sec:NuSTAR}
{\it{NuSTAR}} \, observations of IC\,2497 were obtained on 2016 September 29 (ID 60201017002, 39 ks; program 02041, PI L. Sartori, GO cycle-2). We processed the raw event files using the {\it{NuSTAR}} \,  Data Analysis Software package ({\textsc{NuSTARDAS}} \, version 1.6.0)\footnote{\url{http://heasarc.gsfc.nasa.gov/docs/nustar/analysis/nustar_swguide.pdf}} and the calibration files from \textsc{CALDB} version 20160922. Because of the slightly elevated background event rates around the SAA we processed the data with the parameters \textsc{SAAMODE=strict} and  \textsc{TENTACLE=yes} as recommended by the {\it{NuSTAR}} \, Science Operations Center (SOC). We extracted the source and background spectra using the \textsc{nuproducts} task included in the {\textsc{NuSTARDAS}} \, package with the appropriate response and ancillary files. For each focal plane module we extracted the source spectra using a circular aperture with 45 arcsec radius. We checked that no additional bright sources are detected by {\it{Chandra}} \, within this region. For the background we considered \textcolor{black}{{\rm{two source-free apertures with 60 arcsec radius on the same detector chip as the source}}}. Similar to the {\it{Chandra}} \, spectra, we grouped the {\it{NuSTAR}} \, spectra with a minimum of 3 photons per bin using the {\textsc{HEAsoft}} \, task \textsc{grppha}. The obtained net counts and background counts are listed in Table \ref{tab:nu}. \textcolor{black}{\rm{For the analysis we limited the {\it{NuSTAR}} \, spectrum to 3-24 keV instead of 3-79 keV because of the drop in sensitivity and the heightened background above 24 keV (e.g. \citealt{Wik2014}).}}

\begin{table*}
 \begin{tabular}{lccc}
  \hline
  \hline
  \noalign{\smallskip}
    & FPMA & FPMB & FPMA + FPMB \\
   \noalign{\smallskip}
  \hline
  \noalign{\smallskip}
  Source net counts & $59.9^{+12.4}_{-11.4}$ & $84.6^{+13.2}_{-12.2}$ & $144.5^{+17.7}_{-16.7}$ \\
  \noalign{\smallskip}
  Background counts & $63.1^{+16.6}_{-15.9}$ & $58.4^{+17.8}_{-17.1}$ & $121.5^{+24.1}_{-23.3}$  \\
  \noalign{\smallskip}
  \hline
 \end{tabular}
 \caption{Source net counts and background counts in the 8-24 keV band, for each {\it{NuSTAR}} \, focal plane.}
 \label{tab:nu}
\end{table*}

\subsection{Optical data}

Optical longslit spectra passing through both the nucleus of IC\,2497 and HV were obtained with the double-spectrograph system at the 4.2m William Herschel Telescope (WHT) on La Palma and with the 3m Shane telescope at the Lick Observatory. All the spectra were obtained with a 2\,arcsec slit width and nearly identical position on the sky. Details about observations and data reduction are given in \cite{Lintott2009}.

\textcolor{black}{{\rm{Narrow-band H$\alpha$ images are obtained with the {\it{HST}} using the tunable ramp filters on the Advanced Camera for Surveys (ACS). Details about the data  and data processing are described in \cite{Keel2012b} and \cite{Keel2017}.}}}

\section{Analysis}

The goal of this analysis is to measure the level of obscuration and the current intrinsic bolometric luminosity of the AGN in IC\,2497. A comparison of the obtained current luminosity with the past luminosity inferred from the photoionization state of HV will allow us to probe if the quasar really dropped in luminosity within the last $\sim 100$ kyr as suggested by previous studies, and the magnitude of the drop.

\subsection{X-ray emission}
We fit the obtained spectra using {\textsc{Xspec}} \citep{Arnaud1996} with Cash statistics \citep{Gehrels1986}, and estimate best-fitting parameters and errors using the MCMC routine implemented in {\textsc{Xspec}}. Because of the low signal to noise of the {\it{NuSTAR}} \, observations, we first examine the {\it{Chandra}} \, observations alone. We then perform combined fits considering different models. We verified that no significant variability is present between the observations but allow a $10\%$ offset to account for different responses of the detectors and uncertainties in the cross-calibrations. Errors for the best-fitting parameters are quoted at the $90\%$ confidence level.

\begin{figure*}
\includegraphics[scale=0.5]{spec_v2_2.pdf}
\caption{Broadband X-ray spectrum from {\it{Chandra}} \, (green) and {\it{NuSTAR}} \, (blue). The crosses correspond to the data while the lines show the models. For plotting clarity we rebinned the spectrum with a minimum of 9 counts per bin. {\it{Left}}: The fit corresponds to the best fit to the {\it{Chandra}} \, data only, an unabsorbed power law with an APEC component, with extrapolation at higher energies. This model clearly underpredicts the hard X-ray counts observed with {\it{NuSTAR}}. {\it{Right}}: The fit corresponds to the \textsc{MYTorus} model with $\theta_{\rm i} = 90^{\rm o}$. This model is a better match to the data, and indicates that the source is \textcolor{black}{\rm{obscured. The box shows a zoom-in on the iron \Ka lline region with a different binning (minimum of 3 counts per bin)}}.}
\label{fig:spec}
\end{figure*}

\subsubsection{{\it{Chandra}} \, observations}

The $0.5-8.0$\,keV {\it{Chandra}} \, spectrum is best fitted (\textcolor{black}{\rm{C-stat/dof = 181.9/132}}) by an unabsorbed power law ($\Gamma = 2.14^{+0.27}_{-0.28}$), as expected for an unobscured AGN or for the Thomson-scattered emission of an obscured AGN, plus an APEC component representing the hot gas in the galaxy ($kT = 0.93^{+0.07}_{-0.11}$ keV). \textcolor{black}{\rm{We notice that the power law component may also be partially attributed to emission from star formation and X-ray binaries, as discussed in Section \ref{sec:SF}}}. The total absorption component is fixed to the Galactic value ($N_{\rm H} = 1.31 \times 10^{20}$\,cm$^{-2}$, as given by the Colden Calculator with the NRAO data set\footnote{\url{http://cxc.harvard.edu/toolkit/colden.jsp}}). A similar value for the photon index ($\Gamma = 2.23^{+0.57}_{-0.38}$) is obtained also when  the counts below 1.5 keV in order to avoid contamination by the hot gas. By fitting the $3.0-8.0$\,keV spectrum with a simple power law we obtain a photon index $\Gamma = -0.76^{+1.45}_{-0.05}$. This low value is an indication that the source is heavily obscured and that the power law observed in the soft X-ray spectrum corresponds to the scattered component. However, due to the narrow energy interval considered, adding an obscured component to our fiducial model for the {\it{Chandra}} \, data does not improve the quality of the fit significantly, so that a combined fit with {\it{NuSTAR}} \, is needed to determine the real amount of obscuration.

\subsubsection{Possible contribution from star formation}\label{sec:SF}

\textcolor{black}{\rm{IC 2497 is a luminous infrared galaxy (LIRG) with far-IR luminosity $L_{\rm FIR} = 6.6 \times 10^{44}$\,erg\,s$^{-1}$}. Assuming the FIR - star formation rate relation of \cite{Kennicutt1998}, this corresponds to a star formation rate $SFR \sim 30 M_{\odot}$\,yr$^{-1}$. Using the relations of \cite{Ranalli2003} or \cite{Lehmer2016}, the expected $2-10$ keV luminosity is $L_{\rm 2-10 keV} \sim 1.5 \times 10^{41}$ erg s$^{-1}$, which is comparable (within the large scatter of the correlations) to the observed $2-10$ keV luminosity measured from the {\it{Chandra}} \, and {\it{NuSTAR}} \, spectra (we stress that all the used relations have a large scatter, so that the results have to be considered as an order of magnitude estimation). In addition, the low {\it{W1}}-{\it{W2}} colour, {\it{W1}}-{\it{W2}} = 0.3, suggests that the mid-IR emission is dominated by stellar emission (see Section \ref{sec:mir} for a discussion of the mid-IR emission). Altogether, this points to a scenario where star formation provides an important contribution to the soft part of the X-ray spectrum. We will consider this component in the analysis below.  }

\subsubsection{Phenomenological model} \label{sec:PE}

Applying the unabsorbed power law model described above to the $0.5-24.0$\,keV spectrum significantly underpredicts the hard X-ray counts observed with {\it{NuSTAR}}, especially in the $10.0-20.0$\,keV range \textcolor{black}{\rm{(see Fig. \ref{fig:spec}, left)}}. This again is an indication that the source is highly obscured (e.g. \citealt{Koss2016}; \citealt{Ricci2016}). To infer the intrinsic column density $N_{\rm H}$ and the X-ray luminosity of the source we first fit the $0.5-24.0$\,keV spectrum with a phenomenological model which includes 1) an absorbed power law representing the intrinsic emission of the AGN, where we consider both photoelectric absorption and Compton scattering, 2) \textcolor{black}{\rm{a second unabsorbed power law representing the scattered AGN component and the contribution from star formation}}, 3) the unabsorbed reflected X-ray emission from neutral material above the accretion disc (the PEXRAV model; \citealt{Magdziarz1995}), 4) a Gaussian representing the narrow \Ka lline usually detected in AGN (e.g. \citealt{Nandra1994}; \citealt{Shu2010}; \citealt{Ricci2014}), with fixed energy $E = 6.4$keV and sigma $\sigma = 10$eV, and 5) an APEC model for the hot gas in the galaxy: \\

\textsc{phabs$_{\rm GAL}$} $\times$ (\textsc{zphabs} $\times$ \textsc{cabs} $\times$ \textsc{zpow} $+$ \textsc{pexrav} $+$ \textsc{zpow} $+$ \textsc{zgauss} $+$ \textsc{apec}). \\

\textcolor{black}{\rm{The photon index of the intrinsic power law is fixed at $\Gamma = 1.9$ as expected for an average AGN (e.g. \citealt{Nandra1994}; \citealt{Piconcelli2005}; \citealt{Ricci2015}), while the photon index of the unabsorbed power law is let free to vary. The best fit results in a column density $N_{\rm H} = 1.87^{+1.27}_{-0.21} \times 10^{24}$ cm$^{-2}$ and an intrinsic 2-10 keV luminosity (intended as the luminosity of the intrinsic power law) $L_{\rm 2-10 keV} = 2.3 \times 10^{43}$ erg s$^{-1}$. We note that the reflection parameter $R$ is an upper limit and the best-fitting parameters do not change significantly if we fix $R = -1 \times 10^{-10}$. This suggests that the reflection component is not statistically significant.}}

In order to estimate the equivalent width (EW) of the narrow \Ka lline we fit the {\it{Chandra}} \, and {\it{NuSTAR}} \, data in the $3.0-24.0$\,keV range simultaneously with a simplified model: \textsc{zpow} $+$ \textsc{zgauss}. We obtained EW $ 1.74^{+0.86}_{-0.72}$ keV, which is again consistent with heavy obscuration (e.g. \citealt{Matt1996}).

\subsubsection{BNtorus model} \label{sec:BN}

A more physically motivated model for highly obscured sources is the \textsc{BNtorus} model developed by \cite{Brightman2011}, which considers both absorbed and reprocessed emission of an intrinsic power law from a torus. In order to fit the whole $0.5-24.0$\,keV spectrum simultaneously we add a scattered power law component with same photon index as the intrinsic one and an APEC model in a similar way as described above: \\

\textsc{phabs$_{\rm GAL}$} $\times$ (\textcolor{black}{\textsc{BNtorus}} $+$ \textsc{zpow} $+$ \textsc{apec}).\\

In order to obtain a good fit we need to fix the inclination angle of the torus to the maximum permitted value $\theta_{\rm i} = 87.1^{\rm o}$. Varying the opening angle of the torus between $\theta_{\rm oa} = 40^{\rm o}$ and $\theta_{\rm oa} = 80^{\rm o}$ provides fit parameters which are consistent within the uncertainties and statistically indistinguishable. For $\theta_{\rm oa} = 60^{\rm o}$ the fit implies \textcolor{black}{\rm{$N_{\rm H} = 2.10^{+1.05}_{-0.57} \times 10^{24}$ $\rm {cm}^{-2}$}} and \textcolor{black}{\rm{$L_{\rm 2-10 keV} = 1.0 \times 10^{43}$ erg $\rm s^{-1}$}}.

\subsubsection{MYtorus model} \label{sec:MY}

A second physically motivated model for highly obscured AGN is the \textsc{MYtorus} model developed by \cite{Murphy2009}, which includes three components: the intrinsic absorbed power law continuum (\textsc{MYtorusZ}), the scattered continuum (\textsc{MYtorusS}) and the line emission (\textsc{MYtorusL}). In order to avoid artefacts present in the \textsc{MYtorus} model at low energies we cut the spectrum at 0.65 keV instead of 0.5 keV. As before, we add a scattered power law and an APEC component so that the total model is: \\

\textsc{phabs$_{\rm GAL}$} $\times$ (\textsc{MYtorusZ} $\times$ \textsc{zpow} $+$ \textsc{MYtorusS} $+$ \textsc{MYtorusL} $+$ \textsc{zpow} $+$ \textsc{apec}). \\

\textcolor{black}{\rm{As in the previous model, the photon index of the intrinsic power law is fixed at $\Gamma = 1.9$ while the photon index of the unabsorbed power law is let free to vary. The power law parameters of the MYTorus components are tied together.}} Due to the low signal to noise of the spectrum, the column density and normalizations of the different \textsc{MYtorus} components are considered to be the same \textcolor{black}{\rm{(the available data do not allow us to constrain a different normalization of the scattered and line components with respect to the zeroth-order continuum represented by \textsc{MYtorusZ})}}. In addition, we fix the inclination angle to $\theta_{\rm i} = 90^{\rm o}$ \textcolor{black}{\rm{(Fig. \ref{fig:spec}, right, and Fig. \ref{fig:spec_2})}} and $\theta_{\rm i} = 75^{\rm o}$. The fits return \textcolor{black}{\rm{$N_{\rm H} \sim 2 \times 10^{24}$ cm$^{-2}$}} and \textcolor{black}{\rm{$L_{\rm 2-10 keV} \sim 1.2 \times 10^{43}$ erg $\rm s^{-1}$}} (see Table \ref{tab:res} for the values corresponding to the individual fits).

\begin{figure}
\includegraphics[scale=0.44]{spec_mod_2.pdf}
\caption{\textcolor{black}{\rm{{\it{Top}}: same as Fig. \ref{fig:spec}, right, but for the unfolded spectrum. {\it{Bottom}}: unfolded model components for the best-fitting \textsc{MYTorus} model with $\theta_{\rm i} = 90^{\rm o}$.}}}
\label{fig:spec_2}
\end{figure}

\begin{table*}
 \begin{tabular}{lcc}
  \hline
  \hline
  \noalign{\smallskip}
  \multicolumn{3}{c}{Phenomenological model (Sec. \ref{sec:PE})}\\
  \noalign{\smallskip}
  \hline
  \noalign{\smallskip}
  Column density $N_{\rm H}$ ($10^{24}$ $\rm cm^{-2}$) & \multicolumn{2}{c}{$1.87^{+1.27}_{-0.21}$}\\
  \noalign{\smallskip}
  Reflection parameter $R$ & \multicolumn{2}{c}{$\le 0.02$}\\
  \noalign{\smallskip}
  Plasma temperature $kT$ (keV) & \multicolumn{2}{c}{$0.95^{+0.1}_{-0.08}$}\\
  \noalign{\smallskip}
  Soft $\Gamma$ & \multicolumn{2}{c}{$2.32^{+0.5}_{-0.14}$}\\
  \noalign{\smallskip}
  C-stat/dof & \multicolumn{2}{c}{$246.8/222$}\\
  \noalign{\smallskip}
  Intrinsic 2-10 keV luminosity $L_{\rm 2-10 keV}$ (erg $\rm s^{-1}$) & \multicolumn{2}{c}{$(2.3^{+4.4}_{-0.8}) \times 10^{43}$}\\
  \noalign{\smallskip}
  Bolometric luminosity $L_{\rm bol}$ (erg $\rm s^{-1}$) & \multicolumn{2}{c}{$(4.9^{+14.9}_{-2.1}) \times 10^{44}$}\\
  \noalign{\smallskip}
  \hline
  \hline
  \noalign{\smallskip}
  \multicolumn{3}{c}{BNtorus model (Sec. \ref{sec:BN})}\\
  \noalign{\smallskip}
  \hline
  \noalign{\smallskip}
  Column density $N_{\rm H}$ ($10^{24}$ $\rm cm^{-2}$) & \multicolumn{2}{c}{$2.10^{+1.05}_{-0.57}$}\\
  \noalign{\smallskip}
  Inclination angle $\theta_{\rm i}$$^{\rm (1)}$ & \multicolumn{2}{c}{$87.1^{\rm o}$}\\
  \noalign{\smallskip}
  Torus opening angle $\theta_{\rm {oa}}$$^{\rm (1)}$ & \multicolumn{2}{c}{$60^{\rm o}$}\\
  \noalign{\smallskip}
  Plasma temperature $kT$ (keV) & \multicolumn{2}{c}{$1.00^{+0.15}_{-0.1}$}\\
  \noalign{\smallskip}
  Soft $\Gamma$ & \multicolumn{2}{c}{$2.86^{+0.12}_{-0.33}$}\\
  \noalign{\smallskip}
  C-stat/dof & \multicolumn{2}{c}{$245.3/217$}\\
  \noalign{\smallskip}
  Intrinsic 2-10 keV luminosity $L_{\rm 2-10 keV}$ (erg $\rm s^{-1}$) & \multicolumn{2}{c}{$(1.0^{+0.4}_{-0.4}) \times 10^{43}$}\\
  \noalign{\smallskip}
  Bolometric luminosity $L_{\rm bol}$ (erg $\rm s^{-1}$) & \multicolumn{2}{c}{$(1.9^{+1.0}_{-0.8}) \times 10^{44}$}\\
  \noalign{\smallskip}
  \hline
  \hline
  \noalign{\smallskip}
  \multicolumn{3}{c}{MYtorus model (Sec. \ref{sec:MY})}\\
  \noalign{\smallskip}
  \hline
  \noalign{\smallskip}
   & $\theta_{\rm i} = 90^{\rm o}$ & $\theta_{\rm i} = 75^{\rm o}$ \\
  \noalign{\smallskip}
  Column density $N_{\rm H}$ ($10^{24}$ $\rm cm^{-2}$) & $1.81^{+0.70}_{-0.36}$ & $2.09^{+1.03}_{-0.49}$ \\
  \noalign{\smallskip}
  Plasma temperature $kT$ (keV) & $0.98^{+0.23}_{-0.10}$ & $1.01^{+0.15}_{-0.14}$ \\
  \noalign{\smallskip}
  Soft $\Gamma$ & $2.62^{+0.34}_{-0.32}$ & $2.90^{+0.09}_{-0.33}$ \\
  \noalign{\smallskip}
  C-stat/dof & $239.8/213$ & $241.9/213$ \\
  \noalign{\smallskip}
  Intrinsic 2-10 keV luminosity $L_{\rm 2-10 keV}$ (erg $\rm s^{-1}$) & $(1.3^{+1.1}_{-0.5}) \times 10^{43}$ & $(1.2^{+0.6}_{-0.4}) \times 10^{43}$ \\
  \noalign{\smallskip}
  Bolometric luminosity $L_{\rm bol}$ (erg $\rm s^{-1}$) & $(2.4^{+2.8}_{-1.0}) \times 10^{44}$ & $(2.1^{+1.4}_{-0.8}) \times 10^{44}$ \\
  \noalign{\smallskip}
  \hline
 \end{tabular}
 \caption{Summary of the joint {\it{Chandra}} \, and {\it{NuSTAR}} \, X-ray spectral analysis. Details about the different models are given in the text. $^{\rm (1)}$ Frozen parameters. \textcolor{black}{\rm{The errors on the fit parameters correspond to the $90\%$ confidence range. The intrinsic 2-10 keV luminosities are intended as the luminosity of the intrinsic power law.}}}
 \label{tab:res}
\end{table*}

\subsection{Luminosity estimation}\label{sec:lum}

\textcolor{black}{\rm{ In this section we estimate the current bolometric luminosity from the X-ray fits and from mid-IR data, and the past bolometric luminosity ($\sim 100$ kyr ago) using optical spectroscopy. The results are summarized in Fig. \ref{fig:lum}.}}

\begin{figure*}
\includegraphics[scale=0.6]{lum.png}
\caption{\textcolor{black}{\rm{Past and current bolometric luminosity ranges as estimated from optical (green), X-ray (purple) and mid-IR (blue) data. The past luminosity range corresponds to the range between the minimum and maximum values obtained from optical data. The current bolometric luminosity range from X-ray data corresponds to the range obtained using the $68\%$ (strong purple) and $90\%$ (light purple) confidence intervals on the intrinsic power law normalization. The mid-IR data allows to compute an upper limit for the current bolometric luminosity. See Section \ref{sec:lum} for more details.}}}
\label{fig:lum}
\end{figure*}

\subsubsection{Current bolometric luminosity from the X-ray spectrum}

The simultaneous {\it{Chandra}} \, and {\it{NuSTAR}} \, analysis of the  $0.5-24.0$\,keV spectrum of IC\,2497 suggests that the source is highly obscured, \textcolor{black}{\rm{$N_{\rm H} = 1.8-2.1 \times 10^{24}$ cm$^{-2}$}}, with intrinsic 2-10 keV luminosity \textcolor{black}{\rm{ (luminosity of the intrinsic power law) $L_{\rm 2-10 keV} = 1 - 2.3 \times 10^{43}$ erg s$^{-1}$}}.

\textcolor{black}{\rm{Assuming the luminosity dependent X-ray to bolometric correction from \cite{Marconi2004}, the measured 2-10 keV luminosity corresponds to a bolometric luminosity $L_{\rm bol} = 1.9-4.9 \times 10^{44}$ erg $\rm s^{-1}$.}}

\textcolor{black}{\rm{Fig. \ref{fig:lum} shows the range covered by the estimated bolometric luminosity considering also the uncertainties on the intrinsic power law. For every model we compute the 2-10 keV luminosity for the best-fitting normalization and for the minimum and maximum value allowed from the $90\%$ and $68\%$ confidence regions. We then infer the bolometric luminosity as explained above. In the plot we show the mean value obtained from the best-fit parameters, and the minimum and maximum values for the two confidence levels measured with the different models.}}

\subsubsection{Current bolometric luminosity from mid-IR emission}\label{sec:mir}

The {\it{WISE}} observations \citep{Wright2010} of IC\,2497 allow us to compute an upper limit to the current bolometric luminosity of the obscured quasar. First, the {\it{W1}}-{\it{W2}} colour is low, {\it{W1}}-{\it{W2}} = 0.3\footnote{We checked that no variability is present between {\it{WISE}} and {\it{NEOWISE}} (\citealt{Mainzer2011}).}, which suggests that the mid-IR flux is dominated by stellar emission (e.g. \citealt{Stern2005}; \citealt{Jarrett2011}; \citealt{Stern2012}). Moreover, the emission detected with {\it{W3}} and {\it{W4}} (12 and 22 $\mu$m, respectively) appears effectively not to be absorbed e.g. by dust also in the case of highly obscured quasars (e.g. \citealt{Gandhi2009}; \citealt{Asmus2015}), so that the observed luminosity corresponds to the intrinsic one. Assuming a radio quiet quasar SED template from \cite{Elvis1994} and normalising it to the lowest non absorbed {\it{WISE}} point, i.e. {\it{W3}}, we obtain an upper limit for the bolometric luminosity $L_{\rm bol} = 5.0 \times 10^{44}$ erg $\rm s^{-1}$. This value is consistent with that obtained from the X-ray analysis.

\subsubsection{Past bolometric luminosity from optical emission}

A detailed analysis of the optical longslit spectra, which allow the computation of the lower and upper limits to the AGN output needed to photoionize HV, is presented in \cite{Lintott2009}. Here we provide a brief summary of their findings. A lower limit for the required AGN output is obtained from simple recombination balance, i.e. by assuming that the number of ionizations is equal to the number of recombinations, and that the number of ionizing photons is high enough to power the observed recombination lines. The integrated $\rm H \beta$ luminosity of HV is $1.4 \times 10^{41}$ erg s$^{-1}$. Assuming typical nebular conditions and a flat ionizing continuum ($F_{\nu} \propto \nu^{-1}$) this corresponds to a lower limit for the ionizing luminosity $L_{\rm ion} > 1 \times 10^{45}$\,erg\,s$^{-1}$. On the other hand, the ionization parameter obtained from the [$\mathrm{O}$\textsc{ii}]$\lambda 3727$/[$\mathrm{O}$\textsc{iii}]$\lambda 5007$ ratio of HV is log(U)\,=\,$-2.2$ while the [$\mathrm{S}$\textsc{ii}]$\lambda 6717 /6731$ doublet ratio provides an upper limit for the electron density $n_{\rm e} < 50$\,cm$^{-3}$. This provides an upper limit for the ionizing photon density $n_{\rm phot} < 0.32$ cm$^{-3}$. Assuming a flat ionizing continuum and the mean projected separation between IC\,2497 and HV, $d=20$ kpc, the photon density obtained corresponds to an upper limit to the ionizing luminosity $L_{\rm ion} < 3.2 \times 10^{45}$\,erg\,s$^{-1}$.

\textcolor{black}{\rm{A similar analysis was performed also in \cite{Keel2017}. Based on {\it{HST}} narrow-band H$\alpha$ images and photoionization balance, they infer a lower limit (because of possible dust obscuration and/or a clumpy gas structure) for the isotropic emission rate of ionizing photons $Q_{\rm ion} \sim 3 \times 10^{55}$ photons s$^{-1}$ for almost the whole time range probed by the photoionized cloud, which again corresponds to a lower limit for the ionizing luminosity $L_{\rm ion} > 1.1 \times 10^{45}$\,erg\,s$^{-1}$.}}

\textcolor{black}{\rm{Altogether, the observations point to a past ionizing luminosity in the range $L_{\rm ion} = 1 - 3.2 \times 10^{45}$\,erg\,s$^{-1}$. By scaling a radio quasar template from \citealt{Elvis1994} to match the estimated ionizing luminosity between 1 and 4 Ryd, we infer a past bolometric luminosity $L_{\rm bol} = 0.6 - 2.1 \times 10^{46}$\,erg\,s$^{-1}$.}}

\section{Discussion}
We used new hard X-ray {\it{NuSTAR}} \, observations of IC\,2497, together with archival {\it{Chandra}} \, observations and {\it{WISE}} fluxes, to investigate the current state of its quasar. The data rule out the possibility that the quasar is currently radiating \textcolor{black}{at a rate} $L_{\rm bol} \sim 10^{46}$ erg s$^{-1}$ as required by the photoionization state of Hanny's Voorwerp, and suggest a drop in luminosity by a factor of $\sim$\textcolor{black}{\rm{50}} in the last $\sim$ 100 kyr. In addition, the newly obtained observations suggest that the source is Compton thick with $N_{\rm H} \sim 2 \times$\,{\textcolor{black}{$10^{24}$}}\,cm$^{-2}$. However, the obscuration is not the explanation for the discrepancy between the needed and observed levels of AGN emission; this suggests that significant AGN variability is required to explain the properties of this system. These findings also support the idea first proposed by \cite{Lintott2009} that Hanny's Voorwerp is a quasar ionization echo, although the magnitude of the drop in luminosity is significantly lower than previously thought (factor $\sim$ \textcolor{black}{\rm{50}} instead of $> 100$, \textcolor{black}{\rm{assuming the mean current and past luminosities obtained with our analysis}}).

Assuming a stellar mass = $M_* \sim 1.7 \times 10^{11} M_\odot$ as estimated from the K band magnitude and the $M_* - M_{BH}$ relation in \cite{Jahnke2009}, \textcolor{black}{\rm{as well as the mean current and past luminosities obtained with our analysis}}, the measured drop in bolometric luminosity corresponds to a drop in Eddington ratio (ER) from $\rm \lambda_{Edd} \sim$ \textcolor{black}{\rm{0.35}} to $\rm \lambda_{Edd} \sim$ \textcolor{black}{\rm{0.007}}. The ER found for luminous AGN is typically $\rm \lambda_{Edd} \sim 0.01 - 1$ (e.g. \citealt{Panessa2006}; \citealt{Kollmeier2006}). An ER $\rm \lambda_{Edd} \sim 0.007$ is more consistent with the predictions for a radiative inefficient (RIAF) or advection-dominated accretion flow (ADAF; e.g. \citealt{Narayan1995}; \citealt{Narayan1998}; \citealt{Ho2009}), where most of the energy is radiated in kinetic form. The drop in ER could therefore suggest that IC\,2497 is entering a regime where we might expect a change from a radiatively efficient to a radiatively inefficient accretion state (e.g. \citealt{Done2007}; \citealt{Alexander2012}). This hypothesis is supported also by the fact that IC\,2497 and its quasar show other characteristics common to ADAF sources, such as the presence of a compact radio source or jet, and the prevalence of low-ionization nebular conditions (e.g. \citealt{Terashima2002}; \citealt{Ptak2004}; \citealt{Ho2009}).

Accretion state changes are common in X-ray binaries. In these systems the switch from a radiatively efficient to a radiatively inefficient state is often linked to changes in the jet activity, and provides a source of mechanical feedback (e.g. \citealt{Done2007}; \citealt{Pakull2010}). As observed by \cite{Jozsa2009} and \cite{Rampadarath2010}, IC\,2497 shows a large scale radio jet which connects the galaxy to Hanny's Voorwerp, and a jet hotspot $\sim 300$ pc away from the AGN which may \textcolor{black}{\rm{correspond}} to a jet launched after the quasar started dropping in luminosity. A hypothesis is that the nuclear jet is connected to the current accretion state change, in analogy to what is observed in X-ray binaries, while the large scale radio jet was launched during a previous AGN phase. In fact, as suggested by \cite{Schawinski2015}, AGN may ``flicker" on and off 100-1000 times, with a typical AGN phase lasting $\sim 10^5$ yr, so that the host galaxy may present signatures related to different AGN phases. The idea of the AGN currently doing mechanical feedback on the galaxy is supported also by the presence of a nuclear outflow and a nuclear ring of expanding gas observed in the optical {\it{HST}} images (\citealt{Keel2012b}) and of a bubble expanding in the hot gas around the AGN seen in the {\it{Chandra}} \, data (\citealt{Sartori2016}). The analogy between IC\,2497 and what is observed in X-ray binaries provides further indication of a possible unification of black hole accretion physics from stellar mass BHs to SMBHs. However, it is important to notice that different time-scales have to be taken into account since the dynamical and viscous time-scales increase with black hole mass. If we assume the underlying physics for black hole accretion to be the same, the analogous of a state transition lasting $\sim 1$h for an $\sim 10M_{\odot}$ BH would last $\sim 10^{4-5}$yr for a $\sim 10^{9} M_{\odot}$ SMBH (see \citealt{Schawinski2010b}, \citealt{Sartori2016} and references therein for a more detailed description).

Other examples of AGN variability are given by the so-called {\it{changing-look AGN}}, which present a change in the AGN type (e.g. from Type 1 to Type 1.9) due to broadening or narrowing of the Balmer lines{\footnote{In this \textcolor{black}{\rm{section}} we refer to optical changing-look AGN, while other studies describe systems where the X-ray changing-look behaviour, from Compton thick to Compton thin or viceversa, is due to a change in the obscuration (e.g. \cite{Risaliti2005}; \citealt{Ricci2016b})} (\citealt{Denney2014}; \citealt{LaMassa2015}; \citealt{Ruan2015}; \citealt{MacLeold2016}; \citealt{Runnoe2016}; \citealt{Husemann2016}; \citealt{McElroy2016}; \citealt{Gezari2016}; \textcolor{blue}{Stern et al. submitted}). The appearance or disappearance of broad emission lines is often accompanied by a change in luminosity of a factor $\sim$10 over $\sim$10 yr time-scales. As described above, these time-scales are much shorter compared to time-scales expected for accretion state changes (e.g. \citealt{Sobolewska2011}; \citealt{Hickox2014}), and possible alternative explanations of the changing-look behaviour are variable absorption due to a clumpy torus (e.g. \citealt{Elitzur2012}), transient events, e.g. tidal disruption of a star by the central black hole (e.g. \citealt{Eracloeous1995}; \citealt{Merloni2015}), or major changes in the photoionization balance. The magnitude of the drop in luminosity measured in IC\,2497 is \textcolor{black}{\rm{at least}} a factor $\sim 2$ higher than what has been observed in changing-look AGN. Moreover, the {\it{Chandra}} \, and {\it{NuSTAR}} \, data do not show significant variability, and the upper limits obtained from archival {\it{WISE}}, {\it{NEOWISE}} and {\it{IRAS}} data seem to exclude that the total drop in luminosity happened within the last decades. For these reasons, we argue that the AGN in IC\,2497 should not be classified as a changing-look AGN. On the other hand, we suggest that changing-look AGN correspond to a short-time ($\sim 10-100$ yr) variability which is superimposed on the long-term ($\sim 10^{5-6}$ yr) AGN phases suggested by this work and other observations (e.g. \citealt{Schawinski2015}), high resolution (sub-kpc) simulations (\citealt{Hopkins2010}; \citealt{Novak2011}; \citealt{Bournaud2011}; \citealt{Gabor2013}; \citealt{DeGraf2014}; \citealt{Sijacki2015}) and theoretical models (\citealt{Siemiginowska1997}; \citealt{Sanders1981}; \citealt{DiMatteo2005}; \citealt{Hopkins2005}; \citealt{Springel2005}; \citealt{King2007}; \citealt{King2015}).

\section{Conclusions}
In this paper we present new {\it{NuSTAR}} \, observations of IC\,2497, the galaxy associated with Hanny's Voorwerp. The analysis of these new data, together with archival {\it{Chandra}} \, X-ray observations, optical longslit spectroscopy \textcolor{black}{{\rm{and narrow-band imaging}}}, and {\it{WISE}} \, mid-IR photometry, shows that the AGN in IC\,2497 is Compton thick, and its luminosity dropped by a factor of $\sim$\textcolor{black}{\rm{50}} within the last $\sim 100$ kyr. These findings support the idea that Hanny's Voorwerp is an ionization echo, although the magnitude of the quasar's drop in luminosity is significantly lower than previously thought. On the other hand, the observed change in Eddington ratio, from $\rm \lambda_{Edd} \sim$ \textcolor{black}{\rm{0.35}} to $\rm \lambda_{Edd} \sim$ \textcolor{black}{\rm{0.007}}, suggests that the quasar in IC\,2497 is entering a regime where it is switching from a radiatively efficient to a radiatively inefficient accretion state, in a similar way as observed in X-ray binaries.

\section*{Acknowledgements}

\textcolor{black}{\rm{We thank the anonymous referee for helpful comments which improved the quality of the manuscript. We also thank the {\it{NuSTAR}} \, Cycle 2 TAC for the {\it{NuSTAR}} \, data on which this paper is based.}}
LFS and KS gratefully acknowledge support from Swiss National Science Foundation Grants PP00P2\_138979 and PP00P2\_166159. MJK acknowledges support from the Swiss National Science Foundation through the Ambizione fellowship grant PZ00P2\_154799/1. MJK acknowledges support from NASA through ADAP award NNH16CT03C. ET acknowledge support from: CONICYT-Chile grants Basal-CATA PFB-06/2007 and FONDECYT Regular 1160999. CR acknowledges the China-CONICYT fund, FONDECYT 1141218 and Basal-CATA PFB--06/2007. WPM acknowledges partial support by HST guest observer programme \#14271  under NASA contract NAS 5-26555 (HST) and Chandra Award GO5-16101X under NASA contract NAS8-03060 (CXC). GL acknowledges support from a Herchel Smith Research Fellowship of the University of Cambridge.  This work is based on observations with the {\it{Chandra}} \, satellite and has made use of software provided by the {\it{Chandra}} \, X-ray Center (CXC) in the application package \textsc{CIAO}. This work made use of data from the {\it{NuSTAR}} \, mission, a project led by the California Institute of Technology, managed by the Jet Propulsion Laboratory, and funded by the National Aeronautics and Space Administration. We thank the {\it{NuSTAR}} \, Operations, Software and Calibration teams for support with the execution and analysis of these observations. This research has made use of the {\it{NuSTAR}} \, Data Analysis Software (\textsc{NuSTARDAS}) jointly developed by the ASI Science DataCenter (ASDC, Italy) and the California Institute of Technology (USA). The research has made use of NASA's Astrophysics Data System Bibliographic Service.


\section{Introduction}

The Lambek calculus~\cite{Lambek58} was introduced for describing natural language syntax
by means of type-logical (categorial) grammars. Further research on type-logical grammar
showed that the original system proposed by Lambek appears to be insufficient to cover
intrinsic linguistic phenomena, and so various extensions and modifications of the Lambek
calculus were introduced (see, for example, books by Morrill~\cite{MorrillBook},
Moot and Retor\'{e} \cite{MootRetore} and others). One of
these extensions is the {\em Lambek calculus with brackets} introduced by Moortgat~\cite{Moortgat}.
While the original Lambek calculus is fully associative, brackets block associativity in
specified situations, thus disallowing derivations of ungrammatical phrases like
``the girl whom John loves Mary and Pete loves'' (see~\cite{MorrillBook} for a more detailed
analysis).

The original Lambek grammars generate precisely context-free languages, as shown by Pentus~\cite{PentusCFG}.
For the Lambek calculus extended with brackets, the context-free upper bound was claimed by J\"ager~\cite{Jaeger}.
Unfortunately, J\"ager's argument relied upon a lemma by Versmissen~\cite{Versmissen}, which was
afterwards shown to be incorrect~\cite{FaddaMorrill}\cite{FSCD}.

Recently, however, Kanazawa~\cite{KanazawaArXiv}
returned to this question and presented a new proof, based on a insightful combination of J\"ager's ideas and
the original Pentus' approach. Kanazawa proved this result for two versions of the Lambek calculus with brackets:
not allowing empty antecedents (we denote it by $\mathbf{Lb}$) and allowing them ($\mathbf{Lb}^{\boldsymbol{*}}$).

The third variant of the calculus, $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$, is obtained from $\mathbf{Lb}^{\boldsymbol{*}}$ by adding the {\em multiplicative unit constant,}
$\mathbf{1}$. As noticed by Kanazawa~\cite{KanazawaArXiv}, Pentus-style reasoning is not applicable to the case with
the unit constant (even without brackets). In Pentus' proof, the target context-free grammar rules are essentially
all derivable sequents of a bounded size. For $\mathbf{Lb}$ and $\mathbf{Lb}^{\boldsymbol{*}}$, the set of these sequents is finite; however,
the measure of size used by Pentus ignores occurrences of $\mathbf{1}$, thus this set of sequents becomes infinite, and
doesn't yield a context-free grammar. A workaround for this issue (in the case without brackets) was presented
by Kuznetsov~\cite{FG2011}: the unit constant gets eliminated by a faithful translation of the calculus with
the unit into the original system. The same problem for the case with brackets ({\it i.e.,} whether $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$-grammars
generate precisely context-free languages) was left as an open question in~\cite{KanazawaArXiv}.

In this paper, we extend the construction from~\cite{FG2011} to an embedding of $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$ into $\mathbf{Lb}^{\boldsymbol{*}}$.
Thus, we show that $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$-grammars generate the same class of languages as $\mathbf{Lb}^{\boldsymbol{*}}$-grammars---and, due
to Kanazawa~\cite{KanazawaArXiv}, it is the class of context-free languages.

The translation of $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$ to $\mathbf{Lb}^{\boldsymbol{*}}$ also has an algorithmic application. While the original Lambek calculus,
and, therefore, the Lambek calculus with brackets is NP-complete (Pentus~\cite{PentusNP}), for sequents of
bounded formula complexity and bracket nesting depth Kanovich {\it et al.}~\cite{FSCD}, extending Pentus~\cite{PentusPoly},
present a polynomial-time algorithm for deciding derivability in $\mathbf{Lb}^{\boldsymbol{*}}$. Our translation generalises this algorithm
from $\mathbf{Lb}^{\boldsymbol{*}}$ to $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$.

\section{The Calculi $\mathbf{Lb}^{\boldsymbol{*}}$ and $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$}

In this section and further we use the notation of~\cite{FSCD}; the notation used
in~\cite{KanazawaArXiv} is slightly different.

The syntax of the Lambek calculus with brackets is a bit more involved than a standard
sequential calculus. Formulae of $\mathbf{Lb}^{\boldsymbol{*}}$ are recursively built from variables
($\mathrm{Var} = \{ p_1, p_2, \dots \}$) using three binary connectives,
$\mathop{\backslash}$ (left division), $\mathop{/}$ (right division), $\cdot$ (product) and two unary ones,
$\langle\rangle$ and $[]^{-1}$. The binary connectives come from the original Lambek calculus~\cite{Lambek58};
as one can easily see, the system $\mathbf{Lb}^{\boldsymbol{*}}$ defined below is a conservative extension of the
version of the Lambek calculus which allows empty antecedents~\cite{Lambek61}. The unary
connectives operate brackets that are used to introduce controlled non-associativity.
Sequents of $\mathbf{Lb}^{\boldsymbol{*}}$ are expressions of the form $\Pi \to C$, where $C$ is a formula and
$\Pi$ is a {\em meta-formula.} Meta-formulae are built from formulae using comma (,) and
brackets ([...]). By definition, comma is associative, and the empty meta-formula $\Lambda$
is the unit object w.r.t.\ comma: $\Gamma,\Lambda$ and $\Lambda,\Gamma$ are both considered
graphically equal to $\Gamma$. By $\Delta(\Pi)$ we denote a meta-formula $\Delta$ with
a {\em designated occurrence} of a sub-meta-formula $\Pi$.

Axioms and rules of $\mathbf{Lb}^{\boldsymbol{*}}$ are as follows.

$$
\infer[(\mathrm{ax})]{p_i \to p_i}{}
$$
$$
\infer[(\mathop{\backslash}\to)]{\Delta(\Pi, A \mathop{\backslash} B) \to C}{\Pi \to A & \Delta(B) \to C}
\qquad
\infer[(\to\mathop{\backslash})]{\Pi \to A \mathop{\backslash} B}{A, \Pi \to B}
$$
$$
\infer[(\mathop{/}\to)]{\Delta(B \mathop{/} A, \Pi) \to C}{\Pi \to A & \Delta(B) \to C}
\qquad
\infer[(\to\mathop{/})]{\Pi \to B \mathop{/} A}{\Pi, A \to B}
$$
$$
\infer[(\cdot\to)]{\Delta(A \cdot B) \to C}{\Delta(A, B) \to C}
\qquad
\infer[(\to\cdot)]{\Pi, \Gamma \to A \cdot B}{\Pi \to A & \Gamma \to B}
$$
$$
\infer[(\langle\rangle\to)]{\Delta(\langle\rangle A) \to C}{\Delta([A]) \to C}
\qquad
\infer[(\to\langle\rangle)]{[\Pi] \to \langle\rangle A}{\Pi \to A}
$$
$$
\infer[([]^{-1}\to)]{\Delta([[]^{-1} A]) \to C}{\Delta(A) \to C}
\qquad
\infer[(\to[]^{-1})]{\Pi \to []^{-1} A}{[\Pi] \to A}
$$

Adding the unit constant, $\mathbf{1}$, with the following rules of inference (cf.~\cite{Lambek69}),
$$
\infer[(\mathbf{1}\to)]{\Delta(\mathbf{1}) \to C}{\Delta(\Lambda) \to C}
\qquad
\infer[(\to\mathbf{1})]{\to\mathbf{1}}{}
$$
yields the calculus $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$.

The cut rule of the following form is admissible~\cite{Moortgat}.
$$
\infer[(\mathrm{cut})]{\Delta(\Pi) \to C}{\Pi \to A & \Delta(A) \to C}
$$
Admissibility of cut allows replacing a subformula with an equivalent one preserving
derivability.

\section{Translating $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$ to $\mathbf{Lb}^{\boldsymbol{*}}$}

In this section we present a translation of $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$ formulae that eliminates the unit constant.

Informally, we replace each occurrence of $\mathbf{1}$ by $(q \mathop{\backslash} q)$, where
$q$ is a fresh variable. For classical propositional logic, this would be sufficient, since
$(\varphi \Rightarrow \varphi)$ there is equivalent to the ``{\sc true}'' constant. In order to make
this construction work for the substructural
system $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$, however, we also need to add some extra $(q \mathop{\backslash} q)$'s,
depending on the polarity of specific subformula.

Formally, we define two translations, $\tau^+$ and $\tau^-$, by joint induction:
\begin{align*}
& \tau^+(\mathbf{1}) = \tau^-(\mathbf{1}) = q \mathop{\backslash} q\\
& \tau^+(p_i) = (q \mathop{\backslash} q) \cdot p_i \cdot (q \mathop{\backslash} q) && \tau^-(p_i) = p_i\\
& \tau^+(A \mathop{\backslash} B) = \tau^-(A) \mathop{\backslash} \tau^+(B) && \tau^-(A \mathop{\backslash} B) = \tau^-(A) \mathop{\backslash} \tau^-(B) \\
& \tau^+(B \mathop{/} A) = \tau^+(B) \mathop{/} \tau^+(A) && \tau^-(B \mathop{/} A) = \tau^-(B) \mathop{/} \tau^+(A)\\
& \tau^+(A \cdot B) = \tau^+(A) \cdot \tau^+(B) && \tau^-(A \cdot B) = \tau^-(A) \cdot \tau^-(B)\\
& \tau^+(\langle\rangle A) = (q \mathop{\backslash} q) \cdot \langle\rangle \tau^+(A) \cdot (q \mathop{\backslash} q) &&
\tau^-(\langle\rangle A) = \langle\rangle \tau^-(A)\\
& \tau^+([]^{-1} A) = []^{-1} \tau^-(A) &&
\tau^-([]^{-1} A) = (q \mathop{\backslash} q) \mathop{\backslash} []^{-1} \tau^-(A) \mathop{/} (q \mathop{\backslash} q)
\end{align*}

For metaformulae, we define only $\tau^-$:
\begin{align*}
&\tau^-(\Lambda) = \Lambda\\
&\tau^-(\Gamma, \Delta) = \tau^-(\Gamma), \tau^-(\Delta)\\
&\tau^-([\Gamma]) = [\tau^-(\Gamma)]
\end{align*}

\begin{theorem}\label{Th:main}
For any $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$-sequent $\Pi \to C$ that has no occurrences of $q$,
$\Pi \to C$ is derivable in $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$ if{f} $\tau^-(\Pi) \to \tau^+(B)$ is
derivable in $\mathbf{Lb}^{\boldsymbol{*}}$.
\end{theorem}

In order to make the proof more convenient, we first reformulate $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$: remove $(\mathbf{1}\to)$ and replace some of the other rules with the following ones
(here and further $\mathbf{1}^n$ means $\underbrace{\mathbf{1}, \mathbf{1}, \dots, \mathbf{1}}_{\text{$n$ times}}$).
$$
\infer[(\to\mathbf{1})']
{\mathbf{1}^k \to \mathbf{1}}{}
\qquad
\infer[(\mathrm{ax})']
{\mathbf{1}^k, p_i, \mathbf{1}^m \to p_i}{}
$$
$$
\infer[(\to\langle\rangle)']
{\mathbf{1}^k, [\Pi], \mathbf{1}^m \to \langle\rangle A}
{\Pi \to \langle\rangle A}
\qquad
\infer[([]^{-1}\to)']
{\Delta([\mathbf{1}^k, []^{-1} B, \mathbf{1}^m]) \to C}
{\Delta(B) \to C}
$$
All other rules are left intact.

Let's call the new calculus ${\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}}'$.

\begin{lemma}\label{Lm:LbUp}
A sequent is derivable in ${\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}}'$ if{f} it is derivable in
$\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$.
\end{lemma}

\begin{proof}
Since the new $'$-rules become the original rules of $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$ when $k=m=0$, for the
``if'' part it is sufficient to show that $(\mathbf{1}\to)$ is admissible in ${\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}}'$, {\em i.e.,}
that if $\Psi(\Lambda) \to C$ is derivable in ${\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}}'$, then so is $\Psi(\mathbf{1}) \to C$.

Proceed by induction on derivation.
Let's denote $\Pi$, $\Gamma$, and $\Delta$ (except for the designated part that is affected by the rule)
in the rules of $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$ as the {\em context}. Consider the last rule that derives $\Psi(\Lambda) \to C$.

If this designated occurrence of $\Lambda$ is inside the context of the last rule that derives $\Psi(\Lambda) \to C$,
then we can trace it upwards to the premise (one of the premises) of the rule, replace $\Lambda$ by $\mathbf{1}$ there
(the sequent is still derivable by induction hypothesis), and then apply the rule. In other words, in this case
the two rules are exchangeable, and we can propagate the $(\mathbf{1}\to)$ rule upwards. For $(\mathop{/}\to)$, $(\mathop{\backslash}\to)$,
$(\to\mathop{/})$, $(\to\mathop{\backslash})$, $(\cdot\to)$, $(\to\cdot)$, $(\langle\rangle\to)$, and $(\to[]^{-1})$ this is the only possible
situation: a possible place for inserting $\mathbf{1}$ is always in the context.

If $\Psi(\Lambda) \to C$ is an axiom of the form $(\mathrm{ax})'$, namely, $\mathbf{1}^k, p_i, \mathbf{1}^m \to p_i$, then $\Psi(\mathbf{1}) \to C$ is either
$\mathbf{1}^{k+1}, p_i, \mathbf{1}^m \to p_i$ or $\mathbf{1}^k, p_i, \mathbf{1}^{m+1} \to p_i$; both are again instances of $(\mathrm{ax})'$. The $(\to\mathbf{1})'$ case
is handled in the same way. For the $(\to\langle\rangle)'$ case, the only interesting situation is when $\mathbf{1}$ is added outside
the brackets (the case when it is added into the context $\Pi$, was already considered); in this situation the extra $\mathbf{1}$ again
gets absorbed by the $(\to\langle\rangle)'$ rule. The $([]^{-1}\to)'$ case is dual.

For the ``only if'' part, just notice that each $'$-rule can be represented as a
consequent application of the corresponding rule of $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$ and then $k+m$ times
$(\mathbf{1}\to)$  (or $k$ times for $(\to\mathbf{1})'$).
\end{proof}

\begin{proof}[Proof of Theorem~\ref{Th:main}]
The ``only if'' part: let $\Pi \to C$ be derivable in $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$. Then by Lemma~\ref{Lm:LbUp}
it is derivable in ${\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}}'$. We need to show that $\tau^-(\Pi) \to \tau^+(C)$ is derivable
in $\mathbf{Lb}^{\boldsymbol{*}}$. Proceed by induction on the derivation of $\Pi\to C$ in ${\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}}'$. Applications of rules
without $'$ are translated into $\mathbf{Lb}^{\boldsymbol{*}}$ ``as is.'' The $(\to\mathbf{1})'$ transforms into $(q\mathop{\backslash} q)^k \to q \mathop{\backslash} q$,
which is derivable in $\mathbf{Lb}^{\boldsymbol{*}}$; $(\mathrm{ax})'$ becomes $(q \mathop{\backslash} q)^k, p_i, (q \mathop{\backslash} q)^m \to (q \mathop{\backslash} q) \cdot p_i \cdot (q \mathop{\backslash} q)$,
which is also derivable. Finally, $(\to\langle\rangle)'$ and $([]^{-1}\to)'$ are translated as follows:
$$
\infer[(\to\cdot) \text{ two times}]{(q \mathop{\backslash} q)^k, [\tau^-(\Pi)], (q \mathop{\backslash} q)^m \to (q \mathop{\backslash} q) \cdot \langle\rangle \tau^+(A) \cdot (q \mathop{\backslash} q)}
{(q \mathop{\backslash} q)^k \to q \mathop{\backslash} q & \infer[(\to\langle\rangle)]{[\tau^-(\Pi)] \to \langle\rangle \tau^+(A)}{\tau^-(\Pi) \to \tau^+(A)} & (q \mathop{\backslash} q)^m \to q \mathop{\backslash} q}
$$
$$
\infer[(\mathop{\backslash}\to) \text{ and } (\mathop{/}\to)]
{\Delta([(q \mathop{\backslash} q)^k, (q \mathop{\backslash} q) \mathop{\backslash} []^{-1} \tau^-(B) \mathop{/} (q \mathop{\backslash} q), (q \mathop{\backslash} q)^m]) \to C}
{(q \mathop{\backslash} q)^k \to q \mathop{\backslash} q & \infer[([]^{-1}\to)]{\Delta([[]^{-1} \tau^-(B)]) \to C}{\Delta(\tau^-(B)) \to C} &
(q \mathop{\backslash} q)^m \to q \mathop{\backslash} q}
$$

The ``if'' part is easier: if $\tau^-(\Pi) \to \tau^+(C)$ is derivable in $\mathbf{Lb}^{\boldsymbol{*}}$, then it is also derivable in $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$.
Substitute $\mathbf{1}$ for $q$ (substituting formulae for variables---but not for the unit constant!---preserves derivability).
Since $q$ is a fresh variable, the substitution affects only the $(q \mathop{\backslash} q)$ combinations introduced by $\tau^+$ and $\tau^-$,
and it is easy to see that for any formula $B$ after this substitution $\tau^+(B)$ and $\tau^-(B)$ become equivalent to $B$.
Thus, $\tau^-(\Pi) \to \tau^+(C)$ transforms to a sequent equivalent to $\Pi \to C$, therefore the latter one is derivable in $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$.
\end{proof}

\section{Applications}

The first application of Theorem~\ref{Th:main} is the characterisation of the class of languages generated by
$\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$-grammars.

Let $\Sigma$ be a finite alphabet. An $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$-grammar over $\Sigma$ is a triple $\langle \rhd, H \rangle$, where
$H$ is an $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$-formula (called the {\em target} formula), and $\rhd$ is a finite relation between letters $\Sigma$ and
$\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$-formulae. J\"ager~\cite{Jaeger} gives two definitions of a word $a_1\dots a_n$ being accepted by such a grammar.
The word is {\em s-accepted,} if there exist such formulae $A_1, \dots, A_n$ that $a_i \rhd A_i$ ($i = 1,\dots, n)$ and
the sequent $A_1, \dots, A_n \to H$ is derivable in $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$. The word is {\em t-accepted,} if, again, there exist $A_1, \dots, A_n$
such that $a_i \rhd A_i$, and there also exists a multiformula $\Pi$ such that $\Pi \to H$ is derivable and $A_1, \dots, A_n$
is obtained from $\Pi$ by erasing $[$ and $]$ ({\it i.e.,} it is the {\em yield} of the treelike bracketed structure $\Pi$).
The language t-generated (resp., s-generated) by the grammars defined as the set of all t-accepted (resp., s-accepted) words.

Theorem~\ref{Th:main} immediately yields the following corollary:
\begin{theorem}
The class of languages t-generated (resp., s-generated) by $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$-grammars, coincides with the class of languages
t-generated (resp., s-generated) by $\mathbf{Lb}^{\boldsymbol{*}}$-grammars.
\end{theorem}

\begin{proof}
For the inclusion in the non-trivial direction, apply $\tau^-$ to all formulae associated to letters of
$\Sigma$ by $\rhd$, and $\tau^+$ to the target formula $H$.
\end{proof}

In combination with Kanazawa's result that $\mathbf{Lb}^{\boldsymbol{*}}$-grammars t-generate exactly the class of context-free languages,
this gives a solution for the question left open in~\cite{KanazawaArXiv}: $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$-grammars also t-generate exactly all
context-free languages. We conjecture that our construction also works for the multimodal version of $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$,
solving the second open question from~\cite{KanazawaArXiv}.
The question about the class of s-generated languages appears to be still open, though much less interesting from
the linguistic point of view.

The second application is a fast algorithm for deciding derivability in $\mathbf{Lb}^{\boldsymbol{*}}$. While the full Lambek calculus is
NP-complete~\cite{PentusNP}, Pentus in~\cite{PentusPoly} introduces a depth parameter $d$ (in the product-free
case it is just the Horn depth of $\mathop{\backslash}$ and $\mathop{/}$) and presents an algorithm that works in polynomial time w.r.t.\ $n$
and $2^d$ (where $n$ is the length of the input sequent). Kanovich {\em et al.}~\cite{FSCD} generalise this
algorithm to $\mathbf{Lb}^{\boldsymbol{*}}$, adding the third parameter, $b$, which is the nesting depth of $[...]$, $\langle\rangle$, and $[]^{-1}$. The running
time of this algorithm is bounded by a polynom of $n$, $2^d$, and $n^b$ (if there are brackets, $n^b$ absorbs $n$). Since $\tau^+$ and
$\tau^-$ alter these parameters only linearly, Theorem~\ref{Th:main} yields an algorithm also for $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$ with the same running
time estimation:
\begin{theorem}
There exists an algorithm that decides derivability in $\mathbf{Lb}^{\boldsymbol{*}}_{\mathbf{1}}$ with running time bounded by a polynom of
$n$, $2^d$, and $n^b$.
\end{theorem}

\subsection*{Acknowledgments}

The author
was much inspired by recent and older works of Makoto Kanazawa (to which this note is a modest addition).
The author is also grateful to Lev Beklemishev, Max Kanovich, Glyn Morrill, Vivek Nigam, Mati Pentus, and Andre Scedrov for
in-depth discussions and exchange of ideas in a very warm, friendly, and collaborative atmosphere.


\section{Introduction}

Many features present in galaxy spectra, usually characterized by the intensity of absorption or emission lines,
shape of continuum, etc., can provide important information about their formation and evolution.
In particular, emission lines probe the gaseous and
chemical components of galaxies. The low-ionization nuclear emission-line regions
(LINER) were described by \citep{Heck80} as a class of extragalactic objects
with optical spectra dominated by enhanced low-ionization OI($\lambda 6300$) and
NII($\lambda 6548,6583$) lines. Thus, LINER were defined by intensity ratios of
optical emission lines, namely: (1) I([O II] $\lambda$ 3727)/I([O III] $\lambda$ 5007) $\geq$ 1 where [O II]
$\lambda$ 3727 is  used to designate the [O II] $\lambda$ $\lambda$ 3726, 3729
doublet, and (2) I([O I] $\lambda$ 6300) /I([O III] $\lambda$ 5007) $\geq$ 1/3.

The nature of the ionization source that powers emission lines has not been
determined. In the last three decades different scenarios for the LINER excitation
mechanism have been proposed. The most feasible are: (1) ionization by shock-heated
gas due to direct mechanical energy input of turbulent gas motions, jets or bubbles
\citep{Heck80,DopSut95}; (2) stellar photo-ionization by hot O stars \citep{FilTer92}
or by old post-asymptotic giant branch (AGB) stars \citep{Bin94,Stas08}; and (3)
photo-ionization by an active central black hole \citep{Groves04}, such as an
Active Galactic Nuclei (AGN).

The first mechanism, i.e. gas heated by fast shock waves, is probably the less
favoured scenario because the velocity dispersion of nuclear gas commonly falls
below the value required to explain the observed level of spectral ionization
\citep{Ho03}. However, there is further controversy related to the last two
mechanisms proposed. Observations at radio \citep{nagar05} and X-ray wavelengths
\citep{GM09a} provide strong support for an AGN as the origin of LINER emission.
LINER seem to populate the low luminosity end of the AGN distribution \citep{Kewley06},
where radiatively inefficient accretion flows and external obscuring matter
\citep{GM09b,dudik09} may cause the optical extinction. Thus, LINER are less
luminous than Seyfert galaxies and share similar spectral characteristics, with
the notable exception that LINER show enhanced low-ionization OI($\lambda 6300$)
and NII($\lambda 6548,6583$) lines \citep{Heck80}. Furthermore, there is new
evidence supporting the hypothesis of stellar photo-ionization, as some authors
suggest that hot post-AGB stars and white dwarfs could provide enough
ionization to explain the LINER emission \citep{Bin94,sodre99,Stas08}. These stellar
ionization sources are located in spatially extended regions around the nucleus
as the $H\alpha$ and $H\beta$ brightness profiles do not decrease with $r^{-2}$,
i.e. the radial dependency expected for a central ionization source such as an AGN
\citep{sarzi08,YanB12}.

LINER seems to be a common phenomena and they have been found in about 30\%
of nearby galaxies \citep{Heck80,HoFilSar97}. This percentage increases for galaxies
with early type morphology, reaching 50\% in elliptical passive galaxies
\citep{Fil86,Goud94,Yan06,CapB11}. This early galaxy type is consistent with that
residing in dense galaxy environments. The well known morphology-age-density relation
\citep{Dress80,DML01} shows that red and massive early-type galaxies are associated
to rich galaxy structures preferably located in the cores of groups or clusters
of galaxies. On the other hand, late-type and blue galaxies are typically found
in the outer regions of galaxy groups or in the field.

Since galaxy evolution correlates with different environmental factors, the
study of the neighbourhood of LINER could provide important clues about the
nature of their low ionization emission. If LINER are dominated by star-forming
HII regions we would expect them to follow the morphology-density segregation \citep{Dress80}.
Instead, several recent studies have found that AGN environment do not seem to
follow this relation \citep{PB06,coldwell09,padilla10,coldwell14}. Then, if the
presence of an AGN is the main cause of the LINER emission we could assume that
LINER will follow a similar trend that, for example, Seyfert 2 galaxies.

In addition, it is not clear that LINER constitute a homogeneous class of objects.
Unlike galaxies dominated by star-forming HII regions or Seyferts, whose ionization
sources are clearly identified as young massive stars and AGN, respectively,
LINER can be produced by a wide array of ionization mechanisms. For example, HST
observations have found that both unresolved nuclear and extended H$\alpha$ emission
are present in the majority of nearby LINER \citep{masegosa}, and LINER-like
ratios have long been observed in extended emission line regions \citep{Heck89,sarzi08}.
In addition, \cite{graves07} found that red-sequence galaxies with LINER emission
are younger than their quiescent counterparts, suggesting a connection between the star
formation history and the mechanism generating the low ionization emission. But, there
is no certainty about which mechanism prevails in LINER galaxies.
So the true power source of LINER emission is still under debate and most observations
do not have resolution high enough to detect if the this emission comes from extended
or nuclear regions of the galaxies.

The aim of this paper is to study the relation between LINER objects and their
local density environment in order to shed light on the origin of their main
emission mechanism. The layout of this paper is as follows: in Section 2 we
briefly describe the data selection and LINER classification scheme; in Section 3
we describe the properties of galaxies hosting LINER, the procedure used to
construct the control sample and their correlation with galaxy group parameters.
In Section 4 we investigate the occurrence of LINER in galaxy group and its dependency
with group proximity. Finally, in Section 5 we discuss the results and draw our
conclusions. Throughout this paper, we have assumed a $\Lambda$-dominated
cosmology, with $\Omega_{m} = 0.3$, $\Omega_{\lambda} = 0.7$ and
$H_0 = 100 {\rm km\, s^{-1}} Mpc^{-1}$.

\begin{figure*}
\includegraphics[height=60mm,width=160mm]{fig1.ps}
\caption{BPT diagrams of the selection criteria defined by Kewley {\rm et al. } (2006)
used to classify emission-line galaxies as Seyfert or LINER.
The panels show the line ratios $\log([\rm OIII]/\rm H\beta)$ vs $\log(\rm [NII/H\alpha])$ (left),
$\log([\rm OIII]/\rm H\beta)$ vs $\log(\rm [SII/H\alpha]$ (central) and
$\log([\rm OIII]/\rm H\beta)$ vs $\log(\rm [OI/H\alpha]$ (right).
Seyferts are indicated by black dots and LINER galaxies by grey crosses. The solid lines separate
star-forming galaxies from AGN and the dashed lines represent the Seyfert-LINER demarcation.}
\label{fig1}
\end{figure*}

\section{Data and Sample Selection}
In this paper we used galaxy samples selected from the spectroscopic Sloan
Digital Sky Survey, Data Release 7 (SDSS-DR7, \cite{abazaja09}). The SDSS-DR7
main galaxy sample is essentially a magnitude limited catalog \citep{petro76}
with \textit{$r_{lim}$}$ < 17.77$, consisting of $\sim 700000$ galaxies with
measured spectra and photometry in five optical bands ({\it u,g,r,i,z}), and
with a median redshift of $0.1$ \citep{strauss02}.

The procedures to estimate the physical properties of galaxies used in this analysis
are described by \cite{brinch04},\cite{tremonti04} and \cite{blanton05}. This data is
available from MPA/JHU\footnote{http://www.mpa-garching.mpg.de/SDSS/DR7/} and
NYU\footnote{http://sdss.physics.nyu.edu/vagc/}, including gas-phase
metallicities, stellar masses, indicators of recent major star-bursts, current
total and specific star-formation rates, emission-line fluxes and S\'ersic indices
among others.

\subsection{LINER Selection}
\label{sec:sel}

For the AGN selection we use the publicly available emission-line fluxes, whose
measurement is detailed in \cite{tremonti04}. Additionally, we have corrected
for optical reddening using the Balmer decrement and the \cite{calzetti00} dust
curve. We assume an $R_V=A_V/E(B-V)=3.1$ and an intrinsic Balmer decrement
$(H\alpha/H\beta)_{0}=3.1$ \citep{OM89}. Since the true uncertainties of the
emission-line measurements were underestimated, the signal-to-noise ($S/N$)
of every line was calculated with the flux errors adjusted as suggested by the MPA/JHU
team\footnote{http://www.mpa-garching.mpg.de/SDSS/DR7/raw$\_$data.html}.

The emission-line galaxy sample was restricted to have a redshift range of
$0.04 < z < 0.1$. The lower limit prevents that small fixed-size apertures
affect galaxy properties derived from the fibre spectra, and the upper limit
corresponds to the luminosity completeness limit for the SDSS sample. Furthermore,
we only include galaxies with  $S/N > 2$ for all the lines involved in the three
diagnostic diagrams used to discriminate Seyfert 2 from star-forming and LINER
galaxies. This conservative criteria only moderately reduces the effective sample
size, yet assures a more reliable selection of LINER objects. Thus, from this
refined sample, we separate Seyfert 2, LINER and star-forming galaxies using
the three standard \cite[][: BPT]{BPT81} line-ratio diagrams. The AGN/starburst
separation, as suggested by \cite{Kewley01,Kewley06}, depends on the relative
source location within BPT diagnostic diagrams, and follows the equations

\begin{equation}
\log([\rm OIII]/\rm H\beta) > 0.61/(\log(\rm [NII/H\alpha])-0.47)+1.19,
\label{eqn:bpt1}
\end{equation}

\begin{equation}
\log([\rm OIII]/\rm H\beta) > 0.72/(\log(\rm [SII/H\alpha])-0.32)+1.30,
\label{eqn:bpt2}
\end{equation}

\begin{equation}
\log([\rm OIII]/\rm H\beta) > 0.73/(\log(\rm [OI/H\alpha])+0.59)+1.33,
\label{eqn:bpt3}
\end{equation}

LINER are located below the Seyfert/LINER division line as indicated by equations

\begin{equation}
\label{eqn:bpt4}
\log([\rm OIII]/\rm H\beta) < 1.89\log(\rm [SII/H\alpha])+0.76,
\end{equation}

\begin{equation}
\label{eqn:bpt5}
\log([\rm OIII]/\rm H\beta) < 1.18\log(\rm [OI/H\alpha]+1.30.
\end{equation}

On the grounds of internal consistency, we exclude from analysis ambiguous
galaxies classified as one type of object in one BPT and another in the remaining
two diagrams. Bearing this in mind, we obtain an effective sample of 4778 LINER objects.
The discriminated samples and the selection criteria are shown in the three BPT
diagrams of Figure \ref{fig1}.

\subsection{Galaxy Group Catalogue}
In this work we use the group catalogue constructed by \cite{Zap09}, which were
identified by an improved version of the Huchra \& Geller (1982) friends-of-friends algorithm,
with variable linking lengths of $D_{12}=D_0 R$ and $V_{12}=V_0 R$. These linking lengths are
taken in the direction perpendicular and parallel to the line-of-sight, respectively. The spatial
scaling $R$ takes into account the variation in the space density of galaxies in a
flux-limited catalogue, and is calculated from the ratio of the density of galaxies brighter than
the minimum required to enter the catalogue at the mean distance of the galaxies being linked, to
a characteristic survey depth. As found by \cite{MerZan05}, the algorithm parameters
$D_0=0.239$h$^{-1}$Mpc and $V_0=450$kms$^{-1}$, are selected to produce a reasonably complete
sample (95\%) with low contamination ($<$ 8\%).

In addition, the virial mass of galaxy groups are computed by means of the virial theorem, given by

\begin{equation}
\rm M_{vir} = \frac{\rm 3\sigma_v^2 \rm r_{vir}}{G},
\label{eqn:mvir}
\end{equation}

where $\sigma_v$ is the line-of-sight velocity dispersion and $r_{vir}$ is estimated as in \cite{MerZan05}, i.e.

\begin{equation}
\rm r_{vir} = \frac{\pi}{2}\frac{N_g(N_g-1)}{\sum_{i>j}r_{ij}^{-1}},
\label{eqn:rvir}
\end{equation}

where $N_g$ is the number of galaxy members and $r_{ij}$ are the relative projected distances
between galaxies.

\section{Host galaxy properties}

The morphology of a galaxy can be indicative of evolutionary processes it has
been subject to, and therefore an indicate of its environment. Spiral galaxies
are known to be relatively unevolved galaxies residing typically in low
density environments. In the other hand, elliptical galaxies are known to be
more evolved sources, as a result of dynamical processes such as mergers, and known to inhabit high density
environments like galaxy groups or clusters \citep{toomre, dressler97}, where
they are subject to different effects such as strangulation, ram-pressure and
close-encounters.

The LINER host galaxies have particular characteristics which could give clues
about their history. Figure \ref{dist} shows the normalized distributions of:
$a)$ Redshift, $b)$ (Petrosian) r-band absolute magnitude \citep{blanton03},
extinction and k-corrected bandshifted to $z=0.1$, using the software \texttt{k-correct\_v4.2} of \cite{Blanton2007}.
$c)$ stellar mass in logarithmic scale, $M^{\ast}$,
previously determined by \cite{kauff03b} \footnote{The method relies on spectral
indicators relating to the stellar age and the fraction of stars formed in
recent bursts.}, $d)$ $\rm D_n(4000)$ break index \citep{kauff02} $e)$
S\'ersic index \citep{sersic63}, and $f)$ colour ($M_g-M_r$).

The break index $\rm D_n(4000)$ is defined as the ratio of the average flux
density in the narrow continuum bands ($3850-3950$ and $4000-4100$ $\AA$) and is
suitably correlated to the mean age of the stellar population in a galaxy and can be
used to estimate the star formation rate \citep{brinch04}.
The majority of star formation takes place preferentially in galaxies with low
$\rm D_n(4000)$ values.

The change of surface brightness with distance from the galaxy centre can be well
described by the S\'ersic Law, which has a the form $\ln I(R)= \ln I_{0} -kR^{1/n}$.
Here, $I_{0}$ is the intensity at the centre and the S\'ersic index $n$ controls
the degree of curvature of the profile. Setting $n = 4$ recovers a de Vaucouleurs profile,
i.e. a good description of giant elliptical galaxies; while setting $n = 1$ recovers
the Freeman exponential profile, which is a good description of the light
distribution in disk galaxies and dwarf ellipticals. Most galaxies are well fit by S\'ersic
indices in the range $0.5 < n < 10$.

\begin{figure*}
\includegraphics[width=155mm,height=129mm]{fig2.ps}
\caption{Normalized distributions of galaxy properties for: SDSS-DR7 galaxies
(dotted line), LINER total sample (solid line), and galaxy group members (dashed line).}
\label{dist}
\end{figure*}

Figure \ref{dist} also shows the distributions corresponding to the SDSS-DR7 main
galaxy sample (within the same redshift range). The bi-modality in the distribution of
$\rm D_n(4000)$ and colour for the main galaxy sample is significantly noticeable.
This effect is explained by the presence of two populations of galaxies in
different stages of evolution: a population of young and blue galaxies, and another
composed of old and red galaxies. It is expected that the mean and variance of
both colour distributions are dependent of the luminosity function and stellar mass
\citep{Bernardi03,blanton03b,Hogg04} of the population considered. Moreover,
Fig.\ref{dist} shows that most LINER are hosted by a single population of red
and old galaxies. Consistently, these LINER hosts are significantly more massive
and luminous than the general galaxy sample and have a bulge-type morphology.

These characteristics of LINER are in agreement with the results of \cite{Kewley06}
and \cite{CapB11}, and are consistent with very well known features of galaxy
groups and cluster galaxy members \citep{Dress80, Dress84}.
We include in Figure \ref{dist} the distributions of SDSS-DR7 galaxy groups
as constructed by \cite{Zap09}, noting that
galaxy group members show similar trends as the SDSS main sample, but with a
higher percentage of red and old galaxies. This means that galaxy groups have
a population of blue and young spiral galaxies, probably located in their outer
regions, and another population of bright, red and bulge-type galaxies found
in the centre of the systems. This is clear consequence of the morphological
segregation \citep{DML01} of galaxies in clusters.

In this scenario, LINER galaxies can be thought as galaxies lying in the central
area of rich galaxy groups and clusters. However, after cross-correlating the LINER
catalogue with the catalogue of galaxy groups we found that
only 12.4\% (592 objects) are members of galaxy groups
with more than 10 members, within the redshift interval
considered. Instead of that, 20.4\% of LINER (974 objects) have been found to be a member of
galaxy groups with a range of 4 to 9 members. This low percentage of LINER belonging to galaxy
groups richer than 10 members motivates us
to explore the relation between LINER properties and the group density environment.

\begin{figure*}
\includegraphics[width=175mm,height=135mm]{fig3_1.ps}
\caption{Normalized distributions of galaxy properties for LINER (solid line) and control sample (dashed) belonging
to galaxy groups within $r_p < 2$ $r_{vir} $ and  $\Delta V < 1000$ ${\rm km\, s^{-1}}$. The inner box in panel (b)
corresponds to the luminosity distribution of the brightest member of galaxy groups
matched in  both samples.
}
\label{dist2}
\end{figure*}

\subsection{LINER close to Galaxy Groups and Construction of a control sample}
\label{sec:sel}
Galaxy groups are the most common structures in the universe and the connection
between LINER and high density environments could provide valuable information
about the mechanisms responsible of their low-ionization emission.

Considering the low occurrence rate of LINER in groups, we adopted as LINERS
belonging to galaxy groups as those objects with projected distances $r_p < 2$ $r_{vir} $,
where $r_{vir}$ is the virial radius, and radial velocity difference $\Delta V < 1000$ ${\rm km\, s^{-1}}$,
with respect to the geometric centre of the group.
Thus, with this more relaxed criteria than only to be a galaxy group member  we are
considering galaxies within the dark matter halo of the groups and, also, increasing
the sample to obtain more confident statistical results.
In this way finally we obtained a sample of 3967 LINER galaxies (83\% of the total sample)
hosted by galaxy groups with more than 4 members.

To study the true environmental dependence of LINER properties it is very important
to select an appropriate control sample of galaxies without low-ionization emission
lines that belong to galaxy groups and fulfils the selection criteria as mentioned
before. The type of control samples have been used in a series of papers
\cite{coldwell03,coldwell06} and \cite{coldwell09,coldwell14} to understand the
behaviour of AGN with respect to non-active galaxies. In addition, by using SDSS
mock galaxy catalogues built from the Millennium Simulation, \cite{Perez09} showed
that a suitable control sample for galaxies in pairs should be selected (at least)
with matched distributions of redshift, morphology, stellar mass, and local
density environment. This criteria is also applicable to the case of control
galaxies for LINER. Thus, we construct our control sample by selecting
galaxies without low-ionization emission features from SDSS with matched distributions
of those five parameters.

It is important to mention that the control sample has been taken from a galaxy
sample satisfying the criteria to be closer than $r_p < 2$ $r_{vir} $ and
$\Delta V < 1000$ ${\rm km\, s^{-1}}$ with respect to the galaxy group centre. Moreover, any
galaxy classified as AGN or LINER has not been included in the control sample,
which comprises 3841 galaxies. Figure \ref{dist2} shows the normalized
distribution of properties for LINER host galaxies within the dark matter halos
of galaxy groups, and their respective control sample. We can appreciate that
the control sample suitably reproduces the properties of LINER objects. Both
samples are redder, older and more massive than the rest of the galaxies in the
SDSS catalogue, and have a S\'ersic index arithmetic mean of $~n = 4$.

In addition, we analysed the fraction of both LINER and control samples which correspond to the
central galaxy of groups, identified as the brightest member galaxy. We found a similar percentage of central
galaxies for control ($16.4\%$) and LINER ($13.8\%$) samples. The inner box in Figure 3(b) shows that the
luminosity distribution of the central galaxies in both samples is very similar. Moreover, in Figure \ref{vmass3}
we also compare the virial mass distribution of groups close ($r_p< 2$ $r_{vir} $ and $\Delta V < 1000$ ${\rm km\, s^{-1}}$)
to galaxies of our LINER and control samples. The agreement between both distributions shows that, on average,
LINER and control galaxies belong to galaxy groups of similar virial mass. As we will see in Section 4, this
supports the idea that the difference in their fractions (Fig. 6 and 7) could be driven mainly by the
presence of low-ionization emission in LINER hosts.

\begin{figure}
\includegraphics[width=90mm,height=90mm,]{fig4.ps}
\caption{Virial mass distributions for galaxy groups closer than $r_p < 2$ $r_{vir} $ and $\Delta V < 1000$ ${\rm km\, s^{-1}}$
to LINER (solid line) and control sample (dashed line).}
\label{vmass3}
\end{figure}

\subsection{Host galaxy properties: Dependency with the galaxy group proximity}

\begin{figure}
\includegraphics[width=90mm,height=90mm]{fig5.ps}
\caption{Fraction of red ($M_g-M_r > 0.89$) and old ($D_n(4000)> 1.99$) galaxies
as function of the normalized projected distance to the group centre (left panels)
and to the galaxy group virial mass (right panels). The solid line corresponds to
the LINER sample and the dashed line to the control sample.}
\label{lummass3}
\end{figure}

The similarity between LINER and control sample properties is what enables to explore
how these are related with the galaxy group parameters with the aim to detect some
distinctive feature for LINER. Particularly,  the galaxy colours provide an indirect
constraint on the evolutionary history of galaxies since other galaxy parameters
for example morphology, age, affect galaxy colours, resulting very appropriated for
this study.

As in previous works \citep{coldwell06,coldwell09,coldwell14}, we use colours and
the age indicator to study the fraction of host galaxies redder than
$\rm M_g-M_r > 0.89$ and $D_n(4000) > 1.99$, as a function of the galaxy group
centre distance and the virial mass. These thresholds for the fractions correspond
to mean values of $\rm M_g-M_r$ and  $D_n(4000)$ for the LINER and control
samples (see the distributions in Figure \ref{dist2}). The error bars are calculated with the bootstrap error
resampling technique \citep{barrow84}.

It is possible to observe a slightly higher fraction of red and old galaxies closer to
the geometric group centre, as shown in Figure \ref{lummass3}. In addition, the trends
of LINER and control samples are nearly indistinguishable.
 Moreover, other than a weak excess of red galaxies
(still within the error bars), we do not find any significant difference with
respect to the group virial mass. The conclusion is that the properties of LINER
host galaxies at a given distance or virial mass do not vary respect to control
non-active galaxies, so the location within the galaxy group environment does not
seem to affect colours and/or ages of LINER.

\section{Occurrence of LINER in galaxy groups}

In this section we explore the probability to find LINER respect to the control
sample and the dependence with galaxy groups properties such as the virial mass,
$M_{vir}$ (as derived from the virial theorem using the virial radius and
velocity dispersion of group members). We also considered the group luminosity,
$M_{r4}$, calculated by adding the r-band luminosities of the 4 brightest galaxy
members in each group. This quantity, $M_{r4}$, is a better proxy for the true underlying
group mass than the virial estimate \citep{Eke04,padilla04}. In this way, we
have considered two independent estimates of the group richness for this analysis.

Figure \ref{lummass} shows the fraction of both samples with respect to the
total number of objects for a given bin of virial mass in logarithmic scale
and group luminosity. We observe that the fraction of LINER
drops with the group virial mass and luminosity, while the fraction of control
galaxies increases with these parameters. Precisely for control sample, this is
expected due to the morphology-density relation and the sample distributions
shown in Figure \ref{dist2}. Yet somewhat surprisingly, the relation for LINER
galaxies is in fact the opposite, even when their hosts galaxies are matched to the
control sample.

This tendency is more significant in the left panel of Figure \ref{lummass}, where
the fraction of galaxies in the brightest groups in the control sample is
approximately 2 times the LINER fraction. At low luminosities, the effect is
is reversed as the proportion of LINER reaches its maximum value ($\sim$56\%)
and the difference respect to control falls to $\sim$13\%.
A comparable trend is observed in the right panel of Fig. \ref{lummass}, with
respect to $M_{vir}$, although the excess of control galaxies respect to LINER at
higher values of $M_{vir}$ is somewhat lower than that shown for the group luminosity,
$M_{r4}$.
These results indicate
that LINER galaxies have a preference to populate low mass, low luminosity
environments formed by a higher proportion of blue galaxies with high gas content.

\begin{figure*}
\includegraphics[width=155mm,height=95mm,]{fig6.ps}
\caption{Fraction of LINER (solid line) and the respective control sample (dashed lines) as a function
of the group virial mass $M_{vir}$ (right), and brightness $M_{r4}$ (left).}
\label{lummass}
\end{figure*}

Bearing in mind the fact that several authors have found that the association
between active galaxies and environment do not agree with the expected
morphology-density relation \citep{PB06,coldwell09,padilla10,coldwell14},
these results suggest that the accretion of material inside the
black hole should not be neglected as the main mechanism responsible of
the low-ionization emission.

\begin{figure*}
\includegraphics[width=155mm,height=95mm,]{fig7.ps}
\caption{Fraction of LINER (solid line) and the respective control sample (dashed lines)
as a function of the normalized distance to the centre of galaxy groups with virial
mass $\rm log_{10}(M_{vir}) < 13.8$ (left) and $\rm log_{10}(M_{vir}) \geqslant 13.8$ (right).}
\label{lummass4}
\end{figure*}

Additionally, it is useful to perform this analysis by discriminating groups
with a group mass indicator. Figure \ref{lummass4} shows the fraction of LINER
and control galaxies as a function of the distance to the galaxy group centre,
normalized to the virial radius, $d_{cg}/r_{vir}$, for rich groups ($log_{10}M_{vir} > 13.8M_{\sun}$)
and for poor groups ($log_{10}M_{vir} < 13.8M_{\sun}$). In poor groups we observe a
higher fraction of LINER respect to the control sample, which is not unexpected
given Figure \ref{lummass}. This effect is almost independent of scale with only
a weak increase in the fraction of LINER galaxies toward the inner regions of
poor galaxy groups. In rich groups, we find that the fraction of LINER increases
significantly toward the outer limits of these systems and that the fraction of
LINER objects is lower that in control galaxies.
Therefore, this reinforces the hypothesis suggested, that LINER galaxies prefer
to populate low density environments, where the probability of major mergers rises
due to low galaxy velocity dispersions favouring this type of interactions.
This seems to be in agreement with previous results for AGN found in the literature
\citep{PB06,alonso07,coldwell09}.

\section{Discussion}

From SDSS-DR7 data we derived a large LINER catalogue at redshifts $0.04 < z < 0.1$.
As only a small fraction of LINER are hosted in groups of galaxies,
we adopted LINER objects with $r_p < 2 r_{vir}$ and $\Delta V < 1000 {\rm km\, s^{-1}}$ as belonging
to galaxy groups. In this way, our effective LINER sample comprises 3967 objects,
distributed in 2879 groups. The mean number of LINER per group
is $\sim 1.4$. In order to uncover the true environmental dependence
of the characteristics of galaxies hosting LINER and get a hint
on the process responsible for generating low-ionization emission lines,
we constructed a control sample of non-active galaxies matched in redshift,
luminosity, age, stellar mass, morphology and colour.
We also verified that : (1) the percentage of LINER and control galaxies corresponding to
the central members of groups is quite similar between both samples; (2) galaxy groups nearby
both samples have similar virial masses.

Following previous work, we have explored how host properties depend on mass and the proximity
to galaxy groups. We have calculated the fraction of red ($M_g-M_r > 0.89$) and old
($D_n(4000)> 1.99$) galaxies respect to the normalized projected distance,
$d_{cg}/r_{vir}$, to the geometric group centre, and to the virial mass, $M_{vir}$.
We find no significant difference in the tendencies of LINER and control samples. Therefore,
the properties of hosts in both samples do not seem to vary respect to the location or mass of
the galaxy groups.

In addition, we have carried out a detailed analysis the LINER occurrence within galaxy group
halos and quantified the colour and age dependency with respect to representative
galaxy group parameters such as radius, virial mass and group luminosity. We study
the fraction of LINER with respect to two independent parameters representative
of the galaxy group mass ($M_{vir}$ and $M_{r4}$). The results show a strong
difference between LINER and control galaxies. While the latter fraction
increases for rich groups, the fraction of LINER drops drastically. This
effect is more significant for the highest luminosity groups, where we find about
2 times more control than LINER galaxies, indicating an evident preference of
LINER to inhabit low mass galaxy groups.
Besides, LINER galaxies do not seem to follow the expected morphology-density
relation in the proximity of rich galaxy groups.
Such a behaviour of LINER agrees with
that corresponding to active objects (Seyfert 2, quasars, etc.) analysed in
previous works. So, although we have no certainty about
the degree at which the mechanisms are responsible of the low-ionization
emission, the presence of nuclear activity is certainly
involved and should not be discarded.

In conclusion, we have shown the higher probability of LINER to be found in lower density
environment, such as poor groups, where the galaxy interactions and content of gas are more abundant.
Thus, the global conditions of these environments could favour the presence the LINER emission.

\section{Acknowledgments}
We would like to thanks to anonymous referee for the comments that helped to improve the paper.
This work was supported in part by the Consejo Nacional de
Investigaciones Cient\'ificas y T\'ecnicas de la Rep\'ublica Argentina
(CONICET), the Consejo Nacional de Investigaciones Cient\'ificas, T\'ecnicas y de Creaci\'on Art\'istica de la
Universidad Nacional de San Juan (CICITCA) and the Secretaría de Estado de Ciencia, Tecnolog\'ia e Innovaci\'on
del Gobierno de San Juan (SECITI).

Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan
Foundation, the Participating Institutions, the National Science Foundation,
the U.S. Department of Energy, the National Aeronautics and Space Administration,
the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education
Funding Council for England. The SDSS Web Site is \emph{http://www.sdss.org/}.
The SDSS is managed by the Astrophysical Research Consortium for the Participating
Institutions. The Participating Institutions are the American Museum of Natural
History, Astrophysical Institute Potsdam, University of Basel, University of
Cambridge, Case Western Reserve University, University of Chicago, Drexel University,
Fermilab, the Institute for Advanced Study, the Japan Participation Group, Johns
Hopkins University, the Joint Institute for Nuclear Astrophysics, the Kavli Institute
for Particle Astrophysics and Cosmology, the Korean Scientist Group, the Chinese
Academy of Sciences (LAMOST), Los Alamos National Laboratory, the Max-Planck-Institute
for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New Mexico
State University, Ohio State University, University of Pittsburgh, University of
Portsmouth, Princeton University, the United States Naval Observatory, and the University
of Washington.