Introduction

In most of low-dimensional, direct-gap semiconductor systems (quantum wells and monolayers of transition metal dichalcogenides, quantum wires, and dots), the near band edge electron–hole excitations, the excitonic states, can be classified into two categories depending whether their associated projection of the total angular momentum is jz = ±1 or it is jz = ±2 or 0. According to standard selection rules for optical transitions, the excitons with jz = ±1 can be excited by light and/or can recombine by emitting photons and are therefore referred to as bright excitons. Instead, excitons with jz = ±2 or 0 are termed dark as they do not efficiently couple to the radiation field. Although dark, the jz = ±2 or jz = 0 excitons can be brightened under special conditions, e.g., under application of the magnetic field in an appropriate configuration1,2 and thus possibly initiating the photoluminescence signal (though still hardly observed in absorption-type spectra). Dark excitons can largely alter the dynamics of optical excitations and, as long-lived bosonic quasiparticles, are possible candidates to form a condensate exciton phase3. Those states have been the object of vast investigations in semiconductor quantum wells, dots, and wires, in carbon nanotubes and more recently in monolayer semiconductors4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19. All these studies have been however largely focused on excitonic ground states, whereas the physics of dark excitons associated with the excited excitonic states is a little-explored area20,21,22. This, in particular, applies to two-dimensional semiconductors which often display a characteristic series of Rydberg ns states of bright excitons23,24,25,26, whereas a demonstration of the corresponding series of dark excitonic states is missing so far.

In this paper, we report on low-temperature magneto-optical study of Rydberg series of excitonic states as they appear in a WSe2 monolayer. Notably, in this representative example of a 2D semiconductor, the rich spectrum of excited excitonic states can be relatively easily traced with luminescence experiments25,27,28,29. Profiting on this fact, we brighten dark excitons in our sample by applying an in-plane magnetic field and uncover their ns (n = 1, up to 4) Rydberg series which appears in addition to a Rydberg spectrum of bright excitons24,25. By comparing energy diagrams of dark and bright exciton series, we experimentally derive the single-particle separation between spin–orbit–split conduction band subbands in a WSe2 monolayer. Intriguingly, this separation is found to be of about two times smaller than that commonly assumed from theoretical modeling of electronic bands of monolayer semiconductors.

Results and discussion

