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

We present the results of the search for an isotropic stochastic gravitational wave background (GWB) at nanohertz frequencies using the second data release of the European Pulsar Timing Array (EPTA) for 25 millisecond pulsars and a combination with the first data release of the Indian Pulsar Timing Array (InPTA). A robust GWB detection is conditioned upon resolving the Hellings-Downs angular pattern in the pairwise cross-correlation of the pulsar timing residuals. Additionally, the GWB is expected to yield the same (common) spectrum of temporal correlations across pulsars, which is used as a null hypothesis in the GWB search. Such a common-spectrum process has already been observed in pulsar timing data. We analysed (i) the full 24.7-year EPTA data set, (ii) its 10.3-year subset based on modern observing systems, (iii) the combination of the full data set with the first data release of the InPTA for ten commonly timed millisecond pulsars, and (iv) the combination of the 10.3-year subset with the InPTA data. These combinations allowed us to probe the contributions of instrumental…

We present JADES JWST/NIRSpec spectroscopy of GN-z11, the most luminous candidate z > 10 Lyman break galaxy in the GOODS-North field with M UV = −21.5. We derive a redshift of z = 10.603 (lower than previous determinations) based on multiple emission lines in our low and medium resolution spectra over 0.7 − 5.3 μm. We significantly detect the continuum and measure a blue rest-UV spectral slope of β = −2.4. Remarkably, we see spatially extended Lyman- α in emission (despite the highly neutral intergalactic medium expected at this early epoch), offset 555 km s −1 redwards of the systemic redshift. From our measurements of collisionally excited lines of both low and high ionisation (including [O II ] λ 3727, [Ne III ] λ 3869, and C III ] λ 1909), we infer a high ionisation parameter (log U ∼ −2). We detect the rarely seen N IV ] λ 1486 and N III ] λ 1748 lines in both our low and medium resolution spectra, with other high ionisation lines seen in the low resolution spectrum, such as He II (blended with O III ]) and C IV (with a possible P-Cygni profile). Based on the observed rest-UV…

Context. Data from the Gaia satellite are revolutionising our understanding of the Milky Way. With every new data release, there is a need to update the census of open clusters. Aims. We aim to conduct a blind, all-sky search for open clusters using 729 million sources from Gaia DR3 down to magnitude G ∼ 20, creating a homogeneous catalogue of clusters including many new objects. Methods. We used the Hierarchical Density-Based Spatial Clustering of Applications with Noise (HDBSCAN) algorithm to recover clusters. We validated our clusters using a statistical density test and a Bayesian convolutional neural network for colour-magnitude diagram classification. We inferred basic astrometric parameters, ages, extinctions, and distances for the clusters in the catalogue. Results. We recovered 7167 clusters, 2387 of which are candidate new objects and 4782 of which crossmatch to objects in the literature, including 134 globular clusters. A more stringent cut of our catalogue contains 4105 highly reliable clusters, 739 of which are new. Owing to the scope of our methodology, we are able to…

Pulsar timing arrays offer a probe of the low-frequency gravitational wave spectrum (1–100 nHz), which is intimately connected to a number of markers that can uniquely trace the formation and evolution of the Universe. We present the dataset and the results of the timing analysis from the second data release of the European Pulsar Timing Array (EPTA). The dataset contains high-precision pulsar timing data from 25 millisecond pulsars collected with the five largest radio telescopes in Europe, as well as the Large European Array for Pulsars. The dataset forms the foundation for the search for gravitational waves by the EPTA, presented in associated papers. We describe the dataset and present the results of the frequentist and Bayesian pulsar timing analysis for individual millisecond pulsars that have been observed over the last ~25 yr. We discuss the improvements to the individual pulsar parameter estimates, as well as new measurements of the physical properties of these pulsars and their companions. This data release extends the dataset from EPTA Data Release 1 up to the beginning…

We present rest-frame optical data of the compact z = 5.55 galaxy GS_3073 obtained using the integral field spectroscopy mode of the Near-InfraRed Spectrograph on board the James Webb Space Telescope. The galaxy’s prominent broad components in several hydrogen and helium lines (though absent in the forbidden lines) and v detection of a large equivalent width of He II λ 4686, EW(He II ) ∼20 Å, unambiguously identify it as an active galactic nucleus (AGN). We measured a gas phase metallicity of Z gas / Z ⊙ ∼0.21 −0.04 +0.08 , which is lower than what has been inferred for both more luminous AGN at a similar redshift and lower redshift AGN. We empirically show that classical emission line ratio diagnostic diagrams cannot be used to distinguish between the primary ionisation source (AGN or star formation) for systems with such low metallicity, though different diagnostic diagrams involving He II λ 4686 prove very useful, independent of metallicity. We measured the central black hole mass to be log( M BH / M ⊙ )∼8.2 ± 0.4 based on the luminosity and width of the broad line region of the…

Context . The Medium-Resolution Spectrometer (MRS) provides one of the four operating modes of the Mid-Infrared Instrument (MIRI) on board the James Webb Space Telescope (JWST). The MRS is an integral field spectrometer, measuring the spatial and spectral distributions of light across the 5–28 µm wavelength range with a spectral resolving power between 3700 and 1300. Aims . We present the MRS’s optical, spectral, and spectro-photometric performance, as achieved in flight, and we report on the effects that limit the instrument’s ultimate sensitivity. Methods . The MRS flight performance has been quantified using observations of stars, planetary nebulae, and planets in our Solar System. The precision and accuracy of this calibration was checked against celestial calibrators with well-known flux levels and spectral features. Results . We find that the MRS geometric calibration has a distortion solution accuracy relative to the commanded position of 8 mas at 5 µm and 23 mas at 28 µm. The wavelength calibration is accurate to within 9 km s −1 at 5 µm and 27 km s −1 at 28 µm. The…

We present cosmological parameter constraints using maps from the last Planck data release (PR4). In particular, we detail an upgraded version of the cosmic microwave background likelihood, HiLLiPoP, that is based on angular power spectra and relies on a physical modeling of the foreground residuals in the spectral domain. This new version of the likelihood retains a larger sky fraction (up to 75%) and uses an extended multipole range. Using this likelihood, along with low-ℓ measurements from LoLLiPoP, we derived constraints on ΛCDM parameters that are in good agreement with previous Planck 2018 results, but with smaller uncertainties by 10% to 20%. We demonstrate that the foregrounds can be accurately described in the spectral domain, with a negligible impact on ΛCDM parameters. We also derived constraints on single-parameter extensions to ΛCDM, including AL, ΩK, Neff, and ∑mν. Noteworthy results from this updated analysis include a lensing amplitude value of AL = 1.039 ± 0.052, which is more closely aligned with theoretical expectations within the ΛCDM framework. Additionally, our…

Context. How galaxies form, assemble, and cease their star formation is a central question within the modern landscape of galaxy evolution studies. These processes are indelibly imprinted on the galaxy stellar mass function (SMF), and its measurement and understanding is key to uncovering a unified theory of galaxy evolution. Aims. We present constraints on the shape and evolution of the galaxy SMF, the quiescent galaxy fraction, and the cosmic stellar mass density across 90% of the history of the Universe from z = 7.5 → 0.2 as a means to study the physical processes that underpin galaxy evolution. Methods. The COSMOS survey is an ideal laboratory for studying representative galaxy samples. Now equipped with deeper and more homogeneous near-infrared coverage exploited by the COSMOS2020 catalog, we leverage the large 1.27 deg2 effective area to improve sample statistics and understand spatial variations (cosmic variance) – particularly for rare, massive galaxies – and push to higher redshifts with greater confidence and mass completeness than previous studies. We divide the total…

We present emission-line ratios from a sample of 27 Lyman-break galaxies from z ∼ 5.5 − 9.5 with −17.0 < M1500 < −20.4, measured from ultra-deep JWST/NIRSpec multi-object spectroscopy from the JWST Advanced Deep Extragalactic Survey (JADES). We used a combination of 28 h deep PRISM/CLEAR and 7 h deep G395M/F290LP observations to measure, or place strong constraints on, ratios of widely studied rest-frame optical emission lines including Hα, Hβ, [O II] λλ3726, 3729, [Ne III] λ3869, [O III] λ4959, [O III] λ5007, [O I] λ6300, [N II] λ6583, and [S II] λλ6716, 6731 in individual z > 5.5 spectra. We find that the emission-line ratios exhibited by these z ∼ 5.5 − 9.5 galaxies occupy clearly distinct regions of line-ratio space compared to typical z ∼ 0 − 3 galaxies, instead being more consistent with extreme populations of lower-redshift galaxies. This is best illustrated by the [O III]/[O II] ratio, tracing interstellar medium (ISM) ionisation, in which we observe more than half of our sample to have [O III]/[O II] > 10. Our high signal-to-noise spectra reveal more than an order of…

Context. Gaia has been in operations since 2014, and two full data releases (DR) have been delivered so far: DR1 in 2016 and DR2 in 2018. The third Gaia data release expands from the early data release (EDR3) in 2020, which contained the five-parameter astrometric solution and mean photometry for 1.8 billion sources by providing 34 months of multi-epoch observations that allowed us to systematically probe, characterise, and classify variable celestial phenomena. Aims. We present a summary of the variability processing and analysis of the photometric and spectroscopic time series of 1.8 billion sources carried out for Gaia DR3. Methods. We used statistical and machine learning methods to characterise and classify the variable sources. Training sets were built from a global revision of major published variable star catalogues. For a subset of classes, specific detailed studies were conducted to confirm their class membership and to derive parameters that are adapted to the peculiarity of the considered class. Results. In total, 10.5 million objects are identified as variable in Gaia…

Context. Computational astrophysics nowadays routinely combines grid-adaptive capabilities with modern shock-capturing, high resolution spatio-temporal integration schemes in challenging multidimensional hydrodynamic and magnetohydrodynamic (MHD) simulations. A large, and still growing, body of community software exists, and we provide an update on recent developments within the open-source MPI-AMRVAC code. Aims. Complete with online documentation, the MPI-AMRVAC 3.0 release includes several recently added equation sets and offers many options to explore and quantify the influence of implementation details. While showcasing this flexibility on a variety of hydro-dynamic and MHD tests, we document new modules of direct interest for state-of-the-art solar applications. Methods. Test cases address how higher-order reconstruction strategies impact long-term simulations of shear layers, with and without gas-dust coupling effects, how runaway radiative losses can transit to intricate multi-temperature, multiphase dynamics, and how different flavors of spatio-temporal schemes and/or…

Context. Mid-infrared observations of photodissociation regions (PDRs) are dominated by strong emission features called aromatic infrared bands (AIBs). The most prominent AIBs are found at 3.3, 6.2, 7.7, 8.6, and 11.2 µm. The most sensitive, highest-resolution infrared spectral imaging data ever taken of the prototypical PDR, the Orion Bar, have been captured by JWST. These high-quality data allow for an unprecedentedly detailed view of AIBs. Aims. We provide an inventory of the AIBs found in the Orion Bar, along with mid-IR template spectra from five distinct regions in the Bar: the molecular PDR (i.e. the three H2dissociation fronts), the atomic PDR, and the HIIregion. Methods. We used JWST NIRSpec IFU and MIRI MRS observations of the Orion Bar from the JWST Early Release Science Program, PDRs4All (ID: 1288). We extracted five template spectra to represent the morphology and environment of the Orion Bar PDR. We investigated and characterised the AIBs in these template spectra. We describe the variations among them here. Results. The superb sensitivity and the spectral and spatial…

