Ab initio calculations of oscillator strengths and lifetimes of low-lying states in Mo VI

Ab initio calculations of oscillator strengths and lifetimes of low-lying states in Mo VI


Relativistic coupled-cluster (RCC) calculations have been performed to estimate the electromagnetic forbidden transition probabilities, oscillator strengths and lifetimes of many low-lying states of five times ionized molybdenum (Mo VI). Contributions from the Breit interaction up to the first order of perturbation have been examined. Our results are in good agreement with the available other reported theoretical and experimental results. A long lifetime about of the first excited state, 4d , has been predicted which can be a very useful criteria in the doping process of thin films. Correlations trends from various RCC terms to the transition amplitude calculations are discussed.

I Introduction

Electromagnetic forbidden, both magnetic dipole (M1) and electric quadrupole (E2), transitions of Mo VI are important for temperature and density estimations of tokamak plasmas (1); (2), especially in the collision-radiative model (3); (4). Long lifetimes of metastable states are dominated by these forbidden transitions and these states are generally difficult to observe in the laboratory plasmas due to strong collisions. However, these forbidden transitions of Mo VI have been observed in laboratory in electron spin resonance experiment (5) and therefore, they must be one of the sources of density estimations in astrophysical plasmas where collisions are very low due to high dilute interstellar medium (6). Accurate estimation of abundances of molybdenum in the atmosphere of the evolved stars are important to understand the stellar nucleosynthesis (7).

Hexavelent molybdenum, isoelectronic to rubidium with 4p4d as ground state configuration, is generated by electron impact in the atomic collision process. The electron-impact ionization of multiply charged Mo ions, relevant to astrophysics and laboratory plasma research, have also been investigated (8). Recently, Fisker et al. have given the possibility of the origin of the lightest isotope of molybdenum in proton rich type II supernova (9). The necessity of accurate estimation of allowed dipole transition strengths to find out their mixing these effects in the dipole forbidden transitions in Mo VI is explicitly discussed by T. Yamamoto (10). Again, the transition strength between the fine structure states of 4d lavel can reflect the electronic structure of Mo VI in crystal (11).

A few calculations have been carried out to study the on electric dipole (E1) transitions in Mo VI over the last few decades using the mean-field theory (12); (13). More recently, J. Reader (15) has estimated the E1 transition probabilities among low-lying states by estimating transition strengths in the semiempirical approach with the experimental excitation energies.

For this single reference system, Mo VI, we have performed relativistic coupled-cluster (RCC) calculation with single (S), double (D) and partial triple (T) excitations in the framework of Fock space multi-reference (FSMR). Both the excitation energies and transition probabilities are determined using this RCC method using which lifetimes of many low-lying states are estimated.

Ii Theory and Method of Calculations

ii.1 Theory

The oscillator strength for E1 transition from to is given as


where is the excitation energy between the upper and lower states and is the degeneracy factor of the upper state with total angular momentum .

The single particle reduced matrix elements for the E1, E2 and M1 transition operators are given in ((21)). The emission transition probabilities (in sec) for the E1, E2 and M1 channels from states f to i can be expressed as


where is the transition strength for the coressponding operator (in a.u.) and (in Å) is the corresponding transition wavelength.

The lifetime of a particular excited state can be computed by the reciprocal of the total transition probability, (in sec), arising from all possible states due to spontaneous electromagnetic transitions, i.e.


ii.2 Fock Space Multi-reference RCC theory

The FSMRCC method is one of the most powerful highly correlated many-body approaches due to its all order structure to account the correlation effects (16). The FSMRCC, which is mainly meant for multi-reference systems, is used here for the one valence electron and has been described in details elsewhere (16); (17); (18); (19). Here we present the method briefly.

We first consider the Dirac-Coulomb Hamiltonian for a closed-shell electron system which is given by


with all the standard notations often used.

The theory for a single valence system is based on the concept of common vacuum for both the closed shell - and open shell -electron systems, which allows us to formulate a direct method to determine energy differences (electron attachment energy or negative of the ionization potential). Also, the holes and particles are defined with respect to the common vacuum for both the electron systems. Model space of an (n,m) Fock-space contains determinants of holes and particles distributed within a set of orbitals known as active orbitals. For example, in the present article, we are dealing with (0,1) Fock-space which is a complete model space (CMS) by construction and is given by


