We present C and X-band radio observations of the famous utraluminous X-ray source (ULX) Holmberg IX X-1, previously discovered to be associated with an optical emission line nebula several hundred pc in extent. Our recent infrared study of the ULX suggested that a jet could be responsible for the infrared excess detected at the ULX position. The new radio observations, performed using the Karl G. Jansky Very Large Array (VLA) in B-configuration, reveal the presence of a radio counterpart to the nebula with a spectral slope of -0.56 similar to other ULXs. Importantly, we find no evidence for an unresolved radio source associated with the ULX itself, and we set an upper limit on any 5 GHz radio core emission of 6.6 $\mu$Jy ($4.1\times10^{32}$ erg s$^{-1}$). This is 20 times fainter than what we expect if the bubble is energized by a jet. If a jet exists its core component is unlikely to be responsible for the infrared excess unless it is variable. Strong winds which are expected in super-Eddington sources could also play an important role in inflating the radio bubble. We discuss possible interpretations of the radio/optical bubble and we prefer the jet+winds-blown bubble scenario similar to the microquasar SS 433.

We assess the prospects for detecting the moving lens effect using cosmological surveys. The bulk motion of cosmological structure induces a small-scale dipolar temperature anisotropy of the cosmic microwave radiation (CMB), centered around halos and oriented along the transverse velocity field. We introduce a set of optimal filters for this signal, and forecast that a high significance detection can be made with upcoming experiments. We discuss the prospects for reconstructing the bulk transverse velocity field on large scales using matched filters, finding good agreement with previous work using quadratic estimators.

Interferometric visibilities, reflecting the complex correlations between signals recorded at antennas in an interferometric array, carry information about the angular structure of a distant source. While unknown antenna gains in both amplitude and phase can prevent direct interpretation of these measurements, certain combinations of visibilities called closure phases and closure amplitudes are independent of antenna gains and provide a convenient set of robust observables. However, these closure quantities have subtle noise properties and are generally both linearly and statistically dependent. These complications have obstructed the proper use of closure quantities in interferometric analysis, and they have obscured the relationship between analysis with closure quantities and other analysis techniques such as self calibration. We review the statistics of closure quantities, noting common pitfalls that arise when approaching low signal-to-noise due to the nonlinear propagation of statistical errors. We then develop a strategy for isolating and fitting to the independent degrees of freedom captured by the closure quantities through explicit construction of linearly independent sets of quantities along with their noise covariance in the Gaussian limit, valid for moderate signal-to-noise, and we demonstrate that model fits have biased posteriors when this covariance is ignored. Finally, we introduce a unified procedure for fitting to both closure information and partially calibrated visibilities, and we demonstrate both analytically and numerically the direct equivalence of inference based on closure quantities to that based on self calibration of complex visibilities with unconstrained antenna gains.

Exoplanets can evolve significantly between birth and maturity as their atmospheres, orbits, and structures are shaped by their environment. Young planets ($<$1 Gyr) offer the opportunity to probe these sculpting processes. However, most of the known young planets orbit prohibitively faint stars. We present the discovery of two planets transiting HD 63433 (TOI 1726, TIC 130181866), a young Sun-like ($M_*=0.99\pm0.03$) star. Through kinematics, lithium abundance, and rotation, we confirm that HD 63433 is a member of the Ursa Major moving group ($\tau=414\pm23$ Myr). Based on the TESS light curve and updated stellar parameters, the planet radii are $2.15\pm0.10R_\oplus$ and $2.67\pm0.12R_\oplus$, the orbital periods are 7.11 and 20.55 days, and the orbital eccentricities are lower than abut 0.2. Using HARPS-N velocities, we measure the Rossiter-McLaughlin signal of the inner planet, demonstrating the orbit is prograde. Since the host star is bright (V=6.9), both planets are amenable to transmission spectroscopy, radial velocity measurements of their masses, and more precise determination of the stellar obliquity. This system is therefore poised to play an important role in our understanding of planetary system evolution in the first billion years after formation.

KELT-9 b is an ultra hot Jupiter transiting a rapidly rotating, oblate early-A-type star in a polar orbit. We model the effect of rapid stellar rotation on KELT-9 b's transit light curve using photometry from the Transiting Exoplanet Survey Satellite (\tess) to constrain the planet's true spin-orbit angle and to explore how KELT-9 b may be influenced by stellar gravity darkening. We constrain the host star's equatorial radius to be $1.089\pm0.017$ times as large as its polar radius and its local surface brightness to vary by $\sim38$\% between its hot poles and cooler equator. We model the stellar oblateness and surface brightness gradient and find that it causes the transit light curve to lack the usual symmetry around the time of minimum light. We take advantage of the light curve asymmetry to constrain KELT-9 b's true spin orbit angle (${87^\circ}^{+10^\circ}_{-11^\circ}$), agreeing with \citet{gaudi2017giant} that KELT-9 b is in a nearly polar orbit. We also apply a gravity darkening correction to the spectral energy distribution model from \citet{gaudi2017giant} and find that accounting for rapid rotation gives a better fit to available spectroscopy and yields a more reliable estimate for the star's polar effective temperature.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

The Event Horizon Telescope image of the supermassive black hole in the galaxy M87 is dominated by a bright, unresolved ring. General relativity predicts that embedded within this image lies a thin "photon ring," which is composed of an infinite sequence of self-similar subrings that are indexed by the number of photon orbits around the black hole. The subrings approach the edge of the black hole "shadow," becoming exponentially narrower but weaker with increasing orbit number, with seemingly negligible contributions from high order subrings. Here, we show that these subrings produce strong and universal signatures on long interferometric baselines. These signatures offer the possibility of precise measurements of black hole mass and spin, as well as tests of general relativity, using only a sparse interferometric array.

We report on the first detection of pulsed radio emission from a radio pulsar with the ALMA telescope. The detection was made in the Band-3 frequency range (85-101 GHz) using ALMA in the phased-array mode developed for VLBI observations. A software pipeline has been implemented to enable a regular pulsar observing mode in the future. We describe the pipeline and demonstrate the capability of ALMA to perform pulsar timing and searching. We also measure the flux density and polarization properties of the Vela pulsar (PSR J0835$-$4510) at mm-wavelengths, providing the first polarimetric study of any ordinary pulsar at frequencies above 32 GHz. Finally, we discuss the lessons learned from the Vela observations for future pulsar studies with ALMA, particularly for searches near the supermassive black hole in the Galactic Center, and the potential of using pulsars for polarization calibration of ALMA.

WASP-33b is a retrograde hot Jupiter with a period of 1.2 days orbiting around a rapidly rotating and pulsating A-type star. A previous study found that the transit chord of WASP-33b had changed slightly from 2008 to 2014 based on Doppler tomographic measurements. They attributed the change to orbital precession caused by the non-zero oblateness of the host star and the misaligned orbit. We aim to confirm and more precisely model the precession behavior using additional Doppler tomographic data of WASP-33b obtained with the High Dispersion Spectrograph on the 8.2m Subaru telescope in 2011, as well as the datasets used in the previous study. Using equations of a long-term orbital precession, we constrain the stellar gravitational quadrupole moment $J_{2}=(9.14\pm 0.51)\times 10^{-5}$ and the angle between the stellar spin axis and the line of sight $i_{\star}=96^{+10}_{-14}$ deg. These values update that the host star is more spherical and viewed more equator than the previous study. We also estimate that the precession period is $\sim$840 years. We also find that the precession amplitude of WASP-33b is $\sim$67 deg and WASP-33b transits in front of the host star for only $\sim$20\% of the whole precession period.

CMB-HD is a proposed ultra-deep (0.5 uk-arcmin), high-resolution (15 arcseconds) millimeter-wave survey over half the sky that would answer many outstanding questions in both fundamental physics of the Universe and astrophysics. This survey would be delivered in 7.5 years of observing 20,000 square degrees, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors.

We perform a joint analysis of the counts and weak lensing signal of redMaPPer clusters selected from the Dark Energy Survey (DES) Year 1 dataset. Our analysis uses the same shear and source photometric redshifts estimates as were used in the DES combined probes analysis. Our analysis results in surprisingly low values for $S_8 =\sigma_8(\Omega_{\rm m}/0.3)^{0.5}= 0.65\pm 0.04$, driven by a low matter density parameter, $\Omega_{\rm m}=0.179^{+0.031}_{-0.038}$, with $\sigma_8-\Omega_{\rm m}$ posteriors in $2.4\sigma$ tension with the DES Y1 3x2pt results, and in $5.6\sigma$ with the Planck CMB analysis. These results include the impact of post-unblinding changes to the analysis, which did not improve the level of consistency with other data sets compared to the results obtained at the unblinding. The fact that multiple cosmological probes (supernovae, baryon acoustic oscillations, cosmic shear, galaxy clustering and CMB anisotropies), and other galaxy cluster analyses all favor significantly higher matter densities suggests the presence of systematic errors in the data or an incomplete modeling of the relevant physics. Cross checks with X-ray and microwave data, as well as independent constraints on the observable--mass relation from SZ selected clusters, suggest that the discrepancy resides in our modeling of the weak lensing signal rather than the cluster abundance. Repeating our analysis using a higher richness threshold ($\lambda \ge 30$) significantly reduces the tension with other probes, and points to one or more richness-dependent effects not captured by our model.

We present a catalog of 316 trans-Neptunian bodies detected by the Dark Energy Survey (DES). These objects include 245 discoveries by DES (139 not previously published) detected in $\approx 60,000$ exposures from the first four seasons of the survey ("Y4" data). The survey covers a contiguous 5000 deg$^2$ of the southern sky in the $grizY$ optical/NIR filter set, with a typical TNO in this part of the sky being targeted by $25-30$ Y4 exposures. We describe the processes for detection of transient sources and the linkage into TNO orbits, which are made challenging by the absence of the few-hour repeat observations employed by TNO-optimized surveys. We also describe the procedures for determining detection efficiencies vs. magnitude and estimating rates of false-positive linkages. This work presents all TNOs which were detected on $\ge 6$ unique nights in the Y4 data and pass a "sub-threshold confirmation" test wherein we demand the the object be detectable in a stack of the individual images in which the orbit indicates an object should be present, but was not detected. This eliminates false positives and yields TNO detections complete to $r\lesssim 23.3$ mag with virtually no dependence on orbital properties for bound TNOs at distance $30\,{\rm AU}<d<2500\,{\rm AU}.$ The final DES TNO catalog is expected to yield $>0.3$ mag more depth, and arcs of $>4$ years for nearly all detections.

We present a catalog of 316 trans-Neptunian bodies detected by the Dark Energy Survey (DES). These objects include 245 discoveries by DES (139 not previously published) detected in $\approx 60,000$ exposures from the first four seasons of the survey ("Y4" data). The survey covers a contiguous 5000 deg$^2$ of the southern sky in the $grizY$ optical/NIR filter set, with a typical TNO in this part of the sky being targeted by $25-30$ Y4 exposures. We describe the processes for detection of transient sources and the linkage into TNO orbits, which are made challenging by the absence of the few-hour repeat observations employed by TNO-optimized surveys. We also describe the procedures for determining detection efficiencies vs. magnitude and estimating rates of false-positive linkages. This work presents all TNOs which were detected on $\ge 6$ unique nights in the Y4 data and pass a "sub-threshold confirmation" test wherein we demand the the object be detectable in a stack of the individual images in which the orbit indicates an object should be present, but was not detected. This eliminates false positives and yields TNO detections complete to $r\lesssim 23.3$ mag with virtually no dependence on orbital properties for bound TNOs at distance $30\,{\rm AU}<d<2500\,{\rm AU}.$ The final DES TNO catalog is expected to yield $>0.3$ mag more depth, and arcs of $>4$ years for nearly all detections.

We report the discovery of an intermediate-mass transiting brown dwarf, TOI-503b, from the TESS mission. TOI-503b is the first brown dwarf discovered by TESS and orbits a metallic-line A-type star with a period of $P=3.6772 \pm 0.0001$ days. The light curve from TESS indicates that TOI-503b transits its host star in a grazing manner, which limits the precision with which we measure the brown dwarf's radius ($R_b = 1.34^{+0.26}_{-0.15} R_J$). We obtained high-resolution spectroscopic observations with the FIES, Ond\v{r}ejov, PARAS, Tautenburg, and TRES spectrographs and measured the mass of TOI-503b to be $M_b = 53.7 \pm 1.2 M_J$. The host star has a mass of $M_\star = 1.80 \pm 0.06 M_\odot$, a radius of $R_\star = 1.70 \pm 0.05 R_\odot$, an effective temperature of $T_{\rm eff} = 7650 \pm 160$K, and a relatively high metallicity of $0.61\pm 0.07$ dex. We used stellar isochrones to derive the age of the system to be $\sim$180 Myr, which places its age between that of RIK 72b (a $\sim$10 Myr old brown dwarf in the Upper Scorpius stellar association) and AD 3116b (a $\sim$600 Myr old brown dwarf in the Praesepe cluster). We argue that this brown dwarf formed in-situ, based on the young age of the system and the long circularization timescale for this brown dwarf around its host star. TOI-503b joins a growing number of known short-period, intermediate-mass brown dwarfs orbiting main sequence stars, and is the second such brown dwarf known to transit an A star, after HATS-70b. With the growth in the population in this regime, the driest region in the brown dwarf desert ($35-55 M_J \sin{i}$) is reforesting and its mass range shrinking.

Black hole images present an annular region of enhanced brightness. In the absence of propagation effects, this "photon ring" has universal features that are completely governed by general relativity and independent of the details of the emission. Here, we show that the polarimetric image of a black hole also displays universal properties. In particular, the photon ring exhibits a self-similar pattern of polarization that encodes the black hole spin. We explore the corresponding universal polarimetric signatures of the photon ring on long interferometric baselines, and propose a method for measuring the black hole spin using a sparse interferometric array. These signatures could enable spin measurements of the supermassive black hole in M87, as well as precision tests of general relativity in the strong field regime, via a future extension of the Event Horizon Telescope to space.

Transiting planets in nearby young clusters offer the opportunity to study the atmospheres and dynamics of planets during their formative years. To this end, we focused on K2-25b -- a close-in ($P$=3.48 days), Neptune-sized exoplanet orbiting a M4.5 dwarf in the 650Myr Hyades cluster. We combined photometric observations of K2-25 covering a total of 44 transits and spanning >2 yr, drawn from a mix of space-based telescopes (Spitzer Space Telescope and K2) and ground-based facilities (Las Cumbres Observatory Global Telescope network and MEarth). The transit photometry spanned 0.6--4.5$\mu$m, which enabled our study of K2-25b's transmission spectrum. We combined and fit each dataset at a common wavelength within a Markov Chain Monte Carlo framework, yielding consistent planet parameters. The resulting transit depths ruled out a solar-composition atmosphere for K2-25b for the range of expected planetary masses and equilibrium temperature at a $>4\sigma$ confidence level, and are consistent with a flat transmission spectrum. Mass constraints and transit observations at a finer grid of wavelengths (e.g., from the Hubble Space Telescope) are needed to make more definitive statements about the presence of clouds or an atmosphere of high mean molecular weight. Our precise measurements of K2-25b's transit duration also enabled new constraints on the eccentricity of K2-25's orbit. We find K2-25b's orbit to be eccentric ($e>0.20$) for all reasonable stellar densities and independent of the observation wavelength or instrument. The high eccentricity is suggestive of a complex dynamical history and motivates future searches for additional planets or stellar companions.

We present the first detailed elemental abundances in the ultra-faint Magellanic satellite galaxies Carina II (Car II) and Carina III (Car III). With high-resolution Magellan/MIKE spectroscopy, we determined abundances of nine stars in Car II including the first abundances of an RR Lyrae star in an ultra-faint dwarf galaxy; and two stars in Car III. The chemical abundances demonstrate that both systems are clearly galaxies and not globular clusters. The stars in these galaxies mostly display abundance trends matching those of other similarly faint dwarf galaxies: enhanced but declining [alpha/Fe] ratios, iron-peak elements matching the stellar halo, and unusually low neutron-capture element abundances. One star displays a low outlying [Sc/Fe] = -1.0. We detect a large Ba scatter in Car II, likely due to inhomogeneous enrichment by low-mass AGB star winds. The most striking abundance trend is for [Mg/Ca] in Car II, which decreases from +0.4 to -0.4 and indicates clear variation in the initial progenitor masses of enriching core-collapse supernovae. So far, the only ultra-faint dwarf galaxies displaying a similar [Mg/Ca] trend are likely satellites of the Large Magellanic Cloud. We find two stars with [Fe/H] < -3.5, whose abundances likely trace the first generation of metal-free Population III stars and are well-fit by Population III core-collapse supernova yields. An appendix describes our new abundance uncertainty analysis that propagates line-by-line stellar parameter uncertainties.

