Determination of CaOH and CaOCH$_3$ vibrational branching ratios for direct laser cooling and trapping

Alkaline earth monoalkoxide free radicals (MORs) have molecular properties conducive to direct laser cooling to sub-millikelvin temperatures. Using dispersed laser induced fluorescence (DLIF) measurements from a pulsed supersonic molecular beam source we determine vibrational branching ratios and Franck-Condon factors for the MORs CaOH and CaOCH$_{3}$. With narrow linewidth continuous-wave dye laser excitation, we precisely measure fluorescence branching for both $\tilde{X}-\tilde{A}$ and $\tilde{X}-\tilde{B}$ electronic systems in each molecule. Weak symmetry-forbidden decays to excited bending states with non-zero vibrational angular momentum are observed. Normal mode theoretical analysis combined with ab initio structural calculations are performed and compared to experimental results. Our measurements and analysis pave the way for direct laser cooling of these (and other) complex nonlinear polyatomic molecules. We also describe a possible approach to laser cooling and trapping of molecules with fewer symmetries like chiral species.


Introduction
Detailed quantum mechanical understanding of dynamics, interactions, and reactions for complex molecules requires precisely targeted preparation of internal molecular states and relative kinetic motions in the gas-phase solvent-free environment [1][2][3]. Supersonic molecular beams have allowed groundbreaking achievements in experimental studies of inelastic and reactive molecular interactions, as well as opened the path to coherent control of internal molecular states [4][5][6][7][8][9]. With some crucial exceptions, however, large velocity in the laboratory frame usually leads to short experimental coherence times and averaging over multiple internal and motional quantum states. The emergence of cryogenic buffer-gas beams (CBGBs) led to a significant reduction in the forward velocity but the highest molecular fluxes still maintained 10K kinetic energy in the laboratory frame [10]. Therefore, other experimental techniques have emerged in particular those aimed at reducing the forward velocity of molecular beams to the degree sufficient for eventual three-dimensional confinement inside conservative electromagnetic traps, enabling seconds-long coherence times [11]. Molecular trapping also opens opportunities for accurate collisional studies at very low temperatures.
The use of electric, magnetic, optical and mechanical techniques for molecular beam slowing have all been experimentally demonstrated, with published comprehensive reviews providing excellent survey of the field [11,12]. Recent advances in direct molecular cooling as well as coherent molecule formation from unbound atoms have led to the emergence and rapid development of a new vibrant field of cold and ultracold chemistry [3,13]. The unprecedented opportunity to completely control the internal molecular degrees of freedom, as well as external motion, of diverse molecular species has yielded exciting observations on the role of quantum statistics in molecular reactive collisions as well as confining geometry of the reactants, among other results [3].
OPEN ACCESS RECEIVED However, with few crucial exceptions [6,9,14,15], the studies have been limited to the exploration of diatomic molecules which can be prepared from associated laser-cooled atomic samples.
Extension of the available experimental techniques to production of a chemically diverse set of polyatomic molecules will aid in benchmarking state-of-the-art ab initio calculations as well as shed light on novel reaction mechanisms at work. Given the recent rapid experimental progress on laser cooling and trapping of diatomic molecules [16][17][18][19][20], we consider the extension of these techniques to polyatomic species. The internal quantum complexity of polyatomic molecules grows significantly upon transitioning from polar diatomic radicals with a single vibrational mode to nonlinear polyatomic molecules with multiple rotational and vibrational degrees of freedom. Nevertheless, a large number of different monovalent alkaline earth derivatives have been proposed to be amenable to laser cooling [21,22], with pioneering experimental work on the laser-cooling of the triatomic SrOH molecule [23].
In polyatomic molecules, the absence of strict angular momentum selection rules for transitions between totally symmetric vibrational modes presents a new challenge for photon cycling [24], a necessary ingredient for laser cooling. Additionally, higher-order perturbation mechanisms present only in polyatomic molecules such as Fermi resonance as well as Jahn-Teller and Renner-Teller interactions can lead to Born-Oppenheimer breakdown and coupling between vibrational modes, which can result in loss channels otherwise forbidden by the symmetries of the original unperturbed states. Spin-rotation and hyperfine splittings are in the MHz regime and can be efficiently addressed by adding radio-frequency sidebands to the cycling laser; rotational splittings are in the GHz range and can be coupled to the 'bright' rotational state with microwaves; however, vibrational transitions are in the THz range and would require additional laser systems, significantly increasing experimental complexity. Thus, precise measurements of Franck-Condon factors (FCFs) and vibrational branching ratios (VBRs) are crucial in determining the experimental feasibility of laser cooling.
While there is plenty of theoretical work on the estimation of FCFs and VBRs for both diatomic [25][26][27] and polyatomic molecules [21,22,28], accurate experimental studies have been limited primarily to diatomic species [29][30][31][32] and, recently, SrOH [33]. Accurate theoretical predictions of FCFs for polyatomic molecules remain a challenge because their multidimensional nature not only introduces interactions between the degrees of freedom on each potential surface but also between the coordinates of the two electronic states. The pair of molecules we have chosen to study here represent the simplest (by geometric structure) members of the alkaline earth monoalkoxide radicals (MOR) family [34, 35], a large class of polyatomic molecules that has previously been proposed for laser cooling applications [21]. A simple 'triatomic' model of the MOR molecules has been used in order to estimate branching ratios and indicated that laser cooling of relatively large MOR molecules with up to 15 constituent atoms could be possible. Our precise measurements of VBRs for CaOH and CaOCH 3 described below provide the first experimental validation of the theoretical model of 'shielding' of the R ligand group (e.g. CH 3 ) vibrational modes by the intermediate oxygen atom [21]. Increased density of vibrational states for larger polyatomic molecules could potentially lead to internal vibrational redistribution among multiple vibrational normal modes leading to enhanced loss probabilities and, therefore, necessitating accurate experimental studies.
CaOH is one of the first polyatomic molecules to be proposed suitable for laser cooling and trapping and detailed theoretical ab initio calculations have been performed of the FCFs [22] as well as energy levels [36]. Our work provides an important comparison between the theoretical and experimental results, benchmarking theoretical calculations by Isaev and Berger [22]. Our precise measurements of the branching ratios for a more complex symmetric-top molecule (STM) CaOCH 3 provide, to the best of our knowledge, the first accurate experimental study of the VBRs for a nonlinear radical suitable for direct laser cooling. Particularly, our results indicate the necessity to consider pseudo Jahn-Teller couplings in the excited electronic levels of polyatomic molecules [37]. Perhaps surprisingly, despite the presence of Jahn-Teller interaction in the electronically degenerate E symmetry state of the C 3v symmetric top, using the lowest allowed X A A E «˜electronic excitation will lead to ∼16 scattered photons without any additional repumping lasers. Based on our experimental measurements and theoretical analysis, we propose a feasible route for producing large ultracold samples (∼10 3 -10 6 molecules) of both radicals via laser cooling on the lower spin-orbit (SO) branch of the X Ã -ẽ lectronic transition. Performing a detailed study of inelastic and reactive collisions on increasingly complex members of the CaOR molecule family would provide a unique window into the scattering properties of fundamental ligand groups in organic chemistry.

