\section{Introduction}
\label{sec:intro}
\vspace{-0.15cm}
Acoustic monitoring is a suitable method to monitor and survey avian populations \cite{brandes2008automated,Lee} in their natural habitats. Bird species classification is an important module in such acoustic monitoring systems and directly helps in determining the avian diversity of an area of interest.
One of the major problems in designing such a bird classification framework is the limited availability of labelled training data.
Therefore, it may not be possible to utilize the data intensive classification frameworks such as deep neural networks (DNN) and convolutional neural networks (CNN) to their full potential.
Hence, there is a need to develop classification techniques which could perform well even in scarcity of labelled training data.

Leveraging recent advances in dictionary learning based hashing for nearest neighbour retrieval \cite{compact_codes,robust_hash}, we propose archetypal analysis (AA) based supervised, convex-sparse hashing for species classification. The proposed hashing framework is an alternative to existing audio classification techniques and can easily be extended to other purposes such as audio indexing. The framework is characterized by the utilization of a hashing mechanism for classification and does not require any explicit classifier or reconstruction error calculations used in various dictionary learning frameworks. Generally, hashing based nearest neighbour search can be seen as a variant of template matching where hash codes act as templates for the train/test examples. However, the use of class-specific archetypal dictionaries provides generative modelling characteristics to the proposed framework. Hence, the proposed framework can be seen as a blend of generative modelling and hash-based template matching. Hashing based approaches such as \cite{shazam} have been successful in music information retrieval. However, these approaches cannot handle the large with-in class variabilities which are quite common in bioacoustic signals. The generative nature of the proposed framework helps in effectively overcoming this issue.

The proposed hashing framework builds upon the work in \cite{compact_codes,ccse} and utilizes archetypal analysis (AA) \cite{CVPR_AA} for acoustic modelling. AA requires less amount of training data in comparison to other modelling techniques, such as Gaussian mixture models, for providing effective acoustic modelling \cite{aa_BAD}. The compressed super-frame (CSF) representation \cite{ccse} of bird vocalizations, obtained from the spectrogram (by embedding context to each frame), is given as input to the proposed framework. Each of the input CSF is hashed into a compact binary code designated as \emph{conv-codes} using class-specific AA and Bloom filters \cite{bloom} (see Section \ref{sec:prop}). These conv-codes are clustered in a class-specific manner using K-medoids with Jaccard coefficient \cite{jac} as similarity metric. The cluster centers obtained from this procedure are used for populating the hash table. The proposed framework utilizes the similarity between test conv-codes and these cluster centers for classification. Moreover, both conv-codes and cluster centers are bit strings, hence efficient bit-level operations can be used to calculate the aforementioned similarity. In this work, Jaccard coefficient \cite{jac} is used as a similarity measure for classification.

In addition, this study also proposes a min-hash \cite{min} variant of the proposed framework. This variant approximates the Jaccard coefficient between two conv-codes using min-hashes extracted from convex-sparse representations. The utilization of min-hashing makes it possible to use the direct addressing \cite{direct} instead of searching the hash table for classification (see Section \ref{ssec:min} for details). The direct addressing reduces the hash table search to a constant time ($O(1)$), increasing the computational efficiency of the proposed framework.

The rest of the paper is organized as follows:
Section \ref{sec:back} discuss some of the methods proposed in the literature for bird species classification using acoustic data.
In Section \ref{sec:prop}, the proposed framework and its min-hash variant is described in detail.
Performance analysis and conclusion are provided in Sections \ref{sec:perf} and \ref{sec:con}, respectively.

\vspace{-0.25cm}
\section{Related Works}
\label{sec:back}
\vspace{-0.15cm}
Various methods targeting bird species classification have been proposed in the literature.
In some classical methods, the syllables (the basic unit of birdsong) modelled as frequency and amplitude modulated sinusoidal pulses are used for identifying bird species \cite{Harma,Somervuo,syllable_pair}.
However, the use of sinusoidal modelling based classification is limited to the species producing tonal sounds. Other methods based on support vector machines (SVM), deep neural networks (DNN) and convolutional neural networks (CNN) have also been proposed for this task.
In \cite{Raich}, a convolutional neural network (CNN), having Segnet \cite{segnet} like architecture, is proposed to segment vocalizations from the spectrogram and simultaneously classifies these bird vocalizations.
Hence, this framework overcomes the problem of segmentation of bird vocalizations before classification.
However, this network requires pixel-wise labelling of spectrograms, which is generally not available. An unsupervised and scalable feature learning method based on spherical K-means is proposed in \cite{dan_skmeans}.
These unsupervised features are shown to provide good classification performance for a large number of bird species, using a random forest classifier.
In addition, SVM with various dynamic kernels such as probabilistic sequence kernel (PSK), GMM supervector (GMMSV) kernel, GMM-UBM mean interval (GUMI) kernel, GMM-based intermediate matching kernel (GMM-IMK) and GMM-based pyramid match kernel (PMK) have also been used for bird species classification \cite{deep}.
In a recent study \cite{ccse}, AA is utilized to obtain the convex representations for bird species classification. This framework shows comparable performance to the existing classification frameworks including DNN and SVM.

