Extrasolar Kuiper Belts

Abstract: Extrasolar debris disks are the dust disks found around nearby main sequence stars arising from the break-up of asteroids and comets orbiting the stars. Far-IR surveys (e.g., with Herschel) showed that ~20% of stars host detectable dust levels. While dust temperatures suggest a location at 10s of au comparable with our Kuiper belt, orders of magnitude more dust is required implying a planetesimal population more comparable with the primordial Kuiper belt. High resolution imaging (e.g., with ALMA) has mapped the nearest and brightest disks, providing evidence for structures shaped by an underlying planetary system. Some of these are analogous to structures in our own Kuiper belt (e.g., the hot and cold classical, resonant, scattered disk and cometary populations), while others have no Solar System counterpart. CO gas is seen in some debris disks, and inferred to originate in the destruction of planetesimals with a similar volatile-rich composition to Solar System comets. This chapter reviews our understanding of extrasolar Kuiper belts and of how our own Kuiper belt compares with those of neighbouring stars.

\section{Introduction}
\label{s:introduction}
As the other chapters in this book testify, studies of the Kuiper belt continue to play a
key role in shaping our understanding of the structure and history of the Solar System.
This chapter aims to provide some context to this relatively detailed information on our own
planetary system, by considering how this compares to the planetary systems of other stars,
in particular with regard to the component of their debris disks
which might be considered analogous to the Kuiper belt.
Such comparisons help to ascertain whether our system is typical or atypical, either
in terms of its current architecture, or in terms of its past history.
They also allow a deeper understanding of processes inferred to have occurred in our own history,
since those processes may be ongoing for some stars, furthermore providing evidence for how these
play out in different environments.

By now it is clear that our Solar System is just one of many planetary systems, with observations
of nearby stars showing that approximately half of stars have a planetary system in which at least
one of its planets may be detected in current surveys \citep{Winn2015}.
However, this seemingly abundant information is restricted to the planets that reside close to their
stars (within a few au), and information about the outer regions of planetary systems remains
scarce.
Indeed just a few percent of stars have planets $\gg 5$\,au detectable by direct imaging (which means at
least a few times more massive than Jupiter), although microlensing surveys hint that Neptune-analogues
may be more common.

In contrast, our understanding of the dust content of outer planetary systems is relatively well advanced
\citep[e.g.,][]{Wyatt2008, Hughes2018}.
It is over 30 years since the far-IR satellite {\it IRAS} first discovered circumstellar dust orbiting nearby main sequence
stars \citep{Aumann1984}.
Because this dust is short-lived it was recognised that it must be continually replenished from
larger planetesimals situated in a source region that was inferred from the dust temperature to be
at 10s of au \citep{Backman1993}, a conclusion that was reinforced by imaging of the dust structures \citep{Smith1984}.
Stars with such dust have been called Vega-type (after the first discovery), and the circumstellar dust
is referred to as a debris disk
\footnote{Here {\it debris} refers to the debris left over after planet formation, which
applies to any circumstellar material that is not a planet.
This could include a dust or gas component left over from the protoplanetary disk,
but also planetesimals and the dust and gas resulting from their destruction.},
and often interpreted to arise from an extrasolar analogue to the
Solar System's Kuiper belt \citep[e.g.,][]{Wyatt2003a, MoroMartin2008}.

The extent to which that analogy is appropriate
remains open for interpretation, but what is clear is that $\sim 20$\% of
stars have dust, and so presumably also planetesimals, at somewhere around 20-150\,au from their
stars \citep{Wyatt2008, Hughes2018}.
Given the aforementioned difficulty of detecting planets in these outer regions, it is usually uncertain
whether this dust lies at the outer edge of a planetary system or whether there are other similarities
with our own Kuiper belt, but there are clues in the dust structures.
These often reinforce the Solar System analogy, with the caveat that there is some anthropocentric
bias in this interpretation.

This chapter starts in \S \ref{s:observations} by summarising the various observational methods
used to obtain information about these putative extrasolar Kuiper belts, then in
\S \ref{s:perspective} considers how our Solar System fits within this context.
More detail is given in \S \ref{s:properties} about specific aspects about the physical and
dynamical properties of the extrasolar systems and their comparison with the Solar System,
before concluding in \S \ref{s:conclusions}.

\section{Extrasolar Kuiper Belt Observations}
\label{s:observations}
As noted in \S \ref{s:introduction}, most of the information about extrasolar Kuiper belts comes from observations
of dust created in the destruction of their planetesimals.
These observations can be roughly split into those used to discover the presence of dust using photometry (discussed
in \S \ref{ss:discovery}) and those that provide more detailed characterisation through high resolution
imaging (discussed in \S \ref{ss:characterisation}).
There is also an emerging area of observations of gas in the systems (see \S \ref{ss:gas}).

\subsection{Discovery: Photometry and SED Fitting}
\label{ss:discovery}
Debris disk dust is most readily detected from its thermal emission which manifests itself as the star appearing brighter
than expected from purely photospheric emission at wavelengths that depend on the dust temperature,
although care must be taken to ensure that potential extragalactic contamination has been removed from
the photometry \citep[e.g.,][]{Kennedy2012, Gaspar2014}.
The emission from dust created at 10s of au peaks at far-IR wavelengths and that is where the majority of debris
disks have been discovered.
The most recent surveys were those undertaken by {\it Herschel} \citep{Eiroa2013, Matthews2014, Hughes2018},
including an unbiased survey of the nearest several hundred stars in which dust emission was detected around
20\% of stars \citep{Thureau2014, Sibthorpe2018}.

These photometric surveys are supplemented by photometry at mid-IR and sub-mm wavelengths to build up
the spectral energy distribution (SED) of the dust emission, which generally shows that the emission is dominated
by a single temperature and so dust at a single distance (see e.g. Fig.~\ref{fig:1}).
Given the luminosity of the star this can be translated into a black body distance (i.e., the distance at which
the dust would be if it behaved like a black body), although it is recognised that the small dust that dominates
is an inefficient emitter and so this is likely an underestimate of the true location of the dust by a factor that
can be several depending on the star \citep{Booth2013, Pawellek2014}.
In addition to its black body radius $r_{\rm bb}$, the disk's SED is also characterised by its fractional
luminosity $f$ (i.e., the luminosity of the dust
divided by that of the star), with values of $10^{-6}$ to $10^{-3}$ typical for known disks.

\begin{figure}
\includegraphics[scale=0.51]{wyattfig1}
\caption{Observations of the nearby (7.8\,pc) 400\,Myr A-type main sequence star Fomalhaut.
{\it (Left)} Photometric observations show that the spectral energy distribution is comprised of
two temperature components, one from the star and another from cold dust at $\sim 70$\,K.
{\it (Right)} Optical images in which the stellar emission has been subtracted using coronagraphic techniques
to reveal scattered light from circumstellar dust that is distributed in a narrow ring at $\sim 130$\,au
from the star \citep{Kalas2008}.
The centre of this ring is offset from the star implying an eccentricity of $\sim 0.11$, and a planet-like object
(Fomalhaut-b) is seen close to the inner edge of the ring.}
\label{fig:1}
\end{figure}

For some systems the SED has been characterised in greater detail, in which case more realistic modelling
is performed to constrain the dust size distribution and composition \citep[e.g.,][]{Olofsson2012, Lebreton2012}.
However, the simple shape of the spectrum means that there is still value in considering the disks in terms
of their two observable parameters ($r_{\rm bb}$ and $f$), particularly because there are a number of
observational biases that can be readily understood within this context \citep{Wyatt2008}.
Nevertheless, some systems have an SED that is best characterised as having two temperatures \citep{Chen2006, Kennedy2014}.
In the context of an extrasolar Kuiper belt interpretation, this means that there is an additional
warm component of emission that might be analogous to the Solar System's zodiacal cloud.
However, there remains debate about whether this component originates in an extrasolar asteroid belt
analogue \citep{Su2013}, exocomets scattered in from the extrasolar Kuiper belt \citep{Wyatt2017}, or is a transient
dust component perhaps caused by a recent collision \citep{Jackson2012}.

\subsection{Characterisation: Imaging}
\label{ss:characterisation}
Debris disks around the nearest stars have a spatial extent of $>1$\,arcsec which means that these can be readily
resolved.
The low resolution of far-IR instrumentation precludes this except in a few cases \citep[e.g.,][]{Acke2012},
but this technique has been
successful {at optical and near-IR wavelengths} with {\it HST}
(see Fig.~\ref{fig:1}), and more recently SPHERE on VLT and GPI on Gemini,
where starlight scattered by the smallest sub-$\mu$m
dust can be imaged \citep[e.g.,][]{Schneider2014}.
In the sub-mm with ALMA (Atacama Large Millimeter/submillimeter Array),
thermal emission from mm to cm-sized grains can be
imaged as well (see Fig.~\ref{fig:5}).
Given the size of dust dominating the emission at the different wavelengths, it is generally assumed that
sub-mm observations trace the distribution of the parent planetesimals, while the picture is more complicated
at optical wavelengths because the orbits of small dust grains are significantly modified by radiation forces.
Indeed, different radial and even nonaxisymmetric structures are seen in multiwavelength imaging of the
same disk providing a valuable diagnostic tool for the underlying dynamics.

Such imaging provided the first evidence that the dust is configured in a disk \citep[rather than a spherical distribution,][]{Smith1984},
and has allowed direct measurement of the radial location of the dust, its radial and vertical extents and the presence of gaps
(see \S \ref{ss:radial}),
as well as providing evidence for asymmetries in the form of clumps, eccentricities and warps (see \S \ref{ss:structures}).
All of these give vital clues to the underlying dynamics and the existence of planets.

