%%
%% Beginning of file 'sample.tex'
%%
%% Modified 2005 December 5
%%
%% This is a sample manuscript marked up using the
%% AASTeX v5.x LaTeX 2e macros.

%American Astronomical Society (AAS)	

%% The first piece of markup in an AASTeX v5.x document
%% is the \documentclass command. LaTeX will ignore
%% any data that comes before this command.

%% The command below calls the preprint style
%% which will produce a one-column, single-spaced document.
%% Examples of commands for other substyles follow. Use
%% whichever is most appropriate for your purposes.
%%
%%\documentclass[12pt,preprint]{aastex}
%% manuscript produces a one-column, double-spaced document:
%% preprint2 produces a double-column, single-spaced document:

%\documentclass[preprint2]{aastex} %Use for ApJL submission
\documentclass[apjl]{emulateapj} %Use for ApJ submission or ApJL spacing estimate

% options for emulateapj package
%[iop] - apj style but adjusted due to change in publisher from chicago to iop
%[chicago] - (default) typeset as was done in the University of Chicago Press
%[twocolumn] - (default) two-column mode
%[onecolumn] - main text in one-column mode
%[apj]       - typeset as for main journal
%[apjl]      - (default) typeset as for ApJ Letters 
%[tighten]   - some adjustments to approximate grid typesetting
%[numberedappendix]   - number appendix sections as A, B, etc
%[appendixfloats]  - use separate numbering for floats within appendix
%[twocolappendix]  - make appendix in two-col mode in a two-col paper
%[revtex4]   - force using revtex4

\usepackage{amsmath}  % needed for \tfrac, \bmatrix, etc.
\usepackage{amsfonts} % needed for bold Greek, Fraktur, and blackboard bold
\usepackage{graphicx} % needed for figures
\usepackage{wasysym}
\usepackage{calc}
\usepackage{amssymb}
\usepackage{relsize}
%\usepackage{float}
% \usepackage{pxfonts} %for upright greek alphabets
\usepackage{natbib} %for bibtex

% \usepackage[percent]{overpic} %to overlap images

%% Sometimes a paper's abstract is too long to fit on the
%% title page in preprint2 mode. When that is the case,
%% use the longabstract style option.

%% \documentclass[preprint2,longabstract]{aastex}

%% If you want to create your own macros, you can do so
%% using \newcommand. Your macros should appear before
%% the \begin{document} command.
%%
%% If you are submitting to a journal that translates manuscripts
%% into SGML, you need to follow certain guidelines when preparing
%% your macros. See the AASTeX v5.x Author Guide
%% for information.

%% You can insert a short comment on the title page using the command below.
%\slugcomment{The Astrophysical Journal, 680:L87-L90, 2014 April 11}

%% If you wish, you may supply running head information, although
%% this information may be modified by the editorial offices.
%% The left head contains a list of authors,
%% usually a maximum of three (otherwise use et al.).  The right
%% head is a modified title of up to roughly 44 characters.
%% Running heads will not print in the manuscript style.

%short title length maximum of 44 characters ; 44 now
%\shorttitle{Measuring Temporal Photon Bunching in Black Bodies}
%\shortauthors{P. K. Tan and C. Kurtsiefer et al.}

%% This is the end of the preamble.  Indicate the beginning of the
%% paper itself with \begin{document}.

\begin{document}

%% LaTeX will automatically break titles if they run longer than
%% one line. However, you may use \\ to force a line break if
%% you desire.

\title{Measuring Temporal Photon Bunching in Blackbody Radiation}

%% Use \author, \affil, and the \and command to format
%% author and affiliation information.
%% Note that \email has replaced the old \authoremail command
%% from AASTeX v4.0. You can use \email to mark an email address
%% anywhere in the paper, not just in the front matter.
%% As in the title, use \\ to force line breaks.

\author{P. K. Tan\altaffilmark{1,2,*}}
\author{G. H. Yeo\altaffilmark{2}}
\author{H. S. Poh\altaffilmark{1}}
\author{A. H. Chan\altaffilmark{2}}
\author{C. Kurtsiefer\altaffilmark{1,2,+}}
\altaffiltext{1}{Center for Quantum Technologies, 3 Science Drive 2, Singapore 117543}
\altaffiltext{2}{Department of Physics, National University of Singapore,
  Singapore 117551}
\altaffiltext{*}{email: pengkian@physics.org}
\altaffiltext{+}{email: phyck@nus.edu.sg} 

%% Notice that each of these authors has alternate affiliations, which
%% are identified by the \altaffilmark after each name.  Specify alternate
%% affiliation information with \altaffiltext, with one command per each
%% affiliation.

