active galactic nuclei pose a fundamental problem  how does rotationally dominated gas on kpc scales lose orders of magnitude in angular momentum to fall onto an accretion disk on au scales ( e.g. , * ? ? ?
we do not know the mechanism that dissipates angular momentum and allows gas to accrete .
there are no shortage of ideas , including major or minor mergers ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) , bars and bars within bars ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) , and nuclear spirals @xcite as perhaps indicated by dust lanes @xcite . even looking directly at molecular gas kinematics on pc scales in active galaxies , it is very hard to find clear evidence for the true driver of nuclear activity @xcite .
thanks to the advent of near - infrared integral - field spectrographs ( nir ifu hereafter ) , there has been significant progress in understanding gas flows in the centers of active galactic nuclei ( agn ) in recent years .
the nir observations allow us to penetrate gas and dust in the active nuclei , and , in conjunction with adaptive optics , to probe very near to the nucleus .
in some nearby cases , inflows are directly observed along circumnuclear spirals in ionized gas on @xmath2 pc scales ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
furthermore , there are intriguing hints of kinematic differences in the nuclei of local seyfert galaxies .
for example , @xcite suggest that the ionized gas in the disks of seyfert galaxies ( i.e. , outside of the narrow - line region or nlr ) is more kinematically disturbed than the gas in a control inactive sample . @xcite
also report that seyfert galaxies have more concentrated stellar luminosity profiles , lower stellar velocity dispersions , and elevated h@xmath4  1 - 0 s(1 ) emission within @xmath5 pc relative to an inactive subsample .    in this work ,
we specifically focus on megamaser disk galaxies . in these special systems , luminous water megamasers at 22 ghz trace a geometrically thin molecular disk on sub - pc scales in orbit around the central bh ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
the megamaser spots show near perfect keplerian rotation in many ( e.g. , * ? ? ?
* ) , but not all cases ( e.g. , * ? ? ?
* ; * ? ? ?
thanks to the keplerian rotation in the maser disk , we know the bh mass to a few percent , limited in most cases by the distance to the galaxy ( e.g. , * ? ? ?
direct geometric distances have also been derived by measuring the accelerations of the systemic masers for the few best cases @xcite .
in addition to distance measurements and bh masses , the megamaser disk reveals the spin axis of the accretion disk on sub - pc scales . to get the long path lengths required for masing ,
the disk must be virtually edge - on ; we thus typically know the orientation of the disk on sub - pc scales to within a couple of degrees @xcite .
we have already found that the spin axis of the megamaser disk aligns with the jet on @xmath2 pc scales , while the spin axis of the disk does not appear to align with flattened disk - like structures on @xmath6 pc scales identified in _ hst _  imaging @xcite . in this paper
, we explore the use of nir ifu data to define the gas flows on @xmath7  pc scales , for comparison with the megamaser disk .
we currently have nir ifu observations of the centers of two galaxies in the kuo et al .
sample : ngc 4388 and ngc 1194 .
these two are the only two megamaser disk galaxies with published bh masses where we can come close to resolving the gravitational sphere of influence , in the former case because the galaxy is nearby ( in virgo ) and in the latter because the bh is the most massive known with a maser disk @xcite .
we observed ngc 4388 with sinfoni under program 386.b-0244a on the very large telescope @xcite on april 12 & may 5 2011 in natural seeing conditions and in service mode .
we have been awarded time to observe this galaxy at two spatial resolutions , the higher one involving assistance of the adaptive optics system ( ao ) , however to date these observations are only partially carried out .
thus , in this paper we present the lower spatial resolution data set that covers a wider field of view .
our observations cover the central 8@xmath8 8 , with a spatial sampling of 0125 @xmath8 0250 .
we used the @xmath9-band grating ( 1.95  2.45 @xmath10 m ) that gives a spectral resolution r@xmath113500 ( 6.2   fwhm , as measured on sky lines ) .
for the observations we used standard near - ir nodding techniques . observing blocks consisted of a sequence of object and sky frames ( osoosoosoosoo ) .
each individual integration was 300 sec and the sky fields were offset by 240  to the north and east .
science frames were dithered by 03 in order to reject bad pixels .
there were two observing blocks , with a total on - source integration time of 100 min , and 40 min in total on sky frames .
additionally , after each observing block and at a similar airmass , we observed a b dwarf star to act as a telluric star .
we used the eso sinfoni pipeline v2.3.3 to perform the basic data reduction on each observing block , consisting of 10 object and 4 sky exposures . in brief
, the pipeline extracts the raw data , applies distortion , bad pixel and flat - field corrections and wavelength calibration , and stores the combined sky - subtracted spectra from one observing block in a 3-dimensional data cube .
the final pixel size is @xmath12 .
we reduced the telluric stars in the same way as the science frames . for each telluric star we extracted a one - dimensional spectrum , removed the hydrogen brackett@xmath3 absorption line at @xmath13 m after fitting it with a lorentzian profile , and divided the star spectrum by a black body spectrum with the same temperature as the star . the last step in preparing
the telluric spectrum was to apply small shifts in the spectral direction ( @xmath140.05 pixels or 0.123  ) and scalings to minimize the residuals of the telluric features . to do this , we extracted a central one - dimensional spectrum from each science data cube and cross - correlated and fitted it with the corresponding telluric spectrum .
we derived the wavelength - dependent correction required to match the continuum shape of the telluric star ( without an absolute zero - point ) .
thus , when we divided each individual spaxel in the six galaxy data cubes by the corresponding best - fitting telluric spectrum , we also obtained a relative flux calibration .
based on the seeing monitor , during our observations the seeing varied from @xmath15  for the first observation , and between @xmath16  for the second observation .
note that the seeing monitor takes measurements at 5000   , thus our seeing is typically better than this . from the standard stars ( taken at the end )
we measure a fwhm of 08 for the first observation and 05 for the second .
we aligned the two obs with integer offsets and summed them ; all analysis was done on this final summed data cube .
we observed ngc 1194 using the ifu oh - suppressing infrared imaging spectrograph ( osiris * ? ? ?
* ) assisted by the laser - guide star adaptive optics system @xcite on the 10 m keck ii telescope . the data were taken over the latter two - thirds of a single night on oct .
we observed in the @xmath0band with the kbb filter , for a spectral coverage of @xmath17 and a spectral resolution of @xmath18 .
we used the @xmath19 spatial scale , and oriented the long axis of the ifu along the major axis of the galaxy at a pa of @xmath20 .
we dithered the galaxy by @xmath21 perpendicular to the long axis of the ifu , both to facilitate bad pixel removal and to slightly increase the field - of - view to @xmath22 .
unfortunately , the nucleus of ngc 1194 was too faint to use as a tip - tilt star .
thus , we used a nearby ( 54 ) star ( id=0888 - 0029937 ; @xmath23 ) .
natural seeing conditions varied over the night from @xmath24 to @xmath25 . from our final science data
cube we find that the psf core had a fwhm of @xmath26 as measured from the galaxy core .
we also observed a telluric standard ( a0v ) star approximately once per unit change in airmass .
finally , we interleaved sky and object observations with an object - sky - object pattern , spending 600s per sky and object exposure .
we acquired a total of 170 min on source .
we followed closely the reductions outlined in @xcite , using the osiris data reduction pipeline ( v2.3 ) made available by the instrument team .
we repeat salient details here for completeness .
the pipeline handles sky subtraction , cosmic - ray removal , glitch identification , spectral extraction into a data cube , wavelength calibration , atmospheric dispersion correction , and telluric correction using an a0v star .
we experimented with a scaled sky subtraction method from @xcite , but did not find substantial improvements from the direct method .
finally , we determined the centroid of each exposure by hand , applied these sub - pixel offsets , and averaged all cubes into our final data cube .
here we describe our analysis procedure for the ngc 4388-sinfoni data cube in detail , highlighting any differences in analysis with the ngc 1194 data cube .
we use the direct - pixel  fitting code ppxf of @xcite to fit the stellar continuum . to achieve adequate signal - to - noise ( s / n ) we use the voronoi binning scheme of @xcite , which determines contiguous bins with uniform s /
n , in this case s / n@xmath27 per pixel .
the sinfoni pipeline does not return an error array , and so the noise is determined empirically from line - free regions of the spectrum as the variance in the spectrum once large outliers are clipped . in the case of the osiris data , we use the per - pixel error array generated by the reduction pipeline .
we then fit the coadded spectrum corresponding to each bin with a weighted combination of stellar templates .
each template is shifted to the systemic velocity of the galaxy and convolved with a gaussian line - broadening function .
a polynomial of fourth order is added to account for nonthermal continuum from the agn , and potentially flux calibration errors as well .
example fits are shown in figure [ fig : specn4388 ] .
the measurement uncertainties in both velocity and velocity dispersion are determined with monte carlo simulations . in the case of the ngc
1194 data , the data reduction software generates an error array
. however , for the ngc 4388 data , we use residuals from the best fit ( with the emission lines removed ) to determine the average sigma .
we create 100 perturbed input spectra assuming gaussian errors in the spectra .
we then refit these 100 artificial spectra .
the uncertainty in each parameter is calculated as the values encompassing 68% of the trials .
we experimented with a number of fitting regions , including a fit to the bandhead only ( @xmath28 , short ) , a fit to the full region ( @xmath29 , full ) , and an intermediate region ( @xmath30 , best ) .
we found that the first fit delivered reliable velocity dispersions , but unstable radial velocities because of the narrow wavelength range , and because the bandhead at @xmath31 is often filled in by [ ] @xmath32 emission .
fits to the full spectral region suffered from telluric residuals at @xmath33 .
the just - right spectral region was therefore the third , intermediate case .
the median errors on the velocities are @xmath34  km  s@xmath35  for the best , short , and full regions respectively , while the median errors on the dispersion measurements are @xmath36  km  s@xmath35 . in all cases , we masked high - ew emission lines , including the h@xmath4  lines at @xmath37 @xmath38 and also [ ] @xmath39 , @xmath40 , br@xmath41 , and [ ] .    to bracket uncertainties caused by template mismatch
, we utilize two stellar template libraries in fitting the sinfoni data .
one is a library of stars observed with sinfoni using the same observational set - up . in this case , no correction for instrumental resolution is required .
the native sinfoni templates range from k4iii to m5iii ( plus a k4.5v star ) @xcite . as a check
, we also use the stellar template library of @xcite , observed at higher spectral resolution .
we use this second template set only to test our sensitivity to template mismatch .
the wallace templates cover a wider range in spectral type , including km supergiants and km giants .
we discuss the stellar template fits to the ngc 4388 data in detail in s 4.1 .
even with different templates , we find good agreement in the two sets of stellar velocity dispersion measurements .
taking @xmath42 as the answer based on the sinfoni templates and @xmath43 as the same for the wallace templates , we find @xmath44 ( where the latter number is simply the standard deviation in this ratio ) .
we find that the two measurements agree within 20% ( i.e. , within the observational uncertainties ) and we find no significant systematic offset . in the case of the osiris data ,
our primary fitting uses a single k5iii star , which provides an acceptable fit to our moderate s / n data .
we do not a priori know the location of the black hole , or the precise photometric center of the galaxy .
we could use the continuum to determine the photometric center , but then we would be sensitive to obscuration , which is still significant in the nir ( see figure [ fig : image ] ) .
instead , we use the equivalent width of the co bandhead to determine where the active nucleus peaks @xcite .
we use the index definition from @xcite : the index band is @xmath45 , with a continuum band centered on 2.2900 and a width of 0.0003  @xmath46 , and we also followed their prescription to correct for velocity dispersion although this correction is negligible ( see also * ? ? ? * ) .
these authors show that in the absence of an active nucleus , star - forming galaxies have a very uniform co ew .
a declining co ew towards the galaxy center can be attributed to infill by the nonthermal continuum from the accreting black hole .
we show an example in  [ sec : stellarkin4388 ] .
we create an agn continuum map by taking the co ew at the edge of the cube ( ew=9.5  in ngc 4388 , consistent with stellar - population models ) and assuming that all dilution further inwards results from the agn , such that continuum ( agn ) = [ total continuum ] @xmath8 [ 1 - co ew / co outer ] .
the photometric center of the continuum map falls @xmath47 south of the co ew center .
the optical narrow - line region cone is only visible to the south (  [ sec : nlr ] ) , so it is sensible that the continuum centroid should be shifted towards the low - reddening side of the disk .
the best - fit center also agrees within less than a pixel with the peak h@xmath4  and br@xmath3  emission , although these latter are not as well defined .
we will use this position as the center of the cube throughout .
in the case of ngc 1194 , the co - derived and photo - centers agree within @xmath48 .      while voronoi binning is very useful to create uniform continuum s / n ( e.g. , * ? ? ? * ) , it is more difficult to apply to patchy emission - line maps , where by using large bins we may smear out interesting features in the line emission maps . we therefore create a second data cube using a @xmath49 pixel smoothing in the outer regions and a @xmath50 pixel binning in the inner @xmath50region . with this level of binning , we increase the s / n by a factor of two but still maintain multiple pixels across the psf .
we measure the emission lines from this uniformly gridded cube . in the case of ngc 4388 , because we have such high s / n in the continuum , we again perform continuum - subtraction using ppxf before fitting the emission lines . however , there are no strong stellar absorption features beneath the emission lines of interest , and thus for ngc 1194 , where the s /
n is poor in the continuum , we simply fit a local continuum value .
we fit each emission line independently .
we model the intrinsic line shape as a gaussian , although we allow for up to two gaussian components for emission lines with s / n@xmath51 in ngc 4388 . particularly in the br@xmath3  and
[ si vi ] lines , the line shapes are often non - gaussian and warrant a second component . in @xmath52 of cases within the central @xmath53 the s / n is sufficient to warrant a two - component fit .
we calculate the non - parametric fwhm from the sum of gaussian components , and from that value a corresponding line dispersion .
we also experimented with an empirically determined line - broadening function by fitting a gauss - hermite model to the strong sky emission lines , which caused only very minor changes in our fits .
we derive errors on all the fitted parameters by fitting 100 artificial spectra that are created as a sum of the original spectrum and a gaussian random noise array generated from the extracted errors .
ngc 4388 is an sbb galaxy at a distance of 19 mpc @xcite , corresponding to a physical resolution of @xmath1 pc given our @xmath54 seeing .
the galaxy magnitude is @xmath55 mag , with a b / t@xmath56 based on fits to the sdss data @xcite .
we show an hst / wfc3 image combining f336w ( broad @xmath57 ) , f435w ( broad @xmath58 ) , and f814w ( broad @xmath59 ) in figure [ fig : image ] . with a photometric p.a .
of 90 , the galaxy is very nearly edge - on , which is not true for the majority of the megamaser disk galaxies @xcite .
the inclination determined from the large - scale kinematics is 78@xcite , and the systemic velocity from the same study is @xmath60  km  s@xmath35 .
the gas kinematics show an isophotal twist on arcmin scales , that @xcite successfully model as a bar within the inner 1.5 kpc , with a p.a . of 100  on the sky .
the bar is also apparent in the photometry as a boxy bulge .
ngc 4388 is well - known as a galaxy that is falling into the virgo cluster @xcite .
gas is being stripped from the galaxy , as evidenced by an and ionized gas tail @xcite .
there is even a detection of the photoionized gas in the soft x - rays @xcite .
we are interested in the center of ngc 4388 , in particular because of the megamaser disk on sub - pc scales @xcite .
based on the megamaser rotation curve , the bh mass is found to be @xmath61  @xmath62  @xcite .
combined with the bolometric luminosity estimate from @xcite , the eddington ratio is @xmath63@xmath64 , which is in very rough agreement with our estimate from the [ ]  luminosity @xcite . like
most megamaser disk galaxies @xcite , the active galactic nucleus ( agn ) in ngc 4388 is compton thick and quite bright in hard x - rays ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
ngc 4388 has also been seen to change state in the x - ray @xcite , switching from being compton thick to compton thin .
weak broad h@xmath65  has also been reported at the galaxy center @xcite .
we first examine stellar rotation and dispersion .
we clearly see rotation within our @xmath66 pc field of view that aligns with the disk on larger scales .
following @xcite , we first use ppxf to fit the average spectrum over all data , excluding a 05 radius around the agn ( figure [ fig : specn4388 ] ) . in the fiducial fits , we do not fit higher - order gauss - hermite moments . we will examine @xmath67 and @xmath68 in  [ sec : stelrotcurve ] below .
all of our default fits employ the sinfoni templates . from this high s
/ n spectrum we can derive the best - fit spectral mix , which may be degenerate in lower s / n individual spectra . based on our fit to this average spectrum , we find a dispersion of @xmath69  km  s@xmath35 , which compares well with our previous measurement of @xmath70  km  s@xmath35@xcite .
the best - fit template is composed of 50% k4/5iii , 38% k7iii , and 12% m5 ii / iii stars .
we then fit spatially resolved spectra across the cube , using voronoi binning to ensure comparable s / n in all bins .
we both allow the mix of templates to be a free parameter , and use the best - fit mixture from the high s / n fit to the entire cube .
the latter fit has the advantage that by decreasing the number of free parameters , we increase the fitting stability in the agn - dominated region @xcite .
however , we find excellent agreement between the two measurements , with the scatter in their fractional difference being only 9% .
thus , we will present results with the stellar templates allowed to vary .
the fitted rotation and dispersion fields are shown in figure [ fig : sigdispfield ] .
our best - fit systemic velocity is @xmath71  km  s@xmath35 , consistent with the measurement of @xcite .
we do find a systematic difference of @xmath72  km  s@xmath35  between the absolute velocity derived when we use the sinfoni stellar templates and that from the wallace et al .
we are not sure of the origin of this difference ; it is possible that the sinfoni template stars , while all at rest relative to each other , have not been shifted to a rest velocity .
we thus quote an uncertainty of @xmath73  km  s@xmath35  in the systemic velocity .      to derive the rotation curve we use the _ kinemetry _ formalism of @xcite , which is similar to ellipse fitting of photometry , but operates on both even and odd moments of the kinematic field .
the kinematics are modeled using a sixth order expansion along ellipses with the position angle and ellipticity as free parameters .
we show the resulting rotation curve in figure [ fig : stelrotcurve ] , along with the best - fit ellipse pa , ellipticity @xmath74 , and the amplitudes of the higher - order terms that encapsulate deviations from simple circular rotation .
we clearly measure rotation in the stellar kinematics , with a maximum amplitude of @xmath75  km  s@xmath35 .
the stellar rotation axis of 90  is aligned with the major axis of the kpc - scale disk . on larger scales ,
@xcite measure a stellar rotation of 60 - 80  km  s@xmath35  out to @xmath76  ( @xmath66 pc ) , consistent with ours , which then appears to fall slightly at larger radius .    with @xmath77@xmath78
km  s@xmath35  beyond 1(100 pc ) , we find a @xmath79 .
overall , the stars are dispersion - dominated , but we see signs of kinematic components apart from just a bulge .
firstly , for a @xmath79 , we would expect an ellipticity of @xmath80 for an isotropic oblate rotator @xcite , rather than the @xmath81 that we observe in the isophotes ( figure [ fig : image ] ) .
secondly , within a radius of @xmath82 the dispersion field actually drops ( although the decrement is more pronounced on the east side of the disk ) . on the same scale
, we see a jump in velocity , pa , and ellipticity likely signaling a transition from an inner disk to more bulge - dominated kinematics .
thirdly , we also detect a disk - like structure in the inner 100 pc in the _
hst_/wfc3 f160w image ( figure [ fig : image ] ) .
the nuclear disk appears to have pa@xmath83 , misaligned with the large - scale disk , and a size very similar to the @xmath84drop region .
we detect this same disk component in more detailed two - dimensional fitting of the _ hst_+ground - based data ( r. lsker et al .  in preparation ) .
the most likely explanation of these observations is that there is a disk within @xmath85  ( 180 pc ) embedded in the larger - scale bulge / bar . only on unextincted sight lines
are the kinematics dominated by the disk component , while on heavily extincted sight lines we are seeing the kinematics in the bulge / bar .
the falling rotation curve outside 400 pc may also reflect the increasing dominance of the bulge / bar on larger scales . as additional confirmation of this picture
, we refit the stellar kinematics and measure the higher - order gauss - hermite moments @xmath67 and @xmath68 ( e.g. , * ? ? ?
although individual @xmath67 measurements have only s / n @xmath86 , we do find the classic anti - correlation between @xmath67 and @xmath87 that points to the superposition of a bulge and disk , as seen both in simulations @xcite and data ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
we also find two additional hints that the @xmath0band extinction is higher to the west .
firstly , we find a significantly redder continuum slope to the west side of the disk , corresponding to a differential @xmath88 mag assuming a @xcite reddening law .
secondly , we clearly see a dip in the h@xmath4  emission on the same scale on the west side of the disk ( see  [ sec : gasflux ] ) .
next we will ensure that the @xmath84drop is real and not an artifact of our fitting procedure .
the agn is the primary source of uncertainty in our modeling of @xmath77  in the galaxy center .
contamination from the non - thermal central source dilutes the stellar signal , and is degenerate with @xmath77  ( e.g. , * ? ? ? * ) .
first , we look at the s / n of the dispersion measurements as a function of position . we find that the s / n is always @xmath89 except within @xmath21 of the agn .
this first test suggests that our measurements at @xmath90  ( 100 pc ) should not be dominated by the agn .
next , we re - run ppxf , but instead of fitting the continuum with the default polynomial of order 4 , we only allow a first or second - order polynomial ( figure [ fig : sigdispfield ] ) . while it is clear that , particularly in the case of the first - order polynomial , we achieve very poor fits to the nuclear region ( @xmath91 ) , we still find significant evidence for an asymmetric @xmath84drop .
we then return to the wallace templates ( rather than our default sinfoni templates ) and try restricting the template set to k dwarfs .
we find the same result .
finally , we examine our fits with the template fixed to the mixture derived from the high s / n fit described in  [ sec : stellarkin4388 ] . in all of these cases , although the dispersion profile within 05 is not well determined , we recover the @xmath84drop on 1  ( 100 pc ) scales .
we therefore conclude that the observed drop is real and that the asymmetry results from patchy reddening .
although in general the north side of the disk is more heavily extincted (  [ sec : nlr ] ) , there are also dust lanes extending to the south , which are apparently preferentially affecting the west side of the 100-pc  scale disk .
a falling dispersion towards the center of megamaser disk galaxies is not unexpected .
these so - called `` @xmath92-drops '' are common both in spiral ( e.g , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) , and in early - type galaxies @xcite , in the latter often associated with nuclear star clusters @xcite .
the typical scale is @xmath93 pc , and
the most likely explanation for the falling dispersion is that the galaxy harbors a gas disk that is currently or was recently forming stars .
these bright , dynamically cold stars dominate the central dispersion and cause @xmath77  to fall .
in fact , ngc 2273 , one of the nearer megamaser disk galaxies , has a dramatic @xmath84drop that coincides with a known disk in the inner @xmath6 pc @xcite .
@xcite suggest that low central dispersions are more common in active galaxies .
if the @xmath84drop is caused by a young stellar disk , then we might also expect to see younger stars coinciding with the lower stellar velocity dispersion .
we create a median spectrum over the region with @xmath77@xmath94  km  s@xmath35 , and rerun ppxf to investigate the mixture of stellar populations using the wallace templates , which cover a wider range of spectral type . while there is a slight preference for more supergiant light in the region of the @xmath84drop , we do not find a significant difference in best - fit spectral type as a function of region .
the co ew profile is also quite symmetric and has reached its asymptotic values already in the @xmath84drop region ( figure [ fig : sigdispfield ] ) . with a larger spiral sample ,
@xcite also find similar stellar population ages inside and outside of the @xmath84drop regions .
we propose that the observed stellar kinematics are the superposition of three distinct components . on kpc - scales , the galaxy disk has pa@xmath95 . within our 300 pc aperture ,
the kinematics are dominated by the hotter bulge / bar @xcite .
based on the highest @xmath77  measurements within the cube , presumably occurring where dust obscuration is highest , we measure @xmath77@xmath96  km  s@xmath35  in the bulge / bar . then , within the inner 100 pc , we see evidence in both the kinematics and the nir isophotes for a nuclear disk component at a pa@xmath97 , misaligned by @xmath98  from the large - scale disk and with a @xmath99 .
finally , on pc scales , there is a masing disk at pa@xmath100 , which is misaligned not only from the kpc - scale disk , but also by @xmath101  from the nuclear stellar disk @xcite .      with a best - fit stellar rotation field in hand ,
we turn to the gas kinematics .
we will examine three transitions : h@xmath4@xmath102 1 - 0s(1 ) , br@xmath3 , and the coronal line [ ] @xmath103 .
we note that the [ ca viii]@xmath104 line , also coronal , has been used to trace the nlr ( e.g. , * ? ? ?
* ) , but is blended with the co bandhead . thus we choose to focus on the [ ] line here .
the two - dimensional velocity and line - ratio fields will provide information about the origins and excitation mechanisms of each transition . the fluxes that we present are the integrated flux from our one or two gaussian models as described in  3.3 .
we calculate the fwhm non - parametrically from the sum of both components , and the dispersions are simply calculated as fwhm/2.35 .
we also measure and present the skewness of the lines , which encapsulates the asymmetry in the two - component fits .    in figure
[ fig : emissionflux ] , we show the distribution of fluxes in the h@xmath4 , br@xmath3 , and [ ] emission lines , while in figures [ fig : h2rot ] , [ fig : oiiisivi ] , and [ fig:2dgas ] , we summarize their kinematics .
we see that the overall morphology of the h@xmath4follows the stellar continuum , and is even more flattened than the stellar distribution , suggesting that the emission arises from the nuclear disk .
furthermore , close inspection reveals that the h@xmath4emission is asymmetric , being considerably brighter towards the east .
the brighter side of the h@xmath4  disk is coincident with the observed @xmath84drop in the stellar distribution ; we can directly see the extinction of the h@xmath4  in the region where the @xmath84drop is also extincted .
in contrast , the high - ionization lines are oriented at pa@xmath105 , aligned with the jet on similar scales , and also with the [ ]  emission that traces the narrow - line region ( nlr ; fig .
[ fig : oiiisivi ] ) .
the br@xmath3  mostly follows the high - ionization lines ( fig .
[ fig:2dgas ] ) .
their light distribution strongly suggests that their emission arises in front of the disk in projection , since the [ ] is actually brighter towards the west ( similar to the [ ] ) .
based partially on these flux distributions , as well as the kinematics and line ratios presented below , we will suggest that the h@xmath4  mainly traces star formation in the nuclear disk , while br@xmath3  and [ ] trace the nlr ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) .
we focus first on the h@xmath4@xmath106 kinematics ( figure [ fig : h2rot ] ) . as with the stellar kinematics above
, we fit the h@xmath4  velocity field using kinemetry .
the best - fit rotation curve , pa , and eccentricity are shown in blue in figure [ fig : stelrotcurve ] . with v/@xmath107 , the gas is considerably colder than the stars , but we are not seeing pure rotation .
the rotation field shows an @xmath108-shaped kink .
the twist is apparent as a discontinuity in the h@xmath4  ellipticity @xmath74 and pa in the rotation profiles at @xmath109 ( 150 pc ) , a similar radius to the changing stellar kinematics . in the skewness map , along the high - dispersion kink
, we see a tendency for red asymmetry , which is plausibly attributed to the superposition of two velocity components .
interior to the twist , the kinematic pa of the gas on @xmath110 pc scales appears to align with the megamaser disk .
thus , it is possible that whatever mechanism causes this kink is responsible for supplying material on pc scales .    on the same scales ,
we find a two - armed spiral in the dispersion map , with an inner pa @xmath105 .
we find a corresponding three - arm spiral pattern in the residual map by subtracting the best - fit rotation model from the observed velocity field . to interpret these trends , we recall that any @xmath111 mode perturbation can cause @xmath112 deviations in the even velocity moments , including density and dispersion , and due to geometry @xmath113 ( and @xmath114 ) deviations in the odd velocity moments ( e.g. , * ? ? ? * ) .
we propose a few mechanisms that may explain these kinematic properties .
first , in the study by @xcite of the inner region of ngc 1097 , similar signatures were explained as the result of a two - arm nuclear spiral in the disk plane ( see also * ? ? ?
given that the nuclear disk in ngc4388 is nearly edge - on , such a co - planar spiral density wave seems unlikely .
the second possibility is that we are seeing the superposition of the kpc - scale disk and the 100 pc nuclear disk that we saw in the stellar kinematics .
however , the gas pa@xmath115 , while approximately aligned with the maser disk on sub - pc scales , does not align with the stellar disk of pa@xmath97 .
the final possibility is that the twist we observe here is related to the kink from the bar observed on kpc scales ( e.g. , * ? ? ?
* ) . only with our upcoming ao observations
can we address the @xmath116 pc gas kinematics definitively , and hopefully help determine what is driving the gas inward towards the agn .
the h@xmath4  shows very regular rotation along the orientation of the large - scale galaxy disk outside of the inner 1  ( 100 pc ) .
additional information is encoded in the ratios of h@xmath4  to br@xmath3  ( figure [ fig : h2brg ] ) , and ratios of different h@xmath4transitions ( figure [ fig : h2rat ] ) .
typically , in the hard ionizing field of the nlr , the ratio of h@xmath4  to br@xmath3  takes a value @xmath117h@xmath4/br@xmath3@xmath118 , while it is @xmath119 in star - forming regions @xcite .
finally , in dense spiral arms , even higher ratios of h@xmath4/br@xmath3@xmath120 are observed .
@xcite suggest that the increased h@xmath4  emission arises from the higher gas densities within spiral arms ( see also * ? ? ?
ngc 4388 obeys these trends in general ( figure [ fig : h2brg ] ) . along the disk ,
the h@xmath4  is strong , and the ratio has values @xmath121 , as is commonly seen in spiral arms , while in the galaxy center , the line ratio ranges from @xmath122h@xmath4/br@xmath3@xmath123 , values common for seyfert galaxies .
the only surprising region lies off the disk along a @xmath124angle east of north above and below the galaxy plane , where the br@xmath3emission is strongest . here
, the h@xmath4/br@xmath3  ratio is very low @xmath125 .
such low ratios are reportedly associated with star - forming regions ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
however , in this case , the br@xmath3  and [ si vi ] morphologies match both the orientation of the radio jet and the nlr traced by [ ] @xmath126 in the optical .
thus , we suggest that the low h@xmath4  relative to br@xmath3  in these regions is instead due to very low densities of the nlr @xcite .
we now construct spectra in characteristic regions to examine both the h@xmath4  line ratios and the line kinematics .
we select regions of the map based on the h@xmath4/br@xmath3  line ratio and construct high s / n spectra from the continuum - subtracted cubes .
region 1 is the center , defined as all pixels with co ew@xmath127  that lie within @xmath109 of the center .
we extract two regions in the disk ( h@xmath4/br@xmath3@xmath128 ; regions 2 & 3 ) and two regions in the nlr ( h@xmath4/br@xmath3@xmath129 ; regions 4 & 5 ) . finally , we sample gas that falls beyond the main galaxy and overlaps with the [ ]  ionization cone ( region 8) , as well as two regions north of the disk , towards the northern radio lobe ( regions 6 & 7 ) .
the extracted spectra for the strong transitions ( h@xmath4@xmath130 , br@xmath3 , and [ si vi ] ) from each region are shown in figure [ fig : h2brg ] , the spectra from different h@xmath4  transitions are shown in figure [ fig : h2rat ] , and measured properties from these spectra are presented in table 1 . because we fit the lines using an instrumental profile , the linewidths are naturally corrected for instrumental resolution .
each line is fitted independently .
in fact , as shown in figure [ fig : h2rat ] , the centroids and line shapes for the different transitions can differ , particularly for regions in the nlr ionization cone .    for each region , we can ask about the excitation mechanisms of the h@xmath4 .
there are two ways to excite the warm h@xmath4 : thermal processes or fluorescence via uv radiation .
the thermal processes that can in principle heat the gas include shocks , x - ray emission , or heating by uv emission .    the predicted line ratios of various h@xmath4  transitions provide diagnostics of the primary excitation mechanism .
specifically , if fluorescence dominates the excitation , then the ratio of h@xmath4  2 - 1s(1 ) 2.25@xmath131/1 - 0s(1 ) 2.12@xmath131 is expected to be @xmath132 while the 1 - 0s(2 ) 2.03@xmath131/1 - 0s(0 ) 2.22@xmath131 ratio is expected to be @xmath90 @xcite .
we show this pair of ratios and the expectation for pure thermal emission as the solid line , while the region dominated by fluorescence is shown schematically with the large x at the lower right in figure [ fig : h2rat ] .
there is room for some uv fluorescence , but these line ratios suggest that the bulk of the excitation is thermal .
similar conclusions have been drawn for larger samples of seyfert galaxies @xcite .
our excitation diagram is constructed for the integrated line profile but , in fact the line ratio is clearly a function of velocity , particularly for the line emission arising from the nlr ionization cones .
we do not have sufficient s / n to construct line ratios for the line core and line wings separately , but these different line shapes are pointing to different dominant excitation mechanisms in the line wings . higher s / n is needed to quantitatively examine velocity - dependent line ratios .
we are still left with the question of which thermal process dominates the excitation .
previous work has drawn various conclusions , with @xcite for instance favoring shocks in outflows , @xcite favoring a mix of processes including emission from photo - dissociation regions , @xcite favoring x - ray heating or shocks , and so on .
recent papers have favored x - ray heating in the vicinity of an agn ( e.g. , * ? ? ?
however , in our case , the fact that there is little h@xmath4  coming the nlr while the bulk of the h@xmath4emission on these scales follows the morphology of the starlight ( figure [ fig : emissionflux ] ) strongly suggests that the h@xmath4  is predominantly excited by stellar processes ( e.g. , uv radiation and/or shocks from supernovas ) rather than the agn .
bridging the gap between the accretion flow on @xmath14pc scales and the large - scale ionized gas outflow , we can examine the ionized gas emission on 100 pc scales .
in this region , the radio emission is complex , extends predominantly north - south ( perpendicular to the disk ) and contains two unresolved peaks and a more diffuse `` bubble '' to the north @xcite .
@xcite show that the northern radio source is coincident with the optical nucleus , while the southern radio knot is coincident with a bright region of ionized gas seen in [ ] , strongly suggesting a jet - ism interaction .
high - ionization gas is present to the south of the nucleus in our data cubes as well .
previous near - infrared spectra ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) show extended h@xmath4  coincident with the galaxy disk , and high - ionization gas emission from br@xmath3  and [ ] coincident with the nucleus and ionization cone . with our high spatial resolution and three - dimensional coverage
, we can robustly disentangle the disk and nlr kinematics .
the br@xmath3  and [ si vi ]  emission lines show similar spatial distributions that parallel the nlr emission seen in [ ]  ( figures [ fig : emissionflux ] and [ fig : oiiisivi ] ) .
their emission is distributed at an angle of pa@xmath105 , tracing the general orientation of the radio jet seen on similar scales @xcite . the ionization cone
seen in [ ]  continues to the south - west beyond the sinfoni field - of - view @xcite , and is filled with radio plasma @xcite . to the north , beyond our field of view , the jet opens into another wide angle fan .
it is likely that the southern ionization cone sits in front of the galaxy disk , while the northern one ( which we can see in this map and in the radio but not in the optical ) is extincted by the galaxy disk .
[ tab : region ]    in figure [ fig:2dgas ] , we show the velocity and fwhm fields for br@xmath3  and [ si vi ] .
both the morphology and the kinematics differ from that seen in the h@xmath4 .
specifically , we do not see coherent rotation in either transition .
the skewness plots are also interesting . in both transitions we see redshifted gas towards the south - west and blue - shifted gas towards the north - east .
our interpretation is that we are seeing bulk motion of the gas along the jet , which apparently points towards us to the north and away from us to the south of the galaxy disk . in both cases
we are probing the region of the nlr where the jet has `` broken - out '' of the galaxy disk .
the higher velocity gas is moving along the walls of the bicone , as accelerated by the jet ( e.g. , * ? ? ?
* ; * ? ? ?
the coadded spectra from different regions reveal a similar story . in the disk regions ( regions 2 & 3 ) , where there is little emission from the nlr gas , we find strong and narrow h@xmath4 , and weak , but rather broad , br@xmath3  and [ ] . in the nlr regions ( 4 & 5 ) , br@xmath3  and [ si vi ]
are stronger and slightly broader than in the disk regions .
the line profiles are most complex at positions 6 , 7 , and 8 , regions that are also bright in [ ]  emission and where we see the strongest skewness in the lines ( fig .
[ fig : oiiisivi ] & fig . [ fig:2dgas ] ) .
we clearly see the red wing in the spectra from regions 6 & 8 , and the corresponding blue wing in the spectra from region 7 .     again in line with recent work ( e.g. , * ? ? ?
* ) , we find that the high - ionization lines and hydrogen recombination lines trace `` feedback '' .
that is , the emission in these lines is dominated by extraplanar gas likely excited by the agn , and perhaps entrained by the nuclear jet .
this nlr emission is also aligned roughly perpendicular to the megamaser disk on sub - pc scales .
finally , we ask whether there is a broad component to the br@xmath3 , since @xcite detect a weak broad base to the h@xmath65  line in ngc 4388 , with a fwhm of 3900  km  s@xmath35 . within the central nuclear point source emission ,
all the species have lines that are broad and symmetric ( table 1 ) .
however , we find no evidence for a true broad - line component in br@xmath3 . instead , all the transitions have comparable widths .
neither @xcite nor @xcite detect broad emission from br@xmath3  or br@xmath65 in their long - slit nir spectra either . perhaps the broad component detected in the optical
is scattered broad emission rather than directly transmitted light ( e.g. , * ? ? ?
ngc 1194 is an inclined ( @xmath133 ) s0 galaxy at @xmath134 mpc .
the galaxy pa is 145 , and we aligned the ifu along the major axis . given our psf of @xmath135 , we are resolving 40 pc scales with our ao - assisted observations .
thus we achieve very comparable spatial resolution to the ngc 4388 observations presented above .
we have detected the large - scale disk in ngc 1194 in , and find a circular velocity of @xmath136  km  s@xmath35  @xcite . from sdss imaging ,
we find a bulge - to - total ratio of @xmath137 @xcite .
the bh in ngc 1194 has a mass of @xmath138  @xmath62 , one of the most massive in the megamaser disk sample .
ngc 1194 is in the _ _
iras-__selected 12@xmath131 sample @xcite , and in the swift / bat 22 month sample @xcite .
it has a hard x - ray luminosity of log @xmath139 from 14 - 195 kev , but is not compton thick @xcite .
depending on the bolometric correction , the eddington ratio is @xmath140 @xcite , considerably lower than ngc 4388 .
the galaxy was also part of the _ hst_/ [ ]  snapshot survey of @xcite .
they find extended [ ]  emission with an extent of 700 pc along the major axis of the galaxy , and 470 pc along the minor axis .
the nucleus was detected in first and at 8.4 ghz by @xcite , with a flux density of 0.9 mjy , and an upper limit on the radio core size of 52 pc .
there are three existing @xmath0band spectra of ngc 1194 @xcite on @xmath90  ( 900 pc ) scales . in all cases
the co bandhead is detected , but the line emission is very weak .
davies et al
.  do not detect br@xmath3 , but marginally detect the h@xmath4@xmath141 line .          due to the heavy binning required to achieve adequate s / n
, we can not derive much spatially resolved information for this galaxy . in terms of the gas ,
the br@xmath3  and [ si vi ]  flux distributions are very centrally concentrated , with a fwhm that is consistent with that of the continuum .
the h@xmath4@xmath142 emission is clearly extended , but the distribution is quite patchy .
we therefore construct average spectra in annuli , taking the central spectrum within the fwhm of the psf ( 016 ) .
then we construct six spectra with outer radii of 05 , 1 , and 15 ( 140 , 280 , and 420 pc ) respectively , on either side of the galaxy minor axis .
three examples of these spectra are shown in figure [ fig : lineregions ] .
we can measure the radial velocities and stellar velocity dispersions reliably in these coadded spectra , using ppxf as above .
we do find that the stars are rotating , with an amplitude of @xmath2  km  s@xmath35(figure [ fig : n1194stelsigcurve ] ) .
in contrast , we do not detect a clear rotation curve in the gas .
instead , the velocities are very asymmetric .
asymmetric velocities in high - ionization lines are often attributed to outflow ( e.g. , * ? ? ?
* ; * ? ? ?
* ) , and are sometimes seen in h@xmath4  as well ( e.g. , * ? ? ?
* ) , which may be the explanation here .
another possibility is that we are seeing spiral arms or some other non - axisymmetric gas distribution that is unrelated to the agn .
our limited s / n prevents us from distinguishing between these cases .
neither the [ si vi ]  nor the br@xmath3  is spatially extended enough to measure reliable velocities beyond the central region .
all of the gas emission lines fluxes are weak , and the intrinsic luminosities are also low .
it is perhaps not surprising that this s0 galaxy is gas poor compared to ngc 4388 above , based simply on their respective morphology . on the other hand ,
the center of ngc 1194 is quite dusty , so there must be some associated gas , and indeed we detect atomic hydrogen in this galaxy on large scales @xcite .
the strongest gas emission comes from the high - ionization [ si vi ] line , but as mentioned above it is only marginally spatially resolved .
thus , we suspect that the emission from this transition and br@xmath3  emerge from the inner nlr .
in contrast , the h@xmath4  is clearly spatially resolved . in analogy with ngc 4388 above , these different gas morphologies
may be suggesting that the h@xmath4  is excited by both the agn and stellar processes .
we also measure the spatially resolved ratio of h@xmath4/br@xmath3 , and find that the ratio is lowest in the galaxy center ( having agn - like values ) and then rises outward , perhaps indicating an increased contribution to the h@xmath4  from stellar processes .
these line ratios may suggest that the asymmetric velocities we observe in the h@xmath4  are due to non - circular motions in the inner regions of the galaxy ( e.g. , spiral arms ) rather than outflow .
we have analyzed integral - field @xmath0band observations of two megamaser disk galaxies .
our data probe @xmath1 pc scales in both galaxies , and thus allow us to study the distributions of stars and gas at the centers of these galaxies to investigate agn fueling . in ngc 1194 , an s0 , there is very little gas , and with our s / n ratios , we can not say much about the two - dimensional velocity or dispersion fields .
however , in ngc 4388 we have excellent s / n , and uncover a variety of interesting features :    1 .
the stellar velocity field demonstrates well - ordered rotation aligned with the kpc - scale disk .
although there is net rotation in the stars , still they are kinematically hot , with @xmath79 , likely dominated by the dispersion in the kpc - scale bulge .
2 .   in the inner 1  ( 100 pc )
, we also see evidence for a nuclear disk , offet in pa by 15  from the kpc - scale disk .
the evidence includes disky isophotes in the _
hst_/wfc3 f160w image and a distinct drop in stellar velocity dispersion on one side of the putative disk ( the other side apparently extincted ) .
features like this , so called @xmath84drops , have been seen in the centers of a number of nearby spirals ( e.g. , * ? ? ?
* ; * ? ? ?
* ) , and there is some hint that they are more common among actively accreting galaxies @xcite .
we note that this 100 pc - scale nuclear disk is misaligned by @xmath1135  from the megamaser disk on sub - pc scales .
the h@xmath4  gas also shows regular kinematics on @xmath143 pc scales , with a well - defined rotation field that is aligned with the kpc - scale disk . with a higher rotation amplitude ,
the gas is considerably colder than the stars ( @xmath144 ) .
we also see a kink in the rotation curve on @xmath2 pc scales .
intriguingly , the h@xmath4  rotation within 100 pc appears to align with the megamaser disk .
however , the physical origin of the kinematic twist is not yet clear .
4 .   based on h@xmath4  diagnostic line ratios ,
we conclude that the molecular hydrogen is mainly excited by thermal processes . based on the morphology and kinematics of the h@xmath4  gas
, we speculate that these thermal processes are mostly associated with stars rather than radiation from the agn ( e.g. , * ? ? ?
* ; * ? ? ?
* ) . 5 .
in contrast , the br@xmath3  and high - ionization lines ( particularly [ si vi ] ) have a completely different spatial distribution and kinematics .
they trace the inner narrow - line region , also seen on 500 pc scales in [ ] , as well as the 100 pc - scale jet .    in the future ,
we hope to analyze a larger sample of megamaser disk galaxies in this manner . combined with direct tracers of the cold molecular gas ( e.g. , with alma ) , as well as high - resolution imaging from _ hst _
, we hope to build a multi - phase map of agn fueling from 100 to 0.1 pc @xcite .
the referee provided a very thorough review that improved the quality of this manuscript .
we thank k. gebhardt and r. van den bosch for assistance with the early stages of this project .
j.l.w . has been supported by an nsf astronomy and astrophysics postdoctoral fellowship under award no .
some of the data presented herein were obtained at the w.m .
keck observatory , which is operated as a scientific partnership among the california institute of technology , the university of california and the national aeronautics and space administration .
the observatory was made possible by the generous financial support of the w.m .
keck foundation .
the authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of mauna kea has always had within the indigenous hawaiian community .
we are most fortunate to have the opportunity to conduct observations from this mountain .
keck telescope time was granted by noao , through the telescope system instrumentation program ( tsip ) .
tsip is funded by nsf .


the standard model ( sm ) has been tested at an unprecedented level of precision by various experiments . however , some excesses have not yet been completely resolved .
the first case is the muon anomalous magnetic moment ( muon @xmath3 ) , where the discrepancy between experimental data and the sm prediction currently is @xmath8  @xcite .
the second case is the angular observable @xmath1 of @xmath2  @xcite , where @xmath9 deviation , resulting from an integrated luminosity of 3.0 fb@xmath10 at the lhcb  @xcite , recently confirmed an earlier result with @xmath11 deviations  @xcite .
moreover , the same measurement with @xmath12 deviations was reported by belle  @xcite .
also the other relevant @xmath13 observables are defined in ref .
various possible resolutions to this excess have been widely studied  @xcite .
the third excess is that a resonance with a mass of 750 gev in the diphoton invariant mass spectrum was indicated by the atlas  @xcite and cms  @xcite experiments .
the width and spin of the diphoton resonance are still uncertain .
if the resonance is a fundamental particle , its width prefers to be narrow ; however , if it is a composite state , then its width is wide .
the production cross sections for the diphoton excess at @xmath14 tev in the atlas and cms are respectively given by : @xmath15 where @xmath16 stands for the diphoton resonance .
the challenge of explaining this excess is to obtain the large production cross section shown in eq .
( [ eq : data_di ] ) .
various proposals to explain this excess have been broadly studied  @xcite .    to resolve the excesses in a specific framework
, we propose the extension of the sm by including leptoquarks ( lqs ) and one higgs singlet , where the leptoquarks are colored scalars that simultaneously couple to the leptons and quarks , and the higgs singlet is the 750 gev diphoton resonance candidate .
therefore , the @xmath17 decays can arise from the tree - level lq - mediated feynman diagrams when the muon @xmath3 is induced from lq loops .
although lepton - flavor violating ( lfv ) processes @xmath4 can put strong constraints on some yukawa couplings , it is found that the wilson coefficient @xmath18 from the lq contributions can satisfy the requirements , which result from the global analysis of @xmath2 , and the muon @xmath3 can fit the current data .
in addition , the results are also consistent with the measurement of the rare flavor - changing process @xmath19 .
it is known that the production cross section of eq .
( [ eq : data_di ] ) can not be achieved if only a 750 gev particle is added to the sm without introducing other new effects .
therefore , we need not only some new colored particles to produce the resonance , but also some mechanism by which to enhance the branching ratio ( br ) for the resonance decay to diphoton . intriguingly , lqs meet the requirements for the resonance production and diphoton decay .
some leptoquark explanations of the diphoton excess can be found in refs .  @xcite    the paper is organized as follows :
we introduce the model and discuss the relevant couplings in sec .
ii . in sec .
iii , we study the phenomenon : the sm higgs diphoton decay , lfv processes , wilson coefficients of @xmath20 from lq contributions for @xmath21 decays , the implication of @xmath19 , and the production cross section for @xmath22 in the model .
the conclusion is given in sec .
in this section , we briefly introduce the model and relevant interactions with the lqs . to reconcile the measurements of @xmath23 and @xmath24
, we extend the sm by adding two different representations of lq , which are @xmath25 and @xmath26 under @xmath27 gauge symmetry . in addition
, we also add one real higgs singlet @xmath28 as the candidate for 750 gev diphoton resonance .
the gauge - invariant yukawa interactions of the sm fermions and lqs are written as : @xmath29 where the subscripts @xmath30 are the flavor indices ; @xmath31 and @xmath32 are the @xmath33 lepton and quark doublets , @xmath34 , and @xmath35 are the yukawa couplings . since we do not study the cp violating effects , hereafter , we take all new yukawa couplings as real numbers .
we use the representations of the lqs as : @xmath36 where the superscripts are the electric charges of the particles . the interactions in eq .
( [ eq : int1 ] ) are then expressed as : @xmath37 +   \tilde{k}_{ij } \left [ \bar \ell_{l i }   \ , u_{rj } \phi^{-5/3 } - \bar \nu_{l i}\ , u_{r j } \phi^{-2/3 } \right ]    \nonumber \\   & + y_{ij } \left [ \bar u^c_{li }   \ , \nu_{lj } \delta^{-2/3 } - \frac{1}{\sqrt{2 } } \bar u^c_{l i }   \ , \ell_{l j } \delta^{1/3 } - \frac{1}{\sqrt{2 } } \bar d^c_{l i } \ , \nu_{l j } \delta^{1/3 } - \bar d^c_{l i }   \ , \ell_{l j } \delta^{4/3 } \right ] + h.c . \end{aligned}\ ] ]    since the lqs are colored scalar bosons , they can couple to the sm higgs @xmath38 and higgs singlet @xmath28 via the scalar potential . in order to get the higgs and @xmath28 couplings to the lqs
, we write the gauge - invariant scalar potential as : @xmath39 ^ 2 + \lambda'_\delta tr \left ( \left [ \delta_{1/3}^\dagger \delta_{1/3 } \right]^2 \right ) \nonumber \\ & + \lambda_{h \phi } ( h^\dagger h)\left(\phi_{7/6}^\dagger \phi_{7/6 } \right ) + \lambda_{h \delta }   ( h^\dagger h ) tr\left ( \delta_{1/3}^\dagger \delta_{1/3 } \right ) + \lambda_{\phi \delta } \left(\phi_{7/6}^\dagger \phi_{7/6 } \right )   tr \left ( \delta_{1/3}^\dagger \delta_{1/3 } \right ) , \nonumber \\ & + \mu_{sh } s ( h^\dagger h ) + \mu_{s \phi } s   \left(\phi_{7/6}^\dagger \phi_{7/6 } \right ) + \mu_{s \delta } s \ , tr \left ( \delta_{1/3}^\dagger \delta_{1/3 } \right ) \nonumber \\ & + \lambda_{sh } s^2   ( h^\dagger h ) + \lambda_{s \phi } s^2   \left(\phi_{7/6}^\dagger \phi_{7/6 } \right ) + \lambda_{s \delta } s^2 \ , tr \left ( \delta_{1/3}^\dagger \delta_{1/3 } \right ) \ , .
\end{aligned}\ ] ] as usual , we adopt the representations of the higgs doublet @xmath38 and singlet @xmath28 as : @xmath40 where @xmath41 and @xmath42 are the goldstone bosons ; @xmath43 is the sm higgs field , and @xmath44 is the vacuum expectation value ( vev ) of @xmath38 ( @xmath28 ) .
it is known that the vevs of scalar fields are dictated by the scalar potential .
thus , by the minimal conditions @xmath45 and @xmath46 , the stable vevs are determined by : @xmath47    using the scalar potential of eq .
( [ eq : potential ] ) , the mass - square matrix of @xmath43 and @xmath48 is found by : (    cc , & s    ) (    cc m^2 _ & _ sh v/ + _ sh v v_s + _ sh v/ + _ sh v v_s & m^2_s + _ sh v^2 /2    ) (    c + s    ) [ eq : m2 ] with @xmath49 .
the parameters @xmath50 and @xmath51 lead to the mixture of @xmath43 and @xmath48 . since two neutral scalar bosons are involved in the model , we only need one mixing angle to parametrize the mixing effect .
we formulate the mass eigenstates to be (    c h + h_s    ) = (    cc & - + &    ) (    c + s    ) , [ eq : mixing ] where @xmath52 is the sm - like higgs boson , @xmath53 is the new heavy higgs boson and the new resonance , and their masses are obtained as m^2_h(h_s ) = [ eq : hmass ] with @xmath54 and @xmath55 .
the relationship of the @xmath56 and mixing angle can be expressed by @xmath57 .
since it is known that the higgs diphoton decay measurement has been consistent with the sm prediction , the data from the higgs to diphoton decay can therefore give strict limits on the new physics effects . in order to suppress the contribution from the mixing of @xmath52 and @xmath53
, we set @xmath58 ; thus , the @xmath59 can be simplified as @xmath60 . with @xmath61 gev , @xmath62 gev , and @xmath63 gev
, it can be seen that @xmath64 gev and @xmath65 .
based on the scheme of a small @xmath66 , the contributions from the mixing of @xmath52 and @xmath53 can be ignored .
based on the introduced new interactions , in this section , we study the implications of the higgs diphoton decay , @xmath24 , @xmath4 , the muon @xmath3 , @xmath67 , and the 750 gev diphoton excess . since each of these processes has its own unique characteristics , we discuss these phenomena one by one below .
the higgs measurement is usually described by the signal strength parameter , which is defined as the ratio of observation to the sm prediction and is expressed as : ^f_i = _ i _ f , [ eq : muf ] where @xmath68 stands for the possible channels , and @xmath69 denotes the signal strength of production ( decay ) .
although vector - boson fusion ( vbf ) can also produce the sm higgs , we only consider the gluon - gluon fusion ( ggf ) process because it is the most dominant .
the diphoton higgs decay approached the precision measurement since the 125 gev higgs was observed . therefore , any significant deviation from the sm prediction ( i.e. , @xmath70 ) can imply the new physics effects .
as stated earlier , the sm higgs can couple to the lqs via the scalar potential . from eq .
( [ eq : potential ] ) , it can be seen that after spontaneous symmetry breaking ( ssb ) , the quartic terms @xmath71 and @xmath72 can lead to trilinear couplings of higgs to lqs as : @xmath73 where @xmath74 and @xmath75 . with the couplings in eq .
( [ eq : hlq ] ) , the effective lagrangian for @xmath76 by lq - loop can be formulated as : @xmath77 where @xmath78 and the loop function is given by : @xmath79 with @xmath80 ^ 2 $ ] for @xmath81 .
accordingly , the signal strength of the higgs production and decay to diphoton can be respectively obtained as : @xmath82 where @xmath83 , @xmath84 is the number of colors ; @xmath85 and @xmath86 , and the functions for vector - boson and fermion loops are given by @xmath87\ , , \nonumber   \\
a_1(x )   & = 2 + 3 x + 3(2x -x^2)f(x)\,.\end{aligned}\ ] ]    since the effects of the doublet and triplet lqs are similar , for simplicity , we set @xmath88 and @xmath89 .
the @xmath90 as a function of @xmath91 is presented in fig .
[ fig : muif](a ) , and that of @xmath92 is shown in fig .
[ fig : muif](b ) , where the curves in plot ( a ) are @xmath93 tev , and those in plot ( b ) are @xmath94 tev . for comparison
, we also show the results of atlas  @xcite and cms  @xcite with @xmath95 errors in the plots . from the plots
, it can be clearly seen that with @xmath92 of @xmath96 gev , the lq contributions can significantly shift the @xmath90 away from the sm prediction and that the results are consistent with the current data . on the contrary ,
the @xmath90 approaches the sm result when @xmath92 is of the order of gev .     as a function of ( a ) @xmath91 and ( b ) @xmath92 , where the curves in plots ( a ) and ( b ) denote @xmath97 tev and @xmath98 tev , respectively .
, title="fig:",width=264 ]   as a function of ( a ) @xmath91 and ( b ) @xmath92 , where the curves in plots ( a ) and ( b ) denote @xmath97 tev and @xmath98 tev , respectively . , title="fig:",width=264 ]      in the following analysis , we study the rare lepton - flavor violating processes , e.g. @xmath4 and @xmath100 , muon @xmath3 @xmath101 , and the fcnc process @xmath24 .
we first discuss the radiative lfv processes for @xmath4 . with the couplings in eq .
( [ eq : int2 ] ) , the lq - loop induced decay amplitude for @xmath4 can be written as : @xmath102 \ell_i f^{\mu \nu } \,,\label{eq : lllga }    \end{aligned}\ ] ] where the coefficient @xmath103 is expressed as : @xmath104     \left (   \frac{5}{\delta(m_{t } , m_\phi ) } - \frac{2(1-x)}{\delta(m_\phi , m_{t } ) } \right ) \
, , \nonumber \\
\delta(m_1,m_2 )   & =     x m_1 ^ 2 + ( y+z ) m_2 ^ 2 \ , , \nonumber \\ \int [ dx ] & =   \int dx dy dz \delta(1-x - y - z)\ , , \label{eq : crji}\end{aligned}\ ] ] @xmath105 can be obtained from @xmath106 by exchanging @xmath107 and @xmath108 . in order to balance the chirality of the leptons ,
it is found that the contributions from @xmath109 , @xmath110 , @xmath111 , and @xmath112 are suppressed by the lepton masses .
since the lq @xmath113 can couple to left - handed and right - handed up - type quarks , the chirality flip by the mass insertion in the propagator of the up - type quark can lead to freeing of the lepton masses in the feynman diagrams , which are associated with @xmath114 and @xmath115 .
in addition , the top - quark is much heavier than the @xmath116- and @xmath117-quarks ; therefore , we only present the top - quark contribution in @xmath106 .
straightforwardly , the br for @xmath4 can be expressed as : @xmath118 where @xmath119 for @xmath120 and the brs for @xmath121 in the sm have been applied .
the current experimental upper limits are shown in table  [ tab : cif ] . according to eq .
( [ eq : lllga ] ) , muon @xmath3 can be easily obtained by setting @xmath122 and found as : @xmath123    .
current upper bounds on the brs for the decays @xmath4  @xcite .
[ cols="^,^,^",options="header " , ]     if the photon in @xmath4 is replaced by the higgs , similar feynman diagrams can contribute to @xmath124
. since the upper limit of @xmath125 is of @xmath126 and can give strong constraints on the parameters @xmath127 and @xmath128 , if we set @xmath129 and @xmath130 to be small , then it is apparent that @xmath131 and @xmath132 are much smaller than current upper limits .
hence , we just study the decay @xmath133 .
the one - loop induced effective couplings for @xmath134 are written as : @xmath135 where @xmath136 is expressed as  @xcite : @xmath137 \nonumber \\ & +   \frac { n_c \mu_{h \phi } } { ( 4\pi)^2 } \sum_{i=1\text{-}3 , q = u , d } \left [ m_\mu ( k^\dagger)_{2 i } k_{i 3 } g(m_{q_i},m_\phi ) + m_\tau \tilde k^\dagger_{2 i } \tilde k_{i 3 } \tilde g(m_{q_i},m_\phi ) \right ]   \nonumber \\ & + \frac{n_c m_\tau \mu_{h \delta}}{(4 \pi)^2 } \sum_{i=1\text{-}3,q = u , d } ( y^\dagger)_{2 i } y_{i 3 }   \tilde g(m_{q_i},m_\delta ) \ , ,   \label{eq : cl}\end{aligned}\ ] ] @xmath138 can be obtained from @xmath139 by exchanging @xmath107 and @xmath108 , and the loop functions are given by : @xmath140 \log \left[z+(1-z)r_t - x y
r_h - i \epsilon \right ] \nonumber \\ & + \int_0 ^ 1 dx \log \left [ x + ( 1-x ) r_t - i \epsilon \right]\ , , \nonumber   \\ b(r_t , r_h ) = & \int [ dx ] \frac{xy r_h - r_t } { z+(1-z)r_t - xy r_h } \
, , \nonumber \\
g(m_1,m_2 ) \approx   & \int [ dx ] \frac { z}{- x z m_h^2 + x m_1 ^ 2 + ( y+z ) m_2 ^ 2 } \ , , \nonumber \\
\tilde g(m_1,m_2 ) \approx & \int [ dx ] \frac { y } { - x z m_h^2 + x m_1 ^ 2 + ( y+z ) m_2 ^ 2 } \,.\end{aligned}\ ] ] the @xmath141 in @xmath142 denotes an infinitesimal positive value .
it can be seen that the terms associated with @xmath143 , @xmath144 , and @xmath145 in eq .
( [ eq : cl ] ) are proportional to the lepton masses .
the situation is similar to the @xmath106 in the decays @xmath4 .
although @xmath146 of tev ( @xmath147 can enhance these effects , due to the effects of being related to @xmath148 , their contributions are at least @xmath149 smaller than those from @xmath150 . accordingly , the br for @xmath151 is formulated as : @xmath152 where @xmath153 is the width of the higgs . due to @xmath154 being less than @xmath155 , we use @xmath156 mev in our numerical estimations .
next , we discuss the decays for @xmath17 . in this study
, we do not directly investigate the angular analysis of @xmath24 ; instead , we refer to the results , which were done by using the global analysis to get the best - fit value of @xmath157 for the new physics contributions  @xcite . thus
, we just derive the wilson coefficients of @xmath18 and @xmath158 from the lq contributions . with the yukawa couplings in eq .
( [ eq : int2 ] ) , the effective hamiltonian for @xmath17 mediated by @xmath159 and @xmath160 can be respectively found as : @xmath161 we can decompose the eq .
( [ eq : heff ] ) in terms of the effective operators @xmath162 and @xmath163 , defined as @xmath164 .
the associated wilson coefficients of @xmath165 from the lqs then are found as : @xmath166 where @xmath167 is a scale factor from the sm effective hamiltonian .
it is worth mentioning that the interaction @xmath168 can contribute to @xmath19 , where the brs , measured by lhcb and cms  @xcite and predicted by the sm  @xcite , are : @xmath169 to consider the constraint from @xmath19 , we adopt the expression for the br as  @xcite : @xmath170 with @xmath95 errors , the allowed range for @xmath171 is obtained as @xmath172 .
we use this result to constrain the free parameters .
since the parameters in the decays @xmath4 , @xmath133 , @xmath101 , and @xmath24 are strongly correlated , in the following analysis , we take the current upper limits of @xmath173 shown in table  [ tab : cif ] as the inputs and attempt to find the allowed parameter space , such that the excesses in @xmath101 and @xmath24 can be satisfied , and the @xmath174 can be as large as possible .    from @xmath106 in eq .
( [ eq : crji ] ) , the dominant effects on the radiative lfv processes are from the @xmath113 and the top - quark loop ; thus , there is no possible cancellation in any of the decay amplitudes . with the upper bound of @xmath175
, we see that @xmath176 and @xmath177 have to be very small . in order to explain the excesses of muon @xmath3 and @xmath2 , we set @xmath178 . as a result
, @xmath179 is negligible in this model .
the related parameters for @xmath180 decays are @xmath181 and @xmath182 , respectively .
these parameters simultaneously influence @xmath183 , muon @xmath3 , and @xmath184 ; therefore we have to analyze these processes together to get the allowed parameter space .    since eqs .
( [ eq : crji ] ) , ( [ eq : cl ] ) , and ( [ eq : c9c10 ] ) involve many free parameters , in order to efficiently perform a numerical analysis , we set the ranges of relevant parameters as : @xmath185 \ ; \text{gev},\quad   \mu_{lq } \in [ 1\ , , 100]\ ; \text{gev}\ , , \quad \ { k_{22 } , \tilde k_{22 } , y_{22 } \ }   \in [ -0.1\ , , 0.1]\,,\nonumber \\ &   \ { k_{33},\tilde k_{33 } , y_{33 } \ }   \in [ -0.01,0.01]\ , , \quad \ { k_{23},\tilde k_{23 } , y_{23 } \ }   \in [ -0.1\ , , 0.1]\ , ,   \nonumber \\ & k_{32 } \in \text{sign}(k_{22 } ) [ 0,0.5 ] \ , , \quad \tilde k_{32 } \in [ -0.5\ , , 0.5]\ , , \quad y_{32 } \in -\text{sign}(y_{22 } ) [ 0\ , , 0.5 ] \ , .
\label{eq : range_scanning}\end{aligned}\ ] ] in order to avoid the constraints from @xmath186 ( @xmath187 ) and get @xmath188 , we set ( @xmath189 in eq .
( [ eq : range_scanning ] ) .
additionally , the negative value of @xmath190 can be achieved when @xmath191 and @xmath192 are opposite in sign .
as mentioned earlier , the yukawa couplings in both decays @xmath193 and @xmath194 are the same , we can not remove the constraints from the radiative lfv processes in this model .
the brs for @xmath194 thus are of @xmath195 and much smaller than the current upper limits of @xmath196  @xcite .
one way to escape the constraint from @xmath193 is to add a new lq  @xcite .
since we focus on the excesses of muon @xmath3 and @xmath2 , we leave the more complicated model for further study .    with the chosen ranges of parameters in eq .
( [ eq : range_scanning ] ) , we first show the values for @xmath197 and @xmath171 in fig .
[ fig : scatter](a ) , where the bounds from @xmath193 have been considered ; the horizontal band is from the measurement of @xmath5 ; the vertical band is the range that can explain the excess of @xmath2 , and we used @xmath198 parameter sets and obtained @xmath199 allowed points that satisfy the constraints . it can be seen that the @xmath197 and @xmath171 from the contributions of doublet @xmath200 and triplet @xmath201 lqs can simultaneously satisfy the constraint of @xmath5 and explain the excess in @xmath2 . from eq .
( [ eq : mug2 ] ) , it is known that muon @xmath3 is associated with the yukawa couplings @xmath202 .
although only @xmath203 is related to @xmath197 and @xmath171 , since the yukawa couplings @xmath204 , @xmath205 , and @xmath206 are taken to be the same order of magnitude , we present the correlations of @xmath101 and @xmath190 in fig .  [
fig : scatter](b ) , where only the allowed range of @xmath190 is shown . by plot ( b ) , it can be seen clearly that the excesses in @xmath101 and @xmath2 can be simultaneously fitted in the model .     and @xmath171 using the ranges of parameters in eq .
( [ eq : range_scanning ] ) , where the bound from the @xmath207 and the allowed range of @xmath197 from global analysis for @xmath2 are shown .
( b ) correlation of @xmath101 and @xmath197 , where we only show the values of @xmath197 that can fit the excess in @xmath2 .
[ fig : scatter],title="fig:",width=264 ]   and @xmath171 using the ranges of parameters in eq .
( [ eq : range_scanning ] ) , where the bound from the @xmath207 and the allowed range of @xmath197 from global analysis for @xmath2 are shown .
( b ) correlation of @xmath101 and @xmath197 , where we only show the values of @xmath197 that can fit the excess in @xmath2 .
[ fig : scatter],title="fig:",width=264 ]      as discussed before , the 750 gev scalar boson candidate is a singlet , and it can mix with the sm higgs doublet via the scalar potential .
since the mixing angle of @xmath52 and @xmath53 is small in our scheme , we use @xmath28 instead of the physical state @xmath53 and ignore the mixing effect @xmath208 .
the relevant interactions for the @xmath28 production and decays are from the trilinear and quartic couplings , i.e. , @xmath209 , @xmath210 .
it was discussed earlier that the vev of @xmath28 is smaller than 1 gev , so we thus focus on the trilinear couplings in the numerical analysis .    similar to eq .
( [ eq : hgg ] ) , which is for @xmath76 effective coupling , using the trilinear interactions of @xmath28 and lqs in eq .
( [ eq : potential ] ) and the gauge couplings of gluons and lqs , the loop induced effective coupling for @xmath211 is found as : @xmath212 where @xmath213 and @xmath214 are the trilinear couplings with mass dimension - one .
the partial decay widths for @xmath215 decays via the lq loops are expressed as : @xmath216 due to @xmath217 , it can be seen that @xmath218 is the dominant channel .
although @xmath28 can decay to @xmath219 , their partial decay widths are the same as that of @xmath220 in terms of order of magnitude .
additionally , @xmath221 can be also induced by lq loops , however the decay becomes unimportant when @xmath92 is small .
hence , the br for @xmath222 can be estimated by : @xmath223 it is found that with @xmath224 and @xmath225 , the br for diphoton decay is insensitive to the parameters and becomes @xmath226 .    to calculate the production cross section of @xmath227 using the lq loops , we implement the model to the calchep  @xcite and use the cteq6l parton distribution function ( pdf )  @xcite
thus , the @xmath228 at @xmath14 tev as a function of @xmath91 is presented in fig .
[ fig : xsbr ] , where we use the k - factor @xmath229 for the purpose of including next - to - leading order qcd corrections  @xcite , and the solid , dashed , and dotted lines are for @xmath230 tev , respectively . from the results
, it can be seen that by taking proper values of @xmath231 , @xmath232 fb can be achieved with @xmath7 tev .     at @xmath14 tev as a function of @xmath91 , where the @xmath229 for next - to - leading order qcd corrections is used , and the solid , dashed , and dotted lines are for @xmath233 tev , respectively.,width=264 ]
in order to resolve the excesses of muon @xmath3 and @xmath2 decay , we investigate the extension of the sm by including leptoquarks , in which the particles are colored scalar bosons and can couple to quarks and leptons . in order to accommodate the measurement of @xmath5 and the excess of @xmath2 , we study a model with both doublet and triplet leptoquarks . moreover , we add a real singlet field to the model as the candidate of the 750 gev diphoton resonance , which is indicated by the atlas and cms experiments .
after ssb , the couplings of the sm higgs the lqs are described by @xmath234 .
if @xmath146 is of @xmath96 gev , the signal strength parameter @xmath235 can significantly deviate from the sm prediction and is still consistent with the current higgs measurements .
in this study , lepton - flavor violating processes @xmath4 give strict constraints on the yukawa couplings @xmath236 and @xmath237 .
as a result , the branching ratios for the lepton - flavor violating higgs decays @xmath238 are less than @xmath239 ) . nevertheless , the sizable couplings @xmath240 , @xmath241 , and @xmath242 can still explain the excess of muon @xmath3 and provide the necessary values for the wilson coefficient @xmath18 , such that the excess in @xmath2 can be resolved .
the 750 gev scalar boson @xmath28 can couple to leptoquarks through the trilinear couplings @xmath243 and @xmath244 . in light of these couplings , the @xmath28 production and decays in the @xmath245 collisions can occur through the leptoquark loops . with the @xmath246 of tev , the production cross section for
@xmath247 can reach 5 fb , which fits the current atlas and cms results .
we remark on the production of the leptoquarks at the collider .
the leptoquark pairs can be produced via qcd channels . since the masses of leptoquarks used for the purpose of explaining the excesses are @xmath6 tev
, their production cross sections can be sizable .
thus , our model can be tested by searching for specific signatures , such as @xmath248 , @xmath249 , @xmath250 , and @xmath251 that depend on the electric charges of leptoquarks and their yukawa couplings .
the detailed simulation is beyond the scope of this paper and the relevant study will be presented elsewhere .
this work was partially supported by the ministry of science and technology of taiwan r.o.c .
, under grant most-103 - 2112-m-006 - 004-my3 ( chc ) .
h. o. would like to thank the members of kias for their hospitality during his visit .
a.  abdesselam _ et al . _ [ belle collaboration ] , arxiv:1604.04042 [ hep - ex ] .
j.  matias , f.  mescia , m.  ramon and j.  virto , jhep * 1204 * , 104 ( 2012 ) [ arxiv:1202.4266 [ hep - ph ] ] . s.  descotes - genon , j.  matias and j.  virto , phys .
d * 88 * , 074002 ( 2013 ) [ arxiv:1307.5683 [ hep - ph ] ] .
r.  gauld , f.  goertz and u.  haisch , jhep * 1401 * , 069 ( 2014 ) [ arxiv:1310.1082 [ hep - ph ] ] .
a.  datta , m.  duraisamy and d.  ghosh , phys .
d * 89 * , no . 7 , 071501 ( 2014 ) [ arxiv:1310.1937 [ hep - ph ] ] .
t.  hurth and f.  mahmoudi , jhep * 1404 * , 097 ( 2014 ) [ arxiv:1312.5267 [ hep - ph ] ] .    s.  descotes - genon , l.  hofer , j.  matias and j.  virto , jhep * 1412 * , 125 ( 2014 ) [ arxiv:1407.8526 [ hep - ph ] ] . w.  altmannshofer and d.  m.  straub , eur .
j.  c * 75 * , no . 8 , 382 ( 2015 ) [ arxiv:1411.3161 [ hep - ph ] ] .
s.  descotes - genon , l.  hofer , j.  matias and j.  virto , jhep * 1606 * , 092 ( 2016 ) [ arxiv:1510.04239 [ hep - ph ] ] .
b.  gripaios , m.  nardecchia and s.  a.  renner , jhep * 1505 * , 006 ( 2015 ) [ arxiv:1412.1791 [ hep - ph ] ] .
a.  crivellin , g.  dambrosio and j.  heeck , phys .
lett .   * 114 * , 151801 ( 2015 ) [ arxiv:1501.00993 [ hep - ph ] ] . s.  sahoo and r.  mohanta , phys .
d * 91 * , no .
9 , 094019 ( 2015 ) [ arxiv:1501.05193 [ hep - ph ] ] .
a.  bharucha , d.  m.  straub and r.  zwicky , arxiv:1503.05534 [ hep - ph ] .
d.  becirevic , s.  fajfer and n.  kosnik , phys .
d * 92 * , no . 1 , 014016 ( 2015 ) [ arxiv:1503.09024 [ hep - ph ] ] .
a.  crivellin , l.  hofer , j.  matias , u.  nierste , s.  pokorski and j.  rosiek , phys .
d * 92 * , no . 5 , 054013 ( 2015 ) [ arxiv:1504.07928 [ hep - ph ] ] .
c.  j.  lee and j.  tandean , jhep * 1508 * , 123 ( 2015 ) [ arxiv:1505.04692 [ hep - ph ] ] .
r.  alonso , b.  grinstein and j.  martin camalich , jhep * 1510 * , 184 ( 2015 ) [ arxiv:1505.05164 [ hep - ph ] ] .
s.  sahoo and r.  mohanta , phys .
d * 93 * , no . 3 , 034018 ( 2016 ) [ arxiv:1507.02070 [ hep - ph ] ] . g.  belanger , c.  delaunay and s.  westhoff , phys .  rev .
d * 92 * , 055021 ( 2015 ) [ arxiv:1507.06660 [ hep - ph ] ] .
s.  sahoo and r.  mohanta , phys .
d * 93 * , no .
11 , 114001 ( 2016 ) [ arxiv:1512.04657 [ hep - ph ] ] . c.  w.  chiang , x.  g.  he and g.  valencia , phys .
d * 93 * , no . 7 , 074003 ( 2016 ) [ arxiv:1601.07328 [ hep - ph ] ] .
i.  dorsner , s.  fajfer , a.  greljo , j.  f.  kamenik and n.  kosnik , phys .
rept .   * 641 * , 1 ( 2016 ) [ arxiv:1603.04993 [ hep - ph ] ] .
s.  m.  boucenna , a.  celis , j.  fuentes - martin , a.  vicente and j.  virto , phys .
b * 760 * , 214 ( 2016 ) [ arxiv:1604.03088 [ hep - ph ] ] .
k.  harigaya and y.  nomura , phys .
b * 754 * , 151 ( 2016 ) [ arxiv:1512.04850 [ hep - ph ] ] . y.  mambrini , g.  arcadi and a.  djouadi , phys .
b * 755 * , 426 ( 2016 ) [ arxiv:1512.04913 [ hep - ph ] ] .
m.  backovic , a.  mariotti and d.  redigolo , jhep * 1603 * , 157 ( 2016 ) [ arxiv:1512.04917 [ hep - ph ] ] .
a.  angelescu , a.  djouadi and g.  moreau , phys .
b * 756 * , 126 ( 2016 ) [ arxiv:1512.04921 [ hep - ph ] ] .
y.  nakai , r.  sato and k.  tobioka , phys .
* 116 * , no .
15 , 151802 ( 2016 ) [ arxiv:1512.04924 [ hep - ph ] ] .
d.  buttazzo , a.  greljo and d.  marzocca , eur .
j.  c * 76 * , no . 3 , 116 ( 2016 ) [ arxiv:1512.04929 [ hep - ph ] ] .
s.  di chiara , l.  marzola and m.  raidal , phys .
d * 93 * , no .
9 , 095018 ( 2016 ) [ arxiv:1512.04939 [ hep - ph ] ] .
s.  knapen , t.  melia , m.  papucci and k.  zurek , phys .
d * 93 * , no . 7 , 075020 ( 2016 ) [ arxiv:1512.04928 [ hep - ph ] ] .
a.  pilaftsis , phys .
d * 93 * , no . 1 , 015017 ( 2016 ) [ arxiv:1512.04931 [ hep - ph ] ] .
r.  franceschini _ et al .
_ , jhep * 1603 * , 144 ( 2016 ) [ arxiv:1512.04933 [ hep - ph ] ]
. j.  ellis , s.  a.  r.  ellis , j.  quevillon , v.  sanz and t.  you , jhep * 1603 * , 176 ( 2016 ) [ arxiv:1512.05327 [ hep - ph ] ] .
r.  s.  gupta , s.  jager , y.  kats , g.  perez and e.  stamou , arxiv:1512.05332 [ hep - ph ] .
a.  kobakhidze , f.  wang , l.  wu , j.  m.  yang and m.  zhang , phys .
b * 757 * ( 2016 ) 92 [ arxiv:1512.05585 [ hep - ph ] ] .
a.  falkowski , o.  slone and t.  volansky , jhep * 1602 * , 152 ( 2016 ) [ arxiv:1512.05777 [ hep - ph ] ] .
r.  benbrik , c.  h.  chen and t.  nomura , phys .  rev .
d * 93 * , no . 5 , 055034 ( 2016 ) [ arxiv:1512.06028 [ hep - ph ] ] . f.  wang , l.  wu , j.  m.  yang and m.  zhang , phys .
b * 759 * ( 2016 ) 191 [ arxiv:1512.06715 [ hep - ph ] ] .    p.  s.  b.  dev and d.  teresi , arxiv:1512.07243 [ hep - ph ] .
b.  c.  allanach , p.  s.  b.  dev , s.  a.  renner and k.  sakurai , phys .
d * 93 * , no .
11 , 115022 ( 2016 ) [ arxiv:1512.07645 [ hep - ph ] ]
. k.  cheung , p.  ko , j.  s.  lee , j.  park and p.  y.  tseng , arxiv:1512.07853 [ hep - ph ] .
f.  wang , w.  wang , l.  wu , j.  m.  yang and m.  zhang , arxiv:1512.08434 [ hep - ph ] . c.  w.  chiang , m.  ibe and t.  t.  yanagida , jhep * 1605 * , 084 ( 2016 ) [ arxiv:1512.08895 [ hep - ph ] ] .
x.  j.  huang , w.  h.  zhang and y.  f.  zhou , phys .
d * 93 * , 115006 ( 2016 ) [ arxiv:1512.08992 [ hep - ph ] ] .
s.  kanemura , k.  nishiwaki , h.  okada , y.  orikasa , s.  c.  park and r.  watanabe , arxiv:1512.09048 [ hep - ph ] . t.  nomura and h.  okada , phys .
b * 755 * , 306 ( 2016 ) [ arxiv:1601.00386 [ hep - ph ] ] .
ko , y.  omura and c.  yu , jhep * 1604 * , 098 ( 2016 ) [ arxiv:1601.00586 [ hep - ph ] ] .
ko and t.  nomura , phys .
b * 758 * , 205 ( 2016 ) [ arxiv:1601.02490 [ hep - ph ] ] .
i.  dorsner , s.  fajfer and n.  kosnik , phys .
d * 94 * , no . 1 , 015009 ( 2016 ) [ arxiv:1601.03267 [ hep - ph ] ] .
t.  nomura and h.  okada , arxiv:1601.04516 [ hep - ph ] .
x.  f.  han , l.  wang and j.  m.  yang , phys .
b * 757 * , 537 ( 2016 ) [ arxiv:1601.04954 [ hep - ph ] ] .
g.  belanger and c.  delaunay , arxiv:1603.03333 [ hep - ph ] .
m.  bauer and m.  neubert , phys .
d * 93 * , no .
11 , 115030 ( 2016 ) [ arxiv:1512.06828 [ hep - ph ] ] . c.  w.  murphy , phys
.  lett .
b * 757 * , 192 ( 2016 ) [ arxiv:1512.06976 [ hep - ph ] ] .
w.  chao , arxiv:1512.08484 [ hep - ph ] . c.  hati , phys .
d * 93 * , no . 7 , 075002 ( 2016 ) [ arxiv:1601.02457 [ hep - ph ] ] .
f.  f.  deppisch , s.  kulkarni , h.  pas and e.  schumacher , phys .
d * 94 * , no . 1 , 013003 ( 2016 ) [ arxiv:1603.07672 [ hep - ph ] ] .
u.  k.  dey , s.  mohanty and g.  tomar , arxiv:1606.07903 [ hep - ph ] .
a.  di iura , j.  herrero - garcia and d.  meloni , arxiv:1606.08785 [ hep - ph ] .
s.  baek and k.  nishiwaki , phys .
d * 93 * , no .
1 , 015002 ( 2016 ) [ arxiv:1509.07410 [ hep - ph ] ] .
s.  baek , t.  nomura and h.  okada , arxiv:1604.03738 [ hep - ph ] .
v.  khachatryan _ et al . _
[ cms and lhcb collaborations ] , nature * 522 * , 68 ( 2015 ) [ arxiv:1411.4413 [ hep - ex ] ] .
a.  belyaev , n.  d.  christensen and a.  pukhov , comput .
commun .   * 184 * , 1729 ( 2013 ) [ arxiv:1207.6082 [ hep - ph ] ] .
m.  nadolsky , h.  l.  lai , q.  h.  cao , j.  huston , j.  pumplin , d.  stump , w.  k.  tung and c .- p .
yuan , phys .
d * 78 * , 013004 ( 2008 ) [ arxiv:0802.0007 [ hep - ph ] ] .
a.  djouadi , phys .
* 457 * , 1 ( 2008 ) [ hep - ph/0503172 ] .