\section{Introduction}

The total cross section $\sigma_{tot}(s)$ is a fundamental quantity in collisions of strongly interacting particles. At present one of the main theoretical approaches for the
description of $\sigma_{tot}$ is the QCD-inspired formalism \cite{durand,luna001,luna009,giulia001}. In this approach the energy dependence of the total cross section
is obtained from QCD via an eikonal formulation.
The high-energy dependence of $\sigma_{tot}$ is driven by the rapid increase in gluon density at small-$x$. In this work we explore the
non-perturbative dynamics of QCD in order to describe, in both $pp$ and $\bar{p}p$ channels, the total cross sections $\sigma_{tot}^{pp,\bar{p}p}(s)$ and the ratios of the real to
imaginary part of the forward scattering amplitude, $\rho^{pp,\bar{p}p}(s)$. In our calculations we introduce an energy-dependent dipole form factor which represents the overlap
density for the partons at impact parameter $b$.
The behavior of these forward quantities is derived from the QCD parton model using standard parton-parton elementary processes and an updated set of gluon distribution,
namely the CTEQ6L1 set. In our analysis we have included recent measurements of $pp$ total cross sections at the LHC by the TOTEM Collaboration \cite{TOTEM}. This work is organized
as follows: in the next section we introduce our QCD-based model where the onset of the dominance of semihard partons is managed by the dynamical gluon mass. We investigate a form
factor which the spatial distribution of gluons changes with the energy and introduce integral dispersion relations to connect the real and imaginary parts of eikonals. In
Section III we present our results and in Section IV we draw our conclusions.

\section{The QCD-inspired eikonal model}

