Wednesday, November 25, 2020

Day 1781

Monday, Tuesday.  Wednesday.


2011.10553
Mapping Gaia parallax systematic errors over the sky with faint Milky Way stars
Fardal, et al

Parallaxes measured by the Gaia mission have huge significance for astronomy, but parallaxes in Gaia DR2 are known to have systematic errors that depend on the source position and other quantities. We use the abundant information in faint Milky Way stars, along with the GOG simulation of the Gaia catalog, to probe the spatial dependence of Gaia DR2 parallax systematic errors in an empirical way. The parallax signal, concentrated in thick disk turnoff stars with magnitude G ~ 17, is sufficient to construct maps of the parallax systematic error over the majority of the sky. These maps show a locally regular "waffle pattern" on ~1 degree scales following Gaia scan directions, stronger linear "scar" features, and coherent variations on larger scales. The parallax bias maps also retain traces of astrophysical effects such as dust clouds. The waffle pattern, known from earlier maps of the Magellanic Clouds, extends over the entire sky; its local rms amplitude averages 15 microarcsec and varies by about a factor of two. The strength of this pattern increases by a factor ~6 from magnitude G = 13 to G = 20. Correlations with parallaxes of quasars and of stars with independent distance estimates support our bias estimates. Using similar methods, we map systematic errors in the proper motion and examine the relationship with the parallax systematics. We provide a code package to access and query our bias maps. Similar tests on the general stellar population should be useful in quantifying systematic errors in future Gaia releases.


2011.10578
Completeness of the Gaia-verse III: using hidden states to infer gaps, detection efficiencies and the scanning law from the DR2 light curves
Boubert, et al

