The leading and next-to-leading order QCD predictions for Higgs boson pair production at hadron colliders suffer from a significant mass renormalisation scheme uncertainty related to the choice of the top-quark mass. The functional dependence of the result on the value of the intermediate quark mass can be understood in the high-energy limit using the Method of Regions and the tools of Soft-Collinear Effective Theory. In this work, we study the origin of the sizeable logarithmic mass corrections in the $gg \to HH$ amplitudes at leading and next-to-leading power in the limit $s,|t|,|u| \gg m_t^2 \gg m_H^2$. We argue that the mass corrections follow a predictable factorised pattern that can be exploited to simplify their computation. We present results for the leading power leading logarithmic corrections, our analysis leads to a significant reduction in the theoretical uncertainty of the double Higgs production amplitudes at high-energy due to the top-quark mass scheme.

Dense neutrino gases can exhibit collective flavor instabilities, triggering large flavor conversions that are driven primarily by neutrino-neutrino refraction. One broadly distinguishes between fast instabilities that exist in the limit of vanishing neutrino masses, and slow ones, that require neutrino mass splittings. In a related series of papers, we have shown that fast instabilities result from the resonant growth of flavor waves, in the same way as turbulent electric fields in an unstable plasma. Here we extend this framework to slow instabilities, focusing on the simplest case of an infinitely homogeneous medium with axisymmetric neutrino distribution. The relevant length and time scales are defined by three parameters: the vacuum oscillation frequency $\omega_E=\delta m^2/2E$, the scale of neutrino-neutrino refraction energy $\mu=\sqrt{2}G_F(n_\nu+n_{\overline\nu})$, and the ratio between lepton and particle number $\epsilon=(n_\nu-n_{\overline\nu})/(n_\nu+n_{\overline\nu})$. We distinguish between two very different regimes: (i) For $\omega_E\ll \mu \epsilon^2$, instabilities occur at small spatial scales of order $(\mu\epsilon)^{-1}$ with a time scale of order $\epsilon \omega_E^{-1}$. This novel branch of slow instability arises from resonant interactions with neutrinos moving along the axis of symmetry. (ii) For $\mu \epsilon^2\ll \omega_E\ll \mu$, the instability is strongly non-resonant, with typical time and length scales of order $1/\sqrt{\omega_E \mu}$. Unstable modes interact with all neutrino directions at once, recovering the characteristic scaling of the traditional studies of slow instabilities. In the inner regions of supernovae and neutron-star mergers, the first regime may be more likely to appear, meaning that slow instabilities in this region may have an entirely different character than usually envisaged.

Making inferences about physical properties of the Universe requires knowledge of the data likelihood. A Gaussian distribution is commonly assumed for the uncertainties with a covariance matrix estimated from a set of simulations. The noise in such covariance estimates causes two problems: it distorts the width of the parameter contours, and it adds scatter to the location of those contours which is not captured by the widths themselves. For non-Gaussian likelihoods, an approximation may be derived via Simulation-Based Inference (SBI). It is often implicitly assumed that parameter constraints from SBI analyses, which do not use covariance matrices, are not affected by the same problems as parameter estimation with a covariance matrix estimated from simulations. We investigate whether SBI suffers from effects similar to those of covariance estimation in Gaussian likelihoods. We use Neural Posterior and Likelihood Estimation with continuous and masked autoregressive normalizing flows for density estimation. We fit our approximate posterior models to simulations drawn from a Gaussian linear model, so that the SBI result can be compared to the true posterior. We test linear and neural network based compression, demonstrating that neither methods circumvent the issues of covariance estimation. SBI suffers an inflation of posterior variance that is equal or greater than the analytical result in covariance estimation for Gaussian likelihoods for the same number of simulations. The assumption that SBI requires a smaller number of simulations than covariance estimation for a Gaussian likelihood analysis is inaccurate. The limitations of traditional likelihood analysis with simulation-based covariance remain for SBI with a finite simulation budget. Despite these issues, we show that SBI correctly draws the true posterior contour given enough simulations.

Transmission spectroscopy is a prime method to study the atmospheres of extrasolar planets. We obtained a high-resolution spectral transit time series of the hot Jupiter WASP-121 b with CRIRES+ to study its atmosphere via transmission spectroscopy of the He I λ10833 triplet lines. Our analysis shows a prominent He I λ10833 absorption feature moving along with the planetary orbital motion, which shows an observed, transit-averaged equivalent width of approximately 30 mÅ, a slight redshift, and a depth of about 2%, which can only be explained by an atmosphere overflowing its Roche lobe. We carried out 3D hydrodynamic modeling to reproduce the observations, which favors asymmetric mass loss with a more pronounced leading tidal tail, possibly also explaining observational evidence for additional absorption stationary in the stellar rest frame. A trailing tail is not detectable. From our modeling, we derived estimates of ≥2 × 10¹³ g s^‑1 for the stellar and 5.4 × 10¹² g s^‑1 for the planetary mass loss rate, which is consistent with X-ray and extreme-ultraviolet (XUV) driven mass loss in WASP-121 b.

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

The impact of the dynamical state of gas-rich satellite galaxies at the early moments of their infall into their host systems and the relation to their quenching process are not completely understood at the low-mass regime. Two such nearby systems are the infalling Milky Way (MW) dwarfs Leo T and Phoenix located near the MW virial radius at 414 kpc (1.4R_vir), both of which present intriguing offsets between their gaseous and stellar distributions. Here we present hydrodynamic simulations with RAMSES to reproduce the observed dynamics of Leo T: its 80 pc stellar-HI offset and the 35 pc offset between its older (≳5 Gyr) and younger (∼200 ‑ 1000 Myr) stellar population. We considered internal and environmental properties such as stellar winds, two HI components, cored and cuspy dark matter profiles, and different satellite orbits considering the MW circumgalactic medium. We find that the models that best match the observed morphology of the gas and stars include mild stellar winds that interact with the HI generating the observed offset, and dark matter profiles with extended cores. The latter allow long oscillations of the off-centred younger stellar component, due to long mixing timescales (≳200 Myr), and the slow precession of near-closed orbits in the cored potentials; instead, cuspy and compact cored dark matter models result in the rapid mixing of the material (≲200 Myr). These models predict that non-equilibrium substructures, such as spatial and kinematic offsets, are likely to persist in cored low-mass dwarfs and to remain detectable on long timescales in systems with recent star formation.

We construct the equation of state of hypernuclear matter and study the structure of neutron stars employing a chiral hyperon-nucleon interaction of the Jülich--Bonn group tuned to femtoscopic $\Lambda p$ data of the ALICE collaboration, and $\Lambda\Lambda$ and $\Xi$N interactions determined from Lattice QCD calculations by the HAL QCD collaboration that reproduce the femtoscopic $\Lambda\Lambda$ and $\Xi^-p$ data. We employ the ab-initio microscopic Brueckner--Hartree--Fock theory extended to the strange baryon sector. A special focus is put on the uncertainties of the hyperon interactions and how they are effectively propagated to the composition, equation of state, and mass-radius relation of neutron stars. To such end, we consider the uncertainty due to the experimental error of the femtoscopic $\Lambda p$ data used to fix the chiral hyperon-nucleon interaction and the theoretical uncertainty, estimated from the residual cut-off dependence of this interaction. We find that the final maximum mass of a neutron star with hyperons is in the range $1.3-1.4$ $M_\odot$, in agreement with previous works. The hyperon puzzle, therefore, remains still an open issue if only two-body hyperon-nucleon and hyperon-hyperon interactions are considered.

We critically reconsider the argument based on 't Hooft anomaly matching that aims at proving chiral symmetry breaking in confining four-dimensional quantum chromodynamics-like theories with <inline-formula><mml:math display="inline"><mml:msub><mml:mi>N</mml:mi><mml:mi>c</mml:mi></mml:msub><mml:mo>></mml:mo><mml:mn>2</mml:mn></mml:math></inline-formula> colors and <inline-formula><mml:math display="inline"><mml:msub><mml:mi>N</mml:mi><mml:mi>f</mml:mi></mml:msub></mml:math></inline-formula> flavors. The main line of reasoning relies on a property of the solutions of the anomaly matching and persistent mass equations called <inline-formula><mml:math display="inline"><mml:msub><mml:mi>N</mml:mi><mml:mi>f</mml:mi></mml:msub></mml:math></inline-formula>-independence. In previous works, the validity of <inline-formula><mml:math display="inline"><mml:msub><mml:mi>N</mml:mi><mml:mi>f</mml:mi></mml:msub></mml:math></inline-formula>-independence was assumed based on qualitative arguments, but it was never proven rigorously. We provide a detailed proof and clarify under which (dynamical) conditions it holds. Our results are valid for a generic spectrum of massless composite fermions including baryons and exotics.

We carried out long-term monitoring of the LIGO/Virgo/KAGRA binary black hole (BBH) merger candidate S230922g in search of electromagnetic emission from the interaction of the merger remnant with an embedding active galactic nuclei (AGN) accretion disk. Using a dataset primarily composed of wide-field imaging from the Dark Energy Camera and supplemented by additional photometric and spectroscopic resources, we searched <inline-formula><mml:math display="inline"><mml:mo>∼</mml:mo><mml:mn>70</mml:mn><mml:mo>%</mml:mo></mml:math></inline-formula> of the sky area probability for transient phenomena and discovered six counterpart candidates. One especially promising candidate—AT 2023aagj—exhibited temporally varying asymmetric components in spectral broad line regions, a feature potentially indicative of an off-center event such as a BBH merger. This represents the first live search and multiwavelength, photometric, and spectroscopic monitoring of a gravitational wave BBH optical counterpart candidate in the disk of an AGN.

The projected sensitivity of the effective electron neutrino-mass measurement with the KATRIN experiment is below 0.3 eV (90 % CL) after 5 years of data acquisition. The sensitivity is affected by the increased rate of the background electrons from KATRIN's main spectrometer. A special shifted-analysing-plane (SAP) configuration was developed to reduce this background by a factor of two. The complex layout of electromagnetic fields in the SAP configuration requires a robust method of estimating these fields. We present in this paper a dedicated calibration measurement of the fields using conversion electrons of gaseous <inline-formula id="IEq1"><mml:math id="IEq1_Math"><mml:mmultiscripts><mml:mrow></mml:mrow><mml:mrow></mml:mrow><mml:mtext>83m</mml:mtext></mml:mmultiscripts></mml:math></inline-formula>Kr, which enables the neutrino-mass measurements in the SAP configuration.

Context. Increasing evidence shows that warped disks are common, challenging the methods used to model their velocity fields. Molecular line emission of these disks is characterized by a twisted pattern, similar to the signal from radial flows, complicating the study of warped disk kinematics. Previous attempts to model these features have encountered difficulties in distinguishing between the underlying kinematics of different disks. Aims. This study aims to advance gas kinematics modeling capabilities by extending the Extracting Disk Dynamics (eddy) package to include warped geometries and radial flows. We assess the performance of eddy in recovering input parameters for scenarios involving warps, radial flows, and combinations of the two. Additionally, we provide a basis to break the visual degeneracy between warped disks and radial flow, establishing a criterion to distinguish them. Methods. We extended the eddy package to handle warped geometries by including a parametric prescription of a warped disk and a ray-casting algorithm to account for the surface self-obscuration arising from the 3D to 2D projection. The effectiveness of the tool was tested using the radiative transfer code RADMC3D, generating synthetic models for disks with radial flows, warped disks, and warped disks with radial flows. Results. We demonstrate the efficacy of our tool in accurately recovering the geometrical parameters of systems, particularly in data with sufficient angular resolution. Importantly, we observe minimal impact from thermal noise levels typical in Atacama Large Millimeter/submillimeter Array (ALMA) observations. Furthermore, our findings reveal that fitting an incorrect model type produces characteristic residual signatures, which serve as kinematic criteria for disk classification. Conclusions. Characterizing gas kinematics requires careful consideration of twisted motions. While our model provides insights into disk geometries, caution is needed when interpreting parameters in regions with complex kinematics or low-resolution data. Future ALMA baseline observations should help clarify warped disk kinematics.

We present a detailed survey of electric, magnetic, and quadrupole form factors of light and heavy spin-1 vector mesons. It complements our analogous analysis of the electromagnetic form factors of pseudoscalar and scalar mesons reported earlier. Our formalism is based upon the Schwinger-Dyson equations treatment of a vector <inline-formula><mml:math display="inline"><mml:mo>×</mml:mo></mml:math></inline-formula> vector contact interaction and the Bethe-Salpeter equation description of relativistic two-body bound states. We compute the form factors, associated moments, and charge radii, comparing these quantities to earlier theoretical studies and experimental results if and when possible. We also investigate the quark-mass dependence of the charge radii and find the anticipated hierarchy such that it decreases with increasing dressed quark masses. In addition, our analysis shows that the magnetic moment is independent of the mass of the light and heavy mesons. Our results agree with most measurements reported earlier, finding a negative quadrupole moment, implying the charge distribution is oblate.

Real-world systems are shaped by both their complex internal interactions and the changes in their noisy environments. In this work, we study how a shared active bath affects the statistical dependencies between two interacting Brownian particles by evaluating their mutual information. We decompose the mutual information into three terms: information stemming from the internal interactions between the particles; information induced by the shared bath, which encodes environmental changes; a term describing information interference that quantifies how the combined presence of both internal interactions and environment either masks (destructive interference) or boosts (constructive interference) information. By studying exactly the case of linear interactions, we find that the sign of information interference depends solely on that of the internal coupling. However, when internal interactions are described by a nonlinear activation function, we show that both constructive and destructive interference appear depending on the interplay between the timescale of the active environment, the internal interactions, and the environmental coupling. Finally, we show that our results generalize to hierarchical systems where asymmetric couplings to the environment mimic the scenario where the active bath is only partially accessible to one particle. This setting allows us to quantify how this asymmetry drives information interference. Our work underscores how information and functional relationships in complex multiscale systems are fundamentally shaped by the environmental context.

Studying the orbital motion of stars around Sagittarius A* in the Galactic center provides a unique opportunity to probe the gravitational potential near the supermassive black hole at the heart of our Galaxy. Interferometric data obtained with the GRAVITY instrument at the Very Large Telescope Interferometer (VLTI) since 2016 has allowed us to achieve unprecedented precision in tracking the orbits of these stars. GRAVITY data have been key to detecting the in-plane, prograde Schwarzschild precession of the orbit of the star S2 that is predicted by general relativity. By combining astrometric and spectroscopic data from multiple stars, including S2, S29, S38, and S55 – for which we have data around their time of pericenter passage with GRAVITY – we can now strengthen the significance of this detection to an approximately 10σ confidence level. The prograde precession of S2's orbit provides valuable insights into the potential presence of an extended mass distribution surrounding Sagittarius A*, which could consist of a dynamically relaxed stellar cusp comprising old stars and stellar remnants, along with a possible dark matter spike. Our analysis, based on two plausible density profiles – a power-law and a Plummer profile – constrains the enclosed mass within the orbit of S2 to be consistent with zero, establishing an upper limit of approximately 1200 M_⊙ with a 1σ confidence level. This significantly improves our constraints on the mass distribution in the Galactic center. Our upper limit is very close to the expected value from numerical simulations for a stellar cusp in the Galactic center, leaving little room for a significant enhancement of dark matter density near Sagittarius A*.

