Multiply lensed sources experience a relative time delay in the arrival of photons. This effect can be used to measure absolute distances and the Hubble constant ($H_0$) and is known as time-delay cosmography. The methodology is independent of the local distance ladder and early-universe physics and provides a precise and competitive measurement of $H_0$. With upcoming observatories, time-delay cosmography can provide a 1% precision measurement of $H_0$ and can decisively shed light on the current reported 'Hubble tension'. This paper presents the theoretical background and the current techniques applied for time-delay cosmographic studies and the measurement of the Hubble constant. The paper describes the challenges and systematics in the different components of the analysis and strategies to mitigate them. The current measurements are discussed in context and the opportunities with the anticipated data sets in the future are laid out.

We study the effect of super-sample covariance (SSC) on the power spectrum and higher-order statistics: bispectrum, halo mass function, and void size function. We also investigate the effect of SSC on the cross-covariance between the statistics. We consider both the matter and halo fields. Higher-order statistics of the large-scale structure contain additional cosmological information beyond the power spectrum and are a powerful tool to constrain cosmology. They are a promising probe for ongoing and upcoming high precision cosmological surveys such as DESI, PFS, Rubin Observatory LSST, Euclid, SPHEREx, SKA, and Roman Space Telescope. Cosmological simulations used in modeling and validating these statistics often have sizes that are much smaller than the observed Universe. Density fluctuations on scales larger than the simulation box, known as super-sample modes, are not captured by the simulations and in turn can lead to inaccuracies in the covariance matrix. We compare the covariance measured using simulation boxes containing super-sample modes to those without. We also compare with the Separate Universe approach. We find that while the power spectrum, bispectrum and halo mass function show significant scale- or mass-dependent SSC, the void size function shows relatively small SSC. We also find significant SSC contributions to the cross-covariances between the different statistics, implying that future joint-analyses will need to carefully take into consideration the effect of SSC.

The early release science results from the $\textit{James Webb Space Telescope (JWST)}$ have yielded an unexpected abundance of high-redshift luminous galaxies that seems to be in tension with current theories of galaxy formation. However, it is currently difficult to draw definitive conclusions form these results as the sources have not yet been spectroscopically confirmed. It is in any case important to establish baseline predictions from current state-of-the-art galaxy formation models that can be compared and contrasted with these new measurements. In this work, we use the new large-volume ($L_\mathrm{box}\sim 740 \, \mathrm{cMpc}$) hydrodynamic simulation of the MillenniumTNG project to make predictions for the high-redshift ($z\gtrsim8$) galaxy population and compare them to recent $\textit{JWST}$ observations. We show that the simulated galaxy population is broadly consistent with observations until $z\sim10$. From $z\approx10-12$, the observations indicate a preference for a galaxy population that is largely dust-free, but is still consistent with the simulations. Beyond $z\gtrsim12$, however, our simulation results underpredict the abundance of luminous galaxies and their star-formation rates by almost an order of magnitude. This indicates either an incomplete understanding of the new $\textit{JWST}$ data or a need for more sophisticated galaxy formation models that account for additional physical processes such as Population~III stars, variable stellar initial mass functions, or even deviations from the standard $\Lambda$CDM model. We emphasise that any new process invoked to explain this tension should only significantly influence the galaxy population beyond $z\gtrsim10$, while leaving the successful galaxy formation predictions of the fiducial model intact below this redshift.

Winds in protoplanetary disks play an important role in their evolution and dispersal. However, what physical process is driving the winds is still unclear (i.e. magnetically vs thermally driven), and can only be understood by directly confronting theoretical models with observational data. We use hydrodynamic photoevaporative disk-wind models and post-process them with a thermo-chemical model to produce synthetic observables for the o-H$_2$ at 2.12 micron and [OI] at 0.63 micron spectral lines and directly compare the results to a sample of observations. Our photoevaporative disk-wind model is consistent with the observed signatures of the blueshifted narrow low-velocity component (NLVC), which is usually associated with slow disk winds, for both tracers. Only for one out of seven targets that show blueshifted NLVCs does the photoevaporative model fail to explain the observed line kinematics. Our results also indicate that interpreting spectral line profiles by simple methods, such as the thin-disk approximation, to determine the line emitting region can yield misleading conclusions. The photoevaporative disk-wind models are largely consistent with the studied observational data set, but it is not possible to clearly discriminate between different wind-driving mechanisms. Further improvements to the models, such as consistent modelling of the dynamics and chemistry and detailed modelling of individual targets would be beneficial. Furthermore, a direct comparison of magnetically driven disk-wind models to the observational data set is necessary in order to determine whether or not spatially unresolved observations of multiple wind tracers are sufficient to discriminate between theoretical models.

Under some assumptions on the hierarchy of relevant energy scales, we compute the nonrelativistic QCD (NRQCD) long-distance matrix elements (LDMEs) for inclusive production of $J/\psi$, $\psi(2S)$, and $\Upsilon$ states based on the potential NRQCD (pNRQCD) effective field theory. Based on the pNRQCD formalism, we obtain expressions for the LDMEs in terms of the quarkonium wavefunctions at the origin and universal gluonic correlators, which do not depend on the heavy quark flavor or the radial excitation. This greatly reduces the number of nonperturbative unknowns and substantially enhances the predictive power of the nonrelativistic effective field theory formalism. We obtain improved determinations of the LDMEs for $J/\psi$, $\psi(2S)$, and $\Upsilon$ states thanks to the universality of the gluonic correlators, and obtain phenomenological results for cross sections and polarizations at large transverse momentum that agree well with measurements at the LHC.

Dark matter (DM) with self-interactions is a promising solution for the small-scale problems of the standard cosmological model. Here we perform the first cosmological simulation of frequent DM self-interactions, corresponding to small-angle DM scatterings. The focus of our analysis lies in finding and understanding differences to the traditionally assumed rare DM (large-angle) self-scatterings. For this purpose, we compute the distribution of DM densities, the matter power spectrum, the two-point correlation function, and the halo and subhalo mass functions. Furthermore, we investigate the density profiles of the DM haloes and their shapes. We find that overall large-angle and small-angle scatterings behave fairly similarly with a few exceptions. In particular, the number of satellites is considerably suppressed for frequent compared to rare self-interactions with the same cross-section. Overall, we observe that while differences between the two cases may be difficult to establish using a single measure, the degeneracy may be broken through a combination of multiple ones. For instance, the combination of satellite counts with halo density or shape profiles could allow discriminating between rare and frequent self-interactions. As a by-product of our analysis, we provide - for the first time - upper limits on the cross-section for frequent self-interactions.

Nuclear weak decays provide important probes to fundamental symmetries in nature. A precise description of these processes in atomic nuclei requires comprehensive knowledge on both the strong and weak interactions in the nuclear medium and on the dynamics of quantum many-body systems. In particular, an observation of the hypothetical double beta decay without emission of neutrinos (0 νββ) would unambiguously demonstrate the Majorana nature of neutrinos and the existence of the lepton-number-violation process. It would also provide unique information on the ordering and absolute scale of neutrino masses. The next-generation tonne-scale experiments with sensitivity up to 10²⁸ years after a few years of running will probably provide a definite answer to these fundamental questions based on our current knowledge on the nuclear matrix element (NME), the precise determination of which is a challenge to nuclear theory. Beyond-mean-field approaches have been frequently adapted for the study of nuclear structure and decay throughout the nuclear chart for several decades. In this review, we summarize the status of beyond-mean-field calculations of the NMEs of 0 νββ decay assuming the standard mechanism of an exchange of light Majorana neutrinos. The challenges and prospects in the extension and application of beyond-mean-field approaches for 0 νββ decay are discussed.

Majoron-like bosons would emerge from a supernova (SN) core by neutrino coalescence of the form $\nu\nu\to\phi$ and $\bar\nu\bar\nu\to\phi$ with 100 MeV-range energies. Subsequent decays to (anti)neutrinos of all flavors provide a flux component with energies much larger than the usual flux from the "neutrino sphere." The absence of 100 MeV-range events in the Kamiokande II and IMB signal of SN 1987A implies that $\lesssim0.03$ of the total energy was thus emitted and provides the strongest constraint on the majoron-neutrino coupling of $g\lesssim 10^{-9}\,{\rm MeV}/m_\phi$ for $100~{\rm eV}\lesssim m_\phi\lesssim100~{\rm MeV}$. It is straightforward to extend our new argument to other hypothetical feebly interacting particles.

The Sunyaev-Zeldovich (SZ) effect is a powerful tool in modern cosmology. With future observations promising ever improving SZ measurements, the relativistic corrections to the SZ signals from galaxy groups and clusters are increasingly relevant. As such, it is important to understand the differences between three temperature measures: (a) the average relativistic SZ (rSZ) temperature, (b) the mass-weighted temperature relevant for the thermal SZ (tSZ) effect, and (c) the X-ray spectroscopic temperature. In this work, we compare these cluster temperatures, as predicted by the BAHAMAS & MACSIS, ILLUSTRISTNG, MAGNETICUM, and THE THREE HUNDRED PROJECT simulations. Despite the wide range of simulation parameters, we find the SZ temperatures are consistent across the simulations. We estimate a $\simeq 10{{\ \rm per\ cent}}$ level correction from rSZ to clusters with Y ≃ 10^-4 Mpc^-2. Our analysis confirms a systematic offset between the three temperature measures; with the rSZ temperature $\simeq 20{{\ \rm per\ cent}}$ larger than the other measures, and diverging further at higher redshifts. We demonstrate that these measures depart from simple self-similar evolution and explore how they vary with the defined radius of haloes. We investigate how different feedback prescriptions and resolutions affect the observed temperatures, and discover the SZ temperatures are rather insensitive to these details. The agreement between simulations indicates an exciting avenue for observational and theoretical exploration, determining the extent of relativistic SZ corrections. We provide multiple simulation-based fits to the scaling relations for use in future SZ modelling.

In compact astrophysical objects, the neutrino density can be so high that neutrino-neutrino refraction can lead to fast flavor conversion of the kind $\nu_e \bar\nu_e \leftrightarrow \nu_x \bar\nu_x$ with $x=\mu,\tau$, depending on the neutrino angle distribution. Previously, we have shown that in a homogeneous, axisymmetric two-flavor system, these collective solutions evolve in analogy to a gyroscopic pendulum. In flavor space, its deviation from the weak-interaction direction is quantified by a variable $\cos\vartheta$ that moves between $+1$ and $\cos\vartheta_{\rm min}$, the latter following from a linear mode analysis. As a next step, we include collisional damping of flavor coherence, assuming a common damping rate $\Gamma$ for all modes. Empirically we find that the damped pendular motion reaches an asymptotic level of pair conversion $f=A+(1-A)\cos\vartheta_{\rm min}$ (numerically $A\simeq 0.370$) that does not depend on details of the angular distribution (except for fixing $\cos\vartheta_{\rm min}$), the initial seed, nor $\Gamma$. On the other hand, even a small asymmetry between the neutrino and antineutrino damping rates strongly changes this picture and can even enable flavor instabilities in otherwise stable systems. Furthermore, we establish a formal connection with a stationary and inhomogeneous neutrino ensemble, showing that our findings also apply to this system.

We study the allowed parameter space of the scalar sector in the superweak extension of the standard model (SM). The allowed region is defined by the following conditions: (i) stability of the vacuum, (ii) perturbativity up to the Planck scale, and (iii) the pole mass of the Higgs boson falling into its experimentally measured range. We employ renormalization group equations and quantum corrections at two-loop accuracy. We study the dependence on the Yukawa couplings of the sterile neutrinos at selected values. We also check the exclusion limit set by the precise measurement of the mass of the W boson. Our method for constraining the parameter space using two-loop predictions can also be applied to simpler models such as the singlet scalar extension of the SM in a straightforward way.

We explore the effective field theory of a vector field X^μ that has a Stückelberg mass. The absence of a gauge symmetry for X^μ implies Lorentz-invariant operators are constructed directly from X^μ. Beyond the kinetic and mass terms, allowed interactions at the renormalizable level include X_μX^μH^†H , (X_μX^μ)², and X_μj^μ, where j^μ is a global current of the SM or of a hidden sector. We show that all of these interactions lead to scattering amplitudes that grow with powers of √{s }/m_X, except for the case of X_μj^μ where j^μ is a nonanomalous global current. The latter is well known when X is identified as a dark photon coupled to the electromagnetic current, often written equivalently as kinetic mixing between X and the photon. The power counting for the energy growth of the scattering amplitudes is facilitated by isolating the longitudinal enhancement. We examine in detail the interaction with an anomalous global vector current X_μj_anom^μ, carefully isolating the finite contribution to the fermion triangle diagram. We calculate the longitudinally-enhanced observables Z →X γ (when m_X<m_Z), f f ¯→X γ , and Z γ →Z γ when X couples to the baryon number current. Introducing a "fake" gauge-invariance by writing X^μ=A^μ-∂^μπ /m_X, the would-be gauge anomaly associated with A_μj_anom^μ is canceled by j_anom^μ∂_μπ /m_X; this is the four-dimensional Green-Schwarz anomaly-cancellation mechanism at work. Our analysis demonstrates there is a much larger set of possible interactions that an EFT with a Stückelberg vector field can have, revealing scattering amplitudes that grow with energy. The growth of these amplitudes can be tamed by a dark Higgs sector, but this requires dark Higgs boson interactions (and reintroduces fine-tuning in the dark Higgs sector) that can be separated from X interactions only in the limit g ≪1 .

