These lectures discuss multi-particle production in QCD and in gravity at ultrarelativistic energies, their double copy relations, and strong parallels in emergent shockwave dynamics. Dispersive techniques are applied to derive the BFKL equation for multi-gluon production in Regge asymptotics. Identical methods apply in gravity and are captured by a gravitational Lipatov equation. The building blocks in both cases are Lipatov vertices and reggeized propagators satisfying double copy relations; in gravity, Weinberg's soft theorem is recovered as a limit of the Lipatov framework. BFKL evolution in QCD generates wee parton states of maximal occupancy characterized by an emergent semi-hard saturation scale. A Color Glass Condensate EFT, and renormalization group equations in this framework, describe wee parton correlations, and their rapidity evolution. A shockwave picture of deeply inelastic scattering and hadron-hadron collisions follows, with multi-particle production described by Cutkosky's rules in strong time-dependent fields. Gluon radiation in this framework has a double copy in gravitational shockwave collisions, with a similar correspondence applicable to gluon and graviton shockwave propagators. Possible extensions of this double copy are outlined for computing multi-particle production in gravitational shockwave collisions, self-force and tidal contributions, and classical and quantum noise in the focusing of geodesics.

We discuss the cold dark matter plus massive neutrinos simulations of the MillenniumTNG project, which aim to improve understanding of how well ongoing and future large-scale galaxy surveys will measure neutrino masses. Our largest simulations, <inline-formula><tex-math>$3000\, {\rm Mpc}$</tex-math></inline-formula> on a side, use <inline-formula><tex-math>$10240^3$</tex-math></inline-formula> particles of mass <inline-formula><tex-math>$m_{\mathrm{ p}} = 6.66\times 10^{8}\, h^{-1}{\rm M}_\odot$</tex-math></inline-formula> to represent cold dark matter, and <inline-formula><tex-math>$2560^3$</tex-math></inline-formula> to represent a population of neutrinos with summed mass <inline-formula><tex-math>$M_\nu = 100$</tex-math></inline-formula> meV. Smaller volume runs with <inline-formula><tex-math>${\sim} 630\, {\rm Mpc}$</tex-math></inline-formula> also include cases with <inline-formula><tex-math>$M_\nu = 0\, \textrm {and}\, 300\, {\rm meV}$</tex-math></inline-formula>. All simulations are carried out twice using the paired-and-fixed technique for cosmic variance reduction. We evolve the neutrino component using the particle-based <inline-formula><tex-math>$\delta f$</tex-math></inline-formula> importance sampling method, which greatly reduces shot noise in the neutrino density field. In addition, we modify the GADGET-4 code to account both for the influence of relativistic and mildly relativistic components on the expansion rate and for non-Newtonian effects on the largest represented simulation scales. This allows us to quantify accurately the impact of neutrinos on basic statistical measures of non-linear structure formation, such as the matter power spectrum and the halo mass function. We use semi-analytic models of galaxy formation to predict the galaxy population and its clustering properties as a function of summed neutrino mass, finding significant (<inline-formula><tex-math>${\sim} 10~{{\ \rm per\ cent}}$</tex-math></inline-formula>) impacts on the cosmic star formation rate history, the galaxy mass function, and the clustering strength. This offers the prospect of identifying combinations of summary statistics that are optimally sensitive to the neutrino mass.

We introduce an extension of the evolution mapping framework to cosmological models that include massive neutrinos. The original evolution mapping framework exploits a degeneracy in the linear matter power spectrum when expressed in Mpc units, which compresses its dependence on cosmological parameters into those that affect its shape and a single extra parameter <inline-formula><mml:math><mml:msub><mml:mi>σ</mml:mi><mml:mn>12</mml:mn></mml:msub></mml:math></inline-formula>, defined as the rms linear variance in spheres of radius 12 Mpc. We show that by promoting the scalar amplitude of fluctuations, <inline-formula><mml:math><mml:msub><mml:mi>A</mml:mi><mml:mi>s</mml:mi></mml:msub></mml:math></inline-formula>, to a shape parameter, we can additionally describe the suppression due to massive neutrinos at any redshift to sub-0.01% accuracy across a wide range of masses and for different numbers of mass eigenstates. This methodology has been integrated into the public COMET package, enhancing its ability to emulate predictions of state-of-the-art perturbative models for galaxy clustering, such as the effective field theory (EFT) model. Additionally, the updated software now accommodates a broader cosmological parameter space for the emulator, enables the simultaneous generation of multiple predictions to reduce computation time, and incorporates analytic marginalization over nuisance parameters to expedite posterior estimation. Finally, we explore the impact of different infrared resummation techniques on galaxy power spectrum multipoles, demonstrating that any discrepancies can be mitigated by EFT counterterms without impacting the cosmological parameters.

We use the narrow [Ne v] $λ$3427 emission line detected in the recently published JWST spectra of two galaxies, at z = 6.9 and 5.6, to study the key properties of the active galactic nuclei (AGN) and the supermassive black holes in their centers. Using a new empirical scaling linking the [Ne v] line emission with AGN accretion-driven (continuum) emission, derived from a highly complete low-redshift AGN sample, we show that the [Ne v] emission in the two z > 5 galaxies implies total (bolometric) AGN luminosities of order L_bol~(4-8)x10^45 erg/s. Assuming that the radiation emitted from these systems is Eddington limited, the (minimal) black hole masses are of order M_BH>10^7 M_sun. Combined with the published stellar masses of the galaxies, estimated from dedicated fitting of their spectral energy distributions, the implied BH-to-stellar mass ratios are of order M_BH/M_host~0.1-1. This is considerably higher than what is found in the local Universe, but is consistent with the general trend seen in some other z > 5 AGN. Given the intrinsic weakness of the [Ne v] line and the nature of the [Ne v]-to-L_bol scaling, any (rare) detection of the [Ne v] $λ$3427 line at z > 5 would translate to similarly high AGN luminosities and SMBH masses, thus providing a unique observational path for studying luminous AGN well into the epoch of reionization, including obscured sources.

In the build-up of galactic discs gas infall is an important ingredient and it produces radial gas inflows as a physical consequence of angular momentum conservation, since the infalling gas on to the disc at a specific radius has lower angular momentum than the circular motions of the gas at the point of impact. NGC300 is a well studied isolated, bulge-less, and low-mass disc galaxy ideally suited for an investigation of galaxy evolution with radial gas inflows. To investigate the effects of radial gas inflows on the physical properties of NGC300, a chemical evolution model for NGC300 is constructed by assuming its disc builds up progressively by infalling of metall-free gas and outflowing of metal-enriched gas. Radial gas inflows are also considered in the model. Our model including the radial gas inflows and an inside-out disc formation scenario can simultaneously reproduce the present-day observed radial profiles of HI gas mass surface density, SFR surface density, sSFR, gas-phase and stellar metallicity. We find that, although the value of radial gas inflow velocity is as low as -0.1 km/s, the radial gas inflows steepen the present-day radial profiles of HI gas mass surface density, SFR surface density, and metallicity, but flatten the radial sSFR profile. Incorporating radial gas inflows significantly improves the agreement between our model predicted present-day sSFR profile and the observations of NGC300. It predicts a significant flattening of the metallicity gradient with cosmic time. We also find that the model predicted star formation has been more active recently, indicating that the radial gas inflows may be help to sustain star formation in local spirals, at least in NGC300.

The δN formalism has been the major computational tool to study the superhorizon evolution of the scalar type perturbation sourced by scalar fields. Recently, this formalism was generalized to compute an arbitrary scalar, vector, and tensor type perturbations, including the gravitational waves (GWs), sourced by an arbitrary bosonic fields. In this paper, we explain how to use the generalized δN formalism (the gδN formalism), considering a model with U(1) gauge fields as a concrete example. Several new findings on this model and prospects on future gravitational wave experiments are also discussed, including the condition for the two linear polarizations of GWs to have different amplitudes. This paper provides a detailed explanation of our previous paper published in Physical Review Letters. We also discuss the Weinberg's adiabatic mode for an anisotropic background, showing a qualitative difference from the one for the FLRW background.

A flagship application of quantum computers is the simulation of other quantum systems, including quantum field theories. In this article, we show how quantum computers can be employed to naturally calculate Feynman diagrams and their interferences in Quantum Chromodynamics (QCD). We simulate the colour parts of the interactions directly on the quantum computer, while the kinematic parts are for now pre-computed classically. For processes where some of the external particles are identical, we find the first hints of a potential quantum advantage. We validate our techniques using simulated quantum computers. Furthermore, for toy examples we also demonstrate our algorithms on a 56-qubit trapped-ion quantum computer. The work constitutes a further key step towards a full quantum simulation of generic perturbative QCD processes.

We investigate whether the effective theory for isolated, massive, and weakly interacting spin-$3/2$ particles is compatible with causality and unitarity-i.e., the positivity of scattering amplitudes. We find no solution to positivity constraints, except when gravitons are also present and couple in a (nearly) supersymmetric way. Gravity is thus bootstrapped from $S$-matrix consistency conditions for the longitudinal and transverse polarizations of massive spin-$3/2$ states. For two such particles forming a $U(1)$-charged state, a (gravi)photon gauging the symmetry is also required, with couplings characteristic of supergravity and consistent with both the no global symmetry and weak gravity conjectures. We further explore the EFT-hedron associated with the longitudinal polarizations, the Goldstinos, through novel $t$-$u$ symmetric dispersion relations. We identify the extremal UV models that lie at the corners of the allowed parameter space, recovering familiar models of supersymmetry breaking and uncovering new ones.

We present new observations that densely sample the microwave (4-360 GHz) continuum spectra from eight young systems in the Taurus region. Multi-component, empirical model prescriptions were used to disentangle the contributions from their dust disks and other emission mechanisms. We found partially optically thick, free-free emission in all these systems, with positive spectral indices (median $α_{\rm c} \approx 1$ at 10 GHz) and contributing 5-50% of the 43 GHz fluxes. There is no evidence for synchrotron or spinning dust grain emission contributions for these targets. The inferred dust disk spectra all show substantial curvature: their spectral indices decrease with frequency, from $α_{\rm d} \approx 2.8$-4.0 around 43 GHz to 1.7-2.1 around 340 GHz. This curvature suggests that a substantial fraction of the (sub)millimeter ($\gtrsim$ 200 GHz) dust emission may be optically thick, and therefore the traditional metrics for estimating dust masses are flawed. Assuming the emission at lower frequencies (43 GHz) is optically thin, the local spectral indices and fluxes were used to constrain the disk-averaged dust properties and estimate corresponding dust masses. These masses are roughly an order of magnitude higher ($\approx 1000 \, M_\oplus$) than those found from the traditional approach based on (sub)millimeter fluxes. These findings emphasize the value of broad spectral coverage - particularly extending to lower frequencies ($\sim$cm-band) - for accurately interpreting dust disk emission; such observations may help reshape our perspective on the available mass budgets for planet formation.

We present the ALMA-CRISTAL survey, an ALMA Cycle 8 Large Program designed to investigate the physical properties of star-forming galaxies at 4 ≲ z ≲ 6 through spatially resolved, multiwavelength observations. This survey targets 19 star-forming main-sequence galaxies selected from the ALPINE survey, using ALMA Band 7 observations to study [C II] 158 μm line emission and dust continuum, complemented by JWST/NIRCam and HST imaging to map stellar and UV emission. The CRISTAL sample expanded to 39 after including newly detected galaxies in the CRISTAL fields, archival data, and pilot study targets. The resulting dataset provides a detailed view of gas, dust, and stellar structures on kiloparsec scales at the end of the era of reionization. The survey reveals diverse morphologies and kinematics, including rotating disks, merging systems, [C II] emission tails from potential interactions, and clumpy star formation. Notably, the [C II] emission in many cases extends beyond the stellar light seen in HST and JWST imaging. Scientific highlights include CRISTAL-10, exhibiting an extreme [C II] deficit similar to Arp 220, and CRISTAL-13, where feedback from young star-forming clumps likely causes an offset between the stellar clumps and the peaks of [C II] emission. CRISTAL galaxies exhibit global [C II]/FIR ratios that decrease with increasing FIR luminosity, similar to trends seen in local galaxies but shifted to higher luminosities, likely due to their higher molecular gas content. CRISTAL galaxies also span a previously unexplored range of global FIR surface brightness at high-redshift, showing that high-redshift galaxies can have elevated [C II]/FIR ratios. These elevated ratios are likely influenced by factors such as lower-metallicity gas, the presence of significant extraplanar gas, and contributions from shock-excited gas.

The gas-phase metallicity affects heating and cooling processes in the star-forming galactic interstellar medium (ISM) as well as ionising luminosities, wind strengths, and lifetimes of massive stars. To investigate its impact, we conduct magnetohydrodynamic simulations of the ISM using the FLASH code as part of the SILCC project. The simulations assume a gas surface density of 10 M$_\odot$ pc$^{-2}$ and span metallicities from 1/50 Z$_\odot$ to 1 Z$_\odot$. We include non-equilibrium thermo-chemistry, a space- and time-variable far-UV background and cosmic ray ionisation rate, metal-dependent stellar tracks, the formation of HII regions, stellar winds, type II supernovae, and cosmic ray injection and transport. With the metallicity decreasing over the investigated range, the star formation rate decreases by more than a factor of ten, the mass fraction of cold gas decreases from 60% to 2.3%, while the volume filling fraction of the warm gas increases from 20% to 80%. Furthermore, the fraction of H$_\mathrm{2}$ in the densest regions drops by a factor of four, and the dense ISM fragments into approximately five times fewer structures at the lowest metallicity. Outflow mass loading factors remain largely unchanged, with values close to unity, except for a significant decline at the lowest metallicity. Including the major processes that regulate ISM properties, this study highlights the strong impact of gas phase metallicity on the star-forming ISM.

The M-theoretic emergence proposal claims that in an isotropic decompactification limit to M-theory the full effective action is generated via quantum effects by integrating out only the light towers of states of the theory. In the BPS particle sector, these include transversally wrapped M2- and M5-branes possibly carrying Kaluza-Klein momentum. This implies that a longitudinally wrapped M5-brane, i.e. a wrapped D4-brane, is not to be included in emergence computations. In this work we collect explicit evidence supporting this point by examining an F4 gauge coupling in six dimensions, making use of the duality between heterotic string theory on T4 and strongly coupled type IIA on K3. In this instance, the M-theoretic emergence proposal can be viewed as a tool for making predictions for the microscopic behavior of string theoretic amplitudes.

We investigate the high-ionization, narrow [Ne v] $λ$3427 emission line in a sample of over 340 ultrahard X-ray (14-195 keV) selected Active Galactic Nuclei (AGN) drawn from the BASS project. The analysis includes measurements in individual and stacked spectra, and considers several key AGN properties such as X-ray luminosity, supermassive black hole (SMBH) mass, Eddington ratios, and line-of-sight column density. The [Ne v] $λ$3427 line is robustly detected in ~43% (146/341) of the AGN in our sample, with no significant trends between the detection rate and key AGN/SMBH properties. In particular, the detection rate remains high even at the highest levels of obscuration (>70% for log[N_H/cm^-2] > 23). On the other hand, even some of our highest signal-to-noise spectra (S/N > 50) lack a robust [Ne v] detection. The typical (median) scaling ratios between [Ne v] line emission and (ultra-)hard X-ray emission in our sample are log L[Ne v]/L(14-150 keV) = -3.75 and log L[Ne v]/L(2-10 keV) = -3.36. The scatter on these scaling ratios, of ~0.5 dex, is comparable to, and indeed smaller than, what is found for other commonly used tracers of AGN radiative outputs (e.g., [O III] $λ$5007). Otherwise, we find no significant relations between the (relative) strength of [Ne v] and the basic AGN/SMBH properties under study, in contrast with simple expectations from models of SMBH accretion flows. Our results reaffirm the usability of [Ne v] as an AGN tracer even in highly obscured systems, including dual AGN and high redshift sources.

