Day 1631

Wednesday.  Thursday.

1909.07976
Cosmological simulations of galaxy formation
Vogelsberger, et al

Over the last decades, cosmological simulations of galaxy formation have been instrumental for advancing our understanding of structure and galaxy formation in the Universe. These simulations follow the non-linear evolution of galaxies modeling a variety of physical processes over an enormous range of scales. A better understanding of the physics relevant for shaping galaxies, improved numerical methods, and increased computing power have led to simulations that can reproduce a large number of observed galaxy properties. Modern simulations model dark matter, dark energy, and ordinary matter in an expanding space-time starting from well-defined initial conditions. The modeling of ordinary matter is most challenging due to the large array of physical processes affecting this matter component. Cosmological simulations have also proven useful to study alternative cosmological models and their impact on the galaxy population. This review presents a concise overview of the methodology of cosmological simulations of galaxy formation and their different applications.

1909.07986

Cosmic dissonance: new physics are systematics behind a short sound horizon?

Arendse, et al

Persistent tension between low-redshift and the Cosmic Microwave Background (CMB) measurements, in terms of the sound horizon and the Hubble-Lema\^itre constant, suggests new physics beyond the Standard Model, departures from concordance cosmology, or residual systematics. Assessing the likelihood of new physics devised to resolve the tension requires thorough consistency tests of several independent and high-precision distance calibrations. We examine recent updated distance calibrations from Cepheids, gravitational lensing time-delay observations, and the Tip of the Red Giant Branch. Through selected cosmographic methods for combining observations of the Baryon Acoustic Oscillations (BAO), Type Ia supernovae and local distance calibrators, we obtain robust measurements of the sound horizon, independent of cosmological models and the CMB. For our different models and probes, the sound horizon scale varies between $r_{\rm s}=(135\pm3)$ and $r_{\rm s}=(140\pm3)$~Mpc. The combined tension in the sound horizon and the Hubble-Lema\^itre constant ranges between $3\sigma$ and $5\sigma$, depending on whether distance ladder calibrations are used, and is independent of possible changes to the late-Universe expansion history. Early-Universe extensions with $\textrm{N}_{\textrm{eff}}=3.27\pm0.15$ are allowed by the data, unless a Cepheid distance calibration is included. Some late-dark-energy models, as constrained by the relative supernova distance moduli, do not resolve the current tension. Results from time-delay lenses are consistent with those from distance-ladder calibrations. With upcoming lens samples, the tension in the inferred sound horizon will be completely independent from local calibrations. New proposals to resolve the tension should examine CMB and late-Universe constraints separately, and the inference on both $H_0$ and $r_s.$

Day 1630

Tuesday.

Nature
Water vapor in the atmosphere of the habitable-zone eight-Earth-mass planet K2-18 b
Tsiaras, et al

In the past decade, observations from space and the ground have found water to be the most abundant molecular species, after hydrogen, in the atmospheres of hot, gaseous extrasolar planets. Being the main molecular carrier of oxygen, water is a tracer of the origin and the evolution mechanisms of planets. For temperate, terrestrial planets, the presence of water is of great importance as an indicator of habitable conditions. Being small and relatively cold, these planets and their atmospheres are the most challenging to observe, and therefore no atmospheric spectral signatures have so far been detected. Super-Earths—planets lighter than ten Earth masses—around later-type stars may provide our first opportunity to study spectroscopically the characteristics of such planets, as they are best suited for transit observations. Here, we report the detection of a spectroscopic signature of water in the atmosphere of K2-18 b—a planet of eight Earth masses in the habitable zone of an M dwarf—with high statistical confidence (Atmospheric Detectability Index = 5.0, ~3.6σ (refs.8,9)). In addition, the derived mean molecular weight suggests an atmosphere still containing some hydrogen. The observations were recorded with the Hubble Space Telescope/Wide Field Camera 3 and analysed with our dedicated, publicly available, algorithms. Although the suitability of M dwarfs to host habitable worlds is still under discussion, K2-18 b offers an unprecedented opportunity to gain insight into the composition and climate of habitable-zone planets.

1909.06748
Planet formation around super massive black holes i the active galactic nuclei
Wada, et al

As a natural consequence of the elementary processes of dust growth, we discovered that a new class of planets can be formed around supermassive black holes (SMBHs). We investigated a growth path from sub-micron sized icy dust monomers to Earth-sized bodies outside the "snow line'', located several parsecs from SMBHs in low luminosity active galactic nuclei (AGNs). In contrast to protoplanetary disks, the "radial drift barrier'' does not prevent the formation of planetesimals. In the early phase of the evolution, low collision velocity between dust particles promotes sticking; therefore, the internal density of the dust aggregates decreases with growth. When the porous aggregate's size reaches 0.1--1 cm, the collisional compression becomes effective, and the decrease in internal density stops. Once 10--100 m sized aggregates are formed, they are decoupled from gas turbulence, and the aggregate layer becomes gravitationally unstable, leading to the formation of planets by the fragmentation of the layer, with ten times the mass of the earth. The growth time scale depends on the turbulent strength of the circumnuclear disk and the black hole mass $M_{BH}$, and it is comparable to the AGN's lifetime ($\sim 10^8$ yr) for low mass ($M_{BH} \sim 10^6 M_\odot$) SMBHs.

Day 1629

Monday.

1909.05898
The CFHT Large Area U-band Deep Survey (CLAUDS)
Sawicki, et al

The CFHT Large Area U-band Deep Survey (CLAUDS) uses data taken with the MegaCam mosaic imager on CFHT to produce images of 18.60 deg2 with median seeing of FWHM=0.92 arcsec and to a median depth of U = 27.1 AB (5 sigma in 2 arcsec apertures), with selected areas that total 1.36 deg2 reaching a median depth of U=27.7 AB. These are the deepest U-band images assembled to date over this large an area. These data are located in four fields also imaged to comparably faint levels in grizy and several narrowband filters as part of the Hyper Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP). These CFHT and Subaru datasets will remain unmatched in their combination of area and depth until the advent of the Large Synoptic Survey Telescope (LSST). This paper provides an overview of the scientific motivation for CLAUDS and gives details of the observing strategy, observations, data reduction, and data merging with the HSC-SSP. Three early applications of these deep data are used to illustrate the potential of the dataset: deep U-band galaxy number counts, z~3 Lyman break galaxy (LBG) selection, and photometric redshifts improved by adding CLAUDS U to the Subaru HSC grizy photometry.

