Spin-dependent magneto-thermopower of narrow-gap lead chalcogenide quantum wells

A semi-classical analysis of magneto-thermopower behaviour, namely, the Seebeck and Nernst effect (NE) in quantum wells of IV-VI lead salts with significant extrinsic Rashba spin-orbit coupling (RSOC) is performed in this report. In addition to the spin-dependent Seebeck effect that has been observed before, we also theoretically predict a similar spin-delineated behavior for its magneto-thermal analog, the spin-dependent NE. The choice of lead salts follows from a two-fold advantage they offer, in part, to their superior thermoelectric properties, especially PbTe, while their low band gaps and high spin-orbit coupling make them ideal candidates to study RSOC governed effects in nanostructures. The calculations show a larger longitudinal magneto-thermopower for the spin-up electrons while the transverse components are nearly identical. In contrast, for a magnetic field free case, the related power factor calculations reveal a significantly higher contribution from the spin-down ensemble and suffer a reduction with an increase in the electron density. We also discuss qualitatively the limitations of the semi-classical approach for the extreme case of a high magnetic field and allude to the observed thermopower behaviour when the quantum Hall regime is operational. Finally, techniques to modulate the thermopower are briefly outlined.

. The four corners of the outer and embedded inner squares in the two-paneled figure are symbolically marked by specific gradients (potential, in case of charge (C) and spin (S) and temperature for a thermal (Q) bias) and their current counterparts in presence of a magnetic field. The arrow-denoted linkages show several processes that may arise appearing as a flow of charge (j C ), spin (j S ), or thermal (j Q ) currents. Some processes lead to the creation of a potential or temperature gradient. The left figure (a) sketches the spin-driven transport processes; for instance, the spin Seebeck effect (SSE) is the generation of a spin voltage by a temperature gradient while the reciprocal spin Peltier effect (SPE) indicates a finite heat current arising from a spin injection. Similarly, the spin Nernst effect (SNE) describes a transverse spin current created on account of a longitudinal temperature gradient. The spin Hall effect (SHE) is the production of a finite spin current density that flows perpendicular to a charge current. Inverse SHE (ISHE) is the reciprocal process of SHE. The spin conductivity (SC) is spin current that flows for a finite spin accumulation. The more conventional effects, such as Nernst (NE), Ettingshausen (EE), Righi-Leduc (RLE), and Hall (HE) are observed for a finite H (the right panel (b) illustrates these collective effects) and constitute classical magneto-thermoelectrics. The Seebeck effect (SE) and its reciprocal, the Peltier process illustrate the interplay between charge and heat currents when the magnetic field is absent. The Nernst and Seebeck effects are described in this work; a quantitative determination is made by connecting the potential and temperature gradient via a tensor matrix. Lastly, the valley (V) is an additional degree of freedom and serves as the origin of an anomalous current (j V ) and the valley Hall effect (VHE) for a finite Berry curvature. The Berry curvature is the momentum-space analogue of a real-space magnetic field. A difference in carrier population establishes a valley voltage.
Scientific RePoRTS | (2018) 8:5972 | DOI: 10.1038/s41598-018-23511-2 discernible difference. The derived results, beyond the predicted RSOC-governed variance, reveal a direct relationship to the dominant energy-dependent scattering mechanisms operational in the film. Remarkably though, for processes that are energy independent, the transverse part of thermopower vanishes for a finite magnetic field. We also point out that a similar cessation of the transverse part can happen when either the product of the cyclotron frequency and the transit time approaches zero or for the case of a very high magnetic field. We bring to attention though the flaw in such a prediction of vanishing transverse thermopower for a high magnetic field by noting that the semi-classical approach in this case is an inaccurate portrayal of electron motion and ignores the quantization arising from formation of Landau levels (LLs). We qualitatively discuss the quantum aspect of this problem and the more accurate oscillating thermopower behaviour observed with the onset of the quantum Hall regime (that owes its genesis to the formation of LLs) for high magnetic fields. From the generalized magneto-thermopower results, the magnetic field free power-factor (the product of the square of the Seebeck coefficient and electric conductivity) of PbTe thin films is easily obtainable and show a larger contribution arising from the electrons with RSOC-induced spin-down polarization. Finally, we indicate possible ways to modulate thermopower including changes to the dimensional confinement and application of strain -effects that are mirrored in the effective mass of the electrons.