The obliquity between the stellar spin axis and the planetary orbit, detected via the Rossiter-McLaughlin (RM) effect, is a tracer of the formation history of planetary systems. While obliquity measurements have been extensively applied to hot Jupiters and short-period planets, they remain rare for cold and long-period planets due to observational challenges, particularly their long transit durations. We report the detection of the RM effect for the 19-hour-long transit of HIP 41378 f, a temperate giant planet on a 542-day orbit, observed through a worldwide spectroscopic campaign. We measure a slight projected obliquity of 21 $\pm$ 8 degrees and a significant 3D spin-orbit angle of 52 $\pm$ 6 degrees, based on the measurement of the stellar rotation period. HIP 41378 f is part of a 5-transiting planetary system with planets close to mean motion resonances. The observed misalignment likely reflects a primordial tilt of the stellar spin axis relative to the protoplanetary disk, rather than dynamical interactions. HIP 41378 f is the first non-eccentric long-period (P>100 days) planet observed with the RM effect, opening new constraints on planetary formation theories. This observation should motivate the exploration of planetary obliquities across a longer range of orbital distances through international collaboration.

Cosmic microwave background (CMB) photons are deflected by large-scale structure through gravitational lensing. This secondary effect introduces higher-order correlations in CMB anisotropies, which are used to reconstruct lensing deflections. This allows mapping of the integrated matter distribution along the line of sight, probing the growth of structure, and recovering an undistorted view of the last-scattering surface. Gravitational lensing has been measured by previous CMB experiments, with $\textit{Planck}$'s $42\,σ$ detection being the current best full-sky lensing map. We present an enhanced $\textit{LiteBIRD}$ lensing map by extending the CMB multipole range and including the minimum-variance estimation, leading to a $49$ to $58\,σ$ detection over $80\,\%$ of the sky, depending on the final complexity of polarized Galactic emission. The combination of $\textit{Planck}$ and $\textit{LiteBIRD}$ will be the best full-sky lensing map in the 2030s, providing a $72$ to $78\,σ$ detection over $80\,\%$ of the sky, almost doubling $\textit{Planck}$'s sensitivity. Finally, we explore different applications of the lensing map, including cosmological parameter estimation using a lensing-only likelihood and internal delensing, showing that the combination of both experiments leads to improved constraints. The combination of $\textit{Planck}$ + $\textit{LiteBIRD}$ will improve the $S_8$ constraint by a factor of 2 compared to $\textit{Planck}$, and $\textit{Planck}$ + $\textit{LiteBIRD}$ internal delensing will improve $\textit{LiteBIRD}$'s tensor-to-scalar ratio constraint by $6\,\%$. We have tested the robustness of our results against foreground models of different complexity, showing that a significant improvement remains even for the most complex foregrounds.

We derive constraints on the neutrino mass using a variety of recent cosmological datasets, including DESI BAO, the full-shape analysis of the DESI matter power spectrum and the one-dimensional power spectrum of the Lyman-$α$ forest (P1D) from eBOSS quasars as well as the cosmic microwave background (CMB). The constraints are obtained in the frequentist formalism by constructing profile likelihoods and applying the Feldman-Cousins prescription to compute confidence intervals. This method avoids potential prior and volume effects that may arise in a comparable Bayesian analysis. Parabolic fits to the profiles allow one to distinguish changes in the upper limits from variations in the constraining power $σ$ of the different data combinations. We find that all profiles in the $Λ$CDM model are cut off by the $\sum m_ν\geq 0$ bound, meaning that the corresponding parabolas reach their minimum in the unphysical sector. The most stringent 95% C.L. upper limit is obtained by the combination of DESI DR2 BAO, Planck PR4 and CMB lensing at 53 meV, below the minimum of 59 meV set by the normal ordering. Extending $Λ$CDM to non-zero curvature and $w_0w_\mathrm{a}$CDM relaxes the constraints past 59 meV again, but only $w_0w_\mathrm{a}$CDM exhibits profiles with a minimum at a positive value. Using a combination of DESI DR1 full-shape, BBN and eBOSS Lyman-$α$ P1D, we successfully constrain the neutrino mass independently of the CMB. This combination yields $\sum m_ν\leq 285$ meV (95% C.L.). The addition of DESI full-shape or Lyman-$α$ P1D to CMB and DESI BAO results in small but noticeable improvement of the constraining power of the data. Lyman-$α$ free-streaming measurements especially improve the constraint. Since they are based on eBOSS data, this sets a promising precedent for upcoming DESI data.

The LiteBIRD satellite mission aims at detecting Cosmic Microwave Background $B$ modes with unprecedented precision, targeting a total error on the tensor-to-scalar ratio $r$ of $δr \sim 0.001$. Operating from the L2 Lagrangian point of the Sun-Earth system, LiteBIRD will survey the full sky across 15 frequency bands (34 to 448 GHz) for 3 years.The current LiteBIRD baseline configuration employs 4508 detectors sampling at 19.1 Hz to achieve an effective polarization sensitivity of $ 2 μ\mathrm{K-arcmin}$ and an angular resolution of 31 arcmin (at 140 GHz).We describe the first release of the official LiteBIRD simulations, realized with a new simulation pipeline developed using the LiteBIRD Simulation Framework, see github.com/litebird/litebird_sim . This pipeline generates 500 full-sky simulated maps at a Healpix resolution of nside=512. The simulations include also one year of Time Ordered Data for approximately one-third of LiteBIRD's total detectors.

LiteBIRD, the Lite (Light) satellite for the study of $B$-mode polarization and Inflation from cosmic background Radiation Detection, is a space mission focused on primordial cosmology and fundamental physics. In this paper, we present the LiteBIRD Simulation Framework (LBS), a Python package designed for the implementation of pipelines that model the outputs of the data acquisition process from the three instruments on the LiteBIRD spacecraft: LFT (Low-Frequency Telescope), MFT (Mid-Frequency Telescope), and HFT (High-Frequency Telescope). LBS provides several modules to simulate the scanning strategy of the telescopes, the measurement of realistic polarized radiation coming from the sky (including the Cosmic Microwave Background itself, the Solar and Kinematic dipole, and the diffuse foregrounds emitted by the Galaxy), the generation of instrumental noise and the effect of systematic errors, like pointing wobbling, non-idealities in the Half-Wave Plate, et cetera. Additionally, we present the implementation of a simple but complete pipeline that showcases the main features of LBS. We also discuss how we ensured that LBS lets people develop pipelines whose results are accurate and reproducible. A full end-to-end pipeline has been developed using LBS to characterize the scientific performance of the LiteBIRD experiment. This pipeline and the results of the first simulation run are presented in Puglisi et al. (2025).

We assess the computational feasibility of end-to-end Bayesian analysis of the JAXA-led LiteBIRD experiment by analysing simulated time ordered data (TOD) for a subset of detectors through the Cosmoglobe and Commander3 framework. The data volume for the simulated TOD is 1.55 TB, or 470 GB after Huffman compression. From this we estimate a total data volume of 238 TB for the full three year mission, or 70 TB after Huffman compression. We further estimate the running time for one Gibbs sample, from TOD to cosmological parameters, to be approximately 3000 CPU hours. The current simulations are based on an ideal instrument model, only including correlated 1/f noise. Future work will consider realistic systematics with full end-to-end error propagation. We conclude that these requirements are well within capabilities of future high-performance computing systems.

Precision spectroscopy of the electron spectrum of the tritium <inline-formula><mml:math><mml:mi>β</mml:mi></mml:math></inline-formula>-decay near the kinematic endpoint is a direct method to determine the effective electron antineutrino mass. The KArlsruhe TRItium Neutrino (KATRIN) experiment aims to determine this quantity with a sensitivity of better than <inline-formula><mml:math><mml:mrow><mml:mrow><mml:mn>0.3</mml:mn></mml:mrow><mml:mrow><mml:mspace></mml:mspace><mml:mtext>eV</mml:mtext></mml:mrow></mml:mrow></mml:math></inline-formula> (<inline-formula><mml:math><mml:mrow><mml:mn>90</mml:mn><mml:mo>%</mml:mo></mml:mrow></mml:math></inline-formula> C.L.). An inhomogeneous electric potential in the tritium source of KATRIN can lead to distortions of the <inline-formula><mml:math><mml:mi>β</mml:mi></mml:math></inline-formula>-spectrum, which directly impact the neutrino-mass observable. This effect can be quantified through precision spectroscopy of the conversion-electrons of co-circulated metastable <inline-formula><mml:math><mml:mrow><mml:mmultiscripts><mml:mrow></mml:mrow><mml:mrow></mml:mrow><mml:mrow><mml:mn>83</mml:mn><mml:mtext>m</mml:mtext></mml:mrow></mml:mmultiscripts><mml:mtext>Kr</mml:mtext></mml:mrow></mml:math></inline-formula>. Therefore, dedicated, several-weeks long measurement campaigns have been performed within the KATRIN data taking schedule. In this work, we infer the tritium source potential observables from these measurements, and present their implications for the neutrino-mass determination.

Galaxy evolution is driven by spatially distributed processes with varying time-scales. Integral field spectroscopy provides spatially resolved information about these processes. Nevertheless, disentangling these processes, which are related to both the underlying stellar populations and the interstellar medium can be challenging. We present a case study on NGC 613, observed with MUSE (Multi-Unit Spectroscopic Explorer) for the TIMER (Time Inference with MUSE in Extragalactic Rings) project, a local barred galaxy, which shows several gas ionization mechanisms and is rich in both large and inner-scale stellar structures. We develop a set of steps to overcome fundamental problems in the modelling of emission lines with multiple components, together with the characterization of the stellar populations. That results in the disentanglement of the gas ionization mechanisms and kinematics, along with an optimal parametrization for star formation history recovery. Our analysis reveals evidence of gas inflows, which are associated with the bar dust lanes traced with Hubble Space Telescope. In addition, we show the gas kinematics in a central biconical outflow, which is aligned with a radio jet observed with Very Large Array. The emission line provides estimates of electron density, gas-phase metallicity, and the mass outflow rate, allowing us to distinguish intertwined ionization mechanisms and to identify a part of the multiphase gas cycle in NGC 613. It traces the gas kinematics from the bar lanes to inner scale gas reservoirs, where it can eventually trigger star formation or AGN activity, as observed in the outflow.

We present primordial non-Gaussianity predictions from a new high-precision code for simulating axion-U(1) inflation on a discrete lattice. We measure the primordial scalar curvature power spectrum and bispectrum from our simulations, determining their dependence on both scale and axion-gauge coupling strength. Both the gauge-sourced power spectrum and the bispectrum exhibit a strong blue tilt due to our choice of an $α$-attractor inflaton potential. We provide fitting functions for the power spectrum and bispectrum that accurately reproduce these statistics across a wide range of scales and coupling strengths. While our fitting function for the bispectrum has a separable form, results from high-resolution simulations demonstrate that the full shape is not separable. Thus, our simulations generate realizations of primordial curvature perturbations with nontrivial correlators that cannot be generated using standard techniques for primordial non-Gaussianity. We derive bounds on the axion-gauge coupling strength based on the bispectrum constraints from the cosmic microwave background, demonstrating a new method for constraining inflationary primordial non-Gaussianity by simulating the nonlinear dynamics.

Cosmological correlators are important observables in cosmology. They are often approximated by de Sitter space correlators. In this paper, we give a first precise diagrammatical computation of higher loop diagrams to all orders for a conformally coupled scalar in four dimensions. We show that, in contrast to flat space, diagrams of necklace topology do not resum using the natural application of de Sitter-invariant regularization and are thus hard to evaluate. We propose a modification to the UV-regularization of the loops, compatible with de Sitter invariance, but much easier to work with. The modified diagrams can be resummed to give a glimpse at non-perturbative effects for de Sitter correlators. Furthermore, they fit nicely with recently proposed cosmological dressing rules.

Context. Planetary nebulae (PNe) and their luminosity function (PNLF) in galaxies have been used as a cosmic distance indicator for decades, yet a fundamental understanding is still lacking to explain the universality of the PNLF among different galaxies. So far, models for the PNLF have generally assumed near-solar metallicities and employed simplified stellar populations. Aims. In this work, we investigate how metallicity and helium abundances affect the resulting PNe and PNLF as well as the importance of the initial-to-final mass relation (IFMR) and circumnebular extinction in order to resolve the tension between PNLF observations and previous models. Methods. We introduce PICS (PNe In Cosmological Simulations), a PN model framework that takes into account the stellar metal-licity and is applicable to realistic stellar populations obtained from both cosmological simulations and observations. The framework combines current stellar evolution models with post-AGB tracks and PN models to obtain the PNe from the parent stellar population. Results. We find that metallicity plays an important role for the resulting PNe, as old metal-rich populations can harbor much brighter PNe than old metal-poor populations. In addition, we show that the helium abundance is a vital ingredient at high metallicities, and we explored the impact on the PNLF of a possible saturation of the helium content at higher metallicities. We present PNLF grids for different stellar ages and metallicities, where the observed PNLF bright end can be reached even for old stellar populations of 10 Gyr at high metallicities. Finally, we find that the PNLFs of old stellar populations are extremely sensitive to the IFMR, potentially allowing for the production of bright PNe. Conclusions. With PICS, we have laid the groundwork for studying how models and assumptions relevant to PNe affect the PNe and PNLF. Two of the central ingredients for the PNe and PNLF are the metallicity and helium abundance. Future applications of PICS include self-consistent modeling of PNe in a cosmological framework to explain the origin of the universality of the PNLF bright-end cutoff and using it as a diagnostic tool for galaxy formation.