1909.05940

The Urban Observatory: A multi-model imaging platform for the study of dynamics in complex urban systems

Dobler, et al

We describe an "Urban Observatory" facility designed for the study of complex urban systems via persistent, synoptic, and granular imaging of dynamical processes in cities. An initial deployment of the facility has been demonstrated in New York City and consists of a suite of imaging systems - both broadband and hyperspectral - sensitive to wavelengths from the visible (~400 nm) to the infrared (~13 micron) operating at cadences of ~0.01 - 30 Hz (characteristically ~0.1 Hz). Much like an astronomical survey, the facility generates a large imaging catalog from which we have extracted observables (e.g., time-dependent brightnesses, spectra, temperatures, chemical species, etc.), collecting them in a parallel source catalog. We have demonstrated that, in addition to the urban science of cities as systems, these data are applicable to a myriad of domain-specific scientific inquiries related to urban functioning including energy consumption and end use, environmental impacts of cities, and patterns of life and public health. We show that an Urban Observatory facility of this type has the potential to improve both a city's operations and the quality of life of its inhabitants.

1909.05947

NASA's Meteoroid Engineering Model (MEM) 3 and its ability to replicate spacecraft impact rates

Moorhead, et al

Meteoroids pose one of the largest risks to spacecraft outside of low Earth orbit. In order to correctly predict the rate at which meteoroids impact and damage spacecraft, environment models must describe the mass, directionality, velocity, and density distributions of meteoroids. NASA's Meteoroid Engineering Model (MEM) is one such model; MEM 3 is an updated version of the code that better captures the correlation between directionality and velocity and incorporates a bulk density distribution. This paper describes MEM 3 and compares its predictions with the rate of large particle impacts seen on the Long Duration Exposure Facility (LDEF) and the Pegasus II and III satellites.

1909.06190

How to measure galaxy-galaxy-galaxy lensing with higher precision and accuracy

Linke, Simon, Schneider, Hilbert

Galaxy-galaxy-galaxy lensing (G3L) is a powerful tool for constraining the three-point correlation between the galaxy and the matter field and thereby models of galaxy evolution. We propose three improvements to current measurements of G3L, designed to improve the precision and the accuracy by using the galaxies' redshifts and removing biases of the estimator. We further show how to account for lens galaxy magnification by the cosmic large-scale structure and how to convert the G3L signal from angular to physical scales. The improvements are tested on simple mock data and simulated data based on the Millennium Run with an implemented semi-analytic model of galaxies. Our improvements increase the signal-to-noise ratio by on average $35~\%$ at angular scales between $0.1'$ and $10'$ and physical scales between $0.02$ and $2 \, h^{-1}\,\textrm{Mpc}$. They also remove the bias of the G3L estimator at angular scales below $1'$, which was originally up to $40\,\%$. The signal due to lens magnification is approximately $10\,\%$ of the total signal.

Day 1628

Thursday.  Friday.

1909.04663
The ultra-diffuse dwarf galaxies NGC 1052-DR2 and 1052-DF4 are in conflict with standard cosmology
Haslbauer, et al

Recently van Dokkum et al. (2018b) reported that the galaxy NGC 1052-DF2 (DF2) lacks dark matter if located at $20$ Mpc from Earth. In contrast, DF2 is a dark-matter-dominated dwarf galaxy with a normal globular cluster population if it has a much shorter distance near $10$ Mpc. However, DF2 then has a high peculiar velocity wrt. the cosmic microwave background of $886$ $\rm{km\,s^{-1}}$, which differs from that of the Local Group (LG) velocity vector by $1298$ $\rm{km\,s^{-1}}$ with an angle of $117 \, ^{\circ}$. Taking into account the dynamical $M/L$ ratio, the stellar mass, half-light radius, peculiar velocity, motion relative to the LG, and the luminosities of the globular clusters, we show that the probability of finding DF2-like galaxies in the lambda cold dark matter ($\Lambda$CDM) TNG100-1 simulation is at most $1.0\times10^{-4}$ at $11.5$ Mpc and is $4.8\times10^{-7}$ at $20.0$ Mpc. At $11.5$ Mpc, the peculiar velocity is in significant tension with the TNG100-1, TNG300-1, and Millennium simulations, but occurs naturally in a Milgromian cosmology. At $20.0$ Mpc, the unusual globular cluster population would challenge any cosmological model. Estimating that precise measurements of the internal velocity dispersion, stellar mass, and distance exist for $100$ galaxies, DF2 is in $2.6\sigma$ ($11.5$ Mpc) and $4.1\sigma$ ($20.0$ Mpc) tension with standard cosmology. Adopting the former distance for DF2 and assuming that NGC 1052-DF4 is at $20.0$ Mpc, the existence of both is in tension at $\geq4.8\sigma$ with the $\Lambda$CDM model. If both galaxies are at $20.0$ Mpc the $\Lambda$CDM cosmology has to be rejected by $\geq5.8\sigma$.

1909.04664

EDGE: the origin of scatter in ultra-faint dwarf stellar masses and surface brightnesses

Rey, et al

