Magnetic fields in central stars of planetary nebulae?

Abstract: Most of the planetary nebulae (PN) have bipolar or other non-spherically symmetric shapes. The presence of a magnetic field in the central star may be the reason for this lack of symmetry, but observational works published in the literature have so far reported contradictory results.We try to correlate the presence of a magnetic field with the departures from the spherical geometry of the envelopes of planetary nebulae.We determine the magnetic field from spectropolarimetric observations of ten central stars of planetary nebulae. The results of the analysis of the observations of four stars was previously presented and discussed in the literature, while the observations of six stars, plus additional measurements for a star previously observed, are presented here for the first time.All our determinations of magnetic field in the central planetary nebulae are consistent with null results. Our field measurements have a typical error bar of 150-300 G. Previous spurious field detections obtained with FORS were probably due to the use of different wavelength calibration solutions for frames obtained at different position angles of the retarder waveplate.Currently, there is no observational evidence for the presence of magnetic fields with a strength of the order of hundreds Gauss or higher in the central stars of planetary nebulae.

\section{Introduction}
There is still no conclusive theory why more than 80\,\%\ of the known planetary
nebulae (PNe) have bipolar and non-spherically symmetric structures
\citep{Zuckerman-Aller:86,Stanghellini-etal:93,CorradiSchwarz:95}. An overview
of possible mechanisms that shape PNe is given by \cite{BalickFrank02}. Several
of these processes imply the presence of magnetic fields which deflect the
outflow of the matter along the magnetic field lines. Recently,
\cite{ThirumalaiHeyl:10} have published model calculations for AGB stars
incorporating both magnetism and stellar winds with dust grains.

Such magnetic fields could be either fossil remnants from previous stages of
stellar evolution, or could be generated by a dynamo at the interface between a
rapidly rotating stellar core and a more slowly rotating envelope.
\cite{Blackman-etal:01} argue that some remnant field anchored in the core will
survive even without a convection zone, although the convective envelope may not
be removed completely.

The idea that magnetic fields are important ingredients shaping PNe has been
supported by the detection of SiO, H$_2$O, and OH MASER emission in
circumstellar envelopes of AGB stars pointing at milligauss fields in these
nebula
\citep{Kemball-Diamond:97,Szymczak-Cohen:97,Miranda-etal:01,Vlemmings-etal:02,Vlemmings-etal:05,Vlemmings-etal:06,Herpin-etal:06,Sabin-etal:07,Kemball-etal:09,
Gomez:09,Vlemmings:11}.  Moreover, using an idea of \cite{Pascoli:85},
\cite{HugginsManley:05} connected the extreme filamentary structures seen in
high-resolution optical images of certain planetary nebulae to magnetic fields
consistent with those measured in the MASER from  the precursor circumstellar
envelopes.

For the first time, and with the help of optical circular
spectropolarimetry carried out with the FORS1 spectrograph of the UT1
(``Antu'') telescope of the Very Large Telescope (VLT) of the European
Southern Observatory (ESO), \cite{Jordan-etal:05} reported on the
detection of magnetic fields in the central stars of the planetary
nebulae NGC\,1360 and LSS\,1362. For the central stars of EGB\,5
and Abell\,36 the existence of a magnetic field was found to be
probable but with less certainty.

\cite{Pascoli:08} pointed out that the magnetic field at the surface of
central stars of planetary nebulae is not necessarily connected to the magnetic
fields in the nebula itself; the latter can be  a fossil component of the
primary field embedded in the AGB star. Nevertheless, the reported detection of magnetic fields in central stars of
planetary nebulae has triggered several additional observational and theoretical
studies on the shaping of planetary nebula, e.g.\,by \cite{Garcia-Diaz-etal:08},
\cite{Tsui:08}, and \cite{PascoliLahoche:10} taking magnetic fields into
account.  On the other hand, \cite{Soker:06} cast strong doubts that magnetic
fields could be the main agent shaping planetary nebulae. He argued  that a
single star cannot supply the energy and angular momentum if the magnetic fields
have a large-scale structure required to shape the outflow from an AGB star.

Recently, the detection of magnetic fields in the central stars of
planetary nebulae (CPNs) was called into question by
\citet{Leone-etal:11}, who re-observed NGC\,1360 and LSS\,1362 with
the FORS2 instrument, and concluded that their effective magnetic
field is null within an uncertainty of $\sim 100$\,G (NGC\,1360), and
$\sim 290$\,G (LSS\,1362). Furthermore, both \citet{Leone-etal:11} and
\citet{Bagnulo-etal:12} re-analysed the observations previously
obtained with FORS1 by \citet{Jordan-etal:05}, and could not confirm
the original detection.

The conclusion based on the observations of four CPNs\ reduced by two
independent groups is that there is no observational evidence for
magnetic fields in the CPNs. The aim of this work is to enlarge the
sample of CPNs\ checked for magnetic field, in order to estimate the
occurrence of the magnetic field in CPNs\ and, if magnetic field is
detected, whether its presence correlates with the asymmetry of their
envelope. To achieve this goal, two of the teams that had presented
(discordant) results on previous FORS measurements have joint their
efforts to present here a more complete survey of magnetic fields of
CPNs\ obtained during a three-night observing run with FORS1 carried
out in 2005.

\begin{table*}
\input{Table_1}
\end{table*}

\section{Observations, data reduction, and magnetic field determinations}

All spectropolarimetric data reported in this paper were taken with
the FORS1 instrument \citep{Appenzeller-etal:98} of the ESO Very Large
Telescope. The polarimetric optics of FORS1, now moved to the twin
instrument FORS2, are based on the principle described by
\cite{Appenzeller:67}. FORS2 now is one of the few optical
spectropolarimeters available for the study of stellar magnetism. Due
to the large aperture of the telescope (8\,m) FORS2 is best suited for
the study of faint stars like white dwarfs
\citep{AznarCuadrado-etal:04,Jordan-etal:07}, sub\-dwarfs
\citep{Jordan-etal:05}, and CPNs\
\citep{Jordan-etal:05}.

The FORS1 dataset previously analysed by \citet{Jordan-etal:05},
\citet{Leone-etal:11}, and \citet{Bagnulo-etal:12} consists of
observations of \object{NGC\,1360}=\object{CD--26\,1340},
\object{EGB\,5}=\object{PN\,G211.9+22.6},
\object{LSS\,1362}=\object{PN\,G273.6+06.1}, and
\object{Abell\,36}=\object{ESO\,577-24}. These observations were
obtained by \cite{Jordan-etal:05}, in service mode between November
2, 2003, and January 27, 2004, using grism 600B, and a 0.8\arcsec\
slit width, for a spectral resolution of about 1000.

Three additional nights of telescope time were obtained in visitor
mode between June 4, 2005, and June 6, 2005, at the UT2 (``Kueyen'')
of the ESO VLT. The instrument setup included again grism 600\,B, but
a slit width of 0.7\arcsec, for a spectral resolution of about 1200.

The criterion for the selection of objects were the optimum visibility during the
observation nights and the brightness of the stars in order to reach an optimum signal-to-noise ratio.
Moreover, we made sure that at least two of the stars were in the  center of a nebula with
almost spherical shape (\object{LSE\,125} and  \object{Hen\,2-194})
so that in case of positive detections the magnetic field can be  correlated with the
topology of the nebulae.

During the visitor run, six new CPNs\
were observed two or three times: \object{HD\,107969}=\object{NGC\,4361},
\object{LSE\,125}=\object{PN\,G335.5+12.4},
\object{Hen\,2-194}=\object{ESO\,392-2},
\object{HD\,154072}=\object{IC\,4637},
\object{HD\,161044}=\object{IC\,1266}, and
\object{WD\,2226-210}=\object{NGC\,7293}.
\object{Abell\,36}=\object{ESO\,577-24}, already checked for magnetic
field in the previous service observing run, was re-observed three
times. Additionally, we took one spectropolarimetric dataset for
\object{HD\,160917}, a non-magnetic B9V comparison star.

Figure\,\ref{Fig_spectra} shows the summed up high-quality  spectra for all
CPNs. Only \object{Hen\,2-194} and \object{HD\,161044} were subject to
significant nebula emission.

We analysed all this old and new observational material
by adopting a method described in detail by \citet{Bagnulo-etal:12}.
Spectra are calibrated and extracted using the ESO FORS pipeline
\citep{Izzo-etal:10}, then combined using the difference
method to obtain both the reduced Stokes $V$ profiles (\ensuremath{P_V}) and the
null profiles (\ensuremath{N_V}), as described by \citet{Bagnulo-etal:12}.  The mean
longitudinal magnetic field \ensuremath{\langle B_z\rangle}\ was then calculated by using a
least-square technique based on the relationship
\begin{equation}
\ensuremath{P_V}(\lambda)= - g_\mathrm{eff} \ \ensuremath{C_Z} \ \lambda^{2} \
                \frac{1}{I(\lambda)} \
                \frac{\mathrm{d}I(\lambda)}{\mathrm{d}\lambda} \
                \ensuremath{\langle B_z\rangle}\;,
\label{Eq_Bz}
\end{equation}
where \ensuremath{P_V}\ is the reduced Stokes $V$ profile, $I(\lambda)$ is the Stokes $I$
profile of a spectral line, $g_\mathrm{eff}$ is the effective Land\'{e} factor, and
\begin{equation}
\ensuremath{C_Z} \simeq 4.67 \times 10^{-13}\,\AA^{-1}\ \mathrm{G}^{-1} \; .
\end{equation}
For more details, see, e.g., \citet{Bagnulo-etal:12}.
Compared with \citet{Bagnulo-etal:12}, in this work we have implemented
a sigma-clipping algorithm in the determination of the magnetic field from
the correlation diagram of circular polarisation against local flux
derivative. We have also calculated the null profiles \ensuremath{N_V}\ and compared their
oscillation about zero with the \ensuremath{P_V}\ error bars, and also measured the
null field \ensuremath{\langle N_z\rangle}, i.e., the magnetic field obtained by applying Eq.~(\ref{Eq_Bz})
to the \ensuremath{N_V}\ profile instead of \ensuremath{P_V}. The statistical significance of the
null field values is extensively discussed in \citet{Bagnulo-etal:12}. Here we
report that all null field values were found consistent with zero within error bars.
We also note that the targets of this survey are relatively faint, and that most of
the observations discussed here are not characterised by a ultra-high signal-to-noise
ratio. The main contributor to the error bars is thus photon-noise and background
subtraction (since our targets are embedded in a circumstellar envelope).

The original data reduction of the observations obtained within
Programme ID 072.D-0089 by \cite{Jordan-etal:05} was based on two
distinct wavelength calibration solutions for the frames obtained at
the two different position angles of the retarder plate adopted for
the science observations: science frames obtained with the retarder
waveplate at position angle $-45\degr$ were calibrated with
calibration frames obtained with the retarder waveplate at $-45\degr$,
and science frames obtained with the retarder waveplate at $+45\degr$
were calibrated with calibration frames obtained at
$+45\degr$. However, \citet{Bagnulo-etal:06} and
\citet{Bagnulo-etal:09} shown that this method prevents wavelength
calibration errors to cancel out, and may lead to spurious
detections. While this problem did not occurr in various tests
\citep[e.g.][]{Bagnulo-etal:02}, in the case of for instance the
075.D-0289 data it would lead to spurious magnetic field
measurements of the order of $\sim 9$\,kG.

