Thermoelectric transport in Weyl semimetals under a uniform concentration of torsional dislocations

In this article, we present an effective continuum model for a Weyl semimetal, to calculate its thermal and thermoelectric transport coefficients in the presence of a uniform concentration of torsional dislocations. We model each dislocation as a cylindrical region of finite radius a, where the corresponding elastic strain is described as a gauge field leading to a local pseudo-magnetic field. The transport coefficients are obtained by a combination of scattering theory, Green's functions and the Kubo formulae in the linear response regime. We applied our theoretical results to predict the electrical and thermal conductivities as well as the Seebeck coefficient for several transition metal monopnictides, i.e. TaAs, TaP, NbAs and NbP.


Introduction
Not long aer being postulated as a theoretical concept, [1][2][3][4][5][6][7] Weyl semimetals (WSMs) were discovered in TaAs crystals. 8WSMs constitute important examples of three-dimensional, gapless materials with non-trivial topological properties, as their band structure displays an even number of Weyl nodes.5][6][7] In particular, each node is a monopolar source of Berry curvature, and hence they are protected from being gaped since their topological charge (chirality) is an invariant. 7This implies that in Weyl fermions, the projection of spin over their momentum direction is preserved, a condition referred to as "spinmomentum locking".Some remarkable properties related to the existence of Weyl nodes in the bulk band structure are the presence of Fermi arcs, 8 the chiral anomaly, and the chiral magnetic effect. 9][12][13][14][15][16] Different estimations in the literature report rst-principles calculations for the optical conductivities in the monopnictide family (TaAs, TaP, NbAs and NbP), [17][18][19][20][21][22] which in the lowfrequency (DC) limit are in the range s xx ∼ 10 4 to 10 6 U −1 cm −1 (see Table 1).Concerning the electronic contribution to the thermal conductivity, including the aforementioned scattering mechanisms, estimations based on rst-principles calculations [23][24][25] report values in the range k (el) xx ∼ 20-100 W K −1 m −1 .The lattice contribution, on the other hand, is strongly dependent on the masses of the nuclei, and hence it varies in a wider range for the different materials k (l) xx ∼ 1-190 W K −1 m −1 (see Table 1).It has been proposed that generic semi-metals may constitute attractive candidates for thermoelectric applications due to their relatively large Seebeck coefficients at room temperature jSj ∼ 10 2 mV K −1 . 26This parameter is very sensitive to the density of carriers through the chemical potential, but different estimations in the literature for the family of transition metal monopnictides report values in the range jSj ∼ 10 2 to 10 3 mV K −1 (ref.23-25) at room temperature.Therefore, the general concept of "Topological Thermoelectrics" has generated a lot of interest in the materials science community, with excellent recent reviews 27 on the subject.
In contrast with results reported in the literature, the effects of mechanical strain and dislocations or disclinations have been theoretically explored to a much lesser extent in the context of electronic and thermal transport properties.Those defects can be modeled in a continuum approximation by gauge elds [28][29][30] in WSMs.More recently, the role of gauge elds has been explored in acoustic crystal realizations of topological materials as well, 31 particularly in their role in representing topological defects. 32ble 1 Values of the DC conductivity s xx , the electronic k (el) xx , and lattice k (l) xx contributions to the thermal conductivities at 300 K reported in the literature In our previous studies, we have studied quasi-ballistic transport through a nano-junction in a WSM with a single torsional dislocation, in combination with an external magnetic eld.For such a system, we obtained the electronic 33,34 and thermoelectric 33,35 transport coefficients, using the Landauer ballistic formalism in combination with a mathematical analysis for the quantum mechanical scattering cross-section. 36ore recently, we considered the case of a diluted, uniform concentration of torsional dislocations and their effects on the electrical conductivity of type I WSMs, 37 by means of the Kubo linear-response formalism.The effect of the random distribution of dislocations, with a concentration n d (per unit area), is incorporated in the form of a disorder-averaged self-energy into the corresponding Dyson's equation for the retarded and advanced Green's functions.Furthermore, as described in ref.
37 a vertex correction obtained as a solution to the Bethe-Salpeter equation was incorporated into the Kubo linear response formulae.
In the present work, our purpose is to further extend this study, in combination with Onsager relations of nonequilibrium thermodynamics, to obtain the electronic component of the thermal conductivity and Seebeck coefficient in these materials, limited by this particular scattering mechanism, as a function of temperature and concentration of dislocations.We remark that this is the single scattering mechanism that we shall focus on this study, since it requires a special modeling as compared to other mechanisms that have been already discussed extensively in the literature.Moreover, as we state in the Discussion, Mathiessen's rule allows one to combine all these different contributions via the overall relaxation time in the estimation of the transport coefficients.
We present explicit evaluations of our analytical expressions for the electrical and thermal conductivity, as well as for the Seebeck coefficient, as a function of temperature and concentration of dislocations n d , for several materials in the family of transition metal monopnictides, i.e.TaAs, TaP, NbAs and NbP, with microscopic parameters estimated from ab initio calculations as reported in the literature. 17,38,39Our calculations show that, although the Wiedemann-Franz law is satised for all such compounds in the low-temperature limit, the Seebeck coefficient leads to a large gure of merit ZT (el) > 2 even at room temperature for TaAs.Therefore, our theoretical results suggest that the transition metal monopnictides may constitute very attractive candidates for thermoelectric applications in energy harvesting.Since our model only captures the electronic contribution to the thermal conductivity, this possibility must be further explored to evaluate in more detail the role of phonon-related scattering effects and lattice thermal conductivity.

