Ab initio theory study of laser cooling of barium monohalides

Abstract

The feasibility of laser cooling barium monohalides BaX (X = F, Cl, Br, I) is investigated using ab initio methods with the inclusion of spin–orbit coupling (SOC) effects. Calculated spectroscopic constants for BaF, BaCl, BaBr and BaI are in very good agreement with the available experimental measurements. The results demonstrate that the calculated electronic structure is accurate and can be used to establish the optical scheme of laser cooling. The highly diagonal Franck–Condon factors (FCFs) (BaF: f₀₀ = 0.980, f₁₁ = 0.939, f₂₂ = 0.894; BaCl: f₀₀ = 0.998, f₁₁ = 0.995, f₂₂ = 0.992) between the X²Σ⁺_1/2 and A²Π_1/2 states are determined, which are found to be in good agreement with previous theoretical results. The radiative lifetimes (BaF: 39.13–41.20 ns; BaCl: 117.99–110.23 ns) of the A²Π_1/2–X²Σ⁺_1/2 transition for the first five vibrational levels show that the A²Π_1/2 is a rather short lifetime state. The current study indicates that BaF and BaCl are two good choices of molecules for laser cooling. Therefore, BaI and BaBr are not promising laser-cooling candidates because the FCFs of the A²Π_1/2–X²Σ⁺_1/2 transition are off-diagonal. We further propose the three-laser cooling schemes based on the A²Π_1/2–X²Σ⁺_1/2 transition for BaF and BaCl.

---

1. Introduction

The transverse laser cooling of SrF^1,2 and YO³ has overthrown the notion that laser-cooling molecules is impracticable. The criteria to be considered candidates for laser cooling are a strong transition with highly-diagonal Franck–Condon factors (FCFs) in a wavelength region accessible to tabletop lasers and strong optical forces resulting from the short radiative lifetime. Diagonal FCFs minimize the number of lasers required to keep the molecule in a closed-loop cooling cycle. Naturally, the experimental demonstrations have initiated the search for other molecules which are appropriate for laser cooling. In 2004, Di Rosa⁴ conducted a brief survey of promising candidates for laser cooling. SrF and YO did not occur on his study and there clearly could be others that have not been investigated by theoretical and experimental studies. In this paper, thus, we focus on the theoretical studies of the laser cooling of barium monohalides, though neither have yet been cooled.

In recent years, the interest in laser cooling of diatomic polar molecules has grown owing to a variety of prospective applications.⁵ Laser cooling of SrF,^1,2 KRb,⁶ YO³ and CaF⁷ has been experimentally performed, and polar molecules (such as RaF,⁸ AlH,⁹ AlF⁹) have been investigated in previous theoretical research. Alkaline-earth metal atoms and molecules are of persistent interest in the investigation of cold atoms and molecules. These molecules like monohydrides MH (M = Be, Mg, Ca, Sr, and Ba)¹⁰ and some monohalides (such as BeF,¹¹ MgF¹²) have been identified as the promising candidates for laser cooling theoretically. However, the study on laser cooling of barium monohalides BaX (X = F, Cl, Br, I) is very limited. In addition to being nominated as laser cooling candidates, BaX (X = F, Cl, Br, I) play an important role in stellar spectroscopy.¹³ Considering that the A²Π → X²Σ⁺ transition of barium monohalides has been suggested by Di Rosa⁴ and experiments have been made for SrF,^1,2 YO³ and CaF,⁷ we believe it is possible to cool barium monohalides with a similar experimental method.

