Emptying Dirac valleys in bismuth using high magnetic fields

# Emptying Dirac valleys in bismuth using high magnetic fields

Zengwei Zhu, Jinhua Wang, Huakun Zuo, Benoît Fauqué, Ross D. McDonald, Yuki Fuseya, and Kamran Behnia (1)Wuhan National High Magnetic Field Center and School of Physics
Huazhong University of Science and Technology, Wuhan 430074, China
(2)MS-E536, NHMFL,Los Alamos National Laboratory,Los Alamos, New Mexico 87545, USA
(3)Laboratoire Physique et Etude de Matériaux (CNRS-UPMC)
ESPCI Paris, PSL Research University, 75005 Paris, France
(4) JEIP, USR 3573 CNRS, Collège de France, PSL Research University, 11, place Marcelin Berthelot, 75231 Paris Cedex 05, France
(5) Department of Engineering Science, University of Electro-Communications, Chofu, Tokyo 182-8585, Japan
July 19, 2019
###### Abstract

The Fermi surface of elemental bismuth consists of three small rotationally equivalent electron pockets, offering a valley degree of freedom to charge carriers. A relatively small magnetic field can confine electrons to their lowest Landau level. This is the quantum limit attained in other dilute metals upon application of sufficiently strong magnetic field. Here, we report on the observation of another threshold magnetic field never encountered before in any other solid. Above this field, , one or two valleys become totally empty. Drying up a Fermi sea by magnetic field in the Brillouin zone leads to a manyfold enhancement in electric conductance. We trace the origin of the large drop in magnetoresistance across to transfer of carriers between valleys with highly anisotropic mobilities. The non-interacting picture of electrons with field-dependent mobility explains most results. Coulomb interaction may play a role in shaping the fine details.

## I Introduction

Like spin and unlike momentum, the valley degree of freedom is a discrete quantum number for Bloch waves in a crystal, with a potential for storing information. ValleytronicsRycerz () has been initially explored in AlAs-based two-dimensional electron gas Shayegan () and subsequently in a variety of systems such as grapheneXiao2007 (), transition-metal dichalcogenidesXu2014 (), diamondIsberg () and siliconRenard (). Magnetic field can modulate occupation of different valleysGokmen (); Arora () or their contribution to the total conductivityZhu2012 (); Collaudin (); Jo2014 (). In bismuth, the latter effect is visible even at room temperature and in presence of a field as weak as 0.5 TZhu2012 (); Collaudin ().

The large orbital magnetoresistance of bismuth has been known for several decadesKapitza (). Following the recent observation of a large magnetoresistance in WTeAli2014 (), the amplitude of magnetoresistance in compensated semi-metals has attracted renewed attentionAli2014 (); Liang2015 (); Shekhar2015 (); Novak2015 (). These studies have brought to foreground a fundamental yet unanswered questions: What sets the magnitude and field dependence of orbital magnetoresistance in the high-field regime where the scattering time exceeds by far the cyclotron periodSong2015 (); Abrikosov ()?

A longstanding open question in condensed-matter physics is the fate of the three-dimensional electron gas pushed beyond the quantum limitHalperin (). When electrons are confined to their lowest Landau levels, the Coulomb interaction, neglected in the single-particle treatment of electrons, is expected to play a significant role. While many dilute metals are pushed beyond the quantum limit, graphite remains the only 3D system in which electronic instabilitiesFauque2013 (), believed to be charge-density-wave statesYaguchi2009 (), have been detected. A unified picture of this regime and the consequences of multiple valleys are yet to be drawn.