Scattering by a uniform concentration of dislocations
Let us start with an effective continuum model for a type I WSM, in the presence of a uniform concentration n d = N d /A (per unit transverse surface) of identical cylindrical dislocations of nite radius a, as depicted in Fig. 1.The spatial distribution of such defects is represented by the density function where X j is the position of the j th -dislocation's axis.We model this system using the Hamiltonian 37 where for s = (s x , s y , s z ) the vector of Pauli matrices represents the free-particle Hamiltonian at each of the Weyl nodes K ± = ±b/2, labeled by their corresponding chirality index x = ±, and v F is the Fermi velocity.In addition, V represents the scattering potential due to the presence of the random distribution of dislocations, given that the contribution from a single dislocation defect is given by [34][35][36][37] Ĥx Here, r j = jx − X j j is the distance from the center of each dislocation 1 # j # N d (see Fig. 1), and f = (−sin f,cos f,0) the azimuthal unit vector in polar coordinates.Eqn (5) contains the interaction with each cylindrical dislocation of radius a, where Fig. 1 Random distribution of torsional dislocations, as seen from a plane perpendicular to the cylinder axis.Each dislocation is modeled as a cylinder of radius a, whose central axis is located at the vector X j (on the perpendicular plane x-y).The position of an arbitrary point P on the plane is x = (x, y), and we define r j = jx − X j j as its relative distance from the j th dislocation axis.
torsional strain is described as a pseudo-magnetic eld B x in its interior r < a, [34][35][36][37] as described by the Heaviside function Q(a − r).The corresponding lattice mismatch effect at the boundary r = a is represented by a repulsive delta barrier 35,37 with strength V 0 .In this formalism, the pseudo-magnetic eld B x = V × A x representing strain is proportional to the torsional angle q (in degrees), a relation that is convenient to express in terms of its ux through the circular cross-section of each cylindrical region: jB x ja 2 = 1.36qf0. 35Here, we dened the modied ux quantum representing the dislocations in these materials by f0 h ħv F e z 330 T Å 2 . 35 described in detail in ref. 37, we include the effect of disorder by taking the congurational average over the statistical distribution of dislocations, dened as where f(X j ) is any function of the dislocations' positions and P(X j ) is their statistical distribution function in the sample.In particular, for a uniform distribution we have P(X j ) = 1/A, where A is the area of the plane normal to each cylinder's axis.
As we shall present in the next section, for the calculation of the thermal and thermoelectric transport coefficients in this material, we are interested in the disorder-averaged retarded Green's function Here E ðxlÞ k ¼ xlħv F jkj is the energy spectrum of the "free" massless Weyl fermions, with l = ±1 the band index and x = ±1 their chirality.In addition, the retarded self-energy has the form As usual, the real part of the self-energy renormalizes the single-particle energy spectrum, while the imaginary part represents the scattering relaxation time s (xl) (k) through the relation The advanced self-energy is given by the complex conjugate of the retarded self-energy, i.e., S (xl) A where the elements of the T-matrix are given by T ðxlÞ k 0 k k k ¼ hJ k 0 ;l TðxlÞ ðEÞ J k;l i and k = (k ‖ , k z ), with k ‖ = (k x , k y ) the momentum on the plane perpendicular to the cylindrical dislocation axis.As we showed in ref. 37, the real part of the selfenergy is expressed in terms of the phase shis d m (k) for each angular momentum component m ˛Z This innite series over highly oscillatory terms converges to zero, and therefore no contribution arises from the real part of the self-energy.On the other hand, the imaginary part of the self-energy in eqn (9) gives the scattering relaxation time in terms of the phase shis and we can see that it is a positive denite quantity, inversely proportional to the concentration of dislocations s The phase shis d m (k) for this system were calculated in ref. 35  and their analytical expression is given in eqn (1) of the ESI.† Also, the explicit expression for the T-matrix elements in terms of these phase shis is given in eqn (2) of the ESI.†

