\section{Introduction}

The standard model (SM) of particle physics and the general relativity (GR) cannot explain everything in the universe, such as dark energy and that what occurs in the vicinity of a black hole. Thus, on the very high energy scales, one reconsiders combining the SM with GR in a unified theory, i.e., "quantum gravity". The standard model extension (SME), proposed by Kostelesk\'{y} and collaborators \cite{Kostelecky1,Kostelecky2,Kostelecky3,Kostelecky4,Kostelecky5,Kostelecky6,Kostelecky7,Kostelecky8}, is this effective field theory combining GR and SM at low energy scales, and incorporates all possible background fields that violate the fundamental symmetries existent in nature, the Lorentz invariance, happening on the high energy scales (about Planck scale). Studying the Lorentz violation (LV) is a useful approach toward investigating the foundations of modern physics. Besides SME, string theory \cite{kostelecky198939}, noncommutative field theories \cite{carroll,carroll2,carroll3}, spacetime-varying fields \cite{bertolami69,bertolami692,bertolami693}, loop quantum gravity theory \cite{gambini,gambini2}, brane world scenarios \cite{burgess03,burgess032,burgess033}, massive gravity \cite{fernando} and Einstein-aether theory \cite{jacobson,jacobson2} are other proposals of Lorentz violation.

SME can be used to calculate a number of predictions which can be tested in modern high-precision experiments \cite{Kostelecky8,tasson}. The primary LV in the gravity sector of SME is the form of $s^{\mu\nu}R_{\mu\nu}$, where $s^{\mu\nu}$ is a tensor field possessing a nonzero background configuration and defines preferred frames, $R_{\mu\nu}$ is the Ricci tensor. The simplest form of $s^{\mu\nu}$ is the bumblebee field $B^\mu$, which is expected to produce new and maybe even eccentric phenomena through the coupling term $B^{\mu}B^{\nu}R_{\mu\nu}$. The name of "bumblebee" comes from the alleged remark that bumblebees should not be allowed to fly since we did not know how their wings produce lift. The surprising property of this bumblebee gravity is that, unlike the absence of $U(1)$ gauge symmetry, it does not forbid the propagation of massless vector modes \cite{bluhm2008}. So one expects to reveal a variety of physical relics which may be of interest in studies of dark energy and dark matter due to the appearance of Nambu-Goldstone and massive Higgs in theories with the spontaneous Lorentz symmetry breaking \cite{Kostelecky6,Kostelecky7,bluhm}.

In this Einstein-bumblebee gravity theory, Bertolami and P\'aramos studied the 4-dimensional vacuum solutions including the purely radial, the radial/temporal or the axial/temporal LV \cite{bertolami}. In recent years, Casana {\it et al} \cite{casana} obtained an exact Schwarzschild-like black hole solution; Xu {\it et al} studied more general static black hole solutions and their properties \cite{xu}. In 2020, we obtained an exact Kerr-like solution and studied its black hole shadow \cite{ding2020}. Though this solution does not seem to satisfy the bumblebee field equations, Liu {\it et al} found that it can still be satisfied under some certain conditions, i.e., considered as an approximate solution of the bumblebee field equations \cite{liu}. After that we derived a slowly rotating black hole solution which satisfies all field equations \cite{ding2021}. Lately, we derived a black hole solution and a cosmological solution of Einstein-Gauss-Bonnet gravity coupled to bumblebee field \cite{ding2022}, found that the Guass-Bonnet term and the bumblebee field can both act as a form of dark energy. In  a high dimensional spacetime, we obtained an exact  AdS-like black hole solution \cite{ding2023}, found that the conceptions of black hole horizon area/entropy and volume inside horizon should be renewed due to its anisotropy.

In this paper, we would seek a rotating black hole solution in 3-dimensional spacetime. The first black hole solution is BTZ(Ba\~{n}ados, Teitelboim and Zanelli) solution \cite{banados1,banados2}, which is asymptotically anti-de Sitter(AdS) and has on curvature singularity. However, it is a genuine black hole solution due to the presence of event horizon, Hawking radiation, entropy,  and playing a significant role to understand physical properties in higher dimensions by using many of toy models \cite{bengtsson1,bengtsson2}.  The black hole thermodynamics has raised some challenging questions: a statistical derivation of black hole entropy and an account of its microstates. A promised idea is the holographic principle: there is a nontrivial match between features of 2-dimensional conformal field theory (CFT) and features of black holes \cite{castro}. The microscopic degrees of freedom of the black hole are described by CFT living in the boundary. Quantum studies around BTZ black holes can help better understand AdS/CFT correspondence, T-duality and U-duality to classes of asymptotically flat black strings \cite{horowitz1,horowitz2}.

We will derive a 3-dimensional  rotating black hole solution and study the central charges of the dual CFT in the theory of Einstein gravity coupled to the bumblebee fields. The rest paper is organized as follows.   In Sec. II we give the background for the Einstein-bumblebee theory. In Sec. III, we give the  black hole solution by solving the gravitational and bumblebee field equations. In Sec. IV, we study its central charges  and find some effects of the Lorentz breaking constant $\ell$. Sec. V is for a summary.

\section{Einstein-bumblebee gravity in 3-dimensional spacetime}

In the bumblebee gravity model, one introduces the bumblebee vector field $B_{\mu}$ which has a nonzero vacuum expectation value, to lead a spontaneous Lorentz symmetry breaking in the gravitational sector via a given potential. In the three dimensional spacetime, the
 action of Einstein-bumblebee gravity is \cite{ding2021},
\begin{eqnarray}
\mathcal{S}=
\int d^3x\sqrt{-g}\Big[\frac{R-2\Lambda}{2\kappa}+\frac{\varrho}{2\kappa} B^{\mu}B^{\nu}R_{\mu\nu}-\frac{1}{4}B^{\mu\nu}B_{\mu\nu}
-V(B_\mu B^{\mu}\mp b^2)+\mathcal{L}_M\Big], \label{action}
\end{eqnarray}
where $R$ is Ricci scalar and $\Lambda$ is the cosmological constant. $\kappa=4\pi G/c^4$ for the three dimensions, where $G$ is the Newtonian constant.

The coupling constant $\varrho$ dominates the non-minimal gravity interaction to bumblebee field $B_\mu$. The term $\mathcal{L}_M$ represents possible interactions with matter or external currents.
The constant $b$ is a real positive constant. The potential $V(B_\mu B^{\mu}\mp b^2)$ triggers Lorentz and/or $CPT$ (charge, parity and time) violation. It gives a nonzero vacuum expectation value (VEV) for bumblebee field $B_{\mu}$ indicating that the vacuum of this model obtains a prior direction in the spacetime. This potential has a minimum at $B^{\mu}B_{\mu}\pm b^2=0$ and $V'(b_{\mu}b^{\mu})=0$ to ensure the destroying the $U(1)$ symmetry, where the field $B_{\mu}$ acquires a nonzero VEV, $\langle B^{\mu}\rangle= b^{\mu}$. Another vector $b^{\mu}$ is a function of the spacetime coordinates and has a constant value $b_{\mu}b^{\mu}=\mp b^2$, where $\pm$ signs mean that $b^{\mu}$ is timelike or spacelike, respectively.
The bumblebee field strength is
\begin{eqnarray}
B_{\mu\nu}=\partial_{\mu}B_{\nu}-\partial_{\nu}B_{\mu}.
\end{eqnarray}
This antisymmetry of $B_{\mu\nu}$ implies the constraint \cite{bluhm}
\begin{eqnarray}
\nabla ^\mu\nabla^\nu B_{\mu\nu}=0.
\end{eqnarray}

Varying the action (\ref{action}) with respect to the metric yields the gravitational field equations
\begin{eqnarray}\label{einstein0}
G_{\mu\nu}+\Lambda g_{\mu\nu}=\kappa T_{\mu\nu}^B+\kappa T_{\mu\nu}^M,
\end{eqnarray}
where $G_{\mu\nu}=R_{\mu\nu}-g_{\mu\nu}R/2$, and the bumblebee energy momentum tensor $T_{\mu\nu}^B$ is
\begin{eqnarray}\label{momentum}
&&T_{\mu\nu}^B=B_{\mu\alpha}B^{\alpha}_{\;\nu}-\frac{1}{4}g_{\mu\nu} B^{\alpha\beta}B_{\alpha\beta}- g_{\mu\nu}V+
2B_{\mu}B_{\nu}V'\nonumber\\
&&+\frac{\varrho}{\kappa}\Big[\frac{1}{2}g_{\mu\nu}B^{\alpha}B^{\beta}R_{\alpha\beta}
-B_{\mu}B^{\alpha}R_{\alpha\nu}-B_{\nu}B^{\alpha}R_{\alpha\mu}\nonumber\\
&&+\frac{1}{2}\nabla_{\alpha}\nabla_{\mu}(B^{\alpha}B_{\nu})
+\frac{1}{2}\nabla_{\alpha}\nabla_{\nu}(B^{\alpha}B_{\mu})
-\frac{1}{2}\nabla^2(B^{\mu}B_{\nu})-\frac{1}{2}
g_{\mu\nu}\nabla_{\alpha}\nabla_{\beta}(B^{\alpha}B^{\beta})\Big].
\end{eqnarray}
The prime denotes differentiation with respect to the argument,
\begin{eqnarray}
V'=\frac{\partial V(x)}{\partial x}\Big|_{x=B^{\mu}B_{\mu}\pm b^2}.
\end{eqnarray}
Varying instead with respect to the the bumblebee field generates the bumblebee equations of motion (supposing that there is no coupling between the bumblebee field and $\mathcal{L}_M$),
\begin{eqnarray}\label{motion}
\nabla ^{\mu}B_{\mu\nu}=2V'B_\nu-\frac{\varrho}{\kappa}B^{\mu}R_{\mu\nu}.
\end{eqnarray}

The contracted Bianchi identities ($\nabla ^\mu G_{\mu\nu}=0$) lead to conservation of the total energy-momentum tensor
\begin{eqnarray}\label{}
\nabla ^\mu T_{\mu\nu}=\nabla ^\mu\big( T^B_{\mu\nu}+T^M_{\mu\nu}\big)=0.
\end{eqnarray}

We suppose that there is no matter field and the bumblebee field is frosted at its VEV like in Refs \cite{casana,bertolami}, i.e., it is
\begin{eqnarray}\label{bbu}
B_\mu=b_\mu.
\end{eqnarray}
And the potential has a linear function form \begin{eqnarray}
V=\frac{\lambda}{2}(B_\mu B^\mu-b^2)^2,
\end{eqnarray}
where $\lambda$ is a constant. So this potential is $V=0$ and  its derivative $V'=0$ under the condition (\ref{bbu}).  Then the first two terms in Eq. (\ref{momentum}) are like those of the electromagnetic field, the only distinctive are the coupling items to Ricci tensor. Under this condition,  Eq. (\ref{einstein0}) leads to gravitational field equations \cite{ding2021}
\begin{eqnarray}\label{bar}
G_{\mu\nu}=\kappa (\lambda b_\mu b_\nu+b_{\mu\alpha}b^{\alpha}_{\;\nu}-\frac{1}{4}g_{\mu\nu} b^{\alpha\beta}b_{\alpha\beta})+\varrho\Big(\frac{1}{2}
g_{\mu\nu}b^{\alpha}b^{\beta}R_{\alpha\beta}- b_{\mu}b^{\alpha}R_{\alpha\nu}
-b_{\nu}b^{\alpha}R_{\alpha\mu}\Big)
+\bar B_{\mu\nu},
\end{eqnarray}
with
\begin{eqnarray}\label{barb}
&&\bar B_{\mu\nu}=\frac{\varrho}{2}\Big[
\nabla_{\alpha}\nabla_{\mu}(b^{\alpha}b_{\nu})
+\nabla_{\alpha}\nabla_{\nu}(b^{\alpha}b_{\mu})
-\nabla^2(b_{\mu}b_{\nu})-g_{\mu\nu}\nabla_\alpha\nabla_\beta(b^\alpha b^\beta)\Big].
\end{eqnarray}
\section{Rotating BTZ-like black hole solution}
The static spherically symmetric black hole metric in the four dimensional spacetime  has the form
\begin{eqnarray}\label{metric}
&&ds^2=-e^{2\Phi(r)}dt^2+e^{2\psi(r)}dr^2+r^2\big[\Omega(r) dt+d\phi\big]^2,
\end{eqnarray}
where