Experimental configuration
Since the details of the experimental apparatus have previously been described in other publications [33], we provide only a brief account of the experimental configuration employed. As a source of cold molecules we used a pulsed supersonic molecular beam with argon as a carrier gas. Both CaOH and CaOCH 3 beams were created by flowing argon over a container filled with liquid methanol (CH 3 OH) maintained at room temperature. Room temperature vapor pressure of methanol (∼10 kPa) together with the backing pressure of argon (∼4000 kPa) was enough to seed a sufficient number of gas-phase molecules into the carrier gas to observe large, stable yields of both molecular radicals of interest following laser ablation of a calcium metal target. The peak fluorescence signal per single rovibrational level was approximately of similar magnitude for CaOH and CaOCH 3 , which allowed for accurate studies of both radicals using the same beam source. Molecules in a specific rotational state were excited using a single frequency cw-dye laser. Dispersed fluorescence data was collected with a spectrometer that included 0.67 m focal length, high-efficiency Czerny-Turner-type monochromator with a low-dispersion grating and a cooled, gated intensified charge-coupled device camera. The relative sensitivity as a function of wavelength for the spectrometer was precisely calibrated beforehand. The wavelength resolution of the emission spectrum was controlled by adjusting the width of the monochromator input aperture and absolute calibration was performed with the argon emission lamp. The measured emission wavelengths agreed to 0.1 nm with theoretical predictions which allowed for unambiguous identification of the vibrational loss channels. The narrow linewidth (1 MHz) of the excitation dye laser allowed accurate rotationally-resolved excitation. A Burleigh high-resolution wavemeter was used as a rough frequency reference (<1 GHz), additionally confirmed by co-recording with some Doppler I 2 spectrum [38]. For each dispersed laser induced fluorescence (DLIF) measurement, three separate datasets were taken: (i) signal dataset with ablation Nd:YAG laser on and dye laser on, (ii) scattered background dataset with no Nd:YAG laser but dye laser on, and (iii) spurious ablation glow dataset with Nd:YAG laser present but no dye laser. The resulting measurement was obtained by subtracting datasets (ii) and (iii) from (i) in order to remove the scattered light offset and reduce the background ablation glow (e.g. metastable calcium emission).

