\section{Introduction}\label{intro}
The production of elements with $Z\ge38$ occurs via the rapid ($r$-) and slow ($s$-) neutron-capture processes. While the site of the $s$-process is
well-understood \citep{gallino2005,karakas2010,frebel-rev2018}, a number of sites for the $r$-process\ have been proposed over the last few decades. A key requirement is the ability to provide the strong neutron flux needed for rapid neutron capture to occur. This has been predicted to happen in extreme sites like neutron star binary mergers (NSBM) (e.g., \citealt{lattimer1974}), magneto-rotationally driven jets (e.g., \citealt{winteler2012}), or collapsar disk winds (e.g., \citealt{siegel2019}).

Metal-poor stars enhanced in $r$-process\ elements have been found in both the Milky Way halo (e.g., \citealt{rpa1}, \citealt{rpa2} and \citealt{roederer2018a}) as well as dwarf galaxies, most notably the $r$-process ultra-faint dwarf galaxy Reticulum\,II (Ret\,II) \citep{ji2016Nat,roederer2016}. These stars formed from gas that was previously enriched by a nucleosynthesis event during which heavy $r$-process elements were made. Altogether, they thus provide the best accessible evidence for the early operation of the $r$-process, as many $r$-process elements can be detected in their spectra, and the corresponding chemical abundances be measured.

In 2017, the LIGO/Virgo gravitational wave observatory discovered the first NSBM event GW170817 \citep{abbott2017}, after which a kilonova, likely powered by the decay of newly synthesized of $r$-process\ elements, was observed (e.g., \citealt{pian2017}). Together with chemical-evolution models of the Galaxy that implement NS-NS mergers \citep{cescutti2015,wehmeyer2015,cote2017} and the existence of Reticulum\,II, these lines of evidence support the hypothesis that NS mergers are likely a dominant site of $r$-process\ nucleosynthesis. However, improved chemical-evolution models by \cite{cote2018b}, more recently suggested that NSBM may not in fact be the only sources of $r$-process\ elements in the Galaxy. Additional sources, such as magneto-rotational (``jet'') supernovae \citep{nishimura2017} could be needed to reproduce the observed $r$-process\ patterns in, for e.g., moderately enhanced $r$-process\ ($r$-I: $+0.3<\mathrm{[Eu/Fe]}\leq +1.0$) metal-poor stars. Alternative sites for highly $r$-process-enhanced ($r$-II: [Eu/Fe] $> +1.0$) stars have also been suggested, such as collapsars \citep{siegel2019}.
This conundrum underscores the need for studying large numbers of $r$-process-enhanced stars, to better interpret and place constraints on heavy-element formation through the $r$-process, and all of its associated astrophysical site(s) throughout cosmic history.

Motivated by the importance of $r$-process-enhanced halo stars, the $R$-Process Alliance (RPA) collaboration aims at significantly increasing the numbers of such stars, as, until recently, only a relatively small number of them were known. Moreover, a major goal is to combine $r$-process\ observations, $r$-process\ theory efforts, and results from chemical-evolution simulations to eventually produce a more complete understanding of the origin of the $r$-process-enhanced star population in the Galactic halo.  To this end, two extensive ``data release'' papers have already been published by \citet{rpa1} (hereafter RPA-1) and \citet{rpa2} (hereafter RPA-2) which report on many newly discovered $r$-process-enhanced stars. This paper (hereafter RPA-3) describes the third data release, based on spectra collected with the 6.5\,m Magellan Clay telescope.

With the emergence of a large $r$-process-enhanced star sample, and thus statistically meaningful measures of chemical abundances, kinematics, frequencies, etc. of these stars, questions about their formation and any prior enrichment of their birth gas clouds can soon be addressed more rigorously. This is also of interest with respect to
establishing possible connections between the $r$-process\ signature observed in the Ret\,II stars and that of $r$-process\ stars found in the Galactic halo. Given their nearly identical abundance patterns, the halo stars likely originated in ancient dwarf galaxy analogs similar to Ret\,II \citep{roederer2018}. Establishing the origin story of halo $r$-process-enhanced stars will provide important insides into the environment in which the earliest $r$-process events took place.

 From a cosmological perspective, the $r$-process-enhanced stars also offer another advantage, since they can be dated from their abundance ratios of long-lived isotopes of Th and U relative to abundances of stable $r$-process\ elements such as Eu (e.g., \citealt{cayrel2001,frebel2007,placco2017,holmbeck2018}). With half-lives of 4.5 Gyrs and 14.05 Gyrs, respectively, $^{238}$U and $^{232}$Th decay over cosmologically relevant timescales.
Such ages could be used, for example, to place lower limits on the time of their associated $r$-process\ production event(s). We note, however, that such ratios can not be reliably used for actinide-boost $r$-process\ stars \citep{schatz2002}, which are known to exhibit unusually high Th and U abundances as compared to stable elements, such as Eu. Interestingly, the actinide elements ratio, Th/U (if measured simultaneously in a star), seems to remain a robust tool to determine the ages of actinide-boost $r$-process\ stars \citep{holmbeck2019}.

This paper is outlined as follows.  In Section\,\ref{sec:obs}, we describe the observations, data reduction, and radial-velocity measurements. In Section\,\ref{sec:stell_param}, we present the stellar parameter determinations of the sample using both
1D (1-Dimensional), LTE (Local Thermodynamic Equilibrium) and 1D, NLTE (Non-Local Thermodynamic Equilibrium) methods. Section\,\ref{sec:abund} presents the chemical abundances of the light elements and the $\alpha$-elements, as well as the Fe-peak and neutron-capture elements determined for the stars. We then  discuss the results and conclude in Sections\,\ref{sec:disc} and \ref{sec:conc}, respectively.

\section{Observations}\label{sec:obs}
\subsection{Target Selection and Observations}
Metal-poor, candidate $r$-process-enhanced target stars were selected from various metal-poor surveys, including the RAdial Velocity Experiment (RAVE; \citealp{RAVE2006}),  Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST) (C. Liu, private communication), Best \& Brightest (B\&B; \citealp{schlaufman2014,placco2019}), SkyMapper\footnote{The SkyMapper stars are taken from the extremely metal poor (EMP; \AB{Fe}{H} $<-3$) sample discussed in \citet{dacosta2019}; specifically it includes all stars that satisfy $g_{\mathrm{SkyMapper}}$ $\leq$ 14.1 and --3.0 $\leq$ [Fe/H]$_{\mathrm{fitter}}$ $\leq$ --2.0, where [Fe/H]$_{\mathrm{fitter}}$ is the abundance estimated from the low-resolution spectra (see \citealt{dacosta2019} for details).} \citep{skymapper2018}, Melendez \& Placco (M\&P; \citealp{melendez2016}), and the Hamburg/ESO Survey (HES; \citealp{H-ESO2008}). Most of our targets were selected from the RAVE fourth and fifth data releases (DR4 and DR5; \citealp{kordopatis2013,kunder2017}). For those latter stars, their metal-poor nature was vetted by \citet{placco2018}, who derived estimates of their stellar parameters based on medium-resolution spectroscopy acquired on a number of different telescopes. Additional details on the selection criteria for these stars is presented in \citet{rpa1} and \citet{placco2018,placco2019}.

All targets were observed at high spectral resolution using the Magellan Inamori Kryogenic Echelle (MIKE) spectrograph \citep{bernstein2003} on the Clay (Magellan\,II) telescope at Las Campanas Observatory, Chile. Overall, a total of 148 stars with visual magnitudes ($V$) between 10.0 and 14.5 were observed. The spectra were collected over the course of three years, from January 2015 through December 2018. Either the 0\farcs7 or 1\farcs0 slit with $2\times2$ binning was used, yielding nominal resolving powers of $R\sim$28,000 or $R\sim$22,000 in the red, and $R\sim35,000$ or $R\sim28,000$ in the blue, respectively. All spectra cover the wavelength range from $\sim3860$\,{\AA} to $\sim9000$\,{\AA}.

Exposure times depended on magnitude and weather conditions, but were typically between 900\,s for brighter stars to 1800\,s for fainter ones, aiming for Signal-to-Noise (S/N) values of 25\,pixel$^{-1}$ at 4000\,{\AA}. When needed, more than one exposure was taken to increase the S/N to the desired level. The quality of these spectra enables the detection or a meaningful upper limit on the \euii\ line at 4129\,{\AA}, which serves as the major $r$-process\ diagnostic. Many other elements can also be detected, including C, Sr, and Ba. The Two Micron All Sky Survey (2MASS) IDs, right ascension (R.A.) and declination (Decl.), visual magnitudes ($V$), exposure times, Modified Julian Dates (MJD), slit widths, S/N ratios at 4000\,{\AA} and survey sources are listed in Table\,\ref{tab:ident}.

\subsection{Data Reduction and Radial-Velocity Measurements}