Context . Detailed astrochemical models are a key component to interpret the observations of interstellar and circumstellar molecules since they allow important physical properties of the gas and its evolutionary history to be deduced. Aims . We update one of the most widely used astrochemical databases to reflect advances in experimental and theoretical estimates of rate coefficients and to respond to the large increase in the number of molecules detected in space since our last release in 2013. Methods . We present the sixth release of the UMIST Database for Astrochemistry (UDfA), a major expansion of the gas-phase chemistry that describes the synthesis of interstellar and circumstellar molecules. Since our last release, we have undertaken a major review of the literature which has increased the number of reactions by over 40% to a total of 8767 and increased the number of species by over 55% to 737. We have made a particular attempt to include many of the new species detected in space over the past decade, including those from the QUIJOTE and GOTHAM surveys, as well as providing…

X-ray quasi-periodic eruptions (QPEs) represent a recently discovered example of extreme X-ray variability associated with supermassive black holes. These are high-amplitude bursts recurring every few hours that are detected in the soft X-ray band from the nuclei of nearby galaxies whose optical spectra lack the broad emission lines typically observed in unobscured active galaxies. The physical origin of this new X-ray variability phenomenon is still unknown and several theoretical models have been presented. However, no attempt has been made so far to account for the varying QPE recurrence time and luminosity in individual sources, nor for the diversity of the QPE phenomenology in the different known erupters. We present a semi-analytical model based on an extreme mass-ratio inspiral (EMRI) system where the secondary intersects, along its orbit, a rigidly precessing accretion disc surrounding the primary. We assume that QPEs result from emission from an adiabatically expanding, initially optically thick gas cloud expelled from the disc plane at each impact. We produced synthetic…

Context. The astrometric discovery of sub-stellar mass companions orbiting stars is exceedingly hard due to the required sub-milliarcsecond precision, limiting the application of this technique to only a few instruments on a target-per-target basis and to the global astrometry space missions H IPPARCOS and Gaia . The third Gaia data release ( Gaia DR3) includes the first Gaia astrometric orbital solutions whose sensitivity in terms of estimated companion mass extends down to the planetary-mass regime. Aims. We present the contribution of the exoplanet pipeline to the Gaia DR3 sample of astrometric orbital solutions by describing the methods used for fitting the orbits, the identification of significant solutions, and their validation. We then present an overview of the statistical properties of the solution parameters. Methods. Using both a Markov chain Monte Carlo and a genetic algorithm, we fitted the 34 months of Gaia DR3 astrometric time series with a single Keplerian astrometric-orbit model that had 12 free parameters and an additional jitter term, and retained the solutions…

Context. Gaia Data Release 3 (DR3) provides a number of new data products that complement the early DR3 made available two years ago. Among these is the first Gaia catalogue of eclipsing-binary candidates containing 2 184 477 sources with brightnesses from a few magnitudes to 20 mag in the Gaia G -band and covering the full sky. Aims. We present the catalogue, describe its content, provide tips for its use, estimate its quality, and show illustrative samples. Methods. Candidate selection is based on the results of variable object classification performed within the Gaia Data Processing and Analysis Consortium. Candidates are then further filtered using eclipsing-binary-tailored criteria based on the G -band light curves. To find the orbital period, a large ensemble of trial periods is first acquired using three distinct period-search methods applied to the cleaned G light curve of each source. The G light curve is then modelled with up to two Gaussians and a cosine for each trial period. The best combination of orbital period and geometric model is finally selected using Bayesian…

The exoplanet population characterized by relatively short orbital periods (P < 100 d) around solar-type stars is dominated by super-Earths and sub-Neptunes. However, these planets are missing in our Solar System and the reason behind this absence is still unknown. Two theoretical scenarios invoke the role of Jupiter as the possible culprit: Jupiter may have acted as a dynamical barrier to the inward migration of sub-Neptunes from beyond the water iceline; alternatively, Jupiter may have considerably reduced the inward flux of material (pebbles) required to form super-Earths inside that iceline. Both scenarios predict an anti-correlation between the presence of small planets and that of cold Jupiters in exoplanetary systems. To test that prediction, we homogeneously analyzed the radial-velocity measurements of 38 Kepler and K2 transiting small planet systems gathered over nearly ten years with the HARPS-N spectrograph, as well as publicly available radial velocities collected with other facilities. We used Bayesian differential evolution Markov chain Monte Carlo techniques, which in…

The Gaia mission has provided an invaluable wealth of astrometric data for more than a billion stars in our Galaxy. The synergy between Gaia astrometry, photometry, and spectroscopic surveys gives us comprehensive information about the Milky Way. Using the Bayesian isochrone-fitting code StarHorse , we derive distances and extinctions for more than 10 million unique stars listed in both Gaia Data Release 3 and public spectroscopic surveys: 557 559 in GALAH+ DR3, 4 531 028 in LAMOST DR7 LRS, 347 535 in LAMOST DR7 MRS, 562 424 in APOGEE DR17, 471 490 in RAVE DR6, 249 991 in SDSS DR12 (optical spectra from BOSS and SEGUE), 67 562 in the Gaia -ESO DR5 survey, and 4 211 087 in the Gaia RVS part of the Gaia DR3 release. StarHorse can increase the precision of distance and extinction measurements where Gaia parallaxes alone would be uncertain. We used StarHorse for the first time to derive stellar ages for main-sequence turnoff and subgiant branch stars, around 2.5 million stars, with age uncertainties typically around 30%; the uncertainties drop to 15% for subgiant-branch-only stars,…

Context. Gaia DR3 contains 1.8 billion sources with G -band photometry, 1.5 billion of which with G BP and G RP photometry, complemented by positions on the sky, parallax, and proper motion. The median number of field-of-view transits in the three photometric bands is between 40 and 44 measurements per source and covers 34 months of data collection. Aims. We pursue a classification of Galactic and extra-galactic objects that are detected as variable by Gaia across the whole sky. Methods. Supervised machine learning (eXtreme Gradient Boosting and Random Forest) was employed to generate multi-class, binary, and meta-classifiers that classified variable objects with photometric time series in the G , G BP , and G RP bands. Results. Classification results comprise 12.4 million sources (selected from a much larger set of potential variable objects) and include about 9 million variable stars classified into 22 variability types in the Milky Way and nearby galaxies such as the Magellanic Clouds and Andromeda, plus thousands of supernova explosions in distant galaxies, 1 million active…

Context. LBQS 0302−0019 is a blue quasar (QSO) at z ∼ 3.3 that hosts powerful outflows and resides in a complex environment consisting of an obscured active galactic nucleus (AGN) candidate and multiple companions, all within 30 kpc in projection. Aims. We aim to characterise this complex system using JWST NIRSpec Integral Field Spectrograph (IFS) observations obtained as part of the NIRSpec IFS GTO programme “Galaxy Assembly with NIRSpec IFS” (GA-NIFS); these data cover the QSO rest-frame optical emission lines with a spatial resolution of ∼0.1″ and a sampling of 0.05″ (∼380 pc) over a contiguous sky area of ∼3″ × 3″ (23 × 23 kpc 2 ). Methods. We developed a procedure to correct for the spurious oscillations (or “wiggles”) in NIRSpec single-spaxel spectra caused by the spatial under-sampling of the point spread function. We performed a QSO–host decomposition with the QD EBLEND 3D tools. We used multi-component kinematic decomposition of the optical emission line profiles to infer the physical properties of the emitting gas in the QSO environment. Results. The QSO–host decomposition…

Unusually high N/O abundance ratios were recently reported for a very compact, intensively star-forming object GN-z11 at z = 10.6 from JWST/NIRSpec observations. We present an empirical comparison with the C, N, and O abundance ratios in Galactic globular clusters (GCs) over a large metallicity range. We show that hot hydrogen-burning nucleosynthesis within supermassive stars (SMS) formed through runaway collisions can consistently explain the observed abundances ratio in GN-z11 and in GCs. This suggests that a proto-globular cluster hosting a SMS could be at the origin of the strong N-enrichment in GN-z11. Our model predicts the behavior of N/O, C/O, and Ne/O ratios as a function of metallicity, which can be tested if high- z objects similar to GN-z11 are detected with JWST in the future. Further studies and statistics will help differentiate the proto-GC scenario from the Wolf-Rayet scenario that we quantify with a population synthesis model, and shed more light on this peculiar object.

Contact. The star formation rate (SFR) in high-redshift galaxies is expected due to competing physical processes. This stochastic variability might boost the luminosity of galaxies and might explain the over-abundance seen at z ≳ 10 by the James Webb Space Telescope. Aims. We quantify the amplitude and timescales of this variability and identify the key physical processes. Methods. We selected 245 z = 7.7 galaxies with stellar mass 5 × 10 6 ≲ M ⋆ / M ⊙ ≲ 5 × 10 10 from SERRA , which is a suite of high-resolution radiation-hydrodynamic cosmological simulations. After fitting the average SFR trend, ⟨SFR⟩, we quantified the time-dependent variation, δ ( t )≡log[SFR/⟨SFR⟩], for each system and performed a periodogram analysis to search for periodicity modulations. Results. We find that δ ( t ) is distributed as a zero-mean Gaussian, with standard deviation σ δ ≃ 0.24 (corresponding to a UV magnitude s.d. σ UV ≃ 0.61) that is independent of M ⋆ . However, the modulation timescale increases with stellar mass: t δ ∼ (9, 50, 100) Myr for M ⋆ ∼ (0.1, 1, 5)×10 9 M ⊙ , respectively. These…

Our position inside the Galactic disc has previously prevented us from establishing an accurate rotation curve (RC). The advent of Gaia and its third data release ( Gaia DR3) made it possible to specify the RC up to twice the optical radius. We aim to establish a new RC of the Galaxy from the Gaia DR3 by drastically reducing systematic uncertainties. Our goal is to provide a new estimate of the mass of the Galaxy. We compared different estimates, established a robust assessment of the systematic uncertainties, and addressed differences in methodologies, particularly regarding distance estimates. We find a sharply decreasing RC for the Milky Way; the decrease in velocity between 19.5 and 26.5 kpc is approximately 30 km s −1 . We identify, for the first time, a Keplerian decline of the RC, starting at ∼19 kpc and ending at ∼26.5 kpc from the Galaxy centre, while a flat RC is rejected with a significance of 3 σ . The total mass is revised downwards to 2.06 −0.13 +0.24 × 10 11 M ⊙ , which is in agreement with the absence of a significant mass increase at radii larger than 19 kpc. We…

Context. Recent James Webb Space Telescope (JWST) spectroscopic observations of the z = 10.6 galaxy GN-z11 have revealed a very peculiar UV spectrum exhibiting intense emission lines of nitrogen, which are not typically detected in galaxy spectra. This observation indicates a super-solar N/O abundance ratio at low metallicity, which only resembles the abundances seen in globular cluster (GC) stars. This discovery suggests that we might be seeing proto-GCs in formation or possibly even signatures of supermassive stars. Aims. To examine whether other objects with strong N IV and/or N III emission lines (N-emitters, hereafter) exist and to better understand their origin and nature, we have examined the available JWST spectra and data from the literature. Methods. Using the NIRSpec/JWST observations from CEERS, we found an extreme N-emitter, CEERS-1019 at z = 8.6782, showing intense N IV] λ1486 and N III] λ1750 emission. From the observed rest-UV and optical lines, we conclude that it is compatible with photoionization from stars and we have determined accurate abundances for C, N, O,…

Strong-lensing time delays enable the measurement of the Hubble constant ( H 0 ) independently of other traditional methods. The main limitation to the precision of time-delay cosmography is mass-sheet degeneracy (MSD). Some of the previous TDCOSMO analyses broke the MSD by making standard assumptions about the mass density profile of the lens galaxy, reaching 2% precision from seven lenses. However, this approach could potentially bias the H 0 measurement or underestimate the errors. For this work, we broke the MSD for the first time using spatially resolved kinematics of the lens galaxy in RXJ1131−1231 obtained from the Keck Cosmic Web Imager spectroscopy, in combination with previously published time delay and lens models derived from Hubble Space Telescope imaging. This approach allowed us to robustly estimate H 0 , effectively implementing a maximally flexible mass model. Following a blind analysis, we estimated the angular diameter distance to the lens galaxy D d = 865 −81 +85 Mpc and the time-delay distance D Δt = 2180 −271 +472 Mpc, giving H 0 = 77.1 −7.1 +7.3 km s −1 Mpc −1…