We present the first quantitative spectral analysis of blue supergiant stars in the nearby galaxy NGC 2403. Out of a sample of 47 targets observed with the LRIS spectrograph at the Keck I telescope we have extracted 16 B- and A-type supergiants for which we have data of sufficient quality to carry out a comparison with model spectra of evolved massive stars and infer the stellar parameters. The radial metallicity gradient of NGC 2403 that we derive has a slope of -0.14 (+/- 0.05) dex/r_e, and is in accordance with the analysis of H II region oxygen abundances. We present evidence that the stellar metallicities that we obtain in extragalactic systems in general agree with the nebular abundances based on the analysis of the auroral lines, over more than one order of magnitude in metallicity. Adopting the known relation between stellar parameters and intrinsic luminosity we find a distance modulus m-M = 27.38 +/- 0.08 mag. While this can be brought into agreement with Cepheid-based determinations, it is 0.14 mag short of the value measured from the tip of the red giant branch. We update the mass-metallicity relation secured from chemical abundance studies of stars in resolved star-forming galaxies.

The impact of axion-like particles on the light polarization around the horizon of suppermassive black hole (SMBH) is discussed in the light of the latest polarization measurement of the Event Horizon Telescope (EHT). We investigate different sources of the polarization due to axion interaction with photons and the magnetic field of SMBH. These can modify the linear and circular polarization parameters of the emitted light. We have shown that a significant circular polarization can be produced via the photon scattering from the background magnetic field with axions as off-shell particles. This can further constrain the parameter space of ultralight axion-like particles and their couplings with photons. The future precise measurements of circular polarization can probe the features of ultralight axions in the near vicinity of SMBH.

While primordial black holes (PBHs) with masses M_PBH≳10^-11 M_⊙ cannot comprise the entirety of dark matter, the existence of even a small population of these objects can have profound astrophysical consequences. A subdominant population of PBHs will efficiently accrete dark matter particles before matter-radiation equality, giving rise to high-density dark matter spikes. We consider here the scenario in which dark matter is comprised primarily of weakly interacting massive particles (WIMPs) with a small subdominant contribution coming from PBHs, and revisit the constraints on the annihilation of WIMPs in these spikes using observations of the isotropic gamma-ray background (IGRB) and the cosmic microwave background (CMB), for a range of WIMP masses, annihilation channels, cross sections, and PBH mass functions. We find that the constraints derived using the IGRB have been significantly overestimated (in some cases by many orders of magnitude), and that limits obtained using observations of the CMB are typically stronger than, or comparable to, those coming from the IGRB. Importantly, we show that ∼O (M_⊙) PBHs can still contribute significantly to the dark matter density for sufficiently low WIMP masses and p-wave annihilation cross sections.

Context. The Gl 486 system consists of a very nearby, relatively bright, weakly active M3.5 V star at just 8 pc with a warm transiting rocky planet of about 1.3 R_⊕ and 3.0 M_⊕. It is ideal for both transmission and emission spectroscopy and for testing interior models of telluric planets.
Aims: To prepare for future studies, we aim to thoroughly characterise the planetary system with new accurate and precise data collected with state-of-the-art photometers from space and spectrometers and interferometers from the ground.
Methods: We collected light curves of seven new transits observed with the CHEOPS space mission and new radial velocities obtained with MAROON-X at the 8.1 m Gemini North telescope and CARMENES at the 3.5 m Calar Alto telescope, together with previously published spectroscopic and photometric data from the two spectrographs and TESS. We also performed near-infrared interferometric observations with the CHARA Array and new photometric monitoring with a suite of smaller telescopes (AstroLAB, LCOGT, OSN, TJO). This extraordinary and rich data set was the input for our comprehensive analysis.
Results: From interferometry, we measure a limb-darkened disc angular size of the star Gl 486 at θ_LDD = 0.390 ± 0.018 mas. Together with a corrected Gaia EDR3 parallax, we obtain a stellar radius R_* = 0.339 ± 0.015 R_⊕. We also measure a stellar rotation period at P_rot = 49.9 ± 5.5 days, an upper limit to its XUV (5-920 A) flux informed by new Hubble/STIS data, and, for the first time, a variety of element abundances (Fe, Mg, Si, V, Sr, Zr, Rb) and C/O ratio. Moreover, we imposed restrictive constraints on the presence of additional components, either stellar or sub-stellar, in the system. With the input stellar parameters and the radial-velocity and transit data, we determine the radius and mass of the planet Gl 486 b at R_p = 1.343_−0.062^+0.063 R_⊕ and M_p = 3.00_−0.12^+0.13 M_⊕, with relative uncertainties of the planet radius and mass of 4.7% and 4.2%, respectively. From the planet parameters and the stellar element abundances, we infer the most probable models of planet internal structure and composition, which are consistent with a relatively small metallic core with respect to the Earth, a deep silicate mantle, and a thin volatile upper layer. With all these ingredients, we outline prospects for Gl 486 b atmospheric studies, especially with forthcoming James Webb Space Telescope (Webb) observations.

The observed pattern of linear polarization of the cosmic microwave background photons is a sensitive probe of physics violating parity symmetry under inversion of spatial coordinates. A new parity-violating interaction might have rotated the plane of linear polarization by an angle β as the cosmic microwave background photons have been traveling for more than 13 billion years. This effect is known as "cosmic birefringence." In this paper, we present new measurements of cosmic birefringence from a joint analysis of polarization data from two space missions, P l a n c k and WMAP. This dataset covers a wide range of frequencies from 23 to 353 GHz. We measure β =0.342 °_-0.091^{° +0.094 °} [68% confidence level (CL)] for nearly full-sky data, which excludes β =0 at 99.987% CL. This corresponds to the statistical significance of 3.6 σ . There is no evidence for frequency dependence of β . We find a similar result, albeit with a larger uncertainty, when removing the Galactic plane from the analysis.

The flux of high energy neutrinos and photons produced in a blazar could get attenuated when they propagate through the dark matter spike around the central black hole and the halo of the host galaxy. Using the observation by IceCube of a few high-energy neutrino events from TXS 0506+056, and their coincident gamma ray events, we obtain new constraints on the dark matter-neutrino and dark matter-photon scattering cross sections. Our constraints are orders of magnitude more stringent than those derived from considering the attenuation through the intergalactic medium and the Milky Way dark matter halo. When the cross-section increases with energy, our constraints are also stronger than those derived from the CMB and large-scale structure.

Shape measurements of galaxies and galaxy clusters are widespread in the analysis of cosmological simulations. But the limitations of those measurements have been poorly investigated. In this paper, we explain why the quality of the shape measurement does not only depend on the numerical resolution, but also on the density gradient. In particular, this can limit the quality of measurements in the central regions of haloes. We propose a criterion to estimate the sensitivity of the measured shapes based on the density gradient of the halo and apply it to cosmological simulations of collisionless and self-interacting dark matter. By this, we demonstrate where reliable measurements of the halo shape are possible and how cored density profiles limit their applicability.

We study the effect of magnification in the Dark Energy Survey Year 3 analysis of galaxy clustering and galaxy-galaxy lensing, using two different lens samples: a sample of Luminous red galaxies, redMaGiC, and a sample with a redshift-dependent magnitude limit, MagLim. We account for the effect of magnification on both the flux and size selection of galaxies, accounting for systematic effects using the Balrog image simulations. We estimate the impact of magnification on the galaxy clustering and galaxy-galaxy lensing cosmology analysis, finding it to be a significant systematic for the MagLim sample. We show cosmological constraints from the galaxy clustering auto-correlation and galaxy-galaxy lensing signal with different magnifications priors, finding broad consistency in cosmological parameters in $\Lambda$CDM and $w$CDM. However, when magnification bias amplitude is allowed to be free, we find the two-point correlations functions prefer a different amplitude to the fiducial input derived from the image simulations. We validate the magnification analysis by comparing the cross-clustering between lens bins with the prediction from the baseline analysis, which uses only the auto-correlation of the lens bins, indicating systematics other than magnification may be the cause of the discrepancy. We show adding the cross-clustering between lens redshift bins to the fit significantly improves the constraints on lens magnification parameters and allows uninformative priors to be used on magnification coefficients, without any loss of constraining power or prior volume concerns.

Puffy disc is a numerical model, expected to capture the properties of the accretion flow in X-ray black hole binaries in the luminous, mildly sub-Eddington state. We fit the kerrbb and kynbb spectral models in XSPEC to synthetic spectra of puffy accretion discs, obtained in general relativistic radiative magnetohydrodynamic simulations, to see if they correctly recover the black hole spin and mass accretion rate assumed in the numerical simulation. We conclude that neither of the two models is capable of correctly interpreting the puffy disc parameters, which highlights a necessity to develop new, more accurate, spectral models for the luminous regime of accretion in X-ray black hole binaries. We propose that such spectral models should be based on the results of numerical simulations of accretion.

In star-forming galaxies, the far-infrared (FIR) and radio-continuum luminosities obey a tight empirical relation over a large range of star-formation rates (SFR). To understand the physics, we examine magnetohydrodynamic galaxy simulations, which follow the genesis of cosmic ray (CR) protons at supernovae and their advective and anisotropic diffusive transport. We show that gravitational collapse of the proto-galaxy generates a corrugated accretion shock, which injects turbulence and drives a small-scale magnetic dynamo. As the shock propagates outwards and the associated turbulence decays, the large velocity shear between the supersonically rotating cool disc with respect to the (partially) pressure-supported hot circumgalactic medium excites Kelvin-Helmholtz surface and body modes. Those interact non-linearly, inject additional turbulence and continuously drive multiple small-scale dynamos, which exponentially amplify weak seed magnetic fields. After saturation at small scales, they grow in scale to reach equipartition with thermal and CR energies in Milky Way-mass galaxies. In small galaxies, the magnetic energy saturates at the turbulent energy while it fails to reach equipartition with thermal and CR energies. We solve for steady-state spectra of CR protons, secondary electrons/positrons from hadronic CR-proton interactions with the interstellar medium, and primary shock-accelerated electrons at supernovae. The radio-synchrotron emission is dominated by primary electrons, irradiates the magnetized disc and bulge of our simulated Milky Way-mass galaxy and weakly traces bubble-shaped magnetically loaded outflows. Our star-forming and star-bursting galaxies with saturated magnetic fields match the global FIR-radio correlation (FRC) across four orders of magnitude. Its intrinsic scatter arises due to (i) different magnetic saturation levels that result from different seed magnetic fields, (ii) different radio synchrotron luminosities for different specific SFRs at fixed SFR, and (iii) a varying radio intensity with galactic inclination. In agreement with observations, several 100-pc-sized regions within star-forming galaxies also obey the FRC, while the centres of starbursts substantially exceed the FRC.

The integrated shear 3-point correlation function ζ_± is a higher-order statistic of the cosmic shear field that describes the modulation of the 2-point correlation function ξ_± by long-wavelength features in the field. Here, we introduce a new theoretical model to calculate ζ_± that is accurate on small angular scales, and that allows to take baryonic feedback effects into account. Our model builds on the realization that the small-scale ζ_± is dominated by the non-linear matter bispectrum in the squeezed limit, which can be evaluated accurately using the non-linear matter power spectrum and its first-order response functions to density and tidal field perturbations. We demonstrate the accuracy of our model by showing that it reproduces the small-scale ζ_± measured in simulated cosmic shear maps. The impact of baryonic feedback enters effectively only through the corresponding impact on the non-linear matter power spectrum, thereby permitting to account for these astrophysical effects on ζ_± similarly to how they are currently accounted for on ξ_±. Using a simple idealized Fisher matrix forecast for a DES-like survey we find that, compared to ξ_±, a combined $\xi _{\pm }\ \&\ \zeta _{\pm }$ analysis can lead to improvements of order $20\!-\!40{{\ \rm per\ cent}}$ on the constraints of cosmological parameters such as σ₈ or the dark energy equation of state parameter w₀. We find similar levels of improvement on the constraints of the baryonic feedback parameters, which strengthens the prospects for cosmic shear data to obtain tight constraints not only on cosmology but also on astrophysical feedback models. These encouraging results motivate future works on the integrated shear 3-point correlation function towards applications to real survey data.

We study the evolution of dust in a cosmological volume using a hydrodynamical simulation in which the dust production is coupled with the MUPPI (MUlti Phase Particle Integrator) sub-resolution model of star formation and feedback. As for the latter, we keep as reference the model setup calibrated previously to match the general properties of Milky Way-like galaxies in zoom-in simulations. However, we suggest that an increase of the star formation efficiency with the local dust-to-gas ratio would better reproduce the observed evolution of the cosmic star formation density. Moreover, the paucity of quenched galaxies at low redshift demands a stronger role of active galactic nucleus feedback. We tune the parameters ruling direct dust production from evolved stars and accretion in the interstellar medium to get scaling relations involving dust, stellar mass and metallicity in good agreement with observations. In low-mass galaxies, the accretion process is inefficient. As a consequence, they remain poorer in silicate and small grains than higher mass ones. We reproduce reasonably well the few available data on the radial distribution of dust outside the galactic region, supporting the assumption that the dust and gas dynamics are well coupled at galactic scales.

We study a neural network framework for the numerical evaluation of Feynman loop integrals that are fundamental building blocks for perturbative computations of physical observables in gauge and gravity theories. We show that such a machine learning approach improves the convergence of the Monte Carlo algorithm for high-precision evaluation of multi-dimensional integrals compared to traditional algorithms. In particular, we use a neural network to improve the importance sampling. For a set of representative integrals appearing in the computation of the conservative dynamics for a compact binary system in General Relativity, we perform a quantitative comparison between the Monte Carlo integrators VEGAS and i-flow, an integrator based on neural network sampling.

While the role of local interactions in nonequilibrium phase transitions is well studied, a fundamental understanding of the effects of long-range interactions is lacking. We study the critical dynamics of reproducing agents subject to auto-chemotactic interactions and limited resources. A renormalization group analysis reveals distinct scaling regimes for fast (attractive or repulsive) interactions; for slow signal transduction the dynamics is dominated by a diffusive fixed point. Further, we present a novel nonlinear mechanism that stabilizes the continuous transition against the emergence of a characteristic length scale due to a chemotactic collapse.