We perform a thorough investigation of the universality of the long-distance matrix elements (LDMEs) of nonrelativistic QCD factorization based on a next-to-leading order (NLO) fit of <inline-formula><mml:math><mml:mi>J</mml:mi><mml:mo>/</mml:mo><mml:mi>ψ</mml:mi></mml:math></inline-formula> color octet LDMEs to high transverse momentum <inline-formula><mml:math><mml:msub><mml:mi>p</mml:mi><mml:mi>T</mml:mi></mml:msub></mml:math></inline-formula> <inline-formula><mml:math><mml:mi>J</mml:mi><mml:mo>/</mml:mo><mml:mi>ψ</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math><mml:msub><mml:mi>η</mml:mi><mml:mi>c</mml:mi></mml:msub></mml:math></inline-formula> production data at the LHC. We thereby apply a novel fit-and-predict procedure to systematically take into account scale variations, and predict various observables never studied in this context before. In particular, the LDMEs can well describe <inline-formula><mml:math><mml:mi>J</mml:mi><mml:mo>/</mml:mo><mml:mi>ψ</mml:mi></mml:math></inline-formula> hadroproduction up to the highest measured values of <inline-formula><mml:math><mml:msub><mml:mi>p</mml:mi><mml:mi>T</mml:mi></mml:msub></mml:math></inline-formula>, as well as <inline-formula><mml:math><mml:mrow><mml:mi>ϒ</mml:mi><mml:mo>(</mml:mo><mml:mi>n</mml:mi><mml:mi>S</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> production via potential nonrelativistic QCD based relations. Furthermore, <inline-formula><mml:math><mml:mi>J</mml:mi><mml:mo>/</mml:mo><mml:mi>ψ</mml:mi></mml:math></inline-formula> production in <inline-formula><mml:math><mml:mi>γ</mml:mi><mml:mi>γ</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math><mml:mi>γ</mml:mi><mml:mi>p</mml:mi></mml:math></inline-formula> collisions is surprisingly reproduced down to <inline-formula><mml:math><mml:msub><mml:mi>p</mml:mi><mml:mi>T</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mtext> </mml:mtext><mml:mtext> </mml:mtext><mml:mi>GeV</mml:mi></mml:math></inline-formula>, as long as the region of large inelasticity <inline-formula><mml:math><mml:mi>z</mml:mi></mml:math></inline-formula> is excluded, which may be of significance in future quarkonium studies, in particular at the EIC and the high-luminosity LHC. In addition, our summary reveals an interesting pattern as to which observables still evade a consistent description.

Progress in technology has always been a driver of progress in science. Answering challenging questions requires bespoke technology and strategies in its use. The expanding toolbox of bottom-up synthetic biology allows us to study the organization principles underlying life. This include artificial membrane systems, phase-separating protein droplets, nucleic acids programmed to encode and use information, and many others. To study the details of the underlying physical chemistry, novel biophysical techniques are often needed. Fluorescence techniques allow studying the dynamics of biomolecules with high sensitivity and at high spatial and temporal resolution. [...]

Since the dawn of physics, the task of the experimental physicist has been to design a setup and perform measurements. A new development has become clear: After the data has been measured, modern computation allows us to further extract information. This work presents a collection of new techniques based on numerical methods in statistics and machine learning for the IceCube Neutrino Observatory. It aims to extend the analysis possibilities in particle physics and neutrino astronomy.

Neutron beta decay offers many tests of the Standard Model of particle physics, including the determination of the beta asymmetry and the Fierz interference term. I present systematic studies regarding the recent Perkeo III measurement campaign with the aim to extract the Fierz term and the electron backscatter behavior in both Perkeo III and the new spectrometer PERC. I also present a design for PERC’s secondary detector and a measurement of the neutron depolarization induced by supermirrors.

Supernovae (SNe) are among the most energetic events in the Universe, injecting vast amounts of energy and enriched material into the interstellar medium (ISM). Their ashes, usually referred to as supernova remnants (SNRs) or superbubbles (SBs) in the case SNRs powered by the clustered explosions of many stars, play a central role in shaping the multiphase structure of the ISM, driving turbulence, regulating star formation, and contributing to galactic outflows. Despite their importance, many aspects of their long-term evolution and interaction with their galactic environment remain poorly understood, particularly in the context of a realistic, structured ISM.

In this thesis, I explore the dynamical evolution of SNRs across a range of scales and physical conditions. I begin by developing a novel one-zone blastwave model, based on the thin-shell and sector approximations, that is designed to capture the expansion of SNRs in non-uniform, time-dependent environments, including the effects of gravity, shear induced by differential rotation, and cooling. This model reproduces known analytic limits and serves as a foundation for exploring more complex configurations.

We present a readout electronics system designed for Compton scattering detection using silicon drift detectors (SDDs). This system represents the first component of a Compton telescope to be integrated into a CubeSat module as part of the ComPol project. The instrument is designed to observe Cygnus X-1, a binary system consisting of a black hole and a blue supergiant star, with the aim of measuring its X-ray spectrum and polarization. Through the detection of Compton events, the system will determine the angle and degree of polarization of the source. The payload consists of a Si-based detector to scatter the X-rays and a scintillating cerium bromide crystal (CeBr3) with silicon photon multiplier (SiPM) matrix for the detection of the scattered photons. The analog readout chain of the first detection system based on SDD matrices consists of charge preamplifiers, an analog pulse processor (ASIC), while a field programmable gate array (FPGA) module handles the digital data flow to the PC. Experimental characterisation of the prototype shows an energy resolution of 151 eV full width at half maximum (FWHM) with a shaping time of 3 μs. This paper presents the prototype boards designed for the final flight module, together with the results of tests performed with the prototype system.

In the past, researchers have mostly relied on single-resolution images from individual telescopes to detect gravitational lenses. We present a search for galaxy-scale lenses that, for the first time, combines high-resolution single-band images (in our case from the Hubble Space Telescope, HST) with lower-resolution multiband images (in our case from the Legacy survey, LS) using machine learning. This methodology simulates the operational strategies employed by future missions, such as combining the images of Euclid and the Rubin Observatory's Legacy Survey of Space and Time (LSST). To compensate for the scarcity of lensed galaxy images for network training, we generated mock lenses by superimposing arc features onto HST images, saved the lens parameters, and replicated the lens system in the LS images. We tested four architectures based on ResNet-18: (1) using single-band HST images, (2) using three bands of LS images, (3) stacking these images after interpolating the LS images to HST pixel scale for simultaneous processing, and (4) merging a ResNet branch of HST with a ResNet branch of LS before the fully connected layer. We compared these architecture performances by creating receiver operating characteristic (ROC) curves for each model and comparing their output scores. At a false-positive rate of 10^‑4, the true-positive rate is ∼0.41, ∼0.45, ∼0.51 and ∼0.55, for HST, LS, stacked images and merged branches, respectively. Our results demonstrate that models integrating images from both the HST and LS significantly enhance the detection of galaxy-scale lenses compared to models relying on data from a single instrument. These results show the potential benefits of using both Euclid and LSST images, as wide-field imaging surveys are expected to discover approximately 100 000 lenses.

We used strong gravitational lensing to study the mass distribution of the galaxy cluster MACS J0035.4‑2015, by modeling its total mass distribution. The combination of high-resolution imaging from the Hubble Space Telescope with ground-based spectroscopy from the Multi Unit Spectroscopic Explorer mounted at the Very Large Telescope allowed us to model the observed multiple image positions with ≈0.″3 precision. We find that MACS J0035.4‑2015 can be best described by a combination of an elliptical dark matter halo modeled as an isothermal mass profile, with the brightest cluster galaxy and cluster members each modeled with a spherical truncated isothermal parameterization. With these assumptions, the total mass is estimated to be ≈ 6 × 10¹³ M_⊙ within 100 kpc. The data and mass model presented here form the basis for future cosmological and astrophysical studies of this cluster.

In this work, we significantly enhance masked particle modeling (MPM), a self-supervised learning scheme for constructing highly expressive representations of unordered sets relevant to developing foundation models for high-energy physics. In MPM, a model is trained to recover the missing elements of a set, a learning objective that requires no labels and can be applied directly to experimental data. We achieve significant performance improvements over previous work on MPM by addressing inefficiencies in the implementation and incorporating a more powerful decoder. We compare several pre-training tasks and introduce new reconstruction methods that utilize conditional generative models without data tokenization or discretization. We show that these new methods outperform the tokenized learning objective from the original MPM on a new test bed for foundation models for jets, which includes using a wide variety of downstream tasks relevant to jet physics, such as classification, secondary vertex finding, and track identification.

RNA and proteins are the foundation of life and a natural starting point to explore its origins. However, the prebiotic relationship between the two is asymmetric. While RNA evolved to assemble proteins from amino acids, a significant mirror-symmetric effect of amino acids to trigger the synthesis of RNA was missing. We describe ambient alkaline conditions where amino acids, without additional chemical activators, promote RNA copolymerisation more than 100-fold, starting from prebiotically plausible ribonucleoside-2',3'-cyclic phosphates. The observed effect is explained by acid-base catalysis, with optimal efficiency at pH values near the amine pK_aH. The fold-change in oligomerisation yield is nucleobase-selective, resulting in increased compositional diversity necessary for subsequent molecular evolution and favouring the formation of natural 3'‑5' linkages. The elevated pH offers recycling of oligonucleotide sequences back to 2',3'-cyclic phosphates, providing conditions for high-fidelity replication by templated ligation. The findings reveal a clear functional role of amino acids in the evolution of RNA earlier than previously assumed.

Our understanding of the γ-ray sky has improved dramatically in the past decade, however, the unresolved γ-ray background (UGRB) still has a potential wealth of information about the faintest γ-ray sources pervading the Universe. Statistical cross-correlations with tracers of cosmic structure can indirectly identify the populations that most characterize the γ-ray background. In this study, we analyze the angular correlation between the γ-ray background and the matter distribution in the Universe as traced by gravitational lensing, leveraging more than a decade of observations from the Fermi-Large Area Telescope (LAT) and 3 years of data from the Dark Energy Survey (DES). We detect a correlation at signal-to-noise ratio of 8.9. Most of the statistical significance comes from large scales, demonstrating, for the first time, that a substantial portion of the UGRB aligns with the mass clustering of the Universe as traced by weak lensing. Blazars provide a plausible explanation for this signal, especially if those contributing to the correlation reside in halos of large mass (∼ 10¹⁴ M _⊙) and account for approximately 30–40% of the UGRB above 10 GeV. Additionally, we observe a preference for a curved γ-ray energy spectrum, with a log-parabolic shape being favored over a power-law. We also discuss the possibility of modifications to the blazar model and the inclusion of additional γ-ray sources, such as star-forming galaxies, misalinged active galactic nuclei, or particle dark matter.

Quasi-periodic eruptions (QPEs) are rapid, recurring X-ray bursts from supermassive black holes, believed to result from interactions between accretion disks and surrounding matter. The galaxy SDSS1335+0728, previously stable for two decades, exhibited an increase in optical brightness in December 2019, followed by persistent active galactic nucleus (AGN)-like variability for 5 yr, suggesting the activation of a ~10⁶-M_⊙ black hole. Since February 2024, X-ray emission has been detected, revealing extreme ~4.5-d QPEs with high fluxes and amplitudes, long timescales, large integrated energies and a ~25-d superperiod. Low-significance UV variations are reported, probably related to the long timescales and large radii from which the emission originates. This discovery broadens the possible formation channels for QPEs, suggesting that they are linked not solely to tidal disruption events but more generally to newly formed accretion flows, which we are witnessing in real time in a turn-on AGN candidate.

The R3B experiment at FAIR studies nuclear reactions using high-energy radioactive beams. One key detector in R3B is the CALIFA calorimeter consisting of 2544 CsI(Tl) scintillator crystals designed to detect light charged particles and gamma rays with an energy resolution in the per cent range after Doppler correction. Precise cluster reconstruction from sparse hit patterns is a crucial requirement. Standard algorithms typically use fixed cluster sizes or geometric thresholds. To enhance performance, advanced machine learning techniques such as agglomerative clustering were implemented to use the full multi-dimensional parameter space including geometry, energy and time of individual interactions. An Edge Detection Neural Network exhibited significant differences. This study, based on Geant4 simulations, demonstrates improvements in cluster reconstruction efficiency of more than 30%, showcasing the potential of machine learning in nuclear physics experiments.

We show how a method to construct canonical differential equations for multi-loop Feynman integrals recently introduced by some of the authors can be extended to cases where the associated geometry is of Calabi-Yau type and even beyond. This can be achieved by supplementing the method with information from the mixed Hodge structure of the underlying geometry. We apply these ideas to specific classes of integrals whose associated geometry is a one-parameter family of Calabi-Yau varieties, and we argue that the method can always be successfully applied to those cases. Moreover, we perform an in-depth study of the properties of the resulting canonical differential equations. In particular, we show that the resulting canonical basis is equivalent to the one obtained by an alternative method recently introduced in the literature. We apply our method to non-trivial and cutting-edge examples of Feynman integrals necessary for gravitational wave scattering, further showcasing its power and flexibility.

High-energy physics requires the generation of large numbers of simulated data samples from complex but analytically tractable distributions called matrix elements. Surrogate models, such as normalizing flows, are gaining popularity for this task due to their computational efficiency. We adopt an approach based on flow annealed importance sampling bootstrap (FAB) that evaluates the differentiable target density during training and helps avoid the costly generation of training data in advance. We show that FAB reaches higher sampling efficiency with fewer target evaluations in high dimensions in comparison to other methods.

We explore the potential for improving constraints on gravity by leveraging correlations in the dispersion measure derived from Fast Radio Bursts (FRBs) in combination with cosmic shear. Specifically, we focus on Horndeski gravity, inferring the kinetic braiding and Planck mass run rate from a stage-4 cosmic shear mock survey alongside a survey comprising 104 FRBs. For the inference pipeline, we utilise the Boltzmann code hi_class to predict the linear matter power spectrum in modified gravity scenarios, while non-linear corrections are obtained from the halo-model employed in HMcode, including feedback mechanisms. Our findings indicate that FRBs can disentangle degeneracies between baryonic feedback and cosmological parameters, as well as the mass of massive neutrinos. Since these parameters are also degenerate with modified gravity parameters, the inclusion of FRBs can enhance constraints on Horndeski parameters by up to 40 percent, despite being a less significant measurement. Additionally, we apply our model to current FRB data and use the uncertainty in the DM−z relation to impose limits on gravity. However, due to the limited sample size of current data, constraints are predominantly influenced by theoretical priors. Despite this, our study demonstrates that FRBs will significantly augment the limited set of cosmological probes available, playing a critical role in providing alternative tests of feedback, cosmology, and gravity. All codes used in this work are made publically available.

We calculate two-loop renormalization group equations (RGEs) in the Standard Model Effective Field Theory (SMEFT) with right-handed neutrinos, i.e., the so-called νSMEFT, up to dimension five. Besides the two-loop RGEs of dimension-five (dim-5) operators, we also present those of the renormalizable couplings, including contributions from dim-5 operators. We check consistency relations among the first and second poles of ɛ ≡ (4 - d)/2 with d being the space-time dimension for all renormalization constants and find that those for lepton doublet and right-handed neutrino wave-function renormalization constants, as well as for renormalization constants of charged-lepton and neutrino Yukawa coupling matrices, do not hold. This leads to divergent RG functions for these fields and Yuwawa coupling matrices. We figure out that such infinite RG functions arise from the non-invariance of fields and Yukawa coupling matrices under field redefinitions, considering that flavor transformations are a kind of linear field redefinitions. Those infinite RG functions will disappear once one restores contributions from the derivative of renormalization constants with respect to the Wilson coefficients of redundant operators or, alternatively, considers the RGEs of flavor invariants, which are physical quantities and remain invariant under field redefinitions.

Heavy-ion storage rings have relatively large momentum acceptance which allows for multiple ion species to circulate at the same time. This needs to be considered in radioactive decay measurements of highly charged ions, where atomic charge exchange reactions can significantly alter the intensities of parent and daughter ions. In this study, we investigate this effect using the decay curves of ion numbers in the recent <inline-formula><mml:math><mml:mmultiscripts><mml:mrow></mml:mrow><mml:mrow></mml:mrow><mml:mn>205</mml:mn></mml:mmultiscripts></mml:math></inline-formula>Tl<inline-formula><mml:math><mml:mmultiscripts><mml:mrow></mml:mrow><mml:mrow></mml:mrow><mml:mrow><mml:mn>81</mml:mn><mml:mo>+</mml:mo></mml:mrow></mml:mmultiscripts></mml:math></inline-formula> bound-state beta decay experiment conducted using the Experimental Storage Ring at GSI Darmstadt. To understand the intricate dynamics of ion numbers, we present a set of differential equations that account for various atomic and nuclear reaction processes—bound-state beta decay, atomic electron recombination and capture, and electron ionization. By incorporating appropriate boundary conditions, we develop a set of differential equations that accurately simulate the decay curves of various simultaneously stored ions in the storage ring: <inline-formula><mml:math><mml:mmultiscripts><mml:mrow></mml:mrow><mml:mrow></mml:mrow><mml:mn>205</mml:mn></mml:mmultiscripts></mml:math></inline-formula>Tl<inline-formula><mml:math><mml:mmultiscripts><mml:mrow></mml:mrow><mml:mrow></mml:mrow><mml:mrow><mml:mn>81</mml:mn><mml:mo>+</mml:mo></mml:mrow></mml:mmultiscripts></mml:math></inline-formula>, <inline-formula><mml:math><mml:mmultiscripts><mml:mrow></mml:mrow><mml:mrow></mml:mrow><mml:mn>205</mml:mn></mml:mmultiscripts></mml:math></inline-formula>Pb<inline-formula><mml:math><mml:mmultiscripts><mml:mrow></mml:mrow><mml:mrow></mml:mrow><mml:mrow><mml:mn>81</mml:mn><mml:mo>+</mml:mo></mml:mrow></mml:mmultiscripts></mml:math></inline-formula>, <inline-formula><mml:math><mml:mmultiscripts><mml:mrow></mml:mrow><mml:mrow></mml:mrow><mml:mn>205</mml:mn></mml:mmultiscripts></mml:math></inline-formula>Pb<inline-formula><mml:math><mml:mmultiscripts><mml:mrow></mml:mrow><mml:mrow></mml:mrow><mml:mrow><mml:mn>82</mml:mn><mml:mo>+</mml:mo></mml:mrow></mml:mmultiscripts></mml:math></inline-formula>, <inline-formula><mml:math><mml:mmultiscripts><mml:mrow></mml:mrow><mml:mrow></mml:mrow><mml:mn>200</mml:mn></mml:mmultiscripts></mml:math></inline-formula>Hg<inline-formula><mml:math><mml:mmultiscripts><mml:mrow></mml:mrow><mml:mrow></mml:mrow><mml:mrow><mml:mn>79</mml:mn><mml:mo>+</mml:mo></mml:mrow></mml:mmultiscripts></mml:math></inline-formula>, and <inline-formula><mml:math><mml:mmultiscripts><mml:mrow></mml:mrow><mml:mrow></mml:mrow><mml:mn>200</mml:mn></mml:mmultiscripts></mml:math></inline-formula>Hg<inline-formula><mml:math><mml:mmultiscripts><mml:mrow></mml:mrow><mml:mrow></mml:mrow><mml:mrow><mml:mn>80</mml:mn><mml:mo>+</mml:mo></mml:mrow></mml:mmultiscripts></mml:math></inline-formula>. Through a quantitative comparison between simulations and experimental data, we provide insights into the detailed reaction mechanisms governing stored heavy ions within the storage ring. Our approach effectively models charge-changing processes, reduces the complexity of the experimental setup, and provides a simpler method for measuring the decay half-lives of highly charged ions in storage rings.

We study the femtoscopic correlation functions of meson-baryon pairs in the strangeness <inline-formula><mml:math><mml:mi>S</mml:mi><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:math></inline-formula> sector, employing unitarized s-wave scattering amplitudes derived from the chiral Lagrangian up to next-to-leading order. For the first time, we deliver predictions on the <inline-formula><mml:math><mml:msup><mml:mi>π</mml:mi><mml:mo>-</mml:mo></mml:msup><mml:mi>Λ</mml:mi></mml:math></inline-formula> and <inline-formula><mml:math><mml:msup><mml:mi>K</mml:mi><mml:mo>+</mml:mo></mml:msup><mml:msup><mml:mi>Ξ</mml:mi><mml:mo>-</mml:mo></mml:msup></mml:math></inline-formula> correlation functions which are feasible to be measured at the Large Hadron Collider. We also demonstrate that the employed model is perfectly capable of reproducing the <inline-formula><mml:math><mml:msup><mml:mi>K</mml:mi><mml:mo>-</mml:mo></mml:msup><mml:mi>p</mml:mi></mml:math></inline-formula> correlation function data measured by the same collaboration, without the need to modify the coupling strength to the <inline-formula><mml:math><mml:msup><mml:mover><mml:mi>K</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mn>0</mml:mn></mml:msup><mml:mi>n</mml:mi></mml:math></inline-formula> channel, as has been recently suggested. In all cases, the effects of the source size on the correlation are tested. In addition, we present detailed analysis of the different coupled-channel contributions, together with the quantification of the relative relevance of the different terms in the interaction. These calculations require the knowledge of the so-called production weights, for which we present two novel methods to compute them.

Context. Heavy chemical elements such as iron in the intra-cluster medium (ICM) of galaxy clusters are a signpost of the interaction between the gas and stellar components. Observations of the ICM metallicity in present-day massive systems, however, pose a challenge to the underlying assumption that the cluster galaxies have produced the amount of iron that enriches the ICM. Aims. We evaluate the iron share between ICM and stars within simulated galaxy clusters with the twofold aim of investigating the origin of possible differences with respect to observational findings and of shedding light on the observed excess of iron on the ICM with respect to expectations based on the observed stellar population. Methods. We evaluated the iron mass in gas and stars in a sample of 448 simulated systems with masses M_tot,500>10¹⁴ M_⊙ at z = 0.07. These were extracted from the high-resolution (352 h^‑1 cMpc)³ volume of the MAGNETICUM cosmological hydrodynamical simulations. We compared our results with observational data of low-redshift galaxy clusters. Results. The iron share in simulated clusters features a shallow dependence on the total mass, and its value is close to unity on average. In the most massive simulated systems, the iron share is thus smaller than observational values by almost an order of magnitude. The dominant contribution to this difference is related to the stellar component, whereas the chemical properties of the ICM agree well overall with the observations. We find larger stellar mass fractions in simulated massive clusters, which in turn yield higher stellar iron masses, than in observational data. Conclusions. Consistently with the modelling, we confirm that the stellar content within simulated present-day massive systems causes the metal enrichment in the ICM. It will be crucial to alleviate the stellar mass discrepancy between simulations and observations to definitely assess the iron budget in galaxy clusters.

We present one of the largest uniform optical spectroscopic surveys of X-ray selected sources to date that were observed as a pilot study for the Black Hole Mapper (BHM) survey. The BHM program of the Sloan Digital Sky Survey (SDSS)-V is designed to provide optical spectra for hundreds of thousands of X-ray selected sources from the SRG/eROSITA all-sky survey. This significantly improves our ability to classify and characterise the physical properties of large statistical populations of X-ray emitting objects. Our sample consists of 13 079 sources in the eROSITA eFEDS performance verification field, 12 011 of which provide reliable redshifts from 0 ≲ z ≤ 5.8. The vast majority of these objects were detected as point-like sources (X-ray flux limit F_{0.5 ‑ 2 keV} ≳ 6.5 × 10^‑15 erg/s/cm²) and were observed for about 20 years with fibre-fed SDSS spectrographs. After including all available redshift information for the eFEDS sources from the dedicated SDSS-V plate programme and archival data, we visually inspected the SDSS optical spectra to verify the reliability of these redshift measurements and the performance of the SDSS pipeline. The visual inspection allowed us to recover reliable redshifts (for 99% of the spectra with a signal-to-noise ratio of > 2) and to assign classes to the sources, and we confirm that the vast majority of our sample consists of active galactic nuclei (AGNs). Only ∼3% of the eFEDS/SDSS sources are Galactic objects. We analysed the completeness and purity of the spectroscopic redshift catalogue, in which the spectroscopic completeness increases from 48% (full sample) to 81% for a cleaner, brighter (r_AB < 21.38) sample that we defined by considering a high X-ray detection likelihood, a reliable counterpart association, and an optimal sky coverage. We also show the diversity of the optical spectra of the X-ray selected AGNs and provide spectral stacks with a high signal-to-noise ratio in various sub-samples with different redshift and optical broad-band colours. Our AGN sample contains optical spectra of (broad-line) quasars, narrow-line galaxies, and optically passive galaxies. It is considerably diverse in its colours and in its levels of nuclear obscuration.

The ultra-hot Jupiter (UHJ) TOI-2109b marks the lower edge of the equilibrium temperature gap between 3500 and 4500 K, an unexplored thermal regime that separates KELT-9b, the hottest planet yet discovered, from all other currently known gas giants. To study the thermochemical structure of TOI-2109b's atmosphere, we obtained high-resolution emission spectra of both the planetary day- and nightsides with CAHA/CARMENES and VLT/CRIRES⁺. By applying the cross-correlation technique to the high-resolution spectra, we identified the emission signatures of Fe I (S/N = 4.3) and CO (S/N = 6.3), as well as a thermal inversion layer in the dayside atmo-sphere; no significant H₂O signal was detected from the dayside. None of the analyzed species were detectable from the nightside atmosphere. We applied a Bayesian retrieval framework that combines high-resolution spectroscopy with photometric measurements to constrain the dayside atmospheric parameters and derive upper limits for the nightside hemisphere. The dayside thermal inversion extends from approximately 3200 to 4600 K, with an atmospheric metallicity consistent with that of the host star (0.36 dex). Only weak constraints could be placed on the C/O ratio, with a lower limit of 0.15. The retrieved spectral line broadening is consistent with tidally locked rotation, indicating the absence of strong dynamical processes in the atmosphere. An upper temperature limit of approximately 2400 K and a maximum atmospheric temperature gradient of about 700 K/log bar could be derived for the planetary nightside. Comparison of the retrieved dayside temperature-pressure profile with theoretical models, the absence of strong atmospheric dynamics, and significant differences in the thermal constraints between the day- and nightside hemispheres suggest a limited heat transport efficiency across the planetary atmosphere. Overall, our results place TOI-2109b in a transitional regime between the UHJs below the thermal gap, which show both CO and H₂O emission lines, and KELT-9b, where molecular features are largely absent.

We present GalSBI, a phenomenological model of the galaxy population for cosmological applications using simulation-based inference. The model is based on analytical parametrizations of galaxy luminosity functions, morphologies and spectral energy distributions. Model constraints are derived through iterative Approximate Bayesian Computation, by comparing Hyper Suprime-Cam deep field images with simulations which include a forward model of instrumental, observational and source extraction effects. We developed an emulator trained on image simulations using a normalizing flow. We use it to accelerate the inference by predicting detection probabilities, including blending effects and photometric properties of each object, while accounting for background and PSF variations. This enables robustness tests for all elements of the forward model and the inference. The model demonstrates excellent performance when comparing photometric properties from simulations with observed imaging data for key parameters such as magnitudes, colors and sizes. The redshift distribution of simulated galaxies agrees well with high-precision photometric redshifts in the COSMOS field within 1.5σ for all magnitude cuts. Additionally, we demonstrate how GalSBI's redshifts can be utilized for splitting galaxy catalogs into tomographic bins, highlighting its potential for current and upcoming surveys. GalSBI is fully open-source, with the accompanying Python package, galsbi (https://cosmo-docs.phys.ethz.ch/galsbi/), offering an easy interface to quickly generate realistic, survey-independent galaxy catalogs.

We present magnetohydrodynamic simulations of star formation in the multiphase interstellar medium (ISM) to quantify the impact of non-ionizing far-ultraviolet (FUV) radiation within the SILCC PROJECT simulation framework. Our study incorporates the radiative transfer of ionizing radiation and self-consistent modelling of variable FUV radiation from star clusters, advancing beyond previous studies using static or simplified FUV fields. This enables a more accurate capture of the dynamic interaction between radiation and the evolving ISM alongside other stellar feedback channels. The interstellar radiation field (ISRF) near young star clusters can reach <inline-formula><tex-math>$G_0 \approx 10^4$</tex-math></inline-formula> (in Habing units), far exceeding the solar neighbourhood value of <inline-formula><tex-math>$G_0 = 1.7$</tex-math></inline-formula>. Despite these high intensities, FUV radiation minimally impacts the integrated star formation rate compared to ionizing radiation, stellar winds, and supernovae. A slight reduction in star formation burstiness is linked to increased photoelectric (PE) heating efficiency by the variable FUV field. Dust near star-forming regions can be heated up to 60 K via the PE effect, with a broad temperature distribution. PE heating rates in variable FUV models exhibit higher peaks but lower averages than static ISRF models. Simulations under solar neighbourhood conditions without stellar winds or ionizing radiation but with supernovae yield unexpectedly high star formation rates of <inline-formula><tex-math>$\sim 0.1~\mathrm{M_\odot ~yr^{-1}~kpc^{-2}}$</tex-math></inline-formula>. Our analysis reveals increased cold neutral medium volume-filling factors (VFF) outside stellar clusters, reduced thermally unstable gas, and sharper warm–cold gas separation. The variable FUV field also promotes a cold diffuse gas phase with a molecular component, exhibiting a VFF of <inline-formula><tex-math>$\sim 5{-}10$</tex-math></inline-formula> per cent.

Slow flavor evolution (defined as driven by neutrino masses and not necessarily "slow") is receiving fresh attention in the context of compact astrophysical environments. In Part I of this series, we have studied the slow-mode dispersion relation following our recently developed analogy to plasma waves. The concept of resonance between flavor waves in the linear regime and propagating neutrinos is the defining feature of this approach. It is best motivated for weak instabilities, which probably is the most relevant regime in self-consistent astrophysical environments because these will try to eliminate the cause of instability. We here go beyond the dispersion relation alone (which by definition applies to infinite media) and consider the group velocities of unstable modes that determines whether the instability relaxes within the region where it first appears (absolute), or away from it (convective). We show that all weak instabilities are convective so that their further evolution is not local. Therefore, studying their consequences numerically in small boxes from given initial conditions may not always be appropriate.

We investigate the phases of a strongly coupled QCD-like theory at finite baryon chemical potential using s-confining supersymmetric QCD deformed by anomaly-mediated supersymmetry breaking. Focusing on the case of three colors and four flavors, we identify novel phases including spontaneous breaking of baryon number and/or parity. Both first-order and second-order phase transitions are observed as the baryon chemical potential is varied. These findings may offer insights into possible phases of real QCD at intermediate baryon densities.

Primordial black holes (PBHs) with masses below $10^9\,\rm{g}$ are typically assumed to have negligible cosmological impact due to their rapid evaporation via Hawking radiation. However, the 'memory burden' effect, which is a quantum suppression of PBH evaporation, can dramatically alter their decay dynamics. In this work, we revisit early-Universe constraints on ultralight PBHs in this mass range, demonstrating that memory burden significantly alters previous constraints. We compute new cosmological bounds from BBN that strongly limit the presence of ultralight PBHs in the early Universe. We report that the PBHs in the mass range $10^0$-$10^2\,\rm{g}$ for $k=2$ are unconstrained by observations.

The properties of satellite halos provide a promising probe for dark matter (DM) physics. Observations motivate current efforts to explain surprisingly compact DM halos. If DM is not collisionless but has strong self-interactions, halos can undergo gravothermal collapse, leading to higher densities in the central region of the halo. However, it is challenging to model this collapse phase from first principles. To improve on this, we seek to better understand numerical challenges and convergence properties of self-interacting dark matter (SIDM) N-body simulations in the collapse phase. Especially we aim for a better understanding of the evolution of satellite halos. To do so, we run SIDM N-body simulations of a low mass halo in isolation and within an external gravitational potential. The simulation setup is motivated by the perturber of the stellar stream GD-1. We find that the halo evolution is very sensitive to energy conservation errors, and a too large SIDM kernel size can artificially speed up the collapse. Moreover, we demonstrate that the King model can describe the density profile at small radii for the late stages that we have simulated. Furthermore, for our highest-resolved simulation (N = 5x10^7) we make the data public. It can serve as a benchmark. Overall, we find that the current numerical methods do not suffer from convergence problems in the late collapse phase and provide guidance on how to choose numerical parameters, e.g. that the energy conservation error is better kept well below 1%. This allows to run simulations of halos becoming concentrated enough to explain observations of GD-1 like stellar streams or strong gravitational lensing systems.

The precision measurement of the tritium <inline-formula><mml:math><mml:mrow><mml:mtext>β</mml:mtext></mml:mrow></mml:math></inline-formula>-decay spectrum performed by the KATRIN experiment provides a unique way to search for general neutrino interactions (GNIs). All theoretically allowed GNI terms at dimension 6 involving neutrinos are incorporated into a low-energy effective field theory, and can be identified by specific signatures in the measured tritium <inline-formula><mml:math><mml:mrow><mml:mtext>β</mml:mtext></mml:mrow></mml:math></inline-formula> spectrum. In this Letter an effective description of the impact of GNIs on the <inline-formula><mml:math><mml:mrow><mml:mtext>β</mml:mtext></mml:mrow></mml:math></inline-formula> spectrum is formulated and the first constraints on the effective GNI parameters are derived based on the <inline-formula><mml:math><mml:mn>4</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn>10</mml:mn><mml:mn>6</mml:mn></mml:msup></mml:math></inline-formula> electrons collected in the second measurement campaign of KATRIN in 2019. In addition, constraints on selected types of interactions are investigated, thereby exploring the potential of KATRIN to search for more specific new physics cases, including a right-handed W boson, a charged Higgs boson, or leptoquarks.

We describe the recent experiments which claimed an observation of a tetra-neutron signal. Production reactions like transfer, knockout, fragmentation or photodisintegration have been used at very different experiments and facilities to form systems just made of neutrons. As a possible explanation of the partly contradicting results we suggest that some observed the bound ground state and some an unbound but still correlated state of the four neutrons at different exitation energy. We also refer to some of the theoretical works.

The current generation of cryogenic solid state detectors used in direct dark matter and CE\textnu NS searches typically reach energy thresholds of $\mathcal{O}$(10)$\,$eV for nuclear recoils. For a reliable calibration in this energy regime a method has been proposed, providing mono-energetic nuclear recoils at low energies $\sim\,$100$\,$eV$\,$-$\,$1$\,$keV. In this work we report on the observation of a peak at (1113.6$^{+6.5}_{-6.5}$)$\,$eV in the data of an Al$_{2}$O$_{3}$ crystal in CRESST-III, which was irradiated with neutrons from an AmBe calibration source. We attribute this mono-energetic peak to the radiative capture of thermal neutrons on $^{27}$Al and the subsequent de-excitation via single $γ$-emission. We compare the measured results with the outcome of Geant4 simulations and investigate the possibility to make use of this effect for the energy calibration of Al$_{2}$O$_{3}$ detectors at low energies. We further investigate the possibility of a shift in the expected energy scale of this effect caused by the creation of defects in the target crystal.

We present cosmological constraints from 8 strongly lensed quasars (hereafter, the TDCOSMO-2025 sample). Building on previous work, our analysis incorporated new deflector stellar velocity dispersions measured from spectra obtained with the James Webb Space Telescope (JWST), the Keck Telescopes, and the Very Large Telescope (VLT), utilizing improved methods. We used integrated JWST stellar kinematics for 5 lenses, VLT-MUSE for 2, and resolved kinematics from Keck and JWST for RX J1131-1231. We also considered two samples of non-time-delay lenses: 11 from the Sloan Lens ACS (SLACS) sample with Keck-KCWI resolved kinematics; and 4 from the Strong Lenses in the Legacy Survey (SL2S) sample. We improved our analysis of line-of-sight effects, the surface brightness profile of the lens galaxies, and orbital anisotropy, and corrected for projection effects in the dynamics. Our uncertainties are maximally conservative by accounting for the mass-sheet degeneracy in the deflectors' mass density profiles. The analysis was blinded to prevent experimenter bias. Our primary result is based on the TDCOSMO-2025 sample, in combination with $Ω_{\rm m}$ constraints from the Pantheon+ Type Ia supernovae (SN) dataset. In the flat $Λ$ Cold Dark Matter (CDM), we find $H_0=71.6^{+3.9}_{-3.3}$ km s$^{-1}$ Mpc$^{-1}$. The SLACS and SL2S samples are in excellent agreement with the TDCOSMO-2025 sample, improving the precision on $H_0$ in flat $Λ$CDM to 4.6%. Using the Dark Energy Survey SN Year-5 dataset (DES-SN5YR) or DESI-DR2 baryonic acoustic oscillations (BAO) likelihoods instead of Pantheon+ yields very similar results. We also present constraints in the open $Λ$CDM, $w$CDM, $w_0w_a$CDM, and $w_ϕ$CDM cosmologies. The TDCOSMO $H_0$ inference is robust and consistent across all presented cosmological models, and our cosmological constraints in them agree with those from the BAO and SN.

The dilaton, a pseudo-Nambu-Goldstone boson (pNGB) of broken scale invariance, is an appealing ultralight dark matter (DM) candidate. Its mass is protected by conformal invariance and it can be searched for in tabletop experiments. However, contrary to standard pNGBs of internal symmetries, the dilaton generically has a large non-derivative self-coupling, leading to radiative contributions to its mass of the order of its decay constant. Hence typical ultralight dilatons should also have sub-eV decay constants, which would incur significant deviations from standard DM behavior at structure formation times, in severe tension with observations. Therefore, a fine-tuning is required to generate a hierarchy between the mass and the decay constant. In this work, we consider whether supersymmetry (SUSY) can be used to protect this hierarchy from quantum corrections. To ensure an ultralight dilaton mass robust against realistic SUSY-breaking contributions, we must consider a novel dilaton stabilization mechanism. The observed DM abundance can be produced by the misalignment mechanism for dilaton masses ranging from $10^{-11}$ to $1$ eV. Unfortunately, irreducible SUSY-breaking corrections due to gravity restrict the couplings between the dilaton and the Standard Model to be extremely small, beyond the reach of any current or proposed experiments. Our work demonstrates that constructing a consistent model of ultralight dilaton DM is quite involved.

We calculate the $θ$ dependence in a cousin of QCD, where the vacuum structure can be analyzed exactly. The theory is $\mathcal{N}=2$ $SU(2)$ gauge theory with $N_F=0,1,2,3$ flavors of fundamentals, explicitly broken to $\mathcal{N}=1$ via an adjoint superpotential, and coupled to anomaly mediated supersymmetry breaking (AMSB). The hierarchy $m_{AMSB}\ll μ_{\mathcal{N}=1}\ll Λ$ ensures the validity of our IR analysis. As expected from ordinary QCD, the vacuum energy is a function of $θ$ which undergoes 1st order phase transitions between different vacua where the various dyons condense. For $N_F=0$ we find the expected phase transition at $θ=π$, while for $N_F=1,2,3$ we find phase transitions at fractional values of $π$.

Surprisingly compact substructures in galaxies and galaxy clusters, but also field halos, have been observed by gravitational lensing. They could be difficult to explain with collisionless dark matter (DM). To explain those objects, recent studies focused on the gravothermal collapse that halos consisting of self-interacting dark matter (SIDM) can undergo. However, simple models of elastic scattering could face problems explaining those compact objects during very later stages of the collapse and the post-collapse phase, where a black hole may have formed from DM. We aim to explain compact halos while avoiding the gravothermal catastrophe to which typical SIDM models are subject. Therefore, we investigate the evolution of a DM halo for an SIDM model consisting of two species with unequal masses, which features only interactions between the different species but not within themselves. Employing $N$-body simulations, we study the effect of unequal-mass SIDM models on the evolution of an isolated DM halo. In particular, the late stages of its evolution with high central densities are simulated. We find that our two-species SIDM models can produce density cores with their size depending on the mass ratio of the two species. Moreover, mass segregation caused by the unequal particle masses leads to a finite final density state or at least a slowly growing density, which depends on the mass ratio and the mass fraction of the two DM species. SIDM models consisting of two DM species can simultaneously explain DM halos with density cores, as well as systems that are denser in their centre than expected from collisionless DM, while avoiding the gravothermal catastrophe. They are a compelling alternative to single-species models, offering a rich phenomenology.

We argue that the Standard Model is accompanied by a new pseudoscalar degree of freedom, <inline-formula><mml:math><mml:msub><mml:mi>η</mml:mi><mml:mi>w</mml:mi></mml:msub></mml:math></inline-formula> meson, which cancels the topological susceptibility of the electroweak vacuum and gets its mass from this effect. The prediction is based on the analyticity properties of the Chern-Simons correlator combined with the basic features of gravity. Depending on the quality level of the <inline-formula><mml:math><mml:mrow><mml:mi>U</mml:mi><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:msub><mml:mrow><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mi>B</mml:mi><mml:mo>+</mml:mo><mml:mi>L</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> symmetry, the <inline-formula><mml:math><mml:msub><mml:mi>η</mml:mi><mml:mi>w</mml:mi></mml:msub></mml:math></inline-formula> emerges as a <inline-formula><mml:math><mml:mi>B</mml:mi><mml:mo>+</mml:mo><mml:mi>L</mml:mi></mml:math></inline-formula> pseudo-Goldstone boson or as a Stückelberg 2-form of the electroweak gauge redundancy. An intriguing scenario of the first category is the emergence of <inline-formula><mml:math><mml:msub><mml:mi>η</mml:mi><mml:mi>w</mml:mi></mml:msub></mml:math></inline-formula> in the form of the phase of a <inline-formula><mml:math><mml:mi>U</mml:mi><mml:mo>(</mml:mo><mml:mn>1</mml:mn><mml:msub><mml:mo>)</mml:mo><mml:mrow><mml:mi>B</mml:mi><mml:mo>+</mml:mo><mml:mi>L</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula>-violating fermion condensate triggered by the instantons, somewhat similar to the <inline-formula><mml:math><mml:msup><mml:mi>η</mml:mi><mml:mo>'</mml:mo></mml:msup></mml:math></inline-formula> meson in QCD. Regardless of its origin, the presence of the <inline-formula><mml:math><mml:msub><mml:mi>η</mml:mi><mml:mi>w</mml:mi></mml:msub></mml:math></inline-formula> meson in the theory appears to be a matter of consistency.

Nonlinearities in King plots (KP) of isotope shifts (IS) can reveal the existence of beyond-standard-model (BSM) interactions that couple electrons and neutrons. However, it is crucial to distinguish higher-order standard model (SM) effects from BSM physics. We measure the IS of the transitions <inline-formula><mml:math><mml:mrow><mml:msub><mml:mrow><mml:mmultiscripts><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:mmultiscripts></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>→</mml:mo><mml:msub><mml:mrow><mml:mmultiscripts><mml:mrow><mml:mi>P</mml:mi></mml:mrow><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:mmultiscripts></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> in <inline-formula><mml:math><mml:mrow><mml:msup><mml:mrow><mml:mi>Ca</mml:mi></mml:mrow><mml:mrow><mml:mn>14</mml:mn><mml:mo>+</mml:mo></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math><mml:mrow><mml:msub><mml:mrow><mml:mmultiscripts><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:mmultiscripts></mml:mrow><mml:mrow><mml:mn>1</mml:mn><mml:mo>/</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>→</mml:mo><mml:mmultiscripts><mml:mrow><mml:msub><mml:mrow><mml:mi>D</mml:mi></mml:mrow><mml:mrow><mml:mn>5</mml:mn><mml:mo>/</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:mmultiscripts></mml:mrow></mml:math></inline-formula> in <inline-formula><mml:math><mml:msup><mml:mi>Ca</mml:mi><mml:mo>+</mml:mo></mml:msup></mml:math></inline-formula> with sub-Hz precision as well as the nuclear mass ratios with relative uncertainties below <inline-formula><mml:math><mml:mn>4</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn>10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn>11</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula> for the five stable, even isotopes of calcium (<inline-formula><mml:math><mml:mrow><mml:mmultiscripts><mml:mrow><mml:mi>Ca</mml:mi></mml:mrow><mml:mrow><mml:mn>40</mml:mn><mml:mo>,</mml:mo><mml:mn>42</mml:mn><mml:mo>,</mml:mo><mml:mn>44</mml:mn><mml:mo>,</mml:mo><mml:mn>46</mml:mn><mml:mo>,</mml:mo><mml:mn>48</mml:mn></mml:mrow></mml:mmultiscripts></mml:mrow></mml:math></inline-formula>). Combined, these measurements yield a calcium KP nonlinearity with a significance of <inline-formula><mml:math><mml:mo>∼</mml:mo><mml:msup><mml:mn>10</mml:mn><mml:mn>3</mml:mn></mml:msup><mml:mi>σ</mml:mi></mml:math></inline-formula>. Precision calculations show that the nonlinearity cannot be fully accounted for by the expected largest higher-order SM effect, the second-order mass shift, and identify the little-studied nuclear polarization as the only remaining SM contribution that may be large enough to explain it. Despite the observed nonlinearity, we improve existing KP-based constraints on a hypothetical Yukawa interaction for most of the new boson masses between <inline-formula><mml:math><mml:mrow><mml:mn>10</mml:mn><mml:mtext> </mml:mtext><mml:mtext> </mml:mtext><mml:mi>eV</mml:mi><mml:mo>/</mml:mo><mml:msup><mml:mrow><mml:mi>c</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math><mml:mrow><mml:msup><mml:mrow><mml:mn>10</mml:mn></mml:mrow><mml:mrow><mml:mn>7</mml:mn></mml:mrow></mml:msup><mml:mtext> </mml:mtext><mml:mtext> </mml:mtext><mml:mi>eV</mml:mi><mml:mo>/</mml:mo><mml:msup><mml:mrow><mml:mi>c</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>.

We formulate the BCOV theory of deformations of complex structures as a pull-back to the super moduli space of the worldline of a spinning particle. In this approach the appearance of a non-local kinetic term in the target space action has the same origin as the mismatch of pictures in the Ramond sector of super string field theory and is resolved by the same type of auxiliary fields in shifted pictures. The BV-extension is manifest in this description. A compensator for the holomorphic 3-form can be included by resorting to a description in the large Hilbert space.

We present a first step toward field-level cosmological inference beyond the standard $Λ$CDM model, focusing on optimizing precision tests in the nonlinear regime of large-scale structure (LSS). As an illustrative case, we study the model-independent ``bootstrap'' coefficient of the second-order perturbation theory (PT) kernel for matter in real space, which we use as a proxy for new physics effects in the nonlinear sector. We discuss in details the ultraviolet (UV) cutoff dependence induced by discretizing fields on a grid, which requires proper renormalization to eliminate grid artifacts. We formulate a Wilsonian perturbative framework in which the evolution from a UV theory defined at a high cutoff $Λ_\text{uv}$ down to lower cutoffs is computed analytically, even beyond the validity of a derivative expansion. Within this framework, we develop an extended version of the GridSPT code incorporating the bootstrap parameterization and demonstrate how cutoff-independent predictions can be achieved through the inclusion of appropriate counterterms. We validate our approach at third- and fifth-order in PT, emphasizing the importance of higher-derivative contributions for unbiased parameter extraction. Our framework is readily extendable to biased tracers and redshift-space distortions.

Observations and high-resolution hydrodynamical simulations indicate that massive star clusters form through a complex hierarchical assembly. We use simulations including post-Newtonian dynamics (the BIFROST code) and stellar evolution (the SEVN module) to investigate this collisional assembly. With a full initial stellar mass function, we study the effect of initial binary, triple and massive single stars (450 $M_\odot$) on the assembly, structure, and kinematics of massive ($M_\mathrm{cl}\sim10^6 M_\odot$, $N=1.8 \times 10^6$) star clusters. Simultaneously, intermediate mass black holes (IMBHs), potential seeds for supermassive black holes, can form and grow in our models by stellar collisions, tidal disruption events (TDEs) and black hole (BH) mergers. At a fixed cluster mass, stellar multiplicity or a high mass limit increase the numbers (up to $\sim$ 10) and masses (up to $10^4 M_\odot$) of the formed IMBHs within the first 10 Myr of cluster evolution. The TDE rates peak at $Γ_\mathrm{tde}\sim 5 \times 10^{-5}$ yr$^{-1}$ after IMBH formation at $\sim 2$ Myr. In all simulations, we find gravitational wave driven mergers involving stellar BHs and IMBHs. Initial multiplicity or a high mass limit also result in IMBH-IMBH mergers. The IMBH masses correlate with the initial cluster masses, surface densities and velocity dispersions approximately as $M_\bullet \propto M_\mathrm{cl}$, $M_\bullet\proptoΣ_\mathrm{h}^\mathrm{3/2}$ and $M_\bullet\proptoσ^\mathrm{3}$. Our results suggest the dense $z\sim10$ star clusters recently observed by the James Webb Space Telescope host IMBHs with masses above $M_\bullet \gtrsim 10^4 M_\odot$.

Spatially resolved stellar kinematics has become a key ingredient in time-delay cosmography to break the mass-sheet degeneracy in the mass profile and, in turn, provide a precise constraint on the Hubble constant and other cosmological parameters. In this paper, we present the first measurements of 2D resolved stellar kinematics for the lens galaxy in the quadruply lensed quasar system \lensname, using integral field spectroscopy from the JWST's Near-Infrared Spectrograph (NIRSpec), marking the first such measurement conducted with the JWST. In extracting robust kinematic measurements from this first-of-its-kind dataset, we made methodological improvements both in the data reduction and kinematic extraction. In our kinematic extraction procedure, we performed joint modeling of the lens galaxy, the quasar, and its host galaxy's contributions in the spectra to deblend the lens galaxy component and robustly constrain its stellar kinematics. Our improved methodological frameworks are released as software pipelines for future use: \textsc{squirrel} for extracting stellar kinematics, and \textsc{RegalJumper} for JWST-NIRSpec data reduction, incorporating additional artifact cleaning beyond the standard JWST pipeline. We compared our measured stellar kinematics from the JWST NIRSpec with previously obtained ground-based measurements from the Keck Cosmic Web Imager integral field unit and find that the two datasets are statistically consistent at a $\sim$1.1$σ$ confidence level. Our measured kinematics will be used in a future study to improve the precision of the Hubble constant measurement.

In this work, we study the impact of an imperfect knowledge of the instrument bandpasses on the estimate of the tensor-to-scalar ratio $r$ in the context of the next-generation LiteBIRD satellite. We develop a pipeline to include bandpass integration in both the time-ordered data (TOD) and the map-making processing steps. We introduce the systematic effect by having a mismatch between the ``real'', high resolution bandpass $τ$, entering the TOD, and the estimated one $τ_s$, used in the map-making. We focus on two aspects: the effect of degrading the $τ_s$ resolution, and the addition of a Gaussian error $σ$ to $τ_s$. To reduce the computational load of the analysis, the two effects are explored separately, for three representative LiteBIRD channels (40 GHz, 140 GHz and 402 GHz) and for three bandpass shapes. Computing the amount of bias on $r$, $Δr$, caused by these effects on a single channel, we find that a resolution $\lesssim 1.5$ GHz and $σ\lesssim 0.0089$ do not exceed the LiteBIRD budget allocation per systematic effect, $Δr < 6.5 \times 10^{-6}$. We then check that propagating separately the uncertainties due to a resolution of 1 GHz and a measurement error with $σ= 0.0089$ in all LiteBIRD frequency channels, for the most pessimistic bandpass shape of the three considered, still produces a $Δr < 6.5 \times 10^{-6}$. This is done both with the simple deprojection approach and with a blind component separation technique, the Needlet Internal Linear Combination (NILC). Due to the effectiveness of NILC in cleaning the systematic residuals, we have tested that the requirement on $σ$ can be relaxed to $σ\lesssim 0.05$. (Abridged)

Strong gravitationally lensed supernovae (LSNe), though rare, are exceptionally valuable probes for cosmology and astrophysics. Upcoming time-domain surveys like the Vera Rubin Observatory's Legacy Survey of Space and Time (LSST) offer a major opportunity to discover them in large numbers. Early identification is crucial for timely follow-up observations. We develop a deep learning pipeline to detect LSNe using multi-band, multi-epoch image cutouts. Our model is based on a 2D convolutional long short-term memory (ConvLSTM2D) architecture, designed to capture both spatial and temporal correlations in time-series imaging data. Predictions are made after each observation in the time series, with accuracy improving as more data arrive. We train the model on realistic simulations derived from Hyper Suprime-Cam (HSC) data, which closely matches LSST in depth and filters. This work focuses exclusively on Type Ia supernovae (SNe Ia). LSNe Ia are injected onto HSC luminous red galaxies (LRGs) at various phases of evolution to create positive examples. Negative examples include variable sources from HSC Transient Survey (including unclassified transients), and simulated unlensed SNe Ia in LRG and spiral galaxies. Our multi-band model shows rapid classification improvements during the initial few observations and quickly reaches high detection efficiency: at a fixed false-positive rate (FPR) of $0.01\%$, the true-positive rate (TPR) reaches $\gtrsim 60\%$ by the 7th observation and exceeds $\gtrsim 70\%$ by the 9th. Among the negative examples, SNe in LRGs remain the primary source of FPR, as they can resemble their lensed counterparts under certain conditions. The model detects quads more effectively than doubles and performs better on systems with larger image separations. Although trained and tested on HSC-like data, our approach applies to any cadenced imaging survey, particularly LSST.

To give more credence to the M-theoretic Emergence Proposal it is important to show that also classical kinetic terms in a low energy effective action arise as a quantum effect from integrating out light towers of states. We show that for compactifications of type IIA on Calabi-Yau manifolds, the classical weak coupling Yukawa couplings, which are the triple intersection numbers of the Calabi-Yau threefold, can be obtained from the 1/2-BPS protected one-loop Schwinger integral over $D2$-$D0$ bound states, after employing a novel regularization for the final infinite sum of Gopakumar-Vafa invariants. Approaching the problem in a consecutive manner from 6D decompactification over emergent string to the ultimate M-theory limits, we arrive at a mathematically concrete regularization that involves finite distance degeneration limits of Calabi-Yau threefolds in an intriguing way. We test and challenge this proposal by the concrete determination of the periods around such degeneration points for threefolds with one Kähler modulus and the two examples $\mathbb P_{1,1,1,6,9}[18]$ and $\mathbb P_{1,1,2,2,6}[12]$.

We present COSMOS2025, the COSMOS-Web catalog of photometry, morphology, photometric redshifts and physical parameters for more than 700,000 galaxies in the Cosmic Evolution Survey (COSMOS) field. This catalog is based on our \textit{James Webb Space Telescope} 255\,h COSMOS-Web program, which provides deep near-infrared imaging in four NIRCam (F115W, F150W, F277W, F444W) and one MIRI (F770W) filter over the central $\sim 0.54 {\, \rm deg}^2$ ($\sim 0.2 {\, \rm deg}^2$ for MIRI) in COSMOS. These data are combined with ground- and space-based data to derive photometric measurements of NIRCam-detected sources using both fixed-aperture photometry (on the space-based bands) and a profile-fitting technique on all 37 bands spanning 0.3-8 micron. We provide morphology for all sources from complementary techniques including profile fitting and machine-learning classification. We derive photometric redshifts, physical parameters and non-parametric star formation histories from spectral energy distribution (SED) fitting. The catalog has been extensively validated against previous COSMOS catalogs and other surveys. Photometric redshift accuracy measured using spectroscopically confirmed galaxies out to $z\sim9$ reaches $σ_{\rm MAD} = 0.012$ at $m_{\rm F444W}<28$ and remains at $σ_{\rm MAD} \lesssim 0.03$ as a function of magnitude, color, and galaxy type. This represents a factor of $\sim 2$ improvement at 26 AB mag compared to COSMOS2020. The catalog is approximately 80\% complete at $\log(M_{\star}/{\rm M}_{\odot}) \sim 9$ at $z \sim 10$ and at $\log(M_{\star}/{\rm M}_{\odot}) \sim 7$ at $z \sim 0.2$, representing a gain of 1\,dex compared to COSMOS2020. COSMOS2025 represents the definitive COSMOS-Web catalog. It is provided with complete documentation, together with redshift probability distributions, and it is ready for scientific exploitation today.

A commonly employed method to detect protoclusters in the young universe is the search for overdensities of massive star forming galaxies, such as submillimeter galaxies (SMGs), around high-mass halos, including those hosting quasars. In this work, we study the Megaparsec environment surrounding nine physically associated quasar pairs between $z=2.45$ and $z=3.82$ with JCMT/SCUBA-2 observations at 450 $μ$m and 850 $μ$m covering a field of view of roughly 13.7 arcmin in diameter (or 32 Mpc$^2$ at the median redshift) for each system. We identify a total of 170 SMG candidates and 26 non-SMG and interloper candidates. A comparison of the underlying 850 $μ$m source models recovered with Monte Carlo simulations to the blank field model reveals galaxy overdensities in all fields, with a weighted average overdensity factor of $δ_{\rm cumul} = 3.4 \pm 0.3$. From this excess emission at 850 $μ$m, we calculate a star formation rate density of $1700 \pm 100$ M$_{\odot}$ yr$^{-1}$ Mpc$^{-3}$, consistent with predictions from protocluster simulations and observations. Compared to fields around single quasars, those surrounding quasar pairs have higher excess counts and more centrally peaked star formation, further highlighting the co-evolution of SMGs and quasars. We do not find preferential alignment of the SMGs with the quasar pair direction or their associated Ly$α$ nebulae, indicating that cosmic web filaments on different scales might be traced by the different directions. Overall, this work substantiates the reliability of quasar pairs to detect overdensities of massive galaxies and likely sites of protocluster formation. Future spectroscopic follow-up observations are needed to confirm membership of the SMG candidates with the physically associated quasar pairs and definitively identify the targeted fields as protoclusters.

We develop a framework to study the relation between the stellar mass of a galaxy and the total mass of its host dark matter halo using galaxy clustering and galaxy-galaxy lensing measurements. We model a wide range of scales, roughly from $\sim 100 \; {\rm kpc}$ to $\sim 100 \; {\rm Mpc}$, using a theoretical framework based on the Halo Occupation Distribution and data from Year 3 of the Dark Energy Survey (DES) dataset. The new advances of this work include: 1) the generation and validation of a new stellar mass-selected galaxy sample in the range of $\log M_\star/M_\odot \sim 9.6$ to $\sim 11.5$; 2) the joint-modeling framework of galaxy clustering and galaxy-galaxy lensing that is able to describe our stellar mass-selected sample deep into the 1-halo regime; and 3) stellar-to-halo mass relation (SHMR) constraints from this dataset. In general, our SHMR constraints agree well with existing literature with various weak lensing measurements. We constrain the free parameters in the SHMR functional form $\log M_\star (M_h) = \log(εM_1) + f\left[ \log\left( M_h / M_1 \right) \right] - f(0)$, with $f(x) \equiv -\log(10^{αx}+1) + δ[\log(1+\exp(x))]^γ/ [1+\exp(10^{-x})]$, to be $\log M_1 = 11.559^{+0.334}_{-0.415}$, $\log ε= -1.689^{+0.333}_{-0.220}$, $α= -1.637^{+0.107}_{-0.096}$, $γ= 0.588^{+0.265}_{-0.220}$ and $δ= 4.227^{+2.223}_{-1.776}$. The inferred average satellite fraction is within $\sim 5-35\%$ for our fiducial results and we do not see any clear trends with redshift or stellar mass. Furthermore, we find that the inferred average galaxy bias values follow the generally expected trends with stellar mass and redshift. Our study is the first SHMR in DES in this mass range, and we expect the stellar mass sample to be of general interest for other science cases.