Aims. The nanohertz gravitational wave background (GWB) is expected to be an aggregate signal of an ensemble of gravitational waves emitted predominantly by a large population of coalescing supermassive black hole binaries in the centres of merging galaxies. Pulsar timing arrays (PTAs), which are ensembles of extremely stable pulsars at approximately kiloparsec distances precisely monitored for decades, are the most precise experiments capable of detecting this background. However, the subtle imprints that the GWB induces on pulsar timing data are obscured by many sources of noise that occur on various timescales. These must be carefully modelled and mitigated to increase the sensitivity to the background signal. Methods. In this paper, we present a novel technique to estimate the optimal number of frequency coefficients for modelling achromatic and chromatic noise, while selecting the preferred set of noise models to use for each pulsar. We also incorporated a new model to fit for scattering variations in the Bayesian pulsar timing package temponest. These customised noise models…

We reconstructed the star formation history of the Sco-Cen OB association using a novel high-resolution age map of the region. We developed an approach to produce robust ages for Sco-Cen’s recently identified 37 stellar clusters using the SigMA algorithm. The Sco-Cen star formation timeline reveals four periods of enhanced star formation activity, or bursts, remarkably separated by about 5 Myr. Of these, the second burst, which occurred about 15 million years ago, is by far the dominant one, and most of Sco-Cen’s stars and clusters were in place by the end of this burst. The formation of stars and clusters in Sco-Cen is correlated but not linearly, implying that more stars were formed per cluster during the peak of the star formation rate. Most of the clusters that are large enough to have supernova precursors were formed during the second burst around 15 Myr ago. Star and cluster formation activity has been continuously declining since then. We have clear evidence that Sco-Cen formed from the inside out and that it contains 100-pc long chains of contiguous clusters exhibiting…

We present new astrometric and polarimetric observations of flares from Sgr A* obtained with GRAVITY, the near-infrared interferometer at ESO’s Very Large Telescope Interferometer (VLTI), bringing the total sample of well-covered astrometric flares to four and polarimetric flares to six. Of all flares, two are well covered in both domains. All astrometric flares show clockwise motion in the plane of the sky with a period of around an hour, and the polarization vector rotates by one full loop in the same time. Given the apparent similarities of the flares, we present a common fit, taking into account the absence of strong Doppler boosting peaks in the light curves and the EHT-measured geometry. Our results are consistent with and significantly strengthen our model from 2018. First, we find that the combination of polarization period and measured flare radius of around nine gravitational radii (9 R g ≈ 1.5 R ISCO , innermost stable circular orbit) is consistent with Keplerian orbital motion of hot spots in the innermost accretion zone. The mass inside the flares’ radius is consistent…

We present ten novel [OIII]λ4363 auroral line detections up to z ∼ 9.5 measured from ultra-deep JWST/NIRSpec MSA spectroscopy from the JWST Advanced Deep Extragalactic Survey (JADES). We leverage the deepest spectroscopic observations taken thus far with NIRSpec to determine electron temperatures and oxygen abundances using the direct Te method. We directly compare these results against a suite of locally calibrated strong-line diagnostics and recent high-z calibrations. We find the calibrations fail to simultaneously match our JADES sample, thus warranting a self-consistent revision of these calibrations for the high-z Universe. We find a weak dependence between R2 and O3O2 with metallicity, thus suggesting these line ratios are inefficient in the high-z Universe as metallicity diagnostics and degeneracy breakers. We find R3 and R23 are still correlated with metallicity, but we find a tentative flattening of these diagnostics, thus suggesting future difficulties when applying these strong line ratios as metallicity indicators in the high-z Universe. We also propose and test an…

Quasar feedback in the form of powerful outflows is invoked as a key mechanism to quench star formation in galaxies, although direct observational evidence is still scarce and debated. Here we present Early Release Science JWST NIRSpec IFU observations of the z = 1.59 prototypical obscured Active Galactic Nucleus (AGN) XID2028: This target represents a unique test case for studying quasar feedback at the peak epoch of AGN-galaxy co-evolution because extensive multi-wavelength coverage is available and a massive and extended outflow is detected in the ionised and molecular components. With the unprecedented sensitivity and spatial resolution of the JWST, the NIRSpec dataset reveals a wealth of structures in the ionised gas kinematics and morphology that were previously hidden in the seeing-limited ground-based data. In particular, we find evidence of an interaction between the interstellar medium of the galaxy and the quasar-driven outflow and radio jet that produces an expanding bubble from which the fast and extended wind detected in previous observations emerges. The new…

Context. High-mass X-ray binaries (HMXBs) are a particular class of high-energy sources that require multi-wavelength observational efforts to be properly characterised. New identifications and the refinement of previous measurements are regularly published in the literature by independent teams of researchers and might, when they are collected in a catalogue, offer a tool for facilitating further studies of HMXBs. Aims. We update previous instances of HMXB catalogues in the Galaxy and provide the community easy access to the most complete set of observables on Galactic HMXBs. In addition to the fixed version that is available in Vizier, we also aim to host and maintain a dynamic version that can be updated upon request from users. Any modification will be logged in this version. Methods. Using previous HMXB catalogues supplemented by listings of hard X-ray sources detected in the past 20 yr, we produced a base set of HMXBs and candidates by means of identifier and sky coordinate cross matches. We queried in Simbad for unreferenced HMXBs. We searched for as many hard X-ray, soft…

The CRyogenic InfraRed Echelle Spectrograph (CRIRES) Upgrade project CRIRES + extended the capabilities of CRIRES. It transformed this VLT instrument into a cross-dispersed spectrograph to increase the wavelength range that is covered simultaneously by up to a factor of ten. In addition, a new detector focal plane array of three Hawaii 2RG detectors with a 5.3 μm cutoff wavelength replaced the existing detectors. Amongst many other improvements, a new spectropolarimetric unit was added and the calibration system has been enhanced. The instrument was installed at the VLT on Unit Telescope 3 at the beginning of 2020 and successfully commissioned and verified for science operations during 2021, partly remotely from Europe due to the COVID-19 pandemic. The instrument was subsequently offered to the community from October 2021 onwards. This article describes the performance and capabilities of the upgraded instrument and presents on sky results.

The escape fraction of Lyman-continuum (LyC) photons ( f esc ) is a key parameter for determining the sources of cosmic reionization at z ≥ 6. At these redshifts, owing to the opacity of the intergalactic medium, the LyC emission cannot be measured directly. However, LyC leakers during the epoch of reionization could be identified using indirect indicators that have been extensively tested at low and intermediate redshifts. These include a high [O III ]/[O II ] flux ratio, high star-formation surface density, and compact sizes. In this work, we present observations of 29 4.5 ≤ z ≤ 8 gravitationally lensed galaxies in the Abell 2744 cluster field. From a combined analysis of JWST-NIRSpec and NIRCam data, we accurately derived their physical and spectroscopic properties: our galaxies have low masses (log( M ⋆ )∼8.5), blue UV spectral slopes ( β ∼ −2.1), compact sizes ( r e ∼ 0.3 − 0.5 kpc), and high [O III ]/[O II ] flux ratios. We confirm that these properties are similar to those characterizing low-redshift LyC leakers. Indirectly inferring the fraction of escaping ionizing photons,…

Aims. Using the direct-imaging technique, we searched for low-mass companions around the star AF Lep, which presents a significant proper-motion anomaly (PMa) signal obtained from the comparison of H IPPARCOS and Gaia eDR3 catalogs. Methods. We observed AF Lep in two epochs with VLT/SPHERE using its subsystems IFS and IRDIS in the near-infrared, covering wavelengths ranging from the Y to the K spectral bands (between 0.95 and 2.3 μm). We then reduced the data using the high-contrast imaging techniques angular differential imaging (ADI) and spectral differential imaging in order to be able to retrieve the signal from low-mass companions of the star. Results. A faint companion was retrieved at a separation of ~0.335″ from the star and with a position angle of ~70.5° in the first epoch and with a similar position in the second epoch. This corresponds to a projected separation of ~9 au. The extracted photometry allowed us to estimate a mass for the companion of between 2 and 5 M Jup . This mass is in good agreement with astrometric measurements of the dynamic mass of the companion,…

The current discrepancy between the Hubble constant, H 0 , derived from the local distance ladder and from the cosmic microwave background is one of the most crucial issues in cosmology, as it may possibly indicate unknown systematics or new physics. Here, we present a novel non-parametric method to estimate the Hubble constant as a function of redshift. We establish independent estimates of the evolution of Hubble constant by diagonalizing the covariance matrix. From type Ia supernovae, baryon acoustic oscillation data and the observed Hubble parameter data, a decreasing trend in the Hubble constant with a significance of a 5.6 σ confidence level is found. At low redshift, its value is dramatically consistent with that measured from the local distance ladder and it drops to the value measured from the cosmic microwave background at high redshift. Our results may relieve the Hubble tension, with a preference for recent solutions, especially with respect to novel physics.

We explored the impact of the latest equation of state (EOS) for dense hydrogen–helium mixtures, which takes into account the interactions between hydrogen and helium species during the evolution of very low-mass stars and brown dwarfs (BDs). These interactions modify the thermodynamic properties of the H/He mixture, notably the entropy, a quantity of prime importance for these fully convective bodies, but also the onset and the development of degeneracy throughout the body. This translates into a faster cooling rate, that is, cooler isentropes for a given mass and age, and thus larger BD masses and smaller radii for a given effective temperature and luminosity than the models based on previous EOSs. This means that objects of a given mass and age in the rangeM ≲ 0.1 M⊙,τ ≳ 108yr will have cooler effective temperatures and fainter luminosities. Confronting these new models with several observationally determined BD dynamical masses, we show that this improves the agreement between evolutionary models and observations and resolves at least part of the observed discrepancy between the…

We present the first analysis of cosmic shear measured in DES Y3 that employs the entire range of angular scales in the data. To achieve this, we built upon recent advances in the theoretical modelling of weak lensing provided by a combination of N -body simulations, physical models of baryonic processes, and neural networks. Specifically, we used BACCOemu to model the linear and non-linear matter power spectrum including baryonic physics, allowing us to robustly exploit scales smaller than those used by the DES Collaboration. We show that the additional data produce cosmological parameters that are tighter but consistent with those obtained from larger scales, while also constraining the distribution of baryons. In particular, we measured the mass scale at which haloes have lost half of their gas, log M c = 14.38 −0.56 +0.60 log( h −1 M ⊙ ), and a parameter that quantifies the weighted amplitudes of the present-day matter inhomogeneities, S 8 = 0.799 −0.015 +0.023 . Our constraint on S 8 is statistically compatible with that inferred from the Planck satellite data at the 0.9 σ…

Context. Quasi-periodic eruptions (QPEs) are repeating thermal X-ray bursts associated with accreting massive black holes, the precise underlying physical mechanisms of which are still unclear. Aims. We present a new candidate QPE source, AT 2019vcb (nicknamed Tormund by the Zwicky Transient Facility Collaboration), which was found during an archival search for QPEs in the XMM-Newton archive. It was first discovered in 2019 as an optical tidal disruption event (TDE) at z = 0.088, and its X-ray follow-up exhibited QPE-like properties. Our goals are to verify its robustness as QPE candidate and to investigate its properties to improve our understanding of QPEs. Methods. We performed a detailed study of the X-ray spectral behaviour of this source over the course of the XMM-Newton archival observation. We also report on recent Swift and NICER follow-up observations to constrain the source’s current activity and overall lifetime, as well as an optical spectral follow-up. Results. The first two Swift detections and the first half of the 30 ks XMM-Newton exposure of Tormund displayed a…

