We present three cosmological radiation-hydrodynamic zoom simulations of the progenitor of a Milky Way-mass galaxy from the MEGATRON suite. The simulations combine on-the-fly radiative transfer with a detailed non-equilibrium thermochemical network (81 ions and molecules), resolving the cold and warm gas in the circumgalactic medium (CGM) on spatial scales down to 20 pc and on average 200 pc at cosmic noon. Comparing our full non-equilibrium calculation with local radiation to traditional post-processed photoionization equilibrium (PIE) models assuming a uniform UV background (UVB), we find that non-equilibrium physics and local radiation fields fundamentally impact the thermochemistry of the CGM. Recombination lags and local radiation anisotropy shift ions away from their PIE+UVB values and modify covering fractions (for example, HI damped Ly$\alpha$ absorbers differ by up to 40%). In addition, a resolution study with cooling-length refinement allows us to double the resolution in the cold and warm CGM gas, reaching 120 pc on average. When refining on cooling length, the mass of the lightest cold clumps decreases tenfold to $\approx 10^4\,M_\odot$, their boundary layers develop sharper ion stratification, and the warm gas is better resolved, boosting the abundance of warm gas tracers such as CIV and OIII. Together, these results demonstrate that non-equilibrium thermochemistry coupled to radiative transfer, combined with physically motivated resolution criteria, is essential to predict circumgalactic absorption and emission signatures and to guide the design of targeted observations with existing and upcoming facilities.

Stellar surface inhomogeneities such as spots and faculae introduce Doppler variations that challenge exoplanet detection via the radial velocity method. While their impact on disc-integrated spectra is well established, detailed studies of the underlying local line profiles have so far been limited to the Sun. We present an observational campaign targeting the active star WASP-85 A during transits of its hot Jupiter companion. The transits span two stellar rotation periods, allowing us to probe the evolution of active regions. From ground-based photometry we identify seven active regions, six containing dark spots. Using simultaneous ESPRESSO transit spectroscopy, we spatially resolve these regions on the stellar surface by using the planet as a probe. We detect significant bisector shape changes, line broadening, and net redshifts during spot occultations, with velocity shifts of 108-333 m/s (mean uncertainty 50 m/s). The observed broadening is consistent with the Zeeman effect, implying magnetic field strengths (Stokes $I$) $B$ = 2.7-4.4 kG (mean uncertainty 0.6 kG), comparable to solar umbrae. Combined with our photometric spot model, this yields lower limits to the disc-integrated field $Bf = 16 \pm 3$ G and $61 \pm 9$ G for the two hemispheres probed -- at least three times higher than Sun-as-a-star values. We also measure centre-to-limb variations in FWHM, line depth, equivalent width, and convective blueshift, which broadly agree with solar observations and 3D MHD models. This work demonstrates a new way to characterise the surfaces of exoplanet host stars, paving the way for future analyses incorporating synthetic line profiles from 3D MHD simulations.

Dwarf galaxies provide powerful laboratories for studying galaxy formation physics. Their early assembly, shallow gravitational potentials, and bursty, clustered star formation histories make them especially sensitive to the processes that regulate baryons through multi-phase outflows. Using high-resolution, cosmological zoom-in simulations of a dwarf galaxy from \textit{the Pandora suite}, we explore the impact of stellar radiation, magnetic fields, and cosmic ray feedback on star formation, outflows, and metal retention. We find that our purely hydrodynamical model without non-thermal physics - in which supernova feedback is boosted to reproduce realistic stellar mass assembly - drives violent, overly enriched outflows that suppress the metal content of the host galaxy. Including radiation reduces the clustering of star formation and weakens feedback. However, the additional incorporation of cosmic rays produces fast, mass-loaded, multi-phase outflows consisting of both ionized and neutral gas components, in better agreement with observations. These outflows, which entrain a denser, more temperate ISM, exhibit broad metallicity distributions while preserving metals within the galaxy. Furthermore, the star formation history becomes more bursty, in agreement with recent JWST findings. These results highlight the essential role of non-thermal physics in galaxy evolution and the need to incorporate it in future galaxy formation models.

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.

