Symmetry-enforced topological nodal planes at the Fermi surface of a chiral magnet

Despite recent efforts to advance spintronics devices and quantum information technology using materials with non-trivial topological properties, three key challenges are still unresolved1–9. First, the identification of topological band degeneracies that are generically rather than accidentally located at the Fermi level. Second, the ability to easily control such topological degeneracies. And third, the identification of generic topological degeneracies in large, multisheeted Fermi surfaces. By combining de Haas–van Alphen spectroscopy with density functional theory and band-topology calculations, here we show that the non-symmorphic symmetries10–17 in chiral, ferromagnetic manganese silicide (MnSi) generate nodal planes (NPs)11,12, which enforce topological protectorates (TPs) with substantial Berry curvatures at the intersection of the NPs with the Fermi surface (FS) regardless of the complexity of the FS. We predict that these TPs will be accompanied by sizeable Fermi arcs subject to the direction of the magnetization. Deriving the symmetry conditions underlying topological NPs, we show that the 1,651 magnetic space groups comprise 7 grey groups and 26 black-and-white groups with topological NPs, including the space group of ferromagnetic MnSi. Thus, the identification of symmetry-enforced TPs, which can be controlled with a magnetic field, on the FS of MnSi suggests the existence of similar properties—amenable for technological exploitation—in a large number of materials.

