\subsection{Nodal ARPES data and power-law analysis with $k$-independent self-energy}

In Fig. \ref{pll} we show the imaginary part of the nodal self-energy of (Pb,Bi)$_{2}$Sr$_{2-x}$La$_x$CuO$_{6+\delta}$  ~over a large range in frequency, doping and temperature.
Assuming the bare dispersion to be linear, $\varepsilon({k}) = v_{F}({k}-k_{F})$, with $v_{F}$ being the Fermi velocity and $k_{F}$ the Fermi wave vector, the commonly made assumption of negligible $k$ dependence - $\Sigma({k},\omega) \simeq \Sigma(\omega)$ - reduces the single-particle spectral function to a Lorentzian lineshape as a function of momentum ${k}$ at each fixed frequency $\omega$, i.e., $L(k)= \frac{W}{\pi}\frac{\Gamma/2}{(k-k_{*})^{2}+(\Gamma/2)^{2}}$.
Here $W(\omega)$ is the intensity, $k_*(\omega)$ the peak position, and $\Gamma(\omega)$ its width.
For the results presented in the rest of this paper, $W(\omega)$ is not a key parameter and will not be discussed further.
In practice, the peaks are better fit using a Voigt lineshape, the Gaussian part of which accounts for experimental resolution. In this commonly adopted framework, the imaginary part of the self-energy is then proportional to the width of the Lorentzian component as $\Sigma^{''}(\omega) = v_{F}\Gamma(\omega)/2$ \cite{Valla1999}.
\\

\begin{figure}[h!]
\includegraphics[width=17cm]{Figures/PLL.png}
\caption{ {\bf Nodal self-energy of single-CuO$_2$-layer (Pb,Bi)-2201.} {\bf a}-{\bf e}: Temperature-dependent self-energies from ARPES for five different doping levels, extracted using symmetric Voigt fits to MDCs, plotted using colour-coded solid lines for temperature. The dashed, black lines are results to two-dimensional (in $T$ and $\omega$) fits employing three parameters $(\alpha,\beta,\lambda)$, using the power-law-liquid formalism introduced in \cite{Reber2019} and given in Equation 1. Fitting parameters used are indicated in each data panel, and gathered together with the analogous parameters for Bi-2212 \cite{Reber2019} in panel {\bf h}. The red dashed lines indicate the marginal Fermi liquid ($\alpha$ = 0.5 at optimal doping), which is shown above the canonical Presland dome {\bf i} used to determine the doping level of the measured samples. Panel \textbf{g} shows a typical measured ARPES dataset, containing both the positive-$k$ and negative-$k$ nodal branches, along the $k$-space direction indicated in the schematic Fermi-surface in \textbf{f}.  Table S1 in the SI lists all temperatures measured for each doping level.}
\label{pll}
\end{figure}

Confirming and extending the results in \cite{Reber2019}, the ($\omega$,$T$)-dependent self-energy extracted in this way and shown in Fig.\ref{pll}
can be very well described by a remarkably simple phenomenological model dubbed the Power Law Liquid (PLL) \cite{Reber2019}, with three dimensionless parameters $\alpha$, $\beta$ and $\lambda$:
\begin{equation}\label{Reber_pll}
\Gamma = \frac{2\Sigma^{''}(\omega, T)}{v_F} = G_{0}(\omega,T) + \lambda\frac{ [(\hbar \omega)^{2} + (\beta k_{B} T)^{2} ]^{\alpha}}{(\hbar \omega_{N})^{2 \alpha - 1}}~.
\end{equation}
Here $\lambda$ is a coupling constant describing the strength of the interaction, normalized to an energy scale $\hbar \omega_{N} = 0.5$eV for all dopings, and the parameter $\beta$ sets the balance between the relative influence of temperature and frequency.
$G_{0}(\omega,T)$ is an extra term, combining a self-energy contribution from impurity scattering, as well as electron-phonon coupling, described in the SI (S1).
The electron-phonon contribution is most clearly seen around an energy of 70 meV \cite{Bogdanov2000,Cuk2005}, and is prominent because the normal state ARPES data presented here access much lower temperatures than those from Bi-2212 \cite{Reber2019}.
Note that Equation~(\ref{Reber_pll}) describes the marginal Fermi-liquid \cite{Varma1989} with a power of unity ($\alpha = 1/2$) at optimal doping, and the quadratic temperature and frequency behaviour of a Fermi liquid emerges when $\alpha = 1$.

Figs.\ref{pll}b-f show the experimentally extracted self-energies (coloured lines), together with the result of a two-dimensional ($\omega$,$T$) fit to the data using Equation~(\ref{Reber_pll}) (dashed lines).
An overview of the parameters extracted from the fits is given in Fig.~\ref{Reber_pll}g, together with the parameters for Bi-2212 \cite{Reber2019}. Compared to the latter, the ARPES data presented here cover a complementary doping range from near optimal doping ($p=0.14$) to such overdoping that superconductivity disappears.
Across this doping interval, the power increases smoothly from 1.04 ($\alpha=0.52$) to 1.68 ($\alpha= 0.84$), meaning that even on overdoping out of the superconducting dome, the quadratic power of the Fermi-liquid ($\alpha=1$) is not yet reached.

The continuity in $\alpha$ values evident in Fig.~\ref{Reber_pll}g emphasises the remarkably good agreement between the nodal self-energy behaviour in both single- and bi-layer cuprates, also illustrating that these self-energies near $E_F$ are well described by Equation ~\ref{Reber_pll} for a temperature range spanning 200K.

\subsection{Momentum-dependent power-law exponent and asymmetric ARPES MDCs}
The above analysis does hinge on the aforementioned assumption of negligible $k$ dependence: such that the spectral function part of the MDC lineshape can be well fit using a symmetric Lorentzian for each peak \cite{Gaebel2004}.
In Fig. \ref{MDC_fits}a-d we show the Lorentzian fits in detail.
Panel a shows that a Lorentzian spectral function gives a very good fit \emph{close to $E_F$}, and yields a small residual.
However, at higher binding energy, such as at 200 and 300meV, the foot of the peaks clearly show that the data (black crosses) are not captured completely by the two symmetric peaks in the fit (red lines).
The experimental bottom line is that the MDC peaks in Fig. \ref{MDC_fits}b \& \ref{MDC_fits}c are asymmetric: showing more spectral weight at large momenta $|{k}| > |k_{*}(\omega)|$, compared to $|{k}| < |k_{*}(\omega)|$ (leading to residual values of differing signs inside and outside the MDC peak pair).
This behaviour cannot be captured by a fit based on Lorentzians and the power law liquid.
The analysis and discussion covered by Figs. S2-S4 of the SI excludes that the MDC asymmetry comes from the non-zero curvature of an underlying `bare band', from inhomogeneous detector response, or from strongly structured background signals.
A simple yet profound possible explanation of the MDC asymmetry is then a $k$ dependency of the electronic self-energy itself.

\begin{figure}[h!]
\includegraphics[width=16cm]{Figures/MDCfits.png}
\caption
{\textbf{Closer look at MDC lineshapes for UD32 nodal ARPES data.} \textbf{a}-\textbf{c}: a trio of MDCs at energies indicated, with symmetric peak fits in red. The residuals grow as the binding energy grows. Here \textbf{d} shows the resulting $2\Sigma''(\omega)/v_F = \Gamma(\omega)$ from the symmetric fits. \textbf{e}-\textbf{g}: The same three MDCs now fit in purple using our model given by Equation \ref{Mod_lor}, including the holographically-predicted $k$ dependence, with $V$ set to -1.
The residuals are clearly superior to the symmetric fit for energies further from $E_F$. \textbf{h:}  The imaginary part of the self-energy  $2\Sigma^{''}(k, \omega)/v_F = \Gamma_H(k, \omega)$ at $k_*(\omega)$ (red), which includes the $k$-dependent self-energy (blue), and the free fitting parameter $G_0$ (green).}
\label{MDC_fits}
\end{figure}

Inclusion of a $k$-dependent self-energy, which is missing in the PLL ansatz of Equation \eqref{Reber_pll}, provides a significant experimental extension to the data-based fitness tests for aspiring theories for these materials.    One model-independent approach would be to linearly expand the measured self-energy around $k_{F}$, reported in the SI (Fig. S8), which significantly enlarges the set of fitting parameters, by means of an additional frequency and temperature dependent function.
Instead, we can rely on the theoretical input provided by AdS/CFT, and adjust the fitting ansatz without {\em any} expansion in the number of parameters.

\subsection{Semi-holographic theoretical description}

Physically, a semi-holographic model describes an electron that interacts with a CFT (accounting for a quantum critical state deformed by non-zero $T$ and chemical potential) via linear coupling to a fermionic operator ${\cal O}$ that has a unique scaling dimension. As a result, the self-energy becomes proportional to the correlation function $\langle {\cal O O}^\dagger \rangle (k,\omega,T)$ in the CFT. The latter can be determined from the holographic dictionary \cite{Zaanen2015, hartnoll2016holographic}, and automatically inherits certain scaling properties from the `critical' CFT. To be applicable to the cuprate strange metals, the model must have a dynamical critical exponent $z \rightarrow \infty$ in order to recover the power law liquid when $k$ is close to $k_F$.
From this imposition of local quantum criticality, there follows a fundamental condition from the theory side that the scaling exponents \emph{have to be momentum dependent} \cite{Faulkner2011a}.

Here, for the CFT we use an Einstein-Maxwell-Dilaton model of holography, specifically the Gubser-Rocha model \cite{Gubser2010, Gubser2012}, as it offers an analytical treatment of the gravitational spacetime.
In the long-wavelength limit, within the framework of an emergent particle-hole symmetry, our semi-holography model then gives at $T=0$:
\begin{subequations} \label{holo_spectralfunction}
\begin{eqnarray}
    \frac{2 \Sigma^{''}(k, \omega)}{v_F} &=& \lambda\frac{ [(\hbar \omega)^{2} ]^{\alpha(k)}}{(\hbar \omega_{N})^{2 \alpha(k) - 1}}~,\\
    \alpha(k) &=& \alpha \left( 1-\frac{{k}-k_{F}}{k_{F}} \right)~.
\end{eqnarray}
\end{subequations}

The semi-holographic self-energy can also be generalized to non-zero temperature \cite{Faulkner2011a, Gubser2012}, and under our conditions is well approximated by replacing $(\hbar \omega)^2$ with $(\hbar \omega)^2 + (\beta k_B T)^2$ in Equation (\ref{holo_spectralfunction}a).
This also eases comparison to the PLL in Equation (1), and highlights the key new insight that as frequency increases, and $k \simeq k_F - \omega/v_F$ departs from $k_F$, a $k$ dependence emerges in the \emph{exponent} describing the ($\omega,T$)-dependence of the self-energy.
For low frequencies, where $k \simeq k_F$, one returns to Equation~\eqref{Reber_pll} of the PLL.
The behaviour predicted by holography should leave a clear experimental fingerprint, namely that the ARPES MDCs are asymmetric: the data shown in Fig. \ref{MDC_fits} show this is the case.

Fig.~\ref{Holo_spectral}a shows a simulated spectral function, generated using the semi-holographic self-energy shown in Figs.~\ref{Holo_spectral}c \& ~\ref{Holo_spectral}d.
The resulting MDC asymmetry at non-zero energy is illustrated in Fig.~\ref{Holo_spectral}b, visible as an increased intensity at the $|k|>|k_{*}|$ side of the peak maximum, just like in the experimental data of Fig. \ref{MDC_fits}c.
The next step is to test whether this holographic approach can yield superior quantitative fits to the ARPES data.

\begin{figure}[h!]
\includegraphics[width=9cm]{Figures/Sim_Kdep_rev.png}
\caption{{\bf Simulated spectral function at optimal doping, using the $k$-dependent self-energy from semi-holography.} \textbf{a}: Spectral function with a linear bare band (v$_F=$ 4 eV\AA, k$_{F}$=0.45\AA), convoluted with the experimental $k$ resolution and \textbf{b} MDC at non-zero binding energy generated using the full $\omega$ and $k$-dependent self-energy plotted in panels \textbf{c} and \textbf{d} from the semi-holographic model at $T$=20 K. Note that \textbf{c} includes the imaginary part of a $k$-independent phonon self-energy.}
\label{Holo_spectral}
\end{figure}

\subsection{Fitting the semi-holographic theoretical model to the ARPES data}

