Hydrodynamical cosmological simulations are increasing their level of realism by considering more physical processes and having greater resolution or larger statistics. However, usually either the statistical power of such simulations or the resolution reached within galaxies are sacrificed. Here, we introduce the NewHorizon project in which we simulate at high resolution a zoom-in region of $\sim(16\,\rm Mpc)^3$ that is larger than a standard zoom-in region around a single halo and is embedded in a larger box. A resolution of up to 34 pc is reached within galaxies; this allows the simulation to capture the multi-phase nature of the interstellar medium and the clumpy nature of the star formation process in galaxies. In this introductory paper, we present several key fundamental properties of galaxies and their black holes, including the galaxy mass function, cosmic star formation rate, galactic metallicities, the Kennicutt-Schmidt relation, the stellar-to-halo mass relation, galaxy sizes, stellar kinematics and morphology, gas content within galaxies and its kinematics, and the black hole mass and spin properties over time. The various scaling relations are broadly reproduced by NewHorizon with some differences with the standard observables. Owing to its exquisite spatial resolution, NewHorizon captures the inefficient process of star formation in galaxies, which evolve over time from being more turbulent, gas rich, and star bursting at high redshift. These high-redshift galaxies are also more compact, and they are more elliptical and clumpier until the level of internal gas turbulence decays enough to allow for the formation of discs. The NewHorizon simulation gives access to a broad range of galaxy formation and evolution physics at low-to-intermediate stellar masses, which is a regime that will become accessible in the near future through surveys such as the LSST.

The InterGalactic Magnetic Field (IGMF), which could permeate the cosmic voids but was never detected so far, is considered a relic of the early Universe. Constraints on its strength $B$ can be derived from its influence on time-delayed very-high-energy photons from Gamma-Ray Bursts (GRBs) in the electromagnetic cascades along their path to the Earth. The present lower limit achieved on its intensity is $10^{-18}\;\mathrm{G}$. In this work, we simulate data from the Cherenkov Telescope Array Observatory (CTAO), accounting for realistic observational constraints, and we apply a joint spectral and temporal fit to characterize the IGMF. GRBs 190114C and 221009A are used as test cases to assess the sensitivity of CTAO. They demonstrate that a broad range of IGMF strengths can be probed with a lower bound as high as $10^{-15}\;\mathrm{G}$. Notably, we show that observations by the CTAO first large telescope, LST-1, already allow us to exclude field strengths up to $3\times 10^{-17}\;\mathrm{G}$.

