Day 1577

Thursday.

1905.12496
The Hubble Constant determined through an inverse distance ladder including quasar time delays and type Ia supernovae
Taubenberger, et al

Context. The precise determination of the present-day expansion rate of the Universe, expressed through the Hubble constant $H_0$, is one of the most pressing challenges in modern cosmology. Assuming flat $\Lambda$CDM, $H_0$ inference at high redshift using cosmic-microwave-background data from Planck disagrees at the 4.4$\sigma$ level with measurements based on the local distance ladder made up of parallaxes, Cepheids and Type Ia supernovae (SNe Ia), often referred to as "Hubble tension". Independent, cosmological-model-insensitive ways to infer $H_0$ are of critical importance. Aims. We apply an inverse-distance-ladder approach, combining strong-lensing time-delay-distance measurements with SN Ia data. By themselves, SNe Ia are merely good relative distance indicators, but by anchoring them to strong gravitational lenses one can obtain an $H_0$ measurement that is relatively insensitive to other cosmological parameters. Methods. A cosmological parameter estimate is performed for different cosmological background models, both for strong-lensing data alone and for the combined lensing + SNe Ia data sets. Results. The cosmological-model dependence of strong-lensing $H_0$ measurements is significantly mitigated through the inverse distance ladder. In combination with SN Ia data, the inferred $H_0$ consistently lies around 73-74 km s$^{-1}$ Mpc$^{-1}$, regardless of the assumed cosmological background model. Our results agree nicely with those from the local distance ladder, but there is a >2$\sigma$ tension with Planck results, and a ~1.5$\sigma$ discrepancy with results from an inverse distance ladder including Planck, Baryon Acoustic Oscillations and SNe Ia. Future strong-lensing distance measurements will reduce the uncertainties in $H_0$ from our inverse distance ladder.

1905.12505

Anisotropies of different mass compositions of cosmic rays

Qiao, et al

The spectral hardenings of cosmic ray nuclei above $\sim 200$ GV followed by softenings around 10 TV, the knee of the all-particle spectrum around PeV energies, as well as the pattern change of the amplitude and phase of the large-scale anisotropies around 100 TeV indicate the complexities of the origin and transportation of Galactic cosmic rays. It has been shown that nearby source(s) are most likely to be the cause of such spectral features of both the spectra and the anisotropies. In this work, we study the anisotropy features of different mass composition (or mass groups) of cosmic rays in this nearby source model. We show that even if the spectral features from the nearby source component is less distinctive compared with the background component from e.g., the population of distant sources, the anisotropy features are more remarkable to be identified. Measurements of the anisotropies of each mass composition (group) of cosmic rays by the space experiments such as DAMPE and HERD and the ground-based experiments such as LHAASO in the near future are expected to be able to critically test this scenario.

Day 1576

Wednesday.

1905.11402
Lyman continuum observations across cosmic time: recent developments, future requirements
McCandliss, et al

Quantifying the physical conditions that allow radiation emitted shortward of the hydrogen ionization edge at 911.7 {\AA} to escape the first collapsed objects and ultimately reionize the universe is a compelling problem for astrophysics. The escape of LyC emission from star-forming galaxies and AGN is intimately tied to the emergence and sustenance of the metagalactic ionizing background that pervades the universe to the present day and in turn is tied to the emergence of structure at all epochs. JWST was built in part to search for the source(s) responsible for reionization, but it cannot observe LyC escape directly, because of the progressive increase in the mean transmission of the intergalactic medium towards the epoch of reionization. Remarkable progress has been made to date in directly detecting LyC leaking from star-forming galaxies using space-based and the ground-based observatories, but there remain significant gaps in our redshift coverage of the phenomenon. Ongoing projects to measure LyC escape at low- and intermediate-z will provide guidance to JWST investigations by analyzing the robustness of a set of proposed LyC escape proxies, and also provide a closeup examination of the physical conditions that favor LyC escape. However, currently available facilities are inadequate for deeply probing LyC escape at the faint end of the galaxy luminosity function. Doing so will require facilities that can detect LyC emission in the restframe to limiting magnitudes approaching 28 $< m^*_{(1+z)900} <$ 32 for $M^*_{(1+z)1500}$ galaxies. The goal of acquiring statistically robust samples for determining LyC luminosity functions across cosmic time will require multi-object spectroscopy from spacebased flagship class and groundbased ELT class telescopes along with ancillary panchromatic imaging and spectroscopy spanning the far-UV to the mid-IR.

1905.11410

Brown dwarf atmospheres as the potentially most detectable and abundant sites for life

Lingam, Loeb

We show that the total habitable volume in the atmospheres of cool brown dwarfs with effective temperatures of $\sim 250$-$350$ K is possibly larger by two orders of magnitude than that of Earth-like planets. We also study the role of aerosols, nutrients and photosynthesis in facilitating life in brown dwarf atmospheres. Our predictions might be testable through searches for spectral edges in the near-infrared and chemical disequilibrium in the atmospheres of nearby brown dwarfs that are either free-floating or within $\sim 10$ AU of stars. For the latter category, we find that the James Webb Space Telescope (JWST) may be able to achieve a signal-to-noise ratio of $\sim 5$ after a few hours integration per source for the detection of biogenic spectral features in $\sim 10^3$ cool brown dwarfs.

1905.11636

Effects of baryons on weak lensing peak statistics

Weiss, et al

Upcoming weak-lensing surveys have the potential to become leading cosmological probes provided all systematic effects are under control. Recently, the ejection of gas due to feedback energy from active galactic nuclei (AGN) has been identified as major source of uncertainty, challenging the success of future weak-lensing probes in terms of cosmology. In this paper we investigate the effects of baryons on the number of weak-lensing peaks in the convergence field. Our analysis is based on full-sky convergence maps constructed via light-cones from $N$-body simulations, and we rely on the baryonic correction model of Schneider et al. (2019) to model the baryonic effects on the density field. As a result we find that the baryonic effects strongly depend on the Gaussian smoothing applied to the convergence map. For a DES-like survey setup, a smoothing of $\theta_{\kappa}\gtrsim8$ arcmin is sufficient to keep the baryon signal below the expected statistical error. Smaller smoothing scales lead to a significant suppression of high peaks with $\kappa> 0.2$, while lower peaks are not affected. The situation is more severe for a Euclid-like setup, where a smoothing of $\theta_{\kappa}\gtrsim16$ arcmin is required to keep the baryonic suppression signal below the statistical error. Smaller smoothing scales require a full modelling of baryonic effects since both low and high peaks are strongly affected by baryonic feedback.

1905.11721

The large-scale general-relativistic correction for Newtonian mocks

Adamek, Fidler

We clarify the subtle issue of finding the correct mapping of Newtonian simulations to light-cone observables at very large distance scales. A faithful general-relativistic interpretation specifies a gauge, i.e. a chart that relates the simulation data to points of the space-time manifold. It has already been pointed out that the implicit gauge choice of Newtonian simulations is indeed different from the Poisson gauge that is commonly adopted for relativistic calculations, the difference being most significant at large scales. It is therefore inconsistent, for example, to predict weak-lensing observables from simulations unless this gauge issue is properly accounted for. Using perturbation theory as well as fully relativistic N-body simulations we quantify the systematic error introduced this way, and we discuss several solutions that would render the calculations relativistically self-consistent.

1905.11886

Towards a model-independent measurement of the halo mass function with observables

Dong, et al

In the CDM paradigm, the halo mass function is a sensitive probe of the cosmic structure. In observations, halo mass is typically estimated from its relation with other observables. The resulting halo mass function is subject to systematic bias, such as the Eddington bias, due to the scatter or uncertainty in the observable - mass relation. Exact correction for the bias is not easy, as predictions for the observables are typically model-dependent in simulations. In this paper, we point out an interesting feature in the halo mass function of the concordence $\Lambda$CDM model: the total halo mass within each evenly-spaced logarithmic mass bin is approximately the same over a large mass range. We show that this property allows us to construct an almost bias-free halo mass function using only an observable (as a halo mass estimator) and stacked weak lensing measurements as long as the scatter between the true halo mass and the observable-inferred mass has a stable form in logarithmic units. The method is not sensitive to the form of the mass-observable relation. We test the idea using cosmological simulations, and show that the method performs very well for realistic observables.

Day 1575

Tuesday.