where ’s are the expansion coefficients of , and ’s are the model space configurations made of DF orbitals. The dynamical electron correlation effects are introduced through the valence-universal wave-operator (16); (17)




At this juncture, it is convenient to single out the core-cluster amplitudes and call them . The rest of the cluster amplitudes will henceforth be called . Since is in normal ordered, we can rewrite Eq. (8) as


The “valence-universal” wave-operator in Eq. (10) is parameterized in such a way that the states generated by its action on the reference space satisfy the Fock-space Bloch equation




Here, P is the projection operator of model space. Eq. (11) is valid for all (k,l) starting from k=l=0 (i.e., the core problem) to some desired parent model space, with k=m, l=n. In this present calculation, we truncate Eq. (9) at and .

In this work, single () and double () excitations are considered for and clusters operator. Therefore, the total correlated wavefunction of the system with single valence orbital v, can be written as


Important triple excitations, correspond to the correlation to the valence orbitals, are included in the open shell FSMRCC-SD calculations by an approximation that is similar in spirit to FSMRCC-SD(T) (20). The approximate valence triple excitation amplitudes are given by


where are the amplitudes corresponding to the simultaneous excitations from core orbitals and valence to virtual orbitals , and , respectively. and are the connected composites involving and , and and , respectively, where is the two electron Coulomb ()integral and ’s are the orbital energies.

The transition matrix element due to any operator can be expressed as

Iii Results and Discussions

In the present calculation, the radial wavefunctions of DF orbitals of closed shell Mo VII are obtained using Gaussian type orbitals (GTO) basis with finite nuclear size as discussed in our earlier paper (22). We have used universal basis set, where the exponent is related with two parameters and , same for all the symmetries, expressed as


We have considered and as 0.00625 and 2.72, respectively, after obtaining best fit of the bound orbital energies and evaluating the expectation values of different radial functions ( , , ) generated with GTOs and GRASP2 (23). In the DF calculations, we have taken 22, 20, 17, 15, and 12 number of GTOs for s, p, d, f, and g type symmetries, respectively, to generate the atomic orbitals. In Fig. 1, we have given the relative errors obtained for different orbitals in the calculations of these quantities using the above chosen parameters. Since these errors are very small, it shows that there is a good agreement between results. We assume that both the bound and continuum orbitals generated using the above parameters will describe well both inside and outside of the nucleus. Therefore, we have considered all the orbitals obtained using GTOs for the rest of the calculations.

The number of the DF orbitals for different symmetries used in the present calculation is based on convergent criteria of core correlation energy of Mo VII for which it satisfies numerical completeness. The number of DF orbitals considered for s, p, d, f, and g type symmetries in the RCC calculations are 12, 11, 10, 9 and 8, respectively; and among them 9, 8, 7, 5, and 5 are bound orbitals, respectively, including all the core orbitals. The amplitudes are first determined by solving the closed shell RCC equations for the closed-shell system (Mo VII), then amplitudes are solved from the open-shell equations for the single-valence states of Mo VI.

Figure 1: The relative energies and expectation values of , and of the DF GTO orbitals to the DF GRASP orbitals.

Table I summarizes the calculated excitation energies (EE) and fine structure splitting (FS) of low-lying excited states and their comparison with the recent experimental results (15). The average deviation is around 0.5% for EE. We have also presented the contribution from the partial triple excitations to the EE (E), which is around 0.3% to the total EE. We have examined the first order excitation energy corrections due to Breit interaction using large scale relativistic CI calculations. Maximum contribution is coming for 5s state, which is around , whereas, contributions to 5p and 4f states are around and ; respectively. For 5d state, it is as small as 54 consistent with the result obtained by Pan and Beck (14). Since, the contributions due to Breit interaction are relatively small, we do not consider them here selfconsistently to evaluate wavefunctions.

