\section{Introduction}
Alkaline-earth and alkaline-earth-like atoms, such as strontium and cadmium, have attracted considerable recent interest and research in the field of high-precision metrology and inertial sensing due to a number of favorable characteristics. For example, the presence of two valence electrons provides access to narrow intercombination transitions which have been utilized for e.g. high-precision optical atomic clocks and for atom interferometry. Strontium has received considerable attention, especially in the atomic clock community~\cite{Ushijima_2015,Campbell_2017,Oelker_2019}, but also more recently for atomic interferometry, with recent demonstrations on the narrow $^1S_0\to{}^3P_1$ intercombination~\cite{delAguila_2018,Rudolph_2020} and ultra-narrow $^1S_0\to{}^3P_0$ clock transitions~\cite{Hu_2017,Hu_2019}, as well as being suggested as the test species for both ground-based~\cite{MAGIS_2018,AION_2020} and satellite missions~\cite{SAGE_2019,AEDGE_2020} for gravitational wave and dark matter detection. Cadmium is presently less well developed, though it has attracted renewed interest due to its low susceptibility to background blackbody radiation~\cite{Yamaguchi_2019,Porsev_2020}. Strontium and cadmium also possess favorable features which would be beneficial when operating as a dual-species atom interferometer used to search for possible violations of the weak equivalence principle or other extensions to general relativity~\cite{Tinsley_2019}. The frequencies of all of the $^1S_0\to{}^1P_1$, $^1S_0\to{}^3P_0$ and $^1S_0\to{}^3P_1$ transitions of strontium and cadmium are close to being in a 1:2 ratio (Sr:Cd), allowing for the momentum-transfer and Rabi frequencies in a dual-Bragg interferometer to be almost equal, which is important for ensuring common-mode rejection of noise and other systematic effects~\cite{Barrett_2015}. Additionally, the presence of both fermionic and bosonic isotopes of both atoms allows for potential spin-gravity couplings to be investigated~\cite{Tarallo_2014} and using a Sr-Cd dual interferometer would provide good sensitivity to both standard model extensions and dilaton models~\cite{Damour_2012,Hohensee_2013}.

Crucial to enabling all of these experiments is the laser cooling, trapping and manipulation of atomic samples. A key transition in this respect is the $^1S_0\to{}^1P_1$ transition, with respective wavelengths of 460.86~nm and 228.89~nm for strontium and cadmium. For example, the $^1S_0\to{}^1P_1$ transition is key in the preparation of cold strontium samples, being employed in Zeeman slowers and magneto-optical traps~\cite{Courtillot_2003,Xu_2003}, as well as being used to drive Bragg transitions in atomic interferometers~\cite{Mazzoni_2015}. Similarly for cadmium the deep ultraviolet (DUV) transition at 228.89~nm is used for cooling and trapping, with the DUV produced by frequency-doubling blue light~\cite{Yamaguchi_2019,Brickman_2007,Kaneda_2016}. There are, however, a limited number of options for producing laser light in the region 450-470~nm with the fulfilment of all of the required properties of narrow linewidth, high power and single-mode operation with good mode quality. The impact of this lack of coverage is compounded by an unfortunate lack of similar sources in the corresponding near-infrared (NIR) region (900-940~nm), which could otherwise be used as a master source for the production of the necessary blue light via second-harmonic generation (SHG).

Solutions for producing blue light, either directly or via SHG, with all of these properties do exist. For example, there are high-power, tunable sources in the 900-940~nm region based upon e.g. Ti:sapphire lasers~\cite{Fengqin_2015} or neodymium-doped fibre~\cite{Rota-Rodrigo_2017}, the latter of which can produce Watt-level light in the region 915-937~nm. These systems, however, are complex and typically necessitate a large footprint and expense. Therefore, experiments often rely on tapered amplifier systems which, while capable of producing usable power in the region 500-1000~mW following SHG~\cite{Barbiero_2020}, suffer from bad ageing effects, moderate intensity noise and, especially, poor and variable beam quality, limiting their practical use. More recently, efforts at producing blue laser light directly from diodes have advanced and very recently a commercial system for producing 1~W directly at 460.86~nm has been made available. To achieve competitive powers these systems typically require a master-and-slave configuration with an extended-cavity diode laser (ECDL) used to inject a second diode, with long-term operation limited by the maintenance of good injection and mode quality suitable for only 65\% coupling to a single-mode fibre\cite{Schkolnik_2020}. Due to the relatively new nature of this technology there remain some open questions, however, regarding for example the typical lifetime of the diodes.

In this paper, we instead develop a vertical-external-cavity surface-emitting laser (VECSEL)~\cite{Kuznetsov_1997}, which is tunable in the region between around 910-928~nm, with Watt-level output. VECSELs have the advantage of producing high-power, high-brightness, narrow-linewidth and low-noise tunable single-mode light in a relatively simple configuration which minimizes the mechanical footprint of the laser (19 $\times$ 32~cm$^2$). With respect to its low-noise properties, VECSELs are so-called class-A lasers, with dynamics governed by the photon lifetime and therefore don't exhibit relaxation oscillations. Moreover, they also don't suffer from the amplified spontaneous emission present in amplified semiconductor laser diode systems. Due to this combination of desirable features, VECSELs have been employed recently in quantum technology applications, for example in the demonstration of full-scale functionality for the manipulation of Mg~\cite{Burd_2016} and Be~\cite{Burd_2020} ions, and have been recognised as an attractive solution for other ion and neutral atom systems, for example pertaining to the intercombination transition of strontium~\cite{Moriya_2020} and for the cooling and trapping of cadmium~\cite{Kaneda_2016}. Here we use a common design for an external resonant frequency-doubling, bow-tie cavity capable of producing up to 1~W of stable SHG, which when combined with a suitable doubling crystal effectively extends the tunability and performance of the VECSEL itself into the blue region of the visible spectrum. Specifically, we have developed a system which is suitable to act as a master source for the laser cooling and trapping on the $^1S_0\to{}^1P_1$ transition of both strontium and cadmium. We anticipate good long-term operation for this system, as the VECSEL in the 900-nm range is made with a similar reliable material system to that used in 980-nm diode pump lasers for telecommunication fibre amplifiers, and the optical pumping of the VECSEL is performed with a reliable 808-nm diode, commonly employed to pump neodymium-doped yttrium aluminium garnet (Nd:YAG) lasers.

In addition, we use the VECSEL as a master laser to perform spectroscopy upon a beam of atomic cadmium, identifying the $^1S_0\to{}^1P_1$ transitions from all eight naturally occurring isotopes and measuring their relative and absolute transition frequencies. Whilst absolute measurements of the $^1S_0\to{}^1P_1$ transition of other atoms has received attention~\cite{Laupretre_2020}, measurements on cadmium are more limited and often decades old~\cite{Burns_1956,Moore_1958,NIST_Spectra}. Moreover, some of these data were determined from interferometers operating in air and therefore requiring the refractive index at the transition wavelength of 228.89~nm to be well known. Whilst empirical formulae exist for calculating the refractive index of air as a function of wavelength~\cite{Peck_1972}, the refractive index also depends upon other environmental factors, for example the ambient temperature, pressure and humidity.

Complete measurement of the isotope and hyperfine shifts of the $^1S_0\to{}^3P_1$ intercombination transition of cadmium at 326~nm have been performed~\cite{Maslowski_2009}, but comparable data for the $^1S_0\to{}^1P_1$ transition is lacking, being only available in fragmentary form~\cite{Kelly_1961} or inferable from the published data of the frequency dependence of trapping in magneto-optical traps~\cite{Brickman_2007,Kaneda_2016}. Furthermore, precision measurement of isotope shifts provides important data for theoretical calculations of atomic structure~\cite{Miyake_2019}, as well as presenting the opportunity to search for and constrain potential extensions to the standard model of elementary particles and interactions, such as new long-range or Higgs-like forces~\cite{Delaunay_2017,Berengut_2018}.

\section{Laser Source Setup}
\subsection{Master Laser: NIR VECSEL}
\begin{figure}[t]
\centering\includegraphics[width=0.91\textwidth]{figureVecselGain.pdf}
\caption{a) VECSEL gain mirror structure as depicted with refractive index (black line) on left-hand axis. The simulated standing-wave electric field modulus (red line) is a critical design parameter for the thickness of the layers and is plotted on the right-hand axis for reference. b) Photograph of the processed flip-chip gain mirror on a $3 \times 3 \times 0.3$~mm$^3$ diamond heatspreader. c) Schematic of the gain mirror cooling assembly incorporating the gain mirror on the heatspeader, a cooling mount, a thermoelectric element (TEC), and a water-cooled heatsink.\label{fig:figureVECSELChip}}
\end{figure}
The fundamental laser source is a VECSEL based on a novel gain mirror designed for high-power emission in the 900-930-nm emission band, operating with an 808-nm pump diode source. The monolithic semiconductor multi-layer structure is depicted in Fig.~\ref{fig:figureVECSELChip}a (left axis) and was grown on a GaAs substrate using solid-source molecular beam epitaxy (MBE). The structure incorporates an active region with GaInP window (first-grown layer) and multiple GaInAs quantum wells (QWs) embedded in GaAsP barriers, and a distributed Bragg reflector (DBR) formed by alternating $\lambda$/4-thick layers of AlAs and GaAs. As the carrier confinement of a single low-indium GaInAs QW near the GaAs band-gap $\sim$870~nm is low, a large number of QWs are used to mitigate possible carrier leakage. A total of 24 (12~$\times$~2) QWs are placed in the antinodes of the simulated cavity standing wave; this is a relatively high number compared to typical gain structure operating at the 1-$\mu$m wavelength range (approx. double). The choice is governed by the need to ensure a good carrier collection in the QWs and hence sufficient single-pass gain. For the same reason GaAsP barriers are used~\cite{Zhang_2009}, instead of the more commonly used GaAs barriers, to increase the QW confinement while maintaining a good level of absorption for the 808-nm pump photons and ensuring good optical material quality.

The grown and diced gain mirror is flip-chip bonded~\cite{Kuznetsov_1997} from the mirror back surface onto a high thermal conductance $3 \times 3 \times 0.3$~mm$^3$ diamond heatspreader, after which the GaAs substrate is removed with selective wet-etching, thus revealing the GaInP window/etchstop layer (see photograph in Fig.~\ref{fig:figureVECSELChip}b). The new air/semiconductor surface is antireflection coated using ion-beam sputtering (IBS) to minimize the reflection of the focused pump light (AOI $\sim$30$^\circ$), to optimize the reflection of the intracavity light (AOI 0$^\circ$), and to protect the semiconductor surface. The remaining semiconductor structure on diamond is only $\sim$5~$\mu$m thick and serves as a disk-shaped gain element with an $\sim$80 aspect ratio compared with the used $\sim$400~$\mu$m pump diameter. The disk-shaped geometry enables efficient heat extraction of the pump laser load via the heatspeader diamond, to an actively TEC-stabilized cooling mount, and further to a heatsink cooled with a low-vibration chiller, which helps maintain the whole laser housing at a near constant temperature.

A $\sim$12-cm long air-cavity, similar to~\cite{Burd_2016}, is formed between the gain mirror and an external curved output coupling mirror, ensuring that the pump spot diameter and the simulated fundamental mode diameter are matched. A temperature-stabilized etalon and birefringent filter are inserted inside the cavity to enable single-frequency operation and tuning (see Fig.~\ref{fig:figureSetup}a). The cavity is placed inside a sealed and temperature-stabilized casing for stable operation in changing atmospheric and temperature conditions. The casing is further purged with dry nitrogen to prevent potential deleterious effects from absorption by remaining atmospheric water vapor. Selection of the optical frequency is coarsely performed via the temperature control of the etalon and the temperature and angular control of the birefringent filter. Continuous frequency control across the full $\sim$1.2 GHz free-spectral range and locking of the laser frequency is enabled by a piezoelectric transducer (PZT) on which the output coupling mirror is mounted, providing a measured bandwidth of $\sim$50~kHz.