In their analysis of the data obtained in 072.D-0089,
\cite{Jordan-etal:05} did compare the Stokes profiles obtained using a
common wavelength calibration with those obtained using two distinct
solutions, and concluded that both methods would produce similar
Stokes profiles. However, at that time they did not notice that,
although the profiles looked similar, field measurements were
different: after the analysis of the profiles obtained using two
wavelength calibration frames, a field was firmly detected, while
using a unique wavelength calibration frame, field detections would
disappear. We conclude that for that reason, the relatively high
magnetic fields reported by \cite{Jordan-etal:05} are spurious.

A similar conclusion was drawn by \citet{Bagnulo-etal:12} and
\citet{Landstreet-etal:12} for the measurements of magnetic fields in
subdwarfs by \citet{OToole-etal:05}: the kG magnetic fields are also
disappearing when using a single wavelength calibration frame for the
entire science dataset. We note instead that the detections of weak
magnetic fields in white dwarfs by \citet{AznarCuadrado-etal:04} were
basically confirmed by \citet{Bagnulo-etal:12}.

\section{Results}\label{Sect:Results}
\begin{figure}[t]
\includegraphics[width=0.5\textwidth,angle=0]{Fig_spectra.pdf}\\
\caption{\label{Fig_spectra}
Normalised spectra of all central stars of planetary nebulae observed in the  075.D-0289 campaign ordered by increasing effective temperature
from top to the bottom. The 072.D-0089 are shown in Fig.\,1 of \cite{Jordan-etal:05}.
}
\end{figure}

\begin{figure}[t]
\includegraphics[width=0.39\textwidth,angle=270]{Fig_Abell36.pdf}\\
\caption{\label{Fig_Abell36}
The observations of Abell\,36 were obtained with FORS1 on 2005-06-03.  The
top panel shows the observed flux $F$ (black solid line, in arbitrary
units, and not corrected for the instrument response), the
$\ensuremath{P_V} = V/I$ profile (red solid line centred about 0), and the null profile
\ensuremath{N_V}\ (blue solid line, offset by $-0.75$\,\% for display purpose). The null
profile is expected to be centred about zero and scattered according
to a Gaussian with $\sigma$ given by the \ensuremath{P_V}\ error bars. The
\ensuremath{P_V}\ error bars are represented with light blue bars centred about
$-0.75$\,\%.  The slope of the interpolating lines in the bottom panels
gives the mean longitudinal field from \ensuremath{P_V}\ (left bottom panel) and
from the null profile (right bottom panel) both calculated using the H
Balmer and metal lines. The corresponding \ensuremath{\langle B_z\rangle}\ and \ensuremath{\langle N_z\rangle}\ values are $-33
\pm 158$\,G and $-420 \pm 161$\,G, respectively.
}
\end{figure}

Table\,\ref{Tab_Observations} lists our field determinations from all
new FORS observations. We have also included new field determinations
from the observations carried out in service mode by
\citet{Jordan-etal:05}.  Note that the field estimates for the
observations obtained in 072.D-0089 slightly differ from those
published by \citet{Bagnulo-etal:12} because of the implementation of
the sigma-clipping algorithm, and also a slightly different choice of
the spectral points considered for field
determination. Figure\,\ref{Fig_Abell36} shows an example of
reduced data.

\section{Discussion}
Our analysis of ten CPNs\ does not show significant evidence for the
existence of longitudinal magnetic fields above a few hundred
Gauss. This is in contradiction with the result of
\cite{Jordan-etal:05}, who determined magnetic field strengths of
several kG in their analysis of four CPNs. Our current results are
consistent with the investigation of the central star of
NGC\,1360 by \citet{Leone-etal:11}, who determined
an upper limit for
the magnetic field of 300\,G \citep[while][reported a longitudinal
magnetic field of up to 3\,kG]{Jordan-etal:05}. For the central star of
LSS\,1362  \citet{Leone-etal:11} obtained an upper limit of 600\,G.
We conclude that a non-optimal wavelength calibration method has led
\citet{Jordan-etal:05} to a spurious detection of magnetic fields.

Our re-reduction of the 072.D-0089 data leads to a weighted mean
magnetic field (see Table\,\ref{Tab_Observations}) for NGC\,1360 of
$244\pm162$\,G, assuming that a possible (and weak) magnetic field
appears constant with time.

With our sample of ten stars (Abell\,36 has been observed in both
observation campaigns) we conclude that there is no confirmed case for
a magnetic field in the central star of a planetary nebuale at a kG
level. Magnetic fields of the order to 100-300\,G, however, cannot be
discarded. Indirect evidence for the existence of mG fields in
proto-planetary nebulae could still support an influence of magnetic
fields on the shape of PN. However, these fields need no be connected
to a field of the central star \citep{Pascoli:08}.

\acknowledgements{We thank the staff of the ESO VLT for carrying out
  the service observations of programme 072.D-0089 and supporting our
  visitor mode observations of programme 075.D-0289 (PI=S. Jordan for
  both programmes).}

\bibliographystyle{aa}

Solar and Stellar Astrophysics (astro-ph.SR)


A Census of Young Stars and Brown Dwarfs in IC 348 and NGC 1333

Abstract: We have obtained optical and near-infrared spectra of candidate members of the star-forming clusters IC 348 and NGC 1333. We classify 100 and 42 candidates as new members of the clusters, respectively, which brings the total numbers of known members to 478 and 203. We also have performed spectroscopy on a large majority of the previously known members of NGC 1333 in order to provide spectral classifications that are measured with the same scheme that has been applied to IC 348 in previous studies. The new census of members is nearly complete for Ks<16.8 at Aj<1.5 in IC 348 and for Ks<16.2 at Aj<3 in NGC 1333, which correspond to masses of <=0.01 Msun for ages of 3 Myr according to theoretical evolutionary models. The faintest known members extend below these completeness limits and appear to have masses of ~0.005 Msun. In extinction-limited samples of cluster members, NGC 1333 exhibits a higher abundance of objects at lower masses than IC 348. It would be surprising if the initial mass functions of these clusters differ significantly given their similar stellar densities and formation environments. Instead, it is possible that average extinctions are lower for less massive members of star-forming clusters, in which case extinction-limited samples could be biased in favor of low-mass objects in the more heavily embedded clusters like NGC 1333. In the H-R diagram, the median sequences of IC 348 and NGC 1333 coincide with each other for the adopted distances of 300 and 235 pc, which would suggest that they have similar ages. However, NGC 1333 is widely believed to be younger than IC 348 based on its higher abundance of disks and protostars and its greater obscuration. Errors in the adopted distances may be responsible for this discrepancy.

\section{Introduction}
\label{sec:intro}

A thorough census of young stars and brown dwarfs in nearby star-forming
regions is important for measuring the global properties of young stellar
populations (e.g., initial mass functions, disk fractions) and for providing
well-defined samples of targets for a variety of studies of star and planet
formation.
The Perseus molecular cloud is one of the nearest and richest of these regions
\citep[$\sim300$~pc,][]{sch14}.
It contains several hundred young stars, most of which reside in two clusters,
IC~348 and NGC~1333 \citep{her08,wal08}.
Age estimates for IC~348 based on evolutionary models range from
2 to 6~Myr \citep{luh03,bel13}.
NGC~1333 appears to be younger than IC~348 based on its greater obscuration
and higher abundance of circumstellar disks and protostars \citep{mue07,gut08}.

Candidate members of IC~348 have been identified via
proper motions \citep{fre56,sch99},
H$\alpha$ emission \citep{her54,her98},
X-ray emission \citep{pre96,pre01,pre02,pre04,for11,ste12},
outflow signatures \citep{wal06,hat09},
optical and near-infrared (IR) photometry
\citep{str74b,lad95,luh98,luh03,luh99,naj00,mai03,bur09,mue03,alv13},
mid-IR excess emission \citep{lad06,jor06,jor07,cie07,mue07,reb07,cur09,eva09,you15},
variability \citep{fla12,fla13,cod14},
and kinematics \citep{cot15}.
Optical and near-IR spectroscopy has been used to measure spectral types
and confirm membership for many of those candidates
\citep{str74b,har54,her98,luh98,luh99,luh03,luh05flam,mue07,alv13}.
Similar diagnostics of youth and membership have been applied
to NGC~1333, including
X-ray emission \citep{pre97,pre03,get02,win10,for11}, outflows \citep{hat09},
optical and near-IR photometry
\citep{str76,asp94,lad96,asp03,wil04,gre07,ots08,sch09},
mid-IR excesses \citep{jor06,jor07,cie07,reb07,gut08,eva09,you15},
kinematics \citep{fos15}, and variability \citep{reb15b}.
As in IC~348, spectroscopic classification has been performed on many
of the resulting candidates
\citep{asp03,wil04,gre07,win09,win10,sch09,sch12a,sch12b}.

The current census of IC~348 is incomplete in the outer portions of the
cluster and among the least massive brown dwarfs.
The census of NGC~1333 has significant incompleteness as well,
particularly at substellar masses.
In addition, the methods of spectral classification that have been applied
to NGC~1333 are less uniform than in IC~348.
To address these issues, we have performed a survey for new
members down to $\sim0.005$~$M_\odot$ across the full extent of each cluster,
and we have measured spectral types for a large fraction of the known members
of NGC~1333 with the classification scheme that we have previously applied
to IC~348 and other nearby star-forming regions.
In our presentation of this work, we begin by compiling lists of all known
members of these clusters from previoius studies (Section~\ref{sec:previous}).
We then select candidate cluster members based on X-ray emission,
mid-IR excess emission, optical and near-IR color-magnitude diagrams, and
proper motions (Section~\ref{sec:select}) and use optical and near-IR
spectra to measure their spectral types and determine whether they are members
(Section~\ref{sec:spectra}). We also measure new spectral types for a
large number of the known members of NGC~1333. We conclude by analyzing
several aspects of the new samples of members of IC~348 and NGC~1333, which
include their completeness, ages, mass functions, disk fractions, and spatial
distributions (Section~\ref{sec:analysis}).

\section{Census from Previous Studies}
\label{sec:previous}