Experimental determination of branching ratios
In order to eliminate rotational branching during the photon cycling process the use of J J 1   ¢type angular momentum transitions has been proposed [39,40] and for the first time demonstrated with diatomic SrF in [41] and linear triatomic SrOH in [42]. Rotationally-resolved high-resolution spectroscopy of low-J rotational levels of CaOH and CaOCH 3 has been previously performed and assigned for both X Ã -˜and X B -ẽ lectronic bands [43][44][45][46][47][48][49]. The availability of prior spectroscopic measurements, low rotational temperatures of the skimmed supersonic molecular beam as well as the narrow linewidth of the cw dye laser, allowed us to deliberately address only the rotational transitions which can be used for optical cycling and laser cooling. Additionally, the absence of state-changing collisions in the probing region eliminated any possible systematic errors.
While both CaOH and CaOCH 3 have spin-rotation and hyperfine splittings arising from the unpaired electron and hydrogen spins, correspondingly, for the purposes of the experiments performed here we label the states using rotational quantum numbers N K , ñ | for non-degenerate quantum states (i.e. Σ + states for linear molecules and A 1 states for STMs) and J K , ñ | states for degenerate states (i.e. Π electronic states for linear molecules and E states for STMs) with K=0 for linear molecules. Upon the electronic excitation to a specific rotational level, the spontaneous emission rate for a dipole-allowed transition is governed by the Einstein A coefficient [50]: under the Born-Oppenheimer separation tot el vib rot y y y y = . Therefore, the intensity of different vibrational bands will be proportional to since the same electronic transition dipole moment R e | | and Hönl-London rotational factor are shared by all the vibronic emission bands. Additionally, the rotational branching is limited by the J 0, 1  =  selection rule and all J branches contribute identically for all the vibrational bands. For polyatomic molecules, the FCF is multidimensional and given as By comparing the integrated areas under the dispersed fluorescence emission peaks, we can determine the relative branching ratios for decays from a given electronic, vibrational and rotational state.

