\section{Introduction}
\label{sec:intro}
In the past few years, the accuracy of face recognition networks has significantly improved \cite{schroff2015facenet,taigman2014deepface,ranjan2019fast,deng2018arcface,bansal2018deep}. These improvements have led to the usage of face recognition systems in a large number of applications. This has raised concerns about bias against protected categories such as age, gender or race. A recent study performed by NIST \cite{grother2019face} found evidence that characteristics such as gender and ethnicity impact the verification and matching performance of existing algorithms. Similarly, \cite{buolamwini2018gender} showed that most face-based gender classifiers perform better on male faces than female faces.\\

Several works \cite{wang2019racial,amini2019uncovering,krishnapriya2020issues,vangara2019characterizing,nagpal2019deep,article,cavazos2019accuracy,georgopoulos2020investigating} have recently analyzed and proposed techniques to mitigate bias against race and skintone in face recognition. However, the issue of gender bias has not been widely explored. It is often assumed that gender imbalance in face recognition training datasets is a cause of gender bias in face recognition. However, \cite{albiero2020does} show that to obtain similar verification performance on male and female faces, we need to find the appropriate percentage of male and female identities (which may not be equal) in the training dataset. Finding such an appropriate mixture is not scalable.
\cite{albiero2020face} show that minimizing the effect of gendered hairstyles improves the matching scores for genuine female pairs but does not improve that of impostor female pairs.\\
\begin{figure}
\centering
{\includegraphics[width=\linewidth]{images/gb_teaser.pdf}}
\caption{\small We present a framework called AGENDA, which learns to reduce gender information from face descriptors $f_{in}$. After training, we can transform the face descriptors using AGENDA and generate their gender agnostic representations $f_{out}$, which can then be used for face recognition/verification.}
\label{fig:gbteaser}
\vspace{-0.5cm}
\end{figure}

Recent studies \cite{dhar2019attributes, hill2019deep} have also shown that face recognition networks encode gender information while being trained for identity classification. This implies that face descriptors \footnote{Face descriptors refer to the features extracted from the penultimate layer of a previously trained face recognition network.} extracted from such networks can be trained to predict face gender. This can be viewed as privacy leakage since a malicious agent can learn private gender information about an individual without receiving prior authorization. \\

Privacy leakage in face representations is an important issue, as it can allow for the unauthorized extraction of private sensitive attributes (such as race, gender, and age). This issue of estimating soft biometric attributes  in raw facial images has been studied in  \cite{fu2010age,lu2004ethnicity,makinen2008evaluation}. To prevent gender leakage, \cite{othman2014privacy} propose a face-morphing strategy to suppress gender in face images, while preserving identity information. However, there does not exist such a study with respect to post-hoc utilization of face descriptors extracted from previously trained networks. In this regard, we investigate solutions to the following problems:
\begin{itemize}[leftmargin=3em]
    \item Privacy leakage in face descriptors.
    \item Bias in face verification.
\end{itemize}

In our preliminary experiments (in Sec. \ref{subsec:prednbias}), we find that face descriptors from which the gender is more difficult to predict generally demonstrate lower gender bias in face verification tasks. Therefore, we hypothesize that \textit{reducing the ability to predict gender in face descriptors will reduce gender bias in face verification tasks}. To test this hypothesis, we present an \textbf{A}dversarial \textbf{Gen}der \textbf{D}e-biasing \textbf{a}lgorithm (AGENDA) (Fig. \ref{fig:gbteaser}) that transforms face descriptors so that they can be used to accurately classify identity but not gender. Reducing the predictability of gender in face descriptors will impede malicious gender classifiers, thus reducing the possibility of gender leakage. Although in this work we focus on gender leakage and gender bias, with sufficient training data, AGENDA could be modified to reduce the leakage of other sensitive attributes like race and age.

We define gender bias (in Eq. \ref{eq:bias}) as the difference in the face verification performance on male and female faces. As a result of reducing gender predictability (either by using a baseline or AGENDA), we find that gender bias in face verification decreases considerably. Our technique can be used as a post-hoc de-biasing measure for face recognition networks that might exhibit biases. Moreover, our proposed technique does not rely on the gender mixture (i.e. \% of male and female identities) of the training dataset.  To summarize, the contributions of our paper are as follows:
\begin{itemize}[leftmargin=*]
    \item To the best of our knowledge, we are the first to examine verification performance for gender agnostic face descriptors and analyze the subsequent reduction in gender bias.
    \item We experimentally verify that face descriptors with low gender predictability generally demonstrate lower gender bias in face verification.
    \item To further reduce gender predictability in face descriptors, we propose a method, AGENDA, that unlearns gender information in descriptors while training them for classification. Once trained, AGENDA can then be used to generate gender de-biased representations of face descriptors. When such representations are used for face verification, we find that AGENDA significantly outperforms the baseline with respect to bias reduction. Finally, we analyze the bias versus accuracy trade-off for face verification in male and female faces.
\end{itemize}

\section{Problem statement}
Our goal is to reduce gender information in face descriptors so that the ability of a classifier to predict the gender from these descriptors is reduced. From our initial experiments (section \ref{subsec:prednbias}) we find that face descriptors that show low gender predictability demonstrate relatively lower gender bias in face recognition/verification. From this observation, we hypothesize that reducing gender predictability of face descriptors will considerably reduce the gender bias demonstrated by the descriptors. At this point, we quantitatively describe gender bias in the context of face verification.  We define gender bias, at a given false positive rate (FPR) as the absolute difference between the verification performance for male-male and female-female pairs.
\begin{equation}
    \text{Bias}^{(F)} = |\text{TPR}_{m}^{(F)} - \text{TPR}_{f}^{(F)}|
\label{eq:bias}
\end{equation}
where $\text{TPR}_{m}^{(F)}$ and $\text{TPR}_{f}^{(F)}$ denote the true positive rate for the verification of male-male and female-female pairs respectively at a given FPR $F$. We present a framework that generates de-biased descriptors that, when used for face verification, obtains low bias value and reasonable TPR\textsubscript{f} and TPR\textsubscript{m} at any given FPR.

\begin{figure}
\centering
\hspace*{-0.75cm}\includegraphics[width=1.15\linewidth]{images/gb_pipeline_aaai.pdf}
\caption{\small \textbf{AGENDA architecture}. Face descriptors $f_{in}$ are extracted from a previously trained network $P$ and are fed to a model $M$. $M$ consists of a single linear layer with PReLU \cite{he2015delving} activation that outputs transformed face descriptor $f_{out}$. This is then fed to classifier $C$ and ensemble $E$. The arrows indicate the dataflow at various training stages. In stage 1, $M$ and $C$ are initialized and trained to classify identity using the gradients of $L_{class}$. In stage 2, $E$ is initialized and trained to classify gender using gradients of $L_g$. In stage 3, $M$ and $C$ are trained using the gradients of $L_{br}$ to gender-debias $f_{out}$ while simultaneously being able to classify identity. In stage 4, one member of ensemble $E$ is trained to classify gender from $f_{out}$ using the gradients of $L_g$. Stages 3 and 4 are repeated in alternating fashion, where the ensemble member of $E$ being trained in stage 4 changes at each iteration.}
\label{fig:gbpipeline}
\end{figure}
\section{Proposed approach}
\label{sec:approach}

The key idea in our proposed approach - AGENDA, is to train a model to classify identities while discouraging it to predict gender. Firstly, for a given image $I$, we extract a face descriptor $f_{in}$ using a pre-trained network $P$.
\begin{equation}
    f_{in} = P(I)
\end{equation}
The AGENDA architecture (Fig. \ref{fig:gbpipeline}) is composed of three components:\\
(1) \textbf{Generator model $M$}: A model that takes in face descriptor $f_{in}$ from a pre-trained network $P$, and generates a lower dimensional descriptor $f_{out} \in \mathbb{R}^{256}$. $M$ consists of a single linear layer with 256 units, followed by a PReLU layer. The weights of $M$ are denoted as $\phi_M$. \\
(2) \textbf{Classifier} $C$: A classifier that takes in the output of $M$($f_{out}$) and generates a prediction vector for identity classification. The weights of $C$ are denoted as $\phi_C$.\\
(3) \textbf{Ensemble of discriminators $E$}: An ensemble of $K$ gender prediction models represented as $E_1, E_2 \ldots E_K$. Each of these models is a simple MLP network with an input layer size of 256 units, a SELU activated \cite{klambauer2017self} linear layer with 128 units, and a sigmoid activated output layer with 2 units followed by softmax. We collectively denote the weights of all the models in $E$ as $\phi_E$ and weights of $k^{th}$ model $E_k$ as $\phi_{E_k}$.

We now interpret AGENDA as an adversarial approach. $M$ is a generator that should ideally generate gender agnostic descriptors $f_{out}$. $f_{out}$ is fed to the ensemble $E$ of gender prediction models which acts as a discriminator. $E$ should be able to accurately predict the gender using $f_{out}$. Our aim is to train $M$ to generate descriptors $f_{out}$ that can fool $E$ in terms of gender prediction, and can also be used to classify identities. Therefore, we need to impose two constraints on $f_{out}$: (i.) a penalty term that preserves the identity information, and (ii.) a penalty term that minimizes the gender predictability of $f_{out}$. To this end, we propose a bias reducing classification loss $L_{br}$, which is explained in section~\ref{sec:bias}.
\subsection{Bias reducing classification loss $L_{br}$}
\label{sec:bias}
\setlength{\belowdisplayskip}{1pt}
\setlength{\abovedisplayskip}{1.2pt}
After extracting the descriptor $f_{in}$ from a pre-trained face recognition network, we pass it through $M$ to obtain a lower dimensional descriptor $f_{out}$.
\begin{equation}
    f_{out} = M(f_{in}, \phi_M)
