We investigate wormhole solutions within the framework of the semi-classical Einstein equations in the presence of the conformal anomaly (or trace anomaly). These solutions are sourced by a stress-energy tensor (SET) derived from the trace anomaly, and depend on two positive coefficients, $\alpha$ and $\lambda$, determined by the matter content of the theory and on the degrees of freedom of the involved quantum fields. For a Type B anomaly ($\alpha=0$), we obtain wormhole geometries assuming a constant redshift function and show that the SET components increase with the parameter $\lambda$. In the case of a Type A anomaly ($\lambda=0$), we generalize previously known solutions, yielding a family of geometries that includes Lorentzian wormholes, naked singularities, and the Schwarzschild black hole. Using isotropic coordinates, we identify parameter choices that produce traversable wormhole solutions. Extending to the full trace-anomaly contribution, we solve the differential equation near the throat to obtain the redshift function and demonstrate that both the Ricci and Kretschmann scalars remain finite at the throat. We further analyze the trajectories of null and timelike particles, showing that the height and width of the effective potential for null geodesics increase monotonically with $\alpha$, while the innermost stable circular orbit (ISCO) radius also grows with larger $\alpha$. These results illustrate the rich interplay between trace anomaly effects and the structure and dynamics of wormhole spacetimes.

The distribution of close-in exoplanets is shaped by the interplay between atmospheric and dynamical processes. The Neptunian Desert, Ridge, and Savanna illustrate the sensitivity of these worlds to such processes, making them ideal to disentangle their roles. Determining how many Neptunes were brought close-in by early disk-driven migration (DDM; maintaining primordial spin-orbit alignment) or late high-eccentricity migration (HEM; generating large misalignments) is essential to understand how much atmosphere they lost. We propose a unified view of the Neptunian landscape to guide its exploration, speculating that the Ridge is a hot spot for evolutionary processes. Low-density Neptunes would mainly undergo DDM, getting fully eroded at shorter periods than the Ridge, while denser Neptunes would be brought to the Ridge and Desert by HEM. We embark on this exploration via ATREIDES, which relies on spectroscopy and photometry of 60 close-in Neptunes, their reduction with robust pipelines, and their interpretation through internal structure, atmospheric, and evolutionary models. We carried out a systematic RM census with VLT/ESPRESSO to measure the distribution of 3D spin-orbit angles, correlate its shape with system properties and thus relate the fraction of aligned-misaligned systems to DDM, HEM, and atmospheric erosion. Our first target, TOI-421c, lies in the Savanna with a neighboring sub-Neptune TOI-421b. We measured their 3D spin-orbit angles (Psib = 57+11-15 deg; Psic = 44.9+4.4-4.1 deg). Together with the eccentricity and possibly large mutual inclination of their orbits, this hints at a chaotic dynamical origin that could result from DDM followed by HEM. ATREIDES will provide the community with a wealth of constraints for formation and evolution models. We welcome collaborations that will contribute to pushing our understanding of the Neptunian landscape forward.

The standard $\Lambda$CDM paradigm of the physical Universe suffers from well-known conceptual problems and is challenged by observational data. Alternative models exist in the literature, both phenomenological and physically motivated, many of them suffering from similar or new problems. We propose a method to mechanically generate alternative models in a procedure informed by data and tuned to mitigate specific problems. We implement a computational framework, dubbed CosmoGen, that integrates evolutionary algorithms for symbolic regression, with computation of structure formation and background cosmological quantities that are used to guide the evolutionary process. As a proof-of-concept we apply the procedure to the specific case of dark energy fluid models, and ask the framework to generate models capable of alleviating the cosmological tensions $S_8$ and $H_0$. The system generates models with high fitness values, and through a Bayesian analysis of an illustrative model, we show that the model indeed alleviates the tensions.

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.

The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance.

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.