Data were reduced using the latest versions of the  \texttt{Carnegie Python Distribution}\footnote{https://code.obs.carnegiescience.edu/mike} \citep{kelson2003}. Each order of each spectrum was then normalized and merged into a final spectrum, which was also
radial-velocity shifted by cross-correlation with the \ion{Ca}{ii} line near 8542\,{\AA}
of the spectrum of HD\,140283, using the \texttt{Spectroscopy Made Hard} (\texttt{SMH}) analysis software (first described in \citealt{casey2014}). Heliocentric velocity (RV$_{\mathrm{helio}}$) corrections were then determined with the \texttt{rvcorrect} package in \texttt{IRAF}. We compare our RV$_{\mathrm{helio}}$ values to those from the $Gaia$ second data release (DR2; \citealp{gaia2018}), whenever available (110 stars). We present the comparisons in Figure\,\ref{fig:rad_vel}.
The radial velocities for eleven of the stars differ significantly from the $Gaia$ velocities, up to 60\,km\,s$^{-1}$. Excluding these outliers, our velocities agree with those from {\it Gaia} within a mean difference of 1.93\,km\,s$^{-1}$, with a standard deviation of $\sigma=3.10$\,km\,s$^{-1}$. We also compute the biweight-scale standard deviation\footnote{https://docs.astropy.org/en/stable/api/astropy.\\{\hspace*{0.55cm}stats.biweight.biweight\_scale.html}} of the full sample (including outliers). We find that $\sigma_{\mathrm{biweight}}=3.73$ agrees well with our reported $\sigma$.
Seven of those are found to be double-lined spectroscopic binaries (SB2), as their spectra exhibit clear contamination from a companion; these are labeled as such in Table\,\ref{tab:ident}. No spectral disentanglement of the two components was attempted for these stars, and thus they were not further analyzed with the rest of the sample.  We note that it is highly unlikely that the source of the $r$-process\ enhancement is due to binarity, as it is unable to produce the neutron flux density ($n_{n}>10^{22}$\,cm$^{-3}$) needed to induce the $r$-process\ \citep{thielemann2017}.

\begin{figure}[ht!]
\begin{center}
\vspace*{-0.4cm}
\hspace*{-0.9cm}
\includegraphics[scale=0.2]{rvs_comp_helio_gaia.pdf}
\caption{\label{fig:rad_vel}Histogram showing radial-velocity differences of our targets, as measured in this study and by Gaia DR2. Eleven stars display significant differences $>1\sigma = 3.09$\,km\,s$^{-1}$\ (value determined excluding outliers). Seven of those stars were found to be SB2 binaries, and are labeled as such in Table\,\ref{tab:ident}.}
\end{center}
\end{figure}

\section{Stellar Atmospheric Parameters}\label{sec:stell_param}

We determined stellar atmospheric parameters for the sample stars, namely the effective temperature ($T_{\rm eff}$), surface gravity ($\log g$), metallicity (\AB{Fe}{H}), and microturbulent velocity ($\xi_{t}$), from the abundances of \fei\ and \feii\ lines as determined from equivalent width (EW) measurements. The EWs of the lines were obtained by fitting Gaussian line profiles to the spectral absorption features using \texttt{SMH}. The \fei\ and \feii\ abundances were derived with 1D stellar atmospheric models, using three different spectroscopic methods: (i) assuming LTE, (ii) assuming LTE with $T_{\rm eff}$\ corrected to a photometric scale, following \citet{frebel2013}, as well as (iii) assuming NLTE. The atomic line-by-line properties for the \fei\ and \feii\ lines, including their EW and abundance measurements, are shown in Table\,\ref{tab:linelist}. The linelist is compiled from several data sources, and the $\log\,gf$ values for our lines are adopted from the up-to-date list compiled by \citet{roederer2018a} (see their Table\,2 and the references therein).

\subsection{LTE Stellar Parameters}\label{sec:lte}
The LTE stellar atmospheric parameters were estimated from the abundances of \fei\ and \feii\ lines using the 2017 version of the LTE radiative transfer code \texttt{MOOG} \citep{sneden1973}, including Rayleigh scattering treatment following \citet{sobeck2011}\footnote{https://github.com/alexji/moog17scat}.  Stellar 1D, LTE atmospheric models from \citet{castelli2004}, with standard $\alpha$-element enhancement of [$\alpha$/Fe] = +0.4 were used.
$T_{\rm eff}$\ values were determined by enforcing excitation equilibrium of the abundance of \fei\ lines as a function of excitation potential, $\chi$. The $\log g$\ values were determined by enforcing ionization equilibrium between the abundances inferred from the \fei\ and \feii\ lines. The $\xi_{t}$\ parameter was determined by ensuring no \fei\ abundance trends with reduced equivalent widths existed (REW; $\log(\mathrm{EW}/\lambda)$). The \AB{Fe}{H}\ abundances were determined from the average abundances of the \fei\ and \feii\ lines.
The number of \fei\ and \feii\ lines used, as well as the LTE stellar parameters derived for each star, are listed in Table\,\ref{tab:stell_param}.
Stellar  parameter  uncertainties  are  determined assuming systematic uncertainties following the analysis in \citet{ji2016b}, and they are estimated to be up to 150\,K, 0.3\,dex, and 0.2\,km\,s$^{-1}$\ for $T_{\rm eff}$, $\log g$, and $\xi_{t}$, respectively. The \AB{Fe}{H}\ uncertainties are determined from the standard deviations of the \fei\ and \feii\ abundances.

\subsection{Corrected LTE Stellar Parameter Estimates}\label{sec:lte_corr}
Spectroscopic stellar parameters derived using the assumption of LTE have long been known to suffer from systematic uncertainties, especially for cooler metal-poor stars where atomic collisions are not sufficiently frequent to ensure that actual excitation and ionization equilibrium of the atomic populations is achieved.
Determinations of $T_{\rm eff}$\ via photometric methods have been shown to be less prone to uncertainties compared to the LTE spectroscopic approach \citep{casagrande2010}. However, these require known absolute magnitudes and reddening toward the stars.
As a middle ground, we thus also determined the $T_{\rm eff}$\ of our target stars by ``correcting'' the LTE spectroscopic temperatures to a photometric scale, following the empirical calibration introduced in \citet{frebel2013}, hereafter denoted as $T_{\rm eff}$(FR13 corr.):

\begin{displaymath}
    {T_{\mathrm{eff}}\text{(FR13 corr.)} = {0.9\times} {T_{\mathrm{eff}}\text{(LTE)}} + 670}
\end{displaymath}

For comparison to the ``strict'' LTE $T_{\rm eff}$\ results, we divide the stars into three categories covering different $T_{\rm eff}$\ ranges: $T_{\rm eff}$\ $\leq4500$\,K, $4500 < $ $T_{\rm eff}$\ $< 5000$\,K and $T_{\rm eff}$\ $\geq5000$\,K. We show the differences obtained relative to the strict LTE results in Figure \ref{fig:nlte_corr}, for each $T_{\rm eff}$\ range.
By design, the differences in the temperatures are more pronounced for cooler stars ($T_{\rm eff}$\ $<4500$\,K), by up to 250\,K, while smaller shifts were obtained for the hotter stars ($T_{\rm eff}$$>5000$\,K), up to 150\,K.

After fixing $T_{\rm eff}$\ to the FR13 corrected scale, we then also derived $\log g$, $\xi_{t}$, and \AB{Fe}{H}\ as explained in Section~\ref{sec:lte}, by forcing both ionization equilibrium and no abundance trends with line strengths. We compare the ``corrected'' $\log g$, $\xi_{t}$, and \AB{Fe}{H}\ to those determined strictly in LTE, as again shown in the boxplots in Figure\,\ref{fig:nlte_corr}. For $\log g$, the differences are presented for three surface gravity sub-categories: $\log g$\ $\leq1.5$, $1.5 < $ $\log g$\ $ < 3.0$ and $\log g$\ $\geq3.0$. The shifts from LTE are largest for the giant stars with $\log g$\ $\leq1.5$, up to $\sim +1.0$\,dex, while slightly smaller differences, up to +0.7\,dex, were obtained for the giants and sub-giants with $1.5 <$ $\log g$\ $< 3.0$. For stars with $\log g$\ $\geq 3.0$, the resulting differences were smaller still, up to +0.2\,dex, and lie within the LTE $\log g$\ uncertainties.

We also note that the differences obtained for \AB{Fe}{H}\ increase with decreasing metallicities, and are the largest for the extremely metal-poor stars (EMP) with \AB{Fe}{H}\ $\leq-3.0$ (Figure\,\ref{fig:nlte_corr}), reaching up to +0.5\,dex. On average, the corrections are $\sim +0.2$\,dex.  The corrections for $\xi_{t}$\ are negative relative to the uncorrected LTE case, and are decreasing with increasing $\xi_{t}$, reaching $-0.5$\,km\,s$^{-1}$\ for $\xi_{t}$\  $\geq3.0$\,km\,s$^{-1}$.

\begin{figure*}[ht!]
\begin{center}
\hspace*{-1.3cm}
\includegraphics[scale=0.2]{rpa_nlte_corr.pdf}
\caption{\label{fig:nlte_corr}Boxplots showing the stellar parameter changes for $T_{\rm eff}$, $\log g$, \AB{Fe}{H}, and $\xi_{t}$, when comparing LTE results with those from applying the \citet{frebel2013} correction (labeled as FR13 Corr.), and from using a NLTE approach (labeled NLTE).
The boxplots are color-coded for the ranges of stellar parameters (in LTE), as indicated on top of each panel. For $\log g$, we also show the differences arising from using $Gaia$ DR2 parallaxes and our uncorrected LTE values. Red uncertainty bars represent typical NLTE parameter uncertainties, corresponding to 100\,K for $T_{\rm eff}$, 0.3\,dex for $\log g$, 0.2\,dex for \AB{Fe}{H}, and 0.2\,km\,s$^{-1}$\ for $\xi_{t}$.}
\end{center}
\end{figure*}

\begin{figure*}
\begin{center}
\vspace*{-0.4cm}
\hspace*{-1cm}
\includegraphics[scale=0.2]{isoch_nlte.pdf}
\caption{\label{fig:isoch}Left panel: NLTE stellar parameters of the target stars ($T_{\rm eff}$\ and $\log g$\ color-coded as a function of \AB{Fe}{H}), over-plotted with Dartmouth isochrones at 12.5\,Gyrs for four different metallicities, $\AB{Fe}{H}\ =-1.00$, $\AB{Fe}{H}\ =-1.50$, $\AB{Fe}{H}\ =-2.00$, and $\AB{Fe}{H}\ =-2.50$, respectively \citep{dartmouth2008}. Right panel: NLTE $\log g$\ versus $\xi_{t}$\ for the target stars, color-coded as a function of \AB{Fe}{H}. The polynomial fit to the data is shown by the gray dotted line. Also shown, for comparison, are the stars (show by the small star markers) from \citet{cayrel2004} and \citet{barklem2005}, as well as the polynomial fit to their $\log g$-$\xi_{t}$\ (black dotted line).}
\end{center}
\end{figure*}

\subsection{NLTE Stellar Parameter Estimates}\label{sec:nlte}

Atmospheric stellar parameter determinations for metal-poor stars from LTE spectroscopic methods are affected by unaccounted-for departures from statistical equilibrium that can introduce significant systematic uncertainties, since line formation and populations of non-dominant species (in this case \fei) can potentially deviate from the Saha–Boltzmann equilibrium assumed in LTE \citep{lind2012,amarsi2016,Ezzeddine2017}. To account for such departures, especially in giant and warm stars, it is necessary to investigate the formation of iron lines  in NLTE.
Therefore, we also determine stellar parameters for our target stars using 1D, NLTE radiative transfer models.

\subsubsection{NLTE Methods}

NLTE abundances were computed for \fei\ and \feii\ lines from their EWs using the radiative transfer code \texttt{MULTI} in its 2.3 version \citep{Carlsson1986,carlsson1992}, and employing 1D \texttt{MARCS} model atmospheres \citep{gustafsson1975,gustafsson2008} interpolated to the corresponding parameters. Blanketing from background opacities, excluding Fe lines, were employed from the \texttt{MARCS} opacity package (B. Plez, private communication).

The \fei/\feii\ atomic model used in the NLTE calculations is described in more detail in \citet{Ezzeddine2017}. This model was built adopting up-to-date atomic data. Quantum computations of inelastic collision rates of \fei\ with neutral hydrogen atoms were implemented from \citet{barklem2018}. These collisions play an important dominant role (over electrons) in NLTE calculations of cool stars.

The NLTE stellar parameters were derived following the procedure in \citet{Ezzeddine2017}, where first approximations of the parameters were derived simultaneously by multi-dimensional fitting of the observed \fei\ and \feii\ EW of the target stars to theoretical NLTE EW in a pre-computed dense grid of $T_{\rm eff}$, $\log g$, $\xi_{t}$, and \AB{Fe}{H}\ parameter space, using the Levenberg-Marquardt algorithm. The excitation and ionization equilibrium, as well as the trends with reduced EW, were then checked for each star, ensuring no trends.

\subsubsection{NLTE Results}
The derived NLTE stellar parameters are shown in Table\,\ref{tab:stell_param}. For $T_{\rm eff}$, $\log g$, and $\xi_{t}$, we adopt typical total uncertainties arising from the uncertainties in the atomic data and abundance slopes versus $\chi$ and REW, of 100\,K in $T_{\rm eff}$, 0.3\,dex in $\log g$, and 0.2\,km\,s$^{-1}$\ in $\xi_{t}$. \AB{Fe}{H}\ uncertainties were adopted from the standard deviations of the line-by-line abundances and are shown in  Table\,\ref{tab:stell_param}, which vary from 0.08 to 0.18\,dex.
The NLTE gravities, as a function of $T_{\rm eff}$\ and \AB{Fe}{H}, are shown on the left panel of Figure~\ref{fig:isoch}, over-plotted on four 12.5\,Gyr Dartmouth isochrones \citep{dartmouth2008}, with [$\alpha$/Fe] = +0.4 and \AB{Fe}{H}\ $=-1.00$, \AB{Fe}{H}\ $=-1.50$, \AB{Fe}{H}\ =$-2.00$, and \AB{Fe}{H}\ $=-2.50$, respectively.
 The stars mostly occupy the upper portion of the giant branch in the $T_{\rm eff}$-$\log g$\ diagram. This likely stems from our target selection bias toward cool bright giants for the search for $r$-process stars by the RPA \citep{placco2018,placco2019}. We note that, despite this bias, the RPA sample in this work, as well as other RPA data releases, are expected to be fairly representative, as all initial “snapshot” spectra of the target stars were exposed to determine a reliable r-process enhancement measurement, by adjusting exposure times accordingly to obtain
S/N ratios $> 30$ (per pixel) in the blue around the \euii\ line at 4129\,{\AA} (a representative line of $r$-process\ enhancement, see Section\,\ref{sec:rproc_abund} below).
Four of the stars are found to be hot dwarfs ($T_{\rm eff}$\ $>5500$\,K, $\log g$\ $>4.0$), while another two stars are field horizontal-branch stars.  Additionally, one of the stars, 2MASS~J09291557+0838002, is a newly identified ultra-metal-poor (UMP) star at \AB{Fe}{H}(NLTE)$=-4.02$. Four $r$-process-enhanced stars at the $r$-I, $r/s$, and limited-$r$\ levels are detected at \AB{Fe}{H}\ $>-1.5$. We also plot, on the right panel of Figure\,\ref{fig:isoch}, $\log g$\ versus $\xi_{t}$\ for our target stars. For comparison, we also show the results from previous studies by \citet{cayrel2004} and \citet{barklem2005} for their very metal-poor star samples. We fit the $\log g$-$\xi_{t}$\ dependence for both our target stars, as well as those from previous studies with third degree polynomials\footnote{The fitting coefficients for the third degree polynomial,  y = ax$^3$ + bx$^2$ + cx + d (where x and y correspond to $\log g$\ and $\xi_{t}$, respectively), obtained for our target stars are a = $-0.02$, b = 0.37, c= $-1.98$, and d = 4.73.}.
We find that $\xi_{t}$\ increases for giant stars with lower $\log g$\ as compared to the dwarfs. A similar trend is obtained from previous studies, as shown by the black star markers and their third degree polynomial fit on the same Figure, although the NLTE dependence from our study shows a more severe increase for $\xi_{t}$\ for the giant stars. These positive NLTE $\xi_{t}$\ corrections are nominally driven by the abundance differences obtained between LTE and NLTE from the stronger saturated \fei\ lines (with $EW>50$\,m{\AA}), as explained in the detailed NLTE study by \citet{lind2012}.

We also compare the NLTE stellar parameters from our sample and from the RPA-1 and RPA-2 papers \citep{rpa1,rpa2} to their corresponding LTE values (Figure\,\ref{fig:nlte_corr}).
We find that the NLTE temperatures are higher than those determined in LTE, and the shifts are more pronounced for cooler stars, reaching up to +200\,K. On average, the $T_{\rm eff}$\ shifts between NLTE and LTE are +150\,K for $T_{\rm eff}$\ $\leq4500$\,K, +120\,K for $4500 <$ $T_{\rm eff}$\ $\leq5000$\,K and +50\,K for $T_{\rm eff}$\ $>5000$\,K. The NLTE corrections ($T_{\rm eff}$(NLTE)$-$ $T_{\rm eff}$(LTE)) are, however, lower than those obtained by using the FR13 corrections, by $\sim 100$\,K, on average.

The NLTE corrections in $\log g$, defined by $\log g$(NLTE) $-$ $\log g$(LTE), are also positive, and found to be increasing toward lower gravities, up to +1.2\,dex. These corrections agree with those derived for the  stars in \citet{rpa2}. On average, the corrections are +1.2\,dex for $\log g$\ $\,\leq1.5$, +1.0\,dex for $1.5 <$ $\log g$\,$ < 3.0$, and +0.5\,dex for $\log g$\ $\geq\,$ 3.0. The corrections in \AB{Fe}{H}\ are also increasing toward lower metallicities, in line with other NLTE predictions \citep{mashonkina2011,lind2012,amarsi2016}, and can reach up to +0.6 for \AB{Fe}{H}\ $\leq-3.0$, +0.45 for $-3.0 < $ \AB{Fe}{H}\ $ <-2.0$, and +0.3 for \AB{Fe}{H}\ $\geq-2.0$.
We note that these metallicity NLTE corrections are expected to be larger, on average, than those nominally obtained by fixing the $T_{\rm eff}$\ and $\log g$\ values prior to determining the LTE and NLTE \AB{Fe}{H}\ (e.g., \citealt{lind2012,mashonkina2011}). Therefore, the \AB{Fe}{H}\ corrections obtained in our study also reflect the abundance changes due to the $T_{\rm eff}$\ and $\log g$\ NLTE corrections obtained for the target stars.
Finally, the NLTE $\xi_{t}$\ corrections are negative, and agree on average with those obtained by applying the FR13 corrections. The corrections are, however, slightly lower in NLTE for $\xi_{t}$\ $\geq$ 3.0\,km\,s$^{-1}$, and can reach $-0.6$\,km\,s$^{-1}$.

We note that any differences in the atmospheric models used in the NLTE and LTE calculations (\texttt{MARCS} vs. \texttt{ATLAS}) are expected to have negligible effects on the resulting stellar parameters in LTE. \citet{roederer2014b} showed in their analysis of a large number of metal-poor stars, including both giants and subgiants, that differences obtained between both models are less than 8\,K for $T_{\rm eff}$, 0.005\,dex for $\log g$, 0.002\,km\,s$^{-1}$\ for $\xi_{t}$, and 0.006\,dex for \AB{Fe}{H}. These values are negligible as compared to the differences obtained for the stellar parameters between NLTE and LTE.

\subsection{Estimates of $log\,g$ from $Gaia$ DR2 Parallaxes}\label{sec:gaia}
For $\log g$, we also derive surface gravities from $Gaia$ DR2 parallax-based distances \citep{bailer-jones2018}. To calculate distance-based surface gravities, we adopt a relation between $\log g$,
stellar mass, $V_{\mathrm{mag}}$, and bolometric corrections.
For the dwarfs and giants in our sample, a mass of 0.7 and 0.8 solar masses was adopted, respectively. The RAVE $V_{\mathrm{mag}}$ were used, and the bolometric corrections were adopted from \citet{casagrande2014}.

The derived $\log g$$^{Gaia}$ values are shown in the last column of Table~\ref{tab:stell_param}. For comparison, we also show the differences  between $\log g$$^{Gaia}$ and those derived in LTE in Figure~\ref{fig:nlte_corr}. The gravity corrections are more pronounced for the giant stars with $\log g$\ $\leq1.5$, with an average shift of $\sim +1.2$\,dex. Slightly lower shifts are obtained for stars with $1.5 < $ $\log g$\,$ < 3.0$, of up to +0.6\,dex, and for stars with $\log g$\ $\geq\,$ 3.0, up to +0.1\,dex. The $Gaia$ gravities agree with those obtained in NLTE within uncertainties of $\sim$0.3\,dex, on average.

We note that 26 of our sample stars have  {\it Gaia} DR2 parallax-based distances larger than 6\,kpc \citep{bailer-jones2018}. Those stars are likely to lead to larger uncertainties in the derived $\log g$$^{Gaia}$ ($>0.2$\,dex), due to the larger systematic shifts of their {\it Gaia} DR2 parallaxes \citep{arenou2018}. We flag these stars as such in Table\,\ref{tab:stell_param}, and warn against the use of their {\it Gaia} $\log g$.

Given that the NLTE $\log g$\ values agree overall with those from non-spectroscopic ones ($Gaia$), and that the NLTE stellar parameters are considered to be better representatives of the stellar atmospheres of metal-poor stars than LTE (e.g., see review by \citealt{asplund2005}), we decided to adopt the NLTE stellar parameters for our target stars, and use them to determine the chemical abundances of our stars.
Nevertheless, our analysis clearly shows that in the absence of a full NLTE approach, the FRE13 values (calibrated to photometric $T_{\rm eff}$) provide reasonable stellar parameters, lying within the uncertainty levels of the corresponding NLTE values, where typical values of 100\,K for $T_{\rm eff}$, 0.3\,dex for $\log g$, 0.2\,dex for \AB{Fe}{H}, and 0.2\,km\,s$^{-1}$\ for $\xi_{t}$\ are shown by the red error bars on Figure\,\ref{fig:nlte_corr}.

\section{Chemical Abundances}\label{sec:abund}

We derive abundances for the light elements, $\alpha$ elements, and Fe-peak elements for the target stars, including C, O, Na, Mg, Al, Si, K, Ca, Sc, \tii, \tiii, V, Cr, Mn, Co, Ni, and Zn using \texttt{MOOG}. Additionally, we derive the abundances for the neutron-capture elements Eu, Ba, and Sr, in order to sort our target stars into the sub-classifications suggested by the reviews of \citet{Beers2005} and \citet{frebel-rev2018}. The abundance ratios [X/H] are calculated relative to the solar abundances presented by \citet{asplund2009}. The abundances (except for C, O, Al, and the neutron-capture elements) were derived from a line-by-line analysis using EWs of the lines and a curve-of-growth (COG) method. Sufficiently weak lines lying on the linear part of the COG were chosen with REW $\leq -4.5$. The remaining abundances were derived using spectrum synthesis of their corresponding lines, taking into account isotopic ratios.  Linelists for each region of interest are generated with \texttt{linemake}\footnote{https://github.com/vmplacco/linemake} \citep{sneden2008}.  The Hyper-fine structure (HFS) was considered for the Fe-peak elements, including Sc, V, Mn, and Co, when necessary.
For carbon, an average isotopic shift ratio of  $^{12}$C/$^{13}$C = 4/1 was used. This was determined for each star individually, by fitting the lines around the 4313\,{\AA} CH feature, specifically, those at 4217\,{\AA}, 4225\,{\AA}, as well as 4302\,{\AA}. Upper limits were determined by matching the noise levels in the observed spectral lines with the corresponding synthetic spectral lines. Our line-by-line atomic data, EW and abundances are listed in Table\,\ref{tab:linelist}. The average abundances are also shown in Table\,\ref{tab:abund_fe_peak}.

\subsection{Neutron-Capture Element Abundances}\label{sec:rproc_abund}
Strontium abundances were derived from the \srii\ lines at 4077\,{\AA}, 4161\,{\AA}, and 4215\,{\AA}. With the line at 4077\,{\AA} being saturated in most of the stars, and the line 4161\,{\AA} being very weak, especially in EMP stars, the Sr abundances were primarily derived from the line at 4215\,{\AA}.  We note that the strong 4077 and 4215\,{\AA}
\srii\ lines may be affected by NLTE effects, however, \citet{short2006} reported an offset of $<0.1$\,dex for these lines for stellar parameters similar to those of our sample of stars. Barium abundances were derived mainly from the two \baii\ lines at 5853\,{\AA} and 6141\,{\AA}, which are not saturated, unlike the line at 4554\,{\AA}. For Eu, the abundances were derived from several \euii\ lines at 3907\,{\AA}, 4129\,{\AA}, 4205\,{\AA}, 4435\,{\AA}, 6437\,{\AA}, and 6645\,{\AA}, following \citet{rpa1}. The bluer line at 3907\,{\AA} was used when the blue S/N of the spectra was high enough ($\geq80$). However, the strongest Eu line, located at 4129\,{\AA}, was was the primary line used for most stars. Examples of observed Sr, Ba, and Eu lines at 4215\,{\AA}, 5853\,{\AA}, and 4129\,{\AA}, respectively, for representative $r$-II, $r$-I, and limited-$r$\ stars, as well as their abundance syntheses, are shown in Figure\,\ref{fig:obs_spec}.

The abundance ratios [Sr/Fe], [Ba/Fe] and [Eu/Fe], calculated with NLTE \AB{Fe}{H}\ values, are listed in Table\,\ref{tab:rproc_clas}. They
were also used to classify our stars into their sub-classes following \citet{frebel-rev2018}:

\begin{itemize}
    \item $r$-II:  [Eu/Fe] $ > +1.0$ and [Ba/Eu] $ < 0.0$
    \item $r$-I: $+0.3 \leq $ [Eu/Fe] $ \leq +1.0$ and [Ba/Eu] $ < 0.0$
    \item limited-$r$: [Eu/Fe] $ < +0.3$, [Sr/Ba] $ > +0.5$, and [Sr/Eu] $ > 0.0$
    \item $s$-process: [Ba/Fe] $ > +1.0$ and [Ba/Eu] $ > +0.5$
    \item $r/s$: $0.0 < $ [Ba/Eu] $ < +0.5$
    \item CEMP-$r$: [C/Fe] $> +0.7$, [Eu/Fe] $ > +0.3$, and [Ba/Eu]$ < 0.0$
    \item CEMP-no: [C/Fe] $ >+0.7$ and [Ba/Fe] $ < 0.0$
    \item Non-$r$-process\ enhanced (non-RPE):  [Eu/Fe] $ < +0.3$, [Sr/Ba] $ \leq 0.5$, and [Sr/Eu] $ < 0.0$
\end{itemize}

\begin{figure*}[ht!]
\begin{center}
\hspace*{-1.0cm}
\includegraphics[scale=0.2]{rproc_syn.pdf}
\caption{\label{fig:obs_spec}Observed spectral lines (dotted lines) and their abundance syntheses (gray lines) of \srii\ at 4215\,{\AA}, \baii\ at 5853\,{\AA},  and \euii\ at 4129\,{\AA} for $r$-II\ (J00101758$-$1735387, top panel), $r$-I\ (J12233047$-$4216365, middle panel), and limited-$r$\ (J18405985$-$4841353, bottom panel) stars discovered in our sample. Representative uncertainties of $\pm 0.2$\,dex are shown for each lines by the filled syntheses colors.}
\end{center}
\end{figure*}

\begin{figure*}[ht!]
\begin{center}
\hspace*{-1.3cm}
\includegraphics[scale=0.42]{rproc_hist.pdf}
\caption{\label{fig:hist}Cumulative histograms showing the number of $r$-I\ (filled gray), $r$-II\ (filled red), and limited-$r$\ (filled blue) stars discovered in this work, as a function of NLTE \AB{Fe}{H} (upper panels) and V magnitudes (lower panels). Also shown, for comparison, are all the stars discovered and analyzed in the RPA-1 \citep{rpa1} and RPA-2 \citep{rpa2} samples, as well as the total targets (dashed gray histogram) analyzed in this work (including non-RPE).}
\end{center}
\end{figure*}

\begin{figure}[ht!]
\begin{center}
\hspace*{-1.0cm}
\includegraphics[scale=0.42]{rproc_eu_abund.pdf}
\caption{\label{fig:abund_eufeh} Derived [Eu/Fe] abundances, as a function of [Fe/H], for the $r$-I, $r$-II, limited-$r$, $r/s$, as well as non-RPE stars. Also shown are the stars from the RPA-1 \citep{rpa1} and RPA-2 \citep{rpa2} samples, as well as $r$-I\ and $r$-II\ stars from the literature (extracted from JINAbase \citealt{jinabase2018}).
The arrows display upper limit measurements. Statistics of the RPA-3 stars are shown by marginal red histograms for both the [Eu/Fe] and [Fe/H] ranges, and those of RPA-1 and RPA-2 stars are shown by marginal gray histograms. $r$-I\ and $r$-II\ stars from the  literature are shown by open black histograms. Literature values are taken from \citet{westin2000}, \citet{johnson2002}, \citet{cayrel2004}, \citet{honda2004}, \citet{barklem2005}, \citet{ivans2006}, \citet{mashonkina2010}, \citet{SCM2014}, and \citet{roederer2014b}. Upper limits from the literature data have been suppressed for clarity.}
\end{center}
\end{figure}

\begin{figure}[ht!]
\begin{center}
\hspace*{-0.5cm}
\includegraphics[scale=0.42]{rproc_baeu_abund.pdf}
\caption{\label{fig:abund_baeu} Same as in Figure~\ref{fig:abund_eufeh}, but for [Ba/Eu] abundance ratios. Dividing dotted lines indicating differences between the $s$-process-enhanced stars ([Ba/Eu] $ > +$0.5), $r/s$-enhanced stars ($0.0 < $ [Ba/Eu] $ <$ +0.5), and $r$-process-enhanced stars ([Ba/Eu] $ < $ 0.0) abundance ratios are also shown.
}
\end{center}
\end{figure}

\begin{figure}[ht!]
\begin{center}
\hspace*{-1.0cm}
\includegraphics[scale=0.42]{rproc_eu_abs_abund.pdf}
\caption{\label{fig:eu_abs} Absolute Eu abundances, as a function of [Fe/H], for the $r$-II, $r$-I, limited-$r$, $r/s$, and $s$-process-enhanced stars identified in this work.}
\end{center}
\end{figure}

\begin{figure*}[ht!]
\begin{center}
\hspace*{-2.0cm}
\includegraphics[scale=0.225]{rpa_rel-abund.pdf}
\caption{\label{fig:abund_rel}Light, $\alpha$ and fe-peak abundance ratios (relative to Fe) of the target stars, as a function of [Fe/H], for the different $r$-process\ sub-classifications. Literature studies (shown as gray points) by \citet{cayrel2004}, \citet{barklem2005}, and \citet{cohen2013}, and \citet{roederer2014b} are also shown for comparison. Stars with distinguished abundance patterns (see discussion in Section\,\ref{sec:elem_abund}), 2MASS~J00101758$-$1735387 and 2MASS~J05241392$-$0336543, are also shown by the black star and yellow circle, respectively.}
\end{center}
\end{figure*}

\subsection{CO and the Light-Element Abundances}\label{sec:light-el_abund}

We derive the C abundances via spectrum synthesis of the molecular CH $G$-band at 4313\,{\AA} \citep{masseron2014}. As most of our stars lie on the red giant branch (see Figure\,\ref{fig:isoch}), the carbon abundances are assumed to be depleted by the CN cycle having operated during the stars' evolution. Relative [C/Fe] ratios of the target stars were therefore corrected for evolutionary status following \citet{placco2014}. We show both uncorrected and corrected values in Table\,\ref{tab:rproc_clas}. Additionally, we derive abundances of the light elements O and Na, using the forbidden \oi\ lines at 6300 and 6363\,{\AA} for O (whenever detectable), and the \nai\ doublet at 5889 and 5895\,{\AA} for Na, if the lines are sufficiently weak. We note that \citet{lind2011} reported non-significant NLTE abundance corrections of $<0.1$\,dex for the \nai\ doublet at our metallicity range. The [O/H] and [Na/H] abundances are listed in Table\,\ref{tab:abund_light}.

\subsection{The $\alpha$-element, K, and Al abundances}\label{sec:alpha-el_abund}
We derive $\alpha$-element abundances for \mgi, \cai, \tii, and \tiii\ from the EWs of their corresponding lines.  Potassium abundances were determined when the \ki\ lines at 7664 and 7698\,{\AA} were found to be unpolluted by telluric lines.
For Al, we determine the abundances by synthesizing the blue resonance \ali\ line at 3961\,{\AA}, taking into account the C abundance.  We caution that this line can be prone to significant 3D and NLTE corrections, up to 0.6\,dex for our range of metallicities, as reported by \citet{nordlander2017}.

\subsection{Iron-Peak Element Abundances}\label{sec:fepeak_abund}
We also derive the Fe-peak abundances for our target stars using \scii, \vii, \crii, \mni, \coi, \nii\ and \zni\ from their EWs using a multitude of lines, whenever available. We did not derive abundances from \vi\ and \cri\ lines, as they are prone to larger NLTE effects than their corresponding ionized species, \vii\ and \crii\ \citep{bergemann2010}.  We note that recent results by \citet{bergemann2019} showed that \mni\ lines can be prone to large NLTE and 3D effects, which they found to be increasingly larger for lower metallicities. For our sample of stars in this work, with [Fe/H] ranging typically between $-3.0<\mathrm{[Fe/H]}<-2.0$, the authors report 3D, NLTE corrections between 0.3 and 0.6\,dex for \mni\ lines. We do not include these corrections in this paper, but anticipate to include them in future papers once the correction grids have been made publicly available.

\section{Discussion}\label{sec:disc}

\subsection{$R$-Process-Enhanced Stars}
Following the neutron-capture sub-classifications defined in \citet{Beers2005} and \citet{frebel-rev2018}, we identify in this work 10 new $r$-II, 62 new $r$-I, and 17 new limited-$r$\ stars. Additionally, we find one $s$-process-enhanced and five $r/s$\ stars.
Overall, 46 stars from our target sample of 141 stars are non-$r$-process-enhanced (hereafter non-RPE, following the definition in \citealt{rpa2}). This shows that our overall selection procedure works well. One of our non-RPE stars, 2MASS~J00452379$-$2112161, is a rediscovery of the ultra metal-poor star CD$-38\,^{\circ}245$ \citep{bessell1984}.
Seventeen of the 141 stars were found to be Carbon-enhanced metal-poor (CEMP) stars, one enhanced at the $r$-II\ level, four at the $r$-I\ level, and three are limited-$r$\ stars. Additionally, three out of five of the $r/s$-enhanced stars are CEMP stars, as well as the $s$-process-enhanced star. Five stars were also found to be CEMP-no stars.

This brings up the total number of $r$-process-enhanced stars (at the $r$-II, $r$-I, limited-$r$, and $r/s$\ levels) identified by the RPA efforts from this work, \citet{rpa1}, \citet{rpa2}, \citet{roederer2018a}, \citet{holmbeck2018}, \citet{cain2018}, \citet{gull2018}, and \citet{placco2017} to 253 out of a total of 381 stars. Of those, a total of 28 stars are enhanced at the $r$-II\ level (7.4\%), 163 stars at the $r$-I\ level (42.7\%), 56 stars at the limited-$r$\ level (14.7\%), and 6 stars are enhanced at the $r/s$\ level (1.6\%). We caution that these values should not be used as physically interpretable estimates of the percentage of r-process enhanced stars.  They likely represent upper bounds as some stars were selected for publication on account of their high level of r-process enhancement \citep{placco2017,cain2018,gull2018,holmbeck2018}, and other sample selection biases have not yet been fully quantified.  The RPA intends to report more definitive estimates of the percentage of r-process enhanced stars in future work.

Figure~\ref{fig:hist} shows the cumulative histograms of newly identified $r$-I, $r$-II, and limited-$r$\ stars in this work, as a function of metallicity (upper panels) and V magnitudes (lower panels). None of those were previously identified in the RPA-1 and RPA-2 samples (open histograms; \citealp{rpa1,rpa2}).
Similar to the results found in the RPA-1 and RPA-2 samples (also shown on the plots), we find  $r$-I\ and $r$-II\ stars over a wide metallicity range, from \AB{Fe}{H}\ = $-3.5$ to \AB{Fe}{H}\  = $-1.2$. The limited-$r$\ stars are found primarily between \AB{Fe}{H}\ $ =-2.8$ and \AB{Fe}{H} =$-1.5$.
Our $r$-I\ stars cover a wider range of metallicities than those of the RPA-1 and RPA-2 samples, with one $r$-I\ star identified at \AB{Fe}{H}\ $ = -3.4$, and another at \AB{Fe}{H} $ = -1.2$. A similar cutoff in the $r$-I\ stars distribution is found around \AB{Fe}{H}$=-3.0$.
The mean metallicity of our newly identified $r$-II\ stars is \AB{Fe}{H}\ $\sim-2.3$, which agrees with the RPA-2 results. The \AB{Fe}{H}\ distribution from the RPA-1 sample is slightly shifted toward more metal-poor values, which can be explained by their implementation of an LTE analysis as compared to the 1D, NLTE and $<$3D$>$, NLTE analyses implemented in this work and RPA-2 papers respectively.
Our sample of stars also shift toward higher V magnitudes than those analyzed by the RPA-1 and RPA-2 efforts (Figure\,\ref{fig:hist}, lower panels). This shift correlates with the sizes of the telescopes used in this work and those by RPA-1 and RPA-2 efforts. Our $r$-II\ stars are particularly extended up to $V=14.2$ dex.

We also show the [Eu/Fe] abundance distribution, as a function of metallicity, in Figure~\ref{fig:abund_eufeh}. For comparison, we also show the stars from the RPA-1 and RPA-2 samples, as well as those from literature studies by \citet{westin2000}, \citet{johnson2002}, \citet{cayrel2004}, \citet{honda2004}, \citet{barklem2005}, \citet{ivans2006}, \citet{mashonkina2010}, \citet{SCM2014}, and \citet{roederer2014b}. The histograms on the top and right axes represent the number of stars covering each abundance range. Our results broadly overlap with the literature values, where the [Eu/Fe] ratios also cover a wide range of metallicities.
This reiterates the interpretation by \citet{sneden2000}, who suggested multiple formation sites for $r$-process-enhanced stars and a need for investigating $r$-process\ production and its site(s). We also plot the absolute Eu abundances derived for our target stars, as a function of \AB{Fe}{H}\, in Figure\,\ref{fig:eu_abs}. A clear increase of Eu with Fe is seen for the $r$-II, $r$-I, limited-$r$, and non-RPE stars. A similar trend was found by the RPA-1 survey \citep{rpa1}.

We additionally show the [Ba/Eu] ratios, as a function of \AB{Fe}{H}, in Figure\,\ref{fig:abund_baeu}. This can be used as a diagnostic for the level of $r$- versus $s$-process origin of the neutron-capture
elements measured in stars \citep{frebel-rev2018}. We also plot the  RPA-1 and RPA-2 stars on the same plot for comparison. Just as in previous work, the majority of newly identified $r$-II\ stars strictly have [Ba/Eu] $ < -0.5$, indicating a pure $r$-process\ origin, while the $r$-I\ and limited-$r$\ stars show larger variations in [Ba/Eu].

\subsection{Other Elemental Abundances}\label{sec:elem_abund}

The abundance ratios of the light elements, $\alpha$-elements, and Fe-peak elements relative to Fe are shown in Figure\,\ref{fig:abund_rel}, compared with literature values from \citet{cayrel2004}, \citet{barklem2005}, and \citet{cohen2013}. Classes of $r$-II, $r$-I, limited-$r$, and non-RPE star are distinguished to investigate possible systematic differences.

[C/Fe] and [O/Fe] ratios exhibit increasing trends toward lower metallicities.  [Mg/Fe], [Si/Fe], and [Ca/Fe] ratios show typical [$\alpha$/Fe] enhancements, of $\sim +0.4$, on average, while the [Al/Fe] ratios exhibit sub-solar values, except for one limited-$r$\ and four $r$-I\ stars. We note that [Ti/Fe] ratios  derived from the \tiii\ lines are shown, as they are less prone to NLTE effects than \tii\ lines \citep{sitnova2016}.

Fe-peak element abundance ratios mostly agree with those from the literature.  [Co/Fe] and [Zn/Fe] abundances exhibit increasing trends toward lower metallicities, in agreement with literature results (e.g., \citealt{Ezzeddine2019}), whereas [Mn/Fe] abundances show decreasing trends toward lower metallicities.

The abundances of the majority of our stars agree with typical Galactic field metal-poor star abundances, with the exception of the $r$-II\ star 2MASS~J00101758$-$1735387 (shown by black star symbol in Figure\,\ref{fig:abund_rel}),
which exhibits sub-solar and lower abundance ratios than the rest of the sample stars for C, Na, Mg,  Al, K, Ca, Sc, Ti, V, Cr, Mn, Co, Ni, and Zn. This star exhibits similar abundance patterns to other halo $r$-I\ and $r$-II\ stars identified in \citet{sakari2019} and \citet{xing2019}. These stars exhibit similar abundance patters to those identified in Milky Way (MW) dwarf galaxies \citep{mcwilliam2018}. The former studies suggested that such halo stars have an accretion origin from a dwarf galaxy previously enriched by a NS merger event, which enriched the gas from which 2MASS~J00101758$-$1735387 was formed. Additional investigation of its full $r$-process\ abundance pattern is underway, and will be presented in a future paper.

Another $r$-II\ star, \mbox{2MASS~J05241392$-$0336543}
exhibits a substantially enhanced [Na/Fe] ratio ($+$1.46).
As shown in Figure~\ref{fig:abund_rel} by the yellow circle marker, the [Mg/Fe] ratio in this star ([Mg/Fe]~$= +$0.09)
is slightly lower than most other stars in the sample, but all other light-element abundance ratios in \mbox{2MASS~J05241392$-$0336543} are typical
for its metallicity (\AB{Fe}{H} $= -2.50$).

The atypical [Na/Fe] and [Mg/Fe] abundance ratios appear unrelated to the light-element abundance variations typically observed in globular cluster stars, and discussion of its abundance anomalies and their origin is deferred to a future paper.

\section{Summary and Conclusions}\label{sec:conc}

In this paper we have presented high-resolution spectroscopic observations for 148 newly studied metal-poor stars. We identify 32 new stars with $-2.0\leq$\,[Fe/H]\,$<-1.0$, 103  stars with $-3.0\leq$\,[Fe/H]\,$<-2.0$, 5 stars with $-4.0\leq$\,[Fe/H]\,$<-3.0$, as well as one new ultra metal-poor star at [Fe/H]\,$<-4.0$. One of the stars is a rediscovery of the  previously observed ultra metal-poor star,  CD$-38\,^{\circ}245$. This is part of the ongoing RPA effort to identify new $r$-process-enhanced stars in the Galaxy. Seven stars were found to be double-lined spectroscopic binaries, and were thus excluded from further study, leaving the sample to contain 141 newly studied stars.

We determine 1D, NLTE atmospheric stellar parameters for all stars, and derive abundances for the light elements, including C and O, the $\alpha$-elements, the Fe-peak elements, and the neutron-capture elements, whenever suitable spectral lines are available. Among the 141 stars, we find 10 new $r$-II, 62 new $r$-I, and 17 new limited-$r$\ stars, based on their Eu, Ba, and Sr abundances. Additionally, we identified one new $s$-process-enhanced star and five new $r/s$\ stars.

Of the ten $r$-II\ stars, one is a CEMP star,  another is found to exhibit sub-solar elemental abundance ratios relative to Fe for many of the light, $\alpha$, and Fe-peak elements, consistent with accretion from a dwarf galaxy. All of the $r$-II\ stars exhibit consistently low ratios of [Ba/Eu] ($<-0.5$), which is a signature of pure $r$-process\ enhancement.
The $r$-II\ and $r$-I\ stars span a range of metallicities, from $-3.0$ to $-1.5$ in [Fe/H], with one $r$-I\ star detected at \AB{Fe}{H}\,$=-3.4$ and another at \AB{Fe}{H}\,$=-1.2$.

The neutron-capture and other elemental abundance ratios were compared to literature values, including the results from the RPA-1 and RPA-2 pilot samples of \citet{rpa1} and \citet{rpa2}, respectively. The metallicity distribution of our stars derived using 1D, NLTE agrees with those derived for the RPA-2 sample \citep{rpa2} using $<$3D$>$, NLTE. The RPA-1 \AB{Fe}{H}\ distribution, determined using 1D, LTE, is slightly shifted to the metal-poor end as compared to our sample, due to NLTE effects. Our sample is also shifted toward less bright stars as compared to the RPA-1 and RPA-2 sample. Our $r$-II\ stars cover a range of magnitudes, centered roughly around $V=13$.
The results from the three RPA releases from the Southern Hemisphere (this work, and \citealp{rpa1}) and the Northern Hemisphere \citep{rpa2} have significantly increased the numbers of known $r$-II, $r$-I, and limited-$r$\ stars, roughly doubling the previously known samples to
a total of 28 $r$-II\ stars, 163 $r$-I\ stars, and 56 limited-$r$\ stars, as well as 6 stars enhanced at the $r/s$\ level.
The RPA has already obtained snapshot spectroscopy for $\sim 1500$ bright metal-poor stars, and new snapshot spectra for $\sim 1000$ additional candidate stars are presently being obtained with a variety of telescopes in the Northern and Southern Hemispheres.  These observations are planned to be publicly released, as the data are reduced and analyzed, over the course of the next several years.

The full sample obtained by the RPA will be used to establish what fraction of accreted dwarf galaxies may have contributed metal-poor stars to the halo of the Milky Way; some partial information is already feasible to obtain \citep{roederer2018}.  In addition, detailed information on the nature of the environment in which the $r$-process events associated with individual stars occurred is being explored by consideration of the association of $r$-process-enhanced stars with the so-called ``Dynamically Tagged Groups" (DTGs) of very metal-poor stars (Yuan et al., submitted; Gudin et al., in preparation), following the pioneering work of \citet{roederer2018}. When ultimately combined with new accelerator results on the fundamental nuclear properties for isotopes specific to the $r$-process\ (using, e.g., the Facility for Rare Isotope Beams; FRIB, expected 2022 \citealt{FRIB2019}),
stellar abundances will provide key missing information required to answer many of the outstanding questions regarding the nature and astrophysical sites of the $r$-process.

\acknowledgements
We thank Gary Da Costa for providing the SkyMapper targets.
R.E., A.F., V.M.P, T.C.B., and I.U.R. acknowledge support from JINA-CEE (Joint Institute for Nuclear Astrophysics - Center for the Evolution of the Elements), funded by the NSF under Grant No. PHY-1430152. I.U.R. acknowledges partial support from NSF grants AST-1613536 and AST-1815403. A.F. is partially supported by NSF-CAREER grant AST-1255160 and NSF grant 1716251. JM thanks FAPESP (2014/18100-4).
This work made use of NASA's Astrophysics Data System Bibliographic Services, and the SIMBAD database, operated at CDS, Strasbourg, France \citep{simbad}. This work has made use of data from the European Space Agency (ESA) mission
{\it Gaia} (\url{https://www.cosmos.esa.int/gaia}), processed by the {\it Gaia}
Data Processing and Analysis Consortium (DPAC,
\url{https://www.cosmos.esa.int/web/gaia/dpac/consortium}). Funding for the DPAC
has been provided by national institutions, in particular the institutions
participating in the {\it Gaia} Multilateral Agreement.

\facilities{Magellan-Clay (MIKE, \citealt{bernstein2003})}

\software{IRAF~\citep{iraf}, matplotlib~\citep{matplotlib}, MOOG~\citep{sneden1973,sobeck2011}, MULTI~\citep{Carlsson1986,carlsson1992}, MARCS~\citep{gustafsson1975,gustafsson2008}}


\section{Introduction}

  It is a well-known fact that there exists no continuous function $f:[0,1]\to \mathbb{R} $ such that every point $x\in \mathbb{R}$ has either two or zero preimages. A somewhat less well-known fact is that there exist continuous functions $f:[0,1]\to \mathbb{R} $ where the preimage $f^{-1}(x)$ has an even cardinality for all $x\in \mathbb{R}$, in fact, one can construct such a function where the set of cardinalities of preimages is $\{0,2,4\}$.

  \medskip

 \hspace{4cm}  \includegraphics[width=0.5\textwidth]{waveexamplepic}

  In this paper, we are interested in the question of what other sets may occur as the set of preimage cardinalities for a continuous function. Given a continuous function $f:[a,b]\to \mathbb{R}, a < b$ we let $$\Omega _f = \{n\in \mathbb{N}\cup \{\infty \} \ : \ \exists x\in \mathbb{R} \ s.t. \ |f^{-1}(x)| = n\}.$$ Our question is what subsets of $\mathbb{N}$ may occur as $\Omega _f$ where $f:[a,b]\to \mathbb{R}$ is a continuous function. Notice that by boundedness of a continuous function in the compact interval, it is necessary that $0\in \Omega _f$. We also study the case when $f$ belongs to a sub-class of the class of continuous functions.

   \medskip

  One can also ask what subsets $S$ of $\mathbb{N}\cup \{\infty \} $ occur as $\Omega _f $ but in this question if we let $\infty \in S$ (otherwise, we fall back to the previous question) then it becomes easy: for any non-empty subset $A\subseteq \mathbb{N}$ where $0\in A$, the set $S = A\cup \{\infty \}$ does occur as $\Omega _f$ for some continuous $f:[a,b]\to \mathbb{R}$.

  \medskip

  By re-scaling the domain, we may just consider functions over the interval $[0,1]$ in this question (instead of an arbitrary interval). For an arbitrary interval $[a,b]$, $C[a,b]$ will denote the space of all continuous functions $f:[a,b]\to \mathbb{R}$. In this space, we define a subset  $$\mathcal{F}([a,b]) := \{f\in C[a,b] \ : \ f^{-1}(x) \ \mathrm{is \ finite \ for \ all} \ x\in \mathbb{R}\}.$$ In the case of the interval $[0,1]$, we will just write $\mathcal{F}$ instead of $\mathcal{F}([0,1])$. For functions $f\in \mathcal{F}$, we can also write $$\Omega _f = \{n\in \mathbb{N} \ : \ \exists x\in \mathbb{R} \ s.t. \ |f^{-1}(x)| = n\}.$$ If $\Omega _f = S$, then we also say that $f$ is an $S$-function. Notice that a function $f\in \mathcal{F}$ has finitely many local maximums and local minimums, moreover, the end points 0 and 1 are either a local maximum or a local minimum.

  \medskip

  Given a subset $S\subseteq \mathbb{N}$ with $|S|\geq 2$, we let $$m_1(S)  = \min (S\backslash \{0\}) \  \mathrm{and} \ m_2(S) = \max S.$$ Here, we use a convention that the maximum of an unbounded subset of natural numbers equals $\infty $. Our main theorem is the following

  \medskip

  {\bf Theorem 1.} Let $S\subseteq \mathbb{N}$ such that $0\in S$ and $|S|\geq 2$. Then there exists  $f\in C[0,1]$ with $\Omega _f = S$ if and only if $m_2 \geq 2m_1 - 1$.

  \medskip

  With slight care, one can indeed make all these functions from the class $C^{\infty }$. An important tool in this construction is the use of functions which "wiggles" infinitely many times. Hence this construction would not work for real analytic functions (by our definition, a function $f:[0,1]\to \mathbb{R}$ is real analytic if it can be extended holomorphically to some domain of $\mathbb{C}$ which includes the unit real segment $[0,1]$; in the sequel, instead of "real analytic" we will simply use the term "analytic"). Motivated by this fact, we also study the same question for all analytic functions, in particular, for polynomials. In fact, by far, not every subset of type $\Omega _f$ for a continuous function occurs as $\Omega _p$ where $p$ is a polynomial. First of all, for a polynomial $p(x)$, we necessarily have $\Omega _p$ is finite (same is true for analytic functions $f:[0,1]\to \mathbb{R}$). But even among finite subsets, for example, the set $\{0, 2, 4\}$ already does not occur as $\Omega _p$. By the following theorem, we achieve a complete characterization of the sets $\Omega _p$ as well where $p(x)$ is an analytic function. To state this theorem, we need to introduce some notations. Let $$A_{+} = \{(n, n+k, n+2k)\in \mathbb{N}^3 \ : \ n\geq 0, k > 0\},$$ \ $$A_{-} = \{(n, n-k, n-2k)\in \mathbb{N}^3 \ : \ n\geq 0, 0 < k\leq \frac{n}{2}\},$$ \
  $$B_{+} = \{(n, n+k, n+2k-1)\in \mathbb{N}^3 \ : \ n\geq 0, k \geq 0\},$$ \ $$ B_{-} = \{(n, n-k, n-2k-1)\in \mathbb{N}^3 \ : \ n\geq 0, 0 \leq k\leq \frac{n-1}{2}\},$$ \
  $$C_{+} = \{(n, n+k, n+2k-2)\in \mathbb{N}^3 \ : \ n\geq 0, k \geq 0\},$$ \ $$C_{-} = \{(n, n-k, n-2k-2)\in \mathbb{N}^3 \ : \ n\geq 0, 0 \leq k \leq \frac{n}{2}-1\},$$ \ $$A = A_{+}\cup A_{-}, B = B_{+}\cup B_{-}, C = C_{+}\cup C_{-}. $$

  \bigskip

  {\bf Theorem 2.} Let $S\subseteq \mathbb{N}$ such that $0\in S$. Then there exists an analytic function $f:[0,1]\to \mathbb{R}$ with $\Omega _f = S$ if and only if there exists a finite sequence $(x_0, x_1, x_2, \dots , x_{2m-1}, x_{2m})$ of natural numbers for some $m\geq 1$ such that the following conditions hold:

   i) $x_i\in S$ for all $i\in \{0, \dots, 2m\}$, and for all $z\in S$ there exists $i\in \{0, \dots, 2m\}$ such that $z = x_i$,

   ii) $x_0 = 0$ and $x_{2m} = 0$,

   iii) for all $0\leq i\leq m-1, (x_{2i}, x_{2i+1}, x_{2i+2})\in A\cup B\cup C$,

   iv) either $|\{i \ : \ 0\leq i\leq m-1, (x_{2i}, x_{2i+1}, x_{2i+2})\in B\}| = 2$ or $|\{i \ : \ 0\leq i\leq m-1, (x_{2i}, x_{2i+1}, x_{2i+2})\in C\}| = 1$,

   v) if any triple belongs to $B$, then no triple belongs to $C$.

   \medskip

   {\bf Remark 1.} If $S$ satisfies conditions i)-v), then the analytic function $p$ with $\Omega _p = S$ is constructed by interpolation method, thus in the statement of the theorem, "analytic" can be replaced with "polynomial".

  \section{Proof of Theorem 1}

  In this section we will assume that $S$ is a subset of $\mathbb{N}$ such that $0\in S$ and $|S|\geq 2$. The symbols $I, I_1, I_2, J, J_1, J_2$ will always denote a set $\{1,\dots , n\}$ for some $n\geq 1$ or one of the sets $\mathbb{Z}_{+}, \mathbb{Z}$. Such sets will be called {\em index sets}. We will first consider the special case of Theorem 1 when $1\in S$. In this special case Theorem 1 states the following

  \medskip

  {\bf Proposition 1.} Let $S\subseteq \mathbb{N}$ such that $0, 1\in S$. Then there exists  $f\in C[0,1]$ with $\Omega _f = S$.

  \medskip

   To prove Proposition 1, we introduce some terminology.

   {\bf Definition 1.} Let $n\geq 1$. A function $f\in C([a,b])$ is called an {\em $n$-wave function} if there exists real numbers $m < M$ and a sequence $a = x_0 < x_1 < \dots x_{n-1} < x_n = b$ such that $f$ is strictly increasing in the interval $[x_i, x_{i+1}]$ for all even $i\in \{0, \dots , n-1\}$ and strictly decreasing in the interval $[x_i, x_{i+1}]$ for all odd $i\in \{0, \dots , n-1\}$, moreover, $f(x_i) = m$ for all even $i\in \{1, \dots , n\}$ and $f(x_i) = M$ for all odd $i\in \{0, \dots , n\}$. The restriction $f|_{[x_{i-1}, x_i]}$ is called the $i$-th leaf of $f$; the interval $[x_{i-1}, x_i]$ is called the support of this leaf of $f$.

   \medskip

   Along with $n$-wave functions we will also use the notion of an $n$-stair function.

   \medskip

    {\bf Definition 2.} Let $n\geq 2$ be a positive odd integer. A function $f\in C([a,b])$ is called an {\em $n$-stair function} if there exists a sequence $a = z_0 < z_1 < \dots z_{n-1} < z_n = b$ such that for all $i\in \{0, \dots , n-1\}$ the restriction $f|_{[z_i, z_{i+1}]}$ is a $k_i$-wave function for some positive odd integer $k_i\geq 1$, and for all $0\leq i < j\leq n-1$, $x\in [z_i, z_{i+1}], y\in [z_j, z_{j+1}]$, $f(x) \leq f(y)$. The sequence $(k_1, \dots , k_n)$ will be called {\em the signature of $f$}. The interval $[z_i, z_{i+1}]$ is called the $(i+1)$-st block of $f$.

   \medskip

    \hspace{4cm}  \includegraphics[width=0.5\textwidth]{5-stairfunction.png}

    In the above picture we have 5-stair function with signature $(5,1,3,7,1)$; so this function has 5 blocks, and the restriction to these blocks are 5-wave, 1-wave, 3-wave, 7-wave and 1-wave functions.

 \medskip

    We also need to define infinite stair functions. This is done similarly as follows

    {\bf Definition 3.} A function $f\in C([a,b])$ is called an {\em $\infty $-stair function} if there exists a two-sided sequence $a < \dots < z_{-1} < z_0 < z_1 < \dots  < b$ with $\displaystyle \mathop{\lim }_{n\to \infty }z_{-n} = a$ and $\displaystyle \mathop{\lim }_{n\to \infty }z_n = b$ such that for all $i\in \mathbb{Z}$ the restriction $f|_{[z_i, z_{i+1}]}$ is a $k_i$-wave function for some positive odd integer $k_i\geq 1$, and for all $i < j$, $x\in [z_i, z_{i+1}], y\in [z_j, z_{j+1}]$, $f(x) \leq f(y)$. The sequence $(k_i)_{i\in \mathbb{Z}}$ will be called {\em the signature of $f$}. If $k_i = k$ for all $i\in \mathbb{Z}$ (i.e. if the sequence $(k_i)_{i\in \mathbb{Z}}$ is constant), then we also say that $f$ has a signature $k$. If, for some positive odd $p, q$ and $n_0\in \mathbb{Z}$, $k_i = p$ for all  $i\leq n_0$ and $k_i = q$ for all  $i > n_0$ then we also say that $f$ has a signature $(p,q)$. The intervals $[z_i, z_{i+1}]$ are called {\em the blocks} of $f$.

    \medskip

     \hspace{4cm}  \includegraphics[width=0.5\textwidth]{signature}

     In the above figure an $\infty $-stair function of signature $(3,5)$ is represented.

     \medskip

   {\bf Definition 4.} Let $a\leq c < d\leq b$, $f\in C([a,b]), g\in C([c,d])$ such that $f(c) = g(c)$ and $f(d) = g(d)$. The function $h:[a,b]\to \mathbb{R}$ defined as $$h(x) = \left\{
	\begin{array}{ll}
		f(x)  & \mbox{if } x \in [a,c]\cup [d,b] \\
		g(x) & \mbox{if } x \in [c,d]
	\end{array}
\right.$$ is called {\em the amendment of $f$ by $g$}. More generally, if $g_i:K_i\to \mathbb{R}, i\in I$ is a collection of continuous functions such that the intervals $K_i = [c_i,d_i]$ have mutually disjoint interiors and $g_i(c_i) = f(c_i), g_i(d_i) = f(d_i)$ for all $i\in I$, then {\em the amendment of $f$ by the collection $(g_i)_{i\in I}$} is defined as $$h(x) = \left\{
	\begin{array}{ll}
		f(x)  & \mbox{if } x \in [a,b]\backslash \displaystyle \mathop{\cup }_{i\in I}K_i \\
		g_i(x) & \mbox{if } x \in K_i, i\in I
	\end{array}
\right.$$
   If $[c,d] = [a,b]$ then the amendment is also called {\em a replacement.}

   \medskip

   Notice that the amendment of continuous functions is still continuous. We will be interested in amending $n$-stair functions with $r$-stair functions for various values of $n, r$ (including infinite values) at different sub-intervals.

   \bigskip

   {\bf Proof of Proposition 1.} In the case $S = \{0,1\}$ the claim is obvious so we may assume that $\{0,1\}$ is a proper subset of $S$. Let $r = \min (S\backslash \{0,1\})$, $(S\backslash \{0,1\})\cap 2\mathbb{N} = \{r_i : i\in I_1\}$ and $(S\backslash \{0,1\})\cap (2\mathbb{N}+1) = \{s_i : i\in I_2\}$. We also define the functions $f:[0,1]\to [0,1], h:[0,1]\to [0,1]$ as $f(x) = x, \forall x\in [0,1]$ and $h(x) = \left\{
	\begin{array}{ll}
		2x  & \mbox{if } x \in [0,\frac{1}{2}] \\
		-x+\frac{3}{2} & \mbox{if } x \in [\frac{1}{2},1]
	\end{array}
\right.$  We consider the following cases.

   \medskip

   {\em Case 1.}  $(S\backslash \{0,1\})\cap 2\mathbb{N} = \emptyset $ (i.e. $I_1 = \emptyset $).

   Let $C_i, i\in I_2$ be mutually disjoint open intervals in $(0,1)$ such that $\displaystyle \mathop{\cup }_{i\in I_2}\overline{C_i}=[0,1]$. For all $i\in I_2$, we let $g_i:C_i\to C_i$ be a surjective $\infty $-stair function with signature $s_i$. Then we take the identity function $f:[0,1]\to [0,1]$ and for each $i\in I_2$ we amend $f$ with $g_i$ in the interval $C_i$. The amendment $\overline{f}:[0,1]\to \mathbb{R}$ satisfies the conditions $\mathrm{Range}(\overline{f}) = [0,1]$,  $\overline{f}(0) = 0, \overline{f}(1) = 1$ and $\Omega _{\overline{f}} = S$.

   \medskip

   {\em Case 2.} $2\in S$.

 Let $B_i, i\in I_2$ and $C_i, i\in I_1$ be mutually disjoint open intervals in $(0,1)$ such that $\displaystyle \mathop{\cup }_{i\in I_2}\overline{B_i}=[0,\frac{1}{4}]$ and $\displaystyle \mathop{\cup } _{i\in I_1}\overline{C_i}=[\frac{1}{4}, \frac{1}{2}]$.

  For all $i\in I_2$ we let $g_i:B_i\to 2B_i$ be a surjective $\infty $-stair function with signature $s_i$ and for all $i\in I_1$ we let $g_i:C_i\to 2C_i$ be a surjective $\infty $-stair function with signature $r_i-1$. Then in the intervals $B_i, i\in I_2$ and $C_i, i\in I_1$ we amend $h$ with $g_i$. The amendment $\overline{h}:[0,1]\to \mathbb{R}$ satisfies the conditions  $\mathrm{Range}(\overline{h}) = [0,1]$, $\overline{h}(0) = 0, \overline{h}(\frac{1}{2}) = 1, \overline{h}(1) = \frac{1}{2}, \overline{h}^{-1}(0) = \{0\}, \overline{h}^{-1}(1) = \{\frac{1}{2}\}, \overline{h}^{-1}(\frac{1}{2}) = \{\frac{1}{4}, 1\}$ and $\Omega _{\overline{h}} = S$.

   \medskip

   {\em Case 3.} $r$ is odd, $I_1\neq \emptyset , I_2\neq \emptyset $.

   Let $\phi :[0,\frac{1}{2}]\to [0,1]$ be an $\infty $-stair function with signature $(r,r-2)$ such that every block supporting an $r$-wave function lies in $[0,\frac{1}{4}]$ and every block supporting an $(r-2)$-wave function lies in $[\frac{1}{4}, \frac{1}{2}]$, moreover, $|\phi ^{-1}(\frac{1}{2})| = r-1$. Let $B_i, i\in J_2$ and $C_i, i\in J_1$ be mutually disjoint open blocks of $\phi $ such that $\displaystyle \mathop {\cup }_{i\in J_2}\overline{B_i}=[0,\frac{1}{4}]$ and $\displaystyle \mathop{\cup } _{i\in J_1}\overline{C_i}=[\frac{1}{4}, \frac{1}{2}]$.

   \medskip

   Let $u_i, i\in J_2$ be a leaf of $\phi $ in the block $B_i$, and $v_i, i\in J_1$ be a leaf of $\phi $ in the block $C_i$. Let also $a:J_1\to I_1$ and $b:J_2\to I_2$ be surjective maps (such maps exist because the sets $J_1$ and $J_2$ are both infinite). We replace $u_i, i\in J_2, v_i, i\in J_1$ respectively with $\infty $-stair functions $\overline{u_i}$ of signature $s_{b(i)}-r+1$ and $\overline{v_i}$ of signature $r_{a(i)}-r$. Then we amend $f$ with functions $\overline{u_i}, i\in J_2$ and $\overline{v_i}, i\in J_1$.

   \medskip

   {\em Case 4.} $r$ is even, $r\geq 4$, $I_1\neq \emptyset , I_2\neq \emptyset $.

   This case is very similar to the previous case with the only difference that $\phi :[0,1]\to [0, \frac{1}{2}]$ is assumed to be an $\infty $-stair function with signature $(r+1,r-3)$ such that every block supporting an $(r+1)$-wave function lies in $[0,\frac{1}{4}]$ and every block supporting an $(r-3)$-wave function lies in $[\frac{1}{4}, \frac{1}{2}]$, $|\phi ^{-1}(\frac{1}{2})| = r-1$, and the signatures of $\overline{u_i}, \overline{v_i}$ are $s_{b(i)}-r$ and $r_{a(i)}-r+3$ respectively.

   \medskip

   {\em Case 5.} $(S\backslash \{0,1\})\cap (2\mathbb{N}+1) = \emptyset $ (i.e. $I_2 = \emptyset $).

   \medskip

   Since $S\backslash \{0,1\}$ does not contain an odd element we have that $r$ is even. Let $H:[-\frac{1}{2},1]\to [0,1]$ defined as $h(x) = \left\{
	\begin{array}{ll}
		-x  & \mbox{if } x \in [-\frac{1}{2},0] \\
		h(x) & \mbox{if } x \in [0,1]
	\end{array}
\right.$

  and $\psi :[0,\frac{1}{2}]\to [0,1]$ be an $\infty $-stair function with signature $(r-1,r-3)$ such that $|\psi ^{-1}(\frac{1}{2})| = r-2.$ Let also $D_i, i\in J$ be mutually disjoint open blocks of $\psi $ such that $\displaystyle \mathop{\cup }_{i\in J}\overline{D_i}=[\frac{1}{4},\frac{1}{2}]$. Let $w_i, i\in J$ be a leaf of $\psi $ in the block $D_i$, and $c:J\to I_1$ be a surjective function. We replace $w_i, i\in J$ with $\infty $-stair function $\overline{w_i}$ of signature $r_{c(i)}-r+3$. Then we amend $H$ with functions $\overline{w_i}, i\in J$. The amendment $\overline{H}:[-\frac{1}{2},1]\to \mathbb{R}$ satisfies the conditions $\mathrm{Range}(\overline{H}) = [0,1]$, $\overline{H}(-\frac{1}{2}) = \overline{H}(1) = \frac{1}{2}, \overline{H}(0) = 0, \overline{H}(\frac{1}{2}) = 1$ and $\Omega _{\overline{H}} = S$. By re-scaling the domain of $\overline{H}$ we can obtain a function $F:[0,1]\to \mathbb{R}$ with $\Omega _F = S$.  \ $\square $

  \bigskip

  For the proof of Theorem 1, we need to introduce the following

  {\bf Definition 5.} Let $S\subseteq \mathbb{N}$ and $0\in S$. $S$ is called of type I if for all even $x\in S$ there exists an odd $y\in S$ such that $y > x$. If $S$ is not of type I then we say that $S$ is of type II.

  \medskip

   The above definition implies that if $S$ is finite, then $S$ is of type II if and only if $\max S$ is even.

   \bigskip

   Now we are ready for the proof of the main theorem.

   \medskip

   {\bf Proof of Theorem 1.} First we will prove the "only if part".

   \medskip

     If $(0,1)$ does not contain any local maximum point of $f$ then either $f$ is a strictly monotone function or $f$ has only one local minimum in $(0,1)$ and local maximums at 0 and 1. In the latter case $\Omega _f = \{0,1,2\}$ whereas in the former case $\Omega _f = \{0,1\}$; the inequality  $m_2\geq 2m_1-1$ holds in both cases.

    \medskip

     Thus we may assume that $(0,1)$ does contain a local maximum point. Similarly, we may assume that $(0,1)$ contains a local minimum point. Then, without loss of generality we may also assume that at least one of the ends 0, 1 is not a local maximum. Let $y_1$ be the biggest value of $f$ in $[0,1]$  such that $f^{-1}(y_1)$ contains a local maximum in $(0,1)$. If we have $l$ local maximums in $f^{-1}(y_1)\cap (0,1)$, then $|f^{-1}(y_1)| \leq l+1$, and for a sufficiently small $\epsilon > 0$, $|f^{-1}(y_1-\epsilon )| \geq 2l$; moreover, if $|f^{-1}(y_1)| = l+1$, then $|f^{-1}(y_1-\epsilon )| \geq 2l+1$. Thus $|f^{-1}(y_1-\epsilon )| \geq 2|f^{-1}(y_1)|-1$. Hence we obtain that $m_2\geq 2m_1-1$.

    The "if part" is proved in Proposition 1 when $1\in S$. So we will assume that $1\notin S$.  Let $S\cap 2\mathbb{Z} = \{r_i \ : \ i\in I_1\}$ and $S\cap (2\mathbb{Z}+1) = \{s_i \ : \ i\in I_2\}$. We also let $p = \min (S\backslash \{0\})$. Then $p\geq 2$.

     We will consider the following cases.

    {\em Case 1.} $S$ is of type I.

    We can find an odd $q\in S$ such that $q\geq 2p-1$. Let $C_i, i\in \mathbb{Z}$ be mutually disjoint open intervals such that $\displaystyle \mathop{\cup }_{i\in \mathbb{Z}}\overline{C_i} = [\frac{1}{4}, \frac{3}{4}]$ and $C_i$ is on the left of $C_j$ for all $i<j$ (i.e. for all $x\in C_i, y\in C_j, x < y$).

    We let $f:[\frac{1}{4}, \frac{3}{4}]\to [\frac{1}{4}, \frac{3}{4}]$ be a surjective $\infty $-stair function of signature $(\dots , 1, 2p-1, 1, 2p-1, 1, 2p-1, \dots )$ such that every block supported on $C_i$ for even $i\in \mathbb{Z}$ is a 1-wave function, and every block supported on $C_i$ for odd $i\in \mathbb{Z}$ is a $(2p-1)$-wave function. We also let $f_1:[0, \frac{1}{4}]\to [0, \frac{1}{4}]$ and $f_2:[\frac{3}{4},1]\to [\frac{3}{4},1]$ be surjective $(2p-1)$-wave functions.

    Let $u_j$ be leaf of $f_j, j\in \{1,2\}.$ Let also $v_i$ be a leaf of $f$ supported on the block $C_{2i+1}, i\in \mathbb{Z}$ and $w_i = f|_{C_{2i}}, i\in \mathbb{Z}$. So $w_i$ is a 1-wave function for all $i\in \mathbb{Z}$.
    For all $i\in \mathbb{Z}$, we replace $v_i$ with an $\infty $-stair function $\overline{v_i}$ of signature $(q-2p+2)$; for all $j\in \{1,2\}$ we also replace $u_j$ with an $\infty $-stair function $\overline{u_j}$ of signature $(q-2p+2)$. Furthermore, for all $i\in 2\mathbb{Z}$ we replace $w_i$ with an $\infty $-stair function $\overline{w_i}$ of signature $(r_{a(i)}-1,r_{a(i)}+1)$, and for all $i\in (2\mathbb{Z}+1)$ we replace $w_i$ with an $\infty $-stair function $\overline{w_i}$ of signature $(s_{b(i)})$ where $a:2\mathbb{Z}\to I_1, b:2\mathbb{Z}+1\to I_2$ are surjective maps.
    The blocks of $\overline{w_i}, i\in 2\mathbb{Z}$ are either $(r_{a(i)}-1)$-wave functions or $(r_{a(i)}+1)$-wave functions. For each $i\in 2\mathbb{Z}$, in the former case we replace a leaf of it with $\infty $-stair function of signature $(q-r_{a(i)}+2)$ whereas in the latter case e replace a leaf of it with $\infty $-stair function of signature $(q-r_{a(i)})$. Let $\overline{\overline{w_i}}$ be this amendment of $\overline{w_i}$.

    \hspace{4cm}  \includegraphics[width=0.5\textwidth]{Case1}

    Now we amend $F:[0,1]\to [0,1]$ with functions $\overline{u_j}, j\in \{1,2\}$ and $\overline{v_i}, i\in \mathbb{Z},  \overline{w_i}, i\in 2\mathbb{Z}+1$ and  $\overline {\overline{w_i}}, i\in 2\mathbb{Z}$ where $F$ is defined as $F(x) = \left\{
	\begin{array}{lll}
		f_1(x)  & \mbox{if } x \in [0, \frac{1}{4}] \\
		f(x) & \mbox{if } x \in [\frac{1}{4},\frac{3}{4}] \\
		f_2(x)  & \mbox{if } x \in [\frac{3}{4},1]
	\end{array}
\right.$  (The above picture represents $F$ when $p=3$).

  For the amendment $\overline{F}:[0,1]\to \mathbb{R}$ we have $\Omega _{\overline{F}} = S$.

  \medskip

   {\em Case 2:} $S$ is of type II.

   Let the intervals $C_i, i\in \mathbb{Z}$ be as in Case 1 and $a:2\mathbb{Z}\to I_1, b:2\mathbb{Z}+1\to I_2$ are surjective maps. We can find an even $q\in S$ such that $q-1\geq 2p-1$.

   We let $h:[\frac{1}{4}, \frac{3}{4}]\to [\frac{1}{4}, \frac{3}{4}]$ be a surjective $\infty $-stair function of signature $(\dots , 1, 2p-1, 1, 2p-1, 1, 2p-1, \dots )$ such that every block supported on $C_i$ for even $i\in \mathbb{Z}$ is a 1-wave function, and every block supported on $C_i$ for odd $i\in \mathbb{Z}$ is a $(2p-1)$-wave function. We also let $h_1:[0, \frac{1}{4}]\to [0, \frac{1}{4}]$ be a surjective $(2p-1)$-wave function and $h_2:[\frac{3}{4},\frac{7}{8}]\to [\frac{3}{4},\frac{7}{8}]$ be a surjective $(2p-2)$-wave function with $h_2(\frac{7}{8}) = \frac{3}{4}$.

   Let $u_j$ be the first leaf of $h_j, j\in \{1,2\}$, and $v_i$ be a leaf of $h$ supported on the block $C_{2i+1}, i\in \mathbb{Z}$ and $w_i = h|_{C_{2i}}, i\in \mathbb{Z}$.

   Let the last two leafs of $h_2$ be supported on the interval $D = (d_1, d_2)$. We take a function $\xi : D\to [\frac{3}{4},\frac{7}{8}]$ such that $\xi ^{-1}(\frac{3}{4}) = \{d_1,d_2\}, \xi ^{-1}(\frac{7}{8}) = \{d_3, d_4\}$ for some $d_3, d_4\in D$, (so $|\xi ^{-1}(\frac{3}{4})| = |\xi ^{-1}(\frac{7}{8})| = 2$) and $|\xi ^{-1}(y)| = 4$ for all $y\in (\frac{3}{4},\frac{7}{8})$. Then we amend $h_2$ with $\xi $; let $h_3$ be this amendment.

    For all $i\in \mathbb{Z}$, we replace $v_i$ with an $\infty $-stair function $\overline{v_i}$ of signature $(q-2p+1)$. We also amend $u_1$ with an $\infty $-stair function $\overline{u_1}$ of signature $(q-2p+1)$, and $u_2$ with an $\infty $-stair function $\overline{u_2}$ of signature $(q-2p+1)$. Finally, for all $i\in 2\mathbb{Z}$ we replace $w_i$ with an $\infty $-stair function $\overline{w_i}$ of signature $(s_{a(i)}-2,s_{a(i)})$, and for all $i\in (2\mathbb{Z}+1)$ we replace $w_i$ with an $\infty $-stair function $\overline{w_i}$ of signature $(r_{b(i)}-1)$.

    The blocks of $\overline{w_i}, i\in 2\mathbb{Z}$ are either $(s_{a(i)}-2)$-wave functions or $(s_{a(i)})$-wave functions. For each $i\in 2\mathbb{Z}$, in the former case we replace a leaf of it with $\infty $-stair function of signature $(q-s_{a(i)}+2)$ whereas in the latter case e replace a leaf of it with $\infty $-stair function of signature $(q-s_{a(i)})$. Let $\overline{\overline{w_i}}$ be this amendment of $\overline{w_i}$.

    Now let $H:[0,1]\to [0,1]$ be defined as $H(x) = \left\{
	\begin{array}{llll}
		h_1(x)  & \mbox{if } x \in [0, \frac{1}{4}] \\
		h(x) & \mbox{if } x \in [\frac{1}{4},\frac{3}{4}] \\
		h_3(x)  & \mbox{if } x \in [\frac{3}{4},\frac{7}{8}] \\
		-6x+6 & \mbox{if } x \in [\frac{7}{8},1]
	\end{array}
\right.$

  We amend $H$  with functions $\overline{u_j}, j\in \{1,2\}$ and $\overline{v_i}, \overline{w_i}, i\in 2\mathbb{Z}+1$ and  $\overline {\overline{w_i}}, i\in 2\mathbb{Z}$. For the amendment $\overline{H}:[0,1]\to \mathbb{R}$ we have $\Omega _{\overline{H}} = S$.

  \ $\square $

    \bigskip

    {\bf Remark 2.} In the proof of Proposition 1 and Theorem 1, we use  $\infty $-stair functions $\sigma :[a,b]\to [c,d]$. We can allow "the heights of the waves" to be small enough and fast converging to zero so we can build such a function $\sigma _1:[a,b]\to [c,c+\delta ]$ of class $C^{\infty }$, where $\delta $ is some (possibly small) positive number and $\sigma ^{(n)}(a) = \sigma ^{(n)}(b) = 0$ for all $n\geq 1$. Then we re-scale the range by defining $\sigma _2(x) = \frac{c}{\delta }(\delta + c - d)+\frac{d-c}{\delta }\sigma _1(x)$ thus make $\sigma _2$ of class $C^{\infty }$. This allows to make the function $f$ in the claim of Proposition 1 and Theorem 1 from the $C^{\infty }$ class.

    \bigskip

    \section{Proof of Theorem 2}

    Let us first recall a well-known fact that if an analytic function has infinitely many zeros on a compact interval then it is identically zero. Thus an analytic function $p:[0,1]\to \mathbb{R}$ has finitely many zeros and  finitely many critical points. It also belongs to $\mathcal{F}$. Then there exists a finite sequence $y_1 < y_2 < \dots < y_{m-1}$ such that any local maximum or local minimum point of $p$ belongs to $\displaystyle \mathop{\sqcup }_{1\leq i\leq m-1}p^{-1}(y_i)$, and each of the subsets $p^{-1}(y_i)$ contains at least one local maximum or local minimum. Let $y_0 = -\infty $ and $y_m = \infty $, and let $\phi :\mathbb{R} \to \mathbb{N}$ be defined as $\phi (y) = |p^{-1}(y)|$. Then in any interval $(y_i, y_{i+1}), 0\leq i\leq m-1$, $\phi $ is constant. Let $z_i$ be a point in $(y_i, y_{i+1}), 0\leq i\leq m-1$.

    \medskip

    We let $x_{2i} = |p^{-1}(z_i)|, 0\leq i\leq m$,  and $x_{2i+1} = |p^{-1}(y_i)|, 1\leq i\leq m$. Then for all $0\leq i\leq m-1$,

    if $p^{-1}(y_i)$ does not contain end points, contains $l$ local minimums and $r$ local maximums, and $r < l$, then  $(x_{2i}, x_{2i+1}, x_{2i+2})\in A_{+}$;

    if $p^{-1}(y_i)$ does not contain end points, contains $l$ local minimums and $r$ local maximums, and $r > l$, then  $(x_{2i}, x_{2i+1}, x_{2i+2})\in A_{-}$;

    if $p^{-1}(y_i)$ contains one end point as a local minimum, $l$ local minimums and $r$ local maximums, and $r \leq l$, then  $(x_{2i}, x_{2i+1}, x_{2i+2})\in B_{+}$;

    if $p^{-1}(y_i)$ contains one end point as a local minimum, $l$ local minimums and $r$ local maximums, and $r > l$, then  $(x_{2i}, x_{2i+1}, x_{2i+2})\in B_{-}$;

    if $p^{-1}(y_i)$ contains one end point as a local maximum, $l$ local minimums and $r$ local maximums, and $r \leq l$, then  $(x_{2i}, x_{2i+1}, x_{2i+2})\in B_{+}$;

    if $p^{-1}(y_i)$ contains one end point as a local maximum, $l$ local minimums and $r$ local maximums, and $r > l$, then  $(x_{2i}, x_{2i+1}, x_{2i+2})\in B_{-}$;

      if $p^{-1}(y_i)$ contains two end points of opposite types, $l$ local minimums and $r$ local maximums, and $r \leq l$, then  $(x_{2i}, x_{2i+1}, x_{2i+2})\in C_{+}$;

      if $p^{-1}(y_i)$ contains two end points of opposite types, $l$ local minimums and $r$ local maximums, and $r > l$, then  $(x_{2i}, x_{2i+1}, x_{2i+2})\in C_{-}$;

      if $p^{-1}(y_i)$ contains two end points both of which are local maximums, $l$ local minimums and $r$ local maximums, and $r > l+2$, then  $(x_{2i}, x_{2i+1}, x_{2i+2})\in C_{-}$;

      if $p^{-1}(y_i)$ contains two end points both of which are local maximums, $l$ local minimums and $r$ local maximums, and $r \leq l+2$, then  $(x_{2i}, x_{2i+1}, x_{2i+2})\in C_{+}$;

      if $p^{-1}(y_i)$ contains two end points both of which are local minimums, $l$ local minimums and $r$ local maximums, and $l \geq r+2$, then  $(x_{2i}, x_{2i+1}, x_{2i+2})\in C_{+}$;

      if $p^{-1}(y_i)$ contains two end points both of which are local minimums, $l$ local minimums and $r$ local maximums, and $l > r+2$, then  $(x_{2i}, x_{2i+1}, x_{2i+2})\in C_{-}$;

      in all other cases $x_{2i} = x_{2i+1} = x_{2i+2}$.

      \medskip

      Thus we proved the "only if part" of the theorem. To prove the "if part", we first observe the following simple fact.

      \medskip

      {\bf Lemma 1.} Let $f\in \mathcal{F}([a,b])$ and $y_1 < y_2 < y_3$ such that for all $y \in [y_1, y_3]\backslash \{y_2\}$, the set $f^{-1}(y)$ does not contain any local maximum or local minimum, and $|f^{-1}(y_1)| = n, |f^{-1}(y_2)| = n+k, |f^{-1}(y_3)| = n+2k$ for some $n\geq 0, k\in \mathbb{Z}\backslash \{0\}$. Let also $r, s\in \mathbb{Z}\backslash \{0\}$ such that $r+s = k$. Then $f$ can be amended in the set $f^{-1}([y_1, y_3])$ such that the amendment $F$ belongs to $\mathcal{F}([a,b])$, moreover, there exists $z_1, z_2 \in (y_1, y_3)$ such that $z_1 < y_2 < z_2$, for all $y \in [y_1, y_3]\backslash \{z_1, z_2, y_2\}$, the set $F^{-1}(y)$ does not contain any local maximum or local minimum, and $|F^{-1}(y_1)| = n, |F^{-1}(z_1)| = n+r, |F^{-1}(y_2)| = n+2r, |F^{-1}(z_2)| = n+2r+s, |F^{-1}(y_3)| = n+2k$.  $\square $

      \medskip

      Now, given any sequence $(x_0, x_1, x_2, \dots , x_{2m-1}, x_{2m})$ satisfying conditions i)-v) we have two cases: 1) None of the triples $(x_{2i}, x_{2i+1}, x_{2i+2}), 0\leq i\leq m-1$ belongs to $C$ or 2) Exactly one of the triples $(x_{2i}, x_{2i+1}, x_{2i+2}), 0\leq i\leq m-1$ belongs to $C$.

      In the former case we have $0\leq p < q\leq m-1$ such that $$(x_{2p}, x_{2p+1}, x_{2p+2}), (x_{2q}, x_{2q+1}, x_{2q+2})\in B$$ and in the latter case we have $0\leq k\leq m-1$ such that $(x_{2k}, x_{2k+1}, x_{2k+2})\in C$.

      We can take a sequence $(0, 1, \dots , m-1)$ and inductively build a piecewise linear function $f:[0,1]\to \mathbb{R}$ (i.e. $f$ is continuous, has finitely many singularity points, and in between these singularities it is defined as a polynomial of degree at most 1, such that

      (i) if $c$ is a local maximum or local minimum, then $f(c)\in \{0, 1, \dots , m-1\}$

      (ii) $f(a) = p, f(b) = q$ in Case 1), and $f(a) = f(b) = k$ in Case 2).

      (iii) For all $i\in \{0, 1, \dots , 2m\}$, $|f^{-1}(i)| = x_i$.

       \medskip

       Here, we carry the induction on the quantity $m$. For $m\leq 4$ we can see by a direct check that the claim holds. If $m\geq 5$ then there exists $i\in \{0, \dots , m-2\}$ such that none of the triples $(x_{2i}, x_{2i+1}, x_{2i+2}), (x_{2i+2}, x_{2i+3}, x_{2i+4})$ belongs to $B\cup C$. Then using Lemma 1 we reduce the question to the case of a sequence of length $2m-2$.

       \medskip

       Now, conditions (i)-(ii) imply that $\Omega _f = S$. But, by interpolation formula [1], for any piecewise linear function one can find a polynomial $p:[0,1]\to \mathbb{R}$ with same range and the same set of local maximums and local minimums. Then $\Omega _p = \Omega _f.$ \ $\square $

    \medskip

    {\bf Acknowledgement:} I thank Azer Akhmedov for his encouragements and for so many helpful comments and suggestions on the original draft.

    \vspace{1cm}

    REFERENCES:

    \bigskip

    [1] \ Whittaker, E. T. and Robinson, G. "Lagrange's Formula of Interpolation." §17 in The Calculus of Observations: A Treatise on Numerical Mathematics, 4th ed. New York: Dover, pp. 28-30, 1967.

    \vspace{1cm}

    E-mail: 147414@west-fargo.k12.nd.us

  \end{document}