Tuesday, April 30, 2019

Day 1558

Tuesday.


1904.11987
Disruption of giant molecular clouds and formation of bound star clusters under they nfluence of momentum stellar feedback
Li, Vogelsberger, Marinacci, Gnedin

Energetic feedback from star clusters plays a pivotal role in shaping the dynamical evolution of giant molecular clouds (GMCs). To study the effects of stellar feedback on the star formation efficiency of the clouds and the dynamical response of embedded star clusters, we perform a suite of isolated GMC simulations with star formation and momentum feedback subgrid models using the moving-mesh hydrodynamics code \textsc{Arepo}. The properties of our simulated GMCs span a wide range of initial mass, radius, and velocity configurations. We find that the ratio of the final stellar mass to the total cloud mass, $\epsilon_{\rm int}$, scales strongly with the initial cloud surface density and momentum feedback strength. This correlation is explained by an analytic model that considers force balancing between gravity and momentum feedback. For all simulated GMCs, the stellar density profiles are systematically steeper than that of the gas at the epochs of the peaks of star formation, suggesting a centrally-concentrated stellar distribution. We also find that star clusters are always in a sub-virial state with virial parameter $\sim0.6$ prior to gas expulsion. Both the sub-virial dynamical state and steeper stellar density profiles prevent clusters from dispersal during the gas removal phase of their evolution. The final cluster bound fraction is a continuously increasing function of $\epsilon_{\rm int}$. GMCs with star formation efficiency smaller than 0.5 are still able to form clusters with large bound fractions.


1904.12029
The origin of [CII] 158um emission toward the HII region complex S235
Anderson, et al

Although the 2P3/2-2P1/2 transition of [CII] at 158um is known to be an excellent tracer of active star formation, we still do not have a complete understanding of where within star formation regions the emission originates. Here, we use SOFIA upGREAT observations of [CII] emission toward the HII region complex Sh2-235 (S235) to better understand in detail the origin of [CII] emission. We complement these data with a fully-sampled Green Bank Telescope radio recombination line map tracing the ionized hydrogen gas. About half of the total [CII] emission associated with S235 is spatially coincident with ionized hydrogen gas, although spectroscopic analysis shows little evidence that this emission is coming from the ionized hydrogen volume. Velocity-integrated [CII] intensity is strongly correlated with WISE 12um intensity across the entire complex, indicating that both trace ultra-violet radiation fields. The 22um and radio continuum intensities are only correlated with [CII] intensity in the ionized hydrogen portion of the S235 region and the correlations between the [CII] and molecular gas tracers are poor across the region. We find similar results for emission averaged over a sample of external galaxies, although the strength of the correlations is weaker. Therefore, although many tracers are correlated with the strength of [CII] emission, only WISE 12um emission is correlated on small-scales of the individual HII region S235 and also has a decent correlation at the scale of entire galaxies. Future studies of a larger sample of Galactic HII regions would help to determine whether these results are truly representative.


1904.12071
Forecasting super-sample covariance in future weak lensing surveys with SuperSCRAM
Digman, McEwen, Hirata

The observable universe contains density perturbations on scales larger than any finite volume survey. Perturbations on scales larger than a survey can measure degrade its power to constrain cosmological parameters. The dependence of survey observables such as the weak lensing power spectrum on these long-wavelength modes results in super-sample covariance. Accurately forecasting parameter constraints for future surveys requires accurately accounting for the super-sample effects. If super-sample covariance is in fact a major component of the survey error budget, it may be necessary to investigate mitigation strategies that constrain the specific realization of the long-wavelength modes. We present a Fisher matrix based formalism for approximating the magnitude of super-sample covariance and the effectiveness of mitigation strategies for realistic survey geometries. We implement our formalism in the public code SuperSCRAM: Super-Sample Covariance Reduction and Mitigation. We illustrate SuperSCRAM with an example application, where the modes contributing to super-sample covariance in the WFIRST weak lensing survey are constrained by the low-redshift galaxy number counts in the wider LSST footprint. We find that super-sample covariance increases the volume of the error ellipsoid in 7D cosmological parameter space by a factor of 4.5 relative to Gaussian statistical errors only, but our simple mitigation strategy more than halves the contamination, to a factor of 2.0.


1904.12186
Three years of Sun-as-a-star radial-velocity observations on the approach to solar minimum
Cameron, et al

The time-variable velocity fields of solar-type stars limit the precision of radial-velocity determinations of their planets' masses, obstructing detection of Earth twins. Since 2015 July we have been monitoring disc-integrated sunlight in daytime using a purpose-built solar telescope and fibre feed to the HARPS-N stellar radial-velocity spectrometer. We present and analyse the solar radial-velocity measurements and cross-correlation function (CCF) parameters obtained in the first 3 years of observation, interpreting them in the context of spatially-resolved solar observations. We describe a Bayesian mixture-model approach to automated data-quality monitoring. We provide dynamical and daily differential-extinction corrections to place the radial velocities in the heliocentric reference frakme, and the CCF shape parameters in the sidereal frame. We achieve a photon-noise limited radial-velocity precision better than 0.43 m s$^{-1}$ per 5-minute observation. The day-to-day precision is limited by zero-point calibration uncertainty with an RMS scatter of about 0.4 m s$^{-1}$. We find significant signals from granulation and solar activity. Within a day, granulation noise dominates, with an amplitude of about 0.4 m s$^{-1}$ and an autocorrelation half-life of 15 minutes. On longer timescales, activity dominates. Sunspot groups broaden the CCF as they cross the solar disc. Facular regions temporarily reduce the intrinsic asymmetry of the CCF. The radial-velocity increase that accompanies an active-region passage has a typical amplitude of 5 m s$^{-1}$ and is correlated with the line asymmetry, but leads it by 3 days. Spectral line-shape variability thus shows promise as a proxy for recovering the true radial velocity.

Monday, April 29, 2019

Day 1557

Monday.


1904.11494
Simulating the diverse instabilities of dust in magnetized gas
Hopkins, et al

Recently Squire & Hopkins showed that charged dust grains moving through magnetized gas under the influence of any external force (e.g. radiation pressure, gravity) are subject to a spectrum of instabilities. Qualitatively distinct instability families are associated with different Alfvenic or magnetosonic waves and drift or gyro motion. We present a suite of simulations exploring these instabilities, for grains in a homogeneous medium subject to an external acceleration. We vary parameters such as the ratio of Lorentz-to-drag forces on dust, plasma $\beta$, size scale, and acceleration. All regimes studied drive turbulent motions and dust-to-gas fluctuations in the saturated state, can rapidly amplify magnetic fields into equipartition with velocity fluctuations, and produce instabilities that persist indefinitely (despite random grain motions). Different parameters produce diverse morphologies and qualitatively different features in dust, but the saturated gas state can be broadly characterized as anisotropic magnetosonic or Alfvenic turbulence. Quasi-linear theory can qualitatively predict the gas turbulent properties. Turbulence grows from small to large scales, and larger-scale modes usually drive more vigorous gas turbulence, but dust velocity and density fluctuations are more complicated. In many regimes, dust forms structures (clumps, filaments, sheets) that reach extreme over-densities (up to $\gg 10^{9}$ times mean), and exhibit substantial sub-structure even in nearly-incompressible gas. These can be even more prominent at lower dust-to-gas ratios. In other regimes, dust self-excites scattering via magnetic fluctuations that isotropize and amplify dust velocities, producing fast, diffusive dust motions.


1904.11500
North--South asymmetry in solar activity and solar cycle prediction, IV: prediction for lengths of upcoming solar cycles
Javaraiah

We analyzed the daily sunspot-group data reported by the Greenwich Photoheliographic Results (GPR) during the period 1874-1976 and Debrecen Photoheliographic Data (DPD) during the period 1977-2017 and studied North-South asymmetry in the maxima and minima of the Solar Cycles 12-24. We derived the time series of the 13-month smoothed monthly mean corrected whole-spot areas of the sunspot groups in the Sun's whole sphere (WSGA), northern hemisphere (NSGA), and southern hemisphere (SSGA). From these smoothed time series we obtained the values of the maxima and minima, and the corresponding epochs, of the WSGA, NSGA, and SSGA Cycles 12-24. We find that there exists a 44-66 years periodicity in the North-South asymmetry of the minimum. A long periodicity (130-140 years) may exist in the asymmetry of the maximum. A statistically significant correlation exists between the maximum of SSGA Cycle n and the rise time of WSGA Cycle n + 2. A reasonably significant correlation also exists between the maximum of WSGA Cycle n and the decline time of WSGA Cycle n + 2. These relations suggest that the solar dynamo carries memory over at least three solar cycles. Using these relations we obtained the values 11.7 + or - 0.15 years, 11.2 + or - 0.2 years, and 11.45 + or - 0.3 years for the lengths of WSGA Cycles 24, 25, and 26, respectively, and hence, July 2020, October 2031, and March 2043 for the minimum epochs (start dates) of WSGA Cycles 25, 26, and 27, respectively. We also obtained May 2025 and March 2036 for the maximum epochs of WSGA Cycles 25 and 26, respectively. It seems during the late Maunder minimum sunspot activity was absent around the epochs of the maxima of the NSGA-cycles and the minima of the SSGA-cycles, and some activity was present at the epochs of the maxima of some SSGA-cycles and the minima of some NSGA-cycles.