Results
The thermopower calculations begin with the pair of basic equations that connect heat and charge current. In its most general form, we can write 16 In Eq. 1, we identify the electric field as F, the heat flux is denoted by q, the electric resistivity is ρ while Q(Π) is the Seebeck/thermopower (Peltier) coefficient. The thermal conductivity is identified by the tensor, χ. These tensor quantities are simple constants for an isotropic model. However, an isotropic system can also acquire a tensorial character in presence of a magnetic field, H. For an isotropic system, the pair of equations (Eq. 1) is modified to To explain the additional terms H × j and H × ∇T, note that H is a pseudo-vector while F, ∇T, and q are polar vectors; the change to the heat and electric current, therefore, must be of the form given in Eq. 2 to be in conformity with the general relation that the vector product of a pseudo-vector and polar vector must be a polar vector.
The presence of a magnetic field, in addition, to the usual longitudinal thermopower can also aid in setting up an electromotive force (emf) or a thermal gradient. Two familiar examples in this regard are the Ettingshausen and Nernst effects; while the former describes the generation of a transverse temperature gradient (for a magnetic field directed out-of-plane and an in-plane current), the second phenomenon ( Fig. 2) produces an emf perpendicular to the direction of a temperature difference. As a quantitative illustration of the two aforementioned thermoelectric effects, we begin by writing out the heat equation (the second in the pair of equations in Eq. 2) in component form and set ∂T/∂x = 0; further, assuming that no heat and electric current flows along the y-axis (q y = 0, j y = 0), we have ν χ ∂ ∂ = ′ T y Hj . This is the Ettingshausen effect: A temperature gradient transverse to the flow of current exists in presence of a magnetic field (assumed here directed along the out-of-plane z-axis). The Ettingshausen coefficient is ν χ ′ / 1 . Likewise, for a finite temperature gradient, say, along the x-axis (∂T/∂x ≠ 0; ∂T/∂y = 0) and vanishing charge currents (j x , j y ), a potential drop (E y = ν 2 H∂T/∂x). develops. This is the Nernst effect and the corresponding coefficient is ν 2 H.
A few other notable examples of the influence of a magnetic field on thermoelectric behaviour are the Righi-Leduc effect and the adiabatic Nernst and Hall phenomena 17 . The quantitative determination of Figure 2. The Nernst effect which describes an electromotive force (emf) transverse to a temperature gradient in presence of a magnetic field (directed out-of-plane in this figure) is schematically shown here. The transverse emf as we show later is adjustable via electron dopant density and the transport scattering time.
The integrals in Eq. 3 can be simplified using the Sommerfeld expansion (low-temperature regime) to yield a more compact form: The derivative in Eq. 4 is evaluated at the Fermi energy. To use this compact form in presence of a magnetic field which introduces anisotropy, we must change to a tensor notation. The longitudinal and transverse thermopower expressions are therefore: Here, ρ = 1/σ, where σ is the magneto-electric conductivity tensor. The differentiation is indicated by the primed (superscript) notation. Likewise, for the transverse part, we have xx xy xy yy 2 2 In light of the above thermopower expressions, it is evident that we need to determine the magneto-electric conductivity which can be derived rigorously within the Boltzmann approximation; for a proof, see Chap. 8 in Ref. 19 . However, a single band model with parabolic energy bands should essentially offer an identical result. We simply quote it (and also for its reciprocal, the resistivity, ρ) here 20 In Eq. 6, the relaxation time is τ, the effective mass is indicated by m * , and ω c = eB/m * denotes the cyclotron frequency. The electron density is n in appropriate unit. The external magnetic field is included through the vector B = μ 0 H, where μ 0 is the magnetic permeability. Inserting the magneto-conductivity tensor expression in Eq. 5a and b, the two desired thermopower quantities (Q xx , Q xy ) can be computed. For an estimate of the relaxation time, we employ a semi-empirical approach that applies to non-degenerate semiconductors; this method varies the relaxation time as τ = λε s . Here, ε is the energy while λ and the exponent s are scattering mechanism specific material constants. For cases where this semi-empirical estimation of the relaxation time is reasonably accurate, the thermopower expressions can be readily obtained. Substituting for τ as λε s and evaluating the matrix product [ρ][σ′] whose elements appear in Eqs. 5a and b, we have In writing Eq. 7 that involves the derivative of the magneto-conductivity tensor, we used the relation Here, g(ε) is the two-dimensional density of states (see Eq. 19, Methods Section). The magneto-thermopower tensor components in Eq. 7 are defined at ε f , the Fermi level. In passing, observe that while the thermopower components seem independent of electron density, they however, are modulated via the n-dependent ε f , which is Eq. 22 derived in the Methods section. Further, it is easy to show using an analogous form of Eq. 5b that Q xy = −Q yx .
Before we present numerical estimates for thermopower, it is instructive to enunciate on few key aspects of the expressions (for the low B-field regime, where the classical determination of the magneto-conductivity tensor remains valid) derived in Eq. 7. First of all, note that the presence of a magnetic field (directed out-of-plane, along z-axis) induces a transverse heat gradient and the thermopower (Q), has a finite off-diagonal component, Q xy . Further, for impurity scattering events that typically do not show a strong functional dependence on energy, a possible scenario when electrons are located close to the conduction band minimum, the exponent s can be approximated to zero. A proof of this for the case of delta-doped non-interacting impurity scattering appears in the Methods section. Setting s = 0 in Eq. 7, we recover thermopower expressions under a vanishing magnetic field 21,22 . For the other extreme case of high-field limit, such that ω c τ → ∞, it is easy to verify that Q xy ceases to exist, much like when the magnetic field is turned off. However, in this case, we must exercise caution (also see the discussion presented in the summary) that the expressions in Eq. 7 may be incorrect as regards to their formulation goes at high magnetic fields. The quantization of the electron orbit via formation of quantum mechanical Landau levels beginning at high magnetic fields and the possibility of the onset of quantum hall features radically alters the physical description and therefore any analysis (similar to the one presented) in this regime that has its genesis in classical arguments would be mostly erroneous. Further, note that expressions for the thermopower tensor components in a magnetic field can be reasonably well-estimated if accurate material-specific values such as the effective mass and relaxation time (τ) are at hand. The thermopower can be modulated by tailoring the Scientific RePoRTS | (2018) 8:5972 | DOI:10.1038/s41598-018-23511-2 relaxation time to a variety of scattering events, reflected here through adjustments of the exponent, s. As an illustration of the process that fixes the value of s, we take up the specific case of non-interacting impurities and compute the relaxation time using a self-consistent Born approximation (SCBA). The SCBA (see Eq. 24 in the Methods section and the accompanying discussion that follows) furnishes the relaxation time via the imaginary part of the self-energy of electrons whose motion is disturbed by the impurity scattering potential 23 . In this regard, a calculation of the relaxation time via SCBA and other required parameters must proceed by setting up an appropriate Hamiltonian for the thermoelectric material accounting for dimensional-effects and composition of the bulk/nanostructure. The numerical diagonalization of the pertinent Hamiltonian supplies the dispersion relation useful in the evaluation of the effective mass, assigning the relevant energy scale of the problem, and an assessment of the part played by scattering that alters carrier mobility. The Methods section outlines the Hamiltonian (Eq. 16) for quantum wells of lead salts, the target thermoelectric material in this report.
As for the semi-empirical approach, it is reasonable to state using the form of Eq. 7 that an enhanced thermopower is achievable if the exponent (s) acquires a higher number; in fact, suggestions 24 have been made to dope PbTe, a well-studied thermoelectric such that s is reinforced via impurity-scattering. This scheme, however, works insofar as there is no obvious degradation of electron mobility., To see this, recall that the efficient generation of thermoelectric currents is normally typified via the figure-of-merit (ZT) which relies on electric conductivity and can suffer a considerable reduction in an attempt to augment the scattering. However, it is remarkable that a boost to scattering (and consequently a higher s) led to an observable increase in ZT for bismuth structures 25 . In this context, also note how for scattering events which do not have an energy dependence (s = 0) or simply a constant relaxation time, the transverse magneto-thermopower or the so-called Nernst-Ettingshausen effect vanishes. In practice, scattering events are always energy-dependent phenomena; however, as an approximation, Bloch electrons close to the conduction band minimum which are scattered by constant spherical delta potentials of impurities, the relaxation time can be shown to be a fixed number. The relevant relaxation time calculation which proves this constancy is presented in the Methods section (Eq. 24 and the concomitant analytic results).
Finally, observe from Eq. 7 the possibility of carrier ensembles displaying a thermopower response to temperature and potential gradients tied to their respective Fermi levels. A common occurrence of distinguishable Fermi levels within a material can be readily recognized for spin-split Bloch electrons, where the splitting emerges either from an intrinsic or the Rashba spin-orbit coupling (soc) (RSOC). In general, the intrinsic soc is significant in heavy metals with a large nucleus. For nanostructures such as quantum wells carved out of compounds of heavy metals which additionally possess a narrow band gap, a pronounced RSOC splitting has been observed. This splits the Fermi surface into two groups of spin-polarized electrons. A simple realization of this can be seen (see Fig. 3) by considering the minimal quadratic Hamiltonian with the linear Rashba contribution 26 where λ > 0 R is the Rashba coefficient and determines the robustness of the splitting. The related dispersion is of the form : ε = ħ 2 k 2 /2m * ± λ R k; the splitting is therefore 2λ R k, where k is the in-plane wave vector given by + k k x y 2 2 . The Fermi surfaces for the spin-up and spin-down carriers are therefore non-degenerate and the difference in their energies determined by λ R . A straightforward calculation connects the dependence of the surface Fermi energy (ε f ) to RSOC via λ R for a preset electron density (n). We quote the result here (see Methods section for a proof). The upper (lower) sign is for the spin-up (down) carriers.
The upper (lower) sign is for spin-up (down) electrons and α = ħ 2 /2m * . The Fermi energy (for a specific n) is evidently governed by the strength of RSOC. The corresponding vector (for later use) is The upper (lower) sign in Eq. 10 is for the spin-up (down) ensemble. A more elaborate set of remarks on this aspect and numerical calculations in context of realistic quantum well systems is included in the Methods section.
We can now insert the Rashba-split Fermi energy expressions (Eq. 9) in Eq. 7 and conceivably expect a spin-dependent magneto-thermopower. This could be essentially thought of as the magneto-thermal analog of the spin-dependent Seebeck phenomenon and may likewise be referred to as the spin-dependent Nernst-Ettingshausen effect. In what follows, we employ this framework to examine the spin-delineated magneto-thermoelectric tensor quantities and the related power factor (for a zero B field case) of lead-salt quantum wells.
Lead salt films. We alluded to the utility of tellurides above as desirable thermoelectric materials. For numerical calculations we select PbTe as the candidate telluride, the choice of which is primarily driven by its well-regarded thermoelectric performance and a narrow band gap coupled to a strong spin-orbit coupling. The latter two attributes makes it suitable to observe spin-dependent thermoelectrics. The derived spin-resolved magneto-thermopower results here are generalized and can be applied to a wide variety of materials by selecting the case-specific effective mass and the correct form of the RSOC Hamiltonian. It is worth noting that the lead chalcogenides, PbS and PbSe, are iso-structural with PbTe, crystallizing in a rock-salt crystal arrangement and are expected to show similar thermoelectric behaviour. Briefly, the lead salts crystallize in the rock-salt structure (a two-atom basis) which consists of two inter-penetrating face-centred-cubic (fcc) lattices displaced by (0.5, 0.5, 0.5) × a, where a is the lattice constant 27   substrate surface (see Fig. 4) while the rest three (referred to as the oblique-valleys) are tilted by θ = 70.5°. Quantum confinement in a film, however, removes the degeneracy at the L point of the Brillouin zone pushing the three oblique valleys to a energy beyond that of their sole longitudinal counterpart. In this report, we ignore such distinctions and work with the longitudinal valley tacitly assuming that the other three equally contribute to magneto-thermopower.
For a numerical evaluation of the magneto-thermopower in thin films of lead salts, we start by fixing certain relevant quantities/material parameters. The carrier (electron) concentration 28  For this work, we extrapolate to believe that doping densities considered for AlGaAs FETs will continue to hold good for the current set of materials. Note from Eq. 9 the presence of two spin-polarized Fermi surfaces for a definite concentration. The variable parameters are the carrier effective mass (m * ), the magnetic field strength (H), and the choice of scattering mechanism (s) that influences the relaxation time. As we qualitatively discuss later, m * is calibrated by selecting different film thicknesses and alloying (x is the mole fraction) with SnTe to yield Pb x Sn 1−x Te. The implication of using Pb x Sn 1−x Te in addition to a ready source of composition-controlled tunable energy-levels lies in a band gap closing feature beyond a threshold tin concentration. This is discussed more fully in the closing section. For definiteness, we consider pristine PbTe films and plot the spin-resolved longitudinal and transverse magneto-thermopower of conduction electrons in a 6.0 nm thick film which is modeled as a quantum well (QW). The effective mass of the longitudinal L-valley conduction electrons in Eq. 8 is calculated starting from a 4 × 4 k.p Hamiltonian adapted for a quantum well. The necessary steps are presented (Eqs. 14 and 16) in the Methods section. The external Rashba spin-splitting (RSOC) in each case is determined assuming an asymmetry-induced out-of-plane electric field (F) operates; setting F = 10 6 V/m, we obtain the Rashba coefficient (Eq. 15 in the Methods section) and the spin-polarized energy states. With these simulation parameters in mind, the quantities Q xx and Q xy are shown for several magnetic field strengths in Fig. 5. For the 6.0 nm PbTe thin film structure (QW), the transverse effective mass of the conduction electrons and the direct band gap were determined to be 0.0565*m 0 (the free electron mass is m 0 ) and 0.2131 eV.
We now comment on the nature of the plots and choice of other variables used: First of all, observe that the longitudinal thermopower does not change sign as the magnetic field is reversed while the transverse component is odd in B, a result which readily follows from the nature of the conductivity expressions in Eq. 6. The longitudinal conductivity and resistivity are even while their transverse counterparts are odd functions of B bestowing an even and odd character to each term in Q xx and Q xy , respectively. The relaxation time in magneto-thermopower expressions was set to 0.1 ns and the scattering exponent s = 0.7. In this very approximate scheme of relaxation time determination, we are heuristically guided by standard empirical formulations that predict scattering exponents of s = 0.5 and s = 2 for neutral and ionized impurities, respectively 19 . The scattering of conduction electrons initiated by acoustic phonons is modeled by letting s = −0.5. Since most phonon modes would be suppressed at low temperatures, we choose an intermediate value of s = 0.7 which lies between the assigned exponents (0.5 < s < 2) for neutral and charged impurity scattering. Lastly, observe that for the operational Rashba spin splitting, the transverse magneto-thermopower is closely matched for both spin ensembles in contrast to the longitudinal case where the spin-up conduction electrons exhibit a higher value. This result for the longitudinal case is in agreement with the expression derived in Eq. 7 that shows an inverse relationship to carrier energy, and thereby translating into a higher magneto-thermopower for the spin-up electrons. Note that the spin-up electrons occupy an energetically higher Fermi surface in comparison to their oppositely polarized counterpart. In context of Fig.  5, the QW conduction electron concentration was assumed to be n = 10 12 cm −2 from which the RSOC splitting energy (2λ R k) is 12.0 meV for |k| = 1.0 Å. In the context of RSOC-induced splitting, we briefly digress here to emphasize on a key PbTe material parameter. Firstly, note that the Rashba parameter (λ R ) is given by the relation λ R = λ 0 E z , where λ 0 is a material-specific number (see Eq. 15, Methods section) and E z is the magnitude of a structural inversion asymmetry (SIA) generated 29 out-of-plane electric field. For estimating E z , supposing that it arises only from the charge density in the QW, it is given as E = en s /ε, where n s is the charge density and ε is the static dielectric constant; for PbTe, the static dielectric constant is ≈420. The implications of this large dielectric constant are profound and briefly touched upon in the following sub-section. Before closing, to be clear again, the spin-dependent thermopower expressions are a single-particle phenomenon where two channels (the spin-up and spin-down) independently carry heat and spin, unlike, say, for example, the spin-Seebeck effect which is believed to be a magnon-driven spin current affair 6,30 . Power factor of Pb-based thin films. An important marker in gauging thermoelectric processes is the power factor defined (for zero magnetic fields) as the product of Q 2 and electric conductivity. The longitudinal and transverse Seebeck coefficients (for the zero magnetic field case) and can be easily obtained by setting the cyclotron frequency (ω c ) to zero in Eq. 7. Straightforwardly, the transverse thermopower vanishes while the longitudinal component (Q xx ) reduces to: The utility of the Seebeck coefficient lies in the formulation of the figure of merit generally expressed as = , where σ is the electric conductivity and κ is the combined lattice and electronic contribution to thermal conductivity. Evidently, for a higher ZT, which is desirable for improved thermoelectrics, the product Q 2 σ must be maximized while ensuring a low thermal conductivity. A concurrent fulfillment of this dual set of conditions is purportedly a difficult proposition, however, lead salts offer much promise and multiple pathways in achievement of this goal 31,32 . Expanding on this point, first notice that electric conductivity σ is a charge density dependent quantity (σ = neμ) which in case of Pb-salts is augmented by the valley degeneracy (g = 4) of the L-valley high-symmetry point in the Brillouin zone. The conductivity also receives more augmentation from the mobility (μ) which is large for low effective masses. Typically, low effective masses are an outcome of the strong intrinsic spin-orbit coupling, a feature, as we have mentioned before is a defining hallmark of lead compounds. Unsurprisingly, the effective masses are indeed low for lead salts. In alliance with the favourable microscopic arrangement, the highly polarizable bands in lead salts give rise to a large static dielectric constant 33 leading to a significant Bohr radius (a B ). A large a B serves to conceal impurities and imperfections that may lower the overall mobility and therefore supports a high electric conductivity. It is then tenable to expect a substantial power factor for thin films (QW) of lead salts. In the following, we describe a set of steps to estimate it; as before, PbTe is the representative material.
For a semi-classical approximation of the conductivity, a linearized form of the Boltzmann equation within the relaxation time approach allows us to write . The Fermi velocity is v f = ħk f /m * , where k f is given by Eq. 10. Note that in expressing the velocity operator, we have ignored the Rashba contribution and only retained the parabolic part of the Hamiltonian. The relaxation time can be gauged by considering impurity scattering and excluding any phonon-assisted disruption. The relaxation time for such a condition can be written as The quantity n i is the impurity density and v i the corresponding strength. Putting Eq. 11, the zero-temperature form of Eq. 12, and the expression for τ together, the power factor (Q 2 σ) assumes the form where the upper (lower) sign is, as usual, for the spin-up (down) ensemble. We have also approximated the DOS discarding the Rashba contribution and subsituted for D(ε) as m * /2πħ 2 . For a quantitative prediction of the power factor, we use a 6.0 nm wide PbTe QW. The relaxation time (τ) is approximated by setting the impurity density to n i = 2.2 × 10 9 cm −2 while the corresponding potential is assigned the value, v i = 0.1 eVA 2 . As before, the L-valley conduction electron effective mass for this QW is m * = 0.0565m 0 . The relaxation time using the relation noted above for the selected impurity parameters and effective mass is roughly 1.0ns. The power factor for a 6.0 nm PbTe QW is displayed in Fig. 6. The low-temperature power factor calculation presented here represents the optimal case; impurities and related surface scattering events tend to lower the conductivity. The experimental power factor is therefore expected to be reduced vis-à-vis the current theoretical estimate in this report. In fact, the highest reported power factor in PbTe QWs (n ≈ 10 24 cm −3 ) at room temperature is 130 μW/cmK 2 31 . Further, notice that the effective mass appears in Eq. 13 pointing to the distinct possibility of reduction in QW dimension or greater confinement, for instance, in a nanowire to modulate the power factor. Indeed, such observations have been extensively reported to tune the power factor; for another line of reasoning on the aspect of dimensional analysis of this problem that also takes into account the thermal de Broglie wave length, see Ref. 34 .