The completeness of the Gaia catalogues heavily depends on the status of that space telescope through time. Stars are only published with each of the astrometric, photometric and spectroscopic data products if they are detected a minimum number of times. If there is a gap in scientific operations, a drop in the detection efficiency or Gaia deviates from the commanded scanning law, then stars will miss out on potential detections and thus be less likely to make it into the Gaia catalogues. We lay the groundwork to retrospectively ascertain the status of Gaia throughout the mission from the tens of individual measurements of the billions of stars, by developing novel methodologies to infer both the orientation and angular velocity of Gaia through time and gaps and efficiency drops in the detections. We have applied these methodologies to the Gaia DR2 variable star epoch photometry -- which are the only publicly available Gaia time-series at the present time -- and make the results publicly available. We accompany these results with a new Python package scanninglaw (https://github.com/gaiaverse/scanninglaw) that you can use to easily predict Gaia observation times and detection probabilities for arbitrary locations on the sky.


2011.11604
HSTCosmicrays: A Python package for analyzing cosmic rays in HST calibration data
Miles, et al

HSTCosmicrays is a python-based pipeline designed to find and characterize cosmic rays found in dark frames (exposures taken with the shutter closed). Dark exposures are obtained routinely by all the Hubble Space Telescope (HST) instruments for calibration. The main processing pipeline runs locally or in the cloud on AWS. To date, we have characterized more than 1.2 billion cosmic rays in ~76,000 dark frames obtained with CCDs from the four active instruments ACS/HRC, ACS/WFC, STIS, WFC3/UVIS, and the legacy instrument WFPC2.


2011.11613
Strong detection of the CMB lensingxgalaxy weak lensing cross-correlation from ACT-DR4, Planck Legacy and KiDS-1000
Robertson, et al

We measure the cross-correlation between galaxy weak lensing data from the Kilo Degree Survey (KiDS-1000, DR4) and cosmic microwave background (CMB) lensing data from the Atacama Cosmology Telescope (ACT, DR4) and the Planck Legacy survey. We use two samples of source galaxies, selected with photometric redshifts, $(0.1<z_{\rm B}<1.2)$ and $(1.2<z_{\rm B}<2)$, which produce a combined detection significance of the CMB lensing/weak galaxy lensing cross-spectrum of $7.7\sigma$. With the lower redshift galaxy sample, for which the cross-correlation is detected at a significance of $5.3\sigma$, we present joint cosmological constraints on the matter density parameter, $\Omega_{\rm m}$, and the matter fluctuation amplitude parameter, $\sigma_8$, marginalising over three nuisance parameters that model our uncertainty in the redshift and shear calibration, and the intrinsic alignment of galaxies. We find our measurement to be consistent with the best-fitting flat $\Lambda$CDM cosmological models from both Planck and KiDS-1000. We demonstrate the capacity of CMB-weak lensing cross-correlations to set constraints on either the redshift or shear calibration, by analysing a previously unused high-redshift KiDS galaxy sample $(1.2<z_{\rm B}<2)$, with the cross-correlation detected at a significance of $7\sigma$. This analysis provides an independent assessment for the accuracy of redshift measurements in a regime that is challenging to calibrate directly owing to known incompleteness in spectroscopic surveys.


2011.11648
Unraveling the origin of magnetic fields in galaxies
Martin-Alvarez, et al

Despite their ubiquity, there are many open questions regarding galactic and cosmic magnetic fields. Specifically, current observational constraints cannot rule out if magnetic fields observed in galaxies were generated in the Early Universe or are of astrophysical nature. Motivated by this we use our magnetic tracers algorithm to investigate whether the signatures of primordial magnetic fields persist in galaxies throughout cosmic time. We simulate a Milky Way-like galaxy in four scenarios: magnetised solely by primordial magnetic fields, magnetised exclusively by SN-injected magnetic fields, and two combined primordial + SN magnetisation cases. We find that once primordial magnetic fields with a comoving strength $B_0 >10^{-12}$ G are considered, they remain the primary source of galaxy magnetisation. Our magnetic tracers show that, even combined with galactic sources of magnetisation, when primordial magnetic fields are strong, they source the large-scale fields in the warm metal-poor phase of the simulated galaxy. In this case, the circumgalactic and intergalactic medium can be used to probe $B_0$ without risk of pollution by magnetic fields originated in the galaxy. Furthermore, whether magnetic fields are primordial or astrophysically-sourced can be inferred by studying local gas metallicity. As a result, we predict that future state-of-the-art observational facilities of magnetic fields in galaxies will have the potential to unravel astrophysical and primordial magnetic components of our Universe.

Monday, November 23, 2020

Day 1780

Wednesday.  Thursday.  Friday.


2011.08385
First constraints on the intrinsic CMB dipole and our velocity with Doppler and aberration
Ferreira, Quartin

The CMB dipole is usually assumed to be completely due to the relative velocity between the solar system and the CMB restframe. We test this hypothesis by measuring independently the Doppler and aberration effects on the CMB using Planck 2018 data. We remove the contributions arising solely from the conversion of intensity into temperature and arrive at measurements which are independent from the CMB dipole itself. Doppler and aberration are measured independently in the TT and EE maps for both SMICA and NILC. Combining these new measurements with the dipole one we get the first constraints on the intrinsic CMB dipole. Neglecting a dipolar lensing contribution we can put an upper limit on its amplitude: 3.5 mK (95% CI). We also obtain the estimate of the peculiar velocity of the solar system which does not assume a negligible intrinsic dipole contribution: $v = (300 \pm 99)$ km/s with $(l,b) = (276 \pm 32, 51 \pm 19)^\circ$ [SMICA] and $v = (300 \pm 100)$ km/s with $(l,b) = (280 \pm 32, 50 \pm 20)^\circ$ [NILC] with negligible systematic contributions. These values are consistent with the peculiar velocity hypothesis of the dipole.


2011.09252
Effects of the COVID-19 lockdown on urban light emissions: ground and satellite comparison
Bustamante-Calabria, et al

'Lockdown' periods in response to COVID-19 have provided a unique opportunity to study the impacts of economic activity on environmental pollution (e.g. NO$_2$, aerosols, noise, light). The effects on NO$_2$ and aerosols have been very noticeable and readily demonstrated, but that on light pollution has proven challenging to determine. The main reason for this difficulty is that the primary source of nighttime satellite imagery of the earth is the SNPP-VIIRS/DNB instrument, which acquires data late at night after most human nocturnal activity has already occurred and much associated lighting has been turned off. Here, to analyze the effect of lockdown on urban light emissions, we use ground and satellite data for Granada, Spain, during the COVID-19 induced confinement of the city's population from March 14 until May 31, 2020. We find a clear decrease in light pollution due both to a decrease in light emissions from the city and to a decrease in anthropogenic aerosol content in the atmosphere which resulted in less light being scattered. A clear correlation between the abundance of PM10 particles and sky brightness is observed, such that the more polluted the atmosphere the brighter the urban night sky. An empirical expression is determined that relates PM10 particle abundance and sky brightness at three different wavelength bands.


2011.09835
Multi-CCD point spread function modelling
Liaudat, et al

Galaxy imaging surveys observe a vast number of objects that are affected by the instrument's Point Spread Function (PSF). Weak lensing missions, in particular, aim at measuring the shape of galaxies, and PSF effects represent an important source of systematic errors which must be handled appropriately. This demands a high accuracy in the modelling as well as the estimation of the PSF at galaxy positions. Sometimes referred to as non-parametric PSF estimation, the goal of this paper is to estimate a PSF at galaxy positions, starting from a set of noisy star image observations distributed over the focal plane. To accomplish this, we need our model to first of all, precisely capture the PSF field variations over the Field of View (FoV), and then to recover the PSF at the selected positions. This paper proposes a new method, coined MCCD (Multi-CCD PSF modelling), that creates, simultaneously, a PSF field model over all of the instrument's focal plane. This allows to capture global as well as local PSF features through the use of two complementary models which enforce different spatial constraints. Most existing non-parametric models build one model per Charge Coupled Device (CCD), which can lead to difficulties in capturing global ellipticity patterns. We first test our method on a realistic simulated dataset comparing it with two state-of-the-art PSF modelling methods (PSFEx and RCA). We outperform both of them with our proposed method. Then we contrast our approach with PSFEx on real data from CFIS (Canada France Imaging Survey) that uses the CFHT (Canada-France-Hawaii Telescope). We show that our PSF model is less noisy and achieves a 22% gain on pixel Root Mean Squared Error (RMSE) with respect to PSFEx. We present, and share the code of, a new PSF modelling algorithm that models the PSF field on all the focal plane that is mature enough to handle real data.

Tuesday, November 17, 2020

Day 1779

Monday, Tuesday.


2011.07722
Bayesian forward modeling of cosmic shear data
Porqueres, et al

We present a Bayesian hierarchical modelling approach to infer the cosmic matter density field, and the lensing and the matter power spectra, from cosmic shear data. This method uses a physical model of cosmic structure formation to infer physically plausible cosmic structures, which accounts for the non-Gaussian features of the gravitationally evolved matter distribution and light-cone effects. We test and validate our framework with realistic simulated shear data, demonstrating that the method recovers the unbiased matter distribution and the correct lensing and matter power spectrum. While the cosmology is fixed in this test, and the method employs a prior power spectrum, we demonstrate that the lensing results are sensitive to the true power spectrum when this differs from the prior. In this case, the density field samples are generated with a power spectrum that deviates from the prior, and the method recovers the true lensing power spectrum. The method also recovers the matter power spectrum across the sky, but as currently implemented, it cannot determine the radial power since isotropy is not imposed. In summary, our method provides physically plausible inference of the dark matter distribution from cosmic shear data, allowing us to extract information beyond the two-point statistics and exploiting the full information content of the cosmological fields.

Sunday, November 15, 2020

Day 1778

Wednesday, Thursday, Friday.


2011.04671
The limits of cosmology: role of the Moon
Silk

The lunar surface allows a unique way forward in cosmology, to go beyond current limits. The far side provides an unexcelled radio-quiet environment for probing the dark ages via 21 cm interferometry to seek elusive clues on the nature of the infinitesimal fluctuations that seeded galaxy formation. Far-infrared telescopes in cold and dark lunar polar craters will probe back to the first months of the Big Bang and study associated spectral distortions in the CMB. Optical and IR megatelescopes will image the first star clusters in the universe and seek biosignatures in the atmospheres of unprecedented numbers of nearby habitable zone exoplanets. The goals are compelling and a stable lunar platform will enable construction of telescopes that can access trillions of modes in the sky, providing the key to exploration of our cosmic origins.


2011.04759
International coordination and support for SmallSat-enabled space weather activities
Nieves-Chinchilla, et al

Advances in space weather science and small satellite (SmallSat) technology have proceeded in parallel over the past two decades, but better communication and coordination is needed among the respective worldwide communities contributing to this rapid progress. We identify six areas where improved international coordination is especially desirable, including: (1) orbital debris mitigation; (2) spectrum management; (3) export control regulations; (4) access to timely and low-cost launch opportunities; (5) inclusive data policies; and (6) education. We argue the need for internationally coordinated policies and programs to promote the use of SmallSats for space weather research and forecasting while realizing maximum scientific and technical advances through the integration of these two increasingly important endeavors.


2011.05270
Identification and mitigation of a vibrational telescope systematic with application to Spitzer
Challener, et al

We observed Proxima Centauri with the Spitzer Space Telescope InfraRed Array Camera (IRAC) five times in 2016 and 2017 to search for transits of Proxima Centauri b. Following standard analysis procedures, we found three asymmetric, transit-like events that are now understood to be vibrational systematics. This systematic is correlated with the width of the point-response function (PRF), which we measure with rotated and non-rotated Gaussian fits with respect to the detecor array. We show that the systematic can be removed with a novel application of an adaptive elliptical-aperture photometry technique, and compare the performance of this technique with fixed and variable circular-aperture photometry, using both BiLinearly Interpolated Subpixel Sensitivity (BLISS) maps and non-binned Pixel-Level Decorrelation (PLD). With BLISS maps, elliptical photometry results in a lower standard deviation of normalized residuals, and reduced or similar correlated noise when compared to circular apertures. PLD prefers variable, circular apertures, but generally results in more correlated noise than BLISS. This vibrational effect is likely present in other telescopes and Spitzer observations, where correction could improve results. Our elliptical apertures can be applied to any photometry observations, and may be even more effective when applied to more circular PRFs than Spitzer's.


2011.05297
Euclid: Identification of asteroid streaks in simulated images using StreakDet software
Pöntinen, et al

The ESA Euclid space telescope could observe up to 150 000 asteroids as a side product of its primary cosmological mission. Asteroids appear as trailed sources, that is streaks, in the images. Owing to the survey area of 15 000 square degrees and the number of sources, automated methods have to be used to find them. Euclid is equipped with a visible camera, VIS (VISual imager), and a near-infrared camera, NISP (Near-Infrared Spectrometer and Photometer), with three filters. We aim to develop a pipeline to detect fast-moving objects in Euclid images, with both high completeness and high purity. We tested the StreakDet software to find asteroids from simulated Euclid images. We optimized the parameters of StreakDet to maximize completeness, and developed a post-processing algorithm to improve the purity of the sample of detected sources by removing false-positive detections. StreakDet finds 96.9% of the synthetic asteroid streaks with apparent magnitudes brighter than 23rd magnitude and streak lengths longer than 15 pixels ($10\,{\rm arcsec\,h^{-1}}$), but this comes at the cost of finding a high number of false positives. The number of false positives can be radically reduced with multi-streak analysis, which utilizes all four dithers obtained by Euclid. StreakDet is a good tool for identifying asteroids in Euclid images, but there is still room for improvement, in particular, for finding short (less than 13 pixels, corresponding to 8$\,{\rm arcsec\,h^{-1}}$) and/or faint streaks (fainter than the apparent magnitude of 23).


2011.05345
Cosmology dependency of halo masses and concentrations in hydrodynamic simulations
Ragagnin, et al

We employ a set of Magneticum cosmological hydrodynamic simulations that span over $15$ different cosmologies, and extract masses and concentrations of all well-resolved haloes between $z=0-1$ for critical over-densities $\Delta_\texttt{vir}, \Delta_{200c}, \Delta_{500c}, \Delta_{2500c}$ and mean overdensity $\Delta_{200m}.$ We provide the first mass-concentration (Mc) relation and sparsity relation (i.e. $M_{\Delta1} - M_{\Delta2}$ mass conversion) of hydrodynamic simulations that is modelled by mass, redshift and cosmological parameters $\Omega_m, \Omega_b, \sigma_8, h_0$ as a tool for observational studies. We also quantify the impact that the Mc relation scatter and the assumption of NFW density profiles have on the uncertainty of the sparsity relation. We find that converting masses with the aid of a Mc relation carries an additional fractional scatter ($\approx 4\%$) originated from deviations from the assumed NFW density profile. For this reason we provide a direct mass-mass conversion relation fit that depends on redshift and cosmological parameters. We release the package hydro\_mc, a python tool that perform all kind of conversions presented in this paper.


2011.05347
INSPIRE: INvestigating Stellar Populations In RElics -- I. Survey presentation and pilot program
Spiniello, et al

Massive ETGs are thought to form through a two-phase process. At early times, an intense and fast starburst forms blue and disk-dominated galaxies. After quenching, the remaining structures become red, compact and massive, i.e., 'red nuggets'. Then, a time-extended second phase which is dominated by mergers, causes structural evolution and size growth. Given the stochastic nature of mergers, a small fraction of red nuggets survives, without any interaction, massive and compact until today: relic galaxies. Since this fraction depends on the processes dominating the size growth, counting relics at low-z is a valuable way to disentangle between different galaxy evolution models. In this paper, we introduce the INvestigating Stellar Population In RElics (INSPIRE) Project, that aims at spectroscopically confirming and fully characterizing a large number of relics at 0.1<z<0.5. We focus here on the first results based on a pilot program targeting three systems, representative of the whole sample. For these, we extract 1D optical spectra over an aperture comprising ~30 % of the galaxies light, and obtain line-of-sight integrated stellar velocity and velocity dispersion. We also infer the stellar [$\alpha$/Fe] abundance from line-index measurements and mass-weighted age and metallicity from full-spectral fitting with single stellar population models. Two galaxies have large integrated stellar velocity dispersion values, confirming their massive nature. They are populated by stars with super-solar metallicity and [$\alpha$/Fe]. Both objects have formed >80 % of their stellar mass within a short (0.5 - 1.0 Gyrs) initial star formation episode occurred only ~1 Gyr after the Big Bang. The third galaxy has a more extended star formation history and a lower velocity dispersion. Thus we confirm two out of three candidates as relics.


2011.05751
The impact of the mass spectrum of lenses in quasar microlensing studies.  Constraints on a mixed population of primordial black holes and stars
Esteban-Gutiérrez, et al

We show that quasar microlensing magnification statistics induced by a population of point microlenses distributed according to a mass-spectrum can be very well approximated by that of a single-mass, "monochromatic", population. When the spatial resolution (physically defined by the source size) is small as compared with the Einstein radius, the mass of the monochromatic population matches the geometric mean of the mass-spectrum. Otherwise, the best-fit mass can be larger. Taking into account the degeneracy with the geometric mean, the interpretation of quasar microlensing observations under the hypothesis of a mixed population of primordial black holes and stars, makes the existence of a significant population of intermediate mass black holes ($\sim$ 100$M_\odot$) unlikely but allows, within a two-$\sigma$ confidence interval, the presence of a large population ($\gtrsim 40\%$ of the total mass) of substellar black holes ($\sim$ 0.01$M_\odot$).


2011.05876
Testing the equivalence principle on cosmological scales using the odd multipoles of galaxy cross-power spectrum and bispectrum
Umeh, Koyama, Crittenden

One of the cornerstones of general relativity is the equivalence principle. However, the validity of the equivalence principle has only been established on solar system scales for standard matter fields; this result cannot be assumed to hold for the non-standard matter fields that dominate the gravitational dynamics on cosmological scales. Here we show how the equivalence principle may be tested on cosmological scales for non-standard matter fields using the odd multipoles of the galaxy cross-power spectrum and bispectrum. This test makes use of the imprint on the galaxy cross-power spectrum and bispectrum by the parity-violating general relativistic deformations of the past-light cone, and assumes that galaxies can be treated as test particles that are made of baryons and cold dark matter. This assumption leads to a non-zero galaxy-baryon relative velocity if the equivalence principle does not hold between baryons and dark matter. We show that the relative velocity can be constrained to be less than 25% of the galaxy velocity using the cross-power spectrum of the HI intensity mapping/Halpha galaxy survey and the bispectrum of the Halpha galaxy survey.


2011.06469
Consistency of cosmic shear analyses in harmonic and real space
Doux, et al

Recent cosmic shear studies have reported discrepancies of up to $1\sigma$ on the parameter ${S_{8}=\sigma_{8}\sqrt{\Omega_{\rm m}/0.3}}$ between the analysis of shear power spectra and two-point correlation functions, derived from the same shear catalogs. It is not a priori clear whether the measured discrepancies are consistent with statistical fluctuations. In this paper, we investigate this issue in the context of the forthcoming analyses from the third year data of the Dark Energy Survey (DES-Y3). We analyze DES-Y3 mock catalogs from Gaussian simulations with a fast and accurate importance sampling pipeline. We show that the methodology for determining matching scale cuts in harmonic and real space is the key factor that contributes to the scatter between constraints derived from the two statistics. We compare the published scales cuts of the KiDS, Subaru-HSC and DES surveys, and find that the correlation coefficients of posterior means range from over 80% for our proposed cuts, down to 10% for cuts used in the literature. We then study the interaction between scale cuts and systematic uncertainties arising from multiple sources: non-linear power spectrum, baryonic feedback, intrinsic alignments, uncertainties in the point-spread function, and redshift distributions. We find that, given DES-Y3 characteristics and proposed cuts, these uncertainties affect the two statistics similarly; the differential biases are below a third of the statistical uncertainty, with the largest biases arising from intrinsic alignment and baryonic feedback. While this work is aimed at DES-Y3, the tools developed can be applied to Stage-IV surveys where statistical errors will be much smaller.

Friday, November 13, 2020

Day 1777

Monday, Tuesday, Wednesday.


2011.03369
Asteroid resource utilization: ethical concerns and progress
Rivkin, et al

As asteroid mining moves toward reality, the high bar to entering the business may limit participation and increase inequality, reducing or eliminating any benefit gained by marginalized people or developing nations. Consideration of ethical issues is urgently needed, as well as participation in international, not merely multilateral, solutions.


2011.03407
Dark Energy Survey Year 3 Results: Photometric Data Set for Cosmology
Sevilla-Narbe, et al

2011.03408
Dark Energy Survey Year 3 Results: Weak Lensing Shape Catalogue
Gatti, et al

2011.03409
Dark Energy Survey Year 3 Results: Point-spread function modeling
Jarvis, Bernstein, et al

We introduce a new software package for modeling the point-spread function (PSF) of astronomical images, called Piff (PSFs In the Full FOV), which we apply to the first three years (known as Y3) of the Dark Energy Survey (DES) data. We describe the relevant details about the algorithms used by Piff to model the PSF, including how the PSF model varies across the field of view (FOV). Diagnostic results show that the systematic errors from the PSF modeling are very small over the range of scales that are important for the DES Y3 weak lensing analysis. In particular, the systematic errors from the PSF modeling are significantly smaller than the corresponding results from the DES year one (Y1) analysis. We also briefly describe some planned improvements to Piff that we expect to further reduce the modeling errors in future analyses.

2011.03410
Dark Energy Survey internal consistency tests of the joint cosmological probes analysis with posterior predictive distributions
Doux, et al

2011.03411
Dark Energy Survey Year 3 Results: Optimizing the lens sample in combined galaxy clustering and galaxy-galaxy lensing analysis
Porredon, et al


2011.03985
The Earths long-term climate changes and ice ages: a derivation of Milankovitch cycles from first principles
Rainey

Long-term changes in the tilt of the Earths axis, relative to the plane of its orbit, are of great significance to long-term climate change, because they control the size of the arctic and antarctic circles. These Milankovitch cycles have generally been calculated by numerical integration of Newtons equations of motion, and there is controversy over the results because they are sensitive to numerical drift over the very long computer simulations involved. In this paper the cycles are calculated from first principles, without any reliance on computer simulation. The problem is one of planetary precession, and is solvable by the methods used to study the precession of a spinning top. It is shown that the main component of Milankovitch cycles has a period of 41,000 years and is due to one of the modes of precession of the Earth-Venus system. The other mode of this system produces a component of period 29,500 years, and a third component of period 54,000 years results from the influence of the precession of the orbits of Jupiter and Saturn. These results agree closely with several of the numerical simulations in the literature, and strongly suggest that other different results in the literature are wrong.

Wednesday, November 11, 2020

Day 1776

Tuesday.  Wednesday.  Thursday.  Friday.


2011.01820
Optical to NIR magnitude measurements of the Starlink LEO Darksat satellite and effectiveness of the darkening treatment
Tregloan-Reed, et al

Four observations of Starlink's LEO communication satellites, Darksat and STARLINK-1113 were conducted on two nights with two telescopes. The Chakana 0.6\,m, telescope at the Ckoirama observatory (Chile) observed both satellites on 2020-03-05 (UTC) and 2020-03-07 (UTC) using a Sloan {\it r'} and Sloan {\it i'} filter, respectively. The ESO Visible and Infrared Survey Telescope for Astronomy (VISTA) 4.0\,m telescope with the VISTA InfraRed CAMera (VIRCAM) observed both satellites on 2020/03/05 (UTC) and 2020/03/07 (UTC) in the NIR {\it J}-band and NIR {\it Ks}-band, respectively. The calibration, image processing, and analysis of the Darksat images give r ~5.6\,mag, i ~5.0\,mag, J ~4.2\,mag, and Ks ~4.0\,mag when scaled to a range of 550\,km (airmass $=1$) and corrected for the solar incidence and observer phase angles. In comparison the STARLINK-1113 images give r ~4.9\,mag, i ~4.4\,mag, J ~3.8\,mag, and Ks ~3.6\,mag when corrected for range, solar incidence and observer phase angles. The data and results presented in this work, show that the special darkening coating used by Starlink for Darksat has darkened the Sloan {\it r'} magnitude by 50\,\%, Sloan {\it i'} magnitude by 42\,\%, NIR J magnitude by 32\,\%, NIR Ks magnitude by 28\,\%. The results show that both satellites increase in reflective brightness with increasing wavelength and that the effectiveness of the darkening treatment reduces at longer wavelengths. This shows that the mitigation strategies being developed by Starlink and other LEO satellite operators  This work highlights the continued importance of obtaining multi-wavelength observations of many different LEO satellites in order to characterise their reflective properties, to aid the community in developing impact simulations and develop mitigation tools.


2011.01836
Two-point statistics without bins: a continuous-function generalization of the correlation function estimator for large-scale structure
Storey-Fisher, Hogg

The two-point correlation function (2pcf) is the key statistic in structure formation; it measures the clustering of galaxies or other density field tracers. Estimators of the 2pcf, including the standard Landy-Szalay (LS) estimator, evaluate the 2pcf in hard-edged separation bins, which is scientifically inappropriate and results in a poor trade-off between bias and variance. We present a new 2pcf estimator, the Continuous-Function Estimator, which generalizes LS to a continuous representation and obviates binning in separation or any other pair property. Our estimator, inspired by the mathematics of least-squares fitting, replaces binned pair counts with projections onto basis functions; it outputs the best linear combination of basis functions to describe the 2pcf. The choice of basis can take into account the expected form of the 2pcf, as well as its dependence on pair properties other than separation. We show that the Continuous-Function Estimator with a cubic-spline basis better represents the shape of the 2pcf compared to LS. We also estimate directly the baryon acoustic scale, using a small number of physically-motivated basis functions. Critically, this leads to a reduction in the number of mock catalogs required for covariance estimation, which is currently the limiting step in many 2pcf analyses. We discuss further applications of the Continuous-Function Estimator, including determination of the dependence of clustering on galaxy properties and searches for potential inhomogeneities or anisotropies in large-scale structure.


2011.01945
Reconciling galaxy cluster shapes, measured by theorists vs observers
Harvey, et al

If properly calibrated, the shapes of galaxy clusters can be used to investigate many physical processes: from feedback and quenching of star formation, to the nature of dark matter. Theorists frequently measure shapes using moments of inertia of simulated particles'. We instead create mock (optical, X-ray, strong- and weak-lensing) observations of the twenty-two most massive ($\sim10^{14.7}\,M_\odot$) relaxed clusters in the BAHAMAS simulations. We find that observable measures of shape are rounder. Even when moments of inertia are projected into 2D and evaluated at matched radius, they overestimate ellipticity by 56\% (compared to observable strong lensing) and 430\% (compared to observable weak lensing). Therefore, we propose matchable quantities and test them using observations of eight relaxed clusters from the {\emph Hubble Space Telescope} and {\emph Chandra X-Ray Observatory}. We also release our HST data reduction and lensing analysis software to the community. In real clusters, the ellipticity and orientation angle at all radii are strongly correlated. In simulated clusters, the ellipticity of inner ($<r_{\mathrm{vir}}/20$) regions becomes decoupled: for example with greater misalignment of the central cluster galaxy. This may indicate overly efficient implementation of feedback from active galactic nuclei. Future exploitation of cluster shapes as a function of radii will require better understanding of core baryonic processes. Exploitation of shapes on any scale will require calibration on simulations extended all the way to mock observations.


2011.02377
On the halo-mass and radial scale dependence of the lensing is low effect
Lange, et al

The canonical $\Lambda$CDM cosmological model makes precise predictions for the clustering and lensing properties of galaxies. It has been shown that the lensing amplitude of galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) is lower than expected given their clustering properties. We present new measurements and modelling of galaxies in the BOSS LOWZ sample. We focus on the radial and stellar mass dependence of the lensing amplitude mis-match. We find an amplitude mis-match of around $35\%$ when assuming $\Lambda$CDM with Planck Cosmological Microwave Background (CMB) constraints. This offset is independent of halo mass and radial scale in the range $M_{\rm halo}\sim 10^{13.3} - 10^{13.9} h^{-1} M_\odot$ and $r=0.1 - 60 \, h^{-1} \mathrm{Mpc}$ ($0.05\ h/{\rm Mpc} \lesssim k \lesssim 20 \ h/{\rm Mpc}$). The observation that the offset is both mass and scale independent places important constraints on the degree to which astrophysical processes (baryonic effects, assembly bias) can fully explain the effect. This scale independence also suggests that the "lensing is low" effect on small and large radial scales probably have the same physical origin. Resolutions based on new physics require a nearly uniform suppression, relative to $\Lambda$CDM predictions, of the amplitude of matter fluctuations on these scales. The possible causes of this are tightly constrained by measurements of the CMB and of the low-redshift expansion history.