\end{equation}
\textbf{First constraint:} To make $f_{out}$ proficient at classifying identities we provide it to classifier $C$ and use cross-entropy classification loss $L_{class}$ to train both $C$ and $M$.
\begin{equation}
    \mathbf{\hat{y}_{id}} = C(f_{out}, \phi_C)
\end{equation}
\begin{equation}
\label{eq:lclass}
    L_{class}(\phi_M, \phi_C) = -\mathbf{y_{id}} . \text{log}(\mathbf{\hat{y}_{id}})
\end{equation}
$\mathbf{y_{id}}$ is a one hot identity label and $\mathbf{\hat{y}_{id}}$ is the corresponding softmaxed output of classifier $C$.
\\
\textbf{Training discriminators:} $M$ generates $f_{out}$ which is fed to ensemble $E$. Each of the gender prediction models in $E$, denoted $E_k$, are used for computing cross entropy loss $L^{(E_k)}_g$ for gender classification. $L_g$ is computed as the sum of cross-entropy losses for each $E_k$.
\begin{equation}
    L^{(E_{k})}_g(\phi_M, \phi_{E_{k}}) = - y_{g} \text{log } y_{g}^{(k)} - (1-y_{g}) \text{log } (1 - y_{g}^{(k)})
\end{equation}
\begin{equation}
\label{eq:gc}
L_g(\phi_M, \phi_E)=\sum^{K}_{k=1}L^{(E_{k})}_g
\end{equation}
$y_g$ is the binary gender label for the input face descriptor, and $y^{(k)}_g$ represents the respective softmaxed outputs of $E_k$ in the ensemble.\\
\textbf{Training generator (second constraint):} After training $E$, $M$ is trained to transform $f_{in}$ into gender agnostic descriptor $f_{out}$. We then provide $f_{out}$ to each model in ensemble $E$:
\begin{equation}
    o^{male}_{k}, o^{female}_{k} = E_{k}(f_{out}, \phi_{E_{k}})~~\text{for}~ k=1\ldots K
\end{equation}
The outputs $o_k^{*}$ represent the gender probability scores. If an optimal classifier operating on $f_{out}$ were to always produce a posterior probability of 0.5 for both male and female classes then this implies that no gender information is present in the descriptor. To this end, we define the adversarial loss $L^{(E_{k})}_a$ for the $k^{th}$ model in $E$ to be:
\begin{equation}
    L^{(E_{k})}_a(\phi_M,\phi_{E_{k}}) = - (0.5*\text{log}(o^{male}_{k}) + 0.5*\text{log}(o^{female}_{k}))
\end{equation}
Here, we use an ensemble of gender prediction models instead of a single model because we want $f_{out}$ to be gender agnostic with respect to multiple gender predictors. This approach was motivated by the work of \cite{wu2018towards} to solve `the $\forall$ challenge'. After computing the adversarial loss for model $M$ with respect to all the models in $E$, we select the one for which the loss is maximum. We term this loss as debiasing loss $L_{deb}$.
\begin{equation}
\label{eq:ldeb}
    L_{deb}(\phi_M, \phi_E) = \text{max}\{L^{(E_{k})}_a(\phi_M, \phi_{E_k}) |^{K}_{k=1} \}
\end{equation}
The idea is that we would like to penalize $M$ with respect to the strongest gender predictor for which it was not able to fool. This approach was introduced in \cite{wu2018towards}.  $L_{deb}$ is then combined with $L_{class}$ to compute a bias reducing classification loss $L_{br}$ as follows:
\begin{equation}
\label{eq:lbr}
    L_{br}(\phi_C,\phi_M,\phi_E) = L_{class}(\phi_C, \phi_M)+\lambda L_{deb}(\phi_M, \phi_E)
\end{equation}
Here, $\lambda$ is used to weight the de-biasing loss.
\subsection{Stage-wise Training}
\label{sec:stagewise}
We now explain the various stages of training AGENDA. \\
\textbf{Stage 1 - Initializing and training $M$ and $C$}: Using input descriptors $f_{in}$ from a pre-trained network, we train $M$ and $C$ from scratch for $T_{fc}$ iterations using $L_{class}$ (Eq. \ref{eq:lclass}). \\
\textbf{Stage 2 - Initializing and training $E$}: Once $M$ is trained to perform classification, we feed the outputs $f_{out}$ of $M$ to an ensemble $E$ of $K$ gender prediction models. $E$ is trained from scratch to classify gender for $T_{gtrain}$ iterations using $L_g$ (Eq. \ref{eq:gc}). $\phi_M, \phi_C$ remain unchanged in this stage.\\
\textbf{Stage 3 - Update model $M$ and classifier $C$}: Here, $M$ is trained to generate descriptors $f_{out}$ that are proficient in classifying identities and are relatively gender-agnostic. $f_{out}$ is fed to the ensemble $E$ and the classifier $C$, the outputs of which result in $L_{deb}$ (Eq. \ref{eq:ldeb}) and $L_{class}$ (Eq. \ref{eq:lclass}) respectively. We combine them to compute $L_{br}$ (Eq. \ref{eq:lbr}) for training $M$ and $C$ for $T_{deb}$ iterations, while $\phi_E$ remains locked.  While computing $L_{br}$, the gradient updates for $L_{deb}$ are propagated to $\phi_M$ and those for $L_{class}$ are propagated to both $\phi_M$ and $\phi_C$.\\
\textbf{Stage 4 - Update ensemble $E$ (discriminator)}: In stage 3, $M$ is trained to generate gender-debiased descriptors $f_{out}$ to fool the models in $E$, whereas in stage 4, models in $E$ are re-trained to classify gender using $f_{out}$. Therefore, we run stages 3 and 4 alternatively, for $T_{ep}$ episodes, after which we re-initialize and re-train all the models in $E$ (as done in stage 2). Here, one episode indicates an instance of running stages 3 and 4 consecutively. In stage 4, we heuristically choose one of the models in $E$, and train it for $T_{plat}$ iterations or until it reaches an accuracy of $G_{thresh}$ on the validation set. $\phi_M$ and $\phi_C$ remain locked in this stage. The motivation for training a single model in $E$ at a time is deferred to discussion in the supplementary material. The full details of training are described in Algorithm \ref{alg:AGENDA}.\\
\vspace{-0.5cm}
\begin{algorithm}[H]
\algsetup{linenosize=\scriptsize}
 \scriptsize
\floatname{algorithm}{Algorithm}
\caption{AGENDA}
\label{alg:AGENDA}
\begin{algorithmic}[1]
\STATE \textbf{Required}: $N_{ep}$: Number of training episodes
\STATE \textbf{Required}: Hyperparams: $\lambda, K, T_{fc}, G_{thresh}, T_{deb}, T_{gtrain}, T_{plat}, T_{ep}$\\
\STATE \textbf{Required} Learning rates: $\alpha_1, \alpha_2, \alpha_3$
\FOR {$i$ in \textbf{range}($N_{ep}$)}
\STATE Begin \textbf{Stage 1} (initial training of $M$ and $C$)
\IF{$i$ == 0}
\STATE Initialize $\phi_M$ and $\phi_C$ with random weights
\FOR {$n$ in \textbf{range}($T_{fc}$)}
\STATE $\phi_M\longleftarrow \phi_M - \alpha_1\nabla_{\phi_M}L_{class}(\phi_M,\phi_C)$
\STATE $\phi_C\longleftarrow \phi_C - \alpha_1\nabla_{\phi_C}L_{class}(\phi_M,\phi_C)$
\ENDFOR
\ENDIF
\STATE Begin \textbf{Stage 2} (initial training of $E$)
\IF {$i$ mod $T_{ep}$ == 0}
\STATE Initialize $\phi_E$ with random weights
\FOR{$n$ in \textbf{range}($T_{gtrain}$)}
\STATE $\phi_E \longleftarrow \phi_E - \alpha_2\nabla_{\phi_E}L_g(\phi_M, \phi_E)$
\ENDFOR
\ENDIF
\STATE Begin \textbf{Stage 3} (update $M$ and $C$)
\FOR{$n$ in \textbf{range}($T_{deb}$)}
\STATE $\phi_M \longleftarrow \phi_M - \alpha_3 \nabla_{\phi_M}L_{br}(\phi_C, \phi_M, \phi_E)$
\STATE $\phi_C \longleftarrow \phi_C - \alpha_3 \nabla_{\phi_C}L_{br}(\phi_C, \phi_M, \phi_E)$
\ENDFOR
\STATE Begin \textbf{Stage 4} (update $E_k$)
\STATE $k$ = $i$ mod $K$
\FOR{$n$ in \textbf{range}($T_{plat}$)}
\STATE Compute validation gender prediction accuracy $A$ of $E_k$
\IF{$A>G_{thresh}$}
\STATE break
\ENDIF
\STATE $\phi_{E_{k}} \longleftarrow \phi_{E_{k}} - \alpha_2 \nabla_{\phi_{E_k}} L^{(E_k)}_{g}(\phi_M, \phi_{E_{k}})$
\ENDFOR
\ENDFOR
\end{algorithmic}
\end{algorithm}
\section{Experiments}
\subsection{Pre-trained networks and evaluation dataset}
We separately evaluate the face descriptors, obtained from the penultimate layer of following two pre-trained networks:\\
\textbf{Arcface} : Resnet-101 trained on MS1MV2 \footnote{https://github.com/deepinsight/insightface/wiki/Dataset-Zoo} with Additive Angular margin (Arcface) loss \cite{deng2018arcface}. There are 59,563 males and 22,499 females in this dataset.\\
\textbf{Crystalface} : Resnet-101 trained on a mixture of UMDFaces\footnote{\label{note2}http://umdfaces.io/}\cite{bansal2017umdfaces}, UMDFaces-Videos\textsuperscript{\ref{note2}}\cite{bansal2017s} and MS-Celeb-1M \cite{guo2016ms}, with crystal loss \cite{ranjan2019fast}. There are 39,712 males and 18,308 females in this dataset.\\
For evaluation, we use the aligned faces in the IJB-C dataset, and follow the 1:1 face verification protocol defined in \cite{maze2018iarpa}. The alignment is done using \cite{ranjan2017all}. However, instead of verifying all the pairs, we only verify male-male and female-female pairs. There are 6.4 million male-male and 2.1 million female-female pairs defined in the protocol. The gender labels are provided in the dataset. Using Arcface and Crystalface we extract 512 dimensional face descriptors for the aligned faces in the IJB-C dataset which are then used for verification of male-male and female-female pairs. The gender-wise verification plots for Crystalface and Arcface descriptors are provided in Fig. \ref{fig:predcomp}.
\subsection{Relation between gender predictability and bias}
\label{subsec:prednbias}
We compare the predictability (i.e. ability to classify gender) of face descriptors extracted from Arcface and Crystalface networks. We train a logistic regression classifier on 60k IJB-C face descriptors (30k males and females) to classify gender and test it on 20k IJB-C descriptors (10k males and females). The images for training and testing are selected randomly, and the face descriptors are extracted using the pre-trained networks (Arcface or Crystalface). Note that there is no overlap between identities in train and test split. From the  results in Table \ref{tab:predcomp}, we find that descriptors extracted using Arcface show relatively low gender predictability than Crystalface descriptors. In Table \ref{tab:predcomp}, we also find that the gender bias (computed using Eq. \ref{eq:bias} and plots shown in Fig \ref{fig:predcomp}) is lower in most FPRs when Arcface descriptors are used for face verification, as comapred to Crystalface descriptors. This shows that \textit{face descriptors with low gender predictability appear to demonstrate lower gender bias in face verification}, thus forming the basis of our initial hypothesis (mentioned in Sec. \ref{sec:intro}). Therefore, we propose techniques to reduce the predictability of gender in face descriptors while making them proficient in identity classification.
\begin{table}
\caption{\small Gender-wise verification results and corresponding bias. Acc. refers to gender prediction accuracy of the logistic regression classifier}
\scriptsize
\centering
\hskip-0.38cm\begin{tabular}{cc|ccc|ccc|ccc|ccc}
\toprule
FPR  & &\multicolumn{3}{c|}{$10^{-6}$}& \multicolumn{3}{c|}{$10^{-5}$} & \multicolumn{3}{c|}{$10^{-4}$} & \multicolumn{3}{c}{$10^{-3}$} \\
\midrule
Network&Acc.& TPR\textsubscript{m} & TPR\textsubscript{f}& Bias & TPR\textsubscript{m} & TPR\textsubscript{f}& Bias & TPR\textsubscript{m} & TPR\textsubscript{f}& Bias & TPR\textsubscript{m} & TPR\textsubscript{f}& Bias\\
\midrule
Arcface& \textbf{76.01}& 0.82&0.74& 0.08 & 0.92 & 0.90 & \textbf{0.02} &0.96 & 0.93 & \textbf{0.03} &0.97&0.96&\textbf{0.01} \\
Crystalface& 80.50& 0.67 &0.63 & \textbf{0.04} & 0.84 &0.8 & 0.04 & 0.92& 0.87 & 0.05&0.96&0.93&0.03\\
\bottomrule
\end{tabular}
 \label{tab:predcomp}
\end{table}
\begin{figure}
\centering
{\includegraphics[width=1\linewidth]{images/combined_rg1arc_stretched2.pdf}}
\caption{\small Gender-wise verification plots for Arcface and Crystalface}
\vspace{-0.4cm}
\label{fig:predcomp}
\end{figure}
\subsection{Baseline: CorrPCA}
\label{sec:exp1}
 Since our hypothesis involves removing gender specific information from face descriptors, we propose a naive approach, termed as `Correlation-based PCA' (CorrPCA) for this task. We first compute the eigenspace of the descriptors and isolate the eigenvectors that encode gender information. After this, we remove these eigenvectors and transform the test face descriptors using the remaining subspace. A similar approach was used in the early nineties \cite{turk1991face} to reduce the impact of illumination on PCA features.\\ \textbf{Isolating gender specific components:} We first randomly select 80k (40k males and females) aligned face images from MS-Celeb-1M dataset. Then, we extract the 512-dimensional descriptors for these images, using a given pre-trained network (Arcface/Crystalface). We then compute the eigenspace $S\in \mathbb{R}^{512\times512}$ of the descriptors $X \in \mathbb{R}^{80k \times 512}$ using PCA. Using each eigenvector in $S$, we transfrom the original descriptors $X$ as follows :

 \begin{equation}
     v_s = X\cdot s
 \end{equation}
  Here, $s$ is a row in $S$ (i.e. an eigenvector of $X$), and $v_s$ is the corresponding component of each descriptor. Now, for all the 80k images, we  have a vector $v_s \in \mathbb{R}^{80k\times1}$ and a label vector $\ell$, with gender labels for all the sampled images. After this, we compute the Spearmann correlation coefficient between $v_s$ and $\ell$. We select the eigenvectors in $S$ for which this correlation is lower than $\delta$ and denote them collectively as a subspace $R$. Thus, $R$ is a subspace that has relatively lower gender information. We find that the number of eigenvectors in the subspace $R$ is 504 and 487 for the descriptors of Arcface and Crystalface respectively.\\
 \textbf{Transforming test face descriptors using remaining components:} We then extract the 512-dimensional descriptors for the IJB-C dataset using the given pre-trained network (Arcface/Crystalface), to obtain the transformed descriptors $X_{ijbc}$. Finally, we transform $X_{ijbc}$ using the subspace spanned by $R$ into a new feature space $T_{ijbc}$. We use $\delta=0.1$ in this protocol, for both networks. We evaluate the gender leakage and gender bias of the transformed descriptor sets for IJB-C dataset in Sec. \ref{subsec:res}.\\
