1 Introduction
The fields of physical chemistry or chemical physics have seen astonishing strides towards the creation of cold and ultracold molecules in electronic and rovibrational ground state in recent two decades. Such researches have fostered a wealth of interdisciplinary explorations, such as many-body quantum physics and chemistry [
1,
2], quantum computation [
3,
4], quantum simulation [
5–
7], cold and ultracold chemistry [
8–
13], and precision measurement [
14–
20]. Although a couple of diatomic molecules have already been successfully cooled to ultracold regime, these simple molecules are lack of potential to explore rich chemistry in a general sense. In comparison, free radicals involving in many crucial chemical reactions, but yet to be cooled to submillikelvin temperature. For example, methyl radical CH
3, the simplest organic polyatomic radical, is one of the most important and fundamental intermediates in hydrocarbon chemistry. It plays a key role in various reactions including combustion, atmospheric and interstellar chemistry. Creating ultracold CH
3 would help to understand the quantum mechanisms related to many elementary reactions. For example, at very low temperature, two types of reactions could happen, one is barrierless reaction CH
3 + OH
CH
2O + H
2, and the other is tunnelling process, CH
3 + H
2 CH
4 + H [
21,
22]. However, their reaction rate and branching ratio are still ambiguous [
23]. Understanding these reactions would give a thrust to the advancement of cold chemistry.
CH
3 molecule has an unpaired electron, and has a linear Zeeman shift in strong magnetic fields because its spin is decoupled from the molecular axis and preferentially oriented relative to the external field. An efficient cooling method for CH
3 is translational deceleration and trapping by a well-designed time-varying magnetic [
24,
25], which can be followed by a second-stage cooling [
26]. However, one of the fundamental requirement for the second-stage cooling to proceed is a much larger molecular density than what can be achieved by the usual decelerator. In a recent work, CH
3 molecules with a translational temperature of 200 mK are obtained and are trapped for more than 1 s, but the estimated density is on the order of 5.0×10
7 cm
−3 [
27]. This density is not large enough to ensure further cooling of CH
3 molecules to ultracold regime even if assuming favourable ratio of elastic to inelastic collision cross section between CH
3 molecules.
2 Co-deceleration scheme
In order to improve the CH
3 density, and also inspired by the recent studies on the collision of the O
2–Li mixture [
28], here we propose to simultaneously decelerate CH
3 molecule with Li atom using the moving-magnetic trap decelerator [
29,
30]. This method has two unique advantages. It genuinely has a larger deceleration efficiency over conventional Zeeman deceleration due to larger phase space acceptance for small final velocity, and inherently smoothing deceleration with true three-dimensional trapping potentials, meanwhile it offers simultaneous deceleration of CH
3 and Li since this deceleration scheme does not rely on the ratio of mass to magnetic dipole moment. The simple electronic configuration and light mass of Li atom make study of atom–molecule collisions much less complicated. Thus the decelerated mixture can be used as an ideal test system for sympathetic cooling of CH
3 molecules.
The scheme is shown in Fig. 1, which is based on the moving magnetic trapping technique as demonstrated in Ref. [
29]. The supersonic beam of CH
3 molecules and Li atoms can be created from a pulsed valve of Even–Lavie type, whose temperature can be continuously adjusted through a liquid nitrogen flowing jacket [
27], therefore enabling smoothly velocity tuning of the beam. Using a mixture of 15% CH
4 seeded in Kr, the supersonic beam of CH
3 can be created by DC electric discharge, and have a mean velocity of 330 m/s with a standard deviation of 18% [
24,
27]. Lithium atoms can be entrained into the beam post-nozzle by picking up laser-ablated lithium [
28] or effusive lithium atoms from a heated oven [
31,
32]. As the CH
3–Kr mixture expands over the supersonic expansion, the lithium mixes and assumes the temperature and spatial profile of the CH
3–Kr mixture. After passing through a 5 mm skimmer downstream, the beam enters the decelerator.
In the molecular beam experiments [
24,
27], the typical dimensions of the pulsed methyl radical beam is quite long in beam propagation direction and roughly expand over 4 mm ×4 mm × 3 cm. In order to create a trapping region which can load the central part of the beam, we model our trap from two coils with a inner diameter of 4 mm and a center-to-center distance of 10 mm, which can be wrapped around a thin-walled glass tube used as the deceleration channel of the mixture. The coil geometry is designed as 16 turns (4×4) of a 25 gauge (AWG) Cu wire for the first coil whereas the second has 8 (4×2). The decelerator consists of 198 overlapping quadrupole traps overall thus extending over 997.7 mm.
In Fig. 1 the red arrows represent the current flow, and the current in yellow coils flow in opposite direction to the current in orange coils, thus creating a quadrupole magnetic trap in this anti-Helmholtz configuration. Such a pair of coils are arranged in series along the atomic/molecular beam axis, so the quadrupole magnetic traps created by neighbouring pairs are spatially overlapping with each other. Considering two neighbouring pairs, current pulse of first pair has a half sine shape profile with a pulse width of and the delay time between these two successive pulses is , which means the second current pulse is send through the second pair when the current flowing through the first pair reaches its maximum. Afterwards the current decreases in first pair and increases in the second pair. Thus their current pulse sequences are temporally overlapping, creating a three dimensional moving magnetic trap.
We adapt sinusoidal current pulses as [
29]. They can be described using the form of
, where subscript
n represents the
n-th pair of anti-Helmholtz coils, and modulating frequency
ω sets how fast the trap moves in time. Figure 2 illustrates the resulting spatial magnetic field distribution at
,
,
, and
during the moving of the magnetic trap from first pair to second pair, where
corresponds to the time of activation for the first pair.
The beam will move at constant speed with fixed ω, while the deceleration is accomplished by chirping ω. The deceleration process depends on the frequency chirping rate , which in turn determines current pulse timings and consequent time-varying magnetic fields. The timing of pulse sequence for each pair of anti-Helmholtz coils is calculated according to the initial and final velocity of the particle before and after the decelerator. Assuming a linear frequency chirp of the current pulse, the time-dependent phase of n-th pair can be approximated as , where a is the deceleration value, L is the distance between neighbouring pairs, v0 is the velocity reaching the n-th pair, vi and vf corresponds to the velocity of the beam before and after the whole decelerator, respectively. Therefore, the full chirped current profile can be given as .
3 Deceleration dynamics
In this section, we show the simulation of the one-dimensional simultaneous deceleration of CH3 and Li. The deceleration process in 1D can be approximated by a fictitious time-independent conservative force to a first approximation. It generates a scalar potential:
which tilts the magnetic field potential in the decelerating frame, such that it lowers the front barrier and increases the back barrier, as seen from Fig. 3. The dynamics of the particle inside the trap during deceleration can be approximated by
where is the trap depth, , and ε is a factor proportional to the deceleration value.
In the following, we first show how the magnetic potential is modified according to different acceleration value, based on which we give the Monte-Carlo simulation of the deceleration dynamic process, resulting in the simulated time-of-flight spectrum. Then we obtain the dependence of deceleration efficiency on the final target velocity.
We use a current of 500 A which provides a field magnitude of 1.8 T at the front barrier and 0.9 T magnitude at the back. Figure 3 shows the modified longitudinal field due to deceleration for several deceleration values. In the case of CH3 deceleration with 55.3 km/s2 the fictitious force adds a negative 48 T/m tilt to the 415 T/m of the initial magnetic potential, reducing the height of the front barrier to 1.28 T, which is equivalent to a trap depth of 860 mK.
With these calculation, we model the motion of the particles inside the decelerator. We have calculated their trajectories using Monte-Carlo simulations without any free parameters. In the simulations a million particles are assumed to be Gaussian distributed at the beam origin position where CH
3 are produced by the discharge for both CH
3 and Li. Based on previous experiments [
24,
27], nearly all the produced radicals are populated in the lowest rotational state
and
, which belong to ortho and para type according to nuclear spin statistics, respectively. Assuming the same production efficiency of 7.36% as NH from NH
3 [
35], together with the seeding ratio of 15%, we have 1.1% of CH
3 in the beam. If we use a lower bound 5% [
31] as the entraining efficiency of lithium atoms into the supersonic beam, then the concentration ratio of lithium to CH
3 is about 5:1.
In current simulation, we take the concentration ratio as 1:1 for convenience. Assuming rotational population of the CH
3 radical on the lowest rotational states of each nuclear spin isomers, and population of lithium atom on the lowest six hyperfine states follow Boltzmann distribution, we have performed simulations for six final velocities, 250 m/s, 190 m/s, 140 m/s, 100 m/s, 65 m/s, and 35 m/s. The resulting time-of-flight spectrum from one-dimensional simulation are presented in Fig. 4, which consists of seven time-of-flight traces including the free flight (dashed line) and decelerated packet at various final velocities for both CH
3 (blue line) and Li (red inverted line). We notice CH
3 molecule and Li atom arrive at roughly the same time, which clearly indicate they have been simultaneously decelerated. The observed small difference is due to a longer deceleration time needed to decelerate to the final velocity as pointed out previously by Ref. [
29].
We also obtain the relative number of decelerated CH3 molecules and Li atoms as a function of the final velocity from the simulation, which is shown in left panel of Fig. 5. Both curves were normalized to the final velocity of 255 m/s. For both species the relative number show a monotonic dependence on the final velocity, the smaller the final velocity, the smaller the relative number of decelerated particles. Another feature is larger deceleration efficiency for lithium atoms than for CH3 radicals due to larger magnetic moment to mass ratio for lithium atom. A relative sharp decrease of deceleration efficiency when final velocity is lower than 65 m/s is because of the dependence of the effective magnetic potential on the deceleration value. The effective trapping potentials in the moving frame of reference with a final velocity of 100 m/s for both species are plotted in right panel of Fig. 5. The trap depth of CH3 is more than 0.3 larger is again due to smaller mass to magnetic moment ratio.
4 CH3–Li collision properties
Study collisions between CH
3 and Li would give an important hint for sympathetic cooling of CH
3 molecules by Li atoms. Tscherbul
et al. theoretically studied cold collisions of CH
3 molecules and
3He using both unmodified and strongly anisotropic interaction potentials for He–CH
2 [
36], which gives the ratio of the rate constant for elastic scattering and spin relaxation 9.8×10
12 and 2.8×10
7 at
and
, indicating CH
3 an promising candidate for sympathetic cooling experiments using cold
3He gas. Compared with
3He, lithium atom can be easily laser cooled, and is also an excellent coolant atoms according to recent theoretical calculations [
37–
39]. In a magnetic trap, Tscherbul
et al. have shown the inelastic cross sections for interspecies collisions between
molecular radicals and alkali-metal atoms are strongly suppressed due to the weakness of the spin–rotation interaction in
molecules [
37], and the spin–relaxation collisions would probably be suppressed between spin-stretched CH
3 and Li, thus sympathetic cooling of methyl radical with laser-cooled lithium atoms is likely to be successful.
Here we use quantum diffractive scattering [
40] to model the scattering between co-trapped CH
3 and Li. Assuming the interaction between CH
3 and Li is dominated by long-range van der Waals force, then the interaction potential can be model by an ideal Lennard–Jones potential,
, where the value
C6 can be approximated by
using London dispersion force. Here,
and
are ionization energy of CH
3 and Li, respectively.
and
is the polarizability of CH
3 and Li, respectively. The scattering wavefunction and scattering amplitude can be expanded in terms of the Legendre polynomials
and
respectively, where k is collision wave vector and is the phase shift of the lth partial wave.
The determination of the scattering amplitude and resultant collision cross section
requires finding the partial wave phase shifts, which can be obtained by numerical integration of the radial Schrodinger equation
The solution to the radial equation for each partial wave l is independently computed using the logarithmic-derivative method.
In order to study the diffractive scattering between CH3 and Li, and the prospects for further sympathetic cooling, we assume lithium atoms in the trap are laser cooled to a temperature of 500 μK following magnetic trapping. Figure 6 shows the theoretically computed total cross section for the CH3–Li collisions, which is averaged over a normal velocity distribution at 200 mK. The inset is the partial cross section
as a function of the partial wave value l, which exhibit a universal shape between and and core dependent oscillations below. Figure 7 is a plot of the velocity-averaged loss rate constant versus collision energy for the CH3–Li collisions.
More accurate calculation of scattering cross sections would require not only highly accurate potential energy surfaces constructed by high-level
ab initio electronic structure calculations such as coupled cluster method with single, double, and perturbative triple excitations [CCSD(T)], but also multi-channel scattering calculation [
38], where scattering cross sections between levels
i and
f are given by
5 Conclusion
We have demonstrated the capability of co-deceleration of lithium atom and CH3 molecule using Monte-Carlo simulation of the deceleration process in the moving trap decelerator, and have characterized their deceleration by comparing their deceleration efficiencies, revealing the dependence of the deceleration efficiency on the deceleration value. Our scheme offer several advantages over previous experiments: larger density of decelerated CH3 molecules, and co-trapping of title molecule and atom providing the possibility of study the collision properties of between them, thus open the door for investigating the prospects of sympathetic cooling. Many polyatomic free radicals in the doublet state have similar linear Zeeman effect as CH3 molecule since their spin-rotation interaction is typically smaller than the rotational spacing, thus can be Zeeman decelerated in the same way as we propose here for CH3.
With our ongoing collisional study between lithium atoms and CH
3 molecules, many promising applications will be enabled. For example, after loading them into a magnetic trap, the possibility of creating ultracold CH
3 molecules by sympathetic cooling with ultracold lithium atoms can be stringently tested if elastic collision cross section, inelastic collision cross section, and reactive cross section between them are measured. Bimolecular collisions can also be studied inside such a trap, which has been shown for oxygen molecules [
26,
41], opening an new avenue to investigate the possibility of evaporative cooling. Besides, study of cold reactions between excited lithium atoms and CH
3 molecules are also possible according to previous theoretical calculations [
33,
34]. With the ability of continuously changing the collision energy by tuning the trap depth, we can measure the reaction kinetics between them at the very low temperatures.