In the current study, we focus on investigating the feasibility of direct laser cooling of BaX (X = F, Cl, Br, I), then proposing the scheme to perform laser cooling of BaX (X = F, Cl, Br, I). Karthikeyan et al.¹⁴ have evaluated the FCFs for electronic transitions A–X of BaF molecule by employing a reliable numerical integration procedure. The calculated FCF (f₀₀) of BaF is 0.951, and the laser excitation wavelength is about 857 nm. For BaF, Kang et al.¹⁵ also reported results about FCFs by ab initio. The radiative lifetime obtained for the A²Π state of the BaF molecule is 56.0 ± 0.9 ns using laser spectroscopy.¹⁶ The first seven electronic states of the BaCl molecule have been studied by laser induced fluorescence, combined with high resolution Fourier spectroscopy.¹⁷ The zero pressure radiative lifetime of BaCl has been measured by laser-induced fluorescence in the near-infrared wavelength region.¹⁸ Husain et al.¹⁹ have determined the spectroscopic parameters of the X-state and low lying excited states of BaBr, and calculated the FCFs using a simplified Morse oscillator model following Tuckett²⁰ (A–X: f₀₀ = 0.6903). The A²Π state of BaBr has been analyzed by Ernst et al.²¹ Miliordos et al.²² have been examined the electronic structure of BaI by ab initio multi-reference configuration interaction (MRCI) and coupled cluster (RCCSD(T)) methods. As mentioned above, BaX (X = F, Cl, Br, I) has been investigated extensively in the past. However, there is a lack of systematic studies of laser cooling of BaX (X = F, Cl, Br, I) up to now.

The paper is organized as follows. The theoretical methods and the basis sets used in the ab initio calculations for the included electronic states of BaX (X = F, Cl, Br, I) are briefly described in Section 2. We present the results and data discussions, and propose the laser cooling schemes for BaX (X = F, Cl, Br, I) in Section 3. In Section 4, the conclusions of this work are given.

2. Computational details

To investigate the feasibility of direct laser cooling of barium monohalides, we first calculate the electronic structure for the two lowest doublet electronic states (A²Π and X²Σ⁺), then determine the radiative properties including radiative lifetimes, FCFs and diode laser excitation wavelengths of the A²Π → X²Σ⁺ transitions. All the present ab initio calculations are carried out in the C_2ν point group using the MOLPRO²³ package. The spectroscopic constants, including equilibrium bond distance (R_e), harmonic frequency (w_e), anharmonic vibrational frequency (w_eχ_e), rotational constant (B_e) and electronic transition energy (T_e) for the A²Π and X²Σ⁺ states of BaX (X = F, Cl, Br, I) are evaluated by solving the nuclear Schrödinger equation by Le Roy's LEVEL 8.0 program.²⁴ With the potential energy curves (PECs) and transition dipole moments of different electronic states, this program can also obtain the FCFs and radiative lifetimes of the various vibrational levels. The FC factor is one of the parameters which controls the intensity distribution in the emission of molecular bands. The square of the overlap integral is termed as FC factor (f_ν′ν),²⁵where ψ_ν′ and ψ_ν are the vibrational wave functions for the upper and lower states, respectively. The radiative lifetime τ_ν′ is calculated by the following formula,where A_ν′ν is the Einstein A coefficient for the rate of spontaneous emission from initial rovibrational level (ν′, J′) into final rovibrational level (ν, J),^26,27in the above formula, S(J′,J) is the Hönl–London rotational intensity factor, M(r) is the transition dipole function, ψ_ν′,J′ and ψ_ν,J are normalized radial wave functions.

The complete active space self-consistent field (CASSCF)^28,29 and MRCI plus Davidson correction (MRCI + Q)^30–32 methods have been performed for the ab initio electronic structure calculations of BaX (X = F, Cl, Br, I). Scalar relativistic effects are accounted for using the Douglas–Kroll–Hess (DHK)^33,34 transformation of the relativistic Hamiltonian. For the heavier Ba atom, the spin–orbit coupling (SOC) effects are important for the electronic structure. In this work, SOC effects are evaluated with one- and two-electrons Breit–Pauli operators³⁵ at the MRCI level of theory. It is worth noting, we can neglect the rotational degrees of freedom in the first approximation, since we focus on investigating the vibrational cooling.