Disk winds are an important mechanism for accretion and disk evolution around young stars. The accreting intermediate-mass T-Tauri star RY Tau has an active jet and a previously known disk wind. Archival optical and new near-infrared observations of the RY Tau system show two horn-like components stretching out as a cone from RY Tau. Scattered light from the disk around RY Tau is visible in near-infrared but not seen at optical wavelengths. In the near-infrared, dark wedges that separates the horns from the disk, indicating we may see the scattered light from a disk wind. We use archived ALMA and SPHERE/ZIMPOL I-band observations combined with newly acquired SPEHRE/IRDIS H-band observations and available literature to build a simple geometric model of the RY Tau disk and disk wind. We use Monte Carlo radiative transfer modelling \textit{MCMax3D} to create comparable synthetic observations that test the effect of a dusty wind on the optical effect in the observations. We constrain the grain size and dust mass needed in the disk wind to reproduce the effect from the observations. A model geometrically reminiscent of a dusty disk wind with small micron to sub-micron size grains elevated above the disk can reproduce the optical effect seen in the observations. The mass in the obscuring component of the wind has been constrained to $1\times10^{-9} M_{\odot} \leq M \leq 5\times10^{-8} M_{\odot}$ which corresponds to a lower limit mass loss rate in the wind of about $\sim 1\times10^{-8}M_{\odot}\mathrm{yr}^{-1}$. While an illuminate dust cavity cannot be ruled out without measurements of the gas velocity, we argue that a magnetically launched disk wind is the most likely scenario.

Current best limits on the 21-cm signal during reionization are provided at large scales ($\gtrsim$100 Mpc). To model these scales, enormous simulation volumes are required which are computationally expensive. We find that the primary source of uncertainty at these large scales is sample variance, which decides the minimum size of simulations required to analyse current and upcoming observations. In large-scale structure simulations, the method of `fixing' the initial conditions (ICs) to exactly follow the initial power spectrum and `pairing' two simulations with exactly out-of-phase ICs has been shown to significantly reduce sample variance. Here we apply this `fixing and pairing' (F\&P) approach to reionization simulations whose clustering signal originates from both density fluctuations and reionization bubbles. Using a semi-numerical code, we show that with the traditional method, simulation boxes of $L\simeq 500$ (300) Mpc are required to model the large-scale clustering signal at $k$=0.1 Mpc$^{-1}$ with a precision of 5 (10) per cent. Using F\&P, the simulation boxes can be reduced by a factor of 2 to obtain the same precision level. We conclude that the computing costs can be reduced by at least a factor of 4 when using the F\&P approach.

To use strong gravitational lenses as an astrophysical or cosmological probe, models of their mass distributions are often needed. We present a new, time-efficient automation code for uniform modeling of strongly lensed quasars with GLEE, a lens modeling software, for high-resolution multi-band data. By using the observed positions of the lensed quasars and the spatially extended surface brightness distribution of the lensed quasar host galaxy, we obtain a model of the mass distribution of the lens galaxy. We apply this uniform modeling pipeline to a sample of nine strongly lensed quasars with HST WFC 3 images. The models show in most cases well reconstructed light components and a good alignment between mass and light centroids. We find that the automated modeling code significantly reduces the user input time during the modeling process. The preparation time of required input files is reduced significantly. This automated modeling pipeline can efficiently produce uniform models of extensive lens system samples which can be used for further cosmological analysis. A blind test through a comparison with the results of an independent automated modeling pipeline based on the modeling software Lenstronomy reveals important lessons. Quantities such as Einstein radius, astrometry, mass flattening and position angle are generally robustly determined. Other quantities depend crucially on the quality of the data and the accuracy of the PSF reconstruction. Better data and/or more detailed analysis will be necessary to elevate our automated models to cosmography grade. Nevertheless, our pipeline enables the quick selection of lenses for follow-up monitoring and further modeling, significantly speeding up the construction of cosmography-grade models. This is an important step forward to take advantage of the orders of magnitude increase in the number of lenses expected in the coming decade.

We provide improved Standard Model theory predictions for the exclusive rare semimuonic processes B → K^(*)μ⁺μ⁻ and B_s → ϕμ⁺μ⁻. Our results are based on a novel parametrization of the non-local form factors, which manifestly respects a recently developed dispersive bound. We critically compare our predictions to those obtained in the framework of QCD factorization. Our predictions provide, for the first time, parametric estimates of the systematic uncertainties due to non-local contributions. Comparing our predictions within the Standard Model to available experimental data, we find a large tension for B → Kμ⁺μ⁻. A simple model-independent analysis of potential effects beyond the Standard Model yields results compatible with other approaches, albeit with larger uncertainties for the B → K^*μ⁺μ⁻ and B_s → ϕμ⁺μ⁻ decays. Our approach yields systematically improvable predictions, and we look forward to its application in further analyses beyond the Standard Model.

In this work we report the realization of the first low-threshold cryogenic detector that uses diamond as absorber for astroparticle physics applications. We tested two 0.175 g CVD diamond samples, each instrumented with a W-TES. The sensors showed transitions at about 25 mK. We present the performance of the diamond detectors and we highlight the best performing one, where we obtained an energy threshold as low as 16.8 eV. This promising result lays the foundation for the use of diamond for different fields of applications where low threshold and excellent energy resolution are required, as i.e. light dark matter searches and BSM physics with coherent elastic neutrino nucleus scattering.

We present an analysis of the kinematics of the Radcliffe Wave, a 2.7 kpc long sinusoidal band of molecular clouds in the solar neighborhood recently detected via 3D dust mapping. With Gaia DR2 astrometry and spectroscopy, we analyze the 3D space velocities of ~1500 young stars along the Radcliffe Wave in action-angle space, using the motion of the wave's newly born stars as a proxy for its gas motion. We find that the vertical angle of young stars-corresponding to their orbital phase perpendicular to the Galactic plane-varies significantly as a function of position along the structure, in a pattern potentially consistent with a wavelike oscillation. This kind of oscillation is not seen in a control sample of older stars from Gaia occupying the same volume, disfavoring formation channels caused by long-lived physical processes. We use a "wavy midplane" model to try to account for the trend in vertical angles seen in young stars, and find that while the best-fit parameters for the wave's spatial period and amplitude are qualitatively consistent with the existing morphology defined by 3D dust, there is no evidence for additional velocity structure. These results support more recent and/or transitory processes in the formation of the Radcliffe Wave, which would primarily affect the motion of the wave's gaseous material. Comparisons of our results with new and upcoming simulations, in conjunction with new stellar radial velocity measurements in Gaia DR3, should allow us to further discriminate between various competing hypotheses.

The Standard Model (SM) does not contain by definition any new physics (NP) contributions to any observable but contains four CKM parameters which are not predicted by this model. We point out that if these four parameters are determined in a global fit that includes processes which are infected by NP, the resulting SM contributions to rare decay branching ratios cannot be considered as true SM contributions to the latter. On the other hand true SM predictions, that are free from the CKM dependence, can be obtained for suitable ratios of the $K$ and $B$ rare decay branching ratios to $\Delta M_s$, $\Delta M_d$ and $|\varepsilon_K|$, all calculated within the SM. These three observables contain by now only small hadronic uncertainties and are already well measured so that rather precise true SM predictions for the ratios in question can be obtained. In this context the rapid test of NP infection in the $\Delta F=2$ sector is provided by a $|V_{cb}|-\gamma$ plot that involves $\Delta M_s$, $\Delta M_d$, $|\varepsilon_K|$, and the mixing induced CP-asymmetry $S_{\psi K_S}$. As with the present hadronic matrix elements this test turns out to be negative, assuming negligible NP infection in the $\Delta F=2$ sector and setting the values of these four observables to the experimental ones, allows to obtain SM predictions for all $K$ and $B$ rare decay branching ratios that are most accurate to date and as a byproduct to obtain the full CKM matrix on the basis of $\Delta F=2$ transitions alone. Using this strategy we obtain SM predictions for 26 branching ratios for rare semileptonic and leptonic $K$ and $B$ decays with the $\mu^+\mu^-$ pair or the $\nu\bar\nu$ pair in the final state. Most interesting turn out to be the anomalies in the low $q^2$ bin in $B^+\to K^+\mu^+\mu^-$ ($5.1\sigma$) and $B_s\to \phi\mu^+\mu^-$ ($4.8\sigma$).

In this work we derive constraints on interacting dark matter-dark radiation models from a full-shape analysis of BOSS-DR12 galaxy clustering data, combined with Planck legacy cosmic microwave background (CMB) and baryon acoustic oscillation (BAO) measurements. We consider a set of models parameterized within the effective theory of structure formation (ETHOS), quantifying the lifting of the $S_8$ tension in view of KiDS weak-lensing results. The most favorable scenarios point to a fraction $f\sim 10-100\%$ of interacting dark matter as well as a dark radiation temperature that is smaller by a factor $\xi\sim 0.1-0.15$ compared to the CMB, leading to a reduction of the tension to the $\sim 1\sigma$ level. The temperature dependence of the interaction rate favored by relaxing the $S_8$ tension is realized for a weakly coupled unbroken non-Abelian $SU(N)$ gauge interaction in the dark sector. To map our results onto this $SU(N)$ model, we compute higher-order corrections due to Debye screening. We find a lower bound $\alpha_d\equiv g_d^2/(4\pi)\gtrsim 10^{-8} (10^{-9})$ for dark matter mass $1000 (1)$ GeV for relaxing the $S_8$ tension, consistent with upper bounds from galaxy ellipticities and compatible with self-interactions relevant for small-scale structure formation.

First-order cosmological phase transitions in the early Universe source sound waves and, subsequently, a background of stochastic gravitational waves. Currently, predictions of these gravitational waves rely heavily on simulations of a Higgs field coupled to the plasma of the early Universe, the former providing the latent heat of the phase transition. Numerically, this is a rather demanding task since several length scales enter the dynamics. From smallest to largest, these are the thickness of the Higgs interface separating the different phases, the shell thickness of the sound waves, and the average bubble size. In this work, we present an approach to perform Higgsless simulations in three dimensions, producing fully nonlinear results, while at the same time removing the hierarchically smallest scale from the lattice. This significantly reduces the complexity of the problem and contributes to making our approach highly efficient. We provide spectra for the produced gravitational waves for various choices of wall velocity and strength of the phase transition, as well as introduce a fitting function for the spectral shape.

Context. Recent observations found that observed cluster member galaxies are more compact than their counterparts in ΛCDM hydrodynamic simulations, as indicated by the difference in their strong gravitational lensing properties, and they reported that measured and simulated galaxy-galaxy strong lensing events on small scales are discrepant by one order of magnitude. Among the possible explanations for this discrepancy, some studies suggest that simulations with better resolution and implementing different schemes for galaxy formation could produce simulations that are in better agreement with the observations.
Aims: In this work, we aim to assess the impact of numerical resolution and of the implementation of energy input from AGN feedback models on the inner structure of cluster sub-haloes in hydrodynamic simulations.
Methods: We compared several zoom-in re-simulations of a sub-sample of cluster-sized haloes obtained by varying mass resolution and softening the length and AGN energy feedback scheme. We studied the impact of these different setups on the sub-halo (SH) abundances, their radial distribution, their density and mass profiles, and the relation between the maximum circular velocity, which is a proxy for SH compactness
Results: Regardless of the adopted numerical resolution and feedback model, SHs with masses of M_SH ≲ 10¹¹ h⁻¹ M_⊙, the most relevant mass range for galaxy-galaxy strong lensing, have maximum circular velocities ∼30% smaller than those measured from strong lensing observations. We also find that simulations with less effective AGN energy feedback produce massive SHs (M_SH ≳ 10¹¹ h⁻¹ M_⊙) with higher maximum circular velocity and that their V_max − M_SH relation approaches the observed one. However, the stellar-mass number count of these objects exceeds the one found in observations, and we find that the compactness of these simulated SHs is the result of an extremely over-efficient star formation in their cores, also leading to larger than observed SH stellar mass.
Conclusions: Regardless of the resolution and galaxy formation model adopted, simulations are unable to simultaneously reproduce the observed stellar masses and compactness (or maximum circular velocities) of cluster galaxies. Thus, the discrepancy between theory and observations that emerged previous works. It remains an open question as to whether such a discrepancy reflects limitations of the current implementation of galaxy formation models or the ΛCDM paradigm.

The high-$x$ data from the ZEUS Collaboration are used to extract parton density distributions of the proton deep in the perturbative regime of QCD. The data primarily constrain the up-quark valence distribution and new results are presented on its $x$-dependence as well as on the momentum carried by the up-quark. The results were obtained using Bayesian analysis methods which can serve as a model for future parton density extractions.

The dividing line between galaxies that are quenched by reionization ("relics") and galaxies that survive reionization (i.e. continue forming stars) is commonly discussed in terms of a halo mass threshold. We probe this threshold in a physically more complete and accurate way than has been possible to date, using five extremely high resolution ($M_\mathrm{target}=4M_\odot$) cosmological zoom-in simulations of dwarf galaxies within the halo mass range $1-4\times10^9M_\odot$. The employed LYRA simulation model features resolved interstellar medium physics and individual, resolved supernova explosions. In our results, we discover an interesting intermediate population of dwarf galaxies close to the threshold mass but which are neither full reionization relics nor full reionization survivors. These galaxies initially quench at the time of reionization but merely remain quiescent for ~500Myr. At $z\approx5$ they recommence star formation in a synchronous way, and remain star-forming until the present day. These results demonstrate that the halo mass at $z=0$ is not a good indicator of survival close to the threshold. While the star formation histories we find are diverse, we show that they are directly related to the ability of a given halo to retain and cool gas. Whereas the latter is most strongly dependent on the mass (or virial temperature) of the host halo at the time of reionization, it also depends on its growth history, the UV background (and its decrease at late times) and the amount of metals retained within the halo.