Recent JWST observations have revealed a puzzling population of optically red and compact galaxies with peculiar "V"-shaped spectra at high redshift, known as "Little Red Dots" (LRDs). Until now, most spectroscopically confirmed LRDs are found at z > 4 and it has been speculated that LRDs are tracing the early stages of black hole evolution. We report an independent rediscovery of a broad-line active galactic nucleus (AGN), SDSS J102530.29+140207.3, at z = 0.1, which shows spectral features matching those of LRDs seen in the early Universe, including the V-shaped spectrum, broad Balmer lines (with widths of 1000-2000 km/s), and deep Balmer absorption. We present a new GTC observation of this LRD, which reveals an optical continuum similar to those of G-to-K giant stars including an unambiguous G-band absorption originating from the CH molecule. In addition, this local LRD shows a series of absorption lines potentially related to low-ionization ions or atoms but are deeper than what is observed in empirical stellar templates. We further identify a series of [FeII] emission lines indicative of low-ionization gas, which we find also present in an JWST-selected LRD at z = 2.26. We find small but statistically significant variability in H$\alpha$ consistent with previous findings. Finally, with the new X-ray observation from NuSTAR, we confirm the extreme X-ray weakness of this LRD, which might imply Compton-thick gas obscuration with $N_{\rm H}>10^{24}~{\rm cm^{-2}}$. All evidence suggests SDSS J102530.29+140207.3 has a complex gaseous environment and the strong ionic, atomic, and molecular absorptions are hard to explain with typical stellar and AGN models.

Resolved high-redshift galaxy gas kinematics is a rapidly evolving field driven by increasingly powerful instrumentation. However, the resolution and sensitivity still impose constraints on interpretation. We investigate the uncertainties inherent to high-$z$ galaxy kinematical analysis by modelling a suite of rotating disk galaxies, generating synthetic interferometric ALMA observations, and fitting them with the 3D-kinematical tools 3DBarolo, GalPaK3D, and Qubefit. We present the recovered 3D-fitted kinematical parameters to assess their reliability, quantify the range of values possible for individual source studies, and establish the systematic biases present for observed samples. The $V/\sigma_{\rm V}$ ratio, which indicates how dynamically cold a system is, is of particular importance and depends on the choice of 3D-fitting tool. On average, 3DBarolo and Qubefit slightly overestimates $V/\sigma_{\rm V}$ (<1\sigma) and GalPaK3D underestimates it (<2\sigma). Therefore, all three tools are reliable for kinematical studies of averages of high-redshift galaxy samples. The value range possible for individual sources is significant, however, even more so for samples of not purely rotation dominated sources. To determine whether an observed galaxy is rotation dominated enough to be fitted with a 3D-kinematical tool, $V/\sigma_{\rm V}$ can be extracted directly from the observed data cube, with some caveats. We recommend that the median offsets, value ranges, and tool-dependent biases presented in this paper are taken into account when interpreting 3D-fitted kinematics of observed high-redshift galaxies.

Wide-field slitless spectroscopy (WFSS) is a powerful tool for studying large samples of galaxies across cosmic times. With the arrival of JWST, and its NIRCAM grism mode, slitless spectroscopy can reach a medium spectral resolution of $(R\sim 1600)$, allowing it to spatially resolve the ionised-gas kinematics out to $z\sim 9$. However, the kinematic information is convolved with morphology along the dispersion axis, a degeneracy that must be modelled to recover intrinsic properties. We present the Grism Emission-line Kinematics tOol ($\texttt{geko}$), a Python package that forward-models NIRCam grism observations and infers emission-line morphologies and kinematics within a Bayesian framework. $\texttt{geko}$ combines Sérsic surface-brightness models with arctangent rotation curves, includes full point-spread function (PSF) and line-spread function (LSF) convolution, and leverages gradient-based sampling via $\texttt{jax}$/$\texttt{numpyro}$ for efficient inference. It recovers parameters such as effective radius, velocity dispersion, rotational velocity, rotational support, and dynamical mass, with typical run times of $\sim$20 minutes per galaxy on GPUs. We validate performance using extensive mock data spanning position angle, S/N, and morphology, quantifying where degeneracies limit recovery. Finally, we demonstrate applications to real FRESCO H$\alpha$ emitters at $z\approx 4-6$, recovering both rotation- and dispersion-dominated systems. $\texttt{geko}$ opens the way to statistical studies of galaxy dynamics in the early Universe and is publicly available at this https URL.