Monday, November 2, 2020

Day 1775

Monday.


2010.16416
KiDS-1000 Cosmology: constraints beyond flat $\Lambda$CDM
Tröster, et al

We present constraints on extensions to the flat $\Lambda$CDM cosmological model by varying the spatial curvature $\Omega_K$, the sum of the neutrino masses $\sum m_\nu$, the dark energy equation of state parameter $w$, and the Hu-Sawicki $f(R)$ gravity $f_{R0}$ parameter. With the combined $3\times2$pt measurements of cosmic shear from the Kilo-Degree Survey (KiDS-1000), galaxy clustering from the Baryon Oscillation Spectroscopic Survey (BOSS), and galaxy-galaxy lensing from the overlap between KiDS-1000, BOSS, and the spectroscopic 2-degree Field Lensing Survey (2dFLenS), we find results that are fully consistent with a flat $\Lambda$CDM model with $\Omega_K=0.011^{+0.054}_{-0.057}$, $\sum m_\nu<1.76$ eV (95% CL), and $w=-0.99^{+0.11}_{-0.13}$. The $f_{R0}$ parameter is unconstrained in our fully non-linear $f(R)$ cosmic shear analysis. Considering three different model selection criteria, we find no clear preference for either the fiducial flat $\Lambda$CDM model or any of the considered extensions. Besides extensions to the flat $\Lambda$CDM parameter space, we also explore restrictions to common subsets of the flat $\Lambda$CDM parameter space by fixing the amplitude of the primordial power spectrum to the Planck best-fit value, as well as adding external data from supernovae and lensing of the CMB. Neither the beyond-$\Lambda$CDM models nor the imposed restrictions explored in this analysis are able to resolve the $\sim 3\sigma$ tension in $S_8$ between the $3\times2$pt constraints and Planck, with the exception of $w$CDM, where the $S_8$ tension is resolved. The tension in the $w$CDM case persists, however, when considering the joint $S_8$-$w$ parameter space. The joint flat $\Lambda$CDM CMB lensing and $3\times2$pt analysis is found to yield tight constraints on $\Omega_{\rm m}=0.307^{+0.008}_{-0.013}$, $\sigma_8=0.769^{+0.022}_{-0.010}$, and $S_8=0.779^{+0.013}_{-0.013}$.