\subsection{Gas}
\label{ss:gas}
Debris disks used to be considered as gas-free, in contrast to the protoplanetary disks found around
young stars ($<10$\,Myr).
Recent discoveries with ALMA show that this can no longer be considered the case, since there are now of order
10 main sequence stars known to exhibit CO gas emission.
Many of the first gas detections are around the more massive A-type stars \citep{Kospal2013, Greaves2016, Moor2017},
though there are a number of observational biases that facilitate the detection of gas around early-type stars,
and gas is now seen around later spectral types as well \citep{Marino2016, Matra2019a}.
Most of these observations can be explained within the context of a model in which the gas in debris disks,
like the dust, is the product of of the destruction of planetesimals that are then required to be volatile-rich
\citep{Kral2017}.
An open question is whether the youngest ($<50$\,Myr) stars retain gas from the protoplanetary disk phase
\citep{Kospal2013, Kral2018}.
Observations of debris disk gas not only allow the mass and composition of gas to be measured (see \S \ref{ss:composition}),
but also spatially resolved
with additional kinematic information, resulting in a potentially powerful probe of the disk structure and dynamics
\citep[see \S \ref{ss:structures},][]{Dent2014}.

\section{An Extrasolar Perspective of the Kuiper Belt}
\label{s:perspective}
It is important to note that the Solar System's debris disk is fainter than any known debris disk around another star.
If one of the stars in the debris disk surveys hosted an exact replica of the Solar System, it would only have
been possible to detect the photospheric emission from the star and not that of circumstellar dust (see Figs.~\ref{fig:2}
and \ref{fig:3}).
The challenge with detecting a true Kuiper belt analogue lies in the fact that its thermal emission is outshined by more
than an order of magnitude at all wavelengths by the Sun.

\begin{figure}
\includegraphics[scale=0.49]{wyattfig2}
\caption{The Solar System's Kuiper belt.
{\it (Left)} Face-on images of the number density distribution of Kuiper belt objects that have been debiased from the
observed distributions.
To the left is shown the full Kuiper belt which appears as a narrow ring dominated by the classical Kuiper belt
objects \citep{Matthews2016}, and to the right the contribution of resonant Kuiper belt objects that shows two clumps
at $\pm 90^\circ$ from Neptune (shown with a circle) due to the Plutino population \citep{Lawler2014}.
{\it (Right)} The evolution of the 70\,$\mu$m emission from dust in the Kuiper belt (shown as its level compared with
that from the Sun) expected in the model of \citet{Gomes2005}, and compared with the dust emission from nearby
Sun-like stars \citep{Booth2009}.
In this model the massive primordial Kuiper belt is comparable in brightness to those observed around nearby stars
(shown with asterisks),
but is depleted to a non-detectable level at $\sim 800$\,Myr by the Late Heavy Bombardment.
}
\label{fig:2}
\end{figure}

In fact the level of thermal emission from dust in the Kuiper belt is not well known.
Its emission has not been detected since it is masked by that of the zodiacal cloud which is much closer to the
Earth, although the far-IR all-sky surveys like {\it IRAS}, {\it COBE} and {\it Planck}
do provide upper limits on the level of emission
that can be present \citep[e.g.,][]{Backman1995}.
Direct detection of dust grains in the outer Solar System gives an indication of the number of dust grains present
\citep{Piquette2019}, but this must
be combined with dynamical models for dust production and evolution to understand the whole dust structure \citep{Vitense2014}.
In doing so, the best estimate is that Kuiper belt emission peaks at $\sim 70$\,$\mu$m at a flux 1\% that of the Sun \citep{Vitense2012}.
Even for bright stars their far-IR fluxes cannot be predicted with such accuracy, which combined with calibration uncertainties
in photometric measurements, means that $\gg 10$\% excesses are required for a confident detection.

If it were possible to resolve the emission from Kuiper belt dust, the structures that would be observed
at long wavelengths would resemble the structures known to be present in the distribution of Kuiper belt objects
\citep[e.g.,][]{Lawler2014};
i.e., there would be a prominent ring at $\sim 40$\,au from the classical Kuiper belt
(see Fig.~\ref{fig:2} left),
with some non-axisymmetric structure caused by the resonant Kuiper Belt Objects (see Fig.~\ref{fig:2} middle),
and with emission extending out to larger radii from the scattered disk (which is present but hardly noticeable in
Fig.~\ref{fig:2} left).
Given the low density of the dust distribution, the smallest dust (as traced by the shorter far-IR wavelengths) would
migrate in by Poynting-Robertson drag and so fill in the 40\,au hole.
However, very little dust would make it past Jupiter or even Saturn, meaning that this tenuous dust distribution
would have drops in density associated with these planets \citep{Liou1999, MoroMartin2002}.

There are, however, indications that the Kuiper belt was more massive in the past by several orders of
magnitude.
At such an epoch there would have been correspondingly more dust, resulting in detectable levels of
emission \citep[see Fig.~\ref{fig:2} right,][]{Booth2009}.
Thus the known debris disks could be analogues to the primordial Kuiper belt, and thus representative of
systems that either never will, or have yet to undergo, large-scale instabilities.
Note though that the exact timing of the depletion in the Solar System remains uncertain \citep[e.g.,][]{Gomes2005, Morbidelli2018}.

\section{Extrasolar Kuiper Belt Properties (and Comparison with Solar System)}
\label{s:properties}

\subsection{Mass and Radius Distribution}
\label{ss:mass}
Far-IR surveys from the past few decades have provided large quantities of information
on the incidence of debris disks, as measured by the excess emission above photospheric
levels at different wavelengths (typically at 24, 70, 100 and 160\,$\mu$m).
Surveys have been carried out for nearby stars with a variety of spectral types, distances and
ages.
While main sequence star ages are usually quite inaccurate, there is sufficient information, particularly
when comparing associations of young stars for which ages can be better determined, to show that
disks tend to get fainter with age \citep{Rieke2005, Su2006}.
As such the surveys are usually interpreted within the context of population models that take
collisional erosion into account (see \S \ref{ss:evolution}).

\begin{figure}
\hspace{0.5in}
\includegraphics[scale=0.75]{wyattfig3}
\caption{Distribution of debris disk luminosities and radii \citep[figure adapted from][]{Sibthorpe2018}.
The colour scale shows the fraction of a sample of $\sim 300$ Sun-like (FGK) stars from \citet{Phillips2010}
that have dust emission as a function of their measured fractional luminosities and black body radii.
The approximate locations of asteroid belt and Kuiper belt dust are shown in yellow, noting that the black body
radius is a factor of a few smaller than the true radius because dust emits hotter than black body.
The contours show the fraction of stars in the sample for which dust emission could have been detected,
which is used to correct the observed incidence of debris to get the fraction plotted in the colour scale.
The blue dashed lines and annotation explain why the known disks (which are present around $\sim 20$\% of stars)
lie mostly in a well defined region of parameter space.
}
\label{fig:3}
\end{figure}

Fig.~\ref{fig:3} summarises the state of our knowledge of the debris disk population seen around nearby
Sun-like (i.e., F, G or K-type) stars.
For each of the $\sim 300$ stars in the sample, the various observations (be they detections or non-detections
of dust emission) have been combined to determine the best fit fractional luminosity and black body radius to their disk
(see Fig.~\ref{fig:1}),
or the upper limit on the dust luminosity for a given temperature to allow a non-detection.
Particularly where disks are only detected at a single wavelength the fitted parameters can have significant
uncertainty, and so the disks are placed on this plot accordingly
\footnote{More specifically, uncertainties in the photometric data mean that analysis of the
spectral energy distribution can only determine a probability that the disk is present in a given
part of fractional luminosity vs black body radius parameter space.}, and the colour scale shows the fraction of stars
in the sample predicted to have a disk in different regions of parameter space \citep[for more details on
the method, see][]{Sibthorpe2018}.
Since observational biases mean that disks of a given level could only have been detected around a fraction
of the stars in the sample as shown in the contours \citep{Wyatt2008}, this
is taken into account to get an estimate of the underlying distribution.

As noted in the figure, the distribution of known disks is bounded on the bottom edge by an
observational limitation, while the top edges are set by physical considerations.
For disks close to the star ($<40$\,au in black body radius) the maximum fractional luminosity is set by collisional erosion,
since most nearby stars have the same roughly Gyr age and so their disks would be expected to have undergone similar
levels of collisional processing and would have undergone more depletion the closer they are to the star
\citep{Wyatt2007b}.
Some disks may exhibit higher fractional luminosities if there happens to be a young star in the sample, or if there is
a transient dust component \citep{Wyatt2007a}.
For disks further from the star, the upper edge is set by a constraint on the mass required to
produce the given fractional luminosity.

While Fig.~\ref{fig:3} gives a quick visual impression of the distribution of debris disk properties, population models
that also address the age dependence of the phenomenon can be used to get a more accurate parameterisation
for the distributions of initial masses and radii of the disks.
Such models show that it is possible to fit the observations by assuming that all stars started out with a planetesimal
belt, with a radius and initial mass chosen from given distributions, that mass then evolving with time due to collisional erosion
\citep{Wyatt2007b, Lohne2008, Gaspar2013, Sibthorpe2018}, although it is possible that a few
of the brightest disks may require unrealistic initial masses \citep[$\sim 1000M_\oplus$,][]{Krivov2018}.
There is evidence that the disks of A-type stars tend to be brighter than their Sun-like star counterparts \citep{Greaves2003},
even accounting for the various biases (such as the on average younger ages of the more massive stars), which may be because of a difference
in their planetesimal belt properties (e.g., in having more massive disks, or one made of smaller planetesimals, or lower levels of
stirring).
Given that many of the well known bright debris disks discussed in this review are found around A stars (like $\beta$ Pic),
it may be worth bearing this caveat in mind when making comparisons with the Solar System.
There are few disks detected around the less massive M stars, but this is likely an observational bias and there is no
evidence as yet to say that their disks are any different to those of Sun-like stars \citep[e.g.,][]{Morey2014}.

Despite the success of these models, it is worth noting that they cannot make any strong conclusions about the 80\% of stars
that do not have detectable disks, since these would be based on assumptions about the shape of the distributions.
For example, a common assumption is that the mass distribution is log-normal, with the non-detections arising from
close-in disks that undergo rapid depletion and so are not detectable in the nearby star population, but it could be that the
mass distribution is bimodal, with 80\% of stars having no disk.
Also, more recently a correlation has been found between planetesimal belt radii determined from
high resolution imaging of the sub-mm dust emission and the luminosity of the star \citep{Matra2018b}.
If this holds up as more disks are resolved, this would suggest a preferential location for the formation of
planetesimal belts that could be linked to ice-lines in the protoplanetary disk for example, which would provide an
additional constraint on the population models (as well as significant insight into the origin of the belts).

