\section{introduction}

The cosmic microwave background (CMB) temperature anisotropy contains a wealth of cosmological information and has played a pivotal role in our understanding of the Universe. Besides the primordial fluctuations, various secondary anisotropies,  e.g.\ gravitational lensing, the thermal Sunyaev-Zel'dovich effect,  the kinetic Sunyaev-Zel'dovich effect, as well as the integrated Sachs-Wolfe effect, are playing an increasingly important role in constraining cosmological constituents and dynamics.

Among the secondary effects imprinted on the CMB gravitational lensing is of great importance. The projected gravitational lensing potential is a line-of-sight probe which contains information about the geometric distance traversed by CMB photons and time-dependent gravitational potentials. As such it is very sensitive to  late universe parameters, such as the sum of neutrino masses, the dark energy equation of state and spatial curvature. Since the projected gravitational lensing potential contains both geometric and structure growth information, it effectively breaks the angular diameter distance degeneracy~\cite{Smith:2008an}. Gravitational lensing measurements can also be used to de-lens the $B$-mode polarization of the CMB~\cite{Smith:2010gu}, enabling us to learn about primordial gravitational waves~\cite{Kamionkowski:1996zd} and the energy scale of inflation.

Tentative CMB weak lensing searches have been done with WMAP-7 year data sets~\cite{Smidt:2010by,Feng:2011jx} using non-Gaussian statistics. However, WMAP-7 alone cannot detect weak lensing of the CMB because WMAP temperature maps have insufficient sensitivity~\cite{Feng:2011jx}. Recently, the Atacama Cosmology Telescope ~\cite{Das:2011ak} and South Pole Telescope (SPT)~\cite{vanEngelen:2012va} have performed the first internal lensing reconstruction detections using non-Gaussianity. In addition, Atacama Cosmology Telescope and SPT also measured the gravitational lensing signal from the smoothing effects of the acoustic peaks on the CMB temperature power spectrum~\cite{Reichardt:2008ay,Das:2010ga,Keisler:2011aw}. As the experimental sensitivity improves, internal measurements, either from the power spectrum or the trispectrum, will become more precise in the near future.

The correlation between lensing and large scale structure arises from large scale structure, which deflects CMB photons in the late universe. The signal-to-noise ratio of lensing measurements can be enhanced if the CMB maps are cross correlated with highly sensitive large scale structure tracers, such as luminous red galaxies (LRGs) (which cover the redshift range $0.2<z<0.7$), quasars (which covers the redshift range $z<2.7$) from the Sloan Digital Sky Survey (SDSS),  or the NVSS of radio galaxies which has a higher mean redshift ($z\sim1$) than the LRGs and quasars. Hirata \textit{et al}.~\cite{Hirata:2004rp} used the cross correlation between WMAP-1 and LRGs and quasars from SDSS imaging and found no statistically significant signal. Then Smith \textit{et al}.~\cite{Smith:2007rg} used the cross correlation between WMAP-3 and NVSS, and found a $3.4\sigma$ signal, including systematics. Using a slightly less optimal estimator than Ref.~\cite{Smith:2007rg}, Hirata \textit{et al}.~\cite{Hirata:2008cb} obtained results consistent with, though at slightly lower significance than, Ref.~\cite{Smith:2007rg} for WMAP-3 with LRGs (0.95$\sigma$), WMAP-3 with quasars (1.64$\sigma$), and WMAP-3 with NVSS (2.13$\sigma$) respectively. Recently, SPT found a greater than $4\sigma$ cross correlation between the SPT convergence field and the galaxy survey from the Blanco Cosmology Survey, the Wide-field Infrared Survey Explorer, and Spitzer~\cite{Bleem:2012gm}. In this work, we use WMAP data released in years 1, 3, 5, 7, with NVSS to probe the lensing-galaxy correlation. We follow the methods developed in  Smith \textit{et al}.~\cite{Smith:2007rg} and Hirata \textit{et al}.~\cite{Hirata:2008cb}  using all of WMAP's datasets and compare our results to these earlier analyses.

The structure of the paper is as follows. We introduce the data sets in Sec.~\ref{sec:data}. Gravitational lensing effects on the CMB as well as the lensing extraction technique  are reviewed in Sec.~\ref{sec:lensing}. We describe the cross correlation estimators in Sec.~\ref{sec:cross}, and the forecast for Planck in Sec.~\ref{sec:forcast}. We discuss our results in Sec.~\ref{sec:con}.

\begin{figure*}
\includegraphics[width=6cm]{1.eps}\qquad
\includegraphics[width=6cm]{2.eps}
\caption{WMAP Kp0 mask (left) with $f_{\rm sky}=0.77$ and NVSS mask (right) with $f_{\rm sky}=0.573$.}\label{mask}
\end{figure*}
\section{Data Sets:  WMAP and NVSS}\label{sec:data}

The CMB data we use are from WMAP's Q-, V-, and W-band raw differencing assemblies (DAs). All of these DAs are masked by the Kp0 mask (Fig.~\ref{mask}) to remove bright sources and the galactic plane leaving $77\%$ of the sky.

The input for the galaxy distribution is the NVSS of radio galaxies. The NVSS~\cite{Condon:1998iy} team provides the software ``NVSSlist" to convert its raw catalog to a deconvolved one which is corrected for known biases and systematic errors. We use the deconvolved catalog to extract the galaxy count map. We use this software here, without specifying  either a minimum or maximum flux cut. The NVSS map is pixelized with a  {\sc HEALPix}\ pixelization scheme with $N_{\rm side}=256$. We remove the galactic plane ($|b|<10^\circ$) and the part of the sky unobserved by the survey ($\delta<-36.87^\circ$).  We also carefully remove bright sources with flux $> 1\,\rm Jy$ and mask out a disk of radius $1^{\circ}$ around them, forming the NVSS mask shown in Fig.~\ref{mask}. The resulting galaxy count map has $1,224,990$ sources, a sky fraction $f_{\rm sky}=0.573$, a mean number of sources per pixel $\bar n=2.72$, and a surface density of 170,249 galaxies per steradian. This agrees well with previous studies~\cite{Liu:2010re}.

\section{Gravitational lensing of the CMB}
\label{sec:lensing}

The effect of lensing on the CMB's primordial temperature $\tilde{T}$ in direction ${\bf n}$ can be represented by
\begin{equation}
T(\mathbf{n})=\widetilde{T}(\mathbf{n}+\mathbf{d}(\mathbf{n})),
\label{cmblensing}
\end{equation}
where $T$ is the lensed temperature and the deflection angle field $\mathbf{d}(\mathbf{n}) = \mathbf{\nabla} \phi$, with $\phi$ being the lensing potential. The operator $\nabla$ is the covariance derivative on the sphere with respect to
the angular position ${\bf n}$. We use Gaussian natural units with $\hbar=c=1$ throughout this paper.