1904.11508
A profile in FIRE: resolving the radial distributions of satellite galaxies in the Local Group with simulations
Samuel, Wetzel, et al

While many tensions between Local Group (LG) satellite galaxies and LCDM cosmology have been alleviated through recent cosmological simulations, the spatial distribution of satellites remains an important test of physical models and physical versus numerical disruption in simulations. Using the FIRE-2 cosmological zoom-in baryonic simulations, we examine the radial distributions of satellites with Mstar > 10^5 Msun around 8 isolated Milky Way- (MW) mass host galaxies and 4 hosts in LG-like pairs. We demonstrate that these simulations resolve the survival and physical destruction of satellites with Mstar >~ 10^5 Msun. The simulations broadly agree with LG observations, spanning the radial profiles around the MW and M31. This agreement does not depend strongly on satellite mass, even at distances <~ 100 kpc. Host-to-host variation dominates the scatter in satellite counts within 300 kpc of the hosts, while time variation dominates scatter within 50 kpc. More massive host galaxies within our sample have fewer satellites at small distances, because of enhanced tidal destruction of satellites via the baryonic disks of host galaxies. Furthermore, we quantify and provide fits to the tidal depletion of subhalos in baryonic relative to dark matter-only simulations as a function of distance. Our simulated profiles imply observational incompleteness in the LG even at Mstar >~ 10^5 Msun: we predict 2-10 such satellites to be discovered around the MW and possibly 6-9 around M31. To provide cosmological context, we compare our results with the radial profiles of satellites around MW analogs in the SAGA survey, finding that our simulations are broadly consistent with most SAGA systems.


1904.11716
The role of N2 as a geo-biosignature for the detection and characterization of Earth-like habitats
Lammer, et al

Since the Archean, N2 has been a major atmospheric constituent in Earth's atmosphere. Nitrogen is an essential element in the building blocks of life, therefore the geobiological nitrogen cycle is a fundamental factor in the long term evolution of both Earth and Earth-like exoplanets. We discuss the development of the Earth's N2 atmosphere since the planet's formation and its relation with the geobiological cycle. Then we suggest atmospheric evolution scenarios and their possible interaction with life forms: firstly, for a stagnant-lid anoxic world, secondly for a tectonically active anoxic world, and thirdly for an oxidized tectonically active world. Furthermore, we discuss a possible demise of present Earth's biosphere and its effects on the atmosphere. Since life forms are the most efficient means for recycling deposited nitrogen back into the atmosphere nowadays, they sustain its surface partial pressure at high levels. Also, the simultaneous presence of significant N2 and O2 is chemically incompatible in an atmosphere over geological timescales. Thus, we argue that an N2-dominated atmosphere in combination with O2 on Earth-like planets within circumstellar habitable zones can be considered as a geo-biosignature. Terrestrial planets with such atmospheres will have an operating tectonic regime connected with an aerobe biosphere, whereas other scenarios in most cases end up with a CO2-dominated atmosphere. We conclude with implications for the search for life on Earth-like exoplanets inside the habitable zones of M to K-stars.


1904.11725
Nitrogen fractionation in external galaxies
Viti, et al

In star forming regions in our own Galaxy, the 14N/15N ratio is found to vary from $\sim$ 100 in meteorites, comets and protoplanetary disks up to $\sim$ 1000 in pre-stellar and star forming cores, while in external galaxies the very few single-dish large scale measurements of this ratio lead to values of 100-450. The extent of the contribution of isotopic fractionation to these variations is, to date, unknown. In this paper we present a theoretical chemical study of nitrogen fractionation in external galaxies in order to determine the physical conditions that may lead to a spread of the 14N/15N ratio from the solar value of $\sim$440 and hence evaluate the contribution of chemical reactions in the ISM to nitrogen fractionation. We find that the main cause of ISM enrichment of nitrogen fractionation is high gas densities, aided by high fluxes of cosmic rays.


1904.11831
ASIME 2018 White Paper.  In-space utilization of asteroids: asteroid composition -- answers to questions from the asteroid miners
Graps, et al

In keeping with the Luxembourg government's initiative to support the future use of space resources, ASIME 2018 was held in Belval, Luxembourg on April 16-17, 2018. The goal of ASIME 2018: Asteroid Intersections with Mine Engineering, was to focus on asteroid composition for advancing the asteroid in-space resource utilisation domain. What do we know about asteroid composition from remote-sensing observations? What are the potential caveats in the interpretation of Earth-based spectral observations? What are the next steps to improve our knowledge on asteroid composition by means of ground-based and space-based observations and asteroid rendez-vous and sample return missions? How can asteroid mining companies use this knowledge? ASIME 2018 was a two-day workshop of almost 70 scientists and engineers in the context of the engineering needs of space missions with in-space asteroid utilisation. The 21 Questions from the asteroid mining companies were sorted into the four asteroid science themes: 1) Potential Targets, 2) Asteroid-Meteorite Links, 3) In-Situ Measurements and 4) Laboratory Measurements. The Answers to those Questions were provided by the scientists with their conference presentations and collected by A. Graps or edited directly into an open-access collaborative Google document or inserted by A. Graps using additional reference materials. During the ASIME 2018, first day and second day Wrap-Ups, the answers to the questions were discussed further. New readers to the asteroid mining topic may find the Conversation boxes and the Mission Design discussions especially interesting.


Friday, April 26, 2019

Day 1556

Friday.



1904.10965
Astrometric requirements for strong lensing time-delay cosmography
BIrrer, Treu

The time delay between the arrival of photons of multiple images of time variable sources can be used to constrain absolute distances in the Universe (Refsdal 1964), and in turn obtain a direct estimate of the Hubble constant and other cosmological parameters. To convert the time delay into distances, it is well known that the gravitational potential of the main deflector and the contribution of the matter along the line-of-sight need to be known to a sufficient level of precision. In this paper, we discuss a new astrometric requirement that is becoming important as time-delay cosmography improves in precision and accuracy with larger samples, and better data and modelling techniques. We derive an analytic expression for the propagation of astrometric uncertainties on the multiple image positions into the inference of the Hubble constant and derive requirements depending on image separation and relative time delay. We note that this requirement applies equally to the image position measurements and to the accuracy of the model in reproducing them. To illustrate the requirement, we discuss some example lensing configurations and highlight that, especially for time delays of order 10 days or shorter, the relative astrometric requirement is of order milli-arcseconds. With current optical infrared technology, astrometric uncertainties may be the dominant limitation for strong lensing cosmography in the small image-separation regime.


1904.10972
The X-ray halo scaling relations of supermassive black holes
Gaspari, et al

We study the correlations between (direct) masses of supermassive black holes (SMBHs) and X-ray hot halo properties, by using a Bayesian analysis of archival datasets and theoretical models. We analyze fundamental and composite X-ray variables (plasma temperature, luminosity, density, pressure, and gas/total masses) from galactic to cluster scales. We show novel key scalings, with the tightest relation being the $M_\bullet - T_{\rm x}$, followed by $M_\bullet - L_{\rm x}$ (scatter 0.2-0.3 dex). The tighter scatter and larger correlation coefficient of the X-ray halo scalings compared with the optical counterparts (including the $M_\bullet-\sigma_\ast$), together with the multivariate analysis, suggest that the plasma atmospheres play a more central role than the stellar component in the growth of SMBHs (and ultramassive BHs), in particular accounting for the group/cluster core halo. The derived gas mass scalings also correlate better with $M_\bullet$ than dark matter mass. We provide key insights on environmental features, relic galaxies, and coronae. The comparison of the optical and X-ray fundamental planes shows that, while stars can be described mainly via the virial theorem, X-ray halos are better described by univariate scalings with deviations from self-similar collapse due to feedback processes. We test 3 major channels for BH growth: hot gas accretion, chaotic cold accretion (CCA), and hierarchical BH mergers. Hot/Bondi-like models are ruled out by the data, showing inconsistent anti-correlation with X-ray halos and too low feeding. Cosmological simulations show that binary SMBH mergers are a sub-dominant channel over most of the cosmic time and too rare to induce a central-limit-theorem effect. The scalings are consistent with the predictions of CCA, the rain of matter condensing out of the turbulent X-ray halos, sustaining a self-regulated feedback loop throughout cosmic time.


1904.10992
Half-mass radii for ~7,000 galaxies at 1.0<z<2.5: most of the evolution in the mass-die relation is due to color gradients
Suess, et al

Radial mass-to-light ratio gradients cause the half-mass and half-light radii of galaxies to differ, potentially biasing studies that use half-light radii. Here we present the largest catalog to date of galaxy half-mass radii at z > 1: 7,006 galaxies in the CANDELS fields at 1.0 < z < 2.5. The sample includes both star-forming and quiescent galaxies with stellar masses 9.0 < log(M_* / M_\odot) < 11.5. We test three methods for calculating half-mass radii from multi-band PSF-matched HST imaging: two based on spatially-resolved SED modeling, and one that uses a rest-frame color profile. All three methods agree, with scatter <~0.3 dex. In agreement with previous studies, most galaxies in our sample have negative color gradients (the centers are redder than the outskirts, and r_e,mass < r_e,light). We find that color gradient strength has significant trends with increasing stellar mass, half-light radius, U-V color, and stellar mass surface density. These trends have not been seen before at z>1. Furthermore, color gradients of star-forming and quiescent galaxies show a similar redshift evolution: they are flat at z>~2, then steeply decrease as redshift decreases. This affects the galaxy mass-size relation. The normalizations of the star-forming and quiescent r_mass-M_* relations are 10-40% smaller than the corresponding r_light-M_* relations; the slopes are ~0.1-0.3 dex shallower. Finally, the half-mass radii of star-forming and quiescent galaxies at M_* = 10^{10.5}M_\odot only grow by ~1%$ and ~8% between z~2.25 and z~1.25. This is significantly less than the ~37% and ~47% size increases found when using the half-light radius.