%% Mark off your abstract in the ``abstract'' environment. In the manuscript
%% style, abstract will output a Received/Accepted line after the
%% title and affiliation information. No date will appear since the author
%% does not have this information. The dates will be filled in by the
%% editorial office after submission.

%abstract length maximum of 250 words, single paragraph ; 81 words now
\begin{abstract}
Light from thermal black body radiators such as stars exhibits photon bunching
behaviour at sufficiently short time-scales. However, with available detector
bandwidths, this bunching signal is difficult to be directly used for intensity
interferometry with sufficient statistics in astronomy. Here we present an 
experimental technique to increase the photon bunching signal in blackbody
radiation via spectral filtering of the light source. Our measurements reveal
strong temporal photon bunching in light from blackbody radiation, including
the Sun. Such filtering techniques may revive the interest in intensity
interferometry as a tool in astronomy. 
\end{abstract}

%% Keywords should appear after the \end{abstract} command. The uncommented
%% example has been keyed in ApJ style. See the instructions to authors
%% for the journal to which you are submitting your paper to determine
%% what keyword punctuation is appropriate.

%max of 6 keywords, refer to list for keywords allowed at
% http://aas.org/authors/astronomical-subject-keywords-update-august-2013
\keywords{instrumentation: interferometers --- radiation mechanisms: thermal --- stars: fundamental parameters --- Sun: fundamental parameters --- techniques: interferometric}

%% From the front matter, we move on to the body of the paper.
%% In the first two sections, notice the use of the natbib \citep
%% and \citet commands to identify citations.  The citations are
%% tied to the reference list via symbolic KEYs. The KEY corresponds
%% to the KEY in the \bibitem in the reference list below. We have
%% chosen the first three characters of the first author's name plus
%% the last two numeral of the year of publication as our KEY for
%% each reference.

%% Authors who wish to have the most important objects in their paper
%% linked in the electronic edition to a data center may do so by tagging
%% their objects with \objectname{} or \object{}.  Each macro takes the
%% object name as its required argument. The optional, square-bracket 
%% argument should be used in cases where the data center identification
%% differs from what is to be printed in the paper.  The text appearing 
%% in curly braces is what will appear in print in the published paper. 
%% If the object name is recognized by the data centers, it will be linked
%% in the electronic edition to the object data available at the data centers  
%%
%% Note that for sources with brackets in their names, e.g. [WEG2004] 14h-090,
%% the brackets must be escaped with backslashes when used in the first
%% square-bracket argument, for instance, \object[\[WEG2004\] 14h-090]{90}).
%%  Otherwise, LaTeX will issue an error. 

\section{Background}

Statistical intensity fluctuations from stationary thermal light sources
formed the basis of the landmark experiments by \citet{hbt:54, hbt:57}, 
where correlation measurements between a pair of
photodetection signals were used to assess both temporal
and spatial coherence properties of light fields. The explanation for these
fluctuations as photon bunching phenomenon, as compared to a Poissonian
distribution for uncorrelated detection processes was a
major contribution for the development of the field of quantum optics \citep{glauber:63} 

\citet{hbt:56, hbt:74} developed their observation into a spatial intensity
interferometer which allowed them to infer stellar angular radii
from the coherence properties of star light over relatively large distances of 188\,m. Since a key quality of any telescope system is its angular resolution, given by its real or virtual aperture, they managed to
carry out these measurements with an unprecedented accuracy.

While stellar amplitude interferometry has evolved to become a standard technique to
increase the effective angular resolution of radio telescope arrays, the
angular resolution of the Hanbury-Brown--Twiss approach is still unmatched by
conventional diffraction-limited optical telescopes. It is also very
challenging to reproduce with Michelson interferometry due to the required
sub-wavelength optical path length stability \citep{davis:99}. Such setups involve remote locations for
observatory sites to reduce atmospheric turbulence, as well as expensive and
complicated adaptive optics to augment the measurable sensitivity by
suppressing optical noise \citep{millour:08}. Yet, the comparatively simple
optical intensity interferometry technique has not really found widespread
adoption in the observation of celestial objects. 

\pagebreak

\section{The problem}
Temporal photon bunching, i.e., the tendency for photons to cluster together
on a short time-scale, has been consistently observed from narrow-band
quasi-thermal light sources such as low pressure glow discharge lamps
\citep{martienssen:64, morgan:66}, and 
laser beams \citep{asakura:70} which are passed through a rotating ground
glass to randomize the optical phase in time.

Photon bunching in time for stationary fields is quantified by a second order
correlation function $g^{(2)}(\tau)$, which is the probability of observing a
second photodetection event at a time $\tau$ after a first one, normalized to the
probability of an uncorrelated detection. Optical modes in thermal states are
the textbook example where $g^{(2)}=2$ \citep{glauber:63, mandel:95,
  salehteich:07}. Blackbody radiation would 