We demonstrate how the least luminous galaxies in the Universe, ultra-faint dwarf galaxies, are sensitive to their dynamical mass at the time of cosmic reionization. We select a low-mass ($\sim \text{1.5} \times 10^{9} \, \text{M}_{\odot}$) dark matter halo from a cosmological volume, and perform zoom hydrodynamical simulations with multiple alternative histories using "genetically modified" initial conditions. Earlier forming ultra-faints have higher stellar mass today, due to a longer period of star formation before their quenching by reionization. Our histories all converge to the same final dynamical mass, demonstrating the existence of extended scatter ($\geq$ 1 dex) in stellar masses at fixed halo mass due to the diversity of possible histories. One of our variants builds less than 2 \% of its final dynamical mass before reionization, rapidly quenching in-situ star formation. The bulk of its final stellar mass is later grown by dry mergers, depositing stars in the galaxy's outskirts and hence expanding its effective radius. This mechanism constitutes a new formation scenario for highly diffuse ($\text{r}_{1 /2} \sim 820 \, \text{pc}$, $\sim 32 \, \text{mag arcsec}^2$), metal-poor ($\big[ \mathrm{Fe}\, / \mathrm{H} \big]= -2.9$), ultra-faint ($\mathcal{M}_V= -5.7$) dwarf galaxies within the reach of next-generation low surface brightness surveys.

1909.04670

A new sample of (wandering) massive black holes in dwarf galaxies from high-resolution radio observations

Reines, et al

We present a sample of nearby dwarf galaxies with radio-selected accreting massive black holes (BHs), the majority of which are non-nuclear. We observed 111 galaxies using sensitive, high-resolution observations from the Karl G. Jansky Very Large Array (VLA) in its most extended A-configuration at X-band (~8-12 GHz), yielding a typical angular resolution of ~0.25" and rms noise of ~15 uJy. Our targets were selected by cross matching galaxies with stellar masses M_stellar < 3 x 10^9 M_sun and redshifts z<0.055 in the NASA-Sloan Atlas with the VLA Faint Images of the Radio Sky at Twenty centimeters (FIRST) Survey. With our new high-resolution VLA observations, we detect compact radio sources towards 39 galaxies and carefully evaluate possible origins for the radio emission including thermal HII regions, supernova remnants, younger radio supernovae, background interlopers, and AGNs in the target galaxies. We find that 13 dwarf galaxies almost certainly host active massive BHs despite the fact that only one object was previously identified as having optical signatures of an AGN. We also identify a candidate dual radio AGN in a more massive galaxy system. The majority of the radio-detected BHs are offset from the center of the host galaxies with some systems showing signs of interactions/mergers. Our results indicate that massive BHs need not always live in the nuclei of dwarf galaxies, confirming predictions from simulations. Moreover, searches attempting to constrain BH seed formation using observations of dwarf galaxies need to account for such a population of "wandering" BHs.

1909.04672

The MUSE Ultra Deep Field (MUDF). II. Survey design and the gaseous properties of galaxy groups at 0.5<z<1.5

Fossati, et al

We present the goals, design, and first results of the MUSE Ultra Deep Field (MUDF) survey, a large programme using the Multi Unit Spectroscopic Explorer (MUSE) instrument at the ESO Very Large Telescope. The MUDF survey is collecting ~ 150 hours on-source of integral field optical spectroscopy in a 1.5 x 1.2 square arcmin region which hosts several astrophysical structures along the line of sight, including two bright z ~ 3.2 quasars with close separation (~ 500 kpc). Following the description of the data reduction procedures, we present the analysis of the galaxy environment and gaseous properties of seven groups detected at redshifts 0.5 < z < 1.5, spanning a large dynamic range in halo mass, log(Mh/Msun) ~ 11 - 13.5. For four of the groups, we find associated MgII absorbers tracing cool gas in high-resolution spectroscopy of the two quasars, including one case of correlated absorption in both sightlines at distance ~ 480 kpc. The absorption strength associated with the groups is higher than what has been reported for more isolated galaxies of comparable mass and impact parameters. We do not find evidence for widespread cool gas giving rise to strong absorption within these groups. Combining these results with the distribution of neutral and ionised gas seen in emission in lower-redshift groups, we conclude that gravitational interactions in the group environment strip gas from the galaxy haloes into the intragroup medium, boosting the cross section of cool gas and leading to the high fraction of strong MgII absorbers that we detect.

1909.04679

The biggest splash

Belokurov, et al

Using a large sample of bright nearby stars with accurate Gaia Data Release 2 astrometry and auxiliary spectroscopy we map out the properties of the principle Galactic components such as the "thin" and "thick" discs and the halo. We show that in the Solar neighborhood, there exists a large population of metal-rich ([Fe/H]>-0.7) stars on highly eccentric orbits. By studying the evolution of elemental abundances, kinematics and stellar ages in the plane of azimuthal velocity v_phi and metallicity [Fe/H], we demonstrate that this metal-rich halo-like component, which we dub the Splash, is linked to the alpha-rich (or "thick") disc. Splash stars have little to no angular momentum and many are on retrograde orbits. They are predominantly old, but not as old as the stars deposited into the Milky Way in the last major merger. We argue, in agreement with several recent studies, that the Splash stars may have been born in the Milky Way's proto-disc prior to the massive ancient accretion event which drastically altered their orbits. We can not, however, rule out other (alternative) formation channels. Taking advantage of the causal connection between the merger and the Splash, we put constraints of the epoch of the last massive accretion event to have finished 9.5 Gyr ago. The link between the local metal-rich and metal-poor retrograde stars is confirmed using a large suite of cutting-edge numerical simulations of the Milky Way's formation.

1909.04765

The separation distribution of ultra-wide binaries across galactic populations

Tian, et al

