KICC papers

arXiv:2409.11464v3 Announce Type: replace Abstract: For the first time, we present a systematic search for galaxies with extended emission line and potential outflow features using \textit{JWST} medium-band images in the GOODS-S field. This is done by comparing the morphology in medium-band images to adjacent continuum and UV bands. We look for galaxies that have a maximum extent 50\% larger, an excess area 30\% greater, or an axis ratio difference of more than 0.3 in the medium band compared to the reference bands. After visual inspection, we find 326 candidate galaxies at $1.4 < z < 8.4$, with a peak in the population near cosmic noon, benefiting from the good coverage of the medium-band filters. By fitting their SEDs, we find that the candidate galaxies are at least 20\% more bursty in their star-forming activity and have 50\% more young stellar populations compared to a control sample selected based on the continuum band flux. Additionally, these candidates exhibit a significantly higher production rate of ionizing photons. We further find that candidates hosting known AGN produce extended emission that is more anisotropic compared to non-AGN candidates. A few of our candidates have been spectroscopically confirmed to have prominent outflow signatures through NIRSpec observations, showcasing the robustness of the photometric selection. Future spectroscopic follow-up will better help verify and characterize the kinematics and chemical properties of these systems.

arXiv:2505.01512v1 Announce Type: new Abstract: The discovery of inhabited exoplanets hinges on identifying biosignature gases. JWST is revealing potential biosignatures in exoplanet atmospheres, though their presence is yet to provide strong evidence for life. The central challenge is attribution: how to confidently identify biogenic sources while ruling out, or deeming unlikely, abiotic explanations? Attribution is particularly difficult for individual planets, especially regarding system-scale stochastic processes that could set atmospheric conditions. To address this, we here propose a comparative multi-planet approach: comparing atmospheric compositions across multiple planets within a system to empirically define the 'abiotic baseline'. This baseline serves as a reference point for biosignatures, and enables marginalisation over inaccessible, shared abiotic parameters. This is possible because planets within a system are linked by their birth in the same natal disk, having been irradiated by the same evolving star, and having a related dynamical history. Observations aligning with the abiotic baseline, where the locally informed abiotic models demonstrate high out-of-sample predictive accuracy, are likely non-biological. Deviations from the baseline -- potentially biotic anomalies -- suggest an alternative origin. We present the application of Bayesian leave-one-out cross-validation to evaluate the performance of geochemical- and biogeochemical-climate models in explaining these anomalies, using the expected log pointwise predictive density as a diagnostic. When biogeochemical models outperform their abiotic counterparts, the anomaly may be shaped by life, and constitutes a comparative biosignature. If both models perform poorly, the anomaly is flagged as an "unknown unknown" -- a signature of either unrecognised abiotic chemistry, or life as we don't yet know it.

arXiv:2410.19905v2 Announce Type: replace Abstract: Energetic feedback processes associated with accreting supermassive black holes can expel gas from massive haloes and significantly alter various measures of clustering on ~Mpc scales, potentially biasing the values of cosmological parameters inferred from analyses of large-scale structure (LSS) if not modelled accurately. Here we use the state-of-the-art FLAMINGO suite of cosmological hydrodynamical simulations to gauge the impact of feedback on large-scale structure by comparing to Planck + ACT stacking measurements of the kinetic Sunyaev-Zel'dovich (kSZ) effect of SDSS BOSS galaxies. We make careful like-with-like comparisons to the observations, aided by high precision KiDS and DES galaxy-galaxy lensing measurements of the BOSS galaxies to inform the selection of the simulated galaxies. In qualitative agreement with several recent studies using dark matter only simulations corrected for baryonic effects, we find that the kSZ effect measurements prefer stronger feedback than predicted by simulations which have been calibrated to reproduce the gas fractions of low redshift X-ray-selected groups and clusters. We find that the increased feedback can help to reduce the so-called S8 tension between the observed and CMB-predicted clustering on small scales as probed by cosmic shear (although at the expense of agreement with the X-ray group measurements). However, the increased feedback is only marginally effective at reducing the reported offsets between the predicted and observed clustering as probed by the thermal SZ (tSZ) effect power spectrum and tSZ effect--weak lensing cross-spectrum, both of which are sensitive to higher halo masses than cosmic shear.