To improve the accuracy of potential energy curve calculations, the effects of the core–valence correction and the relativistic correlation are taken into account. This was done by adopting the small-core scalar relativistic effective core potential (ECP) ECP10MWB,³⁶ ECP10MDF,³⁷ ECP28MDF,³⁸ and ECP46MDF³⁹ together with the corresponding valence basis sets for Cl, Br, I and Ba, respectively, and, the AVQZ⁴⁰ all-electron basis sets for the F atom. Within the C_2ν point group symmetry, all molecular orbitals were labeled by their symmetry (a₁, b₁, b₂, a₂). In our calculations of BaF molecule, eight molecular orbitals are selected as active spaces, which correspond to the 2s2p shells of the F atom and 5p6s shells of the Ba atom. That is to say, ten electrons are distributed in a (4, 2, 2, 0) active space for BaF molecule. For BaCl, BaBr and BaI, these active spaces are also (4, 2, 2, 0). The 3s3p shells of the Cl atom, 4s4p shells of the Br atom and 5s5p shells of the Br atom are put in the active space, respectively.

3. Results and discussions

3.1 Potential energy curves (PECs) and spectroscopic constants

Because the A²Π → X²Σ⁺ transition has been suggested by experiments for SrF, YO and CaF, we only carry out ab initio calculations to the two states (A²Π and X²Σ⁺). The PECs of A²Π and X²Σ⁺ states of BaX (X = F, Cl, Br, I) are presented in Fig. 1(a) and (b). Fig. 1 shows that the PECs of A²Π and X²Σ⁺ states look surprisingly similar. It implies that highly diagonal FCFs for the A²Π → X²Σ⁺ transition are possible.

Fig. 1 (Color online) Potential energy curves of X²Σ⁺ (a) and A²Π (b) of BaX (X = F, Cl, Br, I) at MRCI + Q/ECP46MDF(Ba), AVQZ(F), ECP10MWB(Cl), ECP10MDF(Br), ECP28MDF(I).

For the heavier Ba atom, SOC effects are important for the electronic structure. Concerning the effects of SOC on the electronic structure of BaX (X = F, Cl, Br, I), the A²Π state splits into A²Π_1/2 and A²Π_3/2. Our present calculated spectroscopic constants of A²Π_1/2 and X²Σ⁺_1/2 with available experimental values are tabulated in Table 1. The reasons why we focus on A²Π_1/2 and X²Σ⁺_1/2 are explained as follows. Many experimental studies provided some accurate and important spectroscopic data for A²Π_1/2 and X²Σ⁺_1/2 states. Earlier theoretical researches^14,15 reported BaF has highly diagonal FCFs between A²Π_1/2 and X²Σ⁺_1/2.

Table 1 Spectroscopic constants of X²Σ⁺_1/2 and A²Π_1/2 states of BaX (X = F, Cl, Br, I) calculated at MRCI + Q/ECP46MDF(Ba), AVQZ(F), ECP10MWB(Cl), ECP10MDF(Br), ECP28MDF(I). T_e is the adiabatic relative electronic energy referred to the ground state, R_e is the equilibrium internuclear distance, w_e and w_eχ_e are the harmonic and anharmonic vibrational frequency, B_e is the rotational constant

Molecule	States	Te (cm^−1)	Re (Å)	we (cm^−1)	weχe (cm^−1)	Be (cm^−1)	Ref.
BaF	X^2Σ^+1/2	0	2.186	460.72	1.79	0.2114	This work
		0	2.161^b	469.41^a	1.84^a	—	Expt. (^aRef. 41; ^bRef. 42)
	A^2Π1/2	11 893	2.199	436.58	1.72	0.2089	This work
		11647^a	2.183^c	435.50^a	1.68^a	—	Expt. (^aRef. 41; ^cRef. 43)
BaCl	X^2Σ^+1/2	0	2.728	267.38	0.79	0.0814	This work
		0	—	279.89	0.81	0.0840	Expt. (Ref. 17)
	A^2Π1/2	10 657	2.683	261.33	0.80	0.0810	This work
		10 679	—	257.30	0.81	0.0811	Expt. (Ref. 17)
BaBr	X^2Σ^+1/2	0	2.927	175.14	0.39	0.0390	This work
		0	2.847	192.30	0.41	—	Expt. (Ref. 19)
	A^2Π1/2	10 171	2.911	169.98	0.40	—	This work
		—	2.899	177.13	0.40	—	Expt. (Ref. 19)