As non-Lorentzian MDCs force a major departure from a cornerstone of ARPES data analysis methodology, we now describe how to deal with this in the analysis of real data from (Pb,Bi)$_{2}$Sr$_{2-x}$La$_x$CuO$_{6+\delta}$ .
The non-zero-$T$ version of the self-energy with the $k$-dependence in Equation (\ref{holo_spectralfunction}b), suggests a modified fit-function $L_{H}(k)$ at the fixed $\omega$ and $T$ relevant for each MDC:
\begin{subequations} \label{Mod_lor}
\begin{eqnarray}
    L_{H}(k) &=& \frac{W}{\pi} \frac{\frac{\Gamma_H}{2}}{(k-k_{*})^2 + \left(\frac{\Gamma_H}{2}\right)^2}~,\\
    \Gamma_H(k) &=& G_{0} + \lambda \frac{[(\hbar\omega)^2+(\beta k_{B}T)^2]^{{\alpha({k})}}}{\hbar\omega_N^{2{\alpha({k})}-1}}~,\\
    \alpha({k})&=& \alpha \left(1 +V \left[\frac{{k}-k_{F}}{k_{F}}\right] \right)~.
\end{eqnarray}
\end{subequations}

$\Gamma_{H}(k)$ in Equation (\ref{Mod_lor}b) captures both the peak width and its asymmetry via momentum dependence built into the exponent $\alpha({k})$ in Equation (\ref{Mod_lor}c).
By fixing $\alpha$, $\beta$, $\lambda$ and $\omega_N$ to the PLL values at low energies, only $G_{0}$ and the asymmetry parameter $V$ remain free to vary in the fitting process for each MDC.

The Gubser-Rocha holographic model used here, however, actually $\emph{requires}$ that $V=-1$ for all frequencies.
Thus, in this case the resulting fit function has exactly the same number of free parameters as the PLL, and is used to fit the experimental MDCs as shown using the purple lines in Fig. \ref{MDC_fits}e-g.

Comparing Figs. \ref{MDC_fits}b and \ref{MDC_fits}c (PLL, $k$-independent) with Figs. \ref{MDC_fits}f and \ref{MDC_fits}g ($k$-dependent scaling exponents), it is clear that on adopting Equation \eqref{Mod_lor}, the residuals in the panels f and g for energies well below $E_F$ drop, becoming as low as they were for the PLL at $E_F$.
In particular, the zoom to Figs. \ref{MDC_fits}c and \ref{MDC_fits}g illustrates clearly that the asymmetry of the MDC peaks is now captured almost perfectly.
Fig. \ref{MDC_fits}h shows the total self-energy along the loci of the MDC peak-maxima in red, with the holographic, k-dependent part in blue, and the free fitting parameter G$_0 (\omega,T)$ in green. The latter can be seen to automatically take on the combined form of an offset (impurity scattering) plus a step function centered at the phonon energy of 70 meV.
\vspace{5mm}

\subsection{Testing the semi-holographic prediction}

Upon determining all the free parameters in the theory (see Methods), two experimental tests of the semi-holographic results can then be performed. The first is a mild one, unconnected to the observed MDC asymmetry, and concerns the value of $\beta$, shown in Fig. \ref{doping_V}a.
The experimental values (blue squares) all reside in an interval between $\beta =3-4$ across the doping levels studied.
Those coming from our model lie between $\beta =2-3$, and show a doping-dependent (generally upward) trend very similar to that of the experimental values.
A detailed derivation of parameters such as $\beta$ falls outside the experimental focus of this paper, and will be presented in a separate publication.

\begin{figure}[h!]
\includegraphics[width=14cm]{Figures/V_Beta.png}
\caption{{\bf Testing the predictions from semi-holography}
\textbf{a}: Comparison between the coefficient $\beta$, determining the amount of temperature dependence in the self-energy, obtained from the PLL fit to the ARPES data (blue squares) and predicted by the holographic model (black line), as a function of the doping-dependent power $\alpha$, which is translated into doping on the upper x-axis. {\bf b:} Asymmetry parameter $V$ of the momentum-dependent scaling exponent extracted from the experimental ARPES data over the full doping range measured at 8K, obtained by performing MDC fits to the data using Equation \eqref{Mod_lor}, where the leading scaling exponent $\alpha$, the coupling constant $\lambda$ and the quantity $(\hbar\omega)^{2} + (\beta k_{B} T)^{2}$ are fixed. A value $V = 0$ means there is no asymmetry, and $V = -1$ is the prediction from the semi-holographic Gubser-Rocha model. The high temperature fits are shown in the SI section S6.}
\label{doping_V}
\end{figure}

The second, more stringent test is given by the experimentally observed MDC asymmetry.
The non-zero value of $V$ is a non-trivial result, which is built into holography, but must be accounted for by other theoretical approaches aiming to describe the spectral function of the nodal charge carriers in cuprates. For the Gubser-Rocha model adopted here, the asymmetry parameter should have a frequency-independent value of $V=-1$.
We test this by performing a second fit of the ARPES data, now leaving $V$ as a free parameter for each $\omega$.
Fig. \ref{doping_V}b shows that - with no guiding restraint at all applied to $V(\omega)$ - the experimental data yield for all doping levels, temperatures (see SI) and frequencies a value of $V \approx -1$, close to the semi-holographic prediction.

Fig. \ref{doping_V} shows that there is still some room for improvement as regards to both $\beta$ and $V$ in the current approach.
Thus, these ARPES data and their parametrization presented here are an invitation to the many existing $z=\infty$ holographic models besides the analytical Gubser-Rocha model - or to completely different theoretical methods - to take the next step in capturing the observed $\beta$ values and modest frequency dependence of $V$ seen in experiment.

Before providing a discussion of other contexts in which k-dependent self energies have been discussed, we re-iterate the main results of this paper.
The $\omega$ and $T$ dependence of the electronic self-energy of nodal carriers in the normal state of the single-layer cuprate (Pb,Bi)-2201 is shown in ARPES experiments to follow a single power law with a $k$-\textit{dependent scaling exponent}, yielding asymmetric, non-Lorentzian MDCs for energies away from $E_F$.
From fits to the ARPES data from optimal doping across the whole overdoped portion of the phase diagram that work well over wide ranges in energy and temperature of more then 200 meV and 250 K, respectively, at $k=k_F$ the $k$-dependent power-law exponent is found to take a nearly marginal Fermi-liquid value of unity ($\alpha = 0.51$) at optimal doping, growing to 1.68 ($\alpha = 0.84$) for the most overdoped, non-superconducting system studied. For $k=k_F$, these results connect smoothly to those of \cite{Reber2019}.
However, for values $k\neq k_F$ these powers change.
Such a $k$ dependence is a key prediction from (semi)holographic models of strange metals as a quantum critical phase.
The form of the $k$ dependence and the asymmetry parameter describing its magnitude are all successfully fitted using semi-holographic single-fermion spectral functions.
Importantly, any competing theoretical model for strange metals must be able to account for the observed ARPES MDC peak-asymmetry that is well described here by a momentum-dependent power-law self-energy.

There are few studies of the $k$ dependence of the self-energy in solids.
In the cuprate context, \cite{Randeria2004} and \cite{Veenstra2011} present theory/modelling of possible $k$ dependence in the self-energy resulting in a change in the dispersion velocity or arising from strong electron-phonon coupling, respectively.
In a different condensed-matter realisation, $k$-dependent power laws have recently been observed in the zero-bias conduction anomaly in transport spectroscopy of nanowires \cite{Jin2019}, and linked to non-linear Luttinger liquid theory \cite{Imambekov2009}.
One might therefore wonder if there is a link between the 1D physics in the nanowire case, and the carriers in the nodal $k$-space direction of the quasi-2D (Pb,Bi)-2201 system that seem well described by the 1D CFT encoded semi-holographically in the AdS$_2$ infrared geometry of the $z=\infty$ models.

The results presented here herald future comparisons of the experimentally-tested, theoretical spectral function to other experimental probes such as optical and DC conductivities \cite{Jacobs2015}. Alternatively, many-body condensed-matter theories can be guided by both the experimental self-energies presented here and the results from the AdS/CFT analogue, opening up possible pathways out of phenomenology and into a microscopic understanding of strange metals.

\section{Methods}

\subsection{ARPES measurements}

All ARPES data presented here were recorded at beamline I05 of Diamond Light Source, using (horizontally) linearly polarized light at a photon energy of 28 eV.
Due to this chosen photon energy we can combine excellent energy resolution (12meV) with enough range in $k$ space to measure both the negative-$k$ and positive-${k}$ nodal branches in a single $I({k},\omega)$ image (see Fig. \ref{pll}A), which is very helpful in distinguishing the subtle effects at play from possible complicating factors such as signal background.
The energy resolution, together with the Fermi energy position - was confirmed by means of reference data from an amorphous Au film held in electrical contact with the sample.
To account for resolution in the angular direction, both the Lorentzian MDC fit function used for the PLL analysis, and the  L$_H$ MDC fit function from Equation\eqref{Mod_lor} are convolved with a Gaussian of width 0.01 \AA$^{-1}$. In the symmetric PLL case this yields a Voigt lineshape as indicated in the main text and in Fig. \ref{MDC_fits}.

All data were recorded in swept mode, to ensure any detector-response inhomogeneity is averaged out in the energy direction, and were also confirmed to be free of detector non-linearity effects \cite{Reber2014}. All data were recorded in the X$\Gamma$X direction in the Brillouin zone.

High quality single crystals of diffraction replica-free (Pb,Bi)$_{2}$Sr$_{2-x}$La$_x$CuO$_{6+\delta}$ , or (Pb,Bi)-2201), were grown using floating-zone techniques.
Individual crystals were annealed in varying atmospheres for varying lengths of time, so as to change the oxygen content and with it the hole doping, controlling the carrier concentration and $T_c$.
The critical temperatures were determined either via resistivity or AC-susceptibility measurements and the doping level was read off using $T_{c}$ and the Presland formula \cite{Presland1991, ando2000}. For the sample OD0K, $T_c\lesssim2K$ $\lesssim0.05 T_c^{max}$.

\vspace{5mm}
\subsection{Holography and semi-holography}
The Anti-de Sitter/Conformal Field Theory correspondence (AdS/CFT) links a classical gravitational theory with one additional spatial dimension to a strongly interacting quantum theory on the boundary of the extended space. In practice, it provides a systematic way to compute the response functions for a quantum theory with strong interactions, even at non-zero temperature and in the presence of a chemical potential, by means of solving a (relatively speaking) simpler problem, namely solving classical gravitational equations.
For our purposes, therefore, the curved space is three dimensional and the quantum theory should describe the correlated physics of the maximally entangled fermions in the quasi-2D high-$T_c$ cuprates \cite{Zaanen2015, hartnoll2016holographic}. To be precise, holography describes the physics of a composite fermion $\cal O$ with a large number of degrees of freedom, i.e., in the so-called large-$N$ limit. However, as we here are interested in testing theoretical predictions against spectral functions obtained through ARPES data, we ultimately want to compute the response function for an elementary fermion $\psi$. The solution is found within the semi-holographic framework \cite{Faulkner:2010tq}, where the fermion $\psi$ is linearly coupled to the composite fermion $\cal O$ and the Green's function describing the dynamics of $\psi$ then acquires the form
\begin{equation}\label{green_semi}
  G_{\psi\psi}(\omega, k) = \frac{1}{\omega - v_F (k-k_F) - g_s G^{-1}(\omega, k)} \text{ ,}
\end{equation}
where $G^{-1}(\omega, k)$, assuming the role of the self-energy, is the two-point function of the composite fermionic operator $\cal O$, and $G_{\psi\psi}$ can be shown to be properly normalized as the Green's function of an elementary fermion \cite{Gursoy2012}.
The most important piece of the puzzle, the qualitative form of the complex self-energy at energies near the Fermi level, thus comes from  holography. It is in this context that we compute the inverse Green's function $G^{-1}(\omega, k)$ appearing in Eq. \eqref{green_semi}. This is done by solving the Dirac equation on the curved gravitational background spacetime \cite{Iqbal2011}.