We present an extensive sample of ultra-wide binary stars in the solar neighborhood, focusing on separations of $0.01<s/{pc}<1$. Using data from Gaia DR2, we define kinematic sub-populations via the systems' tangential velocities, i.e., disk-like ($v_{\perp, tot}<40$km/s), intermediate ($v_{\perp,tot}=(40-85)$km/s), and halo-like ($v_{\perp,tot}>85$km/s) samples, presuming that these velocity cuts represent a rough ordering in the binaries' age and metallicity. Through stringent cuts on astrometric precision, we can obtain pure binary samples with thousands of binaries in each sample. For all three populations, the distribution of binary separations extends smoothly to 1pc, displaying neither strong truncation nor bimodality. Fitting a smoothly-broken power law for the separation distribution, we find that its slope at separations $s=10^{2.5-4}$AU is the same for all sub-populations, $p(s)\propto s^{\gamma}$, with $\gamma\sim-1.54$. However, the logarithmic slope of $p(s)$ steepens at $s \gtrsim 10^4$AU. We find some evidences that the degree of steepening increases with the binaries' age, with a slope-change of only $\Delta\gamma\sim0.5$ for disk-like stars, but $\Delta\gamma>1$ for halo-like stars. This trend is contrary to what might be expected if steepening at wide separations were due to gravitational perturbations by molecular clouds or stars, which would preferentially disrupt disk binaries. If we were to interpret steepening at $s\gtrsim 10^4$AU as a consequence of disruption by massive halo objects, we would have to invoke a MACHO population inconsistent with other constraints. As a more plausible alternative, we propose a simple model to predict the separation distribution of wide binaries formed in dissolving star clusters. This model generically predicts $\gamma\sim-1.5$ as observed, with steepening at larger separations due to the finite size of binaries' birth clusters.

1909.05277

Cosmological parameters from the BOSS galaxy power spectrum

Ivanov, et al

We present cosmological parameter measurements from the publicly available Baryon Oscillation Spectroscopic Survey (BOSS) data on anisotropic galaxy clustering in Fourier space. Compared to previous studies, our analysis has two main novel features. First, we use a complete perturbation theory model that properly takes into account the non-linear effects of dark matter clustering, short-scale physics, galaxy bias, redshift-space distortions, and large-scale bulk flows. Second, we employ a Markov-Chain Monte-Carlo technique and consistently reevaluate the full power spectrum likelihood as we scan over different cosmologies. Assuming a minimal $\Lambda$CDM cosmology with massive neutrinos, fixing the primordial power spectrum tilt, and imposing the big bang nucleosynthesis (BBN) prior on the physical baryon density $\omega_b$, we find the following late-Universe parameters: Hubble constant $H_0=(67.89\pm 1.06)$ km$\,$s$^{-1}$Mpc$^{-1}$, matter density fraction $\Omega_m=0.295\pm 0.010$, and the mass fluctuation amplitude $\sigma_8=0.721\pm 0.043$. These parameters were measured directly from the BOSS data and independently of the Planck cosmic microwave background observations. Our constraints remain unchanged if we replace the BBN prior on $\omega_b$ by a tighter Planck constraint. Finally, we discuss the information content of the BOSS power spectrum and show that it is dominated by the location of the baryon acoustic oscillations and the power spectrum shape. We argue that the contribution of the Alcock-Paczynski effect is marginal in $\Lambda$CDM, but becomes important for non-minimal cosmological models.

1909.05852

A gravitational lensing detection of filamentary structures connecting luminous red galaxies

Xia, et al

We present a weak lensing detection of filamentary structures in the cosmic web, combining data from the Kilo-Degree Survey, the Red Cluster Sequence Lensing Survey and the Canada-France-Hawaii Telescope Lensing Survey. The line connecting luminous red galaxies with a separation of $3 - 5\, h^{-1}\text{Mpc}$ is chosen as a proxy for the location of filaments. We measure the average weak lensing shear around $\sim$11,000 candidate filaments selected in this way from the Sloan Digital Sky Survey. After nulling the shear induced by the dark matter haloes around each galaxy, we report a $3.3\,\sigma$ detection of an anisotropic shear signal from the matter that connects them. Adopting a filament density profile, motivated from $N$-body simulations, the average density at the centre of these filamentary structures is found to be $15 \pm 4$ times the critical density.

Day 1627

Tuesday.  Wednesday.

1909.03107
Cosmological evidence modelling: a new simulation-based approach to constrain cosmology on non-linear scales
Lange, et al

Extracting accurate cosmological information from galaxy-galaxy and galaxy-matter correlation functions on non-linear scales ($\lesssim 10 h^{-1} \mathrm{Mpc}$) requires cosmological simulations. Additionally, one has to marginalise over several nuisance parameters of the galaxy-halo connection. However, the computational cost of such simulations prohibits naive implementations of stochastic posterior sampling methods like Markov chain Monte Carlo (MCMC) that would require of order $\mathcal{O}(10^6)$ samples in cosmological parameter space. Several groups have proposed surrogate models as a solution: a so-called emulator is trained to reproduce observables for a limited number of realisations in parameter space. Afterwards, this emulator is used as a surrogate model in an MCMC analysis. Here, we demonstrate a different method called Cosmological Evidence Modelling (CEM). First, for each simulation, we calculate the Bayesian evidence marginalised over the galaxy-halo connection by repeatedly populating the simulation with galaxies. We show that this Bayesian evidence is directly related to the posterior probability of cosmological parameters. Finally, we build a physically motivated model for how the evidence depends on cosmological parameters as sampled by the simulations. We demonstrate the feasibility of CEM by using simulations from the Aemulus simulation suite and forecasting cosmological constraints from BOSS CMASS measurements of redshift-space distortions. Our analysis includes an exploration of how galaxy assembly bias affects cosmological inference. Overall, CEM has several potential advantages over the more common approach of emulating summary statistics, including the ability to easily marginalise over highly complex models of the galaxy-halo connection and greater accuracy, thereby reducing the number of simulations required.

1909.03138

The impact of type Ia supernovae in quiescent galaxies: I. Formation of the multiphase interstellar medium

Li, et al