BaI	X^2Σ^+1/2	0	3.126	149.92	0.27	0.0262	This work
		0	3.088	152.16	0.27	—	Expt. (Ref. 44–46)
	A^2Π1/2	9753	3.167	140.47	0.27	0.0255	This work
		9605	3.139	141.75	0.28	—	Expt. (Ref. 45 and 46)

---

For the X²Σ⁺_1/2 state of BaF, the corresponding percentage errors in R_e, w_e and w_eχ_e are 1.16%, 1.85% and 2.72% with respect to experiment,^41,42 respectively. It is also encouraging to see that our calculated T_e, R_e, w_e and w_eχ_e of the A²Π_1/2 state are in good agreement with experiment^41,43 and the relative errors are only 2.11%, 0.73%, 0.25% and 2.38%, respectively. In the case of the X²Σ⁺_1/2 state of BaCl, our present results compare very well with the experimental data¹⁷ (w_e ∼12.51 cm⁻¹, w_eχ_e ∼0.02 cm⁻¹, B_e ∼0.0026 cm⁻¹). And for the A²Π_1/2 of BaCl, our spectroscopic constants calculated at MRCI level are also in reasonable agreement with experiment¹⁷ (T_e ∼22 cm⁻¹, w_e ∼4.03 cm⁻¹, w_eχ_e ∼0.01 cm⁻¹, B_e ∼0.0001 cm⁻¹). Concerning the X²Σ⁺_1/2 of BaBr, the difference of R_e, w_e and w_eχ_e values in the experiment¹⁹ and present work is not signficant, and the relative errors are 2.81%, 8.92% and 4.88%, respectively. While for the A²Π_1/2 of BaBr, Table 1 shows that our present results accord with experiment¹⁹ (R_e ∼0.012 Å, w_e ∼7.15 cm⁻¹, w_eχ_e ∼0.00 cm⁻¹). For BaI, present values of R_e, w_e and w_eχ_e for the X²Σ⁺_1/2 state are 3.126 Å, 140.47 cm⁻¹ and 0.27 cm⁻¹, which are very close to the experimental data.^44–46 Miliordos et al.²² provided a even larger R_e (3.251 Å) and a even smaller w_eχ_e (0.24 cm⁻¹) compared to the experimental data^44–46 by MRCI method. As for the A²Π_1/2 of BaI, the calculated R_e and T_e results are only 0.028 Å and 148 cm⁻¹ larger, while the w_e and w_eχ_e results are only 1.28 cm⁻¹ and 0.01 cm⁻¹ smaller than the observed experimental data.^45,46 The differences between the previous theoretical value²¹ and experimental data^45,46 are bigger (T_e ∼1536 cm⁻¹, R_e ∼0.194 Å, w_e ∼12.75 cm⁻¹, w_eχ_e ∼0.31 cm⁻¹). In summary, the results demonstrate that the electronic structure of BaX (X = F, Cl, Br, I) calculated at MRCI + Q/ECP46MDF(Ba), AVQZ(F), ECP10MWB(Cl), ECP10MDF(Br), ECP28MDF(I) level is accurate and can be used to establish the optical scheme of laser cooling.

3.2 FCFs

The first criterion to be nominated as candidates for laser cooling is highly diagonal FCFs. As listed in the literature, most FCFs f₀₀ of experimental cooling molecules are greater than 0.9. For example, the FCF f₀₀ of SrF reaches 0.98. A larger FCF has significant benefits for limiting the number of lasers required to keep the molecule in a closed-loop cooling cycle. Of course, some cooling schemes for molecules with smaller FCFs have also been proposed.⁴⁷ The FCF values for the A–X transitions of BaX are tabulated in Table 2 and 3. Our computed FCF values of BaF are in good agreement with the recent results.^14,15 For BaBr, the present results at the MRCI level are a little bigger than the values reported by Husain et al.¹⁹ Husain et al.¹⁹ calculated the FCFs using a simplified Morse oscillator model following Tuckett. From Table 2, it could be clearly seen that the Δν = 0 bands of BaX have the strongest transition probability. That is to say, the highly diagonal FCFs f₀₀ (BaF: 0.980 and BaCl: 0.998) for the A–X transition are determined. To visually demonstrate the overlap of the vibrational wave functions for A–X transition, we have plotted the FCFs of spontaneous radiative transitions in Fig. 2. Fig. 2 also shows that BaF and BaCl have large diagonal FCFs.