Onsager coefficients in the linear response regime
In the present work, our purpose is to study the thermal and thermoelectric transport coefficients in these topological materials in the presence of a nite concentration of dislocations n d as the single scattering mechanism.For this purpose, in what follows we shall apply the basic principles of nonequilibrium thermodynamics.Associated with the particle current ĵ, dened by the operator and the heat current operator j (x) the macroscopic currents are given by the corresponding ensemble averages The entropy production rate is expressed in terms of the macroscopic currents and gradients as follows 41,42

Paper Nanoscale Advances
Let us introduce the Onsager coefficients by means of the tensor notation The transport coefficients can also be expressed in terms of these tensors, by applying the corresponding denition as follows: rst, let us assume that VT = 0 and Vm = 0, such that the electrical conductivity tensor is then given by On the other hand, we remark that the thermal conductivity, by denition, is measured under conditions such that no electric current ows through the material J = 0.Then, combining eqn (17a) and (17b) we conclude that the thermal conductivity tensor is given by the expression while the Seebeck coefficient (thermopower) is given by The Onsager coefficients can be expressed in terms of the Kubo formulae in the linear response regime.From the entropy production rate in eqn ( 16), we have where J 1 = J, J 2 = J Q , and In eqn (21), F(t) is usually termed the "dissipation function". 41e shall apply the Luttinger formalism 43 for the evaluation of the Onsager coefficients.For this purpose, we begin by expressing the Kubo formulae for the different currents (i = 1, 2) in the form where s is a positive quantity that guarantees the adiabatic switching-on of the perturbation that drives the system out of equilibrium, and the limit s / 0 + is taken at the end of the calculation.In eqn (23), we also dened the equilibrium density operator where is the grand-canonical partition function.
When inserting eqn (21) into eqn (23), we obtain Then, we conclude In these expressions, the electric eqn ( 13) and heat current eqn ( 14) operators, respectively, can be combined into a single denition (for i = 1, 2) Finally, as we show in Section 3 of the ESI, † we obtain the corresponding Onsager coefficients by considering the spatial average of the corresponding tensors.This is equivalent, in Fourier space, to take the limit of the momentum q / 0 in each of these coefficients where we have taken into account the vertex corrections G RA (k, E), as described in Section 2 of the ESI.† At low temperatures, a closed analytical solution is possible since the derivative of the Fermi distribution takes a compact support at the Fermi energy.Therefore, we can evaluate the vertex function at the Fermi momentum k x F , to obtain for the bulk Onsager coefficients the simplied expressions (for i, j = 1, 2)