Final Remarks
We have presented a semi-classical analysis of magneto-thermopower of a two-dimensional electron gas system in a PbTe quantum well. The calculations were done in the low-temperature regime and moderate magnetic fields extending up to B = 3.0T on either side of the zero mark. Quantum mechanical effects were incorporated by extracting a realistic effective mass for the conduction electrons of the PbTe quantum well from an appropriate k.p Hamiltonian. However, it is important to clarify that at higher magnetic fields than those considered here, there exists a sharp departure from classical behaviour. The electrons begin to circulate in Landau levels (LL) and may further quantize to reach the integer quantum Hall (IQH) regime. We have ignored the quantum mechanical coupling of the magnetic field to electron motion and the attendant changes to longitudinal and transverse magneto-conductance, especially when the IQH effect sets in. Succinctly, for a pre-defined Fermi level positioned between two consecutive LLs -at a high magnetic field and low-temperature -the transverse conductivity is quantized as σ xy = νe 2 /h, where ν is the filling factor while the longitudinal tensor component ceases to exist (σ xx = 0). By sweeping the Fermi level between energetically higher LLs, the well-known stair-like behaviour indicating a larger ν is observed. For situations where the Fermi level is placed between two successive LLs, the constant conductivity, by a simple application of Eq. 7 leads to a vanishing transverse magneto-thermopower. Simultaneously, the longitudinal part (Q xx ) is zero by virtue of σ xx = 0 in the IQH regime. In the other scenario, when the Fermi level is adjusted to align with a certain LL, both the transverse and longitudinal thermopower are restored as neither σ xx vanishes nor does σ xy display any quantization effects. Briefly, the successive vanishing and renewed establishment of thermopower as the Fermi level is progressively adjusted between successive LLs give an overall oscillatory pattern. An oscillatory plot that is in conformity with the reasoning expressed here was presented by Jonson and Girvin in Ref. 35 for two-dimensional electron gas in GaAs in the IQH regime.
In summary, we have theoretically analyzed the magneto-thermopower of lead salts that crystallize in the rocksalt structure. Tellurides of Pb-salts in various forms have been recognized as thermoelectrics with high efficiency, which can be further enhanced through alloying, dimensional confinement, tuning the carrier effective mass, and amending their large dielectric constant. Beyond these physical quantities and design parameters, the narrow band gap and strong intrinsic spin-orbit coupling contribute to a very sizable Rashba spin-splitting; within the scope of the presented model this leads to quantitatively different thermopower and power factor numbers. It is also fitting to note here that quantum wells of Pb-salts, in particular, PbTe alloyed with SnTe show topological crystalline insulator behaviour; a quantum state of matter that has been theoretically predicted (and Figure 6. The spin-resolved power factor (pf) of a 6.0 nm wide PbTe quantum well is plotted for T = 1 K as a function of carrier density. The pf is higher for the spin-down ensemble. This is also reflected in the enhanced charge conductivity (inset) of the spin-down branch calculated using Eq. 12. No electrostatic screening was included in the numerical estimation of the power factor and conductivity. The impurity scattering time (τ) for this calculation was set to 1.0 ns.
Scientific RePoRTS | (2018) 8:5972 | DOI:10.1038/s41598-018-23511-2 observed in few experimental demonstrations) to offer a much higher ZT than hitherto possible with conventional materials. A future work with alloyed PbTe in topological crystalline insulator phase and the attendant magneto-thermopower may uncover unexplored behaviour. As a final observation, these calculations can be transferred (with a modified minimal Hamiltonian of the type used in Eq. 8) to another well-known class of thermoelectrics -the rhombohedrally crystallizing binary tetradymites -that carry similar properties as the lead salts including a large spin-orbit coupling (this drives the creation of topologically non-trivial surface bands), low band gap, and high dielectric constant. A notable tetradymite is Bi 2 Te 3 that complies with all the aforementioned PbTe attributes, however, the presence of topological surface states in thin films complicates the measurement and interpretation of experimental data. The Bi 2 Te 3 structures are also beset by a large concentration of impurities that enter the growth process modifying the intrinsic character.

