The 100 most-cited articles in Astronomy and Astrophysics, 2020

We present cosmological parameter results from the final full-mission Planck measurements of the cosmic microwave background (CMB) anisotropies, combining information from the temperature and polarization maps and the lensing reconstruction. Compared to the 2015 results, improved measurements of large-scale polarization allow the reionization optical depth to be measured with higher precision, leading to significant gains in the precision of other correlated parameters. Improved modelling of the small-scale polarization leads to more robust constraints on many parameters, with residual modelling uncertainties estimated to affect them only at the 0.5 σ level. We find good consistency with the standard spatially-flat 6-parameter ΛCDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted “base ΛCDM” in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis gives dark matter density Ω c h 2 = 0.120 ± 0.001, baryon density Ω b h 2 = 0.0224 ± 0.0001, scalar spectral index n s = 0.965 ± 0.004, and optical…

Context. We present the early installment of the third Gaia data release, Gaia EDR3, consisting of astrometry and photometry for 1.8 billion sources brighter than magnitude 21, complemented with the list of radial velocities from Gaia DR2. Aims. A summary of the contents of Gaia EDR3 is presented, accompanied by a discussion on the differences with respect to Gaia DR2 and an overview of the main limitations which are present in the survey. Recommendations are made on the responsible use of Gaia EDR3 results. Methods. The raw data collected with the Gaia instruments during the first 34 months of the mission have been processed by the Gaia Data Processing and Analysis Consortium and turned into this early third data release, which represents a major advance with respect to Gaia DR2 in terms of astrometric and photometric precision, accuracy, and homogeneity. Results. Gaia EDR3 contains celestial positions and the apparent brightness in G for approximately 1.8 billion sources. For 1.5 billion of those sources, parallaxes, proper motions, and the ( G BP − G RP ) colour are also…

We describe the legacy Planck cosmic microwave background (CMB) likelihoods derived from the 2018 data release. The overall approach is similar in spirit to the one retained for the 2013 and 2015 data release, with a hybrid method using different approximations at low ( ℓ < 30) and high ( ℓ ≥ 30) multipoles, implementing several methodological and data-analysis refinements compared to previous releases. With more realistic simulations, and better correction and modelling of systematic effects, we can now make full use of the CMB polarization observed in the High Frequency Instrument (HFI) channels. The low-multipole EE cross-spectra from the 100 GHz and 143 GHz data give a constraint on the ΛCDM reionization optical-depth parameter τ to better than 15% (in combination with the TT low- ℓ data and the high- ℓ temperature and polarization data), tightening constraints on all parameters with posterior distributions correlated with τ . We also update the weaker constraint on τ from the joint TEB likelihood using the Low Frequency Instrument (LFI) channels, which was used in 2015 as part…

Context. Gaia Early Data Release 3 ( Gaia EDR3) contains results for 1.812 billion sources in the magnitude range G = 3–21 based on observations collected by the European Space Agency Gaia satellite during the first 34 months of its operational phase. Aims. We describe the input data, the models, and the processing used for the astrometric content of Gaia EDR3, as well as the validation of these results performed within the astrometry task. Methods. The processing broadly followed the same procedures as for Gaia DR2, but with significant improvements to the modelling of observations. For the first time in the Gaia data processing, colour-dependent calibrations of the line- and point-spread functions have been used for sources with well-determined colours from DR2. In the astrometric processing these sources obtained five-parameter solutions, whereas other sources were processed using a special calibration that allowed a pseudocolour to be estimated as the sixth astrometric parameter. Compared with DR2, the astrometric calibration models have been extended, and the spin-related…

Aims. Solar Orbiter, the first mission of ESA’s Cosmic Vision 2015–2025 programme and a mission of international collaboration between ESA and NASA, will explore the Sun and heliosphere from close up and out of the ecliptic plane. It was launched on 10 February 2020 04:03 UTC from Cape Canaveral and aims to address key questions of solar and heliospheric physics pertaining to how the Sun creates and controls the Heliosphere, and why solar activity changes with time. To answer these, the mission carries six remote-sensing instruments to observe the Sun and the solar corona, and four in-situ instruments to measure the solar wind, energetic particles, and electromagnetic fields. In this paper, we describe the science objectives of the mission, and how these will be addressed by the joint observations of the instruments onboard. Methods. The paper first summarises the mission-level science objectives, followed by an overview of the spacecraft and payload. We report the observables and performance figures of each instrument, as well as the trajectory design. This is followed by a summary…

Context. Gaia Early Data Release 3 ( Gaia EDR3) gives trigonometric parallaxes for nearly 1.5 billion sources. Inspection of the EDR3 data for sources identified as quasars reveals that their parallaxes are biased, that is, they are systematically offset from the expected distribution around zero, by a few tens of microarcseconds. Aims. We attempt to map the main dependences of the parallax bias in EDR3. In principle, this could provide a recipe for correcting the EDR3 parallaxes. Methods. Quasars provide the most direct way for estimating the parallax bias for faint sources. In order to extend this to brighter sources and a broader range of colours, we used differential methods based on physical pairs (binaries) and sources in the Large Magellanic Cloud. The functional forms of the dependences were explored by mapping the systematic differences between EDR3 and DR2 parallaxes. Results. The parallax bias is found to depend in a non-trivial way on (at least) the magnitude, colour, and ecliptic latitude of the source. Different dependences apply to the five- and six-parameter…

eROSITA (extended ROentgen Survey with an Imaging Telescope Array) is the primary instrument on the Spectrum-Roentgen-Gamma (SRG) mission, which was successfully launched on July 13, 2019, from the Baikonour cosmodrome. After the commissioning of the instrument and a subsequent calibration and performance verification phase, eROSITA started a survey of the entire sky on December 13, 2019. By the end of 2023, eight complete scans of the celestial sphere will have been performed, each lasting six months. At the end of this program, the eROSITA all-sky survey in the soft X-ray band (0.2–2.3 keV) will be about 25 times more sensitive than the ROSAT All-Sky Survey, while in the hard band (2.3–8 keV) it will provide the first ever true imaging survey of the sky. The eROSITA design driving science is the detection of large samples of galaxy clusters up to redshifts z > 1 in order to study the large-scale structure of the universe and test cosmological models including Dark Energy. In addition, eROSITA is expected to yield a sample of a few million AGNs, including obscured objects,…

Context. Gaia Early Data Release 3 ( Gaia EDR3) contains astrometry and photometry results for about 1.8 billion sources based on observations collected by the European Space Agency Gaia satellite during the first 34 months of its operational phase. Aims. In this paper, we focus on the photometric content, describing the input data, the algorithms, the processing, and the validation of the results. Particular attention is given to the quality of the data and to a number of features that users may need to take into account to make the best use of the Gaia EDR3 catalogue. Methods. The processing broadly followed the same procedure as for Gaia DR2, but with significant improvements in several aspects of the blue and red photometer (BP and RP) preprocessing and in the photometric calibration process. In particular, the treatment of the BP and RP background has been updated to include a better estimation of the local background, and the detection of crowding effects has been used to exclude affected data from the calibrations. The photometric calibration models have also been updated to…

We present a joint cosmological analysis of weak gravitational lensing observations from the Kilo-Degree Survey (KiDS-1000), with redshift-space galaxy clustering observations from the Baryon Oscillation Spectroscopic Survey (BOSS) and galaxy-galaxy lensing observations from the overlap between KiDS-1000, BOSS, and the spectroscopic 2-degree Field Lensing Survey. This combination of large-scale structure probes breaks the degeneracies between cosmological parameters for individual observables, resulting in a constraint on the structure growth parameter S 8 = σ 8 √(Ω m /0.3) = 0.766 −0.014 +0.020 , which has the same overall precision as that reported by the full-sky cosmic microwave background observations from Planck . The recovered S 8 amplitude is low, however, by 8.3 ± 2.6% relative to Planck . This result builds from a series of KiDS-1000 analyses where we validate our methodology with variable depth mock galaxy surveys, our lensing calibration with image simulations and null-tests, and our optical-to-near-infrared redshift calibration with multi-band mock catalogues and a…

We present cosmological constraints from a cosmic shear analysis of the fourth data release of the Kilo-Degree Survey (KiDS-1000), which doubles the survey area with nine-band optical and near-infrared photometry with respect to previous KiDS analyses. Adopting a spatially flat standard cosmological model, we find S 8 = σ 8 (Ω m /0.3) 0.5 = 0.759 −0.021 +0.024 for our fiducial analysis, which is in 3 σ tension with the prediction of the Planck Legacy analysis of the cosmic microwave background. We compare our fiducial COSEBIs (Complete Orthogonal Sets of E/B-Integrals) analysis with complementary analyses of the two-point shear correlation function and band power spectra, finding the results to be in excellent agreement. We investigate the sensitivity of all three statistics to a number of measurement, astrophysical, and modelling systematics, finding our S 8 constraints to be robust and dominated by statistical errors. Our cosmological analysis of different divisions of the data passes the Bayesian internal consistency tests, with the exception of the second tomographic bin. As…

We analyse the Planck full-mission cosmic microwave background (CMB) temperature and E -mode polarization maps to obtain constraints on primordial non-Gaussianity (NG). We compare estimates obtained from separable template-fitting, binned, and optimal modal bispectrum estimators, finding consistent values for the local, equilateral, and orthogonal bispectrum amplitudes. Our combined temperature and polarization analysis produces the following final results: f NL local = −0.9 ± 5.1; f NL equil = −26 ± 47; and f NL ortho = −38 ± 24 (68% CL, statistical). These results include low-multipole (4 ≤ ℓ < 40) polarization data that are not included in our previous analysis. The results also pass an extensive battery of tests (with additional tests regarding foreground residuals compared to 2015), and they are stable with respect to our 2015 measurements (with small fluctuations, at the level of a fraction of a standard deviation, which is consistent with changes in data processing). Polarization-only bispectra display a significant improvement in robustness; they can now be used…

The star S2 orbiting the compact radio source Sgr A* is a precision probe of the gravitational field around the closest massive black hole (candidate). Over the last 2.7 decades we have monitored the star’s radial velocity and motion on the sky, mainly with the SINFONI and NACO adaptive optics (AO) instruments on the ESO VLT, and since 2017, with the four-telescope interferometric beam combiner instrument GRAVITY. In this Letter we report the first detection of the General Relativity (GR) Schwarzschild Precession (SP) in S2’s orbit. Owing to its highly elliptical orbit ( e = 0.88), S2’s SP is mainly a kink between the pre-and post-pericentre directions of motion ≈±1 year around pericentre passage, relative to the corresponding Kepler orbit. The superb 2017−2019 astrometry of GRAVITY defines the pericentre passage and outgoing direction. The incoming direction is anchored by 118 NACO-AO measurements of S2’s position in the infrared reference frame, with an additional 75 direct measurements of the S2-Sgr A* separation during bright states (“flares”) of Sgr A*. Our 14-parameter model…

