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

In the original version, the bounds given in Eqs. (87a) and (87b) on the contribution to the early-time optical depth, (15,30), contained a numerical error in deriving the 95th percentile from the Monte Carlo samples. The corrected 95% upper bounds are: τ(15,30) < 0:018 (lowE, flat τ(15, 30), FlexKnot), (1) τ(15, 30) < 0:023 (lowE, flat knot, FlexKnot): (2) These bounds are a factor of 3 larger than the originally reported results. Consequently, the new bounds do not significantly improve upon previous results from Planck data presented in Millea & Bouchet (2018) as was stated, but are instead comparable. Equations (1) and (2) give results that are now similar to those of Heinrich & Hu (2021), who used the same Planck 2018 data to derive a 95% upper bound of 0.020 using the principal component analysis (PCA) model and uniform priors on the PCA mode amplitudes.

Context. The chemical composition of the Sun is a fundamental yardstick in astronomy, relative to which essentially all cosmic objects are referenced. As such, having accurate knowledge of the solar elemental abundances is crucial for an extremely broad range of topics. Aims. We reassess the solar abundances of all 83 long-lived elements, using highly realistic solar modelling and state-of-the-art spectroscopic analysis techniques coupled with the best available atomic data and observations. Methods. The basis for our solar spectroscopic analysis is a three-dimensional (3D) radiative-hydrodynamical model of the solar surface convection and atmosphere, which reproduces the full arsenal of key observational diagnostics. New complete and comprehensive 3D spectral line formation calculations taking into account of departures from local thermodynamic equilibrium (non-LTE) are presented for Na, Mg, K, Ca, and Fe using comprehensive model atoms with reliable radiative and collisional data. Our newly derived abundances for C, N, and O are based on a 3D non-LTE analysis of permitted and…

The orbital observatory Spectrum-Roentgen-Gamma (SRG), equipped with the grazing-incidence X-ray telescopes Mikhail Pavlinsky ART-XC and eROSITA, was launched by Roscosmos to the Lagrange L2 point of the Sun–Earth system on July 13, 2019. The launch was carried out from the Baikonur Cosmodrome by a Proton-M rocket with a DM-03 upper stage. The German telescope eROSITA was installed on SRG under an agreement between Roskosmos and the DLR, the German Aerospace Agency. In December 2019, SRG started to perform its main scientific task: scanning the celestial sphere to obtain X-ray maps of the entire sky in several energy ranges (from 0.2 to 8 keV with eROSITA, and from 4 to 30 keV with ART-XC). By mid-June 2021, the third six-month all-sky survey had been completed. Over a period of four years, it is planned to obtain eight independent maps of the entire sky in each of the energy ranges. The sum of these maps will provide high sensitivity and reveal more than three million quasars and over one hundred thousand massive galaxy clusters and galaxy groups. The availability of eight sky maps…

We report the detection for the first time in space of three new pure hydrocarbon cycles in TMC-1: c-C3HCCH (ethynyl cyclopropenylidene), c-C5H6 (cyclopentadiene), and c-C9H8 (indene). We derive a column density of 3.1 × 1011 cm−2 for the first cycle and similar values, in the range (1−2) × 1013 cm−2, for the second and third. This means that cyclopentadiene and indene, in spite of their large size, are exceptionally abundant, only a factor of five less abundant than the ubiquitous cyclic hydrocarbon c-C3H2. The high abundance found for these two hydrocarbon cycles together with the high abundance previously found for the propargyl radical (CH2CCH) and other hydrocarbons, such as vinyl and allenyl acetylene (Agúndez et al. 2021, A&A, 647, L10; Cernicharo et al. 2021a, A&A, 647, L2; Cernicharo et al. 2021b, A&A, 647, L3), start to allow us to quantify the abundant content of hydrocarbon rings in cold dark clouds and to identify the intermediate species that are probably behind the in situ bottom-up synthesis of aromatic cycles in these environments. While c-C3HCCH is most likely…

Stars orbiting the compact radio source Sgr A* in the Galactic Center serve as precision probes of the gravitational field around the closest massive black hole. In addition to adaptive optics-assisted astrometry (with NACO/VLT) and spectroscopy (with SINFONI/VLT, NIRC2/Keck and GNIRS/Gemini) over three decades, we have obtained 30–100 μas astrometry since 2017 with the four-telescope interferometric beam combiner GRAVITY/VLTI, capable of reaching a sensitivity of mK = 20 when combining data from one night. We present the simultaneous detection of several stars within the diffraction limit of a single telescope, illustrating the power of interferometry in the field. The new data for the stars S2, S29, S38, and S55 yield significant accelerations between March and July 2021, as these stars pass the pericenters of their orbits between 2018 and 2023. This allows for a high-precision determination of the gravitational potential around Sgr A*. Our data are in excellent agreement with general relativity orbits around a single central point mass, M• = 4.30 × 106 M⊙, with a precision of…

Context. The multiplicity fraction of stars, down to the substellar regime, is a parameter of fundamental importance for stellar formation, evolution, and planetology. The census of multiple stars in the solar neighborhood is however incomplete. Aims. Our study is aimed at detecting companions of H IPPARCOS catalog stars from the proper motion anomaly (PMa) they induce on their host star, namely, the difference between their long-term H IPPARCOS - Gaia and short-term Gaia proper motion vectors. We also aim to detect resolved, gravitationally bound companions of the H IPPARCOS catalog stars (117 955 stars) and of the Gaia EDR3 stars closer than 100 pc (542 232 stars). Methods. Using the H IPPARCOS and EDR3 data, we revised the PMa catalog for the H IPPARCOS stars. In order to identify gravitationally bound visual companions of our sample, we searched the Gaia EDR3 catalog for common proper-motion (CPM) candidates. Results. The detection of tangential velocity anomalies with a median accuracy of σ (Δ v T ) = 26 cm s −1 per parsec of distance is demonstrated with the EDR3. This…

We present a catalogue of 362 million stellar parameters, distances, and extinctions derived from Gaia ’s Early Data Release (EDR3) cross-matched with the photometric catalogues of Pan-STARRS1, SkyMapper, 2MASS, and AllWISE. The higher precision of the Gaia EDR3 data, combined with the broad wavelength coverage of the additional photometric surveys and the new stellar-density priors of the StarHorse code, allows us to substantially improve the accuracy and precision over previous photo-astrometric stellar-parameter estimates. At magnitude G = 14 (17), our typical precisions amount to 3% (15%) in distance, 0.13 mag (0.15 mag) in V -band extinction, and 140 K (180 K) in effective temperature. Our results are validated by comparisons with open clusters, as well as with asteroseismic and spectroscopic measurements, indicating systematic errors smaller than the nominal uncertainties for the vast majority of objects. We also provide distance- and extinction-corrected colour-magnitude diagrams, extinction maps, and extensive stellar density maps that reveal detailed substructures in the…

Hydrodynamical cosmological simulations are increasing their level of realism by considering more physical processes and having greater resolution or larger statistics. However, usually either the statistical power of such simulations or the resolution reached within galaxies are sacrificed. Here, we introduce the N EW H ORIZON project in which we simulate at high resolution a zoom-in region of ∼(16 Mpc) 3 that is larger than a standard zoom-in region around a single halo and is embedded in a larger box. A resolution of up to 34 pc, which is typical of individual zoom-in, up-to-date resimulated halos, is reached within galaxies; this allows the simulation to capture the multi-phase nature of the interstellar medium and the clumpy nature of the star formation process in galaxies. In this introductory paper, we present several key fundamental properties of galaxies and their black holes, including the galaxy mass function, cosmic star formation rate, galactic metallicities, the Kennicutt–Schmidt relation, the stellar-to-halo mass relation, galaxy sizes, stellar kinematics and…

