1 Introduction

Near future research programs in physics are strongly related to new developments in muon physics and thus mainly defined by advances in muon beams accelerator, storage and beam shaping techniques [1]. Any new peculiarities to be registered at relativistic muon interaction in a media may become a powerful instrument for various applications. However, even well-known phenomena applied at new conditions might be of growing interest for planning future experiments.

In our recent work [2] in terms of quantum electrodynamics we have modified well-known expression for Cherenkov radiation (ChR), which differs from the standard one by a factor that takes into account the media dispersion (see e.g. [3,4,5]). It was shown that a newly derived formula represents one of the limiting approximations of more general expression, while other approach in calculations leads us to another known radiation, so-called channeling radiation (CR). In those both cases the large emission angles have been considered. We have shown that the electromagnetic radiation at large angles by relativistic channeled particles appears to be a new type of optical radiation – Cherenkov-channeling radiation (ChCR).

Cherenkov radiation is one of the recognised phenomena utilised for muons registration (e.g. so-called Cherenkov detectors) [6]. Hence, hopefully ChCR could be also efficient for detecting muons, and muons states as well. The last is rather important in studying various processes, for instance, in the process of muon-to-electron transition, which takes place in a muon atom [7].

On the other hand, spontaneous CR in optical region becomes feasible in crystals with the refractive index \(n > 1\) (see [8]). It can be observed at large angles to the longitudinal momentum of a channeled particle, close to Cherenkov angle for a given projectile [9,10,11,12]. Due to the mass-factor, CR for heavy projectiles typically does not represent notable interest as a radiation source [13], while CR at large angles might be very useful for a beam diagnostics. For instance, it would be interesting to apply theoretical discoveries of this work to muon channeling, which recently was successfully realised at Fermilab muon collider [14, 15].

Previously we have shown that the intensity of ChCR by light particles, namely, by electrons, might be much higher than ChR [2]. In this report, in view of the importance for developing novel radiation-based tools, we present our results in studying both ChCR and CR by much heavier particles, i.e. by muons. We also analyse the properties of ultraviolet ChCR by relativistic muons in optically transparent crystals that takes into account the crystal dispersion \(n = n(\omega )\).

2 ChR by muons for non-zero media dispersion

As recently shown [2], channeled relativistic particles emit ChCR at the angles close to Cherenkov ones. The expression for standard ChR derived from a general definition of ChCR after the summation over all quantum states of the channeling motion in aligned crystal

$$\begin{aligned} dN_{ChR} \simeq \frac{\alpha }{c\hbar }\beta \left( 1-\frac{1}{n^2\beta ^2}\right) G(\omega ) \end{aligned}$$
(1)

under the condition \(\cos \theta _{Ch} = 1/(n\beta )\) can be simplified to the following expression

$$\begin{aligned} dN_{ChR} \simeq \frac{\alpha }{c\hbar }\beta \sin ^2\theta _{Ch} G(\omega ), \end{aligned}$$
(2)

where \(\alpha \) is the fine structure constant, c is the speed of light in a vacuum, \(\hbar \) is the reduced Planck constant, \(\beta =v/c\) is the normalised particle velocity, \(\theta _{Ch}\) is the Cherenkov angle, and \(G(\omega ) = 1 + (\omega /n(\omega ))(\partial n/\partial \omega )\) is the function, which takes into account non-zero dispersion function for a given media.

Equations (1) and (2) at \(G(\omega ) = 1\) exactly match the ChR angular distribution reduced within standard ChR theory (see e.g. [3,4,5, 16, 17]). To analyse these equations as functions of the factor \(G(\omega )\), we consider two optically transparent crystals:

  • a C (diamond) crystal as one of the most unique optically transparent crystal for the frequency range from \(\sim 220\) nm (\(\sim 5.6\) eV) up to microwaves [18];

  • a Si crystal with one of the largest value of the refractive index \(n(\omega )\) and its derivative \(\partial n(\omega )/\partial \omega \) as well [18]. Ultra thin Si crystal (thickness \(L < 1\,\upmu \)m) is optically transparent.