We use $N$-body simulations to study halo assembly bias (i.e., the dependence of halo clustering on properties beyond total mass) in the density and primordial non-Gaussianity (PNG) linear bias parameters $b_1$ and $b_\phi$, respectively. We consider concentration, spin and sphericity as secondary halo properties, for which we find a clear detection of assembly bias for $b_1$ and $b_\phi$. At fixed total mass, halo spin and sphericity impact $b_1$ and $b_\phi$ in a similar manner, roughly preserving the shape of the linear $b_\phi(b_1)$ relation satisfied by the global halo population. Halo concentration, however, drives $b_1$ and $b_\phi$ in opposite directions. This induces significant changes to the $b_\phi(b_1)$ relation, with higher concentration halos having higher amplitude of $b_\phi(b_1)$. For $z=0.5$ and $b_1 \approx 2$ in particular, the population comprising either all halos, those with the $33\%$ lowest or those with the $33\%$ highest concentrations have a PNG bias of $b_\phi \approx 3$, $b_\phi \approx -1$ and $b_\phi \approx 9$, respectively. Varying the halo concentration can make $b_\phi$ very small and even change its sign. These results have important ramifications for galaxy clustering constraints of the local PNG parameter $f_{\rm NL}$ that assume fixed forms for the $b_\phi(b_1)$ relation. We illustrate the significant impact of halo assembly bias in actual data using the BOSS DR12 galaxy power spectrum: assuming that BOSS galaxies are representative of all halos, the $33\%$ lowest or the $33\%$ highest concentration halos yields $\sigma_{f_{\rm NL}} = 44, 165, 19$, respectively. Our results suggest taking host halo concentration into account in galaxy selection strategies to maximize the signal-to-noise on $f_{\rm NL}$. They also motivate more simulation-based efforts to study the $b_\phi(b_1)$ relation of halos and galaxies.

The experimental detection of the CEν NS allows the investigation of neutrinos and neutrino sources with all-flavor sensitivity. Given its large content in neutrons and stability, Pb is a very appealing choice as target element. The presence of the radioisotope 21⁰Pb (T1_/2∼ 22 yrs) makes natural Pb unsuitable for low-background, low-energy event searches. This limitation can be overcome employing Pb of archaeological origin, where several half-lives of 21⁰Pb have gone by. We present results of a cryogenic measurement of a 15 g PbWO₄ crystal, grown with archaeological Pb (older than ∼2000 yrs) that achieved a sub-keV nuclear recoil detection threshold. A ton-scale experiment employing such material, with a detection threshold for nuclear recoils of just 1 keV would probe the entire Milky Way for SuperNovae, with equal sensitivity for all neutrino flavors, allowing the study of the core of such exceptional events.

Substructures are known to be good tracers for the dynamical states and recent accretion histories of the most massive collapsed structures in the Universe, galaxy clusters. Observations find extremely massive substructures in some clusters, especially Abell 2744, which are potentially in tension with the $\Lambda$CDM paradigm since they are not found in simulations directly. However, the methods to measure substructure masses strongly differ between observations and simulations. Using the fully hydrodynamical cosmological simulation suite Magneticum Pathfinder we develop a method to measure substructure masses in projection from simulations, similar to the observational approach. We identify a simulated Abell 2744 counterpart that not only has eight substructures of similar mass fractions but also exhibits similar features in the hot gas component. This cluster formed only recently through a major merger together with at least 6 massive minor merger events since z=1, where prior the most massive component had a mass of less than $1\times10^{14}M_\odot$. We show that the mass fraction of all substructures and of the eighth substructure separately are excellent tracers for the dynamical state and assembly history for all galaxy cluster mass ranges, with high fractions indicating merger events within the last 2Gyr. Finally, we demonstrate that the differences between subhalo masses measured directly from simulations as bound and those measured in projection are due to methodology, with the latter generally 2-3 times larger than the former. We provide a predictor function to estimate projected substructure masses from SubFind masses for future comparison studies between simulations and observations.

We present MUSE spectroscopy, Megacam imaging, and Chandra X-ray emission for SPT-CL J0307-6225, a z=0.58 major merging galaxy cluster with a large BCG-SZ centroid separation and a highly disturbed X-ray morphology. The galaxy density distribution shows two main overdensities with separations of 0.144 and 0.017 arcmin to their respective BCGs. We characterize the central regions of the two colliding structures, namely 0307-6225N and 0307-6225S, finding velocity derived masses of M_{200, N} = 2.44 ± 1.41 × 10¹⁴ M_⊙ and M_{200, S} = 3.16 ± 1.88 × 10¹⁴ M_⊙, with a line-of-sight velocity difference of |Δv| = 342 km s^-1. The total dynamically derived mass is consistent with the SZ derived mass of 7.63 h$_{70}^{-1}$ ± 1.36 × 10¹⁴ M_⊙. We model the merger using the Monte Carlo Merger Analysis Code, estimating a merging angle of 36$^{+14}_{-12}$ degrees with respect to the plane of the sky. Comparing with simulations of a merging system with a mass ratio of 1:3, we find that the best scenario is that of an ongoing merger that began 0.96$^{+0.31}_{-0.18}$ Gyr ago. We also characterize the galaxy population using Hδ and [OII] λ3727 Å lines. We find that most of the emission-line galaxies belong to 0307-6225S, close to the X-ray peak position, with a third of them corresponding to red-cluster sequence galaxies, and the rest to blue galaxies with velocities consistent with recent periods of accretion. Moreover, we suggest that 0307-6225S suffered a previous merger, evidenced through the two equally bright BCGs at the center with a velocity difference of ~674 km s^-1.

We introduce the star formation and Supernova (SN) feedback model of the SATIN (Simulating AGNs Through ISM with Non-Equilibrium Effects) project to simulate the evolution of the star forming multi-phase interstellar medium (ISM) of entire disk galaxies. This galaxy-wide implementation of a successful ISM feedback model naturally covers an order of magnitude in gas surface density, shear and radial motions. It is implemented in the adaptive mesh refinement code RAMSES at a peak resolution of 9 pc. New stars are represented by star cluster (sink) particles with individual SN delay times for massive stars. With SN feedback, cooling and gravity, the galactic ISM develops a realistic three-phase structure. The star formation rates naturally follow observed scaling relations for the local Milky Way gas surface density. SNe drive additional turbulence in the warm (300 K < $T$ < 10$^4$ K) gas and increase the kinetic energy of the cold gas, cooling out of the warm phase. The majority of the gas leaving the galactic ISM is warm and hot with mass loading factors of $3 \le \eta \le 10$. While the hot gas is leaving the system, the warm and cold gas falls back onto the disc in a galactic fountain flow.

Shape measurements of galaxies and galaxy clusters are widespread in the analysis of cosmological simulations. But the limitations of those measurements have been poorly investigated. In this paper, we explain why the quality of the shape measurement does not only depend on the numerical resolution, but also on the density gradient. In particular, this can limit the quality of measurements in the central regions of haloes. We propose a criterion to estimate the sensitivity of the measured shapes based on the density gradient of the halo and apply it to cosmological simulations of collisionless and self-interacting dark matter. By this, we demonstrate where reliable measurements of the halo shape are possible and how cored density profiles limit their applicability.

Thermal bombs are a widely used method to artificially trigger explosions of core-collapse supernovae (CCSNe) to determine their nucleosynthesis or ejecta and remnant properties. Recently, their use in spherically symmetric (1D) hydrodynamic simulations led to the result that {56,57}Ni and 44Ti are massively underproduced compared to observational estimates for Supernova 1987A, if the explosions are slow, i.e., if the explosion mechanism of CCSNe releases the explosion energy on long timescales. It was concluded that rapid explosions are required to match observed abundances, i.e., the explosion mechanism must provide the CCSN energy nearly instantaneously on timescales of some ten to order 100 ms. This result, if valid, would disfavor the neutrino-heating mechanism, which releases the CCSN energy on timescales of seconds. Here, we demonstrate by 1D hydrodynamic simulations and nucleosynthetic post-processing that these conclusions are a consequence of disregarding the initial collapse of the stellar core in the thermal-bomb modelling before the bomb releases the explosion energy. We demonstrate that the anti-correlation of 56Ni yield and energy-injection timescale vanishes when the initial collapse is included and that it can even be reversed, i.e., more 56Ni is made by slower explosions, when the collapse proceeds to small radii similar to those where neutrino heating takes place in CCSNe. We also show that the 56Ni production in thermal-bomb explosions is sensitive to the chosen mass cut and that a fixed mass layer or fixed volume for the energy deposition cause only secondary differences. Moreover, we propose a most appropriate setup for thermal bombs.

We present the first results of a comprehensive supernova (SN) radiative-transfer (RT) code-comparison initiative (StaNdaRT), where the emission from the same set of standardized test models is simulated by currently-used RT codes. A total of ten codes have been run on a set of four benchmark ejecta models of Type Ia supernovae. We consider two sub-Chandrasekhar-mass ($M_\mathrm{tot} = 1.0$ M$_\odot$) toy models with analytic density and composition profiles and two Chandrasekhar-mass delayed-detonation models that are outcomes of hydrodynamical simulations. We adopt spherical symmetry for all four models. The results of the different codes, including the light curves, spectra, and the evolution of several physical properties as a function of radius and time, are provided in electronic form in a standard format via a public repository. We also include the detailed test model profiles and several python scripts for accessing and presenting the input and output files. We also provide the code used to generate the toy models studied here. In this paper, we describe in detail the test models, radiative-transfer codes and output formats and provide access to the repository. We present example results of several key diagnostic features.

LiteBIRD is a JAXA-led strategic large-class satellite mission designed to measure the polarization of the cosmic microwave background and Galactic foregrounds from 34 to 448 GHz across the entire sky from L2 in the late 2020s. The scientific payload includes three telescopes which are called the low-, mid-, and high-frequency telescopes each with their own receiver that covers a portion of the mission's frequency range. The low frequency telescope will map synchrotron radiation from the Galactic foreground and the cosmic microwave background. We discuss the design, fabrication, and characterization of the low-frequency focal plane modules for low-frequency telescope, which has a total bandwidth ranging from 34 to 161 GHz. There will be a total of 4 different pixel types with 8 overlapping bands to cover the full frequency range. These modules are housed in a single low-frequency focal plane unit which provides thermal isolation, mechanical support, and radiative baffling for the detectors. The module design implements multi-chroic lenslet-coupled sinuous antenna arrays coupled to transition edge sensor bolometers read out with frequency-domain mulitplexing. While this technology has strong heritage in ground-based cosmic microwave background experiments, the broad frequency coverage, low optical loading conditions, and the high cosmic ray background of the space environment require further development of this technology to be suitable for LiteBIRD. In these proceedings, we discuss the optical and bolometeric characterization of a triplexing prototype pixel with bands centered on 78, 100, and 140 GHz.

Over the last years, large (sub-)millimetre surveys of protoplanetary disks have well constrained the demographics of disks, such as their millimetre luminosities, spectral indices, and disk radii. Additionally, several high-resolution observations have revealed an abundance of substructures in the disks dust continuum. The most prominent are ring like structures, likely due to pressure bumps trapping dust particles. The origins and characteristics of these bumps, nevertheless, need to be further investigated. The purpose of this work is to study how dynamic pressure bumps affect observational properties of protoplanetary disks. We further aim to differentiate between the planetary- versus zonal flow-origin of pressure bumps. We perform one-dimensional gas and dust evolution simulations, setting up models with varying pressure bump features. We subsequently run radiative transfer calculations to obtain synthetic images and the different quantities of observations. We find that the outermost pressure bump determines the disks dust size across different millimetre wavelengths. Our modelled dust traps need to form early (< 0.1 Myr), fast (on viscous timescales), and must be long lived (> Myr) to obtain the observed high millimetre luminosities and low spectral indices of disks. While the planetary bump models can reproduce these observables irrespectively of the opacity prescription, the highest opacities are needed for the zonal flow bump model to be in line with observations. Our findings favour the planetary- over the zonal flow-origin of pressure bumps and support the idea that planet formation already occurs in early class 0-1 stages of circumstellar disks. The determination of the disks effective size through its outermost pressure bump also delivers a possible answer to why disks in recent low-resolution surveys appear to have the same sizes across different millimetre wavelengths.

The emergence of functional oligonucleotides on early Earth required a molecular selection mechanism to screen for specific sequences with prebiotic functions. Cyclic processes such as daily temperature oscillations were ubiquitous in this environment and could trigger oligonucleotide phase separation. Here, we propose sequence selection based on phase separation cycles realized through sedimentation in a system subjected to the feeding of oligonucleotides. Using theory and experiments with DNA, we show sequence-specific enrichment in the sedimented dense phase, in particular of short 22-mer DNA sequences. The underlying mechanism selects for complementarity, as it enriches sequences that tightly interact in the condensed phase through base-pairing. Our mechanism also enables initially weakly biased pools to enhance their sequence bias or to replace the most abundant sequences as the cycles progress. Our findings provide an example of a selection mechanism that may have eased screening for the first auto-catalytic self-replicating oligonucleotides.

In the context of false vacuum decay at zero temperature, it is well known that bubbles expand with uniform acceleration in the rest frame of nucleation. We show that this uniformly accelerating expansion suffers from an instability. This can be observed as a tachyonic mode in the spectrum of fluctuations for the energy functional in the reference frame in which the uniformly accelerating bubble wall appears static. In such a frame, arbitrary small perturbations cause an amplifying departure from the static wall solution. This implies that the nucleated bubble is not a critical point of the energy functional in the rest frame of nucleation but becomes one in the accelerating frame. The aforementioned instability for vacuum bubbles can be related to the well-known instability for the nucleated critical static bubbles during finite-temperature phase transitions in the rest frame of the plasma. It is proposed that zero-temperature vacuum decays as seen from accelerating frames have a dual description in terms of finite-temperature phase transitions.