We present a study of magnetoresistance of bismuth up to 90T, and a comprehensive angle-dependence up to 65 T, with implications for our understanding of magnetoresistance in a multi-valley electron system far beyond the quantum limit. Previous studies demonstrated that by lifting the degeneracy of anisotropic valleys, magnetic field affects their contribution to the total conductivityZhu2012 (); Collaudin (); Jo2014 (). When the quantum limit is approached, even the carrier population in each valley becomes differentZhu2011 (). However prior to these measurements magnetic fields high enough to wipe out a valley evaporating its Fermi sea had not been applied. Our study detects a large anisotropic drop in magnetoresistance. Comparing the data with theoretical calculations, we argue that this drop happens because one or two Fermi seas, depending on the orientation of magnetic field, totally dry up. Thus, bismuth becomes the only solid in which the amplitude of the magnetic field to attain 100 percent valley polarization, is known and attainable. The semi-classic and multi-valleyAubrey () picture of non-interacting electrons provides a qualitative explanation of the observation. The large drop in magnetoresistance is a consequence of transfer of carriers from a high-mobility to a low-mobility valley. A quantitative model for the drop of resistance at would require a full consideration of the role of Coulomb interactions.

## Ii Results

As reminded in panels a-c of Fig. 1, bismuth is a compensated semi-metal with electrons at L-point on top of a very small gap and with a quasi-linear dispersion, giving rise to three Dirac valleys, which are elongated Fermi surface ellipsoids located at the boundary of the Brillouin zone (For a review seeFuseya2015a ()). The extreme anisotropy of these ellipsoids make them quasi-two-dimensional (Q2D) metals. In real space, Bloch waves of each three Q2D metal can be easily displaced by electric field along the trigonal axis and one of the three binary axes, but not along the third perpendicular orientation, one of the three bisectrix axes. The rhombohedral structure of the bismuth crystal can be seen as a distortion of two interpenetrating Face-Centered Cubic (FCC) structures (Fig. 1c). For each valley index, the preferential charge flow along a binary axis favors one out of the three bonds linking atoms of one FCC network (in yellow) to two of its nearest neighbors belonging to the other FCC network (in black)(Fig. 1c). According to the present findings, a sufficiently strong magnetic field can completely freeze one or two of these three equivalent channels.

Our principal experimental result is presented in Fig. 1d. For both binary and bisectrix orientations of magnetic field, the magnetoresistance stops to increase and begins to drop when the field exceeds 35 T. When the field is along the binary axis, the drop is more drastic and starts at a slightly higher magnetic field. As we will see below, both these features find a natural explanation if the drop is due to the total evacuation of one or two electron pockets. The angle dependence of the magnetoresistance presented in a polar plot in Fig. 1e clarifies the link between this result and previous studies of angle-dependent magnetoresistance in bismuthZhu2012 (); Collaudin (). As seen in the figure, at 5 T, angular oscillations of magnetoresistance are visible, arising from the anisotropy of the mobility. For each ellipsoid the magnetoresistance is larger when the field is along its longer axis. As the field increases, the relative magnitude of the oscillations are damped because the relative contribution of the hole pocket to the total conductivity is enhanced. At 55 T, however, pronounced angular oscillations appear again with an amplitude much larger than what seen at low fields. Note also the loss of threefold symmetryZhu2012 (); Collaudin (); Kuechler (), which leads to the slight but visible different in for nominally equivalent binary and bisectrix axes. The origin of this field-induced ’nematicty’ is a subject of ongoing research. We also measured magnetoresistance of another bismuth crystal up to 90.5 T with magnetic field oriented along the binary axis and found that the system remains metallic with significant implications for the Landau sepctrum as discusses below.