We have searched previous studies of IC~348 and NGC~1333 for objects
that exhibit evidence of membership. For IC~348, we began with the
census of 288 members compiled by \citet{luh03}. In that tabulation, pairs
of objects with separations of $<1\arcsec$ appeared as single entries.
We adopt the same approach in this work.
\citet{luh03} noted that LRL~1434\footnote{When referring to objects
in IC~348, we use the number identifications from our previous studies and
from this work, which are found in the second column of Table~\ref{tab:mem348}.}
was unusually faint for a cluster member
near its spectral type, which would indicate that it is either a field star
or a member that is detected in scattered light, as in the case of an edge-on
disk. Because its H$\alpha$ emission and Na~I absorption in a low
signal-to-noise (S/N) spectrum seemed to indicate that it was young,
\citet{luh03} adopted it as a member. However, it was not detected in
subsequent mid-IR images of the cluster \citep{lad06,mue07}, indicating that
it does not have a disk, and hence is unlikely to be seen in scattered light.
Therefore, we omit LRL~1434 from our sample of members.
In a low S/N spectrum from \citet{luh03}, LRL~624 appeared
to have the weak Na~I absorption that is expected for a young, low-gravity
cluster member. However, we now treat it as a field star based on its proper
motion, which is inconsistent with membership in IC~348 (Section~\ref{sec:pm}).
\citet{luh05wfpc} measured spectral types for a likely companion in the census
from \citet{luh03}, LRL~78~B, and for the candidate companion LRL~166~B.
Because LRL~166~B and its primary are separated by less than $1\arcsec$,
they appear as a single entry in our tabulation of members.
We have added to our sample the 16, 45, and 16 members identified by
\citet{luh05flam}, \citet{mue07}, and \citet{alv13}, respectively.
The latter study presented spectra for several additional objects whose
membership was uncertain. Based on our inspection of those spectra, we
have classified five of those candidates as likely members
(Section~\ref{sec:class}), consisting of sources 3, 5, 14, 20, and 31
from \citet{alv13} (LRL~670, LRL~5209, LRL~10378, LRL~22443, LRL~54229).
We include in our census the protostar HH~211-IR, the candidate protostellar
brown dwarf IC~348-SMM2E \citep{pal14}, and LRL 1898, LRL 54361, LRL 54362,
LRL 54419, LRL 54459, LRL 54460, LRL 55400, and LRL 57025 \citep{mue07}.
The latter eight objects are probable protostars based on their spectral
energy distributions and their proximity to millimeter cores and
known protostars. Through the above steps, we arrived at a sample
of 378 known members of IC~348 based on previous studies.
Later in this work, we describe the identification of 100 new members,
resulting in a total of 478 known members.
We present the full sample of members in Table~\ref{tab:mem348}.
The new members can be identified by the presence of a spectral type from
this study alone with the exception of LRL~60~B and LRL~187~B.
Although the latter stars lacked classifications prior
to our spectroscopy, they are counted as previously known members rather
than new ones since they appeared within the sample of members from
\citet{luh03}.

To construct a census of known members of NGC~1333, we assessed the evidence
of membership for all objects proposed to be members in previous surveys
of the cluster. The evidence consisted of signatures of youth in the form
of strong emission lines, Li absorption, X-ray emission,
mid-IR excess emission, and the shape of the gravity-sensitive steam bands.
We also examined whether the proposed members exhibited radial velocities
\citep{fos15} and proper motions (Section~\ref{sec:pm}) that are
consistent with those of the larger population of objects that show evidence
of youth.
The radial velocity of 2MASS~03290289+3116010 and the proper motion of source
38 from \citet{sch12a} indicate that they are unlikely to be cluster members.
The previously reported members that we find have sufficient evidence of
membership in NGC~1333 are listed in Table~\ref{tab:mem1333}.
We also include in that tabulation the new members found in our study.
This census contains a total of 203 members. For 42 sources, our
new spectral classifications are the only ones available, some of which
were identified as candidate members in previous studies.

In our tabulations of members of IC~348 and NGC~1333, we have included
all known young stars and brown dwarfs within the fields encompassed by the
maps in Figure~\ref{fig:map}. However, one could argue that the
objects at the largest distances from the centers of the clusters
should instead be assigned to the distributed population that
is present across the Perseus cloud \citep{you15}.

\section{Identification of Candidate Members}
\label{sec:select}

To improve the completeness of the census of known members of IC~348
and NGC~1333, we have obtained spectra of candidate members that
have been identified through several signatures of cluster membership.
In this section, we describe the selection of these candidates.

\subsection{Survey Fields}
\label{sec:fields}

We have made use of several imaging surveys of IC~348 and NGC~1333
for identifying candidate members. In Figure~\ref{fig:map}, we show maps
of the positions of the known members of the clusters and the boundaries
of the fields in those surveys.
IC~348 was observed with the imaging array of the Advanced CCD Imaging
Spectrometer (ACIS-I) on the {\it Chandra} X-ray Observatory
\citep{pre01,pre02,for11,ste12}, CFH12K on the
Canada-France-Hawaii Telescope (CFHT) \citep[$IZ$,][]{luh03},
and WIRCam and MegaCam on the CFHT
\citep[$z\arcmin$$JHK_s$,][]{alv13}.
The MegaCam data encompass the entire field surrounding IC~348 in
Figure~\ref{fig:map}.
Data in $ZYJHK$ are also available for all of IC~348 from Data Release 10
of the United Kingdom Infrared Telescope Infrared Deep Sky Survey
\citep[UKIDSS,][]{law07}.
NGC~1333 was observed with ACIS-I \citep{get02,win10,for11} and
Suprime-Cam on Subaru Telescope \citep[$i\arcmin$$z\arcmin$,][]{sch09}.
Unpublished images from WIRCam are also available for NGC~1333.
We reduced all of the WIRCam images in $JHK_s$ that are available
for IC~348 and NGC~1333 from the CFHT archive, which were obtained through
programs 07BH20 (K. Allers), 07BH12 (B. Biller), 09BD95 (L. Albert),
O6BF23, 08BF98, and 09BF50 (J. Bouvier). The data from O6BF23 were those
in IC~348 analyzed by \citet{alv13}. We also reduced the Suprime-Cam images
from \citet{sch09}, which were retrieved from the Subaru Telescope data archive.

Both clusters have been observed on many occasions at mid-IR wavelengths
with the Infrared Array Camera \citep[IRAC;][]{faz04} on the
{\it Spitzer Space Telescope} \citep{wer04}.
Through the c2d {\it Spitzer} Legacy project \citep{eva03}, shallow IRAC images
were obtained for much of the Perseus cloud, including all of
IC~348 and NGC~1333 \citep{jor06,eva09}.
Deeper images were taken for the larger IRAC fields indicated in
Figure~\ref{fig:map} \citep{lad06,mue07,gut08}. The smaller IRAC
fields in Figure~\ref{fig:map} were monitored for approximately one
month \citep{fla13,reb15b}. In the latest IRAC observations,
most of each cluster was imaged at an additional epoch through program
90071 (A. Kraus) to facilitate the identification of candidate members
via proper motions. The IRAC observations prior to May 2009 were performed with
bands at 3.6, 4.5, 5.8, and 8.0~\micron, denoted as [3.6], [4.5], [5.8],
and [8.0], respectively. The later images were collected only in the
[3.6] and [4.5] bands because of the depletion of the liquid helium coolant.
The Multiband Imaging Photometer for {\it Spitzer} \citep[MIPS;][]{rie04}
also has been used to fully map each cluster at 24~\micron\
\citep{lad06,reb07,gut08,cur09}.

In addition to the known members and the survey fields, we also show in
Figure~\ref{fig:map} the positions of candidate disk-bearing stars that were
identified by the c2d survey and that have not been observed with spectroscopy
to confirm membership. Only a few of these candidates are present within
the map of IC~348 and none are found in the map of NGC~1333,
indicating that the distributions of known members trace the full extent of
the clusters. The distribution for IC~348 has a radius of $\sim14\arcmin$,
which is consistent with previous estimates of 10--$15\arcmin$
for the cluster radius \citep{sch99,mue03}.
We wish to achieve a census of the members of these clusters
within fields that have well-defined boundaries, are large enough to
encompass most of the members, and are small enough to be covered by as
many imaging surveys as possible. Given these considerations, we have
searched for new members primarily within a radius of $14\arcmin$ from the B5
star BD+$31\arcdeg$643 in IC~348 and within the $18\arcmin\times18\arcmin$
field in NGC~1333 that was observed by ACIS-I.
We will assess the completeness of our new census within these fields in
Section~\ref{sec:completeness}.

\subsection{X-ray Emission}

Because young stars are bright in X-rays, one can search for
members of star-forming regions via their X-ray emission \citep{fei87,wal88}.
The X-ray studies of IC~348 and NGC~1333 (see Section~\ref{sec:intro})
have done so by checking for X-ray sources that have optical and near-IR
data that are consistent with those expected for cluster members.
Most of the resulting candidate members have been observed with the
spectroscopy that is needed to confirm membership.
We have pursued spectroscopy of the remaining candidates that lack spectra.
For this sample, we selected X-ray sources identified in {\it Chandra}
images of IC~348 and NGC~1333 by \citet{ste12} and K. Getman (in preparation),
respectively, that are not rejected as non-members by our color-magnitude
diagrams (Section~\ref{sec:cmd}). Those two studies have generated
catalogs of sources found in all available ACIS-I images of the two clusters.
In Table~\ref{tab:mem348}, the members of IC~348 that have X-ray detections
can be identified by the presence of source names from \citet{ste12}.
For the members of NGC~1333, we indicate in Table~\ref{tab:mem1333}
the X-ray detections from the new catalog of K. Getman under the column
for membership evidence.

\subsection{Mid-IR Excess Emission}
\label{sec:midir}

When a star is born, it is surrounded by an accretion disk and an infalling
envelope.
The stars in a young cluster that still retain these structures can be
identified via mid-IR emission in excess above that expected from a stellar
photosphere.
The {\it Spitzer} images described in Section~\ref{sec:fields}
have been previously used to search for new members of IC~348 and NGC~1333 in
that manner \citep{mue07,gut08,eva09,you15}.
As with the X-ray candidates, we have sought spectroscopy for the small
fraction of mid-IR candidates that lack previous spectral classifications.
When assembling this sample of candidates, we rejected those
that appear to be knots of extended emission rather than stars based
on visual inspection of the IRAC and MIPS images, which consist of
sources 171, 216, 222, and 227 from \citet{eva09} and \citet{you15}.
If a young star is seen primarily in scattered light, as in
the case of an edge-on disk, it is likely to appear unusually faint for
its color compared to other cluster members. As a result, it would be prone
to rejection as a field star in optical and near-IR color-magnitude diagrams.
Therefore, we have retained mid-IR candidates in our spectroscopic sample
regardless of their locations in the color-magnitude diagrams in
Section~\ref{sec:cmd}.
Active galactic nuclei and stars on the asymptotic giant branch are
common types of contaminants in a sample of this kind that is selected via
red mid-IR colors.
The absence or presence of mid-IR excess emission is indicated for
each known member of IC~348 and NGC~1333 in Tables~\ref{tab:mem348}
and \ref{tab:mem1333}, except for a few of the faintest members that
lack sufficiently accurate mid-IR photometry.

\subsection{Optical and Near-IR Color-Magnitude Diagrams}
\label{sec:cmd}

In the Hertzsprung-Russell (H-R) diagram, the members of a young, coeval
stellar population appear along the main sequence at higher masses and diverge
above the main sequence at lower masses, which is manifested in color-magnitude
diagrams as a band that becomes redder at fainter magnitudes.
For a nearby young cluster, that stellar sequence is brighter than most
foreground and background stars. As a result, color-magnitude diagrams can
be used to select a sample of candidate cluster members that has relatively
little contamination from field stars.