Table 2 The calculated FCFs f_ν′ν of BaF and BaCl for A²Π_1/2(ν′) → X²Σ⁺_1/2(ν) transitions at MRCI + Q/ECP46MDF(Ba), AVQZ(F), ECP10MWB(Cl) (numbers in bold – published FCF values). Numbers in parentheses indicate the power of 10

Molecules		ν′ = 0	1	2	3
BaF	ν = 0	0.9802	0.0197	0.7824(−4)	0.9904(−6)
	Ref. 14 (^aRef. 15)	0.9510, 0.9810^a	0.0490, 0.0189^a	0.0000, 0.3960(−3)^a	0.0000^a
	1	0.0194	0.9389	0.0414	0.0030(−1)
	Ref. 14 (^aRef. 15)	0.0480	0.8540, 0.9400^a	0.0960	0.0030
	2	0.0041(−1)	0.0401	0.8937	0.0650
	Ref. 14 (^aRef. 15)	0.0020	0.0930	0.7580, 0.8960^a	0.1410
	3	0.3362(−5)	0.0013	0.0620	0.8452
	Ref. 14	0.0000	0.0050	0.1350	0.6660
BaCl	ν = 0	0.9982	0.0017	0.4883(−4)	0.1043(−6)
	1	0.0023	0.9949	0.0032	0.1440(−3)
	2	0.8568(−4)	0.0031	0.9920	0.0045
	3	0.2042(−5)	0.2423(−3)	0.0042	0.9894

---

Table 3 The calculated FCFs f_ν′ν of BaI and BaBr for A²Π_1/2(ν′) → X²Σ⁺_1/2(ν) transitions at MRCI + Q/ECP46MDF(Ba), ECP10MDF(Br), ECP28MDF(I). Numbers in bold – published FCF values (^aRef. 19)

Molecules	f00 f01 f02 f03
	f10 f11 f12 f13
BaI	0.779 0.197 0.023 0.002
	0.192 0.436 0.307 0.059
BaBr	0.791(0.690)^a 0.189(0.216)^a 0.020 0.001
	0.188(0.260)^a 0.440(0.272)^a 0.302 0.061

---

Fig. 2 (Color online) The calculated FCFs of (a) BaF and (b) BaCl for the lowest vibrational levels of the cooling transition A²Π_1/2 → X²Σ⁺_1/2 at MRCI + Q/ECP46MDF(Ba), AVQZ(F), ECP10MWB(Cl).

However for the A–X transition of BaI, the diagonal term f₁₁(0.436) is calculated, and the off-diagonal term f₁₀(0.192), f₁₂(0.307), f₁₃(0.059) are also obtained (see Table 3). These results show that the A–X transition of BaI does not have highly diagonal FCFs. The calculated f₀₀ of 0.779 for BaI is not large enough to be potentially viable for direct laser cooling. The relative probabilities from A²Π_1/2(ν′ = 0) to X²Σ⁺_1/2(ν) are governed by the FCFs which are about 78% for ν = 0, 20% for ν = 1, 2.3% for ν = 2, 0.2% for ν = 3, 0.007% for ν = 4, and negligibly small for all ν > 4. Five cooling lasers are required for limiting the inefficiency or loss of BaI molecule in the cooling cycle. Cycling transitions requiring one or two repump lasers are common in atomic systems, experimentally. For the A–X transition of BaI, the number of lasers required is not practical. The case of BaBr is similar to the case of BaI (see Table 3). So BaI and BaBr are not good candidate molecules for direct laser cooling. Thus we do not need to discuss BaI and BaBr molecules in the following section.

3.3 Spontaneous radiative lifetime and radiative width