Fig.~\ref{fig:3} also includes the location of the Solar System's Kuiper belt and Asteroid belt
(see \S \ref{s:perspective}).
This shows that roughly 20\% of stars have disks that are comparable in radius to the Kuiper belt,
in having dust at temperatures that give a black body radius of 10-100\,au, but that have at least
an order of magnitude more dust at fractional luminosities of $10^{-6}$ to $10^{-3}$.
For the reasons noted above it is not possible to quantify where our own Kuiper belt fits into the
distribution.

\subsection{Evolution: Collisional vs Dynamical Erosion}
\label{ss:evolution}
As noted in \S \ref{ss:mass}, the evidence that there is any evolution in debris disk properties comes
not from seeing those changes in any one system, but rather from seeing the change in properties for samples
of stars with different ages.
Inevitably then there remains uncertainty in how individual disks evolve.
However, since it is necessary for collisions to occur to replenish the dust that is observed, and that
dust is lost from the system relatively quickly (compared with the age of the star), some level of collisional
erosion is unavoidable.
That evolution has been studied extensively both analytically and numerically, in ways that have also been applied
to small body populations in the Solar System.
In the Solar System evidence for past collisional evolution can be found in the detailed shape of the asteroid and
Kuiper belt size distributions \citep[e.g.,][]{Bottke2005, Fraser2009, Singer2019}, whereas in extrasolar systems the evolution of the
size distribution in debris disks
is manifested in the evolution of disk brightness which must be inferred from the population studies.

There are several unknowns in the problem, such as the initial size distribution of the planetesimals,
their level of stirring (i.e., their collision velocity), and their size-dependent dispersal threshold.
With these factors known, the collisional evolution would be well understood \citep[e.g.,][]{Wyatt2011},
and when coupled with models for dust optical properties that determine both its emission properties and loss due to
radiation forces, this can be readily translated into an evolution in the disk's observables.
In principle this means that these fundamental physical parameters can be extracted from the observations,
although in practise there are degeneracies between the parameters which make it hard to extract unambiguous
information, other than that the vast majority of observations are consistent with purely steady state
collisional erosion.

Indeed it is easier to identify the disk properties or processes that are incompatible with the observations,
at least on a population level.
For example, a slow decline in disk incidence with age that is consistent with collisional erosion means
that large scale dynamical depletions at 100s of Myr must be relatively rare \citep{Booth2009}, otherwise there would be
fewer bright disks at several Gyr ages.
Thus either the Kuiper belt is anomalous in this respect \citep[i.e., if its depletion is responsible for the
Late Heavy Bombardment,][]{Gomes2005}, or the depletion actually occurred much earlier on \citep{Morbidelli2018}.
It is also possible to conclude that the majority of debris disks are not broad planetesimal disks \citep[even if each is only
bright at one radius at any given time, e.g.,][]{Kenyon2010},
since broad disks are hot early on and cold later on in a manner that is incompatible
with the observed 24 and 70\,$\mu$m brightness evolution \citep[][see also \S \ref{ss:radial}]{Kennedy2010}.

However, for comparison with the Kuiper belt it is interesting to consider what has been concluded about the
largest planetesimals in debris disks and the origin and level of stirring within them.
While observations are only sensitive to dust up to $\sim $cm-size, the short collisional lifetime of such dust means that
it must be continually replenished.
Extrapolating the size distribution and equating the resulting collisional lifetime with age of the star shows that
debris disks are replenished by planetesimals at least several km in size \citep{Wyatt2002}, with population models
pointing to sizes up to 100\,km \citep{Gaspar2013}.
It has been suggested that some disks may be comprised of much smaller planetesimals, say 10-100\,m, that are colliding
at very low velocity so that the collisions are barely disruptive \citep{Heng2010, Krivov2013}, explaining the lack of small
grains inferred to some systems \citep{Eiroa2011}, although chance alignment of background galaxies could be confusing the
interpretation of these stars \citep{Gaspar2014}.
Detailed observations of the halo of small grains seen to extend beyond a debris disk due to radiation pressure
has the potential to constrain the level of stirring, with eccentricities inferred to be $<1$\%
in one system \citep{Thebault2008}.
For some disks the vertical structure has been measured with inclinations inferred to be of order a few percent
\citep[][]{Matra2019b, Daley2019}.

The origin of the stirring is not well constrained.
It could arise from large bodies embedded in the disk that excite relative velocities to around their escape
velocity \citep{Kenyon2010}.
However, the growth of Pluto-sized objects at 100s of au in the disk takes too long to explain large disks
seen around Gyr-old stars.
One possibility is that these form directly from the protoplanetary disk (rather than from the collisional growth
of smaller planetesimals), or it could be that a massive disk of $\sim 200$\,km planetesimals can stir itself
without the need for larger bodies \citep{Krivov2018b}.
If the belts do mark the outer edge of the planetary system then those planets could stir the disk through their
gravitational perturbations.
This could be through overlapping resonances from a planet at the inner edge \citep{Quillen2006}, but
could also arise from the secular perturbations of a more distant eccentric or inclined planet \citep{Mustill2009}.
It is also possible that the planetesimals were born on stirred orbits, e.g. if these formed in the early
phases of the protoplanetary disk \citep{Booth2016}.

To illustrate the range of interpretations for debris disks that do not encompass a traditional Kuiper belt-like
analogy, note that in one model the planetesimals form during protoplanetary disk dispersal from the
dust that is swept up by a photoevaporating gas disk \citep{Carrera2017}.
While such a model does not recover the observed distribution of debris disk radii, it is worth bearing in mind that
our understanding of protoplanetary disk dispersal processes is incomplete,
in particular regarding the fate of the $>1M_\oplus$ of mm-sized dust that
resides in the outer regions of the system \citep{Wyatt2015}.
As such it is not inconceivable that our understanding of the origin and evolution of debris disks could need
revision.

\subsection{Planetesimal Composition}
\label{ss:composition}
A standard way of determining the composition of circumstellar material is to look for solid state
features in the emission spectrum, such as the silicate features near 10\,$\mu$m.
This has been successful for the few systems for which warm dust is present \citep[e.g.,][]{Lisse2009},
however for the majority that are too cold to have significant mid-IR emission, this leaves far-IR features
of which only a few were measured with {\it Herschel}.
Nevertheless, when applied to the disk of $\beta$~Pic, the 69\,$\mu$m feature was used to show that 3.6\%
of the dust mass is in crystalline olivine grains, moreover constraining them to be relatively Mg-rich
\citep{deVries2012}.

As noted in \S \ref{ss:gas}, a new technique has been developed to probe the volatile content of the planetesimals
by studying gas in the systems.
Molecular gas has a short and well defined lifetime in optically thin environments since it is destroyed by interstellar
radiation, e.g., CO has a lifetime of just $\sim 120$\,years.
Thus if a disk is observed to have a certain mass of CO, the rate at which CO is being replenished can be readily inferred.
Since observations of the dust disk can also be used to determine the rate at which the refractory component is being
replenished, these two rates can be used to find the volatile fraction of the planetesimals feeding the cascade.
While there are some assumptions in that calculation, such as that both gas and dust are being replenished from
steady state grinding, and the uncertainties in the CO mass measurement can be large if excitation conditions are not
known \citep[e.g.,][]{Matra2015}, this at least provides an order of magnitude estimate.

\begin{figure}
\hspace{1.4cm} \includegraphics[scale=0.50]{wyattfig4}
\caption{Planetesimal volatile fraction inferred from observations of gas in debris disks compared with
the same fractions for Solar System comets \citep{Matra2017}.
Volatile fraction here is defined as the fraction of the total mass that is made up of CO and CO$_2$.}
\label{fig:4}
\end{figure}

For now this has been achieved for a handful of disks with conclusions being that the fraction of the planetesimals
that is made up of CO and CO$_2$ is similar to that of Solar System comets
\citep[see Fig.~\ref{fig:4},][]{Marino2016, Matra2017}.
Moreover the same model also explains those systems for which CO is not detected \citep{Kral2017}.
Thus it is plausible that all debris disks are made of volatile-rich planetesimals, and so that all have associated
secondary gas disks, but that only a fraction of these have CO at detectable levels.
While no other gas molecules have been detected in debris disks, this is to be expected given the lower abundance
relative to CO and moreover the even shorter lifetime of most molecules \citep{Matra2018a}.

Atomic gas in the form of CI, CII and OI has also been detected \citep[e.g.,][]{RiviereMarichalar2014, Higuchi2017},
where it is inferred that this gas is the photodissociation product of CO and CO$_2$.
These are much longer lived and survive until the atomic gas disk has undergone viscous diffusion to be accreted
onto the star \citep{Kral2016}.
The observed ratio of carbon to CO thus constrains the viscous diffusion timescale, with implications for the
underlying physics of that accretion \citep{Kral2016b}.
For $\beta$ Pic, this resulted in a prediction for the amount of OI that was too low compared with that observed,
which was used to infer that the planetesimals were also water-rich, similar to Solar System comets, since that
would explain the detection \citep{Kral2016}.

While there is an emerging paradigm that has been outlined above, it is worth noting that there remain uncertainties.
For example, the model does not explain the distribution of CI observed in $\beta$ Pic \citep{Cataldi2018}.
It was also suggested that the $<50$\,Myr systems with protoplanetary disk levels of CO have a component of primordial
gas that has yet to disperse \citep{Kospal2013, Moor2017}.
However, it was recently shown that such high levels of CO could be explained in a secondary production scenario,
since if the gas production rate is particularly high, or the viscous diffusion timescale is too long, then
sufficient carbon can accumulate to shield the CO from photodissociating, which then can also start to self-shield
\citep{Kral2018}.

In consideration of the Solar System, the points to take away are that the volatile-rich composition of our Kuiper belt
may be similar for other stars, and also that there may have been a tenuous atomic gas disk present early on.

