 charge unit ($\alpha=e^2/(4\pi)$, $e<0$) are used in the paper.

Hyperfine splitting in muonic atoms {#s:ked}
===================================

The ground-state hyperfine splitting in muonic atoms can be written in the form: $$\label{E:STS}
\varDelta E=\varDelta E_{\mathrm{NS}}+\varDelta E_{\mathrm{BW}}+\varDelta E_{\mathrm{QED}},$$ where $\varDelta E_{\mathrm{NS}}$ is the hyperfine splitting value incorporating the relativistic and nuclear charge distribution (“nuclear size”) effects, $\varDelta E_{\mathrm{BW}}$ is the BW contribution, $\varDelta E_{\mathrm{QED}}$ is the QED correction. The $\varDelta E_{\mathrm{NS}}$ value can be calculated by the formula: $$\label{E:NS}
\varDelta E_{\mathrm{NS}}=-\alpha\frac{4}{3}\frac{\mu}{\mu_{N}}\frac{1}{m_p}\frac{(2I+1)}{2I}
\int\limits_{0}^{\infty}dr \: g(r)f(r),$$ where $\alpha$ is the fine structure constant, $\mu$ is the nuclear magnetic moment, $\mu_N$ is the nuclear magneton, $m_p$ is the proton mass, and $I$ is the nuclear spin. $g(r)$ and $f(r)$ are the radial parts of the Dirac wave function: $$\Psi(\textbf{r})=
\begin{pmatrix}
    g(r)\varOmega_{\kappa m}(\textbf{n}) \\
    \mathit{i}f(r)\varOmega_{-\kappa m}(\textbf{n})
\end{pmatrix},$$ which are determined by solving the Dirac equation with the Fermi distribution of the nuclear charge ($4\pi\int dr r^{2}\rho(r)=1$): $$\rho (r)=\frac{\rho_{0}}{1+\mathrm{exp}(\frac{r-c}{a})}.$$ Here $c$ is the half-density radius and $a$ is related to the skin thickness $t$ by $t=(4\log 3)a$, defined as the distance over which the charge density falls from 90% to 10% of its maximum value.

The individual contributions to $\varDelta E$ for muonic atoms of $^{203}$Tl,
charge unit ($ \alpha = e^2/(4\pi)$, $ e<0 $) are used in the paper. 

 Hyperfine splitting in muonic atom { # south: ked } 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

 The ground - state hyperfine splitting in muonic atoms can be write in the form: $ $ \label{E: STS } 
 \varDelta E=\varDelta E_{\mathrm{NS}}+\varDelta E_{\mathrm{BW}}+\varDelta E_{\mathrm{QED}},$$ where $ \varDelta E_{\mathrm{NS}}$ is the hyperfine splitting value incorporating the relativistic and nuclear bang distribution (“ nuclear size ”) effect, $ \varDelta E_{\mathrm{BW}}$ is the BW contribution, $ \varDelta E_{\mathrm{QED}}$ is the QED correction. The $ \varDelta E_{\mathrm{NS}}$ value can be calculate by the formula: $ $ \label{E: NS } 
 \varDelta E_{\mathrm{NS}}=-\alpha\frac{4}{3}\frac{\mu}{\mu_{N}}\frac{1}{m_p}\frac{(2I+1)}{2I } 
 \int\limits_{0}^{\infty}dr \: g(r)f(r),$$ where $ \alpha$ is the fine social organization constant, $ \mu$ is the nuclear magnetic moment, $ \mu_N$ is the nuclear magneton, $ m_p$ is the proton mass, and $ I$ is the nuclear spin. $ g(r)$ and $ f(r)$ are the radial character of the Dirac wave function: $ $ \Psi(\textbf{r})= 
 \begin{pmatrix } 
     g(r)\varOmega_{\kappa m}(\textbf{n }) \\ 
     \mathit{i}f(r)\varOmega_{-\kappa m}(\textbf{n }) 
 \end{pmatrix},$$ which are determined by solving the Dirac equation with the Fermi distribution of the nuclear charge ($ 4\pi\int dr r^{2}\rho(r)=1 $ ): $ $ \rho (r)=\frac{\rho_{0}}{1+\mathrm{exp}(\frac{r - c}{a})}.$$ Here $ c$ is the half - concentration radius and $ a$ is related to the skin thickness $ t$ by $ t=(4\log 3)a$, define as the distance over which the charge density fall from 90% to 10% of its maximum value. 

 The individual contribution to $ \varDelta E$ for muonic atoms of $ ^{203}$Tl,
 chwrge unit ($\alpha=e^2/(4\pi)$, $e<0$) are used in the pakee.

Hyperhine spmitting kn muonic atoms {#s:ked}
===================================

The grouid-stqte htperfine splitting in ouonic atlms can ve wcitten in the focj: $$\label{C:FTS}
\vzvDeltc X=\varDelta E_{\mathtm{NS}}+\varDelta E_{\mathrm{BW}}+\varDalga E_{\mathrm{QED}},$$ where $\varDelta E_{\mathrm{NF}}$ is thr jyperfine splijting dalus incorporating the relativistic ahd nucltar charge distribition (“nuclear size”) effects, $\varFelta E_{\mathrm{BW}}$ is the BW conjdibtrion, $\varDeltx E_{\mathrm{QTD}}$ is the QES correction. The $\varDelta E_{\mathfm{NS}}$ ralue can bg dapwulated by vhe fogmula: $$\label{E:NS}
\varDeltd E_{\mathtm{NS}}=-\alpha\frac{4}{3}\fvac{\mu}{\ku_{N}}\drac{1}{m_p}\frac{(2I+1)}{2I}
\int\limits_{0}^{\mnfty}dr \: g(r)f(r),$$ where $\wlpha$ is dhz fine structure consrabt, $\mu$ is dhe vycldar megnstic mlmeit, $\mu_N$ is tge nuclear nagneton, $m_p$ is the kroein mass, and $I$ is thq guclear spin. $g(r)$ and $f(r)$ are the radial pdrta of the Dirac wave funxtion: $$\Psi(\textbf{r})=
\begin{kmatrix}
    g(r)\varOmega_{\kappa m}(\textbf{n}) \\
    \mathit{i}f(r)\varOmega_{-\ka[pa m}(\vebtby{k})
\cnd{poqtgix},$$ which are determined by solving the Dirac sqianion with the Ferii distribuyiln jf the nucleat chargz ($4\pj\int dr r^{2}\rho(r)=1$): $$\rho (g)=\frac{\rhj_{0}}{1+\mathem{exp}(\frac{w-c}{a})}.$$ Nere $c$ is the half-density rqdius and $a$ ps rwlated to the skin thickness $c$ by $t=(4\kog 3)a$, defined as the distancz over which the fharge dehrity falls from 90% to 10% mf its maximum value.

The inqividual rontrnbutions to $\farDelea E$ for mkonic atoms of $^{203}$Tl,
charge unit ($\alpha=e^2/(4\pi)$, $e<0$) are used in Hyperfine in muonic {#s:ked} =================================== The atoms be written in form: $$\label{E:STS} \varDelta E_{\mathrm{NS}}+\varDelta E_{\mathrm{BW}}+\varDelta E_{\mathrm{QED}},$$ where $\varDelta E_{\mathrm{NS}}$ the hyperfine splitting value incorporating the relativistic and nuclear charge distribution (“nuclear size”) $\varDelta E_{\mathrm{BW}}$ is the BW contribution, $\varDelta E_{\mathrm{QED}}$ is the QED correction. The E_{\mathrm{NS}}$ can calculated the formula: $$\label{E:NS} \varDelta E_{\mathrm{NS}}=-\alpha\frac{4}{3}\frac{\mu}{\mu_{N}}\frac{1}{m_p}\frac{(2I+1)}{2I} \int\limits_{0}^{\infty}dr \: g(r)f(r),$$ where $\alpha$ is the fine structure constant, $\mu$ the nuclear magnetic moment, $\mu_N$ is the nuclear $m_p$ is the proton and $I$ is the nuclear $g(r)$ $f(r)$ are radial of Dirac wave function: \begin{pmatrix} g(r)\varOmega_{\kappa m}(\textbf{n}) \\ \mathit{i}f(r)\varOmega_{-\kappa m}(\textbf{n}) \end{pmatrix},$$ which are determined by solving the Dirac equation with the distribution of charge ($4\pi\int r^{2}\rho(r)=1$): (r)=\frac{\rho_{0}}{1+\mathrm{exp}(\frac{r-c}{a})}.$$ $c$ is the and $a$ is related to the by $t=(4\log 3)a$, defined as the distance over the charge falls from 90% to 10% of maximum value. The individual contributions to $\varDelta E$ muonic atoms of $^{203}$Tl,
 charge unit ($\alpha=e^2/(4\pi)$, $e<0$) are useD in the papeR.

HypeRfiNe sPlIttiNg in Muonic atoms {#s:keD}
===================================

the gRound-state hyperfine splIttinG iN MuonIC aToms cAn be wriTTeN IN thE fOrM: $$\laBeL{e:StS}
\vardelTa E=\varDElta E_{\mathrM{NS}}+\VaRDelta E_{\mathrM{bW}}+\VarDelta E_{\mAthRm{QED}},$$ where $\vaRDeLta E_{\maThRm{Ns}}$ Is the HypErfinE splitTIng valUe incorpoRaTIng the RElativiSTIc And nUclear charge distrIBuTIon (“nuclear size”) EffectS, $\vARDELTa E_{\MatHrm{BW}}$ is the bW ContrIBution, $\vARDELTA E_{\mAThrm{QED}}$ is the QeD correctioN. the $\VarDelTa e_{\maTHrm{NS}}$ vAlue cAn BE caLculated by tHe foRmula: $$\labeL{E:NS}
\vaRdelta E_{\mAThrm{NS}}=-\aLpha\frAc{4}{3}\fRac{\Mu}{\mu_{n}}\FrAc{1}{M_p}\fRaC{(2i+1)}{2I}
\iNT\lImiTS_{0}^{\inFty}dr \: g(r)f(R),$$ wHeRe $\alpHa$ is THE FIne sTruCturE consTant, $\mu$ is the nuCleAr maGNetIc momEnt, $\mu_n$ is tHe NucleAr magnEton, $m_P$ iS the proton mass, aNd $I$ iS the nucleAr sPiN. $g(r)$ AnD $f(r)$ arE The radIal ParTs of the dirac waVE fuNcTION: $$\PSi(\textbf{r})=
\begin{pmatRiX}
    G(R)\vArOmega_{\kAppa m}(\tEXtBf{N}) \\
    \Mathit{i}f(R)\vArOMega_{-\KAPpa m}(\tExtbF{N})
\eNd{pmatriX},$$ which ARe DeTermineD bY solviNg The dirAc equATion With thE Fermi diStribUTion of the nucleAR charge ($4\pi\int dR R^{2}\rHO(R)=1$): $$\rHO (r)=\frAc{\rHo_{0}}{1+\mathrm{exp}(\Frac{R-C}{a})}.$$ HeRe $c$ iS ThE haLF-densIty raDiUS aND $a$ is related to the skiN tHickneSs $t$ by $T=(4\log 3)a$, defined aS the distanCE OVer which The cHArGE density falls fRom 90% to 10% Of its maximUM value.

ThE indiVidual coNtributioNS To $\varDelTa E$ For MuoNic ATOmS of $^{203}$Tl,
 charge unit ($\alpha=e^2/ (4\pi)$, $ e<0$)  ar e u se d in  the  paper.

Hyper f inesplitting in muonic at oms { #s : ked} == ===== ======= = == = = === == == === == = == =

Th e g round-s tate hyper fin esplitting in mu onic atoms  ca n be written  in  the f or m:$ $\lab el{ E:STS }
\var D elta E =\varDelt aE _{\mat h rm{NS}} + \ va rDel ta E_{\mathrm{BW} } +\ v arDelta E_{\ma thrm{Q ED } }, $ $  wh ere  $\varDelt aE_{\m a thrm{NS } }$ i s  th e  hyperfine sp litting val u e i ncorpo ra tin g  the r elati vi s tic  and nuclea r ch arge dist ributi o n (“nuc l ear siz e”) ef fec ts,  $\v a rD el taE_ { \ma t hr m{B W }}$  is theBW  c ontri buti o n , $\va rDe ltaE_{\m athrm{QED}}$istheQ EDcorre ction . Th e$\var DeltaE_{\m at hrm{NS}}$ value  can  be calcu lat ed  by  t he fo r mula:$$\ lab el{E:NS }
\varD e lta  E _ { \ ma thrm{NS}}=-\alpha\ fr a c {4 }{3}\fra c{\mu} { \m u_ { N}}\frac {1 }{m _p}\ f r ac{(2 I+1) } {2 I}
\int\ limits _ {0 }^ {\infty }d r \: g (r )f( r), $$ wh e re $ \alpha $ is the  fine structure cons t ant, $\mu$ is th e nu c lear  ma gnetic mome nt,$ \mu_ N$ i s  t hen uclea r mag ne t on ,  $m_p$ is the proto nmass,and $ I$ is the nuc lear spin. $ g (r)$ and  $f( r )$ are the radial  part s of the D i rac wave  func tion: $$ \Psi(\tex t b f{r})=
\ beg in{ pma tri x } 
   g(r)\varOm e g a_{\ ka ppa m}( \te xtbf{n} ) \ \
  \ma th it{i}f(r) \varOmeg a_ {- \k ap pam}(\t e xtbf{n}) 
\ end {p mat rix}, $ $ whic h are  det er mi n edby solv i ng t he D ir ac  equ ati on  with  the Fer mi dist ributionoft he n uc le ar char ge ($4\pi\int  d r r^{2}\rh o( r)= 1$): $ $ \ rho (r)= \frac{\rho_{0}}{1+\math r m{exp}( \fr ac{r- c}{a })}.$$ He re$c$ is  th e  half- densit y rad iu s a n d  $a$i s  r ela te d to the s k i n t hickn es s $t $ by $t =(4\log 3)a$, defi n edas the distan ceover w hi cht he cha rg e  de n s ity falls from90% to 10%  o f  i ts maximum val ue .

Theindivid ual c o ntribut ions to $ \varDelta  E $ fo r muo nic atomsof $^{20 3}$Tl,
 charge_unit ($\alpha=e^2/(4\pi)$,_$e<0$) are used in_the paper.

Hyperfine_splitting_in muonic_atoms_{#s:ked}
===================================

The ground-state hyperfine_splitting in muonic_atoms can be written_in the form:_$$\label{E:STS}
\varDelta_E=\varDelta E_{\mathrm{NS}}+\varDelta E_{\mathrm{BW}}+\varDelta E_{\mathrm{QED}},$$ where $\varDelta E_{\mathrm{NS}}$ is the hyperfine splitting value incorporating the_relativistic_and nuclear_charge_distribution_(“nuclear size”) effects, $\varDelta E_{\mathrm{BW}}$_is the BW contribution, $\varDelta_E_{\mathrm{QED}}$ is_the QED correction. The $\varDelta E_{\mathrm{NS}}$ value can_be_calculated by the_formula: $$\label{E:NS}
\varDelta E_{\mathrm{NS}}=-\alpha\frac{4}{3}\frac{\mu}{\mu_{N}}\frac{1}{m_p}\frac{(2I+1)}{2I}
\int\limits_{0}^{\infty}dr \: g(r)f(r),$$ where $\alpha$ is the_fine structure constant, $\mu$ is the_nuclear magnetic moment,_$\mu_N$_is_the nuclear magneton, $m_p$_is the proton mass, and $I$_is the nuclear spin. $g(r)$ and_$f(r)$ are the radial parts of the_Dirac wave function: $$\Psi(\textbf{r})=
\begin{pmatrix}
  _ g(r)\varOmega_{\kappa m}(\textbf{n}) \\
 _ _\mathit{i}f(r)\varOmega_{-\kappa m}(\textbf{n})
\end{pmatrix},$$ which are determined_by solving the_Dirac equation_with the Fermi_distribution of the nuclear charge ($4\pi\int_dr r^{2}\rho(r)=1$): $$\rho_(r)=\frac{\rho_{0}}{1+\mathrm{exp}(\frac{r-c}{a})}.$$ Here $c$ is the half-density_radius_and $a$ is_related_to_the skin_thickness $t$ by_$t=(4\log_3)a$, defined_as_the distance over which the charge_density_falls from 90% to 10% of its_maximum value.

The individual contributions_to_$\varDelta E$ for muonic_atoms of $^{203}$Tl,