1904.10994
The life cycle of dust
Sadavoy, et al

Dust offers a unique probe of the interstellar medium (ISM) across multiple size, density, and temperature scales. Dust is detected in outflows of evolved stars, star-forming molecular clouds, planet-forming disks, and even in galaxies at the dawn of the Universe. These grains also have a profound effect on various astrophysical phenomena from thermal balance and extinction in galaxies to the building blocks for planets, and changes in dust grain properties will affect all of these phenomena. A full understanding of dust in all of its forms and stages requires a multi-disciplinary investigation of the dust life cycle. Such an investigation can be achieved with a statistical study of dust properties across stellar evolution, star and planet formation, and redshift. Current and future instrumentation will enable this investigation through fast and sensitive observations in dust continuum, polarization, and spectroscopy from near-infrared to millimeter wavelengths.


1904.11014
The nature of ionized gas in the Milky Way galactic fountain
Werk, et al

We address the spatial scale, ionization structure and metal content of gas at the Milky Way disk-halo interface detected as absorption in the foreground of seven halo blue horizontal branch stars (BHBs). The closely-spaced, high-latitude BHB sightlines reach to different heights above the disk, ranging from z = 3 - 14 kpc in the direction of the Northern Galactic Pole. With our high-resolution, S/N > 14, HST/COS and Keck/HIRES spectra, we constrain the distances, velocities, sizes, and carbon content of infalling, intermediate-velocity gaseous structures (IV gas; -90 < v_LSR [km s^-1] < -25 ). We detect transitions that trace multiple ionization states (e.g. CaII, FeII, SiIV, CIV) with column densities that remain constant with height from the disk, indicating that the gas most likely lies within z < 3.4 kpc. The intermediate ionization state gas traced by CIV and SiIV is strongly correlated over the full range of transverse separations probed by our sightlines, indicating large, coherent structures greater than 500 pc in size. The low ionization state material traced by CaII and FeII neither exhibits a correlation with N_HI nor transverse separation, implying cloudlets or clumpiness on scales less than 10 pc. We find that the observed ratio log(N_SiIV/ N_CIV), with a median value of -0.69+/-0.04, is sensitive to the total carbon content of the ionized gas under the assumption of either photoionization or collisional ionization. Our data directly support a physical scenario in which the Milky Way is fed by complex, multiphase processes at its disk-halo interface that involve kpc-scale ionized envelopes or streams containing pc-scale, cool clumps


1904.11070
RascalC: a Jackknife approach to estimating single and multi-tracer galaxy covariance matrices
Philcox, et al

To make use of clustering statistics from large cosmological surveys, accurate and precise covariance matrices are needed. We present a new code to estimate large scale galaxy correlation function covariances in arbitrary survey geometries that produces results comparable to a suite of $10^6$ mocks in $\sim 100$ CPU-hours, orders-of-magnitude faster than pre-existing codes. As in previous works, non-Gaussianity is encapsulated via a shot-noise rescaling, with calibrations performed by comparing models to jackknifed survey data. The approach requires data from only a single dataset (without an input correlation function model), and the deviations between large scale model covariances from a mock survey and those from a large suite of mocks are found to be be indistinguishable from noise. In addition, the choice of input mock are shown to be irrelevant for desired noise levels below $\sim 10^5$ mocks. Coupled with its generalization to multi-tracer data-sets, this shows the algorithm to be an excellent tool for analysis, reducing the need for large numbers of mock simulations to be computed.


1904.11185
Dust production scenarios in galaxies at z~6-8.3
Lesniewska, Michalowski

The mechanism of dust formation in galaxies at high redshift is still unknown. Asymptotic giant branch (AGB) stars and explosions of supernovae (SNe) are possible dust producers, and non-stellar processes may substantially contribute to dust production, for example grain growth in the interstellar medium (ISM). Our aim is to determine the contribution to dust production of AGB stars and SNe in nine galaxies at z ~ 6-8.3, for which observations of dust have been recently attempted. In order to determine the origin of the observed dust we have determined dust yields per AGB star and SN required to explain the total amounts of dust in these galaxies. We find that AGB stars were not able to produce the amounts of dust observed in the galaxies in our sample. In order to explain these dust masses, SNe would have to have maximum efficiency and not destroy the dust which they formed. Therefore, the observed amounts of dust in the galaxies in the early universe were formed either by efficient supernovae or by a non-stellar mechanism, for instance the grain growth in the interstellar medium.

Thursday, April 25, 2019

Day 1555

Thursday.



1904.10459
Ionization bias and the ghost proximity effect near $z\gtrsim6$ quasars in the shadow of proximate absorption systems
Davies

The larger-than-expected scatter in the opacity of the Ly$\alpha$ forest suggests that the metagalactic ionizing background is strongly fluctuating at $z > 5.5$. Models for ionizing background fluctuations predict a strong positive bias on large scales, so the environments of massive $>10^{12}\,{\rm M}_\odot$ dark matter halos, e.g. $z\sim6$ quasar hosts, would be ideal laboratories to constrain the sources of ionizing photons. While the quasars themselves should overwhelm any plausible ionizing photon contribution from neighboring galaxies, proximate damped Ly$\alpha$ absorbers (DLAs) have recently been discovered in the foreground of $z\sim6$ quasars, and the Ly$\alpha$ forest in the shadow of these DLAs could probe the local ionization environment. Using Gpc$^3$ simulations of $z=6$ ionizing background fluctuations, we show that while the Ly$\alpha$ forest signal from ionization bias around a quasar host halo should be strong, it is likely suppressed by the associated intergalactic matter overdensity. We also show that the quasar itself may still overwhelm the clustering signal via a "ghost" of the proximity effect from the quasar radiation causing a large-scale bias in the ionizing photon mean free path. This ghost proximity effect is sensitive to the lifetime and geometry of quasar emission, potentially unlocking a new avenue for constraining these fundamental quasar properties. Finally, we present observations of a $z\sim6$ quasar with a proximate DLA which shows a strong excess in Ly$\alpha$ forest transmission at the predicted location of the ghost proximity effect.


1904.10521
Accounting for correlations when fitting extra cosmological parameters
Huang, Addison, Bennett

Current cosmological tensions motivate investigating extensions to the standard $\Lambda$CDM model. Additional model parameters are typically varied one or two at a time, in a series of separate tests. The purpose of this paper is to highlight that information is lost by not also examining the correlations between these additional parameters, which arise when their effects on model predictions are similar, even if the parameters are not varied simultaneously. We show how these correlations can be quantified with simulations and Markov Chain Monte Carlo (MCMC) methods. As an example, we assume that $\Lambda$CDM is the true underlying model, and calculate the correlations expected between the phenomenological lensing amplitude parameter, $A_L$, the running of the spectral index, $n_{\rm run}$, and the primordial helium mass fraction, $Y_P$, when these parameters are varied one at a time along with the $\Lambda$CDM parameters in fits to the $\textit{Planck}$ 2015 temperature power spectrum. These correlations are not small, ranging from 0.31 ($A_L-n_{\rm run}$) to $-0.93$ ($n_{\rm run}-Y_P$). We find that the values of these three parameters from the $\textit{Planck}$ data are consistent with $\Lambda$CDM expectations within $0.9\sigma$ when the correlations are accounted for. This does not explain the 1.8-2.7$\sigma$ $\textit{Planck}$ preference for $A_L>1$, but provides an additional $\Lambda$CDM consistency test. For example, if $A_L>1$ was a symptom of an underlying systematic error or some real but unknown physical effect that also produced spurious correlations with $n_{\rm run}$ or $Y_P$ our test might have revealed this. We recommend that future cosmological analyses examine correlations between additional model parameters in addition to investigating them separately, one a time.


1904.10526
Upper limits on very-high-energy Gamma-ray emission from Core-collapse supernovae observed with H.E.S.S
H.E.S.S. collaboration, et al

Young core-collapse supernovae with dense-wind progenitors may be able to accelerate cosmic-ray hadrons beyond the knee of the cosmic-ray spectrum, and this may result in measurable gamma-ray emission. We searched for gamma-ray emission from ten supernovae observed with the High Energy Stereoscopic System (H.E.S.S.) within a year of the supernova event. Nine supernovae were observed serendipitously in the H.E.S.S. data collected between December 2003 and December 2014, with exposure times ranging from 1.4 hours to 53 hours. In addition we observed SN 2016adj as a target of opportunity in February 2016 for 13 hours. No significant gamma-ray emission has been detected for any of the objects, and upper limits on the $>1$ TeV gamma-ray flux of the order of $\sim$10$^{-13}$ cm$^{-2}$s$^{-1}$ are established, corresponding to upper limits on the luminosities in the range $\sim$2 $\times$ 10$^{39}$ erg s$^{-1}$ to $\sim$1 $\times$ 10$^{42}$ erg s$^{-1}$. These values are used to place model-dependent constraints on the mass-loss rates of the progenitor stars, implying upper limits between $\sim$2 $\times 10^{-5}$ and $\sim$2 $\times 10^{-3}$M$_{\odot}$yr$^{-1}$ under reasonable assumptions on the particle acceleration parameters.