CaOH measurements
Two electronic transitions have been studied for CaOH that have potential to be used in the laser cooling and trapping applications. Previously, laser cooling of the isoelectronic CaF and SrOH molecules has been demonstrated using either the X A -˜electronic transition and, therefore, we have performed detailed measurements on both bands for CaOH as well. Additionally, in order to increase the photon cycling rate it is favorable to avoid coupling multiple excitation lasers to the same vibronic level, thus requiring separation of the main cycling and repumping lasers to address different electronic levels (e.g. main cycling on X Ã -˜while repumping on X B -˜or vise versa). As was previously demonstrated both theoretically [51] and experimentally [52], when multiple ground state sub-levels n vg from different vibrational levels v are coupled to the same excited state sub-levels n 0e , the effective scattering rate for molecules is reduced by a factor of n n n 2 e e v vg 0 0 + å ( ) , resulting in a lower radiation pressure force (F rad ). Therefore, the scattering rate and F rad is maximized for decoupling all repumpers from the main cycling laser, resulting in a maximal possible scattering rate of R n n n e e g scat,max sp 0 0 0 ) vibrational state is symmetry allowed and observed in our measurements ( figure 1(b)), but the more precise determination of this band's intensity is limited by the contamination signal arising due to the 657 nm emission from the metastable atomic calcium created during the laser ablation process. Additionally, it is important to note that for both A X ˜and B X ˜emission bands, decay to vibrational levels with v2 is significantly suppressed to below 0.3%. This rapid decrease in the intensity of offdiagonal vibrational decays should enable experimentally accessible laser cooling and trapping.

CaOCH 3 measurements
The presence of K 0 ¹ states for nonlinear molecules requires additional considerations for achieving effective photon cycling beyond those examined previously for diatomic and linear polyatomic molecules. Reference [21] outlined how to use K K 1 0  =  ¢ = transition for the perpendicular X Ã -˜transition and K K 0 0  =  ¢ = transition for the parallel X B -˜to achieve rotationally closed excitations. Figures 1 and 3 provide a comparison of the DLIF for CaOH and CaOCH 3 . Both of the molecules were excited on the X Ã -ẽ lectronic transitions around 630nm and in a rotationally resolved manner to J 1 2 ¢ = state. As can be clearly seen from the data presented, both CaOH and CaOCH 3 have a very small number of vibrational decay channels despite significantly increased structural complexity in going from a triatomic to a hexatomic molecule. While  CaOCH 3 has eight distinct vibrational modes (4 of a 1 symmetry and 4 degenerate pairs of e symmetry), only two are optically active at the three part per thousand level as seen in the figure. Our measurements, performed on a cold molecular beam with a narrow-band dye laser exciting only the rotational level suitable for laser cooling, confirm that even nonlinear CaOR molecules behave effectively like triatomics, showing great promise for laser cooling and trapping.
In order to account for possible systematic effects associated with multiple-photon cycling in our measurements, we have also used excitation to the upper SO component of the Ã state ( 2 3 2 P for CaOH and E 2 3 2 for CaOCH 3 ) which does not support rotationally-closed excitations. For all measurements performed, we did not add any additional laser sidebands, thus addressing only a single spin-rotation component. In the presence of optical saturation of the transitions where more than one photon is scattered per single molecule, we anticipate to see 3 5 » less signal in the 0 0 0 emission band when using a non-rotationally closed transition (e.g. R or Q branch). However, for both molecules the measured branching ratios for excitation to either the lower or the upper SO component of the Ã state are essentially consistent with each other (see tables 5 and 6), providing strong evidence that the measured branch intensities are dominated by the actual FCF values and setting an upper limit on the associated systematic error of 1%. In order to confirm our measurements for the X B t ransitions we used R-type branch to excite molecules in rotational levels N″>1, thus preventing optical cycling. Additional data for such measurements is provided in the appendix. We have also used dispersed LIF to explore the possibility of using the X B -˜electronic transition in CaOCH 3 around 566 nm for laser cooling. Figure 4 shows the DLIF results indicating two additional features not present in the A X ˜emission spectrum: (i) emission to 1 quanta of the Ca-O-C bending mode and (ii) a new feature around 1750 cm −1 away from the excitation band. Since B is the second excited electronic level, there is a possibility of mixing with the excited vibrational levels of the Ã state. The proposed intensity borrowing mechanism is shown in figure 5, indicating the resulting decays nominally forbidden by the symmetry arguments. As clearly seen in the data presented in figure 4, a relatively strong emission to the band offset by 1751 cm −1 from the diagonal vibronic transition is observed in the B X ˜decay of CaOCH 3 , making laser cooling on this electronic transition more challenging compared to X A E 2 1 2

