\section{Introduction}\label{sec1}
Many published papers proposing new m-valued circuits justified their proposal by interconnection issues in binary circuits. Just to present one among a lot of quotes, "One of the main problems in binary logic is the high volume of interconnections which can increase the chip area and power consumption" \cite{Sha}. It is quite evident that m-valued signals carry more information than binary signals and that less m-valued signals would be needed than binary ones to carry the same amount of information. However, if the objective is to reduce the interconnect issues of binary circuits, 4-valued circuits are the best solution as they allow a simple interface between binary and 4-valued circuits.
There are two different options:
\begin{itemize}
\item Using binary circuits and using 4-valued signaling between binary circuits, as shown in Figure \ref{4Vint}. In that case, the number of interconnects between the binary parts is divided by 2.
\item Using 4-valued circuits in a binary environment, as shown in Figure \ref{4VC}. In that case, the 4-valued circuits are based on 4-valued operators with  4-valued interconnects.

\end{itemize}

\begin{figure}[htbp]
\centerline{\includegraphics  [width =7 cm]{4Vint.png}}
\caption{4-valued interconnects between binary circuits}
\label{4Vint}
\end{figure}
\begin{figure}[htbp]
\centerline{\includegraphics  [width =4 cm]{4VC.png}}
\caption{4-valued  circuits in a binary environment}
\label{4VC}
\end{figure}

\section{Methodology}
In the last decade, most proposed designs of m-valued circuits use simulations of CNTFET technologies. The CNTFET circuitry is similar to the CMOS one while presenting a significant advantage: the threshold voltage varies as the inverse function of the diameter of the tube. Different diameters provide the different threshold voltages that are needed in m-valued circuits.

 However, this technology  is far from being mature. In 2019, a 16-bit RISC microprocessor has been built with 14,000 CNFET transistors \cite{Hills}. While this is a significant advance for CNTFET technology, we may observe that the Intel 8086 CPU, which was a 16-bit microprocessor, has been launched in 1978 with 29,000 transistors, more than 40 years ago! As of 2019, the largest  transistor count in a commercially available microprocessor was 39.54 billion MOSFETs, in AMD's Zen 2 fabricated using TSMC's 7 nm FinFET semiconductor manufacturing process. As of 2020, the largest transistor count in a GPU (NVidia Ampere) was 54 billion transistors with the same 7 nm process. The CNTFET 16-bit microprocessor manufacturing process had 5 metal layers, while the number of metal layers in nano-CMOS technologies ranges from 8 to 15, with a trade-off between integration and cost. So, comparing m-valued circuits using CNTFET technologies with the most recent FinFET or FDSOI  technologies would be a nonsense.

As the point is to compare the interconnections between binary and 4-valued circuits, we use the CNTFET technology as a basis for all the comparisons. Directly comparing interconnects is impossible. While most 4-valued proposals present the electrical scheme, they generally neither provide any layout information nor even the electrical scheme of the corresponding binary circuits. While interconnects are represented on a plane in the electrical scheme, they are implemented according to several layers in the IC layout.  Figure \ref{Layout} presents the layout and the circuit scheme of a 2-input multiplexer \cite{Huang}. All the connections are not shown in the electrical scheme. The connections are implemented as Metal 2 (blue), Metal 1 (violet) and polysilicon (pink).  However, considering the connections at the transistor level gives insight on the interconnection at the layout levels. There are few possibilities that more connections at the transistor levels lead to less interconnections and less chip area at the layout level. If there are significant differences between number of connections at transistor level between m-valued and binary circuits, we can easily assume that the situation is similar at the layout level. It is also very doubtful that many more transistors could lead to less interconnects. The ratio m-valued transistor count/binary transistor count therefore gives some indication.

\begin{figure}[htbp]
\centerline{\includegraphics  [width =8 cm]{2-1MUX.png}}
\caption{2:1 CNTFET multiplexer. Layout (left) and circuit (right)}
\label{Layout}
\end{figure}

\section{4-valued interconnects between binary circuits}
\subsection{Telecom applications}
There is one significant application of this approach. Reducing the number of
interconnects with multiple levels is used in amplitude
modulation: for instance, PAM-4 coding \cite{PAM4} (Figure \ref{PAM4}) that uses 4
levels to code 2 bits is adopted for high-speed data
transmission (IEEE802.3bs).

\begin{figure}[htbp]
\centerline{\includegraphics  [width = 8 cm]{PAM4.png}}
\caption{PAM4  encoding}
\label{PAM4}
\end{figure}

However, telecom industry and computer industry are quite different.
\subsection{Digital electronics and computer industry}
\subsubsection{Dividing by 2 the number of interconnects}
Different papers have been published in which 4-valued signaling is used
between binary circuits in different technologies: TTL \cite{etiTTL}, ECL  \cite
{etiECL}, CMOS, CNTFET, etc.  2-valued to 4-valued encoder and 4-valued to 2-valued
decoder circuits must be designed for each technology.

With CNTFET technology, a possible implementation is now presented. The truth table of the decoder using Gray code is presented in Table \ref{NQI}. The threshold detector circuits implementing functions NQI, IQI and PQI are shown in Figure \ref{Tdetectors}. The corresponding 4 to 2 decoder with 14 T is presented in Figure \ref {QBGray}. A possible implementation of the 2 to 4 encoder is presented in Figure \ref{BQGray}. It has 12 T.

\begin{table}
\centering
\caption{Gray 4 to 2 decoder}
\begin{tabular}{|c||c|c|c||c|c||}
  \hline
 Q&NQI&IQI&PQI&X&Y\\
\hline
 0&3&3&3&1&0\\
 1&0&3&3&1&1\\
2&0&0&3&0&1\\
3&0&0&0&0&0\\
  \hline
\end{tabular}
\label {NQI}
\end{table}

\begin{figure}[htbp]
\centerline{\includegraphics  [width =8 cm]{4to2decoder.png}}
\caption{Threshold detector circuits presented in \cite{Ebrahimi}}
\label{Tdetectors}
\end{figure}

\begin{figure}[htbp]
\centerline{\includegraphics  [width =6 cm]{QBGray.png}}
\caption{Quaternary to binary decoder (Gray code)}
\label{QBGray}
\end{figure}

\begin{figure}[htbp]
\centerline{\includegraphics  [width =4 cm]{BQGray.png}}
\caption{Binary to quaternary encoder (Gray code)}
\label{BQGray}
\end{figure}

\subsubsection{Doubling the number of interconnects}
 It turns out that 4-valued interconnects have not been used by IC designers.
For "long distance" or "high-speed" transmissions, the solution has been to use differential signals instead of single end signals, as shown in Fig. \ref{Dif}. Differential signaling has advantages over single-ended such as:
\begin{itemize}
\item Better noise margins, as the signal swing is the difference between positive and negative signals
\item Self-reference as the threshold is $(positive signal + negative signal)/2$
\item Reduced switching time due to lower voltage swing.
\item Subtraction between positive and negative signals rejects common-mode noise.
\end{itemize}

\begin{figure}[htbp]
\centerline{\includegraphics  [width = 7 cm]{DIFSIG.png}}
\caption{Differential signaling}
\label{Dif}
\end{figure}

While differential links are used at physical level, different buses and links have been defined and used in computer architectures. There are many examples: SuperSpeed USB  3.0 \cite{USB}, PCIe \cite{PCI}, XAUI \cite{XAUI}], InfiniBand \cite{InfiniBand}, RapidIO \cite{RapidIO} , and SATA  \cite{SATA}. We only give some more details for two examples:
\begin{itemize}
\item PCI Express: It is a high-speed serial bus. A PCI express link (or interconnect) between two devices consists in a number of lanes, each lane being composed of two differential signaling pairs, one for receiving data and the other one for transmitting. Successive versions have been defined. PCI Express 1.0 has been introduced in 2003 and PCI Express 5.0 is expected in 2022.
\item NVLink: To improve the interconnection bandwidth between CPU and GPUs, NVidia introduced a new interconnect architecture called NVLink \cite{NVlink}. A single NVLink is a bidirectional interface incorporating 32 wires forming eight differential pairs in each direction. The first implementation NVLink1 was introduced in the Tesla P100 GPU (2016). Each link has a 40 GB/s bidirectional bandwidth. NVLink1 is also used by IBM in the Power8 microarchitecture. NVlink2 was introduced with GV100 GPU (2018). Bidirectional bandwidth is now 50 GB/s.
\end{itemize}

\section{4-valued interconnects within 4-valued circuits}
In this section, we examine whether 4-valued circuits reduce the number of interconnects compared to the corresponding binary circuits. According to Shannon theory of information, when N events have the same probability to occur, the corresponding amount of information is $ I=log_2(N)$ bits (or shannons). When N=2, I=1 bit and when N=4, I=2 bits. A quaternary wire carries two times the amount of information of a binary one. It means that the information ratio is 2. A quaternary circuit will be more efficient than a binary one if its complexity is not greater than 2.
Comparing quaternary circuits and the corresponding binary ones is not easy. Which technology do we consider? How to evaluate complexity?

\subsection{Quaternary logic gates}
Inverters are the first gate to consider. Figure \ref{4VINV} shows a 4-valued inverter presented in \cite {Sharifi}. This scheme allows a quick comparison with binary circuits. A 4-valued inverter carries the same amount of information than two binary inverters: 2 bits. Right part of Figure 2 includes two binary inverters, which means that the six supplementary transistors are just the overhead of the 4-valued inverter compared to two binary inverters. All the connections between these transistors are overhead, which doesn't exist for binary implementation. It should also be noticed that the 4-valued implementation uses three power supply rails ($V_{dd}, 2V_{dd}/3, V_{dd}/3$) instead of one. Power supply distribution at the layout level is much more complicated. The quaternary inverter presented in \cite{Ebrahimi} has also 10 transistors.

\begin{figure}[htbp]
\centerline{\includegraphics  [width = 5 cm]{4VINV.png}}
\caption{A 4-valued CNTFET inverter}
\label{4VINV}
\end{figure}