In the present study, we pay attention to that the bumblebee field has a radial vacuum energy expectation because that the spacetime curvature has a strong radial variation, on the contrary that the temporal changes are very slow. So the bumblebee field is supposed to be spacelike($b_\mu b^\mu=$ positive constant) as that
\begin{eqnarray}\label{bu}
b_\mu=\big(0,b_0e^{\psi(r)},0\big),
\end{eqnarray}
where $b_0$ is a positive constant.
Then the bumblebee field strength is
\begin{eqnarray}
b_{\mu\nu}=\partial_{\mu}b_{\nu}-\partial_{\nu}b_{\mu},
\end{eqnarray}
whose components are all zero. And their divergences are all zero, i.e.,
\begin{eqnarray}
\nabla^{\mu}b_{\mu\nu}=0.
\end{eqnarray}
From the equation of motion (\ref{motion}), we have
\begin{eqnarray}
b^{\mu}R_{\mu\nu}=0\label{motion2}.
\end{eqnarray}
The gravitational field equations (\ref{bar}) become
\begin{eqnarray}\label{}
G_{\mu\nu}+\Lambda g_{\mu\nu}=\bar B_{\mu\nu}.
\end{eqnarray}

 For the metric (\ref{metric}), the nonzero components of Einstein tensor $G_{\mu\nu}$ and the bumblebee tensor $\bar B_{\mu\nu}$ are shown in the appendix.