Aside from the large diagonal FCFs, sufficiently short lifetime is desirable for rapid laser cooling. A shorter lifetime of the transition could produce a strong Doppler force, which will produce a significant rate of optical cycling. To produce larger spontaneous scattering forces, the rate of optical cycling (10⁵–10⁸ s⁻¹) is suitable. From eqn (1), the radiative lifetime of vibrational level ν′ for a given state is given by the inverse of the total Einstein coefficient. The radiative width (Γ) is related to the lifetime of the state by Γ = 1/(2πcτ). The computed radiative lifetimes and radiative width of BaF and BaCl are listed in Table 4 for the first five vibrational levels (ν′ = 0–5). As shown in Table 4, radiative lifetimes increase with increasing ν′, and radiative widths decrease with increasing ν′. We found that our radiative lifetimes for ν′ = 0 agree well with the experimental results (BaF: (56.0 ± 0.9) ns;¹⁶ BaCl: (107 ± 4) ns¹⁸). The radiative lifetimes (BaF: 39.13–41.20 ns; BaCl: 117.99–110.23 ns) of the A–X transition for the first five vibrational levels show that the A²Π_1/2 is a rather short lifetime state. Thus, BaF and BaCl meet the second criterion for direct laser cooling.

Table 4 The calculated radiative lifetimes (ns) and radiative width (cm⁻¹) (in bold) for the A²Π_1/2 state at MRCI + Q/ECP46MDF(Ba), AVQZ(F), ECP10MWB(Cl)

Molecule	State	ν′ = 0	1	2	3	4
BaF	A^2Π1/2	49.13	49.63	50.15	50.66	51.20
		1.08(−4)	1.07(−4)	1.06(−4)	1.05(−4)	1.04(−4)
BaCl	A^2Π1/2	117.99	115.95	113.99	112.09	110.23
		4.49(−5)	4.57(−5)	4.64(−5)	4.73(−5)	4.81(−5)

---

3.4 Laser cooling scheme

As shown in Fig. 3 (a), we have suggested a three-laser cooling scheme for the BaF molecule. The main cycling laser can drive the X²Σ⁺_1/2(ν = 0) → A²Π_1/2(ν′ = 0) transition at wavelength λ₀₀ = 841.8 nm. Our computed wavelength λ₀₀ agrees well with the result (822.0 nm) given by Kang et al.¹⁵ However, the decays to the X²Σ⁺_1/2(ν = 1) state, to the X²Σ⁺_1/2(ν = 2) state, and the X²Σ⁺_1/2(ν ≥ 3) state from the A²Π_1/2 state occur with probability f₀₁ (= 2%), f₀₂ (≈ 0.008%) and f⁺₀₃(<10⁻⁶), respectively. To enhance the cooling effect, using the X²Σ⁺_1/2(ν = 1) → A²Π_1/2 (ν′ = 0) transition as the first vibrational pump and the X²Σ⁺_1/2(ν = 2) → A²Π_1/2(ν′ = 1) transition for the second pump are required. With this scheme, a BaF molecule can scatter about N_scat = 1/f⁺₀₃ > 10⁶ photons before decaying into the higher vibrational levels (ν ≥ 3) in X²Σ⁺_1/2. Therefore, two repumping lasers with λ₁₀ = 812.2 nm and λ₂₁ = 814.0 nm are needed. The required cooling wavelengths of λ₀₀ = 841.8 nm, λ₁₀ = 812.2 nm and λ₂₁ = 814.0 nm are located in the near-infrared range, where the continuous wavelength laser radiation can be easily generated. One problem with the above laser cooling scheme is the presence of a metastable state A′²Δ between the X²Σ⁺_1/2 and A²Π_1/2.⁴¹ It is notable that the decay rate γ ∝ d²w³. Here the transition dipole moment of A–X (d_AX) is much larger than that of the A–A′ (d_AA′): d_AX/d_AA′ ≈ 68;¹⁵ the frequency of A–X (w_AX) is much larger than that of the A–A′ (w_AA′): w_AX/w_AA′ ≈ 12.⁴¹ Thus γ_AX/γ_AA′ ≈ 10⁶. So the presence of the intervening A′²Δ state leak is unlikely to limit the laser cooling transition A–X significantly. The laser cooling of BaCl is similar to that of BaF. In order to avoid or reduce cumbersome problems, here we only show the proposed laser–driven transitions (solid red) and spontaneous decay (dotted line) with calculated f_ν′ν for BaCl in Fig. 3(b).