We present a spectroscopic analysis of MACS J0138$-$2155, at $z=0.336$, the first galaxy cluster hosting two strongly-lensed supernovae (SNe), Requiem and Encore, providing us with a chance to obtain a reliable $H_0$ measurement from the time delays between the multiple images. We take advantage of new data from the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope, covering a central $1 \rm \, arcmin^2$ of the lensing cluster, for a total depth of 3.7 hours, including 2.9 hours recently obtained by our Target of Opportunity programme. Our new spectroscopic catalogue contains reliable redshifts for 107 objects, including 50 galaxy cluster members with secure redshift values in the range $0.324 < z < 0.349$, and 13 lensed multiple images from four background sources between $0.767\leq z \leq 3.420$, including four images of the host galaxy of the two SNe. We exploit the MUSE data to study the stellar kinematics of 14 bright cluster members and two background galaxies, obtaining reliable measurements of their line-of-sight velocity dispersion. Finally, we combine these results with measurements of the total magnitude of the cluster members in the Hubble Space Telescope F160W band to calibrate the Faber-Jackson relation between luminosity and stellar velocity dispersion ($L \propto \sigma^{1/\alpha}$) for the early-type cluster member galaxies, measuring a slope $\alpha=0.25^{+0.05}_{-0.05}$. A pure and complete sample of cluster member galaxies and a reliable characterisation of their total mass structure are key to building accurate total mass maps of the cluster, mitigating the impact of parametric degeneracies, which is necessary for inferring the value of $H_0$ from the measured time delays between the lensed images of the two SNe.

In modern collider experiments, the quest to explore fundamental interactions between elementary particles has reached unparalleled levels of precision. Signatures from particle physics detectors are low-level objects (such as energy depositions or tracks) encoding the physics of collisions (the final state particles of hard scattering interactions). The complete simulation of them in a detector is a computational and storage-intensive task. To address this computational bottleneck in particle physics, alternative approaches have been developed, introducing additional assumptions and trade off accuracy for speed. The field has seen a surge in interest in surrogate modeling the detector simulation, fueled by the advancements in deep generative models. These models aim to generate responses that are statistically identical to the observed data. In this paper, we conduct a comprehensive and exhaustive taxonomic review of the existing literature on the simulation of detector signatures from both methodological and application-wise perspectives. Initially, we formulate the problem of detector signature simulation and discuss its different variations that can be unified. Next, we classify the state-of-the-art methods into five distinct categories based on their underlying model architectures, summarizing their respective generation strategies. Finally, we shed light on the challenges and opportunities that lie ahead in detector signature simulation, setting the stage for future research and development.

A new and detailed measurement of the cross section for hard exclusive neutral-pion muoproduction on the proton was performed in a wide kinematic region, with the photon virtuality $Q^2$ ranging from 1 to 8 (GeV/$c$)$^{\rm\, 2}$ and the Bjorken variable $x_{\rm Bj}$ ranging from 0.02 to 0.45. The data were collected at COMPASS at CERN using 160 GeV/$c$ longitudinally polarised $\mu^+$ and $\mu^-$ beams scattering off a 2.5 m long liquid hydrogen target. From the average of the measured $\mu^+$ and $\mu^-$ cross sections, the virtual-photon--proton cross section is determined as a function of the squared four-momentum transfer between the initial and final state proton in the range 0.08 (GeV/$c$)$^{\rm\, 2}$$< |t| <$ 0.64 (GeV/$c$)$^{\rm\, 2}$. From its angular distribution, the combined contribution of transversely and longitudinally polarised photons are determined, as well as transverse--transverse and longitudinal--transverse interference contributions. They are studied as functions of four-momentum transfer $|t|$, photon virtuality $Q^2$ and virtual-photon energy $\nu$. The longitudinal--transverse interference contribution is found to be compatible with zero. The significant transverse--transverse interference contribution reveals the existence of a dominant contribution by transversely polarized photons. This provides clear experimental evidence for the chiral-odd GPD $\overline{E}_T$. In addition, the existence of a non-negligible contribution of longitudinally polarized photons is suggested by the $|t|$-dependence of the cross section at $x_{\rm Bj} < $ 0.1 . Altogether, these results provide valuable input for future modelling of GPDs and thus of cross sections for exclusive pseudo-scalar meson production. Furthermore, they can be expected to facilitate the study of next-to-leading order corrections and higher-twist contributions.

The Intra Cluster Medium (ICM) of Galaxy Clusters (GCs) is a highly dynamic environment. It is shaped by mergers and bulk motions on large scales. Small scales are dominated by turbulence, connected to large scales via a turbulent cascade. Both observations and simulations predict this turbulence to be subsonic. This poses numerical challenges for calculating the turbulent cascade and requires careful numerical treatment of hydrodynamics.

Many different numerical methods have been developed and applied to this specific problem. They can be divided according to their discretization approach into grid-based volume-discretization methods, such as stationary and moving meshes, and mass-discretization methods such as Smoothed Particle Hydrodynamics (SPH). More recently, Meshless Finite Mass (MFM) has been developed. The gas is discretized by mass, but fluxes between neighbors are calculated, thus combining the advantages of SPH with grid-based methods.

In this work, we present a new implementation of MFM in the cosmological simulation code OpenGadget3. It is based on the implementation in the Gandalf code but has been extended to allow for cosmological applications. One main goal is the application to subsonic turbulence in the ICM of GCs and a detailed and fair comparison with other hydrodynamical methods. [...]

We calculate the four-graviton scattering amplitude in Type II superstring theory at one-loop up to seventh order in the low-energy expansion through the recently developed iterated integral formalism of Modular Graph Functions (MGFs). We propose a new assignment of transcendental weight to the numbers that appear in the amplitude, which leads to a violation of uniform transcendentality. Furthermore, the machinery of the novel method allows us to propose a general form of the amplitude, which suggests that the expansion is expressible in terms of single-valued multiple zeta values and logarithmic derivatives of the Riemann zeta function at positive and negative odd integers.

This chapter reviews the three-dimensional structure, age, kinematics, and chemistry of the Milky Way (MW) region within ~2 kpc from its center (hereafter referred to as the 'bulge') from an observational perspective. While not exhaustive in citations, this review provides historical context and discusses the main controversies and limitations in the current consensus. The nuclear bulge region, within $\sim$200 pc from the Galactic center, has been excluded from this review. This very complex region, hosting dense molecular clouds and active star formation, would deserve a dedicated paper.

Accretion-induced collapse (AIC) or merger-induced collapse (MIC) of white dwarfs (WDs) in binary systems is an interesting path to neutron star (NS) and magnetar formation, alternative to stellar core collapse and NS mergers. Such events could add a population of compact remnants in globular clusters, they are expected to produce yet unidentified electromagnetic transients including gamma-ray and radio bursts, and to act as sources of trans-iron elements, neutrinos, and gravitational waves. Here we present the first long-term (>5 s post bounce) hydrodynamical simulations in axi-symmetry (2D), using energy- and velocity-dependent three-flavor neutrino transport based on a two-moment scheme. Our set of six models includes initial WD configurations for different masses, central densities, rotation rates, and angular momentum profiles. Our simulations demonstrate that rotation plays a crucial role for the proto-neutron star (PNS) evolution and ejecta properties. We find early neutron-rich ejecta and an increasingly proton-rich neutrino-driven wind at later times in a non-rotating model, in agreement with electron-capture supernova models. In contrast to that and different from previous results, our rotating models eject proton-rich material initially and increasingly more neutron-rich matter as time advances, because an extended accretion torus forms around the PNS and feeds neutrino-driven bipolar outflows for many seconds. AIC and MIC events are thus potential sites of r-process element production, which may imply constraints on their occurrence rates. Finally, our simulations neglect the effects of triaxial deformation and magnetic fields, serving as a temporary benchmark for more comprehensive future studies.

Chemical clocks based on [s-process elements/alpha-elements] ratios are widely used to estimate ages of Galactic stellar populations. However, the [s/alpha] vs. age relations are not universal, varying with metallicity, location in the Galactic disc, and specific s-process elements. Current Galactic chemical evolution models struggle to reproduce the observed [s/alpha] increase at young ages. We provide chemical evolution models for the Milky Way disc to identify the conditions required to reproduce the observed [s/H], [s/Fe], and [s/alpha] vs. age relations. We adopt a multi-zone chemical evolution model including state-of-the-art nucleosynthesis prescriptions for neutron-capture elements (AGB stars, rotating massive stars, neutron star mergers, magneto-driven supernovae). We explore variations in gas infall, AGB yield dependencies on progenitor stars, and rotational velocity distributions for massive stars. Results are compared with open cluster data from the Gaia-ESO survey. A three-infall scenario for disc formation captures the rise of [s/alpha] with age in the outer regions but fails in the inner ones, especially for second s-process peak elements. Ba production in the last 3 Gyr of chemical evolution would need to increase by half to match observations. S-process contributions from low-mass AGB stars improve predictions but require increases not supported by nucleosynthesis calculations, even with potential i-process contribution. Variations in the metallicity dependence of AGB yields show inconsistent effects across elements. Distributions of massive star rotational velocities fail to improve results due to balanced effects on elements. We confirm that there is no single relationship [s/alpha] vs. age, but that it varies along the MW disc. Current prescriptions for neutron-capture element yields cannot fully capture the complexity of evolution, particularly in the inner disc.

Superclusters are the most massive structures in the universe. To what degree they are actually bound against an accelerating expansion of the background is of significant cosmological and astrophysical interest. In this study, we introduce a cross matched set of superclusters from the SLOW constrained simulations of the local (z<0.05) universe. Identifying the superclusters provides estimates on the efficacy of the constraints in reproducing the local large-scale structure accurately. The simulated counterparts can help identifying possible future observational targets containing interesting features such as bridges between pre-merging and merging galaxy clusters and collapsing filaments and provide comparisons for current observations. By determining the collapse volumes for the superclusters we further elucidate the dynamics of cluster-cluster interactions in those regions. Using catalogs of local superclusters and the most massive simulated clusters, we search for counterparts of supercluster members of six regions. We evaluate the significance of these detections by comparing their geometries to supercluster regions in random simulations. We then run an N-body version of the simulation into the far future and determine which of the member clusters are gravitationally bound to the host superclusters. Furthermore we compute masses and density contrasts for the collapse regions. We demonstrate the SLOW simulation of the local universe to accurately reproduce local supercluster regions in mass of their members and three-dimensional geometrical arrangement. We furthermore find the bound regions of the local superclusters consistent in size and density contrast with previous theoretical studies. This will allow to connect future numerical zoom-in studies of the clusters to the large scale environments and specifically the supercluster environments these local galaxy clusters evolve in.

We update the Standard Model (SM) predictions for the lifetimes of the $B^+$, $B_d$ and $B_s$ mesons within the heavy quark expansion (HQE), including the recently determined NNLO-QCD corrections to non-leptonic decays of the free $b$-quark. In addition, we update the HQE predictions for the lifetime ratios $\tau (B^+)/\tau (B_d)$ and $\tau (B_s)/\tau (B_d)$, and provide new results for the semileptonic branching fractions of the three mesons entirely within the HQE. We obtain a considerable improvement of the theoretical uncertainties, mostly due to the reduction of the renormalisation scale dependence when going from LO to NNLO, and for all the observables considered, we find good agreement, within uncertainties, between the HQE predictions and the corresponding experimental data. Our results read, respectively, $\Gamma (B^+) = 0.587^{+0.025}_{-0.035}~{\rm ps}^{-1}$, $\Gamma (B_d) = 0.636^{+0.028}_{-0.037}~{\rm ps}^{-1}$, $\Gamma (B_s) = 0.628^{+0.027}_{-0.035}~{\rm ps}^{-1}$, for the total decay widths, $\tau (B^+)/\tau (B_d) = 1.081^{+0.014}_{-0.016}$, $\tau (B_s)/\tau (B_d) = 1.013^{+0.007}_{-0.007}$, for the lifetime ratios, and ${\cal B}_{\rm sl} (B^+) = (11.46^{+0.47}_{-0.32}) \%$, ${\cal B}_{\rm sl} (B_d) = (10.57^{+0.47}_{-0.27}) \%$, ${\cal B}_{\rm sl} (B_s) = (10.52^{+0.50}_{-0.29}) \%$, for the semileptonic branching ratios. Finally, we also provide an outlook for further improvements of the HQE determinations of the $B$-meson decay widths and of their ratios.

We construct defects describing the transition between different phases of gauged linear sigma models with higher rank abelian gauge groups, as well as defects embedding these phases into the GLSMs. Our construction refers entirely to the sector protected by B-type supersymmetry, decoupling the gauge sector. It relies on an abstract characterization of such transition defects and does not involve an actual perturbative analysis. It turns out that the choices that are required to characterize consistent transition defects match with the homotopy classes of paths between different phases. Our method applies to non-anomalous as well as anomalous GLSMs, and we illustrate both cases with examples. This includes the GLSM associated to the resolution of the $A_N$ singularity and one describing the entire parameter space of $N = 2$ minimal models, in particular, the relevant flows between them. Via fusion with boundary conditions, the defects we construct yield functors describing the transport of D-branes on parameter space. We find that our results match with known results on D-brane transport.

Strongly lensed Type Ia supernovae (LSNe Ia) are a promising probe with which to measure the Hubble constant (H₀) directly. To use LSNe Ia for cosmography, a time-delay measurement between multiple images, a lens-mass model, and a mass reconstruction along the line of sight are required. In this work, we present the machine-learning network LSTM-FCNN, which is a combination of a long short-term memory network (LSTM) and a fully connected neural network (FCNN). The LSTM-FCNN is designed to measure time delays on a sample of LSNe Ia spanning a broad range of properties, which we expect to find with the upcoming Rubin Observatory Legacy Survey of Space and Time (LSST) and for which follow-up observations are planned. With follow-up observations in the i band (cadence of one to three days with a single-epoch 5σ depth of 24.5 mag), we reach a bias-free delay measurement with a precision of around 0.7 days over a large sample of LSNe Ia. The LSTM-FCNN is far more general than previous machine-learning approaches such as the random forest (RF) one, whereby an RF has to be trained for each observational pattern separately, and yet the LSTM-FCNN outperforms the RF by a factor of roughly three. Therefore, the LSTM-FCNN is a very promising approach to achieve robust time delays in LSNe Ia, which is important for a precise and accurate constraint on H₀.

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

Submillimeter single-dish telescopes offer two key advantages compared to interferometers: they can efficiently map larger portions of the sky and recover larger spatial scales. Nonetheless, fluctuations in the atmosphere limit the accurate retrieval of signals from astronomical sources. Therefore, we introduce a user-friendly simulator named maria to optimize scanning strategies and instrument designs to efficiently reduce atmospheric noise and filtering effects. We further use this tool to produce synthetic time streams and maps from hydrodynamical simulations, enabling a fair comparison between theory and reality. maria has implemented a suite of telescope and instrument designs intended to mimic current and future facilities. To generate synthetic time-ordered data, each mock observatory scans through the atmosphere in a configurable pattern over the celestial object. We generate evolving and location-and-time-specific weather for each of the fiducial sites using a combination of satellite and ground-based measurements. While maria is a generic virtual telescope, this study specifically focuses on mimicking broadband bolometers observing at 100 GHz. To validate our virtual telescope, we compare the mock time streams with real MUSTANG-2 observations and find that they are quantitatively similar by conducting a k-sample Anderson-Darling test resulting in p<0.001. Subsequently, we image the time-ordered data to create noise maps and mock observations of clusters of galaxies for both MUSTANG-2 and an instrument concept for the 50m Atacama Large Aperture Submillimeter Telescope (AtLAST). Furthermore, using maria, we find that a 50m dish provides the highest levels of correlation of atmospheric signals across adjacent detectors compared to smaller apertures (e.g., 42-cm and 6-m survey experiments), facilitating removal of atmospheric signal on large scales.