A cool phase of the interstellar medium has been observed in many giant elliptical galaxies, but its origin remains unclear. We propose that uneven heating from Type Ia supernovae (SNe Ia), together with radiative cooling, can lead to the formation of the cool phase. The basic idea is that since SNe Ia explode randomly, gas parcels which are not directly heated by SN shocks will cool, forming multiphase gas. We run a series of idealized high-resolution numerical simulations, and find that cool gas develops even when the overall SNe heating rate $H$ exceeds the cooling rate $C$ by a factor as large as 1.4. We also find that the time for multiphase gas development depends on the gas temperature. When the medium has a temperature $T = 3\times 10^6$ K, the cool phase forms within one cooling time \tc; however, the cool phase formation is delayed to a few times \tc\ for higher temperatures. The main reason for the delay is turbulent mixing. Cool gas formed this way would naturally have a metallicity lower than that of the hot medium. For constant $H/C$, there is more turbulent mixing for higher temperature gas. We note that this mechanism of producing cool gas cannot be captured in cosmological simulations, which usually fail to resolve individual SN remnants.

1909.03996

The formation of the Martian moons

Rosenblatt, et al

Almost all the planets of our solar system have moons. Each planetary system has however unique characteristics. The Martian system has not one single big moon like the Earth, not tens of moons of various sizes like for the giant planets, but two small moons: Phobos and Deimos. How did form such a system? This question is still being investigated on the basis of the Earth-based and space-borne observations of the Martian moons and of the more modern theories proposed to account for the formation of other moon systems. The most recent scenario of formation of the Martian moons relies on a giant impact occurring at early Mars history and having also formed the so-called hemispheric crustal dichotomy. This scenario accounts for the current orbits of both moons unlike the scenario of capture of small size asteroids. It also predicts a composition of disk material as a mixture of Mars and impactor materials that is in agreement with remote sensing observations of both moon surfaces, which suggests a composition different from Mars. The composition of the Martian moons is however unclear, given the ambiguity on the interpretation of the remote sensing observations. The study of the formation of the Martian moon system has improved our understanding of moon formation of terrestrial planets: The giant collision scenario can have various outcomes and not only a big moon as for the Earth. This scenario finds a natural place in our current vision of the early solar system when conditions were favorable for giant collisions to occur. The next step in exploration of Martian moon is a sample return mission to test the giant collision scenario for their origin, and to provide tests of models of early solar system dynamics since Mars may retain material exchanged between the inner and outer solar system.

Day 1626

Monday.

1909.02615
Game of cones: a nulling strategy for modelling lensing convergence in cones with large deviation theory
Barthelemy, et al

The distribution of the cosmic convergence field is modelled from first principles using a large deviation principle. The geometry of the past light cone is accounted for by constructing the total weak-lensing signal from contributions of the matter density in thin disk slices. The prediction of this model is successfully tested against numerical simulation with ray tracing, and found to be accurate within at least 5 per cent in the tails at redshift 1 and opening angle of 10 arcmin and even more so with increasing source redshift and opening angle. An accurate analytical approximation to the theory is also provided for practical implementation. The lensing kernel that mixes physical scales along the line-of-sight tends to reduce the domain of validity of this theoretical approach compared to the three dimensional case of cosmic densities in spherical cells. This effect is shown to be avoidable if a nulling procedure is implemented in order to localise the lensing line-of-sight integrations in a tomographic analysis. Accuracy in the tails is thus achieved within a percent for source redshifts between 0.5 and 1.5 and an opening angle of 10 arcmin. Applications to future weak-lensing surveys like Euclid and the specific issue of shape noise are discussed.

1909.02908

Study on the 3D anisotropic propagation of Galactic cosmic rays

Liu, et al

Conventional cosmic-ray propagation models usually assume an isotropic diffusion coefficient to account for the random deflection of cosmic rays by the turbulent interstellar magnetic field. Such a picture is very successful in explaining a lot of observational phenomena related to the propagation of galactic cosmic rays, such as broken power-law energy spectra, secondary-to-primary ratios, etc. However, the isotropic diffusion presupposition is facing severe challenges from observations. In particular, recent observations on the large-scale anisotropy of TeV cosmic rays show that the dipole direction differs from the prediction of the conventional model. One possible reason is that the large-scale regular magnetic field, which leads to an anisotropic diffusion of cosmic rays, has not been included in the model provided by the public numerical packages. In this work, we propose a pseudo source method to solve the 3-dimensional anisotropic transport equation. Based on a realistic configuration of Galactic magnetic field, both proton energy spectrum and B/C ratio are reproduced with proper transport parameters.

1909.02928

Oceanographic constrations on exoplanet life

Olson, et al

Liquid water oceans are at the center of our search for life on exoplanets because water is a strict requirement for life as we know it. However, oceans are dynamic habitats--and some oceans may be better hosts for life than others. In Earth's ocean, circulation transports essential nutrients such as P and controls the distribution and productivity of life. Of particular importance is upwelling due to wind-driven divergence in surface layers, which returns essential nutrients that tend to accumulate at depth via the biological pump back to the sunlit regions where life is concentrated. Atmospheric buildup of biogenic gases like CH$_4$ that are produced in the subsurface ocean also depends on upwelling. Ocean circulation is likely to impose constraints on the activity and atmospheric expression of life in exo-oceans as well, but we lack an understanding of how ocean dynamics may differ on other planets. We address this issue by using ROCKE-3D, a fully coupled ocean-atmosphere GCM, to investigate ocean dynamics on a diversity of habitable planets. Our results suggest that planets that rotate slower and have higher surface pressure than Earth may be attractive targets for remote life detection because upwelling is enhanced on these planets. Seasonal deepening of the mixed layer on high obliquity planets may also enhance nutrient return from depth into the surface mixed layer. Efficient nutrient recycling favors greater biological activity, more biosignature production, and more detectable life. We also find that an absence of continents may be problematic for sustaining a globally active, remotely detectable ocean biosphere due to limited nutrient fluxes from coastal upwelling and continental weathering. Our results demonstrate the importance of considering oceanographic phenomena for exoplanet life detection and motivate further contributions to the emerging field of exo-oceanography.

Day 1625

Friday.

1909.02007
Resolving the metallicity distribution of the stellar halo with the H3 survey
Conroy, et al

