Constraining electroweak and strongly charged long-lived particles with CheckMATE

Abstract: Long-lived particles have become a new frontier in the exploration of physics beyond the Standard Model. In this paper, we present the implementation of four types of long-lived particle searches, viz. displaced leptons, disappearing track, displaced vertex (together with muons or with missing energy), and heavy charged tracks. These four categories cover the signatures of a large range of physics models. We illustrate their potential for exclusion and discuss their mutual overlaps in mass-lifetime space for two simple phenomenological models involving either a U(1)-charged or a coloured scalar.

\section{Introduction}

The Large Hadron Collider (LHC) started operations over a decade ago and a large number of searches for physics beyond the Standard Model (BSM) have been performed.  In particular, the general-purpose experiments ATLAS and CMS have published over a thousand papers each, minutely testing the predictions of the Standard Model and of more exotic theories. For the simplest cases of new physics consisting of new strongly charged particles, corresponding searches already place limits on the BSM masses at about 2-3 TeV \cite{Aad:2020aze}, closing in on the maximum reach achievable at a 13 TeV proton collider.  We turn our attention therefore to the equally well-motivated but experimentally more challenging cases of new physics where BSM particles are long-lived.  A long lifetime may be the natural consequence of a compressed phase space (e.g.\ in particular dark matter models \cite{Ibe:2006de,Cirelli:2005uq,FileviezPerez:2008bj,Cirelli:2009uv,Buckley:2009kv,Johansen:2010ac,Mahbubani:2017gjh}), a suppressed connection to light SM decay products caused by heavy mediators \cite{Giudice:2004tc}, or direct but feeble couplings to SM particles (possibly resulting from an approximate symmetry) \cite{Hall:2009bx,ArkaniHamed:2008qp,Graham:2012th,Zwane:2015bra}.  These new kinds of scenarios have raised considerable interest in the community~\cite{Alimena:2019zri} and making the results of associated searches available to reinterpretation in terms of other theory models has been a frequent request to the {\tt CheckMATE} collaboration.

With the end of Run 2 of the LHC, both ATLAS and CMS have turned their attention to searches for exotic physics based on new particles with long lifetimes. The signatures of a long-lived  particle (LLP) depend on its charges as well as its  decay modes and actual lifetime, and they can be fairly complicated to systematically characterize.  Several possibilities for LLP signatures indeed exist, based on the LLP decay modes:  \begin{itemize}
    \item neutral LLP $\rightarrow$ invisible or neutral stable particles $\Rightarrow$ Missing momentum
    \item neutral LLP $\rightarrow$ charged leptons $\Rightarrow$ leptons with large impact parameter (i.e. ``displaced'' leptons)
    \item neutral LLP $\rightarrow$ coloured particles $\Rightarrow$ displaced vertices, or emerging jets
    \item stable, charged LLP $\Rightarrow$ charged track (with its time of flight dependent on mass and boost)
    \item charged LLP $\rightarrow$ invisible $\Rightarrow$ ``disappearing'' track
    \item charged LLP $\rightarrow$ other charged stable object(s) $\Rightarrow$ kink-track or displaced vertex.
\end{itemize}
There is a built-in complementarity in different searches simply because particle decay follows an exponential distribution. For example, charged LLPs with intermediate lifetimes will be visible in both disappearing track and heavy charged track searches.  Similarly, a neutral particle decaying into quarks mostly in the electromagnetic or hadronic calorimeter will also likely appear as a smaller, simultaneous signal in the displaced vertex (decay in the tracker) and emerging jet (decay in the hadronic calorimeter) searches.
 Furthermore, the lifetime of a particle in the lab-frame also depends on its boost, which means that the production mechanism can also significantly alter where in the detector the particle decays.  The same particle may accordingly result in different signal distributions depending on whether it is produced ``directly'' or in the decays of a much heavier particle (i.e.\ with a higher boost).  Finally, several decay modes may be open to a single new particle resulting in sensitivity in multiple searches.  The identification of the underlying physics therefore requires a full coverage in terms of the lifetime of new particles.

As we can see, the identification of an LLP is highly complicated and so far, there are no standard algorithms like those available for the identification of standard objects, such as leptons, $b$- or $\tau$-tagged jets, etc. Consequently, it is not always clear how the results of a dedicated LLP search can be  ``reinterpreted'' for a physics model that differs from the tested one, though it displays a priori similar signatures.  A detailed study of models capable of LLP signatures, the reinterpretation struggles and recommendations have been detailed in the community study~\cite{Alimena:2019zri}.  In this present work, we use the signal efficiencies published by the experiments  in order to implement five searches in the {\tt CheckMATE} reinterpretation package.

{\tt CheckMATE} \cite{Drees:2013wra,Dercks:2016npn} is a public tool that allows the reinterpretation of a wide variety of ATLAS and CMS results for new physics models in a coherent and cohesive manner. It consists of an engine written in C++ that runs each analysis cut-by-cut in order to assess the final number of expected events satisfying the requirements of the corresponding analysis.  The engine is also capable of using external libraries like Madgraph~\cite{Alwall:2014hca} and Pythia8~\cite{Sj_strand_2015}\footnote{Internal running of Monte Carlo processes or use of Madgraph interface in {\tt CheckMATE} is currently compatible only with the 8.2 series of Pythia8.} in order to generate events, while the detector simulation is performed by  Delphes~\cite{de_Favereau_2014}. The User Interface and the statistical analyses are provided by a collection of Python scripts (including the AnalysisManager~\cite{Kim:2015wza} that guides the users through the implementation of their own analyses).

In section~2, we briefly summarize the main ingredients of the LLP recast, referring the reader to the appendix for a more complete description. Then, in section~3, we illustrate the performance of the implemented searches in two simple models with LLPs, and discuss their complementarity. Conclusions and a brief outlook are proposed in section~4. The appendix consists of a short guide for the user, as well as a more detailed presentation of the implemented LLP searches.

\section{Implementation of long-lived particle searches}

 Below we offer a brief description of the implemented LLP searches and a comparison with experimentally published results.  This includes the 8 and 13 TeV versions of the CMS displaced lepton search \cite{Khachatryan_2015,CMS-PAS-EXO-16-022}, two different displaced vertex searches \cite{Aaboud:2017iio,Aad:2020srt}, the 13 ATLAS disappearing track \cite{Aaboud:2017mpt} and heavy charged particle track \cite{Aaboud:2019trc} searches.  Together, these searches are capable of probing a wide range of parameter space. Details of the implementation are available in the appendix.  Technically, each analysis is encapsulated in a detector-specific ``analysis handler'' class which provides special functions and efficiencies specific to the detector in question.  We deliberately separate the analysis handlers for long-lived particle searches from those used in prompt searches: this accounts for the fact that the implemented prompt searches do not in fact use any decay length information and all particles denoted stable by the Monte Carlo generator are then clustered into jets based on their kinematics only.\footnote{For this reason, we advise user discretion when applying prompt search limits to models with LLPs.  It may be possible to find a conservative limit from prompt searches by e.g.\ removing decay products of LLPs from the Monte Carlo events beforehand. However each case needs to be evaluated separately and we do not provide a built-in solution for that reason.}

\subsection{Displaced Lepton searches \label{sec:DLS}}

The displaced lepton searches \cite{Khachatryan_2015,CMS-PAS-EXO-16-022} look for two high-$p_T$, isolated leptons ($\ell$) with large impact parameter relative to the primary vertex.  The benchmark used for this search is motivated by R-parity violating (RPV) supersymmetry \cite{Barbier:2004ez,Graham:2012th} where a top-squark ($\tilde{t}_1$) decays via the lepton-number-violating LQD operator as $\tilde t_1 \rightarrow \ell b$.  The leptons thus produced have large $p_T$ and are well isolated.  The two searches implemented here correspond to 8 TeV \cite{Khachatryan_2015} and 13 TeV \cite{CMS-PAS-EXO-16-022} versions of the CMS displaced supersymmetry search.  The identification and fiducial acceptances are provided on generator-level events.  We therefore reproduce the Monte-Carlo production process for validation of the search. Corresponding details are provided in appendix~\ref{app:A}.

\begin{figure}[t]
\centering
\begin{subfigure}{.49\linewidth}
  \includegraphics[width=\textwidth]{DL_exclusion_8.pdf}
  \caption{8 TeV}
  \label{fig:8TeVDLex}
\end{subfigure}
\begin{subfigure}{.49\linewidth}
  \includegraphics[width=\textwidth]{DL_exclusion_13.pdf}
  \caption{13 TeV}
  \label{fig:13TeVDLex}
\end{subfigure}
\caption{A Comparison of the exclusion limits on the Displaced Lepton search provided by CMS  with those obtained from {\tt CheckMATE} (left: 8 TeV, 19.7 fb$^{-1}$; right: 13 TeV, 2.6 fb$^{-1}$). }
\label{fig:DLexclusion}
\end{figure}

The event selection for both 8 TeV and 13 TeV analyses was performed in two stages.  The first stage (i.e.``pre-selection'') selects events with exactly one electron and one muon with opposite electric charges, each expected from the decay of a different top squark.  Further selection cuts and isolation requirements are then applied.  In the second stage, the events are classified into three signal regions (SR) corresponding to increasing ranges of the leptonic impact parameter $d_0$.

The validation results are shown in Fig.~\ref{fig:DLexclusion}, with the experimental exclusion limits displayed in blue, while the recast produces the black exclusion bound. A reasonable agreement is observed at $8$~TeV in Fig.~\ref{fig:8TeVDLex}. The situation at $13$~TeV is somewhat more subtle.