We report the discovery and characterization of three transiting giant planets in the TIC118798035 system. The three planets were identified as transiting candidates from data of the TESS mission, and confirmed with ground-based photometric transit observations along with radial velocity variations obtained with FEROS, HARPS and ESPRESSO. The three planets present transit timing variations (TTVs). We performed a N-body orbital fitting to the TTVs and radial velocities finding that TIC118798035 b is as warm low-density Neptune with a mass of 0.0250$\pm$0.0023 $M_J$, a radius of 0.655$\pm$0.018 $R_J$, and an orbital period of 11.507 d; TIC118798035 c is a warm Saturn with a mass of 0.403$\pm$0.024 $M_J$, a radius of 0.973$\pm$0.023 $R_J$, and an orbital period of 22.564 d; and TIC118798035 d is a warm Jupiter with a mass of 0.773$\pm$0.052 $M_J$, a radius of 0.923$\pm$0.044 $R_J$, and an orbital period of 48.925 d. The bulk metallicities of the three planets don't fully follow the mass-metallicity correlation found for the giant planets of the solar system, which hints at a somewhat different formation history for the planets of the TIC118798035 system. TIC118798035 is the only system having more than two transiting planets larger than 0.5 $R_J$ with a precise orbital and physical characterization, amenable for future atmospheric studies.

Recent JWST/NIRSpec observations have revealed high-$z$ star-forming galaxies depart from the Fundamental Metallicity Relation (FMR), yet the $z = 0$ FMR has not been well-characterized in the low-mass regime ($\rm log(M_{\star}/M_{\odot}) \lesssim 9$) for an appropriate comparison of low- and high-$z$ systems. We attempt to rectify this limitation through a meta-analysis, providing a local, observational comparison for future high-$z$ FMR studies. We analyzed common FMR fitting methods for $\sim 700$ [OIII]$\lambda 4363$ emitters with $\rm log(M_{\star}/M_{\odot}) \lesssim 9$ at $z \sim 0$. We find no evidence of the FMR below $\rm log(M_{\star}/M_{\odot}) \lesssim 9$ through any method, suggesting that slowly-evolving, quasi-steady state gas reservoirs are not yet established. We simultaneously find a weak positive correlation between metallicity and star formation, and that these systems are gas-rich with substantial diversity in effective yields ($y_{\rm eff}$) spanning $\rm \sim 3~dex$. We demonstrate increasing $y_{\rm eff}$ correlates with decreasing FMR offsets, which in the context of the analytical and non-equilibrium gas models of Dalcanton et al. (2007), indicates a scenario where star formation bursts rapidly return and eject metals from the ISM before subsequent gas-balancing. Pristine infall diluting the ISM metal-content cannot lead to the $y_{\rm eff}$ diversity we measure, and thus is not the primary process behind FMR deviations. Our results suggest low-$\rm M_{\star}$ systems, regardless of redshift, depart from a steady-state gas reservoir shaping the canonical FMR, in which metallicity variations are primarily driven by star formation and enriched outflows. With this characterization, we demonstrate $z \gtrsim 3$ [OIII]$\lambda 4363$ systems are indeed more metal-poor than $z \sim 0$ counterparts ($\rm \Delta 12+log(O/H) = 0.3~dex$) at fixed $\rm M_{\star}$.

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.

DES-5Y supernovae, combined with DESI BAO, appear to favour Chevallier-Polarski-Linder $(w_0, w_a)$ dynamical dark energy over $\Lambda$CDM. arXiv:2408.07175 suggested that this is driven by a systematic in the DES pipeline, which particularly affects the low-redshift supernovae brought in from legacy surveys. It is difficult to investigate these data in isolation, however, as the complicated supernovae pipelines must properly account for selection effects. In this work, we discover that the Bayesian evidence previously found for flexknot dark energy (arXiv:2503.17342) is beaten by a magnitude offset between the low- and high-redshift supernovae. In addition, we find that the possible tension between DES-5Y and DESI is significantly reduced by such an offset. We also take the opportunity to trial Nested Bridge Sampling with Sequential Monte Carlo as an alternative method for calculating Bayes factors.