To identify candidate members of IC~348, we have used the diagrams
of $I$ versus $R-I$ and $m_{791}$ versus $m_{791}-m_{850}$
from \citet{luh99} and \citet{luh05wfpc}, respectively.
We also have constructed the following extinction-corrected diagrams:
$K_s$ versus $I-Z$, $I-K_s$, and $Z-K_s$ based on the $I$ and $Z$ data from
\citet{luh03}, $K_s$ versus $Z-Y$, $Z-K_s$, and $Y-K_s$ based on the $Z$ and $Y$
data from UKIDSS, and $K_s$ versus $Z-K_s$ based on the $Z$ data that we have
measured from MegaCam images. The $K_s$ (or $K$) measurements are from the
Two Micron Point Source Catalog \citep[2MASS,][]{skr06}, \citet{mue03},
UKIDSS, and the WIRCam images described in Section~\ref{sec:fields}.
The extinctions of stars in these diagrams were estimated in the manner
described by \citet{luh03tau} using the extinction law from \citet{car89}.
In Figure~\ref{fig:cmd348}, we show the extinction-corrected diagrams
for the known members of IC~348 (including the new ones from this work)
and all other sources with the exception of those that have been
spectroscopically classified as field stars in this work and in previous
studies. In each diagram, we have marked a boundary along the lower envelope
of the locus of known members, which we have used for selecting candidate
members for spectroscopy.
The small number of known members that appear below
the cluster sequence in some of the color-magnitude diagrams are likely seen
in scattered light. They consist of LRL 435, LRL 725, LRL 904, LRL 1287, and
LRL 4011. \citet{luh03} previously noted that LRL 425 and LRL 725 were
unusually faint for their spectral types and colors.

For NGC~1333, we have constructed diagrams of $K_s$ versus $B-K_s$, $R-K_s$, and
$I-K_s$ where $B$, $R$, and $I$ are photographic data from
the USNO-B1.0 Catalog \citep{mon03} and extinction-corrected diagrams
of $K_s$ versus $i\arcmin-z\arcmin$, $i\arcmin-K_s$, and $z\arcmin-K_s$ based
on the $i\arcmin$ and $z\arcmin$ data from Suprime-Cam.
The $K_s$ data are from 2MASS, UKIDSS, and WIRCam.
These diagrams are shown in Figure~\ref{fig:cmd1333} for the known members
of NGC~1333 and all other detected sources, excluding known field stars.
As with IC~348, we have plotted a boundary that follows the lower envelope
of the cluster sequence in each diagram for use in selecting candidate members.
The known members that are below those boundaries consist of
2MASS J03291228+3123065 and sources 39, 66, 92, 102, 110, 113, and 122 from
\citet{gut08}. All of these objects exhibit mid-IR excess emission, so it is
plausible that they are seen in scattered light because of occulting disks.

For each cluster, a source is considered a candidate member
if it is above a boundary in any diagram and is not below a boundary
in any diagram. As an exception to the latter criterion, we do not
reject candidates identified based on mid-IR excesses, as mentioned
in the previous section.

\subsection{Proper Motions}
\label{sec:pm}

Nearby young clusters (150--300~pc) have proper motions of
10--20~mas~yr$^{-1}$, and the members of a given cluster exhibit dispersions
of $\sim1$~mas~yr$^{-1}$ ($\sim1$~km~s$^{-1}$). In comparison, foreground
and background stars typically have much larger and smaller motions,
respectively.
As a result, precise measurements of proper motions can be used to identify
possible cluster members. To measure proper motions in IC~348 and NGC~1333, we
have made use of the multiple epochs of IRAC observations that are available
for the clusters, which were performed through programs 6 (G. Fazio), 36
(G. Fazio), 178 (N. Evans), 30516 (L. Looney), 50596 (G. Rieke), 60160
(J. Muzerolle), 61026 (J. Stauffer), 80174 (K. Flaherty), and 90071 (A. Kraus).
We measured proper motions for all sources in these images by applying the
astrometric techniques and distortion corrections from \cite{esp16}.
Because the different epochs of astrometry have been registered using
the stars that are in common among them (most of which are
background stars), our analysis has produced relative proper motions.
We have ignored proper motion measurements for sources with median values
of S/N in the final epoch of exposures at 4.5~\micron\ that are below
4 and 7.5 and that have errors larger than 6 and 8~mas~yr$^{-1}$
for IC~348 and NGC~1333, respectively. We were able to adopt a lower threshold
of S/N for IC~348 because more epochs of IRAC images are available for it,
which provides lower proper motion errors at a given S/N.
We also inspected the IRAC images of the members that exhibited discrepant
motions compared to the bulk of the population to identify and exclude
measurements that were erroneous due to extended emission or blending of stars.
Two previously identified members, LRL 624 in IC~348
and source 38 from \citet{sch12a} in NGC~1333, appear to have reliable
proper motion measurements that are inconsistent with membership, and hence
have been rejected from the sample of members for each cluster, as mentioned
in Section~\ref{sec:previous}.

Among the 478 and 203 adopted members of IC~348 and NGC~1333, 405 and 141
sources (85 and 69\%) have useful proper motions, respectively.
A larger fraction of the stellar population in IC~348 has measured motions
because of its greater number of IRAC epochs. The proper motions
measured for known members are included in Tables~\ref{tab:mem348} and
\ref{tab:mem1333}. Those data have median values of
($\mu_{\alpha}, \mu_{\delta}=1.9, -2.1$~mas~yr$^{-1}$) for IC~348
and ($\mu_{\alpha}, \mu_{\delta}=2.3, -3.0$~mas~yr$^{-1}$) for NGC~1333.
The motions for the known members and all other sources with measured
motions are plotted in Figure~\ref{fig:pm}.
For each cluster, most of the members are within 3~mas~yr$^{-1}$ of the
median value, and nearly all members are within 6~mas~yr$^{-1}$ of it.
Therefore, we consider objects within 6~mas~yr$^{-1}$ of the median motions
to be candidates, and those within 3~mas~yr$^{-1}$ are the most promising ones.
Sources with motions that differ by $>6$~mas~yr$^{-1}$ are rejected
as non-members (although the previously known members that are only slightly
beyond that threshold are retained as members).
It is evident from Figure~\ref{fig:pm} that the known members overlap with
a large number of field stars in their proper motion measurements.
As a result, the sample of proper motion candidates has significant
contamination from field stars, and it would be inefficient to pursue
spectroscopy of candidates identified based on these measurements alone.
However, these proper motions are valuable for refining the samples of
candidates selected with the other methods that we have already described.
If a candidate found with other diagnostics appeared to be a non-member
based on its motion, we inspected its IRAC images prior to rejection
to check for blends and extended emission that might lead to an erroneous
proper motion measurement.
Many of our spectra of candidates were obtained before we measured the
IRAC proper motions. As a result, some of those objects in our spectroscopic
sample would have been rejected by proper motions if they had been available,
as indicated in Table~\ref{tab:non}.

\section{Spectroscopy}
\label{sec:spectra}

\subsection{Observations}
\label{sec:obs}

We have obtained spectra of 152 candidate members of IC~348 from the analysis
in the previous section and 16 previously known members.
We observed 130 sources with the near-IR spectrograph SpeX
\citep{ray03} at the NASA Infrared Telescope Facility (IRTF).
They consisted of 122 candidates (80 new members, 42 non-members) and eight
known members. Three of the latter objects (LRL~464, LRL~659, LRL~10111)
were candidates at the time of our spectroscopy and were independently
identified as members by \citet{alv13}. One known member in
the SpeX sample (LRL~233) was placed in the slit during the observation
of a candidate at a separation of $3\arcsec$ (LRL~3171).
The spectrum of the candidate was not useful because of low S/N, but
we later successfully observed it with Gemini North. Another known member
observed with SpeX (LRL~62) has a discrepant radial velocity relative to other
cluster members \citep{cot15}, so we sought to verify the evidence of youth
previously found through optical spectroscopy.
The final three members observed with SpeX consist of the companions LRL~60~B
and LRL~187~B, which lacked previous spectral classifications, and LRL~60~A.
The SpeX data were collected in the prism mode with the $0\farcs8$ slit
(0.8--2.5~\micron, $R=150$).
We performed optical spectroscopy on 15 candidates (8 new members,
7 non-members) and six known members (LRL~141, LRL~174, LRL~294, LRL~334,
LRL~366, LRL~10094) with the Inamori Magellan Areal Camera and Spectrograph
(IMACS) on the Magellan I telescope at Las Campanas Observatory. Those data
were taken during the multi-slit observations described by \citet{mue07} for
their sample of IR-selected candidates. The IMACS spectra spanned from
6300--8900~\AA\ and exhibited a resolution of 3~\AA. One of the known members
observed with IMACS (LRL~10094) was originally selected as a candidate member
and was subsequently classified as a member by \citet{alv13}.
We used the near-IR camera on the Keck~I telescope \citep[NIRC,][]{mat94}
to obtain low-resolution ($\sim$100) spectra of one candidate (LRL~6005) and
two known members from \citet{luh05flam} that had uncertain spectral types
(LRL~1050, LRL~2103). They were observed with the gr120 and gr150 grisms, which
together provided coverage from 1--2.5~\micron.
We observed two and 12 candidates with the
Gemini Near-Infrared Imager \citep[NIRI,][]{hod03}
and the Gemini Near-Infrared Spectrograph \citep[GNIRS,][]{eli06},
respectively, at the Gemini North telescope. NIRI was operated with
the $H$-band grism and the $0\farcs75$ slit ($R=500$).
The GNIRS data were collected in the cross-dispersed mode with the
31.7~l~mm$^{-1}$ grating and the $0\farcs67$ slit (1--2.5~\micron, $R=800$).
The 152 candidates observed by these spectrographs consist of 100 new members
and 52 non-members based on the classifications from the next section.

In NGC~1333, we have performed spectroscopy on 55 candidate members and 124
known members with SpeX and 14 candidates with GNIRS. The instrument
configurations were the same as those employed for the targets in IC~348
except for BD+$30\arcdeg$547, for which we used the SXD mode of SpeX
with the $0\farcs8$ slit (0.7--2.5~\micron, $R=750$) because its earlier
type required higher resolution for classification.
We included a large number of known members in our spectroscopic sample
because we wish to maximize the number of members that are classified in
the same manner as the members of IC~348.
In the next section, we find that 42 and 26 candidates are members
and non-members, respectively, and that one of the candidates
(J03284883+3117537) has an uncertain classification because of low S/N.

The SpeX data were reduced with the Spextool package \citep{cus04} and
corrected for telluric absorption in the manner described by \citet{vac03}.
For the IMACS and NIRC spectra, we used routines within IRAF
to apply bias subtraction and flat fielding to the two-dimensional images,
extract spectra of the targets from those images, and perform
wavelength calibration on the extracted data with spectra of arc lamps.
Each NIRC spectrum was also corrected for telluric absorption using the
spectrum of an A star that was observed at a similar airmass.

