Kramers Weyl Semimetals as Quantum Solenoids and Their Applications in Spin-Orbit Torque Devices

TL;DR

Kramers Weyl semimetals exhibit a unique longitudinal magnetoelectric effect enabling their use as quantum solenoids and in spin-orbit torque devices with electric control of magnetization.

Contribution

This work reveals the chiral lattice symmetry-induced longitudinal magnetoelectric effect in Kramers Weyl semimetals and proposes their application in nanoscale spintronic devices.

Findings

Kramers Weyl semimetals act as quantum solenoids with parallel charge current and magnetization.

They enable electric control of magnetization switching in spintronic devices.

The materials exhibit a unique longitudinal magnetoelectric effect due to chiral symmetry.

Abstract

Kramers Weyl semimetals are Weyl semimetals that have Weyl points pinned at the time reversal invariant momenta. Recently it has been discovered that all chiral crystals host Weyl points at time reversal invariant momenta, so metals with chiral lattice symmetry all belong to the category of Kramers Weyl semimetals. In this work, we show that due to the chiral lattice symmetry, Kramers Weyl semimetals have the unique longitudinal magnetoelectric effect in which the charge current induced spin and orbital magnetization is parallel to the direction of the current. This feature allows Kramers Weyl semimetals to act as nanoscale quantum solenoids with both orbital and spin magnetization. As the moving electrons of Kramers Weyl semimetal can generate longitudinal magnetization, Kramers Weyl semimetals can be used for new designs of spin-orbit torque devices with all electric control of magnetization switching for magnets with perpendicular magnetic anisotropy.

Figures

Tables

Table 1. Table 1: of Magnetoelectric susceptibility pseudotensor α 𝛼 \alpha for the chiral crystals in the 11 chiral point groups. α i ​ j subscript 𝛼 𝑖 𝑗 \alpha_{ij} with i , j = x , y , z formulae-sequence 𝑖 𝑗 𝑥 𝑦 𝑧 i,j=x,y,z are in general the elements in α 𝛼 \alpha . In point group with symmetry, the 9 elements in α i ​ j subscript 𝛼 𝑖 𝑗 \alpha_{ij} is no longer independent. α 0 subscript 𝛼 0 \alpha_{0} means α 0 = α x ​ x = α y ​ y = α z ​ z subscript 𝛼 0 subscript 𝛼 𝑥 𝑥 subscript 𝛼 𝑦 𝑦 subscript 𝛼 𝑧 𝑧 \alpha_{0}=\alpha_{xx}=\alpha_{yy}=\alpha_{zz} in T and O point group. In the {C 3 , C 4 , C 6 , D 3 , D 4 , D 6 } group α x ​ x = α y ​ y subscript 𝛼 𝑥 𝑥 subscript 𝛼 𝑦 𝑦 \alpha_{xx}=\alpha_{yy} is then denoted as α ∥ = α x ​ x = α y ​ y subscript 𝛼 parallel-to subscript 𝛼 𝑥 𝑥 subscript 𝛼 𝑦 𝑦 \alpha_{\parallel}=\alpha_{xx}=\alpha_{yy} . α − superscript 𝛼 \alpha^{-} means the antisymmetric elements as α − = − α x ​ y = α y ​ x superscript 𝛼 subscript 𝛼 𝑥 𝑦 subscript 𝛼 𝑦 𝑥 \alpha^{-}=-\alpha_{xy}=\alpha_{yx} in group {C 3 , C 4 , C 6 }. The principal axis of the crystal is set along z 𝑧 z .

Point group	$\alpha$	Point group	$\alpha$
O		T
D2		D3
D4		D6
C1		C2
C3		C4
C6

Table 2. Table S1: List of Magnetoelectric suceptibility pseudotensor α i ​ j subscript 𝛼 𝑖 𝑗 \alpha_{ij} and the SOC pseudotensor v i ​ j subscript 𝑣 𝑖 𝑗 v_{ij} for crystals in the 18 gyrotropic point groups. Under the crystal symmetry transformation, the α i ​ j subscript 𝛼 𝑖 𝑗 \alpha_{ij} and v i ​ j subscript 𝑣 𝑖 𝑗 v_{ij} have the same behavior. The elements in α i ​ j subscript 𝛼 𝑖 𝑗 \alpha_{ij} and v i ​ j subscript 𝑣 𝑖 𝑗 v_{ij} is no longer independent. In T and O point group, α i ​ j = α 0 ​ δ i ​ j subscript 𝛼 𝑖 𝑗 subscript 𝛼 0 subscript 𝛿 𝑖 𝑗 \alpha_{ij}=\alpha_{0}\delta_{ij} , v i ​ j = v ​ δ i ​ j subscript 𝑣 𝑖 𝑗 𝑣 subscript 𝛿 𝑖 𝑗 v_{ij}=v\delta_{ij} . In the group { C 3 , C 4 , C 6 , D 3 , D 4 , D 6 } subscript C 3 subscript C 4 subscript C 6 subscript D 3 subscript D 4 subscript D 6 \left\{\textrm{C}_{3},\textrm{C}_{4},\textrm{C}_{6},\textrm{D}_{3},\textrm{D}_{4},\textrm{D}_{6}\right\} , α ∥ subscript 𝛼 parallel-to \alpha_{\parallel} and v ∥ subscript 𝑣 parallel-to v_{\parallel} is used to denote α x ​ x = α y ​ y = α ∥ subscript 𝛼 𝑥 𝑥 subscript 𝛼 𝑦 𝑦 subscript 𝛼 parallel-to \alpha_{xx}=\alpha_{yy}=\alpha_{\parallel} , v x ​ x = v y ​ y = v ∥ subscript 𝑣 𝑥 𝑥 subscript 𝑣 𝑦 𝑦 subscript 𝑣 parallel-to v_{xx}=v_{yy}=v_{\parallel} respectively. The superscript d and - are used to denote the symmetric and anti-symmetric element respectively. The coordinate is set to have principal axis along z 𝑧 z . δ i ​ j subscript 𝛿 𝑖 𝑗 \delta_{ij} is the Kronecker delta function with i , j = x , y , z formulae-sequence 𝑖 𝑗 𝑥 𝑦 𝑧 i,j=x,y,z .