Nonlinearities in King plots (KP) of isotope shifts (IS) can reveal the existence of beyond-Standard-Model (BSM) interactions that couple electrons and neutrons. However, it is crucial to distinguish higher-order Standard Model (SM) effects from BSM physics. We measure the IS of the transitions ${{}^{3}P_{0}~\rightarrow~{}^{3}P_{1}}$ in $\mathrm{Ca}^{14+}$ and ${{}^{2}S_{1/2} \rightarrow {}^{2}D_{5/2}}$ in $\mathrm{Ca}^{+}$ with sub-Hz precision as well as the nuclear mass ratios with relative uncertainties below $4\times10^{-11}$ for the five stable, even isotopes of calcium (${}^{40,42,44,46,48}\mathrm{Ca}$). Combined, these measurements yield a calcium KP nonlinearity with a significance of $\sim 900 \sigma$. Precision calculations show that the nonlinearity cannot be fully accounted for by the expected largest higher-order SM effect, the second-order mass shift, and identify the little-studied nuclear polarization as the only remaining SM contribution that may be large enough to explain it. Despite the observed nonlinearity, we improve existing KP-based constraints on a hypothetical Yukawa interaction for most of the new boson masses between $10~\mathrm{eV/c^2}$ and $10^7~\mathrm{eV/c^2}$.

The existence of light QCD axions, whose mass depends on an additional free parameter, can lead to a new ground state of matter, where the sourced axion field reduces the nucleon effective mass. The presence of the axion field has structural consequences, in particular, it results in a thinner (or even prevents its existence) heat-blanketing envelope, significantly altering the cooling patterns of neutron stars. We exploit the anomalous cooling behavior to constrain previously uncharted regions of the axion parameter space by comparing model predictions with existing data from isolated neutron stars. Notably, this analysis does not require the light QCD axion to be the dark matter candidate.

We perform an all-order analysis of double-logarithmic corrections to the so-called soft-overlap contribution to heavy-to-light transition form factors at large hadronic recoil. Specifically, we study $B_c \to \eta_c$ transitions within a perturbative non-relativistic framework, treating both the bottom and charm quarks as heavy with the hierarchy $m_b \gg m_c \gg\Lambda_{\rm QCD}$. Our diagrammatic analysis shows that double-logarithmic corrections arise from two distinct sources: Exponentiated soft-gluon effects described by standard Sudakov factors, and rapidity-ordered soft-quark configurations, leading to implicit integral equations, which so far have only been studied in the context of energetic muon-electron backward scattering. We find that the all-order structure of the double logarithms is governed by a novel type of coupled integral equations, which encode the non-trivial interplay between these two effects. Whereas a closed-form solution to these equations is currently unknown, we present useful iteration formulas, and derive the asymptotic behaviour of the soft-overlap form factor for infinitely large recoil energies, showing that the Sudakov suppression is somewhat weakened by the intertwined soft-quark and soft-gluon corrections. In a broader context, our findings shed light onto the physical origin and mathematical structure of endpoint divergences arising from soft-collinear factorization and the related Feynman mechanism for power-suppressed hard exclusive processes.

We study the femtoscopic correlation functions of meson-baryon pairs in the strangeness $S=-1$ sector, employing unitarized s-wave scattering amplitudes derived from the chiral Lagrangian up to next-to-leading order. For the first time, we deliver predictions on the $\pi^-\Lambda$ and $K^+\Xi^-$ correlation functions which are feasible to be measured at the Large Hadron Collider. We also demonstrate that the employed model is perfectly capable of reproducing the $K^-p$ correlation function data measured by the same collaboration, without the need to modify the coupling strength to the $\bar{K}^0n$ channel, as has been recently suggested. In all cases, the effects of the source size on the correlation are tested. In addition, we present detailed analysis of the different coupled-channel contributions, together with the quantification of the relative relevance of the different terms in the interaction.

Rotation matters for the life of a star. It causes a multitude of dynamical phenomena in the stellar interior during a star's evolution, and its effects accumulate until the star dies. All stars rotate at some level, but most of those born with a mass higher than 1.3 times the mass of the Sun rotate rapidly during more than 90% of their nuclear lifetime. Internal rotation guides the angular momentum and chemical element transport throughout the stellar interior. These transport processes change over time as the star evolves. The cumulative effects of stellar rotation and its induced transport processes determine the helium content of the core by the time it exhausts its hydrogen isotopes. The amount of helium at that stage also guides the heavy element yields by the end of the star's life. A proper theory of stellar evolution and any realistic models for the chemical enrichment of galaxies must be based on observational calibrations of stellar rotation and of the induced transport processes. In the last few years, asteroseismology offers such calibrations for single and binary stars. We review the current status of asteroseismic modelling of rotating stars for different stellar mass regimes in an accessible way for the non-expert. While doing so, we describe exciting opportunities sparked by asteroseismology for various domains in astrophysics, touching upon topics such as exoplanetary science, galactic structure and evolution, and gravitational wave physics to mention just a few. Along the way we provide ample sneak-previews for future 'industrialised' applications of asteroseismology to slow and rapid rotators from the exploitation of combined Kepler, Transiting Exoplanet Survey Satellite (TESS), PLAnetary Transits and Oscillations of stars (PLATO), Gaia, and ground-based spectroscopic and multi-colour photometric surveys. We end the review with a list of takeaway messages and achievements of asteroseismology that are of relevance for many fields of astrophysics.

Aims. Detecting diffuse synchrotron emission from the cosmic web is still a challenge for current radio telescopes. We aim to make predictions about the detectability of cosmic web filaments from simulations. Methods. We present the first cosmological magnetohydrodynamic simulation of a 500 h^‑1 c Mpc volume with an on-the-fly spectral cosmic ray (CR) model. This allows us to follow the evolution of populations of CR electrons and protons within every resolution element of the simulation. We modeled CR injection at shocks, while accounting for adiabatic changes to the CR population and high-energy-loss processes of electrons. The synchrotron emission was then calculated from the aged electron population, using the simulated magnetic field, as well as different models for the origin and amplification of magnetic fields. We used constrained initial conditions, which closely resemble the local Universe, and compared the results of the cosmological volume to a zoom-in simulation of the Coma cluster, to study the impact of resolution and turbulent reacceleration of CRs on the results. Results. We find a consistent injection of CRs at accretion shocks onto cosmic web filaments and galaxy clusters. This leads to diffuse emission from filaments of the order S_ν ≈ 0.1 μJy beam^‑1 for a potential LOFAR observation at 144 MHz, when assuming the most optimistic magnetic field model. The flux can be increased by up to two orders of magnitude for different choices of CR injection parameters. This can bring the flux within a factor of ten of the current limits for direct detection. We find a spectral index of the simulated synchrotron emission from filaments of α ≈ ‑1.0 to –1.5 in the LOFAR band.

Estimates of seismic wave speeds in the Earth (seismic velocity models) are key input parameters to earthquake simulations for ground motion prediction. Owing to the non-uniqueness of the seismic inverse problem, typically many velocity models exist for any given region. The arbitrary choice of which velocity model to use in earthquake simulations impacts ground motion predictions. However, current hazard analysis methods do not account for this source of uncertainty. We present a proof-of-concept ground motion prediction workflow for incorporating uncertainties arising from inconsistencies between existing seismic velocity models. Our analysis is based on the probabilistic fusion of overlapping seismic velocity models using scalable Gaussian process (GP) regression. Specifically, we fit a GP to two synthetic 1-D velocity profiles simultaneously, and show that the predictive uncertainty accounts for the differences between the models. We subsequently draw velocity model samples from the predictive distribution and estimate peak ground displacement using acoustic wave propagation through the velocity models. The resulting distribution of possible ground motion amplitudes is much wider than would be predicted by simulating shaking using only the two input velocity models. This proof-of-concept illustrates the importance of probabilistic methods for physics-based seismic hazard analysis.

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

Context. Computing reliable photometric redshifts (photo-z) for active galactic nuclei (AGN) is a challenging task, primarily due to the complex interplay between the unresolved relative emissions associated with the supermassive black hole and its host galaxy. Spectral energy distribution (SED) fitting methods, while effective for galaxies and AGN in pencil-beam surveys, face limitations in wide or all-sky surveys with fewer bands available, lacking the ability to accurately capture the AGN contribution to the SED, hindering reliable redshift estimation. This limitation is affecting the many tens of millions of AGN detected in existing datasets, such as those AGN clearly singled out and identified by SRG/eROSITA. Aims. Our goal is to enhance photometric redshift performance for AGN in all-sky surveys while simultaneously simplifying the approach by avoiding the need to merge multiple data sets. Instead, we employ readily available data products from the 10th Data Release of the Imaging Legacy Survey for the Dark Energy Spectroscopic Instrument, which covers >20 000 deg² of extragalactic sky with deep imaging and catalog-based photometry in the ɡriɀW1-W4 bands. We fully utilize the spatial flux distribution in the vicinity of each source to produce reliable photo-z. Methods. We introduce PICZL, a machine-learning algorithm leveraging an ensemble of convolutional neural networks. Utilizing a cross-channel approach, the algorithm integrates distinct SED features from images with those obtained from catalog-level data. Full probability distributions are achieved via the integration of Gaussian mixture models. Results. On a validation sample of 8098 AGN, PICZL achieves an accuracy σ_NMAD of 4.5% with an outlier fraction η of 5.6%. These results significantly outperform previous attempts to compute accurate photo-z for AGN using machine learning. We highlight that the model's performance depends on many variables, predominantly the depth of the data and associated photometric error. A thorough evaluation of these dependencies is presented in the paper. Conclusions. Our streamlined methodology maintains consistent performance across the entire survey area, when accounting for differing data quality. The same approach can be adopted for future deep photometric surveys such as LSST and Euclid, showcasing its potential for wide-scale realization. With this paper, we release updated photo-z (including errors) for the XMM-SERVS W-CDF-S, ELAIS-S1 and LSS fields.

The Standard Model extended by a real scalar singlet $S$ with an approximate $\mathbb{Z}_2$ symmetry offers a minimal framework for realizing electroweak baryogenesis (EWBG) during a first-order electroweak phase transition. In this work, we explore a novel mechanism where spontaneous $\mathbb{Z}_2$ breaking enables EWBG via domain walls separating two distinct phases of the $S$ field. These domain walls feature restored (or weakly broken) EW symmetry in their cores and sweep through space, generating the baryon asymmetry below the temperature of EW symmetry breaking. We identify the key conditions for the existence of EW-symmetric domain wall cores and analyze the dynamics required for wall propagation over sufficient spatial volumes. Additionally, we outline the CP-violating sources necessary for baryogenesis under different regimes of domain wall evolution. The parameter space accommodating this mechanism spans singlet masses from sub-eV to 15 GeV, accompanied by a non-vanishing mixing with the Higgs boson. Unlike the standard realization of EWBG in the minimal singlet-extended SM, which is notoriously difficult to test, our scenario can be probed by a wide range of existing and upcoming experiments, including fifth force searches, rare meson decays, and EDM measurements.

Our current understanding is that intermediate- to high-mass stars form in a way similar to low-mass stars, that is, through disk accretion. However, the physical conditions that play a role in disk formation, evolution, and the possibility of (sub)stellar companion formation, are significantly different. We search for the mm counterparts of four intermediate- to high-mass (4-10 Solar mass) young stellar objects (YSOs) in the giant Hii region M17 at a distance of 1.7 kpc. These objects expose their photospheric spectrum such that their location on the pre-main-sequence (PMS) is well established. They have a circumstellar disk that is likely remnant of the formation process. With ALMA we have detected, for the first time, these four YSOs in M17, in Band 6 and 7, as well as four other serendipitous objects. Besides the flux measurements, the source size and spectral index provide important constraints on the physical mechanism(s) producing the observed emission. We apply different models to estimate the dust and gas mass contained in the disks. All our detections are spatially unresolved, constraining the source size to <120 au, and have a spectral index in the range 0.5-2.7. The derived (upper limits on the) disk dust masses are on the order of a few Earth masses and estimations of the upper limits on the gas mass vary between $10^{-5}$ and $10^{-3}$ Solar mass. In two objects (B331 and B268) free-free emission indicates the presence of ionized material around the star. The four serendipitous detections are likely (low-mass) YSOs. We compare the derived disk masses of our M17 targets to those obtained for YSOs in low-mass star-forming regions (SFRs) and Herbig stars, as a function of stellar mass, age, luminosity, and outer disk radius. The M17 sample, though small, is both the most massive and the youngest sample, yet has the lowest mean disk mass.

Star-galaxy separation is a crucial step in creating target catalogues for extragalactic spectroscopic surveys. A classifier biased towards inclusivity risks including high numbers of stars, wasting fibre hours, while a more conservative classifier might overlook galaxies, compromising completeness and hence survey objectives. To avoid bias introduced by a training set in supervised methods, we employ an unsupervised machine learning approach. Using photometry from the Wide Area VISTA Extragalactic Survey (WAVES)-Wide catalogue comprising nine-band <inline-formula><tex-math id="TM0001" notation="LaTeX">$u - K_s$</tex-math></inline-formula> data, we create a feature space with colours, fluxes, and apparent size information extracted by PROFOUND. We apply the non-linear dimensionality reduction method UMAP (Uniform Manifold Approximation and Projection) combined with the classifier HDBSCAN (Hierarchical Density-Based Spatial Clustering of Applications with Noise) to classify stars and galaxies. Our method is verified against a baseline colour and morphological method using a truth catalogue from Gaia, SDSS (Sloan Digital Sky Survey), GAMA (Galaxy And Mass Assembly), and DESI (Dark Energy Spectroscopic Instrument). We correctly identify 99.75 per cent of galaxies within the AB magnitude limit of <inline-formula><tex-math id="TM0002" notation="LaTeX">$Z=21.2$</tex-math></inline-formula>, with an F1 score of <inline-formula><tex-math id="TM0003" notation="LaTeX">$0.9971 \pm 0.0018$</tex-math></inline-formula> across the entire ground truth sample, compared to <inline-formula><tex-math id="TM0004" notation="LaTeX">$0.9879 \pm 0.0088$</tex-math></inline-formula> from the baseline method. Our method's higher purity (<inline-formula><tex-math id="TM0005" notation="LaTeX">$0.9967 \pm 0.0021$</tex-math></inline-formula>) compared to the baseline (<inline-formula><tex-math id="TM0006" notation="LaTeX">$0.9795 \pm 0.0172$</tex-math></inline-formula>) increases efficiency, identifying 11 per cent fewer galaxy or ambiguous sources, saving approximately 70 000 fibre hours on the 4MOST (4-m Multi-Object Spectroscopic Telescope) instrument. We achieve reliable classification statistics for challenging sources including quasars, compact galaxies, and low surface brightness galaxies, retrieving 92.7 per cent, 84.6 per cent, and 99.5 per cent of them, respectively. Angular clustering analysis validates our classifications, showing consistency with expected galaxy clustering, regardless of the baseline classification.

