Theoretical Study of the Feasibility of Laser Cooling the - American

Mar 2, 2018 - College of Optoelectronic Technology, Chengdu University of Information Technology, Chengdu 610225, China. Key Laboratory of High ...
0 downloads 0 Views 1MB Size
Subscriber access provided by Kaohsiung Medical University

Article 24

35

Theoretical Study of The Feasibility of Laser Cooling The Mg Cl Molecule Including Hyperfine Structure and Branching Ratios Quan-Shun Yang, Shichang Li, You Yu, and Tao Gao J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/acs.jpca.7b11047 • Publication Date (Web): 02 Mar 2018 Downloaded from http://pubs.acs.org on March 4, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry A is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

24

Theoretical Study of The Feasibility of Laser Cooling The Mg35Cl Molecule Including Hyperfine Structure and Branching Ratios

Quan-Shun Yang, Shi-Chang Li, You Yu, Tao Gao∗ Institute of Atomic and Molecular Physics, Sichuan University, 610065, Chengdu, China College of Optoelectronic Technology, Chengdu University of Information Technology, Chengdu 610225, China Key Laboratory of High Energy Density Physics and Technology of Ministry of Education, Sichuan University, Chengdu, 610064, China

Abstract The possibility of laser cooling the 24Mg35Cl molecule is investigated using the electronic, rovibrational and hyperfine structure. Twelve low-lying Λ-S electronic states of the 24

Mg35Cl molecule have been calculated at the multi-reference configuration interaction

level of theory. The spin-orbit coupling effects are taken into account in the electronic structure calculations. Spectroscopic constants agree well with previously obtained theoretical and experimental values. Based on the potential energy curves and transition dipole moments, the highly diagonally distributed Franck-Condon factors for the A2Π → X2Σ+ transition and short radiative lifetime of the A2Π state are determined. Then, employing a quantum effective Hamiltonian approach, we investigate the hyperfine manifolds of the X2Σ+ state and obtain the zero-field hyperfine spectrum with the errors relative to the experimental data not exceeding 8 kHz ∼ 20 kHz. Finally, we design a laser cooling scheme with one cooling main laser beam and two repumping laser beams with modulated sidebands, which is sufficient for the implementation of efficient laser slowing and cooling of the 24Mg35Cl molecule. Moreover, it is important to note that the dissociation energy (2.2593 eV) of the B2Σ+ state is obtained for the first time at the multi-



Corresponding authors. Tel.:+86 028 85405234. E-mail: [email protected] (T. Gao) 1 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

reference configuration interaction level. We hope that this can provide a helpful reference for experimental observation.

1. Introduction The combination of laser cooling and trapping has produced ultracold temperature atoms, resulting in revolutionary advances in atomic clocks1, quantum information processing2 and simulation of condensed-matter systems3,4. With their additional types of internal motion (such as vibration and rotation), ultracold polar molecules yield even richer phenomena and applications in quantum simulations of condensed-matter system5, quantum computation6 and chemical dynamics7,8. However, cooling of molecules to cold and ultracold temperature has not been realized before 2010. One of the main difficulties in cooling molecules is the lack of a closed two-level transition cycle. This is due to the absence of rigorous selection rules controlling the branching ratios for the decay of an excited electronic state into different vibrational levels of the ground state. Fortunately, the direct laser cooling of SrF9 has paved the way to a new approach to this problem. In this case, the branching ratios are governed by the molecule’s Franck–Condon factors (FCFs). The completion realization of a full three-dimensional MOT for SrF was employed10, and a further improvement based on a MOT beam polarization scheme was derived from a detailed calculation by Tarbutt11. Then, one- and two-dimensional transverse laser cooling and MOT of YO12 and the longitudinal laser cooling of CaF13 molecule was successfully demonstrated experimentally. Finally, the YO molecule was cooled to 2 mK via Doppler cooling with the addition of magneto-optical trapping forces. Recently, the molecules had been cooled to below 1 mk by some group14-18. As early as 2004, Di Rosa presented a brief survey19 of candidate molecules for laser cooling. Therefore, theoretically, a large number of polar molecules20-23 have also been investigated for the feasibility of their use in laser cooling. As with atomic systems19, a potential laser cooling candidate molecule should meet certain conditions: a) highly diagonal Franck-Condon factors (FCFs) that limit the vibrational branching and the number of lasers, b) the excited state should exhibit a short radiative lifetime for a high photon scattering rate, and c) the absence of an intermediate electronic state between the

2 ACS Paragon Plus Environment

Page 2 of 28

Page 3 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

ground state and the first excited state or the intermediate state having no effect on the laser-cooling cycle. The properties of 24Mg35Cl have been studied theoretically24-27 and experimentally2835

by many researchers. For instance, the electronic structures for low-lying states of

24

Mg35Cl have been elucidated through the first-order configuration interaction (FOCI)

and second-order configuration interaction (SOCI) method by Yarkony et al.26 Then, the rotational spectra for the X2Σ+ state of the 24Mg35Cl molecule in various isotropic forms and vibrational states were observed in the 130-290 GHz frequency range by Bogey et al.29 Hirao et al.30 studied the A2Π–X2Σ+ emission band system of the 24Mg35Cl molecule by high-resolution Fourier transform spectroscopy in 2002. The emission bands systems of the 24Mg35Cl and 24Mg37Cl molecules were further studied by high resolution Fourier Transform Spectroscopy (FTS) by Gutterres et al.31 We exploit

24

Mg35Cl for the

following reasons: (i) Recently, laser cooling of 24Mg35Cl has been studied theoretically20. However, some problems are still not solved well, such as the relevant hyperfine structure (HFS) of

24

Mg35Cl, the modulating frequency covering all of the hyperfine

energy levels by the proposed sideband frequency distributions. (ii) Cold molecules have found great utility in precision measurements36. Thus, in the present work, we seek to perform a systematic theoretical investigation of the feasibility of laser cooling the 24

Mg35Cl molecule including hyperfine structure and branching ratios. The paper is organized as follows. In Sec. 2, the ab initio methods and basis set for

twelve low-lying Λ-S electronic states calculations of

24

Mg35Cl molecule are briefly