The discovery of high-redshift galaxies exhibiting a steep spectral UV downturn potentially indicative of two-photon continuum emission marks a turning point in our search for signatures of top-heavy star formation in the early Universe. We develop a photometric search method for identifying further nebular-dominated galaxy candidates, whose nebular continuum dominates over the starlight, due to the high ionising photon production efficiencies $\xi_\mathrm{ion}$ associated with massive star formation. We utilise the extensive medium-band imaging from JADES, which enables the identification of Balmer jumps across a wide range of redshifts (1.5 < z < 8.5), through the deficit in rest-frame optical continuum level. As Balmer jumps are a general recombination feature of young starbursts ($\lesssim 3$~Myr), we further demand a high observed \log\, (\xi_\mathrm{ion, obs}/\mathrm{(Hz\ erg^{-1})}) > 25.60 to power the strong nebular continuum, together with a relatively non-blue UV slope indicating a lack of stellar continuum emission. Our nebular-dominated candidates, constituting ${\sim}$10% of galaxies at $z \sim 6$ (decreasing to ${\sim}$3% at $z \sim 2$, not completeness-corrected) are faint in the rest-frame optical (median $M_\mathrm{opt} = -17.95$) with extreme line emission (median $\mathrm{EW}_\mathrm{H\alpha,rest} = 1567$ Å, $\mathrm{EW}_\mathrm{[O\ III] + H\beta,rest} = 2244$ Å). However, hot H II region temperatures, collisionally-enhanced two-photon continuum emission, and strong UV lines are expected to accompany top-heavy star formation. Thus nebular-dominated galaxies do not necessarily exhibit the biggest Balmer jumps, nor the largest $\xi_\mathrm{ion, obs}$ or reddest UV slopes. Hence continuum spectroscopy is ultimately required to establish the presence of a two-photon downturn in our candidates, thus advancing our understanding of primordial star formation and AGN.

There is an ever-growing need in the gravitational wave community for fast and reliable inference methods, accompanied by an informative error bar. Nested sampling satisfies the last two requirements, but its computational cost can become prohibitive when using the most accurate waveform models. In this paper, we demonstrate the acceleration of nested sampling using a technique called posterior repartitioning. This method leverages nested sampling's unique ability to separate prior and likelihood contributions at the algorithmic level. Specifically, we define a `repartitioned prior' informed by the posterior from a low-resolution run. To construct this repartitioned prior, we use a $\beta$-flow, a novel type of conditional normalizing flow designed to better learn deep tail probabilities. $\beta$-flows are trained on the entire nested sampling run and conditioned on an inverse temperature $\beta$. Applying our methods to simulated and real binary black hole mergers, we demonstrate how they can reduce the number of likelihood evaluations required for convergence by up to an order of magnitude, enabling faster model comparison and parameter estimation. Furthermore, we highlight the robustness of using $\beta$-flows over standard normalizing flows to accelerate nested sampling. Notably, $\beta$-flows successfully recover the same posteriors and evidences as traditional nested sampling, even in cases where standard normalizing flows fail.

Galaxies exhibit a tight correlation between their star-formation rate and stellar mass over a wide redshift range known as the star-forming main sequence (SFMS). With JWST, we can now investigate the SFMS at high redshifts down to masses of $\sim10^6$ M$_{\odot}$, using sensitive star-formation rate tracers such as H$\alpha$ emission -- which allow us to probe the variability in star formation histories. We present inferences of the SFMS based on 316 H$\alpha$-selected galaxies at $z\sim4$-$5$ with $\log(\rm M_\star/M_\odot) = 6.4$ -$10.6$. These galaxies were identified behind the Abell 2744 lensing cluster with NIRCam grism spectroscopy from the ``All the Little Things'' (ALT) survey. At face value, our data suggest a shallow slope of the SFMS (SFR $\propto \mathrm{M}_\star^\alpha$, with $\alpha=0.45$). After correcting for the H$\alpha$-flux limited nature of our survey using a Bayesian framework, the slope steepens to $\alpha = 0.59^{+0.10}_{-0.09}$, whereas current data on their own are inconclusive on the mass dependence of the scatter. These slopes differ significantly from the slope of $\approx1$ expected from the observed evolution of the galaxy stellar mass function and from simulations. When fixing the slope to $\alpha=1$, we find evidence for a decreasing intrinsic scatter with stellar mass (from $\approx 0.5$ dex at M$_\star=10^8$ M$_\odot$ to $0.4$ dex at M$_\star=10^{10}$ M$_\odot$). This tension might be explained by a (combination of) luminosity-dependent SFR(H$\alpha$) calibration, a population of (mini)-quenched low-mass galaxies, or underestimated dust attenuation in high-mass galaxies. Future deep observations across facilities can quantify these processes, enabling better insights into the variability of star formation histories.