In general, as it has been done for the model in this paper, the gravitational equations needed to compute the holographic Green's function can only be solved numerically. Nonetheless, the qualitative behavior of the low-energy Green's function of such a fermionic operator can be obtained analytically and it is completely determined by the infrared properties of the theory, captured in the duality by the inner geometry of the gravitational spacetime.
The electron self-energy obtained in this a way is the one we use for the nodal holes near the Fermi surface within the framework of an emergent particle-hole symmetry, and it possesses the desired ($\omega,T$)-scaling properties, i.e., power-law exponents matching the experimental ARPES data.
In particular, this implies that in the class of Einstein-Maxwell-Dilaton geometries characterized by the dynamical critical exponent $z = \infty$, and in the presence of a Fermi surface, the inverse Green's function at $T = 0$ (the low-energy behavior can be easily generalized to small non-zero temperature \cite{Gubser2012, Faulkner2011}) has to assume a low-energy form near $k_F$ \cite{hartnoll2016holographic} given by
\begin{equation}\label{green_composite}
  G^{-1}(\omega, k) = g \omega (-\omega^2)^{\nu_k-1/2} + \dots \text{ ,}
\end{equation}
where $g$ is constant in the limit $\omega, k-k_F \rightarrow 0$.
It is important to note that the $k$-dependent power $\nu_k$ is solely dictated by the IR of the theory and thus is independent of any of the UV details. The lifetime of the fermionic excitations near the Fermi surface are then described by
\begin{equation}
  \Sigma^{''}(\omega, k) = -\frac{\lambda v_F}{2} (\omega^2)^{\nu_k} + \dots~.
\end{equation}
Theoretically, the constant $\alpha \equiv \nu_{k_F}$ in the exponent in Eq. 2b, depends on the two (dimensionless) parameters in the Dirac equation, the mass $m$ and the charge $q$, that encode certain defining properties of the conformal field theory \cite{hartnoll2016holographic}. We keep the mass fixed close to the limit of $m = -1/2$ of the range allowed for by semi-holography ($-1/2 < m < 1/2$) \cite{Gursoy2012}. Then, varying only $q$ as a function of the hole doping $p$, captures the doping dependence of the exponent $\alpha(p)$ seen in the experimental data, using the numerically obtained linear relationship $\alpha \approx 1.9 q$.
Finally, the value of $\lambda$ (see Eq. 2a) leading to a match with experiment at each doping level can be obtained by adjusting the strength of the coupling $g_s$ between the fermion $\psi$ and the conformal field theory.