Point group	$v_{ij}$	${\alpha }_{ij}$	Point group	$v_{ij}$	${\alpha }_{ij}$
C1			C6v
C2			D2d
C3			S4
C4			D2
C6			D3
C1v			D4
C2v			D6
C3v			T
C4v			O

Equations

$$H_0\left(\mathbf{k}\right)=\sum _{i,j}\frac{{\hslash }^2}{2m_{ij}}{k}_i{k}_j+{\sigma }_i\hslash v_{ij}{k}_j.$$

$$H_0\left(\mathbf{k}\right)=\frac{{\hslash }^2}{2m}{\mathbf{k}}^2+\hslash v\mathbf{k}\cdot \mathbf{\sigma },$$

$$H_0\left(\mathbf{k}\right)=\sum _{i=x,y,z}\frac{{\hslash }^2}{2m_i}{k}_i^2+\hslash v_i{\sigma }_i{k}_i.$$

$$M_i=\sum _{i,j}{\alpha }_{ij}{E}_j,$$

$$\mathcal{H}=\sum _{\nu ,{\nu }^{\prime },\mathbf{k}}{c}_{\nu ,\mathbf{k}}^{†}{H}_{0,\nu {\nu }^{\prime }}\left(\mathbf{k}\right)c_{{\nu }^{\prime },\mathbf{k}},$$

$${\alpha }_{ij}=-\tau \frac{e}{\hslash }\frac{1}{{\left(2\pi \right)}^d}{\int }_{\text{BZ}}d\mathbf{k}\sum _n{M}_{n\mathbf{k},i}{v}_{n\mathbf{k},j}\frac{df\left(E_{n\mathbf{k}}\right)}{d{E}_{n\mathbf{k}}}.$$

$${\mathbf{S}}_{{\mathbf{k}}_{\text{F}}\pm }=\pm \frac{g}{2}\frac{e\hslash }{2m_{\text{e}}}{\hat{\mathbf{k}}}_{\text{F}\pm },\text{and}{\mathbf{m}}_{{\mathbf{k}}_{\text{F}}\pm }=\frac{ev}{2}\frac{{\hat{\mathbf{k}}}_{\text{F}\pm }}{|{\mathbf{k}}_{\text{F}\pm }|}.$$

$${\alpha }_0=-\frac{e^{2}v\tau }{6{\hslash }^2{\pi }^2}\sqrt{2m{E}_{\text{F}}+m^2{v}^2}\left(g\frac{m}{m_{\text{e}}}-1\right).$$

$$\alpha =\det \left(\hat{R}\right)\hat{R}\alpha {\hat{R}}^{\text{T}}.$$

$${\mathbf{S}}_{{\mathbf{k}}_{\text{F}\pm }}=\pm \frac{g}{2}\frac{e\hslash }{2m_{\text{e}}}{\hat{\mathbf{k}}}_{\text{F}},{\mathbf{m}}_{{\mathbf{k}}_{\text{F}}\pm }=\pm \frac{1}{2}ev\frac{{\hat{\mathbf{k}}}_{\text{F}}}{|{\mathbf{k}}_{\text{F}}|},$$

$$N_{\pm }\left(E_{\text{F}}\right)=\frac{m}{2{\pi }^2{\hslash }^3}\left(\sqrt{2m{E}_{\text{F}}+m^2{v}^2}+\frac{m{v}^2}{\sqrt{\frac{2E_{\text{F}}}{m}+v^2}}\mp 2mv\right),$$

$${\mathbf{v}}_{{\mathbf{k}}_{\text{F}\pm }}=\sqrt{\frac{2E_{\text{F}}}{m}+v^2}{\hat{\mathbf{k}}}_{\text{F}\pm }.$$

$${\alpha }_0=-\frac{e^{2}v\tau }{6{\hslash }^2{\pi }^2}\sqrt{2m{E}_{\text{F}}+m^2{v}^2}\left(g\frac{m}{m_{\text{e}}}-1\right).$$

$${\mathcal{H}}_{\text{imp}}=\int d\mathbf{k}d{\mathbf{k}}^{\prime }{\psi }_{{\mathbf{k}}^{\prime }}^{†}{V}_{\mathbf{k}-{\mathbf{k}}^{\prime }}{\psi }_{\mathbf{k}},$$

$${\hat{V}}_{\mathbf{k}-{\mathbf{k}}^{\prime }}=V_{\mathbf{k}-{\mathbf{k}}^{\prime }}{U}_{{\mathbf{k}}^{\prime }}^{†}{U}_{\mathbf{k}},$$

$$U_{\mathbf{k}}=\left(\begin{array}{cc}\cos \frac{\theta }{2}{e}^{-i\varphi }& \sin \frac{\theta }{2}{e}^{-i\varphi }\\ \sin \frac{\theta }{2}& -\cos \frac{\theta }{2}\end{array}\right),$$

$${\hat{V}}_{\mathbf{k}-{\mathbf{k}}^{\prime },+}=$$

$${\hat{V}}_{\mathbf{k}-{\mathbf{k}}^{\prime },-}=$$