State CC Exp. CC Exp.
4d 0 0 0
4d 2670.55 2583.50 2670.55 2583.50 572.36
5s 118536.63 119725.62 396.72
5p 181795.01 182404.47 220.12
5p 186825.31 187331.19 5030.30 4926.72 235.92
4f 269411.28 267047.22 492.50
4f 269653.95 267456.84 242.67 409.62 621.98
5d 282598.52 282825.59 697.69
5d 283396.42 283610.94 797.90 785.35 816.42
6s 315632.41 313806.81 857.47
6p 342531.49 340570.78 690.67
6p 344853.69 342562.44 2322.21 1991.66 697.91

Ref. (15).

Table 1: Excitation energies (EE) and fine structure splitting (FS) (in cm) of Mo VI and their comparison with experimental results. Contribution from absolute value of partial triple excitation to the EE (E) are also presented.

Since the transition rate is proportional to the square of the transition amplitude, therefore precise description of the wavefunction is necessary due to one order higher dependance on wavefunctions than energy. In Table II, we compare the E1 transition amplitude in both length and velocity gauges for few transitions. We find a good agreement between them, which is one of the characteristics to judge the accuracy of the wavefunctions.

Upper Lower D D
5p 4d 0.9851 0.8341
5p 5s 1.7554 1.7522
5p 4d 0.4288 0.3619
5p 5s 2.4885 2.5678
4f 4d 1.3928 1.5438
4f 4d 0.4272 0.4474
4f 4d 2.6555 2.5330
5d 5p 2.7486 2.6942
5d 5p 1.2593 1.2278
5d 5p 3.9600 3.7900
6s 5p 0.9771 0.8873
6s 5p 1.4711 1.1963
6p 4d 0.1984 0.1686
6p 5d 2.5559 2.2828
6p 6s 3.4776 3.1844
6p 4d 0.2845 0.2605
6p 5d 1.1055 0.9949
6p 6s 4.9149 4.6861
Table 2: Absolute values of E1 transition amplitude in length (D) and velocity (D) gauges for Mo VI.

In Table III, we compare our ab initio oscillator strength values correspond to E1 transitions with the recent semi-empirical calculations by Pan and Beck (14) and by Reader (15). Reader has obtained wavefunctions using fitting parameters by comparing calculated and experimental energies, whereas, Pan and Beck have used relativistic CI method for the available transitions. Our calculated values of oscillator strength for the 4d 4f transition at the DF level given in the table, agrees well with similar calculations by Zilitis (13), 1.023.

Term Present calculations Other calculations
Upper Lower DF CC (15) (14)
4d 4f 1.0099 0.3967 0.3226 0.2896
4d 4f 0.0484 0.0246 0.0153 0.0139
Table 3: Oscillator strengths for E1 transitions in length form and their comparison with earlier results.

Weighted oscillator strengths correspond to E1 transitions are presented in Table IV. Here we have used the length gauge values of E1 transition amplitudes and our calculated wavelengths. All these transitions, fall in ultraviolet and visible regions, are useful for astrophysical observations and may be for laboratory researches. According to Cowan (24), if the initial states are dominated by and final states are dominated by , the oscillator strength ratio of , and transitions in a given multiplet are 6:5:1, which we find the same for 4d 5p transitions.

Upper Lower gf
5p 4d 54.82 0.5354
5P 5s 159.54 0.5920
5p 4d 53.38 0.1042
5p 4d 54.12 0.9676
5p 5s 147.91 1.2844
4f 4d 37.44 1.5868
4f 4d 37.81 0.1478
4f 4d 37.75 5.7164
5d 5p 99.58 2.3133
5d 5p 104.71 0.4613
5d 4f 633.77 0.3446
5d 5p 103.86 4.5998
5d 4f 603.72 0.0260
5d 4f 619.03 0.5054
6s 5p 76.10 0.3881
6s 5p 79.06 0.8468
6p 4d 29.36 0.0409
6p 5s 45.28 0.0415
6p 5d 173.17 1.1894
6p 6s 373.63 0.9882
6p 4d 29.19 0.0092
6p 4d 29.41 0.0841
6p 5s 44.87 0.0559
6p 5d 167.40 0.2311
6p 5d 169.63 2.0902
6p 6s 347.75 2.1442
Table 4: Transition wavelengths (in nm) and weighted oscillator strengths (gf) corresponding to electric dipole (E1) transitions of Mo VI.