Nanoscale Advances
Paper where the total transport relaxation time, s (xl) tr (k x F ), is given by We remark that, along with the scattering relaxation time derived directly from the self-energy in eqn (12), the transport relaxation time in eqn (31) is inversely proportional to the concentration of dislocations s (xl) tr . The details of its derivation, as well as the computation of the integrals in eqn (30), are described in detail in the ESI.† In terms of these integrals, we nally obtain closed analytical formulae for the Onsager coefficients: and where in all those expressions, Li s (z) stands for the polylogarithm function of order s.The electrical conductivity is obtained from eqn ( 18) and ( 32) and then, we obtain Similarly, the electronic thermal conductivity is obtained from eqn (19)

Results
In this section, we shall evaluate our analytical expressions to estimate the transport coefficients of several WSMs in the family of transition metal monopnictides.For this purpose, we shall consider the microscopic/atomistic parameters obtained from rst-principles calculations, as reported in ref. 39 and 17.We shall also take into account the anisotropies reported by ref. 17 in the Fermi velocities and density of charge carriers at different Weyl nodes (x = ±) and bands (l = ±), respectively.These two references computed values for the Fermi energy with respect the position of each of the Weyl nodes as presented in Table 2.
We shall assume that the z-direction is aligned with the crystallographic direction of the defect axes, while the temperature and/or voltage gradients are imposed parallel to the xyplane.Therefore, we shall employ the x-and y-components of the Fermi velocities v (xl) F,a averaged from the reported values given in ref. 17 and 39 (see Table 3), for the conduction (l = +1) and valence (l = −1) bands, as well as for each of the chiral Weyl nodes (x = ±1), respectively.From the energies presented in Table 2 and the Fermi velocities given in Table 3, we can compute the Fermi momenta at each node using the formula given in ref. 17.
where k x W is the wave-vector location of each Weyl node in momentum space.The computed values of the Fermi momenta are shown in Table 4.
In order to estimate the geometric and structural parameters involved in the model, we follow the analysis presented in our previous work. 35Therefore, we assume that the cylindrical regions representing the dislocations have a radius a = 15 nm.From the proportionality relation between the torsional angle q (in degrees) and the pseudo-magnetic eld representing strain jB x ja 2 = 1.36qf0, 35 the modied ux quantum associated with the dislocations in these materials is approximately 1:5 300 $4:14 Â 10 5 T Å 2 z 330 T Å 2 .
Moreover, for deniteness, in this work we have chosen a torsion angle q = 15°.Finally, for the parameter a = V 0 /ℎv F that captures the effect of the delta barrier representing the lattice