1904.10857
Extracting the global signal from 21-cm fluctuations: the multi-tracer approach
Fialkov, Barkana, Jarvis

The multi-tracer technique employs a ratio of densities of two differently biased galaxy samples that trace the same underlying matter density field, and was proposed to alleviate the cosmic variance problem. Here we propose a novel application of this approach, applying it to two different tracers one of which is the 21-cm signal of neutral hydrogen from the epochs of reionization and comic dawn. The second tracer is assumed to be a sample of high-redshift galaxies, but the approach can be generalized and applied to other high-redshift tracers. We show that the anisotropy of the ratio of the two density fields can be used to measure the sky-averaged 21-cm signal, probe the spectral energy distribution of radiative sources that drive this signal, and extract large-scale properties of the second tracer, e.g., the galaxy bias. Using simulated 21-cm maps and mock galaxy samples, we find that the method works well for an idealized galaxy survey. However, in the case of a realistic galaxy survey which only probes highly biased luminous galaxies, the inevitable Poisson noise makes the reconstruction far more challenging. This difficulty can be mitigated with the greater sensitivity of future telescopes along with larger survey volumes.


1904.10879
Accuracy requirements for empirically-measured selection functions
Farr

I give formulas for the accuracy to which a selection function must be measured via Monte-Carlo injections in order to have un-biased population inference. The number of found injections scales linearly with the number of objects in the population; the coefficient in front of the linear term depends on both the distribution of injections and the inferred population distribution.

Tuesday, April 23, 2019

Day 1554

Wednesday.



1904.09987
The minimum metallicity of globular clusters and its physical origin -- implications for the galaxy mass-metallicity relation and observations of proto-globular clusters at high redshift
Kruijssen

In the local Universe, globular clusters (GCs) with metallicities $[{\rm Fe}/{\rm H}]<-2.5$ are extremely rare. In this Letter, the close connection between GC formation and galaxy evolution is used to show that this GC metallicity `floor' results from the galaxy mass-metallicity relation of ultra low-luminosity galaxies (ULLGs) at high redshift, where the most metal-poor GCs must have formed. Galaxies with metallicities $[{\rm Fe}/{\rm H}]\lesssim-2.5$ have too low masses to form GCs with initial masses $M_{\rm i}\gtrsim10^5~{\rm M}_\odot$, needed to survive for a Hubble time. This translates the galaxy mass-metallicity relation into a maximum initial cluster mass-metallicity relation for $[{\rm Fe}/{\rm H}]\lesssim-1.8$, which naturally leads to the observed colour-magnitude relation of metal-poor GCs at $z=0$ (the `blue tilt'). Its strength traces the slope of the gas phase mass-metallicity relation of ULLGs. Based on the observed blue tilt of GCs in the Virgo and Fornax Clusters, the galaxy mass-metallicity relation is predicted to have a slope of $\alpha=0.4\pm0.1$ for $10^5\lesssim M_\star/{\rm M}_\odot\lesssim10^7$ at $z\gtrsim2$. The GC metallicity floor implies a minimum host galaxy mass and a maximum redshift for GC formation. Any proto-GCs that may be detected at $z>9$ are most likely to end up in galaxies presently more massive than the Milky Way, whereas GCs in low-mass galaxies such as the Fornax dSph ($M_\star\approx4\times10^7~{\rm M}_\odot$) formed at $z\lesssim3$.


1904.10140
Graphene molecule compared with fullerene C60 as circumstellar carbon dust of planetary nebula
Ota

It had been understood that astronomically observed infrared spectrum of carbon rich planetary nebula as like Tc 1 and Lin 49 comes from fullerene (C60). Also, it is well known that graphene is a raw material for synthesizing fullerene. This study seeks some capability of graphene based on the quantum-chemical DFT calculation. It was demonstrated that graphene plays major role rather than fullerene. We applied two astrophysical conditions, which are void creation by high speed proton and photo-ionization by the central star. Model molecule was ionized void-graphene (C23) having one carbon pentagon combined with hexagons. By molecular vibrational analysis, we could reproduce six major bands from 6 to 9 micrometer, large peak at 12.8, and largest peak at 19.0. Also, many minor bands could be reproduced from 6 to 38 micrometer. Also, deeply void induced molecules (C22) and (C21) could support observed bands.


1904.10220
Galactic cosmic rays after AMS-02
Evoli, et al

The unprecedented quality of the data collected by the AMS-02 experiment onboard the International Space Station allowed us to address subtle questions concerning the origin and propagation of cosmic rays. Here we discuss the implications of these data for the injection spectrum of elements with different masses and for the diffusion coefficient probed by cosmic rays through their journey from the sources to the Earth. We find that the best fit to the spectra of primary and secondary nuclei requires 1) a break in the energy dependence of the diffusion coefficient at energies $\sim 300$ GV; 2) an injection spectrum that is the same for all nuclei heavier than helium, and different injections for both protons and helium. Moreover, if to force the injection spectrum of helium to be the same as for heavier nuclei, the fit to oxygen substantially worsens. Accounting for a small, $X_{s}\sim 0.4~\rm g~cm^{-2}$, grammage accumulated inside the sources leads to a somewhat better fit to the B/C ratio but makes the difference between He and other elements even more evident. The statistic and systematic error bars claimed by the AMS collaboration exceed the error that is expected from calculations once the uncertainties in the cross sections of production of secondary nuclei are taken into account. In order to make this point more quantitative, we present a novel parametrization of a large set of cross sections, relevant for cosmic ray physics, and we introduce the uncertainty in the branching ratios in a way that its effect can be easily grasped.


1904.10275
Projection effects in galaxy cluster samples: insights from X-ray redshifts
Ramos-Ceja, et al

Up to now, the largest sample of galaxy clusters selected in X-rays comes from the ROSAT All-Sky Survey (RASS). Although there have been many interesting clusters discovered with the RASS data, the broad point spread function (PSF) of the ROSAT satellite limits the amount of spatial information of the detected objects. This leads to the discovery of new cluster features when a re-observation is performed with higher resolution X-ray satellites. Here we present the results from XMM-Newton observations of three clusters: RXCJ2306.6-1319, ZwCl1665 and RXCJ0034.6-0208, for which the observations reveal a double or triple system of extended components. These clusters belong to the extremely expanded HIghest X-ray FLUx Galaxy Cluster Sample (eeHIFLUGCS), which is a flux-limited cluster sample ($f_\textrm{X,500}\geq 5\times10^{-12}$ erg s$^{-1}$ cm$^{-2}$ in the $0.1-2.4$ keV energy band). For each structure in each cluster, we determine the redshift with the X-ray spectrum and find that the components are not part of the same cluster. This is confirmed by an optical spectroscopic analysis of the galaxy members. Therefore, the total number of clusters is actually 7 and not 3. We derive global cluster properties of each extended component. We compare the measured properties to lower-redshift group samples, and find a good agreement. Our flux measurements reveal that only one component of the ZwCl1665 cluster has a flux above the eeHIFLUGCS limit, while the other clusters will no longer be part of the sample. These examples demonstrate that cluster-cluster projections can bias X-ray cluster catalogues and that with high-resolution X-ray follow-up this bias can be corrected.


1904.10288
The heavy-element content red of planets: a tracer of their formation sites
Hasegawa, et al

Identification of the main planet formation site is fundamental to understanding how planets form and migrate to the current locations. We consider the heavy-element content trend of observed exoplanets derived from improved measurements of mass and radius, and explore how this trend can be used as a tracer of their formation sites. Using gas accretion recipes obtained from detailed hydrodynamical simulations, we confirm that the disk-limited gas accretion regime is most important for reproducing the heavy-element content trend. Given that such a regime is specified by two characteristic masses of planets, we compute these masses as a function of the distance ($r$) from the central star, and then examine how the regime appears in the mass-semimajor axis diagram. Our results show that a plausible solid accretion region emerges at $r \simeq 0.6$ au and expands with increasing $r$, using the conventional disk model. Given that exoplanets that possess the heavy-element content trend distribute currently near their central stars, our results imply the importance of planetary migration that would occur after solid accretion onto planets might be nearly completed at $r \geq 0.6$ au. Self-consistent simulations would be needed to verify the predictions herein.


1904.10438
Mini-survey of the northern sky to Dec <+30
Capak, et al

We propose an extension of the LSST survey to cover the northern sky to DEC < +30 (accessible at airmass <1.8). This survey will increase the LSST sky coverage by ~9,600 square degrees from 18,900 to 28,500 square degrees (a 50% increase) but use only 0.6-2.5% of the time depending on the synergies with other surveys. This increased area addresses a wide range of science cases that enhance all of the primary LSST science goals by significant amounts. The science enabled includes: increasing the area of the sky accessible for follow-up of multi-messenger transients including gravitational waves, mapping the milky way halo and halo dwarfs including discovery of RR Lyrae stars in the outer galactic halo, discovery of z>7 quasars in combination Euclid, enabling a second generation DESI and other spectroscopic surveys, and enhancing all areas of science by improving synergies with Euclid, WFIRST, and unique northern survey facilities. This white paper is the result of the Tri-Agency Working Group (TAG) appointed to develop synergies between missions and presents a unified plan for northern coverage. The range of time estimates reflects synergies with other surveys. If the modified DESC WFD survey, the ecliptic plane mini survey, and the north galactic spur mini survey are executed this plan would only need 0.6% of the LSST time, however if none of these are included the overall request is 2.5% of the 10 year survey life. In other words, the majority of these observations are already suggested as part of these other surveys and the intent of this white paper is to propose a unified baseline plan to carry out a broad range of objectives to facilitate a combination of multiple science objectives. A companion white paper gives Euclid specific science goals, and we support the white papers for southern extensions of the LSST survey.