Primordial black holes (PBHs) with masses below $10^9\,\rm{g}$ are typically assumed to have negligible cosmological impact due to their rapid evaporation via Hawking radiation. However, the 'memory burden' effect, which is a quantum suppression of PBH evaporation, can dramatically alter their decay dynamics. In this work, we revisit early-Universe constraints on ultralight PBHs in this mass range, demonstrating that memory burden significantly alters previous constraints. We compute new cosmological bounds from BBN that strongly limit the presence of ultralight PBHs in the early Universe. We report that the PBHs in the mass range $10^0$-$10^2\,\rm{g}$ for $k=2$ are unconstrained by observations.

Quantum simulators offer great potential for investigating dynamical properties of quantum field theories. However, preparing accurate non-trivial initial states for these simulations is challenging. Classical Euclidean-time Monte-Carlo methods provide a wealth of information about states of interest to quantum simulations. Thus, it is desirable to facilitate state preparation on quantum simulators using this information. To this end, we present a fully classical pipeline for generating efficient quantum circuits for preparing the ground state of an interacting scalar field theory in 1+1 dimensions. The first element of this pipeline is a variational ansatz family based on the stellar hierarchy for bosonic quantum systems. The second element of this pipeline is the classical moment-optimization procedure that augments the standard variational energy minimization by penalizing deviations in selected sets of ground-state correlation functions (i.e., moments). The values of ground-state moments are sourced from classical Euclidean methods. The resulting states yield comparable ground-state energy estimates but exhibit distinct correlations and local non-Gaussianity. The third element of this pipeline is translating the moment-optimized ansatz into an efficient quantum circuit with an asymptotic cost that is polynomial in system size. This work opens the way to systematically applying classically obtained knowledge of states to prepare accurate initial states in quantum field theories of interest in nature.