$$\frac{1}{{\tau }_{\pm }}=\frac{m{k}_{\text{F}\pm }{n}_{\text{imp}}}{4{\pi }^2{\hslash }^3}{\int }_0^{\pi }d\vartheta {\int }_0^{2\pi }d\phi {\hat{V}}_{\mathbf{k}-{\mathbf{k}}^{\prime },\pm }^2\left(1-\cos \vartheta \right)\sin \vartheta ,$$

$${\int }_0^{\pi }d\vartheta {\int }_0^{2\pi }d\phi {\hat{V}}_{\mathbf{k}-{\mathbf{k}}^{\prime },\pm }^2\left(1-\cos \vartheta \right)\sin \vartheta =$$

$${\int }_0^{\pi }d\vartheta {\int }_0^{2\pi }d\phi V_{\mathbf{k}-{\mathbf{k}}^{\prime },\pm }^2\left(1-\cos \vartheta \right)\sin \vartheta =$$

Full text

I Introduction

Weyl semimetals are nodal topological materials characterised by isolated band touching points, called Weyl points, in 3D momentum space Xiangang ; Balents1 ; Balents2 ; Weng . Due to the Weyl points, which act as the monopoles of the Berry curvature in the momentum space Xiangang ; Balents1 ; Balents2 ; Weng , Weyl semimetals exhibit many exotic properties such as the chiral magnetic effect Ninomiya ; Yamamoto ; Genfu ; Chenglong , the presence of topologically protected Fermi arcs states Xiangang ; Balents1 ; Balents2 ; Weng ; Suyang ; Lv , unconventional quantum oscillations Potter , and novel optical phenomenon Chris1 ; Chris2 ; Moore0 . Recently, a new type of Weyl semimetals called Kramers Weyl semimetals (KWS) in chiral crystals have been discovered Bradlyn ; Shoucheng ; Guoqing2 ; Guoqing ; Hasan ; DingHong ; Yulin ; Ando ; Felser1 ; Felser2 . Chiral crystals are crystals which lack inversion, mirror and improper rotation symmetries. As a result, a chiral crystal has a definite handedness and can be described by the 11 chiral point groups. It was shown that generally, due to the low lattice symmetry, band splittings appear away from time-reversal invariant points in momentum space and result in Kramers Weyl points pinned at time-reversal invariant momenta. However, it is not clear how the properties of KWS are distinct from that of other Weyl semimetals with non-chiral point group symmetry.

In this work, we point out that the chiral lattice symmetry in KWS brings about a unique property: an electric field applied along the principal symmetry axis of the crystal would induce spin and orbital magnetization which is parallel to the applied electric field. This is in sharp contrast to the case of all other noncentrosymmetric Weyl semimetals which give zero magnetoelectric response if the electric field is applied along the principal symmetry axis. This distinctive longitudinal magnetoelectric response involves both the spin Edelstein and orbital magnetization Moore ; Jing , and it arises from the special form of the spin-orbit coupling (SOC) of KWS. Importantly, for the representative KWS K2Sn2O3 and RhSn in the T point group we considered, the induced magnetization at a given electric field can be two to three orders of magnitude larger than the magnetization induced in materials with the strong Rashba spin-orbit couplings such as in Au (111) surfaces and Bi/Ag bilayers Mertig1 ; Mertig2 . Therefore, if the electrons which carry both orbital and spin angular momentum in the KWS are injected into a ferromagnetic layer, the torque induced by the electrons can cause magnetization switching in the ferromagnetic layer. As a result, we propose a KWS based magnetization switching device that is different from the magnetic tunneling device based on spin transfer torque Manchon ; Sinova and the conventional spin orbit torque devices Manchon ; Sinova ; Ralph ; Jungwirth . The new KWS based devices allow all electric control of magnetic switching for ferromagnets with perpendicular magnetic anisotropy, which is important for high density magnetic memories.

II Results

II.1 Effective Hamiltonians for Kramers Weyl Semimetals

In chiral crystals which respect time reversal symmetry, the energy bands are at least doubly degenerate at time-reversal invariant momenta due to the Kramers theorem. In the absence of inversion, mirror and improper rotation symmetries in chiral crystals, and away from the time-reversal invariant points, SOC would lift the Kramers degeneracy in momentum space to create Kramers Weyl points Guoqing . To be more specific, in the spin $\frac{1}{2}$ basis $\psi_{\bm{k}}=\left[\phi_{\bm{k},\uparrow},\phi_{\bm{k},\downarrow}\right]^{\textrm{T}}$, which satisfies the relation $\Theta\psi_{\bm{k}}=i\sigma_{y}\psi_{-\bm{k}}$ under time-reversal operation $\Theta=i\sigma_{y}K$, the effective Hamiltonian $H_{0}\left(\bm{k}\right)$ can be obtained through standard $\bm{k}\cdot\bm{p}$ method Guoqing . Up to second order near a time-reversal invariant momentum $\bm{k}_{0}$, the general form of the Kramers Weyl Hamiltonian can be written as

[TABLE]

Here, $\bm{k}$ is measured from $\bm{k}_{0}$, $i,j=x,y,z$, $m_{ij}$ is the effective mass tensor, $\sigma_{i}$ are the Pauli matrices in spin space, and $v_{ij}$ is the SOC pseudotensor. In chiral crystals, the little group at $\bm{k}_{0}$ is isomorphic to a chiral point group which guarantees $\det\left(v\right)\neq 0$ so that the Kramers Weyl point emerges at $\bm{k}=0$.

