As the statistical precision of cosmological measurements increases, the accuracy of the theoretical description of these measurements needs to increase correspondingly in order to infer the underlying cosmology that governs the Universe. To this end, we have created the Cosmology Likelihood for Observables in Euclid (CLOE), which is a novel cosmological parameter inference pipeline developed within the Euclid Consortium to translate measurements and covariances into cosmological parameter constraints. In this first in a series of six papers, we describe the theoretical recipe of this code for the Euclid primary probes. These probes are composed of the photometric 3x2pt observables of cosmic shear, galaxy-galaxy lensing, and galaxy clustering, along with spectroscopic galaxy clustering. We provide this description in both Fourier and configuration space for standard and extended summary statistics, including the wide range of systematic uncertainties that affect them. This includes systematic uncertainties such as intrinsic galaxy alignments, baryonic feedback, photometric and spectroscopic redshift uncertainties, shear calibration uncertainties, sample impurities, photometric and spectroscopic galaxy biases, as well as magnification bias. The theoretical descriptions are further able to accommodate both Gaussian and non-Gaussian likelihoods and extended cosmologies with non-zero curvature, massive neutrinos, evolving dark energy, and simple forms of modified gravity. These theoretical descriptions that underpin CLOE will form a crucial component in revealing the true nature of the Universe with next-generation cosmological surveys such as Euclid.

Tidal disruption events (TDEs) are expected to release much of their energy in the far-ultraviolet (UV), which we do not observe directly. However, infrared (IR) observations can observe re-radiation of the optical/UV emission from dust, and if this dust is observed in the process of sublimation, we can infer the un-observed UV radiated energy. TDEs have also been predicted to show spectra shallower than a blackbody in the IR, but this has not yet been observed. We present near/mid-IR observations of the TDE AT 2019azh spanning from -3 d before peak until >1750 d after. We evaluate these observations for consistency with dust emission or direct emission from the TDE. We fit the IR data with a modified blackbody associated with dust emission. The UV+optical+IR data are compared with simulated spectra produced from general relativistic radiation magnetohydrodynamics simulations of super-Eddington accretion. We model the data at later times (> 200 d) as an IR echo. The IR data at the maximum light can not be self-consistently fit with dust emission. Instead, the data can be better fit with a reprocessing model, with the IR excess arising due to the absorption opacity being dominated by free-free processes in the dense reprocessing envelope. We infer a large viewing angle of $\sim$60$^{\circ}$, consistent with previously reported X-ray observations, and a tidally disrupted star with mass > 2 M$_{\odot}$. The IR emission at later times is consistent with cool dust emission. We model these data as an IR echo and find that the dust is distant (0.65 pc), and clumpy, with a low covering factor. We show that TDEs can have an IR excess not arising from dust and that IR observations at early times can constrain the viewing angle for the TDE in the unified model. Near-IR observations are therefore essential to distinguish between hot dust and a non-thermal IR excess.

Long-period radio transients (LPTs) are an emerging group of radio transients that show periodic polarised radio bursts with periods varying from a few minutes to a few hours. Fewer than a dozen LPTs have been detected so far, and their origin (source and emission mechanism) remains unclear. Here, we report the discovery of a 1.5 hr LPT, ASKAP J144834-685644, adding to the current sample of sources. ASKAP J144834-685644 is one of the very few LPTs that has been detected from X-rays to radio. It shows a steep radio spectrum and polarised radio bursts, which resemble the radio emission in known LPTs. In addition, it also shows highly structured and periodic narrow-band radio emission. Multi-wavelength properties suggest that the spectral energy distribution (SED) peaks at near ultraviolet wavelengths, indicating the presence of a hot magnetic source. Combining multi-wavelength information, we infer that ASKAP J144834-685644 may be a near edge-on magnetic white dwarf binary (MWD), although we can not fully rule out ASKAP J144834-685644 being an isolated white dwarf pulsar or even a transitional millisecond pulsar (despite the lack of radio pulsations). If ASKAP J144834-685644 is an MWD binary, the observed broadband spectral energy distribution can be explained by emission from an accretion disk. This hints that some fraction of optically bright LPTs may be accreting binaries with the radio period being the orbital period. It might further suggest a connection between optically bright non-accreting synchronized LPTs, such as polars, and non-accreting asynchronous WD pulsars, such as AR Sco and J1912-4410.