Figure 1 shows measured dependencies for the real part of refractive index versus photon frequency (energy) for chosen C and Si crystals [18] as well as the dependences of Cherenkov angle \(\theta _{Ch} = \arccos [1/(n\beta )]\) as function of photon energy. The use of these data in combination with Eqs. (1) and (2) allows the number of ChR photons emitted per unit path of relativistic particle to be easily calculated for a whole Cherenkov cone.

Fig. 1
figure 1

Refractive indexes \(n = n(\omega )\) (dot-solid lines related to the left ordinate) and Cherenkov angles versus ChR photons energy (solid and dashed lines related to the right ordinate) for both C and Si crystals. Dots indicate the experimental data [18]

The results of calculations for 10 GeV muons are presented in Fig. 2. To compare, the dashed curve indicates similar results obtained by the standard formula at \(G(\omega ) = 1\), i.e. without taking into account the crystal dispersion. In a Si crystal the maximum number of ChR photons, \(dN^{max}_{ChR}\simeq 1740\) ph/(eV cm) at \(\hbar \omega \simeq 3.2\) eV is about 5 times greater compared to the number of ChR photons calculated using the standard expression, while in a C crystal the maximum, \(dN^{max}_{ChR} \simeq 1010\) ph/(eV cm), observed at \(\hbar \omega \simeq 7\) eV shows near triple increase.

Fig. 2
figure 2

Spectral distributions of ChR photons per unit path and unit energy calculated for both C (diamond) and Si crystals with the dispersion taken into account (solid line) and without it (at \(G(\omega ) = 1\), dashed line)

As mentioned in [19], the fact that Cherenkov angle \(\theta _{Ch} \) is strictly dependent on the particle velocity \(\beta \), is of fundamental importance for various applications [20,21,22,23]. This reliance sets the design constrains on experimental setups for ChR registration. Figure 3 demonstrates the function \(\theta _{Ch}\) versus \(\beta \) for non-zero crystal dispersion.

Fig. 3
figure 3

Dependency of Cherenkov angles \(\theta _{Ch}\) on particle velocity \(\beta = v/c\) for both C (diamond) and Si crystals

These dependences confirm a well-known feature of Cherenkov angle \(\theta _{Ch}\) in a media with \(n(\omega )\) to tend to the limit value [19], becoming independent of the particle velocity \(\beta \). Taking into account the total internal reflection at the air-crystal interface, we can conclude that the discussed behaviour of ChR in crystals is not suitable for muon detectors in a wide range of energies (e.g. ring-imaging Cherenkov detectors). This finding is valid also for other relativistic particles.

However, the use of optically transparent crystals with high refraction indexes may drastically change situation providing various applications with reasonable ChR intensities.

3 Spontaneous CR by muons

The radiation intensity by various projectiles in crystals may reveal an essential increase for different spectral and angular intervals at specific channeling conditions [8, 10]. Typically this postulate is true for radiation emitted under small angles to the projectile momentum (“forward radiation”). However, bound channeling motion can contribute effectively to other known radiation phenomena observed at large angles [2]. This fact, mostly proved for light particles, should be also valid for relativistic muons independently of notable mass factor. Below in brief we remind properties of a new type of radiation, mixed ChCR, but applied for the case of muon planar channeling.

The number of ChCR photons within the solid angle \(\varDelta o\) emitted by a channeled muon \(\mu \) per unit length and unit photon energy during its transition \(i \rightarrow f\) between the quantum states of transverse motion (at \(i \ne f\)) can be estimated by the following expression [11, 12]

$$\begin{aligned} dN= & {} \sum _{i,f}\frac{1}{c\hbar }\frac{1}{\hbar \omega }\frac{dI_{if}}{d\omega \varDelta o}, \end{aligned}$$
(3)
$$\begin{aligned} \frac{dI_{if}}{d\omega \varDelta o}= & {} \frac{e^2 x_{if}^2\varOmega _{if}^2\omega }{4 c^3\pi \beta ^3 n^2(\omega )}P_i(\theta _0)G(\omega )\varTheta (W_{\varDelta })\varTheta (W_0)\end{aligned}$$
(4)
$$\begin{aligned}&\times \left[ \varDelta ^{+}\left( 1-2\frac{\omega }{\omega _m} + \left( \frac{\omega }{{\omega _m}}\right) ^2\right) + n^2(\omega )\beta ^2\varDelta ^{-}\right] ,\nonumber \\ \end{aligned}$$
(5)