The dereddened near-IR spectra of the new members of IC~348 and
NGC~1333 are presented in Figures~\ref{fig:sp348a}--\ref{fig:sp348f} and
Figures~\ref{fig:sp1333a}--\ref{fig:sp1333g}, respectively.
We also include all of our SpeX data for previously known members from this
work and previous studies \citep{luh05wfpc,mue07} with the exception of the
four featureless spectra in IC~348 that were presented by \citet{mue07}.
Within the lists of all known members of the clusters in
Tables~\ref{tab:mem348} and \ref{tab:mem1333}, we indicate the dates
and instruments for our new spectra and the previous SpeX data.
The candidates from our spectroscopic sample that we have classified as
non-members are found in Table~\ref{tab:non}.

\subsection{Spectral Classification}
\label{sec:class}

We have used the spectra that we have collected to measure spectral types of
our targets and to help determine whether they are members of IC~348 and
NGC~1333. Seven of our targets of optical spectroscopy lack the M-type
spectral features (TiO, VO)
expected if they were members of IC~348, indicating that they are likely
early-type field stars or giants that are behind the cluster.
The remaining 14 objects with optical spectra do exhibit M-type
features, all of which are young based on their Na~I and K~I
absorption lines, which are sensitive to surface gravity.
We measured spectral types from these data through comparison to the
averages of dwarf and giant standards \citep{luh99}.
In previous studies \citep[e.g.,][]{mue07}, we have measured spectral
types from near-IR spectra of young late-type objects via comparison
to spectra of individual young late-type objects that we had classified at
optical wavelengths. Because of the relatively large number of young
objects that now have both optical types and near-IR spectra, we have
recently combined the spectra of several such objects for each subclass
(K. Luhman, in preparation).  We have used the resulting spectra as
standards for classifying the near-IR spectra in IC~348 and NGC~1333
that exhibit steam absorption ($\gtrsim$M0) and
evidence of youth \citep[e.g., triangular $H$-band continuum][]{luc01}.
For the remaining IR data, spectral types
were measured with spectra of standard dwarfs and giants from our previous
studies and from \citet{cus05} and \citet{ray09}.
In addition, we have revised our previous classifications of SpeX data for
members of IC~348 \citep{luh05wfpc,mue07} using the new standard spectra.
The resulting changes are $\leq0.25$~subclass for most objects.
We also measured spectral types from the near-IR spectra of candidate
members of IC~348 presented by \citet{alv13}.
For one of those objects, LRL~5209, the S/N appeared to be lower than expected
for its magnitude, so we performed our own reduction of the raw data.
We have included the new version of the spectrum with the data
that we have collected in Figure~\ref{fig:sp348f}. Most of our classifications
are similar to those from \citet{alv13}. For those cases, we adopt the
types from \citet{alv13}.
We do present our new spectral types for LRL 670, LRL 5209, and LRL 22443
since they differ noticeably from the previous measurements.
In addition, whereas \citet{alv13} classified the membership status of LRL 670,
LRL~5209, LRL 10378, LRL 22443, and LRL 54229 as uncertain, we find that the
spectra of these objects from \citet{alv13} do show sufficient evidence of
membership in the form of the gravity sensitive features. As a result, we
have included those objects in our census of members, as mentioned in
Section~\ref{sec:previous}.

For the members of IC~348 and NGC~1333 that have measured spectral types
and that have SpeX data from this work and our previous studies, we
have estimated extinctions by comparing the observed spectral slopes at
1~\micron\ to the slopes our young standards and adopting
the extinction law of \citet{car89}.
If a spectral type could not be measured from a SpeX spectrum, but a
classification was available from another source (e.g., optical spectrum),
then we adopted that classification when estimating the extinction.
For a spectrum with low S/N at 1~\micron, we estimated the extinction from the
slope at longer wavelengths if a $K$-band excess was unlikely to be present
based on an absence of a mid-IR excess.
The resulting extinctions were used for dereddening the spectra in
Figures~\ref{fig:sp348a}--\ref{fig:sp1333g}.
Among those dereddened spectra, most of the objects with high extinctions
can be identified by their lower S/N at shorter wavelengths.

The spectral types that we have measured for members of IC~348 and NGC~1333
are included within the lists of all known members in Tables~\ref{tab:mem348}
and \ref{tab:mem1333}.
The errors for the optical and IR types are $\pm0.25$ and
0.5~subclass, respectively, unless indicated otherwise. The uncertainties
in the IR types tend to be large at $\gtrsim$M9 because of a degeneracy between
spectral type and reddening in near-IR spectra with low resolution and
low-to-moderate S/N. For instance, a young M9 object with $A_V\sim3.5$ can
appear quite similar to an unreddened L3 object in data of this kind.
Our classifications for non-members are presented in Table~\ref{tab:non}.
Most of these non-members are giants and early-type stars, which are
difficult to distinguish from K and early-M cluster members in color-magnitude
diagrams. Meanwhile, cooler members have more distinctive colors, so
the yield of confirmed members can be higher at fainter magnitudes when
the appropriate bands of photometry are available. For instance,
13 of the 14 candidate low-mass members of NGC~1333 observed with GNIRS
were confirmed as such.
Finally, we note that several of the candidate members of IC~348 that we
selected from our color-magnitude diagrams and confirmed with spectroscopy
were previously identified as photometric candidates by \citet{alv13},
consisting of LRL 6005, LRL 5231, LRL 1254, LRL 1824, and LRL 22778.

\subsection{Comparison to Previous Work}
\label{sec:compare}

We have compared our membership lists and spectral types for IC~348
and NGC~1333 to those from previous studies.
We have presented spectra for 100 and 42 members of these clusters,
respectively, that have not been previously classified, which we refer
to as new members. Two additional stars, LRL~60~B and LRL~187~B, also
lack previous classifications, but they are treated as previously
known members, as mentioned in Section~\ref{sec:previous}.
To illustrate the magnitudes and spectral types at which the census of
these clusters has expanded, we show in Figure~\ref{fig:histonew} the
distributions of extinction-corrected $M_K$ and spectral types for
previously known and new members. We have adopted the extinctions
estimated from our IR spectra in Section~\ref{sec:class} when available.
Otherwise, extinctions are estimated with photometry in the
manner described in Section~\ref{sec:hr}.
Members that lack measured spectral types
(and hence extinctions) are absent from Figure~\ref{fig:histonew}, which
consist of protostars with featureless spectra.
The new members of IC~348 are predominantly low-mass stars at M4--M6 in
the outskirts of the cluster, but they also include several objects that are
the faintest known members. In NGC~1333, the new members are
distributed more uniformly with magnitude and spectral type among the
low-mass stars and brown dwarfs.

\citet{reb15b} compiled a sample of 130 members of NGC~1333, 17 of which
are absent from our census of the cluster for the following reasons.
Seven of these missing objects are among the candidate members that we
have identified in Section~\ref{sec:select}, but they lack spectroscopic
confirmation of their membership. They consist of 2MASS J03291532+3129346
and NGC 1333 IRS J03284883+3117537, J03285358+3112147, J03285508+3114163,
J03291565+311911, J03285709+3121250, and J03291317+3119495.
Sources 12, 45, 113, and 174 from \citet{win10} have been detected in X-rays,
but are rejected as field stars by our color-magnitude diagrams and have
no other evidence of membership. Using the identifiers from \citet{gut08},
sources 38 and 76 are field stars based on our spectroscopy,
sources 8 and 38 are an outflow lobe and a galaxy, respectively
\citep{arn12}, source 95 has uncertain membership (Section~\ref{sec:comments}),
and source 26 is extended at both near- and mid-IR wavelengths, so it is
unclear whether it a star.

\citet{reb15b} identified five new candidate members of NGC~1333 based on
their mid-IR variability.
One of these candidates, SSTYSV J032903.46+311617.9, is in our spectroscopic
sample. Based on its featureless near-IR spectrum and red color in
$[3.6]-[4.5]$ (it is blended with a brighter star at longer wavelengths),
it is probably a protostar, so we have included it in our census of members.
\citet{reb15b} discussed at length a second of their candidates,
SSTYSV J032911.86+312155.7. They classified it as a possible protostar based
on its flux at 24~\micron\ relative to shorter wavelengths. However, we find
that no detection is apparent in the 24~\micron\ images.
Bright extended emission associated with the nearby star SVS~3 prevents
detections in the {\it Spitzer} bands longward of 4.5~\micron.
Nevertheless, given its close proximity to other members like SVS~3 and
the excess at 4.5~\micron\ relative to bands at shorter wavelengths,
it is a promising candidate member.
We include it in our list of candidates that have not been observed with
spectroscopy in Section~\ref{sec:completeness}.
Two other objects from \citet{reb15b}, SSTYSV J032918.65+312021.8 and
SSTYSV J032907.24+312409.7, are also in our sample of candidates based
on their variability, mid-IR excess emission, and their positions
in our color-magnitude diagrams.
The final candidate from \citet{reb15b}, SSTYSV J032836.43+312856.7,
does not exhibit mid-IR excess emission and is rejected by our color-magnitude
diagrams, so we do not consider it to be a candidate member.

As mentioned in Section~\ref{sec:obs}, we have performed spectroscopy on
a large number of the previously known members of NGC~1333 with the goal
of obtaining spectral classifications that are derived with the same
scheme applied to IC~348 in our previous studies and this work.
Spectral types are available from both our work and previous
studies for 90 objects. For a majority of the previous classifications,
there is not a systematic difference from our measurements.
However, we do find that our spectral types are an average
of $\sim1$~subclass earlier than those from \citet{sch09,sch12a,sch12b},
which corresponds to $\sim50$\% higher
mass estimates when combined with evolutionary models \citep[e.g.,][]{bar98}.
There remain several (very faint) members from \citet{sch09,sch12a,sch12b}
for which we have not obtained spectra and measured spectral types.

\subsection{Comments on Individual Sources}
\label{sec:comments}

{\it LRL~62 and LRL~155}. The optical spectrum of LRL~62 from \citet{luh99}
exhibits a spectral type of M4.5 and evidence of youth in the form of weak K~I
and Na~I lines. However, its radial velocity of $\sim$5~km~s$^{-1}$ differs
significantly from the mean value of 15.4~km~s$^{-1}$
($\sigma=0.7$~km~s$^{-1}$) for a sample of known members \citep{cot15}.
We obtained a near-IR spectrum to verify its youth through additional
gravity-sensitive features. Our new classification based on that spectrum
is consistent with the optical result. Like LRL~62,
LRL~155 also has a discrepant radial velocity ($\sim$23~km~s$^{-1}$).
Its near-IR spectrum indicates a spectral type of M1, but it is not sensitive
to signatures of low surface gravity for this type.
Both stars have higher extinctions ($A_J=0.7$ and 1) than
expected for foreground dwarfs ($A_J<0.1$) and are too bright for background
dwarfs. Therefore, we treat them as members of IC~348.
The stars may have been ejected from the cluster through dynamical interactions
with other members \citep{kro98,wei11}.

{\it IC 348 IRS J03442484+3213482}.
It was identified as a candidate protostar by \citet{eva09} and \citet{you15}
(source 405 in those studies). Its near-IR spectrum from SpeX is featureless.
Protostars can have featureless spectra if significant veiling is present.
However, such spectra are normally very red (see Fig.~\ref{fig:sp1333a}),
whereas the spectrum of this object is much bluer ($H-K_s\sim0.5$).
In addition, it is far from the area in the southern part
of the cluster where most of the known protostars are found.
As a result, we conclude that it is more likely to be a galaxy than
a protostar, and we classify it as a non-member for the purposes of this work.