ESA's Euclid cosmology mission relies on the very sensitive and accurately calibrated spectroscopy channel of the Near-Infrared Spectrometer and Photometer (NISP). With three operational grisms in two wavelength intervals, NISP provides diffraction-limited slitless spectroscopy over a field of $0.57$ deg$^2$. A blue grism $\text{BG}_\text{E}$ covers the wavelength range $926$--$1366$\,nm at a spectral resolution $R=440$--$900$ for a $0.5''$ diameter source with a dispersion of $1.24$ nm px$^{-1}$. Two red grisms $\text{RG}_\text{E}$ span $1206$ to $1892$\,nm at $R=550$--$740$ and a dispersion of $1.37$ nm px$^{-1}$. We describe the construction of the grisms as well as the ground testing of the flight model of the NISP instrument where these properties were established.

We present a cosmological analysis using the second and third moments of the weak lensing mass (convergence) maps from the first three years of data (Y3) data of the Dark Energy Survey (DES). The survey spans an effective area of 4139 square degrees and uses the images of over 100 million galaxies to reconstruct the convergence field. The second moment of the convergence as a function of smoothing scale contains information similar to standard shear 2-point statistics. The third moment, or the skewness, contains additional non-Gaussian information. The data is analysed in the context of the $\Lambda$CDM model, varying 5 cosmological parameters and 19 nuisance parameters modelling astrophysical and measurement systematics. Our modelling of the observables is completely analytical, and has been tested with simulations in our previous methodology study. We obtain a 1.7\% measurement of the amplitude of fluctuations parameter $S_8\equiv \sigma_8 (\Omega_m/0.3)^{0.5} = 0.784\pm 0.013$. The measurements are shown to be internally consistent across redshift bins, angular scales, and between second and third moments. In particular, the measured third moment is consistent with the expectation of gravitational clustering under the $\Lambda$CDM model. The addition of the third moment improves the constraints on $S_8$ and $\Omega_{\rm m}$ by $\sim$15\% and $\sim$25\% compared to an analysis that only uses second moments. We compare our results with {\it Planck} constraints from the Cosmic Microwave Background (CMB), finding a $2.2$ \textendash $2.8\sigma$ tension in the full parameter space, depending on the combination of moments considered. The third moment independently is in $2.8\sigma$ tension with {\it Planck}, and thus provides a cross-check on analyses of 2-point correlations.

Mass loss is a key aspect of stellar evolution, particularly in evolved massive stars, yet episodic mass loss remains poorly understood. To investigate this, we need evolved massive stellar populations across various galactic environments. However, spectral classifications are challenging to obtain in large numbers, especially for distant galaxies. We addressed this by leveraging machine-learning techniques. We combined \textit{Spitzer} photometry and Pan-STARRS1 optical data to classify point sources in 26 galaxies within 5 Mpc, and a metallicity range 0.07-1.36 Z$_\odot$. \textit{Gaia} DR3 astrometry was used to remove foreground sources. Classifications are derived using a machine-learning model developed by Maravelias et al. (2022). We report classifications for 1,147,650 sources, with 276,657 sources ($\sim24\%$) being robust. Among these are 120,479 Red Supergiants (RSGs; $\sim11\%$). The classifier performs well even at low metallicities ($\sim0.1$ Z$_\odot$) and distances under 1.5 Mpc, with a slight decrease in accuracy beyond $\sim3$ Mpc due to \textit{Spitzer}'s resolution limits. We also identified 21 luminous RSGs ($\textrm{log}(L/L_\odot)\ge5.5$), 159 dusty Yellow Hypergiants in M31 and M33, as well as 6 extreme RSGs ($\textrm{log}(L/L_\odot)\ge6$) in M31, challenging observed luminosity limits. Class trends with metallicity align with expectations, though biases exist. This catalog serves as a valuable resource for individual-object studies and \textit{James Webb} Space Telescope target selection. It enables follow-up on luminous RSGs and Yellow Hypergiants to refine our understanding of their evolutionary pathways. Additionally, we provide the largest spectroscopically confirmed catalog of massive stars and candidates to date, comprising 5,273 sources (including $\sim330$ other objects).