We can see that Eq. \eqref{green_composite} has the same qualitative form as the semi-holographic response \eqref{green_semi} for $\omega, k-k_F \rightarrow 0$.
Semi-holography thus only changes the UV of the theory, i.e. the constant $g$, but it does not change the emergent behavior $\Sigma'' \propto (\omega^2)^{\nu_k}$. In this sense, while we showed that in our `bottom-up' approach we possess enough tuning knobs so as to get a semi-holographic spectral function that closely resembles the one observed in the ARPES experiments, we do not know whether the values used do indeed correspond to a consistent theory of gravity.
Nonetheless, our main point that the self-energy contains a momentum-dependent power is a universal feature of $z = \infty$ holographic calculations, independent of the particular UV completion applied.
In fact, while the theoretical model proposed here is rooted in holography, any quantum critical theory with infinite dynamical exponent and a large-$N$ suppression of higher-order correlation functions will be able to give rise to an effective correlation function of the same form as in Eq. \eqref{green_semi} \cite{hartnoll2016holographic}.

\newpage


\section{Introduction}

In the last years state-of-art results in natural language processing (NLP) are often obtained with Transformer-like architectures based on the self-attention mechanism \cite{vaswani2017attention} such as BERT \citep{devlin-etal-2019-bert}, GPT-3 \cite{brown2020language}, T5 \cite{raffel2020exploring}, which could have billions of parameters. Due to many parameters, these architectures require lots of time and hardware resources to be trained.

Curriculum learning (CL) is one of the popular methods to reduce training time and increase the resulting quality of the model. Inspired by the importance of adequately ordering information when teaching humans \cite{avrahami1997teaching}, curriculum learning increases the difficulty of training samples shown to the model over time \cite{elman1993learning}. Previous studies have demonstrated that curriculum learning significantly impacts training time and quality in different machine learning domains, such as computer vision \cite{Soviany2020CurriculumLW} and reinforcement learning \cite{narvekar2020curriculum}. In NLP, some results hint that CL might be beneficial \cite{platanios-etal-2019-competence,xu2020curriculum,kocmi-bojar-2017-curriculum}; however, these results are not as optimistic as in reinforcement learning setup.

We suggest dividing recent research in curriculum learning into two main categories: {\em task-driven} curriculum and {\em data-driven} curriculum. The idea of the task-driven curriculum was inspired by human behavior. First, the model learns how to solve a simple task, and then the difficulty is gradually increased. This type of curriculum proposed by \citet{bengio2009curriculum} is considered to be classical, and a majority of curriculum-related results are obtained in this framework. Alternatively to the task-driven curriculum, some curricula try to use some form of filtering or sorting of training data that could facilitate learning a model on a given task. We suggest calling these curricula {\em data-driven} and distinguishing them from the classical task-based approach.

This paper attempts to understand when data-driven curriculum learning works for transformer-based language models. Generally, data-driven curriculum learning is organized in two steps: first, estimating the complexity for the elements that comprise the dataset; second, designing a  sampling strategy, thus forming a curriculum. In the first part of the paper, we list potentially useful natural language processing complexity measures. The second part discusses possible sampling strategies that might apply to corresponding complexity measures. We run extensive experiments with different metrics and sampling strategies on three classes of NLP tasks: unsupervised learning with masked language modeling, text classification, and machine translation. Our experiments show that data-driven curriculum learning does not give quality increase or time reduction on all metric-sampling strategy setups and often makes results even worse.

\begin{figure*}[h!]
	\centering
	\begin{subfigure}{.4\textwidth}
		\centering
    	\includegraphics[scale=0.22]{s140_SM_res.png}
    	\caption{Sentiment140 with sort-merge sampler for all complexity measures.}
    	\label{fig:s140_SM_all_metrics}
	\end{subfigure}
	\begin{subfigure}{.4\textwidth}
		\centering
    	\includegraphics[scale=0.22]{s140_max_wf_rank.png}
    	\caption{Sentiment140 with max word rank complexity measure for all samplers.}
    	\label{fig:s140_max_wf_rank_all_samplers}
	\end{subfigure}\\
	\begin{subfigure}{.4\textwidth}
		\centering
    	\includegraphics[scale=0.22]{hnd_SS_res.png}
    	\caption{Hyperpartisan News with sort-shuffle samples for all complexity measures}
    	\label{fig:hnd_SS_all_metrics}
	\end{subfigure}
	\begin{subfigure}{.4\textwidth}
		\centering
        \includegraphics[scale=0.22]{hnd_max_wf_rank.png}
        \caption{Hyperpartisan News with max word rank complexity measure for all samplers.}
        \label{fig:hnd_max_wf_rank_all_samplers}
	\end{subfigure}
	\caption{Pre-trained BERT fine-tuned on Sentiment140 and Hyperpartisan News Detection datasets. Accuracy of the classifier as a function of the number of training steps.}
	\label{fig:text_classification}
\end{figure*}

\section{Metrics}

The first important part of the curriculum learning pipeline is measuring the complexity of samples for a given dataset. Texts could have a complex structure, and one can measure their complexity in different ways. A variety of heuristically motivated methods is accompanied by several metrics based on specific aspects of information theory. For a review of heuristic text complexity measures such as length of TF-IDF \cite{aizawa2003information} we address the reader to Appendix \ref{a:heu}. In this paper, we also explore the metrics initially proposed by \citet{ay2006unifying} to measure the complexity of finite systems and try to see if one could apply these metrics to NLP tasks.

\citet{ay2006unifying} observes that for finite systems, a set of parts impacts the complexity of the system as well as inter-dependencies of the parts. In the context of NLP, this means that text is more than just a bag of words. The authors propose four different metrics to estimate the complexity of a system. However, one of these metrics maximizes on single-letter texts, such as "Aaaaaaaaa," while the second was created to measure cyclic sequences and does not apply to texts. Thus we experiment with two other metrics, namely, Tononi, Sporns, and Edelman (TSE) \cite{tononi1994measure} and excess entropy (EE), and adapt them to the complexity of texts. For the calculation of TSE and EE for NLP we address the reader to Appendix \ref{a:cal}.

\section{Samplers}

The second important part of curriculum learning is the sampling strategy (or sampler) - the algorithm deciding which samples should be shown to the model at which moment. Let us observe existing curricula and suggest some new ones.

\textbf{Competence-based. CB}\\
A competence-based curriculum, offered by \citet{platanios-etal-2019-competence}, uniformly samples data from increasing dataset's prefix. Competence is a function $c(t)$, which defines the size of the dataset prefix.
$$c(t) = \min \left( 1, \sqrt{t\frac{1 - c_0^2}{T} + c_0^2}\right)$$
Where $T$ - total number of steps, $t$ - current step, $c_0$ - hyperparameter set to $0.01$.

\textbf{Hyperbolic. HYP} \\
The main idea of this sampler is to increase average batch complexity through time. All samples are split by complexity into $N$ sequential buckets with equal size. Training time is divided into $N$ epochs and the probability of sampling the element from the $j$-th bucket on the $i$-th epoch is proportional to the distance between $j$ and $i$.
$$Pr_i(j) = \frac{c}{|j - i|^{0.5}}$$
Where $Pr_i(j)$ - probability to sample from $j$-th bucket on the $i$-th epoch, $c$ - constant to guarantee that sum of all probabilities equals to $1$.

\textbf{Difficulty-based. DB} \\
This sampler is a reversed version of the competence-based one. A difficulty-based sampler takes elements from a linearly decreasing suffix instead of sampling from a gradually increasing prefix.

\textbf{Sort-shuffle. SS}\\
All previously described samplers do not guarantee that the model would see each element in the training data. Sort-shuffle samples each element exactly once, randomly splitting the data into batches and sorting by average complexity.

\textbf{Sort-merge. SM}\\
Many complexity estimates correlate with the length of the text. The main idea of a sort-merge sampler is to remove this correlation and train the model on stable length distribution. This algorithm consists of four main steps: sort dataset by length; sequentially split into buckets; sort each bucket by a complexity metric; form $i$-th batch from $i$-th elements from each bucket. Like a sequential one, the sort-merge sampler shows each element to the model exactly once.

Equipped with the list of metrics and curriculum samplers, we can discuss our experimental results.

\section{Experiments}
We perform our experiments on three NLP tasks: text classification, machine translation (NMT), and masked language modeling (MLM). Here we discuss the first task of classification in detail. The extensive results of the experiments are available in Appendix \ref{a:ex}. All the experiments are performed with the HuggingFace library \cite{wolf-etal-2020-transformers}, which provides the models with their setups, such as hyperparameters and tokenizers. We did not change default parameters in our experiment unless specifically stated otherwise. Thus, the dataset and the model specify every experiment. We use the base version of the BERT model \cite{devlin-etal-2019-bert} for MLM and classification, and the small version of the T5 model \cite{raffel2020exploring} for machine translation. Experiments were performed on BooksCorpus\footnote{\url{https://huggingface.co/datasets/bookcorpus}} dataset for MLM, Sentiment140\footnote{\url{https://www.kaggle.com/kazanova/sentiment140}} and Hyperpartisan News Detection\footnote{\url{https://huggingface.co/datasets/hyperpartisan_news_detection}} for classification, and WMT16-en-de\footnote{\url{https://huggingface.co/datasets/wmt16}} for machine translation. To estimate the curriculum's convergence speed, we calculate the average number of steps to reach a threshold that is 10\% lower than the resulting saturation quality metric for every problem.

\subsection{Text Classification}

 Figure~\ref{fig:text_classification} summarizes the experiments with BERT for text classification. Neither different samplers nor complexity measures improve a BERT-based classifier's resulting accuracy.

\subsection{Masked Language Modelling}

 Figure~\ref{fig:BooksCorpus_graphs} shows the results of MLM pretraining of BERT on BooksCorpus. Irrespective of sampling, the complexity measures have similar ranking in terms of their performance on MLM: length, likelihood, TSE, EE, TF-IDF, maximum word rank. Since sorted sampler takes length into account by design, it is not included in the corresponding plots. Data-based curricula show inferior results in comparison with the baseline.

 \begin{figure*}[t]
	\centering
	\begin{subfigure}{.25\textwidth}
		\centering
		\includegraphics[scale=0.2]{BooksCorpus_CB}
		\caption{Competence-based}
		\label{fig:BooksCorpus_graphs_CB}
	\end{subfigure}
	\begin{subfigure}{.25\textwidth}
		\centering
		\includegraphics[scale=0.2]{BooksCorpus_DB}
		\caption{Difficulty-based}
		\label{fig:BooksCorpus_graphs_DB}
	\end{subfigure}
	\begin{subfigure}{0.25\textwidth}
		\centering
		\includegraphics[scale=0.2]{BooksCorpus_Hyp}
		\caption{Hyperbolic}
		\label{fig:BooksCorpus_graphs_Hyp}
	\end{subfigure}\\
	\begin{subfigure}{.25\textwidth}
		\centering
		\includegraphics[scale=0.2]{BooksCorpus_SS}
		\caption{Sort-shuffle}
		\label{fig:BooksCorpus_graphs_SS}
	\end{subfigure}
	\begin{subfigure}{.25\textwidth}
		\centering
		\includegraphics[scale=0.2]{BooksCorpus_SM}
		\caption{Sort-merge}
		\label{fig:BooksCorpus_graphs_SM}
	\end{subfigure}
	\begin{subfigure}{.25\textwidth}
		\centering
		\includegraphics[scale=0.5]{BooksCorpus_legend}

		\label{fig:BooksCorpus_graphs_legend}
	\end{subfigure}
	\caption{Loss function dependency on the number of training steps on MLM for BooksCorpus dataset during the first 40k steps of training. Every plot depicts results for six different complexity estimates combined with a specific sampler.}
	\label{fig:BooksCorpus_graphs}
\end{figure*}

\subsection{Neural Machine Translation}

Table~\ref{table:MT_quality} shows the experiments with T5 model \cite{raffel2020exploring} for machine translation and various curricula. We use the BLEU metric to estimate the quality of the resulting models. We calculate the average BLEU score over ten validations at saturation. Once again, curriculum learning does not give any notable benefits.

\begin{table}[h!]
    \caption{The average BLEU score from 50k to 100k steps on WMT16 dataset. Results better than the baseline are highlighted. '-' denotes the cases when complexity measure and sampler are not compatible.}
	\label{table:MT_quality}
	\centering
	\begin{tabular}{lcccccc}
		\hline
		Metrics & & \multicolumn{3}{c}{Samplers} &\\
		\hline
		& CB & DB & Hyp & SS & SM & \\
		\hline
		baseline & & \multicolumn{3}{c}{18.3} & \\
		\hline
		length &  10.1 & 17.4 & 16.3 & - & -\\
		TSE &  10.3 & \bf 18.4 &  16.8 &  13.8 &  14.8\\
		EE &  10.2 &  18.2 &  16.9 &  13.3 &  15.0\\
		\hline
	\end{tabular}
\end{table}

\section{Discussion}

We try to interpret obtained results cautiously. Though \citet{platanios-etal-2019-competence} report that competence-based sampling is beneficial for recurrent neural networks, we could not reproduce this result in transformer-based architectures. We also run experiments to check whether data-based curricula could work on non-transformer architectures. The results do not look encouraging; see Appendix \ref{a:plat}.

Curriculum learning depends on subtle factors, for example, a correct choice of hyperparameters. It is hard to check all possible values of hyperparameters, yet to the best of our capabilities, we address this issue in Appendix \ref{a:hyp}. The results do not seem to depend on the learning rate, and once again, curriculum learning shows no benefits.

At this point, we can only conclusively say two things: (1) a deeper investigation of the underlying information theoretic principles that stand behind curriculum learning is badly needed; (2) until we better understand these principles, data-based curriculum learning is a gamble with very low odds to gain either speed or resulting performance.

\section{Conclusion}
In this work, we ran extensive experiments with curriculum learning for transformer-based architectures on three NLP tasks: masked language modeling, text classification, and machine translation. We demonstrate that curricula do not help in the standard training setting and sometimes even worsen results.

\section{Acknowledgments}
The publication was supported by the grant for research centers in the field of AI provided by the Analytical Center for the Government of the Russian Federation (ACRF) in accordance with the agreement on the provision of subsidies (identifier of the agreement 000000D730321P5Q0002) and the agreement with HSE University  No. 70-2021-00139. This research was supported in part through computational resources of HPC facilities at HSE University \cite{kostenetskiy2021hpc}


\section{Introduction}
Recently, magnesium monofluoride (MgF) has been identified as a promising candidate molecule for laser cooling and magneto-optical trapping experiments \cite{Kang2015, Xu2016, Chae2021}. Compared to other group II monofluoride molecules that have been laser-cooled so far\cite{Barry2014,Truppe2017,Anderegg2017a}, MgF is lighter and has a stronger optical cycling transition in the ultraviolet. These properties allow for exerting a large radiation force to rapidly slow and cool the molecules and produce a magneto-optical trap with a high capture velocity. The predicted low off-diagonal vibrational branching and the simple hyperfine structure of MgF reduces the complexity of the optical setup significantly.

Rotationally resolved optical spectra of MgF have been recorded in absorption \cite{Jenkins1934, Fowler1941, Barrow1967, Walker1968,Novikov1971} and emission \cite{Novikov1971}. The vibration-rotation emission in the electronic ground state has also been studied \cite{Barber1995}. Precise hyperfine and rotational constants of MgF in its X$^2\Sigma^+$  electronic ground state for vibrational states $v=0,1,2$ and 3 were determined from its millimeter-wave spectrum \cite{Anderson1994,Anderson1994a}. Recently, Xu et al. \cite{Xu2019} recorded optical absorption spectra of the A$^2\Pi \leftarrow $X$^2\Sigma^+$  transition, resolving a prior debate in the literature about the sign and value of the spin-orbit coupling constant. However, this study suffered from a large systematic frequency offset of about \SI{-4.1}{\giga\hertz} and the Hamiltonian that was used to fit to the experimental data did not account for the presence of $\Lambda$-doubling and hyperfine structure (\textit{vide infra}) in the excited state\footnote[2]{A further study from the same group \cite{Gu2021} was submitted during the preparation of this article. We find a systematic frequency offset of about $+$\SI{2.3}{GHz} in the reported line centers.}. Recently, optical cycling experiments have been performed \cite{Xia2021}, a first step towards laser cooling experiments.

MgF has also been studied theoretically using ab initio methods \cite{Walker1970, Kang2015, Xu2016}. Pelegrini \textit{et al.} \cite{Pelegrini2005} calculated various properties of MgF as part of a wider study of group II monofluorides. These predictions show good agreement with the available experimental data for CaF. For MgF, they predicted a radiative lifetime of \SI{7.16}{\nano\second} for the A$^2\Pi$ ,$v'=0$ level, and a decay probability of 1.4\% to the X$^2\Sigma^+$ , $v''=1$ first vibrationally excited state. They also predicted the permanent electric dipole moments for the ground and excited states to be \SI{2.67}{D} and \SI{4.23}{D} respectively.

Here, we present hyperfine-resolved UV laser-induced fluorescence (LIF) spectra of MgF produced in a cryogenic buffer gas molecular beam. The large frequency calibration error present in the previous study by Xu \textit{et al.}\cite{Xu2019} is corrected by calibrating our wavemeter with known transition frequencies in Yb. The eigenvalues of an appropriate Hamiltonian are fitted to the measured hyperfine energy levels to derive precise spectroscopic constants for the A$^2\Pi$  state. We record and analyse transition isotope shifts between the two bosonic isotopologues and compare to predictions from mass scaling arguments. The spectral width of isolated lines is measured with high accuracy to determine the radiative lifetime of the A$^2\Pi$  state. The electric dipole moments of the ground and excited states are deduced from the Stark shifting of individual rotational lines in electric fields of up to \SI{10.6}{\kilo\volt\per\centi\meter}. We then determine how opposite parity levels in the excited state mix in an external electric field. This effect can result in large losses from the optical cycling scheme, if stray electric fields are not well-controlled.

The group II metal monofluorides are interesting candidates for laser cooling because of the single, unpaired, and metal-centered electron that is polarized away from the fluorine atom. The internuclear distance and potential energy curves of the ground and first excited states are very similar. This leads to a very diagonal Franck-Condon matrix which reduces the number of vibrational repump lasers required to scatter a large number of photons. Measurements of the fine and hyperfine structure and the electric dipole moments provide information about the spin density at the fluorine nucleus and the charge distribution. This provides information to better understand the bonding structure in these molecules \cite{Kaendler1989}. We compare our results to the other group II monofluoride molecules, CaF, SrF and BaF, which have been studied in detail.

\section{Hamiltonian}
We use the following effective Hamiltonian, which operates in a given vibrational state with energy $E_0$:
\begin{equation}\label{eq:hamil}
    \begin{aligned}[b]
        H&=B\mathbf{N^2}-D\mathbf{N^4}+A \mathbf{L}\cdot \mathbf{S}+\gamma\mathbf{N}\cdot\mathbf{S}\\
        &-\frac{1}{2}p(N_+S_++N_-S_-)+\frac{1}{2}q(N_+^2+N_-^2)\\
        &+a L_zI_z+b_F\mathbf{S}\cdot\mathbf{I}+\frac{1}{3}c(3S_zI_z-\mathbf{S}\cdot\mathbf{I})\\
        &-\frac{1}{2}d(S_+I_++S_-I_-)
    \end{aligned}
\end{equation}
 It describes rotation ($B$, $D$), spin-orbit $A$, spin-rotation $\gamma$, $\Lambda$-doubling ($p$, $q$) and magnetic hyperfine interactions ($a$, $b_F$, $c$ and $d$). The ground state of MgF is a X$^2\Sigma^+$  state, for which $\Lambda=0$ and therefore $A=p=q=a=d=0$.

 Fig. \ref{fig:biglevelscheme} shows the relevant levels and transitions in absence of hyperfine structure. For a given value of $J$, the level with the lowest energy is labeled $\mathcal{F}_1$. We use $\Delta J_{\mathcal{F'}\mathcal{F}''}(N'')$ to label the transitions. The energies of the $J$-levels in the ground state can be calculated using,
 \begin{equation}
    \begin{aligned}[b]
    E(N)&=E_0+BN(N+1)+\gamma N/2 &  \textrm{for}    &&  J&=N+1/2\\
    E(N)&=E_0+BN(N+1)-\gamma (N+1)/2 & \textrm{for} &&  J&=N-1/2 ,
    \end{aligned}
    \label{eq:XstateEnergies}
\end{equation}
and in the A$^2\Pi$ state using the following formula,
 \begin{equation}
    \begin{aligned}[b]
        E(J)&=E_0-\frac{1}{2}\gamma+B\left(J(J+1
        -\frac{3}{4}\right)\\
        &\mp\sqrt{\left(\frac{1}{2}(A+\gamma)-B\right)^2+
    \left(B-\frac{1}{2}\gamma\right)^2\left(J(J+1)-\frac{3}{4}\right)} \hspace{0.2cm}.
    \end{aligned}
    \label{eq:AstateEnergies}
\end{equation}
Here, the minus (plus) sign applies for the $\mathcal{F}_1(\mathcal{F}_2)$ levels respectively.
In the case of $J'=1/2$, expression \eqref{eq:AstateEnergies} reduces to
\begin{equation}
    E\left(J'=\frac{1}{2}\right)=E_0-\frac{1}{2}\gamma-\frac{1}{2}(A+\gamma)
    B.
\end{equation}

\noindent When comparing transition frequencies of different isotopologues, an explicit form for the energy $E_0$ is required. Here we assume that

\begin{equation}
    E_0 = T_e + \omega_e(v+1/2) - \omega_ex_e(v+1/2)^2,
    \label{eqn:E0definition}
\end{equation}

\noindent with $T_e$ being the potential energy minimum of an electronic state, and the remaining terms describing the vibrational energy up to second order in the vibrational quantum number $v$.

\begin{figure}
    \centering
    \includegraphics[]{FIG1_Levelscheme_review.eps}
    \caption{Schematic energy level diagram for a regular $^2\Pi -^2\Sigma^+$ system, where $p$ and $q$ are both positive. Positive parity states are shown in red while negative parity states are shown blue (not to scale). $N'$ is not well defined for low-$J'$. For clarity, where the spacings are small, arrowheads extending above (below) a doublet pair mark a transition to the upper (lower) doublet.}
    \label{fig:biglevelscheme}
\end{figure}

\section{Experimental setup}
\label{experimental}

The spectrometer used for this study is similar to the one described previously \cite{Doppelbauer2021,Hofsass2021}. MgF molecules are produced in a cryogenic helium buffer gas cell that is cooled to \SI{2.7}{K} using a closed-cycle helium cryocooler. The cell's geometry is based on the design of Truppe \textit{et al.} \cite{Truppe2018}; it has a length of \SI{40}{mm} with a bore diameter of \SI{10}{mm}, and an aperture of \SI{4}{mm}. A Mg rod is ablated by \SI{20}{mJ} of pulsed Nd:YAG laser (Continuum Minilite II, 1064 nm) radiation focused to a waist diameter of 0.4 mm. The hot Mg atoms react with NF$_3$ gas (\SI{0.001}{sccm} flow rate, 100 K) to form MgF molecules. The molecules are cooled by collisions with the cryogenic He buffer gas which is flowing into the cell continuously at a rate of \SI{1}{sccm}. This helium flow also extracts the pulse of molecules from the cell into a molecular beam with a rotational temperature of about 4 K and a mean forward velocity of typically \SI{160}{\meter\per\second}. The forward velocity increases over several thousand ablation shots, but can be restored to its original value by cleaning the cell. The molecules are detected by laser-induced fluorescence (LIF) \SI{44}{cm} downstream from the buffer gas cell aperture \cite{Doppelbauer2021,Hofsass2021}. The transverse velocity spread of the molecular beam is reduced to about \SI{1}{\meter\per\second} by a 2x\SI{2}{\square\milli\meter} square aperture placed at the entrance to the LIF detector. Here, a continuous wave (CW) \SI{359}{nm} laser beam from the second harmonic of a titanium sapphire laser intersects the molecular beam perpendicularly. We ensure that any Doppler shift arising from misalignment of the probe beam is below \SI{10}{MHz} by measuring spectra with and without retroreflecting the laser. The LIF is imaged onto a photomultiplier tube (PMT, Hamamatsu R928) and the resulting photo-current is amplified to give a time-dependent fluorescence signal. We measure and stabilize the fundamental wavelength of the titanium sapphire laser using a wavemeter (HighFinesse WS8-10), which has an absolute accuracy of \SI{20}{\mega\hertz} and a measurement resolution of \SI{0.4}{\mega\hertz}. The wavemeter is calibrated using a temperature-stabilized HeNe laser (SIOS), whose absolute frequency is known to within \SI{5}{\mega\hertz}. Additionally, we determine the $(6s6p)^1P_1\leftarrow(6s^2)^1S_0$ transition frequency and isotope shifts of Yb by applying high-resolution laser spectroscopy to a pulsed buffer gas beam of Yb atoms. These frequencies are known with an absolute and relative uncertainty of better than \SI{1}{\mega\hertz} \cite{Kleinert2016}. Our experimental spectrum is presented in Fig. \ref{fig:Ybspectrum}. The lines under the spectrum show the transition frequencies measured by Kleinert \textit{et al.}\cite{Kleinert2016}. In our spectrum, the line-centers are determined from a fit to multiple Lorentzian lineshapes and listed in Table \ref{tab:Ybtable}. We reproduce the absolute transition frequencies within \SI{10}{MHz} and the relative frequencies (isotope shifts) to within \SI{1}{MHz} over the range of several GHz. Since the deviation between our measured line centers and the published values is within the absolute accuracy of the wavemeter, we use the calibration from the HeNe laser without further correction. We found it was necessary to add a constant flow of dry nitrogen gas through the doubling cavity, in order to avoid absorption of the fundamental light by water vapor. This effect was particularly pronounced near the Q$_1$(0) fundamental
frequency of \SI{834.3255}{\tera\hertz}.

\begin{figure}
    \centering
    \includegraphics[]{FIG2_Yb.eps}
    \caption{Laser induced fluorescence spectrum of the $(6s6p)^1P_1\leftarrow(6s^2)^1S_0$ transition in Yb used to verify the accuracy of our wavemeter. The black dots represent the data recorded in this study and the blue curve shows a fit using Lorentzian lineshapes. The red sticks represent the line positions obtained by Kleinert \textit{et al.}\cite{Kleinert2016} with an absolute accuracy of better than \SI{1}{\mega\hertz}. The inset shows the spectral line from the $^{174}$Yb isotope in more detail.}
    \label{fig:Ybspectrum}
\end{figure}

\begin{table}
\caption{\label{tab:Ybtable} The measured $(6s6p)^1P_1\leftarrow(6s^2)^1S_0$ transition frequencies of Yb relative to the transition frequency of $^{174}$Yb as determined by Kleinert \textit{et al.} \cite{Kleinert2016}. For the isotopes with a nuclear spin of $I\neq0$, $F'$ is given in brackets. The last column gives the absolute frequency differences of the isotope shifts between the two measurements with a standard deviation (SD) of \SI{1.3}{\mega\hertz}.}
\begin{ruledtabular}
\begin{tabular}{lccc}
Isotope&Isotope shift (MHz)\cite{Kleinert2016}&This study (MHz)&$\Delta f$\\
\hline
176 & -508.89$\pm$0.09 & -502.11$\pm$0.13 & 6.78\\
173 ($F'=5/2$) & -250.78$\pm$0.33 & -243.79$\pm$0.43 & 6.99\\
174 & 0 & 4.16$\pm$0.05 & 4.16 \\
172 & 531.11$\pm$0.09 & 536.51$\pm$0.08 & 5.40\\
173 ($F'=7/2$) & 589.75$\pm$0.24 & 595.14$\pm$0.24 & 5.39\\
171 ($F'=3/2$) & 835.19$\pm$0.20 & 839.57$\pm$0.20 & 4.38\\
171 ($F'=1/2$) & 1153.68$\pm$0.25 & 1160.96$\pm$0.54 & 7.28\\
170 & 1190.36$\pm$0.49 & 1196.78$\pm$0.63 & 6.42\\
168 & 1888.80$\pm$0.11 & 1892.44$\pm$0.11 & 3.64\\
\end{tabular}
\end{ruledtabular}
\end{table}

To determine the electric dipole moment of MgF in the A$^2\Pi$  and X$^2\Sigma^+$  states, we install transparent copper mesh electrodes below and above the molecular beam to apply electric fields to the molecules inside the LIF detector. The distance between the electrodes is measured to be \SI{9.0(3)}{\milli\meter}. The voltage on the electrodes is supplied by a high-voltage power supply (Spellman SL1200) and measured with a calibrated high-voltage probe and multimeter with a combined relative accuracy of 10$^{-4}$.

\section{Isotope shifts, spectroscopic constants, hyperfine structure, and $\Lambda$-doubling}
The vibrational, rotational, fine and hyperfine structure of the ground electronic state of MgF is well known. To improve the spectroscopic parameters for the rotational, fine and hyperfine structure of the A$^2\Pi$  state, we record low-$J$ rotational lines of the A$^2\Pi \leftarrow $X$^2\Sigma^+$  transition. The isotope shifts between the $^{26}$MgF and $^{24}$MgF are discussed in section \ref{sec:isotopeShifts}. From section \ref{sec:hyperfineStructure} onwards, we focus on the most abundant $^{24}$MgF, summarizing the results of our measurements, and point out some important differences with the other group II monofluorides.

\subsection{Isotope shifts}
\label{sec:isotopeShifts}

We use a Mg metal ablation target with a natural isotopic abundance of 79\%, 10\% and 11\% for $^{24}$Mg, $^{25}$Mg, and $^{26}$Mg, respectively. $^{24}$Mg and $^{26}$Mg are bosons with a nuclear spin $I(^{24}\textrm{Mg})=0$ and $I(^{26}\textrm{Mg})=0$ whereas $^{25}$Mg is a fermion with a nuclear spin $I(^{25}\textrm{Mg})=5/2$. Fluorine has one stable isotope with a nuclear spin of $I(^{19}\textrm{F})=1/2$.

Fig. \ref{fig:R22isotopes} shows a typical spectrum when exciting the R$_2$/Q$_{21}$(1) line. In this example we scan over the three MgF isotopologues, and observe isotope shifts of \SI{-3.35}{GHz} ($^{25}$MgF) and \SI{-6.48}{GHz} ($^{26}$MgF), relative to $^{24}$MgF. The inset of Fig. \ref{fig:R22isotopes} shows the more complex hyperfine structure of $^{25}$MgF. Each rotational line is split into two groups of hyperfine lines separated by about 1GHz, which are both further split by a few hundred MHz. The larger splitting arises from the Fermi interaction between the electron and $^{25}$Mg nuclear spin in the ground electronic state\cite{Anderson1994}. We did not analyse the excited state in detail in this study.

\begin{figure}
    \centering
    \includegraphics{FIG3_Isotopes_2.eps}
    \caption{Isotope shift of the R$_2$/Q$_{21}$(1) line of MgF. The $^{26}$Mg and $^{24}$Mg isotopes are bosons with $I=0$, while $^{25}$Mg is a fermion with $I=5/2$ and has a more complex hyperfine structure. A zoom-in on the $^{25}$MgF isotope is shown as an inset.}
    \label{fig:R22isotopes}
\end{figure}

The shift in transition frequencies between isotopologues is characteristic to a molecular species. This molecular isotope shift can be used as an additional means of identification, and to reveal small deviations from the Born-Oppenheimer approximation. Since $^{24}$MgF and $^{26}$MgF exhibit the same hyperfine structure, the shift in the gravity center of a rotational line can be found straightforwardly by comparing the positions of equivalent hyperfine peaks. We define the transition isotope shift, $\delta\nu_{i}$, as,

\begin{equation}
\delta\nu_{i} = \nu_i(^{26}\mathrm{MgF}) - \nu_i(^{24}\mathrm{MgF})
\end{equation}

\noindent with $\nu_i(j)$ being the absolute frequency of the optical transition $i$ in isotopologue $j$. Within the Born-Oppenheimer approximation, the isotope shift comes about through changes to the relevant reduced masses in the molecular system. For the rovibrational constants, the relevant reduced mass is that computed from the two atomic masses, $m_{\mathrm{Mg}}$ and $m_{\mathrm{F}}$,
\begin{equation}
m_{mol} = \frac{m_{\mathrm{Mg}}m_{\mathrm{F}}}{m_{\mathrm{Mg}}+m_{\mathrm{F}}}\hspace{0.3cm}.
\end{equation}
For the electronic contribution $T_e$, it is that of the valence electron mass, $m_e$, and the remaining molecular mass,

\begin{equation}
 m_{el} = \frac{m_e(m_{\mathrm{MgF}}-m_e)}{m_{\mathrm{MgF}}} \hspace{0.3cm}.
\end{equation}

\noindent We define $\rho = \sqrt{m_{mol}(^{24}\mathrm{MgF})/m_{mol}(^{26}\mathrm{MgF})}$, which has the value of about $0.983$, and $\rho_{el} = m_{el}(^{24}\mathrm{MgF})/m_{el}(^{26}\mathrm{MgF})$, noting that $1-\rho_{el} = 5.67\times 10^{-7}$. To predict the isotope shifts, we use equations \eqref{eq:XstateEnergies}, \eqref{eq:AstateEnergies} and \eqref{eqn:E0definition} to calculate the energy differences, applying the following relations,
\begin{equation}
\begin{gathered}
B^* = \rho^2 B,  \hspace{0.1cm}\\
\omega_e^* = \rho \omega_e, \hspace{0.1cm}
\omega_ex_e^* = \rho^2 \omega_ex_e, \\
T_e^* = \rho_{el}^{-1}T_e, \hspace{0.2cm}
\label{eqn:isotopeShiftRelations}
\end{gathered}
\end{equation}

\noindent where the asterisks refer to the constants for $^{26}\mathrm{MgF}$. The constants $\omega_e, \omega_ex_e$ are taken from the work of Novikov and Gurvich \cite{Novikov1971}\footnote[2]{While Barber \textit{et al.}\cite{Barber1995} provided more accurate values for the ground state, it is only the difference in ground and excited state constants that matters.}. The remaining values are taken from tables \ref{tab:gdstateparams} and \ref{tab:exstateparams} in the subsequent sections of this article. From the difference in vibrational constants we expect an isotope shift of \SI{-7.33}{GHz}, while the change in $T_e$ contributes $+$\SI{0.47}{GHz}.

\begin{figure}
	\includegraphics[]{FIG5_MgF_isotope_shifts.eps}
	\caption{Isotope shifts of the A$^2\Pi \leftarrow $X$^2\Sigma^+$  transition in MgF. The blue joined points are calculated as described in the text, and the red joined points are from our dataset. (a) The A$^2\Pi_{1/2} \leftarrow $X$^2\Sigma^+$  branches. (b) The A$^2\Pi_{3/2} \leftarrow $X$^2\Sigma^+$  branches.}
	\label{fig:isotopeShifts}
\end{figure}

In Fig. \ref{fig:isotopeShifts} a and b, we plot the calculated $\delta\nu_i$ values for the four branches of the A$^2\Pi_{1/2} \leftarrow $X$^2\Sigma^+$  and A$^2\Pi_{3/2} \leftarrow $X$^2\Sigma^+$  transitions. We plot the predicted values up to $N=7$, and compare with those available from our measurements. We observe a systematic difference between the measured and calculated values of $450-500$~MHz; the mean of the differences is 470~MHz and their standard deviation is 30~MHz. Thus, while the rotational dependence of the isotope shift is well described by relations \eqref{eqn:isotopeShiftRelations}, the combined shift of the electronic and vibrational terms is not. An inaccuracy of 1.8~cm$^{-1}$ in any one of the vibrational constants would be required to account for our observations, which seems unlikely. According to Hougen\cite{Hougen1970}, we should also include a term $B\expect{L_{\perp}^2}$ in the definition of $E_0$, which accounts for the component of the electronic orbital angular momentum perpendicular to the internuclear axis. A non-zero value of $\expect{L_{\perp}^2}$ would increase the discrepancy, by up to 540~MHz. It is therefore likely that the specific mass shift contribution to $T_e$, the field shift of the Mg nucleus, or deviations from the Born-Oppenheimer approximation are responsible for the additional shift. These can only be derived by more sophisticated calculations, and our values provide an important benchmark in this regard.

\subsection{Spectroscopic constants of the A$^2\Pi$  state}
\label{sec:hyperfineStructure}
To determine the spin-orbit, rotational, spin-rotation, $\Lambda$-doubling, and hyperfine constants of the A$^2\Pi$  state, we record 25 hyperfine-resolved rotational lines of the A$^2\Pi_{1/2} \leftarrow $X$^2\Sigma^+$  and A$^2\Pi_{3/2} \leftarrow $X$^2\Sigma^+$  transitions. The lines are slightly broadened by a small, uncompensated ambient magnetic field of 0.8 G in the LIF detector, and the effect of optical pumping between hyperfine and rotational states (discussed in section \ref{sec:lifetime}); they are fitted using a sum of Lorentzian lineshapes. Using a Voigt profile did not change the fit residuals significantly. In this way, we determine the centers of 56 hyperfine lines for $N''\leq 4$, the observed frequencies of which are listed in Table \ref{tab:linelist} of the Appendix. These line centers were used in a least-squares fit to determine spectroscopic parameters for the A$^2\Pi$  state, and we list the best fitted values in Table \ref{tab:exstateparams}. In our analysis, we fixed the ground state parameters to those determined by Anderson \textit{et al.}\cite{Anderson1994}, and reproduce these in Table \ref{tab:gdstateparams} for reference. The A$^2\Pi$  state is well approximated by a Hund's case (a) coupling scheme for the angular momenta. Relevant details regarding the Hamiltonian are provided in Appendix \ref{Appendix:Hamiltonian}. The $\Lambda$-splitting is determined by the linear combination $p+2q$, and the hyperfine splittings are determined by the parameters $a$, $b_F+2c/3$ and $d$ \cite{Frosch1952a}\footnote[2]{\noindent Note that $b_F$ is related to $b$ and $c$ in equation 6.5 of Frosch and Foley \cite{Frosch1952a} by $b_F=b+c/3$.}. To independently measure $p$ and $q$, or $b_F$ and $c$, requires exciting to higher $J$ levels where the Hund's case (a) approximation breaks down. These parameter pairs are otherwise strongly correlated when fitted separately and so we state the linear combinations in Table \ref{tab:exstateparams}. The same reasoning applies for the parameters $A$ and $\gamma$.

\begin{table}
\caption{\label{tab:gdstateparams} Spectroscopic parameters of the X$^2\Sigma^+$  state of MgF, reproduced from the work of Anderson \textit{et al.}\cite{Anderson1994}. In the original article, the hyperfine structure parameter of the fluorine nucleus $b(\textrm{F})=b_F(\textrm{F})-c(\textrm{F})/3$ was used.}
\begin{ruledtabular}
\begin{tabular}{ld}
Parameter &  \multicolumn{1}{c}{Value (MHz)}  \\
\hline
 $E_0(X^2\Sigma^+)$              &          0.0       \\
 $B              $ &      15496.8125    \\
 $     D         $ &          0.03238   \\
 $\gamma         $ &         50.697     \\
 $b_F(\textrm{F})         $&        214.2       \\
 $c(\textrm{F})           $ &        178.5       \\
\end{tabular}
\end{ruledtabular}
\end{table}

\begin{table}
\caption{\label{tab:exstateparams} Experimentally determined spectroscopic constants of the A$^2\Pi$  state of MgF and their standard deviation (SD). The $\Lambda$-doubling parameters $p$ and $q$, and the spin-orbit ($A$) and spin rotation ($\gamma$) constants are strongly correlated. We state their linear combinations $p+2q$ and $A+\gamma$ that are well constrained by the fit.}
\begin{ruledtabular}
\begin{tabular}{lddd}
Parameter & \multicolumn{1}{c}{Value (MHz)}  & \multicolumn{1}{c}{SD} & \multicolumn{1}{c}{SD$\sqrt{Q}$\footnote{SD$\sqrt{Q}$ includes correlations
between parameters\cite{Watson1977}.}}\\
\hline
 $E_0(A^2\Pi)-\frac{1}{2}\gamma    $ &
 834855315.2&         3 &       4\\
 $A+\gamma              $ &    1091346 &         3 &      3 \\
 $B              $ &      15788.2 &      0.3 &         0.4\\
 $\gamma         $ &        -53 &         16 &       136\\
 $p+2q           $ &         15 &         2 &       2\\
 $a($F$)           $ &        109 &          6 &         7\\
 $b_F($F$)+2c/3    $ &        -52 &         14 &        16 \\
 $d($F$)           $ &        135  &          7  &         7\\

\end{tabular}
\end{ruledtabular}
\end{table}

\subsection{Hyperfine structure}

Fig. \ref{fig:levelschemeandspectra} a and b show hyperfine-resolved spectra of the P$_{1}$/Q$_{12}$(1) and Q$_{1}$(0) lines of the A$^2\Pi_{1/2} \leftarrow $X$^2\Sigma^+$  transition respectively. These lines originate from different rotational levels in the X$^2\Sigma^+$  ground state and reach opposite parity levels in the same $J'=1/2$ level of the excited state. The structure of the P$_{1}$/Q$_{12}$(1) line is dominated by the ground state fine and hyperfine interactions, and the excited state hyperfine interaction is not resolved. However, in the negative parity $\Lambda$-doublet, the excited state splitting is \SI{179}{\mega\hertz}, and we resolve both the ground and excited state hyperfine structure in our spectra of the Q$_{1}$(0) line (Fig. \ref{fig:levelschemeandspectra} b). This is caused by a dependence of the magnetic hyperfine interaction on the sign of $\Lambda$, and therefore a difference in the linear combinations of $\Lambda$ states \cite{Frosch1952a}. The magnetic hyperfine constant $d($F$)$ encapsulates this effect, and its influence on the hyperfine splittings is illustrated in the level scheme of Fig. \ref{fig:levelschemeandspectra} c.\\

The hyperfine parameters $a,b_F,c$ and $d$ are related to properties of the electronic wavefunction in the molecule. A first-order approximation was initially described by Frosch and Foley\cite{Frosch1952a} and then subsequently simplified and corrected by Dousmanis\cite{Dousmanis1955}\footnote[3]{In particular, we note the correction to the value of $d$, acknowledged by Frosch and Foley.}. The interaction between the electron and nuclear magnetic moment can be split into two parts: one part is sensitive to the electron density at the nucleus, contained in $b_F$ and the other is sensitive to the orbital and spin wavefunction away from the nucleus, which determines $a, c$ and $d$. The coordinates of the electron relative to the interacting nucleus are expressed in the form $(r_1,\chi)$, where $r_1$ is the distance from the nucleus and $\chi$ is the opening angle subtended with respect to the internuclear axis. According to Dousmanis\footnote[4]{Here we use SI units, and our equations relate to the CGS units of Dousmanis by a factor $\mu_0/(4\pi)$.},

\begin{equation}
\label{eq:Dousmanis}
\begin{aligned}
 a  &= \frac{\mu_0}{2\pi}\mu_B\mu(\mathrm{F}) \expect{1/r_1^{3}},\\
 b_F  &= b + c/3 = \frac{4}{3}\mu_0\mu_B\mu(\mathrm{F})\psi^2(0),\\
 c &= \frac{3}{4\pi}\mu_0\mu_B\mu(\mathrm{F})\expect{\frac{3\cos^2\chi-1}{r_1^3}},\\
 d &= \frac{3}{4\pi}\mu_0\mu_B\mu(\mathrm{F})\expect{\sin^2\chi/r_1^3}
\end{aligned}
\end{equation}

\noindent Here, $\mu(\mathrm{F}) = 5.25\mu_N$ is the magnetic moment of the fluorine nucleus \cite{Stone2005}, $\mu_N$ is the nuclear magneton, and $\mu_B$ is the Bohr magneton. The angled brackets denote expectation values of the A$^2\Pi$  electronic wavefunction, and $\psi^2(0)$ represents its probability density at the fluorine nucleus. From the value of $a$ we find a typical radius $\tilde{r} = \expect{r_1^{-3}}^{-1/3} = \SI{0.88}{\angstrom}$, roughly half the internuclear equilibrium separation of \SI{1.75}{\angstrom}. Approximating $\expect{\sin^2\chi/r_1^3}\approx \expect{\sin^2\chi}/\tilde{r}^3$ leads to a typical value $\tilde{\chi} = \ang{65.3}$. A wavefunction uniformly distributed over $\chi$ has $\tilde{\chi}=\ang{54.5}$. These observations suggest that, much like in the ground electronic state, there is appreciable electron density between the two nuclei.\\
This is in stark contrast with the other group II monofluorides CaF, SrF and BaF, where there is no resolvable interaction with the fluorine spin \cite{Bernath1980,Bernath1981,Kaendler1989,Chen2016a,Xu2017}, and upper limits of a few MHz have been inferred for $d($F$)$. It is only for the fermionic isotopologues of these molecules, where the metal nucleus has non-zero spin, that structure has been observed, and $d(^{87}\textrm{Sr})$, $d(^{135}\textrm{Ba})$, and $d(^{137}\textrm{Ba})$ could be determined \cite{Le2009,Steimle2011}.\\
Finally, we note that the hyperfine interaction causes mixing between the $X^2\Sigma^+(N=3,F=2)$ and $X^2\Sigma^+(N=1,F=2)$ states, and between the $A^2\Pi(J=1/2,F=1)$ and $A^2\Pi(J=3/2,F=1)$ states. This mixing results in losses from the P$_{11}/$Q$_{12}$(1) optical cycling transition. For $^{24}$MgF we calculate a branching ratio of $1.6\times 10^{-6}$ for $^{24}$MgF to the $N=3$ states. For $^{25}$MgF these losses are estimated to be an order of magnitude larger. A more detailed analysis of this effect is given in Appendix A.

\subsection{$\Lambda$-doubling}

In their discussion of the $\Lambda$-doubling in MgF, Walker and Richards derived values of $p$ and $q$ for MgF by extrapolation of the $\Lambda$-splitting observed at $J'+1/2 > 16$ \cite{Walker1970}. Their analysis gives $p+2q= \SI{-50.3}{MHz}$, which disagrees with our result in both the sign of the interaction and the magnitude. This discrepancy is likely due to the omission of the spin-rotation interaction in their analysis of the $^2\Sigma^+ - ^2\Sigma^+$ bands of MgF, and due to insufficient resolution in the A$^2\Pi \leftarrow \mathrm{X}^2\Sigma^+$ absorption spectra. We here provide an improved measurement and update their discussion of the origin of $\Lambda$-doubling in the group II monofluorides. In Table \ref{tab:vanVleckComp}, we compare measured values of $p + 2q$ with the values obtained according to Van Vleck's pure precession approximation \cite{VanVleck1929,Hinkley1972}. The approximation is valid when the $\Lambda$-doubling is dominated by the interaction with a single $^2\Sigma$ state, whose $\sigma$ molecular orbital is mainly derived from an atomic $p$ orbital \cite{Mulliken1931}. Under this assumption,
\begin{equation}
\begin{aligned}
    &p= p_{\textrm{vv}}= \frac{2ABl(l+1)}{E_{\Pi}-E_{\Sigma}} \hspace{0.3cm},\\
    &q= q_{\textrm{vv}} = \frac{2B^2l(l+1)}{E_{\Pi}-E_{\Sigma}} \hspace{0.3cm}.
    \label{eq:vv}
\end{aligned}
\end{equation}

\noindent Here, $l=1$ is the orbital angular momentum of the unpaired electron, and $E_{\Pi}-E_{\Sigma}$ is the energy difference between the interacting states. Equations \eqref{eq:vv} apply for the interaction with a $\Sigma^+$ state, with the sign reversing when the interaction is with a $\Sigma^-$ state. In the table we give $p_{\mathrm{vv}} + 2q_\mathrm{{vv}}$ values for pure precession with the nearby  B$^2\Sigma^+$ states which are higher in energy, and this appears to work well for the heavier monofluorides CaF, SrF and BaF. The trend is primarily due to the decreasing spin-orbit interaction moving up the group. In the case of MgF, the interaction changes sign and is an order of magnitude smaller than expected from the interaction with the B$^2\Sigma^+$ state. Therefore, $\Lambda$-doubling in MgF is more complex than for the heavier group II monofluorides and may comprise interactions with many $\Sigma$-states.
\textit{Ab initio} calculations by Kang \textit{et al.} \cite{Kang2015} imply that the B$^2\Sigma^+$ state of MgF is of mixed character, which may explain the marked difference in $\Lambda$-doubling compared to the heavier group II monofluorides.

\begin{table}
\caption{\label{tab:vanVleckComp} A comparison of the experimentally obtained $\Lambda$-doubling constants of the A$^2\Pi$  state $p+2q$, with the value of $p_{\textrm{vv}}+2q_{\textrm{vv}}$ as discussed in the text, for the group II monofluorides. All values are stated in \SI{}{MHz}, and for MgF we use values obtained in this study. Sources of spectroscopic data are shown in the column headings where needed.}
\begin{ruledtabular}
\begin{tabular}{lcccc}
 & \multicolumn{1}{c}{MgF \cite{Barrow1967}} & \multicolumn{1}{c}{CaF \cite{Kaledin1999}} & \multicolumn{1}{c}{SrF \cite{Steimle1993,Steimle1977}} & \multicolumn{1}{c}{BaF \cite{Barrow1967a}} \\
\hline
   $10^{-3}(p+2q)$ (Exp.) &    \SI{0.015}{}$^a$  & \SI{-1.32}{}  & $\SI{-3.90}{}^b$  & $\SI{-7.67}{}^b$ \\
    $10^{-3}(p_{\textrm{vv}} + 2q_{\textrm{vv}})$ & \SI{-0.25}{}    &   \SI{-1.31}{} &   \SI{-4.17}{} &     \SI{-6.66}{} \\
\end{tabular}
\end{ruledtabular}
\footnotetext{This work.}
\footnotetext{In these studies, only $p$ is reported, and $q$ is assumed to be zero.}
\end{table}

\begin{figure*}
    \centering
    \includegraphics[]{FIG5_Hyperfine_proof.eps}
    \caption{Hyperfine resolved spectra of the P$_1$/Q$_{12}$(1) (a) and Q$_1$(0) (b) lines. Experimental data are shown as black dots, and the line positions simulated from the fitted spectroscopic constants are shown as colored bars. The dark blue lines are Lorentzian fits to the experimental data to determine the line-center. (c) Energy level scheme of the hyperfine states involved in the transitions shown in (a) and (b). Energy differences are given in MHz. The two hyperfine states in the positive parity doublet in A$^2\Pi_{1/2}$ , $J=1/2$ are not resolved; the calculation from spectroscopic constants determines the spacing to \SI{0.4}{\mega\hertz}. (d) Dependence of the linewidth of the three hyperfine lines of the Q$_1(0)$ spectrum as a function of the peak laser intensity. The hyperfine lines broaden differently with increasing intensity because of hyperfine and rotational pumping. Straight lines are linear fits to the data. The relative slopes agree well with simulations using rate equations.}
    \label{fig:levelschemeandspectra}
\end{figure*}

\section{Radiative lifetime of the A$^2\Pi$ ,$v'=0$ level}
\label{sec:lifetime}
So far, the radiative lifetime of the A$^2\Pi$  state is known only theoretically \cite{Pelegrini2005}. We determine the lifetime experimentally by measuring the Q$_1$(0) spectral line shape (Fig. \ref{fig:levelschemeandspectra} b) at low laser intensity. This line is convenient because the hyperfine structure in both ground and excited states is fully resolved, and the ground $N=0$ rotational level contains the largest population of slow molecules for which the residual Doppler broadening is smallest. For these measurements, we reduced the expected Doppler broadening to \SI{1}{MHz} by replacing the 2x\SI{2}{\square\milli\meter} square molecular beam aperture in the detector with a \SI{1}{mm} slit. We measure spectra at several probe laser intensities, and extract the full width at half-maximum (FWHM) for each line with a Lorentzian fitting function. Fig. \ref{fig:levelschemeandspectra} d plots the linewidths for the three hyperfine components of the Q$_1$(0) line against the peak laser intensity together with linear fits (solid lines). Linear fits to the data show that each line broadens differently with increasing laser intensity. This broadening occurs at laser intensities well below the predicted two-level saturation intensity, $I_s=\pi h c\Gamma/3\lambda^3=$ \SI{62}{\milli\watt\per\centi\meter\squared}, and is the result of optical pumping between rotational and hyperfine states of the molecule, which we discuss in the following paragraph.

In general, the spectrum of an open transition will broaden when the number of photon scattering events is sufficient to optically pump the molecule to a state not addressed by the laser. This broadening can occur at an arbitrarily low laser intensity $I$, provided that the interaction time $t_i$ is large enough. In the absence of hyperfine structure, molecules are pumped to $N''=2$ on the Q$_{1}(0)$ line after an average of three scattering events. With the inclusion of hyperfine structure, molecules are pumped to both $N''=2$ and also between hyperfine levels of $N''=0$, further reducing the number of scattering events. For our experiments where the typical laser interaction time is $t_i=\SI{10}{\micro\second}$, this effect becomes significant even when $I\sim 10^{-3}I_{s}$. To verify our understanding, we simulate the interaction with the laser using rate equations, the measured hyperfine splittings, and the branching ratios for each hyperfine decay channel. We assume all three polarization components are excited with equal probability; this is a reasonable approximation given the magnetic field in the detector mixes the ground states by spin precession during the interaction time with the laser. We fit the simulated spectra with the same Lorentzian model, and find that the model predicts the relative broadening rates within the experimental uncertainties. A complete quantitative treatment of this effect requires detailed information about the laser profile and the collection optics, which is beyond the requirements of this paper. We refer the interested reader to the work of Wall \textit{et al.}\cite{Wall2008} as an example.

To estimate an upper bound for the true Lorentzian linewidth $\Gamma/(2\pi)$, where $\Gamma=1/\tau_0$ and $\tau_0$ is the radiative lifetime of the A$^2\Pi$ ,$v'=0$ level, we use the spectrum taken at a peak intensity of \SI{80}{\micro\watt\per\centi\meter\squared}, shown in Fig. \ref{fig:levelschemeandspectra} b. Here the FWHM of the individual hyperfine lines are consistent within the uncertainties, indicating that the effect of optical pumping is small. We fit the data to a sum of three Voigt profiles, finding a Lorentzian FWHM of $\Gamma/(2\pi)=\SI{22.0(5)}{MHz}$ and a Gaussian FWHM of $3.8(1.8)$\SI{}{MHz}. The ground state splitting can be estimated from our experimental spectra, and compared to the precise measurements of Anderson \textit{et al.} to estimate an uncertainty on the linearity of the laser scan. We find agreement within $\pm \SI{1}{MHz}$, consistent with the Yb measurements presented in Section \ref{experimental}. The Gaussian contribution to the lineshape arises from the Doppler effect, residual Zeeman shifts and laser frequency instability. To estimate its systematic uncertainty, we measured the $(3s^27s) ^2S_{1/2} \leftarrow (3s^23p) ^2P_{1/2}$ narrow transition in a buffer gas beam of atomic aluminum. The atoms are produced in the same beam machine, with a forward velocity similar to the MgF molecules. The probe-light at \SI{225.8}{nm} is generated from the fourth harmonic of the same titanium sapphire laser, which increases the laser frequency noise by at least a factor two. The sensitivity to Doppler shifts arising from the shorter probe wavelength is increased by a factor 1.6 relative to the Q$_{1}(0)$ line in MgF. The Al transition has a natural linewidth of \SI{2.5}{\mega\hertz} \cite{Paschen1932}, and using this value we arrive at an upper bound of the Gaussian FWHM contribution of \SI{9.2(2)}{\mega\hertz} at the Al detection wavelength. We therefore vary the Gaussian contribution in the final fit between 1 and \SI{4.6}{\mega\hertz}, and find that this changes the fitted value of $\Gamma/(2\pi)$ by at most \SI{1}{MHz}. From this, we estimate a lower bound for the radiative lifetime of the A$^2\Pi$ ,$v'=0$ state to be $\tau_0 = 1/\Gamma = (7.23\pm0.16_{stat}\pm0.33_{sys})$ ns, in agreement with the theoretical prediction by Pelegrini \textit{et al.} \cite{Pelegrini2005} of \SI{7.16}{\nano\second}. The uncertainty is dominated by the systematic uncertainty in the Gaussian contribution to the spectral lineshape.

\vspace{-2mm}
\section{Electric dipole moment measurements}

The application of an external electric field $\boldsymbol{E}$ in the detector introduces an additional term,

\begin{equation}
    H_{Stark}=-\boldsymbol{\mu}_i\cdot\boldsymbol{E},
    \label{eqn:HStark}
\end{equation}

\noindent to the Hamiltonian given in Equation \eqref{eq:hamil}. Here, $\boldsymbol{\mu}_i$ is the vector dipole moment operator in electronic state $i$. From the Stark splitting and shifting of the LIF spectra we can determine the magnitude of the dipole moments $|\boldsymbol{\mu}_A|$ and $|\boldsymbol{\mu}_X|$.

The Stark Hamiltonian \eqref{eqn:HStark} couples nearby molecular states of opposite parity having the same total angular momentum projection $M_F$ onto an axis parallel to $\boldsymbol{E}$. To calculate the energies under an applied field, we diagonalize the Hamiltonian matrix constructed using the relevant Hund's case basis functions, including all levels with $J\leq 9/2$. In the electronic ground state the Stark Hamiltonian mainly mixes states separated by one unit of the quantum number $N''$. For the electric fields applied in this study of up to \SI{10.6}{kV\per\centi\meter}, the Stark shift is quadratic. For the excited states, the dominant interaction is between the closely spaced $\Lambda$-doublet levels. The interaction is strong compared to this splitting and results in a linear Stark shift even for low electric field strengths. The Stark effect overcomes the hyperfine interaction at modest electric field strengths of about \SI{0.1}{kV\per\centi\meter}, after which the levels separate by their value of the angular momentum projection $M_{J'}$.

Fig. \ref{fig:starkshifts} a and b show spectra of the  R$_{21}$(0) and Q$_{1}$(0) lines, respectively, under field-free conditions and with an electric field of \SI{4.44}{\kilo\volt\per\centi\meter} applied in the detector. The total span of each spectrum is determined by the Stark effect of the excited state, and the small shift in the gravity center is due to the ground state Stark shift. With the laser polarization oriented perpendicular to the electric field, we excite and observe all four $M_{J'}$ components of the R$_{21}$(0) lines in Fig. \ref{fig:starkshifts} a. We measure spectra at various applied fields and fit the modeled spectra to extract best fit values for the dipole moments, obtaining $\mu_A=3.20\pm0.01_{stat}\pm0.22_{sys}$~D and $\mu_X=2.88\pm0.03_{stat}\pm0.20_{sys}$~D. The systematic uncertainty is dominated by the determination of the electric field between the mesh electrodes. We assume a measurement uncertainty of \SI{0.3}{mm} in the mesh separation of \SI{9}{\milli\meter}, and with finite element modelling we estimate a reduction of $2\%$ in the electric field strength relative to infinite plate electrodes of the same separation. The combination of these effects far exceeds the statistical uncertainty from the fitting procedure. Simulated spectra using the best fit values are shown inverted in Fig. \ref{fig:starkshifts} a and b, demonstrating good quantitative agreement with the measurements. The FWHM of the lines measured at low and high field are consistent within \SI{1}{MHz}, from which we deduce that spatial inhomogeneity of the electric field across the probe beam is below $0.1\%$.

To examine the behavior of the energy levels in MgF in electric fields in more detail, we show the simulated Stark shifts for low-$J$ levels at high, moderate and small electric field strengths in Fig. \ref{fig:starkshifts} c. Between 10 and \SI{100}{kV/cm} (left panel), different rotational levels interact significantly; this leads to avoided crossings between the excited states, indicated by circles in the Fig.\ref{fig:starkshifts}. At intermediate fields up to \SI{10}{kV/cm} (center panel), the Stark effect of the ground state becomes comparable to the excited state. At electric fields below \SI{0.25}{kV/cm} (right panel), mixing in the ground state is negligible, whereas in the excited state the nearby $\Lambda$-states mix significantly.

The X$^2\Sigma^+$  ground state dipole moment has been predicted theoretically  \cite{Langhoff1986,Fowler1991,Kobus2000,Pelegrini2005,Wu2015,Hou2017}. The values obtained by different methods range from \SI{2.67}{D}\cite{Pelegrini2005} to \SI{3.126}{D}\cite{Wu2015}. We note that our measured value of $\mu_X$ is in good agreement with the theoretical value obtained by Fowler and Sadlej\cite{Fowler1991}, $\mu_X =\SI{2.8611}{D}$, using the complete active space self-consistent field (CASSCF) method. Pelegrini \textit{et al.} also calculated $\mu_A=4.23$ D, which is in poor agreement with our experimental results. Table \ref{tab:dipolemoments} compares our experimental values to those of other group II monofluorides. Steimle \textit{et al.}\cite{Steimle2011} reported two different values of the electric dipole moment for the A$^2\Pi_{1/2}$  and A$^2\Pi_{3/2}$  states of BaF, which can be explained by interactions with other states. Within the  statistical uncertainty of 0.3\%, our Stark shifted spectra of MgF are well described by a single dipole moment for the A$^2\Pi$  state. This suggests very little influence of nearby perturbing states on the dipole moment, in contrast to BaF.

\begin{figure*}
    \centering
    \includegraphics[]{FIG6_Stark.eps}
   \caption{Stark effect measurement of the R$_{21}$(0) (a) and Q$_{1}$(0) (b) lines of the A$^2\Pi \leftarrow $X$^2\Sigma^+$  transition in MgF. The field-free spectra are shown in blue and when an electric field of 4.44 kV/cm is applied in red. The simulated spectra assuming $\mu_X=2.88(20)$ D and $\mu_A=3.20(22)$ D for the ground and excited state dipole moments, respectively, are shown in black. (c) Calculated Stark shifts of the $X^2\Sigma^+, N''=0, J''=1/2$, $A^2\Pi_{1/2}, J'=1/2$, and $A^2\Pi_{3/2}, J'=3/2$ levels. The left panel shows the shifts of up to \SI{100}{\kilo\volt\per\centi\meter}, the center shows the region of $0-$\SI{10}{\kilo\volt\per\centi\meter} and the right panel shows a zoom-in to $0-$\SI{0.25}{\kilo\volt\per\centi\meter}. Avoided crossings are marked with black circles. The different $|M_F|$ components are color-coded. If the Stark curves of two states with the same $|M_F|$ overlap at the resolution of the figure, they are shown as dashed lines. The electric field used in the spectra (a and b) is indicated with a dashed vertical line.}
    \label{fig:starkshifts}
\end{figure*}

\section{Electric field induced rotational branching}
The laser cooling scheme for MgF relies on the parity ($P$) and angular momentum ($J$) selection rules of electric dipole transitions. In zero electric field, optical cycling is possible on the P$_1/$Q$_{12}(1)$ transition, which excites molecules from the $N''=1$ ground states to the $J'=1/2, P'=+$ excited states.
However, a small electric field $E$ results in an exchange of population, $\epsilon^2$, between excited states of opposite parity. The relative transition strength from $N''=1$ becomes  $1-\epsilon^2$ when exciting the states of mostly positive character, and these states decay to $N''=0,2$ with probability $\epsilon^2$. In addition, excitation to the states with mostly negative parity character becomes weakly allowed, with a relative transition strength $\epsilon^2$. These states decay to $N''=0,2$ with near unit probability. This second loss channel cannot be neglected, because efficient optical cycling requires frequency sidebands to address the ground state hyperfine splitting, causing the weak transitions to be driven near resonance. As a result, we estimate the loss probability from the cooling cycle due to uncontrolled electric fields, to lowest order in $\epsilon$, as $\mathcal{P}_{\mathrm{loss}}\approx 2\epsilon^2$.

In Fig. \ref{fig:lossrate}, we plot $\mathcal{P}_{\mathrm{loss}}$ versus the electric field strength for the two possible values of $F'$, using our measured spectroscopic constants. We fit the data below 5 V/cm to $\mathcal{P}_{\mathrm{loss}}=\alpha_{F'}|E|^2$, finding that $\alpha_0=4\times10^{-5}\SI{}{cm^2/V^2}$ and $\alpha_1=2\times10^{-4}\SI{}{cm^2/V^2}$. For the P$_1$/Q$_{12}$(1) cooling transition, electric fields of \SI{18}{\volt\per\centi\meter} and \SI{9}{\volt\per\centi\meter} lead to a mixing of the parity eigenstates of about \SI{1.4}{\%}, for $F'=0$ and $1$, respectively. At this level, losses from the optical cycle due to parity mixing match losses to the $v''=1$ vibrational manifold in the X$^2\Sigma^+$  state predicted by Ref. \cite{Pelegrini2005}. We note that such electric field induced losses have already been observed experimentally in SrF \cite{Norrgard2016} and AlF \cite{Hofsass2021}.
\begin{table}
\caption{\label{tab:dipolemoments} Reported dipole moments of group-II monofluorides, in debye. Values for MgF are from this study.}
\begin{ruledtabular}
\begin{tabular}{lcccc}
Parameter & MgF & CaF & SrF & BaF  \\
\hline
 $\mu(X^2\Sigma^+)$              &          2.88(20)     & 3.07(7)\cite{Childs1984} &  3.4963(4)\cite{Ernst1985} & 3.170(3)\cite{Ernst1986} \\
  $\mu(A^2\Pi_{1/2})$              &          3.20(22)     &  &  & 1.50(2)\cite{Steimle2011} \\
   $\mu(A^2\Pi_{3/2})$              &          3.20(22)    & 2.45(6)\cite{Ernst1989} & 2.064(50)\cite{Kaendler1989} & 1.31(2)\cite{Steimle2011}  \\

\end{tabular}
\end{ruledtabular}
\end{table}

\begin{figure}
    \centering
    \includegraphics[]{FIG7_Mixing_review_arxiv.eps}
   \caption{Calculated population mixing of opposite parity levels, $\epsilon^2$, for the different hyperfine levels of the A$^2\Pi$ , $J'=1/2$ states in an electric field. This results in unwanted rotational branching to $N''=0$ and 2 in the ground state when driving the P$_1$/Q$_{12}(1)$ laser cooling transition. The $\left|M_{F'}\right|=0,1$ levels in $F'=1$ are shown in blue, the $F'=0, \left|M_{F'}\right|=0$ level is shown in red. The dashed line indicates the level of expected losses due to vibrational branching to $v''=1$. The quadratic loss coefficients are $\alpha_0=4\times10^{-5}\SI{}{cm^2/V^2}$ and $\alpha_1=2\times10^{-4}\SI{}{cm^2/V^2}$.}
    \label{fig:lossrate}
\end{figure}

\section{Conclusion}
We recorded hyperfine-resolved CW LIF spectra of 25 low-$J$ lines of the A$^2\Pi \leftarrow $X$^2\Sigma^+$  transition in MgF, and analyzed the $^{24}$MgF isotopologue in detail. By fitting the eigenvalues of the effective Hamiltonian to the measured line positions, we determined the spectroscopic parameters of the A$^2\Pi$  state that are relevant for laser cooling experiments: rotational, fine- and hyperfine structure constants. We calibrated our wavemeter using the precisely known Yb $(6s6p)^1P_1\leftarrow(6s^2)^1S_0$ transition frequencies \cite{Kleinert2016}, and correct a \SI{-4.1}{GHz} systematic error in the line frequencies presented by Xu \textit{et al.} \cite{Xu2019}. Transition isotope shifts between the $^{24}$MgF and $^{26}$MgF isotopologues were recorded, and we observe an unexplained \SI{470}{MHz} transition frequency shift which may indicate deviations from the Born-Oppenheimer approximation. We studied the broadening of the hyperfine lines due to optical pumping and recorded high-resolution spectra of the Q$_1(0)$ line to determine the radiative lifetime of the A$^2\Pi$ , $v'=0$ level to be $\tau_0=\SI{7.23(36)}{ns}$.

By studying the fluorescence spectra under an applied electric field, we experimentally determined the dipole moments of the X$^2\Sigma^+$  and A$^2\Pi$  states. Our value for the ground state, $\mu_X=2.88\pm0.03_{stat}\pm0.20_{sys}$~D is in good agreement with the value predicted by \textit{ab initio} calculations \cite{Kobus2000,Langhoff1986,Fowler1991,Wu2015,Hou2017}. Using our value of  $\mu_A=3.20\pm0.01_{stat}\pm0.22_{sys}$~D, we predict the electric field strength at which parity-mixing in the excited states limits optical cycling on the P$_1$/Q$_{12}$(1) line. We find that \SI{9}{\volt\per\centi\meter} is sufficient for unwanted rotational branching to match the expected vibrational branching. To scatter more than $10^4$ photons, stray electric fields have to be controlled to below the \SI{1}{\volt\per\centi\meter} level. Coincidentally, the hyperfine structure in the $J'=1/2,P'=+1$ level is less than \SI{1}{\mega\hertz}, which simplifies the laser cooling scheme significantly. On the contrary, the hyperfine splitting in the $J'=1/2,P'=-1$, level is large, which increases the separation between opposite parity hyperfine levels. This reduces the sensitivity of the optical cycling scheme to stray electric fields substantially.

There are a number of notable differences between MgF and the heavier group II monofluorides CaF, SrF and BaF. First, the sign of the $\Lambda$-splitting is inverted in MgF, and its magnitude is about 100 times smaller. Second, MgF has the largest dipole moment in the A$^2\Pi$  state, whereas it has the smallest dipole moment in its ground state. Third, the interaction of the electronic angular momentum with the fluorine nuclear spin leads to a resolvable hyperfine structure in the excited state of MgF, and from our measured hyperfine constants we can infer that the electronic wavefunction has significant probability density between the nuclei. This supports the conclusions of Anderson \textit{et al.} \cite{Anderson1994} regarding the greater covalency of the chemical bonding in MgF.

Finally, we note that our measurements form a stringent set of benchmarks for precise quantum chemical calculations on MgF. The hyperfine constants and dipole moment measurements are strong benchmarks for molecular orbital calculations of the X$^2\Sigma^+$  and A$^2\Pi$  states, while the transition isotope shifts presented constrain vibrational constants and deviations from the Born-Oppenheimer approximation.

\section*{Acknowledgments} \noindent The authors acknowledge the mechanical and electronic workshops of the Fritz Haber Institute for their expert technical assistance. This project received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 Research and Innovation Programme (CoMoFun, Grant Agreement No. 949119).

\section*{Conflict of Interest}
\noindent The authors have no conflicts to disclose.
\section*{Data Availability}
\noindent The data that support the findings of this study are openly available in Zenodo.com at https://doi.org/10.5281/zenodo.6110184, and may be used under the Creative Commons Attribution 4.0 International license.