In the KWS, the specific forms of the SOC that creates the Weyl point are determined by the point group symmetry as $v=\det\left(\hat{R}\right)\hat{R}v\hat{R}^{\textrm{T}}$, with $\hat{R}$ the symmetry transformation matrix. In materials within the cubic point group $\left\{\textrm{T},\textrm{O}\right\}$ such as K2Sn2O3 Guoqing , RhSi Hasan , CoSi Hasan ; DingHong ; Ando , AlPt Yulin , PtGa Felser1 , and PdGa Felser2 of cubic chiral B20 structure Burkov , at $\bm{k}_{0}$ (such as the $\Gamma$ point) the little group isomorphic to $\left\{\textrm{T},\textrm{O}\right\}$ can give rise to the isotropic Weyl Hamiltonian

[TABLE]

where the high symmetry cubic point group $\left\{\textrm{T},\textrm{O}\right\}$ forces $v_{ij}$ to be proportional to the identity matrix. In those materials within the point group $\left\{\textrm{T},\textrm{O}\right\}$, the high symmetry may also enable extra multi-fold point degeneracies to appear at $\bm{k}_{0}$ Bradlyn ; Shoucheng ; Guoqing2 . In the materials CsCuBr3 Guoqing , elemental Te, Se Miyake ; Itou ; Takeshi ; Tetsuaki , etc., the dihedral point group there has lower crystal symmetry and allows the anisotropy to show up in the Weyl Hamiltonian

[TABLE]

In the materials belonging to the cyclic point group, as the crystal symmetry is further reduced, the constraints on the SOC is further reduced. The complete forms of the SOC pseudotensor $v_{ij}$ in the KWS Hamiltonian are summarised in the Supplementary Table 1 for all the chiral point groups.

The SOC in the KWS creates the Kramers Weyl points at the time reversal invariant momenta and allows the coupling between the spin and momentum. Under an electric field, the SOC enables the charge carriers to have net magnetization and such magnetoelectric response respects the same crystal symmetry present in the SOC. As shown below, for materials with cubic point groups $\left\{\textrm{T},\textrm{O}\right\}$, the simple form of the isotropic Weyl Hamiltonian in Eq. 2 allows us to calculate the magnetoelectric susceptibility analytically. For the KWS in dihedral and cyclic point group, the magnetoelectric responses are calculated numerically for the selected materials in Supplementary Note 1 and Supplementary Note 2 respectively.

II.2 Magnetoelectric Pseudotensors and Their Symmetry Properties

In magnetoelectric effects, induced magnetization $\bm{M}$ and the applied electric field $\bm{E}$ are related by the magnetoelectric pseudotensor $\alpha$ such that:

[TABLE]

where $i,j=x,y,z$ and $\alpha_{ij}$ are elements of the magnetoelectric pseudotensor $\alpha$. For a generic Hamiltonian

[TABLE]

where $c^{\dagger}_{\nu,\bm{k}}\left(c_{\nu,\bm{k}}\right)$ is the creation (annihilation) operator, $H_{0,\nu\nu^{\prime}}\left(\bm{k}\right)$ is the element of the Hamiltonian matrix $H_{0}\left(\bm{k}\right)$, $\alpha_{ij}$ can be obtained from the linear response theory as Moore ; Jing

[TABLE]

In Eq.6, $f\left(E_{n\bm{k}}\right)$ is the Fermi Dirac distribution function, $E_{n\bm{k}}$ is the energy dispersion of band $n$ from the Hamiltonian $H_{0}\left(\bm{k}\right)$, $v_{n\bm{k},j}=\frac{\partial E_{n\bm{k}}}{\partial k_{j}}$, $d$ is the dimension of the system, $\tau$ is the effective scattering time and $i,j=x,y,z$ denote the spatial components. The total magnetic moment $\bm{M}_{n\bm{k}}=\bm{S}_{n\bm{k}}+\bm{m}_{n\bm{k}}$ carried by the Bloch electrons consists of both the spin magnetic moment $\bm{S}_{n\bm{k}}=\bra{\phi_{\bm{k},n}}\frac{1}{2}g\mu_{\textrm{b}}\bm{\sigma}\ket{\phi_{\bm{k},n}}$ and the orbital magnetic moment $\bm{m}_{n\bm{k}}=\frac{ie}{2\hbar}\bra{\partial_{\bm{k}}\phi_{\bm{k},n}}\times\left[H_{0}\left(\bm{k}\right)-E_{n\bm{k}}\right]\ket{\partial_{\bm{k}}\phi_{\bm{k},n}}$. Here, $\mu_{\textrm{b}}=\frac{e\hbar}{2m_{\textrm{e}}}$ is the Bohr magneton, $g$ is the Lande $g$ factor which is set to be 2 in our calculations and $\ket{\phi_{\bm{k},n}}$ denotes a Bloch state. As we will show explicitly below, the orbital magnetization is related to the Berry curvature of the Bloch states which has the form $\bm{\Omega}_{n\bm{k}}=i\bra{\partial_{\bm{k}}\phi_{\bm{k},n}}\times\ket{\partial_{\bm{k}}\phi_{\bm{k},n}}$  Niuqian .

The linear response theory applies to generic Hamiltonians. However, to shed light on the general properties of KWS, we note that the form of $\alpha$ can be determined by point group symmetries which is independent of the details of the Hamiltonian. The group theory analysis of $\alpha$ is elaborated in the Method Section as well as in Ref. Wenyu2 , and the general form for the chiral point groups is provided in Table 1. From the group theory point of view, KWS can be classified into three sub-classes. For KWS belonging to the cubic point groups $\left\{\textrm{T},\textrm{O}\right\}$ , $\alpha$ is proportional to the identity matrix as shown in Table 1. This implies that the induced magnetization is always parallel to the direction of the applied electric field. Therefore, these KWS can behave as classical solenoids in all electric field directions without the need to fabricate any spiral structures.