Fig. 3 (Color online) Proposed laser cooling schemes for (a) BaF and (b) BaCl using the A²Π_1/2 (ν′) → X²Σ⁺_1/2(ν) (solid red) transition. The decay pathways with calculated f_ν′ν are shown as dotted line.

4. Conclusions

In the present work, accurate ab initio calculations including the SOC effects are performed on barium monohalides at the MRCI + Q level of theory. The AVQZ basis set (F), pseudopotentials ECPnMDF (n = 46, 10, 28) (Ba, Br, I), ECP10MWB(Cl) and corresponding basis sets are used in the present calculations. The PECs for the A²Π and X²Σ⁺ states of BaX (X = F, Cl, Br, I) are presented and analyzed. Our calculated spectroscopic constants (R_e, w_e, w_eχ_e, B_e, T_e) of the X²Σ⁺_1/2 and A²Π_1/2 states are in good agreement with available experimental data. Using the PECs and transition dipole moments, we also obtain the FCFs and radiative lifetimes for the X²Σ⁺_1/2 → A²Π_1/2 transition. The results indicate that BaF and BaCl both have highly-diagonal FCFs and short radiative lifetime. Thus BaF and BaCl are promising candidates. The three-laser cooling schemes that drive the A²Π_1/2 → X²Σ⁺_1/2 transitions have been proposed for BaF and BaCl. However for BaI and BaBr, the A–X transitions do not have highly diagonal FCFs, and the calculated f₀₀ (BaI: 0.779; BaBr: 0.791) is not large enough to be potentially viable for direct laser cooling. BaI and BaBr are not identified as very promising laser-cooling candidates.

Conflicts of interest

There are no conflicts to declare.

Acknowledgements

This work was supported by the National Natural Science Foundation of China (Grant No. 11704052).

References