The Cosmological Principle (CP) -- the notion that the Universe is spatially isotropic and homogeneous on large scales -- underlies a century of progress in cosmology. It is conventionally formulated through the Friedmann-Lema\^itre-Robertson-Walker (FLRW) cosmologies as the spacetime metric, and culminates in the successful and highly predictive $\Lambda$-Cold-Dark-Matter ($\Lambda$CDM) model. Yet, tensions have emerged within the $\Lambda$CDM model, most notably a statistically significant discrepancy in the value of the Hubble constant, $H_0$. Since the notion of cosmic expansion determined by a single parameter is intimately tied to the CP, implications of the $H_0$ tension may extend beyond $\Lambda$CDM to the CP itself. This review surveys current observational hints for deviations from the expectations of the CP, highlighting synergies and disagreements that warrant further study. Setting aside the debate about individual large structures, potential deviations from the CP include variations of cosmological parameters on the sky, discrepancies in the cosmic dipoles, and mysterious alignments in quasar polarizations and galaxy spins. While it is possible that a host of observational systematics are impacting results, it is equally plausible that precision cosmology may have outgrown the FLRW paradigm, an extremely pragmatic but non-fundamental symmetry assumption.

We introduce SpectraPyle, a versatile spectral stacking pipeline developed for the Euclid mission's NISP spectroscopic surveys, aimed at extracting faint emission lines and spectral features from large galaxy samples in the Wide and Deep Surveys. Designed for computational efficiency and flexible configuration, SpectraPyle supports the processing of extensive datasets critical to Euclid's non-cosmological science goals. We validate the pipeline using simulated spectra processed to match Euclid's expected final data quality. Stacking enables robust recovery of key emission lines, including Halpha, Hbeta, [O III], and [N II], below individual detection limits. However, the measurement of galaxy properties such as star formation rate, dust attenuation, and gas-phase metallicity are biased at stellar mass below log10(M*/Msol) ~ 9 due to the flux-limited nature of Euclid spectroscopic samples, which cannot be overcome by stacking. The SFR-stellar mass relation of the parent sample is recovered reliably only in the Deep survey for log10(M*/Msol) > 10, whereas the metallicity-mass relation is recovered more accurately over a wider mass range. These limitations are caused by the increased fraction of redshift measurement errors at lower masses and fluxes. We examine the impact of residual redshift contaminants that arises from misidentified emission lines and noise spikes, on stacked spectra. Even after stringent quality selections, low-level contamination (< 6%) has minimal impact on line fluxes due to the systematically weaker emission of contaminants. Percentile-based analysis of stacked spectra provides a sensitive diagnostic for detecting contamination via coherent spurious features at characteristic wavelengths. While our simulations include most instrumental effects, real Euclid data will require further refinement of contamination mitigation strategies.

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.

The Euclid satellite will provide data on the clustering of galaxies and on the distortion of their measured shapes, which can be used to constrain and test the cosmological model. However, the increase in precision places strong requirements on the accuracy of the theoretical modelling for the observables and of the full analysis pipeline. In this paper, we investigate the accuracy of the calculations performed by the Cosmology Likelihood for Observables in Euclid (CLOE), a software able to handle both the modelling of observables and their fit against observational data for both the photometric and spectroscopic surveys of Euclid, by comparing the output of CLOE with external codes used as benchmark. We perform such a comparison on the quantities entering the calculations of the observables, as well as on the final outputs of these calculations. Our results highlight the high accuracy of CLOE when comparing its calculation against external codes for Euclid observables on an extended range of operative cases. In particular, all the summary statistics of interest always differ less than $0.1\,\sigma$ from the chosen benchmark, and CLOE predictions are statistically compatible with simulated data obtained from benchmark codes. The same holds for the comparison of correlation function in configuration space for spectroscopic and photometric observables.

Euclid is expected to establish new state-of-the-art constraints on extensions beyond the standard LCDM cosmological model by measuring the positions and shapes of billions of galaxies. Specifically, its goal is to shed light on the nature of dark matter and dark energy. Achieving this requires developing and validating advanced statistical tools and theoretical prediction software capable of testing extensions of the LCDM model. In this work, we describe how the Euclid likelihood pipeline, Cosmology Likelihood for Observables in Euclid (CLOE), has been extended to accommodate alternative cosmological models and to refine the theoretical modelling of Euclid primary probes. In particular, we detail modifications made to CLOE to incorporate the magnification bias term into the spectroscopic two-point correlation function of galaxy clustering. Additionally, we explain the adaptations made to CLOE's implementation of Euclid primary photometric probes to account for massive neutrinos and modified gravity extensions. Finally, we present the validation of these CLOE modifications through dedicated forecasts on synthetic Euclid-like data by sampling the full posterior distribution and comparing with the results of previous literature. In conclusion, we have identified in this work several functionalities with regards to beyond-LCDM modelling that could be further improved within CLOE, and outline potential research directions to enhance pipeline efficiency and flexibility through novel inference and machine learning techniques.