In Table V, we present M1 and E2 transition probabilities and their corresponding wavelengths. However most of the transitions come in ultraviolet region, there are few transitions fall in infrared region. Though these transitions produce weak lines but they are important parameters in astrophysical studies. As expected, transition probability for E2 transitions come greater in value than that of M1 transitions except for the transitions fall in infrared region.

Upper Lower A A
4d 4d 3870.71 3.937210 2.00610
5s 4d 83.52 6.991910 2.867710
5s 4d 85.36 9.344110
5p 5p 2029.74 4.240210 1.1401
4f 5p 118.14 4.538810
4f 5p 125.44 9.858010 1.619110
4f 5p 124.80 4.225810
4f 4f 24412.87 7.555510 1.319310
5d 4d 35.35 1.122110 6.479510
5d 4d 35.68 4.843510 2.838210
5d 5s 61.31 1.209010 5.762710
5d 4d 35.25 3.111110 2.5858
5d 4d 35.58 1.274610 1.603410
5d 5s 61.01 1.211210
5d 5d 12733.17 2.979010 5.386010
6s 4d 31.86 1.729710 1.597010
6s 4d 32.13 8.983610
6s 5s 51.52 2.9841
6s 5d 322.77 2.645610 2.671710
6s 5d 331.17 3.563810
6p 5p 63.22 8.445410
6p 5p 65.25 7.491110 4.007310
6p 4f 136.01 2.274910
6p 5p 62.44 3.608910 1.377610
6p 5p 64.42 3.580710 1.6161
6p 4f 132.42 3.535710 3.152610
6p 4f 131.34 2.050510
6p 6p 5020.93 9.961910 1.118510
Table 5: Transition wavelengths (in nm) and transition probabilities corresponding to electric quadrupole (E2) and magnetic dipole transitions (M1) (in sec) of Mo VI.

In Table VI, the contributions due to the different correlation terms like core-correlation, pair-correlation, core-polarization and important two-body contributions are estimated for few transitions to highlight the effect of correlations. Significant correlation contributions from the higher order core-polarization, like , are noticeable compared to the lowest order contributions. Also, contributions from two-body correlations are almost comparable for most of the cases. For the M1 transition, 4d 4d , there are cancelation observed among different correlation effects. The contribution comes from core-polarization term is more compared to other term. Even in the case of 5s 4d M1 transition, the Dirac-Fock contribution is almost canceled by the lowest order of core-polarization, which makes the core-corelation effect more dominant to the total value of transition matrix element.

Dirac-Fock Core-corr. Pair-corr. Core-polar. Core-polar. Two-body Total
(lowest) (higher) contr.
5p(1/2) 4d(3/2) (E1) -1.0898 -2.1536E-3 2.4778E-2 8.1783E-2 -8.4547E-3 -8.5320E-4 -0.9851
5p(3/2) 4d(3/2) (E1) -0.4704 -8.5909E-4 1.1185E-2 3.1035E-2 -3.4575E-3 -3.6146E-4 -0.4288
4d(5/2) 4d(3/2) (E2) 1.4451 -2.5900E-2 -3.1118E-2 -1.3272E-1 3.6817E-3 -1.1045E-3 1.2432
5s(1/2) 4d(3/2) (E2) 2.4193 2.8071E-3 -6.8921E-2 -3.1504E-2 2.1471E-2 -2.6223E-4 2.3176
4d(5/2) 4d(3/2) (M1) 1.5488 -1.3535E-2 -1.4813E-4 1.5097E-4 1.2398E-2 -1.8252E-3 1.5284
5s(1/2) 4d(3/2) (M1) -2.5086E-5 9.4303E-6 -3.1876E-7 2.4685E-5 -4.8217E-6 6.6800E-9 3.5805E-6
Table 6: Explicit contributions from the CCSD(T) calculations to the absolute magnitude of transition amplitudes.

Table VII summarizes the calculated lifetime of the low-lying excited states. Recent calculations of the lifetimes of 5p and 5p states by Zilitis (13) are also compared here. Here we can see that the lifetime of the 4d state comes in the order of second, which suggest that Mo VI can be used for uniform doping in thin film. The lifetime of the 5s state is found to be of the order of microsecond due to only forbidden transition contributions. Lifetime of the states are larger than the states due to strong contributions from the allowed dipole transitions 4f 4d .