1904.10439
Enhancing LSST science with Euclid synergy
Capak, et al

This white paper is the result of the Tri-Agency Working Group (TAG) appointed to develop synergies between missions and is intended to clarify what LSST observations are needed in order to maximally enhance the combined science output of LSST and Euclid. To facilitate LSST planning we provide a range of possible LSST surveys with clear metrics based on the improvement in the Dark Energy figure of merit (FOM). To provide a quantifiable metric we present five survey options using only between 0.3 and 3.8% of the LSST 10 year survey. We also provide information so that the LSST DDF cadence can possibly be matched to those of \emph{Euclid} in common deep fields, SXDS, COSMOS, CDFS, and a proposed new LSST deep field (near the Akari Deep Field South). Co-coordination of observations from the Large Synoptic Survey Telescope (LSST) and Euclid will lead to a significant number of synergies. The combination of optical multi-band imaging from LSST with high resolution optical and near-infrared photometry and spectroscopy from \emph{Euclid} will not only improve constraints on Dark Energy, but provide a wealth of science on the Milky Way, local group, local large scale structure, and even on first galaxies during the epoch of reionization. A detailed paper has been published on the Dark Energy science case (Rhodes et al.) by a joint LSST/Euclid working group as well as a white paper describing LSST/Euclid/WFIRST synergies (Jain et al.), and we will briefly describe other science cases here. A companion white paper argues the general science case for an extension of the LSST footprint to the north at airmass < 1.8, and we support the white papers for southern extensions of the LSST survey.

Day 1553

Tuesday.



1904.09326
Disentangling nature from nurture: tracing the origin of seed black holes
Natarajan, et al

The origin and properties of black hole seeds that grow to produce the detected population of supermassive black holes are unconstrained at present. Despite the existence of several potentially feasible channels for the production of initial seeds in the high redshift universe, since even actively growing seeds are not directly observable at these epochs, discriminating between models remains challenging. Several new observables that encapsulate information about seeding have been proposed in recent years, and these offer exciting prospects for truly unraveling the nature of black hole seeds in the coming years. One of the key challenges for this task lies in the complexity of the problem, the required disentangling of the confounding effects of accretion physics and mergers, as mergers and accretion events over cosmic time stand to erase these initial conditions. Nevertheless, some unique signatures of seeding do survive and still exist in: local scaling relations between black holes and their galaxy hosts at low-masses; in high-redshift luminosity functions of accreting black holes; and in the total number and mass functions of gravitational wave coalescence events from merging binary black holes. One of the clearest discriminants for seed models are these high redshift gravitational wave detections of mergers from space detectable in the milliHertz range. These predicted event rates offer the most direct constraints on the properties of initial black hole seeds. Improving our theoretical understanding of black hole dynamics and accretion will also be pivotal in constraining seeding models in combination with the wide range of multi-messenger data.


1904.09898
Mass loss via Solar wind and coronal mass ejections during solar cycle 23 and 24
Mischra, et al

Similar to the Sun, other stars shed mass and magnetic flux via ubiquitous quasi-steady wind and episodic stellar coronal mass ejections (CMEs). We investigate the mass loss rate via solar wind and CMEs as a function of solar magnetic variability represented in terms of sunspot number and solar X-ray background luminosity. We estimate the contribution of CMEs to the total solar wind mass flux in the ecliptic and beyond, and its variation over different phases of the solar activity cycles. The study exploits the number of sunspots observed, coronagraphic observations of CMEs near the Sun by SOHO/LASCO, in situ observations of the solar wind at 1 AU by WIND, and GOES X-ray flux during solar cycle 23 and 24. We note that the X-ray background luminosity, occurrence rate of CMEs and ICMEs, solar wind mass flux, and associated mass loss rates from the Sun do not decrease as strongly as the sunspot number from the maximum of solar cycle 23 to the next maximum. Our study confirms a true physical increase in CME activity relative to the sunspot number in cycle 24. We show that the CME occurrence rate and associated mass loss rate can be better predicted by X-ray background luminosity than the sunspot number. The solar wind mass loss rate which is an order of magnitude more than the CME mass loss rate shows no obvious dependency on cyclic variation in sunspot number and solar X-ray background luminosity. These results have implications to the study of solar-type stars.

Sunday, April 21, 2019

Day 1552

Monday.



1904.08929
The MASSIVE survey XIII -- spatially resolved stellar kinematics in the central 1 kpc of 20 massive elliptical galaxies with the GMOS-North Integral-field spectrograph
Ene, et al

We use observations from the GEMINI-N/GMOS integral-field spectrograph (IFS) to obtain spatially resolved stellar kinematics of the central $\sim 1$ kpc of 20 early-type galaxies (ETGs) with stellar masses greater than $10^{11.7} M_\odot$ in the MASSIVE survey. Together with observations from the wide-field Mitchell IFS at McDonald Observatory in our earlier work, we obtain unprecedentedly detailed kinematic maps of local massive ETGs, covering a scale of $\sim 0.1-30$ kpc. The high ($\sim 120$) signal-to-noise of the GMOS spectra enable us to obtain two-dimensional maps of the line-of-sight velocity, velocity dispersion $\sigma$, as well as the skewness $h_3$ and kurtosis $h_4$ of the stellar velocity distributions. All but one galaxy in the sample have $\sigma(R)$ profiles that increase towards the center, whereas the slope of $\sigma(R)$ at one effective radius ($R_e$) can be of either sign. The $h_4$ is generally positive, with 14 of the 20 galaxies having positive $h_4$ within the GMOS aperture and 18 having positive $h_4$ within $1 R_e$. The positive $h_4$ and rising $\sigma(R)$ towards small radii are indicative of a central black hole and velocity anisotropy. We demonstrate the constraining power of the data on the mass distributions in ETGs by applying Jeans anisotropic modeling (JAM) to NGC~1453, the most regular fast rotator in the sample. Despite the limitations of JAM, we obtain a clear $\chi^2$ minimum in black hole mass, stellar mass-to-light ratio, velocity anisotropy parameters, and the circular velocity of the dark matter halo.


1904.08941
The supersonic project: shining light on SIGOs - a new formation channel for globular clusters
Chiou, et al

Supersonically induced gas objects (SIGOs) with little to no dark matter component are predicted to exist in patches of the Universe with non-negligible relative velocity between baryons and the dark matter at the time of recombination. Using {\sc arepo} hydrodynamic simulations we find that the gas densities inside these objects are high enough to allow stars to form. An estimate of the luminosity of the first star clusters formed within these SIGOs suggests that they may be observed at high redshift using future HST and JWST observations. Furthermore, our simulations indicate that SIGOs lie in a distinct place in the luminosity-radius parameter space, which can be used observationally to distinguish SIGOs from dark-matter hosting gas systems. Finally, as a proof-of-concept, we model star formation before reionization and evolve these systems to current times. We find that SIGOs occupy a similar part of the magnitude-radius parameter space as globular clusters. These results suggest that SIGOs may be linked with present-day local globular clusters. Since the relative velocity between the baryons and dark matter is coherent over a few Mpc scales, we predict that if this is the dominant mechanism for the formation of globular clusters, their abundance should vary significantly over these scales.

Friday, April 19, 2019

Day 1551

Friday.



1904.08406
Revealing the differences in the SMBH accretion rate distributions of starburst and non-starburst galaxies
Grimmett, et al

We infer and compare the specific X-ray luminosity distributions for a sample of massive (i.e. $\log_{10} (M*/M\odot) > 10.5$) galaxies split according to their far-infrared-derived star-forming properties (i.e., starburst and non-starburst) and redshift. We model each distribution as a power-law with an upper and lower turnover, and adopt a maximum likelihood method to include information from non-detections in the form of upper limits. When we use our inferred distributions to calculate the ratios of high to low sLx AGN (corresponding to above and below $0.1\lambda_{\text{Edd}}$, respectively) we find that starbursts have significantly higher proportions of high sLx AGN compared to their non-starburst counterparts. These findings help explain the increase in average X-ray luminosity in bins of increasing SFR reported by previous studies.


1904.08431
Origin of the Golden Mass of galaxies and black holes
Dekel, et al