{\it BD+$30\arcdeg$547}. \citet{pre97} described it as a likely foreground star,
but its proper motion is consistent with membership in NGC~1333
(E. Mamajek, private communication) and it exhibits X-ray emission.
Some previous studies have concluded that it has IR excess emission at
24~\micron\ \citep{eva09,you15,reb15b}, but that appears to be due to
contamination from a protostar at a separation of $3\farcs4$.

{\it 2MASS J03290575+3116396}. Our SpeX type (M0--M2) is much later than
previous types (late-G, A3).

{\it [SVS76] NGC 1333 7}. It has been classified as an early B giant
through optical spectroscopy \citep{tur80}. Meanwhile, several studies have
reported IR excess emission for this star based on images
from {\it Spitzer}. However, it is detected with only low S/N at 8~\micron\ and
is not detected at 24~\micron\ because of bright extended emission.
Among the bands at shorter wavelengths where better photometry is available,
only the [5.8] band exhibits a significant excess.
Given the excess in that band and its location near the center of the cluster,
the star is a promising candidate member, but we exclude it from
our census of confirmed members because of the uncertainty in its
spectral classification.

{\it NGC 1333 IRS J03290347+3116179}. It has a red, featureless near-IR
spectrum. It is too close to a brighter star (2MASS J03290375+3116039)
to be detected by {\it Spitzer} at $>5$~\micron, but its
$[3.6]-[4.5]$ color is indicative of a protostar \citep{reb15b}, which
would be consistent with the appearance of its spectrum.
Its membership is further supported by its detection in X-rays and its close
proximity to known members of the cluster.

{\it 2MASS J03290895+3122562}. Its near-IR spectrum is red and featureless,
which is consistent with the protostellar nature implied by its
IRAC colors \citep{gut08,eva09}.

{\it 2MASS J03294415+3119478}. The strength of the steam bands for this
object imply a type of M7--M8, but the overall slope of the SpeX data
does not agree with that of any reddened standard. In addition,
it is fainter than most members near its type in color-magnitude diagrams.
With these characteristics, it is similar to 2MASS~J04381486+2611399, which
is a low-mass member of Taurus that is seen in scattered light from an edge-on
disk \citep{luh07edge}.

{\it 2MASS J03283695+3123121}.
\citet{pre97} suggested that it is a foreground star based on its proper
motion from \citet{her83}, but the motion that we have measured with
IRAC is consistent with membership in NGC~1333.

{\it NGC 1333 IRS J03284883+3117537}.
The S/N of our spectrum of this object is too low for classification.
It is a candidate member based on its location on color-magnitude diagrams,
proper motion, and excesses in the IRAC bands. It also may
be detected at low S/N at 24~\micron\ in images from MIPS, which would
further support the presence of excess emission. It is included in
our sample of remaining candidate members that lack classifications
(Section~\ref{sec:completeness}).

{\it 2MASS 03302246+3132403}. \citet{eva09} and \citet{you15} identified it
as a possible protostar based on its mid-IR excess emission.
However, it is detected in the optical bands of the Digitized Sky Survey,
whereas most protostars are too heavily reddened for detections in those data.
In addition, it is not near known protostars or high column densities of gas.
Therefore, we conclude that it is probably a galaxy.

\section{Analysis of New Census}
\label{sec:analysis}

\subsection{Completeness}
\label{sec:completeness}

As discussed in Section~\ref{sec:fields}, we have focused our survey
for members of IC~348 and NGC~1333 within a radius of $14\arcmin$ from
BD+$31\arcdeg$643 in the former and within the ACIS-I images of the latter,
which cover a field with a size of $18\arcmin\times18\arcmin$
(see Fig.~\ref{fig:map}).
To characterize the completeness of our new census for each field,
we use a color-magnitude diagram in two bands that can
detect objects at both low masses and high extinctions.
Given the available data, the best options for these bands are $H$ and $K_s$.
In Figure~\ref{fig:hk}, we plot diagrams of $K_s$ versus $H-K_s$ for
all known members of IC~348 and NGC~1333, the remaining candidate members
identified in Section~\ref{sec:select} that lack spectra and are
within the survey fields, and all other objects within those fields that
are detected in $H$ and $K_s$ and that are not rejected as field stars
by any of the color-magnitude diagrams that we used in selecting candidates.
Thirteen and 20 members of IC~348 and NGC~1333, respectively, are absent from
those diagrams, which consist of companions that are unresolved from brighter
stars and protostars that are extended\footnote{Some of these protostars have
measurements in the 2MASS Point Source Catalog, but are found to be
dominated by extended emission in the higher resolution images from WIRCam.}
or are not detected in $H$ or $K_s$. In Figure~\ref{fig:hk}, there are few
remaining objects in the survey fields with undetermined membership status down
to rather faint magnitudes for low-to-moderate levels of extinction.
Specifically, our census appears to be nearly complete for extinction-corrected
magnitudes of $K_s<16.8$, 15.8, and 15.3 in the IC~348 field and for
$K_s<17.3$, 16.2, and 15.3 in the NGC~1333 field for $A_J<1.5$, 3, and 5,
respectively.

We can use our new census of IC~348 and NGC~1333 to examine the completeness
of previous surveys for members. \citet{luh03} estimated that their census
of a $16\arcmin\times14\arcmin$ field in IC~348 was nearly complete for
extinction-corrected magnitudes of $H<15.5$ ($\lesssim$M8) for $A_V<4$
($A_J<1.13$). Subsequent studies have not uncovered any additional members
in their survey field and in that range of magnitudes and extinctions.
\citet{alv13} noted that three of their new members were within the field
from \citet{luh03}. However, one of those objects, LRL~659, actually
falls slightly outside of that field. The other two members, LRL~2050 and
22528 (M8 and M9), are fainter than the completeness limit from \citet{luh03}.
Meanwhile, \citet{alv13} concluded that their survey for brown dwarfs
($\gtrsim$M6.5) was complete down to $\sim0.013$~$M_\odot$ ($\lesssim$M9)
for $A_V\leq4$ for a field that encompasses the entire cluster.
However, we have found several new members at M6.5--M9 with $A_V<4$,
consisting of LRL~1254, LRL~1824, LRL~5103, LRL~10256, LRL~22185, LRL~22191,
LRL~22317, and LRL~22778.
\citet{sch09,sch12a,sch12b} searched for members of NGC~1333 within a
$30\arcmin\times30\arcmin$ field covering the entire cluster using images
that exhibited completeness limits of $J=20.8$ and $K=18$. We have
identified 32 additional members above those limits (13 at $\leq$M6,
19 at $>$M6). The incompleteness in their census
and the systematic offset between their spectral classification scheme
and that applied to other regions like IC~348 (Section~\ref{sec:class})
cast doubt on the validity of the statements by
\citet{sch09,sch12a,sch12b,sch13} regarding the mass function in NGC~1333.

We also can characterize the fraction of members in our census that have
been detected in X-rays by ACIS-I on {\it Chandra}.
\citet{ste12} analyzed the four existing ACIS-I observations of IC~348,
arriving at list of 290 detected sources. They found X-ray counterparts
for 187 of the 316 members from their adopted cluster census that were
within the ACIS-I images. Using our updated census, 388 known members
were observed by ACIS-I, 197 of which were detected.
The members with X-ray detections can be identified in Table~\ref{tab:mem348}
via the presence of source designations from \citet{ste12}.
K. Getman (in preparation) has performed a similar analysis for the
existing ACIS-I data in NGC~1333.
Among the 186 known members within those images, 98 have counterparts
in their catalog of ACIS-I sources, as indicated in the column for evidence
of membership in Table~\ref{tab:mem1333}.
In Figure~\ref{fig:histox}, we plot the distributions of extinction-corrected
$M_K$ and spectral types for all known members of IC~348 and
NGC~1333 within the ACIS-I images and for the members detected in those data.
As expected, the fraction of members with X-ray detections decreases with
fainter magnitudes and later spectral types, quickly approaching zero
at $M_K>6$ and $>$M7.

In Tables~\ref{tab:cand348} and \ref{tab:cand1333}, we present our remaining
candidate members that lack spectra and that are within the $14\arcmin$ radius
field in IC~348 and within the ACIS-I field in NGC~1333. Although it has been
previously observed with spectroscopy, [SVS76] NGC 1333 7 is included with
these candidates since its membership is uncertain (Section~\ref{sec:comments}).
The probability of membership varies substantially among these candidates.
Those identified via both color-magnitude diagrams and proper motions are
promising while those selected by proper motions alone (i.e., they lack
the optical data needed for the color-magnitude diagrams) are much less likely
to be members (Section~\ref{sec:pm}).
Based on their positions in the color-magnitude diagrams in
Figure~\ref{fig:hk}, most of the candidates should have spectral types
of $\gtrsim$M6 if they are members. To check whether they have the near-IR
colors expected for those types, we have included diagrams of $J-H$ versus
$H-K_s$ for each cluster in Figure~\ref{fig:hk}.
Some of the candidates at $H-K_s<1.5$ do resemble known late-type members
in both $J-H$ and $H-K_s$, but most of the candidates at $H-K_s>1.5$ have
colors indicative of earlier types, and thus are likely to be background stars.

\subsection{H-R Diagrams}
\label{sec:hr}

We have constructed H-R diagrams for the known members of IC~348 and NGC~1333
in terms of $M_K$ versus spectral type. We use absolute magnitude
and spectral type instead of bolometric luminosity and effective temperature
to avoid uncertainties in bolometric corrections and conversions between
spectral type and temperature.
We choose $K_s$ for the band of the absolute magnitude because it is long
enough in wavelength that extinctions are relatively low for most objects
while short enough in wavelength that the fluxes are likely to be dominated
by stellar photospheres rather than circumstellar disks. In addition,
given the sensitivities of the available images, those at $K_s$ detect the
largest fraction of the clusters members. The analysis in this section
was also performed with the $J$ and $H$ bands, which produced identical
results to those from $K_s$.
We estimated extinctions for all known members of IC~348 and NGC~1333 that
have measured spectral types in a similar manner as done by \citet{fur11}
for members of Taurus. For the members that we observed with IR spectroscopy,
we have adopted the extinctions derived during the spectral classifications.
For each of the remaining objects, we calculated the extinction from the
excess in $J-H$ relative to the color expected for a young stellar
photosphere at the spectral type in question \citep{luh10tau}.
After correcting the $K_s$ measurements for extinction, we converted
them to absolute magnitudes using distances of 300~pc for IC~348
\citep{her08} and 235~pc for NGC~1333 \citep{hir08}.
The resulting values of $M_K$ are plotted as a function of spectral type in
Figure~\ref{fig:hr}. For comparison, we have included data for the members
of the Upper Sco association compiled by \citet{luh12usco}. We have
adopted a distance of 145~pc for Upper Sco \citep{pm08} and have estimated
extinctions with the same methods that were applied to IC~348 and NGC~1333.

