Identifying non-contaminated sample of high-redshift galaxies with escaping Lyman continuum (LyC) flux is important for understanding the sources and evolution of cosmic reionization. We present CLAUDS $u$-band photometry of the COSMOS field to probe LyC radiation from spectroscopically confirmed galaxies at $z\geq3.5$ and outside the standard Lyman-break galaxy colour selection expectations. Complementary to the CLAUDS data, we use Subaru multi-filter photometry, \textit {Hubble Space Telescope} ($HST$) multi-filter imaging, and the spectroscopic surveys D10K, VUDS and 3D-HST. We present a sample of Lyman continuum galaxy (LCG) candidates in the redshift range $3.5\lesssim z\lesssim5.1$. Here, we introduce 5 LCG candidates, where two are flagged quality 1 and three quality 2. The estimated $f_{\rm esc}^{\rm abs}$ for quality 1 candidates are in the range $\sim5\% - 73\%$ and $\sim30\% - 93\%$. These estimates are based on our derived parameters from individual galaxies as inputs to a range of BPASS models as well as mean intergalactic medium (IGM) and maximal intergalactic and circumgalactic media (IGM+CGM) transmission. We conclude that our search for LCGs is most likely biased to lines of sight with low HI densities or free from Lyman limit systems. Our two best LCG candidates have EW (Ly$\alpha)\leq50$\AA\ and we find no correlation or anti-correlation between EW (Ly$\alpha$), $f_{\rm esc}^{\rm abs}$, and $R_{\rm obs}$, the ratio of ionizing to non-ionizing observed flux in the measured passbands. Stacking candidates without solid LyC detections ($S/N<3$) results in an estimated $f_{\rm esc}^{\rm abs}$ from galaxies not greater than $1\%$.

We present the first results of our pilot study of 8 photometrically selected Lyman continuum (LyC) emitting galaxy candidates from the COSMOS field and focus on their optical emission line ratios. Observations were performed in the H and K bands using the Multi-Object Spectrometer for Infra-Red Exploration (MOSFIRE) instrument at the Keck Observatory, targeting the [OII], H$\beta$, and [OIII] emission lines. We find that photometrically selected LyC emitting galaxy candidates have high ionization parameters, based on their high [OIII]/[OII] ratios (O32), with an average ratio for our sample of 2.5$\pm$0.2. Preliminary results of our companion Low Resolution Imaging Spectrometer (LRIS) observations, targeting LyC and Ly$\alpha$, show that those galaxies with the largest O32 are typically found to also be Ly$\alpha$ emitters. High O32 galaxies are also found to have tentative non-zero LyC escape fractions ($f_{esc}(LyC)$) based on $u$ band photometric detections. These results are consistent with samples of highly ionized galaxies, including confirmed LyC emitting galaxies from the literature. We also perform a detailed comparison between the observed emission line ratios and simulated line ratios from density bounded H$_{\textrm{II}}$ regions modeled using the photoionization code MAPPINGS V. Estimates of $f_{esc}(LyC)$ for our sample fall in the range from 0.0-0.23 and suggest possible tension with published correlations between O32 and $f_{esc}(LyC)$, adding weight to dichotomy of arguments in the literature. We highlight the possible effects of clumpy geometry and mergers that may account for such tension.