Obviously, the situation is similar for NOR, NAND and more complex gates. Table \ref{QNAND-TC} compare the transistor count of two different 4-valued NAND with the  equivalent binary NAND. For processing 2 bits of information, the quaternary Nands have at least two times more transistors. The results are the same for NOR gates.
\begin{table}
\centering
\caption{2-input Nand gate transistor count}
\label{QNAND-TC}
\begin{tabular}{|c|c|c|c|c|c|}
\hline
2-input Nand  & 4-V Nand \cite{Sharifi}& 4-V Nand\cite{Ebrahimi}&Binary Nand  \\
\hline
Transistor count  &20&16&4 \\
\hline
4V/2V TC ratio &5&4&1\\
\hline
\end{tabular}
\end{table}
\subsection{Quaternary full adders} \label{ENCDEC}
A detailed examination of the implementation of quaternary adders has been presented in \cite{etiQADD}. They are based on the following versions:
\begin{itemize}
\item V1: Interfacing a 2-bit binary adder with a 4-to-2 decoder and a 2-to-4 encoder.
\item V2: Direct implementation of Table \ref{T1}.  A function f(inputs) is decompose into f(inputs) = 3.f3 + 2.f2 + 1.f1 where f3, f2 and f1 are respectively the binary functions of the inputs for which the functions have values 3, 2 and 1. f3, f2 and f1 includes the NQI, IQI and PQI functions of input variables (Table \ref{NQI}). This approach is used in the adder presented in \cite{Ebrahimi}.
\item V3:  This intermediate approach uses multiplexers to implement subfunctions that can be derived from Table \ref{T1}. An example of subfunction is the successor function: When A = 1 and Ci = 0 then QS = (B+1) mod. 4. This approach is used in adders designed in \cite {Moaiyeri} and \cite {Roosta}. Results  \cite {Roosta} are presented with 1 and 3 power supplies.
\end{itemize}

\begin{table}
\centering
\caption{Truth table of a quaternary adder}
\begin{tabular}{|c|c|c||c|c|c|c|c|c||c|c|}
  \hline
A&B&Ci&QS&QC& &A&B&Ci&QS&QC\\
\hline
 0&0&0&0&0&&0&0&1&1&0\\
0&1&0&1&0&&0&1&1&2&0\\
0&2&0&2&0&&0&2&1&3&0\\
0&3&0&3&0&&0&3&1&0&1\\

 1&0&0&1&0&& 1&0&1&2&0\\
1&1&0&2&0&&1&1&1&3&0\\
1&2&0&3&0&&1&2&1&0&1\\
1&3&0&0&1&&1&3&1&1&1\\

 2&0&0&2&0&& 2&0&1&3&0\\
2&1&0&3&0&&2&1&1&0&1\\
2&2&0&0&1&&2&2&1&1&1\\
2&3&0&1&1&&2&3&0&2&1\\

 3&0&0&3&0&&3&0&1&0&1\\
3&1&0&0&1&&3&1&1&1&1\\
3&2&0&1&1&&3&2&1&2&1\\
3&3&0&3&1&&3&3&1&3&1\\
  \hline
\end{tabular}
\label {T1}
\end{table}

\begin{table}
\centering
\caption{NQ1, IQI and PQI funtions}
\begin{tabular}{|c||c|c|c||c|c|}
  \hline
 Q&NQI&IQI&PQI\\
\hline
 0&3&3&3\\
 1&0&3&3\\
2&0&0&3\\
3&0&0&0\\
  \hline
\end{tabular}
\label {NQI}
\end{table}

For a fair comparison, we assume that the binary full adder uses the conventional 28 T design (Figure \ref{28T}). The binary XOR gates used in V1 have 10T, which is typical of standard cell library implementation.
\begin{figure}[htbp]
\centerline{\includegraphics  [width = 7 cm]{28TBFA.png}}
\caption{28T binary full adder}
\label{28T}
\end{figure}

Table \ref{Adders} present the transistor count of the three different 4-valued full adders and the transistor count ratio with the conventional binary full adder. The TC ratio is always greater than 2, which is the information ratio. Except the quaternary adder  presented \ in \cite{Roosta} with 3 power supplies, the TC ratio is greater or equal to 4. The direct interfacing of 2-bit binary adder with 4 to 2 decoders and 2 to 4 encoder is close to the best quaternary versions. We outline that the results for version 1 are somewhat pessimistic (conventional implementation of XOR and 1-bit full adder) while version 3  \cite{Roosta} are somewhat optimistic. As detailed in \cite{etiQADD},when reducing the high fan out supported by some inverters, the TC raises to 100 (1 power supply) and 148 (3 power supplies).

\begin{table}
\centering
\caption{Adder transistor count and 4V/2V transistor ratios}
\label{Adders}
\begin{tabular}{|c|c|c|c|c|c|c|}
\hline
4-V adders&V1&V1& V2\cite{Ebrahimi}&V3\cite{Moaiyeri} &V3\cite{Roosta}&V3\cite{Roosta} \\
&3  PS&1 PS&&&3  PS&1 PS\\
\hline
TC &112 &112&111&154&82&130\\
\hline
TC ratio&4 &4&4&5.5&2.9&4.65\\
\hline
\end{tabular}
\end{table}

N-digit quaternary ripple carry adders are implementing by concatenating 1-digit quaternary adders. A N-digit quaternary adder must be compared with a 2N-bit binary adder. As the TC ratio is always greater than 2, the quaternary adders always use  two times more transistors for all published quaternary adders (except \cite{Roosta}) than the corresponding binary ones. More transistors mean more interconnects. Version 1 is a good example: the decoder and encoder circuits are just overhead over the binary implementation.
As shown in \cite{etiQADD}, quaternary carry computation of a 4-digit adder is slightly more efficient than binary carry computation of 8-bit adder, either for Carry-Look Ahead Adders or Carry-Save Adder. However the gain is too small to compensate  the transistor count difference between quaternary and binary full adders.

\subsection{Quaternary multipliers}
MVL designers generally only propose the implementation of a 1-digit quaternary multiplier. Sometimes, N*N quaternary multipliers are proposed \cite{Hajare}. In both cases, no comparison is done with the corresponding binary multipliers.
As detailed in \cite{etiQMUL}, the best implementation of a quaternary 1-digit multiplier uses the technique quoted in version 3.
The corresponding transistor count ranges from 54 T up to 76 T. It depends on the fan-out supported by some inverters. The corresponding 1-bit multiplier is implemented by an AND gate (6 T). However, a direct comparison cannot be done as an NxN combinational multiplier uses both 1-digit multipliers and 1-digit adders.
The simplest technique to implement NxN multipliers use the Wallace technique. For instance, Figure \ref{88WT} shows the Wallace reduction tree for a 8x8 bit multiplier. The situation is more complicated for the reduction of a 4x4 quaternary multiplier, as shown in Figure \ref{WQ44}:
\begin{itemize}
\item There are 16 1-digit multipliers versus 64 1-bit multipiers. However, the 1-digit multiplier is far more complicated than a simple AND gate.
\item Each line to be reduced has only 4 partial products instead of 8
\item However, a 1-digit multiplier delivers a quaternary product term and a ternary product term. With base 4, 3x3= 21, which means that the carry values are 0, 1 or 2. In Figure \ref{WQ44}, 3/2/1 corresponds respectively to quaternary/ternary/binary values. It means that different types of quaternary adders and half adders must be used in the reduction tree. They are named $Q_{332}$, $Q_{322}$,$ QHA_{32}$ and $QHA_{31}$ in Figure \ref{WQ44}. The corresponding designs are detailed  in \cite{etiQMUL}.
\end{itemize}
.
\begin{figure}[htbp]
\centerline{\includegraphics  [width =8 cm]{88WT.png}}
\caption{8*8 bit multiplication with a Wallace tree}
\label{88WT}
\end{figure}

\begin{figure}[htbp]
\centerline{\includegraphics  [width =8 cm]{WQ44.png}}
\caption{4x4 quaternary multiplier with  a Wallace Tree}
\label{WQ44}
\end{figure}

A 4x4 quaternary multiplier may also be designed by interfacing a 8x8 binary multiplier with the 4 to 2 decoder and 2 to 4 encoder circuits used for the V1 quaternary adder (subsection \ref{ENCDEC}).
With conventional binary circuits and the optimistic version of \cite{Roosta}, the transistor count comparison for 8x8 binary multipliers and two 4x4 quaternary multipliers is given in Table \ref{Multipliers}

\begin{table}
\centering
\caption{Quaternary multiplier transistor counts}
\label{Multipliers}
\begin{tabular}{|c|c|c|c|c|c|c|}
\hline
8x8 bit multiplier&4x4 Q multiplier&4x4 Q multiplier  \\&encoder and decoder&direct implementation\\\hline
1892 &2032 &2888\\
\hline
\end{tabular}
\end{table}

Significant conclusions can be derived from Table \ref{Multipliers}:
\begin{itemize}
\item Interfacing a binary multiplier with 4 to 2 decoder and 2 to 4 encoder circuits provides a better quaternary multiplier than the best direct implementation of a quaternary multiplier.
\item Compared to the binary implementation of the multiplier, the encoder and decoder circuits are just an overhead that is needed to implement the best multiplier. It means the best quaternary multiplier has more interconnects that the corresponding binary one.
\end{itemize}

\subsection{Conclusion for logic gates and arithmetic circuits}
For basic logic gates, adders and multipliers, the best quaternary circuits have always more than 2 times the transistor count of the corresponding binary ones. If a quaternary wire carries two times the information of a binary one, it doesn't mean that quaternary circuits reduce the number of interconnects. The quaternary circuits are too much more complicated and actually augment the number of interconnects.

\section{Concluding remarks}
For more than 50 years, multi-valued circuits have been proposed as a potential solution to overcome interconnection issues. 4-valued circuits could divide by two the number of interconnects. It turns out that this assumption has been proved being false. Ironically, doubling the number of interconnects by using differential signaling pairs has been the solution that has been largely used to overcome the most significant interconnect issues.
The lack of success of multivalued interconnects results from misunderstandings about both binary and multivalued circuits:
\begin{itemize}
\item The main issue of binary integrated circuits is not interconnections. It is power dissipation. While "Heat/Power wall" or "Memory wall" are well-known concepts for IC designers, "Interconnection wall" is not quoted. Increasing the number of metal layers in modern FinFET technologies reduces the interconnection issues. Reducing power consumption means using as small power supplies as possible, in the 0.8 V to 1V range. Within this range, there is no room for more than two different voltage levels.
\item As shown in this paper for quaternary circuits, m-valued circuits are far more complicated than binary circuits in terms of transistor count and thus interconnects. The increased complexity is always greater than the increased computational power, expressed as $ log_2 (m)/log_2 (2)$. Multipliers are the most spectacular example.
\end{itemize}


