Double Bubbles Minimize

The classical isoperimetric inequality in R^3 states that the surface of smallest area enclosing a given volume is a sphere. We show that the least area surface enclosing two equal volumes is a double bubble, a surface made of two pieces of round spheres separated by a flat disk, meeting along a single circle at an angle of 120 degrees.Comment: 57 pages, 32 figures. Includes the complete code for a C++ program as described in the article. You can obtain this code by viewing the source of this articl

unWISE tomography of Planck CMB lensing

MB lensing tomography, or the cross-correlation between CMB lensing maps and large-scale structure tracers over a well-defined redshift range, has the potential to map the amplitude and growth of structure over cosmic time, provide some of the most stringent tests of gravity, and break important degeneracies between cosmological parameters. In this work, we use the unWISE galaxy catalog to provide three samples at median redshifts $z \sim 0.6, 1.1$ and 1.5, fully spanning the Dark Energy dominated era, together with the most recent Planck CMB lensing maps. We obtain a combined cross-correlation significance $S/N = 79.3$ over the range of scales $100 < \ell < 1000$. We measure the redshift distribution of unWISE sources by a combination of cross-matching with the COSMOS photometric catalog and cross-correlation with BOSS galaxies and quasars and eBOSS quasars. We also show that magnification bias must be included in our analysis and perform a number of null tests. In a companion paper, we explore the derived cosmological parameters by modeling the non-linearities and propagating the redshift distribution uncertainties.Comment: 51 pages, 22 figures. Comments welcome! Revisions reflect version accepted by JCA

On Galactic density modeling in the presence of dust extinction

Inferences about the spatial density or phase-space structure of stellar populations in the Milky Way require a precise determination of the effective survey volume. The volume observed by surveys such as Gaia or near-infrared spectroscopic surveys, which have good coverage of the Galactic mid-plane region, is highly complex because of the abundant small-scale structure in the three-dimensional interstellar dust extinction. We introduce a novel framework for analyzing the importance of small-scale structure in the extinction. This formalism demonstrates that the spatially-complex effect of extinction on the selection function of a pencil-beam or contiguous sky survey is equivalent to a low-pass filtering of the extinction-affected selection function with the smooth density field. We find that the angular resolution of current 3D extinction maps is sufficient for analyzing Gaia sub-samples of millions of stars. However, the current distance resolution is inadequate and needs to be improved by an order of magnitude, especially in the inner Galaxy. We also present a practical and efficient method for properly taking the effect of extinction into account in analyses of Galactic structure through an effective selection function. We illustrate its use with the selection function of red-clump stars in APOGEE using and comparing a variety of current 3D extinction maps.Comment: Code available at https://github.com/jobovy/mwdust and at https://github.com/jobovy/apogee-map

The Real Estate Investment Trust: State Law Problems

We present a map of the three-dimensional (3D) distribution of dust in the Orion complex. Orion is the closest site of high-mass star formation, making it an excellent laboratory for studying the interstellar medium and star formation. We used data from the Gaia-TGAS catalogue combined with photometry from 2MASS and WISE to get the distances and extinctions of individual stars in the vicinity of the Orion complex. We use a Gaussian process and adopt a non-parametric method to infer the probability distribution function of the dust densities at arbitrary points throughout the region. We map the dust distribution towards different parts of the Orion complex. We find that the distance and depth of the cloud are compatible with other recent works, which show that the method can be applicable to local molecular clouds to map their 3D dust distribution. We also demonstrate the danger of only using colours of stars to derive their extinctions without considering further physical constraints, such as the colour-magnitude diagram (CMD)

The stellar population structure of the Galactic disk

The spatial structure of stellar populations with different chemical abundances in the Milky Way contains a wealth of information on Galactic evolution over cosmic time. We use data on 14,699 red-clump stars from the APOGEE survey, covering 4 kpc <~ R <~ 15 kpc, to determine the structure of mono-abundance populations (MAPs)---stars in narrow bins in [a/Fe] and [Fe/H]---accounting for the complex effects of the APOGEE selection function and the spatially-variable dust obscuration. We determine that all MAPs with enhanced [a/Fe] are centrally concentrated and are well-described as exponentials with a scale length of 2.2+/-0.2 kpc over the whole radial range of the disk. We discover that the surface-density profiles of low-[a/Fe] MAPs are complex: they do not monotonically decrease outwards, but rather display a peak radius ranging from ~5 kpc to ~13 kpc at low [Fe/H]. The extensive radial coverage of the data allows us to measure radial trends in the thickness of each MAP. While high-[a/Fe] MAPs have constant scale heights, low-[a/Fe] MAPs flare. We confirm, now with high-precision abundances, previous results that each MAP contains only a single vertical scale height and that low-[Fe/H], low-[a/Fe] and high-[Fe/H], high-[a/Fe] MAPs have intermediate (h_Z~300 to 600 pc) scale heights that smoothly bridge the traditional thin- and thick-disk divide. That the high-[a/Fe], thick disk components do not flare is strong evidence against their thickness being caused by radial migration. The correspondence between the radial structure and chemical-enrichment age of stellar populations is clear confirmation of the inside-out growth of galactic disks. The details of these relations will constrain the variety of physical conditions under which stars form throughout the MW disk.Comment: Code available at https://github.com/jobovy/apogee-map

An Empirical Relation between Sodium Absorption and Dust Extinction

Dust extinction and reddening are ubiquitous in astronomical observations and are often a major source of systematic uncertainty. We present here a study of the correlation between extinction in the Milky Way and the equivalent width of the NaI D absorption doublet. Our sample includes more than 100 high resolution spectra from the KECK telescopes and nearly a million low resolution spectra from the Sloan Digital Sky Survey (SDSS). We measure the correlation to unprecedented precision, constrain its shape, and derive an empirical relation between these quantities with a dispersion of order 0.15 magnitude in E(B-V). From the shape of the curve of growth we further show that a typical sight line through the Galaxy, as seen within the SDSS footprint, crosses about three dust clouds. We provide a brief guide on how to best estimate extinction to extragalactic sources such as supernovae, using the NaI D absorption feature, under a variety of circumstances.Comment: MNRAS accepte

A 3D Dust Map Based on Gaia, Pan-STARRS 1 and 2MASS

We present a new three-dimensional map of dust reddening, based on Gaia parallaxes and stellar photometry from Pan-STARRS 1 and 2MASS. This map covers the sky north of a declination of -30 degrees, out to a distance of several kiloparsecs. This new map contains three major improvements over our previous work. First, the inclusion of Gaia parallaxes dramatically improves distance estimates to nearby stars. Second, we incorporate a spatial prior that correlates the dust density across nearby sightlines. This produces a smoother map, with more isotropic clouds and smaller distance uncertainties, particularly to clouds within the nearest kiloparsec. Third, we infer the dust density with a distance resolution that is four times finer than in our previous work, to accommodate the improvements in signal-to-noise enabled by the other improvements. As part of this work, we infer the distances, reddenings and types of 799 million stars. We obtain typical reddening uncertainties that are ~30% smaller than those reported in the Gaia DR2 catalog, reflecting the greater number of photometric passbands that enter into our analysis. Our 3D dust map can be accessed at https://doi.org/10.7910/DVN/2EJ9TX or through the Python package "dustmaps," and can be queried interactively at http://argonaut.skymaps.info. Our catalog of stellar parameters can be accessed at https://doi.org/10.7910/DVN/AV9GXO.Comment: 43 pages, 19 figure