Dating the ages and weighting the stellar populations in galaxies are essential steps when studying galaxy formation through cosmic times. Evolutionary population synthesis models with different input physics are used for this purpose. Moreover, the contribution from the thermally pulsing asymptotic giant branch (TP-AGB) stellar phase, which peaks for intermediate-age 0.6-2 Gyr, has been debated for decades. Here we report the detection of strong cool-star signatures in the rest-frame near-infrared spectra of three young (~1Gyr), massive (~10^10Msun) quiescent galaxies at large look-back time, z=1-2, using JWST/NIRSpec. The coexistence of oxygen- and carbon-type absorption features, spectral edges and features from rare species, such as vanadium and possibly zirconium, reveal a strong contribution from TP-AGB stars. Population synthesis models with a significant TP-AGB contribution reproduce the observations better than those with a weak TP-AGB, which are commonly used. These findings call for revisions of published stellar population fitting results, as they point to populations with lower masses and younger ages and have further implications for cosmic dust production and chemical enrichment. New generations of improved models are needed, informed by these and future observations.

Geospatial Foundation Models (GeoFMs) are transforming Earth Observation (EO), but evaluation lacks standardized protocols. GEO-Bench-2 addresses this with a comprehensive framework spanning classification, segmentation, regression, object detection, and instance segmentation across 19 permissively-licensed datasets. We introduce ''capability'' groups to rank models on datasets that share common characteristics (e.g., resolution, bands, temporality). This enables users to identify which models excel in each capability and determine which areas need improvement in future work. To support both fair comparison and methodological innovation, we define a prescriptive yet flexible evaluation protocol. This not only ensures consistency in benchmarking but also facilitates research into model adaptation strategies, a key and open challenge in advancing GeoFMs for downstream tasks. Our experiments show that no single model dominates across all tasks, confirming the specificity of the choices made during architecture design and pretraining. While models pretrained on natural images (ConvNext ImageNet, DINO V3) excel on high-resolution tasks, EO-specific models (TerraMind, Prithvi, and Clay) outperform them on multispectral applications such as agriculture and disaster response. These findings demonstrate that optimal model choice depends on task requirements, data modalities, and constraints. This shows that the goal of a single GeoFM model that performs well across all tasks remains open for future research. GEO-Bench-2 enables informed, reproducible GeoFM evaluation tailored to specific use cases. Code, data, and leaderboard for GEO-Bench-2 are publicly released under a permissive license.

The Euclid mission aims to measure the positions, shapes, and redshifts of over a billion galaxies to provide unprecedented constraints on the nature of dark matter and dark energy. Achieving this goal requires a continuous reassessment of the mission's scientific performance, particularly in terms of its ability to constrain cosmological parameters, as our understanding of how to model large-scale structure observables improves. In this study, we present the first scientific forecasts using CLOE (Cosmology Likelihood for Observables in Euclid), a dedicated Euclid cosmological pipeline developed to support this endeavour. Using advanced Bayesian inference techniques applied to synthetic Euclid-like data, we sample the posterior distribution of cosmological and nuisance parameters across a variety of cosmological models and Euclid primary probes: cosmic shear, angular photometric galaxy clustering, galaxy-galaxy lensing, and spectroscopic galaxy clustering. We validate the capability of CLOE to produce reliable cosmological forecasts, showcasing Euclid's potential to achieve a figure of merit for the dark energy parameters $w_0$ and $w_a$ exceeding 400 when combining all primary probes. Furthermore, we illustrate the behaviour of the posterior probability distribution of the parameters of interest given different priors and scale cuts. Finally, we emphasise the importance of addressing computational challenges, proposing further exploration of innovative data science techniques to efficiently navigate the Euclid high-dimensional parameter space in upcoming cosmological data releases.

ESASky is a science-driven discovery portal to explore the multi-wavelength sky and visualise and access multiple astronomical archive holdings. The tool is a web application that requires no prior knowledge of any of the missions involved and gives users world-wide simplified access to the highest-level science data products from multiple astronomical space-based astronomy missions plus a number of ESA source catalogues. The first public release of ESASky features interfaces for the visualisation of the sky in multiple wavelengths, the visualisation of query results summaries, and the visualisation of observations and catalogue sources for single and multiple targets. This paper describes these features within ESASky, developed to address use cases from the scientific community. The decisions regarding the visualisation of large amounts of data and the technologies used were made in order to maximise the responsiveness of the application and to keep the tool as useful and intuitive as possible.