\vspace{-0.2cm}
\subsection{Training details for AGENDA}
\label{sec:agendahp}
For training AGENDA, we use a combination of UMDFaces, UMDFaces-Videos and MS-Celeb-1M datasets. The face alignment and gender labels are obtained using \cite{ranjan2017all}. We perform our experiments using input face descriptors $f_{in}$ from Arcface and Crystalface, extracted for these datasets. The hyperparameter information for AGENDA is provided below:\\
\textbf{Stage 1}: $T_{fc}=66000$ iterations, learning rate $\alpha_1 = 10^{-5}$. \\
\textbf{Stage 2}: $T_{gtrain}=30000$ iterations, learning rate $\alpha_2 = 10^{-3}$. Here, we use $K=1$ and $5$ for Arcface and Crystalface, respectively.\\
\textbf{Stage 3}: $T_{deb} = 1200$ iterations,learning rate $\alpha_3 = 10^{-4}$. We compute $L_{br}$  using $\lambda=10$ and $1$, when using $f_{in}$ from Arcface and Crystalface respectively.\\
\textbf{Stage 4}: $T_{plat}=2000$ iterations, learning rate $\alpha_2 = 10^{-3}$ (same as stage 2). $G_{thresh}=0.90$ and $0.80$ for Arcface and Crystalface, respectively.\\
In all the aforementioned training stages, we use an Adam optimizer and a batch size of 400, and we ensure that each batch is balanced in terms of gender.
\begin{figure}
{\centering
\subfloat[]{\includegraphics[height = 0.39\linewidth]{images/combined_arcface.pdf}}
\subfloat[]{\includegraphics[height = 0.39\linewidth]{images/combined_rg1_trimmed.pdf}}
\caption{\small Gender-wise IJB-C verification results for (a.) Arcface and (b.) Crystalface using face descriptors from original pre-trained networks and their AGENDA counterparts. AGENDA reduces gender bias at low FPRs.}
\vspace{-0.4cm}
\label{fig:bias_AGENDA}
}
\end{figure}