For KWS with point groups Dn, a pure longitudinal magnetization parallel to the electric field is also obtained when the electric field is applied along the direction of any of the symmetry axes. For KWS with cyclic point groups, in general, magnetization with both components parallel and perpendicular to the direction of the applied electric field are generated. Interestingly, for all other Weyl semimetals without chiral point group symmetry, the magnetoelectric response is zero if the electric field is applied along the principal axis, which can be obtained from the Supplementary Table 1. Therefore, the longitudinal magnetoelectric response along the principal symmetry axis is a very distinctive feature of KWS due to the special spin texture of KWS which determines the spin and orbital magnetization.

Combined with the effective Hamiltonian for isotropic KWS in Eq. 2, the magnetoelectric susceptibility in Eq. 6 can be explicitly calculated and it shows how the SOC strength and band dispersion will influence the magnetoelectric response.

II.3 Longitudinal Magnetoelectric Response in KWS

We first consider an effective Hamiltonian which describes a Kramers Weyl point near the $\Gamma$ point in chiral crystals with point groups $\left\{\textrm{T},\textrm{O}\right\}$, where the isotropic Weyl Hamiltonian is $H_{0}\left(\bm{k}\right)=\frac{\hbar^{2}}{2m}\bm{k}^{2}+\hbar v\bm{k}\cdot\bm{\sigma}$. At the Fermi energy $E_{\textrm{F}}=\frac{\hbar^{2}}{2m}\bm{k}_{\textrm{F}}^{2}\pm\hbar v|\bm{k}_{\textrm{F}}|$, there are two spherical Fermi surfaces with corresponding wave vectors $\bm{k}_{\textrm{F}\pm}=k_{\textrm{F}\pm}\hat{\bm{k}}$, where $k_{\textrm{F}\pm}=\frac{1}{\hbar}\sqrt{2mE_{\textrm{F}}+m^{2}v^{2}}\mp\frac{mv}{\hbar}$ and $\hat{\bm{k}}=\frac{\bm{k}}{|\bm{k}|}$. The spin and orbital magnetic moments of the two Fermi surfaces with Fermi momenta $\bm{k}_{\textrm{F}\pm}$ can be written as

[TABLE]

It is important to note that the orbital magnetic moment $\bm{m}_{\bm{k}_{\textrm{F}}\pm}$ is proportional to the Berry curvature generated by the Weyl point inside Fermi surfaces which is $\bm{\Omega}_{\bm{k}_{\textrm{F}\pm}}=\mp\frac{\hat{\bm{k}}_{\textrm{F}\pm}}{2\bm{k}^{2}_{\textrm{F}}\pm}$. The spin texture on a Fermi surface is schematically shown in Fig.1a. It is clear that without breaking time-reversal symmetry, the total magnetic moment of all the electrons is zero. By applying an electric field, the steady state distribution of the electronic state can generate a net magnetization as indicated in Eq. 6. With this special form of spin texture of an isotropic KWS, at the Fermi energy $E_{\textrm{F}}$, we obtain the isotropic longitudinal magnetoelectric susceptibility $\alpha_{0}$ as (see the method section for detailed derivation)

[TABLE]

In $\alpha_{0}$, the first and second terms are the spin and the orbital contributions respectively. It is important to note that, for hole bands with negative effective mass, the spin and orbital contribution will always add together to enhance the magnetoelectric response. From Table 1, the magnetoelectric response of an isotropic KWS can be written as $\bm{M}=\alpha_{0}\bm{E}$ which indicates that the magnetization induced is parallel to the applied electric field as is schematically shown in Fig.1b. From Eq. 8, it is clear that strong Weyl SOC $v$, long scattering time $\tau$ and large effective mass $m$ can give large magnetoelectric response.

To acquire a strong magnetoelectric response, it is also preferable to have long scattering time as shown in Eq.8. In KWS, since the Kramers Weyl points are pinned at the time reversal invariant momenta, the electrons on opposite sides of the Weyl point have opposite spin. As a result, for elastic back scattering from scalar impurities, the intra Weyl pocket scattering is suppressed by the SOC, similar to the case in the surface states of topological insulators Yazdani ; Xue ; Shoucheng2 . For an ideal KWS that only has Fermi surfaces enclosing the Weyl point at time reversal invariant momenta, the Weyl SOC can enhance the scattering time by a factor of 3 as shown in the method section using Born approximation. Therefore, we choose a scattering time $\tau=1$ps, which is in the range of the inter Weyl point scattering time in TaAs Mertig2 ; Ramshaw ; ShuangJia . With the effective mass $m=1.4m_{e}$, Weyl SOC strength $\hbar v=-70$meV$\cdot$nm, the KWS have the band dispersion shown in Fig. 1c. The isotropic magnetoelectric susceptibility $\alpha_{0}$ as a function of Fermi energy $E_{\textrm{F}}$ is then evaluated as shown in Fig. 1d. For the Fermi pockets enclosing a single Kramers Weyl point, the isotropic magnetoelectric susceptibility $\alpha_{0}$ increases with the square root of the Fermi energy $E_{\textrm{F}}$.