We present the discoveries of KELT-25b (TIC 65412605, TOI-626.01) and KELT-26b (TIC 160708862, TOI-1337.01), two transiting companions orbiting relatively bright, early A-stars. The transit signals were initially detected by the KELT survey, and subsequently confirmed by \textit{TESS} photometry. KELT-25b is on a 4.40-day orbit around the V = 9.66 star CD-24 5016 ($T_{\rm eff} = 8280^{+440}_{-180}$ K, $M_{\star}$ = $2.18^{+0.12}_{-0.11}$ $M_{\odot}$), while KELT-26b is on a 3.34-day orbit around the V = 9.95 star HD 134004 ($T_{\rm eff}$ =$8640^{+500}_{-240}$ K, $M_{\star}$ = $1.93^{+0.14}_{-0.16}$ $M_{\odot}$), which is likely an Am star. We have confirmed the sub-stellar nature of both companions through detailed characterization of each system using ground-based and \textit{TESS} photometry, radial velocity measurements, Doppler Tomography, and high-resolution imaging. For KELT-25, we determine a companion radius of $R_{\rm P}$ = $1.64^{+0.039}_{-0.043}$ $R_{\rm J}$, and a 3-sigma upper limit on the companion's mass of $\sim64~M_{\rm J}$. For KELT-26b, we infer a planetary mass and radius of $M_{\rm P}$ = $1.41^{+0.43}_{-0.51}$ $M_{\rm J}$ and $R_{\rm P}$ = $1.940^{+0.060}_{-0.058}$ $R_{\rm J}$. From Doppler Tomographic observations, we find KELT-26b to reside in a highly misaligned orbit. This conclusion is weakly corroborated by a subtle asymmetry in the transit light curve from the \textit{TESS} data. KELT-25b appears to be in a well-aligned, prograde orbit, and the system is likely a member of a cluster or moving group.

Multi-messenger astrophysics is a fast-growing, interdisciplinary field that combines data, which vary in volume and speed of data processing, from many different instruments that probe the Universe using different cosmic messengers: electromagnetic waves, cosmic rays, gravitational waves and neutrinos. In this Expert Recommendation, we review the key challenges of real-time observations of gravitational wave sources and their electromagnetic and astroparticle counterparts, and make a number of recommendations to maximize their potential for scientific discovery. These recommendations refer to the design of scalable and computationally efficient machine learning algorithms; the cyber-infrastructure to numerically simulate astrophysical sources, and to process and interpret multi-messenger astrophysics data; the management of gravitational wave detections to trigger real-time alerts for electromagnetic and astroparticle follow-ups; a vision to harness future developments of machine learning and cyber-infrastructure resources to cope with the big-data requirements; and the need to build a community of experts to realize the goals of multi-messenger astrophysics.

The Neptune desert is a feature seen in the radius-mass-period plane, whereby a notable dearth of short period, Neptune-like planets is found. Here we report the {\it TESS} discovery of a new short-period planet in the Neptune desert, orbiting the G-type dwarf TYC\,8003-1117-1 (TOI-132). {\it TESS} photometry shows transit-like dips at the level of $\sim$1400 ppm occurring every $\sim$2.11 days. High-precision radial velocity follow-up with HARPS confirmed the planetary nature of the transit signal and provided a semi-amplitude radial velocity variation of $\sim$11.5 m s$^{-1}$, which, when combined with the stellar mass of $0.97\pm0.06$ $M_{\odot}$, provides a planetary mass of 22.83$^{+1.81}_{-1.80}$ $M_{\oplus}$. Modeling the {\it TESS} high-quality light curve returns a planet radius of 3.43$^{+0.13}_{-0.14}$ $R_{\oplus}$, and therefore the planet bulk density is found to be 3.11$^{+0.44}_{-0.450}$ g cm$^{-3}$. Planet structure models suggest that the bulk of the planet mass is in the form of a rocky core, with an atmospheric mass fraction of 4.3$^{+1.2}_{-2.3}$\%. TOI-132 b is a {\it TESS} Level 1 Science Requirement candidate, and therefore priority follow-up will allow the search for additional planets in the system, whilst helping to constrain low-mass planet formation and evolution models, particularly valuable for better understanding the Neptune desert.

We present a simulated cosmology analysis using the second and third moments of the weak lensing mass (convergence) maps. The analysis is geared towards the third year (Y3) data from the Dark Energy Survey (DES), but the methodology can be applied to other weak lensing data sets. The second moment, or variances, of the convergence as a function of smoothing scale contains information similar to standard shear 2-point statistics. The third moment, or the skewness, contains additional non-Gaussian information. We present the formalism for obtaining the convergence maps from the measured shear and for obtaining the second and third moments of these maps given partial sky coverage. We estimate the covariance matrix from a large suite of numerical simulations. We test our pipeline through a simulated likelihood analyses varying 5 cosmological parameters and 10 nuisance parameters. Our forecast shows that the combination of second and third moments provides a 1.5 percent constraint on $S_8 \equiv \sigma_8 (\Omega_{\rm m}/0.3)^{0.5}$ for DES Y3 data. This is 20 percent better than an analysis using a simulated DES Y3 shear 2-point statistics, owing to the non-Gaussian information captured by the inclusion of higher-order statistics. The methodology developed here can be applied to current and future weak lensing datasets to make use of information from non-Gaussianity in the cosmic density field and to improve constraints on cosmological parameters.

High-resolution imaging of supermassive black holes is now possible, with new applications to testing general relativity and horizon-scale accretion and relativistic jet formation processes. Over the coming decade, the EHT will propose to add new strategically placed VLBI elements operating at 1.3mm and 0.87mm wavelength. In parallel, development of next-generation backend instrumentation, coupled with high throughput correlation architectures, will boost sensitivity, allowing the new stations to be of modest collecting area while still improving imaging fidelity and angular resolution. The goal of these efforts is to move from imaging static horizon scale structure to dynamic reconstructions that capture the processes of accretion and jet launching in near real time.

We present constraints on extensions of the minimal cosmological models dominated by dark matter and dark energy, $\Lambda$CDM and $w$CDM, by using a combined analysis of galaxy clustering and weak gravitational lensing from the first-year data of the Dark Energy Survey (DES Y1) in combination with external data. We consider four extensions of the minimal dark energy-dominated scenarios: 1) nonzero curvature $\Omega_k$, 2) number of relativistic species $N_{\rm eff}$ different from the standard value of 3.046, 3) time-varying equation-of-state of dark energy described by the parameters $w_0$ and $w_a$ (alternatively quoted by the values at the pivot redshift, $w_p$, and $w_a$), and 4) modified gravity described by the parameters $\mu_0$ and $\Sigma_0$ that modify the metric potentials. We also consider external information from Planck CMB measurements; BAO measurements from SDSS, 6dF, and BOSS; RSD measurements from BOSS; and SNIa information from the Pantheon compilation. Constraints on curvature and the number of relativistic species are dominated by the external data; when these are combined with DES Y1, we find $\Omega_k=0.0020^{+0.0037}_{-0.0032}$ at the 68% confidence level, and $N_{\rm eff}<3.28\, (3.55)$ at 68% (95%) confidence. For the time-varying equation-of-state, we find the pivot value $(w_p, w_a)=(-0.91^{+0.19}_{-0.23}, -0.57^{+0.93}_{-1.11})$ at pivot redshift $z_p=0.27$ from DES alone, and $(w_p, w_a)=(-1.01^{+0.04}_{-0.04}, -0.28^{+0.37}_{-0.48})$ at $z_p=0.20$ from DES Y1 combined with external data; in either case we find no evidence for the temporal variation of the equation of state. For modified gravity, we find the present-day value of the relevant parameters to be $\Sigma_0= 0.43^{+0.28}_{-0.29}$ from DES Y1 alone, and $(\Sigma_0, \mu_0)=(0.06^{+0.08}_{-0.07}, -0.11^{+0.42}_{-0.46})$ from DES Y1 combined with external data, consistent with predictions from GR.

We present constraints on extensions of the minimal cosmological models dominated by dark matter and dark energy, $\Lambda$CDM and $w$CDM, by using a combined analysis of galaxy clustering and weak gravitational lensing from the first-year data of the Dark Energy Survey (DES Y1) in combination with external data. We consider four extensions of the minimal dark energy-dominated scenarios: 1) nonzero curvature $\Omega_k$, 2) number of relativistic species $N_{\rm eff}$ different from the standard value of 3.046, 3) time-varying equation-of-state of dark energy described by the parameters $w_0$ and $w_a$ (alternatively quoted by the values at the pivot redshift, $w_p$, and $w_a$), and 4) modified gravity described by the parameters $\mu_0$ and $\Sigma_0$ that modify the metric potentials. We also consider external information from Planck CMB measurements; BAO measurements from SDSS, 6dF, and BOSS; RSD measurements from BOSS; and SNIa information from the Pantheon compilation. Constraints on curvature and the number of relativistic species are dominated by the external data; when these are combined with DES Y1, we find $\Omega_k=0.0020^{+0.0037}_{-0.0032}$ at the 68% confidence level, and $N_{\rm eff}<3.28\, (3.55)$ at 68% (95%) confidence. For the time-varying equation-of-state, we find the pivot value $(w_p, w_a)=(-0.91^{+0.19}_{-0.23}, -0.57^{+0.93}_{-1.11})$ at pivot redshift $z_p=0.27$ from DES alone, and $(w_p, w_a)=(-1.01^{+0.04}_{-0.04}, -0.28^{+0.37}_{-0.48})$ at $z_p=0.20$ from DES Y1 combined with external data; in either case we find no evidence for the temporal variation of the equation of state. For modified gravity, we find the present-day value of the relevant parameters to be $\Sigma_0= 0.43^{+0.28}_{-0.29}$ from DES Y1 alone, and $(\Sigma_0, \mu_0)=(0.06^{+0.08}_{-0.07}, -0.11^{+0.42}_{-0.46})$ from DES Y1 combined with external data, consistent with predictions from GR.

We describe the design and assembly of the LUX-ZEPLIN experiment, a direct detection search for cosmic WIMP dark matter particles. The centerpiece of the experiment is a large liquid xenon time projection chamber sensitive to low energy nuclear recoils. Rejection of backgrounds is enhanced by a Xe skin veto detector and by a liquid scintillator Outer Detector loaded with gadolinium for efficient neutron capture and tagging. LZ is located in the Davis Cavern at the 4850' level of the Sanford Underground Research Facility in Lead, South Dakota, USA. We describe the major subsystems of the experiment and its key design features and requirements.

We explore a simple spherical model of optically thin accretion on a Schwarzschild black hole, and study the properties of the image as seen by a distant observer. We show that a dark circular region in the center --- a shadow --- is always present. The outer edge of the shadow is located at the photon ring radius $b_{\rm ph} \equiv \sqrt{27}r_g$, where $r_g=GM/c^2$ is the gravitational radius of the accreting mass $M$. The location of the shadow edge is independent of the inner radius at which the accreting gas stops radiating. The size of the observed shadow is thus a signature of the spacetime geometry and it is hardly influenced by accretion details. We briefly discuss the relevance of these results for the Event Horizon Telescope image of the supermassive black hole in M87.

We present the discovery of KELT J072709+072007 (HD 58730), a very low mass ratio ($q \equiv M_2/M_1 \approx 0.08$) eclipsing binary (EB) identified by the Kilodegree Extremely Little Telescope (KELT) survey. We present the discovery light curve and perform a global analysis of four high-precision ground-based light curves, the Transiting Exoplanets Survey Satellite (TESS) light curve, radial velocity (RV) measurements, Doppler Tomography (DT) measurements, and the broad-band spectral energy distribution (SED). Results from the global analysis are consistent with a fully convective ($M_2 = 0.253^{+0.021}_{-0.017}\ M_{\odot})$ M star transiting a late-B primary ($M_1 = 3.348^{+0.057}_{-0.082}\ M_{\odot};\ T_{\rm eff,1} = 12000^{+580}_{-530}\ {\rm K}$). We infer that the system is younger than 272 Myr ($3\sigma$), and the M star mass and radius are consistent with values from a pre-main sequence isochrone of comparable age. We separately and analytically fit for the variability in the out-of-eclipse TESS phase curve, finding good agreement between these results and those from the global analysis. Such systems are valuable for testing theories of binary star formation and understanding how the environment of a star in a close-but-detached binary affects its physical properties. In particular, we examine how a star's properties in such a binary might differ from the properties it would have in isolation.

We investigate the nonlinear dynamics of cold atom systems that can in principle serve as quantum simulators of false vacuum decay. The analog false vacuum manifests as a metastable vacuum state for the relative phase in a two-species Bose-Einstein condensate (BEC), induced by a driven periodic coupling between the two species. In the appropriate low energy limit, the evolution of the relative phase is approximately governed by a relativistic wave equation exhibiting true and false vacuum configurations. In previous work, a linear stability analysis identified exponentially growing short-wavelength modes driven by the time-dependent coupling. These modes threaten to destabilize the analog false vacuum. Here, we employ numerical simulations of the coupled Gross-Pitaevski equations (GPEs) to determine the non-linear evolution of these linearly unstable modes. We find that unless a physical mechanism modifies the GPE on short length scales, the analog false vacuum is indeed destabilized. We briefly discuss various physically expected corrections to the GPEs that may act to remove the exponentially unstable modes. To investigate the resulting dynamics in cases where such a removal mechanism exists, we implement a hard UV cutoff that excludes the unstable modes as a simple model for these corrections. We use this to study the range of phenomena arising from such a system. In particular, we show that by modulating the strength of the time-dependent coupling, it is possible to observe the crossover between a second and first order phase transition out of the false vacuum.

Operating 1.5 million km from Earth at the Sun-Earth L2 Lagrange point, the Origins Space Telescope equipped with a slightly modified version of its HERO heterodyne instrument could function as a uniquely valuable node in a VLBI network. The unprecedented angular resolution resulting from the combination of Origins with existing ground-based millimeter/submillimeter telescope arrays would increase the number of spatially resolvable black holes by a factor of a million, permit the study of these black holes across all of cosmic history, and enable new tests of general relativity by unveiling the photon ring substructure in the nearest black holes.

We present the discovery of KELT-24 b, a massive hot Jupiter orbiting a bright (V=8.3 mag, K=7.2 mag) young F-star with a period of 5.6 days. The host star, KELT-24 (HD 93148), has a $T_{\rm eff}$ =$6509^{+50}_{-49}$ K, a mass of $M_{*}$ = $1.460^{+0.055}_{-0.059}$ $M_{\odot}$, radius of $R_{*}$ = $1.506\pm0.022$ $R_{\odot}$, and an age of $0.78^{+0.61}_{-0.42}$ Gyr. Its planetary companion (KELT-24 b) has a radius of $R_{\rm P}$ = $1.272\pm0.021$ $R_{\rm J}$, a mass of $M_{\rm P}$ = $5.18^{+0.21}_{-0.22}$ $M_{\rm J}$, and from Doppler tomographic observations, we find that the planet's orbit is well-aligned to its host star's projected spin axis ($\lambda$ = $2.6^{+5.1}_{-3.6}$). The young age estimated for KELT-24 suggests that it only recently started to evolve from the zero-age main sequence. KELT-24 is the brightest star known to host a transiting giant planet with a period between 5 and 10 days. Although the circularization timescale is much longer than the age of the system, we do not detect a large eccentricity or significant misalignment that is expected from dynamical migration. The brightness of its host star and its moderate surface gravity make KELT-24b an intriguing target for detailed atmospheric characterization through spectroscopic emission measurements since it would bridge the current literature results that have primarily focused on lower mass hot Jupiters and a few brown dwarfs.

The Event Horizon Telescope (EHT) recently produced the first horizon-scale image of a supermassive black hole. Expanding the array to include a 3-meter space telescope operating at >200 GHz enables mass measurements of many black holes, movies of black hole accretion flows, and new tests of general relativity that are impossible from the ground.