We present an interface between PineAPPL and Matrix, which allows fully differential cross sections to be calculated in the form of interpolation grids, accurate at next-to-next-to-leading order (NNLO) in QCD and next-to-leading order in electroweak (EW) theory. This interface is the first publicly available tool to calculate interpolation grids at NNLO QCD accuracy for a wide set of processes. Interpolation grids provide the functionality to compute predictions for arbitrary parton distribution functions (PDFs) as well as PDF uncertainties without the need to repeat the actual calculation. Another important application of the these grids is to perform global analyses of PDFs using exact NNLO calculations instead of $K$-factors, which have several drawbacks. This exact treatment of NNLO corrections is also an important prerequisite for fitting PDFs at next-to-next-to-next-to-leading order level with reliable uncertainties. The new version of the Matrix code interfaced to PineAPPL, as well as the grids produced for this publication, are available on the Matrix website and on PloughShare, respectively.

The mass accretion rate of galaxy clusters is a key factor in determining their structure, but a reliable observational tracer has yet to be established. We present a state-of-the-art machine learning model for constraining the mass accretion rate of galaxy clusters from only X-ray and thermal Sunyaev–Zeldovich observations. Using idealized mock observations of galaxy clusters from the MillenniumTNG simulation, we train a machine learning model to estimate the mass accretion rate. The model constrains 68% of the mass accretion rates of the clusters in our data set to within 33% of the true value without significant bias, a ∼58% reduction in the scatter over existing constraints. We demonstrate that the model uses information from both radial surface brightness density profiles and asymmetries.