We address the origin of the golden mass and time for galaxy formation and the onset of rapid black-hole growth. The preferred dark-halo mass of ~$10^{12}M_\odot$ is translated to a characteristic epoch, z~2, at which the typical forming halos have a comparable mass. We put together a coherent picture based on existing and new simple analytic modeling and cosmological simulations. We describe how the golden mass arises from two physical mechanisms that suppress gas supply and star formation below and above the golden mass, supernova feedback and virial shock heating of the circum-galactic medium (CGM), respectively. Cosmological simulations reveal that these mechanisms are responsible for a similar favored mass for the dramatic events of gaseous compaction into compact star-forming "blue nuggets", caused by mergers, counter-rotating streams or other mechanisms. This triggers inside-out quenching of star formation, to be maintained by the hot CGM, leading to today's passive early-type galaxies. The blue-nugget phase is responsible for transitions in the galaxy structural, kinematic and compositional properties, e.g., from dark-matter to baryon central dominance and from prolate to oblate shape. The growth of the central black hole is suppressed by supernova feedback below the critical mass, and is free to grow once the halo is massive enough to lock the supernova ejecta by its deep potential well and the hot CGM. A compaction near the golden mass makes the black hole sink to the galactic center and triggers a rapid black-hole growth. This ignites feedback by the Active Galactic Nucleus that helps keeping the CGM hot and maintaining long-term quenching.


1904.08434
The Herschel Dwarf Galaxy Survey: II. Physical conditions, origin of [CII] emission, and porosity of the multiphase low-metallicity ISM
Cormier, et al

The sensitive infrared telescopes, Spitzer and Herschel, have been used to target low-metallicity star-forming galaxies, allowing us to investigate the properties of their interstellar medium (ISM) in unprecedented detail. Interpretation of the observations in physical terms relies on careful modeling of those properties. We have employed a multiphase approach to model the ISM phases (HII region and photodissociation region) with the spectral synthesis code Cloudy. Our goal is to characterize the physical conditions (gas densities, radiation fields, etc.) in the ISM of the galaxies from the Herschel Dwarf Galaxy Survey. We are particularly interested in correlations between those physical conditions and metallicity or star-formation rate. Other key issues we have addressed are the contribution of different ISM phases to the total line emission, especially of the [CII]157um line, and the characterization of the porosity of the ISM. We find that the lower-metallicity galaxies of our sample tend to have higher ionization parameters and galaxies with higher specific star-formation rates have higher gas densities. The [CII] emission arises mainly from PDRs and the contribution from the ionized gas phases is small, typically less than 30% of the observed emission. We also find correlation - though with scatter - between metallicity and both the PDR covering factor and the fraction of [CII] from the ionized gas. Overall, the low metal abundances appear to be driving most of the changes in the ISM structure and conditions of these galaxies, and not the high specific star-formation rates. These results demonstrate in a quantitative way the increase of ISM porosity at low metallicity. Such porosity may be typical of galaxies in the young Universe.


1904.08446
Mars obliquity history constrained by elliptic crater orientations
Holo, et al

The dynamics of Mars' obliquity are believed to be chaotic, and the historical ~3.5 Gyr (late-Hesperian onward) obliquity probability density function (PDF) is high uncertain and cannot be inferred from direct simulation alone. Obliquity is also a strong control on post-Noachian Martian climate, enhancing the potential for equatorial ice/snow melting and runoff at high obliquities (> 40{\deg}) and enhancing the potential for desiccation of deep aquifers at low obliquities (< 25{\deg}). We developed a new technique using the orientations of elliptical craters to constrain the true late-Hesperian-onward obliquity PDF. To do so, we developed a forward model of the effect of obliquity on elliptic crater orientations using ensembles of simulated Mars impactors and ~3.5 Gyr-long Mars obliquity simulations. In our model, the inclinations and speeds of Mars crossing objects bias the preferred orientation of elliptic craters which are formed by low-angle impacts. Comparison of our simulation predictions with a validated database of elliptic crater orientations allowed us to invert for best-fitting obliquity history. We found that since the onset of the late-Hesperian, Mars' mean obliquity was likely low, between ~10{\deg} and ~30{\deg}, and the fraction of time spent at high obliquities > 40{\deg} was likely < 20%.

Thursday, April 18, 2019

Day 1550

Thursday.



1904.07874
On the detectability of visible-wavelength line emission from the local circumgalactic and intergalactic medium
Lokhorst, et al

We describe a new approach to studying the intergalactic and circumgalactic medium in the local Universe: direct detection through narrow-band imaging of ultra-low surface brightness visible-wavelength line emission. We use the hydrodynamical cosmological simulation EAGLE to investigate the expected brightness of this emission at low redshift ($z$ $\lesssim$ 0.2). H$\alpha$ emission in extended halos (analogous to the extended Ly$\alpha$ halos/blobs detected around galaxies at high redshifts) has a surface brightness of $\gtrsim700$ photons cm$^{-2}$ sr$^{-1}$ s$^{-1}$ out to $\sim$100 kpc. Mock observations show that the Dragonfly Telephoto Array, equipped with state-of-the-art narrow-band filters, could directly image these structures in exposure times of $\sim$10 hours. H$\alpha$ fluorescence emission from this gas can be used to place strong constraints on the local ultra-violet background, and on gas flows around galaxies. Detecting H$\alpha$ emission from the diffuse intergalactic medium (the "cosmic web") is beyond current capabilities, but would be possible with a hypothetical 1000-lens Dragonfly array.


1904.07905
Weak lensing in the Horizon-AGN simulation light cone.  Small scale baryonic effects
Gouin, et al

Context. Accurate model predictions including the physics of baryons are required to make the most of the upcoming large cosmological surveys devoted to gravitational lensing. The advent of hydrodynamical cosmological simulations enables such predictions on sufficiently sizeable volumes. Aims. Lensing quantities (deflection, shear, convergence) and their statistics (convergence power spectrum, shear correlation functions, galaxy-galaxy lensing) are computed in the past lightcone built in the Horizon-AGN hydrodynamical cosmological simulation, which implements our best knowledge on baryonic physics at the galaxy scale in order to mimic galaxy populations over cosmic time. Methods. Lensing quantities are generated over a one square degree field of view by performing multiple-lens plane ray-tracing through the lightcone, taking full advantage of the 1 kpc resolution and splitting the line of sight over 500 planes all the way to redshift z~7. Two methods are explored (standard projection of particles with adaptive smoothing, and integration of the acceleration field) to assert a good implementation. The focus is on small scales where baryons matter most. Results. Standard cosmic shear statistics are impacted at the 10% level by the baryonic component for angular scales below a few arcmin. The galaxy-galaxy lensing signal, or galaxy-shear correlation function, is consistent with measurements for the redshift z~0.5 massive galaxy population. At higher redshift z>1, the impact of magnification bias on this correlation is relevant for separations greater than 1 Mpc. Conclusions. This work is pivotal for all current and upcoming weak lensing surveys and represents a first step towards building a full end-to-end generation of lensed mock images from large cosmological hydrodynamical simulations.


1904.07941
Ionising the intergalactic medium by star clusters ? The first empirical evidence
Vanzella, et al

We present a VLT/X-Shooter spectroscopy of the Lyman continuum (LyC) emitting galaxy Ion2 at z=3.2121 and compare it to that of the recently discovered strongly lensed LyC-emitter at z=2.37, known as the "Sunburst" arc. Three main results emerge from the X-Shooter spectrum: (a) the Lya has three distinct peaks with the central one at the systemic redshift, indicating a ionised tunnel through which both Lya and LyC radiation escape; (b) the large O32 oxygen index ([OIII]5007 / [OII]3727) of 9.18_{-1.32}^{+1.82} is compatible to those measured in local (z~0.4) LyC leakers; (c) there are narrow nebular high-ionisation metal lines with \sigma_v < 20 km/s, which confirms the presence of young hot, massive stars. The HeII 1640 appears qualitatively broad, consistent with a young stellar component including Wolf-Rayet stars. Similarly, the Sunburst LyC-emitter shows a triple-peaked Lya profile and spectral features arising from young hot and massive stars. The strong lensing magnification, (\mu > 30-50), makes this exceptional object the very likely first case of gravitationally-bound star cluster ever observed at a cosmological distance, with a stellar mass M ~< 10^7 Msun and an effective radius smaller than 20 pc. Intriguingly, sources like Sunburst but without lensing magnification might appear as Ion2-like galaxies, in which spatially unresolved massive star clusters dominate the ultraviolet emission. This work supports the idea that dense young star clusters can contribute to the ionisation of the IGM through holes created by stellar feedback.


1904.07947
Fast radio bursts
Petroff, et al

The discovery of radio pulsars over a half century ago was a seminal moment in astronomy. It demonstrated the existence of neutron stars, gave a powerful observational tool to study them, and has allowed us to probe strong gravity, dense matter, and the interstellar medium. More recently, pulsar surveys have led to the serendipitous discovery of fast radio bursts (FRBs). While FRBs appear similar to the individual pulses from pulsars, their large dispersive delays suggest that they originate from far outside the Milky Way and hence are many orders-of-magnitude more luminous. While most FRBs appear to be one-off, perhaps cataclysmic events, two sources are now known to repeat and thus clearly have a longer-lived central engine. Beyond understanding how they are created, there is also the prospect of using FRBs -- as with pulsars -- to probe the extremes of the Universe as well as the otherwise invisible intervening medium. Such studies will be aided by the high implied all-sky event rate: there is a detectable FRB roughly once every minute occurring somewhere on the sky. The fact that less than a hundred FRB sources have been discovered in the last decade is largely due to the small fields-of-view of current radio telescopes. A new generation of wide-field instruments is now coming online, however, and these will be capable of detecting multiple FRBs per day. We are thus on the brink of further breakthroughs in the short-duration radio transient phase space, which will be critical for differentiating between the many proposed theories for the origin of FRBs. In this review, we give an observational and theoretical introduction at a level that is accessible to astronomers entering the field.