1905.10522
An extended catalog of galaxy-galaxy strong gravitational lenses discovered in DES using convolutional neural networks
Jacobs, et al

We search Dark Energy Survey (DES) Year 3 imaging for galaxy-galaxy strong gravitational lenses using convolutional neural networks, extending previous work with new training sets and covering a wider range of redshifts and colors. We train two neural networks using images of simulated lenses, then use them to score postage stamp images of 7.9 million sources from the Dark Energy Survey chosen to have plausible lens colors based on simulations. We examine 1175 of the highest-scored candidates and identify 152 probable or definite lenses. Examining an additional 20,000 images with lower scores, we identify a further 247 probable or definite candidates. After including 86 candidates discovered in earlier searches using neural networks and 26 candidates discovered through visual inspection of blue-near-red objects in the DES catalog, we present a catalog of 511 lens candidates.

1905.10739

Connecting galaxy structure and star formation: the role of environment in formation of S0 galaxies

Mishra, et al

In this work, we investigate the reason behind the increased occurrence of S0 galaxies in high density environments. Our sample comprises of $\sim$ 2500 spiral and $\sim$ 2000 S0 galaxies spanning a wide range of environments. Dividing the galaxies into categories of classical and pseudobulge hosting spiral and S0 galaxies, we have studied their properties as a function of the environment. We find that the fraction of pseudobulge hosting disc galaxies decreases with increase in density. The classical bulge hosting spirals and S0 galaxies follow a similar trend in less dense environments but towards higher densities, we observe an increase in the fraction of classical bulge host S0 galaxies at the expense of spirals. Comparing the structural and the star formation properties of galaxies on the size-mass and $NUV-r$ colour-mass planes respectively, we infer that classical bulge hosting spirals are likely to get transformed into S0 morphology. We notice a trend of galaxy structure with environment such that the fraction of classical bulge hosting spiral galaxies is found to increase with environment density. We also find that among classical bulge hosting spirals, the fraction of quenched galaxies increases in denser environments. We surmise that the existence of more classical bulge hosting spirals galaxies and more efficient quenching leads to the observed increased occurrence of S0 galaxies in high density environments. The relation between galaxy structure and environment also exists for the disc galaxies irrespective of their visual morphology, which is driven mainly by halo mass.

1905.10779

The star formation history n the solar neighborhood as told by massive white dwarfs

Isern

White dwarfs are the remnants of low and intermediate mass stars. Because of electron degeneracy, their evolution is just a simple gravothermal process of cooling. Recently, thanks to Gaia data, it has been possible to construct the luminosity function of massive (0.9 < M/Msun < 1.1) white dwarfs in the solar neighborhood (d < 100 pc). Since the lifetime of their progenitors is very short, the birth times of both, parents and daughters, are very close and allow to reconstruct the (effective) star formation rate. This rate started growing from zero during the early Galaxy and reached a maximum 6-7 Gyr ago. It declined and ~5 Gyr ago started to climb once more reaching a maximum 2 - 3 Gyr in the past and decreased since then. There are some traces of a recent star formation burst, but the method used here is not appropriate for recently born white dwarfs.

Day 1574

Monday.

1905.09839
A fundamental test for stellar feedback recipes in galaxy simulations
Fujimoto, et al

Direct comparisons between galaxy simulations and observations that both reach scales < 100 pc are strong tools to investigate the cloud-scale physics of star formation and feedback in nearby galaxies. Here we carry out such a comparison for hydrodynamical simulations of a Milky Way-like galaxy, including stochastic star formation, HII region and supernova feedback, and chemical post-processing at 8 pc resolution. Our simulation shows excellent agreement with almost all kpc-scale and larger observables, including total star formation rates, radial profiles of CO, HI, and star formation through the galactic disc, mass ratios of the ISM components, both whole-galaxy and resolved Kennicutt-Schmidt relations, and giant molecular cloud properties. However, we find that our simulation does not reproduce the observed de-correlation between tracers of gas and star formation on < 100 pc scales, known as the star formation 'uncertainty principle', which indicates that observed clouds undergo rapid evolutionary lifecycles. We conclude that the discrepancy is driven by insufficiently-strong pre-supernova feedback in our simulation, which does not disperse the surrounding gas completely, leaving star formation tracer emission too strongly associated with molecular gas tracer emission, inconsistent with observations. This result implies that the cloud-scale de-correlation of gas and star formation is a fundamental test for feedback prescriptions in galaxy simulations, one that can fail even in simulations that reproduce all other macroscopic properties of star-forming galaxies.

1905.09908

Kinematics of multiple stellar populations in globular clusters with Gaia

Cordoni, et al

The internal dynamics of multiple stellar populations in Globular Clusters (GCs) provide unique constraints on the physical processes responsible for their formation. Specifically, the present-day kinematics of cluster stars, such as rotation and velocity dispersion, seem to be related to the initial configuration of the system. In a recent work, we analyzed for the first time the kinematics of the different stellar populations in NGC 0104 over a large field of view, exploiting the Gaia Data Release 2 proper motions combined with multi-band ground-based photometry. In this paper, we extend this analysis to six GCs, namely NGC 0288, NGC 5904, NGC 6121, NGC 6752 and NGC 6838 and further explore NGC 0104. Among the analyzed clusters only NGC 0104 and NGC 5904 show significant rotation in the plane of the sky. By separating our sample in 1G and 2G stars we find that overall these two populations exhibit a similar rotation pattern in NGC 0104. However, some hints of different rotations between 1G and 2G stars are observed in the external regions of this cluster. Interestingly, 1G and 2G stars in NGC 5904 exhibit different rotation curves, with distinct phases. The radial components of the motion of 1G and 2G stars show different radial trends, in contrast with what is observed in most of the other clusters. There is no evidence for rotation among the selected 1G and 2G stars of the remaining clusters. The analysis of the velocity-dispersion profiles of multiple populations confirms that 2G stars of NGC 0104 show stronger anisotropy than the 1G.

1905.09920

Multiwavelength cluster mass estimates and machine learning

Cohn, Battaglia

One emerging application of machine learning methods is the inference of galaxy cluster masses. Often cluster mass predictions are made from observables by fitting or deriving scaling relations; if multiwavelength measurements are available, these scaling relation based estimates are combined into a likelihood. Here, machine learning is used in a simulation to instead directly combine five multiwavelength measurements to obtain cluster masses. Comparisons of the contributions of each observable to the accuracy of the resulting mass measurement are made using model-agnostic Importance Permutation values, as well as by brute force comparison of different combinations of observables. As machine learning relies upon the accuracy of the training set in capturing the observables, their correlations, and the observational selection function, and the training set originates from simulations, a few ways of testing whether a simulation and observations are consistent are explored as well.

1905.10312

Evidence for disc regulation in the lowest-mass stars of the young stellar cluster NGC 2264

Orcajo, et al

In the pre-main-sequence stage, star-disc interactions have been shown to remove stellar angular momentum and regulate the rotation periods of stars with M2 and earlier spectral types. Whether disc regulation also extends to stars with later spectral types still remains a matter of debate. Here we present a star-disc interaction study in a sample of over 180 stars with spectral types M3 and later (corresponding to stellar masses $\leq 0.3 M_\odot$) in young stellar cluster NGC 2264. Combining rotation periods from the literature, new and literature spectral types, and newly presented deep Spitzer observations, we show that stars with masses below 0.3 $M_\odot$ with discs also rotate slower than stars without a disc in the same mass regime. Our results demonstrate that disc-regulation still operates in these low-mass stars, although the efficiency of this process might be lower than in higher-mass objects. We confirm that stars with spectral types earlier and later than M2 have distinct period distributions and that stars with spectral types M5 and later rotate even faster M3 and M4-type stars.

1905.10359

Stellar systems following the $R^{1/M}$ luminosity law.  III.  Photometric, intrinsic and dynamical properties for all S\'sic indices

Baes, Ciotti