Context.The large astrometric and photometric survey performed by theGaiamission allows for a panoptic view of the Galactic disc and its stellar cluster population. Hundreds of stellar clusters were only discovered after the latestGaiadata release (DR2) and have yet to be characterised. Aims.Here we make use of the deep and homogeneousGaiaphotometry down toG= 18 to estimate the distance, age, and interstellar reddening for about 2000 stellar clusters identified withGaiaDR2 astrometry. We use these objects to study the structure and evolution of the Galactic disc. Methods.We relied on a set of objects with well-determined parameters in the literature to train an artificial neural network to estimate parameters from theGaiaphotometry of cluster members and their mean parallax. Results.We obtain reliable parameters for 1867 clusters. Our catalogue confirms the relative lack of old stellar clusters in the inner disc (with a few notable exceptions). We also quantify and discuss the variation of scale height with cluster age, and we detect the Galactic warp in the distribution of old…

Context. ESPRESSO is the new high-resolution spectrograph of ESO’s Very Large Telescope (VLT). It was designed for ultra-high radial-velocity (RV) precision and extreme spectral fidelity with the aim of performing exoplanet research and fundamental astrophysical experiments with unprecedented precision and accuracy. It is able to observe with any of the four Unit Telescopes (UTs) of the VLT at a spectral resolving power of 140 000 or 190 000 over the 378.2 to 788.7 nm wavelength range; it can also observe with all four UTs together, turning the VLT into a 16 m diameter equivalent telescope in terms of collecting area while still providing a resolving power of 70 000. Aims. We provide a general description of the ESPRESSO instrument, report on its on-sky performance, and present our Guaranteed Time Observation (GTO) program along with its first results. Methods. ESPRESSO was installed on the Paranal Observatory in fall 2017. Commissioning (on-sky testing) was conducted between December 2017 and September 2018. The instrument saw its official start of operations on October 1, 2018,…

Context. The Extreme Ultraviolet Imager (EUI) is part of the remote sensing instrument package of the ESA/NASA Solar Orbiter mission that will explore the inner heliosphere and observe the Sun from vantage points close to the Sun and out of the ecliptic. Solar Orbiter will advance the “connection science” between solar activity and the heliosphere. Aims. With EUI we aim to improve our understanding of the structure and dynamics of the solar atmosphere, globally as well as at high resolution, and from high solar latitude perspectives. Methods. The EUI consists of three telescopes, the Full Sun Imager and two High Resolution Imagers, which are optimised to image in Lyman- α and EUV (17.4 nm, 30.4 nm) to provide a coverage from chromosphere up to corona. The EUI is designed to cope with the strong constraints imposed by the Solar Orbiter mission characteristics. Limited telemetry availability is compensated by state-of-the-art image compression, onboard image processing, and event selection. The imposed power limitations and potentially harsh radiation environment necessitate the use…

Aims. The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. The estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on various individual methodologies and numerical implementations, which were developed for different observational probes and for the combination thereof. In this paper we present validated forecasts, which combine both theoretical and observational ingredients for different cosmological probes. This work is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. Methods. We describe in detail the methods adopted for Fisher matrix forecasts, which were applied to galaxy clustering, weak lensing, and the combination thereof. We estimated the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on…

Context. The third Gaia data release is published in two stages. The early part, Gaia EDR3, gives very precise astrometric and photometric properties for nearly two billion sources together with seven million radial velocities from Gaia DR2. The full release, Gaia DR3, will add radial velocities, spectra, light curves, and astrophysical parameters for a large subset of the sources, as well as orbits for solar system objects. Aims. Before the publication of the catalogue, many different data items have undergone dedicated validation processes. The goal of this paper is to describe the validation results in terms of completeness, accuracy, and precision for the Gaia EDR3 data and to provide recommendations for the use of the catalogue data. Methods. The validation processes include a systematic analysis of the catalogue contents to detect anomalies, either individual errors or statistical properties, using statistical analysis and comparisons to the previous release as well as to external data and to models. Results. Gaia EDR3 represents a major step forward, compared to Gaia DR2, in…

The H0LiCOW collaboration inferred via strong gravitational lensing time delays a Hubble constant value of H 0 = 73.3 −1.8 +1.7 km s −1 Mpc −1 , describing deflector mass density profiles by either a power-law or stars (constant mass-to-light ratio) plus standard dark matter halos. The mass-sheet transform (MST) that leaves the lensing observables unchanged is considered the dominant source of residual uncertainty in H 0 . We quantify any potential effect of the MST with a flexible family of mass models, which directly encodes it, and they are hence maximally degenerate with H 0 . Our calculation is based on a new hierarchical Bayesian approach in which the MST is only constrained by stellar kinematics. The approach is validated on mock lenses, which are generated from hydrodynamic simulations. We first applied the inference to the TDCOSMO sample of seven lenses, six of which are from H0LiCOW, and measured H 0 = 74.5 −6.1 +5.6 km s −1 Mpc −1 . Secondly, in order to further constrain the deflector mass density profiles, we added imaging and spectroscopy for a set of 33 strong…

Context. Sky surveys produce enormous quantities of data on extensive regions of the sky. The easiest way to access this information is through catalogues of standardised data products. XMM-Newton has been surveying the sky in the X-ray, ultra-violet, and optical bands for 20 years. Aims. The XMM-Newton Survey Science Centre has been producing standardised data products and catalogues to facilitate access to the serendipitous X-ray sky. Methods. Using improved calibration and enhanced software, we re-reduced all of the 14 041 XMM-Newton X-ray observations, of which 11 204 observations contained data with at least one detection and with these we created a new, high quality version of the XMM-Newton serendipitous source catalogue, 4XMM-DR9. Results. 4XMM-DR9 contains 810 795 detections down to a detection significance of 3 σ , of which 550 124 are unique sources, which cover 1152 degrees 2 (2.85%) of the sky. Filtering 4XMM-DR9 to retain only the cleanest sources with at least a 5 σ detection significance leaves 433 612 detections. Of these detections, 99.6% have no pileup.…

All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported following the O1/O2 runs have near-zero effective spins. There are only three potential explanations for this. If the BH spin magnitudes are large, then: (i) either both BH spin vectors must be nearly in the orbital plane or (ii) the spin angular momenta of the BHs must be oppositely directed and similar in magnitude. Then there is also the possibility that (iii) the BH spin magnitudes are small. We consider the third hypothesis within the framework of the classical isolated binary evolution scenario of the BH-BH merger formation. We test three models of angular momentum transport in massive stars: a mildly efficient transport by meridional currents (as employed in the Geneva code), an efficient transport by the Tayler-Spruit magnetic dynamo (as implemented in the MESA code), and a very-efficient transport (as proposed by Fuller et al.) to calculate natal BH spins. We allow for binary evolution to increase the BH spins through accretion and account for the potential spin-up of stars through tidal interactions.…

The processing of raw data from modern astronomical instruments is often carried out nowadays using dedicated software, known as pipelines, largely run in automated operation. In this paper we describe the data reduction pipeline of the Multi Unit Spectroscopic Explorer (MUSE) integral field spectrograph operated at the ESO Paranal Observatory. This spectrograph is a complex machine: it records data of 1152 separate spatial elements on detectors in its 24 integral field units. Efficiently handling such data requires sophisticated software with a high degree of automation and parallelization. We describe the algorithms of all processing steps that operate on calibrations and science data in detail, and explain how the raw science data is transformed into calibrated datacubes. We finally check the quality of selected procedures and output data products, and demonstrate that the pipeline provides datacubes ready for scientific analysis.

Context. The AGN bolometric correction is a key element for understanding black hole (BH) demographics and computing accurate BH accretion histories from AGN luminosities. However, current estimates still differ from each other by up to a factor of two to three, and rely on extrapolations at the lowest and highest luminosities. Aims. Here we revisit this fundamental question by presenting general hard X-ray ( K X ) and optical ( K O ) bolometric corrections, computed by combining several AGN samples spanning the widest (about 7 dex) luminosity range ever used for this type of studies. Methods. We analysed a total of ∼1000 type 1 and type 2 AGN for which we performed a dedicated SED-fitting. Results. We provide a bolometric correction separately for type 1 and type 2 AGN; the two bolometric corrections agree in the overlapping luminosity range. Based on this we computed for the first time a universal bolometric correction for the whole AGN sample (both type 1 and type 2). We found that K X is fairly constant at log( L BOL / L ⊙ ) < 11, while it increases up to about one order of…

Aims. We produce a clean and well-characterised catalogue of objects within 100 pc of the Sun from the Gaia Early Data Release 3. We characterise the catalogue through comparisons to the full data release, external catalogues, and simulations. We carry out a first analysis of the science that is possible with this sample to demonstrate its potential and best practices for its use. Methods. Theselection of objects within 100 pc from the full catalogue used selected training sets, machine-learning procedures, astrometric quantities, and solution quality indicators to determine a probability that the astrometric solution is reliable. The training set construction exploited the astrometric data, quality flags, and external photometry. For all candidates we calculated distance posterior probability densities using Bayesian procedures and mock catalogues to define priors. Any object with reliable astrometry and a non-zero probability of being within 100 pc is included in the catalogue. Results. We have produced a catalogue of 331 312 objects that we estimate contains at least 92% of stars…

The Solar Orbiter mission seeks to make connections between the physical processes occurring at the Sun or in the solar corona and the nature of the solar wind created by those processes which is subsequently observed at the spacecraft. The mission also targets physical processes occurring in the solar wind itself during its journey from its source to the spacecraft. To meet the specific mission science goals, Solar Orbiter will be equipped with both remote-sensing and in-situ instruments which will make unprecedented measurements of the solar atmosphere and the inner heliosphere. A crucial set of measurements will be provided by the Solar Wind Analyser (SWA) suite of instruments. This suite consists of an Electron Analyser System (SWA-EAS), a Proton and Alpha particle Sensor (SWA-PAS), and a Heavy Ion Sensor (SWA-HIS) which are jointly served by a central control and data processing unit (SWA-DPU). Together these sensors will measure and categorise the vast majority of thermal and suprathermal ions and electrons in the solar wind and determine the abundances and charge states of…

We present a new set of solar metallicity atmosphere and evolutionary models for very cool brown dwarfs and self-luminous giant exoplanets, which we term ATMO 2020. Atmosphere models are generated with our state-of-the-art 1D radiative-convective equilibrium code ATMO , and are used as surface boundary conditions to calculate the interior structure and evolution of 0.001–0.075 M ⊙ objects. Our models include several key improvements to the input physics used in previous models available in the literature. Most notably, the use of a new H–He equation of state including ab initio quantum molecular dynamics calculations has raised the mass by ~1−2% at the stellar–substellar boundary and has altered the cooling tracks around the hydrogen and deuterium burning minimum masses. A second key improvement concerns updated molecular opacities in our atmosphere model ATMO , which now contains significantly more line transitions required to accurately capture the opacity in these hot atmospheres. This leads to warmer atmospheric temperature structures, further changing the cooling curves and…