\subsection{Radial and Vertical Structure}
\label{ss:radial}

\begin{figure}
\includegraphics[scale=1.04]{wyattfig5}
\caption{Compilation of ALMA images of the distribution of mm-cm-sized dust in the debris disks of 15 nearby stars.
Individual images are from: $\beta$ Pic \citep{Matra2019b}, HD107146 \citep{Marino2018b}, q$^1$ Eri (J. Lovell, private communication),
HD181327 \citep{Marino2016}, HD92945 \citep{Marino2019}, $\eta$ Crv \citep{Marino2017a}, HR4796 \citep{Kennedy2018},
HD98800 \citep{Kennedy2019}, Fomalhaut \citep{MacGregor2017}, $\epsilon$ Eri \citep{Booth2017},
AU Mic \citep{Daley2019}, HR8799 \citep{BoothM2016}, HD32297 and HD61005 \citep{Schneider2014, MacGregor2018},
49 Cet \citep{Hughes2017}.
}
\label{fig:5}
\end{figure}

As Fig.~\ref{fig:1} shows, debris disks can be radially very narrow, with a fractional radial width (i.e., $\Delta r/r$)
of $\sim 10$\% in the example shown of the Fomalhaut disk.
That width can be measured for disks that have been imaged, but more specifically requires that disks
be resolved at mm wavelengths to
determine the width of the planetesimal belt, since the smaller dust probed at shorter wavelengths will be spread out
from the parent belt by radiation forces.
There is then some bias in the sense that a broad disk will have lower surface brightness than a narrow one, and so will
be harder to detect in high resolution imaging.
Nevertheless, Fig.~\ref{fig:5} shows that it is generally found that the disks are perhaps better described as rings or belts,
with the fractional radial width ranging from 10-100\%.

There are some notable exceptions to this description in that some disks are broad \citep[e.g.,][]{MacGregor2016}.
For example, the 61~Vir disk extends from 30\,au to beyond 150\,au
\citep{Wyatt2012, Marino2017b}, a property which means that its surface brightness
is so low that it is barely noticeable in ALMA imaging (but detected in the visibility data), with a radial profile that may be explained by
a broad disk of $\sim 10$\,km-sized planetesimals undergoing collisional erosion \citep{Marino2017b}.
Other examples of broad disks include HD107146 and HD92945, which also show an annular gap in the radial profile
\citep[see Fig.~\ref{fig:5},][]{Marino2018b, Marino2019} that may be related to planets either embedded in the disk or nearby
\citep{Yelverton2018}.

The inner edges of the debris belts have received some attention, since it is noted that if the belt is truncated by
dynamical interactions with a planet located at the inner edge, then the sharpness of the belt edge is indicative
of the mass of the planet \citep[e.g.,][]{Chiang2009}.
A more massive planet stirs up the disk to a greater extent resulting in a shallower slope.
Analogy with the Solar System shows that the sharpness and location of the disk's outer edge can also
be an important indicator of the system's past history \citep{Gladman2001}.

The outer edge of a debris disk is sometimes quoted as the location at which the disk becomes too faint
to detect.
However, it is expected that for planetesimals confined to a belt, there would be
a halo of small $\mu$m-sized dust that extends beyond the belt due to radiation pressure \citep{Strubbe2006}.
In this case the outer edge of the planetesimal belt would be evident
as the location where there is a drop in the disk's surface brightness.
However, it is now becoming clear that
the extended dust components in some systems that were thought to be a halo of small dust
also contain mm-sized dust \citep{MacGregor2018}.
Since this dust is unlikely to be affected by radiation pressure, this calls into question whether
the planetesimal disk is in fact broad with simply a step change in density at what had been assumed to be the outer
edge of a radially confined planetesimal belt.
The presence of a faint outer belt of debris in one system could support this view \citep{Marino2016}.
Alternatively the extended halo of mm-sized grains could represent an analogue of the scattered disk,
i.e., made of planetesimals (and the dust derived from their destruction) that are on eccentric orbits
undergoing scattering with a planet located near the inner edge.
Indeed the favoured interpretation of the radial structure of the HR8799 disk is that it has a low eccentricity disk
(i.e., a classical Kuiper belt analogue) as well as a high eccentricity population of comets in a scattered disk analogue
\citep[see Fig.~\ref{fig:6},][]{Geiler2019}, an interpretation strengthened by the fact that multiple planets interior to the belt are
known in this system \citep{Marois2010}.

\begin{figure}
\includegraphics[scale=0.46]{wyattfig6}
\caption{Distributions of orbital elements inferred from ALMA images of debris disks that are analogous to dynamical
populations in the Kuiper belt.
{\it (Left)} Vertical structure in the $\beta$ Pic disk (see Fig.~\ref{fig:5}) is best modelled as the the sum of two
Gaussians that imply a distribution of inclinations reminiscent of the hot and cold classical Kuiper belt \citep{Matra2019b}.
{\it (Right)} Radial structure in the HR8799 disk (see Fig.~\ref{fig:5}) includes a halo of mm-sized grains that is
best modelled with a high eccentricity population potentially analogous to the scattered disk \citep{Geiler2019}.
}
\label{fig:6}
\end{figure}

It was already noted in \S \ref{ss:evolution} that the vertical structure of some disks has now been resolved.
Furthermore, for $\beta$ Pic it was found that the vertical structure was not well described by a single Gaussian, rather it
is fit well as the sum of two Gaussians that correspond to components with mean inclinations of $1^\circ$ and $9^\circ$
\citep{Matra2019b}.
While this is a few times flatter than the Solar System, there is a clear analogy with the cold and hot components
of the classical Kuiper belt (see Fig.~\ref{fig:6}).

\subsection{Dynamical Structures}
\label{ss:structures}
Some of the radial structures discussed in \S \ref{ss:radial} can be considered dynamical structures, in that
the shape of inner and outer disk edges, and gaps within the disks, can all be related to the presence of
planets.
However, such an interpretation is potentially ambiguous, because this may alternatively be indicative of the initial
distribution of planetesimals, in which radial structure may have been imprinted by processes ongoing
during the protoplanetary disk phase (e.g., the formation of planetesimals at ice-lines);
i.e., no planets may be required.
Non-axisymmetric structure is a clearer indicator that additional perturbations are ongoing, since these
would be needed to maintain such structures over long timescales.

Several of the disks that have been imaged are seen to be eccentric.
This was manifested originally by a brightness asymmetry in mid-IR images of the HR4796 disk \citep{Telesco2000},
which was interpreted as arising from the pericentre side of the disk being closer to the star and so
hotter and brighter \citep{Wyatt1999}.
That interpretation is supported by ALMA imaging which showed that the opposite side of the disk is brighter
\citep{Kennedy2018} as expected due to the slower orbital velocity at apocentre \citep{Pan2016}.
Subsequently scattered light imaging showed that the star is not exactly at the centre of the disk
\citep{Milli2017}, although the direction of the offset is not exactly that expected suggesting that
there may be some additional density variations around the ring \citep[see also][]{Lohne2017}.
Fomalhaut's disk (see Fig.~\ref{fig:1}) is perhaps a better example, being brighter and larger on the sky so that
the offset and asymmetries are more clearly seen \citep{Kalas2005, MacGregor2017}, with additional azimuthal structure
also evident \citep{Kalas2013}.

The origin of these eccentricities, which are of order 0.05-0.15, have not been identified, but would arise naturally
from the secular perturbations of a planet on an eccentric orbit.
While a planet candidate was detected orbiting near the inner edge of
the Fomalhaut disk \citep[see Fig.~\ref{fig:1},][]{Kalas2008} its eccentricity
was subsequently found to put it on a disk-crossing orbit \citep{Kalas2013} that would significantly disrupt the disk
suggesting that the planet is less massive than originally thought and/or has only been put on this orbit
recently \citep{Beust2014, Tamayo2014, Lawler2015}.
The presence of an additional planet in the system that is responsible for the eccentricity is thus inferred.

A more secure identification of a disk asymmetry that is caused by planetary perturbations is the warp in the
$\beta$ Pic disk \citep{Mouillet1997}.
The change in the plane of symmetry seen in scattered light images at 80\,au in this edge-on disk \citep[see][]{Apai2015}
was inferred to arise from the secular perturbations of a planet on an inclined orbit at $\sim 10$\,au.
Since that planet has subsequently been found with direct imaging \citep{Lagrange2010}, this supports the idea that
disk structures can be used to identify unseen planets.
That said, there remains some uncertainty in the interpretation of the vertical structure of the $\beta$ Pic disk,
not least because the additional clumpy structure mentioned below is situated at the same location as the warp.

Both the Vega and $\epsilon$~Eri debris disks were inferred to be clumpy based on sub-mm imaging
\citep{Holland1998, Greaves1998}.
Subsequent imaging gave conflicting views on the presence of the clumps \citep{Greaves2005, Marsh2006, Hughes2012, Holland2017},
but it seems likely that at least some of the clumpiness was attributable to background galaxies.
The most secure evidence for clumpiness in a disk comes from $\beta$ Pic for which mid-IR images show a prominent
clump at $\sim 50$\,au in projected separation from the star \citep{Telesco2005}.
The same clump was later imaged in CO \citep[see Fig.~\ref{fig:7},][]{Dent2014}, for which the additional kinematic information showed that
the clump is both radially and azimuthally broad and situated at $\sim 85$\,au distance from the star \citep{Matra2017b}.
Two models were proposed to explain the clump, either a recent collision between two Mars-sized protoplanets
\citep{Jackson2014}, or the trapping of planetesimals into the 2:1 and 3:2 resonances of a $\sim 30M_\oplus$ planet
that migrated out from 31-57\,au over 23\,Myr \citep{Wyatt2003, Matra2019b}.
The radial breadth of the clump disfavours a collisional origin leaving the favoured interpretation as a structure
analogous to the resonant Kuiper belt objects (the Plutinos and Twotinos).
Comparison of Figs.~\ref{fig:2} and \ref{fig:7} show there are some notable differences,
however, such as that the tight clump observed in the $\beta$ Pic disk would require the planetesimals to have
a smaller libration width than resonant objects in the Kuiper belt.
This could be explained within the context of this model
by the stochasticity of the migration and/or the initial orbital eccentricities \citep[e.g.,][]{Reche2008, Nesvorny2016}.