\subsection{Results}
\label{subsec:res}
\textbf{Evaluating gender leakage}: We follow the same steps as in section ~\ref{subsec:prednbias} for computing the gender classification accuracy (Table \ref{tab:predcomp}) of raw face descriptors. We present the accuracy of the logistic regression classifier trained on face descriptors obtained from CorrPCA and AGENDA in Table \ref{tab:logacc}. We find that for both Arcface, the Crystalface, the classification accuracy goes down when the face descriptors are transformed using AGENDA or CorrPCA framework, which indicates that gender leakage from the descriptors is likely reduced.\\
\textbf{Evaluating verification performance and bias}: We first evaluate the 1:1 verification performance obtained by CorrPCA-transformed IJB-C face descriptors. From Table \ref{tab:biasres}, we can infer that this method helps to reduce the gender bias in Crystalface, whereas the bias and performance in Arcface remains mostly unchanged. Following that, we evaluate the verification performance and corresponding bias obtained after using AGENDA. After AGENDA training, we feed the 512 dimensional $f_{in}$ (extracted from Arcface/Crystalface) of aligned IJB-C images to the trained model $M$ (Fig. \ref{fig:gbpipeline}), that generates 256 dimensional gender de-biased descriptors $f_{out}$. We then use $f_{out}$ to perform gender-wise IJB-C 1:1 face verification. From Fig. \ref{fig:bias_AGENDA}. We find that when using $f_{out}$ from Arcface, the bias is reduced (especially in low FPRs), without catastrophically losing verification performance. This is done by improving female-female verification at low FPRs, while slightly decreasing male-male verification. Similarly, for $f_{out}$ from Crystalface, we find that the gender bias is reduced at low FPRs. The bias is especially close to 0 after FPR $10^{-5}$.\\
From Table \ref{tab:biasres},  we can infer that AGENDA consistently outperforms CorrPCA in terms of bias reduction at almost all the FPRs under consideration.  For FPRs not reported in Table \ref{tab:biasres}, the gender bias in verification is exactly same for CorrPCA, AGENDA and the original pre-trained network.\\
\newcommand{\tabarc}{
\begin{tabular}{cccc|ccc|ccc|ccc}
\toprule
FPR & \multicolumn{3}{c|}{$10^{-6}$}& \multicolumn{3}{c|}{$10^{-5}$} & \multicolumn{3}{c|}{$10^{-4}$} & \multicolumn{3}{c}{$10^{-3}$}\\
\midrule
Method& TPR\textsubscript{m} & TPR\textsubscript{f}& Bias & TPR\textsubscript{m} & TPR\textsubscript{f}& Bias & TPR\textsubscript{m} & TPR\textsubscript{f}& Bias & TPR\textsubscript{m} & TPR\textsubscript{f}& Bias \\
\midrule
Original& 0.82&0.74& 0.08 & 0.92 & 0.90 & 0.02 &0.96 & 0.93 & 0.03 &0.97&0.96&0.01 \\
CorrPCA& 0.82& 0.76&\cellcolor{yellow} 0.06 &0.92 & 0.90 & 0.02 &0.96 & 0.93 & 0.03&0.97&0.96&0.01\\
AGENDA & 0.81&0.79 & \cellcolor{green}\textbf{0.02}&0.90& 0.89& \cellcolor{green}\textbf{0.01} & 0.94 & 0.93 & \cellcolor{green}\textbf{0.01} &0.96&0.95&0.01\\
\bottomrule
\end{tabular}}
\newcommand{\tabrg}{
\begin{tabular}{cccc|ccc|ccc|ccc}
\toprule
FPR & \multicolumn{3}{c|}{$10^{-6}$}& \multicolumn{3}{c|}{$10^{-5}$} & \multicolumn{3}{c|}{$10^{-4}$} & \multicolumn{3}{c}{$10^{-3}$}\\
\midrule
Method & TPR\textsubscript{m} & TPR\textsubscript{f}& Bias & TPR\textsubscript{m} & TPR\textsubscript{f}& Bias & TPR\textsubscript{m} & TPR\textsubscript{f}& Bias & TPR\textsubscript{m} & TPR\textsubscript{f}& Bias\\
\midrule
Original & 0.67 &0.63 & 0.04 & 0.84 &0.8 & 0.04 & 0.92& 0.87 & 0.05&0.96&0.93&0.03\\
CorrPCA & 0.67&0.65&\cellcolor{yellow}0.02&0.84&0.81&\cellcolor{yellow}0.03&0.91&0.87&\cellcolor{yellow}0.04&0.95&0.93&\cellcolor{green}0.02 \\
AGENDA &0.32&0.32& \cellcolor{green}\textbf{0.0}&0.67&0.67&\cellcolor{green}\textbf{0.0}&0.83&0.81&\cellcolor{green}\textbf{0.02}&0.92&0.90&\cellcolor{green}\textbf{0.02} \\
\bottomrule
\end{tabular}}
\newcommand{\tabrgtpe}{
\begin{tabular}{cccc|ccc|ccc|ccc}

\toprule
FPR & \multicolumn{3}{c|}{$10^{-6}$}& \multicolumn{3}{c|}{$10^{-5}$} & \multicolumn{3}{c|}{$10^{-4}$} & \multicolumn{3}{c| }{$10^{-3}$}\\
\midrule
Method & TPR\textsubscript{m} & TPR\textsubscript{f}& Bias & TPR\textsubscript{m} & TPR\textsubscript{f}& Bias & TPR\textsubscript{m} & TPR\textsubscript{f}& Bias & TPR\textsubscript{m} & TPR\textsubscript{f}& Bias\\
\midrule
Original & 0.80&0.69& 0.11 & 0.88 & 0.84 & 0.04 &0.93 & 0.89 & 0.04 &0.96&0.94&0.02\\
CorrPCA & 0.79&0.69&\cellcolor{yellow}0.10&0.88&0.84&0.04&0.93&0.90&\cellcolor{yellow}0.02&0.96&0.94&0.02 \\
AGENDA &0.68&0.60& \cellcolor{green}\textbf{0.08}&0.82&0.80&\cellcolor{green}\textbf{0.02}&0.90&0.87&\cellcolor{green}\textbf{0.03}&0.95&0.93&\cellcolor{green}\textbf{0.02} \\

\bottomrule
\end{tabular}}
\begin{table}
\caption{\small Performance of logistic regression classifier trained using descriptors extracted from original networks, transformed using CorrPCA and extracted using AGENDA.}
\scriptsize
\centering
\begin{tabular}{c|ccc|ccc}
Network & \multicolumn{3}{c|}{Arcface} &  \multicolumn{3}{c}{Crystalface}\\
\hline
 Method & Original & CorrPCA & AGENDA& Original & CorrPCA & AGENDA \\
\hline
Gender classif\textsuperscript{n} acc. & 76.01 & 72.33 &64.01 & 80.50 & 75.74 & 67.25 \\
Females misclassified (\%) & 39.15 & 36.08 & 51.10 & 36.63 & 37.70 & 43.82 \\
Males miscalssified (\%) & 8.83 & 20.68 & 20.88 & 4.3& 10.8 & 23.31\\
\end{tabular}
 \label{tab:logacc}
\end{table}
\begin{table}%
 \caption{\small Gender-wise IJB-C verification results and their corresponding gender bias for (a.) Arcface and (b.) Crystalface. AGENDA shows lower bias than CorrPCA at most FPRs. Green: least bias, Yellow: second least bias.}%
  \scriptsize
  \centering
  \subfloat[][Arcface]{\tabarc}%
  \qquad
  \subfloat[][Crystalface]{\tabrg}

  \label{tab:biasres}%
  \vspace{-0.6cm}
\end{table}

\begin{figure*}
{\centering
\subfloat[Arcface, $\lambda=10$]{\includegraphics[width=0.25\linewidth]{images/arcK.pdf}}
\subfloat[Crystalface, $\lambda=1$]{\includegraphics[width=0.25\linewidth]{images/rg1K.pdf}}
\subfloat[Arcface, $K=5$]{\includegraphics[width=0.25\linewidth]{images/netA_lam.pdf}}
\subfloat[Crystalface, $K=5$]{\includegraphics[width=0.25\linewidth]{images/netB_lam.pdf}}
\caption{\small \textbf{Ablation experiments:} TPR (at FPR=$10^{-5}$) versus Bias used to train AGENDA on Arcface and Crystalface with different number of gender prediction models $K$ in $E$ (a, b) and with different weights $\lambda$ on $L_{deb}$ (c, d). Evaluation is done on IJBC dataset}
\vspace{-0.3cm}
\label{fig:ablation}
}
\end{figure*}
\vspace{-0.3cm}
\subsection{Ablation study - AGENDA}
Here, we evaluate two hyperparameters used for training the AGENDA framework on Arcface and Crystalface : (a.) the number of gender prediction models $K$ in the ensemble $E$ used to compute $L_{deb}$ (Eq. \ref{eq:ldeb}). (b.) the weight for $L_{deb}$ defined in Eq. \ref{eq:lbr}. We analyze how changing these hyperparameters vary the resultant bias reduction and verification performance at a fixed $\text{FPR}=10^{-5}$ in IJB-C dataset. \\
\textbf{Varying K} :We experiment with $K=1,2,3,4,5$ and $10$. Here, we fix all the other hyperparameters and use the same values specified in Sec. \ref{sec:agendahp}. In Fig. \ref{fig:ablation}(a) and (b), we find that in Arcface, changing $K$ does not have much effect on gender bias or verification TPR at FPR $10^{-5}$. However, for Crystalface, we find that as we increase $K$, the gender bias keeps decreasing which in turn leads to drop in verification performance at FPR $10^{-5}$.  We find that at $K=4$, the bias drops to 0, and as we further keep increasing $K$, the verification performance decreases. \\
\textbf{Varying $\lambda$: } We fix $K=5$ and evaluate $\lambda = 0.1, 1, 10$ for training the AGENDA framework using $f_{in}$. All the other hyperparameters use the same values specified in in Sec. \ref{sec:agendahp}. The results are presented in Fig. \ref{fig:ablation}(c) and (d). For both Arcface and Crystalface, as we keep on increasing the value of $\lambda$, the gender bias keeps generally decreasing and the verification TPR keeps  decreasing.
\begin{figure}
    \centering
    \includegraphics[width=0.4\textwidth]{images/corrAB_stretched.pdf}
  \caption{\small Distrubution of gender information in the feature space}
  \vspace{-0.6cm}
  \label{fig:cp}
\end{figure}
\subsection{Analysis of performance drop}
\label{sec:res2}
Gender is an important face attribute that helps deep networks recognize faces. So, minimizing its predictability by using AGENDA is expected to decrease the overall performance as a side-effect. As seen in Fig. \ref{fig:bias_AGENDA}, the verification performance in Crystalface reduces considerably as compared to Arcface, after AGENDA is applied. To understand this behavior, we analyze the distribution of the feature space of both Arcface and Crystalface. For this, we first randomly select 80k images (40k males and females) from the IJB-C dataset, and extract their face descriptors using a given pre-trained network (Arcface or Crystalface). Using PCA, we compute the eigenspace of this set of face descriptors. As done in the implementation of our CorrPCA baseline, we compute the gender correlation of each of the 512 eigenvectors in the eigenspace (Fig. \ref{fig:cp}).

 The top-256 eigenvectors encode more identity information than the bottom ones, since the networks are trained to classify identity. We find that the gender correlation of these identity encoding (top-256) eigenvectors of Crystalface is generally higher than Arcface. This seems to indicate that gender and identity are more entangled in face descriptors from Crystalface than Arcface. Therefore, the drop in verification performance, when the descriptors are gender de-biased is also expected to be more for Crystalface, since the verification performance depends on identity information encoded in the descriptors.