Paper
Nanoscale Advances mismatch at the edge of the cylindrical dislocation, we follow our previous estimations based on Frank's law 35 by setting a = 3p/4.As clearly seen in eqn ( 36)- (38), our analytical expressions for the electronic transport coefficients depend on the total transport relaxation time due to the scattering with the dislocations at the Fermi energy.In Table 5, we present calculated values of such relaxation times, for different transition metal monopnictides, assuming the structural parameters in our model.For an estimation of the concentration of defects n d in real crystal systems, ref. 38 reports a native concentration of dislocations in the range n d ∼ 10 5 to 10 7 cm −2 for the materials TiO 2 and SrTiO 3 .These concentrations can be enhanced using different treatments up to 10 13 cm −2 , close to the rendering amorphous limit.Also, as is pointed out in ref. 44, the maximal practical density of screw dislocations detected in materials using electronic microscopy is in the range 10 11 to 10 12 cm −2 .Assuming then that a realistic concentration of dislocations would be in the range n d ∼ 10 9 to 10 11 cm −2 , an important aspect to check is if the ratio n d /(pk F 2 ) < 1 for the materials involved, in order for our approximations to be valid.From the values for the Fermi momenta reported in Table 4, we see that the four materials satisfy k ± F > 0.01 nm −1 .Therefore, for the aforementioned range of concentrations, the ratio n d /(pk F 2 ) ∼ 10 −3 to 10 −1 , and hence our approximations are well justied for all four materials analyzed in this study.
The DC conductivity and electronic thermal conductivity, as a function of the concentration of dislocations, are displayed in Fig. 2 and 3, respectively, where a temperature T = 5 K was assumed.We see that both transport coefficients exhibit an inverse proportionality, i.e.In particular, as seen in Table 5, for a concentration of dislocations n d = 2 × 10 9 cm −2 the relaxation times in all four materials are on the order of s tr ∼ 10 −12 s.As displayed in Fig. 4, the electrical conductivity of all materials depends on temperature, and assuming   Nanoscale Advances Paper a concentration of defects of n d = 2 × 10 9 cm −2 , it is in the range of s xx ∼ 10 4 to 10 5 U −1 cm −1 , with NbP being the better conductor.A similar hierarchy among the four materials is observed in Fig. 5 for the electronic thermal conductivity, which displays a nearly linear dependence up to room temperature.
Results computed from eqn (36) for the DC conductivity along xand y-directions at zero and room temperatures are presented in Table 6.We observe anisotropy between the x-and y-directions, due to the anisotropy in the components of the Fermi velocity, as can be appreciated in the values displayed in Table 3.
As seen in Fig. 5, the room temperature electronic thermal conductivity in all four compounds, assuming the same concentration of defects n d = 2 × 10 9 cm −2 , is on the order k (el) xx ∼ 10 to 10 2 W K −1 m −1 .The reason why NbP exhibits higher values of electrical and thermal conductivity as compared to the Fig. 3 Electronic thermal conductivity k (el) xx versus the concentration of dislocations computed from eqn (37) at T = 5 K, for the transition metal monopnictides TaAs, TaP, NbAs and NbP.Table 6 Values of the s xx and s yy DC conductivities (in the units of 10 3 U −1 cm −1 ) at 0 K and 300 K for each material.The result was computed from eqn (36) other materials in the monopnictide family is because its k F is almost an order of magnitude larger, as shown in Table 4.The effect becomes dominant due to the presence of k F in the exponential of the argument of the polylogarithmic functions in the analytical expressions for the electric and thermal conductivities.In physical terms, since such an exponential factor arises from the Fermi-Dirac distribution at nite temperature, it implies a higher population of chiral Weyl fermions available for transport as compared with the other materials in the same transition monopnictide family.Results calculated using eqn (37) for the electronic contribution to the thermal conductivity along x-and y-directions at room temperatures are presented in Table 7. Again, we observe anisotropy between the x-and y-directions, as a consequence of the anisotropy in the Fermi velocity components in these materials.
From the expressions for the electrical conductivity in eqn (36), the electronic thermal conductivity in eqn (37), and the Seebeck coefficient in eqn (38) In Fig. 7, we represent the Lorenz number, calculated from eqn ( 40) for all four materials as a function of temperature.Remarkably, L(T / 0) / L 0 = (p 2 /3)(k B /e) 2 , and hence the Wiedemann-Franz law is indeed satised in the limit of very low temperatures, a common feature for normal metallic systems, that is however also veried in these semimetal compounds.
Table 7 Values of the k (el) xx and k (el) yy electronic thermal conductivities (in the units of W K −1 m −1 ) at 300 K for each material.The result was computed from eqn (37)

Nanoscale Advances Paper
In Fig. 6, we present the Seebeck coefficient as a function of temperature, calculated from eqn (38) for the different materials.The negative sign of the Seebeck coefficient is consistent with the choice of a positive chemical potential, where the charge carriers are therefore electrons (instead of holes).For all materials, the Seebeck coefficient at room temperature is on the order of jSj ∼ 10 2 mV K −1 , and its absolute value grows at lower temperatures.This is consistent with different estimations in the literature for the family of transition metal monopnictides, reporting values in the range jSj ∼ 10 2 to 10 3 mV K −1 (ref.23-25)  at room temperature.
Finally, in Fig. 8 we present the gure of merit ZT (el) calculated from eqn (41), for all different materials as a function of temperature.As both the DC conductivity and the electronic thermal conductivity are inversely proportional to the concentration of dislocations, i.e.
parameter cancels in their ratio in eqn (41), and hence ZT (el) turns out to be independent of n d .However, a weak dependence on the presence of dislocations remains, since the scattering relaxation time is still a function of such defects through the scattering phase shis d m (k), as seen in eqn (31).Nevertheless, we could check that this effect also tends to cancel with the relaxation time upon taking the ratio leading to ZT (el) , and in practice this value becomes nearly independent on the presence of dislocations.Since according to Fig. 7 all four materials satisfy the Wiedemann-Franz law at very low temperatures L (el) (T / 0) / L 0 = 2.44 × 10 −8 V 2 K −2 , the low temperature limit of the gure of merit depends only on the Seebeck coef-cient, ZT (el) (T / 0) ∼ S 2 (T)/L 0 , and hence it decreases to zero as T / 0, and increases with temperature as seen in Fig. 8 for all four materials.Remarkably, near room temperature, TaAs presents ZT (el) > 2, which suggests that it could be an excellent candidate for thermoelectric applications.These ndings are compatible with previous studies that proposed generic semimetals for thermoelectric applications due to their relatively large Seebeck coefficients at room temperature jSj ∼ 10 2 mV K −1 , 26 in agreement with the order of magnitude of our current estimations for the monopnictides, as well as with independent estimations for these compounds, reported in the literature [23][24][25] to be in the range jSj ∼ 10 2 to 10 3 mV K −1 at room temperature.