We present a comparison of our results from ground-based observations of asteroid (21) Lutetia with imaging data acquired during the flyby of the asteroid by the ESA Rosetta mission. This flyby provided a unique opportunity to evaluate and calibrate our method of determination of size, 3-D shape, and spin of an asteroid from ground-based observations. We present our 3-D shape-modeling technique KOALA which is based on multi-dataset inversion. We compare the results we obtained with KOALA, prior to the flyby, on asteroid (21) Lutetia with the high-spatial resolution images of the asteroid taken with the OSIRIS camera on-board the ESA Rosetta spacecraft, during its encounter with Lutetia. The spin axis determined with KOALA was found to be accurate to within two degrees, while the KOALA diameter determinations were within 2% of the Rosetta-derived values. The 3-D shape of the KOALA model is also confirmed by the spectacular visual agreement between both 3-D shape models (KOALA pre- and OSIRIS post-flyby). We found a typical deviation of only 2 km at local scales between the profiles from KOALA predictions and OSIRIS images, resulting in a volume uncertainty provided by KOALA better than 10%. Radiometric techniques for the interpretation of thermal infrared data also benefit greatly from the KOALA shape model: the absolute size and geometric albedo can be derived with high accuracy, and thermal properties, for example the thermal inertia, can be determined unambiguously. We consider this to be a validation of the KOALA method. Because space exploration will remain limited to only a few objects, KOALA stands as a powerful technique to study a much larger set of small bodies using Earth-based observations.

Synthetic Aperture Radar (SAR) data enables large-scale surveillance of maritime vessels. However, near-real-time monitoring is currently constrained by the need to downlink all raw data, perform image focusing, and subsequently analyze it on the ground. On-board processing to generate higher-level products could reduce the data volume that needs to be downlinked, alleviating bandwidth constraints and minimizing latency. However, traditional image focusing and processing algorithms face challenges due to the satellite's limited memory, processing power, and computational resources. This work proposes and evaluates neural networks designed for real-time inference on unfocused SAR data acquired in Stripmap and Interferometric Wide (IW) modes captured with Sentinel-1. Our results demonstrate the feasibility of using one of our models for on-board processing and deployment on an FPGA. Additionally, by investigating a binary classification task between ships and windmills, we demonstrate that target classification is possible.

We present a catalog of 717 candidate galaxies at $z > 8$ selected from 125 square arcminutes of NIRCam imaging as part of the JWST Advanced Deep Extragalactic Survey (JADES). We combine the full JADES imaging dataset with data from the JEMS and FRESCO JWST surveys along with extremely deep existing observations from HST/ACS for a final filter set that includes fifteen JWST/NIRCam filters and five HST/ACS filters. The high-redshift galaxy candidates were selected from their estimated photometric redshifts calculated using a template fitting approach, followed by visual inspection from seven independent reviewers. We explore these candidates in detail, highlighting interesting resolved or extended sources, sources with very red long-wavelength slopes, and our highest redshift candidates, which extend to $z_{phot} = 18$. We also investigate potential contamination by stellar objects, and do not find strong evidence from SED fitting that these faint high-redshift galaxy candidates are low-mass stars. Over 93\% of the sources are newly identified from our deep JADES imaging, including 31 new galaxy candidates at $z_{phot} > 12$. Using 42 sources in our sample with measured spectroscopic redshifts from NIRSpec and FRESCO, we find excellent agreement to our photometric redshift estimates, with no catastrophic outliers and an average difference of $\langle \Delta z = z_{phot}- z_{spec} \rangle= 0.26$. These sources comprise one of the most robust samples for probing the early buildup of galaxies within the first few hundred million years of the Universe's history.