In our model the increase of $\sigma_{tot}^{pp,\bar{p}p}(s)$ is associated with elementary semihard processes in the hadrons. In the eikonal representation $\sigma_{tot}(s)$
and $\rho(s)$ are given by
\begin{eqnarray}
\sigma_{tot}(s)   =  4\pi   \int_{_{0}}^{^{\infty}}   \!\!  b\,   db\,
[1-e^{-\chi_{_{R}}(s,b)}\cos \chi_{_{I}}(s,b)]
\label{eq01}
\end{eqnarray}
and
\begin{eqnarray}
\rho(s) = \frac{-\int_{_{0}}^{^{\infty}}   \!\!  b\,
db\, e^{-\chi_{_{R}}(s,b)}\sin \chi_{_{I}}(s,b)}{\int_{_{0}}^{^{\infty}}   \!\!  b\,
db\,[1-e^{-\chi_{_{R}}(s,b)}\cos \chi_{_{I}}(s,b)]} ,
\label{eq03}
\end{eqnarray}
respectively, where $s$ is the square of the total center-of-mass energy, $b$ is the impact parameter, and
$\chi(s,b)=\textnormal{Re}\, \chi(s,b) + i\textnormal{Im}\, \chi(s,b) \equiv\chi_{_{R}}(s,b)+i\chi_{_{I}}(s,b)$ is the (complex) eikonal function. The eikonal functions for $pp$ and
$\bar{p}p$ scatterings are the sum of the soft and semihard (SH) parton interactions in the hadron-hadron collision,
\begin{eqnarray}
\chi(s,b) = \chi_{_{soft}}(s,b) + \chi_{_{SH}}(s,b).
\end{eqnarray}
It follows from the QCD parton model that the real part of the eikonal, $\textnormal{Re}\,\chi_{_{SH}}(s,b)$, can be factored as
\begin{eqnarray}
\textnormal{Re}\,\chi_{_{SH}}(s,b) = \frac{1}{2}\, W_{\!\!_{SH}}(b)\,\sigma_{_{QCD}}(s),
\end{eqnarray}
where $W_{\!\!_{SH}}(b)$ is an overlap density for the partons at impact parameter space $b$,
\begin{eqnarray}
W_{\!\!_{SH}}(b) &=& \int d^{2}b'\, \rho_{A}(|{\bf b}-{\bf b}'|)\, \rho_{B}(b'),
\label{ref009}
\end{eqnarray}
and $\sigma_{_{QCD}}(s)$ is the usual QCD cross section
\begin{eqnarray}
\sigma_{_{QCD}}(s) &=& \sum_{ij} \frac{1}{1+\delta_{ij}} \int_{0}^{1}\!\!dx_{1}
\int_{0}^{1}\!\!dx_{2} \int_{Q^{2}_{min}}^{\infty}\!\!d|\hat{t}|
\frac{d\hat{\sigma}_{ij}}{d|\hat{t}|}(\hat{s},\hat{t}) \nonumber \\
 &\times & f_{i/A}(x_{1},|\hat{t}|)f_{j/B}(x_{2},|\hat{t}|)\, \Theta \! \left( \frac{\hat{s}}{2} - |\hat{t}| \right),
\label{eq08}
\end{eqnarray}
with $|\hat{t}|\equiv Q^{2}$, $\hat{s}=x_{1}x_{2}s$ and $i,j=q,\bar{q},g$.

The eikonal functions for $pp$ and $\bar{p}p$ scatterings are written as $\chi_{pp}^{\bar{p}p}(s,b) = \chi^{+} (s,b) \pm \chi^{-} (s,b)$, with
$\chi^{+}(s,b) = \chi^{+}_{_{soft}}(s,b) + \chi^{+}_{_{SH}}(s,b)$ and $\chi^{-}(s,b) = \chi^{-}_{_{soft}}(s,b) + \chi^{-}_{_{SH}}(s,b)$. Since in the parton model
$\chi^{-}_{_{SH}}(s,b)$ decreases rapidly with increasing $s$, the difference between $pp$ and $\bar{p}p$ cross sections is due only to the different weighting of the quark-antiquark
(valence) annihilation cross sections in the two channels. Thus the crossing-odd eikonal $\chi^{-}(s,b)$ receives no contribution from semihard processes at high energies, and it
is sufficient to take $\chi_{_{SH}}(s,b)=\chi^{+}_{_{SH}}(s,b)$ and, consequently, $\chi^{-}(s,b) = \chi^{-}_{_{soft}}(s,b)$. The connection between the real and imaginary parts of
the crossing-even eikonal $\chi^{+}(s,b)$ is obtained by means of dispersion relation
\begin{eqnarray}
\textnormal{Im}\,\chi^{+}(s,b) = -\frac{2s}{\pi}\, {\cal P}\!\! \int_{0}^{\infty}ds'\,
\frac{\textnormal{Re}\,\chi^{+}(s',b)}{s^{\prime 2}-s^{2}} ,
\label{idr001}
\end{eqnarray}
valid at $s \gg 1$ GeV$^{2}$. The soft term is given by
\begin{eqnarray}
\chi^{+}_{_{soft}}(s,b) = \frac{1}{2}\, W^{+}_{\!\!_{soft}}(b;\mu^{+}_{_{soft}})\, \left[ A' +\frac{B'}{(s/s_{0})^{\gamma}}\, e^{i\pi\gamma/2}
-i C'\left[ \ln\left(\frac{s}{s_{0}}\right) -i\frac{\pi}{2} \right] \right] ,
\label{soft01}
\end{eqnarray}
where $\sqrt{s_{0}}\equiv 5$ GeV and $A'$, $B'$, $C'$, $\gamma$ and $\mu^{+}_{_{soft}}$ are fitting parameters. The odd eikonal $\chi^{-}(s,b)$ is given by
\begin{eqnarray}
\chi^{-}(s,b) &=& \frac{1}{2}\, W^{-}_{\!\!_{soft}}(b;\mu^{-}_{_{soft}})\,D'\, \frac{e^{-i\pi/4}}{\sqrt{s}},
\label{softminus}
\end{eqnarray}
where $\mu^{-}_{_{soft}}\equiv 0.5$ GeV and $D'$ is also a fitting parameter.

In the expression (\ref{eq08}) the term $d\hat{\sigma}_{ij}/d|\hat{t}|$ is the differential cross section for $ij$ scattering ($i,j=q,\bar{q},g$). We select parton-parton processes
containing at least one gluon in the initial state, i.e., we consider the processes
\begin{eqnarray}
\frac{d\hat{\sigma}}{d\hat{t}}(gg\to gg)=\frac{9\pi\bar{\alpha}^{2}_{s}}{2\hat{s}^{2}}\left(3 -\frac{\hat{t}\hat{u}}{\hat{s}^{2}}-
\frac{\hat{s}\hat{u}}{\hat{t}^{2}}-\frac{\hat{t}\hat{s}}{\hat{u}^{2}} \right) ,
\label{pp001}
\end{eqnarray}
\begin{eqnarray}
\frac{d\hat{\sigma}}{d\hat{t}}(qg\to qg)=\frac{\pi\bar{\alpha}^{2}_{s}}{\hat{s}^{2}}\, (\hat{s}^{2}+\hat{u}^{2}) \left(
\frac{1}{\hat{t}^{2}}-\frac{4}{9\hat{s}\hat{u}} \right) ,
\label{pp002}
\end{eqnarray}
\begin{eqnarray}
\frac{d\hat{\sigma}}{d\hat{t}}(gg\to \bar{q}q)=\frac{3\pi\bar{\alpha}^{2}_{s}}{8\hat{s}^{2}}\, (\hat{t}^{2}+\hat{u}^{2}) \left(
\frac{4}{9\hat{t}\hat{u}}-\frac{1}{\hat{s}^{2}} \right) ,
\label{pp003}
\end{eqnarray}
where
\begin{eqnarray}
\bar{\alpha}_{s} (Q^{2})= \frac{4\pi}{\beta_0 \ln\left[
(Q^{2} + 4M_g^2(Q^{2}) )/\Lambda^2 \right]},
\label{eq26}
\end{eqnarray}
\begin{eqnarray}
M^2_g(Q^{2}) =m_g^2 \left[\frac{ \ln
\left(\frac{Q^{2}+4{m_g}^2}{\Lambda ^2}\right) } {
\ln\left(\frac{4{m_g}^2}{\Lambda ^2}\right) }\right]^{- 12/11} ;
\label{mdyna}
\end{eqnarray}
here $\bar{\alpha}_{s} (Q^{2})$ is a non-perturbative QCD effective charge, obtained by Cornwall through the use of the pinch technique in order to derive a gauge invariant
Schwinger-Dyson equation for the gluon propagator \cite{cornwall}.

In the case of semihard gluons we consider the possibility of a ``broadening'' of the spatial distribution. Our assumption suggests an increase of the average gluon radius
when $\sqrt{s}$ increases, and can be properly implemented using the energy-dependent dipole form factor
\begin{eqnarray}
G^{(d)}_{\!\!_{SH}}(s,k_{\perp};\nu_{\!\!_{SH}})=\left( \frac{\nu_{_{SH}}^{2}}{k_{\perp}^{2}+\nu_{_{SH}}^{2}} \right)^{2},
\end{eqnarray}
where $\nu_{_{SH}}= \nu_{1}-\nu_{2}\ln ( \frac{s}{s_{0}} )$, with $\sqrt{s_{0}}\equiv 5$ GeV. Here $\nu_{1}$ and $\nu_{2}$ are constants to be fitted. From equation (\ref{ref009}) we have
\begin{eqnarray}
W_{\!\!_{SH}}(s,b;\nu_{_{SH}}) &=& \frac{1}{2\pi}\int_{0}^{\infty}dk_{\perp}\, k_{\perp}\, J_{0}(k_{\perp}b)\,[G^{(d)}_{\!\!_{SH}}(s,k_{\perp};\nu_{\!\!_{SH}}]^{2} \nonumber \\
 &=& \frac{\nu^{2}_{_{SH}}}{96\pi} (\nu_{_{SH}} b)^{3} K_{3}(\nu_{_{SH}} b).
\end{eqnarray}
Since semihard interactions dominate at high energies, we consider an energy-dependence only in the case of $W_{\!\!_{SH}}(s,b;\nu_{_{SH}})$. In this way the soft overlap densities
$W^{+}_{\!\!_{soft}}(b;\mu^{+}_{_{soft}})$ and $W^{-}_{\!\!_{soft}}(b;\mu^{-}_{_{soft}})$ are static, i.e., $\mu^{+}_{_{soft}}$ and $\mu^{-}_{_{soft}}$ are not energy dependent.

\section{Results}

First, in order to determine the model parameters, we set the value of the gluon mass scale to $m_{g} = 400$ MeV. This choice for the mass scale is not only
consistent to our LO procedure, but is also the one usually obtained in other calculations of strongly interacting
processes \cite{luna001,luna009,luna010,luna2014}. Next, we carry out a global fit to high-energy forward $pp$ and
$\bar{p}p$ scattering data above $\sqrt{s} = 10$ GeV, namely the total cross section $\sigma_{tot}^{pp,\bar{p}p}$ and
the ratio of the real to imaginary part of the forward scattering amplitude $\rho^{pp,\bar{p}p}$. We use data sets
compiled and analyzed by the Particle Data Group \cite{PDG} as well as the recent data at LHC from the TOTEM
Collaboration, with the statistic and systematic errors added in quadrature. We include in the dataset the first estimate for the $\rho$ parameter made by the TOTEM
Collaboration in their $\rho$-independent measurement at $\sqrt{s}=7$ TeV, namely $\rho^{pp}=0.145\pm0.091$ \cite{TOTEM}.

The values of the fitted parameters are given in the legend of Figure 1. The result $\chi^{2}/DOF = 1.062$ for the fit was obtained
for 154 degrees of freedom. The results of the fit to $\sigma_{tot}$ and $\rho$ for both $pp$ and $\bar{p}p$ channels are
displayed in Figure 1, together with the experimental data. The curves depicted in Figure 1 were all calculated using the cutoff $Q_{min}=1.3$ GeV, the value of the CTEQ6 fixed
initial scale $Q_{0}$.

\begin{figure}[!h]
\centering
\hspace*{-0.2cm}\includegraphics[scale=0.60]{bblfigure006b.eps}
\label{fig1}
\caption{The total cross section (upper) and the ratio of the real to imaginary part of the forward scattering amplitude (lower), for $pp$ ($\bullet$) and $\bar{p}p$ ($\circ$).
The values of the fitted parameters are: $\nu_{1} = 2.770\pm0.865$ GeV, $\nu_{2} = (7.860\pm5.444)\times 10^{-2}$ GeV, $A' = 108.9\pm8.6$ GeV$^{-1}$, $B' = 30.19\pm15.78$ GeV$^{-1}$,
$C' = 1.260\pm0.437$ GeV$^{-1}$, $\gamma = 0.719\pm 0.200$, $\mu^{+}_{soft} = 0.457\pm0.209$ GeV and $D' = 21.73\pm2.26$ GeV$^{-1}$.}
\end{figure}

\section{Conclusions}

In this work we have investigated infrared effects on semihard parton-parton interactions. We have studied $pp$ and $\bar{p}p$ scattering in the LHC energy region with the
assumption that the increase of their total cross sections arises exclusively from semihard interactions.
In the calculation of $\sigma_{tot}^{pp,\bar{p}p}$ and $\rho^{pp,\bar{p}p}$ we have considered the phenomenological implication of an energy-dependent form factor for semihard
partons. We introduce integral dispersion relations to connect the real and imaginary parts of eikonals with energy-dependent form factors.
In our analysis we have included the recent data at LHC from the TOTEM Collaboration.
Our results show that very good descriptions of $\sigma_{tot}^{pp,\bar{p}p}$ and
$\rho^{pp,\bar{p}p}$ data are obtained by constraining the value of the cutoff at $Q_{min} = 1.3$ GeV.
The $\chi^{2}/DOF$ for the global fit was 1.062 for 154 degrees of freedom. This good
statistical result shows that our eikonal model, where non-perturbative effects are naturally included via a QCD
effective charge, is well suited for detailed predictions of the forward quantities to be measured at higher energies.
In the semihard sector we have considered a form factor in which the average gluon radius increases
with $\sqrt{s}$. With this assumption we have obtained another form in which the eikonal can be factored into the QCD
parton model, namely
$\textnormal{Re}\,\chi_{_{SH}}(s,b) = \frac{1}{2}\, W_{\!\!_{SH}}(s,b)\,\textnormal{Re}\,\sigma_{_{QCD}}(s)$. The imaginary
part of this {\it semi-factorizable} eikonal was obtained by means of appropriate integral dispersion relations
which take into account eikonals with energy-dependent form factors.

At the moment we are analyzing the effects of different updated sets of parton distributions on the
forward quantities, namely CTEQ6L and MSTW, investigating the uncertainty on these forward observables coming from the uncertainties associated with the dynamical
mass scale and the parton distribution functions, and exploring a new ansatz for the energy-dependent form factor.

\begin {thebibliography}{99}

\bibitem{durand} L. Durand and H. Pi, Phys. Rev. Lett. {\bf 58}, 303 (1987);
Phys. Rev. D {\bf 38}, 78 (1988);
{\bf 40}, 1436 (1989).

\bibitem{luna001} E.G.S. Luna, A.F. Martini, M.J. Menon, A. Mihara, and A.A. Natale, Phys. Rev. D {\bf 72}, 034019
(2005);
E.G.S. Luna, Phys. Lett. B {\bf 641}, 171 (2006);
E.G.S. Luna and A.A. Natale, Phys. Rev. D {\bf 73}, 074019 (2006);
E.G.S. Luna, Braz. J. Phys. {\bf 37}, 84 (2007);
D.A. Fagundes, E.G.S. Luna, M.J. Menon, and A.A. Natale, Nucl. Phys. A {\bf 886}, 48 (2012);
E.G.S. Luna and P.C. Beggio, Nucl. Phys. A {\bf 929}, 230 (2014).

\bibitem{luna009} E.G.S. Luna, A.L. dos Santos, and A.A. Natale, Phys. Lett. B {\bf 698}, 52 (2011).

\bibitem{giulia001} A. Corsetti, A. Grau, G. Pancheri, and Y.N. Srivastava, Phys. Lett. B {\bf 382}, 282 (1996);
A. Grau, G. Pancheri, and Y.N. Srivastava, Phys. Rev. D {\bf 60}, 114020 (1999);
R.M. Godbole, A. Grau, G. Pancheri, and Y.N. Srivastava, Phys. Rev. D {\bf 72}, 076001 (2005);
A. Achilli {\it et al.}, Phys. Lett. B {\bf 659}, 137(2008);
A. Grau, R.M. Godbole, G. Pancheri, and Y.N. Srivastava, Phys. Lett. B {\bf 682}, 55 (2009);
G. Pancheri, D.A. Fagundes, A. Grau, S. Pacetti, and Y.N. Srivastava, arXiv:1301.2925 [hep-ph].

\bibitem{TOTEM} G. Antchev {\it et al.}, EPL {\bf 96}, 21002 (2011);
G. Antchev {\it et al.}, EPL {\bf 101}, 21002 (2013);
G. Antchev {\it et al.}, EPL {\bf 101}, 21004 (2013);
G. Antchev {\it et al.}, Phys. Rev. Lett. {\bf 111}, 012001 (2013).

\bibitem{cornwall} J. M. Cornwall, Phys. Rev. D {\bf 22}, 1452 (1980); D {\bf 26}, 1453 (1982).

\bibitem{luna010} A. Doff, E.G.S. Luna, and A.A. Natale, Phys. Rev. D {\bf88}, 055008 (2013).

\bibitem{luna2014} E.G.S. Luna and A.A. Natale, J. Phys. G {\bf 42}, 015003 (2015).

\bibitem{PDG} K.A. Olive {\it et al.}, Chin. Phys. C {\bf 38}, 090001 (2014).

\end{thebibliography}
\smallskip

\end{document}


\section{INTRODUCTION}

Since the failure of a second reaction wheel, the \textit{Kepler} telescope provides the community with K2 observations of fields along the ecliptic plane \citep{Howell14}, and is therefore observing some previously known WASP exoplanets. One result has been the discovery of two additional transiting planets in the WASP-47 system \citep{Becker15}.

The transiting hot Jupiter WASP-85Ab was discovered by \citet{Brown14} and was observed by K2 in Campaign 1. It has a 2.66-d orbit around its G5, \textit{V}\thinspace =\thinspace 11.2 host star. WASP-85A forms a close visual binary, with an angular separation of 1.5\thinspace arcsec, with a cooler and dimmer K0, \textit{V}\thinspace =\thinspace 11.9 WASP-85B. The WASP photometry showed that one of the two binary components exhibits a 14.6\thinspace $\pm$\thinspace 1.5\thinspace d rotational modulation, which indicates stellar activity \citep{Brown14}.

When a starspot is occulted by a transiting planet the lightcurve exhibits a ``bump'' of temporary brightening \citep{Silva03}. The detection of starspots can provide a measurement or a constraint on the stellar obliquity, i.e. the angle between the stellar rotation and the planet's orbital axis. For example, the recurrence of the same starspot in consecutive transits has revealed that Kepler-17 has a low obliquity \citep{Desert11}. The ability to detect the recurring starspot occultation in a particular system depends on the ratio between the stellar rotational and planet orbital periods, and the lifetime of starspots. Alternatively, if starspots are occulted at certain preferential phases of the transit, this may constrain the obliquity if the transit chord crosses the active stellar latitudes at the observed active phases, as in the case of HAT-P-11 \citep{Sanchis11}.

The majority of hot Jupiters around stars cooler than around 6250\thinspace K are found in aligned orbits, suggestive of a disc migration mechanism. Hotter stars, on the other hand, host hot Jupiters that are in many cases found in misaligned orbits \citep{Albrecht2012}, which cannot be accounted for solely by planet disc migration mechanism. Instead, it is believed that planet--planet scattering, Kozai mechanism and tidal dissipation play an additional role in shaping misaligned orbits of hot Jupiters (e.g. \citealt{Triaud10}).

In this paper, we present a detection of starspots in the K2 lightcurve of WASP-85, which we use to constrain the stellar obliquity. We also refine the system parameters and the stellar rotational period, and search for transit-timing variations, as could be caused by an additional body in the system.

\section{THE K2 OBSERVATIONS}

The K2 Campaign 1 observations cover an 82-d time-span between 2014 May 30 and 2014 August 20 and provide a total of 120\thinspace 660 short-cadence images. We retrieved the 1-min short-cadence target pixel file for WASP-85 via the Minkulski Archive for Space Telescopes (MAST).

Given the degraded pointing accuracy of K2, the photometric precision is reduced because the position of an observed object drifts within the photometric aperture, and because of the non-uniform pixel sensitivity of the on-board detectors. \citet{Vanderburg14} developed a self-flat-fielding (SFF) procedure with which they correct the long-cadence photometry to within a factor of two of the original \textit{Kepler} precision by correlating the measured flux and arclength of the drift. Our K2 data-reduction procedure for short-cadence observations is similar to the procedure presented by \citet{Vanderburg14}.

\begin{figure*}
\includegraphics[width=\textwidth]{f1.png}
\caption{SFF calibration procedure. All the panels correspond to a time interval of 5\thinspace d. 1a: Position of the center of the object given as CCD row vs. CCD column. Spatial outliers are marked with red and flagged as thruster events. 1b: Measured normalized flux vs. arclength of the spacecraft drift. The correlation is fitted with a Gaussian convolution (marked with red). 1c: Lightcurve before the SFF correction with the quality-flagged data points removed. 1d: Lightcurve after the SFF correction with the quality-flagged data points removed, including the thruster events which were identified during the SFF procedure. The positive flux excursions in the corrected lightcurve are the remaining cosmic rays that have avoided being identified and quality-flagged as cosmic rays in the downloaded K2 target pixel file.}
\end{figure*}

We first applied a standard data reduction procedure using PyRAF tools for \textit{Kepler} (PyKE) version 2.6.2 \citep{Still12}. Besides removing all the frames that were originally quality-flagged by K2, we manually adjusted quality flags and removed an additional 100 frames between 2014 June 1 18:04\thinspace UT and 19:41\thinspace UT, since they were obtained while the spacecraft was still in coarse pointing after the attitude tweak. Next, we used a large 100-pixels aperture mask to accommodate the diluted photometric signal from both binary components, WASP-85A and WASP-85B, and to reduce the photometric losses due to the relative drift between the object and the extraction aperture. Photometry of the individual binary components could not be performed because the angular separation (1.5\thinspace arcsec) is smaller than the \textit{Kepler} point spread function (full width half maximum of 6~arcsec).

After the data extraction and removal of low-frequency modulations, we performed our SFF procedure. The procedure is based on the PyKE tool \textsc{kepsff} with two main modifications to adjust and improve the efficiency for the short-cadence data. First, the 1-min short-cadence frames obtained during the thruster firing events are seen to have object positions that are outliers from the mean position drift on the detector. Data points with positional deviations of more than 6$\sigma$ were flagged as thruster firing events and deleted (see Figure~1a). This level was established by trial and error, as rejecting the thruster firing events while removing a minimal number of good data points. A total of 490 data points were quality-flagged as thruster firing events, which represents 0.4 per cent of the entire lightcurve.

\begin{figure}
\includegraphics[width=8.5cm]{f2.png}
\caption{K2 Campaign 1 short-cadence lightcurve of WASP-85 before (shown in black) and after the SFF correction (red). Note the presence of transits and starspot rotational modulation. The corrected lightcurve is shown with an offset of \mbox{$-$20\thinspace 000}\thinspace e\textsuperscript{-}\thinspace s\textsuperscript{-1} for clarity.}
\end{figure}

Second, instead of using the polynomial to fit the flux versus arclength of the positional drift of the spacecraft, we used convolution with a Gaussian kernel, which provided an improved fit of the correlation compared to the polynomial. For each data point the Gaussian convolution calculates the peak of the Gaussian flux distribution of a specified number of data points nearest in arclength to a particular data point. At this stage the deviating data points are typically of astrophysical origin and should be rejected to avoid any bias in the SFF calibration. Trial and error indicated that the flux-versus-drift arclength correlation was best fitted when using a Gaussian kernel with a width of 50 data points and a rejection threshold of 4$\sigma$ (see Figure~1b). After the Gaussian peaks have been calculated for all the data points the flux of every data point was replaced with a corresponding Gaussian peak. Figures~1c and 1d show an example lightcurve of a 5-d calibration time interval, before and after SFF, respectively. The time interval of 5\thinspace d was chosen as the best compromise between including a larger statistical sample of spacecraft drifts and a long-term variable drift behaviour.

With these two modifications we achieved nearly twice the precision improvement compared to the original \textsc{kepsff} PyKE tool for the short-cadence dataset of WASP-85. The median 1-min photometric precision for WASP-85 (\textit{K}$_{p}$\thinspace =\thinspace 10.2) was 778 parts per million (ppm) before SFF and 260\thinspace ppm after SFF, an improvement by a factor 3. For a comparison, the 1-min photometric precision of the original \textit{Kepler} mission for a \textit{K}$_{p}$\thinspace =\thinspace 10.5 star was 202\thinspace ppm \citep{Still13}.

After removal of all the quality-flagged data points, i.e. images affected by the cosmic rays, coarse pointing, spacecraft attitude tweak, thruster firing and reaction wheel de-saturation events, and a downlink data gap near the middle of Campaign 1, the total number of useful data points was reduced from 120\thinspace 660 to 108\thinspace 872. Figure~2 shows the entire short-cadence lightcurve of WASP-85 before and after applying the SFF correction. A total of 30 transits are visible. Also visible is the rotational modulation caused by the presence of starspots.

\section{ANALYSIS AND RESULTS}
\subsection{System parameters}

We obtained the system parameters using the Markov Chain Monte Carlo (MCMC) code presented in \citet{CollierCameron07} and further described in \citet{Pollacco08} and \citet{Anderson15}. In addition to the high-quality K2 photometry presented herein we used the HARPS radial-velocity measurements of WASP-85A from \citet{Brown14}. These measurements are not thought to have been significantly contaminated by WASP-85B. We corrected the K2 lightcurve for dilution by WASP-85B using a flux ratio of 0.50 (see \citet{Brown14}). Because the presence of starspot occultations may lead to inaccurate determination of system parameters \citep{Oshagh13} we removed the starspot occultation events from the measured lightcurve prior to the MCMC analysis. Limb darkening was accounted for using a four-parameter law, with coefficients calculated for \textit{Kepler} bandpass and interpolated from the tabulations of \citet{Sing10}.

To refine the ephemeris (but not the other parameters), we produced a separate MCMC analysis that also included the ground-based photometry presented in \citet{Brown14}, with four-parameter limb-darkening coefficients from \citet{Claret00,Claret04}, as appropriate for different bandpasses. This extended the photometric baseline from 80\thinspace d to 6.5\thinspace yr and reduced the uncertainty on the period by a factor of 4.

We present the system parameters in Table~1 and superimpose the corresponding transit model on the K2 lightcurve in Figure~3.

There is a discrepancy between the stellar effective temperature derived from the HARPS spectra ($T_{\rm eff} = 5685 \pm 65$\thinspace K; \citealt{Brown14}) and the value obtained from fitting the transits in our MCMC analysis ($T_{\rm eff} = 6112 \pm 27$\thinspace K). The same discrepancy, albeit smaller, is present in the discovery paper, which used only the ground-based transits ($T_{\rm eff} = 5910 \pm 54$\thinspace K from MCMC; \citealt{Brown14}). When fitting the photometry (either ground-based or K2) with a transit model generated using the spectroscopic $T_{\rm eff}$, large systematics are evident during ingress and egress (see upper residual panel in Figure~3), suggestive of a poor modelling of the limb darkening. The $T_{\rm eff}$ discrepancy could be explained in part if light from WASP-85B had contaminated the spectra of WASP-85A, leading to a lower apparent $T_{\rm eff}$ for the primary, though any contamination is expected to be at a low level.

\begin{figure}
\includegraphics[width=8.5cm]{f3.png}
\caption{Best-fit MCMC model and its residuals from the phase-folded K2 lightcurve near the transit. The lightcurve has been corrected for dilution by WASP-85B, using a flux ratio of 0.50. The lower residual panel corresponds to the effective temperature of 6112\thinspace K as fitted by our MCMC analysis. The upper residual panel serves as a comparison and corresponds to the spectroscopic effective temperature of 5685\thinspace K which results in large residuals near ingress and egress.}
\end{figure}

\begin{table*}
\centering
\begin{minipage}{11cm}
\caption{System parameters from MCMC for WASP-85Ab and its host star}
\begin{tabular}{@{}lccccc@{}}
\hline
Parameter&Symbol&Value&Unit\\
\hline
Transit epoch$^{a}$&\textit{t}$_{\rm 0}$&2456847.472856 $\pm$ 0.000014&BJD\\
Orbital period$^{a}$&\textit{P}&2.65567770 $\pm$ 0.00000044&d\\
Area ratio&$(R_{\rm p}/R_{\star})^{2}$&0.01870 $\pm$ 0.00002&...\\
Transit width&\textit{t}$_{\rm T}$&0.10816$^{+0.00006}_{-0.00002}$&d\\
Ingress and egress duration&\textit{t}$_{\rm 12}$, \textit{t}$_{\rm 34}$&0.0130377$^{+0.000008}_{-0.000016}$&d\\
Impact parameter&\textit{b}&0.047$^{+0.005}_{-0.018}$&...\\
Orbital inclination&\textit{i}&89.69$^{+0.11}_{-0.03}$&$^{\circ}$\\
Orbital eccentricity&\textit{e}&0 (adopted)&...\\
Orbital separation&\textit{a}&0.039 $\pm$ 0.001&AU\\
Stellar effective temperature$^{b}$&\textit{T}$_{\rm eff}$&6112 $\pm$ 27&K\\
Stellar mass&\textit{M}$_{\star}$&1.09 $\pm$ 0.08&M$_\odot$\\
Stellar radius&\textit{R}$_{\star}$&0.935 $\pm$ 0.023&R$_\odot$\\
Stellar density&$\rho_{\star}$&1.330 $\pm$ 0.007&$\rho_\odot$\\
Planet equilibrium temperature&\textit{T}$_{\rm p}$&1452 $\pm$ 6&K\\
Planet mass&\textit{M}$_{\rm p}$&1.265 $\pm$ 0.065&M$_{\rm Jup}$\\
Planet radius&\textit{R}$_{\rm p}$&1.24 $\pm$ 0.03&R$_{\rm Jup}$\\
Planet density&$\rho_{\rm p}$&0.66 $\pm$ 0.02&$\rho_{\rm Jup}$\\
\hline
\end{tabular}
\begin{itemize}[leftmargin=0.35cm]
\setlength\itemsep{0cm}
\item[$^{a}$]Epoch and period derived by fitting the photometric datasets from the K2 and all the available ground-based observations.
\item[$^{b}$]This value, obtained from a fit to the transit lightcurve, is discrepant with the spectroscopic value ($T_{\rm eff} = 5685 \pm 65$\thinspace K).
\end{itemize}
\end{minipage}
\end{table*}

\subsection{TTV and TDV}

Transit-timing variations (TTVs) and transit-duration variations (TDVs) can arise from dynamical interactions with additional planets in the system (e.g. \citealt{Becker15} and citations therein). Typical reported TTV amplitudes for perturbed transiting exoplanets are of the order of a few tens of minutes with periods of a few hundred days \citep{Mazeh13}. TDVs on the other hand are expected to be observed with smaller amplitudes and in phase with TTV. The first unequivocal exoplanet TDV detection was measured for system Kepler-88 and indicated a TDV amplitude of 5\thinspace min and a period of 630\thinspace d \citep{Nesvorny13}.

We measured the transit-timings and transit-durations for WASP-85Ab for each of the transits using the MCMC analysis similarly as described in Section~3.1, and show the deviations in Figure~4. The detected starspot occultations have again been removed from the lightcurve prior to the analysis as in Section~3.1. The error bars in Figure~4 do not include any contribution due to the remaining unremoved stellar activity features, which will be contributing to the observed scatter of the TTV and TDV data points. One data point has been excluded from TTV and TDV analysis due to the thruster firing event residual that took place at the beginning of ingress of the fifth transit near BJD~2456821.

\begin{figure}
\includegraphics[width=8.5cm]{f4.png}
\caption{TTV (upper panel) and TDV (lower panel) for WASP-85Ab. The measurements indicate a sinusoidal TTV variability with a small semi-amplitude of 7.8\thinspace s and a period 14.7\thinspace d, likely a signature of stellar activity rather than actual TTVs. The TDV measurements are consistent with the assumption of constant transit-duration. Note that the data point with the largest error bars corresponds to the transit during which a reaction wheel de-saturation took place near the ingress (see transit 17 in Figure~5).}
\end{figure}

The hypothesis of the TDV values being zero gives a $\chi^2$ of 31.8 for 28 degrees of freedom, which means we cannot reject it. The $\chi^2$ for TTV is 90.6 for 28 degrees of freedom which rejects the white noise distribution hypothesis and suggests a sinusoidal TTV signal with a semi-amplitude of 7.8\thinspace s and a period of 14.7\thinspace d. This TTV signal, if real, would to the best of our knowledge have the smallest detected TTV amplitude to date. On the other hand, it is known that stellar activity can induce TTV signals with semi-amplitudes of up to 200\thinspace s \citep{Oshagh13} and that the imperfect removal of some of the starspot occultation events, and the effect of the stellar rotational modulation, could have resulted in the TTV variability. The TTV period of 14.7\thinspace d is also compatible with the measured rotational period of $15.1 \pm 0.6$\thinspace d (see Section~3.4) and it is therefore more likely that the small measured TTV signal results from stellar activity.

The very small TTV amplitude and the absence of statistically significant TDV indicates that a second, non-transiting, massive planet is unlikely, unless at a much longer period.

\subsection{Starspot detection}

The short 1-min cadence of K2 observations and the restored photometric precision using the SFF enabled us to detect individual starspot occultations. We show in Figure~5 the transit lightcurves after subtracting the transit model, thus showing only the residuals. We first looked for starspots by calculating the $\chi^2$ value compared to a straight line. This was done by setting a box width equal to a phase width of 0.09, and sliding that box in time. The box width was determined by examining a few well defined starspot occultation events in the system. To estimate the threshold of significance we used the same box on out-of-transit data. Using this approach we were able to detect the 12 most prominent starspot occultation events above a $\chi^2_{\rm red}$ threshold of 4 (marked with asterisks in Figure~5). Below this threshold the residual lightcurve proved to be too unstable for this occultation detection technique.

To identify additional occultation events we asked four colleagues to mark definite and possible occultation events in the residual lightcurves near transits. The order of lightcurves was shuffled to avoid any bias towards spots that repeat in consecutive transits. All the events which were identified by at least two people as a possible starspot occultation event are listed in Table~2 and marked with red ellipses in Figure~5. We identified a total of 23 starspot occultations. The measured $\chi^2_{\rm red}$ of the individual occultation events are given in Table~2 and starspots are marked with different color intensities in Figure~5 as a rough indicator of their detectability.

\begin{figure*}
\includegraphics[width=\textwidth]{f5.png}
\caption{Starspot occultations in the model-subtracted lightcurve of WASP-85Ab. Vertical dashed lines specify the extent of the transit. Red ellipses mark the 23 identified starspot occultation events. The color intensity denotes the detectability of events and asterisks correspond to events detected by our $\chi^2$ occultation detection technique. Potential starspot occultation pairs are marked with red lines. Blue arrows mark the positions of thruster firing events. Note that transit number 16 lies in the K2 data gap.}
\end{figure*}

Blue arrows are added in Figure~5 to mark the positions of thruster firing events which occasionally introduce artefacts in the lightcurve. Thruster firing artefacts are especially present in the first quarter of the observing campaign due to the somewhat inconsistent spacecraft drifts which might lead to a confusion with occultation events. Nevertheless, a typical potential thruster firing artefact only causes a small step-like deviation in the lightcurve and in the out-of-transit WASP-85 lightcurve never resembles a starspot-like bump. Therefore, it seems unlikely that any of the in-transit starspot-like event would be a thruster firing artefact, despite occasional overlaps.

\begin{table}
\centering
\caption{Starspot occultation events}
\begin{tabular}{@{}cccccc@{}}
\hline
Transit&$\chi^2_{\rm red}$&Phase&Stellar&Pair&Differential\\
number&&&longitude$^a$ [$^{\circ}$]&&longitude [$^{\circ}$]\\
\hline
5&1.57&0.0075&24.9&&\\
6&7.26&$-$0.0077&$-$25.6&6+7&60.7 $\pm$ 3.4\\
7&2.46&$-$0.0005&$-$1.7&&\\
7&4.91&0.0103&35.1&&\\
10&6.65&$-$0.0059&$-$19.4&10+11&63.7 $\pm$ 4.0\\
11&8.73&0.0125&44.3&&\\
13&5.25&$-$0.0037&$-$11.9&13+14&62.8 $\pm$ 4.2\\
14&1.65&0.0139&50.9&&\\
19&5.79&0.0051&16.7&&\\
21&6.21&0.0004&1.2&&\\
23&3.16&$-$0.0028&$-$9.0&&\\
24&8.13&$-$0.0100&$-$33.9&24+25&63.8 $\pm$ 3.8\\
24&3.06&0.0062&20.3&&\\
25&5.86&$-$0.0023&$-$7.5&25+26&63.3 $\pm$ 4.5\\
25&4.21&0.0089&29.9&&\\
26&1.82&0.0148&55.8&&\\
27&2.90&$-$0.0136&$-$49.4&27+28&68.2 $\pm$ 4.2\\
28&3.10&0.0058&18.8&&\\
29&5.63&$-$0.0028&$-$9.0&29+30&62.1 $\pm$ 4.3\\
29&2.31&0.0031&10.0&&\\
30&4.95&$-$0.0064&$-$20.9&30+31&61.3 $\pm$ 3.8\\
30&2.29&0.0143&53.1&&\\
31&1.81&0.0116&40.4&&\\
\hline
\end{tabular}
\begin{itemize}[leftmargin=0.35cm]
\setlength\itemsep{0cm}
\item[$^{a}$]Stellar longitude of zero corresponds to the meridian that runs through the center of the stellar disc.
\end{itemize}
\end{table}

\subsection{Are the starspots recurring?}

The presence of 23 starspots occultation events among 30 transits indicates that WASP-85A is magnetically active, that starspots occur uniformly along the transit chord, and that some starspot occultations are potentially occurring in pairs.

If the stellar rotational period is large compared to the orbital period of the planet and if the rotational and orbital axes are aligned, it is possible to see an occultation of the same starspot in consecutive transits. Figure~5 reveals eight such potential starspot occultation pairs at consistent phase shifts. For example, an occultation event at phase 0.013 during transit 11 was likely caused by the same starspot that was responsible for the occultation event at phase $-$0.006 in transit 10. All eight potential occultation pairs are marked with red lines in Figure~5 and their phase shifts and differential stellar longitudes are given in Table~2.

A detection of eight potential occultation pairs (16 spots) among a total of 23 starspot occultation events coupled with the small scatter among measured differential stellar longitudes indicates that occultation pairs were probably caused by occultations of the recurring starspots. It is therefore reasonable to suggest that planet's orbital and stellar rotational axis are aligned. Using the system parameter values from Section~3.1 we estimate that the obliquity angle has to be $<$14$^{\circ}$. From the mean differential stellar longitude of the starspot pairs of $63.2 \pm 2.3^{\circ}$ (see Table~2) and known orbital period we derive a stellar rotational period of $15.1 \pm 0.6$\thinspace d. This is in agreement with the rotational modulation period (see Section~3.5) and with the spectroscopically determined rotational period of 14.2\thinspace $\pm$\thinspace 4\thinspace d for WASP-85A \citep{Brown14}.

\subsection{Rotational modulation}

The K2 lightcurve shows a clear rotational modulation (see Figure~2). To estimate the period we used an autocorrelation function (ACF). For this task, we first removed the transits from the lightcurve and then applied 10$\sigma$ clipping. The rotational period is manifested as the highest ACF peak near 13.6\thinspace d and its multiples (see Figure~6). Peaks corresponding to half of this period are caused by starspot regions  on opposite sides of the star. The 80-d time-span of the K2 observations covers five cycles of the 13.6-d periodicity. Following the procedure demonstrated by \citet{McQuillan13} we determined the rotational period as a median of the intervals $\Delta\tau$ between the ACF peaks associated with the rotational period. We conservatively estimated the error as the FWHM of the highest ACF peak.

The resulting rotational period of $13.6 \pm 1.6$~d is in agreement with the spectroscopically determined rotational period for WASP-85A, 14.2\thinspace $\pm$\thinspace 4\thinspace d, if the star's rotational axis is orthogonal to the line of sight. But, the value is also in agreement with the estimated rotational period for WASP-85B, 11.7\thinspace $\pm$\thinspace 4\thinspace d \citep{Brown14}. Therefore, it is uncertain which of the two binary components is causing the observed rotational modulation. However, because WASP-85A is active (see Figure~5), it is reasonable to assume that the rotational modulation is caused by the rotation of the spotted surface of WASP-85A. Because rotational modulations depend strongly on starspots coming and going, which causes the phase shifts, the period derived from it is therefore less accurate and less reliable than the period determined via starspot occultations (see Section~3.4).

We used Spot Oscillation And Planet 2.0 (SOAP 2.0, \citet{Dumusque14}) to fit the starspot rotational modulation of WASP-85 and found that the three characteristic ``M''-shaped modulation features  can be well fitted with two main starspot regions on the opposite sides of the rotating stellar surface. However, the change in the lightcurve over time shows that the starspot regions are changing over the K2 observations.

\begin{figure}
\includegraphics[width=8.5cm]{f6.png}
\caption{Autocorrelation function of the K2 Campaign 1 short-cadence lightcurve of WASP-85. Marked with \textit{P} is the first and highest peak corresponding to 13.6\thinspace d rotational period. Additional four peaks are visible at multiples of this period. Peaks corresponding to half period are caused by another starspot region on the opposite side of the star.}
\end{figure}

\subsection{Phase-curve variations}

We folded the K2 lightcurve on the orbital period of the planet. In order to do this we first removed the rotational modulation using our two-region model of the starspots. However, since this model is only an approximation, the folded lightcurve is still dominated by the residuals of the starspot modulation, which do not average out owing to the limited length of the K2 observations (Figure~7).

We did not detect the secondary eclipse nor any other phase-curve variations caused by the orbiting exoplanet. Using the system parameters from Table~1 the expected semi-amplitudes of ellipsoidal, Doppler, and reflection variations are 3.7, 2.1, and (for an albedo of 0.1) 15.4\thinspace ppm, respectively \citep{Pfahl08}. We conservatively estimate an upper limit of 50\thinspace ppm for the semi-amplitude of variations over the planet's phase-curve. Figure~7 shows a comparison between the measured phase-curve of WASP-85 and the simulated reflection modulation with a semi-amplitude of 50\thinspace ppm. The inability to detect the reflection modulation and the secondary eclipse in the measured K2 phase-curve of WASP-85 indicates that the planet's visual geometric albedo is lower than $\sim$0.3.

\begin{figure}
\includegraphics[width=8.5cm]{f7.png}
\caption{Measured phase-curve of WASP-85, binned to 200 bins, (shown in black) and a simulated reflection phase-curve modulation with a semi-amplitude of 50\thinspace ppm (red). The inability to see the reflection modulation and secondary eclipse in the measured phase-curve of similar amplitude as shown in the simulated phase-curve suggests that the actual reflection semi-amplitude is lower than 50\thinspace ppm.}
\end{figure}

\section{CONCLUSION AND DISCUSSION}

We have modified the PyKE self-flat-fielding method to nearly double its performance in restoring the degraded short-cadence photometric precision due to the pointing drift of the K2 spacecraft. We find similar system parameters to the discovery paper \citet{Brown14}, measure a rotational modulation period of 13.6\thinspace d, and reconstruct the ``M''-shaped photometric modulations with two dynamic starspot regions on either side of a magnetically active star. The presence of starspot occultations indicates that the host star WASP-85A is active and suggests that the stellar rotation and the planet's orbital axis are aligned. The phase shifts of potential occultation pairs correspond to a stellar rotational period of $15.1 \pm 0.6$\thinspace d. The agreement between the stellar rotation and transit-timing variation period suggests that its small semi-amplitude of 8\thinspace s is a consequence of stellar activity rather than an inter-planet gravitational interaction. Measured transit-duration variability is consistent with white noise distribution around zero. We do not detect any phase-curve modulations nor the secondary eclipse of the planet.

Starspots provide one method of constraining the sky-projected obliquity of the host star, while the Rossiter-McLaughin (R--M) effect provides an independent estimate.

For example, in WASP-4 the R--M effect implies a projected alignment angle of $\lambda = -4^{\circ +43}_{\,\,\, -34}$ \citep{Triaud10}, while \citet{Sanchis11b} report that on two occasions the same starspot was occulted, from which they constrain the alignment to $\lambda = -1^{\circ +14}_{\,\,\, -12}$. Similarly in WASP-6 the R--M angle is $\lambda = 11^{\circ +14}_{\,\,\, -18}$ \citep{Gillon09}, and \citet{Tregloan15} suggest that the same starspot was occulted in two transits and report $\lambda = 7.2 \pm 3.7^{\circ}$. In WASP-19 the R--M angle is $\lambda = 4.6 \pm 5.2^{\circ}$ \citep{Hellier11}, while \citet{Tregloan13} saw the same starspot occulted on consecutive nights and report $\lambda = 1.0 \pm 1.2^{\circ}$.

HAT-P-11 is misaligned, with an R--M angle of $\lambda =
103^{\circ +26}_{\,\,\, -10}$ \citep{Winn10}. The misalignment is
confirmed by the analysis of starspots by \citet{Sanchis11}.
Similarly, for Kepler-63 \cite{Sanchis13} report $\lambda =
-110^{\circ +22}_{\,\,\, -14}$ from analysing both starspots and the R--M
effect.

There is currently no measurement of the R--M effect in our
system, WASP-85A, but we predict that it would confirm that the system
is aligned. Qatar-2 is a similar case, where starspots indicate alignment \citep{Mancini14} though this is not yet confirmed by an R--M measurement.

There are no reported cases where the R--M angle conflicts with the analysis of starspots, and thus the starspots are a reliable way of measuring alignment in systems too faint for R--M measurements. Examples are Kepler-30 (\textit{V}\thinspace =\thinspace 15.6) where occultations of starspots by both Kepler-30c and Kepler-30d indicate alignment \citep{Sanchis12}, and Kepler-17 (\textit{V}\thinspace =\thinspace 14.3) where starspots also indicate alignment \citep{Desert11}.

The benefit of ongoing K2 observations of known WASP planets is that the systems are generally bright enough to enable good measurements of the alignment angle, $\lambda$, the stellar rotational period, the stellar radius, and the star's projected rotational broadening, $v \sin i$. Given all of these we obtain, not only the projected alignment $\lambda$, but the actual alignment angle $\psi$ (see, e.g., \citealt{Winn05}).

Adopting the rotational period of 15.1\thinspace $\pm$ 0.6\thinspace d, and combining with a stellar radius of 0.94\thinspace R$_{\odot}$ (Table~1) and a $v \sin i$ of 3.41 $\pm$ 0.89 km\thinspace s$^{-1}$ \citep{Brown14} we find that the star's rotational axis must be at $>$50$^{\circ}$ to the line of sight. This leads to the constraint that $\psi<41^{\circ}$.

The measured effective temperature and spin-orbit alignment of WASP-85 system are consistent with the current sample of predominantly aligned systems with stars below the characteristic effective temperature of 6250\thinspace K.

\acknowledgements{We would like to thank the anonymous referee for their comments which led to improving this paper. We gratefully acknowledge the financial support from the Science and Technology Facilities Council (STFC), under grants ST/J001384/1, ST/M001040/1 and ST/M50354X/1. This paper includes data collected by the K2 mission. Funding for the K2 mission is provided by the NASA Science Mission directorate. We thank Amaury Triaud for the use of RV data obtained with HARPS on the ESO 3.6-m telescope at La Silla (089.C-0151). This paper is based in part on observations made with the IO:O camera on the 2.0-m Liverpool Telescope (PL12B13). This work made use of PyKE \citep{Still12}, a software package for the reduction and analysis of \textit{Kepler} data. This open source software project is developed and distributed by the NASA Kepler Guest Observer Office.}

\bibliographystyle{apj}


\section{Problem Statement}
In this version of Colored Bin-Packing, we have $n$ items of unit weight, where each item is one of $x$ colors, where $x\geq 3$. Again, we have an unlimited supply of bins, each with weight limit $L$, in which to pack the items and our goal is to minimize the total number of bins.

\section{Algorithm}

The unit-weight case introduces weight constraints. Furthermore, we know that color constraints, specifically the discrepancy, forces us to use more bins than the weight constraint alone. The first half of the algorithm depends on discrepancy. If discrepancy is not a problem (i.e. is not more than 0), we simply order and split the bins up similar to the zero-weight case. If discrepancy is an issue (i.e. is more than 0), we handle the bins differently depending on if bin capacity is even or odd. In both cases we pack bins by alternating between items of the most frequent color and items of other colors. If bin capacity is even, we may be able to condense bins by using items of other colors that top some bins to merge bins containing only a single item of the most frequent color. If bin capacity is odd, we check if we can eventually reduce discrepancy to zero. This is because we use one additional max-color item for every bin we pack. If discrepancy eventually becomes zero, we are in the case where discrepancy is not a problem, and apply that algorithm. If this is not so, we will end up with bins topped by items of the most frequent color. The final set will be optimal, since all the bins will be topped with items of the this color and therefore cannot be combined.

The \textsc{Unit-Weight} algorithm solves this problem.

We first note that we must use at least $\bigg\lceil{\dfrac{n}{L}}\bigg\rceil$ number of bins.

Let $D$ denote the discrepancy, or the difference between the number of $MaxColor$ items and $OtherColors$ items. By the previous proof, we know that for 0-weight items, we will either need one bin (in the case where $D\leq 0$), or $D$ bins (in the case where $D > 0$).

\begin{algorithm}[ht] \caption{Algorithm \textsc{UnitWeight} ($I$). Input is a string $I$ of length $n$, representing items of different colors; $L$ is the weight limit of each bin}
\label{algunit}
\begin{algorithmic}[1]
\STATE Scan $I$ to count the number of times each color occurs. Also keep track of the color that occurs the most; let $MaxColor$ denote this color and let $MaxCount$ denote the number of items of $MaxColor$
\STATE Let $OtherColors$ denote the set of colors that are not $MaxColor$ and let $OtherCount$ denote the number of items of any color in $OtherColors$.
\STATE Let $D = MaxCount - OtherCount$ //Discrepancy
\STATE $S:$ solution //set of bins; initially empty
\IF{$D \le 0$}
	\STATE {\textsc{Split}}($I, n, L$) //\textsc{Split}() gets solution from \textsc{Zero-Weight()} and splits it into bins of capacity $L$
	\STATE Add the solution from {\textsc{Split}} to $S$.
    \STATE \textbf{return} $S$;
\ELSE
	\IF {$L$ is even}
		\STATE Create new current bin $C$. Start with placing a $MaxColor$ item, alternating with an $OtherColors$ item, until C is full. Add $C$ to $S$.
		\STATE Repeat the process until we run out of $OtherColors$ items. Make sure any partially full bins are topped with a $MaxColor$ item. There will be at most one such bin.
		\STATE Place the remaining $MaxColor$ items into bins of their own. Add these bins to $S$.
		\STATE Call  \textsc{Condense($S$)} to condense all the single-item bins. Add these bins to $S$.
	\ELSIF {$L$ is odd}
    \IF {$D \leqslant   \lceil \frac{OTHER}{\lfloor{L/2}\rfloor} \rceil $}
    \STATE //if there are enough OtherColors items such that alternating will allow the discrepancy to be reduced to zero

				\FOR{1 to $D$}
					\STATE Create new current bin $C$. Start with placing a $MaxColor$ item, alternating with an $OtherColors$ item, until C is full. Add $C$ to $S$.
				\ENDFOR
				\STATE Call {\textsc{Split}}($R$), where $R$ is the remaining items to be placed. //D=0 now so \textsc{Split} may be called
				\STATE Add the solution from {\textsc{Split}} to $S$.
			\ELSE
				\STATE Create new current bin $C$. Start with placing a $MaxColor$ item, alternating with an $OtherColors$ item, until C is full. Add $C$ to $S$.
				\STATE Repeat this until all $OtherColors$ items are used. Make sure any partially full bins are topped with a $MaxColor$ item. There will be at most one such bin.
				\STATE Place the remaining $MaxColor$ items in bins of their own. Add these to $S$.
			\ENDIF
	\ENDIF

\ENDIF
\end{algorithmic}
\end{algorithm}

\begin{algorithm}
\hspace{-1in}
\caption{Algorithm \textsc{Split} ($I, n, L)$}
\begin{algorithmic}[1]
	\STATE $R$ = \textsc{Zero-Weight}($I$):
	\FOR {$i = 1$ to $\bigg\lceil{\dfrac{n}{L}}\bigg\rceil$}
        \STATE Create new bin $C$
		\STATE Remove the first $L$ items (or remaining items if there are fewer than $L$) from $R$ and pack into $C$
        \STATE Add $C$ to $S$
	\ENDFOR
\end{algorithmic}
\end{algorithm}

\begin{algorithm}\caption{Algorithm \textsc{Condense} ($S$). $S$ is a set of bins, both completely and partially filled}
\begin{algorithmic}[1]

	\STATE Separate $S$ into three sets respectively called M-bins, P-bins and F-bins.
\STATE M-bins (MaxColor bins): contain only one item and the item is of MaxColor.
\STATE P-bins (partially full bins): contain both $MaxColor$ and $OtherColors$ items, are topped with a MaxColor item, and can fit at least two more items. There will be at most one of this bin type.
\STATE F-bins (full bins): fully packed bins that are topped with an $OtherColors$ item.
\STATE //Note that since \textsc{Condense} is called only when $L$ is even, all fully packed bins (F-bins) will be topped with an $OtherColors$ item, therefore all bins fall into one of these three categories.

        \IF {P-bins is not empty}
        		\STATE Let $current$ be the bin in P-bins
                \STATE Remove $current$ from P-bins
		\ELSIF {M-bins is not empty}
		\STATE Let $current$ be a bin in M-bins
		\STATE Remove $current$ from M-bins
	\ENDIF
	\WHILE {F-bins is not empty and M-bins is not empty}
		\IF {weight of $current$ $+ 2 > L$}
			\STATE Remove $current$ from its containing set
            \STATE Let $current$ be the next bin in M-bins  //switch to a new bin since $current$ is full
		\ENDIF
		\STATE Let $x$ be the top-most $OtherColors$ item from a bin in F-bins
        \STATE Let $y$ be a $MaxColor$ item from a bin (or the bin??) in M-bins
        \STATE Pack $x$ followed by $y$ into $current$.
		\STATE Remove $x$ and $y$ from their containing bins.
	\ENDWHILE
\end{algorithmic}
\end{algorithm}

\newpage

\section{Proof}
\textbf{Case 1:} $D \leqslant 0$ (first if statement) If the discrepancy is not more than zero, we use the \textsc{Zero-Weight} algorithm to first order the items so they satisfy the color requirements. Then, without breaking this order, we pack the first $L$ items in the first bin, the next $L$ items in the second, and so on, until we have packed all the items.
\\
\\
\\
\\
\noindent \textbf{Example:}\\
\underline{Items:} 4 W, 3 B, 2 Y\\
\underline{Weight Limit:} 3\\
In the zero-weight case, we would sort these items in this order, using only one bin:\\
\begin{center}BWBWBWYWY\\
\end{center}
We use this result to get the unit weight solution:\\
\begin{center}BWB WBW YWY\\
\end{center}
where a space denotes a new bin.
\\
\\
This is based on the two requirements defined by the problem: color and weight. We know from the zero-weight case that the \textsc{Zero-Weight} algorithm satisfies the color restriction when discrepancy is less than or equal to 0. By allocating the items according to the weight limit, we use the minimum amount of bins necessary based on bin capacity. Since both are satisfied, the solution is optimal.
\\
\\
\\
\textbf{Case 2:} $D > 0$ (else statement)
\begin {enumerate}
\item {When $L$ is even.}

\noindent \textbf{Example:}\\
\underline{Items:} 12 W, 3 B, 2 Y 2 G\\
\underline{Weight Limit:} 4\\
\underline{D:}  12 - (3 + 2 + 2) = 12 - 7 = 5\\
We use  {\textsc{Unitweight}} to first pack these into bins,\\
\begin{center} WBWB  WYWG  WBWY  WGW  W  W  W  W\\
\end{center}
where a space denotes a new bin. Note that there are four bins that only contains a single  $MAXCOLOR$ - white - item.
We reduce the number of bins further by calling {\textsc{Condense}}.\\
This becomes,
\begin{center}WBW  WYW  WBWY  WGW  WBW  WGW\\
\end{center}
giving us 6 bins instead of 8.
\\
\\
Note that we use the same definition of N-bins, H-bins and F-bins seen in {\textsc{Condense}}.

With a positive discrepancy, it is easy to see that the color requirement will force us to have more bins than the minimum number of bins enforced by the weight limit. This is because each excess $MaxColors$ item requires an individual bin - these bins form the set of `M-bins'.  \\

Keeping this in mind, we also want to use $OtherColors$ items to buffer as many $MaxColor$ items as possible. We do this to minimize the effect excess $MaxColor$ items will have on the eventual number of bins. The more $MaxColor$ items we buffer, the fewer M-bins we are forced to create. Hence we start every bin with a $MaxColor$ item, and place only one $OtherColors$ item between each $MaxColor$ item. Completely packed bins topped with an $OtherColors$ item form the set called 'F-bins'. \\

When we run out of $OtherColors$ items while packing a bin, we get a potential P-bins. To use up as many $MaxColor$ items as possible, we cover the $OtherColors$ item currently topping the bin, with a $MaxColor$ item. If the bin is able to accept 2 or more items after we do so, we get a P-bins. Because this happens only when we run out of $OtherColors$ items, a P-bins only occurs once.\\

We refer to the packing before the call to \textsc{Condense} as the \textit{initial packing}. After the initial packing, $S$ contains F-bins, M-bins and a possible P-bins. Because $L$ is even and all F-bins start with a $MaxColor$ item, all F-bins must be topped with an $OtherColors$ item. These $OtherColors$ items are used to combine the M-bins, and further reduce our total number of bins using the algorithm {\textsc{Condense}}. \\

To show {\textsc{Condense}}'s optimality, we prove that the number of bins in $S$ cannot be further reduced after {\textsc{Condense}} is run.\\\\
There are several cases, all of which are dependent on {\textsc{Condense}}'s exit condition.

\begin{enumerate}[i.]
\item $S$ contains no F-bins.\\

\item $S$ contains no M-bins.\\

\item $S$ contains neither F-bins nor M-bins.
\end{enumerate}
Since combining bins requires both F-bins and M-bins, any one of these three conditions means that we cannot reduce the total number of bins by combining any bins in $S$.
\item {When $L$ is odd.}

We observe that when the weight limit is odd, if we use as many $MaxColor$ items as possible when packing a bin, we use one more $MaxColor$ than $OtherColors$ item for every bin. To do so, we fill the bin in a similar fashion to if $L$ is even - starting with a $MaxColor$ item and alternating with $OtherColors$ items. However, because the $L$ is odd, we end with a $MaxColor$ item. This is the excess $MaxColor$ item we pack when packing each bin.\\

Next, we note that if $D > 0$, the number of M-bins depends on the number of $OtherColors$ items. This is because since we alternate between $OtherColors$ and $MaxColor$ items when packing the bins, each bin takes exactly  $\lfloor {L/2} \rfloor$ $OtherColors$ items. If $OtherCount$ is the total number of $OtherColors$, we can form $\lceil \frac{OtherCount}{\lfloor{L/2}\rfloor} \rceil$ P-bins and F-bins. The remaining $MaxColor$ items must each be packed in a separate bin which are the M-bins.\\

Thus, if we can form enough P-bins and F-bins to pack the excess $MaxColor$ items, we can reduce $D$ to 0, otherwise we form M-bins. This gives us 2 cases.

\begin{enumerate}
\item{We eventually reduce $D$ to 0.}

\noindent \textbf{Example:}\\
\underline{Items:} 8 W, 3 B, 2 Y 2 G\\
\underline{Weight Limit:} 5\\
\underline{D:}  8 - (3 + 2 + 2) = 8 - 7 = 1\\
We use {\textsc{Unitweight}} to first pack these into bins.\\
After the first bin,
\begin{center} WBWBW\\
\end{center}
where a space denotes a new bin, {\textsc{Discrepancy}} is now zero and we pack the remaining items using {\textsc{Split}}.
We first sort the remaining items,\\
\begin{center}WYWGWBWYWG\\
\end{center}
before splitting the sorted items into bins to give us,
 \begin{center}WYWGW  BWYWG\\
\end{center}
This gives us,
\begin{center}WBWBW  WYWGW  BWYWG\\
\end{center}
as our overall solution.
\\
\\

If $D \leqslant   \lceil \frac{OtherCount}{\lfloor{L/2}\rfloor} \rceil $, we know that $D = 0$ after $D$ bins are filled in Step -- (fill in). From Case 1, we know that {\textsc{Split}} gives us an optimal packing when $D = 0$. Thus, for $S$ to be optimal, the $D$ bins packed during this phase must also be optimal.

To show that the $D$ bins are optimally packed, we note that the initial packing falls into two categories:\\

\begin{enumerate}
\item {Initial packing yields a P-bin:}

If the initial packing yields a P-bin, it means that $D = 1$ when we started filling the P-bin. The phase ends when we place the last excess $MaxColor$ item into the P-bin. This reduces $D$ to 0.\\

In this case, we show that there are no more items left. \\

Since $D = 0$, we know the number of $MaxColor$ and $OtherColors$ items must be equal. Any additional $OtherColors$ item would mean the P-bins would not exist, as the last item placed into the bin would be an $OtherColors$ item. This is a contradiction. Since there are no $OtherColors$ items remaining there must be no $MaxColor$ items remaining. Since there are no items remaining,  we cannot combine any of the bins, and the solution set is optimal.\\

\item {Initial packing does not yield any P-bins.}

Note that we will not have any M-bins. When $D=0$, there are two cases:

\begin{enumerate}
\item $MaxCount = OtherCount = 0$:

We won't need any additional bins so clearly there are no M-bins.

\item $MaxCount = OtherCount > 0$:
Since there are no excess $MaxColor$ items, there will be no bins containing only a $MaxColor$ item.

Here, the solution set contains only F-bins. Since we have no M-bins, bin combination cannot happen, and the set is optimal. \\
 \end{enumerate}

 \item{$D$ remains greater than 0.}

\noindent \textbf{Example:}\\
\underline{Items:} 15 W, 3 B, 2 Y 2 G\\
\underline{Weight Limit:} 5\\
\underline{D:}  15 - (3 + 2 + 2) = 15 - 7 = 6\\
We use {\textsc{Unitweight}} to pack these into bins, knowing that $D$ will not eventually reduce to zero.
\begin{center} WBWBW  WYWGW  WBWYW  WGW  W  W  W  W\\
\end{center}
where a space denotes a new bin. Though there are bins containing only one items, we cannot condense them since all the non-$MAXCOLOR$ items need them in their bins. The solution is optimal.\\
\\
\\

This means we did not have enough P-bins and/or F-bins to pack excess $MaxColor$ items and therefore cannot reduced $D$ to 0. $S$ contains at least one M-bins. However, because all the bins are topped by $MaxColor$ items, we cannot combine any bins to reduce the total number of bins. $S$ is optimal.
\end{enumerate}
\end{enumerate}

\end{enumerate}

\section{Time Complexity}
We analyze the time complexity based on the various cases.
If $D$ is less than or equal to 0, the algorithm first calls \textsc{Split} which runs \textsc{Zero-Weight} to order the items based on color, then packs the items based on weight, using at $O(n)$ time in total where $n$ is the number of items.

If $D$ is more than 0, we consider two sub-cases: $L$ is odd or $L$ is even.
If $L$ is odd, we check if we can eventually reduce the $D$ to 0. If we can, we pack items until $D$ is 0. This takes at most $O(n)$, since we consider each item only
once. When $D$ is 0 we call \textsc{Split} which takes $O(n)$ time.  If we cannot
reduce discrepancy to 0, we pack items by alternating colors. Since each item is considered exactly once, this takes $O(n)$ time.

If $L$ is even, again we pack by alternating colors then call \textsc{Condense}. The runtime of \textsc{Condense} depends on the
number of bins and since there will be $O(n)$ bins, \textsc{Condense} runs in $O(n)$ time.

\bigskip

\section{Bin Coloring Problem with 0-weight items of more than 2 colors, reordering allowed}

In this section, we consider a version of the Bin Coloring Problem where items are of more than 2 colors, all items have a weight of 0, and reordering is allowed. More formally, the problem is defined as follows:

\begin{itemize}
\item The input to the problem is a string. Each character $x$ represents an item with the color $x$.
\item All items have a weight of 0.
\item No adjacent items in the bin can be of the same color.
\item Reordering is allowed. For example, for $i$ and $j$ such that $i>j$, the $i$th item in the string can be packed either before or after the $j$th item in the string.
\end{itemize}

We provide an algorithm \textsc{zero-weight} that is optimal for this problem.

Since the only constraint in the 0-weight situation is the color constraint (that no two adjacent items can have the same color), we pay special attention to the most frequent color, i.e. the color with the most number of items. In the first case, if this color has no more items than all the other colors combined, we can pack all the items in one bin, using other color items as buffers between items of the most frequent color. Specifically, we first alternate between items of other colors until the we have exactly enough items of other colors to buffer the items of the most frequent color. Then, we alternate between items of other colors and items of the most frequent color. If on the other hand, the most frequent color occurs more than all the other colors combined, we are forced to use all the other color items as buffers. Even so, we will be running out of the other color items in the end and left with several items of the most frequent because there are too many of them. For each of the items of the  most frequent color that are remaining, we are forced to open one new bin, since they are all of the same color.

\begin{algorithm}\caption{Algorithm \textsc{zero-weight}. Input is a string $I$, representing items of different colors. }
\label{algo1}
\begin{algorithmic} [1]
\STATE Scan $I$ to count the number of times each color occurs. Also keep track of the color that occurs the most; let $MaxColor$ denote this color and let $MaxCount$ denote the number of items of $MaxColor$
\STATE Let $OtherColors$ denote the set of colors that are not $MaxColor$ and let $OtherCount$ denote the number of items of any color in $OtherColors$.
\STATE $Discrepancy = MaxCount - OtherCount$
\IF {$Discrepancy \le 0$}
\STATE Create a new bin $B$ whose size is equal to the length of $I$.
\STATE Alternate between an item of $MaxColor$ and an item of any arbitrary $OtherColors$ until $B$ is full and all items are packed. Only one bin is needed in this case.
\STATE \textbf{return} $B$
\ELSE
\STATE Create an empty set $S$ that will be the solution
\STATE Create a bin $B_0$ of size $OtherCount * 2 + 1$. Add $B_0$ to $S$.
\STATE Alternate between an item of $MaxColor$ and an item of any arbitrary $OtherColors$ until $B_0$ is full. When $B_0$ is full, there will be no items of $OtherColors$ left.
\STATE Create a bin $B_x$ for each of the remaining $x$ items of $MaxColor$ and pack each of these items in one of these bins. In this case, the total number of bins is $Discrepancy$. Add each $B_x$ to $S$
\STATE \textbf{return} $S$
\ENDIF
\end{algorithmic}
\end{algorithm}

\begin{proof}

First note that from a previous paper (Bin Packing with colors), if there are only two colors, then the minimum number of bins necessary is the discrepancy of the two colors. For example, if there are 5 Black items and 3 White items, the number of bins necessary is $5-3=2$ bins. This is done by starting with an item of the color with more items ($B$) and alternating with items of the other color ($W$) until $W$ runs out. We then add one more item of color $B$ into this first bin. If $|B$ and $|W|$ denote the number of black and white items, respectively, we need $|B|-|W|-1$ more bins to fit the rest of the items of color $W$ for a total of $|B|-|W|$ bins.

For more than 2 colors, we consider two cases:

\underline{Case 1}: $Discrepancy \le 0$:
See below for an example.
In this case there are at least as many $OtherColors$ items as there are $MaxColor$ items. In this case, we need only one bin, which is the fewest number of bins any algorithm can get. To do this, we alternate between items of $OtherColors$ until the sum of the unused $OtherColors$ items is one less than the number of $MaxColors$ items. We then add one $MaxColor$ item to the bin. From here, we alternate between any $OtherColors$ item and a $MaxColor$ item until all items are used. The last item will be one of $MaxColor$.

Note that it is always possible to alternate the $OtherColors$ items until $OtherCount$ is one less than the number of $MaxColor$ items. The only case when we can no longer alternate is when there are items of only one color, $X$, remaining, and because of the color restriction, we can no longer put these items into the bin. We claim that at this point, $OtherCount$, which in this case is the number of $X$ items, will be one less than $MaxCount$. We prove this by way of contradiction. Suppose when we stop alternating, $OtherCount$ is equal or greater than $MaxCount$. Note that an $X$ item must be on top of the bin, otherwise we have not completed the alternating step since we can pack another $X$ item on top of the bin. There must be at least $MaxCount$ $X$ items remaining, so there were at least $MaxCount + 1$ total $X$ items, which contradicts our assumption that $MaxColor$ is the color with the maximum number of items.

After these alternations, we have one more $MaxColor$ items than all $OtherColors$ items, and it is easy to see that we can alternate between them and put all items in the same bin.

\underline{Case 2}: $Discrepancy > 0 $. In this case, all other items can be regarded as some uniform color \textit{non-MaxColor} that is not equal to $MaxColor$, which makes this the black/white bin packing problem. This is because there are so many $MaxColor$ items that we must pack each \textit{non-MaxColor} item between two $MaxColor$ items to use as few bins as possible. The first bin $B_0$ will start with a \textit{non-MaxColor} item, alternate $MaxColor$ and \textit{non-MaxColor} items, and end with a \textit{non-MaxColor} item, so will have $OtherCount$*2+1 items in total. Each additional bin will have an item of $MaxColor$. Then the minimum number of bins necessary is still the difference between the number of $MaxColor$ and \textit{non-MaxColor} items, i.e. the discrepancy, as is shown in the Black/White bin packing problem.

\end{proof}

The first step of the algorithm requires a scan of the input string to determine $OtherCount$ and $MaxCount$. The second step requires a scan to separate the items into in two sets: items of $OtherColors$ and items of $MaxColor$.  Packing each item takes constant time, so the remainder of the algorithm runs linear to the number of items. Therefore the algorithm has total runtime $O(3n)$ where $n$ is the number of items.

\bigskip
\bigskip

\noindent \textsc {Example: }

\medskip

\noindent \underline {Items:} WWWWWWWWBBYY (8W, 2B, 2Y)

\medskip

\noindent In this case, $MaxColor$ is $W$. $MaxCount$ equals $8$ and $OtherCount$ equals $4$. $Discrepancy > 0$.

\medskip

\noindent According to our algorithm, we first create a bin $B_0$ of size $4*2+1=9$. Then we alternate between $W$ and any $OtherColor$ to fill the first bin.

\medskip

\begin{center} WBWBWYWYW
\end{center}

\medskip

\noindent Now, we have used all the $OtherColors$, and we still have 3 $W$'s remaining. Therefore, we create $3$ new bins, each holding one of the remaining $W$'s.

\medskip

\noindent Therefore, we get the following final result:

\medskip

\begin{center} WBWBWYWYW / W / W / W
\end{center}

\end{document}


\section{Introduction}
The Anderson impurity model is one of the most important and fundamental
model Hamiltonians in condensed-matter physics~\cite{Anderson:1961}.  It
was originally invented to study localized magnetic impurities in
metallic environments and there exist three distinct parameter regimes
called empty orbital regime, mixed valence regime, and local moment
regime~\cite{Hewson:1993}.  Of particular interest is the local moment
regime, where the electrical resistivity grows logarithmically toward
the low temperature as a consequence of the celebrated Kondo
effect~\cite{Kondo:1964,Schrieffer:1966}.  By extending the Anderson
impurity model, orbital degeneracy can also be
incorporated~\cite{Coqblin:1969} and periodic arrangement of magnetic
impurities is considered to be relevant to heavy fermion
physics~\cite{Coleman:2007}.

Further application of the Anderson impurity model is possible by
introducing additional degrees of freedom corresponding to left and
right leads to study tunneling of electrons through a quantum dot both
in and out of equilibrium~\cite{Meir:1991,Meir:1993}.  In fact, one
degree of freedom can be decoupled by a canonical transformation and the
other interacting degree of freedom turns out to be described by the
original Anderson impurity model~\cite{Glazman:1988}.  Therefore, the
Kondo effect emerges again in quantum dot systems but now as the
logarithmic growth of the electrical conductance by lowering the
temperature~\cite{Glazman:1988,Ng:1988}, where a number of beautiful
observations of the Kondo effect have been
made~\cite{Grobis:2007,Wiel:2010}.

The purpose of this Rapid Communication is to show that all the rich
physics associated with the Anderson impurity model can be simulated
with ultracold atoms by employing standard techniques such as
magnetic-field-induced Feshbach resonances, species-selective optical
lattices, and laser couplings of atomic hyperfine
states~\cite{Lewenstein:2012}.  In particular, we place special emphasis
on the transport measurement of the Kondo effect analogous to quantum
dot experiments, which should be of great importance in ultracold atom
experiments because transport is usually difficult to study with the
exception of recent progress made in
Refs.~\cite{Brantut:2012,Stadler:2012,Krinner:2013,Brantut:2013,Krinner:2015,Lee:2015,Husmann:2015}.
Future realization of the Kondo effect and its transport measurement
with ultracold atoms will pave the way for developing new insights into
yet unresolved aspects of Kondo physics such as the formation and
dynamics of the Kondo screening cloud~\cite{Kouwenhoven:2001} and the
quantum criticality in heavy fermion systems~\cite{Coleman:2007}.

\section{Setup and measurement protocol}
Our study is based on the simple and versatile scheme to realize the
orbital Kondo effect with ultracold atoms~\cite{Nishida:2013} (see
Refs.~\cite{Recati:2002,Falco:2004,Duan:2004,Paredes:2005,Gorshkov:2010,Lal:2010,Bauer:2013,Kuzmenko:2015,Sundar:2015,Nakagawa:2015,Patton:2015}
for other proposals).  The system consists of a Fermi sea of
spin-polarized $\up$ fermions of species $A$ interacting with a spinless
impurity atom of different species $B$ which is loaded into a ground
state of an isotropic potential.  By tuning the interspecies attraction
with an $s$-wave Feshbach resonance, the impurity atom and a
spin-polarized fermion can form a bound molecule that occupies a
degenerate orbital of the confinement potential with orbital angular
momentum $\ell\geq1$~\cite{Nishida:2010}.  In particular, when the total
energy of the bound molecule coincides with the scattering threshold of
the $A_\up$ and $B$ atoms, an $\ell$th partial-wave resonance is induced
and low-energy physics in its vicinity is described by a two-channel
Hamiltonian:
\begin{align}\label{eq:H_up}
 H_\up &= \int\!\frac{d\k}{(2\pi)^3}\,\ek\,\psi_{A\up}^\+(\k)\psi_{A\up}(\k)
 + \sum_{m=-\ell}^\ell\delta_m\phi_m^\+\phi_m \notag\\
 & + \sum_{m=-\ell}^\ell\int\!\frac{d\k}{(2\pi)^3}
 \left[V_\ell^m(\k)\psi_{A\up}^\+(\k)\psi_B^\+\phi_m + \text{H.c.}\right].
\end{align}
Here $\psi_{A\up}^\+(\k)$ creates an $A_\up$ atom with energy
$\ek=\hbar^2\k^2/(2M)$, while $\psi_B^\+$ creates the impurity $B$ atom
in the ground state of the confinement potential whose energy is chosen
to be zero.  The bound molecule is created by $\phi_m^\+$ in one of the
degenerate orbitals labeled by the magnetic quantum number
$|m|\leq\ell$.  Its coupling to the $A_\up$ and $B$ atoms is assumed to
have the harmonic form of
\begin{align}\label{eq:coupling}
 V_\ell^m(\k) = v_m|\k|^\ell Y_\ell^m(\hat\k)
 \exp\!\left[-\frac{\k^2}{2\Lambda^2}\right]
\end{align}
with the wave-number cutoff $\Lambda$ set by an inverse characteristic
extent of the confined $B$ atom.  Because only one $B$ atom is confined,
the particle number operators of the localized $B$ atom and bound
molecule are constrained by
\begin{align}\label{eq:constraint}
 N_B = \psi_B^\+\psi_B + \sum_{m=-\ell}^\ell\phi_m^\+\phi_m = 1.
\end{align}
When the rotational symmetry is exact, we have an equal detuning
$\delta_m$ and coupling $v_m$ for all $m$ and thus the degeneracy is
($2\ell+1$)-fold, while we shall develop a general formulation so that
it is also applicable to study the effect of symmetry breaking later.

Interestingly, the low-energy effective Hamiltonian (\ref{eq:H_up})
naturally realizable with ultracold atoms is nothing but the
infinite-$U$ Anderson impurity model in the slave-particle
representation~\cite{Barnes:1976,Coleman:1984} with its fictitious
degrees of freedom corresponding to our real atom and molecule as in
Eq.~(\ref{eq:constraint}).  The empty orbital, mixed valence, and local
moment regimes of the original Anderson impurity model are thus
translated into atomic regime ($\<\psi_B^\+\psi_B\>\simeq1$), resonant
regime ($0\lesssim\<\psi_B^\+\psi_B\>\lesssim1$), and molecular regime
($\<\psi_B^\+\psi_B\>\simeq0$), respectively, in the language of
ultracold atoms.  In particular, the orbital Kondo effect emergent in
the molecular limit was elaborated for $\ell=1$ in
Ref.~\cite{Nishida:2013}, while the analysis therein can be
straightforwardly generalized for an arbitrary $\ell$ to find that the
Kondo temperature in the SU($2\ell+1$) symmetric case has a universal
leading exponent given by
\begin{align}\label{eq:T_kondo}
 \TK \propto \TF\exp\!\left[-\frac{\pi}{(2\ell+1)a_\ell\kF^{2\ell+1}}\right]
\end{align}
with $a_\ell\ll\kF^{-2\ell-1}$ being the $\ell$th partial-wave
scattering length.  Because the Kondo effect in condensed-matter systems
manifests itself most sharply in their transport measurements, it is
highly desired although challenging to establish an analogous transport
signature of the orbital Kondo effect in ultracold atom experiments.

In what follows, we indeed show that the conductance measurement of the
Kondo effect in quantum dot experiments can be equivalently performed
with ultracold atoms by adopting the idea from Ref.~\cite{Knap:2012}.
To this end, we introduce another spin $\down$ component of fermionic
species $A$,
\begin{align}\label{eq:H_down}
 H_\down = \int\!\frac{d\k}{(2\pi)^3}\,\ek\,\psi_{A\down}^\+(\k)\psi_{A\down}(\k),
\end{align}
as well as the intercomponent coupling driven by a resonant laser field
with the Rabi frequency $\Omega$,
\begin{align}\label{eq:H_up-down}
 H_{\up\down} = i\frac{\hbar\Omega}{2}\int\!\frac{d\k}{(2\pi)^3}
 \left[\psi_{A\up}^\+(\k)\psi_{A\down}(\k) - \psi_{A\down}^\+(\k)\psi_{A\up}(\k)\right],
\end{align}
which are expressed in the rotating frame.  It is legitimate to assume
that interactions of $A_\down$ atoms with $A_\up$ atoms and with the
impurity $B$ atom are both negligible in the dilute limit because they
are generally off-resonance when the interaction of $A_\up$ atoms with
the impurity $B$ atom is tuned close to a confinement-induced resonance.
It is also important that the confinement-induced resonance can be
turned on and off without changing the magnetic field but by controlling
the potential strength acting on the impurity $B$
atom~\cite{Nishida:2010}.  With all these setups, we are ready to
propose a simple conductance measurement with ultracold atoms consisting
of the following three protocols.

(i) Preparation: We first stay away from any confinement-induced
resonances so that both $A_\up$ and $A_\down$ atoms negligibly interact
with the impurity $B$ atom.  After introducing $N_-$ and $N_+$ numbers
of $A_\up$ and $A_\down$ atoms, we apply the intercomponent coupling
(\ref{eq:H_up-down}) for a duration of $\pi/(2\Omega)$, which transforms
the two spin components into two superposition states according to
$|A_\up\>\to|A_-\>\equiv(|A_\up\>-|A_\down\>)/\sqrt2$ and
$|A_\down\>\to|A_+\>\equiv(|A_\up\>+|A_\down\>)/\sqrt2$.  Consequently,
$N_-$ and $N_+$ numbers of $A_-$ and $A_+$ atoms are prepared at a
common temperature $T$ and their corresponding chemical potentials are
denoted by $\mu_-\equiv\mu-\Delta\mu/2$ and
$\mu_+\equiv\mu+\Delta\mu/2$, respectively.

(ii) Transport: We now turn on the confinement-induced $\ell$th
partial-wave resonance of $A_\up$ atoms with the impurity $B$ atom.
Because $A$ atoms are prepared on the $\pm$ basis, it is appropriate to
express the total Hamiltonian composed of Eqs.~(\ref{eq:H_up}) and
(\ref{eq:H_down}) in terms of the corresponding creation operators
$\psi_{A\pm}^\+(\k)\equiv[\psi_{A\up}^\+(\k)\pm\psi_{A\down}^\+(\k)]/\sqrt2$
as
\begin{align}\label{eq:transport}
 & H_\up + H_\down = \sum_{\sigma=\pm}\int\!\frac{d\k}{(2\pi)^3}\,\ek\,
 \psi_{A\sigma}^\+(\k)\psi_{A\sigma}(\k)
 + \sum_{m=-\ell}^\ell\delta_m\phi_m^\+\phi_m \notag\\
 & + \sum_{\sigma=\pm}\sum_{m=-\ell}^\ell\int\!\frac{d\k}{(2\pi)^3}
 \left[\frac{V_\ell^m(\k)}{\sqrt2}\psi_{A\sigma}^\+(\k)\psi_B^\+\phi_m
 + \text{H.c.}\right],
\end{align}
where $A_-$ and $A_+$ atoms equally interact with the impurity $B$
atom.  Remarkably, the resulting Hamiltonian is identical to the
Anderson impurity model extended to study tunneling of electrons through
a quantum dot~\cite{Meir:1991,Meir:1993} with roles of left and right
leads played by our two superposition states.  The chemical potential
imbalance $\Delta\mu\neq0$ thus causes a transport of fermions from
majority to minority superposition states through scatterings with the
impurity $B$ atom~\cite{transport}.  After a period of $\Delta t$, the
numbers of $A_-$ and $A_+$ atoms change into $N_-+I\Delta t$ and
$N_+-I\Delta t$ with $I$ being the steady-state current.

(iii) Measurement: We then turn off the confinement-induced resonance to
terminate the transport.  By applying the intercomponent coupling
(\ref{eq:H_up-down}) again but for a duration of $3\pi/(2\Omega)$, the
two superposition states are transformed back into the original two spin
components as $|A_-\>\to-|A_\up\>$ and $|A_+\>\to-|A_\down\>$ whose
numbers are now $N_-+I\Delta t$ and $N_+-I\Delta t$, respectively.
Finally, from the measured numbers of $A_\up$ and $A_\down$ atoms
compared to their initial values, the transported fermion number
$I\Delta t$ can be extracted to determine the conductance $I/\Delta\mu$
of the transport Hamiltonian (\ref{eq:transport}).  This constitutes the
ultracold atom equivalent of the conductance measurement in quantum dot
systems.

\section{Conductance and the Kondo effect}
In order to provide quantitative guides on how the orbital Kondo effect
emerges in the proposed conductance measurement with ultracold atoms, we
study in detail the linear conductance
$G=\lim_{\Delta\mu\to0}I/\Delta\mu$ of the transport Hamiltonian
(\ref{eq:transport}).  With the aid of the Meir-Wingreen formula for the
steady-state current $I=(dN_-/dt-dN_+/dt)/2$~\cite{Meir:1992}, the
linear conductance is given by
\begin{align}\label{eq:conductance}
 G = \frac1\hbar\sum_{m=-\ell}^\ell\int\!\frac{d\k}{(2\pi)^3}
 \frac{|V_\ell^m(\k)|^2}{2}f'_T(\ek-\mu)\,\Im\,\G_m^R(\ek-\mu),
\end{align}
where $f_T(z)=1/(e^{\beta z}+1)$ is the Fermi-Dirac distribution
function at an inverse temperature $\beta=1/(\kB T)$ and
\begin{align}\label{eq:retarded}
 \G_m^R(\eps) = -\frac{i}{\hbar}\int_0^\infty\!dt\,e^{i\eps t/\hbar}
 \<\{\psi_B^\+(t)\phi_m(t),\phi_m^\+(0)\psi_B(0)\}\>
\end{align}
is the retarded Green's function with the expectation value taken with
respect to the equilibrium state at $\Delta\mu=0$.  Therefore, $A_-$ and
$A_+$ atoms now have the equal chemical potential $\mu_\pm=\mu$ and thus
it is advantageous to express the transport Hamiltonian
(\ref{eq:transport}) on the original spin basis so that it is decoupled
into Eqs.~(\ref{eq:H_up}) and (\ref{eq:H_down})~\cite{decoupling}.
Because $A_\down$ atoms do not interact with the impurity $B$ atom, the
expectation value is simply evaluated as
$\<\cdots\>=\Tr[e^{-\beta(H_\up-\mu N_\up)}\cdots]/Z$ along with the
constraint (\ref{eq:constraint}), where
\begin{align}
 N_\up = \int\!\frac{d\k}{(2\pi)^3}\psi_{A\up}^\+(\k)\psi_{A\up}(\k)
 + \sum_{m=-\ell}^\ell\phi_m^\+\phi_m
\end{align}
is the particle number operator of $A_\up$ atoms.  While the equilibrium
Green's function (\ref{eq:retarded}) can be computed by means of various
methods~\cite{Georges:1996}, we here employ the so-called noncrossing
approximation~\cite{Bickers:1987,Hewson:1993}, which is known to be
reliable for large degeneracy and not too low temperature, and thus
adequate for our purpose.

\begin{figure}[t]
 \includegraphics[width=\columnwidth]{noncrossing}
 \caption{Noncrossing approximation for self-energies of the localized
 (a) $B$ atom [Eq.~(\ref{eq:self-energy_B})] and (b) bound molecule
 [Eq.~(\ref{eq:self-energy_m})] as well as (c) the Matsubara Green's
 function corresponding to Eq.~(\ref{eq:retarded}).  Solid, dashed, and
 doubled lines represent the propagators of $A_\up$ atom, $B$ atom, and
 bound molecule, respectively.  \label{fig:noncrossing}}
\end{figure}

The propagators of the localized $B$ atom and bound molecule are
sometimes called resolvents and denoted by
$\R_B(z)=[z-\Sigma_B(z)]^{-1}$ and
$\R_m(z)=[z+\mu-\delta_m-\Sigma_m(z)]^{-1}$, respectively.  Their
self-energies within the noncrossing approximation are depicted in
Fig.~\ref{fig:noncrossing} and determined self-consistently according to
\begin{align}\label{eq:self-energy_B}
 \Sigma_B(z) = \sum_{\ell=-m}^m\int\!\frac{d\k}{(2\pi)^3}\,
 |V_\ell^m(\k)|^2\,f_T(\ek-\mu)\,\R_m(z+\ek-\mu)
\end{align}
and
\begin{align}\label{eq:self-energy_m}
 \Sigma_m(z) = \int\!\frac{d\k}{(2\pi)^3}\,
 |V_\ell^m(\k)|^2\,f_T(\mu-\ek)\,\R_B(z-\ek+\mu).
\end{align}
In terms of the corresponding spectral densities of
$\rho_B(\eps)=-(1/\pi)\,\Im\,\R_B(\eps+i0^+)$ and
$\rho_m(\eps)=-(1/\pi)\,\Im\,\R_m(\eps+i0^+)$, the imaginary part of the
retarded Green's function (\ref{eq:retarded}) is expressed as
\begin{align}\label{eq:imaginary}
 -\frac1\pi\,\Im\,\G_m^R(\eps)
 = \frac{1+e^{-\beta\eps}}{Z_B}\int_{-\infty}^\infty\!dz\,
 e^{-\beta z}\rho_B(z)\rho_m(z+\eps)
\end{align}
with
\begin{align}
 Z_B = \int_{-\infty}^\infty\!dz\,e^{-\beta z}
 \left[\rho_B(z)+\sum_{m=-\ell}^\ell\rho_m(z)\right]
\end{align}
being the impurity partition function~\cite{Bickers:1987,Hewson:1993}.
The substitution of Eq.~(\ref{eq:imaginary}) into
Eq.~(\ref{eq:conductance}) now allows us to compute the linear
conductance numerically for a given set of parameters.

Besides the chemical potential $\mu$ and temperature $T$ of $A_\up$
atoms, the linear conductance (\ref{eq:conductance}) depends on the
detuning $\delta_m$ and coupling $v_m$ as well as the wave-number cutoff
$\Lambda$ through $V_\ell^m(\k)$ defined in Eq.~(\ref{eq:coupling}).  In
order to make contact with ultracold atom experiments, the bare
parameters $\delta_m$ and $v_m$ should be expressed in terms of physical
parameters such as the scattering length $a_m$ and the resonance range
$r_m$ characterizing low-energy scatterings in the $\ell$th partial-wave
channel.  They can be related by matching the two-body scattering
$\T$ matrix in the vacuum computed from the two-channel Hamiltonian
(\ref{eq:H_up}) with the standard form of
\begin{align}\label{eq:T-matrix}
 \T_\ell(k) = \frac{8\pi^2\hbar^2}{M}\sum_{m=-\ell}^\ell
 \frac{k^{2\ell}\,Y_\ell^m(\hat\k_\mathrm{out})\bar{Y}_\ell^m(\hat\k_\mathrm{in})}
 {ik^{2\ell+1}+1/a_m+r_mk^2+O(k^4)},
\end{align}
where we find
\begin{align}
 \frac{1}{a_m} = -\frac{8\pi^2\hbar^2\delta_m}{Mv_m^2}
 + \frac{\Gamma(\ell+\frac12)\Lambda^{2\ell+1}}{\pi}
\end{align}
and
\begin{align}
 r_m = \frac{1}{a_m\Lambda^2} + \frac{4\pi^2\hbar^4}{M^2v_m^2}
 + \frac{\Gamma(\ell-\frac12)\Lambda^{2\ell-1}}{\pi}.
\end{align}
In the vicinity of the confinement-induced resonance
$1/a_m\ll\Lambda^{2\ell+1}$, the two-body scattering $\T$ matrix
(\ref{eq:T-matrix}) has a pole in terms of the scattering energy
$\eps=\hbar^2k^2/(2M)$ at $\eps_m=-\hbar^2/(2Ma_mr_m)$, which is the
physical molecular energy and tunable in ultracold atom experiments.
These dimensionful quantities are customarily normalized with respect to
the Fermi wave number $\kF$ defined through the particle number density
of $A_\up$ atoms as $\kF^3/(6\pi^2)=\int\!d\k/(2\pi)^3f_T(\ek-\mu)$.
Consequently, the linear conductance (\ref{eq:conductance}) is
parameterized by $\eps_m/\eF$, $r_m/\kF^{2\ell-1}$, $\Lambda/\kF$, and
$T/\TF$ with $\eF=\hbar^2\kF^2/(2M)$ and $\TF=\eF/\kB$ being the Fermi
energy and temperature, respectively.

\begin{figure}[t]
 \includegraphics[width=0.92\columnwidth]{SU3-R1_arinv}\bigskip\\
 \includegraphics[width=0.92\columnwidth]{SU3-R1_logT}
 \caption{Linear conductance $G$ in the SU(3) symmetric case with
 $\ell=1$, $\Lambda=10\kF$, and $r_m=\Lambda$ as a function of the
 threefold degenerate molecular energy $\eps_m=-\hbar^2/(2Ma_mr_m)$ for
 selected $T/\TF=0.01$, 0.02, 0.05, 0.1, 0.2, 0.5, and 1 from highest to
 lowest curves (upper panel) and as a function of the temperature $T$
 for selected $\eps_m/\eF=0.6$, 0.65, 0.7, 0.75, 0.8, 0.85, and 0.9 from
 rightmost to leftmost curves (lower panel).  All curves in the lower
 panel are plotted in the same temperature range of
 $0.01\TF\leq T\leq\TF$ but normalized individually by $\TK$ for each
 $\eps_m/\eF$ so as to best fit to a common empirical form
 (\ref{eq:empirical}) of $hG(T)\approx(9/4)/[1+67.(T/\TK)^2]^{0.26}$
 indicated by the dashed curve.  \label{fig:SU3-R1}}
\end{figure}

While the formulation developed so far is general, we now focus on the
confinement-induced $p$-wave resonance in a dilute system with $\ell=1$
and $\Lambda=10\kF$.  Figure~\ref{fig:SU3-R1} shows the computed linear
conductance $G$ in units of the Planck constant $h=2\pi\hbar$ in the
SU(3) symmetric case with the resonance range $r_m=\Lambda$ chosen
corresponding to the case where the confined $B$ atom is lighter than
$A$ atoms~\cite{Nishida:2010}.  In the upper panel, $hG$ is plotted as a
function of the threefold degenerate molecular energy
$0\leq\eps_m/\eF\leq2$ for selected temperatures $T/\TF=0.01$, 0.02,
0.05, 0.1, 0.2, 0.5, and 1 accessible in ultracold atom experiments.
When the molecular energy is above the Fermi energy
($\eps_m\gtrsim\eF$), the conductance decreases toward the low
temperature until it eventually reaches a single Lorentzian curve with
its width set by the molecular lifetime.  On the other hand, when the
molecular energy is below the Fermi energy ($\eps_m\lesssim\eF$), the
conductance exhibits a strikingly distinct behavior, i.e., logarithmic
growth by lowering the temperature, which is the hallmark of the Kondo
effect.  Actually, the most remarkable feature of the Kondo effect is
the fact that the conductance turns out to be the universal function of
the temperature characterized by a single quantity $\TK$ in which all
microscopic details are encoded~\cite{Costi:1994}.  Such universality is
demonstrated in the lower panel of Fig.~\ref{fig:SU3-R1} where $hG$ for
selected molecular energies $\eps_m/\eF=0.6$, 0.65, 0.7, 0.75, 0.8,
0.85, and 0.9 are plotted in the same temperature range of
$0.01\TF\leq T\leq\TF$ but normalized individually by
$\TK\approx0.0072\TF$, $0.012\TF$, $0.019\TF$, $0.029\TF$, $0.044\TF$,
$0.063\TF$, and $0.081\TF$ so that all curves come together to make up a
single universal curve as much as possible.  The obtained universal
function is found to be well described by the empirical
form~\cite{Goldhaber-Gordon:1998} of
\begin{align}\label{eq:empirical}
 G(T) = \frac{G_0}{[1+(3^{1/\gamma}-1)(T/\TK)^2]^\gamma},
\end{align}
where $G_0=(2\ell+1)\sin^2[\pi/(2\ell+1)]/h$ is the zero temperature
conductance for the SU($2\ell+1$) Kondo effect in our
realization~\cite{Carmi:2011}, the Kondo temperature is defined by
$G(\TK)=G_0/3$, and the exponent $\gamma\approx0.26$ is chosen so as to
achieve the best fit to the numerical data.  While the Kondo temperature
decreases exponentially toward the molecular limit
$\eps_m/\eF\to-\infty$ as in Eq.~(\ref{eq:T_kondo}), we find a
reasonable range of molecular energy $0.5\lesssim\eps_m/\eF\lesssim1$
where the universal logarithmic growth of the conductance is observable
in ultracold atom experiments.

\begin{figure}[t]
 \includegraphics[width=0.92\columnwidth]{SU3-R5_arinv}\bigskip\\
 \includegraphics[width=0.92\columnwidth]{SU2-R1_arinv}
 \caption{Same as the upper panel of Fig.~\ref{fig:SU3-R1} with $\ell=1$
 and $\Lambda=10\kF$ but for a narrower resonance with $r_m=5\Lambda$
 (upper panel) and in a symmetry broken case to SU(2) with twofold
 degenerate $\eps_m=\eps_{\pm1}$ and decoupled $\eps_0\to\infty$ (lower
 panel).  \label{fig:SU3-R5_SU2-R1}}
\end{figure}

The linear conductance $G$ as a function of the molecular energy is also
computed for different sets of parameters as shown in
Fig.~\ref{fig:SU3-R5_SU2-R1} with $\ell=1$ and $\Lambda=10\kF$ retained.
The upper panel is in the SU(3) symmetric case as before but now for the
larger resonance range of $r_m=5\Lambda$ corresponding to the case where
the confined $B$ atom is heavier than $A$ atoms~\cite{Nishida:2010}.
While the Kondo effect is still observable in the plotted temperature
range, the Kondo temperature gets lower and the peak structure gets
narrower.  On the other hand, the lower panel is aimed at elucidating
the effect of symmetry breaking by setting the resonance range back to
$r_m=\Lambda$.  Here we suppose that the isotropic confinement potential
acting on the impurity $B$ atom is strongly deformed to a uniaxial one
so that the $p$-orbital degeneracy is reduced to twofold with
$\eps_m=\eps_{\pm1}$ and the nondegenerate molecular state is decoupled
with its energy $\eps_0\to\infty$.  The resulting peak structure of the
conductance in the lower panel of Fig.~\ref{fig:SU3-R5_SU2-R1} is found
to remain almost unchanged from that in the upper panel of
Fig.~\ref{fig:SU3-R1} and thus the observability of the Kondo effect is
not impaired by the symmetry breaking from SU(3) to SU(2).

\section{Concluding remarks}
In this Rapid Communication, we proposed and elaborated a simple and
versatile scheme to perform the conductance measurement with ultracold
atoms by employing spin superposition states, which can be implemented,
for example, in a Fermi gas of lithium atoms with impurity ytterbium
atoms~\cite{Hara:2014}.  In particular, we showed that a
confinement-induced $p$-wave or higher partial-wave resonance leads to
the universal logarithmic growth of the conductance toward the low
temperature, which is within reach of observation and thus provides a
clear evidence of the orbital Kondo effect in ultracold atom
experiments.  Not only the proposed transport measurement is applicable
both in and out of equilibrium, but our system is highly tunable and can
be easily extended to a dense Kondo lattice~\cite{Nishida:2013}, which
is difficult to realize in quantum dot systems in spite of its
importance to heavy fermion physics.  It will be interesting to further
incorporate a strong attraction between two components of fermions so
that they form Cooper pairs exhibiting the BCS-BEC crossover.  The
possibility to study such a rich variety of Kondo physics with ultracold
atoms is now opened up.

\acknowledgments
The author thanks Tomosuke Aono, Mikio Eto, Toshimasa Fujisawa, Ryotaro
Inoue, Norio Kawakami, Mikio Kozuma, and Yoshiro Takahashi for valuable
discussions.  This work was supported by JSPS KAKENHI Grant No.~15K17727
and MEXT KAKENHI Grant No.~15H05855.  Part of the numerical calculations
were carried out at the YITP computer facility in Kyoto University.