The S\'ersic or $R^{1/m}$ model has become the de facto standard model to describe the surface brightness profiles of early-type galaxies and the bulges of spiral galaxies. The photometric, intrinsic and dynamical properties of this model have been investigated, but mainly for fairly large S\'ersic indices $m$. For small values of $m$, appropriate for low-mass and dwarf ellipticals, a detailed investigation of these properties is still lacking. In this study, we use a combination of numerical and analytical techniques to investigate the S\'ersic model over the entire range of S\'ersic parameters, focusing on the small $m$ regime, where a number of interesting and surprising properties are found. For all values $m<1$, the model is characterised by a finite central luminosity density, and for $m<\tfrac12$, even a central depression in the luminosity density profile. This behaviour translates to the dynamical properties: we show that all S\'ersic models with $m \geqslant\tfrac12$ can be supported by an isotropic velocity dispersion tensor, and that these isotropic models are stable to both radial and non-radial perturbations. The models with $m < \tfrac12$, on the other hand, cannot be supported by an isotropic velocity dispersion tensor.

Day 1573

Friday.

1905.09288
$\bar T$: A new cosmological parameter?
Yoo, et al

The background photon temperature $\bar T$ is one of the fundamental cosmological parameters. Despite its significance, $\bar T$ has never been allowed to vary in the data analysis, owing to the precise measurement of the comic microwave background (CMB) temperature by COBE FIRAS. However, even in future CMB experiments, $\bar T$ will remain unknown due to the unknown monopole contribution $\Theta_0$ at our position to the observed (angle-averaged) temperature $\langle T\rangle^{\rm obs}$. By fixing $\bar T\equiv\langle T\rangle^{\rm obs}$, the standard analysis underestimates the error bars on cosmological parameters, and the best-fit parameters obtained in the analysis are biased in proportion to the unknown amplitude of $\Theta_0$. Using the Fisher formalism, we find that these systematic errors are smaller than the error bars from the $Planck$ satellite. However, with $\bar T\equiv\langle T\rangle^{\rm obs}$, these systematic errors will always be present and irreducible, and future cosmological surveys might misinterpret the measurements.

1905.09353

Physical correlations of the scatter between galaxy mass, stellar content, and halo mass

Bradshaw, et al

We use the UniverseMachine to analyze the source of scatter between the central galaxy mass, the total stellar mass in the halo, and the dark matter halo mass. We also propose a new halo mass estimator, the cen+N mass: the sum of the stellar mass of the central and the N most massive satellites. We show that, when real space positions are perfectly known, the cen+N mass has scatter competitive with that of richness-based estimators. However, in redshift space, the cen+N mass suffers less from projection effects in the UniverseMachine model. The cen+N mass is therefore a viable low scatter halo mass estimator, and should be considered an important tool to constrain cosmology with upcoming spectroscopic data from DESI. We analyze the scatter in stellar mass at fixed halo mass and show that the total stellar mass in a halo is uncorrelated with secondary halo properties, but that the central stellar mass is a function of both halo mass and halo age. This is because central galaxies in older halos have had more time to grow via accretion. If the UniverseMachine model is correct, accurate galaxy-halo modeling of mass selected samples therefore needs to consider halo age in addition to mass.

1905.09664

Effects of dark matter pressure on the ellipticity of cosmic voids

Rezaei

The dark matter in or around the cosmic voids affects their shapes. The thermodynamical properties of dark matter can alter the ellipticity of cosmic voids. Here, applying the dark matter equation of state from the pseudo-isothermal density profile of galaxies, we explore the shapes of cosmic voids with the non zero pressure dark matter in different cosmological models. For this purpose, the linear growth of density perturbation in the presence of dark matter pressure is calculated. In addition, the matter transfer function considering the dark matter pressure, as well as the linear matter power spectrum in the presence of the dark matter pressure are presented. Employing these results, the probability density distribution for the ellipticity of cosmic voids with the non zero pressure dark matter is calculated. Our calculations verify that the dark matter pressure leads to more spherical shapes for the cosmic voids.

1905.09781

Model-independent determination of $H_0$ and $\Omega_{K0}$ from strong lensing and type Ia supernovae

Collett, et al

We present the first determination of the Hubble constant $H_0$ from strong lensing time delay data and type Ia supernova luminosity distances that is independent of the cosmological model. We also determine the spatial curvature model-independently. We assume that light propagation over long distances is described by the FLRW metric and geometrical optics holds, but make no assumption about the contents of the Universe or the theory of gravity on cosmological scales. We find $H_0=75.7^{+4.5}_{-4.4}$ km/s/Mpc and $\Omega_{K0}=0.12^{+0.27}_{-0.25}$. This is a 6\% determination of $H_0$. A weak prior from the cosmic microwave background on the distance to the last scattering surface improves this to $H_0=76.8^{+4.2}_{-3.8}$ km/s/Mpc and $\Omega_{K0}=0.18^{+0.25}_{-0.18}$. Assuming zero spatial curvature, we get $H_0=74.2^{+3.0}_{-2.9}$ km/s/Mpc, a precision of $4\%$. The measurements also provide a consistency test of the FLRW metric: we find no evidence against it.

1905.09800

Gaussbock: fast parallel-iterative cosmological parameter estimation with Bayesian nonparametrics

Moews, Zuntz

We present and apply Gaussbock, a new embarrassingly parallel iterative algorithm for cosmological parameter estimation designed for an era of cheap parallel computing resources. Gaussbock uses Bayesian nonparametrics and truncated importance sampling to accurately draw samples from posterior distributions with an orders-of-magnitude speed-up in wall time over alternative methods. Contemporary problems in this area often suffer from both increased computational costs due to high-dimensional parameter spaces and consequent excessive time requirements, as well as the need for fine tuning of proposal distributions or sampling parameters. Gaussbock is designed specifically with these issues in mind. We explore and validate the performance and convergence of the algorithm on a fast approximation to the Dark Energy Survey Year 1 (DES Y1) posterior, finding reasonable scaling behavior with the number of parameters. We then test on the full DES Y1 posterior using large-scale supercomputing facilities, and recover reasonable agreement with previous chains, although the algorithm can underestimate the tails of poorly-constrained parameters. In addition, we provide the community with a user-friendly software tool for accelerated cosmological parameter estimation based on the methodology described in this paper.

Day 1572

Thursday.

1905.08799
Revealing the galaxy-halo connection in IllustrisTNG
Bose, Eisenstein, et al

We use the IllustrisTNG (TNG) simulations to explore the galaxy-halo connection as inferred from state-of-the-art cosmological, magnetohydrodynamical simulations. With the high mass resolution and large volume achieved by combining the 100 Mpc (TNG100) and 300 Mpc (TNG300) volumes, we establish the mean occupancy of central and satellite galaxies and their dependence on the properties of the dark matter haloes hosting them. We derive best-fitting HOD parameters from TNG100 and TNG300 for target galaxy number densities of $\bar{n}_g = 0.032\,h^3$Mpc$^{-3}$ and $\bar{n}_g = 0.016\,h^3$Mpc$^{-3}$, respectively, corresponding to a minimum galaxy stellar mass of $M_\star\sim1.9\times10^9\,{\rm M}_\odot$ and $M_\star\sim3.5\times10^9\,{\rm M}_\odot$, respectively, in hosts more massive than $10^{11}\,{\rm M}_\odot$. Consistent with previous work, we find that haloes located in dense environments, with low concentrations, later formation times, and high angular momenta are richest in their satellite population. At low mass, highly-concentrated haloes and those located in overdense regions are more likely to contain a central galaxy. The degree of environmental dependence is sensitive to the definition adopted for the physical boundary of the host halo. We examine the extent to which correlations between galaxy occupancy and halo properties are independent and demonstrate that HODs predicted by halo mass and present-day concentration capture the qualitative dependence on the remaining halo properties. At fixed halo mass, concentration is a strong predictor of the stellar mass of the central galaxy, which may play a defining role in the fate of the satellite population. The radial distribution of satellite galaxies, which exhibits a universal form across a wide range of host halo mass, is described accurately by the best-fit NFW density profile of their host haloes.

1905.08803

The Galactic mid plane is not a plane: implications for dynamical analysis with Gaia data and beyond

Beane, et al