Context. In the current era of multi-wavelength and multi-messenger astronomy, international organisations are actively working on the definition of new standards for the publication of astronomical data, and substantial effort is devoted to make them available through public archives. Aims. We present a set of tools that allow user-friendly access and basic scientific analysis of observations in Hierarchical Progressive Survey (HiPS) format, and we use them to gauge the quality of representative skymaps at ultraviolet, optical, and infrared wavelengths. Methods. We apply a fully-automatic procedure to derive aperture photometry in 10 different bands for the 323 nearby galaxies in the Herschel Reference Sample (HRS), and compare its results with the rigorous analyses involving specialised knowledge and human intervention carried out by the HRS team. Results. Our experiment shows that 9 of the 10 skymaps considered preserve the original quality of the data, and the photometric fluxes returned by our pipeline are consistent with the HRS measurements within a few per cent. In the case of Herschel PACS maps at 100 {\mu}m, we uncovered a systematic error that we ascribe to an inconsistent combination of data products with different spatial resolution. For the remaining skymaps, the estimated statistical uncertainties provide a realistic indication of the differences with respect to the HRS catalogue. Conclusions. In principle, the currently available HiPS skymaps in Flexible Image Transport System (FITS) format allow to carry out broadband photometric analyses with an accuracy of the order of a few percent, but some level human intervention is still required. In addition to assessing data quality, we also propose a series of recommendations to realise the full potential of the HiPS format for the scientific analysis of large astronomical data sets.

We discuss the methods used to compile a high signal-to-noise dataset representative of both the instrumental and cosmic background signal measured at high galactic latitude by the XMM-Newton EPIC cameras. The characteristics of the EPIC background are described and the potential applications of the derived dataset in general science analysis are outlined. In the case of the cosmic X-ray background, the transition between a hard power-law spectrum (due to the integrated emission of unresolved, largely extragalactic, point sources) and a softer thermal spectrum (produced by hot plasma associated with the Galactic plane and halo) is unambiguously detected around ~1keV. We derive a value for the intensity of the power-law component of 2.15 (+/- 0.26) e-11 erg/sq cm/s/sq deg in the 2-10 keV band (Normalisation at 1keV of 11.1 photons /sq cm/s/sr/keV). The implication is that recent, very deep Chandra observations have resolved ~70 - 90% of the 2-10 keV background into discrete sources. Our measurement is towards the higher end of the range of quoted background normalisations.

We report here on observations of a tidal disruption event, XMMSL2 J1404-2511, discovered in an XMM-Newton slew, in a quiescent galaxy at z=0.043. X-ray monitoring covered the epoch when the accretion disc transitioned from a thermal state, with kT~80 eV, to a harder state dominated by a warm, optically-thick corona. The bulk of the coronal formation took place within 7 days and was coincident with a temporary drop in the emitted radiation by a factor 4. After a plateau phase of ~100 days, the X-ray flux of XMMSL2 J1404-2511 decayed by a factor 500 within 230 days. We estimate the black hole mass in the galaxy to be $M_{BH}=4\pm{2}\times10^{6}$ solar masses and the peak X-ray luminosity $L_{X}\sim6\times10^{43}$ ergs/s. The optical/UV light curve is flat over the timescale of the observations with $L_{opt}\sim 2\times10^{41}$ ergs/s. We find that TDEs with coronae are more often found in an X-ray sample than in an optically-selected sample. Late-time monitoring of the optical sample is needed to test whether this is an intrinsic property of TDEs or is due to a selection effect. From the fast decay of the X-ray emission we consider that the event was likely due to the partial stripping of an evolved star rather than a full stellar disruption, an idea supported by the detection of two further re-brightening episodes, two and four years after the first event, in the SRG/eROSITA all-sky survey.

We investigate the composition of the solid-state materials in the winds around S-type AGB stars. The S stars produce dust in their wind that bears a resemblance to the dust produced in some O-rich AGB stars. However, the reported resemblance is mostly based on IRAS/LRS spectra with limited spectral resolution, sensitivity, and wavelength coverage. We investigate the dust composition around S stars using ISO/SWS data that surpass the previous studies in terms of spectral resolution and wavelength coverage. We compare the dust spectra from the 9 sources with the O-rich AGB spectra and a subset of M super-giants. We constructed average dust emission spectra of the different categories. We report the discovery of several previously unreported dust emission features in the S star spectra. The long wavelength spectra of W Aql and pi1 Gru exhibit the "30" micrometer feature attributed to MgS. Two sources exhibit a series of emission bands between 20 and 40 micrometer that we tentatively ascribe to Diopside. We show that the 10-20 micrometer spectra of the S stars are significantly different from the O-rich AGB stars. The O-rich stars exhibit a structured emission feature that is believed to arise from amorphous silicate and aluminium-oxide. The S stars lack the substructure found in the O-rich stars. Instead they show a smooth peak with a varying peak-position from source to source. We suggest that this feature is caused by a family of related material, whose exact composition determines the peak position. The observed trend mimics the laboratory trend of non-stoichiometric silicates. In this scenario the degree of non-stoichiometry is related to the Mg to SiO4 ratio, in other words, to the amount of free O available during the dust grain growth.

In meteor science, the identification of meteor showers is a crucial and complex problem. The most common method is to perform a systematic search of a database of observed orbits using an orbit dissimilarity criterion (D-criterion) and an algorithm. D-criteria compare the result of an orbit dissimilarity function (D-function) and a threshold. These D-functions associate one value to two meteoroids orbits. If this value is lower than the threshold, the meteoroids' orbits are considered similar. Group of meteors are thus formed. However, not all D-criteria have been evaluated, and their high number makes it hard to know which should be prioritised. This paper presents a review of each D-function, the tests they passed, the threshold choice, and the algorithms they are used with. We show what methods are currently used in the search for meteor showers, presenting statistics based on papers justifying the existence of established meteor showers. We review D-functions from eight different papers. We describe how thresholds and clustering algorithms are usually chosen. We also analyse tests that were performed on D-criteria. We discover that most of those criteria were not properly tested, and that some have been criticised for their theoretical background. Thus, we recommend performing a post-search analysis of the groups found, both in a statistical sense (to make sure the groups formed could not have been formed randomly) and an orbital dynamics sense (to check whether the group could indeed come from a singular parent body), to present the findings as potential meteor showers.

Observational astronomy of tidal disruption events (TDEs) began with the detection of X-ray flares from quiescent galaxies during the ROSAT all-sky survey of 1990-1991. The flares complied with theoretical expectations, having high peak luminosities ($L_{\rm x}$ up to $\ge4\times 10^{44}$ erg/s), a thermal spectrum with $kT\sim$few$\times10^5$ K, and a decline on timescales of months to years, consistent with a diminishing return of stellar debris to a black hole of mass $10^{6-8}$ solar masses. These measurements gave solid proof that the nuclei of quiescent galaxies are habitually populated by a super-massive black hole. Beginning in 2000, XMM-Newton, Chandra and Swift have discovered further TDEs which have been monitored closely at multiple wavelengths. A general picture has emerged of, initially near-Eddington accretion, powering outflows of highly-ionised material, giving way to a calmer sub-Eddington phase, where the flux decays monotonically, and finally a low accretion rate phase with a harder X-ray spectrum indicative of the formation of a disk corona. There are exceptions to this rule though which at the moment are not well understood. A few bright X-ray TDEs have been discovered in optical surveys but in general X-ray TDEs show little excess emission in the optical band, at least at times coincident with the X-ray flare. X-ray TDEs are powerful new probes of accretion physics down to the last stable orbit, revealing the conditions necessary for launching jets and winds. Finally we see that evidence is mounting for nuclear and non-nuclear intermediate mass black holes based on TDE flares which are relatively hot and/or fast.

Quasars powered by massive black holes (BHs) with mass estimates above a billion solar masses have been identified at redshift 6 and beyond. The existence of such BHs requires almost continuous growth at the Eddington limit for their whole lifetime, of order of one billion years. In this paper, we explore the possibility that positively skewed scale-dependent non-Gaussian primordial fluctuations may ease the assembly of massive BHs. In particular, they produce more low-mass halos at high redshift, thus altering the production of metals and ultra-violet flux, believed to be important factors in BH formation. Additionally, a higher number of progenitors and of nearly equal-mass halo mergers would boost the mass increase provided by BH-BH mergers and merger-driven accretion. We use a set of two cosmological simulations, with either Gaussian or scale-dependent non-Gaussian primordial fluctuations to perform a proof-of-concept experiment to estimate how BH formation and growth are altered. We estimate the BH number density and the fraction of halos where BHs form, for both simulations and for two popular scenarios of BH formation (remnants of the first generation of stars and direct collapse in the absence of metals and molecular hydrogen). We find that the fractions of halos where BHs form are almost identical, but that non-Gaussian primordial perturbations increase the total number density of BHs for the both BH formation scenarios. We also evolve BHs using merger trees extracted from the simulations and find that non-Gaussianities increase both the BH masses and the number of the most massive BHs.

Observations of planetary material polluting the atmospheres of white dwarfs are an important probe of the bulk composition of exoplanetary material. Medium- and high-resolution optical and ultraviolet spectroscopy of seven white dwarfs with known circumstellar dust and gas emission are presented. Detections or meaningful upper limits for photospheric absorption lines are measured for: C, O, Na, S, P, Mg, Al, Si, Ca, Ti, Cr, Fe, and Ni. For 16 white dwarfs with known observable gaseous emission discs (and measured photospheric abundances), there is no evidence that their accretion rates differ, on average, from those without detectable gaseous emission. This suggests that, typically, accretion is not enhanced by gas drag. At the effective temperature range of the white dwarfs in this sample (16,000-25,000K) the abundance ratios of elements are more consistent than absolute abundances when comparing abundances derived from spectroscopic white dwarf parameters versus photometric white dwarf parameters. Crucially, this highlights that the uncertainties on white dwarf parameters do not prevent white dwarfs from being utilised to study planetary composition. The abundances of oxygen and silicon for the three hydrogen-dominated white dwarfs in the sample with both optical and ultraviolet spectra differ by 0.62 dex depending on if they are derived from the optical or ultraviolet spectra. This optical/ultraviolet discrepancy may be related to differences in the atmospheric depth of line formation; further investigations into the white dwarf atmospheric modelling are needed to understand this discrepancy.

We present results of ultra-deep ISOCAM observations through a cluster-lens at 7 and 15 micron with the Infrared Space Observatory (ISO) satellite. These observations reveal a large number of luminous Mid-Infrared (MIR) sources. Cross-identification in the optical and Near-Infrared (NIR) wavebands shows that about half of the 7 micron sources are cluster galaxies. The other 7 micron and almost all 15 micron sources are identified as lensed distant galaxies. Thanks to the gravitational amplification they constitute the faintest MIR detected sources, allowing us to extend the number counts in both the 7 and 15 micron bands. In particular, we find that the 15 micron counts have a steep slope alpha_15 = -1.5 +/- 0.3 and are large, with N_15 (>30 microJy}) = 13 +/- 5 per square arcmin. These numbers rule out non-evolutionary models and favour very strong evolution. Down to our counts limit, we found that the resolved 7 and 15 microns background radiation intensity is respectively (2 +/-0.5) 10^(-9) and (5 +/-1) 10^(-9) W m^(-2) sr^(-1).

We have performed a comprehensive spectral and timing analysis of the first NuSTAR observation of the high-mass X-ray binary 4U 1538-522. The observation covers the X-ray eclipse of the source, plus the eclipse ingress and egress. We use the new measurement of the mid-eclipse time to update the orbital parameters of the system and find marginally-significant evolution in the orbital period, with $\dot{P}_{\rm orb}/P_{\rm orb} = \left(-0.95 \pm 0.37\right) \times 10^{-6}$ yr$^{-1}$. The cyclotron line energy is found approximately 1.2 keV higher than RXTE measurements from 1997--2003, in line with the increased energy observed by Suzaku in 2012 and strengthening the case for secular evolution of 4U 1538-522's CRSF. We additionally characterize the behavior of the iron fluorescence and emission lines and line-of-sight absorption as the source moves into and out of eclipse.

Nongravitational accelerations in the absence of observed activity have recently been identified on NEOs, opening the question of the prevalence of anisotropic mass-loss in the near-Earth environment. Motivated by the necessity of nongravitational accelerations to identify 2010 VL$_{65}$ and 2021 UA$_{12}$ as a single object, we investigate the problem of linking separate apparitions in the presence of nongravitational perturbations. We find that nongravitational accelerations on the order of $10^{-9}$ au/d$^2$ can lead to a change in plane-of-sky positions of $\sim10^3$ arcsec between apparitions. Moreover, we inject synthetic tracklets of hypothetical nongravitationally-accelerating NEOs into the Minor Planet Center orbit identification algorithms. We find that at large nongravitational accelerations ($|A_i|\geq10^{-8}$ au/d$^2$) these algorithms fail to link a significant fraction of these tracklets. We further show that if orbits can be determined for both apparitions, the tracklets will be linked regardless of nongravitational accelerations, although they may be linked to multiple objects. In order to aid in the identification and linkage of nongravitationally accelerating objects, we propose and test a new methodology to search for unlinked pairs. When applied to the current census of NEOs, we recover the previously identified case but identify no new linkages. We conclude that current linking algorithms are generally robust to nongravitational accelerations, but objects with large nongravitational accelerations may potentially be missed. While current algorithms are well-positioned for the anticipated increase in the census population from future survey missions, it may be possible to find objects with large nongravitational accelerations hidden in isolated tracklet pairs.