In any imaging survey, measuring accurately the astronomical background light is crucial to obtain good photometry. This paper introduces BKGnet, a deep neural network to predict the background and its associated error. BKGnet has been developed for data from the Physics of the Accelerating Universe Survey (PAUS), an imaging survey using a 40 narrow-band filter camera (PAUCam). Images obtained with PAUCam are affected by scattered light: an optical effect consisting of light multiply that deposits energy in specific detector regions contaminating the science measurements. Fortunately, scattered light is not a random effect, but it can be predicted and corrected for. We have found that BKGnet background predictions are very robust to distorting effects, while still being statistically accurate. On average, the use of BKGnet improves the photometric flux measurements by 7% and up to 20% at the bright end. BKGnet also removes a systematic trend in the background error estimation with magnitude in the i-band that is present with the current PAU data management method. With BKGnet, we reduce the photometric redshift outlier rate

We present a new infrared survey covering the three Euclid deep fields and four other Euclid calibration fields using Spitzer's Infrared Array Camera (IRAC). We have combined these new observations with all relevant IRAC archival data of these fields in order to produce the deepest possible mosaics of these regions. In total, these observations represent nearly 11% of the total Spitzer mission time. The resulting mosaics cover a total of approximately 71.5deg$^2$ in the 3.6 and 4.5um bands, and approximately 21.8deg$^2$ in the 5.8 and 8um bands. They reach at least 24 AB magnitude (measured to sigma, in a 2.5 arcsec aperture) in the 3.6um band and up to ~ 5 mag deeper in the deepest regions. The astrometry is tied to the Gaia astrometric reference system, and the typical astrometric uncertainty for sources with 16<[3.6]<19 is <0.15 arcsec. The photometric calibration is in excellent agreement with previous WISE measurements. We have extracted source number counts from the 3.6um band mosaics and they are in excellent agreement with previous measurements. Given that the Spitzer Space Telescope has now been decommissioned these mosaics are likely to be the definitive reduction of these IRAC data. This survey therefore represents an essential first step in assembling multi-wavelength data on the Euclid deep fields which are set to become some of the premier fields for extragalactic astronomy in the 2020s.

We characterise the selection cuts and clustering properties of a magnitude-limited sample of bright galaxies that is part of the Bright Galaxy Survey (BGS) of the Dark Energy Spectroscopic Instrument (DESI) using the ninth data release of the Legacy Imaging Surveys (DR9). We describe changes in the DR9 selection compared to the DR8 one as explored in Ruiz-Macias et al. (2021). We also compare the DR9 selection in three distinct regions: BASS/MzLS in the north Galactic Cap (NGC), DECaLS in the NGC, and DECaLS in the south Galactic Cap (SGC). We investigate the systematics associated with the selection and assess its completeness by matching the BGS targets with the Galaxy and Mass Assembly (GAMA) survey. We measure the angular clustering for the overall bright sample (r $\leq$ 19.5) and as function of apparent magnitude and colour. This enables to determine the clustering strength and slope by fitting a power-law model that can be used to generate accurate mock catalogues for this tracer. We use a counts-in-cells technique to explore higher-order statistics and cross-correlations with external spectroscopic data sets in order to check the evolution of the clustering with redshift and the redshift distribution of the BGS targets using clustering-redshifts. While this work validates the properties of the BGS bright targets, the final target selection pipeline and clustering properties of the entire DESI BGS will be fully characterised and validated with the spectroscopic data of Survey Validation.