Orbital properties of stars, computed from their six-dimensional phase space measurements and an assumed Galactic potential, are used to understand the structure and evolution of the Galaxy. Stellar actions, computed from orbits, have the attractive quality of being invariant under certain assumptions and are therefore used as quantitative labels of a star's orbit. We report a subtle but important systematic error that is induced in the actions as a consequence of local midplane variations expected for the Milky Way. This error is difficult to model because it is non-Gaussian and bimodal, with neither mode peaking on the null value. An offset in the vertical position of the Galactic midplane of $\sim15\,\text{pc}$ for a thin disk-like orbit or $\sim 120\,\text{pc}$ for a thick disk-like orbit induces a $25\%$ systematic error in the vertical action $J_z$. In FIRE simulations of Milky Way-mass galaxies, these variations are on the order of $\sim100\,\text{pc}$ at the solar circle. From observations of the mean vertical velocity variation of $\sim5\text{--}10\,\text{km}\,\text{s}^{-1}$ with radius, we estimate that the Milky Way midplane variations are $\sim60\text{--}170\,\text{pc}$, consistent with three-dimensional dust maps. Action calculations and orbit integrations, which assume the global and local midplanes are identical, are likely to include this induced error, depending on the volume considered. Variation in the local standard of rest or distance to the Galactic center causes similar issues. The variation of the midplane must be taken into account when performing dynamical analysis across the large regions of the disk accessible to Gaia and future missions.

1905.08892

From Centaurs to comets - 40 years

Pixinho, et al

In 1977, while Apple II and Atari computers were being sold, a tiny dot was observed in an inconvenient orbit. The minor body 1977 UB, to be named (2060) Chiron, with an orbit between Saturn and Uranus, became the first Centaur, a new class of minor bodies orbiting roughly between Jupiter and Neptune. The observed overabundance of short-period comets lead to the downfall of the Oort Cloud as exclusive source of comets and to the rise of the need for a Trans-Neptunian comet belt. Centaurs were rapidly seen as the transition phase between Kuiper Belt Objects (KBOs), also known as Trans-Neptunian Objects (TNOs) and the Jupiter-Family Comets (JFCs). Since then, a lot more has been discovered about Centaurs: they can have cometary activity and outbursts, satellites, and even rings. Over the past four decades since the discovery of the first Centaur, rotation periods, surface colors, reflectivity spectra and albedos have been measured and analyzed. However, despite such a large number of studies and complementary techniques, the Centaur population remains a mystery as they are in so many ways different from the TNOs and even more so from the JFCs.

1905.08951

Globular cluster formation from colliding substructure

Madau, et al

We investigate a scenario where the formation of Globular Clusters (GCs) is triggered by high-speed collisions between infalling atomic-cooling subhalos during the assembly of the main galaxy host, a special dynamical mode of star formation that operates at high gas pressures and is intimately tied to LCDM hierarchical galaxy assembly. The proposed mechanism would give origin to "naked" globulars, as colliding dark matter subhalos and their stars will simply pass through one another while the warm gas within them clashes at highly supersonic speed and decouples from the collisionless component, in a process reminiscent of the Bullet galaxy cluster. We find that the resulting shock-compressed layer cools on a timescale that is tipically shorter than the crossing time, first by atomic line emission and then via fine-structure metal-line emission, and is subject to gravitational instability and fragmentation. Through a combination of kinetic theory approximation and high-resolution N-body simulations, we show that this model may produce: (a) a GC number-halo mass relation that is linear down to dwarf galaxy scales and agrees with the trend observed over five orders of magnitude in galaxy mass; (b) a population of old globulars with a median age of 12 Gyr and an age spread similar to that observed; (c) a spatial distribution that is biased relative to the overall mass profile of the host. This is because, in an inelastic collision, the splash remnant will lose orbital energy and fall deeper into the Galactic potential rather than sharing the orbits of the progenitor subhalos; and (d) a bimodal metallicity distribution with a spread similar to that observed in massive galaxies. Additional, hydrodynamic simulations of subhalo-subhalo high-speed impacts should be performed to further validate a collision-driven scenario for the formation of GCs.

1905.08991

Weak lensing measurement of filamentary structure with the SDSS BOSS and Subr Hyper Suprime-Cam Data

Kondo, et al

We report the weak lensing measurement of filaments between Sloan Digital Sky Survey (SDSS) III/Baryon Oscillation Spectroscopic Survey (BOSS) CMASS galaxy pairs at $z\sim0.55$, using the Subaru Hyper Suprime-Cam (HSC) first-year galaxy shape catalogue. Despite of the small overlap of $140$ deg$^2$ between these surveys we detect the filament lensing signal at 3.9$\sigma$ significance, which is the highest signal-to-noise lensing measurement of filaments between galaxy-scale halos at this redshift range. We derive a theoretical prediction and covariance using mock catalogues based on full-sky ray-tracing simulations. We find that the intrinsic scatter of filament properties and the fluctuations in large scale structure along the line-of-sight are the primary component of the covariance and the intrinsic shape noise from source galaxies no longer limits our lensing measurement. This fact demonstrates the statistical power of the HSC survey due to its deep observations and high number density of source galaxies. Our result is consistent with the theoretical prediction and supports the "thick" filament model. As the HSC survey area increases, we will be able to study detailed filament properties such as the dark matter distributions and redshift evolution of filaments.

1905.08993

Morphological diversity of spiral galaxies originating int the cold gas inflow from cosmic webs

Noguchi

Spiral galaxies comprise three major structural components; thin discs, thick discs, and central bulges. Relative dominance of these components is known to correlate with the total mass of the galaxy, and produces a remarkable morphological variety of spiral galaxies. Although there are many formation scenarios regarding individual components, no unified theory exists which explains this systematic variation. The cold-flow hypothesis predicts that galaxies grow by accretion of cold gas from cosmic webs (cold accretion) when their mass is below a certain threshold, whereas in the high-mass regime the gas that entered the dark matter halo is first heated by shock waves to high temperatures and then accretes to the forming galaxy as it cools emitting radiation (cooling flow). This hypothesis also predicts that massive galaxies at high redshifts have a hybrid accretion structure in which filaments of cold inflowing gas penetrate surrounding hot gas. In the case of Milky Way, the previous study suggested that the cold accretion created its thick disc in early times and the cooling flow formed the thin disc in later epochs. Here we report that extending this idea to galaxies with various masses and associating the hybrid accretion with the formation of bulges reproduces the observed mass-dependent structures of spiral galaxies: namely, more massive galaxies have lower thick disc mass fractions and higher bulge mass fractions. The proposed scenario predicts that thick discs are older in age and poorer in iron than thin discs, the trend observed in the Milky Way (MW) and other spiral galaxies.

Day 1571

Wednesday.

1905.08262
Cosmological information contents on the light-cone
Yoo, et al

We develop a theoretical framework to describe the cosmological observables on the past light cone such as the luminosity distance, weak lensing, galaxy clustering, and the cosmic microwave background anisotropies. We consider that all the cosmological observables include not only the background quantity, but also the perturbation quantity, and they are subject to cosmic variance, which sets the fundamental limits on the cosmological information that can be derived from such observables, even in an idealized survey with an infinite number of observations. To quantify the maximum cosmological information content, we apply the Fisher information matrix formalism and spherical harmonic analysis to cosmological observations, in which the angular and the radial positions of the observables on the light cone carry different information. We discuss the maximum cosmological information that can be derived from five different observables: (1) type Ia supernovae, (2) cosmic microwave background anisotropies, (3) weak gravitational lensing, (4) local baryon density, and (5) galaxy clustering. We compare our results with the cosmic variance obtained in the standard approaches, which treat the light cone volume as a cubic box of simultaneity. We discuss implications of our formalism and ways to overcome the fundamental limit.

1905.08754

Towards determining the neutrino mass hierarchy: weak lensing and galaxy clustering forecasts with baryons and intrinsic alignments

Copeland, Taylor, Hall

The capacity of Stage IV lensing surveys to measure the neutrino mass sum and differentiate between the normal and inverted mass hierarchies depends on the impact of nuisance parameters describing small-scale baryonic astrophysics and intrinsic alignments. For a Euclid-like survey, we perform the first combined weak lensing and galaxy clustering Fisher analysis with baryons, intrinsic alignments, and massive neutrinos for both hierarchies. We use a matter power spectrum generated from a halo model that captures the impact of baryonic feedback and adiabatic contraction. For weak lensing, we find that baryons cause severe degradation to forecasts of the neutrino mass sum, $\Sigma$, approximately doubling $\sigma_{\Sigma}$. We show that including galaxy clustering constraints from Euclid and BOSS, and cosmic microwave background (CMB) Planck priors, can reduce this degradation to $\sigma_{\Sigma}$ to 9% and 16% for the normal and inverted hierarchies respectively. The combined forecasts, $\sigma_{\Sigma_{\rm{NH}}}=0.034\, \rm{eV}$ and $\sigma_{\Sigma_{\rm{IH}}}=0.034\, \rm{eV}$, preclude a meaningful distinction of the hierarchies but could be improved upon with future CMB priors on $n_s$ and information from neutrinoless double beta decay to achieve a 2$\sigma$ distinction. The effect of intrinsic alignments on forecasts is shown to be minimal, with $\sigma_{\Sigma}$ even experiencing mild improvements due to information from the intrinsic alignment signal. We find that while adiabatic contraction and intrinsic alignments will require careful calibration to prevent significant biasing of $\Sigma$, there is less risk presented by feedback from energetic events like AGN and supernovae.