Several members of IC~348 and NGC~1333 are unusually faint for their
spectral types, appearing below the cluster sequences in Figure~\ref{fig:hr}.
These objects include LRL~276, LRL~435, LRL~621, LRL~622, LRL~725, and
LRL~4011 in IC~348 and sources 39, 92, 99, 110, and 122 from \citet{gut08}
in NGC~1333. The first five sources in IC~348 exhibited similar positions
in the H-R diagram from \citet{luh03}. As noted in Section~\ref{sec:midir},
stars that are occulted by circmumstellar disks often are observed
primarily in scattered light, which results in underestimates of their
luminosities. All of these stars do show evidence of disks in the form
of mid-IR excess emission, so it is plausible that they are occulted by
edge-on disks.

To compare the ages of IC~348, NGC~1333, and Upper Sco, we have computed the
median values of $M_K$ as a function of spectral type for each population.
This was done by applying local linear quantile regression with the
function {\tt lprq} in the {\it quantreg}
package \citep{koe16}
within R \citep{R} using a bandpass of one spectral type.
The resulting median sequences are plotted together in Figure~\ref{fig:hr}.
We do not include the median for $<$M1 in NGC~1333 because of the small
number of members at those types.
The sequences for IC~348 and NGC~1333 are not offset vertically
from each other, which would suggest that they have similar ages.
However, NGC~1333 exhibits clear evidence of a younger age in the form of
a greater abundance of protostars and
circumstellar disks \citep[][Section~\ref{sec:disks}]{mue07,gut08} and
higher extinction (Figure~\ref{fig:hk}).
In order for the H-R diagram to produce a younger age for NGC~1333, we would
need to adopt a larger distance for it \citep{her83} or a smaller distance
for IC~348 \citep{rip14}, e.g., if the clusters have similar distances.
Indeed, one would not expect the cluster distances to differ as much
as we have assumed (65~pc) given that their projected separation is only
$\sim$17~pc, although it is possible that they reside in separate clouds
along the line of site rather than a single cloud \citep{bal08}.
Meanwhile, the median sequences for IC~348 (at 300~pc) and NGC~1333 (at 235~pc)
are 0.4~mag brighter than the sequence for Upper Sco (at 145~pc),
which corresponds to an age difference of 0.25~dex based on evolutionary
models \citep[e.g.,][]{bar98,bar15}.
If Upper Sco has an age of 11~Myr \citep{pec12}\footnote{A younger
age of 4--5~Myr has also been proposed for Upper Sco based on its low-mass
stars \citep{deg89,pre02b,sle06,her15}.}, then IC~348 and NGC~1333 would have
ages of 6~Myr. The latter agrees with the value derived for
IC~348 by \citet{bel13} from color-magnitude diagrams and evolutionary models.
However, a distance of 250~pc was adopted for IC~348
in that study. Using that distance, the sequences in Figure~\ref{fig:hr}
would indicate similar ages for IC~348 and Upper Sco, which would be difficult
to reconcile with the fact that IC~348 has higher abundances of disks and
protostars and, unlike Upper Sco, is still associated with a molecular cloud.
The {\it Gaia} mission \citep{per01} should soon help isolate the sources
of these discrepancies by providing accurate parallactic distances for IC~348,
NGC~1333, and Upper Sco, as well as other nearby clusters and associations.

\subsection{Initial Mass Functions}

Previous studies have estimated the initial mass functions (IMFs) in
IC~348 and NGC~1333 based on earlier samples of spectroscopically confirmed
members \citep{luh98,luh03,gre07,alv13,sch09,sch12a,sch12b}\footnote{The IMFs
in these clusters also have been constrained via IR luminosity functions
\citep{lad95,lad96,mue03}}.
Typically, the masses of individual objects were derived by combining estimates
of bolometric luminosities and effective temperatures with the values predicted
by evolutionary models.
As a result, the IMFs depended on the adopted bolometric corrections,
temperature scales, and models.
To avoid those dependencies, we examine the IMFs in IC~348 and NGC~1333
in terms of observational parameters that should be roughly correlated with
stellar mass, spectral type and extinction-corrected $K_s$. Note that these
parameters still depend on the methods adopted for measuring spectral
types and extinctions.
As done in our previous studies of IMFs in star-forming regions, we
attempt to construct a sample of members in each cluster that is representative
and unbiased in terms of mass by considering all known members within
a field and an extinction threshold for which the current census has
a high level of completeness.
Guided by the analysis of completeness in Section~\ref{sec:completeness},
we select extinction thresholds that are high enough to encompass large
numbers of members while low enough that that the completeness extends to low
masses, arriving at $A_J<1.5$ for the $14\arcmin$ radius field in IC~348
and $A_J<3$ for the ACIS-I field in NGC~1333.
In Section~\ref{sec:completeness}, we found that the census of IC~348
and NGC~1333 within these fields and extinction limits should be
nearly complete for extinction-corrected magnitudes of $K_s<16.8$ and 16.2,
respectively. These samples contain 341 and 120 members, respectively,
which correspond to 71\% and 59\% of the known members.

The distributions of spectral types and extinction-corrected $M_K$
for our extinction-limited samples of members of IC~348 and NGC~1333 are
plotted in Figure~\ref{fig:histoav}.
Relative to IC~348, NGC~1333 exhibits a surplus of objects with late
spectral types and faint magnitudes.
For instance, N($\geq$M6.5)/N($<$M6.5)=54/287=0.188$^{+0.025}_{-0.02}$
and 42/78=0.538$\pm0.056$ and
N($M_K\geq6.5$)/N($M_K<6.5$)=50/291=0.172$^{+0.025}_{-0.019}$
and 40/80=0.50$^{+0.056}_{-0.054}$ for IC~348 and NGC~1333, respectively.
There are multiple possible explanations for the differences in these ratios.
They could reflect a variation in the IMF between the two clusters, although
a significant variation would be surprising given that the clusters
have similar stellar densities and environments, and indeed have arisen
from the same cloud (or at least related clouds).
Because most of the spectral types in NGC~1333 have been measured with
IR spectra, and the resulting classifications tend to have larger uncertainties
than the optical types that are more frequently available in IC~348, one
would expect the distribution of spectral types to be somewhat broader
in NGC~1333, which would result in a higher abundance of late spectral types.
However, this effect would not explain the difference
between the two clusters in their distributions of $M_K$.
Another possibility is that the average extinctions of members of star-forming
regions vary with stellar mass, in which case extinction-limited samples
would not be representative of the stellar populations in those clusters.
For instance, if members at lower masses tend to have less extinction
in the most embedded clusters like NGC~1333, then an extinction-limited sample
could capture most of the brown dwarfs but miss many of the stars.
In a cluster like IC~348 that has dispersed more of its natal cloud,
the dependence of extinction on stellar mass could be smaller, leading to
a larger, more representative ratio of stars to brown dwarfs in an
extinction-limited sample.
This scenario could account for the differences in the
extinction-limited samples for IC~348 and NGC~1333.
Finally, we note that both samples contain several known members that
are fainter than the completeness limits, and the degree of incompleteness
beyond those limits may differ between the two samples, which could
somewhat inflate the perceived abundance of low-mass objects in one cluster
relative to the other.
To determine whether these last two issues are
responsible for the surplus of low-mass objects in the extinction-limited
sample for NGC~1333 relative to the sample in IC~348, it will be necessary
to obtain additional data (e.g., photometry, spectroscopy, proper motions)
that can extend the completeness limits of the census to higher extinctions
and lower masses in both clusters.

Our work has provided new constraints on the minimum masses of the IMFs
in IC~348 and NGC~1333. In each cluster, members are present down to
and below the completeness limits, and thus the minimum of the IMF has
not been detected. The faintest known members have $M_K=10.4$ and 11.2,
which correspond to masses of $\sim$0.004--0.006 and 0.003--0.005~$M_\odot$,
respectively, for ages of 1--3~Myr according to the evolutionary models of
\citet{bur97} and \citet{cha00}.

\subsection{Disk Fractions}
\label{sec:disks}

We can combine our census of IC~348 and NGC~1333 with the previous mid-IR
imaging of these clusters to measure the fractions of members that have
circumstellar disks.
A disk is present in the first three stages of a young stellar object, which
consist of classes~0 and I (protostar+disk+infalling envelope)
and class~II \citep[star+disk,][]{lw84,lad87,and93,gre94}.
A star that has fully cleared its primordial disk is in the class~III stage.
A disk fraction can be defined as either N(I+II)/N(II+III) or N(II)/N(II+III)
(class~0 objects are rare enough that their contribution is usually negligible).
Because the ages measured for young clusters (with H-R diagrams)
often apply to the class~II
and III sources, we choose the latter definition, as done in \citet{luh10tau}.
Since protostars normally have heavily veiled, featureless spectra
(Fig.~\ref{fig:sp1333a}), we can exclude them from our calculations of disk
fractions by considering only members that have measured spectral types.
Members of star-forming regions are often discovered based on mid-IR excess
emission from disks. As a result, samples of members can be biased in favor
of disks, making it difficult to measure disk fractions that are representative
of the stellar populations. However, IC~348 and NGC~1333 have been thoroughly
surveyed for members using a variety of methods, and we have shown that
the current membership samples for the $14\arcmin$ radius field in IC~348
and the ACIS-I field in NGC~1333 are nearly complete for a wide range
of masses and extinctions (Section~\ref{sec:completeness}).
We have assigned the presence or absence of mid-IR excess emission for each
member based on the results of previous disk surveys in these clusters
with {\it Spitzer}
\citep{luh05frac,lad06,mue07,gut08,cur09,eva09,arn12,reb15b,you15}.
A few of the faintest brown dwarfs lack sufficiently accurate photometry for
determining whether excess emission in present; they are excluded from our
calculations of disk fractions.
In Table~\ref{tab:disks} and Figure~\ref{fig:disks}, we list and plot
the fraction of members that have excess emission as a function of spectral
type for the $14\arcmin$ radius field in IC~348 and for the ACIS-I field in
NGC~1333.
The latter has a higher disk fraction, which agrees with previous
analysis of near-IR photometry \citep{lad95,lad96} and the {\it Spitzer} data
\citep{lad06,mue07,gut08}.
The higher disk fraction in NGC~1333 is consistent with its higher
abundance of protostars \citep{mue07,gut08} and the younger age implied by
its greater obscuration. In each cluster, the disk fraction is
roughly constant for the full range of spectral types.

\subsection{Spatial Distributions}