Aims. The Spectrometer Telescope for Imaging X-rays (STIX) on Solar Orbiter is a hard X-ray imaging spectrometer, which covers the energy range from 4 to 150 keV. STIX observes hard X-ray bremsstrahlung emissions from solar flares and therefore provides diagnostics of the hottest (⪆10 MK) flare plasma while quantifying the location, spectrum, and energy content of flare-accelerated nonthermal electrons. Methods. To accomplish this, STIX applies an indirect bigrid Fourier imaging technique using a set of tungsten grids (at pitches from 0.038 to 1 mm) in front of 32 coarsely pixelated CdTe detectors to provide information on angular scales from 7 to 180 arcsec with 1 keV energy resolution (at 6 keV). The imaging concept of STIX has intrinsically low telemetry and it is therefore well-suited to the limited resources available to the Solar Orbiter payload. To further reduce the downlinked data volume, STIX data are binned on board into 32 selectable energy bins and dynamically-adjusted time bins with a typical duration of 1 s during flares. Results. Through hard X-ray diagnostics, STIX…

We present the NPIPE processing pipeline, which produces calibrated frequency maps in temperature and polarization from data from the Planck Low Frequency Instrument (LFI) and High Frequency Instrument (HFI) using high-performance computers. NPIPE represents a natural evolution of previous Planck analysis efforts, and combines some of the most powerful features of the separate LFI and HFI analysis pipelines. For example, following the LFI 2018 processing procedure, NPIPE uses foreground polarization priors during the calibration stage in order to break scanning-induced degeneracies. Similarly, NPIPE employs the HFI 2018 time-domain processing methodology to correct for bandpass mismatch at all frequencies. In addition, NPIPE introduces several improvements, including, but not limited to: inclusion of the 8% of data collected during repointing manoeuvres; smoothing of the LFI reference load data streams; in-flight estimation of detector polarization parameters; and construction of maximally independent detector-set split maps. For component-separation purposes, important improvements…

A new realization of the International Celestial Reference Frame (ICRF) is presented based on the work achieved by a working group of the International Astronomical Union (IAU) mandated for this purpose. This new realization follows the initial realization of the ICRF completed in 1997 and its successor, ICRF2, adopted as a replacement in 2009. The new frame, referred to as ICRF3, is based on nearly 40 years of data acquired by very long baseline interferometry at the standard geodetic and astrometric radio frequencies (8.4 and 2.3 GHz), supplemented with data collected at higher radio frequencies (24 GHz and dual-frequency 32 and 8.4 GHz) over the past 15 years. State-of-the-art astronomical and geophysical modeling has been used to analyze these data and derive source positions. The modeling integrates, for the first time, the effect of the galactocentric acceleration of the solar system (directly estimated from the data) which, if not considered, induces significant deformation of the frame due to the data span. The new frame includes positions at 8.4 GHz for 4536 extragalactic…

The magnetometer instrument on the Solar Orbiter mission is designed to measure the magnetic field local to the spacecraft continuously for the entire mission duration. The need to characterise not only the background magnetic field but also its variations on scales from far above to well below the proton gyroscale result in challenging requirements on stability, precision, and noise, as well as magnetic and operational limitations on both the spacecraft and other instruments. The challenging vibration and thermal environment has led to significant development of the mechanical sensor design. The overall instrument design, performance, data products, and operational strategy are described.

The Atacama Large Millimeter Array (ALMA) Large Program to INvestigate [CII] at Early times (ALPINE) targets the [CII] 158 μ m line and the far-infrared continuum in 118 spectroscopically confirmed star-forming galaxies between z = 4.4 and z = 5.9. It represents the first large [CII] statistical sample built in this redshift range. We present details regarding the data processing and the construction of the catalogs. We detected 23 of our targets in the continuum. To derive accurate infrared luminosities and obscured star formation rates (SFRs), we measured the conversion factor from the ALMA 158 μ m rest-frame dust continuum luminosity to the total infrared luminosity ( L IR ) after constraining the dust spectral energy distribution by stacking a photometric sample similar to ALPINE in ancillary single-dish far-infrared data. We found that our continuum detections have a median L IR of 4.4 × 10 11 L ⊙ . We also detected 57 additional continuum sources in our ALMA pointings. They are at a lower redshift than the ALPINE targets, with a mean photometric redshift of 2.5 ± 0.2. We…

Context. It has become clear in recent years that binarity plays a crucial role in many aspects of planetary nebulae (PNe), particularly with regard to the striking morphologies they exhibit. To date, there are nearly 60 known binary central stars of PNe (bCSPNe). However, both theory and observation indicate that this figure represents only the tip of the iceberg, with the Galactic PN population hosting orders of magnitude more stars. Aims. We are involved in a search for new bCSPNe with the aim of enhancing the statistical validation of the key role of binarity in the formation and shaping of PNe. New discoveries of bCSPNe and their characterization carry important implications not only for understanding PN evolution, but also for studying binary evolution and the common-envelope phase, which is still poorly understood. Methods. We used data from the TESS satellite to search for variability in the eight CSPNe that belong to the two-minute cadence of preselected targets in Cycle 1, with their available pipeline-extracted light curves. We identified strong periodicities and analysed…

Planck 2018 resultsIn this volume, we proudly present a special feature on the results from the ESA Planck mission, based on data released by ESA and the Planck Collaboration in July 2018.This 2018 data release has significantly lower systematic residuals for both Planck instruments, LFI and HFI, and a more accurate photometric calibration for HFI.These calibration improvements are most significant over the largest angular scales and for the polarized emission.The resulting frequency maps were used to separate the diffuse sky emission into maps of the cosmic microwave background, the thermal Sunyaev-Zeldovich effect, the cosmic infrared background, and the Galactic thermal dust, carbon monoxide, anomalous dust, free-free, and synchrotron emission.The 12 articles in the special issue describe the released data products and present scientific results extracted by the Planck Collaboration from this data.The six-parameter ΛCDM model continues to provide an excellent fit to the cosmic microwave background data at high and low redshift.Planck measures five of those six parameters to…

Context. After years of scientific progress, the origin of stellar binary black holes is still a great mystery. Several formation channels for merging black holes have been proposed in the literature. As more merger detections are expected with future gravitational-wave observations, population synthesis studies can help to distinguish between them. Aims. We study the formation of coalescing binary black holes via the evolution of isolated field binaries that go through the common envelope phase in order to obtain the combined distributions of observables such as black-hole spins, masses and cosmological redshifts of mergers. Methods. To achieve this aim, we used a hybrid technique that combines the parametric binary population synthesis code COMPAS with detailed binary evolution simulations performed with the MESA code. We then convolved our binary evolution calculations with the redshift- and metallicity-dependent star-formation rate and the selection effects of gravitational-wave detectors to obtain predictions of observable properties. Results. By assuming efficient angular…

Aims. We study the evolution of the rotation and the high energy X-ray, extreme ultraviolet (EUV), and Ly-α emission for F, G, K, and M dwarfs, with masses between 0.1 and 1.2 M⊙, and provide a freely available set of evolutionary tracks for use in planetary atmosphere studies. Methods. We develop a physical rotational evolution model constrained by observed rotation distributions in young stellar clusters. Using rotation, X-ray, EUV, and Ly-α measurements, we derive empirical relations for the dependences of high energy emission on stellar parameters. Our description of X-ray evolution is validated using measurements of X-ray distributions in young clusters. Results. A star’s X-ray, EUV, and Ly-α evolution is determined by its mass and initial rotation rate, with initial rotation being less important for lower mass stars. At all ages, solar mass stars are significantly more X-ray luminous than lower mass stars and stars that are born as rapid rotators remain highly active longer than those born as slow rotators. At all evolutionary stages, habitable zone planets receive higher…

The Low Frequency Array (LOFAR) is an ideal instrument to conduct deep extragalactic surveys. It has a large field of view and is sensitive to large-scale and compact emission. It is, however, very challenging to synthesize thermal noise limited maps at full resolution, mainly because of the complexity of the low-frequency sky and the direction dependent effects (phased array beams and ionosphere). In this first paper of a series, we present a new calibration and imaging pipeline that aims at producing high fidelity, high dynamic range images with LOFAR High Band Antenna data, while being computationally efficient and robust against the absorption of unmodeled radio emission. We apply this calibration and imaging strategy to synthesize deep images of the Boötes and Lockman Hole fields at ~150 MHz, totaling ~80 and ~100 h of integration, respectively, and reaching unprecedented noise levels at these low frequencies of ≲30 and ≲23μJy beam−1in the inner ~3 deg2. This approach is also being used to reduce the LOTSS-wide data for the second data release.

The ALMA-ALPINE [CII] survey is aimed at characterizing the properties of a sample of normal star-forming galaxies (SFGs). The ALMA Large Program to INvestigate (ALPINE) features 118 galaxies observed in the [CII]-158 μ m line and far infrared (FIR) continuum emission during the period of rapid mass assembly, right after the end of the HI reionization, at redshifts of 4 < z < 6. We present the survey science goals, the observational strategy, and the sample selection of the 118 galaxies observed with ALMA, with an average beam minor axis of about 0.85″, or ∼5 kpc at the median redshift of the survey. The properties of the sample are described, including spectroscopic redshifts derived from the UV-rest frame, stellar masses, and star-formation rates obtained from a spectral energy distribution (SED) fitting. The observed properties derived from the ALMA data are presented and discussed in terms of the overall detection rate in [CII] and FIR continuum, with the observed signal-to-noise distribution. The sample is representative of the SFG population in the main sequence at these…

The SpHere INfrared Exoplanet (SHINE) project is a 500-star survey performed with SPHERE on the Very Large Telescope for the purpose of directly detecting new substellar companions and understanding their formation and early evolution. Here we present an initial statistical analysis for a subsample of 150 stars spanning spectral types from B to M that are representative of the full SHINE sample. Our goal is to constrain the frequency of substellar companions with masses between 1 and 75 MJup and semimajor axes between 5 and 300 au. For this purpose, we adopt detection limits as a function of angular separation from the survey data for all stars converted into mass and projected orbital separation using the BEX-COND-hot evolutionary tracks and known distance to each system. Based on the results obtained for each star and on the 13 detections in the sample, we use a Markov chain Monte Carlo tool to compare our observations to two different types of models. The first is a parametric model based on observational constraints, and the second type are numerical models that combine advanced…

Accurate distances to local molecular clouds are critical for understanding the star and planet formation process, yet distance measurements are often obtained inhomogeneously on a cloud-by-cloud basis. We have recently developed a method that combines stellar photometric data withGaiaDR2 parallax measurements in a Bayesian framework to infer the distances of nearby dust clouds to a typical accuracy of ∼5%. After refining the technique to target lower latitudes and incorporating deep optical data from DECam in the southern Galactic plane, we have derived a catalog of distances to molecular clouds in Reipurth (2008, Star Formation Handbook, Vols. I and II) which contains a large fraction of the molecular material in the solar neighborhood. Comparison with distances derived from maser parallax measurements towards the same clouds shows our method produces consistent distances with ≲10% scatter for clouds across our entire distance spectrum (150 pc−2.5 kpc). We hope this catalog of homogeneous distances will serve as a baseline for future work.