Discussion
Our results for the electrical conductivity can be compared with independent estimations reported in the literature for the monopnictide family (TaAs, TaP, NbAs and NbP), where different scattering mechanisms where considered, particularly the electron-phonon interaction but not the dislocations studied in our work.The reported values, summarized in Table 1, show that the electrical conductivity is in the range s xx ∼ 10 4 to 10 6 U −1 cm −1 , in agreement in the order of magnitude with our results in Fig. 4 for an estimated concentration of dislocations n d = 2 × 10 9 cm −2 .Our theoretical model is concerned with the role of scattering with the quenched distribution of dislocations, but it does not include other possible mechanisms, particularly the electron-phonon scattering.However, its contribution to the transport relaxation time may be signicant as temperature increases enough to excite the relevant phonon modes.Moreover, the phonon spectrum itself can develop interesting topological features that may generate novel electron-phonon scattering mechanisms in WSMs, as discussed for instance in ref. 46.The latter is an entirely different mechanism, whose detailed analysis requires a separate model beyond the scope of the present work.The combination and competition between both scattering mechanisms can be estimated using Mathiessen's rule, such that the overall relaxation time including electron-phonon scattering would be 47 For instance, an estimation of the electron-phonon contribution is reported in ref. 18 and 25 as computed from rst principles for TaAs and NbAs, and NbP respectively.The reported electron-phonon relaxation times at 300 K are on the order s e-p ∼ 10 −13 seconds for all materials. 18,25In addition, ref.
23 reports an estimated value of the overall relaxation time (including electron-phonon, impurity and piezoelectric xx and (b) shows ZT (el) yy .
Paper Nanoscale Advances scattering) of s el = 3.01 × 10 −13 at 300 K, 23 which is still on the same order.In contrast, the calculated transport relaxation times for the scattering mechanism considered in this work, assuming a concentration of dislocations of n d = 2 × 10 9 cm −2 and displayed in Table 5, are larger, on the order of 10 −12 seconds.However, as shown explicitly in Fig. 2 and 3 for the transport coefficients, the corresponding values for the relaxation time are inversely proportional to the concentration of dislocations, and hence their relative importance in comparison with other possible scattering mechanisms is strongly determined by this sample-dependent parameter.On the other hand, electron-phonon scattering is strongly dependent on temperature, and hence as an estimation we can interpolate it from its reported value s e-ph (300) at 300 K, using the common Bloch-Gruneisen expression, 47 with Q D the Debye temperature and the function 47 In Fig. 9, we represent the transport relaxation time s tr (solid blue line) for the scattering with dislocations studied in this work, as a function of their concentration n d , exhibiting the expected inverse proportionality.For the sake of comparison, we also present as solid horizontal lines the values for the electron-phonon scattering relaxation time at three different temperatures s e-p (T), estimated from eqn (43).As can be seen in Fig. 9, at 300 K the scattering due to dislocations dominates over electron-phonon at concentrations n d > 3 × 10 11 cm −2 , while the corresponding concentration threshold is given by n d > 2 × 10 10 cm −2 at T = 100 K, and n d > 6 × 10 10 cm −2 at T = 50 K, respectively.
Concerning the thermal conductivity, for the sake of comparison, ref. 23 reports a rst principles calculation for the electronic contribution to the thermal conductivity of TaAs, including electron-phonon scattering but no dislocations as in this work.Their result at 300 K is k (el) xx = 56.87W K −1 m −1 , which is within the range of our calculated values as displayed in Fig. 5 and in Table 7, even when taking into account the anisotropy in the x-and y-directions already discussed.
Since the contribution from the overall relaxation time tends to cancel when taking the ratio of the transport coefficients in eqn (40), the Lorenz number obtained when other scattering mechanisms are present should still be close to our calculation.Indeed, ref. 23 also reports a value for the Lorenz number L (el) xx = 2.27 × 10 −8 V 2 K −2 , which is a small deviation from the Wiedemann-Franz law, in agreement with our results in Fig. 7.
Finally, in Fig. 8, we display the gure of merit ZT (el) calculated from eqn (41), for all different materials as a function of temperature.At room temperature, NbAs and TaAs exhibit a comparatively large gure of merit, with ZT (el) xx ∼ 2.5 for TaAs, suggesting that they could be excellent candidates for thermoelectric applications.We remark that, upon including the phonon contribution to the total thermal conductivity k Tot = k (el) xx + k (l) xx , this value will decrease.Indeed, a crude estimation of this effect may be introduced using the formula where we dened the correction factor due to the presence of the lattice conductivity by However, as can be appreciated in Table 1, at 300 K the lattice thermal conductivities reported in the literature are in general smaller than their electronic counterparts, with the exception of TaP, and thus the correction factor is not far from unity for most cases.Considering our calculated values for the electronic thermal conductivity (assuming n d = 2 × 10 9 cm −1 ) at 300 K presented in Table 7, along with the values for the lattice thermal conductivity reported in the literature at the same temperature as displayed in Table 1, we have that for TaAs C l = (1 + 36.06/26.17)−1 = 0.42, while for NbAs C l = 0.96, indicating that even when including the phonon effects, the gure of merit for those two compounds ZT ∼ C l × ZT (el) $ 1 is still comparatively large, thus suggesting that they could be very attractive for thermoelectric applications.Moreover, as the presence of the torsional dislocations studied in this work will affect the mechanical properties of the lattice, thus enhancing phononphonon scattering, the lattice thermal conductivity will decrease in the presence of such defects as compared with the literature values quoted in Table 1.Therefore, our values for the gure of merit may actually be closer to reality than our corrected estimations here, but a detailed analysis of phonon effects is a matter for a separate study.As a nal comment, we Fig. 9 The figure shows the logarithm of the relaxation time s (in seconds) vs. the logarithm of the concentration of dislocations n d (in cm −2 ) for the material TaAs.Horizontal lines correspond to the values of s e-ph (T) for three different temperatures.The value at T = 300 K is s e-ph (300) ∼ 10 −13 seconds. 18The values at T = 50 K and T = 100 K were computed using eqn (43), where Q D = 352 K is the experimental Debye temperature reported for TaAs. 45