We explore a novel approach to compute the force between a static quark-antiquark pair with the gradient flow algorithm on the lattice. The approach is based on inserting a chromoelectric field in a Wilson loop. The renormalization issues, associated with the finite size of the chromoelectric field on the lattice, can be solved with the use of gradient flow. We compare numerical results for the flowed static potential to our previous measurement of the same observable without a gradient flow.

We present TimeEvolver, a program for computing time evolution in a generic quantum system. It relies on well-known Krylov subspace techniques to tackle the problem of multiplying the exponential of a large sparse matrix iH, where H is the Hamiltonian, with an initial vector v. The fact that H is Hermitian makes it possible to provide an easily computable bound on the accuracy of the Krylov approximation. Apart from effects of numerical roundoff, the resulting a posteriori error bound is rigorous, which represents a crucial novelty as compared to existing software packages such as Expokit[1]. On a standard notebook, TimeEvolver allows to compute time evolution with adjustable precision in Hilbert spaces of dimension greater than 10⁶. Additionally, we provide routines for deriving the matrix H from a more abstract representation of the Hamiltonian operator.

We present a photometric analysis of star and star cluster (SC) formation in a high-resolution simulation of a dwarf galaxy starburst that allows the formation of individual stars to be followed. Previous work demonstrated that the properties of the SCs formed in the simulation are in good agreement with observations. In this paper, we create mock spectral energy distributions and broad-band photometric images using the radiative transfer code SKIRT 9. We test several observational star formation rate (SFR) tracers and find that 24 $\mu$m, total infrared and Hα trace the underlying SFR during the (post)starburst phase, while UV tracers yield a more accurate picture of star formation during quiescent phases prior to and after the merger. We then place the simulated galaxy at distances of 10 and 50 Mpc and use aperture photometry at Hubble Space Telescope resolution to analyse the simulated SC population. During the starburst phase, a hierarchically forming set of SCs leads inaccurate source separation because of crowding. This results in estimated SC mass function slopes that are up to ~0.3 shallower than the true slope of ~-1.9 to -2 found for the bound clusters identified from the particle data in the simulation. The masses of the largest clusters are overestimated by a factor of up to 2.9 due to unresolved clusters within the apertures. The aperture-based analysis also produces a relation between cluster formation efficiency and SFR surface density that is slightly flatter than that recovered from bound clusters. The differences are strongest in quiescent SF environments.

SN 2020cxd is a representative of the family of low-energy, underluminous Type IIP supernovae (SNe), whose observations and analysis were recently reported by Yang et al. (2021). Here, we re-evaluate the observational data for the diagnostic SN properties by employing the hydrodynamic explosion model of a 9 M_⊙ red supergiant progenitor with an iron core and a pre-collapse mass of 8.75 M_⊙. The explosion of the star was obtained by the neutrino-driven mechanism in a fully self-consistent simulation in three dimensions (3D). Multiband light curves and photospheric velocities for the plateau phase are computed with the one-dimensional radiation-hydrodynamics code STELLA, applied to the spherically averaged 3D explosion model as well as sphericized radial profiles in different directions of the 3D model. We find that the overall evolution of the bolometric light curve, duration of the plateau phase, and basic properties of the multiband emission can be well reproduced by our SN model with its explosion energy of only 0.7 × 10⁵⁰ erg and an ejecta mass of 7.4 M_⊙. These values are considerably lower than the previously reported numbers, but they are compatible with those needed to explain the fundamental observational properties of the prototype low-luminosity SN 2005cs. Because of the good compatibility of our photospheric velocities with line velocities determined for SN 2005cs, we conclude that the line velocities of SN 2020cxd are probably overestimated by up to a factor of about 3. The evolution of the line velocities of SN 2005cs compared to photospheric velocities in different explosion directions might point to intrinsic asymmetries in the SN ejecta.

LiteBIRD is a JAXA-led strategic large-class satellite mission designed to measure the polarization of the cosmic microwave background and Galactic foregrounds from 34 to 448 GHz across the entire sky from L2 in the late 2020s. The scientific payload includes three telescopes which are called the low-, mid-, and high-frequency telescopes each with their own receiver that covers a portion of the mission's frequency range. The low frequency telescope will map synchrotron radiation from the Galactic foreground and the cosmic microwave background. We discuss the design, fabrication, and characterization of the low-frequency focal plane modules for low-frequency telescope, which has a total bandwidth ranging from 34 to 161 GHz. There will be a total of 4 different pixel types with 8 overlapping bands to cover the full frequency range. These modules are housed in a single low-frequency focal plane unit which provides thermal isolation, mechanical support, and radiative baffling for the detectors. The module design implements multi-chroic lenslet-coupled sinuous antenna arrays coupled to transition edge sensor bolometers read out with frequency-domain mulitplexing. While this technology has strong heritage in ground-based cosmic microwave background experiments, the broad frequency coverage, low optical loading conditions, and the high cosmic ray background of the space environment require further development of this technology to be suitable for LiteBIRD. In these proceedings, we discuss the optical and bolometeric characterization of a triplexing prototype pixel with bands centered on 78, 100, and 140 GHz.

We investigate the formation and evolution of "primordial" dusty rings occurring in the inner regions of protoplanetary discs, with the help of long-term, coupled dust-gas, magnetohydrodynamic simulations. The simulations are global and start from the collapse phase of the parent cloud core, while the dead zone is calculated via an adaptive $\alpha$ formulation by taking into account the local ionization balance. The evolution of the dusty component includes its growth and back reaction on to the gas. Previously, using simulations with only a gas component, we showed that dynamical rings form at the inner edge of the dead zone. We find that when dust evolution as well as magnetic field evolution in the flux-freezing limit are included, the dusty rings formed are more numerous and span a larger radial extent in the inner disc, while the dead zone is more robust and persists for a much longer time. We show that these dynamical rings concentrate enough dust mass to become streaming unstable, which should result in rapid planetesimal formation even in the embedded phases of the system. The episodic outbursts caused by the magnetorotational instability have significant impact on the evolution of the rings. The outbursts drain the inner disc of grown dust, however, the period between bursts is sufficiently long for the planetesimal growth via streaming instability.The dust mass contained within the rings is large enough to ultimately produce planetary systems with the core accretion scenario. The low mass systems rarely undergo outbursts and thus, the conditions around such stars can be especially conducive for planet formation.

We present extended Lyα emission out to 800 kpc of 1034 [O III]-selected galaxies at redshifts 1.9 < z < 2.35 using the Hobby-Eberly Telescope Dark Energy Experiment. The locations and redshifts of the galaxies are taken from the 3D-HST survey. The median-stacked surface brightness profile of the Lyα emission of the [O III]-selected galaxies agrees well with that of 968 bright Lyα-emitting galaxies (LAEs) at r > 40 kpc from the galaxy centers. The surface brightness in the inner parts (r < 10 kpc) around the [O III]-selected galaxies, however, is 10 times fainter than that of the LAEs. Our results are consistent with the notion that photons dominating the outer regions of the Lyα halos are not produced in the central galaxies but originate outside of them.

We present [C II] synthetic observations of smoothed particle hydrodynamics (SPH) simulations of a dwarf galaxy merger. The merging process varies the star formation rate (SFR) by more than three orders of magnitude. Several star clusters are formed, the feedback of which disperses and unbinds the dense gas through expanding H II regions and supernova (SN) explosions. For galaxies with properties similar to the modeled ones, we find that the [C II] emission remains optically thin throughout the merging process. We identify the warm neutral medium ( $3\lt \mathrm{log}{T}_{\mathrm{gas}}\lt 4$ with χ _HI > 2χ _H2) to be the primary source of [C II] emission (~58% contribution), although at stages when the H II regions are young and dense (during star cluster formation or SNe in the form of ionized bubbles), they can contribute ≳50% to the total [C II] emission. We find that the [C II]/far-IR (FIR) ratio decreases owing to thermal saturation of the [C II] emission caused by strong far-UV radiation fields emitted by the massive star clusters, leading to a [C II] deficit medium. We investigate the [C II]-SFR relation and find an approximately linear correlation that agrees well with observations, particularly those from the Dwarf Galaxy Survey. Our simulation reproduces the observed trends of [C II]/FIR versus Σ_SFR and Σ_FIR, and it agrees well with the Kennicutt relation of SFR-FIR luminosity. We propose that local peaks of [C II] in resolved observations may provide evidence for ongoing massive cluster formation.

Cosmological shock waves are essential to understanding the formation of cosmological structures. To study them, scientists run computationally expensive high-resolution 3D hydrodynamic simulations. Interpreting the simulation results is challenging because the resulting data sets are enormous, and the shock wave surfaces are hard to separate and classify due to their complex morphologies and multiple shock fronts intersecting. We introduce a novel pipeline, Virgo, combining physical motivation, scalability, and probabilistic robustness to tackle this unsolved unsupervised classification problem. To this end, we employ kernel principal component analysis with low-rank matrix approximations to denoise data sets of shocked particles and create labeled subsets. We perform supervised classification to recover full data resolution with stochastic variational deep kernel learning. We evaluate on three state-of-the-art data sets with varying complexity and achieve good results. The proposed pipeline runs automatically, has only a few hyperparameters, and performs well on all tested data sets. Our results are promising for large-scale applications, and we highlight now enabled future scientific work.

RES-NOVA is a newly proposed experiment for detecting neutrinos from astrophysical sources, mainly Supernovae, using an array of cryogenic detectors made of PbWO₄ crystals produced from archaeological Pb. This unconventional material, characterized by intrinsic high radiopurity, enables low-background levels in the region of interest for the neutrino detection via Coherent Elastic neutrino-Nucleus Scattering (CEν NS). This signal lies at the detector energy threshold, O(1 keV), and it is expected to be hidden by naturally occurring radioactive contaminants of the crystal absorber. Here, we present the results of a radiopurity assay on a 0.84 kg PbWO₄ crystal produced from archaeological Pb operated as a cryogenic detector. The crystal internal radioactive contaminations are: 23²Th <40 μ Bq/kg, 23⁸U <30 μ Bq/kg, 22⁶Ra 1.3 mBq/kg and 21⁰Pb 22.5 mBq/kg. We also present a background projection for the final experiment and possible mitigation strategies for further background suppression. The achieved results demonstrate the feasibility of realizing this new class of detectors.

The understanding of the accretion process has a central role in the understanding of star and planet formation. We aim to test how accretion variability influences previous correlation analyses of the relation between X-ray activity and accretion rates, which is important for understanding the evolution of circumstellar disks and disk photoevaporation. We monitored accreting stars in the Orion Nebula Cluster from November 24, 2014, until February 17, 2019, for 42 epochs with the Wendelstein Wide Field Imager in the Sloan Digital Sky Survey u'g'r' filters on the 2 m Fraunhofer Telescope on Mount Wendelstein. Mass accretion rates were determined from the measured ultraviolet excess. The influence of the mass accretion rate variability on the relation between X-ray luminosities and mass accretion rates was analyzed statistically. We find a typical interquartile range of ~ 0.3 dex for the mass accretion rate variability on timescales from weeks to ~ 2 years. The variability has likely no significant influence on a correlation analysis of the X-ray luminosity and the mass accretion rate observed at different times when the sample size is large enough. The observed anticorrelation between the X-ray luminosity and the mass accretion rate predicted by models of photoevaporation-starved accretion is likely not due to a bias introduced by different observing times.

Elliptic modular graph forms (eMGFs) are non-holomorphic modular forms depending on a modular parameter $\tau$ of a torus and marked points $z$ thereon. Traditionally, eMGFs are constructed from nested lattice sums over the discrete momenta on the worldsheet torus in closed-string genus-one amplitudes. In this work, we develop methods to translate the lattice-sum realization of eMGFs into iterated integrals over modular parameters $\tau$ of the torus with particular focus on cases with one marked point. Such iterated-integral representations manifest algebraic and differential relations among eMGFs and their degeneration limit $\tau \rightarrow i\infty$. From a mathematical point of view, our results yield concrete realizations of single-valued elliptic polylogarithms at arbitrary depth in terms of meromorphic iterated integrals over modular forms and their complex conjugates. The basis dimensions of eMGFs at fixed modular and transcendental weights are derived from a simple counting of iterated integrals and a generalization of Tsunogai's derivation algebra.

We describe a systematic approach for the evaluation of Witten diagrams for multi-loop scattering amplitudes of a conformally coupled scalar ϕ⁴-theory in Euclidean AdS₄, by recasting the Witten diagrams as flat space Feynman integrals. We derive closed form expressions for the anomalous dimensions for all double-trace operators up to the second order in the coupling constant. We explain the relation between the flat space unitarity methods and the discontinuities of the short distance expansion on the boundary of Witten diagrams.

We present a general framework for calculating post-Minskowskian, classical, conservative Hamiltonians for $N$ non-spinning bodies in general relativity from relativistic scattering amplitudes. Novel features for $N>2$ are described including the subtraction of tree-like iteration contributions and the calculation of non-trivial many-body Fourier transform integrals needed to construct position space potentials. A new approach to calculating these integrals as an expansion in the hierarchical limit is described based on the method of regions. As an explicit example, we present the $\mathcal{O}\left(G^2\right)$ 3-body momentum space potential in general relativity as well as for charged bodies in Einstein-Maxwell. The result is shown to be in perfect agreement with previous post-Newtonian calculations in general relativity up to $\mathcal{O}\left(G^2 v^4\right)$. Furthermore, in appropriate limits the result is shown to agree perfectly with relativistic probe scattering in multi-center extremal black hole backgrounds and with the scattering of slowly-moving extremal black holes in the moduli space approximation.