described. Sec. 3 shows the computational results. The hyperfine structure and hyperfine structure branching ratio calculation procedures are presented in Sec. 4 and Sec. 5. A brief summary is given in Sec. 6.

2. METHOD 3 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

For

24

Page 4 of 28

Mg35Cl, we calculate the potential energy curves (PECs) of twelve low-lying

Λ-S states. The complete active space self-consistent field (CASSCF)37, 38 mothed is used for the calculations. Afterward, the multi-reference configuration interaction (MRCI)39, 40 plus Davidson correction (+Q)41 calculations with the CASSCF wave functions as a zeroorder reference function are carried out. The calculations are performed in the C2v symmetry, where all molecular orbitals are labeled by (a1, b1, b2, a2) irreducible representations, and the electronic states are labeled by (A1, B1, B2, A2) irreducible representations. In the CASSCF calculations, the nine active (valence) electrons representing the Mg 3s2 and Cl 3s23p5 shells are placed into the active space. The twenty electrons in the 1s22s22p6 shell of the Mg and Cl atoms are placed into the closed spaces but are still optimized. The closed spaces are denoted (6,2,2,0) and the total orbital space is denoted (10,4,4,0). In this case, the active space is referred to as CAS (9, 8). Additionally, state averaged were computed using the CASSCF method which contained all possible electronic states (7A1, 5B1, 5B2, 1A2). We choose the same orbital space and closed spaces in the subsequent MRCI calculations. The correlation consistent polarized valence quadruple zeta cc-PVQZ-DK (=VQZDK)42 and augmented correlation consistent polarized valence quadruple zeta aug-ccPVQZ-DK (=AVQZ-DK)43 all-electron basis sets are used for the Mg and Cl atoms in the calculations of the Λ-S states, respectively. The scalar relativistic effects are accounted for using the Douglas–Kroll–Hess (DKH)44,

45

transformation of the relativistic

Hamiltonian. The spin-orbit coupling (SOC) effects are also considered using the SOC operator Hso defined within the Breit-Pauli approximation46. The same active space and basis sets are used in the SOC calculations at the MRCI theory level. All calculations are implemented with the MOLPRO47 software. The spectroscopic constants (Re, Te, ωe, ωeχe, Be, De) are obtained by using Le Roy’s LEVEL 8.048 program. Using the PECs and transition dipole moments (TDMs), the FCFs and radiative lifetimes of the various vibrational levels can also be determined from LEVEL

calculations48.

8.0

4 ACS Paragon Plus Environment

Page 5 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

FIG. 1. PECs of the low-lying double electronic states of the Σ+, Π and ∆ of 24Mg35Cl molecule at the MRCI level of theory.

FIG. 2. PECs of the

low-lying

quadruple

electronic states of the

Σ+ and Π of

molecule at the MRCI

level of theory.

24

Mg35Cl

3. Results and discussion 3.1. PECs and spectroscopic constants In Figs. 1 and 2, we present the PECs of twelve low-lying Λ-S states at the MRCI level of theory. The corresponding spectroscopic constants of these bound states calculated with the LEVEL 8.048 program are listed in Table 1. Previous theoretical and experimental values for these bound states are also tabulated in Table 1. Table 1. Spectroscopic constants of bound states of 24Mg35Cl calculated at the MRCI level of theory. state X2Σ+

Te(cm-1) 0 0 0

Re(Å) 2.19 2.20

ωe(cm-1) 483.20 462.12 467.53

ωeχe(cm-1) 2.24 2.10

5 ACS Paragon Plus Environment

Be(cm-1) 0.25 0.26 0.24

De(eV) 3.42 3.30

This work Expa Ref20

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

A2Π

22Π B2Σ+(int) B2Σ+(ext) C2Σ+ a Ref. 35. b c

0 0 0 26958.71 26062.04 26310.80 26741.33 26739.91 32363.35 31945.56 38613.09 30867.66 42918.81

2.17 2.20 2.17 2.18 2.16 2.17 2.52 2.55 2.15 3.66 2.37

469.90 466.08 462.10 515.92 492.33 487.60 490.80

1.95 2.00

0.25 0.25 0.25 0.25

0.55 0.54

22.98

0.25 0.25 0.19 0.18 0.26

2.07 2.01 0.10

0.21 16.22

0.09 0.21

2.26 0.77

2.20 1.27 2.18

681.20 622.72 540.39 552.00 179.32 705.16

Page 6 of 28

14.29

3.37

Ref 27 Expb Expc This work Ref 20 Ref 27 Expb Expc This work Ref 20 This work Expa This work This work

Ref. 34.

Ref. 31.

For the X2Σ+ state, Wan et al.20 found a slightly larger equilibrium distance Re of 2.2025 Å compared to the experimental data. Wan et al.20 did not report their result for the anharmonicity constant ωeχe for the X2Σ+ state. However, our calculations yield a shorter equilibrium distance (2.1940 Å) and a slightly larger anharmonicity constant (2.2378 cm-1) for the X2Σ+ state compared to the observed values (2.1960 Å and 2.1000 cm-1). The percentage errors for our calculated equilibrium distance and the anharmonicity constant are 0.09% and 6.56%, respectively. The Be values for the X2Σ+ and A2Π states are 0.2462 cm-1 and 0.2520 cm-1 in our calculation, in excellent agreement with the experimental values of 0.2456 cm-1 and 0.2517 cm-1. As shown in Table 1, the spectroscopic parameters for the A2Π state are also close to the experimental data, for example, the anharmonicity constant ωeχe of 2.1989 cm-1 is in good agreement with the experimental value of 2.1800 cm-1. Moreover, it is important to note that the dissociation energy of the B2Σ+ state is obtained here for the first time at the multi-reference configuration interaction level. The calculated harmonic constant ωe (540.3859 cm-1) for the B2Σ+(int) state is in accordance with the experimental data (552.0000 cm-1). These values will be very helpful as reference data for experimental observation. Overall, the spectroscopic constants derived from our present work are in good agreement with previous theoretical and experimental values.

6 ACS Paragon Plus Environment

Page 7 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Table 2. Spectroscopic constants of the first six Ω states of 24Mg35Cl at MRCI level.

A2Π3/2

Te(cm-1) 0 0 26888.87

A2Π1/2

26947.78

state X2Σ+1/2

2

2 Π3/2