\begin{figure}[t]
\centering\includegraphics[width=\textwidth]{figureSetup.pdf}
\caption{a) A tunable VECSEL is employed as a master laser source, with the wavelength selectable via the combination of the angle of a birefrigent filter (BRF), the temperature of an etalon, and the position of the outcoupling mirror (OC). After passing through an optical isolator, the output of a VECSEL is mode-matched to a resonant bow-tie cavity with an LBO crystal to produce blue light for experiments with either strontium or cadmium. When performing experiments with Cd, a second Brewster-cut BBO frequency-doubling cavity is employed. The DUV output of this second doubling stage is used to perform spectroscopy on a beam of atoms whilst a reference cavity monitors the change in frequency. For more details, see the text. b) A photograph of the SHG cavity without the lid, showing the cavity dimensions in mm.\label{fig:figureSetup}}
\end{figure}

\subsection{SHG Production}
The output of the VECSEL is directed first through a 60~dB optical isolator before being mode-matched with a single biconvex, spherical lens to a frequency-doubling, bow-tie cavity which uses a lithium triborate (LiB$_3$O$_5$) crystal as the doubling medium (Fig.~\ref{fig:figureSetup}), the main design for which is inspired by a system at 399~nm~\cite{Pizzocaro_2014}. Two separate cavities, which share a common geometrical and mechanical and optical design, are installed, allowing for the fast switching between producing light for either strontium-based or cadmium-based experiments, with the power splitting controlled by a polarising beamsplitter and half-wave plate combination. The cavity geometry is designed to work across the range of the output of the VECSEL~\cite{Freegarde_2001}, with the only substantial difference for switching between the two blue wavelengths being the necessary phase matching angle of the crystal. The reflectivities of the mirrors are suitably high for the wavelengths required for both strontium and cadmium and are thus interchangeable.

An outline of the SHG cavity parameters and a photograph of the opened cavity are shown in Figs.~\ref{fig:figureSetup}a and~b, with distances kept small so as to minimize the size and footprint. Considering the case at 460.86~nm for strontium, two mirrors of radius of curvature 75~mm and with an angle of incidence of 9$^{\circ}$, provide a focus of 40~$\mu$m and 42~$\mu$m, in the tangential and sagittal dimensions respectively, at the centre of the crystal, with the calculated Boyd-Kleinman factor of 0.28 near to the maximum value for this wavelength and crystal length~\cite{Boyd_1968,Freegarde_1997}. Two plane mirrors fold the beam and support a secondary waist of 220~$\mu$m and 244~$\mu$m at their mid-point, again in the tangential and sagittal dimensions respectively. Light is coupled into this cavity via one of these mirrors, which has a transmission of 1\%, selected heuristically so as to optimize the power of the SHG~\cite{Pizzocaro_2014}. The total optical path length is 362~mm and the free-spectral range 829~MHz. The optical finesse of this cavity is around 400, as measured from the linewidth and by inducing sidebands at 10~MHz with an electro-optical modulator to act as frequency reference. The difference to these values when considering the case of the cadmium-specific 457.49~nm are negligible.

The cavity is monolithic in design for enhanced thermal and mechanical stability, being produced from a single block of aluminium and with the mirror mounts being screwed directly into this bulk material. Good phase matching of the crystal is ensured by mounting it on a five-axis aligner (Thorlabs) and providing active temperature control via a thermoelectric cooler. One of the plane mirrors is mounted onto a PZT which controls the optical path length of the cavity and allows it to be stabilized via the H\"ansch-Couillaud locking method~\cite{Hansch_1980}. This mirror is only 6.3~mm in diameter, compared to a diameter of 12.7~mm for the other mirrors, so as to minimize the mechanical load and thereby extend the PZT control bandwidth, with a calculated resonance frequency of order 100~kHz.

Finally, the whole cavity can be placed under vacuum with the lid being sealed by a rubber O-ring and a valve. Without this vacuum element, it is significantly more difficult to achieve a strong and stable lock due to the high circulating powers leading to localized heating within the cavity. The associated change in optical path length of the cavity, arising from the negative thermo-optic coefficient, means that the symmetry with respect to a change of the fundamental laser frequency or optical path length of the cavity is broken. For example, when changing the frequency with the PZT of the cavity, the sweep is effectively slowed when the PZT is pulling and vice versa when the PZT is pushing~\cite{Kaneda_2016}. Operating the cavity under vacuum removes this problem in addition to providing some protection from acoustic noise.

\section{Optical Characterization}
\begin{figure}[t]
\centering\includegraphics[width=\textwidth]{figureVECSEL.pdf}
\caption{a) Power as a function of wavelength with 24~A of pump laser current, close to saturation. b) The power output of the VECSEL at 915~nm as a function of the pump laser current. c) Spectral density of the VECSEL frequency noise when free-running (solid red line) as determined from the error signal of H\"ansch-Couillaud locking system. The free-running linewidth as determined from these measurements is shown in the legend -- see text for details. The background (dashed black line) is measured with the laser running but out of resonance with the cavity. \label{fig:figureVECSEL}}
\end{figure}

In the NIR, the VECSEL is capable of producing an output of greater than 1~W in the region of 908-915~nm and greater than 2~W between 915 and 928~nm, with a peak output of around 3~W at 925~nm (Fig.~\ref{fig:figureVECSEL}a). Moreover, the spatial mode is of a high quality, being nearly circular and having a measured M$^2$~<~1.1. This represents a clear improvement over tapered amplifier systems which, as mentioned above, suffer from poor mode quality. The practical benefits of this would include e.g. highly efficient coupling into a single-mode fibre, but here it is manifest as excellent coupling into the SHG cavity of 90\% or greater, even when using only a single spherical lens to perform the mode matching. This is advantageous as it allows nearly all of the power in the NIR to be available for frequency doubling. In comparison, although it is possible to generate a similar amount of power using an ECDL and tapered amplifier system in this regime, much lower coupling into the SHG cavity is typically achievable ($\sim$50\%), despite usually utilising more complicated optical coupling systems. Furthermore, this efficiency can vary greatly depending upon the exact tapered amplifier chip employed due to variable beam quality (M$^2$~<~1.2-1.7) and also on the operating current of the amplifier which affects the beam shape.

Fig.~\ref{fig:figureSHG}a shows the results of the doubling efficiency from the two resonant bow-tie cavities, which are capable of producing greater than 1~W at the wavelength required for the $^1S_0 \to {}^1P_1$ of strontium at 460.86~nm. Slightly less power is produced at the cadmium-specific wavelength of 457.49~nm, despite a slightly higher observed conversion efficiency, due to the reduced power of the VECSEL at this wavelength (Fig.~\ref{fig:figureVECSEL}a). Additionally, a lower coupling efficiency closer to 70\% was operational during these measurements, limited only by the use of a different set of optics. Fundamentally, there are no limitations which would prevent the possibility of also reaching 1~W at this wavelength. The data is fit using the least-squares method, with the linear losses of the cavity as the only parameter, to a theoretical prediction of the SHG based upon the cavity parameters and the measured powers~\cite{Pizzocaro_2014}. The extracted finesse values of 387$\pm$6 and 367$\pm$7, for the cadmium and strontium cavities respectively, are slightly lower than the values determined by measuring the linewidth, suggesting that even greater SHG powers may be possible with further optimization. The cavity and the laser can remain locked for long periods and the power of the SHG over one hour, as measured with an amplified, biased photodiode, shows variation contained within $\sim\pm$2\%, whilst the fluctuations of the output from the VECSEL over the same period is $\pm$0.5\% (Fig.~\ref{fig:figureSHG}b).

The relative intensity noise (RIN) of the visible output and the NIR input are measured with a small-area, amplified photodiode with a bandwidth of $\sim$100~MHz. Fig.~\ref{fig:figureSHG}c shows that at low frequencies, the RIN of the SHG approaches the limit of the noise of the fundamental laser at around -100~dBc in the region 10-1000~Hz. At higher frequencies, the RIN of the SHG clearly exceeds that of the NIR input to the cavity, due to the limited bandwidth of the locking electronics currently employed. The performance of the VECSEL itself, however, approaches the shot-noise limit at high frequencies for the measurement power sent to the photodiode (8.5~mW).

In addition to providing a long-term, stable lock, the output of the H\"ansch-Couillaud system can be used to estimate the linewidth of the VECSEL output, when free-running and without feedback. By calibrating the slope of error signal with the measured linewidth of the cavity, it is possible to extract the power spectral density of the laser frequency fluctuations ($S_\nu\left(f\right)$), which can in turn be used to provide an estimate of the laser linewidth. There are multiple definitions of the laser linewidth and here we evaluate the linewidth based upon the phase noise spectral density ($\int_{\Delta\nu/2}^{f_{max}} S_\phi\left(f\right)df=1$, $S_\phi\left(f\right)=S_\nu\left(f\right)/f^2$), a so-called fast linewidth, and the root-mean-squared value evaluated at various timescales ($\Delta\nu_{rms}=\int_{1/t}^{f_{max}} S_\nu\left(f\right)df$). The upper limit of the integrals ($f_{max}$) is set by the maximum Fourier frequency of the collected data ($\sim$5~MHz). For the free-running laser, the measured fast linewidth is 8~kHz and $\Delta\nu_{rms}$ is 4~kHz at 100~$\mu$s and 104~kHz at 100~ms timescales. These linewidths are considerably below the natural linewidths of the transitions in question, which are $2\pi\times$30.5~MHz and $2\pi\times$90.1~MHz for the strontium and cadmium $^1S_0 \to {}^1P_1$ transitions, respectively~\cite{Nagel_2005,Xu_2004}. Moreover, this linewidth allows for easy and efficient coupling to the frequency-doubling cavity, which as well as being practically preferable, also leads to a reduction in the intensity noise of the SHG, which is a function not just of the intensity noise of the fundamental source, but also its frequency noise via the lock to the cavity. It may additionally facilitate the direct locking to higher finesse cavities for applications demanding narrower linewidths.

\begin{figure}[t]
\centering\includegraphics[width=\textwidth]{figureSHG.pdf}
\caption{a) SHG power as a function of the power of the fundamental light coupled into the two cavities (blue square -- Cd, red diamonds -- Sr). The data is fit by varying the linear losses of the cavity, with all other parameters fixed by theory. Errors bars represent the systematic error from the uncertainty of the calibration of the power meters ($\pm$3\%), whilst the shaded region represents the uncertainty of the fit. b) Power measurement of the SHG over the period of an hour. c) Relative intensity noise of the fundamental light (red line) incident upon the frequency-doubling cavity and of the SHG (blue line). The dashed black line is the shot-noise limit for the measurement of the fundamental (8.5~mW). The background trace (grey line) is collected with the incident light on the measurement photodiode blocked.\label{fig:figureSHG}}
\end{figure}

\section{Spectroscopy of the $^1S_0 \to {}^1P_1$ Cd Transition}
\begin{figure}[t]
\centering\includegraphics[width=\textwidth]{figureCd.pdf}
\caption{a) The fluorescence from the cadmium beam is imaged onto a photomultiper tube. See text for more details. b) The top panel shows a sample fluorescence signal (solid black line) fit to Voigt profiles of equal width and fixed relative amplitudes given by the known natural abundances (solid red line). The dashed lines show the contributions of the individual transitions. The dotted line shows the expected Doppler-broadened profile for an ideal cadmium gas at the same temperature as the beam. The lower panel shows the mean result of multiple fits to different fluorescence data sets, with the error bars smaller than the shown data points. The separation between the points for $^{112}$Cd and $^{113}$Cd-B are also indistinguishable at this scale. For the fermionic isotopes, the labels A and B respectively refer to the $F=1/2\to F'=1/2$ and $F=1/2\to F'=3/2$ hyperfine transitions. See Table~\ref{tab:tabRef} for precise values.\label{fig:figureCd}}
\end{figure}
The suitability of this laser for atomic physics applications is demonstrated by spectroscopy performed on the 228.89~nm cadmium $^1S_0 \to {}^1P_1$ transition, using the setup shown in Fig.~\ref{fig:figureCd}a. This requires an additional stage of frequency doubling with the visible output of the first bow-tie cavity being sent to a second bow-tie cavity, this time using a beta-barium borate crystal (BBO, $\beta$-BaB$_2$O$_4$) cut at Brewster's angle as the doubling medium~\cite{Kaneda_2016,Wilson_2011,Hannig_2018}. Mode matching to this cavity is performed by a cylindrical telescope, to correct for the walk-off from the LBO cavity, and a spherical telescope and the observed coupling is approximately 70\%. Accounting for this mode-matching efficiency, we estimate that this second cavity (Agile Optic) should allow for the production of 250~mW of DUV assuming a finesse of 350 and an input-coupling mirror of 1\% transmission~\cite{Hannig_2018}.