We present a simple parton-shower model that replaces the explicit angular ordering of the coherent branching formalism with a differentially accurate simulation of soft-gluon radiation by means of a non-trivial dependence on azimuthal angles. We introduce a global kinematics mapping and provide an analytic proof that it satisfies the criteria for next-to leading logarithmic accuracy. In the new algorithm, initial and final state evolution are treated on the same footing. We provide an implementation for final-state evolution in the numerical code Alaric and present a first comparison to experimental data.

We present for the first time Next-to-Leading (NLO) QCD renormalization group (RG) evolution matrices for nonleptonic Δ F =2 transitions in the Standard Model effective field theory (SMEFT). To this end we transform first the known two-loop QCD anomalous dimension matrices (ADMs) of the BSM (Beyond the SM) operators in the so-called Buras Misiak Urban basis into the ones in the common weak effective theory (WET) basis (the so-called Jenkins Manohar Stoffer basis) for which tree-level and one-loop matching to the SMEFT are already known. This subsequently allows us to find the two-loop QCD ADMs for the SMEFT nonleptonic Δ F =2 operators in the Warsaw basis. Having all these ingredients we investigate the impact of these NLO QCD effects on the QCD RG evolution of SMEFT Wilson coefficients for nonleptonic Δ F =2 transitions from the new physics scale Λ down to the electroweak scale μ_ew. The main benefit of these new contributions is that they allow one to remove renormalization scheme dependences present in the one-loop matchings both between the WET and SMEFT and also between SMEFT and a chosen UV completion. But the Next-to-Leading (NLO) QCD effects, calculated here in the Naive dimensional regularisation minimal subtraction scheme, turn out to be small, in the ballpark of a few percent but larger than one-loop Yukawa top effects when only the Δ F =2 operators are considered. The more complicated class of nonleptonic Δ F =1 decays will be presented soon in another publication.

Planes of satellites are observed around many galaxies. However, these observations are still considered a point of tension for the $\Lambda$CDM paradigm. We use the fully hydrodynamical cosmological $\Lambda$CDM state-of-the-art simulation Magneticum Pathfinder to investigate the existence of such planes over a large range of haloes, from Milky Way to galaxy cluster masses. To this end, we develop the Momentum in Thinnest Plane (MTP) method to identify planes and quantify the properties of their constituent satellites, considering both position and momentum. We find that thin planes ($20\%\cdot R_\mathrm{halo}$) containing at least $50\%$ of the total number of satellites can be found in almost all systems. In Milky Way mass-like systems, around 86\% of such planes are even aligned in momentum ($90\%$ of the total satellite momentum), where the fraction is smaller if more satellites are required to be inside the plane. We further find a mass dependency, with more massive systems exhibiting systematically thicker planes. This may point towards the change from continuous accretion of small objects along filaments and sheets for less massive haloes to the accretion of large objects (e.g., major mergers) dominating the growth of more massive haloes. There is no correlation between the existence of a plane and the main galaxy's morphology. Finally, we find a clear preference for the minor axes of the satellite planes and the host galaxy to be aligned, in agreement with recent observations.

Event reconstruction is a central step in many particle physics experiments, turning detector observables into parameter estimates; for example estimating the energy of an interaction given the sensor readout of a detector. A corresponding likelihood function is often intractable, and approximations need to be constructed. In our work, we first show how the full likelihood for a many-sensor detector can be broken apart into smaller terms, and secondly how we can train neural networks to approximate all terms solely based on forward simulation. Our technique results in a fast, flexible, and close-to-optimal surrogate model proportional to the likelihood and can be used in conjunction with standard inference techniques. We illustrate our technique for parameter inference in neutrino telescopes based on maximum likelihood and Bayesian posterior sampling. Given its great flexibility, we also showcase our method for detector optimization, and apply it to simulation of a ton-scale water-based liquid scintillator detector.

We update the Standard Model (SM) predictions for $B$-meson lifetimes within the heavy quark expansion (HQE). Including for the first time the contribution of the Darwin operator, SU(3)$_F$ breaking corrections to the matrix element of dimension-six four-quark operators and the so-called eye-contractions, we obtain for the total widths $\Gamma (B^+) = (0.58^{+0.11}_{-0.07}) \, \mbox{ps}^{-1}$, $\Gamma (B_d) = (0.63^{+0.11}_{-0.07}) \, \mbox{ps}^{-1}$, $\Gamma (B_s) = (0.63^{+0.11}_{-0.07}) \, \mbox{ps}^{-1}$, and for the lifetime ratios $\tau (B^+) / \tau (B_d) = 1.086 \pm 0.022$, $\tau (B_s) / \tau (B_d) = 1.003 \pm 0.006 \, (1.028 \pm 0.011)$. The two values for the last observable arise from using two different sets of input for the non-perturbative parameters $\mu_\pi^2(B_d)$, $\mu_G^2(B_d)$, and $\rho_D^3(B_d)$ as well as from different estimates of the SU(3)$_F$ breaking in these parameters. Our results are overall in very good agreement with the corresponding experimental data, however, there seems to emerge a tension in $\tau (B_s)/\tau (B_d)$ when considering the second set of input parameters. Specifically, this observable is extremely sensitive to the size of the parameter $\rho_D^3 (B_d)$ and of the SU(3)$_F$ breaking effects in $\mu_\pi^2$, $\mu_G^2$ and $\rho_D^3$; hence, it is of utmost importance to be able to better constrain all these parameters. In this respect, an extraction of $\mu_\pi^2 (B_s), \mu_G^2 (B_s), \rho_D^3 (B_s)$ from future experimental data on inclusive semileptonic $B_s$-meson decays or from direct non-perturbative calculations, as well as more insights about the value of $\rho_D^3 (B)$ extracted from fit, would be very helpful in reducing the corresponding theory uncertainties.

We study the structure of atomic hydrogen (H I) in the host galaxy of GRB 171205A/SN 2017iuk at z = 0.037 through H I 21 cm emission line observations with the Karl G. Jansky Very Large Array. These observations reveal unusual morphology and kinematics of the H I in this otherwise apparently normal galaxy. High column density, cold H I is absent from an extended North-South region passing by the optical center of the galaxy, but instead is extended toward the South, on both sides of the galaxy. Moreover, the H I kinematics do not show a continuous change along the major axis of the galaxy as expected in a classical rotating disk. We explore several scenarios to explain the H I structure and kinematics in the galaxy: feedback from a central starburst and/or an active galactic nucleus, ram-pressure stripping, accretion, and tidal interaction from a companion galaxy. All of these options are ruled out. The most viable remaining explanation is the penetrating passage of a satellite through the disk only a few Myr ago, redistributing the H I in the GRB host without yet affecting its stellar distribution. It can also lead to the rapid formation of peculiar stars due to a violent induced shock. The location of GRB 171205A in the vicinity of the distorted area suggests that its progenitor star(s) originated in extreme conditions that share the same origin as the peculiarities in H I. This could explain the atypical location of GRB 171205A in its host galaxy.

We present GIGANTES, the most extensive and realistic void catalog suite ever released-containing over 1 billion cosmic voids covering a volume larger than the observable universe, more than 20 TB of data, and created by running the void finder VIDE on QUIJOTE's halo simulations. The GIGANTES suite, spanning thousands of cosmological models, opens up the study of voids, answering compelling questions: Do voids carry unique cosmological information? How is this information correlated with galaxy information? Leveraging the large number of voids in the GIGANTES suite, our Fisher constraints demonstrate voids contain additional information, critically tightening constraints on cosmological parameters. We use traditional void summary statistics (void size function, void density profile) and the void autocorrelation function, which independently yields an error of 0.13 eV on ∑ m _ν for a 1 h ^-3 Gpc³ simulation, without cosmic microwave background priors. Combining halos and voids we forecast an error of 0.09 eV from the same volume, representing a gain of 60% compared to halos alone. Extrapolating to next generation multi-Gpc³ surveys such as the Dark Energy Spectroscopic Instrument, Euclid, the Spectro-Photometer for the History of the Universe and Ices Explorer, and the Roman Space Telescope, we expect voids should yield an independent determination of neutrino mass. Crucially, GIGANTES is the first void catalog suite expressly built for intensive machine-learning exploration. We illustrate this by training a neural network to perform likelihood-free inference on the void size function, giving a ~20% constraint on Ω_m. Cosmology problems provide an impetus to develop novel deep-learning techniques. With GIGANTES, machine learning gains an impressive data set, offering unique problems that will stimulate new techniques.

Using one of the largest volumes of the hydrodynamical cosmological simulation suit Magneticum, we study the evolution of protoclusters identified at redshift = 4, with properties similar to SPT2349-56. We identify 42 protoclusters in the simulation, as massive and equally rich in substructures as observed, confirming that these structures are already virialized. The dynamics of the internally fast rotating member galaxies within these protoclusters resembles observations, merging rapidly to form the cores of the BCGs of the assembling clusters. Half of the gas reservoir of these structures is in a hot phase, with the metal-enrichment at a very early stage. These systems show a good agreement with the observed amount of cold star-forming gas, largely enriched to solar values. We predict that some of the member galaxies are already quenched at z = 4, rendering them undetectable through measurements of their gas reservoir. Tracing the evolution of protoclusters reveals that none of the typical mass indicators at high redshift are good tracers to predict the present-day mass of the system. We find that none of the simulated protoclusters with properties as SPT2349-56 at z = 4.3, are among the top ten most massive clusters at redshift z = 0, with some barely reaching masses of M = 2 x 10^14Msun. Although the average star-formation and mass-growth rates in the simulated galaxies match observations at high redshift reasonably well, the simulation fails to reproduce the extremely high total star-formation rates within observed protoclusters, indicating that the sub-grid models are lacking the ability to reproduce higher star-formation efficiency (or lower depletion timescales).

We provide a numerical framework with which spatially and spectrally accurate representations of X-ray binary populations can be studied from hydrodynamical cosmological simulations. We construct average spectra accounting for a hot gas component and verify the emergence of observed scaling relations between galaxy wide X-ray luminosity ($L_{X}$) and stellar mass ($M_{\star}$) as well as star-formation rate (SFR). Using simulated galaxy halos extracted from the $(48\,h^{-1} \mathrm{cMpc})^3$ volume of the Magneticum Pathfinder cosmological simulations at $z = 0.07$ we generate mock spectra with the X-ray photon-simulator Phox. We extend the Phox code to account for the stellar component in the simulation and study the resulting contribution in composite galactic spectra. Average X-ray luminosity functions are perfectly reproduced up to the one-photon luminosity limit. Comparing our resulting $L_{X}-\mathrm{SFR}-M_{\star}$ relation for X-ray binaries with recent observations of field galaxies in the Virgo galaxy cluster we find significant overlap. Invoking a metallicity dependent model for high-mass X-ray binaries yields an anti-correlation between mass-weighted stellar metallicity and SFR normalized luminosity. The spatial distribution of high-mass X-ray binaries coincides with star-formation regions of simulated galaxies while low-mass X-ray binaries follow the stellar mass surface density. X-ray binary emission is the dominant contribution in the 2-10 keV band in the absence of an actively accreting central super-massive black hole with 50% contribution in the 0.5-2 keV band rivaling the hot gas component. Our modelling remains consistent with observations despite uncertainties connected to our approach. The predictive power and easily extendable framework hold great value for future investigations of galactic X-ray spectra.

Small grains play an essential role in astrophysical processes such as chemistry, radiative transfer, gas/dust dynamics. The population of small grains is mainly maintained by the fragmentation process due to colliding grains. An accurate treatment of dust fragmentation is required in numerical modelling. However, current algorithms for solving fragmentation equation suffer from an over-diffusion in the conditions of 3D simulations. To tackle this challenge, we developed a Discontinuous Galerkin scheme to solve efficiently the non-linear fragmentation equation with a limited number of dust bins.

In this paper we present COMET, a Gaussian process emulator of the galaxy power spectrum multipoles in redshift-space. The model predictions are based on one-loop perturbation theory and we consider two alternative descriptions of redshift-space distortions: one that performs a full expansion of the real- to redshift-space mapping, as in recent effective field theory models, and another that preserves the non-perturbative impact of small-scale velocities by means of an effective damping function. The outputs of COMET can be obtained at arbitrary redshifts (up to $z \sim 3$), for arbitrary fiducial background cosmologies, and for a large parameter space that covers the shape parameters $\omega_c$, $\omega_b$, and $n_s$, as well as the evolution parameters $h$, $A_s$, $\Omega_K$, $w_0$, and $w_a$. This flexibility does not impair COMET's accuracy, since we exploit an exact degeneracy between the evolution parameters that allows us to train the emulator on a significantly reduced parameter space. While the predictions are sped up by at least two orders of magnitude, validation tests reveal an accuracy of $0.1\,\%$ for the monopole and quadrupole ($0.3\,\%$ for the hexadecapole), or alternatively, better than $0.25\,\sigma$ for all three multipoles in comparison to statistical uncertainties expected for the Euclid survey with a tenfold increase in volume. We show that these differences translate into shifts in mean posterior values that are at most of the same size, meaning that COMET can be used with the same confidence as the exact underlying models. COMET is a publicly available Python package that also provides the tree-level bispectrum multipoles in redshift-space and Gaussian covariance matrices.

This work demonstrates that nontopological solitons with large global charges and masses, even above the Planck scale, can form in the early universe and dominate the dark matter abundance. In solitosynthesis, solitons prefer to grow as large as possible under equilibrium dynamics when an initial global charge asymmetry is present. Their abundance is set by when soliton formation via particle fusion freezes out, and their charges are set by the time it takes to accumulate free particles. This work improves the estimation of both quantities, and in particular shows that much larger-charged solitons form than previously thought. The results are estimated analytically and validated numerically by solving the coupled Boltzmann equations. Without solitosynthesis, phase transitions can still form solitons from particles left inside false-vacuum pockets and determine their present-day abundance and properties. Even with zero charge asymmetry, solitons formed in this way can have very large charges on account of statistical fluctuations in the numbers of (anti)particles inside each pocket.