therefore be the obvious choice for observing photon bunching, if it were
not for the short coherence time associated with its spectral properties.

In a scenario where spatial coherence effects are neglected, the second order
temporal correlation function for light fields can be written as
\begin{equation}
g^{(2)}(\tau)=1+\left|\gamma(\tau)\right|^2\,,
\end{equation}
where $\gamma(\tau)$ is the complex degree of coherence, or normalized
self-coherence of the stationary light field with $\gamma(0)=1$, which itself is proportional to
the Fourier transform of the spectral density $S(\nu)$ of the light field
\citep{mandel:95}. For well-defined spectral densities and coherence
functions, a meaningful coherence time $\tau_c$ can be defined, which is
reciprocal to a characteristic width $\Delta \nu$ of the spectral distribution
$S(\nu)$. 

For blackbody light sources in the optical regime, the coherence time
$\tau_{c}$ given by Planck's law is on the order of $10^{-14}$\,s, which is much
shorter than the best timing resolution $\tau_{t}$ of a few
$10^{-11}$\,s of existing photodetectors. 
As shown by \citet{scarl:68}, the observable excess from uncorrelated
pair detections, $g^{(2)}(\tau = 0)-1$, scales as $\tau_{c}/\tau_{t}$
and is thus only significant when the coherence time of the
source is of the order of the detection timing response.

Therefore, empirical observation of temporal photon bunching from blackbody
light sources has been extremely difficult, with only recent breakthroughs
using two-photon absorption in semiconductors \citep{boitier:09}, and in
ghost imaging research using a narrowband Faraday anomalous dispersion optical
filter \citep{karmakar:12, wu:14}.

In this paper, we present a spectral filtering technique that significantly
increases the the observed photon bunching signature from a blackbody light
source, such that intensity interferometry techniques for
observing celestial objects may become attractive again.


\section{Experiment}
As with previous attempts to observe photon bunching from blackbody radiation
\citep{karmakar:12, wu:14}, we employ narrowband optical filtering of the
light source in order to increase 
the coherence time $\tau_c$, and use single photon detectors with low
timing uncertainties. The experimental setup is shown in
Fig.~\ref{fig:setup}.

\begin{figure}
\centering
%\includegraphics[width=\hsize]{apj_hbt.eps}
\includegraphics[width=\hsize]{figure1.eps}
\caption{Setup to determine the temporal correlation $g^{{(2)}}(\tau)$ for
  wideband thermal light. The initial spectrum gets filtered to a narrow
  optical bandwidth such that the temporal decay of the second order
  correlation function can be observed with conventional single photon
  detectors in a Hanbury-Brown--Twiss experiment. We employ a combination of
  a grating and a temperature-tuned etalon as a spectral filter, and ensure
  spatial coherence in the setup by using single mode optical fibers (SMF).}
\label{fig:setup}
\end{figure}

Firstly, light from different sources under investigation is coupled into a
single mode optical fiber to enforce spatial coherence, because the
photon bunching signature $g^{(2)}(\tau=0)-1\propto1/M$, where $M$ is the
number of transverse modes \citep{glauber:63}. 

The light is then collimated with an aspheric lens and sent through the first
spectral filter, which is either an interferometric bandpass filter with a full width at half maximum (FWHM) of 3\,nm centered at 546.1\,nm to match the
$6p7s\,^3S_1 - 6p6s\,^3P^0_2$ transition in a Mercury (Hg) discharge lamp, or a grating
monochromator for broadband continuous spectrum. While the interference filter
is sufficient to select the spectral lines of the different mercury isotopes,
the grating monochromator had to be used for the blackbody sources.

\begin{figure}
\centering
%\includegraphics[width=\hsize]{transmission.eps}
\includegraphics[width=\hsize]{figure2.eps}
\caption{(Top) Transmission profile of grating monochromator. The solid line
  is a fit to a Gaussian profile to measurements for different grating angles
  (symbols), resulting in a full width at half maximum (FWHM) of
  0.122$\pm$0.002\,nm. (Bottom) Transmission of Hg light near 546\,nm through the
  temperature-tuned etalon with 2\,GHz bandwidth. The hyperfine structure for
  different isotopes is partially resolved, as identified in
  \citep{sansonetti:10}.}
\label{fig:transmission}
\end{figure}