Based on two decades of radial velocity (RV) observations using Keck/HIRES and McDonald/Tull, and more recent observations using the Automated Planet Finder, we found that the nearby star HR 5183 (HD 120066) hosts a 3$M_J$ minimum mass planet with an orbital period of $74^{+43}_{-22}$ years. The orbit is highly eccentric (e$\simeq$0.84), shuttling the planet from within the orbit of Jupiter to beyond the orbit of Neptune. Our careful survey design enabled high cadence observations before, during, and after the planet's periastron passage, yielding precise orbital parameter constraints. We searched for stellar or planetary companions that could have excited the planet's eccentricity, but found no candidates, potentially implying that the perturber was ejected from the system. We did identify a bound stellar companion more than 15,000 au from the primary, but reasoned that it is currently too widely separated to have an appreciable effect on HR 5183 b. Because HR 5183 b's wide orbit takes it more than 30 au (1") from its star, we also explored the potential of complimentary studies with direct imaging or stellar astrometry. We found that a Gaia detection is very likely, and that imaging at 10 $\mu$m is a promising avenue. This discovery highlights the value of long-baseline RV surveys for discovering and characterizing long-period, eccentric Jovian planets. This population may offer important insights into the dynamical evolution of planetary systems containing multiple massive planets.

We present the design and implementation of an automated data calibration and reduction pipeline for very-long-baseline interferometric (VLBI) observations taken at millimeter wavelengths. These short radio-wavelengths provide the best imaging resolution available from ground-based VLBI networks such as the Event Horizon Telescope (EHT) and the Global Millimeter VLBI Array (GMVA), but require specialized processing due to the strong effects from atmospheric opacity and turbulence as well as the heterogeneous nature of existing global arrays. The pipeline builds upon a calibration suite (HOPS) originally designed for precision geodetic VLBI. To support the reduction of data for astronomical observations, we have developed an additional framework for global phase and amplitude calibration which provides output in a standard data format for astronomical imaging and analysis. We test the pipeline on observations taken at 3.5 mm (86 GHz) by the GMVA joined by the phased Atacama Large Millimeter/submillimeter Array in April 2017, and demonstrate the benefits from the specialized processing of high frequency VLBI data with respect to classical analysis techniques.

Compensated isocurvature perturbations (CIPs) are modulations of the relative baryon and dark matter density that leave the total matter density constant. The best current constraints from the primary cosmic microwave background (CMB) are consistent with CIPs some two orders of magnitude larger in amplitude than adiabatic perturbations, suggesting that there may be a huge gap in our knowledge of the early Universe. However, it was recently suggested by Barreira~et.~al. that CIPs which are correlated with the primordial curvature perturbation, as arises in some versions of the curvaton model, lead to a new observable: scale dependent galaxy bias. Combining a galaxy survey with an unbiased tracer of the density field facilitates a measurement of the amplitude of correlated CIPs that is free from cosmic variance, the main limitation on constraints from the primary CMB. Among the most promising tracers to use for this purpose is the remote dipole field, reconstructed using the technique of kinetic Sunyaev Zel'dovich (kSZ) tomography. In this paper, we evaluate the detection significance on the amplitude of correlated CIPs possible with next-generation CMB and galaxy surveys using kSZ tomography. Our analysis includes all relativistic contributions to the observed galaxy number counts and allows for both CIPs and primordial non-Gaussianity, which also gives rise to a scale dependent galaxy bias. We find that kSZ tomography can probe CIPs of comparable amplitude to the adiabatic fluctuations, representing an improvement of over two orders of magnitude upon current constraints, and an order of magnitude over what will be possible using future CMB or galaxy surveys alone.

This paper proposes a cloud architecture for the correlation of wide bandwidth VLBI data. Cloud correlation facilitates processing of entire experiments in parallel using flexibly allocated and practically unlimited compute resources. This approach offers a potential improvement over dedicated correlation clusters, which are constrained by a fixed number of installed processor nodes and playback units. Additionally, cloud storage offers an alternative to maintaining a fleet of hard-disk drives that might be utilized intermittently. We describe benchmarks of VLBI correlation using the DiFX-2.5.2 software on the Google Cloud Platform to assess cloud-based correlation performance. The number of virtual CPUs per Virtual Machine was varied to determine the optimum configuration of cloud resources. The number of stations was varied to determine the scaling of correlation time with VLBI arrays of different sizes. Data transfer rates from Google Cloud Storage to the Virtual Machines performing the correlation were also measured. We also present an example cloud correlation configuration. Current cloud service and equipment pricing data is used to compile cost estimates allowing an approximate economic comparison between cloud and cluster processing. The economic comparisons are based on cost figures which are a moving target, and are highly dependent on factors such as the utilization of cluster and media, which are a challenge to estimate. Our model suggests that shifting to the cloud is an alternative path for high data rate, low duty cycle wideband VLBI correlation that should continue to be explored. In the production phase of VLBI correlation, the cloud has the potential to significantly reduce data processing times and allow the processing of more science experiments in a given year for the petabyte-scale data sets increasingly common in both astronomy and geodesy VLBI applications.

Observational cosmology has become an important laboratory for testing General Relativity, with searches for modified gravity forming a significant portion of the science case for existing and future surveys. In this paper, we illustrate how future measurements of the Cosmic Microwave Background (CMB) temperature and polarization anisotropies can be combined with large galaxy surveys to improve constraints on modified gravity using the technique of Sunyaev Zel'dovich (SZ) tomography. SZ tomography uses the correlations between the kinetic/polarized SZ contributions to the small-angular scale CMB and the distribution of structure measured in a galaxy redshift survey to reconstruct the remote dipole and quadrupole fields, e.g. the CMB dipole and quadrupole observed throughout the Universe. We compute the effect of a class of modifications of gravity on the remote dipole and quadrupole fields, illustrating that these observables combine a number of the desirable features of existing probes. We then perform a fisher forecast of constraints on a two-parameter class of modifications of gravity for next-generation CMB experiments and galaxy surveys. By incorporating information from the reconstructed remote dipole and quadrupole fields, we find that it is possible to improve the constraints on this model by a factor of $\sim 2$ beyond what is possible with a galaxy survey alone. We conclude that SZ tomography is a promising method for testing gravity with future cosmological datasets.

The imaging fidelity of the Event Horizon Telescope (EHT) is currently determined by its sparse baseline coverage. In particular, EHT coverage is dominated by long baselines, and is highly sensitive to atmospheric conditions and loss of sites between experiments. The limited short/mid-range baselines especially affect the imaging process, hindering the recovery of more extended features in the image. We present an algorithmic contingency for the absence of well-constrained short baselines in the imaging of compact sources, such as the supermassive black holes observed with the EHT. This technique enforces a specific second moment on the reconstructed image in the form of a size constraint, which corresponds to the curvature of the measured visibility function at zero baseline. The method enables the recovery of information lost in gaps of the baseline coverage on short baselines and enables corrections of any systematic amplitude offsets for the stations giving short-baseline measurements present in the observation. The regularization can use historical source size measurements to constrain the second moment of the reconstructed image to match the observed size. We additionally show that a characteristic size can be derived from available short-baseline measurements, extrapolated from other wavelengths, or estimated without complementary size constraints with parameter searches. We demonstrate the capabilities of this method for both static and movie reconstructions of variable sources.

This white paper envisions a revolutionary post-DESI, post-LSST dark energy program based on intensity mapping of the redshifted 21cm emission line from neutral hydrogen at radio frequencies. The proposed intensity mapping survey has the unique capability to quadruple the volume of the Universe surveyed by optical programs, provide a percent-level measurement of the expansion history to $z \sim 6$, open a window to explore physics beyond the concordance $\Lambda$CDM model, and to significantly improve the precision on standard cosmological parameters. In addition, characterization of dark energy and new physics will be powerfully enhanced by cross-correlations with optical surveys and cosmic microwave background measurements. The rich dataset obtained by the proposed intensity mapping instrument will be simultaneously useful in exploring the time-domain physics of fast radio transients and pulsars, potentially in live "multi-messenger" coincidence with other observatories. The core dark energy/inflation science advances enabled by this program are the following: (i) Measure the expansion history of the universe over $z=0.3-6$ with a single instrument, extending the range deep into the pre-acceleration era, providing an unexplored window for new physics; (ii) Measure the growth rate of structure in the universe over the same redshift range; (iii) Observe, or constrain, the presence of inflationary relics in the primordial power spectrum, improving existing constraints by an order of magnitude; (iv) Observe, or constrain, primordial non-Gaussianity with unprecedented precision, improving constraints on several key numbers by an order of magnitude. Detailed mapping of the enormous, and still largely unexplored, volume of cosmic space will thus provide unprecedented information on fundamental questions of the vacuum energy and early-universe physics.

Wide field surveys for transiting planets are well suited to searching diverse stellar populations, enabling a better understanding of the link between the properties of planets and their parent stars. We report the discovery of HAT-P-69b (TOI 625.01) and HAT-P-70b (TOI 624.01), two new hot Jupiters around A stars from the HATNet survey which have also been observed by the Transiting Exoplanet Survey Satellite (TESS). HAT-P-69b has a mass of 3.58 +0.58/-0.58 MJup and a radius of 1.676 +0.051/-0.033 RJup, residing in a prograde 4.79-day orbit. HAT-P-70b has a radius of 1.87 +0.15/-0.10 RJup and a mass constraint of < 6.78 (3 sigma) MJup, and resides in a retrograde 2.74-day orbit. We use the confirmation of these planets around relatively massive stars as an opportunity to explore the occurrence rate of hot Jupiters as a function of stellar mass. We define a sample of 47,126 main-sequence stars brighter than Tmag=10 that yields 31 giant planet candidates, including 18 confirmed planets, 3 candidates, and 10 false positives. We find a net hot Jupiter occurrence rate of 0.41+/-0.10 % within this sample, consistent with the rate measured by Kepler for FGK stars. When divided into stellar mass bins, we find the occurrence rate to be 0.71+/-0.31% for G stars, 0.43+/-0.15% for F stars, and 0.26+/-0.11% for A stars. Thus, at this point, we cannot discern any statistically significant trend in the occurrence of hot Jupiters with stellar mass.

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

The kinetic Sunyaev Zel'dovich (kSZ) effect, cosmic microwave background (CMB) temperature anisotropies induced by scattering of CMB photons from free electrons, forms the dominant blackbody component of the CMB on small angular scales. Future CMB experiments on the resolution/sensitivity frontier will measure the kSZ effect with high significance. Through the cross-correlation with tracers of structure, it will be possible to reconstruct the remote CMB dipole field (e.g. the CMB dipole observed at different locations in the Universe) using the technique of kSZ tomography. In this paper, we derive a quadratic estimator for the remote dipole field using the Cosmic Infrared Background (CIB) as a tracer of structure. We forecast the ability of current and next-generation measurements of the CMB and CIB to perform kSZ tomography. The total signal-to-noise of the reconstruction is order unity for current datasets, and can improve by a factor of up to $10^3$ for future datasets, based on statistical error only. The CIB-based reconstruction is highly correlated with a galaxy survey-based reconstruction of the remote dipole field, which can be exploited to improve constraints on cosmological parameters and models for the CIB and distribution of baryons.

The K2 mission has enabled searches for transits in crowded stellar environments very different from the original Kepler mission field. We describe here the reduction and analysis of time series data from K2's Campaign 0 superstamp, which contains the 150 Myr open cluster M35. We report on the identification of a substellar transiting object orbiting an A star at the periphery of the superstamp. To investigate this transiting source, we performed ground based follow-up observations, including photometry with the Las Cumbres Observatory telescope network and high resolution spectroscopy with Keck/HIRES. We confirm that the host star is a hot, rapidly rotating star, precluding precision radial velocity measurements. We nevertheless present a statistical validation of the planet or brown dwarf candidate using speckle interferometry from the WIYN telescope to rule out false positive stellar eclipsing binary scenarios. Based on parallax and proper motion data from Gaia Data Release 2 (DR2), we conclude that the star is not likely to be a member of M35, but instead is a background star around 100 pc behind the cluster. We present an updated ephemeris to enable future transit observations. We note that this is a rare system as a hot host star with a substellar companion. It has a high potential for future follow-up, including Doppler tomography and mid-infrared secondary transit observations.

Measurements of the Cosmic Microwave Background (CMB) temperature anisotropies on large angular scales have uncovered a number of anomalous features of marginal statistical significance, such as a hemispherical power asymmetry, lack of power on large angular scales, and features in the power spectrum. Because the primary CMB temperature has been measured at the cosmic variance limit, determining if these anomalies are hints of new physics as opposed to foregrounds, systematics, or simply statistical flukes, requires new observables. In this paper, we highlight the potential contribution that future measurements of the kinetic Sunyaev-Zel'dovich effect (kSZ) and the polarized Sunyaev Zel'dovich effect (pSZ) could make in determining the physical nature of several CMB anomalies. The kSZ and pSZ effects, temperature and polarization anisotropies induced by scattering from free electrons in the reionized Universe, are the dominant blackbody contribution to the CMB on small angular scales. Using the technique of SZ tomography, measurements of kSZ and pSZ effects can be combined with galaxy surveys to reconstruct the remote CMB dipole and quadrupole fields, providing a 3-dimensional probe of large scale modes inside our Hubble volume. Building on previous work, we forecast the additional constraining power that these observables might offer for a representative set of anomaly models. We find that the remote CMB dipole and quadrupole contain a similar amount of information on anomaly models as the primary CMB polarization. The information from CMB temperature, polarization, and the remote dipole and quadrupole fields is complementary, and the full set of observables can improve constraints on anomaly models by a factor of $\sim 2-4$ using next-generation CMB experiments and galaxy surveys. This could be sufficient to definitively establish the physical origin of several CMB anomalies.

A millimeter-wave survey over half the sky, that spans frequencies in the range of 30 to 350 GHz, and that is both an order of magnitude deeper and of higher-resolution than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. By providing such a deep, high-resolution millimeter-wave survey (about 0.5 uK-arcmin noise and 15 arcsecond resolution at 150 GHz), CMB-HD will enable major advances. It will allow 1) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10/hMpc), which probes dark matter particle properties. It will also allow 2) measurements of the thermal and kinetic Sunyaev-Zel'dovich effects on small scales to map the gas density and gas pressure profiles of halos over a wide field, which probes galaxy evolution and cluster astrophysics. In addition, CMB-HD would allow us to cross critical thresholds in fundamental physics: 3) ruling out or detecting any new, light (< 0.1eV), thermal particles, which could potentially be the dark matter, and 4) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe. Such a survey would also 5) monitor the transient sky by mapping the full observing region every few days, which opens a new window on gamma-ray bursts, novae, fast radio bursts, and variable active galactic nuclei. Moreover, CMB-HD would 6) provide a census of planets, dwarf planets, and asteroids in the outer Solar System, and 7) enable the detection of exo-Oort clouds around other solar systems, shedding light on planet formation. CMB-HD will deliver this survey in 5 years of observing half the sky, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. The telescopes will field about 2.4 million detectors (600,000 pixels) in total.

Young exoplanets are snapshots of the planetary evolution process. Planets that orbit stars in young associations are particularly important because the age of the planetary system is well constrained. We present the discovery of a transiting planet larger than Neptune but smaller than Saturn in the 45 Myr Tucana-Horologium young moving group. The host star is a visual binary, and our follow-up observations demonstrate that the planet orbits the G6V primary component, DS Tuc A (HD 222259A, TIC 410214986). We first identified transits using photometry from the Transiting Exoplanet Survey Satellite (TESS; alerted as TOI 200.01). We validated the planet and improved the stellar parameters using a suite of new and archival data, including spectra from SOAR/Goodman, SALT/HRS and LCO/NRES; transit photometry from Spitzer; and deep adaptive optics imaging from Gemini/GPI. No additional stellar or planetary signals are seen in the data. We measured the planetary parameters by simultaneously modeling the photometry with a transit model and a Gaussian process to account for stellar variability. We determined that the planetary radius is $5.70\pm0.17$ Earth radii and that the orbital period is 8.1 days. The inclination angles of the host star's spin axis, the planet's orbital axis, and the visual binary's orbital axis are aligned within 15 degrees to within the uncertainties of the relevant data. DS Tuc Ab is bright enough (V=8.5) for detailed characterization using radial velocities and transmission spectroscopy.