\section{Introduction}
Precision measurements of symmetry-violating electromagnetic moments provide stringent tests of physics beyond the Standard Model (BSM)~\cite{Bernreuther1991, Pospelov2005, Engel2013, DeMille2015, Nakai2017, Cesarotti2019}. Polar molecules containing a high-$Z$ nucleus have emerged as ideal probes because they may combine intrinsic sensitivity to CP-violating moments with structural features that aid high-precision experiments~\cite{DeMille2015, Meyer2006, Meyer2008, Hutzler2020, Hudson2011}. For example, experiments using a variety of (diatomic) polar molecules have set increasingly tight limits on the value of an electron electric dipole moment (eEDM), ultimately constraining CP-violating BSM physics at energy scales $\gtrsim$10~TeV~\cite{Hudson2011, ACME2014, ACME2018, Cairncross2017}. Molecules with enhanced sensitivity to a nuclear magnetic quadrupole moment (nMQM) have been proposed for similar tests~\cite{Maison2019, Denis2020MQM}. Concurrent with these advances in precision measurement, molecular laser cooling has recently been extended to both diatomic~\cite{Lim2018} and triatomic~\cite{Augenbraun2020YbOH} species containing the heavy Yb nucleus. Measurements combining properly chosen molecular species with laser cooling to the $\mu$K regime are predicted to provide several orders of magnitude increased sensitivity to both the eEDM and the nMQM, potentially allowing experiments to probe PeV-scale BSM physics~\cite{Kozyryev2017PolyEDM}.

A number of laser-coolable molecules have been proposed for next-generation eEDM measurements, with linear species garnering the most attention due to their structural simplicity~\cite{Kozyryev2017PolyEDM, Lim2018, Isaev2017, Gaul2019, Denis2019, Nakhate2019}. The degenerate vibrational bending mode in YbOH is a promising venue for future eEDM or nMQM measurements due to its small parity doublet, analogous to the $\Omega$-doublets that have allowed strong systematic error rejection in eEDM experiments~\cite{Meyer2008, Eckel2013, ACME2014, ACME2018, Cairncross2017}. The recent laser cooling of YbOH to $\ll$1~mK in one dimension represents an important first step toward eEDM measurements with laser-cooled polyatomic molecules~\cite{Augenbraun2020YbOH}.
However, in YbOH the energy of the first bending excitation is relatively high ($\omega_2 \sim 330$~cm$^{-1}${}) and  anharmonic~\cite{Mengesha2020}, limiting feasible coherence times to $\lesssim 800$~ms~\cite{LanChangPrivate}.