Dispersion calculations. A minimal Hamiltonian that captures the basic band dispersion of the conduction electrons in a lead chalcogenide quantum well is
RS R x y y x 2 where we have included the linear Rashba spin-orbit interaction. The Rashba coupling constant is λ R > 0 and is particularly strong in narrow band gap materials with strong intrinsic spin-orbit coupling, such as the lead chalcogenides. In particular, the strength of the Rashba coupling coefficient is λ R = λ 0 〈E(z)〉, where 〈E(z)〉 serves as the average electric field. The material-dependent λ 0 is given as 36 In Eq. 15, the fundamental direct band gap at is E g , the spin-orbit splitting is Δ, and m * is the effective mass. Note that the Dresselhaus coupling term is absent since bulk inversion symmetry is preserved for crystals which have rock salt crystal arrangement. The effective mass of band carriers in a lead chalcogenide quantum well can be derived from a continuum bulk 4 × 4 k.p model valid in the vicinity of the high-symmetry L-valley. In the spin-resolved (↑↓) basis set of ↑ − L 6 , ↓ − L 6 , ↑ + L 6 , and ↓ + L 6 , where ± L 6 are the Bloch functions of the conduction (−) and valence band (+) edge, it takes the form 15 vector serves as the z-axis, the k x , k y , and k z terms must be replaced by performing the operation Pk, where P is the transformation matrix and k is the k-vector triad. The transformation matrix, P for one of the equivalent tilted valleys, say, [111] can be written noting that its axis can be made parallel to that of the longitudinal valley, [111], through an anti-clockwise rotation of θ = cos −1 (1/3). Note that the orthonormal coordinate system in this case is = − x i j ( )/ 2, = − + + y i j k ( 2 )/ 6 , and = + − z i j k ( ) / 3. There exists, in addition to [111], another set of tilted equivalent valleys, namely, [111] and [111].
Relation between Fermi energy and Rashba coupling parameter. An outcome of the inclusion of the spin-orbit interaction Hamiltonians is that Fermi energy now depends on the strength of the Rashba spin-orbit coupling parameter. We quoted the results in the manuscript (Eqs 9 and 10), a quantitative calculation is given here. To prove this, we first derive the density of states (DOS) beginning with standard expression D(E) = π 1 4 2 ∫d 2 kδ(E − ε(k)). Note that the Rashba-split conduction electron energy is of the form E = αk 2 ± λ R k, where α = ħ 2 /2m * . To evaluate the δ(⋅), we recall the identity δ(g(k)) = δ − | ′ | k k g k ( ) ( ) i i , where k i is the non-degenerate root of g(k). For our case, g(k) = E − αk 2  λ R k, where the  differentiates the spin-up function from the spin-down version. The parameters α and β are always positive. We retain the positive root for g(k) in each case giving  it is easy to see that the spin-up branch has a lower DOS than its oppositely spin-polarized description. Finally, when α λ  k k R 2 (or large energies), the DOS approaches the value = πα πħ ⁎ m 1 4 2 2 . This is the standard DOS expression for a two-dimensional system without spin-degeneracy. We next calculate the electron density, n, at Fermi energy (zero temperature), which by definition is Employing the standard relation : x i x 1 1 , we arrive at i i 2 2 1 2