To seek a large magnetoelectric response from realistic materials, we study two representative materials in the T point group: the K2Sn2O3 Guoqing , the RhSn Burkov . The K2Sn2O3 is a small gap insulator with a Kramers Weyl point at H near the conduction band bottom. Around the Kramers Weyl point, the conduction bands of K2Sn2O3 have large spin splitting shown in Fig.2a. In the slight n-doped state, the electrons will occupy the Fermi pockets enclosing the single Kramers Weyl point at H so that the isotropic Weyl Hamiltonian can be effectively described by Eq.2, with the effective mass $m=1.4m_{\textrm{e}}$ and SOC strength $\hbar v=-70$meV$\cdot$nm. To validate the calculations using an effective Hamiltonian in the form of Eq.2, a realistic tight-binding model for K2Sn2O3 is further constructed (see the method section) to calculate the magnetoelectric susceptibility $\alpha_{0}$, as is shown in Fig. 2b. Below $E=0.7$eV, only the Weyl bands enclosing the Kramers Weyl point at H are involved so that the numerically calculated magnetoelectric susceptibility matches well with the analytical formula in Eq.8. It is interesting to note that there is a four-fold fermion around the energy $E=0.75$ eV. As the Fermi surfaces enclosing multi-fold fermion also contribute to the magnetoelectric susceptibility $\alpha_{0}$ Felix , the four-fold fermion at H point induce a sudden increase of $\alpha_{0}$ when the chemical potential is above 0.7eV. Assuming that the chemical potential lies at $E=0.6$eV (about 80meV above the Weyl point), an electric field of $10^{5}$V/m is enough to generate a magnetization of $0.005\mu_{\textrm{b}}$ per unit cell.

On the other hand, RhSn, which has the same B20 structure as the recently studied KWS RhSi Hasan , is semimetallic. The SOC generates the Kramers Weyl points at $\Gamma$ as well as other multi-fold fermion band crossing points  Bradlyn ; Shoucheng ; Guoqing2 at the time reversal invariant momenta as shown in Fig. 3a. A four-fold fermion band crossing point near the Fermi energy is highlighted by a red circle in Fig. 3a. With the realistic tight binding model constructed from Wannier functions (see the method section), the magnetoelectric susceptibility $\alpha_{0}$ is calculated in the energy range including all the bands. The results are shown in Fig. 3b. At the Fermi energy, which is near the multi-fold fermion Weyl point, an electric field of $10^{5}$V/m can generate a magnetization of 0.007$\mu_{\textrm{b}}$ per unit cell . Compared with the magnetoelectric response of Au(111) surface and Bi/Ag bilayers which have large Rashba SOC Mertig1 ; Mertig2 , the K2Sn2O3 and RhSn can generate magnetizations which are two to three orders of magnitude larger.

For KWS with dihedral (Dn) point groups, such as CsCuBr3 Guoqing , and elemental Te, Se Miyake ; Itou ; Takeshi ; Tetsuaki , a pure longitudinal magnetoelectric response can also be obtained when the electric field is applied along the direction of the symmetry axes. On the other hand, for KWS with cyclic (Cn) point groups, such as Ca2B5Os3, a longitudinal magnetoelectric response is generally accompanied with a transverse response. The calculations of magnetoelectric tensors for CsCuBr3 with D2 point group and Ca2B5Os3 with C2 point group can be found in the Supplementary Note 1 and Supplementary Note 2 respectively.

III Discussion

We note that the current induced magnetization in Weyl semimetals was first studied by Johansson et al. in TaAs which belong to the point group C4v with mirror planes Mertig2 . Therefore, the Weyl points appear at general $\bm{k}$ points and the resulting spin polarization induced by an electric field is perpendicular to the direction of the electric field. This transverse magnetoelectric effect is similar to the case with Rashba SOC Edelstein which is a property of the polar point group as shown in Supplementary Table 1. If the applied electric field is along the principal axis, however, the magnetoelectric response is zero in TaAs as determined by the C4v point group symmetry.

On the other hand, the longitudinal magnetoelectric response in Weyl semimetals was first studied by Yoda et al. Murakami ; Yoda . However, in their models, helical hopping textures Murakami ; Yoda are needed for electrons to hop in a spiral manner, imitating the movement of electrons in a classical solenoid. As a result, an orbital magnetization parallel to the direction of an applied electric field would be generated and the longitudinal magnetization is present even without spin-orbit coupling. Unfortunately, no realistic materials that possess such helical hopping textures are identified.

In this work, the finite longitudinal magnetoelecrtric response comes from the current induced net magnetic moments accumulation on Fermi surfaces of the KWS. The net magnetic moments on Fermi surfaces basically have two origins: 1) the SOC induced spin splitting; 2) the finite orbital magnetic moment distribution on Fermi surfaces. The two-fold Kramers Weyl point we focused on in this work contributes to both the spin and orbital magnetic moments, but for a general KWS of chiral crystal symmetry, the complicated band structure diversifies the origin of magnetic moments on Fermi surfaces.

As we have seen, the Fermi pockets enclosing multi-fold fermion contribute to the magnetoelectric susceptibility as well. It is due to the fact that multi-fold fermion endows orbital magnetic moments to the Fermi surfaces outside. Along with the SOC induced spin splitting, the Fermi pockets enclosing multi-fold fermion can also give rise to spin and orbital magnetization driven by charge current Felix . In a KWS, the magnetoelectric response from the multi-fold fermion and two-fold Kramers Weyl point is the same type as the response tensor form is determined by the chiral crystal symmetry of the KWS. Since the Fermi surfaces with a multi-fold fermion inside have different magnetic moment, group velocity and density of states from that with a two-fold Kramers Weyl point, the magnetoelectric susceptibility can differ in the value for the two cases.

Besides the two-fold Kramers Weyl point and multi-fold fermion at high symmetry points in the Brillouin zone, a Weyl point at a general $\bm{k}$ point can generate finite longitudinal magnetoelectric susceptibility as long as the crystal symmetry does not reduce it to zero. This is the case of the local maximal magnetoelectric susceptibility $\alpha_{0}$ of RhSn at $E=-0.46$eV shown in Fig. 3 b. There is a Weyl point in $\Gamma-R$ line around the energy of $E=-0.46$eV which enhances the magnetoelectric susceptibility at that energy. In a brief summary, the longitudinal magnetoelectric susceptibility in KWS of chiral crystal symmetry originates from SOC and Weyl points in the band structure. The magnitude of longitudinal magnetoelectric response at given current density relies on the Fermi surface magnetic moment, Fermi velocity, effective scattering time and density of states at the Fermi energy of KWS.