The Galactic stellar halo is predicted to have formed at least partially from the tidal disruption of accreted dwarf galaxies. This assembly history should be detectable in the orbital and chemical properties of stars. The H3 Survey is obtaining spectra for 200,000 stars, and, when combined with Gaia data, is providing detailed orbital and chemical properties of Galactic halo stars. Unlike previous surveys of the halo, the H3 target selection is based solely on magnitude and Gaia parallax; the survey therefore provides a nearly unbiased view of the entire stellar halo at high latitudes. In this paper we present the distribution of stellar metallicities as a function of Galactocentric distance and orbital properties for a sample of 4232 kinematically-selected halo giants to 100 kpc. The stellar halo is relatively metal-rich, [Fe/H]=-1.2, and there is no discernable metallicity gradient over the range $6<R_{\rm gal}<100$ kpc. However, the halo metallicity distribution is highly structured including distinct metal-rich and metal-poor components at $R_{\rm gal}<10$ kpc and $R_{\rm gal}>30$ kpc, respectively. Metal-poor stars with [Fe/H]$<-2$ are a small population of the halo at all distances and orbital categories. We associate the "in-situ" stellar halo with stars displaying thick-disk chemistry on halo-like orbits; such stars are confined to $|z|<10$ kpc. The majority of the stellar halo is resolved into discrete features in orbital-chemical space, suggesting that the bulk of the stellar halo formed from the accretion and tidal disruption of dwarf galaxies. (ABRIDGED)

1909.02561

On the road to percent accuracy III: non-linear reaction of the matter power spectrum to massive neutrinos

Cataneo, et al

We analytically model the non-linear effects induced by massive neutrinos on the total matter power spectrum using the halo model reaction framework of Cataneo et al. 2019. In this approach the halo model is used to determine the relative change to the matter power spectrum caused by new physics beyond the concordance cosmology. Using standard fitting functions for the halo abundance and the halo mass-concentration relation, the total matter power spectrum in the presence of massive neutrinos is predicted to percent-level accuracy, out to $k=10 \, h \, {\rm Mpc}^{-1}$. We find that refining the prescriptions for the halo properties using $N$-body simulations improves the recovered accuracy to better than 1%. This paper serves as another demonstration for how the halo model reaction framework, in combination with a single suite of standard $\Lambda$CDM simulations, can recover percent-level accurate predictions for beyond-$\Lambda$CDM matter power spectra, well into the non-linear regime.

Day 1624

Thursday.

1909.01344
Galaxy disc scaling realtions: a tight linear galaxy -- halo connection challenges abundance matching
Posti, et al

In $\Lambda$CDM cosmology, to first order, galaxies form out of the cooling of baryons within the virial radius of their dark matter halo. The fractions of mass and angular momentum retained in the baryonic and stellar components of disc galaxies put strong constraints on our understanding of galaxy formation. In this work, we derive the fraction of angular momentum retained in the stellar component of spirals, $f_j$, the global star formation efficiency $f_M$, and the ratio of the asymptotic circular velocity ($V_{\rm flat}$) to the virial velocity $f_V$, and their scatter, by fitting simultaneously the observed stellar mass-velocity (Tully-Fisher), size-mass, and mass-angular momentum (Fall) relations. We compare the goodness of fit of three models: (i) where the logarithm of $f_j$, $f_M$, and $f_V$ vary linearly with the logarithm of the observable $V_{\rm flat}$; (ii) where these values vary as a double power law; and (iii) where these values also vary as a double power law but with a prior imposed on $f_M$ such that it follows the expectations from widely used abundance matching models. We conclude that the scatter in these fractions is particularly small ($\sim 0.07$ dex) and that the linear model is by far statistically preferred to that with abundance matching priors. This indicates that the fundamental galaxy formation parameters are small-scatter single-slope monotonic functions of mass, instead of being complicated non-monotonic functions. This incidentally confirms that the most massive spiral galaxies should have turned nearly all the baryons associated with their haloes into stars. We call this the failed feedback problem.

1909.01405

Studying black holes on horizon scales with space-VLBI

Haworth, et al

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.

1909.01478

Trans-Neptunian objects found in the first four years of the Dark Energy Survey

Bernardinelli, Berstein, Sako, et al

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.

1909.01565

How do bound star clusters form?

Krumholz, McKee

Gravitationally-bound clusters that survive gas removal represent an unusual mode of star formation in the Milky Way and similar spiral galaxies. While forming, they can be distinguished observationally from unbound star formation by their high densities, virialised velocity structures, and star formation histories that accelerate toward the present, but extend multiple free-fall times into the past. In this paper we examine several proposed scenarios for how such structures might form and evolve, and carry out a Bayesian analysis to test these models against observed distributions of protostellar age, counts of young stellar objects relative to gas, and the overall star formation rate of the Milky Way. We show that models in which the acceleration of star formation is due either to a large-scale collapse or a time-dependent increase in star formation efficiency are unable to satisfy the combined set of observational constraints. In contrast, models in which clusters form in a "conveyor belt" mode where gas accretion and star formation occur simultaneously, but the star formation rate per free-fall time is low, can match the observations.

1909.01592

A space mission to map the entire observable universe using the CMB as a backlight

Basu, et al

This Science White Paper, prepared in response to the ESA Voyage 2050 call for long-term mission planning, aims to describe the various science possibilities that can be realized with an L-class space observatory that is dedicated to the study of the interactions of cosmic microwave background (CMB) photons with the cosmic web. Our aim is specifically to use the CMB as a backlight -- and survey the gas, total mass, and stellar content of the entire observable Universe by means of analyzing the spatial and spectral distortions imprinted on it. These distortions result from two major processes that impact on CMB photons: scattering by electrons (Sunyaev-Zeldovich effect in diverse forms, Rayleigh scattering, resonant scattering) and deflection by gravitational potential (lensing effect). Even though the list of topics collected in this White Paper is not exhaustive, it helps to illustrate the exceptional diversity of major scientific questions that can be addressed by a space mission that will reach an angular resolution of 1.5 arcmin (goal 1 arcmin), have an average sensitivity better than 1 uK-arcmin, and span the microwave frequency range from roughly 50 GHz to 1 THz. The current paper also highlights the synergy of our BACKLIGHT mission concept with several upcoming and proposed ground-based CMB experiments.