32255.48

22Π1/2

32396.04

B2Σ+1/2(int) B2Σ+1/2(ext) a Ref. 35.

38616.28 30923.27

Re(Å) 2.19 2.20 2.17 2.18 2.17 2.18 2.17 2.52 2.54 2.54 2.57 2.16 3.66

ωe(cm-1) 483.55 467.52 518.27 492.79 519.18 495.04 491.60 681.20 616.72 720.03 635.54 525.75 179.96

ωeχe(cm-1) 2.17 2.44 2.65 2.54 14.29 17.12 22.46 0.21

Be(cm-1) 0.25 0.24 0.25 0.25 0.25 0.25 0.25 0.19 0.18 0.18 0.18 0.26 0.09

De(eV) 3.38 3.26 0.52 0.51 0.57 0.56 2.07 2.03 2.07 2.01 0.11 2.22

FIG. 3. PECs of six Ω states of 24Mg35Cl molecule at MRCI+Q level.

7 ACS Paragon Plus Environment

This work Ref 20 This work Ref 20 This work Ref 20 Expa This work Ref 20 This work Ref 20 This work This work

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

FIG. 4. The PDMs (a) of the double electronic states of the Σ+, Π, and ∆, (b) of the quadruple electronic states of the Σ+, and Π of 24Mg35Cl molecule.

The SOC effects are also taken into account for the 24Mg35Cl molecule. The PECs of six Ω states are shown in Fig. 3. The spectroscopic parameters together with previous theoretical values for the six Ω states are collected in Table 2. Comparison of the spectroscopic constants presented in Table 1 (with the SOC) and in Table 2 (without the SOC) shows that the equilibrium distance Re, anharmonicity constant ωexe and rotation constant Be of the X2Σ+1/2, B2Σ+1/2(int), B2Σ+1/2(ext), and 22Π3/2 states are almost identical to those of their X2Σ+, B2Σ+(int), B2Σ+(ext), and 22Π parents. For the harmonic constant ωe, the differences of X2Σ+, B2Σ+(ext), A2Π, and 22Π with X2Σ+1/2 B2Σ+1/2(ext), A2Π1/2, A2Π3/2, and 22Π3/2 are 0.3494 cm-1, 0.6448 cm-1, 3.2655 cm-1, 2.3559 cm-1 and 0.0033 cm-1, respectively. Therefore, the SOC did not show a significant impact for the spectroscopic parameters. The permanent dipole moments (PDMs) of the twelve low-lying electronic states and the TDMs of the 2Π → 2Σ+ transitions are presented in Figs. 4 and 5. As shown in Fig. 4, the calculated PDMs are decreasing as the internuclear distance R increases for X2Σ+, reaching a minimum (2.74 a.u.) at R > Re and finally tending to zero. The evaluated PDM data (3.167 Debye) at Re are in accordance with the previous theoretical value (3.381 Debye)28. For the A2Π state, the magnitude of the PDMs gradually decreases with increasing R and reaches a minimum (−1.82 a.u.) at R=2.30 Å. Subsequently, it starts to rise, reaches a maximum (0.95 a.u.) at R=2.65 Å, and drops thereafter, in reasonable agreement with the results obtained by Wan et al.20 8 ACS Paragon Plus Environment

Page 8 of 28

Page 9 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

FIG. 5. TDMs for the transition from A2Π, 22Π and 32Π states (a) to the ground state X2Σ+, (b) to the B2Σ+ state and (c) to the C2Σ+ state.

3.2. FCFs and radiative lifetime We obtain the FCFs of the A2Π → X2Σ+ transition and spontaneous radiative lifetimes of the A2Π state using the LEVEL 8.0 program48. The FCFs describe the overlap of the vibrational wave functions, and diagonal FCFs can restrain the vibrational branching of a state. The corresponding FCF data along with the previous theoretical values for the 24

Mg35Cl molecule are listed in Table 3 for the transition discussed here. The FCF of

~0.932 for the A2Π → X2Σ+ transition is slightly larger than the previously calculated value of ~0.928 obtained by Wan et al.20 and is in extremely close agreement with the value of ~0.932 predicted by Kumaran et al.49 using the closed-form approximation method. However, the FCFs are not the only important parameters for laser cooling molecules. A potential candidate molecule should also exhibit a rate of optical cycling that produces a large spontaneous scattering force for laser cooling. Examination of the 9 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 28

data presented in Table 4 shows that the A2Π (v′=0) state has a sufficiently short radiative lifetime of γ-1 = 12.28 ns, allowing for fast optical cycling. In addition, there is no intermediate state that can affect the laser cooling cycle of the

24

Mg35Cl molecule.

Therefore, the 24Mg35Cl molecule meets the three conditions of a potential candidate for laser cooling.

FIG. 6. Proposed laser-cooling scheme for 24Mg35Cl using the X2Σ+ → A2Π transition (solid lines) and spontaneous decay (dotted lines) with calculated fv'v″ for 24Mg35Cl. Solid short line indicates relevant upper and lower vibrational structure in 24Mg35Cl, respectively.

10 ACS Paragon Plus Environment

Page 11 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Table 3. The calculated FCFs fv'v″ along with theoretical values and wavelength (unit is nm) λν'ν″ of the A2Π → X2Σ+ transition. The diagonal FCFs are in bold. Numbers in parentheses indicate the power of 10.

Present work

The closed-form approximationa

Referenceb a

f00 f01 f02 f03 9.32(10-1) 6.50(10-2) 2.62(10-3) 9.51(10-5) 9.32(10-1) 6.53(10-2) 2.28(10-3) 6.43(10-5) 9.28(10-1) 6.90(10-2) 3.00(10-3) 1.40(10-4)

f10 f11 f12 f13 6.55(10-2) 8.00(10-1) 1.26(10-1) 8.00(10-3)

f20 f21 f22 f23 2.16(10-3) 1.28(10-1) 6.69(10-1) 1.84(10-1)

f30 f31 f32 f33 6.01(10-5) 7.09(10-3) 1.86(10-1) 5.42(10-1)

λ00

λ10

λ21

371.30

378.02

377.52

8.10(10-1) 1.21(10-1) 6.54(10-3) 7.00(10-2) 7.92(10-1) 1.27(10-1) 1.40(10-4)