\begin{figure}
\hspace{2cm} \includegraphics[scale=0.60]{wyattfig7}
\caption{Dynamical structure in the debris disk of $\beta$ Pic that is analogous to that of resonant Kuiper belt
objects.
While the disk is edge-on, velocity information from CO observations allows its face-on distribution to be
inferred showing that this includes two clumps of unequal brightness on opposite sides of the star \citep{Dent2014}.
The best explanation is that outward migration of a planet close to the inner edge of the disk trapped volatile-rich
planetesimals into its 2:1 and 3:2 resonances causing over-densities in these regions where collisions are
most frequent \citep{Wyatt2003}.
The CO lifetime is short compared with the orbital timescale meaning that it does not travel far from where
it is created before being photodissociated.
}
\label{fig:7}
\end{figure}

\subsection{Connection from Outer to Inner System}
\label{ss:connection}
The fraction of stars with excess 10\,$\mu$m emission indicative of warm dust in the region where there
may be habitable planets was long recognised to be low \citep{Aumann1991}, and was quantified more recently
from the WISE photometric survey to find that excesses above 10\% of the stellar emission only
occur to 1:$10^4$ nearby stars \citep{Kennedy2013},
although this rate is higher if considering a sample of young $<100$\,Myr stars.
This result is explained by collisional evolution \citep[see also Fig.~\ref{fig:3} and \S \ref{ss:mass},][]{Wyatt2007a},
which suggests that the rare bright excesses are transient populations of warm dust, perhaps originating in recent
collisions.
However, ongoing interferometric surveys have shown that excesses at 0.1-1\% levels are more common, and present
to $\sim 20$\% of stars \citep{Ertel2018}, and moreover are more prevalent (though not exclusively) around
stars known to host cold outer disks \citep{Mennesson2014}.

This does not prove that the warm dust must originate somehow in the outer belt, since it could be that the
conditions of a protoplanetary disk that leads to the formation of an outer cold belt also favour the formation
of a massive asteroid belt \citep[e.g.,][]{Geiler2017}.
However, there are two mechanisms by which outer belts can replenish a warm dust population.
The first is by Poynting-Robertson drag \citep{Burns1979}.
While the dust population dragged in from the known outer belts is expected to be significantly depleted by mutual
collisions \citep[which is not the case in the more tenuous belt in the Solar System,][]{Wyatt2005},
the levels of dust expected to make it in to $\sim 1$\,au are comparable to those inferred from
observations in some cases \citep{Mennesson2014, Kennedy2015},
although the dust population may be further depleted if the dust has to cross the path of any planets
\citep{Bonsor2018}.
The second is by comet scattering.
That is, an intervening planetary system could scatter planetesimals into the inner regions where they could
sublimate or disintegrate to replenish the warm dust, much as comets in the Solar System replenish the zodiacal
cloud \citep{Nesvorny2010}.

\begin{figure}
\includegraphics[scale=0.66]{wyattfig8}
\caption{Evidence for inward scattering of exocomets in the debris disk of the $\sim 1$\,Gyr F star $\eta$ Crv.
{\it (Left)} ALMA image of its outer debris belt
(and the stellar photospheric emission at the centre),
with inset (bottom right) of the mid-IR emission showing warm dust close to the
star, and annotation to illustrate the origin of that warm dust in exocomet scattering \citep{Marino2017a}.
{\it (Top right)} The profile of CO emission showing this peaks at $\sim 20$\,au close to where the exocomets
may be expected to cross the H$_2$O or CO$_2$ ice lines \citep{Marino2017a}.
{\it (Bottom right)} The spectrum of mid-IR emission can be fitted by a model for the dust composition that
suggests it is water- and carbon-rich \citep{Lisse2012}.
}
\label{fig:8}
\end{figure}

For $\eta$ Crv shown in Fig.~\ref{fig:8}, its $\sim 10$\% excess from dust at
$\sim 1$\,au \citep{Smith2009, Defrere2015}
cannot be explained by steady state collisions in an asteroid belt given its $\sim 1$\,Gyr age \citep{Wyatt2007a}.
In this case an origin in the outer belt seen at 150\,au \citep{Wyatt2005b} is most likely, and is supported by the
composition of the warm dust inferred from its infrared spectrum \citep{Lisse2012}.
The warm dust is too bright to be explained by P-R drag, but inward scattering of comets is the favoured
scenario, particularly following the tentative detection of CO gas at $\sim 20$\,au from the star
\citep{Marino2017a}.
The detection of gas at a location where there is no dust is surprising, but is readily explained if this is the
location where comets start sublimating on their way in before depositing their dust at $\sim 1$\,au.

The requirement to sustain a high rate of comet scattering sets constraints on the intervening planetary system,
e.g. since the comets should not be ejected before reaching the inner regions \citep{Wyatt2017}, but it is found
that a plausible planetary system architecture is one with a chain of $3-30M_\oplus$ planets \citep{Marino2017a}.
More generally, some level of comet scattering is expected in all systems with an outer belt, with the
rate determined by the processes that perturb planetesimals out of their orbits in the belt, and the
architecture of their planetary system \citep{Marino2018a}.
Since these planetesimals may impact any planets and so deliver water or strip their atmospheres,
this has significant interest for those seeking to understand the prevalence of conditions that favour the
development of life.

\section{Conclusions}
\label{s:conclusions}
The Sun is not unique in having either a planetary system or in having belts of planetesimals.
Other stars also have what can be considered as extrasolar analogues to the Kuiper belt.
However, the extent of that analogy depends on how strict the similarities must be.
Certainly other stars have belts of planetesimals at distances from their stars similar to our own
Kuiper belt, but these belts contain orders of magnitude more mass, perhaps comparable to
that in our primordial Kuiper belt.
Since we cannot yet detect a true Kuiper belt analogue around another star
we cannot yet tell how ours compares within the population, except to say that it is not in
the top 20\%.

There are, however, many other similarities to the Kuiper belt in the debris disks observed
around other stars.
For example, several have volatile-rich compositions inferred to be similar to Solar System
comets (from gas observations).
In terms of structures, there is evidence for analogues to the hot and cold classical belt
(from vertically resolved edge-on disks), to the resonant Kuiper belt object population (from clumpy disks),
and to the scattered disk (from haloes of mm-sized grains).
Many disks are seen to be radially narrow with sharp inner and outer edges (from resolved imaging).
It is also inferred that in some systems exocomets are scattered into the inner regions
(from observations of close-in dust and gas).

It is nevertheless worth bearing in mind that these similarities are based in some cases on
a relatively small number of what must inevitably be considered the extremes of the debris disk
population (in that they are bright enough to have been studied in sufficient detail to identify
any similarity).
Also, the similarities that are highlighted in the literature (and in this review)
reflect to some extent our inherently anthropocentric view (e.g., some of the aforementioned
features may have other interpretations).
Moreover, there are also examples of disks that do not fit into a Kuiper belt-like
mold (e.g., the broad disks, or the eccentric rings).

These differences in the extrasolar debris disk population are to be expected, because
just as we have learned by comparing the Solar System to
the exoplanet population, planet formation has a diverse range of outcomes of which
the Solar System is just one example.
By studying and comparing our Kuiper belt to extrasolar debris disks, we have
a different perspective on those outcomes, one which is biased towards processes that
occur in the outer regions of the systems that are hard to access through
exoplanet observations.
What the similarities do show is that, perhaps, we share some dynamical processes
in common with other stars, and by considering the Solar System as one outcome we
may achieve a deeper understanding of the planet formation process and the uniqueness of our own situation.

\bibliographystyle{spbasicHBexo}

Earth and Planetary Astrophysics (astro-ph.EP)

Extrasolar Kuiper Belts


Creation of Quantum-Degenerate Gases of Ytterbium in a Compact 2D-/3D-MOT Setup

Abstract: We report on the first experimental setup based on a 2D-/3D-MOT scheme to create both Bose-Einstein condensates and degenerate Fermi gases of several ytterbium isotopes. Our setup does not require a Zeeman slower and offers the flexibility to simultaneously produce ultracold samples of other atomic species. Furthermore, the extraordinary optical access favors future experiments in optical lattices. A 2D-MOT on the strong 1S0-1P1 transition captures ytterbium directly from a dispenser of atoms and loads a 3D-MOT on the narrow 1S0-3P1 intercombination transition. Subsequently, atoms are transferred to a crossed optical dipole trap and cooled evaporatively to quantum degeneracy.

\section{Introduction}
Owing to their unique properties, alkaline-earth-like atoms have become the subject of intense experimental and theoretical research in recent years. Ultracold gases of these two-electron atoms have been proposed for the realization and study of a variety of novel quantum systems including heavy fermion materials and the Kondo insulator\cite{Foss-Feig_2010A, Foss-Feig_2010B, Foss-Feig_2011}, SU(N)-symmetric Hamiltonians\cite{Gorshkov_2010}, quantum information processing\cite{Daley_2008, Gorshkov_2009} and artificial gauge fields\cite{Gerbier_2010, Dalibard_2011, Szirmai_2011}. Quantum degenerate gases have been reported by only a few groups for isotopes of ytterbium\cite{Takasu_2003, Fukuhara_2007A, Fukuhara_2007B, Fukuhara_2009, Taie_2010, Sugawa_2011, Hansen_2011} (Yb), strontium\cite{Stellmer_2009, deEscobar_2009, DeSalvo_2010, Tey_2010, Stellmer_2010}, and calcium\cite{Kraft_2009, Halder_2012}. However, all of these experiments rely on Zeeman slowers as the initial cooling stage.