We develop a new approach to Vlasov Perturbation Theory (VPT) that solves for the hierarchy of cumulants of the phase-space distribution function to arbitrarily high truncation order in the context of cosmological structure formation driven by collisionless dark matter. We investigate the impact of higher cumulants on density and velocity power spectra as well as the bispectrum, and compare to scale-free <inline-formula><mml:math><mml:mi>N</mml:mi></mml:math></inline-formula>-body simulations. While there is a strong difference between truncation at the first cumulant, i.e., standard perturbation theory (SPT), and truncation at the second (i.e., including the velocity dispersion tensor), the third cumulant has a small quantitative impact and fourth and higher cumulants only have a minor effect on these summary statistics at weakly non-linear scales. We show that spurious exponential growth is absent in vector and tensor modes if scalar-mode constraints on the non-Gaussianity of the background distribution function that results from shell crossing are satisfied, guaranteeing the screening of UV modes for all fluctuations of any type, as expected physically. We also show analytically that loop corrections to the power spectrum are finite within VPT for any initial power spectra consistent with hierarchical clustering, unlike SPT. Finally, we discuss the relation to and contrast our predictions with effective field theory (EFT), and discuss how the advantages of VPT and EFT approaches could be combined.

Due to their inherent capabilities of capturing non-local dependencies, Transformer neural networks have quickly been established as the paradigmatic architecture for large language models and image processing. Next to these traditional applications, machine learning (ML) methods have also been demonstrated to be versatile tools in the analysis of image-like data of quantum phases of matter, e.g. given snapshots of many-body wave functions obtained in ultracold atom experiments. While local correlation structures in image-like data of physical systems can reliably be detected, identifying phases of matter characterized by global, non-local structures with interpretable ML methods remains a challenge. Here, we introduce the correlator Transformer (CoTra), which classifies different phases of matter while at the same time yielding full interpretability in terms of physical correlation functions. The network's underlying structure is a tailored attention mechanism, which learns efficient ways to weigh local and non-local correlations for a successful classification. We demonstrate the versatility of the CoTra by detecting local order in the Heisenberg antiferromagnet, and show that local gauge constraints in one- and two-dimensional lattice gauge theories can be identified. Furthermore, we establish that the CoTra reliably detects non-local structures in images of correlated fermions in momentum space (Cooper pairs) and that it can distinguish percolating from non-percolating images.