\subsection{Effect of Triplet Probabilistic Embedding}
In \cite{ranjan2019fast}, the face descriptors from Crystalface are not directly used for verification. Instead, the descriptors undergo triplet probabilistic embedding (TPE) \cite{Swami_2016_triplet} for generating a template representation of a given identity. TPE is an embedding learned to generate more discriminative, low-dimensional representations of given input descriptors, that have been shown to achieve better verification results. We apply TPE on the descriptors obtained using Crystalface and find that \textit{TPE improves the overall verification performance, but it also increases bias at all FPRs}. For comparison, we apply TPE after transforming the Crystalface descriptors with AGENDA. From Table \ref{tab:tpeb}, we can infer that the gender bias in the verification results obtained after applying TPE on AGENDA-transformed descriptors is lower than when TPE is applied on original face descriptors of Crystalface. The details of training TPE are provided in the supplementary material.
\begin{table}
\caption{\small IJB-C 1:1 verification results after applying TPE on face descriptors from Crystalface and its AGENDA counterpart}
\scriptsize
\hskip-0.3cm\begin{tabular}{cccc|ccc|ccc|ccc}
\toprule
FPR & \multicolumn{3}{c|}{$10^{-6}$}& \multicolumn{3}{c|}{$10^{-5}$} & \multicolumn{3}{c|}{$10^{-4}$} & \multicolumn{3}{c}{$10^{-3}$}\\
\midrule
Method& TPR\textsubscript{m} & TPR\textsubscript{f}& Bias & TPR\textsubscript{m} & TPR\textsubscript{f}& Bias & TPR\textsubscript{m} & TPR\textsubscript{f}& Bias & TPR\textsubscript{m} & TPR\textsubscript{f}& Bias\\
\midrule
Orig. + TPE& 0.80&0.69& 0.11 & 0.88 & 0.84 & 0.04 &0.93 & 0.89 & 0.04 &0.96&0.94&0.02 \\
AGENDA +TPE & 0.57&0.51 & \textbf{0.06}&0.75& 0.73& \textbf{0.02} & 0.88 & 0.85 & \textbf{0.03} &0.93&0.91&0.02\\
\bottomrule
\end{tabular}
\vspace{-0.5cm}
\label{tab:tpeb}
\end{table}
\section{Conclusion}
Implicit encoding of gender in face descriptors during training may lead to privacy leakage as such descriptors can be trained to classify face gender. Moreover, the expression of gender also appears to contribute to gender bias when such descriptors are used for face verification.  We address the issue of annonymizing the gender of face descriptors. Our initial experiments show that face descriptors showing lower gender predictability generally demonstrate lower gender bias in face verification. Motivated by this finding and the need for anonymizing the gender of face representations, we propose a framework - AGENDA, that adversarially reduces gender information from face descriptors, while training them to classify identities. The results of our experiments with AGENDA and CorrPCA provide evidence in support of our hypothesis that reduction of gender predictability in face descriptors will decrease gender bias in face verification. AGENDA significantly outperforms our PCA-based baseline in terms of bias reduction. The efficacy of AGENDA is evaluated on face descriptors from two SOTA recognition networks. However, gender is an important facial attribute, and we find that reducing gender information leads to slight decrease in verification accuracy. In the near future, we intend to apply and modify CorrPCA and AGENDA to reduce the information of other attributes like age and race in face recognition features; and apply a combination of several de-biasing losses to reduce the strength of multiple attributes simultaneously.
\section*{Acknowledgement}
The authors would like to thank P. Jonathon Phillips, Ankan Bansal, Aniket Roy and Rajeev Ranjan for their helpful suggestions. This research is based upon work supported by the Office of the Director of National Intelligence (ODNI), Intelligence Advanced Research Projects Activity (IARPA), via IARPA R\&D Contract No. 2019-022600002. The views and conclusions contained herein are those of the authors and should not be interpreted as necessarily representing the official policies or endorsements, either expressed or implied, of the ODNI, IARPA, or the U.S. Government. The U.S. Government is authorized to reproduce and distribute reprints for Governmental purposes notwithstanding any copyright annotation thereon.
\bibstyle{aaai21}