From the early days of laser cooling\cite{Pritchard_1987} it is well known that Zeeman\cite{Phillips_1982} or chirp\cite{Ertmer_1985} slowing of a thermal beam enhances the loading rate of a magneto-optical trap (MOT) by many orders of magnitude as compared to loading from a background vapor directly. For alkali metals the concept of a 2D-MOT\cite{Dieckmann_1998} has been developed as a powerful alternative which results in a similar flux of cold atoms, but in a much more compact and versatile setup. Here, atoms are transversely cooled from a background vapor into a well-collimated beam of cold atoms. Because longitudinally slow atoms are cooled preferentially, the resulting velocity distribution peaks at several tens of meters per second allowing them to be captured by the 3D-MOT. Often, a ``pushing'' beam or moving molasses is used to enhance the yield of the cold atom source further. 2D-MOTs are extremely efficient, yet uniquely simple and highly tunable. In particular, a single setup can readily be used to cool multiple elements, whereas multispecies Zeeman slowers require sophisticated designs. 2D-MOTs are therefore inherently well suited for producing and studying ultracold mixtures.\cite{Ospelkaus_2006} Furthermore, high-precision spectroscopy experiments, which are prone to frequency shifts from black-body radiation, can benefit from the fact that there is no direct line of sight from the source of thermal atoms to the cold or ultracold sample.

The traditional concept of a 2D-MOT is limited to elements with sufficiently large vapor pressures near room temperature. Alkaline-earth-like elements, except for mercury\cite{Petersen_2008}, have high melting temperatures $T_l$ (e.g. $T_l = 1097\ \mathrm{K}$ for Yb) and low vapor pressures (e.g. less than $10^{-15}\ \mathrm{Pa}$ at $T\approx300\ \mathrm{K}$ for Yb{\cite{Habermann_1964}}). At room temperature this effectively prevents any 2D-MOT from being loaded from a background gas. An alternative concept has been demonstrated for lithium\cite{Tiecke_2009}, where the 2D-MOT is side-loaded from an oven instead of background vapor and high loading rates of up to $10^9\ \mathrm{s}^{-1}$ have been reported.

In this paper we report on the first successful implementation of a 2D-MOT loaded via a dispenser for Yb and the subsequent generation of quantum degenerate gases of both bosonic ${}^{174}\mathrm{Yb}$ and fermionic ${}^{173}\mathrm{Yb}$. Based on numerical simulations, we have designed a very compact ``glass-cell'' setup for a 2D- and 3D-MOT that allows extraordinarily large optical access for future experiments based on optical lattices.

\section{Experimental Concept}
\subsection{Atomic Properties of Yb}
\begin{figure}
	\centering
	\includegraphics{10_LevelScheme.pdf}
	\caption{Partial Level Scheme of Yb. Only the lowest states in the triplet and singlet manifolds, which are relevant to laser cooling and precision spectroscopy, are shown.\cite{Golub_1988,Takasu_2004,Porsev_2004, Das_2005} The hyperfine structure of the ${}^1P_1$ state is given for the fermionic isotope ${}^{173}\mathrm{Yb}$.}
	\label{fig:levelScheme}
\end{figure}
Yb is a rare-earth element with as many as seven stable isotopes. Its five bosonic isotopes (${}^{168}\mathrm{Yb}$, ${}^{170}\mathrm{Yb}$, ${}^{172}\mathrm{Yb}$, ${}^{174}\mathrm{Yb}$, ${}^{176}\mathrm{Yb}$) are spinless ($I=0$), while the remaining two fermionic isotopes have nuclear spins of $I=1/2$ (${}^{171}\mathrm{Yb}$) and $I=5/2$ (${}^{173}\mathrm{Yb}$). The electronic structure of Yb is dominated by its two valence electrons and hence very similar to alkaline-earth elements. As shown in Fig.~\ref{fig:levelScheme}, its level scheme consists of states with singlet and triplet electronic spin.\cite{Note1} The intercombination transitions connecting singlet and triplet states are only allowed due to imperfect $L$-$S$-coupling and thus are very weak. In particular, the transition from the ground state ${}^1S_0$ to the triplet state ${}^3P_0$ is even forbidden as a single-photon transition for the bosonic isotopes, although small transition amplitudes may be created artificially by magnetically-induced state mixing\cite{Taichenachev_2006} or via multi-photon schemes\cite{Santra_2005, Hong_2005, Ovsiannikov_2007}. In fermionic isotopes, hyperfine interaction yields a non-zero natural linewidth of the transition on the order of $10\ \mathrm{mHz}$. With a lifetime of several seconds\cite{Porsev_2004}, it is metastable on experimentally relevant timescales and gives rise to a series of unique and intriguing features.\cite{Note2}

The ultranarrow ${}^1S_0-{}^3P_0$ transition is not only crucial for novel optical frequency standards\cite{Takamoto_2005, Barber_2006, Poli_2008, Lemke_2009}, its fascinating properties also make it ideal for the spectroscopic analysis of quantum gases. Energy levels of atoms confined to a site of an optical lattice can be probed locally in the optical domain, and due to the coherent nature of the probe even controlled manipulation of quantum states is possible. The wavelength of the lattice acts as a tool to tune the differential AC-Stark shifts of the two states of the spectroscopy transition. In this way, different scenarios can be realized, e.g.~a ``magic'' lattice\cite{Katori_2002} where the main order differential light-shift induced by the trapping potential on the transition is cancelled, thus enabling high-precision spectroscopy. Moreover, state-selective lattices at wavelengths where the polarizability of either ${}^1S_0$ or ${}^3P_0$ vanishes have been proposed\cite{Daley_2008} and allow an unprecedented degree of control of the system.

For laser cooling, two complementary transitions are of interest. The broad principal transition ${}^1S_0\rightarrow{}^1P_1$ allows an extremely strong radiation pressure. In bosonic isotopes, the lack of substructure in the ground state ${}^1S_0$ prevents any Sisyphus cooling on this transition.\cite{Note3} On the other hand, the intercombination transition ${}^1S_0\rightarrow{}^3P_1$ is very narrow and has a \textit{Doppler cooling} limit of merely $T_{D}=4.4\ \mathrm{\mu{}K}$, which provides ideal starting conditions for evaporative cooling. Since the ${}^1S_0$ ground state does not allow magnetic trapping, evaporation and further experimental steps require optical dipole traps. In combination with optical lattices, a setup with wide optical access is therefore desirable for Yb experiments.

\subsection{Side-loaded 2D-MOT}
\begin{figure*}
	\centering
	\includegraphics{20_SchematicView.pdf}
	\caption{Sketch of the 2D-/3D-MOT system. \textbf{(a)} The 3D-MOT using the intercombination transition ${}^1S_0\rightarrow{}^3P_1$ is loaded from a 2D-MOT operated close to the principal transition ${}^1S_0\rightarrow{}^1P_1$ in a separate glass cell. Both cells are mounted to a central vacuum chamber and connected by a dual differential pumping stage. A pushing beam enhances the loading rate of the 3D-MOT. \textbf{(b)} Top-view of \textbf(a). The 2D-MOT is loaded transversely from the beam of atoms emitted by a dispenser.}
	\label{fig:2d3dsketch}
\end{figure*}
\begin{figure*}
	\centering
	\includegraphics{30_2DMotParameters.pdf}
	\caption{Loading rate of a 3D-MOT of ${}^{174}\mathrm{Yb}$ as a function of magnetic field gradient and detuning of the 2D-MOT. Numerical simulations \textbf{(a)} for a capture velocity of $10\ \mathrm{m}/\mathrm{s}$ of the 3D-MOT and corresponding experimental results \textbf{(b)}. Optimal loading is achieved at a gradient of about $55\ \mathrm{G}/\mathrm{cm}$ and a detuning of $\Delta=-1.2\Gamma$.}
	\label{fig:2dmot_sim_exp}
\end{figure*}
\begin{figure}
	\centering
	\includegraphics{35_CaptureVelocity.pdf}
	\caption{Loading rate of a ${}^{174}\mathrm{Yb}$ 3D-MOT as a function of its capture velocity. Results of numerical simulations are shown for magnetic field gradients of $40\ \mathrm{G}/\mathrm{cm}$ ($\blacklozenge$), $55\ \mathrm{G}/\mathrm{cm}$ ($\bullet$) and $70\ \mathrm{G}/\mathrm{cm}$ ($\blacksquare$) at detunings of $-1.4\Gamma$, $-1.2\Gamma$ and $-1.1\Gamma$, respectively, yielding optimal maximum loading rates.}
	\label{fig:sim_capturevelocity}
\end{figure}
With two nearly closed cooling transitions, Yb is a straightforward candidate for a dual MOT scheme. The principal transition ${}^1S_0\rightarrow{}^1P_1$ is used for cooling in the 2D-MOT to achieve a strong cooling force, while the subsequent 3D-MOT operates on the intercombination transition ${}^1S_0\rightarrow{}^3P_1$ that allows for a final temperature in the microkelvin range. The central concept of our 2D-/3D-MOT setup to overcome the low vapor pressure of Yb is to use a side-loaded 2D-MOT\cite{Tiecke_2009}, but use dispensers instead of an oven. As shown in Fig.~\ref{fig:2d3dsketch}, the 2D-MOT is formed by a perpendicular pair of counterpropagating, circularly polarized laser beams and a 2D magnetic quadrupole field generated by rectangular coils. It is loaded transversely from a dispenser that has been installed within a glass cell. Using the narrow intercombination transition, however, the velocity capture range of the 3D-MOT typically is only a few meters per second. As detailed further in section \ref{section:2d3dmot}, we use laser broadening\cite{Kuwamoto_1999} to enhance the loading rate substantially.

