Introduction {#Sec1}
============

Bismuth-doped glasses are the subject of great research interest due to their potential applications for the development of new fiber amplifiers and tunable laser sources. Indeed, some Bismuth-doped glasses, depending on their compositions, show remarkably wide near infrared (NIR) photoluminescence (PL) bands in the spectral range 1100--1800 nm, covering all optical telecommunication windows where the traditional rare-earths-based glasses have a rather narrow operating spectral range^[@CR1]--[@CR4]^. However, even though NIR emission has been observed in a variety of Bismuth-doped glasses, silica based^[@CR5]^, germanate^[@CR6]^, chalcogenide^[@CR7]^, fluoride^[@CR8]^, and chloride^[@CR9]^, the exact origin of this luminescence is still unknown. Because the identification of the nature of this NIR PL is crucial for further development of efficient light sources, research efforts have been devoted to the elucidation of the nature of active NIR-luminescent centre and a number of hypotheses have been proposed to explain the origin of NIR PL in Bismuth-doped materials. A detailed overview can be found in the literature^[@CR10],\ [@CR11]^. Note that in spite of the significant number of experimental and theoretical results, there is no consensus among researchers regarding the nature of NIR PL in the various Bismuth-doped glasses. Furthermore, some experimental observations appear controversial^[@CR11]^. Because the conventional optical spectroscopy not always allows to draw unambiguous conclusions, new advanced experiments are highly desirable to clarify the nature of the luminescent centres in Bismuth-doped glasses.

Recently^[@CR12]^, we used magnetic circular polarization of luminescence (MCPL) to investigate Bismuth-doped silica glass without other co-dopant. This technique enables the study of excited states (ES). Our experiments have shown that the lasing ES is a magnetic multiplet (NIR PL band 1400 nm). Furthermore, it was demonstrated that the magnetic field (MF) and temperature dependences of the MCPL signal correspond to an isolated non-Kramers doublet of the even-electron system. This important finding completely excludes any odd-electron system as a possible origin of the NIR PL. Unfortunately, based on the MCPL data alone, we were unable to reveal the magnetic multiplicity of the ground state (GS) of lasing optical centre. On the other hand, the well known red emitting centre, which was identified in all anterior investigations as Bi^2+^ ion^[@CR3],\ [@CR11],\ [@CR13]^, was not studied in our recent work^[@CR12]^. Moreover, to the best our knowledge, experiments directly proving the correctness of the assigning of red PL to the Bi^2+^ ion have so far not been presented. On the contrary, the absence of the electron paramagnetic resonance (EPR) signal in Bismuth-doped silica glass is likely to contradict this assumption. In this regard, the goal of the present paper was to clarify the spin multiplicity of the GS of both centres, the lasing one (responsible for the NIR PL) and the red PL emitting centre, by the magnetic circular dichroism (MCD).

The MCD is based on the Zeeman effect and provides valuable information about splitting of the electronic states in magnetic and crystal (CF) fields. Unlike the MCPL technique, in the MCD experiments we measure not the luminescence polarization but the difference in absorption of the right (σ^+^) and left (σ^−^) circularly polarized light in the external magnetic field. In general, MCD contains physically and mathematically separable contributions of three terms: $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${\mathscr{A}}$$\end{document}$, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ {\mathcal B} $$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${\mathscr{C}}$$\end{document}$ ^[@CR14],\ [@CR15]^. The most interesting paramagnetic $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${\mathscr{C}}$$\end{document}$-term is non-zero if the ground state of the optical centre is degenerate. This term is magnetic field and temperature dependent. The field-induced mixing of the closely lying states produces a $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ {\mathcal B} $$\end{document}$-term which depends only on the MF magnitude. The third, diamagnetic $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${\mathscr{A}}$$\end{document}$-term, is due to the excited state degeneracy and it is temperature independent. Due to these MCD term properties, this technique allows to investigate the nature of the GS of paramagnetic centre by measuring the temperature and magnetic field dependences of the MCD signal.

In the present paper, for the first time, we investigate the Bismuth-doped silica glass without other co-dopant by the MCD technique in the wavelength range 300--1600 nm. We show that GS of the lasing centre is a magnetic singlet *S* = 0. This result allows us to conclude that the attribution of the NIR PL to Bi^+^ ion in this material is not contradictory. The second very interesting and somewhat unexpected result is that the GS of centre responsible for the red PL in this glass is not just a magnetic multiplet corresponding to the half-integer effective spin, but it is also split in a zero magnetic field. The very high energy of this splitting indicates that it is due to the crystal field. It turns out that this centre is also the even-electron system. Thus, we conclude that the assignment of the red PL in this material to Bi^2+^ ion that was made in all previous works was wrong. This means that the Bi^2+^ ions are simply absent in this material.

Results and Discussion {#Sec2}
======================

In Fig. [1](#Fig1){ref-type="fig"} we show the PL spectra recorded at 1.4 K in the visible and NIR spectral regions under excitation at 375 and 450 nm. Three intense PL bands, labeled as PL1, PL2, PL3 with peaks at 1440, 830 and 585 nm, are well known and correspond to those reported earlier^[@CR16]^ (see also for the completeness)^[@CR3],\ [@CR17]^. It is worth noting that the narrow PL2 band can be observed under both excitation wavelengths, though it is much stronger and red-shifted under excitation at 375 nm, which can be explained by the inhomogeneous broadening and partial superposition of the absorption (excitation) spectra corresponding to different centres as it can be seen in Fig. [2(a,b)](#Fig2){ref-type="fig"}.Figure 1Spectra of luminescence in visible and NIR ranges under different excitation wavelengths. Figure 2(**a**) Absorption spectrum recorded at room temperature. (**b**) Excitation spectra of PL1, PL2 and PL3 luminescence bands recorded at 1.4 K. (**c**) MCD spectra recorded at various temperatures in magnetic field B = 6 T.

In Fig. [2(a)](#Fig2){ref-type="fig"} we show the absorption spectrum, recorded in the range 300--550 nm. This spectrum consists of two complex bands with maxima near 380 and 420 nm. Each of these bands can be decomposed into two sub-bands as it is suggested by the Gaussian fit also shown in Fig. [2(a)](#Fig2){ref-type="fig"}. Labels D1, B1, B2 and D2 correspond to bands with peaks at 460, 410, 375 and 350 nm, respectively. In Fig. [2(b)](#Fig2){ref-type="fig"} the PL excitation spectra recorded at 1440, 830 and 585 nm are shown. There are two bands in the excitation spectrum of the luminescence band PL2 recorded in the spectral region 300--550 nm with the intensity peaks at 415 and 375 nm. These bands can be put in the direct correspondence to the bands B1 and B2 in the Gaussian deconvolution of absorption spectrum. The excitation spectrum of the NIR PL1 band also consists of two bands centred at 370 and 420 nm. The excitation spectrum of PL3 band also reveals two bands with maxima at 465 and 350 nm. Obviously, they should be put in correspondence to the absorption bands D1 and D2. It is seen from comparison of the excitation spectra that the luminescence bands PL1 and PL2 belong to the same centre, while the band PL3 to another one.

The MCD spectra recorded at fixed MF 6 T and different temperatures in the spectral range 300--550 nm are shown in Fig. [2(c)](#Fig2){ref-type="fig"}. They consist of two bands, one positive and one negative, whose position and spectral width are very close to the absorption and excitation bands D1 and D2, respectively. In the magnetic field **B** of 6 T at 1.45 K we measured Δ*A* = −0.63 × 10^−3^ and 1.27 × 10^−3^ at 350 and 465 nm, respectively. The MCD spectra are temperature dependent, which definitely indicates the spin multiplicity of the corresponding ground state. The MCD bands exhibit the same temperature dependence as expected for two different bands of the same centre. The absence of the MCD signal in B1 and B2 absorption bands indicates that the GS of corresponding centre is a magnetic singlet with effective spin *S* = 0. It should be also noted that we could not detect any MCD signal in the whole spectral range between 550 and 1600 nm. This is not surprising taking into account previous reports^[@CR3],\ [@CR17],\ [@CR18]^ where it was shown that the NIR PL and absorption bands at 1440, 830, 375 (B1) and 415 nm (B2) belong to the same centre. This observation, in principle, does not contradict the assumption that the Bi^+^ ion is responsible for the NIR PL in this material^[@CR19]^.

The most valuable information about the MCD active centre can be obtained from the analysis of saturation of magnetization. Magnetic field dependences of MCD were recorded at different fixed temperatures in two spectral points corresponding to the band peaks: 465 and 352 nm (bands D1 and D2, respectively). Because qualitatively two measured curve subsets were identical, in Fig. [3(a)](#Fig3){ref-type="fig"} we show only the one corresponding to the band D1. It is seen that the isotherms do not superimpose but exhibit a small degree of nesting. The curves tend towards the saturation at high MF and low temperature. Temperature dependences of MCD recorded at fixed magnitudes of the external magnetic field plotted as a function of 1/*T* are shown in Fig. [3(b)](#Fig3){ref-type="fig"}. The MCD signal increases with decreasing temperature and it becomes saturated at high magnetic field and low temperature.Figure 3Magnetic field (**a**) and temperature (**b**) dependences of MCD at 465 nm (MCD and absorption band D1). Solid lines in (**b**) represent the fit by Eq. ([2](#Equ2){ref-type=""}).

The analysis of saturation of magnetization data was performed in terms of the spin Hamiltonian with axial (*D*) and rhombic (*E*) components of the crystal field^[@CR20],\ [@CR21]^:$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$ {\mathcal H} =D\{{S}_{z}^{2}-\frac{1}{3}S(S+\mathrm{1)}\}+E({S}_{x}^{2}-{S}_{y}^{2})+{g}_{\parallel }{\mu }_{B}B{S}_{z}\,\cos \,\theta ,$$\end{document}$$where *S* is the effective spin, *θ* is the angle between the applied magnetic field *B* and the symmetry axis (z-axis), g~\|\|~ is the axial component of g-tensor and *μ* ~*B*~ is the Bohr magneton. In the frame of this model, the odd- and even-electron systems exhibit qualitatively different behaviour. In particular, the crystal field cannot remove completely the degeneracy in the case of a half-integer effective spin *S* = 1/2, 3/2, ... (Kramers doublet), as it is shown in Fig. [4(a)](#Fig4){ref-type="fig"} for the simplest case of *S* = 1/2. Kramers doublet can be split only by the MF, and when it is applied parallel to z-axis (*θ* = 0), the gap between spin sub-levels increases linearly with MF. In the case of integer spin, the crystal field potential can completely remove the degeneracy of the state as it is shown in Fig. [4(b)](#Fig4){ref-type="fig"} for the system with *S* = 1. In order to reach the lowermost energy position for the non-Kramers doublet *M* ~*s*~ = ±1, the axial component of the crystal field must be negative (*D* \< 0). The rhombic component of the crystal field removes the residual degeneracy, producing the so-called zero-field splitting (ZFS) with corresponding energy gap between *M* ~*s*~ = ±1 sub-levels Δ = 2\|E\|. MCD can be observed only in the presence of MF because of the peculiarities of \|±〉 wave functions^[@CR21]^.Figure 4Energy level diagrams corresponding to (**a**) *S* = 1/2 and (**b**) *S* = 1 effective spin in the assumption of the Hamiltonian (1).

Quantitative information about the studied system was obtained by fitting the saturation experimental data to the following equation^[@CR21],\ [@CR22]^:$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${\rm{\Delta }}A={A}^{sat}{\int }_{0}^{1}\frac{2g{\mu }_{B}B\,{\cos }^{2}\,\theta }{\sqrt{{{\rm{\Delta }}}^{2}+{\mathrm{(2}g{\mu }_{B}B\cos \theta )}^{2}}}\times \,\tanh (\frac{\sqrt{{{\rm{\Delta }}}^{2}+{\mathrm{(2}g{\mu }_{B}B\cos \theta )}^{2}}}{2kT})d(\cos \,\theta ),$$\end{document}$$where *A* ^*sat*^ is the saturation value of MCD and *g* is the effective *g*-factor, which depends on the value of the effective spin *S*. The simultaneous fit of saturation curves shown in Fig. [3(b)](#Fig3){ref-type="fig"} results in low *g*-value *g*  *=* 1.08 ± 0.02 and Δ = 1.22 ± 0.11 cm^−1^ in the assumption that the GS is triplet, *S* = 1. The non-zero value of ZFS energy Δ definitely indicates that the optical centre responsible for the observed MCD in Bismuth-doped silica glass is a system with even number of electrons. As a result, any odd-electron system, and Bi^2+^ ions in particular, should be excluded as possible origin of the MCD bands. Moreover, the even-electron systems with a non-degenerate GS, like Bi^5+^ and Bi^+^ ions, also cannot cause the MCD. Amorphous silica glass has a very large family of the point defects^[@CR23],\ [@CR24]^. However, there are no suitable among the known defects, which would have a triplet ground state and the absorption bands in the range 300--600 nm. In this regard, we put forward the hypothesis that the MCD bands Di and red luminescence are related to the unknown defect in the glass network. Probably, Bismuth doping leads to the formation of some new kind of defect or it changes the properties of the existing one. We would also like to note that similar conclusions about the nature of the red luminescence were made by us in the detailed magneto-optical investigation of Bismuth-doped aluminosilicate glass^[@CR25]^.

Our experiments show that the observation of the conventional EPR in Bismuth-doped silica glass is possible. However, attempts to detect with conventional spectrometers an EPR signal, that could be assigned to a Bismuth-related centre, failed^[@CR10]^. The first and main reason for this is that all the EPR experiments, at the best of our knowledge, were performed with the microwave frequencies in the X- band (9.5 ± 0.5 GHz), which is too low to match measured here $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${\rm{ZFS}}\simeq 36.6$$\end{document}$ GHz. The second reason, but much less important, is that the microwave transitions $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$|-\rangle \iff |+\rangle $$\end{document}$ are allowed in the configuration when the microwave MF is parallel to the static magnetic field **B**, whereas in the most conventional EPR spectrometers these fields are perpendicular to each other. Finally, in the glass materials all the paramagnetic centres are randomly oriented. This causes an extremely large broadening of the resonance line so that its intensity may be below the limit of sensitivity of the spectrometer.

Conclusion {#Sec3}
==========

In conclusion, magnetic circular dichroism has been studied in a Bismuth-doped silica glass without other co-dopant. It was shown that the ground state of the optical centre responsible for the red photoluminescence is the spin multiplet. Magnetic field and temperature dependences of MCD revealed the zero-field splitting in the GS. This fact implies that the red luminescence is caused by the even-electron system. This interpretation, in turn, denies the assumption that the Bi^2+^ ions are the source of red luminescence in the studied material. Thus, the magneto-optical methods do not reveal any trace of the Bi^2+^ ion in the studied material. Also, our investigation made prediction of the experimental conditions for observing the conventional EPR signal. Finally, we conclude that the lasing centre responsible for the NIR PL has a non-degenerate GS (effective spin *S* = 0), that, in principle, does not contradict the hypothesis that this photoluminescence is caused by the electronic transitions in Bi^+^ ions.

Materials and Methods {#Sec4}
---------------------

The investigated bulk SiO~2~:Bi samples with the dimensions 2 mm × 4 mm × 5 mm were similar to that reported by us previously^[@CR16]^. The only difference was the enhanced atomic ratio Bi/Si, which was estimated ≈400 ppm using the electron probe microanalysis (EPMA).

In Fig. [5](#Fig5){ref-type="fig"} we show the setup of the MCD spectrometer. The sample was placed in a helium closed-cycle magneto-optical cryostat SpectromagPT (Oxford Instruments). The system provides magnetic field up to 7 T with the homogeneity of 0.01% over 10 mm diameter spherical volume and temperature stabilization in 1.4--300 K range. As a light source we used a quartz tungsten or UV enhanced Xenon lamp. Light from the source was dispersed by a monochromator with a spectral resolution 2.5 and 5 nm in the range 350--900 and 900--1600 nm, respectively. After the collimating lens and exchangeable long wave pass filter (LWPF) the light was linearly polarized with a large aperture Glan-Thompson polarizer. The photoelastic modulator (PEM, II/FS20 optical head, Hinds Instruments) acted as a dynamic wave plate oriented at 45° to the polarization plane of the incident light beam. In order to transform the linear polarization of light to the circular, σ^+^ or σ^−^, PEM periodically introduced the retardation $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\delta ={\delta }_{0}\,\sin \,\mathrm{(2}\pi {f}_{PEM}t)$$\end{document}$, where *δ* ~0~ = *λ*/4 is the peak retardation, *f* ~*PEM*~ = 20 KHz is the modulator frequency. Then the light beam was chopped at low frequency $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$${f}_{CH}\simeq 300$$\end{document}$ Hz. After passing through the sample it was detected by the amplified photodetector (Si or Ge photodiode). In the case of a small dichroism Δ*A* \< 0.1 the MCD signal is Δ*A* = *V* ~*PEM*~/*V* ~*CH*~, where *V* ~*PEM*~ and *V* ~*CH*~ are the signals at PEM's and chopper's frequencies, respectively. The phase-sensitive technique with a multi-frequency lock-in amplifier (MFLI, Zurich Instruments) was used to record *V* ~*PEM*~ and *V* ~*CH*~ signals simultaneously.Figure 5Setup for MCD measurements. L1/L2 - UV Xenon or quartz tungsten lamp, LWPF - long wave pass filter; GTP - Glan-Thompson polarizer; PEM - photoelastic modulator; CH - chopper; SpMag - cryomagnetic system; PD - Si or Ge amplified photo-detector; LIA - multi-frequency lock-in amplifier.

**Publisher\'s note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

O.V.L. and I.M.R. were supported by "Agence Nationale de la Recherche" via ANR grant "BOATS" 12BS04-0019-01, and "Fonds Européen de Développement Economique des Régions" (FEDER) through the "Contrat de Projets Etat Région" (CPER Photonics for Society P4S) 2015--2020.

O.V.L. conducted the experiments, analysed the data. H.E.H. and M.B. prepared samples. V.B.A. synthesized the precursors and discussed the results with I.M.R. and O.V.L. I.M.R. conceived the experiment and wrote the manuscript. All authors reviewed the manuscript.

Competing Interests {#FPar1}
===================

The authors declare that they have no competing interests.