The core of the monochromator is a blazed reflective diffraction grating
with 1200\,lines/mm mounted on a rotation stage. Light from the
optical fiber (NA = 0.13) is collimated with an achromatic doublet ($f=30\,$mm),
illuminating about $10^4$ lines of the grating. The first diffraction order is
imaged with another achromat of ($f=30\,$mm) into another single mode optical
fiber. This arrangement acts as a tunable narrowband filter in the
visible regime, allowing to adjust for peak emissions
of blackbody light sources. A test of the transmission of a narrowband laser
source ($\lambda$ = 532\,nm) is shown in Fig.~\ref{fig:transmission} (top
part). The fit to a Gaussian transmission profile leads to an optical
bandwidth of 0.122$\pm$0.002\,nm (FWHM). We typically observe a transmission
of narrowband light from single mode fiber to single mode fiber of
$\approx$15\%.

In a second step, the light passed through a solid etalon (thickness
$d$=0.5\,mm) with flat surfaces. The etalon has a free spectral range of
$FSR=c/(2nd)=$ 205\,GHz, with a refractive index $n=1.46$ for the fused
silica material (Suprasil311), and the speed of light $c$ in vacuum. At
$\lambda$ = 546\,nm, the free spectral range corresponds to
$\Delta\lambda=\lambda^2/c\cdot FSR\approx$ 0.2\,nm, so that only one
transmission order should fall within the transmission width of the
monochromator.
Both surfaces of the etalon have dielectric coatings with a
reflectivity of $R\approx$ 97\% over a wavelength range from 390\,nm to
810\,nm, resulting in a transmission bandwidth (FWHM) of
\begin{equation}
  \Delta\nu \approx FSR\cdot{1-R\over\pi\sqrt{R}}\approx2\,\mathrm{GHz}\,.
\end{equation}

The central transmission frequency of the etalon is tuned by temperature, with
a scaling rate of -4.1\,GHz/K due to thermal expansion and temperature
dependence of the refractive index. A temperature stability of $\pm$5\,mK over the
measurement duration of several hours ensures that the etalon stays on
a specific wavelength. Temperature tuning allows the etalon to be always used
under normal incidence, which avoids frequency walk-off
\citep{green:80} and spatial distortion \citep{park:05}. To demonstrate the
resolving power of this device, we tuned the peak transmission across the
various spectral lines from a Hg gas discharge lamp at 546.1\,nm, revealing a
partially resolved hyperfine structure for different isotopes
% with non-zero nuclear spins
(see Fig.~\ref{fig:transmission}, bottom part).

A Glan-Taylor polarizer (PBS) following the etalon selects
a well-defined linear polarization mode for the photons, with the same
rationale of minimizing the number of modes ending up in the temporal
correlation setup. A half wave plate ($\lambda/2$) transforms the
polarization emerging from the grating/fiber combination such that there is
maximal transmission through the PBS. 

A second Glan-Taylor polarizer, rotated by $\approx45^\circ$  with respect to
the polarization direction of the first, distributes the light onto two
silicon avalanche photodetectors (APD) that provide detection timing signals
to a correlator. The relative orientation of the two polarizers allows for the
balancing of the photodetection rates to optimize the integration time, and
also helps to suppress  the optical crosstalk 
between the two APDs due to hot carrier recombination in the devices upon
detection which would otherwise introduce false coincidence events
\citep{kurtsiefer:01a}.

The temporal correlation was carried out with a digital sampling
oscilloscope (4\,GHz analog bandwidth, 40\,Gsamples/sec, LeCroy) that can
measure and histogram detection timing differences between two APD events with a
timing uncertainty below 10\,ps. 

\begin{figure}
%\centerline{\includegraphics[width=\hsize]{apd_combined.eps}}
\centerline{\includegraphics[width=\hsize]{figure3.eps}}
\caption{Comparison of timing jitter between two actively quenched thin avalanche
  photodiodes (PDF), and two passively quenched deep avalanche photodiodes
  (C90302S) by photopair detection from light generated by spontaneous
  parametric down conversion (SPDC) with a coherence time below 10\,ps.}
\label{fig:jitter}
\end{figure}

To assess the overall timing uncertainty with the detector/correlator combination,
we exposed two different avalanche photodetector sets to photon pairs
generated by spontaneuous parametric down conversion (SPDC) with a coherence time below
10\,ps. The resulting temporal correlations shown in
Fig.~\ref{fig:jitter} reflects the overall timing resolution of the system.
For APDs with a deep conversion zone (C30902S, Perkin-Elmer, passive
quenching), the coincidence distribution appears roughly Gaussian, with 1.2\,ns
FWHM. Since the photon bunching signature $g^{(2)}(\tau)$ of the light field
under investigation needs to be convoluted with this detector response, we
would observe a significant reduction of the photon bunching peak with
these detectors. A second set of APDs, based on thin conversion zones (type
PDF, Micro Photon Devices, active quenching), leads to a much better localized
coincidence distribution with 40\,ps FWHM, although a significant fraction of
the distribution shows a tail with timing differences up to 500\,ps.