We apply the framework of Vlasov Perturbation Theory (VPT) to the two-loop matter power spectrum within $Λ$CDM cosmologies. The main difference to Standard Perturbation Theory (SPT) arises from taking the velocity dispersion tensor into account, and the resulting screening of the backreaction of UV modes renders loop integrals cutoff-independent. VPT is informed about non-perturbative small scale dynamics via the average value of the dispersion generated by shell-crossing, which impacts the evolution of perturbations on weakly non-linear scales. When using an average dispersion from halo models, the VPT power spectrum agrees with the one from the simulation, up to differences from missing three-loop contributions. Alternatively, treating the average dispersion as free parameter we find a remarkably stable prediction of the matter power spectrum from collisionless dynamics at percent level for a wide range of the dispersion scale. We quantify the impact of truncating the Vlasov hierarchy for the cumulants of the phase-space distribution function, finding that the two-loop matter power spectrum is robust to neglecting third and higher cumulants. Finally, we introduce and validate a simplified fast scheme fVPT that can be easily incorporated into existing codes and is as numerically efficient as SPT.

Most young stars and therefore planetary systems form in high-mass star forming regions and are exposed to ultraviolet radiation, affecting the protoplanetary disk. These regions are located at large distances and only now with JWST become accessible to study the inner disks surrounding young stars. We present the eXtreme UV Environments (XUE) program, which provides the first detailed characterization of the physical and chemical properties of the inner disks around young intermediate-mass stars exposed to external irradiation from nearby massive stars. We present high signal to noise MIRI-MRS spectroscopy of 12 disks located in three sub-clusters of the high-mass star-forming region NGC 6357. Based on their mid-infrared spectral energy distribution, we classify the XUE sources into Group I and II based on the Meeus scheme. We analyze their molecular emission features, and compare their spectral indices and 10 $μ$m silicate emission profiles to those of nearby Herbig and intermediate T Tauri disks. Despite being more massive, the XUE stars host disks with molecular richness comparable to isolated T Tauri systems. The 10 $μ$m silicate features show lower F$_{11.3}$/F$_{9.8}$ ratios at a given F$_{\mathrm{peak}}$, but current uncertainties prevent conclusions about their inner disk properties. Most disks display water emission from the inner disk, suggesting that even in these extreme environments rocky planets can form in the presence of water. The absence of strong line fluxes and other irradiation signatures suggests that the XUE disks have been truncated by external UV photons. However, this truncation does not appear to significantly impact the chemical richness of their inner regions. These findings indicate that even in extreme environments, IMTT disks can retain the ingredients necessary for rocky planet formation.

Individual stars in the Milky Way (MW) and its satellites have been shown to trace galaxy stellar mass dependent sequences in the $α$-abundance ([$α$/Fe]) vs metallicity ([Fe/H]) plane. Testing the universality of such sequences has been elusive as deep absorption-line spectra required for [$α$/Fe] and [Fe/H] measurements beyond the local group are mostly limited to integrated light from nearby, relatively high-mass, early-type galaxies. However, analogous to [$α$/Fe] vs [Fe/H] for stars, we now have log(O/Ar) vs 12+log(Ar/H) for the integrated nebular light of star-forming galaxies (SFGs). From Sloan-Digital Sky Survey (SDSS) observations of $\sim3000$ SFGs out to z$\sim0.3$, where we directly determined O & Ar abundances, we obtain for the first time the distribution of an ensemble of SFGs in the log(O/Ar) vs 12+log(Ar/H) plane. We show that higher (<M$\rm_{*}$>$\sim2.6\times10^9$M$_{\odot}$) and lower mass (<M$\rm_{*}$>$\sim1.7\times10^7$M$_{\odot}$) SFGs clearly trace distinct mass dependent sequences in this plane, qualitatively consistent with the mass dependence of chemical enrichment sequences observed for the stars in the MW and its satellites. Such sequences are consistent with expectations from galaxy chemical evolution (GCE) models that are driven primarily by the interplay of core-collapse and Type Ia supernovae.

We present a very deep CO(3-2) observation of a massive, gas-rich, main sequence, barred spiral galaxy at $z\approx1.52$. Our data were taken with the IRAM-NOEMA interferometer for a 12-antenna equivalent on-source integration time of $\sim$ 50 hours. We fit the major axis kinematics using forward modelling of a rotating disk, and then subtract the two-dimensional beam convolved best-fit model revealing signatures of planar non-circular motions in the residuals. The inferred in-plane radial velocities are remarkably large, of the order of $\approx60$ km/s. Direct comparisons with a high-resolution, simulated, gas-rich, barred galaxy, obtained with the moving mesh code AREPO and the TNG sub-grid model, show that the observed non-circular gas flows can be explained as radial flows driven by the central bar, with an inferred net inflow rate of the order of the SFR. Given the recent evidence for a higher-than-expected fraction of barred disk galaxies at cosmic noon, our results suggest that rapid gas inflows due to bars could be important evolutionary drivers for the dominant population of star-forming galaxies at the peak epoch of star and galaxy formation.

The liquid-vapor transition is a classic example of a discontinuous (first-order) phase transition. Such transitions underlie many phenomena in cosmology, nuclear and particle physics, and condensed-matter physics. They give rise to long-lived metastable states, whose decay can be driven by either thermal or quantum fluctuations. Yet, direct experimental observations of how these states collapse into a stable phase remain elusive in the quantum regime. Here, we use a trapped-ion quantum simulator to observe the real-time dynamics of ``bubble nucleation'' induced by quantum fluctuations. Bubbles are localized domains of the stable phase which spontaneously form, or nucleate, and expand as the system is driven across a discontinuous quantum phase transition. Implementing a mixed-field Ising spin model with tunable and time-dependent interactions, we track the microscopic evolution of the metastable state as the Hamiltonian parameters are varied in time with various speeds, bringing the system out of equilibrium. Site-resolved measurements reveal the emergence and evolution of finite-size quantum bubbles, providing direct insight into the mechanism by which the metastable phase decays. We also identify nonequilibrium scaling behavior near the transition, consistent with a generalized Kibble-Zurek mechanism. Our results demonstrate the power of quantum simulators to probe out-of-equilibrium many-body physics, including quantum bubble nucleation, a key feature of discontinuous quantum phase transitions, with application to studies of matter formation in the early universe.

Heavy chemical elements such as iron in the intra-cluster medium (ICM) of galaxy clusters are a signpost of the interaction between the gas and stellar components. Observations of the ICM metallicity in present-day massive systems, however, pose a challenge to the underlying assumption that the cluster galaxies have produced the amount of iron that enriches the ICM. We evaluate the iron share between ICM and stars within simulated galaxy clusters with the twofold aim of investigating the origin of possible differences with respect to observational findings and of shedding light on the observed excess of iron on the ICM with respect to expectations based on the observed stellar population. We evaluated the iron mass in gas and stars in a sample of 448 simulated systems with masses M500 > 1e14 Msun at z=0.07. These were extracted from the high-resolution (352 cMpc/h)^3 volume of the Magneticum cosmological hydrodynamical simulations. We compared our results with observational data of low-redshift galaxy clusters. The iron share in simulated clusters features a shallow dependence on the total mass, and its value is close to unity on average. In the most massive simulated systems, the iron share is thus smaller than observational values by almost an order of magnitude. The dominant contribution to this difference is related to the stellar component, whereas the chemical properties of the ICM agree well overall with the observations. We find larger stellar mass fractions in simulated massive clusters, which in turn yield higher stellar iron masses, than in observational data. Consistently with the modelling, we confirm that the stellar content within simulated present-day massive systems causes the metal enrichment in the ICM. It will be crucial to alleviate the stellar mass discrepancy between simulations and observations to definitely assess the iron budget in galaxy clusters.

We present analysis of the plateau and late-time phase properties of a sample of 39 Type II supernovae (SNe II) that show narrow, transient, high-ionization emission lines (i.e., "IIn-like") in their early-time spectra from interaction with confined, dense circumstellar material (CSM). Originally presented by Jacobson-Galán et al 2024a, this sample also includes multicolor light curves and spectra extending to late-time phases of 35 SNe with no evidence for IIn-like features at <2 days after first light. We measure photospheric phase light-curve properties for the distance-corrected sample and find that SNe II with IIn-like features have significantly higher luminosities and decline rates at +50 days than the comparison sample, which could be connected to inflated progenitor radii, lower ejecta mass, and/or persistent CSM interaction. However, we find no statistical evidence that the measured plateau durations and $^{56}$Ni masses of SNe II with and without IIn-like features arise from different distributions. We estimate progenitor zero-age main sequence (ZAMS) masses for all SNe with nebular spectroscopy through spectral model comparisons and find that most objects, both with and without IIn-like features, are consistent with progenitor masses <12.5 M$_{\odot}$. Combining progenitor ZAMS masses with CSM densities inferred from early-time spectra suggests multiple channels for enhanced mass loss in the final years before core collapse such as a convection-driven chromosphere or binary interaction. Finally, we find spectroscopic evidence for ongoing ejecta-CSM interaction at radii $>10^{16}$ cm, consistent with substantial progenitor mass-loss rates of $\sim 10^{-4}$--$10^{-5}$ M$_{\odot}$ yr$^{-1}$ ($v_w < 50$ km/s) in the final centuries to millennia before explosion.

Context. We present simulations of a massive young star cluster using \textsc{Nbody6++GPU} and \textsc{MOCCA}. The cluster is initially more compact than previously published models, with one million stars, a total mass of $5.86 \times 10^5~\mathrm{M}_{\odot}$, and a half-mass radius of $0.1~\mathrm{pc}$. Aims. We analyse the formation and growth of a very massive star (VMS) through successive stellar collisions and investigate the subsequent formation of an intermediate-mass black hole (IMBH) in the core of a dense star cluster. Methods. We use both direct \textit{N}-body and Monte Carlo simulations, incorporating updated stellar evolution prescriptions (SSE/BSE) tailored to massive stars and VMSs. These include revised treatments of stellar radii, rejuvenation, and mass loss during collisions. While the prescriptions represent reasonable extrapolations into the VMS regime, the internal structure and thermal state of VMSs formed through stellar collisions remain uncertain, and future work may require further refinement. Results. We find that runaway stellar collisions in the cluster core produce a VMS exceeding $5 \times 10^4~\mathrm{M}_{\odot}$ within 5 Myr, which subsequently collapses into an IMBH. Conclusions. Our model suggests that dense stellar environments may enable the formation of very massive stars and massive black hole seeds through runaway stellar collisions. These results provide a potential pathway for early black hole growth in star clusters and offer theoretical context for interpreting recent JWST observations of young, compact clusters at high redshift.

Virtually all extragalactic use cases of the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) require the use of galaxy redshift information, yet the vast majority of its sample of tens of billions of galaxies will lack high-fidelity spectroscopic measurements thereof, instead relying on photometric redshifts (photo-$z$) subject to systematic imprecision and inaccuracy best encapsulated by photo-$z$ probability density functions (PDFs). We present the version 1 release of Redshift Assessment Infrastructure Layers (RAIL), an open source Python library for at-scale probabilistic photo-$z$ estimation, initiated by the LSST Dark Energy Science Collaboration (DESC) with contributions from the LSST Interdisciplinary Network for Collaboration and Computing (LINCC) Frameworks team. RAIL's three subpackages provide modular tools for end-to-end stress-testing, including a forward modeling suite to generate realistically complex photometry, a unified API for estimating per-galaxy and ensemble redshift PDFs by an extensible set of algorithms, and built-in metrics of both photo-$z$ PDFs and point estimates. RAIL serves as a flexible toolkit enabling the derivation and optimization of photo-$z$ data products at scale for a variety of science goals and is not specific to LSST data. We thus describe to the extragalactic science community, including and beyond Rubin the design and functionality of the RAIL software library so that any researcher may have access to its wide array of photo-$z$ characterization and assessment tools.

Direct detection experiments have established the most stringent constraints on potential interactions between particle candidates for relic, thermal dark matter and Standard Model particles. To surpass current exclusion limits a new generation of experiments is being developed. The upcoming upgrade of the CRESST experiment will incorporate $\mathcal{O}$(100) detectors with different masses ranging from $\sim$2g to $\sim$24g, aiming to achieve unprecedented sensitivity to sub-GeV dark matter particles with a focus on spin-independent dark matter-nucleus scattering. This paper presents a comprehensive analysis of the planned upgrade, detailed experimental strategies, anticipated challenges, and projected sensitivities. Approaches to address and mitigate low-energy excess backgrounds $-$ a key limitation in previous and current sub-GeV dark matter searches $-$ are also discussed. In addition, a long-term roadmap for the next decade is outlined, including other potential scientific applications.

Effective Field Theory (EFT) modeling is expected to be a useful tool in the era of future higher-redshift galaxy surveys such as DESI-II and Spec-S5 due to its robust description of various large-scale structure tracers. However, large values of EFT bias parameters of higher-redshift galaxies could jeopardize the convergence of the perturbative expansion. In this paper we measure the bias parameters and other EFT coefficients from samples of two types of star-forming galaxies in the state-of-the-art MilleniumTNG and Astrid hydrodynamical simulations. Our measurements are based on the field-level EFT forward model that allows for precision EFT parameter measurements by virtue of cosmic variance cancellation. Specifically, we consider approximately representative samples of Lyman-break galaxies (LBGs) and Lyman-alpha emitters (LAEs) that are consistent with the observed (angular) clustering and number density of these galaxies at $z=3$. Reproducing the linear biases and number densities observed from existing LAE and LBG data, we find quadratic bias parameters that are roughly consistent with those predicted from the halo model coupled with a simple halo occupation distribution model. We also find non-perturbative velocity contributions (Fingers of God) of a similar size for LBGs to the familiar case of Luminous Red Galaxies. However, these contributions are quite small for LAEs despite their large satellite fraction values of up to $\sim 30\%$. Our results indicate that the effective momentum reach $k_{\rm{Max}}$ at $z=3$ for LAEs (LBGs) will be in the range $0.3-0.6 ~h\rm{Mpc}^{-1}$ ($0.2-0.8~h\rm{Mpc}^{-1}$), suggesting that EFT will perform well for high redshift galaxy clustering. This work provides the first step toward obtaining realistic simulation-based priors on EFT parameters for LAEs and LBGs.