Despite recent efforts to advance spintronics devices and quantum information technology using materials with non-trivial topological properties, three key challenges are still unresolved 1-9 . First, the identification of topological band degeneracies that are generically rather than accidentally located at the Fermi level. Second, the ability to easily control such topological degeneracies. And third, the identification of generic topological degeneracies in large, multisheeted Fermi surfaces. By combining de Haas-van Alphen spectroscopy with density functional theory and band-topology calculations, here we show that the non-symmorphic symmetries [10][11][12][13][14][15][16][17] in chiral, ferromagnetic manganese silicide (MnSi) generate nodal planes (NPs) 11,12 , which enforce topological protectorates (TPs) with substantial Berry curvatures at the intersection of the NPs with the Fermi surface (FS) regardless of the complexity of the FS. We predict that these TPs will be accompanied by sizeable Fermi arcs subject to the direction of the magnetization. Deriving the symmetry conditions underlying topological NPs, we show that the 1,651 magnetic space groups comprise 7 grey groups and 26 black-and-white groups with topological NPs, including the space group of ferromagnetic MnSi. Thus, the identification of symmetry-enforced TPs, which can be controlled with a magnetic field, on the FS of MnSi suggests the existence of similar properties-amenable for technological exploitation-in a large number of materials.
Nearly a century ago Wigner, von Neumann and Herring 1,2 addressed the conditions under which Bloch states form degenerate band crossings, but their topological character and technological relevance has been recognized only recently 3-5 . To be useful 4-9 , tiny changes of a control parameter must generate a large response, underscoring the lack of control over the band filling as the unresolved key challenge in materials with band crossings known so far. This raises the question whether topological band crossings exist that are (1) generically located at the Fermi level, (2) separated sufficiently in the Brillouin zone (BZ) and (3) easy to control.
Natural candidates are systems with non-symmorphic symmetriesfor example, screw rotations-that generate positions in reciprocal space at which band-crossings are symmetry-enforced. The associated key characteristics include 10-17 : (1) the crossings are due to symmetry alone, that is, they occur on all bands independent of details such as chemical composition; (2) pairs of band crossings with opposite chirality are separated in k-space by about half a reciprocal lattice vector; (3) the band crossings may be enforced on entire planes 11,12 , forming so-called nodal planes (NPs) with non-zero topological charge; and (4) their existence may be controlled by means of symmetry breaking. Thus, if in a material the Fermi surfaces (FSs) cross such topological NPs, they enforce pairwise FS degeneracies with large Berry curvatures. The topology of these FS degeneracies, which we refer to as topological protectorates (TPs), will be independent of material-specific details and, moreover, may be controlled by symmetry breaking. The putative existence of topological NPs has been studied in phononic metamaterials [18][19][20] , and mentioned in a study of non-magnetic chiral systems focusing on Kramers-Weyl fermions 21 .
To demonstrate the formation of symmetry-enforced TPs at the intersection of NPs with the FS, we decided to study the ferromagnetic state of manganese silicide (MnSi), which has attracted great interest for its itinerant-electron magnetism 22 , helimagnetism, skyrmion lattice 23 and quantum phase transition 24 . Crystallizing in space group (SG) 198, MnSi is a magnetic sibling of non-magnetic RhSi (ref. 25 ), CoSi (ref. 26 ) and PdGa (ref. 27 ), in which sizeable Fermi arcs and multifold fermions were recently inferred from angle-resolved photoemission spectroscopy. MnSi is ideally suited for our study, as magnetic fields exceeding around 0.7 T stabilize ferromagnetism with magnetic screw-rotation symmetries enforcing NPs.
The topological charge ν of this duo of NPs (Fig. 1b) may be determined with the fermion doubling theorem 29 , which states that ν summed over all band crossings must be zero. We note that besides the NPs, there is an odd number of symmetry-enforced band crossings on the Y 1 -Γ-Y and R 1 -U-R lines forming Weyl points (ν = ±1) and four-fold points (ν = ±2), respectively (Fig. 1c, d, Extended Data Fig. 2, Supplementary Note 1). Moreover, due to the effective mirror symmetries, accidental Weyl points away from these high-symmetry lines must form pairs or quadruplets with the same ν. As the sum over ν of all of these Weyl and four-fold points is odd, the duo of NPs must carry a non-zero topological charge to satisfy the fermion doubling theorem. Hence, the duo of NPs at the BZ boundary is the topological partner of a single Weyl point on the Y 1 -Γ-Y line (Fig. 1b). This is a counter-example to Weyl semimetals, in which Weyl points occur always in pairs. Shown in Fig. 1d is the band structure of a generic tight-binding model satisfying SG 19.27 (Supplementary Note 2), where pairs of bands form NPs on the BZ boundaries k x = ±π and k z = ±π, whereas on the Y 1 -Γ-Y and R 1 -U-R lines there are Weyl and four-fold points, respectively. Explicit calculation of the Chern numbers shows that all of these band crossings, including those at the NPs, exhibit non-zero topological charges as predicted above. In turn, all of the FSs carry substantial Berry curvatures. The numerical analysis shows that these Berry curvatures become extremal at the NPs and close to the four-fold and Weyl points (Extended Data Fig. 3). By the bulk-boundary correspondence 3,4 , the non-trivial topology of these band crossings generates large Fermi arcs on the surface, which extend over half of the BZ of the surface (Extended Data Fig. 4). These arguments may be extended to 254 of the 1,651 magnetic SGs, of which 33 have NPs whose topological charges are enforced to be non-zero by symmetry alone (Supplementary Note 3).  5). Ten bands are found to cross the Fermi level (Fig. 1e). In agreement with our symmetry analysis and the tight-binding model (Fig. 1d), we find the same generic band crossings, namely: (1) NPs on the BZ boundaries k x = ±π and k z = ±π;