Very-long-baseline interferometry (VLBI) at frequencies above 230 GHz with Earth-diameter baselines gives spatial resolution finer than the ${\sim}50 \mu$as "shadow" of the supermassive black hole at the Galactic Center, Sagittarius A* (Sgr A*). Imaging static and dynamical structure near the "shadow" provides a test of general relativity and may allow measurement of black hole parameters. However, traditional Earth-rotation synthesis is inapplicable for sources (such as Sgr A*) with intra-day variability. Expansions of ground-based arrays to include space-VLBI stations may enable imaging capability on time scales comparable to the prograde innermost stable circular orbit (ISCO) of Sgr A*, which is predicted to be 4-30 minutes, depending on black hole spin. We examine the basic requirements for space-VLBI, and we develop tools for simulating observations with orbiting stations. We also develop a metric to quantify the imaging capabilities of an array irrespective of detailed image morphology or reconstruction method. We validate this metric on example reconstructions of simulations of Sgr A* at 230 and 345 GHz, and use these results to motivate expanding the Event Horizon Telescope (EHT) to include small dishes in Low Earth Orbit (LEO). We demonstrate that high-sensitivity sites such as the Atacama Large Millimeter/Submillimeter Array (ALMA) make it viable to add small orbiters to existing ground arrays, as space-ALMA baselines would have sensitivity comparable to ground-based non-ALMA baselines. We show that LEO-enhanced arrays sample half of the diffraction-limited Fourier plane of Sgr A* in less than 30 minutes, enabling reconstructions of near-horizon structure with normalized root-mean-square error $\lesssim0.3$ on sub-ISCO timescales.

We introduce a new picture of vacuum decay which, in contrast to existing semiclassical techniques, provides a real-time description and does not rely on classically-forbidden tunneling paths. Using lattice simulations, we observe vacuum decay via bubble formation by generating realizations of vacuum fluctuations and evolving with the classical equations of motion. The decay rate obtained from an ensemble of simulations is in excellent agreement with existing techniques. Future applications include bubble correlation functions, fast decay rates, and decay of non-vacuum states.

We present HI 21 cm spectroscopy from the GBT for the host galaxies of 31 nearby AGNs with direct M$_{\textrm{BH}}$ measurements from reverberation mapping. These are the first published HI detections for 12 galaxies, and the spectral quality is generally an improvement over archival data for the remainder of the sample. We present measurements of emission-line fluxes, velocity widths, and recessional velocities from which we derive HI mass, total gas mass, and redshifts. Combining M$_{\textrm{GAS}}$ with constraints on M$_{\textrm{STARS}}$ allows exploration of the baryonic content of these galaxies. We find a typical M$_{\textrm{GAS}}$/M$_{\textrm{STARS}}$ fraction of 10%, with a few reaching $\sim$30-50%. We also examined several relationships between M$_{\textrm{STARS}}$, M$_{\textrm{GAS}}$, M$_{\textrm{BH}}$, baryonic mass, and morphological type. We find a weak preference for galaxies with larger M$_{\textrm{GAS}}$ to host more massive black holes. We also find gas-to-stellar fractions to weakly correlate with later types in unbarred spirals, with an approximately constant fraction for barred spirals. Consistent with previous studies, we find declining M$_{\textrm{GAS}}$/M$_{\textrm{STARS}}$ with increasing M$_{\textrm{STARS}}$, with a slope suggesting the gas reservoirs have been replenished. Finally, we find a clear relationship for M$_{\textrm{BH}}$-M$_{\textrm{BARY}}$ with a similar slope as M$_{\textrm{BH}}$-M$_{\textrm{STARS}}$ reported by Bentz & Manne-Nicholas (2018). The dwarf Seyfert NGC 4395 appears to follow this relationship as well, even though it has a significantly higher gas fraction and smaller M$_{\textrm{BH}}$ than the remainder of our sample.

Our aim is to investigate the nature and formation of brown dwarfs by adding a new well-characterised object to the small sample of less than 20 transiting brown dwarfs. One brown dwarf candidate was found by the KESPRINT consortium when searching for exoplanets in the K2 space mission Campaign 16 field. We combined the K2 photometric data with a series of multi-colour photometric observations, imaging and radial velocity measurements to rule out false positive scenarios and to determine the fundamental properties of the system. We report the discovery and characterisation of a transiting brown dwarf in a 5.17 day eccentric orbit around the slightly evolved F7V star EPIC 212036875. We find a stellar mass of 1.15+/-0.08 M$_\odot$, a stellar radius of 1.41+/-0.05 R$_\odot$, and an age of 5.1+/-0.9 Gyr. The mass and radius of the companion brown dwarf are 51+/-2 MJ and 0.83+/-0.03 RJ, respectively, corresponding to a mean density of 108+15-13 g cm-3. EPIC 212036875 b is a rare object that resides in the brown dwarf desert. In the mass-density diagram for planets, brown dwarfs and stars, we find that all giant planets and brown dwarfs follow the same trend from ~0.3 MJ to the turn-over to hydrogen burning stars at ~73 MJ. EPIC 212036875 b falls close to the theoretical model for mature H/He dominated objects in this diagram as determined by interior structure models, as well as the empirical fit. We argue that EPIC 212036875 b formed via gravitational disc instabilities in the outer part of the disc, followed by a quick migration. Orbital tidal circularisation may have started early in its history for a brief period when the brown dwarf's radius was larger. The lack of spin-orbit synchronisation points to a weak stellar dissipation parameter which implies a circularisation timescale of >23 Gyr, or suggests an interaction between the magnetic and tidal forces of the star and the brown dwarf.

We announce the discovery of KELT-23Ab, a hot Jupiter transiting the relatively bright ($V=10.3$) star BD+66 911 (TYC 4187-996-1), and characterize the system using follow-up photometry and spectroscopy. A global fit to the system yields host-star properties of $T_{eff}=5900\pm49 K$, $M_*=0.945^{+0.060}_{-0.054} M_{\odot}$, $R_*=0.995\pm0.015 R_{\odot}$, $L_*=1.082^{+0.051}_{-0.048} L_{\odot}$, log$g_{*}=4.418^{+0.026}_{-0.025}$ (cgs), and $\left[{\rm Fe}/{\rm H}\right]=-0.105\pm0.077$. KELT-23Ab is a hot Jupiter with mass $M_P=0.938^{+0.045}_{-0.042} M_{\rm J}$, radius $R_P=1.322\pm0.025 R_{\rm J}$, and density $\rho_P=0.504^{+0.038}_{-0.035}$ g cm$^{-3}$. Intense insolation flux from the star has likely caused KELT-23Ab to become inflated. The time of inferior conjunction is $T_0=2458149.40776\pm0.00091~\rm {BJD_{TDB}}$ and the orbital period is $P=2.255353^{+0.000031}_{-0.000030}$ days. There is strong evidence that KELT-23A is a member of a long-period binary star system with a less luminous companion, and due to tidal interactions, the planet is likely to spiral into its host within roughly a Gyr. This system has one of the highest positive ecliptic latitudes of all transiting planet hosts known to date, placing it near the Transiting Planet Survey Satellite and James Webb Space Telescope continuous viewing zones. Thus we expect it to be an excellent candidate for long-term monitoring and follow-up with these facilities.

The CHANG-ES galaxy sample consists of 35 nearby edge-on galaxies that have been observed using the VLA at 1.6 GHz and 6.0 GHz. Here we present the 3rd data release of our sample, namely the B-configuration 1.6 GHz sample. In addition, we make available the {\it band-to-band} spectral index maps between 1.6 GHz and 6.0 GHz, the latter taken in the matching resolution C-configuration. The images can be downloaded from https://www.queensu.ca/changes. These are our highest resolution images ($\approx$ 3 arcsec) and we examine the possible presence of low luminosity active galactic nuclei in the sample as well as some in-disk structure. New features can be seen in the spectral index maps that are masked in the total intensity emission, including hidden spiral arms in NGC~3448 and two previously unknown radio lobes on either side of the nucleus of NGC~3628. Our AGN detection rate, using only radio criteria, is 55\% which we take as a lower limit because some weaker embedded AGNs are likely present which could be revealed at higher resolution. Archival XMM-Newton data were used to search for further fingerprints of the AGNs in the studied sample. In galaxy disks, discrete regions of flat spectral index are seen, likely due to a thermal emission fraction that is higher than the global average.

We present the first cosmological parameter constraints using measurements of type Ia supernovae (SNe Ia) from the Dark Energy Survey Supernova Program (DES-SN). The analysis uses a subsample of 207 spectroscopically confirmed SNe Ia from the first three years of DES-SN, combined with a low-redshift sample of 122 SNe from the literature. Our "DES-SN3YR" result from these 329 SNe Ia is based on a series of companion analyses and improvements covering SN Ia discovery, spectroscopic selection, photometry, calibration, distance bias corrections, and evaluation of systematic uncertainties. For a flat LCDM model we find a matter density Omega_m = 0.331 +_ 0.038. For a flat wCDM model, and combining our SN Ia constraints with those from the cosmic microwave background (CMB), we find a dark energy equation of state w = -0.978 +_ 0.059, and Omega_m = 0.321 +_ 0.018. For a flat w0waCDM model, and combining probes from SN Ia, CMB and baryon acoustic oscillations, we find w0 = -0.885 +_ 0.114 and wa = -0.387 +_ 0.430. These results are in agreement with a cosmological constant and with previous constraints using SNe Ia (Pantheon, JLA).

We present the first cosmological parameter constraints using measurements of type Ia supernovae (SNe Ia) from the Dark Energy Survey Supernova Program (DES-SN). The analysis uses a subsample of 207 spectroscopically confirmed SNe Ia from the first three years of DES-SN, combined with a low-redshift sample of 122 SNe from the literature. Our "DES-SN3YR" result from these 329 SNe Ia is based on a series of companion analyses and improvements covering SN Ia discovery, spectroscopic selection, photometry, calibration, distance bias corrections, and evaluation of systematic uncertainties. For a flat LCDM model we find a matter density Omega_m = 0.331 +_ 0.038. For a flat wCDM model, and combining our SN Ia constraints with those from the cosmic microwave background (CMB), we find a dark energy equation of state w = -0.978 +_ 0.059, and Omega_m = 0.321 +_ 0.018. For a flat w0waCDM model, and combining probes from SN Ia, CMB and baryon acoustic oscillations, we find w0 = -0.885 +_ 0.114 and wa = -0.387 +_ 0.430. These results are in agreement with a cosmological constant and with previous constraints using SNe Ia (Pantheon, JLA).

We describe catalog-level simulations of Type Ia supernova (SN~Ia) light curves in the Dark Energy Survey Supernova Program (DES-SN), and in low-redshift samples from the Center for Astrophysics (CfA) and the Carnegie Supernova Project (CSP). These simulations are used to model biases from selection effects and light curve analysis, and to determine bias corrections for SN~Ia distance moduli that are used to measure cosmological parameters. To generate realistic light curves, the simulation uses a detailed SN~Ia model, incorporates information from observations (PSF, sky noise, zero point), and uses summary information (e.g., detection efficiency vs. signal to noise ratio) based on 10,000 fake SN light curves whose fluxes were overlaid on images and processed with our analysis pipelines. The quality of the simulation is illustrated by predicting distributions observed in the data. Averaging within redshift bins, we find distance modulus biases up to 0.05~mag over the redshift ranges of the low-z and DES-SN samples. For individual events, particularly those with extreme red or blue color, distance biases can reach 0.4~mag. Therefore, accurately determining bias corrections is critical for precision measurements of cosmological parameters. Files used to make these corrections are available at https://des.ncsa.illinois.edu/releases/sn.

We describe catalog-level simulations of Type Ia supernova (SN~Ia) light curves in the Dark Energy Survey Supernova Program (DES-SN), and in low-redshift samples from the Center for Astrophysics (CfA) and the Carnegie Supernova Project (CSP). These simulations are used to model biases from selection effects and light curve analysis, and to determine bias corrections for SN~Ia distance moduli that are used to measure cosmological parameters. To generate realistic light curves, the simulation uses a detailed SN~Ia model, incorporates information from observations (PSF, sky noise, zero point), and uses summary information (e.g., detection efficiency vs. signal to noise ratio) based on 10,000 fake SN light curves whose fluxes were overlaid on images and processed with our analysis pipelines. The quality of the simulation is illustrated by predicting distributions observed in the data. Averaging within redshift bins, we find distance modulus biases up to 0.05~mag over the redshift ranges of the low-z and DES-SN samples. For individual events, particularly those with extreme red or blue color, distance biases can reach 0.4~mag. Therefore, accurately determining bias corrections is critical for precision measurements of cosmological parameters. Files used to make these corrections are available at https://des.ncsa.illinois.edu/releases/sn.

The astrophysical $S$-factor for the radiative proton capture reaction on $^7$Be ($S_{17}$) at low energies is affected by the $s$-wave scattering lengths. We report the measurement of elastic and inelastic scattering cross sections for the $^7$Be+p system in the center-of-mass energy range 0.474 - 2.740 MeV and center-of-mass angular range of 70$^\circ$- 150$^\circ$. A radioactive $^7$Be beam produced at Oak Ridge National Laboratory's (ORNL) Holifield Radioactive Ion Beam Facility was accelerated and bombarded a thin polypropylene (CH$_{2}$)$_\text n$ target. Scattered ions were detected in the segmented Silicon Detector Array. Using an $\textit{R}$-matrix analysis of ORNL and Louvain-la-Neuve cross section data, the $s$-wave scattering lengths for channel spins 1 and 2 were determined to be 17.34$^{+1.11}_{-1.33}$ and -3.18$^{+0.55}_{-0.50}$ fm, respectively. The uncertainty in the $s$-wave scattering lengths reported in this work is smaller by a factor of 5-8 compared to the previous measurement, which may reduce the overall uncertainty in $S_{17}$ at zero energy. The level structure of $^8$B is discussed based upon the results from this work. Evidence for the existence of 0$^+$ and 2$^+$ levels in $^8$B at 1.9 and 2.21 MeV, respectively, is observed.

We describe the first public data release of the Dark Energy Survey, DES DR1, consisting of reduced single epoch images, coadded images, coadded source catalogs, and associated products and services assembled over the first three years of DES science operations. DES DR1 is based on optical/near-infrared imaging from 345 distinct nights (August 2013 to February 2016) by the Dark Energy Camera mounted on the 4-m Blanco telescope at Cerro Tololo Inter-American Observatory in Chile. We release data from the DES wide-area survey covering ~5,000 sq. deg. of the southern Galactic cap in five broad photometric bands, grizY. DES DR1 has a median delivered point-spread function of g = 1.12, r = 0.96, i = 0.88, z = 0.84, and Y = 0.90 arcsec FWHM, a photometric precision of < 1% in all bands, and an astrometric precision of 151 mas. The median coadded catalog depth for a 1.95" diameter aperture at S/N = 10 is g = 24.33, r = 24.08, i = 23.44, z = 22.69, and Y = 21.44 mag. DES DR1 includes nearly 400M distinct astronomical objects detected in ~10,000 coadd tiles of size 0.534 sq. deg. produced from ~39,000 individual exposures. Benchmark galaxy and stellar samples contain ~310M and ~ 80M objects, respectively, following a basic object quality selection. These data are accessible through a range of interfaces, including query web clients, image cutout servers, jupyter notebooks, and an interactive coadd image visualization tool. DES DR1 constitutes the largest photometric data set to date at the achieved depth and photometric precision.

We describe the first public data release of the Dark Energy Survey, DES DR1, consisting of reduced single epoch images, coadded images, coadded source catalogs, and associated products and services assembled over the first three years of DES science operations. DES DR1 is based on optical/near-infrared imaging from 345 distinct nights (August 2013 to February 2016) by the Dark Energy Camera mounted on the 4-m Blanco telescope at Cerro Tololo Inter-American Observatory in Chile. We release data from the DES wide-area survey covering ~5,000 sq. deg. of the southern Galactic cap in five broad photometric bands, grizY. DES DR1 has a median delivered point-spread function of g = 1.12, r = 0.96, i = 0.88, z = 0.84, and Y = 0.90 arcsec FWHM, a photometric precision of < 1% in all bands, and an astrometric precision of 151 mas. The median coadded catalog depth for a 1.95" diameter aperture at S/N = 10 is g = 24.33, r = 24.08, i = 23.44, z = 22.69, and Y = 21.44 mag. DES DR1 includes nearly 400M distinct astronomical objects detected in ~10,000 coadd tiles of size 0.534 sq. deg. produced from ~39,000 individual exposures. Benchmark galaxy and stellar samples contain ~310M and ~ 80M objects, respectively, following a basic object quality selection. These data are accessible through a range of interfaces, including query web clients, image cutout servers, jupyter notebooks, and an interactive coadd image visualization tool. DES DR1 constitutes the largest photometric data set to date at the achieved depth and photometric precision.