Context. Weak lensing and clustering statistics beyond two-point functions can capture non-Gaussian information about the matter density field, thereby improving the constraints on cosmological parameters relative to the mainstream methods based on correlation functions and power spectra. Aims. This paper presents a cosmological analysis of the fourth data release of the Kilo Degree Survey based on the density split statistics, which measures the mean shear profiles around regions classified according to foreground densities. The latter is constructed from a bright galaxy sample, which we further split into red and blue samples, allowing us to probe their respective connection to the underlying dark matter density. Methods. We use the state-of-the-art model of the density splitting statistics and validate its robustness against mock data infused with known systematic effects such as intrinsic galaxy alignment and baryonic feedback. Results. After marginalising over the photometric redshift uncertainty and the residual shear calibration bias, we measure for the full KiDS-bright sample a structure growth parameter of $S_8 = \sigma_8 \sqrt{\Omega_\mathrm{m}/0.3} = 0.74^{+0.03}_{-0.02}$ that is competitive to and consistent with two-point cosmic shear results, a matter density of $\Omega_\mathrm{m} = 0.28 \pm 0.02$, and a constant galaxy bias of $b = 1.32^{+0.12}_{-0.10}$.

A thermonuclear explosion triggered by a He-shell detonation on a carbon-oxygen white-dwarf core has been predicted to have strong UV line blanketing at early times due to the iron-group elements produced during He-shell burning. We present the photometric and spectroscopic observations of SN 2016dsg, a subluminous peculiar Type I supernova consistent with a thermonuclear explosion involving a thick He shell. With a redshift of 0.04, the i-band peak absolute magnitude is derived to be around -17.5. The object is located far away from its host, an early-type galaxy, suggesting it originated from an old stellar population. The spectra collected after the peak are unusually red, show strong UV line blanketing and weak O I λ7773 absorption lines, and do not evolve significantly over 30 days. An absorption line around 9700-10500 Å is detected in the near-infrared spectrum and is likely from the unburnt He in the ejecta. The spectroscopic evolution is consistent with the thermonuclear explosion models for a sub-Chandrasekhar-mass white dwarf with a thick He shell, while the photometric evolution is not well described by existing models.

We investigate shock structures driven by merger events in high-resolution simulations that result in a galaxy with a virial mass M ~ 1e12 Msol. We find that the sizes and morphologies of the internal shocks resemble remarkably well those of the newly-detected class of odd radio circles (ORCs). This would highlight a so-far overlooked mechanism to form radio rings, shells and even more complex structures around elliptical galaxies. Mach numbers of M = 2-3 for such internal shocks are in agreement with the spectral indices of the observed ORCs. We estimate that ~5 percent of galaxies could undergo merger events which occasionally lead to such prominent structures within the galactic halo during their lifetime, explaining the low number of observed ORCs. At the time when the shock structures are matching the physical sizes of the observed ORCs, the central galaxies are typically classified as early-type galaxies, with no ongoing star formation, in agreement with observational findings. Although the energy released by such mergers could potentially power the observed radio luminosity already in Milky-Way-like halos, our predicted luminosity from a simple, direct shock acceleration model is much smaller than the observed one. Considering the estimated number of candidates from our cosmological simulations and the higher observed energies, we suggest that the proposed scenario is more likely for halo masses around 1e13 Msol in agreement with the observed stellar masses of the galaxies at the center of ORCs.

We present multiple results on the production of loosely-bound molecules in bottomonium annihilations and $e^+e^-$ collisions at $\sqrt{s} = 10.58$ GeV. We perform the first comprehensive test of several models for deuteron production against all the existing data in this energy region. We fit the free parameters of the models to reproduce the observed cross sections, and we predict the deuteron spectrum and production and the cross section for the $e^+e^- \to d\bar{d} + X$ process both at the $\Upsilon(1,2,3S)$ resonances and at $\sqrt{s}=10.58$ GeV. The predicted spectra show differences but are all compatible with the uncertainties of the existing data. These differences could be addressed if larger datasets are collected by the Belle~II experiment. Fixing the source size parameter to reproduce the deuteron data, we then predict the production rates for $H$ dibaryon and hypertriton in this energy region using a simple coalescence model. Our prediction on $H$ dibaryon production rate is below the limits set by the direct search at the Belle experiment, but in the range accessible to the Belle~II experiment. The systematic effect due to the MC modelling of quarks and gluon fragmentation into baryons is reduced deriving a new tuning of the \pythia MonteCarlo generator using the available measurement of single- and double-particle spectra in $\Upsilon$ decays.

Euclid's photometric galaxy cluster survey has the potential to be a very competitive cosmological probe. The main cosmological probe with observations of clusters is their number count, within which the halo mass function (HMF) is a key theoretical quantity. We present a new calibration of the analytic HMF, at the level of accuracy and precision required for the uncertainty in this quantity to be subdominant with respect to other sources of uncertainty in recovering cosmological parameters from Euclid cluster counts. Our model is calibrated against a suite of N-body simulations using a Bayesian approach taking into account systematic errors arising from numerical effects in the simulation. First, we test the convergence of HMF predictions from different N-body codes, by using initial conditions generated with different orders of Lagrangian Perturbation theory, and adopting different simulation box sizes and mass resolution. Then, we quantify the effect of using different halo-finder algorithms, and how the resulting differences propagate to the cosmological constraints. In order to trace the violation of universality in the HMF, we also analyse simulations based on initial conditions characterised by scale-free power spectra with different spectral indexes, assuming both Einstein--de Sitter and standard $\Lambda$CDM expansion histories. Based on these results, we construct a fitting function for the HMF that we demonstrate to be sub-percent accurate in reproducing results from 9 different variants of the $\Lambda$CDM model including massive neutrinos cosmologies. The calibration systematic uncertainty is largely sub-dominant with respect to the expected precision of future mass-observation relations; with the only notable exception of the effect due to the halo finder, that could lead to biased cosmological inference.

Dark matter (DM) with self-interactions is a promising solution for the small-scale problems of the standard cosmological model. Here we perform the first cosmological simulation of frequent DM self-interactions, corresponding to small-angle DM scatterings. The focus of our analysis lies in finding and understanding differences to the traditionally assumed rare DM (large-angle) self scatterings. For this purpose, we compute the distribution of DM densities, the matter power spectrum, the two-point correlation function and the halo and subhalo mass functions. Furthermore, we investigate the density profiles of the DM haloes and their shapes. We find that overall large-angle and small-angle scatterings behave fairly similarly with a few exceptions. In particular, the number of satellites is considerably suppressed for frequent compared to rare self-interactions with the same cross-section. Overall we observe that while differences between the two cases may be difficult to establish using a single measure, the degeneracy may be broken through a combination of multiple ones. For instance, the combination of satellite counts with halo density or shape profiles could allow discriminating between rare and frequent self-interactions. As a by-product of our analysis, we provide - for the first time - upper limits on the cross-section for frequent self-interactions.

Many processes during the evolution of protoplanetary disks and during planet formation are highly sensitive to the sizes of dust particles that are present in the disk: the efficiency of dust accretion in the disk and volatile transport on dust particles, gravoturbulent instabilities leading to the formation of planetesimals, or the accretion of pebbles onto large planetary embryos to form giant planets are typical examples of processes that depend on the sizes of the dust particles involved. Furthermore, radiative properties like absorption or scattering opacities depend on the particle sizes. To interpret observations of dust in protoplanetary disks, a proper estimate of the dust particle sizes is needed. We present DustPy: a Python package to simulate dust evolution in protoplanetary disks. DustPy solves gas and dust transport including viscous advection and diffusion as well as collisional growth of dust particles. DustPy is written with a modular concept, such that every aspect of the model can be easily modified or extended to allow for a multitude of research opportunities.

We report the discovery and characterization of two small transiting planets orbiting the bright M3.0V star TOI-1468 (LSPM J0106+1913), whose transit signals were detected in the photometric time series in three sectors of the TESS mission. We confirm the e planetary nature of both of them using precise radial velocity measurements from the CARMENES and MAROON-X spectrographs, and supplement them with ground-based transit photometry. A joint analysis of all these data reveals that the shorter-period planet, TOI-1468 b ($P_{\rm b}$ = 1.88 d), has a planetary mass of $M_{\rm b} = 3.21\pm0.24$ $M_{\oplus}$ and a radius of $R_{\rm b} =1.280^{+0.038}_{-0.039} R_{\oplus}$, resulting in a density of $\rho_{\rm b} = 8.39^{+ 1.05}_{- 0.92}$ g cm$^{-3}$, which is consistent with a mostly rocky composition. For the outer planet, TOI-1468 c ($P_{\rm c} = 15.53$ d), we derive a mass of $M_{\rm c} = 6.64^{+ 0.67}_{- 0.68}$ $M_{\oplus}$, a radius of $R_{\rm c} = 2.06\pm0.04\,R_{\oplus}$, and a bulk density of $\rho_{c} = 2.00^{+ 0.21}_{- 0.19}$ g cm$^{-3}$, which corresponds to a rocky core composition with a H/He gas envelope. These planets are located on opposite sides of the radius valley, making our system an interesting discovery as there are only a handful of other systems with the same properties. This discovery can further help determine a more precise location of the radius valley for small planets around M dwarfs and, therefore, shed more light on planet formation and evolution scenarios.

Context. The rapidly improving quality and resolution of both low surface brightness observations and cosmological simulations of galaxies enables an increasingly thorough investigation of the imprints of the formation history in the outer, unrelaxed regions of galaxies, and a direct comparison to another tracer of galaxy formation, the internal kinematics. Aims. Using the state-of-the-art hydrodynamical cosmological simulation Magneticum Pathfinder, we identify tidal tails, shells, streams, and satellite planes, and connect their existence to the amount of rotational support and the formation histories of the host galaxies. Methods. Tidal features are visually classified from a three-dimensional rendering of the simulated galaxies by several scientists independently. Only features that were identified by at least half of the participating individuals are considered as existing features. The results are compared to observations of the MATLAS survey. Results. Shells are preferentially found around kinematically slowly rotating galaxies in both simulations and observations, while streams can be found around all kind of galaxies with a slightly higher probability to be present around less rotationally supported galaxies. Tails and satellite planes, however, appear independently of the internal kinematics of the central galaxy, indicating that they are formed through processes that have not (yet) affected the internal kinematics. Conclusions. As shells are formed through radial merger events while streams are remnants of circular merger infall, this suggests that the orbital angular momentum of the merger event could play a more crucial role in transforming the host galaxy than previously anticipated. The existence of a shell around a given slow rotator can further be used to distinguish the radial merger formation scenario from other formation pathways of slow rotators.

After playing a pivotal role in the birth of the Standard Model in the 70s, the study of charm physics has undergone a revival during the last decade, triggered by a wealth of precision measurements from the charm and B factories, and from the CDF and especially the LHCb experiments. In this paper, we sum up how the unique phenomenology of charmed hadrons can be used to test the Standard Model and we review the latest measurements performed in this field by the LHCb experiment. These include the historic first observations of CP violation and of a nonzero mass difference between the charmed neutral-meson mass eigenstates, the most precise determination of their decay-width difference to date, and a search for time-dependent CP violation reaching the unprecedented precision of 10−4. These results challenge our comprehension of nonperturbative strong interactions, and their interpretation calls for further studies on both the theoretical and experimental sides. The upcoming upgrades of the LHCb experiment will play a leading role in this quest.

Negative feedback from active galactic nuclei (AGN) is the leading mechanism for the quenching of massive galaxies in the vast majority of modern galaxy evolution models. However, direct observational evidence that AGN feedback causes quenching on a population scale is lacking. Studies have shown that luminous AGN are preferentially located in gas-rich and star-forming galaxies, an observation that has sometimes been suggested to be in tension with a negative AGN feedback picture. We investigate three of the current cosmological simulations (IllustrisTNG, EAGLE, and SIMBA) along with post-processed models for molecular hydrogen gas masses and perform similar tests to those used by observers. We find that the simulations predict: (i) no strong negative trends between L_bol and $f_{\mathrm{ H}_2}$ or specific star formation rate (sSFR); (ii) both high-luminosity ($L_{\rm {bol}} \ge 10^{44}\rm {\, erg\, s^{-1}}$) and high Eddington ratio (λ_Edd $\ge 1{{\ \rm per\ cent}}$) AGN are preferentially located in galaxies with high molecular gas fractions and sSFR; and (iii) that the gas-depleted and quenched fractions of AGN host galaxies are lower than a control sample of non-active galaxies. These three findings are in qualitative agreement with observational samples at z = 0 and z = 2 and show that such results are not in tension with the presence of strong AGN feedback, which all simulations we employ require to produce realistic massive galaxies. However, we also find quantifiable differences between predictions from the simulations, which could allow us to observationally test the different subgrid feedback models.

The 1.33-mm survey of protoplanetary discs in the Taurus molecular cloud found annular gaps and rings to be common in extended sources (≳ 55AU), when their 1D visibility distributions were fit parametrically. We first demonstrate the advantages and limitations of non-parametric visibility fits for data at the survey's 0.12-arcsec resolution. Then we use the non-parametric model in Frankenstein (frank) to identify new substructure in three compact and seven extended sources. Among the new features, we identify three trends: a higher occurrence rate of substructure in the survey's compact discs than previously seen, underresolved (potentially azimuthally asymmetric) substructure in the innermost disc of extended sources, and a 'shoulder' on the trailing edge of a ring in discs with strong depletion at small radii. Noting the shoulder morphology is present in multiple discs observed at higher resolution, we postulate it is tracing a common physical mechanism. We further demonstrate how a superresolution frank brightness profile is useful in motivating an accurate parametric model, using the highly structured source DL Tau in which frank finds two new rings. Finally, we show that sparse (u, v) plane sampling may be masking the presence of substructure in several additional compact survey sources.