To find the optimal 2D-MOT parameters we have performed detailed numerical simulations. Using the geometry as described above, except for the absence of a pushing beam, we numerically integrated the possible trajectories of Yb atoms emitted from the dispenser for different velocities assuming the mean radiation pressure in the 2D-MOT. We assume that Yb atoms emitted from the dispenser stick to the glass cell at their first contact, which is justified by the high desorption temperature. Thus, atoms perform a ``single pass'' through the 2D-MOT, in strong contrast to 2D-MOTs loading from a thermal background vapor. Spontaneous emissions and collisions have been ignored in this calculation. We sampled the velocity distribution in steps of $1\ \mathrm{m}/\mathrm{s}$ per axis up to a maximum of $250\ \mathrm{m}/\mathrm{s}$ and weighted each trajectory with the flux of the source in the associated velocity bin assuming a thermal Maxwell-Boltzmann distribution. The 2D-MOT beams have an elliptic cross-section of $1\times4\ \mathrm{cm}^2$ in diameter and a homogeneous intensity of $I=0.6I_{sat}$ where $I_{sat}=60\ \mathrm{mW}$ is the saturation intensity of the transition ${}^1S_0\rightarrow{}^1P_1$. From these results we then calculated the rate of atoms entering the 3D-MOT volume at a velocity smaller than the capture velocity. Examples of simulation results for a capture velocity $v_c \approx 10\ \mathrm{m}/\mathrm{s}$ are shown in Fig.~\ref{fig:2dmot_sim_exp}a. Evidently, the 2D-MOT efficiently decelerates and cools a substantial amount of atoms emitted by the dispenser, and the resulting beam of atoms is indeed slow enough to subsequently be captured by the 3D-MOT. The loading rate exhibits a distinct maximum at a magnetic quadrupole field gradient of $55\ \mathrm{G}/\mathrm{cm}$ and a detuning $\Delta_{2D} = -1.2\Gamma$. As expected due to the large linewidth of the cooling transition, the optimal gradient is significantly higher than for typical alkali 2D-MOTs. Experimentally obtained loading rates of our 2D-/3D-MOT system for the same range of magnetic field gradients and detunings are shown in Fig.~{\ref{fig:2dmot_sim_exp}b}, and they are in excellent agreement with the simulation. In particular, the optimal set of parameters is reproduced exactly within the resolution of the scans.

Fig.~{\ref{fig:sim_capturevelocity}} shows the expected loading rates as a function of the capture velocity for several 2D-MOT configurations. It illustrates a strong dependence of the optimal 2D-MOT parameters on the capture velocity around $10\ \mathrm{m}/\mathrm{s}$ and a saturation of the loading rate at values beyond some $20\ \mathrm{m}/\mathrm{s}$.

\section{Experimental setup}
\subsection{Apparatus}
\begin{figure}
	\centering
	\includegraphics{40_ApparatusCAD.pdf}
	\caption{Sketch of the CAD model of the 2D-/3D-MOT apparatus.}
	\label{fig:vacuumSystem}
\end{figure}
The vacuum system of our experimental apparatus consists of two separate glass cells that are mounted to a single stainless steel vacuum chamber. Fig.~{\ref{fig:vacuumSystem}} shows a CAD section view of the apparatus. Both cells are made from synthetic fused silica glass, and broadband anti-reflective coatings have been applied to their outer surfaces. While the 2D-MOT is produced in the upper cell, the 3D-MOT is located in the ``science'' cell below, and all subsequent experimental steps take place in this cell as well.

The central chamber connects the 2D-MOT cell to the science cell via two differential pumping stages, equipped with graphite tubes, along its central axis. These are necessary to maintain the low pressure required to achieve sufficiently long lifetimes in the science cell ($p_{3D}\sim10^{-11}\ \mathrm{mbar}$) even when the 2D-MOT is in operation ($p_{2D}\sim10^{-9}\ \mathrm{mbar}$). Each cell is pumped by a $55\ \mathrm{l}/\mathrm{s}$ ion getter pump, and an additional titanium sublimation pump has been installed for the science cell. The chamber itself has been manufactured from solid, non-ferromagnetic steel with internal tubes to connect the pumps and other equipment to the cells. Its compact designs reduces the distance between both MOTs to $30\ \mathrm{cm}$, which is essential for an efficient transfer of atoms from the 2D-MOT into the 3D-MOT.

The 2D-MOT cell ($50\times{}50\times{}120\ \mathrm{mm}^3$) contains two Alvasource\cite{Note4} AS-4-Yb-500-S Yb dispenser sources that generate the atomic beams loading the 2D-MOT when heated. Both dispensers are connected to a high-power feedthrough and have been mounted upright on an insulating ceramic ring with the emitter slit facing towards the center of the cell. To prevent the Yb atoms from coating any critical surface of the 2D-MOT cell, aluminium plates with rectangular apertures have been installed in front of each dispenser. We typically operate the dispensers at currents of about $6.5\ \mathrm{A}$, which yields a sufficiently large flux of Yb. During operation, the pressure measured close to the 2D-MOT cell rises from some $5\times{}10^{-10}\ \mathrm{mbar}$ to a few $10^{-9}\ \mathrm{mbar}$, likely due to outgassing of the dispensers, as emitted Yb directly sticks to the walls. Two identical pairs of external rectangular coils generate the magnetic quadrupole field for the 2D-MOT. The coils are water-cooled and produce the necessary gradients of more than $50\ \mathrm{G}/\mathrm{cm}$ without radiating a significant amount of heat.

The science cell ($26\times36\times77\ \mathrm{mm}^3$) provides excellent optical access to the atomic samples. A pair of Helmholtz coils generates the magnetic quadrupole field for the 3D-MOT, but it has been designed for large a variety of tasks, such as magnetic traps for another species or bias fields for Feshbach resonances in mixtures. It is capable of producing gradients of up to $200\ \mathrm{G}/\mathrm{cm}$ or homogeneous fields of up to $900\ \mathrm{G}$ and does not reduce the optical access significantly.

\subsection{Optical Setup}
The 2D-MOT is formed by two elliptical beams of about $1\ \mathrm{cm}$ by $4\ \mathrm{cm}$ in diameter that are retroreflected in a cat's eye configuration. We use a commercial diode laser system with a tapered amplifier that is frequency-doubled in a bow-tie cavity and provides about $180\ \mathrm{mW}$ of laser power ($I_0 = 1.9 I_{sat}$) to each beam at a wavelength of $399\ \mathrm{nm}$. The laser is offset-locked to a master ECDL using an additional $\mathrm{GaN}$ laser diode and can easily be scanned across a broad frequency range. The master laser itself is stabilized by fluorescence spectroscopy of an atomic beam. Additional light for pushing and imaging beams is generated by an injection-locked slave laser diode of the same type.

The 3D-MOT has been set up in a retroreflected three-beam-configuration (Fig.~{\ref{fig:2d3dsketch}}). A commercial fiber laser is frequency-doubled in a waveguide and provides $6\ \mathrm{mW}$ ($I_0 = 45 I_{sat}$) and $2 \times 16\ \mathrm{mW}$ ($I_0 = 130 I_{sat}$) for the axial and diagonal beams with a diameter of $15\ \mathrm{mm}$. This laser is stabilized via Doppler-free fluorescence spectroscopy of the atomic beam mentioned above, and its inherent narrow linewidth of less than $50\ \mathrm{kHz}$ is particularly useful for operating a MOT at small detunings.

\section{A 2D-/3D-MOT System for Yb}
\subsection{Basic 2D-/3D-MOT Configuration\label{section:2d3dmot}}
\begin{figure}
	\centering
	\includegraphics{50_3DMotBroadening.pdf}
	\caption{3D-MOT loading rate of ${}^{174}\mathrm{Yb}$ measured as a function of the FWHM $\xi$ of the broadened laser spectrum. The high frequency edge of the spectrum is kept at a constant detuning from the ${}^1S_0\rightarrow{}^3P_1$ resonance frequency.}
	\label{fig:motBroadening}
\end{figure}
Fluorescence in the 2D-MOT cell indicates that atoms are captured in a cigar-shaped region along the axis of the magnetic field. However, the emerging flux of cold atoms cannot be captured efficiently in a single-frequency 3D-MOT setup operated close to the ${}^1S_0\rightarrow{}^3P_1$ intercombination transition, because the narrow linewidth addresses only a very small velocity range. The capture velocity can be increased by actively broadening the frequency spectrum of the laser, as has already been reported for Yb setups using Zeeman slowers\cite{Kuwamoto_1999}. This temporarily introduces additional frequency components, preferably further detuned to the red, at the expense of temperature, as detailed in section \ref{section:motperformance}. Due to the low saturation intensity of the transition, there is typically enough laser power to cover a broad frequency range at a reasonable power spectral density. We use the RF signal driving the AOM that controls the power in the 3D-MOT beams to modulate the laser frequency across several megahertz at a modulation frequency of typically $200\ \mathrm{kHz}$. To keep the minimal red detuning roughly at $-3.2\ \mathrm{MHz}$, we also increase the center detuning accordingly. Fig.~\ref{fig:motBroadening} shows the MOT loading rate as a function of the FWHM $\xi$ of the broadened laser spectrum for the bosonic isotope ${}^{174}\mathrm{Yb}$; we attribute the reduced loading rate above $7\ \mathrm{MHz}$ to insufficient spectral power density.

For optimal loading from the 2D-MOT, we typically operate the 3D-MOT at a center detuning of $\Delta_{3D}=-6.7\ \mathrm{MHz}$, a broadening of $\xi = 7\ \mathrm{MHz}$, and a magnetic field gradient of $2\ \mathrm{G}/\mathrm{cm}$. As shown in Fig.~\ref{fig:2dmot_sim_exp}b, the 2D-MOT is operated at a magnetic field gradient of $55\ \mathrm{G}/\mathrm{cm}$ and a detuning of $\Delta_{2D}=-1.2\Gamma$. Because of the large linewidth of the ${}^1S_0\rightarrow{}^1P_1$ transition used for the 2D-MOT, higher laser power may still increase attainable loading rates significantly, however our experimental setup is currently limited by available laser power and power dissipation in the coils (ca. $400\ \mathrm{W}$). When reducing the laser power for the 2D-MOT by a factor of $0.5$, the loading rate of the 3D-MOT is decreased by a factor of $0.4$.

\subsection{Pushing Beam}
\begin{figure}
	\centering
	\includegraphics{60_PushingBeam.pdf}
	\caption{3D-MOT loading rate of ${}^{174}\mathrm{Yb}$ measured as a function of power in the pushing beam. Loading rates are normalized to the value without a pushing beam.}
	\label{fig:pushingBeam}