By contrast, nonlinear symmetric top molecules possess long-lived parity-doubled states due to rigid body rotation about the molecular symmetry axis (``$K$-doublets").\footnote{The quantum number $K$ denotes the molecule-frame projection of $N$, the total angular momentum excluding spin. States with $\lvert K \rvert > 0$ are doubly degenerate. Note that, necessarily, $N \geq \lvert K \rvert$.} For instance, the $K^{\dprime} = 1$ level in YbOCH$_3$ is expected to have a spontaneous lifetime much longer than 100~s due to its low energy ($\sim$10~cm$^{-1}${}) and proton spin statistics. This could lead to $> 30 \times$ improved statistical sensitivity to BSM physics if full advantage were taken of the feasible coherence times. Moreover, the very small $K$-doubling of metastable rotational levels enables full polarization and internal co-magnetometry, two features that have enabled robust rejection of systematic errors in recent eEDM measurements~\cite{Eckel2013, ACMELongPaper}. To fully leverage the lifetime of the eEDM-sensitive $K$-doublet in a neutral species requires cooling to ultracold temperatures. Laser cooling of molecules is a proven route toward this end~\cite{Barry2014, truppe2017CaF, Ding2020}. Importantly, the recent 1D laser cooling of CaOCH$_3$ to $<$1~mK shows that the complicated vibrational and rotational
structure in symmetric top molecules does not adversely affect the laser cooling process~\cite{Mitra2020}. Ytterbium monomethoxide (YbOCH$_3$), isoelectronic to CaOCH$_3$, is thus an intriguing candidate for eEDM/nMQM measurements using laser-cooled molecules, but to our knowledge there are no previously reported observations of this species.

The alkaline-earth monomethoxides, MOCH$_3$ (M = Ca, Sr, Ba), have received considerable spectroscopic attention~\cite{Wormsbecher1982, Crozet2002, Forthomme2011, OBrien1988, Namiki1998, Paul2019, brazier1986laser}. These molecules comprise an alkaline-earth atom ionically and monovalently bonded to a relatively electronegative methoxy ligand~\cite{Ellis2001, Bernath1997}. They have linear M-O-C backbones with three off-axis H atoms ($C_{3v}$ symmetry) in both the ground ($\tilde{X}\,\,{}^2 \! A_1$) and low-lying excited ($\tilde{A}\,\,{}^2 \! E$ and $\tilde{B} \, ^2A_1$) states~\cite{brazier1986laser, Crozet2005}. The unpaired valence electron is in a metal-centered, hybrid non-bonding orbital~\cite{Namiki1998} resulting in ``diagonal" vibrational branching ratios (VBRs) and short radiative lifetimes, properties highly favorable for laser cooling~\cite{Augenbraun2020ATM, kozyryev2016MOR}. It is not \textit{a priori} obvious that Yb-containing analogues will share all of these properties, e.g. due to previously observed perturbations in YbF and YbOH~\cite{Lim2017, Mengesha2020}.

Here, we present an initial experimental characterization of the previously unobserved lanthanide-containing symmetric top molecule YbOCH$_3$. This is, to date, the highest-mass neutral symmetric top molecule suitable for precision measurement experiments that has been experimentally characterized. An intense electronic transition near 579~nm has been detected and assigned to the origin band of an $\tilde{A}\,\,{}^2 \! E_{1/2} \leftarrow \tilde{X} \, ^2 A_1$ electronic transition. Dispersed laser-induced fluorescence (DLIF) measurements for YbOCH$_3$ and the isotopologue YbOCD$_3$ are used to characterize the vibrational structure in the ground and excited electronic states. In addition, we have recorded a portion of the high-resolution spectrum of the $0_0^0 \, \tilde{A}\,\,{}^2 \! E_{1/2} \leftarrow \tilde{X} \, ^2 A_1$  band. We complement these measurements with \textit{ab initio} electronic structure calculations and find excellent agreement between theory and experiment. Finally, we discuss the prospects for direct laser cooling of YbOCH$_3$ and its use in future precision measurements.

\section{Experimental Setup}
Molecular beams of YbOCH$_3$ are produced in a setup similar to that used in recent studies of YbOH~\cite{Nakhate2019, Steimle2019, Mengesha2020}. Briefly, a rotating ytterbium rod is ablated at $\sim$20~Hz with a short pulse of 532~nm radiation ($\sim10$~ns, $\sim5$ mJ). The ablation plume is entrained in and reacted with a gas mixture of methanol vapor and Ar in a supersonic expansion. The gas mixture is produced by passing Ar at $\sim 4000$~kPa over room-temperature liquid methanol (CH$_3$OH; vapor pressure $\sim10$~kPa). Typical pulse widths for the molecular beam are $\sim 50$~$\mu$s. We use a number of complementary spectroscopic methods including two-dimensional (2D) spectroscopy, DLIF spectroscopy, radiative decay, and high-resolution excitation spectroscopy to conclusively confirm detection of YbOCH$_3$ and provide initial spectroscopic characterization.

For the 2D spectroscopy~\cite{Reilly2006, Gascooke2011, Kokkin2014} and DLIF measurements, the free-jet expansion is probed approximately 10~cm downstream from the source using radiation from an excimer-pumped, tunable, pulsed dye laser ($\sim$10~ns pulse width, $\sim$3~cm$^{-1}${} linewidth). The molecular laser-induced fluorescence (LIF) is focused into a  0.67~m, high-efficiency ($f$-number = 6.0) Czerny-Turner-type monochromator with a low-dispersion grating (300~lines/mm). The DLIF from the grating is imaged on a cooled and temporally gated intensified charge-coupled device (ICCD). The CCD array is binned to produce an array of intensities versus emission wavelength. The ICCD can be gated with resolution $<$1~ns, which enables temporal separation of the LIF from background light due to the dye laser and ablation plume. The ICCD gate width is typically set to 200~ns and delayed 10~ns after the pulsed laser. The wavelength calibration and relative sensitivity of the spectrometer were calibrated prior to data collection using an argon pen lamp and blackbody source.

High-resolution measurements are performed in a separate apparatus, as in Ref.~\cite{Steimle2019}. The molecular beam is produced as described above, but skimmed in order to reduce Doppler broadening. A single-frequency, continuous-wave (cw) dye laser (linewidth $\sim 1$~MHz) is used to excite the molecules, and the resulting fluorescence is detected on a cooled photomultiplier tube (PMT). Laser powers of approximately 20~mW are used, resulting in comparable power- and Doppler-broadened linewidths, typically $\sim0.001$~cm$^{-1}${}.

\section{Observations}

\subsection{2D Spectroscopy}

\begin{figure*}[t]
\includegraphics[width=2\columnwidth]{./figures/Fig1.png}
\caption{(a) Pulsed-dye laser 2D spectrum in the vicinity of the YbOCH$_3$ $\tilde{A}\,\,{}^2 \! E_{1/2} \leftarrow \tilde{X}\,\,{}^2 \! A_1$ bandhead. Integration at fixed fluorescence wavelength can be used to obtain excitation spectra. Color scale indicates number of LIF counts detected. Integration at fixed excitation wavenumber would yield DLIF spectra, but we describe higher resolution DLIF data below. (b) Excitation spectrum obtained by vertical integration in the region labeled ``off-resonance." Strong decay to the bending mode is observed from approximately 10~cm$^{-1}${} to the blue of the origin frequency (see text). (c) Excitation spectrum by vertical integration in the region labeled ``on-resonance." A strong bandhead associated with the $0_0^0 \, \tilde{A}\,\,{}^2 \! E_{1/2} \leftarrow \tilde{X}\,\,{}^2 \! A_1$ vibronic transition is observed, characteristic of a $^2E - ^2A_1$ transition.}
\label{fig:2Dspec}
\end{figure*}

\begin{table}[tb]
\setlength{\tabcolsep}{9pt}
\begin{centering}
\renewcommand{\arraystretch}{1}
\caption{Vibrational normal mode numbering and symmetry labels for MOCH$_3$ molecules under transformations of the $C_{3v}$ point group.}
\begin{tabular}{ccl}
\hline \hline
Mode & Symmetry & Description \\
\hline
 $\nu_1$ & $a_1$ & CH$_3$ symmetric stretch \\
 $\nu_2$ & $a_1$ & CH$_3$ symmetric bend \\
 $\nu_3$ & $a_1$ & O-C stretch \\
 $\nu_4$ & $a_1$ & M-O stretch \\
 $\nu_5$ & $e$ & CH$_3$ asymmetric stretch \\
 $\nu_6$ & $e$ & CH$_3$ asymmetric bend \\
 $\nu_7$ & $e$ & O-CH$_3$ wag \\
 $\nu_8$ & $e$ & M-O-C bend\\

\hline \hline

\end{tabular}
\label{tab:VibLabels}
\end{centering}
\end{table}

Initial survey scans used 2D (excitation vs. DLIF) spectroscopy to search for YbOCH$_3$ fluorescence. A $75$~nm-wide fluorescence spectral window was monitored while scanning the excitation laser wavelength. To facilitate data analysis, this window tracked the excitation laser as it was scanned. The pulsed dye laser was scanned over a range of about 400~cm$^{-1}${} near the YbOH $\tilde{A} \, ^2 \Pi_{1/2} \leftarrow \tilde{X} \, ^2 \Sigma^+$ origin band (17323~cm$^{-1}${}~\cite{Steimle2019}). A strong and broad fluorescence signal was observed around 17271~cm$^{-1}${}. The 2D spectrum in the vicinity of this signal is shown in Fig.~\ref{fig:2Dspec}(a). The ground state vibrational frequencies observed in emission closely matched those expected of YbOCH$_3$, leading to an initial assignment of this as the YbOCH$_3$ $ 0_0^0 \, \tilde{A}\,\,{}^2 \! E_{1/2} \leftarrow \tilde{X}\,\,{}^2 \! A_1$ band. (For convenience, we include the vibrational normal mode numbering and symmetry labeling in Tab.~\ref{tab:VibLabels}.) Ytterbium monohydroxide is also produced in the ablation reaction, and by comparing the intensities of the two species's origin bands we estimate that we produced about an order of magnitude more YbOCH$_3$ than YbOH.

We obtain the excitation spectra shown in Fig.~\ref{fig:2Dspec}(b) by vertically integrating $\pm2$~nm slices of the the off-resonance fluorescence. Similarly, Fig.~\ref{fig:2Dspec}(c) is obtained by vertical integration of the on-resonance fluorescence. Horizontal integration at fixed excitation wavenumber can be used to obtain dispersed fluorescence traces, although these have lower resolution than the DLIF measurements described in Sec.~\ref{sec:ObsDLIF}.

The on-resonance detected excitation spectrum [Fig.~\ref{fig:2Dspec}(c)] exhibits an intense, sharp, blue-degraded band near 17271~cm$^{-1}${} which is assigned as the $0^0_0 \, \tilde{A}\,^2E_{1/2} \rightarrow \tilde{X}\,\,{}^2 \! A_1$ emission. The excitation spectrum is quite compact, implying that the geometry changes little upon excitation and that the electronic orbital angular momentum in the $\tilde{A} \,^2E$ state is largely unquenched (i.e., $\zeta_e \approx 1$)~\cite{Brazier1989, Marr1996}. The off-resonance excitation spectrum, shown in Fig.~\ref{fig:2Dspec}(b), is obtained by vertical integration over a $\pm$2~nm range centered $\sim$130~cm$^{-1}${} to the red of the excitation wavelength (i.e., by monitoring the Stokes-shifted emission). The extracted excitation spectrum exhibits both the sharp band near 17271~cm$^{-1}${} and a weaker, broader, and unstructured fluorescence feature near 17281~cm$^{-1}${}. The latter is assigned as excitation from $\tilde{X}\,\,{}^2 \! A_1(v=0)$ to a state of unknown character that we simply designate as [17.28]. The relative intensities of the bands at 17271~cm$^{-1}${} and 17281~cm$^{-1}${} for the on- and off-resonance detection suggest that the bending mode in the [17.28] state is quite active.

Also evident in the 2D spectrum are Stokes-shifted signals centered about 400~cm$^{-1}${} and 535~cm$^{-1}${} to the red of the excitation. These are assigned as $4^0_1\, \tilde{A}\,^2E_{1/2}\rightarrow \tilde{X}\,\,{}^2 \! A_1$ and $4^0_18^0_1 \tilde{A}\,\,{}^2 \! E_{1/2} \rightarrow \tilde{X}\,\,{}^2 \! A_1$, respectively (with similar assignments for the [17.28] bands). Again, the relative intensities of these bands suggest higher bending mode activity in the [17.28] state. It is noteworthy that in YbOH there is also an unassigned vibronic state approximately 10~cm$^{-1}${} above the $\tilde{A}\,^2\Pi_{1/2}(000)$ state~\cite{Mengesha2020}. A key difference is that in YbOCH$_3$ this higher-energy state couples preferentially to the Yb-O-C bending mode while in YbOH the analogous level exhibits a DLIF spectrum nearly identical to the diagonal origin band.

\begin{figure*}[t]
\centering
\includegraphics[scale=0.50]{./figures/Fig2.png}
\caption{(a) Pulsed-dye laser 2D spectrum in the vicinity of the $4^1_0 \, \tilde{A}\,\,{}^2 \! E_{1/2} \leftarrow \tilde{X}\,\,{}^2 \! A_1$ bands. Integration at fixed fluorescence wavelength can be used to obtain excitation spectra. The DLIF features near 17641~cm$^{-1}${} (marked ``(i)") and 17683~cm$^{-1}${} (marked ``(ii)") are shown in Fig.~\ref{fig:DLIFnu4}. The weak feature marked by an asterisk is due to the YbOH [17.68] band. (b) Pulsed-dye laser 2D spectrum in the vicinity of the $0^0_0 \, \tilde{A}\,\,{}^2 \! E_{3/2}\leftarrow \tilde{X}\,\,{}^2 \! A_1$ band.}
\label{fig:2Dnu4E32}
\end{figure*}

We also recorded 2D spectra in the vicinity of several bands involving excited vibrational levels in the $\tilde{X}\,\,{}^2 \! A_1$ and $\tilde{A}\,\,{}^2 \! E$ states. A weak feature was observed about $130$~cm$^{-1}${} to the red of the origin band near the expected $8^0_1\, \tilde{A}\,\,{}^2 \! E_{1/2} \leftarrow \tilde{X} \, ^2 A_1$ vibronic transition. We observed two bands near the expected position of the $4^1_0 \, \tilde{A}\,^2E_{1/2} \leftarrow \tilde{X}\,\,{}^2 \! A_1$ transition, approximately 400~cm$^{-1}${} above the origin band [Fig.~\ref{fig:2Dnu4E32}(a)]. These bands, near 17641~cm$^{-1}${} (designated as the [17.64] state) and 17681~cm$^{-1}${} (designated as the [17.68] state), are blue-shifted by 370~cm$^{-1}${} and 410~cm$^{-1}${} relative to the origin band. Based on the long observed progressions in the Yb-O stretching mode, these levels both appear to have strong $v_4=1$ character (see Sec.~\ref{sec:ObsDLIF}). This observation is similar to YbF, where two closely-spaced excited states (the [557] and [561] states) have substantial $v = 1$ character~\cite{Lim2017, Smallman2014}. Fluorescence from these states do not show the strong activity in the $\nu_8$ mode that was observed on the origin band.

In addition, we explored a region approximately 1300~cm$^{-1}${} above the origin band, where the $\tilde{A}\,^2E_{3/2} \leftarrow \tilde{X}\,\,{}^2 \! A_1$ transition would be expected. As shown in Fig.~\ref{fig:2Dnu4E32}(b), a band was observed near 18540~cm$^{-1}${} with weak off-diagonal decays at frequencies matching the YbOCH$_3$ $\nu_4$ and $\nu_8$ modes. This implies an effective spin-orbit constant of $a \zeta_e d \approx$1270~cm$^{-1}${}, similar to that of YbF~\cite{Dunfield1995} and YbOH~\cite{SteimleTBD}, indicating that this the spin-orbit interaction is not significantly quenched.

The $\tilde{A}\,\,{}^2 \! E_{1/2}(v=0)$ vibronic state radiative lifetime was measured by fixing the dye laser wavelength to the $0^0_0 \, \tilde{A}\,\,{}^2 \! E_{1/2} \leftarrow \tilde{X} \, ^2 A_1$ bandhead and recording the DLIF spectrum at variable time delay after the pulsed laser excitation. The ICCD gate (200~ns wide) was delayed after the pulsed dye laser in steps of $\sim3$~ns and the resulting fluorescence fit to an exponential to determine the excited state lifetime. The radiative lifetime of the $\tilde{A} \,^2E_{1/2} (v=0)$ level is determined to be 37(4)~ns, somewhat longer than the value of 22(2)~ns for the $\tilde{A}\,^2\Pi_{1/2}(000)$ state of YbOH~\cite{Mengesha2020}.

\subsection{Dispersed Fluorescence}
\label{sec:ObsDLIF}

\begin{figure*}[t]
\centering
\includegraphics[width=2\columnwidth]{./figures/Fig3.png}
\caption{DLIF spectrum resulting from excitation of the  $0^0_0 \, \tilde{A}\,\,{}^2 \! E_{1/2} \leftarrow \tilde{X} \, ^2 A_1$ bandhead. Numbers above each peak indicate transition assignments using the vibrational labeling of Tab.~\ref{tab:VibLabels}. (a) YbOCH$_3$, excitation at 17271.1~cm$^{-1}${}. Inset: expanded vertical axis to emphasize the weakest decay features observed. (b) YbOCD$_3$, excitation at 17249.7~cm$^{-1}${}. The small shoulder marked by an asterisk is attributed to impurity YbOH/YbOD.}
\label{fig:DLIF}
\end{figure*}

Higher resolution DLIF spectra were recorded by fixing the excitation laser wavelength and reducing the entrance slit of the monochromator. Typically 5,000 ablation pulses were co-added to achieve sensitive measurement of the vibrational frequencies and FCFs. Background traces were recorded separately to eliminate scattered light due to either the probe or ablation laser.

The DLIF spectrum resulting from pulsed dye laser excitation at the YbOCH$_3$ bandhead ($\sim$17271~cm$^{-1}${}) is shown in Fig.~\ref{fig:DLIF}(a). The measured vibrational intervals are listed in Tab.~\ref{tab:VibFreqs} and vibrational frequencies are extracted as discussed in Sec.~\ref{sec:DiscXState}. The DLIF is consistent with a relatively diagonal FCF matrix, as expected. A short progression in $\nu_4$ (the Yb-O stretching mode) is observed, with relative intensities similar to those of YbOH~\cite{Mengesha2020}. We observed prominent decay to $\nu_8$ (Yb-O-C bending mode), which was unexpected because the $8^0_1 \, \tilde{A}\,^2E_{1/2} \rightarrow \tilde{X}\,\,{}^2 \! A_1$ transition is nominally forbidden by symmetry. Vibronic coupling due to the (pseudo)-Jahn-Teller interaction has previously been invoked to explain similar decays in CaOCH$_3$~\cite{kozyryevCaOCH3}. We do not observe fluorescence at $2\nu_8$, which is symmetry allowed. Weak features associated with the $\nu_3$ (C-O stretch) are also observed, and the frequencies and branching ratios are consistent with values observed in CaOCH$_3$~\cite{kozyryevCaOCH3, Paul2019}.

Isotopic studies aided in the assignment of the ligand-centered vibrations. Figure~\ref{fig:DLIF}(b) displays the DLIF spectrum following excitation in YbOCD$_3$ and vibrational frequencies extracted from these data are listed in Tab.~\ref{tab:VibFreqs}. The origin band of YbOCD$_3$ was found near 17250~cm$^{-1}${}. This $\approx 20$~cm$^{-1}${} isotope shift is similar to the large shift ($\approx 10$~cm$^{-1}${}) observed upon deuteration of CaOCH$_3$~\cite{Crozet2005}. Interestingly, the branching ratios observed for YbOCD$_3$ differ strikingly from those observed in YbOCH$_3$, although the frequencies agree well with the expected isotope shifts.

We have also measured DLIF spectra following excitation to bands near the expected position of the $\tilde{A}\,^2E_{1/2}(v_4=1)$ state. As described above, two such states were found near 17641~cm$^{-1}${} and 17683~cm$^{-1}${}, which we call the [17.64] and [17.68] states, respectively. DLIF measurements from these levels are shown in Figs.~\ref{fig:DLIFnu4}(a) and (b). These spectra both show strong decays to a progression in $\nu_4$, confirming large $v_4=1$ character of each excited state. Weak features associated with decay to $v_8=1$ are also observed, although with significantly lower intensity than observed on the origin band. The presence of two states with significant $v_4=1$ character is likely due to mixing of the $\tilde{A}\,^2E_{1/2} (v_4=1)$ vibronic level with vibronic levels of states associated with $f$-orbital vacancies, e.g. ([Xe]$4f^{13}6s^2$)$_\text{Yb$^+$}$, similar to the case of YbF~\cite{Lim2017, Smallman2014}. This spectral region is complicated by the presence of YbOH bands within a few cm$^{-1}$ of both the [17.64] and [17.68] states. We attribute a few small background features to these YbOH bands, which fluoresce at a range of wavelengths including $\sim$577~nm and $\sim$582~nm~\cite{Mengesha2020}.

\begin{figure*}[tb]
\centering
\includegraphics[width=2\columnwidth]{./figures/Fig4.png}
\caption{DLIF spectrum resulting from excitation at (a) 17641.75~cm$^{-1}${} and (b) 17683.40~cm$^{-1}${}. The numbers above the spectral features are the measured shifts (in cm$^{-1}$) relative to the excitation wavelength (``Exc."). Features marked by an asterisk are due to YbOH impurity fluorescence as described in the text.}
\label{fig:DLIFnu4}
\end{figure*}

\begin{table*}[t]
\setlength{\tabcolsep}{9pt}
\begin{centering}
\renewcommand{\arraystretch}{1}
\caption{Observed and calculated (\textit{ab initio}) vibrational intervals in $\tilde{X}\,\,{}^2 \! A_1$ for YbOCH$_3$ and YbOCD$_3$. Calculated values are in the harmonic approximation. Previously measured values for alkaline-earth monomethoxides are included for comparison. All frequencies are in cm$^{-1}$. Typical measurement error bars are $\pm$5~cm$^{-1}${} depending upon signal-to-noise ratio and proximity to argon emission calibration lines.}
\begin{tabular}{ccccccc}
\hline \hline
Mode & YbOCH$_3$ (meas.) & YbOCH$_3$ (calc.) & YbOCD$_3$ (meas.) & YbOCD$_3$ (calc.) & SrOCH$_3$ $^{a}$& CaOCH$_3$ $^{b}$  \\
\hline

 $\nu_3$ &   1151 & 1154 & - & 1085 & 1138 & 1156 \\
 $\nu_4$ &  400 & 390 & 370 & 369 & 405 & 487   \\
 $\nu_8$ & 130 & 134 & 120 & 124 & 135 &   142\\
 $2\nu_8$ & 260 & 268 & - & 248 & - &   -\\
 $\nu_4 + \nu_8$ &  533 & 524 & 500 & 493 &- & 625 \\
 $2\nu_4$ &  805 & 780 & - & 738 & 806&973 \\
 $2\nu_4 + \nu_8$ &  933 & 914 & 880 & 862 & -&- \\
 $\nu_4 + 2\nu_8$ & - & 658 & 680 & 617 &- &- \\
 $\nu_3 + \nu_8$ & 1287 &  1288 & 1270& 1262 & -&-\\
 $3\nu_4$ & 1197 & 1170 & - & - & - & -\\
 $4\nu_4$ & 1576 & 1560 & - & - & - & -\\
 $5\nu_4$ & 1972 & 1950 & - & - & - & -\\

\hline \hline
\multicolumn{5}{l}{$^{a}$ Measured values from Refs.~\cite{OBrien1988} and \cite{brazier1986laser}}\\
\multicolumn{5}{l}{$^{b}$ Measured values from Ref.~\cite{brazier1986laser}}\\
\end{tabular}
\label{tab:VibFreqs}
\end{centering}
\end{table*}

\subsection{CW Excitation and DLIF Spectra}
The high-resolution spectrum recorded near the $0^0_0 \, \tilde{A}\,\,{}^2 \! E_{1/2} \leftarrow \tilde{X}\,\,{}^2 \! A_1$ bandhead identified above is shown in Fig.~\ref{fig:HiRes}(a). The spectrum is quite congested due to the presence of several isotopes of Yb with relatively high natural abundance, including $^{174}$Yb (32\%), $^{172}$Yb (22\%), $^{173}$Yb (16\%), $^{171}$Yb (14\%), and $^{176}$Yb (13\%). A full analysis of the high-resolution spectra will be the focus of a future publication, where lower rotational temperature and isotope-selective enhanced production~\cite{Bernath1997, Jadbabaie2020} is used to greatly reduce the complexity of the spectrum. We make some qualitative observations here, showing the spectrum is consistent with a $^2E_{1/2} \leftarrow ^2A_1$ transition.

\begin{figure*}[t]
\centering
\includegraphics[width=2\columnwidth]{./figures/Fig5.png}
\caption{(a) High-resolution excitation spectrum near the YbOCH$_3$ $0_0^0 \, \tilde{A}\,^2E_{1/2} \leftarrow \tilde{X} \, ^2 A_1$ bandhead. Features marked by an asterisk were selected for cw-excitation DLIF measurements. Inset: Details of spectrum near 17273.15~cm$^{-1}${}, with isotopic and $K$-stack structure visible. (b) DILF spectrum resulting from rotationally-resolved, cw-dye laser excitation at 17270.970~cm$^{-1}${}. Inset: Zoom near signal baseline to show the relative intensity sensitivity of the measurement.}
\label{fig:HiRes}
\end{figure*}

The rotational energy level structure of a $^2A_1$ state is that of an open-shell prolate symmetric top molecule near the Hund's case (b) limit. For an MOCH$_3$ species, the $K_R^{\dprime}=1$ stack of levels begins approximately 6~cm$^{-1}${} ($\approx A^{\dprime}$) above the $K_R^{\dprime}=0$ stack. The $\lvert K_R^{\dprime} \rvert =1$ levels obey different nuclear spin statistics and will also be populated significantly and show prominently in the spectrum. The energy level pattern of a $^2E_{1/2}$ state is that of an open-shell prolate symmetric top near the Hund's case (a) limit. The parity doubling in the $K_R^\prime=0, K^\prime=1$ level is expected to be large and hence the dominant $\tilde{A}\,\,{}^2 \! E_{1/2}(K_R^\prime=0,K^\prime=1) \leftarrow \tilde{X} \,^2A_1(K_R^{\dprime}=0)$ band is expected to have the same general appearance as the $\tilde{A}\,^2\Pi_{1/2} \leftarrow \tilde{X}\,^2\Sigma^+$ band of YbOH. Indeed, we observe that the high-resolution spectrum is reminiscent of the $\tilde{A} \, ^2\Pi_{1/2} \leftarrow \tilde{X} \, ^2\Sigma^+$ origin band of YbOH due to the sharply peaked, blue-degraded bandhead and the widely spaced peaks in the $P$/$R$ branches~\cite{Steimle2019, Mengesha2020}. Using the expected YbOCH$_3$ rotational constant, $B$, determined by either empirically scaling parameters for YbOH or by \textit{ab initio} calculations, we estimate that these peaks are spaced by approximately $4B$. This branch spacing implies an excited state parity doubling $\epsilon_1 - 2h_1 \approx B$~\cite{HerzbergVol1}.\footnote{In this expression, $\epsilon_1$ is a higher-order spin-rotation parameter and $h_1$ is a Jahn-Teller parameter~\cite{Crozet2005}.} Here, $\epsilon_1 - 2h_1$ plays the role of the usual $p+2q$ ($\Lambda$-type doubling) parameter in a linear $^2\Pi$ state~\cite{Dick2006SrCH3}. Note that a similar branch spacing of about $4B$ is observed in the YbOH $\tilde{A} \, ^2\Pi_{1/2} \leftarrow \tilde{X} \, ^2\Sigma^{+}$ origin band~\cite{Steimle2019}.

We recorded several DLIF spectra using single-frequency, cw-dye laser excitation on selected rotational lines. These measurements were used to eliminate the possibility that the strong symmetry-forbidden decays observed above originate from the pulsed laser simultaneously exciting many quantum states and/or isotopologues. We selected a few rotationally resolved transitions for cw-excitation DLIF spectra. The selected transitions, marked by asterisks in Fig.~\ref{fig:HiRes}(a), were chosen based on their relative isolation and preliminary assignment of $K^{\prime}=0$ and $K^{\prime}=1$ character. A representative cw-excitation DLIF spectrum is shown in Fig.~\ref{fig:HiRes}(b). The sensitivity of the cw-dye laser measurements is lower than that of the pulsed dye laser measurements due to decreased signal and an inability to temporally exclude the excitation light. This loss in sensitivity required enlarging the entrance slit on the monochromator, lowering the spatial resolution. DLIF spectra recorded on other rotationally-resolved transitions yielded VBRs consistent with those of Fig.~\ref{fig:HiRes}(b) and with the pulsed-laser excitation spectra. These measurements confirmed that the VBRs determined from the pulsed-dye laser excitation are not strongly dependent on the rotational state selected.

\section{Analysis}
\subsection{Vibrational Structure}
The wavelength axes of the DLIF spectra were calibrated and linearized prior to data collection. Ground-state vibrational intervals were determined from the DLIF spectra by locally fitting each peak and taking appropriate differences. The measured frequencies and assignments are listed in Table~\ref{tab:VibFreqs} along with comparisons to the alkaline-earth monomethoxides SrOCH$_3$ and CaOCH$_3$. The observed ground state vibrational energies were fit to the phenomenological expression~\cite{HerzbergVol3}

\begin{equation}
    E(v) = \sum_{i} \omega_i \left(v_i + \frac{d_i}{2} \right) + \sum_{i,k} x_{ik} \left(v_i + \frac{d_i}{2} \right)^2,
    \label{eq:VibFreqs}
\end{equation}

\noindent where mode $i$ has degeneracy $d_i$, frequency $\omega_i$, and $v_i$ quanta of excitation. The constants $x_{ik}$ denote the leading anharmonic corrections. Only the observed modes ($\nu_3$, $\nu_4$, and $\nu_8$)  were included in the fit. The fit incorporated DLIF traces with pulsed-dye laser excitation to the $\tilde{A}(v=0)$, [17.64], and [17.68] states. The results of this fit are presented in Sec.~\ref{sec:DiscXState}.

\subsection{Vibrational Branching Ratios}
Vibrational branching ratios (VBRs), $b_{v^\prime, v^{\dprime}}$, were determined by the ratios of integrated areas under each observed peak. These are related to the FCFs by~\cite{nguyen2018fluorescence}

\begin{equation}
    \begin{split}
    b_{v^\prime,v^{\dprime}} &= \frac{I_{iv^\prime,fv^{\dprime}}}{\sum_{fv^{\dprime}} I_{iv^\prime,fv^{\dprime}}} \\ & \approx \frac{\text{FCF}_{iv^\prime,fv^{\dprime}} \times \nu^3_{iv^\prime,fv^{\dprime}}}{\sum_{fv^{\dprime}}\text{FCF}_{iv^\prime,fv^{\dprime}} \times \nu^3_{iv^\prime,fv^{\dprime}}},
    \end{split}
\label{eq:VBRFCF}
\end{equation}

\noindent where $i$ and $f$ label initial and final states, respectively. Table~\ref{tab:VBRs} lists the observed vibrational branching ratios and FCFs determined for the YbOCH$_3$ $\tilde{A}\,\,{}^2 \! E_{1/2}-\tilde{X}\,\,{}^2 \! A_1$, [17.64]$-\tilde{X}\,\,{}^2 \! A_1$, and [17.68]$-\tilde{X}\,\,{}^2 \! A_1$ bands. The spectrometer and ICCD relative intensity sensitivities were calibrated prior to data collected. The typical noise level of the DLIF spectra corresponds to $\sim$100 counts per pixel, equivalent to $\sim$1 part in $10^3$ relative to the dominant peak. Systematic errors associated with the calibration lead to $\sim$5\% relative uncertainty in peak height near 580~nm and $\sim$10\% relative uncertainty near 630~nm.

\subsection{Electronic Structure Calculations}
To aid assignment of the DLIF spectra, \textit{ab initio} calculations were performed for YbOCH$_3$ and YbOCD$_3$ using the ORCA quantum chemistry program~\cite{neese2012orca, neese2020orca}. Molecular orbitals, optimized geometries, normal modes, and vibrational frequencies for the ground electronic state were calculated at the level of unrestricted Kohn-Sham (UKS) density-functional theory (DFT) using the B3LYP functional~\cite{Bykov2015, Stephens1994}. We carefully studied the accuracy of various basis set and functional choices. Final calculations employed a quadruple-$\zeta$ basis set, including two sets of polarization functions and diffuse orbitals on the heavy atoms. A 28-electron small-core pseudopotential was used for the Yb atom~\cite{Dolg1989, ORourke2019, isaev2015polyatomic}. Time-dependent DFT (TD-DFT) was used to model the $\tilde{A}\,\,{}^2 \! E$ state, using the same basis and functional as were used for the ground electronic state. This allowed computation of excitation energies and geometry optimization using analytic Hessians. Vibrational frequencies were determined via numerical differentiation. The theoretical methods were validated by performing test calculations for YbOH and CaOCH$_3$. See App.~\ref{app:AbInitio} for comparison between different basis sets, validation against previous measurements, and additional computational outputs.

\begin{table*}[tb]
\setlength{\tabcolsep}{9pt}
\begin{centering}
\renewcommand{\arraystretch}{1}
\caption{Branching ratios and FCFs from the $\tilde{A}\,^2E_{1/2} (v=0)$, [17.64], and [17.68] states of YbOCH$_3$ to the $\tilde{X}\,^2A_1$ levels. Values are determined from the pulsed-dye laser excitation DLIF spectra. After accounting for the detection noise and systematic uncertainties due to calibration, the noise levels of the smallest VBRs are $\pm0.7\times10^{-3}$ for the $0^0_0$ band and $\pm2\times10^{-3}$ for the [17.64] and [17.68] bands.}
\begin{tabular}{ccccccc}
\hline \hline
 & \multicolumn{2}{c}{$\tilde{A}\,^2E_{1/2}(v=0)$} & \multicolumn{2}{c}{[17.64]} & \multicolumn{2}{c}{[17.68]} \\
 & VBR & FCF & VBR & FCF & VBR & FCF \\ \hline
$0_0$ & 0.7184 & 0.729 & 0.132 & 0.122 & 0.526 & 0.504 \\
$8_1$ & 0.1653 & 0.164 & - & - & - & - \\
$4_1$ & 0.0566 & 0.054 & 0.624 & 0.617 & 0.375 & 0.385\\
$4_1 8_1$ & 0.0413 & 0.038 & 0.0181 & 0.0183 & 0.044 & 0.046 \\
$4_2$ & $1.1\times10^{-3}$ & $9.89\times 10^{-4}$ & 0.201 & 0.214 & 0.0268 & 0.0293 \\
$4_2 8_1$ & $7.9\times10^{-3}$ & $6.8 \times 10^{-3}$ & - & - & - & - \\
$3_1$ & $6.7\times10^{-3}$ & $5.49 \times 10^{-3}$ & - & - & - & - \\
$3_1 8_1$ & $2.65\times10^{-3}$ & $2.13 \times 10^{-3}$ & - & - & - & - \\
$4_3$ & - & - & 0.021 & 0.0244 & $5.13\times10^{-3}$ & $6.02 \times 10^{-3}$\\
$4_4$ & - & - & $3.6\times10^{-3}$ & $4.36 \times 10^{-3}$ & 0.015 & 0.0194 \\
$4_5$ & - & - & - &                          & $7.6\times10^{-3}$ & $1.0 \times 10^{-2}$ \\

\hline \hline

\end{tabular}
\label{tab:VBRs}
\end{centering}
\end{table*}

\section{Discussion}
The primary purpose of this investigation was to observe YbOCH$_3$ and determine its suitability for laser cooling and precision tests of fundamental physics. The short radiative lifetime and reasonably diagonal FCFs indicate that YbOCH$_3$ is a promising candidate for these applications, similar to YbF~\cite{Sauer1996,Sauer1999,Lim2017} and YbOH~\cite{Nakhate2019}. The data provide information necessary to begin further studies of optical cycling and vibrational repumping. As discussed below, YbOCH$_3$ also presents an interesting platform for future studies of (pseudo-)Jahn-Teller interaction in the presence of strong spin-orbit coupling.

\subsection{The $\tilde{X}\,\,{}^2 \! A_1$ State}
\label{sec:DiscXState}
Assignments of the observed vibrational decays were relatively straightforward using the \textit{ab initio} predictions, isotopic data, and comparison to alkaline-earth monomethoxides. For the origin band, the peaks redshifted by 130~cm$^{-1}${} and 400~cm$^{-1}${} are assigned to  $8^0_1$ and $4^0_1$, respectively. The feature redshifted by 1150~cm$^{-1}${} is near the expected frequency of either $\nu_3$ (1$\times$ O-C stretch) or $3\nu_4$ (3$\times$ Yb-O stretch). Because this decay has a stronger intensity than $4^0_2$, we favor the assignment to $3^0_1$. Furthermore, the DLIF spectra from the [17.64] and [17.68] states show a feature at slightly higher frequency ($\approx 1170$~cm$^{-1}${}) that is assigned to $3\nu_4$. The weak feature redshifted by $\sim$1280~cm$^{-1}${} is assigned to the decay $3^0_1 8^0_1$ based on the predicted frequencies for various combination bands. In addition, the relative intensity of this feature to the $3^0_1$ peak is roughly consistent with the strengths of other combination bands involving $\nu_8$.

The vibrational frequencies are determined by fitting observed vibrational energy intervals to Eq.~\ref{eq:VibFreqs}. We find $\omega_3 = 1152(5)$~cm$^{-1}${}, $\omega_4 = 405(3)$~cm$^{-1}${}, $\omega_8 = 130(3)$~cm$^{-1}${}, and $x_{44} = -2.0(6)$~cm$^{-1}${}. These fitted frequencies agree well with the computed harmonic frequencies and with the corresponding values for CaOCH$_3$ and SrOCH$_3$ after accounting for differences in reduced mass~\cite{brazier1986laser}. The anharmonic contribution to the Yb-O stretching mode is similar to that observed in YbOH in both magnitude and sign~\cite{Mengesha2020}. Finally, the fitted value $\omega_3$ is similar to that predicted for the methoxy anion, consistent with a highly ionic metal-ligand bond~\cite{Weil1985}.

Our \textit{ab initio} calculations predict that the $\tilde{X}\,\,{}^2 \! A_1$ state has a linear Yb-O-C bond and symmetric, off-axis H atoms with $C_{3v}$ symmetry, just like the alkaline-earth monomethoxides. The calculated values $r_{\rm OC} = 1.39$~\AA~and $r_{\rm CH} = 1.10$~\AA~for YbOCH$_3$ are quite similar to those measured in CaOCH$_3$ ($r_{\rm OC} = 1.411(7)$~\AA~and $r_{\rm CH}=1.094$~\AA)~\cite{Crozet2005}, while the calculated value of $r_{\rm YbO} = 2.04$~\AA~for YbOCH$_3$ is very close to the measured value for YbOH ($r_{\rm YbO} = 2.0397$~\AA)~\cite{Nakhate2019}. In addition, our \textit{ab initio} calculations predict a molecule-frame dipole moment of $\mu(\tilde{X}) = 2.2$~D in the $\tilde{X}\,\,{}^2 \! A_1$ ground state. This is nearly identical to the ground-state dipole moment of YbOH (1.9(2)~D~\cite{Steimle2019}) and comparable to that of CaOCH$_3$ (1.58(8)~D~\cite{Namiki1998}). These calculations also report the $\tilde{X}\,\,{}^2 \! A_1$ electronic state arises primarily from Yb-centered $6s\sigma$ (65\%) and $6p\sigma$ (15\%) orbitals. This is quite similar to YbF where experimental analysis of the hyperfine structure showed that the $X \, ^2\Sigma^+$ state had 57\% $6s$ character~\cite{Steimle2007}. The resulting polarization of the valence electron away from the Yb-O bond is consistent with the predicted molecule-frame dipole moment being significantly smaller than that of an electrostatic model with point charges located near the Yb$^{+}$ and OCH$_3^{-}$ moieties~\cite{Marr1996}. Together, these results help confirm the expectation that YbOCH$_3$ in the $\tilde{X}\,\,{}^2 \! A_1$ ground state conforms to many of the patterns established by previous studies of alkaline-earth monomethoxides. The predicted vibrational frequencies for YbOCH$_3$ and YbOCD$_3$ show excellent agreement with our measurements, as can be seen in Tab.~\ref{tab:VibFreqs}. Although we only observed activity in three of the vibrational modes, we report the full results of the vibrational structure in App.~\ref{app:AbInitio}.

\subsection{The Excited States}
\label{sec:DiscAState}
We have observed a number of excited states assigned to YbOCH$_3$, associated with bands near 17271~cm$^{-1}${}, 17281~cm$^{-1}${}, 17641~cm$^{-1}${}, 17683~cm$^{-1}${}, and 18542~cm$^{-1}${}. The band at 17271~cm$^{-1}${} is assigned to terminate on the $\tilde{A}\,\,{}^2 \! E_{1/2} (v=0)$ level. As discussed above, both the [17.64] and [17.68] states seem to have strong $v_4=1$ (Yb-O stretching) character, but on the basis of the FCF measurements we preliminarily assign the [17.64] state as dominantly $\tilde{A}\,\,{}^2 \! E_{1/2}(v_4=1)$, while the [17.68] state appears to be a perturbing state similar to those observed in both YbF~\cite{Lim2017} and YbOH~\cite{Mengesha2020}. The band near 18542~cm$^{-1}${} is tentatively assigned as $0^0_0 \, \tilde{A}\,^2E_{3/2} \leftarrow \tilde{X}\,\,{}^2 \! A_1$. This implies an effective spin-orbit splitting $a \zeta_e d \approx 1270$~cm$^{-1}${}, similar to that of YbF~\cite{Dunfield1995} and YbOH~\cite{SteimleTBD}. Future high-resolution studies will be useful to provide further confirmation of these assignments.

Observation of the relatively intense $8^0_1$ feature in the $\tilde{A}\,^2E_{1/2}(v=0)$ DLIF spectrum (Fig.~\ref{fig:DLIF}) was at first unexpected because, assuming both the $\tilde{X}$ and $\tilde{A}$ states have $C_{3v}$ symmetry, the $\nu_8$  mode is of $e$ symmetry and this decay is symmetry-forbidden within the Born-Oppenheimer (BO) approximation. One possible explanation for the strength of this transition is that YbOCH$_3$ possesses a bent geometry ($C_s$ symmetry) in the $\tilde{A}$ state. We disfavor this explanation for several reasons. First, we do \textit{not} observe emission from $\tilde{A}\,\,{}^2 \! E_{1/2}$ to $\tilde{X}\,\,{}^2 \! A_1 (v_8=2)$, which would be expected if there were significant geometrical change upon electronic excitation. Similarly, we do not observe strong decays to $\tilde{X}\,\,{}^2 \! A_1 (v_8=1)$ from the [17.64] or [17.68] states. Second, a bent excited state would likely be accompanied by significant quenching of the electronic orbital angular momentum, but the above assignment of the $\tilde{A}\,^2E_{3/2}$ state is consistent with $a \zeta_e d \approx 1270$~cm$^{-1}${}, indicating little quenching. Third, we did not observe rotational bands associated with the $K_a$ structure expected for a bent molecule; the contour of the high-resolution spectra could be well simulated using a symmetric top Hamiltonian. A more detailed discussion of the high-resolution spectrum will be the subject of a separate publication.

Instead, we attribute the relatively strong $8^0_1$ signal to spin-orbit vibronic coupling. Following Ref.~\cite{Paul2019}, we consider the possibility of Jahn-Teller coupling within $\tilde{A}\,\,{}^2 \! E$, characterized by parameter $k_8$, and pseudo-Jahn-Teller coupling between $\tilde{A}\,\,{}^2 \! E$ and nearby states of $A_1$ symmetry, characterized by parameter $\lambda_8$. Here, we assume that $\nu_8$ is the only mode with significant Jahn-Teller activity. Spin-orbit coupling between $\tilde{A}\,\,{}^2 \! E$ and $\tilde{B}\,^2A_1$ can interfere with (pseudo-)Jahn-Teller terms and lead to a second-order, effective coupling between the $\tilde{A}\,\,{}^2 \! E(v_8=1)$ and $\tilde{A}\,\,{}^2 \! E(v_8=0)$ vibronic levels. Given the intensity of the $8^0_1$ feature, we infer an effective matrix element $$\langle \tilde{A}\,\,{}^2 \! E_{1/2},v_8=1 \rvert H_{JT} \lvert \tilde{A}\,\,{}^2 \! E_{1/2},v_8=0 \rangle \approx 50~\text{cm$^{-1}${}}.$$ Again following Ref.~\cite{Paul2019}, we may also express this matrix element as $k_8 + 2C_{AB}\lambda_8$. The first term represents direct coupling between $\tilde{A}\,\,{}^2 \! E(v_8=0,1)$ while the second term represents the spin-orbit vibronic coupling via an intermediate $^2A_1$ state. $C_{AB}$ is the mixing induced by spin-orbit coupling between $\tilde{A}\,\,{}^2 \! E$ and $\tilde{B}\,^2A_1$, approximated using the measured spin-orbit parameter of YbOCH$_3$ and the estimated $\tilde{B}\,^2A_1$ energy scaled from YbF. It has been shown for CaOCH$_3$ that the $k_8$ term is negligible compared to the $\lambda_8$ term. Assuming the same holds for YbOCH$_3$, we estimate $\lambda_8 \approx 55$~cm$^{-1}${}. This is approximately three times smaller than the value of $\lambda_8$ computed by \textit{ab initio} methods for CaOCH$_3$, which is consistent with the expected trend as spin-orbit coupling strength increases~\cite{Barckholtz1998}. Though preliminary, this estimate shows that the strong $8^0_1$ feature can be explained reasonably by the spin-orbit vibronic coupling mechanism.

Recent studies of CaOCH$_3$ by our group lend support to these quantitative estimates. In those measurements, which will be discussed in more detail elsewhere, we observed the analogous vibronically-induced $8^0_1$ feature in CaOCH$_3$, albeit with weaker VBR $\approx 3 \times 10^{-3}$. Note that this agrees quite well with the \textit{ab initio} prediction of Ref.~\cite{Paul2019}, where the decay was predicted but not observed. Applying the same analysis as used above, we determine a linear pseudo-Jahn-Teller coefficient for CaOCH$_3$ of $\lambda_{8}^\text{CaOCH$_3$} \approx 145$~cm$^{-1}${}. This is in good agreement with the \textit{ab initio} value of $150$~cm$^{-1}${}~\cite{Paul2019} and provides a satisfactory check of the estimation methods used above. Nonetheless, \textit{ab initio} treatment of the vibronic coupling problem in YbOCH$_3$ would be highly desirable. In addition, we note that other interpretations are possibly valid, such as the $k_8$ parameter being nonnegligible or vibronic coupling between $\tilde{A}\,\,{}^2 \! E$ and nearby electronic states arising from excitations of Yb $4f$ electrons. The presence of a strong perturbation makes it very challenging to use standard semi-empirical methods~\cite{Mengesha2020,kozyryevCaOCH3, Augenbraun2020ATM} to predict the FCFs in YbOCH$_3$. This reinforces the need for careful experimental studies to guide future experiments with this molecule.

TD-DFT calculations targeting the $\tilde{A}\,\,{}^2 \! E$ state predict an excitation energy $T_e = 18060$~cm$^{-1}${}. This is in good agreement with the measured value of $T_e \approx 17900$~cm$^{-1}${} (accounting for the spin-orbit coupling). The $\tilde{A}\,\,{}^2 \! E$ state is predicted to have $C_{3v}$ symmetry, with bond lengths/angles $r_{\rm YbO} = 2.010$~\AA, $r_{\rm OC} = 1.396$~\AA, $r_{\rm CH} = 1.095$~\AA, and $\theta_\text{OCH} = 111.283^\circ$. The calculations predict a linear Yb-O-C structure. The predicted change in Yb-O bond length is quite similar to that observed for YbOH and consistent with the relatively diagonal branching ratios~\cite{Steimle2019}. The calculation predicts a vibrational frequency $\omega_4 = 410$~cm$^{-1}${}, which is in reasonable agreement with the splitting observed between the $\tilde{A}\,\,{}^2 \! E$ and [17.64]/[17.68] states. See App.~\ref{app:AbInitio} for a full set of predicted excited-state vibrational frequencies.

\subsection{Prospects for Laser Cooling and Trapping}

Our measurements confirm the viability of achieving rapid photon cycling with YbOCH$_3$. The short radiative lifetime of the $\tilde{A}$ state (corresponding to $\Gamma \approx  27 \times 10^6$~s$^{-1}$) will enable large optical forces, comparable to those of previously laser-cooled molecules. Based on our VBR measurements, a feasible laser cooling scheme requiring about 8 distinct laser wavelengths will achieve vibrational closure below a part in $10^3$.  Additional frequencies required for repumping rotational and spin-rotation substructure can be added as laser frequency sidebands using electro- or acousto-optic modulators. Assuming a cryogenic buffer-gas beam containing $10^9$ YbOCH$_3$ molecules and a loss probability per photon scatter of $<1 \times 10^{-3}$ (the smallest observed VBR) to unaddressed vibrational states, this scheme would leave $\sim 10^4 - 10^5$ YbOCH$_3$ molecules in a ``bright" state following $\sim$10,000 photon scatters. This is comparable to the experimental complexity required to laser cool much lighter (and simpler) species~\cite{Baum2020PRA, Mitra2020, kozyryevCaOCH3}. Laser-enhanced molecule production such as that demonstrated for YbOH~\cite{Jadbabaie2020} could increase the molecule number by at least an order of magnitude. The presence of multiple excited electronic states with strong electric-dipole allowed transitions to excited vibrational states in $\tilde{X}$ provides a convenient method of decoupling the repumping pathways, similar to the scheme used for YbF~\cite{Lim2018}.

\subsection{Suitability for Precision Measurements}
\label{sec:DiscPrecisionMeasurements}
To the best of our knowledge, YbOCH$_3$ is the highest-$Z$ symmetric top molecule amenable to optical cycling to be experimentally characterized. It is a promising system for precision tests of parity and/or time-reversal symmetry violation. In addition to the high-$Z$ nucleus that leads to a valence electron experiencing strong relativistic effects, the presence of easily polarized, internal co-magnetometer states in the $\lvert K^{\dprime} \rvert \neq 0$ manifolds allow full access to a large internal effective electric field and rejection of potential systematic errors~\cite{Kozyryev2017PolyEDM, Hutzler2020}. Consider, for example, laser-cooled and trapped YbOCH$_3$ used for future tests of the eEDM. Under the reasonable assumptions that the internal effective electric field of YbOCH$_3$ matches that of YbOH ($\mathcal{E}_\text{eff} \approx 25$~GV/cm)~\cite{Denis2019, Gaul2019, Prasannaa2019} and that achievable experimental lifetimes approach $10-100$~s~\cite{OHara1999, Bernon2013, Wang2018}, $10^5$ trapped YbOCH$_3$ molecules and one week of averaging could provide a statistical sensitivity to the eEDM four orders of magnitude beyond the current limit~\cite{ACME2018}. This is about an order of magnitude more sensitive than an equivalent measurement in YbOH.

Even with a much smaller photon budget, transverse laser cooling and high-fidelity readout enabled by rapid photon cycling could lead to new limits on precision tests of $CP$-violation~\cite{Ho2020}. Consider an experiment using $^{173}$YbOCH$_3$ in a molecular beam to search for a nMQM. Combining the demonstrated production of YbOCH$_3$, enhanced chemical production~\cite{Jadbabaie2020}, and transverse laser cooling using $<$150 photons~\cite{Mitra2020}, bright beams with $>10^7$ molecules per pulse in the science state appear possible. Forward velocities around 50~m/s and laser cooling to $< 1$~mK in the transverse direction~\cite{Augenbraun2020YbOH, Mitra2020} would make coherence times on the order of 10~ms feasible. Under these conditions, a sensitivity on the order of 10~$\mu$Hz would be achievable with one week of averaging. Assuming (reasonably) that the nMQM sensitivity parameter is similar to that of YbOH, this would be near or below the level required to set new limits on $T$-violating BSM physics, e.g., the QCD $\bar{\theta}$ parameter or the difference of quark chromo-EDMs~\cite{Maison2019, Denis2020MQM}. Note that symmetric top molecules offer some intrinsic advantages over linear species for such measurements such as requiring smaller laboratory electric fields for full polarization and that nuclear spin statistics ease the task of state preparation by naturally populating the $\lvert K^{\dprime}\rvert=1$ science state~\cite{Kozyryev2017PolyEDM, Yu2020}.

YbOCH$_3$ is also an interesting candidate for possible tests of $P$-violation. Recently, the degenerate bending modes in open-shell linear triatomic molecules have been proposed for tests of nuclear spin-dependent parity violation~\cite{Norrgard2019, Hao2020}. Similar searches could be conducted in symmetric top molecules that meet certain structural conditions. One key requirement of such an experiment is that the ground-state parity doubling is not smaller than all spin-rotation/hyperfine structure, such that opposite-parity states with the same total angular momentum may be Zeeman-tuned into near resonance. In $^{171}$YbOCH$_3$, we expect that hyperfine and $K$-doublet splittings will be of the same order of magnitude, $\sim$0.3~MHz~\cite{Namiki1998}. The size of the $K$-doubling increases with $N$, so that by working in, e.g., the $N^{\dprime}=2$ or $3$ state it may be possible to tune the parity-doubling to be larger than the hyperfine structure. Future measurements with hyperfine resolution will be required to determine the exact splittings. Alternatively, chiral variants such as YbOCHDT may be useful for studies of parity violation arising from $P$-odd cosmic fields, e.g., in searches for axion-like particles~\cite{Gaul2020a, Gaul2020b}. Chiral species are predicted to be particularly sensitive to oscillating pseudovector fields that cannot be easily probed in existing spin-precession experiments~\cite{Gaul2020a}. Because such effects scale rapidly with $Z$~\cite{Gaul2020b}, YbOCHDT may provide orders-of-magnitude increased sensitivity over previously considered chiral probes. Production and characterization of chiral ytterbium monomethyl, YbCHDT, may be especially promising by bringing the heavy nucleus closer to the chiral center.

\section{Conclusion}
We have observed ytterbium monomethoxide, YbOCH$_3$, and performed initial measurements to assess its feasibility for laser cooling and precision measurements of $CP$-violating electromagnetic moments. The measured vibrational frequencies, vibrational branching ratios, and radiative lifetimes are compatible with a direct laser cooling experiment. Electronic structure calculations have aided the vibrational assignments. We have interpreted some of the vibrational branching ratios using a simple linear vibronic coupling model. Despite having stronger off-diagonal decays than the isoelectronic species YbF and YbOH, the vibrational branching ratios in YbOCH$_3$ converge rapidly enough to permit efficient laser cooling. Even in the absence of full 3D laser cooling, YbOCH$_3$ would allow optical cycling of $\sim$150 photons with about 3 vibrational repumping lasers. This scale enables efficient state preparation, transverse cooling, and/or unit-efficiency readout for precision measurements in a molecular beam.

Further data are needed to sharpen the developing picture of the $\tilde{A}\,\,{}^2 \! E_{1/2} - \tilde{X}\,\,{}^2 \! A_1$ band. The $\tilde{A} \, ^2\Pi$ states in YbF and YbOH are quite complex due to a series of perturbing states nearby in energy~\cite{Mengesha2020, Lim2017}, and it appears that analogous states exist in YbOCH$_3$. Careful studies of the $\tilde{A}$ state and nearby perturbing levels will be necessary to develop a laser cooling scheme capable of scattering $\gg 10^3$ photons with YbOCH$_3$. Detailed understanding of vibronic interactions will also inform the feasibility of laser cooling increasingly complex polyatomic molecules~\cite{Dickerson2020}. Finally, vibronic interactions with matrix elements off-diagonal in $K$ may lead to leakage from the optical cycle at rates high enough to require repumping. Characterizing these effects using high-resolution spectroscopy will also be important prior to laser cooling YbOCH$_3$.

\section{Acknowledgments}
We are grateful to Nick Hutzler, Ivan Kozyryev, and Phelan Yu for constructive feedback on this manuscript. We thank Lan Cheng and the PolyEDM collaboration for many useful discussions. This work was supported by the W. M. Keck Foundation and the Heising-Simons Foundation. BLA acknowledges financial support from the NSF GRFP. ZDL was supported by the Center for Fundamental Physics (Fundamental Physics Grant) and the Templeton Foundation. HS was supported by the Ezoe Memorial Recruit Foundation. The computations in this paper were run on the FASRC Cannon cluster supported by the FAS Division of Science Research Computing Group at Harvard University.