Using a new sample of extremely metal poor systems, the EMPRESS survey has recently reported a primordial helium abundance that is $3\sigma$ smaller than the prediction from the Standard BBN scenario. This measurement could be interpreted as a hint for a primordial lepton asymmetry in the electron neutrino flavor. Motivated by the EMPRESS results, we present a comprehensive analysis of the lepton asymmetry using measurements of the abundances of primordial elements, along with CMB data from Planck. Assuming that there is no dark radiation in our Universe, we find an electron neutrino chemical potential $\xi_{\nu_e} = 0.037 \pm 0.013$, which deviates from zero by $2.8\sigma$. If no assumption is made on the abundance of dark radiation in the Universe, the chemical potential is $\xi_{\nu_e} = 0.037 \pm 0.020$, which deviates from zero by $1.9\sigma$. We also find that this result is rather insensitive to the choice of nuclear reaction rates. If the true helium abundance corresponds to the EMPRESS central value, future CMB observations from the Simons Observatory and CMB-S4 will increase the significance for a non-zero lepton asymmetry to $4\sigma$ and $5\sigma$ respectively, assuming no dark radiation, or to $3\sigma$ when no assumption is made on the abundance of dark radiation.

We investigate shock structures driven by merger events in high-resolution simulations that result in a galaxy with a virial mass M ~ 1e12 Msol. We find that the sizes and morphologies of the internal shocks resemble remarkably well those of the newly-detected class of odd radio circles (ORCs). This would highlight a so-far overlooked mechanism to form radio rings, shells and even more complex structures around elliptical galaxies. Mach numbers of M = 2-3 for such internal shocks are in agreement with the spectral indices of the observed ORCs. We estimate that ~5 percent of galaxies could undergo merger events which occasionally lead to such prominent structures within the galactic halo during their lifetime, explaining the low number of observed ORCs. At the time when the shock structures are matching the physical sizes of the observed ORCs, the central galaxies are typically classified as early-type galaxies, with no ongoing star formation, in agreement with observational findings. Although the energy released by such mergers could potentially power the observed radio luminosity already in Milky-Way-like halos, our predicted luminosity from a simple, direct shock acceleration model is much smaller than the observed one. Considering the estimated number of candidates from our cosmological simulations and the higher observed energies, we suggest that the proposed scenario is more likely for halo masses around 1e13 Msol in agreement with the observed stellar masses of the galaxies at the center of ORCs.

We present a Hubble Space Telescope (HST) weak gravitational lensing study of nine distant and massive galaxy clusters with redshifts $1.0 \lesssim z \lesssim 1.7$ ($z_\mathrm{median} = 1.4$) and Sunyaev Zel'dovich (SZ) detection significance $\xi > 6.0$ from the South Pole Telescope Sunyaev Zel'dovich (SPT-SZ) survey. We measured weak lensing galaxy shapes in HST/ACS F606W and F814W images and used additional observations from HST/WFC3 in F110W and VLT/FORS2 in $U_\mathrm{HIGH}$ to preferentially select background galaxies at $z\gtrsim 1.8$, achieving a high purity. We combined recent redshift estimates from the CANDELS/3D-HST and HUDF fields to infer an improved estimate of the source redshift distribution. We measured weak lensing masses by fitting the tangential reduced shear profiles with spherical Navarro-Frenk-White (NFW) models. We obtained the largest lensing mass in our sample for the cluster SPT-CLJ2040$-$4451, thereby confirming earlier results that suggest a high lensing mass of this cluster compared to X-ray and SZ mass measurements. Combining our weak lensing mass constraints with results obtained by previous studies for lower redshift clusters, we extended the calibration of the scaling relation between the unbiased SZ detection significance $\zeta$ and the cluster mass for the SPT-SZ survey out to higher redshifts. We found that the mass scale inferred from our highest redshift bin ($1.2 < z < 1.7$) is consistent with an extrapolation of constraints derived from lower redshifts, albeit with large statistical uncertainties. Thus, our results show a similar tendency as found in previous studies, where the cluster mass scale derived from the weak lensing data is lower than the mass scale expected in a Planck $\nu\Lambda$CDM (i.e. $\nu$$\Lambda$ Cold Dark Matter) cosmology given the SPT-SZ cluster number counts.

A key component of the Dark Energy Spectroscopic Instrument (DESI) survey validation (SV) is a detailed visual inspection (VI) of the optical spectroscopic data to quantify key survey metrics. In this paper we present results from VI of the quasar survey using deep coadded SV spectra. We show that the majority (~70%) of the main-survey targets are spectroscopically confirmed as quasars, with ~16% galaxies, ~6% stars, and ~8% low-quality spectra lacking reliable features. A non-negligible fraction of the quasars are misidentified by the standard DESI spectroscopic pipeline but we show that the majority can be recovered using post-pipeline "afterburner" quasar-identification approaches. We combine these "afterburners" with our standard pipeline to create a modified pipeline to improve the overall quasar completeness. At the depth of the main DESI survey both pipelines achieve a good-redshift purity (reliable redshifts measured within 3000 km/s) of ~99%; however, the modified pipeline recovers ~94% of the visually inspected quasars, as compared to just ~86% from the standard pipeline. We demonstrate that both pipelines achieve an overall redshift precision and accuracy of ~100 km/s and ~70 km/s, respectively. We constructed composite spectra to investigate why some quasars are missed by the standard spectroscopic pipeline and find that they are more host-galaxy dominated and/or dust reddened than the standard-pipeline quasars. We also show example spectra to demonstrate the overall diversity of the DESI quasar sample and provide strong-lensing candidates where two targets contribute to a single DESI spectrum.

Cosmological shock waves are essential to understanding the formation of cosmological structures. To study them, scientists run computationally expensive high-resolution 3D hydrodynamic simulations. Interpreting the simulation results is challenging because the resulting data sets are enormous, and the shock wave surfaces are hard to separate and classify due to their complex morphologies and multiple shock fronts intersecting. We introduce a novel pipeline, Virgo, combining physical motivation, scalability, and probabilistic robustness to tackle this unsolved unsupervised classification problem. To this end, we employ kernel principal component analysis with low-rank matrix approximations to denoise data sets of shocked particles and create labeled subsets. We perform supervised classification to recover full data resolution with stochastic variational deep kernel learning. We evaluate on three state-of-the-art data sets with varying complexity and achieve good results. The proposed pipeline runs automatically, has only a few hyperparameters, and performs well on all tested data sets. Our results are promising for large-scale applications, and we highlight now enabled future scientific work.

Using the Schwinger-Keldysh-formalism, reformulated in [1] as an effective field theory in Euclidean anti-de Sitter, we evaluate the one-loop cosmological four-point function of a conformally coupled interacting scalar field in de Sitter. Recasting the Witten cosmological correlator as flat space Feynman integrals, we evaluate the one-loop cosmological four-point functions in de Sitter space in terms of single-valued multiple polylogarithms. From it we derive anomalous dimensions and OPE coefficients of the dual conformal field theory at space-like, future infinity. In particular, we find an interesting degeneracy in the anomalous dimensions relating operators of neighboring spins.

We report the discovery of a multiplanetary system transiting the M0 V dwarf HD 260655 (GJ 239, TOI-4599). The system consists of at least two transiting planets, namely HD 260655 b, with a period of 2.77 d, a radius of R_b = 1.240 ± 0.023 R_⊕, a mass of M_b = 2.14 ± 0.34 M_⊕, and a bulk density of ρ_b = 6.2 ± 1.0 g cm⁻³, and HD 260655 c, with a period of 5.71 d, a radius of {R_c} = 1.533 _{- 0.046} ^{+ 0.051}{R_ \oplus }, a mass of M_c = 3.09 ± 0.48 M_⊕, and a bulk density of {ρ _c} = 4.7 _{- 0.8} ^{+ 0.9}{{g}} g cm⁻³. The planets have been detected in transit by the Transiting Exoplanet Survey Satellite (TESS) mission and confirmed independently with archival and new precise radial velocities obtained with the HIRES and CARMENES instruments since 1998 and 2016, respectively. At a distance of 10 pc, HD 260655 has become the fourth closest known multitransiting planet system after HD 219134, LTT 1445 A, and AU Mic. Due to the apparent brightness of the host star (J = 6.7 mag), both planets are among the most suitable rocky worlds known today for atmospheric studies with the James Webb Space Telescope, both in transmission and emission.

The long-standing controversy about the isospin dependence of the effective Dirac mass in ab initio calculations of asymmetric nuclear matter is clarified by solving the relativistic Brueckner-Hartree-Fock equations in the full Dirac space. The symmetry energy and its slope parameter at the saturation density are E_sym(ρ₀) =33.1 MeV and L =65.2 MeV, in agreement with empirical and experimental values. Further applications predict the neutron star radius R_{1.4 M ⊙}≈12 km and the maximum mass of a neutron star M_max≤2.4 M_⊙ .

Biological processes operate in a spatially and temporally ordered manner to reliably fulfill their function. This is achieved by pattern formation, which generally involves many different spatial and temporal scales. The resulting multiscale patterns exhibit complex dynamics for which it is difficult to find a simplified description at large scales while preserving information about the patterns at small scales. Here, we introduce an approach for mass-conserving reaction-diffusion systems that is based on a linear theory and therefore conceptually simple to apply. We investigate multiscale patterns of the Min protein system and show that our approach enables us to explain and predict the intricate dynamics from the large-scale mass redistribution of the total protein densities.

Using the potential non-relativistic quantum chromodynamics (pNRQCD) effective field theory, we derive a Lindblad equation for the evolution of the heavy-quarkonium reduced density matrix that is accurate to next-to-leading order (NLO) in the ratio of the binding energy of the state to the temperature of the medium. The resulting NLO Lindblad equation can be used to more reliably describe heavy-quarkonium evolution in the quark-gluon plasma at low temperatures compared to the leading-order truncation. For phenomenological application, we numerically solve the resulting NLO Lindblad equation using the quantum trajectories algorithm. To achieve this, we map the solution of the three-dimensional Lindblad equation to the solution of an ensemble of one-dimensional Schrödinger evolutions with Monte-Carlo sampled quantum jumps. Averaging over the Monte-Carlo sampled quantum jumps, we obtain the solution to the NLO Lindblad equation without truncation in the angular momentum quantum number of the states considered. We also consider the evolution of the system using only the complex effective Hamiltonian without stochastic jumps and find that this provides a reliable approximation for the ground state survival probability at LO and NLO. Finally, we make comparisons with our prior leading-order pNRQCD results and experimental data available from the ATLAS, ALICE, and CMS collaborations.

Like most spiral galaxies, the Milky Way contains a population of blue, metal-poor globular clusters and another of red, metal-rich ones. Most of the latter belong to the bulge, and therefore they are poorly studied compared to the blue (halo) ones because they suffer higher extinction and larger contamination from field stars. These intrinsic difficulties, together with a lack of low-mass bulge globular clusters, are reasons to believe that their census is not complete yet. Indeed, a few new clusters have been confirmed in the last few years. One of them is VVV CL001, the subject of the present study. We present a new spectroscopic analysis of the recently confirmed globular cluster VVV CL001, made by means of MUSE@VLT integral field data. Individual spectra were extracted for stars in the VVV CL001 field. Radial velocities were derived by cross-correlation with synthetic templates. Coupled with proper motions from the VVV (VISTA Variables in the Vía Láctea) survey, these data allow us to select 55 potential cluster members, for which we derive metallicities using the public code THE CANNON. The mean radial velocity of the cluster is V_helio = -324.9 ± 0.8 km s^-1, as estimated from 55 cluster members. This high velocity, together with a low metallicity [Fe/H] = -2.04 ± 0.02 dex, suggests that VVV CL001 could be a very old cluster. The estimated distance is d = 8.23 ± 0.46 kpc, placing the cluster in the Galactic bulge. Furthermore, both its current position and the orbital parameters suggest that VVV CL001 is most probably a bulge globular cluster.

In star-forming galaxies, the far-infrared (FIR) and radio-continuum luminosities obey a tight empirical relation over a large range of star-formation rates (SFR). To understand the physics, we examine magneto-hydrodynamic galaxy simulations, which follow the genesis of cosmic ray (CR) protons at supernovae and their advective and anisotropic diffusive transport. We show that gravitational collapse of the proto-galaxy generates a corrugated accretion shock, which injects turbulence and drives a small-scale magnetic dynamo. As the shock propagates outwards and the associated turbulence decays, the large velocity shear between the supersonically rotating cool disc with respect to the (partially) pressure-supported hot circumgalactic medium excites Kelvin-Helmholtz surface and body modes. Those interact non-linearly, inject additional turbulence and continuously drive multiple small-scale dynamos, which exponentially amplify weak seed magnetic fields. After saturation at small scales, they grow in scale to reach equipartition with thermal and CR energies in Milky Way-mass galaxies. In small galaxies, the magnetic energy saturates at the turbulent energy while it fails to reach equipartition with thermal and CR energies. We solve for steady-state spectra of CR protons, secondary electrons/positrons from hadronic CR-proton interactions with the interstellar medium, and primary shock-accelerated electrons at supernovae. The radio-synchrotron emission is dominated by primary electrons, irradiates the magnetised disc and bulge of our simulated Milky Way-mass galaxy and weakly traces bubble-shaped magnetically-loaded outflows. Our star-forming and star-bursting galaxies with saturated magnetic fields match the global FIR-radio correlation (FRC) across four orders of magnitude. Its intrinsic scatter arises due to (i) different magnetic saturation levels that result from different seed magnetic fields, (ii) different radio synchrotron luminosities for different specific SFRs at fixed SFR and (iii) a varying radio intensity with galactic inclination. In agreement with observations, several 100-pc-sized regions within star-forming galaxies also obey the FRC, while the centres of starbursts substantially exceed the FRC.