High-resolution, low-noise observations of the cosmic microwave background (CMB) planned for the near-future will enable new cosmological probes based on re-scattered CMB photons -- the secondary CMB. At the same time, enormous galaxy surveys will map out huge volumes of the observable Universe. Using the technique of kinetic Sunyaev Zel'dovich (kSZ) tomography these new probes can be combined to reconstruct the remote dipole field, the CMB dipole as observed from different vantage points in our Universe. The volume accessible to future galaxy surveys is large enough that general relativistic corrections to the observed distribution of galaxies must be taken into account. These corrections are interesting probes of gravity in their own right, but can also obscure potential signatures of primordial non-Gaussianity. In this paper, we demonstrate that correlations between the reconstructed remote dipole field and the observed galaxy number counts can in principle be used to detect general relativistic corrections. We show that neglecting general relativistic corrections leads to an $\mathcal{O}(1)$ bias on the inferred amplitude of primordial non-Gaussianity, $f_{\rm NL}$. In addition, we demonstrate that the reconstructed remote dipole field can provide useful constraining power on various bias parameters appearing in the galaxy number counts, and can significantly mitigate the effects of alignment bias.

We present two-temperature, radiative general relativistic magnetohydrodynamic simulations of Magnetically Arrested Discs (MAD) that launch powerful relativistic jets. The mass accretion rates of our simulations are scaled to match the luminosity of the accretion flow around the supermassive black hole in M87. We consider two sub-grid prescriptions for electron heating: one based on a Landau-damped turbulent cascade, and the other based on heating from trans-relativistic magnetic reconnection. The simulations produce jets with power on the order of the observed value for M87. Both simulations produce spectra that are consistent with observations of M87 in the radio, millimetre, and submillimetre. Furthermore, the predicted image core-shifts in both models at frequencies between 15 GHz and 86 GHz are consistent with observations. At 43 and 86 GHz, both simulations produce wide opening angle jets consistent with VLBI images. Both models produce 230~GHz images with distinct black hole shadows that are resolvable by the Event Horizon Telescope (EHT), although at a viewing angle of 17 degrees, the 230 GHz images are too large to match EHT observations from 2009 and 2012. The 230 GHz images from the simulations are dynamic on time-scales of months to years, suggesting that repeated EHT observations may be able to detect the motion of rotating magnetic fields at the event horizon.

Astrometry, the measurement of positions and motions of the stars, is one of the oldest disciplines in Astronomy, extending back at least as far as Hipparchus' discovery of the precession of Earth's axes in 190 BCE by comparing his catalog with those of his predecessors. Astrometry is fundamental to Astronomy, and critical to many aspects of Astrophysics and Geodesy. In order to understand our planet's and solar system's context within their surroundings, we must be able to to define, quantify, study, refine, and maintain an inertial frame of reference relative to which all positions and motions can be unambiguously and self-consistently described. It is only by using this inertial reference frame that we are able to disentangle our observations of the motions of celestial objects from our own complex path around our star, and its path through the galaxy, and the local group. Every aspect of each area outlined in the call for scientific frontiers in astronomy in the era of the 2020-2030 timeframe will depend on the quality of the inertial reference frame. In this white paper, we propose support for development of radio Very Long Baseline Interferometry (VLBI) capabilities, including the Next Generation Very Large Array (ngVLA), a radio astronomy observatory that will not only support development of a next generation reference frame of unprecedented accuracy, but that will also serve as a highly capable astronomical instrument in its own right. Much like its predecessors, the Very Long Baseline Array (VLBA) and other VLBI telescopes, the proposed ngVLA will provide the foundation for the next three decades for the fundamental reference frame, benefitting astronomy, astrophysics, and geodesy alike.

Having discovered that Earth-sized planets are common, we are now embarking on a journey to determine if Earth-like planets are also common. Finding Earth-like planets is one of the most compelling endeavors of the 21st century - leading us toward finally answering the question: Are we alone? To achieve this forward-looking goal, we must determine the masses of the planets; the sizes of the planets, by themselves, are not sufficient for the determination of the bulk and atmospheric compositions. Masses, coupled with the radii, are crucial constraints on the bulk composition and interior structure of the planets and the composition of their atmospheres, including the search for biosignatures. Precision radial velocity is the most viable technique for providing essential mass and orbit information for spectroscopy of other Earths. The development of high quality precision radial velocity instruments coupled to the building of the large telescope facilities like TMT and GMT or space-based platforms like EarthFinder can enable very high spectral resolution observations with extremely precise radial velocities on minute timescales to allow for the modeling and removal of radial velocity jitter. Over the next decade, the legacy of exoplanet astrophysics can be cemented firmly as part of humankind's quest in finding the next Earth - but only if we can measure the masses and orbits of Earth-sized planets in habitable zone orbits around Sun-like stars.

We recommend an intensive effort to survey and understand the obliquity distribution of small close-in extrasolar planets over the coming decade. The orbital obliquities of exoplanets--i.e., the relative orientation between the planetary orbit and the stellar rotation--is a key tracer of how planets form and migrate. While the orbital obliquities of smaller planets are poorly explored today, a new generation of facilities coming online over the next decade will make such observations possible en masse. Transit spectroscopic observations with the extremely large telescopes will enable us to measure the orbital obliquities of planets as small as $\sim2R_{\oplus}$ around a wide variety of stars, opening a window into the orbital properties of the most common types of planets. This effort will directly contribute to understanding the formation and evolution of planetary systems, a key objective of the National Academy of Sciences' Exoplanet Science Strategies report.

The identification of the electromagnetic counterpart of the gravitational wave event, GW170817, and discovery of neutrinos and gamma-rays from TXS 0506+056 heralded the new era of multi-messenger astrophysics. As the number of multi-messenger events rapidly grow over the next decade, the cyberinfrastructure requirements to handle the increase in data rates, data volume, need for event follow up, and analysis across the different messengers will also explosively grow. The cyberinfrastructure requirements to enhance multi-messenger astrophysics will both be a major challenge and opportunity for astronomers, physicists, computer scientists and cyberinfrastructure specialists. Here we outline some of these requirements and argue for a distributed cyberinfrastructure institute for multi-messenger astrophysics to meet these challenges.

Opening up a new window of millimeter-wave observations that span frequency bands in the range of 30 to 500 GHz, survey half the sky, and are both an order of magnitude deeper (about 0.5 uK-arcmin) and of higher-resolution (about 10 arcseconds) than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. In particular, such a survey would allow for major advances in measuring the distribution of dark matter and gas on small-scales, and yield needed insight on 1.) dark matter particle properties, 2.) the evolution of gas and galaxies, 3.) new light particle species, 4.) the epoch of inflation, and 5.) the census of bodies orbiting in the outer Solar System.

We present cosmological results from a combined analysis of galaxy clustering and weak gravitational lensing, using 1321 deg$^2$ of $griz$ imaging data from the first year of the Dark Energy Survey (DES Y1). We combine three two-point functions: (i) the cosmic shear correlation function of 26 million source galaxies in four redshift bins, (ii) the galaxy angular autocorrelation function of 650,000 luminous red galaxies in five redshift bins, and (iii) the galaxy-shear cross-correlation of luminous red galaxy positions and source galaxy shears. To demonstrate the robustness of these results, we use independent pairs of galaxy shape, photometric redshift estimation and validation, and likelihood analysis pipelines. To prevent confirmation bias, the bulk of the analysis was carried out while blind to the true results; we describe an extensive suite of systematics checks performed and passed during this blinded phase. The data are modeled in flat $\Lambda$CDM and $w$CDM cosmologies, marginalizing over 20 nuisance parameters, varying 6 (for $\Lambda$CDM) or 7 (for $w$CDM) cosmological parameters including the neutrino mass density and including the 457 $\times$ 457 element analytic covariance matrix. We find consistent cosmological results from these three two-point functions, and from their combination obtain $S_8 \equiv \sigma_8 (\Omega_m/0.3)^{0.5} = 0.783^{+0.021}_{-0.025}$ and $\Omega_m = 0.264^{+0.032}_{-0.019}$ for $\Lambda$CDM for $w$CDM, we find $S_8 = 0.794^{+0.029}_{-0.027}$, $\Omega_m = 0.279^{+0.043}_{-0.022}$, and $w=-0.80^{+0.20}_{-0.22}$ at 68% CL. The precision of these DES Y1 results rivals that from the Planck cosmic microwave background measurements, allowing a comparison of structure in the very early and late Universe on equal terms. Although the DES Y1 best-fit values for $S_8$ and $\Omega_m$ are lower than the central values from Planck ...

We present cosmological results from a combined analysis of galaxy clustering and weak gravitational lensing, using 1321 deg$^2$ of $griz$ imaging data from the first year of the Dark Energy Survey (DES Y1). We combine three two-point functions: (i) the cosmic shear correlation function of 26 million source galaxies in four redshift bins, (ii) the galaxy angular autocorrelation function of 650,000 luminous red galaxies in five redshift bins, and (iii) the galaxy-shear cross-correlation of luminous red galaxy positions and source galaxy shears. To demonstrate the robustness of these results, we use independent pairs of galaxy shape, photometric redshift estimation and validation, and likelihood analysis pipelines. To prevent confirmation bias, the bulk of the analysis was carried out while blind to the true results; we describe an extensive suite of systematics checks performed and passed during this blinded phase. The data are modeled in flat $\Lambda$CDM and $w$CDM cosmologies, marginalizing over 20 nuisance parameters, varying 6 (for $\Lambda$CDM) or 7 (for $w$CDM) cosmological parameters including the neutrino mass density and including the 457 $\times$ 457 element analytic covariance matrix. We find consistent cosmological results from these three two-point functions, and from their combination obtain $S_8 \equiv \sigma_8 (\Omega_m/0.3)^{0.5} = 0.783^{+0.021}_{-0.025}$ and $\Omega_m = 0.264^{+0.032}_{-0.019}$ for $\Lambda$CDM for $w$CDM, we find $S_8 = 0.794^{+0.029}_{-0.027}$, $\Omega_m = 0.279^{+0.043}_{-0.022}$, and $w=-0.80^{+0.20}_{-0.22}$ at 68% CL. The precision of these DES Y1 results rivals that from the Planck cosmic microwave background measurements, allowing a comparison of structure in the very early and late Universe on equal terms. Although the DES Y1 best-fit values for $S_8$ and $\Omega_m$ are lower than the central values from Planck ...

The observation of gravitational waves from a binary neutron star merger by LIGO/VIRGO and the associated electromagnetic counterpart provides a high precision test of orbital dynamics, and therefore a new and sensitive probe of extra forces and new radiative degrees of freedom. Axions are one particularly well-motivated class of extensions to the Standard Model leading to new forces and sources of radiation, which we focus on in this paper. Using an effective field theory (EFT) approach, we calculate the first post-Newtonian corrections to the orbital dynamics, radiated power, and gravitational waveform for binary neutron star mergers in the presence of an axion. This result is applicable to many theories which add an extra massive scalar degree of freedom to General Relativity. We then perform a detailed forecast of the potential for Advanced LIGO to constrain the free parameters of the EFT, and map these to the mass $m_a$ and decay constant $f_a$ of the axion. At design sensitivity, we find that Advanced LIGO can potentially exclude axions with $m_a \lesssim 10^{-11} \ {\rm eV}$ and $f_a \sim (10^{14} - 10^{17}) \ {\rm GeV}$. There are a variety of complementary observational probes over this region of parameter space, including the orbital decay of binary pulsars, black hole superradiance, and laboratory searches. We comment on the synergies between these various observables.

We investigate how the spectral properties of atomic (HI) and molecular (H$_2$) gas, traced by CO(2-1), are related in M33 on $80$ pc scales. We find the HI and CO(2-1) velocity at peak intensity to be highly correlated, consistent with previous studies. By stacking spectra aligned to the velocity of HI peak intensity, we find that the CO line width ($\sigma_{\rm HWHM}=4.6\pm0.9$ km s$^{-1}$; $\sigma_{\rm HWHM}$ is the effective Gaussian width) is consistently smaller than the HI line width ($\sigma_{\rm HWHM}=6.6\pm0.1$ km s$^{-1}$), with a ratio of ${\sim}0.7$, in agreement with Druard et al. (2014). The ratio of the line widths remains less than unity when the data are smoothed to a coarser spatial resolution. In other nearby galaxies, this line width ratio is close to unity which has been used as evidence for a thick, diffuse molecular disk that is distinct from the thin molecular disk dominated by molecular clouds. The smaller line width ratio found here suggests that M33 has a marginal thick molecular disk. From modelling individual lines-of-sight, we recover a strong correlation between HI and CO line widths when only the HI located closest to the CO component is considered. The median line width ratio of the line-of-sight line widths is $0.56\pm0.01$. There is substantial scatter in the HI--CO(2-1) line width relation, larger than the uncertainties, that results from regional variations on $<500$ pc scales, and there is no significant trend in the line widths, or their ratios, with galactocentric radius. These regional line width variations may be a useful probe of changes in the local cloud environment or the evolutionary state of molecular clouds.

We report on a first census of Galactic black hole X-ray binary (BHXRB) properties with the second data release (DR2) of {\em Gaia}, focusing on dynamically confirmed and strong candidate black hole transients. DR2 provides five-parameter astrometric solutions including position, parallax and proper motion for 11 of a sample of 24 systems. Distance estimates are tested with parallax inversion as well as Bayesian inference. We derive an empirically motivated characteristic scale length of $L$=2.17$\pm$0.12 kpc for this BHXRB population to infer distances based upon an exponentially decreasing space density prior. Geometric DR2 parallaxes provide new, independent distance estimates, but the faintness of this population in quiescence results in relatively large fractional distance uncertainties. Despite this, DR2 estimates generally agree with literature distances. The most discrepant case is BW Cir, for which detailed studies of the donor star have suggested a distant location at >~25 kpc. A large DR2 measured parallax and relatively high proper motion instead prefer significantly smaller distances, suggesting that the source may instead be amongst the nearest of XRBs. However, both distances create problems for interpretation of the source, and follow-up data are required to resolve its true nature. DR2 also provides a first distance estimate to one source, MAXI J1820+070, and novel proper motion estimates for 7 sources. Peculiar velocities relative to Galactic rotation exceed $\sim$ 50 km s$^{-1}$ for the bulk of the sample, with a median system kinetic energy of peculiar motion of $\sim$ 5 $\times$ 10$^{47}$ erg. BW Cir could be a new high-velocity BHXRB if its astrometry is confirmed. A putative anti-correlation between peculiar velocity and black hole mass is found, as expected in mass-dependent BH kick formation channels, but this trend remains weak in the DR2 data.

The Galactic Center supermassive black hole Sagittarius A* (Sgr A*) is one of the most promising targets to study the dynamics of black hole accretion and outflow via direct imaging with very long baseline interferometry (VLBI). At 3.5 mm (86 GHz), the emission from Sgr A* is resolvable with the Global Millimeter VLBI Array (GMVA). We present the first observations of Sgr A* with the phased Atacama Large Millimeter/submillimeter Array (ALMA) joining the GMVA. Our observations achieve an angular resolution of ~87{\mu}as, improving upon previous experiments by a factor of two. We reconstruct a first image of the unscattered source structure of Sgr A* at 3.5 mm, mitigating effects of interstellar scattering. The unscattered source has a major axis size of 120 $\pm$ 34{\mu}as (12 $\pm$ 3.4 Schwarzschild radii), and a symmetrical morphology (axial ratio of 1.2$^{+0.3}_{-0.2}$), which is further supported by closure phases consistent with zero within 3{\sigma}. We show that multiple disk-dominated models of Sgr A* match our observational constraints, while the two jet-dominated models considered are constrained to small viewing angles. Our long-baseline detections to ALMA also provide new constraints on the scattering of Sgr A*, and we show that refractive scattering effects are likely to be weak for images of Sgr A* at 1.3 mm with the Event Horizon Telescope. Our results provide the most stringent constraints to date for the intrinsic morphology and refractive scattering of Sgr A*, demonstrating the exceptional contribution of ALMA to millimeter VLBI.

We present Karl G. Jansky Very Large Array radio observations of the long gamma-ray burst GRB 161219B ($z=0.147$) spanning $1-37$ GHz. The data exhibit unusual behavior, including sharp spectral peaks and minutes-timescale large-amplitude variability centered at $20$ GHz and spanning the full frequency range. We attribute this behavior to scattering of the radio emission by the turbulent ionized Galactic interstellar medium (ISM), including both diffractive and refractive scintillation. However, the scintillation is much stronger than predicted by a model of the Galactic electron density distribution (NE2001); from the measured variability timescale and decorrelation bandwidth we infer a scattering measure of $SM\approx {(8-70)\times 10^{-4}}$ kpc m$^{-20/3}$ (up to ${25}$ times larger than predicted in NE2001) and a scattering screen distance of $d_{\rm scr}\approx {0.2-3}$ kpc. We infer an emission region size of $\theta_s \approx {0.9-4}$ $\mu$as ($\approx {(1-4)}\times 10^{16}$ cm) at $\approx4$ days, and find that prior to 8 days the source size is an order of magnitude smaller than model predictions for a uniformly illuminated disk or limb-brightened ring, indicating a slightly off-axis viewing angle or significant substructure in the emission region. Simultaneous multi-hour broadband radio observations of future GRB afterglows will allow us to characterize the scintillation more completely, and hence to probe the observer viewing angle, the evolution of the jet Lorentz factor, the structure of the afterglow emission regions, and ISM turbulence at high Galactic latitudes.