Calculated electronic structure
(2) an odd number of Weyl points along Y 1 -Γ-Y; and (3) an odd number of four-fold points along R 1 -U-R.
The calculated FSs as matched to experiment are shown in Fig. 1f, highlighting the NPs at the BZ boundaries at k x = ±π and k z = ±π (see Extended Data Table 1 for key parameters and Extended Data Fig. 5). Eight FS sheets centred at Γ comprise two small isolated hole pockets (sheets 1 and 2), two intersecting hole pockets with avoided crossings and magnetic breakdown due to spin-orbit coupling (sheets 3 and 4) and two pairs of jungle-gym-type sheets (sheets 5 and 6, and sheets 7 and 8). Sheets 9 and 10 are centred at R, comprising eight threefingered electron pockets around the [111] axes and a tiny electron pocket, respectively. The sheet pairs (5, 6), (7,8) and (9,10) Videos 1 and 2). For instance, rotating the moments within the x-y plane away from [010] breaks the magnetic screw rotation θC x 2 , but keeps θC z 2 intact. In turn, the TPs gap out on the k y = ±π and k x = ±π planes, whereas they remain degenerate at the k z = ±π planes (Extended Data Fig. 1

Experimental results
To experimentally prove the mechanism causing generic TPs at the intersection of the FS with symmetry-enforced NPs and their dependence on the direction of the magnetization, we mapped out the FS by means of the de Haas-van Alphen (dHvA) effect using capacitive cantilever magnetometry (Methods, Extended Data  30 . Figure 2b shows the oscillatory high-field part of the torque τ(1/B) at temperature T = 35 mK with the low-frequency components removed for clarity. To extract the dHvA frequencies, a fast Fourier transform (FFT) analysis of τ(1/B) was carried out, where the effects of demagnetizing fields and the unsaturated magnetization were taken into account (Methods). The FFT frequencies correspond to extremal FS cross-sections in low effective fields of about 0.7−1.9 T (Methods). Typical dHvA frequencies and FFT amplitudes, shown for φ = 82.5° in Fig. 2c, show five different regimes of dHvA frequencies labelled I to V. They comprise over 40 dHvA frequencies corresponding to different extremal FS orbits, as denoted by Greek letters (Fig. 2c, Extended Data Table 2). In our data analysis, we delineated artefacts due to the finite FFT window, such as the side lobes between κ 2 and 2κ 1 , or 3κ 2 and ξ 1 (Methods). Fitting the temperature dependence of the FFT amplitudes within Lifshitz-Kosevich theory 31 , the effective masses for each of the orbits were deduced ranging from m* = 0.4m e to m* = 17m e , where m e is the bare electron mass (Fig. 2d).
To relate the dHvA frequencies to the calculated FS orbits, the torque amplitude was inferred from the DFT band structure by means of the Lifshitz-Kosevich formalism, using small rigid band shifts of the order of 10 meV to improve the matching following convention (Methods, Extended Data Table 1). The assignment to experiment was based on the consistency between dHvA frequency, angular dispersion, strength of torque signal, field dependence of the dHvA frequencies, effective masses and presence of magnetic breakdown, as explained in Methods, Extended Data Table 2, Extended Data Figs. 6, 7, Supplementary Note 5. Figure 3a shows an intensity map of the experimental data of the (001) plane as a function of φ, where the theoretical dHvA branches are depicted by coloured lines (colours correspond to the FS sheets in Fig. 1). For comparison, Figure 3b shows an intensity map of the calculated dHvA spectra, where the experimental frequencies are marked by grey crosses.
For regimes I to IV, featuring contributions of the large FS sheets (5, 6) and (7,8), all frequencies may be assigned unambiguously (Extended Data Figs. 6, 7, Supplementary Note 5). Namely, regime I contains the loop orbits around U associated with pair (5, 6) (blue and orange) and the neck orbit of sheet 8 (yellow). Regime II exhibits the dHvA branches originating from neck orbits around Γ-Y-Γ on sheet 7 (purple). The neck orbits of sheet 8 (yellow), which evades detection because of the large slope of the dispersion, its high mass and the suppression of the magnetic torque near the [010] high-symmetry direction, is consistent with an anomalous frequency splitting at the expected crossing with the loop orbits of pair (5, 6) (blue and orange) around 6.5 kT (Supplementary Note 5). Regime III arises from both pairs (5, 6) and (7,8), that is, neck orbits around Γ-Y-Γ of (5, 6) and loop orbits around U of (7,8). The remaining cascade of frequencies in regime III reflects breakdown orbits (translucent yellow) arising from avoided crossings between sheets 3 and 4 (red and green). Regime IV is, finally, dominated by sheet 2 of the isolated hole pocket and the first harmonic of sheet 2.  )  35  294  500  800  1,600  2, As the magnetic torque generically vanishes at high-symmetry directions, which corresponds to the ⟨100⟩ axes in regimes I to IV, the associated FS sheets are centred at the Γ point. Likewise, the lowest frequency in regime V corresponds to a Γ-centred FS sheet, which can be assigned to the small hole pocket of sheet 1. In stark contrast, for regime V above about 0.05 kT, the high-symmetry directions correspond to the ⟨111⟩ axes, whereas the torque for the ⟨100⟩ axes is finite (see also Fig. 3b, Extended Data Fig. 5g). Hence, regime V is related to FS pockets in the vicinity of the R point that may be assigned to FS sheets (9,10). This allows for a basic estimate of the size and the effective mass of FS sheets (9, 10) without the need for a detailed account of their shape, completing the assignment. The calculations demonstrate the presence of symmetry-enforced crossings of sheets (9, 10) if they intersect the NPs (Fig. 1).
To confirm that we observed the entire FS, we calculated the Sommerfeld coefficient of the specific heat from the density of states at the Fermi level as rescaled by the measured mass enhancements (Extended Data Table 1). Excellent agreement is observed within a few percent of experiment 32 , γ ≈ 28 mJ mol −1 K −2 at B = 12 T. This analysis reveals, that sheets (5, 6), (7,8) and (9,10), which form TPs, contribute 86% to the total density of states at the Fermi level.

Topological NPs
Spectroscopic evidence of the symmetry-enforced topological band degeneracies at the BZ boundaries may be inferred from FS sheets (5, 6). Identical characteristics are observed for FS sheets (7,8) (Extended Data Fig. 7, Supplementary Note 5). We note that the dHvA cyclotron orbits are perpendicular to the NPs for fundamental reasons, piercing through them at specific points of the TPs. As shown in Fig. 4a, a magnetic field parallel to [010] leads to extremal cross-sections for FS sheets (5, 6), supporting cyclotron orbits in the vicinity of the U and the Y 1 points on planes depicted by blue and green shading, respectively. Centred with respect to the U point are possible cyclotron orbits comprising different segments of FS sheets 5 and 6, which interact at TP1 to TP4 with the BZ boundaries at k x = ±π and k z = ±π. In the absence of the non-symmorphic symmetries, these intersections would exhibit anticrossing and magnetic breakdown, leading to several orbits with different cross-sections and hence several dHvA frequencies. Instead, the behaviour is distinctly different to magnetic breakdown or Klein tunnelling 33,34 .
As the BZ boundaries at k x = ±π and k z = ±π represent symmetryenforced NPs, the crossing points of sheets 5 and 6 at TP1 to TP4 are, hence, protected band degeneracies at which the wavefunctions are orthogonal, that is, TP1 to TP4 are part of the TPs that suppress transitions between orbits (we call orbits containing at least one TP 'topological orbits'). In turn, two independent topological orbits (topological orbits 1 and 2) with identical areas and hence the same dHvA frequencies are expected (Fig. 4b, top). This is in excellent agreement with experiment, which shows a single dHvA frequency for field parallel [010] (φ = 90° in Fig. 4c). Rotating the direction of the magnetic field within the x-y plane away from [010], the NP at k x = ±π gaps out, whereas the NP at k z = ±π remains protected. Thus, the associated loop orbits around U (Fig. 4b, bottom) continue to include two points on the FS at k z = ±π (TP3 and TP4), leading to two additional topological orbits (topological orbits 3 and 4) of identical cross-section with the same dHvA frequency, in perfect agreement with the observed spectra (Fig. 4c).
Comparing the extremal cross-sections of the neck orbits around Γ-Y 1 -Γ with those around Γ-X-Γ, the latter crosses an NP whereas the former does not. With respect to Γ-X-Γ, there would be two extremal cross-sections with identical areas, positioned symmetrically with respect to X (Fig. 4d, top left), whereas for the cross-sections with respect to Γ-Y 1 -Γ there are two extremal orbits with different areas positioned asymmetrically with respect to Y 1 (Fig. 4d, bottom left). Thus, within our symmetry analysis and our DFT calculations, we expect  (Fig. 1f) from which the dHvA branch originates. The first harmonic (1sth) of the branches originating from sheet 2 is also labelled for clarity. A line cut of this colour map for fixed field angle φ = 82.5° is shown in Fig. 2c. More than 40 dHvA branches were observed as listed in Extended Data Table 2. b, Torque amplitudes of the dHvA spectra inferred from the ab initio band structure (Methods), as a function of f and φ, with the experimental frequencies of the dHvA branches indicated by crosses. To obtain a quantitative matching between theoretical and experimental dHvA branches, small rigid energy shifts to the ab initio bands were applied, as summarized in Extended Data Table 1, Supplementary Note 4. The detailed procedure how the experimental and theoretical dHvA branches were matched is described in the main text and in Supplementary Note 5.
Article a single dHvA branch for neck orbits parallel to a NP compared with two dHvA branches for neck orbits that are not parallel to a NP (Fig. 4d,  top right). Keeping in mind that only neck orbits around Y 1 are accessible experimentally, we clearly observe two branches, giving strong evidence that there are no NPs on the k y = ±π BZ boundary (Fig. 4e).

Concluding remarks
The  Fig. 4). These Fermi arcs reflect the presence of duos of NPs. Analogous Fermi arcs will not exist in non-magnetic materials with SG 198 [25][26][27] , which support trios of NPs (Supplementary Note 1).
In systems with symmetry-enforced NPs and TPs, tiny changes of the direction of the magnetization will control the topological band crossing in the bulk and the Fermi arcs, causing massive changes of Berry curvature that may be exploited technologically. The formation of TPs irrespective of the complexity of the FS raises the question of whether they affect the transport properties 38 and enable exotic states of matter 39 . Extending the analysis presented here to all 1,651 magnetic SGs, we find that there is a large number of candidate materials, such as CoNb 3 S 6 (ref. 40

Online content
Any methods, additional references, Nature Research reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at https://doi.org/10.1038/s41586-021-03543-x.

Fig. 4 | Extremal orbits and spectroscopic signatures of NPs and TPs.
Identical features presented here for sheet pair (5, 6) are also observed for FS sheet pair (7, 8) (Extended Data Fig. 7, Supplementary Note 5). a, FS sheet pair (5, 6) for a field (B) along the [010] direction (for an alternative colour shading see Extended Data Fig. 6d1). Planes illustrating loop-and neck-type orbits around the U point and the Γ-Y 1 -Γ line are indicated by blue and green shading, respectively. Loop orbits with respect to the U point intersect at TP1 to TP4 with the NPs on the k x = π and k z = π BZ boundaries. The NPs enforce degeneracies at TP1 to TP4, where the wave functions are orthogonal. b, Instead of anticrossing and magnetic breakdown, topological orbits stabilize. Top: cross-sectional areas under field along [010] at φ = 90°. Bottom: schematic cross-sectional areas under rotated field for φ = 80°. c, Intensity map of dHvA spectra in the regime of loop-and neck-type orbits around the U point (Fig. 3). The spectra are in excellent agreement with the topological orbits. No evidence for independent orbits of FS sheets 5 and 6 are observed. d, Top left: symmetrical positions of extremal orbits 1 and 2 in a plot of the FS cross-sectional area a along Γ-X-Γ with respect to the NP at the X point. The orbits give rise to identical dHvA frequencies. Note that these orbits are not accessible experimentally. Bottom left: asymmetrical position of extremal FS cross-sections along Γ-Y 1 -Γ with respect to the BZ boundary at k y = ±π. Top right: the associated orbits give rise to different dHvA frequencies. e, Intensity map of dHvA spectra in the regime of neck-type orbits around the Y 1 point (Fig. 3). Spectra are in excellent agreement with two orbits as shown in d (bottom left and top right), that is, no NP at the BZ boundary at k y = ±π containing Y 1 .

Sample preparation
For our study, two MnSi samples were prepared from a high-quality single crystalline ingot obtained by optical float-zoning 42 . The samples were oriented by X-ray Laue diffraction and cut into 1 × 1 × 1 mm 3 cubes with faces perpendicular to [

Experimental methods
Quantum oscillations of the magnetization, that is, the dHvA effect, was measured by means of cantilever magnetometry measuring the magnetic torque τ = m × B. The double-beam type cantilevers sketched in Extended Data Fig. 5e were obtained from CuBe foil by standard optical lithography and wet-chemical etching. The cantilever position was read out in terms of the capacitance between the cantilever and a fixed counter electrode using an Andeen-Hagerling AH2700A capacitance bridge, similar to the design described in refs. 43,44 .
Angular rotation studies were performed in a 3 He insert with a manual rotation stage at a base temperature T = 280 mK under magnetic fields up to 15 T. In addition, the effective charge carrier mass was determined using a dilution refrigerator insert with fixed sample stage under magnetic fields up to 14 T (16 T using a Lambda stage) at temperatures down to 35 mK.
We discuss partial rotations in the (001) and (110) crystallographic planes. The angle φ is measured from [100] in the (001) plane and the angle θ is measured from [001] in the (110) plane. Corresponding data are shown in Fig. 3a and Extended Data Fig. 5g. Owing to the topology of the FS and the simple cubic BZ, the (001) plane rotation shows most of the extremal orbits and is already sufficient for an assignment to the FS sheets. For this reason, the discussion of the dHvA data in the main text focuses on the rotation in the (001) plane.
The response of the cantilever was calibrated by means of the electrostatic displacement, taking into account the cantilever bending line obtained from an Euler-Bernoulli approach 45 . Applying a d.c. voltage, U, to the capacitance C 0 = ε 0 A/d 0 , defined by the area A, the plate distance d 0 and the vacuum permittivity ε 0 , leads to an electrostatic force F = C 0 U 2 /2d 0 . This force is equivalent to a torque τ = βFL, where L is the effective beam length and β = 0.78 is a geometry-dependent prefactor accounting for the different mechanical response of a bending beam to a torque and force, respectively. From this, the calibration constant K(C) = τ/ΔC quantifying the capacitance change ΔC in response to the torque was obtained for different values of C. Changes in K(C) up to 10% were recorded during magnetic field sweeps. The torque was calculated using (1)

Evaluation of the dHvA signal
The dependence of the capacitance, C(B ext ), was converted into torque and corrected as described below, where B ext is the applied magnetic field. An exemplary torque curve obtained at T = 280 mK and φ = 82.5° is shown in Fig. 2a. In the regime below B ≈ 0.7 T the transitions from helical to conical and field-polarized state generated a strongly hysteretic behaviour. At higher fields, magnetic quantum oscillations on different amplitude and frequency scales could be readily resolved. The first low-frequency components appeared at magnetic fields as low as B ≈ 4 T, whereas several high-frequency components, corresponding to larger extremal cross-sections, could only be resolved in high fields (Fig. 2b). Consequently, the data acquisition and evaluation was optimized by treating low-and high-frequency components separately.
To eliminate the non-oscillatory component of the signal, low-order polynomial fits or curves obtained by adjacent averaging over suitable field intervals were subtracted from the data, producing consistent results. FFTs of τ(1/B) were used to determine the frequency components contained in the signal. Field sweeps were performed from 0 T to 15 T at 0.03-0.04 T min −1 and from 15 T to 10 T at 0.008 T min −1 . FFTs over the range 4 T to 15 T (10 T to 15 T) were performed to evaluate frequency components below (above) f = 350 T for measurements in the 3 He insert and from 10 T to 14 T (11 T to 16 T with Lambda stage) in the dilution refrigerator. The values correspond to the applied field before taking into account demagnetization. Rectangular FFT windows were chosen to maximize the ability to resolve closely spaced frequency peaks. See Supplementary Note 4 for details.

Internal magnetic field and dHvA frequency f(B) in a weak itinerant magnet
MnSi is a weak ferromagnet with an unsaturated magnetization up to the largest magnetic fields studied. This results in two different peculiarities concerning the observed dHvA frequencies.
(1) The field governing the quantum oscillations is the internal field 31 , where M exp is the low-field value of the magnetization in the fieldpolarized phase determined experimentally. The applied field was corrected by this value. The field dependence of the magnetic moment yields only a minor correction of the internal field that may be neglected.
(2) The effect of the unsaturated magnetization on the Fermi surface is more prominent and may be described in a good approximation as a rigid Stoner exchange splitting that scales with the magnitude of the magnetization. Consequently, FS cross-sectional areas are enlarged with increasing B for the majority electron orbits and minority hole orbits. Cross-sectional areas shift downwards for majority hole and minority electron orbits. This change in cross-sectional area is not directly proportional to the change in the observed dHvA frequencies f, that is, the dHvA frequencies deviate from the field-dependent frequency where I is the Stoner exchange parameter, m b is the band mass, the ± is for electron and hole orbits, respectively, s = ±1 is the spin index and f 0 is the hypothetical frequency without exchange splitting. Note, that this model is only meaningful in the field-polarized regime B ≳ 0.7 T. Using the experimental M(B) curve of MnSi 32 , we estimate that the frequencies f(B) in the windows used for f > 350 T defined above with centre fields B average = 2B high B low /(B low + B high ) ranging from 11.8 T to 13.2 T correspond to the extremal cross-sections at B ≈ 1.7−1.9 T (Extended Data Fig. 5f).
For the window used for frequencies f < 350 T, it is B average = 6.5 T and f(B) corresponds to the extremal cross-sections at B ≈ 0.7 T. Thus, even under large magnetic fields, the experimental frequency values correspond to a field-polarized state in a low field.

Quantum oscillatory torque and Lifshitz-Kosevich equation
Evaluation and interpretation of the quantum oscillatory torque magnetization was performed using the Lifshitz-Kosevich formalism 31 . The components of M parallel (∥) and perpendicular (⊥) to the field are given by: where V is the sample volume, p is the harmonic index, A″ is the curvature of the cross-sectional area parallel to B, and f is the dHvA frequency observed (see comments above). The phase γ = 1/2 corresponds to a parabolic band. In general, the phase includes also contributions due to Berry phases when the orbit encloses topologically non-trivial structures in k-space. The ± holds for maximal and minimal cross-sections, respectively. The torque amplitude is given by τ osc = M osc,⊥ B. The torque thus vanishes in high-symmetry directions where f(φ) is stationary. This feature of τ may be used to infer additional information about the symmetry properties of a dHvA branch. R T describes the temperature dependence of the oscillations * describes the influence of a finite scattering time τ. Here, ω c =eB/m* is the cyclotron frequency.

DFT calculations
The band structure and FS sheets of MnSi in the field-polarized phase were calculated using DFT. The calculations included the effect of spin-orbit coupling. In all calculations, the magnetic part of the exchange-correlation terms was scaled 49 to match the experimental magnetic moment of 0.41μ B per Mn atom at low fields. As input for the DFT calculations, the experimental crystal structure of MnSi was used, that is, space group P2 1 3 (198) with an experimental lattice constant a = 4.558 Å. Both Mn and Si occupy Wyckoff positions 4a with coordinates (u, u, u), (−u + 1/2, −u, u + 1/2), (−u, u + 1/2, −u + 1/2), (u + 1/2, −u + 1/2, −u) where u Mn = 0.137 and u Si = 0.845 (Extended Data Fig. 5a). Calculations were carried out using WIEN2k 50 , ELK 51 and VASP 52,53 using different versions of the local spin density approximation. The results are consistent within the expected reproducibility of current DFT codes 54 . The remaining uncertainties motivate a comprehensive experimental FS determination as reported in this study. In the main text, we focus on the results obtained with WIEN2k, using the local spin density approximation parametrization of Perdew and Wang 55 and a sampling of the full BZ with a 23 × 23 × 23 Γ-centred grid. The results of Extended Data Figs. 1, 2, 4 were obtained using VASP with the PBE functional 56 and a BZ sampling with a 15 × 15 × 15 k-mesh centred around Γ.
Bands used for the determination of the Fermi surface were calculated with WIEN2k on a 50 × 50 × 50 k-mesh. Owing to the presence of spin-orbit coupling, but the absence of both inversion and time-reversal symmetry, band structure data had to be calculated for different directions of the spin quantization axis. For a given experimental plane of rotation, calculations were performed in angular steps of 10°. The bands were then interpolated k-point-wise using third-order splines to obtain band structure information in 1° steps.
For the prediction of the dHvA branches from the DFT results, the Supercell k-space Extremal Area Finder (SKEAF) 57 was used on interpolated data corresponding to 150 × 150 × 150 k-points in the full BZ. The theoretical torque amplitudes shown in Fig. 3b were calculated directly from the prefactors in equations (4) and (5) convoluted with a suitable distribution function.
To compute the surface states of MnSi in the field-polarized phase (Extended Data Fig. 4), we first constructed a DFT-derived tight-binding model using the maximally localized Wannier function method as implemented in Wannier90 58 . Using this tight-binding model, we computed the momentum-resolved surface density of states by means of an iterative Green's function method, using WannierTools 59 . The symmetry eigenvalues of the DFT bands were computed from expectation values using VASP pseudo wavefunctions, as described in ref. 60 .

Magnetic breakdown
The probabilities for magnetic breakdown at a junction i is given by The probability for no breakdown to occur is thus q i = 1 − p i . The breakdown fields B 0 were calculated from Chamber's formula where k g is the gap in k-space and a and b are the curvatures of the trajectories at the breakdown junction 31 . In our study of MnSi, we observed magnetic breakdown in particular between sheets 3 and 4, which exhibit up to eight junctions depending on the magnetic field direction and between FS sheet pairs touching the BZ surfaces on which the NP degeneracy is lifted. Only breakdown orbits that are closed after one cycle are considered in the analysis. Further details can be found in the Supplementary Note 5.

Assignment of dHvA orbits and rigid band shifts
The assignment of the experimental dHvA branches to the corresponding extremal FS cross-sections was based on the following criteria: (1) dHvA frequency-determining sheet size in terms of the cross-sectional area; (2) angular dispersion-relating to sheet shape, topology and symmetry; (3) torque signal strength-relating to sheet shape and symmetry; (4) direction of f(B) shift-relating to spin orientation and charge carrier type; (5) effective mass-relating to the temperature dependence; (6) magnetic breakdown behaviour-relating to proximity of neighbouring sheets. The majority of the observed dHvA branches could be related directly to the FS as calculated. In addition, we used the well-established procedure of small rigid band shifts to optimize the matching. While this procedure is, in general, neither charge nor spin conserving, it results in a very clear picture of the experimental FS. One has to bear in mind, however, that the deviations between the true FS and the calculated FS are not due to a rigid band shift (this might be justified, for example, in case of unintentional doping, which we rule out here). Rather, it may be attributed to differences in the band dispersions that originate in limitations of our DFT calculations (for example, neglecting electronic correlations and the coupling to the spin fluctuation spectrum).