Day 1570

Monday.  Tuesday.

1905.06954
Prospects for observing the cosmic web in Lyman-{\alpha} emission
Wistok, et al

Mapping the intergalactic medium (IGM) in Lyman-$\alpha$ emission would yield unprecedented tomographic information on the large-scale distribution of baryons and potentially provide new constraints on the UV background and various feedback processes relevant for galaxy formation. Here, we examine the Lyman-$\alpha$ emission of the moderate-density IGM due to collisional excitations and recombinations in the presence of a UV background in the Sherwood simulation suite. We focus on large-scale-structure filaments in which Lyman-$\alpha$ radiative transfer effects are expected to be moderate. At low density the emission is primarily due to fluorescent re-emission of the ionising UV background due to recombinations, while collisional excitations dominate at higher densities. We discuss prospects of current and future observational facilities to detect this emission and focus on VLT/MUSE for a more detailed sensitivity analysis. We construct mock MUSE observations resembling the MUSE Hubble Deep Field South from our simulations and show that our predictions are consistent with recent analyses of diffuse Lyman-$\alpha$ emission using MUSE at redshifts $3 < z < 6$. We find that it should be most feasible to detect the Lyman-$\alpha$ emission from filaments in the IGM in overdense regions, somewhat surprisingly towards the high-redshift end ($z \gtrsim 4.5$) accessible by MUSE, and if narrowband widths close to the MUSE spectral resolution limit are used. This is due to the higher densities and lower temperatures in the IGM at higher redshift. High-redshift protoclusters therefore appear to be the ideal environment to observe filamentary structures in the IGM in Lyman-$\alpha$ emission.

1905.07024

Probing galaxy assembly bias in BOSS galaxies using void probabilities

Walsh, Tinker

We measure the void probability function (VPF) of galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS). The VPF provides complementary information to standard two-point statistics in that it is sensitive to galaxy bias in the most extreme underdensities in the cosmic web. Thus the VPF is ideal for testing whether halo occupation of galaxies depends on large-scale density, an effect known as galaxy assembly bias. We find that standard HOD model---one parameterized by halo mass only---fit only to the two-point function, accurately predicts the VPF. Additionally, for HOD models where density dependence is explicitly incorporated, the best-fit models fit to the combination of the correlation function and the VPF have zero density dependence. Thus galaxy assembly bias is not a strong source of systematic uncertainty when modeling the clustering of massive galaxies.

1905.07070
The detectability and characterization of the TRAPPIST-1 exoplanet atmospheres with JWST
Lustig-Yaeger, et al

The James Webb Space Telescope (JWST) will offer the first opportunity to characterize terrestrial exoplanets with sufficient precision to identify high mean molecular weight atmospheres, and TRAPPIST-1's seven known transiting Earth-sized planets are particularly favorable targets. To assist community preparations for JWST, we use simulations of plausible post-ocean-loss and habitable environments for the TRAPPIST-1 exoplanets, and test simulations of all bright object time series spectroscopy modes and all MIRI photometry filters to determine optimal observing strategies for atmospheric detection and characterization using both transmission and emission observations. We find that transmission spectroscopy with NIRSpec Prism is optimal for detecting terrestrial, CO2 containing atmospheres, potentially in fewer than 10 transits for all seven TRAPPIST-1 planets, if they lack high altitude aerosols. If the TRAPPIST-1 planets possess Venus-like H2SO4 aerosols, up to 12 times more transits may be required to detect atmospheres. We present optimal instruments and observing modes for the detection of individual molecular species in a given terrestrial atmosphere and an observational strategy for discriminating between evolutionary states. We find that water may be prohibitively difficult to detect in both Venus-like and habitable atmospheres due to its presence lower in the atmosphere where transmission spectra are less sensitive. Although the presence of biogenic O2 and O3 will be extremely challenging to detect, abiotically produced oxygen from past ocean loss may be detectable for all seven TRAPPIST-1 planets via O2-O2 collisionally-induced absorption at 1.06 and 1.27 microns, or via NIR O3 features for the outer three planets. Our results constitute a suite of hypotheses on the nature and detectability of highly-evolved terrestrial exoplanet atmospheres that may be tested with JWST.

1905.07343

Does the evolution of complex life depend on the stellar spectral energy distribution?

Haqq-Misra

This paper presents the proportional evolutionary time hypothesis, which posits that the mean time required for the evolution of complex life is a function of stellar mass. The "biological available window" is defined as the region of a stellar spectrum between 200 to 1200 nm that generates free energy for life. Over the $\sim$4 Gyr history of Earth, the total energy incident at the top of the atmosphere and within the biological available window is $\sim$10$^{34}$ J. The hypothesis assumes that the rate of evolution from the origin of life to complex life is proportional to this total energy, which would suggest that planets orbiting other stars should not show signs of complex life if the total energy incident on the planet is below this energy threshold. The proportional evolutionary time hypothesis predicts that late K- and M-dwarf stars (M < 0.7 M$_{\odot}$) are too young to host any complex life at the present age of the universe. F-, G-, and early K-dwarf stars (M > 0.7 M$_{\odot}$) represent the best targets for the next generation of space telescopes to search for spectroscopic biosignatures indicative of complex life.

1905.07410

Cosmic inference: constraining parameters with observations and highly limited number of simulations

Takhtaganov, et al

Cosmological probes pose an inverse problem where the measurement result is obtained through observations, and the objective is to infer values of model parameters which characterize the underlying physical system -- our Universe. Modern cosmological probes increasingly rely on measurements of the small-scale structure, and the only way to accurately model physical behavior on those scales, roughly 65 Mpc/h or smaller, is via expensive numerical simulations. In this paper, we provide a detailed description of a novel statistical framework for obtaining accurate parameter constraints by combining observations with a very limited number of cosmological simulations. The proposed framework utilizes multi-output Gaussian process emulators that are adaptively constructed using Bayesian optimization methods. We compare several approaches for constructing multi-output emulators that enable us to take possible inter-output correlations into account while maintaining the efficiency needed for inference. Using Lyman alpha forest flux power spectrum, we demonstrate that our adaptive approach requires considerably fewer --- by a factor of a few in Lyman alpha P(k) case considered here --- simulations compared to the emulation based on Latin hypercube sampling, and that the method is more robust in reconstructing parameters and their Bayesian credible intervals.

1905.07557

Testing halo assembly bias using galaxy clusters

Sunayama, More

We critically examine the methodology behind the claimed observational detection of halo assembly bias using optically selected galaxy clusters by Miyatake et al. (2016) and More et al. (2016). We mimic the optical cluster detection algorithm and apply it to two different mock catalogs generated from the Millennium simulation galaxy catalog, one in which halo assembly bias signal is present, while the other in which the assembly bias signal has been expressly erased. We split each of these cluster samples into two using the average cluster-centric distance of the member galaxies to measure the difference in the clustering strength of the subsamples with respect to each other. We observe that the subsamples split by cluster-centric radii show differences in clustering strength, even in the catalog where the true assembly bias signal was erased. We show that this is a result of contamination of the member galaxy sample from interlopers along the line-of-sight. This undoubtedly shows that the particular methodology adopted in the previous studies cannot be used to claim a detection of the assembly bias signal. We figure out the tell-tale signatures of such contamination, and show that the observational data also shows similar signatures. Furthermore, we also show that projection effects in optical galaxy clusters can bias the inference of the 3-dimensional edges of galaxy clusters (splashback radius), so appropriate care should be taken while interpreting the splashback radius of optical clusters.

1905.07803

The matter fluctuation amplitude inferred from the weak lensing power spectrum and correlation function in CHFTLenS data