The existence of a radius valley in the Kepler size distribution stands as one of the most important observational constraints to understand the origin and composition of exoplanets with radii between those of Earth and Neptune. In this work we provide insights into the existence of the radius valley, first from a pure formation point of view and then from a combined formation-evolution model. We run global planet formation simulations including the evolution of dust by coagulation, drift, and fragmentation, and the evolution of the gaseous disc by viscous accretion and photoevaporation. A planet grows from a moon-mass embryo by either silicate or icy pebble accretion, depending on its position with respect to the water ice line. We include gas accretion, type I–II migration, and photoevaporation driven mass-loss after formation. We perform an extensive parameter study evaluating a wide range of disc properties and initial locations of the embryo. We find that due to the change in dust properties at the water ice line, rocky cores form typically with ∼3 M ⊕ and have a maximum mass…

Time-delay cosmography of lensed quasars has achieved 2.4% precision on the measurement of the Hubble constant, H 0 . As part of an ongoing effort to uncover and control systematic uncertainties, we investigate three potential sources: 1- stellar kinematics, 2- line-of-sight effects, and 3- the deflector mass model. To meet this goal in a quantitative way, we reproduced the H0LiCOW/SHARP/STRIDES (hereafter TDCOSMO) procedures on a set of real and simulated data, and we find the following. First, stellar kinematics cannot be a dominant source of error or bias since we find that a systematic change of 10% of measured velocity dispersion leads to only a 0.7% shift on H 0 from the seven lenses analyzed by TDCOSMO. Second, we find no bias to arise from incorrect estimation of the line-of-sight effects. Third, we show that elliptical composite (stars + dark matter halo), power-law, and cored power-law mass profiles have the flexibility to yield a broad range in H 0 values. However, the TDCOSMO procedures that model the data with both composite and power-law mass profiles are informative.…

Context. Clouds are ubiquitous in exoplanet atmospheres and they represent a challenge for the model interpretation of their spectra. When generating a large number of model spectra, complex cloud models often prove too costly numerically, whereas more efficient models may be overly simplified. Aims. We aim to constrain the atmospheric properties of the directly imaged planet HR 8799e with a free retrieval approach. Methods. We used our radiative transfer code petitRADTRANS for generating the spectra, which we coupled to the PyMultiNest tool. We added the effect of multiple scattering which is important for treating clouds. Two cloud model parameterizations are tested: the first incorporates the mixing and settling of condensates, the second simply parameterizes the functional form of the opacity. Results. In mock retrievals, using an inadequate cloud model may result in atmospheres that are more isothermal and less cloudy than the input. Applying our framework on observations of HR 8799e made with the GPI, SPHERE, and GRAVITY, we find a cloudy atmosphere governed by disequilibrium…

Context. Recent years have seen building evidence that planet formation starts early, in the first ~0.5 Myr. Studying the dust masses available in young disks enables us to understand the origin of planetary systems given that mature disks are lacking the solid material necessary to reproduce the observed exoplanetary systems, especially the massive ones. Aims. We aim to determine if disks in the embedded stage of star formation contain enough dust to explain the solid content of the most massive exoplanets. Methods. We use Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 (1.1–1.3 mm) continuum observations of embedded disks in the Perseus star-forming region together with Very Large Array (VLA) Ka -band (9 mm) data to provide a robust estimate of dust disk masses from the flux densities measured in the image plane. Results. We find a strong linear correlation between the ALMA and VLA fluxes, demonstrating that emission at both wavelengths is dominated by dust emission. For a subsample of optically thin sources, we find a median spectral index of 2.5 from which we derive…

Context. Open clusters are key targets for studies of Galaxy structure and evolution, and stellar physics. Since the Gaia data release 2 (DR2), the discovery of undetected clusters has shown that previous surveys were incomplete. Aims. Our aim is to exploit the Big Data capabilities of machine learning to detect new open clusters in Gaia DR2, and to complete the open cluster sample to enable further studies of the Galactic disc. Methods. We use a machine-learning based methodology to systematically search the Galactic disc for overdensities in the astrometric space and identify the open clusters using photometric information. First, we used an unsupervised clustering algorithm, DBSCAN, to blindly search for these overdensities in Gaia DR2 ( l , b , ϖ , μ α * , μ δ ), and then we used a deep learning artificial neural network trained on colour–magnitude diagrams to identify isochrone patterns in these overdensities, and to confirm them as open clusters. Results. We find 582 new open clusters distributed along the Galactic disc in the region | b | < 20°. We detect substructure in…

Ensemble studies of red-giant stars with exquisite asteroseismic (Kepler), spectroscopic (APOGEE), and astrometric (Gaia) constraints offer a novel opportunity to recast and address long-standing questions concerning the evolution of stars and of the Galaxy. Here, we infer masses and ages for nearly 5400 giants with availableKeplerlight curves and APOGEE spectra using the codePARAM, and discuss some of the systematics that may affect the accuracy of the inferred stellar properties. We then present patterns in mass, evolutionary state, age, chemical abundance, and orbital parameters that we deem robust against the systematic uncertainties explored. First, we look at age-chemical-abundances ([Fe/H] and [α/Fe]) relations. We find a dearth of young, metal-rich ([Fe/H] > 0.2) stars, and the existence of a significant population of old (8−9 Gyr), low-[α/Fe], super-solar metallicity stars, reminiscent of the age and metallicity of the well-studied open cluster NGC 6791. The age-chemo-kinematic properties of these stars indicate that efficient radial migration happens in the thin disc. We…

Aims. We present the detailed characterisation of a sample of 56 sources serendipitously detected in ALMA band 7 as part of the ALMA Large Program to INvestigate CII at Early Times (ALPINE). These sources, detected in COSMOS and ECDFS, have been used to derive the total infrared luminosity function (LF) and to estimate the cosmic star formation rate density (SFRD) up to z ≃ 6. Methods. We looked for counterparts of the ALMA sources in all the available multi-wavelength (from HST to VLA) and photometric redshift catalogues. We also made use of deeper UltraVISTA and Spitzer source lists and maps to identify optically dark sources with no matches in the public catalogues. We used the sources with estimated redshifts to derive the 250 μ m rest-frame and total infrared (8–1000 μ m) LFs from z ≃ 0.5 to 6. Results. Our ALMA blind survey (860 μ m flux density range: ∼0.3–12.5 mJy) allows us to further push the study of the nature and evolution of dusty galaxies at high- z , identifying luminous and massive sources to redshifts and faint luminosities never probed before by any far-infrared…

Context. We describe the atomic and molecular data that were used for the abundance analyses of FGK-type stars carried out within the Gaia -ESO Public Spectroscopic Survey in the years 2012 to 2019. The Gaia -ESO Survey is one among several current and future stellar spectroscopic surveys producing abundances for Milky-Way stars on an industrial scale. Aims. We present an unprecedented effort to create a homogeneous common line list, which was used by several abundance analysis groups using different radiative transfer codes to calculate synthetic spectra and equivalent widths. The atomic data are accompanied by quality indicators and detailed references to the sources. The atomic and molecular data are made publicly available at the CDS. Methods. In general, experimental transition probabilities were preferred but theoretical values were also used. Astrophysical gf -values were avoided due to the model-dependence of such a procedure. For elements whose lines are significantly affected by a hyperfine structure or isotopic splitting, a concerted effort has been made to collate the…

Following the recent discovery of X-ray quasi-periodic eruptions (QPEs) coming from the nucleus of the galaxy GSN 069, here we report on the detection of QPEs in the active galaxy named RX J1301.9+2747. QPEs are rapid and recurrent increases of the X-ray count-rate by more than one order of magnitude with respect to a stable quiescent level. During a XMM-Newton observation lasting 48 ks that was performed on 30 and 31 May 2019, three strong QPEs lasting about half an hour each were detected in the light curves of RX J1301.9+2747. The first two QPEs are separated by a longer recurrence time (about 20 ks) compared to the second and third (about 13 ks). This pattern is consistent with the alternating long-short recurrence times of the GSN 069 QPEs, although the difference between the consecutive recurrence times is significantly smaller in GSN 069. Longer X-ray observations will better clarify the temporal pattern of the QPEs in RX J1301.9+2747 and will allow a detailed comparison with GSN 069 to be performed. The X-ray spectral properties of QPEs in the two sources are remarkably…

Aims. We aim to study vertical settling and radial drift of dust in protoplanetary disks from a different perspective: an edge-on view. An estimation of the amplitude of settling and drift is highly relevant to understanding planet formation. Methods. We analyze a sample of 12 HST-selected edge-on protoplanetary disks (i.e., seen with high inclinations) for which the vertical extent of the emission layers can be constrained directly. We present ALMA high angular resolution continuum images (~0.1′′) of these disks at two wavelengths, 0.89 and 2.06 mm (respectively ALMA bands 7 and 4), supplemented with archival band 6 data (1.33 mm) where available. Results. Several sources show constant brightness profiles along their major axis with a steep drop at the outer edges. Two disks have central holes with additional compact continuum emission at the location of the central star. For most sources, the millimeter continuum emission is more compact than the scattered light, both in the vertical and radial directions. Six sources are resolved along their minor axis in at least one millimetric…

The isotropy of the late Universe and consequently of the X-ray galaxy cluster scaling relations is an assumption greatly used in astronomy. However, within the last decade, many studies have reported deviations from isotropy when using various cosmological probes; a definitive conclusion has yet to be made. New, effective and independent methods to robustly test the cosmic isotropy are of crucial importance. In this work, we use such a method. Specifically, we investigate the directional behavior of the X-ray luminosity-temperature ( L X – T ) relation of galaxy clusters. A tight correlation is known to exist between the luminosity and temperature of the X-ray-emitting intracluster medium of galaxy clusters. While the measured luminosity depends on the underlying cosmology through the luminosity distance D L , the temperature can be determined without any cosmological assumptions. By exploiting this property and the homogeneous sky coverage of X-ray galaxy cluster samples, one can effectively test the isotropy of cosmological parameters over the full extragalactic sky, which is…

Aims. Understanding the link between the galaxy properties and the dark matter halos they reside in and their coevolution is a powerful tool for constraining the processes related to galaxy formation. In particular, the stellar-to-halo mass relation (SHMR) and its evolution throughout the history of the Universe provides insights on galaxy formation models and allows us to assign galaxy masses to halos in N -body dark matter simulations. To address these questions, we determine the SHMR throughout the entire cosmic history from z ∼ 4 to the present. Methods. We used a statistical approach to link the observed galaxy stellar mass functions on the COSMOS field to dark matter halo mass functions up to z ∼ 4 from the ΛCDM DUSTGRAIN - pathfinder simulation, which is complete for M h > 10 12.5 M ⊙ , and extended this to lower masses with a theoretical parameterization. We propose an empirical model to describe the evolution of the SHMR as a function of redshift (either in the presence or absence of a scatter in stellar mass at fixed halo mass), and compare the results with several…

Context. Molecular gas is a necessary fuel for star formation. The CO (1−0) transition is often used to deduce the total molecular hydrogen but is challenging to detect in low-metallicity galaxies in spite of the star formation taking place. In contrast, the [C II ] λ 158 μ m is relatively bright, highlighting a potentially important reservoir of H 2 that is not traced by CO (1−0) but is residing in the C + -emitting regions. Aims. Here we aim to explore a method to quantify the total H 2 mass ( M H 2 ) in galaxies and to decipher what parameters control the CO-dark reservoir. Methods. We present Cloudy grids of density, radiation field, and metallicity in terms of observed quantities, such as [O I ], [C I ], CO (1−0), [C II ], L TIR , and the total M H 2 . We provide recipes based on these models to derive total M H 2 mass estimates from observations. We apply the models to the Herschel Dwarf Galaxy Survey, extracting the total M H 2 for each galaxy, and compare this to the H 2 determined from the observed CO (1−0) line. This allows us to quantify the reservoir of H 2 that is…