The neglect of modelling both stellar and nebular emission significantly affects the derived physical properties of galaxies, particularly those with high star formation rates. While this issue has been studied, it has not been established a clear threshold for a significant impact on the estimated physical properties of galaxies due to accounting for both stellar and nebular emission. We analyse galaxies from SDSS-DR7 across a wide range of star-forming activity levels, comparing the results obtained from two spectral fitting tools: FADO (which considers both stellar and nebular continuum) and STARLIGHT (only considers the stellar continuum). A strong linear correlation is found between the rest-frame H$\alpha$ and H$\beta$ equivalent widths (EWs) and the optical nebular contribution, identifying these as reliable tracers. The results show that when the nebular contribution exceeds 8% (corresponding to EW(H$\alpha$)$\simeq$500 Å and EW(H$\beta$)$\simeq$110 Å), there is a significant impact on the estimation of galaxy properties, namely stellar mass, age and metallicity. Our results highlight the importance of taking into account both the stellar and nebular continuum when analysing the optical spectra of star-forming galaxies. In particular, this is a fundamental aspect for galaxies with a rest-frame EW(H$\alpha$)$\gtrsim$500 Å (or the scaled value of 375 Å for pseudo-continuum measures). At low redshifts, this mostly impacts extreme emission line galaxies, while at higher redshifts it becomes a dominant aspect given the higher star-forming activity in the younger Universe. In light of current JWST observations and future instruments designed for high-redshift observations, such as MOONS, this reveals as a critical issue to take into consideration.

We present the discovery and characterization of two warm mini-Neptunes transiting the K3V star TOI-815 in a K-M binary system. Analysis of the spectra and rotation period reveal it to be a young star with an age of $200^{+400}_{-200}$Myr. TOI-815b has a 11.2-day period and a radius of 2.94$\pm$0.05$\it{R_{\rm\mathrm{\oplus}}}$ with transits observed by TESS, CHEOPS, ASTEP, and LCOGT. The outer planet, TOI-815c, has a radius of 2.62$\pm$0.10$\it{R_{\rm\mathrm{\oplus}}}$, based on observations of three non-consecutive transits with TESS, while targeted CHEOPS photometry and radial velocity follow-up with ESPRESSO were required to confirm the 35-day period. ESPRESSO confirmed the planetary nature of both planets and measured masses of 7.6$\pm$1.5 $\it{M_{\rm \mathrm{\oplus}}}$ ($\rho_\mathrm{P}$=1.64$^{+0.33}_{-0.31}$gcm$^{-3}$) and 23.5$\pm$2.4$\it{M_{\rm\mathrm{\oplus}}}$ ($\rho_\mathrm{P}$=7.2$^{+1.1}_{-1.0}$gcm$^{-3}$) respectively. Thus, the planets have very different masses, unlike the usual similarity of masses in compact multi-planet systems. Moreover, our statistical analysis of mini-Neptunes orbiting FGK stars suggests that weakly irradiated planets tend to have higher bulk densities compared to those suffering strong irradiation. This could be ascribed to their cooler atmospheres, which are more compressed and denser. Internal structure modeling of TOI-815b suggests it likely has a H-He atmosphere constituting a few percent of the total planet mass, or higher if the planet is assumed to have no water. In contrast, the measured mass and radius of TOI-815c can be explained without invoking any atmosphere, challenging planetary formation theories. Finally, we infer from our measurements that the star is viewed close to pole-on, which implies a spin-orbit misalignment at the 3$\sigma$ level.

We compare the performance of the flat-sky approximation and Limber approximation for the clustering analysis of the photometric galaxy catalogue of Euclid. We study a 6 bin configuration representing the first data release (DR1) and a 13 bin configuration representative of the third and final data release (DR3). We find that the Limber approximation is sufficiently accurate for the analysis of the wide bins of DR1. Contrarily, the 13 bins of DR3 cannot be modelled accurately with the Limber approximation. Instead, the flat-sky approximation is accurate to below $5\%$ in recovering the angular power spectra of galaxy number counts in both cases and can be used to simplify the computation of the full power spectrum in harmonic space for the data analysis of DR3.