Ref. 49.

b

Ref. 20.

A three-laser cooling scheme for the laser cooling of 24Mg35Cl is presented in Fig. 6. The vibrational and rotational branching loss must be addressed for the laser cooling cycle. Therefore, the FCF (0.932) of the X2Σ+ (v″=0) → A2Π (v ′ =0) transition is suitable as the main pump. The first vibrational repump (X2Σ+ (v″=0) → A2Π (v′=1)) and second vibrational repump (X2Σ+ (v″=1) → A2Π (v ′ =2)) are also required to obtain more than 104 scattered photons. Table 4. Spontaneous radiative lifetimes (ns) of A2Π → X2Σ+ transition of 24Mg35Cl molecule Radiative lifetime transition

A 2Π

v′=0

v′=1

v′=2

v′=3

v′=4

12.28

12.35

12.43

12.51

12.65

11 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 28

4. HYPERFINE STRUCTURE OF THE 24Mg35Cl RADICAL In laser cooling of molecules, it is important to consider the hyperfine structure (HF) of the ground electronic state for the elimination of dark states. Thus, all hyperfine levels should be pumped simultaneously in the case of the molecules accumulating into one state, and then, nearly closed optical transitions are implement. In this work, we have estimated the HFs of the ground state X2Σ+ of the

24

Mg35Cl molecule using the matrix

diagonalization technique. In a diatomic molecular system, the effective Hamiltonian is a sum of terms representing the various inter-coupling of angular momenta within the molecule50. Therefore, the effective Hamiltonian for the 24Mg35Cl ground state is given by H hf = γvT 1 ( Sˆ )T 1 ( Nˆ ) + bF T 1 ( Iˆ)T 1 ( Sˆ ) + cvTq1= 0 ( Iˆ)Tq1= 0 ( Sˆ ) − eT 2 (Q)T 2 (∇E )

(1)

where γv is the spin-rotational constant, bF is the Fermi contact constant, and cv is the dipole-dipole constant. All parameters36 are listed in Table 5. Table 5. Rotational constants and hyperfine constants for 24Mg35Cl (in MHz) v=0 B 7339.1014 D 8.1640×10-3 γ 66.547 bF 36.520 cv 36.463 eQq -11.622

A typical Hund’s case (b) case molecule first S (electron spin) couples with N (the total angular momentum excluding electron spin) to form J (total angular momentum). Then, J couples with I (total nuclear spin) to form the total angular momentum F (total angular momentum of molecule). The angular momentum states are written as | ((S, N) J, I) F)〉, and this is known as the bβJ basis. Based on the bβJ basis, the corresponding matrix elements for each term of the hyperfine structure are tabulated as follows: η′, Λ′, N ′, S , J ′, I , F | γvT 1 ( Sˆ )T 1 ( Nˆ ) | η, Λ, N , S , J , I , F S N J  = δJJ ′δFF ′ γv ( −1) N + J + S [ S ( S + 1)(2 S + 1)]1 / 2 [ N ( N + 1)(2 N + 1)]1 / 2    N S 1 12 ACS Paragon Plus Environment

(2)

Page 13 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

η′, Λ′, N ′, S , J ′, I , F | bF T 1 ( Iˆ)T 1 ( Sˆ ) | η, Λ, N , S , J , I , F = δFF ′bF ( −1) J ′ + F + I + N + J +1+ S [( 2 J ′ + 1)( 2 J + 1)]1 / 2 [ S ( S + 1)(2 S + 1)]1 / 2 [ N ( N + 1)( 2 N + 1)]1 / 2 (3)

I J ′ F  J S N  × ×    J I 1   S J′ 1

η′, Λ′, N ′, S , J ′, I , F | cvTq1= 0 ( Iˆ)Tq1= 0 ( Sˆ ) | η, Λ, N , S , J , I , F = δFF ′C (−1) J ′ + F + I + N [(2 J ′ + 1)(2 J + 1)]1 / 2 [ S ( S + 1)(2 S + 1)]1 / 2 [ N ( N + 1)(2 N + 1)]1 / 2 (2 N + 1) J N 2 N   ×  N ×  0 0 0   S 

J ′ 1  I J ′ F  N 2 ×   J I 1   S 1

(4)

η′, Λ′, N ′, S , J ′, I , F | −eT 2 (Q )T 2 (∇E ) | η, Λ, N , S , J , I , F  I J ′ F  ( I + 1)( 2 I + 1)(2 I + 3) eQq ( −1) J ′ + F + I  (−1) J ′ + N + S [( 2 J ′ + 1)( 2 J + 1)]1 / 2  2 J I 4 ( 2 − 1 ) I I   (5) J′ S N 2 N   N (2 N + 1) ( −1) ×  N 2  0 0 0 

= δ JJ ′ δ F F ′ N × J

where δii' is defined as δii' =0 for i ≠ i' and δii' =1 for i = i', and eQq is the nuclear quadrupole coupling constant of

35

Cl (ICl = 3/2). By numerically diagonalizing the

constructed matrix, we obtain the state eigenvectors and the energy eigenvalues. Then, using the energy eigenvalues of different rotational hyperfine energy levels, we obtain the transition frequencies among the rotational hyperfine levels. The experimental observational spectral data and our theoretically calculated results for the X2Σ+ (v″=0) state of the 24Mg35Cl molecule are presented in Table 6. The fourth column presents the theoretical calculation spectrum data, with the errors relative to the experimental data not exceeding 8 kHz ∼ 20 kHz. Therefore, the theoretical methods and results used in this work are reliable. However, the A2Π state is best described by Hund’s case (a); the parity of the rotational ladder in the A2Π states is given by half-integral J, levels with parities (1) J-1/2 or (-1) J+1/2 are designated + or -, respectively. The hyperfine splitting between F = 0 and F = 1 for the Hund’s case (a) state A2Π is still unresolved.

13 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 28

Table 6. The experimental spectral data and our theoretically computed results and their comparisons for 24Mg35Cl. The first three columns are the permissible transitions between the rotational hyperfine levels of the X2Σ+ state. The forth column is observational spectrum data in Ref. 36. Our theoretical results are included in the fifth column and the differences between experimental data and theoretically calculated results are given in the last column. N″→N′