Sunday, November 1, 2020

Day 1774

Wednesday, Thursday, Friday.


2010.13845
Disk dominated galaxies retain their shapes below $z=1.0$
Hoffmann, Laigle, Chisari, et al

The high abundance of disk galaxies without a large central bulge challenges predictions of current hydrodynamic simulations of galaxy formation. We aim to shed light on the formation of these objects by studying the redshift and mass dependence of their 3D shape distribution in the COSMOS galaxy survey. This distribution is inferred from the observed distribution of 2D shapes, using a reconstruction method which we test using hydrodynamic simulations. We find a moderate bias for the inferred average disk circularity and relative thickness with respect to the disk radius, but a large bias on the dispersion of these quantities. Applying the 3D shape reconstruction method on COSMOS data, we find no significant dependence of the inferred 3D shape distribution on redshift. The relative disk thickness shows a significant mass dependence which can be accounted for by the scaling of disk radius with galaxy mass. We conclude that the shapes of disk dominated galaxies are overall not subject to disruptive merging or feedback events below redshift $z=1.0$. This favours a scenario where these disks form early and subsequently undergo a tranquil evolution in isolation. In addition, our study shows that the observed 2D shapes of disk dominated galaxies can be well fitted using an ellipsoidal model for the galaxy 3D morphology combined with a Gaussian model for the 3D axes ratio distribution, confirming findings from similar work reported in the literature. Such an approach allows to build realistic mock catalogs with intrinsic galaxy shapes that will be essential for the study of intrinsic galaxy alignment as a contaminant of weak lensing surveys.