Lu, Haiman

Based on the cosmic shear data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), Kilbinger et al. (2013) obtained a constraint on the amplitude of matter fluctuations of $\sigma_8({\Omega_\mathrm{m}}/0.27)^{0.6}=0.79\pm0.03$ from the two-point correlation function (2PCF). This is $\approx3\sigma$ lower than the value $0.89\pm0.01$ derived from Planck data on cosmic microwave background (CMB) anisotropies. On the other hand, based on the same CFHTLenS data, but using the power spectrum, and performing a different analysis, Liu et al. (2015) obtained the higher value of $\sigma_8({\Omega_\mathrm{m}}/0.27)^{0.64}=0.87^{+0.05}_{-0.06}$. We here investigate the origin of this difference, by performing a fair side-by-side comparison of the 2PCF and power spectrum analyses on CFHTLenS data. We find that these two statistics indeed deliver different results, even when applied to the same data in an otherwise identical procedure. We identify excess power in the data on small scales ($\ell>5,000$) driving the larger values inferred from the power spectrum. We speculate on the possible origin of this excess small-scale power.

1905.08017

Occurrence of great magnetic storms on 6-8 March 1582

Hattori, Hayakawa, Ebihara

Although knowing the occurrence frequency of severe space weather events is important for a modern society, it is insufficiently known due to the lack of magnetic or sunspot observations, before the Carrington event in 1859 known as one of the largest events during the last two centuries. Here, we show that a severe magnetic storm occurred on 8 March 1582 based on auroral records in East Asia. The equatorward boundary of auroral visibility reached 28.8{\deg} magnetic latitude. The equatorward boundary of the auroral oval is estimated to be 33.0{\deg} invariant latitude (ILAT), which is comparable to the storms on 25/26 September 1909 (~31.6{\deg} ILAT, minimum Dst of -595 nT), 28/29 August 1859 (~36.5{\deg} ILAT), and 13/14 March 1989 (~40{\deg} ILAT, minimum Dst of -589 nT). Assuming that the equatorward boundary is a proxy for the scale of magnetic storms, we presume that the storm on March 1582 was severe. We also found that the storm on March 1582 lasted, at least, for three days by combining European records. The auroral oval stayed at mid-latitude for the first two days and moved to low-latitude (in East Asia) for the last day. It is plausible that the storm was caused by a series of ICMEs (interplanetary coronal mass ejections). We can reasonably speculate that a first ICME could have cleaned up interplanetary space to make the following ICMEs more geo-effective, as probably occurred in the Carrington and Halloween storms.

Day 1569

Wednesday.  Thursday.  Friday.

1905.05184
A Neutrino beacon
Jackson

Observational SETI has concentrated on using electromagnetism as the carrier , namely radio waves and laser radiation. Michael Hippke [2] has pointed out that it may be possible to use neutrinos or gravitational waves as signals. Gravitational waves demand the command of the generation of very large scale amounts of energy, Jackson and Benford [3]. This paper describes a beacon that uses beamed neutrinos as the signal. Neutrinos, like gravitational waves, have the advantage of extremely low extinction in the interstellar medium. To make use of neutrinos an advanced civilization can use a gravitational lens as a focus and amplifier. The lens can be a neutron star or a black hole. Using wave optics one can calculate the advantage of gravitational lensing for amplification of a beam and along the optical axis it is exceptionally large. Even though the amplification is very large the dimeter of the beam is quite small, less that a centimeter. This implies that a large constellation of neutrino transmitters would have to enclose the local neutron star or black hole to cover the sky. This means that such a beacon would have to be built by a Kardashev Type II civilization.

1905.05199

Origins of scaling relations of globular cluster systems

Choksi, Gnedin

Globular cluster (GC) systems demonstrate tight scaling relations with the properties of their host galaxies. In previous work, we developed an analytic model for GC formation in a cosmological context and showed that it matches nearly all of the observed scaling relations across 4 orders of magnitude in host galaxy mass. Motivated by the success of this model, we investigate in detail the physical origins and evolution of these scaling relations. The ratio of the combined mass in GCs $M_{\rm GC}$ to the host dark matter halo mass $M_h$ is nearly constant at all redshifts, but its normalization evolves by a factor of $\sim$10 from birth to $z=0$. The relation is steeper than linear at halo masses $M_h \lesssim 10^{11.5} M_{\odot}$, primarily due to non-linearity in the stellar mass-halo mass relation. The near constancy of the ratio $M_{\rm GC}/M_h$, combined with the shape of the stellar mass-halo mass relation, sets the characteristic $U-$shape of the GC specific frequency as a function of host galaxy mass. The contribution of accreted satellite galaxies to the buildup of GC systems is a strong function of the host galaxy mass, ranging from $\approx$0% at $M_h \approx 10^{11} M_{\odot}$ to 80% at $M_h \approx 10^{15} M_{\odot}$. The metal-poor clusters are significantly more likely to form ex-situ relative to the metal-rich clusters, but a substantial fraction of metal-poor clusters still form in-situ in lower mass galaxies. Similarly, the fraction of red clusters increases from $\approx 10$% at $M_h = 10^{11} M_{\odot}$ to $\approx 60$% at $M_h \approx 10^{13} M_{\odot}$, and flattens at higher $M_h$. Clusters formation occurs essentially continuously at high redshift, while at low redshift galactic mergers become increasingly important for cluster formation.

1905.05202

Newly discovered bright z~9-10 galaxies and improved constraints on their prevalence using the full CANDELS area

Bouwens, et al

We report the results of an expanded search for z~9-10 candidates over the ~883 arcmin^2 CANDELS+ERS fields. This study adds 147 arcmin^2 to the search area we consider over the CANDELS COSMOS, UDS, and EGS fields, while expanding our selection to include sources with bluer J_{125}-H_{160} colors than our previous J_{125}-H_{160}>0.5 mag selection. In searching for new z~9-10 candidates, we make full use of all available HST, Spitzer/IRAC, and ground-based imaging data. As a result of our expanded search and use of broader color criteria, 3 new candidate z~9-10 galaxies are identified. We also find again the z=8.683 source previously confirmed by Zitrin+2015. This brings our sample of probable z~9-11 galaxy candidates over the CANDELS+ERS fields to 19 sources in total, equivalent to 1 candidate per 47 arcmin^2 (1 per 10 WFC3/IR fields). To be comprehensive, we also discuss 28 mostly lower likelihood z~9-10 candidates, including some sources that seem to be reliably at z>8 using the HST+IRAC data alone, but which the ground-based data show are much more likely at z<4. One case example is a bright z~9.4 candidate COS910-8 which seems instead to be at z~2. Based on this expanded sample, we obtain a more robust LF at z~9 and improved constraints on the volume density of bright z~9 and z~10 galaxies. Our improved z~9-10 results again reinforce previous findings for strong evolution in the UV LF at z>8, with a factor of ~10 evolution seen in the luminosity density from z~10 to z~8.

1905.05409

Probing the assembly of dwarf galaxies through cosmic time with damped Lyman-$\alpha$ absorption spectroscopy

Jeon, Besla, Bromm

We investigate the absorption features associated with a gas-rich dwarf galaxy using cosmological hydrodynamics simulations. Our goal is to explore whether the progenitors of the lowest mass dwarf galaxies known to harbor neutral hydrogen today (M_star~10^6 solar mass, M_halo=4x10^9 solar mass) could possibly be detected as Damped Lyman-alpha Absorbers (DLAs) over cosmic time. We trace the evolution of a single dwarf galaxy, pre-selected to contain DLAs, from the era of the first metal-free, so-called Population~III (Pop~III), stars, down to z=0, thus allowing us to study the metal enrichment history of DLAs associated with the simulated galaxy. We find that the progenitors of the simulated dwarf are expected to be seen for most of their evolution as DLAs that are contaminated by normal, Population~II, stars. The time period during which DLAs are only metal-enriched by Pop~III stars, on the other hand, is likely very brief, confined to high redshifts, z~6. The susceptibility of the dwarfs to the external UV radiation background allows them to preserve neutral gas only at the centre (a few ~100 pc). This results in a small probability that the simulated dwarf would be observed as a DLA. This study suggests that DLAs are unlikely to be hosted in the lowest mass dwarfs that can harbor neutral gas (M_halo~ 4x10^9 solar mass), below which neutral gas is unlikely to exist. However, this study does illustrate that, when detected, absorption lines provide a powerful method for probing ISM conditions inside the smallest dwarf galaxies at intermediate to high redshifts.