Concerning the applications of KWS, due to its unique longitudinal magnetoelectric response, the KWS/ferromagnet heterostructures can be used for new designs of spin-orbit torque devices as shown in Fig.4a. It is interesting to note that a KWS can cause current induced magnetization switching in the ferromagnetic layer by two effects. First, as depicted in Fig.4a, the KWS can inject electrons which carry both orbital and spin angular momentum into the ferromagnetic layer. The injection of angular momentum can cause magnetization switching similar to the spin transfer torque induced magnetization switching in magnetic tunnelling junctions as shown in Fig.4b. The important difference is that the electrons injected by KWS can carry both orbital and spin angular momentum. Particularly, close to the Weyl points, the orbital magnetization carried by the electrons can be significant. Therefore, KWS can work as a source of spin and orbital angular momentum for magnetization switching. Second, when a current is passed along the principal symmetry axis of the KWS, a magnetization is induced at interface between the KWS and the ferromagnet. This current induced magnetization can cause the magnetization switching of the ferromagnet layer through ferromagnetic coupling between the KWS and the ferromagnet.

In recent years, tremendous progress has been made for the study of spin-orbit torque devices as depicted in Fig.4c. For example, current induced magnetization switching due to spin-orbit torques, in the presence of an in-plane magnetic field, has been realized experimentally in heavy metal/ferromagnet heterostructures Gambardella ; Mingzhong . It was also demonstrated that charge currents in multi-layer WTe2 can induce out-of-plane magnetizations Ralph and WTe2 can also be used for magnetization switching for ferromagnets with in-plane anisotropy Hyunsoo . However, an all electric control of magnetization switching for ferromagnets with perpendicular magnetic anisotropy (PMA) through spin-orbit torques, which is important for high density magnetic memories, has not been experimentally realized. Therefore, due to the unique longitudinal response of KWS, they allow new designs of spintronic devices for magnetization switching of ferromagnets with perpendicular magnetic anisotropy as depicted in Fig.4a.

Concerning the candidate materials of KWS, the cubic chiral B20 structured materials Burkov , such as the recent experimentally studied RhSi Hasan , CoSi Hasan ; DingHong ; Ando , AlPt Yulin , PtGa Felser1 , PdGa Felser2 , and RhGe, RhSn all belong to the T point group that has the isotropic longitudinal magnetoelectric response. Particularly, the PtGa Felser1 that has the strongest SOC reported so far in chiral crystals is expected to produce strong longitudinal magnetoelectric effect. Besides, elemental Se and Te also have chiral lattice structure Miyake ; Itou . Interestingly, a few superconducting materials with chiral lattice symmetry and strong spin-orbit coupling such as Li2Pt3B Salamon , Li2Pd3B Hirata , Mo3Al2C Prozorov , TaRh2B2 and NbRh2B2 Cava have been experimentally studied. In their normal state, these superconducting materials can also have spintronic applications.

In conclusion, the magnetoelectric effect of KWS are studied in this work using three different methods: 1. Group theory analysis; 2. Analytical calculations through linear response theory; and 3. Realistic tight-binding calculations are constructed through first principle calculations and the magnetoelectric response of realistic materials are calculated. A new design for current-controlled magnet memory device is proposed and new candidates of KWS are suggested.

IV Methods

IV.1 Symmetry Analysis for the Magnetoelectric Response in Chiral Crystals

In the chiral crystals with magnetoelectric effect $M_{i}=\sum_{i,j}\alpha_{ij}E_{j}$ with $i,j=x,y,z$, under the crystal symmetry the magnetization transforms as $M_{i}\rightarrow\det\left(\hat{R}\right)\hat{R}_{ij}M_{j}$ while the electric field transforms as $E_{i}\rightarrow\hat{R}_{ij}E_{j}$, where $\hat{R}$ represents the symmetry transformation operator and is an orthogonal matrix. As a result, the magnetoelectric susceptibility $\alpha_{ij}$ under the crystal symmetry respects

[TABLE]

The chiral point groups, which do not allow improper rotations, can be divided into three sub-classes: the cubic point groups $\left\{\textrm{T},\textrm{O}\right\}$, the dihedral point groups $\left\{\textrm{D}_{n}\right\}$ with $n=2,3,4,6$ and the cyclic point groups Cn with $n=1,2,3,4,6$. In the cubic point groups $\left\{\textrm{T},\textrm{O}\right\}$, the multiple high order rotation axes along different directions would force the magnetoelectric susceptibility to be proportional to the identity matrix. In the dihedral point groups $\left\{\textrm{D}_{n}\right\}$ with $n=2,3,4,6$, the Cn rotation axis along $z$ and the in-plane C2 rotation axis along the $x$-axis would eliminate all the off diagonal elements and leave only the diagonal elements in the magnetoelectric susceptibility $\alpha=\textrm{diag}\left\{\alpha_{xx},\alpha_{yy},\alpha_{zz}\right\}$. In $\left\{\textrm{D}_{3},\textrm{D}_{4},\textrm{D}_{6}\right\}$, the principal axis would further make $\alpha_{xx}=\alpha_{yy}$. In the cyclic point groups $\left\{\textrm{C}_{n}\right\}$ with $n=1,2,3,4,6$, the lower symmetry would allow the off-diagonal elements to coexist with the diagonal elements. The explicit forms of the magnetoelectric susceptibility pseudotensor $\alpha_{ij}$ from the symmetry analysis is shown in Table 1 for the 11 chiral point groups.