Aims. We introduce a novel sub-resolution prescription to correct for the unresolved dynamical friction (DF) onto black holes (BHs) in cosmological simulations, to describe BH dynamics accurately, and to overcome spurious motions induced by numerical effects. Methods. We implemented a sub-resolution prescription for the unresolved DF onto BHs in the OpenGadget3 code. We carried out cosmological simulations of a volume of (16 comoving Mpc)³ and zoomed-in simulations of a galaxy group and of a galaxy cluster. We assessed the advantages of our new technique in comparison to commonly adopted methods for hampering spurious BH displacements, namely repositioning onto a local minimum of the gravitational potential and ad hoc boosting of the BH particle dynamical mass. We inspected variations in BH demography in terms of offset from the centres of the host sub-halos, the wandering population of BHs, BH–BH merger rates, and the occupation fraction of sub-halos. We also analysed the impact of the different prescriptions on individual BH interaction events in detail. Results. The newly introduced DF correction enhances the centring of BHs on host halos, the effects of which are at least comparable with those of alternative techniques. Also, the correction becomes gradually more effective as the redshift decreases. Simulations with this correction predict half as many merger events with respect to the repositioning prescription, with the advantage of being less prone to leaving substructures without any central BH. Simulations featuring our DF prescription produce a smaller (by up to ~50% with respect to repositioning) population of wandering BHs and final BH masses that are in good agreement with observations. Regarding individual BH–BH interactions, our DF model captures the gradual inspiraling of orbits before the merger occurs. By contrast, the repositioning scheme, in its most classical renditions, describes extremely fast mergers, while the dynamical mass misrepresents the dynamics of the black holes, introducing numerical scattering between the orbiting BHs. Conclusions. The novel DF correction improves the accuracy if tracking BHs within their hosts galaxies and the pathway to BH- BH mergers. This opens up new possibilities for better modeling the evolution of BH populations in cosmological simulations across different times and different environments.

Within $Z^\prime$ models, neutral meson mixing severely constrains beyond the Standard Model (SM) effects in flavour changing neutral current (FCNC) processes. However, in certain regions of the $Z^\prime$ parameter space, the contributions to meson mixing observables become negligibly small even for large $Z^\prime$ couplings. While this a priori allows for significant new physics (NP) effects in FCNC decays, we discuss how large $Z^\prime$ couplings in one neutral meson sector can generate effects in meson mixing observables of other neutral mesons, through correlations stemming from $\text{SU(2)}_L$ gauge invariance and through Renormalization Group (RG) effects in the SM Effective Field Theory~(SMEFT). This is illustrated with the example of $B_s^0-\bar B_s^0$ mixing, which in the presence of both left- and right-handed $Z^\prime bs$ couplings $\Delta_L^{bs}$ and $\Delta_R^{bs}$ remains SM-like for $\Delta_R^{bs}\approx 0.1\,\Delta_L^{bs}$. We show that in this case, large $Z^\prime bs$ couplings generate effects in $D$ and $K$ meson mixing observables, but that the $D$ and $K$ mixing constraints and the relation between $\Delta_R^{bs}$ and $\Delta_L^{bs}$ are fully compatible with a lepton flavour universality~(LFU) conserving explanation of the most recent $b\to s\ell^+\ell^-$ experimental data without violating other constraints like $e^+ e^-\to\ell^+\ell^-$ scattering. Assuming LFU, invariance under the $\text{SU(2)}_L$ gauge symmetry leads then to correlated effects in $b\to s\nu\bar\nu$ observables presently studied intensively by the Belle~II experiment, which allow to probe the $Z^\prime$ parameter space that is opened up by the vanishing NP contributions to $B_s^0-\bar B_s^0$ mixing. In this scenario the suppression of $B\to K(K^*)\mu^+\mu^-$ branching ratios implies {\em uniquely} enhancements of $B\to K(K^*)\nu\bar\nu$ branching ratios up to $20\%$.

Stable 205Tl ions have the lowest known energy threshold for capturing electron neutrinos (ve) of Eve≥50.6  keV. The Lorandite Experiment (LOREX), proposed in the 1980s, aims at obtaining the longtime averaged solar neutrino flux by utilizing natural deposits of Tl-bearing lorandite ores. To determine the ve capture cross section, it is required to know the strength of the weak transition connecting the ground state of 205Tl and the 2.3 keV first excited state in 205Pb. The only way to experimentally address this transition is to measure the bound-state beta decay (ßb) of fully ionized 205Tl81+ ions. After three decades of meticulous preparation, the half-life of the ßb decay of 205Tl81+ has been measured to be 291−27+33   days using the Experimental Storage Ring (ESR) at GSI, Darmstadt. The longer measured half-life compared to theoretical estimates reduces the expected signal-to-noise ratio in the LOREX, thus challenging its feasibility.

We discuss the relation between the Koonin-Pratt femtoscopic correlation function (CF) and invariant mass distributions from production experiments. We show that the equivalence is total for a zero source-size and that a Gaussian finite-size source provides a form-factor for the virtual production of the particles. Motivated by this remarkable relationship, we study an alternative method to the Koonin-Pratt formula, which connects the evaluation of the CF directly with the production mechanisms. The differences arise mostly from the <inline-formula><mml:math display="inline"><mml:mi>T</mml:mi></mml:math></inline-formula>–matrix quadratic terms and increase with the source size. We study the case of the <inline-formula><mml:math display="inline"><mml:msup><mml:mi>D</mml:mi><mml:mn>0</mml:mn></mml:msup><mml:msup><mml:mi>D</mml:mi><mml:mrow><mml:mo>*</mml:mo><mml:mo>+</mml:mo></mml:mrow></mml:msup></mml:math></inline-formula> and <inline-formula><mml:math display="inline"><mml:msup><mml:mi>D</mml:mi><mml:mo>+</mml:mo></mml:msup><mml:msup><mml:mi>D</mml:mi><mml:mrow><mml:mo>*</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula> correlation functions of interest to unravel the dynamics of the exotic <inline-formula><mml:math display="inline"><mml:msub><mml:mi>T</mml:mi><mml:mrow><mml:mi>c</mml:mi><mml:mi>c</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mn>3875</mml:mn><mml:msup><mml:mo stretchy="false">)</mml:mo><mml:mo>+</mml:mo></mml:msup></mml:math></inline-formula>, and find that these differences become quite sizable already for 1 fm sources. We nevertheless conclude that the lack of coherence in high-multiplicity-event reactions and in the creation of the fire-ball source that emits the hadrons certainly make much more realistic the formalism based on the Koonin-Pratt equation. We finally derive an improved Lednicky-Lyuboshits (LL) approach, which implements a Lorentz ultraviolet regulator that corrects the pathological behavior of the LL CF in the punctual source-size limit.

Inspired by earlier results on recursions for open-string tree-level amplitudes, and by a result of Brown and Dupont relating open- and closed-string tree-level amplitudes via single-valued periods, we identify a recursive relation for closed-string tree-level amplitudes. We achieve this by showing that closed-string analogues of Selberg integrals satisfy the Knizhnik-Zamolodchikov equation for a suitable matrix representation of the free Lie algebra on two generators, and by identifying the limits at z=1 and z=0, which are related by the Deligne associator, with N-point and (N-1)-point closed-string amplitudes, respectively.

A quantitative spectroscopic study of blue supergiant stars in the Hubble constant anchor galaxy NGC 4258 is presented. The non-LTE analysis of Keck I telescope LRIS spectra yields a central logarithmic metallicity (in units of the solar value) of [Z] = ‑0.05 ± 0.05 and a very shallow gradient of ‑(0.09 ± 0.11) r/r ₂₅ with respect to galactocentric distance in units of the isophotal radius. Good agreement with the mass–metallicity relationship of star-forming galaxies based on stellar absorption line studies is found. A comparison with H II region oxygen abundances obtained from the analysis of strong emission lines shows reasonable agreement when the M. Pettini & B. E. J. Pagel calibration is used, while the D. Zaritsky et al. calibration yields values that are 0.2–0.3 dex larger. These results allow us to put the metallicity calibration of the Cepheid period–luminosity relation in this anchor galaxy on a purely stellar basis. Interstellar reddening and extinction are determined using Hubble Space Telescope and JWST photometry. Based on extinction-corrected magnitudes, combined with the stellar effective temperatures and gravities we determine, we use the flux-weighted gravity–luminosity relationship to estimate an independent spectroscopic distance. We obtain a distance modulus m ‑ M = 29.38 ± 0.12 mag, in agreement with the geometrical distance derived from the analysis of the water maser orbits in the galaxy's central circumnuclear disk.

HD 54879 is the most recently discovered magnetic O-type star. Previous studies ruled out a rotation period shorter than 7 yr, implying that HD 54879 is the second most slowly rotating known magnetic O-type star. We report new high-resolution spectropolarimetric measurements of HD 54879, which confirm that a full stellar rotation cycle has been observed. We derive a stellar rotation period from the longitudinal magnetic field measurements of <inline-formula> </inline-formula> days (about 7.02 yr). The radial velocity of HD 54879 has been stable over the last decade of observations. We explore equivalent widths and longitudinal magnetic fields calculated from lines of different elements, and conclude the atmosphere of HD 54879 is likely chemically homogeneous, with no strong evidence for chemical stratification or lateral abundance nonuniformities. We present the first detailed magnetic map of the star, with an average surface-magnetic-field strength of 2954 G, and a strength for the dipole component of 3939 G. There is a significant amount of magnetic energy in the quadrupole components of the field (23%). Thus, we find HD 54879 has a strong magnetic field with a significantly complex topology.

Large angular scale surveys in the absence of atmosphere are essential for measuring the primordial B-mode power spectrum of the Cosmic Microwave Background (CMB). Since this proposed measurement is about three to four orders of magnitude fainter than the temperature anisotropies of the CMB, in-flight calibration of the instruments and active suppression of systematic effects are crucial. We investigate the effect of changing the parameters of the scanning strategy on the in-flight calibration effectiveness, the suppression of the systematic effects themselves, and the ability to distinguish systematic effects by null-tests. Next-generation missions such as LiteBIRD, modulated by a Half-Wave Plate (HWP), will be able to observe polarisation using a single detector, eliminating the need to combine several detectors to measure polarisation, as done in many previous experiments and hence avoiding the consequent systematic effects. While the HWP is expected to suppress many systematic effects, some of them will remain. We use an analytical approach to comprehensively address the mitigation of these systematic effects and identify the characteristics of scanning strategies that are the most effective for implementing a variety of calibration strategies in the multi-dimensional space of common spacecraft scan parameters. We verify that LiteBIRD's standard configuration yields good performance on the metrics we studied. We also present Falcons.jl, a fast spacecraft scanning simulator that we developed to investigate this scanning parameter space.

We present a novel multi-messenger approach for probing nonstandard neutrino properties through the detection of gravitational waves (GWs) from collapsing stellar cores and associated supernova explosions. We show that neutrino flavor conversion inside the proto-neutron star (PNS), motivated by physics Beyond the Standard Model (BSM), can significantly boost PNS convection. This effect leads to large-amplitude GW emission over a wide frequency range during an otherwise relatively quiescent GW phase shortly after core bounce. Such a signal provides a promising new avenue for exploring nonstandard neutrino phenomena and other BSM physics impacting PNS convection.

Even in the absence of neutrino masses, a neutrino gas can exhibit a homogeneous flavor instability that leads to a periodic motion known as the fast flavor pendulum. A well-known necessary condition is a crossing of the angular flavor lepton distribution. However, in contrast to an earlier finding by two of us, the Nyquist criterion inspired by plasma physics, while being a more restrictive necessary condition, is not always sufficient. The question depends on the unstable branch of the dispersion relation being bounded by critical points that both lie under the light cone (points with subluminal phase velocity), in which case the Nyquist criterion is sufficient. While the lepton-number angle distribution, assumed to be axially symmetric, easily allows one to determine the real branch of the dispersion relation and to recognize if instead superluminal critical points exist, this graphical method does not translate into a simple instability condition. We discuss the dispersion relation for the homogeneous mode in the more general context of modes with arbitrary wave number and stress that it plays no special role on this continuum, except for its regular but fragile long-term behavior owed to its many symmetries.

We present ultraviolet, optical, and near-infrared photometric and optical spectroscopic observations of the luminous fast blue optical transient (LFBOT) CSS 161010:045834–081803 (CSS 161010). The transient was found in a low-redshift (z = 0.033) dwarf galaxy. The light curves of CSS 161010 are characterized by an extremely fast evolution and blue colors. The V-band light curve shows that CSS 161010 reaches an absolute peak of <inline-formula> </inline-formula> mag in 3.8 days from the start of the outburst. After maximum, CSS 161010 follows a power-law decline ∝t ^‑2.8±0.1 in all optical bands. These photometric properties are comparable to those of well-observed LFBOTs such as AT 2018cow, AT 2020mrf, and AT 2020xnd. However, unlike these objects, the spectra of CSS 161010 show a remarkable transformation from a blue and featureless continuum to spectra dominated by very broad, entirely blueshifted hydrogen emission lines with velocities of up to 10% of the speed of light. The persistent blueshifted emission and the lack of any emission at the rest wavelength of CSS 161010 are unique features not seen in any transient before CSS 161010. The combined observational properties of CSS 161010 and its M _* ∼ 10⁸ M _⊙ dwarf galaxy host favor the tidal disruption of a star by an intermediate-mass black hole as its origin.

We assess the capabilities of the LiteBIRD mission to map the hot gas distribution in the Universe through the thermal Sunyaev-Zeldovich (SZ) effect. Our analysis relies on comprehensive simulations incorporating various sources of Galactic and extragalactic foreground emission, while accounting for the specific instrumental characteristics of the LiteBIRD mission, such as detector sensitivities, frequency-dependent beam convolution, inhomogeneous sky scanning, and 1/f noise. We implement a tailored component-separation pipeline to map the thermal SZ Compton y-parameter over 98% of the sky. Despite lower angular resolution for galaxy cluster science, LiteBIRD provides full-sky coverage and, compared to the Planck satellite, enhanced sensitivity, as well as more frequency bands to enable the construction of an all-sky thermal SZ y-map, with reduced foreground contamination at large and intermediate angular scales. By combining LiteBIRD and Planck channels in the component-separation pipeline, we also obtain an optimal y-map that leverages the advantages of both experiments, with the higher angular resolution of the Planck channels enabling the recovery of compact clusters beyond the LiteBIRD beam limitations, and the numerous sensitive LiteBIRD channels further mitigating foregrounds. The added value of LiteBIRD is highlighted through the examination of maps, power spectra, and one-point statistics of the various sky components. After component separation, the 1/f noise from LiteBIRD's intensity channels is effectively mitigated below the level of the thermal SZ signal at all multipoles. Cosmological constraints on S ₈ = σ ₈ (Ω_m /0.3)^0.5 obtained from the LiteBIRD-Planck combined y-map power spectrum exhibits a 15 % reduction in uncertainty compared to constraints derived from Planck alone. This improvement can be attributed to the increased portion of uncontaminated sky available in the LiteBIRD-Planck combined y-map.

