Ming-jie Wan*ab,
Guo-sen Wangc,
Xing-yong Huang
ab,
Duo-hui Huang*ab and
Kang-lin Weid
aFaculty of Science, Yibin University, Yibin, China. E-mail: wanmingjie1983@sina.com; hdhzhy912@163.com
bComputational Physics Key Laboratory of Sichuan Province, Yibin University, Yibin, China
cSchool of Mathematics and Physics, Chengdu University of Technology, Chengdu, China
dIntelligent Manufacturing Institute of Yibin University, Yibin, China
First published on 28th November 2022
In this paper, the potential energy curves of 22 Λ–S states as well as 51 Ω states were calculated using the internally contracted multiconfiguration interaction and Davidson correction method. Through the obtained transition data, the spectroscopy data of the low excitation bound state are fitted and compared with the same main group ions. The phenomenon of avoided crossing that occurs in the Ω state is analyzed, and finally it is concluded that this phenomenon mainly occurs in the energy region between 20000 cm−1 and 40
000 cm−1. The potential laser cooling transition cycle in the Ω state is analyzed. The Franck–Condon factor, radiative lifetime and Einstein coefficient between
are calculated. In this paper, we argue that direct laser cooling of SeBr+ is not feasible. The content of our study provides a theoretical basis for subsequent calculations to explore the properties of SeBr+ spectrum.
In 2010, Shuman12 at Yale University reported an experiment using three laser beams to prepare SrF molecules with temperatures close to supercooling, which first confirmed the feasibility of using lasers to cool molecules directly in experimental. The success of experimental laser cooling of neutral molecules has led many researchers to turn their attention to anions and cations, and people have begun to investigate whether cations are also capable of laser cooling, so that the spectroscopic and transition properties of cations can be studied. The spectroscopic properties of molecules with the same electronic structure as SeBr+ have also been extensively studied. NCl,13 PCl,14 AsCl,15,16 SbCl,17 AsI,18,19 SbBr,20 SbI21 all have the same electronic state structure as SeBr+, which means they may have similar spectroscopic properties. For SCl+ of the same main group, Lu et al.22 calculated the spectroscopic data of the ground and low excited states of this ion, detailed potential energy curves in the Ω state and the potential laser-cooled transition in the Ω state. Wang et al.23,24 calculated the spectroscopic data of SeCl+ and TeCl+ in the Λ–S state as well as in the Ω state, compared in detail the spectroscopic properties in the presence of the same main group of ions, and concluded from detailed calculations that direct laser cooling of SeCl+ and TeCl+ is not feasible. Meanwhile, molecules with the same electronic state structure and the same potential energy curve trend as SeBr+ have also been extensively studied, and Belinassi25 investigated the transition properties between the ground state to the lower excited state of the SeI molecule. In addition, the spectrum properties of AsBr,18,26 a molecule with the same electronic number and electronic state structure as SeBr+, have also been reported. However, the spectroscopic properties and transition properties of SeBr+ are rarely reported, so this paper will investigate the spectroscopic properties and transition properties of SeBr+ and discuss the feasibility of using laser cooling of SeBr+. Discuss the phenomenon of avoided crossing in the Ω states, which means that the potential energy curves of two states with the same symmetry cannot intersect in a diatomic molecule. Since the relative atomic masses of both selenium and bromine are large, the spin–orbit coupling effect is not negligible. The phenomena of avoided crossing in the Ω state will also be more complicated. The obtained calculation results will be compared with ions of the same main group to analyze the similarities and differences of the spectroscopic properties. This would be the first report on SeBr+ molecular spectroscopy data, bridging the experimental and theoretical void.
In this paper, we present the research work in the following sections: the calculation details are shown in Section 2, the results and discussions are shown in Section 3, and the conclusions are shown in Section 4.
In the calculations, the initial guessed molecular orbitals and wave functions of the ground state of the SeBr+ ion were first obtained using the Hartree–Fock (HF) method.28,29 On this basis, we used the state-averaged complete active space self-consistent field (SA-CASSCF)30,31 for optimization to obtain the state-averaged multistate wave functions. The internally contracted multi-reference configuration interaction method (icMRCI)32,33 was then used to calculate the correlation energy, and the Davidson correction (+Q) was included for optimization. To improve the accuracy of the potential energy curves for the Λ–S and Ω states, the aug-cc-pVQZ-PP basis set34 was used for both Se+ and Br. The choice of these basis sets implies that a relativistic effective core potential is used in this paper when considering scalar-relativistic effects. In the CASSCF calculation, for SeBr+, we chose 8 molecular orbitals (4, 2, 2, 0) as the active orbitals, denoted as CAS(12,8). These 8 molecular orbitals include the 4s4p orbital of Se+ and the 4s4p orbital of Br. In the calculation of MRCI+Q, we consider the CV correlation effect, where the 3d orbitals of Se+ are used as the correlation level. When considering the spin–orbit coupling effect, the spin–orbit matrix elements and eigenstates are calculated using the spin–orbit pseudopotential.
With our obtained Λ–S and Ω state potential energy curves and calculated transition dipole moments, the one-dimensional radial Schrödinger equation is solved using the LEVEL 8.0 program35 to obtain the equilibrium nuclear bond lengths (Re), harmonic frequency (ωe), first anharmonic correction (ωeχe), rotation constants (Be) and vertical transition energies (Te) of the bound states.
Atomic states | Molecular states | This work | Exp.38,39 |
---|---|---|---|
Se+(4Su) + Br(2Pu) | X3Σ−, 13Π, 15Σ−, 15Π | 0 | 0 |
Se+(4Du) + Br(2Pu) | 11Σ+, 11Σ−, 21Σ−, 11Π, 21Π, 31Π, 11Δ, 21Δ, 11Φ, 13Σ+, 23Σ−, 33Σ−, 23Π, 33Π, 43Π, 13Δ, 23Δ, 13Φ | 13![]() |
13![]() |
Based on the obtained potential energy curves, we fitted the spectroscopic data of the bound states and show some of them in Table 2. (The spectrum data of all electronic states are shown in Table SI1.†) We will focus on the analysis of the first three electronic states in comparison with the electronic states corresponding to the ions of the same main group. Table 3 shows the spectrum constants of the low excited electronic states of ions in the same main group as SeBr+. The equilibrium bond lengths of SCl+,22 SeCl+,23 TeCl+24 show an increasing trend, while the rest of the spectrum constants show a decreasing trend. The equilibrium bond lengths of the low excited electronic states of SeBr+ are larger than those of SeCl+ and smaller than those of TeCl+. This phenomenon is consistent with the fact that the relative molecular mass of SeBr+ is between SeCl+ and TeCl+. Besides, the rest of the spectrum constants of SeBr+ are smaller than those of SCl+, SeCl+ and TeCl+. Comparing the vertical transition energies of 11Δ and 11Σ+, we found that the overall vertical transition energies of SeCl+ are smaller compared with SCl+, SeCl+, TeCl+. For some neutral molecules with similar electronic state structures as SeBr+, the vertical transition energies from the ground state (X3Σ−) to the first (11Δ) and second (11Σ+) excited states of SbBr36 and AsBr26 are smaller than those of SeBr+, a property similar to that of ions of the same main group. Besides, the spectrum properties of 13Δ and 13Σ+ are very similar, and the difference between the values of various spectrum data is very small.
Λ–S states | Te/cm−1 | Re/Å | ωe/cm−1 | Be/cm−1 | ωeχe/cm−1 |
---|---|---|---|---|---|
X3Σ− | 0 | 2.1802 | 378.5 | 0.0905 | 1.0977 |
11Δ | 5733 | 2.1776 | 381.7 | 0.0907 | 1.0549 |
11Σ+ | 10![]() |
2.1890 | 366.1 | 0.0897 | 1.2493 |
15Π | 19![]() |
2.9854 | 130.5 | 0.0482 | 1.1235 |
11Σ− | 20![]() |
2.4452 | 257.8 | 0.0719 | 0.8183 |
13Π | 20![]() |
2.9313 | 104.4 | 0.0500 | 0.8814 |
13Δ | 21![]() |
2.4542 | 254.7 | 0.0714 | 0.8267 |
Λ–S states | Monovalent molecules | Te/cm−1 | Re/Å | ωe/cm−1 | Be/cm−1 | ωeχe/cm−1 |
---|---|---|---|---|---|---|
X3Σ− | SCl+22 | 0 | 1.8779 | 686.1 | 0.2839 | 3.12 |
SeCl+23 | 0 | 2.0351 | 520.8 | 0.1686 | 1.9586 | |
TeCl+24 | 0 | 2.2287 | 448.9 | 0.1237 | 1.5727 | |
SeBr+ | 0 | 2.1802 | 378.5 | 0.0905 | 1.0977 | |
11Δ | SCl+22 | 6334 | 1.8763 | 691.0 | 0.2844 | 2.89 |
SeCl+23 | 6415 | 2.0310 | 527.7 | 0.1693 | 1.8346 | |
TeCl+24 | 6183 | 2.2244 | 455.0 | 0.1241 | 1.4745 | |
SeBr+ | 5733 | 2.1776 | 381.7 | 0.0907 | 1.0549 | |
11Σ+ | SCl+22 | 11![]() |
1.8818 | 675.1 | 0.2828 | 3.22 |
SeCl+23 | 11![]() |
2.0374 | 513.3 | 0.1682 | 2.1417 | |
TeCl+24 | 11![]() |
2.2289 | 446.6 | 0.1236 | 1.6931 | |
SeBr+ | 10![]() |
2.1890 | 366.1 | 0.0897 | 1.2493 |
The main electronic configurations at the equilibrium bond lengths are listed in Table 4. (The electronic configuration of all electronic states at the equilibrium nuclear spacing is shown in Table SI2.†) The electronic configuration of the ground state (X3Σ−) in the equilibrium nuclear bond length attachment is mainly 11σ212σ04π45π46π22δ2, accounting for 83.35%. Similar to the ground state electronic configuration is 11Δ, with a percentage of 79.55%. The electronic configuration of 13Δ and 13Σ+ in the equilibrium nuclear bond length attachment is 11σ212σ04π45π36π32δ2 with the percentages of 87.54% and 87.49%, respectively, a difference of only 0.27%. The reason for this phenomenon is due to the fact that the two states are very close in energy and have similar properties at the equilibrium nuclear bond lengths. The electronic configuration of 33Π and 43Π is also very similar, but the difference in the percentage is 3.59%.
Λ–S states | Main CSFs at Re (%) | Λ–S states | Main CSFs at Re (%) |
---|---|---|---|
X3Σ− | 11σ212σ04π45π46π22δ2 (83.35) | 21Δ | 11σ212σ04π45π36π32δ2 (70.75) |
11σ112σ14π45π46π22δ2 (8.96) | |||
11Δ | 11σ212σ04π45π46π22δ2 (79.55) | 23Π | 11σ112σ04π45π46π32δ2 (86.0) |
11Σ+ | 11σ212σ04π45π46π22δ2 (74.12) | 33Π | 11σ212σ14π45π36π22δ2 (81.29) |
11σ212σ04π45π36π32δ2 (9.0) | |||
13Π | 11σ212σ14π45π46π12δ2 (66.73) | 43Π | 11σ212σ14π45π36π22δ2 (84.70) |
11σ212σ14π45π36π22δ2 (15.59) |
Atomic state | Ω states | ΔE | |
---|---|---|---|
This work | Exp.38,39 | ||
a At the nuclear distance R = 7 Å. | |||
Se+(4S3/2) + Br(2P3/2) | 3, 2, 2, 1, 1, 1, 0+, 0+, 0−, 0− | 0.0 | 0.0 |
Se+(4S3/2) + Br(2P1/2) | 2, 1, 1, 0+, 0− | 3392.3 | 3685.2 |
Se+(4D3/2) + Br(2P3/2) | 3, 2, 2, 1, 1, 1, 0+, 0+, 0−, 0− | 13![]() |
13![]() |
Se+(4D5/2) + Br(2P3/2) | 4, 3, 3, 2, 2, 2, 1, 1, 1, 1, 0+, 0+, 0−, 0− | 13![]() |
13![]() |
Se+(4D3/2) + Br(2P1/2) | 2, 1, 1, 0+, 0− | 16![]() |
17![]() |
Se+(4D5/2) + Br(2P1/2) | 3, 2, 2, 1, 1, 0+, 0− | 17![]() |
17![]() |
The potential energy curves of the electronic states corresponding to these dissociation channels are shown in Fig. 4. As can be seen from the figure, there is a complex phenomenon of avoided crossing in the Ω state. For the Ω = 0+ state, the phenomenon of avoided crossing occurs mainly in the region between 20000 cm−1 and 35
000 cm−1 of energy. The phenomenon of avoided crossing between (3)0+ and (4)0+ is most evident at the nuclear bond length of 2.6 Å. The phenomenon of avoided crossing for the electronic states with Ω = 0− occur mainly between (2)0− and (3)0−, and (7)0− and (8)0−. For the electronic state with Ω = 1, there is a very complex phenomenon of avoided crossing, which occurs mainly around 25
000 cm−1 to 30
000 cm−1. This phenomenon occurs mainly between nuclear bond lengths of 2.5–3 Å and is more pronounced between some electronic states, for example, between (5)1–(9)1. The phenomenon of avoided crossing between (3)2 and (4)2, (5)2 and (6)2 is very obvious and occurs near the nuclear bond length of 3 Å. For the rest of the electronic states with Ω = 2, it is not easy to observe the phenomenon of avoided crossing, so it will not be discussed in detail. For the state Ω = 3, 4, there is no obvious phenomenon of avoided crossing, so we will not discuss it in detail. Based on the obtained potential energy curves, we obtained the spectroscopic data of some of the lower excited states by solving the one-dimensional radial Schrödinger equation and show them in Table 6. The vertical transition energy from the first excited state(X3Σ−1(1)1) to the ground state
is 418 cm−1, compared with 13 cm−1, 191 cm−1, and 965 cm−1 for SF+,37 SeCl+,23 and TeCl+24 of the same main group, respectively. It can be noticed that this value becomes larger and larger as the relative atomic mass increases. The vertical transition energy of 11Δ2(1)2 also changes somewhat, with a value of 5733 cm−1 in the Λ–S state and 6074 cm−1 in the Ω state. The change in this value indicates that the spin–orbit coupling effect has some influence on SeBr+. The spectroscopic data of the ground state (X3Σ−) did not change significantly after considering the spin–orbit coupling effect. For 11Δ, the vertical transition energy increases and the rest of the spectrum data do not change significantly. The changes before and after the 11Σ+ splitting are mainly reflected in the increase of vertical transition energy and first anharmonic correction, while the equilibrium nuclear bond length and the rest of the spectral data do not change much. In this paper, the SO matrix elements between the SeBr+ electronic states are calculated and the strength of the spin–orbit coupling interaction is determined by their magnitudes, and the possible 29 SO matrix elements between the calculated electronic states are listed in Table SI3.† SO1 and SO2 denote the SO matrix elements between the two components of 11Π and 13Σ+ so they are equal in size. Similarly, SO3 = SO4, SO5 = SO6, SO7 = SO8, SO9 = SO10, SO11 = SO12, SO13 = SO14, SO15 = SO16, SO17 = SO18, SO19 = SO20, SO21 = SO22, SO23 = SO24, SO25 = SO26, SO27 = SO28. The variation of the obtained SO matrix elements with the nuclear bond length is shown in Fig. 5. As a whole, the values of SO1, SO3, SO19, SO21 and SO25 are relatively small, while the rest of the values other than these are relatively large, with the maximum even reaching 1200 cm−1, which is sufficient to show that the rest of the states in the spin–orbit coupling effect are more affected. To ensure the accuracy of the calculation results, the possibility of laser cooling of SeBr+ under the effect of spin–orbit coupling will be discussed in this paper.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d2ra06580k |
This journal is © The Royal Society of Chemistry 2022 |