\section{Results}

\begin{figure}
\centering
%\includegraphics[width=\hsize]{spectrum_hgarsun.eps}
\includegraphics[width=\hsize]{figure4.eps}
\caption{Spectrum of a Hg glow discharge lamp, an Ar arc lamp as a simulator of a stellar light source, and the Sun.}
\label{fig:spectrum}
\end{figure}

Measurements for $g^{{(2)}}(\tau)$ were performed on three light sources: a
low pressure Mercury (Hg) discharge lamp as a quasi-thermal light source, a
radio frequency-driven Argon (Ar) arc discharge lamp to produce a plasma with a
blackbody temperature $\approx$ 6000\,K as the thermal light source, and the
Sun as a celestial thermal light source. Their respective spectral distributions are shown in Fig.~\ref{fig:spectrum}. The optical filters were fixed to
maximize transmission for the spectral lines around 546.1\,nm in the Hg
spectrum. The resulting timing difference histograms are shown in
Fig.~\ref{fig:g2}.

\begin{figure}
\centering
%\includegraphics[width=\hsize]{g2_combined.eps}
\includegraphics[width=\hsize]{figure5.eps}
\caption{Intensity correlation function $g^{{(2)}}(\tau)$ from the three
  thermal light sources after spectral filtering. All light sources show
  a significant bunching at $\tau=0$, and a decay time that is
  well-resolved with the avalanche photodetectors. Error bars indicate
  uncertainties due to Poissonian counting statistics.}
\label{fig:g2}
\end{figure}

To ensure proper normalization of the resulting 
$g^{(2)}(\tau)$, we independently recorded pair events  with
a resolution of about 300\,ps with a time stamp unit that, unlike the
oscilloscope, does not exhibit a significant dead time, and found that
$g^{(2)}(|\tau|\gtrsim2\,\textrm{ns})\rightarrow1$. To extract a meaningful
set of model parameters, we assumed a Lorentzian line shape for the spectral
power distribution; this is mostly to reflect the approximate transmission
function of the etalon. A corresponding model function for the obtained pair
rates in the histograms is
\begin{equation}
N(\tau)=a+b\cdot e^{-\left|2\tau/\tau_c\right|}\,
\end{equation}
where $a,b$ capture the absolute event number and the degree of bunching,
respectively, while $\tau_c$ represents a coherence time. The $g^{(2)}$ scale
on the right side of diagrams in Fig.~\ref{fig:g2} refers to the absolute
number of events after division by fit parameter $a$. 

For the Hg emission line at 546\,nm from the low pressure glow discharge lamp,
coincidence data was obtained over an
integration time of 68 hours, with single detector rates of 115000\,s$^{-1}$.
The relatively large discrepancy between expected coincidences with the single
rates and the observed coincidences is due to a dead time of the oscilloscope
in this measurement.   
From the numerical fit, we find $g^{{(2)}}(\tau=0)=1.79\pm0.01$
and an exponential decay with a characteristic coherence time of $\tau_{{c}}$
= 0.436$\pm$0.006\,ns. The photon bunching signature comes very close to the
expected value of 2 for a thermal state, with the difference possibly due to
the residual timing uncertainty of our photodetectors.
While the Doppler-broadened line width for individual transitions in
Hg vapour at temperatures around 400\,K is on the order of 235\,MHz, the
contributions of different hyperfine transitions from different isotopes
result in a complex overall line shape of several GHz width, out of
which we select a reasonable fraction with the etalon. 

For the Ar arc discharge lamp, a wavelength of 540\,nm was
chosen to avoid contributions of a prominent peak at 546\,nm, possibly caused by mercury traces in the plasma. We find that
$g^{{(2)}}(0) =1.45\pm0.03$ and a shorter coherence time
$\tau_{c}=0.31\pm0.01$\,ns with a reduced $\chi^{2}_{r}=1.066$.
The wider spectrum of the cut-out blackbody radiation makes the decay time
shorter, so a peak reduction due to the detector response becomes more
prominent. Furthermore, residual transmission of the grating monochromator or
stray light into other transmission peaks of the etalon may contribute to
light that is not coherent with the main peak on observable time-scales. 

