Context. The Kepler high-precision planetary sample has revealed a `radius valley' separating compact super-Earths from sub-Neptunes with lower densities. Super-Earths are generally assumed to be rocky planets that were probably born in situ, while the composition and formation of sub-Neptunes remains debated. Numerous statistical studies have explored planetary and stellar properties and their correlations to provide observational clues. However, no conclusive result on the origin of the radius valley or the composition of sub-Neptunes has been derived to date.
Aims: To provide more constraints, our aim is to investigate the distributions of the orbital spacing of sub-Neptunes and super-Earth planets in Kepler systems and compare their distributions with theoretical predictions of planet pairs of different formation pathways and compositions in synthetic planetary systems.
Methods: Based on the Kepler planetary sample, we derived the distributions of period ratios of sub-Neptune and super-Earth planet pairs. Using synthetic planetary systems generated by the Generation III Bern Model, we also obtained theoretical predictions of period ratio distributions of planet pairs of different compositions and origins.
Results: We find that Kepler sub-Neptune pairs show a significant preference to be near first-order mean motion resonances by a factor of 1.7<sub>−0.3</sub><sup>+0.3</sup>. This is smaller than the model predictions for `water-rich' pairs but larger than that of `water-poor' pairs by confidence levels of ~2σ. Kepler super-Earth pairs show no significant preference for mean motion resonances from a random distribution. The derived normalised fraction of near first-order resonances of actual Kepler super-Earth pairs is consistent with that of synthetic water-poor planet pairs but significantly (≳3σ) smaller than that of synthetic water-rich planet pairs.
Conclusions: The orbital migration has been more important for sub-Neptunes than for super-Earths, suggesting a partial ex situ formation of the former and an origin of the radius valley caused in part by distinct formation pathways. However, the model comparisons also show that sub-Neptunes in Kepler multiple systems are not likely (~2σ) to all be water-rich planets born ex situ but a mixture of the two (in situ and ex situ) pathways. Whereas, Kepler super-Earth planets are predominantly composed of water-poor planets that were born inside the ice line, likely through a series of giant impacts without large-scale migration.

Context. Gas in protoplanetary disks mostly cools via thermal accommodation with dust particles. Thermal relaxation is thus highly sensitive to the local dust size distributions and the spatial distribution of the grains. So far, the interplay between thermal relaxation and gas turbulence has not been dynamically modeled in hydrodynamic simulations of protoplanetary disks with dust.
Aims: We aim to study the effects of the vertical shear instability (VSI) on the thermal relaxation times, and vice versa. We are particularly interested in the influence of the initial dust grain size on the VSI and whether the emerging turbulence is sustained over long timescales.
Methods: We ran three axisymmetric hydrodynamic simulations of a protoplanetary disk including four dust fluids that initially resemble MRN size distributions of different initial grain sizes. From the local dust densities, we calculated the thermal accommodation timescale of dust and gas and used the result as the thermal relaxation time of the gas in our simulation. We included the effect of dust growth by applying the monodisperse dust growth rate and the typical growth limits.
Results: We find that the emergence of the VSI is strongly dependent on the initial dust grain size. Coagulation also counteracts the emergence of hydrodynamic turbulence in our simulations, as shown by others before. Starting a simulation with larger grains (100 μm) generally leads to a less turbulent outcome. While the inner disk regions (within ∼70 au) develop turbulence in all three simulations, we find that the simulations with larger particles do not develop VSI in the outer disk.
Conclusions: Our simulations with dynamically calculated thermal accommodation times based on the drifting and settling dust distribution show that the VSI, once developed in a disk, can be sustained over long timescales, even if grain growth is occurring. The VSI corrugates the dust layer and even diffuses the smaller grains into the upper atmosphere, where they can cool the gas. Whether the instability can emerge for a specific stratification depends on the initial dust grain sizes and the initial dust scale height. If the grains are initially ≳100 μm and if the level of turbulence is initially assumed to be low, we find no VSI turbulence in the outer disk regions.

We describe an attempt of string theoretic derivation of the Gibbons-Hawking entropy. Despite not admitting a de Sitter vacuum, the string theory, by the power of open-close correspondence, captures the Gibbons-Hawking entropy as the entropy of Chan-Paton species on a de Sitter-like state obtained via $D$-branes. Moreover, this derivation sheds a new light at the origin of the area-form, since the equality takes place for a critical 't Hooft coupling for which the species entropy of open strings saturates the area-law unitarity bound.

Observations and high-resolution hydrodynamical simulations indicate that massive star clusters assemble hierarchically from subclusters with a universal power-law cluster mass function. We study the consequences of such assembly for the formation of intermediate-mass black holes (IMBHs) at low metallicities (Z = 0.01 Z_⊙) with our updated N-body code BIFROST based on the hierarchical fourth-order forward integrator. BIFROST integrates few-body systems using secular and regularized techniques including post-Newtonian equations of motion up to order PN3.5 and gravitational-wave recoil kicks for BHs. Single stellar evolution is treated using the fast population synthesis code SEVN. We evolve three cluster assembly regions with N_tot = 1.70-2.35 × 10⁶ stars following a realistic IMF in ~1000 subclusters for t = 50 Myr. IMBHs with masses up to m_• ~ 2200 M_⊙ form rapidly mainly via the collapse of very massive stars (VMSs) assembled through repeated collisions of massive stars followed by growth through tidal disruption events and BH mergers. No IMBHs originate from the stars in the initially most massive clusters. We explain this by suppression of hard massive star binary formation at high velocity dispersions and the competition between core collapse and massive star lifetimes. Later the IMBHs form subsystems resulting in gravitational-wave BH-BH, IMBH-BH, and IMBH-IMBH mergers with a m_• ~ 1000 M_⊙ gravitational-wave detection being the observable prediction. Our simulations indicate that the hierarchical formation of massive star clusters in metal poor environments naturally results in formation of potential seeds for supermassive black holes.

We present a machine-learning-based model for the total density profiles of subhaloes with masses $M \gtrsim 7\times 10^8\, h^{-1}{\rm M}_\odot$ in the IllustrisTNG100 simulation. The model is based on an interpretable variational encoder (IVE) which returns the independent factors of variation in the density profiles within a low-dimensional representation, as well as the predictions for the density profiles themselves. The IVE returns accurate and unbiased predictions on all radial ranges, including the outer region profile where the subhaloes experience tidal stripping; here its fit accuracy exceeds that of the commonly used Einasto profile. The IVE discovers three independent degrees of freedom in the profiles, which can be interpreted in terms of the formation history of the subhaloes. In addition to the two parameters controlling the normalization and inner shape of the profile, the IVE discovers a third parameter that accounts for the impact of tidal stripping on to the subhalo outer profile; this parameter is sensitive to the mass loss experienced by the subhalo after its infall on to its parent halo. Baryonic physics in the IllustrisTNG galaxy formation model does not impact the number of degrees of freedom identified in the profile compared to the pure dark matter expectations, nor their physical interpretation. Our newly proposed profile fit can be used in strong lensing analyses or other observational studies which aim to constrain cosmology from small-scale structures.

Lambda- Cold Dark Matter (LambdaCDM) has been successful at explaining the large-scale structures in the universe but faces severe issues on smaller scales when compared to observations. Introducing self-interactions between dark matter particles claims to provide a solution to the small-scale issues in the LambdaCDM simulations while being consistent with the observations at large scales. The existence of the energy region in which these self-interactions between dark matter particles come close to saturating the S-wave unitarity bound can result in the formation of dark matter bound states called darkonium. In this scenario, all the low energy scattering properties are determined by a single parameter, the inverse scattering length gamma. In this work, we set bounds on gamma by studying the impact of darkonium on the observations at direct detection experiments using data from CRESST-III and XENON1T. The exclusion limits on gamma are then subsequently converted to exclusion limits on the self-interaction cross-section and compared with the constraints from astrophysics and N-body simulations.

RR Lyrae stars are excellent tracers of stellar populations for old, metal-poor components in the Milky Way Galaxy and the Local Group. Their luminosities have a metallicity-dependence, but determining spectroscopic [Fe/H] metallicities for RR Lyrae stars, especially at distances outside the solar neighbourhood, is challenging. Using 40 RRLs with metallicities derived from both Fe(II) and Fe(I) abundances, we verify the calibration between the [Fe/H] of RR Lyrae stars from the Calcium triplet. Our calibration is applied to all RR Lyrae stars with Gaia RVS spectra in Gaia DR3 as well as to 80 stars in the inner Galaxy from the BRAVA-RR survey. The co-added Gaia RVS RR Lyrae spectra provide RR Lyrae metallicities with an uncertainty of 0.25~dex, which is a factor of two improvement over the Gaia photometric RR Lyrae metallicities. Within our Galactic bulge RR Lyrae star sample, we find a dominant fraction with low energies without a prominent rotating component. Due to the large fraction of such stars, we interpret these stars as belonging to the $in-situ$ metal-poor Galactic bulge component, although we can not rule out that a fraction of these belong to an ancient accretion event such as Kraken/Heracles.

Context. Weak gravitational lensing offers a powerful method to investigate the projected matter density distribution within galaxy clusters, granting crucial insights into the broader landscape of dark matter on cluster scales.
Aims: In this study, we make use of the large photometric galaxy cluster data set derived from the publicly available Third Data Release of the Kilo-Degree Survey, along with the associated shear signal. Our primary objective is to model the peculiar sharp transition in the cluster profile slope, that is what is commonly referred to as the splashback radius. The data set under scrutiny includes 6962 galaxy clusters, selected by AMICO (an optimised detection algorithm of galaxy clusters) on the KiDS-DR3 data, in the redshift range of 0.1 ≤ z ≤ 0.6, all observed at a signal-to-noise ratio greater than 3.5.
Methods: Employing a comprehensive Bayesian analysis, we model the stacked excess surface mass density distribution of the clusters. We adopt a model from recent results on numerical simulations that capture the dynamics of both orbiting and infalling materials, separated by the region where the density profile slope undergoes a pronounced deepening.
Results: We find that the adopted profile successfully characterizes the cluster masses, consistent with previous works, and models the deepening of the slope of the density profiles measured with weak-lensing data up to the outskirts. Moreover, we measure the splashback radius of galaxy clusters and show that its value is close to the radius within which the enclosed overdensity is 200 times the mean matter density of the Universe, while theoretical models predict a larger value consistent with a low accretion rate. This points to a potential bias of optically selected clusters preferentially characterized by a high density at small scales compared to a pure mass-selected cluster sample.

Liquid-liquid phase separation is a widespread organizing principle of live cells, including, for example, the spatiotemporal organization of bacterial chromatin. The biophysics of phase-separating systems is often studied in vitro to avoid the complexity of the live-cell environment and facilitate application of advanced biophysical methods. One attractive method for measuring, e.g., partition coefficients, in such systems is fluorescence correlation spectroscopy (FCS). FCS circumvents some of the limitations of the widespread confocal laser scanning microscopy image-based measurements. Here, we describe how to perform partition coefficient measurements in biological phase-separating systems. Our protocol details typical workflows for the preparation of in vitro reconstituted condensates, FCS data acquisition, and subsequent data analysis, including corrections of some common artifacts. Our recommendations should help avoid many pitfalls of partition coefficient determination in these challenging systems.

One of science’s greatest challenges is determining how life can spontaneously emerge from a mixture of molecules. A complicating factor is that life and its molecules are inherently unstable—RNA and proteins are prone to hydrolysis and denaturation. For the de novo synthesis of life or to better understand its emergence at its origin, selection mechanisms are needed for unstable molecules. Here we present a chemically fuelled dynamic combinatorial library to model RNA oligomerization and deoligomerization and shine new light on selection and purification mechanisms under kinetic control. In the experiments, oligomers can only be sustained by continuous production. Hybridization is a powerful tool for selecting unstable molecules, offering feedback on oligomerization and deoligomerization rates. Moreover, we find that templation can be used to purify libraries of oligomers. In addition, template-assisted formation of oligomers within coacervate-based protocells changes its compartment’s physical properties, such as their ability to fuse. Such reciprocal coupling between oligomer production and physical properties is a key step towards synthetic life.

We review observational, experimental and theoretical results related to Dark Matter.

. We study the excitation spectrum of light and strange mesons in diffractive scattering. We identify different hadron resonances through partial wave analysis, which inherently relies on analysis models. Besides statistical uncertainties, the model dependence of the analysis introduces dominant sys tematic uncertainties. We discuss several of their sources for the π⁻π⁻π⁺ and K⁰_S K^¯ final states and present methods to reduce them. We have developed a new approach exploiting a-priori knowledge of signal continuity over adjacent final-state-mass bins to stably fit a large pool of partial-waves to our data, al lowing a clean identification of very small signals in our large data sets. For two-body final states of scalar particles, such as K⁰_S K^¯, mathematical ambigui ties in the partial-wave decomposition lead to the same intensity distribution for different combinations of amplitude values. We will discuss these ambiguities and present solutions to resolve or at least reduce the number of possible so lutions. Resolving these issues will allow for a complementary analysis of the a_J-like resonance sector in these two final states.

Charged and quasi-neutral beams propagating through an unmagnetised plasma are subject to numerous collisionless instabilities on the small scale of the plasma skin depth. The electrostatic two-stream instability, driven by longitudinal and transverse wakefields, dominates for dilute beams. This leads to modulation of the beam along the propagation direction and, for wide beams, transverse filamentation. A three-dimensional spatiotemporal two-stream theory for warm beams with a finite extent is developed. Unlike the cold beam limit, diffusion due to a finite emittance gives rise to a dominant wavenumber, and a cut-off wavenumber above which filamentation is suppressed. Particle-in-cell simulations give excellent agreement with the theoretical model. This work provides deeper insights into the effect of diffusion on filamentation of finite beams, crucial for comprehending plasma-based accelerators in laboratory and cosmic settings.

Radio emission has been detected from tens of white dwarfs, in particular in accreting systems. Additionally, radio emission has been predicted as a possible outcome of a planetary system around a white dwarf. We searched for 3 GHz radio continuum emission in 846 000 candidate white dwarfs previously identified in Gaia using the Very Large Array Sky Survey (VLASS) Epoch 1 Quick Look Catalogue. We identified 13 candidate white dwarfs with a counterpart in VLASS within 2 arcsec. Five of those were found not to be white dwarfs in follow-up or archival spectroscopy, whereas seven others were found to be chance alignments with a background source in higher resolution optical or radio images. The remaining source, WDJ204259.71+152108.06, is found to be a white dwarf and M-dwarf binary with an orbital period of 4.1 d and long-term stochastic optical variability, as well as luminous radio and X-ray emission. For this binary, we find no direct evidence of a background contaminant, and a chance alignment probability of only ≈2 per cent. However, other evidence points to the possibility of an unfortunate chance alignment with a background radio and X-ray emitting quasar, including an unusually poor Gaia DR3 astrometric solution for this source. With at most one possible radio emitting white dwarf found, we conclude that strong (≳1-3 mJy) radio emission from white dwarfs in the 3 GHz band is virtually non-existent outside of interacting binaries.