Before discussing the origin of the drop in magnetoresistance, let us consider the Landau spectrum of bismuth in the present configuration. Our data extends the maps of Landau spectrum in bismuth previously drawn up to 12 TZhu2011 () and 28 TZhu2012a () based on studies of the Nernst effect, to 65 T. Fig. 2 shows the theoretical evolution of the Landau levels, the Fermi energy (a and d), the carrier concentration (b and e) and its distribution among valleys (c and f) with magnetic field (See more about the theoretical details in Supplementary Note 2). With magnetic field perpendicular to the trigonal axis and exceeding 15 T, electrons are confined to the lowest spin-polarized Landau sub-level(0Zhu2011 ()) and holes occupy their three (each doubly degenerate) lowest Landau levels. With increasing magnetic field, three events occur at successive magnetic fields. The 2 and 1 Landau levels cross the Fermi level and last occupied electron LL for one (when bisectrix) or two (when binary) pockets becomes empty. The Fermi energy is significantly affected by a modest magnetic field. This is because the carriers are light (almost a thousandth of the free electron mass along bisectrixZhu2011 ()) and therefore the cyclotron energy is large. Charge neutrality imposes equality between electron and hole concentrations. Since electron degeneracy rises faster than hole degeneracy, the Fermi energy shifts downward with magnetic field up to 20 T. The evacuation of each hole Landau level produces a deceleration kink in this downward shift. Note the non-monotonic field-dependence of the 0 level for different valleys. It depends on the magnitude of the interband coupling and the the g-factor corrections (See more about the theoretical details in Supplementary Note 2) and sets the field at which a valley becomes empty.

As seen in panels b and e, the carrier concentration is drastically modified by magnetic field. It increases first, before saturating and dropping at higher fields. The distribution of carriers among electron pockets (panels c and f) is also strongly affected by magnetic field. It ends up by total evacuation of one or two electron valleys. Thus, magnetic field can weigh enough on the balance of energy to dry a Fermi sea without causing a metal-insulator transition. The total number of remaining electrons is still matching the number of holes respecting charge neutrality and keeping the system metallic.

Fig. 3a presents a color map of the angle-dependent magnetoresistance. The angular evolution of the second derivative is presented in Fig. 3b and 3c. The evacuation of each hole Landau level occurs at a field which does not strongly vary much with azimuthal angle. The Zeeman splitting of the hole Landau levels vanishes when the field is perfectly perpendicular to the trigonal axisFuseya (), but becomes finite in case of slight misalignment. With these features in mind, one can identify the horizontal bright lines of Fig.3b as hole Landau levels and their bifurcation at intermediate angles as the result of finite Zeeman splitting arising from imperfect alignment. The experimental data in panel c is to be compared with the theoretical angle-dependent Landau spectrum of panel d (which takes into account the misalignment(See more about the theoretical details in Supplementary Note 2). The excellent agreement validates parameters used for the simulation. The drop in magnetoresistance occurs when the 0 electron Landau level of a valley crosses the chemical potential, leaving no carriers in the valley.

## Iii Discussion

The large drop in magnetoresistance cannot be due to the evacuation of the 1 hole level, as suggested decades ago by authors observing the onset of this dropMiura1983 (); Miura1994 () and using a theoretical model contradicted by the absence of metal-insulator transition evidenced by Fig.1f (See more about the theoretical details in Supplementary Note 2) for details). This interpretation becomes implausible due to the sharp angular dependence of the drop resolved here. Why is the drop much larger than any other quantum oscillation arising from the evacuation of a hole or electron Landau level (Fig. 1 d, f) and has a different temperature dependence(See the experimental results in Supplementary Note 1)? Why it is more than twice larger and significantly wider when the field is along the binary axis? Such questions find straightforward answers as soon as one considers the consequences of emptying an electron valley by magnetic field.

As seen in Fig. 1d, magnetoresistance in bismuth is approximately linear in field over an extended window of magnetic field up to 25 T. Its angle dependence in presence of moderate magnetic field can be quantitatively explained in a semi-classical approach taking carrier mobility as a tensorCollaudin (). The non-trivial field dependence of magnetoresistance (its deviation from the quadratic behavior expected for a compensated semi-metal) reflects the fact that both carrier density and the components of the mobility tensor vary with magnetic fieldSong2015 (). In this context, a change in total carrier concentration, a change in the occupation of different valleys or a change in the relevant mobility component are three distinct possibilities to cause a change in magnetoresistance. Let us examine each of these possibilities.