1905.05517

Star-forming rings in lenticular galaxies: origin of the gas

Proshina, et al

Rings in S0s are enigmatic features which can however betray the evolutionary paths of particular galaxies. We have undertaken long-slit spectroscopy of five lenticular galaxies with UV-bright outer rings. The observations have been made with the Southern African Large Telescope (SALT) to reveal the kinematics, chemistry, and the ages of the stellar populations and the gas characteristics in the rings and surrounding disks. Four of the five rings are also bright in the H-alpha emission line, and the spectra of the gaseous rings extracted around the maxima of the H-alpha equivalent width reveal excitation by young stars betraying current star formation in the rings. The integrated level of this star formation is 0.1-0.2 solar mass per year, with the outstanding value of 1 solar mass per year in NGC 7808. The difference of chemical composition between the ionized gas of the rings which demonstrate nearly solar metallicity and the underlying stellar disks which are metal-poor implies recent accretion of the gas and star formation ignition; the star formation history estimated by using different star formation indicators implies that the star formation rate decreases with e-folding time of less than 1 Gyr. In NGC 809 where the UV-ring is well visible but the H-alpha emission line excited by massive stars is absent, the star formation has already ceased.

1905.05645
Radiation-pressure-driven dust transport ot galaxy haloes at $z\sim 10$
Hirashita, Inoue

The origin of dust in galaxy halos or in the circum-galactic medium (CGM) is still a mystery. We investigate if the radiation pressure in high-redshift ($z\sim 10$) galaxies can efficiently transport dust to halos. To clarify the first dust enrichment of galaxy halos in the early Universe, we solve the motion of a dust grain considering radiation pressure, gas drag, and gravity in the vertical direction of the galactic disc. Radiation pressure is estimated in a consistent manner with the stellar spectra and dust extinction. As a consequence, we find that dust grains with radii $a\sim 0.1~\mu$m successfully escape from the galactic disc if the ongoing star formation episode converts more than 15 per cent of the baryon content into stars and lasts $\gtrsim 30$ Myr, while larger and smaller grains are trapped in the disc because of gravity and gas drag, respectively. We also show that grain charge significantly enhances gas drag at a few--10 scale heights of the galactic disc, where the grain velocities are suppressed to $\sim 1$ km s$^{-1}$. There is an optimum dust-to-gas ratio ($\sim 10^{-3}$) in the galactic disc and an optimum virial mass $\sim 10^{10}$--$10^{11}$ M$_{\odot}$ for the transport of $a\sim 0.1~\mu$m grains to the halo. We conclude that early dust enrichment of galaxy halos at $z\gtrsim 10$ is important for the origin of dust in the CGM.

1905.05719

Kepler-62f: Kepler's first small planet in the habitable zone, but is it real?

Borucki, et al

Kepler-62f is the first exoplanet small enough to plausibly have a rocky composition orbiting within the habitable zone (HZ) discovered by the Kepler Mission. The planet is 1.4 times the size of the Earth and has an orbital period of 267 days. At the time of its discovery, it had the longest period of any small planet in the habitable zone of a multi-planet system. Because of its long period, only four transits were observed during Kepler's interval of observations. It was initially missed by the Kepler pipeline, but the first three transits were identified by an independent search by Eric Agol, and it was identified as a planet candidate in subsequent Kepler catalogs. However in the latest catalog of exoplanets (Thompson et al., 2018), it is labeled as a false positive. Recent exoplanet catalogues have evolved from subjective classification to automatic classifications of planet candidates by algorithms (such as `Robovetter'). While exceptionally useful for producing a uniform catalogue, these algorithms sometimes misclassify planet candidates as a false positive, as is the case of Kepler-62f. In particularly valuable cases, i.e., when a small planet has been found orbiting in the habitable zone (HZ), it is important to conduct comprehensive analyses of the data and classification protocols to provide the best estimate of the true status of the detection. In this paper we conduct such analyses and show that Kepler-62f is a true planet and not a false positive. The table of stellar and planet properties has been updated based on GAIA results.

1905.05766

Predicting the density profiles of the first haloes

Delos, et al

The first dark matter halos form by direct collapse from peaks in the matter density field, and evidence from numerical simulations and other analyses suggests that the dense inner regions of these objects largely persist today. These halos would be the densest dark matter structures in the Universe, and their abundance can probe processes that leave imprints on the primordial density field, such as inflation or an early matter-dominated era. They can also probe dark matter through its free-streaming scale. The first halos are qualitatively different from halos that form by hierarchical clustering, as evidenced by their $\rho\propto r^{-3/2}$ inner density profiles. In this work, we present and tune models that predict the density profiles of these halos from properties of the density peaks from which they collapsed. These models predict the coefficient $A$ of the $\rho=Ar^{-3/2}$ small-radius asymptote of the density profile along with the maximum circular velocity $v_\mathrm{max}$ and associated radius $r_\mathrm{max}$. These models are universal: they can be applied to any cosmology, and we confirm this by validating the models using six $N$-body simulations carried out in wildly disparate cosmological scenarios. With their connection to the primordial density field established, the first dark matter halos will serve as probes of the early Universe and the nature of dark matter.

1905.05821

A chromaticity analysis and PSF subtraction techniques for SCExAO/CHARIS data

Gerard, et al

We present an analysis of instrument performance using new observations taken with the Coronagraphic High Angular Resolution Imaging Spectrograph (CHARIS) instrument and the Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) system. In a correlation analysis of our datasets (which use the broadband mode covering J through K band in a single spectrum), we find that chromaticity in the SCExAO/CHARIS system is generally worse than temporal stability. We also develop a point spread function (PSF) subtraction pipeline optimized for the CHARIS broadband mode, including a forward modelling-based exoplanet algorithmic throughput correction scheme. We then present contrast curves using this newly developed pipeline. An analogous subtraction of the same datasets using only the H band slices yields the same final contrasts as the full JHK sequences; this result is consistent with our chromaticity analysis, illustrating that PSF subtraction using spectral differential imaging (SDI) in this broadband mode is generally not more effective than SDI in the individual J, H, or K bands. In the future, the data processing framework and analysis developed in this paper will be important to consider for additional SCExAO/CHARIS broadband observations and other ExAO instruments which plan to implement a similar integral field spectrograph broadband mode.

1905.05826

The fundamental metallicity relation emerges from the local anti-correlation between star formation rate and gas-phase metallicity existing ind disk galaxies

Almeida, Sanchez-Menguiano

The fundamental metallicity relation (FMR) states that galaxies of the same stellar mass but larger star formation rate (SFR) tend to have smaller gas-phase metallicity (<Zg>). It is thought to be fundamental because it naturally arises from the stochastic feeding of star-formation from external metal-poor gas accretion, a process extremely elusive to observe but essential according the cosmological simulations of galaxy formation. In this letter, we show how the FMR emerges from the local anti-correlation between SFR surface density and Zg recently observed to exist in disk galaxies. We analytically derive the global FMR from the local law, and then show that both relations agree quantitatively when considering the star-forming galaxies of the MaNGA survey. Thus, understanding the FMR becomes equivalent to understanding the origin of the anti-correlation between SFR and metallicity followed by the set of star-forming regions of any typical galaxy. The correspondence between local and global laws is not specific of the FMR, so that a number of local relations should exist associated with known global relations.

1905.06066

Revisiting the wave optics effect on primordial black hole constraints from optic micro lensing search

Sugiyama, Kurita, Takada