Time-delay cosmography is a powerful technique to constrain cosmological parameters, particularly the Hubble constant ($H_{0}$). The TDCOSMO collaboration is performing an ongoing analysis of lensed quasars to constrain cosmology using this method. In this work, we obtain constraints from the lensed quasar~WGD$\,$2038$-$4008~using new time-delay measurements and previous mass models by TDCOSMO. This is the first TDCOSMO lens to incorporate multiple lens modeling codes and the full time-delay covariance matrix into the cosmological inference. The models are fixed before the time delay is measured, and the analysis is performed blinded with respect to the cosmological parameters to prevent unconscious experimenter bias. We obtain $D_{\Delta t} = 1.68^{+0.40}_{-0.38}$ Gpc using two families of mass models, a power-law describing the total mass distribution, and a composite model of baryons and dark matter, although the composite model is disfavored due to kinematics constraints. In a flat $\Lambda$CDM cosmology, we constrain the Hubble constant to be $H_{0} = 65^{+23}_{-14}\, \rm km\ s^{-1}\,Mpc^{-1}$. The dominant source of uncertainty comes from the time delays, due to the low variability of the quasar. Future long-term monitoring, especially in the era of the Vera C. Rubin Observatory's Legacy Survey of Space and Time, could catch stronger quasar variability and further reduce the uncertainties. This system will be incorporated into an upcoming hierarchical analysis of the entire TDCOSMO sample, and improved time delays and spatially-resolved stellar kinematics could strengthen the constraints from this system in the future.

Theoretical models of structure formation predict the presence of a hot gaseous atmosphere around galaxies. While this hot circum-galactic medium (CGM) has been observationally confirmed through UV absorption lines, the detection of its direct X-ray emission remains scarce. We investigate theoretical predictions of the intrinsic CGM X-ray surface brightness (SB) using simulated galaxies and connect them to their global properties such as gas temperature, hot gas fraction and stellar mass. We select a sample of galaxies from the ultra-high resolution ($48\ \rm{cMpc\, h^{-1}}$) cosmological volume of the Magneticum Pathfinder set of hydrodynamical cosmological simulations. We classify them as star-forming (SF) or quiescent (QU) based on their specific star-formation rate. For each galaxy we generate X-ray mock data using the X-ray photon simulator PHOX, from which we obtain SB profiles out to the virial radius for different X-ray emitting components, namely gas, active galactic nuclei and X-ray binaries (XRBs). We fit a $\beta$-profile to each galaxy and observe trends between its slope and global quantities of the simulated galaxy. We find marginal differences between the average total SB profile of the CGM in SF and QU galaxies, with the contribution from hot gas being the largest ($>50\%$) at radii $r>0.05\,R_{\rm{vir}}$. The contribution from X-ray binaries (XRBs) equals the gas contribution for small radii and is non-zero for large radii. The galaxy population shows positive correlations between global properties and normalization of the SB profile. The slope of fitted $\beta$-profiles correlates strongly with the total gas luminosity, which in turn shows strong connections to the current accretion rate of the central super-massive black hole (SMBH).