arXiv:2505.00095v1 Announce Type: new Abstract: Metallicity ($Z$) estimates based on ultraviolet (UV) emission lines from the narrow-line regions (NLRs) of active galactic nuclei (AGNs) have been found to differ from those derived from optical lines. However, the origin of this discrepancy ($ZR$) remains poorly understood. To investigate the source of $ZR$, we compiled from the literature the fluxes of narrow near-UV ($1000 < \lambda(\angstrom) < 2000)$ and optical ($3000 < \lambda(\angstrom) < 7000)$ emission line measurements for a sample of 11 AGNs (9 at $z<0.4$ and 2 at $z\sim2.4$). Metallicity values for our sample were derived using a semi-empirical calibration based on the $C43$=log[(\ion{C{iv}$\lambda$1549+\ion{C{iii}]$\lambda$1909)/\ion{He}{ii}$\lambda$1640] emission-line ratio and compared with those obtained via direct measurement of the electron temperature ($T_{\rm e}$-method) and via calibrations based on optical emission-lines. The source of the discrepancy was investigated in terms of the ionization parameter ($U$), electron density ($N_{\rm e}$), and carbon abundance (C/H). We found a weak correlation between $ZR$, $U$ and $N_{\rm e}$. However, a moderate correlation was observed between $ZR$ and direct estimates of C/H, suggesting that the previously assumed (C/O)-$Z$ relations in photoionization models used to derive UV carbon-line calibrations may not be valid for AGNs. By combining a large set of abundance estimates for local star-forming regions with those of our AGN sample, we derived a new (C/O)-$Z$ relation. Comparisons between the results of photoionization models that assume this new abundance relation and the UV observational data of our sample produce $Z$ values derived from the $C43$ index that are consistent with those obtained using the $T_{\rm e}$-method.

arXiv:2409.20379v2 Announce Type: replace Abstract: A precise measurement of photometric redshifts (photo-z) is key for the success of modern photometric galaxy surveys. Machine learning (ML) methods show great promise in this context, but suffer from covariate shift (CS) in training sets due to selection bias where interesting sources are underrepresented, and the corresponding ML models show poor generalisation properties. We present an application of the StratLearn method to the estimation of photo-z, validating against simulations where we enforce the presence of CS to different degrees. StratLearn is a statistically principled approach that relies on splitting the source and target datasets into strata based on estimated propensity scores (i.e. the probability for an object to be in the source set given its observed covariates). After stratification, two conditional density estimators are fit separately to each stratum, then combined via a weighted average. We benchmark our results against the GPz algorithm, quantifying the performance of the two codes with a set of metrics. Our results show that the StratLearn-z metrics are only marginally affected by the presence of CS, while GPz shows a significant degradation of performance in the photo-z prediction for fainter objects. For the strongest CS scenario, StratLearn-z yields a reduced fraction of catastrophic errors, a factor of 2 improvement for the RMSE and one order of magnitude improvement on the bias. We also assess the quality of the conditional redshift estimates with the probability integral transform (PIT). The PIT distribution obtained from StratLearn-z features fat fewer outliers and is symmetric, i.e. the predictions appear to be centered around the true redshift value, despite showing a conservative estimation of the spread of the conditional redshift distributions. Our julia implementation of the method is available at https://github.com/chiaramoretti/StratLearn-z.

arXiv:2504.21087v1 Announce Type: new Abstract: The Einstein Telescope (ET), along with other third-generation gravitational wave (GW) detectors, will be a key instrument for detecting GWs in the coming decades. However, analyzing the data and estimating source parameters will be challenging, especially given the large number of expected detections - of order $10^5$ per year - which makes current methods based on stochastic sampling impractical. In this work, we use Dingo-IS to perform Neural Posterior Estimation (NPE) of high-redshift events detectable with ET in its triangular configuration. NPE is a likelihood-free inference technique that leverages normalizing flows to approximate posterior distributions. After training, inference is fast, requiring only a few minutes per source, and accurate, as corrected through importance sampling and validated against standard Bayesian inference methods. To confirm previous findings on the ability to estimate parameters for high-redshift sources with ET, we compare NPE results with predictions from the Fisher information matrix (FIM) approximation. We find that FIM underestimates sky localization errors substantially for most sources, as it does not capture the multimodalities in sky localization introduced by the geometry of the triangular detector. FIM also overestimates the uncertainty in luminosity distance by a factor of $\sim 3$ on average when the injected luminosity distance is $d^{\mathrm{inj}}_{\mathrm{L}} > 10^5~$Mpc, further confirming that ET will be particularly well suited for studying the early Universe.

arXiv:2504.20144v1 Announce Type: cross Abstract: We report the existence of two exotic compact objects in the leading relativistic model of modified Newtonian dynamics, namely aether-scalar-tensor theory. This model is consistent with precision cosmology and gravitational wave constraints on tensor speed. Black hole mimickers could subtly change observations: gravitational waves from their mergers might show unusual echoes or altered ringdown patterns, and images of their horizon-scale shadows might be slightly different from those of a true black hole. Shapiro-free lenses are massless objects that deflect light without any gravitational time delay, producing distinctive lensing events. These predictions connect to ongoing and future gravitational-wave searches, horizon-scale imaging, and time-domain lensing surveys.

arXiv:2503.17342v2 Announce Type: replace Abstract: We present an updated reconstruction of the dark energy equation of state, $w(a)$, using the newly released DESI DR2 Baryon Acoustic Oscillation (BAO) data in combination with Pantheon+ and DES5Y Type Ia supernovae measurements, respectively. Building on our previous analysis in arXiv:2503.08658, which employed a nonparametric flexknot reconstruction approach, we examine whether the evidence for dynamical dark energy persists with the improved precision of the DESI DR2 dataset. We find that while the overall qualitative structure of $w(a)$ remains consistent with our earlier findings, the statistical support for dynamical dark energy is reduced when considering DESI DR2 data alone, particularly for more complex flexknot models with higher numbers of knots. However, the evidence for simpler dynamical models, such as $w$CDM and CPL (which correspond to $n=1$ and $n=2$ knots respectively), increases relative to $\Lambda$CDM with DESI DR2 alone, with CPL being the preferred dynamical model, consistent with previous DESI analyses. When combined with Pantheon+ data, the conclusions remain broadly consistent with our earlier work, but when instead combined with DES5Y supernovae data, there is an increased preference for flexknot models for all values of $n$ considered. This results in all such models being preferred over $\Lambda$CDM, with the CPL model being the most favoured by a Bayes factor of $\sim 2.3$ relative to $\Lambda$CDM.

arXiv:2504.20992v1 Announce Type: new Abstract: We analyse the large-scale angular clustering of quasars in the Gaia-unWISE quasar catalog, Quaia, and their cross-correlation with maps of the lensing convergence of the Cosmic Microwave Background (CMB), to constrain the level of primordial non-Gaussianity (PNG). Specifically, we target the scale-dependent bias that would be induced by PNG on biased tracers of the matter inhomogeneities on large scales. The Quaia sample is particularly well suited for this analysis, given the large effective volume covered, and our ability to map out the main potential sources of systematic contamination and mitigate their impact. Using the universality relation to characterise the response of the quasar overdensity to PNG ($p_\phi=1$), we report constraints on the local-type PNG parameter $f_{\rm NL}$ of $f_{\rm NL}=-20.5^{+19.0}_{-18.1}$ (68\% C.L.) by combining the quasar auto-correlation and its cross-correlation with CMB lensing in two tomographic redshift bins (or $f_{\rm NL}=-28.7^{+26.1}_{-24.6}$ if assuming a lower response for quasars, $p_\phi=1.6$). Using the CMB lensing cross-correlations alone, we find $f_{\rm NL}=-13.8^{+26.7}_{-25.0}$. These are the tightest constraints on $f_{\rm NL}$ to date from angular clustering statistics and cross-correlations with CMB lensing.

arXiv:2504.20038v1 Announce Type: new Abstract: We present the tightest cosmic microwave background (CMB) lensing constraints to date on the growth of structure by combining CMB lensing measurements from the Atacama Cosmology Telescope (ACT), the South Pole Telescope (SPT) and \textit{Planck}. Each of these surveys individually provides lensing measurements with similarly high statistical power, achieving signal-to-noise ratios of approximately 40. The combined lensing bandpowers represent the most precise CMB lensing power spectrum measurement to date with a signal-to-noise ratio of 61 and an amplitude of $A_\mathrm{lens}^\mathrm{recon} = 1.025 \pm 0.017$ with respect to the theory prediction from the best-fit CMB \textit{Planck}-ACT cosmology. The bandpowers from all three lensing datasets, analyzed jointly, yield a $1.6\%$ measurement of the parameter combination $S_8^\mathrm{CMBL} \equiv \sigma_8\,(\Omega_m/0.3)^{0.25} = 0.825^{+0.015}_{-0.013}$. Including Dark Energy Spectroscopic Instrument (DESI) Baryon Acoustic Oscillation (BAO) data improves the constraint on the amplitude of matter fluctuations to $\sigma_8 = 0.829 \pm 0.009$ (a $1.1\%$ determination). When combining with uncalibrated supernovae from \texttt{Pantheon+}, we present a $4\%$ sound-horizon-independent estimate of $H_0=66.4\pm2.5\,\mathrm{km\,s^{-1}\,Mpc^{-1}} $. The joint lensing constraints on structure growth and present-day Hubble rate are fully consistent with a $\Lambda$CDM model fit to the primary CMB data from \textit{Planck} and ACT. While the precise upper limit is sensitive to the choice of data and underlying model assumptions, when varying the neutrino mass sum within the $\Lambda\mathrm{CDM}$ cosmological model, the combination of primary CMB, BAO and CMB lensing drives the probable upper limit for the mass sum towards lower values, comparable to the minimum mass prior required by neutrino oscillation experiments.

arXiv:2412.05363v3 Announce Type: replace Abstract: Recent observations by JWST have revealed supersolar $^{14}$N abundances in galaxies at very high redshift. On the other hand, these galaxies show subsolar metallicity. The observed N/O ratios are difficult to reproduce in the framework of chemical evolution models for the Milky Way. Our aim is to reproduce these high N/O ratios with chemical evolution models assuming different histories of star formation triggering galactic winds coupled with detailed nucleosynthesis prescriptions for $^{14}$N, $^{12}$C, $^{16}$O and $^{56}$Fe. We compute several models for small galaxies ($10^{9}\text{ - }10^{10}\text{ M}_{\odot}$) with high star formation efficiency and strong galactic winds. These winds are assumed to be differential, carrying out mainly the products of the explosion of core-collapse supernovae. We find that only models with high star formation rates, normal initial mass function, and differential galactic winds can reproduce the observed chemical abundances. We also find that with the same assumptions about star formation and galactic winds, but with a very rapid formation resulting from fast gas infall, we can also reproduce the estimated ages of these objects. We find no necessity to invoke peculiar nucleosynthesis from Population III stars, very massive stars and supermassive stars.

arXiv:2504.18637v1 Announce Type: new Abstract: The Sunyaev-Zeldovich (SZ) effect is a window into the astrophysical processes of galaxy clusters, and relativistic corrections (the "rSZ") promise to provide a global census of the gas feedback within clusters. Upcoming wide-field millimeter-wave surveys such as the Simons Observatory (SO), Fred Young Submillimeter Telescope, and CMB-S4 will make increasingly precise measurements of the SZ effect and its relativistic corrections. We present simulated full-sky maps of the rSZ effect and a fast code to generate it, for use in the development of analysis techniques and pipelines. As part of the websky simulation suite, our mock observations have semi-realistic cross-correlations with other large-scale structure tracers, offering insights into the formation and evolution of galaxy clusters and large-scale structure. As a demonstration of this, we examine what an SO-like experiment can learn from the rSZ effect. We find that high significance detections will be possible, provided that the instrumental systematics are under control, and that the evolution of cluster temperatures with mass and redshift can be probed in a manner complementary to X-ray measurements.

arXiv:2504.18384v1 Announce Type: new Abstract: While mounting observational evidence suggests that intermediate mass black holes (IMBHs) may be important in shaping the properties of dwarf galaxies both at high redshifts and in the local Universe, our theoretical understanding of how these IMBHs grow is largely incomplete. To address this, we perform high-resolution simulations of an isolated dwarf galaxy with a virial mass of $10^{10}~{\rm M}_{\odot}$ harbouring a $10^4~{\rm M}_{\odot}$ IMBH at its centre at a peak spatial resolution of $\lesssim 0.01$ pc. Within the fully multi-phase interstellar medium (ISM), we incorporate explicit sampling of stars from the initial mass function, photo-ionization, photoelectric heating, individual supernovae (SNe), as well as a Shakura-Sunyaev accretion disc model to track the evolution of BH mass and spin. We find that a nuclear star cluster (NSC) effectively captures the ISM gas and promotes formation of a circumnuclear disc (CND) on scales of $\lesssim 7$ pc. Simultaneously, gravitational torques from the NSC reduce CND angular momentum on (sub-)parsec scales, circularizing the gas onto the $\alpha$-accretion disc and promoting sustained IMBH growth at $\sim 0.01$ of the Eddington rate. While in the innermost regions ($\lesssim 0.5$ pc), star formation is highly suppressed, the CND is susceptible to fragmentation, leading to the formation of massive, young stars. Interestingly, despite an in-situ SN rate of $0.3~{\rm Myr}^{-1}$, the dense CND persists, sustaining BH accretion and leading to its net spin-up. Our study demonstrates the complexity of IMBH accretion within a multi-phase ISM, and paves the way for next-generation studies where IMBH growth in a fully cosmological context can be captured.

arXiv:2504.18245v1 Announce Type: new Abstract: The multi-tracer (MT) technique has been shown to outperform single-tracer analyses in the context of galaxy clustering. In this paper, we conduct a series of Fisher analyses to further explore MT information gains within the framework of non-linear bias expansion. We examine how MT performance depends on the bias parameters of the subtracers, showing that directly splitting the non-linear bias generally leads to smaller error bars in $A_s$, $h$, and $\omega_{\rm cdm}$ compared to a simple split in $b_1$. This finding opens the door to identifying subsample splits that do not necessarily rely on very distinct linear biases. We discuss different total and subtracer number density scenarios, as well as the possibility of splitting into more than two tracers. Additionally, we consider how different Fingers-of-God suppression scales for the subsamples can be translated into different $k_{\rm max}$ values. Finally, we present forecasts for ongoing and future galaxy surveys.

arXiv:2407.18349v2 Announce Type: replace Abstract: Understanding the impact of baryonic physics on cosmic structure formation is crucial for accurate cosmological predictions, especially as we usher in the era of large galaxy surveys with the Rubin Observatory as well as the Euclid and Roman Space Telescopes. A key process that can redistribute matter across a large range of scales is feedback from accreting supermassive black holes. How exactly these active galactic nuclei (AGN) operate from sub-parsec to Mega-parsec scales however remains largely unknown. To understand this, we investigate how different AGN feedback models in the Fable simulation suite affect the cosmic evolution of the matter power spectrum (MPS). Our analysis reveals that AGN feedback significantly suppresses clustering at scales $k \sim 10\,h\,cMpc^{-1}$, with the strongest effect at redshift $z = 0$ causing a reduction of $\sim 10\%$ with respect to the dark matter-only simulation. This is due to the efficient feedback in both radio (low Eddington ratio) and quasar (high Eddington ratio) modes in our fiducial Fable model. We find that variations of the quasar and radio mode feedback with respect to the fiducial Fable model have distinct effects on the MPS redshift evolution, with the radio mode being more effective on larger scales and later epochs. Furthermore, MPS suppression is dominated by AGN feedback effects inside haloes at $z = 0$, while for $z \gtrsim 1$ the matter distribution both inside and outside of haloes shapes the MPS suppression. Hence, future observations probing earlier cosmic times beyond $z \sim 1$ will be instrumental in constraining the nature of AGN feedback.

arXiv:2504.17549v1 Announce Type: new Abstract: The co-evolution of massive black holes (BHs) and their host galaxies is well-established within the hierarchical galaxy formation paradigm. Large-scale cosmological simulations are an ideal tool to study the repeated BH mergers, accretion and feedback that conspire to regulate this process. While such simulations are of fundamental importance for understanding the complex and intertwined relationship between BHs and their hosts, they are plagued with numerical inaccuracies at the scale of individual BH orbits. To quantify this issue, taking advantage of the $(100 \, h^{-1}\,\text{cMpc})^3$ FABLE simulation box, we track all individual BH mergers and the corresponding host galaxy mergers as a function of cosmic time. We demonstrate that BH mergers frequently occur prematurely, well before the corresponding merger of the host galaxies is complete, and that BHs are sometimes erroneously displaced from their hosts during close galaxy encounters. Correcting for these artefacts results in substantial macrophysical delays, spanning over several Gyrs, which are additional to any microphysical delays arising from unresolved BH binary hardening processes. We find that once the macrophysical delays are accounted for, high-mass BH merger events are suppressed, affecting the predictions for the BH population that may be observable with LISA and pulsar timing arrays. Furthermore, including these macrophysical delays leads to an increase in the number of observable dual active galactic nuclei, especially at lower redshifts, with respect to FABLE. Our results highlight the pressing need for more accurate modelling of BH dynamics in cosmological simulations of galaxy formation as we prepare for the multi-messenger era.

arXiv:2504.16988v1 Announce Type: cross Abstract: Precision observations of orbital systems have recently emerged as a promising new means of detecting gravitational waves and ultra-light dark matter, offering sensitivity in new regimes with significant discovery potential. These searches rely critically on precise modeling of the dynamical effects of these signals on the observed system; however, previous analyses have mainly only relied on the secularly-averaged part of the response. We introduce here a fundamentally different approach that allows for a fully time-resolved description of the effects of oscillatory metric perturbations on orbital dynamics. We find that gravitational waves and ultra-light dark matter can induce large oscillations in the orbital parameters of realistic binaries, enhancing the sensitivity to such signals by orders of magnitude compared to previous estimates.

arXiv:2504.17007v1 Announce Type: new Abstract: The James Webb Space Telescope (JWST) discovered 79 transients out to $z$$\sim$4.8 through the JADES Transient Survey (JTS), but the JTS did not find any $z$$>$5 transients. Here, we present the first photometric evidence of a $z$$>$5 transient/variable source with JWST. The source, AT 2023adya, resides in a $z_{\mathrm{spec}}$$=$5.274 galaxy in GOODS-N, which dimmed from $m_{\rm F356W}$$=$26.05$\pm$0.02 mag to 26.24$\pm$0.02 mag in the rest-frame optical over approximately two rest-frame months, producing a clear residual signal in the difference image ($m_{\rm F356W}$$=$28.01$\pm$0.17 mag; SN$_\mathrm{var}$$=$6.09) at the galaxy center. Shorter-wavelength bands (F090W/F115W) show no rest-frame ultraviolet brightness change. Based on its rest-frame V-band absolute magnitude of M$_\mathrm{V}$$=$$-$18.48 mag, AT 2023adya could be any core-collapse supernova (SN) subtype or an SN Ia. However, due to low SN Ia rates at high redshift, the SN Ia scenario is unlikely. Alternatively, AT 2023adya may be a variable active galactic nucleus (AGN). However, the JWST NIRCam/Grism spectrum shows no broad H$\alpha$ emission line (FWHM$=$130$\pm$26 km s$^{-1}$), disfavoring the variable AGN scenario. It is also unlikely that AT 2023adya is a tidal disruption event (TDE) because the TDE models matching the observed brightness changes have low event rates. Although it is not possible to determine AT 2023adya's nature based on the two-epoch single-band photometry alone, this discovery indicates that JWST can push the frontier of transient/variable science past $z$$=$5 and towards the epoch of reionization.

arXiv:2504.16791v1 Announce Type: new Abstract: Radiometers are crucial instruments in radio astronomy, forming the primary component of nearly all radio telescopes. They measure the intensity of electromagnetic radiation, converting this radiation into electrical signals. A radiometer's primary components are an antenna and a Low Noise Amplifier (LNA), which is the core of the ``receiver'' chain. Instrumental effects introduced by the receiver are typically corrected or removed during calibration. However, impedance mismatches between the antenna and receiver can introduce unwanted signal reflections and distortions. Traditional calibration methods, such as Dicke switching, alternate the receiver input between the antenna and a well-characterised reference source to mitigate errors by comparison. Recent advances in Machine Learning (ML) offer promising alternatives. Neural networks, which are trained using known signal sources, provide a powerful means to model and calibrate complex systems where traditional analytical approaches struggle. These methods are especially relevant for detecting the faint sky-averaged 21-cm signal from atomic hydrogen at high redshifts. This is one of the main challenges in observational Cosmology today. Here, for the first time, we introduce and test a machine learning-based calibration framework capable of achieving the precision required for radiometric experiments aiming to detect the 21-cm line.

arXiv:2502.02647v2 Announce Type: replace Abstract: The field of high redshift galaxy formation has been revolutionised by JWST, which is yielding unprecedented insights on galaxy assembly at early times. Our key aim is to study the physical mechanisms that can explain the unexpected abundance of bright galaxies at $z \geq 11$, as well as their metal enrichment and spectral properties. We also use recent data to determine the key sources of reionisation. To do so, we implement cold gas fractions and star formation efficiencies derived from the SPHINX20 high-resolution radiation-hydrodynamics simulation into DELPHI, a semi-analytic model that tracks the assembly of dark matter halos and their baryonic components from $z \sim 4.5-40$. In addition, we explore two different methodologies to boost galaxy luminosities at $z \geq 11$: a stellar initial mass function (IMF) that becomes increasingly top-heavy with decreasing metallicity and increasing redshift (eIMF model), and star formation efficiencies that increase with increasing redshift (eSFE model). Our key findings are: (i) both the eIMF and eSFE models can explain the abundance of bright galaxies at $z \geq 11$; (ii) dust attenuation plays an important role for the bright-end of the UV LF at $z \leq 11$; (iii) the mass-metallicity relation is in place as early as $z \sim 17$ in all models although its slope is model-dependent; (iv) within the spread of both models and observations, all of our models are in good agreement with current estimates of $\beta$ slopes at $z \sim 5-17$ and Balmer break strengths at $z \sim 6-10$; (v) in the eIMF model, galaxies at $z\geq12$ or with $\rm{M_{UV}}\geq-18$ show values of $\xi_{\rm{ion}} \sim 10^{25.55}~{\rm [Hz~erg^{-1}]}$, twice larger than in other models; (vi) star formation in galaxies below $10^{9}\rm{M_{\odot}}$ is the key driver of reionisation, providing the bulk ($\sim 85\%$) of ionising photons down to its midpoint at $z \sim 7$.