We finally collected sunlight into a single mode optical fiber with a
coelostat and transferred it into the lab. At the same working wavelength of
546.1\,nm, we find $g^{{(2)}}(0)=1.3\pm0.1$ and a coherence
time $\tau_{c}= 0.26\pm0.05$\,ns, with $\chi^{2}_{r}=1.0244$ over an
integration time of 45 minutes. At the moment, we cannot explain
% hopefully neither Ginglish nor Singlish now :)
the reduction of the visibility of the photon
bunching, but seeing-induced atmospheric fluctuations in optical path
length differences \citep{dravins:07a} may be a reason.

\section{Summary and Outlook}

We have successfully developed and demonstrated an experimental technique to
resolve and explicitly measure the temporal photon bunching from thermal
blackbody light sources, including both discrete spectrum and blackbody
spectrum such as the Sun. We observed a photon bunching signature, i.e., the
excess of $g^{(2)}(0)$ above 1, that exceeds the recently reported records
of $g^{(2)}(0)=1.03$ \citep{karmakar:12} and $g^{(2)}(0)=1.04$ \citep{wu:14}
by about an order of magnitude. 

Compared to photon bunching observations using two-photon absorption
\citep{boitier:09}, our method would work at astronomical
intensity levels, and not be constrained to cases with the simultaneous arrivals of two photons at
the same detector.

We believe that with better timing control of the photodetectors and a better
match of the etalon with the grating monochromator as a first stage, it should
not be too technically challenging to reach the theoretical limit $g^{(2)}(0)=2$
for blackbody radiators.

The substantial increase in the observable photon bunching signature may
reopen considerations of intensity interferometry as a useful tool in
astronomical applications, where large baselines can be implemented with
significantly less effort than that in first order (amplitude) stellar
interferometers. Current atomic clock technology should 
even allow for the development of very large scale optical intensity
interferometers with kilometric baselines, leading to micro-arcsecond
resolution. With this, one could investigate stellar limb darkening coefficients, resolve spectroscopic binaries, measure stellar oblateness to better understand the internal dynamics, and also determine the pulsations of cepheid variables to better derive the Hubble constant \citep{barbieri:09, borra:13, dravins:07, dravins:13, foellmi:09, lebohec:08, naletto:09, nunez:12, ofir:06, rou:12, wittkowski:04}.

Beyond the traditional application of estimating the angular size of a
celestial object via its spatial coherence length, this technique may also be
used to measure photon decoherence of starlight across extended time-scales
and distances, which are impossible to create in Earth-based
laboratory conditions. Such experiments may provide complementary
observational evidence that help to set constraints in theoretical models of
quantum gravity \citep{lieu:03, maziashvili:09, ng:03, ragazzoni:03}. 

%% If you wish to include an acknowledgments section in your paper,
%% separate it off from the body of the text using the \acknowledgments
%% command.

%% Included in this acknowledgments section are examples of the
%% AASTeX hypertext markup commands. Use \url without the optional [HREF]
%% argument when you want to print the url directly in the text. Otherwise,
%% use either \url or \anchor, with the HREF as the first argument and the
%% text to be printed in the second.

\acknowledgments
We acknowledge the support of this work by the National Research Foundation and the Ministry of Education in Singapore.

%% To help institutions obtain information on the effectiveness of their
%% telescopes, the AAS Journals has created a group of keywords for telescope
%% facilities. A common set of keywords will make these types of searches
%% significantly easier and more accurate. In addition, they will also be
%% useful in linking papers together which utilize the same telescopes
%% within the framework of the National Virtual Observatory.
%% See the AASTeX Web site at http://aastex.aas.org/
%% for information on obtaining the facility keywords.

%% After the acknowledgments section, use the following syntax and the
%% \facility{} macro to list the keywords of facilities used in the research
%% for the paper.  Each keyword will be checked against the master list during
%% copy editing.  Individual instruments or configurations can be provided 
%% in parentheses, after the keyword, but they will not be verified.

%{\it Facilities:} \facility{Nickel}, \facility{HST (STIS)}, \facility{CXO (ASIS)}.

%% Appendix material should be preceded with a single \appendix command.
%% There should be a \section command for each appendix. Mark appendix
%% subsections with the same markup you use in the main body of the paper.

%% Each Appendix (indicated with \section) will be lettered A, B, C, etc.
%% The equation counter will reset when it encounters the \appendix
%% command and will number appendix equations (A1), (A2), etc.

%% The reference list follows the main body and any appendices.
%% Use LaTeX's thebibliography environment to mark up your reference list.
%% Note \begin{thebibliography} is followed by an empty set of
%% curly braces.  If you forget this, LaTeX will generate the error
%% "Perhaps a missing \item?".
%%
%% thebibliography produces citations in the text using \bibitem-\cite
%% cross-referencing. Each reference is preceded by a
%% \bibitem command that defines in curly braces the KEY that corresponds
%% to the KEY in the \cite commands (see the first section above).
%% Make sure that you provide a unique KEY for every \bibitem or else the
%% paper will not LaTeX. The square brackets should contain
%% the citation text that LaTeX will insert in
%% place of the \cite commands.