Context. The High Accuracy Radial velocity Planet Searcher (HARPS) spectrograph has been mounted since 2003 at the ESO 3.6 m telescope in La Silla and provides state-of-the-art stellar radial velocity (RV) measurements with a precision down to ∼1 m s −1 . The spectra are extracted with a dedicated data-reduction software (DRS), and the RVs are computed by cross-correlating with a numerical mask. Aims. This study has three main aims: (i) Create easy access to the public HARPS RV data set. (ii) Apply the new public SpEctrum Radial Velocity AnaLyser (SERVAL) pipeline to the spectra, and produce a more precise RV data set. (iii) Determine whether the precision of the RVs can be further improved by correcting for small nightly systematic effects. Methods. For each star observed with HARPS, we downloaded the publicly available spectra from the ESO archive and recomputed the RVs with SERVAL. This was based on fitting each observed spectrum with a high signal-to-noise ratio template created by coadding all the available spectra of that star. We then computed nightly zero-points (NZPs) by…

Here we present a publicly available database of opacities for molecules of astrophysical interest named ExoMolOP that has been compiled for over 80 species, and is based on the latest line list data from the ExoMol, HITEMP, and MoLLIST databases. These data are generally suitable for characterising high-temperature exoplanet or cool stellar and substellar atmospheres, and have been computed at a variety of pressures and temperatures, with a few molecules included at room temperature only from the HITRAN database. The data are formatted in different ways for four different exoplanet atmosphere retrieval codes; ARCiS, TauREx, NEMESIS, and petitRADTRANS, and include both cross sections (at R = λ /Δ λ = 15000) and k -tables (at R = λ /Δ λ = 1000) for the 0.3–50 μ m wavelength region. Opacity files can be downloaded and used directly for these codes. Atomic data for alkali metals Na and K are also included, using data from the NIST database and the latest line shapes for the resonance lines. Broadening parameters have been taken from the literature where available, or have been…

We combine high-resolution spectroscopic data from APOGEE-2 survey Data Release 16 (DR16) with broad-band photometric data from several sources as well as parallaxes from Gaia Data Release 2 (DR2). Using the Bayesian isochrone-fitting code StarHorse, we derived the distances, extinctions, and astrophysical parameters for around 388 815 APOGEE stars. We achieve typical distance uncertainties of ∼6% for APOGEE giants, ∼2% for APOGEE dwarfs, and extinction uncertainties of ∼0.07 mag, when all photometric information is available, and ∼0.17 mag if optical photometry is missing. StarHorse uncertainties vary with the input spectroscopic catalogue, available photometry, and parallax uncertainties. To illustrate the impact of our results, we show that thanks to Gaia DR2 and the now larger sky coverage of APOGEE-2 (including APOGEE-South), we obtain an extended map of the Galactic plane. We thereby provide an unprecedented coverage of the disc close to the Galactic mid-plane (|ZGal| < 1 kpc) from the Galactic centre out to RGal ∼ 20 kpc. The improvements in statistics as well as distance and…

The [C II ] 158 μ m line is one of the strongest IR emission lines, which has been shown to trace the star formation rate (SFR) of galaxies in the nearby Universe, and up to z ∼ 2. Whether this is also the case at higher redshift and in the early Universe remains debated. The ALPINE survey, which targeted 118 star-forming galaxies at 4.4 < z < 5.9, provides a new opportunity to examine this question with the first statistical dataset. Using the ALPINE data and earlier measurements from the literature, we examine the relation between the [C II ] luminosity and the SFR over the entire redshift range from z ∼ 4 − 8. ALPINE galaxies, which are both detected in [C II ] and in dust continuum, show good agreement with the local L ([CII])–SFR relation. Galaxies undetected in the continuum by ALMA are found to be over-luminous in [C II ] when the UV SFR is used. After accounting for dust-obscured star formation, by an amount of SFR(IR) ≈ SFR(UV) on average, which results from two different stacking methods and SED fitting, the ALPINE galaxies show an L ([CII])–SFR relation comparable to the…

We present a new catalogue of ∼2400 optically selected quasars with spectroscopic redshifts and X-ray observations from either Chandra or XMM–Newton . The sample can be used to investigate the non-linear relation between the ultraviolet (UV) and X-ray luminosity of quasars as well as to build a Hubble diagram up to a redshift of z ∼ 7.5. We selected sources that are neither reddened by dust in the optical and UV nor obscured by gas in the X-rays, and whose X-ray fluxes are free from flux-limit-related biases. After checking for any possible systematics, we confirm, in agreement with our previous works, that the X-ray to UV relation provides distance estimates matching those from supernovae up to z ∼ 1.5, and its slope shows no redshift evolution up to z ∼ 5. We provide a full description of the methodology for testing cosmological models, further supporting a trend whereby the Hubble diagram of quasars is well reproduced by the standard flat cold dark matter model up to z ∼ 1.5–2, but strong deviations emerge at higher redshifts. Since we have minimised all non-negligible systematic…

We present a combined tomographic weak gravitational lensing analysis of the Kilo Degree Survey (KV450) and the Dark Energy Survey (DES-Y1). We homogenize the analysis of these two public cosmic shear datasets by adopting consistent priors and modeling of nonlinear scales, and determine new redshift distributions for DES-Y1 based on deep public spectroscopic surveys. Adopting these revised redshifts results in a 0.8σreduction in the DES-inferred value forS​8, which decreases to a 0.5σreduction when including a systematic redshift calibration error model from mock DES data based on the MICE2 simulation. The combined KV450+DES-Y1 constraint onS8= 0.762−0.024+0.025is in tension with thePlanck2018 constraint from the cosmic microwave background at the level of 2.5σ. This result highlights the importance of developing methods to provide accurate redshift calibration for current and future weak-lensing surveys.

Context. WASP-121 b is a hot Jupiter that was recently found to possess rich emission (day side) and transmission (limb) spectra, suggestive of the presence of a multitude of chemical species in the atmosphere. Aims. We survey the transmission spectrum of WASP-121 b for line-absorption by metals and molecules at high spectral resolution and elaborate on existing interpretations of the optical transmission spectrum observed with the Hubble Space Telescope (HST). Methods. We applied the cross-correlation technique and direct differential spectroscopy to search for sodium and other neutral and ionised atoms, TiO, VO, and SH in high-resolution transit spectra obtained with the HARPS spectrograph. We injected models assuming chemical and hydrostatic equilibrium with a varying temperature and composition to enable model comparison, and employed two bootstrap methods to test the robustness of our detections. Results. We detect neutral Mg, Na, Ca, Cr, Fe, Ni, and V, which we predict exists in equilibrium with a significant quantity of VO, supporting earlier observations by HST/WFC3.…

Context. Star formation takes place in giant molecular clouds, resulting in mass-segregated young stellar clusters composed of Sun-like stars, brown dwarfs, and massive O-type(50–100 M ⊙ ) stars. Aims. We aim to identify candidate hub-filament systems (HFSs) in the Milky Way and examine their role in the formation of the highest mass stars and star clusters. Methods. The Herschel survey HiGAL has catalogued about 10 5 clumps. Of these, approximately 35 000 targets are detected at the 3 σ level in a minimum of four bands. Using the DisPerSE algorithm we detect filamentary skeletons on 10′ × 10′ cut-outs of the SPIRE 250 μ m images (18′′ beam width) of the targets. Any filament with a total length of at least 55′′ (3 × 18′′) and at least 18′′ inside the clump was considered to form a junction at the clump. A hub is defined as a junction of three or more filaments. Column density maps were masked by the filament skeletons and averaged for HFS and non-HFS samples to compute the radial profile along the filaments into the clumps. Results. Approximately 3700 (11%) are candidate HFSs, of…

Context. The relevance of M dwarfs in the search for potentially habitable Earth-sized planets has grown significantly in the last years. Aims. In our on-going effort to comprehensively and accurately characterise confirmed and potential planet-hosting M dwarfs, in particular for the CARMENES survey, we have carried out a comprehensive multi-band photometric analysis involving spectral energy distributions, luminosities, absolute magnitudes, colours, and spectral types, from which we have derived basic astrophysical parameters. Methods. We have carefully compiled photometry in 20 passbands from the ultraviolet to the mid-infrared, and combined it with the latest parallactic distances and close-multiplicity information, mostly from Gaia DR2, of a sample of 2479 K5 V to L8 stars and ultracool dwarfs, including 2210 nearby, bright M dwarfs. For this, we made extensive use of Virtual Observatory tools. Results. We have homogeneously computed accurate bolometric luminosities and effective temperatures of 1843 single stars, derived their radii and masses, studied the impact of…

Aims. Mapping the interstellar medium in 3D provides a wealth of insights into its inner working. The Milky Way is the only galaxy for which detailed 3D mapping can be achieved in principle. In this paper, we reconstruct the dust density in and around the local super-bubble. Methods. The combined data from surveys such as Gaia, 2MASS, PANSTARRS, and ALLWISE provide the necessary information to make detailed maps of the interstellar medium in our surrounding. To this end, we used variational inference and Gaussian processes to model the dust extinction density, exploiting its intrinsic correlations. Results. We reconstructed a highly resolved dust map, showing the nearest dust clouds at a distance of up to 400 pc with a resolution of 1 pc. Conclusions. Our reconstruction provides insights into the structure of the interstellar medium. We compute summary statistics of the spectral index and the 1-point function of the logarithmic dust extinction density, which may constrain simulations of the interstellar medium that achieve a similar resolution.

We present weak lensing shear catalogues from the fourth data release of the Kilo-Degree Survey, KiDS-1000, spanning 1006 square degrees of deep and high-resolution imaging. Our ‘gold-sample’ of galaxies, with well-calibrated photometric redshift distributions, consists of 21 million galaxies with an effective number density of 6.17 galaxies per square arcminute. We quantify the accuracy of the spatial, temporal, and flux-dependent point-spread function (PSF) model, verifying that the model meets our requirements to induce less than a 0.1σchange in the inferred cosmic shear constraints on the clustering cosmological parameterS8= σ8√Ωm/0.3.. Through a series of two-point null-tests, we validate the shear estimates, finding no evidence for significant non-lensingB-mode distortions in the data. The PSF residuals are detected in the highest-redshift bins, originating from object selection and/or weight bias. The amplitude is, however, shown to be sufficiently low and within our stringent requirements. With a shear-ratio null-test, we verify the expected redshift scaling of the…

The Radio and Plasma Waves (RPW) instrument on the ESA Solar Orbiter mission is described in this paper. This instrument is designed to measure in-situ magnetic and electric fields and waves from the continuous to a few hundreds of kHz. RPW will also observe solar radio emissions up to 16 MHz. The RPW instrument is of primary importance to the Solar Orbiter mission and science requirements since it is essential to answer three of the four mission overarching science objectives. In addition RPW will exchange on-board data with the other in-situ instruments in order to process algorithms for interplanetary shocks and type III langmuir waves detections.