1904.07999
Improved photomeson model for interactions of cosmic ray nuclei
Morejon, et al

Photon-hadronic interactions are important for the sources and the transport of Ultra-High Energy Cosmic Rays (UHECRs). Current state-of-the-art cosmic ray transport simulations handle nuclear disintegration at energies of the Giant Dipole Resonance at a more sophisticated level, as well as the photohadronic interactions of nucleons in the high-energy regime above the pion production threshold. However, the interactions of nuclei above the pion production threshold are commonly modeled by treating the nucleus as a superposition of free nucleons -- ignoring the effect of the nuclear medium. We construct an improved, inclusive model for the photomeson regime for nuclei with $A \leq 56$ by employing more accurate, data-driven parametrizations of the interaction cross section, the fragmentation of the primary nucleus and the inclusive pion production cross section that directly affects the production of astrophysical neutrinos. We apply our results to two multi-messenger scenarios (Tidal Disruption Events and Gamma-Ray Bursts) in which photonuclear interactions in the photomeson regime are the dominant cooling process for the highest energy cosmic rays. While we find moderate changes to the mass composition of UHECRs, the astrophysical neutrino fluxes exhibit a significant (factor of a few) reduction compared to the na\"ive superposition of free nucleons for sources of UHECR nuclei with a populated cascade. The numerical code implementing the model has been made publicly available, which facilitates the integration of our results in similar frameworks.


1904.08160
Cosmic ray models
Kachelriess, Semikoz

We review progress in high-energy cosmic ray physics focusing on recent experimental results and models developed for their interpretation. Emphasis is put on the propagation of charged cosmic rays, covering the whole range from $\sim (20-50)$\,GV, i.e.\ the rigidity when solar modulations can be neglected, up to the highest energies observed. We discuss models aiming to explain the anomalies in Galactic cosmic rays, the knee, and the transition from Galactic to extragalactic cosmic rays.


1904.08210
Fitting B/C cosmic-ray dta in the AMS-02 era: a cookbook
Derome, et al

AMS-02 on the ISS has been releasing data of unprecedented accuracy. This poses new challenges for their interpretation. We refine the methodology to get a statistically sound determination of the cosmic-ray propagation parameters. We inspect the numerical precision of the model calculation, nuclear cross-section uncertainties and energy correlations in data systematic errors. We use the 1D diffusion model in USINE. Our $\chi^2$ analysis includes a covariance matrix of errors for AMS-02 systematics and nuisance parameters to account for cross-section uncertainties. Mock data are used to validate some of our choices. We show that any mis-modelling of nuclear cross-section values or the energy correlation length of the covariance matrix of errors biases the analysis, also making good models ($\chi^2_{min}/{dof}\approx1$) appear as excluded ($\chi^2_{min}/{dof}\gg1$). We provide a framework to mitigate these effects (AMS-02 data are interpreted in a companion paper). New production cross-section data as well as the publication, by the AMS-02 collaboration, of a covariance matrix of errors for each data set would be an important step towards an unbiased view of cosmic-ray propagation in the Galaxy.


1904.08344
Detection of Cosmic Rays from ground: an introduction
Di Sciascio

Cosmic rays are the most outstanding example of accelerated particles. They are about 1\% of the total mass of the Universe, so that cosmic rays would represent by far the most important energy transformation process of the Universe. Despite large progresses in building new detectors and in the analysis techniques, the key questions concerning origin, acceleration and propagation of the radiation are still open. One of the reasons is that there are significant discrepancies among the different results obtained by experiments located at ground probably due to unknown systematic errors affecting the measurements. In this note we will focus on detection of Galactic CRs from ground with EAS arrays. This is not a place for a complete review of CR physics (for which we recommend, for instance \cite{spurio,gaisser,grieder,longair,kampert,blasi,kachelriess}) but only to provide elements useful to understand the basic techniques used in reconstructing primary particle characteristics (energy, mass and arrival direction) from ground, and to show why indirect measurements are difficult and results still conflicting.

Wednesday, April 17, 2019

Day 1549

Wednesday.



1904.07237
Next generation cosmography with strong lensing and stellar dynamics
Yildirim, Suyu, Halkola

We present a joint strong lensing and stellar dynamical framework for future time-delay cosmography purposes. Based on a pixelated source reconstruction and the axisymmetric Jeans equations, we are capable of constraining cosmological distances and hence the current expansion rate of the Universe ($H_0$) to the few percent level per lens, when high signal-to-noise integral field unit (IFU) observations from the next generation of telescopes become available. For illustrating the power of this method, we mock up IFU stellar kinematic data of the prominent lens system RXJ1131-1231, given the specifications of the James Webb Space Telescope. Our analysis shows that the time-delay distance ($D_{\Delta t}$) can be constrained with 2.3% uncertainty at best, if future IFU stellar kinematics are included in the fit and if the set of candidate model parameterisations contains the true lens potential. These constraints would translate to a 2.4% precision measurement on $H_0$ from the single lens RXJ1131-1231, assuming flat $\Lambda$CDM cosmology, and can be expected to yield an $H_0$ measure with 1.4% uncertainty, when only three such distance measurements are combined. In addition, the angular diameter distance ($D_\mathrm{d}$) to RXJ1131-1231 can be constrained with 1.8% precision, providing two distance measurements from a single lens system, which is extremely powerful to further constrain the matter density ($\Omega_\mathrm{m}$). The measurement accuracy of $D_\mathrm{d}$, however, is highly sensitive to any systematics in the measurement of the stellar kinematics. For both distance measurements, we strongly advise to probe a large set of physically motivated lens potentials in the future, to minimise the systematic errors associated with the lens mass parameterisation.


1904.07239
Generalised model-independent characterization of strong gravitational lenses VI: the origin of the formalism intrinsic degeneracies and their influence on $H_0$
Wagner

We give a physical interpretation of the formalism intrinsic degeneracies of the gravitational lensing formalism that we derived on a mathematical basis in part IV of this series. We find that all degeneracies occur due to the partition of the mass density along the line of sight. Usually, it is partitioned into a background (cosmic) density and a foreground deflecting object. The latter can be further partitioned into a main deflecting object and perturbers. Weak deflecting objects along the line of sight are also added, either to the deflecting object or as a correction of the angular diameter distances, perturbing the cosmological background density. A priori, this is an arbitrary choice of reference frame and partition. They can be redefined without changing the lensing observables which are sensitive to the integrated deflecting mass density along the entire line of sight. Reformulating the time delay equation such that this interpretation of the degeneracies becomes easily visible, we note that the source can be eliminated from this formulation, which simplifies reconstructions of the deflecting mass distribution or the inference of the Hubble constant, $H_0$. Subsequently, we list necessary conditions to break the formalism intrinsic degeneracies and discuss ways to break them by model choices or including non-lensing observables like velocity dispersions along the line of sight with their advantages and disadvantages. We conclude with a systematic summary of all formalism intrinsic degeneracies and possibilities to break them.


1904.07248
Machine learning in Astronomy: a practical overview
Baron

Astronomy is experiencing a rapid growth in data size and complexity. This change fosters the development of data-driven science as a useful companion to the common model-driven data analysis paradigm, where astronomers develop automatic tools to mine datasets and extract novel information from them. In recent years, machine learning algorithms have become increasingly popular among astronomers, and are now used for a wide variety of tasks. In light of these developments, and the promise and challenges associated with them, the IAC Winter School 2018 focused on big data in Astronomy, with a particular emphasis on machine learning and deep learning techniques. This document summarizes the topics of supervised and unsupervised learning algorithms presented during the school, and provides practical information on the application of such tools to astronomical datasets. In this document I cover basic topics in supervised machine learning, including selection and preprocessing of the input dataset, evaluation methods, and three popular supervised learning algorithms, Support Vector Machines, Random Forests, and shallow Artificial Neural Networks. My main focus is on unsupervised machine learning algorithms, that are used to perform cluster analysis, dimensionality reduction, visualization, and outlier detection. Unsupervised learning algorithms are of particular importance to scientific research, since they can be used to extract new knowledge from existing datasets, and can facilitate new discoveries.


1904.07254
Linking gas and galaxies at high redshift: MUSE surveys the environments of six damped Lyman alpha galaxies at z~3
Mackenzie, et al

We present results from a survey of galaxies in the fields of six z>3 Damped Lyman alpha systems (DLAs) using the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope (VLT). We report a high detection rate of up to ~80% of galaxies within 1000 km/s from DLAs and with impact parameters between 25 and 280 kpc. In particular, we discovered 5 high-confidence Lyman alpha emitters associated with three DLAs, plus up to 9 additional detections across five of the six fields. The majority of the detections are at relatively large impact parameters (>50 kpc) with two detections being plausible host galaxies. Among our detections, we report four galaxies associated with the most metal-poor DLA in our sample (Z/Z_sun = -2.33), which trace an overdense structure resembling a filament. By comparing our detections with predictions from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) cosmological simulations and a semi-analytic model designed to reproduce the observed bias of DLAs at z>2, we conclude that our observations are consistent with a scenario in which a significant fraction of DLAs trace the neutral regions within halos with a characteristic mass of 10^11-10^12 M_sun, in agreement with the inference made from the large-scale clustering of DLAs. We finally show how larger surveys targeting ~25 absorbers have the potential of constraining the characteristic masses of halos hosting high-redshift DLAs with sufficient accuracy to discriminate between different models.