The goal of this paper is to make a connection between tropical geometry, representations of quantum affine algebras, and scattering amplitudes in physics. The connection allows us to study important and difficult questions in these areas: (1)We give a systematic construction of prime modules (including prime non-real modules) of quantum affine algebras using tropical geometry. We also introduce new objects which generalize positive tropical Grassmannians.(2)We propose a generalization of Grassmannian string integrals in physics, in which the integrand is a product indexed by prime modules of a quantum affine algebra. We give a general formula of u-variables using prime tableaux (corresponding to prime modules of quantum affine algebras of type A) and Auslander–Reiten quivers of Grassmannian cluster categories.(3)We study limit g-vectors of cluster algebras. This is another way to obtain prime non-real modules of quantum affine algebras systematically. Using limit g-vectors, we construct new examples of non-real modules of quantum affine algebras. We give a systematic construction of prime modules (including prime non-real modules) of quantum affine algebras using tropical geometry. We also introduce new objects which generalize positive tropical Grassmannians. We propose a generalization of Grassmannian string integrals in physics, in which the integrand is a product indexed by prime modules of a quantum affine algebra. We give a general formula of u-variables using prime tableaux (corresponding to prime modules of quantum affine algebras of type A) and Auslander–Reiten quivers of Grassmannian cluster categories. We study limit g-vectors of cluster algebras. This is another way to obtain prime non-real modules of quantum affine algebras systematically. Using limit g-vectors, we construct new examples of non-real modules of quantum affine algebras.

This report summarizes the work of the EMMI Rapid Reaction Task Force on "Real and Virtual Photon Production at Ultra-Low Transverse Momentum and Low Mass at the LHC". We provide an overview of the soft-photon puzzle, i.e., of the long-standing discrepancy between experimental data and predictions based on Low's soft-photon theorem, also referred to as "anomalous" soft photon production, and we review the current theoretical understanding of soft radiation and soft theorems. We also focus on low-mass dileptons as a tool for determining the electrical conductivity of the medium produced in high-energy nucleus–nucleus collisions. We discuss how both topics can be addressed with the planned ALICE 3 detector at the LHC.

The stellar haloes of dwarf galaxies are becoming an object of interest in the extragalactic community due to their detection in some recent observations. Additionally, new cosmological simulations of very high resolution were performed, allowing their study. These stellar haloes could help shed light on our understanding of the assembly of dwarf galaxies and their evolution, and allow us to test the hierarchical model for the formation of structures at small scales. We aim to characterise the stellar haloes of simulated dwarf galaxies and analyse their evolution and accretion history. We use a sample of 17 simulated galaxies from the Auriga Project with a stellar mass range from 3.28x10^8 Msun to 2.08x10^10 Msun. We define the stellar halo as the stellar material located outside an ellipsoid with semi-major axes equal to 4 times the half light radius (Rh) of each galaxy. We find that the inner regions of the stellar halo (4 to 6 times the Rh) are dominated by in-situ material. For the less massive simulated dwarfs (M*<=4.54x10^8 Msun), this dominance extends to all radii. We find that this in-situ stellar halo is mostly formed in the inner regions of the galaxies and then ejected into the outskirts during interactions and merger events. In ~50% of the galaxies, the stripped gas from satellites contributed to the formation of this in-situ halo. The stellar haloes of the galaxies more massive than M*>=1x10^9 Msun are dominated by the accreted component beyond 6 Rh. We find that the more massive dwarf galaxies accrete stellar material until later times (t90~4.44 Gyr ago, being t90 the formation time) than the less massive ones (t90~8.17 Gyr ago), impacting on the formation time of the accreted stellar haloes. The galaxies have a range of 1 to 7 significant progenitors contributing to their accreted component but there is no correlation between this quantity and the galaxies' accreted mass.

A novel coalescence process is shown to take place in plasma fluid simulations, leading to the formation of large-scale magnetic islands that become dynamically important in the system. The parametric dependence of the process on the plasma $\beta$ and the background magnetic shear is studied, and the process is broken down at a fundamental level, allowing to clearly identify its causes and dynamics. The formation of magnetic-island-like structures at the spatial scale of the unstable modes is observed quite early in the non-linear phase of the simulation for most cases studied, as the unstable modes change their structure from interchange-like to tearing-like. This is followed by a slow coalescence process that evolves these magnetic structures towards larger and larger scales, adding to the large-scale tearing-like modes that already form by direct coupling of neighbouring unstable modes, but remain sub-dominant without the contribution from the smaller scales through coalescence. The presence of the cubic non-linearities retained in the model is essential in the dynamics of this process. The zonal fields are key actors of the overall process, acting as mediators between the competitive mechanisms from which Turbulence Driven Magnetic Islands can develop. The zonal current is found to slow down the formation of large-scale magnetic islands, acting as an inhibitor, while the zonal flow is needed to allow the system to transfer energy to the larger scales, acting as a catalyst for the island formation process.

Brown dwarfs are the bridge between low-mass stars and giant planets. One way of shedding light on their dominant formation mechanism is to study them at the earliest stages of their evolution, when they are deeply embedded in their parental clouds. Several works have identified pre- and proto-brown dwarfs candidates using different observational approaches. The aim of this work is to create a database with all the objects classified as very young substellar candidates in the litearature in order to study them in an homogeneous way. We have gathered all the information about very young substellar candidates available in the literature until 2020. We have retrieved their published photometry from the optical to the centimeter regime, and we have written our own codes to derive their bolometric temperatures and luminosities, and their internal luminosities. We have also populated the database with other parameters extracted from the literature, like e.g. the envelope masses, their detection in some molecular species, and presence of outflows. The result of our search is the SUCANES database, containing 174 objects classified as potential very young substellar candidates in the literature. We present an analysis of the main properties of the retrieved objects. Since we have updated the distances to several star forming regions, this has allowed us to reject some candidates based on their internal luminosities. We have also discussed the derived physical parameters and envelope masses for the best substellar candidates isolated in SUCANES. As an example of a scientific exploitation of this database, we present a feasibility study for the detection of radiojets with upcoming facilities: the ngVLA and the SKA interferometers. The SUCANES database is accessible through a Graphical User Interface and it is open to any potential user.

While supervised neural networks have become state of the art for identifying the rare strong gravitational lenses from large imaging data sets, their selection remains significantly affected by the large number and diversity of non-lens contaminants. This work evaluates and compares systematically the performance of neural networks in order to move towards a rapid selection of galaxy-scale strong lenses with minimal human input in the era of deep, wide-scale surveys. We used multiband images from PDR2 of the Hyper-Suprime Cam (HSC) Wide survey to build test sets mimicking an actual classification experiment, with 189 securely-identified strong lenses from the literature over the HSC footprint and 70 910 non-lens galaxies in COSMOS covering representative lens-like morphologies. Multiple networks were trained on different sets of realistic strong-lens simulations and non-lens galaxies, with various architectures and data preprocessing, mainly using the deepest gri-bands. Most networks reached excellent area under the Receiver Operating Characteristic (ROC) curves on the test set of 71 099 objects, and we determined the ingredients to optimize the true positive rate for a total number of false positives equal to zero or 10 (TPR₀ and TPR₁₀). The overall performances strongly depend on the construction of the ground-truth training data and they typically, but not systematically, improve using our baseline residual network architecture presented in Paper VI (Cañameras et al., A&A, 653, L6). TPR₀ tends to be higher for ResNets (≃ 10–40%) compared to AlexNet-like networks or G-CNNs. Improvements are found when (1) applying random shifts to the image centroids, (2) using square-root scaled images to enhance faint arcs, (3) adding z-band to the otherwise used gri-bands, or (4) using random viewpoints of the original images. In contrast, we find no improvement when adding g – αi difference images (where α is a tuned constant) to subtract emission from the central galaxy. The most significant gain is obtained with committees of networks trained on different data sets, with a moderate overlap between populations of false positives. Nearly-perfect invariance to image quality can be achieved by using realistic PSF models in our lens simulation pipeline, and by training networks either with large number of bands, or jointly with the PSF and science frames. Overall, we show the possibility to reach a TPR₀ as high as 60% for the test sets under consideration, which opens promising perspectives for pure selection of strong lenses without human input using the Rubin Observatory and other forthcoming ground-based surveys.

Metals in the diffuse, ionized gas at the boundary between the Milky Way's interstellar medium (ISM) and circumgalactic medium, known as the disk–halo interface (DHI), are valuable tracers of the feedback processes that drive the Galactic fountain. However, metallicity measurements in this region are challenging due to obscuration by the Milky Way ISM and uncertain ionization corrections that affect the total hydrogen column density. In this work, we constrain ionization corrections to neutral hydrogen column densities using precisely measured electron column densities from the dispersion measures of pulsars that lie in the same globular clusters as UV-bright targets with high-resolution absorption spectroscopy. We address the blending of absorption lines with the ISM by jointly fitting Voigt profiles to all absorption components. We present our metallicity estimates for the DHI of the Milky Way based on detailed photoionization modeling of the absorption from ionized metal lines and ionization-corrected total hydrogen columns. Generally, the gas clouds show a large scatter in metallicity, ranging between 0.04 and 3.2 Z _⊙, implying that the DHI consists of a mixture of gaseous structures having multiple origins. We estimate the inflow and outflow timescales of the DHI ionized clouds to be 6–35 Myr. We report the detection of an infalling cloud with supersolar metallicity that suggests a Galactic fountain mechanism, whereas at least one low-metallicity outflowing cloud (Z < 0.1 Z _⊙) poses a challenge for Galactic fountain and feedback models.

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

Multi-messenger observations of astrophysical transients provide powerful probes of the underlying physics of the source as well as beyond the Standard Model effects. We explore transients that can occur in the vicinity of supermassive black holes at the center of galaxies, including tidal disruption events (TDEs), certain types of blazars, or even supernovae. In such environments, the dark matter (DM) density can be extremely high, resembling a dense spike or core. We study a novel effect of neutrino diffusion sustained via frequent scatterings off DM particles in these regions. We show that for transients occurring within DM spikes or cores, the DM-neutrino scattering can delay the arrival of neutrinos with respect to photons, but this also comes with a suppression of the neutrino flux and energy loss. We apply these effects to the specific example of TDEs, and demonstrate that currently unconstrained parameter space of DM-neutrino interactions can account for the sizable $O$(days) delay of the tentative high-energy neutrinos observed from some TDEs.

Cosmic shear, galaxy clustering, and the abundance of massive halos each probe the large-scale structure of the universe in complementary ways. We present cosmological constraints from the joint analysis of the three probes, building on the latest analyses of the lensing-informed abundance of clusters identified by the South Pole Telescope (SPT) and of the auto- and cross-correlation of galaxy position and weak lensing measurements (3$\times$2pt) in the Dark Energy Survey (DES). We consider the cosmological correlation between the different tracers and we account for the systematic uncertainties that are shared between the large-scale lensing correlation functions and the small-scale lensing-based cluster mass calibration. Marginalized over the remaining $\Lambda$CDM parameters (including the sum of neutrino masses) and 52 astrophysical modeling parameters, we measure $\Omega_\mathrm{m}=0.300\pm0.017$ and $\sigma_8=0.797\pm0.026$. Compared to constraints from Planck primary CMB anisotropies, our constraints are only 15% wider with a probability to exceed of 0.22 ($1.2\sigma$) for the two-parameter difference. We further obtain $S_8\equiv\sigma_8(\Omega_\mathrm{m}/0.3)^{0.5}=0.796\pm0.013$ which is lower than the Planck measurement at the $1.6\sigma$ level. The combined SPT cluster, DES 3$\times$2pt, and Planck datasets mildly prefer a non-zero positive neutrino mass, with a 95% upper limit $\sum m_\nu<0.25~\mathrm{eV}$ on the sum of neutrino masses. Assuming a $w$CDM model, we constrain the dark energy equation of state parameter $w=-1.15^{+0.23}_{-0.17}$ and when combining with Planck primary CMB anisotropies, we recover $w=-1.20^{+0.15}_{-0.09}$, a $1.7\sigma$ difference with a cosmological constant. The precision of our results highlights the benefits of multiwavelength multiprobe cosmology.

We present the Astrophysical Multimessenger Modeling (AM ³ ) software. AM ³ is a documented open-source software (source code at gitlab.desy.de/am3/am3; user guide and documentation at am3.readthedocs.io/en/latest/) that efficiently solves the coupled integro-differential equations describing the temporal evolution of the spectral densities of particles interacting in astrophysical environments, including photons, electrons, positrons, protons, neutrons, pions, muons, and neutrinos. The software has been extensively used to simulate the multiwavelength and neutrino emission from active galactic nuclei (including blazars), gamma-ray bursts, and tidal disruption events. The simulations include all relevant nonthermal processes, namely synchrotron emission, inverse Compton scattering, photon–photon annihilation, proton–proton and proton–photon pion production, and photo-pair production. The software self-consistently calculates the full cascade of primary and secondary particles, including nonlinear feedback processes and predictions in the time domain. It also allows the user to track separately the particle densities produced by means of each distinct interaction process, including the different hadronic channels. With its efficient hybrid solver combining analytical and numerical techniques, AM ³ combines efficiency and accuracy at a user-adjustable level. We describe the technical details of the numerical framework and present three examples of applications to different astrophysical environments.

We extend the evolution-mapping approach, introduced in the first paper of this series to describe non-linear matter density fluctuations, to statistics of the cosmic velocity field. This framework classifies cosmological parameters into shape parameters, which determine the shape of the linear matter power spectrum, <inline-formula><tex-math id="TM0001" notation="LaTeX">$P_{\rm L}(k, z)$</tex-math></inline-formula>, and evolution parameters, which control its amplitude at any redshift. Evolution-mapping leverages the fact that density fluctuations in cosmologies with identical shape parameters but different evolution parameters exhibit similar non-linear evolutions when expressed as a function of clustering amplitude. We analyse a suite of N-body simulations sharing identical shape parameters but spanning a wide range of evolution parameters. Using a method for estimating the volume-weighted velocity field based on the Voronoi tessellation of simulation particles, we study the non-linear evolution of the velocity divergence power spectrum, <inline-formula><tex-math id="TM0002" notation="LaTeX">$P_{\theta \theta }(k)$</tex-math></inline-formula>, and its cross-power spectrum with the density field, <inline-formula><tex-math id="TM0003" notation="LaTeX">$P_{\delta \theta }(k)$</tex-math></inline-formula>. We demonstrate that the evolution-mapping relation applies accurately to <inline-formula><tex-math id="TM0004" notation="LaTeX">$P_{\theta \theta }(k)$</tex-math></inline-formula> and <inline-formula><tex-math id="TM0005" notation="LaTeX">$P_{\delta \theta }(k)$</tex-math></inline-formula>. While this breaks down in the strongly non-linear regime, deviations can be modelled in terms of differences in the suppression factor, <inline-formula><tex-math id="TM0006" notation="LaTeX">$g(a) = D(a)/a$</tex-math></inline-formula>, similar to those for the density field. Such modelling describes the differences in <inline-formula><tex-math id="TM0007" notation="LaTeX">$P_{\theta \theta }(k)$</tex-math></inline-formula> between models with the same linear clustering amplitude to better than 1 per cent accuracy at all scales and redshifts considered. Evolution-mapping simplifies the description of the cosmological dependence of non-linear density and velocity statistics, streamlining the sampling of large cosmological parameter spaces for cosmological analysis.