Dust attenuation in galaxies has often been used as a proxy for the extinction of point sources, such as supernovae, even though this approach ignores fundamental differences between the two cases. We present an analysis of the impact of geometric effects and scattering within dusty media on recovered galaxy dust properties. We use SKIRT, a radiative transfer code, to simulate observations of point sources embedded in dust clouds, as well as spiral and elliptical galaxies. We examine various galaxy morphologies, inclinations, and instrument apertures. We find that in galaxies the scattering of light into the line of sight and the presence of sources at different depths within the galaxy make attenuation fundamentally different from extinction. For a medium with intrinsic extinction slope Rv=3.068, we recover effective attenuation slopes Rv_e ranging from 0.5 to 7, showing that the two quantities are not analogous, even for local resolved observations. We find that Rv_e greatly depends on dust density, galaxy morphology, and inclination, the latter being the most significant. A single simulated galaxy, viewed from different angles, can reproduce the well-known relation between attenuation strength Av_e and Rv_e observed for star-forming galaxy samples. An increase in dust density leads to higher Rv_e across all inclinations, which, assuming a correlation between stellar mass and dust density, explains the increase in Rv_e with mass observed in star-forming galaxies. However, we are unable to explain the differences in Rv_e between star-forming and quiescent high-mass galaxies. We conclude that highly attenuated regions of simulated face-on galaxies yield Rv_e within 10% of the intrinsic extinction slope of the medium, allowing for the distinction of different dust types. For edge-on spirals, however, the median Rv_e for low Av_e regions appears to better approximate the extinction slope.

In this paper, we explore higher-dimensional asymptotically flat wormhole geometries in the framework of Gauss-Bonnet (GB) gravity and investigate the effects of the GB term, by considering a specific radial-dependent redshift function and by imposing a particular equation of state. This work is motivated by previous assumptions that wormhole solutions were not possible for the $k=1$ and \alpha &lt; 0 case, where $k$ is the sectional curvature of an $(n-2)$-dimensional maximally symmetric space, and $\alpha$ is the Gauss-Bonnet coupling constant. However, we emphasize that this discussion is purely based on a nontrivial assumption that is only valid at the wormhole throat, and cannot be extended to the entire radial-coordinate range. In this work, we provide a counterexample to this claim, and find for the first time specific solutions that satisfy the weak energy condition throughout the entire spacetime, for $k=1$ and \alpha &lt; 0. In addition to this, we also present other wormhole solutions which alleviate the violation of the WEC in the vicinity of the wormhole throat.

Analytical computations in relativistic cosmology can be split into two sets: time evolution relating the initial conditions to the observer's light-cone and light propagation to obtain observables. Cosmological perturbation theory in the FLRW coordinates constitutes an efficient tool for the former task, but the latter is dramatically simpler in light-cone-adapted coordinates that trivialize the light rays towards the observer world-line. Here we point out that time evolution and observable reconstruction can be combined into a single computation that relates directly initial conditions to observables. This is possible if one works uniquely in such light-cone coordinates, thus completely bypassing the FLRW "middle-man" coordinates. We first present in detail these light-cone coordinates, extending and generalizing the presently available material in the literature, and construct a particularly convenient subset for cosmological perturbation theory. We then express the Einstein and energy-momentum conservation equations in these coordinates at the fully non-linear level. This is achieved through a careful 2+1+1 decomposition which leads to relatively compact expressions and provides good control over the geometrical interpretation of the involved quantities. Finally, we consider cosmological perturbation theory to linear order, paying attention to the remaining gauge symmetries and consistently obtaining gauge-invariant equations. Moreover, we show that it is possible to implement statistical homogeneity on stochastic fluctuations, despite the fact that the coordinate system privileges the observer world-line.

We use three domain wall simulations from the radiation era to the late time dark energy domination era based on the PRS algorithm to calculate the energy-momentum tensor components of domain wall networks in an expanding universe. Unequal time correlators in the radiation, matter and cosmological constant epochs are calculated using the scaling regime of each of the simulations. The CMB power spectrum of a network of domain walls is determined. The first ever quantitative constraint for the domain wall surface tension is obtained using a Markov chain Monte Carlo method; an energy scale of domain walls of 0.93 MeV, which is close but below the Zel'dovich bound, is determined.