where e is the electron charge (the charge of a muon \(\mu \)), \(x_{if} = \langle f \mid x e^{{-\mathbf {i}}\kappa _x x} \mid i \rangle \) is the matrix element for the muon \(\mu \) transition \(i \rightarrow f\) with the x-projection of the emitted photon wave vector \(k_x\), \(\hbar \varOmega _{if} = \varepsilon _i - \varepsilon _f\) is the photon energy in a rest system, \(\omega _m\) is the maximum observed radiation frequency.

Such simplification (at \(i \ne f\)) is suitable when ChCR photon energy is essentially less than the energy of channeled relativistic muon [2]. In this approximation the matrix elements \(x_{if}\) can be calculated in a dipole approximation, while the number of ChCR photons for intrazone transitions (\(i = f\)) has to be calculated utilising the expression of Eq. (1).

The solid angle \(\varDelta o\) corresponds to the photon detector elementary area, \(P_i(\theta _0)\) denotes to the initial population of the i-th transverse quantum state as a function of the incidence angle \(\theta _0 = \arctan [p_x/{p_z}]\) with respect to the channeling planes. Other notations in Eq. (5) are as follows

$$\begin{aligned} \varDelta ^{\pm }=(\varDelta \varphi \pm \cos (\varDelta \varphi + 2\varphi ) \sin \varDelta \varphi ), \end{aligned}$$
(6)

and the arguments of the Heaviside’s function \(\varTheta (\ldots )\) -

$$\begin{aligned} W_{\varDelta }= & {} \frac{\varOmega _{if}}{1 - n(\omega )\beta \cos (\theta + \varDelta \theta )} - \omega , \nonumber \\ W_0= & {} \omega - \frac{\varOmega _{if}}{1 - n(\omega )\beta \cos \theta } \end{aligned}$$
(7)

are determined by the detuning from the resonance frequencies.

4 ChCR by planar channeled GeV muons

The analysis of ChCR is performed for 10 GeV muons, \(\mu ^{-}\) and \(\mu ^{+}\), channeled along (110) planes in both Si and C (diamond) crystals. All calculations for the angular and spectral distributions of the ChCR photons number dN are carried out using the Eqs. (1), (3) and (5) for the photon beam collimator \(\varDelta \theta = \varDelta \varphi = 0.3\) mrad that is chosen as typical angular size of the area unit cell for the photon detector.

Fig. 4
figure 4

ChCR angular distributions at the energies corresponding to the maximum values of \(G(\omega )\) (see Fig. 2). ChCR photons are emitted by channeled 10 GeV muons (\(\mu ^{-}\) and \(\mu ^{+}\)): (right) \(\hbar \omega \sim 3.3\) eV in Si (\(\theta _{Ch}^{max} \simeq 81.4^{\circ }\)); (left) \(\hbar \omega \simeq 7.2\) eV in diamond (C) (\(\theta _{Ch}^{max} \simeq 73.0^{\circ }\))

The angular distributions of ChCR-photons number dN at fixed energies, i.e. \(\hbar \omega \sim 3.3\) eV in Si with \(\theta _{Ch}^{max} \simeq 81.4^{\circ }\) and \(\hbar \omega \simeq 7.2\) eV in C (diamond) with \(\theta _{Ch}^{max} \simeq 73.0^{\circ }\) that correspond to the extrema in ChR photons spectra (Fig. 2), are shown in Fig. 4. In both cases we can reveal the influence of channeling on the radiation spectra. ChCR by planar-channeled \(\mu ^{+}\) is more intense than those by \(\mu ^{-}\), and, moreover, in the spectra for \(\mu ^{+}\) we can clearly resolve a few side peaks. These results confirm the admitted difference for channeling of positively and negatively charged particles, namely, a positively charged projectile at planar channeling displays longer life time in various quantum states.