2010.13909
Modeling of magneto-rotational stellar evolution I. Method and first applications
Takahashi, Langer

While magnetic fields have long been considered to be important for the evolution of magnetic non-degenerate stars and compact stars, it has become clear in recent years that actually all of the stars are deeply affected. This is particularly true regarding their internal angular momentum distribution, but magnetic fields may also influence internal mixing processes and even the fate of the star. We propose a new framework for stellar evolution simulations, in which the interplay between magnetic field, rotation, mass loss, and changes in the stellar density and temperature distributions are treated self-consistently. For average large-scale stellar magnetic fields which are symmetric to the axis of rotation of the star, we derive 1D evolution equations for the toroidal and poloidal components from the mean-field MHD equation by applying Alfven's theorem, and a conservative form of the angular momentum transfer due to the Lorentz force is formulated. We implement our formalism into a numerical stellar evolution code and simulate the magneto-rotational evolution of 1.5 M$_\odot$ stars. The Lorentz force aided by the $\Omega$ effect imposes torsional Alfven waves propagating through the magnetized medium, leading to near-rigid rotation within the Alfven timescale. Our models with different initial spins and B-fields can reproduce the main observed properties of Ap/Bp stars. Calculations continued to the red-giant regime show a pronounced core-envelope coupling, which reproduces the core and surface rotation periods determined by asteroseismic observations.