As seen in Fig. 2, following the evacuation of 1 Landau level, carrier concentration decreases. This would lead to an enhancement (and not a decrease) in magnetoresistance. Whatever causes the drop in magnetoresistance should counter and outweigh the enhancement due to the decreased carrier number. On the other hand, the change in the fraction of electrons occupying each valley (shown in panels c and f of Fig. 2) would generate a drop in magnetoresistance according to what we know of the mobility anisotropy of each valley. When the field is parallel to the longer axis of one electron ellipsoid, the relevant component of the mobility tensor for electrons in this ellipsoid is large. When this valley is emptied and its electrons move to the other two, the relevant component of the mobility tensor becomes lower and therefore magnetoresistance drops. When the magnetic field is along the binary axis, two valleys become empty with consequences more drastic than emptying one, hence a larger drop in magnetoresistance. As seen in panels c and f of Fig. 2, the field window for carrier transfer between valleys is remarkably close to the width of the experimentally-observed drop in magnetoresistance.

We can express our interpretation in quantitative terms. In the high-field limit (), the conductivity for charge current applied along the trigonal axis, ,can be approximated in the following way (See the detailed derivation in Supplementary Note 3):

 σbin33≃enhB2(faμa+2fbμb) (1)
 σbis33≃enhB2(faμ′a+2fbμ′b) (2)

The indexes bin and bis refer to orientation of the magnetic field., , and are carrier densities of the hole pocket and the electron pockets , b, and c (For binary and bisectrix, ). are the fraction of carriers residing in a valley. Parameters , , and are effective mobilities, each combinations of the components of the electron and hole mobility tensor(See the detailed analysis in Supplementary Note 3). If we assume that the effective mobilities remain constant, the change in conductivity depends on the variation of , and . The decrease in pulls down conductivity. But since and , the change in and leads to an opposite effect. Because one or two valleys are emptied, and change more drastically than (See Fig. 2). Therefore, the increase in conductivity caused by the transfer of carriers from one valley to another easily outweighs the decrease expected due to the change in total carrier density.

We note that the effect discussed here bears a similarity to the Gunn effectGunn (). In a semiconductor such as GaAs, a sufficiently strong electric field can displace carriers from one band to another causing a change in differential conductivityBlakemore (). The effect we are observing here is also a consequence of carrier transfer between branches with different mobilities caused by magnetic (and not electric) field.

We also measured the variation of the Hall resistivity, , up to 60 T. The results are presented in Fig. 4. As seen in the figure, for both configurations, shows a jump more drastic than the one seen in the longitudinal resistivity, . In a compensated semi-metal, a finite Hall resistivity is expected if the mobility of hole-like and electron-like carriers differ, which is the case of bismuth. In semi-classical picture of multi-valley transportAubrey (), the Hall conductivity is expected to vanish in the high-field limit () when the current is injected along symmetry axes. However, this picture assumes that carriers are not scattered from one valley to another. The finite Hall conductivity observed here indicates the existence of such scattering events. Interestingly, the Hall response becomes vanishingly small above for the binary configuration (and not for bisectrix). Only in the former configuration electrons are confined to a single valley. Therefore, the magnitude of at 60T for the two configurations is in qualitative agreement with our picture.

A quantitative explanation of the magnitude of the drops in and requires an accurate knowledge of the evolution of the components of the mobility tensor across . In the supplement(See the more analysis in Supplementary Note 4), we employ the semi-classical transport theory to fit the experimental data over a wide field window. This narrows down the possible scenarios for the evolution of mobility with magnetic field and leads us to conclude that the non-interacting picture gives a reasonable account of the data. However, the absence of a perfect fit leaves room for a possible role of Coulomb interaction. The sharpness, the amplitude and the anisotropy of the drop in magnetoresistance cannot be reproduced without invoking an abrupt change (by roughly a factor of two) in and , the components of the electron mobility tensor in the vicinity at . This leaves room for speculation on the role played by interaction.