J″→J′

0→1

1/2→1/2

0→1

1/2→3/2

F″→F′ 2→2 1→2 2→3 2→2 1→2 1→0

vexp (MHz) 14570.719 14643.751 14708.397 14690.968 14764.000 14720.436

vcal (MHz) 14570.701 14643.741 14708.378 14690.952 14763.992 14720.426

vexp-vcal (MHz) 0.018 0.010 0.019 0.019 0.008 0.010

Based on the calculated data, we designed the closure of the rotational structures for 24

Mg35Cl, as presented in Fig. 7. Clearly, we sought to cool on the N=0 rotational state

for the ground state (X2Σ+). However, due to parity and angular momentum selection rules, molecules can decay back to the N=0 and N=2 levels, and rotational branching opens up. For

24

Mg35Cl, this loss of N = 2 and N = 0 still exist. Therefore, we should

eliminate this loss by using the microwaves resonant with both |N = 0; J = 1/2; F = 0,1〉 ↔ |N = 1; J = 1/2; F =1,0〉 transitions. The rotational state (N=1) is excited to the |A, J=1/2, +〉 state with the wavelength of 371.30 nm. Moreover, six hyperfine energy levels (|J=3/2, F=3〉, |J=3/2, F =2〉, |J=3/2, F=1〉, |J=3/2, F=0〉, |J=1/2, F=2〉, and |J=1/2, F=1〉) for N=1 are obtained due to the I=3/2 nuclear spin of Cl, which splits J into F=J±1/2 levels through magnetic and electric hyperfine interaction. Meanwhile, each of the population levels requires a repump laser. However, experimentally, it is challenging to use so many lasers. Fortunately, we can use an electro-optical modulator (EOM) with the laser detuning of δ = −17.5 MHz and the modulation frequency of fMod = 35 MHz to cover all hyperfine energy levels by the proposed sideband frequency distributions. So far, we have discussed all of the conditions for the laser cooling of the

24

Mg35Cl molecule.

A2Π shows a short lifetime and highly diagonal FCFs in the vibrational level. We eliminated the rotational branching by parity selection rules and microwave remixing. Additionally, the sideband modulation will cover all six hyperfine levels of |X, N=1〉

14 ACS Paragon Plus Environment

Page 15 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

FIG. 7. The relevant hyperfine energy-level structure, laser frequencies, and spontaneous radiation (dashed lines) for 24Mg35Cl. Due to the angular momentum and parity selection rules a spontaneous decay back to the odd parity rotational state N=1 of X2Σ+ is allowed.

Table 7. The J-mixing phenomena in N=1 rotational states of the ground X state. N

1

Nominal label |J=3/2, F=3〉 |J=3/2, F=2〉 |J=1/2, F=2〉 |J=1/2, F=1〉 |J=3/2, F=1〉 |J=3/2, F=0〉

Actual label |J=3/2, F=3〉 0.3711|J=1/2, F=2〉+0.9286|J=3/2, F=2〉 0.9286|J=1/2, F=2〉-0.3711|J=3/2, F=2〉 0.9354|J=1/2, F=1〉+0.3535|J=3/2, F=1〉 -0.3535|J=1/2, F=1〉+0.9354|J=3/2, F=1〉 |J=3/2, F=0〉

5. BRANCHING RATIOS We calculate the branching ratios, which reflect the distributions of the transition strengths for all possible hyperfine decays from |A, J=1/2, +〉 to |X, N=1, −〉. These calculations are done in the Hund’s case (a) basis. Thus, we must transform the basis sets |A, J, +〉 and |X, N, J, F, −〉 to the Hund’s case (a) basis (|Λ, S, Σ, J, Ω, I, F, mF〉). We also should consider the mixing among the J states for the same N before deriving the matrix elements for the transition in A-X. The mixing is summarized in Table 7. The conversion for the pure J state between the Hund’s case (b) and the Hund’s case (a) is given by

15 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

N, J = N +

1 1  1 1 1 1  Σ = ,Ω = = − Σ = − ,Ω = −  2 2 2 2 2 2

N, J = N −

1 1  1 1 1 1  Σ = ,Ω = = + Σ = − ,Ω = −  2 2 2 2 2 2

Page 16 of 28

   

(6)

   

(7)

while for the A state, it is given by Π, J =

1 1  1 1 1 1 1 1 1 1  Λ = 1, S = , Σ = − , J = , Ω = ,+ = − Λ = −1, S = , Σ = , J = , Ω = −  2 2 2 2 2 2 2 2 2 2

   

(8)

Here, we calculate the transition dipole moments d = Ψe T p1 (dˆ ) Ψg

(9) ^

where |Ψe〉 and |Ψg〉 are Zeeman sublevels under the Hund’s case (a) basis. T1(d) is the first rank tensor denoting the electric dipole moment, and p denotes the polarization of light. Employing the Wigner-Eckart theorem, we obtain d = Λ′, S , Σ′, Ω′, J ′, I , F ′, m′F T p1 ( dˆ ) Λ, S , Σ, Ω, J , I , F , mF  F′ 1 F  J F I  (−1) F + J ′+ I +1 (2 F ′ + 1)( 2 F + 1)  = (−1) F ′−m    F ′ J ′ 1  − m′F p mF  F

(10)

× Λ′, S , Σ′, Ω′, J ′ T 1 (dˆ ) Λ, S , Σ, Ω, J

We can further reduce the final term Λ′, S , Σ ′, Ω′, J ′ T 1 ( dˆ ) Λ, S , Σ, Ω, J 1

=

∑( −1)

J ′− Ω′

q = −1

 J′ 1 J   × Λ′, S , Σ ′ T 1 ( dˆ ) Λ, S , Σ ( 2 J ′ + 1)(2 J + 1)   − Ω′ q Ω  ^

where the term 〈Λ′, S, Σ′ |T1(d)|Λ, S, Σ〉 is common to all branches.

16 ACS Paragon Plus Environment

(11)

Page 17 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Table 8. Branching ratios for the A2Π1/2, J′=1/2 → X2Σ+, N=1, i.e. the main optical cycling transition. Columns are denoted by the excited state’s hyperfine numbers (F′, m′F) and the rows are denoted by the ground states numbers (J, F, mF). J 3/2 3/2