The spatial distributions of the stellar populations in IC~348 and NGC~1333
has been previously studied through analysis of probable members detected
in near- and mid-IR imaging \citep{lad95,lad96,mue03,gut08}.
The current census of each cluster now offers confirmation of membership,
measurements of spectral types for most members, and a high level of
completeness for most locations, masses, and extinctions.
These features allow us to examine the spatial distributions of spectral types
and offsets in $M_K$ from the median cluster sequence at a given spectral
type ($\Delta M_K$), which
serve as proxies for stellar masses and ages, respectively. We can also
measure the spatial dependence of disk fractions using the mid-IR excess data
compiled in the previous section. We define $\Delta M_K$ as $M_K$(median
sequence at a star's spectral type)$-M_K$(star), i.e., higher
values of $\Delta M_K$ correspond to younger implied ages.
We consider only members within the $14\arcmin$ radius field in IC~348
and the ACIS-I field in NGC~1333 because of the well-defined completeness
in those areas. The class~0 and I objects exhibit distinct spatial distributions
compared to members in the more evolved classes, so we exclude them from our
analysis by considering only members that have measured spectral types.
We also omit members that
are later than M9 since their types tend to have large uncertainties.
For each of these samples of members, we have computed surface density
as a function of position using the {\tt kde2d} function in the R package
{\it MASS} \citep{ven02},
which performs a two dimensional kernel density estimation with a bivariate
normal kernel. We then identified the density contours that would divide the
sample for a given cluster into three subsets that have equal numbers.
We selected that number of sections to allow coarse measurements
of the variations of median spectral type, $\Delta M_K$, and disk fraction
with surface density while also providing good number statistics within
each section. The resulting contours are plotted in
Figure~\ref{fig:map2} with the locations of all known members of each cluster.
The average surface densities in these sections are 0.3, 1.0, and
3.0~arcmin$^{-2}$ in IC~348 and 0.2, 0.9, and 2.6~arcmin$^{-2}$ in NGC~1333.
As found in previous studies, IC~348 is more centrally concentrated than
NGC~1333, which exhibits a double cluster morphology \citep{lad96}.

For each of the three sections within IC~348 and NGC~1333 in
Figure~\ref{fig:map2}, we have computed the median spectral type,
the median $\Delta M_K$, and the disk fraction.
The errors in the medians were estimated with bootstrapping.
The resulting values are plotted for each section in Figure~\ref{fig:map3}.
None of these parameters exhibit significant variations among the
sections in either cluster.
Through analysis of the near-IR luminosity function of IC~348, \citet{mue03}
found a higher abundance of solar-mass stars in the core relative to the
outskirts of the cluster.
That surplus is also detected when distributions of spectral types are
compared between the inner and outer portions of the cluster,
although it is does not have a noticeable effect on the median types
because low-mass stars are the dominant component of the stellar population
throughout the cluster.
Using X-ray and near-IR photometry, \citet{get14} have detected
age gradients in NGC~2024 and the Orion Nebula Cluster in which younger
stars are found in the cores of the clusters.
A trend of that kind is not present in the median values of $\Delta M_K$
for IC~348 and NGC~1333. For perspective, an age gradient like that
reported for Orion (1.2--1.9~Myr) should correspond to a difference of
0.3~mag in $M_K$ according to evolutionary models of low-mass stars.
The errors in median $\Delta M_K$ are larger for NGC~1333
than for IC~348 because the former has a broader sequence at a given
spectral type.

We note that a variety of more sophisticated methods are available for
characterizing the spatial distributions of members of star-forming clusters
\citep{gut09,kuh14} and searching for evidence of mass segregation
\citep{sag88,hil98,all09,mas11}.
The optimum approach for measuring the latter has been a subject of
debate in recent years \citep{asc09,olc11,par15}.

\subsection{Candidate Binary Systems}

Our census of IC~348 and NGC~1333 may contain resolved components of
multiple systems. In Table~\ref{tab:pairs}, we have compiled
all pairs of objects from our census that have separations less than
$6\arcsec$. We have omitted binaries that have been resolved only
in high-resolution imaging and that lack spectral classifications of both
components \citep{duc99}. For some of these pairs, both components have
spectral types of late M, making them candidates for wide binary brown
dwarfs \citep{luh04bin,luh09fu}. The numbers of pairs are 23 and 8
for IC~348 and NGC~1333, respectively. To roughly estimate the fraction
of these pairs that comprise binary systems, we performed
a Monte Carlo simulation of the projected separations of unrelated cluster
members using surface density maps of the known members.
In $\sim$90\% of the realizations, the number of chance alignments with
separations of $<6\arcsec$ is between 5--15 for IC~348 and between 2--8
for NGC~1333. Thus, a significant fraction of the candidate binaries
could consist of unrelated cluster members.

\section{Conclusions}

We have sought to improve the completeness of the census of stars and
brown dwarfs in IC~348 and NGC~1333
and the accuracies of spectral types of known members of the latter.
The results of this study are summarized as follows:

\begin{enumerate}

\item
We have obtained optical and near-IR spectra of candidate members of
IC~348 and NGC~1333 that have been selected based on X-ray emission, mid-IR
excess emission, positions in color-magnitude diagrams, and proper motions.
We have classified 100 and 42 of the candidates as new members of IC~348
and NGC~1333, respectively.
The total numbers of known members are now 478 and 203.
Two stars in IC~348, LRL~62 and LRL~155, have radial velocities that
differ significantly from that of the bulk of the cluster,
but they exhibit other evidence of membership, which suggests that
they may have been ejected via dynamical interactions.

\item
We have searched for new members primarily within a radius of $14\arcmin$ from
the B5 star BD+$31\arcdeg$643 in IC~348 and within the area in NGC~1333 that
was observed by ACIS-I on {\it Chandra} ($18\arcmin\times18\arcmin$).
These fields are large enough to encompass most or all members of the clusters.
The new census is nearly complete for extinction-corrected magnitudes of
$K_s<16.8$, 15.8, and 15.3 in the IC~348 field and for
$K_s<17.3$, 16.2, and 15.3 in the NGC~1333 field
for $A_J<1.5$, 3, and 5, respectively. For perspective, $K_s=15$ and 17
correspond to masses of $\sim$0.025 and 0.008~$M_\odot$, respectively, for
an age of 3~Myr according to evolutionary models \citep{bur97,cha00,bar15}.
IC~348 and NGC~1333 now have two of the most complete membership lists
among star-forming clusters.

\item
The known members of IC~348 and NGC~1333 extend down to (and below)
the completeness limits of the current census. As a result, we have not
yet detected the low-mass cutoffs in the mass functions of these clusters.
The faintest known members have $M_K=10.4$ and 11.2, which imply masses of
$\sim$0.004--0.006 and 0.003--0.005~$M_{\odot}$, respectively, for ages of
1--3~Myr based on evolutionary models.

\item
In addition to the candidate members, we have performed spectroscopy
on a large fraction (77\%) of the previously known members of NGC~1333.
These data provide greater uniformity in the spectral types among members
of this cluster and relative to stars in other young clusters like IC~348.

\item
To estimate the IMFs in IC~348 and NGC~1333, we have attempted to select
a sample of members in each cluster that is unbiased in terms of mass.
We have constructed extinction-limited samples for this purpose,
which should have a high level of completeness down to low masses.
The resulting sample for NGC~1333 has a higher abundance of low-mass objects
than the sample for IC~348.
For instance, N($\geq$M6.5)/N($<$M6.5)=0.188$^{+0.025}_{-0.02}$ in IC~348
and 0.538$\pm0.056$ in NGC~1333. Similar fractions are found when the clusters
are compared in terms of $M_K$.  A variation in the IMF between the clusters
would be surprising given their similar densities and environments.
Instead, it is possible that average extinctions are lower for
objects at lower masses, in which case extinction-limited samples may be
biased in favor of low-mass objects in heavily embedded clusters like NGC~1333.
To test that explanation, the completeness limits of the
census of IC~348 and NGC~1333 need to be extended to higher extinctions.

\item
We have constructed H-R diagrams for IC~348 and NGC~1333 in terms of spectral
type and $M_K$. For the adopted distances of 300 and 235~pc, the median
sequences of the clusters coincide, which suggests that they have similar ages.
In contrast, NGC~1333 shows strong evidence of a younger age in the form of
higher abundances of disk-bearing stars and protostars and greater obscuration.
This discrepancy may indicate that IC~348 is closer or NGC~1333 is more
distant than we have assumed. The {\it Gaia} mission should soon test this
explanation by measuring parallactic distances for the clusters.

\item
Based on mid-IR photometry, the fraction of members that have circumstellar
disks is higher NGC~1333 than in IC~348 (N(II)/N(II+III)$\sim$0.6 and 0.4),
which agrees with results for earlier samples of members.
In each cluster, the disk fraction is roughly constant across the entire
range of stellar masses (0.01--3~$M_\odot$).

\item
For each cluster, we have examined the spatial distribution of stellar masses
and ages by computing the median spectral types and median offsets in $M_K$
from the median cluster sequence in three sections with differing stellar
densities. We also have measured the disk fraction in each of these areas.
None of these parameters exhibit significant variations with stellar density.

\end{enumerate}

\acknowledgements

This work was supported by grant AST-1208239 from the NSF.
We thank Konstantin Getman for providing his X-ray catalog of
NGC~1333 and Catarina Alves de Oliveira for providing her near-IR spectra of
candidate members of IC~348. We also thank Cameron Bell, Catarina Alves
de Oliveira, Eric Feigelson, Konstantin Getman, Charles Lada,
and Eric Mamajek for helpful discussions and comments on the manuscript.
The IRTF is operated by the University of Hawaii under contract
NNH14CK55B with NASA. The Gemini data were obtained through programs
GN-2008B-Q-21, GN-2014B-Q-55, GN-2015B-Q-43, and GN-2015B-FT-10.
Gemini Observatory is operated by AURA under a cooperative agreement with
the NSF on behalf of the Gemini partnership: the NSF (United States), the NRC
(Canada), CONICYT (Chile), the ARC (Australia),
Minist\'{e}rio da Ci\^{e}ncia, Tecnologia e Inova\c{c}\~{a}o (Brazil) and
Ministerio de Ciencia, Tecnolog\'{i}a e Innovaci\'{o}n Productiva (Argentina).
2MASS is a joint project of the University of Massachusetts and IPAC
at Caltech, funded by NASA and the NSF.  This work used data from
the NASA/IPAC Infrared Science Archive, operated by JPL under contract
with NASA, and the SIMBAD database, operated at CDS, Strasbourg, France.
The Digitized Sky Survey was produced at the Space Telescope Science
Institute under U.S. Government grant NAG W-2166. The images of these
surveys are based on photographic data obtained using the Oschin Schmidt
Telescope on Palomar Mountain and the UK Schmidt Telescope. The plates
were processed into the present compressed digital form with the
permission of these institutions.
WIRCam is a joint project of CFHT, Taiwan, Korea, Canada, and France.
MegaCam is a joint project of CFHT and CEA/DAPNIA.
CFHT is operated by the NRC of Canada, the Institute
National des Sciences de l'Univers of the Centre National de la Recherche
Scientifique of France, and the University of Hawaii.
Subaru Telescope is operated by National Astronomical Observatory of Japan.
The W.M. Keck Observatory is operated as a scientific partnership among
Caltech, the University of California, and NASA. The Observatory was made
possible by the generous financial support of the W.M. Keck Foundation.
The Center for Exoplanets and Habitable Worlds is supported by the
Pennsylvania State University, the Eberly College of Science, and the
Pennsylvania Space Grant Consortium.

Solar and Stellar Astrophysics (astro-ph.SR); Astrophysics of Galaxies (astro-ph.GA)