{'fraction_non_alphanumeric': 0.0302025262860772, 'fraction_numerical': 0.020958295109731143, 'mean_word_length': 4.5359375, 'pattern_counts': {'":': 0, '<': 12, '<?xml version=': 0, '>': 0, 'https://': 2, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 3, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The Third Data Release of the KODIAQ Survey

We present and make publicly available the third data release (DR3) of the Keck Observatory Database of Ionized Absorption toward Quasars (KODIAQ) survey. KODIAQ DR3 consists of a fully-reduced sample of 727 quasars at 0.1<z<6.4 observed with ESI at moderate resolution (4000<R<10000). DR3 contains 872 spectra available in flux calibrated form, representing a sum total exposure time of approximately 2.8 megaseconds. These co-added spectra arise from a total of 2753 individual exposures of quasars taken from the Keck Observatory Archive (KOA) in raw form and uniformly processed using a data reduction package made available through the XIDL distribution. DR3 is publicly available to the community, housed as a higher level science product at the KOA and in the igmspec database.


INTRODUCTION
Studies of absorption lines towards quasars have impacted a wide range of fields within astrophysics and cosmology. From studies of individual systems to determine the primordial D/H ratio (e.g., Cooke et al. 2018), to studies of tens of high resolution, high signal-to-noise spectra to constrain the temperature-density relation in the IGM (Hiss et al. 2019), to a determination of the mean free path of ionizing radiation from a stack of thousands of lines of sight (Prochaska & Wolfe 2009), quasar lines of sight have become one of the most valuable tools available to understand the universe at cosmological distances.
Large quasar and galaxy surveys such as SDSS (Ross et al. 2020) provide spectra in bulk in a ready-foranalysis state. By contrast, quasar spectra from larger ground-based facilities are often not available to the community in a science-ready format, if at all, and with little information on the provenance of the data. To help remedy the situation, a number of groups have provided the community with pipeline-reduced data from large aperture telescopes (e.g., Murphy 2018). As a part of these endeavours, a primary goal of the Keck Observatory Database of Ionized Absorption toward Quasars (KODIAQ, Lehner et al. 2014) is to uniformly reduce and make publicly available spectra of quasars obtained with instruments on the Keck telescopes over the last two-plus decades. In the first (O'Meara et al. 2015) and second (O'Meara et al. 2017) data releases (DR1 and DR2, respectively), KODIAQ has provided spectra from Keck's high-resolution echelle spectrograph HIRES. The DR2 sample encompasses DR1 and adds additional spectra, bringing the sample to a total of 300 unique quasar sightlines. In this work, we significantly increase the number of quasars with the third data release (DR3) from KODIAQ, this time adding 727quasar spectra obtained with the Echellette Sepctrograph and Imager (ESI, Sheinis et al. 2002) instrument on the Keck II telescope. Although lower resolution than HIRES, ESI spectra are of significant resolution to both serve statistical studies of the Lyman-α forest as well as metal line absorption, along with studies of individual absorbers such as damped Lyman Alpha (DLA) systems (e.g. Rafelski et al. 2012).

THE DATA
The ESI data in DR3 were obtained by individual PIs between December, 1999 (shortly after instrument commissioning) and January, 2015. ESI is an echellette system, with ten echelle orders spanning the range 3930 < λ < 10930Å. Cosmetic defects in the ESI CCD impact some pixels in the wavelength regions 4460 < λ < 4565Å. Depending on seeing conditions, slit width, and desired wavelength, ESI has a resolution of 4000 < R < 10000. Observations in DR3 are made using the full range of ESI slit widths, ranging between 0.3 and 6.0 arcseconds, with the majority of observations made with the 0.5, 0.75, and 1.0 arcsecond slits. The observations are summarized in Table 1. The coordinates and quasar emission redshifts presented in Table 1 were obtained by passing the coordinates from the raw data file first through the Simbad 1 databases and then, if not found, the NED 2 database, selecting the quasar closest on the sky. If an object was not found within thirty arcseconds of a quasar in either database, the coordinates were taken directly form the data header. We caution that in the case of close quasar pairs, the coordinates, chosen to correspond to the nearest object, may be inaccurate and refer to the other pair member.
2.1. Data Reduction 1 http://simbad.u-strasbg.fr/simbad/ 2 https://ned.ipac.caltech.edu Once the data were downloaded and reviewed for calibrations, they were uniformly reduced using the ESIRedux 3 code. ESIRedux is a suite of IDL routines within the XIDL 4 package of astronomical utilities. ESIRedux is functionally very similar to HIRedux, which is described in O'Meara et al. (2015), but with a few important differences which we outline below. The workflow of ESIredux is as follows: First, the raw data is grouped into setups. Unlike for HIRES, ESI is essentially fixed in wavelength coverage and echelle order placement on the detector, so setups are governed only by slit width and detector binning. Next, for each setup, a bias level and gain is determined, and a flat field image is constructed from a combination of internal lamps, dome flats, and if present, twilight sky flats. The flat field images serve the traditional role of removing pixel-to-pixel variations, but also provide the location of the slit edges for each echelle order. Twilight flats are used primarily in the bluest orders, where the internal and dome flats have lower flux. A two-dimensional wavelength image for the detector is then created using a combination of CuAr, HgNe, and Xe lamps. The wavelength solutions are compared to archived databases, and the residuals to the fit are most commonly found to be sub-pixel. Cosmic ray rejection is then performed using standard frame comparison techniques. As with HIRedux, sky subtraction and object definition and extraction is performed using the methods outlined in Bochanski et al. (2009) on an order-by-order basis, resulting in a spectrum of each order with counts as a function of wavelength. Finally, for objects with multiple exposures, the individual orders are weighted mean combined.
A key difference between the ESI data in KODIAQ DR3 and the HIRES data from DR2 is that the data in DR3 have been flux calibrated instead of continuum normalized. For expediency, all exposures in DR3 were flux calibrated using a single archival sensitivity function derived from a median of sensitivity functions from multiple flux standard stars. Examples of DR3 ESI QSO spectra are shown in Fig. 1.

Issues of note
A number of issues should be noted before using data from DR3 for precision analysis. First, the reddest echelle order, with a central wavelength of approximately λ = 9800Å was often poorly extracted. In these cases, that order has been rejected, and the data trimmed. The rightmost column in Table 1 indicates the wavelength range of the spectra in DR3. Second, the detector defects at λ ≈ 4500Å are present in the reduced data, and caution should be used near those wavelengths. Finally, some spectra suffer discontinuities at the regions in wavelength where the echelle orders overlap. These discontinuities are the result of orderto-order differences in flux calibration. Users that intend an analysis of the intrinsic quasar spectra should first inspect the data for these occurrences. We note that all intermediate data reduction steps for the data in DR3 are provided in the Keck Observatory Archive (KOA), 5 so that users may adjust specific steps, such as flux calibration, to match their needs. Finally, Table 1 shows that a number of quasars were observed on multiple dates, and often by multiple PIs. We did not attempt to combine the data across these observations, given frequent differences in slit width and potential inaccuracies in flux calibration. 5 https://koa.ipac.caltech.edu/Datasets/KODIAQ/ 2.2. KODIAQ at the KOA and igmspec As with DR2, the data products of DR3 are available for community download from the KOA. As DR3 stems from a different instrument, KOA hosts DR3 as a separate contributed data set from DR2. Nevertheless, a number of quasars have data in both DR3 and DR2. Table 2 lists those quasars. Note that in a few cases, the HIRES DR2 name does not exactly match the ESI DR3 name. Those cases are highlighted in Table 2. As with DR2, within DR3 at KOA, users can search for and download individual quasar data, or the full DR3 data set at once. DR3 will make available all intermediate data reduction products available, grouped by observing run. Spectra from DR3 will also be ingested into the igmspec database 6 .

PROPERTIES OF KODIAQ DR3
The KODIAQ DR3 consists of ESI spectra for 727 quasars. As in previous data releases, quasars are named according to their J2000 R.A./decl. coordinates as resolved by SIMBAD or NED, unless the quasar is not listed in these databases, in which case we use the position on the sky as provided by the raw data headers. In Table 1 we list the quasar redshift as provided by SIMBAD or NED. If the quasar is unlisted, a crude estimate is made for the redshift by using the Lyman−α and other emission lines. The full DR3 is comprised of 2753 individual exposures, grouped into 872 spectral coadds. The aggregate exposure time of the DR3 is ∼ 2.8 megaseconds.

General Properties
Fig. 2 shows the quasar redshift distribution for the DR3. The median redshift for the sample isz = 3.21 and the sample spans the range in redshift of 0.10 < z < 6.44. The distribution of DR3 quasars on the sky is shown in Figure 3.  In Fig. 4, we show median signal to noise per pixel, derived from the pipeline optimal extraction values of the flux and 1-σ errors, within ±5Å for the 872 in DR3 at 3 rest wavelengths λ r where λ r = λ obs /(1 + z em ) and λ obs is the observed wavelength in the spectrum.
The three values of λ r are chosen as markers for typical wavelengths of interest for IGM/CGM studies, namely λ r = 1450Å for studies of heavy element absorption free of Lyman-α forest contamination, λ r = 1170Å for stud-ies of the forest itself, and λ r = 1450Å for determination of the H I column density (N HI ) in Lyman limit systems from the shape of the Lyman break.

Cosmological Properties
Given the large number of quasars in DR3, and with a significant fraction of the spectra having moderate to high signal to noise, the sample can lend itself well to a variety of statistical analysis in cosmology and elsewhere. Here, we highlight a number of statistical properties of DR3 with respect to its ability to be used for studied of neutral hydrogen and heavy element absorption. In Fig. 5 we show the rest wavelength coverage provided by the quasars in DR3 in the quasar rest frame.
Vertical lines denote the rest frame wavelengths of commonly studied lines. Fig. 5 illustrates that the number of quasars that can be used for statistical studies of H I and heavy element absorption numbers well into the hundreds. Fig. 6 shows this large sample in a different way, namely showing as a function of redshift how many quasars can be used for specific ions. We can further subdivide the DR3 sample by exploring how many quasars provide coverage for specific ions as a function of redshift and signal to noise. Fig. 7 provides this summary for four ions: H I Lyman-α, C IV, O VI, and Mg II. Finally, Fig. 8 provides the redshift sensitivity function g(z) for the DR3 sample. The function g(z) is calculated in the same manner as in O'Meara et al. (2015) and O'Meara et al. (2017). g(z) is presented in Fig. 8 for 4 values of the signal-to-noise ratio, and for H I Lyman-α and C IV as representative ions. In all the above metrics, DR3 represents a significant sample for statistical studies of hydrogen and heavy element absorption, even at moderate to high signal-to-noise, which is a requirement for studies of weak absorption.

SUMMARY AND FUTURE
We present here a summary of and make publicly available the 727 quasar sample in the KODIAQ DR3. As with KODIAQ DR1 and DR2, all spectra and in-termediate data reduction products are made available at the KOA. Unlike DR1 and DR2, the data are roughly uniform in wavelength coverage (and nearly always larger per quasar in wavelength coverage than the HIRES spectra in previous releases), but like DR1 and DR2, they sample a large range of quasar redshift, and provide a wide range of signal-to-noise.
When using data products from DR3, we ask that publications provide the standard KOA acknowledgement, acknowledgment of the original PIs that obtained the data, and the following text: "Some or all of the data presented in this work were obtained from the Keck Observatory Database of Ionized Absorbers toward QSOs (KODIAQ), which was funded through NASA ADAP grants NNX10AE84G and NNX16AF52G, along with a citation to this paper. ACKNOWLEDGEMENTS Support for this research was made by NASA through the Astrophysics Data Analysis Program (ADAP) grants NNX10AE84G and NNX16AF52G, along with NSF grant award number 1516777. This research has made use of the Keck Observatory Archive (KOA), which is operated by the W. M. Keck Observatory and the NASA Exoplanet Science Institute (NExScI), under contract with the National Aeronautics and Space Administration. The data presented herein were obtained at the W.M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support of the W.M. Keck Foundation.
The authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of Maunakea has always had within the indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain. The authors wish to recognize and sincerely appreciate the work of the entire WMKO staff over the last two-plus decades, and to the efforts of the team at the NASA Exoplanet Science Institute (NExScI) who are responsible for maintaining the KOA.          Table 1 continued     Table 1 continued