ERRATUM This article is an erratum for: [https://doi.org/10.1051/0004-6361/202039657] ​​​​​​​

Context.Previous theoretical works on planet formation around low-mass stars have often been limited to large planets and individual systems. As current surveys routinely detect planets down to terrestrial size in these systems, models have shifted toward a more holistic approach that reflects their diverse architectures. Aims.Here, we investigate planet formation around low-mass stars and identify differences in the statistical distribution of modeled planets. We compare the synthetic planet populations to observed exoplanets and we discuss the identified trends. Methods.We used the Generation III Bern global model of planet formation and evolution to calculate synthetic populations, while varying the central star from Solar-like stars to ultra-late M dwarfs. This model includes planetary migration,N-body interactions between embryos, accretion of planetesimals and gas, and the long-term contraction and loss of the gaseous atmospheres. Results.We find that temperate, Earth-sized planets are most frequent around early M dwarfs (0.3M⊙–0.5M⊙) and that they are more rare for Solar-type…

Context. The explosion of observational data on exoplanets gives many constraints on theoretical models of planet formation and evolution. Observational data probe very large areas of the parameter space and many different planet properties. Aims. Comparing theoretical models with observations allows one to take a key step forward towards understanding planetary systems. It however requires a model able to (i) predict all the necessary observable quantities (not only masses and orbits, but also radii, luminosities, magnitudes, or evaporation rates) and (ii) address the large range in relevant planetary masses (from Mars mass to super-Jupiters) and distances (from stellar-grazing to wide orbits). Methods. We have developed a combined global end-to-end planetary formation and evolution model, the Generation III Bern model, based on the core accretion paradigm. This model solves as directly as possible the underlying differential equations for the structure and evolution of the gas disc, the dynamical state of the planetesimals, the internal structure of the planets yielding their…

The GRAVITY instrument on the ESO VLTI pioneers the field of high-precision near-infrared interferometry by providing astrometry at the 10−100 μas level. Measurements at this high precision crucially depend on the control of systematic effects. We investigate how aberrations introduced by small optical imperfections along the path from the telescope to the detector affect the astrometry. We develop an analytical model that describes the effect of these aberrations on the measurement of complex visibilities. Our formalism accounts for pupil-plane and focal-plane aberrations, as well as for the interplay between static and turbulent aberrations, and it successfully reproduces calibration measurements of a binary star. The Galactic Center observations with GRAVITY in 2017 and 2018, when both Sgr A* and the star S2 were targeted in a single fiber pointing, are affected by these aberrations at a level lower than 0.5 mas. Removal of these effects brings the measurement in harmony with the dual-beam observations of 2019 and 2020, which are not affected by these aberrations. This also…

Determining the architecture of multi-planetary systems is one of the cornerstones of understanding planet formation and evolution. Resonant systems are especially important as the fragility of their orbital configuration ensures that no significant scattering or collisional event has taken place since the earliest formation phase when the parent protoplanetary disc was still present. In this context, TOI-178 has been the subject of particular attention since the first TESS observations hinted at the possible presence of a near 2:3:3 resonant chain. Here we report the results of observations from CHEOPS, ESPRESSO, NGTS, and SPECULOOS with the aim of deciphering the peculiar orbital architecture of the system. We show that TOI-178 harbours at least six planets in the super-Earth to mini-Neptune regimes, with radii ranging from 1.152−0.070+0.073 to 2.87−0.13+0.14 Earth radii and periods of 1.91, 3.24, 6.56, 9.96, 15.23, and 20.71 days. All planets but the innermost one form a 2:4:6:9:12 chain of Laplace resonances, and the planetary densities show important variations from planet to…

As the number of merging binary black holes observed with ground-based gravitational-wave detectors grows, increasingly accurate theoretical models are required to compare them to the observed sample and disentangle contributions from multiple channels. In formation models involving isolated binary stars, important uncertainties remain regarding the stability of mass transfer and common-envelope evolution. To study some of these uncertainties, we have computed binary simulations using the MESA code consisting of a 30 M ⊙ star in a low metallicity ( Z ⊙ /10) environment with a black-hole companion. We have developed an updated prescription to compute mass transfer rates including the possibility of outflows from outer Lagrangian points, as well as a method to self-consistently determine the core-envelope boundary in cases where there is common-envelope evolution. We find that binaries survive common-envelope evolution only if unstable mass transfer happens after the formation of a deep convective envelope, resulting in a narrow range (0.2 dex) in period for successful envelope…

We report the detection, for the first time in space, of a new non-functionalised hydrocarbon cycle in the direction of TMC-1: o -C 6 H 4 (ortho-benzyne). We derive a column density for this hydrocarbon cycle of (5.0 ± 1.0) × 10 11 cm −2 . The abundance of this species is around 30 times lower than that of cyclopentadiene and indene. We compare the abundance of benzyne with that of other pure hydrocarbons, cycles or chains, and find that it could be formed from neutral-radical reactions such as C 2 H + CH 2 CHCCH and C + C 5 H 5 , and possibly through C 4 H + C 2 H 4 , C 3 H + CH 2 CCH 2 , and C 3 H 2 + C 3 H 3 . Hence, the rich content of hydrocarbon cycles observed in TMC-1 could arise through a bottom-up scenario involving reactions of a few radicals with the abundant hydrocarbons recently revealed by the QUIJOTE line survey.

At least 30% of main sequence stars host planets with sizes of between 1 and 4 Earth radii and orbital periods of less than 100 days. We useN-body simulations including a model for gas-assisted pebble accretion and disk–planet tidal interaction to study the formation of super-Earth systems. We show that the integrated pebble mass reservoir creates a bifurcation between hot super-Earths or hot-Neptunes (≲15M⊕) and super-massive planetary cores potentially able to become gas giant planets (≳15M⊕). Simulations with moderate pebble fluxes grow multiple super-Earth-mass planets that migrate inwards and pile up at the inner edge of the disk forming long resonant chains. We follow the long-term dynamical evolution of these systems and use the period ratio distribution of observed planet-pairs to constrain our model. Up to ~95% of resonant chains become dynamically unstable after the gas disk dispersal, leading to a phase of late collisions that breaks the original resonant configurations. Our simulations naturally match observations when they produce a dominant fraction (≳95%) of unstable…

Context. The heating of the solar corona by small heating events requires an increasing number of such events at progressively smaller scales, with the bulk of the heating occurring at scales that are currently unresolved. Aims. The goal of this work is to study the smallest brightening events observed in the extreme-UV quiet Sun. Methods. We used commissioning data taken by the Extreme Ultraviolet Imager (EUI) on board the recently launched Solar Orbiter mission. On 30 May 2020, the EUI was situated at 0.556 AU from the Sun. Its High Resolution EUV telescope (HRI EUV , 17.4 nm passband) reached an exceptionally high two-pixel spatial resolution of 400 km. The size and duration of small-scale structures was determined by the HRI EUV data, while their height was estimated from triangulation with simultaneous images from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory mission. This is the first stereoscopy of small-scale brightenings at high resolution. Results. We observed small localised brightenings, also known as ‘campfires’, in a quiet Sun region…

Context. Water is a key molecule in the physics and chemistry of star and planet formation, but it is difficult to observe from Earth. The Herschel Space Observatory provided unprecedented sensitivity as well as spatial and spectral resolution to study water. The Water In Star-forming regions with Herschel (WISH) key program was designed to observe water in a wide range of environments and provide a legacy data set to address its physics and chemistry. Aims. The aim of WISH is to determine which physical components are traced by the gas-phase water lines observed with Herschel and to quantify the excitation conditions and water abundances in each of these components. This then provides insight into how and where the bulk of the water is formed in space and how it is transported from clouds to disks, and ultimately comets and planets. Methods. Data and results from WISH are summarized together with those from related open time programs. WISH targeted ~80 sources along the two axes of luminosity and evolutionary stage: from low- to high-mass protostars (luminosities from <1 to > 10 5…

Recent observations of extrasolar gas giants suggest super-stellar C/O ratios in planetary atmospheres, while interior models of observed extrasolar giant planets additionally suggest high heavy element contents. Furthermore, recent observations of protoplanetary disks revealed super-solar C/H ratios, which are explained by inward drifting and evaporating pebbles enhancing the volatile content of the disk. We investigate in this work how the inward drift and evaporation of volatile-rich pebbles influences the atmospheric C/O ratio and heavy element content of giant planets growing by pebble and gas accretion. To achieve this goal, we perform semi-analytical 1D models of protoplanetary disks, including the treatment of viscous evolution and heating, pebble drift, and simple chemistry to simulate the growth of planets from planetary embryos to Jupiter-mass objects by the accretion of pebbles and gas while they migrate through the disk. Our simulations show that the composition of the planetary gas atmosphere is dominated by the accretion of vapor that originates from inward drifting…

We study the impact of mass-transfer physics on the observable properties of binary black hole populations that formed through isolated binary evolution. We used the POSYDON framework to combine detailed MESA binary simulations with the COSMIC population synthesis tool to obtain an accurate estimate of merging binary black hole observables with a specific focus on the spins of the black holes. We investigate the impact of mass-accretion efficiency onto compact objects and common-envelope efficiency on the observed distributions of the effective inspiral spin parameter χeff, chirp mass Mchirp, and binary mass ratio q. We find that low common envelope efficiency translates to tighter orbits following the common envelope and therefore more tidally spun up second-born black holes. However, these systems have short merger timescales and are only marginally detectable by current gravitational-wave detectors as they form and merge at high redshifts (z ∼ 2), outside current detector horizons. Assuming Eddington-limited accretion efficiency and that the first-born black hole is formed with a…

We present constraints on the tensor-to-scalar ratio r using Planck data. We use the latest release of Planck maps, processed with the NPIPE code, which produces calibrated frequency maps in temperature and polarisation for all Planck channels from 30 GHz to 857 GHz using the same pipeline. We computed constraints on r using the BB angular power spectrum, and we also discuss constraints coming from the TT spectrum. Given Planck ’s noise level, the TT spectrum gives constraints on r that are cosmic-variance limited (with σ r = 0.093), but we show that the marginalised posterior peaks towards negative values of r at about the 1.2 σ level. We derived Planck constraints using the BB power spectrum at both large angular scales (the ‘reionisation bump’) and intermediate angular scales (the ‘recombination bump’) from ℓ = 2 to 150 and find a stronger constraint than that from TT , with σ r = 0.069. The Planck BB spectrum shows no systematic bias and is compatible with zero, given both the statistical noise and the systematic uncertainties. The likelihood analysis using B modes yields the…

We present redshift distribution estimates of galaxies selected from the fourth data release of the Kilo-Degree Survey over an area of ∼1000 deg2(KiDS-1000). These redshift distributions represent one of the crucial ingredients for weak gravitational lensing measurements with the KiDS-1000 data. The primary estimate is based on deep spectroscopic reference catalogues that are re-weighted with the help of a self-organising map (SOM) to closely resemble the KiDS-1000 sources, split into five tomographic redshift bins in the photometric redshift range 0.1 < zB ≤ 1.2. Sources are selected such that they only occupy that volume of nine-dimensional magnitude-space that is also covered by the reference samples (‘gold’ selection). Residual biases in the mean redshifts determined from this calibration are estimated from mock catalogues to be ≲0.01 for all five bins with uncertainties of ∼0.01. This primary SOM estimate of the KiDS-1000 redshift distributions is complemented with an independent clustering redshift approach. After validation of the clustering-zon the same mock catalogues and a…

Context. The CARMENES exoplanet survey of M dwarfs has obtained more than 18 000 spectra of 329 nearby M dwarfs over the past five years as part of its guaranteed time observations (GTO) program. Aims. We determine planet occurrence rates with the 71 stars from the GTO program for which we have more than 50 observations. Methods. We use injection-and-retrieval experiments on the radial-velocity time series to measure detection probabilities. We include 27 planets in 21 planetary systems in our analysis. Results. We find 0.06 −0.03 +0.04 giant planets (100 M ⊕ < M pl sin i < 1000 M ⊕ ) per star in periods of up to 1000 d, but due to a selection bias this number could be up to a factor of five lower in the whole 329-star sample. The upper limit for hot Jupiters (orbital period of less than 10 d) is 0.03 planets per star, while the occurrence rate of planets with intermediate masses (10 M ⊕ < M pl sin i < 100 M ⊕ ) is 0.18 −0.05 +0.07 planets per star. Less massive planets with 1 M ⊕ < M pl sin i < 10 M ⊕ are very abundant, with an estimated rate of 1.32 −0.31 +0.33 planets per star…

The majority of massive stars live in binary or multiple systems and will interact with a companion during their lifetimes, which helps to explain the observed diversity of core-collapse supernovae. Donor stars in binary systems can lose most of their hydrogen-rich envelopes through mass transfer. As a result, not only are the surface properties affected, but so is the core structure. However, most calculations of the core-collapse properties of massive stars rely on single-star models. We present a systematic study of the difference between the pre-supernova structures of single stars and stars of the same initial mass (11–21 M ⊙ ) that have been stripped due to stable post-main-sequence mass transfer at solar metallicity. We present the pre-supernova core composition with novel diagrams that give an intuitive representation of the isotope distribution. As shown in previous studies, at the edge of the carbon-oxygen core, the binary-stripped star models contain an extended gradient of carbon, oxygen, and neon. This layer remains until core collapse and is more extended in mass for…

Aims. We perform a comprehensive determination of the systemic proper motions of 74 dwarf galaxies and dwarf galaxy candidates in the Local Group based on Gaia early data release 3. The outputs of the analysis for each galaxy, including probabilities of membership, will be made publicly available. The analysis is augmented by a determination of the orbital properties of galaxies within 500 kpc. Methods. We adopt a flexible Bayesian methodology presented in the literature, which takes into account the location of the stars on the sky, on the colour-magnitude diagram, and on the proper motion plane. We applied some modifications, in particular to the way the colour-magnitude diagram and spectroscopic information are factored in, for example, by including stars in several evolution phases. The bulk motions were integrated in three gravitational potentials: two where the Milky Way was treated in isolation and has a mass 0.9 & 1.6 × 10 12 M ⊙ , and a time-varying potential, which includes the infall of a massive Large Magellanic Cloud (LMC). Results. We were able to determine bulk proper…

Context.Planetary formation and evolution is a combination of multiple interlinked processes. Constraining the mechanisms observationally requires statistical comparison to a large diversity of planetary systems. Aims.We want to understand global observable consequences of different physical processes (accretion, migration, and interactions) and initial properties (like disc masses and metallicities) on the demographics of the planetary population. We also want to study the convergence of our scheme with respect to one initial condition, the initial number of planetary embryo in each disc. Methods.We selected distributions of initial conditions that are representative of known protoplanetary discs. Then, we used the Generation III Bern model to perform planetary population synthesis. We synthesise five populations with each a different initial number of Moon-mass embryos per disc: 1, 10, 20, 50, and 100. The last is our nominal population consisting of 1000 stars (systems) that was used for an extensive statistical analysis of planetary systems around 1M⊙stars. Results.The…

We identify stellar structures in the PHANGS sample of 74 nearby galaxies and construct morphological masks of sub-galactic environments based on Spitzer 3.6 μ m images. At the simplest level, we distinguish five environments: centres, bars, spiral arms, interarm regions, and discs without strong spirals. Slightly more sophisticated masks include rings and lenses, which are publicly released but not explicitly used in this paper. We examine trends with environment in the molecular gas content, star formation rate, and depletion time using PHANGS–ALMA CO(2–1) intensity maps and tracers of star formation. The interarm regions and discs without strong spirals clearly dominate in area, whereas molecular gas and star formation are quite evenly distributed among the five basic environments. We reproduce the molecular Kennicutt–Schmidt relation with a slope compatible with unity within the uncertainties and without significant slope differences among environments. In contrast to what has been suggested by early studies, we find that bars are not always deserts devoid of gas and star…

Diffuse cluster radio sources, in the form of radio halos and relics, reveal the presence of cosmic rays and magnetic fields in the intracluster medium (ICM). These cosmic rays are thought to be (re)accelerated through the ICM turbulence and shock waves generated by cluster merger events. Here we characterize the presence of diffuse radio emission in known galaxy clusters in the HETDEX Spring Field, covering 424 deg 2 . For this, we developed a method to extract individual targets from LOFAR observations processed with the LoTSS DDF-pipeline software. This procedure enables improved calibration as well as the joint imaging and deconvolution of multiple pointings of selected targets. The calibration strategy can also be used for LOFAR low-band antenna and international-baseline observations. The fraction of Planck PSZ2 catalog clusters with any diffuse radio emission apparently associated with the ICM is 73 ± 17%. We detect a total of ten radio halos and twelve candidate halos in the HETDEX Spring Field. Of these ten radio halos, four are new discoveries, two of which are located in…

We present the first results of the Galaxy Activity, Torus, and Outflow Survey (GATOS), a project aimed at understanding the properties of the dusty molecular tori and their connection to the host galaxy in nearby Seyfert galaxies. Our project expands the range of active galactic nuclei (AGN) luminosities and Eddington ratios covered by previous surveys of Seyferts conducted by the Atacama Large Millimeter Array (ALMA), allowing us to study the gas feeding and feedback cycle in a combined sample of 19 Seyferts. We used ALMA to obtain new images of the emission of molecular gas and dust using the CO(3–2) and HCO + (4–3) lines as well as their underlying continuum emission at 870 μm with high spatial resolutions (0.1″ ∼ 7 − 13 pc) in the circumnuclear disks (CND) of ten nearby ( D < 28 Mpc) Seyfert galaxies selected from an ultra-hard X-ray survey. Our new ALMA observations detect 870 μm continuum and CO line emission from spatially resolved disks located around the AGN in all the sources. The bulk of the 870 μm continuum flux can be accounted for by thermal emission from dust in the…

Using the astrometry and integrated photometry from the Gaia Early Data Release 3, we map the density variations in the distribution of young upper main sequence (UMS) stars, open clusters, and classical Cepheids in the Galactic disc within several kiloparsecs of the Sun. We derive maps of relative over- and under-dense regions for UMS stars in the Galactic disc using both bivariate kernel density estimators and wavelet transformations. The resulting overdensity maps exhibit large-scale arches that extend in a clumpy but coherent way over the entire sampled volume, indicating the location of the spiral arm segments in the vicinity of the Sun. Peaks in the UMS overdensity are well matched by the distribution of young and intrinsically bright open clusters. By applying a wavelet transformation to a sample of classical Cepheids, we find that their overdensities possibly extend the spiral arm segments on a larger scale (≃10 kpc from the Sun). While the resulting map based on the UMS sample is generally consistent with previous models of the Sagittarius-Carina spiral arm, the geometry of…

Context. Outflows accelerated by active galactic nuclei (AGN) are commonly observed in the form of coherent, mildly collimated high-velocity gas directed along the AGN ionisation cones and kinetically powerful (≳10 44 − 45 erg s −1 ) jets. Recent works found that outflows can also be accelerated by low-power (≲10 44 erg s −1 ) jets, and the most recent cosmological simulations indicate that these are the dominant source of feedback on sub-kiloparsec scales, but little is known about their effect on the galaxy host. Aims. We study the relation between radio jets and the distribution and kinematics of the ionised gas in IC 5063 , NGC 5643 , NGC 1068, and NGC 1386 as part of our survey of nearby Seyfert galaxies called Measuring Active Galactic Nuclei Under MUSE Microscope (MAGNUM). All these objects host a small-scale (≲1 kpc) low-power (≲10 44 erg s −1 ) radio jet that has small inclinations (≲45°) with respect to the galaxy disc. Methods. We employed seeing-limited optical integral field spectroscopic observations from the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large…

Context. The nearest stars provide a fundamental constraint for our understanding of stellar physics and the Galaxy. The nearby sample serves as an anchor where all objects can be seen and understood with precise data. This work is triggered by the most recent data release of the astrometric space mission Gaia and uses its unprecedented high precision parallax measurements to review the census of objects within 10 pc. Aims. The first aim of this work was to compile all stars and brown dwarfs within 10 pc observable by Gaia and compare it with the Gaia Catalogue of Nearby Stars as a quality assurance test. We complement the list to get a full 10 pc census, including bright stars, brown dwarfs, and exoplanets. Methods. We started our compilation from a query on all objects with a parallax larger than 100 mas using the Set of Identifications, Measurements, and Bibliography for Astronomical Data database (SIMBAD). We completed the census by adding companions, brown dwarfs with recent parallax measurements not in SIMBAD yet, and vetted exoplanets. The compilation combines astrometry and…

Context. Open clusters (OCs) trace the evolution of the Galactic disc with great accuracy. Gaia and large ground-based spectroscopic surveys make it possible to determine their properties and study their kinematics with unprecedented precision. Aims. We study the kinematical behaviour of the OC population over time. We take advantage of the latest age determinations of OCs to investigate the correlations of the 6D phase-space coordinates and orbital properties with age. The phase-space distribution, age-velocity relation, and action distribution are compared to those of field stars. We also investigate the rotation curve of the Milky Way traced by OCs, and we compare it to that of other observational or theoretical studies. Methods. We gathered nearly 30 000 radial velocity (RV) measurements of OC members from both Gaia -RVS data and ground-based surveys and catalogues. We computed the weighted mean RV, Galactic velocities, and orbital parameters of 1382 OCs. We investigated their distributions as a function of age and by comparison to field stars. Results. We provide the largest RV…

Astronomical Roentgen Telescope – X-ray Concentrator (ART-XC) is the hard X-ray instrument with grazing incidence imaging optics on board the Spektr-Roentgen-Gamma (SRG) observatory. The SRG observatory is the flagship astrophysical mission of the Russian Federal Space Program, which was successively launched into orbit around the second Lagrangian point (L2) of the Earth-Sun system with a Proton rocket from the Baikonur cosmodrome on 13 July 2019. The ART-XC telescope will provide the first ever true imaging all-sky survey performed with grazing incidence optics in the 4–30 keV energy band and will obtain the deepest and sharpest map of the sky in the energy range of 4–12 keV. Observations performed during the early calibration and performance verification phase as well as during the ongoing all-sky survey that started on 12 December 2019 have demonstrated that the in-flight characteristics of the ART-XC telescope are very close to expectations based on the results of ground calibrations. Upon completion of its four-year all-sky survey, ART-XC is expected to detect approximately…

Context. Reliable predictions of mass-loss rates are important for massive-star evolution computations. Aims. We aim to provide predictions for mass-loss rates and wind-momentum rates of O-type stars, while carefully studying the behaviour of these winds as functions of stellar parameters, such as luminosity and metallicity. Methods. We used newly developed steady-state models of radiation-driven winds to compute the global properties of a grid of O-stars. The self-consistent models were calculated by means of an iterative solution to the equation of motion using full non-local thermodynamic equilibrium radiative transfer in the co-moving frame to compute the radiative acceleration. In order to study winds in different galactic environments, the grid covers main-sequence stars, giants, and supergiants in the Galaxy and both Magellanic Clouds. Results. We find a strong dependence of mass-loss on both luminosity and metallicity. Mean values across the grid are Ṁ~L*2.2 and Ṁ~L*0.95; however, we also find a somewhat stronger dependence on metallicity for lower luminosities. Similarly,…

Aims.The eROSITA Final Equatorial-Depth Survey has been carried out during the performance verification phase of the Spectrum-Roentgen-Gamma/eROSITA telescope and was completed in November 2019. This survey is designed to provide the first eROSITA-selected sample of clusters and groups and to test the predictions for the all-sky survey in the context of cosmological studies with clusters of galaxies. Methods.In the area of ~140 square degrees covered by eFEDS, 542 candidate clusters and groups of galaxies were detected as extended X-ray sources with theeSASSsource detection algorithm. We performed imaging and spectral analysis of the 542 cluster candidates with eROSITA X-ray data and studied the properties of the sample. Results.We provide the catalog of candidate galaxy clusters and groups detected by eROSITA in the eFEDS field down to a flux of ~10–14erg s–1cm–2in the soft band (0.5–2 keV) within 1’. The clusters are distributed in the redshift range ɀ=[0.01, 1.3] with a median redshift ɀmedian= 0.35. With eROSITA X-ray data, we measured the temperature of the intracluster medium…

Context. Since July 2014, the Gaia mission has been engaged in a high-spatial-resolution, time-resolved, precise, accurate astrometric, and photometric survey of the entire sky. Aims. We present the Gaia Science Alerts project, which has been in operation since 1 June 2016. We describe the system which has been developed to enable the discovery and publication of transient photometric events as seen by Gaia . Methods. We outline the data handling, timings, and performances, and we describe the transient detection algorithms and filtering procedures needed to manage the high false alarm rate. We identify two classes of events: (1) sources which are new to Gaia and (2) Gaia sources which have undergone a significant brightening or fading. Validation of the Gaia transit astrometry and photometry was performed, followed by testing of the source environment to minimise contamination from Solar System objects, bright stars, and fainter near-neighbours. Results. We show that the Gaia Science Alerts project suffers from very low contamination, that is there are very few false-positives. We…

Context. The distribution of member stars in the surroundings of an open cluster (OC) can shed light on the process of its formation, evolution, and dissolution. The analysis of structural parameters of OCs as a function of their age and position in the Galaxy constrains theoretical models of cluster evolution. The Gaia catalog is very appropriate for finding members of OCs at large distance from their centers. Aims. We revisit the membership lists of OCs from the solar vicinity, in particular, by extending these membership lists to the peripheral areas through Gaia EDR3. We then take advantage of these new member lists to study the morphological properties and the mass segregation levels of the clusters. Methods. We used the clustering algorithm HDBSCAN on Gaia parallaxes and proper motions to systematically search for members up to 50 pc from the cluster centers. We fit a King’s function on the radial density profile of these clusters and a Gaussian mixture model (GMM) on their two-dimensional member distribution to study their shape. We also evaluated the degree of mass…

MeerKAT’s large number (64) of 13.5 m diameter antennas, spanning 8 km with a densely packed 1 km core, create a powerful instrument for wide-area surveys, with high sensitivity over a wide range of angular scales. The MeerKAT Galaxy Cluster Legacy Survey (MGCLS) is a programme of long-track MeerKAT L -band (900−1670 MHz) observations of 115 galaxy clusters, observed for ∼6−10 h each in full polarisation. The first legacy product data release (DR1), made available with this paper, includes the MeerKAT visibilities, basic image cubes at ∼8″ resolution, and enhanced spectral and polarisation image cubes at ∼8″ and 15″ resolutions. Typical sensitivities for the full-resolution MGCLS image products range from ∼3−5 μJy beam −1 . The basic cubes are full-field and span 2° × 2°. The enhanced products consist of the inner 1.2° × 1.2° field of view, corrected for the primary beam. The survey is fully sensitive to structures up to ∼10′ scales, and the wide bandwidth allows spectral and Faraday rotation mapping. Relatively narrow frequency channels (209 kHz) are also used to provide H I…

We present constraints on extensions to the standard cosmological model of a spatially flat Universe governed by general relativity, a cosmological constant (Λ), and cold dark matter (CDM) by varying the spatial curvature ΩK, the sum of the neutrino masses ∑mν, the dark energy equation of state parameter w, and the Hu-Sawicki f(R) gravity fR0 parameter. With the combined 3 × 2 pt measurements of cosmic shear from the Kilo-Degree Survey (KiDS-1000), galaxy clustering from the Baryon Oscillation Spectroscopic Survey (BOSS), and galaxy-galaxy lensing from the overlap between KiDS-1000, BOSS, and the spectroscopic 2-degree Field Lensing Survey, we find results that are fully consistent with a flat ΛCDM model with ΩK = 0.011−0.057+0.054, ∑mν < 1.76 eV (95% CL), and w = −0.99−0.13+0.11. The fR0 parameter is unconstrained in our fully non-linear f(R) cosmic shear analysis. Considering three different model selection criteria, we find no clear preference for either the fiducial flat ΛCDM model or any of the considered extensions. In addition to extensions to the flat ΛCDM parameter space,…

We present a new calibration of the peak absolute magnitude of Type Ia supernovae (SNe Ia) based on the surface brightness fluctuations (SBF) method, aimed at measuring the value of the Hubble constant. We build a sample of calibrating anchors consisting of 24 SNe hosted in galaxies that have SBF distance measurements. Applying a hierarchical Bayesian approach, we calibrate the SN Ia peak luminosity and extend the Hubble diagram into the Hubble flow by using a sample of 96 SNe Ia in the redshift range 0.02 < z < 0.075, which was extracted from the Combined Pantheon Sample. We estimate a value of H0 = 70.50 ± 2.37 (stat.) ± 3.38 (sys.) km s−1 Mpc−1 (i.e., 3.4% stat., 4.8% sys.), which is in agreement with the value obtained using the tip of the red giant branch calibration. It is also consistent, within errors, with the value obtained from SNe Ia calibrated with Cepheids or the value inferred from the analysis of the cosmic microwave background. We find that the SNe Ia distance moduli calibrated with SBF are on average larger by 0.07 mag than those calibrated with Cepheids. Our…

Over the past decade, several works have used the ratio between total (rest 8−1000 μ m) infrared and radio (rest 1.4 GHz) luminosity in star-forming galaxies ( q IR ), often referred to as the infrared-radio correlation (IRRC), to calibrate the radio emission as a star formation rate (SFR) indicator. Previous studies constrained the evolution of q IR with redshift, finding a mild but significant decline that is yet to be understood. Here, for the first time, we calibrate q IR as a function of both stellar mass ( M ⋆ ) and redshift, starting from an M ⋆ -selected sample of > 400 000 star-forming galaxies in the COSMOS field, identified via ( NUV − r )/( r − J ) colours, at redshifts of 0.1 < z < 4.5. Within each ( M ⋆ , z ) bin, we stacked the deepest available infrared/sub-mm and radio images. We fit the stacked IR spectral energy distributions with typical star-forming galaxy and IR-AGN templates. We then carefully removed the radio AGN candidates via a recursive approach. We find that the IRRC evolves primarily with M ⋆ , with more massive galaxies displaying a systematically…

Context. The solar telescope connected to HARPS-N has been observing the Sun since the summer of 2015. Such a high-cadence, long-baseline data set is crucial for understanding spurious radial-velocity signals induced by our Sun and by the instrument. On the instrumental side, this data set allowed us to detect sub- m s −1 systematics that needed to be corrected for. Aims. The goals of this manuscript are to (i) present a new data reduction software for HARPS-N, (ii) demonstrate the improvement brought by this new software during the first three years of the HARPS-N solar data set, and (iii) release all the obtained solar products, from extracted spectra to precise radial velocities. Methods. To correct for the instrumental systematics observed in the data reduced with the current version of the HARPS-N data reduction software (DRS version 3.7), we adapted the newly available ESPRESSO DRS (version 2.2.3) to HARPS-N and developed new optimised recipes for the spectrograph. We then compared the first three years of HARPS-N solar data reduced with the current and new DRS. Results. The…

The hypothesis that the late Universe is isotropic and homogeneous is adopted by most cosmological studies, including studies of galaxy clusters. The cosmic expansion rate H 0 is thought to be spatially constant, while bulk flows are often presumed to be negligible compared to the Hubble expansion, even at local scales. The effects of bulk flows on the redshift–distance conversion are hence usually ignored. Any deviation from this consensus can strongly bias the results of such studies, and thus the importance of testing these assumptions cannot be understated. Scaling relations of galaxy clusters can be effectively used for this testing. In previous works, we observed strong anisotropies in cluster scaling relations, whose origins remain ambiguous. By measuring many different cluster properties, several scaling relations with different sensitivities can be built. Nearly independent tests of cosmic isotropy and large bulk flows are then feasible. In this work, we make use of up to 570 clusters with measured properties at X-ray, microwave, and infrared wavelengths to construct ten…

Context. A key piece of information to understand the origin and role of protoplanetary disk substructures is their dust content. In particular, disk substructures associated with gas pressure bumps can work as dust traps, accumulating grains and reaching the necessary conditions to trigger the streaming instability. Aims. In order to shed some light on the origin and role that disk substructures play in planet formation, we aim to characterize the dust content of substructures in the disk of TW Hya. Methods. We present Atacama Large Millimeter Array (ALMA) observations of TW Hya at 3.1 mm with ~50 milliarcsecond resolution. These new data were combined with archival high angular resolution ALMA observations at 0.87, 1.3, and 2.1 mm. We analyze these multiwavelength data to infer a disk radial profile of the dust surface density, maximum particle size, and slope of the particle size distribution. Results. Most previously known annular substructures in the disk of TW Hya are resolved at the four wavelengths. Inside the inner 3 au cavity, the 2.1 and 3.1 mm images show a compact…

Context. The Mid-Infrared Instrument (MIRI) onboard the James Webb Space Telescope (JWST) will provide imaging, coronagraphy, low-resolution spectroscopy, and medium-resolution spectroscopy at unprecedented sensitivity levels in the mid-infrared wavelength range. The Medium Resolution Spectrometer (MRS) of MIRI is an integral field spectrograph that provides diffraction-limited spectroscopy between 4.9 and 28.3 μm, within a field of view (FOV) varying from ∼13 to ∼56 arcsec square. The design for MIRI MRS conforms with the goals of the JWST mission to observe high redshift galaxies and to study cosmology as well as observations of galactic objects, and stellar and planetary systems. Aims. From ground testing, we calculate the physical parameters essential for general observers and calibrating the wavelength solution and resolving power of the MRS which is critical for maximizing the scientific performance of the instrument. Methods. We have used ground-based observations of discrete spectral features in combination with Fabry-Perot etalon spectra to characterize the wavelength…

We investigate the inner regions of the Milky Way using data from APOGEE and Gaia EDR3. Our inner Galactic sample has more than 26 500 stars within |XGal|< 5 kpc, |YGal|< 3.5 kpc, |ZGal|< 1 kpc, and we also carry out the analysis for a foreground-cleaned subsample of 8000 stars that is more representative of the bulge–bar populations. These samples allow us to build chemo-dynamical maps of the stellar populations with vastly improved detail. The inner Galaxy shows an apparent chemical bimodality in key abundance ratios [α/Fe], [C/N], and [Mn/O], which probe different enrichment timescales, suggesting a star formation gap (quenching) between the high- and low-α populations. Using a joint analysis of the distributions of kinematics, metallicities, mean orbital radius, and chemical abundances, we can characterize the different populations coexisting in the innermost regions of the Galaxy for the first time. The chemo-kinematic data dissected on an eccentricity–|Z|max plane reveal the chemical and kinematic signatures of the bar, the thin inner disc, and an inner thick disc, and a broad…

Context. The full third Gaia data release will provide, for the first time, the calibrated spectra obtained with the blue and red Gaia slitless spectrophotometers (BP and RP, respectively). Gaia is a very complex mission and cannot be considered as a single instrument, but rather as many instruments. The two lines of sight with wide fields of view introduce strong variations of the observations across the large focal plane with more than one hundred different detectors. The main challenge when facing Gaia spectral calibration is that no lamp spectra or flat fields are available during the mission. Also, the significant size of the line spread function with respect to the dispersion of the prisms produces alien photons contaminating neighbouring positions of the spectra. This makes the calibration special and different from standard approaches. Aims. This work gives a detailed description of the internal calibration model for the spectrophotometric data used to obtain the content of the Gaia catalogue. The main purpose of the internal calibration is to bring all the epoch spectra…

The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will scan thousands of square degrees of the northern sky with a unique set of 56 filters using the dedicated 2.55 m Javalambre Survey Telescope (JST) at the Javalambre Astrophysical Observatory. Prior to the installation of the main camera (4.2 deg2field-of-view with 1.2 Gpixels), the JST was equipped with the JPAS-Pathfinder, a one CCD camera with a 0.3 deg2field-of-view and plate scale of 0.23 arcsec pixel−1. To demonstrate the scientific potential of J-PAS, the JPAS-Pathfinder camera was used to perform miniJPAS, a ∼1 deg2survey of the AEGIS field (along the Extended Groth Strip). The field was observed with the 56 J-PAS filters, which include 54 narrow band (FWHM ∼ 145 Å) and two broader filters extending to the UV and the near-infrared, complemented by theu, g, r, iSDSS broad band filters. In this miniJPAS survey overview paper, we present the miniJPAS data set (images and catalogs), as we highlight key aspects and applications of these unique spectro-photometric data and describe how to access…

Context. The efficiency of the different processes responsible for the evolution of interstellar dust on the scale of a galaxy are, to date, very uncertain, spanning several orders of magnitude in the literature. Yet, precise knowledge of the grain properties is key to addressing numerous open questions about the physics of the interstellar medium and galaxy evolution. Aims. This article presents an empirical statistical study, aimed at quantifying the timescales of the main cosmic dust evolution processes as a function of the global properties of a galaxy. Methods. We modeled a sample of ≃800 nearby galaxies, spanning a wide range of metallicities, gas fractions, specific star formation rates, and Hubble stages. We derived the dust properties of each object from its spectral energy distribution. Through an additional level of analysis, we inferred the timescales of dust condensation in core-collapse supernova ejecta, grain growth in cold clouds, and dust destruction by shock waves. Throughout this paper, we have adopted a hierarchical Bayesian approach, resulting in a single large…

Aims. The 11-year solar cycle (Schwabe cycle) is the dominant pattern of solar magnetic activity reflecting the oscillatory dynamo mechanism in the Sun’s convection zone. Solar cycles have been directly observed since 1700, while indirect proxy data suggest their existence over a much longer period of time, but generally without resolving individual cycles and their continuity. Here we reconstruct individual solar cycles for the last millennium using recently obtained 14C data and state-of-the-art models. Methods. Starting with the 14C production rate determined from the so far most precise measurements of radiocarbon content in tree rings, solar activity was reconstructed in the following three physics-based steps: (1) correction of the 14C production rate for the changing geomagnetic field; (2) computation of the open solar magnetic flux; and (3) conversion into sunspot numbers outside of grand minima. All known uncertainties, including both measurement and model uncertainties, were straightforwardly accounted for by a Monte-Carlo method. Results. Cyclic solar activity is…

AU Microscopii (AU Mic) is a young, active star whose transiting planet was recently detected. Here, we report our analysis of its TESS light curve, where we modeled the BY Draconis type quasi-periodic rotational modulation by starspots simultaneously to the flaring activity and planetary transits. We measured a flare occurrence rate in AU Mic of 6.35 flares per day for flares with amplitudes in the range of 0.06% < f max < 1.5% of the star flux. We employed a Bayesian Markov chain Monte Carlo analysis to model the five transits of AU Mic b observed by TESS, improving the constraints on the planetary parameters. The measured planet-to-star effective radius ratio of R p ∕ R ⋆ = 0.0496 ± 0.0007 implies a physical radius of 4.07 ± 0.17 R ⊕ and a planet density of 1.4 ± 0.4 g cm −3 , confirming that AU Mic b is a Neptune-size moderately inflated planet. While a single feature possibly due to a second planet was previously reported in the former TESS data, we report the detection of two additional transit-like events in the new TESS observations of July 2020. This represents substantial…

Context. Until recently, the 3D shape, and therefore density (when combining the volume estimate with available mass estimates), and surface topography of the vast majority of the largest ( D ≥ 100 km) main-belt asteroids have remained poorly constrained. The improved capabilities of the SPHERE/ZIMPOL instrument have opened new doors into ground-based asteroid exploration. Aims. To constrain the formation and evolution of a representative sample of large asteroids, we conducted a high-angular-resolution imaging survey of 42 large main-belt asteroids with VLT/SPHERE/ZIMPOL. Our asteroid sample comprises 39 bodies with D ≥ 100 km and in particular most D ≥ 200 km main-belt asteroids (20/23). Furthermore, it nicely reflects the compositional diversity present in the main belt as the sampled bodies belong to the following taxonomic classes: A, B, C, Ch/Cgh, E/M/X, K, P/T, S, and V. Methods. The SPHERE/ZIMPOL images were first used to reconstruct the 3D shape of all targets with both the ADAM and MPCD reconstruction methods. We subsequently performed a detailed shape analysis and…

We present new colour–effective temperature ( T eff ) transformations based on the photometry of the early third data release (EDR3) of the ESA/ Gaia mission. These relations are calibrated on a sample of about 600 dwarf and giant stars for which T eff has previously been determined with the infrared flux method from dereddened colours. The 1 σ dispersion of the transformations is of 60–80 K for the pure Gaia colours (BP−RP) 0 , (BP− G ) 0 , and ( G −RP) 0 , improving to 40–60 K for colours including the 2MASS K s -band, namely (BP− K s )0 , (RP− K s ) 0 , and (G− K s ) 0 . We validate these relations in the most challenging case of dense stellar fields, where the Gaia EDR3 photometry could be less reliable, providing guidance for the safe use of Gaia colours in crowded environments. We compare the T eff from the Gaia EDR3 colours with those obtained from standard ( V − K s ) 0 colours for stars in three Galactic globular clusters of different metallicity, namely NGC 104, NGC 6752, and NGC 7099. The agreement between the two estimates of T eff is excellent, with mean differences of…

Aims. The catalog of Stars With ExoplanETs (SWEET-Cat) was originally introduced in 2013. Since then many more exoplanets have been confirmed, increasing significantly the number of host stars listed there. A crucial step toward a comprehensive understanding of these new worlds is the precise and homogeneous characterization of their host stars. Better spectroscopic stellar parameters along with new results from Gaia eDR3 provide updated and precise parameters for the discovered planets. A new version of the catalog, whose homogeneity in the derivation of the parameters is key to unraveling star–planet connections, is available to the community. Methods. We made use of high-resolution spectra for planet-host stars, either observed by our team or collected through public archives. The spectroscopic stellar parameters were derived for the spectra following the same homogeneous process using ARES and MOOG (ARES+MOOG) as for the previous SWEET-Cat releases. We re-derived parameters for the stars in the catalog using better quality spectra and/or using the most recent versions of the…

We present the OBELISK project, a cosmological radiation-hydrodynamics simulation that follows the assembly and reionization of a protocluster progenitor during the first two billion years after the big bang, down to z = 3.5. The simulation resolves haloes down to the atomic cooling limit and tracks the contribution of different sources of ionization: stars, active galactic nuclei, and collisions. The OBELISK project is specifically designed to study the coevolution of high-redshift galaxies and quasars in an environment favouring black hole growth. In this paper, we establish the relative contribution of these two sources of radiation to reionization and their respective role in establishing and maintaining the high-redshift ionizing background. Our volume is typical of an overdense region of the Universe and displays star formation rate and black hole accretion rate densities similar to those of high-redshift protoclusters. We find that hydrogen reionization happens inside-out, is completed by z ∼ 6 in our overdensity, and is predominantly driven by galaxies, while accreting black…

We determined effective temperatures, surface gravities, and metallicities for a sample of 343 M dwarfs observed with CARMENES, the double-channel, high-resolution spectrograph installed at the 3.5 m telescope at Calar Alto Observatory. We employed S TE P AR S YN , a Bayesian spectral synthesis implementation particularly designed to infer the stellar atmospheric parameters of late-type stars following a Markov chain Monte Carlo approach. We made use of the BT-Settl model atmospheres and the radiative transfer code turbospectrum to compute a grid of synthetic spectra around 75 magnetically insensitive Fe I and Ti I lines plus the TiO γ and ϵ bands. To avoid any potential degeneracy in the parameter space, we imposed Bayesian priors on T eff and log g based on the comprehensive, multi-band photometric data available for the sample. We find that this methodology is suitable down to M7.0 V, where refractory metals such as Ti are expected to condense in the stellar photospheres. The derived T eff , log g , and [Fe/H] range from 3000 to 4200 K, 4.5 to 5.3 dex, and −0.7 to 0.2 dex,…

Context.GaiaEarly Data Release 3 (GaiaEDR3) provides accurate astrometry for about 1.6 million compact (QSO-like) extragalactic sources, 1.2 million of which have the best-quality five-parameter astrometric solutions. Aims.The proper motions of QSO-like sources are used to reveal a systematic pattern due to the acceleration of the solar systembarycentre with respect to the rest frame of the Universe. Apart from being an important scientific result by itself, the acceleration measured in this way is a good quality indicator of theGaiaastrometric solution. Methods.Theeffect of the acceleration was obtained as a part of the general expansion of the vector field of proper motions in vector spherical harmonics (VSH). Various versions of the VSH fit and various subsets of the sources were tried and compared to get the most consistent result and a realistic estimate of its uncertainty. Additional tests with theGaiaastrometric solution were used to get a better idea of the possible systematic errors in the estimate. Results.Our best estimate of the acceleration based onGaiaEDR3 is (2.32 ±…

We report the detection of the sulfur-bearing species NCS, HCCS, H 2 CCS, H 2 CCCS, and C 4 S for the first time in space. These molecules were found towards TMC-1 through the observation of several lines for each species. We also report the detection of C 5 S for the first time in a cold cloud through the observation of five lines in the 31–50 GHz range. The derived column densities are N (NCS) = (7.8 ± 0.6) × 10 11 cm −2 , N (HCCS) = (6.8 ± 0.6) × 10 11 cm −2 , N (H 2 CCS) = (7.8 ± 0.8) × 10 11 cm −2 , N (H 2 CCCS) = (3.7 ± 0.4) × 10 11 cm −2 , N (C 4 S) = (3.8 ± 0.4) × 10 10 cm −2 , and N (C 5 S) = (5.0 ± 1.0) × 10 10 cm −2 . The observed abundance ratio between C 3 S and C 4 S is 340, that is to say a factor of approximately one hundred larger than the corresponding value for CCS and C 3 S. The observational results are compared with a state-of-the-art chemical model, which is only partially successful in reproducing the observed abundances. These detections underline the need to improve chemical networks dealing with S-bearing species.

In recent years, the advent of a new generation of radial velocity instruments has allowed us to detect planets with increasingly lower mass and to break the one Earth-mass barrier. Here we report a new milestone in this context by announcing the detection of the lowest-mass planet measured so far using radial velocities: L 98-59 b, a rocky planet with half the mass of Venus. It is part of a system composed of three known transiting terrestrial planets (planets b–d). We announce the discovery of a fourth nontransiting planet with a minimum mass of 3.06−0.37+0.33 M⊕and an orbital period of 12.796−0.019+0.020days and report indications for the presence of a fifth nontransiting terrestrial planet. With a minimum mass of 2.46−0.82+0.66 M⊕and an orbital period 23.15−0.17+0.60days, this planet, if confirmed, would sit in the middle of the habitable zone of the L 98-59 system. L 98-59 is a bright M dwarf located 10.6ṗc away. Positioned at the border of the continuous viewing zone of theJames WebbSpace Telescope, this system is destined to become a corner stone for comparative…

Planet-forming disks are not isolated systems. Their interaction with the surrounding medium affects their mass budget and chemical content. In the context of the ALMA-DOT program, we obtained high-resolution maps of assorted lines from six disks that are still partly embedded in their natal envelope. In this work, we examine the SO and SO 2 emission that is detected from four sources: DG Tau, HL Tau, IRAS 04302+2247, and T Tau. The comparison with CO, HCO + , and CS maps reveals that the SO and SO 2 emission originates at the intersection between extended streamers and the planet-forming disk. Two targets, DG Tau and HL Tau, offer clear cases of inflowing material inducing an accretion shock on the disk material. The measured rotational temperatures and radial velocities are consistent with this view. In contrast to younger Class 0 sources, these shocks are confined to the specific disk region impacted by the streamer. In HL Tau, the known accreting streamer induces a shock in the disk outskirts, and the released SO and SO 2 molecules spiral toward the star in a few hundred years.…

The cosmic radio dipole is of fundamental interest to studies of cosmology. Recent works have put forth open questions about the nature of the observed cosmic radio dipole. In the current work, we use simulated source count maps to test a linear and a quadratic estimator for possible biases in the estimated dipole amplitude with respect to the masking procedure. We find a superiority on the part of the quadratic estimator, which we used to analyse the TGSS-ADR1, WENSS, SUMSS, and NVSS radio source catalogues, spread over a decade of frequencies. We applied the same masking strategy to all four surveys to produce comparable results. In order to address the differences in the observed dipole amplitudes, we cross-matched the two surveys located at both ends of the analysed frequency range. For the linear estimator, we identified a general bias in the estimated dipole directions. The positional offsets of the quadratic estimator to the cosmic microwave background (CMB) dipole for skies with 107 simulated sources is found to be below one degree and the absolute accuracy of the estimated…

Star formation rate (SFR) measurements at z > 4 have relied mostly on the rest-frame far-ultraviolet (FUV) observations. The corrections for dust attenuation based on the IRX- β relation are highly uncertain and are still debated in the literature. Hence, rest-frame far-infrared (FIR) observations are necessary to constrain the dust-obscured component of the SFR. In this paper, we exploit the rest-frame FIR continuum observations collected by the ALMA Large Program to INvestigate [CII] at Early times (ALPINE) to directly constrain the obscured SFR in galaxies at 4.4 < z < 5.9. We used stacks of continuum images to measure average infrared luminosities taking both detected and undetected sources into account. Based on these measurements, we measured the position of the main sequence of star-forming galaxies and the specific SFR (sSFR) at z ∼ 4.5 and z ∼ 5.5. We find that the main sequence and sSFR do not significantly evolve between z ∼ 4.5 and z ∼ 5.5, as opposed to lower redshifts. We developed a method to derive the obscured SFR density (SFRD) using the stellar masses or…

Aims.Thanks to the high angular resolution, sensitivity, image fidelity, and frequency coverage of ALMA, we aim to improve our understanding of star formation. One of the breakthroughs expected from ALMA, which is the basis of our Cycle 5 ALMA-IMF Large Program, is the question of the origin of the initial mass function (IMF) of stars. Here we present the ALMA-IMF protocluster selection, first results, and scientific prospects. Methods.ALMA-IMF imaged a total noncontiguous area of ~53 pc2, covering extreme, nearby protoclusters of the Milky Way. We observed 15 massive (2.5 −33 × 103M⊙), nearby (2−5.5 kpc) protoclusters that were selected to span relevant early protocluster evolutionary stages. Our 1.3 and 3 mm observations provide continuum images that are homogeneously sensitive to point-like cores with masses of ~0.2M⊙and ~0.6M⊙, respectively, with a matched spatial resolution of ~2000 au across the sample at both wavelengths. Moreover, with the broad spectral coverage provided by ALMA, we detect lines that probe the ionized and molecular gas, as well as complex molecules. Taken…

Context. Recent observations by the Extreme Ultraviolet Imager (EUI) on board Solar Orbiter have characterized prevalent small-scale transient brightenings in the corona above the quiet Sun termed campfires. Aims. In this study we search for comparable brightenings in a numerical model and then investigate their relation to the magnetic field and the processes that drive these events. Methods. We used the MURaM code to solve the 3D radiation magnetohydrodynamic equations in a box that stretches from the upper convection zone to the corona. The model self-consistently produces a supergranular network of the magnetic field and a hot corona above this quiet Sun. For the comparison with the model, we synthesized the coronal emission as seen by EUI in its 174 Å channel, isolated the seven strongest transient brightenings, and investigated the changes of the magnetic field in and around these in detail. Results. The transients we isolated have a lifetime of about 2 min and are elongated loop-like features with lengths around 1 Mm to 4 Mm. They tend to occur at heights of about 2 Mm to 5…

Context. The analysis of the latest release of the Apache Point Observatory Galactic Evolution Experiment project (APOGEE DR16) data suggests the existence of a clear distinction between two sequences of disc stars at different Galactocentric distances in the [ α /Fe] versus [Fe/H] abundance ratio space: the so-called high- α sequence, classically associated with an old population of stars in the thick disc with high average [ α /Fe], and the low- α sequence, which mostly comprises relatively young stars in the thin disc with low average [ α /Fe]. Aims. We aim to constrain a multi-zone two-infall chemical evolution model designed for regions at different Galactocentric distances using measured chemical abundances from the APOGEE DR16 sample. Methods. We performed a Bayesian analysis based on a Markov chain Monte Carlo method to fit our multi-zone two-infall chemical evolution model to the APOGEE DR16 data. Results. An inside-out formation of the Galaxy disc naturally emerges from the best fit of our two-infall chemical-evolution model to APOGEE-DR16: Inner Galactic regions are…

The Cluster HEritage project withXMM-Newton– Mass Assembly and Thermodynamics at the Endpoint of structure formation (CHEX-MATE) is a three-mega-second Multi-Year Heritage Programme to obtain X-ray observations of a minimally-biased, signal-to-noise-limited sample of 118 galaxy clusters detected byPlanckthrough the Sunyaev–Zeldovich effect. The programme, described in detail in this paper, aims to study the ultimate products of structure formation in time and mass. It is composed of a census of the most recent objects to have formed (Tier-1: 0.05 < z < 0.2; 2 × 1014 M⊙ < M500 < 9 × 1014 M⊙), together with a sample of the highest mass objects in the Universe (Tier-2:z < 0.6;M500 > 7.25 × 1014 M⊙). The programme will yield an accurate vision of the statistical properties of the underlying population, measure how the gas properties are shaped by collapse into the dark matter halo, uncover the provenance of non-gravitational heating, and resolve the major uncertainties in mass determination that limit the use of clusters for cosmological parameter estimation. We will acquire X-ray…

We investigate the strongly lensed ( μ ≃ ×10 − 100) Lyman continuum (LyC) galaxy, dubbed Sunburst, at z = 2.37, taking advantage of a new accurate model of the lens. A characterization of the intrinsic (delensed) properties of the system yields a size of ≃3 sq. kpc, a luminosity of M UV = −20.3, and a stellar mass of M ≃ 10 9 M ⊙ ; 16% of the ultraviolet light is located in a 3 Myr old gravitationally bound young massive star cluster (YMC), with an effective radius of ∼8 pc (corresponding to 1 milliarcsec without lensing) and a dynamical mass of ∼10 7 M ⊙ (similar to the stellar mass) – from which LyC radiation is detected ( λ < 912 Å). The star formation rate and stellar mass surface densities for the YMC are Log 10 (Σ SFR [ M ⊙ yr −1 kpc −2 ]) ≃ 3.7 and Log 10 (Σ M [ M ⊙ pc −2 ]) ≃ 4.1, with sSFR > 330 Gyr −1 , consistent with the values observed in local young massive star clusters. The inferred outflowing gas velocity (> 300 km s −1 ) exceeds the escape velocity of the cluster. The resulting relative escape fraction of the ionizing radiation emerging from the entire galaxy is…

Context. The physical processes driving the formation of Galactic spiral arms are still under debate. Studies using open clusters favour the description of the Milky Way spiral arms as long-lived structures following the classical density wave theory. Current studies comparing the Gaia DR2 field stars kinematic information of the solar neighbourhood to simulations, find a better agreement with short-lived arms with a transient behaviour. Aims. Our aim is to provide an observational, data-driven view of the Milky Way spiral structure and its dynamics using open clusters as the main tracers, and to contrast it with simulation-based approaches. We used the most complete catalogue of Milky Way open clusters, with astrometric Gaia EDR3 updated parameters, estimated astrophysical information, and radial velocities, to revisit the nature of the spiral pattern of the Galaxy. Methods. We used a Gaussian mixture model to detect overdensities of open clusters younger than 30 Myr that correspond to the Perseus, Local, Sagittarius, and Scutum spiral Arms, respectively. We used the birthplaces of…

We report the detection of two isomers of ethynyl cyclopentadiene ( c -C 5 H 5 CCH), namely 1- and 2-ethynyl-1,3-cyclopentadiene, in the direction of TMC-1. We derive column densities of (1.4 ± 0.2) × 10 12 cm −2 and (2.0 ± 0.4) × 10 12 cm −2 , respectively, for these two cyclopentadiene derivatives, which imply that they are about ten times less abundant than cyclopentadiene. We also report the tentative detection of ethynyl benzene (C 6 H 5 CCH), for which we estimate a column density of (2.5 ± 0.4) × 10 12 cm −2 . We derived abundances for the corresponding cyano derivatives of cyclopentadiene and benzene and found values significantly lower than previously reported. The rotational temperature of the ethynyl and cyano derivatives of these cycles is about 9 K, that is, very close to the gas kinetic temperature of the cloud. The abundance ratio of the 1- and 2-isomers of ethynyl cyclopentadiene is 1.4 ± 0.5, while for the two isomers of cyano cyclopentadiene it is 2.4 ± 0.6. The relative abundances of CCH over CN derivatives is 7.7 ± 2.2 for cyclopentadiene, which probably reflects…

Context. The LOw Frequency ARray (LOFAR) is the only radio telescope that is presently capable of high-sensitivity, high-resolution (i.e. < 1 mJy beam −1 and < 15″) observations at ultra-low frequencies (< 100 MHz). To utilise these capabilities, the LOFAR Surveys Key Science Project is undertaking a large survey to cover the entire northern sky with Low Band Antenna (LBA) observations. Aims. The LOFAR LBA Sky Survey (LoLSS) aims to cover the entire northern sky with 3170 pointings in the frequency range between 42 − 66 MHz, at a resolution of 15″ and at a sensitivity of 1 mJy beam −1 (1 σ ). In this work, we outline the survey strategy, the observational status, and the calibration techniques. We also briefly describe several of our scientific motivations and present the preliminary public data release. Methods. The preliminary images were produced using a fully automated pipeline aimed at correcting all direction-independent effects in the data. Whilst the direction-dependent effects, such as those from the ionosphere, have not yet been corrected, the images presented in this work…

Aims.We aim to demonstrate the scientific potential of theGaiaEarly Data Release 3 (EDR3) for the study of different aspects of the Milky Way structure and evolution and we provide, at the same time, a description of several practical aspects of the data and examples of their usage. Methods.We used astrometric positions, proper motions, parallaxes, and photometry from EDR3 to select different populations and components and to calculate the distances and velocities in the direction of the anticentre. In this direction, theGaiaastrometric data alone enable the calculation of the vertical and azimuthal velocities; also, the extinction is relatively low compared to other directions in the Galactic plane. We then explore the disturbances of the current disc, the spatial and kinematical distributions of early accreted versus in situ stars, the structures in the outer parts of the disc, and the orbits of open clusters Berkeley 29 and Saurer 1. Results.With the improved astrometry and photometry of EDR3, we find that: (i) the dynamics of the Galactic disc are very complex with oscillations…

Aims.The complexity of star formation at the physical scale of molecular clouds is not yet fully understood. We investigate the mechanisms regulating the formation of stars in different environments within nearby star-forming galaxies from the Physics at High Angular resolution in Nearby GalaxieS (PHANGS) sample. Methods.Integral field spectroscopic data and radio-interferometric observations of 18 galaxies were combined to explore the existence of the resolved star formation main sequence (Σstellarversus ΣSFR), resolved Kennicutt–Schmidt relation (Σmol. gasversus ΣSFR), and resolved molecular gas main sequence (Σstellarversus Σmol. gas), and we derived their slope and scatter at spatial resolutions from 100 pc to 1 kpc (under various assumptions). Results.All three relations were recovered at the highest spatial resolution (100 pc). Furthermore, significant variations in these scaling relations were observed across different galactic environments. The exclusion of non-detections has a systematic impact on the inferred slope as a function of the spatial scale. Finally, the scatter…

We present CO(2−1) and adjacent continuum observations of seven nearby radio-quiet type-2 quasars (QSO2s) obtained with ALMA at ∼0.2″ resolution (370 pc at z ∼ 0.1). These QSO2s are luminous ( L [OIII] > 10 8.5 L ⊙ ∼ M B < −23), and their host galaxies massive ( M * ∼ 10 11 M ⊙ ). The CO morphologies are diverse, including disks and interacting systems. Two of the QSO2s are red early-type galaxies with no CO(2–1) detected. In the interacting galaxies, the central kiloparsec contains 18–25% of the total cold molecular gas, whereas in the spirals it is only ∼5–12%. J1010+0612 and J1430+1339 show double-peaked CO flux maps along the major axis of the CO disks that do not have an optical counterpart at the same angular resolution. Based on our analysis of the ionized and molecular gas kinematics and millimeter continuum emission, these CO morphologies are most likely produced by active galactic nucleus (AGN) feedback in the form of outflows, jets, and/or shocks. The CO kinematics of the QSO2s with CO(2−1) detections are dominated by rotation but also reveal noncircular motions.…

VANDELS is an ESO Public Spectroscopic Survey designed to build a sample of high-signal-to-noise ratio, medium-resolution spectra of galaxies at redshifts between 1 and 6.5. Here we present the final Public Data Release of the VANDELS Survey, comprising 2087 redshift measurements. We provide a detailed description of sample selection, observations, and data reduction procedures. The final catalogue reaches a target selection completeness of 40% at iAB = 25. The high signal-to-noise ratio of the spectra (above 7 in 80% of the spectra) and the dispersion of 2.5 Å allowed us to measure redshifts with high precision, the redshift measurement success rate reaching almost 100%. Together with the redshift catalogue and the reduced spectra, we also provide optical mid-infrared photometry and physical parameters derived through fitting the spectral energy distribution. The observed galaxy sample comprises both passive and star forming galaxies covering a stellar mass range of 8.3 < Log(M*/M⊙) < 11.7.

Intermediate-mass black holes (IMBHs) span the approximate mass range 100−105 M⊙, between black holes (BHs) that formed by stellar collapse and the supermassive BHs at the centers of galaxies. Mergers of IMBH binaries are the most energetic gravitational-wave sources accessible by the terrestrial detector network. Searches of the first two observing runs of Advanced LIGO and Advanced Virgo did not yield any significant IMBH binary signals. In the third observing run (O3), the increased network sensitivity enabled the detection of GW190521, a signal consistent with a binary merger of mass ∼150 M⊙ providing direct evidence of IMBH formation. Here, we report on a dedicated search of O3 data for further IMBH binary mergers, combining both modeled (matched filter) and model-independent search methods. We find some marginal candidates, but none are sufficiently significant to indicate detection of further IMBH mergers. We quantify the sensitivity of the individual search methods and of the combined search using a suite of IMBH binary signals obtained via numerical relativity, including…

Upcoming studies of extrasolar gas giants will give precise insights into the composition of planetary atmospheres, with the ultimate goal of linking it to the formation history of the planet. Here, we investigate how drifting and evaporating pebbles that enrich the gas phase of the disk influence the chemical composition of growing and migrating gas giants. To achieve this goal, we perform semi-analytical 1D models of protoplanetary disks, including viscous evolution, pebble drift, and evaporation, to simulate the growth of planets from planetary embryos to Jupiter-mass objects by the accretion of pebbles and gas while they migrate through the disk. The gas phase of the protoplanetary disk is enriched due to the evaporation of inward drifting pebbles crossing evaporation lines, leading to the accretion of large amounts of volatiles into the planetary atmosphere. As a consequence, gas-accreting planets are enriched in volatiles (C, O, N) compared to refractories (e.g., Mg, Si, Fe) by up to a factor of 100, depending on the chemical species, its exact abundance and volatility, and…

Strong lensing time delays can measure the Hubble constant H0 independently of any other probe. Assuming commonly used forms for the radial mass density profile of the lenses, a 2% precision has been achieved with seven Time-Delay Cosmography (TDCOSMO) lenses, in tension with the H0 from the cosmic microwave background. However, without assumptions on the radial mass density profile – and relying exclusively on stellar kinematics to break the mass-sheet degeneracy – the precision drops to 8% with the current data obtained using the seven TDCOSMO lenses, which is insufficient to resolve the H0 tension. With the addition of external information from 33 Sloan Lens ACS (SLACS) lenses, the precision improves to 5% if the deflectors of TDCOSMO and SLACS lenses are drawn from the same population. We investigate the prospect of improving the precision of time-delay cosmography without relying on mass profile assumptions to break the mass-sheet degeneracy. Our forecasts are based on a previously published hierarchical framework. With existing samples and technology, 3.3% precision on H0 can…

We have discovered a 300 kpc-wide giant Lyman-α (Lyα) nebula centered on the massive galaxy group RO-1001 at z = 2.91 in the Cosmic Evolution Survey field. Keck Cosmic Web Imager observations reveal three cold gas filaments converging into the center of the potential well of its ∼4 × 1013 M⊙ dark matter halo, hosting 1200 M⊙ yr−1 of star formation as probed by Atacama Large Millimeter Array and NOrthern Extended Millimeter Array observations. The nebula morphological and kinematics properties and the prevalence of blueshifted components in the Lyα spectra are consistent with a scenario of gas accretion. The upper limits on active galactic nuclei activity and overall energetics favor gravity as the primary Lyα powering source and infall as the main source of gas flows to the system. Although interpretational difficulties remain, with outflows and likely also photoionization with ensuing recombination still playing a role, this finding provides arguably an ideal environment to quantitatively test models of cold gas accretion and galaxy feeding inside an actively star-forming massive…

Context. Results from observations report a growing number of metal-poor stars showing an abundance pattern midway between the s - and r -processes. These so-called r/s-stars raise the need for an intermediate neutron capture process ( i -process), which is thought to result from the ingestion of protons in a convective helium-burning region, but whose astrophysical site is still largely debated. Aims. We investigate whether an i -process during the asymptotic giant branch (AGB) phase of low-metallicity low-mass stars can develop and whether it can explain the abundances of observed r/s-stars. Methods. We computed a 1 M ⊙ model at [Fe/H] = −2.5 with the stellar evolution code STAREVOL, using a nuclear network of 1091 species (at maximum) coupled to the transport processes. The impact of the temporal and spatial resolutions on the resulting abundances was assessed. We also identified key elements and isotopic ratios that are specific to i -process nucleosynthesis and carried out a detailed comparison between our model and a sample of r/s-stars. Results. At the beginning of the AGB…

Context. The majority of massive stars are part of binary systems. In about a quarter of these, the companions are so close that mass transfer occurs while they undergo core hydrogen burning, first on the thermal and then on the nuclear timescale. The nuclear timescale mass transfer leads to observational counterparts: the semi-detached so-called massive Algol binaries. These systems may provide urgently needed tests of the physics of mass transfer. However, comprehensive model predictions for these systems are sparse. Aims. We use a large grid of detailed evolutionary models of short-period massive binaries and follow-up population synthesis calculations to derive probability distributions of the observable properties of massive Algols and their descendants. Methods. Our results are based on ∼10 000 binary model sequences calculated with the stellar evolution code MESA, using a metallicity suitable for the Large Magellanic Cloud (LMC), covering initial donor masses between 10 M ⊙ and 40 M ⊙ and initial orbital periods above 1.4 d. These models include internal differential rotation…

Submillimeter/millimeter observations of dusty star-forming galaxies with the Atacama Large Millimeter/submillimeter Array (ALMA) have shown that dust continuum emission generally occurs in compact regions smaller than the stellar distribution. However, it remains to be understood how systematic these findings are. Studies often lack homogeneity in the sample selection, target discontinuous areas with inhomogeneous sensitivities, and suffer from modest u v coverage coming from single array configurations. GOODS-ALMA is a 1.1 mm galaxy survey over a continuous area of 72.42 arcmin 2 at a homogeneous sensitivity. In this version 2.0, we present a new low resolution dataset and its combination with the previous high resolution dataset from the survey, improving the u v coverage and sensitivity reaching an average of σ = 68.4 μJy beam −1 . A total of 88 galaxies are detected in a blind search (compared to 35 in the high resolution dataset alone), 50% at S / N peak ≥ 5 and 50% at 3.5 ≤ S / N peak ≤ 5 aided by priors. Among them, 13 out of the 88 are optically dark or faint sources ( H -…

The International LOFAR Telescope is an interferometer with stations spread across Europe. With baselines of up to ~2000 km, LOFAR has the unique capability of achieving sub-arcsecond resolution at frequencies below 200 MHz. However, it is technically and logistically challenging to process LOFAR data at this resolution. To date only a handful of publications have exploited this capability. Here we present a calibration strategy that builds on previous high-resolution work with LOFAR. It is implemented in a pipeline using mostly dedicated LOFAR software tools and the same processing framework as the LOFAR Two-metre Sky Survey (LoTSS). We give an overview of the calibration strategy and discuss the special challenges inherent to enacting high-resolution imaging with LOFAR, and describe the pipeline, which is publicly available, in detail. We demonstrate the calibration strategy by using the pipeline on P205+55, a typical LoTSS pointing with an 8 h observation and 13 international stations. We perform in-field delay calibration, solution referencing to other calibrators in the field,…

We present the first identification in interstellar space of the propargyl radical (CH 2 CCH). This species was observed in the cold dark cloud TMC-1 using the Yebes 40 m telescope. The six strongest hyperfine components of the 2 0, 2 –1 0, 1 rotational transition, lying at 37.46 GHz, were detected with signal-to-noise ratios from 4.6 σ to 12.3 σ . We derived a column density of 8.7 × 10 13 cm −2 for CH 2 CCH, which translates to a fractional abundance relative to H 2 of 8.7 × 10 −9 . This radical has a similar abundance as methyl acetylene, with an abundance ratio CH 2 CCH/CH 3 CCH close to one. The propargyl radical is thus one of the most abundant radicals detected in TMC-1, and it is probably the most abundant organic radical with a certain chemical complexity ever found in a cold dark cloud. We constructed a gas-phase chemical model and find calculated abundances that agree with, or fall two orders of magnitude below, the observed value depending on the poorly constrained low-temperature reactivity of CH 2 CCH with neutral atoms. According to the chemical model, the propargyl…

Massive starburst galaxies in the early Universe are estimated to have depletion times of ∼100 Myr and thus be able to convert their gas very quickly into stars, possibly leading to a rapid quenching of their star formation. For these reasons, they are considered progenitors of massive early-type galaxies (ETGs). In this paper, we study two high- z starbursts, AzTEC/C159 ( z ≃ 4.57) and J1000+0234 ( z ≃ 4.54), observed with ALMA in the [C II ] 158- μ m emission line. These observations reveal two massive and regularly rotating gaseous discs. A 3D modelling of these discs returns rotation velocities of about 500 km s −1 and gas velocity dispersions as low as ≈ 20 km s −1 , leading to very high ratios between regular and random motion ( V / σ ≳ 20), at least in AzTEC/C159. The mass decompositions of the rotation curves show that both galaxies are highly baryon-dominated with gas masses of ≈10 11 M ⊙ , which, for J1000+0234, is significantly higher than previous estimates. We show that these high- z galaxies overlap with z = 0 massive ETGs in the ETG analogue of the stellar-mass…

Aims. We aim to demonstrate that the presence and mass of an exoplanet can now be effectively derived from the astrometry of another exoplanet. Methods. We combined previous astrometry of β Pictoris b with a new set of observations from the GRAVITY interferometer. The orbital motion of β Pictoris b is fit using Markov chain Monte Carlo simulations in Jacobi coordinates. The inner planet, β Pictoris c, was also reobserved at a separation of 96 mas, confirming the previous orbital estimations. Results. From the astrometry of planet b only, we can (i) detect the presence of β Pictoris c and (ii) constrain its mass to 10.04 −3.10 +4.53 M Jup . If one adds the astrometry of β Pictoris c, the mass is narrowed down to 9.15 −1.06 +1.08 M Jup . The inclusion of radial velocity measurements does not affect the orbital parameters significantly, but it does slightly decrease the mass estimate to 8.89 −0.75 +0.75 M Jup . With a semimajor axis of 2.68 ± 0.02 au, a period of 1221 ± 15 days, and an eccentricity of 0.32 ± 0.02, the orbital parameters of β Pictoris c are now constrained as precisely…

Context. Data from ground-based, high-resolution solar telescopes can only be used for science with calibrations and processing, which requires detailed knowledge about the instrumentation. Space-based solar telescopes provide science-ready data, which are easier to work with for researchers whose expertise is in the interpretation of data. Recently, data-processing pipelines for ground-based instruments have been constructed. Aims. We aim to provide observers with a user-friendly data pipeline for data from the Swedish 1-meter Solar Telescope (SST) that delivers science-ready data together with the metadata needed for proper interpretation and archiving. Methods. We briefly describe the CHROMospheric Imaging Spectrometer (CHROMIS) instrument, including its (pre)filters, as well as recent upgrades to the CRisp Imaging SpectroPolarimeter (CRISP) prefilters and polarization optics. We summarize the processing steps from raw data to science-ready data cubes in FITS files. We report calibrations and compensations for data imperfections in detail. Misalignment of Ca II data due to…

Context. A large-scale magnetic field permeates our Galaxy and is involved in a variety of astrophysical processes, such as star formation and cosmic ray propagation. Dust polarization has been proven to be one of the most powerful observables for studying the field properties in the interstellar medium (ISM). However, it does not provide a direct measurement of its strength. Different methods have been developed that employ both polarization and spectroscopic data in order to infer the field strength. The most widely applied method was developed by Davis (1951, Phys. Rev., 81, 890) and Chandrasekhar & Fermi (1953, ApJ, 118, 1137), hereafter DCF. The DCF method relies on the assumption that isotropic turbulent motions initiate the propagation of Alfvén waves. Observations, however, indicate that turbulence in the ISM is anisotropic and that non-Alfvénic (compressible) modes may be important. Aims. Our goal is to develop a new method for estimating the field strength in the ISM that includes the compressible modes and does not contradict the anisotropic properties of turbulence.…

Aims. Planet-forming discs are believed to be very weakly turbulent in the regions outside of 1 AU. For this reason, it is now believed that magnetised winds could be the dominant mechanism driving accretion in these systems. However, currently, no self-consistent approach can describe discs that are subject to a magnetised wind in a way similar to the α disc model. In this article, I explore the parameter space of wind-driven protoplanetary discs in a systematic manner and present scaling laws that can be used in reduced models in a similar way to α disc models. Methods. I computed a series of self-similar wind solutions, assuming that the disc is dominated by ambipolar and Ohmic diffusion. These solution were obtained by searching for stationary solutions in the finite-volume code PLUTO using a relaxation method and continuation. Results. Self-similar solutions are obtained for values of plasma β ranging from 10 2 to 10 8 for several Ohmic and ambipolar diffusion strengths. Mass accretion rates of about 10 −8 M ⊙ yr −1 are obtained for the poloidal field strength β = O (10 4 ) or…

Aims. As a part of the CALYPSO Large Program, we aim to constrain the properties of protostellar jets and outflows by searching for corresponding emission in a sample of 21 Class 0 protostars covering a broad range of internal luminosities (Lint from 0.035 L⊙ to 47 L⊙). Methods. We analyzed high angular (~0.′′5–1.′′0) resolution IRAM Plateau de Bure Interferometer (PdBI) observations in three typical outflow and jet tracers, namely: CO (2−1), SO (56−45), and SiO (5−4). We obtained the detection rate, spatial distribution, kinematics, and collimation of the outflow and jet emission in the three lines. Molecular column densities and abundances, the jet mass-loss and momentum rates, and mechanical luminosity are estimated from the integrated line intensities. Results. Blue- and red-shifted emission in CO (2−1), which probes outflowing gas, was detected in all the sources in the sample and detected for the first time in the cases of SerpS-MM22 and SerpS-MM18b. Collimated high-velocity jets in SiO (5−4) were detected in 67% of the sources (for the first time in IRAS4B2, IRAS4B1,…

We present new Atacama Large Millimeter/Submillimeter Array observations of the [C II ] 158 μm transition and the dust continuum in HZ4, a typical star-forming galaxy when the Universe was only ∼1 Gyr old ( z ≈ 5.5). Our high ≈0.3″ spatial resolution allows us to study the relationships between [C II ] line emission, star formation rate, and far-infrared emission on spatial scales of ∼2 kpc. In the central ∼4 kpc of HZ4, the [C II ]/FIR is ∼3 × 10 −3 on global scales as well as on spatially resolved scales of ∼2 kpc, comparable to the ratio observed in local moderate starburst galaxies such as M 82 or M 83. For the first time in an individual normal galaxy at this redshift, we find evidence for outflowing gas from the central star-forming region in the direction of the minor axis of the galaxy. The projected velocity of the outflow is ∼400 km s −1 , and the neutral gas-mass outflow rate is ∼3 − 6 times higher than the star formation rate in the central region. Finally, we detect a diffuse component of [C II ] emission, or [C II ] halo, that extends beyond the star-forming disk and…

Context. Blue straggler stars are exotic objects present in all stellar environments whose nature and formation channels are still partially unclear. They seem to be particularly abundant in open clusters (OCs), thus offering a unique chance to tackle these problems statistically. Aims. We aim to build up a new and homogeneous catalogue of blue straggler stars (BSS) in Galactic OCs using Gaia to provide a more solid assessment of the membership of these stars. We also aim to explore possible relationships of the straggler abundance with the parent cluster’s structural and dynamical parameters. As a by-product, we also search for possible yellow straggler stars (YSS), which are believed to be stragglers in a more advanced evolution stage. Methods. We employed photometry, proper motions, and parallaxes extracted from Gaia DR2 for 408 Galactic star clusters and searched for stragglers within them after performing a careful membership analysis. Results. The number of BBS emerging from our more stringent, selection criteria turns out to be significantly smaller than in previous versions…

Context. The interstellar medium is the locus of physical processes affecting the evolution of galaxies which drive or are the result of star formation activity, supermassive black hole growth, and feedback. The resulting physical conditions determine the observable chemical abundances that can be explored through molecular emission observations at millimeter and submillimeter wavelengths. Aims. Our goal is to unveiling the molecular richness of the central region of the prototypical nearby starburst galaxy NGC 253 at an unprecedented combination of sensitivity, spatial resolution, and frequency coverage. Methods. We used the Atacama Large Millimeter/submillimeter Array (ALMA), covering a nearly contiguous 289 GHz frequency range between 84.2 and 373.2 GHz, to image the continuum and spectral line emission at 1.6″(∼28 pc) resolution down to a sensitivity of 30 − 50 mK. This article describes the ALMA Comprehensive High-resolution Extragalactic Molecular Inventory (ALCHEMI) large program. We focus on the analysis of the spectra extracted from the 15″ (∼255 pc) resolution ALMA Compact…

The treatment and criteria for development of unstable Roche lobe overflow (RLOF) that leads to the common envelope (CE) phase have hindered the area of evolutionary predictions for decades. In particular, the formation of black hole-black hole (BH-BH), black hole-neutron star (BH-NS), and neutron star-neutron star (NS-NS) merging binaries depends sensitively on the CE phase in classical isolated binary evolution model. All these mergers are now reported as LIGO/Virgo sources or source candidates. CE is even considered by some as a mandatory phase in the formation of BH-BH, BH-NS, or NS-NS mergers in binary evolution models. At the moment, there is no full first-principles model for the development of the CE. We employed the StarTrack population synthesis code to test the current advancements in studies on the stability of RLOF for massive donors to assess their effect on the LIGO/Virgo source population. In particular, we allowed for more restrictive CE development criteria for massive donors ( M > 18 M ⊙ ). We also tested a modified condition for switching between different types…

Context. Herbig Ae/Be stars (HAeBes) have so far been studied based on relatively small samples that are scattered throughout the sky. Their fundamental stellar and circumstellar parameters and statistical properties were derived with heterogeneous approaches before Gaia. Aims. Our main goal is to contribute to the study of HAeBes from the largest sample of such sources to date, for which stellar and circumstellar properties have been determined homogeneously from the analysis of the spectral energy distributions (SEDs) and Gaia EDR3 parallaxes and photometry. Methods. Multiwavelength photometry was compiled for 209 bona fide HAeBes for which Gaia EDR3 distances were estimated. Using the Virtual Observatory SED Analyser (VOSA), photospheric models were fit to the optical SEDs to derive stellar parameters, and the excesses at infrared (IR) and longer wavelengths were characterized to derive several circumstellar properties. A statistical analysis was carried out to show the potential use of such a large dataset. Results. The stellar temperature, luminosity, radius, mass, and age were…

Sub-Chandrasekhar mass carbon-oxygen white dwarfs with a surface helium shell have been proposed as progenitors of Type Ia supernovae (SNe Ia). If true, the resulting thermonuclear explosions should be able to account for at least some of the range of SNe Ia observables. To study this, we conducted a parameter study based on three-dimensional simulations of double detonations in carbon-oxygen white dwarfs with a helium shell, assuming different core and shell masses. An admixture of carbon to the shell and solar metallicity are included in the models. The hydrodynamic simulations were carried out using the A REPO code. This allowed us to follow the helium shell detonation with high numerical resolution, and this improves the reliability of predicted nucleosynthetic shell detonation yields. The addition of carbon to the shell leads to a lower production of 56 Ni, while including solar metallicity increases the production of intermediate mass elements. The production of higher mass elements is further shifted to stable isotopes at solar metallicity. Moreover, we find different core…

Using VLTI/GRAVITY and SINFONI data, we investigate the subparsec gas and dust structure around the nearby type 1 active galactic nucleus (AGN) hosted by NGC 3783. The K-band coverage of GRAVITY uniquely allows simultaneous analysis of the size and kinematics of the broad line region (BLR), the size and structure of the near-infrared(near-IR)-continuum-emitting hot dust, and the size of the coronal line region (CLR). We find the BLR, probed through broad Brγ emission, to be well described by a rotating, thick disc with a radial distribution of clouds peaking in the inner region. In our BLR model, the physical mean radius of 16 light-days is nearly twice the ten-day time-lag that would be measured, which closely matches the ten-day time-lag that has been measured by reverberation mapping. We measure a hot dust full-width at half-maximum (FWHM) size of 0.74 mas (0.14 pc) and further reconstruct an image of the hot dust, which reveals a faint (5% of the total flux) offset cloud that we interpret as an accreting or outflowing cloud heated by the central AGN. Finally, we directly measure…

Context.Molecular filaments and hubs have received special attention recently thanks to new studies showing their key role in star formation. While the (column) density and velocity structures of both filaments and hubs have been carefully studied, their magnetic field (B-field) properties have yet to be characterized. Consequently, the role of B-fields in the formation and evolution of hub-filament systems is not well constrained. Aims.We aim to understand the role of the B-field and its interplay with turbulence and gravity in the dynamical evolution of the NGC 6334 filament network that harbours cluster-forming hubs and high-mass star formation. Methods.We present new observations of the dust polarized emission at 850μm toward the 2 pc × 10 pc map of NGC 6334 at a spatial resolution of 0.09 pc obtained with theJames Clerk MaxwellTelescope (JCMT) as part of the B-field In STar-forming Region Observations (BISTRO) survey. We study the distribution and dispersion of the polarized intensity (PI), the polarization fraction (PF), and the plane-of-the-sky B-field angle (χB_POS) toward…

Astronomy & Astrophysics has the great pleasure of publishing a special issue on the early installment of the third data release (EDR3) of the Gaia space mission which has been obtaining exquisite astrometric and photometric measurements since July 2014.Similar to the full third data release planned for 2022, Gaia EDR3 is based on the analysis on the first 34 months of the mission.It provides parallaxes, proper motions, and three-color photometry for 1.5 billion sources, which will all be unchanged in the full third release that will only add information to EDR3.

Context.The recent early third data release (EDR3) from theGaiamission has produced parallaxes for 1.468 × 109sources with better quality than those reported in the previous data release. Nevertheless, there are calibration issues with the data that require corrections to the published values and uncertainties. Aims.We want to properly characterize the behavior of the random and systematic uncertainties of theGaiaEDR3 parallaxes in order to maximize the precision of the derived distances without compromising their accuracy. We also aim to provide a step-by-step procedure for the calculation of distances to stars and stellar clusters when using these parallaxes. Methods.We reanalyzed some of the data presented in the calibration papers for quasar and Large Magellanic Cloud (LMC) parallaxes and combine these results with measurements for six bright globular clusters. We calculated the angular covariance of EDR3 parallaxes at small separations (up to a few degrees) based on the LMC results and combined it with the results for larger angles using quasars to obtain an approximate…