We present an overview of the James Webb Space Telescope (JWST) Advanced Deep Extragalactic Survey (JADES), an ambitious program of infrared imaging and spectroscopy in the GOODS-S and GOODS-N deep fields, designed to study galaxy evolution from high redshift to cosmic noon. JADES uses about 770 hours of Cycle 1 guaranteed time largely from the Near-Infrared Camera (NIRCam) and Near-Infrared Spectrograph (NIRSpec) instrument teams. In GOODS-S, in and around the Hubble Ultra Deep Field and Chandra Deep Field South, JADES produces a deep imaging region of ~45 arcmin$^2$ with an average of 130 hrs of exposure time spread over 9 NIRCam filters. This is extended at medium depth in GOODS-S and GOODS-N with NIRCam imaging of ~175 arcmin$^2$ with an average exposure time of 20 hrs spread over 8-10 filters. In both fields, we conduct extensive NIRSpec multi-object spectroscopy, including 2 deep pointings of 55 hrs exposure time, 14 medium pointings of ~12 hrs, and 15 shallower pointings of ~4 hrs, targeting over 5000 HST and JWST-detected faint sources with 5 low, medium, and high-resolution dispersers covering 0.6-5.3 microns. Finally, JADES extends redward via coordinated parallels with the JWST Mid-Infrared Instrument (MIRI), featuring ~9 arcmin$^2$ with 43 hours of exposure at 7.7 microns and twice that area with 2-6.5 hours of exposure at 12.8 microns For nearly 30 years, the GOODS-S and GOODS-N fields have been developed as the premier deep fields on the sky; JADES is now providing a compelling start on the JWST legacy in these fields.

A majority of JWST/NIRSpec/IFU studies at high redshifts to date have focused on UV-bright or massive objects, while our understanding of low-mass galaxies at early cosmic times remains limited. In this work, we present NIRSpec/IFS high-resolution observations of two low-mass (M_* < 10^9 \ M_\odot), low-metallicity ([12 + \log(\text{O/H})] < 8) galaxies at $z \sim 7.66$ with evidence of hosting AGN. Using spatially-resolved maps of rest-frame optical emission lines, we find flat metallicity profiles, indicative of ISM redistribution by outflows or past merging. We identify kinematical components decoupled from galactic rotation with velocities of $\sim 250 - 500 \ \text{km} \ \text{s}^{-1}$. We argue that these components are likely tracing outflows, possibly AGN-driven, for which we infer outflow rates of $\sim 21 - 40 \ M_\odot \ \text{yr}^{-1}$, suggesting they may suppress future star formation. We compare our observational results to those from the new large-volume AESOPICA simulations, which fully incorporate different models of black hole growth and AGN feedback. We find that our observational results of $v_\text{out}/v_\text{esc}$ and $\dot{M}_\text{out}$/SFR are consistent with the AGN population in these simulations, hinting that AGN-driven feedback may contribute to quenching both in our systems and in a wider population of low-mass galaxies in the early Universe. This novel study demonstrates the necessity of deep IFU observations to decompose the complex kinematics and morphology of high-$z$ galaxies, trace outflows, and constrain the effect of feedback in these low-mass systems.

