Charge disproportionation triggers bipolar doping in FeSb2-xSnxSe4

2 days ago - In this study, we demonstrate the ability to tune magnetic ordering in the ... crystal lattice leading to large separation between magnet...
0 downloads 0 Views 6MB Size
Subscriber access provided by UNIV OF CAMBRIDGE

Article

Charge disproportionation triggers bipolar doping in FeSb2-xSnxSe4 ferromagnetic semiconductors enabling temperature-induced Lifshitz transition Honore Djieutedjeu, Juan S. Lopez, Ruiming Lu, Brandon Buchanan, Xiaoyuan Zhou, Hang Chi, Kulugammana G. S. Ranmohotti, Ctirad Uher, and Pierre F. P. Poudeu J. Am. Chem. Soc., Just Accepted Manuscript • DOI: 10.1021/jacs.9b01884 • Publication Date (Web): 10 May 2019 Downloaded from http://pubs.acs.org on May 12, 2019

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Charge disproportionation triggers bipolar doping in FeSb2-xSnxSe4 ferromagnetic semiconductors enabling temperature-induced Lifshitz transition Honore Djieutedjeu1, Juan S. Lopez1, Ruiming Lu1, Brandon Buchanan1, Xiaoyuan Zhou2, Hang Chi2, Kulugammana G. S. Ranmohotti1,3, Ctirad Uher2 and Pierre F. P. Poudeu1,* 1) Laboratory for Emerging Energy and Electronic Materials (LE3M), Department of Materials Science and Engineering, University of Michigan, Ann Arbor, MI, 48109, USA. 2) Department of Physics, University of Michigan, Ann Arbor, MI, 48109, USA 3) Division of Science, Mathematics and Technology, Governors State University, University Park, IL, 60484, USA.

*Corresponding Author: Pierre Ferdinand Poudeu Poudeu E-mail: [email protected]; Fax: +1-734-763-4788

ACS Paragon Plus Environment

1

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 44

Abstract Ferromagnetic semiconductors (FMSs) featuring high Curie transition temperature (Tc) and strong correlation between itinerant carriers and localized magnetic moments are of tremendous importance for the development of practical spintronic devices. The realization of such materials hinges on the ability to generate and manipulate a high density of itinerant spinpolarized carriers and the understanding of their responses to external stimuli. In this study, we demonstrate the ability to tune magnetic ordering in the p-type FMS FeSb2-xSnxSe4 (0 ≤ x ≤ 0.20) through carrier density engineering. We found that the substitution of Sb by Sn FeSb2xSnxSe4

increases ordering of metal atoms within the selenium crystal lattice leading to large

separation between magnetic centers. This results in a decrease in the Tc from 450 K for samples with x ≤ 0.05 to 325 K for samples with 0.05 < x ≤ 0.2. In addition, charge disproportionation arising from the substitution of Sb3+ by Sn2+ triggers partial oxidation of Sb3+ to Sb5+, which is accompanied by the generation of both electrons and holes. This leads to a drastic decrease of the electrical resistivity and thermopower, simultaneously with a large increase in the magnetic susceptibility and saturation magnetization upon increasing Sn content. The observed bipolar doping induces a very interesting temperature-induced quantum electronic transition (Lifshitz transition), which is manifested by the presence of an anomalous peak in the resistivity curve simultaneously with a reversal of the sign of majority charge carriers from hole-like to electron-like at the temperature of maximum resistivity. This study suggests that while there is a strong correlation between the overall magnetic moment and free carriers spin in FeSb2-xSnxSe4 FMSs, the magnitude of the Curie temperature strongly depend on the spatial separation between localized magnetic centers rather than the concentration of magnetic atoms or the density of itinerant carriers. Keywords:

Quantum

electronic

transition;

Bipolar

doping;

Spintronics;

High-Tc

ferromagnetism; Transition metal chalcogenides

ACS Paragon Plus Environment

2

Page 3 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Introduction Manipulating magnetic ordering in various magnetic semiconductors by chemical (doping, solid solution) or physical (external stimuli such as applied bias voltage or magnetic field) means has been intensively investigated,1-3 with the goal to demonstrate the ability to use both the charge and spin degrees of freedom of free carriers in spintronics devices.1,4 Important milestones towards realizing this goal are the ability to control the density and transport of spin-polarized charge carriers as well as understanding their interactions with localized magnetic moments in magnetic semiconductors, and particularly in high temperature ferromagnetic semiconductors (FMSs).5-9 In this respect, significant progress has been made over the past three decades on dilute ferromagnetic semiconducting (FMS) material systems.4,10-13 For instance, nonmagnetic conventional semiconductors, such as III-V, II-VI, and IV-IV compounds can be turned into magnetic semiconductors (MS) by introducing sizable concentrations of magnetic ions, such as Mn, Fe, Ni and Co into the crystal lattice. In addition, the electronic properties of such conventional semiconductors are modulated by the incorporation of the magnetic impurities. This approach has been used for several decades with III-V, II-VI and IV-IV semiconductors.1,4,14,15 While III-V semiconductors, such as GaAs are already being used in a large number of electronic devices, the integration of the magnetically doped materials (Ga1xMnx)As

into practical devices remains restricted by the low magnetic transition temperature

(173 K) and low solubility of the magnetic ion into the host semiconductor.4,13 Nevertheless, such multifunctional materials could find application in devices where the magnetic properties can be controlled through optically activated charge carriers,16 such as the spin-based field effect transistors (FETs), spin-polarized lasers, light emitting diodes (LEDs), non volatile memory and quantum computing. In magnetically doped conventional semiconductors, extrinsic carriers (electrons or holes) generated within the valence band (VB) or the conduction band (CB) are coupled to the concentration of the magnetic impurity (Mn, Fe, etc.), which serves both as a doping agent (donor or acceptor) and as the magnetically active species. Fascinating collective carrierinduced or carrier-mediated ferromagnetic phenomena, which incorporate the effects of phonons, electrons and spins on the overall magnetic and transport behavior, have been reported in magnetic semiconductors and diluted magnetic semiconductors.2,17-19 However, the mechanism of the exchange interaction between the spin of charge carriers and the localized moment of magnetic atoms remains controversial.10,13 In the case of carrier-induced ferromagnetism, critical carrier concentration has been reported for the onset of ferromagnetism and for magnetic transitions such as AFM to FM ordering and FM to PM disorder behavior. For

ACS Paragon Plus Environment

3

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 44

instance, it was found in [(PbTe)1-x(SnTe)x]1-y(MnTe)y that with a very small value of y, the FM onset hole concentration was p ~ 2 x 1020 cm-3, while the maximum Tc value was observed with samples exhibiting hole concentration of p ~ 7 x 1020 cm-3.20 Ferromagnetic ordering in this system was rationalized by considering that the number of effective carriers taking part in the RKKY exchange interaction mechanism increases when the electrons from L valence band are four times more mobile than those in Σ band.21,22 In material systems with s-d and p-d orbital hybridization, such as p-type magnetic semiconductors with large hole density (Mn-doped GaAs), it is widely accepted that the FM ordering of adjacent magnetic centers (Mn – Mn) within the host material is mediated by very strong spin-spin interaction involving free carriers. However, the intrinsic contribution of individual material parameters (carrier density, carrier mobility, fraction of magnetic centers (concentration of Mn) or the distance between magnetic centers) to the change in the Curie temperature (Tc) remains unclear.10,13,23 In addition, the exact band picture resulting from the substitution of Ga by Mn in (Ga,Mn)As remains controversial.10

24

This limited understanding

of the mechanism of ferromagnetism in Mn-doped FMSs is due to the fact that both the density of holes and the concentration of magnetic centers in FMSs are modified simultaneously when Mn ions substitute for Ga in (Mn,Ga)As. In addition, it is challenging to control the Mn – Mn separation within the host material as well as the exact concentration of Mn atoms contributing to the ferromagnetic order in Mn-doped FMSs. This is due to the fact that some of the Mn ions enter into interstitial sites rather than occupying the crystallographic position of the substituted atoms.25-27 This limited understanding of the mechanism by which holes mediate ferromagnetic ordering and affect Tc in Mn-doped FMSs severely hampers our ability to use spin-dependent electronic transport to manipulate magnetic ordering and to develop efficient ferromagnetic semiconductors suitable for spin-based devices at practical temperatures.

ACS Paragon Plus Environment

4

Page 5 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

To achieve a greater level of understanding of key parameters (carrier concentration, concentration of magnetic atoms, distance between magnetic atoms, etc.) controlling the correlations

between

electronic transport and magnetic

behavior

in

magnetic semiconductors,

it

necessary

design

compounds symmetry

to with

is low

crystal

structure in which the substructure responsible for magnetic behavior is decoupled Figure 1: Schematic illustration of the crystal structure of FeSb2-xSnxSe4 projected along the b-axis. The relevant semiconducting and magnetic substructural units are highlighted. The M1 and M2 sites (blue) are mixed occupied by Sb 96% and Fe(4%); whereas the M4 site (cyan) is mixed occupied by Fe (84%) and Sb(16%). The M3 site (green) is fully occupied by Fe.

from

the

substructure responsible for

electronic

conduction. Such low symmetry

structure

enables flexibilities in independently tuning the carrier

concentration