Very-metal-poor stars ([Fe/H] < ‑2) are important laboratories for testing stellar models and reconstructing the formation history of our galaxy. Asteroseismology is a powerful tool to probe stellar interiors and measure ages, but few asteroseismic detections are known in very-metal-poor stars and none have allowed detailed modeling of oscillation frequencies. We report the discovery of a low-luminosity Kepler red giant (KIC 8144907) with high signal-to-noise ratio oscillations, [Fe/H] = ‑2.66 ± 0.08 and [α/Fe] = 0.38 ± 0.06, making it by far the most metal-poor star to date for which detailed asteroseismic modeling is possible. By combining the oscillation spectrum from Kepler with high-resolution spectroscopy, we measure an asteroseismic mass and age of 0.79 ± 0.02(ran) ± 0.01(sys) M _⊙ and 12.0 ± 0.6(ran) ± 0.4(sys) Gyr, with remarkable agreement across different codes and input physics, demonstrating that stellar models and asteroseismology are reliable for very-metal-poor stars when individual frequencies are used. The results also provide a direct age anchor for the early formation of the Milky Way, implying that substantial star formation did not commence until redshift z ≈ 3 (if the star formed in situ) or that the Milky Way has undergone merger events for at least ≈12 Gyr (if the star was accreted by a dwarf satellite merger such as Gaia-Enceladus).

Strongly gravitationally lensed supernovae (LSNe) are promising probes for providing absolute distance measurements using gravitational-lens time delays. Spatially unresolved LSNe offer an opportunity to enhance the sample size for precision cosmology. We predict that there will be approximately three times as many unresolved as resolved LSNe Ia in the Legacy Survey of Space and Time (LSST) by the Rubin Observatory. In this article, we explore the feasibility of detecting unresolved LSNe Ia from a pool of preclassified SNe Ia light curves using the shape of the blended light curves with deep-learning techniques. We find that ∼30% unresolved LSNe Ia can be detected with a simple 1D convolutional neural network (CNN) using well-sampled rizy-band light curves (with a false-positive rate of ∼3%). Even when the light curve is well observed in only a single band among r, i, and z, detection is still possible with false-positive rates ranging from ∼4 to 7% depending on the band. Furthermore, we demonstrate that these unresolved cases can be detected at an early stage using light curves up to ∼20 days from the first observation with well-controlled false-positive rates, providing ample opportunity to trigger follow-up observations. Additionally, we demonstrate the feasibility of time-delay estimations using solely LSST-like data of unresolved light curves, particularly for doubles, when excluding systems with low time delays and magnification ratios. However, the abundance of such systems among those unresolved in LSST poses a significant challenge. This approach holds potential utility for upcoming wide-field surveys, and overall results could significantly improve with enhanced cadence and depth in the future surveys.

Context. A low-mass companion potentially in the brown dwarf mass regime was discovered on a ~12 yr orbit (~5.5 au) around HD 167665 using radial velocity (RV) monitoring. Joint RV–astrometry analyses confirmed that HD 167665B is a brown dwarf with precisions on the measured mass of ~4–9%. Brown dwarf companions with measured mass and luminosity are valuable for testing formation and evolutionary models. However, its atmospheric properties and luminosity are still unconstrained, preventing detailed tests of evolutionary models. Aims. We further characterize the HD 167665 system by measuring the luminosity and refining the mass of its companion and reassessing the stellar age. Methods. We present new high-contrast imaging data of the star and of its close-in environment from SPHERE and GRAVITY, which we combined with RV data from CORALIE and HIRES and astrometry from HIPPARCOS and Gaia. Results. The analysis of the host star properties indicates an age of 6.20 ± 1.13 Gyr. GRAVITY reveals a point source near the position predicted from a joint fit of RV data and HIPPARCOS–Gaia proper motion anomalies. Subsequent SPHERE imaging confirms the detection and reveals a faint point source of contrast of ∆H2 = 10.95 ± 0.33 mag at a projected angular separation of ~180 mas. A joint fit of the high-contrast imaging, RV, and HIPPARCOS intermediate astrometric data together with the Gaia astrometric parameters constrains the mass of HD 167665B to ~1.2%, 60.3 ± 0.7 M_J. The SPHERE colors and spectrum point to an early or mid-T brown dwarf of spectral type T4_‑2⁺¹. Fitting the SPHERE spectrophotometry and GRAVITY spectrum with synthetic spectra suggests an effective temperature of ~1000–1150 K, a surface gravity of ~5.0–5.4 dex, and a bolometric luminosity log(L/L_⊙)=‑4.892_‑0.028^+0.024 dex. The mass, luminosity, and age of the companion can only be reproduced within 3σ by the hybrid cloudy evolutionary models of Saumon & Marley (2008, ApJ, 689, 1327), whereas cloudless evolutionary models underpredict its luminosity.

The high-precision measurements of exoplanet transit light curves that are now available contain information about the planet properties, their orbital parameters, and stellar limb darkening (LD). Recent 3D magnetohydrodynamical (MHD) simulations of stellar atmospheres have shown that LD depends on the photospheric magnetic field, and hence its precise determination can be used to estimate the field strength. Among existing LD laws, the uses of the simplest ones may lead to biased inferences, whereas the uses of complex laws typically lead to a large degeneracy among the LD parameters. We have developed a novel approach in which we use a complex LD model but with second derivative regularization during the fitting process. Regularization controls the complexity of the model appropriately and reduces the degeneracy among LD parameters, thus resulting in precise inferences. The tests on simulated data suggest that our inferences are not only precise but also accurate. This technique is used to re-analyse 43 transit light curves measured by the NASA Kepler and Transiting Exoplanet Survey Satellite missions. Comparisons of our LD inferences with the corresponding literature values show good agreement, while the precisions of our measurements are better by up to a factor of 2. We find that 1D non-magnetic model atmospheres fail to reproduce the observations while 3D MHD simulations are qualitatively consistent. The LD measurements, together with MHD simulations, confirm that Kepler-17, WASP-18, and KELT-24 have relatively high magnetic fields (<inline-formula><tex-math id="TM0001" notation="LaTeX">$\gt 200$</tex-math></inline-formula> G). This study paves the way for estimating the stellar surface magnetic field using the LD measurements.

The energy released by active galactic nuclei (AGN) has the potential to heat or remove the gas of the ISM, thus likely impacting the cold molecular gas reservoir of host galaxies at first, with star formation following as a consequence on longer timescales. Previous works on high-z galaxies, which compared the gas content of those without identified AGN, have yielded conflicting results, possibly due to selection biases and other systematics. To provide a reliable benchmark for galaxy evolution models at cosmic noon (z = 1 ‑ 3), two surveys were conceived: SUPER and KASHz, both targeting unbiased X-ray-selected AGN at z > 1 that span a wide bolometric luminosity range. In this paper we assess the effects of AGN feedback on the molecular gas content of host galaxies in a statistically robust, uniformly selected, coherently analyzed sample of AGN at z = 1 ‑ 2.6, drawn from the KASHz and SUPER surveys. By using targeted and archival ALMA data in combination with dedicated SED modeling, we retrieve CO and far-infrared (FIR) luminosity as well as M_* of SUPER and KASHz host galaxies. We selected non-active galaxies from PHIBBS, ASPECS, and multiple ALMA/NOEMA surveys of submillimeter galaxies in the COSMOS, UDS, and ECDF fields. By matching the samples in redshift, stellar mass, and FIR luminosity, we compared the properties of AGN and non-active galaxies within a Bayesian framework. We find that AGN hosts at given FIR luminosity are on average CO depleted compared to non-active galaxies, thus confirming what was previously found in the SUPER survey. Moreover, the molecular gas fraction distributions of AGN and non-active galaxies are statistically different, with the distribution of AGN being skewed to lower values. Our results indicate that AGN can indeed reduce the total cold molecular gas reservoir of their host galaxies. Lastly, by comparing our results with predictions from three cosmological simulations (TNG, Eagle, and Simba) filtered to match the properties of observed AGN, AGN hosts, and non-active galaxies, we confirm already known discrepancies and highlight new discrepancies between observations and simulations.

Aims. Our goal is twofold. First, to detect new clusters we apply the newest methods for the detection of clustering with the best available wide-field sky surveys in the mid-infrared because they are the least affected by extinction. Second, we address the question of cluster detection's completeness, for now limiting it to the most massive star clusters. Methods. This search is based on the mid-infrared Galactic Legacy Infrared Mid Plane Survey Extraordinaire (GLIMPSE), to minimize the effect of dust extinction. The search Ordering Points To Identify the Clustering Structure (OPTICS) clustering algorithm is applied to identify clusters, after excluding the bluest, presumably foreground sources, to improve the cluster-to-field contrast. The success rate for cluster identification is estimated with a semi-empirical simulation that adds clusters, based on the real objects, to the point source catalog, to be recovered later with the same search algorithm that was used in the search for new cluster candidates. As a first step, this is limited to the most massive star clusters with a total mass of 104 $M_\odot$. Results. Our automated search, combined with inspection of the color-magnitude diagrams and images yielded 659 cluster candidates; 106 of these appear to have been previously identified, suggesting that a large hidden population of star clusters still exists in the inner Milky Way. However, the search for the simulated supermassive clusters achieves a recovery rate of 70 to 95%, depending on the distance and extinction toward them. Conclusions. The new candidates, if confirmed, indicate that the Milky Way still harbors a sizeable population of still unknown clusters. However, they must be objects of modest richness, because our simulation indicates that there is no substantial hidden population of supermassive clusters in the central region of our Galaxy.

Future cosmic microwave background (CMB) experiments are primarily targeting a detection of the primordial $B$-mode polarisation. The faintness of this signal requires exquisite control of systematic effects which may bias the measurements. In this work, we derive requirements on the relative calibration accuracy of the overall polarisation gain ($\Delta g_\nu$) for LiteBIRD experiment, through the application of the blind Needlet Internal Linear Combination (NILC) foreground-cleaning method. We find that minimum variance techniques, as NILC, are less affected by gain calibration uncertainties than a parametric approach, which requires a proper modelling of these instrumental effects. The tightest constraints are obtained for frequency channels where the CMB signal is relatively brighter (166 GHz channel, $\Delta {g}_\nu \approx 0.16 \%$), while, with a parametric approach, the strictest requirements were on foreground-dominated channels. We then propagate gain calibration uncertainties, corresponding to the derived requirements, into all frequency channels simultaneously. We find that the overall impact on the estimated $r$ is lower than the required budget for LiteBIRD by almost a factor $5$. The adopted procedure to derive requirements assumes a simple Galactic model. We therefore assess the robustness of obtained results against more realistic scenarios by injecting the gain calibration uncertainties, according to the requirements, into LiteBIRD simulated maps and assuming intermediate- and high-complexity sky models. In this case, we employ the so-called Multi-Clustering NILC (MC-NILC) foreground-cleaning pipeline and obtain that the impact of gain calibration uncertainties on $r$ is lower than the LiteBIRD gain systematics budget for the intermediate-complexity sky model. For the high-complexity case, instead, it would be necessary to tighten the requirements by a factor $1.8$.

Jet observables at hadron colliders feature ''super-leading'' logarithms, double-logarithmic corrections resulting from a breakdown of color coherence due to complex phases in hard-scattering amplitudes. While these effects only arise in high orders of perturbation theory and are suppressed in the large-$N_c$ limit, they formally constitute leading logarithmic corrections to the cross sections. We present the first analysis of the corresponding contributions to a hadronic cross section, including all partonic channels and interference effects. Interestingly, some interference terms in partonic $q\bar q\to q\bar q$ scattering are only linearly suppressed in $1/N_c$. Our results for the $pp\to 2$ jets gap-between-jets cross section demonstrate the numerical importance of super-leading logarithms for small values of the veto scale $Q_0$, showing that these contributions should be accounted for in precision studies of such observables.

We present a machine learning search for high-redshift (5.0 < z < 6.5) quasars using the combined photometric data from the Dark Energy Spectroscopic Instrument (DESI) Imaging Legacy Surveys and the Wide-field Infrared Survey Explorer survey. We explore the imputation of missing values for high-redshift quasars, discuss the feature selections, compare different machine learning algorithms, and investigate the selections of class ensemble for the training sample, then we find that the random forest model is very effective in separating the high-redshift quasars from various contaminators. The 11 class random forest model can achieve a precision of 96.43% and a recall of 91.53% for high-redshift quasars for the test set. We demonstrate that the completeness of the high-redshift quasars can reach as high as 82.20%. The final catalog consists of 216,949 high-redshift quasar candidates with 476 high probable ones in the entire Legacy Surveys DR9 footprint, and we make the catalog publicly available. Using Multi Unit Spectroscopic Explorer (MUSE) and DESI early data release (EDR) public spectra, we find that 14 true high-redshift quasars (11 in the training sample) out of 21 candidates are correctly identified for MUSE, and 20 true high-redshift quasars (11 in the training sample) out of 21 candidates are correctly identified for DESI-EDR. Additionally, we estimate photometric redshift for the high-redshift quasar candidates using a random forest regression model with a high precision.

In this work we revisit the problem of the dynamical stability of hierarchical triple systems with applications to circumbinary planetary orbits. We derive critical semimajor axes based on simulating and analyzing the dynamical behavior of 3 × 10⁸ binary star–planet configurations. For the first time, three-dimensional and eccentric planetary orbits are considered. We explore systems with a variety of binary and planetary mass ratios, binary and planetary eccentricities from 0 to 0.9, and orbital mutual inclinations ranging from 0° to 180°. Planetary masses range between the size of Mercury and the lower fusion boundary (approximately 13 Jupiter masses). The stability of each system is monitored over 10⁶ planetary orbital periods. We provide empirical expressions in the form of multidimensional, parameterized fits for two borders that separate dynamically stable, unstable, and mixed zones. In addition, we offer a machine learning model trained on our data set as an alternative tool for predicting the stability of circumbinary planets. Both the empirical fits and the machine learning model are tested for their predictive capabilities against randomly generated circumbinary systems with very good results. The empirical formulae are also applied to the Kepler and TESS circumbinary systems, confirming that many planets orbit their host stars close to the stability limit of those systems. Finally, we present a REST application programming interface with a web-based application for convenient access to our simulation data set.