˜.
Any nonlinear molecule (e.g. CaOCH 3 ) in an orbitally degenerate electronic state (e.g. E 2 state) will always distort in a way such as to lower the symmetry and remove degeneracy. This so-called Jahn-Teller effect (JTE) nullifies the v 2, 4, i  =   ¼selection rule for non-symmetric vibrations in electronic transitions [50]. JTE in the Ã electronic state as well as near-degeneracy between the excited combination band in the Ã state and the ground vibrational level of the B electronic state lead to Born-Oppenheimer approximation breakdown.
is responsible for mixing of the zeroth-order Born-Oppenheimer vibrational and electronic wavefunctions. Careful theoretical understanding of vibronic coupling in the excited B state is currently in progress and beyond the scope of the present paper.

Comparison with theoretical predictions
CaOH three-dimensional potential energy surfaces have been previously mapped out using ab initio calculations for both the ground [36] as well as excited electronic states [57]. While table 1 provides a comparison between such ab initio results from different studies, we further focus our comparison on the calculations performed by Isaev and Berger [22] since, to our knowledge, they provide the only previous study of CaOH FCFs relevant to laser cooling prospects. Even though there are no previous ab initio calculations of FCFs for CaOCH 3 , the use of isotopically substituted methoxy groups has allowed accurate experimental determination of the radical's geometry in both X and Ã states [45]. Because of the highly diagonal FCF matrices for both CaOH and CaOCH 3 molecules, only a few excited vibrational modes were present in our DLIF data as seen in figures 1-4. Other than the unexpected peak seen in the 1751 nm overtone band for the B X ˜emission, the frequencies of other vibrational modes have been previously measured and tabulated which allowed for unambiguous vibrational character assignment. Theoretical analysis of the multidimensional FCFs for all vibrational modes-including those not observed in the experiment-depends on the frequencies of the corresponding vibrational motions. We therefore initialized  our calculations with ab initio molecular geometries and vibrational frequencies using the ORCA quantum chemistry program, which is discussed in [59]. Tables 1-3 provide a summary of the calculated normal vibrational modes for CaOH and CaOCH 3 . We observe excellent agreement between our ORCA calculations with previous ab initio calculations for CaOH by Isaev and Berger [22] as well as previously experimentally determined values for the vibrational normal modes. In the case of MOCH 3 molecules, there are 12 normal modes, 4 of a 1 symmetry and 4 degenerate pairs of e 1 symmetry. Thus, while there are 12 internal coordinates, we reduced the problem to 8 symmetry coordinates listed in table 4. The details of the numerical computation of CaOH and CaOCH 3 FCFs are described in the appendix. The VBRs were computed by weighting each FCF by the cube of the transition frequency from the vibrationless excited electronic state to the vibrational state in the ground electronic state and normalizing. Results of these calculations are shown in tables 5 and 6. We obtain excellent agreement between the measured and calculated VBRs for CaOH and a reasonable agreement for a more complex CaOCH 3 molecule. Incorporating possible couplings between different vibrational modes as well as inclusion of the anharmonic terms in the vibrational potential could resolve the slight discrepancy between the measured and predicted VBRs for CaOCH 3 . Notice that our theoretical predictions for FCFs are much closer to the measured values we have observed compared to the ab initio values predicted by Isaev and Berger [22]. By using available experimental vibrational frequencies to determine the force constants as well as incorporating experimentally measured bond lengths and structure into our analysis, we were able to circumvent some of the challenges associated with predictions of FCFs with