%% We have used macros to produce journal name abbreviations.
%% AASTeX provides a number of these for the more frequently-cited journals.
%% See the Author Guide for a list of them.

%% Note that the style of the \bibitem labels (in []) is slightly
%% different from previous examples.  The natbib system solves a host
%% of citation expression problems, but it is necessary to clearly
%% delimit the year from the author name used in the citation.
%% See the natbib documentation for more details and options.

\begin{thebibliography}{}

\bibitem[{Asakura(1970)}]{asakura:70}
Asakura, T. 1970, Opto-electronics, 2, 115

\bibitem[{Barbieri {et~al.}(2009)Barbieri, Daniel, de~Wit, Dravins, Jensen,
  Kervella, Bohec, Malbet, Nunex, Ralston, \& Ribak}]{barbieri:09}
Barbieri, C., Daniel, M.~K., de~Wit, W.~J., {et~al.} 2009, Astro2010: The
  Astronomy and Astrophysics Decadal Survey, 2010

\bibitem[{Boitier {et~al.}(2009)Boitier, Godard, Rosencher, \&
  Fabre}]{boitier:09}
Boitier, F., Godard, A., Rosencher, E., \& Fabre, C. 2009, Nature Physics, 5,
  267

\bibitem[{Borra(2013)}]{borra:13}
Borra, E.~F. 2013, Mon. Not. R. Astron. Soc.

\bibitem[{Davis {et~al.}(1999)Davis, Tango, Booth, ten Brummelaar, \& andS.
  M.~Owens}]{davis:99}
Davis, J., Tango, W.~J., Booth, A.~J., ten Brummelaar, T.~A., \& andS.
  M.~Owens, R. A.~M. 1999, Mon. Not. R. Astron. Soc., 303, 773

\bibitem[{Dravins(2007)}]{dravins:07a}
Dravins, D. 2007, High Time Resolution Astrophysics, Astrophysics and Space
  Science Library (Springer)

\bibitem[{Dravins \& LeBohec(2007)}]{dravins:07}
Dravins, D., \& LeBohec, S. 2007, Proc. SPIE, 6986, 698609

\bibitem[{Dravins {et~al.}(2013)Dravins, LeBohec, Jensen, Nu{\~n}ez, \& for~the
  CTA~Consortium}]{dravins:13}
Dravins, D., LeBohec, S., Jensen, H., Nu{\~n}ez, P.~D., \& for~the
  CTA~Consortium. 2013, Astroparticle Physics, 43, 331, elsevier

\bibitem[{Foellmi(2009)}]{foellmi:09}
Foellmi, C. 2009, Astronomy \& Astrophysics, 507, 1719, eSO

\bibitem[{Glauber(1963)}]{glauber:63}
Glauber, R. 1963, Phys. Rev., 131, 2766

\bibitem[{Green(1980)}]{green:80}
Green, J.~M. 1980, J. Phys. E: Sci. Instrum., 13, 1302, iOP

\bibitem[{Hanbury-Brown \& Twiss(1954)}]{hbt:54}
Hanbury-Brown, R., \& Twiss, R.~Q. 1954, Phil. Mag., 45, 663

\bibitem[{Hanbury-Brown \& Twiss(1956)}]{hbt:56}
---. 1956, Nature, 178, 1046

\bibitem[{Hanbury-Brown \& Twiss(1957)}]{hbt:57}
---. 1957, Proc. Roy. Soc., A242, 300

\bibitem[{Hanbury-Brown \& Twiss(1974)}]{hbt:74}
---. 1974, Mon. Not. R. Astron. Soc., 167, 121

\bibitem[{Karmakar {et~al.}(2012)Karmakar, Meyers, \& Shih}]{karmakar:12}
Karmakar, S., Meyers, R., \& Shih, Y. 2012, Proc. SPIE, 8518,
  doi:10.1117/12.929157

\bibitem[{Kurtsiefer {et~al.}(2001)Kurtsiefer, Zarda, Mayer, \&
  Weinfurter}]{kurtsiefer:01a}
Kurtsiefer, C., Zarda, P., Mayer, S., \& Weinfurter, H. 2001, J. Mod. Opt.,
  48(13), 2039

\bibitem[{LeBohec {et~al.}(2008)LeBohec, Barbieri, de~Wit, Dravins, Feautrier,
  Foellmi, Glindemann, Hall, Holder, Holmes, Kervella, Kieda, Coarer, Lipson,
  Malbet, Morel, Nuñez, Ofir, Ribak, Saha, Schoeller, Zhilyaev, \&
  Zinnecker}]{lebohec:08}