Recently several long-period radio transients have been discovered, with strongly polarised coherent radio pulses appearing on timescales between tens to thousands of seconds [1,2]. In some cases the radio pulses have been interpreted as coming from rotating neutron stars with extremely strong magnetic fields, known as magnetars; the origin of other, occasionally periodic and less well-sampled radio transients, is still debated [3]. Coherent periodic radio emission is usually explained by rotating dipolar magnetic fields and pair production mechanisms, but such models do not easily predict radio emission from such slowly-rotating neutron stars and maintain it for extended times. On the other hand, highly magnetic isolated white dwarfs would be expected to have long spin periodicities, but periodic coherent radio emission has not yet been directly detected from these sources. Here we report observations of a long-period (21 minutes) radio transient, which we have labeled GPMJ1839-10. The pulses vary in brightness by two orders of magnitude, last between 30 and 300 seconds, and have quasi-periodic substructure. The observations prompted a search of radio archives, and we found that the source has been repeating since at least 1988. The archival data enabled constraint of the period derivative to &lt;3.6\times10^{-13}s s$^{-1}$, which is at the very limit of any classical theoretical model that predicts dipolar radio emission from an isolated neutron star.

We describe distinctive stellar features indicating the presence of hyperons in neutron stars. A strongly negative curvature of the mass-radius relation R(M) is characteristic of hyperons, which can be determined from measurements of neutron stars with three different masses. Similarly, a reduced second derivative of the tidal deformability as function of mass $\lambda(M)$ points to hyperonic degrees of freedom in NS matter. The slopes of such curves R(M) and \lambda(M) can distinguish a hyperonic equation of state from purely nucleonic models if they appear increased (decreased for \lambda(M)) relative to the maximum mass of neutron stars.

Neutron stars are the dense and highly magnetic relics of supernova explosions of massive stars. The quest to constrain the Equation of State (EoS) of ultra-dense matter and thereby probe the behavior of matter inside neutron stars, is one of the core goals of modern physics and astrophysics. A promising method involves investigating the long-term cooling of neutron stars, and comparing theoretical predictions with various sources at different ages. However, limited observational data, and uncertainties in source ages and distances, have hindered this approach. In this work, re-analyzing XMM-Newton and Chandra data from dozens of thermally emitting isolated neutron stars, we have identified three sources with unexpectedly cold surface temperatures for their young ages. To investigate these anomalies, we conducted magneto-thermal simulations across diverse mass and magnetic fields, considering three different EoS. We found that the "minimal" cooling model, failed to explain the observations, regardless the mass and the magnetic field, as validated by a machine learning classification method. The existence of these young cold neutron stars suggests that any dense matter EoS must be compatible with a fast cooling process at least in certain mass ranges, eliminating a significant portion of current EoS options according to recent meta-modelling analysis.

Ultraluminous X-ray binaries have challenged our assumptions of extreme accretion rates in X-ray binaries, and impact other subfields of astronomy, such as cosmology, gravitational wave sources and supernovæ. Our understanding of ULXs has changed tremendously over the last 35 years, and we now know that ULXs can be powered by accreting neutron stars as well as black holes, and can be found in a wide range of stellar environments. In this chapter, we introduce the observational techniques used to discover and characterize ULXs, and discuss our current understanding of their unique accretion physics and formation channels.