As part of an on-going effort to identify, understand and correct for astrophysics biases in the standardization of Type Ia supernovae (SN Ia) for cosmology, we have statistically classified a large sample of nearby SNe Ia into those that are located in predominantly younger or older environments. This classification is based on the specific star formation rate measured within a projected distance of 1 kpc from each SN location (LsSFR). This is an important refinement compared to using the local star formation rate directly, as it provides a normalization for relative numbers of available SN progenitors and is more robust against extinction by dust. We find that the SNe Ia in predominantly younger environments are ΔY = 0.163 ± 0.029 mag (5.7σ) fainter than those in predominantly older environments after conventional light-curve standardization. This is the strongest standardized SN Ia brightness systematic connected to the host-galaxy environment measured to date. The well-established step in standardized brightnesses between SNe Ia in hosts with lower or higher total stellar masses…

Context. Complex organic molecules are detected in many sources in the warm inner regions of envelopes surrounding deeply embedded protostars. Exactly how these species form remains an open question. Aims. This study aims to constrain the formation of complex organic molecules through comparisons of their abundances towards the Class 0 protostellar binary IRAS 16293–2422. Methods. We utilised observations from the ALMA Protostellar Interferometric Line Survey of IRAS 16293–2422. The species identification and the rotational temperature and column density estimation were derived by fitting the extracted spectra towards IRAS 16293–2422 A and IRAS 16293–2422 B with synthetic spectra. The majority of the work in this paper pertains to the analysis of IRAS 16293–2422 A for a comparison with the results from the other binary component, which have already been published. Results. We detect 15 different complex species, as well as 16 isotopologues towards the most luminous companion protostar IRAS 16293–2422 A. Tentative detections of an additional 11 isotopologues are reported. We also…

Common-envelope (CE) evolution in massive binary systems is thought to be one of the most promising channels for the formation of compact binary mergers. In the case of merging binary black holes (BBHs), the essential CE phase takes place at a stage when the first BH is already formed and the companion star expands as a supergiant. We aim to decipher the kinds of BH binaries with supergiant companions that could potentially evolve through and survive a CE phase. To this end, we compute envelope binding energies from detailed massive stellar models at different evolutionary stages and metallicities. We make multiple physically extreme choices of assumptions that favor easier CE ejection as well as account for recent advancements in mass-transfer stability criteria. We find that even with the most optimistic assumptions, a successful CE ejection in BH binaries is only possible if the donor is a massive convective-envelope giant, namely a red supergiant (RSG). The same is true for neutron-star binaries with massive companions. In other words, pre-CE progenitors of BBH mergers are BH…

Context. Water is one of the key chemical elements in planetary structure modelling. Due to its complex phase diagram, equations of state often only cover parts of the pressure-temperature space needed in planetary modelling. Aims. We aim to construct an equation of state of H 2 O spanning a very wide range, from 0.1 Pa to 400 TPa and 150 to 10 5 K, which can be used to model the interior of planets. Methods. We combined equations of state valid in localised regions to form a continuous equation of state spanning over the above-mentioned pressure and temperature range. Results. We provide tabulated values for the most important thermodynamic quantities: the density, adiabatic temperature gradient, entropy, internal energy, and bulk speed of sound of water over this pressure and temperature range. For better usability we also calculated density-temperature and density-internal energy grids. We discuss further the impact of this equation of state on the mass radius relation of planets compared to other popular equations of state like ANEOS and QEOS. Conclusions. AQUA is a combination…

Massive binaries that merge as compact objects are the progenitors of gravitational-wave sources. Most of these binaries experience one or more phases of mass transfer, during which one of the stars loses all or part of its outer envelope and becomes a stripped-envelope star. The evolution of the size of these stripped stars is crucial in determining whether they experience further interactions and understanding their ultimate fate. We present new calculations of stripped-envelope stars based on binary evolution models computed with MESA. We use these to investigate their radius evolution as a function of mass and metallicity. We further discuss their pre-supernova observable characteristics and potential consequences of their evolution on the properties of supernovae from stripped stars. At high metallicity, we find that practically all of the hydrogen-rich envelope is removed, which is in agreement with earlier findings. Only progenitors with initial masses below 10 M ⊙ expand to large radii (up to 100 R ⊙ ), while more massive progenitors remain compact. At low metallicity, a…

The era of large transient surveys, gravitational-wave observatories, and multi-messenger astronomy has opened up new possibilities for our understanding of the evolution and final fate of massive stars. Most massive stars are born in binary or higher-order multiple systems and exchange mass with a companion star during their lives. In particular, the progenitors of a large fraction of compact-object mergers, and Galactic neutron stars (NSs) and black holes (BHs) have been stripped of their envelopes by a binary companion. Here, we study the evolution of single and stripped binary stars up to core collapse with the stellar evolution code M ESA and their final fates with a parametric supernova (SN) model. We find that stripped binary stars can have systematically different pre-SN structures compared to genuine single stars and thus also different SN outcomes. These differences are already established by the end of core helium burning and are preserved up to core collapse. Consequently, we find that Case A and B stripped stars and single and Case C stripped stars develop qualitatively…

Context. Ultra-hot Jupiters are excellent laboratories for the study of exoplanetary atmospheres. WASP-121b is one of the most studied; many recent analyses of its atmosphere report interesting features at different wavelength ranges. Aims. In this paper we analyze one transit of WASP-121b acquired with the high-resolution spectrograph ESPRESSO at VLT in one-telescope mode, and one partial transit taken during the commissioning of the instrument in four-telescope mode. Methods. We take advantage of the very high S/N data and of the extreme stability of the spectrograph to investigate the anomalous in-transit radial velocity curve and study the transmission spectrum of the planet. We pay particular attention to the removal of instrumental effects, and stellar and telluric contamination. The transmission spectrum is investigated through single-line absorption and cross-correlation with theoretical model templates. Results. By analyzing the in-transit radial velocities we were able to infer the presence of the atmospheric Rossiter–McLaughlin effect. We measured the height of the…

Context. Spectroscopic surveys of massive galaxy clusters reveal the properties of faint background galaxies thanks to the magnification provided by strong gravitational lensing. Aims. We present a systematic analysis of integral-field-spectroscopy observations of 12 massive clusters, conducted with the Multi Unit Spectroscopic Explorer (MUSE). All data were taken under very good seeing conditions (∼0″​​.6) in effective exposure times between two and 15 h per pointing, for a total of 125 h. Our observations cover a total solid angle of ∼23 arcmin 2 in the direction of clusters, many of which were previously studied by the MAssive Clusters Survey, Frontier Fields (FFs), Grism Lens-Amplified Survey from Space and Cluster Lensing And Supernova survey with Hubble programmes. The achieved emission line detection limit at 5 σ for a point source varies between (0.77–1.5) × 10 −18 erg s −1 cm −2 at 7000 Å. Methods. We present our developed strategy to reduce these observational data, detect continuum sources and line emitters in the datacubes, and determine their redshifts. We constructed…

CONTEXT: Yebes 40m radio telescope is the main and largest observing instrument at Yebes Observatory and it is devoted to Very Long Baseline Interferometry (VLBI) and single dish observations since 2010. It has been covering frequency bands between 2 GHz and 90 GHz in discontinuous and narrow windows in most of the cases, to match the current needs of the European VLBI Network (EVN) and the Global Millimeter VLBI Array (GMVA). AIMS: Nanocosmos project, a European Union funded synergy grant, opened the possibility to increase the instantaneous frequency coverage to observe many molecular transitions with single tunnings in single dish mode. This reduces the observing time and maximises the output from the telescope. METHODS: We present the technical specifications of the recently installed 31.5 - 50GHz (Q band) and 72 - 90.5 GHz (W band) receivers along with the main characteristics of the telescope at these frequency ranges. We have observed IRC+10216, CRL 2688 and CRL 618, which harbour a rich molecular chemistry, to demonstrate the capabilities of the new instrumentation for…

Context. This work is part of the Gaia Data Processing and Analysis Consortium papers published with the Gaia Early Data Release 3 (EDR3). It is one of the demonstration papers aiming to highlight the improvements and quality of the newly published data by applying them to a scientific case. Aims. We use the Gaia EDR3 data to study the structure and kinematics of the Magellanic Clouds. The large distance to the Clouds is a challenge for the Gaia astrometry. The Clouds lie at the very limits of the usability of the Gaia data, which makes the Clouds an excellent case study for evaluating the quality and properties of the Gaia data. Methods. The basis of our work are two samples selected to provide a representation as clean as possible of the stars of the Large Magellanic Cloud (LMC) and the Small Magellanic Cloud (SMC). The selection used criteria based on position, parallax, and proper motions to remove foreground contamination from the Milky Way, and allowed the separation of the stars of both Clouds. From these two samples we defined a series of subsamples based on cuts in the…

Context. Many stars do not live alone, but instead have one or more stellar companions. Observations show that these binaries, triples, and higher-order multiples are common. While the evolution of single stars and binaries have been studied extensively, the same is not true for the evolution of stellar triples. Aims. To fill in this gap in our general understanding of stellar lives, we aim to systematically explore the long-term evolution of triples and to map out the most common evolutionary pathways that triples go through. We quantitatively study how triples evolve, which processes are the most relevant, and how this differs from binary evoluion. Methods. We simulated the evolution of several large populations of triples with a population synthesis approach. We made use of the triple evolution code TRES to simulate the evolution of each triple in a consistent way, including three-body dynamics (based on the secular approach), stellar evolution, and their mutual influences. We simulated the evolution of the system up until mass transfer starts, the system becomes dynamically…

We present cosmological constraints from a joint cosmic shear analysis of the Kilo-Degree Survey (KV450) and the Dark Energy Survey (DES-Y1), which were conducted using Complete Orthogonal Sets of E/B-Integrals (COSEBIs). With COSEBIs, we isolated any B-modes that have a non-cosmic shear origin and demonstrate the robustness of our cosmological E-mode analysis as no significant B-modes were detected. We highlight how COSEBIs are fairly insensitive to the amplitude of the non-linear matter power spectrum at high k-scales, mitigating the uncertain impact of baryon feedback in our analysis. COSEBIs, therefore, allowed us to utilise additional small-scale information, improving the DES-Y1 joint constraints on S8 = σ8(Ωm/0.3)0.5 and Ωm by 20%. By adopting a flat ΛCDM model we find S8 = 0.755−0.021+0.019, which is in 3.2σ tension with the Planck Legacy analysis of the cosmic microwave background.

Context. Complex organic molecules (COMs) have been detected in a few Class 0 protostars but their origin is not well understood. While the usual picture of a hot corino explains their presence as resulting from the heating of the inner envelope by the nascent protostar, shocks in the outflow, disk wind, the presence of a flared disk, or the interaction region between envelope and disk at the centrifugal barrier have also been claimed to enhance the abundance of COMs. Aims. Going beyond studies of individual objects, we want to investigate the origin of COMs in young protostars on a statistical basis. Methods. We use the CALYPSO survey performed with the Plateau de Bure Interferometer of the Institut de Radioastronomie Millimétrique to search for COMs at high angular resolution in a sample of 26 solar-type protostars, including 22 Class 0 and four Class I objects. We derive the column densities of the detected molecules under the local thermodynamic equilibrium approximation and search for correlations between their abundances and with various source properties. Results. Methanol is…