Contamination-free assessments of the interstellar medium and star formation in quasar host galaxies, particularly based on the far-infrared, offer insights into the role of supermassive black holes in galaxy evolution. Motivated by predictions of quasar heating of dust on both nuclear and galaxy-wide scales, we perform two-component (host galaxy + point source) modeling of high-resolution (~0.1") ALMA observations of the FIR continuum in Band 5 (lambda_rest~500 um) of three highly luminous quasars (L_bol~10^47 erg/s), powered by supermassive black holes having M_BH~10^9 M_sun, at z=2. We include Band 9 (lambda_rest~154 um; 0.06" and 0.3") data at high S/N which places further constraints on the unresolved nuclear component in two cases. To break the degeneracy between quasar and stellar heating, we use CO (J=5-4), observed in Band 5, to gauge the expected contribution of star formation to the infrared luminosity. We find very good agreement between the strength and spatial distribution of the extended continuum component and its prediction based on CO (J=5-4). This is supported by the location of our three quasars along the L_(CO 5-4)-L_(IR, SFR) luminosity relation for inactive star-forming galaxies. As a consequence, there is no evidence for additional continuum emission on extended scales which could be attributed to quasar-heated dust. As expected, the nuclear (i.e., torus) contribution is present and subdominant (12% in Band 9 for one quasar with a typical star-forming host) or non-existent (<8% in Band 9 for the starbursting host). Based on the continuum and CO, the presence of substantial levels of ongoing star formation agrees with previous estimates from unresolved ALMA continuum observations which finds SFRs consistent with star-forming main-sequence galaxies. Therefore, our results do not provide evidence for a quasar-mode feedback, even for the most luminous cases at z=2.

Weakened magnetic braking (WMB) was originally proposed in 2016 to explain anomalously rapid rotation in old field stars observed by the Kepler mission. The proximate cause was suggested to be a transition in magnetic morphology from larger to smaller spatial scales. In a series of papers over the past five years, we have collected spectropolarimetric measurements to constrain the large-scale magnetic fields for a sample of stars spanning this transition, including a range of spectral types from late F to early K. During this time, we gradually improved our methods for estimating the wind braking torque in each of our targets, and for evaluating the associated uncertainties. Here, we reanalyze the entire sample with a focus on uniformity for the relevant observational inputs. We supplement the sample with two additional active stars to provide more context for the evolution of wind braking torque with stellar Rossby number (Ro). The results demonstrate unambiguously that standard spin-down models can reproduce the evolution of wind braking torque for active stars, but WMB is required to explain the subsequent abrupt decrease in torque as Ro approaches a critical value for dynamo excitation. This transition is seen in both the large-scale magnetic field and the X-ray luminosity, indicating weakened coronal heating. We interpret these transitions as evidence of a rotational threshold for the influence of Coriolis forces on global convective patterns and the resulting inefficiency of the global stellar dynamo.

We present a statistical study of spatially resolved chemical enrichment in 18 main-sequence galaxies at $z=4$--6, observed with \jwst/NIRSpec IFU as part of the ALPINE-CRISTAL-\jwst\ survey. Performing an optimized reduction and calibration procedure, including local background subtraction, light-leakage masking, stripe removal, and astrometry refinement, we achieve robust emission-line mapping on kiloparsec scales. Although line-ratio distributions vary across galaxies in our sample, we generally find mild central enhancements in [O\,\textsc{iii}]/H$\beta$, [O\,\textsc{ii}]/[O\,\textsc{iii}], [S\,\textsc{ii}]$_{6732}$/[S\,\textsc{ii}]$_{6718}$, H$\alpha$/H$\beta$, and $L_{\rm H\alpha}/L_{\rm UV}$, consistent with elevated electron density, dust obscuration, and bursty star formation accompanied by reduced metallicity and ionization parameter. These features point to inside-out growth fueled by recent inflows of pristine gas. Nevertheless, the median metallicity gradient is nearly flat over a few kpc scale, $\Delta \log({\rm O/H}) = 0.02 \pm 0.01$ dex kpc$^{-1}$, implying efficient chemical mixing through inflows, outflows, and mergers. From pixel-by-pixel stellar and emission-line characterizations, we further investigate the resolved Fundamental Metallicity Relation (rFMR). Metallicity is described by a fundamental plane with stellar mass and SFR surface densities, but with a stronger dependence on $\Sigma_{\rm SFR}$ than seen in local galaxies. Our results indicate that the regulatory processes linking star formation, gas flows, and metal enrichment were already vigorous $\sim$1 Gyr after the Big Bang, producing the nearly flat metallicity gradient and a stronger coupling between star formation and metallicity than observed in evolved systems in the local universe.

Protoclusters, galaxy clusters' high redshift progenitors, hold the keys to understanding the formation and evolution of clusters and their member galaxies. However, their cosmological distances and spatial extensions (tens of Mpc) have inhibited complete mapping of their structure and constituent galaxies, which is key to robustly linking protoclusters to their descendants. Here we report the discovery of the Bigfoot, a tridimensional structure at $z = 3.98$ including 11 subgroups traced by 55 (700) spectroscopic (photometric) redshifts with JWST, extending over $15\times 37$ $\times 49{\rm{cMpc^3}}$ in the PRIMER-UDS field. Bigfoot's large-scale and mass function of member galaxies closely match constrained simulations' predictions for the progenitors of today's most massive clusters (${M_0} > 10^{15} {M_{_ \odot }}$). All subgroups with ${M_{\rm{h}}} > {10^{12.5}}{M_{_ \odot }}$ exhibit enhanced fractions of massive galaxies ($>{10^{10.0} {M_{_ \odot }}}$) compared to lower-mass halos and the field, demonstrating the accelerated formation of massive galaxies in massive halos. The presence of this massive protocluster with a large central halo (${10^{13.0} {M_{_ \odot }}}$) in a JWST deep field bears important cosmological implication that favors high ${\sigma _8}$ of PLANCK cosmology over low-redshift probes.

Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster medium. We estimate the sensitivity of the Cherenkov Telescope Array (CTA) to detect diffuse gamma-ray emission from the Perseus galaxy cluster. We perform a detailed spatial and spectral modelling of the expected signal for the DM and the CRp components. For each, we compute the expected CTA sensitivity. The observing strategy of Perseus is also discussed. In the absence of a diffuse signal (non-detection), CTA should constrain the CRp to thermal energy ratio within the radius $R_{500}$ down to about $X_{500}<3\times 10^{-3}$, for a spatial CRp distribution that follows the thermal gas and a CRp spectral index $\alpha_{\rm CRp}=2.3$. Under the optimistic assumption of a pure hadronic origin of the Perseus radio mini-halo and depending on the assumed magnetic field profile, CTA should measure $\alpha_{\rm CRp}$ down to about $\Delta\alpha_{\rm CRp}\simeq 0.1$ and the CRp spatial distribution with 10% precision. Regarding DM, CTA should improve the current ground-based gamma-ray DM limits from clusters observations on the velocity-averaged annihilation cross-section by a factor of up to $\sim 5$, depending on the modelling of DM halo substructure. In the case of decay of DM particles, CTA will explore a new region of the parameter space, reaching models with \tau_{\chi}&gt;10^{27}s for DM masses above 1 TeV. These constraints will provide unprecedented sensitivity to the physics of both CRp acceleration and transport at cluster scale and to TeV DM particle models, especially in the decay scenario.

For the first time in nearly a decade, a new, bright transient was detected in the central parsec (pc) of the Galaxy. MAXI J1744-294 was never observed in outburst prior to January 2025. We present the results of a broadband, multi-wavelength study of this enigmatic source, including data from the NuSTAR, Chandra, XMM-Newton, Swift, and NICER X-ray telescopes, as well as complementary radio and near-infrared observations. We find that MAXI J1744-294 remained in the bright/soft state throughout the first months of 2025. Spectral hardening was observed in April 2025, followed by a decline in flux. Based on the spectral and temporal characteristics of the source, we identify MAXI J1744-294 as a candidate black hole (BH) low-mass X-ray binary (LMXB) $-$ the fourth candidate BH transient discovered within a (projected) distance of one pc from the Galactic supermassive black hole Sgr A*. This discovery provides further evidence for a cusp of BH-LMXBs in the central pc of our Galaxy, as argued for in previous observational work and suggested by analytical and theoretical work. Our multi wavelength study, involving a complementary range of observatories and spanning different outburst states, can serve as a model for future time domain astrophysics research.

H2 is the most abundant molecule in the interstellar medium and is a useful tool to study photodissociation regions, where radiative feedback from massive stars on molecular clouds is dominant. The James Webb Space Telescope, with its high spatial resolution, sensitivity, and wavelength coverage provides unique access to the detection of most of H2 lines and the analysis of its spatial morphology. Our goal is to use H2 line emission detected with the JWST in the Horsehead nebula to constrain the physical parameters (e.g., extinction, gas temperature, thermal pressure) throughout the PDR and its geometry. The study of H2 morphology reveals that FUV-pumped lines peak closer to the edge of the PDR than thermalized lines. From H2 lines, we estimate the value of extinction throughout the PDR. We find that AV is increasing from the edge of the PDR to the second and third H2 filaments. Then, we study the H2 excitation in different regions across the PDR. The temperature profile shows that the observed gas temperature is quite constant throughout the PDR, with a slight decline in each of the dissociation fronts. This study also reveals that the OPR is far from equilibrium. We observe a spatial separation of para and ortho rovibrational levels, indicating that efficient ortho-para conversion and preferential ortho self-shielding are driving the spatial variations of the OPR. Finally, we derive a thermal pressure in the first filament around P > 6x10$^6$ K cm$^{-3}$, about ten times higher than that of the ionized gas. We highlight that template stationary 1D PDR models cannot account for the intrinsic 2D structure and the very high temperature observed in the Horsehead nebula. We argue the highly excited, over-pressurized H2 gas at the edge of the PDR interface could originate from the mixing between the cold and hot phase induced by the photo-evaporation of the cloud.

This white paper explores the potential for strategic synergies between the JWST and the Ariel telescopes, two flagship observatories poised to revolutionise the study of exoplanet atmospheres. Both telescopes have the potential to address common fundamental questions about exoplanets-especially concerning their nature and origins-and serve a growing scientific community. With their operations now anticipated to overlap, starting from 2030, there is a unique opportunity to enhance the scientific outputs of both observatories through coordinated efforts. In this report, authored by the Ariel-JWST Synergy Working Group, part of the Ariel Consortium Science Team, we summarise the capabilities of JWST and Ariel; we highlight their key differences, similarities, synergies, and distinctive strengths. Ariel is designed to conduct a broad survey of exoplanet atmospheres but remains highly flexible, allowing the mission to integrate insights from JWST's discoveries. Findings from JWST, including data from initiatives shaped by NASA's decadal survey priorities and community-driven research themes, will inform the development of Ariel's core survey strategy. Conversely, Ariel's ability to perform broad-wavelength coverage observations for bright targets provides complementary avenues for exoplanet researchers, particularly those interested in time-domain observations and large-scale atmospheric studies. This paper identifies key pathways for fostering JWST-Ariel synergies, many of which can be initiated even before Ariel's launch. Leveraging their complementary designs and scopes, JWST and Ariel can jointly address fundamental questions about the nature, formation, and evolution of exoplanets. Such strategic collaboration has the potential to maximise the scientific returns of both observatories and lay the foundation for future facilities in the roadmap to exoplanet exploration.

We report on coronagraphic observations of the nearest solar-type star, $\alpha$ Cen A, using the MIRI instrument on the James Webb Space Telescope. With three epochs of observation (August 2024, February 2025, and April 2025), we achieve a sensitivity sufficient to detect $T_{\rm eff}\approx$ 225-250 K (1-1.2 $R_{\rm Jup}$) planets between 1"-2" and exozodiacal dust emission at the level of $>$5-8$\times$ the brightness of our own zodiacal cloud. The lack of exozodiacal dust emission sets an unprecedented limit of a few times the brightness of our own zodiacal cloud$-$a factor of $\gtrsim$10 more sensitive than measured toward any other stellar system to date. In August 2024, we detected a F$_\nu$(15.5 $\mu$m) = 3.5 mJy point source, called $S1$, at a separation of 1.5" from $\alpha$ Cen A. Because the August 2024 epoch had only one successful observation at a single roll angle, it is not possible to unambiguously confirm $S1$ as a bona fide planet. Our analysis confirms that $S1$ is neither a background nor a foreground object. $S1$ is not recovered in the February and April 2025 epochs. However, if $S1$ is the counterpart of the object, $C1$, seen by the VLT/NEAR program in 2019, we find that there is a 52% chance that the $S1+C1$ candidate was missed in both follow-up JWST/MIRI observations due to orbital motion. Incorporating constraints from the non-detections, we obtain families of dynamically stable orbits for $S1+C1$ with periods between 2-3 years. These suggest that the planet candidate is on an eccentric ($e \approx 0.4$) orbit significantly inclined with respect to $\alpha$ Cen AB orbital plane ($i_{\rm mutual} \approx 50^\circ$, or $\approx 130^\circ$). Based on the photometry and orbital properties, the planet candidate could have a temperature of 225 K, a radius of $\approx$1-1.1 $R_{\rm Jup}$ and a mass between 90-150 $M_{\rm Earth}$, consistent with RV limits.

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.

Cross-correlations techniques offer an alternative method to search for molecular species in JWST observations of exoplanet atmospheres. In a previous article, we applied cross-correlation functions for the first time to JWST NIRSpec/G395H observations of exoplanet atmospheres, resulting in a detection of CO in the transmission spectrum of WASP-39b and a tentative detection of CO isotopologues. Here we present an improved version of our cross-correlation technique and an investigation into how efficient the technique is when searching for other molecules in JWST NIRSpec/G395H data. Our search results in the detection of more molecules via cross-correlations in the atmosphere of WASP-39b, including $\rm H_{2}O$ and $\rm CO_{2}$, and confirms the CO detection. This result proves that cross-correlations are a robust and computationally cheap alternative method to search for molecular species in transmission spectra observed with JWST. We also searched for other molecules ($\rm CH_{4}$, $\rm NH_{3}$, $\rm SO_{2}$, $\rm N_{2}O$, $\rm H_{2}S$, $\rm PH_{3}$, $\rm O_{3}$ and $\rm C_{2}H_{2}$) that were not detected, for which we provide the definition of their cross-correlation baselines for future searches of those molecules in other targets. We find that that the cross-correlation search of each molecule is more efficient over limited wavelength regions of the spectrum, where the signal for that molecule dominates over other molecules, than over broad wavelength ranges. In general we also find that Gaussian normalization is the most efficient normalization mode for the generation of the molecular templates.

Based on the DUSTGRAIN-pathfinder suite of simulations, we investigate observational degeneracies between nine models of modified gravity and massive neutrinos. Three types of machine learning techniques are tested for their ability to discriminate lensing convergence maps by extracting dimensional reduced representations of the data. Classical map descriptors such as the power spectrum, peak counts and Minkowski functionals are combined into a joint feature vector and compared to the descriptors and statistics that are common to the field of digital image processing. To learn new features directly from the data we use a Convolutional Neural Network (CNN). For the mapping between feature vectors and the predictions of their underlying model, we implement two different classifiers; one based on a nearest-neighbour search and one that is based on a fully connected neural network. We find that the neural network provides a much more robust classification than the nearest-neighbour approach and that the CNN provides the most discriminating representation of the data. It achieves the cleanest separation between the different models and the highest classification success rate of 59% for a single source redshift. Once we perform a tomographic CNN analysis, the total classification accuracy increases significantly to 76% with no observational degeneracies remaining. Visualising the filter responses of the CNN at different network depths provides us with the unique opportunity to learn from very complex models and to understand better why they perform so well.

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.