\begin{figure*}[t]
	\centering
	\includegraphics[scale=0.43]{hashing_new.pdf}
	\caption{Illustration of the proposed convex sparse hashing approach for two
species. A dictionary of 20 archetypes is used for obtaining
convex representations. First ten archetypes ((1-10) correspond to
species 1 while remaining ten archetypes (11-20) correspond to species
2.
The cluster centers obtained by class-specific K-medoids clustering
(T=2) are used to obtain the hash table having 4 entries (2 per
species).}
	\label{fig:prop}
\end{figure*}
\vspace{-0.35cm}
\section{Proposed Framework}
\label{sec:prop}
\vspace{-0.15cm}
In this section, we describe the convex-sparse audio hashing in detail. First, we explain the process for obtaining archetypal dictionaries and the hash table. Then, the process to classify a given audio recording is described.
\vspace{-0.25cm}
\subsection{Preprocessing and dictionary learning}
\vspace{-0.1cm}
The audio recordings are converted into compressed super-frame (CSF) \cite{ccse} representation by concatenating $W$ neighbouring frames of the spectrogram and projecting the concatenated vector on a random Gaussian matrix. The concatenation helps in effectively capturing the frequency-temporal modulations that give species-specific nature to bird vocalizations, while random projections help in compressing the concatenated vectors ($Wm$-dimensional to $K$-dimensional, where $m$ is the number of frequency bins in a frame and $K<<Wm$). More details about CSF representation can be found in \cite{ccse}. Only CSFs corresponding to the vocalization regions of an audio are used for learning dictionaries.
A semi-supervised segmentation method proposed in \cite{thakur_mlsp} is used here for segmenting bird vocalizations.

In this work, archetypal analysis (AA) \cite{CVPR_AA} is employed to obtain the class-specific dictionaries from the CSF representation. The CSFs obtained from vocalizations of a species are pooled together to form a matrix $\mathbf{X} \in \mathbb{R}^{K \times l}$, where $K$ is the dimensionality of the CSFs and $l$ is the number of pooled CSFs.
This matrix $\mathbf{X}$ is factorized as: $\mathbf{X=DA}$.
The dictionary $\mathbf{D} \in \mathbb{R}^{K\times d}$ has $d$ archetypes that model the convex hull of the data and characterize the extremal rather than the average behaviour as provided by other acoustic modelling techniques such as Gaussian mixture models. These archetypes are the convex combination of input data points such that $\mathbf{D=XB}$. The archetypal dictionary, $\mathbf{D}$ can be obtained by solving the following optimization problem:
\vspace{-0.25cm}
\begin{equation}
\begin{aligned}
 \argmin_{\substack{\mathbf{B},\mathbf{A}\\ \mathbf{b}_j \in \Delta_l,\mathbf{a}_i \in \Delta_d}} &\Vert \mathbf{X}\mathbf{-DA} \Vert_{F}^2= \argmin_{\substack{\mathbf{B},\mathbf{A}\\ \mathbf{b}_j \in \Delta_l,\mathbf{a}_i \in \Delta_d}}\Vert \mathbf{X}\mathbf{-XBA} \Vert_{F}^2,\\
 \Delta_l \triangleq [\mathbf{b}_j\succeq 0, &\Vert\mathbf{b}_j\Vert_1=1], \Delta_d \triangleq [\mathbf{a}_i\succeq 0, \Vert\mathbf{a}_i\Vert_1=1].
\label{eq:aa}
\end{aligned}
\end{equation}

Here $\mathbf{a}_i$ and $\mathbf{b}_j$ are columns of $\mathbf{A} \in \mathbb{R}^{d\times l}$ and $\mathbf{B} \in \mathbb{R}^{l\times d}$, respectively. This optimization objective can be solved using the block-coordinate descent scheme. More details about the implementation of AA can be found in \cite{CVPR_AA}. The final dictionary $\mathbf{D}_f$ is obtained by concatenating all the class-specific dictionaries: $\mathbf{D}_f=[\mathbf{D}_{1} \mathbf{D}_{2} \ldots\mathbf{D}_{q}]$, where $\mathbf{D}_q$ is the archetypal dictionary of the $q$th class.

\vspace{-0.3cm}
\subsection{Generating conv-codes and populating hash table}
\label{ssec:conv_codes}
\vspace{-0.1cm}

A convex sparse representation ($\mathbf{y}_i$) is obtained by projecting a CSF ($\mathbf{x}_i$) on a simplex corresponding to the dictionary, $\mathbf{D}_f \in \mathbb{R}^{K \times qd}$ ($qd$ is the number of archetypes from all classes):
\vspace{-0.35cm}

\begin{equation}
\mathbf{y}_i=\argmin_{\substack{\mathbf{y}_i\\ \mathbf{y}_i \in \Delta_{qd}}} \Vert \mathbf{x}_i\mathbf{-D}_f\mathbf{y}_i \Vert^2_{2}
\label{eq:active_Set}
\end{equation}
where $\Delta_{qd} \triangleq [\mathbf{y}_i\succeq 0, \Vert\mathbf{y}_i\Vert_1=1]$.

In this work, the active-set algorithm proposed in \cite{CVPR_AA} is used to solve equation \ref{eq:active_Set}. The convexity constraints on the decomposition leads to sparsity such that only few coefficients of $\mathbf{y}_i$ are significant \cite{CVPR_AA}. The archetypes corresponding to these significant coefficients have maximum contribution in representing the CSF. The set formed by indices of the top $Z$ coefficients having maximum magnitudes in convex sparse representation  is termed as \emph{effective-set}. An effective-set is computed for each CSF extracted from training audio recordings. Ideally, a characteristic or fixed set of archetypes contribute to the possible effective-sets for vocalizations of each bird species. Hence, effective-sets obtained for different bird species should be different, making them a suitable representation for the species classification task.

The set operations, required to work with effective-sets, are usually computationally expensive. On the other hand, the bit-level operations or manipulations in binary strings are relatively computationally efficient.
Hence, these effective sets are converted into compact binary strings or conv-codes using Bloom filters \cite{bloom}. Bloom filters are designed to facilitate the set membership queries in an efficient manner by storing the information about the presence/absence of given elements in compact binary strings. To obtain a conv-code for a given effective set using Bloom filters, initially all bits in conv-code are set to 0. Then, the non-cryptographic hash functions of Bloom filters hashes each element of an effective set to one of the locations of the conv-code, where the corresponding bit is flipped to 1, marking the presence of that element in the given effective-set. These hash functions are a sequence of bit-level operations applied on an input element to produce a digest that corresponds to one of the location on the conv-code.

In this work, conv-codes extracted from the training audio recordings are used for populating the hash table. The conv-codes of each class are clustered using K-medoids (with Jaccard coefficient as the similarity metric) to obtain $T$ cluster centers. In K-medoids clustering, an input data point is the medoid or the cluster center. Hence, the cluster centers obtained using K-medoids are subset of the training conv-codes. These cluster centers are used as keys, pointing to the respective species or class label, in the hash table. Thus, the hash table has $qT$ entries, where $q$ is the number of classes. It must be noted that $qT$ is significantly less than the number of input CSFs. Figure \ref{fig:prop} depicts the process of creating the hash table from training input audio recordings of two bird species.

\vspace{-0.3cm}
\subsection{Classification}
\label{ssec:test}
\vspace{-0.1cm}
In this section, we discuss the classification strategy employed for a test audio recording.
Initially, a test input audio recording is processed to segment the bird vocalizations regions from the background.
A segmented bird vocalization is represented by a set of CSFs, $\mathcal{X}=[\mathbf{x}^t_1\mathbf{x}^t_2\ldots \mathbf{x}^t_n \ldots \mathbf{x}^t_s]$, where $\mathbf{x}^t_n$ is the $n$th CSF in the test vocalization.
A convex-sparse representation ($\mathbf{y}^t_n$) is obtained for CSF ($\mathbf{x}^t_n$) using equation \ref{eq:active_Set}.
This CSF, $\mathbf{x}^t_n$, is represented by an effective-set obtained by choosing the indicies of the top $Z$ coefficients having maximum magnitudes in $\mathbf{y}^t_n$.
Finally, $\mathbf{x}^t_n$ is converted into a conv-code ($\mathbf{C}^t_n$) using the bloom filter.
During classification, this conv-code is matched with all the conv-codes used for populating the hash tables during training. The hash table entry that exhibits maximum similarity with $\mathbf{C}^t_n$ is regarded as its nearest neighbour and $\mathbf{C}^t_n$ is assigned the species label pointed by the respective hash table entry.
In this work, we have used Jaccard coefficient \cite{jac} as a metric for this matching. The Jaccard coefficient between the test conv-code and the hash table entry compute the extent of intersection between their corresponding effective-sets.
As discussed earlier, the possible effective-sets for each class consist of a particular set of archetypes and ideally, this set of archetypes is unique for each bird species.
Hence, two effective-sets belonging to same species exhibits maximum intersection and their corresponding conv-codes exhibit maximum Jaccard similarity.

Each test CSF in $\mathcal{X}$ is classified into one of the available bird classes. A voting rule is applied on these CSF wise decisions to get the final species label of the bird vocalization represented by $\mathcal{X}$.

\vspace{-0.4cm}
\subsection{Reducing computation using min-hashing}
\label{ssec:min}
\vspace{-0.1cm}
The retrieval time required to find a best matched conv-code from the hash table is $O(qT)$, where $qT$ is the number of hash table entries. Although bit-manipulation operations are fast, the required computation to find the best matched conv-code can be reduced further. To decrease the retrieval time, we propose a variant of the proposed framework that utilizes min-hashing \cite{shrivastava2012fast,min} on the convex-sparse representations. To obtain the min-hash, convex-sparse is randomly permuted and the indices of $Z$ most significant coefficients in this permuted representations are noted. The first of these indices act as min-hash for the convex-sparse representation under processing. It has been shown in the literature that the probability of two sets having same min-hash is equal to the expected Jaccard similarity between these sets \cite{mining, shrivastava2012fast}. Hence, it can be inferred that two sets having high similarity are most likely to produce the same min-hash.

Based on this, we propose to use direct addressing \cite{direct} for creating a hash table whose keys are min-hashes and values are corresponding label information. An array (a data-structure) can be seen as over-simplified example of the direct addressing. The index of an array is a key and the value corresponding to this key can be obtained by directly accessing the memory corresponding to this index. In this work, min-hashes are equivalent to indices of array and value at each index is the corresponding species label.
Since min-hashes are indices of coefficients of the convex-sparse representation, the size of this array is equal to the number of archetypes in the dictionary. For testing a CSF, min-hash is obtained for its convex representation. The hash table is directly accessed at the location of this min-hash to obtain the species label. This hash table look-up requires $O(1)$ computation only.

\vspace{-0.4cm}
\section{Performance Analysis}
\label{sec:perf}
\vspace{-0.2cm}
\subsection{Dataset Used}
\vspace{-0.1cm}
The proposed framework is evaluated on a dataset\footnotemark \footnotetext{We are in process of making the dataset public.} containing audio recordings of 50 different bird species, obtained from three different sources.
The recordings of 26 bird species were obtained from the Great Himalayan national park (GHNP), in north India. These recordings were used in \cite{deep} for performance comparison. The recordings of 7 bird species were obtained from the bird audio database maintained by the Art \& Science centre, UCLA \cite{data2}.
The remaining 17 bird species audio recordings were obtained from the Macaulay Library \cite{mac}.
All the recording used here are 16-bit mono, sampled at 44.1 kHz and are of durations varying from $15$ seconds to $3$ minutes.
The information about these 50 species along with the total number of recordings and vocalizations per species is available at \url{https://goo.gl/z6UEQa}.

\begin{figure*}[t]
	\centering
	\includegraphics[trim={0cm 6.5cm 1cm 0cm},clip,scale=0.42]{hash_exp1.eps}
	\caption{Classification performances of different methods obtained with (A) 33.33\% and (B) 20\% of the training data (Averaged across three and five folds respectively).}
	\label{fig:exp1}
\end{figure*}
\vspace{-0.3cm}
\subsection{Experimental setup}
\vspace{-0.1cm}
The spectrograms are obtained from each input audio recording by using the STFT with 512 FFT points on a frame rate of 20 ms with 50\% overlap.
The window size $W=5$ is used to obtain super-frame representations which results in 1285-d ($1285=(257 \times 5)$) representation.
Random projections are used to compress this 1285-d representation to a 500-d representation.
A 25-atom archetypal dictionary is learned for each class and thus the final dictionary having 1250 (\textit{50}$\times $\textit{25}) atoms is used to obtain the convex sparse representations.
The parameter $Z=4$ is used for obtaining effective-sets from these convex representations.
A Bloom filter with 2 different non cryptographic hash functions i.e. Murmur3 \cite{murmur} and SpookyHash \cite{spooky} is used to convert the effective sets into 1024-bits conv-codes. K-medoids with 10 cluster centers (i.e. $T=10$) is applied on training conv-codes of each class to obtain the hash table entries. After training, the hash table containing 500 distinct entries ($10\times50$) is obtained.  All the parameters mentioned in this section are obtained empirically on a validation set.

\vspace{0.15cm}
\noindent\textbf{Comparative methods:}
The classification performance of the proposed framework and its min-hash variant is compared with various existing bird species classification methods such as CCSE framework \cite{ccse}, SVM with dynamic kernels (IMK, PMK, GUMI and PSK), spherical K-means with random forest based framework proposed in \cite{dan_skmeans} and a DNN based approach proposed in \cite{deep}.
CCSE framework is also based on super-frames and class-wise archetypal analysis, but it uses reconstruction error for classification.
For SVM and DNN based classification schemes, MFCC using delta and acceleration coefficients are used as the feature representation.
In addition, we have also compared the proposed framework with incoherent dictionary hashing (IDH) \cite{compact_codes} for the application of species classification, using nearest neighbour retrieval.
For a fair comparison, an incoherent dictionary of 1250 atoms (same as the one used in the proposed framework) is used in IDH. Apart from that, the performance of the proposed framework is also compared with its min-hash variant (MHV).

\vspace{0.1cm}
\noindent\textbf{Train/test data distribution:}
The 20\% of the vocalizations present in the audio recordings of each class are used for the validation purposes. These vocalizations are not used for training or testing.
The remaining vocalizations are used for performance comparison using a three-fold cross validation scheme. In each fold, 33.33\% of the vocalizations (per class) are used for training while the remaining are used for testing.
The results presented here are averaged across three folds.
\vspace{-0.3cm}

\subsection{Results and Discussion}
\vspace{-0.1cm}
The classification performance of the proposed framework along with other comparative methods is shown in Fig.~\ref{fig:exp1}(A).
Following inferences can be made from this figure:
\begin{itemize}[leftmargin=*]
\item The classification performance of the proposed framework is comparable to the DNN i.e. DNN shows a relative improvement of only 2.1\%.
\item Similar classification performances are observed for CCSE and the proposed framework, with CCSE showing a relative improvement of 0.9\% over the proposed approach.
This observation indicates that maximum intersection between effective-sets of two super-frames is same as the minimum class/dictionary specific reconstruction error used for classification, as in the CCSE approach.

\item The proposed framework outperforms IMK, PMK, GUMI, PSK and SKM by a relative improvement of 3.44\%, 2.76\%, 4.13\%, 2.3\% and 0.6\%, respectively.

\item The incoherent dictionary learning used in IDH could not perform up-to the level of AA used in the proposed framework.
This is due to the fact that the proposed framework utilizes class-specific AA dictionaries which provide more discriminative representations.

\item Min-hash variant shows comparable performance to SVM based frameworks. However, on comparison with the proposed framework, min-hash variant shows a drop of 3.84\% in classification accuracy. As discussed earlier, this drop in classification accuracy occurs due to the loss of information incurred while representing effective-sets by one index i.e. min-hash only. However, using min-hashing decreases the required classification time by approximately three times.

\end{itemize}
\vspace{-0.1cm}

Although the proposed method performs similar to CCSE framework, however it is computationally efficient. The minimum reconstruction error classification used in CCSE is computationally expensive in comparison to the indexing. In our experiments, it has been observed that average time taken for classifying a test bird vocalization using the proposed hashing framework is $0.13$ seconds as opposed to $0.4$ seconds in CCSE. Similarly, for the min-hash variant, this average classification time is $0.05$ seconds. Hence, min-hash variant makes up for the drop in accuracy by providing rapid classification. The running times mentioned here are the average of ten runs for classifying approximately $3200$ vocalizations.

\vspace{0.1cm}
\noindent\textbf{Performance in low data conditions:} To analyse the performance of the proposed framework and other comparative methods in low data conditions, we conducted the same experiment, with the same parameter setting but with different train/test data distribution.
Here, we have used only 20\% of the available recordings per class for training and the rest for testing.
In this experiment, a $5$-fold cross-validation is used for performance comparison.
The results of this experiment are illustrated in Fig.~\ref{fig:exp1}(B).
It can be observed that all the methods exhibits a drop in average performance when compared with the previous experiment (using 33.3\% recordings for the training).
However, this drop is more significant for methods that uses classifiers such as DNN, SVM and random forest.
The relative drops of 15.34\%, 14.4\%, 16.05\%, 14.1\%, 10.69\% and 16.03\% are observed for IMK, PMK, GUMI, PSK, SKM and DNN, respectively.
On the contrary, dictionary learning methods (CCSE, IDH, the proposed framework and min-hash variant) observe a relatively low drop in average classification performance.
The relative drops of 8.13\%, 8.04\%, 7.93\% and 8.18\% are observed for CCSE, IDH, the proposed framework and min-hash variant. This justifies the claim that the proposed framework requires less amount of training data to provide effective classification as compared to the SVM and DNN based methods.

\vspace{-0.5cm}
\section{Conclusion}
\label{sec:con}
\vspace{-0.17cm}
In this work, we proposed an AA based convex sparse hashing framework for bird species classification.
Bloom filters are used to convert convex representations into binary conv-codes, that are further used for indexing and classification.
The experiments demonstrate that the classification performance of the proposed framework is similar to state-of-the art methods.
Our experimentation has validated that this framework requires less training data in comparison to the SVM and DNN based classification frameworks.
No extra effort has been made to make the conv-codes more compact and to make the retrieval time more efficient.
In future, we will work on this aspect to further decrease the time required for classification.

\bibliographystyle{IEEEbib}


\section{Introduction}

In noncentrosymmetric metals, spin and momentum of electrons are coupled such that the deformation of the Fermi surface due to dissipative current can polarize electron spins and this is called the magnetoelectric effect or the Edelstein effect~\cite{Levitov1, Edelstein1, Fiebig}. More interestingly, in some noncentrosymmetric superconductors, supercurrents can also give rise to a static magetization without dissipation~\cite{Fiebig, Levitov2, Edelstein2, Edelstein3, Yip, Fujimoto}. And conversely, a static magnetization can drive local supercurrents~\cite{Yip, Samokhin, Agterberg, Feigelman, Sigrist1, Mironov}. As pointed out by Edelstein~\cite{Edelstein2, Edelstein3}, in superconductors with a polar axis $\bm{c}$ and Rashba spin-orbit coupling (SOC), the magnetization $\bm{M}$ induced by supercurrents can be expressed as $\bm{M}\propto\bm{c}\times\bm{J}^{\textrm{S}}$~\cite{Edelstein2} where $\bm{J}^{\textrm{S}}$ is the supercurrent density. In such Rashba superconductors, the induced magnetization is always perpendicular to the direction of supercurrent and the polar axis, as is shown in Fig. \ref{figure1} (a).

In recent years, many noncentrosymmetric superconductors such as superconductors with chiral lattice structures~\cite{Salamon, Hirata, Prozorov, Carnicom}, and superconducting transition metal dichalcogenides (TMDs)~\cite{Jianting, Saito, Mak, Cava, Sajadi, Sanfeng, Guangtong} are discovered. Since the newly emergent noncentrosymmetric superconductors have different crystal symmetries, the SOC in these materials have different forms. For example, superconductors with point group symmetry T and O, which belong to chiral (or enantiomorphic) point groups, have isotropic SOC in the form of $v\bm{p}\cdot\bm{\sigma}$, where $\bm{p}$ is momentum, $\bm{\sigma}$ is the Pauli matrices and $v$ is the SOC strength~\cite{Guoqing}. On the other hand, multilayer 1Td-structure WTe$_2$ and MoTe$_2$ possess anisotropic SOC~\cite{Guangtong}. For atomically thin monolayer 2H-MoS$_2$ and 2H-NbSe$_2$,  Ising SOC~\cite{Jianting, Saito, Mak, Noah, Benjamin, Wenyu} is present which pins electron spins to the out-of-plane directions. These different forms of SOC are expected to cause unconventional magnetoelectric effects different from the one in Fig. \ref{figure1} (a) for Rashba superconductors. However, a general understanding of the magnetoelectric effects of all these noncentrosymmetric superconductors is lacking. In this work, through linear response theory and group theory analysis, we provide a general and powerful way to understand magnetoelectric effects in noncentrosymmetric superconductors. Importantly, we point out that in superconducting chiral crystals and superconducting quasi-2D TMDs, applying supercurrent can give novel magnetization shown in Fig. \ref{figure1} (b) and (c) respectively.

First of all, we note that among the superconductors within the 21 noncentrosymmetric point groups, the magnetoelectric effect is generally non-zero only for the ones belonging to the 18 gyrotropic point groups~\cite{gyrotropy, Moore}, we call these superconductors gyrotropic superconductors.  The explicit forms of the magnetoelectric pseudotensors for the 18 gyrotropic point groups are listed in Table \ref{magnetoelectric_pseudotensor}. According to Table \ref{magnetoelectric_pseudotensor}, many of the novel magnetoelectric response of noncentrosymmetric superconductors can be identified immediately.

For example, for materials belonging to the T and O point groups, such as Li$_2$Pd$_3$B, Li$_2$Pt$_3$B~\cite{Salamon, Hirata}, and Mo$_3$Al$_2$C~\cite{Prozorov}, the magnetoelectric effect is purely longitudinal, meaning that the induced spin magnetization is always parallel to the direction of the supercurrent, as is shown in Fig. \ref{figure1} (b). This is a quantum analogue of classical solenoids in which the induced magnetic field is parallel to the current directions. We further point out that, for quasi-two-dimensional materials, the only symmetries which allow to have induced spin magnetization perpendicular to the atomic plane are the ones which have the polar axis lying inside the atomic plane. The newly discovered multilayer 1Td-structure WTe$_2$~\cite{Sajadi, Sanfeng} and MoTe$_2$~\cite{Guangtong} naturally fulfils this condition so that there exists the in-plane supercurrent induced out of plane magnetization seen from Fig. \ref{figure1} (c).

On the other hand, Ising superconductors such as 2H-structure MoS$_2$ and NbSe$_2$, the magnetoelectric effect is indeed zero due to their D$_{3h}$ point group, even though the SOC in these materials are particularly strong. These 2H-TMDs are interesting examples of noncentrosymmetric superconductors which give rise to zero magnetoelectric response. Importantly, under uniaxial strain which reduce D$_{3h}$ to C$_{2v}$, a purely out-of-plane magnetization can be induced by a supercurrent in 2H-TMD. Therefore, our results, as summerized in Table \ref{magnetoelectric_pseudotensor}, can be used to provide guiding principles to generate and manipulate magnetoelectric responses in superconductors.

\begin{figure*}
\centering
\includegraphics[width=\textwidth]{Figure1.pdf}
\caption{The supercurrent driven magnetization in (a) polar superconductors with Rashba SOC, (b) superconducting chiral crystals belonging to T and O point groups, and (c) quasi-2D superconducting TMDs with polar axis lying inside the atomic plane. In polar superconductors with Rashba SOC, the supercurrent driven magnetization $\bm{M}$ is the transverse response as $\bm{M}\propto\bm{c}\times\bm{J}^{\textrm{S}}$. In superconducting chiral crystals of T and O point group, the magnetization arising from supercurrent is parallel to the supercurrent direction and the magnetoelectric effect becomes a longitudinal response. In quasi-2D superconducting TMDs with C$_1$, C$_2$, C$_{1v}$ or C$_{2v}$ point group symmetry, the polar axis $\bm{c}$ lies inside the atomic plane and allows the in-plane supercurrent to induce the out of plane magnetization.}\label{figure1}
\end{figure*}

The rest of the paper is organised as follows. First, we construct the microscopic model for the magnetoelectric effect and show that the induced spin magnetization is related to the supercurrent by a rank-two pseudotensor $T_{ij}$. Second, we analyse the symmetry properties of $T_{ij}$ under different point group operations and list the general form of $T_{ij}$ in Table \ref{magnetoelectric_pseudotensor}. Third, we discuss the application of Table \ref{magnetoelectric_pseudotensor} in the understanding of several interesting superconductors with gyrotropic point group, including the superconductor with chiral lattice symmetry and superconducting TMDs. Finally, we demonstrate how an unconventional magnetoelectric response can be generated by strain in non-gyrotropic superconductors using 2H-structure TMD as an example.

\section{Results}

\subsection{Theory for the superconducting magnetoelectric effect}
To study a noncentrosymmetric superconductor which carries a supercurrent, we consider the Bogliubov de Gennes Hamiltonian
\begin{align}
H\left(\bm{p}, \bm{q}\right)=&\begin{pmatrix}
\xi_{\bm{p}+\frac{1}{2}\bm{q}}+\bm{g}_{\bm{p}+\frac{1}{2}\bm{q}}\cdot\bm{\sigma} & i\Delta_{\bm{q}}\sigma_y\\
-i\Delta_{\bm{q}}^\ast\sigma_y & -\xi_{-\bm{p}+\frac{1}{2}\bm{q}}-\bm{g}_{-\bm{p}+\frac{1}{2}\bm{q}}\cdot\bm{\sigma}^\ast
\end{pmatrix}.
\end{align}
Here the Nambu basis is $\left[c^\dagger_{\bm{p}+\frac{1}{2}\bm{q}, \uparrow}, c^\dagger_{\bm{p}+\frac{1}{2}\bm{q}, \downarrow}, c_{-\bm{p}+\frac{1}{2}\bm{q}, \uparrow}, c_{-\bm{p}+\frac{1}{2}\bm{q}, \downarrow}\right]$ with $c^{\left(\dagger\right)}_{\bm{p}', s}$ the annihilation (creation) operator at the momentum $\bm{p}'$ and spin index $s=\uparrow/\downarrow$.  The kinetic energy of an electron with momentum $\bm{p}$ is described by $\xi_{\bm{p}}$. The vector $\bm{g}_{\bm{p}}$ describes the SOC of the material and its specific form is determined by the lattice symmetry, and $\bm{\sigma}$ denotes the Pauli matrices. From the definition of the basis, it is clear that the $\Delta_{\bm{q}}\sigma_{y}$ term pair electrons with net momentum $\bm{q}$ and opposite spin. Therefore, the Hamiltonian describes a superconductor which sustains a supercurrent when $\bm{q}$ is finite. With the Hamiltonian, we can calculate the spin magnetization $\bm{M}$ as ~\cite{Gorkov}
\begin{align}\label{Magnetization_formula}
M_k=-\frac{1}{2\beta V}\sum_{\bm{p}, \bm{q}, n}\frac{\partial}{\partial h_k}\textrm{tr}\log G^{-1}\left(\bm{p}, \bm{q}, i\omega_n\right)|_{\bm{h}=0}.
\end{align}
Here $\beta$ reads $\beta^{-1}=k_bT$ with $k_b$ the Boltzmann constant, $k=x, y, z$ denotes the component of $\bm{M}$, and the Matsubara Green's function is defined as
\begin{align}
G^{-1}\left(\bm{p}, \bm{q}, i\omega_n\right)=i\omega_n-H\left(\bm{p}, \bm{q}\right)-\frac{1}{2}g\mu_{\textrm{B}}\begin{pmatrix}
\bm{h}\cdot\bm{\sigma} & 0 \\
0 & -\bm{h}\cdot\bm{\sigma}^\ast
\end{pmatrix},
\end{align}
with $\bm{h}$ describing a Zeeman field, $\mu_{\textrm{B}}$ the Bohr magneton, and $g$ the Lande $g$ factor. Since the supercurrent associated finite momentum $\bm{q}$ couples with the velocity operator as $\nabla\left(\xi_{\bm{p}}+\bm{g}_{\bm{p}}\cdot\bm{\sigma}\right)\cdot\frac{1}{2}\bm{q}$ at the linear order, the spin magnetic moment is induced through the rank-two pseudotensor $J_{ij}\left(\bm{p}\right)=\partial_{p_j}g_i\left(\bm{p}\right)$ with $i, j = x, y, z$. As a result, in the presence of supercurrent, the coupling term $\sigma_iJ_{ij}\frac{1}{2}q_j$ acts as an effective Zeeman field that generates  the spin magnetization.

In order to obtain an analytical form for the supercurrent induced magnetization $\bm{M}$, we expand $\textrm{tr}\log G^{-1}\left(\bm{p}, \bm{q}, i\omega_n\right)$ in terms of the Zeeman field $\bm{h}$ and the pairing order parameters $\Delta_{\bm{q}}$ and $\Delta^\ast_{\bm{q}}$. Then by taking the derivative on the Zeeman field $\bm{h}$ at $\bm{h}=0$ and summing over $\bm{p}$, $\bm{q}$, the Matsubara frequencies $\omega_n$, we obtain the analytical form for the magnetization $M$, which can be written as $M_k=T_{kj}\tilde{J}^{\textrm{S}}_{j}$, where  $\bm{J}^{\textrm{S}}$ is the supercurrent density. We denote $T_{kj}$ as a component of the magnetoelectric pseudotensor $T$~\cite{SM} with the form:
\begin{widetext}
\begin{align}\label{Key}
T_{kj}&=\left\{\frac{\Delta_0^2}{\bm{g}^2_{\bm{p}_{\textrm{F}}}}f\left(\frac{|\bm{g}_{\bm{p}_{\textrm{F}}}|}{\pi k_bT_c}\right)\left\langle J_{kj}\left(\bm{p}_{\textrm{F}}\right)\right\rangle+\left[\frac{7\Delta_0^2\zeta\left(3\right)}{8\pi^2k_b^2T_c^2}-\frac{\Delta_0^2}{\bm{g}^2_{\bm{p}_{\textrm{F}}}}f\left(\frac{|\bm{g}_{\bm{p}_{\textrm{F}}}|}{\pi k_bT_c}\right)\right]\left\langle\hat{\bm{g}}_{\bm{p}_{\textrm{F}}k}\hat{\bm{g}}_{\bm{p}_{\textrm{F}}i}J_{ij}\left(\bm{p}_{\textrm{F}}\right)\right\rangle\right\}\frac{2\mu_{\textrm{B}}g}{3\sqrt{3}\xi}N\left(E_{\textrm{F}}\right).
\end{align}
\end{widetext}
Here,  $\bm{p}_{\textrm{F}}$ represents the Fermi momentum, $N\left(E_{\textrm{F}}\right)$ is the density of states at the Fermi level, $\Delta_0$ is the homogeneous pairing order parameter magnitude, $\xi$ is the superconducting coherence length, and $\tilde{\bm{J}}^{\textrm{S}}=\bm{J}^{\textrm{S}}/J^{\textrm{S}}_{\textrm{max}}$ is the supercurrent density normalized by the maximum supercurrent density $J^{\textrm{S}}_{\textrm{max}}=2e\Delta_0^2\frac{2\hbar}{3\sqrt{3}m^\ast\xi}$ with $m^\ast$ being the effective pairing mass~\cite{Tinkham}. We have $\left\langle\right\rangle=\int\frac{d\Omega_{\bm{p}_{\textrm{F}}}}{4\pi}$ as the angle average at the Fermi surface and define the function $f\left(\rho\right)$ as
\begin{align}
f\left(\rho\right)=\textrm{Re}\sum_{n=0}^\infty\left(\frac{1}{2n+1}-\frac{1}{2n+1+i\rho}\right).
\end{align}
The magnetoelectric susceptibility obtained in Eq. \ref{Key} applies for superconductors with a pair of spin splitted bands. For noncentrosymmetric superconductor with multiple bands, the total magnetoelectric susceptibility is the sum of contribution from all the paired Fermi pockets as is discussed in the Supplementary Materials~\cite{SM}. As the magnetoelectric pseudotensor $T$ is constructed from the SOC pseudovector $\bm{g}_{\bm{p}_{\textrm{F}}}$ at the Fermi surface, its general form is determined by the crystal point group symmetry as shown in the next section.

\subsection{Symmetry Analysis for Magnetoelectric Pseudotensor in Three Dimensions}
The SOC pseudovector $\bm{g}_{\bm{p}}$ under the point group operation $\hat{R}$ respects $\bm{g}\left(\bm{p}\right)=\det\left(\hat{R}\right)\hat{R}\bm{g}\left(\hat{R}^{-1}\bm{p}\right)$~\cite{Samokhin2, Smidman}. Hence inheriting from the pseudovector $\bm{g}_{\bm{p}}$, we show that the magnetoelectric pseudotensor $T_{kj}$ under the crystal symmetry is subject to the constraints~\cite{SM}
\begin{align}\label{Transformation}
T=\det\left(\hat{R}\right)\hat{R}T\hat{R}^{\textrm{T}},
\end{align}
where $\hat{R}$ is the orthogonal matrix of the point group transformation. The general derivation of Eq. \ref{Transformation} for all the noncentrosymmetric point groups is given in the Materials and Methods section. All the nonzero components of the magnetoelectric pseudotensor $T$ of the 18 gyrotropic point groups are listed in Table \ref{magnetoelectric_pseudotensor}.

\newcommand{\tabincell}[2]{\begin{tabular}{@{}#1@{}}#2\end{tabular}}
\begin{table}[ht]
\caption{List of Magnetoelectric pseudotensor $T$ for the noncentrosymmetric superconductors in the 18 gyrotropic point groups. $T_{ij}$ with $i, j=x, y, z$ are in general the elements in $T$. $T_{\parallel}$ denotes the elements $T_{xx}$, $T_{yy}$ as $T_{\parallel}=T_{xx}=\pm T_{yy}$  and $T_0$ denotes the diagonal elements in T and O point group when  $T_{xx}=T_{yy}=T_{zz}$. $T_{\textrm{d}}$ and $T^-$ denote the symmetric and anti-symmetric off diagonal elements in $T$ respectively.The principal axis is along the $z$-direction.}
\centering
\begin{tabular}{c c c c c c}
\hline\hline
Point group & $T_{ij}$ & Point group & $T_{ij}$\\ [0.5ex]
\hline
C$_1$ & $\begin{pmatrix}
T_{xx} & T_{xy} & T_{xz} \\
T_{yx} & T_{yy} & T_{yz} \\
T_{zx} & T_{zy} & T_{zz}
\end{pmatrix}$ & C$_2$ & $\begin{pmatrix}
T_{xx} & T_{xy} & 0 \\
T_{yx} & T_{yy} & 0 \\
0 & 0 & T_{zz}
\end{pmatrix}$  \\
C$_3$ & $\begin{pmatrix}
T_{\parallel} & -T^- & 0 \\
T^- & T_{\parallel} & 0 \\
0 & 0 & T_{zz}
\end{pmatrix}$ & C$_4$ & $\begin{pmatrix}
T_{\parallel} & -T^- & 0 \\
T^- & T_{\parallel} & 0 \\
0 & 0 & T_{zz}
\end{pmatrix}$ \\
C$_6$ & $\begin{pmatrix}
T_{\parallel} & -T^- & 0 \\
T^- & T_{\parallel} & 0 \\
0 & 0 & T_{zz}
\end{pmatrix}$ & C$_{1v}$ & $\begin{pmatrix}
0 & T_{xy} & 0 \\
T_{yx} & 0 & T_{yz} \\
0 & T_{zy} & 0
\end{pmatrix}$ \\
C$_{2v}$ & $\begin{pmatrix}
0 & T_{xy} & 0 \\
T_{yx} & 0 & 0 \\
0 & 0 & 0
\end{pmatrix}$ & C$_{3v}$ & $\begin{pmatrix}
0 & -T^- & 0 \\
T^- & 0 & 0 \\
0 & 0 & 0
\end{pmatrix}$ \\
C$_{4v}$ & $\begin{pmatrix}
0 & -T^- & 0 \\
T^- & 0 & 0 \\
0 & 0 & 0
\end{pmatrix}$ &  C$_{6v}$ & $\begin{pmatrix}
0 & -T^- & 0 \\
T^- & 0 & 0 \\
0 & 0 & 0
\end{pmatrix}$ \\
D$_{2d}$ & $\begin{pmatrix}
T_{\parallel} & 0 & 0 \\
0 & -T_{\parallel} & 0 \\
0 & 0 & 0
\end{pmatrix}$ & S$_4$ & $\begin{pmatrix}
T_{\parallel} & T_{\textrm{d}} & 0 \\
T_{\textrm{d}} & -T_{\parallel} & 0 \\
0 & 0 & 0
\end{pmatrix}$ \\
D$_{2}$ & $\begin{pmatrix}
T_{xx} & 0 & 0 \\
0 & T_{yy} & 0 \\
0 & 0 & T_{zz}
\end{pmatrix}$ & D$_3$ & $\begin{pmatrix}
T_{\parallel} & 0 & 0 \\
0 & T_{\parallel} & 0 \\
0 & 0 & T_{zz}
\end{pmatrix}$ \\
D$_4$ & $\begin{pmatrix}
T_{\parallel} & 0 & 0 \\
0 & T_{\parallel} & 0 \\
0 & 0 & T_{zz}
\end{pmatrix}$ & D$_6$ & $\begin{pmatrix}
T_{\parallel} & 0 & 0 \\
0 & T_{\parallel} & 0 \\
0 & 0 & T_{zz}
\end{pmatrix}$ \\
T & $\begin{pmatrix}
T_0 & 0 & 0 \\
0 & T_0 & 0 \\
0 & 0 & T_0
\end{pmatrix}$ & O & $\begin{pmatrix}
T_0 & 0 & 0 \\
0 & T_0 & 0 \\
0 & 0 & T_0
\end{pmatrix}$ \\[1ex]
\hline
\end{tabular}
\label{magnetoelectric_pseudotensor}
\end{table}

\subsection{Unconventional Magnetoelectric Effects in Gyrotropic Superconductors}
With the general form of $T_{ij}$ in Table \ref{magnetoelectric_pseudotensor}, some unconventional and novel magnetoelectric responses of gyrotropic superconductors can be identified immediately. One particular interesting case is for materials with point group symmetry T and O. In this case, the tensor $T_{ij}$ is proportional to the identity matrix and implies that the induced magnetization is parallel to the supercurrent direction shown in Fig. \ref{figure1} (b), namely, $\bm{M} \propto\tilde{\bm{J}}^{\textrm{S}}$. This is a quantum analogue of classical solenoids but without the need to fabricate any helical structures and the longitudinal response is induced by SOC. It is interesting to note that a few superconductors with relatively high $T_c$ have point group O such Li$_2$Pt$_3$B, Li$_2$Pd$_3$B and Mo$_3$Al$_2$C~\cite{Salamon, Hirata, Prozorov}. Another interesting examples are superconducting TaRh$_2$B$_2$ and NbRh$_2$B$_2$ which had been newly discovered by Carnicom et al.~\cite{Carnicom}. These superconductors have point group C$_3$. From Table \ref{magnetoelectric_pseudotensor}, we immediately realize that a supercurrent along the $z$-direction (the principal symmetry axis direction) will generate a pure longitudinal response such that the supercurrent-induced magnetization is parallel to the direction of the supercurrent. On the other hand, a supercurrent perpendicular to the $z$-direction generates a magnetization in the $x$-$y$ plane.

Importantly, in the normal state, these materials with chiral lattice symmetries are Kramers Weyl Semimetals which has Weyl points pinned at time-reversal invariant momenta~\cite{Guoqing}. Therefore, using Table \ref{magnetoelectric_pseudotensor}, one can immediately identify the novel magnetoelectric properties of a superconducting Kramers Weyl semimetal. More noncentrosymmetric superconductors with gyrotropic point groups are listed in the Supplementary Materials~\cite{SM}.

Another interesting result from Table \ref{magnetoelectric_pseudotensor} is that, in quasi-two-dimensions, only C$_1$, C$_{1v}$, C$_{2v}$ and C$_2$ symmetries with polar axis lying inside the atomic plane allow an out-of-plane magnetization induced by an in-plane supercurrent, as is present in Fig. \ref{figure1} (c). Interestingly, recently discovered few layer 1Td-WTe$_2$~\cite{Sajadi, Sanfeng} and 1Td-MoTe$_2$~\cite{Guangtong} have such a low symmetry: C$_{1v}$. Moreover, due to the large SOC in these materials, the magnetoelectric effect is expected to be strong. With C$_{1v}$ symmetry, we show below how the magnetization direction can also be controlled by the direction of the supercurrent. The magnetoelectric effect of 2H-structure TMDs are also discussed in the next section.

\begin{figure}
\centering
\includegraphics[width=3.6in]{Figure2.pdf}
\caption{The magnetization direction (red arrows) can be controlled by the in-plane supercurrent directions $\tilde{\bm{J}}^{\textrm{S}}$ as calculated using realistic parameters for bilayer 1Td-MoTe$_2$. The in-plane mirror invariant line is defined along the $\tilde{\bm{J}}_x^{\textrm{S}}$ direction.}\label{figure2}
\end{figure}

\subsection{Unconventional magnetoelectric effects of superconducting TMDs}
In this section, we apply the general theory obtained above to two interesting TMDs, namely, bilayer or multilayer 1Td-MoTe$_2$ and strained monolayer 2H-NbSe$_2$. For the case of bilayer 1Td-MoTe$_2$, the crystal structure belongs to the C$_{1v}$ point group and only has one in-plane mirror symmetry with $y\rightarrow-y$. The superconductivity occurs at the critical temperature $T_c=4$K and it has the in-plane anisotropic upper critical field $H_{c2}$ which exceeds the Pauli limit $H_p$ due to its anisotropic SOC $\bm{g}\left(p_x, p_y\right)=\left(v_xp_y, v_yp_x, v_zp_y\right)$. From the angle dependent in-plane $H_{c2}\left(\phi\right)$, the anisotropic SOC strength at the Fermi level is estimated to be $\bm{g}\left(\phi\right)=\left(16.9\sin\phi, 23.5\cos\phi, 23.1\sin\phi\right)$meV~\cite{Guangtong}. With the Fermi wave vector $\left(p_{\textrm{F}x}, p_{\textrm{F}y}\right)=\left(1.009, 1.050\right)\textrm{nm}^{-1}$ and the Fermi energy density of states $N\left(E_{\textrm{F}}\right)=4.3776\textrm{eV}^{-1}\textrm{nm}^{-2}$ from the first principle calculation~\cite{Guangtong}, we can further use the pairing gap $\Delta_0=1.76k_bT_c$ at zero temperature and the coherence length $\xi=20.79$nm ~\cite{Guangtong} to estimate the magnetoelectric pseudotensor $T$ as $T_{xy}=2.52\times10^{-4}\mu_{\textrm{B}}/\textrm{nm}^2$, $T_{yx}=2.64\times10^{-4}\mu_{\textrm{B}}/\textrm{nm}^2$ and $T_{zy}=3.43\times10^{-3}\mu_{\textrm{B}}/\textrm{nm}^2$ while all other elements are zero.

From the form of $T$, we note that the supercurrent induced magnetization $\left(M_x, M_y, M_z\right)$ is a function of the in-plane supercurrent direction shown in Fig.\ref{figure2}. When the supercurrent flows along $y$ direction, the $z$ direction magnetization $M_z$ gets the optimised value. For a supercurrent close to the critical value, the maximum $M_z$ can reach $M_z=7.6\times10^{-4}\mu_{\textrm{B}}$ per unit cell. The strength of supercurrent induced magnetization is in the same order as Bi/Ag bilayers with strong Rashba SOC~\cite{Johansson1, Johansson2} which has been observed \cite{Fert}. Therefore, we expect that 1Td-MoTe$_2$ is an excellent candidate to observe the superconducting magnetoelectric effect which has yet to be observed experimentally. Moreover, our prediction of the supercurrent induced magnetization also applies to the recently realized superconducting 1Td-WTe$_2$~\cite{Sajadi, Sanfeng}.

For the case of 2H-NbSe$_2$ which is also known as an Ising superconductor~\cite{Jianting, Saito, Mak, Benjamin, Wenyu}, the materials have point group symmetry D$_{3h}$~\cite{Noah}.  Due to the in-plane mirror invariant line at $x=0$ and the three-fold rotation symmetry, the SOC takes the form $g_z\left(\bm{p}\right)\propto p_x\left(p_x^2-3p_y^2\right)$ around the $\Gamma$ pocket and $g_z\left(\bm{p}\mp K\right)=\pm\beta_{\textrm{so}}$ around the $K\left(K'\right)$ pockets\cite{Wenyu}. According to the general analysis above, the three-fold rotation symmetry forces the magnetoelectric pseudotensor $T$ to be zero and there is no magnetoelectric effect. This is an interesting example of a noncentrosymmetric superconductors with zero magnetoelectric response. Interestingly, a uniaxial strain that breaks the three-fold symmetry and reduce the point group symmetry to C$_{2v}$ with the polar axis lying inside the atomic plane. As a result, SOC which is linearly proportional to $p_x$, namely, $g_z\left(\bm{p}\right)\propto p_x$, can be induced. Therefore, under uniaxial strain, the magnetoelectric pseudotensor $T$ will have nonzero element $T_{zx}$ and the $x$ component of the supercurrent will generate the out-of-plane magnetization.  This provides a novel way to use strain to generate and manipulate spin polarizations in non-gyrotropic superconductors.

\section{Discussion}
In this work, we presented the general form of the magnetoelectric pseudotensor for gyrotropic superconductors as summerized in Table \ref{magnetoelectric_pseudotensor}. Our theory provides a powerful tool for the search of unconventional magnetoelectric effects (or the Edelstein effect) in materials. Guided by the general theory, we demonstrated novel ways of generating and manipulating spin polarizations in noncentrosymmetric superconductors such as 1Td-MoTe$_2$ and 2H-NbSe$_2$. In particularly, we predict that 1Td-MoTe$_2$ is an excellent candidate for the first realization of the Edelstein effect in superconductors.

Importantly, Table \ref{magnetoelectric_pseudotensor} can be used to identify the magnetoelectric effects of a large number of superconductors with chiral lattice symmetry such as Li$_2$Pt$_3$B, Li$_2$Pd$_3$B and Mo$_3$Al$_2$C with O point group \cite{Salamon, Hirata, Prozorov}, and newly discovered superconducting chiral crystals such as TaRh$_2$B$_2$ and NbRh$_2$B$_2$ with C$_{3}$ point group~\cite{Carnicom}. The supercurrent driven longitudinal magnetization in those chiral crystals resembles the current-induced magnetization of solenoids at microscopic scale. These superconducting chiral crystals can have potential applications for new designs superconducting spintronic devices~\cite{Linder}.

It is also important to note that the Edelstein effect (current-induced spin magnetization effect) had been observed in several materials with Rashba SOC in the normal state~\cite{Gambardella}. However, Edelstein effect has not been observed in superconducting materials. In this work, we identified several superconducting materials, such as 1Td-MoTe$_2$ and strained NbSe$_2$, which possess strong Edelstein effect and the proposed Edelstein effect can be observed through magneto-optical Kerr effect measurements~\cite{Fai}.

\section*{Acknowledgement}
W.-Y. He and K. T. Law are thankful for the support of HKRGC through C6026-16W, 16324216, 16307117 and 16309718. K. T. Law is further supported by the Croucher Foundation and the Dr. Tai-chin Lo Foundation.