2010.13992
Precision of the Chinese Space Station Telescope (CSST) stellar radial velocities
Sun, Deng, Yuan

The Chinese Space Station Telescope (CSST) spectroscopic survey plans to deliver high-quality low-resolution ($R > 200$) slitless spectra for hundreds of millions of targets down to a limiting magnitude of about 21 mag, covering a large survey area (17500 deg$^2$) and a wide wavelength range (255-1000 nm by 3 bands GU, GV, and GI). In this work, we use empirical spectra of the Next Generation Spectral Library to simulate the CSST stellar spectra at $R = 250$, and investigate their capabilities in measuring radial velocities. We find that velocity uncertainties depend strongly on effective temperature, weakly on metallicity for only FGK stars, and hardly on surface gravity. It is possible to deliver stellar radial velocities to a precision of about $3 \,\mathrm{km}\,\mathrm{s}^{-1}$ for AFGKM stars, and about $10 \,\mathrm{km}\,\mathrm{s}^{-1}$ for OB stars, at signal-to-noise ratio (SNR) of 100. Velocity uncertainties using single GU/GV/GI band spectra are also explored. Given the same SNR, the GU band performs best, the GV band the second best, and then the GI band. The effects of spectral normalization and imperfect template on velocity measurements are investigated and found to be very weak. The uncertainties caused by wavelength calibration are considered and found to be moderate. Given the possible precision of radial velocities, the CSST spectroscopic survey can enable interesting science such as searching for hyper-velocity stars. Limitations of our results are also discussed.