Laser cooling and trapping prospects
Our experimental measurements and theoretical calculations confirm the possibility of achieving direct laser cooling for CaOH and CaOCH 3 with a few repumping lasers. Furthermore, our measurements provide the quantitative guidance to choosing the most experimentally efficient route to identifying the main cooling and all the repumping transition bands. Using the measured values for the branching ratios, we have determined the optimal laser cooling schemes for CaOH and CaOCH 3 depicted in figures 6 and 7, respectively. Because of the relatively light mass, CaOH can be slowed down to a complete stop from 100 m s −1 CBGB using only 9000 photons which can be obtained with three repumping lasers indicated. Thus, the technical complexity of laser slowing CaOH is comparable to experiments with diatomic molecules SrF and CaF. While additional structural complexity for CaOCH 3 results in less diagonal FCFs, with the repumping scheme in figure 7 and assuming ∼10 9 molecules in a single rovibrational quantum state at 100 m s −1 , one can obtain between 10 4 and 10 7 slowed molecules in the trapping region, assuming 1 part per 1000 or 1 part per 3000 loss probability to unaddressed vibrational levels.

Beyond symmetric top molecules
Because of the strongly ionic nature of the Ca-O bond in CaOR molecules, the electronic transitions involve non-bonding Ca + orbitals perturbed by an RO − ligand even for asymmetric rotor molecules (ARMs). Since the Ca-O-C bond angle remains 180 • , the local symmetry near the calcium metal is linear and therefore many characteristics of electronic transitions carry over from linear CaOH to more complex CaOR symmetric top and asymmetric rotor radicals [35,63]. For example, there is 'SO' splitting of about the same size, and the electronic orbital angular momentum remains essentially unquenched. However, the symmetry reduction in going from C v ¥ (CaOH) and C 3v (CaOCH 3 ) to molecules with lower symmetry properties like C S and no symmetry has potential to impact the VBRs in a significant manner. For linear (C v ¥ ) and symmetric top (C 3v ) molecules, Ca-O-C bending vibrations have π or e symmetry, correspondingly, and therefore only v 2, 4,  =   ¼ branches appear in the electronic decays under the Born-Oppenheimer approximation 6 . Such a bending mode becomes non-degenerate for molecules with lower symmetries than C 3v , therefore leading to allowed decays of comparable intensity to both one quanta in the Ca-O stretching and Ca-O-C bending modes as has been previously observed for calcium monoalkoxides [35]. While this presents a technical challenge, requiring an additional laser for repumping lost molecules from such a decay channel, it does not pose a fundamental threat to achieving photon cycling in ARMs like CaOC 2 H 5 or CaOCH 2 CH 3 .
Depending on the orientation of the transition dipole moment relative to the principal axes of the molecule, vibrational and electronic spectra of ARMs are characterized as a-type, b-type or c-type transitions with the selection rules [50]:   Therefore, as can be seen from the selection rules above, for near-oblate or near-prolate ARMs, the overall structure of the emission bands resembles that of parallel or perpendicular transitions for STMs. Choosing either a-type or ctype transitions and driving either K a a  =  ¢ = , correspondingly will resemble optical cycling transitions for STMs as described above. Additionally, using quantum states such as K 0 a  = or K 0 c  = for a-type and c-type transitions, correspondingly, will restrict J selection rules to J 1  =  only (similarly to N 1  =  transitions used for laser cooling diatomic and linear molecules). Thus, while reduced symmetry for ARMs leads to novel challenges in repumping vibrational levels, carefully chosen rotational transitions can significantly limit (or even eliminate) rotational branching. As a first step in achieving photon cycling with ARMs, it is necessary to perform high resolution spectroscopy of the relevant electronic transitions and identify suitable bands for achieving optical cycling. While we have already completed such work for CaOCH 3 as described above, similar spectroscopy on more complex CaOR radicals (e.g. CaOCH 2 CH 3 ) can be performed in the future. A natural chiral candidate for laser cooling is a chiral analog of calcium monomethoxide, CaOCHDT. A recent theoretical and experimental work by Liu and coworkers [66] on calcium ethoxide (CaOC 2 H 5 ) confirms that Ca-O stretching mode is the dominant vibrational loss channel for this molecule and indicates its potential suitability for optical cycling and laser cooling using a scheme outlined here for ARMs.