We use a recent census of the Milky Way (MW) satellite galaxy population to constrain the lifetime of particle dark matter (DM). We consider two-body decaying dark matter (DDM) in which a heavy DM particle decays with lifetime $\tau$ comparable to the age of the Universe to a lighter DM particle (with mass splitting $\epsilon$) and to a dark radiation species. These decays impart a characteristic "kick velocity," $V_{\mathrm{kick}}=\epsilon c$, on the DM daughter particles, significantly depleting the DM content of low-mass subhalos and making them more susceptible to tidal disruption. We fit the suppression of the present-day DDM subhalo mass function (SHMF) as a function of $\tau$ and $V_{\mathrm{kick}}$ using a suite of high-resolution zoom-in simulations of MW-mass halos, and we validate this model on new DDM simulations of systems specifically chosen to resemble the MW. We implement our DDM SHMF predictions in a forward model that incorporates inhomogeneities in the spatial distribution and detectability of MW satellites and uncertainties in the mapping between galaxies and DM halos, the properties of the MW system, and the disruption of subhalos by the MW disk using an empirical model for the galaxy--halo connection. By comparing to the observed MW satellite population, we conservatively exclude DDM models with \tau &lt; 18\ \mathrm{Gyr} ($29\ \mathrm{Gyr}$) for $V_{\mathrm{kick}}=20\ \mathrm{km}\, \mathrm{s}^{-1}$ ($40\ \mathrm{km}\, \mathrm{s}^{-1}$) at $95\%$ confidence. These constraints are among the most stringent and robust small-scale structure limits on the DM particle lifetime and strongly disfavor DDM models that have been proposed to alleviate the Hubble and $S_8$ tensions.

We observed the periodic radio transient GLEAM-X J162759.5-523504.3 (GLEAM-X J1627) using the Chandra X-ray Observatory for about 30-ks on January 22-23, 2022, simultaneously with radio observations from MWA, MeerKAT and ATCA. Its radio emission and 18-min periodicity led the source to be tentatively interpreted as an extreme magnetar or a peculiar highly magnetic white dwarf. The source was not detected in the 0.3-8 keV energy range with a 3-sigma upper-limit on the count rate of 3x10^{-4} counts/s. No radio emission was detected during our X-ray observations either. Furthermore, we studied the field around GLEAM-X J1627 using archival ESO and DECam data, as well as recent SALT observations. Many sources are present close to the position of GLEAM-X J1627, but only two within the 2" radio position uncertainty. Depending on the assumed spectral distribution, the upper limits converted to an X-ray luminosity of L_{X}<6.5x10^{29} erg/s for a blackbody with temperature kT=0.3 keV, or L_{X}<9x10^{29} erg/s for a power-law with photon index Gamma = 2 (assuming a 1.3 kpc distance). Furthermore, we performed magneto-thermal simulations for neutron stars considering crust- and core-dominated field configurations. Based on our multi-band limits, we conclude that: i) in the magnetar scenario, the X-ray upper limits suggest that GLEAM-X J1627 should be older than ~1 Myr, unless it has a core-dominated magnetic field or has experienced fast-cooling; ii) in the white dwarf scenario, we can rule out most binary systems, a hot sub-dwarf and a hot magnetic isolated white dwarf (T>10.000 K), while a cold isolated white dwarf is still compatible with our limits.

The nature of two recently discovered radio emitters with unusually long periods of 18min (GLEAM-X J1627-52) and 21min (GPM J1839-10) is highly debated. Their bright radio emission resembles that of radio magnetars, but their long periodicities and lack of detection at other wavelengths challenge the neutron-star interpretation. In contrast, long rotational periods are common in white dwarfs but, although predicted, dipolar radio emission from isolated magnetic white dwarfs has never been unambiguously observed. In this work, we investigate these long-period objects as potential isolated neutron-star or white-dwarf dipolar radio emitters and find that both scenarios pose significant challenges to our understanding of radio emission via pair production in dipolar magnetospheres. We also perform population-synthesis simulations based on dipolar spin-down in both pictures, assuming different initial-period distributions, masses, radii, beaming fractions, and magnetic-field prescriptions, to assess their impact on the ultra-long pulsar population. In the neutron-star scenario, we do not expect a large number of ultra-long period pulsars under any physically motivated (or even extreme) assumptions for the period evolution. On the other hand, in the white-dwarf scenario, we can easily accommodate a large population of long-period radio emitters. However, no mechanism can easily explain the production of such bright coherent radio emission in either scenarios.