LeBohec, S., Barbieri, C., de~Wit, W.-J., {et~al.} 2008, Proc. SPIE, 7013:
  Optical and Infrared Interferometry, doi:10.1117/12.787443, optical and
  Infrared Interferometry

\bibitem[{Lieu \& Hillman(2003)}]{lieu:03}
Lieu, R., \& Hillman, L.~W. 2003, The Astrophysical Journal, 585, L77

\bibitem[{Mandel \& Wolf(1995)}]{mandel:95}
Mandel, L., \& Wolf, E. 1995, Optical Coherence and Quantum Optics (Cambridge
  University Press)

\bibitem[{Martienssen \& Spiller(1964)}]{martienssen:64}
Martienssen, W., \& Spiller, E. 1964, American Journal of Physics, 32,
  doi:10.1119/1.1970023

\bibitem[{Maziashvili(2009)}]{maziashvili:09}
Maziashvili, M. 2009, Astroparticle Physics, 31, 344

\bibitem[{Millour(2008)}]{millour:08}
Millour, F. 2008, New Astronomy Reviews, 52, 177

\bibitem[{Morgan \& Mandel(1966)}]{morgan:66}
Morgan, B.~L., \& Mandel, L. 1966, Physical Review Letters, 16, 1012

\bibitem[{Naletto {et~al.}(2009)Naletto, Barbieri, Occhipinti, Capraro, Paola,
  Facchinetti, Verroi, Zoccarato, Anzolin, Belluso, Billotta, Bolli, Bonanno,
  Deppo, Fornasier, Germanà, Giro, Marchi, Messina, Pernechele, Tamburini,
  Zaccariotto, \& Zampieri}]{naletto:09}
Naletto, G., Barbieri, C., Occhipinti, T., {et~al.} 2009, Astronomy \&
  Astrophysics, 508, 531, eSO

\bibitem[{Ng {et~al.}(2003)Ng, Christiansen, \& van Dam}]{ng:03}
Ng, Y.~J., Christiansen, W.~A., \& van Dam, H. 2003, The Astrophysical Journal,
  591, L87

\bibitem[{Nu{\~n}ez(2012)}]{nunez:12}
Nu{\~n}ez, P.~D. 2012, Master's thesis, University of Utah

\bibitem[{Ofir \& Ribak(2006)}]{ofir:06}
Ofir, A., \& Ribak, E.~N. 2006, Proc. SPIE, 6268: Advances in Stellar
  Interferometry, doi:10.1117/12.689123

\bibitem[{Park {et~al.}(2005)Park, Ko, \& Park}]{park:05}
Park, S., Ko, H., \& Park, M.-H. 2005, Optical Engineering, 44(4), 048001

\bibitem[{Ragazzoni {et~al.}(2003)Ragazzoni, Turatto, \&
  Gaessler}]{ragazzoni:03}
Ragazzoni, R., Turatto, M., \& Gaessler, W. 2003, The Astrophysical Journal,
  587, L1

\bibitem[{Rou {et~al.}(2012)Rou, Nunez, Kieda, \& LeBohec}]{rou:12}
Rou, J., Nunez, P.~D., Kieda, D., \& LeBohec, S. 2012, Mon. Not. R. Astron.
  Soc., 000, 1

\bibitem[{Saleh \& Teich(2007)}]{salehteich:07}
Saleh, B. E.~A., \& Teich, M.~C. 2007, Fundamentals of Photonics (John Wiley \&
  Sons)

\bibitem[{Sansonetti \& Veza(2010)}]{sansonetti:10}
Sansonetti, C.~J., \& Veza, D. 2010, J. Phys. B: At. Mol. Opt. Phys., 43,
  205003

\bibitem[{Scarl(1968)}]{scarl:68}
Scarl, D.~B. 1968, Phys.Rev, 175, 1661, issue 5

\bibitem[{Wittkowski {et~al.}(2004)Wittkowski, Aufdenberg, \&
  Kervella}]{wittkowski:04}
Wittkowski, M., Aufdenberg, J.~P., \& Kervella, P. 2004, Astronomy \&
  Astrophysics, 413, 711

\bibitem[{Wu {et~al.}(2014)Wu, Zhai, \& Liu}]{wu:14}
Wu, L.~A., Zhai, G.~J., Liu, X.~F., Chen, X.~H., Yao, X.~R., Yu, W.~K. Yu,
Optics Letters (in press)

\end{thebibliography}

%\bibliographystyle{apj} 
%\bibliography{PK}{}

\end{document}

%%
%% End of file `sample.tex'.