Context. Expanding the sample of directly imaged companions to nearby, young stars that are amenable to detailed astrometric and spectroscopic studies is critical for the continued development and validation of theories of their evolution and atmospheric processes. Aims. The recent release of the Gaia astrometric catalog allows us to efficiently search for these elusive companions by targeting those stars that exhibit the astrometric reflex motion induced by an orbiting companion. The nearby (27 pc), young (24 Myr) star AF Leporis (AF Lep) was targeted because of its significant astrometric acceleration measured between the H IPPARCOS and Gaia astrometric catalogs, consistent with a wide-orbit planetary-mass companion detectable with high-contrast imaging. Methods. We used the SPHERE instrument on the VLT to search for faint substellar companions in the immediate vicinity of AF Lep. We used observations of a nearby star interleaved with those of AF Lep to efficiently subtract the residual point spread function. This provided sensitivity to faint planetary-mass companions within 1″…

Context. Debris from past merger events is expected and also known, to some extent, to populate the stellar halo near the Sun. Aims. We aim to identify and characterise such merger debris using Gaia DR3 data supplemented with metallicity and chemical abundance data from LAMOST LRS and APOGEE for halo stars within 2.5 kpc from the Sun. Methods. We utilised a single linkage-based clustering algorithm to identify over-densities in the integrals of motion space that could be due to merger debris. Combined with metallicity information and chemical abundances, we characterised these statistically significant over-densities. Results. We find that the local stellar halo contains seven main dynamical groups, with some of them shown to be in situ and some of accreted origin, most of which are already known. We report the discovery of a new substructure, which we dubbed ED-1. In addition, we find evidence for 11 independent smaller clumps, 5 of which are new: ED-2, 3, 4, 5, and 6, and typically rather tight dynamically. We identify their narrow range of metallicities, along with their…

Aims. We aim to derive a new constraint on the expansion history of the Universe by applying the cosmic chronometers method in the VANDELS survey, studying the age evolution of high-redshift galaxies with a full-spectral-fitting approach. Methods. We selected a sample of 39 massive (log( M ⋆ / M ⊙ ) > 10.8) and passive (log(sSFR/yr −1 ) < −11) galaxies from the fourth data release of the VANDELS survey at 1 < z < 1.5. To minimise the potential contamination by star-forming outliers, we selected our sample by combining different selection criteria, considering both photometric and spectroscopic information. The analysis of the observed spectral features provides direct evidence of an age evolution with redshift and of mass-downsizing, with more massive galaxies presenting stronger age-related features. To estimate the physical properties of the sample, we performed full spectral fitting with the code BAGPIPES , jointly analysing spectra and photometry of our sources without any cosmological assumption regarding the age of the population. Results. The derived physical properties of…

We report the discovery of a remarkable Lyα emitting galaxy at z = 7.2782, JADES-GS+53.16746−27.7720 (shortened to JADES-GS-z7-LA), with rest-frame equivalent width, EW0(Lyα) = 388.0 ± 88.8 Å and UV magnitude −17.0. The spectroscopic redshift is confirmed via rest-frame optical lines [O II], Hβ and [O III] in its JWST/NIRSpec Micro-Shutter Assembly (MSA) spectrum. The Lyα line is detected in both lower resolution (R ∼ 100) PRISM as well as medium resolution (R ∼ 1000) G140M grating spectra. The line spread function-deconvolved Lyα full width at half maximum in the grating is 383.9 ± 56.2 km s−1 and the Lyα velocity offset compared to the systemic redshift is 113.3 ± 80.0 km s−1, indicative of very little neutral gas or dust within the galaxy. We estimate the Lyα escape fraction to be > 70%. JADES-GS-z7-LA has a [O III]/[O II] ratio (O32) of 11.1 ± 2.2 and a ([O III] + [O II])/Hβ ratio (R23) of 11.2 ± 2.6, consistent with low metallicity and high ionization parameters. Deep NIRCam imaging also revealed a close companion source (separated by 0.23″), which exhibits similar photometry…

The distribution of close-in exoplanets is shaped by a complex interplay between atmospheric and dynamical processes. The Desert-Rim Exoplanets Atmosphere and Migration (DREAM) program aims at disentangling those processes through the study of the hot Neptune desert, whose rim hosts planets that are undergoing, or survived, atmospheric evaporation and orbital migration. In this first paper, we use the Rossiter-McLaughlin revolutions (RMR) technique to investigate the orbital architecture of 14 close-in planets ranging from mini-Neptune to Jupiter-size and covering a broad range of orbital distances. While no signal is detected for the two smallest planets, we were able to constrain the sky-projected spin-orbit angle of six planets for the first time, to revise its value for six others, and, thanks to constraints on the stellar inclination, to derive the 3D orbital architecture in seven systems. These results reveal a striking three-quarters of polar orbits in our sample, all being systems with a single close-in planet but of various stellar and planetary types. High-eccentricity…

We present a new catalogue of low-mass X-ray binaries (LMXBs) in the Galaxy. The catalogue contains source names, coordinates, source types, fluxes, distances, system parameters, and other characteristic properties of 349 LMXBs, including systems that have been newly discovered or reclassified since the most recently reported LMXB catalogues. The aim of this catalogue is to provide a list of all currently known Galactic objects identified as LMXBs with some basic information on each system (including X-ray and optical/IR properties where possible). Literature published before May 2023 has been taken into account where possible when compiling this information. References for all reported properties as well as object-finding charts in several energy bands are provided as part of the catalogue. We plan to update the catalogue regularly, in particular to reflect new objects discovered in the ongoing large-scale surveys such as Gaia and eROSITA.

We present JWST/Near Infrared Spectrograph (NIRSpec) integral field spectroscopy (IFS) of a lensed Population III candidate stellar complex (dubbed Lensed And Pristine 1, LAP1), with a lensing-corrected stellar mass of ≲10 4 M ⊙ and an absolute luminosity of M UV > −11.2 ( m UV > 35.6), confirmed at redshift 6.639 ± 0.004. The system is strongly amplified ( μ ≳ 100) by straddling a critical line of the Hubble Frontier Field galaxy cluster MACS J0416. Although the stellar continuum is currently not detected in the Hubble and JWST/Near Infrared Camera (NIRCam) and Near Infrared Imager and Slitless Spectrograph (NIRISS) imaging, arclet-like shapes of Lyman and Balmer lines, Ly α , H γ , H β and H α are detected with NIRSpec IFS with signal-to-noise ratios (S/N) of approximately 5 − 13 and large equivalent widths (> 300 − 2000 Å), along with a remarkably weak [O III ] λ λ 4959, 5007 at S / N ≃ 4. LAP1 shows a large ionizing photon production efficiency, log( ξ ion [erg Hz −1 ]) > 26. From the metallicity indexes R23 = ([O III ] + [O II ])/H β ≲ 0.74 and R3 = ([O III ]/H β ) = 0.55 ±…

Aims.Integral field spectroscopy (IFS) with JWST NIRSpec will significantly improve our understanding of the first quasars, by providing spatially resolved, infrared spectroscopic capabilities that cover key rest-frame optical emission lines that have been previously unobservable. Methods.Here we present our results from the first twoz > 6 quasars observed as a part of the Galaxy Assembly with NIRSpec IFS (GA-NIFS) GTO programme, with DELS J0411–0907 atz = 6.82 and VDES J0020–3653 atz = 6.86. Results.By observing the Hβ, [O III]λλ4959, 5007, and Hαemission lines in these high-zquasars for the first time, we measured accurate black hole masses,MBH = 1.85−0.8+2 × 109 M⊙and 2.9−1.3+3.5 × 109M⊙, corresponding to Eddington ratios ofλEdd = 0.8−0.4+0.7and 0.4−0.2+0.3for DELS J0411–0907 and VDES J0020–3653, respectively. These provide a key comparison for existing estimates from the more uncertain Mg IIline. We performed quasar–host decomposition using models of the quasars’ broad lines to measure the underlying host galaxies. We also discovered multiple emission line regions surrounding…

The second data release of the LOFAR Two-Metre Sky Survey (LoTSS) covers 27% of the northern sky, with a total area of ~5700 deg 1 . The high angular resolution of LOFAR with Dutch baselines (6 arcsec) allows us to carry out optical identifications of a large fraction of the detected radio sources without further radio followup; however, the process is made more challenging by the many extended radio sources found in LOFAR images as a result of its excellent sensitivity to extended structure. In this paper we present source associations and identifications for sources in the second data release based on optical and near-infrared data, using a combination of a likelihood-ratio cross-match method developed for our first data release, our citizen science project Radio Galaxy Zoo: LOFAR, and new approaches to algorithmic optical identification, together with extensive visual inspection by astronomers. We also present spectroscopic or photometric redshifts for a large fraction of the optical identifications. In total 4 116 934 radio sources lie in the area with good optical data, of…

Super-Eddington accretion is one scenario that may explain the rapid assembly of ∼109 M⊙ supermassive black holes (BHs) within the first billion year of the Universe. This critical regime is associated with radiatively inefficient accretion and accompanied by powerful outflows in the form of winds and jets. By means of hydrodynamical simulations of BH evolution in an isolated galaxy and its host halo with 12 pc resolution, we investigate how super-Eddington feedback affects the mass growth of the BH. It is shown that super-Eddington feedback efficiently prevents BH growth within a few Myr. The super-Eddington accretion events remain relatively mild with typical rates of about 2−3 times the Eddington limit, because of the efficient regulation by jets in that regime. We find that these jets are powerful enough to eject gas from the centre of the host galaxy all the way up to galactic scales at a few kpc, but do not significantly impact gas inflows at those large scales. By varying the jet feedback efficiency, we find that weaker super-Eddington jets allow for more significant BH…

Context. Isotope abundance ratios provide a powerful tool for tracing stellar nucleosynthesis, evaluating the composition of stellar ejecta, and constraining the chemical evolution of the Milky Way. Aims. We aim to measure the 12C/13C, 32S/34S, 32S/33S, 32S/36S, 34S/33S, 34S/36S, and 33S/36S isotope ratios across the Milky Way. Methods. With the IRAM 30 meter telescope, we performed observations of the J = 2−1 transitions of CS, C33S, C34S, C36S, 13CS, 13C33S, and 13C34S as well as the J = 3−2 transitions of C33S, C34S, C36S, and 13CS toward a large sample of 110 high-mass star-forming regions. Results. We measured the 12C/13C, 32S/34S, 32S/33S, 32S/36S, 34S/33S, 34S/36S, and 33S/36S abundance ratios with rare isotopologs of CS, thus avoiding significant saturation effects. With accurate distances obtained from parallax data, we confirm previously identified 12C/13C and 32S/34S gradients as a function of galactocentric distance. In the central molecular zone, 12C/13C ratios are higher than suggested by a linear fit to the disk values as a function of galactocentric radius. While…

The nucleosynthetic isotope dichotomy between carbonaceous (CC) and non-carbonaceous (NC) meteorites has been interpreted as evidence for spatial separation and the coexistence of two distinct planet-forming reservoirs for several million years in the solar protoplanetary disk. The rapid formation of Jupiter’s core within one million years after the formation of calcium-aluminium-rich inclusions (CAIs) has been suggested as a potential mechanism for spatial and temporal separation. In this scenario, Jupiter’s core would open a gap in the disk and trap inward-drifting dust grains in the pressure bump at the outer edge of the gap, separating the inner and outer disk materials from each other. We performed simulations of dust particles in a protoplanetary disk with a gap opened by an early-formed Jupiter core, including dust growth and fragmentation as well as dust transport, using the dust evolution software DustPy . Our numerical experiments indicate that particles trapped in the outer edge of the gap rapidly fragment and are transported through the gap, contaminating the inner disk…