For these spectroscopic experiments, however, the power is kept low, ensuring the delivered optical intensity is below the saturation intensity ($I/I_{sat}\sim$0.2, $I_{sat}=$0.99~W/cm$^2$), reducing the impact of power broadening and also minimising power-declining ageing effects in the crystal~\cite{Kaneda_2016}. Therefore only approximately 1~mW of collimated DUV output is sent transverse to a beam of natural cadmium in a retro-reflected configuration. By sweeping the PZT of the VECSEL with a triangular wave and keeping the two doubling cavities on resonance, it is possible to continuously tune over 4~GHz in the DUV. The fluorescence signal is measured with a photomultiplier tube and the input polarization is set to be horizontal to maximize the amount of the dipole emission collected. A lens (f=40~mm) and a curved mirror (R=30~mm), placed inside the vacuum chamber vertically offset from the atomic beam, collect the atomic fluorescence, whilst a second lens and narrow optical filter, outside the chamber, focus the light onto the PMT and remove background light (Fig.~\ref{fig:figureCd}a).

The cadmium atomic beam is based upon an oven, with a similar design to one previously used for strontium~\cite{Schioppo_2012}, with capillaries of length 10~mm and internal diameter 0.23~mm used to ensure a collimated beam. The ambient pressure when the oven is at room temperature is $\sim$10$^{-10}$~mbar. Under typical operating conditions of a temperature of around 100$^\circ$C and with a measured flow rate of $\sim$10$^{10}$~atoms/s, the divergence of the atomic beam is suitably low that it allows for different isotopes to be easily identified, as shown in Fig.~\ref{fig:figureCd}b. The observed width of the emission peaks is close to that of the natural linewidth, meaning that saturated absorption spectroscopy need not be used, with the advantage that much lower intensities may be employed. For comparison, when using an atomic vapor of cadmium at the same temperature, a single Doppler profile would cover the transitions of all the isotopes.

In addition to being coupled into the second-stage SHG cavity, a small fraction of the blue light is matched to a reference bow-tie cavity with a free spectral range of 204.2$\pm$0.8~MHz. By fitting the resonance peaks of this cavity to a parabola it provides a calibration that allows the relative frequency of the blue light with time to be determined. Furthermore, the NIR light is simultaneously coupled into a wavemeter (Bristol Instruments 621) which provides an absolute reference at an accuracy of 0.2~ppm, which corresponds to 262~MHz at the transition frequency. Finally, some of the DUV is sent to a reference photodiode, which allows for any power variations across the spectrum to be corrected for.

\begin{table}[t]
  \begin{center}
    \caption{Summary of measured isotope shift relative to the $^{114}$Cd isotope. Error values for this work sequentially represent the error from the fitting procedure and the systematic error arising from the frequency conversion -- see text for details. The two different methods from reference~\cite{Kelly_1961} are presented along with the values determined in this work. The first of these methods distinguishes neither $^{110}$Cd nor $^{112}$Cd from their respective nearby fermionic transitions, $^{111}$Cd and $^{113}$Cd, respectively, which are therefore quoted here with the same value. The second method attempts to extract the contributions of just the bosonic isotopes.}
    \label{tab:tabRef}
    \begin{tabular}{l c c r}
      \hline
      \textbf{Transition} & \textbf{Determination 1~\cite{Kelly_1961}} & \textbf{Determination 2~\cite{Kelly_1961}}  & \textbf{This work / MHz}\\
      \hline

      {$^{106}$Cd} & - & - & 1747.5~$\pm$~4.3~$\pm$~9.7 \\
      $^{108}$Cd & - & - & 1256.0~$\pm$~4.3~$\pm$~7.0 \\
      $^{110}$Cd & 878~$\pm$~17 & 905~$\pm$~35 & 833.2~$\pm$~2.8~$\pm~$4.6  \\
      $^{111}$Cd - F$'$=1/2 & - & - & 587.0~$\pm$~3.1~$\pm$~3.3 \\
      $^{111}$Cd - F$'$=3/2  & 878~$\pm$~17 & - &  871.5~$\pm$~2.9~$\pm$~4.8  \\
      $^{112}$Cd  & 375~$\pm$~15  & 395~$\pm$~30 & 399.2~$\pm$~2.8~$\pm$~2.2  \\
      $^{113}$Cd - F$'$=1/2  & - & - & 150.0~$\pm$~2.9~$\pm$~0.8   \\
      $^{113}$Cd - F$'$=3/2  & 375~$\pm$~15 & - & 401.1~$\pm$~4.0~$\pm$~2.2   \\

      $^{116}$Cd  & - & - & -299.0~$\pm$~2.6~$\pm$~1.7  \\

      \hline
    \end{tabular}
  \end{center}
\end{table}

In order to extract the relative frequency spacing of the different transitions, we fit the collected fluorescence signal to a sum of Voigt profiles, whose relative amplitudes are fixed based upon the known natural isotope abundances and, in the case of the fermionic isotopes $^{111}$Cd and $^{113}$Cd, the relative strength of the different hyperfine transitions. The width of the Voigt profile for the individual transitions is variable but equal for all transitions and the Lorentzian contribution set to the known natural linewidth, neglecting decoherence due to e.g. collisions. Neglecting the contributions of power broadening (6.3~MHz) and time-of-flight broadening (170~kHz), this also allows for the extraction of the atomic beam divergence from the Gaussian contribution to the Voigt profile~\cite{Schioppo_2012}, resulting in an estimated atomic beam divergence (half angle) of approximately 40~mrad. The relative frequency of each transition can be extracted from the fit results and a total of ten readings each whilst the frequency of the laser is being increased and decreased. The initial values provided to the fitting procedure of each isotope shift are selected from a uniform pseudorandom distribution with a width of 20~MHz, centred on values chosen using previous measurements~\cite{Kaneda_2016,Brickman_2007,Kelly_1961}, to reduce possible initial-value biasing effects. Using both positive and negative frequency gradients helps to average out any remaining systematic error from imperfect frequency calibration with the reference cavity.
\begin{table}[!t]
  \begin{center}
    \caption{Absolute frequency measurement for $^1S_0\to{}^1P_1$ transition of $^{114}$Cd. Error bars for this work are sequentially the statistical and systematic errors. The two values from~\cite{Burns_1956} represent values taken with two different types of lamp: a Michelson lamp; and an electrodeless lamp. Also shown are quoted values from natural cadmium.}
    \label{tab:tabAbs}
    \begin{threeparttable}
    \begin{tabular}{l l}
      \hline
      \textbf{Reference} & \textbf{Transition Frequency / MHz}\\
      \hline

      $^{114}$Cd, This Work\tnote{a} & 1309864506 $\pm$ 4 $\pm$ 262  \\
      $^{114}$Cd, Michelson~\cite{Burns_1956}\tnote{b}  & 1309864420 $\pm$ 86 \\
      $^{114}$Cd, Electrodeless~\cite{Burns_1956}\tnote{b}  & 1309864600 $\pm$ 86 \\
      Natural Cd~\cite{Burns_1956}\tnote{b}  & 1309864720 $\pm$ 86 \\
      Natural Cd~\cite{NIST_Spectra}\tnote{c} & 1309864720 $\pm$ 218 \\
      Natural Cd~\cite{Moore_1958}\tnote{c}  & 1309867418 \\

      \hline
    \end{tabular}
    \begin{tablenotes}\footnotesize
    \item[a]Derived from NIR wavemeter.
    \item[b]Derived from wavenumber. Errors are calculated from the quoted wavelengths and their uncertainty of $\pm$0.015~pm, by calculating the refractive index at the transition wavelength ($n-1$=3.084$\times10^{-4}$)~\cite{Peck_1972}.
    \item[c]Derived from wavenumber.
    \end{tablenotes}
    \end{threeparttable}
  \end{center}
\end{table}

The determined frequency separations, relative to the $^{114}$Cd isotope transition, are presented in Fig.~\ref{fig:figureCd}b and summarized in Table~\ref{tab:tabRef}. These values are the mean, weighted by the error given by the fitting procedure, for the total of 20 data sets. The two error values respectively represent the error from fitting procedure and the systematic error which comes from the uncertainty in the frequency conversion of the temporal data. The error arising from the fitting procedure is the quadrature summation of two factors: the standard error of the mean from the fits to multiple data sets; and the mean error from all the fitting procedure of each data set, as parametrized by the covariance matrix, which is not assumed to scale inversely with repeated measurements. The error in the frequency conversion is dominated by the uncertainty in the free spectral range of the reference cavity. Errors such as the Doppler shift, detailed below, are here assumed to be common to all isotopes and therefore ignored for these relative frequency determinations. All the extracted values from the fit are consistent with the previously reported values within 3$\sigma$~\cite{Kelly_1961}.

To measure the absolute frequency of the transition of the $^{114}$Cd isotope, the wavemeter rather than the reference cavity is employed to determine the frequency axis. In this case, multiple data sets are recorded with the laser sweeping across the $^{114}$Cd transition and the traces, which are noisy due to the inherent fluctuations of the wavemeter, are fit using the same procedure as above. The mean and standard deviation of these fits give a measured value of 1309864506$\pm$4$\pm$262~MHz, where the first error is the statistical standard error of the mean and the second the systematic error. In addition to including the known accuracy of the wavemeter, the final systematic error is calculated from adding this error in quadrature to the errors arising from the first- and second-order Doppler shifts and the recoil energy. The largest of these additional systematic errors is the first-order Doppler shift, which arises from imperfect orthogonality between the atomic beam and the probe laser beam. We estimate this angle to deviate from perpendicularity by $<4$~mrad, giving an additional error of 5~MHz, meaning that the accuracy of the wavemeter remains the dominant contribution. This value of the transition frequency is consistent with references~\cite{Burns_1956,NIST_Spectra}, but not with reference~\cite{Moore_1958} (Table~\ref{tab:tabAbs}).

\section{Summary \& Outlook}
A tunable, narrow-linewidth VECSEL system has been used to generate Watt-level light in the blue region of the visible spectrum via SHG in a resonant bow-tie cavity, which is stable to the $\sim$2\% level on the timescales of one hour. The relative simplicity, narrow linewidth and good spatial mode make this system an attractive proposition for applications in both the NIR and SHG regions. The parameters obtained address the need for laser cooling and trapping of strontium and cadmium atoms, where higher powers allow for faster sample preparation times. This is an important requirement for increasing the duty cycle of experiments, which is a key target for increasing the sensitivity of many experiments.

The laser has been used as a master source in conjunction with a beam of cadmium atoms to measure the frequency of the $^{114}$Cd $^1S_0 \to {}^1P_1$ transition. New values for the isotope shifts of this transition have also been provided for all naturally occurring isotopes. To achieve a better determination of some of the individual isotopes, such as distinguishing the $^{112}$Cd transition from the $F=1/2\to F'=3/2$ $^{113}$Cd transition, it may be required to use enriched samples. When used in conjunction with the isotope shifts of other transitions of cadmium, especially when the accuracy attains the level of $\sim$kHz, such measurements could be used to make the King plots used to probe new physics~\cite{King_1963,Counts_2020}. This is therefore of particular interest for the narrow intercombination and optical clock transitions~\cite{Miyake_2019}.

The developed atomic cadmium beam itself has the potential to act as a tool for further, precision spectroscopic experiments, or to act as the atomic source for e.g. loading into a magneto-optical trap. In comparison to the vapor-based loading employed elsewhere~\cite{Yamaguchi_2019,Kaneda_2016,Brickman_2007}, this has the advantage of allowing for an improved vacuum in the trap region and improved signal-to-noise from fluorescence measurements due to reduced ambient scattering.

\section*{Funding}
This work has been supported by the European Research Council, Grant No.772126 (TICTOCGRAV).

\section*{Acknowledgments}
We thank Jacopo Catani and Stefan Truppe for fruitful discussions. We thank Gabriele Santambrogio for the loaning of the NIR wavemeter employed in the absolute frequency measurements. We thank Fabio Corti for work on electronic systems and Tommaso Rossetti for assistance with data collection and analysis. This work has been supported by the European Research Council, Grant No.772126 (TICTOCGRAV). J.-P.P acknowledges the support of the Jenny and Antti Wihuri Foundation, the Walter Ahlstr\"om Foundation and the Finnish Foundation for Technology Promotion