The Large Synoptic Survey Telescope (LSST) will provide for unbiased sampling of variability properties of objects with $r$ mag $<$ 24. This should allow for those objects whose variations reveal their orbital periods ($P_{orb}$), such as low mass X-ray binaries (LMXBs) and related objects, to be examined in much greater detail and with uniform systematic sampling. However, the baseline LSST observing strategy has temporal sampling that is not optimised for such work in the Galaxy. Here we assess four candidate observing strategies for measurement of $P_{orb}$ in the range 10 minutes to 50 days. We simulate multi-filter quiescent LMXB lightcurves including ellipsoidal modulation and stochastic flaring, and then sample these using LSST's operations simulator (OpSim) over the (mag, $P_{orb}$) parameter space, and over five sightlines sampling a range of possible reddening values. The percentage of simulated parameter space with correctly returned periods ranges from $\sim$23 %, for the current baseline strategy, to $\sim$70 % for the two simulated specialist strategies. Convolving these results with a $P_{orb}$ distribution, a modelled Galactic spatial distribution and reddening maps, we conservatively estimate that the most recent version of the LSST baseline strategy will allow $P_{orb}$ determination for $\sim$18 % of the Milky Way's LMXB population, whereas strategies that do not reduce observations of the Galactic Plane can improve this dramatically to $\sim$32 %. This increase would allow characterisation of the full binary population by breaking degeneracies between suggested $P_{orb}$ distributions in the literature. Our results can be used in the ongoing assessment of the effectiveness of various potential cadencing strategies.

We present angular diameter distance measurements obtained by locating the BAO scale in the distribution of galaxies selected from the first year of Dark Energy Survey data. We consider a sample of over 1.3 million galaxies distributed over a footprint of 1318 deg$^2$ with $0.6 < z_{\rm photo} < 1$ and a typical redshift uncertainty of $0.03(1+z)$. This sample was selected, as fully described in a companion paper, using a color/magnitude selection that optimizes trade-offs between number density and redshift uncertainty. We investigate the BAO signal in the projected clustering using three conventions, the angular separation, the co-moving transverse separation, and spherical harmonics. Further, we compare results obtained from template based and machine learning photometric redshift determinations. We use 1800 simulations that approximate our sample in order to produce covariance matrices and allow us to validate our distance scale measurement methodology. We measure the angular diameter distance, $D_A$, at the effective redshift of our sample divided by the true physical scale of the BAO feature, $r_{\rm d}$. We obtain close to a 4 per cent distance measurement of $D_A(z_{\rm eff}=0.81)/r_{\rm d} = 10.75\pm 0.43 $. These results are consistent with the flat $\Lambda$CDM concordance cosmological model supported by numerous other recent experimental results. All data products are publicly available here: https://des.ncsa.illinois.edu/releases/y1a1/bao

Gravitational potentials which change in time induce fluctuations in the observed cosmic microwave background (CMB) temperature. Cosmological structure moving transverse to our line of sight provides a specific example known as the moving lens effect. Here we explore how the observed CMB temperature fluctuations combined with the observed matter over-density can be used to infer the transverse velocity of cosmological structure on large scales. We show that near-future CMB surveys and galaxy surveys will have the statistical power to make a first detection of the moving lens effect, and we discuss applications for the reconstructed transverse velocity.

We argue that the exclusion of the Galactic Plane and Bulge from the uniform wide-fast-deep (WFD) LSST survey cadence is fundamentally inconsistent with two of the main science drivers of LSST: Mapping the Milky Way and Exploring the Transient Optical Sky. We outline the philosophical basis for this claim and then describe a number of important science goals that can only be addressed by WFD-like coverage of the Plane and Bulge.

We present spectroscopy from the first three seasons of the Dark Energy Survey Supernova Program (DES-SN). We describe the supernova spectroscopic program in full: strategy, observations, data reduction, and classification. We have spectroscopically confirmed 307 supernovae, including 251 type Ia supernovae (SNe Ia) over a redshift range of $0.017 < z < 0.85$. We determine the effective spectroscopic selection function for our sample, and use it to investigate the redshift-dependent bias on the distance moduli of SNe Ia we have classified. We also provide a full overview of the strategy, observations, and data products of DES-SN, which has discovered 12,015 likely supernovae during these first three seasons. The data presented here are used for the first cosmology analysis by DES-SN ('DES-SN3YR'), the results of which are given in DES Collaboration (2018a).

We present spectroscopy from the first three seasons of the Dark Energy Survey Supernova Program (DES-SN). We describe the supernova spectroscopic program in full: strategy, observations, data reduction, and classification. We have spectroscopically confirmed 307 supernovae, including 251 type Ia supernovae (SNe Ia) over a redshift range of $0.017 < z < 0.85$. We determine the effective spectroscopic selection function for our sample, and use it to investigate the redshift-dependent bias on the distance moduli of SNe Ia we have classified. We also provide a full overview of the strategy, observations, and data products of DES-SN, which has discovered 12,015 likely supernovae during these first three seasons. The data presented here are used for the first cosmology analysis by DES-SN ('DES-SN3YR'), the results of which are given in DES Collaboration (2018a).

The Event Horizon Telescope (EHT), a global 230 GHz VLBI array, achieves an angular resolution of ${\approx} 20\,\mu{\rm as}$, sufficient to resolve the supermassive black hole Sagittarius A* (Sgr A*). This resolution may soon enable measurement of the black hole "shadow" size and asymmetry, predicted to be ${\approx}50\,\mu$as and ${\lesssim} 3\,\mu$as, respectively. Measurements that depart from these values could indicate a violation of the "no-hair theorem." However, refractive scattering by the turbulent ionized interstellar medium distorts the image of Sgr A*, affecting its apparent size and asymmetry. In this paper, we present a general analytic approach to quantify the expected image wander, distortion, and asymmetry from refractive scattering. If the turbulence in the scattering material of Sgr A* is close to Kolmogorov, we estimate the mean refractive image wander, distortion, and asymmetry to be 0.53 $\mu$as, 0.72 $\mu$as, and 0.52 $\mu$as at 230 GHz. However, alternative scattering models with flatter power spectra can yield larger values, up to 2.1 $\mu$as, 6.3 $\mu$as, and 5.0 $\mu$as, respectively. We demonstrate that these effects can be reduced by averaging images over multiple observations. For a small number of observations, the effects of scattering can be comparable to or greater than those from black hole spin, and they determine a fundamental limit for testing general relativity via images of Sgr A*.

Kinetic Sunyaev Zel'dovich (kSZ) tomography provides a powerful probe of the radial velocity field of matter in the Universe. By cross-correlating a high resolution CMB experiment like CMB S4 and a galaxy survey like DESI or LSST, one can measure the radial velocity field with very high signal to noise over a large volume of the universe. In this paper we show how this measurement can be used to improve constraints on primordial non-Gaussianities of the local type. The velocity field provides a measurement of the unbiased matter perturbations on large scales, which can be cross-correlated with the biased large-scale galaxy density field. This results in sample variance cancellation for a measurement of scale-dependent bias due to a non-zero $f_{NL}$. Using this method we forecast that CMB S4 and LSST combined reach a sensitivity $\sigma_{f_{NL}} \sim 0.5$, which is a factor of three improvement over the sensitivity using LSST alone (without internal sample variance cancellation). We take into account critical systematics like photometric redshifts, the kSZ optical depth degeneracy, and systematics affecting the shape of the galaxy auto-power spectrum and find that these have negligible impact, thus making kSZ tomography a robust probe for primordial non-Gaussianities. We also forecast the impact of mass binning on our constraints. The techniques proposed in this paper could be an important component of achieving the theoretically important threshold of $\sigma_{f_{NL}} \lesssim 1$ with future surveys.

Several statistics have been proposed for measuring the kSZ effect by combining the small-scale CMB with galaxy surveys. We review five such statistics, and show that they are all mathematically equivalent to the optimal bispectrum estimator of type $\langle ggT \rangle$. Reinterpreting these kSZ statistics as special cases of bispectrum estimation makes many aspects transparent, for example optimally weighting the estimator, or incorporating photometric redshift errors. We analyze the information content of the bispectrum and show that there are two observables: the small-scale galaxy-electron power spectrum $P_{ge}(k_S)$, and the large-scale galaxy-velocity power spectrum $P_{gv}(k)$. The cosmological constraining power of the kSZ arises from its sensitivity to fluctuations on large length scales, where its effective noise level can be much better than galaxy surveys.

The kinetic Sunyaev Zel'dovich (kSZ) and polarized Sunyaev Zel'dovich (pSZ) effects are temperature and polarization anisotropies induced by the scattering of CMB photons from structure in the post-reionization Universe. In the case of the kSZ effect, small angular scale anisotropies in the optical depth are modulated by the cosmic microwave background (CMB) dipole field, i.e. the CMB dipole observed at each spacetime point, which is sourced by the primordial dipole and especially the local peculiar velocity. In the case of the pSZ effect, similar small-scale anisotropies are modulated by the CMB quadrupole field, which receives contributions from both scalar and tensor modes. Statistical anisotropies in the cross correlations of CMB temperature and polarization with tracers of the inhomogeneous distribution of electrons provide a means of isolating and reconstructing the dipole and quadrupole fields. In this paper, we present a set of unbiased minimum variance quadratic estimators for the reconstruction of the dipole and quadrupole fields, and forecast the ability of future CMB experiments and large scale structure surveys to perform this reconstruction. Consistent with previous work, we find that a high fidelity reconstruction of the dipole and quadrupole fields over a variety of scales is indeed possible, and demonstrate the sensitivity of the pSZ effect to primordial tensor modes. Using a principle component analysis, we estimate how many independent modes could be accessed in such a reconstruction. We also comment on a few first applications of a detection of the dipole and quadrupole fields, including a reconstruction of the primordial contribution to our locally observed CMB dipole, a test of statistical homogeneity on large scales from the first modes of the quadrupole field, and a reconstruction technique for the primordial potential on the largest scales.

The detection and characterization of primordial gravitational waves through their impact on the polarization anisotropies of the cosmic microwave background (CMB) is a primary science goal of current and future observations of the CMB. An ancillary dataset that will become accessible with the great leaps in sensitivity of CMB experiments is the polarized Sunyaev Zel'dovich (pSZ) effect, small-scale CMB polarization anisotropies induced by scattering from free electrons in the post-reionization Universe. The cross correlation of the pSZ effect with galaxy surveys, a technique known as pSZ tomography, can be used to reconstruct the remote quadrupole field: the CMB quadrupole observed from different locations in the Universe. Primordial gravitational waves leave a distinct imprint on the remote quadrupole field, making pSZ tomography a potential new method to characterize their properties. Building on previous work, we explore the utility of the full set of correlations between the primary CMB and the reconstructed remote quadrupole field to both provide exclusion limits on the amplitude of primordial gravitational waves, as well as to provide constraints on several phenomenological models of the tensor sector: axion gauge field inflation, general models with chiral tensors, and models with modified late-time decay of tensors. We find that relatively futuristic experimental requirements are necessary to provide competitive exclusion limits compared with the primary CMB. However, pSZ tomography can be a powerful probe of the late-time evolution of tensors and, through cross-correlations with the primary CMB, can provide mild improvements on parameter constraints in various models with chiral primordial gravitational waves.

Millimeter polarimetry of Sgr A* probes the linearly polarized emission region on a scale of $\sim 10$ Schwarzschild radii ($R_S$) as well as the dense, magnetized accretion flow on scales out to the Bondi radius ($\sim 10^5 R_S$) through Faraday rotation. We present here multi-epoch ALMA Band 6 (230 GHz) polarimetry of Sgr A*. The results confirm a mean rotation measure, ${\rm RM} \approx -5 \times 10^5 {\rm\ rad\ m^{-2}}$, consistent with measurements over the past 20 years and support an interpretation of the RM as originating from a radiatively inefficient accretion flow (RIAF) with $\dot{M} \approx 10^{-8} { \rm\ M_{\odot}\ y^{-1} }$. Variability is observed for the first time in the RM on time scales that range from hours to months. The long-term variations may be the result of changes in the line of sight properties in a turbulent accretion flow. Short-term variations in the apparent RM are not necessarily the result of Faraday rotation and may be the result of complex emission and propagatation effects close to the black hole, some of which have been predicted in numerical modeling. We also confirm the detection of circular polarization at a mean value of $-1.1 \pm 0.2 \%$. It is variable in amplitude on time scales from hours to months but the handedness remains unchanged from that observed in past centimeter- and millimeter-wavelength detections. These results provide critical constraints for the analysis and interpretation of Event Horizon Telescope data of Sgr A*, M87, and similar sources.

We perform a joint analysis of the auto and cross-correlations between three cosmic fields: the galaxy density field, the galaxy weak lensing shear field, and the cosmic microwave background (CMB) weak lensing convergence field. These three fields are measured using roughly 1300 sq. deg. of overlapping optical imaging data from first year observations of the Dark Energy Survey and millimeter-wave observations of the CMB from both the South Pole Telescope Sunyaev-Zel'dovich survey and Planck. We present cosmological constraints from the joint analysis of the two-point correlation functions between galaxy density and galaxy shear with CMB lensing. We test for consistency between these measurements and the DES-only two-point function measurements, finding no evidence for inconsistency in the context of flat $\Lambda$CDM cosmological models. Performing a joint analysis of five of the possible correlation functions between these fields (excluding only the CMB lensing autospectrum) yields $S_{8}\equiv \sigma_8\sqrt{\Omega_{\rm m}/0.3} = 0.782^{+0.019}_{-0.025}$ and $\Omega_{\rm m}=0.260^{+0.029}_{-0.019}$. We test for consistency between these five correlation function measurements and the Planck-only measurement of the CMB lensing autospectrum, again finding no evidence for inconsistency in the context of flat $\Lambda$CDM models. Combining constraints from all six two-point functions yields $S_{8}=0.776^{+0.014}_{-0.021}$ and $\Omega_{\rm m}= 0.271^{+0.022}_{-0.016}$. These results provide a powerful test and confirmation of the results from the first year DES joint-probes analysis.

We perform a joint analysis of the auto and cross-correlations between three cosmic fields: the galaxy density field, the galaxy weak lensing shear field, and the cosmic microwave background (CMB) weak lensing convergence field. These three fields are measured using roughly 1300 sq. deg. of overlapping optical imaging data from first year observations of the Dark Energy Survey and millimeter-wave observations of the CMB from both the South Pole Telescope Sunyaev-Zel'dovich survey and Planck. We present cosmological constraints from the joint analysis of the two-point correlation functions between galaxy density and galaxy shear with CMB lensing. We test for consistency between these measurements and the DES-only two-point function measurements, finding no evidence for inconsistency in the context of flat $\Lambda$CDM cosmological models. Performing a joint analysis of five of the possible correlation functions between these fields (excluding only the CMB lensing autospectrum) yields $S_{8}\equiv \sigma_8\sqrt{\Omega_{\rm m}/0.3} = 0.782^{+0.019}_{-0.025}$ and $\Omega_{\rm m}=0.260^{+0.029}_{-0.019}$. We test for consistency between these five correlation function measurements and the Planck-only measurement of the CMB lensing autospectrum, again finding no evidence for inconsistency in the context of flat $\Lambda$CDM models. Combining constraints from all six two-point functions yields $S_{8}=0.776^{+0.014}_{-0.021}$ and $\Omega_{\rm m}= 0.271^{+0.022}_{-0.016}$. These results provide a powerful test and confirmation of the results from the first year DES joint-probes analysis.