\section{Related work}
\textbf{Bias in face recognition:} Several empirical studies \cite{grother2019face, buolamwini2018gender, drozdowski2020demographic} have shown that many publicly available face recognition systems demonstrate bias towards attributes such as race and gender. \cite{wang2019racial, wang2020mitigating} highlight the issue of racial bias in face recognition, and propose strategies to mitigate the same. In the context of gender bias, most experiments show that the face recognition performance on females is lower than that of males. Use of cosmetics by females has been assumed to play a major role in the resulting gender bias \cite{cook2019fixed, klare2012face}. However, \cite{albiero2020analysis} show that cosmetics only play a minor role in the gender gap.\\
\textbf{Building fairer training datasets:} It has been speculated that the unequal percentage of male and female identities in training datasets might lead to gender bias in face recognition. However, \cite{albiero2020does} show that the bias is not mitigated when equal number of male and female identities are used for training. \cite{sattigeri2018fairness} uses a GAN to generate a (Celeb-A \cite{liu2015faceattributes} like) dataset which is less biased with respect to gender, for predicting attractiveness. However, such a method cannot be used to generate unbiased versions of large `in the wild datasets' like MS-Celeb-1M. In contrast to these approaches, our proposed bias mitigation algorithm does not rely on the gender mixture or quality of the training dataset.\\
\textbf{Adversarial techniques to suppress attributes:} \cite{wu2018towards} introduce an approach to anonymize identity and private attributes in a given video, while performing activity recognition. \cite{wang2019balanced} present an adversarial approach to minimize gender leakage in object classification to mitigate bias. Similarly, \cite{alvi2018turning} propose a technique to adversarially minimize the predictability of multiple attributes - gender, pose and ancestral origin, while performing age classification. \cite{li2019deepobfuscator} propose techniques to adversarially minimize gender predictability to reduce gender bias while predicting smile and presence of high-cheekbones. In some of the aforementioned experiments, the attribute under consideration is ephemeral to the target task. For example, in \cite{wu2018towards}, an action is not specific to an identity. Similarly, the presence of smile in \cite{li2019deepobfuscator} may not be unique to a single gender. In contrast, attributes like gender and race may not be ephemeral to face recognition. A given identity can be generally tied to a single gender or race. Therefore, because of the high level of entanglement between identity and gender or race, disentangling them is more involved.\\
\textbf{Gender privacy:} \cite{mirjalili2018gender,mirjalili2017soft} introduce techniques to synthesize perturbed face images using an adversarial approach so that gender classifiers are confounded, but the performance of a commercial face-matcher (in terms of similarity score) is preserved. However, such perturbations have not been demonstrated for pre-trained face descriptors.  It has been already shown in \cite{dhar2019attributes,hill2019deep} that face descriptors implicitly encode gender information during training. Therefore, a  classifier can be easily trained to predict the face gender, using these descriptors as input. Inspired by existing adversarial methods to remove certain attributes, we propose a framework to reduce the gender information in pre-trained face descriptors, while making them efficient for the task of identity classification.  \\


\section{Introduction}\label{sec:a}

Systems containing chiral fermions are very diverse and posses a number of intriguing properties that are  revealed when the topological charge is induced in them \cite{Zhitnitsky:2012ej,Kharzeev:2007tn,Kharzeev:2015znc,Li:2014bha,Gorbar:2018nmg,Marsh:2015xka,Klinkhamer:2004hg}. In particular, in the presence of the topological charge the Green's functions of the electromagnetic field acquire $P$ and $CP$-odd components \cite{Carroll:1989vb} which manifest as novel fermion interactions.

This paper considers scattering of electrically charged particles in chiral systems with finite topological charge. The scattering is mediated by electromagnetic field excitations some of which are $P$ and $CP$ even and others are odd. The latter induce novel terms in the scattering cross section that can be studied experimentally. They also affect the transport cross sections and thus the kinetic coefficients of chiral systems which are also of certain practical interest. The results of this paper pertain to any chiral system, however, for the sake of semantic simplicity it is convenient to refer to  light charged fermions as ``electrons" and to heavy charged fermions as ``ions".

The new scattering effects in chiral matter have the most clear interpretation when ``ions" are non-relativistic.  More precisely, we consider scattering of ``electrons" on heavy ``ions" assuming that the momentum transfer is much smaller than the ion mass $M$, which will be referred to as the static limit.  It is actually not difficult to obtain fully relativistic results, however they are quite bulky and hardly illuminating.

In the Born approximation, the scattering amplitude is proportional to the potential energy in the momentum space. The latter is encoded in the photon propagator which can be derived using the Maxwell-Chern-Simons theory \cite{Wilczek:1987mv,Sikivie:1984yz} that couples electrodynamics to the topological charge induced by external sources. The corresponding term in the Lagrangian is
\cite{Wilczek:1987mv,Carroll:1989vb, Sikivie:1984yz}
\ball{i1}
\mathcal{L}_A =-\frac{c_A}{4}\theta F_{\mu\nu}\tilde F^{\mu\nu}\,,
\gal
where $c_A$ is the  chiral anomaly coefficient of QED  \cite{Adler:1969gk,Bell:1969ts} and the dimensionless pseudoscalar field $\theta$ describes the topological charge. In many systems $\theta$ is believed to be a slowly varying function of coordinates and time \cite{Kharzeev:2015znc}. This is the approximation also assumed in this work. In particular, we treat the first derivative $\partial\theta$ as constant and adopt a fairly standard notation of the constant vector $b^\mu= (b_0,-\b b)=c_A\partial^\mu \theta= c_A(\dot \theta, -\b \nabla \theta)$; $b_0$ is also known as the chiral conductivity $\sigma_\chi$  \cite{Kharzeev:2009pj,Fukushima:2008xe}.
The photon propagator in chiral matter reads \cite{Carroll:1989vb,Lehnert:2004hq}
\ball{b10}
D_{\mu\nu}(q)= -i \frac{q^2 g_{\mu\nu}+i\epsilon_{\mu\nu\rho \sigma}b^\rho q^\sigma+b_\mu b_\nu}{q^4+b^2 q^2-(b\cdot q)^2}\,.
\gal
 Since we are interested in the static limit, it is convenient to introduce a notation $D_{\mu\nu}(\b q)= \lim_{q^0\to 0} D_{\mu\nu}(q)$ and the corresponding expression in the configuration space $D_{\mu\nu}(\b x)$. Throughout the paper the bold face font distinguishes the three-dimensional vectors. The potential induced by a stationary current $J^\nu(\b x)$ can be computed as
\ball{b11}
A^\mu (\b x) = -i\int d^3x' D^{\mu\nu}(\b x-\b x')J_\nu(\b x')= -i\int \frac{d^3q}{(2\pi)^3} e^{i\b q\cdot \b x}D^{\mu\nu}(\b q) J_\nu (\b q)\,.
\gal

In the forthcoming sections we consider electron scattering off potential $A^\mu$ at the leading order of the perturbation theory.  Following \cite{Qiu:2016hzd} it is instructive to consider two different types of chiral matter: (i) homogenous matter with $b_0\neq 0$, $\b b=0$ and (ii) stationary matter with $b_0=0$, $\b b\neq 0$. The paper is organized accordingly: \sec{sec:b} deals with the homogeneous chiral matter in which case the scattering cross section \eq{d5} is found to have a resonance at momentum transfer $q^2=- b_0^2$. It appears due to the periodic variation of the vector potential with the wavenumber $b_0$\cite{Qiu:2016hzd,Tuchin:2014iua} and is intimately related to the chiral instability of the electromagnetic field \cite{Carroll:1989vb}. Therefore, at $T\gg b_0$ the transport cross section $\sigma_T$ is enhanced as can be seen in \eq{d8}. As a consequence, the mean free path, which in a dilute gas of density $n$ can be estimated as  $\ell \sim 1/n\sigma_T$, is shorter than at $b_0=0$ by a factor of $\sim M^2/T^2$. This implies suppression of transport coefficients at high temperatures and, in particular, of the ratio of the shear viscosity to the entropy density is $\eta/s\sim T\ell \aver{v}$, which may be a reason for the smallness of this ratio in Quark Gluon Plasma \cite{Heinz:2013th}. Indeed, the gluon propagator has exactly same form as \eq{b10} apart from the color factor.

The stationary chiral matter is discussed in \sec{sec:d} and the corresponding scattering and transport cross sections given by \eq{c8} and \eq{c10} reflect the axial symmetry with respect of the vector $\b b$.  In \sec{sec:j} the effect of the new terms in the transport cross section is illustrated by computing the electrical conductivity  using the classical transport theory. In homogeneous matter the result is displayed in \fig{fig9} which shows suppression of the conductivity at $T\gg b_0$. In stationary matter the applied external electric field induces electric current along its direction and in the direction of $\b b$ with the corresponding conductivities $\sigma$ and $\sigma'$. Their temperature dependence is shown in \fig{fig5}. The discussion and conclusions are presented in \sec{sec:s}.

\section{Homogenous matter}\label{sec:b}

\subsection{Potential}

In homogeneous chiral matter with $\b b=0$, $b_0\neq 0$, the components of the propagator \eq{b10} read in the static limit \cite{Qiu:2016hzd}
\begin{subequations}\label{b12}
\bal
&D_{00}(\b q)= \frac{i}{\b q^2}\,,\label{b12a}\\
&D_{0i}(\b q)= D_{0i}(\b q)= 0\,,\label{b12b}\\
&D_{ij}(\b q)=-\frac{i\delta_{ij}}{\b q^2-b_0^2}-\frac{\epsilon_{ijk}q^k}{b_0(\b q^2-b_0^2)}+\frac{\epsilon_{ijk}q^k}{b_0\b q^2}\,.\label{b12c}
\gal
\end{subequations}
The current density of the static point source (the ``ion") of charge $e'$ is $J^\nu(\b x) =e' \delta\indices{^\nu_0}\delta(\b x)$. It induces the Coulomb potential
\ball{b15}
A^0(\b q)= e'/\b q^2\,,\qquad \b A(\b q)=0
\gal
implying that the scattering cross section off the point charge is given by the Rutherford formula and is not affected by the anomaly (in the static limit).

A non-trivial contribution comes about if the ``ion" is in a state $\psi$ with a finite expectation value of the magnetic moment $\b \mu$. Indeed, the spin current associated with such a state is $\b\nabla \times \psi^* \b \mu \psi$. In the point particle limit the spin current can be written as  $\b J(\b x) = \b \nabla \times (\b \mu\delta(\b x))$. It represents the first non-vanishing  multipole moment of the vector potential. Altogether the electrical current of ion is
\ball{b16}
J^0(\b x) =e'\delta(\b x)\,,\qquad
\b J(\b x) = \b \nabla \times (\b \mu\delta(\b x))\,,
\gal
which in momentum space reads
\ball{b17}
J^0(\b q)=e'\,,\qquad \b J(\b q)= i\b q\times \b \mu\,.
\gal
According to \eq{b11} and  \eq{b12} it  produces the potential
\begin{subequations}\label{b19}
\bal
A^\ell(\b q)&=-iD^{\ell i}(\b q)J_i(\b q)=  -\epsilon_{ijk}\mu^kq^j\left( -\frac{i\delta_{\ell i}}{\b q^2-b_0^2}-\frac{\epsilon_{\ell i r}q^r}{b_0(\b q^2-b_0^2)}+\frac{\epsilon_{\ell i r}q^r}{b_0\b q^2}\right)\label{b19a}\\
&
= -\frac{1}{\b q^2-b_0^2}\left[ i(\b \mu\times \b q)^\ell + \frac{b_0}{\b q^2}(\b \mu\cdot \b q q^\ell-\b q^2 \mu^\ell)\right]\,, \label{b19b}
\gal
\end{subequations}
while the time component is still given by the first equation of \eq{b15}. Transformation to the configuration space is accomplished using the integral
\ball{b20}
\int \frac{e^{i\b q\cdot \b x}}{|\b q|^2-b_0^2-i0}d^3q=
 \frac{2\pi^2}{|\b x|}e^{i|\b x|b_0}\,.
\gal
The potential,   or, more precisely, the zero-frequency component of the vector potential, in the configuration space reads
\ball{b22}
\b A(\b x)=&\frac{\b\mu \times \b x}{4\pi |\b x|^3}\left( 1-ib_0|\b x|\right)e^{ib_0|\b x|}-\frac{\b\mu}{2\pi b_0|\b x|^3}\left( 1-e^{ib_0|\b x|}+ib_0|\b x|e^{ib_0|\b x|}\right)\nonumber\\
&+\frac{(\b\mu\times \b x)\times \b x}{4\pi b_0 |\b x|^5}\left[3\left(e^{ib_0|\b x|}-1\right) -ib_0|\b x|(3-ib_0|\b x|)e^{ib_0|\b x|}\right]\,.
\gal
 In the anomaly-free limit  $b_0\to 0$ \eq{b22} reduces to the classical result
\ball{b24}
\b A(\b x)= \frac{1}{4\pi}\frac{\b \mu\times\b x}{|\b x|^3}= \frac{1}{4\pi}\int \frac{\b \nabla'\times (\b \mu\delta(\b x'))}{|\b x-\b x'|}d^3x'\,, \qquad b^\mu=0\,.
\gal
Since $\mu\sim e/M$, the magnetic contribution \eq{b22} is a relativistic correction to the Coulomb potential.

 The oscillatory behavior of the potential \eq{b22} stems from the non-dissipative nature of the anomalous current \cite{Kharzeev:2011ds}. Its imaginary part indicates the radiative instability  with respect to the pair-production \cite{Tuchin:2018sqe}. An important feature of the photon propagator  in chiral matter \eq{b10} is the emergence of the unstable modes that produce the chiral instability of the electromagnetic field. This instability originates in the momentum interval $|\b q|<b_0$ and appears as the pole in the upper-half of the complex $q^0$-plane \cite{Carroll:1989vb,Tuchin:2014iua}. In the static limit $q^0\to 0$ there is a single mode $|\b q|=b_0$ that causes the chiral instability and it appears as the singularity in \eq{b19b}.

\subsection{Cross sections}

The scattering cross section of ``electron" of charge $e$ off an ``ion" of charge $e'$ reads
\ball{c3}
\frac{d\sigma}{d\Omega'} = \frac{e^2}{8\pi^2}\left[ p\cdot A^*(\b q) \, p'\cdot A(\b q)-p\cdot p' |A(\b q)|^2+p\cdot A(\b q)\, p'\cdot A^*(\b q)+m^2|A(\b q)|^2\right]\,,
\gal
where $p=(E,\b p)$ and $p'=(E,\b p')$ are particle momentum before and after scattering and $\b q=\b p'-\b p$ is the momentum transfer. In terms of the scattering angle $\vartheta$ the momentum transfer is $|\b q|= 2|\b p|\sin\frac{\vartheta}{2}$. To be sure, the electron wave function also gets anomalous contributions, however those external leg contributions are irrelevant for the scattering problem.

 The static limit which is used in derivation of \eq{c3} requires that the energy transfer be negligible compared to (i) the momentum transfer, i.e.\ $q^0\ll |\b q|$ and (ii) the anomaly parameter, i.e.\  $q^0\ll b_0$ as can be seen by examining the denominator of \eq{b10}. Using  $q_0= \b q^2/2M$ the first condition  implies  that $|\b q|\ll M$, while  the second condition imposes a stronger constraint $|\b q|\ll \sqrt{b_0M}$ (since $b_0\ll M$). In the static limit the interaction time $1/q_0$ is much longer than the time $1/b_0$ it takes the chiral instability to develop.

Substituting the vector potential from\eq{b19b} and the scalar potential from \eq{b15} into \eq{c3} yields
\ball{d1}
\frac{d\sigma}{d\Omega'}= \frac{e^2}{4\pi^2}
\left\{
\left[ \frac{e'E}{\b q^2}+\frac{b_0(\b \mu\cdot \b q\, \b p\cdot \b q- \b p\cdot \b \mu \, \b q^2)}{(\b q^2-b_0^2)\b q^2}\right]^2 + \frac{[\b \mu\cdot (\b p\times \b p')]^2}{(\b q^2-b_0^2)^2}
-\frac{\b q^2}{4}\left[ \frac{e'^2}{\b q^4}-\frac{(\b \mu\times \b q)^2}{(\b q^2-b_0^2)^2}\left(1+\frac{b_0^2}{\b q^2}\right)\right]
\right\}\,.
\gal
Since the magnetic moments are usually randomly oriented,  Eq.~\eq{d1} needs to be averaged over its directions. Using $\aver{\mu_i\mu_j}= \mu^2\delta_{ij}/3$ one finds
\ball{d5}
\aver{\frac{d\sigma}{d\Omega'}}= \frac{e^2}{8\pi^2}\left\{
\frac{2E^2e'^2}{\b q^4}\left( 1-\frac{\b q^2}{4E^2}\right)+\frac{2\mu^2}{3(\b q^2-b_0^2)^2}\left(1+\frac{b_0^2}{\b q^2}\right) \left[ (\b p\times \b q)^2+\frac{\b q^4}{2}\right]
\right\}\,.
\gal
The first term in \eq{d5} corresponds to scattering off the Coulomb potential, while the second one involves a contribution of the magnetic moment. It is the last contribution that is dependent on the anomaly parameter $b_0$. It produces the resonant behavior at momentum transfers $\b q^2= b_0^2$. The origin of this behavior can be seen in the coordinate space, where the vector potential \eq{b22} oscillates with the wave number $b_0$. As have been already noted, $|\b q|= b_0$ is the runaway mode responsible for the chiral instability.

The transport cross section is defined as
\ball{c4}
\sigma_T= \int (1-\cos\vartheta)d\sigma = \frac{1}{2\b p^2}\int \b q^2 d\sigma \,.
\gal
Integration on the right-hand-side runs over the directions of momentum $\b p'$.
Plugging \eq{d5} into \eq{c4} and using $(\b p\times \b q)^2 = \b p^2\b q^2-\b q^4/4$ we obtain for the transport cross section
\ball{d7}
\sigma_T=
 \frac{e^2}{16\pi \b p^4}\int_0^{(2p)^2}\left\{
\frac{2E^2e'^2}{\b q^2}\left( 1-\frac{\b q^2}{4E^2}\right)+\frac{2\mu^2}{3}\frac{\b q^2}{(\b q^2-b_0^2)^2+\Gamma^4}\left(1+\frac{b_0^2}{\b q^2}\right) \left( \b p^2 \b q^2+\frac{\b q^4}{4}\right)
\right\}d\b q^2
\gal
where $\Gamma$   is related to the photon decay rate $w$ in the chiral medium as $\Gamma^2= q_0w$  where $w\sim e^2b_0$ \cite{Tuchin:2018sqe}. Actually, any process that tames the chiral instability also contributes to $w$.

Integration in \eq{d7} yields
\ball{d8}
\sigma_T=\frac{e^2}{16\pi \b p^4}\left( 4E^2e'^2L+ \frac{2\mu^2}{3}4\b p^4\mathcal{I}\right)\,,
\gal
where $L$ stands for the Landau logarithm  (the first, Coulomb term, in \eq{d7} is computed with the logarithmic accuracy) and we defined
\ball{d9}
\mathcal{I}= 1+\frac{1}{\epsilon}\left[ 2a(1+a)-\epsilon^2\right]\left( \arctan\frac{1-a}{\epsilon}+\arctan\frac{a}{\epsilon}\right)+\frac{1}{2}(1+3a)\ln \left(1+\frac{1-2a}{a^2+\epsilon^2}\right)\,,
\gal
with
\ball{d10}
a= \frac{b_0^2}{4\b p^2}\,,\qquad \epsilon= \frac{\Gamma^2}{4\b p^2}\,.
\gal
This function is displayed in \fig{fig3}.
At large momenta  \eq{d9} can be expanded at small $a$ while keeping the ratio $a/\epsilon$ fixed. With the logarithmic accuracy this gives $\mathcal{I}\approx \ln a^{-1}$. If $|\b p|\gg M$, then the anomalous term dominates and the transport cross section which reads\footnote{A more precise expression is obtained by keeping sub-logarithmic corrections in \eq{d9} which is accounted for by replacing in \eq{d12} $\ln(4\b p^2/b_0^2)\to \ln(4\b p^2/b_0^2)+ 1+\pi y+2y \arctan y-(1/2)\ln(1+1/y^2)$ with $y=b_0^2/\Gamma^2$. \label{foot}}
\ball{d12}
\sigma_T\approx \frac{e^2\mu^2}{6\pi } \ln\frac{4\b p^2}{b_0^2} \,,\qquad |\b p|\gg b_0\,, \Gamma\,.
\gal
In the opposite limit $|\b p|\ll b_0,\Gamma$, $\mathcal{I}\approx 6\b p^2b_0^2/(b_0^4+\Gamma^4)$ so that the Coulomb term dominates the transport cross section. It is noteworthy that the transport cross section is only logarithmically sensitive to the chiral conductivity $b_0$ and does not depend on the external parameter $\Gamma$ in two asymptotic regimes. The transition between these regimes at moderate momenta is the most sensitive to their values.
\begin{figure}[ht]
      \includegraphics[height=5cm]{i1.pdf}
  \caption{Function $\mathcal{I}(a)$ at $\epsilon=0.03$.}
\label{fig3}
\end{figure}

\section{Stationary matter}\label{sec:d}

In stationary matter with $\b b\neq 0$, $b_0= 0$, the propagator components are \cite{Qiu:2016hzd}
\begin{subequations}\label{b33}
\bal
&D_{00}(\b q)=\frac{i\b q^2}{\b q^4+\b b^2\b q^2-(\b b\cdot \b q)^2}\,,\label{b34}\\
&D_{0i}(\b q)=-D_{i0}(\b q)= \frac{\epsilon_{ijk}b^jq^k}{\b q^4+\b b^2\b q^2-(\b b\cdot \b q)^2}\,,\label{b35}\\
&D_{ij}(\b q)=-i\frac{\b q^2\delta_{ij}+b_ib_j}{\b q^4+\b b^2\b q^2-(\b b\cdot \b q)^2}\,.\label{b36}
\gal
\end{subequations}
The anomalous contribution arises already from the leading multipole moment sourced by $J^0=e'\delta(\b x)$. The corresponding  potential $ A^\mu (\b q)= -ie'D^{\mu 0}(\b q)$ has components
\ball{b40}
A^0(\b q) = \frac{e'\b q^2}{\b q^4+\b b^2\b q^2-(\b b\cdot \b q)^2}\,,\qquad
\b A(\b q) = \frac{ie'\b b\times \b q}{\b q^4+\b b^2\b q^2-(\b b\cdot \b q)^2}\,.
\gal
The expressions in the configuration space can be found in \cite{Qiu:2016hzd}.  This configuration is axially symmetric with respect to the vector $\b b$.

Employing \eq{b40} in \eq{c3} one derives the cross section
\ball{c8}
\frac{d\sigma}{d\Omega'} = \frac{(ee')^2}{16\pi^2}\frac{4E^2\b q^4 +4[\b p\cdot (\b b\times \b p')]^2-\b q^2[\b q^4-(\b b\times \b q)^2]}{ \b q^4(\b q+\unit q\times \b b)^4}\,,
\gal
where $\unit q= \b q/|\b q|$. Two features are the most noteworthy in \eq{c8} as compared to the Rutherford formula:  breaking of the axial symmetry by $\b b$ and disappearance of the small scattering angle divergence. The latter, explicitly seen  in \eq{c9},  is particularly significant in the kinetic theory.  We note the lack of the resonant  behavior, which is closely related to the absence of the chiral instability at $b_0=0$.

The transport cross section is obtained by substitution of \eq{c8} into \eq{c4}.  It is convenient to set up a Cartesian reference frame with the $z$-axis along the momentum $\b p$ and vector $\b b$ in the $yz$-plane so that the the coordinates of all relevant vectors can be parameterized as
\ball{c6}
\b p = |\b p|(0,0,1)\,,\qquad \b p'= |\b p|(\sin\vartheta\cos\phi, \sin\vartheta\sin\phi, \cos\vartheta)\,,\qquad \b b= |\b b|(0,\sin\chi, \cos\chi)\,
\gal
and the solid angle $d\Omega' = d\cos\vartheta d\phi$.
The angular integration  can be done analytically in an important case $|\b p|\gg |\b b|$ when the cross section is dominated by the small angle scattering. We have
\begin{subequations}
\bal
\sigma_T\approx &\frac{(ee')^2}{16\pi^2}\int_0^{2\pi}d\phi  \int_0^{\sim 1}d\vartheta^2\frac{E^2\vartheta^2+b^2\sin^2\chi\cos^2\phi}{(p^2\vartheta^2+b^2-b^2\sin^2\chi\sin^2\phi)^2} \label{c9}\\
\approx & \frac{(ee')^2}{16\pi^2}\frac{E^2}{\b p^4}\int_0^{2\pi}\ln \frac{\b p^2}{\b b^2(1-\sin^2\chi\sin^2\phi)}d\phi=  \frac{(ee')^2}{4\pi}\frac{E^2}{\b p^4}\ln \frac{\b p^2}{\b b^2(1+|\cos\chi|)}\,.\label{c10}
\gal
\end{subequations}
In the limit $\b b\to 0$ the denominator in the argument of the logarithm is replaced by the Debye mass $m_D$ squared. \fig{fig1} shows the transport cross section calculated with \eq{c8} and an approximation \eq{c10}.

\begin{figure}[ht]
\begin{tabular}{cc}
      \includegraphics[height=5cm]{sigmat1.pdf} &
      \includegraphics[height=5cm]{sigmat2.pdf}
      \end{tabular}
  \caption{The transport cross section at $\b b\neq 0$ and $b_0=0$ versus $E/b$ at $\chi=1$ (left panel) and versus $\chi$ at  $E/b=100$ (right panel). Solid line: exact result, dashed line: approximation \eq{c10}. }
\label{fig1}
\end{figure}

\section{Electrical conductivity}\label{sec:j}

The resonance in the scattering cross section at $|\b q|\approx  b_0\neq 0$ enhances the transport cross section \eq{d8} at $|\b p|\gtrsim b_0$ as seen in \fig{fig3}. Since the mean free path is inversely proportional to the transport cross section, at least in the dilute gas limit, the transport coefficients get suppressed at $T \gtrsim b_0$. In this section we employ the classical transport theory to illustrate this and other effects of the chiral anomaly on the electrical conductivity. We neglect the quantum statistical effects so that the equilibrium distribution of particles of energy $E$ at temperature $T= \beta^{-1}$  is proportional to $e^{-\beta E}$. The terms in the kinetic equations that involve the Berry curvature \cite{Son:2012zy} are likewise neglected as they are not essential for the purpose of this illustration. The main goal here is to estimate the effect of the anomaly on the collision integral.

In the linear response approximation the electric current induced by constant external electric field $\bm{\mathcal{E}}$ is given by
\ball{k1}
\b j = -e\int \b v \delta f d^3p = \frac{e^2}{T}\int f_0\nu^{-1} (\bm{\mathcal{E}}\cdot \b v) \b v d^3p\,,
\gal
where $\nu = n' v_\text{rel}\sigma_T$ is the collision rate of electrons with ions of density $n'$, $v_\text{rel}$ is the relative velocity and $\b v = \partial E/\partial \b p$.  In the chiral limit, the equilibrium distribution of electrons $f_0$ is the one of free massless fermions because of the peculiar form of their dispersion relation $(E\pm b_0)^2=(\b p\mp \b b)^2$ (see Appendix). Calculations in this section are performed in the chiral limit and furthermore it is assumed that $e=e'$ and $n=n'$ for notational simplicity. We proceed by computing the current \eq{k1} in two limiting cases.

\subsection{Homogeneous matter}

Electrical conductivity  at $\b b=0$ can be computed using the standard formula
\ball{k21}
\sigma = \frac{e^2}{3T}\int f_0 \frac{1}{n \sigma_T}d^3p \,,
\gal
where the transport cross section is given by \eq{d8}. As explained there, at  low energies/temperatures the Coulomb term dominates the transport cross section so that one recovers the textbook result $\sigma= 16\pi T/e^2 L$, where $L=\ln(T/m_D)$. At high temperatures one obtains using \eq{d12}
\ball{k23}
\sigma = \frac{\pi}{\mu^2 T \ln(T/b_0)}\,,\qquad T\gg b_0,\,\Gamma\,.
\gal
Remarkably,  this formula does not depend on the resonance width $\Gamma$ (which is essentially a cutoff) and only weakly depends on $b_0$. It shows that the electrical conductivity of homogenous chiral matter with $b_0\neq 0$ is suppressed at high temperatures by a factor of  $e^2/8\mu^2T^2\sim M^2/2T^2$, where $M$ is the ion  mass. This is exhibited in \fig{fig9}, neglecting the temperature dependance of the Landau logarithm. The suppression happens because of the resonance in the scattering cross section at the scattering angle $\vartheta= b_0/|\b p|$, which in turn can be traced back to the oscillatory behavior of the potential as explained in \sec{sec:b}.

\begin{figure}[ht]
      \includegraphics[height=5cm]{condb0.pdf}
  \caption{Electrical conductivity \eq{k21} at $\Gamma= 0.1b_0$, $M= 5b_0$ and $L=5$  (solid line). Dashed line: the Coulomb limit, dotted line: the anomalous contribution \eq{k23} with the correction mentioned in footnote \ref{foot}. In relativistic heavy-ion collisions  $b_0$ is conjectured to be on the order of 1--10 MeV \cite{Kharzeev:2015znc}. }
\label{fig9}
\end{figure}

\subsection{Stationary matter}

At finite $\b b$ (and $b_0=0$) the current is no longer isotropic. It follows from \eq{k1} that the most general relationship between the current and the electric field has form
\ball{k3}
j_i=\Sigma_{ij}\mathcal{E}_j\,,
\gal
where $\Sigma_{ij}$ is a symmetric tensor. Since the only available symmetric tensor are $\delta_{ij}$ and $\hat b^i\hat b^j$  it can be expanded as
\ball{k5}
\Sigma_{ij}= \sigma \delta_{ij}+ \sigma' \hat b^i\hat b^j\,.
\gal
These terms describe the current flowing in the direction of the electric field and in the direction of the $\theta$-field gradient $\b b$:
\ball{k6}
\b j = \sigma \bm{\mathcal{E}}+ \sigma'(\unit b \cdot \bm{\mathcal{E}})\unit b\,.
\gal
The  two conductivities read
\bal
&\sigma= \frac{1}{2}(\Sigma\indices{^k_k}-\Sigma^{k\ell}\hat b_k\hat b_\ell)= \frac{e^2}{2T}\int f_0\nu^{-1} [v^2 -(\b v\cdot \unit b)^2] d^3p\,,\label{k7}\\
&\sigma'= \frac{1}{2}(3\Sigma^{k\ell}\hat b_k\hat b_\ell-\Sigma\indices{^k_k})= \frac{e^2}{2T}\int f_0\nu^{-1} [3(\b v\cdot \unit b)^2-v^2] d^3p\,.\label{k8}
\gal
In the case $\b b=0$, the integrands of \eq{k7} and \eq{k8} are spherically symmetric so that $\aver{v_iv_j}=v^2\delta_{ij}/3$, $\sigma'=0$ and $\sigma$ is given by \eq{k21}.

In the limit $T\gg |\b b|$, the main contribution to these integrals arises from $|\b p|\gg |\b b|$ region where the integrals can be evaluated analytically. Substituting \eq{c10} into \eq{k7} and \eq{k8} we get in the chiral limit
\ball{k11}
\sigma= \frac{16\pi T}{e^2\ln (T/|\b b|)}\,,\qquad \sigma'= \frac{4\pi T}{e^2\ln^2 (T/|\b b|)}\,.
\gal
We observe the logarithmic suppression of electric current in the $\b b$ direction as compared to the electric field direction. The temperature dependence of the conductivities is displayed in \fig{fig5}. The suppression effect that the resonance at finite $b_0$ has on the magnitude of the electric conductivity can be appreciated by comparing the left panel of \fig{fig5} and \fig{fig9} at $T\gg b_0$ and $T\gg |\b b|$ respectively.

\begin{figure}[ht]
\begin{tabular}{cc}
      \includegraphics[height=5cm]{cond1.pdf} &
      \includegraphics[height=5cm]{primecond1.pdf}
      \end{tabular}
  \caption{Electrical conductivities \eq{k7} and \eq{k8} (solid lines) and their high temperature approximations \eq{k11} (dashed lines).  }
\label{fig5}
\end{figure}

\section{Discussion and Conclusions}\label{sec:s}

The main feature of the electron-ion scattering differential cross section at finite $b_0$ is the resonance, which appears as the pole of the photon propagator \eq{b10} and is intimately related to the chiral instability of the electromagnetic field. It exists  at any collision energy but its particular form depends of course on the collision kinematics. In the case of electron scattering on a heavy ion,  the resonance scattering angle is found to be  $\vartheta=b_0/|\b p|$ assuming that $\b b=0$. The same result also holds in the ultrarelativistic limit when $E\approx E'$ (or, more precisely, $|q^0|\ll 2|\b p|\sin\frac{\vartheta}{2}$).
At finite $\b b$ the denominator of \eq{b10}  is proportional to $\b q^2+\b b^2\sin^2\alpha-b_0^2$ where $\alpha$ is the angle between $\b b$ and $\b q$. Clearly, the resonance emerges at any finite $b_0$ regardless of the magnitude of $\b b$, though at large $\b b$ it is confined to small values of $\sin\alpha$. The magnitude of the cross section at the resonance depends on the mechanism that tames the chiral instability of the electromagnetic field. Nevertheless, the transport cross section at high momenta given by \eq{d12}, is not sensitive to it  and is only logarithmically sensitive to $b_0$.

It was pointed out in \cite{Khaidukov:2013sja,Qiu:2016hzd} that the third term in \eq{b12c} is topological in the sense that the corresponding contribution to the interaction energy of two current loops is proportional to their Gauss linking number. The topological term is certainly present in the scattering amplitude describing interaction of the electron current with the ion magnetic moment. However, once the polarization sum is performed to obtain \eq{c3} the topological contribution can not be easily identified and is, actually, not specially interesting in the scattering problem. Nevertheless, it is of interest in bound states, which can also be addressed using similar methods.

In conclusion, we computed the scattering and transport cross sections in chiral matter. The scattering is resonant at the scattering angles $\vartheta\approx b_0/|\b p|$, which implies suppression of the transport coefficients at temperatures higher than the chiral conductivity $T>b_0$.  At finite $\b b$ the transport is anisotropic as it acquires a component in the $\b b$-direction.  These observations provide a new avenue for experimental investigation of the chiral anomaly in chiral matter. In particular, their generalization to QCD may explain  the unusually small values of the shear viscosity to the entropy density ratio of the Quark-Gluon Plasma.

\acknowledgments
This work  was supported in part by the U.S. Department of Energy under Grant No.\ DE-FG02-87ER40371.