The Last Journey is a large-volume, gravity-only, cosmological N-body simulation evolving more than 1.24 trillion particles in a periodic box with a side-length of 5.025Gpc. It was implemented using the HACC simulation and analysis framework on the BG/Q system, Mira. The cosmological parameters are chosen to be consistent with the results from the Planck satellite. A range of analysis tools have been run in situ to enable a diverse set of science projects, and at the same time, to keep the resulting data amount manageable. Analysis outputs have been generated starting at redshift z~10 to allow for construction of synthetic galaxy catalogs using a semi-analytic modeling approach in post-processing. As part of our in situ analysis pipeline we employ a new method for tracking halo sub-structures, introducing the concept of subhalo cores. The production of multi-wavelength synthetic sky maps is facilitated by generating particle lightcones in situ, also beginning at z~10. We provide an overview of the simulation set-up and the generated data products; a first set of analysis results is presented. A subset of the data is publicly available.

We describe the first major public data release from cosmological simulations carried out with Argonne's HACC code. This initial release covers a range of datasets from large gravity-only simulations. The data products include halo information for multiple redshifts, down-sampled particles, and lightcone outputs. We provide data from two very large LCDM simulations as well as beyond-LCDM simulations spanning eleven w0-wa cosmologies. Our release platform uses Petrel, a research data service, located at the Argonne Leadership Computing Facility. Petrel offers fast data transfer mechanisms and authentication via Globus, enabling simple and efficient access to stored datasets. Easy browsing of the available data products is provided via a web portal that allows the user to navigate simulation products efficiently. The data hub will be extended by adding more types of data products and by enabling computational capabilities to allow direct interactions with simulation results.

We perform a joint analysis of intrinsic alignments and cosmology using tomographic weak lensing, galaxy clustering and galaxy-galaxy lensing measurements from Year 1 (Y1) of the Dark Energy Survey. We define early- and late-type subsamples, which are found to pass a series of systematics tests, including for spurious photometric redshift error and point spread function correlations. We analyse these split data alongside the fiducial mixed Y1 sample using a range of intrinsic alignment models. In a fiducial Nonlinear Alignment Model (NLA) analysis, assuming a flat \lcdm~cosmology, we find a significant difference in intrinsic alignment amplitude, with early-type galaxies favouring $A_\mathrm{IA} = 2.38^{+0.32}_{-0.31}$ and late-type galaxies consistent with no intrinsic alignments at $0.05^{+0.10}_{-0.09}$. We find weak evidence of a diminishing alignment amplitude at higher redshifts in the early-type sample. The analysis is repeated using a number of extended model spaces, including a physically motivated model that includes both tidal torquing and tidal alignment mechanisms. In multiprobe likelihood chains in which cosmology, intrinsic alignments in both galaxy samples and all other relevant systematics are varied simultaneously, we find the tidal alignment and tidal torquing parts of the intrinsic alignment signal have amplitudes $A_1 = 2.66 ^{+0.67}_{-0.66}$, $A_2=-2.94^{+1.94}_{-1.83}$, respectively, for early-type galaxies and $A_1 = 0.62 ^{+0.41}_{-0.41}$, $A_2 = -2.26^{+1.30}_{-1.16}$ for late-type galaxies. In the full (mixed) Y1 sample the best constraints are $A_1 = 0.70 ^{+0.41}_{-0.38}$, $A_2 = -1.36 ^{+1.08}_{-1.41}$. For all galaxy splits and IA models considered, we report cosmological parameter constraints that are consistent with the results of Troxel et al. (2017) and Dark Energy Survey Collaboration (2017).