Table 1 collects values of the polar angles \(\theta \) that correspond to the maximum number \(dN^{max}\) of ChCR-photons that come to the photon detector within the angular cone \(\varDelta \theta =\varDelta \varphi =0.3\) mrad. These data prove general conclusion that the angles for maxima ChCR are slightly over those for ordinary ChR under other similar conditions.

Table 1 Maximum number of ChCR photons (hitting the photon detector) emitted by 10 GeV muons planar channeled in both C (110) and Si (110) crystals and related polar angles

Further calculations are devoted to the spectral distributions of the ChCR photon number for the polar angles \(\theta \) given in Table 1. The results are shown in Fig. 5. From the plots we can see that the ChCR spectral distribution simulated for some fixed angle forms by contribution of various photons with different energies appearing in known angular distributions (Fig. 4). This peculiarity is most pronounced for both types of muons (\(\mu ^{+}\) and \(\mu ^{-}\)) in a C crystal. The sharp maxima are observed near the radiation energy \(\hbar \omega \sim 10\) eV.

This analysis on ChCR spectral-angular distributions manifests once more that the maxima of ChCR by their photon both energies and numbers are determined by the \(G(\omega )\) extrema.

Another interesting feature of ChCR relates to the fact that the channeling mode allows expanding the range of angles, at which this type of radiation is allowed [2]. To study similar observation for muons, we have calculated the angular distributions of ChCR for fixed photon energies 7.2 and 5.6 eV emitted by positively charged muons \(\mu ^{+}\) channeled in a C crystal. Figure 6 presents those distributions for the angles greater than \(\theta > 73^{\circ }\).

Fig. 5
figure 5

Spectral distribution of ChCR by channeled 10 GeV muons (\(\mu ^{-}\) and \(\mu ^{+}\)) at the polar angle \(\theta \) (from Table 1): \(\theta \simeq 81.77^{\circ }\) (\(\mu ^{-}\)) and \(\theta \simeq 81.82^{\circ }\) (\(\mu ^{+}\)) in Si; \(\theta \simeq 73.45^{\circ }\) (\(\mu ^{-}\)) and \(\theta \simeq 73.64^{\circ }\) (\(\mu ^{+}\)) in C (diamond)

Fig. 6
figure 6

ChCR angular distributions for two chosen energies 7.2 eV and 5.6 eV emitted by 10 GeV \(\mu ^{+}\) channeled in C crystal (for angles greater than \(\theta > 73^{\circ }\))

Table 2 Maximum number of ChCR photons hitting the photon detector of the angular size \(\varDelta \theta = \varDelta \varphi = 0.3\) mrad emitted by \(\mu ^{-}\) and \(\mu ^{+}\) muons in both C and Si crystals under the azimuth angle \(\theta \)

We would like to underline here that the number of ChCR photons dN for these energies almost coincides with the number \(dN_{ChR}\) of Cherenkov photons for the same energies (see below in Table 3). However, ChR photons by quasi free projectiles are not emitted under such large angles. The maximum value of Cherenkov angle for ChR photons emitted by relativistic 10 GeV particles in C crystals equals to \(\theta ^{max}_{Ch} \simeq 73.46^{\circ }\) and corresponds to the photon energy \(\hbar \omega \simeq 11.3\) eV.

As known, a C crystal is optically transparent for photon energies \(\hbar \omega \le 5.6\) eV (220 nm), while for a ultrathin Si crystal of \(L < 1\, \mu \)m thickness we get \(T > 75 \% \) transparency for photon energies \(\hbar \omega \le 2.1\) eV (600 nm) [18]. Due to such high transmission characteristics we have plotted the angular distributions of ChCR in those crystals at fixed photon energies (Fig. 7). The maximum numbers \(dN^{max}\) of ChCR-photons hitting the photon detector with solid angle \(\varDelta o\) and corresponding polar angles are shown in Table 2, separately for C and Si crystals as well as for \(\mu ^{+}\) and \(\mu ^{-}\).

Fig. 7
figure 7

ChCR angular distributions for the frequencies of optical crystal transparency 2.07 eV for Si and 5.6 eV for C