Context. Planets orbiting members of open or globular clusters offer a great opportunity to study exoplanet populations systematically, as stars within clusters provide a mostly homogeneous sample, at least in chemical composition and stellar age. However, even though there have been coordinated efforts to search for exoplanets in stellar clusters, only a small number of planets have been detected. One successful example is the seven-year radial velocity (RV) survey `Search for giant planets in M 67' of 88 stars in the open cluster M 67, which led to the discovery of five giant planets, including three close-in (P < 10 days) hot-Jupiters.
Aims: In this work, we continue and extend the observation of stars in M 67, with the aim being to search for additional planets.
Methods: We conducted spectroscopic observations with the Habitable Planet Finder (HPF), HARPS, HARPS-North, and SOPHIE spectrographs of 11 stars in M 67. Six of our targets showed a variation or long-term trends in their RV during the original survey, while the other five were not observed in the original sample, bringing the total number of stars to 93.
Results: An analysis of the RVs reveals one additional planet around the turn-off point star S1429 and provides solutions for the orbits of stellar companions around S2207 and YBP2018. S1429 b is a warm-Jupiter on a likely circular orbit with a period of [77.48_-0.19^+0.18] days and a minimum mass of M sin i = 1.80 ± 0.2 M_J. We update the hot-Jupiter occurrence rate in M 67 to include the five new stars, deriving [4.2_-2.3^+4.1 %] when considering all stars, and [5.4_-3.0^+5.1 %] if binary star systems are removed.

A growing number of supernovae (SNe) are now known to exhibit evidence for significant interaction with a dense, pre-existing, circumstellar medium (CSM). SNe Ibn comprise one such class that can be characterized by both rapidly evolving light curves and persistent narrow He I lines. The origin of such a dense CSM in these systems remains a pressing question, specifically concerning the progenitor system and mass-loss mechanism. In this paper, we present multiwavelength data of the Type Ibn SN 2020nxt, including HST/STIS ultraviolet spectra. We fit the data with recently updated CMFGEN models designed to handle configurations for SNe Ibn. The UV coverage yields strong constraints on the energetics and, when combined with the CMFGEN models, offer new insight on potential progenitor systems. We find the most successful model is a ≲4 M_⊙ helium star that lost its $\sim 1\, {\rm M}_\odot$ He-rich envelope in the years preceding core collapse. We also consider viable alternatives, such as a He white dwarf merger. Ultimately, we conclude at least some SNe Ibn do not arise from single, massive (>30 M_⊙) Wolf-Rayet-like stars.

Context. Disc winds and planet-disc interactions are two crucial mechanisms that define the structure, evolution, and dispersal of protoplanetary discs. While winds are capable of removing material from discs, eventually leading to their dispersal, massive planets can shape their disc by creating sub-structures such as gaps and spiral arms.
Aims: We studied the interplay between an X-ray photoevaporative disc wind and the sub-structures generated due to planet-disc interactions to determine how their mutual interactions affect the disc's and the planet's evolution.
Methods: We performed 3D hydrodynamic simulations of viscous discs (α = 6.9 × 10⁻⁴) that host a Jupiter-like planet and undergo X-ray photoevaporation. We traced the gas flows within the disc and wind and measured the rate of accretion onto the planet, as well as the gravitational torque that is acting on it.
Results: Our results show that the planetary gap removes the wind's pressure support, allowing wind material to fall back into the gap. This opens new pathways for material from the inner disc (and part of the outer disc) to be redistributed through the wind towards the gap. Consequently, the gap becomes shallower and the flow of mass across the gap in both directions is significantly increased, as is the planet's mass-accretion rate (by factors of ≈5 and ≈2, respectively). Moreover, the wind-driven redistribution results in a denser inner disc and a less dense outer disc, which, combined with the recycling of a significant portion of the inner wind, leads to longer lifetimes for the inner disc, contrary to the expectation in a planet-induced photoevaporation scenario that has been proposed in the past.

The SRG/eROSITA All-Sky Survey (eRASS) is expected to contain ~100 quasars that emitted their light when the universe was less than a billion years old, i.e. at z>5.6. By selection, these quasars populate the bright end of the AGN X-ray luminosity function and their count offers a powerful demographic diagnostic of the parent super-massive black hole population. Of the >~ 400 quasars that have been discovered at z>5.6 to date, less than 15 % have been X-ray detected. We present a pilot survey to uncover the elusive X-ray luminous end of the distant quasar population. We have designed a quasar selection pipeline based on optical, infrared and X-ray imaging data from DES DR2, VHS DR5, CatWISE2020 and the eRASS. The core selection method relies on SED template fitting. We performed optical follow-up spectroscopy with the Magellan/LDSS3 instrument for the redshift confirmation of a subset of candidates. We have further obtained a deeper X-ray image of one of our candidates with Chandra ACIS-S. We report the discovery of five new quasars in the redshift range 5.6 < z < 6.1. Two of these quasars are detected in eRASS and are by selection X-ray ultra-luminous. These quasars are also detected at radio frequencies. The first one is a broad absorption line quasar which shows significant X-ray dimming over 3.5 years, i.e. about 6 months in the quasar rest frame. The second radio-detected quasar is a jetted source with compact morphology. We show that a blazar configuration is likely for this source, making it the second most distant blazar known to date. With our pilot study, we demonstrate the power of eROSITA as a discovery machine for luminous quasars in the epoch of reionization. The X-ray emission of the two eROSITA detected quasars are likely to be driven by different high-energetic emission mechanisms a diversity which will be further explored in a future systematic full-hemisphere survey.

Context. Radio jets are present in a diverse sample of AGN. However, the mechanisms of jet powering are not fully understood, and it remains unclear to what extent they obey mass-invariant scaling relations similar to those found for the triggering and fuelling of X-ray-selected AGN.
Aims: We use the multi-wavelength data in the eFEDS field observed by eROSITA/Spectrum-Roentgen-Gamma (SRG) and LOFAR to study the incidence of X-ray and radio AGN as a function of several stellar mass (M_*)-normalised AGN power indicators.
Methods: From the LOFAR - eFEDS survey, we defined a new sample of radio AGN, with optical counterparts from Legacy Survey DR9, according to a radio-excess relative to their host star formation rate. We further divided the sample into compact and complex radio morphologies. In this work, we used the subset matching to the well-characterised, highly complete spectroscopic GAMA09 galaxies (0 < z < 0.4). We release this value-added LOFAR - eFEDS catalogue^*. We calculated the fraction of GAMA09 galaxies hosting radio, X-ray, and both radio and X-ray AGN as functions of the specific black hole kinetic (λ_Jet) and radiative (λ_Edd) power.
Results: Despite the soft-X-ray eROSITA-selected sample, the incidence of X-ray AGN as a function of λ_Edd shows the same mass-invariance and power law slope (−0.65) as that found in previous studies once corrected for completeness. Across the M_* range probed, the incidence of compact radio AGN as a function of λ_Jet is described by a power law with constant slope, showing that it is not only high mass galaxies hosting high power jets and vice versa. This slope is steeper than that of the X-ray incidence, which has a value of around −1.5. Furthermore, higher-mass galaxies are more likely to host radio AGN across the λ_Jet range, indicating some residual mass dependence of jet powering. Upon adding complex radio morphologies, including 34 FRIIs, three of which are giant radio galaxies, the incidence not only shows a larger mass dependence but also a jet power dependence, being clearly boosted at high λ_Jet values. Importantly, the latter effect cannot be explained by such radio AGN residing in more dense environments (or more massive dark matter haloes). The similarity in the incidence of quiescent and star-forming radio AGN reveals that radio AGN are not only found in "red and dead" galaxies. Overall, our incidence analysis reveals some fundamental statistical properties of radio AGN samples, but highlights open questions regarding the use of a single radio luminosity-jet power conversion. We explore how different mass and accretion rate dependencies of the incidence can explain the observed results for varying disk-jet coupling models.

The source catalogue is available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (ftp://130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/686/A43 or on the LOFAR Surveys DR website: https://lofar-surveys.org/efeds.html

The exploration of outflows in protobinary systems presents a challenging yet crucial endeavour, offering valuable insights into the dynamic interplay between protostars and their evolution. In this study, we examine the morphology and dynamics of jets and outflows within the IRAS 4A protobinary system. This analysis is based on ALMA observations of SiO(5-4), H₂CO(3_{0, 3}-2_{0, 3}), and HDCO(4_{1, 4}-3_{1, 3}) with a spatial resolution of ~150 au. Leveraging an astrochemical approach involving the use of diverse tracers beyond traditional ones has enabled the identification of novel features and a comprehensive understanding of the broader outflow dynamics. Our analysis reveals the presence of two jets in the redshifted emission, emanating from IRAS 4A1 and IRAS 4A2, respectively. Furthermore, we identify four distinct outflows in the region for the first time, with each protostar, 4A1 and 4A2, contributing to two of them. We characterize the morphology and orientation of each outflow, challenging previous suggestions of bends in their trajectories. The outflow cavities of IRAS 4A1 exhibit extensions of 10 and 13 arcsec with position angles (PA) of 0° and -12°, respectively, while those of IRAS 4A2 are more extended, spanning 18 and 25 arcsec with PAs of 29° and 26°. We propose that the misalignment of the cavities is due to a jet precession in each protostar, a notion supported by the observation that the more extended cavities of the same source exhibit lower velocities, indicating they may stem from older ejection events.

A bright (m _F150W,AB = 24 mag), z = 1.95 supernova (SN) candidate was discovered in JWST/NIRCam imaging acquired on 2023 November 17. The SN is quintuply imaged as a result of strong gravitational lensing by a foreground galaxy cluster, detected in three locations, and remarkably is the second lensed SN found in the same host galaxy. The previous lensed SN was called "Requiem," and therefore the new SN is named "Encore." This makes the MACS J0138.0‑2155 cluster the first known system to produce more than one multiply imaged SN. Moreover, both SN Requiem and SN Encore are Type Ia SNe (SNe Ia), making this the most distant case of a galaxy hosting two SNe Ia. Using parametric host fitting, we determine the probability of detecting two SNe Ia in this host galaxy over a ∼10 yr window to be ≈3%. These observations have the potential to yield a Hubble constant (H ₀) measurement with ∼10% precision, only the third lensed SN capable of such a result, using the three visible images of the SN. Both SN Requiem and SN Encore have a fourth image that is expected to appear within a few years of ∼2030, providing an unprecedented baseline for time-delay cosmography.

We estimate the efficiency of mitigating the lensing B-mode polarization, the so-called delensing, for the LiteBIRD experiment with multiple external data sets of lensing-mass tracers. The current best bound on the tensor-to-scalar ratio, r, is limited by lensing rather than Galactic foregrounds. Delensing will be a critical step to improve sensitivity to r as measurements of r become more and more limited by lensing. In this paper, we extend the analysis of the recent LiteBIRD forecast paper to include multiple mass tracers, i.e., the CMB lensing maps from LiteBIRD and CMB-S4-like experiment, cosmic infrared background, and galaxy number density from Euclid- and LSST-like survey. We find that multi-tracer delensing will further improve the constraint on r by about 20%. In LiteBIRD, the residual Galactic foregrounds also significantly contribute to uncertainties of the B-modes, and delensing becomes more important if the residual foregrounds are further reduced by an improved component separation method.

We present initial results from the Dark Energy Spectroscopic Instrument (DESI) complete calibration of the colour-redshift relation (DC3R2) secondary target survey. Our analysis uses 230 k galaxies that overlap with KiDS-VIKING ugriZYJHK_s photometry to calibrate the colour-redshift relation and to inform photometric redshift (photo-z) inference methods of future weak lensing surveys. Together with emission line galaxies (ELGs), luminous red galaxies (LRGs), and the Bright Galaxy Survey (BGS) that provide samples of complementary colour, the DC3R2 targets help DESI to span 56 per cent of the colour space visible to Euclid and LSST with high confidence spectroscopic redshifts. The effects of spectroscopic completeness and quality are explored, as well as systematic uncertainties introduced with the use of common Self-Organizing Maps trained on different photometry than the analysis sample. We further examine the dependence of redshift on magnitude at fixed colour, important for the use of bright galaxy spectra to calibrate redshifts in a fainter photometric galaxy sample. We find that noise in the KiDS-VIKING photometry introduces a dominant, apparent magnitude dependence of redshift at fixed colour, which indicates a need for carefully chosen deep drilling fields, and survey simulation to model this effect for future weak lensing surveys.

We study the possibility of using the LiteBIRD satellite B-mode survey to constrain models of inflation producing specific features in CMB angular power spectra. We explore a particular model example, i.e. spectator axion-SU(2) gauge field inflation. This model can source parity-violating gravitational waves from the amplification of gauge field fluctuations driven by a pseudoscalar "axionlike" field, rolling for a few e-folds during inflation. The sourced gravitational waves can exceed the vacuum contribution at reionization bump scales by about an order of magnitude and can be comparable to the vacuum contribution at recombination bump scales. We argue that a satellite mission with full sky coverage and access to the reionization bump scales is necessary to understand the origin of the primordial gravitational wave signal and distinguish among two production mechanisms: quantum vacuum fluctuations of spacetime and matter sources during inflation. We present the expected constraints on model parameters from LiteBIRD satellite simulations, which complement and expand previous studies in the literature. We find that LiteBIRD will be able to exclude with high significance standard single-field slow-roll models, such as the Starobinsky model, if the true model is the axion-SU(2) model with a feature at CMB scales. We further investigate the possibility of using the parity-violating signature of the model, such as the TB and EB angular power spectra, to disentangle it from the standard single-field slow-roll scenario. We find that most of the discriminating power of LiteBIRD will reside in BB angular power spectra rather than in TB and EB correlations.

We explore the capability of measuring lensing signals in LiteBIRD full-sky polarization maps. With a 30 arcmin beam width and an impressively low polarization noise of 2.16 μK-arcmin, LiteBIRD will be able to measure the full-sky polarization of the cosmic microwave background (CMB) very precisely. This unique sensitivity also enables the reconstruction of a nearly full-sky lensing map using only polarization data, even considering its limited capability to capture small-scale CMB anisotropies. In this paper, we investigate the ability to construct a full-sky lensing measurement in the presence of Galactic foregrounds, finding that several possible biases from Galactic foregrounds should be negligible after component separation by harmonic-space internal linear combination. We find that the signal-to-noise ratio of the lensing is approximately 40 using only polarization data measured over 80% of the sky. This achievement is comparable to Planck's recent lensing measurement with both temperature and polarization and represents a four-fold improvement over Planck's polarization-only lensing measurement. The LiteBIRD lensing map will complement the Planck lensing map and provide several opportunities for cross-correlation science, especially in the northern hemisphere.

In this work, we propose a novel approach to the homotopy transfer procedure starting from a set of homotopy data such that the first differential complex is a differential graded module over the second one. We show that the module structure may be used to induce an $A_\infty$-algebra on the second differential complex, constructed in a similar fashion to the homotopy transfer $A_\infty$-algebra. We prove that, under certain conditions, the $A_\infty$-algebras obtained with this procedure are quasi-isomorphic to the homotopy transfer one. On the other hand, when the side conditions do not hold, we find that there are cases where the existence of an $A_\infty$-quasi-isomorphism with the homotopy transfer $A_\infty$-algebra is obstructed. In other words, we obtain a new $A_\infty$-algebra on the second complex, inequivalent to the homotopy transfer one. Lastly, we prove that these $A_\infty$-algebras are not infinitesimal Hochschild deformations of the homotopy transfer $A_\infty$-algebra.

Simulating high-resolution detector responses is a computationally intensive process that has long been challenging in Particle Physics. Despite the ability of generative models to streamline it, full ultra-high-granularity detector simulation still proves to be difficult as it contains correlated and fine-grained information. To overcome these limitations, we propose Intra-Event Aware Generative Adversarial Network (IEA-GAN). IEA-GAN presents a Transformer-based Relational Reasoning Module that approximates an event in detector simulation, generating contextualized high-resolution full detector responses with a proper relational inductive bias. IEA-GAN also introduces a Self-Supervised intra-event aware loss and Uniformity loss, significantly enhancing sample fidelity and diversity. We demonstrate IEA-GAN's application in generating sensor-dependent images for the ultra-high-granularity Pixel Vertex Detector (PXD), with more than 7.5 M information channels at the Belle II Experiment. Applications of this work span from Foundation Models for high-granularity detector simulation, such as at the HL-LHC (High Luminosity LHC), to simulation-based inference and fine-grained density estimation.

Testing gravity and the concordance model of cosmology, $\Lambda$CDM, at large scales is a key goal of this decade's largest galaxy surveys. Here we present a comparative study of dark matter power spectrum predictions from different numerical codes in the context of three popular theories of gravity that induce scale-independent modifications to the linear growth of structure: nDGP, Cubic Galileon and K-mouflage. In particular, we compare the predictions from full $N$-body simulations, two $N$-body codes with approximate time integration schemes, a parametrised modified $N$-body implementation and the analytic halo model reaction approach. We find the modification to the $\Lambda$CDM spectrum is in $2\%$ agreement for $z\leq1$ and $k\leq 1~h/{\rm Mpc}$ over all gravitational models and codes, in accordance with many previous studies, indicating these modelling approaches are robust enough to be used in forthcoming survey analyses under appropriate scale cuts. We further make public the new code implementations presented, specifically the halo model reaction K-mouflage implementation and the relativistic Cubic Galileon implementation.

Axion quark nuggets (AQNs) are hypothetical objects with a mass greater than a few grams and sub-micrometer size, formed during the quark-hadron transition. Originating from the axion field, they offer a possible resolution of the similarity between visible and dark components of the Universe. These composite objects behave as cold dark matter, interacting with ordinary matter and resulting in pervasive electromagnetic radiation throughout the Universe. This work aims to predict the electromagnetic signature in large-scale structures from the AQN-baryon interaction, accounting for thermal and non-thermal radiations. We use Magneticum hydrodynamical simulations to describe the distribution and dynamics of gas and dark matter at cosmological scales. We calculate the electromagnetic signature from radio, starting at $\nu \sim$ 1 GHz, up to a few keV X-ray energies. We find that the AQNs signature is characterized by monopole and fluctuation signals. The amplitude of both signals strongly depends on the average AQN mass and the ionization level of the baryonic environment. We identify a most optimistic scenario with a signal often near the sensitivity limit of existing instruments, such as FIRAS and the South Pole Telescope for high-resolution. Fluctuations in the Extra-galactic Background Light caused by the AQN can be tested with space-based imagers Euclid and James Webb Space Telescope. We also identify a minimal configuration, still out of reach of existing instruments, but future experiments might be able to pose constraints on the AQN model. We conclude that this is a viable dark matter model, which does not violate the canons of cosmology, nor existing observations. The best chances for testing this model reside in 1) ultra-deep IR and optical surveys, 2) spectral distorsions of the CMB and 3) low-frequency (1 GHz < $\nu $ < 100 GHz) and high-resolution ($\ell > 10^4$) observations.

We use the reconstructed properties of the Amaterasu particle, the second-highest energy cosmic ray ever detected, to map out three-dimensional constraints on the location of its unknown source. We highlight possible astrophysical sources that are compatible with these regions and requirements. Among these, M82, a powerful starburst galaxy, stands out as a strong candidate due to its position and proximity. To derive our constraints, we use CRPropa 3 to model all relevant propagation effects, including deflections in the Galactic and extra-Galactic magnetic fields. We consider key input quantities such as source distance, position, energy, and the strength and coherence length of the extra-Galactic magnetic field as free parameters. We then infer constraints on these parameters by applying approximate Bayesian computation. We present our results, demonstrating the impact of different assumptions for the arrival mass of the Amaterasu particle and the systematic uncertainties on the energy scale.

Context. In this work, we test feedback propagation models on the test case of 2MASS 0918+2117 (2M0918), a z = 0.149 X-ray variable AGN that shows tentative evidence for nuclear ultra-fast outflows (UFOs) in a 2005 XMM-Newton observation. We also investigate whether UFOs can be related to the observed X-ray variability.
Aims: We observed 2M0918 with XMM-Newton and NuSTAR in 2020 to confirm the presence of and characterize the UFOs. We performed a kinematic analysis of the publicly available 2005 SDSS optical spectrum to reveal and measure the properties of galaxy-scale ionized outflows. Furthermore, we constructed 20-year-long light curves of observed flux, line-of-sight column density, and intrinsic accretion rate from the spectra of the first four SRG/eROSITA all-sky surveys and archival observations from Chandra and XMM-Newton.
Methods: We detect UFOs with v ∼ 0.16c and galaxy-scale ionized outflows with velocities of ∼700 km s⁻¹. We also find that the drastic X-ray variability (factors > 10) can be explained in terms of variable obscuration and variable intrinsic luminosity.
Results: Comparing the energetics of the two outflow phases, 2M0918 is consistent with momentum-driven wind propagation. 2M0918 expands the sample of AGN with both UFOs and ionized gas winds from 5 to 6 and brings the sample of AGN hosting multiscale outflows to 19, contributing to a clearer picture of feedback physics. From the variations in accretion rate, column density, and ionization level of the obscuring medium, we propose a scenario that connects obscurers, an accretion enhancement, and the emergence of UFOs.

Global magnetic fields of early-type stars are commonly characterised by the mean longitudinal magnetic field <B_z> and the mean field modulus , derived from the circular polarisation and intensity spectra, respectively. Observational studies often report a root mean square (rms) of <B_z> and an average value of . In this work, I used numerical simulations to establish statistical relationships between these cumulative magnetic observables and the polar strength, B_d, of a dipolar magnetic field. I show that in the limit of many measurements randomly distributed in rotational phase, <B_z>_rms = 0.179_−0.043^+0.031 B_d and _avg = 0.691_−0.023^+0.020 B_d. The same values can be recovered with only three measurements, provided that the observations are distributed uniformly in the rotational phase. These conversion factors are suitable for ensemble analyses of large stellar samples, where each target is covered by a small number of magnetic measurements.

Context. The rapidly improving quality and resolution of both low surface brightness observations and cosmological simulations of galaxies enable us to address the important question of how the formation history is imprinted in the outer unrelaxed regions of galaxies, and to inspect the correlations of these imprints with another tracer of galaxy formation, the internal kinematics.
Aims: Using the hydrodynamical cosmological simulation called Magneticum Pathfinder, we identified tidal tails, shells, streams, and satellite planes, and connected them to the amount of rotational support and the formation histories of the host galaxies. This presents the first combined statistical census considering all these four types of features in hydrodynamical cosmological simulations.
Methods: Tidal features were visually classified from a three-dimensional rendering of the simulated galaxies by several scientists independently. Only features that were identified by at least half of the participating individuals were considered to be existing features. The data on satellite planes and kinematic properties of the simulated galaxies were taken from previous work. The results were compared to observations, especially from the MATLAS survey.
Results: Generally, prominent features are much more common around elliptical than around disk galaxies. Shells are preferentially found around kinematically slowly rotating galaxies in both simulations and observations, while streams can be found around all types of galaxies, with a slightly higher probability to be present around less rotationally supported galaxies. Tails and satellite planes, however, appear independently of the internal kinematics of the central galaxy, indicating that they are formed through processes that have not (yet) affected the internal kinematics. Prolate rotators have the overall highest probability to exhibit tidal features, but the highest likelihood for a specific type of feature is found for galaxies with kinematically distinct cores (KDCs), nearly 20% of which exhibit streams.
Conclusions: As shells are formed through radial merger events while streams are remnants of circular merger infall, this suggests that the orbital angular momentum of the merger event plays a more crucial role in transforming the host galaxy than previously anticipated. The existence of a shell around a given slow rotator furthermore is a sign of a radial merger formation for this particular slow rotator because one-third of the galaxies with a shell were transformed into slow rotators by the merger event that also caused the shells. The appearance of a stream around a KDC is a direct indicator for the multiple merger formation pathway of that KDC as opposed to the major merger pathway.

Context. Recent observations of protostellar cores suggest that most of the material in the protostellar phase is accreted along streamers. Streamers in this context are defined as velocity coherent funnels of denser material potentially connecting the large-scale environment to the small scales of the forming accretion disk.
Aims: Using simulations that simultaneously resolve the driving of turbulence on the filament scale as well as the collapse of the core down to protostellar disk scales, we aim to understand the effect of the turbulent velocity field on the formation of overdensities in the accretion flow.
Methods: We performed a three-dimensional numerical study on a core collapse within a turbulent filament using the RAMSES code and analysed the properties of overdensities in the accretion flow.
Results: We find that overdensities are formed naturally by the initial turbulent velocity field inherited from the filament and subsequent gravitational collimation. This leads to streams that are not really filamentary but show a sheet-like morphology. Moreover, they have the same radial infall velocities as the low density material. As a main consequence of the turbulent initial condition, the mass accretion onto the disk does not follow the predictions for solid body rotation. Instead, most of the mass is funneled by the overdensities to intermediate disk radii.

We developed a grid of stellar rotation models for low-mass and solar-type classical T Tauri stars (CTTS; 0.3M _⊙ < M _* < 1.2M _⊙). These models incorporate the star–disk interaction and magnetospheric ejections to investigate the evolution of the stellar rotation rate as a function of the mass of the star M _*, the magnetic field (B _*), and stellar wind ( ). We compiled and determined stellar parameters for 208 CTTS, such as projected rotational velocity , mass accretion rate , stellar mass M _*, ages, and estimated rotational periods using Transiting Exoplanet Survey Satellite (TESS) data. We also estimated a representative value of the mass-loss rate for our sample using the [O I] λ 6300 spectral line. Our results confirm that measurements in CTTS agree with the rotation rates provided by our spin models in the accretion-powered stellar winds picture. In addition, we used the approximate Bayesian computation technique to explore the connection between the model parameters and the observational properties of CTTS. We find that the evolution of with age might be regulated by variations in (1) the intensity of B _* and (2) the fraction of the accretion flow ejected in magnetic winds, removing angular momentum from these systems. The youngest stars in our sample (∼1 Myr) show a median branching ratio and median B _* ∼ 2000 G, in contrast to ∼0.01 and 1000 G, respectively, for stars with ages ≳3 Myr.

We study two-loop corrections to the scattering amplitude of four massive leptons in quantum electrodynamics. These amplitudes involve previously unknown elliptic Feynman integrals, which we compute analytically using the differential equation method. In doing so, we uncover the details of the elliptic geometry underlying this scattering amplitude and show how to exploit its properties to obtain compact, easy-to-evaluate series expansions that describe the scattering of four massive leptons in QED in the kinematical regions relevant for Bhabha and Møller scattering processes.

In this work we demonstrate that significant gains in performance and data efficiency can be achieved in High Energy Physics (HEP) by moving beyond the standard paradigm of sequential optimization or reconstruction and analysis components. We conceptually connect HEP reconstruction and analysis to modern machine learning workflows such as pretraining, finetuning, domain adaptation and high-dimensional embedding spaces and quantify the gains in the example usecase of searches of heavy resonances decaying via an intermediate di-Higgs system to four b-jets. To our knowledge this is the first example of a low-level feature extraction network finetuned for a downstream HEP analysis objective.

Context. There are several symbiotic stars (e.g., BF Cyg, Z And, and FN Sgr) in which periodic signals of tens of minutes have been detected. These periods have been interpreted as the spin period of magnetic white dwarfs that accrete through a magnetic stream originating from a truncated accretion disc.
Aims: To shed light on the origin of magnetic symbiotic stars, we investigated the system FN Sgr in detail. We searched for a reasonable formation pathway to explain its stellar and binary parameters including the magnetic field of the accreting white dwarf.
Methods: We used the MESA code to carry out pre-CE and post-CE binary evolution and determined the outcome of CE evolution assuming the energy formalism. For the origin and evolution of the white dwarf magnetic field, we adopted the crystallization scenario.
Results: We found that FN Sgr can be explained as follows. First, a non-magnetic white dwarf is formed through CE evolution. Later, during post-CE evolution, the white dwarf starts to crystallize and a weak magnetic field is generated. After a few hundred million years, the magnetic field penetrates the white dwarf surface and becomes detectable. Meanwhile, its companion evolves and becomes an evolved red giant. Subsequently, the white dwarf accretes part of the angular momentum from the red giant stellar winds. As a result, the white dwarf spin period decreases and its magnetic field reaches super-equipartition, getting amplified due to a rotation- and crystallization-driven dynamo. The binary then evolves into a symbiotic star, with a magnetic white dwarf accreting from an evolved red giant through atmospheric Roche-lobe overflow.
Conclusions: We conclude that the rotation- and crystallization-driven dynamo scenario, or any age-dependent scenario, can explain the origin of magnetic symbiotic stars reasonably well. This adds another piece to the pile of evidence supporting this scenario. If our formation channel is correct, our findings suggest that white dwarfs in most symbiotic stars formed through CE evolution might be magnetic, provided that the red giant has spent ≳3 Gyr as a main-sequence star.

Efficient algorithms are being developed to search for strong gravitational lens systems owing to increasing large imaging surveys. Neural networks have been successfully used to discover galaxy-scale lens systems in imaging surveys such as the Kilo Degree Survey, Hyper-Suprime Cam (HSC) Survey and Dark Energy Survey over the last few years. Thus, it has become imperative to understand how some of these networks compare, their strengths and the role of the training datasets as most of the networks make use of supervised learning algorithms. In this work, we present the first-of-its-kind systematic comparison and benchmarking of networks from four teams that have analysed the HSC Survey data. Each team has designed their training samples and developed neural networks independently but coordinated apriori in reserving specific datasets strictly for test purposes. The test sample consists of mock lenses, real (candidate) lenses and real non-lenses gathered from various sources to benchmark and characterise the performance of each of the network. While each team's network performed much better on their own constructed test samples compared to those from others, all networks performed comparable on the test sample with real (candidate) lenses and non-lenses. We also investigate the impact of swapping the training samples amongst the teams while retaining the same network architecture. We find that this resulted in improved performance for some networks. These results have direct implications on measures to be taken for lens searches with upcoming imaging surveys such as the Rubin-Legacy Survey of Space and Time, Roman and Euclid.

We consider the strong coupling limit of lattice QCD with massless staggered quarks and study the resource requirements for quantum simulating the theory in its Hamiltonian formulation. The bosonic Hilbert space of the color-singlet degrees of freedom grows quickly with the number of quark flavors, making it a suitable testing ground for resource considerations across different platforms. In particular, in addition to the standard model of computation with qubits, we consider mapping the theory to qudits $(d>2)$ and qumodes, as used on trapped-ion systems and photonic devices, respectively.

Active and sterile neutrinos may acquire "dark moments" via one-loop diagrams with a massless dark photon and new light particles in the loop. Due to the kinetic mixing between the dark photon and the Standard Model photon, neutrinos would obtain effective electromagnetic moments. This mechanism allows for enhanced electromagnetic moments that can evade constraints on particles directly charged under electromagnetism. We show that in a wide region of parameter space, the model features testable predictions for the anapole and magnetic moment of active and sterile neutrinos with dark matter direct detection experiments sensitive to the solar neutrino flux.

The mass loss mechanism of red supergiant stars is not well understood, even though it has crucial consequences for their stellar evolution and the appearance of supernovae that occur upon core-collapse. We argue that outgoing shock waves launched near the photosphere can support a dense chromosphere between the star's surface and the dust formation radius at several stellar radii. We derive analytic expressions for the time-averaged density profile of the chromosphere, and we use these to estimate mass loss rates due to winds launched by radiation pressure at the dust formation radius. These mass loss rates are similar to recent observations, possibly explaining the upward kink in mass loss rates of luminous red supergiants. Our models predict that low-mass red supergiants lose less mass than commonly assumed, while high-mass red supergiants lose more. The chromospheric mass of our models is $$0.01 solar masses, most of which lies within a few stellar radii. This can help explain the early light curves and spectra of type-II P supernovae without requiring extreme pre-supernova mass loss. We discuss implications for stellar evolution, type II-P supernovae, SN 2023ixf, and Betelgeuse.

We compute the suppression of bottomonium in the quark-gluon plasma using the three-loop QCD static potential. The potential describes the spin-averaged bottomonium spectrum below threshold with a less than 1% error. Within potential nonrelativistic quantum chromodynamics and an open quantum systems framework, we compute the evolution of the bottomonium density matrix. The values of the quarkonium transport coefficients are obtained from lattice QCD measurements of the bottomonium in-medium width and thermal mass shift; we additionally include for the first time a vacuum contribution to the dispersive coefficient γ . Using the three-loop potential and the values of the heavy quarkonium transport coefficients, we find that the resulting bottomonium nuclear modification factor is consistent with experimental observations, while at the same time reproducing the lattice measurements of the in-medium width.

The static QCD force from the lattice can be used to extract Λ_{MS ¯}, which determines the running of the strong coupling. Usually, this is done with a numerical derivative of the static potential. However, this introduces additional systematic uncertainties; thus, we use another observable to measure the static force directly. This observable consists of a Wilson loop with a chromoelectric field insertion. We work in the pure SU(3) gauge theory. We use gradient flow to improve the signal-to-noise ratio and to address the field insertion. We extract Λ_MS^{¯ n_f=0} from the data by exploring different methods to perform the zero-flow-time limit. We obtain the value √{8 t₀ }Λ_MS^{¯ n_f=0}=0.62 9_-26⁺²² , where t₀ is a flow-time reference scale. We also obtain precise determinations of several scales: r₀/r₁, √{8 t₀ }/r₀, √{8 t₀ }/r₁, and we compare these to the literature. The gradient flow appears to be a promising method for calculations of Wilson loops with chromoelectric and chromomagnetic insertions in quenched and unquenched configurations.

Off-lightcone Wilson-line operators are constructed using local operators connected by time-like or space-like Wilson lines, which ensure gauge invariance. Off-lightcone Wilson-line operators have broad applications in various contexts. For instance, space-like Wilson-line operators play a crucial role in determining quasi-distribution functions (quasi-PDFs), while time-like Wilson-line operators are essential for understanding quarkonium decay and production within the potential non-relativistic QCD (pNRQCD) framework. In this work, we establish a systematic approach for calculating the matching from the gradient-flow scheme to the MS ¯ scheme in the limit of small flow time for off-lightcone Wilson-line operators. By employing the one-dimensional auxiliary-field formalism, we simplify the matching procedure, reducing it to the matching of local current operators. We provide one-loop level matching coefficients for these local current operators. For the case of hadronic matrix element related to the quark quasi-PDFs, we show at one-loop level that the finite flow time effect is very small as long as the flow radius is smaller than the physical distance z, which is usually satisfied in lattice gradient flow computations. Applications include lattice gradient flow computations of quark/gluon quasi-PDFs, gluonic correlators related to quarkonium decay and production in pNRQCD, and spin-dependent potentials in terms of chromoelectric and chromomagnetic field insertions into a Wilson loop.

We scrutinize the Sommerfeld enhancement in dark matter pair annihilation for p-wave and higher-ℓ partial waves. For the Yukawa potential these feature a super-resonant Breit-Wigner peak in their velocity-dependence close to Sommerfeld resonances as well as a universal scaling with velocity for all ℓ ≥ 1 that differs from the s-wave case. We provide a quantum mechanical explanation for these phenomena in terms of quasi-bound states sustained by the centrifugal barrier of the partial-wave potential, and give approximate WKB expressions capturing the main effects. The impact of quasi-bound states is exemplified for wino dark matter and models with light mediators, with a focus on indirect detection signals. We note that quasi-bound states can also explain similar peaks in the bound-state formation and self-scattering cross sections.

The reduction and distortion of quantum correlations in the presence of classical noise leads to varied levels of inefficiency in the availability of entanglement as a resource for quantum information processing protocols. While generically minimizing required entanglement for mixed quantum states remains challenging, a class of many-body Gaussian quantum states (denoted N IC ) is here identified that exhibits two-mode bipartite entanglement structure, resembling that of pure states, for which the logarithmic negativity entanglement measure remains invariant upon inclusion of the classical correlations and optimal entanglement resources can be clearly quantified. This subclass is found to be embedded within a broader class of many-body Gaussian states (denoted N -SOL) that retain two-mode entanglement structure for detection processes. These two entanglement classes are relevant in theoretical and experimental applications from the scalar field vacuum to the local axial motional modes of trapped ion chains. Utilizing the subspace that heralds inseparability in response to partial transposition, a minimum noise filtering process is designed to be necessary, sufficient, and computable for determining membership in these classes of entanglement structure. Application of this process to spacelike regions of the free scalar field vacuum is found to improve resource upper bounds, providing new understanding of the entanglement required for the quantum simulation of quantum fields as observed by arrays of local detectors.

We study the evolution of the two-pion correlation function parameters with collision energy in the context of relativistic heavy-ion collisions within the NICA energy range. To this end, we perform UrQMD simulations in the cascade mode to produce samples of pions from <inline-formula id="IEq1"><mml:math><mml:mrow><mml:mn>5</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn>10</mml:mn><mml:mn>6</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> Bi+Bi collisions for each of the studied energies. The effects of the quantum-statistical correlations are introduced using the correlation afterburner code CRAB. We fit the correlation function using Gaussian, exponential and symmetric Lévy shapes and show that for all collision energies the latter provides the best fit. We separate the sample into pions coming from primary processes and pions originating from the decay of long-lived resonances, and show that the source size for the latter is significantly larger than for the former. The source size for the secondaries, is similar but in general larger than the size for the whole pion sample. To further characterize the pion source, we also simulate the effects of a non-ideal detector introducing a momentum smearing parameter, representing the minimum pair momentum and thus a maximum source size that can be resolved. The values of the correlation function intercept parameter are therefore modified from the values they attain for the perfect detector case. Using the core-halo picture of the source, we show that the values of the intercept parameter are influenced by the presence of a significant fraction of core pions coming from the decay of long-lived but slow-moving resonances. These findings serve as a benchmark to compare with future Monte Carlo studies that consider an Equation of State and thus allow for a phase transition within the studied energy domain.

During the public Kaggle competition "IceCube – Neutrinos in Deep Ice", thousands of reconstruction algorithms were created and submitted, aiming to estimate the direction of neutrino events recorded by the IceCube detector. Here we describe in detail the three ultimate best, award-winning solutions. The data handling, architecture, and training process of each of these machine learning models is laid out, followed up by an in-depth comparison of the performance on the Kaggle datatset. We show that on cascade events in IceCube above 10 TeV, the best Kaggle solution is able to achieve an angular resolution of better than 5<inline-formula id="IEq1"><mml:math><mml:mmultiscripts><mml:mrow></mml:mrow><mml:mrow></mml:mrow><mml:mo>∘</mml:mo></mml:mmultiscripts></mml:math></inline-formula>, and for tracks correspondingly better than 0.5<inline-formula id="IEq2"><mml:math><mml:mmultiscripts><mml:mrow></mml:mrow><mml:mrow></mml:mrow><mml:mo>∘</mml:mo></mml:mmultiscripts></mml:math></inline-formula>. These results indicate that the Kaggle solutions perform at a level comparable to the current state-of-the-art in the field, and that they may even be able to outperform existing reconstruction resolutions for certain types of events.

Feynman integrals are solutions to linear partial differential equations with polynomial coefficients. Using a triangle integral with general exponents as a case in point, we compare D-module methods to dedicated methods developed for solving differential equations appearing in the context of Feynman integrals, and provide a dictionary of the relevant concepts. In particular, we implement an algorithm due to Saito, Sturmfels, and Takayama to derive canonical series solutions of regular holonomic D-ideals, and compare them to asymptotic series derived by the respective Fuchsian systems.

We present a method for obtaining unbiased signal estimates in the presence of a significant unknown background, eliminating the need for a parametric model for the background itself. Our approach is based on a minimal set of conditions for observation and background estimators, which are typically satisfied in practical scenarios. To showcase the effectiveness of our method, we apply it to simulated data from the planned dielectric axion haloscope MADMAX.

Tidal features provide signatures of recent mergers and offer a unique insight into the assembly history of galaxies. The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will enable an unprecedentedly large survey of tidal features around millions of galaxies. To decipher the contributions of mergers to galaxy evolution it will be necessary to compare the observed tidal features with theoretical predictions. Therefore, we use cosmological hydrodynamical simulations NEWHORIZON, EAGLE, ILLUSTRISTNG, and MAGNETICUM to produce LSST-like mock images of z ~ 0 galaxies (z ~ 0.2 for NEWHORIZON) with $M_{\scriptstyle \star ,\text{ 30 pkpc}}\ge 10^{9.5}$ M$_{\scriptstyle \odot }$. We perform a visual classification to identify tidal features and classify their morphology. We find broadly good agreement between the simulations regarding their overall tidal feature fractions: $f_{{\small NewHorizon}}=0.40\pm 0.06$, $f_{{\small EAGLE}}=0.37\pm 0.01$, $f_{{\small TNG}}=0.32\pm 0.01$, and $f_{{\small Magneticum}}=0.32\pm 0.01$, and their specific tidal feature fractions. Furthermore, we find excellent agreement regarding the trends of tidal feature fraction with stellar and halo mass. All simulations agree in predicting that the majority of central galaxies of groups and clusters exhibit at least one tidal feature, while the satellite members rarely show such features. This agreement suggests that gravity is the primary driver of the occurrence of visually identifiable tidal features in cosmological simulations, rather than subgrid physics or hydrodynamics. All predictions can be verified directly with LSST observations.

In this study of the 'Resolving supermAssive Black hole Binaries In galacTic hydrodynamical Simulations' (RABBITS) series, we investigate the orbital evolution of supermassive black holes (SMBHs) during galaxy mergers. We simulate both disc and elliptical galaxy mergers using the KETJU code, which can simultaneously follow galaxy (hydro-)dynamics and small-scale SMBH dynamics with post-Newtonian corrections. With our SMBH binary subgrid model, we show how active galactic nuclei (AGNs) feedback affects galaxy properties and SMBH coalescence. We find that simulations without AGN feedback exhibit excessive star formation, resulting in merger remnants that deviate from observed properties. Kinetic AGN feedback proves more effective than thermal AGN feedback in expelling gas from the centre and quenching star formation. The different central galaxy properties, which are a result of distinct AGN feedback models, lead to varying rates of SMBH orbital decay. In the dynamical friction phase, galaxies with higher star formation and higher SMBH masses possess denser centres, become more resistant to tidal stripping, experience greater dynamical friction, and consequently form SMBH binaries earlier. As AGN feedback reduces gas densities in the centres, dynamical friction by stars dominates over gas. In the SMBH hardening phase, compared to elliptical mergers, disc mergers exhibit higher central densities of newly formed stars, resulting in accelerated SMBH hardening and shorter merger time-scales (i.e. $\lesssim 500$ Myr versus $\gtrsim 1$ Gyr). Our findings highlight the importance of AGN feedback and its numerical implementation in understanding the SMBH coalescing process, a key focus for low-frequency gravitational wave observatories.

Context. It has been claimed for more than a decade that energies other than orbital and thermodynamic internal are required to explain post-common envelope (CE) binaries with sufficiently long orbital periods (≳1 d) hosting AFGK-type main-sequence stars (∼0.5 − 2.0 M_⊙) paired with oxygen-neon white dwarfs (≳1.1 M_⊙). This would imply a completely different energy budget during CE evolution for these post-CE binaries in comparison to the remaining systems hosting M dwarfs and/or less massive white dwarfs.
Aims: In this first in a series of papers related to long-period post-CE binaries, we investigated whether extra energy is required to explain the currently known post-CE binaries with sufficiently long orbital periods consisting of oxygen-neon white dwarfs with AFGK-type main-sequence star companions.
Methods: We carried out binary population simulations with the BSE code adopting empirically derived inter-correlated main-sequence binary distributions for the initial binary population and assuming that the only energy, in addition to orbital, that help to unbind the CE is thermal energy. We also searched for the formation pathways of the currently known systems from the zero-age main-sequence binary to their present-day observed properties.
Results: Unlike what has been claimed for a long time, we show that all such post-CE binaries can be explained by assuming inefficient CE evolution, which is consistent with results achieved for the remaining post-CE binaries. There is therefore no need for an extra energy source. We also found that for CE efficiency close to 100%, post-CE binaries hosting oxygen-neon white dwarfs with orbital periods as long as one thousand days can be explained. For all known systems we found formation pathways consisting of CE evolution triggered when a highly evolved (i.e. when the envelope mass is comparable to the core mass), thermally pulsing, asymptotic giant branch star fills its Roche lobe at an orbital period of several thousand days. Due to the sufficiently low envelope mass and sufficiently long orbital period, the resulting post-CE orbital period can easily be several tens of days.
Conclusions: We conclude that the known post-CE binaries with oxygen-neon white dwarfs and AFGK-type main-sequence stars can be explained without invoking any energy source other than orbital and thermal energy. Our results strengthen the idea that the most common formation pathway of the overall population of post-CE binaries hosting white dwarfs is through inefficient CE evolution.

We present the results of a search for Miras and long-period variables (LPVs) in M33 using griJHK_S archival observations from the Canada-France-Hawaii Telescope. We use multiband information and machine learning techniques to identify and characterize these variables. We recover ~1300 previously discovered Mira candidates and identify ~13 000 new Miras and LPVs. We detect for the first time a clear first-overtone pulsation sequence among Mira candidates in this galaxy. We use O-rich, fundamental-mode Miras in the LMC and M33 to derive a distance modulus for the latter of μ = 24.629 ± 0.046 mag.

The inverted pendulum is a mechanical system with a rapidly oscillating pivot point. Using techniques similar in spirit to the methodology of effective field theories, we derive an effective Lagrangian that allows for the systematic computation of corrections to the so-called Kapitza equation. The derivation of the effective potential of the system requires non-trivial matching conditions, which need to be determined order by order in the power-counting of the problem. The convergence behavior of the series is investigated on the basis of high-order results obtained by this method. The results from this analysis can be used to determine the regions of parameter space, in which the inverted position of the pendulum is stable or unstable to high precision.

Within the framework of non-relativistic QCD (NRQCD) effective field theory, we study the leptoproduction of $J/{\psi}$ at next-to-leading order in perturbative QCD for both unpolarized and polarized electron-ion collisions. We demonstrate that the $J/{\psi}$-tagged deep inelastic scattering in the future Electron-Ion Collider can be served as a golden channel for the reasons including constraining NRQCD long distance matrix elements, probing the nuclear gluon distribution functions, as well as investigating the gluon helicity distribution inside a longitudinal polarized proton.

We present forecasts for constraints on the Hu and Sawicki f (R ) modified gravity model using realistic mock data representative of future cluster and weak lensing surveys. We create mock thermal Sunyaev-Zel'dovich effect selected cluster samples for SPT-3G and CMB-S4 and the corresponding weak gravitational lensing data from next-generation weak-lensing (ngWL) surveys like Euclid and Rubin. We employ a state-of-the-art Bayesian likelihood approach that includes all observational effects and systematic uncertainties to obtain constraints on the f (R ) gravity parameter log₁₀|f_{R 0}|. In this analysis we vary the cosmological parameters [Ω_m,Ω_νh²,h ,A_s,n_s,log₁₀|f_{R 0}|], which allows us to account for possible degeneracies between cosmological parameters and f (R ) modified gravity. The analysis accounts for f (R ) gravity via its effect on the halo mass function which is enhanced on cluster mass scales compared to the expectations within general relativity (GR). Assuming a fiducial GR model, the upcoming cluster dataset SPT -3 G ×ngWL is expected to obtain an upper limit of log₁₀|f_{R 0}|<-5.95 at 95% credibility, which significantly improves upon the current best bounds. The CMB -S 4 ×ngWL dataset is expected to improve this even further to log₁₀|f_{R 0}|<-6.23 . Furthermore, f (R ) gravity models with log₁₀|f_{R 0}|≥-6 , which have larger numbers of clusters, would be distinguishable from GR with both datasets. We also report degeneracies between log₁₀|f_{R 0}| and Ω_m as well as σ₈ for log₁₀|f_{R 0}|>-6 and log₁₀|f_{R 0}|>-5 respectively. Our forecasts indicate that future cluster abundance studies of f (R ) gravity will enable substantially improved constraints that are competitive with other cosmological probes.

A relativistic microscopic optical model potential for nucleon-nucleus scattering is developed based on the ab initio relativistic Brueckner-Hartree-Fock (RBHF) theory with the improved local density approximation, which is abbreviated as the RBOM potential. Both real and imaginary parts of the single-particle potentials in symmetric and asymmetric nuclear matter at various densities are determined uniquely in the full Dirac space. The density distributions of the target nuclei are calculated by the covariant energy density functional theory with the density functional PC-PK1. The central and spin-orbit terms of the optical potentials are quantitatively consistent with the relativistic phenomenological optical potentials. The performance of the RBOM potential is evaluated by considering proton scattering with incident energy E ≤200 MeV on five target nuclei, ²⁰⁸Pb, ¹²⁰Sn, ⁹⁰Zr, ⁴⁸Ca, and ⁴⁰Ca. Scattering observables including the elastic scattering angular distributions, analyzing powers, spin rotation functions, and reaction cross sections are analyzed. Theoretical predictions show good agreements with the experimental data and the results derived from phenomenological optical potentials. We anticipate that the RBOM potential can provide reference for other phenomenological and microscopic optical model potentials, as well as reliable descriptions for nucleon scattering on exotic nuclei in the era of rare-isotope beams.

We report on the SRG/eROSITA detection of ultra-soft ($kT=47^{+5}_{-5}$ eV) X-ray emission (L_X =$2.5^{+0.6}_{-0.5} \times 10^{43}$ erg s^-1) from the tidal disruption event (TDE) candidate AT 2022dsb ~14 d before peak optical brightness. As the optical luminosity increases after the eROSITA detection, then the 0.2-2 keV observed flux decays, decreasing by a factor of ~39 over the 19 d after the initial X-ray detection. Multi-epoch optical spectroscopic follow-up observations reveal transient broad Balmer emission lines and a broad He II 4686 Å emission complex with respect to the pre-outburst spectrum. Despite the early drop in the observed X-ray flux, the He II 4686 Å complex is still detected for ~40 d after the optical peak, suggesting the persistence of an obscured hard ionizing source in the system. Three outflow signatures are also detected at early times: (i) blueshifted H α emission lines in a pre-peak optical spectrum, (ii) transient radio emission, and (iii) blueshifted Ly α absorption lines. The joint evolution of this early-time X-ray emission, the He II 4686 Å complex, and these outflow signatures suggests that the X-ray emitting disc (formed promptly in this TDE) is still present after optical peak, but may have been enshrouded by optically thick debris, leading to the X-ray faintness in the months after the disruption. If the observed early-time properties in this TDE are not unique to this system, then other TDEs may also be X-ray bright at early times and become X-ray faint upon being veiled by debris launched shortly after the onset of circularization.

The arrival of the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST), Euclid-Wide and Roman wide-area sensitive surveys will herald a new era in strong lens science in which the number of strong lenses known is expected to rise from $\mathcal {O}(10^3)$ to $\mathcal {O}(10^5)$. However, current lens-finding methods still require time-consuming follow-up visual inspection by strong lens experts to remove false positives which is only set to increase with these surveys. In this work, we demonstrate a range of methods to produce calibrated probabilities to help determine the veracity of any given lens candidate. To do this we use the classifications from citizen science and multiple neural networks for galaxies selected from the Hyper Suprime-Cam survey. Our methodology is not restricted to particular classifier types and could be applied to any strong lens classifier which produces quantitative scores. Using these calibrated probabilities, we generate an ensemble classifier, combining citizen science, and neural network lens finders. We find such an ensemble can provide improved classification over the individual classifiers. We find a false-positive rate of 10^-3 can be achieved with a completeness of 46 per cent, compared to 34 per cent for the best individual classifier. Given the large number of galaxy-galaxy strong lenses anticipated in LSST, such improvement would still produce significant numbers of false positives, in which case using calibrated probabilities will be essential for population analysis of large populations of lenses and to help prioritize candidates for follow-up.

Polarization of the cosmic microwave background (CMB) can help probe the fundamental physics behind cosmic inflation via the measurement of primordial B modes. As this requires exquisite control over instrumental systematics, some next-generation CMB experiments plan to use a rotating half-wave plate (HWP) as polarization modulator. However, the HWP non-idealities, if not properly treated in the analysis, can result in additional systematics. In this paper, we present a simple, semi-analytical end-to-end model to propagate the HWP non-idealities through the macro-steps that make up any CMB experiment (observation of multi-frequency maps, foreground cleaning, and power spectra estimation) and compute the HWP-induced bias on the estimated tensor-to-scalar ratio, r. We find that the effective polarization efficiency of the HWP suppresses the polarization signal, leading to an underestimation of r. Laboratory measurements of the properties of the HWP can be used to calibrate this effect, but we show how gain calibration of the CMB temperature can also be used to partially mitigate it. On the basis of our findings, we present a set of recommendations for the HWP design that can help maximize the benefits of gain calibration.

Strongly lensed systems with peculiar configurations allow us to probe the local properties of the deflecting lens mass while simultaneously testing general profile assumptions. The quasar HE0230−2130 is lensed by two galaxies at similar redshifts (Δz ∼ 0.003) into four observed images. Using modeled quasar positions from fitting the brightness of the quasar images in ground-based imaging data from the Magellan telescope, we find that lens-mass models where each of these two galaxies is parametrized with a singular power-law (PL) profile predict five quasar images. One of the predicted images is unobserved despite it being distinctively offset from the lensing galaxies and likely bright enough to be observable. This missing image gives rise to new opportunities to study the mass distribution of these galaxies. To interpret the quad configuration of the system, we tested 12 different profile assumptions with the aim of obtaining lens-mass models that correctly predict only four observed images. We tested the effects of adopting: cored profiles for the lensing galaxies; external shear; and additional profiles to represent a dark matter clump. We find that half of our model classes can produce the correct image multiplicity. By comparing the Bayesian evidence of different model parametrizations, we favor two model classes: (i) one that incorporates two singular PL profiles for the lensing galaxies and a cored isothermal sphere in the region of the previously predicted fifth image (rNIS profile), and (ii) one with a bigger lensing galaxy parametrized by a singular PL profile and the smaller galaxy by a cored PL profile with external shear. We estimated the mass of the rNIS clump for each candidate model of our final Markov chain Monte Carlo sample, and find that only 2% are in the range of 10⁶ M_⊙ ≤ M_rNIS ≤ 10⁹ M_⊙, which is the predicted mass range of dark matter subhalos in cold dark matter simulations, or the mass of dark-matter-dominated and low-surface-brightness galaxies. We therefore favor the models with a cored mass distribution for the lens galaxy close to the predicted fifth image. Our study further demonstrates that lensed quasar images are sensitive to the dark matter structure in the gravitational lens. We are able to describe this exotic lensing configuration with relatively simple models, which demonstrates the power of strong lensing for studying galaxies and lens substructure.

We report the identification of 64 new candidates of compact galaxies, potentially hosting faint quasars with bolometric luminosities of L_bol = 10⁴³-10⁴⁶ erg s⁻¹, residing in the reionization epoch within the redshift range of 6 ≲ z ≲ 8. These candidates were selected by harnessing the rich multiband datasets provided by the emerging JWST-driven extragalactic surveys, focusing on COSMOS-Web, as well as JADES, UNCOVER, CEERS, and PRIMER. Our search strategy includes two stages: applying stringent photometric cuts to catalog-level data and detailed spectral energy distribution fitting. These techniques effectively isolate the quasar candidates while mitigating contamination from low-redshift interlopers, such as brown dwarfs and nearby galaxies. The selected candidates indicate physical traits compatible with low-luminosity active galactic nuclei, likely hosting ≈10⁵-10⁷ M_⊙ supermassive black holes (SMBHs) living in galaxies with stellar masses of ≈10⁸-10¹⁰ M_⊙. The SMBHs selected in this study, on average, exhibit an elevated mass compared to their hosts, with the mass ratio distribution slightly higher than those of galaxies in the local Universe. As with other high-z studies, this is at least in part due to the selection method for these quasars. An extensive Monte Carlo analysis provides compelling evidence that heavy black hole seeds from the direct collapse scenario appear to be the preferred pathway to mature this specific subset of SMBHs by z ≈ 7. Notably, most of the selected candidates might have emerged from seeds with masses of ∼10⁵ M_⊙, assuming a thin disk accretion with an average Eddington ratio of f_Edd = 0.6 ± 0.3 and a radiative efficiency of ϵ = 0.2 ± 0.1. This work underscores the significance of further spectroscopic observations, as the quasar candidates presented here offer exceptional opportunities to delve into the nature of the earliest galaxies and SMBHs that formed during cosmic infancy.

FITS files and full Table B.1 are available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (ftp://130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/685/A25

Context. The features in the light curves and spectra of many Type I and Type II supernovae (SNe) can be understood by assuming an interaction of the SN ejecta with circumstellar matter (CSM) surrounding the progenitor star. This suggests that many massive stars may undergo various degrees of envelope stripping shortly before exploding, and may therefore produce a considerable diversity in their pre-explosion CSM properties.
Aims: We explore a generic set of about 100 detailed massive binary evolution models in order to characterize the amount of envelope stripping and the expected CSM configurations.
Methods: Our binary models were computed with the MESA stellar evolution code, considering an initial primary star mass of 12.6 M_⊙ and secondaries with initial masses of between ∼12 M_⊙ and ∼1.3 M_⊙, and focus on initial orbital periods above ∼500 d. We compute these models up to the time of iron core collapse in the primary.
Results: Our models exhibit varying degrees of stripping due to mass transfer, resulting in SN progenitor models ranging from fully stripped helium stars to stars that have not been stripped at all. We find that Roche lobe overflow often leads to incomplete stripping of the mass donor, resulting in a large variety of pre-SN envelope masses. In many of our models, the red supergiant (RSG) donor stars undergo core collapse during Roche lobe overflow, with mass transfer and therefore system mass-loss rates of up to 0.01 M_⊙ yr⁻¹ at that time. The corresponding CSM densities are similar to those inferred for Type IIn SNe, such as <ASTROBJ>SN 1998S</ASTROBJ>. In other cases, the mass transfer becomes unstable, leading to a common-envelope phase at such late time that the mass donor explodes before the common envelope is fully ejected or the system has merged. We argue that this may cause significant pre-SN variability, as witnessed for example in <ASTROBJ>SN 2020tlf</ASTROBJ>. Other models suggest a common-envelope ejection just centuries before core collapse, which may lead to the strongest interactions, as observed in superluminous Type IIn SNe, such as <ASTROBJ>SN 1994W</ASTROBJ> and <ASTROBJ>SN 2006gy</ASTROBJ>.
Conclusions: Wide massive binaries exhibit properties that may not only explain the diverse envelope stripping inferred in Type Ib, IIb, IIL, and IIP SNe, but also offer a natural framework to understand a broad range of hydrogen-rich interacting SNe. On the other hand, the flash features observed in many Type IIP SNe, such as <ASTROBJ>SN 2013fs</ASTROBJ>, may indicate that RSG atmospheres are more extended than currently assumed; this could enhance the parameter space for wide binary interaction.

Among the well-known methods to approximate derivatives of expectancies computed by Monte-Carlo simulations, averages of pathwise derivatives are often the easiest one to apply. Computing them via algorithmic differentiation typically does not require major manual analysis and rewriting of the code, even for very complex programs like simulations of particle-detector interactions in high-energy physics. However, the pathwise derivative estimator can be biased if there are discontinuities in the program, which may diminish its value for applications. This work integrates algorithmic differentiation into the electromagnetic shower simulation code HepEmShow based on G4HepEm, allowing us to study how well pathwise derivatives approximate derivatives of energy depositions in a sampling calorimeter with respect to parameters of the beam and geometry. We found that when multiple scattering is disabled in the simulation, means of pathwise derivatives converge quickly to their expected values, and these are close to the actual derivatives of the energy deposition. Additionally, we demonstrate the applicability of this novel gradient estimator for stochastic gradient-based optimization in a model example.

IRAS04368+2557 in L1527 is a Class 0/I protostar with a clear disk-envelope system revealed by previous Atacama Large Millimeter/submillimeter Array (ALMA) observations. In this paper, we discuss the flared structure of this source with observed sulfur-bearing molecules included in the FAUST ALMA large program. The analyses of molecular distributions and kinematics have shown that CS, SO, and OCS trace different regions of the disk-envelope system. To evaluate the temperature across the disk, we derive rotation temperature with the two observed SO lines. The temperature profile shows a clear, flared "butterfly" structure with the higher temperature being ∼50 K and the central lower temperature region (<30 K) coinciding with the continuum peak, suggesting dynamically originated heating rather than radiation heating from the central protostar. Other physical properties, including column densities, are also estimated and further used to demonstrate the vertical structure of the disk-envelope system. The "warped" disk structure of L1527 is confirmed with our analyses, showing that sulfur-bearing molecules are not only effective material probes but also sufficient for structural studies of protostellar systems.

In cosmological simulations of large-scale structure star formation and feedback in galaxies are modelled by so-called sub-grid models, that represent a physically motivated approximation of processes occurring below the resolution limit. However, when additional physical processes are considered in these simulations, for instance, magnetic fields or cosmic rays, they are often not consistently coupled within the descriptions of the underlying sub-grid star formation models. Here, we present a careful study on how one of the most commonly used sub-grid models for star formation in current large-scale cosmological simulations can be modified to self consistently include the effects of non-thermal components (e.g., magnetic fields) within the fluid. We demonstrate that our new modelling approach, that includes the magnetic pressure as an additional regulation on star formation, can reproduce global properties of the magnetic field within galaxies in a setup of an isolated Milky Way-like galaxy simulation, but is also successful in reproducing local properties such as the anti-correlation between the local magnetic field strength with the local star formation rate as observed in galaxies (i.e. NGC 1097). This reveals how crucial a consistent treatment of different physical processes is within cosmological simulations and gives guidance for future simulations.

Stellar winds of massive ($\gtrsim 9\, \mathrm{M_\odot }$) and very massive ($\gtrsim 100\, \mathrm{M_\odot }$) stars may play an important role in the metal-enrichment during the formation of star clusters. With novel high-resolution hydrodynamical GRIFFIN-project simulations, we investigate the rapid recycling of stellar wind-material during the formation of massive star clusters up to $M_\mathrm{cluster}\sim 2\times 10^5\, \mathrm{M_\odot }$ in a low-metallicity dwarf galaxy starburst. The simulation realizes new stars from a stellar initial mass function (IMF) between $0.08$ and $\sim 400\, \mathrm{M_\odot }$ and follows stellar winds, radiation and supernova-feedback of single massive stars with evolution tracks. Star clusters form on time-scales less than ~5 Myr, and their supernova-material is very inefficiently recycled. Stellar wind-material, however, is trapped in massive clusters resulting in the formation of stars self-enriched in Na, Al, and N within only a few Myr. Wind-enriched (second population, 2P) stars can be centrally concentrated in the most massive clusters ($\gtrsim 10^4\, \mathrm{M_\odot }$) and the locked wind-material increases approximately as $M_\mathrm{cluster}^{2}$. These trends resemble the characteristics of observed 2P stars in globular clusters (GCs). We fit scaling relations to the lognormal distributed wind-mass fractions and extrapolate to possible GC progenitors of $M_\mathrm{cluster}=10^7\, \mathrm{M_\odot }$ to investigate whether a dominant 2P could form. This can only happen if the IMF is well-sampled, single massive stars produce at least a factor of a few more enriched winds, for example, through a top-heavy IMF, and a significant fraction of the first population (unenriched) stars is lost during cluster evolution.

The forward-modelling of galaxy surveys has recently gathered interest as one of the primary methods to achieve the precision on the estimate of the redshift distributions required by stage IV surveys. One of the key aspects of forward-modelling is the connection between the physical properties of galaxies and their intrinsic spectral energy distributions (SEDs), achieved through stellar population synthesis (SPS) codes, e.g. FSPS. However, SPS requires many detailed assumptions about the galaxy constituents, for which the model choice or parameters are currently uncertain. In this work, we perform a sensitivity study of the impact that the SED modelling choices variations have on the mean and scatter of the tomographic galaxy redshift distributions. We use the Prospector-$\beta$ model and its SPS parameters to build observed magnitudes of a fiducial sample of galaxies. We then build new samples by varying one SED modelling choice at a time. We model the colour-redshift relation of these galaxy samples using the KiDS-VIKING remapped version (McCullough et al. 2023) of the Masters et al. (2015) SOM. We place galaxies in the SOM cells according to the simulated galaxy colours. We then build color-selected tomographic bins and compare each variant's binned redshift distributions against the estimates obtained for the fiducial model. We find that the SED components related to the IMF, AGN, gas physics, and attenuation law substantially bias the mean and the scatter of the tomographic redshift distributions with respect to those estimated with the fiducial model. For the uncertainty of these choices currently present in the literature, and regardless of any stellar mass function based reweighting strategy applied, the bias in the mean and the scatter of the tomographic redshift distributions is larger than the precision requirements set by Stage IV galaxy surveys, e.g. LSST and Euclid.

Context. Resolved observations at near-infrared (near-IR) and millimeter wavelengths have revealed a diverse population of planet-forming disks. In particular, near-IR scattered light observations usually target close-by, low-mass star-forming regions. However, disk evolution in high-mass star-forming regions is likely affected by the different environment. Orion is the closest high-mass star-forming region, enabling resolved observations to be undertaken in the near-IR.
Aims: We seek to examine planet-forming disks, in scattered light, within the high-mass star-forming region of Orion in order to study the impact of the environment in a higher-mass star-forming region on disk evolution.
Methods: We present SPHERE/IRDIS H-band data for a sample of 23 stars in the Orion star-forming region observed within the DESTINYS (Disk Evolution Study Through Imaging of Nearby Young Stars) program. We used polarization differential imaging in order to detect scattered light from circumstellar dust. From the scattered light observations we characterized the disk orientation, radius, and contrast. We analysed the disks in the context of the stellar parameters and the environment of the Orion star-forming region. We used ancillary X-shooter spectroscopic observations to characterize the central stars in the systems. We furthermore used a combination of new and archival ALMA mm-continuum photometry to characterize the dust masses present in the circumstellar disks.
Results: Within our sample, we detect extended circumstellar disks in ten of 23 systems. Of these, three are exceptionally extended (V351 Ori, V599 Ori, and V1012 Ori) and show scattered light asymmetries that may indicate perturbations by embedded planets or (in the case of V599 Ori) by an outer stellar companion. Our high-resolution imaging observations are also sensitive to close (sub)stellar companions and we detect nine such objects in our sample, of which six were previously unknown. We find in particular a possible substellar companion (either a very low-mass star or a high-mass brown dwarf) 137 au from the star RY Ori. We find a strong anticorrelation between disk detection and multiplicity, with only two of our ten disk detections located in stellar multiple systems. We also find a correlation between scattered light contrast and the millimeter flux. This trend is not captured by previous studies of a more diversified sample and is due to the absence of extended, self-shadowed disks in our Orion sample. Conversely, we do not find significant correlations between the scattered light contrast of the disks and the stellar mass or age. We investigate the radial extent of the disks and compare this to the estimated far-ultraviolet (FUV) field strength at the system location. While we do not find a direct correlation, we notice that no extended disks are detected above an FUV field strength of ~300 G₀.

The tip of the red giant branch (TRGB) based distance method in the I band is one of the most efficient and precise techniques for measuring distances to nearby galaxies (D ≲ 15 Mpc). The TRGB in the near-infrared (NIR) is 1–2 mag brighter relative to the I band, and has the potential to expand the range over which distance measurements to nearby galaxies are feasible. Using Hubble Space Telescope (HST) imaging of 12 fields in eight nearby galaxies, we determine color-based corrections and zero-points of the TRGB in the Wide Field Camera 3 IR (WFC3/IR) F110W and F160W filters. First, we measure TRGB distances in the I band equivalent Advanced Camera System (ACS) F814W filter from resolved stellar populations with the HST. The TRGB in the ACS F814W filter is used for our distance anchor and to place the WFC3/IR magnitudes on an absolute scale. We then determine the color dependence (a proxy for metallicity/age) and zero-point of the NIR TRGB from photometry of WFC3/IR fields that overlap with the ACS fields. The new calibration is accurate to ∼1% in distance relative to the F814W TRGB. Validating the accuracy of the calibrations, we find that the distance modulus for each field using the NIR TRGB calibration agrees with the distance modulus of the same field as determined from the F814W TRGB. This is a JWST preparatory program, and the work done here will directly inform our approach to calibrating the TRGB in JWST NIRCam and NIRISS photometric filters.

Lipids can spontaneously assemble into vesicle-forming membranes. Such vesicles serve as compartments for even the simplest living systems. Vesicles have been extensively studied for constructing synthetic cells or as models for protocells—the cells hypothesized to have existed before life. These compartments exist almost always close to equilibrium. Life, however, exists out of equilibrium. In this work, we studied vesicle-based compartments regulated by a non-equilibrium chemical reaction network that converts activating agents. Specifically, we use activating agents to condense carboxylates and phosphate esters into acylphosphate-based lipids that form vesicles. These vesicles can only be sustained when condensing agents are present, and without them, they decay. We demonstrate that the chemical reaction network can operate on prebiotic activating agents, opening the door to prebiotically plausible, self-sustainable protocells that compete for resources. In future work, such protocells should be endowed with a genotype, for example, based on self-replicating RNA structures that affect the protocell behavior to enable Darwinian evolution in a prebiotically plausible chemical system.

We present a novel cryogenic VUV spectrofluorometer designed to characterize wavelength shifters (WLS) crucial for experiments based on liquid argon (LAr) scintillation light detection. Wavelength shifters like 1,1,4,4-tetraphenyl-1,3-butadiene (TPB) or polyethylene naphthalate (PEN) are used in these experiments to shift the VUV scintillation light to the visible region. Precise knowledge of the optical properties of the WLS at liquid argon's temperature (87 K) and LAr scintillation wavelength (128 nm) is necessary to model and understand the detector response. The cryogenic VUV spectrofluorometer was commissioned to measure the emission spectra and relative wavelength shifting efficiency (WLSE) of samples between 300 K to 87 K for VUV (120 nm to 190 nm) and UV (310 nm) excitation. New mitigation techniques for surface effects on cold WLS were established. As part of this work, the TPB-based wavelength shifting reflector (WLSR) featured in the neutrinoless double-beta decay experiment LEGEND-200 was characterized. The WLSE was observed to increase by (54 ± 5) % from room temperature (RT) to 87 K. PEN installed in LEGEND-200 was also characterized, and a first measurement of the relative WLSE and emission spectrum at RT and 87 K is presented. The WLSE of amorphous PEN was found to be enhanced by at least (37 ± 4) % for excitation with 128 nm and by (52 ± 3) % for UV excitation at 87 K compared to RT.

Over the last decade, evidence has accumulated that massive stars do not typically evolve in isolation but instead follow a tumultuous journey with a companion star on their way to core collapse. While Roche-lobe overflow appears instrumental for the production of a large fraction of Type Ib and Ic supernovae (SNe), variations in the initial orbital period, P_init, of massive interacting binaries may also produce a wide diversity of case B, BC, or C systems, with pre-SN stars endowed from minute to massive H-rich envelopes. Focusing here on the explosion of the primary donor star, originally 12.6 M_⊙, we used radiation hydrodynamics and nonlocal thermodynamic equilibrium time-dependent radiative transfer to document the gas and radiation properties of such SNe, covering Types Ib, IIb, II-L, and II-P. Variations in P_init are the root cause of the wide diversity of our SN light curves, which present single-peak, double-peak, fast-declining, or plateau-like morphologies in the V band. The different ejecta structures, expansion rates, and relative abundances (e.g., H, He, and ⁵⁶Ni) can lead to a great deal of diversity in terms of spectral line shapes (absorption versus emission strength and width) and evolution. We emphasize that Hα is a key tracer of these modulations, and that He I 7065 Å is an enduring optical diagnostic for the presence of He. Our grid of simulations fares well against representative Type Ib, IIb, and II-P SNe, but interaction with circumstellar material, which is ignored in this work, is likely at the origin of the tension between our Type II-L SN models and observations (e.g., of SN 2006Y). Remaining discrepancies in the rise time to bolometric maximum of our models call for a proper account of both small-scale and large-scale structures in core-collapse SN ejecta. Discrepant Type II-P SN models, with a high plateau brightness but small spectral line widths, can be fixed by adopting more compact red-supergiant star progenitors.

Context. Very high-energy (VHE, E > 100 GeV) observations of the blazar Mrk 501 with MAGIC in 2014 provided evidence for an unusual narrow spectral feature at about 3 TeV during an extreme X-ray flaring activity. The one-zone synchrotron-self Compton scenario, widely used in blazar broadband spectral modeling, fails to explain the narrow TeV component.
Aims: Motivated by this rare observation, we propose an alternative model for the production of narrow features in the VHE spectra of flaring blazars. These spectral features may result from the decay of neutral pions (π⁰ bumps) that are in turn produced via interactions of protons (of tens of TeV energy) with energetic photons, whose density increases during hard X-ray flares.
Methods: We explored the conditions needed for the emergence of narrow π⁰ bumps in VHE blazar spectra during X-ray flares reaching synchrotron energies ∼100 keV using time-dependent radiative transfer calculations. We focused on high-synchrotron peaked (HSP) blazars, which comprise the majority of VHE-detected extragalactic sources.
Results: We find that synchrotron-dominated flares with peak energies ≳100 keV can be ideal periods for the search of π⁰ bumps in the VHE spectra of HSP blazars. The flaring region is optically thin to photopion production, its energy content is dominated by the relativistic proton population, and the inferred jet power is highly super-Eddington. Application of the model to the spectral energy distribution of Mrk 501 on MJD 56857.98 shows that the VHE spectrum of the flare is described well by the sum of a synchrotron self-Compton (SSC) component and a distinct π⁰ bump centered at 3 TeV. Spectral fitting of simulated SSC+π⁰ spectra for the Cherenkov Telescope Array (CTA) show that a π⁰ bump could be detected at a 5σ significance level with a 30-min exposure.
Conclusions: A harder VHE γ-ray spectrum than the usual SSC prediction or, more occasionally, a distinct narrow bump at VHE energies during hard X-ray flares, can be suggestive of a relativistic hadronic component in blazar jets that otherwise would remain hidden. The production of narrow features or spectral hardenings due to π⁰ decay in the VHE spectra of blazars is testable with the advent of CTA.

We present the results from a complex study of an eclipsing O-type binary (Aa+Ab) with the orbital period of P _A = 3.2254367 days that forms part of a higher-order multiple system in a configuration of (A+B)+C. We derived masses of the Aa+Ab binary of M ₁ = 19.02 ± 0.12 and M ₂ = 17.50 ± 0.13 M _⊙, the radii of R ₁ = 7.70 ± 0.05 and R ₂ = 6.64 ± 0.06 R _⊙, and temperatures of T ₁ = 34,250 ± 500 K and T ₂ = 33,750 ± 500 K. From the analysis of the radial velocities, we found a spectroscopic orbit of A in the outer A+B system with P _A+B = 195.8 days (P _A+B/P _A ≈ 61). In the O ‑ C analysis, we confirmed this orbit and found another component orbiting the A+B system with P _AB+C = 2550 days (P _AB+C/P _A+B ≈ 13). From the total mass of the inner binary and its outer orbit, we estimated the mass of the third object, M _B ≳ 10.7 M _⊙. From the light travel time effect fit to the O ‑ C data, we obtained the limit for the mass of the fourth component, M _C ≳ 7.3 M _⊙. These extra components contribute about 20%–30% (increasing with wavelength) to the total system light. From the comparison of model spectra with the multiband photometry, we derived a distance modulus of 18.59 ± 0.06 mag, a reddening of 0.16 ± 0.02 mag, and an R_V of 3.2. This work is part of our ongoing project, which aims to calibrate the surface brightness–color relation for early-type stars. *Based on observations collected at the European Southern Observatory, Chile. ^† This paper includes data gathered with the 6.5 m Magellan Clay Telescope at Las Campanas Observatory, Chile.

We compute differential distributions for Drell-Yan processes at the LHC and the Tevatron colliders at next-to-next-to-leading order in perturbative QCD, including fiducial cuts on the decay leptons in the final state. The comparison of predictions obtained with four different codes shows excellent agreement, once linear power corrections from the fiducial cuts are included in those codes that rely on phase-space slicing subtraction schemes. For $Z$-boson production we perform a detailed study of the symmetric cuts on the transverse momenta of the decay leptons. Predictions at fixed order in perturbative QCD for those symmetric cuts, typically imposed in experiments, suffer from an instability. We show how this can be remedied by an all-order resummation of the fiducial transverse momentum spectrum, and we comment on the choice of cuts for future experimental analyses.

We investigate experimentally the dynamic phase transition of a two-dimensional active nematic layer interfaced with a passive liquid crystal. Under a temperature ramp that leads to the transition of the passive liquid into a highly anisotropic lamellar smectic-A phase, and in the presence of a magnetic field, the coupled active nematic reorganizes its flow and orientational patterns from the turbulent into a quasilaminar regime aligned perpendicularly to the field. Remarkably, while the phase transition of the passive fluid is known to be continuous, or second order, our observations reveal intermittent dynamics of the order parameter and the coexistence of aligned and turbulent regions in the active nematic, a signature of discontinuous, or first order, phase transitions, similar to what is known to occur in relation to flocking in dry active matter. Our results suggest that alignment transitions in active systems are intrinsically discontinuous, regardless of the symmetry and momentum-damping mechanisms.

Cryogenic phonon detectors with transition-edge sensors achieve the best sensitivity to sub-GeV/c² dark matter interactions with nuclei in current direct detection experiments. In such devices, the temperature of the thermometer and the bias current in its readout circuit need careful optimization to achieve optimal detector performance. This task is not trivial and is typically done manually by an expert. In our work, we automated the procedure with reinforcement learning in two settings. First, we trained on a simulation of the response of three Cryogenic Rare Event Search with Superconducting Thermometers (CRESST) detectors used as a virtual reinforcement learning environment. Second, we trained live on the same detectors operated in the CRESST underground setup. In both cases, we were able to optimize a standard detector as fast and with comparable results as human experts. Our method enables the tuning of large-scale cryogenic detector setups with minimal manual interventions.

The recently reported observation of VFTS 243 is the first example of a massive black-hole binary system with negligible binary interaction following black-hole formation. The black-hole mass (≈10 M_⊙) and near-circular orbit (e ≈0.02 ) of VFTS 243 suggest that the progenitor star experienced complete collapse, with energy-momentum being lost predominantly through neutrinos. VFTS 243 enables us to constrain the natal kick and neutrino-emission asymmetry during black-hole formation. At 68% confidence level, the natal kick velocity (mass decrement) is ≲10 km /s (≲1.0 M_⊙ ), with a full probability distribution that peaks when ≈0.3 M_⊙ were ejected, presumably in neutrinos, and the black hole experienced a natal kick of 4 km /s . The neutrino-emission asymmetry is ≲4 %, with best fit values of ∼0 - 0.2 % . Such a small neutrino natal kick accompanying black-hole formation is in agreement with theoretical predictions.

The Euclid photometric survey of galaxy clusters stands as a powerful cosmological tool, with the capacity to significantly propel our understanding of the Universe. Despite being subdominant to dark matter and dark energy, the baryonic component of our Universe holds substantial influence over the structure and mass of galaxy clusters. This paper presents a novel model that can be used to precisely quantify the impact of baryons on the virial halo masses of galaxy clusters using the baryon fraction within a cluster as a proxy for their effect. Constructed on the premise of quasi-adiabaticity, the model includes two parameters, which are calibrated using non-radiative cosmological hydrodynamical simulations, and a single large-scale simulation from the Magneticum set, which includes the physical processes driving galaxy formation. As a main result of our analysis, we demonstrate that this model delivers a remarkable 1% relative accuracy in determining the virial dark matter-only equivalent mass of galaxy clusters starting from the corresponding total cluster mass and baryon fraction measured in hydrodynamical simulations. Furthermore, we demonstrate that this result is robust against changes in cosmological parameters and against variation of the numerical implementation of the subresolution physical processes included in the simulations. Our work substantiates previous claims regarding the impact of baryons on cluster cosmology studies. In particular, we show how neglecting these effects would lead to biased cosmological constraints for a Euclid-like cluster abundance analysis. Importantly, we demonstrate that uncertainties associated with our model arising from baryonic corrections to cluster masses are subdominant when compared to the precision with which mass-observable (i.e. richness) relations will be calibrated using Euclid and to our current understanding of the baryon fraction within galaxy clusters.

Narrowband galaxy surveys have recently gained interest as a promising method to achieve the necessary accuracy on the photometric redshift estimate of individual galaxies for stage-IV cosmological surveys. One key advantage is the ability to provide higher spectral resolution information about galaxies that should allow a more accurate and precise estimation of galaxy stellar population properties. However, the impact of adding narrow-band photometry on the stellar population properties estimate is largely unexplored. The scope of this work is two-fold: on one side, leveraging the predictive power of broad-band and narrow-band data to infer galaxy physical properties such as stellar masses, ages, star formation rates and metallicities. On the other hand, evaluating the improvement of performance in estimating galaxy properties when we use narrow-band data instead of broad-band. In this work we measure the stellar population properties of a sample of galaxies in the COSMOS field for which both narrowband and broadband data are available. In particular, we employ narrowband data from PAUS and broad-band data from CFHTLS. We use two different spectral energy distribution fitting codes to measure galaxy properties, namely CIGALE and Prospector. We find that the increased spectral resolution of narrow-band photometry does not yield a substantial improvement on constraining galaxy properties using spectral energy distribution fitting. Still we find that we obtain a more diverse distribution of metallicities and dust optical depths with cigale when employing the narrowband data. The effect is not as prominent as expected, which we relate this to the low narrowband SNR of a majority of the galaxies, the respective drawbacks of both codes as well as the coverage only in the optical regime. The measured properties are afterwards compared to the COSMOS2020 catalogue, showing good agreement.

We calculate directly in position space the one-loop renormalization kernels of the soft operators O_γ and O_g that appear in the soft-quark contributions to, respectively, the subleading-power γγ → h and gg → h form factors mediated by the b-quark. We present an IR/rapidity divergence-free definition for O_g and demonstrate that with a correspondent definition of the collinear function, a consistent factorization theorem is recovered. Using conformal symmetry techniques, we establish a relation between the evolution kernels of the leading-twist heavy-light light-ray operator, whose matrix element defines the B-meson light-cone distribution amplitude (LCDA), and O_γ to all orders in perturbation theory. Application of this relation allows us to bootstrap the kernel of O_γ to the two-loop level. We construct an ansatz for the kernel of O_g at higher orders. We test this ansatz against the consistency requirement at two-loop and find they differ only by a particular constant.

We study the effects of exceptionally light QCD axions on the stellar configuration of white dwarfs. At finite baryon density, the nonderivative coupling of the axion to nucleons displaces the axion from its in-vacuum minimum, which implies a reduction of the nucleon mass. This dramatically alters the composition of stellar remnants. In particular, the modifications of the mass-radius relationship of white dwarfs allow us to probe large regions of unexplored parameter space without requiring that axions are dark matter.

We have performed a systematic search for galaxy-scale strong lenses using Hyper Suprime-Cam imaging data, focusing on lenses in overdense environments. To identify these lens candidates, we exploit our neural network from HOLISMOKES VI, which is trained on realistic gri mock-images as positive examples, and real images as negative examples. Compared to our previous work, we lower the i-Kron radius limit to >0.5". This results in an increase by around 73 million sources to more than 135 million images. During our visual multi-stage grading of the network candidates, we now inspect simultaneously larger stamps (80"x80") to identify large, extended arcs cropped in the 10"x10" cutouts, and classify additionally their overall environment. Here we also reinspect our previous lens candidates and classify their environment. Using these 546 visually identified lens candidates, we further define various criteria by exploiting extensive and complementary photometric redshift catalogs, to select the candidates in overdensities. In total, we identified 24 grade-A and 138 grade-B candidates with either spatially-resolved multiple images or extended, distorted arcs in the new sample. Furthermore, with our different techniques, we identify in total 237/546 lens candidates in a cluster-like or overdense environment, containing only 49 group- or cluster-scale re-discoveries. These results demonstrate the feasibility of downloading and applying network classifiers to hundreds of million cutouts, necessary in the upcoming era of big data from deep, wide-field imaging surveys like Euclid and the Rubin Observatory Legacy Survey of Space and Time, while leading to a sample size that can be inspected by humans. These networks, with false-positive rates of ~0.01%, are very powerful tools to identify such rare galaxy-scale strong lensing systems, while also aiding in the discovery of new strong lensing clusters.

Neutrino-nucleus scatterings in the detector could induce electron ionization signatures due to the Migdal effect. We derive prospects for a future detection of the Migdal effect via coherent elastic solar neutrino-nucleus scatterings in liquid xenon detectors, and discuss the irreducible background that it constitutes for the Migdal effect caused by light dark matter-nucleus scatterings. Furthermore, we explore the ionization signal induced by some neutrino electromagnetic and non-standard interactions on nuclei. In certain scenarios, we find a distinct peak on the ionization spectrum of xenon around 0.1 keV, in clear contrast to the Standard Model expectation.

Context. The mass and spin of massive black holes (BHs) at the centre of galaxies evolve due to gas accretion and mergers with other BHs. Besides affecting the evolution of relativistic jets, for example, the BH spin determines the efficiency with which the BH radiates energy.
Aims: Using cosmological, hydrodynamical simulations, we investigate the evolution of the BH spin across cosmic time and its role in controlling the joint growth of supermassive BHs and their host galaxies.
Methods: We implemented a sub-resolution prescription that models the BH spin, accounting for both BH coalescence and misaligned accretion through a geometrically thin, optically thick disc. We investigated how BH spin evolves in two idealised setups, in zoomed-in simulations and in a cosmological volume. The latter simulation allowed us to retrieve statistically robust results for the evolution and distribution of BH spins as a function of BH properties.
Results: We find that BHs with M_BH ≲ 2 × 10⁷ M_⊙ grow through gas accretion, occurring mostly in a coherent fashion that favours spin-up. Above M_BH ≳ 2 × 10⁷ M_⊙, the gas angular momentum directions of subsequent accretion episodes are often uncorrelated with each other. The probability of counter-rotating accretion and hence spin-down increases with BH mass. In the latter mass regime, BH coalescence plays an important role. The spin magnitude displays a wide variety of histories, depending on the dynamical state of the gas feeding the BH and the relative contribution of mergers and gas accretion. As a result of their combined effect, we observe a broad range of values of the spin magnitude at the high-mass end. Reorientation of the BH spin direction occurs on short timescales (≲ 10 Myr) only during highly accreting phases (ƒ_Edd ≳ 0.1). Our predictions for the distributions of BH spin and spin-dependent radiative efficiency as a function of BH mass are in very good agreement with observations.

Movie associated to Fig. 7 is available at https://www.aanda.org

Aims: We present the first multiwavelength study of Mrk 501 that contains simultaneous very-high-energy (VHE) γ-ray observations and X-ray polarization measurements from the Imaging X-ray Polarimetry Explorer (IXPE).
Methods: We used radio-to-VHE data from a multiwavelength campaign carried out between March 1, 2022, and July 19, 2022 (MJD 59639 to MJD 59779). The observations were performed by MAGIC, Fermi-LAT, NuSTAR, Swift (XRT and UVOT), and several other instruments that cover the optical and radio bands to complement the IXPE pointings. We characterized the dynamics of the broadband emission around the X-ray polarization measurements through its multiband fractional variability and correlations, and compared changes observed in the polarization degree to changes seen in the broadband emission using a multi-zone leptonic scenario.
Results: During the IXPE pointings, the VHE state is close to the average behavior, with a 0.2-1 TeV flux of 20%-50% of the emission of the Crab Nebula. Additionally, it shows low variability and a hint of correlation between VHE γ-rays and X-rays. Despite the average VHE activity, an extreme X-ray behavior is measured for the first two IXPE pointings, taken in March 2022 (MJD 59646 to 59648 and MJD 59665 to 59667), with a synchrotron peak frequency > 1 keV. For the third IXPE pointing, in July 2022 (MJD 59769 to 59772), the synchrotron peak shifts toward lower energies and the optical/X-ray polarization degrees drop. All three IXPE epochs show an atypically low Compton dominance in the γ-rays. The X-ray polarization is systematically higher than at lower energies, suggesting an energy stratification of the jet. While during the IXPE epochs the polarization angles in the X-ray, optical, and radio bands align well, we find a clear discrepancy in the optical and radio polarization angles in the middle of the campaign. Such results further support the hypothesis of an energy-stratified jet. We modeled broadband spectra taken simultaneous to the IXPE pointings, assuming a compact zone that dominates in the X-rays and the VHE band, and an extended zone stretching farther downstream in the jet that dominates the emission at lower energies. NuSTAR data allow us to precisely constrain the synchrotron peak and therefore the underlying electron distribution. The change between the different states observed in the three IXPE pointings can be explained by a change in the magnetization and/or the emission region size, which directly connects the shift in the synchrotron peak to lower energies with the drop in the polarization degree.

The MWL data are available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (ftp://130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/685/A117

Supernova (SN) SN H0pe is a gravitationally lensed, triply imaged, Type Ia SN (SN Ia) discovered in James Webb Space Telescope imaging of the PLCK G165.7+67.0 cluster of galaxies. Well-observed multiply imaged SNe provide a rare opportunity to constrain the Hubble constant (H ₀), by measuring the relative time delay between the images and modeling the foreground mass distribution. SN H0pe is located at z = 1.783 and is the first SN Ia with sufficient light-curve sampling and long enough time delays for an H ₀ inference. Here we present photometric time-delay measurements and SN properties of SN H0pe. Using JWST/NIRCam photometry, we measure time delays of Δt _ab = <inline-formula> <tex-math> $-{116.6}_{-9.3}^{+10.8}$ </tex-math> </inline-formula> observer-frame days and Δt _cb = <inline-formula> <tex-math> $-{48.6}_{-4.0}^{+3.6}$ </tex-math> </inline-formula> observer-frame days relative to the last image to arrive (image 2b; all uncertainties are 1σ), which corresponds to a ∼5.6% uncertainty contribution for H ₀ assuming 70 km s^‑1 Mpc^‑1. We also constrain the absolute magnification of each image to μ _a = <inline-formula> <tex-math> ${4.3}_{-1.8}^{+1.6}$ </tex-math> </inline-formula>, μ _b = <inline-formula> <tex-math> ${7.6}_{-2.6}^{+3.6}$ </tex-math> </inline-formula>, μ _c = <inline-formula> <tex-math> ${6.4}_{-1.5}^{+1.6}$ </tex-math> </inline-formula> by comparing the observed peak near-IR magnitude of SN H0pe to the nonlensed population of SNe Ia.

The Draco Dwarf spheroidal (dSph) galaxy is one of the nearest and the most dark-matter-dominated satellites of the Milky Way. We obtained multiepoch near-infrared (NIR, JHK _s ) observations of the central region of Draco dSph covering a sky area of ∼21' × 21' using the WIRCam instrument at the 3.6 m Canada–France–Hawaii Telescope. Homogeneous JHK _s time-series photometry for 212 RR Lyrae (173 fundamental-mode, 24 first-overtone, and 15 mixed-mode variables) and five Anomalous Cepheids in Draco dSph are presented and used to derive their period–luminosity relations at NIR wavelengths for the first-time. The small scatter of ∼0.05 mag in these empirical relations for RR Lyrae stars is consistent with those in globular clusters and suggests a very small metallicity spread, up to ∼0.2 dex, among these centrally located variables. Based on empirically calibrated NIR period–luminosity–metallicity relations for RR Lyrae in globular clusters, we determined a distance modulus to Draco dSph of μ _RRL = 19.557 ± 0.026 mag. The calibrated K _s -band period–luminosity relations for Anomalous Cepheids in the Draco dSph and the Large Magellanic Cloud exhibit statistically consistent slopes but systematically different zero points, hinting at possible metallicity dependence of ∼ ‑ 0.3 mag dex^‑1. Finally, the apparent magnitudes of the tip of the red-giant branch in I and J bands also agree well with their absolute calibrations with the adopted RR Lyrae distance to Draco. Our recommended ∼1.5% precise RR Lyrae distance, D _Draco = 81.55 ± 0.98(statistical) ± 1.17(systematic) kpc, is the most accurate and precise distance to Draco dSph galaxy.

We compute all helicity amplitudes for the scattering of five partons in two-loop QCD in all the relevant flavor configurations, retaining all contributing color structures. We employ tensor projection to obtain helicity amplitudes in the 't Hooft-Veltman scheme starting from a set of primitive amplitudes. Our analytic results are expressed in terms of massless pentagon functions, and are easy to evaluate numerically. These amplitudes provide important input to investigations of soft-collinear factorization and to studies of the high-energy limit.

Supernovae (SNs) are an important source of energy in the interstellar medium. Young remnants of supernovae (SNRs) exhibit peak emission in the X-ray region, making them interesting objects for X-ray observations. In particular, the supernova remnant SN1006 is of great interest due to its historical record, proximity, and brightness. Thus, it has been studied with a number of X-ray telescopes. Improving X-ray imaging of this and other remnants is an important but challenging task, as it often requires multiple observations with different instrument responses to image the entire object. Here, we use Chandra observations to demonstrate the capabilities of Bayesian image reconstruction using information field theory (IFT). Our objective is to reconstruct denoised, deconvolved, and spatio-spectral resolved images from X-ray observations and to decompose the emission into different morphologies, namely, diffuse and point-like. Further, we aim to fuse data from different detectors and pointings into a mosaic and quantify the uncertainty of our result. By utilizing prior knowledge on the spatial and spectral correlation structure of the diffuse emission and point sources, this method allows for the effective decomposition of the signal into these two components. In order to accelerate the imaging process, we introduced a multi-step approach, in which the spatial reconstruction obtained for a single energy range is used to derive an informed starting point for the full spatio-spectral reconstruction. We applied this method to 11 Chandra observations of SN1006 from 2008 and 2012, providing a detailed, denoised, and decomposed view of the remnant. In particular, the separated view of the diffuse emission ought to provide new insights into the complex, small-scale structures in the center of the remnant and at the shock front profiles. For example, our analysis reveals sharp X-ray flux increases by up to two orders of magnitude at the shock fronts of SN1006.

Stars can be tidally disrupted when passing near a black hole, and the debris can induce a flux of high-energy neutrinos. It has been discussed that there are hints in IceCube data of high-energy neutrinos produced in Tidal Disruption Events. The emitting region of neutrinos and photons in these astrophysical events is likely to be located in the vicinity of the central black hole, where the dark matter density might be significantly larger than in the outer regions of the galaxy. We explore the potential attenuation of the emitted neutrino and photon fluxes due to interactions with dark matter particles around the supermassive black hole of the host galaxies of AT2019dsg, AT2019fdr and AT2019aalc, and study the implications for some well-motivated models of dark matter-neutrino and dark matter-photon interactions. Furthermore, we discuss the complementarity of our constraints with values of the dark matter-neutrino scattering cross section proven to alleviate some cosmological tensions.

We present photometric measurements of 88 Cepheid variables in the core of the Small Magellanic Cloud (SMC), the first sample obtained with the Hubble Space Telescope (HST) and Wide Field Camera 3, in the same homogeneous photometric system as past measurements of all Cepheids on the SH0ES distance ladder. We limit the sample to the inner core and model the geometry to reduce errors in prior studies due to the non-trivial depth of this Cloud. Without crowding present in ground-based studies, we obtain an unprecedentedly low dispersion of 0.102 mag for a Period-Luminosity relation in the SMC, approaching the width of the Cepheid instability strip. The new geometric distance to 15 late-type detached eclipsing binaries in the SMC offers a rare opportunity to improve the foundation of the distance ladder, increasing the number of calibrating galaxies from three to four. With the SMC as the only anchor, we find H$_0\!=\!74.1 \pm 2.1$ km s$^{-1}$ Mpc$^{-1}$. Combining these four geometric distances with our HST photometry of SMC Cepheids, we obtain H$_0\!=\!73.17 \pm 0.86$ km s$^{-1}$ Mpc$^{-1}$. By including the SMC in the distance ladder, we also double the range where the metallicity ([Fe/H]) dependence of the Cepheid Period-Luminosity relation can be calibrated, and we find $\gamma = -0.22 \pm 0.05$ mag dex$^{-1}$. Our local measurement of H$_0$ based on Cepheids and Type Ia supernovae shows a 5.8$\sigma$ tension with the value inferred from the CMB assuming a $\Lambda$CDM cosmology, reinforcing the possibility of physics beyond $\Lambda$CDM.

We present UV/optical/NIR observations and modeling of supernova (SN) 2024ggi, a type II supernova (SN II) located in NGC 3621 at 7.2 Mpc. Early-time ("flash") spectroscopy of SN 2024ggi within +0.8 days of discovery shows emission lines of H I, He I, C III, and N III with a narrow core and broad, symmetric wings (i.e., IIn-like) arising from the photoionized, optically-thick, unshocked circumstellar material (CSM) that surrounded the progenitor star at shock breakout. By the next spectral epoch at +1.5 days, SN 2024ggi showed a rise in ionization as emission lines of He II, C IV, N IV/V and O V became visible. This phenomenon is temporally consistent with a blueward shift in the UV/optical colors, both likely the result of shock breakout in an extended, dense CSM. The IIn-like features in SN 2024ggi persist on a timescale of $t_{\rm IIn} = 3.8 \pm 1.6$ days at which time a reduction in CSM density allows the detection of Doppler broadened features from the fastest SN material. SN 2024ggi has peak UV/optical absolute magnitudes of $M_{\rm w2} = -18.7$ mag and $M_{\rm g} = -18.1$ mag that are consistent with the known population of CSM-interacting SNe II. Comparison of SN 2024ggi with a grid of radiation hydrodynamics and non-local thermodynamic equilibrium (nLTE) radiative-transfer simulations suggests a progenitor mass-loss rate of $\dot{M} = 10^{-2}$M$_{\odot}$ yr$^{-1}$ ($v_w$ = 50 km/s), confined to a distance of $r < 5\times 10^{14}$ cm. Assuming a wind velocity of $v_w$ = 50 km/s, the progenitor star underwent an enhanced mass-loss episode in the last ~3 years before explosion.

The emergence of biopolymer building blocks is a crucial step during the origins of life. However, all known formation pathways rely on rare pure feedstocks and demand successive purification and mixing steps to suppress unwanted side reactions and enable high product yields. Here we show that heat flows through thin, crack-like geo-compartments could have provided a widely available yet selective mechanism that separates more than 50 prebiotically relevant building blocks from complex mixtures of amino acids, nucleobases, nucleotides, polyphosphates and 2-aminoazoles. Using measured thermophoretic properties, we numerically model and experimentally prove the advantageous effect of geological networks of interconnected cracks that purify the previously mixed compounds, boosting their concentration ratios by up to three orders of magnitude. The importance for prebiotic chemistry is shown by the dimerization of glycine, in which the selective purification of trimetaphosphate (TMP) increased reaction yields by five orders of magnitude. The observed effect is robust under various crack sizes, pH values, solvents and temperatures. Our results demonstrate how geologically driven non-equilibria could have explored highly parallelized reaction conditions to foster prebiotic chemistry.

We present the results from a complex study of an eclipsing O-type binary (Aa+Ab) with the orbital period $P_{A}=3.2254367$ days, that forms part of a higher-order multiple system in a configuration (A+B)+C. We derived masses of the Aa+Ab binary $M_{1}= 19.02 \pm 0.12 \,M_\odot$, $M_{2}= 17.50 \pm 0.13 \,M_\odot$, radii $R_{1}= 7.70 \pm 0.05 \,R_\odot$, $R_{2}= 6.64 \pm 0.06 \,R_\odot$, and temperatures $T_1 = 34250 \pm 500 $ K, $T_2 = 33750 \pm 500 $ K. From the analysis of radial velocities, we found a spectroscopic orbit of A in the outer A+B system with $P_{A+B}=195.8$ days ($P_{A+B}/P_{A}\approx 61$). In the O-C analysis, we confirmed this orbit and found another component orbiting the A+B system with $P_{AB+C}=2550$ days ($P_{AB+C}\,/P_{A+B}\approx 13$). From the total mass of the inner binary and its outer orbit, we estimated the mass of the third object, $M_B \gtrsim 10.7 M_\odot$. From the light-travel time effect fit to the O-C data, we obtained the limit for the mass of the fourth component, $M_C \gtrsim 7.3 M_\odot$. These extra components contribute to about 20% to 30% (increasing with wavelength) of the total system light. From the comparison of model spectra with the multiband photometry, we derived a distance modulus of 18.59 $\pm$ 0.06 mag, a reddening of 0.16 $\pm$ 0.02 mag, and an $R_V$ of $3.2$. This work is part of our ongoing project, which aims to calibrate the surface brightness-color relation for early-type stars.