Indeed, as efficiencies have been provided by the experimental collaboration for the $8$~TeV, but not for the $13$~TeV search---in particular, the considered ranges of $p_T$ and $d_0$ do not match and the modelling of efficiencies in this latter case appears as an important assumption in the recast. An attempt to validate the 13 TeV search using the 8 TeV efficiencies results in poor agreement with the numbers in the signal region at large impact parameters, and therefore in a much stronger expected limit for high values of the LLP lifetime $c\tau$ (see the dashed black curve in Fig.~\ref{fig:13TeVDLex}).  A simple linear interpolation performs rather poorly as well.  We therefore make a conservative estimate of this detector effect by adding a single bin for the $d_0$ range (20 mm - 100 mm) and determine the associated efficiency via a $\chi^2$-fit of the expected number of events in all three signal regions.  The outcome of this procedure is displayed in Fig.~\ref{fig:13TeVDLex} as a solid black line and shows a considerably improved agreement  with the expected exclusion limits.   Exact numbers in each of the signal regions are produced in Appendix~\ref{app:A}.
\FloatBarrier

\subsection{Displaced Vertex searches --- DV +
MET}\label{sec:1710.04901}

This ATLAS search \cite{Aaboud:2017iio} looks for high-mass displaced vertices (DV's), reconstructed from five or more tracks. Large missing transverse momentum is also required. The outcome of $32.8~\text{fb}^{-1}$ of $13$~TeV collision data is a yield consistent with the expected background.

The template considered by the ATLAS collaboration consists in the (strong) production of a pair of long-lived gluinos ($\tilde{g}$), then decaying into light quarks and stable neutralinos. Heavy squark mediators result in suppressed gluino decay widths. Recast instructions were provided in \cite{DVrecast} and include a pre-selection at generator-level, followed by the application of parametrized efficiencies. Previously to our implementation, this strategy has been applied with success by the publically available codes \cite{Cottincode,Lessacode} (see also Contribution~$22$ in \cite{Brooijmans:2018xbu}). Details on the implementation are provided in appendix~\ref{app:C}.

For the validation, we considered the limits on the gluino production cross-sections presented in the ATLAS paper. In a first scenario, the long-lived gluino and the neutralino LLP are separated by a wide mass-gap, with the neutralino fixed at $100$~GeV while the gluino takes mass of $m_{\tilde{g}}=1.4$~TeV or $2$~TeV. The LLP lifetime is varied between $\tau_{\tilde{g}}=0.003$~ns and $50$~ns.
The results from the recast search are shown in Fig.~\ref{fig:validation_1710_04901}, using two statistical approaches: the simplified evaluation of {\tt CheckMATE} defining a ratio $r$ (blue curve)---see Eq.~(1) of \cite{Dercks:2016npn}---and the full $p$-value analysis (red curve); both return very similar limits.
The statistical uncertainty in the simulation is below percent level ($10^6$~events are generated at each point), except for the end point $\tau_{\tilde{g}}=0.003$~ns, where it reaches $3-4\%$. The $95\%$\,CL limits from the experimental analysis is shown in black.
We observe a general qualitative agreement. In the lower row of plots of Fig.~\ref{fig:validation_1710_04901}, the experimental observed limits are normalized to the limiting cross-sections of the recast procedure (with $r$-approach in blue and $p$-values in red). Quantitatively, we find that the bounds agree within $20\%$ of the cross-section value, with outliers at up to $50\%$ discrepancy for small lifetimes.

Nevertheless, this apparent success of the recast strategy with efficiencies applied on truth-level objects needs to be tempered as it seems to perform worse in the case of a compressed spectrum. This was confirmed to us by the authors of  \cite{Cottincode,Lessacode}. A more detailed comparison is provided in the appendix.

\begin{figure}[tbhp]
\begin{center}
\includegraphics[width=0.49\textwidth]{Limits_1400.pdf}
\includegraphics[width=0.49\textwidth]{Limits_2000.pdf}\\ \vspace{-0.2cm}
\includegraphics[width=0.49\textwidth]{Limits_1400_ratio.pdf}
\includegraphics[width=0.49\textwidth]{Limits_2000_ratio.pdf}
\end{center}
\caption{Validation of the DV+MET search in the scenario with large mass-splitting for two different benchmarks (left $m_{\tilde{g}}=1.4$~TeV, right: $m_{\tilde{g}}=2$~TeV.). The bottom panel in both cases shows the ratio of the published ATLAS exclusion to the one obtained by {\tt CheckMATE}.  There is no observable difference between using the full CLs method and an exclusion based on the ratio ($r$) of cross section from {\tt  CheckMATE} of events passing all cuts to the 95\% upper limits on cross sections published by ATLAS. }
\label{fig:validation_1710_04901}
\end{figure}
\FloatBarrier

\subsection{Displaced Vertex searches --- DV + $\mu$}\label{sec:dvplusmu}

In this section we discuss a search for massive, long-lived particles decaying to final states with a DV and an energetic muon~\cite{Aad:2020srt}. The search analyzed $139~\text{fb}^{-1}$ of data collected by ATLAS at the centre of mass energy 13~TeV.

The benchmark process considered by the experiment was pair production of top squarks followed by the RPV decay into a light quark and a muon. Other physics scenarios, for example, long-lived lepto-quarks, right-handed neutrinos or long-lived electroweakinos in RPV, could result in similar signals including a DV and a muon. In Section~\ref{sec:performance} we apply this search to sbottom pair-production followed by the RPV decay.

The event selection defines two mutually exclusive trigger-based signal regions: $E_\mathrm{T}^\text{miss}$ Trigger SR and Muon Trigger SR. The former requires significant missing transverse momentum ($>180$~GeV), while the latter is recorded with the muon trigger and has low ($<180$~GeV ) transverse momentum. Additionally, at least one displaced vertex is required to be present in the fiducial region. There is no explicit requirement for a signal muon to originate from the reconstructed vertex.

The search was validated using a benchmark RPV-supersymmetric (SUSY) model for the process $pp \to \tilde{t}_1 \tilde{t}_1$, $\tilde{t}_1 \to \mu\, q$. In Figure~\ref{fig:validation_2003.11958} we show a comparison of the ATLAS result and {\tt CheckMATE} recasting in the stop lifetime-mass plane, $\tau_{\tilde{t}}$--$m_{\tilde{t}}$. The yellow band shows a 2-sigma range of the ATLAS expected exclusion limit, the blue solid line is the ATLAS observed exclusion while the blue dashed the ATLAS expected exclusion, and the black solid line shows an exclusion line obtained with {\tt CheckMATE}. Generally a good agreement is observed, however in a range of lifetimes $0.01$--$0.1$~ns, the recast exclusion is significantly weaker, though within the 2-sigma band. Further details can be found in Appendix~\ref{app:E}.

\begin{figure}[tbh!]
\begin{center}
\includegraphics[width=0.7\textwidth]{c1.pdf}
\end{center}
\caption{Comparison of exclusion limits in the stop mass and stop lifetime plane, $m_{\tilde{t}}$--$\tau_{\tilde{t}}$, reported by ATLAS: expected -- dashed blue; observed -- solid blue and {\tt CheckMATE} -- black solid. The yellow band shows a 2-sigma range of the ATLAS expected exclusion limit. \label{fig:validation_2003.11958}}
\end{figure}
\FloatBarrier

\subsection{Heavy Charged Particles searches} \label{sec:1902.01636}

In this section we focus on the search for heavy charged long-lived particles performed by the ATLAS experiment using a data sample of $36.1~\text{fb}^{-1}$ of collisions at 13 TeV \cite{Aaboud:2019trc}. In our implementation we cover searches for long-lived charginos and sleptons.

The ATLAS collaboration reported no significant excess of observed data events above the expected background in this search. Thus, the collaboration have published upper limits at 95\% confidence level on the cross-sections for stau and chargino production for specific benchmark models. These limits have been obtained applying the CLs prescription \cite{Read:2002hq}.

For the validation of our implementation in  {\tt CheckMATE} we have employed HistFitter \cite{Baak:2014wma} to estimate the CLs while running $10^5$ toy-experiments, given the low backgrounds for the considered signal regions and assuming a 10\% signal uncertainty.

The validation has been performed through the comparison with the observed upper cross-section limits reported by the ATLAS collaboration, as is depicted in Fig.~\ref{fig:validation_1902}. A very good qualitative agreement is visible between the ATLAS (red) and the {\tt CheckMATE}-derived (dashed blue) limits for both chargino and stau scenarios. More details can be found in Appendix~\ref{app:heavychargedparticle}.

\begin{figure}[tbh!]
\begin{center}
 \includegraphics[scale=0.45]{chargino.pdf}
 \includegraphics[scale=0.45]{stau.pdf}
\end{center}
	\caption{Validation of the HCP search -- observed upper cross-section limits using the best expected signal region for chargino pair production (left) and stau pair production (right). ATLAS results are shown as a solid red line while {\tt CheckMATE} predictions are in form of a dashed blue line. The theory prediction is displayed as a dot-dashed black line.}
	\label{fig:validation_1902}
\end{figure}
\FloatBarrier

\subsection{Disappearing track searches}
\label{subsect:disappearingTrack}
The ATLAS collaboration presented a search\footnote{HepData accessed from https://doi.org/10.17182/hepdata.78375.v5}~\cite{Aaboud:2017mpt} for direct electroweak (EW) gaugino or gluino pair production with wino-like electroweakinos (hence near-degenerate charged and neutral $SU(2)$-triplet fermions). The chargino decays via $\tilde \chi^+ \rightarrow  \pi^+ \tilde{\chi}^0$; the neutralino is stable. The experimental collaboration analysed data based on the integrated luminosity of 36.1 fb$^{-1}$ recorded between 2015 and 2016. Due to the small mass difference between the two states (which is of the order of the pion mass), the chargino is long-lived and the decay products are entirely invisible to the detector.  Thus the chargino appears as a ``disappearing'' track, i.e.\ a track that does not reach the outer edges of the tracker detector but stops before.  In order to define a trigger isolating the signal from large SM backgrounds, the search additionally demands a large momentum jet from initial-state radiation or four jets originating from the gluino decay. The observed number of events is consistent with the SM expectations and constraints for wino-like charginos with a lifetime of 0.2 ns,  and mass up to 460 GeV are derived. In the strong production channel, where the chargino emerges from the decay of a gluino (colour octet fermion), limits on gluino masses up to 1.65 TeV are reported, under the assumption of a chargino mass of 460 GeV and lifetime of 0.2 ns.

In order to reinterpret this search, we follow all procedures regarding production of signal events described in the original paper as closely as possible.  A detailed description can be found in appendix~\ref{app:B}.  The kinematic cuts for both signal regions are summarised in Table \ref{disappearingTrackSR}.  ATLAS further applies quality requirements, as that the tracklet is required to have hits in all four pixel layers, and a disappearance condition is demanded for each event, as the number of SCT hits associated with the tracklet must be zero. Although, the two latter cuts are impossible to simulate in a phenomenological study, ATLAS provides efficiency maps for the tracklets for EW SR and strong SR, respectively.  In addition, the collaboration provides a transverse momentum smearing function for the chargino.  We also use the benchmark SLHA files for the EW and strong scenarios, the pseudo analysis code, and all relevant data made publicly available at \textsc{HEPData} \cite{DisSearchData}.

\begin{table}[tbh]
\begin{center}
  \begin{tabular}{|l | l |}
    \hline
    EW SR & strong SR\\
    \hline
    at least one jet with $p_T>140$ GeV & $p_T(j_1)>100$ GeV, $p_T(j_2)>50$ GeV, $p_T(j_3)>50$ GeV\\
    $E_{\text{T}}^{\text{miss}}>140$ GeV & $E_{\text{T}}^{\text{miss}}>150$ GeV\\
    $\Delta\Phi(E_{\text{T}}^{\text{miss}}$, jets($p_T>50$ GeV)) $>1.0$ & $\Delta\Phi(E_{\text{T}}^{\text{miss}}$, jets($p_T>50$ GeV)) $>0.4$\\
    cuts on charged LLP with $p_T>100$ GeV & cuts on charged LLP with $p_T>100$ GeV\\
    \hline
  \end{tabular}
  \caption{Summary of the selection criteria for signal events for direct electroweakino production and the strong channel channel where the chargino is produced in gluino decays.}
  \label{disappearingTrackSR}
\end{center}
\end{table}
\begin{table}[h]
\begin{center}
  \begin{tabular}{|l | c | c | c | c | c |}
    \hline
    & CM $\tilde\chi_1^\pm\tilde\chi_1^\pm$  & CM $\tilde\chi_1^+\tilde\chi_1^0$ & CM $\tilde\chi_1^-\tilde\chi_1^0$ & CM all channels & ATLAS \\
    \hline
    \hline
    Trigger & 445.1 & 624.0 & 274,4 & 1343.5 & 1276 \\
    \hline
    Lepton Veto & 423.4 & 608.5 & 267.3 & 1308.2 & 1181 \\
    \hline
    MET and jet requirements & 164.2 & 229.6 & 101.0 & 494.8 & 579\\
    \hline
    EW SR & 5.2 & 4.4 & 1.6 & 11.2 & 13.5\\
    \hline
  \end{tabular}
  \caption{Cutflow comparison for a chargino produced in direct electroweak production with ($m_{\tilde \chi^\pm_1}$, $\tau_{\tilde \chi^\pm_1}$) = (400 GeV, 0.2 ns).}
    \label{disappearingTrackCutflowEW}
  \end{center}
\end{table}

\begin{table}[tbh!]
\begin{center}
  \begin{tabular}{| l | c | c |}
    \hline
    & CM & ATLAS \\
    \hline
    \hline
    Trigger & 289 & 285 \\
    \hline
    Lepton Veto & 277 & 278 \\
    \hline
    MET and jet requirements & 216 & 202\\
    \hline
    strong SR & 11 & 11\\
    \hline
  \end{tabular}
  \caption{Cutflow comparison for a chargino produced in strong production channel with ($m_{\tilde g}$, $m_{\tilde \chi^\pm_1}$, $\tau_{\tilde \chi^\pm_1}$) = (1600 GeV, 500 GeV, 0.2 ns) in the high--MET region.}
  \label{disappearingTrackCutflowStrong}
\end{center}
\end{table}

We use the ATLAS benchmark points as the test case scenarios which correspond to $ (m_{\tilde \chi^\pm_1}, \tau_{\tilde \chi^\pm_1}) = (400\,\rm{GeV},\, 0.2\,\rm{ns})$ for the EW signal region and $
(m_{\tilde g}, m_{\tilde \chi^\pm_1}, \tau_{\tilde \chi^\pm_1}) = (1600\,\rm{GeV},\, 500\,\rm{GeV},\, 0.2\,\rm{ns})$ in the strong signal region. The results for the EW and strong SR are summarized in Table~\ref{disappearingTrackCutflowEW} and \ref{disappearingTrackCutflowStrong}, respectively and our recast results show satisfactory agreement with the public ATLAS results.

We did not validate the disappearing track search in a grid scan with ATLAS exclusions, since the event generation requires matched events with two additional partons in the final state in order to reproduce the ATLAS cutflows. As only a few events would pass all selection cuts, such a scan would be costly to perform and we thus opted against it.
\FloatBarrier

\section{Performance and interplay in LLP scenarios\label{sec:performance}}
In this section, we consider two simple LLP scenarios that put forward the complementarity of the implemented searches.

\subsection{Electroweak LLP}

\begin{figure}
    \centering
    \includegraphics[scale=0.5]{EW_excl.pdf}
    \includegraphics[scale=0.5]{EW_excl_yuk_test.pdf}
    \caption{Upper limits in the lifetime-mass plane (left) and coupling-mass plane (right) of the LLP scalar $\phi$. We do not get limits from the disappearing track search due to a veto on leptons in the event. \label{fig:EWLLP_bounds}
    }
\end{figure}

We first consider a model addressing the case of electroweakly produced LLPs.  It extends the SM with a scalar ($\phi$), charged under $U(1)_Y$ and sharing the gauge quantum numbers of a right-handed lepton, and a SM-singlet Dirac fermion ($\chi$).  An extra $Z_2$-symmetry under which both $\phi$ and $\chi$ are charged further constrains their interactions.  The scalar is produced in pairs via the Drell-Yan process $pp \rightarrow \phi^* \phi$. After that it decays according to $\phi \rightarrow \ell \chi$, mediated by the Yukawa coupling $y_\ell \phi \bar \ell_R \chi + \mathrm{h.c}$. The singlet fermion is assumed to be stable.  The model has been implemented in Pythia8.  The Lagrangian is simply
\begin{equation}
    \mathcal{L} = \frac{1}{2}D_\mu \phi D^\mu \phi - \frac{1}{2}m_{\phi}^2 \phi^2+ \bar \chi (i\gamma^\mu \partial _\mu) \chi - m \bar \chi \chi - \sum_\ell (y_\ell \phi \bar \chi \ell_R + \mathrm{h.c.}) \, .
\end{equation}

For small values of the Yukawa coupling, the scalar $\phi$ is long-lived and could be visible in the charged track searches.  In the case where the lifetime is too short to leave the required hits in the tracker modules, the displaced lepton search might detect the products of the decay.  To examine the complementarity of these searches, we set $y_e = y_\mu$, ensuring that the targets of the displaced lepton search are actually produced. As the lifetime is determined by the smallness of the Yukawa coupling and is fairly independent of the mass of $\chi$, we set $m_\chi = 10$ GeV.  These assumptions of course mean a focus on a very specific scenario.

The exclusion results obtained with {\tt CheckMATE} in the lifetime-mass plane are shown in Fig.~\ref{fig:EWLLP_bounds}. The CMS displaced-lepton searches appear sensitive to the considered scenario for LLP lifetimes up to $\sim10^2$\,ns and masses in the range $100-500$\,GeV. The Heavy Charged Particle search from ATLAS impacts the high-lifetime area $c\tau>1$\,m, in a comparable range of masses. The disappearing track search does not show any sensitivity, due mainly to cross-sections much weaker in the considered model than the targeted range of the search: the scalar production cross section is indeed suppressed compared to the fermionic one. Moreover, the original disappearing track final state consists of multiple final states (all wino final state configurations), e.g. the NLO cross section for scalar leptons of mass 125 GeV is about 0.045 pb whereas the combined wino pair production cross section is about 1 pb.
Finally, the cross section is further suppressed since the disappearing track search basically triggers on monojet events which requires a hard recoil of the scalar pair against a hard jet. As a consequence, the current model cannot be probed with the disappearing track search which is focused on mass degenerate wino--like electroweakinos.
Still, it is obvious that two LLP searches aiming at quite different signatures can interplay and lead to complementary exclusion bounds. Current prompt limits on pair production of scalars that decay to electron or muon with missing energy are at 250 GeV with full run 2 data~\cite{Aad:2019qnd}.

\subsection{Strongly-interacting LLP}
Just as EW-charged LLPs manifest in the form of track-based signatures, strongly charged LLPs result in jets originating from a secondary vertex.  A possible example is provided by the minimal supersymmetric standard model, when small R-parity violating couplings \cite{Barbier:2004ez} open up decay channels of coloured R-odd particles.  We shall focus in particular on the $LQD$ coupling $\lambda'_{223}$ where the subscripts refer to an interaction between the second lepton doublet superfield, the second left-handed quark doublet superfield and the right-handed bottom superfield.  This results in various possible decay channels like $\tilde b_1 \rightarrow \mu c$ or $\tilde b_1 \rightarrow \nu s$, where $\tilde b_1$ is the lightest  bottom squark, which we assume to (nearly) coincide with the long-lived right-handed sbottom.

To test this scenario, we employ a simplified supersymmetric model where all the new-physics particles take a mass at the $10$~TeV scale, beyond the discovery reach of the LHC, with the exception of the lightest bottom squark of right-handed type $\tilde{b}_1$, whose mass is scanned over in the electroweak-TeV range.  The lifetime of $\tilde{b}_1$ is then entirely determined by the size of the coupling $\lambda'_{223}$ and can be varied freely. We generate samples of $10^6$ events for $pp \rightarrow \tilde{b}_1 \tilde{b}^*_1$ using Pythia8.  The production cross section is normalised to  the NNLO$_{\text{approx}}$ + NNLL values provided by the LHC cross-section Working Group \cite{SQprod}.

Several LLP searches are potentially sensitive to this scenario. The decay channel muon + jet is an obvious target for the DV+muon analysis, while the decay into neutrino + jet enters the scope of the DV+MET search. Finally, the bottom squark could be detected as a heavy long-lived charged particle. However, we do not compute limits from the corresponding search due to known uncertainties in the hadronization of such long-lived strongly charged scalars.
In addition, the disappearing track search is insensitive here because the displaced jets produced from the sbottom decays are hard, as a general rule.

  The limits obtained with {\tt CheckMATE} are presented in Fig.~\ref{fig:strongLLP_bounds} in the plane corresponding to the lifetime and mass of the bottom squark. We observe that, in this configuration where muon and MET productions are set equal, the DV+muon search proves slightly more competitive than the DV+MET analysis (which is consistent with the respective limits placed on the cross-sections by these searches in their respective benchmark scenario). The most constraining limits are placed on lifetimes $\sim1$~ns and exclude sbottom masses up to $\sim1.6$~TeV. This is more competitive than limits from the R-parity conserving scenario with $\tilde b \rightarrow b \tilde \chi$, which places a limit of 1270 GeV for a massless $\tilde \chi$~\cite{Aad:2021jmg}.
  %

  %
  Two types of prompt searches would also constrain the same model in a different parameter space searches for promptly-decaying leptoquarks, or usual sbottom limits (in the presence of some light LSP with the decay $\tilde b \rightarrow c \tilde \chi^+ \rightarrow c \bar \ell \nu \tilde \chi^0 $).  There have not been direct searches for sbottoms using this topology.  Limits on squarks in the $2\ell + 2\mathrm{jets} + \mathrm{MET}$ or in $ 2\mathrm{jets} + \mathrm{MET}$ are both expected to be much smaller than the standard topology ($\tilde b \rightarrow b \tilde \chi^0$). Limits on standard sbottom decay with full run 2 data are currently between 600-1270~GeV~\cite{Aad:2021jmg} depending on the mass of the final invisible particle. Although there have been searches for third generation leptoquarks, they focus on decays into tops~\cite{Aad:2021rrh, Sirunyan:2020zbk} (requiring e.g. b-tagged jets) and therefore do not apply directly to our model.

\begin{figure}[t]
    \centering
    \includegraphics[scale=0.5]{strong_excl.pdf}
    \includegraphics[scale=0.5]{strong_excl_test.pdf}
    \caption{Limits in the mass-lifetime (left) and mass-$\lambda$ $(=\lambda^{'}_{223})$ plane of the LLP sbottoms.
    \label{fig:strongLLP_bounds}}
\end{figure}
\FloatBarrier

\section{Comments and Outlook}

In this paper, we present the implementation of a new class of analyses in the {\tt CheckMATE} package, dedicated to long-lived particle searches, which can be used to reinterpret experimental limits on new physics models. We demonstrate the interplay of these searches in detecting strongly or weakly charged LLPs and obtain limits comparable to prompt limits in certain ranges of lifetime. In addition, these searches provide a way to probe couplings that can even be much smaller than those currently observable via low-energy intensity-frontier experiments, e.g.\ by measuring meson decays.

For our two scalar models, we find that the electroweak model can be constrained up to a mass of 480 GeV with the CMS dilepton search for lifetimes under 10 cm.  The charged track search is able to set bounds up to about 400 GeV for large lifetimes, greater than 10m.

We find an important gap in the search for electroweakly charged LLPs decaying to leptons in the intermediate lifetime range.  Since the disappearing track search employs a lepton veto, it is not possible currently to understand whether  models with decays into leptons can be observed (as the leptons emerge with potentially very large impact parameters and not from the primary vertex).  It would be fruitful to both improve the experimental search criteria as well as to provide efficiency for lepton identification based on $p_T$ and $d_0$ to remedy this situation.

For the strongly charged LLP, our limits are stronger than typical SUSY searches for particles with the same quantum numbers because of the much smaller backgrounds in these exotic searches.  An important gap here is that interaction of these particles with the detector material results in the efficiencies being dependent on some more unknown parameters than just the mass and lifetime. This is highlighted in Appendix C.  Also, this same issue prevents us from naively applying the charged track search to strongly interacting LLPs where charge exchange with the detector material becomes important.

Due to the absence of standardized detector objects in such searches, each re-interpretation heavily depends on the parametrized efficiencies published by the experiments. The validation of the five searches considered here has shown that this method gives relatively good results, which should allow the recast of the experimental limits for a very diverse range of models sharing an LLP as a common feature. However, since the efficiencies do rely on identifying the right truth-level particle via a user input and PDG code (possibly assigned ad hoc by the user for new particles), we do advise user vigilance when using these results.

Users can implement their own versions of an LLP analysis using the AnalysisManager, which simplifies the setting up of detector parameters, stores the expected and observed events reported by the experiment in the correct format for future use in statistical calculations and provides a skeleton C++ analysis code with access to all detector objects.  Currently, searches based on ionisation, like those for monopoles or multi-charged objects, are not possible to implement because the relevant efficiencies are not publicly available.  However, if the efficiency tables based on mass, charge and momentum were to be provided by the experiments in the future, it would be possible to implement them as well, without further difficulties.  The latest version of {\tt CheckMATE} can always be downloaded from
\url{https://github.com/CheckMATE2}.  Documentation and validation notes can be found at  \hbox{\url{https://checkmate.hepforge.org/}}.

\acknowledgments

We thank Giovanna Cottin and Andre Lessa for useful discussions.
R.~RdA acknowledges partial funding/support from the Elusives ITN (Marie Sk\l{}odowska-Curie grant agreement No 674896), the ``SOM Sabor y origen de la Materia" (FPA 2017-85985-P).
Z.~S.~W. is supported partly by the Ministry of Science and Technology (MoST) of Taiwan with grant number MoST-109-2811-M-007-509, and partly by the Ministry of Science, ICT \& Future Planning of Korea, the Pohang City Government, and the Gyeongsangbuk-do Provincial Government through the Young Scientist Training Asia Pacific Economic Cooperation program of the Asia Pacific Center for Theoretical Physics. Z.~S.~W. also thanks the Witwatersrand University for hospitality where part of this work was initiated. KR is supported by the National Science Centre, Poland under grants: 2016/23/G/ST2/04301, 2015/18/M/ST2/00518, 2018/31/B/ST2/02283, 2019/35/B/ST2/02008, and the Norwegian Financial Mechanism for years 2014-2021,
grant DEC-2019/34/H/ST2/00707. F.D.~acknowledges support of the
BMBF Verbundprojekt 05H2018 and the DFG grant SFB CRC-110 {\em Symmetries and the Emergence of Structure in QCD}.

High Energy Physics - Phenomenology (hep-ph)


Transport of magnetic flux from the canopy to the internetwork

Abstract: Recent observations have revealed that 8% of linear polarization patches in the internetwork quiet Sun are fully embedded in downflows. These are not easily explained with the typical scenarios for the source of internetwork fields which rely on flux emergence from below. We explore using radiative MHD simulations a scenario where magnetic flux is transported from the magnetic canopy overlying the internetwork into the photosphere by means of downward plumes associated with convective overshoot. We find that if a canopy-like magnetic field is present in the simulation, the transport of flux from the canopy is an important process for seeding the photospheric layers of the internetwork with magnetic field. We propose that this mechanism is relevant for the Sun as well, and it could naturally explain the observed internetwork linear polarization patches entirely embedded in downflows.

\section{Introduction}

 High spatial resolution and high signal-to-noise spectropolarimetric
 observations, such as those from the spectropolarimeter (SP, \citealt*{HinSP2001}) on-board the
 Hinode satellite \citep{Hinode,Tsuneta+others2008,Suematsu+others2008,Shimizu+others2008} and the full
 Stokes magnetograph IMaX \citep{MartinezPillet+others2010} on-board
 the balloon-borne telescope Sunrise \citep{Barthol2010,
   Solanki+others2010,Berkefeld+others2010, Gandorfer+others2010}, have made it possible to characterize in
 greater detail the photospheric magnetic fields in the internetwork
 (INW) regions of the Sun. For a recent review on small-scale magnetic
 fields see \citet{deWijn+others2009}.

It is generally proposed that the INW fields come from below (e.g., \citealt{WedemeyerBohm+others2009}),
emerging from the high plasma-$\beta$ ($8 \pi \rho / B^{2}$) region where the magnetic
fields can be inductively amplified, into the
higher layers where the plasma-$\beta$ is low.
Possibilities include a source associated with the global dynamo (recycling of flux from decaying
active regions and inductive amplification, cf. \citealt{Ploner+others2001,Hagenaar+others2003,TrujilloBueno+others2004}) or local dynamo action \citep{Petrovay+Szakaly1993}.

 What all these mechanisms above have in common is that the flux
 originates mainly from below the surface. Therefore, one would expect
 the appearance of INW magnetic fields to be associated with flux
 emergence. This has been observed in several instances, e.g.,
 \citet{Centeno+others2007,MartinezGonzalez+others2007,MartinezGonzalez+BellotRubio2009}. Recent
 observations from the IMaX instrument on-board the Sunrise
 balloon-borne telescope allowed \cite{Danilovic+others2010} to study
 in detail linear polarization patches in the INW. IMaX observed the
 full Stokes vector of the photospheric 525.0 nm \ion{Fe}{1} line at
 high signal-to-noise ratio and high spatial resolution. Of the nearly
 5000 identified linear polarization patches 8 \% were entirely
 embedded in downflows over their whole lifetime. This rules out
 emergence from below the solar surface as their origin. These
 observations together with magnetoconvection simulations prompted us
 to explore a new mechanism for transporting flux into the INW
 photosphere. This mechanism would produce naturally linear
 polarization patches entirely embedded in downflows. In the proposed
 mechanism magnetic field is dragged down from a magnetic canopy
 \citep{Gabriel1976} by downflows due to convective overshoot.

\section{MHD simulation and studying the topology}

 The magnetoconvection simulation used in this paper is the same one
 as in \citet{Pietarila+others2010}, but run for a more extended
 period of solar time. We use the MuRAM code (MPS/University of
 Chicago Radiative MHD, \citealt{Voegler+others2005}) which takes into
 account effects of full compressibility, open boundary conditions as
 well as non-gray radiative transfer and partial ionization. The
 simulation domain is a 24Mm$\times$24Mm$\times$1.68Mm
 (576$\times$576$\times$120 grid points) box of which approximately
 700 km is above the $log(\tau_{500nm})=0$ level. The boundary
 conditions are periodic in the horizontal direction. The upper
 boundary condition for the magnetic field is potential (vertical at
 infinity) and the velocity obeys an open boundary condition both at
 the top and bottom of the domain. The simulation was started from a
 snapshot of a non-magnetic convection run that had reached a statistical steady state, on to which we added a
 magnetic field. The initial configuration of the field was a strip of
 magnetic field with a squared Gaussian profile located in the center
 of the domain in the x-direction and aligned along the y-axis. We
 chose a peak value of 300 G for the strip. This way the initial
 magnetic field is energetically important above the surface but much
 less so below the surface. The full width at half maximum was chosen
 to be 3.157 Mm, so that the patch of magnetic field spans several
 granules. This setup is reminiscent of an isolated, unipolar network
 lane. The simulation was run for $\approx$ 1 hour of solar
 time. Since we are interested in the region outside the initial strip
 of magnetic field, we take advantage of the periodic boundary
 conditions in the horizontal directions and place the strip at the
 edges of the domain in the figures instead of having it in the center
 as in \citet{Pietarila+others2010} (i.e., from now on x=0 and x=24 Mm
 refer to the center of the initial strip).

\section{Results}

 In the beginning of the simulation (left panels in
 Fig.~\ref{fig:bigsnap}) the vertical magnetic field is confined to
 the initial magnetic field strip. The potential field upper boundary
 condition allows the magnetic field to expand near the top of the
 domain and a canopy-like structure composed of nearly horizontal
 field forms rather quickly over the initially field-free
 regions. (From here on the term ``canopy'' is used for this
 structure.) Since the initial field is unipolar and we have periodic
 boundary conditions in the horizontal directions, the canopy field
 expands in a ``wine glass'' shape: the field is vertical at $x=12$ Mm
 where the field from two same-polarity network lanes meets. The
 spreading out of the field is apparent in the magnetic field azimuth:
 it is predominantly positive in the left ($x=0-12$ Mm) and negative
 in the right ($x=12-24$ Mm) halves of the box. The effect of
 granulation on the spatial distribution of the magnetic field is
 visible early on: the field is gathered in the intergranular downflow
 lanes. As the simulation progresses the magnetic field fills up the
 entire domain (right panels of Fig.~\ref{fig:bigsnap}). At all
 heights the granulation pattern is visible in the spatial
 distribution of the magnetic field: field is found predominantly in
 the downflow lanes. The plasma-$\beta$ is well above unity ($\ge$ 10)
 in most of the simulation domain, only close to the very top does it
 decrease to below unity. The field azimuth has the same spatial
 asymmetry (goes from positive to negative with {\it x}) at all
 heights.

\begin{figure*
\includegraphics[width=14cm]{Fig1.eps}
\caption{Magnetic field and vertical velocity close to the beginning
  ($t=1.5$ min, panels on the left) and end ($t=45$ min, panels on the
  right) of the simulation. Sub-panels in top panel from top to bottom
  show x-y slices of the vertical magnetic field strength at 3
  different heights (400, 200 and 0 km) and a x-z slice of the
  simulation domain. $z=0$ km is defined as where the average
  temperature (over the entire simulation domain) at $t=0$ is
  $\approx$ 6000 K. The dashed lines in the bottom sub-panel mark the heights of the x-y
  slices. Initial distribution of field corresponds to unipolar (positive polarity) field at the left and right boundaries of the computational domain with field-free area between them. The two panels below show the same for the horizontal field
  and the vertical velocity, respectively.}
\label{fig:bigsnap}
\end{figure*}

 The magnetic field is transported to the initially nearly field-free
 regions of the photosphere by two mechanisms, namely by horizontal
 diffusion (random walk) from the initial strip and by a process where
 the downflows due to the convective overshoot pull the field down
 from the canopy. With the simulation resolution, $\simeq$41 km in the
 horizontal direction, we do not achieve magnetic Reynolds numbers
 high enough to have local dynamo action. Therefore, we can exclude it
 as a source for the field in the initially nearly field-free
 regions. This does not exclude stretching and inductive amplification
 of the unsigned flux as a source in the later stages of the
 simulation. As we will later show when discussing the time evolution
 of the magnetic field topology this is an important process in
 supplying field to photospheric regions outside the initial strip. The prevalence
 of the spatial azimuth asymmetry throughout the simulation, both in
 height and time, points to a strong connection with the canopy which has the same asymmetry. This
 connection is also seen in the time evolution of the magnetic field
 (Fig.~\ref{fig:absb}). In the initially field-free regions, the
 magnetic field appears at a given height nearly simultaneously at all
 locations in the horizontal plane. The fine structure of the field
 has a spatial scale of about 1 Mm, corresponding to the granulation
 scale. The field fills the volume starting from the top of the domain
 and reaches the surface after $\approx 10$ minutes, i.e, after a few
 granulation timescales. The horizontal diffusion of the field from
 the initial strip is seen as the slight increase with time of the width
 of the strong field region at $x \approx 0$  and $x \approx 24$. Its
 effect in filling the volume with field at distances greater than a
 couple of Mm from the strip is negligible compared to the flux
 transported down from the canopy. The nature of the downflows and canopy-like magnetic fields will be discussed in Section ~\ref{sec:disc}.

The dominance of downflows over horizontal diffusion in providing flux close to the surface is also seen in histograms of the magnetic field
 strength outside the initial strip taken at
 different times during the simulation (two first panels in Fig.~\ref{fig:fields}). They show the strengthening of the horizontal and vertical fields. The vertical field is not unipolar as
 one would expect if the main transportation mechanism was horizontal
 diffusion. In fact, the signed vertical flux is significantly smaller
 than the unsigned flux at all times (right panel in Fig.~\ref{fig:fields}). There are indications of
 saturation of the downflow mechanism after $\approx$ 600 seconds when
 the signed flux increase becomes faster than the increase of the
 unsigned flux.

\begin{figure*}
\includegraphics[width=14cm]{Fig2.eps}
\caption{Time evolution of the mean (averaged along y-direction)
  vertical (left) and horizontal (right) magnetic field at z=400, 200
  and 0 km during the first 18 minutes of the simulation. }
\label{fig:absb}
\end{figure*}

\begin{figure*}
\includegraphics[width=14cm]{Fig3.eps}
\caption{Histograms of the horizontal (left) and vertical (middle)
  magnetic fields at $z=220$ km for each time step of the simulation
  (blacks correspond to the beginning and reds to the end of the
  simulation) in a strip at $x=4-10$ Mm of the
  simulation domain (i.e., outside the initial strip of magnetic
  field). Right plot shows the time evolution of the signed (in black)
  and unsigned (in red) vertical magnetic flux in the strip. }
\label{fig:fields}
\end{figure*}

 To quantify the apparent speed at which the horizontal magnetic field is
 advected down from the canopy we look at when the average (in
 y-direction) total field strength at a given height becomes larger
 than 1 G or the unsigned vertical field exceeds 0.33 G or the
 horizontal field 0.66 G (Fig.~\ref{fig:canvel}). The limit was chosen based
 on images similar to Fig.~\ref{fig:absb} where 1 G outlines well the
 downward moving canopy. The thus computed average downward  speed of the canopy's base is found to be $\approx$1 km
 s$^{-1}$. The simulation domain
 extends $\approx$ 700 km above $\log(\tau) \approx 0$ in height so an
 average speed of 1 km s$^{-1}$ is consistent with field from the canopy reaching
 the solar surface in $\approx$ 10 minutes.

\begin{figure*}
\includegraphics[width=14cm]{Fig4.eps}
\caption{Transport of magnetic fields by downflows. Left and middle
  panels show the time when the base of the canopy (vertical and horizontal components) reaches a given height. The
  canopy base is considered to have reached a height when at that height
  the average (over a 400 km wide strip) vertical field strength is
  above 0.33 G or  the horizontal field exceeds 0.66 G. The different
  colors are for different positions along the x-axis (shown
  in the color bar). The right panel shows the average downward velocity of the
   canopy base as a function of position along the $x$-axis. The velocity is derived from a linear fit between z=$0$ and $500$ km to the curves in the two other panels. Dotted line is for horizontal, dashed for vertical and dash-dotted for total field.  }
\label{fig:canvel}
\end{figure*}

 If the magnetic field is dragged down by downflows one would expect a
 prevalence of U-loops in the field reaching the surface,
 especially in the beginning of the simulation before the field
 topology becomes more complex as convection moves the dragged-down
  field in different directions. A 3-dimensional rendering of the
 magnetic field confirms this as is shown in Fig.~\ref{fig:loops}. Color
 coding of the field lines represents vertical velocity (with red signifying downflows).
 It shows clearly that the U-loop is associated
 with a downflow.

\begin{figure*}
\includegraphics[width=14cm]{Fig5.eps}
\caption{Left: Volume rendering of magnetic field lines in a subvolume
  of the simulation domain ($x=6-10$ Mm, $y=2-6$ Mm and $z=-0.5-0.7$ Mm). The
  snapshot is at $t=14.7$ min. Right: zoom into the region outlined by the
  light blue box. The surface is the ${\rm B_{x}}$ component of the
  magnetic field at $z=0$ km. Color of the field lines in the right panel
  shows the vertical velocity (red is a downflow, blue an
  upflow. Color scale saturates at $\pm$1 km s$^{-1}$). The volume rendering was made using the Visualization and Analysis Platform for Ocean, Atmosphere, and Solar Researchers (VAPOR, www.vapor.ucar.edu)}
\label{fig:loops}
\end{figure*}

 Since we are interested in the interaction of the canopy and the
 convective motions, we also performed an analysis to characterize the
 flux in the lower to mid photospheric layers as belonging to U-loops,
 $\Omega$-loops, or field lines going from below the surface to above
 the higher photospheric layers. For this purpose we define lower and
 upper surfaces at $z=-294$ and $z=406$ km. For each grid point
 $(x,y,z)$ within this domain the field through the point is traced in
 both directions until it leaves the subdomain.  We label all the
 points along the field line as U-type if both ends of the field line
 exit the subdomain through the upper plane; the points are labeled as
 $\Omega$-type if both ends of the field line exit through the bottom
 plane; the field lines whose one end of which passes through the
 bottom plane and other end through the top plane are defined as being
 either positive or negative polarity depending on the orientation of
 the field.

The resulting classification of the magnetic field at z=$224$~km is
shown for four timesteps in Fig.~\ref{fig:loopcount}. The first
snapshot is at 4.2 minutes, i.e., when the downflowing canopy
reaches this height for the first time. At this time the downflows
have had sufficient time to pull parts of the horizontal canopy down
into the lower photosphere. The INW flux is located in the junctions
of the intergranular lanes and it is nearly entirely in the form of
U-loops or positive polarity field passing through the domain. The
positive polarity field is due to flux expulsion: as 4.2 minutes is
comparable to the lifetime of the granulation there has been enough
time for flux expulsion to concentrate the very weak field in the
region from $x=5$ to $20$ Mm into a number of isolated elements which
thread the subdomain from bottom to top (and hence appear black in the
topology maps). At later times, 14.1, 27.3 min, and 45 min, there are
still U-loops present in the INW, but also a large amount of
$\Omega$-loops are seen as well as positive polarity field.  The
U-loops appear as small bipoles in the bottom sub-panels showing the
vertical field, with a vertical magnetic field strength which is
generally stronger than in the $\Omega$-loops and in the vertical INW
field. Because we do not have a supergranular scale flow in the
simulation, flux which was originally in the unipolar network features
begins to get concentrated into granular downflow lanes while at the
same time gradually diffusing outward from its original location, as
can be seen in the snapshot from 45 minutes.

Figure~\ref{fig:time_evol} displays the time evolution of the average
unsigned flux density in the loops and single polarity field.
{\black{Throughout the simulation most of the increase in unsigned flux
    in the volume within the height range $-294$~km~$\le z \le 406$~km
    is from field lines which enter the volume from $z=-294$~km and
    leave through $z=406$~km. In between these two heights the
    unsigned field lines may be rather tangled and form small loops
    which contribute to the net unsigned flux. Since the field lines
    do not again leave the subdomain through the bottom they are not
    classified as loops. Unsigned flux of this type can increase
    through interaction with the near-surface turbulent eddies, and
    represents the pulling of the lowest lying flux above the surface
    into the convection zone. Higher lying flux is then pulled down,
    forming U-loops both ends of which pass through the plane
    $z=406$~km.} \black} During the first 10 minutes of the simulation
    U-loops dominate over $\Omega$-loops.
They remain an important factor
also in the later stages of the simulation.
As the field in the originally non-magnetic region builds up, field lines
with different topologies can interact and reconnect. The reconnection occurs
preferentially in the higher levels where the lower pressures allow the
flux tubes to expand horizontally and, therefore, more easily come into contact
with other flux tubes. For this reason the reconnection tends to decrease the depth
to which U-loops penetrate. Furthermore, magnetic tension
can more easily pull the shallow, reconnected, U-loops out of the simulation domain. This may explain the decrease
of flux associated with U-loops in Figure~\ref{fig:time_evol}. Small-scale inductive
stretching of the INW field lines over several turn-over times of the
granulation {\black can pull the now tangled field lines below the
  $z=-294$~km plane explaining \black} the appearance of the
$\Omega$-loops (green in Figs.~\ref{fig:loopcount} and
~\ref{fig:time_evol}).  {\black The average unsigned flux densities shown
  in Figure~\ref{fig:time_evol} are expected to depend strongly on the
  amount of signed flux in the network lane. The values are also
  likely affected by the choice of a potential field boundary
  condition and the height at which this boundary is placed.
  Furthermore, the amount of tangling within the photosphere in the
  simulations depends on the magnetic Reynolds number, with
  small-scale dynamo action occurring for sufficiently high values
  \citep{Vogler+Schussler2007,PietarilaGraham+others2010}. For these
  reasons, while it is qualitatively clear that the canopy can supply
  a significant amount of flux to the photosphere, it is difficult to
  estimate the amount quantitatively from the present simulations. \black}

\begin{figure*}
\includegraphics[width=7cm]{Fig6a.eps}\includegraphics[width=7cm]{Fig6b.eps}

\includegraphics[width=7cm]{Fig6c.eps}\includegraphics[width=7cm]{Fig6d.eps}
\caption{Snapshots from the simulation at a height of 224~km at times
  t=4.2 (upper left), 14.1 (upper right), 27.3 (lower left), and 45
  minutes (lower right).  For each time the upper part of the panel
  shows the vertical component of the velocity saturated at $-2$ km
  s$^{-1}$ (blue) and $2$ km s$^{-1}$ (red).  The middle parts are
  maps of the distribution of the different types of loops.  U-type
  loops (field lines both ends of which leave through the top of the
  subdomain $-294$ km $\le z \le 406$ km are shown in red),
  $\Omega$-loops (both ends of which leave through the bottom of the
  subdomain) in green, field lines passing from the bottom to the top
  of the subdomain in black, and in the opposite direction in
  grey. The types of loops are only shown where the field exceeds
  $1$~G, all other regions are colored white.  The bottom third of
  each snapshot shows the vertical component of the magnetic field,
  saturated at $\pm$ 50 G.  }
\label{fig:loopcount}
\end{figure*}

\begin{figure*}
\includegraphics[width=7cm]{Fig7a.eps}\includegraphics[width=7cm]{Fig7b.eps}

\caption{ The time evolution of the different types of spatially averaged
  flux densities at height
  224~km. The classification into U, $\Omega$ or straight-through
  topology is based on the subdomain $-294 \le z\le 406$~km.  The
  amount of \black unsigned \black flux associated with U-type loops is shown in red,
  $\Omega$-loops in green, field lines passing from the bottom to the
  top of the subdomain in black, and in the opposite direction in
  grey.  The total amount of signed flux passing through the box has
  been subtracted from the black curve, \black so that the curve gives the increase in
      unsigned flux due to turbulent stretching and tangling of the
      field within the subdomain\black. The right panel shows a zoom
  into the lower left corner of the left panel.  }
\label{fig:time_evol}
\end{figure*}

 Figure ~\ref{fig:snaps} shows the time evolution of the horizontal
 and vertical magnetic fields and the vertical velocity of the U-loops
 shown in Fig.~\ref{fig:loops}. In the beginning there is only a
 little magnetic field in the junction of the intergranular lanes, but
 as the downflow plume reaches this height the field becomes
 stronger. The area with magnetic field becomes larger with time,
 consistent with the U-loop becoming more submerged. The U-loop is
 stretched and distorted as the field is further dragged into the
 downflow lanes. This is clearly seen in the fanning out of the
 horizontal field along the downflow lanes. After some time as the
 granulation evolves the field occupies a large fraction of the
 downflow lanes and is no longer as strong or as localized. The
 horizontal field strength in the feature is over 80 G at its
 strongest. The vertical field strengths are slightly lower, around 60
 G. The size of the feature before it is spread out and weakened is
 less than an arcsecond across. This is a typical example of a submerging U-loop.

\begin{figure*}
\includegraphics[width=14cm]{Fig8.eps}
\caption{Time evolution of magnetic field and velocity at $z=220$ km in
  a region with magnetic flux appearing in the form of a U-loop
  from above. Top row: Horizontal magnetic field strength, middle row:
  vertical magnetic field strength, bottom row: vertical
  velocity (red corresponds to downflows and blue to upflows). Times for snapshots are given in the bottom of the top row
  images. The region is centered at $x=13.5$ Mm and $y=3.7$ Mm. }
\label{fig:snaps}
\end{figure*}

\section{Discussion}
\label{sec:disc}
 In this paper we have looked at a scenario for transporting magnetic
 flux to the quiet Sun INW surface: already existing magnetic field is
 dragged down from a canopy-like configuration (consisting of mostly
 horizontal magnetic field) by downflows related to convective
 overshoot. An observational signature of the mechanism would be
 patches of linear polarization entirely embedded in downflows. Such
 patches were found to be common (8\% of all linear polarization
 patches identified in the study) in the quiet-Sun INW by
 \cite{Danilovic+others2010}.

\nocite{Abbett2007}

A similar process has been reported by Abbett (2007, in particular
Fig.10 and the text describing it). The current study focuses on the canopy fields from
network elements whereas Abbett studied convective dynamo action in the
presence of an initially uniform field. While the rates and
importance of the processes will differ in the two cases, the occurrence
of the process in the two simulations suggests that it is
robust. Importantly we note that Abbett (2007) simulations extended to
coronal heights with a simpler treatment of the radiative losses. It
should also be noted that the proposed mechanism is closely related to
that of downward pumping as described by
\cite{Thomas+others2002,Brummell+others2008}, however, now applied to
the INW surrounded by network magnetic field and overlaid by a canopy.

In the simulation nearly all of the horizontal field is in a high plasma-$\beta$ regime. If the field was stronger, it would appear more rigid to the downflows and, thus, would also less likely be dragged down. The horizontal field near the top of the
simulation domain (our canopy) should be thought of as weaker flux lying beneath the actual chromospheric canopy as inferred from observations. This weaker canopy is there because the chromospheric canopy's lower boundary is not sharp
but is rather disturbed by shocks, overshooting convection, propagating waves and turbulence.

Several aspects of our simulation are idealized, especially the potential field boundary condition imposed at a relatively
low height. The mechanism we have investigated, however, is likely to take place in the real Sun as
well if there is field associated with the overlying canopy at heights reached by the overshooting convection. The field dragged down below the surface by the downward plumes could then act as an additional source for the
generation of more field.

\cite{Danilovic+others2010} present statistics of linear polarization
patches in an INW region observed with the IMaX instrument on
Sunrise. They find that 8\% of the observed linear polarization
patches are associated purely with signatures of
submergence. Moreover, 16\% of all the observed patches begin in a
downflow. These features are not consistent with emergence of flux
from below, where an upflow would initially be expected. However, they
are compatible with the scenario proposed here. Furthermore, the sizes
and field strengths of the observed linear polarization patches are
comparable with the features seen in the simulation. Of course, we
cannot rule out other possible explanations for the linear
polarization features embedded in downflows.\black For example, the
local dynamo simulations suggest the existence of very small-scale
loops with footpoints embedded in intergranular lanes
\citep{Vogler+Schussler2007,Schussler+Vogler2008}. Their polarimetric
signature can vary as they are redistributed by the turbulent flows
\citep{Danilovic+others2010b}. In addition, $\Omega$-loops, emerged or
generated by reconnection above the solar surface, can be pulled back
down by convective motions \citep{Chae+others2004,Stein+Nordlund2006} producing the
linear polarization features associated with downflows as a result.\black

It is difficult for a number of reasons to conclusively identify the mechanism proposed here in observations of the real Sun. We expect a
 mixture of processes to transport magnetic field to the INW and only a small portion of the observed linear
 polarization patches are embedded in downflows. In the simulation the cleanest examples are seen in the
 beginning when there is only little flux outside the initial strip. As the simulation evolves, the magnetic
 field configuration becomes more complex. One possible observational
 test would be to look for spatial patterns in the azimuth of INW
 linear polarization patches embedded in downflows. This
 would, however, require solving the 180 degree ambiguity correctly.

\acknowledgements{This work was partly supported by the WCU grant
  No. R31-10016 from the Korean Ministry of Education, Science and
  Technology.}

\bibliographystyle{aa}

Solar and Stellar Astrophysics (astro-ph.SR)


Qualitative Framing of Financial Incentives - A Case of Emotion Annotation

Abstract: Online labor platforms, such as the Amazon Mechanical Turk, provide an effective framework for eliciting responses to judgment tasks. Previous work has shown that workers respond best to financial incentives, especially to extra bonuses. However, most of the tested incentives involve describing the bonus conditions in formulas instead of plain English. We believe that different incentives given in English (or in qualitative framing) will result in differences in workers' performance, especially when task difficulties vary. In this paper, we report the preliminary results of a crowdsourcing experiment comparing workers' performance using only qualitative framings of financial incentives. Our results demonstrate a significant increase in workers' performance using a specific well-formulated qualitative framing inspired by the Peer Truth Serum. This positive effect is observed only when the difficulty of the task is high, while when the task is easy there is no difference of which incentives to use.

\section{Introduction}
In online labor markets, the quality of the answers returned by inexpert, remote workers presents an important concern for the usefulness of such crowdsourcing process. The~standard techniques for quality assurance, such as combining answers from different workers and posterior removal of non-satisfactory answers \cite{ipeirotis2010quality}, assume the existence of a correct answer.
Yet, for judgment tasks---tasks involving human subjective judgment such as collecting opinions, preferences, relevance estimations, ratings, or emotion categorization---it is difficult to establish the validity of a specific answer based only on the~(dis)agreement with the majority label.

Previous research showed the potential of motivating workers by introducing additional incentives, such as bonuses for good-quality answers \cite{ho2015incentivizing,harris2011you}. Yet, these works require automatic estimation of the answers' correctness.
The schemas without this requirement would be preferable for judgment tasks.
In this direction, other works suggest that giving bonus based on the consistency with peer answers is most advantageous \cite{Shaw:2011:DII:1958824.1958865,huang2013enhancing,faltings2014d}.
Such peer-oriented schemas originate from a game-theoretic approach to incentivize truthful answers \cite{prelec2004}, and thus they imply the mathematical formulation of bonus computation. However, a lay person from Amazon Mechanical Turk (MTurk) is unlikely to fully understand the implications of such computations, especially in case of non-mathematical judgment tasks (such as emotion labeling).
To overcome this, we suggest to study alternative simplified formulations of the peer-oriented schemas suitable for application in judgment tasks without any mathematical formulas. We refer to this approach as qualitative framing.
We report in this paper the results of a crowdsourcing experiment conducted on the MTurk, where we compare several qualitatively framed peer-oriented incentives with the ones inspired by previous research on social incentives in crowdsourcing \cite{Shaw:2011:DII:1958824.1958865} across two difficulty levels in the context of emotion annotation.

\section{Experiment Design}

We use in our experiment a specific judgment task---annotating emotions in tweets.
We ask each worker to label $10$ tweets after completing a mandatory tutorial. For each tweet, the worker should  indicate the dominant emotion felt by the author while writing the tweet by selecting one of $20$ emotion categories, \textit{No~emotion}, or \textit{Other~emotion} from the Geneva Emotion Wheel (GEW, version~2.0) \cite{scherer2013grid}.
 For example, for the tweet \textit{``Woooo! It's a good day"}, the worker is likely to select the emotion \textit{Happiness}.
The worker is also asked to provide the excerpts of the tweet indicating the presence of the chosen emotion (\textit{``Woooo"} and \textit{``good day"} in the example) as well as additional emotion indicators (not from the tweet text). The same task was presented and studied in our previous work \cite{sintsova-musat-pu:2013:WASSA}. We also included optional demographic and feedback questionnaires.
For the sake of brevity, we report here only the results for the emotion labels.

We randomly assigned each worker to an incentive treatment condition employing a specific bonus formulation.
To discover whether the effects of incentives differ depending on the difficulty of the annotation, we assigned workers at random to one of two difficulty levels.

The base payment was set to $\$0.5$ USD, and there was information about a chance to obtain an additional bonus of $\$0.1$ USD. A specific description of the bonus (dependent on the assigned incentive condition) was presented to a worker right before starting the tweets' annotation:

\noindent\textit{``You can qualify for the additional $\$0.1$ bonus if your answers [...]" }

\vspace{1pt}\noindent
\textbf{Normative: }
\textit{``... demonstrate an additional effort." }

\vspace{1pt}\noindent
\textbf{Experts:}
\textit{``... are extremely accurate according to our experts." }

\vspace{1pt}\noindent
\textbf{Professors: }
\textit{``... are approved by our professors."}

\vspace{1pt}\noindent
\textbf{Peer Agreement: }
\textit{``... agree with those of other workers." }

\vspace{1pt}\noindent
\textbf{Peer Truth Serum 1 (PTS1): }
\textit{``... are more surprisingly common with other workers than collectively predicted."} \\
\textbf{Peer Truth Serum 2 (PTS2): }
\textit{``... both agree with those of other workers and at the same time novel."}

\vspace{4pt}\noindent
The second variable was the difficulty of the tweets to label:

\vspace{2pt}\noindent
\textbf{Easy dataset} The first, ``easy'' dataset consisted of $70$ tweets manually chosen such that they were explicitly emotional and not difficult to interpret, judged by us as requiring less effort to annotate.
The example tweet is \textit{``You said it would be different but like usual nothing has changed"}.

\vspace{2pt}\noindent
\textbf{Difficult dataset} The second, ``difficult'' dataset comported $70$ emotional tweets that were less obvious to interpret and believed to require more attention to annotate, because they contained negated emotional terms, or were of a sarcastic nature, or simply confusing. The example tweets are \textit{``I do not regret a single thing being born as Dusun."} and  \textit{``Best part about rush hour is driving into it going Chicago!!! \#not"}.

\section{Results}
We ran the experiment on MTurk between May 13th and June 3rd  2015, with $1,135$ workers completing the full task.
The random assignment of workers to treatment conditions resulted in a fairly even split of the non-spamming workers across all conditions ($94.6 \pm 3.2$ workers per group).

We estimate agreement and correctness of the returned emotion labels to evaluate the workers' performance. We report average per-worker scores.
The computation of each metric for one worker's answer on a tweet is as follows.
\textit{Category agreement} of a worker's emotion category label is computed as the percentage of agreed peer labels for the same tweet.
 To compute \textit{category correctness}, we first obtain the ground-truth categories by aggregating the workers' labels from all incentive conditions.
  For each tweet, we extract the majority label and all the categories with a relatively high assignment number (we use the threshold of $0.5$ minimum ratio from the majority label count). The worker's label is considered to be correct if it is within the set of extracted ground-truth labels for the tweet.

\begin{figure}[t]
\centering
\includegraphics[width=\linewidth]{figures/category_agreement_diff-cut.pdf}
\caption{The effect of incentives on average \textit{category agreement} of workers assigned to the difficult dataset. Error bars indicate a standard error.
}
\label{fig:diffCat}
\end{figure}

\vspace{2pt}
Because the two-way ANOVA revealed a significant interaction between the two studied factors,
we decided to resort on analyzing the simple effects of incentives at each level of difficulty.
We found \textit{no significant effect  of the incentives' framing on the quality of emotion labels on the easy dataset} (p-value of one-way ANOVA is $0.183$ for category agreement and $0.571$ for category correctness on the easy dataset).
Nevertheless, the incentives significantly affect worker's performance on the difficult dataset
Indeed, p-value of one-way ANOVA is~$<10^{-7}$ for category agreement and $0.033$ for category correctness.
Figure \ref{fig:diffCat} visualizes the differences  in category agreement between incentives' conditions on the difficult dataset.
Pairwise comparison of the incentives reveals the\textit{ advantage of PTS2 incentive on the difficult dataset}: it is significantly better than any other incentive in terms of category agreement (the highest p-value is $0.045$ when comparing with PTS1), and it leads to the highest category correctness while significantly outperforming the Normative condition (with p-value~$=0.023$).
The Normative condition, which is appealing only to workers' honesty, in its turn results in the lowest category agreement on the difficult dataset, significantly different from PTS1, PTS2, and Professors incentives (with p-values~$\le0.034$).

The workers' answers to the feedback questionnaire verify that the difficult dataset is indeed perceived as more difficult than the easy dataset and that it requires additional cognitive effort to label.
This implicitly approves our methodology for selecting tweets for the difficult and easy datasets.
The further statistical analysis revealed no effect of the incentives' framing on the workers' perceived task comprehension, enjoyment, easiness, or cognitive effort.

\section{Discussion}

The results of our experiment revealed the differences in the effects of the incentives' framing depending on the difficulty of the emotion annotation task.  When the task is easy, no incentive results in significantly improved performance compared to the other incentives. However, when difficulty of the task increases, our findings suggest that it is beneficial to use a particular incentive: with a well-formulated version of the  Peer Truth Serum (PTS) bonus, workers tend to output more correct and agreeing emotion labels.
 This finding is in accordance with benefits of this schema in its mathematical formulation for non-judgment tasks \cite{faltings2014d}. We study two different qualitative formulations of the Peer Truth Serum (\textit{``surprisingly common"} vs. \textit{``agreeing and novel"}) and found that only the second formulation can lead to significantly better answers  showing the importance of better qualitative framing of the incentives.

The fact that the qualitative framing of incentives affects the quality of answers at least for one fixed bonus amount encourages the future investigation, testing, and deployment of more advantageous qualitative bonus formulations.
Our  findings could be of potential value to the researchers and practitioners aiming to design incentives mechanisms for judgment tasks.


Detecting Low-Mass Supermassive Black Holes

Abstract: We demonstrate the feasibility of uncovering supermassive black holes in late-type, quiescent spiral galaxies by detecting signs of very low-level nuclear activity. We use a combination of x-ray selection and multi-wavelength follow-up. Here, we apply this technique to NGC 3184 and NGC 5457, both of type Scd, and show that strong arguments can be made that both host AGNs.

\section{Introduction}

 As discussed in Mathur et al.\ (these Proceedings) and in \cite{gmff08},
 the only way to detect low-mass supermassive black holes (SMBHs) may
 be by their accretion activity. We choose to use x-ray selection to
 identify candidate AGNs for the following reasons: First, x-ray
 emission is a hallmark of AGNs. Second, x-rays can penetrate
 obscuring material which may be hiding the line emitting
 regions. Third, there are fewer sources of x-rays in a galaxy than
 there are of optical and UV emission and so dilution of the AGN
 signature by host galaxy light is less of a problem, even if the AGN
 is moderately obscured.
Fourth, even if, as expected
in some theories \citep{n00,l03} AGNs that have
luminosities or accretion rates below a cut-off value do not have
broad-line regions, they should still be detectable in x-rays.
The disadvantage is that x-ray
observations by themselves cannot always distinguish between AGNs and
other x-ray sources, such as x-ray binaries (XRBs) and ultraluminous
x-ray sources (ULXs). Multi-wavelength data are needed to determine
the type of source.

As examples, we present here two late-type spiral galaxies, NGC 3184
and NGC 5457 (M101). Their nuclear optical spectra show no signs of
AGNs \cite{hfs97-3}. But using archival \textit{Chandra\/}\ data and
multi-wavelength data from the literature, we show that the sum total
of evidence strongly suggests that the galaxies host AGNs.

\section{NGC 3184}

This galaxy is of type Scd, classified as having an \textsc{Hii}
nucleus by \citet{hfs97-3}, and is at a distance of 8.7 Mpc
\citep{t88}. The argument for the presence of an AGN is based on the
following: (a) \emph{The nucleus is an x-ray source}. Observations
with \textit{Chandra\/}\ totaling 60 ks detect a point-like source with 36 counts
and a more extended component with 117 counts
(Fig.~\ref{fig:n31n54}). Assuming a power-law spectrum with slope
$\Gamma \sim 2$ and correcting for Galactic absorption only, the point
source has luminosity $L(0.3-8\mathrm{keV}) \sim 2\times 10^{37}$
erg s$^{-1}$. Luminosities in this range are typically associated with XRBs
rather than AGNs, but does not rule out an AGN, especially since
intrinsic absorption is unknown. (b) \emph{Infrared line ratios are
not that of an \textsc{Hii} nucleus but instead suggest the presence
of an AGN component}. A stronger argument that the source is an AGN
derives from infrared data. This galaxy is part of the
\textit{Spitzer} Infrared Nearby Galaxies Survey
\citep{kea03}. \citet{dea06} have used the equivalent width of the PAH
feature at $6.2 \mu$m and the fluxes in a mix of high- and
low-ionization IR lines
of S, Ne, and O to create diagnostic diagrams that
distinguish between AGN and star-forming galaxies. This nucleus falls
into the ``transition'' region between AGNs and \textsc{Hii}
regions. This suggests that the emission has an AGN component that
was diluted because of the large aperture used ($\sim 20''$)
to extract the fluxes. (c) \emph{Mid-infrared colors are redder than
those of normal galaxies}. The observed IRAC colors, $[3.6]-[4.5] =
+0.20$ $\pm 0.16$ and $[5.8]-[8.0] = +1.24$ $\pm 0.16$ (magnitudes in
Vega system) (D.~A.~Dale, priv.\ comm.), are redder than more than
80\% of normal late-type galaxies (Fig.~\ref{fig:n3184}; also see
\cite{aea08}). This is expected if there is an AGN, as AGN power-law
emission falls off more slowly than galactic emission in the
MIR. Thus, the IR line ratios and IRAC colors strongly argue that the
source is an AGN.

\begin{figure}
  \includegraphics[height=0.2\textheight]{n31n54_gs.ps}
\caption{\textit{Chandra\/}\ images of the nuclei of NGC 3184 (left) and NGC 5457
  (right). Both panels are $10''$ on a side and have north up and east
  to the left. \textit{NGC 3184\/}: The nucleus consists of a point-like
  source and an extended component to the south of it. The bar on the
  lower left corresponds to a distance of $\sim\!80$ pc. \textit{NGC
  5457\/}: The nucleus is resolved into two sources, marked ``N'' and
  ``S'' in the figure. Source N is the active nucleus and source S is
  a star cluster. The bar on the lower left corresponds to a distance
  of $\sim\! 70$ pc.}
\label{fig:n31n54}
\end{figure}

\begin{figure}
  \includegraphics[height=.3\textheight]{IRAC_color5.eps}
  \caption{The unusually red mid-infrared color of the nucleus of NGC
  3184 is marked on a Spitzer IRAC color-color diagram. The contours
  show the distribution of $z\sim 0$ galaxies, separated by factors of
  1.8 in galaxy density. C1...C4 stand for the 3.6 $\mu$m, 4.5 $\mu$m,
  5.8 $\mu$m, 8.0 $\mu$m bands respectively. The left- and right-hand
  concentrations are spiral and elliptical galaxies respectively.  A
  typical error bar is shown in the lower right. (R.~J.~Assef 2007,
  priv. comm. Also see \cite{aea08}). \textit{Reprinted with
  permission of the author.}}
\label{fig:n3184}
\end{figure}

\section{NGC 5457}

NGC 5457 is a galaxy of type Scd at a distance of approximately 7 Mpc
\citep{ssea98}. The nucleus was classified as \textsc{Hii} by
\citet{hfs97-3}. NGC 5457 was observed multiple times by \textit{Chandra\/}\ for a
total observing time of approximately 1 Ms
(Fig.~\ref{fig:n31n54}). For this galaxy the x-ray data by themselves
provide evidence suggesting the presence of an AGN. (a) \emph{The
nucleus is an x-ray source}. (b) Even more importantly, \emph{the
nuclear source is variable}. The nucleus varies by about a factor of
10 over nine months (Fig.~\ref{fig:n5457}), with luminosity ranging
from $L(0.3-8\mathrm{keV}) \sim 3 \times 10^{37}$ to $3\times 10^{38}$
erg s$^{-1}$. The variability rules out multiple XRBs as the origin of the
x-ray emission, as they would have to be varying in phase. It is
possible for a single XRB to produce the observed x-ray emission, but
XRBs with luminosities this high are uncommon. (c) \emph{The x-ray
spectrum} shows an absorbed power-law with slope $\Gamma \approx 2$,
similar to what is seen in AGNs. The best-fit intrinsic absorption is
consistent with zero when the source is faint. (d) \emph{The ratios of
0.3--2 keV, 2--5 keV, and 5--8 keV counts} for this nucleus are
similar to those seen in Compton-thick AGNs \citep{lea06}. Points (c)
and (d) together suggest a highly obscured AGN where, like in the
Seyfert galaxy NGC 1068, the dominant observed emission is radiation
that has been scattered into our line of sight. Finally, (e) as can be
seen in Fig.~\ref{fig:n5457}, \emph{the behavior of the nucleus is
markedly different from that of a known star cluster} of comparable
x-ray luminosity which is $\sim\! 110$ pc away from the nucleus. This
suggests the nuclear source is not simply another star cluster.

\begin{figure}
\includegraphics[angle=-90,scale=0.75]{vari_ctsonly.ps}
\caption{X-ray variability of the nucleus of NGC 5457 as seen by
\textit{Chandra\/}. The filled squares show the nuclear count rate, and the open
squares the count rate from a nearby star cluster. The nucleus
varies by almost a factor of ten between March and November, 2004.}
\label{fig:n5457}
\end{figure}

\section{Conclusion}

There are compelling, though not conclusive, arguments that NGC 3184
and NGC 5457 actually host AGNs even though neither shows signs of
AGNs in their optical spectra. These AGNs are excellent candidates for
being low-mass SMBHs, as they reside in low-mass bulges. A similar
approach can be used to uncover more candidate low-mass SMBHs.

\begin{theacknowledgments}
We are grateful to D.~A.~Dale for kindly providing \textit{Spitzer}
fluxes for the nucleus of NGC 3184 prior to publication. Support for
this work was provided by the National Aeronautics and Space
Administration through Chandra Award Number GO7-8111X issued by the
Chandra X-ray Observatory Center, which is operated by the Smithsonian
Astrophysical Observatory for and on behalf of the National
Aeronautics Space Administration under contract NAS8-03060.

\end{theacknowledgments}

\bibliographystyle{aipproc}

\input{confproc.bbl}

\end{document}