\section*{Disclosures}

The authors declare no conflicts of interest.


\section{Introduction}

In the first paper of this series \citep{Mac21a} spacetimes with a
local spatial rotational invariance were (re)-investigated. Here the
corresponding issues for local boost invariance are studied. The issues
arising and the methods to be used are set out in Sections 1 and 2 of
the previous paper. Only those key points required to make this paper
reasonably self-contained will be repeated here.

In studying local rotational symmetry of spacetime (LRS), \cite{Ell67}
introduced three definitions which he showed to be equivalent for dust
spacetimes. The results were later extended to spacetimes with a
perfect fluid and electromagnetic field \citep{SteEll68}. Here Ellis's
definition (A${}_m$) will be studied. It reads:
\begin{quote}
  (A${}_m$) At each point $P$ in an open neighbourhood
$U$ of a point $P_o$, there exists a nondiscrete subgroup $g$ of the
Lorentz group in the tangent space $T_P$ which leaves invariant the
curvature tensor and all its covariant derivatives to the $m$-th
order.
\end{quote}
Implicit in this definition are conditions on the smoothness of the
manifold and the correspondence between the $g$ at separated
points. This paper considers the case where the group $g$ of linear
isotropies contains boosts. This is called local boost invariance
(LBI), `local' meaning that the same $g$ applies throughout $U$.
`Spacetime' here just means a four-dimensional Lorentzian
manifold. The field equations of general relativity will not be used,
but in cases where the Ricci tensor takes the form that would be
implied by specific matter content in general relativity that
interpretation is referred to.

The starting conjecture (based on a claim by \cite{Sik76} now known to
be false in general) is that (A${}_1$) is sufficent to imply Ellis's
definition (C), i.e.
\begin{quote}There exists a local group of motions $G_r$ in an open
neighbourhood $W$ of a point $P_o$ which is multiply transitive on some
$q$[-dimensional] surface through each point $P$ of $W$.
\end{quote}
With LBI this would imply that the spacetime had local
boost symmetry (LBS), i.e.\ that the group $G_r$ contained, for every
$P\in W$, a subgroup of boosts leaving $P$ fixed. A theorem of
\cite{Hal89} implies, under appropriate topological and smoothness
conditions, that $W$ is a region of a manifold in which the same $G_r$
acts globally.

Only Petrov type D and conformally flat spacetimes have a Weyl tensor
that can satisfy (A${}_0$) (or (A${}_m$) for larger $m$) with a group
$g$ containing boosts.  There are rather more Ricci tensor types that
can satisfy (A${}_0$) for boosts. For (A${}_0$) to apply to the whole
Riemann tensor, the Weyl and Ricci tensors must of course be
appropriately aligned, and there may be a nonzero Ricci scalar.
The tracefree part of the Ricci tensor can be characterized by its
Segre type. The possible invariance groups of the Ricci tensor were
listed by Segre type in Table 5.2 of \cite{SteKraMac03}. Table 1 here
lists those with nontrivial invariance groups $\hat{I}_0$ which include
a boost.
\begin{table}[ht] \tabcolsep0.7ex
  \centering\caption{Nontrivial invariance groups containing boosts,
    by Ricci tensor type.}
  \vspace{2ex}
\label{tab:5.2}
\begin{tabular}{@{}ll@{}} \noalign{\vspace{3\p@ plus2\p@ minus1\p@}} \hline\hline \noalign{\vspace{3\p@ plus2\p@ minus1\p@}}
Invariance group & Segre type of the Ricci tensor \\ \noalign{\vspace{3\p@ plus2\p@ minus1\p@}}\hline \noalign{\vspace{3\p@ plus2\p@ minus1\p@}}
\noalign{\vspace{3\p@ plus2\p@ minus1\p@}}
Boosts (only) & $[11(1,1)]$ \\
Boosts and rotations & $[(11)(1,1)]$ \\
$SO(2,1)$: three-dimensional Lorentz group & $[1(11,1)]$ \\
Full Lorentz group & $[(111,1)]$ \\ \noalign{\vspace{3\p@ plus2\p@ minus1\p@}}
\hline\hline
\end{tabular}
\end{table}
\cite{CahDef68} showed that for Petrov type D spacetimes with boost
(or spatial rotation) invariance and any compatible Ricci tensor,
(A${}_2$) was a sufficient criterion for the spacetime to be LRS or
LBS. Subsequently \cite{GooWai86} gave criteria for the LRS Petrov
type D case in terms of the spin coefficients and curvature expressed
in a Newman-Penrose (NP) null tetrad. These criteria were shown in
\cite{Mac21a} to be equivalent to (A${}_1$). The discrepancy with
Cahen and Defrise's use of (A${}_2$) is shown in the Appendix of
\cite{Mac21a} to be due to a less suitable choice of frame in the
calculations.

Here it is shown in Section \ref{sec:LBS}, by arguments parallel to those
of \cite{GooWai86}, that (A${}_1$) is also sufficient for LBS in
Petrov type D spacetimes. In Section \ref{CFcases}, the corresponding
question for conformally flat spacetimes is studied. In both these
sections the detailed arguments are closely related to those of
\cite{Mac21a} by the asterisk operation of the GHP formalism
\citep{GerHelPen73}. A recent preprint addresses local boost
invariance of higher dimensional manifolds \citep{McNColWyl19}.

As in the previous paper, the Cartan-Karlhede procedure\ for characterizing spacetimes and
testing their equivalence, as outlined in Section 2 of \cite{Mac21a},
is used. It relies on the computation of ``Cartan invariants'', the
components of the Riemann tensor and its covariant derivatives in
canonically chosen frames.  The implementation used here employs the
Newman-Penrose formalism as set out in Chapter 7 of
\cite{SteKraMac03}. The ``Newman-Penrose equations'' (the Ricci
equations) and Bianchi identities [(7.21a)-(7.21r) and (7.32a)-(7.32k)
  in \cite{SteKraMac03}] will be referred to below as (NPa)--(NPr) and
(Ba)--(Bk).