Nanoscale Advances Paper
remark that no inter-valley scattering is involved in our analysis, for two main reasons.First, since the two valleys are well separated in momentum space, in order to couple them the potential scattering term must involve momentum exchange at least of this order of magnitude.On the other hand, in our formulation of the scattering problem across a single dislocation, the scattering term arises from the elastic gauge eld connection at the same valley as the spinor state being scattered, and hence this constitutes a kinematic constraint preventing inter-valley scattering.Nonetheless, the possibility for inter-valley scattering cannot be ruled out completely, since other mechanisms besides the ones considered in this work may be in place.

Conclusions
Along this article, we have presented a theoretical analysis for thermoelectric transport coefficients in the family of transition metal monopnictides, when the sole scattering mechanism considered is the presence of a uniform, diluted concentration of torsional dislocation defects.Our approach is based on a combination of Green's functions with a statistical average over the random distribution of defects, leading to a Dyson equation with a self-energy in the non-crossing approximation, enhanced with vertex corrections.Moreover, from the analytical expressions for the retarded and advanced Green's functions, by means of general Onsager relations in non-equilibrium thermodynamics and the Luttinger formalism to implement the Kubo formulae, we obtained explicit analytical expressions for the electrical conductivity, thermal conductivity (electronic contribution) and Seebeck coefficient.Our analytical expressions are fairly general, and complemented with geometrical and microscopic parameters obtained from ab initio calculations, we could evaluate them to estimate the corresponding values of those transport coefficients for each material as a function of temperature and concentration of dislocations.This work provides a rst step towards a theoretical analysis of transport in these systems, and hence we shall not delve into details concerning the experimental challenges posed, for instance, by the control of the sign of the charge carriers by a specic doping mechanism.Indeed, experimental and ab initio studies suggest that both electron and hole pockets will in general participate in transport. 48However, recent experimental studies 49 on WSM lms reveal that the presence of grain boundaries (that may play a similar role to our dislocation edges here) favours spontaneous predominance of holes (positive charge carriers) with a very high mobility, thus suggesting that doping may be achieved by defect engineering in these materials.
As presented in the Results and discussions sections, our analytical results predict values for the transport coefficients which are close to those reported in the literature, where different scattering mechanisms than this one are considered, particularly the electron-phonon interaction.For this particular case, we provided quantitative estimations of the range of temperatures and concentrations where each mechanism may become dominant.Finally, we remark that our results indicate that a high gure of merit is expected for at least two compounds in the family of transition metal monopnictides, i.e.TaAs and NbAs, even though our analytical expressions do not include the contribution from the lattice thermal conductivity.Nevertheless, based on reported values for this parameter in the literature, we estimate the order of the correction to our theoretical results to include the lattice effects, leading us to conclude that both TaAs and NbAs could be attractive candidates for thermoelectric applications.In this direction, we remark that the scattering mechanism analyzed in this work, i.e. the presence of torsional dislocation defects, introduces a moderate (as compared to electron-phonon) effect on electronic transport, whereas it may generate a strong phonon scattering mechanism decreasing the lattice thermal conductivity.Therefore, to engineer the concentration of such torsional dislocation defects in these materials, by decreasing the lattice thermal conductivity while nearly preserving the values of the electrical conductivity and Seebeck coefficient, may lead to even higher gures of merit than predicted here.However, an accurate assessment of the phonon transport mechanisms involved in this case goes beyond the scope of the present work and is a matter of outgoing research.
ab involves the linear superposition of band and chiral components.In particular, in the limit of low concentrations n d /pk F 2 < 1, the Onsager coefficients (for i, j = 1, 2) are given by s xx ∼ n d −1 and k xx ∼ n d −1 , since the transport relaxation time dened in eqn (31) is itself inversely proportional to the concentration of defects s (xl) tr ∼ n d −1 .

Fig. 2
Fig. 2 Electrical (DC) conductivity s xx versus the concentration of dislocations computed from eqn (37) at T = 5 K, for the transition metal monopnictides TaAs, TaP, NbAs and NbP.

Fig. 4
Fig. 4 The figure shows the DC electrical conductivity for the transition metal monopnictides TaAs, TaP, NbAs and NbP: (a) displays s xx and (b) displays s yy .Here we assume a concentration of dislocations of n d = 2 × 10 9 cm −2 .

Fig. 6
Fig. 6 Seebeck coefficient versus temperature computed from eqn (38) for the transition metal monopnictides TaAs, TaP, NbAs and NbP: (a) shows S xx and (b) shows S yy .

Table 2
Values of DE x F ¼ E F À E xW from ab initio computations reported in the literature.Here, E F is the Fermi level and E x W is the energy at the Weyl node x.We use the average of the two values

Table 3
Values of the Fermi velocity v (xl) F,a in the units of 10 5 m s −1 .In the valence band (l = −1) they correspond to hole velocities.We use the average of the reported values in ref. 17 and 39.Tables with the reported values are given in Section 4 of the ESI

Table 4
Values of k x F computed from the average DE x F given in Table2and the Fermi velocities given in Table3using eqn(39)