When metals are magnetized, emulsions phase separate, or galaxies cluster, domain walls and patterns form and irremediably coarsen over time. Such coarsening is universally driven by diffusive relaxation toward equilibrium. Here, we discover an inertial counterpart - wave coarsening - in active elastic media, where vibrations emerge and spontaneously grow in wavelength, period, and amplitude, before a globally synchronized state called a time crystal forms. We observe wave coarsening in one- and two-dimensional solids and capture its dynamical scaling. We further arrest the process by breaking momentum conservation and reveal a far-from-equilibrium nonlinear analogue to chiral topological edge modes. Our work unveils the crucial role of symmetries in the formation of time crystals and opens avenues for the control of nonlinear vibrations in active materials.

Dust is a fundamental component of the interstellar medium within galaxies, as dust grains are highly efficient absorbers of ultraviolet (UV) and optical photons. Accurately quantifying this obscuration is crucial for interpreting galaxy spectral energy distributions (SEDs). The extinction curves in the Milky Way (MW) and Large Magellanic Cloud exhibit a strong feature known as the 2175 Å UV bump, most often attributed to small carbonaceous dust grains. This feature was recently detected in faint galaxies out to z=7.55, suggesting rapid formation channels. Here, we report the detection of a strong UV bump in a luminous Lyman-break galaxy at z_prism=7.11235, GNWY-7379420231, through observations taken as part of the NIRSpec Wide GTO survey. We fit a dust attenuation curve that is consistent with the MW extinction curve within 1{\sigma}, in a galaxy just ~700 Myr after the Big Bang. From the integrated spectrum, we infer a young mass-weighted age (t~22-59 Myr) for this galaxy, however spatially resolved SED fitting unveils the presence of an older stellar population (t~252 Myr). Furthermore, morphological analysis provides evidence for a potential merger. The underlying older stellar population suggests the merging system could be pre-enriched, with the dust illuminated by a merger-induced starburst. Moreover, turbulence driven by stellar feedback in this bursty region may be driving polycyclic aromatic hydrocarbon formation through top-down shattering. The presence of a UV bump in GNWY-7379420231 solidifies growing evidence for the rapid evolution of dust properties within the first billion years of cosmic time

Gas-phase metallicity is a fundamental parameter that helps constrain the star-forming history and chemical evolution of a galaxy. Measuring electron temperature through auroral-to-strong line ratios is a direct approach to deriving metallicity. However, there is a longstanding discrepancy between metallicity measured through the direct method and that based on the photoionization models. This paper aims to verify and understand the discrepancies. We bin ~ 1.5 million spaxels from SDSS-IV MaNGA according to metallicity and ionization parameters derived from theoretical strong-line calibrations. We stack the spectra of spaxels within each bin and measure the flux of strong lines and faint auroral lines. Auroral lines for [OII], [SII], [OIII], and [SIII] are detected in the stacked spectra of most bins, and the [NII] auroral line is detected in fewer bins. We apply an empirical method to correct dust attenuation, which makes more realistic corrections for low ionization lines. We derive electron temperatures for these five ionic species and measure the oxygen and sulfur abundances using the direct method. We present the resulting abundance measurements and compare them with those model-calibrated strong-line abundances. The chemical abundances measured with the direct method are lower than those derived from the photoionization model, with a median of 0.09 dex. This discrepancy is smaller compared to the results based on other metallicity calibrations previously reported. However, we notice that the direct method could not account for the variation in ionization parameters, indicating that the precise calibration of metallicity using the direct method has yet to be fully realized. We report significant discrepancies between data and the photoionization model, which illustrates that the one-dimensional photoionization model is incapable of representing the complexity of real situations.

Recent cosmological surveys have provided unprecedented datasets that can be used to reconstruct the history of the dark energy equation of state. In this work, a free-form "flexknot'' parameterisation is employed to represent $w(a)$ as a linear spline between free-moving nodes, the number of which may vary. By combining DESI Baryon Acoustic Oscillation measurements with Pantheon+ or DES5Y supernovae, the functional posteriors of $w(a)$ reveal an unexpected W-shaped structure. While the Bayesian evidence may still favour $\Lambda$CDM, the robustness of these results suggests the structure is indeed present in the data. The tension $R$-statistic and suspiciousness have been marginalised over models, and demonstrate that while the reconstructions from DESI and Pantheon+ agree, DESI and DES5Y do not. We conclude that, while there is no smoking gun for dynamical dark energy, the structure unearthed in this work is generally too complex to be captured by the restrictive $w$CDM or CPL parameterisations.