The Kilodegree Extremely Little Telescope (KELT) project has been conducting a photometric survey for transiting planets orbiting bright stars for over ten years. The KELT images have a pixel scale of ~23"/pixel---very similar to that of NASA's Transiting Exoplanet Survey Satellite (TESS)---as well as a large point spread function, and the KELT reduction pipeline uses a weighted photometric aperture with radius 3'. At this angular scale, multiple stars are typically blended in the photometric apertures. In order to identify false positives and confirm transiting exoplanets, we have assembled a follow-up network (KELT-FUN) to conduct imaging with higher spatial resolution, cadence, and photometric precision than the KELT telescopes, as well as spectroscopic observations of the candidate host stars. The KELT-FUN team has followed-up over 1,600 planet candidates since 2011, resulting in more than 20 planet discoveries. Excluding ~450 false alarms of non-astrophysical origin (i.e., instrumental noise or systematics), we present an all-sky catalog of the 1,128 bright stars (6<V<10) that show transit-like features in the KELT light curves, but which were subsequently determined to be astrophysical false positives (FPs) after photometric and/or spectroscopic follow-up observations. The KELT-FUN team continues to pursue KELT and other planet candidates and will eventually follow up certain classes of TESS candidates. The KELT FP catalog will help minimize the duplication of follow-up observations by current and future transit surveys such as TESS.

We present the discovery and confirmation of two new transiting giant planets from the Kepler extended mission K2. K2-260 b is a hot Jupiter transiting a $V=12.7$ F6V star in K2 Field 13, with a mass and radius of $M_{\star}=1.39_{-0.06}^{+0.05} M_{\odot}$ and $R_{\star}=1.69 \pm 0.03 R_{\odot}$. The planet has an orbital period of $P=2.627$ days, and a mass and radius of $M_P=1.42^{+0.31}_{-0.32} M_J$ and $R_P=1.552^{+0.048}_{-0.057} R_J$. This is the first K2 hot Jupiter with a detected secondary eclipse in the Kepler bandpass, with a depth of $71 \pm 15$ ppm, which we use to estimate a geometric albedo of $A_g\sim0.2$. We also detected a candidate stellar companion at 0.6" from K2-260; we find that it is very likely physically associated with the system, in which case it would be an M5-6V star at a projected separation of $\sim400$ AU. K2-261 b is a warm Saturn transiting a bright ($V=10.5$) G7IV/V star in K2 Field 14. The host star is a metal-rich ([Fe/H]$=0.36 \pm 0.06$), mildly evolved $1.10_{-0.02}^{+0.01} M_{\odot}$ star with $R_{\star}=1.65 \pm 0.04 R_{\odot}$. Thanks to its location near the main sequence turn-off, we can measure a relatively precise age of $8.8_{-0.3}^{+0.4}$ Gyr. The planet has $P=11.633$ days, $M_P=0.223 \pm 0.031 M_J$, and $R_P=0.850^{+0.026}_{-0.022} R_J$, and its orbit is eccentric ($e=0.39 \pm 0.15$). Its brightness and relatively large transit depth make this one of the best known warm Saturns for follow-up observations to further characterize the planetary system.

We present two galaxy shape catalogues from the Dark Energy Survey Year 1 data set, covering 1500 square degrees with a median redshift of $0.59$. The catalogues cover two main fields: Stripe 82, and an area overlapping the South Pole Telescope survey region. We describe our data analysis process and in particular our shape measurement using two independent shear measurement pipelines, METACALIBRATION and IM3SHAPE. The METACALIBRATION catalogue uses a Gaussian model with an innovative internal calibration scheme, and was applied to $riz$-bands, yielding 34.8M objects. The IM3SHAPE catalogue uses a maximum-likelihood bulge/disc model calibrated using simulations, and was applied to $r$-band data, yielding 21.9M objects. Both catalogues pass a suite of null tests that demonstrate their fitness for use in weak lensing science. We estimate the 1$\sigma$ uncertainties in multiplicative shear calibration to be $0.013$ and $0.025$ for the METACALIBRATION and IM3SHAPE catalogues, respectively.

We present two galaxy shape catalogues from the Dark Energy Survey Year 1 data set, covering 1500 square degrees with a median redshift of $0.59$. The catalogues cover two main fields: Stripe 82, and an area overlapping the South Pole Telescope survey region. We describe our data analysis process and in particular our shape measurement using two independent shear measurement pipelines, METACALIBRATION and IM3SHAPE. The METACALIBRATION catalogue uses a Gaussian model with an innovative internal calibration scheme, and was applied to $riz$-bands, yielding 34.8M objects. The IM3SHAPE catalogue uses a maximum-likelihood bulge/disc model calibrated using simulations, and was applied to $r$-band data, yielding 21.9M objects. Both catalogues pass a suite of null tests that demonstrate their fitness for use in weak lensing science. We estimate the 1$\sigma$ uncertainties in multiplicative shear calibration to be $0.013$ and $0.025$ for the METACALIBRATION and IM3SHAPE catalogues, respectively.

We present galaxy-galaxy lensing measurements from 1321 sq. deg. of the Dark Energy Survey (DES) Year 1 (Y1) data. The lens sample consists of a selection of 660,000 red galaxies with high-precision photometric redshifts, known as redMaGiC, split into five tomographic bins in the redshift range $0.15 < z < 0.9$. We use two different source samples, obtained from the Metacalibration (26 million galaxies) and Im3shape (18 million galaxies) shear estimation codes, which are split into four photometric redshift bins in the range $0.2 < z < 1.3$. We perform extensive testing of potential systematic effects that can bias the galaxy-galaxy lensing signal, including those from shear estimation, photometric redshifts, and observational properties. Covariances are obtained from jackknife subsamples of the data and validated with a suite of log-normal simulations. We use the shear-ratio geometric test to obtain independent constraints on the mean of the source redshift distributions, providing validation of those obtained from other photo-$z$ studies with the same data. We find consistency between the galaxy bias estimates obtained from our galaxy-galaxy lensing measurements and from galaxy clustering, therefore showing the galaxy-matter cross-correlation coefficient $r$ to be consistent with one, measured over the scales used for the cosmological analysis. The results in this work present one of the three two-point correlation functions, along with galaxy clustering and cosmic shear, used in the DES cosmological analysis of Y1 data, and hence the methodology and the systematics tests presented here provide a critical input for that study as well as for future cosmological analyses in DES and other photometric galaxy surveys.

We present galaxy-galaxy lensing measurements from 1321 sq. deg. of the Dark Energy Survey (DES) Year 1 (Y1) data. The lens sample consists of a selection of 660,000 red galaxies with high-precision photometric redshifts, known as redMaGiC, split into five tomographic bins in the redshift range $0.15 < z < 0.9$. We use two different source samples, obtained from the Metacalibration (26 million galaxies) and Im3shape (18 million galaxies) shear estimation codes, which are split into four photometric redshift bins in the range $0.2 < z < 1.3$. We perform extensive testing of potential systematic effects that can bias the galaxy-galaxy lensing signal, including those from shear estimation, photometric redshifts, and observational properties. Covariances are obtained from jackknife subsamples of the data and validated with a suite of log-normal simulations. We use the shear-ratio geometric test to obtain independent constraints on the mean of the source redshift distributions, providing validation of those obtained from other photo-$z$ studies with the same data. We find consistency between the galaxy bias estimates obtained from our galaxy-galaxy lensing measurements and from galaxy clustering, therefore showing the galaxy-matter cross-correlation coefficient $r$ to be consistent with one, measured over the scales used for the cosmological analysis. The results in this work present one of the three two-point correlation functions, along with galaxy clustering and cosmic shear, used in the DES cosmological analysis of Y1 data, and hence the methodology and the systematics tests presented here provide a critical input for that study as well as for future cosmological analyses in DES and other photometric galaxy surveys.

We report measurements of the sky-projected spin--orbit angles for three transiting hot Jupiters: two of which are in nearly polar orbits, WASP-100b and WASP-109b, and a third in a low obliquity orbit, WASP-72b. We obtained these measurements by observing the Rossiter--McLaughlin effect over the course of the transits from high resolution spectroscopic observations made with the CYCLOPS2 optical fiber bundle system feeding the UCLES spectrograph on the Anglo-Australian Telescope. The resulting sky-projected spin--orbit angles are $\lambda = {-7^{\circ}}^{+11^{\circ}}_{-12^{\circ}}$, $\lambda = {79^{\circ}}^{+19^{\circ}}_{-10^{\circ}}$, and $\lambda = {99^{\circ}}^{+10^{\circ}}_{-9^{\circ}}$ for WASP-72b, WASP-100b, and WASP-109b, respectively. These results suggests that WASP-100b and WASP-109b are on highly inclined orbits tilted nearly $90^{\circ}$ from their host star's equator while the orbit of WASP-72b appears to be well-aligned. WASP-72b is a typical hot Jupiter orbiting a mid-late F star (F7 with $T_\mathrm{{eff}}=6250\pm120$K). WASP-100b and WASP-109b are highly irradiated bloated hot Jupiters orbiting hot early-mid F stars (F2 with $T_\mathrm{{eff}}=6900\pm120$K and F4 with $T_\mathrm{{eff}}=6520\pm140$K), making them consistent with the trends observed for the majority of stars hosting planets on high-obliquity orbits.

We measure the clustering of DES Year 1 galaxies that are intended to be combined with weak lensing samples in order to produce precise cosmological constraints from the joint analysis of large-scale structure and lensing correlations. Two-point correlation functions are measured for a sample of $6.6 \times 10^{5}$ luminous red galaxies selected using the \textsc{redMaGiC} algorithm over an area of $1321$ square degrees, in the redshift range $0.15 < z < 0.9$, split into five tomographic redshift bins. The sample has a mean redshift uncertainty of $\sigma_{z}/(1+z) = 0.017$. We quantify and correct spurious correlations induced by spatially variable survey properties, testing their impact on the clustering measurements and covariance. We demonstrate the sample's robustness by testing for stellar contamination, for potential biases that could arise from the systematic correction, and for the consistency between the two-point auto- and cross-correlation functions. We show that the corrections we apply have a significant impact on the resultant measurement of cosmological parameters, but that the results are robust against arbitrary choices in the correction method. We find the linear galaxy bias in each redshift bin in a fiducial cosmology to be $b(z$=$0.24)=1.40 \pm 0.08$, $b(z$=$0.38)=1.61 \pm 0.05$, $b(z$=$0.53)=1.60 \pm 0.04$ for galaxies with luminosities $L/L_*>$$0.5$, $b(z$=$0.68)=1.93 \pm 0.05$ for $L/L_*>$$1$ and $b(z$=$0.83)=1.99 \pm 0.07$ for $L/L_*$$>1.5$, broadly consistent with expectations for the redshift and luminosity dependence of the bias of red galaxies. We show these measurements to be consistent with the linear bias obtained from tangential shear measurements.

We measure the clustering of DES Year 1 galaxies that are intended to be combined with weak lensing samples in order to produce precise cosmological constraints from the joint analysis of large-scale structure and lensing correlations. Two-point correlation functions are measured for a sample of $6.6 \times 10^{5}$ luminous red galaxies selected using the \textsc{redMaGiC} algorithm over an area of $1321$ square degrees, in the redshift range $0.15 < z < 0.9$, split into five tomographic redshift bins. The sample has a mean redshift uncertainty of $\sigma_{z}/(1+z) = 0.017$. We quantify and correct spurious correlations induced by spatially variable survey properties, testing their impact on the clustering measurements and covariance. We demonstrate the sample's robustness by testing for stellar contamination, for potential biases that could arise from the systematic correction, and for the consistency between the two-point auto- and cross-correlation functions. We show that the corrections we apply have a significant impact on the resultant measurement of cosmological parameters, but that the results are robust against arbitrary choices in the correction method. We find the linear galaxy bias in each redshift bin in a fiducial cosmology to be $b(z$=$0.24)=1.40 \pm 0.08$, $b(z$=$0.38)=1.61 \pm 0.05$, $b(z$=$0.53)=1.60 \pm 0.04$ for galaxies with luminosities $L/L_*>$$0.5$, $b(z$=$0.68)=1.93 \pm 0.05$ for $L/L_*>$$1$ and $b(z$=$0.83)=1.99 \pm 0.07$ for $L/L_*$$>1.5$, broadly consistent with expectations for the redshift and luminosity dependence of the bias of red galaxies. We show these measurements to be consistent with the linear bias obtained from tangential shear measurements.

Radio images of the Galactic Center supermassive black hole, Sagittarius A* (Sgr A*), are dominated by interstellar scattering. Previous studies of Sgr A* have adopted an anisotropic Gaussian model for both the intrinsic source and the scattering, and they have extrapolated the scattering using a purely $\lambda^2$ scaling to estimate intrinsic properties. However, physically motivated source and scattering models break all three of these assumptions. They also predict that refractive scattering effects will be significant, which have been ignored in standard model fitting procedures. We analyze radio observations of Sgr A* using a physically motivated scattering model, and we develop a prescription to incorporate refractive scattering uncertainties when model fitting. We show that an anisotropic Gaussian scattering kernel is an excellent approximation for Sgr A* at wavelengths longer than 1cm, with an angular size of $(1.380 \pm 0.013) \lambda_{\rm cm}^2\,{\rm mas}$ along the major axis, $(0.703 \pm 0.013) \lambda_{\rm cm}^2\,{\rm mas}$ along the minor axis, and a position angle of $81.9^\circ \pm 0.2^\circ$. We estimate that the turbulent dissipation scale is at least $600\,{\rm km}$, with tentative support for $r_{\rm in} = 800 \pm 200\,{\rm km}$, suggesting that the ion Larmor radius defines the dissipation scale. We find that the power-law index for density fluctuations in the scattering material is $\beta < 3.47$, shallower than expected for a Kolmogorov spectrum ($\beta=11/3$), and we estimate $\beta = 3.38^{+0.08}_{-0.04}$ in the case of $r_{\rm in} = 800\,{\rm km}$. We find that the intrinsic structure of Sgr A* is nearly isotropic over wavelengths from 1.3mm to 1.3cm, with a size that is roughly proportional to wavelength. We discuss implications for models of Sgr A*, for theories of interstellar turbulence, and for imaging Sgr A* with the Event Horizon Telescope.

We report the discovery and characterization of an eclipsing M5V dwarf star, orbiting a slightly evolved F7V main sequence star. In contrast to previous claims in the literature, we confirm that the system does not belong to the galactic open cluster Ruprecht 147. We determine its fundamental parameters combining K2 time-series data with spectroscopic observations from the McDonald Observatory, FIES@NOT, and HIRES@KECK. The very precise photometric data from the K2 mission allows us to measure variations caused by the beaming effect (relativistic doppler boosting), ellipsoidal variation, reflection, and the secondary eclipse. We determined the radial velocity using spectroscopic observations and compare it to the radial velocity determined from the beaming effect observed in the photometric data. The M5V star has a radius of $0.200 \substack{+0.007 \\ -0.008}$ $R_{\odot}$ and a mass of $0.187 \substack{+0.012 \\ -0.013}$ $M_{\odot}$. The primary star has radius of $1.518 \substack{+0.038 \\ -0.049}$ $R_{\odot}$ and a mass of $1.008 \substack{+0.081 \\ -0.097}$ $M_{\odot}$. The orbital period is $ 5.441995 \pm 0.000007$ days. The system is one of the few eclipsing systems with observed beaming effect and spectroscopic radial velocity measurements and it can be used as test case for the modelling of the beaming effect. Current and forthcoming space missions such as TESS and PLATO might benefit of the analysis of the beaming effect to estimate the mass of transiting companions without the need for radial velocity follow up observations, provided that the systematic sources of noise affecting this method are well understood.

We present new L-band ($1\mbox{--}2$ GHz) observations of the nearby spiral galaxy M33 with 80 pc resolution obtained with the Karl G. Jansky Very Large Array. The HI observations, combined with HI measurements from the Green Bank Telescope, improve the spectral resolution and sensitivity ($2.8$ K rms noise in a $0.2$ km s$^{-1}$ channel) compared to previous observations. We find individual profiles are usually non-Gaussian, harbouring line wings, multiple components, and asymmetries. Given this spectral complexity, we quantify the motions in the atomic ISM through moment analysis of the spectra and fits to aligned, stacked profiles. The measured value of the HI line width depends strongly on the method used, with the velocity stacked profiles aligned to the peak velocity giving the minimum value of $\sigma = 7$ km s$^{-1}$ and all other methods giving higher values ($\sigma\sim10$ km s$^{-1}$). All measurements of the line width show a shallow radial trend, with $\sigma$ decreasing by $\sim2$ km s$^{-1}$ from $R_{\rm gal}=0$ to $R_{\rm gal}=8$ kpc. We consider a number of energy sources that might maintain the line width against turbulent dissipation, but no single source is adequate. We find excess emission relative to a Gaussian in the stacked profile line wings, ranging from 9% to 26% depending on how the spectra are aligned. By splitting the line wings into symmetric and asymmetric components, we find that the lagging rotational disk accounts for one-third of the line wing flux. We also find emission far from the rotation axis of the galaxy in multiple discrete HI clouds, including a filament with a projected length of $\sim8$ kpc.