Context. As a result of their formation via massive single and binary stellar evolution, the masses of stellar-remnant black holes (BH) are subjects of great interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star merger events. Aims. In this work, we present new developments in the stellar-remnant formation and related schemes of the current N -body evolution program NBODY7 . We demonstrate that the newly implemented stellar-wind and remnant-formation schemes in the stellar-evolutionary sector or BSE of the NBODY7 code, such as the “rapid” and the “delayed” supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produce neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main-sequence (ZAMS) mass and metallicity, with those from the widely recognised StarTrack population-synthesis program. We also demonstrate the new, recipe-based implementations of various widely debated mechanisms of natal kicks on NSs and BHs, such as…

Aims.The main objective of this project is to characterise chromospheric activity of FGK stars from the HARPS archive. We start, in this first paper, by presenting a catalogue of homogeneously determined chromospheric emission (CE), stellar atmospheric parameters, and ages for 1674 FGK main sequence (MS), subgiant, and giant stars. The analysis of CE level and variability is also performed. Methods.We measured CE in the Ca IIH&K lines using more than 180 000 high-resolution spectra from the HARPS spectrograph, as compiled in the AMBRE project, obtained between 2003 and 2019. We converted the fluxes to bolometric and photospheric corrected chromospheric emission ratio,RHK′. Stellar atmospheric parametersTeff, logg, and [Fe/H] were retrieved from the literature or determined using a homogeneous method.M⋆,R⋆, and ages were determined from isochrone fitting. Results.We show that our sample has a distribution of CE for MS stars that is consistent with an unbiased sample of solar-neighbour MS stars. We analysed the CE distribution for the different luminosity classes and spectral types…

We present dust attenuation properties of spectroscopically confirmed star forming galaxies on the main sequence at a redshift of ∼4.4 − 5.8. Our analyses are based on the far infrared continuum observations of 118 galaxies at rest-frame 158 μ m obtained with the Atacama Large Millimeter Array (ALMA) Large Program to INvestigate [CII] at Early times (ALPINE). We study the connection between the ultraviolet (UV) spectral slope ( β ), stellar mass ( M ⋆ ), and infrared excess (IRX = L IR / L UV ). Twenty-three galaxies are individually detected in the continuum at > 3.5 σ significance. We perform a stacking analysis using both detections and nondetections to study the average dust attenuation properties at z ∼ 4.4 − 5.8. The individual detections and stacks show that the IRX– β relation at z ∼ 5 is consistent with a steeper dust attenuation curve than typically found at lower redshifts ( z < 4). The attenuation curve is similar to or even steeper than that of the extinction curve of the Small Magellanic Cloud. This systematic change of the IRX– β relation as a function of redshift…

Context.Methods used to detect giant exoplanets can be broadly divided into two categories: indirect and direct. Indirect methods are more sensitive to planets with a small orbital period, whereas direct detection is more sensitive to planets orbiting at a large distance from their host star. This dichotomy makes it difficult to combine the two techniques on a single target at once. Aims.Simultaneous measurements made by direct and indirect techniques offer the possibility of determining the mass and luminosity of planets and a method of testing formation models. Here, we aim to show how long-baseline interferometric observations guided by radial-velocity can be used in such a way. Methods.We observed the recently-discovered giant planetβPictoris c with GRAVITY, mounted on the Very Large Telescope Interferometer. Results.This study constitutes the first direct confirmation of a planet discovered through radial velocity. We find that the planet has a temperature ofT = 1250 ± 50 K and a dynamical mass ofM = 8.2 ± 0.8 MJup. At 18.5 ± 2.5 Myr, this putsβPic c close to a ‘hot start’…

Context. The changing-look phenomenon observed in a growing number of active galaxies challenges our understanding of the accretion process close to a black hole. Aims. We propose a simple explanation for the sources where multiple semi-periodic outbursts are observed, and the sources are operating close to the Eddington limit. Methods. The outburst are caused by the radiation pressure instability operating in the narrow ring between the standard gas-dominated outer disk and the hot optically thin inner advection-dominated accretion flow. The corresponding limit cycle is responsible for periodic outbursts, and the timescales are much shorter than the standard viscous timescale due to the narrowness of the unstable radial zone. Results. Our toy model gives quantitative predictions and works well for multiple outbursts like those observed in NGC 1566, NGC 4151, NGC 5548, and GSN 069, although the shapes of the outbursts are not yet well modeled, and further development of the model is necessary.

Context.The intriguing binary LS V +22 25 (LB-1) has drawn much attention following claims of it being a single-lined spectroscopic binary with a 79-day orbit comprising a B-type star and a ≈70 M⊙black hole – the most massive stellar black hole reported to date. Subsequent studies demonstrated a lack of evidence for a companion of such great mass. Recent analyses have implied that the primary star is a stripped He-rich star with peculiar sub-solar abundances of heavy elements, such as Mg and Fe. However, the nature of the secondary, which was proposed to be a black hole, a neutron star, or a main sequence star, remains unknown. Aims.Based on 26 newly acquired spectroscopic observations secured with the HERMES and FEROS spectrographs covering the orbit of the system, we perform an orbital analysis and spectral disentangling of LB-1 to elucidate the nature of the system. Methods.To derive the radial velocity semi-amplitudeK2of the secondary and extract the spectra of the two components, we used two independent disentangling methods: the shift-and-add technique and Fourier…

The structure of magnetic flux ropes injected into the solar wind during reconnection in the coronal atmosphere is explored with particle-in-cell simulations and compared with in situ measurements of magnetic “switchbacks” from the Parker Solar Probe. We suggest that multi-x-line reconnection between open and closed flux in the corona injects flux ropes into the solar wind and that these flux ropes convect outward over long distances before eroding due to reconnection. Simulations that explore the magnetic structure of flux ropes in the solar wind reproduce the following key features of the switchback observations: a rapid rotation of the radial magnetic field into the transverse direction, which is a consequence of reconnection with a strong guide field; and the potential to reverse the radial field component. The potential implication of the injection of large numbers of flux ropes in the coronal atmosphere for understanding the generation of the solar wind is discussed.

Context. Directly imaged planets and substellar companions are key targets for the characterization of self-luminous atmospheres. Their photometric appearance at 4–5 μ m is sensitive to the chemical composition and cloud content of their atmosphere. Aims. We aim to systematically characterize the atmospheres of directly imaged low-mass companions at 4–5 μ m. We want to homogeneously process the data, provide robust flux measurements, and compile a photometric library at thermal wavelengths of these mostly young, low-gravity objects. In this way, we want to find trends related to their spectral type and surface gravity by comparing with isolated brown dwarfs and predictions from atmospheric models. Methods. We used the high-resolution, high-contrast capabilities of NACO at the Very Large Telescope (VLT) to directly image the companions of HIP 65426, PZ Tel, and HD 206893 in the NB4.05 and/or M ′ filters. For the same targets, and additionally β Pic, we also analyzed six archival VLT/NACO datasets which were taken with the NB3.74, L ′, NB4.05, and M ′ filters. The data processing and…

We present the methodology for a joint cosmological analysis of weak gravitational lensing from the fourth data release of the ESO Kilo-Degree Survey (KiDS-1000) and galaxy clustering from the partially overlapping Baryon Oscillation Spectroscopic Survey (BOSS) and the 2-degree Field Lensing Survey (2dFLenS). Cross-correlations between BOSS and 2dFLenS galaxy positions and source galaxy ellipticities have been incorporated into the analysis, necessitating the development of a hybrid model of non-linear scales that blends perturbative and non-perturbative approaches, and an assessment of signal contributions by astrophysical effects. All weak lensing signals were measured consistently via Fourier-space statistics that are insensitive to the survey mask and display low levels of mode mixing. The calibration of photometric redshift distributions and multiplicative gravitational shear bias has been updated, and a more complete tally of residual calibration uncertainties was propagated into the likelihood. A dedicated suite of more than 20 000 mocks was used to assess the performance of…

Context. The census of open clusters in the Milky Way is in a never-before seen state of flux. Recent works have reported hundreds of new open clusters thanks to the incredible astrometric quality of the Gaia satellite, but other works have also reported that many open clusters discovered in the pre Gaia era may be associations. Aims. We aim to conduct a comparison of clustering algorithms used to detect open clusters, attempting to statistically quantify their strengths and weaknesses by deriving the sensitivity, specificity, and precision of each as well as their true positive rate against a larger sample. Methods. We selected DBSCAN, HDBSCAN, and Gaussian mixture models for further study, owing to their speed and appropriateness for use with Gaia data. We developed a preprocessing pipeline for Gaia data and developed the algorithms further for the specific application to open clusters. We derived detection rates for all 1385 open clusters in the fields in our study as well as more detailed performance statistics for 100 of these open clusters. Results. DBSCAN was sensitive to…

Aims. Previous high-precision studies of abundances of elements in solar twin stars are extended to a wider metallicity range to see how the trends of element ratios with stellar age depend on [Fe/H]. Methods. HARPS spectra with signal-to-noise ratios S/N ≳ 600 at λ​ ∼ ​6000 Å were analysed with MARCS model atmospheres to obtain 1D LTE abundances of C, O, Na, Mg, Al, Si, Ca, Ti, Cr, Fe, Ni, Sr, and Y for 72 nearby solar-type stars with metallicities in the range of −0.3 ≲ [Fe/H] ≲ +0.3 and ASTEC stellar models were used to determine stellar ages from effective temperatures, luminosities obtained via Gaia DR2 parallaxes, and heavy element abundances. Results. The age-metallicity distribution appears to consist of the following two distinct populations: a sequence of old stars with a steep rise of [Fe/H] to ∼​ + 0.3 dex at an age of ∼​7 Gyr and a younger sequence with [Fe/H] increasing from about −0.3 dex to ∼​ + 0.2 dex over the last 6 Gyr. Furthermore, the trends of several abundance ratios, [O/Fe], [Na/Fe], [Ca/Fe], and [Ni/Fe], as a function of stellar age, split into two…

Context. Gravitational microlensing is sensitive to compact-object lenses in the Milky Way, including white dwarfs, neutron stars, or black holes, and could potentially probe a wide range of stellar-remnant masses. However, the mass of the lens can be determined only in very limited cases, due to missing information on both source and lens distances and their proper motions. Aims. Our aim is to improve the mass estimates in the annual parallax microlensing events found in the eight years of OGLE-III observations towards the Galactic Bulge with the use of Gaia Data Release 2 (DR2). Methods. We use Gaia DR2 data on distances and proper motions of non-blended sources and recompute the masses of lenses in parallax events. We also identify new events in that sample which are likely to have dark lenses; the total number of such events is now 18. Results. The derived distribution of masses of dark lenses is consistent with a continuous distribution of stellar-remnant masses. A mass gap between neutron star and black hole masses in the range between 2 and 5 solar masses is not favoured by…