The two-point correlation function of the temperature field is~\cite{Okamoto:2003zw}:
\begin{eqnarray}
\left\langle T_{lm} T_{l'm'} \right\rangle = \tilde{C}^{TT}_l\delta_{ll'}\delta_{m-m'} (-1)^m + \hspace{15mm} \nonumber \\
\sum_{LM}(-1)^M \left(
\begin{array}{ccc}
l   &  l'  & L\\
m & m'&-M
\end{array} \right) f^{TT}_{lLl'}\ \phi_{LM}, \label{lensingcouplemode}
\end{eqnarray}
where the second term encodes the effects of lensing with the weighting factor $f^{TT}_{lLl'}$ given by
\begin{eqnarray}
f^{TT}_{lLl'} = \tilde{C}_l^{TT} {} _0F_{l'Ll} + \tilde{C}_{l'}^{TT}
{} _0F_{lLl'}. \label{weighting}
\end{eqnarray}
We use the standard spherical harmonic decomposition
\begin{equation}
T({\bf n})=\displaystyle\sum_{lm}T_{lm}Y_{lm}({\bf n}),
\end{equation}
which defines the temperature modes $T_{lm}$. We use a similar notation for all other quantities defined on a sphere, e.g.\
 $\phi_{LM}$ are the modes of the lensing potential, etc.

Here $\tilde{C}_l^{TT}$ is the unlensed temperature power spectrum, and
\begin{eqnarray}
_0F_{lLl'} = \sqrt{\frac{(2l+1)(2l'+1)(2L+1)}{4\pi}} \times \hspace{15mm} \nonumber\\
\frac{1}{2} [L(L+1)+l'(l'+1)-l(l+1)] \left(
\begin{array}{ccc}
l  &  L &l'\\
0 & 0 & 0
\end{array} \right)
\end{eqnarray}
are proportional to the Wigner $3j$ symbols. Equation~(\ref{lensingcouplemode}) provides a way to extract $\phi_{LM}$ from the $TT$ correlations.

In the late universe, the Poisson equation relates the lensing potential $\Phi(\mathbf{k})$ to the density contrast $\delta(\mathbf{k})$,
\begin{equation}
k^2\Phi(\mathbf{k})=\frac{3H_0^2\Omega_m}{2a}\delta(\mathbf{k}),
\label{possion}
\end{equation}
where $\Omega_m$ is the matter fraction, $a$ is the scale factor and $H_0$ is the Hubble constant. Using the definition $D(\chi)=-2(1/\chi-1/\chi_*)$,  the projected lensing potential can be expressed as an integral along the line-of-sight,
\begin{equation}
\phi({\bf n})=\int_0^{\chi_*}
d\chi\ \Phi(\chi{\bf n})D(\chi),\label{realPhi}
\end{equation}
and it is integrated from 0 to the comoving distance at the last scattering surface $\chi_*$.
Here $\chi(z)$ is the comoving distance at redshift $z$. The galaxy overdensity is also given by a line-of-sight integration as
\begin{equation}
g({\bf n})=\frac{\int d\chi\ b_g\, \mathcal {N}(\chi)\delta(\chi{\bf
n})}{\int d\chi\ \mathcal {N}(\chi)}.\label{realG}
\end{equation}

To better understand the projected galaxy overdensity, $\mathcal {N}(\chi)=dN/d\chi$ is the comoving distance distribution of the galaxies. For NVSS galaxies,  there is a lack of accurate photometric redshifts so approximations to the redshift distributions are made~\cite{Xia:2011hj,Ho:2008bz,Schiavon:2012fc}. Following ~\cite{Smith:2007rg}, we use a Gaussian distribution
\begin{equation}
\frac{dN}{dz} \propto e^{-\frac{(z-z_0)^2}{2\sigma^2}}
\label{redshiftdistri}
\end{equation}
where $\sigma=0.8$ for $z<z_0 (=1.1)$, and $\sigma=0.3$ for $z>z_0$, and the comoving distribution $dN/d\chi$ is easily derived from Eq. (\ref{redshiftdistri}).

The parameter $b_g(z)$ is the redshift-dependent galaxy bias. To keep the model as simple as possible we treat the galaxy bias as a constant which can be determined from a fit to the data shown in Fig.~\ref{nvsscl}. This is different from~\cite{Hirata:2008cb,Ho:2008bz} which used the cross correlation of NVSS with the SDSS and with sources from the 2-Micron All Sky Survey to determine the galaxy bias. The galaxy bias is of great importance because the cross correlation can be directly translated into primordial non-Gaussianity~\cite{NbodyfNL} and may enable general relativity to be tested on cosmological scales~\cite{GR}.

Equations (\ref{realPhi}) and (\ref{realG}) give the general definitions for the lensing potential and galaxy overdensity which are both Gaussian fields, characterized by their variances, (i.e.\ the power spectra) $C_l^{\phi\phi}$ and $C_l^{gg}$. From Eq. (\ref{realPhi}) and (\ref{realG}), one sees that the lensing-galaxy cross correlation is built on the relation described by Eq. (\ref{possion}). Using the Limber approximation $k\sim l/\chi$ we calculate the theoretical galaxy auto-power spectrum and cross-power spectrum in Eqs.~(\ref{clgg}) and (\ref{clphig}).

The galaxy-galaxy power spectrum ( $\langle g_{lm}g^{\ast}_{lm}\rangle$ ) is
\begin{eqnarray}
C_l^{gg}&\simeq&\left(\frac{1}{\int d\chi \mathcal {N}(\chi)}\right)^2 \int d\chi\ b_g^2\,
\frac{\mathcal
{N}^2(\chi)}{\chi^2}P_{\delta}\left(\frac{l}{\chi}\right)\label{clgg}
\end{eqnarray}
and it will be used later for determining the galaxy bias and for simulating galaxy maps. This power spectrum shows the galaxy clustering strength on different angular scales.  We calculate the power spectrum of the NVSS overdensity map using two independent methods: a pseudo-$C_l$ method and a spherical harmonic estimation, as described in Refs.~\cite{Xia:2011hj, Blake:2004dg} respectively. We find {that both methods agree very well (Fig.~\ref{nvsscl}). As a final check, we have computed the NVSS galaxy power spectrum using the NVSS galaxy map in both equatorial and galactic coordinates and find a negligible difference, as expected, since the galaxy clustering is an intrinsic property of the Universe and does not depend on the choice of coordinate system. The galaxy power spectrum we use, obtained using the pixelized map in equatorial coordinates, is plotted in Fig.~\ref{nvsscl}.

The lensing-galaxy cross-power spectrum is
\begin{eqnarray}
C_l^{\phi g}&\simeq&\frac{3H_0^2\Omega_m}{2}\frac{1}{\int d\chi\
\mathcal {N}(\chi)}\nonumber\\&\times&\int d\chi\
b_g(1+z)D(\chi) \frac{\mathcal
{N}(\chi)}{k^2 \chi^2}P_{\delta}\left(\frac{l}{\chi}\right)
\label{clphig}
\end{eqnarray}
which shows the mutual influence between the gravitational potential and the galaxy clustering in the late universe on different angular scales. $P_{\delta}(k)$ is the matter power spectrum defined using the same convention as Ref.~\cite{pk}.  The primordial scalar curvature perturbations are evaluated at the pivot scale $k_0 = 0.002\, \textrm{Mpc}^{-1}$. The cross-power spectrum $C_l^{\phi g}$ will be used to simulate the correlated galaxy maps and will also be fit to data to determine the detection significance.

\begin{figure}
\includegraphics[bb=50 195 465 480,width=8.5cm]{3.eps}
\caption{The NVSS galaxy auto-power spectrum. The $1\sigma$ error bars are from 1000 Monte Carlo simulations of the NVSS galaxy map with galaxy bias $b_g=1$. For both sets of data points, the red points are from the pseudo-$C_l$ method~\cite{Xia:2011hj}, and the green are from the spherical harmonic estimation~\cite{Blake:2004dg}. The theoretical galaxy auto-power spectrum is fit to the red data points derived from the pseudo-$C_l$ method. Both methods show a consistent galaxy auto-power spectrum from the NVSS data. The first bin of the red data points largely deviates from the theoretical curve due to systematic effects.}
\label{nvsscl}
\end{figure}

\section{Cross Correlation Estimation}
\label{sec:cross}

\begin{table*}
\caption{The 6-parameter $\Lambda$CDM model used for the simulations of the temperature, galaxy and lensing potential. The derived parameter $\sigma_8$, based on the 6-parameter model, is shown in column eight. Using these parameters 1000 galaxy simulations with $b_g=1$ were performed to get the reconstructed galaxy biases as well as the $1\sigma$ error bars. From column nine, we see that all the reconstructed galaxy biases are consistent with the input value $b_g=1$. Furthermore, the galaxy biases of the real data are calculated based on the simulations and shown in column ten. Two independent methods were used to calculate the galaxy auto-power spectra, as specified in the footnotes. }
\begin{tabular}{c|c|c|c|c|c|c|c|c|c}
\hline Data
set&$\Omega_bh^2$&$\Omega_{\rm CDM}h^2$&$H_0$&$A_s$&$n_s$&$\tau$&$\sigma_8$&$b^{\rm sim}_g$&$b^{\rm data}_g$\\
\hline
\multirow{2}{*}{}
WMAP-1\footnotemark[3]&0.0226&0.1104&72&$2.212\times10^{-9}$&0.96&0.117&0.76&$0.98\pm0.12$\footnotemark[1]&$2.00\pm0.12$\footnotemark[1]\\
&&&&&&&&$0.96\pm0.10$\footnotemark[2]&$1.95\pm0.10$\footnotemark[2]\\
\hline
\multirow{2}{*}{}WMAP-3\footnotemark[4]&0.02186&0.1105&70.4&$2.393\times10^{-9}$&0.947&0.073&0.77&$0.98\pm0.11$\footnotemark[1]&$1.97\pm0.11$\footnotemark[1]\\
&&&&&&&&$0.96\pm0.09$\footnotemark[2]&$1.92\pm0.09$\footnotemark[2]\\
\hline
\multirow{2}{*}{}WMAP-5\footnotemark[5]&0.02305&0.1182&69.7&$2.484\times10^{-9}$&0.969&0.094&0.85&$0.98\pm0.10$\footnotemark[1]&$1.85\pm0.10$\footnotemark[1]\\
&&&&&&&&$0.96\pm0.08$\footnotemark[2]&$1.81\pm0.08$\footnotemark[2]\\
\hline
\multirow{2}{*}{}WMAP-7\footnotemark[6]&0.02258&0.1109&71&$2.43\times10^{-9}$&0.963&0.088&0.80&$0.98\pm0.11$\footnotemark[1]&$1.91\pm0.11$\footnotemark[1]\\
&&&&&&&&$0.96\pm0.09$\footnotemark[2]&$1.86\pm0.09$\footnotemark[2]\\
\hline
\end{tabular}
\footnotetext[1]{Pseudo-$C_l$ method~\cite{Xia:2011hj}, }
\footnotetext[2]{Spherical harmonic estimation~\cite{Blake:2004dg}}
\footnotetext[3]{WMAP1+CBI+ACBAR+2dFGRS+Ly$\alpha$, ~\cite{wmap1para}}
\footnotetext[4]{WMAP3 ALL, \url{http://lambda.gsfc.nasa.gov/product/map/dr2/parameters.cfm}}
\footnotetext[5]{WMAP5+BAO+SNALL+Ly$\alpha$POST, \url{http://lambda.gsfc.nasa.gov/product/map/dr3/parameters.cfm}}
\footnotetext[6]{WMAP7, \url{http://lambda.gsfc.nasa.gov/product/map/dr4/parameters.cfm}} \label{paratable}
\end{table*}

Monte Carlo simulations are used to estimate the cross correlation between the CMB and the galaxy distribution. The variances ${\tilde C}_l^{TT}$, $C_l^{TT}$, $C_l^{\phi\phi}$ are computed using {\sc CAMB}~\cite{Lewis:1999bs} with the cosmological parameters listed in Table~\ref{paratable}. In addition to these, we derive $C_l^{gg}$ and $C_l^{\phi g}$ from Eq. (\ref{clgg}) and Eq. (\ref{clphig}) with the parameters listed in Table~\ref{paratable}. Simulated CMB temperature modes, ${\tilde a}_{lm}$, are drawn from Gaussian distributions with zero means and variances ${\tilde C}_l^{TT}$. In this work, we will use two sets of cosmological parameters because we want to check the consistency of our results with a previous study~\cite{Smith:2007rg} and also because we want to explore the impact of using the newest parameters from WMAP-7 on the lensing-galaxy cross correlations.

We convert these ${\tilde a}_{lm}$ to an unlensed temperature map, ${\tilde T}({\bf n})$, on which we do a cubic interpolation to precisely implement Eq. (\ref{cmblensing}). This produces a lensed temperature map $T({\bf n})$ that is converted back to harmonic space  to give the lensed modes $a_{lm}$. Then each DA's beam and pixel window transfer function (the pixel window transfer function has negligible effects on the cross-power spectra) from WMAP are multiplied by these modes which are subsequently transformed into a temperature map containing the lensing signal.

We then simulate Gaussian noise in map space where the pixel noise is assumed to be uncorrelated and Gaussian distributed with zero mean and pixel-independent variance determined from ${\sigma_0}/{\sqrt{N_{\rm obs}}}$.  Here, both $\sigma_0$ and $N_{\rm obs}$ are supplied by the WMAP team for different DAs. We add this noise map to the signal map and apply the Kp0 mask to get a simulated WMAP DA made in the same way as the real WMAP maps were produced. The entire procedure can be summarized by Eq.~(\ref{simWMAP}) in which $a_{lm}$ is the lensed CMB mode,  $n({\bf n})$ the white noise, ${\rm M}^{\rm WMAP}({\bf n})$ the mask map, $\nu$ the index of the DA channel, $p_l$ the pixel window transfer function, $b_l$ the beam transfer function
\begin{eqnarray}
T^{(\nu)}({\bf n})&=&{\rm M}^{\rm WMAP}({\bf n})\biggl[\sum_{lm}p_lb^{(\nu)}_la_{lm}Y_{lm}({\bf n})\nonumber\\
&&\qquad + \left(\frac{\sigma_0}{\sqrt{N_{\rm obs}({\bf
n})}}\right)^{(\nu)}n({\bf n})\biggr]\,.
\label{simWMAP}
\end{eqnarray}

To maximize the signal-to-noise ratio, we compute a single harmonic mode ${\hat a}_{lm}$ from eight Q, V, W-band DAs. This reduction step is expressed as~\cite{Smith:2007rg}
\begin{eqnarray}
\mathbf{\hat a}&=&(\mathbf{S}+\mathbf{N})^{-1}\mathbf{a}\nonumber\\
&=&\mathbf{S}^{-1/2}\mathbf{A}^{-1}\mathbf{S}^{1/2}\mathbf{N}^{-1}\mathbf{a}.\label{DAto1}
\end{eqnarray}
Here ${\bf a}$ is the vector of $a_{lm}$s, ${\bf S}$ the signal covariance matrix, ${\bf N}$ the noise covariance matrix, and ${\bf A}=\mathbf{I}+\mathbf{S}^{1/2}\mathbf{N}^{-1}\mathbf{S}^{1/2}$.  We use the second form of Eq.~(\ref{DAto1})  and filter the raw CMB modes using a multigrid-preconditioned-complex conjugate gradient method. The master equation that has to be solved is
\begin{equation}
\mathbf{A}\mathbf{x}=\mathbf{y},
\label{eq13}
\end{equation}
where $\mathbf{x}=\mathbf{S}^{1/2}\mathbf{\hat a}$, and
$\mathbf{y}=\mathbf{S}^{1/2}\mathbf{N}^{-1}\mathbf{a}$.
Equation~(\ref{eq13}) is better for numerical computations because ${\bf A}$ is close to the unit matrix. Appendix~\ref{app:mpcg} gives details of the numerical calculation of Eq.~(\ref{eq13}). We solve Eq.~(\ref{eq13}) with $\hat a_{lm}$ for both the temperature ($\hat T_{lm}$) and the galaxy ($\hat g_{lm}$).

We use the standard quadratic estimator to reconstruct a noisy lensing potential map $\hat\phi_{lm}$ in harmonic space~\cite{Hu:2001fa,Hu:2001tn},
\begin{equation}
\displaystyle\sum_{lm}\hat\phi_{lm}Y_{lm}({\bf n})=\nabla^i({\
}_0A^T({\bf n})\nabla_i{\ }_0B^T({\bf n})),\label{noisypotential}
\end{equation}
where \begin{equation} {}_0A^T(\mathbf{n}) =
\displaystyle\sum_{lm}\hat T_{lm} Y_{lm}(\mathbf{n})\end{equation} and
\begin{equation}
_0B^T({\bf{n}}) = \displaystyle \sum_{lm} \tilde{C}^{TT}_l\hat
T_{lm}  Y_{lm}({\bf{n}}).
\end{equation}
In the above equations, $\nabla_i$ is the gradient operator on a sphere and $\nabla^i=g^{ij}\nabla_j$ . Here $g_{ij}$ is the metric of a sphere.

We also use Monte Carlo simulations for the NVSS galaxy maps. The simulated galaxy modes $g_{lm}$ are drawn from a Gaussian distribution and transformed into a galaxy overdensity map $g({\bf n})$ at  {\sc HEALPix}\ resolution $N_{\rm side}=1024$. The galaxy modes must satisfy the correct galaxy-galaxy auto-correlation and lensing-galaxy cross correlation. From these two constraints the simulated galaxy mode must be
\begin{equation}
g_{lm}=\frac{C_l^{\phi
g}}{C_l^{\phi\phi}}\phi_{lm}+\sqrt{C_l^{gg}-\frac{(C_l^{\phi
g})^2}{C_l^{\phi\phi}}}G_{lm},
\end{equation}
where $G_{lm}$ is a complex Gaussian random variable, and $\phi_{lm}$ is inherited from the deflection field in Eq. (\ref{cmblensing}). From this equation we see that the lensing-galaxy correlation is encoded in the first term. A NVSS map is generated where the galaxy number count in each pixel is drawn from a Poisson distribution with mean
\begin{equation}
\lambda({\bf n})=\bar n(1+g({\bf n})).
\end{equation}
The galaxy count map $\lambda({\bf n})$ is used to generate a simulated galaxy overdensity
map $g^{(\text{sim})} ({\bf n})$,
\begin{equation}
g^{(\text{sim})} ({\bf n}) ={\rm{M}}^{\rm{NVSS}}({\bf n})\left[\frac{\lambda({\bf
n})}{\bar n}-1\right],
\end{equation}
where ${\rm{M}}^{\rm{NVSS}}$ is the NVSS mask shown in Fig.~\ref{mask}. $g^{(\text{sim})} ({\bf n})$ automatically contains the shot-noise with the variance $N_l^{gg}={1}/{\bar n}$ for the galaxy overdensity map. We degrade this map to resolution $N_{\rm side}=256$ i.e.\ the same as the real NVSS data. The harmonic mode $g_{lm}^{(\text{sim})}$, which contains the shot-noise, is obtained from $g^{(\text{sim})} ({\bf n})$ and is further filtered using the same procedure as in Eq.~(\ref{DAto1}),
\begin{equation}
\hat g_{lm} =({\bf S}+{\bf N})^{-1}g_{lm}^{(\text{sim})} .\label{rawgmap}
\end{equation}
Here ${\bf S}$ represents the primordial galaxy covariance and ${\bf N}$ is the shot-noise covariance.

We show the noisy reconstruction of the potential maps and the filtered galaxy map in Fig.~\ref{tracermap}, using the measured WMAP and NVSS data.
\begin{figure*}
\includegraphics[width=8cm]{4.eps}
\vspace{0.5cm}
\includegraphics[width=8cm]{5.eps}
\caption{The noisy reconstruction of the lensing potential map
(Eq.(\ref{noisypotential}) from WMAP-7) band-pass filtered from $20\le l\le 40$ (left). The analogous map from NVSS galaxy data [Eq.(\ref{rawgmap})] within the band $20\le l\le 40$ (right). } \label{tracermap}
\end{figure*}

The lensing-galaxy cross-power spectrum is the observable which will be compared with the counterpart from data. The estimator of the lensing-galaxy cross correlation is expressed as
\begin{equation}
C^{\phi g}_b=\frac{1}{\mathcal
{F}_b}\displaystyle\sum_{\substack{l\in b\\ -l\le m\le l}} (\hat
\phi_{lm}-\langle\hat\phi_{lm}\rangle)^{\ast}\hat
{g}_{lm},\label{cross-estimator}
\end{equation}
where $\mathcal {F}_b$ is the normalization factor at band $b$. It is shown in Ref.~\cite{Smith:2006ud} that the normalization factor can be calculated by either a
direct or a fast method for the full-sky coverage and that the fast method converges in $O(10^2)$ simulations. When there is a sky-cut these methods account for the sky-cut effect very well and a constant $f_{\rm sky}$ is often used. The factor $f_{\rm sky}$ is actually a function of $l$~\cite{Verde:2003ey} so a simple constant approximation may potentially bias the cross-spectra reconstruction. Therefore an end-to-end simulation~\cite{Smith:2007rg} is the best way to get the exact normalization accounting for the sky-cut and that is done here.

\begin{table*}
\caption{The two sets of cosmological parameters used in this work:  we choose two sets of parameters (labeled ``{\bf Set-3}" and ``{\bf Set-7}") to do the cross correlation calculations in this work. In order to compare our results with those from the previous studies~\cite{Smith:2007rg}, we use the parameters they used {\bf Set-3}\ from WMAP-3's cosmological parameters ( row ``WMAP-3" in Table \ref{paratable}) combined with the corresponding galaxy bias of Smith \textit{et al}.~\cite{Smith:2007rg}.  Based on the newest cosmological parameters from WMAP-7 ( row ``WMAP-7" in Table \ref{paratable}) we construct a new parameter set, {\bf Set-7}\ with the corresponding galaxy bias shown in Table \ref{paratable}. }
\begin{tabular}{c|c|c|c|c|c|c|c|c}
\hline Data
set&$\Omega_bh^2$&$\Omega_{\rm CDM}h^2$&$H_0$&$A_s$&$n_s$&$\tau$&$\sigma_8$&$b_g$\\
\hline
{\bf Set-3} &0.02186&0.1105&70.4&$2.393\times10^{-9}$&0.947&0.073&0.77&$1.70$\\
\hline
{\bf Set-7} &0.02258&0.1109&71&$2.43\times10^{-9}$&0.963&0.088&0.80&$1.91$\\
\hline
\end{tabular}
\label{paratable1}
\end{table*}

\begin{table*}
\centering \caption{Measure of lensing-galaxy cross correlation
$\mathcal{C}$ and its significance $\mathcal{C} / \Delta
\mathcal{C}$ using  {\bf Set-3}. For five columns of this table: the second column shows the simulation results, the third column is the case without gradient stripes removed, the fourth column is the case with gradient stripes removed (this column shows the statistical results of the lensing-galaxy cross correlations). The fifth column is the case by setting the NVSS map in equatorial coordinates as a complementary diagnostic test. }
\begin{tabular}{c|cc|cc|cc|c}
\hline\hline
Data set & $\mathcal {C}^{\rm{sim}}$ & $\mathcal{C} / \Delta \mathcal{C}$& $\mathcal {C}$\footnotemark[1]& $\mathcal{C} / \Delta \mathcal{C}$&$\mathcal {C}$\footnotemark[2] & $\mathcal{C} / \Delta \mathcal{C}$& $\mathcal {C}$\footnotemark[3]\\
\hline
WMAP-1$\times$NVSS& $1.00\pm0.47$ & $2.13\sigma$& $1.25\pm0.47$ & $2.66\sigma$&$1.24\pm0.47$&$2.64\sigma$&0.26\\
\hline
WMAP-3$\times$NVSS& $1.00\pm0.35$ & $2.86\sigma$& $1.20\pm0.35$ & $3.43\sigma$&$1.26\pm0.35$&$3.60\sigma$&0.17\\
\hline
WMAP-5$\times$NVSS& $1.00\pm0.31$ & $3.23\sigma$& $1.24\pm0.31$ & $4.00\sigma$&$1.27\pm0.31$&$4.10\sigma$&0.23\\
\hline
WMAP-7$\times$NVSS& $1.00\pm0.30$ & $3.33\sigma$& $1.14\pm0.30$ & $3.80\sigma$&$1.16\pm0.30$&$3.87\sigma$&0.15\\
\hline
\end{tabular} \\
\footnotetext[1]{Without gradient stripes removed.}
\footnotetext[2]{With gradient stripes removed.}
\footnotetext[3]{Galaxy map in equatorial coordinate.}
\label{lensingtable1}
\end{table*}

\begin{table*}
\centering \caption{Measure of lensing-galaxy cross correlation
$\mathcal{C}$ and its significance $\mathcal{C} / \Delta
\mathcal{C}$ using {\bf Set-7}. The format of this table is the same as Table \ref{lensingtable1}. }
\begin{tabular}{c|cc|cc|cc|c}
\hline\hline
Data set & $\mathcal {C}^{\rm{sim}}$ & $\mathcal{C} / \Delta \mathcal{C}$& $\mathcal {C}$\footnotemark[1]& $\mathcal{C} / \Delta \mathcal{C}$&$\mathcal {C}$\footnotemark[2] & $\mathcal{C} / \Delta \mathcal{C}$& $\mathcal {C}$\footnotemark[3]\\
\hline
WMAP-1$\times$NVSS& $1.00\pm0.41$ & $2.44\sigma$& $1.01\pm0.41$ & $2.46\sigma$&$1.00\pm0.41$&$2.44\sigma$&0.20\\
\hline
WMAP-3$\times$NVSS& $1.00\pm0.31$ & $3.23\sigma$& $0.96\pm0.31$ & $3.10\sigma$&$1.01\pm0.31$&$3.26\sigma$&0.13\\
\hline
WMAP-5$\times$NVSS& $1.00\pm0.28$ & $3.57\sigma$& $0.98\pm0.28$ & $3.50\sigma$&$1.01\pm0.28$&$3.61\sigma$&0.18\\
\hline
WMAP-7$\times$NVSS& $1.00\pm0.26$ & $3.85\sigma$& $0.92\pm0.26$ & $3.54\sigma$&$0.93\pm0.26$&$3.58\sigma$&0.11\\
\hline
\end{tabular} \\
\footnotetext[1]{Without gradient stripes removed.}
\footnotetext[2]{With gradient stripes removed.}
\footnotetext[3]{Galaxy map in equatorial coordinate.}
\label{lensingtable2}
\end{table*}

\begin{table}
\centering \caption{Fisher matrix analysis for WMAP$\times$NVSS cross
correlation. The $1\sigma$ error bars are determined from Eq. (\ref{fisherErrorbar}). We calculate two sets of the optimal bounds for this work, based on two sets of parameters: {\bf Set-3}\ (column two); {\bf Set-7}\ (column three). }
\begin{tabular}{c|cc|cc}
\hline
Data set &$\mathcal {C}^{\textrm{optimal}}$\,\footnotemark[1]&$\mathcal{C} / \Delta \mathcal{C}$&$\mathcal {C}^{\textrm{optimal}}$\,\footnotemark[2]&$\mathcal{C} / \Delta \mathcal{C}$\\
\hline
WMAP-1$\times$NVSS&$1\pm0.46$&$2.17\sigma$&$1\pm0.39$&$2.56\sigma$\\
\hline
WMAP-3$\times$NVSS&$1\pm0.29$&$3.45\sigma$&$1\pm0.25$&$4.00\sigma$\\
\hline
WMAP-5$\times$NVSS&$1\pm0.25$&$4.00\sigma$&$1\pm0.21$&$4.76\sigma$\\
\hline
WMAP-7$\times$NVSS&$1\pm0.22$&$4.55\sigma$&$1\pm0.19$&$5.26\sigma$\\
\hline
\end{tabular}\\
\footnotetext[1]{WMAP-3 year cosmological parameters and $b_g=1.70$.}
\footnotetext[2]{WMAP-7 year cosmological parameters and
$b_g=1.91$.}
\label{opboundWMAPNVSS}
\end{table}

\begin{table*}
\centering \caption{Results of the curl null tests for WMAP$\times$NVSS cross correlation. The curl null tests are performed based on two sets of parameters: {\bf Set-3}\ (column two); {\bf Set-7}\ (column three).  }
\begin{tabular}{c|cc|cc}
\hline
Data set &$\mathcal {C}$\footnotemark[1]&$\mathcal{C} / \Delta \mathcal{C}$&$\mathcal {C}$\footnotemark[2]&$\mathcal{C} / \Delta \mathcal{C}$\\
\hline
WMAP-1$\times$NVSS&$-0.11\pm0.47$&$-0.23\sigma$&$-0.03\pm0.41$&$-0.07\sigma$\\
\hline
WMAP-3$\times$NVSS&$0.00\pm0.35$&$0.00\sigma$&$0.04\pm0.31$&$0.13\sigma$\\
\hline
WMAP-5$\times$NVSS&$0.05\pm0.31$&$0.16\sigma$&$0.07\pm0.28$&$0.25\sigma$\\
\hline
WMAP-7$\times$NVSS&$-0.05\pm0.30$&$0.17\sigma$&$-0.03\pm0.26$&$-0.12\sigma$\\
\hline
\end{tabular}\\
\footnotetext[1]{WMAP-3 year cosmological parameters and $b_g=1.70$.}
\footnotetext[2]{WMAP-7 year cosmological parameters and
$b_g=1.91$.}
\label{curltable}
\end{table*}

\begin{table*}
\centering \caption{Gaussianity diagnostics for the probability distribution of $\{\mathcal{C}\}$ which is constructed from 1000 Monte Carlo simulations. The second column is the Kolmogorov-Smirnov test, and the critical value is 0.04 at 5\% confidence level. The Kolmogorov-Smirnov test requires the maximum deviation be $<0.04$ to validate the distribution is Gaussian. The third column is the skewness of $\{\mathcal{C}\}$, and the fourth column is the kurtosis of $\{\mathcal{C}\}$. The upper values in the cells are the results for {\bf Set-3}, the lower values for {\bf Set-7}.  For a Gaussian distribution, the skewness should be 0 and the kurtosis should be 3. As can be seen, all the probability distribution functions pass the Kolmogorov-Smirnov test and are consistent with being Gaussian-distributed.  }
\begin{tabular}{c|c|c|c}
\hline
Data set &maximum distance[$<$0.04]&skewness[$\sim$0]&kurtosis[$\sim$3]\\
\hline
\multirow{2}{*}
{WMAP-1$\times$NVSS}&0.02&0.02&2.89\\
&0.02&-0.05&2.77\\
\hline
\multirow{2}{*}{WMAP-3$\times$NVSS}&0.02&-0.14&2.62\\
&0.03&-0.17&2.58\\
\hline
\multirow{2}{*}{WMAP-5$\times$NVSS}&0.03&-0.17&2.53\\
&0.03&-0.23&2.43\\
\hline
\multirow{2}{*}{WMAP-7$\times$NVSS}&0.03&-0.21&2.44\\
&0.03&-0.21&2.36\\
\hline
\end{tabular}\\
\label{kstable}
\end{table*}

As a systematic check we note that the lensing signal consists of a gradient and a curl component~\cite{curl2005}. The curl component estimator $\psi_{lm}$ is
defined by
\begin{equation}
\displaystyle\sum_{lm}\psi_{lm}Y_{lm}({\bf n})=\displaystyle
\sum_{ij} \epsilon^{ij}\nabla_i({\ }_0A^T({\bf n})\nabla_j{\
}_0B^T({\bf n}))\label{curl-estimator}
\end{equation}
and should vanish because lensing does not generate vorticity. Similar to Eq. (\ref{cross-estimator}), the curl-galaxy
cross correlation diagnostic is calculated by
\begin{equation}
C^{\psi g}_b=\frac{1}{\mathcal
{F}_b}\displaystyle\sum_{\substack{l\in b\\ -l\le m\le l}} (\psi_{lm}-\langle\psi_{lm}\rangle)^{\ast}\hat
{g}_{lm}\label{curl-cross}
\end{equation}
which should also vanish.

The amplitude of the cross correlation is determined using
\begin{equation}
\mathcal {C} = \frac{\sum_{AB} C^{(\textrm{th})}_A
\mathbf{C}^{-1}_{AB} C^{(\textrm{obs})}_{B}}
{\sum_{AB}C^{(\textrm{th})}_A\mathbf{C}^{-1}_{AB}C^{(\textrm{th})}_{B}}.\label{CL}
\end{equation}
${\bf C}_{AB}$ is the covariance matrix for the band powers and $A$ and $B$ stand for the band power index. We find that the off-diagonal correlations of ${\bf C}_{AB}$ are negligible, and the covariance matrix elements can be simply replaced by the band power variance $\sigma_A^2$, i.e.\ ${\bf C}_{AB}=\sigma_A^2\delta_{AB}$.

We have described the procedures used to perform analysis on simulated or measured data. Now we summarize the analysis of the real WMAP and NVSS data.
We fit the theoretical galaxy auto-power spectrum to the NVSS data in Fig.~\ref{nvsscl} and determine the galaxy biases (Table \ref{paratable}) using two methods. The error bars are determined from 1000 simulated galaxy maps with galaxy bias $b_g=1$. Then we choose two sets of parameters (labeled ``{\bf Set-3}" and ``{\bf Set-7}") in Table~\ref{paratable1} to do the cross correlation calculations in this work.  In order to compare our results with those from the previous studies~\cite{Smith:2007rg}, we use the parameters they used {\bf Set-3}\ from WMAP-3's cosmological parameters ( row ``WMAP-3" in Table \ref{paratable}) combined with the corresponding galaxy bias of Smith \textit{et al}.~\cite{Smith:2007rg}.  Based on the newest cosmological parameters from WMAP-7 ( row ``WMAP-7" in Table \ref{paratable}) we construct a new parameter set {\bf Set-7}\ with the corresponding galaxy bias shown in Table \ref{paratable}. For each of the parameter sets we calculate four lensing-galaxy cross correlations from WMAP-1 to WMAP-7.

We carefully treat the known systematics of NVSS, i.e.\ the gradient stripes which are generated by the declination-dependence of the galaxy overdensity field due to the low-flux calibration issue~\cite{Condon:1998iy}. We first make a gradient stripe map only using $m=0$ modes and then subtract it from the galaxy map. We calculate the lensing-galaxy cross correlations for two cases: without the gradient stripes removed and with the gradient stripes removed. We find that this systematic effect does not affect the lensing-galaxy cross correlations as seen from column $\mathcal C^a$ and column $\mathcal C^b$ in Table~\ref{lensingtable1} and Table~\ref{lensingtable2}. The statistical results are those with the gradient stripes removed which are shown in Figs.~\ref{crosspower} and ~\ref{crosspowerpara2} for the Kp0$+$NVSS mask combination. From the two figures, we find that the lensing-galaxy cross-power spectra are consistent with the theoretical predictions and the uncertainty of the cross-power spectrum is decreasing as the year of WMAP increases. All the error bars are calculated from 1000 Monte Carlo simulations, which we confirmed to be sufficient for convergences. As a complementary diagnostic test, we keep the NVSS galaxy overdensity map in equatorial coordinates and calculate the cross-power spectra and all the amplitudes are shown in column $\mathcal C^c$ in Table~\ref{lensingtable1} and Table~\ref{lensingtable2}. As can be seen, they are negligible. All the cross correlation amplitudes are summarized in Table~\ref{lensingtable1} and Table~\ref{lensingtable2}. From the results of  WMAP-3$\times$NVSS in Table~\ref{lensingtable1}: for the statistical results, we get lensing detection significance level of $3.60\sigma$ and ~\cite{Smith:2007rg} got $4.02\sigma$. Both analyses agree quite well. We find the cross-power spectra from WMAP-5 and WMAP-7 clearly and firmly show the lensing-galaxy correlation at $>3\sigma$ level for both cases. All the results are within the optimal bounds shown in Table~\ref{opboundWMAPNVSS}.

Assuming there is no cosmological parity violation the curl-galaxy correlation should be consistent with zero. We show the results of the curl null tests in Figs.~\ref{curlpower} and \ref{curlpowerpara2}. As expected, all the correlations are consistent with zero. The amplitude as well as the significance are given in Table~\ref{curltable}.

We pixelized the NVSS catalog with different  {\sc HEALPix}\ resolutions (e.g. $N_{\rm side}=512, 1024$)  to probe the possible pixel artifacts that could afflict the cross-power spectra and because we want to examine the impact of possible long range spatial correlations. However, we find that different NVSS pixelization resolutions do not affect the cross correlation.

We also use the diagonal elements of the covariance matrix to do the analysis to check the consistency with previous studies. In this case, the estimator has a larger variance as pointed out by Smith \textit{et al}.~\cite{Smith:2009jr}. This contributes to the difference in significance levels between $4\sigma$~\cite{Smith:2007rg} and $2\sigma$~\cite{Hirata:2008cb}.

\begin{figure*}
\includegraphics[bb=50 195 465 480,width=8.5cm]{6.eps}\qquad
\includegraphics[bb=50 195 465 480,width=8.5cm]{7.eps} \\
\vspace{0.5cm}
\includegraphics[bb=50 195 465 480,width=8.5cm]{8.eps}\qquad
\includegraphics[bb=50 195 465 480,width=8.5cm]{9.eps}
\caption{({\bf Set-3}) The lensing-galaxy cross-power spectra for WMAP$\times$ NVSS are calculated from Eq. (\ref{cross-estimator}). The Kp0 mask is used to remove the contaminated regions of the WMAP data. WMAP's data are provided from two Q bands, two V bands and four W bands. The NVSS mask is applied to the galaxy map to remove bright sources and unobserved regions. The theoretical cross-power spectra are shown in blue solid lines, and they are the same for all of the four panels. The real data are shown in the red scattered points. The statistical amplitude for WMAP-1$\times$NVSS is $1.24\pm0.47$, for WMAP-3$\times$NVSS is $1.26\pm0.35$, for WMAP-5$\times$NVSS is $1.27\pm0.31$, for WMAP-7$\times$NVSS is $1.16\pm0.30$. All the error bars are determined from 1000 Monte Carlo simulations. We find that the lensing-galaxy cross-power spectra are consistent with the theoretical predictions and the uncertainty of the cross-power spectrum is decreasing as the year of WMAP increases. }
\label{crosspower}
\end{figure*}
\begin{figure*}
\includegraphics[bb=50 195 465 480,width=8.5cm]{10.eps}\qquad
\includegraphics[bb=50 195 465 480,width=8.5cm]{11.eps} \\
\vspace{0.5cm}
\includegraphics[bb=50 195 465 480,width=8.5cm]{12.eps}\qquad
\includegraphics[bb=50 195 465 480,width=8.5cm]{13.eps}
\caption{({\bf Set-7}) The lensing-galaxy cross-power spectra for WMAP$\times$ NVSS are calculated from Eq. (\ref{cross-estimator}). See Fig.~\ref{crosspower} for detailed descriptions. The statistical amplitude for WMAP-1$\times$NVSS is $1.00\pm0.41$, for WMAP-3$\times$NVSS is $1.01\pm0.31$, for WMAP-5$\times$NVSS is $1.01\pm0.28$, for WMAP-7$\times$NVSS is $0.93\pm0.26$. }
\label{crosspowerpara2}
\end{figure*}

\begin{figure*}
\includegraphics[bb=50 195 465 480,width=8.5cm]{14.eps}\qquad
\includegraphics[bb=50 195 465 480,width=8.5cm]{15.eps} \\
\vspace{0.5cm}
\includegraphics[bb=50 195 465 480,width=8.5cm]{16.eps}\qquad
\includegraphics[bb=50 195 465 480,width=8.5cm]{17.eps}
\caption{({\bf Set-3}) The curl null tests for WMAP$\times$ NVSS are calculated from Eq. (\ref{curl-cross}). The Kp0 mask is used to remove the contaminated regions of the WMAP data. WMAP data are provided from two Q bands, two V bands and four W bands. The NVSS mask is applied to the galaxy map to remove bright sources and unobserved regions. The theoretical lensing-galaxy cross-power spectra with both WMAP and NVSS in galactic coordinates are shown in blue solid lines for comparison, and they are the same for all of the four panels. The curl amplitude for WMAP-1$\times$NVSS is $-0.11\pm0.47$, for WMAP-3$\times$NVSS is $0.00\pm0.35$, for WMAP-5$\times$NVSS is $0.05\pm0.31$, for WMAP-7$\times$NVSS is $-0.05\pm0.30$. As can be seen, all cross-power spectra for the curl null test are consistent with zero (black dotted line). }
\label{curlpower}
\end{figure*}
\begin{figure*}
\includegraphics[bb=50 195 465 480,width=8.5cm]{18.eps}\qquad
\includegraphics[bb=50 195 465 480,width=8.5cm]{19.eps} \\
\vspace{0.5cm}
\includegraphics[bb=50 195 465 480,width=8.5cm]{20.eps}\qquad
\includegraphics[bb=50 195 465 480,width=8.5cm]{21.eps}
\caption{({\bf Set-7}) The curl null tests for WMAP$\times$ NVSS are calculated from Eq. (\ref{curl-cross}). See Fig.~\ref{curlpower} for detailed descriptions. The curl amplitude for WMAP-1$\times$NVSS is $-0.03\pm0.41$, for WMAP-3$\times$NVSS is $0.04\pm0.31$, for WMAP-5$\times$NVSS is $0.07\pm0.28$, for WMAP-7$\times$NVSS is $-0.03\pm0.26$. As can be seen, all cross-power spectra for the curl null test are consistent with zero (black dotted line).}
\label{curlpowerpara2}
\end{figure*}

\section{Forecast for Future Experiments}
\label{sec:forcast}

\begin{figure}
\includegraphics[bb=50 195 465 480,width=8.5cm]{22.eps}
\caption{The signal-to-noise ratio for the lensing-galaxy cross correlation between Planck and NVSS as a function of the maximum multipole used in the analysis.} \label{planck-nvss}
\end{figure}
The revealed cross correlation between WMAP and NVSS hints that the detection significance would be further enhanced if the precision of the CMB data were improved. The upcoming Planck data will improve upon WMAP, so we expect that the cross correlation between Planck and NVSS will be more significant. To predict the optimal bound on the detection signal-to-noise ratio for lensing-galaxy cross correlation we first calculate the equivalent noise $N_l^{\phi g}$ from the following equation
\begin{equation}
N_l^{\phi g}=\left[N_l^{\phi\phi}N_l^{gg}\right]^{1/2}.\label{crossNoise}
\end{equation}
where $N_l^{\phi\phi}$ is the lensing reconstruction noise~\cite{Okamoto:2003zw} and $N_l^{gg}$ is the galaxy shot-noise. The efficient algorithm for calculating $N_l^{\phi\phi}$ is given in Refs.~\cite{Smith:2010gu,Feng:2011jx}. This reconstruction noise can be minimized by combining different CMB channels and the minimum noise is
\begin{equation}
N_l^{\rm min, \phi\phi}=\frac{1}{\sum_{\nu}\left[N_l^{\nu,
\phi\phi}\right]^{-1}},
\end{equation}
Both of the noise spectra effectively propagate the uncertainty $\Delta C_l^{\phi g}$ into the cross-power spectrum $C_l^{\phi g}$. Specifically, we express it as
\begin{equation}
\Delta C_l^{\phi g}=\sqrt{\frac{2}{(2l+1)f_{\rm sky}}}\ (C_l^{\phi
g}+N_l^{\phi g}).\label{fisherErrorbar}
\end{equation}
The optimal bound is then determined from
\[
\left[ \sum_l \left(\frac{C_l^{\phi g}}{\Delta C_l^{\phi
g}}\right)^2\right]^{1/2}\,.
\]

The redshift distribution Eq.~(\ref{redshiftdistri}) was used and the galaxy bias was set equal to $b_g=1$. The instrumental properties for Planck are given in Refs.~\cite{Miller:2008zi,Planck:2006aa}. We show the signal-to-noise ratio for Planck with NVSS as a function of $l_{\text{max}}$ in Fig.~\ref{planck-nvss}. We find that the highest signal-to-noise ratio, i.e.\ $15\sigma$, saturates at $l_{\text{max}}=2000$. Since the lensing-galaxy cross-power spectrum scales as $C_l^{\phi g}\propto b_g$ as illustrated by Eq. (\ref{clphig}), the amplitude of this cross-power spectra is degenerate with the galaxy bias and the signal-to-noise for the cross-power spectrum can also serve as a prediction of the detection significance for the galaxy bias. Thus, we see Planck can detect $b_g$ with high precision which will lead to a better understanding of the correlation between the baryonic matter distribution and the dark matter distribution.

\section{Conclusion}
\label{sec:con}

We have calculated the lensing-galaxy cross-power spectra using WMAP and NVSS and the full covariance matrix to filter the data sets. Specifically, we performed a thorough analysis of WMAP-1, WMAP-3, WMAP-5 and WMAP-7 raw DAs. The cross correlations between WMAP-5, -7's 8 DAs (2Q-bands$+$2V-bands$+$4W-bands) with NVSS clearly and firmly show signals at $>3\sigma$ level. We took the effects of gradient stripes into account for the NVSS data, and determined the significance without and with gradient stripes removed. The major effects caused by the stripes can be seen from the first bin of either the galaxy auto-power spectrum or the lensing-galaxy cross-power spectrum; the first bin decreases if the gradient stripes are marginalized over. However, gradient stripes do not affect the lensing-galaxy correlation (compare Refs.~\cite{Schiavon:2012fc,Nolta:2003uy,MNR:MNR9764}). We have explicitly shown all these results in Tables \ref{lensingtable1} and \ref{lensingtable2}. In these two tables, column $\mathcal{C}^a$ are the results without the gradient stripes removed and column $\mathcal{C}^b$ are our main results with the gradient stripes removed. In order to validate the lensing-galaxy cross correlations, we produced a NVSS galaxy map in equatorial coordinates directly from the NVSS catalog and cross correlated it with the WMAP DA which is in galactic coordinates, we find that all the lensing-galaxy cross correlation amplitudes are negligible.

We investigated the impact of different NVSS pixelization resolutions and found no effect. We compared the sensitivities of the estimators both with the full and diagonal covariance matrix and found that the former more effectively reduces the variance, which is mainly caused by the sky-cut and the inhomogeneous instrumental noise. However, the former scheme involves the inversion of a large matrix which is computationally challenging.

We predicted the detection significance for the lensing-galaxy cross correlation or the galaxy bias for the upcoming Planck data with NVSS and found the detection significance will be improved by a factor of 5.

The minimum variance of the estimator assumes that the CMB and galaxy overdensity modes are Gaussian. However, if the CMB contains gravitational lensing, the bispectrum $\langle TTg\rangle$ is not zero; it induces an additional variance as indicated in Eq.~(\ref{NG}). We analytically and numerically confirm that this variance is actually not noticeable for WMAP and NVSS as pointed out in Ref.~\cite{Smith:2011we}. Furthermore, being aware of the potential non-Gaussian shape of the probability distribution function (PDF)~\cite{Smith:2012ty}, we specifically investigate the PDF of the cross-power spectrum amplitude $\mathcal {C}$ in terms of Kolmogorov-Smirnov test, the skewness and the kurtosis. The diagnostic tests are shown in Table~\ref{kstable}. All the PDFs pass the Kolmogorov-Smirnov tests. All the PDFs are consistent with being Gaussian-distributed (Figs.~\ref{like1} and \ref{like2}).

The lensing-galaxy cross correlations effectively link the early universe to the late universe and the CMB is served as a back light casting the dark cosmic web (which is formed by the dark matter) throughout the major expansion history of the universe. The gravitational lensing is a powerful tool to decode the information of dark matter distribution from the CMB and the lensing-galaxy cross-correlations further unveil the relationship between baryonic matter and dark matter.

\begin{figure*}
\includegraphics[bb=50 195 465 480,width=8.5cm]{23.eps}
\includegraphics[bb=50 195 465 480,width=8.5cm]{24.eps}
\includegraphics[bb=50 195 465 480,width=8.5cm]{25.eps}
\includegraphics[bb=50 195 465 480,width=8.5cm]{26.eps}
\caption{({\bf Set-3}) Probability distribution function
for the lensing-galaxy cross correlation. The likelihood functions are normalized to 1. From 1000 simulations, a set of $\{\mathcal{C}\}$ is generated for each one of the subfigures, then by counting the frequency of $\mathcal{C}$ within a bin, a step-like function (red) is plotted. For comparison, Gaussian likelihood (green) is plotted using the mean and the variance of the set $\{\mathcal{C}\}$. } \label{like1}
\end{figure*}
\begin{figure*}
\includegraphics[bb=50 195 465 480,width=8.5cm]{27.eps}
\includegraphics[bb=50 195 465 480,width=8.5cm]{28.eps}
\includegraphics[bb=50 195 465 480,width=8.5cm]{29.eps}
\includegraphics[bb=50 195 465 480,width=8.5cm]{30.eps}
\caption{({\bf Set-7}) Probability distribution function
for the lensing-galaxy cross correlation. See Fig.~\ref{like1} for detailed descriptions. } \label{like2}
\end{figure*}

\begin{acknowledgments}
We would like to acknowledge helpful discussions with Sudeep Das, Duncan Hanson, Christian Reichardt, Meir Shimon, and Amit Yadav. We acknowledge the use of {\sc CAMB},  {\sc HEALPix}\footnotemark[1]~\cite{healpixpaper}, and  {\sc LAPACK}\ software packages and the LAMBDA archive. The computational resources required for this work were accessed via the GlideinWMS~\cite{Sfiligoi2009} on the Open Science Grid~\cite{Pordes2008}. We are indebted to Frank Wuerthwein, Igor Sfiligoi, Terrence Martin, and Robert Konecny for their insight and support.
\end{acknowledgments}

\footnotetext[1]{\url{http://healpix.jpl.nasa.gov/}}


\section{Introduction}

In cataclysmic variable stars (CVs), a white dwarf primary star accretes
matter by way of Roche lobe overflow from a binary companion,
which resembles a main-sequence star.   The variety of
CV behaviors leads to a complicated taxonomy \citep{warner95}.
Many CVs undergo
dwarf nova outbursts, thought to be caused by
accretion disk instabilities which greatly increase the rate at which
matter moves inward in the disk.  Other CVs, called novalike
variables, remain in a bright state for years at a time.
In still others, called AM Her stars or polars, the matter
that is transferred becomes entrained in a strong white-dwarf
magnetic field, and is funneled directly onto the
white dwarf's magnetic pole.

The main driver of CV evolution is thought to be a gradual
loss of orbital angular momentum. This causes the Roche critical
lobe of the secondary star to shrink, leading to a shortening
of the orbital period $P_{\rm orb}$, and driving mass
transfer on long (Gyr) timescales.  The
mechanisms by which angular momentum is lost are not fully
understood.  It is often supposed that magnetic
braking of the secondary star predominates at longer
periods ($>$ 3 hr), and that magnetic braking
becomes inefficient at short period, so that
gravitational radiation predominates.
Around $P_{\rm orb}$ = 70 min, the
secondary becomes degenerate and its radius begins to {\it increase}
with mass, leading to a slow {\it increase} in the orbital
radius.  This turnaround is often called the {\it period bounce},
even though it is thought to take place very slowly.

The histogram of CV orbital periods shows a significant dip
at roughly 2 hr $< P_{\rm orb} <$ 3 hr, known as the {\it gap}
\citep{kolbstillthere}.  This is often explained as follows.
As the secondary loses mass, its thermal timescale increases
to become comparable to the time for the orbit to evolve, with the
result that the secondary exceeds its equilibrium radius.
At about three hours, the secondary becomes fully convective,
reducing the efficiency of magnetic braking.  As the orbital
evolution slows, the secondary detaches from its Roche lobe,
shutting down mass transfer. The detached system continues
to evolve to shorter periods, crossing the gap and eventually
re-establishing contact with the Roche lobe near
$P_{\rm orb} = 2$ hr.  While the mechanism by which this
happens is somewhat speculative, \citet{kniggedonor} shows
that a discontinuity in the secondary stars' radii occurs
across the gap.

In a steady state, the number of stars with a given $P_{\rm orb}$
should be inversely proportional (roughly) to $\dot P$, the
rate at which the period changes.  If
$\dot P$ really is very slow at short periods, then there
should be a large population of short-period CVs, and the
gradual turnaround at the period bounce around 70 minutes should
lead to a `spike' in the distribution
\citep{pattlate,barker03,unveils}.

Efforts to confront theories such as this with observation
have often been frustrated, because the sample of known
CVs is incomplete in ways that are difficult to
quantify. The discovery
channels for CVs include the following:
(1) Dwarf nova outbursts are conspicuous -- they last for days or
weeks and typically have amplitudes of several magnitudes.
(2) Nearly all CVs have unusual colors compared to normal stars,
most conspicuously ultraviolet excesses arising from accretion
processes or (in some cases) the underlying white dwarf.
(3) With the exception of some novalike variables and dwarf novae near the
peak of outburst, nearly all CVs show emission lines, especially
in the Balmer sequence; these can be strong
enough to be noticed in surveys such as IPHAS \citep{withamiphas}.
(4) A great many other CVs have been discovered as optical counterparts
of X-ray sources.  (5) CVs at minimum light can be rather faint
($M > 10$), and a small handful have turned up in proper
motion surveys.

New, large samples of CVs with consistent selection
criteria are potentially useful for
clarifying issues such as the space density and orbital
period distribution of these objects.  Because the colors
of CVs overlap those of quasars, the Sloan Digital Sky Survey
(SDSS) turned up a large number of spectroscopically confirmed
CVs (\citealt{szkodyi,szkodyii,szkodyiii,szkodyiv,szkodyv,
szkodyvi,szkodyvii,szkodyviii}; hereafter referred to as
SzkodyI-VIII).
\citet{unveils} compiled orbital periods for 137
of these; with the SDSS sample, they were
finally able to discern the long-predicted period spike.

Like the SDSS, the Catalina Real Time Survey \citep[CRTS; ][]{drakecrtts}
has discovered many new CVs.  While the
SDSS CVs were originally selected by color --
they were chosen for spectroscopy largely because their
colors overlap those of quasars -- the CRTS
selects entirely by variability.  Briefly, the CRTS
surveys most of the accessible sky at Galactic latitudes
$|b| > 10^{\circ}$ and declination $\delta > -30^{\circ}$
every lunation, using a 0.7 m Schmidt telescope
in the Catalina mountains near Tucson, Arizona.
They search for variability using
a master catalog that reaches to $m \sim 22$.
Objects that show abrupt outbursts of $> 2$ mag amplitude
lasting less than few weeks
are classified as likely CVs, making the CRTS a
prolific source of dwarf novae in particular.

The CRTS maintains a catalog of ``Confirmed/Likely''
CVs on the World Wide Web
\footnote{available at http://nesssi.cacr.caltech.edu/catalina/AllCV.html}.
We downloaded this catalog on 2012 March 7, when it
contained 584 objects, and used this data set for the
present analysis.

The CRTS CV sample is of great interest because of its
depth and selection criteria, so it is a natural
choice for follow-up studies similar to
\citet{unveils}.  \citet{ww12} describe high-speed
photometry of 20 objects, mostly from the
CRTS, and give orbital periods for 15, including
two eclipsing dwarf novae and two superhumpers.
Only two of their sample had periods longer than the
2-3 hour period gap.  The period distribution of their
sample also showed the spike just above the turnaround.
In addition, they found that dwarf novae with more
recorded outbursts tended to have longer orbital
periods than those with fewer outbursts.

Here, we report on followup spectroscopy and photometry
of the CRTS CVs listed in Table \ref{tab:star_info}.
Most of the CRTS CVs are too faint for us to follow up,
so we selected this sample based largely on apparent
brightness.  Operational considerations, such as the
ease with which observations could be interleaved
other programs, the nearness of the target
to opposition at midnight, and (for radial velocity targets)
the strength and tractability of emission and absorption
spectra also entered into target selection.

We find spectroscopic periods for 8 objects, one of which
was independently measured by \citet{ww12}.  In
addition, we obtained exploratory spectra of
28 others, largely to assess their feasibility for
radial-velocity studies.  We also obtained
standardized magnitudes for 37 objects (17 of which
also have spectra).

{\it Plan of this paper.}
We describe our equipment and techniques in
Section \ref{sec:techniques}.
Table \ref{tab:spectro_summary} summarizes all of our
spectroscopy. In
Section \ref{sec:explore} we describe the spectra of the 28
objects for which we have only a only a quick spectrum.
The eight objects for which we have radial velocity time series are
discussed in Section \ref{sec:rvs};
Table \ref{tab:velocities} gives the velocities, and
Table \ref{tab:parameters} lists parameters for
sinusoidal fits to the velocity time series.
Table \ref{tab:photometry} gives magnitudes and colors
of objects for which we have standardized photometry.
Finally, in Section \ref{sec:statistics}
we consider how the CRTS CV samples overlaps with other
lists of CVs, discuss the apparent selection biases
of CRTS and the implications for the CV population
in general.

\section{Equipment and Techniques}
\label{sec:techniques}

All our data are from MDM Observatory on Kitt Peak,
Ariznoa.  Nearly all are from the 2.4m Hiltner telescope,
but a single photometric measurement was
kindly taken by J. Halpern at the 1.3m McGraw-Hill telescope.

The spectroscopic observing setups were as follows.

{\it Modular spectrograph.} For most of the spectra
we used the modular spectrograph and a 600 line mm$^{-1}$
grating.  The detector was either `Templeton', a
1024$^2$ thinned SITe CCD that gave 2 \AA\ pixel$^{-1}$
from 4600 to 6700 \AA , or `Nellie', a thick
2048$^2$ CCD that gave 1.7 \AA\ pixel$^{-1}$ from
4460 to 7770 \AA , with vignetting toward the ends
of the range.  The choice of detector was
dictated by availability during a series of controller upgrades.
With the modspec, targets are centered using the
image reflected from the polished slit jaws.  For this,
we used new (2010) slit viewing optics and a
self-contained Andor Ikon CCD camera unit; this has
greatly improved the acquisition of faint targets with
this instrument.

All of our radial-velocity time series were taken with
the modular spectograph.  For wavelength calibration we
obtained comparison lamps in twilight and used the
night-sky $\lambda$ 5577 line to track spectrograph
flexure during the night.  We also observed flux standards
in twilight when the sky was clear, and bright
O and B-type stars in order to correct approximately
for telluric absorption.

{\it Ohio State Multi-Object Spectrograph (OSMOS).}
This versatile instrument \citep{martini11} images the focal plane onto
a 4 k $\times$ 4 k CCD, at a scale of 0.273 arcsec
pixel$^{-1}$.  Filters can be placed in the parallel
beam of the reducing camera for direct imaging, or
volume-phase-holographic grism disperser
can be inserted for spectroscopy.  In spectroscopic
mode, one can insert slits in two different locations
in the focal plane, which yield different wavelength
coverage;
we used the `inner' slit, which gives coverage from
3960 to 6875 \AA\ at 0.75 \AA\ pixel$^{-1}$.  OSMOS
has no slit-viewing camera, so targets are acquired
by taking a direct image (without slit and disperser)
and then moving the telescope to align the target with
the known location of the slit.  We inserted
a $V$ filter for the acquisition exposures, and
therefore obtained a rough $V$ magnitude (without
a color transformation) for our targets.
Some humid weather
in 2011 September intermittently caused
condensation in the center of the detector window;
fortunately, the
important H$\alpha$ feature was outside the
affected region.

For most of the spectroscopic reduction we used
IRAF\footnote{
IRAF is distributed by the National Optical Astronomy
Observatory, which is operated by the Association of
Universities for Research in Astronomy (AURA) under
cooperative agreement with the National Science Foundation.}
routines.  To extract one-dimensional spectra from the
modspec images, we used a local implementation of the
optimal-extraction algorithm of \citet{horne}.

We computed synthetic $V$ magnitudes for our spectra
using the passband tabulated by \citet{bessell}.  Our
slit was usually 1.1 arcsec wide, which means that
an unknown fraction of the light was lost; also,
many of our spectra were taken through thin cloud.
Experience suggests that our synthetic magnitudes
are accurate to $\sim 0.2$ mag in good conditions;
in poor conditions, they can be considerably too faint.

We selected eight of our targets for time-series
radial velocity observations, with the aim of determining
orbital periods.  For these, we pushed observations to large
hour angle in order to avoid daily cycle count aliases.
To measure emission-line radial velocities we used the
convolution algorithms described by \citet{sy80}, in
which an antisymmetric function is convolved with the line
profile and the zero of the convolution is taken as the
line center.   We used methods described by \citet{shafter83} to
tune the convolution function for maximum signal.  We estimated
uncertainties in the velocities by propagating the
counting-statistics errors in the spectral channels; these
estimates do not include possible systematic effects.
For absorption-line velocities, we used the \citet{tondav}
convolution algorithms as implemented by \citet{kurtzmink}.
To search for periods we used a `residual-gram' method
described by \citet{tpst}.  Once we had established
a period, we fit the time series with sinusoids
of the form $v(t) = \gamma + K \sin[2 \pi(t - T_0) / P]$.
Uncertainties in the fit parameters were estimated from the
scatter around the best fit using the procedure described
by \citet{cash79}.

For direct imaging, we mostly used the `Nellie' CCD,
which gave 0.24 \AA\ per pixel.  This detector is
insensitive in the ultraviolet, so we used
$BVRI$ filters (although $R$ and
$B$ were sometimes skipped).  We derived photometric
transformations from observations of \citet{landolt92}
standard stars.  The scatter in the transformations
was  generally $< 0.05$ mag.  As noted earlier, we also
have some direct images from OSMOS.

\section{Exploratory Spectra}
\label{sec:explore}

Table \ref{tab:spectro_summary} shows that we
obtained only brief exposures of most of our
spectroscopic targets, generally taking one or
two exposures in a single visit.  Our main purpose
was to verify that the objects showed spectra typical of
CVs.  Figs.~\ref{fig:idspecs1} and \ref{fig:idspecs2} show
these exploratory spectra.  Of the 28 objects surveyed,
27 appear to be bona fide CVs, the exception being
CSS0350+35.  Here are brief descriptions.

{\it CSS0051+20.}  Our one spectrum has modest signal-to-noise,
but shows clear H$\alpha$ emission and confirms that this is a
CV.  The continuum has a shape that hints at an
M-dwarf contribution, but we cannot be sure this is real.

{\it CSS0208+37.}  The Balmer emission lines are strong and
appear double-peaked, suggesting a low orbital inclination.
A prominent blue continuum may be from an accretion disk, or
may be from a white dwarf photosphere; on the other hand, the
Balmer decrement appears extreme, with H$\beta$ much
stronger than H$\alpha$, suggesting that an
instrumental effect might be enhancing the blue end
of the spectrum.  Dwarf novae declining from outburst
can sometimes show blue continua, but our synthesized and
measured magnitudes
($V = 18.3$ and $V = 18.23$, respectively) are both
a bit fainter than the CRTS minimum magnitude of
17.9, ruling out this interpretation; the filter
photometry was obtained only two days before the
spectrum.

{\it CSS0328+28.}  This is an unremarkable dwarf
nova spectrum.

{\it CSS0350+35.} The spectrum is a good match to an
M1 dwarf, and shows only a narrow H$\alpha$ line
similar to a dMe star.   While this might be a
CV in which mass transfer has stopped, there is
no sign of a white dwarf in the spectrum, so
the spectrum is consistent with a flare star.
The light curve at the CRTS does not rule out
either a CV or flare star, so the flare star
classification appears likely.
This is, notably, the
{\it only} object in our spectroscopic sample
that appears not to be a CV.

{\it CSS0401+08.} Our spectrum is consistent with a
dwarf nova in outburst, with weak H$\alpha$
emission and weak, broad H$\beta$ absorption
(with a hint of a central reversal) on a blue
continuum.  The magnitude synthesized from the
spectrum ($V = 16.5$) is also well above the
18.5 mag minimum listed in the CRTS.

{\it CSS0411+23.} The emission lines are relatively
narrow (H$\alpha$ has a FWHM near 300 km s$^{-1}$),
suggesting a relatively low orbital inclination.
Weak helium emission is present at $\lambda \lambda
6678, 5876, 4921$, and 5015.

{\it CSS0440+02.}  The emission lines are broad
($\sim 1900$ km s$^{-1}$ FWHM) and just show
double peaks at our signal to noise ratio.  The
system is probably close to edge-on.

{\it CSS0447+09.} This shows evidence for a
K-type contribution, in the form of an
absorption dip near 5170 \AA , and absorption
in the Na D lines.  The signal-to-noise
ratio is not good enough to warrant a detailed
decomposition of the spectrum, but the absorption
features suggest that about half the light comes
from a mid-K star, probably in the range
K4 to K6.  This suggests that the
orbital period is likely to be $\sim 6$ h
or greater \citep{kniggedonor}, but a much shorter
period is possible if the secondary has lost
much of its mass \citep{thorstenseneiuma}.

{\it CSS0505+08.}  We detect H$\alpha$ at
rather low signal-to-noise ratio, confirming
the CV nature of the object.

{\it CSS0514+08.} Both H$\alpha$ and
H$\beta$ are detected in emission.  The
FWHM of H$\alpha$ is around 1000 km s$^{-1}$,
suggesting an intermediate orbital inclination.
No secondary star is seen at our signal-to-noise
ratio.

{\it CSS0518$-$02.} This object shows an
unusual spectrum, with a blue continuum and
a narrow ($\sim 300$ km s$^{-1}$), relatively
weak H$\alpha$ emission line.  The NaD lines
are detected in absorption, but the continuum
does not show any other convincing late-type
features, so the NaD may be interstellar.
The spectrum's synthetic magnitude
($V = 16.9$) suggests that the object was somewhat
brighter than minimum light ($m = 17.5$), but
not dramatically so.  The CRTS light curve shows
many outbursts of $\sim 2$ magnitudes.  This may
be a Z Cam-type dwarf nova, with a persistent
plateau state near $m = 17.5$, or could be
some other kind of novalike variable.

{\it CSS0545+02.} The red Digital Sky Survey
and the SDSS finding chart (linked from the
CRTS site) show a nebulosity around this object,
strongest along a northeast-southwest axis,
extending to a radius of about 20 arcsec from
the object.  In our spectrum, the apparent
sharp absorption features at H$\alpha$
(6563 \AA), [NII] ($\lambda\lambda$ 6548 and
6584), and [SII] ($\lambda\lambda$ 6717 and
6730) are strong nebular emission features in the sky
that have been over-subtracted.  The long-slit
spectrum shows extended weak H$\alpha$, H$\beta$,
[NII] and [SII] emission along 2.6-arcmin
length of the slit, and [OIII] $\lambda\lambda$ 4959 and
5007 emission extending to $\sim 18$ arcsec
from the star.
The stellar spectrum shows H$\alpha$
emission with a width of $\sim 2000$ km s$^{-1}$,
very much resembling a CV.  The CRTS light curve
hangs mostly at $m = 16.5$ or so, but the
SDSS has $g = 18.7$, in good agreement with our
$V = 18.75$ and $V = 18.9$ synthesized from
our spectrum; the nebula may be affecting the
CRTS measurement.  If this really is a CV, the
nebulosity is especially intriguing; it is also
possible that this is a planetary nebula central
star, in which case the broad H$\alpha$ might arise
from a stellar wind.  The CRTS light curves shows some
faint upper limits, so the system may eclipse.

{\it CSS0558+00.}  Our single, low
signal-to-noise spectrum shows weak,
broad H$\alpha$ on a blue continuum.  The
synthesized magnitude ($V = 16.6$) is well above
the minimum $m = 19.0$ found by CRTS.  It appears
we caught this in a relatively rare outburst;
the CRTS light curve shows only three outbursts,
even though the object is well-observed.  Our
filter photometry, obtained 9 days earlier than the
spectrum, shows the object much fainter, at
$V = 20.45$.  This is a dwarf nova.

{\it CSS0905+12.} This appears to be an
ordinary dwarf nova observed near minimum
light. H$\alpha$ has a FWHM of 1200 km s$^{-1}$,
suggesting an intermediate orbital inclination.

{\it CSS1055+09.}  The spectrum shows a
contribution from an M dwarf which is
presumably the secondary star.  Because of
our limited spectral coverage and signal-to-noise,
we can only roughly estimate the secondary's spectral type
to be M3 $\pm 2$.  The secondary appears to contribute at
least half the continuum light in the region of
our spectrum.  The easy visibility of the secondary
suggests that the orbital period lies above the 2-3 hour
gap.

{\it CSS1139+45.} Our spectrum is very poor, but
shows broad H$\alpha$ emission, confirming that this
is a CV.  The synthesized magnitude and spectrum
both indicate it was at minimum light.

{\it CSS1211$-$08.} This relatively bright
object shows strong Balmer and HeI lines.
H$\alpha$ has a FWHM of 800 km s$^{-1}$, indicating
a fairly low orbital inclination, and the spectrum
and synthesized magnitude indicate that the spectrum
was taken near minimum light.

{\it CSS1556$-$08.} Our spectrum has poor signal-to-noise,
but does show broad H$\alpha$ and H$\beta$ emission.
The synthesized
$V = 19.3$ is fainter than the CRTS minimum, but
conditions were partly cloudy at the time.

{\it CSS1616$-$18.}  The spectrum is typical
of dwarf novae in outburst, with narrow H$\alpha$
emission, broad H$\beta$ absorption, and a blue
continuum.  The synthesized $V = 15.9$, while the
CRTS lists variation between $m = 14.9$ and $m = 17.8$.

{\it CSS1649+03.} Broad H$\alpha$ emission is
clearly detected, but the poor signal-to-noise ratio
precludes further analysis.  Our photometric measurements
and the synthesized magnitudes both agree well with the
magnitude at minimum listed in the CRTS.

{\it CSS1720+18.}  Strong, broad H$\alpha$ and H$\beta$
confirm the CV status.  Our spectrum has
poor signal-to-noise, which is unsurprising given the
object's faintness.

{\it CSS1727+13.} H$\alpha$ has a FWHM of 1200
km s$^{-1}$, indicating an intermediate inclination.
The spectrum is typical of dwarf novae at minimum
light, and has a synthesized $V = 19.0$.  It was taken
only 3 days after photometric measurements showing
$V = 16.96$, indicating that the star had
faded quickly following an outburst.

{\it CSS1735+15.} Our spectrum was taken in partly
cloudy conditions, but does show H$\alpha$ in
emission.  In addition, NaD absorption is
present, and the K-star features near
$\lambda 5168$ are just detected.  The orbital
period is probably well longward of the 2-3
hr gap.  We think that the hump in the continuum
near 5300 \AA\ is probably an artefact.

{\it CSS1752+29.} The H$\alpha$ emission has a
modest equivalent width, and the continuum is
blue, suggesting that the system was in outburst.
Magnitudes from the OSMOS acquisition image and the
synthesized spectrum agree nicely at $V \sim 17.3$,
while the CRTS lists $m = 18.3$ for minimum light, so
the system was likely declining from outburst.
Again, the continuum hump near 5300 \AA\ is probably
an artefact.

{\it CSS2029$-$15 = SY Cap.}  \citet{kato09}
found superhumps
with $P_{\rm sh}$ = 0.063759(22) d = 91.8 min.
Using the $P_{\rm orb}$-$P_{\rm sh}$ relation derived
by \citet{unveils}, we predict $P_{\rm orb} = 89.8$ min.
Our spectrum shows narrow H$\alpha$  (FWHM $\sim$
400 km s$^{-1}$), indicating a low orbital inclination.

{\it CSS2059$-$09.} The H$\alpha$ line has a
FWHM of 1500 km s$^{-1}$, suggesting a fairly low
orbital inclination.
The CRTS light curve shows a
gradual brightening of the minimum magnitude over
time, with outbursts superposed.  Our spectrum and
acquisition image were taken in partly cloudy conditions.

{\it CSS2213+17.} The H$\alpha$ portion of our
spectrum was unaffected by a dewar condensation
problem, and shows a broad emission line, confirming
that this is a CV.

{\it CS2227+28.} This spectrum was also affected by
the condensation, but H$\alpha$ was again in the
clear portion, and is nicely detected with a somewhat
triangular profile, indicating a moderate orbital
inclination.

\section{Radial Velocity Studies}
\label{sec:rvs}

We obtained time series spectroscopy for eight
targets.  Figs.~\ref{fig:montage1} and \ref{fig:montage2}
show the average spectra, periodograms, and folded
velocity curves, and Table~\ref{tab:parameters}
gives the parameters of the sinusoidal fits.
We discuss the objects in order of RA.

\subsection{CSS0501+20}

The spectrum is typical of dwarf
novae at minimum light.  The lines are
single-peaked, and  H$\alpha$ has a
FHWM of 1100 km s$^{-1}$.  We adopt
$P_{\rm orb} = 107.7 \pm 0.2$ min. An
alternate choice of daily cycle count
gives 116.7 min, but the Monte Carlo test
described by \citet{tf85} assigns this
a probability below 1 per cent.  At this
orbital period, this is likely to be an
SU UMa star showing superhumps and superoutbursts.

\subsection{CSS0519+15}

The equivalent widths of the emission lines are
rather smaller than in most dwarf novae at minimum
light, suggesting that we caught the system somewhat
above minimum.  The lines are single-peaked and
the FWHM of H$\alpha$ is $\sim 800$ km s$^{-1}$,
indicating a fairly low inclination.  Our
$P_{\rm orb}, 122.3 (4)$ min, places this in the
period range of the SU UMa stars.

\subsection{CSS0647+49}

The spectrum shows a conspicuous contribution
from a late-type star.  Using spectra of
stars classified by \citet{keenan89}, we
find that the companion's spectral type to
be K4.5 $\pm$ 1 subclass, and that its contribution
to the spectrum is equivalent to $V = 18.0 \pm 0.4$.
In Fig.~\ref{fig:montage1}, the upper spectral
trace shows the average spectrum, and the lower
shows the result of subtracting a scaled K4 star from
the average.
Most of our observations are from
2011 March, but we also have velocities from January
and September.  An unambiguous cycle count over the
whole interval yields $P_{\rm orb} = 8.9160 \pm 0.0005$
hr.   The emission- and absorption-line velocities
are both modulated; the emission-line modulation
is $0.510 \pm 0.007$ cycles out of phase from the
absorption, consistent with the half-cycle offset
expected if the emission lines trace the white dwarf
motion.  If we {\it assume} that this is the case, then
the mass ratio (secondary to white dwarf) is
$0.73 \pm 0.04$.  With this mass ratio and a
typical white dwarf mass of $M_{\rm wd} = 0.8$ M$_\odot$,
the orbital inclination $i$ would be near 35$^\circ$.
While this is only meant to be illustrative, the
rather low secondary velocity amplitude ($K_2 = 117 \pm 4$
km s$^{-1}$) implies an orbital inclination
low enough that eclipses are very unlikely.

We estimate a distance using the secondary's
contribution as follows.  If we fix $P_{\rm orb}$ at its
measured value and assume that the secondary
star fills its Roche critical lobe,
then the secondary's radius,
$R_2$ depends almost entirely on its mass $M_2$;
furthermore, the dependence of $R_2$ on $M_2$ is weak.
Evolutionary calculations by \citet{bk00} suggest that a
K4 star in an 9-hour orbit should have
$M_2 = 0.7 \pm 0.2$ M$_\odot$.  As a check,
taking our mass ratio at face value and assuming
a broadly typical $M_1 = 0.8$ M$_\odot$ for the
white dwarf gives $M_2 = 0.58$ M$_\odot$, in reasonable
agreement.   With this mass range, the approximation
given in eqn.~1 of
\citet{beuermann2ndry} constrains
$R_2$ to be $0.85 \pm 0.10$ R$_\odot$.
From data in \citet{beuermann99}, we estimate the
surface brightness of a K4.5 dwarf to be such that
it would have $M_V = 6.3 \pm 0.5$ if it had
$R = 1$ R$_\odot$.  Scaling this to the secondary's
radius yields $M_V = 6.6 \pm 0.6$.  The secondary's
synthetic $V$ magnitude is 18.0, but this is probably
a little too faint, because the synthetic magnitude
of the average spectrum ($V = 17.2$) is
fainter than the $V = 16.84$ we find from
the more accurate filter photometry.  Discrepancies
this large are to be expected (Section
\ref{sec:techniques}).
Correcting for these losses gives
$V \sim 17.6$ for the secondary contribution.
At this celestial location
($l = 166.19^\circ, b = 19.79^\circ$),
\citet{schlegel98} estimate the total Galactic
extinction to be $E(B-V) = 0.11$, which taking
$R = 3.3$ makes $A_V = 0.36$, assuming the star
lies outside the dust layer.  Putting all this
together yields $(m - M)_0 = 10.9 \pm 0.7$,
or a distance $d = 1300 (+500,-400)$ pc.
Notice that we did not assume that the secondary
star follows a main-sequence mass-radius relation,
but rather combined the Roche size constraint with the
surface brightness.

\subsection{CSS0814$-$00}

This is an SU UMa-type dwarf nova, and was observed in superoutburst
by \citet{kato09}, who detected superhumps.  The superhump
period $P_{\rm sh}$ was not determined cleanly, but
appeared to be  0.0763 d.  Photometry obtained
by \citet{ww10} gave $P_{\rm orb} = 1.796$ h, or
0.0748 d, in rough accordance with expectations based
on the superhump period.  Our spectroscopic period is
essentially identical, at 0.07485(5) d.  The emission
lines are barely double-peaked, suggesting that the
inclination is not small, but \citet{ww10} make no
mention of an eclipse.  As far as we know, this is
the only system studied here which has a period
determination in the literature.

\subsection{CSS0902$-$11}

Like CSS0647+49, this object also has a strong contribution
from the secondary star.  By comparing the spectrum
with stars classified by \citet{keenan89},
we estimate the secondary's spectral type to be K7 $\pm$ 1
subtype, and find that the secondary's contribution is
nominally equivalent to $V = 19.1$.   However, many of the
exposures used in the average spectrum were taken in
partly cloudy conditions and mediocre seeing.
Our average spectrum has synthesized $V = 18.4$, but
our best exposures have $V = 17.8$.   The CRTS light
curve shows that the source is fairly steady when not
in outburst, at
$m = 17.5$, consistent with our best exposures, so
it appears that the mean spectrum from which the secondary
magnitude was derived is about 0.6 mag too faint.
We therefore adopt $V = 19.1 - 0.6 = 18.5$ for the secondary
contribution.

The emission-line radial velocities
did not yield a period, but the absorption spectrum
showed an unambiguous modulation at $6.62 \pm 0.01$ hr.
The relatively small velocity amplitude of the secondary
($K_2 = 100 \pm 6$ km s$^{-1}$) constrains the inclination
to be fairly low for any realistic white dwarf mass, so
eclipses are not expected.

We can once again estimate a distance using the secondary's
contribution, following the same procedure we used
for CSS0647+49.  Using \citet{bk00} as a guide, we estimate
$M_2 = 0.55 \pm 0.15$ M$_\odot$, from which we infer
$R_2 = 0.68 \pm 0.08$ R$_\odot$ at this period.
The surface brightness for K7 $\pm$ 1 is equivalent to
$M_V = 7.25 \pm 0.5$ for a 1 R$_\odot$ star \citep{beuermann2ndry}.
At this location ($l = 239.86^\circ,
b = 22.33^\circ$) the \citet{schlegel98} extinction map
gives $E(B-V) = 0.05$ mag.  Combining this with the
apparent magnitude of the secondary gives
$d = 1100 (+350, -260)$ pc.

\subsection{CSS0912$-$03}

The emission lines are notably broad -- the FWHM of H$\alpha$
is nearly 1700 km s$^{-1}$ -- and the lines show
incipient double peaks, so the inclination is likely to
be fairly high.  Even so, the radial velocity
amplitude is modest.  We detect a modulation in the
emission line velocity, but with the  available
the choice of daily cycle count is ambiguous, so
we give two possible periods in Table \ref{tab:parameters}.
Both possible periods are well below the 2-3 hour gap, so
it is likely that this will prove to be an SU UMa-type
dwarf nova with superoutbursts and superhumps.

\subsection{CSS1706+14}

The spectrum is typical of dwarf novae at minimum
light.  We obtained radial velocities on a single
night in 2011 June, after which the star went into
outburst, suppressing the H$\alpha$ emission and
ending the measurements.  In 2012 May we
obtained velocities on three consecutive nights.
From the combined
data we adopt a period near 0.0582 d, but
periods near 0.111 d and 0.0552 d (the latter being
a daily cycle count alias of our adopted period) are
not completely ruled out.  The $\sim$ 345 d gap in the
time series created fine-scale ringing in the
periodogram.  To derive the period uncertainty in Table
\ref{tab:parameters}, we shifted the 2012 May
data back in time by an integer number of periods,
removing the gap.  The exact period used in this
artificial shift has essentially no effect on the
resulting period uncertainty.

\subsection{CSS1729+22}

The spectrum shows weak M-dwarf features, in particular
the extra continuum around 5950-6170 \AA , and the
band head near $\lambda 6180$.  The features are too
weak to derive good constraints on the spectral type,
but an M1.5 dwarf contributing around 25 per cent of the
flux at 6500 \AA\ gives a reasonably good match.  As
one might expect from the spectrum, the period is
relatively long; the best-fitting $P_{\rm orb}$ is
7.12(3) hr, or 3.37 cycle d$^{-1}$; however, we cannot
entirely rule out an alias at 4.45 cycle d$^{-1}$, or 5.40 hr.

\section{The CRTS Sample and the Cataclysmic Variable Population}
\label{sec:statistics}

As noted earlier, the CRTS CV sample is
of great interest in characterizing the CV population.
In this section, we consider what the CRTS sample can
tell us about the completeness of the available
CV sample.

The number of non-CVs included in the CRTS sample
appears to be small.  In the sample of 36 objects
for which we obtained spectra, we found only a single
apparent non-CV.  The fainter objects in the CRTS
sample were beyond our magnitude limit,
and it is possible that the fainter end of the sample
includes a greater fraction of interlopers.
However, the selection criterion
-- outbursts of more than 2 mag --
should be fairly robust even for faint objects.
We assume, then,
that essentially all CRTS CVs are real CVs.

\subsection{Other Samples Used for Comparison}

We compared the CRTS list to several other samples of
CVs, which we describe here.

{\it The SDSS Sample.} SzkodyI-VIII list 286 CVs.
Although the SDSS Data Release 8 covers some
areas at low galactic latitude that CRTS does not \citep{sdssviii},
all of the SDSS CVs
lie within the nominal sky coverage of the CRTS, so
for practical purposes the SDSS coverage is entirely
contained within the CRTS coverage.  SzkodyI-VIII
do not tabulate the subtypes of their objects, though
they do give limited information on this.  To enable
more detailed comparisons, we classified the objects
in SzkodyI-VIII, primarily on the basis of their
spectra, supplemented by the information in the
text of the papers.  Our classification scheme was
as follows:

\begin{description}

\item[DN]  This class was used for objects known to be
dwarf novae, and for objects whose spectra resemble those
of dwarf novae.
The spectra classified this way tended to have strong,
broad Balmer emission
(with H$\alpha$ equivalent width usually greater than 30 \AA\ ),
relatively flat disk continua (in $F_\lambda$ vs.~$\lambda$), and
weak or absent HeII $\lambda$4686.
Objects showing blue continua and white-dwarf
absorption wings around H$\beta$ were classified as
DN-W; if a K or M-type secondary was present, we
classified the object as DN-2.  In some cases the SDSS
spectrum shows the object in outburst.  Dwarf novae in
outburst can be difficult to distinguish from novalike
variables, but in these cases the flux level in the
spectrum will usually be much greater than expected
based on the imaging data.

\item[NL] This class included spectra showing blue
continua, without white-dwarf absorption,  and relatively
weak emission lines, or
stronger emission lines and substantial HeII $\lambda$4686
(typically half the strength of H$\beta$ in those cases).
The Balmer absorption wings in a novalike variable can
superficially resemble white-dwarf absorption, but with
experience the distinctive white-dwarf line profile can
usually be distinguished from the disk absorption lines
seen in UX-UMa type variables.

\item[AM] These showed HeII $\lambda$4686 similar in
strength to H$\beta$, or other evidence of a
strong magnetic field such as cyclotron humps.

\item[NCV] We assigned 14 objects to this ``Non-CV" class.
This heterogeneous group include objects whose spectra resembled
reflection-effect white dwarf systems, subdwarf B stars, and
chromospherically active M dwarfs.  One,
SDSS J1023+00, has proven to be a binary
containing a millisecond radio pulsar \citep{archibald,wang,
thorj1023}.

\end{description}

Removing 14 apparent non-CVs from the
SDSS sample leaves us with 272 objects.
Because of the limited information --
especially the lack of long-term light curves -- the
classifications should be considered somewhat rough.
They are given in Table \ref{tab:sdssclasses}.

{\it The Ritter-Kolb Catalog.} \citet{ritterkolb}
have maintained a catalog (hereafter RKcat) of
CVs and related object with known or suspected orbital
periods.  Some of these were discovered in the CRTS,
but entries that do not match with CRTS objects
clearly were not.  In our comparisons, we
used the cataclysmic
binary list, `cbcat', from version 7.16 of the
Ritter-Kolb catalog (hereafter RKcat)\footnote{available at
http://www.mpa-garching.mpg.de/RKcat/}, which
contains 926 objects, of which 582 are in the
CRTS survey area.
RKcat provides subclassifications similar to the ones
we invented for the SDSS sample.

{\it The General Catalog of Variable Stars.}
We downloaded the General Catalog of Variable Stars
\citep{gcvs}; the 2012 January version contains 45678
entries, of which 753 are classified as some kind of
CV (types N, NA, NB, NC, NL, UG, UGSS, UGSU, or AM).
Of these, 287 are in the CRTS survey area.

\subsection{Comparison with Other Samples.}

Table \ref{tab:samplestats}
shows the number of cross-matches, and non-matches
between the CRTS sample and the lists detailed in the
previous section.

There are only 44 cross-matches between the
CRTS and the SDSS CV samples.
Fig.~\ref{fig:sdsscrtshist} shows histograms of the
SDSS sample and the CRTS objects that lie in the
SDSS footprint.  Note the following:

\begin{enumerate}

\item
At minimum light, most of
the CRTS objects are too faint to have been identified
as CVs by SDSS.  The fainter CRTS objects are for the most part
detected in the SDSS imaging data (provided they were
in the SDSS footprint),
but were too faint to be selected for spectra.

\item Among the brighter CRTS objects, somewhat more than
half are not in the SDSS sample.

\item The top panel shows that
only a rather small fraction of SDSS objects are recovered
by CRTS.

\end{enumerate}

Why does CRTS miss so many CVs?  As
noted earlier, one expects CRTS to
preferentially select dwarf novae, and to
be less sensitive to other CV subtypes.
Consistent with this, 40 of the 44
cross-matches between SDSS and CRTS are objects
that we classified as `DN'.  The RKcat and
GCVS comparisons (Table \ref{tab:samplestats})
also show this tendency for CRTS to select dwarf novae;
of the 134 CRTS objects that are listed in SDSS,
RKcat, or GCVS, 123 of them are classified as
dwarf novae in at least one of these catalogs, and only
11 are other kinds of CV.  We therefore confirm that,
as one might expect, CVs that are not dwarf novae
are largely passed over by CRTS.

The objects recovered by CRTS are mostly dwarf novae,
but are the objects {\it not} recovered by CRTS
{\it not} dwarf novae?
In the latter part of the table, we give the numbers
of CVs that lie in the nominal CRTS survey area, but which
are {\it not} recovered by CRTS.
As expected, dwarf novae constitute a smaller portion of
the unrecovered objects; the aggregate dwarf nova fraction
is 348/602, rather than 123/134.

However, these figures also imply that over half the
{\it unrecovered} objects actually {\it are}
dwarf novae.  Somehow, 348 dwarf novae that we know of
in the CRTS survey area have slipped through its net.
How can we account for this?

Some of these non-recoveries (or `misses')
are to be expected, because dwarf novae from
one subclass -- the WZ Sge stars -- erupt very infrequently,
in some cases on timescales of decades or more.  An
even more extreme example, the star
GD 552 \citep{undana-gd552}, appears identical to a
short-period dwarf nova at minimum light, but it has
{\it no} observed outbursts.  Some fraction of the
WZ Sge stars will have been missed simply because they
have not outburst during the $\sim 5$ years that
CRTS has operated.

To explore the effect of outburst interval on CRTS selection,
we used data from RKcat, which
gives an average outburst interval (which they denote $T1$)
for some of the listed dwarf novae.
We arbitrarily chose 700 days, roughly half the time span
of the CRTS survey, as the dividing line between
short and long outburst intervals.  In addition, we assigned to the
long-interval group any
dwarf nova subclassified as `WZ'.
Many dwarf novae could not be assigned because
because $T1$ was not given (and they
were not classed as WZ), but sufficient information
existed to classify 137 dwarf novae from
the CRTS footprint; 95 of these were in the
short-interval group, and 42 in the long-interval group.
The objects that {\it were}
recovered in CRTS included 18 from the short-interval group, and
9 from the long-interval group;
among RK dwarf novae that were {\it not}
recovered by CRTS, there were 77 short-interval systems, and 33 long-interval
systems.  The fact that so many short-interval systems are
{\it not} recovered by CRTS argues strongly that long outburst intervals
are not the main reason CRTS is missing so many dwarf novae (although
it must account for some cases).  Indeed, the ratio of
short to long interval dwarf novae for CRTS/RK matches is
remarkably similar to the ratio in the group of RK dwarf novae
that are not matched to CRTS.

The main reason for the incompleteness must therefore lie elsewhere.
Dwarf nova outbursts will be missed if they occur between observations.
Perusal of CRTS light curves suggests that some parts of the
sky are covered rather infrequently.
A good number of dwarf nova outbursts must
therefore `slip through the cracks' in the coverage.  This
is exacerbated by the 2-magnitude criterion; not only must
the object be caught in outburst, but during that part of the
outburst where it is more than 2 magnitudes above minimum.
To see whether outburst amplitude might be an important
selection factor, we again turned to RKcat, this time
finding the outburst amplitude $\Delta m$.  For objects
that did not have superoutburst magnitudes listed, we
let $\Delta m = mag1 - mag3$ (in their notation), but
if superoutburst magnitudes were available, we used
$mag1 - mag4$.  Fig.~\ref{fig:outbcdf} shows the cumulative
distribution functions of $\Delta m$ for the recovered and
unrecovered dwarf novae.  As expected, CRTS is
nearly blind to objects with $\Delta m < 2$; more
interestingly, there is a significant bias against
smaller outburst amplitudes extending all the way
up to $\Delta m = 6$.  This effect probably arises
because, in any given snapshot, a large-amplitude object
will have a greater likelihood of being caught at
$\Delta m > 2$ than a small-amplitude object, which
would have to be fairly close to its peak brightness
in order to exceed the survey threshold.

It seems, then, that the CRTS survey is biased toward
large outburst amplitudes.  How might this affect other
quantities?  In Fig.~\ref{fig:amplvsper}
we plot $\Delta m$ against $P_{\rm orb}$,
for those dwarf novae in the CRTS footprint that have
both quantities tabulated in RKcat.  While there
is a great deal of scatter at any given $P_{\rm orb}$,
there is a clear trend for short-period dwarf novae
to have greater outburst amplitude, and hence a greater
likelihood of being discovered by CRTS.  The CRTS sample
should therefore be biased toward shorter orbital periods.
This must contribute at some level to the preponderance of
short periods found in this paper and by \citet{ww12}.

The fact that so many of the CRTS CV sample are new discoveries,
and its continuing effectiveness in finding new ones, both show
that a great many CVs remain undiscovered.  Future
synoptic surveys with faster cadence, and
less-stringent variability criteria should discover many
more CVs.

\section{Summary}

We obtained spectra of 36 CRTS CV candidates, and confirmed
that all save one appear to be {\it bona fide} CVs.  For eight
of the objects we obtained spectroscopic periods, and found that
three of them had $P_{\rm orb}$ longward of the 2-3 hour gap.
In addition, we examined the overlap between the CRTS CV sample
and others previously-existing samples.  Most CRTS CVs are new
discoveries, but CRTS has not recovered even a majority of the known
dwarf novae in its footprint.   This suggests that a great many
CVs remain undiscovered.  Analysis of the recovered and
unrecovered samples shows that the CRTS sample is biased toward
large outburst amplitudes, which in turn biases it toward
shorter orbital periods.

\acknowledgments
We gratefully acknowledge support from NSF grants
AST-0708810 and AST-1008217, and thank Dartmouth undergraduates
Jason Spellmire and Erin Dauson for conscientious and
cheerful assistance at the telescope.

\clearpage