Substituting the equation $G_{22}+\Lambda g_{22}=\bar B_{22}$ into $G_{02}+\Lambda g_{02}=\bar B_{02}$, one can obtain
\begin{eqnarray}\label{2202}
(1+\ell)\big[(r\Omega''+3\Omega')-r\Omega'(\Phi'+\psi')\big]=0,
\end{eqnarray}
where $\ell=\varrho b_0^2$ and the symbol $\prime$ means the derivative with their argument.
By using the motion equation (\ref{motion2}), one can obtain
\begin{eqnarray}\label{r11}
R_{11}=\frac{\psi'}{r}-(\Phi''+\Phi'^2-\Phi'\psi')+\frac{1}{2}e^{-2\Phi}r^2\Omega'^2=0.\end{eqnarray}
Substituting the equations $G_{22}+\Lambda g_{22}=\bar B_{22}$ and (\ref{r11}) into $G_{00}+\Lambda g_{00}=\bar B_{00}$, one can obtain
\begin{eqnarray}\label{2200}
(1+\ell)\left[r\Omega(r\Omega''+3\Omega')+\frac{1}{4}r^2\Omega'^2
-e^{2\Phi}\frac{\psi'}{r}-r^2\Omega\Omega'(\Phi'+\psi')\right]+e^{2\Phi+2\psi}\Lambda=0.
\end{eqnarray}
Substituting the equation (\ref{r11}) into $G_{11}+\Lambda g_{11}=\bar B_{11}$, one can obtain
\begin{eqnarray}\label{r11rr}
(1+\ell)\left[e^{2\Phi}\frac{\Phi'}{r}+\frac{1}{4}r^2\Omega'^2\right]+e^{2\Phi+2\psi}\Lambda=0.
\end{eqnarray}
The above two equations (\ref{2200}) and (\ref{r11rr}) can give that
\begin{eqnarray}\label{1100}
r\Omega\big[(r\Omega''+3\Omega')-r\Omega'(\Phi'+\psi')\big]-\frac{1}{r}e^{2\Phi}(\Phi'+\psi')=0.
\end{eqnarray}
Then the Eqs. (\ref{2202}) and (\ref{1100}) can give the following equations
\begin{eqnarray}
&&\Phi'+\psi'=0,\label{pp}\\&&r\Omega''+3\Omega'=0\label{omega}.
\end{eqnarray}
Eq. (\ref{omega}) can give
\begin{eqnarray}
\Omega=-\frac{j}{2r^2},
\end{eqnarray}
where $j$ is an integral constant relating to its angular momentum.
From the Eq. (\ref{pp}), one can assume that $e^{2\Phi}=f(r)$ and  $e^{2\psi}=C/f(r)$, where $C$ is a constant to be determined. Substituting it into the equations $G_{11}+\Lambda g_{11}=\bar B_{11}$ and $G_{22}+\Lambda g_{22}=\bar B_{22}$, one can obtain that
\begin{eqnarray}
(1+\ell)^2\big(2f'+\frac{j^2}{r^3}\big)+4C\Lambda r=0,
\end{eqnarray}
which can give that
\begin{eqnarray}
f(r)= -m-\frac{C\Lambda r^2}{(1+\ell)^2}+\frac{j^2}{4r^2},
\end{eqnarray}
where $m$ is also an integral constant relating to its mass.
In order to get a BTZ-like solution, we choose the constant $C=(1+\ell)^2$, then the last rotating BTZ-like black hole solution is
\begin{eqnarray}\label{metricb}
ds^2=-f(r)dt^2+\frac{(1+\ell)^2}{f(r)}dr^2+r^2(d\phi-\frac{j}{2r^2}dt)^2,\qquad f(r)= -m-\Lambda r^2+\frac{j^2}{4r^2}.
\end{eqnarray}
This metric represents a purely radial LV black hole solution in a 3D spacetime. Its Kretschmann scalar is
\begin{eqnarray}
R_{\mu\nu\rho\tau}R^{\mu\nu\rho\tau}=\frac{12\Lambda^2}{(1+\ell)^4},
\end{eqnarray}
which is a finite constant in the whole spacetime as like as the original BTZ black hole, so there is also no curvature singularity at the origin. Its total energy-momentum tensor of the metric (\ref{metricb}) is
\begin{equation}\label{totalmom}
\!\!\!\!\!\!\!\! T_\nu^\mu=\frac{1}{\kappa}\left(
\begin{array}{ccc}
-\epsilon  & 0  & 0 \\
0  & p_r &0    \\
0  & 0  & p_t
\end{array}
\right),
\end{equation}
where $\epsilon$ is the energy density, $p_r$ is the radial pressure and $p_t$ is the tangential pressure, which read
\begin{eqnarray}
-\epsilon=p_r=p_t=\frac{-\Lambda}{(1+\ell)^2}.
\end{eqnarray}

It has two horizons: inner(Cauchy) horizon $r_+$ and outer(event) horizon $r_-$ which can be read from $f(r_{\pm})=0$,
\begin{eqnarray}
r_{\pm}=\sqrt{\frac{2}{-\Lambda}}\left(\sqrt{M+\sqrt{-\Lambda}J}\pm\sqrt{M-\sqrt{-\Lambda}J}\right),
\end{eqnarray}
where $M=m/8$ is its ADM mass and $J=j/8$ is its angular momentum. Like Kerr black hole, it has an ergosphere $r_{erg}$ which can be read from the time-time component $g_{tt}=0$,
\begin{eqnarray}
r_{erg}=\sqrt{r_+^2-r_-^2}.
\end{eqnarray}
Its surface gravity of the event horizon $\kappa$ is
\begin{eqnarray}
\kappa=-\frac{1}{1+\ell}\Big(\Lambda r_++\frac{16J^2}{r_+^3}\Big),
\end{eqnarray}
which can be read from the formula \cite{ding2020},
\begin{eqnarray}
\kappa=-\frac{1}{2}\lim_{r\rightarrow r_+}\sqrt{\frac{1}{g_{rr}X}}\frac{dX}{dr}, \qquad X=g_{tt}-\frac{g_{t\phi}^2}{g_{\phi\phi}}.
\end{eqnarray}
It is easy to prove that in the present case, the function $X=-f(r)$.

\section{Central charge of the dual CFT}

According to the holographic theory \cite{witten}, one may expect that each sector of 3-dimensional  gravity which is either asymptotically anti-de Sitter(AdS) or AdS-like, there exists a dual 2-dimensional conformal field theory(CFT) that might live on the boundary of this AdS spacetime. From the black hole thermodynamics aspect, Yekta \cite{yekta} obtained the central charges of the CFT which can be constructed using the thermodynamics of the outer and inner horizons. He found the result was in complete agreement with that via the method of asymptotic symmetry group analysis \cite{liu1,liu2,liu3}. In this section, we study the central charges by using the thermodynamic method of black hole/CFT correspondence has been proposed in \cite{chen1,chen2,chen3}.

From the horizon equations $f(r_\pm)=0$, one can rewrite the black mass as
\begin{eqnarray}\label{}
M=-\frac{1}{8}\Lambda r_\pm^2+\frac{2J^2}{r_\pm^2}.
\end{eqnarray}
The temperature of outer horizon $T_+$ is defined by
\begin{eqnarray}\label{tp}
T_+=-\frac{1}{2\pi(1+\ell)}\Big(\Lambda r_++\frac{16J^2}{r_+^3}\Big)=\frac{-\Lambda}{2\pi(1+\ell)r_+}(r_+^2-r_-^2).
\end{eqnarray}
The temperature of the inner horizon $T_-$ should be a geometrical positive quantity and is constant over the inner horizon \cite{castro}, so it should be,
\begin{eqnarray}\label{tm} T_-=\frac{1}{2\pi(1+\ell)}\Big(\Lambda r_-+\frac{16J^2}{r_-^3}\Big)=\frac{-\Lambda}{2\pi(1+\ell)r_-}(r_+^2-r_-^2).
\end{eqnarray}
There are three thermodynamical quantities: entropy $S_\pm$ of the inner and outer horizons, angular momentum $J$ and pressure $P=-\Lambda/8\pi(1+\ell)^2$, i.e., mass $M$ can be expressed by the function of these quantities $M=M(S,J,P)$. The entropies are as following
\begin{eqnarray}\label{spm}
S_\pm=\int\frac{dM}{T_\pm}=\frac{1}{2}(1+\ell)\pi r_\pm,
\end{eqnarray}
the angular velocities $\Omega_\pm$ and the thermodynamical volumes $V_\pm$ are
\begin{eqnarray}\label{}
\Omega_\pm=\left(\frac{\partial M}{\partial J}\right)_{S,P}=\sqrt{8\pi P}\frac{r_\mp}{r_\pm},\qquad V_\pm=\left(\frac{\partial M}{\partial P}\right)_{S,J}=(1+\ell)^2\pi r_\pm^2.
\end{eqnarray}
So there exists the first law of thermodynamics and the Smarr formula for the outer horizon
\begin{eqnarray}\label{}
dM= T_+ dS_++\Omega_+ dJ+V_+ dP,\qquad 0\cdot M= T_+ S_+-2PV_++\Omega_+ J.
\end{eqnarray}
For the inner horizon, the Killing vector is spacelike inside the event horizon, then one should assign negative energy $-M$ to the inner horizon \cite{castro}, similar to the negative energies within the ergosphere. So  the first law of thermodynamics and the Smarr formula for the inner  horizon are
\begin{eqnarray}\label{}
dM=- T_- dS_-+\Omega_- dJ+V_- dP,\qquad 0\cdot M=- T_- S_--2PV_-+\Omega_- J.
\end{eqnarray}

It is obviously that the equality $T_+S_+=T_-S_-$ is true from the relations (\ref{tp}), (\ref{tm}) and (\ref{spm}). This means that the entropy product of the inner and outer horizon
\begin{eqnarray}\label{}
S_+S_-=(1+\ell)^2\frac{\pi^2J}{\sqrt{-\Lambda}},
\end{eqnarray}
is universal (mass-independent). So the central charges of the left- and right-moving sectors in the dual CFT must be the same,
 \begin{eqnarray}\label{}
c_L=c_R=6\frac{d}{dJ}\left(\frac{S_+S_-}{4\pi^2}\right)=\frac{3(1+\ell)^2}{2\sqrt{-\Lambda}}.
\end{eqnarray} These results of 2-dimensional CFT on the boundary can be used to describe 3-dimensional AdS gravity in the bulk.

Defining the left-moving temperature $T_L$ and right-moving temperature $T_R$ of the dual CFT,
 \begin{eqnarray}\label{}
T_L=\left(\frac{1}{T_+}+\frac{1}{T_-}\right)^{-1}=\frac{-\Lambda}{2\pi(1+\ell)}(r_+-r_-),\qquad T_R=\left(\frac{1}{T_+}-\frac{1}{T_-}\right)^{-1}=\frac{-\Lambda}{2\pi(1+\ell)}(r_++r_-),
\end{eqnarray}
then the entropies of the inner and outer horizon can be represented by
 \begin{eqnarray}\label{}
S_\pm=\frac{\pi^2}{3\sqrt{-\Lambda}}(c_LT_L\pm c_RT_R).
\end{eqnarray}
In addition, defining the energies of the left- and right-moving sectors of the dual CFT as
 \begin{eqnarray}\label{}
E_L=\frac{\pi^2}{6\sqrt{-\Lambda}}c_LT_L^2=\frac{M-\sqrt{-\Lambda}J}{2},\qquad E_R=\frac{\pi^2}{6\sqrt{-\Lambda}}c_RT_R^2=\frac{M+\sqrt{-\Lambda}J}{2},
\end{eqnarray}
then the black hole mass $M$ and angular momentum $J$ can be represented by
 \begin{eqnarray}\label{}
M=E_L+E_R, \qquad J=\frac{1}{\sqrt{-\Lambda}}(E_L-E_R).
\end{eqnarray}
In this way, we use the central charges on the boundary to reappear  the entropy, mass and angular momentum in the bulk.

\section{Summary}

In this paper, we have studied Einstein gravity coupled to a bumblebee field in a 3-dimensional spacetime. We obtain an exactly rotating BTZ-like black hole solution. The bumblebee field doesn't affect the locations of the black hole horizon and ergosphere.

This black hole is different from 4-dimensional Schwarzschild-like \cite{casana} or Kerr-like \cite{ding2020,ding2021} black holes that it cannot be asymptotically flat and has no curvature singularity at the origin. It is asymptotically AdS and has isotropic pressure, but is different from the 4- or higher dimensional AdS-like black holes \cite{maluf,ding2023} whose pressure are anisotropic.

We study the AdS/CFT correspondence of this black hole, find
that the entropy product of its inner and outer horizons is universal. So the central charges of the
dual CFT on the boundary can be obtained via the thermodynamic method, and they can reappear
black hole mass and angular momentum in the bulk.

\begin{acknowledgments} This work was supported by the Scientific Research Fund of the Hunan Provincial Education Department under No. 19A257, the National Natural Science Foundation (NNSFC)
of China (grant No. 11247013), Hunan Provincial Natural Science Foundation of China (grant No. 2015JJ2085 and No. 2020JJ4284).
\end{acknowledgments}


\section{Introduction}

Person re-identification (Re-ID) which aims at retrieving a specific person over a distributed set of non-overlapping cameras has made unprecedented progress~\cite{eom2019learning,Zhai2020ad,zheng2019pyramidal,zheng2019joint}.
Various efforts have been made to solve posture variations, resolution, occlusion, and blur under visible conditions.
Nevertheless, most surveillance systems use near-infrared (NIR) images instead of visible (VIS) images under low-light conditions, as the near-infrared cameras provide an efficient way to obtain high-quality images under poor illumination.
The problem of visible infrared re-identification (VI-ReID) needs to be considered in practical applications.

\begin{figure}[t]
    \centering
    \includegraphics[width=1.0\columnwidth]{fig2.pdf}
    \caption{\textbf{Geometrical analysis of the feature space optimized by softmax loss and spectral-aware softmax loss (SA-Softmax).}
    %
    $v$ and $n$ indicate the visible and near-infrared samples.
    %
    $W_{1}$ refers to the prototype of class 1.
    %
    $\theta_1$ and $\theta_2$ refer to the cosine distance of intra-class and inter-class in the retrieval phase, respectively.
    %
    (a) Since the original softmax loss does not explicitly optimize modality discrepancy, the query feature has to deal with an extra modality bias to match gallery samples from another side (Detailed discussed in Sec. \ref{sec:limitation}).
    %
    (b) The SA-Softmax optimizes the current samples using a prototype of another spectrum with the same identity, which efficiently explored the entire embedding space in the visible-infrared Re-ID.}
    \label{fig:softmax}
\end{figure}

Compared to the conventional Re-ID task, the primary challenge in such an application lies in the significant modality discrepancy between the VIS and NIR images, which leads to large intra-class variations.
Two typical frameworks are proposed to address this challenge from image generation and representation learning, respectively.
Image generation approaches~\cite{wang2018learning,wang2019aligngan,wang2019learning} aim to explore stable transformations across the spectrum so that samples under different spectral can be directly compared in the target domain.
Although with the progress and adaptation of image synthesis methods such as Generative Adversarial Networks (GAN) \cite{goodfellow2014generative,zhu2017unpaired}, this body of works shows impressive visual quality, the synthesized images are still far from photo-realistic.
Moreover, lacking paired cross-spectral images also makes the training much more challenging.
Another line of methods is representation learning approaches~\cite{dai2018cross,ye2018visible,ye2018hierarchical,lu2020cross,wu2020rgb}, which aim to bridge the modality gap in the common latent subspace.
These approaches combine a well-designed loss function with a single-path or multi-path deep neural network.
Benefiting from end-to-end optimization, these approaches are highly effective and often achieve state-of-the-art performance.

Despite the substantial progress, we observe that the softmax loss training paradigm, which is designed for a single modality classification task, still be considered as a basic training strategy in the cross-modality retrieval task.
Even though the softmax loss has made significant progress in single-modality tasks, as shown in Fig.~\ref{fig:softmax} (a), for the cross-spectral condition, the softmax loss considers the VIS and NIR samples from the same identity as a single class in the optimization, ignoring the gap between the two modalities.
The distributions of the same identity under different spectrums show limited overlapping in the softmax-trained feature space.
Hence, even after obtaining a well-trained feature distribution, the query feature still has to face an extra modality bias when finding a gallery sample from another side in the retrieval.
It makes the query feature easily influenced by noise gallery samples and induces a performance penalty.
Several previous works can make sense in this situation~\cite{schroff2015facenet,wen2016discriminative}, but they either rely on the limited inner-batch relations or only focus on intra-class samples, which do not explore the embedding space in its entirety.

In this paper, we propose a spectral-aware softmax (SA-Softmax) loss to enforce a higher similarity between the cross-modality sample pairs.
In general, the weight of the final fully connected layer of a deep neural network (DNN) trained with the softmax loss could be considered as a prototype for each identity~\cite{deng2019arcface,deng2021variational}.
Therefore, the SA-Softmax emphasizes the modality-gap of each modality in every single identity by dividing the prototype of each identity into VIS and NIR prototypes first.
As shown in Fig.~\ref{fig:softmax} (b), to increase the overlapping between the embedding of VIS and NIR distributions, the SA-Softmax utilizes the prototype from another spectrum to optimize the current spectral samples.
Since the feature embedding and the prototype matrix $W$ have different optimization goals in SA-Softmax, we utilize asynchronous optimizing processing for the feature embedding and the prototype matrix $W$.
Aside from the metric learning perspective, if regarding the spectrum-aware prototype matrix $W$ in SA-Softmax as a discriminator and the CNN feature extractor as a generator, the SA-Softmax can be considered as an adversarial learning strategy.

While SA-Softmax provides good eventual distribution, it still shows poor optimization for the feature embedding in two perspectives.
Firstly, although adopting the prototype from another spectrum as the target, the samples in the optimization will also be highly influenced by the prototype from the current spectrum.
This influence may decrease the similarity of sample pairs in the optimization phase (Detailed discussion can refer to the Sec.~\ref{sec:FM}).
To overcome this drawback, we employ a feature mask to exclude the response between the sample and the prototype from the current spectrum for better optimization.
Secondly, we observe that Softmax loss is more concerned with obtaining a high relative results rather than a high similarity for each identity, implying that it has little encouragement to provide a high intra-class similarity (Detailed discussion can be found in Sec.~\ref{sec:AST}).
Unfortunately, the SA-Softmax inherits this weakness as well.
Hence, we add an absolute-similarity term in the SA-Softmax to decrease the variation caused by the modality, pose, and occlusion in the VI-ReID and make the distribution of every identity more compact.

The main contributions of the paper are summarized as follows:
\begin{itemize}
    \item We point out the weakness of the softmax loss training paradigm in the visible-infrared person re-identification problem and propose an asynchronous optimized Spectral-Aware Softmax (SA-Softmax) loss that fully explores the embedding space with the modality information.

    \item Combining with the feature mask and the absolute-similarity term, we present an effective framework based on the SA-Softmax for the VI-ReID task.

    \item Extensive experiments on two public datasets, RegDB and SYSU-MM01, demonstrate the superiority of our SA-Softmax loss.
\end{itemize}

\section{Related Works}
\textbf{Visible-infrared Person Re-ID.}
On account of the recent advances in deep learning, two typical frameworks are proposed to address this challenge from image generation and representation learning respectively.
Image generation approaches attempt to build an efficient transformation among the modalities. Wu et al. \cite{wu2017rgb} propose a zero-padding framework to align the cross-modality images in a common space.
D$^2$RL \cite{wang2019learning} regards each of the samples in a hyperspectral perspective and utilizes the variational autoencoders (VAE) to supply the missing channels for each modality.
AlignGan \cite{wang2019aligngan} employs a cycle-gan framework with efficient constraints to align from the pixel-level and feature-level.
Hi-CMD \cite{choi2020hi} uses a GAN framework to change the pose and illumination attributes for each sample to capture ID-discriminative and color-irrelevant representations.
Ye et al. \cite{ye2021channel}, X-modality \cite{li2020infrared}, and SMCL \cite{ye2018visible} generate an extra modality for jointly learning. On the other hand, representation learning approaches attempt to obtain a modality-irrelevant feature space through module design and constraints.
From the representation learning perspective, Dai et al. \cite{dai2018cross} employ an adversarial learning strategy to tackle the modality discrepancy.
DDAG \cite{ye2020dynamic} uses the graph structure to aggregate representations from local to global.
Wu et al. \cite{wu2021discover} attempt to discover modality-irrelevant nuances in different patterns to match the cross-modality pedestrian samples in a fine-grid way.
Hao et al. \cite{hao2021cross} combine adversarial learning with center aggregation to extract the centralization features with diversity.
Despite significant achievements in the VI-ReID task over the past few years, the softmax loss training paradigm designed for a single modality classification task is still considered as a basic training strategy in this area.
As has been aforementioned, the drawback of softmax loss in such a cross-modality condition is obvious.

\textbf{Deep Metric Learning.}
As the optimization target, the loss function plays an essential role in feature representation learning.
Based on the most widely used Softmax loss, Large Margin Softmax \cite{liu2016large} firstly adds an extra margin to learn a more discriminative feature embedding.
Since the softmax loss does not explicitly optimize the cosine-similarity, Normface~\cite{wang2017normface} applies the $L2$ normalization to the weight matrix and feature.
More work is now being done to investigate a better strategy for adjusting the adding margin in the L2-normalized softmax, such as SphereFace~\cite{liu2017sphereface}, CosFace~\cite{wang2018cosface}, and ArcFace~\cite{deng2019arcface}.
MagFace~\cite{meng2021magface} emphasizes the value of both direction and magnitude for the feature vector and shows the quality of faces through the magnitude of feature embedding.
While the softmax loss strategy mainly focuses on the cosine similarity, some other works consider the Euclidean feature space as well.
Triplet loss explores the optimization for pair-wise distance by the triplet relation in each mini-batch.
Center loss~\cite{wen2016discriminative} adopts the center of each class to obtain a more compact intra-class distribution.
Circle loss~\cite{sun2020circle} provides a unified perspective for maximizing intra-class similarity while minimizing inter-class similarity.
Although these metric learning methods perform well on related tasks such as face recognition or person re-identification, they are all designed for single modality tasks that follows the prior that all class samples come from the same domain.
The straight idea that clusters all samples from the same class together ignores the characteristics of multi-modality data.
Therefore, adopting these learning strategies under cross-modality conditions will face inevitable modality gaps and achieve limited performance.

\section{Method}
\subsection{The limitation of Softmax Loss}
\label{sec:limitation}
Generally, given extracted feature embedding $x_i$, the most widely used loss paradigm in classification, softmax loss, can be formulated as:
\begin{equation}
\begin{split}
\label{eq:orisoftmax}
&\mathcal{L}_{s}(W, x_i) = -log\frac{e^{s_{i}\cos(\theta_{i})}}{\sum^N_{j=1}e^{s_j\cos(\theta_{j})}},\\
& with \quad s_j=\left \|W^{T}_{j} \right \| \left \|x_{i} \right \|, \theta_{j}=\left \langle x_{i}, W_{j} \right \rangle,
\end{split}
\end{equation}
where $W_j \in \mathbb{R}^d$ denotes the $j$-th column of the prototype matrix $W \in \mathbb{R}^{d \times N}$. $d$ is the dimension of embedding feature, $N$ is the number of class, and $x_i \in \mathbb{R}^d$ refers to the $i$-th training sample, which belongs to the $y_i$-th class.

Eq. \ref{eq:orisoftmax} indicates that the original softmax loss aims to obtain an ideal embedding space by enlarging the intra-class similarity while decreasing the inter-class similarity.
However, we can also observe several typical drawbacks of the original softmax that limit its performance in cross-spectral conditions.
\textbf{Firstly}, the prototype matrix $W$ is optimized by the whole identity samples containing both two modalities.
In general, the optimization of softmax loss is considered as two synchronous parts: the prototype matrix $W$ and the embedding feature $x_i$, whose derivatives are described respectively as:
\begin{equation}
\begin{split}
\label{eq:smd}
&\frac{\partial \mathcal{L}_{softmax} }{\partial x_i} = \sum_{j=1}^{N}(p_{ij}-\mathit{1} (j==y_i))W_j,\\
&\frac{\partial \mathcal{L}_{softmax} }{\partial W_i} = \sum_{i=1}^{B}(p_{ij}-\mathit{1} (j==y_i))x_i,\\
&with \qquad p_{ij} = \frac{e^{W^{T}_{y_i}}x_i}{\sum^N_{j=1}e^{W^{T}_{j}x_i}},
\end{split}
\end{equation}
where $y_i$ is the label of $x_i$, and B refers to the batch size.

During the training phase, as indicated in Eq. \ref{eq:smd}, the softmax loss aims to optimize the sample-to-prototype similarity, while the prototype matrix $W$ optimized by the two modality samples is not a suitable target in such a condition.
Specifically, as shown in Fig. \ref{fig:softmax} (a), for the samples $v_1$ and $n_1$ from the same identity and $n_2$ from another identity, although they have achieved a good training result as $\cos (\left \langle v_{1}, W_{1} \right \rangle) > \cos (\left \langle v_{1}, W_{2} \right \rangle)$, $\cos (\left \langle n_{1}, W_{1} \right \rangle) > \cos (\left \langle n_{1}, W_{2} \right \rangle)$, and $\cos (\left \langle n_{2}, W_{2} \right \rangle) > \cos (\left \langle n_{2}, W_{1} \right \rangle)$, a error retrieval result still occurs as $\cos (\left \langle v_{1}, n_{2} \right \rangle) > \cos (\left \langle v_{1}, n_{1} \right \rangle)$.

\textbf{Secondly}, the softmax is not a balanced optimization for every identity.
The first-order derivative and partial derivative of $L(\theta_i)$ are:
\begin{equation}
\begin{split}
\label{eq:ts}
&\mathcal{L}_{s}'(\theta_i)=\frac{e^{s\cos(\theta_j)}s}{e^{s\cos(\theta_i)}+e^{s\cos(\theta_j)}}\sin(\theta_i)>0,\\
&\frac{\partial \mathcal{L}_{s}'(\theta_i)}{\partial \theta_j} =-\frac{e^{s\cos(\theta_i)+s\cos(\theta_j)}s^2\sin(\theta_j)}{(e^{s\cos(\theta_i)}+e^{s\cos(\theta_j)})^2}\sin(\theta_i)<0.\\
\end{split}
\end{equation}

From Eq.\ref{eq:ts}, we can observe that the $\mathcal{L}_{s}'(\theta_i)$ increases monotonically and strictly with the decrease of the $\theta_j$.
It indicates that those distributions with a high $\theta_j$ at the boundary of the training set are less encouraged to learn compact distributions.
Since the training set and the testing set come from different identities, the distributions of the two parts are not well aligned.
Therefore, those distributions at the boundary of the training set also play an important role in learning a more efficient model.

\textbf{Thirdly}, even ignoring the influence of $s$, the optimization is still inefficient. From Eq.\ref{eq:ts}, after $\theta_i < \theta_j$, the $\frac{e^{s\cos(\theta_j)}s}{e^{s\cos(\theta_i)}+e^{s\cos(\theta_j)}}$ in $\mathcal{L}_{s}'(\theta_i)$ will also drop rapidly and show less ability to further decrease the $\theta_i$.

Therefore, we propose a spectral-aware loss to overcome the above drawback of softmax loss in the cross-spectral task.

\begin{figure*}[t]
\centering
\includegraphics[height=8.4cm,width=18cm]{optimization.pdf}
\vspace{-1em}\caption{\textbf{The optimization process of SA-Softmax loss.}
(a) The softmax loss uses a synchronization strategy to optimize the network parameters and prototype matrix.
(b) For SA-Softmax loss, the network and prototype have different goals.
Hence, SA-Softmax adopts an asynchronous strategy to optimize the network parameters and prototype matrix.
(c) The two-step asynchronous framework for SA-Softmax loss.}
\label{fig:optim}
\end{figure*}

\subsection{Spectral-Aware Softmax Loss}
Despite its success in single modality tasks, using a prototype to aggregate two modality samples with the same identity lacks explicit optimization for the modality discrepancy and introduces a typical gap in the final feature space, as discussed above.
Therefore, SA-Softmax inherits the idea of softmax loss and attempts to fully utilize the extra modality information in each sample pair.

Specifically, as shown in Fig.~\ref{fig:optim}, the SA-Softmax firstly utilizes the modality prototype $[W^v, W^n] \in \mathbb{R}^{d \times 2N}$ instead of the identity prototype $W \in \mathbb{R}^{d \times N}$. Herein, $d$ is the dimension of prototypes, and $N$ is the number of classes.
Then, to enforce a high overlapping between the distribution of two modalities, the SA-Softmax uses the prototype from another spectrum to optimize each sample.
To achieve this kind of cross optimization for each sample and corresponding prototype, the SA-Softmax employs an asynchronous optimization strategy for the prototype matrix $[W^v, W^n]$ and sample feature $x_i$, as shown in Fig.~\ref{fig:optim}.
Therefore, the label $y_i$ of $x_i$ is rewritten to $y^W_i$ and $y^F_i$, respectively, to optimize the prototype matrix and feature extractor as:

\begin{equation}
\begin{split}
\label{eq:saslabel}
\begin{cases}
  & y^W_i = y_i, y^F_i = y_i + N \qquad \text{ if } x_i \in VIS \\
  & y^W_i = y_i + N, y^F_i = y_i \qquad \text{ if } x_i \in NIR.
\end{cases}
\end{split}
\end{equation}

For the prototype matrix $[W^v, W^n]$, this process can be formulated as:
\begin{equation}
\begin{split}
\label{eq:sasw}
&\frac{\partial \mathcal{L}^W_{sas} }{\partial [W^v, W^n]} = \sum_{i=1}^{B}(p_{ij}-\mathit{1} (j==y^W_i))x_i,\\
&with \qquad p_{ij} = \frac{e^{[W^v, W^n]^{T}_{y^W_i}x_i}}{\sum^{2N}_{j=1}e^{[W^v, W^n]^{T}_{j}x_i}}.
\end{split}
\end{equation}
And for feature embedding $x_i$, this process can be formulated as:
\begin{equation}
\begin{split}
\label{eq:sasf}
&\frac{\partial \mathcal{L}^F_{sas} }{\partial x_i} = \sum_{j=1}^{N}(p_{ij}-\mathit{1} (j==y^F_i))[W^v, W^n],\\
&with \qquad p_{ij} = \frac{e^{[W^v, W^n]^{T}_{y^F_i}x_i}}{\sum^{2N}_{j=1}e^{[W^v, W^n]^{T}_{j}x_i}}.
\end{split}
\end{equation}
Here, the $\mathcal{L}^W_{sas}$ and $\mathcal{L}^F_{sas}$ denote the softmax loss results under the label $y^W$ and $y^F$, respectively.
In summary, the SA-Softmax could be given by:
\begin{equation}
\begin{split}
\label{eq:sas}
&\mathcal{L}_{sas}(W, x_i) = \mathcal{L}^W_{sas}(W) + \mathcal{L}^F_{sas}(x_i),\\
& \mathcal{L}^W_{sas}(W) = -log\frac{e^{[W^v, W^n]^{T}_{y^W_i}x_i}}{\sum^{2N}_{j=1}e^{[W^v, W^n]^{T}_{j}x_i}},\\
& \mathcal{L}^F_{sas}(x_i) = -log\frac{e^{[W^v, W^n]^{T}_{y^F_i}x_i}}{\sum^{2N}_{j=1}e^{[W^v, W^n]^{T}_{j}x_i}}.
\end{split}
\end{equation}

\subsection{Feature Mask}
\label{sec:FM}
\begin{figure}[t]
    \centering
    \includegraphics[height=4.5cm,width=8.4cm]{fig3.pdf}
    \caption{\textbf{A geometrical interpretation of the feature mask.}
    %
    $d_1$ and $d_2$ refer to the pair-wise distance for same identity before and after the optimization,respectively.
    %
    (a) Ambiguous optimization in the original SA-Softmax.
    %
    The prototype of the current modality will push away the training samples, which may lead to ambiguous conditions where the distance between the target prototype and samples decreases (loss decreases) while the discrepancy between the intra-class sample pairs becomes larger.
    %
    (b) Feature mask ignores the influence from the prototype of the current modality, resulting in more stable training.}
    \label{fig:mask}
\end{figure}
Although the idea of using the prototype from another modality provides a more compact intra-class distribution for the cross-modality samples, the prototype of the current modality of each sample is not insignificant.
As shown in Fig.~\ref{fig:mask} (a), the prototype of the current modality will greatly affect the optimization and may decrease the similarity of pair-wise samples, which would defeat the purpose of SA-Softmax.
To overcome the tendency, we apply a Feature Mask to eliminate the influence of the current modality prototype for each identity.
In Fig.~\ref{fig:mask} (b), we show how the feature mask works.
After combining the feature mask, the $L^F_{sas}$ can be formulated as:
\begin{equation}
\label{eq:sasFM}
\mathcal{L}^F_{sas}(x_i) = -log\frac{e^{[W^v, W^n]^{T}_{y^F_i}x_i}}{\sum^{2N}_{j=1, j \ne y^W_i}e^{[W^v, W^n]^{T}_{j}x_i}}.
\end{equation}

Compared to the original SA-Softmax, adding the feature mask can significantly alleviate the ambiguous optimization, which indicates better learning results.

\subsection{Absolute-Similarity Term}
\label{sec:AST}
As mentioned above, in Eq.~\ref{eq:orisoftmax} and Eq.~\ref{eq:sas}, we can observe that the softmax loss uses a relative strategy to balance the intra/inter similarity.
It indicates that softmax loss is not a balance optimization.
Those outlier identities with lower $D_{inter}$ are less likely to learn a more compact intra-class distribution, making them more vulnerable to changes in pose, illumination, and so on.
On the other hand, owing to its design, even for those with normal identities, the softmax loss shows limited ability to push them to be more similar.
As shown in Eq. \ref{eq:orisoftmax}, the distance $D_{intra}$ can be formulated as:
\begin{equation}
\label{eq:dinter}
D_{intra}(x_i) = e^{W^{T}_{y_i}x_i} = e^{\left \|W_{y_i} \right \| \left \|x_i \right \| cos<W_{y_i}, x_i>},
\end{equation}
where $<W_{y^F_i}, x_i>$ is the angle between $W_{y_i}$ and $x_i$.
Although it aims to increase the cosine similarity between $W_{y^F_i}$ and $x_i$, from Eq.~\ref{eq:dinter}, the $D_{intra}$ also enjoys the increase of $\left \|W_{y_i} \right \|$ and $\left \|x_i \right \|$, resulting in limited results.
As mentioned in several previous works \cite{wu2020rgb}, since the identities during training and testing do not overlap, the feature has different distributions in the training and testing sets.
Therefore, if we consider each identity as a single small domain, it is obvious that an ideal compact distribution with less variance is an equal demand for every identity.

To alleviate this problem, we add an extra Absolute-Similarity Term (AST) to directly enforce that each sample is more similar to its target prototype.
Following the formulation of the softmax loss, the AST can be given by:
\begin{equation}
\label{eq:ast}
\mathcal{L}_{AST}(x)= \sum_{i=1}^{B}\left \|1-cos<W_{y^F_i}, x_i>  \right \|_2.
\end{equation}

In summary, the final loss function could be shown as:
\begin{equation}
\label{eq:final}
\mathcal{L}
=
\mathop{\underline{\alpha \mathcal{L}_{sas} + (1-\alpha) \mathcal{L}_{softmax}}}_{Relative\; Term}
+
\mathop{\underline{\beta \mathcal{L}_{AST}}}_{Absolute\; Term},
\end{equation}
where $\alpha$ and $\beta$ are the hyper-parameters to trade-off among each part.

\section{Experiment}
\subsection{Experimental Setting}
\textbf{Datasets.} We conduct experiments on two publicly available visible-infrared person re-identification datasets SYSU-MM01 \cite{wu2017rgb} and RegDB \cite{nguyen2017person}.
\begin{itemize}
    \item \textbf{SYSU-MM01} is a large-scale dataset captured by four visible cameras and two infrared cameras in both indoor and outdoor environments.
    %
    The training set contains 395 identities with $22,258$ visible images and $11,909$ infrared images, while the testing set includes $96$ identities with $3,803$ infrared images as the query.
    %
    This dataset contains two different search modes, the all-search mode and the indoor-search mode. In the all-search mode, the gallery images are from all the visible cameras. For the indoor-search mode, the source of the gallery set excludes two outdoor cameras.
    \item \textbf{RegDB} dataset is collected by two aligned cameras, one for visible and the other for far-infrared (thermal).
    %
    It contains $412$ identities, each with $10$ visible images and $10$ infrared images. Following the evaluation protocol of previous works \cite{ye2018hierarchical,ye2020cross}, we choose half of the identities at random for training and the other half for testing.
    %
    The results are the average of $10$ repeating.
\end{itemize}

\textbf{Evaluation Protocol.} We follow the evaluation settings in existing VI-ReID methods\cite{wu2021discover,ye2021channel} and adopt the widely used Cumulative Matching Characteristic (CMC) and mean Average Precision (\emph{m}AP) as evaluation metrics.

\textbf{Implementation details.} We use Pytorch to implement our method and finish all the experiments on a single RTX 3090 GPU.
The mini-batch size is set to 128.
For each mini-batch, we randomly select 8 identities, each with 8 visible images and 8 infrared images.
We select the ResNet-50 based PCB \cite{sun2018beyond} model as the baseline, which is a widely used fine-grid part feature learning framework in both Re-ID and visible-infrared Re-ID.
Following previous works \cite{ye2020dynamic,hao2021cross}, we divide the first convolutional layer to tackle the two modalities input.
We resize all of the images to 384 $\times$ 192 and use random flipping and random erasing \cite{zhong2020random,ye2021channel} to augment the data.
To train the model, the SGD optimizer is used with an initial learning rate of 0.01, which is divided by 10 at the 40th and 80th epoch.

\begin{table*}[t]
  \centering
  \renewcommand{\arraystretch}{1.2}
  \caption{\textbf{Comparison with the state-of-the-arts on SYSU-MM01 and RegDB datasets}. R-1, 10, 20 denotes the Rank-1, 10, 20 accuracy.}
  \resizebox{\textwidth}{!}{
  \begin{tabular}{r|c|cccc|cccc|cccc|cccc}
  \hline
   \multirow{3}{*}{Model} & \multirow{3}{*}{Venue} & \multicolumn{8}{c}{RegDB} & \multicolumn{8}{c}{SYSU-MM01} \\
   \cline{3-18}
   & &\multicolumn{4}{c}{Visible to Thermal}&\multicolumn{4}{c}{Thermal to Visible} &\multicolumn{4}{c}{All Search} &\multicolumn{4}{c}{Indoor Search} \\
   \cline{3-18}
   & & R-1 & R-10 & R-20 & \emph{m}AP   & R-1 & R-10 & R-20 & \emph{m}AP & R-1 & R-10 & R-20 & \emph{m}AP & R-1 & R-10 & R-20 & \emph{m}AP\\
   \hline
    Zero-Padding\cite{wu2017rgb}        & ICCV'17   & 17.8 & 34.2 & 44.4 & 18.9 & 16.6 & 34.7 & 44.3 & 17.8 & 14.8 & 54.1 & 71.3 & 15.9 & 20.6 & 68.4 & 85.8 & 26.9 \\
    HCML\cite{ye2018hierarchical}       & AAAI'18   & 24.4 & 47.5 & 56.8 & 20.8 & 21.7 & 45.0 & 55.6 & 22.2 & 14.3 & 53.2 & 69.2 & 16.2 & 24.5 & 73.3 & 86.7 & 30.1 \\
    BDTR\cite{ye2018visible}            & IJCAI'18  & 33.6 & 58.6 & 67.4 & 32.8 & 32.9 & 58.5 & 68.4 & 32.0 & 17.0 & 55.4 & 72.0 & 19.7 & - & - & - & - \\
    cmGAN\cite{dai2018cross}            & IJCAI'18  & -    & -    & -    & -    & -    & -    & -    & -    & 27.0 & 67.5 & 80.6 & 27.8 & 31.7 & 77.2 & 89.2 & 42.2 \\
    AlignGAN\cite{wang2019aligngan}     & ICCV'19   & 57.9 & -    & -    & 53.6 & 56.3 & - & - & 53.4 & 42.4 & 85.0 & 93.7 & 40.7 & 45.9 & 87.6 & 94.4 & 54.3 \\
    X-Modality\cite{li2020infrared}     & AAAI'20   & 62.2 & 83.1 & 91.7 & 60.2 & - & - & - & - & 49.9 & 89.8 & 96.0 & 50.7 & -    & -    & -    & -    \\
    MGE + FMASP\cite{wu2020rgb}         & IJCV'20   & 65.1 & 83.7 & - & 64.5 & - & - & - & - & 43.6 & 86.3 & - & 45.0 & 48.6 & 89.5 & - & 57.5 \\
    SSFT\cite{lu2020cross}             & CVPR'20   & 72.3 & -    & - & 72.9 & 71.0 & - & - & 71.7 & 61.6 & 89.2  & 93.9 & 63.3  & 70.5 & 94.9 & 97.7 & 72.6 \\
    DDAG\cite{ye2020dynamic}            & ECCV'20   & 69.3 & 86.2 & 91.5 & 63.5 & 68.1 & 85.2 & 90.3 & 61.8 & 54.8 & 90.4 & 95.8 & 53.0 & 61.0 & 94.1 & 98.4 & 68.0 \\
    DG-VAE\cite{pu2020dual}             & MM'20     & 73.0 & 86.9 & - & 71.8 & - & - & - & - & 59.5 & 93.8 & - & 58.5 & - & - & - & - \\
    CICL + IAMA\cite{zhao2021joint}     & AAAI'21   & 78.8 & -    & - & 69.4 & 77.9 & - & - & 69.4 & 57.2 & 94.3 & 98.4 & 59.3 & 66.6 & 98.8 & 99.7 & 74.7 \\
    VCD + VML\cite{tian2021farewell}    & CVPR'21   & 73.2 & -    & - & 71.6 & 71.8 & - & - & 70.1 & 60.0 & 94.2 & 98.1 & 58.8 & 66.1 & 96.6 & 99.4 & 73.0 \\
    MPANet\cite{wu2021discover}         & CVPR'21   & 82.8 & -    & - & 80.7 & 83.7 & - & - & 80.9 & 70.6 & 96.2 & 98.8 & 68.2 & 76.7 & 98.2 & 99.6 & 81.0 \\
    MCLNet\cite{hao2021cross}           & ICCV'21   & 80.3 & 92.7 & 96.0 & 73.1 & 75.9 & 90.9 & 94.6 & 69.5 & 65.4 & 93.3 & 97.1 & 62.0 & 72.6 & 97.0 & 99.2 & 76.6 \\
    SMCL\cite{wei2021syncretic}         & ICCV'21   & 83.9 & -    & - & 79.8 & 83.1 & - & - & 78.6 & 67.4 & 92.9 & 96.8 & 61.8 & 68.8 & 96.6 & 98.8 & 75.6 \\
    CAJ\cite{ye2021channel}             & ICCV'21   & 85.0 & 95.5 & 97.5 & 79.1 & 84.8 & 95.3 & 97.5 & 77.8 & 69.9 & 95.7 & 98.5 & 66.9 & 76.3 & 97.9 & 99.5 & 80.4 \\
    \hline
    SA-Softmax                          & \textbf{-} & \textbf{95.2} & \textbf{99.0} & \textbf{99.6} & \textbf{88.7} & \textbf{93.1} & \textbf{98.1} & \textbf{99.1} & \textbf{86.9} & \textbf{71.0} & \textbf{96.9} & \textbf{99.3} & \textbf{68.5} & \textbf{78.8} & \textbf{99.1} & \textbf{99.8} & \textbf{82.2} \\
    \toprule[1pt]
    \end{tabular}}
    \label{Tbl:final}
\end{table*}

\begin{figure*}[t]
\centering
\includegraphics[height=4.5cm,width=18cm]{fig5.pdf}
\vspace{-1em}\caption{\textbf{Impact of the hyper-parameters in terms of CMC (\%) and \emph{m}AP (\%) on SYSU-MM01.}
For both all-search and indoor-search modes, the performance peaks when  $\alpha$  and $\beta$ are set to 0.7 and 1.0, respectively.}
\label{fig:para}
\end{figure*}

\subsection{Comparison with State-of-the-art Methods}
To comprehensively demonstrate the performance of SA-Softmax, we make a comparison between the combination of the PCB model and SA-Softmax with recent state-of-the-art visible-infrared Re-ID methods in both SYSU-MM01 and RegDB.
The results are shown in Tab.~\ref{Tbl:final}.
It is clear that the PCB+SA-Softmax outperforms the existing SOTAs.
Especially in the RegDB dataset, our method achieves an impressive performance, with Rank-1 accuracy of 95.2\% and \emph{m}AP of 88.7\%  for visible to thermal, and Rank-1 accuracy of 93.1\% and \emph{m}AP of 86.9\% for thermal to visible, respectively.
Meanwhile, in SYSU-MM01, our method also outperforms the SOTAs in both all-search and indoor-search modes.
Since the SYSU-MM01 datasets include more divergence in the pose, viewpoints, and background than the RegDB dataset, which focuses primarily on the modality discrepancy, we believe this is why the SA-Softmax achieves fewer improvements in SYSU-MM01 than in RegDB.
Note that SA-Softmax is a plug-and-play module, our network is trained based on the PCB model, and the performance may be boosted with better data augmentations as well as extra well-designed network modules.

\begin{table}[t]
  \centering
  \renewcommand{\arraystretch}{1.2}
  \caption{\textbf{Ablation study of each components in SA-Softmax in terms of CMC (\%) and \emph{m}AP (\%) on SYSU-MM01.} 'SAS' denotes the original SA-Softamx loss. 'FM' denotes the Feature Mask. 'AST' refers to the Absolute-Similarity Term.}
  \resizebox{85mm}{!}{
  \begin{tabular}{ccc|cc|cc}
  \toprule[1pt]
   \multicolumn{3}{c|}{\textbf{Module}} & \multicolumn{4}{c}{\textbf{SYSU-MM01}} \\
   \cline{1-7}
   \multicolumn{3}{c|}{\emph{Setting}} &\multicolumn{2}{c|}{\emph{All Search}} &\multicolumn{2}{c}{\emph{Indoor Search}}\\
   \cline{1-7}
   SAS & FM & AST & R-1 & \emph{m}AP   & R-1 & \emph{m}AP\\
  \hline
               &           &      & 64.5 & 61.3 & 71.6 & 76.1 \\
    \checkmark &           &      & 69.6 & 66.3 & 77.1 & 81.0 \\
    \checkmark &\checkmark &      & 70.3 & 66.8 & 78.3 & 81.7 \\
    \checkmark &\checkmark &\checkmark      & \textbf{71.0} & \textbf{68.5} & \textbf{78.8} & \textbf{82.2}\\
  \bottomrule[1pt]
    \end{tabular}}
    \label{Tbl:ablation}
\end{table}

\subsection{Ablation Study}
To study the effectiveness of the proposed SA-Softmax loss, we conduct ablation experiments on the baseline PCB model.
In Tab.~\ref{Tbl:ablation}, we show the results of SA-Softmax with different components on the SYSU-MM01.
Compared to the baseline, the original SA-Softmax greatly improves the performance of SYSU-MM01 in both all-search and indoor-search settings.
Then, the performance could improve further after adding the feature mask.
Finally, combined with the absolute-similarity terms, the performance achieves Rank-1 accuracy of 71\% and~\emph{m}AP of 68.5\%.
The results demonstrate that all these components in the SA-Softmax contribute consistently to alleviating modality discrepancy or improving discriminability through an effective training process.

\subsection{Discussion}
\textbf{Impact of the hyper-parameters $\alpha$ and $\beta$.}
In Fig.~\ref{fig:para}, we conduct empirical experiments to measure the model performance under different hyper-parameter settings.
For the discussion of $\alpha$, we choose the SA-Softmax with Feature Mask as the baseline.
As shown in the left of Fig.~\ref{fig:para} (a)(b), we observe that the performance increases significantly when attaching the SA-Softmax loss.
Even if we completely replace the softmax loss with the SA-Softmax loss, the performance under the SYSU-MM01 all-search mode still improves by 3.8\% in Rank-1 accuracy and by 4.1\% in mAP, respectively.
Besides, the performance peak is reached when $\alpha$ is set to 0.7.
From the $\beta$ side, as shown in the right of Fig.~\ref{fig:para} (a)(b), after combining the absolute-similarity term (AST) in the training phase, the Rank-1 accuracy and mAP are also improved when $\beta$ is less than 1.5.
It achieves the best performance when $\beta$ is set to 1.0.
Since a larger $\beta$ clearly lets makes the network more concerned about encouraging the intra-class similarity, the inter-class similarity will be largely ignored.
Therefore, the performance degrades when $\beta$ is greater than 1.0.
Compared to the all-search mode, the AST obtains much smaller increments in the indoor-search mode.
It may be mainly due to the variation among the testing samples.
Note that the all-search mode includes more diverse backgrounds and poses, which makes it necessary to have a more compact intra-class distribution.

\begin{table}[t]
  \centering
  \renewcommand{\arraystretch}{1.2}
  \caption{\textbf{Ablation study of weight mask (WM) in SA-Softmax in terms of CMC (\%) and \emph{m}AP (\%) on SYSU-MM01.} Here, we select the SA-Softmax with feature mask (FM) as the baseline.}
  \resizebox{85mm}{!}{
  \begin{tabular}{cc|cc|cc}
  \toprule[1pt]
   \multicolumn{2}{c|}{\textbf{Module}} & \multicolumn{4}{c}{\textbf{SYSU-MM01}} \\
   \cline{1-6}
   \multicolumn{2}{c|}{\emph{Setting}} &\multicolumn{2}{c|}{\emph{All Search}} &\multicolumn{2}{c}{\emph{Indoor Search}}\\
   \cline{1-6}
   SAS + FM & WM & R-1 & \emph{m}AP   & R-1 & \emph{m}AP\\
  \hline
   \checkmark &            & 70.3 & 66.8 & 78.3 & 81.7 \\
   \checkmark & \checkmark & 65.7 & 63.2 & 72.6 & 78.2\\
  \bottomrule[1pt]
    \end{tabular}}
    \label{Tbl:wm}
\end{table}

\begin{figure}[t]
    \centering
    \includegraphics[height=4cm,width=9.0cm]{figure4.pdf}
    \caption{\textbf{The cosine similarity between the modality prototype and the identity prototype with and without the weight mask (WM).} $P_v$ and $P_n$ are prototypes of two modalities in SA-softmax, whereas $P_s$ is the prototype of the original softmax loss.
    %
    (a) Without WM, the prototype could be supported by another prototype from the same identity. This effect preserves discriminative modality prototypes, thus providing efficient training.
    %
    (b) Due to the lack of support from another prototype, either $P_v$ or $P_n$ becomes more and more similar to $P_s$, causing the SA-Softmax loss to degrade to the original softmax loss.}
    \label{fig:wm}
\end{figure}

\textbf{Why not add the weight mask?}
Besides the feature mask, the same worry occurs in the optimization for the prototype matrix.
However, this condition is totally different from the optimization for feature embedding.
Concretely, if following the above strategy and adding a weight mask, the $L^W_{sas}$ will be formulated as:
\begin{equation}
\begin{split}
\label{eq:sasWM}
&\mathcal{L}^W_{sas}(W) = -log\frac{e^{[W^v, W^n]^{T}_{y^W_i}x_i}}{e^{[W^v, W^n]^{T}_{y^W_i}x_i}+D_{inter}(x_i)},\\
&with \quad D_{inter}(x_i) = {\textstyle \sum^{2N}_{j=1, j \ne y^W_i, j \ne y^F_i}e^{[W^v, W^n]^{T}_{j}x_i}},
\end{split}
\end{equation}
where we could approximate consider the $D_{inter}$ to the inter-class distance in the same part of softmax loss.
We can easily observe that $L^W_{sas}$ will degrade to the original softmax loss.
Although the distribution bias from the training set brings some variances between $W^v$ and $W^n$, it provides limited improvement for this task.
For example, in Fig.~\ref{fig:wm}, we visualize the cosine similarity between the two modality prototypes with the identity prototype (trained by original softmax loss), with and without the weight mask.
After adding the weight mask, the paired prototypes can support each other. This kind of support can help each prototype leave away from the identity prototype to the decision boundary. It not only avoids the degradation from the modality prototype to the identity prototype but also keeps the effectiveness of the SA-Softmax.
Furthermore, we conduct a quantitative experiment in the SYSU-MM01 to show the influence of the weight mask.
As shown in Table~\ref{Tbl:wm}, after adding the weight mask, the performance decreased in both rank-1 accuracy and \emph{m}AP.
Although benefitting from the distribution bias of different modalities in each class, the weight mask still leads to performance degradation.

\begin{table}[t]
  \centering
  \renewcommand{\arraystretch}{1.2}
  \caption{\textbf{Comparison of the SA-Softmax and other softmax-based training strategies on the SYSU-MM01 in terms of CMC (\%) and \emph{m}AP (\%).}}
  \resizebox{85mm}{!}{
  \begin{tabular}{l|cc|cc|cc}
  \toprule[1pt]
   \multirow{3}{*}{\textbf{Method}} & \multicolumn{2}{c|}{\textbf{Hyper-parameter}} & \multicolumn{4}{c}{\textbf{SYSU-MM01}} \\
   \cline{2-7}
   & \multicolumn{2}{c|}{\emph{Setting}} &\multicolumn{2}{c}{\emph{All Search}} & \multicolumn{2}{c}{\emph{Indoor Search}}\\
   \cline{2-7}
   & margin(m) & gamma($\gamma$) & R-1 & \emph{m}AP   & R-1 & \emph{m}AP\\
  \hline
    baseline (Softmax)                       & - & - & 64.5 & 61.3 & 71.6 & 76.1 \\
    AM-Softmax \cite{wang2018additive}       & 0.1 & - & 56.4 & 57.5 & 65.9 & 73.1 \\
    AM-Softmax \cite{wang2018additive}       & 0.2 & - & 55.0 & 55.9 & 63.9 & 71.0 \\
    AM-Softmax \cite{wang2018additive}       & 0.3 & - & 45.2 & 47.5 & 49.5 & 59.7 \\
    Circle Loss \cite{sun2020circle}         & - & 32 & 66.7 & 64.5 & 74.9 & 79.6\\
    Circle Loss \cite{sun2020circle}         & - & 64 & 64.6 & 62.5 & 71.4 & 76.5\\
    Circle Loss \cite{sun2020circle}         & - & 128 & 55.2 & 53.5 & 62.5 & 69.2 \\
    SA-Softmax                               & - & - & \textbf{71.0} & \textbf{68.5} & \textbf{78.8} & \textbf{82.2}\\
  \bottomrule[1pt]
    \end{tabular}}
    \label{Tbl:compare}
\end{table}

\textbf{Effect of SA-Softmax loss.} As we have introduced above, the SA-Softmax loss is not the first one to explore more specific softmax-based loss functions.
Hence, in this part, we compare SA-Softmax with two typical metric learning strategies: AM-Softmax \cite{wang2018additive} and Circle Loss \cite{sun2020circle} in the SYSU-MM01 dataset.

For a fair comparison, we use the same implementation detail as training the SA-Softmax.
For the AM-Softmax, we follow the setting provided by~\cite{wei2018person}, in which $m$ and $s$ are set to 0.3 and 15, respectively.
But we further observe that this hyper-parameter does not work well in the VI-ReID task, as a larger margin $m$ destroys the convergence in the early epochs.
Therefore, we evaluate the AM-Softmax with different $m$ of 0.3, 0.2, and 0.1, and show all the results.
Here, we measure the circle loss under the same conditions as in~\cite{sun2020circle}, but with the gamma ranging from 32 to 128 to demonstrate its performance in the VI-ReID task.

The final results are reported in Tab.~\ref{Tbl:compare}.
We can observe that those well-designed deep metric learning methods face completely different application scenarios compared to the single modality retrieval task.
We have to relax the hyper-parameters to ensure they still work.
Interestingly, adding an extra margin may not make sense in a multi-modality condition, as it would make it difficult for the network to converge in the early epochs.
Meanwhile, the SA-Softmax shows significant improvement on the VI-ReID task compared to other deep metric learning strategies.
Specifically, the SA-Softmax outperforms the circle loss by +4.3\% in Rank-1 accuracy and +4.0\% in the \emph{m}AP under the SYSU-MM01 all-search mode.
Besides, in the indoor-search mode, the SA-Softmax also surpasses the circle loss by 3.9\% in Rank-1 accuracy and 2.6\% in \emph{m}AP.
Benefiting from the design which directly emphasizes the modality discrepancy in the VI-ReID task, the SA-Softmax reaches a considerable superiority against other deep metric methods.
As a plug-and-play module, it indicates the great potential of SA-Softmax to become a new benchmark for the VI-ReID task.

\begin{figure*}[t!]
\centering
\includegraphics[height=18cm,width=18cm]{fig7.pdf}
\vspace{-1em}\caption{\textbf{Visualization of learned features space ($1_{st}$, $3_{rd}$ row) and feature distances histogram ($2_{nd}$, $4_{th}$ row) under different training strategies.}
The circles and triangles in the features space denote the features extracted from the visible and infrared modalities, respectively.
To better evaluate differences, the similarity between pairs of samples from the same modality has been ignored in the feature distance histograms.
It is obvious that the SA-Softmax loss utilizes the training set better and obtains a more compact distribution than the original softmax loss.}
\label{fig:tsne}
\end{figure*}

\textbf{Visualization of SA-Softmax loss.} To demonstrate the effectiveness of SA-Softmax, we visualize the well-learned feature spaces in the training and test sets of SYSU-MM01 via t-SNE~\cite{van2008visualizing}.
For comparison, we also visualize the feature space training by the original softmax loss.
The visualization results are shown in Fig.~\ref{fig:tsne}.
From the top column of Fig.~\ref{fig:tsne}, we can easily observe the significant modality discrepancy between the two spectral image samples, which is the main challenge of the VI-ReID task.
The middle column shows the feature space trained by the original softmax loss.
As we have previously discussed, employing the class prototype for cross-modality samples lacks the penalty for the modality discrepancy.
Therefore, the distributions in each class still show limited overlap between the two spectral samples in the training set.
Besides the limited optimization of intra-class similarity, the lack of a compact distribution also results in class prototypes that do not represent the distribution of each class well, which in turn limits the optimization of the inter-class distribution.
Such inexplicit optimizations in the training set can also lead to unsatisfactory performance in the testing set.

In contrast, as shown in the third column of Fig.~\ref{fig:tsne}, the SA-Softmax can be better adapted to the VI-ReID task, with a more specific design for the cross-modality retrieval task.
Due to the use of modality prototypes and the explicit optimization direction for the modality gap, the prototypes can better describe the distribution of each class, enabling a more efficient training process.
In the visualized feature space of the training set, SA-Softmax exhibits an impressive performance in exploiting the entire training set, which in turn leads to significant improvements in the test set.

Furthermore, to quantitatively analyze the original softmax and SA-Softmax, we visualized the cosine distance distributions of intra-class cross-modality pairs and inter-class cross-modality sample pairs in Fig.~\ref{fig:tsne}.
Compared with the softmax loss, the SA-Softmax significantly increases the similarity between intra-class cross-modality pairs and broadens the distribution discrepancy in both training and testing sets.
In general, the visualization results further demonstrate the superior performance of the SA-Softmax in VI-ReID.

\section{Conclusion}
In this paper, we highlight the limitation of the softmax loss paradigm in the visible-infrared person re-identification task: it ignores the optimization of modality discrepancy. %
To this end, we propose a novel Spectral-Aware Softmax loss (SA-Softmax), which exploits the global-wise modality information and emphasizes the modality discrepancy by employing the modality prototype with an asynchronous optimization.
Based on the analysis and observation of the ambiguous optimization in the SA–Softmax loss, we modify the SA–Softmax loss with the feature mask and absolute-similarity term to obtain a more stable result.
By combining the widely-used backbone, the SA-Softmax demonstrates its superiority among softmax-based training strategies and achieves state-of-the-art performance in two publicly available VI-ReID datasets, RegDB, and SYSU-MM01.

\section*{Acknowledgments}
This work was supported by the National Science Fund for Distinguished Young Scholars (No.62025603), the National Natural Science Foundation of China (No.U1705262, No. 62176222, No. 62176223, No. 62176226, No. 62072386, No. 62072387, No. 62072389, No. 62002305, No. 61772443, No. 61802324 and No. 61702136), Guangdong Basic and Applied Basic Research Foundation (No.2019B1515120049), the Natural Science Foundation of Fujian Province of China (No.2021J01002), and the Fundamental Research Funds for the central universities (No. 20720200077, No. 20720200090 and No. 20720200091).

\bibliographystyle{IEEEtran}


\section{Proof of Theorem \ref{mainthm}}
\label{sec:proof}

The following two results show that we need only consider groups of base size at least three in order to prove Theorem \ref{mainthm}, and give an explicit list of cases that we need to investigate.

\begin{lemma}[{\cite[Lemma 2.3]{LMM}}]
\label{no_b(G)=2}
There are no IBIS groups $G$ with non-abelian socle and $b(G)=2$.
\end{lemma}

\begin{proposition}[{\cite[Theorem 1]{MR2684423}}]
\label{b(G)>2}
The almost simple primitive groups with sporadic socle and $b(G) \geq 3$ are known.
\end{proposition}

Let $G$ be an almost simple group with sporadic socle,  acting primitively on the cosets of a maximal subgroup $H$ of $G$. In general, we expect IBIS groups to be rare, so the majority of the techniques that we employ to prove Theorem \ref{mainthm} involve showing that there exists an irredundant base of $G$ of size larger than $b(G)$, so $G$ is not IBIS.
This may be achieved in a variety of ways, some of which we now set out.  The labellings of these techniques will allow us to be more concise later on.

\begin{itemize}
\item[(T1)] We consider $G$ as a permutation group on cosets of $H$, and we construct a partial irredundant base of size $b(G)$ by random search.
\item[(T2)] We observe that the stabiliser in $G$ of a coset $Hx$ is $H^x$. Therefore,  we find a set of conjugates $\{H^{x_1}, H^{x_2},\dots H^{x_{b(G)}}\}$ of $H$ by random search such that, writing $H_i = \cap_{j=1}^i H^{x_i}$, we have a descending chain of subgroups $G>H_1> H_2> \dots > H_{b(G)} \neq 1$. This implies that $\{Hx_1, Hx_2, \dots Hx_{b(G)}\}$ is a partial irredundant base for $G$.
\item[(T3)] We choose a subgroup $K<G$, and show that $K$ has an irredundant base of size larger than $b(G)$ in its action on the cosets of $H$.
\end{itemize}

 We now embark on the proof of Theorem \ref{mainthm}.
\subsection{$\bm{G_0\in  \{\text{M}_{11},\text{M}_{12}, \text{M}_{22}, \text{M}_{23},\text{M}_{24},\text{J}_1, \text{J}_2,\text{J}_3,\text{HS},\text{McL}\}}$}

In all of these cases,  we are able to construct the permutation representation of $G$ acting on cosets of $H$ in {\sc Magma} by first constructing $G$ as a matrix group, finding its maximal subgroups and using the {\tt LMGCosetImage} function.  We prove most groups are not IBIS by applying (T1).
The remaining cases are when $(G,H,b(G))$ is one of $(\M_{11},A_6.2,4)$, $(\M_{12},\M_{11},5)$, $(\M_{22},\text{L}_3(4),5)$, $(\M_{23},\M_{22},6)$ or $(\M_{24},\M_{23},7)$; we prove they are IBIS by enumerating a set of orbit representatives on $b(G)$-tuples, and showing that each is either an irredundant base, or is not a base.

\subsection{$\bm{G_0\in  \{ \text{Co}_3, \text{Co}_2,\text{He}, \text{Suz}, \text{Fi}_{22}, \text{Fi}_{23}, \text{Fi}_{24}', \text{Ru}, \text{O'N}, \text{Co}_1,
 \text{HN}\}$}}
Here at least one of the permutation representations equivalent to $G$ acting on cosets of $H$ is not readily available.  In most cases however, we are able to construct $H$ using the {\tt AtlasSubgroup} command as part of the AtlasRep package of {\sf GAP} \cite{atlasrep,GAP4}.
We then apply (T2) in each case and deduce that $G$ is not IBIS.
\subsection{$\bm{G_0=\text{Ly}}$}
\label{Ly}
In view of Proposition \ref{b(G)>2}, we have $H\in \{G_2(5), 3.\text{McL}:2\}$, where $b(G)=3$ in each case.

Let $H=G_2(5)$. We appeal to (T3) by considering the centraliser $K=C_G(g_5)$ of an element $g_5\in H$ of order 5.  We find $H$ has two classes of elements of order 5 and,  by the centraliser orders, we deduce that both classes are contained in the class labelled 5A in $G$. Now $K=C_G(g_5)$ is maximal in $G$ \cite[p. 174]{ATLAS} and, replacing $g_5$ with a conjugate if necessary,  a construction of $K$ is available as a straight line program in the Online Atlas \cite{onlineATLAS}.
Without loss, we choose $g_5$ to be in class 5A in $H$, so that $|C_H(g_5)| = 2^3.3.5^6$, and find candidates for $C_H(g_5)<K$ using the \texttt{Subgroups} command in {\sc{Magma}}.  In order to construct an partial irredundant base for $K$ of length 3, we first fix the trivial coset $H$ so that we now must find a partial irredundant base for $K_2=C_H(g_5)$ of size 2. We then randomly search $K$ for two elements $x_1, x_2$ such that the following conditions hold:
\begin{enumerate}
\item $x_1x_2^{-1} \notin H$, so that $Hx_1$ and $Hx_2$ are distinct.
\item $K_2> K_2 \cap K_2^{x_1} >1$, and
\item  $K_2 \cap K_2^{x_1}> K_2 \cap K_2^{x_1} \cap K_2^{x_2} >1$ so that $\{H,Hx_1,Hx_2\}$ forms a partial irredundant base for $K$ by (T2).
\end{enumerate}
Choosing $K$, $K_2$ in this way allows us to verify each of these conditions in $K$ rather than $G$, which significantly reduces the computational resources required. We find an appropriate $x_1, x_2$ via random search and therefore prove that $G$ is not IBIS.

If instead $H = 3.\text{McL}:2$, we repeat the same process, instead considering the centraliser of an involution in class 2A of $H$. We deduce that $G$ is not IBIS.
\subsection{$\bm{G=\text{J}_4}$}
By Proposition \ref{b(G)>2}, we only need consider the action of $G$ on right cosets of one of its maximal subgroups $H \in \{2^{11}:\M_{24}, 2^{1+12}:3.\M_{22}.2, 2^{10}:\mathrm{L}_5(2)\}$, where $b(G)=3$ in each case.

Suppose $H=2^{11}:\M_{24}$.  Now $H$ is the stabiliser of a vector in the irreducible 112-dimensional representation of $G$ over $\mathbb{F}_2$ \cite[p.190]{ATLAS}. Generators for $G$ and $H$ in this representation, along with a vector $v$ fixed by $H$, are given in the Online Atlas \cite{onlineATLAS}.  The action of $G$ on cosets of $H$ is
then permutationally isomorphic to the action of $G$ on $v^G$.  We consider a subgroup $K \cong \M_{24}$ of $H$ and find that $v$ along with two further vectors in $v^G$ obtained by random search
have stabiliser order chain $[|M_{24}|,16,2]$ in $K$. The result follows applying (T3).

Now suppose $H=2^{1+12}:3.\M_{22}.2$. We employ technique (T3) with the centraliser of an element $g_5\in H$ of order 5.  There is a single class of elements of order 5 in $G$ and so we have $|C_G(g_5)| = 2^6.3.5.7$, while $|C_H(g_5)| = 2^6.3.5$.  Two elements $x_1,x_2 \in C_G(g_5)$ are representatives for the same coset in $G/H$ if and only if $x_1x_2^{-1} \in C_H(g_5)$. We find appropriate $x_1, x_2$ via random search such that by (T2) we can deduce that $H, Hx_1$ and $Hx_2$ form a partial irredundant base for $C_G(g_5)$, so $G$ is not IBIS.

Finally, suppose $H=2^{10}:\mathrm{L}_5(2)$.  We apply (T3) with a subgroup $K \cong L_5(2)$ of $H$. The construction of such a subgroup $K$ is available in the Online Atlas, and we find that $K$ fixes a non-zero vector $w \in V$.  As before, we consider the action of $G$ on $w^G$.  By random search, we find $h_1, h_2 \in G$ such that $H=G_{w}>G_{w,w^{h_1}}>G_{w,w^{h_1}, w^{h_2}}>1$, so $G$ is not IBIS.

\subsection{$\bm{G=\text{Th}}$}
By Lemma \ref{no_b(G)=2} and Proposition \ref{b(G)>2},  if $G$ is IBIS, then $H \in \{ {}^3\mathrm{D}_4(2):3, 2^5.\mathrm{L}_5(2)\}$ and $b(G)=3$. In each case,  straight line programs for the generators of $H$ in terms of standard generators of $G$ are available in the Online Atlas.
We proceed by applying (T3) in the same fashion as in Section \ref{Ly}.
More explicitly, if $H={}^3\mathrm{D}_4(2):3$, then we apply (T3) by considering the centraliser of an involution $g_2$ lying in the class labelled 2A of $H$. If instead $H=2^5.\mathrm{L}_5(2)$, we again consider the centraliser of a 2A involution in $H$.
\subsection{$\bm{G=\text{M}}$}
By Proposition \ref{b(G)>2}, it remains to check that the action of $G$ on cosets of $H=2.\text{B}$, where $b(G)=3$.
First fix the identity coset, whose stabiliser is $H$. If $G$ is IBIS, then every irredundant base of $H$ on $G/H$ has size 2. Let $g_{17}$ be an element of order 17. We consider $K=C_H(g_{17})$ with the view of employing (T3).  We may construct $C=|C_G(g_{17})|$ in {\sf{GAP}} as a subgroup of its maximal overgroup $(L_3(2)\times S_4(2):2).2< G$ using the Online Atlas \cite[p.234]{ATLAS}. We know from the character table of $H$  \cite[p.210--219]{ATLAS} that $|K| = 136$. Now $C$ has a single conjugacy class of subgroups of order 136, so without loss of generality, we may fix one of them and set it to be $K$.  Now $C$ is divisible by 3 and 7, while $|K|$ is not, so take $g_3, g_7\in C_G(g_{17})$ to be elements of orders 3 and 7 respectively.  Clearly $H \neq Hg_3, Hg_7$, and we must also verify that $g_3g_7^{-1} \notin H$.  We observe that $k\in K$ fixes $Hg_3$ if and only if $k^{g_3} \in H$.  Since $k^{g_3}\in C$, we must have $k^{g_3} \in H\cap C =  K$, and a similar argument holds for $Hg_7$. Therefore, we check directly in our construction of $C$ and $K$ in {\sf{GAP}} that there exist $g_3$, $g_7$ such that: (1) $g_3g_7^{-1} \notin H$ so that $Hg_3$, $Hg_7$ are distinct, (2) $K_1=K^{g_3}\cap K < K$ and (3) $1< K^{g_7}\cap K_1 < K_1$. Together, these demonstrate that $K$ has a partial irredundant base of size 2 on $G/H$, hence so does $H$. Therefore, $G$ is not IBIS.

\subsection{$\bm{G=\text{B}}$}
By Proposition \ref{b(G)>2}, we need to consider
 \[H\in \{2.{}^2.E_6(2).2, 2^{1+22}.\mathrm{Co}_2, \mathrm{Fi}_{23}, 2^{9+16}.\mathrm{PSp}_8(2), \rm{Th}, (2^2\times F_4(2)):2, 2^{2+10+20}.(M_{22}:2 \times S_3)\}\]
 The action of $G$ on cosets of $H$ has base size 3 in all cases except $H=2.{}^2.E_6(2).2$, where $b(G)=4$.

 First let $H=2.{}^2.E_6(2).2$. We can construct $H$ as the centraliser of an element $t_0$ in conjugacy class $2A$ of $G$. The orbit lengths of $H$ on 2A involutions $t\in B$ are given in \cite[p. 216]{ATLAS} and we reproduce them in Table \ref{B_M1_orbs}.
 \begin{table}[h!]
 \begin{tabular}{ccc}
 \toprule
 $i$& Class of $t_0t_i$ & $|C_H(t_i)|$\\
 \midrule
0&1A & $2^{38}.3^9.5^2.7^2.11.13.17.19$\\
 1&2B & $2^{38}.3^6.5.7.11$\\
 2&2C & $2^{26}.3^6.5^2.7^2.13.17$\\
 3&3A & $2^{18}.3^9.5^2.7.11.13$\\
 4&4B & $2^{30}.3^6.5.7$\\
 \bottomrule
 \end{tabular}
 \caption{Orbit lengths of $H$ on $2A$ involutions of $B$. \label{B_M1_orbs}}
 \end{table}

Let $t_i$ be an orbit representative of the orbit labelled $i$ in Table \ref{B_M1_orbs}. We claim that $\{t_0,t_1,t_2, t_3\}$ forms a partial irredundant base for $G$, implying $G$ is not IBIS. Clearly $G_{t_0} = H$ by definition and $|G_{t_0,t_1}|  = |C_H(t_1)| = 2^{38}.3^6.5.7.11$ by Table \ref{B_M1_orbs}.  Choosing $t_2$ and $t_3$ appropriately, we find that the Sylow 2-subgroup of $|G_{t_0,t_1,t_2}|$ has order $2^{26}$, while the Sylow 2-subgroup of $|G_{t_0,t_1,t_2,t_{3}}|$ has order $2^{18}$. Therefore, these groups are distinct and non-trivial, so the claim is proved.

Now let $H=2^{1+22}.\mathrm{Co}_2$. Similar to the last case, $H$ can be constructed as the centraliser of an involution $t_0$ in class 2B of $B$, and we can instead consider the action of $G$ on its class 2B. M\"uller \cite[Table 1]{MR2379937} computed the orbit lengths of $H$ on 2B involutions.  The smallest non-trivial orbit has length 93150 and so the size of a stabiliser $K$ in $H$ of an element $t_1$ of this orbit is larger than the number of cosets of $H$ in $G$. Therefore, it is impossible to construct an irredundant base of size 3 that begins with $\{t_0,t_1\}$, so $G$ is not IBIS.

Now let $H=\rm{Fi}_{23}$. The $H$-orbits on cosets of $H$ were computed by M\"uller et al. \cite[Table 2]{MR2331753}.  The Sylow 7-subgroup of $H$ has order 7, as do the stabilisers $K_2 \cong O_8^+(3):2_2$  and $K_3\cong \mathrm{PSp}_8(2)$ of elements in the non-trivial orbits labelled 2 and 3 in \cite[Table 2]{MR2331753} respectively. Hence, taking appropriate representatives $Hx_2, Hx_3$ of these orbits, we achieve an irredundant partial base $\{H,Hx_2,Hx_3\}$ of size 3 stabilised by an element of order 7.  Since the stabiliser is not trivial, $G$ is not IBIS.

Next suppose $H = 2^{2+10+20}(M_{22}:2\times S_3)$.
Neunh\"offer et al.  \cite[Table 1]{MR2824523} computed a subset of the $H$-suborbits of $B$ acting on cosets of $H$.  Although only the orbit lengths and sizes of the corresponding stabilisers in $H$ are available, this is sufficient to show $G$ is not IBIS.  There is a single class of elements of order 5 in $H$, and two distinct orbits of $H$ on cosets with stabiliser sizes 240 and 120. Therefore,  by choosing representatives of these orbits appropriately, we can construct a partial irredundant base of size three stabilised by an element of $G$ of order 5. Therefore, $G$ is not IBIS.

For the remaining maximal subgroups $H_4=2^{9+16}.\mathrm{PSp}_8(2)$, $H_5=\rm{Th}$ and $H_6= (2^2\times F_4(2)):2$
we adopt a different approach.  In the following construction, we will use $H$ to refer to any one of the groups $H_4,H_5$ or $H_6$. Let $g_6 \in H$ have order 6, and let $g_2=g_6^3$.  We will construct distinct cosets $Hx$ and $Hy$ such that $Hx$ is fixed by $g_2$ and $g_6$, while $Hy$ is fixed by $g_2$, but not $g_6$. This will show that $\{H,Hx,Hy\}$ form a partial irredundant base for $G$ in its coset action on $H$, and so $G$ is not IBIS.

Now,  $g$ fixes $Hx$ if and only if $g^{x^{-1}}=xgx^{-1} \in H$, and similarly for $Hy$. Hence, it is sufficient to exhibit $x$ and $y$ such that $g_6^{x^{-1}}, g_2^{x^{-1}}, g_2^{y^{-1}}\in H$, while $g_6^{y^{-1}} \notin H$.  Notice that we may simplify the problem slightly by taking $x,y\in HC_G(g_2)$. We must then check that $g_6^{x^{-1}} \in H$ $g_6^{y^{-1}} \notin H$ and that $x,y, xy^{-1} \notin H$ to ensure that $H$, $Hx$ and $Hy$ are distinct.  Since $HC_G(g_6) \subseteq HC_G(g_2)$, the probability of finding $x \in C_G(g_2)$ such that $g_6^{x^{-1}} \in H$ is proportional to
\[
\frac{|C_G(g_6)| }{ |C_G(g_2)|} = \frac{|C_G(g_6)|{|C_H(g_2)|}}{|C_G(g_2)|{|C_H(g_6)|}}.
\]
We try to maximise this quantity in our choice of the class of $g_6$, and using the Character Table Library in GAP (including the labelling of classes given there), we choose $g_6$ from classes $6C$, $6B$ and $6Q$ in $H_4$, $H_5$ and $H_6$ respectively.

We now turn to a computational search in {\sc Magma}. The relevant code to reproduce the procedure we describe is available on the author's website. We are able to construct $G$ as a matrix group in a 4370-dimensional representation over $\mathbb{F}_2$, and construct each of the maximal subgroups $H$ using straight-line programs on the standard generators, which are available in \cite{onlineATLAS}. We observe that each maximal subgroup acts reducibly on the underlying vector space $V$, so we compute the submodules preserved by each maximal subgroup.  In each case, we are able to find a relatively small submodule on which $H$ acts faithfully. We project onto this submodule so that we are able to find the conjugacy classes of $H$ using the \texttt{LMGClasses} function.  Once we have found an appropriate element $g_6$, we then find the preimage in the 4370-dimensional representation of $H$, and define $g_2=g_6^3$.  We compute $C_G(g_2)$ using \texttt{CentraliserOfInvolution} function, and search for $x$ and $y$ satisfying the above conditions by random selection.  Since determining membership of a random element of $G$ in $H$ is very time-consuming, we instead equivalently check whether the element preserves a submodule fixed by $H$. We are able to find an appropriate $x$ and $y$ in each case, so infer that $G$ is not IBIS in its coset actions on $H_4$, $H_5$ and $H_6$.

The proof of Theorem \ref{mainthm} is now complete.
\section*{Acknowledgements}
The author would like to thank Eamonn O'Brien for helpful discussions and assistance with the intensive computations required for the proof for the Baby Monster. The author also acknowledges the support of an Australian Research Council Discovery Early Career Researcher Award (project number DE230100579).

\bibliographystyle{abbrv}