The focus of our study can be clearly read from Fig. 1. Figure 1a illustrates a simple generic scheme of the band structure of a WSe2 monolayer, as it is found in a close vicinity of the direct bandgap at the two K + and K points of the Brillouin zone. We consider only three relevant bands (per each K + and K point): two spin–orbit split subbands in the conduction band (weakly separated in energy, by \({\Delta}_{{{{{{{\mathrm{c}}}}}}}} \approx\) tens of meVs) but only one, top valence subband (neglecting the lower one separated by hundreds of meVs from the upper subband). As generally accepted and shown in Fig. 1a, the interband, intravalley transitions which involve the electronic states of the upper conduction subband are optically allowed, whereas those associated with the lower subband are characterized by jz = 0 and thus are not optically active under conventional conditions/configurations7,30. The jz = 0 transitions are weakly allowed only in the configuration for which the electromagnetic wave has its k-vector aligned along the monolayer plane and nonzero out of plane component of the electric field, or in the presence of an in-plane magnetic field. Please, note, that in the case of TMD monolayers, the jz = ±1 selection rules for optical transitions coincide with the ansatz of conservation of electronic spins in optical transitions. The band-to-band image of optical transitions is not, however, sufficient to describe the optical response of semiconductors, and in particular that of the investigated structure: the electron–hole Coulomb interaction and hence the excitonic effects have to be taken into account. As shown in Fig. 1b, one may expect the appearance of two distinct ns Rydberg series of near band edge excitons in our WSe2 monolayer: the routinely observed bright exciton series23,24,25,26 associated with the upper conduction band subband, but also the dark-exciton series associated with the lower conduction band subband, revealed in the present report. In the limit of large n, the energy difference between bright and dark-exciton states provides a measure of the spin–orbit splitting of the conduction band \({\Delta}_{{{{{{{\mathrm{c}}}}}}}} = E_{{{{{{{\mathrm{g}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}} - E_{{{{{{{\mathrm{g}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}\), where \(E_{{{{{{{\mathrm{g}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}}\) and \(E_{{{{{{{\mathrm{g}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}\) denote, correspondingly, the single-particle bandgaps associated with upper and lower conduction band subbands. It is worth noting, that the energy difference \({\Delta}_{{{{{{{{\mathrm{DB}}}}}}}}}\) between the ground states of dark and bright excitons depends also on their respective, \(E_{{{{{{{\mathrm{b}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}\) and \(E_b^B\), binding energies: \({\Delta}_{{{{{{{{\mathrm{DB}}}}}}}}} = {\Delta}_{{{{{{{\mathrm{c}}}}}}}} + \left( {E_{{{{{{{\mathrm{b}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}} - E_b^{{{{{{{\mathrm{B}}}}}}}}} \right)\). Since \(E_{{{{{{{\mathrm{b}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}\) and \(E_b^{{{{{{{\mathrm{B}}}}}}}}\) are a priori different, measuring the \({\Delta}_{{{{{{{{\mathrm{DB}}}}}}}}}\) parameter alone does not infer the amplitude of \({\Delta}_{{{{{{{\mathrm{c}}}}}}}}\).

Fig. 1: Schematic illustrations of bright and dark transitions in WSe2 monolayer.
figure 1

a Schematic representation of the interband transitions in a single-particle picture. The spin-up and spin-down states are shown as light-orange and light-blue parabolas, respectively. Optically active (bright) transitions at K+ and K valleys are represented by the red arrows, while optically inactive (dark) transitions are depicted in dark gray. The light-gray arrow represents the conduction band spin–orbit splitting (\({\Delta}_{{{{{{{\mathrm{c}}}}}}}}\)). b The energy diagram of excitonic transitions: Rydberg series of bright (\(n{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}}\)) and dark (\(n{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}\)) resonances. The horizontal red and dark-gray lines represent bright and dark excitonic levels, respectively. The light-gray arrows represent the bandgaps for bright and dark excitons (\(E_{{{{{{{\mathrm{g}}}}}}}}^{{{{{{{{\mathrm{B/D}}}}}}}}}\)) and the binding energies of the ground states of bright and dark excitons (\(E_{{{{{{{\mathrm{b}}}}}}}}^{{{{{{{{\mathrm{B/D}}}}}}}}}\)). The energy splitting between ground states of bright and dark excitons (\({\Delta}_{{{{{{{{\mathrm{DB}}}}}}}}}\)) and the conduction band spin–orbit splitting (\({\Delta}_{{{{{{{\mathrm{c}}}}}}}}\)) are represented by the green arrows.

The active part of the structure used for the experiments is a WSe2 monolayer extracted from a commercially available tungsten diselenide crystal. The monolayer was encapsulated in between hexagonal boron nitride (hBN) flakes and deposited on a Si substrate, following the conventional methods of mechanical exfoliation and deterministic dry transfer techniques. Experiments consisted of low-temperature (4.2 K) magneto-photoluminescence measurements carried out in magnetic fields up to 30 T, applied along the monolayer plane. As sketched in the inset to Fig. 2, we used the backscattering geometry where both the excitation and the collected light beams were quasi-perpendicular to the sample/monolayer. More details on sample preparation and experimental details can be found in the “Methods” section.

Fig. 2: Magnetic field-dependence of emissions from bright and dark-exciton Rydberg states.
figure 2

Low-temperature photoluminescence (PL) spectra in arbitrary (arb.) units measured at selected in-plane magnetic fields in the energy region of both a ground and b excited states. Spectra were shifted for clarity and normalized to the 1sB feature in the lower energy region and to the 2sB feature in the higher-energy region. Note that the intensity scales differ for these two regions. The inset schematically shows the backscattering geometry used in the experiment, where the excitation and the collected light beams are quasi-perpendicular to the sample surface.

Representative photoluminescence spectra of the investigated WSe2 monolayer are presented in Fig. 2. Inspecting the low spectral range (Fig. 2a), one recognizes a characteristic set of multiple photoluminescence transitions, typical of WSe2 monolayers31,32,33. Possible formations of many different excitonic states (including nontrivial indirect/inter-valley excitons and excitonic complexes such as trions and biexcitons) account for the spectrum complexity which we disregard here, focusing our attention on neutral, bright, and dark intravalley excitons. First, we examine the photoluminescence (PL) spectrum measured in the absence of magnetic field (B), and, in accordance to previous reports24,25,27,28,29, recognize the characteristic sequence of emission peaks related to the Rydberg series of \(1{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}},2{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}}\), \(3{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}}\), up to \(4{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}}\) states of bright excitons. The PL spectrum measured at B = 0 shows also a weak transition related to the ground state, \(1{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}\), of dark exciton6,7,8. Note that no additional emission features can be found between \(1{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}}\) and \(2{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}}\) states (see Supplementary Fig. S3). This is because this low energy state is largely populated (under optical excitation) and because the emission from this state is legitimate in the direction along the monolayer plane, thus weakly visible in our spectrum collected using the objective with a relatively large numerical aperture7,8. Please note, that we employ here a simplified picture of the dark exciton and view it as raising only a single resonance, i.e., we ignore its doublet, dark-gray components split by ≈0.7 meV6,8, the energy which is negligibly small in the context of our further considerations.

Central for our report is the study of the evolution of the PL spectra measured as a function of the in-plane magnetic field (\(B_\parallel\)). Such studies have been already successfully applied to unveil the ground (\(1{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}\)) states of dark excitons in most of TMD monolayers8,9,10,11,34. The application of the in-plane magnetic field induces the Zeeman effect, and thus the mixing of spin–orbit split conduction subbands; dark excitons acquire an oscillator strength (\(\sim B_\parallel ^2\)) for the emission in the direction perpendicular to the monolayer plane and become apparent in the PL spectra. The efficient magnetic brightening of the ground state of the dark exciton in WSe2 monolayer is confirmed by the present experiments. As can be seen in Fig. 2a, the intensity of the \(1{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}\) emission grows progressively as a function of \(B_\parallel\), and eventually dominates the PL spectrum in the limit of high magnetic fields. Importantly, new transitions driven by the application of the in-plane magnetic field appear also in the higher spectral range (see Fig. 2b) in the vicinity of the PL peaks associated to the excited states of bright excitons. These new transitions, labeled as \(2{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}\), \(3{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}\), and \(4{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}\) in Fig. 2, are attributed to the excited Rydberg states of the dark exciton.

As demonstrated in Fig. 3, the magnetic brightening of dark-exciton states, i.e., the increase of the relative intensity \(I_{{{{{{{\mathrm{D}}}}}}}}(n)/I_{{{{{{{\mathrm{B}}}}}}}}\left( n \right)\) of dark \(n{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}\) emission peak with respect to its bright \(n{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}}\) counterpart, roughly follows the expected \(I_{{{{{{{\mathrm{D}}}}}}}}(n)/I_{{{{{{{\mathrm{B}}}}}}}}\left( n \right) = b_n\;B_\parallel ^2\) rule (see Supplementary Note 1 and Supplementary Figs. S1 and S2 for more details on data analysis and estimation of the PL intensities). This, however, is not true in the case of ground states for the in-plane fields \(B_\parallel > 20\,{{{{{{{\mathrm{T}}}}}}}}\), where the relative intensity deviates from the quadratic dependence, possibly due to the shortened (by optical activation) dark-exciton lifetime leading to the decrease of the dark-exciton population. Hence, only the datapoints for \(B_\parallel < 20\,{{{{{{{\mathrm{T}}}}}}}}\) were used for the estimation of b1. The bn coefficients reflect the relative oscillator strength \(r_n = r_{n{{{{{{{\mathrm{s}}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}\)/\(r_{n{{{{{{{\mathrm{s}}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}}\) as well as the relative population \(p_n = N_{n{{{{{{{\mathrm{s}}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}/N_{n{{{{{{{\mathrm{s}}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}}\) (under optical excitation) of corresponding dark \(n{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}\) and bright \(n{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}}\) states: \(b_n =\) \(r_n \cdot p_n\). Comparing the extracted \(b_n\) parameters with theoretically expected values for \(r_n\) (see Supplementary Note 2), we find that dark and bright states are not far from being equally populated for pairs of the exited states (\(p_n \approx 2\) for \(n = 2,3\)). On the other hand, the population of the ground dark state (\(1{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}\)) is estimated to be five orders of magnitude higher than that of the \(1{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}}\) ground bright state (\(p_n \approx 10^5\)). As discussed in more detail in Supplementary Note 2, the population of excitonic states under optical transitions in WSe2 monolayers appears to be highly nonthermal, which needs its explanation with a future comprehensive study of PL dynamics in these systems.

Fig. 3: Brightening of the dark-exciton Rydberg states.
figure 3

Intensities of dark states relative to the intensity of bright states (Rel. int. = \(I_{{{{{{{\mathrm{D}}}}}}}}(n)/I_{{{{{{{\mathrm{B}}}}}}}}\left( n \right)\)) as a function of the in-plane magnetic field for a the 1s states (blue squares), b the 2s states (red circles), and c the 3s states (green triangles) along with quadratic fits (black lines) of the form \(a + bx^2\). Note that in the case of 1s states, only datapoints from 0 T to 20 T were considered for the fitting, as for higher fields, the relative intensity clearly starts to deviate from quadratic dependence. The error bars represent combined standard uncertainties resulting from the uncertainties of the estimation of the photoluminescence intensities of bright and dark-exciton states. More details on this estimation can be found in Supplementary Note 1.

Energy diagrams of both dark and bright exciton series are of particular interest. Central positions of of \(1{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}\) and \(1{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}}\) emission peaks (at \(E_{1{{{{{{{\mathrm{s}}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}} = 1.682\;{{{{{{{\mathrm{eV}}}}}}}}\) and \(E_{1{{{{{{{\mathrm{s}}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}} = 1.725\;{{{{{{{\mathrm{eV}}}}}}}}\), respectively) can be directly read from the measured spectra with reasonable precision, whereas the central energies of transitions associated with excited excitonic states have been extracted by reproducing the high energy part of the spectra with multiple Lorentzian resonances (see Supplementary Fig. S2 and Supplementary Table S1). As shown in Fig. 4, both series of bright and dark excitonic resonances follow the previously evoked (for bright states) semi-empirical formula25:

$$E_n^{{{{{{{{\mathrm{D}}}}}}}}/{{{{{{{\mathrm{B}}}}}}}}} = E_{{{{{{{\mathrm{g}}}}}}}}^{{{{{{{{\mathrm{D}}}}}}}}/{{{{{{{\mathrm{B}}}}}}}}} - \frac{{Ry^{{{{{{{{\mathrm{D}}}}}}}}/{{{{{{{\mathrm{B}}}}}}}}}}}{{\left( {n + \delta } \right)^2}},$$
(1)

where \(E_{{{{{{{\mathrm{g}}}}}}}}^{{{{{{{{\mathrm{D}}}}}}}}/{{{{{{{\mathrm{B}}}}}}}}}\) is the dark/bright bandgap associated to the lower/upper conduction band subband and, \(Ry^{{{{{{{{\mathrm{D}}}}}}}}/{{{{{{{\mathrm{B}}}}}}}}}\) should be, in the first approximation, identified with the effective Rydberg energy for the dark/bright exciton; \(Ry^{{{{{{{{\mathrm{D}}}}}}}}/{{{{{{{\mathrm{B}}}}}}}}} = 13.6\;{{{{{{{\mathrm{eV}}}}}}}}\cdot \mu ^{{{{{{{{\mathrm{D}}}}}}}}/{{{{{{{\mathrm{B}}}}}}}}}/\left( {\varepsilon ^2m_0} \right)\), where \(\mu ^{{{{{{{{\mathrm{D}}}}}}}}/{{{{{{{\mathrm{B}}}}}}}}}\) is the reduced effective mass embracing the valence band hole and electron from upper/lower conduction subbands, ε is the dielectric constant of the surrounding hBN material and \(m_0\) stands for the electron mass. δ parameter, which accounts for the 2D screening length and the hBN-dielectric screening, is found to be \(\delta = - 0.09\), which is notably common for dark and bright series and well matches the previously reported value (\(- 0.08\) for the bright exciton series25). The validity of the above simple formula to account for the energy ladder of 2D excitonic states in TMD monolayers encapsulated in hBN has been largely justified in ref. 25, and, in particular, shown to be equivalent to the optional approach employing the Rytova–Keldysh potential24,27. As extrapolated from plots shown in Fig. 4, the extracted values of dark and bright bandgaps are, correspondingly, \(E_{{{{{{{\mathrm{g}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}} = 1.8802\; \pm 0.0004\;{{{{{{{\mathrm{eV}}}}}}}}\) and \(E_{{{{{{{\mathrm{g}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}} = 1.8940\; \pm 0.0006\;{{{{{{{\mathrm{eV}}}}}}}}\) and we eventually come up with disclosing the amplitude of the conduction band spin–orbit splitting: \({\Delta}_{{{{{{{\mathrm{c}}}}}}}} = E_{{{{{{{\mathrm{g}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}} - E_{{{{{{{\mathrm{g}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}} = 14 \pm 1\;{{{{{{{\mathrm{meV}}}}}}}}\). Notably, this value is to be accounted by theoretical modeling of electronic bands of TMD monolayers35,36,37 which, as yet, overestimates it by a factor of two and a half. It is, however, important to note that the theoretical estimate35,36,37 of \({\Delta}_{{{{{{{\mathrm{c}}}}}}}} \cong 36\;{{{{{{{\mathrm{meV}}}}}}}}\), refers to freestanding WSe2 monolayers, whereas, \({\Delta}_{{{{{{{\mathrm{c}}}}}}}} = 14\;{{{{{{{\mathrm{meV}}}}}}}}\) we report, is, strictly speaking, derived for the case of WSe2 monolayer embedded in hBN. The dielectric screening which is well evidenced to affect the excitonic resonances in encapsulated TMD monolayers, leads also to the renormalization of, single-particle bandgaps38,39. It remains to be clarified with advanced theoretical calculations whether the dielectric screening may differently renormalize our dark and bright bandgaps and whether the value \({\Delta}_{{{{{{{\mathrm{c}}}}}}}} = 14\;{{{{{{{\mathrm{meV}}}}}}}}\) we report here applies also to other structures of WSe2 monolayers.

Fig. 4: Dark and bright excitonic Rydberg series.
figure 4

Experimentally obtained transition energies for the bright (red diamonds, \(n{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}}\)) and dark (gray squares, \(n{{{{{{{\mathrm{s}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}\)) exciton s-series as a function of \(1/\left( {n + \delta } \right)^2\), where \(\delta = - 0.09\). The error bars represent the line widths of the observed emission lines. The solid red and gray lines represent fits to the data with the model described in the text. The dashed red and gray lines denote the bandgaps for dark (\(E_{{{{{{{\mathrm{g}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}} = 1.880\;{{{{{{{\mathrm{eV}}}}}}}}\)) and bright (\(E_{{{{{{{\mathrm{g}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}} = 1.894\;{{{{{{{\mathrm{eV}}}}}}}}\)) excitons, while the red and gray arrows show the binding energies of bright (\(E_{{{{{{{\mathrm{b}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}} = 171\;{{{{{{{\mathrm{meV}}}}}}}}\)) and dark \(\left( {E_{{{{{{{\mathrm{b}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}} = 198\;{{{{{{{\mathrm{meV}}}}}}}}} \right)\) excitons, respectively. The conduction band splitting (\({\Delta}_{{{{{{{\mathrm{c}}}}}}}} = 14\;{{{{{{{\mathrm{meV}}}}}}}}\)) is also indicated.

Binding energies of bright and dark excitons: \(E_{{{{{{{\mathrm{b}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}} = 170.6 \pm 0.4\;{{{{{{{\mathrm{meV}}}}}}}}\) and \(E_{{{{{{{\mathrm{b}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}} = 198.4 \pm 0.8\;{{{{{{{\mathrm{meV}}}}}}}},\) respectively, are two other extracted parameters. Following up our analysis (see Eq. (1) and Fig. 4), the resulting energy ratio \(\rho _{{{{{{{{\mathrm{DB}}}}}}}}} = E_{{{{{{{\mathrm{b}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}/E_{{{{{{{\mathrm{b}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}} = 1.16\) (or difference \(E_{{{{{{{\mathrm{b}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}} - E_{{{{{{{\mathrm{b}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}} = 28\;{{{{{{{\mathrm{meV}}}}}}}}\)) is anticipated to be explained by different, for bright and dark excitons, reduced effective masses \(\left( {E_{{{{{{{\mathrm{b}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}/E_{{{{{{{\mathrm{b}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}} = \mu ^{{{{{{{\mathrm{D}}}}}}}}/\mu ^{{{{{{{\mathrm{B}}}}}}}}} \right)\). Indeed, when assuming the theoretically calculated masses of a valence hole (\(m_h = 0.36\)) and electrons in two conduction subbands (\(m_{{{{{{{\mathrm{e}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}} = 0.40\) and \(m_{{{{{{{\mathrm{e}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}} = 0.29\))37, we find \(\mu ^{{{{{{{\mathrm{D}}}}}}}}/\mu ^{{{{{{{\mathrm{B}}}}}}}} = 1.18\), a value which is very close to the experimentally derived ratio \(E_{{{{{{{\mathrm{b}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}/E_{{{{{{{\mathrm{b}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}} = 1.16\). On the other hand, one may expect that the difference between \(E_{{{{{{{\mathrm{b}}}}}}}}^{{{{{{{\mathrm{D}}}}}}}}\) and \(E_{{{{{{{\mathrm{b}}}}}}}}^{{{{{{{\mathrm{B}}}}}}}}\) arises not only from the direct Coulomb term but also from the exchange term, as bright and dark excitons consist of parallel or antiparallel spin configurations, respectively40,41,42. Whereas we suggest here that the exchange term might be relatively small, this claim has be to taken with caution, since our estimations refer to theoretically calculated effective masses which perhaps are not sufficiently accurate. Notably, the amplitude of the exchange term, as large as \(20\;{{{{{{{\mathrm{meV}}}}}}}}\), has been theoretically predicted40 but for the case of freestanding WSe2 monolayers. Coulomb interaction, being effectively at the origin of the exchange term, is screened in encapsulated monolayers. With very rough estimation one would expect the reduction of the exchange term down to \(20\;{{{{{{{\mathrm{meV/}}}}}}}}4.5 = 4.4\;{{{{{{{\mathrm{meV}}}}}}}}\) in our monolayer encapsulated in hBN with dielectric constant \(\varepsilon = 4.5\). Nonetheless, our findings once again raise a question about the importance of the exchange interaction for the energy difference between the bright and dark excitons binding energies in TMD monolayers.

Conclusions

In this report, we demonstrated the magnetic brightening of the Rydberg ns-series of dark excitons, up to n = 4, in a WSe2 monolayer encapsulated in hBN. The analysis of the bright and dark excitons series allowed us to determine one of the missing band parameters, the amplitude of the spin–orbit splitting in the conduction band of this broadly investigated 2D structure. Its derived value, \({\Delta}_{{{{{{{\mathrm{c}}}}}}}} = 14\;{{{{{{{\mathrm{meV}}}}}}}}\), is significantly lower than that commonly assumed, what calls for revision of theoretical calculations of electronic bands in TMD monolayers. Moreover, our results suggest that the difference between the binding energies of bright and dark excitons can be fully explained by the difference in the masses of electrons in the two spin–orbit-split conduction bands, without referring to exchange interactions.

Methods

Sample fabrication

The sample used in our study is composed of WSe2 monolayer encapsulated in hBN flakes and deposited on a Si substrate. The heterostructure was fabricated using mechanical exfoliation of bulk WSe2 and hBN crystals. The bottom hBN layer was created by nondeterministic exfoliation onto Si substrate, while the WSe2 monolayer and capping hBN flake were transferred using dry stamping technique onto the bottom hBN flake using PDMS stamps.

Experimental setup

The low-temperature magneto-photoluminescence experiments were performed in the Voigt configuration, with the magnetic field applied parallel to the monolayer plane using a free-beam insert placed in a resistive solenoid producing magnetic fields up to 30 T. The sample was placed on top of the x–y–z piezo stage and kept in gaseous helium at T = 4.2 K. Continuous wave and unpolarized excitation of 515 nm wavelength was used, with the excitation power of 200 µW. The excitation beam was focused and the signal was collected by the same microscope objective with a numerical aperture of NA = 0.35. The unpolarized collected light was analyzed by a 0.5-m-long monochromator equipped with a 500 lines/mm grating and CCD camera.