Context. The existence of axion quark nuggets is a potential consequence of the axion field, which provides a possible solution to the charge-conjugation parity violation in quantum chromodynamics. In addition to explaining the cosmological discrepancy of matter-antimatter asymmetry and a visible-to-dark-matter ratio of Ω_dark/Ω_visible ≃ 5, these composite compact objects are expected to represent a potentially ubiquitous electromagnetic background radiation by interacting with ordinary baryonic matter. We conducted an in-depth analysis of axion quark nugget-baryonic matter interactions in the environment of the intracluster medium in the constrained cosmological Simulation of the LOcal Web (SLOW). Aims. Here, we aim to provide upper limit predictions on electromagnetic counterparts of axion quark nuggets in the environment of galaxy clusters by inferring their thermal and non-thermal emission spectrum originating from axion quark nugget-cluster gas interactions. Methods. We analyzed the emission of axion quark nuggets in a large sample of 161 simulated galaxy clusters using the SLOW simulation. These clusters are divided into a sub-sample of 150 galaxy clusters, ordered in five mass bins ranging from 0.8 to 31.7 × 10¹⁴ M_⊙, along with 11 cross-identified galaxy clusters from observations. We investigated dark matter-baryonic matter interactions in galaxy clusters in their present stage at the redshift of z = 0 by assuming all dark matter consists of axion quark nuggets. The resulting electromagnetic signatures were compared to thermal Bremsstrahlung and non-thermal cosmic ray (CR) synchrotron emission in each galaxy cluster. We further investigated individual frequency bands imitating the observable range of the WMAP, Planck, Euclid, and XRISM telescopes for the most promising cross-identified galaxy clusters hosting detectable signatures of axion quark nugget emission. Results. We observed a positive excess in the low- and high-energy frequency windows, where thermal and non-thermal axion quark nugget emission can significantly contribute to (or even outshine) the emission of the intracluster medium (ICM) in frequencies up to ν_T ≲ 3842.19 GHz and ν_T ϵ [3.97, 10.99] × 10¹⁰GHz, respectively. Emission signatures of axion quark nuggets are found to be observable if CR synchrotron emission of individual clusters is sufficiently low. The degeneracy in the parameters contributing to an emission excess makes it challenging to offer predictions with respect to pinpointing specific regions of a positive axion quark nugget excess; however, a general increase in the total galaxy cluster emission is expected based on this dark matter model. Axion quark nuggets constitute an increment of 4.80% of the total galaxy cluster emission in the low-energy regime of ν_T ≲ 3842.19 GHz for a selection of cross-identified galaxy clusters. We propose that the Fornax and Virgo clusters represent the most promising candidates in the search for axion quark nugget emission signatures. Conclusions. The results from our simulations point towards the possibility of detecting an axion quark nugget excess in galaxy clusters in observations if their signatures can be sufficiently disentangled from the ICM radiation. While this model proposes a promising explanation for the composition of dark matter, with the potential to have this outcome verified by observations, we propose further changes that are aimed at refining our methods. Our ultimate goal is to identify the extracted electromagnetic counterparts of axion quark nuggets with even greater precision in the near future.

Confinement prohibits isolation of color charges, e.g., quarks, in nature via a process called string breaking: the separation of two charges results in an increase in the energy of a color flux, visualized as a string, connecting those charges. Eventually, creating additional charges is energetically favored, hence breaking the string. Such a phenomenon can be probed in simpler models, including quantum spin chains, enabling enhanced understanding of string-breaking dynamics. A challenging task is to understand how string breaking occurs as time elapses, in an out-of-equilibrium setting. This work establishes the phenomenology of dynamical string breaking induced by a gradual increase of string tension over time. It, thus, goes beyond instantaneous quench processes and enables tracking the real-time evolution of strings in a more controlled setting. We focus on domain-wall confinement in a family of quantum Ising chains. Our results indicate that, for sufficiently short strings and slow evolution, string breaking can be described by the transition dynamics of a two-state quantum system akin to a Landau-Zener process. For longer strings, a more intricate spatiotemporal pattern emerges: the string breaks by forming a superposition of bubbles (domains of flipped spins of varying sizes), which involve highly excited states. We finally demonstrate that string breaking driven only by quantum fluctuations can be realized in the presence of sufficiently long-ranged interactions. This work holds immediate relevance for studying string breaking in quantum-simulation experiments.

We present a full general relativistic analytic solution for a radiation-pressure-supported equilibrium fluid torus orbiting a rotating neutron star (NS). We applied previously developed analytical methods that include the effects of both the NS's angular momentum and quadrupole moment in the Hartle-Thorne geometry. The structure, size, and shape of the torus are explored, with a particular focus on the critically thick solution – the cusp tori. For the astrophysically relevant range of NS parameters, we examined how our findings differ from those obtained for the Schwarzschild space-time. The solutions for rotating stars display signatures of an interplay between relativistic and Newtonian effects where the impact of the NS angular momentum and quadrupole moment are almost counterbalanced at a given radius. Nevertheless, the space-time parameters still strongly influence the size of tori, which can be shown in a coordinate-independent way. Finally, we discuss the importance of the size of the central NS which determines whether or not a surrounding torus exists. We provide a set of tools in a Wolfram Mathematica code, which establishes a basis for further investigation of the impact of the NSs' super-dense matter equation of state on the spectral and temporal behaviour of accretion tori.

Context. Dust coagulation and fragmentation impact the structure and evolution of protoplanetary disks and set the initial conditions for planet formation. Dust grains dominate the opacities, they determine the cooling times of the gas via thermal accommodation in collisions, they influence the ionization state of the gas, and the available grain surface area is an important parameter for the chemistry in protoplanetary disks. Therefore, dust evolution is an effect that should not be ignored in numerical studies of protoplanetary disks. Available dust coagulation models are, however, too computationally expensive to be implemented in large-scale hydrodynamic simulations. This limits detailed numerical studies of protoplanetary disks, including these effects, mostly to one-dimensional models. Aims. We aim to develop a simple – yet accurate – dust coagulation model that can be easily implemented in hydrodynamic simulations of protoplanetary disks. Our model shall not significantly increase the computational cost of simulations and provide information about the local grain size distribution. Methods. The local dust size distributions are assumed to be truncated power laws. Such distributions can be fully characterized by only two dust fluids (large and small grains) and a maximum particle size, truncating the power law. We compare our model to state- of-the-art dust coagulation simulations and calibrate it to achieve a good fit with these sophisticated numerical methods. Results. Running various parameter studies, we achieved a good fit between our simplified three-parameter model and DustPy, a state-of-the-art dust coagulation software. Conclusions. We present TriPoD, a sub-grid dust coagulation model for the PLUTO code. With TriPoD, we can perform twodimensional, vertically integrated dust coagulation simulations on top of a hydrodynamic simulation. Studying the dust distributions in two-dimensional vortices and planet-disk systems is thus made possible.

Two decades ago the $\chi_{c1}\left(3872\right)$ was discovered in the hadron spectrum with two heavy quarks. The discovery fueled a surge in experimental research, uncovering dozens of so called XYZ exotics states lying outside the conventional quark model, as well as theoretical investigations into new forms of matter, such as quark-gluon hybrids, mesonic molecules, and tetraquarks, with the potential of disclosing new information about the fundamental strong force. Among the XYZs, the $\chi_{c1}\left(3872\right)$ and $T_{cc}^+\left(3875\right)$ stand out for their striking characteristics and unlashed many discussions about their nature. Here, we address this question using the Born--Oppenheimer Effective Field Theory (BOEFT) and show how QCD settles the issue of their composition. Not only we describe well the main features of the $\chi_{c1}\left(3872\right)$ and $T_{cc}^+\left(3875\right)$ but obtain also model independent predictions in the bottomonium sector. This opens the way to systematic applications of BOEFT to all XYZs.

Recent observations of volume-limited samples of magnetic white dwarfs (WD) have revealed a higher incidence of magnetism in older stars. Specifically, these studies indicate that magnetism is more prevalent in WDs with fully or partially crystallized cores than in those with entirely liquid cores. This has led to the recognition of a crystallization-driven dynamo as an important mechanism for explaining magnetism in isolated WDs. However, recent simulations have challenged the capability of this mechanism to generate surface magnetic fields with the typical strengths detected in WDs. In this Letter, we explore an alternative hypothesis for the surface emergence of magnetic fields in isolated WDs. Those with masses ≳0.55 M_⊙ are the descendants of main sequence stars with convective cores capable of generating strong dynamo magnetic fields. This idea is supported by asteroseismic evidence of strong magnetic fields buried within the interiors of red giant branch stars. Assuming that these fields are disrupted by subsequent convective zones, we estimated magnetic breakout times for WDs with carbon-oxygen (CO) cores and masses ranging from 0.57 M_⊙ to 1.3 M_⊙. Due to the significant uncertainties in breakout times stemming from the treatment of convective boundaries and mass-loss rates, we cannot provide a precise prediction for the emergence time of the main sequence dynamo field. However, we can predict that this emergence should occur during the WD phase for those objects with masses ≳0.65 M_⊙. We also find that the magnetic breakout is expected to occur earlier in more massive WDs, which is consistent with observations of volume-limited samples and the well-established fact that magnetic WDs tend to be more massive than non-magnetic ones. Moreover, within the uncertainties of stellar evolutionary models, we find that the emergence of main sequence dynamo magnetic fields can account for a significant portion of the magnetic WDs. Additionally, we estimated magnetic breakout times due to crystallization-driven dynamos in CO WDs; our results suggest that this mechanism cannot explain the majority of magnetic WDs.

Context. The formation and evolution of protoplanetary disks remains elusive. We have numerous astronomical observations of young stellar objects of different ages with their envelopes and/or disks. Moreover, in the last decade, there has been tremendous progress in numerical simulations of star and disk formation. New simulations use realistic equations of state for the gas and treat the interaction of matter and the magnetic field with the full set of nonideal magnetohydrodynamic (MHD) equations. However, it is still not fully clear how a disk forms and whether it happens from inside-out or outside-in. Open questions remain regarding where material is accreted onto the disk and comes from, how dust evolves in disks, and the timescales of appearance of disk's structures. These unknowns limit our understanding of how planetesimals and planets form and evolve. Aims. We attempted to reconstruct the evolutionary history of the protosolar disk, guided by the large amount of cosmochemical constraints derived from the study of meteorites, while using astronomical observations and numerical simulations as a guide to pinpointing plausible scenarios. Methods. Our approach is highly interdisciplinary and we do not present new observations or simulations in this work. Instead, we combine, in an original manner, a large number of published results concerning young stellar objects observations, and numerical simulations, along with the chemical, isotopic and petrological nature of meteorites. Results. We have achieved a plausible and coherent view of the evolution of the protosolar disk that is consistent with cosmochemical constraints and compatible with observations of other protoplanetary disks and sophisticated numerical simulations. The evidence that high-temperature condensates, namely, calcium-aluminum inclusions (CAIs) and amoeboid olivine aggregates (AOAs), formed near the protosun before being transported to the outer disk can be explained in two ways: there could have either been an early phase of vigorous radial spreading of the disk that occurred or fast transport of these condensates from the vicinity of the protosun toward large disk radii via the protostellar outflow. The assumption that the material accreted toward the end of the infall phase was isotopically distinct allows us to explain the observed dichotomy in nucleosynthetic isotopic anomalies of meteorites. It leads us toward intriguing predictions on the possible isotopic composition of refractory elements in comets. At a later time, when the infall of material waned, the disk started to evolve as an accretion disk. Initially, dust drifted inward, shrinking the radius of the dust component to ∼45 au, probably about to about half of the width of the gas component. Next, structures must have emerged, producing a series of pressure maxima in the disk, which trapped the dust on Myr timescales. This allowed planetesimals to form at radically distinct times without significantly changing any of the isotopic properties. We also conclude that there was no late accretion of material onto the disk via streamers. The disk disappeared at about 5 My, as indicated by paleomagnetic data in meteorites. Conclusions. The evolution of the protosolar disk seems to have been quite typical in terms of size, lifetime, and dust behavior. This suggests that the peculiarities of the Solar System with respect to extrasolar planetary systems probably originate from the chaotic nature of planet formation and not from the properties of the parental disk itself.

In recent years the amount of publicly available astronomical data has increased exponentially, with a remarkable example being large-scale multiepoch photometric surveys. This wealth of data poses challenges to the classical methodologies commonly employed in the study of variable objects. As a response, deep learning techniques are increasingly being explored to effectively classify, analyze, and interpret these large datasets. In this paper we use two-dimensional histograms to represent Optical Gravitational Lensing Experiment phasefolded light curves as images. We use a Convolutional Neural Network (CNN) to classify variable objects within eight different categories (from now on labels): Classical Cepheid, RR Lyrae, Long Period Variable, Miras, Ellipsoidal Binary, Delta Scuti, Eclipsing Binary, and spurious class with Incorrect Periods (Rndm). We set up different training sets to train the same CNN architecture in order to characterize the impact of the training. The training sets were built from the same source of labels but different filters and balancing techniques were applied. Namely: Undersampling, Data Augmentation, and Batch Balancing (BB). The best performance was achieved with the BB approach and a training sample size of ~370 000 stars. Regarding computational performance, the image representation production rate is of ~76 images per core per second, and the time to predict is ~60 μs per star. The accuracy of the classification improves from ~92%, when based only on the CNN, to ~98% when the results of the CNN are combined with the period and amplitude features in a two step approach. This methodology achieves comparable results with previous studies but with two main advantages: the identification of miscalculated periods and the improvement in computational time cost.