This paper introduces cosmoDC2, a large synthetic galaxy catalog designed to support precision dark energy science with the Large Synoptic Survey Telescope (LSST). CosmoDC2 is the starting point for the second data challenge (DC2) carried out by the LSST Dark Energy Science Collaboration (LSST DESC). The catalog is based on a trillion-particle, 4.225 Gpc^3 box cosmological N-body simulation, the `Outer Rim' run. It covers 440 deg^2 of sky area to a redshift of z=3 and is complete to a magnitude depth of 28 in the r-band. Each galaxy is characterized by a multitude of properties including stellar mass, morphology, spectral energy distributions, broadband filter magnitudes, host halo information and weak lensing shear. The size and complexity of cosmoDC2 requires an efficient catalog generation methodology; our approach is based on a new hybrid technique that combines data-driven empirical approaches with semi-analytic galaxy modeling. A wide range of observation-based validation tests has been implemented to ensure that cosmoDC2 enables the science goals of the planned LSST DESC DC2 analyses. This paper also represents the official release of the cosmoDC2 data set, including an efficient reader that facilitates interaction with the data.

A challenging requirement posed by next-generation observations is a firm theoretical grasp of the impact of baryons on structure formation. Cosmological hydrodynamic simulations modeling gas physics are vital in this regard. A high degree of modeling flexibility exists in this space making it important to explore a range of methods in order to gauge the accuracy of simulation predictions. We present results from the first cosmological simulation using Conservative Reproducing Kernel Smoothed Particle Hydrodynamics (CRK-SPH). We employ two simulations: one evolved purely under gravity and the other with non-radiative hydrodynamics. Each contains 2x2304^3 cold dark matter plus baryon particles in an 800 Mpc/h box. We compare statistics to previous non-radiative simulations including power spectra, mass functions, baryon fractions, and concentration. We find self-similar radial profiles of gas temperature, entropy, and pressure and show that a simple analytic model recovers these results to better than 40% over two orders of magnitude in mass. We quantify the level of non-thermal pressure support in halos and demonstrate that hydrostatic mass estimates are biased low by 24% (10%) for halos of mass 10^15 (10^13) Msun/h. We compute angular power spectra for the thermal and kinematic Sunyaev-Zel'dovich effects and find good agreement with the low-l Planck measurements. Finally, artificial scattering between particles of unequal mass is shown to have a large impact on the gravity-only run and we highlight the importance of better understanding this issue in hydrodynamic applications. This is the first in a simulation campaign using CRK-SPH with future work including subresolution gas treatments.

Delensing is an increasingly important technique to reverse the gravitational lensing of the cosmic microwave background (CMB) and thus reveal primordial signals the lensing may obscure. We present a first demonstration of delensing on Planck temperature maps using the cosmic infrared background (CIB). Reversing the lensing deflections in Planck CMB temperature maps using a linear combination of the 545 and 857GHz maps as a lensing tracer, we find that the lensing effects in the temperature power spectrum are reduced in a manner consistent with theoretical expectations. In particular, the characteristic sharpening of the acoustic peaks of the temperature power spectrum resulting from successful delensing is detected at a significance of 16$\rm{\sigma}$, with an amplitude of $A_{\rm{delens}} = 1.12 \pm 0.07$ relative to the expected value of unity. This first demonstration on data of CIB delensing, and of delensing techniques in general, is significant because lensing removal will soon be essential for achieving high-precision constraints on inflationary B-mode polarization.

Correlations of galaxy ellipticities with large-scale structure, due to galactic tidal interactions, provide a potentially significant contaminant to measurements of cosmic shear. However, these intrinsic alignments are still poorly understood for galaxies at the redshifts typically used in cosmic shear analyses. For spiral galaxies, it is thought that tidal torquing is significant in determining alignments resulting in zero correlation between the intrinsic ellipticity and the gravitational potential in linear theory. Here, we calculate the leading-order correction to this result in the tidal-torque model from non-linear evolution, using second-order perturbation theory, and relate this to the contamination from intrinsic alignments to the recently-measured cross-correlation between galaxy ellipticities and the CMB lensing potential. On the scales relevant for CMB lensing observations, the squeezed limit of the gravitational bispectrum dominates the correlation. Physically, the large-scale mode that sources CMB lensing modulates the small-scale power and hence the intrinsic ellipticity, due to non-linear evolution. We find that the angular cross-correlation from tidal torquing has a very similar scale dependence as in the linear alignment model, believed to be appropriate for elliptical galaxies. The amplitude of the cross-correlation is predicted to depend strongly on the formation redshift, being smaller for galaxies that formed at higher redshift when the bispectrum of the gravitational potential was smaller. Finally, we make simple forecasts for constraints on intrinsic alignments from the correlation of forthcoming cosmic shear measurements with current CMB lensing measurements. We note that cosmic variance can be significantly reduced in measurements of the difference in the intrinsic alignments for elliptical and spiral galaxies if these can be separated (e.g., using colour).

We measure the cross-correlation between weak lensing of galaxy images and of the cosmic microwave background (CMB). The effects of gravitational lensing on different sources will be correlated if the lensing is caused by the same mass fluctuations. We use galaxy shape measurements from 139 deg$^{2}$ of the Dark Energy Survey (DES) Science Verification data and overlapping CMB lensing from the South Pole Telescope (SPT) and Planck. The DES source galaxies have a median redshift of $z_{\rm med} {\sim} 0.7$, while the CMB lensing kernel is broad and peaks at $z{\sim}2$. The resulting cross-correlation is maximally sensitive to mass fluctuations at $z{\sim}0.44$. Assuming the Planck 2015 best-fit cosmology, the amplitude of the DES$\times$SPT cross-power is found to be $A = 0.88 \pm 0.30$ and that from DES$\times$Planck to be $A = 0.86 \pm 0.39$, where $A=1$ corresponds to the theoretical prediction. These are consistent with the expected signal and correspond to significances of $2.9 \sigma$ and $2.2 \sigma$ respectively. We demonstrate that our results are robust to a number of important systematic effects including the shear measurement method, estimator choice, photometric redshift uncertainty and CMB lensing systematics. Significant intrinsic alignment of galaxy shapes would increase the cross-correlation signal inferred from the data; we calculate a value of $A = 1.08 \pm 0.36$ for DES$\times$SPT when we correct the observations with a simple IA model. With three measurements of this cross-correlation now existing in the literature, there is not yet reliable evidence for any deviation from the expected LCDM level of cross-correlation, given the size of the statistical uncertainties and the significant impact of systematic errors, particularly IAs. We provide forecasts for the expected signal-to-noise of the combination of the five-year DES survey and SPT-3G.

We present the analysis of the gravitational microlensing event OGLE-2013-BLG-0102. The light curve of the event is characterized by a strong short-term anomaly superposed on a smoothly varying lensing curve with a moderate magnification $A_{\rm max}\sim 1.5$. It is found that the event was produced by a binary lens with a mass ratio between the components of $q = 0.13$ and the anomaly was caused by the passage of the source trajectory over a caustic located away from the barycenter of the binary. From the analysis of the effects on the light curve due to the finite size of the source and the parallactic motion of the Earth, the physical parameters of the lens system are determined. The measured masses of the lens components are $M_{1} = 0.096 \pm 0.013~M_{\odot}$ and $M_{2} = 0.012 \pm 0.002~M_{\odot}$, which correspond to near the hydrogen-burning and deuterium-burning mass limits, respectively. The distance to the lens is $3.04 \pm 0.31~{\rm kpc}$ and the projected separation between the lens components is $0.80 \pm 0.08~{\rm AU}$.

We detect a cold, terrestrial planet in a binary-star system using gravitational microlensing. The planet has low mass (2 Earth masses) and lies projected at $a_{\perp,ph}$ ~ 0.8 astronomical units (AU) from its host star, similar to the Earth-Sun distance. However, the planet temperature is much lower, T<60 Kelvin, because the host star is only 0.10--0.15 solar masses and therefore more than 400 times less luminous than the Sun. The host is itself orbiting a slightly more massive companion with projected separation $a_{\perp,ch}=$10--15 AU. Straightforward modification of current microlensing search strategies could increase their sensitivity to planets in binary systems. With more detections, such binary-star/planetary systems could place constraints on models of planet formation and evolution. This detection is consistent with such systems being very common.

We report the discovery of MOA-2013-BLG-220Lb, which has a super-Jupiter mass ratio $q=3.01\pm 0.02\times 10^{-3}$ relative to its host. The proper motion, $\mu=12.5\pm 1\, {\rm mas}\,{\rm yr}^{-1}$, is one of the highest for microlensing planets yet discovered, implying that it will be possible to separately resolve the host within $\sim 7$ years. Two separate lines of evidence imply that the planet and host are in the Galactic disk. The planet could have been detected and characterized purely with follow-up data, which has important implications for microlensing surveys, both current and into the LSST era.

A measurement by microlensing of the planetary mass function of planets with masses ranging from 5M_E to 10M_J and orbital radii from 0.5 to 10 AU was reported recently. A strategy for extending the mass range down to (1-3)M_E is proposed here. This entails monitoring the peaks of a few tens of microlensing events with moderately high magnifications with 1-2m class telescopes. Planets of a few Earth masses are found to produce deviations of ~ 5% to the peaks of microlensing light curves with durations ~ (0.7-3)hr in events with magnification ~ 100 if the projected separation of the planet lies in the annular region (0.85-1.2)r_E. Similar deviations are produced by Earth mass planets in the annular region (0.95-1.05)r_E. It is possible that sub-Earths could be detected very close to the Einstein ring if they are sufficiently abundant, and also planetary systems with more than one low mass planet.