 follows from Eq.(\[eq3\]) in order to determine the numerical value of the hyperfine structure splitting $\Delta E_{hf}(MHz)$ one needs to evaluate the electron triplet density at the atomic nucleus $\rho_T(0)$. Accurate numerical evaluation of the $\rho_T(0)$ value is the main goal of this study. To perform such an evaluation we need to construct the accurate wave functions of the triplet $S-$states of the four-electron atoms and ions. In general, these wave functions are obtained as the solutions of the corresponding Schrödinger equation for the bound atomic state(s), i.e. $H \Psi = E \Psi$, where $H$ is the Hamiltonian operator defined below (see Eq.(\[Hamil\])) and $\Psi$ is the wave function and $E (< 0)$ is the total energy of the bound atomic state. This problem is considered in detail in the next Section.

Hamiltonian and bound state wave functions
==========================================

In the lowest-order approximation upon the fine-structure constant $\alpha$ we can consider the non-relativistic Schrödinger equation. The non-relativistic Hamiltonian $H$ of an arbitrary four-electron atomic system (i.e. atom, or ion) is written in the form [@LLQ] $$\begin{aligned}
   H = -\frac{\hbar^2}{2 m_e} \Bigl[ \sum^{4}_{i=1} \nabla^2_i + \frac{1}{M_A} \nabla^{2}_{5} \Bigr] - \sum^4_{i=1} \frac{Q e^2}{r_{i5}} +  \sum^{3}_{i=1} \sum^{4}_{j=2 (j>i)} 
  \frac{e^2}{r_{ij}} \label{Hamil}\end{aligned}$$ where $\hbar$ is the reduced Planck constant, $m_e$ is the electron mass, $e$ is the absolute value of the electric charge of the electron. Also, in this equation $Q e$ and $M_A$ are the electric charge and mass of the nucleus ($M_n \gg 1$) expressed in $e$ and $m_e$, respectively. Below, we consider the beryllium atom with the infinitely heavy atomic nucleus, i.e. when $M_A = \infty$, or $\frac{1}{M_A
follows from Eq.(\[eq3\ ]) in order to determine the numerical value of the hyperfine social organization separate $ \Delta E_{hf}(MHz)$ one needs to evaluate the electron three concentration at the atomic nucleus $ \rho_T(0)$. Accurate numerical evaluation of the $ \rho_T(0)$ value is the chief goal of this study. To do such an evaluation we necessitate to construct the accurate wave function of the triplet $ S-$states of the four - electron atoms and ion. In general, these wave functions are obtained as the solutions of the corresponding Schrödinger equation for the bound atomic state(s), i.e. $ hydrogen \Psi = E \Psi$, where $ H$ is the Hamiltonian operator defined below (learn Eq.(\[Hamil\ ]) ) and $ \Psi$ is the wave function and $ E (< 0)$ is the total energy of the jump atomic state. This problem is considered in detail in the future Section. 

 Hamiltonian and bound state wave functions 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

 In the lowest - order approximation upon the fine - structure constant $ \alpha$ we can consider the non - relativistic Schrödinger equation. The non - relativistic Hamiltonian $ H$ of an arbitrary four - electron atomic system (i.e. atom, or ion) is written in the form [ @LLQ ] $ $ \begin{aligned } 
    H = -\frac{\hbar^2}{2 m_e } \Bigl [ \sum^{4}_{i=1 } \nabla^2_i + \frac{1}{M_A } \nabla^{2}_{5 } \Bigr ] - \sum^4_{i=1 } \frac{Q e^2}{r_{i5 } } +   \sum^{3}_{i=1 } \sum^{4}_{j=2 (j > i) } 
   \frac{e^2}{r_{ij } } \label{Hamil}\end{aligned}$$ where $ \hbar$ is the shrink Planck constant, $ m_e$ is the electron mass, $ e$ is the absolute value of the electric charge of the electron. Also, in this equality $ Q e$ and $ M_A$ are the electric cathexis and mass of the nucleus ($ M_n \gg 1 $) expressed in $ e$ and $ m_e$, respectively. Below, we consider the beryllium atom with the infinitely heavy nuclear nucleus, i.e. when $ M_A = \infty$, or $ \frac{1}{M_A
 foplows from Eq.(\[eq3\]) in order to determine tkw numecical vzlue of ghe hyperfine structure splivtint $\Deluc E_{hf}(MHz)$ one needs tu evaluatv the elextroi triplet densitb at the atomid nucnxus $\rho_T(0)$. Accuraje numerical evaluation of tfe $\rho_T(0)$ value is the main goal of thif study. Tl perform such an tvajuatjon we need to construct the accurzte wavt functions of the triplet $S-$states of the fokr-elfctron atoms and ilns. In genetzl, eyese wave fuvctions art mbtained aa the solutions of the correspovding Schrödinget sqkdtion for tie bougd atomic state(s), i.e. $V \Psi = R \Psi$, where $H$ ls thx Haniltonian operator dehined below (see Eq.(\[Haiil\])) and $\Pvi$ is the wave funcriin ang $E (< 0)$ is rhe toual ehergy lf vhe bound afomic state. This problem is comsyeered in detajl in ehq next Section.

Hamiltonian and bound staue wabe functions
==========================================

In the lowewt-order approximation upon the fine-structure constant $\alpha$ we can consider the non-rxlxtirlfguc Schrödinger equation. The non-relativistic Hamijfomisn $H$ of an arbltrary four-electrom wtpiic system (i.e. atom, or jon) is written in hhe fori [@LLQ] $$\begin{alidned}
   H = -\frac{\hbar^2}{2 m_e} \Bigl[ \sum^{4}_{i=1} \nabla^2_i + \frab{1}{M_A} \babla^{2}_{5} \Bigr] - \sum^4_{i=1} \fxac{Q e^2}{r_{i5}} +  \sbm^{3}_{i=1} \suk^{4}_{j=2 (j>i)} 
  \frac{e^2}{r_{ij}} \label{Hamil}\eud{alighed}$$ where $\hhar$ is ths reduced Planck zonxtdnt, $m_e$ is the electron masf, $e$ is thx absplute vxlue of thq electric charna of the electron. Wlso, nn thhs equatioj $Q e$ and $M_A$ are the electric cierge and mass ox tve nuclebs ($M_n \ng 1$) expressed ig $e$ and $m_e$, reskectively. Belod, we consiser the beryllium wtom with the lnfinitely hxavy atomyc nycleys, i.e. wfdn $M_A = \infty$, ot $\frac{1}{M_A
follows from Eq.(\[eq3\]) in order to determine value the hyperfine splitting $\Delta E_{hf}(MHz)$ electron density at the nucleus $\rho_T(0)$. Accurate evaluation of the $\rho_T(0)$ value is main goal of this study. To perform such an evaluation we need to the accurate wave functions of the triplet $S-$states of the four-electron atoms and In these functions obtained as the solutions of the corresponding Schrödinger equation for the bound atomic state(s), i.e. $H = E \Psi$, where $H$ is the Hamiltonian defined below (see Eq.(\[Hamil\])) $\Psi$ is the wave function $E 0)$ is total of bound atomic state. problem is considered in detail in the next Section. Hamiltonian and bound state wave functions ========================================== In lowest-order approximation fine-structure constant we consider non-relativistic Schrödinger equation. Hamiltonian $H$ of an arbitrary four-electron atom, or ion) is written in the form $$\begin{aligned} H -\frac{\hbar^2}{2 m_e} \Bigl[ \sum^{4}_{i=1} \nabla^2_i + \nabla^{2}_{5} \Bigr] - \sum^4_{i=1} \frac{Q e^2}{r_{i5}} + \sum^{3}_{i=1} (j>i)} \frac{e^2}{r_{ij}} \label{Hamil}\end{aligned}$$ where $\hbar$ is the reduced Planck constant, $m_e$ is the electron mass, the absolute value of electric charge of electron. in equation e$ and are the electric charge and mass of the nucleus ($M_n \gg expressed in $e$ and $m_e$, respectively. Below, we consider the with infinitely heavy atomic i.e. when $M_A = or
 follows from Eq.(\[eq3\]) in order to dEtermine thE numeRicAl vAlUe of The hYperfine structURe spLitting $\Delta E_{hf}(MHz)$ one nEeds tO eVAluaTE tHe eleCtron trIPlET DenSiTy At tHe AToMic nuCleUs $\rho_T(0)$. ACcurate numEriCaL evaluation oF ThE $\rho_T(0)$ value Is tHe main goal of ThiS study. to PerFOrm suCh aN evalUation WE need tO construcT tHE accurATe wave fUNCtIons Of the triplet $S-$statES oF The four-electroN atoms AnD IoNS. in gEneRal, these waVe FunctIOns are oBTaINED as THe solutions of The correspoNDinG SchröDiNgeR EquatiOn for ThE BouNd atomic staTe(s), i.E. $H \Psi = E \Psi$, Where $H$ IS the HamILtonian OperatOr dEfiNed bELoW (sEe EQ.(\[HAMil\])) ANd $\psi$ IS thE wave funCtIoN and $E (< 0)$ Is thE TOTAl enErgY of tHe bouNd atomic state. thiS proBLem Is conSiderEd in DeTail iN the neXt SecTiOn.

Hamiltonian anD bouNd state waVe fUnCtiOnS
==========================================

In thE Lowest-OrdEr aPproximAtion upON thE fINE-StRucture constant $\alpHa$ WE CaN consideR the noN-ReLaTIvistic SChRödIngeR EQuatiOn. ThE NoN-relativIstic HAMiLtOnian $H$ oF aN arbitRaRy fOur-ElectROn atOmic syStem (i.e. atOm, or iON) is written in thE Form [@LLQ] $$\begin{aLIgNED}
   H = -\FRac{\hBar^2}{2 M_e} \Bigl[ \sum^{4}_{i=1} \nAbla^2_I + \Frac{1}{m_A} \naBLa^{2}_{5} \bigR] - \Sum^4_{i=1} \fRac{Q e^2}{R_{i5}} +  \SUm^{3}_{I=1} \Sum^{4}_{j=2 (j>i)} 
  \frac{e^2}{r_{ij}} \label{haMil}\end{AlignEd}$$ where $\hbar$ is The reduced pLANck constAnt, $m_E$ Is THe electron mass, $E$ is thE absolute vALue of the ElectRic chargE of the eleCTRon. Also, iN thIs eQuaTioN $q E$ aNd $M_A$ are the eleCTRic cHaRge and mAss Of the nuCleUs ($M_N \gg 1$) ExpReSsed in $e$ anD $m_e$, respeCtIvElY. BEloW, we coNSider the BeRylLiUm aTom wiTH the inFinitEly hEaVy ATomIc nucleUS, i.E. WHen $M_a = \iNfTy$, or $\FraC{1}{M_a
 follows from Eq.(\[eq3\])  in orderto de ter min ethenume rical value of thehyperfine structure sp litti ng $\De l ta  E_{h f}(MHz) $  o n e  ne ed stoev a lu ate t heelectro n tripletden si ty at the at o mi c nucleus$\r ho_T(0)$. Ac cur ate nu me ric a l eva lua tionof the $\rho_ T(0)$ val ue is the main go a l  o f th is study. To perf o rm such an evalua tion w en ee d tocon struct the  a ccura t e wavef un c t i ons of the triple t $S-$state s  of  the f ou r-e l ectron  atom sa ndions. In ge nera l, thesewave f u nctions are obt ainedasthe  sol u ti on s o ft hec or res p ond ing Schr öd in ger e quat i o n forthe  bou nd at omic state(s) , i .e.$ H \ Psi =  E \P si$,  w here$H$ is  theHa miltonian opera tordefined b elo w(se eEq.(\ [ Hamil\ ]))  an d $\Psi $ is th e  wa ve f u nc tion and $E (< 0)$  i s th e totalenergy of  t h e boundat omi c st a t e. Th is p r ob lem is c onside r ed  i n detai lin the  n ext  Se ction . 

Ha milton ian andbound state wave fun c tions
======= = == = = == = ==== === =========== ==== = ==== 

In th e l o west- order  a p pr o ximation upon the f in e-stru cture  constant $\a lpha$ we c a n consider  the no n -relativisticSchrö dinger equ a tion. Th e non -relativ istic Ham i l tonian $ H$ofanarb i t ra ry four-elect r o n at om ic syst em(i.e. a tom , o r i on)  i s written  in thefo rm  [ @L LQ]  $$\b e gin{alig ne d}  H= -\f r ac{\hb ar^2} {2 m _e }\ Big l[ \sum ^ {4 } _ {i=1 }\n abla ^2_ i+ \fr ac{1 } {M_ A} \nab la^{2}_{5 } \ B igr]  -  \ sum^4_{ i=1} \frac{Qe^ 2}{r_{i5}}  +   \ sum^{3 } _ {i=1} \s um^{4}_{j=2 (j>i)} 
  \ f rac{e^2 }{r _{ij} } \l abel{Hami l}\ end{al ign e d}$$ w here $ \hbar $ist h e red u c ed  Pl an ck constan t ,  $m _e$ i stheelectro n mass, $e$ is the abs olute value o f t he e l e ct ric ch a rge  o f  th e electron. Also,  in this e qu a ti on $Q e$ a n d $ M_ A$ arethe ele ctric chargeand massof the nu cl eus( $ M_n  \gg 1$) e xpressed  in $e$ a n d $m_ e $,  resp ect ively.  B elo w, we  consi d erthe b erylli um  atomwithth e infini tely heavy atomic nucle us, i. e. wh en$M_A = \i nft y $,or $\frac {1}{ M_A
 follows_from Eq.(\[eq3\])_in order to determine_the numerical_value_of the_hyperfine_structure splitting $\Delta_E_{hf}(MHz)$ one needs_to evaluate the electron_triplet density at_the_atomic nucleus $\rho_T(0)$. Accurate numerical evaluation of the $\rho_T(0)$ value is the main goal_of_this study._To_perform_such an evaluation we need_to construct the accurate wave_functions of_the triplet $S-$states of the four-electron atoms and_ions._In general, these_wave functions are obtained as the solutions of the_corresponding Schrödinger equation for the bound_atomic state(s), i.e._$H_\Psi_= E \Psi$, where_$H$ is the Hamiltonian operator defined_below (see Eq.(\[Hamil\])) and $\Psi$ is_the wave function and $E (< 0)$_is the total energy of the_bound atomic state. This problem_is considered_in detail in the next_Section.

Hamiltonian and bound_state wave_functions
==========================================

In the lowest-order_approximation upon the fine-structure constant $\alpha$_we can consider_the non-relativistic Schrödinger equation. The non-relativistic_Hamiltonian_$H$ of an_arbitrary_four-electron_atomic system_(i.e. atom, or_ion)_is written_in_the form [@LLQ] $$\begin{aligned}
  _H_= -\frac{\hbar^2}{2 m_e} \Bigl[ \sum^{4}_{i=1} \nabla^2_i +_\frac{1}{M_A} \nabla^{2}_{5} \Bigr] -_\sum^4_{i=1}_\frac{Q e^2}{r_{i5}} + _\sum^{3}_{i=1} \sum^{4}_{j=2 (j>i)} 
 _\frac{e^2}{r_{ij}} \label{Hamil}\end{aligned}$$ where $\hbar$ is the_reduced Planck_constant, $m_e$_is the electron mass, $e$ is the absolute value of the_electric charge of the electron. Also,_in this equation $Q_e$ and_$M_A$_are the electric_charge_and mass_of the nucleus ($M_n \gg 1$) expressed_in $e$_and $m_e$, respectively. Below, we consider_the beryllium atom with_the_infinitely heavy atomic nucleus, i.e. when_$M_A = \infty$, or $\frac{1}{M_A