Context. The Spektrum-Roentgen-Gamma (SRG)/extended Roentgen Survey with an Imaging Telescope Array (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, that is, at z > 5.6. By selection, these quasars populate the bright end of the active galactic nuclei (AGN) X-ray luminosity function, and their space density offers a powerful demographic diagnostic of the parent super-massive black hole (SMBH) population. Aims. 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. Methods. 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 (up to its four-pass cumulative version, eRASS:4). 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. Results. 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, therefore, X-ray ultra-luminous by selection. We also report the detection of these quasars at radio frequencies. The first one is a broad absorption line quasar, which shows significant, order-of-magnitude X-ray dimming over 3.5 years, corresponding to six months in the quasar rest frame. The second X-ray detected quasar is a jetted source with compact morphology. We show that a blazar configuration is likely for this source, making it one of the most distant blazars known to date. Conclusions. 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 we will further explore in a future systematic full-hemisphere survey.

Context. Galaxy groups lying between galaxies and galaxy clusters in the mass spectrum of dark matter halos play a crucial role in the evolution and formation of the large-scale structure. Their shallower potential wells compared to clusters of galaxies make them excellent sources to constrain non-gravitational processes such as feedback from the central active galactic nuclei (AGN). Aims. We investigate the impact of feedback, particularly from AGN, on the entropy and characteristic temperature measurements of galaxy groups detected in the SRG/eROSITA's first All-Sky Survey (eRASS1) to shed light on the characteristics of the feedback mechanisms and help guide future AGN feedback implementations in numerical simulations. Methods. We analyzed the deeper eROSITA observations of 1178 galaxy groups detected in the eRASS1. We divided the sample into 271 subsamples based on their physical and statistical properties and extracted average thermodynamic properties, including the electron number density, temperature, and entropy, at three characteristic radii from cores to outskirts along with the integrated temperature by jointly analyzing X-ray images and spectra following a Bayesian approach. Results. We present the tightest constraints with unprecedented statistical precision on the impact of AGN feedback through our average entropy and characteristic temperature measurements of the largest group sample used in X-ray studies, incorporating major systematics in our analysis. We find that entropy shows an increasing trend with temperature in the form of a power-law-like relation at the higher intra-group medium (IGrM) temperatures, while for the low-mass groups with cooler (T < 1.44 keV) IGrM temperatures, a slight flattening is observed on the average entropy. Overall, the observed entropy measurements agree well with the earlier measurements in the literature. Additionally, comparisons with the state-of-the-art cosmological hydrodynamic simulations (MillenniumTNG, Magneticum, OWL) after applying the selection function calibrated for our galaxy groups reveal that observed entropy profiles in the cores are below the predictions of simulations. At the mid-region, the entropy measurements agree well with the Magneticum simulations, whereas the predictions of MillenniumTNG and OWL simulations fall below observations. At the outskirts, the overall agreement between the observations and simulations improves, with Magneticum simulations reproducing the observations the best. Conclusions. These measurements will pave the way for achieving more realistic AGN feedback implementations in numerical simulations. The future eROSITA Surveys will enable the extension of the entropy measurements in even cooler IGrM temperatures below 0. 5 keV, allowing for the testing of the AGN feedback models in this regime.

JWST observations of the young Galactic supernova remnant Cassiopeia A revealed an unexpected structure seen as a green emission feature in colored composite MIRI F1130W and F1280W images—hence dubbed the Green Monster—that stretches across the central parts of the remnant in projection. Combining the kinematic information from NIRSpec and the MIRI Medium Resolution Spectrograph with the multiwavelength imaging from NIRCam and MIRI, we associate the Green Monster with circumstellar material (CSM) that was lost during an asymmetric mass-loss phase. MIRI images are dominated by dust emission, but their spectra show emission lines from Ne, H, and Fe with low radial velocities indicative of a CSM nature. An X-ray analysis of this feature in a companion paper supports its CSM nature and detects significant blueshifting, thereby placing the Green Monster on the nearside, in front of the Cas A supernova remnant. The most striking features of the Green Monster are dozens of almost perfectly circular 1″–3″ sized holes, most likely created by interaction between high-velocity supernova ejecta material and the CSM. Further investigation is needed to understand whether these holes were formed by small 8000–10,500 km s^‑1 N-rich ejecta knots that penetrated and advanced out ahead of the remnant's 5000–6000 km s^‑1 outer blast wave or by narrow ejecta fingers that protrude into the forward-shocked CSM. The detection of the Green Monster provides further evidence of the highly asymmetric mass loss that Cas A's progenitor star underwent prior to its explosion.

We evaluate the three-loop five-point pentagon-box-box massless integral family in the dimensional regularization scheme, via canonical differential equation. We use tools from computational algebraic geometry to enable the necessary integral reductions. The boundary values of the differential equation are determined analytically in the Euclidean region. To express the final result, we introduce a new representation of weight six functions in terms of one-fold integrals over the product of weight-three functions with weight-two kernels that are derived from the differential equation. Our work paves the way to the analytic computation of three-loop multi-leg Feynman integrals.

The ongoing discrepancy in the Hubble constant ($H_0$) estimates obtained through local distance ladder methods and early universe observations poses a significant challenge to the $\Lambda$CDM model, suggesting potential new physics. Type II supernovae (SNe II) offer a promising technique for determining $H_0$ in the local universe independently of the traditional distance ladder approach, opening up a complimentary path for testing this discrepancy. We aim to provide the first $H_0$ estimate using the tailored expanding photosphere method (EPM) applied to SNe II, made possible by recent advancements in spectral modelling that enhance its precision and efficiency. Our tailored EPM measurement utilizes a spectral emulator to interpolate between radiative transfer models calculated with TARDIS, allowing us to fit supernova spectra efficiently and derive self-consistent values for luminosity-related parameters. We apply the method on public data for ten SNe II at redshifts between 0.01 and 0.04. Our analysis demonstrates that the tailored EPM allows for $H_0$ measurements with precision comparable to the most competitive established techniques, even when applied to literature data not designed for cosmological applications. We find an independent $H_0$ value of $74.9\pm1.9$ (stat) km/s/Mpc, which is consistent with most current local measurements. Considering dominant sources of systematic effects, we conclude that our systematic uncertainty is comparable to or less than the current statistical uncertainty. This proof-of-principle study highlights the potential of the tailored EPM as a robust and precise tool for investigating the Hubble tension independently of the local distance ladder. Observations of SNe II tailored to $H_0$ estimation can make this an even more powerful tool by improving the precision and by allowing us to better understand and control systematic uncertainties.

The impact of viscosity in the intracluster medium (ICM) is still an open question in astrophysics. To address this problem, we have run a set of cosmological simulations of three galaxy clusters with a mass larger than M _Vir > 10¹⁵ M _⊙ at z = 0 using the smoothed particle magnetohydrodynamics-code OPENGADGET3. We aim to quantify the influence of viscosity and constrain its value in the ICM. Our results show significant morphological differences at small scales, temperature variations, and density fluctuations induced by viscosity. We observe a suppression of instabilities at small scales, resulting in a more filamentary structure and a larger amount of small structures due to the lack of mixing with the medium. The conversion of kinetic to internal energy leads to an increase of the virial temperature of the cluster of ∼5%–10%, while the denser regions remain cold. The amplitude of density and velocity fluctuations are found to increase with viscosity. However, comparison with observational data indicates that the simulations, regardless of the viscosity, match the observed slope of the amplitude of density fluctuations, challenging the direct constraint of viscosity solely through density fluctuations. Furthermore, the ratio of density to velocity fluctuations remains close to 1 regardless of the amount of viscosity, in agreement with the theoretical expectations. Our results show for the first time in a cosmological simulation of a galaxy cluster the effect of viscosity in the ICM, a study that is currently missing in the literature.

We compute the photon self-energy to three loops in Quantum Electrodynamics. The method of differential equations for Feynman integrals and a complete $\epsilon$-factorization of the former allow us to obtain fully analytical results in terms of iterated integrals involving integration kernels related to a K3 geometry. We argue that our basis has the right properties to be a natural generalization of a canonical basis beyond the polylogarithmic case and we show that many of the kernels appearing in the differential equations, cancel out in the final result to finite order in $\epsilon$. We further provide generalized series expansions that cover the whole kinematic space so that our results for the self-energy may be easily evaluated numerically for all values of the momentum squared. From the local solution at $p^2=0$, we extract the photon wave function renormalization constant in the on-shell scheme to three loops and confirm its agreement with previously obtained results.

Forward modeling the galaxy density within the Effective Field Theory of Large Scale Structure (EFT of LSS) enables field-level analyses that are robust to theoretical uncertainties. At the same time, they can maximize the constraining power from galaxy clustering on the scales amenable to perturbation theory. In order to apply the method to galaxy surveys, the forward model must account for the full observational complexity of the data. In this context, a major challenge is the inclusion of redshift space distortions (RSDs) from the peculiar motion of galaxies. Here, we present improvements in the efficiency and accuracy of the RSD modeling in the perturbative LEFTfield forward model. We perform a detailed quantification of the perturbative and numerical error for the prediction of momentum, velocity and the redshift-space matter density. Further, we test the recovery of cosmological parameters at the field level, namely the growth rate $f$, from simulated halos in redshift space. For a rigorous test and to scan through a wide range of analysis choices, we fix the linear (initial) density field to the known ground truth but marginalize over all unknown bias coefficients and noise amplitudes. With a third-order model for gravity and bias, our results yield $<1\,\%$ statistical and $<1.5\,\%$ systematic error. The computational cost of the redshift-space forward model is only $\sim 1.5$ times of the rest frame equivalent, enabling future field-level inference that simultaneously targets cosmological parameters and the initial matter distribution.

We use galaxy cluster abundance measurements from the South Pole Telescope (SPT) enhanced by Multi-Component Matched Filter (MCMF) confirmation and complemented with mass information obtained using weak-lensing data from Dark Energy Survey Year~3 (DES Y3) and targeted Hubble Space Telescope (HST) observations for probing deviations from the cold dark matter paradigm. Concretely, we consider a class of dark sector models featuring interactions between dark matter (DM) and a dark radiation (DR) component within the framework of the Effective Theory of Structure Formation (ETHOS). We focus on scenarios that lead to power suppression over a wide range of scales, and thus can be tested with data sensitive to large scales, as realized for example for DM$-$DR interactions following from an unbroken non-Abelian $SU(N)$ gauge theory (interaction rate with power-law index $n=0$ within the ETHOS parameterization). Cluster abundance measurements are mostly sensitive to the amount of DR interacting with DM, parameterized by the ratio of DR temperature to the cosmic microwave background (CMB) temperature, $\xi_{\rm DR}=T_{\rm DR}/T_{\rm CMB}$. We find an upper limit $\xi_{\rm DR}<17\%$ at $95\%$ credibility. When the cluster data are combined with Planck 2018 CMB data along with baryon acoustic oscillation (BAO) measurements we find $\xi_{\rm DR}<10\%$, corresponding to a limit on the abundance of interacting DR that is around three times tighter than that from CMB+BAO data alone. We also discuss the complementarity of weak lensing informed cluster abundance studies with probes sensitive to smaller scales, explore the impact on our analysis of massive neutrinos, and comment on a slight preference for the presence of a non-zero interacting DR abundance, which enables a physical solution to the $S_8$ tension.

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

The stellar disc is the dominant luminous component of the Milky Way (MW). Although our understanding of its structure is rapidly expanding due to advances in large-scale stellar surveys, our picture of the MW disc remains substantially obscured by selection functions and incomplete spatial coverage of observational data. In this work, we present the comprehensive chrono-chemo-kinematic structure of the MW disc, recovered using a novel orbit superposition approach combined with data from APOGEE DR 17. We detect periodic azimuthal metallicity variations within 6-8 kpc with an amplitude of 0.05-0.1 dex peaking along the bar major axis. The radial metallicity profile of the MW also varies with azimuth, displaying a pattern typical among other disc galaxies: a decline outside the solar radius and an almost flat profile in the inner region, attributed to the presence of old, metal-poor high-{\alpha} populations, which comprise about 40% of the total stellar mass. The geometrically defined thick disc and the high-{\alpha} populations have comparable masses, with differences in their stellar population content, which we quantify using the reconstructed 3D MW structure. The well-known [{\alpha}/Fe]-bimodality in the MW disc, once weighted by stellar mass, is less pronounced at a given metallicity for the whole galaxy but distinctly visible in a narrow range of galactic radii (5-9 kpc), explaining its relative lack of prominence in external galaxies and galaxy formation simulations. Analysing a more evident double age-abundance sequence, we construct a scenario for the MW disc formation, advocating for an inner/outer disc dichotomy genetically linked to the MW's evolutionary stages. In this picture, the extended solar vicinity is a transition zone that shares chemical properties of both the inner (old age-metallicity sequence) and outer discs (young age-metallicity sequence).

A dense neutrino gas exhibiting angular crossings in the electron lepton number is unstable and develops fast flavor conversions. Instead of assuming an unstable configuration from the onset, we imagine that the system is externally driven toward instability. We use the simplest model of two neutrino beams initially of different flavor that either suddenly appear or one or both slowly build up. Flavor conversions commence well before the putative unstable state is fully attained, and the final outcome depends on how the system is driven. The system generally sticks to the closest state that is linearly stable, a conclusion that we prove for the first time using quasilinear theory. Our results suggest that in an astrophysical setting, one should focus less on flavor instabilities in the neutrino radiation field and more on the external dynamics that leads to the formation of the unstable state.

We report the discovery of the first example of an Einstein zig-zag lens, an extremely rare lensing configuration. In this system, J1721+8842, six images of the same background quasar are formed by two intervening galaxies, one at redshift z₁ = 0.184 and a second one at z₂ = 1.885. Two out of the six multiple images are deflected in opposite directions as they pass the first lens galaxy on one side, and the second on the other side -- the optical paths forming zig-zags between the two deflectors. In this letter, we demonstrate that J1721+8842, previously thought to be a lensed dual quasar, is in fact a compound lens with the more distant lens galaxy also being distorted as an arc by the foreground galaxy. Evidence supporting this unusual lensing scenario includes: 1- identical light curves in all six lensed quasar images obtained from two years of monitoring at the Nordic Optical Telescope; 2- detection of the additional deflector at redshift z₂ = 1.885 in JWST/NIRSpec IFU data; and 3- a multiple-plane lens model reproducing the observed image positions. This unique configuration offers the opportunity to combine two major lensing cosmological probes: time-delay cosmography and dual source-plane lensing since J1721+8842 features multiple lensed sources forming two distinct Einstein radii of different sizes, one of which being a variable quasar. We expect tight constraints on the Hubble constant and the equation of state of dark energy by combining these two probes on the same system. The z₂=1.885 deflector, a quiescent galaxy, is also the highest-redshift strong galaxy-scale lens with a spectroscopic redshift measurement.

Using the decomposition of the D-dimensional space-time into parallel and perpendicular subspaces, we study and prove a connection between Landau and leading singularities for N-point one-loop Feynman integrals by applying the multidimensional theory of residues. We show that if <inline-formula id="IEq1"><mml:math><mml:mrow><mml:mi>D</mml:mi><mml:mo>=</mml:mo><mml:mi>N</mml:mi></mml:mrow></mml:math></inline-formula> and <inline-formula id="IEq2"><mml:math><mml:mrow><mml:mi>D</mml:mi><mml:mo>=</mml:mo><mml:mi>N</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:math></inline-formula>, the leading singularity corresponds to the inverse of the square root of the leading Landau singularity of the first and second type, respectively. We make use of this outcome to systematically provide differential equations of Feynman integrals in canonical forms and the extension of the connection of these singularities at the multi-loop level by exploiting the loop-by-loop approach. Illustrative examples with the calculation of Landau and leading singularities are provided to supplement our results.

We compute the electron self-energy in Quantum Electrodynamics to three loops in terms of iterated integrals over kernels of elliptic type. We make use of the differential equations method, augmented by an ϵ-factorized basis, which allows us to gain full control over the differential forms appearing in the iterated integrals to all orders in the dimensional regulator. We obtain compact analytic expressions, for which we provide generalized series expansion representations that allow us to evaluate the result numerically for all values of the electron momentum squared. As a by-product, we also obtain ϵ-resummed results for the self-energy in the on-shell limit p² = m², which we use to recompute the known three-loop renormalization constants in the on-shell scheme.

We investigate the formation of bound states of non-relativistic dark matter particles subject to long-range interactions through radiative capture. The initial scattering and final bound states are described by Coulomb potentials with different strengths, as relevant for non-abelian gauge interactions or theories featuring charged scalars. For bound states with generic quantum numbers $n$ and $\ell$, we provide closed-form expressions for the bound-state formation (BSF) cross sections of monopole, dipole and quadrupole transitions, and of arbitrary multipole order when $\ell=n-1$. This allows us to investigate in detail a strong enhancement of BSF that occurs for initial states in a repulsive potential. For $\ell=n-1\gg 1$, we show that the BSF cross section for each single bound state violates the perturbative unitarity bound in the vicinity of a certain critical initial velocity, and provide an interpretation in terms of a smooth matching of classical trajectories. When summing the BSF cross section over all possible bound states in the final state, this leads to a unitarity violation below a certain velocity, but within the validity range of the weakly coupled non-relativistic description. We identify an effectively strong interaction as the origin of this unitarity violation, which is caused by an "anomalously" large overlap of scattering and bound-state wave functions in Coulomb potentials of different strength.

We compute the tree-level current for the emission of two soft quark-antiquark pairs in a hard scattering. We also compute the square of this current and discuss the resulting color correlations, featuring dipole correlations and three-parton correlations. This object is essential for analyzing the infrared singularities at next-to-next-to-next-to-next-to-leading-order (N4LO) predictions in perturbative QCD.

In this paper, the first in a series, we present a new theoretical model for the global structure and dissipation of relativistically magnetized collisionless shock waves. Quite remarkably, we find that in contrast to unmagnetized shocks, energy dissipation does not involve collective plasma interactions. Rather, it is a consequence of nonlinear particle dynamics. We demonstrate that the kinetic-scale shock transition can be modeled as a stationary system consisting of a large set of cold beams coupled through the magnetic field. The fundamental mechanism governing shock dissipation relies on the onset of chaos in orbital dynamics within quasiperiodic solitonic structures. We discuss the impact of upstream temperature and magnetization on the shock profile, recovering the magnetic field compression, downstream velocities, and heating expected from the Rankine-Hugoniot jump conditions. We deduce a rate of entropy generation from the spectrum of Lyapunov exponents and discuss the thermalization of the beam distribution. Our model provides a general framework to study magnetized collisionless shock structures.