1/2

3/2

1/2

3/2

F mF 0 0 -1 1 0 1 -1 1 0 1 -2 -1 2 0 1 2 -2 -1 2 0 1 2 -3 -2 -1 3 0 1 2 3

m′F=-1 1/18 0.1088 0.1088 0 0.0162 0.0162 0 0.0029 0.0014 0.0005 0 0 0.4138 0.2069 0.0690 0 0 0 0 0 0 0 0 0

F′=1 m′F=0 1/18 0.1088 0 0.1088 0.0162 0 0.0162 0 0.0014 0.0020 0.0014 0 0 0.2069 0.2759 0.2069 0 0 0 0 0 0 0 0

m′F=1 1/18 0 0.1088 0.1088 0 0.0162 0.0162 0 0 0.0005 0.0014 0.0029 0 0 0.0690 0.2069 0.4138 0 0 0 0 0 0 0

m′F=-2 0 0.0069 0 0 0.3431 0 0 0.1551 0.0775 0 0 0 0.1227 0.0614 0 0 0 1/6 1/18 1/90 0 0 0 0

m′F=-1 0 0.0035 0.0035 0 0.1715 0.1715 0 0.0775 0.0388 0.1163 0 0 0.0614 0.0307 0.0920 0 0 0 1/9 4/45 1/30 0 0 0

F′=2 m′F=0 0 0.0012 0.0046 0.0012 0.0572 0.2287 0.0572 0 0.1163 0 0.1163 0 0 0.0920 0 0.0920 0 0 0 1/15 1/10 1/15 0 0

m′F=1 0 0 0.0035 0.0035 0 0.1715 0.1715 0 0 0.1163 0.0388 0.0775 0 0 0.0920 0.0307 0.0614 0 0 0 1/30 4/45 1/9 0

m′F=2 0 0 0 0.0069 0 0 0.3431 0 0 0 0.0775 0.1551 0 0 0 0.0614 0.1227 0 0 0 0 1/90 1/18 1/6

Due to the angular momentum selection rules for the electric dipole transition, only the decays of |A, J=1/2, +〉 to the N = 1 level in the ground state are allowed. The calculated branching ratios for A2Π1/2, J′=1/2 → X2Σ+, and N=1 are listed in Table 8. As a consistency check, we sum the branching ratios in each column and obtain 1, as expected. The theoretically calculated spectra are then plotted with the branching ratio as the strength of

each peak (see

Fig. 8).

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

FIG. 8. The proposed sideband modulation scheme to the (re)pumping lasers to simultaneously cover all four hyperfine levels of |X, N=1〉. The theoretically calculated spectra (black line) are plotted with the branching ratio from the |A, J=1/2, +〉 to each hyperfine level as the strength of each peak, the calculated energy value in Fig. 7 as the center frequency. The red solid line indicates the sidebands of an EOM with the laser detuning δ = −17.5MHz and the modulated frequency fMod = 35MHz. The units of the x axis are MHz.

6. Conclusion In this work, the PECs and PDMs for the twelve low-lying Λ-S electronic states and TDMs for the 2Π →2Σ+ transition of the

24

Mg35Cl molecule have been assessed at the