A minimal set of Cartan invariants sufficient for the above procedure,
was defined by \cite{MacAma86}.  It consists of totally symmetrized
spinor derivatives of the Newman-Penrose curvature quantities. Here
the shorthand notation for such spinors, as defined in \cite{Mac21a},
will be used. If $Q^{ABC\ldots}{}_{E'F'\ldots}$ is a relevant
curvature quantity then the notation $Q_{AB'}$ denotes the component
of $Q^{(ABC\ldots)}{}_{(E'F'\ldots)}$ in which $A$ of the $m$ unprimed
indices and $B$ of the $n$ primed indices are contracted with the
basis spinors $\iota$ and $\bar{\iota}$ respectively (and the others
with the basis spinors $o$ and $\bar{o}$). $\chi$ is said to have
valence $(m,\,n)$.  The set defined in \cite{MacAma86} consists of the
totally symmetrized derivatives of $\Psi$, $\Phi$ and $\Lambda$, together
with, at order 1, $\Xi_{DEFW'}={\nabla^C}_{W'}\Psi_{CDEF}$ and at
order $q+2$, the d'Alembertians of quantities at order $q$. For a
totally-symmetrized spinor of valence $(m,\,n)$, only components with
$2(A+B)=m+n$ are LBI.

\section{Petrov type D spacetimes with local boost invariance}
\label{sec:LBS}
In the calculations, the boost invariance is assumed to act in the
$(\bm{k},\,\bm{l})$ plane of a Newman-Penrose tetrad adapted to the Petrov
type D Weyl tensor and to leave the Riemann tensor and its first
derivative unchanged\footnote{The results however are independent of
  the choice of frame: they would simply be more difficult to check in
  a randomly chosen frame.}. $\Lambda$ and $\Psi_2 (\neq 0)$ are boost
invariant, so from (A${}_1$) $D\Lambda = \Delta \Lambda =
D\Psi_2=\Delta\Psi_2=0$. In $\Phi_{AB'}$ only the components
$\Phi_{11'}$ and $\Phi_{02'}$ can be nonzero. Boost invariance of the
Cartan invariants $\nabla{\Psi_{AB'}}$ requires that
 $$\nabla{\Psi_{01'}} =  \nabla{\Psi_{11'}} = \nabla{\Psi_{20'}} =
\nabla{\Psi_{31'}} = \nabla{\Psi_{40'}} =\nabla{\Psi_{41'}} =0 $$
and thus
\begin{equation}\label{C1dagvalues}
  \kappa = \sigma = \rho = \mu = \lambda = \nu =0.
\end{equation}
These spacetimes are members of Kundt's class.

Therefore, in Petrov type D, invariance of the Riemann tensor and its
first derivatives under a boost implies that there is a Newman-Penrose
tetrad (a canonical one for Petrov type D, fixed up to a spatial
rotation and boost) in which the following criteria hold.
\begin{eqnarray}
  {\rm (C1}{}^\dag): &&  \kappa=\sigma=\rho=\mu=\lambda=\nu=0.\\
  {\rm (C2}{}^\dag): &&  \Phi_{01}=\Phi_{12}=0.\\
  {\rm (C3}{}^\dag): &&  \Delta \Lambda= D\Lambda = 0.
\end{eqnarray}
{}[The (Cn${}^\dag$) notation is adopted to emphasize the parallel with
the conditions (C1)--(C3) of \cite{GooWai86}.]
The theorem analogous to Theorem 2.1 of Goode and Wainwright is:
\begin{theorem} \label{LBSThm1}
A space-time $($assumed conformally curved$)$ is LBS if
and only if there exists a null tetrad $(\bm{k},\,
  \bm{l},\,\bm{m},\,\bm{\overline{m}})$ in which
{\rm (C1${}^\dag)-($C3${}^\dag$)} hold.
\end{theorem}

To show that conversely (C1${}^\dag$)--(C3${}^\dag$) imply that the
Riemann tensor and its first derivatives are boost invariant, one
first inserts (C1${}^\dag$) into (NPa, b, j, k, m and n), and, using
(C2${}^\dag$), obtains that
$$\Phi_{00'}=\Phi_{22'} =0, \qquad \Psi_0=\Psi_1=\Psi_3=\Psi_4=0,$$ so
that the Riemann tensor is Petrov type D or conformally flat and is
boost invariant. (Be) and (Bf) give $D\Psi_2=\Delta\Psi_2=0$. From
above, (C1${}^\dag$) and (C3${}^\dag$) then guarantee that $\nabla
\Lambda$ and $\nabla \Psi_{AB'}$ are boost invariant.

Using (C1${}^\dag$)--(C3${}^\dag$) in (Bi) and (Bk) gives $D\Phi_{11'}
=\Delta\Phi_{11'}=0$, so the gradient of $\Phi_{11'}$ is boost
invariant.  The components of $\nabla\Phi$ that must vanish for LRSI
are given in full as \eqref{DPhi01B}--\eqref{DPhi23B} in Section
\ref{CFcases}, and the relations between these components of
$\nabla\Phi$ and the Bianchi identities are discussed there. Given
(C1), they all vanish, due to (Bb) and (Bc), so $\nabla\Phi$ is boost
invariant.

To complete the check of the equivalence of the conditions
(C1${}^\dag$)-(C3${}^\dag$) with the assumption that the Riemann
tensor and its first derivatives are boost invariant, one has to show
that once the remaining frame freedom, a spatial rotation, has been
fixed, so that $\Phi_{02'}$ is an invariant,
$D\Phi_{02'}=\Delta\Phi_{02'} =0$, which follows if $\varepsilon$ and
$\gamma$ are real.  (The boost invariance of $\Xi_{AB'}$ is readily
checked.)

These and other restrictions on the spin coefficients analogous to
those in Section 3 of \cite{Mac21a} are now sought, following
analogous steps in \cite{GooWai86}.  They will enable the LBI of
higher derivatives of the curvature to be checked.  From (Bh) and (Bj)
one finds
$$ \delta(-\Psi_2-\Phi_{11'}+\Lambda) = -3\tau\Psi_2 + 2\bar{\pi}\Phi_{11}
+\pi\Phi_{02'}.$$
Applying the $[\bar{\delta},\, D]$ commutator to
$(-\Psi_2-\Phi_{11'}+\Lambda)$ and using (NPc) [which tells us that
  $D\tau=(\varepsilon-\bar{\varepsilon})\tau$], (Bb) and $\Psi_2\neq
0$ one obtains $D\pi= -(\varepsilon-\bar{\varepsilon})\pi$.

Similarly, from (Bg) and the conjugate of (Bj) one obtains
$$ \bar{\delta}(-\Psi_2-\Phi_{11'}+\Lambda) = 3\pi\Psi_2 -
2\bar{\tau}\Phi_{11} - \tau\Phi_{20'}.$$ Applying the
commutator
$[\bar{\delta},\,\Delta]$ to $(-\Psi_2-\Phi_{11'}+\Lambda)$, using (NPi)
[which tells us $\Delta \pi = -(\gamma-\bar{\gamma})\pi$], (Bc) and
$\Psi_2 \neq 0$ one obtains $\Delta\tau = (\gamma-\bar{\gamma})\tau$.

Thus from (NPc), (NPi), and the arguments above:
\begin{eqnarray}\label{dpitau1}
  D\pi= -(\varepsilon-\bar{\varepsilon})\pi, &\quad&
 \quad \Delta\tau =
 (\gamma-\bar{\gamma})\tau \\
 \label{dpitau2}
 D\tau=(\varepsilon-\bar{\varepsilon})\tau, &\quad& \Delta \pi =
-(\gamma-\bar{\gamma})\pi.
\end{eqnarray}
Note that to arrive at \eqref{dpitau1} in conformally flat cases
one would need to derive it by a different argument.

It will now be shown that a position-dependent rotation (the
remaining frame freedom) can always be used to achieve
\begin{equation}\label{aplusb2}
  \varepsilon=\bar{\varepsilon} \quad {\rm and} \quad
  \gamma=\bar{\gamma},
\end{equation}
and hence that
\begin{equation} \label{DDeltas}
  D\Phi_{02'} = \Delta \Phi_{02'}=D\pi = \Delta\pi=D\tau=\Delta\tau=0.
\end{equation}

For non-zero $\tau$ and $\pi$, if either $\pi$ or $\tau$ is real
\eqref{DDeltas} follows immediately from
\eqref{dpitau1}--\eqref{dpitau2}.  If $\pi$ and $\tau$ are both
non-zero $D(\ln(\tau/\pi))=2(\varepsilon-\bar{\varepsilon})$ and
$\Delta(\ln(\tau/\pi))= 2(\gamma-\bar{\gamma})$. If then $\tau/\pi$ is
real, these imply \eqref{aplusb2}. Under a rotation through an angle
$\theta$,
$(\varepsilon-\bar{\varepsilon})^\star=(\varepsilon-\bar{\varepsilon})+iD\theta$
and $(\gamma-\bar{\gamma})^\star=(\gamma-\bar{\gamma})+i\Delta\theta.$
Hence if $\tau/\pi$ is not real, $\theta=\smfrac{1}{2} {\rm Im}(\ln
(\pi/\tau))$ achieves \eqref{aplusb2}.
If just one of $\pi$ or $\tau$ is nonzero, its argument gives a
suitable $\theta$. If $\pi=0=\tau$ one needs to show that the equations
\begin{equation}\label{thetaeqs}
  iD\theta = -(\varepsilon -\bar{\varepsilon})\quad {\rm and} \quad
  i\Delta\theta = -(\gamma -\bar{\gamma}),
\end{equation}
are compatible. This is done by applying the  $[\Delta,\,D]$
commutator to $\theta$. Note that (NPh) implies that
$\Psi_2$ is real in this case; the imaginary part of (NPf) then
shows that the two equations for $\theta$ are indeed compatible.

Finally, (NPe) and the complex conjugate of (NPd), and (NPr)
and the complex conjugate of (NPo), are used to obtain
\begin{equation}\label{aminusb}
  D(\alpha - \bar{\beta}) = 0 = \Delta(\alpha-\bar{\beta}).
\end{equation}
Inserting these results into the higher derivatives, using CLASSI,
shows that the second, third and fourth derivatives of the Riemann
tensor are also boost invariant. This will also be true in conformally
flat spacetimes if \eqref{C1dagvalues}, \eqref{aminusb} and
\eqref{DDeltas} hold, but in the case studied in Section
\ref{SO21case} (A${}_3$) is required to establish \eqref{DDeltas}.

One can show, as in the LRS case, that in Petrov type D one cannot
have $s=1$ and $t_2=1$, $t_1=t_0=0$ which would require (A${}_4$) to
be checked. To eliminate the
possibility, the calculations follow a similar logic to those in
\cite{Mac21a}. Necessarily $\Psi_2 \neq 0$, and $\Phi_{11'}$ and
$\Phi_{02'}$ would be constant. That $\nabla\Psi_{AB'}$ is constant
implies $\tau$ and $\pi$ are constant. If at least one of them is
nonzero, then (NPg) and/or (NPp) imply that $\bar{\alpha}-\beta$ is
constant (possibly zero). Direct calculation (using CLASSI) then shows
$t_2=0$ i.e.\ all terms in the second derivatives are also constant,
so the Cartan-Karlhede procedure\ terminates.  If both $\pi$ and $\tau$ are zero, (NPg)
implies $\Phi_{02'}=0$, and then inspection (using CLASSI) shows that
all first derivatives of the Riemann tensor, and hence all higher
derivatives, are zero, and the Cartan-Karlhede procedure\ terminates at step 1.

Thus for LBS Petrov type D spacetimes it is sufficient to check
(A${}_3$) and the Cartan-Karlhede procedure\ must terminate at the third step or
earlier. This proves the following.

\begin{theorem} \label{LBSDthm2}
  If a spacetime of Petrov type D is such that the Riemann tensor and
  its first derivative are invariant under a local boost invariance,
  then the spacetime is locally boost symmetric and admits a local
  isometry group $G_r~(r\geq 3)$.
\end{theorem}

The converse of Theorem \ref{LBSDthm2} is obvious, and by the
equivalence shown above this proves Theorem \ref{LBSThm1}. Note that
as in Section 3 of \cite{Mac21a} the invariance of the derivatives of
the Ricci tensor has not been used to derive the results, only
checked, and $\Psi_2 \neq 0$ was used only in deriving (C1${}^\dag$)
and \eqref{dpitau1}.  The LBS conclusion depends only on
\eqref{C1dagvalues}, \eqref{aminusb} and \eqref{DDeltas}.

In the following section it is shown that Theorem \ref{LBSDthm2} is
still true with `Petrov type D' replaced by 'conformally flat', unless
the Ricci tensor is of Segre type [1(11,1)] when some cases require
LBI of the curvature and its first three derivatives to ensure
LBS.

\section{ Conformally flat spacetimes with local boost invariance}
\label{CFcases}

The conformally flat cases to be considered are those Ricci tensor
types appearing in Table 1.  (A${}_m$) is assumed to hold with a group
$g$ which contains a boost. By the same argument as in \cite{Mac21a},
Ricci tensors of Segre type [(111,1)] are easily disposed of: the
spacetimes are of constant curvature, the subgroup $g$ in (A${}_m$) is
the trivial one comprising the whole Lorentz group, $s=6$,
$t_p=0=t_0$, and there is a group $G_{10}$ transitive on the whole
spacetime. The Cartan-Karlhede procedure\ terminates at the first step and (A${}_0$)
suffices because it will imply (A${}_1$).

In the rest of this section the actual (A${}_m$) required for local
LBS in conformally flat spacetimes with the various Ricci tensors
which admit a boost invariance, but are of less symmetry than Segre type
[(111,1)], are studied.

\subsection{The first derivatives and Bianchi identities}
\label{CF1stDerivs}

The first step is to impose LBI on the first derivatives of $\Phi$ and
$\Lambda$. Then one can try to derive (C1${}^\dag$), which were
obtained in Petrov type D cases from invariance of $\nabla\Psi$, and
look for an appropriate variant of the rest of the arguments in Section
\ref{sec:LBS}.

With LBI of the curvature only $\Lambda$, $\Phi_{11'}$ and
$\Phi_{02'}=\overline{\Phi_{20'}}$ can be nonzero components in a
suitable canonically chosen frame, and since $\Phi_{11'}$ and
$\Lambda$ are invariant under the remaining allowed changes of frame,
(A${}_1$) requires that $\Delta\Phi_{11'}=D\Phi_{11'}=\Delta \Lambda =
D\Lambda = 0$. (Note that the invariance of curvature is being assumed
here, rather than $\Phi_{00'}=\Psi_{22'}=0$ being deduced from other
assumptions as in Section \ref{sec:LBS}.) The terms in $\nabla\Phi$
that must vanish if (A${}_1$) holds are
\begin{eqnarray}
3\nabla\Phi_{01'}&=&
2(2\kappa\Phi_{11'}+\bar{\kappa}\Phi_{02'}), \label{DPhi01B}\\
3\nabla\Phi_{02'}&=&
D\Phi_{02'}+2(\bar{\rho}-\varepsilon+\bar{\varepsilon})\Phi_{02'}
+4\sigma\Phi_{11'},\label{DPhi02B}\\
9\nabla\Phi_{11'}/2&=& 2(\rho+\bar{\rho})\Phi_{11'}+ \sigma\Phi_{20'}+\bar{\sigma}\Phi_{02'},\label{DPhi11B}\\
  3\nabla\Phi_{13'}&=&\Delta\Phi_{02'} + 2(\bar{\gamma}-\gamma-\mu)\Phi_{02'}-4\bar{\lambda}\Phi_{11'},\label{DPhi13B}\\
  9\nabla\Phi_{22'}/2&=&-2(\mu+\bar{\mu})\Phi_{11'}-\lambda\Phi_{02}
-\bar{\lambda}\Phi_{20'}.\label{DPhi22B}\\
3\nabla\Phi_{23'}&=&
-2(\nu\Phi_{02'}+2\bar{\nu}\Phi_{11'}).\label{DPhi23B}
\end{eqnarray}
In addition $\nabla\Phi_{00'} = \nabla\Phi_{33'}\equiv 0$ here.  Of
the above equations, \eqref{DPhi01B} is equivalent to (Ba),
\eqref{DPhi11B} to (Bi) or the real part of (Be), \eqref{DPhi22B} to
(Bk) or the real part of (Bf), and \eqref{DPhi23B} to the conjugate of
(Bd). The Bianchi identities (Be) and (Bf) give
\begin{equation}\label{BeAndBf}
  \bar{\sigma}\Phi_{02'}+2\rho\Phi_{11'}=0,\quad {\rm and} \quad
  \bar{\lambda}\Phi_{20'}+2\mu\Phi_{11'}=0.
\end{equation}
In the general case the only information from boost invariance
of $\nabla\Phi$ additional to that in the Bianchi identities comes
from \eqref{DPhi02B} and (Bb), and \eqref{DPhi13B} and the conjugate
of (Bc), which give respectively.
\begin{equation}
  \label{DPhiRemB}
  \bar{\rho}\Phi_{02'}+2\sigma\Phi_{11'}=0,\quad {\rm and} \quad
  \bar{\mu}\Phi_{20'}+2\lambda\Phi_{11'}=0.
\end{equation}
If \eqref{DPhiRemB} is satisfied the first derivatives of the
curvature are boost invariant. (In the special case of Segre type
[1(11,1)], where $D\Phi_{02'}=\Delta\Phi_{02'}=0$, \eqref{DPhiRemB} is
equivalent to (Bb) and (Bc), and contains no additional
information. This case is considered in Section \ref{SO21case} below.)

The Bianchi identities (Bg) and (Bh) give
\begin{equation}\label{pitau}
  2(\bar{\pi}+\tau)\Phi_{11'}+(\pi +\bar{\tau})\Phi_{02'}=0.
\end{equation}
The remaining information in the Bianchi identities is
\begin{eqnarray}\label{CFBg}
  Bg: & \delta\Phi_{20'} -2\bar{\delta}\Lambda=&
  (2\bar{\alpha}-2\beta-\bar{\pi})\Phi_{20'}-2\pi\Phi_{11'},\\
  \label{CFBj}
  Bj:& -\delta\Phi_{11'}-\bar{\delta}\Phi_{02'}+3\delta \Lambda = &
  (-2\alpha+2\bar{\beta}+\pi-\bar{\tau})\Phi_{02'}\\
  && \phantom{-2\alpha}+2(\bar{\pi}-\tau)\Phi_{11'}. \nonumber
\end{eqnarray}

\subsection{Ricci tensors of Segre type [11(1,1)]}
\label{ssec:CFG1}

If such spacetimes obey (A${}_\infty$) the Cartan-Karlhede procedure\ must terminate after
at most 3 steps since $s=1$ at every step and at most two steps are
needed to increase $t_p$ to the maximum of 2. So (A${}_3$) would
suffice. As shown next, only (A${}_1$) is actually required. All the
resulting spacetimes admit a $G_r,~r\geq 3$ acting on submanifolds of
dimension at least 2.

The ratio $\Phi_{02'}/\Phi_{11'}=2c$ is fixed under the remaining
frame freedom but may be position-dependent. (Using the remaining
rotational freedom one could make $\Phi_{02'}$ real and positive, and
$c$ real, but this freedom may be needed as in Section \ref{sec:LBS}
to obtain a frame in which \eqref{aplusb2} holds, using
\eqref{dpitau1}--\eqref{dpitau2}.) One must have $|c| \neq 1$ or the
Segre type will be [1(11,1)]. Relations of the form
$2\Phi_{11'}Q+\Phi_{02'}\bar{Q}=0$ imply $Q=0$, since otherwise
$|Q/\bar{Q}|=1\neq |c|=|\Phi_{02'}/\Phi_{11'}|$ (If $|c|=1$ at some
isolated points, continuity implies $Q=0$ there too.) Hence from (Ba)
and (Bd), $\kappa=\nu=0$, and from \eqref{DPhiRemB} and
\eqref{BeAndBf} $\rho=\mu=\sigma=\lambda=0$. So (C1${}^\dag$)
holds. In addition, \eqref{pitau} implies $\bar{\pi}+\tau=0$.

\eqref{dpitau2} holds (being just (NPc) and (NPi)) and \eqref{dpitau1}
then follows simply from it, using $\bar{\pi}+\tau=0$. One can then
obtain \eqref{aplusb2} and \eqref{DDeltas}, and complete the proof
that (A${}_1$) is sufficient to imply LBS in this case as in Section
\ref{sec:LBS}. One may note that since
$D\Phi_{02'}=\Delta\Phi_{02'}=0$, $c$ is constant in the timelike
two-planes determined by the boost.

\subsection{Ricci tensors of Segre type [(11)(1,1)]}

This Segre type is that of the Ricci tensor of a non-null
electromagnetic field.  Here $s_0=2$ and in a canonical frame
$\Phi_{02'}=0$.  Assuming $\Phi_{11'}\neq 0$, (Ba)--(Bd) give
$$\kappa=\sigma=\lambda=\nu=0.$$ If $s=2$ one must have $t_1=0$ and so
the Cartan-Karlhede procedure\ terminates at step 1 since neither $s$ nor $t$ has changed
and as in Section 3 of \cite{Mac21a} this gives the Bertotti-Robinson type
solutions with a $G_6$ transitive on the whole spacetime, and
(A${}_1$) suffices.

For $s=1$ with $\hat{I}_1$ consisting of the boosts,
$D\Lambda=D\Phi_{11'}=\Delta \Lambda = \Delta\Phi_{11'}=0$ and from
(Be) and (Bf) $\rho=\mu=0$ so the conditions of Theorem \ref{LBSThm1}
hold.  As in the previous subsection, one has $\tau+\bar{\pi}=0$, so
the proof that (A${}_1$) with a one-dimensional group of boosts
implies LBS proceeds as in section \ref{sec:LBS}.

\subsection{Ricci tensors of Segre type [1(11,1)]}
\label{SO21case}

Here the invariance group is SO(2,1), generated by null rotations
about $\bm{k}$ and about $\bm{l}$ and a boost in the $(\bm{k},\,\bm{l})$ plane.  The
SO(2,1) group acts in a hyperplane and leaves invariant one direction
in the $(\bm{m},\,\bm{\overline{m}})$ plane. The Ricci tensor represents a
tachyonic fluid, and a canonical form for it which is manifestly null
rotation invariant (as in \cite{Mac20}) about each of the null
directions has only $\Phi_{11'}$ and $\Phi_{02'}$ non-zero with
$2|\Phi_{11'}|=|\Phi_{02'}|$. (This is a specialization of the form
for Segre type [11(1,1)], treated above.) Using the remaining freedom
of spatial rotation in the $(\bm{m},\,\bm{\overline{m}})$ plane one can set
$\Phi_{02'} =2\Phi_{11'}$: the parameters of both null rotations are
then pure imaginary and the vector orthogonal to the hyperplane in
which the SO(2,1) acts is in the direction $\bm{m}+\bm{\overline{m}}$.

The boost rescales the parameters of the null rotations. One might
therefore have invariance under a two-dimensional subgroup of SO(2,1)
generated by the boost and one of the null rotations (see entry R6 in
Table 6.1 of \cite{Hal04}). Among the quantities defined by
\cite{MacAma86}, only $\nabla^k\Phi_{AB'}$ and
$\nabla^k\Lambda_{AB'}$, both of which are Hermitian, and
d'Alembertians thereof, have to be considered. Boost invariance
implies that of the components of $\nabla^k\Phi_{AB'}$, only those
with $A+B=2+k$ can be nonzero. If then $\nabla^k\Phi_{AB'}$ is
invariant under one of the null rotations, then from the Hermitian
symmetry it will also be invariant under the other. Thus
$\nabla^k\Phi_{AB'}$ will be SO(2,1) invariant. The same applies to
$\nabla^k\Lambda_{AB'}$ and the d'Alembertians of spinors of lower
derivative order.

So $s=2$ is impossible and either $s=3$ or $s=1$. If $s=3$ there can
be at most one independent function of position (as there is just one
spacelike direction fixed under $g$). The spacetime admits a $G_6$ (or
in special cases a $G_7$, cf.\ \cite{RebTei92}) acting on timelike
hypersurfaces of constant curvature. The metrics include analogues of
the FLRW metrics for perfect fluids.

The remaining case is where $s=1$ and $g$ in (A${}_m$) is just the boost
invariance\footnote{It might be convenient to carry out the necessary
  calculations using the 3+1 orthonormal tetrad formalism based on a
  spacelike congruence introduced by \cite{Har82}.}. From (Ba) and
(Bd)--(Bf) one has
\begin{equation}\label{CFtacBian}
\kappa+\bar{\kappa}=\nu+\bar{\nu}=\bar{\mu}+\lambda =\bar{\rho}+\sigma=0,
\end{equation}
and then (Bb) and (Bc) give \eqref{aplusb2}.  In
this case the Bianchi identities and
$D\Phi_{11'}=\Delta\Phi_{11'}=D\Lambda=\Delta \Lambda = 0$ ensure (A${}_1$),
\eqref{DPhiRemB} being equivalent to \eqref{BeAndBf} in this case.
From \eqref{pitau} one has
\begin{equation}\label{ptident}
  (\pi+\bar{\pi})+ (\tau+\bar{\tau})=0.
\end{equation}

The real parts of the left sides of \eqref{CFBg} and \eqref{CFBj} are
proportional to the real part of
$\delta(\Phi_{11'}-\Lambda)$.
For the right sides to be
compatible one must have
\begin{equation}
  \label{abpident} 0=(\alpha+\bar{\alpha})-(\beta+\bar{\beta})+(\pi+\bar{\pi}).
\end{equation}

Eliminating between the imaginary parts of \eqref{CFBg} and \eqref{CFBj} yields
\begin{eqnarray}
  \label{delR}(\delta-\bar{\delta})\Lambda&=&[(\bar{\alpha}-\beta)-(\alpha-\bar{\beta})]\Phi_{11'},  \\
  \label{delPhi}
  (\delta-\bar{\delta})\Phi_{11'}&=&[(\bar{\alpha}-\beta)-(\alpha-\bar{\beta})]\Phi_{11'}.
\end{eqnarray}
The perhaps surprising equality of the right sides does not imply that
$\delta(\Phi_{11'}-\Lambda)=0$.

The real parts of (NPc) and (NPi) give $D (\tau+\bar{\tau})=0=
\Delta(\pi+\bar{\pi})$, so
\begin{equation}\label{dpt}
  D (\tau+\bar{\tau})=\Delta(\pi+\bar{\pi}) = \Delta
  (\tau+\bar{\tau})=D(\pi+\bar{\pi})=0.
\end{equation}
From the imaginary part of (NPg) minus the conjugate of (NPh), using
\eqref{CFtacBian} and \eqref{abpident},
\begin{equation}\label{delpi}
  \delta(\pi+\bar{\pi})=  \bar{\delta}(\pi+\bar{\pi}).
\end{equation}

From the conjugate of (NPd) together with (NPe), and from the
conjugate of (NPr) with (NPo), one finds using \eqref{abpident} that
\begin{equation}
  D(\bar{\alpha}-\beta)= \Delta(\bar{\alpha}-\beta)= 0,
\end{equation}
which is consistent with \eqref{abpident} and \eqref{dpt}.

From (NPb) and the complex conjugate of (NPa), and from the conjugate
of (NPj) and (NPn), one obtains
\begin{equation}\label{kident}
  \kappa(2\bar{\alpha}-2\beta+(\pi + \bar{\pi}))=0=
  \nu(2\bar{\alpha}-2\beta+(\pi + \bar{\pi})).
\end{equation}
Since $\kappa$ and $\nu$ are pure imaginary, the imaginary parts of
these equations follow from \eqref{abpident} while the real parts give
$\kappa\Im_{\alpha\beta}=\nu\Im_{\alpha\beta}=0$, where
$\Im_{\alpha\beta} \equiv
(\alpha-\bar{\beta})-(\bar{\alpha}-\beta)$. (From the relations
between null and orthonormal tetrads obtainable using \cite{GooWai86}
and \cite{Har82} one finds that $\Im_{\alpha\beta}$ plays the same
r\^ole here as $\dot{u}_1$ did in the conformally flat perfect fluids
studied in \cite{Mac21a}.) So far only (A${}_1$) has been used.

Calculating $\nabla^2\Phi$ shows that the $00'$, $01'$, $34'$ and
$44'$ components are identically zero. Most of the other components
with $A+B\neq 4$ vanish on simplification with the help of
\eqref{delPhi}, \eqref{ptident}, \eqref{abpident}, and \eqref{kident},
which all follow from (A${}_1$). The remaining components which must
vanish for LBI ($\nabla^2\Phi_{03}$ and $\nabla^2\Phi_{14'}$) yield
$\Im_{\alpha\beta}\sigma= \Im_{\alpha\beta}\mu =0$. For
the spacetime to obey (A${}_2$) thus requires $\Im_{\alpha\beta}$
multiplying each of $\kappa$, $\sigma$, $\mu$ and $\nu$ to be zero.

To satisfy (A${}_2$) there are now two possibilities: either (a)
$\Im_{\alpha\beta}=(\delta-\bar{\delta})\Phi_{11'}
=(\delta-\bar{\delta})\Lambda=0$, showing that $\Phi_{11'}$ and
$\Lambda$ are constant in the hyperplanes defined by the SO(2,1)
action, or (b) $\Im_{\alpha\beta}\neq 0$ and (C1${}^\dag$) holds.

In case (a), we find that $\nabla^k\Phi$ and $\nabla^k\Lambda$ are
invariant under a null rotation about $\bm{k}$ for $k=1\dots3$ and
therefore $s=3$ by the earlier argument that $s_i$ cannot be 2.  We
thus again have at least a $G_6$ if $\Im_{\alpha\beta}=0$. Note that
(A${}_2$) was checked but does not give extra conditions in this case.

In case (b), one now has the conditions (C1${}^\dag$) and
\eqref{aplusb2}. (NPc) and (NPi) give $D\tau=\Delta\pi=0$. To complete
\eqref{DDeltas}, one has to show that $D\pi=\Delta \tau=0$.  Then the
curvature derivatives up to the fourth will have boost invariance, as
in Section \ref{sec:LBS}. (Note that since
$(\delta-\bar{\delta})\Phi_{11'}\neq 0$, $t_0=0$ need not be
considered here.)  One finds $\nabla^3\Phi_{22'}=
-6\Phi_{11'}\Im_{\alpha\beta}D(\pi+\tau)/25$ and $\nabla^3\Phi_{33'}$
gives the same with $\Delta$ replacing $D$. Thus to obtain
\eqref{DDeltas} one needs (A${}_3$) in this case.

\section*{Acknowledgements}
I am grateful to Jan {\AA}man for his work in developing the software
CLASSI and to Filipe Mena, who was kind enough to check and correct some of the
calculations.

\section{Conclusion}

The work in Sections \ref{sec:LBS} and \ref{CFcases} gives
the following analogue of Theorem 3 of \cite{Mac21a}.

\begin{theorem} In spacetimes with a Ricci tensor of Segre type
    $[1(11,1)]$ whose distinguished spacelike eigenvector is not geodesic, local
  boost invariance of the curvature and its derivatives
  up to the third holds if and only if the spacetime is locally
  boost symmetric. In all other cases, local boost
  invariance of the curvature and its first derivatives holds if and
  only if the spacetime is locally boost symmetric.
\end{theorem}


\section{Introduction}\label{intro}
\emph{Mereology} is the study of the relations of part to whole and the relations of part to part within a whole.
\emph{Quantum mereology} studies such issues as the bipartite decomposition of a quantum system into a ``system'' and ``environment'', the relationship of a distinguished ``system'' and an ``observer'', the emergence of space and time from quantum entanglement, and other fundamental questions in quantum mechanics \cite{Carroll,Cao}.
The general scheme is as follows.
The whole is an isolated quantum system%
\footnote{Obviously, in the exact sense, isolated systems do not exist (or they are fundamentally unobservable), with the possible exception of the Universe as a whole.}
in a given pure state undergoing a given unitary (Schr\"odinger) evolution.
Then, in a way chosen according to certain criteria,
the system is decomposed into a tensor product of subsystems.
By reducing the ``universe'' density matrix,  we obtain mixed states for subsystems and can study the energies and entanglement measures associated with subsystems, and their time evolution.
For the corresponding computations, involving rather tedious combinatorics, we develop a model based on a finite formulation of quantum mechanics.
\section{Factorization of a Hilbert Space}\label{factor}
\paragraph{Tensor product of Hilbert spaces}
The (\emph{global}) Hilbert space \Math{\Hspace} of a \Math{K}-component quantum system is the tensor product of the (\emph{local}) Hilbert spaces \Math{\Hspace_k} of the components:
\MathEqLab{\Hspace=\KroneckerProduct_{k=1}^K\Hspace_k\,.}{HspaceN}
If \Math{\dim\Hspace=\mathcal{N}} and \Math{\dim\Hspace_k=d_k}, then \Math{\mathcal{N}=\prod_{k=1}^K{}d_k}.
For any \Math{d}-dimensional Hilbert space, we denote the \Math{i}th orthonormal basis element by \Math{\barket{i}}, that is,
\Math{\barket{0}=\vect{1,0,\ldots}^{\top}},
 \Math{\barket{1}=\vect{0,1,0,\ldots}^{\top}},\Math{\ldots,\barket{d-1}=\vect{0,0,\ldots,1}^{\top}}.
Tensor monomials of local basis elements form an orthonormal basis in the global Hilbert space:
\MathEqLab{\barket{i}=\barket{i_1}\KroneckerProduct\cdots\KroneckerProduct\barket{i_k}\KroneckerProduct\cdots\KroneckerProduct\barket{i_K},}{1to1corresp}
where \Math{\barket{i}\in\Hspace}, \Math{\barket{i_k}\in\Hspace_k} and
\MathEqLab{i=i_1\prod_{m=2}^Kd_m+\ldots+i_k\prod_{m=k+1}^Kd_m+\ldots+i_K.}{1to1numeric}
\paragraph{Tensor decomposition of a Hilbert space}
We can {reverse the procedure}, since \eqref{1to1corresp} is a one-to-one correspondence --- the sequence \Math{i_1,\ldots,i_K} is uniquely recovered from \Math{i} using  formula \eqref{1to1numeric}.
\par
Starting with an orthonormal basis in an \Math{\mathcal{N}}-dimensional Hilbert space \Math{\Hspace} and a decomposition \Math{\mathcal{N}=d_1\cdots{}d_K},
we can construct a particular isomorphism between \Math{\Hspace} and the tensor product of local spaces of the corresponding dimensions.
\par
When constructing an isomorphism, we must take into account the freedom in the choice of bases in Hilbert spaces.
Any two orthonormal bases are related by a unitary transformation.
Using the properties of the tensor product, we can write
\Ali{&U\barket{\psi}=\,U_1\barket{\psi_1}\KroneckerProduct\cdots\KroneckerProduct{}U_K\barket{\psi_K}
=\vect{U_1\KroneckerProduct\cdots\KroneckerProduct{}U_K}\vect{\barket{\psi_1}\KroneckerProduct\cdots\KroneckerProduct\barket{\psi_K}}\\
&\implies\vect{U_1\KroneckerProduct\cdots\KroneckerProduct{}U_K}^{-1}U\barket{\psi}=\barket{\psi_1}\KroneckerProduct\cdots\KroneckerProduct\barket{\psi_K},
}
\!\!where \Math{\barket{\psi}\in\Hspace,\, \barket{\psi_k}\in\Hspace_k;\,U,\,U_k} are unitary transformations in the corresponding spaces.
We see that local transformations can be absorbed by the global transformation \Math{U},%
~so in general we have
\MathEq{U\barket{\psi}=\barket{\psi_1}\KroneckerProduct\cdots\KroneckerProduct\barket{\psi_K}.}
Thus, to specify a factorization of a Hilbert space \Math{\Hspace} we need a \emph{decomposition} of \Math{\dim\Hspace} and a \emph{unitary transformation} that fixes a basis in  \Math{\Hspace}.
\section{Finite Version of Quantum Mechanics}\label{FQM}
We use a version of quantum theory \cite{Kornyak17ISQS,Kornyak17MMCP,Kornyak11MMCP} in which the groups of unitary evolutions are replaced by linear representations of finite groups, and the field of complex numbers is replaced by its dense constructive subfields that naturally arise from the natural numbers and roots of unity.
\paragraph{Permutation Hilbert space}
The fact that \emph{any} linear (hence \emph{unitary}) representation of a finite group is a subrepresentation of some permutation representation implies that the formalism of quantum mechanics can be completely%
\footnote{Modulo empirically insignificant elements of traditional formalism such as infinities of various kinds.} reproduced based on permutations of some set
\MathEqLab{\Omega=\set{e_1,\ldots,e_\mathcal{N}}\cong\set{1,\ldots,\mathcal{N}}}{omega}
of primary (\emph{``ontic''}) objects on which a permutation group \Math{G\leq\SymGN{\mathcal{N}}} acts.
\par
The Hilbert space on the set \Math{\Omega}, necessary for calculating quantum probabilities, can be most economically constructed on the basis of two primitive concepts: (a) \emph{natural numbers}  \Math{\N=\set{0,1,\ldots}}, abstraction of \emph{counting},   and (b) \emph{roots of unity}, abstraction of \emph{periodicity}.
\par
To construct a field \Math{\NF} sufficient for all calculations in the quantum formalism, in particular, for splitting any representation of any subgroup of  \Math{G} into irreducible components, we can proceed as follows.
We extend the semiring \Math{\N} to the ring  \Math{\N\ordset{\runi{\ell}}},
where \Math{\runi{\ell}} is the \Math{\ell}th primitive root of unity, and  \Math{\ell} is the least common multiple of the orders (periods) of the elements of  \Math{G}.
The \emph{algebraic integer} \Math{\runi{\ell}} can be written in complex form as \Math{\displaystyle\runi{\ell}=\e^{{2\pi}\imu/{\ell}}}.
Finally, constructing the \emph{field of fractions} of the ring \Math{\N\ordset{\runi{\ell}}}, we arrive at the \emph{cyclotomic extension} of the rationals \MathEq{\NF=\Q\farg{\e^{{2\pi}\imu/{\ell}}}.}
For \Math{\ell>2}, the field \Math{\NF}, being a dense subfield of \Math{\C}, is physically (empirically) indistinguishable from the field of complex numbers.
\par
Treating the set \Math{\Omega} as a basis, we obtain an \Math{\mathcal{N}}-dimensional Hilbert space \Math{\Hspace_\mathcal{N}} over the field \Math{\NF}.
The action of \Math{G} on \Math{\Omega} determines the \emph{permutation representation} \Math{\mathcal{P}} in  \Math{\Hspace_\mathcal{N}} by the matrices
\MathEq{\mathcal{P}\farg{g}_{i,j}=\delta_{ig,j},}
where \Math{ig} denotes the (right) action of \Math{g\in{}G} on \Math{i\in\Omega}.
\paragraph{Decomposition of permutation representation}
The permutation representation of any group \Math{G} has the \emph{trivial} one-dimensional subrepresentation in the space spanned by the all-ones vector
\MathEq{\barket{\omega}=(\underbrace{1,1,\ldots,1}_{\mathcal{N}})^{\top}.}
The complement to the trivial subrepresentation is called the \emph{standard representation}.
The operator of projection onto the \Math{\vect{\mathcal{N}-1}}-dimensional \emph{standard space} \Math{\Hspace_{\,\!\star}} has the form
\MathEq{\projstd=\idmat_{\mathcal{N}} - \frac{\veco{\omega}{\omega}}{\mathcal{N}}\,.}
Quantum mechanical behavior (interference, etc.) manifests itself precisely in the standard representation.
Banks made a profound observation \cite{Banks} that the projection of classical permutation evolutions in the whole \Math{\Hspace_\mathcal{N}} leads to truly quantum evolutions in the subspace \Math{\Hspace_{\,\!\star}} and demonstrated that the choice \Math{G=\SymGN{\mathcal{N}}}, where \Math{\mathcal{N}} is the number of fundamental (Planck) elements,%
\footnote{The number  \Math{\mathcal{N}} is estimated as \Math
\sim\Exp\farg{\Exp\farg{20}}} and \Math
\sim\Exp\farg{\Exp\farg{123}}} for 1 cm$^3$ of matter and for the entire Universe, respectively.}
makes it possible to reproduce finite dimensional approximations to all known models of theoretical physics.
\paragraph{Quantum states as projections of natural vectors}
\Math{\SymGN{\mathcal{N}}} is a \emph{rational-representation group}, i.e., its every irreducible representation (the standard representation is one of them) is realizable over \Math{\Q}.
This means that to describe evolutions in \Math{\Hspace_{\,\!\star}}, it is sufficient to consider only vectors with rational components.%
\footnote{Complex numbers (nontrivial elements of cyclotomic extensions) can be required only when a representation of a proper subgroup of \Math{\SymGN{\mathcal{N}}} must be split into irreducible components.}
Taking into account that a quantum state is a ray in vector space, it is easy to show that any quantum state in \Math{\Hspace_{\,\!\star}} can be obtained by projection of vectors from \Math{\Hspace_\mathcal{N}} with natural components.
These \emph{natural} vectors are integer points of the nonnegative orthant:
\Math{\barket{n}=\vect{n_1,\ldots,n_\mathcal{N}}^\top\in\N^{\mathcal{N}}\subset\Hspace_\mathcal{N}.}
To build constructive models, one need to select a finite subset in \Math{\N^{\mathcal{N}}}.
Natural vectors whose coordinates belong to the set \Math{\set{0,1,\ldots,m}, m\geq2,} will be called \Math{m}\emph{th order vectors}.
Quantum states described by collinear vectors are equivalent.
To eliminate this redundancy, we can pre-project the natural vectors onto  \Math{\S^{^{\mathcal{N}-1}}_+}\!\!, the part of the unit sphere in the nonnegative orthant.
The total number of quantum states defined by vectors of the \Math{m}th order is not less than
\MathEqLab{2^{\mathcal{N}}-2\,.}{Rm}
The area of \Math{\S^{^{\mathcal{N}-1}}_+} is equal to
\MathEqLab{\frac{\mathcal{N}\pi^{\mathcal{N}/2}}{2^{\mathcal{N}}\Gamma\farg{\mathcal{N}/2+1}}\approx\sqrt{\frac{\mathcal{N}}{\pi}}\vect{\frac{\e\pi}{2\mathcal{N}}}^{\mathcal{N}/2}\!.}{area}
A rough comparison of the exponentially growing --- as can be seen from \eqref{Rm} --- number of quantum states with the rapidly decreasing area \eqref{area} shows that for large \Math{\mathcal{N}}, the \Math{m}th order vectors represent a significant part of quantum states.
\paragraph{Ontic vectors}
We call \Math{2}nd order vectors \emph{ontic vectors}.
These vectors are attractive for both ontological and computational reasons.
The ontic vector \Math{\barket{q}} can be written as a bit string of length \Math{\mathcal{N}}.
Interpreting this string as a \emph{characteristic function}, one can identify the ontic vector with the corresponding nontrivial subset of the set of ontic elements \eqref{omega}:
 \Math{q\subset\Omega}.
The complete set of ontic vectors is \MathEq{Q=2^\Omega\setminus\set{\emptyset,\Omega},~~ \cabs{Q}=2^\mathcal{N}-2\,.}
\indent{}The inner product of ontic vectors \Math{\barket{q}} and \Math{\barket{r}} in the space \Math{\Hspace_\mathcal{N}} has the form \Math{\inner{q}{r}=\popnumb{q\bitand{}r}},
where \Math{\bitand} is the bitwise AND for bit strings, and \Math{\popnumb{\,\cdot\,}} is the number of ones in a bit string (\emph{population number} or \emph{Hamming weight}).
For the inner product of normalized projections of  \Math{\barket{q}} and \Math{\barket{r}} onto \Math{\Hspace_{\,\!\star}}, we have
\MathEq{S\farg{q,r}\equiv\frac{\innerform{q}{\projstd}{r}}{\sqrt{\innerform{q}{\projstd}{q}\innerform{r}{\projstd}{r}}}
=\frac{\mathcal{N}\!\popnumb{q\,\bitand\,r}-\!\popnumb{q}\popnumb{r}}{\sqrt{\popnumb{q}\popnumb{\bitnot\!q}\popnumb{r}\popnumb{\bitnot\!r}}},}
where \Math{\bitnot} denotes bitwise inversion.
The obvious identities
\Math{\popnumb{\bitnot\!a}=\mathcal{N}-\popnumb{a}} and \Math{\popnumb{a\,\bitand\,b}+\popnumb{a\,\bitand\bitnot\!b}=\popnumb{a}}
imply the folowing symmetry with respect to transpositions of subsets of the ontic set and their complements
\MathEq{S\farg{q,r}=-S\farg{\bitnot\!q,r}=-S\farg{q,\bitnot\!r}=S\farg{\bitnot\!q,\bitnot\!r}\,.}
\section{Decomposition of a Quantum System}\label{Decomposition}
Since any mixed state of a quantum system can be obtained from a pure state in a larger Hilbert space by taking a partial trace, it is natural to assume that at a fundamental level the state of an isolated system must be pure.%
\footnote{This belief is expressed by the metaphor ``the Church of the Larger Hilbert Space'' (J.A. Smolin).}
\par
The original permutation basis in the space \Math{\Hspace_{\mathcal{N}}}, i.e., the set \Math{\Omega}, will be referred as the \emph{ontic basis}.
In this basis, the pure density matrix in \Math{\Hspace_{\,\!\star}} associated with the ontic state \Math{\barket{q}\in\Hspace_{\mathcal{N}}} has the form
\MathEqLab{\mathlarger{\rho}_q^o=\frac{\projstd\veco{q}{q}\projstd}{\innerform{q}{\projstd}{q}}=\frac{1}{\mathcal{N}}
\frac{\vect{\barket{q}-\alpha\barket{\omega}}
\vect{\brabar{q}-\alpha\brabar{\omega}}}
{\alpha\vect{1-\alpha}}\,,}{Rho}
where \Math{\displaystyle\alpha=\frac{\popnumb{q}}{\mathcal{N}}} is the \emph{population density}.
There is an obvious duality: the expression for the density matrix \Math{\displaystyle\mathlarger{\rho}_{\bitnot{}q}^o} is obtained from \eqref{Rho} by replacements
\Math{q\rightarrow\,\,\bitnot\hspace{-4pt}q} and \Math{\alpha\rightarrow1-\alpha}.
\paragraph{Energy basis}
In continuous QM, the evolution of an isolated system is described by the one-parameter unitary group \Math{U_t=\e^{-\imu{}Ht}} generated by the Hamiltonian \Math{H} whose eigenvalues are called \emph{energy eigenvalues}.
In finite QM, the evolution is described by a unitary representation of a cyclic group \Math{U\farg{g}^t} generated by an element \Math{g\in{}G}, where \Math{t} is an integer parameter.
In our case, the generator of evolution is some matrix \Math{\mathcal{P}\farg{g}} from the permutation representation of  \Math{G} in \Math{\Hspace_\mathcal{N}}.
We call the \emph{energy basis} an orthonormal basis in which the matrix \Math{\mathcal{P}\farg{g}} is diagonal.
\par
Any permutation can be written as a product of disjoint cycles.
It is easy to show that the total number of cycles of length \Math{\lcycle} in the whole group \Math{\SymGN{\mathcal{N}}} is \Math{\displaystyle\frac{\mathcal{N}!}{\lcycle}},
and, therefore, the expected number of \Math{\lcycle}-cycles in a single permutation is \Math{\displaystyle\frac{1}{\lcycle}}.
That is, high-frequency (high-energy) evolutions are more common.
\par
The \Math{\lcycle}-cycle matrix has the form
\MathEq{C_{\lcycle}=
\begin{bmatrix}
0&1&0&\cdots&0\\
0&0&1&\cdots&0\\
\vdots&\vdots&\vdots&\vdots&\vdots\\
1&0&0&\cdots&0
\end{pmatrix}\,.}
The diagonal form of this matrix is
	\MathEq
	F_{\lcycle}C_{\lcycle}F^{-1}_{\lcycle}=
\begin{bmatrix}
1&0&0&\cdots&0\\
0&\runi{\lcycle}&0&\cdots&0\\
0&0&\runi{\lcycle}^2&\cdots&0\\
\vdots&\vdots&\vdots&\vdots&\vdots\\
0&0&0&\cdots&\runi{\lcycle}^{\lcycle-1}
\end{pmatrix}\,,}
where \Math{\displaystyle\runi{\lcycle}=\e^{{2\pi}\imu/{\lcycle}}} is the \Math{\lcycle}th primitive root of unity, and
\MathEq{F_{\lcycle}=\frac{1}{\sqrt{\lcycle}}
\begin{bmatrix}
1&1&1&\cdots&1\\
1&\runi{\lcycle}^{-1}&\runi{\lcycle}^{-2}&\cdots&\runi{\lcycle}^{-\vect{\lcycle-1}}\\
1&\runi{\lcycle}^{-2}&\runi{\lcycle}^{-4}&\cdots&\runi{\lcycle}^{-2\vect{\lcycle-1}}\\
\vdots&\vdots&\vdots&\cdots&\vdots\\
1&\runi{\lcycle}^{-\vect{\lcycle-1}}&\runi{\lcycle}^{-2\vect{\lcycle-1}}&\cdots&\runi{\lcycle}^{-\vect{\lcycle-1}\vect{\lcycle-1}}
\end{pmatrix}}
is the \emph{Fourier transform} matrix.
\Math{F_{\lcycle}} is both unitary and symmetric, therefore \Math{\displaystyle{}F^{-1}_{\lcycle}=F_{\lcycle}^*.}
In general, the matrix of the permutation representation of an element \Math{g\in{}\SymGN{\mathcal{N}}} is the direct sum of cyclic matrices
\Math{\displaystyle{}\mathcal{P}\farg{g}=\bigoplus\limits_{m=1}^MC_{\lcycle_m}}, and the corresponding diagonalizing matrix is \Math{\displaystyle{}F=\bigoplus\limits_{m=1}^MF_{\lcycle_m}},
which is the transition matrix from the \emph{ontic basis} to the \emph{energy basis}.
The density matrix in the energy basis can be calculated from \eqref{Rho} by the formula
\MathEq{\mathlarger{\rho}_q^e=F\mathlarger{\rho}_q^oF^{-1}.}
\paragraph{Decomposition of a quantum state}
If the basis of a space \Math{\Hspace} is fixed, then the procedure for the decomposition of a quantum system in the state \Math{\mathlarger{\rho}} is reduced to the following:
\begin{enumerate}
	\item The factorization of \Math{\dim\Hspace} is chosen: \Math{\dim\Hspace=d_1\cdot{}\ldots\cdot{}d_K}.
	\item Subsystems are identified with subsets \Math{A} of the set of points \MathEqLab{X=\set{1,\ldots,K}.}{setX}
	~\\[-20pt]
	\item
The basis elements of \Math{\Hspace} are identified with the tensor monomials from the tensor product
\eqref{HspaceN} in accordance with \eqref{1to1corresp} and using formula \eqref{1to1numeric}.
	\item
The density matrix of the subsystem \Math{A} is calculated by taking the trace over the complement:	\Math{\mathlarger{\rho}_A=\tr_{X\setminus{}A}\mathlarger{\rho}}.
\end{enumerate}
Having density matrices of subsystems, it is possible to calculate the physical characteristics describing subsystems, and interactions and quantum correlations between them:
energies, entropies, mutual information and other entanglement measures.
\section{Entanglement Measures}\label{Measures}
Quantitatively, quantum correlations are described by \emph{measures of entanglement}, which are based on the concept of entropy.
The most commonly used in physics is the \emph{von Neumann entropy}
\Math{S_1\farg{\mathlarger{\rho}}=-\tr\farg{\mathlarger{\rho}\log\mathlarger{\rho}}.}
Also often used are entropies from the \emph{R\'enyi family}
\MathEq{S_\alpha\farg{\mathlarger{\rho}}=\frac{1}{1-\alpha}\log\tr\farg{\mathlarger{\rho}^{\alpha}},~~\alpha\geq0,~\alpha\neq1.}
The common feature of the von Neumann and R\'enyi entropies is their additivity on combinations of independent probability distributions determined by the eigenvalues of the density matrices.
The von Neumann entropy is preferred because it additionally satisfies a stronger requirement, the chain rule for conditional entropies.
\par
In our calculations, we use the 2nd R\'enyi entropy (also called the \emph{collision entropy}) \Math{S_2\farg{\mathlarger{\rho}}=-\log\tr\farg{\mathlarger{\rho}^2}} for the following reasons:
\begin{itemize}
	\item
It is easy to calculate: \Math{\tr\farg{\mathlarger{\rho}^2}=\sum\limits_{i=1}^n\mathlarger{\rho}_{ii}^2+2\sum\limits_{i=1}^{n-1}\sum\limits_{j=i+1}^{n}\cabs{\mathlarger{\rho}_{ij}}^2.}
	\item
The value \Math{\tr\farg{\mathlarger{\rho}^2}} coincides with the Born probability: ``the system observes itself.''
	\item
The value \Math{\tr\farg{\mathlarger{\rho}^2}} is the square of the \emph{Frobenius} (\emph{Hilbert-Schmidt})  \emph{norm} of the density matrix.
\end{itemize}
\section{\emph{\emph{Illustrative Calculation}}}\label{Illustrative}
\begin{center}
\includegraphics[width=0.95\textwidth,height=0.3998\textheight]{EntropyPCA21.pdf}
\end{center}
The figure presents the values of entropy \Math{S_2\farg{\mathlarger{\rho}_A}} computed in the ontic basis for the decomposition \Math{\mathcal{N}=2^{12}}.
Data for subsystems of all possible sizes, computed for ten randomly generated ontic vectors, demonstrate the following features:
\begin{itemize}
	\item
Weak dependence on the quantum state: visually, all graphs are almost identical.
Note that this behavior arises for a sufficiently large number of the decomposition components \eqref{setX}.
In this case, \Math{\cabs{X} = 12}.
	\item
Symmetry \Math{S_2\farg{\mathlarger{\rho}_A}=S_2\farg{\mathlarger{\rho}_{X\setminus{}A}}}	is a manifestation of the \emph{Schmidt} bipartite \emph{decomposition} of a pure state:
both matrices \Math{\mathlarger{\rho}_A} and \Math{\mathlarger{\rho}_{X\setminus{}A}} have identical sets of nonzero eigenvalues.
	\item
For \Math{\cabs{A}} noticeably smaller than \Math{\cabs{X}/2}, the reduced state is close to the \emph{maximally mixed} state:	\Math{S_2\farg{\mathlarger{\rho}_{_A}}\approx\cabs{A}\log{}d}.
In our example, \Math{d=2}.
\end{itemize}