We present photometric redshifts for 1 341 559 galaxies from the Physics of the Accelerating Universe Survey (PAUS) over 50.38 ${\rm deg}^{2}$ of sky to $i_{\rm AB}=23$. Redshift estimation is performed using DEEPz, a deep-learning photometric redshift code. We analyse the photometric redshift precision when varying the photometric and spectroscopic samples. Furthermore, we examine observational and instrumental effects on the precision of the photometric redshifts, and we compare photometric redshift measurements with those obtained using a template method-fitting BCNz2. Finally, we examine the use of photometric redshifts in the identification of close galaxy pairs. We find that the combination of samples from W1+W3 in the training of DEEPz significantly enhances the precision of photometric redshifts. This also occurs when we recover narrow band fluxes using broad bands measurements. We show that DEEPz determines the redshifts of galaxies in the prevailing spectroscopic catalogue used in the training of DEEPz with greater precision. For the faintest galaxies ($i_{\rm AB}=21-23$), we find that DEEPz improves over BCNz2 both in terms of the precision (20-50 per cent smaller scatter) and in returning a smaller outlier fraction in two of the wide fields. The catalogues were tested for the identification of close galaxy pairs, showing that DEEPz is effective for the identification of close galaxy pairs for samples with i_{\rm AB} &lt; 22.5 and redshift 0.2 &lt; z &lt; 0.6. In addition, identifying close galaxy pairs common between DEEPz and BCNz2 is a promising approach to improving the purity of the catalogues of these systems.

One of the potential explanations for the existence of very large grains (VLGs) in the inner envelope of low/intermediate-mass Class 0/I Young Stellar Object is the migration of VLGs from the protostellar disk via a protostellar outflow. To understand whether the grain migration is prevented by RAdiative Torque Disruption (RATD), we perform the numerical modeling of RATD in parallel with the grain propagation, using the gas velocity and density structure inside the jet and outflow from an MHD simulation of an intermediate Class 0 protostar. We found that with the bolometric luminosity $\geq 20L_{\odot}$, RATD can destroy aggregate grains of size $1 \sim 500\rm \mu m$ having maximum tensile strength $S_{\rm max} \leq 10^{5} \rm erg cm^{-3}$ inside the jet/outflow base after &lt; 2 yrs. This effect lets sub-micron grains dominate the outflow and partially prevent the migration of large grains from the inner disk to inner envelope. In contrast, RATD cannot prevent the migration of composite VLGs and submillimeter grains having $S_{\rm max}\geq 10^{7} \rm erg cm^{-3}$. Next, we incorporate RATD into POLARIS, assuming grains are not moving relative to the gas. We found that POLARIS works well in describing the disruption for aggregate grains, but overestimates the disruption efficiency for composite grains. The observed polarization degree can be reduced by twice when aggregate grains are removed from the outflow cavity wall and inner envelope by RATD. However, RATD is not an important factor controlling dust polarization properties as iron inclusions do.

We report the discovery and characterization of ZTF J0112+5827, a new magnetic cataclysmic variable with an orbital period of 80.9 minutes. ROSAT observations revealed X-ray emission with an average flux of $(68.4 \pm 15.7) \times 10^{-14}$erg s$^{-1}$cm$^{-2}$ (0.1--2.4 keV). The ZTF light curves show ellipsoidal-like variability in the $g$ band and two prominent humps at phases $\sim$0.0 and $\sim$0.7 in $i$ and $r$ bands. Spectroscopic observations with the Palomar 200-inch telescope revealed cyclotron emission features and strong He II and Balmer emission lines. Doppler tomography shows clear accretion streams with line-of-sight velocities of $\sim$500 km s$^{-1}$, but no accretion disk. Analysis of cyclotron harmonics indicates a magnetic field strength of $38.7^{+1.3}_{-1.1}$ MG, confirming ZTF J0112+5827 as a polar system containing a strongly magnetic white dwarf.