This study investigates the clustering and bias of Luminous Red Galaxies (LRG) in the BOSS-LOWZ, -CMASS, -COMB, and eBOSS samples, using two types of simulated lightcones: (i) high-fidelity lightcones from UCHUUN-body simulation, employing SHAM technique to assign LRG to (sub)haloes, and (ii) 16 000 covariance lightcones from GLAM-UCHUUN-body simulations, including LRG using HOD data from UCHUU. Our results indicate that UCHUU and GLAM lightcones closely replicate BOSS/eBOSS data, reproducing correlation function and power spectrum across scales from redshifts 0.2 to 1.0, from 2 to $150 \,h^{-1}\,\mathrm{Mpc}$ in configuration space, from 0.005 to $0.7\, h\,\mathrm{Mpc}^{-1}$ in Fourier space, and across different LRG stellar masses. Furthermore, comparing with existing MD-PATCHY and EZMOCK BOSS/eBOSS lightcones based on approximate methods, our GLAM-UCHUU lightcones provide more precise clustering estimates. We identify significant deviations from observations within $20 \,h^{-1}\,\mathrm{Mpc}$ scales in MD-PATCHY and EZMOCK, with our covariance matrices indicating that these methods underestimate errors by between 10 per cent and 60 per cent. Lastly, we explore the impact of cosmology on galaxy clustering. Our findings suggest that, given the current level of uncertainties in BOSS/eBOSS data, distinguishing models with and without massive neutrino effects on large-scale structure (LSS) is challenging. This paper highlights the UCHUU and GLAM-UCHUU simulations' robustness in verifying the accuracy of Planck cosmological parameters, providing a strong foundation for enhancing lightcone construction in future LSS surveys. We also demonstrate that generating thousands of galaxy lightcones is feasible using N-body simulations with adequate mass and force resolution.

More than 30% of white dwarfs exhibit atmospheric metals, which are understood to be from recent or ongoing accretion of circumstellar debris. In cool white dwarfs, surface motions should rapidly homogenise photospheric abundances, and the accreted heavy elements should diffuse inward on a timescale much longer than that for surface mixing. The recent discovery of a metal scar on WD 0816-310 implies its B ≈ 140 kG magnetic field has impeded surface mixing of metals near the visible magnetic pole. Here, we report the discovery of a second magnetic, metal-polluted white dwarf, WD 2138-332, which exhibits periodic variability in longitudinal field, metal line strength, and broadband photometry. All three variable quantities have the same period, and show remarkable correlations: the published light curves have a brightness minimum exactly when the longitudinal field and line strength have a maximum, and a maximum when the longitudinal field and line strength have a minimum. The simplest interpretation of the line strength variability is that there is an enhanced metal concentration around one pole of the magnetic field; however, the variable line-blanketing cannot account for the observed multi-band light curves. More theoretical work is required to understand the efficiency of horizontal mixing of the accreted metal atoms, and the origin of photometric variability. Because both magnetic, metal-polluted white dwarfs that have been monitored to date show that metal line strengths vary in phase with the longitudinal field, we suggest that metal scars around magnetic poles may be a common feature of metal-polluted white dwarfs.