1909.01814

Cooperative multi-spacecraft observation of incoming space threats

Nallapu, Thangavelautham

Earth is constantly being bombarded with material from space. Most of the natural material end up being dust grains that litter the surface of Earth, but larger bodies are known to impact every few decades. The most recent large impact was Chelyabinsk which set off a 500-kiloton explosion which was 40 times that of the Hiroshima nuclear explosion. Apart from meteors, there is a growing threat of space assets deorbiting. With these impending space threats, it is critical to have a constellation of satellites to autonomously lookout for meteors and reentering space debris. By using multiple spacecraft, it is possible to perform multipoint observation of the event. Through multipoint observation, it is possible to triangulate the location of the observed event. The detection, tracking, and analysis of these objects all need to be performed autonomously. Our previous work focused on developing several vision algorithms including blob-detection, feature detection, and neural network-based image segment classification. For this multipoint observation to occur, it requires multiple spacecraft to coordinate their actions particularly fixating on the space observation target. Furthermore, communication and coordination are needed for bringing new satellites into observation view and removing other satellites that have lost their view. In this paper, we analyze state-of-the-art observation technology for small satellites and perform detailed design of its implementation. Through this study, we estimate the error estimates on position, velocity, and acceleration. We presume use of low to mid-tier cameras for the spacecraft.

1909.01929

Recognition and classification of the cosmic-ray events in images captured by CMOS/CCD cameras

Niedzwiecki, et al

Muons and other ionizing radiation produced by cosmic rays and radiative decays affect CMOS/CCD sensor. When particles colliding with sensors atoms cause specific kind of noise on images recorded by cameras. We present a concept and preliminary implementation of method for recognizing those events and algorithms for image processing and their classification by machine learning. Our method consists of analyzing the shape of traces present in images recorded by a camera sensor and metadata related to an image like camera model, GPS location of camera, vertical and horizontal orientation of a camera sensor, timestamp of image acquisition, and other events recognized near-by sensors. The so created feature vectors are classified as either a muon-like event, an electron-like event or the other event, possibly noise. For muon-like events our method estimates azimuth of a muon track. Source of the data is database of CREDO (Cosmic-Ray Extremely Distributed Observatory) project and ESO (European Southern Observatory) archives. The telescope dark frames from ESO are analysed. CREDO project collected so far over 2 millions images of events from many kinds of cameralike: smartphones camera, laptop webcams and Internet of Things cameras localised around the globe.

Day 1623

Tuesday.  Wednesday.

1909.00003
Properties of the circumgalactic medium in cosmic ray-dominated galaxy halos
Ji, et al

We investigate the impact of cosmic rays (CRs) on the circumgalactic medium (CGM) in FIRE-2 simulations, for ultra-faint dwarf through Milky Way (MW)-mass halos hosting star-forming (SF) galaxies. Our CR treatment includes injection by supernovae, anisotropic streaming and diffusion along magnetic field lines, collisional and streaming losses, with constant parallel diffusivity $\kappa\sim3\times10^{29}\,\mathrm{cm^2\ s^{-1}}$ chosen to match $\gamma$-ray observations. With this, CRs become more important at larger halo masses and lower redshifts, and dominate the pressure in the CGM in MW-mass halos at $z\lesssim 1-2$. The gas in these ``CR-dominated'' halos differs significantly from runs without CRs: the gas is primarily cool (a few $\sim10^{4}\,$K), and the cool phase is volume-filling and has a thermal pressure below that needed for virial or local thermal pressure balance. Ionization of the ``low'' and ``mid'' ions in this diffuse cool gas is dominated by photo-ionization, with O VI columns $\gtrsim 10^{14.5}\,\mathrm{cm^{-2}}$ at distances $\gtrsim 150\,\mathrm{kpc}$. CR and thermal gas pressure are locally anti-correlated, maintaining total pressure balance, and the CGM gas density profile is determined by the balance of CR pressure gradients and gravity. Neglecting CRs, the same halos are primarily warm/hot ($T\gtrsim 10^{5}\,$K) with thermal pressure balancing gravity, collisional ionization dominates, O VI columns are lower and Ne VIII higher, and the cool phase is confined to dense filaments in local thermal pressure equilibrium with the hot phase.

1909.00283

Solar system chaos and the Paleocene-Eocene boundary age constrained by geology and astronomy

Zeebe, Lourens

Astronomical calculations reveal the solar system's dynamical evolution, including its chaoticity, and represent the backbone of cyclostratigraphy and astrochronology. An absolute, fully calibrated astronomical time scale has hitherto been hampered beyond $\sim$50 Ma, because orbital calculations disagree before that age. Here we present geologic data and a new astronomical solution (ZB18a), showing exceptional agreement from $\sim$58 to 53 Ma. We provide a new absolute astrochronology up to 58 Ma and a new Paleocene-Eocene boundary age (56.01 $\pm$ 0.05 Ma). We show that the Paleocene-Eocene Thermal Maximum (PETM) onset occurred near a 405-kyr eccentricity maximum, suggesting an orbital trigger. We also provide an independent PETM duration (170 $\pm$ 30 kyr) from onset to recovery inflection. Our astronomical solution requires a chaotic resonance transition at $\sim$50 Ma in the solar system's fundamental frequencies.

1909.00285

The Galilean satellites formed slowly from pebbles

Shibaike, et al

