Abstract
Bulk bismuth has a complex Landau spectrum. The small effective masses and the large g-factors are anisotropic. The chemical potential drifts at high magnetic fields. Moreover, twin boundaries further complexify the interpretation of the data by producing extra anomalies in the extreme quantum limit. Here, we present a study of angle dependence of magnetoresistance up to 65 T in bismuth complemented with Nernst, ultrasound, and magneto-optic data. All observed anomalies can be explained in a single-particle picture of a sample consisting of two twinned crystals tilted by 108° and with two adjacent crystals keeping their own chemical potentials despite a shift between chemical potentials as large as 68 meV at 65 T. This implies an energy barrier between adjacent twinned crystals reminiscent of a metal- semiconductor Schottky barrier or a p-n junction. We argue that this barrier is built by accumulating charge carriers of opposite signs across a twin boundary.
Similar content being viewed by others
Introduction
Bismuth occupies a distinct place in the history of solid state physics (For reviews, see1,2,3,4). Its tiny Fermi surface was the first to be experimentally mapped out thanks to the easily detectable quantum oscillations it generates5,6,7,8,9,10. Recent attention focuses on finding topologically non-trivial states both for the surfaces and the hinges of bismuth crystals11,12,13,14,15,16,17,18,19,20. Recent transport studies on macroscopic crystals have revealed a non-trivial variation of the magnetoresistance21 and the Seebeck effect22 with the amplitude and the orientation of the magnetic field. A robust boundary conductance emerges at quantizing magnetic field, which is yet to be understood23. At the atomic scale, a bismuth layer on a NbSe2 substrate was found to be a platform for detecting nanometric Turing patterns24.
The rhombohedral structure of crystalline bismuth at ambient pressure is the outcome of lowering the cubic symmetry by pulling the cube along its body diagonal axis. Since this can be achieved along each of the four body diagonals, twin formation is common in bismuth crystals25,26 and twin boundaries in bismuth crystals were directly observed decades ago with a scanning tunneling microscope27,28.
The uncontrolled presence of such twin boundaries misled researchers seeking the signatures of collective electronic effects in bismuth. In the extreme quantum limit, when all carriers are confined to their lower Landau levels, a variety of electronic instabilities are conceivable29,30. Experiments seeking such instabilities31,32 explored bismuth beyond its quantum limit, achievable by applying a magnetic field exceeding 9 T along the trigonal axis, and found anomalies which, at first sight, could not be explained in the single-particle picture. Taken at face value, they were attributed to instabilities caused by electronic interactions. However, these interpretations in the case of bismuth were rejected subsequently. The computed Landau spectrum33,34,35,36 was found to be in agreement with the experimental data and the hysteresis attributed to a first-order phase transition32 was not reproduced by torque magnetometry37 or magnetostriction38 measurements.
It was later found that the additional anomalies31 were caused by the presence of a secondary twinned crystal tilted by 108° with respect to the main crystal36. Thus, up to 28 T, the single-particle picture proved to be sufficient to explain the experimental data36. The same theoretical frame was then employed to explain how the evacuation of one or two electron pockets by a 40 T magnetic field, applied along the binary or the bisectrix axes, generates a sudden drop in magnetoresistance39 and an anomaly in magnetization40.
Theory has predicted that when carriers are confined to their lowest Landau levels, the three-dimensional electron gas can suffer a variety of instabilities29,30. The success of the single-particle picture in bismuth pushed beyond the quantum limit is to be contrasted with the case of graphite, another low-carrier semi-metal, in which the quantum limit is accessible in a three-dimensional context. A succession of field-induced phase transitions have been detected in graphite41,42,43,44 indicating the capacity of strong Coulomb interaction to generate collective states29,30 in this material. The two elemental semi-metals differ in several aspects. Screening is weaker in graphite because of its lower dielectric constant. Spin-orbit coupling is much stronger in bismuth. Finally, graphite is much more anisotropic than bismuth, and its electron and hole Fermi pocketsâ longer axis are aligned.
Three-dimensional semimetals have been found to display numerous interesting properties in presence of strong magnetic field. Let us cite two examples. In ZrTe5, for example, there is a magnetic freeze-out beyond the quantum limit and spin-polarised electrons produce an anomalous Hall effect45. In micro-structured Cd3As2, current beams can be guided by tilting a magnetic field46.
In the case of bismuth-antimony alloys, possible signatures of field-induced electronic instability were recently observed47. They are to be distinguished from high-field anomalies due to the presence of twin boundaries seen before48.
Here, we present an extensive study of the Shubnikov-de Haas effect, Nernst effect, ultrasound, and magneto-optics in bismuth at low temperatures and a magnetic field as strong as 65 T. We map the angle dependence of the Landau spectra for the whole solid angle and find that all anomalies can be explained in the single-particle picture, invoking the presence of a secondary twinned crystal. This includes the anomaly observed around 38 T when the magnetic field is applied along the trigonal axis of the main crystal49, whose origin has not been hitherto identified. Thus, up to 65 T, there is no signature of instability caused by electron-electron interaction in bismuth. However, the quantitative agreement between theory and experiment is obtained thanks to a crucial assumption: that the two adjacent twinned crystals each have their own chemical potential μ, set by the orientation of the magnetic field respective to their crystalline axes. We will show that without this assumption, there is no agreement between the computed and the measured Landau spectrum. This unavoidable assumption leads us to conclude that electrons residing in two adjacent crystals do not mix up, despite a 68 meV difference in their chemical potentials at 65 T. We argue that this hypothesis points to the existence of a type of p-n junction at the twin boundary of this semi-metal.
Results
Angle-dependent magnetoresistance in three planes
Figure 1a shows the sketch of the Fermi surface of bismuth, which contains a hole pocket at T point and three cigarette-like electron pockets at L points. Figure 1bâd shows magnetoresistance along three planes, namely binary-trigonal (bin.-tri.), bisectrix-trigonal (bis.-tri.), and binary-bisectrix (bin.-bis.) at 1.6 K. Magnetoresistance shows oscillations which correspond to evacuations of the Landau levels with magnetic field. The different orientations for these planes are indicated. When the magnetic field is along the trigonal axis, MR shows a minimum around 9 T, which corresponds to the evacuation of the penultimate Landau level of hole-like carriers50. As the magnetic field is tilted off the trigonal axis, the cross-section of the hole Fermi pocket perpendicular to the magnetic field increases and this minimum shifts to higher fields. In the binary-bisectrix of Fig. 1d, drops around 40 T are due to the emptying of Dirac pockets39.
Figure 2 presents the second derivative of our magnetoresistance next to the Nernst measurements (Sxy) for three high-symmetric planes36. The two types of data are consistent with each other, but the sensitivity of the Nernst results is relatively higher. Each peak in the Nernst response or in the second derivative of magnetoresistance corresponds to an intersection between an electron (or a hole) Landau level and the chemical potential. Note that the data have been normalized to the largest amplitude of each data set. Figure 2aâc shows the data for the magnetic field rotating in the binary-trigonal, bisectrix-trigonal, and in the binary-bisectrix planes, respectively. In the corner of the panel, the schematic diagram shows the measurement configurations for MR. The quantum oscillations pattern is relatively simple and easily understood when the field is along the trigonal axis51. The hole Landau peaks can be easily identified as their amplitudes are the largest, except for a sudden jump around 40 T along the trigonal axis.
Consistency of measurements
We verified the consistency between diverse techniques when the field was rotated in the bisectrix-trigonal plane and the data obtained on different samples. Figure 3a shows a comparison of the Nernst data collected in Grenoble up to 28 T (dashed lines) with the data collected in Tallahassee up to 45 T (solid lines) for several orientations of the magnetic field. The Nernst data shows an anomaly at 40 T, as reported previously49. In various experiments, the peak positions are close to each other. Panels (b) and (c) compare three distinct sets of data obtained in different laboratories at different times and on different samples. The Nernst data, Sxy, the second derivative of magnetoresistance, \(\frac{{d}^{2}\rho }{d{B}^{2}}\), and the change in the sound velocity (Îv), as the ultrasound propagated along the trigonal direction, are shown for two orientations of the magnetic field, when it is oriented almost along the trigonal axis (b) and when it is 15° (c) off the trigonal axis. There is a good consistency between the three experimental probes, indicating similar peak positions. Specifically, there areâ~â40 T anomalies in Sxy, in \(\frac{{d}^{2}\rho }{d{B}^{2}}\), and in Îv. The slight difference in positions at θâââ15° is consistent with the experimental margin in measuring the angle between the magnetic field and the sample in three different experiments. We also notice that the anomalies in sound velocity and ultrasound attenuation occur at slightly different magnetic fields52.
Comparison with theory
The angle-dependent Landau spectrum can be described by a theoretical model that considers electron pockets at the L points of the Brillouin zone in an extended Dirac Hamiltonian and the hole pocket at the T point as carriers with parabolic dispersion and highly anisotropic g-factor36,39,53,54. In order to explain additional and unexpected anomalies, one needs to assume the presence of an additional secondary twinned domain36,39. In the high field limit, when only a single Landau level is occupied, the interband coupling becomes important and plays a significant role in setting the details of the spectrum39.
The twinned sample originates from the distortion by the symmetry-lowering operation due to Peierls transition55. Since the Schmid factor is 0.48 in bismuth56, applying an external load26 and thermal cycling57 could relatively easily induce twins compared to other materials56. Thus, the strain appearing during sample cleaving by a blaze along cleavage basal trigonal (0001) plane in liquid nitrogen could introduce twinned domains on the sample and these twinned domains should mainly be on the surface. This could be a reason that a smaller RRR was found in a sample with a smaller dimension58, since the ratio of the one plane area to the volume is smaller in a smaller sample. When the sample is twinned, besides the existing principal trigonal, a quasi-trigonal axis will be present.
Figure 4 compares the theoretical prediction with experimental peaks for three high symmetrical planes. The experimental peaks depicted by open symbols are compared with the theoretical results, which are presented by lines from both the main and a secondary tilted crystal as a function of log(B) (See the Supplementary Fig. 2 for a linear plot). For clarity, we have omitted the electron peaks below 28 T, which have been identified earlier36, By considering all contributions from both electrons and holes, as well as the main and twin crystal crystals, the agreement between experiment and theory can be further improved (See Supplementary Fig. 3 for all data sets). Beyond 28 T, we included all peaks associated with electrons from both domains. Figure 4a sketches the main and twinned crystals with their Fermi surfaces. The tilt angle between the two crystals is 108°36. The common binary axis is also shown. All theoretical results of the hole and twinned electron Landau level 0eâ are included.
The theoretical model here is fundamentally the same as the one used before35,36,59, except that the tensor, which is the correction to the additional g-factor in the lowest Landau level \(V{\prime}\), was adjusted to fit the experimental Landau spectra better. Equations (14), (22) from the Ref. 35 and (37), (38) from the ref. 59 were used for calculations. \(V{\prime}\) was adjusted to â0.0025 from â0.062539 (See Supplementary Note 4 for a comparison). As seen in Fig. 4, whereas a slight displacement occurs, which may be due to misalignment during measurement, there is a relatively good agreement between the theory and the experiment. Some more peaks may be due to the electron Landau levels, which have not been clearly identified in the figures, see Fig. 4c. We may conclude, therefore, that the Landau spectrum of bismuth corresponds to what is expected in the single-particle theory. No collective effect due to electron-electron interaction is visible up to 65 T. Because of extreme sensitivity, nullification of the Zeeman splitting of holes along the trigonal axis, the experimentally resolved hole peaks are split due to a small residual misalignment in the case of binary-bisectrix plane (see Fig. 4d). The twinned crystal shares one common binary axis with the original crystal, but tilts its trigonal axis by 108° off39. This is demonstrated by Fig. 4d at zero degrees where two sets of hole peaks merge. The additional peaks vanish as the two crystals are in the same configuration regarding the magnetic field and possess an identical Landau spectrum.
The twinning scenario provides an explanation for the previously reported drop in the magnetoresistance for a magnetic field of 38 T along the trigonal axis of the main crystal49. We can identify it as the result of emptying the electron pockets of the twinned domain. As seen in Fig. 4b, c, the magenta dashed line, which represents the intersection between the 0eâ Landau level of the twinned domain and the chemical potential happens to be atâ~â40 T when the magnetic field is along the trigonal axis of the main crystal. Thus, there is no ultra-quantum feature in the main crystal beyond the quantum limit of holes, which occurs at 9 T along the trigonal axis.
Twins observed by magneto-optics
We confirmed the presence of twin boundaries in our bismuth samples with a magneto-optical probe and found that pulsed magnetic fields have no influence on the domain population. The in-situ images are shown in Fig. 5. The crystal was illuminated with polarized light and the image was recorded using a high-speed camera in pulsed magnetic fields up to 30.8 T oriented along the trigonal axis. The field profile is shown in Fig. 5d. Images are taken at the magnetic field at the times represented by red circles. Three typical images taken at 0 T, 30.8 T, and 0 T after the pulse are shown in (Fig. 5aâc), respectively. The binary axis common to the two domains is labeled in Fig. 5a. The twinned domains occupy a large area close to one-third of the surface (Please see Supplementary Note 1 for more similar results on the observation of twin domains in another sample).
Independent or common chemical potentials for main and twinned crystals
Let us now turn our attention to the most striking aspect of our results. Since the charge is conserved, any change in the density of electrons is to be compensated by an identical change in the density of holes. This leads to a shift in the position of the chemical potential with increasing magnetic field7,35,36,39. This field-induced shift strongly depends on the orientation of the applied magnetic field. The main crystal and its twinned domain are tilted by 108 degrees. Therefore, the magnitude of the field-induced chemical potential shift is different for the two twinned crystals. Thus, the Landau spectrum would be very different in two imaginable cases: (a) The two crystals keep their own distinct chemical potential; (b) The two crystals share the same chemical potential. Figure 6a, b compares the experimental data on twinned hole peaks and the theoretically expected Landau spectrum in the two cases. As seen in Fig. 6b, assuming that the two twinned crystals share a common chemical potential set by the main crystal, the theoretical expectation is incompatible with the experimental results, in sharp contrast with the agreement seen between theory and experiment in the former case, see Fig. 6a.
This conclusion implies a discontinuity in chemical potential and in carrier density across the twin boundary, which is reminiscent of a Schottky barrier. Remarkably, when the magnetic field is along the trigonal axis, the difference in chemical potentials increases up to 68 meV at 65 T (see Fig. 6c). This is remarkable given that the twin boundary width is as narrow as 1ânm28,60. Below, we turn to the underlying factors that impedes electrons belonging to adjacent crystals from mixing up.
Discussion
We now propose a possible scenario explaining this effect. It is well known that the bands are bent near the interface between electronically different solids61. The typical examples are shown in Fig. 7a for the p-n junction and (b) for the Schottky junction. The former is a junction between semiconductors, and the latter is a junction between a metal and a semiconductor. In each case, chemical potentials are different in each side when separated (upper panels of Fig. 7). When they are joined, the carriers are forced to move to make μ common. For example, holes move from p-type to n-type, and electrons move from n-type to p-type in the p-n junction. The n-type band is shifted upward and the p-type band downward by the potential gradient due to the carrier transfer across the interface, resulting in the band bending with a common μ as shown in (Fig. 7a).
A similar band bending is also expected in the present twin boundary between Bi crystals. At zero field, μ is shared between the main and twinned crystals. Under a strong magnetic field, μ of the main crystal shifts differently from that of the twinned one, as shown in Fig. 6c. This situation can be expressed by the model where two electronically different semimetals are joined at the twin boundary as depicted in Fig. 7c.
Suppose that μ of the main crystal is lower than that of the twinned one when two semimetals are separated. When they are joined, the carriers are forced to move to make μ common. As a result, both the electron and hole bands are bent, as if two types of p-n junction coexist. By this band bending, the relative levels between μ, the conduction and the valence bands are kept in each crystal as before joining. The Landau level spectrum of this twinned p-n junction is equivalent to that of the case where the independent chemical potentials are employed. Eventually, the mystery of the apparent jump of μ is solved by this scenario of the twinned p-n junction.
For the twinned p-n junction scenario, we assumed an electric potential difference between the main and twinned crystals. One may question whether an inhomogeneous electric field can exist inside a (semi) metal. The standard answer would be âNo" because finite current will screen the potential difference. However, Landauer argued against this simplistic view62,63. Even in a metal, a point-like scattering center inside a metal will create an unscreened dipolar field, i.e., the electric field is inhomogeneous near the point-like defects. In the present case, the twin boundary, a line or plane defect, will play a similar role to the point-like scattering center of Landauerâs theory.
In this scenario, the Landau eigenvalues curve in the vicinity of a twin boundary in the same way as they bend near the edges of a 2D electron gas64. Any time of such Landau level bending results in edge currents, the continuity equation causes the edge charge and currents to remain close to the edge but distributed in real space in a 2D case64. In the present 3D case, the continuity equation strictly confines the currents and any associated charge variations to the vicinity of the twin boundary. The currents flow along the twin boundaries and essentially account for small differences in magnetization between twins, which should be present due to differences in Landau level filling factors related to the difference in chemical potential. The equilibrium distribution of currents and charge ultimately depends on a balance between the Lorentz force of the current and the twin boundary potential that causes the Landau level eigenstates to curve in real space. These currents flowing along the twin boundary result from a novel form of extended bound state associated with the twin junction. This line of inquiry deserves further investigations.
The twinned p-n junction bears a resemblance to the p-n junction in semiconductors. Therefore, rectification effects observed in semiconductor p-n junctions are also expected in twinned p-n junctions. If these rectification effects are observed, it would provide direct evidence of the potential barrier we propose. However, in the present twinned p-n junction, significant rectification effects are not anticipated to be observed because they are semimetals, not semiconductors. This could explain why direct evidence of the twinned p-n junctions has not been obtained thus far. The rectification effects may become observable if one type of carrier is evaporated, e.g., by applying a gate voltage.
The presence of twins in Bi has already been observed using scanning tunneling microscopy27. Their presence could also be detected using scanning photocurrent microscopy in presence of a magnetic field, though, to our knowledge, there has been no report of such observation yet. In either case, the presence of the potential barrier we propose may be challenging to detect when observing in the semimetallic state. For instance, applying a gate voltage and subsequently scanning after modifying electron-hole compensation in one domain might make detecting the potential barrier possible.
In summary, we have carried out angle-dependence of magnetoresistance along three high symmetric planes, and mapped out their Landau spectra up to 65 T. All anomalies including the one observed atâ~â40 T for a magnetic field oriented along the trigonal axis could be explained by assuming the presence of a second twinned crystal and the two crystals keep separate chemical potentials. We suggested an interpretation of this result by a scenario in which a p-n junction builds up at the twin boundary-areas. This implies the existence of an inhomogeneous electric field in the vicinity of twin boundaries, which can be checked by future experiments.
Methods
Samples and measurements description
Bismuth single crystals commercially obtained used in this study had typical dimensions of 1âÃâ2âÃâ0.5âmm3 and a Room-temperature-to-Residual-resistivity-Ratio (RRR) of \(\left.\rho (300{{{\rm{K}}}})/\rho (4.2{{{\rm{K}}}})\right) \sim 30-100\). Magnetoresistance (MR) measurements in pulsed magnetic field were carried out at Wuhan National High Magnetic Field Center (WHMFC) and in National High Magnetic Field Laboratory- Pulsed Field Facility (NHMFL-PFF) in Los Alamos. The samples were rotated in three distinct planes: binary-trigonal, bisectrix-trigonal and binary-bisectrix. For each set of measurements, the current was applied along the axis perpendicular to the rotating plane. Our high-field MR data complemented the Nernst data obtained below 28 T and reported previously35,36. Additional angle-dependent Nernst measurements up to 45 T were also carried out in the hybrid magnet of the NHMFL DC Facility in Tallahassee. Ultrasound measurements have been collected at Laboratoire National des Champs Magnétique Intense (LNCMI-Toulouse) with pulsed magnetic field up to 54 T. Longitudinal ultrasonic waves were generated using commercial LiNbO3 36° Y-cut transducers glued on a polished surface. A standard pulse-echo technique was used to determine the change in the sound velocity and attenuation. Magneto-optical measurements involving polarized light (for more details, see reference65) were conducted at the International MegaGauss Science Laboratory, the Institute for Solid State Physics at the University of Tokyo. The resulting image was captured using a high-speed camera.
Data availability
The data supporting the present work are available from the corresponding authors upon request.
References
Dresselhaus, M. Electronic properties of the group V semimetals. J. Phys. Chem. Solids 32, Suppl.1, 3 (1971).
Edelâman, V. S. Properties of electrons in bismuth. Soviet Phys. Uspekhi 10, 819 (1977).
Issi, J. P. Low temperature transport properties of the group V semimetals. Aust. J. Phys. 32, 585 (1979).
Fuseya, Y., Ogata, M. & Fukuyama, H. Transport properties and diamagnetism of Dirac electrons in bismuth. J. Phys. Soc. Jpn. 84, 012001 (2015).
Dhillon, J. S. & Shoenberg, D. The de Haas-van Alphen effect III. experiments at fields up to 32 kG. Philos. Trans. R. Soc. Lond. Ser. A, Math. Phys. Sci. 248, 1 (1955).
Shoenberg, D. & Kapitza, P. L. The magnetic properties of bismuth, III. further measurements on the de Haas-van Alphen effect. Proc. R. Soc. Lond. Ser. A. Math. Phys. Sci. 170, 341 (1939).
Smith, G. E., Baraff, G. A. & Rowell, J. M. Effective g factor of electrons and holes in bismuth. Phys. Rev. 135, A1118 (1964).
Bhargava, R. N. de Haas-van Alphen and galvanomagnetic effect in Bi and Bi-Pb alloys. Phys. Rev. 156, 785 (1967).
Liu, Y. & Allen, R. E. Electronic structure of the semimetals Bi and Sb. Phys. Rev. B 52, 1566 (1995).
Kumar, Y., Sharma, P., Karn, N. K. & Awana, V. P. S. Shubnikov-de haas (sdh) oscillation in self-flux grown rhombohedral single-crystalline bismuth. J. Superconduct. Novel Magn. 36, 389 (2023).
Fu, L., Kane, C. L. & Mele, E. J. Topological insulators in three dimensions. Phys. Rev. Lett. 98, 106803 (2007).
Drozdov, I. K. et al. One-dimensional topological edge states of bismuth bilayers. Nat. Phys. 10, 664 (2014).
Benia, H. M., StraÃer, C., Kern, K. & Ast, C. R. Surface band structure of Bi1âxSbx(111). Phys. Rev. B 91, 161406 (2015).
Aguilera, I., Friedrich, C. & Blügel, S. Electronic phase transitions of bismuth under strain from relativistic self-consistent GW calculations. Phys. Rev. B 91, 125129 (2015).
Ito, S. et al. Proving nontrivial topology of pure bismuth by quantum confinement. Phys. Rev. Lett. 117, 236402 (2016).
Murani, A. et al. Ballistic edge states in bismuth nanowires revealed by SQUID interferometry. Nat. Commun. 8, 15941 (2017).
Fuseya, Y. & Fukuyama, H. Analytical solutions for the surface states of Bi1âxSbx (0 â¤xââ²â0.1). J. Phys. Soc. Jpn. 87, 044710 (2018).
Schindler, F. et al. Higher-order topology in bismuth. Nat. Phys. 14, 918 (2018).
Nayak, A. K. et al. Resolving the topological classification of bismuth with topological defects. Sci. Adv. 5, eaax6996 (2019).
Hsu, C.-H. et al. Topology on a new facet of bismuth. Proc. Natl. Acad. Sci. 116, 13255 (2019).
Collaudin, A., Fauqué, B., Fuseya, Y., Kang, W. & Behnia, K. Angle dependence of the orbital magnetoresistance in bismuth. Phys. Rev. X 5, 021022 (2015).
Spathelf, F., Fauqué, B. & Behnia, K. Magneto-seebeck effect in bismuth. Phys. Rev. B 105, 235116 (2022).
Kang, W., Spathelf, F., Fauqué, B., Fuseya, Y. & Behnia, K. Boundary conductance in macroscopic bismuth crystals. Nat. Commun. 13, 189 (2022).
Fuseya, Y., Katsuno, H., Behnia, K. & Kapitulnik, A. Nanoscale turing patterns in a bismuth monolayer. Nat. Phys. 17, 1031 (2021).
Savenko, V. S., Uglov, V. V., Ostrikov, O. M. & Khodoskin, A. P. Twinning of bismuth single crystals bombarded by boron ions. Tech. Phys. Lett. 24, 287 (1998).
Bashmakov, V. I. & Chikova, T. S. Evolution of twin-boundary geometry in bismuth crystals with an increase of load. Crystallogr. Rep. 47, 489 (2002).
Edelman, V. S., Sharvin, D. Y., Khlyustikov, I. N. & Troyanovskii, A. M. STM revealing of twin microlayers with quantized width on cleaved bismuth surface. Europhys. Lett. 34, 115 (1996).
Edelâman, V. S. Scanning tunnel microscopy and spectroscopy of an atomically clean bismuth surface. Phys. Usp. 48, 1057 (2005).
Halperin, B. I. Possible states for a three-dimensional electron gas in a strong magnetic field. Jpn. J. Appl. Phys. 26, 1913 (1987).
MacDonald, A. H. & Bryant, G. W. Strong-magnetic-field states of the pure electron plasma. Phys. Rev. Lett. 58, 515 (1987).
Behnia, K., Balicas, L. & Kopelevich, Y. Signatures of electron fractionalization in ultraquantum bismuth. Science 317, 1729 (2007).
Li, L. et al. Phase transitions of Dirac electrons in bismuth. Science 321, 547 (2008).
Alicea, J. & Balents, L. Bismuth in strong magnetic fields: Unconventional Zeeman coupling and correlation effects. Phys. Rev. B 79, 241101 (2009).
Sharlai, Y. V. & Mikitik, G. P. Origin of the peaks in the Nernst coefficient of bismuth in strong magnetic fields. Phys. Rev. B 79, 081102 (2009).
Zhu, Z., Fauqué, B., Fuseya, Y. & Behnia, K. Angle-resolved landau spectrum of electrons and holes in bismuth. Phys. Rev. B 84, 115137 (2011).
Zhu, Z. et al. Landau spectrum and twin boundaries of bismuth in the extreme quantum limit. Proc. Natl. Acad. Sci. 109, 14813 (2012).
Fauqué, B. et al. Hall plateaus at magic angles in bismuth beyond the quantum limit. Phys. Rev. B 79, 245124 (2009).
Kuchler, R. et al. Thermodynamic evidence for valley-dependent density of states in bulk bismuth. Nat. Mater. 13, 461 (2014).
Zhu, Z. et al. Emptying Dirac valleys in bismuth using high magnetic fields. Nat. Commun. 8, 15297 (2017).
Iwasa, A. et al. Thermodynamic evidence of magnetic-field-induced complete valley polarization in bismuth. Sci. Rep. 9, 1672 (2019).
Yaguchi, H. & Singleton, J. A high-magnetic-field-induced density-wave state in graphite. J. Phys.: Condensed Matter 21, 344207 (2009).
Fauqué, B. et al. Two phase transitions induced by a magnetic field in graphite. Phys. Rev. Lett. 110, 266601 (2013).
Wang, J. et al. Critical point for Bose-Einstein condensation of excitons in graphite. Proc. Natl. Acad. Sci. 117, 30215 (2020).
Zhu, Z. et al. Graphite in 90 t: Evidence for strong-coupling excitonic pairing. Phys. Rev. X 9, 011058 (2019).
Gourgout, A. et al. Magnetic freeze-out and anomalous hall effect in ZrTe5. npj Quantum Mater. 7, 71 (2022).
Huang, X. et al. Magnetic electron collimation in three-dimensional semi-metals. npj Quantum Mater. 5, 12 (2020).
Kinoshita, Y. et al. Field-induced reentrant insulator state of a gap-closed topological insulator (Bi1âxSbx) in quantum-limit states. Phys. Rev. B 107, 125140 (2023).
Banerjee, A. et al. Transport anomalies across the quantum limit in semimetallic Bi0.96Sb0.04. Phys. Rev. B 78, 161103 (2008).
Fauqué, B., Vignolle, B., Proust, C., Issi, J.-P. & Behnia, K. Electronic instability in bismuth far beyond the quantum limit. N. J. Phys. 11, 113012 (2009).
Bompadre, S. G., Biagini, C., Maslov, D. & Hebard, A. F. Unambiguous determination of the g factor for holes in bismuth at high B/T. Phys. Rev. B 64, 073103 (2001).
Behnia, K., Méasson, M.-A. & Kopelevich, Y. Oscillating Nernst-Ettingshausen effect in bismuth across the quantum limit. Phys. Rev. Lett. 98, 166602 (2007).
Akiba, K. et al. Quantitative evaluation of Dirac physics in PbTe. Phys. Rev. B 98, 115144 (2018).
Zhu, Z. et al. Nernst quantum oscillations in bulk semi-metals. J. Phys.: Condensed Matter 23, 094204 (2011).
Fuseya, Y. et al. Origin of the large anisotropic g factor of holes in bismuth. Phys. Rev. Lett. 115, 216401 (2015).
Peierls, R. More Surprises in Theoretical Physics (Princeton Univ Press, Princeton, 1991).
Ostrikov, O. M. Effect of carbon ion bombardment and surface oxidation on the twinning rate of single-crystal bismuth. Tech. Phys. 44, 597 (1999).
Ostrikov, O. M. Twinning in bismuth single crystals induced by thermal cycling. Tech. Phys. 46, 1199 (2001).
Galev, V., Kozlov, V., Kolomoets, N., Skipidarov, S. Y. & Tsvetkova, N. Exponential temperature dependence of the coefficient of the Nernst-Ettingshausen transverse effect in bismuth. Sov. Phys.-JETP Lett. 33, 106 (1981).
Zhu, Z., Fauqué, B., Behnia, K. & Fuseya, Y. Magnetoresistance and valley degree of freedom in bulk bismuth. J. Phys. Cond. Matt. 30, 313001 (2018).
Edelman, V. S. Tunneling spectra at terrace boundaries on the bismuth surface. J. Exp. Theor. Phys. 107, 251 (2008).
Kittel, C. Introduction to solid state physics (John Wiley & Sons, Inc., New York, 2004).
Landauer, R. Residual resistivity dipoles. Zeitschrift für Physik B Condensed Matter 21, 247 (1975).
Landauer, R. Spatial variation of currents and fields due to localized scatterers in metallic conduction. IBM J. Res. Dev. 32, 306 (1988).
MacDonald, A. H., Rice, T. M. & Brinkman, W. F. Hall voltage and current distributions in an ideal two-dimensional system. Phys. Rev. B 28, 3648 (1983).
Katakura, I. et al. Development of high-speed polarizing imaging system for operation in high pulsed magnetic field. Rev. Sci. Instr. 81, 043701 (2010).
Acknowledgements
This work was supported by The National Key Research and Development Program of China (Grants No.2022YFA1403500 and 2016YFA0401704), the National Science Foundation of China (Grant No. 12004123, No. 11574097 and No. 51861135104), and by directors funding grant number 20120772 at LANL. N.H. and RMcD acknowledge support from the US-DOE BES âScience of 100Tâ program. B. F. acknowledges support from Jeunes Equipes de lâInstitut de Physique du Collége de France (JEIP). K. B. was supported by the Agence Nationale de la Recherche (ANR-19-CE30-0014- 04). Y. F. was supported by JSPS KAKENHI grants 23H04862, 23H00268, and 18KK0132. Work performed at National High Magnetic Field Laboratory was supported by National Science Foundation Cooperative Agreements DMR-0654118 and DMR-2128556 and the State of Florida. M. T. acknowledges support from the JSPS KAKENHI grant No. 23H04862.
Author information
Authors and Affiliations
Contributions
Z.Z. and K.B. conceived this work; Y.Y., Z.Z., J.W., P.N., L.X., H.Z., N.H. and R.D.M. performed the magnetoresistance experiment. A.Y. and Y.F. preformed the theoretical calculation. Y. K and M. T preformed the magneto-optics measurements. A.V.S, A.A., M.S.N, B.F. and K.B. preformed the Nernst effect measurement up to 45 T. B.F., D.L., and C.P. preformed ultrasound measurement. Z.Z., B.F., Y.F. and K.B. wrote the manuscript with input from all authors.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Ye, Y., Yamada, A., Kinoshita, Y. et al. High-field immiscibility of electrons belonging to adjacent twinned bismuth crystals. npj Quantum Mater. 9, 12 (2024). https://doi.org/10.1038/s41535-024-00625-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41535-024-00625-7