1. E. S. Shuman, J. F. Barry and D. DeMille, Nature, 2010, 467, 820.
2. J. F. Barry, D. J. McCarron, E. B. Norrgard, M. H. Steinecker and D. DeMille, Nature, 2014, 512, 286.
3. M. T. Hummon, M. Yeo, B. K. Stuhl, A. L. Collopy, Y. Xia and J. Ye, Phys. Rev. Lett., 2013, 110, 143001.
4. M. D. Di Rosa, Eur. Phys. J. D, 2004, 31, 395.
5. L. D. Carr, D. DeMille, R. V. Krems and J. Ye, New J. Phys., 2009, 11, 055049.
6. J. Kobayashi, K. Aikawa, K. Oasa and S. Inouye, Phys. Rev. A, 2014, 89, 021401(R).
7. V. Zhelyazkova, A. Cournol, T. E. Wall, A. Matsushima, J. J. Hudson, E. A. Hinds, M. R. Tarbutt and B. E. Sauer, Phys. Rev. A, 2014, 89, 053416.
8. T. A. Isaev, S. Hoekstra and R. Berger, Phys. Rev. A, 2010, 82, 052521.
9. N. Wells and I. C. Lane, Phys. Chem. Chem. Phys., 2011, 13, 19018.
10. Y. F. Gao and T. Gao, Phys. Rev. A, 2014, 90, 052506.
11. I. C. Lane, Phys. Chem. Chem. Phys., 2012, 14, 15078.
12. S. Y. Kang, Y. F. Gao, F. G. Kuang, T. Gao, J. G. Du and G. Jiang, Phys. Rev. A, 2015, 91, 042511.
13. H. Holweger and E. A. Mueller, Sol. Phys., 1974, 39, 19.
14. B. Karthikeyan, K. Balachandrakumar, V. Raja and N. Rajamanickamc, J. Appl. Spectrosc., 2013, 80, 5.
15. S. Y. Kang, F. G. Kuang, G. Jiang and J. G. Du, Mol. Phys., 2016, 114, 810.
16. L. E. Berg, N. Gador, D. Husain, H. Ludwigs and P. Royen, Chem. Phys. Lett., 1998, 287, 89.
17. C. Amiot, M. Hafid and J. Verges, J. Mol. Spectrosc., 1996, 180, 121.
18. A. Derkatch, C. Lundevall, L.-E. Berg and P. Royen, Chem. Phys. Lett., 2000, 332, 278.
19. D. Husain, J. Lei, F. Castano and M. N. Sánchez Rayo, J. Photochem. Photobiol. Chem., 1997, 107, 1.
20. R. P. Tuckett, personal communication, 1992.
21. W. E. Ernst, G. Weller and T. Torring, Chem. Phys. Lett., 1985, 121, 6.
22. E. Miliordos, A. Papakondylis, A. A. Tsekouras and A. Mavridis, J. Phys. Chem. A, 2007, 111, 10002.
23. H.-J. Werner, P. J. Knowles, R. Lindh, F. R. Manby and M. Schütz, et al., MOLPRO, is a package of ab initio programs, version 2009.1, http://www.molpro.net.
24. R. J. Le Roy, LEVEL 8.0, A Computer Program for Solving the Radial Schrödinger Equation for Bound and Quasibound Levels, CPRR-663, University of Waterloo, Waterloo, ON, 2007.
25. D. R. Bates, Mon. Not. Roy. Astron. Soc., 1952, 112, 614.
26. G. Herzberg, Spectra of Diatomic Molecules, Van Nostrand, New York, 1950.
27. P. F. Bernath, Spectra of Atoms and Molecules, Oxford University Press, New York, 2nd ed, 2005.
28. H.-J. Werner and P. J. Knowles, J. Chem. Phys., 1985, 82, 5053.
29. P. J. Knowles and H.-J. Werner, Chem. Phys. Lett., 1985, 115, 259.
30. H.-J. Werner and P. J. Knowles, J. Chem. Phys., 1988, 89, 5803.
31. P. J. Knowles and H.-J. Werner, Chem. Phys. Lett., 1988, 145, 514.
32. S. R. Langhoff and E. R. Davidson, Int. J. Quantum Chem., 1974, 8, 61.
33. N. Douglas and N. M. Kroll, Ann. Phys., 1974, 82, 89.
34. B. A. Hess, Phys. Rev. A, 1986, 33, 3742.
35. A. Berning, M. Schweizer, H.-J. Werner, P. J. Knowles and P. Palmieri, Mol. Phys., 2000, 98, 1823.
36. A. Bergner, M. Dolg, W. Kuechle, H. Stoll and H. Preuss, Mol. Phys., 1993, 80, 1431.
37. K. A. Peterson, D. Figgen, E. Goll, H. Stoll and M. Dolg, J. Chem. Phys., 2003, 119, 11113.
38. K. A. Peterson, B. C. Shepler, D. Figgen and H. Stoll, J. Phys. Chem. A, 2006, 110, 13877.
39. I. S. Lim, H. Stoll and P. Schwerdtfeger, J. Chem. Phys., 2006, 124, 034107.
40. A. K. Wilson, D. E. Woon, K. A. Peterson and T. H. Dunning Jr, J. Chem. Phys., 1999, 110, 7667.
41. R. F. Barrow, A. Bernard, C. Effantin, J. D'incan, G. Fabre, A. E. Hachimi, R. Stringat and J. Verges, Chem. Phys. Lett., 1988, 147, 535.
42. C. Ryzlewicz and T. Torring, Chem. Phys., 1980, 51, 329.
43. K. P. Huber and G. Herzberg, Molecular Spectra and Molecular Structure IV: Constants of Diatomic Molecules, Van Nostrand-Reinhold, 2005.
44. R. F. Gutterres, J. Vergěs and C. Amiot, J. Mol. Spectrosc., 1999, 196, 29.
45. R. F. Gutterres, J. Vergěs and C. Amiot, J. Mol. Spectrosc., 2000, 200, 253.
46. R. F. Gutterres, J. Vergěs and C. Amiot, J. Mol. Spectrosc., 2001, 206, 62.
47. R. Yang, Y. F. Gao, B. Tang and T. Gao, Phys. Chem. Chem. Phys., 2015, 17, 1900.

---