It is generally accepted that the four major (Galilean) satellites formed out of the gas disk that accompanied Jupiter's formation. However, understanding the specifics of the formation process is challenging as both small particles (pebbles) as well as the satellites are subject to fast migration processes. Here, we hypothesize a new scenario for the origin of the Galilean system, based on the capture of several planetesimal seeds and subsequent slow accretion of pebbles. To halt migration, we invoke an inner disk truncation radius, and other parameters are tuned for the model to match physical, dynamical, compositional, and structural constraints. In our scenario it is natural that Ganymede's mass is determined by pebble isolation. Our slow-pebble-accretion scenario then reproduces the following characteristics: (1) the mass of all the Galilean satellites; (2) the orbits of Io, Europa, and Ganymede captured in mutual 2:1 mean motion resonances; (3) the ice mass fractions of all the Galilean satellites; (4) the unique ice-rock partially differentiated Callisto and the complete differentiation of the other satellites. Our scenario is unique to simultaneously reproduce these disparate properties.

Day 1622

Monday.

1908.11380
The quenching and morphological evolution of central galaxies is facilitated by the feedback-driven expulsion of circumgalactic gas
Davies, et al

We examine the connection between the properties of the circumgalactic medium (CGM) and the quenching and morphological evolution of central galaxies in the EAGLE and IllustrisTNG simulations. The simulations yield very different median CGM mass fractions, $f_{\rm CGM}$, as a function of halo mass, $M_{200}$, with low-mass haloes being significantly more gas-rich in IllustrisTNG than in EAGLE. Nonetheless, in both cases scatter in $f_{\rm CGM}$ at fixed $M_{200}$ is strongly correlated with the specific star formation rate and the kinematic morphology of central galaxies. The correlations are strongest for $\sim L^\star$ galaxies, corresponding to the mass scale at which expulsive AGN feedback becomes efficient. This feedback elevates the CGM cooling time, preventing gas from accreting onto the galaxy to fuel star formation, and thus establishing a preference for quenched, spheroidal galaxies to be hosted by haloes with low $f_{\rm CGM}$ for their mass. In both simulations, $f_{\rm CGM}$ correlates negatively with the host halo's intrinsic concentration, and hence with its binding energy and formation redshift, primarily because early halo formation fosters the rapid early growth of the central black hole (BH). This leads to a lower $f_{\rm CGM}$ at fixed $M_{200}$ in EAGLE because the BH reaches high accretion rates sooner, whilst in IllustrisTNG it occurs because the central BH reaches the mass threshold at which AGN feedback is assumed to switch from thermal to kinetic injection earlier. Despite these differences, there is consensus from these state-of-the-art simulations that the expulsion of efficiently-cooling gas from the CGM is a crucial step in the quenching and morphological evolution of central galaxies.

1908.11392

VEXAS: the VISTA Extension to Auxiliary Surveys -- Data Release 1: the Southern Galactic Hemisphere

Spiniello, Agnello

We present the first public data release of the VISTA EXtension to Auxiliary Surveys (VEXAS), comprising of 9 cross-matched multi-wavelength photometric catalogs where each object has a match in at least two surveys. We aim at a spatial coverage as uniform as possible in the multi-wavelength sky, with the purpose of providing the astronomical community with reference magnitudes and colours for various scientific uses, including: object classification (e.g. quasars, galaxies, and stars; high-z galaxies, white dwarfs, etc.); photometric redshifts of large galaxy samples; searches of exotic objects such as, for example, extremely red objects and lensed quasars. We have cross-matched the wide-field VISTA catalogs (the VISTA Hemisphere Survey and the VISTA Kilo Degree Infrared Galaxy Survey) with the AllWISE mid-infrared Survey, requiring that a match exists within 10 arcsec. We have further matched this table with X-Ray and radio data (ROSAT, XMM, SUMSS). We also performed a second cross-match between VISTA and AllWISE, with a smaller matching radius (3"), including WISE magnitudes. We have then cross-matched this resulting table ($\approx138\times10^6$ objects) with three photometric wide-sky optical deep surveys (DES, SkyMapper, PanSTARRS). We finally include matches to objects with spectroscopic follow-up by the SDSS and 6dFGS. To demonstrate the power of all-sky multi-wavelength cross-match tables, we show two examples of scientific applications of VEXAS, in particular using the publicly released tables to discover strong gravitational lenses (beyond the reach of previous searches), and to build a statistically large sample of extremely red objects. The VEXAS catalog is currently the widest and deepest, public, optical-to-IR photometric and spectroscopic database in the Southern Hemisphere.

1908.11440

Radio flares from collisions of neutrino stars with interstellar asteroids

Siraj, Loeb

We propose that collisions between neutron stars and interstellar asteroids, such as `Oumuamua, could power observable radio flares in the Milky Way galaxy. We find the rate of such events at $\sim 1 \mathrm{\; Jy}$ to be $\sim 10 \mathrm{\; day^{-1}}$.

1908.11448

Testing the accuracy of halo occupation distribution modeling using hydrodynamic simulations

Beltz-Mohrmann, et al

Halo models provide a simple and computationally inexpensive way to investigate the connection between galaxies and their dark matter haloes. However, these models rely on the assumption that the role of baryons can be easily parametrized in the modelling procedure. We aim to examine the ability of halo occupation distribution (HOD) modelling to reproduce the galaxy clustering found in two different hydrodynamic simulations, Illustris and EAGLE. For each simulation, we measure several galaxy clustering statistics on two different luminosity threshold samples. We then apply a simple five parameter HOD, which was fit to each simulation separately, to the corresponding dark matter only simulations, and measure the same clustering statistics. We find that the halo mass function is shifted to lower masses in the hydrodynamic simulations, resulting in a galaxy number density that is too high when an HOD is applied to the dark matter only simulation. However, the exact way in which baryons alter the mass function is remarkably different in the two simulations. After applying a correction to the halo mass function in each simulation, the HOD is able to accurately reproduce all clustering statistics for the high luminosity sample of galaxies. For the low luminosity sample, we find evidence that in addition to correcting the halo mass function, including spatial, velocity, and assembly bias parameters in the HOD is necessary to accurately reproduce clustering statistics.