Ground State, Isoelectronic Ions and Low-Lying Excited States of Lithium Atom in Strong Magnetic Field

In the framework of the variational Monte Carlo method, the ground states of the lithium atom and lithium like ions up to Z = 10 in an external strong magnetic field are evaluated. Furthermore, the two low-lying excited states 2 1 2 s s , 2 1 1 2 s p − and 1 2 1 2 3 s p d − − of the lithium atom in strong magnetic field are also investigated. Simple trial wave functions for lithium are used.


Introduction
The properties of atoms in strong magnetic fields are relevant to several domains such as astrophysics, condensed matter physics, and plasma physics [1].In the past, only hydrogen and helium atoms have been studied for magnetic fields of arbitrary strength [2].The influence of a magnetic field on the properties of atoms and molecules is of great interest.The study of strongly magnetized atoms was first treated in a solid-state context.In semiconductors, the absorption of photons can promote electrons from the valence to the conduction band, forming an electron-hole pair.The absorption spectra of these excitons, when subjected to a magnetic field, display Zeeman effects [3].The methods of treating the resulting problems arising from these topics have resulted in few studies of many-body effects in strong magnetic fields, i.e., in the field regime where perturbation theory is no longer applicable.
The motivation in this area arises from several sources.On the one hand, this is due to the astrophysical discovery of strong magnetic fields on the surfaces of white dwarfs (10 2 -10 5 T) and neutron stars (10 7 -10 9 T) [4] [5].Also, the complex properties of atoms under these extreme conditions are of immediate interest from a purely theoretical point of view.In addition, the observations of excitons with small effective masses and large dielectric constants in semiconductors, which result in very large effective magnetic fields, give additional impetus for this subject [6].
It is well-known that the study of the electronic structure of atoms in the external magnetic field is a quite complicated subject, especially for the strong magnetic field.In a weak magnetic field, which can be treated as a perturbation, the wave function can be expanded in terms of spherical harmonic functions.In a very strong magnetic field, where the Lorentz forces dominate, the wave function can be expanded in terms of Landau-like orbitals.But in the strong magnetic field, where the Lorentz and Coulomb forces are of nearly equal importance, neither of them can be treated as a perturbation [7].Hence, it is necessary to develop non-perturbative techniques to solve this problem.The parameter z β can be used to characterize three different regimes of strength: the low (weak, , and Z is the charge of the atomic nucleus.So far, considerable effort has been devoted to the theoretical investigations of atoms in magnetic fields with arbitrary strength.
Rӧsner et al. [8] have done detailed work on the spectrum of the hydrogen atom in magnetic fields up to 4.7 × 10 8 T, which has been successfully applied to the identification of the observed spectra from many magnetic white dwarfs [9] [10].
For the hydrogen atom in magnetic fields up to 9.4 × 10 8 T, Kravchenko et al.
[11] [12] presented accurate results with an accuracy of 10 −12 .In view of this, the issue of the hydrogen atom over a wide range of magnetic fields can be completely solved.Despite this great success, there were also many magnetic white dwarfs whose spectra remain unexplained or could not be completely accounted for with hydrogen atom.It is thus essential to acquire extensive and accurate data of energy levels and resulting transition wavelengths of multi-electron atoms subjected to strong magnetic fields.Hence, the research emphasis of atoms under these extreme conditions focused on the next lightest element, neutral helium.In contrast to the hydrogen atom, the problem of the helium atom is much more intricate because of the occurrence of the electron-electron repulsion.Even so, there exist many investigations on helium atoms in the literature.Thurner et al. [13] gave the energies of several triplets excited states of helium atoms for magnetic field strengths in the range 4.7 × 10 1 to 4.7 × 10 8 T. Furthermore, they presented a consistent correspondence diagram between the low field and high field domains for the first time.
Many studies have been recently presented by Doma et al. using the Variational Monte Carlo (VMC) method for the helium, lithium, and beryllium atoms [14] [15] [16] [17] [18].The success of the VMC method extended also to the light molecules, where similar investigations have been carried out for the hydrogen molecule and its molecular ion by Doma et al. [19] [20] [21].
In fact, the problem of highly accurate calculations of the bound states in the three-electron atomic systems has attracted continuing attention.For lithium atom and its ions up to Z = 10, in the strong magnetic field and using the VMC method, Doma et al. calculated the energies of these systems [16].
In the present paper, we applied the VMC method and constructed three trial wave functions, which take into consideration the correlation functions due to the interactions between the electrons in the atom, to calculate the ground state of the lithium atom and its ions up to Z = 10 in the presence of an external magnetic field regime with 0 ~100 γ = a.u.Furthermore, the two low-lying excited states of the lithium atom in an external magnetic field is also solved.For this purpose, simple trial wave functions with low number of variational parameters for lithium are used.To find the minimum energy eigenvalues we applied the VMC method.We first constructed the Hamiltonian operator for the three-electron system in the absence of the external field.We made use of the assumption of an infinitely heavy nucleus in the unrestricted Hrtree-Fock (HF) approximation.Hence, the nonrelativistic Hamiltonian for the lithium atom, in the absence of the field, under the Born-Oppenheimer approximation of zero order, that is, with the Li nucleus assumed to be of infinite mass, is (in Hartree atomic units) [22] given by: where i ∇ is the 3-vector of the momentum of the ith electron, Z is the nuclear charge (here, Z = 3), i r is the distance between the ith electron and the Li nucleus, and ij r are the interelectron distances.In our calculations we used the form of H in Hylleraas Coordinates [22] [23] as follows: According to the assumption that the nuclear mass is infinite, and the magnetic field is oriented along the z-axis, the non-relativistic Hamiltonian  for the Li atom in the presence of the magnetic field (in atomic units (a.u.)) takes the form [7]: where x y x y x y ρ ρ ρ ρ are the z-components of the total angular momentum and spin respectively, The VMC method is used to compute quantum expectation values of an operator with a given trial wave function.This method is a combination of the will-known variational principle and the Monte Carlo technique of evaluating multi-dimensional integrals.The variational principle states that the expectation value of the Hamiltonian H with respect to the trial wave function T ψ is the varitional energy [24]: where R is the 3N-dimensional vector of the electron coordinates and exact E is the exact value of the energy of a given state.To evaluate VMC E , Equation (2.4), we rewrite it as: where is positive everywhere and interpreted as a probability distribution function and Also, it is important to calculate the standard deviation of the energy [17] [18]: ( ) where L is the ensemble size of random numbers { } , which may be generated using a variety of methods and N is the number of ensembles.These ensembles so generated must reflect the distribution function itself.Open Journal of Microphysics At each of the points { } i  , the weighted average ( ) ter enough evaluations, the VMC estimate of VMC E will be: ( )

The Trial Wave Functions
In the VMC method, a good trial wave function, For the ground state of the lithium atom, 2 1 2 s s and its isoelectronic ions, in the presence of weak magnetic field, we used the following trial wave function [14]: where In the last equation, Z ′ is an effective charge, which is treated as a variational parameter, to minimize the energy.It is seen that 0 ψ is a linear combination of terms constructed from the single-particle hydrogenic wave functions in such a way that the total spin of the system is 1 2 S = and its z-component is 1 2 In Equation (3.2), α is a spin up function and β is a spin down function.
For the correlation function f, we considered the Jastrow correlation function which is a product of functions of the form [15]: It seems very reasonable to enforce the cusp conditions on trial wave functions because they are obeyed by the exact wave function.These constraints lead to smoother behavior in the local energy,

( ) ( ) ( )
. When two particles interacting via the Coulomb potential approach one another, the potential ener-gy diverges, and therefore the exact wave function must have a cusp so that the local kinetic energy, ( ) ( ) ( ) , supplies an equal and opposite divergences.For the electron-nuclear cusp this condition is where Z is the nuclear charge and i r is any single electron-nucleus coordinate.
Imposition of the cusp conditions is in fact very important in VMC calculations because divergences in the local energy led to poor statistical behavior and even instabilities in the calculations.In order that the function f satisfies the cusp conditions we must have: Hence, the ground state trial wave 1 ψ has two variational parameters, namely k and Z ′ .Since the wave functions and energies of the states are strongly depen- dent on the magnetic field strength, the configurations are also affected by the magnetic field.With increasing field strength, the ground state undergoes two transitions involving three different electronic configurations [25].In fact, the configuration 2 1 2 s s describes the ground state of lithium atom only for relatively weak fields.
For intermediate field regime, the ground state configuration changes to , , where ( ) ( ) ( ) is the usual hydrogen like function in nl state.
The factor ( ) ij f r is the correlation function given by Equation (3.3).In the above equation, ( ) s χ is the spin function, which is given by [26] [27]: In Equation (3.5),A is the three-particle anti symmetrizer For the strong field, the ground state is constituted by According to the azimuthal symmetry with respect to the magnetic field axis, cylindrical coordinates are suitable for representing the one-particle basis functions for strong field.Hence, in the case of strong field, the following Hylleraas-Gaussian basis set-trial wave function with cylindrical symmetry is a good candidate for highly precise calculations, for the ground as well as the low-lying excited states of lithium [28]: where , i i α β are positive nonlinear variational parameters.The exponents i n ρ and zi n obey the following restrictions: z are suitable for describing excitations.At high field strengths, however, where the magnetic field destroys the spherical symmetry of the problem, an isotropic basis set method would be inefficient.The distribution of the i α is chosen to be peaked around B/4, whereas the i β are well tempered in a large regime [29].

Results and Discussions
In the present paper, we solved the Schrödinger equation for the lithium atom in a magnetic field under the assumption of an infinitely heavy nucleus.Accordingly, we applied the VMC method to calculate the ground state energies of the lithium atom in the magnetic field regime between 0 a.u. and 100 a.u. as well as the lithium like ions up to 10 Z = .All energies are obtained in atomic units ( 1 e m e = = =  ), with a set of 1 × 10 7 Monte Carlo integration points to make the statistical error as low as possible.Our calculations for the ground-state and the two-low lying excited states of the lithium atom in an external magnetic field are based on using three accurate trial wave functions 1 2 , ψ ψ and 3 ψ .In Table 1, we present the best values of the calculated ground state energy and standard deviation of the lithium atom by using the three trial wave functions.The exact value of the ground state energy of Li is taken from [26].
In the presence of a magnetic field an atom and its physics is subject to a variety of changes.For instance, the conserved quantum numbers are reduced to the total angular momentum z-projection M, the total z-parity z Π , the total spin z-projection Z S and the total spin 2 S .For weak fields up to 0.18 γ = , the ground state arises from the field-free  ψ for strong field.
Our work E, σ E [30] E [31] Our work E, σ E [30] E [31] Our work values of the magnetic field strength γ .The previous values obtained in [30] [31] are also presented for comparison.It is seen from Table 2 that our energy values are in excellent agreement with the corresponding previous works for all values of 0 γ > .It is clear that the energy values of the ground state of the li- thium atom in the presence of magnetic field are closer to the values obtained by Wang and Qiao [30] for the ground as well as the two states 2 1 2 s s and 2 1 1 2 s p − .In the absence of the field, our result is in excellent agreement with the exact value.
Also, we calculated the total energies of the lithium like ions up to Z = 10 as functions of the magnetic field.In Table 3, we present the total energies of the lithium like ions by using the trial wave function 3 ψ .In most cases of Table 3, the obtained results are in good agreement with the exact values.The associated standard deviations have very small values, which vary between 10 −5 and 10 −4 .This is due to the large number of VMC points.It is clear from Table 3 that our result for the ground state of the lithium atom (Z = 3) is −7.69085 (a.u.), which is very close to the value of [32].To demonstrate the behavior of the total energies of the lithium like ions against the change in the magnetic field strength, the energies are calculated at a fixed value of 0.12212 γ = by using the trial wave function 2 ψ with the cor- responding standard deviation σ .The obtained results are listed in Table 4. Figure 1 shows the variation of the ground state total energy versus the magnetic field strength γ (a.u.) for the three states  , the state 2 1 2 s s represents the ground state of lithium atom.It is a doubly tightly bound state, i.e., it involves two tightly bound orbitals of 1s character.The 2s electronic charge distribution along the z axis expands slightly in weak magnetic fields.Also, we remark that the total energy decreases because of the spin Zeeman-effect.On the other hand, as the field gets larger, the system tends to shrink, and by raising the angular momentum of the most exterior electron, the atom can increase the electronic separation.In the intermediate field regime, the crossover from configuration 2 1 2 s s to which represents the ground state of the lithium atom for high field strengths.For this configuration, in which all single-electron binding energies increase unlimited for γ → ∞ , a shrinking process of this distribution in z direction is also visible.In general, one can see that for each separate configuration, the effect of the increasing field strength consists in compressing the electronic distribution towards the z axis and the spherical symmetry of the atom is clearly broken.Moreover, we remark that the total energy decreases because of the spin Zeeman-effect.On the other hand, as the field gets larger, the system tends to shrink,   and by raising the angular momentum of the most exterior electron, the atom can increase the electronic separation.Furthermore, we present in Figure 2 the

Conclusions
In the present paper, we used the well-known VMC method, with a set of 1 × 10 7 Monte Carlo integration points to make the statistical error as low as possible, to study the ground state of the lithium atom and its like ions up to Z = 10 in the presence of a magnetic field.Also, we present the results of calculating the energies of the first two excited states of the lithium atom in the external magnetic field.The method is flexible enough to yield precise results for arbitrary field strengths.Our calculations for the ground and first two excited states are performed for magnetic field strengths ranging from zero up to 100 γ = . The effect of increasing the magnetic field on the ground state is studied.With increasing field strength, the ground state undergoes two transitions involving three different electronic configurations.For small values of γ (weak field) the state 2 For strong field, our calculations are based on using a compact few-parameter trial wave function which can be considered as the most accurate among several existing few-parameter trial wave functions for the lithium atom and its ions.
Our results for the total energies exhibit good accuracy comparing with the existing data in the literature.
The energies were plotted as functions of the magnetic field strengths γ to show graphically the effect of the magnetic field on the behavior of the total energy.We conclude that the presented analysis introduced in this paper ensures that the VMC method can be considered as an efficient tool to study the three-electron system under the effect of the strong magnetic field.Furthermore, our results confirm that the chosen trial wave functions are more efficient in the low and high field regimes.
Zeeman term.For weak field, the magnetic field parameter falls in the interval 0 0.18 γ < ≤ a.u.For intermediate field, γ takes values in the interval 0 much of the same features as does the exact wave function.One possible guideline in choosing the trial wave function is the use of the constraints about the behavior of the wave function when the distance between one electron and the nucleus or two electrons approaches zero.These constraints are called "cusp conditions" and are related to the derivative of the wave function.More details about the trial wave function can be found in[20] [24].
configuration.To represent the atomic configuration of all the elec- trons in intermediate field, for ground state as well as for the low-lying excited states, another type of trial wave function is employed to obtain our results.This trial wave function is constructed from the hydrogen-like wave function as follows:

P
stands for the permutation of 123 into ijk.
function given by (3.6) and f is the correlation function given by (3.3).The function (3.7) is an eigenfunction of the z-component of the orbital angular momentum with an eigenvalue of i m and is an eigenfunction of the z-parity with eigenvalue ( ) the ρ dependence of the lowest Landau state in the field B if we choose i α to be B/4 and thus represents an adjustment to the existence of the magnetic field.Also, the monomials 2 i k ρ and zi n

Figure 1 .
Figure1.The total energy versus the magnetic field strength γ (a.u.) for the states2
total energies of the lithium like ions (at 0.12212 γ =) versus Z, up to Z = 10, by using the trial wave function 2 ψ .

SM
the ground state where for intermediate and strong fields the ground state configuration arises from the 2 transition field strengths are determined and precise values for the field strength parameter γ corresponding to the crossovers of the ground state configurations are obtained.For this purpose, suitable types of trial wave functions are used to gain greater flexibility in representing the different atomic configurations when the magnetic field increases, and the cylindrical symmetry associated with the magnetic field becomes dominant.For the weak field, we used a trial wave function which is a linear combination of terms constructed from the single-particle hydrogenic wave functions in such a way that the total spin of = .For the intermediate field, where the ground state configuration changes to 2 1 1 2 s p − configuration, we used for the ground state as well as for the low-lying excited states another type of trial wave function.This trial wave function is constructed from the hydrogen-like wave function as explained in Section 3.

Table 2 ,
we present our results for the behavior of the total energy at different

Table 2 .
Ground and first excited states of the lithium atom (in a.u.) as functions of the magnetic field strength.Previous results

Table 3 .
The energy values E (a.u.) and the magnetic field strengths γ as the ground state configuration changes for lithium-like ions.The standard deviations, σ, of our results are

Table 4 .
Total energies of the lithium like ions versus Z, up to Z = 10, for a fixed value of 0.12212 γ = a.u.by using the trial wave function 2