Microlensing of stars, e.g. in the Galactic bulge and Andromeda galaxy (M31), is among the most robust, powerful method to constrain primordial black holes (PBHs) that are a viable candidate of dark matter. If PBHs are in the mass range $M_{\rm PBH} \lower.5ex\hbox{$\; \buildrel > \over \sim \;$} 10^{-10}M_\odot$, its Schwarzschild radius ($r_{\rm Sch}$) becomes comparable with or shorter than optical wavelength ($\lambda)$ used in a microlensing search, and in this regime the wave optics effect on microlensing needs to be taken into account. For a lensing PBH with mass satisfying $r_{\rm Sch}\sim \lambda$, it causes a characteristic oscillatory feature in the microlensing light curve, and it will gives a smoking gun evidence of PBH if detected, because any astrophysical object cannot have such a tiny Schwarzshild radius. Even in a statistical study, e.g. constraining the abundance of PBHs from a systematic search of microlensing events for a sample of many source stars, the wave effect needs to be taken into account. We examined the impact of wave effect on the PBH constraints obtained from the $r$-band (6210\AA) monitoring observation of M31 stars in Niikura et al. (2019), and found that a finite source size effect is dominant over the wave effect for PBHs in the mass range $M_{\rm PBH}\simeq[10^{-11},10^{-10}]M_\odot$. We also discuss that, if a denser-cadence (10~sec), $g$-band monitoring observation for a sample of white dwarfs over a year timescale is available, it would allow one to explore the wave optics effect on microlensing light curve, if it occurs, or improve the PBH constraints in $M_{\rm PBH}\lower.5ex\hbox{$\; \buildrel > \over \sim \;$} 10^{-11}M_\odot$ even from a null detection.

1905.06082

Modelling baryonic feedback for survey cosmology

Chisari, et al

Observational cosmology in the next decade will rely on probes of the distribution of matter in the redshift range between $0<z<3$ to elucidate the nature of dark matter and dark energy. In this redshift range, galaxy formation is known to have a significant impact on observables such as two-point correlations of galaxy shapes and positions, altering their amplitude and scale dependence beyond the expected statistical uncertainty of upcoming experiments at separations under 10 Mpc. Successful extraction of information in such a regime thus requires, at the very least, unbiased models for the impact of galaxy formation on the matter distribution, and can benefit from complementary observational priors. This work reviews the current state of the art in the modelling of baryons for cosmology, from numerical methods to approximate analytical prescriptions, and makes recommendations for studies in the next decade, including a discussion of potential probe combinations that can help constrain the role of baryons in cosmological studies. We focus, in particular, on the modelling of the matter power spectrum, $P(k,z)$, as a function of scale and redshift, and of the observables derived from this quantity. This work is the result of a workshop held at the University of Oxford in November of 2018.

1905.06337

How do galaxies trace a large scale structure?: a case study around a massive protocluster at $z=3.13$

Shi, et al

In the hierarchical theory of galaxy formation, a galaxy overdensity is a hallmark of a massive cosmic structure. However, it is less well understood how different types of galaxies trace the underlying large-scale structure. Motivated by the discovery of a z=3.13 protocluster, we examine how the same structure is populated by Ly$\alpha$-emitting galaxies (LAEs). To this end, we have undertaken a deep narrow-band imaging survey sampling Ly$\alpha$ emission at this redshift. Of the 93 LAE candidates within a 36'x36'~(70x70~Mpc^2) field, 21 galaxies form a significant surface overdensity (delta_g=3.3+/-0.9), which is spatially segregated from the Lyman break galaxy (LBG) overdensity. One possible interpretation is that they trace two separate structures of comparable masses (~ 10^{15}M_sun) where the latter is hosted by a halo assembled at an earlier time. We speculate that the dearth of LAEs in the LBG overdensity region may signal the role of halo assembly bias in galaxy formation, which would suggest that different search techniques may be biased accordingly to the formation age or dynamical state of the host halo. The median Ly$\alpha$- and UV luminosity is 30--70\% higher for the protocluster LAEs relative to the field. This difference cannot be explained by the galaxy overdensity alone, and may require a top-heavy mass function, higher star formation efficiency for protocluster halos, or suppression of galaxy formation in low-mass halos. A luminous Ly$\alpha$ blob and an ultramassive galaxy found in this region paint a picture consistent with the expected early growth of galaxies in clusters.

1905.06341

On the formation of density filaments in the turbulent interstellar medium

Xu, Ji, Lazarian

This study is motivated by recent observations on ubiquitous interstellar density filaments and guided by modern theories of compressible magnetohydrodynamic (MHD) turbulence. The interstellar turbulence shapes the observed density structures. As the fundamental dynamics of compressible MHD turbulence, perpendicular turbulent mixing of density fluctuations entails elongated density structures aligned with the local magnetic field, accounting for low-density parallel filaments seen in diffuse atomic and molecular gas. The elongation of low-density parallel filaments depends on the turbulence anisotropy. When taking into account the partial ionization, we find that the minimum width of parallel filaments in the cold neutral medium and molecular clouds is determined by the neutral-ion decoupling scale perpendicular to magnetic field. In highly supersonic MHD turbulence in molecular clouds, both low-density parallel filaments due to anisotropic turbulent mixing and high-density filaments due to shock compression exist.

1905.06454

Cosmic shear covariance matrix in $w$CDM: cosmology matters

Harnois-Deraps, Giblin, Joachimi

We present here the cosmo-SLICS, a new suite of simulations specially designed for the analysis of current and upcoming weak lensing data beyond the standard two-point cosmic shear. We sample the $[\Omega_{\rm m}, \sigma_8, h, w_0]$ parameter space at 25 points organised in a Latin hyper-cube, spanning a range that contains most of the $2\sigma$ posterior distribution from ongoing lensing surveys. At each of these nodes we evolve a pair of $N$-body simulations in which the sampling variance is highly suppressed, and ray-trace the volumes 800 times to further increase the effective sky coverage. We extract a lensing covariance matrix from these pseudo-independent light-cones and show that it closely matches a brute-force construction based on an ensemble of 800 truly independent $N$-body runs. More precisely, a Fisher analysis reveals that both methods yield marginalized two-dimensional constraints that vary by less than 6% in area, a result that holds under different survey specifications and that matches to within 15% the area obtained from an analytical covariance calculation. Extending this comparison with our 25 $w$CDM models, we probe the cosmology dependence of the lensing covariance directly from numerical simulations, reproducing remarkably well the Fisher results from the analytical models at most cosmologies. We demonstrate that varying the cosmology at which the covariance matrix is evaluated in the first place might have an order of magnitude greater impact on the parameter constraints than varying the choice of covariance estimation technique. We present a test case in which we generate fast predictions for both the lensing signal and its associated variance with a flexible Gaussian process regression emulator, achieving an accuracy of a few percent on the former and 10% on the latter.

1905.06835

Constraining a black hole companiton for M87* though imaging by the Event Horizon Telescope

Safarzadeh, Loeb, Reid

The Event Horizon Telescope (EHT), a global very long baseline interferometric array observing at a wavelength of 1.3 mm, detected the first image of the M87 supermassive black hole (SMBH). M87 is a giant elliptical galaxy at the center of Virgo cluster, which is expected to have formed through merging of cluster galaxies. Consequently M87* hosted mergers of black holes through dynamical friction and could have one or multiple binary companions with a low mass ratio at large separations. We show that a long-term monitoring of the M87 SMBH image over $\sim$1 year with absolute positional accuracy of 1$\approx\mu as$ could detect such binary companions and exclude a large parameter space in semi major axis ($a_0$) and mass ratio ($q$), which is currently not constrained. Moreover, the presence of the accretion disk around M87* excludes a binary companion with $a_0\approx$ of order a mili parsec, as otherwise the accretion disk would have been tidally disrupted.

1905.06838

Pushing the limits of the coronagraph occult's on HST/STIS

Debes, Ren, Schneider

The Hubble Space Telescope (HST)/Space Telescope Imaging Spectrograph (STIS) contains the only currently operating coronagraph in space that is not trained on the Sun. In an era of extreme--adaptive-optics--fed coronagraphs, and with the possibility of future space-based coronagraphs, we re-evaluate the contrast performance of the STIS CCD camera. The 50CORON aperture consists of a series of occulting wedges and bars, including the recently commissioned BAR5 occulter. We discuss the latest procedures in obtaining high contrast imaging of circumstellar disks and faint point sources with STIS. For the first time, we develop a noise model for the coronagraph, including systematic noise due to speckles, which can be used to predict the performance of future coronagraphic observations. Further, we present results from a recent calibration program that demonstrates better than $10^{-6}$ point-source contrast at 0.6", ranging to $3\times10^{-5}$ point-source contrast at 0.25". These results are obtained by a combination of sub-pixel grid dithers, multiple spacecraft orientations, and post-processing techniques. Some of these same techniques will be employed by future space-based coronagraphic missions. We discuss the unique aspects of STIS coronagraphy relative to ground-based adaptive-optics--fed coronagraphs.