through doping within the semiconducting substructure, the concentration as well as separation between magnetic atoms through solid solution substitution within the magnetic substructure. Along these lines, our research group has developed over the past decade several complex structures with tremendous compositional flexibilities.5,28-37 These include (M2-xPbx)Pb2Sb4Se10 (M = Mn, Fe);36,38 FePb4Sb6Se14;34 and the homologous series M2Pn4M’N-3SeN+5 (N  3 (M = Mn, Fe; M’ = Pb, Sn; Pn = Sb, Bi).5,30,33,35-37 The fascinating class of ternary compounds MPn2Se4 (Figure 1) corresponds to the first member of the M2Pn4M’N-3SeN+5 series with N = 3. The monoclinic crystal structure (space group C2/m) contains four crystallographically independent metal positions, M1(4i), M2(4i), M3(2d) and M4(2a) that are distributed within two distinct building units alternating along [001] denoted layer A and layer B (Figure 1). The layer A, which contains metal positions M1 and M3, consists of double rods of face-sharing monocapped trigonal prismatic coordination of Se atoms around M1 site that alternate along the

ACS Paragon Plus Environment

5

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 44

a-axis with M(3)Se6 octahedra that shared edges to form a single chain parallel to the b-axis. Adjacent layer A units are separated by the NaCl-type building unit, layer B, that consists of octahedrally coordinated M2 and M4 metal positions. The pnictogen atoms (Pn) predominantly occupy the M1 and M2 sites, whereas the transition metal atoms (M) are located within the M3 and M4 sites. This preferential distribution of Pn and M atoms within the crystal lattice results in the formation of two important substructures. The first is the one dimensional (1D) chain, [MnSe4n+2], of edge-sharing octahedra of Se atoms containing the magnetic transition metal atoms (M3 and M4). This will be thereafter referred to as the magnetic subunit (Figure 1). The second substructure is the three dimensional (3D) framework separating adjacent magnetic subunits and containing the non-magnetic main group metals (M1 and M2) and selenium atoms. This non-magnetic substructure will be referred to as the semiconducting subunit. Our investigation of the MPn2Se4 phases led to the synthesis of several compositions and the study of their structure-properties relationship. This includes the compositions, MnSb2Se4,35 FeSb2Se4,37 Mn1-xSnxBi2Se4,5 FeBi2Se4.30 Fe1-xSnxSb2Se4,32 MnSb2-xSnxSe4,33 and FeSb229

xInxSe4.

We find that the type of magnetic atoms (Mn, Fe) within the [MnSe4n+2] chains

governs the nature of magnetic ordering in these phases. For instance, the FePn2Se4 phases display ferromagnetic behavior,30,37 whereas the Mn analogs (MnPn2Se4) are predominantly antiferromagnetic.5,35 In addition, the type of majority charge carriers in MPn2Se4 is controlled by the nature of the main group element (Bi, Sb) within the semiconducting unit. For example, the conduction type in MPn2Se4 switches from p-type for MSb2Se4 to n-type for MBi2Se4.5,30,35,37 This ability to independently control the magnetic ordering and carrier type within the same crystal lattice through chemical manipulation of the composition of various structural subunits makes it possible to examine the effect of free-carriers on ferromagnetic ordering in spintronic materials. Here, we report on the FeSb2-xSnxSe4 system, a ferromagnetic (FM) analogue of the MnSb2-xSnxSe4 anti-ferromagnetic (AFM) system,33 in which an interesting transition from AFM to FM behavior was observed for an optimum substitution of Sb with Sn (0.05 ≤ x ≤ 0.15).33 The mechanism for ferromagnetic behavior in Sn-doped MnSb2-xSnxSe4 samples can be well understood using the bound magnetic polarons model (BMP).33 The substitution of Sb3+ by Sn2+ in the FeSb2-xSnxSe4 system is anticipated to alter the carrier density while maintaining a constant concentration of the magnetic atoms. This enables a better assessment of the interactions of localized magnetic moments with charge carriers spin and the consequences of these interactions on the electronic and magnetic properties in ferromagnetic semiconductors.

ACS Paragon Plus Environment

6

Page 7 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Polycrystalline powders of FeSb2-xSnxSe4 (0 ≤ x ≤ 0.25) were synthesized by reacting stoichiometric mixtures of the elements at 775 K.33 The synthesized dark gray polycrystalline powder of FeSb2-xSnxSe4 (0 ≤ x ≤ 0.25) was structurally characterized using X-ray powder diffraction and X-ray photoelectron spectroscopy (XPS). Phase purity and thermal stability of various samples was evaluated using Differential scanning calorimetry (DSC) and their electronic and magnetic properties were investigated. Details on the synthesis procedure, structural characterization and measurements of the electronic and magnetic properties were described elsewhere 33 and can also be found in the supporting information. Results and Discussion Crystal structure. Polycrystalline powders of selected FeSb2-xSnxSe4 compositions were synthesized by solid-state reaction of high purity elements in the appropriate ratios. Careful comparison of the powder XRD patterns of all samples with the theoretical pattern calculated using single crystal data of FeSb2Se4 (ICSD 421938)37 indicated an excellent match in peaks positions (Figure 2A) suggesting that the synthesized FeSb2-xSnxSe4 compositions are isostructural with FeSb2Se4 and that all diffraction peaks can be indexed in the space group C2/m (# 12). No obvious extra peak of unknown phase could be detected on the XRD patterns of various compositions up to x = 0.2. Phase purity of all samples was further confirmed through differential scanning calorimetry (DSC), which displayed a single peak of melting with onset temperature between 862 K and 873 K (Figure 3A). The difference between the melting point onset temperature (~10 K for highest value observed) is within the tolerable range to conclude that substituting Sb with Sn in the FeSb2Se4 crystal lattice does not significantly affect the melting temperature (875 K)37 of Fe1-xSb2SnxSe4 system for Sn concentrations x ≤ 0.25. The effect of the substitution between Sb and Sn on the lattice parameters, atomic positions, bond distances as well as the distribution of metal atoms (Fe, Sb, Sn) within the crystal structure of FeSb2-xSnxSe4 compositions was investigated by Rietveld refinement on XRD patterns of the synthesized materials using FullProf.39 The fairly good agreement in both peak position and intensity between the observed and calculated XRD data for all compositions, as shown in Figure 2A as well as in supplementary materials Figure S1, attests to the good quality of the Rietveld refinement results, which is also supported by the relatively low values of the weighted profile R-factors and 2 (Table S1). The refinement result revealed no significant change of the lattice parameters and unit cell volume (Figure 2B and Table S1) upon substituting Sb by Sn, which is surprising since the 22% larger effective ionic radii (IR) of Sn2+ (93 pm) in octahedral

ACS Paragon Plus Environment

7

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 44

coordination compared to Sb3+ (76 pm),41 points to an increase in the lattice parameters and unit cell volumes. Astonishingly, the c parameter slightly contracts with increasing Sn content and the largest drop (c/c = 7%) is observed for the composition x = 0.2 (Figure 2B). In addition, irregular trend in the local distortion of metal coordination environments can be noticed from the variation of both atomic positions and M – Se bond lengths within similar polyhedra in selected compositions (Table S2 and Table S3). Although this unanticipated alteration of the structural parameters can be partially supported by the redistribution of metal atoms within various atomic sites upon increasing Sn content (Figure 2C and 2D), it also suggests a more complex variation in the chemical states (oxidation state) of elements within the crystal structure. The distribution of metal atoms at various positions was refined in the last step of the Rietveld refinement process assuming full occupancy of all atomic positions.

Figure 2: Structural characterization of selected FeSb2-xSnxSe4 compositions. (A) Example of Rietveld refinement result using X-ray diffraction pattern of the sample with x = 0.01; (B) Refined lattice parameters as a function of Sn content; Distribution of Fe (C) and Sn (D) at M1, M2, M3 and M4 sites as a function of Sn content.

ACS Paragon Plus Environment

8

Page 9 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

The refinement results show that the M3(2d) site is exclusively occupied by Fe, regardless of the Sn content, which is consistent with the occupancy of the M3 site in the parent structure, FeSb2Se4.37 The remaining metal sites (M1(4i), M2(4i) and M4(2a)) show mixed occupancy with an interesting redistribution of Fe and Sn in various samples. For instance, starting with the parent structure (x = 0), Sb atoms are preferentially found at M1(Sb:96%/ Fe:4%) and M2(Sb:96%/ Fe:4%) sites, whereas Fe is predominantly at the M4(Sb:16%/ Fe:84%) site.37 Upon substitution of Sb by Sn, the concentration of Fe gradually decreases at the M1 and M2 sites, whereas higher Fe content is found at the M4 sites (Figure 2C). This implies an increased ordering in the distribution of Fe atoms with increasing substitution of Sb by Sn. Full ordering of Fe at M3 and M4 sites is obtained for compositions with x ≥ 0.1. Simultaneously, Sn atoms preferentially occupied the M1 site for compositions with x ≤ 0.04 and started to fill both M1 and M2 sites in samples with x > 0.04 (Figure 2D). This variation in the distribution of Fe and Sn atoms within the crystal structure of samples with increasing Sn content can partly explain the

irregular

alterations

of

the

structural parameters such as lattice constants and bond distances and is also anticipated to have a significantly impact on the magnetic and electronic properties of FeSb2-xSnxSe4 samples. To

further

investigate

the

origin of the unexpected changes observed in the unit cell parameters, we have probed the oxidation state of various cations (Sn, Fe and Sb) in all FeSb2-xSnxSe4 compositions using Xray photoelectron spectroscopy (XPS). The spectra (Figure 3B to 3D) for Sn_3d,

Sb_3d

and

Fe_2p

shell

electrons for various compositions Figure 3: (A) Differential scanning calorimetric heating curves of selected FeSb2-xSnxSe4 compositions revealing a single peak of melting between 862 K and 873 K. X-ray photoelectron spectroscopy spectra of Sn_3d (B), Sb_3d (C), and Fe_2p (D) shell electrons for selected FeSb2xSnxSe4 compositions.

ranging between 0 ≤ x ≤ 0.25 confirmed

the

presence

of

all

elements. The intensity of the Sn_3d5/2 peak at binding energy ~ 486.4 eV,

ACS Paragon Plus Environment

9

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 44

which is associated with Sn2+ in SnSe,42 consistently increases with increasing concentration of Sn in the FeSb2-xSnxSe4 structure (Figure 3B), confirming an effective incorporation of Sn atoms within the FeSb2Se4 crystal lattice at 2+ oxidation state. The Fe_2p3/2 peaks are observed for all samples at binding energy of 711.3 eV along with a broad satellite peak at 716.5 eV for compositions with x ≥ 0.15 (Figure 3D). This is consistent with binding energy associated with Fe2+ in FeO.43-45 The Sb_3d5/2 peak appears at ~ 530.5 eV in the sample with x = 0 and shifts to higher values with increasing concentration of Sn. The observed binding energy for the Sb_3d5/2 peak in FeSb2Se4 is consistent with the reported value in MnSb2-xSnxSe4,33 Sb2Se35,33 or Sb2O3 (530.2 eV)46, implying a +3 oxidation state for Sb atom. Interestingly, the Sb_3d5/2 peak shape broadens with increasing concentration of Sn from x = 0.05 up to 0.25 and the binding energy shifts from ~ 530.5 eV to ~ 530.9 eV. Such peak shifting suggests that a fraction of Sb3+ atoms has been oxidized to Sb5+ upon increasing Sn concentration. The suggested partial oxidation is triggered by the local charge imbalance (excess negative charge) resulting from the substitution of Sb3+ by Sn2+ and the intermixing of these atoms at M1 and M2 positions in the crystal lattice. Such partial oxidation of the larger Sb3+ (IR: 76 pm) to the slightly smaller Sb5+ (IR: 60 pm) compensates for the anticipated increase in the lattice parameter due to the substitution of small Sb3+ (IR: 76 pm) by a larger Sn2+ (IR: 93 pm) leading to nearly constant lattice parameters for samples with x ≤ 0.15. For samples with higher Sn content (x ≥ 0.15), a significantly large fraction of Sb3+(4d105s25p0), which contains stereoactive electron lone-pair is oxidized to Sb5+(4d105s05p0), which has no electron lone-pair, diminishing the impact of the stereoactivity of the Sb lone-pair on the lattice parameter. The active lone-pair, which is oriented along the caxis, hampered lattice contraction along that direction for samples with x ≤ 0.15 and the removal of such stereoctivity through oxidation of Sb3+ to Sb5+ explains well the observed sharp contraction of the c-axis for samples with x ≥ 0.2. A similar effect was observed in MnSb233

xSnxSe4.

Remarkably, the XPS spectrum of Sb_3d5/2 shell electrons for samples with x ≥ 0.2

also showed a shoulder peak at lower binding energy, 529 eV, which matches with the binding energy of elemental Sb.47 This suggests that for samples with x > 0.15, a partial reduction of Sb3+ to elemental Sb0 also takes place simultaneously with the partial oxidation of Sb3+ to Sb5+. The effect of the substitution of Sb3+ by Sn2+ in FeSb2-xSnxSe4 and the resulting change in the chemical state of Sb from Sb3+ to Sb+5 and Sb0 (for samples with x ≥ 0.2) observed from XPS data on the electronic properties of the resulting substitution compound can be rationalized using the chemical equations (1) through (5). The substitution of Sb3+ by Sn2+ in FeSb2-xSnxSe4 leads to the formation of acceptor states accordingly with equation (1).

ACS Paragon Plus Environment

10

Page 11 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

𝑥𝑆𝑏3 +

xSn2+

.

x(Sn2+)’Sb + xh

(1)

The local excess negative charges resulting from the substitution of Sb3+ by Sn2+ triggers partial oxidation of Sb3+ to Sb5+ (for compositions with x ≤ 0.15) or simultaneous oxidation and reduction of Sb3+ to Sb5+ and elemental Sb0 for samples with higher Sn content. These oxidation and reduction reactions lead to the generation of various concentrations of electrons accordingly to the chemical reactions (2) and (3). y (Sb3+  Sb5+ + 2 e)

(2)

z (Sb3+ +3 e Sb0)

(3)

The global equation for this oxidation/reduction reaction is given by equation (4). (y+z) Sb3+  ySb5+ + z Sb0 + (2y-3z) e

(4)

The complete chemical equation of the substitution and the subsequent oxidation and reduction processes taking place within the FeSb2-xSnxSe4 crystal lattice is given by equation (5). Fe2+ + (2-x)Sb3+ + xSn2+ + 4Se2- Fe2+ + (1-y-z)(2-x)Sb3+ + y(2-x)Sb5+ + z(2-x)Sb0 + x(Sn2+)’Sb

.

+ x h + (2y-3z)(2-x) e

(5) with y = x+[3z(2-x)]/2(2-x)

Therefore, the final chemical formula of the substituted composition is

.

(Fe2+)(Sb3+)[(4-3x)/2- 5z(2-x)/2](Sb5+)[x/2+3z(2-x)/2](Sb0)z(2-x)(Sn2+)x(Se2)4(e)[x-3z(2-x)/2](CB) (h )x(VB) One can derive from the above equations, that the substitution of x mole of Sb3+ by an equivalent amount of Sn2+ within the structure of FeSb2-xSnxSe4 induces the formation of [x/2+3z(2-x)/2] moles of Sb5+ and z(2-x) moles of elemental Sb0 along with the generation of both [x-3z(2-x)/2] moles of electrons and x moles of holes within the FeSb2-xSnxSe4 crystal lattice. Therefore, bipolar electronic conduction is anticipated for various FeSb2-xSnxSe4 compositions along with strong change of the magnetic properties. Electronic transport. To fully characterize the impact of the substitution of Sb by Sn on the electronic transport properties of FeSb2-xSnxSe4 samples, electrical resistivity (Figure 4A and 4B) and thermopower (Figure 4C and 4D) data were measured in the temperature range from 300 K to 10 K. Regardless of the temperature, the electrical resistivity and thermopower of various samples decrease drastically with increasing Sn content, reaching the lowest values for x = 0.15 and slightly increase with further increase in Sn content (x > 0.15). For example, the electrical resistivity of the sample with x = 0 is ~ 200 .m at 300 K and decreases drastically to 8 .m for the sample with x = 0.02. The resistivity further drops to 0.1 .m for samples with x

ACS Paragon Plus Environment

11

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 44

≥ 0.1 (Figure 4A and 4B). Likewise, the thermopower at 300 K is 217 μVK-1 for FeSb2Se4 (x = 0) then decreases to ~150 μVK-1 for x = 0.1. Increasing the Sn content to x = 0.15 leads to a drastic drop of the thermopower to 36 μV/K and further increase in Sn content to x = 0.2 induced a slight increase of the thermopower to 50 μV/K (Figure 4C and 4D). The observed drastic drop in the electrical resistivity and thermopower of various FeSb2-xSnxSe4 samples with a small increase in the Sn content can be associated with the large increase in the carrier density (both electrons and holes) arising from the substitution of Sb3+ by Sn2+ and the accompanying oxidation-reduction processes as discussed above. Indeed, the chemical equation (5) indicates the generation of equal amount of holes and electrons, for compositions with x ≤ 0.15, and an increase in the concentration of both carrier types with Sn content (x value). For instance, the generation of ~8.41020 electrons per cm3 and an equivalent density of acceptor states is anticipated for the composition with x = 0.15. However, for compositions with higher x values (x > 0.15), the presence of a Sb_3d5/2 peak at 529 eV in XPS the data, which is associated with elemental Sb0, suggests the generation of a slightly lower density of electrons ([x-3z(2-x)/2] moles, where z(2-x) is the mole fraction of Sb0 formed) compared to holes (x moles). This analysis is in agreement with the alteration of the electrical resistivity and thermopower, where a decrease in the resistivity and thermopower with increasing Sn content up to x = 0.15 and the slight increase in both electrical resistivity and thermpower for the sample with x = 0.2 is observed. The large density of both electrons and holes generated in various FeSb2-xSnxSe4 samples as well as the observed sharp drop in the measured thermopower with a small increase in Sn content imply a bipolar conduction mechanism. Therefore, the overall electrical resistivity and thermopower measured correspond to the contributions of both holes and electrons to the electronic transport. Here, the total electrical conductivity measured is the sum of electron and hole contributions ( = + + -), whereas the thermopower measured at a given temperature is the weighted sum of the positive thermopower (S+) values from hole transport and the negative thermopower (S-) values from electron transport (ST = [+S+ + -S-]/(+ + -)). All samples show positive overall values of the thermopower within the temperature range from 300 K to 60 K, suggesting holes as the overall majority carriers. The observed irregular trend in the variation of the electrical conductivity and thermopower of various FeSb2-xSnxSe4 with increasing Sn content can be attributed to the alteration in the distribution of Fe and Sn atoms in the crystal structure of various samples.

ACS Paragon Plus Environment

12

Page 13 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Figure 4: Temperature dependence of the electronic properties of selected FeSb2-xSnxSe4 samples. (A) and (B) electrical resistivity highlighting the electronic transition at 62 K from p-type semiconducting behavior to n-type metallic behavior for the sample with x = 0.15 (inset of 4B); (C) and (D) thermopower showing a sign reversal at 62 K for the composition with x = 0.15; (E) thermal conductivity; and (F) figure of merit. One interesting feature of the bipolar electronic transport in various FeSb2-xSnxSe4 samples is the observation of multiple discontinuities and transitions on the temperature dependent electrical resistivity and thermopower curves. For example, Figure 4A and 4B revealed a striking variation in the magnitude of the electrical resistivity of various samples upon cooling from 300 K to 50 K. At 300 K, an electrical resistivity of ~ 200 .m was measured for the sample with x = 0. The electrical resistivity gradually increases to ~ 400 .m at 200 K and then suddenly increases by several orders of magnitude exceeding ~1.4 x 106 .m at 50 K. Temperature dependent electrical resistivity curves of Sn-substituted samples with x =

ACS Paragon Plus Environment

13

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 44

0.02 and 0.1 showed a similar trend. However, the magnitude of the electrical resistivity at a given temperature as well as the temperature at which the sharp increase in the electrical resistivity is observed gradually decreases with increasing Sn content. For instance, a close examination of the temperature dependent electrical resistivity curve of samples with x = 0.02 and 0.1 revealed a discontinuity at 180 K and 60 K, respectively, where a small drop in the electrical resistivity is observed. Likewise, a large increase in the overall thermopower was observed for samples with x = 0, 0.02 and 0.1, upon cooling from 300 K to 60 K. The thermopower for the sample with x = 0 increases from ~ 217 μV/K at 300 K to ~ 400 μV/K at 60 K. Likewise, the thermopower of Sn substituted samples (x = 0.02 and 0.10) rapidly increases upon cooling reaching a maximum of 380 μV/K at 90 K for the sample with x = 0.02 and 300 μV/K at 80 K for the sample with x = 0.1. The observed sudden increase in both the electrical resistivity and thermopower of the samples with x = 0, 0.02 and 0.1 could be associated with large localization of charge carriers upon cooling. However, this picture is not consistent with the drop in the thermopower upon cooling below 80 K. Therefore, the observed temperature dependence of the electronic transport in FeSb2-xSnxSe4 samples may also be associated to the onset of an interesting topological electronic transition as discussed below. The electrical resistivity and thermopower curves for the samples with x ≥ 0.15 surprisingly showed a completely different temperature dependence and an abrupt change with increasing Sn content compared to samples with x ≤ 0.1. For example, the electrical resistivity and thermopower at 100 K drastically drops, respectively, from 49  m and 265 μV/K for the sample with x = 0.1 to 0.27 .m and 17.15 μV/K for the sample with x = 0.15. In addition, the sample with x = 0.15 showed only a marginal increase in the electrical resistivity upon cooling from 300 K (0.08 .m) to 100 K (0.27 .m), whereas the electrical resistivity of the sample with x = 0.1 increases by nearly 600% within the same temperature range (Figure 4A). The abrupt drop in both electrical resistivity and thermopower over a wide temperature range upon small increase in the Sn content from x = 0.1 to x = 0.15 suggests a significant change in the underlying electronic structure of FeSb2-xSnxSe4 samples with x ≥ 0.15. Such alteration of the electronic structure can be probed through measurements of macroscopic physical properties such as temperature dependence electrical resistivity and thermopower. As can be seen from Figure 4B, the electrical resistivity initially increases with decreasing temperature from 0.08 .m at 300 K up to the maximum of 0.41Ω m at peak temperature Tp ~ 62 K and sharply decreases to 0.22 Ω.m with further cooling of the sample to 6.7 K. This anomalous switching of the temperature dependent electrical resistivity from the semiconducting behavior from 300 K to 62 K to the

ACS Paragon Plus Environment

14

Page 15 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

metallic behavior in the temperature range below 62 K is a signature of an electronic topological transition known as Lifshitz transition (LT), which is associated with a change of the Fermi surface (FS) topology due to the variation of electronic band structure or change in the Fermi energy.48 A similar electronic transition was also observed at Tp = 19 K for the sample with x = 0.02 that was successfully measured to very low temperatures (Figure S2). A slightly higher electrical resistivity was observed at low temperatures upon increasing the Sn content to x = 0.2. For this sample, the electrical resistivity gradually increases with decreasing temperatures and plateaus to ~700 .m at 20 K without full reversal to metallic behavior. Further confirmation of the alteration in the electronic structure of FeSb2-xSnxSe4 samples with x ≥ 0.15 is obtained through the temperature dependent measurement of the thermopower. Both samples (x = 0.15 and 0.2) initially show a slow decrease in thermopower with decreasing temperature up to a critical temperature (Tc) (from 36 μV/K at 300 K to 30 μV/K at 150 K for x = 0.15 and from 50 μV/K at 300 K to 40 μV/K at 90 K for x = 0.2) beyond which a sharp drop is observed (Figure 4D). This is a surprising trend given that the electrical resistivity,

within

the

same

temperature range, increases with decreasing

temperature.

Interestingly, the sample with x = 0.15 shows a sign reversal of the thermopower at 62 K, which also corresponds to the temperature of maximum

electrical

resistivity

curve (Figure 4B). Beyond this temperature, negative values of the thermopower

(n-type)

were

observed (-2 μV/K at 40 K), which

Figure 5: Schematic illustration of the electronic band structure of FeSb2-xSnxSe4 samples. (A) The substitution of Sb3+ by Sn2+ leads to the formation of two subbands within the band gap, a lower band (LB) of acceptor states and an upper band (UB) of donor states; (B) proposed down shift, upon cooling, of the overall band structure that explains the change of the Fermi surface topology from hole-like to electron-like (Lifshitz transition).

implies a switching of the majority carriers type from hole to electron. Astonishingly, the change in the majority carriers type vanishes upon increasing the Sn content to x = 0.2, as the thermopower remains positive down to 10 K (10.7 μV/K). This further points to the high sensitivity of the electronic structure of FeSb2-xSnxSe4 samples to doping. The observed larger electrical resistivity and thermopower for

ACS Paragon Plus Environment

15

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 44

the sample with x = 0.2 compared to the samples with x = 0.15 is consistent with the chemical equation (5) that predicts the generation of fewer electrons due to the formation of elemental Sb0 accordingly with the XPS data. Origin of Lifshitz transition. The simultaneous change of the conduction behavior from semiconducting to metallic along with the rapid decrease in the thermopower with decreasing temperature and the reversal of the sign of majority charge carriers type from hole to electron observed at 62 K for the sample with x = 0.15 and at 19 K for the sample with x = 0.02 (Figure S2) strongly points toward a rapid change of the Fermi surface topology with doping as well as temperature, which is a manifestation of Lifshitz transition.48 Similar type of unusual transport properties (strong resistivity peak along with a sign reversal of majority carriers) were reported in iron pnictides,49 transition metal pentatellurides (ZrTe5, HfTe5 etc),28,50 high temperature superconductors,51,52 and strongly correlated electron systems such as NaxCoO2.53 However, the precise origin of such anomalous transport properties remains a subject of a longstanding debate. Several mechanisms have been proposed for rational understanding of the observed transport anomalies. These include, bipolar conduction,54 structural phase transition,55 formation of charge density wave,56 polaronic behavior,57 semimetal-semiconductor transition,58 or temperature induced band shift.59-61 Building on current understanding of Lifshitz transition in semiconductors, the anomalous transport properties observed for various FeSb2-xSnxSe4 samples can be rationalized within the context of bipolar conduction coupled with a down-shift of the overall electronic structure with decreasing temperature as detailed below. The substitution of Sb3+ by Sn2+ leads to the creation of acceptor states (extra holes), which are located within the band gap and just above the valence band (VB). The density of these acceptor states increases with increasing Sn content (x values). At sufficiently, high concentrations, x > 0.1, these acceptor states form a hole-like sub-band, hereafter denoted lower band (LB) within the band gap (Figure 5). Simultaneously, the local charge imbalance resulting from the substitution of Sb3+ by Sn2+ triggers partial oxidation of Sb3+ to Sb5+, which is accompanied by the generation of donor states (extra electrons) that are located within the band gap. These donor states, for sufficiently high x values (x > 0.1), also form a sub-band, denoted upper band (UB), just below the conduction band (CB). The resulting electronic band structure for FeSb2-xSnxSe4 samples with x ≥ 0.15 contains two sub-bands (UB (electron-like) and LB (hole-like)) that are separated by a small energy gap (Figure 5) within the narrow band gap (Eg ~ 0.2 eV) of FeSb2Se4 phase.37 Therefore, the topology of the Fermi surface in the resulting band structure and hence the transport behavior is dictated by the relative position of the sub-bands (UB and LB) with respect

ACS Paragon Plus Environment

16

Page 17 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

to the Fermi energy (EF). It should be noted that the relative position of the sub-bands with respect to the Fermi energy level depends on the doping level (the width of the sub-band increases with x values), but most importantly on the temperature, assuming a temperature dependent evolution of the overall band structure. Within this proposed band picture, the semiconducting behavior of the temperature dependent resistivity as wells as the gradual slow decrease of the thermopower observed from 300 K to 62 K for samples with x ≥ 015 implies that at 300 K, the Fermi energy is located within the energy gap (Figure 5A) between the UB and LB, but closer to the LB (hole-like Fermi surface), which lead to bipolar conduction and an overall p-type semiconducting behavior. However, the UB down-shifts gradually upon cooling, touching the Fermi energy level at 62 K for the sample with x = 0.15 (Figure 5B). At this point the topology of the Fermi surface transitions from hole-like to electron-like, which is manifested on the macroscopic properties by a peak on the resistivity curve and the transition from p-type semiconducting to n-type metallic behavior. Further cooling below 62 K results in the UB crossing the Fermi energy level, leading to the observed n-type metallic conductivity at low temperatures. Although the proposed band picture explains very well the anomalies observed in the temperature dependent electrical resistivity and thermopower of FeSb2-xSnxSe4 samples, a direct visualization of changes in their electronic structure with doping level and temperature using high resolution angle-resolved photoemission spectroscopy (ARPES)62 is highly desired to confirm the proposed band structure picture. Temperature dependent evolution of the overall electronic band structure was observed in WTe2,60 and ZrTe5.50 High-resolution angle-resolved photoemission spectroscopy (ARPES) study of WTe2 and ZrTe5 showed that two branches of bands (upper band (UB) and lower band (LB)) with linear dispersion that are separated by an energy gap dominate electronic conductions in both compounds. It was found that the overall electronic structure shifts downward (lower energy) by as much as 70 meV upon cooling from 255 K to 2 K. The energy gap between the upper band (UB) and lower band (LB) decreases but persist down to very low temperature.50 By analogy, one can assimilate the UB and LB bands observed in both WTe2 and ZrTe5 to the sub-bands form in FeSb2-xSnxSe4 samples by donor and acceptor states, respectively. This implies that Lifshitz transition can potentially be switched “on” and “off” on demand in any narrow band gap semiconductor through the engineering of the electronic band structure near the Fermi level via doping so as to create both hole-like (LB) and electron-like (UB) sub-bands within the band gap. Thermal transport and ZT. To fully evaluate the impact of Sb substitution by Sn on the thermoelectric performance of FeSb2-xSnxSe4 samples, thermal conductivity data were collected

ACS Paragon Plus Environment

17

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 44

on hot pressed pellets of selected compositions within the temperature range from 300 K to 2 K. The total thermal conductivity of the pristine sample (x = 0) is ~ 2 WK-1m-1 at 300 K and rapidly increases with decreasing temperature reaching a maximum of ~ 3.4 WK-1m-1 at 180 K (Figure 4E). Further decrease in temperature resulted in a drastic drop of the thermal conductivity to ~ 0.5 WK-1m-1 at 90 K, which remains constant down to 2 K. Interestingly, the magnitude of the total thermal conductivity as well as the temperature at which the maximum total thermal conductivity is observed decreases rapidly with increasing Sn content (Figure 4E). For instance, the total thermal conductivity at 300 K drops to 1.5 WK-1m-1 for the samples with x = 0.02 and further decreases to 1 WK-1m-1 for the samples with x = 0.1. The total thermal conductivity of the sample with x = 0.02 remains nearly constant at 1.5 WK-1m-1 upon cooling from 300 K to 225 K, then rapidly increase with further cooling to a maximum value of 2.5 WK-1m-1 at 150 K and finally drops to 0.5 WK-1m-1 with further down to 50 K. Further increasing the Sn content to x = 0.1 resulted in a drop of the peak thermal conductivity to 1.5 WK-1m-1 at 125 K. The overall thermoelectric figure of merit (ZT) of FeSb2-xSnxSe4 samples is very low (Figure 4F) due to the extremely large electrical resistivity of the pristine sample and the bipolar conduction resulting from the substitution of Sb by Sn. The very high electrical resistivity of the pristine sample (x = 0) leads to a very low ZT value of ~ 7.0 x10-7 at 300 K despite the large thermopower value (~ 217 μV/K) and low total thermal conductivity (2 WK1m-1).

However, partial substitution of Sb atoms by Sn atoms leads to a significant drop in the

electrical resistivity and total thermal conductivity and moderate reduction of the thermopower owing to bipolar doping. This results in marginal improvement of the overall figure of merit of FeSb2-xSnxSe4 samples with increasing Sn content. Magnetism. The magnetic properties of magnetic semiconductors and correlated electrons materials are very sensitive to the alteration of carrier density as well as the separation between magnetic atoms within the crystal. Therefore, changes in the carrier density of FeSb2-xSnxSe4 samples upon substituting Sb by Sn is anticipated to play a significant role on the observed overall magnetic behavior. To investigate the impact of Sb3+ to Sn2+ substitution on the magnetism in FeSb2-xSnxSe4 solid solution, we have measured the field-cooled (FC) and zero-filed cooled (ZFC) temperature dependent magnetic susceptibility of selected compositions under 100 Oe applied magnetic field (Figure 6 and S3) in the temperature range from 2 K to 300 K. Regardless of the temperature, the FC and ZFC susceptibility increase in magnitude with increasing Sn content (x value) from 5% to 20%. For instance, the FC susceptibility at 300 K is 0.36 emu mol-1Oe-1 for

ACS Paragon Plus Environment

18

Page 19 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

the sample with x = 0 and increases to 0.67 emu mol-1Oe-1, 0.70 emu mol-1Oe-1, 1.01 emu mol1Oe-1,

and 1.42 emu mol-1Oe-1 for compositions with x = 0.05, 0.1, 0.15 and 0.2, respectively.

While the change in the FC susceptibility with increasing Sn content is marginal at 300 K, a strong divergence of the FC susceptibility curves for various compositions is observed upon cooling.

Figure 6: Variation of the magnetic susceptibility of selected FeSb2-xSnxSe4 samples with temperature under zero-field-cooled (ZFC) (A and B) and fieldcooled (FC) (C and D) conditions using an external applied field of 100 Oe.

For example the FC susceptibility at 60 K is 0.51 emu mol-1Oe-1 for the pristine sample (x = 0) and increases to 0.87 emu mol-1Oe-1, 1.70 emu mol-1Oe-1, 2.57 emu mol-1Oe-1, and 4.64 emu mol-1Oe-1 for compositions with x = 0.05, 0.1, 0.15 and 0.2, respectively. This result is quite surprising considering that the concentration of the magnetic atom (Fe) is essentially the same in all compositions, which should lead to comparable magnetic susceptibility values. Therefore, one can speculate that the observed large change in the measured magnetic susceptibility of various samples takes into account contributions from free carriers as well as the change in the atomic ordering of Fe atoms in various samples (Figure 2C). Indeed, the increase in the

ACS Paragon Plus Environment

19

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 44

magnitude of the susceptibility observed on various samples with increasing Sn content is in agreement with the observed drop in the electrical resistivity. This further confirms the effect of increase carrier concentration on the overall magnetic susceptibility. Further discussion of a possible mechanism leading to the observed magnetic property is provided below.

Figure 7: Isothermal field dependent magnetization of selected FeSb2-xSnxSe4 samples collected at 2 K (A and B) and 300 K (C and D) showing ferromagnetism up to 300 K. Note the strong dependence of magnetic parameters, such as the coercivity (Hc) and saturation magnetization (Ms) on the Sn content. The ZFC susceptibility for the samples with x ≤ 0.1 gradually increases with decreasing temperature up to 130 K, then sharply drops within a narrow temperature range between 130 K and 110 K, then continues to decrease slightly down to 40 K beyond which another sharp drop is observed with further cooling (Figure 6B). The observed transition in the susceptibility at 130 K for both samples (x = 0.05 and 0.1) is similar to the one previously reported in FeSb2Se437 and FeSb2-xInxSe429 phases and is in agreement with the similarity in the separation between magnetic centers (Fe atom) within the crystal lattice of FeSb2-xSnxSe4 (x ≤ 0.1) and FeSb2xInxSe4

phases. It is believed that such magnetic transition at 130 K is related to the anisotropic

distortion of the [FeSe6] octahedra upon cooling,37 whereas the sharp drop of the susceptibility below 40 K suggests the presence of AFM (anti-ferromagnetic) ordering within the [FenSe4n+2]

ACS Paragon Plus Environment

20

Page 21 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

magnetic chain that competes with the dominating FM exchange interaction between adjacent magnetic chains in the bulk crystal. Interestingly, the shape of the temperature dependence of the ZFC susceptibility curve drastically changes for samples with Sn content x > 0.1. The drop in the ZFC susceptibility below the magnetic transition at 130 K is weak for the composition with x = 0.15 and nearly fades away for the sample with x = 0.2 (Figure 6A). This is a very significant result that further supports our hypothesis that the stereoactivity of the electron lonepair on Sb atom contributes to the sharp transition in the magnetic susceptibility at 130 K.37 Indeed, as observed from the XPS data (Figure 3B), the fraction of such lone-pairs decreases with increasing Sn content due to the oxidation of a large fraction of Sb3+ to Sb5+. This severely diminishes the effect of the lone-pair, which is manifested by a contraction of the c-axis and a slight increase in the a-axis for samples with x ≥ 0.15. This change in the lattice parameters leads to a more regular shape of the [FeSe6] octahedra. Therefore, the observed fading of the magnetic transition at 130 K for FeSb2-xSnxSe4 samples with x > 0.1 implies that the regular shape of the [FeSe6] octahedra remains upon cooling below 130 K, which inhibits the redistribution of spin within the 3d6 orbitals. The analysis of FC susceptibility curves revealed a large difference in both the shape and magnitude compared to the ZFC curve from the same sample. For example, the FC magnetic susceptibility for compositions with x ≤ 0.10 increases with decreasing temperature down to 130 K where the magnetic transition is observed. While the samples with composition, x = 0 and x = 0.05 showed a drop in the susceptibility as anticipated, the magnitude of such change in the susceptibility is significantly smaller than that of the ZFC. Interestingly, the sample with x = 0.1, instead of the anticipated decrease in the magnetic susceptibility, exhibits a surprising large increase in susceptibility within a narrow temperature range from 130 K to 110 K, then plateaus down to 40 K, below which a rapid drop in the susceptibility is observed (Figure 6D). The drop in the susceptibility below 40 K for the sample with x = 0.1 suggests a competition between the AFM intrachain exchange interaction (FenSe4n+2 magnetic chains) and the dominant interchain FM exchange interaction in FeSb2-xSnxSe4. Remarkably, the magnetic transition below 130 K completely vanishes from the FC susceptibility upon increasing the Sn content to x ≥ 0.15. The FC susceptibility curves of these samples show continuous increase of the susceptibility with decreasing temperature from 300 K down to 2 K (Figure 6C). This is a typical behavior for ferromagnetic semiconductors and suggests a dominant interchain FM exchange interaction in FeSb2-xSnxSe4. In addition, the larger magnitude of the FC susceptibility compared to the ZFC data at the same temperature further suggests a strong contribution from spin-polarized free

ACS Paragon Plus Environment

21

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 44

carriers, which act through the formation of bound magnetic polarons (BMPs) with localized magnetic moments on FenSe4n+2 magnetic chains. The observed alteration in the nature of magnetic transition below 130 K could be related to the variation in the distribution of magnetic atom (Fe) within the crystal structure. As shown in Figure 2C, an ordering of Fe atoms exclusively within the M3 and M4 sites is observed for samples with x ≥ 0.1, whereas a small amount of Fe is also found at M1 and M2 sites for composition with x < 0.1. It derives from the above discussion that the magnetic transition at 130 K is related to the distribution of Fe atoms within the crystal lattice, in addition to the contribution from the stereoactivity of Sb3+ electron lone-pair, whereas the magnitude of the susceptibility is controlled by the Sn content, which in turn determines the concentration of holes and electrons generated within the crystal lattice. Additional characterization of the nature of ferromagnetism in the synthesized FeSb2xSnxSe4

samples was achieved through isothermal field-dependent magnetization measurements

at 2 K and 300 K (Figure 7) spanning the applied external magnetic field between -10 kOe and +10 kOe. The S-shape of the magnetization versus applied flied (M vs H) curves for various samples at 2 K and 300 K suggests the persistence of ferromagnetic character in the samples up to room temperature. Increasing the Sn content induces a marginal drop in the coercivity simultaneously with a remarkable increase in the saturation magnetization leading to steeper magnetization curves for compositions with x ≥ 0.15. For instance, the saturation magnetization at 2 K gradually increases with increasing Sn content, starting from 555 emu mol-1 for the pristine sample (x = 0), reaches a maximum of ~ 1400 emu mol-1 for the sample with x = 0.15, and decreases to ~900 emu mol-1 with further increasing the Sn content to x = 0.2. A similar trend was obtained for the magnetization curves at 300 K. Simultaneously, the coercive field (Hc) at 2 K marginally decreases with increasing Sn content from 180 Oe for x = 0 to 160 Oe for x = 0.15 and further drops to 66 Oe for the sample with x = 0.2. This implies that the large increase in the density of carriers (holes and electrons) arising from the substitution of Sb by Sn gradually transforms the FeSb2-xSnxSe4 samples from “semi-hard” magnetic materials to “soft” magnetic materials, which is characterized by a relatively facile spontaneous magnetization and spin reversal upon application of a weak external magnetic field. The observed trend in the alteration of the saturation magnetization and the coercivity with increasing Sn content is consistent with the change in the electrical resistivity with increasing Sn content, discussed above, where the minimum resistivity is achieved for x = 0.15. This further indicates a strong correlation between electronic transport and magnetic ordering in FeSb2-xSnxSe4 samples. Another interesting result in this study is the surprisingly larger values of the coercivity at 300

ACS Paragon Plus Environment

22

Page 23 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

K for various FeSb2-xSnxSe4 samples compared to the coercivity at 2 K (Figure 6B and 6D). For instance, the sample with x = 0 showed a coercivity Hc ~ 255 Oe at 300 K compared to 180 Oe at 2 K. Similarly, the coercivity for the sample with x = 0.15 increases to 300 Oe at 300 K from 160 Oe at 2K. This anomalous increase in the coercivity with rising temperature suggests that the strength of magnetic coupling is mediated by the interactions of localized magnetic moments with free carriers-spin in FeSb2-xSnxSe4 samples, which inhibits the drop in coercivity with rising temperature as is normally anticipated from thermal spin fluctuations.

Figure 8: High temperature zero-field-cooled (ZFC) and field-cooled (FC) magnetic susceptibility of selected FeSb2-xSnxSe4 samples measured under external field of 100 Oe. Zero-field-cooled (ZFC) data were collected after heating the as-synthesized samples from RT to 600 K under zero-applied external field. To probe the effect of Sb substitution by Sn on the Curie transition temperature of FeSb2-xSnxSe4 samples, high temperature magnetic susceptibility data were collected under field-cooled (FC) and zero-field cooled (ZFC) conditions. ZFC data were collected post thermal treatment consisting of heating and cooling the sample between 300 K and 600 K without any applied external field. The susceptibility measurement was performed during heating from 300

ACS Paragon Plus Environment

23

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 44

K to 600 K under the applied field of 100 Oe. Then, the sample was cooled to 300 K under external applied field of 100 Oe and the FC susceptibility data was measured on heating with the same external applied field. The resulting FC and ZFC data for selected samples are shown in Figure 8. It can be noted from the FC data that all samples maintain ferromagnetic-like character well above 300 K. For example, the FC susceptibility curve for the sample with x = 0 shows a rapid drop with rising temperature from 300 K to 325 K followed by a weaker drop in the susceptibility with further increase in temperature up to 450 K above which a constant value of the susceptibility is measured up to 600 K. This general shape of the FC susceptibility curve with the two magnetic transitions at 325 K and 450 K, respectively is also observed for the sample with x = 0.1 (Figure 8A). However, the magnetic transition at 450 K vanishes with further increase in the Sn content to x ≥ 0.15 and only the transition at 325 K is observed on the FC susceptibility curves. This strong alteration in the general shape and magnitude of the temperature dependent susceptibility can be correlated to the distribution of magnetic metal atom (Fe) at various crystallographic sites within the crystal structure of FeSb2Se437 as well as the magnetic interaction of localized magnetic moments with charge carriers spin. In a previous study on Fe1-ySnySb2Se4 phases, we showed that large separation between magnetic centers, achieved through substitution of Fe by Sn, resulted in a large drop of the Curie temperature, Tc, from 450 K for the composition with y = 0 to 325 K for the structurally ordered sample with y = 0.13.32 In the present work, we found that a similar manipulation of the separation between magnetic centers within the structure of FeSb2-xSnxSe4 samples can be achieved through preferential distribution of Sn at M1 and M2 sites whereas Fe atoms are ordered within M3 and M4 sites (Figure 2C). Therefore, the two magnetic transitions on the temperature dependent FC susceptibility of samples with x ≤ 0.1 can be associated with local magnetic ordering within individual magnetic chains at 325 K through interaction with charge carriers and the formation of magnetic polarons and the subsequent long-range collective ordering at 450 K of overlapping adjacent magnetic polarons (BMPs).32 For samples with x ≥ 0.15, the long-range ordering vanishes due to the larger separation between adjacent magnetic polarons. Therefore, only the intrachain magnetic ordering at 325 K (Figure 8B) is observed. Another surprising magnetic behavior in FeSb2-xSnxSe4 samples is the observation of negative values of the ZFC magnetic susceptibility for compositions with x ≤ 0.1 (Figure 8C and 8D). ZFC susceptibility data were collected after a thermal treatment, under zero applied field, that consisted of heating the samples to 600 K followed by cooling to 300K. ZFC susceptibility data measured during heating under applied field of 100 Oe start from negative values at 300 K (-1.3 emu mol-1Oe-1 for x = 0, and ~ -0.1 emu mol-1Oe-1 for x = 0.1) and rapidly

ACS Paragon Plus Environment

24

Page 25 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

increase with temperature up to 325 K where a change in the slope is observed. The magnitude of the susceptibility continues to increase with further increase in temperature, crosses over to positive values at Tcross ~325 K (for x = 0.1) and 410 K (for x = 0) and plateaus at 0.27 emu mol-1Oe-1 (x = 0) or 0.15 emu mol-1Oe-1 (x = 0.1) in the temperature range from 450 K to 600 K. One plausible origin of the observed negative magnetic susceptibility is the formation of trapped internal magnetic field (TIMF) in the sample during cooling under zero applied field from 600 K. It appears that upon heating the samples to 600 K under no external magnetic field, randomization of magnetic moments from BMPs is achieved. During cooling under no external magnetic field, the alignment of magnetic moments is dictated by the strength of the local magnetic field on the individual magnetic sublattice, layer A and layer B (Figure 1). Since the concentration of the magnetic atom (Fe) within the M3 site in layer A is larger than that of the M4 site in layer B, for FeSb2-xSnxSe4 compositions with x ≤ 0.1 (Figure 4), a larger fraction of randomized magnetic moments will align along the direction dictated by the localized moment on M3 sites (layer A). Therefore the TIMF could originate from the non-zero net magnetic moment resulting from the antiferromagnetic exchange interaction between the two ferromagnetic sublattices corresponding to layer A and layer B (Figure 1) as observed in our previous investigation of the FeSb2-xInxSe4 phases.29 The magnitude and orientation of such TIMF could significantly impact on the sign and magnitude of the zero-field cooled magnetic susceptibility data under various applied external magnetic field. For example, large negative magnetic susceptibility will be measured if the magnitude of the external applied field is smaller than that of the TIMF and its orientation is opposite to that of the TIMF. This analysis is consistent with the observed drop in the magnitude of the negative ZFC magnetic susceptibility upon increasing Sn content from -1.3 emu mol-1 Oe-1 at 300 K (for x = 0) to (-) to -0.1 emu mol-1 Oe-1 at 300 K (for x = 0.1). The TIMF vanishes for FeSb2-xSnxSe4 compositions with x > 0.1, since the full occupancy by Fe of both M3 and M4 sites in layer A and layer B, respectively, leads to similar magnetic moments for both sublattices. Conclusion In summary, we have systematically investigated the impact of partial substitution of Sb by Sn on the magnetic and electronic behavior of the p-type ferromagnetic semiconductor FeSb2xSnxSe4.

We found that the charge imbalance arising from the substitution of Sb3+ by Sn2+

triggers partial oxidation Sb3+ to Sb5+ leading to simultaneous hole and electron doping of the FeSb2-xSnxSe4 matrix. This presumably leads to the formation of two subbands within the band gap of FeSb2-xSnxSe4, the upper band formed by donor states and the lower subband formed by

ACS Paragon Plus Environment

25

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 44

acceptor impurities. Such bipolar doping results in a very interesting temperature-induced quantum electronic transition, which is manifested by the presence of anomalous peak on the temperature dependent resistivity curve at 62 K simultaneously with a reversal of the sign of the majority charge carriers from holes (above 62 K) to electrons (below 62 K) for FeSb2-xSnxSe4 sample with x = 0.15. Such electronic transition is reminiscent of the Lifshitz transition, which is associated with the change in the Fermi surface topology from hole-like to electron-like upon cooling from 300 K to 2 K. It was also found that increasing the Sn content leads to a preferential ordering of Fe atoms at both M3 and M4 metal positions for compositions with x > 0.1. This results in an increase in the separation between magnetic centers leading to a decrease in the Curie temperature from 450 K for samples with x ≤ 0.1 to 325 K for compositions with x > 0.1. This result is very significant as it unambiguously settles the debate on the correlation of the Curie transition temperature, Tc, and the concentration of magnetic centers on one hand or the separation between magnetic centers on the other hand. It is found that the separation between magnetic centers rather than the concentration of magnetic atoms dictates the magnitude of the Curie temperature. Interestingly, the difference in the fraction of magnetic atom, Fe, within the sublattices, layer A and layer B leads to the observation of negative ZFC magnetic susceptibility at high temperatures in FeSb2-xSnxSe4 samples with x ≤ 0.1. ASSOCIATED CONTENT Supporting Information The Supporting Information (Experimental section, Table S1 to S3 and Figure S1 to S3) associated with this article can be found in the online version. AUTHOR INFORMATION Corresponding Author [email protected] Notes The authors declare no competing financial interest. ACKNOWLEDGEMENTS This work was supported by the National Science Foundation under Award # DMR-1561008. PFPP and CU also gratefully acknowledge the financial support from Department of Energy,

ACS Paragon Plus Environment

26

Page 27 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Office of Basic Energy Science under Award # DE-SC-0018941 for electronic and thermal transport measurements. JSL also acknowledges financial support from the National Science Foundation Graduate Research Fellowship under Award # DGE 1256260. Magnetic data was measured on a SQUID magnetometer at the University of Michigan, which was purchased using MRI grants from the NSF (DMR- 1428226 and CHE-104008). References (1) Ohno, H. Making Nonmagnetic Semiconductors Ferromagnetic Science 1998, 281, 951. (2) Coey, J. M. D.; Venkatesan, M.; Fitzgerald, C. B. Donor impurity band exchange in dilute ferromagnetic oxides Nat Mater 2005, 4, 173. (3) Chiba, D.; Fukami, S.; K., S.; N., I.; K., K.; T., O. Electrical Control of the ferromagnetism phase transition in cobalt at room temperature Nat Mater 2011, 10, 853. (4) Dietl, T.; Ohno, H.; Matsukura, F.; Cibert, J.; Ferrand, D. Zener model description of ferromagnetism in zinc-blende magnetic semiconductors Science 2000, 287, 1019. (5) Ranmohotti, K. G. S.; Djieutedjeu, H.; Poudeu, P. F. P. Chemical Manipulation of Magnetic Ordering in Mn1-xSnxBi2Se4 Solid-Solutions J Am Chem Soc 2012, 134, 14033. (6) Wolf, S. A.; Awschalom, D. D.; Buhrman, R. A.; Daughton, J. M.; von Molnar, S.; Roukes, M. L.; Chtchelkanova, A. Y.; Treger, D. M. Spintronics: A spin-based electronics vision for the future Science 2001, 294, 1488. (7) Kioseoglou, G.; Hanbicki, A. T.; Sullivan, J. M.; van 't Erve, O. M. J.; Li, C. H.; Erwin, S. C.; Mallory, R.; Yasar, M.; Petrou, A.; Jonker, B. T. Electrical spin injection from an ntype ferromagnetic semiconductor into a III-V device heterostructure Nature Mater 2004, 3, 799. (8) Chappert, C.; Fert, A.; Van Dau, F. N. The emergence of spin electronics in data storage Nature Mater 2007, 6, 813. (9) Awschalom, D. D.; Flatte, M. E. Challenges for semiconductor spintronics Nat Phys 2007, 3, 153. (10) Dobrowolska, M.; Tivakornsasithorn, K.; Liu, X.; Furdyna, J. K.; Berciu, M.; Yu, K. M.; Walukiewicz, W. Controlling the Curie temperature in (Ga,Mn) As through location of the Fermi level within the impurity band Nature Mater 2012, 11, 444. (11) Richardella, A.; Roushan, P.; Mack, S.; Zhou, B.; Huse, D. A.; Awschalom, D. D.; Yazdani, A. Visualizing Critical Correlations Near the Metal-Insulator Transition in Ga1xMnxAs

Science 2010, 327, 665.

ACS Paragon Plus Environment

27

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 44

(12) Sawicki, M.; Chiba, D.; Korbecka, A.; Nishitani, Y.; Majewski, J. A.; Matsukura, F.; Dietl, T.; Ohno, H. Experimental probing of the interplay between ferromagnetism and localization in (Ga, Mn)As Nat Phys 2010, 6, 22. (13) Dietl, T. A ten-year perspective on dilute magnetic semiconductors and oxides Nature Mater 2010, 9, 965. (14) Ando, K. Seeking Room-Temperature Ferromagnetic Semiconductors Science 2006, 312, 1883. (15) Chattopadhyay, A.; Das Sarma, S.; Millis, A. J. Transition Temperature of Ferromagnetic Semiconductors: A Dynamical Mean Field Study Phys Rev Lett 2001, 87, 227202. (16) Jansen, R. Silicon spintronics Nature Mater 2012, 11, 400. (17) Dietl, T.; Ohno, H. Engineering magnetism in semiconductors Mater Today 2006, 9, 18. (18) Kaminski, A.; Das Sarma, S. Polaron percolation in diluted magnetic semiconductors Phys Rev Lett 2002, 88. (19) Kundaliya, D. C.; Ogale, S. B.; Lofland, S. E.; Dhar, S.; Metting, C. J.; Shinde, S. R.; Ma, Z.; Varughese, B.; Ramanujachary, K. V.; Salamanca-Riba, L.; Venkatesan, T. On the origin of high-temperature ferromagnetism in the low-temperature-processed Mn-Zn-O system Nature Mater 2004, 3, 709. (20) Suski, T.; Igalson, J.; Story, T. Ferromagnetism of (Pb,Sn,Mn)Te under high pressure J Magn Magn Mater 1987, 66, 325. (21) Kanomata, T.; Shirakawa, K.; Kaneko, T. Effect of pressure on the Curie temperature of CdCr2Se4 and HgCr2Se4 J Magn Magn Mater 1986, 54–57, 1297. (22) Baltzer, P. K.; Wojtowicz, P. J.; Robbins, M.; Lopatin, E. Exchange Interactions in Ferromagnetic Chromium Chalcogenide Spinels Phys Rev 1966, 151, 367. (23) Matsuda, T.; Munekata, H. Mechanism of photoexcited precession of magnetization in (Ga,Mn)As on the basis of time-resolved spectroscopy Phys Rev B 2016, 93, 075202. (24) Kobayashi, M.; Muneta, I.; Takeda, Y.; Harada, Y.; Fujimori, A.; Krempasky, J.; Schmitt, T.; Ohya, S.; Tanaka, M.; Oshima, M.; Strocov, V. N. Unveiling the impurity band induced ferromagnetism in the magnetic semiconductor (Ga,Mn)As Phys Rev B 2014, 89, 205204. (25) Edmonds, K. W.; Boguslawski, P.; Wang, K. Y.; Campion, R. P.; Novikov, S. N.; Farley, N. R. S.; Gallagher, B. L.; Foxon, C. T.; Sawicki, M.; Dietl, T.; Nardelli, M. B.; Bernholc, J. Mn interstitial diffusion in (Ga,Mn)As Phys Rev Lett 2004, 92, 037201. (26) Blinowski, J.; Kacman, P. Spin interactions of interstitial Mn ions in ferromagnetic GaMnAs Phys Rev B 2003, 67, 121204.

ACS Paragon Plus Environment

28

Page 29 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

(27) Yu, K. M.; Walukiewicz, W.; Wojtowicz, T.; Kuryliszyn, I.; Liu, X.; Sasaki, Y.; Furdyna, J. K. Effect of the location of Mn sites in ferromagnetic Ga1-xMnxAs on its Curie temperature Phys Rev B 2002, 65, 201303. (28) Casamento, J.; Lopez, J. S.; Moroz, N. A.; Olvera, A.; Djieutedjeu, H.; Page, A.; Uher, C.; Poudeu, P. F. P. Crystal Structure and Thermoelectric Properties of the

7,7L

Lillianite

Homologue Pb6Bi2Se9 Inorg Chem 2017, 56, 261. (29) Moroz, N. A.; Lopez, J. S.; Djieutedjeu, H.; Ranmohotti, K. G. S.; Olvera, A.; Ren, P.; Page, A.; Takas, N. J.; Uher, C.; Poudeu, P. F. P. Indium Preferential Distribution Enables Electronic Engineering of Magnetism in FeSb2-xInxSe4 p-Type High-Tc Ferromagnetic Semiconductors Chem Mater 2016, 28, 8570. (30) Ranmohotti, K. G. S.; Djieutedjeu, H.; Lopez, J.; Page, A.; Haldolaarachchige, N.; Chi, H.; Sahoo, P.; Uher, C.; Young, D.; Poudeu, P. F. P. Coexistence of High-T-c Ferromagnetism and n-Type Electrical Conductivity in FeBi2Se4 J Am Chem Soc 2015, 137, 691. (31) Olvera, A.; Shi, G. S.; Djieutedjeu, H.; Page, A.; Uher, C.; Kioupakis, E.; Poudeu, P. F. P. Pb7Bi4Se13: A Lillianite Homologue with Promising Thermoelectric Properties Inorg Chem 2015, 54, 746. (32) Djieutedjeu, H.; Olvera, A.; Page, A.; Uher, C.; Poudeu, P. F. P. High-T-c Ferromagnetism and Electron Transport in p-Type Fe1-xSnxSb2Se4 Semiconductors Inorg Chem 2015, 54, 10371. (33) Djieutedjeu, H.; Zhou, X. Y.; Chi, H.; Haldolaarachchige, N.; Ranmohotti, K. G. S.; Uher, C.; Young, D.; Poudeu, P. F. P. Donor and acceptor impurity-driven switching of magnetic ordering in MnSb2-xSnxSe4 J Mater Chem C 2014, 2, 6199. (34) Poudeu, P. F. P.; Djieutedjeu, H.; Sahoo, P. Crystal Structure of FePb4Sb6Se14 and its Structural Relationship with FePb3Sb4Se10 Z Anorg Allg Chem 2012, 638, 2549. (35) Djieutedjeu, H.; Makongo, J. P. A.; Rotaru, A.; Palasyuk, A.; Takas, N. J.; Zhou, X. Y.; Ranmohotti, K. G. S.; Spinu, L.; Uher, C.; Poudeu, P. F. P. Crystal Structure, Charge Transport, and Magnetic Properties of MnSb2Se4 Eur J Inorg Chem 2011, 3969. (36) Anglin, C.; Takas, N.; Callejas, J.; Poudeu, P. F. P. Crystal structure and physical properties of the quaternary manganese-bearing pavonite homologue Mn1.34Sn6.66Bi8Se20 J Solid State Chem 2010, 183, 1529. (37) Djieutedjeu, H.; Poudeu, P. F. P.; Takas, N. J.; Makongo, J. P. A.; Rotaru, A.; Ranmohotti, K. G. S.; Anglin, C. J.; Spinu, L.; Wiley, J. B. Structural-Distortion-Driven Cooperative Magnetic and Semiconductor-to-Insulator Transitions in Ferromagnetic FeSb2Se4 Angew Chem Int Edit 2010, 49, 9977.

ACS Paragon Plus Environment

29

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 30 of 44

(38) Poudeu, P. F. P.; Djieutedjeu, H.; Ranmohotti, K. G. S.; Makongo, J. P. A. M.; Takas, N. Geometrical Spin Frustration and Ferromagnetic Ordering in (MnxPb2-x)Pb2Sb4Se10 Inorg Chem 2014, 53, 209. (39) Rodriguez carvajal, J. Recent Advances in Magnetic-Structure Determination by Neutron Powder Diffraction Physica B 1993, 192, 55. (40) Rodriguez-Carvajal, J. Recent Developments of the Program FULLPROF, in Commission on Powder Diffraction (IUCr) Newsletter 2001, 26, 12. (41) Shannon, R. Revised effective ionic radii and systematic studies of interatomic distances in halides and chalcogenides Acta Cryst A 1976, 32, 751. (42) Zhang, W.; Yang, Z.; Liu, J.; Zhang, L.; Hui, Z.; Yu, W.; Qian, Y.; Chen, L.; Liu, X. Room temperature growth of nanocrystalline tin (II) selenide from aqueous solution J cryst growth 2000, 217, 157. (43) Biesinger, M. C.; Payne, B. P.; Grosvenor, A. P.; Lau, L. W.; Gerson, A. R.; Smart, R. S. C. Resolving surface chemical states in XPS analysis of first row transition metals, oxides and hydroxides: Cr, Mn, Fe, Co and Ni Appl Surf Sci 2011, 257, 2717. (44) Grosvenor, A.; Kobe, B.; Biesinger, M.; McIntyre, N. Investigation of multiplet splitting of Fe 2p XPS spectra and bonding in iron compounds Surf Interface Anal 2004, 36, 1564. (45) Graat, P. C. J.; Somers, M. A. J. Simultaneous determination of composition and thickness of thin iron-oxide films from XPS Fe 2p spectra Appl Surf Sci 1996, 100, 36. (46) Kim, D. K.; Lee, J. H.; Heo, Y. W.; Lee, H. Y.; Kim, J. J. Co-doping effect of Zn and Sb in SnO2: Valence stabilization of Sb and expanded solubility limit Ceram Int 2011, 37, 2723. (47) Perez, O. E. L.; Sanchez, M. D.; Teijelo, M. L. Characterization of growth of anodic antimony oxide films by ellipsometry and XPS J Electroanal Chem 2010, 645, 143. (48) Lifshitz, I. M. Anomalies of Electron Characteristics of a Metal in the High Pressure Region Sov Phys Jetp-Ussr 1960, 11, 1130. (49) Liu, C.; Kondo, T.; Fernandes, R. M.; Palczewski, A. D.; Mun, E. D.; Ni, N.; Thaler, A. N.; Bostwick, A.; Rotenberg, E.; Schmalian, J.; Bud'ko, S. L.; Canfield, P. C.; Kaminski, A. Evidence for a Lifshitz transition in electron-doped iron arsenic superconductors at the onset of superconductivity Nat Phys 2010, 6, 419. (50) Zhang, Y.; Wang, C. L.; Yu, L.; Liu, G. D.; Liang, A. J.; Huang, J. W.; Nie, S. M.; Sun, X.; Zhang, Y. X.; Shen, B.; Liu, J.; Weng, H. M.; Zhao, L. X.; Chen, G. F.; Jia, X. W.; Hu, C.; Ding, Y.; Zhao, W. J.; Gao, Q.; Li, C.; He, S. L.; Zhao, L.; Zhang, F. F.; Zhang, S. J.; Yang, F.; Wang, Z. M.; Peng, Q. J.; Dai, X.; Fang, Z.; Xu, Z. Y.; Chen, C. T.; Zhou, X. J.

ACS Paragon Plus Environment

30

Page 31 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Electronic evidence of temperature-induced Lifshitz transition and topological nature in ZrTe5 Nat Commun 2017, 8, 15512. (51) Khan, S. N.; Johnson, D. D. Lifshitz Transition and Chemical Instabilities in Ba1-xKxFe2As2 Superconductors Phys Rev Lett 2014, 112, 156401. (52) LeBoeuf, D.; Doiron-Leyraud, N.; Vignolle, B.; Sutherland, M.; Ramshaw, B. J.; Levallois, J.; Daou, R.; Laliberte, F.; Cyr-Choiniere, O.; Chang, J.; Jo, Y. J.; Balicas, L.; Liang, R. X.; Bonn, D. A.; Hardy, W. N.; Proust, C.; Taillefer, L. Lifshitz critical point in the cuprate superconductor YBa2Cu3Oy from high-field Hall effect measurements Phys Rev B 2011, 83, 054506. (53) Okamoto, Y.; Nishio, A.; Hiroi, Z. Discontinuous Lifshitz transition achieved by bandfilling control in NaxCoO2 Phys Rev B 2010, 81, 121102. (54) Shahi, P.; Singh, D. J.; Sun, J. P.; Zhao, L. X.; Chen, G. F.; Lv, Y. Y.; Li, J.; Yan, J. Q.; Mandrus, D. G.; Cheng, J. G. Bipolar Conduction as the Possible Origin of the Electronic Transition in Pentatellurides: Metallic vs Semiconducting Behavior Phys Rev X 2018, 8, 021055. (55) Okada, S.; Sambongi, T.; Ido, M. Giant Resistivity Anomaly in ZrTe5 J Phys Soc Jpn 1980, 49, 839. (56) Okada, S.; Sambongi, T.; Ido, M.; Tazuke, Y.; Aoki, R.; Fujita, O. Negative Evidences for Charge Spin-Density Wave in ZrTe5 J Phys Soc Jpn 1982, 51, 460. (57) Rubinstein, M. HfTe5 and ZrTe5- Possible polaronic conductors Phys Rev B 1999, 60, 1627. (58) McIlroy, N.; Moore, S.; Zhang, D. Q.; Wharton, J.; Kempton, B.; Littleton, R.; Wilson, M.; Tritt, T. M.; Olson, C. G. Observation of a semimetal-semiconductor phase transition in the intermetallic ZrTe5 J Phys-Condens Mat 2004, 16, L359. (59) Manzoni, G.; Sterzi, A.; Crepaldi, A.; Diego, M.; Cilento, F.; Zacchigna, M.; Bugnon, P.; Berger, H.; Magrez, A.; Grioni, M.; Parmigiani, F. Ultrafast Optical Control of the Electronic Properties of ZrTe5 Phys Rev Lett 2015, 115, 207402. (60) Wu, Y.; Jo, N. H.; Ochi, M.; Huang, L.; Mou, D. X.; Bud'ko, S. L.; Canfield, P. C.; Trivedi, N.; Arita, R.; Kaminski, A. Temperature-Induced Lifshitz Transition in WTe2 Phys Rev Lett 2015, 115, 166602. (61) Chi, H.; Zhang, C.; Gu, G. D.; Kharzeev, D. E.; Dai, X.; Li, Q. Lifshitz transition mediated electronic transport anomaly in bulk ZrTe5 New J Phys 2017, 19, 015005. (62) Damascelli, A.; Hussain, Z.; Shen, Z. X. Angle-resolved photoemission studies of the cuprate superconductors Rev Mod Phys 2003, 75, 473.

ACS Paragon Plus Environment

31

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 44

Table of content summary Bipolar transport in p-type FeSb2-xSnxSe4 semiconductors triggered by partial substitution of Sb3+ by Sn2+ leads to Lifshitz-like electronic transition at 62 K for the composition with x = 0.15.

ACS Paragon Plus Environment

32

Page 33 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Table of Content Graphic 118x108mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 1 161x124mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 34 of 44

Page 35 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Figure 2 212x165mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3 134x184mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 36 of 44

Page 37 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Figure 4 164x185mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 5 183x102mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 38 of 44

Page 39 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Figure 6 214x184mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 7 209x164mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 40 of 44

Page 41 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Figure 8 214x179mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure S1 209x174mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 42 of 44

Page 43 of 44 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of the American Chemical Society

Figure S2 249x184mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of the American Chemical Society 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure S3 158x135mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 44 of 44