We present a study of the environments of 17 Lyman- α emitting galaxies (LAEs) in the reionisation-era (5.8 < z < 8) identified by JWST/NIRSpec as part of the JWST Advanced Deep Extragalactic Survey (JADES). Unless situated in sufficiently (re)ionised regions, Lyman- α emission from these galaxies would be strongly absorbed by neutral gas in the intergalactic medium (IGM). We conservatively estimate sizes of the ionised regions required to reconcile the relatively low Lyman- α velocity offsets (Δ v Ly α < 300 km s −1 ) with moderately high Lyman- α escape fractions ( f esc, Ly α > 5%) observed in our sample of LAEs, suggesting the presence of ionised hydrogen along the line of sight towards at least eight out of 17 LAEs. We find minimum physical ‘bubble’ sizes of the order of R ion ∼ 0.1–1 pMpc are required in a patchy reionisation scenario where ionised bubbles containing the LAEs are embedded in a fully neutral IGM. Around half of the LAEs in our sample are found to coincide with large-scale galaxy overdensities seen in FRESCO at z ∼ 5.8–5.9 and z ∼ 7.3, suggesting Lyman- α…

Context. The Juno mission has provided measurements of Jupiter’s gravity field with an outstanding level of accuracy, leading to better constraints on the interior of the planet. Improving our knowledge of the internal structure of Jupiter is key to understanding its formation and evolution but is also important in the framework of exoplanet exploration. Aims. In this study, we investigated the differences between the state-of-the-art equations of state and their impact on the properties of interior models. Accounting for uncertainty on the hydrogen and helium equation of state, we assessed the span of the interior features of Jupiter. Methods. We carried out an extensive exploration of the parameter space and studied a wide range of interior models using Markov chain Monte Carlo simulations. To consider the uncertainty on the equation of state, we allowed for modifications of the equation of state in our calculations. Results. Our models harbour a dilute core and indicate that Jupiter’s internal entropy is higher than what is usually assumed from the Galileo probe measurements. We…

We report on observations of 68 satellites belonging to the SpaceX Starlink constellation with the LOFAR radio telescope. Radiation associated with Starlink satellites was detected at observing frequencies between 110 and 188 MHz, which is well below the 10.7– 12.7 GHz radio frequencies used for the downlink communication signals. A combination of broad-band features, covering the entire observed bandwidth, as well as narrow-band (bandwidth < 12.2 kHz) emission at frequencies of 125, 135, 143.05, 150, and 175 MHz, was observed. The presence and properties of both the narrow- and broad-band features vary between satellites at different orbital altitudes, indicating possible differences between the operational state of, or the hardware used in, these satellites. While the narrowband detections at 143.05 MHz can be attributed to reflections of radar signals from the French GRAVES Space Surveillance Radar, the signal properties of the broad- and narrow-band features at the other frequencies suggest that this radiation is intrinsic to the Starlink satellites and it is seen for 47 out of…

We present a new clustering method, significance mode analysis (SigMA), for extracting co-spatial and co-moving stellar populations from large-scale surveys such as ESAGaia. The method studies the topological properties of the density field in the multidimensional phase space. We validatedSigMAon simulated clusters and find that it outperforms competing methods, especially in cases where many clusters are closely spaced. We applied the new method toGaiaDR3 data of the closest OB association to Earth, Scorpio-Centaurus (Sco-Cen), and find more than 13 000 co-moving young objects, about 19% of which have a substellar mass.SigMAfinds 37 co-moving clusters in Sco-Cen. These clusters are independently validated by their narrow Hertzsprung-Russell diagram sequences and, to a certain extent, by their association with massive stars too bright forGaia, and are hence unknown toSigMA. We compared our results with similar recent work and find that theSigMAalgorithm recovers richer populations, is able to distinguish clusters with velocity differences down to about 0.5 km s−1, and reaches…

Context. The Mid-InfraRed Instrument (MIRI) Medium Resolution Spectrometer (MRS) on board the James Webb Space Telescope (JWST) allows us to probe the inner regions of protoplanetary disks, where the elevated temperatures result in an active chemistry and where the gas composition may dictate the composition of planets forming in this region. The disk around the classical T Tauri star Sz 98, which has an unusually large dust disk in the millimetre with a compact core, was observed with the MRS, and we examine its spectrum here. Aims. We aim to explain the observations and put the disk of Sz 98 in context with other disks, with a focus on the H2O emission through both its ro-vibrational and pure rotational emission. Furthermore, we compare our chemical findings with those obtained for the outer disk from Atacama Large Millimeter/submillimeter Array (ALMA) observations. Methods. In order to model the molecular features in the spectrum, the continuum was subtracted and local thermodynamic equilibrium (LTE) slab models were fitted. The spectrum was divided into different wavelength…