Conclusions and future directions
While the high density of vibronic states for complex polyatomic molecules can potentially inhibit photon cycling using the X B -˜electronic transitions, our measurements provide a strong case that using the lowest excited electronic level (Ã) will enable photon cycling and laser cooling despite the presence of SO and Jahn-Teller interactions. In addition to precise measurements of the VBRs for both X Ã -˜and X B -˜electronic excitations, we have also performed theoretical calculations that agree quite well for the X Ã -˜band, where harmonic oscillator approximation remains valid for calculating the multidimensional FCFs. Using the analytical integral expressions presented by Sharp and Rosenstock [67], we have calculated the relevant FCFs using previously available experimental data for geometries, vibrational frequencies, and corresponding force fields. The presence of strong optically accessible electronic transitions, as well as ability to scatter multiple photons from a single molecule, opens the possibility for internal state molecular manipulation and control of the molecular motion in the laboratory frame. While radiative force slowing and cooling requires multiple thousands of repeated absorptionemission cycles [68], Sisyphus cooling processes as well stimulated emission slowing methods can provide significant gains even with a limited number of optical cycles [69,70]. Stimulated optical forces arising from polychromatic light beams [71,72] might be especially effective for manipulating MOR molecules with reduced symmetries and increased complexity where scattering tens of thousands of photons would require a large number of vibrational repumping lasers. Additionally, rapid optical cycling of even tens of photons can significantly improve the rate of opto-electrical cooling methods for STMs and ARMs, which have so far relied on using vibrational transitions with long lifetimes and therefore low scattering rates [73,74].
Laser cooled polyatomic molecules provide an ideal starting point for producing ultracold fundamental radicals of chemical interest like OH, CH 3 and OCH 3 . While a dedicated theoretical exploration is necessary for identifying the exact route for efficient photodissociation of polyatomic molecules into underlying constituents with minimal energy released, this approach seems feasible with the use of the STIRAP, which is the reverse of the process employed for forming ultracold diatomic molecules from ultracold alkali atoms [75]. Release of as little as 100 nK of additional kinetic energy into the motion of dissociation products has been achieved with ultracold Sr 2 molecules [76], and an efficient pathway for SrOH zero-kinetic-energy dissociation is currently being developed using ab initio molecular potentials [77]. Moreover, modern theoretical approaches to molecular quantum scattering calculations are reaching the regime of larger polyatomic molecules [78] and, therefore, experimental results on collisions would play a crucial role in further developing the field.
Furthermore, our measurements indicate the possibility of multiple photon cycling for other polyatomic species including polyatomic molecular ions. Internal state cooling and control of polyatomic molecular ions can significantly benefit from laser-induced optical cycling. Polyatomic molecular ions ScOH + , YOH + , as well as ScOCH 3 + and YOCH 3 + , are isoelectronic to the corresponding calcium and strontium compounds [79,80] and should be highly suitable for direct optical manipulation and internal state cooling using methods previously demonstrated with diatomic ions [81]. A few particularly interesting candidates to consider are AcOH + and AcOCH 3 + , which should be similar in electronic structure to RaOH and RaOCH 3 and are sensitive to new physics beyond the standard model [22], but with a significantly longer half-life of the actinium-227 nucleus compared to radium-227 [82].
Another interesting class of molecules to explore is mixed hypermettalic neutral and ionic oxides of the MOM′ type, where M and M′ refer to two different metal elements. Recently, quantum state controlled synthesis of BaOCa + has been demonstrated [83] and theoretical work on neutral MOM (same metal) molecules has been performed motivated by the potential to search for fine structure constant α variation [84]. Neutral polyatomic molecules functionalized with optical cycling centers like Ca or Sr could also enable interesting applications in quantum sciences [85,86].