Term Present calculations Other calculations
4d 4.9854
5s 6.121310
5p 1.496810 1.30010
5p 1.415110 1.26010
4f 7.162810
4f 2.946210
5d 2.157510
5d 2.092210

Ref. (13).

Table 7: Radiative lifetimes (in sec.) for different low-lying states of Mo VI.

Iv Conclusion

Forbidden transition probabilities among the low-lying states of Mo VI relevant for astro- and plasma physics are calculated using highly correlated relativistic coupled-cluster method for the first time in literature to the best of our knowledge. The lifetime of the 4d state is found to be around 5 second, which will be useful in many physical processes. Contributions of different correlation terms are discussed and found strong effect from higher order core-polarization. In the near future, present work will motivate experimentalists to verify our results due to its importance in many areas in physics.

V Acknowledgment

We are grateful to Prof B P Das and Dr Rajat K Chaudhuri , Indian Institute of Astrophysics, Bangalore for providing the CC code. One of us (Narendra) would like to recognize the support of Council of Scientific and Industrial Research (CSIR), India.


  1. Feldman U 1981 Phys. Scr. 24 681.
  2. Herter T, Briotta D A, Gull G E, Shure M A and Houck J R 1982 Astrophys. J. 262 164.
  3. Fournier K B, Pacella D, May M J, Finkenthal M and Goldstein W H 1997 Nucl. Fusion 37 825.
  4. Quinet P 2002 J. Phys. B 35 19.
  5. Baur J, Kunzer M, Dombrowski K F, Kaufmann U, Schneider J, Baranov P G, and Mokhov E N 1997 Inst. Phys. Conf. Ser. 155 933.
  6. Charro E and Martin I 2003 Astrophysics J. 585 1191.
  7. Orlov M Ya and Shavrina A V 1988 Astrophysics28 354.
  8. Hathiramani D, Hartenfeller U, Aichele K, Hofmann G, Steidl M, Stenke M, Voelpel R and Salzborn 1998 J. Phys. B 31 2101.
  9. Fisker J L, Hoffman R D and Pruet J 2009 Astrophys. J. 690 L135.
  10. Yamamoto T 2008 X-Ray Spectrom. 37 572.
  11. Szotek Z, Temmerman W M, Svane A, Petit L, and Winter H 2003 Phys. Rev. B 68 104411.
  12. Migdalek J and Baylis W E 1979 J. Quant. Spectrosc. Radiat. Transfer 22 127.
  13. Zilitis V A 2007 Optics and Spectroscopy 103 895.
  14. Pan L and Beck D R 2006 Phys. Scr. 73 607.
  15. Reader J 2010 J. Phys. B 43 074024.
  16. Lindgren I and Morrison J 1985 Atomic Many-body Theory 3 ed. G. E. Lambropoulos and H. Walther (Berlin: Springer).
  17. Lindgren I and Mukherjee D 1987 Phys. Rep. 151 93.
  18. Haque A and Mukherjee D 1984 J. Chem. Phys. 80 5058.
  19. Pal S, Rittby M, Bartlett R J, Sinha D and Mukherjee D 1987 Chem. Phys. Lett. 137 273; 1988 J. Chem. Phys. 88 4357.
  20. Raghavachari K, Trucks G W, Pople J A and Head-Gordon M 1989 Chem. Phys. Lett. 157 479; Urban M, Noga J, Cole S J and Bartlett R J 1985 Chem. Phys. Lett. 83 4041.
  21. Johnson W R, Plante D R and Sapirstein J 1995 Adv. At. Mol. Opt. Phys. 35 255.
  22. Dixit G, Nataraj H S, Sahoo B K, Chaudhuri R K and Majumder S, 2008 J. Phys. B 41 025001; 2008 Phys. Rev. A 77 012718.
  23. Parpia F A, Fischer C F and Grant I P 2006 Comp. Phys. Comm. 175 745.
  24. Cowan R D 1981 The Theory of Atomic Structure and Spectra (Berkeley, CA: University of California Press).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description