1904.07323
Causation of late quaternary rapid-increase radiocarbon anomalies
Brakenridge

Brief (less than 100 years) rapid-increase anomalies in the Earth's atmospheric radiocarbon production have previously been attributed to either gamma photon radiation from supernovae or to cosmic ray particle radiation from exceptionally large solar flares. Analysis of distances and ages of nearby supernovae remnants, the probable gamma emissions, the predicted Earth incident radiation, and the terrestrial radiocarbon record indicates that supernova causation may be the case. Supernovae include Type Ia white dwarf explosions, Type Ib, c, and II core collapse events, and some types of gamma burst objects. All generate significant pulses of atmospheric radiocarbon depending on distances. Surveys of supernova remnants offer a nearly complete accounting for the past 50,000 years. There are 18 events less than or at 1.4 kilo-parsec distance, and brief radiocarbon anomalies with appropriate sizes occurred for each of the closest events. In calendar years before 1950, these are: Vela, 22 per mil del 14C at 12,760; S165, 20 per mil at 7431; Vela Junior, 13 per mil at 2765; HB9, 9 per mil at 5372; Boomerang, 11 per mil at 10,255; and Cygnus Loop (per mil change not calculated) at 14,722. Although uncertainties remain large, the agreements of prediction to observation support a possible causal connection.


1904.07471
The lowest detected stellar Fe abundance: the halo star SMSS J160540.18-14323.1
Nordlander, et al

We report the discovery of SMSS J160540.18-144323.1, a new ultra-metal poor halo star discovered with the SkyMapper telescope. We measure [Fe/H] = -6.2 +- 0.2 (1D LTE), the lowest ever detected abundance of iron in a star. The star is strongly carbon-enhanced, [C/Fe] = 3.9 +- 0.2, while other abundances are compatible with an alpha-enhanced solar-like pattern with [Ca/Fe] = 0.4 +- 0.2, [Mg/Fe] = 0.6 +- 0.2, [Ti/Fe] = 0.8 +- 0.2, and no significant s- or r-process enrichment, [Sr/Fe] < 0.2 and [Ba/Fe] < 1.0 (3{\sigma} limits). Population III stars exploding as fallback supernovae may explain both the strong carbon enhancement and the apparent lack of enhancement of odd-Z and neutron-capture element abundances. Grids of supernova models computed for metal-free progenitor stars yield good matches for stars of about 10 solar mass imparting a low kinetic energy on the supernova ejecta, while models for stars more massive than roughly 20 solar mass are incompatible with the observed abundance pattern.


1904.07524
The mass-richness relation of optically-selected clusters from weak gravitational lensingand abundance with Subaru HSC first-year data
Murata, et al

Constraining the relation between the richness $N$ and the halo mass $M$ over a wide redshift range for optically-selected clusters is a key ingredient for cluster-related science in optical surveys, including the Subaru Hyper Suprime-Cam (HSC) survey. We measure stacked weak lensing profiles around 1747 HSC CAMIRA clusters over a redshift range of $0.1\leq z_{\rm cl}\leq 1.0$ with $N\geq 15$ using the HSC first-year shear catalog covering $\sim$$140$ ${\rm deg^2}$. The exquisite depth and image quality of the HSC survey allow us to measure lensing signals around the high-redshift clusters at $0.7\leq z_{\rm cl}\leq 1.0$ with a signal-to-noise ratio of 19 in the comoving radius range $0.5\lesssim R\lesssim 15 h^{-1}{\rm Mpc}$. We constrain richness-mass relations $P(\ln N|M,z)$ of the HSC CAMIRA clusters assuming a log-normal distribution without informative priors on model parameters, by jointly fitting to the lensing profiles and abundance measurements under both Planck and WMAP cosmological models. We show that our model gives acceptable $p$-values when we add redshift dependent terms which are proportional to $\ln (1+z)$ and $[\ln (1+z)]^{2}$ into the mean and scatter relations of $P(\ln N|M,z)$. Such terms presumably originate from the variation of photometric redshift errors as a function of the redshift. We show that the constraints on the mean relation $\langle M|N \rangle$ are consistent between the Planck and WMAP models, whereas the scatter values $\sigma_{\ln M|N}$ for the Planck model are systematically larger than those for the WMAP model. We also show that the scatter values for the Planck model increase toward lower richness values, whereas those for the WMAP model are consistent with constant values as a function of richness. This result highlights the importance of the scatter in the mass-richness relation for cluster cosmology.


1904.07718
OGLE-2017-BLG-1186: first application of astroseismology and Gaussian processes to microlensing
Li, et al

We present the analysis of the event OGLE-2017-BLG-1186 from the 2017 $Spitzer$ microlensing campaign. This event is remarkable because its source is photometrically variable. We perform an asteroseismic analysis of the source star and find that it is an oscillating red giant with average timescale of $\sim 9$ days. The asteroseismic analysis also provides us source properties including the source angular size ($\sim 27~\mu{\rm as}$) and distance ($\sim 11.5$ kpc), which are essential for inferring the properties of the lens. When fitting the light curve, we test the feasibility of Gaussian Processes (GPs) in handling the correlated noise caused by the variable source. We find, in this event, that the GP model cannot provide better constraints on the parameters of interest because of the worse estimation of blending fluxes compared to the traditional $\chi^2$ minimization method. We note that this event is the first microlensing system for which asteroseismology and GPs have been used to account for the variable source. With both finite-source effect and microlens parallax measured, we find that the lens is likely a $\sim 0.045~M_{\odot}$ brown dwarf at distance $\sim 9.0$ kpc, or a $\sim 0.073 M_{\odot}$ ultracool dwarf at distance $\sim 9.8$ kpc. Combining the estimated lens properties with a Bayesian analysis using a Galactic model, we find a $\sim 35\%$ probability for the lens to be a bulge object and $\sim 65\%$ to be a background disk object.


1904.07727
Using cosmogenic Lithium, Beryllium and Boron to determine the surface ages of icy objects n the outer solar system
Hedman

Given current uncertainties in the cratering rates and geological histories of icy objects in the outer solar system, it is worth considering how the ages of icy surfaces could be constrained with measurements from future landed missions. A promising approach would be to determine cosmic-ray exposure ages of surface deposits by measuring the amounts of cosmogenic Lithium, Beryllium and Boron at various depths within a few meters of the surface. Preliminary calculations show that ice that has been exposed to cosmic radiation for one billion years should contain these cosmogenic nuclei at concentrations of a few parts per trillion, so any future experiment that might attempt to perform this sort of measurement will need to meet stringent sensitivity requirements.


1904.07747
Ly$\alpha$ view around a z=2.84 hyper luminous QSO at a node of the cosmic web
Kikuta, et al

We report on the results of deep and wide-field (1.1 deg$^2$) narrow-band observations with Subaru/Hyper Suprime-Cam (HSC) of a field around a hyperluminous QSO (HLQSO), HS1549+1919, residing in a protocluster at $z=2.84$, to map the large-scale structure of Ly$\alpha$ emitters (LAEs). One HSC pointing enables us to detect 3490 LAEs and 76 extended Ly$\alpha$ blobs (LABs), probing diverse environments from voids to protoclusters. The HLQSO is found to be near the center of the protocluster, which corresponds to the intersection of $\sim$100 cMpc-scale structures of LAEs. LABs are basically distributed along the large-scale structure, with larger ones particularly clustered around the HLQSO, confirming a previously noted tendency of LABs to prefer denser environments. Moreover, the shapes of LABs near the HLQSO appear to be aligned with the large-scale structure. Finally, a deep Ly$\alpha$ image reveals a diffuse Ly$\alpha$ nebula along a filamentary structure with no luminous UV/sub-mm counterpart. We suggest that the diffuse nebula is due to a cold filament with high clumping factor illuminated by the QSO, with a required high clumpiness provided by unresolved residing halos of mass $\leq 10^{9-10}M_\odot$.


1904.07749
Resolving the interstellar medium in the nuclear region of two z=5.78 quasar host galaxies with ALMA
Wang, et al

We present ALMA observations of the [CII] 158 micron fine structure line and dust continuum emission from two quasars, SDSSJ104433.04-012502.2 and SDSSJ012958.51-003539.7, at z=5.78. The ALMA observations at 0.2" resolution map the dust and gas on kpc scales. The spatially resolved emission show a similar trend of increasing [CII]-FIR ratios with decreasing FIR surface brightnesses toward larger radii as was found in the infrared luminous galaxies with intense star formation. We confirm the velocity gradients of [CII] emission found previously in SDSS J0129-0035. No evidence of order motion is detected in SDSS J1044-0125. The velocity maps and PV-diagrams also suggest turbulent gas clumps in both objects. We tentatively detect a [CII] peak offset 4.9 kpc to the East of SDSS J1044-0125. This may be associated with an infalling companion, or node of gas outflows. All these results suggest significant dynamical evolution of the ISM in the nuclear region of these young quasar-starburst systems. We fit the velocity map of the [CII] emission from SDSS J0129-0035 with a rotating disk model. The result suggests a face-on system with an inclination angle of 16+/-20 degree and constrains the lower limit of the host galaxy dynamical mass to be 2.6x10^10 Msun within the [CII] emitting region. It is likely that SDSS J0129-0035, as well as other young quasars with super massive black hole masses on order of 10^7 Msun to 10^8 Msun, falls close to the black hole and host galaxy mass relation defined by local galaxies.