We present one of the largest uniform optical spectroscopic surveys of X-ray selected sources to date that were observed as a pilot study for the Black Hole Mapper (BHM) survey. The BHM program of the Sloan Digital Sky Survey (SDSS)-V is designed to provide optical spectra for hundreds of thousands of X-ray selected sources from the SRG/eROSITA all-sky survey. This significantly improves our ability to classify and characterise the physical properties of large statistical populations of X-ray emitting objects. Our sample consists of 13079 sources in the eROSITA eFEDS performance verification field, 12011 of which provide reliable redshifts from 0<z<5.8. The vast majority of these objects were detected as point-like sources (X-ray flux limit F(0.5-2 keV)>6.5x10^-15 erg/s/cm^2) and were observed for about 20 years with fibre-fed SDSS spectrographs. After including all available redshift information for the eFEDS sources from the dedicated SDSS-V plate programme and archival data, we visually inspected the SDSS optical spectra to verify the reliability of these redshift measurements and the performance of the SDSS pipeline. The visual inspection allowed us to recover reliable redshifts (for 99% of the spectra with a signal-to-noise ratio of >2) and to assign classes to the sources, and we confirm that the vast majority of our sample consists of active galactic nuclei (AGNs). Only ~3% of the eFEDS/SDSS sources are Galactic objects. We also show the diversity of the optical spectra of the X-ray selected AGNs and provide spectral stacks with a high signal-to-noise ratio in various sub-samples with different redshift and optical broad-band colours. Our AGN sample contains optical spectra of (broad-line) quasars, narrow-line galaxies, and optically passive galaxies. It is considerably diverse in its colours and in its levels of nuclear obscuration.

Nuclear fission leads to the splitting of a nucleus into two fragments^1,2. Studying the distribution of the masses and charges of the fragments is essential for establishing the fission mechanisms and refining the theoretical models^3,4. It has value for our understanding of r-process nucleosynthesis^5,6, in which the fission of nuclei with extreme neutron-to-proton ratios is pivotal for determining astrophysical abundances and understanding the origin of the elements⁷ and for energy applications^8,9. Although the asymmetric distribution of fragments is well understood for actinides (elements in the periodic table with atomic numbers from 89 to 103) based on shell effects¹⁰, symmetric fission governs the scission process for lighter elements. However, unexpected asymmetric splits have been observed in neutron-deficient exotic nuclei¹¹, prompting extensive further investigations. Here we present measurements of the charge distributions of fission fragments for 100 exotic fissioning systems, 75 of which have never been measured, and establish a connection between the neutron-deficient sub-lead region and the well-understood actinide region. These new data comprehensively map the asymmetric fission island and provide clear evidence for the role played by the deformed Z = 36 proton shell of the light fragment in the fission of sub-lead nuclei. Our dataset will help constrain the fission models used to estimate the fission properties of nuclei with extreme neutron-to-proton ratios for which experimental data are unavailable.

Star clusters can interact and merge in galactic discs, haloes, or centres. We present direct N-body simulations of binary mergers of star clusters with <inline-formula><tex-math id="TM0001" notation="LaTeX">$M_{\star } = 2.7 \times 10^4 \, \mathrm{M_{\odot }}$</tex-math></inline-formula> each, using the N-body code BIFROST with subsystem regularization and post-Newtonian dynamics. We include 500 <inline-formula><tex-math id="TM0002" notation="LaTeX">$\mathrm{M_{\odot }}$</tex-math></inline-formula> massive black holes (MBHs) in the progenitors to investigate their impact on remnant evolution. The MBHs form hard binaries interacting with stars and stellar black holes (BHs). A few Myr after the cluster merger, this produces sizable populations of runaway stars (<inline-formula><tex-math id="TM0003" notation="LaTeX">$\sim$</tex-math></inline-formula>800 with <inline-formula><tex-math id="TM0004" notation="LaTeX">$v_{\mathrm{ej}} \gtrsim 50 \, \mathrm{km\, s^{-1}}$</tex-math></inline-formula>) and stellar BHs (<inline-formula><tex-math id="TM0005" notation="LaTeX">$\sim$</tex-math></inline-formula>30) escaping within 100 Myr. The remnants lose <inline-formula><tex-math id="TM0006" notation="LaTeX">$\sim 30{{\ \rm per\ cent}}$</tex-math></inline-formula> of their BH population and <inline-formula><tex-math id="TM0007" notation="LaTeX">$\sim 3{{\ \rm per\ cent}}$</tex-math></inline-formula> of their stars, with <inline-formula><tex-math id="TM0008" notation="LaTeX">$\sim$</tex-math></inline-formula>30 stars accelerated to high velocities <inline-formula><tex-math id="TM0009" notation="LaTeX">$\gtrsim 300 \, \mathrm{km\, s^{-1}}$</tex-math></inline-formula>. Comparison simulations of isolated clusters with central hard MBH binaries and cluster mergers without MBHs show that the process is driven by MBH binaries, while those with a single 1000 <inline-formula><tex-math id="TM0010" notation="LaTeX">$\mathrm{M_{\odot }}$</tex-math></inline-formula> MBH in isolated or merging clusters produce fewer runaway stars at lower velocities. Low-eccentricity merger orbits yield rotating remnants (<inline-formula><tex-math id="TM0011" notation="LaTeX">$v_{\mathrm{rot}} \sim 3 \, \mathrm{km\, s^{-1}}$</tex-math></inline-formula>), but probing the presence of MBHs via kinematics alone remains challenging. We expect the binary MBHs to merge within a Hubble time, producing observable gravitational-wave (GW) events detectable by future GW detectors such as the Einstein Telescope and Laser Interferometer Space Antenna. The results suggest that interactions with low-mass MBH binaries formed in merging star clusters are an important additional channel for producing runaway and high-velocity stars, free-floating stellar BHs, and compact objects.

In some scenarios, the dark matter relic abundance is set by the semi-annihilation of two dark matter particles into one dark matter particle and one Standard Model particle. These semi-annihilations might still be occurring today in the galactic center at a significant rate, generating a flux of boosted dark matter particles. We investigate the possible signals of this flux component in direct detection and neutrino experiments for sub-GeV dark matter masses. We show that for typical values of the semi-annihilation cross-section, the sensitivity of current experiments to the spin-independent dark matter-proton scattering cross-section can be several orders of magnitude larger than current constraints from cosmic-ray boosted dark matter. We also argue that the upcoming DARWIN and DUNE experiments may probe scattering cross-sections as low as <mml:math><mml:msup><mml:mrow><mml:mn>10</mml:mn></mml:mrow><mml:mrow><mml:mo>‑</mml:mo><mml:mn>37</mml:mn></mml:mrow></mml:msup><mml:mspace></mml:mspace><mml:msup><mml:mrow><mml:mi>cm</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:math> for masses between 30 MeV and 1 GeV.

Turbulent flows are known to produce enhanced effective magnetic and passive scalar diffusivities, which can fairly accurately be determined with numerical methods. It is now known that, if the flow is also helical, the effective magnetic diffusivity is reduced relative to the nonhelical value. Neither the usual second-order correlation approximation nor the various τ approaches have been able to capture this. Here we show that the helicity effect on the turbulent passive scalar diffusivity works in the opposite sense and leads to an enhancement. We have also demonstrated that the correlation time of the turbulent velocity field increases with the kinetic helicity. This is a key point in the theoretical interpretation of the obtained numerical results. Simulations in which helicity is being produced self-consistently by stratified rotating turbulence resulted in a turbulent passive scalar diffusivity that was found to be decreasing with increasing rotation rate.

We investigated the ability of the Euclid telescope to detect galaxy-scale gravitational lenses. To do so, we performed a systematic visual inspection of the 0.7 deg² Euclid Early Release Observations data towards the Perseus cluster using both the high-resolution I_E band and the lower-resolution Y_E , J_E, and H_E bands. Each extended source brighter than magnitude 23 in I_E was inspected by 41 expert human classifiers. This amounts to 12086 stamps of 10″ × 10″. We found 3 grade A and 13 grade B candidates. We assessed the validity of these 16 candidates by modelling them and checking that they are consistent with a single source lensed by a plausible mass distribution. Five of the candidates pass this check, five others are rejected by the modelling, and six are inconclusive. Extrapolating from the five successfully modelled candidates, we infer that the full 14 000 deg² of the Euclid Wide Survey should contain 100 000_{‑30 000}^{+ 70 000} galaxy-galaxy lenses that are both discoverable through visual inspection and have valid lens models. This is consistent with theoretical forecasts of 170 000 discoverable galaxy-galaxy lenses in Euclid. Our five modelled lenses have Einstein radii in the range 0'.'68 < θ_E < 1″.24, but their Einstein radius distribution is on the higher side when compared to theoretical forecasts. This suggests that our methodology is likely missing small-Einstein-radius systems. Whilst it is implausible to visually inspect the full Euclid dataset, our results corroborate the promise that Euclid will ultimately deliver a sample of around 10⁵ galaxy-scale lenses. ★This paper is published on behalf of the Euclid Consortium.

We present the ALMA-CRISTAL survey, an ALMA Cycle 8 Large Program designed to investigate the physical properties of star-forming galaxies at $4 \lesssim z \lesssim 6$ through spatially resolved, multi-wavelength observations. This survey targets 19 star-forming main-sequence galaxies selected from the ALPINE survey, using ALMA Band 7 observations to study [CII] 158 $μ$m line emission and dust continuum, complemented by JWST/NIRCam and HST imaging to map stellar and UV emission. The CRISTAL sample expanded to 39 after including newly detected galaxies in the CRISTAL fields, archival data, and pilot study targets. The resulting dataset provides a detailed view of gas, dust, and stellar structures on kiloparsec scales at the end of the era of reionization. The survey reveals diverse morphologies and kinematics, including rotating disks, merging systems, [CII] emission tails from potential interactions, and clumpy star formation. Notably, the [CII] emission in many cases extends beyond the stellar light seen in HST and JWST imaging. Scientific highlights include CRISTAL-10, exhibiting an extreme [CII] deficit similar to Arp 220; and CRISTAL-13, where feedback from young star-forming clumps likely causes an offset between the stellar clumps and the peaks of [CII] emission. CRISTAL galaxies exhibit global [CII]/FIR ratios that decrease with increasing FIR luminosity, similar to trends seen in local galaxies but shifted to higher luminosities, likely due to their higher molecular gas content. CRISTAL galaxies also span a previously unexplored range of global FIR surface brightness at high-redshift, showing that high-redshift galaxies can have elevated [CII]/FIR ratios. These elevated ratios are likely influenced by factors such as lower metallicity gas, the presence of significant extraplanar gas, and contributions from shock-excited gas.

Mass transfer in binary systems is the key process in the formation of various classes of objects, including merging binary black holes (BBHs) and neutron stars. Orbital evolution during mass transfer depends on how much mass is accreted and how much angular momentum is lost $-$ two of the main uncertainties in binary evolution. Here, we demonstrate that there is a fundamental limit to how close binary systems can get via stable mass transfer (SMT), that is robust against uncertainties in orbital evolution. Based on detailed evolutionary models of interacting systems with a BH accretor and a massive star companion, we show that the post-interaction orbit is always wider than $\sim10R_{\odot}$, even when extreme shrinkage due to L2 outflows is assumed. Systems evolving towards tighter orbits become dynamically unstable and result in stellar mergers. This separation limit has direct implications for the properties of BBH mergers: long delay times ($\gtrsim1 \rm Gyr$), and no high BH spins from the tidal spin-up of helium stars. At high metallicity, the SMT channel may be severely quenched due to Wolf-Rayet winds. The reason for the separation limit lies in the stellar structure, not in binary physics. If the orbit gets too narrow during mass transfer, a dynamical instability is triggered by a rapid expansion of the remaining donor envelope due to its near-flat entropy profile. The closest separations can be achieved from core-He burning ($\sim8-15R_{\odot}$) and Main Sequence donors ($\sim15-30R_{\odot}$), while Hertzsprung Gap donors lead to wider orbits ($\gtrsim30-50R_{\odot}$) and non-merging BBHs. These outcomes and mass transfer stability are determined by the entropy structures, which are governed by internal composition profiles. The formation of compact binaries is thus sensitive to chemical mixing in stars and may relate to the blue-supergiant problem.

RNA and proteins are the foundation of life and a natural starting point to explore its origins. However, the prebiotic relationship between the two is asymmetric. While RNA evolved to assemble proteins from amino acids, a significant mirror-symmetric effect of amino acids to trigger the synthesis of RNA was missing. We describe ambient alkaline conditions where amino acids, without additional chemical activators, promote RNA copolymerisation more than 100-fold, starting from prebiotically plausible ribonucleoside-2′,3′-cyclic phosphates. The observed effect is explained by acid-base catalysis, with optimal efficiency at pH values near the amine pKaH. The fold-change in oligomerisation yield is nucleobase-selective, resulting in increased compositional diversity necessary for subsequent molecular evolution and favouring the formation of natural 3′−5′ linkages. The elevated pH offers recycling of oligonucleotide sequences back to 2′,3′-cyclic phosphates, providing conditions for high-fidelity replication by templated ligation. The findings reveal a clear functional role of amino acids in the evolution of RNA earlier than previously assumed.

We present a novel way of probing non-gravitational dark matter interactions: dark astronomy, which leverages the dark radiation emitted by dissipative dark sectors. If the mediator of the dark matter self interactions is a dark photon with a small mass that kinetically mixes with the visible photon, the dark radiation flux becomes accessible to underground experiments. We argue that the emission may be dominantly longitudinally polarized, thereby enhancing the sensitivity of direct detection experiments such as XENON and SENSEI to this signal. We introduce a new detection mechanism based on resonant dark-photon-to-photon conversion at the surface of conducting materials, which offers unique directional sensitivity to dark radiation. This mechanism facilitates the development of experiments that combine dark matter detection techniques with methods of traditional astronomy, opening the possibility to map dark radiation sources within our galaxy.

We revisit the theory of background fields constructed on the BRST-algebra of a spinning particle with <inline-formula><mml:math><mml:mi>N</mml:mi></mml:math></inline-formula> = 4 worldline supersymmetry, whose spectrum contains the graviton but no other fields. On a generic background, the closure of the BRST-algebra implies the vacuum Einstein equations with a cosmological constant that is undetermined. On the other hand, in the "vacuum" background with no metric, the cohomology is given by a collection of free scalar- and vector fields. However, we show that only certain combinations of linear excitations, involving a vector field, can be extended beyond the linear level in turn, an Einstein metric in space-time.

A generalization of Wilson line operators at subleading power in the soft expansion has been recently introduced as an efficient building block of gravitational scattering amplitudes for non-spinning objects. The classical limit in this picture corresponds to the strict Regge limit, where the Post-Minkowskian (PM) expansion corresponds to the soft expansion, interpreted as a sum over correlations of soft emissions. Building on the well-studied worldline model with <inline-formula><tex-math>$$\mathcal{N}$$</tex-math></inline-formula> = 1 supersymmetry, in this work we extend the generalized Wilson line (GWL) approach to the case of spinning gravitating bodies. Specifically, at the quantum level we derive from first-principles a representation for the spin 1/2 GWL that is relevant for the all-order factorization of next-to-soft gravitons with fermionic matter, thus generalizing the exponentiation of single-emission next-to-soft theorems. At the classical level, we identify the suitable generalization of Wilson line operators that enables the generation of classical spin observables at linear order in spin. Thanks to the crucial role played by the soft expansion, the map from Grassmann variables to classical spin is manifest. We also comment on the relation between the GWL approach and the Worldline Quantum Field Theory as well as the Heavy Mass Effective Theory formalism. We validate the approach by rederiving known results in the conservative sector at 2PM order.