2010.14305
No phosphine in the atmosphere of Venus
Villanueva, et al

The detection of phosphine (PH3) has been recently reported in the atmosphere of Venus employing mm-wave radio observations (Greaves et at. 2020). We here demonstrate that the observed PH3 feature with JCMT can be fully explained employing plausible mesospheric SO2 abundances (~100 ppbv as per the SO2 profile given in their figure 9), while the identification of PH3 in the ALMA data should be considered invalid due to severe baseline calibration issues. We demonstrate this by independently calibrating and analyzing the ALMA data using different interferometric analysis tools, in which we observe no PH3 in all cases. Furthermore, for any PH3 signature to be produced in either ALMA or JCMT spectra, PH3 needs to present at altitudes above 70 km, in stark disagreement with their photochemical network. We ultimately conclude that this detection of PH3 in the atmosphere of Venus is incorrect and invite the Greaves et al. team to revise their work and consider a correction or retraction of their original report.


2010.15123
Measuring the mass and concentration of dark matter haloes from the velocity dispersion profile of their stars
Bose, Loeb

We use the IllustrisTNG (TNG) cosmological, hydrodynamical simulations of galaxy formation to measure the velocity dispersion profiles of dark matter and star particles in Milky Way-mass, galaxy group, and cluster-scale dark matter halos. The mean profile calculated from both dark and luminous tracers are similar in shape, exhibiting a large degree of halo-to-halo scatter around the average profile. The so-called "splashback" radius demarcates the outer boundary of the halo, and manifests as a kink in the velocity dispersion profile, located on average between $\sim 1.0-1.5r_{200m}$, where $r_{200m}$ is the radius within which the enclosed density of the halo equals 200 times the mean background density of the universe at that redshift. Interestingly, we find that this location may also be identified as the radius at which the (stacked) velocity dispersion profile drops to 60% of its peak value (for line-of-sight motions of stellar and dark matter particles in TNG halos). We further show that the scatter in the velocity dispersion profiles may be attributed to the variations in the assembly history of the host halos. In particular, this segregates the profile into two regimes: one within $\sim0.1r_{200m}$, where the scatter in the velocity dispersion within is set by the early assembly history of the halo, and the other beyond this radius where the scatter in the velocity dispersion is influenced more strongly by its late-time assembly. Finally, we show that a two-parameter model can be used to fit the measured velocity dispersion profiles and the fit parameters can be related directly to two fundamental halo properties: mass and concentration. We describe a simple model which allows us to express the stellar velocity dispersion profile in terms of the mass and concentration of the host halo as the only free parameters.


2010.15188
The statistical reliability of 267 GHz JCMT observations of Venus
Thompson

In the light of the recent announcement of the discovery of the potential biosignature phosphine in the atmosphere of Venus I present an independent reanalysis of the original JCMT data to assess the statistical reliability of the detection. Two line detection methods are explored, low order polynomial fits and higher order multiple polynomial fits. It is found that, similar to other reanalyses of ALMA Venus spectra, the polynomial fitting process results in false positive detections in the JCMT spectrum. Furthermore, a non-parametric bootstrap analysis reveals that neither line detection method is able to recover a statistically significant detection. There is thus no significant evidence for phosphine absorption in the JCMT Venus spectra.