What can cause such a sudden change in mobility in the vicinity of ? We note that there is a correlation between the boundaries of the triangular blue regions of Fig. 3a and the contours of Landau levels in Fig. 3d. As highlighted by a red circle, two electron lines and one hole line meet together where the drop in magnetoresistance occurs and emphasize that such a meeting never occurs for Landau levels with a higher index number. When this meeting occurs( for B 50 T binary), the energy cost of a trion (a charged exciton with two electrons and one hole), which is the kinetic energy of L-point electrons and T-point holes becomes vanishingly small. Trions have been observed in semiconductors subject to frequency-tuned photonsKheng1993 (). The formation of excitons in a semi-metal such as bismuth as a consequence of Coulomb attraction between electrons and holes is a longstanding ideaHalperin (). The angle-dependent Landau spectrum emerging from this study specifies a region where their energy cost becomes very small. For binary, according to our calculations, the difference in energy between electron and hole Landau levels across the Fermi energy becomes as low as 0.5 meV (See Fig. 5). Now, with a mean distance between electrons and holes of 10 nm and a dielectric coefficient of , the exciton binding energy of an electron-hole pair is E= 0.7 meV. At this stage, we can only speculate if it is more than accidental that the experimentally-observed change in magnetoresistance occurs when the apparent energy cost of binding electrons to holes becomes comparable to the estimated magnitude of attraction between them.

## Iv Conclusion

In summary, based on an extensive set of angle-dependent magnetoresistance data, we conclude that the drop in magnetoresistance is not due to the evacuation of a hole Landau level as proposed previouslyMiura1983 (); Miura1994 (), but the consequence of emptying at least one electronic valley. Thus, our study identifies bismuth as the first solid in which cyclotron energy is strong enough to dry a Fermi sea. This happens because of a unique combination of several features. First of all, the system is dilute and therefore the Fermi energy is small. Second, the carriers are light and both the cyclotron and the Zeeman energy are large. Third, unlike graphene, bismuth is not a perfect Dirac system. In a perfect Dirac system, the spin is locked to angular momentum leading to a cancellation of field-induced cyclotron and Zeeman shifts. Therefore, the lowest Landau level does not shift with magnetic field. The non-trivial evolution of the lowest Landau level(See more about the theoretical details in Supplementary Note 2) leads to the total evacuation of a valley in an attainable magnetic field. By shrinking the Fermi surface with Sb dopingFuseya (), one may lower sufficiently to make field-induced hundred-percent valley polarization as accessible as its spin counterpart.

## V Methods

The dimensions of the samples used in this study were 211mm. They were cut by a wire-saw from a big crystal. The RRR was in the range of 30 to 50. Resistance was measured with a standard four-probe technique. A mechanical rotator was used to rotate the sample under pulsed magnetic field reaching 65 T. 25 m-diameter gold wires were attached to four electrodes on the sample using Dupont 4929N silver paint. The angle could be determined from the signal obtained from a coil attached to the back of rotating platform, which monitored the magnitude of the field parallel to its axis. Electric current was applied along the trigonal direction. The results were reproduced in the study of two different samples and was repeated in two different laboratories (NHMFL-Los Alamos and WHMFC-Wuhan).

## Vi Data availability

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

## Vii Acknowledgements

Z. Z. and this work was supported by the 1000 Youth Talents Plan, the National Science Foundation of China (Grant No. 11574097), The National Key Research and Development Program of China (Grant No.2016YFA0401704) and by directors funding grant number 20120772 at LANL. RMcD acknowledges support from the US-DOE BES ’Science of 100T’ program. In France, this work was part of SUPERFIELD and QUANTUM LIMIT projects funded by Agence Nationale de la Recherche. B.F acknowledge support from Jeunes Equipes de l¡¯Institut de Physique du Collége de France (JEIP). K. B. was supported by China High-end foreign expert program, 111 Program and Fonds-ESPCI-Paris. Y. F. was supported by JSPS KAKENHI grants 16K05437, 15KK0155 and 15H02108.

## Viii Contributions

Z.Z. and K.B. designed the research project. Z.Z. and B.F. prepared the samples. Z.Z. J.W. H.Z. R.D.M carried out experiments. Y.F. provided theoretical support and calculations. Z.Z, Y.F. and K.B. wrote the manuscript and all authors commented on the manuscript.