According to the Λ cold dark matter (ΛCDM) cosmology, present-day galaxies with stellar masses M ⋆ > 10 11 M ⊙ should contain a sizable fraction of dark matter within their stellar body. Models indicate that in massive early-type galaxies (ETGs) with M ⋆ ≈ 1.5 × 10 11 M ⊙ , dark matter should account for ∼15% of the dynamical mass within one effective radius (1 R e ) and for ∼60% within 5 R e . Most massive ETGs have been shaped through a two-phase process: the rapid growth of a compact core was followed by the accretion of an extended envelope through mergers. The exceedingly rare galaxies that have avoided the second phase, the so-called relic galaxies, are thought to be the frozen remains of the massive ETG population at z ≳ 2. The best relic galaxy candidate discovered to date is NGC 1277, in the Perseus cluster. We used deep integral field George and Cynthia Mitchel Spectrograph (GCMS) data to revisit NGC 1277 out to an unprecedented radius of 6 kpc (corresponding to 5 R e ). By using Jeans anisotropic modelling, we find a negligible dark matter fraction within 5 R e ( f DM (5…

In the past few years, we have undertaken an extensive investigation of star clusters and their stellar populations in the Large and Small Magellanic Clouds (LMC and SMC) based on archival images collected with the Hubble Space Telescope. We present photometry and astrometry of stars in 101 fields observed with the Wide Field Channel of the Advanced Camera for Surveys and the Ultraviolet and Visual Channel and the Near-Infrared Channel of Wide Field Camera 3. These fields comprise 113 star clusters. We provide differential-reddening maps for those clusters with significant reddening variations across the field of view. We illustrate various scientific outcomes that arise from the early inspection of the photometric catalogs. In particular, we provide new insights into the extended main-sequence turnoff (eMSTO) phenomenon: (i) We detected eMSTOs in two clusters, KMHK 361 and NGC 265, which had no previous evidence of multiple populations. This finding corroborates the conclusion that the eMSTO is a widespread phenomenon among clusters younger than ∼2 Gyr. (ii) The homogeneous…

Over the last few years, many studies have found an empirical relationship between the abundance of a star and its age. Here we estimate spectroscopic stellar ages for 178 825 red-giant stars observed by the APOGEE survey with a median statistical uncertainty of 17%. To this end, we use the supervised machine learning technique XGBoost , trained on a high-quality dataset of 3060 red-giant and red-clump stars with asteroseismic ages observed by both APOGEE and Kepler . After verifying the obtained age estimates with independent catalogues, we investigate some of the classical chemical, positional, and kinematic relationships of the stars as a function of their age. We find a very clear imprint of the outer-disc flare in the age maps and confirm the recently found split in the local age-metallicity relation. We present new and precise measurements of the Galactic radial metallicity gradient in small age bins between 0.5 and 12 Gyr, confirming a steeper metallicity gradient for ∼2 − 5 Gyr old populations and a subsequent flattening for older populations mostly produced by radial…

Context. The James Webb Space Telescope’s (JWST’s) unprecedented combination of sensitivity, spatial resolution, and infrared coverage has enabled a new era of galaxy morphology exploration across most of cosmic history. Aims. We analyze the near-infrared (NIR ∼ 0.8 − 1 μm) rest-frame morphologies of galaxies with log M * / M ⊙ > 9 in the redshift range of 0 < z < 6, compare with previous HST-based results and release the first JWST-based morphological catalog of ∼20 000 galaxies in the CEERS survey. Methods. We classified the galaxies in our sample into four main broad classes: spheroid, disk+spheroid, disk, and disturbed, based on imaging with four filters: F 150 W , F 200 W , F 356 W , and F 444 W . We used convolutional neural networks (CNNs) trained on HST/WFC3 labeled images and domain-adapted to JWST/NIRCam. Results. We find that ∼90% and ∼75% of galaxies at z < 3 have the same early and late and regular and irregular classification, respectively, in JWST and HST imaging when considering similar wavelengths. For small (large) and faint objects, JWST-based classifications tend…

Context. Recently, Doroshenko and collaborators reported a very low-mass compact star, a Central Compact Object named XMMU J173203.3−344518 inside the supernova remnant HESS J1731−347. Its tiny mass is at odds with all calculations of minimum masses of neutron stars generated by iron cores, therefore (and even if not compellingly) it has been suggested to be a “strange star”. In addition to the mass, the radius and surface temperature were extracted from the data, and the whole body of information should ultimately reveal whether this object is truly consistent with an exotic composition. Aims. Our aim is to understand the status of the compact object XMMU J173203.3−344518 in HESS J1731−347 within the existing models of strange stars, including its prompt formation. Methods. The information obtained on the mass, radius and surface temperature are compared to theoretical calculations performed within usual theoretical models using General Relativity as the assumed theory of gravitation and a handful of cooling scenarios. A qualitative discussion showing the consistency of the…

The redshifted 21 cm signal from neutral hydrogen is a direct probe of the physics of the early universe and has been an important science driver of many present and upcoming radio interferometers. In this study we use a single night of observations with the New Extension in Nançay Upgrading LOFAR (NenuFAR) to place upper limits on the 21 cm power spectrum from cosmic dawn at a redshift ofz = 20.3. NenuFAR is a new low-frequency radio interferometer, operating in the 10–85 MHz frequency range, currently under construction at the Nançay Radio Observatory in France. It is a phased array instrument with a very denseuvcoverage at short baselines, making it one of the most sensitive instruments for 21 cm cosmology analyses at these frequencies. Our analysis adopts the foreground subtraction approach, in which sky sources are modeled and subtracted through calibration and residual foregrounds are subsequently removed using Gaussian process regression. The final power spectra are constructed from the gridded residual data cubes in theuvplane. Signal injection tests are performed at each…

Observations of individual massive stars, super-luminous supernovae, gamma-ray bursts, and gravitational wave events involving spectacular black hole mergers indicate that the low-metallicity Universe is fundamentally different from our own Galaxy. Many transient phenomena will remain enigmatic until we achieve a firm understanding of the physics and evolution of massive stars at low metallicity (Z). The Hubble Space Telescope has devoted 500 orbits to observing ∼250 massive stars at low Z in the ultraviolet (UV) with the COS and STIS spectrographs under the ULLYSES programme. The complementary X-Shooting ULLYSES (XShootU) project provides an enhanced legacy value with high-quality optical and near-infrared spectra obtained with the wide-wavelength coverage X-shooter spectrograph at ESO’s Very Large Telescope. We present an overview of the XShootU project, showing that combining ULLYSES UV and XShootU optical spectra is critical for the uniform determination of stellar parameters such as effective temperature, surface gravity, luminosity, and abundances, as well as wind properties…

Context. Massive stars lose a large fraction of their mass to radiation-driven winds throughout their entire life. These outflows impact both the life and death of these stars and their surroundings. Aims. Theoretical mass-loss rates of hot, massive stars are derived to be used in applications such as stellar evolution. The behaviour of these rates in the OB-star regime is analysed, and their effects on massive-star evolution predictions studied. Methods. Dynamically consistent models are computed by solving the spherically symmetric, steady-state equation-of-motion for a large grid of hot, massive stars with different metallicities. The radiative acceleration is derived from non-local thermodynamic equilibrium radiative transfer in the co-moving frame, and all models cover a large spatial range from deep subsonic atmospheric layers into the radiation-driven and highly supersonic wind outflow. The resulting mass-loss rates are used to derive a simple scaling recipe with stellar parameters (luminosity, mass, effective temperature, and metallicity), and the new recipe is used to…

Context . The global context of making numerous data produced by researchers available requires collecting and organising the data, assigning meaningful metadata, and presenting the data in a meaningful and homogeneous way. The BASECOL database, which collects inelastic rate coefficients for application to the interstellar medium and to circumstellar and cometary atmospheres, meets those requirements. Aims . We aim to present the scientific content of the BASECOL2023 edition. Methods . While the previous versions relied on finding rate coefficients in the literature, the current version is populated with published results sent by the producers of data. The paper presents the database, the type of data that can be found, the type of metadata that are used, and the Virtual Atomic and Molecular Data Centre (VAMDC) standards that are used for the metadata. Finally, we present the different datasets species by species. Results . As the BASECOL database, interconnected with the VAMDC e-infrastructure, uses the VAMDC standards, the collisional data can be extracted with tools using VAMDC…

Context. We present a new catalogue of high mass X-ray binaries (HMXBs) in the Galaxy that we call the Galactic High Mass X-ray Binary Catalogue (XRBcats), which improves upon the most recent of such catalogues. We include new HMXBs discovered since previous publications and revise the classification for several objects previously considered HMXBs or HMXB candidates. The catalogue includes basic information (e.g. source names, coordinates, types), other data (e.g. distance and X-ray flux estimates, binary system parameters), and other characteristic properties of the 169 HMXBs catalogued. We also present finding charts in several bands from the infrared to hard X-rays for each object. Aims. The aim of this catalogue is to provide a list of all currently known Galactic HMXBs, including basic information on both compact objects and non-degenerate counterpart properties (where available). We also include objects tentatively classified as HXMBs in the literature and give a brief motivation for the classification in each case. Methods. The catalogue is compiled based on a search of known…

Despite their importance in stellar evolution, little is known about magnetic fields in the interior of stars. The recent seismic detection of magnetic fields in the core of several red giant stars has given measurements of their strength and information on their topology. We revisit the puzzling case of hydrogen-shell burning giants that show deviations from the expected regular period spacing of gravity modes. These stars also tend to have a measured period spacing that is too low compared to their counterparts. In this Letter, we show that these two features are well accounted for by strong magnetic fields in the cores of these stars. For 11 Kepler red giants showing these anomalies, we placed lower limits on the core field strengths ranging from 40 to 610 kG. For one star, the measured field exceeds the critical field above which gravity waves no longer propagate in the core. We find that this star shows mixed mode suppression at low frequency, which further suggests that this phenomenon might be related to strong core magnetic fields.

X-ray quasi-periodic eruptions (QPEs) are intense repeating soft X-ray bursts from the nuclei of nearby galaxies. Their physical origin is still largely unconstrained, and several theoretical models have been proposed ranging from disc instabilities to impacts between an orbiting companion and the existing accretion disc around the primary, or episodic mass transfer at pericentre in an extreme mass-ratio binary. We present here results from a recent XMM-Newton observation of GSN 069, the galactic nucleus where QPEs were first discovered. After about two years of absence, QPEs have reappeared in GSN 069, and we detect two consecutive QPEs separated by a much shorter recurrence time than ever before. Moreover, their intensites and peak temperatures are remarkably different, a novel addition to the QPE phenomenology. We study the QPE spectral properties from all XMM-Newton observations assuming QPEs to either represent an additional emission component superimposed on that from the disc, or the transient evolution of the disc emission itself. In the former scenario, QPEs are consistent…

Recent cosmic shear studies have shown that higher-order statistics (HOS) developed by independent teams now outperform standard two-point estimators in terms of statistical precision thanks to their sensitivity to the non-Gaussian features of large-scale structure. The aim of the Higher-Order Weak Lensing Statistics (HOWLS) project is to assess, compare, and combine the constraining power of ten different HOS on a common set ofEuclid-like mocks, derived fromN-body simulations. In this first paper of the HOWLS series, we computed the nontomographic (Ωm,σ8) Fisher information for the one-point probability distribution function, peak counts, Minkowski functionals, Betti numbers, persistent homology Betti numbers and heatmap, and scattering transform coefficients, and we compare them to the shear and convergence two-point correlation functions in the absence of any systematic bias. We also include forecasts for three implementations of higher-order moments, but these cannot be robustly interpreted as the Gaussian likelihood assumption breaks down for these statistics. Taken…

Context. Some accreting black holes exhibit much stronger variability patterns than the usual stochastic variations. Radiation pressure instability is one of the proposed mechanisms that might account for this effect. Aims. We model luminosity changes for objects with a black hole mass of 10, 10 5 , and 10 7 solar masses, using the time-dependent evolution of an accretion disk that is unstable as a result of the dominant radiation pressure. We concentrate on the outburst timescales. We explore the influence of the hot coronal flow above the cold disk, the inner purely hot flow, and the effect of the magnetic field on the time evolution of the disk-corona system. For intermediate-mass black holes and active galactic nuclei, we also explore the role of the disk outer radius because a disk that is fed by tidal disruption events (TDE) can be quite small. Methods. We used a 1D vertically integrated time-dependent numerical scheme that models the simultaneous evolution of the disk and corona, which is coupled by the vertical mass exchange. We parameterized the strength of the large-scale…

Context . High-resolution spectrographs fed by adaptive optics (AO) provide a unique opportunity to characterize directly imaged exoplanets. Observations with such instruments allow us to probe the atmospheric composition, spin rotation, and radial velocity of the planet, thereby helping to reveal information on its formation and migration history. The recent upgrade of the Cryogenic High-Resolution Infrared Echelle Spectrograph (CRIRES+) at the VLT makes it a highly suitable instrument for characterizing directly imaged exoplanets. Aims . In this work, we report on observations of β Pictoris b with CRIRES+ and use them to constrain the planets atmospheric properties and update the estimation of its spin rotation. Methods . The data were reduced using the open-source pycrires package. We subsequently forward-modeled the stellar, planetary, and systematic contribution to the data to detect molecules in the planet’s atmosphere. We also used atmospheric retrievals to provide new constraints on its atmosphere. Results . We confidently detected water and carbon monoxide in the atmosphere…

Ices are the main carriers of volatiles in protoplanetary disks and are crucial to our understanding of the protoplanetary disk chemistry that ultimately sets the organic composition of planets. The Director’s Discretionary-Early Release Science (DD-ERS) program Ice Age on the James Webb Space Telescope (JWST) follows the ice evolution through all stages of star and planet formation. JWST’s exquisite sensitivity and angular resolution uniquely enable detailed and spatially resolved inventories of ices in protoplanetary disks. JWST/NIRSpec observations of the edge-on Class II protoplanetary disk HH 48 NE reveal spatially resolved absorption features of the major ice components H2O, CO2, and CO, and multiple weaker signatures from less abundant ices NH3, OCN−, and OCS. Isotopologue 13CO2 ice has been detected for the first time in a protoplanetary disk. Since multiple complex light paths contribute to the observed flux, the ice absorption features are filled in by ice-free scattered light. This implies that observed optical depths should be interpreted as lower limits to the total ice…

Context. The z ∼ 0.1 type-2 QSO J1430+1339, known as the “ Teacup ”, is a complex galaxy showing a loop of ionised gas ∼10 kpc in diameter, co-spatial radio bubbles, a compact (∼1 kpc) jet, and outflow activity. Its closeness offers the opportunity to study in detail the intricate interplay between the central supermassive black hole (SMBH) and the material in and around the galaxy, both the interstellar medium (ISM) and circumgalactic medium (CGM). Aims. We characterise the spatially resolved properties and effects of the galactic ionised gas outflow and compare them with those of the radio jet and with theoretical predictions to infer its acceleration mechanism. Methods. We used VLT/MUSE optical integral field spectroscopic observations to obtain flux, kinematic, and excitation maps of the extended (up to ∼100 kpc) ionised gas and to characterise the properties of stellar populations. We built radial profiles of the outflow properties as a function of distance from the active nucleus, from kiloparsec up to tens of kiloparsec scales, at ∼1 kpc resolution. Results. We detect a…

The origin of the elevated C/O ratios in discs around late M dwarfs compared to discs around solar-type stars is not well understood. Here we endeavour to reproduce the observed differences in the disc C/O ratios as a function of stellar mass using a viscosity-driven disc evolution model and study the corresponding atmospheric composition of planets that grow inside the water-ice line in these discs. We carried out simulations using a coupled disc evolution and planet formation code that includes pebble drift and evaporation. We used a chemical partitioning model for the dust composition in the disc midplane. Inside the water-ice line, the disc’s C/O ratio initially decreases to sub-stellar due to the inward drift and evaporation of water-ice-rich pebbles before increasing again to super-stellar values due to the inward diffusion of carbon-rich vapour. We show that this process is more efficient for very low-mass stars compared to solar-type stars due to the closer-in ice lines and shorter disc viscous timescales. In high-viscosity discs, the transition from sub-stellar to…

Context. Constraining the turbulence level and pebble size in protoplanetary disks is an essential initial step in understanding the aerodynamic properties of pebbles, which are crucial for planet formation. Recent laboratory experiments have revealed that destructive collisions of icy dust particles may occur at much lower velocities than previously believed. These low fragmentation velocities push down the maximum grain size in collisional growth models. Aims. Motivated by the smooth radial distribution of pebble sizes inferred from ALMA/VLA multi-wavelength continuum analysis, we propose a concise model to explain this feature and aim to constrain the turbulence level at the midplane of protoplanetary disks. Methods. Our approach is built on the assumption that the fragmentation threshold is the primary barrier limiting pebble growth within pressure maxima. Consequently, the grain size at the ring location can provide direct insights into the turbulent velocity governing pebble collisions and, by extension, the turbulence level at the disk midplane. We validate this method using…

The mass-loss rate (MLR) is one of the most important parameters in astrophysics, because it impacts many areas of astronomy, such as ionizing radiation, wind feedback, star-formation rates, initial mass functions, stellar remnants, supernovae, and so on. However, the most important modes of mass loss are also the most uncertain, as the dominant physical mechanisms that lead to this phenomenon are stull largely unknown. Here we assemble the most complete and clean red supergiant (RSG) sample (2121 targets) so far in the Small Magellanic Cloud (SMC) with 53 different bands of data to study the MLR of RSGs. In order to match the observed spectral energy distributions (SEDs), we created a theoretical grid of 17 820 oxygen-rich models (“normal” and “dusty” grids are half-and-half) using the radiatively driven wind model of the DUSTY code, covering a wide range of dust parameters. We select the best model for each target by calculating the minimal modified chi-square and visual inspection. The resulting MLRs from DUSTY are converted to real MLRs based on the scaling relation, for which a…

Context. Quaoar is a classical trans-Neptunian object (TNO) with an area-equivalent diameter of 1100 km and an orbital semi-major axis of 43.3 astronomical units. Based on stellar occultations observed between 2018 and 2021, an inhomogeneous ring (Q1R, i.e., Quaoar’s first ring) has been detected around this body. Aims. A new stellar occultation by Quaoar was observed on August 9, 2022, with the aim of improving Quaoar’s shape models and the physical parameters of Q1R, while searching for additional material around the body. Methods. The occultation provided nine effective chords across Quaoar, pinning down its size, shape, and astrometric position. Large facilities, such as Gemini North and the Canada-France-Hawaii Telescope (CFHT), were used to obtain high acquisition rates and signal-to-noise ratios. The light curves were also used to characterize the Q1R ring (radial profiles and orbital elements). Results. Quaoar’s elliptical fit to the occultation chords yields the limb with an apparent semi-major axis of 579.5 ± 4.0 km, apparent oblateness of 0.12 ± 0.01, and area-equivalent…

We have conducted a search for open clusters in the vicinity of classical Galactic Cepheids based on high-quality astrometry from the third data release (DR3) of the ESA mission Gaia to improve the calibration of the Leavitt law (LL). Our approach requires no prior knowledge of existing clusters, allowing us to both detect new host clusters and cross-check previously reported associations. Our Gold sample consists of 34 Cepheids residing in 28 open clusters, including 27 fundamental mode and seven overtone Cepheids. Three new bona fide cluster Cepheids are reported (V0378 Cen, ST Tau, and GH Lup) and the host cluster identifications for three others (VW Cru, IQ Nor, and SX Vel) are corrected. The fraction of Cepheids occurring in open clusters within 2 kpc of the Sun is f CC,2 kpc = 0.088 −0.019 +0.029 . Nonvariable cluster members allow us to determine cluster parallaxes to ∼7 μas in the range 12.5 < G < 17 mag, where recent studies found that parallax corrections by Lindegren et al. (2021, A&A, 649, A4, L21) are accurate and require no further offset corrections. By comparing…

Context. Despite their vital importance to understanding galaxy evolution and our own Galactic ecosystem, our knowledge of the physical properties of the hot X-ray emitting phase of the Milky Way is still inadequate. However, sensitive SRG/eROSITA large area surveys are now providing us with the long-sought data needed to mend this state of affairs. Aims. Our aim is to constrain the properties of the Milky Way hot halo emission toward intermediate Galactic latitudes close to the Galactic anti-center. Methods. We analyzed the spectral properties of the integrated soft X-ray emission observed by eROSITA in the relatively deep eFEDS field. Results. We observe a flux of 12.6 and 5.1 × 10 −12 erg cm −2 s −1 deg −2 in the total (0.3–2) and soft (0.3–0.6 keV) band. We measure the temperature and metal (oxygen) abundance of the hot circumgalactic medium (CGM) to be within kT CGM = 0.153–0.178 keV and Z CGM = 0.052–0.072 Z ⊙ , depending on the contribution of solar wind charge exchange (SWCX). Slightly higher CGM abundances Z CGM = 0.05–0.10 Z ⊙ are possible, considering the uncertain…

Aims. The demographics of the production and escape of ionizing photons from UV-faint early galaxies is a key unknown that has hindered attempts to discover the primary drivers of reionization. With the advent of JWST, it is finally possible to observe the rest-frame optical nebular emission from individual sub- L * z > 3 galaxies to measure the production rate of ionizing photons, ξ ion . Methods. Here we study a sample of 370 z ∼ 3 − 7 galaxies spanning −23 < M UV < −15.5 (median M UV ≈ −18) with deep multiband HST and JWST/NIRCam photometry that covers the rest-UV to the optical from the GLASS and UNCOVER JWST surveys. Our sample includes 102 galaxies with Lyman-alpha emission detected in MUSE spectroscopy. We used H α fluxes inferred from NIRCam photometry to estimate the production rate of ionizing photons that do not escape these galaxies, ξ ion (1 − f esc ). Results. We find median log 10 ξ ion (1 − f esc ) = 25.33 ± 0.47, with a broad intrinsic scatter of 0.42 dex, which implies a broad range of galaxy properties and ages in our UV-faint sample. Galaxies detected with…

The first James Webb Space Telescope (JWST) data on the massive colliding cluster El Gordo allow for 23 known families of multiply lensed images to be confirmed and for eight new members of these families to be identified. Based on these families, which have been confirmed spectroscopically by MUSE, we derived an initial lens model. This model guided the identification of 37 additional families of multiply lensed galaxies, among which 28 are entirely new systems, and nine were previously known. The initial lens model determined geometric redshifts for the 37 new systems. The geometric redshifts agree reasonably well with spectroscopic or photometric redshifts when those are available. The geometric redshifts enable two additional models that include all 60 families of multiply lensed galaxies spanning a redshift range 2 < z < 6. The derived dark-matter distribution confirms the double-peak configuration of mass found by earlier work with the southern and northern clumps having similar masses. We confirm that El Gordo is the most massive known cluster at z > 0.8 and has an estimated…

Context. Spatially resolved studies of the kinematics of galaxies provide crucial insights into their assembly and evolution, enabling one to infer the properties of the dark matter halos, derive the impact of feedback on the interstellar medium (ISM), as well as measure and characterize the outflow motions. To date, most of the kinematic studies at z = 0.5 − 3.5 have been obtained using emission lines tracing the warm, ionized gas (e.g., H α , [OII], and [OIII]). However, whether these provide an exhaustive or only a partial view of the dynamics of galaxies and of the properties of the ISM is still debated. Complementary insights into the cold gas kinematics are therefore needed. Aims. We present the Archival Large Program to Advance Kinematic Analysis (ALPAKA), a project aimed at gathering high-resolution observations of CO and [CI] emission lines of star-forming galaxies at z = 0.5 − 3.5 from the Atacama Large Millimeter Array (ALMA) public archive. With ≈147 h of total integration time, ALPAKA assembles ∼0.25″ observations for 28 star-forming galaxies, which is the largest…

In order to investigate the impact of radio jets on the interstellar medium (ISM) of galaxies hosting active galactic nuclei (AGN), we present subarcsecond-resolution Atacama Large Millimeter/submillimeter Array (ALMA) CO(2-1) and CO(3-2) observations of the Teacup galaxy. This is a nearby ( D L = 388 Mpc) radio-quiet type-2 quasar (QSO2) with a compact radio jet ( P jet ≈ 10 43 erg s −1 ) that subtends a small angle from the molecular gas disc. Enhanced emission line widths perpendicular to the jet orientation have been reported for several nearby AGN for the ionised gas. For the molecular gas in the Teacup, not only do we find this enhancement in the velocity dispersion but also a higher brightness temperature ratio ( T 32 / T 21 ) perpendicular to the radio jet compared to the ratios found in the galaxy disc. Our results and the comparison with simulations suggest that the radio jet is compressing and accelerating the molecular gas, and driving a lateral outflow that shows enhanced velocity dispersion and higher gas excitation. These results provide further evidence that the…

Context.GaiaData Release 3 (GaiaDR3) time series data may contain spurious signals related to the time-dependent scan angle. Aims.We aim to explain the origin of scan-angle-dependent signals and how they can lead to spurious periods, provide statistics to identify them in the data, and suggest how to deal with them inGaiaDR3 data and in future releases. Methods.Using realGaia(DR3) data alongside numerical and analytical models, we visualise and explain the features observed in the data. Results.We demonstrated withGaia(DR3) data that source structure (multiplicity or extendedness) or pollution from close-by bright objects can cause biases in the image parameter determination from which photometric, astrometric, and (indirectly) radial velocity time series are derived. These biases are a function of the time-dependent scan direction of the instrument and thus can introduce scan-angle-dependent signals, which due to the scanning-law-induced sampling ofGaiacan result in specific spurious periodic signals. Numerical simulations in which a period search is performed onGaiatime series…

Context. Binary stars play a crucial role in our understanding of the formation and evolution of star clusters and their stellar populations. Aims. We use Gaia Data Release 3 to homogeneously analyze 78 Galactic open clusters and the unresolved binary systems they host, each composed of two main sequence (MS) stars. Methods. We first investigated the structural parameters of these clusters, such as the core radius and the central density, and determined the cluster mass function (MF) and total mass by interpolating the density profile of each cluster. We measured the fraction of binaries with a large mass ratio and the fraction of blue straggler stars (BSSs), and finally investigated possible connections between the populations of binary stars and BSSs with the main parameters of the host cluster. Results. Remarkably, we find that the MFs of 78 analyzed open clusters follow a similar trend and are well reproduced by two single power-law functions, with a change in slope around masses of 1 M ⊙ . The fraction of binary stars ranges from ∼15% to more than ∼60% without significant…

Context. The Tibet ASγ and LHAASO collaborations recently reported the observation of a γ-ray diffuse emission with energy up to the PeV level from the Galactic plane. Aims. We discuss the relevance of non-uniform cosmic-ray transport scenarios and the implications of these results for cosmic-ray physics. Methods. We used the DRAGON and HERMES codes to build high-resolution maps and spectral distributions of that emission for several representative models under the condition that they reproduce a wide set of local cosmic-ray data up to 100 PeV. Results. We show that the energy spectra measured by Tibet ASγ, LHAASO, ARGO-YBJ, and Fermi-LAT in several regions of interest in the sky can all be reasonably described in terms of the emission arising by the Galactic cosmic-ray “sea”. We also show that all our models are compatible with IceTop γ-ray upper limits. Conclusions. We compare the predictions of conventional and space-dependent transport models with those data sets. Although the Fermi-LAT, ARGO-YBJ, and LHAASO preliminary data slightly favor this scenario, due to the still large…

Context. The origin of the stellar initial mass function (IMF) and its relation with the core mass function (CMF) are actively debated issues with important implications in astrophysics. Recent observations in the W43 molecular complex of top-heavy CMFs, with an excess of high-mass cores compared to the canonical mass distribution, raise questions about our understanding of the star formation processes and their evolution in space and time. Aims. We aim to compare populations of protostellar and prestellar cores in three regions imaged in the ALMA-IMF Large Program. Methods. We created an homogeneous core catalogue in W43, combining a new core extraction in W43-MM1 with the catalogue of W43-MM2&MM3 presented in a previous work. Our detailed search for protostellar outflows enabled us to identify between 23 and 30 protostellar cores out of 127 cores in W43-MM1 and between 42 and 51 protostellar cores out of 205 cores in W43-MM2&MM3. Cores with neither outflows nor hot core emission are classified as prestellar candidates. Results. We found a similar fraction of cores which are…

Recent studies of high-redshift galaxies with James Webb Space Telescope (JWST), such as GN-z11 at z = 10.6, show unexpectedly significant amounts of nitrogen (N) in their spectra. As this phenomenology appears to extend to gravitionally lensed galaxies at Cosmic noon such as the Sunburst Arc at z = 2.37, as well as globular clusters overall, we suggest that the common ingredient among them are very massive stars (VMSs) with zero-age main sequence (ZAMS) masses in the range of 100–1000 M ⊙ . The He II in the Sunburst Arc might also be the result of the disproportionally large contribution of VMS to the total stellar contribution. We analyse the pros and cons of the previous suggestions, including classical Wolf–Rayet (cWR) stars and supermassive stars (SMSs), to conclude that only our VMS alternative ticks all the relevant boxes. We discuss the VMS mass-loss history via their peculiar vertical evolution in the HR diagram resulting from a self-regulatory effect of these wind-dominated VMSs and we estimate that the large amounts of N present in star-forming galaxies may indeed result…

We present bolometric luminosities, black hole masses, and Eddington ratios for 42 luminous quasars at z ≳ 6 using high signal-to-noise ratio VLT/X-shooter spectra, acquired as part of the enlarged ESO Large Programme XQR-30. In particular, we derived the bolometric luminosities from the rest-frame 3000 Å luminosities using a bolometric correction from the literature, as well as the black hole masses by modeling the spectral regions around the C IV 1549 Å and the Mg II 2798 Å emission lines, with scaling relations calibrated in the Local Universe. We find that the black hole masses derived from both emission lines are in the same range and the scatter of the measurements agrees with expectations from the scaling relations. The Mg II-derived masses are between ∼(0.8−12) ×109 M⊙ and the derived Eddington ratios are within ∼0.13−1.73, with a mean (median) of 0.84(0.72). By comparing the total sample of quasars at z > 5.8, from this work and from the literature, to a bolometric luminosity distribution-matched sample at z ∼ 1.5, we find that quasars at high redshift host slightly less…

Context.The characterization of exoplanets requires the reliable determination of the fundamental parameters of their host stars. Spectral fitting plays an important role in this process. For the majority of stellar parameters, matching synthetic spectra to the observations provides a robust and unique solution for the fundamental parameters, such as effective temperature, surface gravity, abundances, radial and rotational velocities, among others. Aims.Here, we present a new software package for fitting high-resolution stellar spectra that is easy to use, available for common platforms, and free from commercial licenses. We call it PySME. It is based on the proven Spectroscopy Made Easy package, later referred to as IDL SME or “original” SME. Methods.The IDL (Interactive Data Language) part of the original SME code has been rewritten in Python, but we kept the efficient C++ and FORTRAN code responsible for molecular-ionization equilibrium, opacities, and spectral synthesis. In the process we updated some components of the optimization procedure to offer more flexibility and better…

Context. The latest Gaia data release in July 2022, DR3, in addition to the refinement of the astrometric and photometric parameters from DR2, added a number of important data products to those available in earlier releases, including radial velocity data, information on stellar multiplicity, and XP spectra of a selected sample of stars. Gaia has proved to be an important search tool for white dwarf stars, which are readily identifiable from their absolute G magnitudes as low luminosity objects in the Hertzsprung–Russell (H–R) diagram. Each data release has yielded large catalogs of white dwarfs, containing several hundred thousand objects, far in excess of the numbers known from all previous surveys (∼40 000). While the normal Gaia photometry ( G , G BP , and G RP bands) and astrometry can be used to identify white dwarfs with high confidence, it is much more difficult to parameterize the stars and determine the white dwarf spectral type from this information alone. Observing all stars in these catalogs with follow-up spectroscopy and photometry is also a huge logistical challenge…

We present MEGARA integral field unit (IFU) observations of five local type-2 quasars (QSO2s, z ∼ 0.1) from the Quasar Feedback (QSOFEED) sample. These active galactic nuclei (AGN) have bolometric luminosities of 10 45.5 − 46 erg s −1 and stellar masses of ∼10 11 M ⊙ . The LR-V grating of MEGARA allows us to explore the kinematics of the ionized gas through the [O III] λ 5007 Å emission line. The nuclear spectra of the five QSO2s, extracted in a circular aperture of ∼1.2″ (∼2.2 kpc) in diameter, matching the resolution of these seeing-limited observations, show signatures of high velocity winds in the form of broad (full width at half maximum, 1300 ≤ FWHM ≤ 2240 km s −1 ) and blueshifted components. We found that four out of the five QSO2s present outflows that we can resolve with our seeing-limited data, and they have radii ranging from 3.1 to 12.6 kpc. In the case of the two QSO2s with extended radio emission, we found that it is well aligned with the outflows, suggesting that low-power jets might be compressing and accelerating the ionized gas in these radio-quiet QSO2s. In the…

The cosmological principle is fundamental to the standard cosmological model. It assumes that the Universe is homogeneous and isotropic on very large scales. As the basic assumption, it must stand the test of various observations. In this work, we investigated the properties of the Pantheon+ sample, including redshift distribution and position distribution, and we give its constraint on the flat ΛCDM model: Ωm = 0.36 ± 0.02 and H0 = 72.83 ± 0.23 km s−1 Mpc−1. Then, using the region fitting (RF) method, we mapped the all-sky distribution of cosmological parameters (Ωm and H0) and find that the distribution significantly deviates from isotropy. A local matter underdensity region exists toward (308.4°−48.7+47.6, −18.2°−28.8+21.1) as well as a preferred direction of the cosmic anisotropy (313.4°−18.2+19.6, −16.8°−10.7+11.1) in galactic coordinates. Similar directions may imply that local matter density might be responsible for the anisotropy of the accelerated expansion of the Universe. Results of statistical isotropy analyses including Isotropy and Isotropy with real-data positions…

The Gaia -ESO Survey is a public spectroscopic survey that targeted ≳10 5 stars covering all major components of the Milky Way from the end of 2011 to 2018, delivering its final public release in May 2022. Unlike other spectroscopic surveys, Gaia -ESO is the only survey that observed stars across all spectral types with dedicated, specialised analyses: from O ( T eff ~ 30 000–52 000 K) all the way to K-M (≳3500 K). The physics throughout these stellar regimes varies significantly, which has previously prohibited any detailed comparisons between stars of significantly different types. In the final data release (internal data release 6) of the Gaia -ESO Survey, we provide the final database containing a large number of products, such as radial velocities, stellar parameters and elemental abundances, rotational velocity, and also, for example, activity and accretion indicators in young stars and membership probability in star clusters for more than 114 000 stars. The spectral analysis is coordinated by a number of working groups (WGs) within the survey, each specialised in one or more…

Detecting planets in the early stages of formation is key to reconstructing the history and diversity of fully developed planetary systems. The aim of this work is to identify potential signatures from planet-disc interaction in the circumstellar discs around MWC 480, HD 163296, AS 209, IM Lup, and GM Aur, through the study of molecular lines observed as part of the ALMA large program MAPS. Extended and localised perturbations in velocity, line width, and intensity have been analysed jointly using the DISCMINER modelling framework, in three bright CO isotopologues, 12 CO, 13 CO, and C 18 O J = 2−1, to provide a comprehensive summary of the kinematic and column density substructures that planets might be actively sculpting in these discs. We find convincing evidence for the presence of four giant planets located at wide orbits in three of the discs in the sample: two around HD 163296, one in MWC 480, and one in AS 209. One of the planet candidates in HD 163296, P94, previously associated with velocity signatures detected in lower velocity resolution 12 CO data, is confirmed and…

Context. Retrieving the physical parameters from spectroscopic observations of exoplanets is key to understanding their atmospheric properties. Exoplanetary atmospheric retrievals are usually based on approximate Bayesian inference and rely on sampling-based approaches to compute parameter posterior distributions. Accurate or repeated retrievals, however, can result in very long computation times due to the sequential nature of sampling-based algorithms. Aims. We aim to amortize exoplanetary atmospheric retrieval using neural posterior estimation (NPE), a simulation-based inference algorithm based on variational inference and normalizing flows. In this way, we aim (i) to strongly reduce inference time, (ii) to scale inference to complex simulation models with many nuisance parameters or intractable likelihood functions, and (iii) to enable the statistical validation of the inference results. Methods. We evaluated NPE on a radiative transfer model for exoplanet spectra ( petitRADTRANS ), including the effects of scattering and clouds. We trained a neural autoregressive flow to…

We report on the discovery of the C7N− anion towards the starless core TMC-1 and towards the carbon-rich evolved star IRC +10216. We used the data of the QUIJOTE (Q-band Ultrasensitive Inspection Journey to the Obscure TMC-1 Environment) line survey towards TMC-1 and found six lines in perfect harmonic frequency relation from J = 27 − 26 up to J = 32 − 31. The frequency of the lines can be reproduced with a rotational constant and a distortion constant of B = 582.68490 ± 0.00024 MHz and D = 4.01 ± 0.13 Hz, respectively. The standard deviation of the fit is 4 kHz. Towards IRC +10216, we identify 17 lines from J = 27 − 26 up to J = 43 − 42; their frequencies are also in harmonic relation, providing B = 582.6827 ± 0.00085 MHz and D = 3.31 ± 0.31 Hz. The nearly exact coincidence of the rotational and distortion constants in both sources points unambiguously to a common molecular carrier. Taking into account the chemical peculiarities of both sources, the carrier could be a radical or an anion. The radical can be discarded, as the observed lines belong to a singlet species. Hence, the…

Context. The sample of Solar system objects has dramatically increased over the last decade. The number of measured properties (e.g., diameter, taxonomy, rotation period, thermal inertia, etc.) has expanded even more quickly. However, this wealth of information is spread over a myriad of studies, with different designations reported per object. Aims. We provide a solution to the identification of Solar system objects based on any of their multiple names or designations. We also compile and rationalize their properties to provide an easy access to them. We aim to continuously update the database as new measurements become available. Methods. We built a Web Service, SsODNet , which offers four access points, each corresponding to an identified necessity in the community: name resolution ( quaero ), compilation of a large corpus of properties ( dataCloud ), determination of the best estimate among compiled values ( ssoCard ), and a statistical description of the population ( ssoBFT ). Results. The SsODNet interfaces are fully operational and freely accessible to everyone. The name…

Context. The detection of quasi-periodic oscillations (QPOs) in the light curves of active galactic nuclei (AGNs) can provide insights into the physics of the super-massive black holes (SMBHs) powering these systems and could represent a signature of the existence of SMBH binaries, setting fundamental constraints on SMBH evolution in the Universe. Aims. The identification of long-term QPOs, characterized by periods on the order of several months to years, is particularly challenging and can only be achieved via all-sky monitoring instruments that provide unbiased, continuous light curves of astrophysical objects. The Fermi-LAT satellite, thanks to its monitoring observing strategy, is an ideal instrument for such a goal. Here, we aim to identify QPOs in the γ-ray light curves of the brightest AGNs within the Fermi-LAT catalog. Methods. We analyzed the light curves of the 35 brightest Fermi-LAT AGNs, including data from the beginning of the Fermi mission (August 2008) to April 2021, with energies from 100 MeV to 300 GeV. Two time binnings were investigated: 7 and 30 days. The search…

We present constraints on the flat Λ cold dark matter cosmological model through a joint analysis of galaxy abundance, galaxy clustering, and galaxy-galaxy lensing observables with the Kilo-Degree Survey. Our theoretical model combines a flexible conditional stellar mass function, which describes the galaxy-halo connection, with a cosmological N-body simulation-calibrated halo model, which describes the non-linear matter field. Our magnitude-limited bright galaxy sample combines nine-band optical-to-near-infrared photometry with an extensive and complete spectroscopic training sample to provide accurate redshift and stellar mass estimates. Our faint galaxy sample provides a background of accurately calibrated lensing measurements. We constrain the structure growth parameter to S8 = σ8√Ωm/0.3 =√0.773−0.030+0.028 and the matter density parameter to Ωm = 0.290−0.017+0.021. The galaxy-halo connection model adopted in the work is shown to be in agreement with previous studies. Our constraints on cosmological parameters are comparable to, and consistent with, joint ‘3 × 2pt’…

Context. The Medium Resolution Spectrometer (MRS) is one of the four observing modes of JWST/MIRI. Using JWST in-flight data of unresolved (point) sources, we can derive the MRS absolute spectral response function (ASRF) starting from raw data. Spectral fringing, caused by coherent reflections inside the detector arrays, plays a critical role in the derivation and interpretation of the MRS ASRF. The fringe corrections implemented in the current pipeline are not optimal for non-extended sources, and a high density of molecular features particularly inhibits an accurate correction. Aims. In this paper, we present an alternative way to calibrate the MIRI/MRS data. Firstly, we derive a fringe correction that accounts for the dependence of the fringe properties on the MIRI/MRS pupil illumination and detector pixel sampling of the point spread function. Secondly, we derive the MRS ASRF using an absolute flux calibrator observed across the full 5–28 µm wavelength range of the MRS. Thirdly, we apply the new ASRF to the spectrum of a G dwarf and compare it with the output of the JWST/MIRI…