Context. Radio relics in galaxy clusters are giant diffuse synchrotron sources powered in cluster outskirts by merger shocks. Although the relic–shock connection has been consolidated in recent years by a number of observations, the details of the mechanisms leading to the formation of relativistic particles in this environment are still not well understood. Aims. The diffusive shock acceleration (DSA) theory is a commonly adopted scenario to explain the origin of cosmic rays at astrophysical shocks, including those in radio relics in galaxy clusters. However, in a few specific cases it has been shown that the energy dissipated by cluster shocks is not enough to reproduce the luminosity of the relics via DSA of thermal particles. Studies based on samples of radio relics are required to further address this limitation of the mechanism. Methods. In this paper, we focus on ten well-studied radio relics with underlying shocks observed in the X-rays and calculate the electron acceleration efficiency of these shocks that is necessary to reproduce the observed radio luminosity of the…

Aims. The dynamics of coalescing compact binaries can be affected by the environment in which the systems evolve, leaving detectable signatures in the emitted gravitational signal. In this paper, we investigate the ability of gravitational-wave detectors to constrain the nature of the environment in which compact binaries merge. Methods. We parametrized a variety of environmental effects by modifying the phase of the gravitational signal emitted by black hole and neutron star binaries. We infer the bounds on such effects by current and future generations of interferometers, studying their dependence on the binary’s parameters. Results. We show that the strong dephasing induced by accretion and dynamical friction can constrain the density of the surrounding medium to orders of magnitude below those of accretion disks. Planned detectors, such as LISA or DECIGO, will be able to probe densities typical of those of dark matter.

Aims. The SINFONI survey for Unveiling the Physics and Effect of Radiative feedback (SUPER) aims to trace and characterise ionised gas outflows and their impact on star formation in a statistical sample of X-ray selected active galactic nuclei (AGN) at z ∼ 2. We present the first SINFONI results for a sample of 21 Type 1 AGN spanning a wide range in bolometric luminosity (log L bol = 45.4–47.9 erg s −1 ). The main aims of this paper are to determine the extension of the ionised gas, characterise the occurrence of AGN-driven outflows, and link the properties of such outflows with those of the AGN. Methods. We used adaptive optics-assisted SINFONI observations to trace ionised gas in the extended narrow line region using the [O III ] λ 5007 line. We classified a target as hosting an outflow if its non-parametric velocity of the [O III ] line, w 80 , was larger than 600 km s −1 . We studied the presence of extended emission using dedicated point-spread function (PSF) observations, after modelling the PSF from the Balmer lines originating from the broad line region. Results. We detect…

Context.TheβPictoris moving group is one of the most well-known young associations in the solar neighbourhood and several members are known to host circumstellar discs, planets, and comets. Measuring its age precisely is essential to the study of several astrophysical processes, such as planet formation and disc evolution, which are strongly age-dependent. Aims.We aim to determine a precise and accurate dynamical traceback age for theβPictoris moving group. Methods.Our sample combines the extremely preciseGaiaDR2 astrometry with ground-based radial velocities measured in an homogeneous manner. We use an updated version of our algorithm to determine dynamical ages. The new approach takes into account a robust estimate of the spatial and kinematic covariance matrices of the association to improve the sample selection process and to perform the traceback analysis. Results.We estimate a dynamical age of 18.5−2.4+2.0Myr for theβPictoris moving group. We investigated the spatial substructure of the association at the time of birth and we propose the existence of a core of stars that is…

Solar-like stars (M ≲ 1.3 M⊙) lose angular momentum through their magnetized winds. The resulting evolution of the surface rotation period, which can be directly measured photometrically, has the potential to be an accurate indicator of stellar age, and is constrained by observations of rotation periods of coeval stars, such as members of Galactic open clusters. A prominent observational feature of the mass–rotation period diagrams of open clusters is a sequence of relatively slower rotators. The formation and persistence of this slow-rotator sequence across several billion years imply an approximately coherent spin-down of the stars that belong to it. In particular, the sequence is observed to evolve coherently toward longer periods in progressively older clusters. Recent observations of the ≈700 Myr Praesepe and the 1 Gyr NGC 6811 clusters, however, are not fully consistent with this general pattern. While the stars of 1 M⊙ on the slow-rotator sequence of the older NGC 6811 have longer periods than their counterparts in the younger Praesepe, as expected, the two sequences…

Accurate photometric redshift calibration is central to the robustness of all cosmology constraints from cosmic shear surveys. Analyses of the Kilo-Degree Survey (KiDS) re-weighted training samples from all overlapping spectroscopic surveys to provide a direct redshift calibration. Using self-organising maps we demonstrate that this spectroscopic compilation is sufficiently complete for KiDS, representing 99% of the effective 2D cosmic shear sample. We used the SOM to define a 100% represented “gold” cosmic shear sample, per tomographic bin. Using mock simulations of KiDS and the spectroscopic training set, we estimated the uncertainty on the SOM redshift calibration, and we find that photometric noise, sample variance, and spectroscopic selection effects (including redshift and magnitude incompleteness) induce a combined maximal scatter on the bias of the redshift distribution reconstruction (Δ⟨z⟩ = ⟨z⟩est − ⟨z⟩true) of σΔ⟨z⟩ ≤ 0.006 in all tomographic bins. Photometric noise and spectroscopic selection effects contribute equally to the observed scatter. We show that the SOM…

Context.Until recently, camera networks designed for monitoring fireballs worldwide were not fully automated, implying that in case of a meteorite fall, the recovery campaign was rarely immediate. This was an important limiting factor as the most fragile – hence precious – meteorites must be recovered rapidly to avoid their alteration. Aims.The Fireball Recovery and InterPlanetary Observation Network (FRIPON) scientific project was designed to overcome this limitation. This network comprises a fully automated camera and radio network deployed over a significant fraction of western Europe and a small fraction of Canada. As of today, it consists of 150 cameras and 25 European radio receivers and covers an area of about 1.5 × 106km2. Methods.The FRIPON network, fully operational since 2018, has been monitoring meteoroid entries since 2016, thereby allowing the characterization of their dynamical and physical properties. In addition, the level of automation of the network makes it possible to trigger a meteorite recovery campaign only a few hours after it reaches the surface of the…

The LOFAR Two-metre Sky Survey (LoTSS) will cover the full northern sky and, additionally, aims to observe the LoTSS deep fields to a noise level of ≲10μJy beam−1over several tens of square degrees in areas that have the most extensive ancillary data. This paper presents the ELAIS-N1 deep field, the deepest of the LoTSS deep fields to date. With an effective observing time of 163.7 h, it reaches a root mean square noise level of ≲20μJy beam−1in the central region (and below 30μJy beam−1over 10 square degrees). The resolution is ~6 arcsecs and 84 862 radio sources were detected in the full area (68 square degrees) with 74 127 sources in the highest quality area at less than 3 degrees from the pointing centre. The observation reaches a sky density of more than 5000 sources per square degree in the central region (~5 square degrees). We present the calibration procedure, which addresses the special configuration of some observations and the extended bandwidth covered (115–177 MHz; central frequency 146.2 MHz) compared to standard LoTSS. We also describe the methods used to calibrate…

We present a novel view on the morphology and dynamical state of ten prominent, nearby (≤500 pc), and young (∼30–300 Myr) open star clusters withGaiaDR2:α Per, Blanco 1, IC 2602, IC 2391, Messier 39, NGC 2451A, NGC 2516, NGC 2547, Platais 9, and the Pleiades. We introduce a pioneering member-identification method that is informed by cluster bulk velocities and deconvolves the spatial distribution with a mixture of Gaussians. Our approach enables inferring the true spatial distribution of the clusters by effectively filtering field star contaminants while at the same time mitigating the effect of positional errors along the line of sight. This first application of the method reveals vast stellar coronae that extend for ≳100 pc and surround the cluster cores, which are comparatively tiny and compact. The coronae and cores form intertwined, coeval, and comoving extended cluster populations, each encompassing tens of thousands of cubic parsec and stretching across tens of degrees on the sky. Our analysis shows that the coronae are gravitationally unbound but largely comprise the bulk of…

The molecular gas content of normal galaxies at z > 4 is poorly constrained because the commonly used molecular gas tracers become hard to detect at these high redshifts. We use the [C II ] 158 μ m luminosity, which was recently proposed as a molecular gas tracer, to estimate the molecular gas content in a large sample of main sequence star-forming galaxies at z = 4.4 − 5.9, with a median stellar mass of 10 9.7 M ⊙ , drawn from the ALMA Large Program to INvestigate [C II ] at Early times survey. The agreement between the molecular gas masses derived from [C II ] luminosities, dynamical masses, and rest-frame 850 μ m luminosities extrapolated from the rest-frame 158 μ m continuum supports [C II ] as a reliable tracer of molecular gas in our sample. We find a continuous decline of the molecular gas depletion timescale from z = 0 to z = 5.9, which reaches a mean value of (4.6 ± 0.8) × 10 8 yr at z ∼ 5.5, only a factor of between two and three shorter than in present-day galaxies. This suggests a mild enhancement of the star formation efficiency toward high redshifts. Our estimates also…

Ultra-hot Jupiters offer interesting prospects for expanding our theories on dynamical evolution and the properties of extremely irradiated atmospheres. In this context, we present the analysis of new optical spectroscopy for the transiting ultra-hot Jupiter WASP-121b. We first refine the orbital properties of WASP-121b, which is on a nearly polar (obliquity ψ North = 88.1 ± 0.25° or ψ South = 91.11 ± 0.20°) orbit, and exclude a high differential rotation for its fast-rotating ( P < 1.13 days), highly inclined ( i ⋆ North = 8.1 −2.6 +3.0 ° or i ⋆ South = 171.9 −3.4 +2.5 °) star. We then present a new method that exploits the reloaded Rossiter-McLaughlin technique to separate the contribution of the planetary atmosphere and of the spectrum of the stellar surface along the transit chord. Its application to HARPS transit spectroscopy of WASP-121b reveals the absorption signature from metals, likely atomic iron, in the planet atmospheric limb. The width of the signal (14.3 ± 1.2 km s −1 ) can be explained by the rotation of the tidally locked planet. Its blueshift (−5.2 ± 0.5 km s −1 )…

The ubiquitous presence of Galactic cirri in deep optical images represents a major obstacle to study the low surface brightness features of extragalactic sources. To address this issue, we have explored the optical properties of cirri usingg,r,i, andzbands in the Sloan Digital Sky Survey (SDSS) Stripe82 region. Using state-of-the-art, custom made, image processing techniques, including the modeling and removal of the instrumental scattered light produced by the stars, we managed to isolate the optical diffuse emission by the cirri, allowing their photometric characterization. We find that their optical colors are driven by the dust column density: The cirri become redder as their 100μm emission increases. Remarkably, the optical colors of the Galactic cirri differ significantly from those of extragalactic sources, with a characteristic bluerr-icolor for a giveng-r, allowing one to detect these by using a simple color relation. Our results show the high potential of deep multi-band optical photometry, on its own, identifying the presence of cirri at a higher spatial resolution than…