## Ix Competing interests

The authors declare no competing financial interests.

## X Corresponding author

Correspondence to Zengwei Zhu. *)zengwei.zhu@hust.edu.cn

## References

• (1) Rycerz,A., Tworzydło, J. and Beenakker, C. W. J. Valley filter and valley valve in graphene, Nat. Phys. 3, 172-175 (2007)
• (2) Shayegan, M., De Poortere, E. P., Gunawan, O., Shkolnikov, Y. P., Tutuc, E. and Vakili, K. Two-dimensional electrons occupying multiple valleys in AlAs. Phys. Status Solidi B 243, 3629-3642 (2006)
• (3) Xiao, D., Yao, W. and Niu, Q. Valley-Contrasting Physics in Graphene: Magnetic Moment and Topological Transport Phys. Rev. Lett. 99, 236809 (2007)
• (4) Xu, X., Yao, W., Xiao, D., Heinz, T. F. Spin and pseudospins in layered transition metal dichalcogenides. Nat. Phys., 10, 343-350 (2014)
• (5) Isberg, J., Gabrysch, M., Hammersberg, J., Majdi, S., Kovi, K. K. and Twitchen, D. J. Generation, transport and detection of valley-polarized electrons in diamond. Nature Materials 12, 760-764 (2013)
• (6) Renard, V. T., Piot, B. A., Waintal, X., Fleury, G., Cooper, D., Niida, Y., Tregurtha, D., Fujiwara, A., Hirayama, Y. and Takashina, K. Valley polarization assisted spin polarization in two dimensions. Nature Communications 6, 7230-7237 (2015)
• (7) Gokmen, T., Padmanabhan, M., Gunawan, O., Shkolnikov, Y. P., Vakili, K., De Poortere, E. P. and Shayegan, M. Parallel magnetic-field tuning of valley splitting in AlAs two-dimensional electrons. Phys. Rev. B 78, 233306 ( 2008)
• (8) Arora, A. et al. Valley Zeeman Splitting and Valley Polarization of Neutral and Charged Excitons in Monolayer MoTe at High Magnetic Fields, Nano Lett. 16, 3624-3629 (2016).
• (9) Zhu, Z., Collaudin, A., Fauqué, B., Kang, W., and Behnia, K. Field-induced polarization of Dirac valleys in bismuth. Nat. Phys 8, 89-94 (2012)
• (10) Collaudin, A., Fauqué, B., Fuseya, Y., Kang, W. and Behnia, K. Angle Dependence of the Orbital Magnetoresistance in Bismuth. Physical Review X, 5, 021022 (2015)
• (11) Jo, Y. J., Park, J., Lee, G., Eom, Man Jin, Choi, E. S., Shim, Ji Hoon, Kang, W. and Kim, Jun Sung Valley-Polarized Interlayer Conduction of Anisotropic Dirac Fermions in SrMnBi, Phys. Rev. Lett. 113, 156602 (2014)
• (12) Kapitza, P. The Study of the Specific Resistance of Bismuth Crystals and Its Change in Strong Magnetic Fields and Some Allied Problems, Proc. R. Soc. A 119, 358-443 (1928)
• (13) Ali, M. N., Xiong, J., Flynn, S., Tao, J., Gibson, Q. D., Schoop, L. M., Liang, T., Haldolaarachchige, N., Hirschberger, M., Ong, N. P. and Cava, R. J. Large non-saturating magnetoresistance in WTe, Nature 514, 205-208 (2014)
• (14) Liang, T., Gibson, Q., Ali, M. N., Liu, M., Cava, R. J., Ong N. P., Ultrahigh mobility and giant magnetoresistance in CdAs: protection from backscattering in a Dirac semimetal. Nature Mater. 14, 280-284 (2015)
• (15) Shekhar, C. et al. Extremely large magnetoresistance and ultrahigh mobility in the topological Weyl semimetal candidate NbP, Nat. Phys 11, 645-649 (2015)
• (16) Novak, M., Sasaki, S., Segawa, K., and Ando, Y. Large linear magnetoresistance in the Dirac semimetal TlBiSSe. Phys. Rev. B 91, 041203(R) (2015).
• (17) Song, J. C. W., Refael, G. and Lee, P. A. Linear magnetoresistance in metals: Guiding center diffusion in a smooth random potential. Phys. Rev. B 92, 180204(R)(2015)
• (18) Abrikosov, A. A. Quantum magnetoresistance. Phys. Rev. B 58, 2788 (1998)
• (19) Halperin, B. I. Possible states for a three-dimensional electron gas in a strong magnetic field, Jpn. J. Appl. Phys. Suppl. 26, 1913-1919 (1987)
• (20) Fauqué, B., LeBoeuf, D., Vignolle, B., Nardone, M., Proust, C. and Behnia, K. Two Phase Transitions Induced by a Magnetic Field in Graphite. Phys. Rev. Lett., 110, 266601 (2013)
• (21) Yaguchi, H. and Singleton, J. A high-magnetic-field-induced density-wave state in graphite. J. Phys. Condens. Matter 21, 344207 (2009)
• (22) Fuseya, Y., Ogata, M. and Fukuyama, H. Transport Properties and Diamagnetism of Dirac Electrons in Bismuth. J. Phys. Soc. Jpn. 84, 012001 (2015)
• (23) Zhu, Z., Fauqué, B., Fuseya, Y. and Behnia, K. Angle-resolved Landau spectrum of electrons and holes in bismuth. Phys. Rev. B 84, 115137 (2011)
• (24) Aubrey, J. E. Magnetoconductivity tensor for semimetals, J. Phys. F 1, 493-497 (1971)
• (25) Küchler, R., Steinke, L., Daou, R., Brando, M., Behnia, K. and Steglich, F. Thermodynamic evidence for valley-dependent density of states in bulk bismuth, Nature Mater. 13, 461-465 (2014)
• (26) Zhu, Z., Fauqué, B., Malone, L., Antunes, A. B., Fuseya Y. and Behnia, K. Landau spectrum and twin boundaries of bismuth in the extreme quantum limit, PNAS 109, 14813-14818 (2012)
• (27) Fuseya, Y., Zhu, Z., Fauqué, B., Kang, W., Lenoir, B. and Behnia, K. Origin of the Large Anisotropic g Factor of Holes in Bismuth. Phys. Rev. Lett. 115, 216401 (2015)
• (28) Hiruma, K. and Miura, M. Magnetoresistance Study of Bi and BiSb Alloys in High Magnetic Fields. II. Landau Levels and Semimetal-Semiconductor Transition. J. Phys. Soc. Jpn. 52, 2118-2127 (1983)
• (29) Miura, N., Clark, R. G., Newbury, R., Starrett, R. P. and Skougarevsky, A. V. Anomalous magnetoresistance and Shubnikov-de Haas effect observed for bismuth in pulsed magnetic-fields up to 50 T at temperatures down to 0.3 K. Physica B 194, 1191-1192 (1994)
• (30) Hartman, R. Temperature dependence of the low-field galvanomagnetic coefficients of bismuth. Phys. Rev. 181, 1070-1086 (1969)
• (31) Gunn, J. B. Microwave oscillations of current in III–V semiconductors. Solid State Commun. 88, 883-886 (1963)
• (32) Blakemore, J. S. Semiconducting and other major properties of gallium arsenide. J. Appl. Phys. 53, R123-R181 (1982)
• (33) Kheng, K., Cox, R. T., d¡¯ Aubigné, Merle Y., Bassani, F., Saminadayar, K., and Tatarenko, S. Observation of negatively charged excitons X in semiconductor quantum wells. Phys. Rev. Lett. 71, 1752-1755 (1993)
• (34) Halperin, B. I. and Rice, T. M. Possible anomalies at a semimetal-semiconductor transition, Rev. Mod. Phys. 40, 755-766 (1968)
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters