We introduce a novel approach to estimate the spectral index, $\beta_s$, of polarised synchrotron emission, combining the moment expansion of CMB foregrounds and the constrained-ILC method. We reconstruct the maps of the first two synchrotron moments, combining multi-frequency data, and apply the `T-T plot' technique between two moment maps to estimate the synchrotron spectral index. This approach offers a new technique for mapping the foreground spectral parameters, complementing the model-based parametric component separation methods. Applying this technique, we derive a new constraint on the spectral index of polarised synchrotron emission using QUIJOTE MFI wide-survey 11 and 13 GHz data, Wilkinson Microwave Anisotropy Probe (WMAP) data at K and Ka bands, and Planck LFI 30 GHz data. In the Galactic plane and North Polar Spur regions, we obtain an inverse-variance-weighted mean synchrotron index of $\beta_s = -3.11$ with a standard deviation of $0.21$ due to intrinsic scatter, consistent with previous results based on parametric methods using the same dataset. We find that the inverse-variance-weighted mean spectral index, including both statistical and systematic uncertainties, is $\beta_s^{\rm plane} = -3.05 \pm 0.01$ in the Galactic plane and $\beta_s^{\rm high\text{-}lat} = -3.13 \pm 0.02$ at high latitudes, indicating a moderate steepening of the spectral index from low to high Galactic latitudes. Our analysis indicates that, within the current upper limit on the AME polarisation fraction, our results are not subject to any appreciable bias. Furthermore, we infer the spectral index over the entire QUIJOTE survey region, partitioning the sky into 21 patches. This technique can be further extended to constrain the synchrotron spectral curvature by reconstructing higher-order moments when better-quality data become available.

The ethical and legal imperative to share research data without causing harm requires careful attention to privacy risks. While mounting evidence demonstrates that data sharing benefits science, legitimate concerns persist regarding the potential leakage of personal information that could lead to reidentification and subsequent harm. We reviewed metadata accompanying neuroimaging datasets from six heterogeneous studies openly available on OpenNeuro, involving participants across the lifespan, from children to older adults, with and without clinical diagnoses, and including associated clinical score data. Using metaprivBIDS (this https URL), a novel tool for the systematic assessment of privacy in tabular data, we found that privacy is generally well maintained, with serious vulnerabilities being rare. Nonetheless, minor issues were identified in nearly all datasets and warrant mitigation. Notably, clinical score data (e.g., neuropsychological results) posed minimal reidentification risk, whereas demographic variables (age, sex, race, income, and geolocation) represented the principal privacy vulnerabilities. We outline practical measures to address these risks, enabling safer data sharing practices.

Three black-hole low-mass X-ray binaries (LMXBs) in the Milky Way show rates of period decay which cannot be easily explained by standard mechanisms. Recently, it has been claimed that the anomalous period decays in two of these systems may be explained by dynamical friction due to very high dark matter (DM) densities around the black holes. We critically assess these claims by performing $N$-body simulations of binaries embedded in dense DM ``spikes". We simulate the previously-studied systems XTE J1118+480 and A0620--00, as well as studying the third binary Nova Muscae 1991 for the first time in this context. These simulations show that feedback on the DM distribution plays a crucial role and we rule out previously-claimed shallow DM spikes. We set lower limits on the steepness $\gamma$ of DM density profiles required to explain the period decay in these LMXBs, requiring $\gamma \gtrsim 2.15-2.20$ in XTE J1118+480 and A0620--00 and $\gamma \gtrsim 2.3$ in Nova Muscae 1991. Improved modeling of the long-term evolution of binaries embedded in DM spikes may allow us to exclude even larger densities in future.

We investigate the capability of the J-PAS survey to constrain the primordial power spectrum using a non-parametric Bayesian method. Specifically, we analyze simulated power spectra generated by a local oscillatory primordial feature template motivated by non-standard inflation. The feature is placed within the range of scales where the signal-to-noise ratio is maximized, and we restrict the analysis to $k \in [0.02,0.2] \text{ h} \text{ Mpc}^{-1}$, set by the expected J-PAS coverage and the onset of non-linear effects. Each primordial power spectrum is reconstructed by linearly interpolating $N$ knots in the $\{\log k, \log P_{\mathcal{R}}(k)\}$ plane, which are sampled jointly with the cosmological parameters $\{H_0,\Omega_b h^2, \Omega_c h^2\}$ using PolyChord. To test the primordial features, we apply two statistical tools: the Bayes factor and a hypothesis test that localizes the scales where features are detected. We assess the recovery under different J-PAS specifications, including redshift binning, tracer type, survey area, and filter strategy. Our results show that combining redshift bins and tracers allows the detection of oscillatory features as small as 2\%.

This is the second paper in the HOWLS (higher-order weak lensing statistics) series exploring the usage of non-Gaussian statistics for cosmology inference within \textit{Euclid}. With respect to our first paper, we develop a full tomographic analysis based on realistic photometric redshifts which allows us to derive Fisher forecasts in the ($\sigma_8$, $w_0$) plane for a \textit{Euclid}-like data release 1 (DR1) setup. We find that the 5 higher-order statistics (HOSs) that satisfy the Gaussian likelihood assumption of the Fisher formalism (1-point probability distribution function, $\ell$1-norm, peak counts, Minkowski functionals, and Betti numbers) each outperform the shear 2-point correlation functions by a factor $2.5$ on the $w_0$ forecasts, with only marginal improvement when used in combination with 2-point estimators, suggesting that every HOS is able to retrieve both the non-Gaussian and Gaussian information of the matter density field. The similar performance of the different estimators\inlinecomment{, with a slight preference for Minkowski functionals and 1-point probability distribution function,} is explained by a homogeneous use of multi-scale and tomographic information, optimized to lower computational costs. These results hold for the $3$ mass mapping techniques of the \textit{Euclid} pipeline: aperture mass, Kaiser--Squires, and Kaiser--Squires plus, and are unaffected by the application of realistic star masks. Finally, we explore the use of HOSs with the Bernardeau--Nishimichi--Taruya (BNT) nulling scheme approach, finding promising results towards applying physical scale cuts to HOSs.

We present new ALMA [OIII]$_{88}$ observations of eight previously [CII]$_{158}$-detected galaxies at $6.8 \lesssim z \lesssim 7.7$. Six of our targets -- the primary sample -- are massive, UV-luminous galaxies drawn from the REBELS survey, while the remaining two are UV-fainter galaxies that were previously serendipitously detected through their luminous [CII] lines in the REBELS fields. We detect [OIII]$_{88}$ emission in all eight galaxies at $6.2 - 17.7\sigma$ significance, and find them to be consistent with the local dwarf galaxy relation between $L_\mathrm{[OIII]}$ and star formation rate. Our sample spans [OIII]/[CII] $\approx 1.9 - 9.6$, which is typical for the high-redshift galaxy population. Five of the primary targets benefit from JWST/NIRSpec observations, enabling a direct comparison of the [OIII]/[CII] ratio against rest-optical ISM diagnostics. We supplement our high-redshift sample with eleven $z\approx6-14$ galaxies in the literature for which similar ALMA and JWST observations are available, and furthermore compare to the [OIII]/[CII] ratios measured for local dwarf galaxies. We find that, at fixed metallicity and ionization parameter, $z>6$ galaxies show elevated [OIII]/[CII] ratios compared to local dwarfs. Instead, we find that a large [OIII]$_{4959,5007}$+H$\beta$ equivalent width -- a proxy for burstiness -- is the main driver of the high [OIII]/[CII] ratios seen in the early Universe, which is primarily due to [CII] being suppressed in bursty galaxies. Given the apparent validity of the [OIII]$_{88}$-SFR relation across most of cosmic time, as well as the abundance of young, bursty galaxies at high redshift, [OIII]$_{88}$ is set to remain a powerful ISM tracer at the cosmic dawn.

Future gravitational wave observatories can probe dark matter by detecting the dephasing in the waveform of binary black hole mergers induced by dark matter overdensities. Such a detection hinges on the accurate modelling of the dynamical friction, induced by dark matter on the secondary compact object in intermediate and extreme mass ratio inspirals. In this paper, we introduce NbodyIMRI, a new publicly available code designed for simulating binary systems within cold dark matter `spikes'. Leveraging higher particle counts and finer timesteps, we validate the applicability of the standard dynamical friction formalism and provide an accurate determination of the maximum impact parameter of particles which can effectively scatter with a compact object, across various mass ratios. We also show that in addition to feedback due to dynamical friction, the dark matter also evolves through a `stirring' effect driven by the time-dependent potential of the binary. We introduce a simple semi-analytical scheme to account for this effect and demonstrate that including stirring tends to slow the rate of dark matter depletion and therefore enhances the impact of dark matter on the dynamics of the binary.

Future gravitational wave interferometers such as LISA, Taiji, DECIGO, and TianQin, will enable precision studies of the environment surrounding black holes. In this paper, we study intermediate and extreme mass ratio binary black hole inspirals, and consider three possible environments surrounding the primary black hole: accretion disks, dark matter spikes, and clouds of ultra-light scalar fields, also known as gravitational atoms. We present a Bayesian analysis of the detectability and measurability of these three environments. Focusing for concreteness on the case of a detection with LISA, we show that the characteristic imprint they leave on the gravitational waveform would allow us to identify the environment that generated the signal, and to accurately reconstruct its model parameters.

Neutral gas in galaxies during the Epoch of Reionisation regulates star formation, dust growth, and the escape of ionising photons, making it a key ingredient in understanding both galaxy assembly and reionisation. Yet, direct constraints on the HI content of galaxies at z>6 have been scarce. With JWST, Ly$\alpha$ damping wings in galaxy spectra can now provide a direct probe of this neutral component. We analyse JWST/NIRSpec prism spectra of 12 UV-luminous galaxies from the REBELS-IFU program at z~6.5-7.7, deriving HI column densities by modelling Ly$\alpha$ damping wings. Significant damped Ly$\alpha$ absorption is detected in eight galaxies, with $N_{\mathrm{HI}}\gtrsim10^{21}$ cm$^{-2}$. We use the column densities and sizes derived for these sources to estimate their HI mass and compare with $L_{\mathrm{[CII]}}$-$M_{\mathrm{HI}}$ calibrations. The resulting HI masses show a tentative correlation with those inferred from [CII], although the [CII]-based estimates are systematically larger, suggesting that the HI reservoirs may extend beyond the [CII]-emitting gas. We also combine the DLA-based measurements with FIR-derived dust-to-gas ratios, dust attenuation, and gas-phase metallicities. No correlation is found between DLA-based and FIR-based dust-to-gas ratios, but combining the REBELS-IFU sample with literature samples at lower metallicities reveals a strong correlation between $A_{\mathrm{V}}/N_{\mathrm{HI}}$ and metallicity. These findings suggest that by $z\sim7$ massive galaxies can already host substantial, enriched reservoirs of neutral gas and dust, consistent with $A_{\mathrm{V}}$/$N_{\mathrm{HI}}$-metallicity trends at lower redshift. At the highest redshifts ($z>8$), however, we see tentative evidence for systematically lower $A_{\mathrm{V}}$/$N_{\mathrm{HI}}$ at fixed metallicity, which may point to pristine gas accretion or more efficient dust destruction/expulsion.

In the context of transient constant-roll inflation near a local maximum, we derive the non-perturbative field redefinition that relates a Gaussian random field with the true non-Gaussian curvature perturbation. Our analysis shows the emergence of a new critical amplitude $\zeta_*$, corresponding to perturbations that prevent the inflaton from overshooting the local maximum, thus becoming trapped in the false minimum of the potential. For potentials with a mild curvature at the local maximum (and thus small non-Gaussianity), we recover the known perturbative field redefinition. We apply these results to the formation of primordial black holes, and discuss the cases for which $\zeta_*$ is smaller or of the same order than the critical value for collapse of spherically symmetric overdensities. In the latter case, we present a simple potential for which the power spectrum needs an amplitude 10 times smaller that in the Gaussian case for producing a sizeable amount of primordial black holes.

J-PAS (Javalambre Physics of the Accelerating Universe Astrophysical Survey) will present a groundbreaking photometric survey covering 8500 deg$^2$ of the visible sky from Javalambre, capturing data in 56 narrow band filters. This survey promises to revolutionize galaxy evolution studies by observing $\sim$10$^8$ galaxies with low spectral resolution. A crucial aspect of this analysis involves predicting stellar population parameters from the observed galaxy photometry. In this study, we combine the exquisite J-PAS photometry with state-of-the-art single stellar population (SSP) libraries to accurately predict stellar age, metallicity, and dust attenuation with a neural network (NN) model. The NN is trained on synthetic J-PAS photometry from different SSP librares (E-MILES, Charlot & Bruzual, XSL), to enhance the robustness of our predictions against individual SSP model variations and limitations. To create mock samples with varying observed magnitudes we add artificial noise in the form of random Gaussian variations within typical observational uncertainties in each band. Our results indicate that the NN can accurately estimate stellar parameters for SSP models without evident degeneracies, surpassing a bayesian SED-fitting method on the same test set. We obtain median bias, scatter and percentage of outliers $\mu$ = (0.01 dex, 0.00 dex, 0.00 mag), $\sigma_{NMAD}$ = (0.23 dex, 0.29 dex, 0.04 mag), f$_{o}$ = (17 %, 24 %, 1 %) at $i \sim$17 mag for age, metallicity and dust attenuation, respectively. The accuracy of the predictions is highly dependent on the signal-to-noise (S/N) ratio of the photometry, achieving robust predictions up to $i$ $\sim$ 20 mag.

Gravitational lensing allows the detection of binary black holes (BBH) at cosmological distances with chirp masses that appear to be enhanced by $1+z$ in the range $195\%$), from which we infer a factor $\simeq 5$ higher intrinsic rate of NSBH events than BBH events, reflecting a higher proportion of neutron stars formed by early star formation. We predict a distinctive locus for lensed NSBH events in the observed binary mass plane, spanning $1

The lack of power anomaly is an unexpected feature observed at large angular scales in the CMB maps produced by the COBE, WMAP and Planck satellites. This signature, which consists in a missing of power with respect to that predicted by the LCDM model, might hint at a new cosmological phase before the standard inflationary era. The main point of this paper is taking into account the latest Planck polarisation data to investigate how the CMB polarisation improves the understanding of this feature. With this aim, we apply to the latest Planck data, both PR3 (2018) and PR4 (2020) releases, a new class of estimators capable of evaluating this anomaly by considering temperature and polarisation data both separately and in a jointly way. This is the first time that the PR4 dataset has been used to study this anomaly. To critically evaluate this feature, taking into account the residuals of known systematic effects present in the Planck datasets, we analyse the cleaned CMB maps using different combinations of sky masks, harmonic range and binning on the CMB multipoles. Our analysis shows that the estimator based only on temperature data confirms the presence of a lack of power with a lower-tail-probability (LTP), depending on the component separation method, $\leq 0.33\%$ and $\leq 1.76\%$ for PR3 and PR4, respectively. To our knowledge, the LTP$\leq 0.33\%$ for the PR3 dataset is the lowest one present in the literature obtained from Planck 2018 data, considering the Planck confidence mask. We find significant differences between these two datasets when polarisation is taken into account. Moreover, we also show that for the PR3 dataset the inclusion of the subdominant polarisation information provides estimates that are less likely accepted in a LCDM cosmological model than the only-temperature analysis over the entire harmonic-range considered.

Recent James Webb Space Telescope (JWST) observations have revealed a population of sources with a compact morphology and a `v-shaped' continuum, namely blue at rest-frame \lambda&lt;4000A and red at longer wavelengths. The nature of these sources, called `little red dots' (LRDs), is still debated, since it is unclear if they host active galactic nuclei (AGN) and their number seems to drastically drop at z<4. We utilise the 63 $deg^2$ covered by the quick Euclid Quick Data Release (Q1) to extend the search for LRDs to brighter magnitudes and to lower z than what has been possible with JWST to have a broader view of the evolution of this peculiar galaxy population. The selection is done by fitting the available photometric data (Euclid, Spitzer/IRAC, and ground-based griz data) with two power laws, to retrieve the rest-frame optical and UV slopes consistently over a large redshift range (i.e, z<7.6). We exclude extended objects and possible line emitters, and perform a visual inspection to remove imaging artefacts. The final selection includes 3341 LRD candidates from z=0.33 to z=3.6, with 29 detected in IRAC. Their rest-frame UV luminosity function, in contrast with previous JWST studies, shows that the number density of LRD candidates increases from high-z down to z=1.5-2.5 and decreases at even lower z. Less evolution is apparent focusing on the subsample of more robust LRD candidates having IRAC detections, which is affected by low statistics and limited by the IRAC resolution. The comparison with previous quasar UV luminosity functions shows that LRDs are not the dominant AGN population at z<4. Follow-up studies of these LRD candidates are key to confirm their nature, probe their physical properties and check for their compatibility with JWST sources, since the different spatial resolution and wavelength coverage of Euclid and JWST could select different samples of compact sources.

We report on a search for sub-GeV dark matter (DM) particles interacting with electrons using the DAMIC-M prototype detector at the Modane Underground Laboratory. The data feature a significantly lower detector single $e^-$ rate (factor 50) compared to our previous search, while also accumulating a ten times larger exposure of $\sim$1.3 kg-day. DM interactions in the skipper charge-coupled devices (CCDs) are searched for as patterns of two or three consecutive pixels with a total charge between 2 and 4 $e^-$. We find 144 candidates of 2 $e^-$ and 1 candidate of 4 $e^-$, where 141.5 and 0.071, respectively, are expected from background. With no evidence of a DM signal, we place stringent constraints on DM particles with masses between 1 and 1000 MeV/$c^2$ interacting with electrons through an ultra-light or heavy mediator. For large ranges of DM masses below 1 GeV/c$^2$, we exclude theoretically-motivated benchmark scenarios where hidden-sector particles are produced as a major component of DM in the Universe through the freeze-in or freeze-out mechanisms.

We search for the signature of parity-violating physics in the cosmic microwave background, called cosmic birefringence, using the Planck data release 4. We initially find a birefringence angle of $\beta=0.30\pm0.11$ (68% C.L.) for nearly full-sky data. The values of $\beta$ decrease as we enlarge the Galactic mask, which can be interpreted as the effect of polarized foreground emission. Two independent ways to model this effect are used to mitigate the systematic impact on $\beta$ for different sky fractions. We choose not to assign cosmological significance to the measured value of $\beta$ until we improve our knowledge of the foreground polarization.

The data volume generated by astronomical surveys is growing rapidly. Traditional analysis techniques in spectroscopy either demand intensive human interaction or are computationally expensive. In this scenario, machine learning, and unsupervised clustering algorithms in particular offer interesting alternatives. The Apache Point Observatory Galactic Evolution Experiment (APOGEE) offers a vast data set of near-infrared stellar spectra which is perfect for testing such alternatives. Apply an unsupervised classification scheme based on $K$-means to the massive APOGEE data set. Explore whether the data are amenable to classification into discrete classes. We apply the $K$-means algorithm to 153,847 high resolution spectra ($R\approx22,500$). We discuss the main virtues and weaknesses of the algorithm, as well as our choice of parameters. We show that a classification based on normalised spectra captures the variations in stellar atmospheric parameters, chemical abundances, and rotational velocity, among other factors. The algorithm is able to separate the bulge and halo populations, and distinguish dwarfs, sub-giants, RC and RGB stars. However, a discrete classification in flux space does not result in a neat organisation in the parameters space. Furthermore, the lack of obvious groups in flux space causes the results to be fairly sensitive to the initialisation, and disrupts the efficiency of commonly-used methods to select the optimal number of clusters. Our classification is publicly available, including extensive online material associated with the APOGEE Data Release 12 (DR12). Our description of the APOGEE database can enormously help with the identification of specific types of targets for various applications. We find a lack of obvious groups in flux space, and identify limitations of the $K$-means algorithm in dealing with this kind of data.

Recent JWST observations have revealed an unexpected abundance of massive galaxy candidates in the early Universe, extending further in redshift and to lower luminosity than what had previously been found by sub-millimeter surveys. These JWST candidates have been interpreted as challenging the $\Lambda$CDM cosmology, but, so far, they have mostly relied only on rest-frame ultraviolet data and lacked spectroscopic confirmation of their redshifts. Here we report a systematic study of 36 massive dust-obscured galaxies with spectroscopic redshifts between $z_{\rm spec}=5-9$ from the JWST FRESCO survey. We find no tension with the $\Lambda$CDM model in our sample. However, three ultra-massive galaxies (log$M_{\star}/M_{\odot}$ $\gtrsim11.0$) require an exceptional fraction of 50% of baryons converted into stars -- two to three times higher than even the most efficient galaxies at later epochs. The contribution from an active nucleus is unlikely because of their extended emission. Ultra-massive galaxies account for as much as 17% of the total cosmic star formation rate density at $z\sim5-6$.

We present a set of preliminary simulations of intermediate mass ratio inspirals (IMRIs) inside dark matter (DM) spikes accounting for post-Newtonian corrections the interaction between the two black holes up to the order 2.5 in $c^2$, as well as relativistic corrections to the dynamical friction (DF) force exerted by the DM distribution. We find that, incorporating relativity reduces of a factor $1/2$ the inspiral time, for equivalent initial orbital parameters, with respect to the purely classical estimates. Vice versa, neglecting the DF of the spike systematically yields longer inspiral times.

The origin of small deviations from statistical isotropy in the Cosmic Microwave Background (CMB) - the so-called CMB anomalies - remains an open question in modern cosmology. In this work, we test statistical isotropy in Planck Data Release 4 (PR4) by estimating the temperature and $E$-mode power spectra across independent sky regions. We find that the directions with higher local bandpower amplitudes in intensity are clustered for multipoles between 200 and 2000 with clustering probabilities consistently below 1% for all these scales when compared to end-to-end (E2E) Planck simulations; notably, this range extends beyond that reported in Planck Data Release 3 (PR3). On the other hand, no significant clustering is observed in the polarization $E$-modes. In a complementary analysis, we search for dipolar variations in cosmological parameters fitted using the previously computed power spectra. When combining temperature and polarization power spectra, we identify a potential anomaly in the amplitude of the primordial power spectrum, $A_{s}$, with only 5 out of 600 simulations exhibiting a dipole amplitude as large as that observed in the data. Interestingly, the dipole direction aligns closely with the known hemispherical power asymmetry, suggesting a potential link between these anomalies. All other cosmological parameters remain consistent with $\Lambda\mathrm{CDM}$ expectations. Our findings highlight the need to further investigate these anomalies and understand their nature and potential implications for better understanding of the early Universe.