We use data from the first-year (Y1) observations of the DES collaboration to measure the galaxy angular power spectrum (APS), and search for its BAO feature using a template-fitting method. We test our methodology in a sample of 1800 DES Y1-like mock catalogs. The APS is measured with the pseudo-$C_\ell$ method, using pixelized maps constructed from the mock catalogs and the DES mask. The covariance matrix of the $C_\ell$'s in these tests are also obtained from the mock catalogs. We use templates to model the measured spectra and estimate template parameters firstly from the $C_\ell$'s of the mocks using two different methods, a maximum likelihood estimator and a MCMC, finding consistent results with a good reduced $\chi^2$. Robustness tests are performed to estimate the impact of different choices of settings used in our analysis. After these tests on mocks, we apply our method to a galaxy sample constructed from DES Y1 data specifically for LSS studies. This catalog comprises galaxies within an effective area of 1318 deg$^2$ and $0.6<z<1.0$. We fit the observed spectra with our optimized templates, considering models with and without BAO features. We find that the DES Y1 data favors a model with BAO at the $2.6\,\sigma$ C.L. with a best-fit shift parameter of $\alpha=1.023\pm 0.047$. However, the goodness-of-fit is somewhat poor, with $\chi^2/$(dof) = 1.49. We identify a possible cause of this issue and show that using a theoretical covariance matrix obtained from $C_\ell$'s that are better adjusted to data results in an improved value of $\chi^2/$(dof) = 1.36 which is similar to the value obtained with the real-space analysis. Our results correspond to a distance measurement of $D_A(z_{\rm eff} = 0.81)/r_d = 10.65 \pm 0.49$, consistent with the main DES BAO findings. This is a companion paper to the main DES BAO article showing the details of the harmonic space analysis.

We use data from the first-year (Y1) observations of the DES collaboration to measure the galaxy angular power spectrum (APS), and search for its BAO feature using a template-fitting method. We test our methodology in a sample of 1800 DES Y1-like mock catalogs. The APS is measured with the pseudo-$C_\ell$ method, using pixelized maps constructed from the mock catalogs and the DES mask. The covariance matrix of the $C_\ell$'s in these tests are also obtained from the mock catalogs. We use templates to model the measured spectra and estimate template parameters firstly from the $C_\ell$'s of the mocks using two different methods, a maximum likelihood estimator and a MCMC, finding consistent results with a good reduced $\chi^2$. Robustness tests are performed to estimate the impact of different choices of settings used in our analysis. After these tests on mocks, we apply our method to a galaxy sample constructed from DES Y1 data specifically for LSS studies. This catalog comprises galaxies within an effective area of 1318 deg$^2$ and $0.6<z<1.0$. We fit the observed spectra with our optimized templates, considering models with and without BAO features. We find that the DES Y1 data favors a model with BAO at the $2.6\,\sigma$ C.L. with a best-fit shift parameter of $\alpha=1.023\pm 0.047$. However, the goodness-of-fit is somewhat poor, with $\chi^2/$(dof) = 1.49. We identify a possible cause of this issue and show that using a theoretical covariance matrix obtained from $C_\ell$'s that are better adjusted to data results in an improved value of $\chi^2/$(dof) = 1.36 which is similar to the value obtained with the real-space analysis. Our results correspond to a distance measurement of $D_A(z_{\rm eff} = 0.81)/r_d = 10.65 \pm 0.49$, consistent with the main DES BAO findings. This is a companion paper to the main DES BAO article showing the details of the harmonic space analysis.

Analog condensed matter systems present an exciting opportunity to simulate early Universe models in table-top experiments. We consider a recent proposal for an analog condensed matter experiment to simulate the relativistic quantum decay of the false vacuum. In the proposed experiment, two ultra-cold condensates are coupled via a time-varying radio-frequency field. The relative phase of the two condensates in this system is approximately described by a relativistic scalar field with a potential possessing a series of false and true vacuum local minima. If the system is set up in a false vacuum, it would then decay to a true vacuum via quantum mechanical tunnelling. Should such an experiment be realized, it would be possible to answer a number of open questions regarding non-perturbative phenomena in quantum field theory and early Universe cosmology. In this paper, we illustrate a possible obstruction: the time-varying coupling that is invoked to create a false vacuum for the long-wavelength modes of the condensate leads to a destabilization of shorter wavelength modes within the system via parametric resonance. We focus on an idealized setup in which the two condensates have identical properties and identical background densities. Describing the system by the coupled Gross-Pitaevskii equations (GPE), we use the machinery of Floquet theory to perform a linear stability analysis, calculating the wavenumber associated with the first instability band for a variety of experimental parameters. However, we demonstrate that, by tuning the frequency of the time-varying coupling, it may be possible to push the first instability band outside the validity of the GPE, where dissipative effects are expected to damp any instabilities. This provides a viable range of experimental parameters to perform analog experiments of false vacuum decay.

The kinetic Sunyaev Zel'dovich (kSZ) effect, cosmic microwave background (CMB) anisotropies induced by scattering from free electrons in bulk motion, is a primary target of future CMB experiments. Measurements of the kSZ effect have the potential to address fundamental questions about the structure and evolution of our Universe on the largest scales and at the earliest times. This potential is unlocked by combining measurements of small-scale CMB anisotropies with large-scale structure surveys, a technique known as kSZ tomography. Previous work established a quadratic estimator for the remote dipole field, the CMB dipole observed at different locations in the Universe. This previous work did not include gravitational lensing, redshift space distortions, or non-linear evolution of structure. In this paper, we investigate how well the remote dipole field can be reconstructed in the presence of such effects by using mock data from a suite of simulations. To properly model both large and small scales, we develop a novel box-in-box simulation pipeline, where small-scale information is obtained from N-body simulations, and large-scale information obtained by evolving fields using linear theory and adding the resulting corrections to the N-body particle data. This pipeline allows us to create properly correlated maps of the primary CMB including lensing, as well as the kSZ effect and density maps on the past light cone of an observer. Analyzing an ensemble of mocks, we find that the dipole field can be reconstructed with high fidelity over a range of angular scales and redshift bins, although there is evidence of excess power from nonlinear structure. We also analyze correlations with the primary CMB, investigating the ability of kSZ tomography to reconstruct the intrinsic CMB dipole. Our results constitute a proof-of-principle that kSZ tomography is a promising technique for future datasets.

We present results from spectroscopic follow-up observations of stars identified in the Kepler field and carried out by teams of the Kepler Follow-Up Observation Program. Two samples of stars were observed over six years (2009-2015): 614 standard stars (divided into "platinum" and "gold" categories) selected based on their asteroseismic detections and 2667 host stars of Kepler Objects of Interest (KOIs), most of them planet candidates. Four data analysis pipelines were used to derive stellar parameters for the observed stars. We compare the $T_{\mathrm{eff}}$, $\log$(g), and [Fe/H] values derived for the same stars by different pipelines; from the average of the standard deviations of the differences in these parameter values, we derive error floors of $\sim$ 100 K, 0.2 dex, and 0.1 dex for $T_{\mathrm{eff}}$, $\log$(g), and [Fe/H], respectively. Noticeable disagreements are seen mostly at the largest and smallest parameter values (e.g., in the giant star regime). Most of the $\log$(g) values derived from spectra for the platinum stars agree on average within 0.025 dex (but with a spread of 0.1-0.2 dex) with the asteroseismic $\log$(g) values. Compared to the Kepler Input Catalog (KIC), the spectroscopically derived stellar parameters agree within the uncertainties of the KIC, but are more precise and are thus an important contribution towards deriving more reliable planetary radii.

The accretion flow around the Galactic Center black hole Sagittarius A* (Sgr A*) is expected to have an electron temperature that is distinct from the ion temperature, due to weak Coulomb coupling in the low-density plasma. We present four two-temperature general relativistic radiative magnetohydrodynamic (GRRMHD) simulations of Sgr A* performed with the code KORAL. These simulations use different electron heating prescriptions, motivated by different models of the underlying plasma microphysics. We compare the Landau-damped turbulent cascade model used in previous work with a new prescription we introduce based on the results of particle-in-cell simulations of magnetic reconnection. With the turbulent heating model, electrons are preferentially heated in the polar outflow, whereas with the reconnection model electrons are heated by nearly the same fraction everywhere in the accretion flow. The spectra of the two models are similar around the submillimetre synchrotron peak, but the models heated by magnetic reconnection produce variability more consistent with the level observed from Sgr A*. All models produce 230~GHz images with distinct black hole shadows which are consistent with the image size measured by the Event Horizon Telescope, but only the turbulent heating produces an anisotropic `disc-jet' structure where the image is dominated by a polar outflow or jet at frequencies below the synchrotron peak. None of our models can reproduce the observed radio spectral slope, the large near-infrared and X-ray flares, or the near-infrared spectral index, all of which suggest non-thermal electrons are needed to fully explain the emission from Sgr A*.

We report results from very long baseline interferometric (VLBI) observations of the supermassive black hole in the Galactic center, Sgr A*, at 1.3 mm (230 GHz). The observations were performed in 2013 March using six VLBI stations in Hawaii, California, Arizona, and Chile. Compared to earlier observations, the addition of the APEX telescope in Chile almost doubles the longest baseline length in the array, provides additional {\it uv} coverage in the N-S direction, and leads to a spatial resolution of $\sim$30 $\mu$as ($\sim$3 Schwarzschild radii) for Sgr A*. The source is detected even at the longest baselines with visibility amplitudes of $\sim$4-13% of the total flux density. We argue that such flux densities cannot result from interstellar refractive scattering alone, but indicate the presence of compact intrinsic source structure on scales of $\sim$3 Schwarzschild radii. The measured nonzero closure phases rule out point-symmetric emission. We discuss our results in the context of simple geometric models that capture the basic characteristics and brightness distributions of disk- and jet-dominated models and show that both can reproduce the observed data. Common to these models are the brightness asymmetry, the orientation, and characteristic sizes, which are comparable to the expected size of the black hole shadow. Future 1.3 mm VLBI observations with an expanded array and better sensitivity will allow a more detailed imaging of the horizon-scale structure and bear the potential for a deep insight into the physical processes at the black hole boundary.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

The potential habitability of an exoplanet is traditionally assessed by determining if its orbit falls within the circumstellar `habitable zone' of its star, defined as the distance at which water could be liquid on the surface of a planet (Kopparapu et al., 2013). Traditionally, these limits are determined by radiative-convective climate models, which are used to predict surface temperatures at user-specified levels of greenhouse gases. This approach ignores the vital question of the (bio)geochemical plausibility of the proposed chemical abundances. Carbon dioxide is the most important greenhouse gas in Earth's atmosphere in terms of regulating planetary temperature, with the long term concentration controlled by the balance between volcanic outgassing and the sequestration of CO2 via chemical weathering and sedimentation, as modulated by ocean chemistry, circulation and biological (microbial) productivity. We develop a model incorporating key aspects of Earth's short and long-term biogeochemical carbon cycle to explore the potential changes in the CO2 greenhouse due to variance in planet size and stellar insolation. We find that proposed changes in global topography, tectonics, and the hydrological cycle on larger planets results in proportionally greater surface temperatures for a given incident flux. For planets between 0.5 to 2 R_earth the effect of these changes results in average global surface temperature deviations of up to 20 K, which suggests that these relationships must be considered in future studies of planetary habitability.

We describe the derivation and validation of redshift distribution estimates and their uncertainties for the galaxies used as weak lensing sources in the Dark Energy Survey (DES) Year 1 cosmological analyses. The Bayesian Photometric Redshift (BPZ) code is used to assign galaxies to four redshift bins between z=0.2 and 1.3, and to produce initial estimates of the lensing-weighted redshift distributions $n^i_{PZ}(z)$ for bin i. Accurate determination of cosmological parameters depends critically on knowledge of $n^i$ but is insensitive to bin assignments or redshift errors for individual galaxies. The cosmological analyses allow for shifts $n^i(z)=n^i_{PZ}(z-\Delta z^i)$ to correct the mean redshift of $n^i(z)$ for biases in $n^i_{\rm PZ}$. The $\Delta z^i$ are constrained by comparison of independently estimated 30-band photometric redshifts of galaxies in the COSMOS field to BPZ estimates made from the DES griz fluxes, for a sample matched in fluxes, pre-seeing size, and lensing weight to the DES weak-lensing sources. In companion papers, the $\Delta z^i$ are further constrained by the angular clustering of the source galaxies around red galaxies with secure photometric redshifts at 0.15<z<0.9. This paper details the BPZ and COSMOS procedures, and demonstrates that the cosmological inference is insensitive to details of the $n^i(z)$ beyond the choice of $\Delta z^i$. The clustering and COSMOS validation methods produce consistent estimates of $\Delta z^i$, with combined uncertainties of $\sigma_{\Delta z^i}=$0.015, 0.013, 0.011, and 0.022 in the four bins. We marginalize over these in all analyses to follow, which does not diminish the constraining power significantly. Repeating the photo-z procedure using the Directional Neighborhood Fitting (DNF) algorithm instead of BPZ, or using the $n^i(z)$ directly estimated from COSMOS, yields no discernible difference in cosmological inferences.

We describe the derivation and validation of redshift distribution estimates and their uncertainties for the galaxies used as weak lensing sources in the Dark Energy Survey (DES) Year 1 cosmological analyses. The Bayesian Photometric Redshift (BPZ) code is used to assign galaxies to four redshift bins between z=0.2 and 1.3, and to produce initial estimates of the lensing-weighted redshift distributions $n^i_{PZ}(z)$ for bin i. Accurate determination of cosmological parameters depends critically on knowledge of $n^i$ but is insensitive to bin assignments or redshift errors for individual galaxies. The cosmological analyses allow for shifts $n^i(z)=n^i_{PZ}(z-\Delta z^i)$ to correct the mean redshift of $n^i(z)$ for biases in $n^i_{\rm PZ}$. The $\Delta z^i$ are constrained by comparison of independently estimated 30-band photometric redshifts of galaxies in the COSMOS field to BPZ estimates made from the DES griz fluxes, for a sample matched in fluxes, pre-seeing size, and lensing weight to the DES weak-lensing sources. In companion papers, the $\Delta z^i$ are further constrained by the angular clustering of the source galaxies around red galaxies with secure photometric redshifts at 0.15<z<0.9. This paper details the BPZ and COSMOS procedures, and demonstrates that the cosmological inference is insensitive to details of the $n^i(z)$ beyond the choice of $\Delta z^i$. The clustering and COSMOS validation methods produce consistent estimates of $\Delta z^i$, with combined uncertainties of $\sigma_{\Delta z^i}=$0.015, 0.013, 0.011, and 0.022 in the four bins. We marginalize over these in all analyses to follow, which does not diminish the constraining power significantly. Repeating the photo-z procedure using the Directional Neighborhood Fitting (DNF) algorithm instead of BPZ, or using the $n^i(z)$ directly estimated from COSMOS, yields no discernible difference in cosmological inferences.

Scattering in the ionized interstellar medium is commonly observed to be anisotropic, with theories of magnetohydrodynamic (MHD) turbulence explaining the anisotropy through a preferred magnetic field direction throughout the scattering regions. In particular, the line of sight to the Galactic Center supermassive black hole, Sgr A*, exhibits strong and anisotropic scattering, which dominates its observed size at wavelengths of a few millimeters and longer. Therefore, inferences of the intrinsic structure of \sgra\ at these wavelengths are sensitive to the assumed scattering model. In addition, extrapolations of the scattering model from long wavelengths, at which its parameters are usually estimated, to 1.3 mm, where the Event Horizon Telescope (EHT) seeks to image Sgr A* on Schwarzschild-radius scales, are also sensitive to the assumed scattering model. Past studies of Sgr A* have relied on simple Gaussian models for the scattering kernel that effectively presume an inner scale of turbulence far greater than the diffractive scale; this assumption is likely violated for Sgr A* at 1.3 mm. We develop a physically motivated model for anisotropic scattering, using a simplified model for MHD turbulence with a finite inner scale and a wandering transverse magnetic field direction. We explore several explicit analytic models for this wandering and derive the expected observational properties --- scatter broadening and refractive scintillation --- for each. For expected values of the inner scale, the scattering kernel for all models is markedly non-Gaussian at 1.3 mm but is straightforward to calculate and depends only weakly on the assumed model for the wandering of the magnetic field direction. On the other hand, in all models, the refractive substructure depends strongly on the wandering model and may be an important consideration in imaging Sgr A* with the EHT.