CASSCF/MRCI (+Q) level of theory. The cc-PVQZ-DK (=VQZ-DK) and aug-cc-PVQZDK (=AVQZ-DK) all-electron basis sets are used for the Mg and Cl atoms, respectively, in the calculations of the Λ-S states. The SOC effects are also taken into account at the MRCI+Q level by using the same basis sets and active space. The SOC calculation results show that the influence of SOC on the spectroscopic properties is weak. Moreover, the hyperfine structure and branching ratio are also computed for the X2Σ+ ground state of the 24Mg35Cl molecule. The spectroscopic constants and transition frequencies derived from the present work are in good agreement with experimental values and available theoretical data. Using the potential energy curves and transition dipole moments, the highly diagonally distributed Franck-Condon factors for the A2Π → X2Σ+ transition and radiative lifetimes for the first six vibrational levels (v'=0-5) of the A2Π state are 18 ACS Paragon Plus Environment

Page 18 of 28

Page 19 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

determined. Our results reveal that the radiative lifetime of A2Π (v'=0) is 12.28 ns, which is sufficiently short for rapid laser cooling. There is no intermediate state affecting the laser cooling cycle of the

24

Mg35Cl molecule. Therefore, the three criteria for efficient

and rapid laser cooling are satisfied well. Finally, a three-laser cooling scheme for the laser cooling of 24Mg35Cl is presented. The X2Σ+ (v″=0) → A2Π (v′=0) transition can serve as the main pump due to its suitable FCF (0.932). The first vibrational repump (X2Σ+ (v″=0) → A2Π (v′=1)) and second vibrational repump (X2Σ+ (v″=1) → A2Π (v′ =2)) are also required to obtain more than 104 scattered photons. Each of the hyperfine population levels also requires a pump laser. However, the use of so many lasers will be challenging experimentally. Fortunately, we can use an EOM with the laser detuning of δ = −17.5 MHz and the modulation frequency of fMod = 35 MHz to cover all hyperfine energy levels by the proposed sideband frequency distributions. Overall, our results indicate that the 24Mg35Cl molecule is a promising candidate laser cooling molecule.

Acknowledgments The authors wish to thank the anonymous referee for valuable comments.

(1) Swallows, M. D.; Bishof, M.; Lin, Y.; Blatt, S.; Martin, M. J.; Rey, A. M.; Ye, J. Suppression of collisional shifts in a strongly interacting lattice clock. Science. 2011, 331, 1043-1046. (2) Saffman, M.; Walker, T. G.; Mølmer, K. Quantum information with Rydberg atoms. Reviews of Modern Physics. 2010, 82, 2313-2355. (3) Weber, T.; Herbig, J.; Mark, M.; Nägerl, H. C.; Grimm, R., Bose-Einstein condensation of cesium. Science. 2003, 299, 232-235. (4) Davis, K. B.; Mewes, M. O.; Andrews, M. R.; Van Druten, N. J.; Durfee, D. S.; Kurn, D. M.; Ketterle, W. Bose-Einstein condensation in a gas of sodium atoms. Physical review letters. 1995, 22, 3969-3972.

19 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(5) Pupillo, G.; Micheli, A.; Büchler, H.-P.; Zoller, P. in Cold Molecules: Theory, Experiment, Applications (eds Krems, R.; Friedrich, B.; Stwalley, W. C.), CRC Press, 2009. (6) DeMille, D. Quantum computation with trapped polar molecules. Phys. Rev. Lett. 2002, 88, 1-4. (7) Balakrishnan, N.; Dalgarno, A. Chemistry at ultracold temperatures. Chem. Phys. Lett. 2001, 341, 652–656. (8) Krems, R. V. Cold controlled chemistry. Phys. Chem. Chem. Phys. 2008, 10, 4079– 4092. (9) Shuman E. S.; Barry J. F.; DeMille D. Laser cooling of a diatomic molecule. Nature. 2010, 467, 820–823. (10) Barry, J. F.; McCarron, D. J.; Norrgard, E. B.; Steinecker, M. H.; DeMille, D. Magneto-optical trapping of a diatomic molecule. Nature. 2014, 512, 286-289. (11) Tarbutt, M. R. Magneto-optical trapping forces for atoms and molecules with complex level structures. New Journal of Physics. 2015, 1, 1-12. (12) Hummon, M. T.; Yeo, M.; Stuhl, B. K.; Collopy, A. L.; Xia, Y.; Ye, J. 2D magnetooptical trapping of diatomic molecules. Phys. Rev. Lett. 2013, 14, 1-5. (13) Wall, T. E.; Kanem, J. F.; Hudson, J. J.; Sauer, B. E.; Cho, D.; Boshier, M. G.; Hinds, E. A.; Tarbutt, M. R. Lifetime of the A (v′= 0) state and Franck-Condon factor of the A−X (0−0) transition of CaF measured by the saturation of laser-induced fluorescence. Phys. Rev. A. 2008, 6, 1-9. (14) Prehn, A.; Ibrügger, M.; Glöckner, R.; Rempe, G.; Zeppenfeld, M. Optoelectrical cooling of polar molecules to submillikelvin temperatures. Physical review letters, 2016, 6, 1-5. (15) Kozyryev, I.; Baum, L.; Matsuda, K.; Augenbraun, B. L.; Anderegg, L.; Sedlack, A. P.; Doyle, J. M. Sisyphus laser cooling of a polyatomic molecule. Physical Review Letters, 2017, 17, 1-5. (16) Truppe, S.; Williams, H. J.; Hambach, M.; Caldwell, L.; Fitch, N. J.; Hinds, E. A.; Tarbutt, M. R. Molecules cooled below the Doppler limit. Nature Physics. 2017, 13, 1173-1176.

20 ACS Paragon Plus Environment

Page 20 of 28

Page 21 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(17) Norrgard, E. B.; McCarron, D. J.; Steinecker, M. H.; Tarbutt, M. R.; DeMille, D. Submillikelvin dipolar molecules in a radio-frequency magneto-optical trap. Physical review letters, 2016, 6, 1-5. (18) Cheng, C.; van der Poel, A. P.; Jansen, P.; Quintero-Pérez, M.; Wall, T. E.; Ubachs, W.; Bethlem, H. L. Molecular fountain. Physical review letters, 2016, 25, 1-5. (19) DiRosa, M. D. Laser-cooling molecules. Eur. Phys. J. D. 2004, 31, 395-402. (20) Wan, M.; Shao, J.; Gao, Y.; Huang, D.; Yang, J.; Cao, Q.; Jin, C.; Wang, F. Laser cooling of MgCl and MgBr in theoretical approach. J. Chem. Phys. 2015, 2, 1-7. (21) Gao, Y.; Gao, T. Laser cooling of the alkaline-earth-metal monohydrides: Insights from an ab initio theory study. Phys. Rev. A. 2014, 5, 1-10. (22) Isaev, T.A.; Hoekstra, S.; Berger, R., Laser-cooled RaF as a promising candidate to measure molecular parity violation. Phys. Rev. A. 2010, 5, 1-5. (23) Xu, L.; Yin, Y.; Wei, B.; Xia, Y.; Yin, J. Calculation of vibrational branching ratios and hyperfine structure of

24

Mg19F and its suitability for laser cooling and magneto-

optical trapping. Physical Review A. 2016, 1, 1-10. (24) Montagnani, R.; Riani, P.; Salvetti, O. Theoretical study of the potential energy curves of the diatomic radicals Me II X. III. Application to MgCl, CaF and CaCl radicals and some preliminary remarks. Theoretical Chemistry Accounts: Theory, Computation, and Modeling (Theoretica Chimica Acta), 1983, 1, 13-19. (25) Langhoff, S. R.; Bauschlicher Jr, C. W.; Partridge, H. Theoretical dissociation energies for the alkali and alkaline‐earth monofluorides and monochlorides. The Journal of Chemical Physics, 1986, 3, 1687-1695. (26) Parlant, G.; Rostas, J.; Taieb, G.; Yarkony, D. R. On the electronic structure and dynamical aspects of the predissociation of the A2ΠΩ states of MgCl. A rigorous quantum mechanical treatment incorporating spin–orbit and derivative coupling effects. J. Chem. Phys. 1990, 93, 6403-6418. (27) Langhoff, S. R.; Bauschlicher Jr, C. W.; Partridge, H.; Ahlrichs, R. Theoretical study of the dipole moments of selected alkaline‐earth halides. physics, 1986, 9, 5025-5031.

21 ACS Paragon Plus Environment

The Journal of chemical

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(28) Parker, A. E. Band Systems of MgCl, CaCl and SrCl. Physical Review. 1935, 5, 349358. (29) Bogey, M.; Demuynck, C.; Destombes, J. L. Millimeter wave spectrum of MgCl X2Σ+ and isotopomers in different vibrational states. Determination of mass-invariant parameters. Chemical physics letters. 1989, 3, 265-268. (30) Hirao, T.; Bernath, P. F.; Fellows, C. E.; Gutterres, R. F.; Vervloet, M., HighResolution Fourier Transform Study of MgCl: The A2Π–X2Σ+ Band System. J. Mol. Spectrosc. 2002, 212, 53-56. (31) Gutterres, R. F.; Santos, R. F. D.; Fellows, C. E. Spectroscopic study of the 24Mg35Cl and

24

Mg37Cl A²pi-X²sigma+ band system. Brazilian Journal of Physics. 2003, 4, 886-

891. (32) Morgan, F. Band spectra of MgCl, MgBr and MgI in absorption. Physical Review. 1936, 7, 603-607. (33) Rostas, J.; Shafizadeh, N.; Taieb, G.; Bourguignon, B.; Prisant, M. G., The a 2ΠX2Σ+ system of MgCl. Evidence for predissociation in the a 2Π state of MgCl. Chem. Phys. 1990, 142, 97-109. (34) Huber, K. P.; Herzberg, G., Molecular Spectra and Molecular Structure IV Constants of Diatomic Molecules; Van Nostrand Reinhold, New York, 1979. (35) Ohshima, Y.; Endo, Y. Fourier-transform microwave spectroscopy of 24Mg35Cl generated by laser ablation. Chemical physics letters. 1993, 213, 95-100. (36) Carr, L. D.; DeMille, D.; Krems, R. V.; Ye, J. Cold and ultracold molecules: science, technology and applications. New J. Phys. 2009, DOI: 10.1088/1367-2630/11/5/055049. (37) Werner, H. -J.; Knowles, P. J. A second order multiconfiguration SCF procedure with optimum convergence. J. Chem. Phys. 1985, 11, 5053-5063. (38) Knowles, P. J.; Werner, H. -J. An efficient second-order MC SCF method for long configuration expansions. Chem. Phys. Lett. 1985, 3, 259-267. (39) Werner, H. -J.; Knowles, P. J. An efficient internally contracted multiconfiguration– reference configuration interaction method. J. Chem. Phys. 1988, 9, 5803-5814. (40) Knowles, P. J.; Werner, H. -J. An efficient method for the evaluation of coupling coefficients in configuration interaction calculations. Chem. Phys. Lett. 1988, 6, 514-522.

22 ACS Paragon Plus Environment

Page 22 of 28

Page 23 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(41) Laughoff, S. R.; Davidson, E. R. Configuration interaction calculations on the nitrogen molecule. Int. J. Quantum Chem. 1974, DOI: 10.1002/qua.560080106. (42) Prascher, B.; Woon, D. E.; Peterson, K. A.; Dunning, T. H.; Wilson, A. K. Gaussian basis sets for use in correlated molecular calculations. VII. Valence, core-valence, and scalar relativistic basis sets for Li, Be, Na, and Mg. Theor. Chem. Acc. 2011, DOI: 10.1007/s00214-010-0764-0. (43) Woon, D. E.; Dunning Jr, T. H. Gaussian basis sets for use in correlated molecular calculations. III. The atoms aluminum through argon. J. Chem. Phys. 1993, 98, 13581371. (44) Douglas, N.; Kroll, N. M. Quantum electrodynamical corrections to the fine structure of helium. Ann. Phys. 1974, 1, 89-155. (45) Hess, B. A. Relativistic electronic-structure calculations employing a twocomponent no-pair formalism with external-field projection operators. Phys. Rev. A. 1986, 33, 3742-3748. (46) Berning, A.; Schweizer, M.; Werner, H.-J.; Knowles, P. J.; Palmieri, P. Spin-orbit matrix elements for internally contracted multireference configuration interaction wavefunctions. Mol. Phys. 2000, 98, 1823-1833. (47) Werner, H. J.; Knowles, P.; Knizia, G.; Manby, F. R.; Schütz, M.; Celani, P.; Korona, T.; Lindh, R.; Mitrushenkov, A.; Rauhut, G. Molpro, version 2010.1, a package of ab initio programs; 2010. (48) Le Roy, R. J., Level 8.0: A Computer Program for Solving the Radial Schrödinger Equation for Bound and Quasibound Levels, Research Report CP- 663; University of Chemical Physics: 2007. (49) Nicholls, R. W. Franck-Condon factor formulae for astrophysical and other molecules. J. Suppl. Ser. 1981, 47, 279-290. (50) Brown, J. H.; Carrington, A. Rotational Spectroscopy of Diatomic Molecules; Cambridge University Press, Cambridge, U.K., 2012.

TOC graphic 23 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

24 ACS Paragon Plus Environment

Page 24 of 28

Page 25 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

FIG. 1. PECs of the low-lying double electronic states of the Σ+,  and Δ of 24Mg35Cl molecule at the MRCI level of theory.

FIG. 2. PECs of the low-lying quadruple electronic states of the Σ+ and  of 24Mg35Cl molecule at the MRCI level of theory.

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

FIG. 3. PECs of six Ω states of 24Mg35Cl molecule at MRCI+Q level.

FIG. 4. The PDMs (a) of the double electronic states of the Σ+, , and Δ, (b) of the quadruple electronic states of the Σ+, and  of 24Mg35Cl molecule.

ACS Paragon Plus Environment

Page 26 of 28

Page 27 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

FIG. 5. TDMs for the transition from A2Π, 22Π and 32Π states (a) to the ground state X2Σ+, (b) to the B2Σ+ state and (c) to the C2Σ+ state.

FIG. 6. Proposed laser-cooling scheme for 24Mg35Cl using the X2Σ+ → A2Π transition (solid lines) and spontaneous decay (dotted lines) with calculated fv'v″ for 24Mg35Cl molecule. Solid short line indicates relevant upper and lower vibrational structure in 24Mg35Cl, respectively.

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

FIG. 7. The relevant hyperfine energy-level structure, laser frequencies, and spontaneous radiations for 24Mg35Cl

molecule. Due to the angular momentum and parity selection rules allow a spontaneous decay

only back to the odd parity rotational state N=1 of X2Σ+.

FIG. 8. The proposed sideband modulation scheme to the (re)pumping lasers to simultaneously cover all four hyperfine levels of |X, N=1. The theoretically calculated spectra (black line) are plotted with the branching ratio from the |A, J=1/2, + to each hyperfine level as the strength of each peak, the calculated energy value in Fig. 7 as the center frequency. The vertical solid color line indicates the sidebands of an EOM with the laser detuning δ = −17.5MHz and the modulated frequency fMod = 35MHz. The units of the x axis are MHz.

ACS Paragon Plus Environment

Page 28 of 28