Frequentist parameter inference using profile likelihoods has received increased attention in the cosmology literature recently since it can give important complementary information to Bayesian credible intervals. Here, we give a pedagogical review to frequentist parameter inference in cosmology with particular focus on when the graphical profile likelihood construction gives meaningful constraints, i.e.\ confidence intervals with correct coverage. This construction rests on the assumption of the asymptotic limit of a large data set such as in \textit{Wilks' theorem}. We assess the validity of this assumption in the context of three cosmological models with \textit{Planck} 2018 \texttt{Plik\_lite} data: While our tests for the $\Lambda$CDM model indicate that the profile likelihood method gives correct coverage, $\Lambda$CDM with the sum of neutrino masses as a free parameter appears consistent with a Gaussian near a boundary motivating the use of the boundary-corrected or Feldman-Cousins graphical method; for $w_0$CDM with the equation of state of dark energy, $w_0$, as a free parameter, we find indication of a violation of the assumptions. Finally, we compare frequentist and Bayesian constraints of these models. Our results motivate care when using the graphical profile likelihood method in cosmology.

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

Galaxy chemical enrichment mechanisms have primarily been constrained by [$\alpha$/Fe] and [Fe/H] measurements of individual stars and integrated light from stellar populations. However such measurements are limited at higher redshifts (z>1). Recently, we proposed an analogous diagram of the oxygen-to-argon abundance ratio, log(O/Ar), vs Ar abundance, 12+log(Ar/H), as a new diagnostic window for emission nebulae. In this Letter, using robust line flux measurements including temperature sensitive auroral lines, we present direct determination of O and Ar abundances in nine SFGs from JWST/NIRSPEC spectra at z$\sim$1.3-7.7, and two more with Keck/MOSFIRE spectra at z$\sim$2.2. Utilising their positions on the log(O/Ar) vs 12+log(Ar/H) plane, we present the first inference of galaxy chemical enrichment mechanisms from an ensemble of galaxies. The SFGs at z$\sim$1.3-3.4 are consistent with the solar neighbourhood galactic chemical enrichment models of the Milky Way Galaxy that are driven by core-collapse and Type Ia supernovae. Such enrichment mechanisms thus occur at least out to z$\sim$3.4. However, the highest-redshift SFGs (z$\sim$3.6-7.7) have very low log(O/Ar) values, revealing a different enrichment process at z>3.6. Such low log(O/Ar) values may be caused by a rapid but intermittent star-formation and/or additional sources. The new diagnostic window for SFGs enables us to reveal the unique fingerprints of galaxy chemical enrichment out to cosmic dawn.

Modern cosmological research in large-scale structure has witnessed an increasing number of machine-learning applications. Among them, convolutional neural networks (CNNs) have received substantial attention due to their outstanding performance in image classification, cosmological parameter inference, and various other tasks. However, many models based on CNNs are criticized as "black boxes" due to the difficulties in relating their outputs intuitively and quantitatively to the cosmological fields under investigation. To overcome this challenge, we present the Cosmological Correlator Convolutional Neural Network (C3NN)—a fusion of CNN architecture and cosmological N-point correlation functions (NPCFs). We demonstrate that its output can be expressed explicitly in terms of the analytically tractable NPCFs. Together with other auxiliary algorithms, we can open the "black box" by quantitatively ranking different orders of the interpretable outputs based on their contribution to classification tasks. As a proof of concept, we demonstrate this by applying our framework to a series of binary classification tasks using Gaussian and log-normal random fields and relating its outputs to the NPCFs describing the two fields. Furthermore, we exhibit the model's ability to distinguish different dark energy scenarios (w ₀ = ‑0.95 and ‑1.05) using N-body simulated weak-lensing convergence maps and discuss the physical implications coming from their interpretability. With these tests, we show that C3NN combines advanced aspects of machine learning architectures with the framework of cosmological NPCFs, thereby making it an exciting tool to extract physical insights in a robust and explainable way from observational data.

For decades, the boundary of cosmic filaments has been a subject of debate. In this work, we determine the physically motivated radii of filaments by constructing stacked galaxy number density profiles around the filament spines. We find that the slope of the profile changes with distance to the filament spine, reaching its minimum at approximately 1 Mpc at $z=0$ in both state-of-the-art hydrodynamical simulations and observational data. This can be taken as the average value of the filament radius. Furthermore, we note that the average filament radius rapidly decreases from $z=4$ to 1, and then slightly increases. Moreover, we find that the radius of the filament depends on the length of the filament, the distance from the connected clusters, and the masses of the clusters. These results suggest a two-phase formation scenario of cosmic filaments. The filaments experienced rapid contraction before $z=1$, but their density distribution has remained roughly stable since then. The subsequent mass transport along the filaments to the connected clusters is likely to have contributed to the formation of the clusters themselves.

Context. Recent observations suggest that planet formation starts early, in protostellar disks of ≤10⁵ yr, which are characterized by strong interactions with the environment, such as through accretion streamers and molecular outflows.
Aims: To investigate the impact of such phenomena on the physical and chemical properties of a disk, it is key to understand what chemistry planets inherit from their natal environment.
Methods: In the context of the ALMA large program Fifty AU Study of the chemistry in the disk/envelope system of solar-like protostars (FAUST), we present observations on scales from ∼1500 au to ∼60 au of H₂CO, HDCO, and D₂CO toward the young planet-forming disk IRS 63.
Results: The H₂CO probes the gas in the disk as well as in a large scale streamer (∼1500 au) impacting onto the southeast disk side. We detected for the first time deuterated formaldehyde, HDCO and D₂CO, in a planet-forming disk and HDCO in the streamer that is feeding it. These detections allowed us to estimate the deuterium fractionation of H₂CO in the disk: [HDCO]/[H₂CO] ∼ 0.1 − 0.3 and [D₂CO]/[H₂CO] ∼ 0.1. Interestingly, while HDCO follows the H₂CO distribution in the disk and in the streamer, the distribution of D₂CO is highly asymmetric, with a peak of the emission (and [D]/[H] ratio) in the southeast disk side, where the streamer crashes onto the disk. In addition, D₂CO was detected in two spots along the blue- and redshifted outflow. This suggests that (i) in the disk, HDCO formation is dominated by gas-phase reactions in a manner similar to H₂CO, while (ii) D₂CO is mainly formed on the grain mantles during the prestellar phase and/or in the disk itself and is at present released in the gas phase in the shocks driven by the streamer and the outflow.
Conclusions: These findings testify to the key role of streamers in the buildup of the disk concerning both the final mass available for planet formation and its chemical composition.

In the present study we employ three distinct physically motivated speed of sound bounds in order to construct hybrid models where the high density phase is described by the maximally stiff equation of state. In particular, we consider the bounds related to special relativity, relativistic kinetic theory and conformality. The low density hadronic phase is described by a state-of-the-art microscopic relativistic Brueckner-Hartree-Fock theory. This work aims to access the effect of the different speed of sound constraints on the relevant parameter space of the key parameters of first-order phase transitions by utilising recent astronomical data. This involves a systematic analysis that also includes two distinct schemes for the construction of hybrid models, namely the Maxwell and Gibbs methods. Finally, a relevant discussion is conducted on the possible occurrence of a thermodynamic inconsistency that is related to the stability of the high density phase over hadronic matter at large densities.

The mass ranges allowed for primordial black holes (PBHs) to constitute all of dark matter (DM) are broadly constrained. However, these constraints rely on the standard semiclassical approximation which assumes that the evaporation process is self-similar. Quantum effects such as memory burden take the evaporation process out of the semiclassical regime latest by the time the black hole loses half of its mass. What happens beyond this time is currently not known. However, theoretical evidence based on prototype models indicates that the evaporation slows down, thereby extending the lifetime of a black hole. This modifies the mass ranges constrained, in particular, by big bang nucleosynthesis (BBN) and cosmic microwave background spectral distortions. We show that previous constraints are largely relaxed when the PBH lifetime is extended, making it possible for PBHs to constitute all of DM in previously excluded mass ranges. In particular, this is the case for PBHs lighter than 10⁹ g that enter the memory burden stage before BBN and are still present today as DM.

The photon parton distribution function (PDF) of the proton is crucial for precise comparisons of LHC cross sections with theoretical predictions. However, it was previously affected by very large uncertainties of around ${\cal O}(100\%)$ or dependent upon phenomenologically inspired models. In the paper~\cite{Manohar:2016nzj}, we demonstrated how the photon PDF could be determined using the proton structure functions $F_2$ and $F_L$ measured in electron--proton scattering experiments. We provided an explicit formula for the PDF, which can be systematically improved order by order in perturbation theory. We obtained a photon PDF with errors $\lesssim 2$\% for $10^{-4} < x < 0.1$. Here, we recall the underlying idea and method used to obtain this result, as well as the progress made since then.

The Gaia mission has provided highly accurate observations that have significantly reduced the scatter in the colour-magnitude diagrams of open clusters. As a result of the improved isochrone sequence of the open cluster M67, we have created new stellar models that avoid commonly used simplifications in 1D stellar modelling, such as mass-independent core overshooting and a constant mixing length parameter. This has enabled us to deliver a precise isochrone specifically designed for M67, available for download. We follow a commonly used qualitative approach to adjust the input physics to match the well-defined colour-magnitude sequence, and we test the model-predicted masses against a known eclipsing binary system at the main sequence turnoff of the cluster. Despite using improvements in photometry and stellar physics we cannot match the masses of both binary components with the same theoretical isochrone. A $\chi ^{2}$-based isochrone fitting approach using our preferred input physics results in a cluster age of $3.95^{+ 0.16}_{- 0.15}$ Gyr.

If the gradient of a probability distribution on a landscape of vacua aligns with the variation of some fundamental parameter, the parameter may be likely to take some non-generic value. Such non-generic values can be associated to critical boundaries, where qualitative changes of the landscape properties happen, or an anthropic bound is located. Assuming the standard volume-weighted and the local probability measures, we discuss ordered landscapes which can produce several types of the aligned probability gradients. The resulting values of the gradients are defined by the "closeness" of a given vacuum to the highest- or the lowest-energy vacuum. Using these ingredients we construct a landscape scanning independently the Higgs mass and the cosmological constant (CC). The probability gradient pushes the Higgs mass to its observed value, where a structural change of the landscape takes place, while the CC is chosen anthropically.

The vector meson-baryon interaction in a coupled channel scheme is revisited within the correlation function framework. As illustrative cases to reveal the important role played by the coupled channels, we focus on the $\phi$p and $\rho^0$p systems given their complex dynamics and the presence of quasi-bound states or resonances in the vicinity of their thresholds. We show that the $\phi$p femtoscopic data provide novel information about a $N^*$ state present in the experimental region and anticipate the relevance of a future $\rho^0$p correlation function measurement in order to pin down the $S=0, Q=+1$ vector meson-baryon interaction as well as to disclose the characterizing features of the $N^*(1700)$ state.

We have applied an operator-overloading forward-mode algorithmic differentiation tool to the Monte-Carlo particle simulation toolkit Geant4. Our differentiated version of Geant4 allows computing mean pathwise derivatives of user-defined outputs of Geant4 applications with respect to user-defined inputs. This constitutes a major step towards enabling gradient-based optimization techniques in high-energy physics, as well as other application domains of Geant4. This is a preliminary report on the technical aspects of applying operator-overloading AD to Geant4, as well as a first analysis of some results obtained by our differentiated Geant4 prototype. We plan to follow up with a more refined analysis.

Context. The vast majority of close binaries containing a compact object, including the progenitors of supernovae Ia and at least a substantial fraction of all accreting black holes in the Galaxy, form through common-envelope (CE) evolution. Despite this importance, we struggle to even understand the energy budget of CE evolution. For decades, observed long-period post-CE binaries have been interpreted as evidence of additional energies contributing during CE evolution. We have recently shown that this argument is based on simplified assumptions for all long-period post-CE binaries containing massive white dwarfs (WDs). The only remaining post-CE binary star that has been claimed to require contributions from additional energy sources to understand its formation is KOI 3278.
Aims: Here, we address in detail the potential evolutionary history of KOI 3278. In particular, we investigate whether extra energy sources, such as recombination energy, are indeed required to explain its existence.
Methods: We used the 1D stellar evolution code MESA to carry out binary evolution simulations and searched for potential formation pathways for KOI 3278 that are able to explain its observed properties.
Results: We find that KOI 3278 can be explained if the WD progenitor filled its Roche lobe during a helium shell flash. In this case, the orbital period of KOI 3278 can be reproduced if the CE binding energy is calculated taking into account gravitational energy and thermodynamic internal energy. While the CE evolution that led to the formation of KOI 3278 must have been efficient - that is, most of the available orbital energy must have been used to unbind the CE - recombination energy is not required.
Conclusions: We conclude that currently not a single observed post-CE binary requires one to assume that energy sources other than gravitational and thermodynamic energy are contributing to CE evolution. KOI 3278, however, remains an intriguing post-CE binary as, unlike its siblings, understanding its existence requires highly efficient CE ejection.

Context. Galaxy clusters are the largest gravitating structures in the universe, and their mass assembly is sensitive to the underlying cosmology. Their mass function, baryon fraction, and mass distribution have been used to infer cosmological parameters despite the presence of systematics. However, the complexity of the scaling relations among galaxy cluster properties has never been fully exploited, limiting their potential as a cosmological probe.
Aims: We propose the first machine learning (ML) method using galaxy cluster properties from hydrodynamical simulations in different cosmologies to predict cosmological parameters combining a series of canonical cluster observables, such as gas mass, gas bolometric luminosity, gas temperature, stellar mass, cluster radius, total mass, and velocity dispersion at different redshifts.
Methods: The ML model was trained on mock "measurements" of these observable quantities from Magneticum multi-cosmology simulations to derive unbiased constraints on a set of cosmological parameters. These include the mass density parameter, Ω_m, the power spectrum normalization, σ₈, the baryonic density parameter, Ω_b, and the reduced Hubble constant, h₀.
Results: We tested the ML model on catalogs of a few hundred clusters taken, in turn, from each simulation and found that the ML model can correctly predict the cosmology from where they have been picked. The cumulative accuracy depends on the cosmology, ranging from 21% to 75%. We demonstrate that this is sufficient to derive unbiased constraints on the main cosmological parameters with errors on the order of ~14% for Ω_m, ~8% for σ₈, ~6% for Ω_b, and ~3% for h₀.
Conclusions: This proof-of-concept analysis, though based on a limited variety of multi-cosmology simulations, shows that ML can efficiently map the correlations in the multidimensional space of the observed quantities to the cosmological parameter space and narrow down the probability that a given sample belongs to a given cosmological parameter combination. More large-volume, mid-resolution, multi-cosmology hydro-simulations need to be produced to expand the applicability to a wider cosmological parameter range. However, this first test is exceptionally promising, as it shows that these ML tools can be applied to cluster samples from multiwavelength observations from surveys such as Rubin/LSST, CSST, Euclid, and Roman in optical and near-infrared bands, and eROSITA in X-rays, to the constrain cosmology and effect of baryonic feedback.

We present JWST NIRCam (F356W and F444W filters) and MIRI (F770W) images and NIRSpec Integral Field Unit (IFU) spectroscopy of the young Galactic supernova remnant Cassiopeia A (Cas A) to probe the physical conditions for molecular CO formation and destruction in supernova ejecta. We obtained the data as part of a JWST survey of Cas A. The NIRCam and MIRI images map the spatial distributions of synchrotron radiation, Ar-rich ejecta, and CO on both large and small scales, revealing remarkably complex structures. The CO emission is stronger at the outer layers than the Ar ejecta, which indicates the re-formation of CO molecules behind the reverse shock. NIRSpec-IFU spectra (3–5.5 μm) were obtained toward two representative knots in the NE and S fields that show very different nucleosynthesis characteristics. Both regions are dominated by the bright fundamental rovibrational band of CO in the two R and P branches, with strong [Ar VI] and relatively weaker, variable strength ejecta lines of [Si IX], [Ca IV], [Ca V], and [Mg IV]. The NIRSpec-IFU data resolve individual ejecta knots and filaments spatially and in velocity space. The fundamental CO band in the JWST spectra reveals unique shapes of CO, showing a few tens of sinusoidal patterns of rovibrational lines with pseudocontinuum underneath, which is attributed to the high-velocity widths of CO lines. Our results with LTE modeling of CO emission indicate a temperature of ∼1080 K and provide unique insight into the correlations between dust, molecules, and highly ionized ejecta in supernovae and have strong ramifications for modeling dust formation that is led by CO cooling in the early Universe.

Past studies have long emphasised the key role played by galactic stellar bars in the context of disc secular evolution, via the redistribution of gas and stars, the triggering of star formation, and the formation of prominent structures such as rings and central mass concentrations. However, the exact physical processes acting on those structures, as well as the timescales associated with the building and consumption of central gas reservoirs are still not well understood. We are building a suite of hydro-dynamical RAMSES simulations of isolated, low-redshift galaxies that mimic the properties of the PHANGS sample. The initial conditions of the models reproduce the observed stellar mass, disc scale length, or gas fraction, and this paper presents a first subset of these models. Most of our simulated galaxies develop a prominent bar structure, which itself triggers central gas fuelling and the building of an over-density with a typical scale of 100−1000 pc. We confirm that if the host galaxy features an ellipsoidal component, the formation of the bar and gas fuelling are delayed. We show that most of our simulations follow a common time evolution, when accounting for mass scaling and the bar formation time. In our simulations, the stellar mass of 10¹⁰ M_⊙ seems to mark a change in the phases describing the time evolution of the bar and its impact on the interstellar medium. In massive discs (M_⋆ ≥ 10¹⁰ M_⊙), we observe the formation of a central gas reservoir with star formation mostly occurring within a restricted starburst region, leading to a gas depletion phase. Lower-mass systems (M_⋆ < 10¹⁰ M_⊙) do not exhibit such a depletion phase, and show a more homogeneous spread of star-forming regions along the bar structure, and do not appear to host inner bar-driven discs or rings. Our results seem to be supported by observations, and we briefly discuss how this new suite of simulations can help our understanding of the secular evolution of main sequence disc galaxies.

We present the serendipitous discovery of (1) a large double radio relic associated with the galaxy cluster PSZ2 G277.93 + 12.34 and (2) a new odd radio circle, ORC J1027-4422, both found in the same deep MeerKAT 1.3 GHz wide-band radio continuum image. The angular separation of the two arc-shaped cluster relics is ~16 arcmin or ~2.6 Mpc for a cluster redshift of z ≈ 0.158. The thin southern relic, which shows several ridges/shocks including one possibly moving inwards, has a linear extent of ~1.64 Mpc. In contrast, the northern relic is about twice as wide, twice as bright, but only has a largest linear size of ~0.66 Mpc. Complementary SRG/eROSITA X-ray images reveal extended emission from hot intracluster gas between the two relics and around the narrow-angle tail (NAT) radio galaxy PMN J1033-4335 (z ≈ 0.153) located just east of the northern relic. The radio morphologies of the NAT galaxy and the northern relic, which are also detected with the Australian Square Kilometer Array Pathfinder (ASKAP) at 888 MHz, suggest both are moving in the same outward direction. The discovery of ORC J1027-4422 in a different part of the same MeerKAT image makes it the fourth known single ORC. It has a diameter of ~90 arcsec corresponding to 400 kpc at a tentative redshift of z ≈ 0.3 and remains undetected in X-ray emission. Supported by simulations, we discuss similarities between outward moving galaxy and cluster merger shocks as the formation mechanisms for ORCs and radio relics, respectively.

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.

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.

Context. Magnetic fields generated in the early Universe undergo turbulent decay during the radiation-dominated era. The decay is governed by a decay exponent and a decay time. It has been argued that the latter is prolonged by magnetic reconnection, which depends on the microphysical resistivity and viscosity. Turbulence, on the other hand, is not usually expected to be sensitive to microphysical dissipation, which affects only very small scales.
Aims: We want to test and quantify the reconnection hypothesis in decaying hydromagnetic turbulence.
Methods: We performed high-resolution numerical simulations with zero net magnetic helicity using the PENCIL CODE with up to 2048³ mesh points and relate the decay time to the Alfvén time for different resistivities and viscosities.
Results: The decay time is found to be longer than the Alfvén time by a factor that increases with increasing Lundquist number to the 1/4 power. The decay exponent is as expected from the conservation of the Hosking integral, but a timescale dependence on resistivity is unusual for developed turbulence and not found for hydrodynamic turbulence. In two dimensions, the Lundquist number dependence is shown to be leveling off above values of ≈25 000, independently of the value of the viscosity.
Conclusions: Our numerical results suggest that resistivity effects have been overestimated in earlier work. Instead of reconnection, it may be the magnetic helicity density in smaller patches that is responsible for the resistively slow decay. The leveling off at large Lundquist number cannot currently be confirmed in three dimensions.

The nature of dark matter (DM) and its interaction with the Standard Model (SM) is one of the biggest open questions in physics nowadays. The vast majority of theoretically motivated ultralight-DM (ULDM) models predict that ULDM couples dominantly to the SM strong/nuclear sector. This coupling leads to oscillations of nuclear parameters that are detectable by comparing clocks with different sensitivities to these nature's constants. Vibrational transitions of molecular clocks are more sensitive to a change in the nuclear parameters than the electronic transitions of atomic clocks. Here, we propose the iodine molecular ion, I₂⁺, as a sensitive detector for such a class of ULDM models. The iodine's dense spectrum allows us to match its transition frequency to that of an optical atomic clock (Ca⁺) and perform correlation spectroscopy between the two clock species. With this technique, we project a few-orders-of-magnitude improvement over the most sensitive clock comparisons performed to date. We also briefly consider the robustness of the corresponding "Earth-bound" under modifications of the Z_N -QCD axion model.

For non-relativistic thermal dark matter, close-to-threshold effects largely dominate the evolution of the number density for most of the times after thermal freeze-out, and hence affect the cosmological relic density. A precise evaluation of the relevant interaction rates in a thermal medium representing the early universe includes accounting for the relative motion of the dark matter particles and the thermal medium. We consider a model of dark fermions interacting with a plasma of dark gauge bosons, which is equivalent to thermal QED. The temperature is taken to be smaller than the dark fermion mass and the inverse of the typical size of the dark fermion-antifermion bound states, which allows for the use of non-relativistic effective field theories. For the annihilation cross section, bound-state formation cross section, bound-state dissociation width and bound-state transition width of dark matter fermion-antifermion pairs, we compute the leading recoil effects in the reference frame of both the plasma and the center-of-mass of the fermion-antifermion pair. We explicitly verify the Lorentz transformations among these quantities. We evaluate the impact of the recoil corrections on the dark matter energy density. Our results can be directly applied to account for the relative motion of quarkonia in the quark-gluon plasma formed in heavy-ion collisions. They may be also used to precisely assess thermal effects in atomic clocks based on atomic transitions; the present work provides a first field theory derivation of time dilation for these processes in vacuum and in a medium.

We use functional methods to match the two-Higgs-doublet model with heavy scalars in the nondecoupling regime to the appropriate nonlinear effective field theory, which takes the form of an electroweak chiral Lagrangian (HEFT). The effective Lagrangian is derived to leading order in the chiral counting. This includes the loop induced h →γ γ and h →Z γ local terms, which enter at the same chiral order as their counterparts in the Standard Model. An algorithm is presented that allows us to compute the coefficient functions to all orders in h . Some of the all-orders results are given in closed form. The parameter regimes for decoupling, nondecoupling, and alignment scenarios in the effective field theory context and some phenomenological implications are briefly discussed.

High-frequency gravitational waves (HFGWs) are predicted in various exotic scenarios involving both cosmological and astrophysical sources. These elusive signals have recently sparked the interest of a diverse community of researchers, due to the possibility of HFGW detection in the laboratory through graviton-photon conversion in strong magnetic fields. Notable examples include the redesign of the resonant cavities currently under development to detect the cosmic axion. In this work, we derive the sensitivities of some existing and planned resonant cavities to detect a HFGW background. As a concrete scenario, we consider the collective signals that originate from the merging of compact objects, such as two primordial black holes (PBHs) in the asteroid mass window. Our findings improve over existing work by explicitly discussing and quantifying the loss in the experimental reach due to the actual coherence of the source. We elucidate on the approach we adopt in relation with recent literature on the topic. Most notably, we give a recipe for the estimate of the stochastic background that focuses on the presence of the signal in the cavity at all times and showing that, in the relevant PBH mass region, the signal is dominated by coherent binary mergers.

Flavor-dependent neutrino transport is described by a well-known kinetic equation for occupation-number matrices in flavor space. However, in the context of fast flavor conversion, we identify an unforeseen predicament: the pivotal self-induced exponential growth of small inhomogeneities strongly violates conservation of neutrino-neutrino refractive energy. We prove that it is traded with the huge reservoir of neutrino kinetic energy through gradients of neutrino flavor coherence (the off-diagonal piece of the flavor density matrix) and derive the missing gradient terms. The usual equations remain sufficient to describe flavor evolution, at the cost of renouncing energy conservation, which cannot play any role in explaining the numerically observed final state.

We propose a simple model that can alleviate the H₀ tension while remaining consistent with big bang nucleosynthesis (BBN). It is based on a dark sector described by a standard Lagrangian featuring a S U (N ) gauge symmetry with N ≥3 and a massive scalar field with a quartic coupling. The scalar acts as a dark Higgs leading to spontaneous symmetry breaking S U (N )→S U (N -1 ) via a first-order phase transition à la Coleman-Weinberg. This setup naturally realizes previously proposed scenarios featuring strongly interacting dark radiation (SIDR) with a mass threshold within hot new early dark energy. For a wide range of reasonable model parameters, the phase transition occurs between the BBN and recombination epochs and releases a sufficient amount of latent heat such that the model easily respects bounds on extra radiation during BBN while featuring a sufficient SIDR density around recombination for increasing the value of H₀ inferred from the cosmic microwave background. Our model can be summarized as a natural mechanism providing two successive increases in the effective number of relativistic degrees of freedom after BBN but before recombination Δ N_BBN→Δ N_NEDE→Δ N_IR alleviating the Hubble tension. The first step is related to the phase transition, and the second is related to the dark Higgs becoming nonrelativistic. This setup predicts further signatures, including a stochastic gravitational wave background and features in the matter power spectrum that can be searched for with future pulsar timing and Lyman-α forest measurements.

We examine the physical significance of torsion co-cycles in the cohomology of a projective Calabi-Yau three-fold for the (2,2) superconformal field theory (SCFT) associated to the non-linear sigma model with such a manifold as a target space. There are two independent torsion subgroups in the cohomology. While one is associated to an orbifold construction of the SCFT, the other encodes the possibility of turning on a topologically non-trivial flat gerbe for the NS-NS B-field. Inclusion of these data enriches mirror symmetry by providing a refinement of the familiar structures and points to a generalization of the duality symmetry, where the topology of the flat gerbe enters on the same footing as the topology of the underlying manifold.

In this paper we present a novel method to extract information on hadron-hadron interactions using for the first time femtoscopic data to constrain the low-energy constants of a QCD effective Lagrangian. This method offers a new way to investigate the nonperturbative regime of QCD in sectors where scattering experiments are not feasible, such as the multistrange and charm ones. As an example of its application, we use the very precise K^-Λ correlation function data, recently measured in p p collisions at LHC, to constrain the strangeness S =-2 meson-baryon interaction. The model obtained delivers new insights on the molecular nature of the Ξ (1620 ) and Ξ (1690 ) states.

We present a novel method for identifying transients suitable for both strong signal-dominated and background-dominated objects. By employing the unsupervised machine learning algorithm known as expectation maximization, we achieve computing time reductions of over 10⁴ on a single CPU compared to conventional brute-force methods. Furthermore, this approach can be readily extended to analyze multiple flares. We illustrate the algorithm's application by fitting the IceCube neutrino flare of TXS 0506+056.

In this work we compare the predictions for the scattering length and effective range of the channels <inline-formula id="IEq3"><mml:math><mml:mrow><mml:msup><mml:mi>K</mml:mi><mml:mn>0</mml:mn></mml:msup><mml:msup><mml:mi mathvariant="normal">Σ</mml:mi><mml:mo>+</mml:mo></mml:msup><mml:mo>,</mml:mo><mml:msup><mml:mi>K</mml:mi><mml:mo>+</mml:mo></mml:msup><mml:msup><mml:mi mathvariant="normal">Σ</mml:mi><mml:mn>0</mml:mn></mml:msup><mml:mo>,</mml:mo><mml:msup><mml:mi>K</mml:mi><mml:mo>+</mml:mo></mml:msup><mml:mi mathvariant="normal">Λ</mml:mi></mml:mrow></mml:math></inline-formula> and <inline-formula id="IEq4"><mml:math><mml:mrow><mml:mi>η</mml:mi><mml:mi>p</mml:mi></mml:mrow></mml:math></inline-formula>, assuming the <inline-formula id="IEq5"><mml:math><mml:mrow><mml:msup><mml:mi>N</mml:mi><mml:mo>∗</mml:mo></mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mn>1535</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:math></inline-formula> state as a molecular state of these channels, or an original genuine state, made for instance from three quarks. Looking at very different scenarios, what we conclude is that the predictions of these two pictures are drastically different, to the point that we advise the measurement of these magnitudes, accessible for instance by measuring correlation functions, in order to gain much valuable information concerning the nature of this state.

In this contribution, we aim to summarise the efforts of the Italian SKA scientific community in conducting surveys of star-forming regions within our Galaxy, in the development of astrochemical research on protostellar envelopes and disks, and in studying the planet formation process itself. The objective is dual: Firstly, to investigate the accumulation and development of dust throughout the formation of planets, and secondly, to chemically examine protoplanetary disks and protostellar envelopes by studying heavy molecules, such as chains and rings containing over seven carbon atoms, which exhibit significantly reduced strength at millimeter wavelengths.

The rotational properties of the transfermium nuclei are investigated in the full deformation space by implementing a shell-model-like approach in the cranking covariant density functional theory on a three-dimensional lattice, where the pairing correlations, deformations, and moments of inertia are treated in a microscopic and self-consistent way. The kinematic and dynamic moments of inertia of the rotational bands observed in the transfermium nuclei No 252 , No 254 , Rf 254 , and Rf 256 are well reproduced without any adjustable parameters using a well-determined universal density functional. It is found for the first time that the emergence of the octupole deformation should be responsible for the significantly different rotational behavior observed in No 252 and No 254 . The present results provide a microscopic solution to the long-standing puzzle on the rotational behavior in No isotopes, and highlight the risk of investigating only the hexacontetrapole (β₆₀) deformation effects in rotating transfermium nuclei without considering the octupole deformation.

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

Dark photons have emerged as promising candidates for dark matter, and their search is a top priority in particle physics, astrophysics, and cosmology. We report the first use of a tunable niobium superconducting radio-frequency cavity for a scan search of dark photon dark matter with innovative data analysis techniques. We mechanically adjusted the resonant frequency of a cavity submerged in liquid helium at a temperature of 2 K, and scanned the dark photon mass over a frequency range of 1.37 MHz centered at 1.3 GHz. Our study leveraged the superconducting radio-frequency cavity's remarkably high quality factors of approximately 10¹⁰, resulting in the most stringent constraints to date on a substantial portion of the exclusion parameter space on the kinetic mixing coefficient ε between dark photons and electromagnetic photons, yielding a value of ε <2.2 ×10^-16.

The study of next-to-leading-power (NLP) corrections in soft emissions continues to attract interest both in quantum chromodynamics (QCD) and in quantum electrodynamics (QED). Soft-photon spectra in particular provide a clean case-study for the experimental verification of the Low-Burnett-Kroll (LBK) theorem. In this paper we study the consistency of the LBK theorem in the context of an ambiguity arising from momentum-conservation constraints in the computation of nonradiative amplitudes. We clarify that this ambiguity leads to various possible formulations of the LBK theorem, which are all equivalent up to power-suppressed effects (i.e., beyond the formal accuracy of the LBK theorem). We also propose a new formulation of the LBK theorem with a modified shifted kinematics which facilitates the numerical computation of nonradiative amplitudes with publicly available tools. Furthermore, we present numerical results for soft-photon spectra in the associated production of a muon pair with a photon, both in e⁺e^- annihilation and proton-proton collisions.

Biomolecular condensates help organize the cell cytoplasm and nucleoplasm into spatial compartments with different chemical compositions. A key feature of such compositional patterning is the local enrichment of enzymatically active biomolecules which, after transient binding via molecular interactions, catalyze reactions among their substrates. Thereby, biomolecular condensates provide a spatial template for nonuniform concentration profiles of substrates. In turn, the concentration profiles of substrates, and their molecular interactions with enzymes, drive enzyme fluxes which can enable novel nonequilibrium dynamics. To analyze this generic class of systems, with a current focus on self-propelled droplet motion, we here develop a self-consistent sharp interface theory. In our theory, we diverge from the usual bottom-up approach, which involves calculating the dynamics of concentration profiles based on a given chemical potential gradient. Instead, reminiscent of control theory, we take the reverse approach by deriving the chemical potential profile and enzyme fluxes required to maintain a desired condensate form and dynamics. The chemical potential profile and currents of enzymes come with a corresponding power dissipation rate, which allows us to derive a thermodynamic consistency criterion for the passive part of the system (here, reciprocal enzyme-enzyme interactions). As a first-use case of our theory, we study the role of reciprocal interactions, where the transport of substrates due to reactions and diffusion is, in part, compensated by redistribution due to molecular interactions. More generally, our theory applies to mass-conserved active matter systems with moving phase boundaries.

We report the first experimental results obtained with the new haloscope of the QUAX experiment located at Laboratori Nazionali di Frascati of INFN (LNF). The haloscope is composed of a OFHC Cu resonant cavity cooled down to about 30 mK and immersed in a magnetic field of 8 T. The cavity frequency was varied in a 6 MHz range between 8.831496 and 8.83803 GHz. This corresponds to a previously unprobed mass range between 36.52413 and 36.5511 μ eV . We don't observe any excess in the power spectrum and set limits on the axion-photon coupling in this mass range down to g_{a γ γ}<0.882 ×10^-13 GeV^-1 with the confidence level set to 90%.

We introduce our novel Bayesian parton density determination code, PartonDensity.jl. The motivation for this new code, the framework, and its validation are described. As we show, PartonDensity.jl provides both a flexible environment for the determination of parton densities and a wealth of information concerning the knowledge update provided by the analyzed dataset.

Recent observations of high-energy neutrinos from active galactic nuclei (AGN), NGC 1068 and TXS 0506 +056 , suggest that cosmic rays (CRs) are accelerated in the vicinity of the central supermassive black hole and high-energy protons and electrons can cool efficiently via interactions with ambient photons and gas. The dark matter density may be significantly enhanced near the black hole, and CRs could lose energies predominantly due to scatterings with the ambient dark matter particles. We propose CR cooling in AGN as a new probe of dark matter-proton and dark matter-electron scatterings. Under plausible astrophysical assumptions, our constraints on sub-GeV dark matter can be the strongest derived to date. Some of the parameter space favored by thermal light dark matter models might already be probed with current multimessenger observations of AGN.

Context. Self-interacting dark matter (SIDM) can tackle or alleviate small-scale issues within the cosmological standard model ΛCDM, and diverse flavours of SIDM can produce unique astrophysical predictions, resulting in different possible signatures which can be used to test these models with dedicated observations of galaxy clusters.
Aims: This work aims to assess the impact of dark matter self-interactions on the properties of galaxy clusters. In particular, the goal is to study the angular dependence of the cross section by testing rare (large angle scattering) and frequent (small angle scattering) SIDM models with velocity-dependent cross sections.
Methods: We re-simulated six galaxy cluster zoom-in initial conditions with a dark matter-only run and with full-physics set-up simulations that include a self-consistent treatment of baryon physics. We tested the dark matter-only setup and the full physics setup with the collisionless cold dark matter, rare self-interacting dark matter, and frequent self-interacting dark matter models. We then studied their matter density profiles as well as their subhalo population.
Results: Our dark matter-only SIDM simulations agree with theoretical models, and when baryons are included in simulations, our SIDM models substantially increase the central density of galaxy cluster cores compared to full-physics simulations using collisionless dark matter. SIDM subhalo suppression in full-physics simulations is milder compared to the one found in the dark matter-only simulations because of the cuspier baryonic potential that prevents subhalo disruption. Moreover, SIDM with small-angle scattering significantly suppresses a larger number of subhaloes compared to large-angle scattering SIDM models. Additionally, SIDM models generate a broader range of subhalo concentration values, including a tail of more diffuse subhaloes in the outskirts of galaxy clusters and a population of more compact subhaloes in the cluster cores.

It has been recently suggested that the strong Emergence Proposal is realized in M-theory limits by integrating out all light towers of states with a typical mass scale not larger than the species scale, i.e. the eleventh dimensional Planck mass. Within the BPS sector, these are transverse M2- and M5-branes, that can be wrapped and particle-like, carrying Kaluza-Klein momentum along the compact directions. We provide additional evidence for this picture by revisiting and investigating further the computation of R⁴-interactions in M-theory à la Green-Gutperle-Vanhove. A central aspect is a novel UV-regularization of Schwinger-like integrals, whose actual meaning and power we clarify by first applying it to string perturbation theory. We consider then toroidal compactifications of M-theory and provide evidence that integrating out all light towers of states via Schwinger-like integrals thus regularized yields the complete result for R⁴-interactions. In particular, this includes terms that are tree-level, one-loop and space-time instanton corrections from the weakly coupled point of view. Finally, we comment on the conceptual difference of our approach to earlier closely related work by Kiritsis-Pioline and Obers-Pioline, leading to a correspondence between two types of constrained Eisenstein series.

We present an extension to a Sunyaev-Zel'dovich Effect (SZE) selected cluster catalogue based on observations from the South Pole Telescope (SPT); this catalogue extends to lower signal to noise than the previous SPT-SZ catalogue and therefore includes lower mass clusters. Optically derived redshifts, centres, richnesses, and morphological parameters together with catalogue contamination and completeness statistics are extracted using the multicomponent matched filter (MCMF) algorithm applied to the S/N > 4 SPT-SZ candidate list and the Dark Energy Survey (DES) photometric galaxy catalogue. The main catalogue contains 811 sources above S/N = 4, has 91 per cent purity, and is 95 per cent complete with respect to the original SZE selection. It contains in total 50 per cent more clusters and twice as many clusters above z = 0.8 in comparison to the original SPT-SZ sample. The MCMF algorithm allows us to define subsamples of the desired purity with traceable impact on catalogue completeness. As an example, we provide two subsamples with S/N > 4.25 and S/N > 4.5 for which the sample contamination and cleaning-induced incompleteness are both as low as the expected Poisson noise for samples of their size. The subsample with S/N > 4.5 has 98 per cent purity and 96 per cent completeness and is part of our new combined SPT cluster and DES weak-lensing cosmological analysis. We measure the number of false detections in the SPT-SZ candidate list as function of S/N, finding that it follows that expected from assuming Gaussian noise, but with a lower amplitude compared to previous estimates from simulations.

We investigate the accuracy of the perturbative galaxy bias expansion in view of the forthcoming analysis of the Euclid spectroscopic galaxy samples. We compare the performance of a Eulerian galaxy bias expansion using state-of-the-art prescriptions from the effective field theory of large-scale structure (EFTofLSS) with a hybrid approach based on Lagrangian perturbation theory and high-resolution simulations. These models are benchmarked against comoving snapshots of the flagship I N-body simulation at z = (0.9, 1.2, 1.5, 1.8), which have been populated with Hα galaxies leading to catalogues of millions of objects within a volume of about 58 h⁻³ Gpc³. Our analysis suggests that both models can be used to provide a robust inference of the parameters (h, ω_c) in the redshift range under consideration, with comparable constraining power. We additionally determine the range of validity of the EFTofLSS model in terms of scale cuts and model degrees of freedom. From these tests, it emerges that the standard third-order Eulerian bias expansion - which includes local and non-local bias parameters, a matter counter term, and a correction to the shot-noise contribution - can accurately describe the full shape of the real-space galaxy power spectrum up to the maximum wavenumber of k_max = 0.45 h Mpc⁻¹, and with a measurement precision of well below the percentage level. Fixing either of the tidal bias parameters to physically motivated relations still leads to unbiased cosmological constraints, and helps in reducing the severity of projection effects due to the large dimensionality of the model. We finally show how we repeated our analysis assuming a volume that matches the expected footprint of Euclid, but without considering observational effects, such as purity and completeness, showing that we can get constraints on the combination (h, ω_c) that are consistent with the fiducial values to better than the 68% confidence interval over this range of scales and redshifts.

Context. The production of neutron-rich elements at neutron densities intermediate to those of the s- and r-processes, the so-called i-process, has been identified as possibly being responsible for the observed abundance pattern found in certain carbon-enhanced metal-poor (CEMP) stars. The production site may be low-metallicity stars on the asymptotic giant branch (AGB) where the physical processes during the thermal pulses are not well known.
Aims: We investigate the impact of overshoot from various convective boundaries during the AGB phase on proton ingestion events (PIEs) and the neutron densities as a necessary precondition for the i-process as well as on the structure and continued evolution of the models.
Methods: We therefore analyzed models of a 1.2 M_⊙, Z = 5 × 10⁻⁵ star. A fiducial model without overshoot on the AGB (overshoot was applied during the pre-AGB evolution) serves as a reference. The same model was then run with various overshoot values and the resulting models were compared to one another. Light element nucleosynthesis is also discussed. Additionally, we introduce a new timescale argument to predict PIE occurrence to discriminate between a physical and a numerical reason for a nonoccurrence. A comparison to observations as well as previous studies was conducted before finally presenting the most promising choice of overshoot parameters for the occurrence of the i-process in low-mass, low-metallicity models.
Results: The fiducial model reveals high neutron densities and a persistent split of the pulse-driven convection zone (PDCZ). Overshoot from the PDCZ results in either temporary or permanent remerging of the split PDCZ, influencing the star's structure and evolution. While both overshoot and non-overshoot models exhibit PIEs generating neutron densities suitable for the i-process, they lead to varied C/O and N/O ratios and notable Li enhancements. Comparison with previous studies and observations of CEMP-r/s stars suggests that while surface enhancements in our models may be exaggerated, abundance ratios align well. Though, for high values of overshoot from the PDCZ the agreement becomes worse.

Context. Recent observations of galaxy mergers inside galaxy cluster environments, such as NGC 5291 in the vicinity of Abell 3574, report high star formation rates in the ejected tidal tails, which point towards currently developing tidal dwarf galaxies. This prompts the intriguing question whether these newly formed stellar structures could get stripped from the galaxy potential by the cluster and thus populate it with dwarf galaxies.
Aims: We verify whether environmental stripping of tidal dwarf galaxies from galaxy mergers inside galaxy cluster environments is a possible evolutionary channel to populate a galaxy cluster with low-mass and low surface brightness galaxies.
Methods: We performed three high-resolution hydrodynamical simulations of mergers between spiral galaxies in a cluster environment, implementing a stellar mass ratio of 2:1 with M₂₀₀ = 9.5 × 10¹¹ M_⊙ for the more massive galaxy. Between the three different simulations, we varied the initial orbit of the infalling galaxies with respect to the cluster center.
Results: We demonstrate that cluster environments are capable of stripping tidal dwarf galaxies from the host potential independently of the infall orbit of the merging galaxy pair, without instantly destroying the tidal dwarfs. Starting to evolve separately from their progenitor, these newly formed dwarf galaxies reach total masses of M_tot ≈ 10^{7 − 9} M_⊙ within the limits of our resolution. In the three tested orbit scenarios, we find three, seven, and eight tidal dwarf galaxies per merger, respectively, which survive longer than 1 Gyr after the merger event. Exposed to ram pressure, these gas dominated dwarf galaxies exhibit high star formation rates while also losing gas to the environment. Experiencing a strong headwind due to their motion through the intracluster medium, they quickly lose momentum and start spiraling towards the cluster center, reaching distances on the order of 1 Mpc from their progenitor. About 4 Gyr after the merger event, we still find three and four intact dwarf galaxies in two of the tested scenarios, respectively. The other stripped tidal dwarf galaxies either evaporate in the hostile cluster environment due to their low initial mass, or are disrupted as soon as they reach the cluster center.
Conclusions: The dwarf production rate due to galaxy mergers is elevated when the interaction with a cluster environment is taken into account. Comparing their contribution to the observed galaxy mass function in clusters, our results indicate that ∼30% of dwarf galaxies in clusters could have been formed by stripping from galaxy mergers.

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

We present detailed forecasts for the constraints on the characteristics of primordial magnetic fields (PMFs) generated prior to recombination that will be obtained with the LiteBIRD satellite. The constraints are driven by some of the main physical effects of PMFs on the CMB anisotropies: the gravitational effects of magnetically-induced perturbations; the effects on the thermal and ionization history of the Universe; the Faraday rotation imprint on the CMB polarization spectra; and the non-Gaussianities induced in polarization anisotropies. LiteBIRD represents a sensitive probe for PMFs. We explore different levels of complexity, for LiteBIRD data and PMF configurations, accounting for possible degeneracies with primordial gravitational waves from inflation. By exploiting all the physical effects, LiteBIRD will be able to improve the current limit on PMFs at intermediate and large scales coming from Planck. In particular, thanks to its accurate B-mode polarization measurement, LiteBIRD will improve the constraints on infrared configurations for the gravitational effect, giving B ^{n _B=-2.9} _{1 Mpc}< 0.8 nG at 95% C.L., potentially opening the possibility to detect nanogauss fields with high significance. We also observe a significant improvement in the limits when marginalized over the spectral index, B ^{n _Bmarg} _{1 Mpc}< 2.2 nG at 95 % C.L. From the thermal history effect, which relies mainly on E-mode polarization data, we obtain a significant improvement for all PMF configurations, with the marginalized case, √⟨B ²⟩^marg<0.50 nG at 95 % C.L. Faraday rotation constraints will take advantage of the wide frequency coverage of LiteBIRD and the high sensitivity in B modes, improving the limits by orders of magnitude with respect to current results, B ^{n _B=-2.9} _{1 Mpc} < 3.2 nG at 95 % C.L. Finally, non-Gaussianities of the B-mode polarization can probe PMFs at the level of 1 nG, again significantly improving the current bounds from Planck. Altogether our forecasts represent a broad collection of complementary probes based on widely tested methodologies, providing conservative limits on PMF characteristics that will be achieved with the LiteBIRD satellite.

NGC 1068 is a nearby, widely studied Seyfert II galaxy presenting radio, infrared, X-ray, and γ-ray emission, along with strong evidence for high-energy neutrino emission. Recently, the evidence for neutrino emission was explained in a multimessenger model, whereby the neutrinos originate from the corona of the active galactic nucleus. In this environment, γ-rays are strongly absorbed, so that an additional contribution is necessary, for instance, from the circumnuclear starburst ring. In this work, we discuss whether the radio jet can be an alternative source of the γ-rays between about 0.1 and 100 GeV, as observed by Fermi-LAT. In particular, we include both leptonic and hadronic processes, namely, accounting for inverse Compton emission and signatures from pp as well as pγ interactions. In order to constrain our calculations, we used VLBA and ALMA observations of the radio knot structures, which are spatially resolved at different distances from the supermassive black hole. Our results show that the best leptonic scenario for the prediction of the Fermi-LAT data is provided by the radio knot closest to the central engine. For that to be the case, a magnetic field strength of ∼1 mG is needed as well as a strong spectral softening of the relativistic electron distribution at (1 − 10) GeV. However, we show that neither such a weak magnetic field strength, nor such a strong softening is expected for that knot. A possible explanation for the ∼10 GeV γ-rays could potentially be provided by hadronic pion production in case of a gas density ≳10⁴ cm⁻³. Nonetheless, this process is not found to contribute significantly to the low-energy end of the Fermi-LAT range. We conclude that the emission sites in the jet are not sufficient to explain the γ-rays across the whole Fermi-LAT energy band.

Context. Signposts of early planet formation are ubiquitous in substructured young discs. Dense, hot, and high-pressure regions that formed during the gravitational collapse process, integral to star formation, facilitate dynamical mixing of dust within the protostellar disc. This provides an incentive to constrain the role of gas and dust interaction and resolve potential zones of dust concentration during star and disc formation stages.
Aims: We explore whether the thermal and dynamical conditions that developed during protostellar disc formation can generate gas flows that efficiently mix and transport the well-coupled gas and dust components.
Methods: We simulated the collapse of dusty molecular cloud cores with the hydrodynamics code PLUTO augmented with radiation transport and self-gravity. We used a two-dimensional axisymmetric geometry and followed the azimuthal component of the velocity. The dust was treated as Lagrangian particles that are subject to drag from the gas, whose motion is computed on a Eulerian grid. We considered 1, 10, and 100 µm-sized neutral, spherical dust grains. Importantly, the equation of state accurately includes molecular hydrogen dissociation. We focus on molecular cloud core masses of 1 and 3 M_⊙ and explore the effects of different initial rotation rates and cloud core sizes.
Results: Our study underlines mechanisms for the early transport of dust from the inner hot disc regions via the occurrence of two transient gas motions, namely meridional flow and outflow. The vortical flow fosters dynamical mixing and retention of dust, while the thermal pressure driven outflow replenishes dust in the outer disc. Notably, these phenomena occur regardless of the initial cloud core mass, size, and rotation rate.
Conclusions: Young dynamical precursors to planet-forming discs exhibit regions with complex hydrodynamical gas features and high-temperature structures. These can play a crucial role in concentrating dust for subsequent growth into protoplanets. Dust transport, especially, from sub-au scales surrounding the protostar to the outer relatively cooler parts, offers an efficient pathway for thermal reprocessing during pre-stellar core collapse.

We suggest that filaments in star-forming regions undergo frequent mergers. As stellar nurseries, filaments play a vital role in understanding star formation and mergers could pave the way for understanding the formation of more complex filamentary systems, such as networks and hubs. We compare the physical properties derived from hydrodynamic RAMSES simulations of merging filaments to those obtained from ALMA observations towards the LDN 1641-North region in Orion. We find similarities in the distributions of line-mass, column density, and velocity dispersion. Such common features support the hypothesis of filament mergers shaping the structure of the interstellar medium.

Ultrahot Jupiters are a type of gaseous exoplanet that orbit extremely close to their host star, resulting in significantly high equilibrium temperatures. In recent years, high-resolution emission spectroscopy has been broadly employed in observing the atmospheres of ultrahot Jupiters. We used the CARMENES spectrograph to observe the high-resolution spectra of the dayside hemisphere of MASCARA-1b in both visible and near-infrared. Through cross-correlation analysis, we detected signals of Fe I and Ti I. Based on these detections, we conducted an atmospheric retrieval and discovered the presence of a strong inversion layer in the planet's atmosphere. The retrieved Ti and Fe abundances are broadly consistent with solar abundances. In particular, we obtained a relative abundance of [Ti/Fe] as −1.0 ± 0.8 under the free retrieval and −0.4_−0.8^+0.5 under the chemical equilibrium retrieval, suggesting the absence of significant titanium depletion on this planet. Furthermore, we considered the influence of planetary rotation on spectral line profiles. The resulting equatorial rotation speed was determined to be 4.4_−2.0^+1.6 km s⁻¹, which agrees with the rotation speed induced by tidal locking.

Context. This is the second paper in a series presenting the results from a 500 h⁻¹Mpc large constrained simulation of the local Universe (SLOW). The initial conditions for this cosmological hydro-dynamical simulation are based on peculiar velocities derived from the CosmicFlows-2 catalog. The simulation follows cooling, star formation, and the evolution of super-massive black holes. This allows one to directly predict observable properties of the intracluster medium (ICM) within galaxy clusters, including X-ray luminosity, temperatures, and the Compton-y signal.
Aims: Comparing the properties of observed galaxy clusters within the local Universe with the properties of their simulated counterparts enables us to assess the effectiveness of the initial condition constraints in accurately replicating the mildly nonlinear properties of the largest, collapsed objects within the simulation.
Methods: Based on the combination of several, publicly available surveys we compiled a sample of galaxy clusters within the local Universe, of which we were able to cross-identify 46 of them with an associated counterpart within the SLOW simulation. We then derived the probability of the cross identification based on mass, X-ray luminosity, temperature, and Compton-y by comparing it to a random selection.
Results: Our set of 46 cross-identified local Universe clusters contains the 13 most massive clusters from the Planck SZ catalog as well as 70% of clusters with M₅₀₀ larger than 2 × 10¹⁴ M_⊙. Compared to previous constrained simulations of the local volume, we found in SLOW a much larger amount of replicated galaxy clusters, where their simulation-based mass prediction falls within the uncertainties of the observational mass estimates. Comparing the median observed and simulated masses of our cross-identified sample allows us to independently deduce a hydrostatic mass bias of (1 − b)≈0.87.
Conclusions: The SLOW constrained simulation of the local Universe faithfully reproduces numerous fundamental characteristics of a sizable number of galaxy clusters within our local neighborhood, opening a new avenue for studying the formation and evolution of a large set of individual galaxy clusters as well as testing our understanding of physical processes governing the ICM.

Context. The metallicity distribution function (MDF) of the Galactic bulge is characterized by a multi-peak shape, with a metal-poor peak centered at [Fe/H] ∼ −0.3 dex and a metal-rich peak centered at [Fe/H] ∼ +0.3 dex. The bimodality of the MDF is also reflected in the [α/Fe] versus [Fe/H] abundance ratios, suggesting the presence of different stellar populations in the bulge.
Aims: In this work we aim to reproduce the observed MDF of the Galactic bulge by testing a scenario in which the metal-poor component of the bulge is formed by stars formed in situ, during a strong burst of star formation, while the metal-rich population is formed by stars created in situ during a second burst of star formation and/or stars accreted from the innermost part of the Galactic disk as an effect of a growing bar.
Methods: We adopted a chemical evolution model that is able to follow the evolution of several chemical species with detailed nucleosynthesis prescriptions. In particular, because of the importance of the production of Fe in constraining the MDF, close attention is paid to the production of this element in both Type Ia supernovae and massive stars. In particular, we included yields from rotating massive stars with different rotational velocity prescriptions. Our model also takes the infall and outflow of gas into account, as well as the effect of stellar migration. Results are compared to ∼13 000 stars from the SDSS/APOGEE survey that belong to the region located at a Galactocentric distance R_GC ≤ 3.5 kpc.
Results: We successfully reproduce the observed double-peak shape of the bulge MDF as well as the abundance trends of the α elements relative to Fe by assuming both (i) a multi-burst star formation history with a quenching of the first burst of ∼10² Myr and (ii) migration of stars from the innermost part of the Milky Way disk, as an effect of a growing bar. According to our results, the fraction of the stellar mass of the bulge-bar that belongs to the inner disk is ∼40%. In terms of the nucleosynthesis, we conclude that models that assume either no rotation for massive stars or a distribution of rotational velocities that favors slow rotation at high metallicities best reproduce the observed MDF as well as the [α/Fe] and the [Ce/Fe] versus [Fe/H] abundance patterns.

Traditionally, weak lensing cosmological surveys have been analyzed using summary statistics motivated by their analytically tractable likelihoods, or by their ability to access higher-order information, at the cost of requiring Simulation-Based Inference (SBI) approaches. While informative, these statistics are neither designed nor guaranteed to be statistically sufficient. With the rise of deep learning, it becomes possible to create summary statistics optimized to extract the full data information. We compare different neural summarization strategies proposed in the weak lensing literature, to assess which loss functions lead to theoretically optimal summary statistics to perform full-field inference. In doing so, we aim to provide guidelines and insights to the community to help guide future neural-based inference analyses. We design an experimental setup to isolate the impact of the loss function used to train neural networks. We have developed the sbi_lens JAX package, which implements an automatically differentiable lognormal wCDM LSST-Y10 weak lensing simulator. The explicit full-field posterior obtained using the Hamilotnian-Monte-Carlo sampler gives us a ground truth to which to compare different compression strategies. We provide theoretical insight into the loss functions used in the literature and show that some do not necessarily lead to sufficient statistics (e.g. Mean Square Error (MSE)), while those motivated by information theory (e.g. Variational Mutual Information Maximization (VMIM)) can. Our numerical experiments confirm these insights and show, in our simulated wCDM scenario, that the Figure of Merit (FoM) of an analysis using neural summaries optimized under VMIM achieves 100% of the reference Omega_c - sigma_8 full-field FoM, while an analysis using neural summaries trained under MSE achieves only 81% of the same reference FoM.

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.

The energy injection through Hawking evaporation has been used to put strong constraints on primordial black holes as a dark matter candidate at masses below 10¹⁷ g. However, Hawking's semiclassical approximation breaks down at latest after half-decay. Beyond this point, the evaporation could be significantly suppressed, as was shown in recent work. In this study we review existing cosmological and astrophysical bounds on primordial black holes, taking this effect into account. We show that the constraints disappear completely for a reasonable range of parameters, which opens a new window below 10¹⁰ g for light primordial black holes as a dark matter candidate.

Context. Identifying past wet merger activity in galaxies has been a longstanding issue in extragalactic formation history studies. Gaia's 6D kinematic measurements in our Milky Way (MW) have vastly extended the possibilities for Galactic archaeology, leading to the discovery of a multitude of early mergers in the MW's past. As recent work has established a link between younger globular clusters (GCs; less than about 10-11 Gyr old) and wet galaxy merger events, the MW provides an ideal laboratory for testing which GC properties can be used to trace extragalactic galaxy formation histories.
Aims: To test the hypothesis that GCs trace wet mergers, we relate the measured GC age distributions of the MW and three nearby galaxies, M 31, NGC 1407, and NGC 3115, to their merger histories and interpret the connection with wet mergers through an empirical model for GC formation.
Methods: The GC ages of observed galaxies are taken from a variety of studies to analyze their age distributions side-by-side with the model. For the MW, we additionally cross-match the GCs with their associated progenitor host galaxies to disentangle the connection to the GC age distribution. For the modeled GCs, we take galaxies with similar GC age distributions as observed to compare their accretion histories with those inferred through observations.
Results: We find that the MW GC age distribution is bimodal, mainly caused by younger GCs (10-11 Gyr old associated with Gaia-Sausage/Enceladus (GSE) and in part by unassociated high-energy GCs. The GSE GC age distribution also appears to be bimodal. We propose that the older GSE GCs (12-13 Gyr old) were accreted together with GSE, while the younger ones formed as a result of the merger. For the nearby galaxies, we find that clear peaks in the GC age distributions coincide with active early gas-rich merger phases. Even small signatures in the GC age distributions agree well with the expected wet formation histories of the galaxies inferred through other observed tracers. From the models, we predict that the involved cold gas mass can be estimated from the number of GCs found in the formation burst.
Conclusions: Multimodal GC age distributions can trace massive wet mergers as a result of GCs being formed through them. From the laboratory of our own MW and nearby galaxies we conclude that the ages of younger GC populations of galaxies can be used to infer the wet merger history of a galaxy.

Massive stars are one of the most important and investigated astrophysical production sites of $^{26}$Al, a short-lived radioisotope with $\sim$ 1 Myr half-life. Its short lifetime prevents us from observing its complete chemical history, and only the $^{26}$Al that was recently produced by massive stars can be observed. Hence, it is considered a tracer of star formation rate (SFR). However, important contributions to $^{26}$Al comes from nova systems that pollute the interstellar medium with a large delay, thus partly erasing the correlation between $^{26}$Al and SFR. In this work we describe the 2D distribution of the mass of $^{26}$Al as well as that of massive stars and nova systems in the Milky Way, to investigate their relative contributions to the production of $^{26}$Al. We use a detailed 2D chemical evolution model where the SFR is azimuthally dependent and is required to reproduce the spiral arm pattern observed in the Milky Way. We test two different models, one where the $^{26}$Al comes from massive stars and novae, and one with massive stars only. We then compare the predictions to the $\sim$ 2 M$_{\odot}$ of $^{26}$Al mass observed by the surveys COMPTEL and INTEGRAL. The results show that novae do not trace SFR and, in the solar vicinity, they concentrate in its minima. The effect of novae on the map of the $^{26}$Al mass consists in damping the spiral pattern by a factor of five. Regarding the nucleosynthesis, we find that $\sim$75% of the $^{26}$Al is produced by novae and the $\sim$25% by massive stars. We conclude that novae cannot be neglected as $^{26}$Al producers since the observations can only be reproduced by including their contribution. Moreover, we suggest that bulge novae should eject around six times more material than the disc ones to well reproduce the observed mass of $^{26}$Al.

How does molecular complexity emerge and evolve during the process leading to the formation of a planetary system? Astrochemistry is experiencing a golden age, marked by significant advancements in the observation and understanding of the chemical processes occurring in the inner regions of protostellar systems. However, many questions remain open, such as the origin of the chemical diversity observed in the early evolutionary stages, which may influence the chemical composition of the forming planets. Additionally, astrochemistry provides us with powerful tools to investigate the accretion/ejection processes occurring in the inner regions of young embedded objects, such as jets, winds, accretion streamers, and shocks. In this chapter, we review the observational efforts carried out in recent years to chemically characterize the inner regions of Solar-System analogs. We summarize our current understanding of molecular complexity in planet-forming disks and shed light on the existing limitations and unanswered questions. Finally, we highlight the important role of future radio facilities, like SKAO and ngVLA, in exploring the chemical complexity of the regions where planetary systems are emerging.

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

Context. The structure and kinematics of the old component of the Galactic bulge are still a matter of debate. It is clear that the bulk of the bulge as traced by red clump stars includes two main components, which are usually identified as the metal-rich and metal-poor components. They have different shapes, kinematics, mean metallicities, and alpha-element abundances. It is our current understanding that they are associated with a bar and a spheroid, respectively. On the other hand, RR Lyrae variables trace the oldest population of the bulge. While it would be natural to think that they follow the structure and kinematics of the metal-poor component, the data analysed in the literature show conflicting results.
Aims: We aim to derive a rotation curve for bulge RR Lyrae stars in order to determine that the old component traced by these stars is distinct from the two main components observed in the Galactic bulge.
Methods: This paper combines APOGEE-2S spectra with OGLE-IV light curves, near-infrared photometry, and proper motions from the VISTA Variables in the Vía Láctea survey for 4193 RR Lyrae stars. Six-dimensional phase-space coordinates were used to calculate orbits within an updated Galactic potential and to isolate the stars.
Results: The stars that stay confined within the bulge represent 57% of our sample. Our results show that bulge RR Lyrae variables rotate more slowly than metal-rich red clump stars and have a lower velocity dispersion. Their kinematics is compatible with them being the low-metallicity tail of the metal-poor component. We confirm that a rather large fraction of halo and thick disc RR Lyrae stars pass by the bulge within their orbits, increasing the velocity dispersion. A proper orbital analysis is therefore critical to isolate bona fide bulge variables. Finally, bulge RR Lyrae seem to trace a spheroidal component, although the current data do now allow us to reach a firm conclusion about the spatial distribution.

Full Tables 1 and 2 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/687/A312

Within the relativistic Brueckner-Hartree-Fock theory in the full Dirac space, the tensor-force effects on infinite nuclear matter are elucidated by subtracting the matrix elements of tensor forces from the realistic nucleon-nucleon interaction. The tensor-force effects for the binding energy per particle of symmetric nuclear matter (SNM) as well as the symmetry energy are attractive and are more pronounced around the empirical saturation density, while the tensor forces have little impact on the pure neutron matter. By tuning the tensor-force strength, an infinite (negative) scattering length in the spin-triplet channel is found. This locates the dilute SNM with only the $^3S_1$-$^3D_1$ channel interaction at the unitary limit. Its ground-state energy is found proportional to the energy of a free Fermi gas with a scaling factor 0.38, revealing good universal properties. This work paves the way to study the tensor-force effects in neutron stars as well as finite nuclei from realistic nucleon-nucleon interactions, highlights the role of the tensor force on the deviation of the nuclear physics to the unitary limit, and provides valuable reference for studies of the four-component unitary Fermi gas.

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.

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.

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.