\end{figure}
The transfer of atoms from the 2D-MOT into the 3D-MOT can be enhanced significantly by a weak additional beam along the axis of the 2D-MOT pushing atoms through the differential pumping stage. We have implemented a single pushing beam as shown schematically in Fig.~{\ref{fig:2d3dsketch}} at a detuning of $\Delta_{p}=-20\ \mathrm{MHz}$ from the ${}^1S_0\rightarrow{}^1P_1$ transition. In Fig.~{\ref{fig:pushingBeam}}, the loading rate of the 3D-MOT is shown as a function of power in the pushing beam. Clearly, a power around $375\ \mathrm{\mu{}W}$ yields an optimal and robust loading rate -- the optimal power depends on its size and the details of its alignment, however, and is given here for reference only. Because of the large linewidth of the ${}^1S_0\rightarrow{}^1P_1$ transition, the radiation pressure of the weak pushing beam exceeds that of the 3D-MOT many times. We therefore align the pushing beam such that it is blocked by the differential pumping tube and does not enter the science cell.

\subsection{MOT Performance\label{section:motperformance}}
For the parameters discussed above, we found a typical value for the 3D-MOT steady-state particle number of $N_{174}^{(0)} = 2\times 10^8$ for $\mathrm{{}^{174}Yb}$, and we typically load $N=6\times{}10^7$ bosonic ${}^{174}\mathrm{Yb}$ atoms within $10\ \mathrm{s}$.

At its optimal loading parameters, the 3D-MOT has a temperature of $T=570\ \mathrm{\mu{}K}$. To reduce the temperature of the sample as close to the Doppler limit of $T_D=4.4\ \mathrm{\mu{}K}$ as possible, the laser-broadening amplitude is gradually reduced to zero within several hundreds of milliseconds, while the center detuning is simultaneously shifted closer to resonance, and the power is reduced to $170\ \mathrm{\mu{}W}$ and $60\ \mathrm{\mu{}W}$ in the diagonal beams and beams along the symmetry axis of the coils, respectively. Furthermore, the MOT is compressed by increasing the field gradient to about $7\ \mathrm{G}/\mathrm{cm}$. With these parameters, the MOT temperature is reduced to typically $T=20\ \mathrm{\mu{}K}$, comparable to temperatures published for other Yb MOTs\cite{Kuwamoto_1999,Hansen_2011}.

\section{Realization of Quantum Degenerate Yb Gases}
\subsection{Crossed Optical Dipole Traps}
For evaporative cooling and trapping of ultracold samples, we have set up a crossed optical dipole trap (ODT) formed by a deep horizontal trap and a vertical trap that creates additional confinement in the crossing region. The horizontal trap is created by a beam of up to $9\ \mathrm{W}$ produced by a solid state laser at $\lambda_{ODT} = 532\ \mathrm{nm}$ and focused down to vertical and horizontal $1/e^2$-waists of $w_y = 18\ \mathrm{\mu{}m}$ and $w_x = 28\ \mathrm{\mu{}m}$, respectively. For the vertical trap, part of the laser power is split off by a partially reflecting mirror and typically up to $1.35\ \mathrm{W}$ are focused to waists of $w_z = 29\ \mathrm{\mu{}m}$ and $w_x = 87\ \mathrm{\mu{}m}$. The powers in each beam are monitored by logarithmically amplified photodiodes and regulated using quartz AOMs that also induce a frequency difference of about $2.5\ \mathrm{MHz}$ between the two beams.

\subsection{Trap Loading}
To transfer atoms into the ODT once the MOT has been loaded, we abruptly switch on the dipole traps and gradually increase the magnetic field gradient to $7\ \mathrm{G}/\mathrm{cm}$ to strongly compress the MOT, while its temperature is reduced as discussed above. Both dipole traps are initially operated at their full powers of $P_h^{(i)} = 9\ \mathrm{W}$ and $P_v^{(i)} = 1.35\ \mathrm{W}$, corresponding to potential depths of $U_h^{(i)}=k_B\times{}625\ \mathrm{\mu{}K}$ and $U_v^{(i)}=k_B\times{}20\ \mathrm{\mu{}K}$ for the horizontal and vertical traps, respectively. Typically $1.2\times{}10^7$ ${}^{174}\mathrm{Yb}$ atoms are trapped in the horizontal ODT at the same temperature measured for the single-frequency MOT. The transfer efficiency strongly depends on the final MOT parameters, in particular the magnetic field gradient.

\subsection{Evaporative Cooling and BEC}
\begin{figure*}
	\centering
	\includegraphics{70_QuantumGases.pdf}
	\caption{Absorption images of quantum degenerate gases of $\mathrm{Yb}$ taken after time of flight and density profiles (circles) integrated along the vertical axis along with fits of appropriate model functions (red lines). \textbf{(a)} Partially condensed Bose gas of ${}^{174}\mathrm{Yb}$ ($20\ \mathrm{ms}$ TOF). The bimodal fit is composed of a Bose distribution and a Thomas-Fermi density distribution of a harmonic trap potential for the condensate fraction. \textbf{(b)} Nearly pure Bose-Einstein condensate of ${}^{174}\mathrm{Yb}$ ($20\ \mathrm{ms}$ TOF) with a fitted Thomas-Fermi density distribution of a harmonic trap potential. \textbf{(c)} Degenerate Fermi gas of ${}^{173}\mathrm{Yb}$ ($8\ \mathrm{ms}$ TOF). The fit shows the distribution of a non-interacting Fermi gas with $T/T_F=0.32(3)$.}
	\label{fig:bec}
\end{figure*}
Evaporative cooling is performed in two stages by lowering the depths of the ODT. In the first stage, atoms are predominantly evaporated from the horizontal trap, but concentrated into the crossed region of the dipole trap as the temperature decreases. Therefor we exponentially ramp the horizontal trap down to $20\ \mathrm{\mu{}K}$ ($P_h = 0.30\ \mathrm{W}$) within $3.5\ \mathrm{s}$, but keep the vertical trap at constant depth.

Subsequently, atoms are evaporated from the crossed region in a second stage. Within $1.5\ \mathrm{s}$ we exponentially reduce the depth of the horizontal trap to nominally $3\ \mathrm{\mu{}K}$ ($P_h = 40\ \mathrm{mW}$) and relax confinement further by lowering the power of the vertical trap to some $P_v = 35\ \mathrm{mW}$ ($0.5\ \mathrm{\mu{}K}$). It is crucial that the confinement along the axis of the horizontal trap provided by the vertical beam is lowered during the final evaporation stage. However, the vertical trap has virtually no effect on the overall trap depth, and its final power can be varied across a large range to control the resulting trap geometry.

The absorption images in Fig.~\ref{fig:bec}a,b show the emergence of a Bose-Einstein condensate of ${}^{174}\mathrm{Yb}$ during the final stages of the evaporation ramp. We routinely produce nearly pure condensates of $N = 1\times{}10^5$ ${}^{174}\mathrm{Yb}$ atoms with no discernible thermal fraction left. The crossed dipole trap at its final parameters has an effective depth of several tens of nanokelvins and trap frequencies $\omega_{x,y,z}=2\pi\times(118, 188, 56)\ \mathrm{Hz}$ determined from dipole and quadrupole oscillations of the condensate.

\subsection{Degenerate Fermi Gases}
We have produced degenerate Fermi gases of the spin-$5/2$ isotope ${}^{173}\mathrm{Yb}$ using a similar sequence as for ${}^{174}\mathrm{Yb}$, because its $s$-wave scattering length\cite{Kitagawa_2008} is favorable for evaporative cooling in a spin mixture.

The 2D-MOT uses the ${}^1S_0\rightarrow{}^1P_1$ ($F=5/2\rightarrow{}F'=7/2$) hyperfine transition of ${}^{173}\mathrm{Yb}$, but its operation is complicated by the proximity of the $F=5/2\rightarrow{}F'=3/2$ hyperfine transition as shown in Fig.~{\ref{fig:levelScheme}}. This results in a strong reduction of the observed 3D-MOT loading rate at detunings close to the hyperfine splitting, but efficient operation of a 2D-MOT is still possible within a narrow window around the same optimal detuning as observed for bosonic ${}^{174}\mathrm{Yb}$. Even more than in the case of ${}^{174}\mathrm{Yb}$, the loading rate strongly depends on the intensity of the 2D-MOT beams, however.

Due to its lower abundance and reduced 2D-MOT performance, the loading time of the ${}^{173}\mathrm{Yb}$ 3D-MOT is significantly larger than for ${}^{174}\mathrm{Yb}$. With an earlier version, i.e. reduced laser powers in the 2D- and 3D-MOT beams, we obtained the following results: A spin mixture of about $2 \times 10^6$ ${}^{173}\mathrm{Yb}$ atoms was transferred from a MOT of $1.5\times10^7$ atoms to the dipole trap. Thermalization is faster than in ${}^{174}\mathrm{Yb}$ due to the enhanced scattering rate, and we used an evaporation ramp of similar shape as described above, but steeper than for the bosons. After $4\ \mathrm{s}$ of evaporation $3\times{}10^4$ atoms at a temperature of $T/T_F = 0.3$ were left. A typical absorption image of a degenerate Fermi gas is shown in Fig.~\ref{fig:bec}c.

\section{Conclusion}
In conclusion, we have presented a novel scheme for the efficient creation of quantum gases of Yb and similar alkaline-earth-like atoms based on transverse loading of a 2D-MOT. Loading rates observed in the subsequent 3D-MOT are in agreement with numerical simulations of the 2D-MOT setup. The compact 2D-/3D-MOT apparatus facilitates laser cooling of additional species to study quantum-degenerate mixtures and provides excellent optical access to the combined 3D-MOT and science cell. Therefore, it allows a wide range of experiments to be conducted. We routinely produce Bose-Einstein condensates of ${}^{174}\mathrm{Yb}$ and degenerate Fermi gases of ${}^{173}\mathrm{Yb}$ with particle numbers well suited for further experiments.

\begin{acknowledgments}
This work has been supported financially by the Deutsche Forschungsgesellschaft within the SFB 925 and GRK 1355 and the European Union within the FET-Open Scheme (iSense). The authors thank Thomas R{\"u}tzel, Hans Kessler, and Jan Carstens for experimental support.
\end{acknowledgments}

Quantum Gases (cond-mat.quant-gas); Quantum Physics (quant-ph)