Fig. 8
figure 8

Spectral distributions of ChCR photons emitted at polar angles \(\theta \) corresponding to the threshold of crystals optical transparency (see Table 2)

These data demonstrate again higher emission ability of \(\mu ^{+}\) with respect to \(\mu ^{-}\) that takes place due to the essential difference in the projectile phase redistribution for positively and negatively charged particles in the transverse plane of channeling. The dechanneling length for \(\mu ^{+}\) at planar channeling is over then the length for \(\mu ^{-}\) that makes the spontaneous radiation by the first one intenser [10, 13].

Successfully we have simulated the ChCR spectral distributions for fixed polar angles \(\theta \), at which the angular spectra show the maxima (see Table 2). These spectra are presented in Fig. 8. According to these results at polar angles \(\theta \) corresponding to the threshold of the crystal optical transparency, the ChCR spectral distributions by muons channeled in both C and Si crystals are qualitatively similar. However, ChCR by negative muons \(\mu ^{-}\) reveals a wide range of continuity into the visible range. Indeed, the ChCR spectral distribution at \(75.96^{\circ }\) by negative muons \(\mu ^{-}\) channeled in (110) Si (see Fig. 8a) shows the widest range of continuity into the visible range of spectrum. It may be the main practical advantage of ChCR compared to ChR. It is important to emphasize that the number of ChCR photons for this spectrum is not less than the number of ChR photons.

To compare results for ChCR and ChR photons, in Table 3 the maximum numbers \(dN_{ChR}^{max}\) and Cherenkov angles \(\theta _{Ch}\) of ChR-photons versus the energy \(\hbar \omega \) are presented. Reported data prove that the maximum numbers \(dN_{ChR}^{max}\) of ordinary ChR photons emitted by quasi free negative \(\mu ^{-}\) and positive \(\mu ^{+}\) muons crossing Si and C crystals are much smaller than the numbers \(dN^{max}\) of energy equivalent ChCR photons emitted by channeled muons (see Tables 1 and 2).

Table 3 Maximum number of ChR-photons hitting the photon detector of the angular size \(\varDelta \theta = \varDelta \varphi = 0.3\) mrad emitted by \(\mu ^{-}\) and \(\mu ^{+}\) muons in both C and Si crystals at the angles \(\theta _{Ch}^{max}\)

5 Conclusions

In this work we have presented our first results of theoretical studies on Cherenkov radiation by relativistic muons at channeling in optically transparent monocrystals, comparison of that with ordinary Cherenkov radiation under other equivalent conditions makes evident the advantages of newly proposed radiation.

We have shown for the first time that the energies (as well as the polar angles of radiation) corresponding to the largest number of photons for ChR, and ChCR as well, emitted by muons are determined by extrema of the function \(G(\omega ) = 1 + \frac{\omega }{n(\omega )}\frac{\partial n}{\partial \omega }\), which depends on the media dispersion.

Obtained results have demonstrated that the channeling process may significantly extend the angular fan for ultraviolet frequencies of photon emission by relativistic muons.

We have confirmed the fact that positively charged projectiles at channeling should radiate more intensively than negative ones, valid also for large emission angles. Indeed, due to longer life time in a channeling regime for positive \(\mu ^{+}\) muons with respect to negative \(\mu ^{-}\) ones, i.e. \(\mu ^{+}\) are characterised by larger dechanneling length compared \(\mu ^{-}\) and emit more ChCR photons. This is the result of a strong phase redistribution of the beam in transverse plane of channeling motion at Coulomb scattering in the periodic field of averaged continuous potentials formed by crystal planes or axes (in this work we have examined only planar channeling case).

In our studies we have discovered that some tail of the ChCR spectrum by relativistic muons may end into the visible frequencies. ChCR by 10 GeV negative \(\mu ^{-}\) muons planar channeled in Si (110) reveals the widest range of continuity in the visible spectrum.

And finally, we have to note another interesting phenomenon, out of the scope of this work, related to the visible spectrum of ChCR emitted by relativistic muons channeled in C crystals, studies on which will be the subject for a separate paper.