A.2. Details of FCF calculations
To compute the FCFs for CaOH and CaOCH 3 , we used Wilson's GF matrix method and the Sharp-Rosenstock expansion of the FC overlap integral using a harmonic potential. In depth discussion and derivation of these methods are extensively addressed elsewhere [67,89,90]. Importantly, the Sharp-Rosenstock expansion depends on the vibrational coordinates of the initial and final states, denoted Q¢ and Q respectively. The relationship between these two coordinate systems is: where K is the difference between the equilibrium geometry of the two states in terms of the initial vibrational coordinates, while J accounts for coordinate mixing, also known as the Duchinsky effect. Wilson's GF matrix method [91] is the extension of eigenvalue vibration problems to molecular vibrations, and can be used to find the relationship between Q and Q¢. For  ]. The elements of the G (geometry) matrix are derived from mass-weighting the summed scalar product of each of the 3N−6 internal coordinates over all atoms, while the F (force) matrix elements are derived using a harmonic approximation-each element of the F matrix is the second derivative of the potential energy surface with respect to the internal coordinates.
From the G matrix, it is possible to derive the kinetic energy, and from the F matrix it is possible to derive the potential energy.
Writing the Lagrangian of the system, we end up with a determinant problem: and ν is the frequency of a molecular vibration. Multiplying by the determinant of G, we get a diagonalization problem: where the corresponding eigenvectors are the vibrations (normal modes) in terms of the internal coordinates. Thus, the linear transformation L formed by the eigenvectors maps from vibrational normal mode coordinates to internal coordinates. By repeating the GF analysis for the initial and final electronic states, we end up with L matrices for both states, and thus can define Q and Q¢.
Specifically, J is defined as the product of the inverse of the initial state L matrix and and the final state L, while K is the transformation of the equilibrium difference in internal coordinates to normal coordinates Due to vibrational degeneracies, it will often be preferable to reduce the dimensions of the GF problem by concatenating the internal coordinates into symmetry coordinates shown in table 4. To summarize, these Franck-Condon calculations required knowledge of the molecule geometry, vibrational frequencies, and the potential energy surface in the form of harmonic force constants defined over an internal coordinate system. Suitable internal and symmetry coordinate systems were analytically derived using literature on methanol (HOCH 3 ) [92], and the force constants were fit using PGOPHER software to molecular geometry and spectroscopic data [93].
The force constant fitting was initialized by first reducing the dimensionality of the GF problem to a triatomic situation by treating the CH 3 methyl group as a single atom. Given an ordering of the eigenvalues and neglecting off-diagonal F matrix interactions, the 3×3 GF eigenvalue problem is completely determined without knowledge of the diagonal matrix terms. We determine the correct ordering by comparing calculated FC values with known measurements and prior theory, initializing the PGOPHER force constant fitting using force values for the M-O stretch, O-CH 3 stretch, and M-O-C bend symmetry coordinates. To prevent misidentification of normal modes with frequencies, we separately fit the force constants for vibrations of different symmetries, which should be completely decoupled in such a way that the L matrix is block-diagonal. We then computed the FC overlap integrals using the Sharp-Rosenstock expansion with the matrices J, K, and G, a diagonal matrix consisting of i reduced vibrational frequencies h 4 i 2 p n , for both the initial and final electronic states. The FCFs were derived by calculating and normalizing the overlap integral for each vibration up to the third quantum number.

A.3. Additional supporting measurements
In this Section we provide additional supporting data in figures 8-10. 7 Three degrees of freedom are taken up for position, three for rotation. Accordingly, a linear molecule has 3N−5 degrees of freedom due to rotational invariance around the primary symmetry axis. -ẽ lectronic transition. By using R-branch, we avoid potential systematic error due to accidental optical cycling. As summarized in tables 5 and 6, our measurements are consistent with theoretical predictions and values obtained by exciting the lower spin-orbit Figure 9. Comparison of DLIF results (blue curves) for CaOH (a) and CaOCH 3 (b) excited to the lower spin-orbit component of the X Ã -˜electronic transition with laser scatter background (red curves) provided for reference. While we have already subtracted laser scatter background from the DLIF data, for both CaOH and CaOCH 3 datasets it constituted only a small fraction of the observed signal near the 0 0 0 emission band (as can be seen from the plots here). Figure 10. Comparison of DLIF results for CaOH (a) and CaOCH 3 (b) excited using R-type branches of the X B -˜electronic transitions. By using a R-branch for N>1 we avoid potential systematic error due to accidental optical cycling, confirming our measurements with the rotationally-closed excitations described in the main text. Observed vibrational decay channels are indicated.