IV.2 Magnetoelectric susceptibility in isotropic KWS

In the isotropic KWS, at the Fermi energy $E_{\textrm{F}}=\frac{\hbar^{2}\bm{k}^{2}}{2m}\pm\hbar v|\bm{k}|$, there are two spherical Fermi surfaces with the corresponding wave vectors $\bm{k}_{\textrm{F}\pm}=k_{\textrm{F}\pm}\hat{\bm{k}}$ and $k_{\textrm{F}}=\frac{1}{\hbar}\left(\sqrt{2mE_{\textrm{F}}+m^{2}v^{2}}\mp mv\right)$, $\hat{\bm{k}}=\frac{\bm{k}}{|\bm{k}|}$. On the Fermi surfaces enclosing the Kramers Weyl point, the spin magnetic dipole moment and orbital magnetic dipole moment are respectively

[TABLE]

with $\pm$ denoting the two band branches. The density of states can be obtained as

[TABLE]

and the Fermi velocities there become

[TABLE]

At the Fermi energy, the magnetoelectric susceptibility $\alpha_{0}$ can be approximated through $\sum_{n}\frac{1}{\left(2\pi\right)^{d}}\int_{\textrm{BZ}}...\frac{df\left(E_{n,\bm{k}}\right)}{dE_{n,\bm{k}}}d\bm{k}\rightarrow-\sum_{n}N_{n}\left(E_{\textrm{F}}\right)\int\frac{d\Omega_{\bm{k}_{\textrm{F},n}}}{4\pi}$ with $n=\pm$ and then we can obtain

[TABLE]

IV.3 Isotropic Weyl SOC Suppressed Back Scattering

In the pure scalar potential scattering process

[TABLE]

the scattering potential can be projected onto the band basis as

[TABLE]

with

[TABLE]

and $\theta$, $\phi$ are the azimuth and polar angle of $\bm{k}$. In the presence of large SOC, the scattering is dominated in each band branch $\left(\pm\right)$ and the effective scattering potential becomes

[TABLE]

In the Born approximation, the scattering in each band branch is obtained as

[TABLE]

with $\varphi=\phi^{\prime}-\phi$, $\vartheta=\theta^{\prime}-\theta$. In the spherical Fermi surfaces, the scattering potential $\hat{V}^{2}_{\bm{k}-\bm{k}^{\prime},\pm}$ only depends on the angle between $\bm{k}$ and $\bm{k}^{\prime}$, so

[TABLE]

As a result, the Weyl SOC reduces the scattering potential to about 0.3 and suppresses the back scattering.

IV.4 First-principles Calculation and Wannier Function Construction

We perform first-principles calculation within the density functional theory framework as implemented in Vienna Abinitio Simulation Package(VASP) Kresse . PAW Blochl ; Joubert type of pseudopotential with PBE exchange functional Ernzerhof is used in the calculations and spin orbit coupling is included in the pseudopotentials. After the first-principle calculation is done, tight binding Hamiltonian which perfectly recovers bands near Fermi surface is built through Wannier function construction using package Wannier90 Vanderbilt ; Mostofi . The lattice structures of all the materials we considered are obtained from the ICSD ICSD .

Data Availability

The data generated from our codes that support the findings of this study are available from the corresponding author upon reasonable request.

Author Contributions

K. T. L. and W.-Y. H. conceived the idea and initiated the project. W.-Y. H. performed the symmetry analysis and the theoretical calculation of the magnetoelectric susceptibility. X. Y. X. performed the first principles calculation and Wannier function construction for the selected Kramers Weyl semimetals in this work. All the authors discussed the results and co-wrote the paper.

Acknowledgement

The authors thank Xi Dai, Kin Fai Mak, Qiming Shao and Binghai Yan for inspiring discussions. K. T. Law is 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.

Competing Interests

The authors declare no competing interests.

V Suppelementary Note 1: Anisotropic Longitudinal Magnetoelectric Response in KWS in Dihedral Point Groups

In the KWS with dihedral point group symmetry, the magnetoelectric susceptibility has the nonzero value only in the diagonal elements. This indicates that the magnetization generated is parallel to the electric field if the electric field is applied along the symmetry axes. For example, the KWS CsCuBr3 belongs to the D2 point group and its realistic band structure in Supplementary Figure S1a shows the Kramers Weyl points at the time reversal invariant momenta {$\Gamma$, Y, S, Z}. Near the Kramers Weyl points, the bands splittings are highly anisotropic as allowed by the D2 point group symmetry. The anisotropic longitudinal magnetoelectric susceptibility $\left\{\alpha_{xx},\alpha_{yy},\alpha_{zz}\right\}$, as a function of the Fermi energy, are numerically evaluated based on the realistic band structure, are depicted in Supplementary Figure S1b.

VI Suppelementary Note 2: Coexistence of Longitudinal and Transverse Magnetoelectric response in KWS in Cyclic Point Groups

The KWS in the cyclic point group possess a polar axis that leaves the crystal invariant under rotation. This symmetry allows both the longitudinal and the transverse magnetoelectric response: an electric field along the polar axis can generate pure longitudinal magnetization while the electric field inside the basal plane produces magnetization with nonzero parallel and perpendicular components. In Ca2B5Os3 with a single C2 rotation axis, the spin splitted bands with Kramers Weyl points are shown in Supplementary Figure S1c. Based on the realistic band structure, the nonvanishing magnetoelectric susceptibility elements under the C2 symmetry are numerically calculated at different energy. The results are shown in Supplementary Figure S1d and large magnetoelectric response is also obtained due to the strong spin-orbit coupling of the material.