Spatial proton gradients create energy in biological systems and are likely a driving force for prebiotic systems. Due to the fast diffusion of protons, they are however difficult to create as steady state, unless driven by other non-equilibria such as thermal gradients. Here, we quantitatively predict the heat-flux driven formation of pH gradients for the case of a simple acid-base reaction system. To this end, we (i) establish a theoretical framework that describes the spatial interplay of chemical reactions with thermal convection, thermophoresis, and electrostatic forces by a separation of timescales, and (ii) report quantitative measurements in a purpose-built microfluidic device. We show experimentally that the slope of such pH gradients undergoes pronounced amplitude changes in a concentration-dependent manner and can even be inverted. The predictions of the theoretical framework fully reflect these features and establish an understanding of how naturally occurring non-equilibrium environmental conditions can drive pH gradients.

We present N-body simulations, including post-Newtonian dynamics, of dense clusters of low-mass stars harbouring central black holes (BHs) with initial masses of 50, 300, and 2000 M_⊙. The models are evolved with the N-body code BIFROST to investigate the possible formation and growth of massive BHs by the tidal capture of stars and tidal disruption events (TDEs). We model star-BH tidal interactions using a velocity-dependent drag force, which causes orbital energy and angular momentum loss near the BH. About ~20-30 per cent of the stars within the spheres of influence of the black holes form Bahcall-Wolf cusps and prevent the systems from core collapse. Within the first 40 Myr of evolution, the systems experience 500-1300 TDEs, depending on the initial cluster structure. Most (>95 per cent) of the TDEs originate from stars in the Bahcall-Wolf cusp. We derive an analytical formula for the TDE rate as a function of the central BH mass, density, and velocity dispersion of the clusters ($\dot{N}_{\mathrm{TDE}} \propto M\mathrm{_{BH}}\rho \sigma ^{-3}$). We find that TDEs can lead a 300 M_⊙ BH to reach $\sim 7000 \, \mathrm{{M}_{\odot }}$ within a Gyr. This indicates that TDEs can drive the formation and growth of massive BHs in sufficiently dense environments, which might be present in the central regions of nuclear star clusters.

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.

Recent observations have shown that the atmospheres of ultrahot Jupiters (UHJs) commonly possess temperature inversions, where the temperature increases with increasing altitude. Nonetheless, which opacity sources are responsible for the presence of these inversions remains largely observationally unconstrained. We used LBT/PEPSI to observe the atmosphere of the UHJ KELT-20 b in both transmission and emission in order to search for molecular agents which could be responsible for the temperature inversion. We validate our methodology by confirming a previous detection of Fe I in emission at 16.9σ. Our search for the inversion agents TiO, VO, FeH, and CaH results in non-detections. Using injection-recovery testing we set 4σ upper limits upon the volume mixing ratios for these constituents as low as ~1 × 10^-9 for TiO. For TiO, VO, and CaH, our limits are much lower than expectations from an equilibrium chemical model, while we cannot set constraining limits on FeH with our data. We thus rule out TiO and CaH as the source of the temperature inversion in KELT-20 b, and VO only if the line lists are sufficiently accurate. *Based on data acquired with the Potsdam Echelle Polarimetric and Spectroscopic Instrument (PEPSI) using the Large Binocular Telescope (LBT) in Arizona.

Earth and other rocky objects in the inner Solar system are depleted in carbon compared to objects in the outer Solar system, the Sun, or the ISM. It is believed that this is a result of the selective removal of refractory carbon from primordial circumstellar material. In this work, we study the irreversible release of carbon into the gaseous environment via photolysis and pyrolysis of refractory carbonaceous material during the disc phase of the early Solar system. We analytically solve the one-dimensional advection equation and derive an explicit expression that describes the depletion of carbonaceous material in solids under the influence of radial and vertical transport. We find both depletion mechanisms individually fail to reproduce Solar system abundances under typical conditions. While radial transport only marginally restricts photodecomposition, it is the inefficient vertical transport that limits carbon depletion under these conditions. We show explicitly that an increase in the vertical mixing efficiency, and/or an increase in the directly irradiated disc volume, favours carbon depletion. Thermal decomposition requires a hot inner disc (>500 K) beyond 3 au to deplete the formation region of Earth and chondrites. We find FU Ori-type outbursts to produce these conditions such that moderately refractory compounds are depleted. However, such outbursts likely do not deplete the most refractory carbonaceous compounds beyond the innermost disc region. Hence, the refractory carbon abundance at 1 au typically does not reach terrestrial levels. Nevertheless, under specific conditions, we find photolysis and pyrolysis combined to reproduce Solar system abundances.

The integrated shear 3-point correlation function $\zeta_{\pm}$ measures the correlation between the local shear 2-point function $\xi_{\pm}$ and the 1-point shear aperture mass in patches of the sky. Unlike other higher-order statistics, $\zeta_{\pm}$ can be efficiently measured from cosmic shear data, and it admits accurate theory predictions on a wide range of scales as a function of cosmological and baryonic feedback parameters. Here, we develop and test a likelihood analysis pipeline for cosmological constraints using $\zeta_{\pm}$. We incorporate treatment of systematic effects from photometric redshift uncertainties, shear calibration bias and galaxy intrinsic alignments. We also develop an accurate neural-network emulator for fast theory predictions in MCMC parameter inference analyses. We test our pipeline using realistic cosmic shear maps based on $N$-body simulations with a DES Y3-like footprint, mask and source tomographic bins, finding unbiased parameter constraints. Relative to $\xi_{\pm}$-only, adding $\zeta_{\pm}$ can lead to $\approx 10-25\%$ improvements on the constraints of parameters like $A_s$ (or $\sigma_8$) and $w_0$. We find no evidence in $\xi_{\pm} + \zeta_{\pm}$ constraints of a significant mitigation of the impact of systematics. We also investigate the impact of the size of the apertures where $\zeta_{\pm}$ is measured, and of the strategy to estimate the covariance matrix ($N$-body vs. lognormal). Our analysis solidifies the strong potential of the $\zeta_{\pm}$ statistic and puts forward a pipeline that can be readily used to improve cosmological constraints using real cosmic shear data.

We study gravitational back-reaction within relational time formulations of quantum mechanics by considering two versions of time: a time coordinate, modelled as a global quantum degree of freedom, and the proper time of a given physical system, modelled via an internal degree of freedom serving as a local quantum "clock". We show that interactions between coordinate time and mass-energy in a global Wheeler-DeWitt-like constraint lead to gravitational time dilation. In the presence of a massive object this agrees with time dilation in a Schwarzchild metric at leading order in $G$. Furthermore, if two particles couple independently to the time coordinate we show that Newtonian gravitational interaction between those particles emerges in the low energy limit. We also observe features of renormalization of high energy divergences.

A configurable calorimeter simulation for AI (COCOA) applications is presented, based on the Geant4 toolkit and interfaced with the Pythia event generator. This open-source project is aimed to support the development of machine learning algorithms in high energy physics that rely on realistic particle shower descriptions, such as reconstruction, fast simulation, and low-level analysis. Specifications such as the granularity and material of its nearly hermetic geometry are user-configurable. The tool is supplemented with simple event processing including topological clustering, jet algorithms, and a nearest-neighbors graph construction. Formatting is also provided to visualise events using the Phoenix event display software.

The singlet sector of the $O(N),$ $\phi^4$-model in AdS$_4$ at large-$N$, gives rise to a (non-local) dual conformal field theory on the conformal boundary of AdS$_4$, which is a deformation of the generalized free field. We identify and compute a AdS$_4$ 3-point 1-loop fish diagram that controls the exact large-$N$ dimensions and operator product coefficients (OPE) for all "double trace" operators as a function of the renormalized $\phi^4$-coupling. We find that the space of $\phi^4$-coupling is compact with a boundary at the bulk Landau pole where the lowest OPE coefficient diverges.

We revisit stellar energy-loss bounds on the Yukawa couplings $g_{\rm B,L}$ of baryophilic and leptophilic scalars $\phi$. The white-dwarf luminosity function yields $g_{\rm B}\lesssim 7 \times 10^{-13}$ and $g_{\rm L}\lesssim 4 \times 10^{-16}$, based on bremsstrahlung from ${}^{12}{\rm C}$ and ${}^{16}{\rm O}$ collisions with electrons. In models with a Higgs portal, this also implies a bound on the scalar-Higgs mixing angle $\sin \theta \lesssim 2 \times 10^{-10}$. Our new bounds apply for $m_\phi\lesssim {\rm 1~keV}$ and are among the most restrictive ones, whereas for $m_\phi\lesssim 0.5\,{\rm eV}$ long-range force measurements dominate. Besides a detailed calculation of the bremsstrahlung rate for degenerate and semi-relativistic electrons, we prove with a simple argument that non-relativistic bremsstrahlung by the heavy partner is suppressed relative to that by the light one by their squared-mass ratio. This large reduction was overlooked in previous much stronger bounds on $g_{\rm B}$. In an Appendix, we provide fitting formulas (few percent precision) for the bremsstrahlung emission of baryophilic and leptophilic scalars as well as axions for white-dwarf conditions, i.e., degenerate, semi-relativistic electrons and ion-ion correlations in the ``liquid'' phase.

Absorption features in stellar atmospheres are often used to calibrate photocentric velocities for kinematic analysis of further spectral lines. The Li feature at $\sim$ 6708 Å is commonly used, especially in the case of young stellar objects for which it is one of the strongest absorption lines. However, this is a complex line comprising two isotope fine-structure doublets. We empirically measure the wavelength of this Li feature in a sample of young stars from the PENELLOPE/VLT programme (using X-Shooter, UVES and ESPRESSO data) as well as HARPS data. For 51 targets, we fit 314 individual spectra using the STAR-MELT package, resulting in 241 accurately fitted Li features, given the automated goodness-of-fit threshold. We find the mean air wavelength to be 6707.856 Å, with a standard error of 0.002 Å (0.09 km/s) and a weighted standard deviation of 0.026 Å (1.16 km/s). The observed spread in measured positions spans 0.145 Å, or 6.5 km/s, which is up to a factor of six higher than typically reported velocity errors for high-resolution studies. We also find a correlation between the effective temperature of the star and the wavelength of the central absorption. We discuss how exclusively using this Li feature as a reference for photocentric velocity in young stars could potentially be introducing a systematic positive offset in wavelength to measurements of further spectral lines. If outflow tracing forbidden lines, such as [O i] 6300 Å, are actually more blueshifted than previously thought, this then favours a disk wind as the origin for such emission in young stars.

A new model grid containing 228,016 synthetic red supergiant explosions (Type II supernovae) is introduced. Time evolution of spectral energy distributions from 1 A to 50,000 A (100 frequency bins in a log scale) is computed at each time step up to 500 days after explosion in each model. We provide light curves for the filters of the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST), Zwicky Transient Facility (ZTF), Sloan Digital Sky Servey (SDSS), and the Neil Gehrels Swift Observatory, but light curves for any photometric filters can be constructed by convolving any filter response functions to the synthetic spectral energy distributions. We also provide bolometric light curves and photosphere information such as photospheric velocity evolution. The parameter space covered by the model grid is five progenitor masses (10, 12, 14, 16, and 18 Msun at the zero-age main sequence, solar metallicity), ten explosion energies (0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, and 5.0 x 10^51 erg), nine 56Ni masses (0.001, 0.01, 0.02, 0.04, 0.06, 0.08, 0.1, 0.2, and 0.3 Msun), nine mass-loss rates (1e-5.0, 1e-4.5, 1e-4.0, 1e-3.5, 1e-3.0, 1e-2.5, 1e-2.0, 1e-1.5, and 1e-1.0 Msun/yr with a wind velocity of 10 km/s), six circumstellar matter radii (1, 2, 4, 6, 8, and 10 x 10^14 cm), and ten circumstellar structures (beta = 0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, and 5.0). 56Ni is assumed to be uniformly mixed up to the half mass of a hydrogen-rich envelope. This model grid can be a base for rapid characterizations of Type II supernovae with sparse photometric sampling expected in LSST through a Bayesian approach, for example. The model grid is available at doi.org/10.5061/dryad.pnvx0k6sj.

Neutrinos propagating in a dense neutrino gas, such as those expected in core-collapse supernovae (CCSNe) and neutron star mergers (NSMs), can experience fast flavor conversions on relatively short scales. This can happen if the neutrino electron lepton number ($\nu$ELN) angular distribution crosses zero in a certain direction. Despite this, most of the state-of-the-art CCSN and NSM simulations do not provide such detailed angular information and instead, supply only a few moments of the neutrino angular distributions. In this study we employ, for the \emph{first} time, a machine learning (ML) approach to this problem and show that it can be extremely successful in detecting $\nu$ELN crossings on the basis of its zeroth and first moments. We observe that an accuracy of $\sim95\%$ can be achieved by the ML algorithms, which almost corresponds to the Bayes error rate of our problem. Considering its remarkable efficiency and agility, the ML approach provides one with an unprecedented opportunity to evaluate the occurrence of FFCs in CCSN and NSM simulations \emph{on the fly}. We also provide our ML methodologies on \href{https://github.com/sajadabbar/ML-nu_FFI/tree/main}{GitHub}.

Spiral arms are observed in numerous protoplanetary discs. These spiral arms can be excited by companions, either on bound or unbound orbits. We simulate a scenario where an unbound perturber, i.e. a flyby, excites spiral arms during a periastron passage. We run three-dimensional hydrodynamical simulations of a parabolic flyby encountering a gaseous protoplanetary disc. The perturber mass ranges from $10\, \rm M_J$ to $1\, \rm M_{\odot}$. The perturber excites a two-armed spiral structure, with a more prominent spiral feature for higher mass perturbers. The two arms evolve over time, eventually winding up, consistent with previous works. We focus on analysing the pattern speed and pitch angle of these spirals during the whole process. The initial pattern speed of the two arms are close to the angular velocity of the perturber at periastron, and then it decreases over time. The pitch angle also decreases over time as the spiral winds up. The spirals disappear after several local orbital times. An inclined prograde orbit flyby induces similar disc substructures as a coplanar flyby. A solar-mass flyby event causes increased eccentricity growth in the protoplanetary disc, leading to an eccentric disc structure which dampens over time. The spirals' morphology and the disc eccentricity can be used to search for potential unbound stars or planets around discs where a flyby is suspected. Future disc observations at high resolution and dedicated surveys will help to constrain the frequency of such stellar encounters in nearby star-forming regions.

We study, in the context of the three-dimensional ${\cal N}=6$ Chern-Simons-matter (ABJM) theory, the infrared-finite functions that result from performing $L-1$ loop integrations over the $L$-loop integrand of the logarithm of the four-particle scattering amplitude. Our starting point are the integrands obtained from the recently proposed all-loop projected amplituhedron for the ABJM theory. Organizing them in terms of negative geometries ensures that no divergences occur upon integration if at least one loop variable is left unintegrated. We explicitly perform the integrations up to $L=3$, finding both parity-even and -odd terms. Moreover, we discuss a prescription to compute the cusp anomalous dimension $\Gamma_{\rm cusp}$ of ABJM in terms of the integrated negative geometries, and we use it to reproduce the first non-trivial order of $\Gamma_{\rm cusp}$. Finally, we show that the leading singularities that characterize the integrated results are conformally invariant.

We examine the influence of quadrupole moment of a slowly rotating neutron star (NS) on the oscillations of a fluid accretion disk (torus) orbiting a compact object the spacetime around which is described by the Hartle-Thorne geometry. Explicit formulae for non-geodesic orbital epicyclic and precession frequencies, as well as their simplified practical versions that allow for an expeditious application of the universal relations determining the NS properties, are obtained and examined. We demonstrate that the difference in the accretion disk precession frequencies for NSs of the same mass and angular momentum, but different oblateness, can reach up to tens of percent. Even higher differences can arise when NSs with the same mass and rotational frequency, but different equations of state (EoS), are considered. In particular, the Lense-Thirring precession frequency defined in the innermost parts of the accretion region can differ by more than one order of magnitude across NSs with different EoS. Our results have clear implications for models of the LMXBs variability.

Simulating high-resolution detector responses is a storage-costly and computationally intensive process that has long been challenging in particle physics. Despite the ability of deep generative models to make this process more cost-efficient, ultra-high-resolution detector simulation still proves to be difficult as it contains correlated and fine-grained mutual information within an event. To overcome these limitations, we propose Intra-Event Aware GAN (IEA-GAN), a novel fusion of Self-Supervised Learning and Generative Adversarial Networks. IEA-GAN presents a Relational Reasoning Module that approximates the concept of an ''event'' in detector simulation, allowing for the generation of correlated layer-dependent contextualized images for high-resolution detector responses with a proper relational inductive bias. IEA-GAN also introduces a new intra-event aware loss and a Uniformity loss, resulting in significant enhancements to image fidelity and diversity. We demonstrate IEA-GAN's application in generating sensor-dependent images for the high-granularity Pixel Vertex Detector (PXD), with more than 7.5M information channels and a non-trivial geometry, at the Belle II Experiment. Applications of this work include controllable simulation-based inference and event generation, high-granularity detector simulation such as at the HL-LHC (High Luminosity LHC), and fine-grained density estimation and sampling. To the best of our knowledge, IEA-GAN is the first algorithm for faithful ultra-high-resolution detector simulation with event-based reasoning.

We characterize twelve young stellar objects (YSOs) located in the Lupus I region, spatially overlapping with the Upper Centaurus Lupus (UCL) sub-stellar association. The aim of this study is to understand whether the Lupus I cloud has more members than what has been claimed so far in the literature and gain a deeper insight into the global properties of the region. We selected our targets using the Gaia DR2 catalog based on their consistent kinematic properties with the Lupus I bona fide members. In our sample of twelve YSOs observed by X-shooter, we identified ten Lupus I members. We could not determine the membership status of two of our targets, namely Gaia DR2 6014269268967059840 and 2MASS J15361110-3444473 due to technical issues. We found out that four of our targets are accretors, among them, 2MASS J15551027-3455045, with a mass of ∼0.03 M_⊙, is one of the least massive accretors in the Lupus complex identified to date. Several of our targets (including accretors) are formed in situ and off-cloud with respect to the main filaments of Lupus I; hence, our study may hint that there are diffused populations of M dwarfs around Lupus I main filaments. In this context, we would like to emphasize that our kinematic analysis with Gaia catalogs played a key role in identifying the new members of the Lupus I cloud.

Based on observations collected at the European Southern Observatory at Paranal under program 105.20P9.001.

We present for the first time fully analytic results for multi-loop equal-mass ice cone graphs in two dimensions. By analysing the leading singularities of these integrals, we find that the maximal cuts in two dimensions can be organised into two copies of the same periods that describe the Calabi-Yau varieties for the equal-mass banana integrals. We obtain a conjectural basis of master integrals at an arbitrary number of loops, and we solve the system of differential equations satisfied by the master integrals in terms of the same class of iterated integrals that have appeared earlier in the context of equal-mass banana integrals. We then go on and show that, when expressed in terms of the canonical coordinate on the moduli space, our results can naturally be written as iterated integrals involving the geometrical invariants of the Calabi-Yau varieties. Our results indicate how the concept of pure functions and transcendental weight can be extended to the case of Calabi-Yau varieties. Finally, we also obtain a novel representation of the periods of the Calabi-Yau varieties in terms of the same class of iterated integrals, and we show that the well-known quadratic relations among the periods reduce to simple shuffle relations among these iterated integrals.

The prebiotic replication of DNA and RNA is a complex interplay between chemistry and the environment. Factors that have direct and indirect effects on prebiotic chemistry include temperature, concentration of monovalent and bivalent ions, the pH of water, ultraviolet irradiation and the presence of gaseous CO₂. We discuss various primordial conditions to host the first replication reactions on the early Earth, including heated rock pores, hydrothermal vents, evaporating water ponds, freezing-thawing ice compartments, ultraviolet irradiation and high CO₂ concentrations. We review how the interplay of replication chemistry with the strand separation and length selectivity of non-equilibrium physics can be provided by plausible geo-environments. Fast molecular evolution has been observed over a few hours in such settings when a polymerase protein is used as replicator. Such experimental findings make us optimistic that it will soon also be possible to probe evolution dynamics with much slower prebiotic replication chemistries using RNA. Our expectation is that the unique autonomous evolution dynamics provided by microfluidic non-equilibria make the origin of life understandable and experimentally testable in the near future.

Certain Feynman integrals are associated to Calabi-Yau geometries. We demonstrate how these integrals can be computed with the method of differential equations. The four-loop equal-mass banana integral is the simplest Feynman integral whose geometry is a nontrivial Calabi-Yau manifold. We show that its differential equation can be cast into an ϵ -factorized form. This allows us to obtain the solution to any desired order in the dimensional regularization parameter ϵ . The method generalizes to other Calabi-Yau Feynman integrals. Our calculation also shows that the four-loop banana integral is only minimally more complicated than the corresponding Feynman integrals at two or three loops.

We bootstrap the symbol of the maximal-helicity-violating four-particle form factor for the chiral part of the stress-tensor supermultiplet in planar N =4 super-Yang-Mills theory at two loops. When minimally normalized, this symbol involves only 34 letters and obeys the extended Steinmann relations in all partially overlapping three-particle momentum channels. In addition, the remainder function for this form factor exhibits an antipodal self-duality: It is invariant under the combined operation of the antipodal map defined on multiple polylogarithms—which reverses the order of the symbol letters—and a simple kinematic map. This self-duality holds on a four-dimensional parity-preserving kinematic hypersurface. It implies the antipodal duality recently noticed between the three-particle form factor and the six-particle amplitude in this theory.

We provide analytic results for two-loop four-point master integrals with one massive propagator and one massive leg relevant to single top production. Canonical bases of master integrals are constructed and the Simplified Differential Equations approach is employed for their analytic solution. The necessary boundary terms are computed in closed form in the dimensional regulator, allowing us to obtain analytic results in terms of multiple polylogarithms of arbitrary transcendental weight. We provide explicit solutions of all two-loop master integrals up to transcendental weight six and discuss their numerical evaluation for Euclidean and physical phase-space points.

We consider a dense neutrino gas in the "fast-flavor limit" (vanishing neutrino masses). For the first time, we identify exact solutions of the nonlinear wave equation in the form of solitons. They can propagate with both sub- or superluminal speed, the latter not violating causality. The soliton with infinite speed is a homogeneous solution and coincides with the usual fast-flavor pendulum except that it swings only once instead of being periodic. The subluminal soliton in the static limit corresponds to a one-swing "spatial pendulum". A necessary condition for such solutions to exist is a ``crossed'' neutrino angle distribution. Based on the Nyquist criterion, we derive a new sufficient condition without solving the dispersion relation. The solitons are very fragile: they are as unstable as the homogeneous neutrino gas alone. Moreover, in the presence of matter, only the solution survives that is homogeneous in a frame comoving with the matter current. Generally, the matter effect cannot be eliminated by transformations in flavor space, but instead has a real physical impact.

Spiral arms are observed in numerous protoplanetary discs. These spiral arms can be excited by companions, either on bound or unbound orbits. We simulate a scenario where an unbound perturber, i.e. a flyby, excites spiral arms during a periastron passage. We run three-dimensional hydrodynamical simulations of a parabolic flyby encountering a gaseous protoplanetary disc. The perturber mass ranges from $10\, \rm M_J$ to $1\, \rm {\rm M}_{\odot }$. The perturber excites a two-armed spiral structure, with a more prominent spiral feature for higher mass perturbers. The two arms evolve over time, eventually winding up, consistent with previous works. We focus on analysing the pattern speed and pitch angle of these spirals during the whole process. The initial pattern speed of the two arms are close to the angular velocity of the perturber at periastron, and then it decreases over time. The pitch angle also decreases over time as the spiral winds up. The spirals disappear after several local orbital times. An inclined prograde orbit flyby induces similar disc substructures as a coplanar flyby. A solar-mass flyby event causes increased eccentricity growth in the protoplanetary disc, leading to an eccentric disc structure which dampens over time. The spirals' morphology and the disc eccentricity can be used to search for potential unbound stars or planets around discs where a flyby is suspected. Future disc observations at high resolution and dedicated surveys will help to constrain the frequency of such stellar encounters in nearby star-forming regions.

We take a major step towards computing D-dimensional one-loop amplitudes in general gauge theories, compatible with the principles of unitarity and the color-kinematics duality. For n-point amplitudes with either supersymmetry multiplets or generic non-supersymmetric matter in the loop, simple all-multiplicity expressions are obtained for the maximal cuts of kinematic numerators of n-gon diagrams. At n = 6, 7 points with maximal supersymmetry, we extend the cubic-diagram numerators to encode all contact terms, and thus solve the long-standing problem of simultaneously realizing the following properties: color-kinematics duality, manifest locality, optimal power counting of loop momenta, quadratic rather than linearized Feynman propagators, compatibility with double copy as well as all graph symmetries. Color-kinematics dual representations with similar properties are presented in the half-maximally supersymmetric case at n = 4, 5 points. The resulting gauge-theory integrands and their supergravity counterparts obtained from the double copy are checked to reproduce the expected ultraviolet divergences.

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

We investigate the mass-metallicity relationship of star forming galaxies by analysing the absorption line spectra of $\sim$200,000 galaxies in the Sloan Digital Sky Survey. The galaxy spectra are stacked in bins of stellar mass and a population synthesis technique is applied yielding metallicities, ages and star formation history of the young and old stellar population together with interstellar reddening and extinction. We adopt different lengths of the initial starbursts and different initial mass functions for the calculation of model spectra of the single stellar populations contributing to the total integrated spectrum. We also allow for deviations of the ratio of extinction to reddening RV from 3.1 and determine the value from the spectral fit. We find that burst length and RV have a significant influence on the determination of metallicities whereas the effect of the initial mass function is small. RV values are larger than 3.1. The metallicities of the young stellar population agree with extragalactic spectroscopic studies of individual massive supergiant stars and are significantly higher than those of the older stellar population. This confirms galaxy evolution models where metallicity depends on the ratio of gas to stellar mass and where this ratio decreases with time. Star formation history is found to depend on galaxy stellar mass. Massive galaxies are dominated by stars formed at early times.

Polarization of the cosmic microwave background (CMB) is sensitive to new physics violating parity symmetry, such as the presence of a pseudoscalar "axionlike" field. Such a field may be responsible for early dark energy (EDE), which is active prior to recombination and provides a solution to the so-called Hubble tension. The EDE field coupled to photons in a parity-violating manner would rotate the plane of linear polarization of the CMB and produce a cross-correlation power spectrum of $E$- and $B$-mode polarization fields with opposite parities. In this paper, we fit the $EB$ power spectrum predicted by the photon-axion coupling of the EDE model with a potential $V(\phi)\propto [1-\cos(\phi/f)]^3$ to polarization data from Planck. We find that the unique shape of the predicted $EB$ power spectrum is not favored by the data and obtain a first constraint on the photon-axion coupling constant, $g=(0.04\pm 0.16)M_{\text{Pl}}^{-1}$ (68% CL), for the EDE model that best fits the CMB and galaxy clustering data. This constraint is independent of the miscalibration of polarization angles of the instrument or the polarized Galactic foreground emission. Our limit on $g$ may have important implications for embedding EDE in fundamental physics, such as string theory.

The standard perturbation theory (SPT) approach to gravitational clustering is based on a fluid approximation of the underlying Vlasov-Poisson dynamics, taking only the zeroth and first cumulant of the phase-space distribution function into account (density and velocity fields). This assumption breaks down when dark matter particle orbits cross and leads to well-known problems, e.g., an anomalously large backreaction of small-scale modes onto larger scales that compromises predictivity. We extend SPT by incorporating second and higher cumulants generated by orbit crossing. For collisionless matter, their equations of motion are completely fixed by the Vlasov-Poisson system, and thus we refer to this approach as Vlasov Perturbation Theory (VPT). Even cumulants develop a background value, and they enter the hierarchy of coupled equations for the fluctuations. The background values are in turn sourced by power spectra of the fluctuations. The latter can be brought into a form that is formally analogous to SPT, but with an extended set of variables and linear as well as nonlinear terms, that we derive explicitly. In this paper, we focus on linear solutions, which are far richer than in SPT, showing that modes that cross the dispersion scale set by the second cumulant are highly suppressed. We derive stability conditions on the background values of even cumulants from the requirement that exponential instabilities be absent. We also compute the expected magnitude of averaged higher cumulants for various halo models and show that they satisfy the stability conditions. Finally, we derive self-consistent solutions of perturbations and background values for a scaling universe and study the convergence of the cumulant expansion. The VPT framework provides a conceptually straightforward and deterministic extension of SPT that accounts for the decoupling of small-scale modes.

Context. In recent years, a new hot topic has emerged in the star and planet formation field, namely, the interaction between the circumstellar disk and its birth cloud. The birth environments of young stars leave strong imprints on the star itself and their surroundings. In this context, we present a detailed analysis of the rich circumstellar environment around the young Herbig Ae/Be star T CrA.
Aims: Our aim is to understand the nature of the stellar system and the extended circumstellar structures, as seen in scattered light images.
Methods: We conducted our analysis on the basis of a set of combined archival data and new adaptive optics images at a high contrast and high resolution.
Results: The scattered light images reveal the presence of a complex environment around T CrA, composed of a bright, forward-scattering rim of the disk's surface that is seen at very high inclinations, along with a dark lane of the disk midplane, bipolar outflows, and streamer features that are likely tracing infalling material from the surrounding birth cloud onto the disk. The analysis of the light curve suggests that the star is a binary with a period of 29.6 yr, confirming previous assertions based on spectro-astrometry. The comparison of the scattered light images with the ALMA continuum and ¹²CO (2-1) line emission shows that the disk is in Keplerian rotation and the northern side of the outflowing material is receding, while the southern side is approaching the observer. The overall system lies on different geometrical planes. The orbit of the binary star is perpendicular to the outflows and is seen edge on. The disk is itself seen edge-on, with a position angle of ~7°. The direction of the outflows seen in scattered light is in agreement with the direction of the more distant molecular hydrogen emission-line objects (MHOs) associated with the star. Modeling of the spectral energy distribution using a radiative transfer scheme is in good agreement with the proposed configuration, as well as the hydrodynamical simulation performed using a smoothed particle hydrodynamics code.
Conclusions: We find evidence of streamers of accreting material around T CrA. These streamers connect the filament, along which T CrA is forming along with the outer parts of the disk, suggesting that the strong misalignment between the inner and outer disk is due to a change in the direction of the angular momentum of the material accreting on the disk during the late phase of star formation. This impacts the accretion taking place in the components of the binary, favoring the growth of the primary with respect the secondary, in contrast to the case of aligned disks.

Reduced images are also available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (ftp://130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/671/A82

Using numerical simulations, we investigate the gravitational evolution of filamentary molecular cloud structures and their condensation into dense protostellar cores. One possible process is the so called 'edge effect', the pile-up of matter at the end of the filament due to self-gravity. This effect is predicted by theory but only rarely observed. To get a better understanding of the underlying processes we used a simple analytic approach to describe the collapse and the corresponding collapse time. We identify a model of two distinct phases: The first phase is free fall dominated, due to the self-gravity of the filament. In the second phase, after the turning point, the collapse is balanced by the ram pressure, produced by the inside material of the filament, which leads to a constant collapse velocity. This approach reproduces the established collapse time of uniform density filaments and agrees well with our hydrodynamic simulations. In addition, we investigate the influence of different radial density profiles on the collapse. We find that the deviations compared to the uniform filament are less than 10%. Therefore, the analytic collapse model of the uniform density filament is an excellent general approach.

We generalize existing constraints on primordial black holes to dark objects with extended sizes using the aLIGO design sensitivity. We show that LIGO is sensitive to dark objects with radius O (10 −10³ km ) if they make up more than ∼O (10^-2−10^-3) of dark matter.

Enzyme-enriched condensates can organize the spatial distribution of their substrates by catalyzing nonequilibrium reactions. Conversely, an inhomogeneous substrate distribution induces enzyme fluxes through substrate-enzyme interactions. We find that condensates move toward the center of a confining domain when this feedback is weak. Above a feedback threshold, they exhibit self-propulsion, leading to oscillatory dynamics. Moreover, catalysis-driven enzyme fluxes can lead to interrupted coarsening, resulting in equidistant condensate positioning, and to condensate division.

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 allowing for a consistent treatment of uncertainties. 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 geometry optimization enabling to learn optimal detector designs. Lastly, we apply our method to realistic simulation of a ton-scale water-based liquid scintillator detector.

[Abridged] Aims. We provide an important step toward a better understanding of the magnetorotational instability (MRI)-dust coevolution in protoplanetary disks by presenting a proof of concept that dust evolution ultimately plays a crucial role in the MRI activity. Methods. First, we study how a fixed power-law dust size distribution with varying parameters impacts the MRI activity, especially the steady-state MRI-driven accretion, by employing and improving our previous 1+1D MRI-driven turbulence model. Second, we relax the steady-state accretion assumption in this disk accretion model, and partially couple it to a dust evolution model in order to investigate how the evolution of dust (dynamics and grain growth processes combined) and MRI-driven accretion are intertwined on million-year timescales. Results. Dust coagulation and settling lead to a higher gas ionization degree in the protoplanetary disk, resulting in stronger MRI-driven turbulence as well as a more compact dead zone. On the other hand, fragmentation has an opposite effect because it replenishes the disk in small dust particles. Since the dust content of the disk decreases over million years of evolution due to radial drift, the MRI-driven turbulence overall becomes stronger and the dead zone more compact until the disk dust-gas mixture eventually behaves as a grain-free plasma. Furthermore, our results show that dust evolution alone does not lead to a complete reactivation of the dead zone. Conclusions. The MRI activity evolution (hence the temporal evolution of the MRI-induced $\alpha$-parameter) is controlled by dust evolution and occurs on a timescale of local dust growth, as long as there is enough dust particles in the disk to dominate the recombination process for the ionization chemistry. Once it is no longer the case, it is expected to be controlled by gas evolution and occurs on a viscous evolution timescale.

Entropy production is a necessary ingredient for addressing the overpopulation of thermal relics. It is widely employed in particle physics models for explaining the origin of dark matter. A long-lived particle that decays to the known particles, while dominating the universe, plays the role of the dilutor. We point out the impact of its partial decay to dark matter on the primordial matter power spectrum. For the first time, we derive a stringent limit on the branching ratio of the dilutor to dark matter from large scale structure observation using the sloan digital sky survey data. This offers a novel tool for testing models with a dark matter dilution mechanism. We apply it to the left-right symmetric model and show that it firmly excludes a large portion of parameter space for right-handed neutrino warm dark matter.

We report on our study of SN 2022xxf during the first four months of its evolution. The light curves (LCs) display two humps at similar maximum brightness separated by 75d, unprecedented for a broad-lined Type Ic supernova (SN IcBL). SN~2022xxf is the most nearby SN IcBL to date (in NGC~3705, $z = 0.0037$, 20 Mpc). Optical and NIR photometry and spectroscopy are used to identify the energy source powering the LC. Nearly 50 epochs of high S/N-ratio spectroscopy were obtained within 130d, comprising an unparalleled dataset for a SN IcBL, and one of the best-sampled SN datasets to date. The global spectral appearance and evolution of SN~2022xxf points to typical SN Ic/IcBL, with broad features (up to $\sim14000$ km~s$^{-1}$) and a gradual transition from the photospheric to the nebular phase. However, narrow emission lines (corresponding to $\sim1000-2500$ km~s$^{-1}$) are present from the time of the second rise, suggesting slower-moving circumstellar material (CSM). These lines are subtle, but some are readily noticeable at late times such as in Mg~I $\lambda$5170 and [O~I] $\lambda$5577. Unusually, the near-infrared spectra show narrow line peaks, especially among features formed by ions of O and Mg. We infer the presence of CSM that is free of H and He. We propose that the radiative energy from the ejecta-CSM interaction is a plausible explanation for the second LC hump. This interaction scenario is supported by the color evolution, which progresses to the blue as the light curve evolves along the second hump, and the slow second rise and subsequent rapid LC drop. SN~2022xxf may be related to an emerging number of CSM-interacting SNe Ic, which show slow, peculiar LCs, blue colors, and subtle CSM interaction lines. The progenitor stars of these SNe likely experienced an episode of mass loss shortly prior to explosion consisting of H/He-free material.

We calculate chromoelectric and chromomagnetic correlators in quenched QCD at 1.5 T_c and 10⁴T_c , with the aim to estimate the heavy quark diffusion coefficient at leading order in the inverse heavy quark mass expansion, κ_E , as well as the coefficient of the first mass-suppressed correction, κ_B. We use gradient flow for noise reduction. At 1.5 T_c we obtain 1.70 ≤κ_E/T³≤3.12 and 1.03 <κ_B/T³<2.61 . The latter implies that the mass-suppressed effects in the heavy quark diffusion coefficient are 20% for bottom quarks and 34% for charm quarks at this temperature.

We consider electrically neutral complex vector particles $V$ below the GeV mass scale that, from a low energy perspective, couple to the photon via higher dimensional form factor interactions. We derive ensuing astrophysical constraints by considering the anomalous energy loss from the Sun, Horizontal Branch, and Red Giant stars as well as from SN1987A that arise from vector pair-production in these environments. Under the assumption that the dark states $V$ constitute dark matter, the bounds are then complemented by direct and indirect detection as well as cosmological limits. The relic density from freeze-out and freeze-in mechanisms is also computed. On the basis of a UV-complete model that realizes the considered effective couplings, we also discuss the naturalness of the constrained parameter space, and provide an analysis of the zero mass limit of $V$.

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 Λ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 Λ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.

Cosmology inference of galaxy clustering at the field level with the EFT likelihood in principle allows for extracting all non-Gaussian information from quasi-linear scales, while robustly marginalizing over any astrophysical uncertainties. A pipeline in this spirit is implemented in the \texttt{LEFTfield} code, which we extend in this work to describe the clustering of galaxies in redshift space. Our main additions are: the computation of the velocity field in the LPT gravity model, the fully nonlinear displacement of the evolved, biased density field to redshift space, and a systematic expansion of velocity bias. We test the resulting analysis pipeline by applying it to synthetic data sets with a known ground truth at increasing complexity: mock data generated from the perturbative forward model itself, sub-sampled matter particles, and dark matter halos in N-body simulations. By fixing the initial-time density contrast to the ground truth, while varying the growth rate $f$, bias coefficients and noise amplitudes, we perform a stringent set of checks. These show that indeed a systematic higher-order expansion of the velocity bias is required to infer a growth rate consistent with the ground truth within errors. Applied to dark matter halos, our analysis yields unbiased constraints on $f$ at the level of a few percent for a variety of halo masses at redshifts $z=0,\,0.5,\,1$ and for a broad range of cutoff scales $0.08\,h/\mathrm{Mpc} \leq \Lambda \leq 0.20\,h/\mathrm{Mpc}$. Importantly, deviations between true and inferred growth rate exhibit the scaling with halo mass, redshift and cutoff that one expects based on the EFT of Large Scale Structure. Further, we obtain a robust detection of velocity bias through its effect on the redshift-space density field and are able to disentangle it from higher-derivative bias contributions.

New-generation direct searches for low mass dark matter feature detection thresholds at energies well below 100 eV, much lower than the energies of commonly used X-ray calibration sources. This requires new calibration sources with sub-keV energies. When searching for nuclear recoil signals, the calibration source should ideally cause mono-energetic nuclear recoils in the relevant energy range. Recently, a new calibration method based on the radiative neutron capture on $^{182}$W with subsequent de-excitation via single $\gamma$-emission leading to a nuclear recoil peak at 112 eV was proposed. The CRESST-III dark matter search operated several CaWO$_{4}$-based detector modules with detection thresholds below 100 eV in the past years. We report the observation of a peak around the expected energy of 112 eV in the data of three different detector modules recorded while irradiated with neutrons from different AmBe calibration sources. We compare the properties of the observed peaks with Geant-4 simulations and assess the prospects of using this for the energy calibration of CRESST-III detectors.

Collisionless shock waves in supernova remnants and the solar wind heat electrons less effectively than they heat ions, as is predicted by kinetic simulations. However, the values of T$_e$/T$_p$ inferred from the H alpha profiles of supernova remnant shocks behave differently as a function of Mach number or Alfvén Mach number than what is measured in the solar wind or predicted by simulations. Here we determine T$_e$/T$_p$ for supernova remnant shocks using H alpha profiles, shock speeds from proper motions, and electron temperatures from X-ray spectra. We also improve the estimates of sound speed and Alfvén speed used to determine Mach numbers. We find that the H alpha determinations are robust and that the discrepancies among supernova remnant shocks, solar wind shocks and computer-simulated shocks remain. We discuss some possible contributing factors, including shock precursors, turbulence and varying preshock conditions.

Light (anti-) nuclei are a powerful tool both in collider physics and astrophysics. In searches for new and exotic physics, the expected small astrophysical backgrounds at low energies make these antinuclei ideal probes for, e.g., dark matter. At the same time, their composite structure and small binding energies imply that they can be used in collider experiments to probe the hadronisation process and two-particle correlations. For the proper interpretation of such experimental studies, an improved theoretical understanding of (anti-) nuclei production in specific kinematic regions and detector setups is needed. In this work, we develop a coalescence framework for (anti-) deuteron production which accounts for both the emission volume and momentum correlations on an event-by-event basis. This framework goes beyond the equal-time approximation, which has been commonly assumed in femtoscopy experiments and (anti-) nucleus production models until now. Using PYTHIA~8 as an event generator, we find that the equal-time approximation leads to an error of O(10%) in low-energy processes like $\Upsilon$ decays, while the errors are negligible at LHC energies. The framework introduced in this work paves the way for tuning event generators to (anti-) nuclei measurements.

In Paper I we showed that clumps in high-redshift galaxies, having a high star formation rate density (\Sigma_SFR), produce disks with two tracks in the [Fe/H]-[\alpha/Fe] chemical space, similar to that of the Milky Way's (MW's) thin + thick disks. Here we investigate the effect of clumps on the bulge's chemistry. The chemistry of the MW's bulge is comprised of a single track with two density peaks separated by a trough. We show that the bulge chemistry of an N-body + smoothed particle hydrodynamics clumpy simulation also has a single track. Star formation within the bulge is itself in the high-\Sigma_SFR clumpy mode, which ensures that the bulge's chemical track follows that of the thick disk at low [Fe/H] and then extends to high [Fe/H], where it peaks. The peak at low metallicity instead is comprised of a mixture of in-situ stars and stars accreted via clumps. As a result, the trough between the peaks occurs at the end of the thick disk track. We find that the high-metallicity peak dominates near the mid-plane and declines in relative importance with height, as in the MW. The bulge is already rapidly rotating by the end of the clump epoch, with higher rotation at low [\alpha/Fe]. Thus clumpy star formation is able to simultaneously explain the chemodynamic trends of the MW's bulge, thin + thick disks and the Splash.

We investigate shock structures driven by merger events in high-resolution simulations that result in a galaxy with a virial mass M ≈ 10¹² M _⊙. 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 ${ \mathcal M }$ = 2-3 for such internal shocks are in agreement with the spectral indices of the observed ORCs. We estimate that ~5% 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 10¹³ M _⊙ in agreement with the observed stellar masses of the galaxies at the center of ORCs. Such shocks might be detectable with next-generation X-ray instruments like the Line Emission Mapper (LEM).

More than 50% of solar-mass stars form in multiple systems. It is therefore crucial to investigate how multiplicity affects the star and planet formation processes at the protostellar stage. We report continuum and C$^{18}$O (2-1) observations of the VLA 1623-2417 protostellar system at 50 au angular resolution as part of the ALMA Large Program FAUST. The 1.3 mm continuum probes the disks of VLA 1623A, B, and W, and the circumbinary disk of the A1+A2 binary. The C$^{18}$O emission reveals, for the first time, the gas in the disk-envelope of VLA 1623W. We estimate the dynamical mass of VLA 1623W, $M_{\rm dyn}=0.45\pm0.08$ M$_{\odot}$, and the mass of its disk, $M_{\rm disk}\sim6\times10^{-3}$ M$_{\odot}$. C$^{18}$O also reveals streamers that extend up to 1000 au, spatially and kinematically connecting the envelope and outflow cavities of the A1+A2+B system with the disk of VLA 1623W. The presence of the streamers, as well as the spatial ($\sim$1300 au) and velocity ($\sim$2.2 km/s) offset of VLA 1623W suggest that either sources W and A+B formed in different cores, interacting between them, or that source W has been ejected from the VLA 1623 multiple system during its formation. In the latter case, the streamers may funnel material from the envelope and cavities of VLA 1623AB onto VLA 1623W, thus concurring to set its final mass and chemical content.

Modeling of strong gravitational lenses is a necessity for further applications in astrophysics and cosmology. With the large number of detections in current and upcoming surveys, such as the Rubin Legacy Survey of Space and Time (LSST), it is pertinent to investigate automated and fast analysis techniques beyond the traditional and time-consuming Markov chain Monte Carlo sampling methods. Building upon our (simple) convolutional neural network (CNN), we present here another CNN, specifically a residual neural network (ResNet), that predicts the five mass parameters of a singular isothermal ellipsoid (SIE) profile (lens center x and y, ellipticity e_x and e_y, Einstein radius θ_E) and the external shear (γ_{ext, 1}, γ_{ext, 2}) from ground-based imaging data. In contrast to our previous CNN, this ResNet further predicts the 1σ uncertainty for each parameter. To train our network, we use our improved pipeline to simulate lens images using real images of galaxies from the Hyper Suprime-Cam Survey (HSC) and from the Hubble Ultra Deep Field as lens galaxies and background sources, respectively. We find very good recoveries overall for the SIE parameters, especially for the lens center in comparison to our previous CNN, while significant differences remain in predicting the external shear. From our multiple tests, it appears that most likely the low ground-based image resolution is the limiting factor in predicting the external shear. Given the run time of milli-seconds per system, our network is perfectly suited to quickly predict the next appearing image and time delays of lensed transients. Therefore, we use the network-predicted mass model to estimate these quantities and compare to those values obtained from our simulations. Unfortunately, the achieved precision allows only a first-order estimate of time delays on real lens systems and requires further refinement through follow-up modeling. Nonetheless, our ResNet is able to predict the SIE and shear parameter values in fractions of a second on a single CPU, meaning that we are able to efficiently process the huge amount of galaxy-scale lenses expected in the near future.

The network code is available under https://github.com/shsuyu/HOLISMOKES-public/tree/main/HOLISMOKES_IX

We present nonlinear solutions of Vlasov perturbation theory (VPT), describing gravitational clustering of collisionless dark matter with dispersion and higher cumulants induced by orbit crossing. We show that VPT can be cast into a form that is formally analogous to standard perturbation theory (SPT), but including additional perturbation variables, nonlinear interactions, and a more complex propagation. VPT nonlinear kernels have a crucial decoupling property: for fixed total momentum, the kernels become strongly suppressed when any of the individual momenta cross the dispersion scale into the nonlinear regime. This screening of UV modes allows us to compute nonlinear corrections to power spectra even for cosmologies with very blue power-law input spectra, for which SPT diverges. We compare predictions for the density and velocity divergence power spectra as well as the bispectrum at one-loop order to N -body results in a scaling universe with spectral indices -1 ≤n_s≤+2 . We find a good agreement up to the nonlinear scale for all cases, with a reach that increases with the spectral index n_s. We discuss the generation of vorticity as well as vector and tensor modes of the velocity dispersion, showing that neglecting vorticity when including dispersion would lead to a violation of momentum conservation. We verify momentum conservation when including vorticity, and compute the vorticity power spectrum at two-loop order, necessary to recover the correct large-scale limit with slope n_w=2 . Comparing to our N -body measurements confirms the cross-over from k⁴ to k² scaling at large scales. Our results provide a proof-of-principle that perturbative techniques for dark matter clustering can be systematically improved based on the known underlying collisionless dynamics.

We analyse the full shape of anisotropic clustering measurements from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) quasar sample together with the combined galaxy sample from the Baryon Oscillation Spectroscopic Survey (BOSS). We obtain constraints on the cosmological parameters independent of the Hubble parameter h for the extensions of the ΛCDM models, focusing on cosmologies with free dark energy equation of state parameter w. We combine the clustering constraints with those from the latest CMB data from Planck to obtain joint constraints for these cosmologies for w and the additional extension parameters - its time evolution w_a, the physical curvature density ω_K and the neutrino mass sum ∑m_ν. Our joint constraints are consistent with a flat ΛCDM cosmological model within 68% confidence limits. We demonstrate that the Planck data are able to place tight constraints on the clustering amplitude today, σ₁₂, in cosmologies with varying w and present the first constraints for the clustering amplitude for such cosmologies, which is found to be slightly higher than the ΛCDM value. Additionally, we show that when we vary w and allow for non-flat cosmologies and the physical curvature density is used, Planck prefers a curved universe at 4σ significance, which is ~2σ higher than when using the relative curvature density Ω_K. Finally, when w is varied freely, clustering provides only a modest improvement (of 0.021 eV) on the upper limit of ∑m_ν.

VISIONS is an ESO public survey of five nearby (d < 500 pc) star-forming molecular cloud complexes that are canonically associated with the constellations of Chamaeleon, Corona Australis, Lupus, Ophiuchus, and Orion. The survey was carried out with VISTA, using VIRCAM, and collected data in the near-infrared passbands J, H, and Ks. With a total on-sky exposure time of 49.4 h VISIONS covers an area of 650 deg$^2$, and it was designed to build an infrared legacy archive similar to that of 2MASS. Taking place between April 2017 and March 2022, the observations yielded approximately 1.15 million images, which comprise 19 TB of raw data. The observations are grouped into three different subsurveys: The wide subsurvey comprises shallow, large-scale observations and has visited the star-forming complexes six times over the course of its execution. The deep subsurvey of dedicated high-sensitivity observations has collected data on the areas with the largest amounts of dust extinction. The control subsurvey includes observations of areas of low-to-negligible dust extinction. Using this strategy, the VISIONS survey offers multi-epoch position measurements, is able to access deeply embedded objects, and provides a baseline for statistical comparisons and sample completeness. In particular, VISIONS is designed to measure the proper motions of point sources with a precision of 1 mas/yr or better, when complemented with data from VHS. Hence, VISIONS can provide proper motions for sources inaccessible to Gaia. VISIONS will enable addressing a range of topics, including the 3D distribution and motion of embedded stars and the nearby interstellar medium, the identification and characterization of young stellar objects, the formation and evolution of embedded stellar clusters and their initial mass function, as well as the characteristics of interstellar dust and the reddening law.

The direct detection of core-collapse supernova (SN) progenitor stars is a powerful way of probing the last stages of stellar evolution. However, detections in archival Hubble Space Telescope images are limited to about one per year. Here, we explore whether we can increase the detection rate by using data from ground-based wide-field surveys. Due to crowding and atmospheric blurring, progenitor stars can typically not be identified in pre-explosion images alone. Instead, we combine many pre-SN and late-time images to search for the disappearance of the progenitor star. As a proof of concept, we implement our search for ZTF data. For a few hundred images, we achieve limiting magnitudes of about 23 mag in the g and r band. However, no progenitor stars or long-lived outbursts are detected for 29 SNe within z<0.01, and the ZTF limits are typically several magnitudes less constraining than detected progenitors in the literature. Next, we estimate progenitor detection rates for the Legacy Survey of Space and Time (LSST) with the Vera C. Rubin telescope by simulating a population of nearby SNe. The background from bright host galaxies reduces the nominal LSST sensitivity by, on average, 0.4 mag. Over the ten-year survey, we expect the detection of about 50 red supergiant progenitors and several yellow and blue supergiants. The progenitors of SNe Ib and Ic are detectable if they are brighter than -4.7 mag or -4.0 mag in the LSST i band, respectively. In addition, we expect the detection of hundreds of pre-SN outbursts depending on their brightness and duration.

We present new computations for Feynman integrals relevant to Higgs plus jet production at three loops, including first results for a non-planar class of integrals. The results are expressed in terms of generalised polylogarithms up to transcendental weight six. We also provide the full canonical differential equations, which allows us to make structural observations on the answer. In particular, we find a counterexample to previously conjectured adjacency relations, for a planar integral of the tennis-court type. Additionally, for a non-planar triple ladder diagram, we find two novel alphabet letters. This information may be useful for future bootstrap approaches.

We analyze Hα or CO rotation curves extending out to several galaxy effective radii for 100 massive, large, star-forming disk galaxies (SFGs) across the peak of cosmic galaxy star formation (z ~ 0.6-2.5), more than doubling the previous sample presented by Genzel et al. and Price et al. The observations were taken with SINFONI and KMOS integral-field spectrographs at the ESO-Very Large Telescope, LUCI-LBT, NOEMA-IRAM, and Atacama Large Millimeter/submillimeter Array. We fit the major-axis kinematics with beam-convolved, forward models of turbulent rotating disks with bulges embedded in dark matter (DM) halos, including the effects of pressure support. The fraction of dark to total matter within the disk effective radius (R _e ~ 5 kpc), f _DM(R _e) = V ² _DM(R _e)/V ² _circ(R _e) decreases with redshift: at z ~ 1 (z ~ 2) the median DM fraction is 0.38 ± 0.23 (0.27 ± 0.18), and a third (half) of all galaxies are maximal disks with f _DM(R _e) < 0.28. DM fractions correlate inversely with the baryonic surface density, and the low DM fractions can be explained with a flattened, or cored, inner DM density distribution. At z ~ 2, there is ≈40% less DM mass on average within R _e compared to expected values based on cosmological stellar-mass-halo-mass relations. The DM deficit is more evident at high star formation rate surface densities (≳2.5 M _⊙ yr^-1 kpc²) and galaxies with massive bulges (≥10¹⁰ M _⊙). A combination of stellar or active galactic nucleus feedback, and/or heating due to dynamical friction, may drive the DM from cuspy into cored mass distributions, pointing to an efficient buildup of massive bulges and central black holes at z ~ 2 SFGs.

We present a calculation of all matching coefficients for N-jettiness beam functions at next-to-next-to-next-to-leading order (N³LO) in perturbative quantum chromodynamics (QCD). Our computation is performed starting from the respective collinear splitting kernels, which we integrate using the axial gauge. We use reverse unitarity to map the relevant phase-space integrals to loop integrals, which allows us to employ multi-loop techniques including integration-by-parts identities and differential equations. We find a canonical basis and use an algorithm to establish non-trivial partial fraction relations among the resulting master integrals, which allows us to reduce their number substantially. By use of regularity conditions, we express all necessary boundary constants in terms of an independent set, which we compute by direct integration of the corresponding integrals in the soft limit. In this way, we provide an entirely independent calculation of the matching coefficients which were previously computed in ref. [1].

We discover analytic equations that can infer the value of $\Omega_{\rm m}$ from the positions and velocity moduli of halo and galaxy catalogues. The equations are derived by combining a tailored graph neural network (GNN) architecture with symbolic regression. We first train the GNN on dark matter halos from Gadget N-body simulations to perform field-level likelihood-free inference, and show that our model can infer $\Omega_{\rm m}$ with $\sim6\%$ accuracy from halo catalogues of thousands of N-body simulations run with six different codes: Abacus, CUBEP$^3$M, Gadget, Enzo, PKDGrav3, and Ramses. By applying symbolic regression to the different parts comprising the GNN, we derive equations that can predict $\Omega_{\rm m}$ from halo catalogues of simulations run with all of the above codes with accuracies similar to those of the GNN. We show that by tuning a single free parameter, our equations can also infer the value of $\Omega_{\rm m}$ from galaxy catalogues of thousands of state-of-the-art hydrodynamic simulations of the CAMELS project, each with a different astrophysics model, run with five distinct codes that employ different subgrid physics: IllustrisTNG, SIMBA, Astrid, Magneticum, SWIFT-EAGLE. Furthermore, the equations also perform well when tested on galaxy catalogues from simulations covering a vast region in parameter space that samples variations in 5 cosmological and 23 astrophysical parameters. We speculate that the equations may reflect the existence of a fundamental physics relation between the phase-space distribution of generic tracers and $\Omega_{\rm m}$, one that is not affected by galaxy formation physics down to scales as small as $10~h^{-1}{\rm kpc}$.

The flavor evolution of a neutrino gas can show "slow" or "fast" collective motion. In terms of the usual Bloch vectors to describe the mean-field density matrices of a homogeneous neutrino gas, the slow two-flavor equations of motion (EOMs) are P_˙ω=(ω B +μ P )×P_ω, where ω =Δ m²/2 E , μ =√{2 }G_F(n_ν+n_{ν ¯}), B is a unit vector in the mass direction in flavor space, and P =∫d ω P_ω. For an axisymmetric angle distribution, the fast EOMs are D_˙v=μ (D₀-v D₁)×D_v, where D_v is the Bloch vector for lepton number, v =cos θ is the velocity along the symmetry axis, D₀=∫d v D_v, and D₁=∫d v v D_v. We discuss similarities and differences between these generic cases. Both systems can have pendulumlike instabilities (soliton solutions), both have similar Gaudin invariants, and both are integrable in the classical and quantum case. Describing fast oscillations in a frame comoving with D₁ (which itself may execute pendulumlike motions) leads to transformed EOMs that are equivalent to an abstract slow system. These conclusions carry over to three flavors.

A primary new capability of JWST is the ability to penetrate the dust in star-forming galaxies to identify and study the properties of young star clusters that remain embedded in dust and gas. In this Letter we combine new infrared images taken with JWST with our optical Hubble Space Telescope (HST) images of the starbursting barred (Seyfert2) spiral galaxy NGC 1365. We find that this galaxy has the richest population of massive young clusters of any known galaxy within 30 Mpc, with ~30 star clusters that are more massive than 10⁶ M _⊙ and younger than 10 Myr. Sixteen of these clusters are newly discovered from our JWST observations. An examination of the optical images reveals that 4 of 30 (~13%) are so deeply embedded that they cannot be seen in the Hubble I band (A _V ≳ 10 mag), and that 11 of 30 (~37%) are missing in the HST B band, so age and mass estimates from optical measurements alone are challenging. These numbers suggest that massive clusters in NGC 1365 remain completely obscured in the visible for ~1.3 ± 0.7 Myr and are either completely or partially obscured for ~3.7 ± 1.1 Myr. We also use the JWST observations to gain new insights into the triggering of star cluster formation by the collision of gas and dust streamers with gas and dust in the bar. The JWST images reveal previously unknown structures (e.g., bridges and overshoot regions from stars that form in the bar) that help us better understand the orbital dynamics of barred galaxies and associated star-forming rings. Finally, we note that the excellent spatial resolution of the NIRCAM F200W filter provides a better way to separate barely resolved compact clusters from individual stars based on their sizes.

With the advent of high-cadence, all-sky automated surveys, supernovae (SNe) are now discovered closer than ever to their dates of explosion. However, young premaximum light follow-up spectra of Type Ic SNe (SNe Ic), probably arising from the most-stripped massive stars, remain rare despite their importance. In this Letter, we present a set of 49 optical spectra observed with the Las Cumbres Observatory through the Global Supernova Project for 6 SNe Ic, including a total of 17 premaximum spectra, of which 8 are observed more than a week before V-band maximum light. This data set increases the total number of publicly available premaximum-light SN Ic spectra by 25%, and we provide publicly available SNID templates that will significantly aid in the fast identification of young SNe Ic in the future. We present a detailed analysis of these spectra, including Fe II 5169 velocity measurements, O I 7774 line strengths, and continuum shapes. We compare our results to published samples of stripped SNe in the literature and find one SN in our sample that stands out. SN 2019ewu has a unique combination of features for an SN Ic: an extremely blue continuum, high absorption velocities, a P Cygni-shaped feature almost 2 weeks before maximum light that TARDIS radiative transfer modeling attributes to C II rather than Hα, and weak or nonexistent O I 7774 absorption feature until maximum light.

Heavy QCD axions are well-motivated extensions of the QCD axion that address the quality problem while still solving the strong CP problem. Owing to the gluon coupling, critical for solving the strong CP problem, these axions can be produced in significant numbers in beam dump and collider environments for axion decay constants as large as PeV, relevant for addressing the axion quality problem. In addition, if these axions have leptonic couplings, they can give rise to long-lived decay into lepton pairs, in particular, dominantly into muons above the dimuon threshold and below the GeV scale in a broad class of axion models. Considering existing constraints, primarily from rare meson decays, we demonstrate that current and future neutrino facilities and long-lived particle searches have the potential to probe significant parts of the heavy QCD axion parameter space via dimuon final states.

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 O (G²) 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 O (G²v⁴). 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.

JWST observations of polycyclic aromatic hydrocarbon (PAH) emission provide some of the deepest and highest resolution views of the cold interstellar medium (ISM) in nearby galaxies. If PAHs are well mixed with the atomic and molecular gas and illuminated by the average diffuse interstellar radiation field, PAH emission may provide an approximately linear, high-resolution, high-sensitivity tracer of diffuse gas surface density. We present a pilot study that explores using PAH emission in this way based on Mid-Infrared Instrument observations of IC 5332, NGC 628, NGC 1365, and NGC 7496 from the Physics at High Angular resolution in Nearby GalaxieS-JWST Treasury. Using scaling relationships calibrated in Leroy et al., scaled F1130W provides 10-40 pc resolution and 3σ sensitivity of Σ_gas ~ 2 M _⊙ pc^-2. We characterize the surface densities of structures seen at <7 M _⊙ pc^-2 in our targets, where we expect the gas to be H I-dominated. We highlight the existence of filaments, interarm emission, and holes in the diffuse ISM at these low surface densities. Below ~10 M _⊙ pc^-2 for NGC 628, NGC 1365, and NGC 7496 the gas distribution shows a "Swiss cheese"-like topology due to holes and bubbles pervading the relatively smooth distribution of the diffuse ISM. Comparing to recent galaxy simulations, we observe similar topology for the low-surface-density gas, though with notable variations between simulations with different setups and resolution. Such a comparison of high-resolution, low-surface-density gas with simulations is not possible with existing atomic and molecular gas maps, highlighting the unique power of JWST maps of PAH emission.

We present a high-resolution view of bubbles within the Phantom Galaxy (NGC 628), a nearby (~10 Mpc), star-forming (~2 M _⊙ yr^-1), face-on (i ~ 9°) grand-design spiral galaxy. With new data obtained as part of the Physics at High Angular resolution in Nearby GalaxieS (PHANGS)-JWST treasury program, we perform a detailed case study of two regions of interest, one of which contains the largest and most prominent bubble in the galaxy (the Phantom Void, over 1 kpc in diameter), and the other being a smaller region that may be the precursor to such a large bubble (the Precursor Phantom Void). When comparing to matched-resolution Hα observations from the Hubble Space Telescope, we see that the ionized gas is brightest in the shells of both bubbles, and is coincident with the youngest (~1 Myr) and most massive (~10⁵ M _⊙) stellar associations. We also find an older generation (~20 Myr) of stellar associations is present within the bubble of the Phantom Void. From our kinematic analysis of the H I, H₂ (CO), and H II gas across the Phantom Void, we infer a high expansion speed of around 15 to 50 km s^-1. The large size and high expansion speed of the Phantom Void suggest that the driving mechanism is sustained stellar feedback due to multiple mechanisms, where early feedback first cleared a bubble (as we observe now in the Precursor Phantom Void), and since then supernovae have been exploding within the cavity and have accelerated the shell. Finally, comparison to simulations shows a striking resemblance to our JWST observations, and suggests that such large-scale, stellar-feedback-driven bubbles should be common within other galaxies.

We compute the two-loop corrections to the helicity amplitudes for the coupling of a massive vector boson to a massless quark-antiquark pair and a gluon, accounting for vector and axial-vector couplings of the vector boson and distinguishing isospin non-singlet and singlet contributions. A new four-dimensional basis for the decomposition of the amplitudes into 12 invariant tensor structures is introduced. The associated form factors are then computed up to two loops in QCD using dimensional regularization. After performing renormalization and infrared subtraction, the finite parts of the renormalized non-singlet vector and axial-vector form factors are shown agree with each other, and to reproduce the previously known two-loop amplitudes. The singlet axial-vector amplitude receives a contribution from the axial anomaly from two loops onwards. This amplitude is computed for massless and massive internal quarks. Our results provide the last missing two-loop amplitudes entering the NNLO QCD corrections of vector-boson-plus-jet production at hadron colliders.

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 Letter, 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 to 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.

Large-scale bars can fuel galaxy centers with molecular gas, often leading to the development of dense ringlike structures where intense star formation occurs, forming a very different environment compared to galactic disks. We pair ~0.″3 (30 pc) resolution new JWST/MIRI imaging with archival ALMA CO(2-1) mapping of the central ~5 kpc of the nearby barred spiral galaxy NGC 1365 to investigate the physical mechanisms responsible for this extreme star formation. The molecular gas morphology is resolved into two well-known bright bar lanes that surround a smooth dynamically cold gas disk (R _gal ~ 475 pc) reminiscent of non-star-forming disks in early-type galaxies and likely fed by gas inflow triggered by stellar feedback in the lanes. The lanes host a large number of JWST-identified massive young star clusters. We find some evidence for temporal star formation evolution along the ring. The complex kinematics in the gas lanes reveal strong streaming motions and may be consistent with convergence of gas streamlines expected there. Indeed, the extreme line widths are found to be the result of inter-"cloud" motion between gas peaks; SCOUSEPY decomposition reveals multiple components with line widths of <σ _CO,scouse> ≈ 19 km s^-1 and surface densities of $\langle \,{{\rm{\Sigma }}}_{{{\rm{H}}}_{2},\mathrm{scouse}}\rangle \,\approx \,800\,{M}_{\odot }\,{\mathrm{pc}}^{-2}$ , similar to the properties observed throughout the rest of the central molecular gas structure. Tailored hydrodynamical simulations exhibit many of the observed properties and imply that the observed structures are transient and highly time-variable. From our study of NGC 1365, we conclude that it is predominantly the high gas inflow triggered by the bar that is setting the star formation in its CMZ.

Much like passive materials, active systems can be affected by the presence of imperfections in their microscopic order, called defects, that influence macroscopic properties. This suggests the possibility to steer collective patterns by introducing and controlling defects in an active system. Here we show that a self-assembled, passive nematic is ideally suited to control the pattern formation process of an active fluid. To this end, we force microtubules to glide inside a passive nematic material made from actin filaments. The actin nematic features self-assembled half-integer defects that steer the active microtubules and lead to the formation of macroscopic polar patterns. Moreover, by confining the nematic in circular geometries, chiral loops form. We find that the exact positioning of nematic defects in the passive material deterministically controls the formation and the polarity of the active flow, opening the possibility of efficiently shaping an active material using passive defects.

In the first paper of this series, we proposed a model-independent framework for characterising the architecture of planetary systems at the system level. There are four classes of planetary system architecture: similar, mixed, anti-ordered, and ordered. In this paper, we investigate the formation pathways leading to these four architecture classes. To understand the role of nature versus nurture in sculpting the final (mass) architecture of a system, we apply our architecture framework to synthetic planetary systems - formed via core-accretion - using the Bern model. General patterns emerge in the formation pathways of the four architecture classes. Almost all planetary systems emerging from protoplanetary disks whose initial solid mass was less than one Jupiter mass are similar. Systems emerging from heavier disks may become mixed, anti-ordered, or ordered. Increasing dynamical interactions (planet-planet, planet-disk) tends to shift a system's architecture from mixed to anti-ordered to ordered. Our model predicts the existence of a new metallicity-architecture correlation. Similar systems have very high occurrence around low-metallicity stars. The occurrence of the anti-ordered and ordered classes increases with increasing metallicity. The occurrence of mixed architecture first increases and then decreases with increasing metallicity. In our synthetic planetary systems, the role of nature is disentangled from the role of nurture. Nature (or initial conditions) pre-determines whether the architecture of a system becomes similar; otherwise nurture influences whether a system becomes mixed, anti-ordered, or ordered. We propose the `Aryabhata formation scenario' to explain some planetary systems which host only water-rich worlds. We finish this paper with a discussion of future observational and theoretical works that may support or refute the results of this paper.

Using a new sample of extremely metal poor systems, the EMPRESS survey has recently reported a primordial helium abundance that is 3 σ smaller than the prediction from the standard big bang nucleosynthesis (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 cosmic microwave background (CMB) data from Planck. Assuming that there is no dark radiation in our Universe, we find an electron neutrino chemical potential ξ_νe=0.043 ±0.015 , which deviates from zero by 2.9 σ . If no assumption is made on the abundance of dark radiation in the Universe, the chemical potential is ξ_νe=0.046 ±0.021 , which deviates from zero by 2.2 σ . 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 nonzero lepton asymmetry to 4 σ and 5 σ respectively, assuming no dark radiation, or to 3 σ when no assumption is made on the abundance of dark radiation.

We train graph neural networks on halo catalogs from Gadget N-body simulations to perform field-level likelihood-free inference of cosmological parameters. The catalogs contain ≲5000 halos with masses ≳10¹⁰ h ^-1 M _⊙ in a periodic volume of ${(25\,{h}^{-1}\,\mathrm{Mpc})}^{3}$ ; every halo in the catalog is characterized by several properties such as position, mass, velocity, concentration, and maximum circular velocity. Our models, built to be permutationally, translationally, and rotationally invariant, do not impose a minimum scale on which to extract information and are able to infer the values of Ω_m and σ ₈ with a mean relative error of ~6%, when using positions plus velocities and positions plus masses, respectively. More importantly, we find that our models are very robust: they can infer the value of Ω_m and σ ₈ when tested using halo catalogs from thousands of N-body simulations run with five different N-body codes: Abacus, CUBEP³M, Enzo, PKDGrav3, and Ramses. Surprisingly, the model trained to infer Ω_m also works when tested on thousands of state-of-the-art CAMELS hydrodynamic simulations run with four different codes and subgrid physics implementations. Using halo properties such as concentration and maximum circular velocity allow our models to extract more information, at the expense of breaking the robustness of the models. This may happen because the different N-body codes are not converged on the relevant scales corresponding to these parameters.

We present initial results of an ongoing survey with the Karl G. Jansky Very Large Array targeting the CO($J$ = 1-0) transition in a sample of 30 submillimeter-selected, dusty star-forming galaxies at $z =$ 2-5 with existing mid--$J$ CO detections from ALMA and NOEMA, of which 17 have been fully observed. We detect CO(1-0) emission in 11 targets, along with three tentative ($\sim$1.5-2$\sigma$) detections; three galaxies are undetected. Our results yield total molecular gas masses of 6-23$\times$10$^{10}$ ($\alpha_\mathrm{CO}$/1) M$_\odot$, with gas mass fractions, $f_\mathrm{gas}$=$M_\mathrm{mol}$/($M_*$+$M_\mathrm{mol}$), of 0.1-0.8 and a median depletion time of (140$\pm$70) Myr. We find median CO excitation ratios of $r_{31}$ = 0.75$\pm$0.39 and $r_{41}$ = 0.63$\pm$0.44, with a significant scatter. We find no significant correlation between the excitation ratio and a number of key parameters such as redshift, CO(1-0) line width or $\Sigma_\mathrm{SFR}$. We only find a tentative positive correlation between $r_{41}$ and the star-forming efficiency, but we are limited by our small sample size. Finally, we compare our results to predictions from the SHARK semi-analytical model, finding a good agreement between the molecular gas masses, depletion times and gas fractions of our sources and their SHARK counterparts. Our results highlight the heterogeneous nature of the most massive star-forming galaxies at high-redshift, and the importance of CO(1--0) observations to robustly constrain their total molecular gas content and ISM properties.

Clusters of galaxies are sensitive to the most nonlinear peaks in the cosmic density field. The weak gravitational lensing of background galaxies by clusters can allow us to infer their masses. However, galaxies associated with the local environment of the cluster can also be intrinsically aligned due to the local tidal gradient, contaminating any cosmology derived from the lensing signal. We measure this intrinsic alignment in Dark Energy Survey (DES) Year 1 redMaPPer clusters. We find evidence of a non-zero mean radial alignment of galaxies within clusters between redshift 0.1-0.7. We find a significant systematic in the measured ellipticities of cluster satellite galaxies that we attribute to the central galaxy flux and other intracluster light. We attempt to correct this signal, and fit a simple model for intrinsic alignment amplitude ($A_{\textrm{IA}}$) to the measurement, finding $A_{\textrm{IA}}=0.15\pm 0.04$, when excluding data near the edge of the cluster. We find a significantly stronger alignment of the central galaxy with the cluster dark matter halo at low redshift and with higher richness and central galaxy absolute magnitude (proxies for cluster mass). This is an important demonstration of the ability of large photometric data sets like DES to provide direct constraints on the intrinsic alignment of galaxies within clusters. These measurements can inform improvements to small-scale modeling and simulation of the intrinsic alignment of galaxies to help improve the separation of the intrinsic alignment signal in weak lensing studies.

Fast Radio Bursts (FRBs) are short astrophysical transients of extragalactic origin. Their burst signal is dispersed by the free electrons in the large-scale-structure (LSS), leading to delayed arrival times at different frequencies. Another potential source of time delay is the well known Shapiro delay, which measures the space-space and time-time metric perturbations along the line-of-sight. If photons of different frequencies follow different trajectories, i.e. if the universality of free fall guaranteed by the weak equivalence principle (WEP) is violated, they would experience an additional relative delay. This quantity, however, is not an observable on the background level as it is not gauge independent, which has led to confusion in previous papers. Instead, an imprint can be seen in the correlation between the time delays of different pulses. In this paper, we derive robust and consistent constraints from twelve localised FRBs on the violation of the WEP in the energy range between 4.6 and 6 meV. In contrast to a number of previous studies, we consider our signal to be not in the model, but in the covariance matrix of the likelihood. To do so, we calculate the covariance of the time delays induced by the free electrons in the LSS, the WEP breaking terms, the Milky Way and host galaxy. By marginalising over both host galaxy contribution and the contribution from the free electrons, we find that the parametrised post-Newtonian parameter $\gamma$ characterising the WEP violation must be constant in this energy range to 1 in $10^{13}$ at 68$\;\%$ confidence. These are the tightest constraints to-date on $\Delta\gamma$ in this low energy range.

Mergers of galaxy clusters are promising probes of dark matter (DM) physics. For example, an offset between the DM component and the galaxy distribution can constrain DM self-interactions. We investigate the role of the intracluster medium (ICM) and its influence on DM-galaxy offsets in self-interacting dark matter (SIDM) models. To this end, we employ Smoothed Particle Hydrodynamics + N-body simulations to study idealised setups of equal- and unequal-mass mergers with head-on collisions. Our simulations show that the ICM hardly affects the offsets arising shortly after the first pericentre passage compared to DM-only (DMO) simulations. But later on, e.g. at the first apocentre, the offsets can be amplified by the presence of the ICM. Furthermore, we find that cross-sections small enough not to be excluded by measurements of the core sizes of relaxed galaxy clusters have a chance to produce observable offsets. We found that different DM models affect the DM distribution and also the galaxy and and ICM distribution, including its temperature. Potentially, the position of the shock fronts, combined with the brightest cluster galaxies (BCGs), provides further clues to the properties of DM. Overall our results demonstrate that mergers of galaxy clusters at stages about the first apocentre passage could be more interesting in terms of DM physics than those shortly after the first pericentre passage. This may motivate further studies of mergers at later evolutionary stages.

We use the BOSS DR12 galaxy power spectrum to constrain compensated isocurvature perturbations (CIP), which are opposite-sign primordial baryon and dark matter perturbations that leave the total matter density unchanged. Long-wavelength CIP $\sigma(\vec{x})$ enter the galaxy density contrast as $\delta_g(\vec{x}) \supset b_\sigma\sigma(\vec{x})$, with $b_\sigma$ the linear CIP galaxy bias parameter. We parameterize the CIP spectra as $P_{\sigma\sigma} = A^2P_{\mathcal{R}\mathcal{R}}$ and $P_{\sigma\mathcal{R}} = \xi\sqrt{P_{\sigma\sigma}P_{\mathcal{R}\mathcal{R}}}$, where $A$ is the CIP amplitude and $\xi$ is the correlation with the curvature perturbations $\mathcal{R}$. We find a significance of detection of $Ab_\sigma \neq 0$ of $1.8\sigma$ for correlated ($\xi = 1$) and $3.7\sigma$ for uncorrelated ($\xi = 0$) CIP. Large-scale data systematics have a bigger impact for uncorrelated CIP, which may explain the large significance of detection. The constraints on $A$ depend on the assumed priors for the $b_\sigma$ parameter, which we estimate using separate universe simulations. Assuming $b_\sigma$ values representative of all halos we find $\sigma_A = 145$ for correlated CIP and $\sigma_{|A|} = 475$ for uncorrelated CIP. Our strongest uncorrelated CIP constraint is for $b_\sigma$ representative of the $33\%$ most concentrated halos, $\sigma_{|A|} = 197$, which is better than the current CMB bounds $|A| \lesssim 360$. We also discuss the impact of the local primordial non-Gaussianity parameter $f_{\rm NL}$ in CIP constraints. Our results demonstrate the power of galaxy data to place tight constraints on CIP, and motivate works to understand better the impact of data systematics, as well as to determine theory priors for $b_\sigma$.

Many essential building blocks of life, including amino acids, sugars, and nucleosides, require

aldehydes for prebiotic synthesis. Pathways for their formation under early earth conditions

are therefore of great importance. We investigated the formation of aldehydes by an

experimental simulation of primordial early earth conditions, in line with the metal-sulfur

world theory in an acetylene-containing atmosphere. We describe a pH-driven, intrinsically

autoregulatory environment that concentrates acetaldehyde and other higher molecular

weight aldehydes. We demonstrate that acetaldehyde is rapidly formed from acetylene over a

nickel sulfide catalyst in an aqueous solution, followed by sequential reactions progressively

increasing the molecular diversity and complexity of the reaction mixture. Interestingly,

through inherent pH changes, the evolution of this complex matrix leads to auto-stabilization

of de novo synthesized aldehydes and alters the subsequent synthesis of relevant biomo-

lecules rather than yielding uncontrolled polymerization products. Our results emphasize the

impact of progressively generated compounds on the overall reaction conditions and

strengthen the role of acetylene in forming essential building blocks that are fundamental for

the emergence of terrestrial life.

We investigate the nucleosynthesis and kilonova properties of binary neutron-star (NS) merger models which lead to intermediate remnant lifetimes of ~0.1-1seconds until black-hole (BH) formation and describe all components of material ejected during the dynamical merger phase, NS-remnant evolution, and final viscous disintegration of the BH torus after gravitational collapse. To this end we employ a combination of hydrodynamics, nucleosynthesis, and radiative-transfer tools to achieve a consistent end-to-end modeling of the system and its observables. We adopt a novel version of the Shakura-Sunyaev scheme allowing to vary the approximate turbulent viscosity inside the NS remnant independently of the surrounding disk. We find that asymmetric progenitors lead to shorter remnant lifetimes and enhanced ejecta masses, although the viscosity affects the absolute values of these characteristics. The integrated production of lanthanides and heavier elements in such binary systems is sub-solar, suggesting that the considered scenarios contribute in a sub-dominant fashion to r-process enrichment. One reason is that BH-tori formed after delayed collapse exhibit less neutron-rich conditions than typically found, and often assumed in previous BH-torus models, for early BH formation. The outflows in our models feature strong anisotropy as a result of the lanthanide-poor polar neutrino-driven wind pushing aside lanthanide-rich dynamical ejecta. Considering the complexity of the models, the estimated kilonova light curves show promising agreement with AT2017gfo after times of several days, while the remaining inconsistencies at early times could possibly be overcome in binary configurations with a more dominant neutrino-driven wind relative to the dynamical ejecta.

We discuss the potential of the multi-tracer technique to improve observational constraints of the local primordial non-Gaussianity (PNG) parameter $f_{\rm NL}$ from the galaxy power spectrum. For two galaxy samples $A$ and $B$, we show the constraining power is $\propto |b_1^B b_\phi^A - b_1^A b_\phi^B|$, where $b_1$ and $b_\phi$ are the linear and PNG galaxy bias parameters. This allows for significantly improved constraints compared to the traditional expectation $\propto |b_1^A - b_1^B|$ based on naive universality-like relations where $b_\phi \propto b_1$. Using IllustrisTNG galaxy simulation data, we find that different equal galaxy number splits of the full sample lead to different $|b_1^B b_\phi^A - b_1^A b_\phi^B|$, and thus have different constraining power. Of all of the strategies explored, splitting by $g-r$ color is the most promising, more than doubling the significance of detecting $f_{\rm NL}b_\phi \neq 0$. Importantly, since these are constraints on $f_{\rm NL}b_\phi$ and not $f_{\rm NL}$, they do not require priors on the $b_\phi(b_1)$ relation. For direct constraints on $f_{\rm NL}$, we show that multi-tracer constraints can be significantly more robust than single-tracer to $b_\phi$ misspecifications and uncertainties; this relaxes the precision and accuracy requirements for $b_\phi$ priors. Our results present new opportunities to improve our chances to detect and robustly constrain $f_{\rm NL}$, and strongly motivate galaxy formation simulation campaigns to calibrate the $b_\phi(b_1)$ relation.

In subleading powers of soft-collinear effective theory (SCET), the Lagrangian contains couplings between soft quarks and hard-collinear quarks. Matrix elements of the hard-collinear part of these couplings are radiative jet functions. Although the radiative jet functions are constructed to contain hard-collinear modes, we find that, in order to render the radiative jet functions gauge invariant, it is necessary, in general, to include in the radiative jet functions certain contributions that contain a soft-quark propagator. In gauges other than the Feynman gauge, the soft-quark propagator is canceled by ``gauge terms'' in the gluon propagator, leaving a purely hard-collinear contribution.

We use the BOSS DR12 galaxy power spectrum to constrain compensated isocurvature perturbations (CIP), which are opposite-sign primordial baryon and dark matter perturbations that leave the total matter density unchanged. Long-wavelength CIP $\sigma(\vec{x})$ enter the galaxy density contrast as $\delta_g(\vec{x}) \supset b_\sigma\sigma(\vec{x})$, with $b_\sigma$ the linear CIP galaxy bias parameter. We parameterize the CIP spectra as $P_{\sigma\sigma} = A^2P_{\mathcal{R}\mathcal{R}}$ and $P_{\sigma\mathcal{R}} = \xi\sqrt{P_{\sigma\sigma}P_{\mathcal{R}\mathcal{R}}}$, where $A$ is the CIP amplitude and $\xi$ is the correlation with the curvature perturbations $\mathcal{R}$. We find a significance of detection of $Ab_\sigma \neq 0$ of $1.8\sigma$ for correlated ($\xi = 1$) and $3.7\sigma$ for uncorrelated ($\xi = 0$) CIP. Large-scale data systematics have a bigger impact for uncorrelated CIP, which may explain the large significance of detection. The constraints on $A$ depend on the assumed priors for the $b_\sigma$ parameter, which we estimate using separate universe simulations. Assuming $b_\sigma$ values representative of all halos we find $\sigma_A = 145$ for correlated CIP and $\sigma_{|A|} = 475$ for uncorrelated CIP. Our strongest uncorrelated CIP constraint is for $b_\sigma$ representative of the $33\%$ most concentrated halos, $\sigma_{|A|} = 197$, which is better than the current CMB bounds $|A| \lesssim 360$. We also discuss the impact of the local primordial non-Gaussianity parameter $f_{\rm NL}$ in CIP constraints. Our results demonstrate the power of galaxy data to place tight constraints on CIP, and motivate works to understand better the impact of data systematics, as well as to determine theory priors for $b_\sigma$.

Context: Several observations of the local Universe (LU) point towards the existence of very prominent structures. The presence of massive galaxy clusters and local super clusters on the one hand, but also large local voids and under-densities on the other hand. However, it is highly non trivial to connect such different observational selected tracers to the underlying dark matter (DM) distribution. Methods (abridged): We used a 500 Mpc/h large constrained simulation of the LU with initial conditions based on peculiar velocities derived from the CosmicFlows-2 catalogue and follow galaxy formation physics directly in the hydro-dynamical simulations to base the comparison on stellar masses of galaxies or X-ray luminosity of clusters. We also used the 2668 Mpc/h large cosmological box from the Magneticum simulations to evaluate the frequency of finding such anomalies in random patches within simulations. Results: We demonstrate that haloes and galaxies in our constrained simulation trace the local DM density field very differently. Thereby, this simulation reproduces the observed 50% under-density of galaxy clusters and groups within the sphere of ~100 Mpc when applying the same mass or X-ray luminosity limit used in the observed cluster sample (CLASSIX), which is consistent with a ~1.5$\sigma$ feature. At the same time, the simulation reproduces the observed over-density of massive galaxy clusters within the same sphere, which on its own also corresponds to a ~1.5$\sigma$ feature. Interestingly, we find that only 44 out of 15635 random realizations (i.e. 0.28%) are matching both anomalies, making the LU to be a ~3$\sigma$ environment. We finally compared a mock galaxy catalogue with the observed distribution of galaxies in the LU, finding also a match to the observed factor of two over-density at ~16 Mpc as well as the observed 15% under-density at ~40 Mpc distance.

Free-floating planets (FFPs) can result from dynamical scattering processes happening in the first few million years of a planetary system's life. Several models predict the possibility, for these isolated planetary-mass objects, to retain exomoons after their ejection. The tidal heating mechanism and the presence of an atmosphere with a relatively high optical thickness may support the formation and maintenance of oceans of liquid water on the surface of these satellites. In order to study the timescales over which liquid water can be maintained, we perform dynamical simulations of the ejection process and infer the resulting statistics of the population of surviving exomoons around free-floating planets. The subsequent tidal evolution of the moons' orbital parameters is a pivotal step to determine when the orbits will circularize, with a consequential decay of the tidal heating. We find that close-in ($a \lesssim 25 $R$_{\rm J}$) Earth-mass moons with CO$_2$-dominated atmospheres could retain liquid water on their surfaces for long timescales, depending on the mass of the atmospheric envelope and the surface pressure assumed. Massive atmospheres are needed to trap the heat produced by tidal friction that makes these moons habitable. For Earth-like pressure conditions ($p_0$ = 1 bar), satellites could sustain liquid water on their surfaces up to 52 Myr. For higher surface pressures (10 and 100 bar), moons could be habitable up to 276 Myr and 1.6 Gyr, respectively. Close-in satellites experience habitable conditions for long timescales, and during the ejection of the FFP remain bound with the escaping planet, being less affected by the close encounter.

Comparing Galactic chemical evolution models to the observed elemental abundances in the Milky Way, we show that neutron star mergers can be a leading r-process site only if at low metallicities such mergers have very short delay times and significant ejecta masses that are facilitated by the masses of the compact objects. Namely, black hole-neutron star mergers, depending on the black hole spins, can play an important role in the early chemical enrichment of the Milky Way. We also show that none of the binary population synthesis models used in this Letter, i.e., COMPAS, StarTrack, Brussels, ComBinE, and BPASS, can currently reproduce the elemental abundance observations. The predictions are problematic not only for neutron star mergers, but also for Type Ia supernovae, which may point to shortcomings in binary evolution models.

We present new VLA observations, between 6.8mm and 66mm, of the edge-on Class~I disk IRAS04302+2247. Observations at 6.8mm and 9.2mm lead to the detection of thermal emission from the disk, while shallow observations at the other wavelengths are used to correct for emission from other processes. The disk radial brightness profile transitions from broadly extended in previous ALMA 0.9mm and 2.1mm observations to much more centrally brightened at 6.8mm and 9.2mm, which can be explained by optical depth effects. The radiative transfer modeling of the 0.9mm, 2.1mm, and 9.2mm data suggests that the grains are smaller than 1cm in the outer regions of the disk and allows us to obtain the first lower limit for the scale height of grains emitting at millimeter wavelengths in a protoplanetary disk. We find that the millimeter dust scale height is between 1au and 6au at a radius 100au from the central star, while the gas scale height is estimated to be about 7au, indicating a modest level of settling. The estimated dust height is intermediate between less evolved Class 0 sources, that are found to be vertically thick, and more evolved Class II sources, which show a significant level of settling. This suggests that we are witnessing an intermediate stage of dust settling.

Context. While it is generally assumed that Class II sources evolve largely in isolation from their environment, many still lie close to molecular clouds and may continue to interact with them. This may result in late accretion of material onto the disk that can significantly influence disk processes and planet formation.
Aims: In order to systematically study late infall of gas onto disks, we identify candidate Class II sources in close vicinity to a reflection nebula (RN) that may be undergoing this process.
Methods: First we targeted Class II sources with known kilo-au scale gas structures - possibly due to late infall of material - and we searched for RNe in their vicinity in optical and near-infrared images. Second, we compiled a catalogue of Class II sources associated with RNe and looked for the large-scale CO structures in archival ALMA data. Using the catalogues of protostars and RNe, we also estimated the probability of Class II sources interacting with surrounding material.
Results: All of the sources with large-scale gas structures also exhibit some reflection nebulosity in their vicinity. Similarly, at least five Class II objects associated with a prominent RNe, and for which adequate ALMA observations are available, were found to have spirals or stream-like structures which may be due to late infall. We report the first detection of these structures around S CrA.
Conclusions: Our results suggest that a non-negligible fraction of Class II disks in nearby star-forming regions may be associated with RNe and could therefore be undergoing late accretion of gas. Surveys of RNe and kilo-au scale gas structures around Class II sources will allow us to better understand the frequency and impact of late-infall phenomena.

This work presents the results from extending the long-term monitoring program of stellar motions within the Galactic Center to include stars with separations of 2-7 arcsec from the compact radio source, Sgr A*. In comparison to the well studied inner 2 arcsec, a longer baseline in time is required to study these stars. With 17 years of data, a sufficient number of positions along the orbits of these outer stars can now be measured. This was achieved by designing a source finder to track the positions of ∼2000 stars in NACO/VLT adaptive-optics-assisted images of the Galactic Center from 2002 to 2019. Of the studied stars, 54 exhibit significant accelerations toward Sgr A*, most of which have separations of between 2 and 3 arcsec from the black hole. A further 20 of these stars have measurable radial velocities from SINFONI/VLT stellar spectra, which allows for the calculation of the orbital elements for these stars, thus increasing the number of known orbits in the Galactic Center by ∼40%. With orbits, we can consider which structural features within the Galactic Center nuclear star cluster these stars belong to. Most of the stars have orbital solutions that are consistent with the known clockwise rotating disk feature. Further, by employing Monte Carlo sampling for stars without radial velocity measurements, we show that many stars have a subset of possible orbits that are consistent with one of the known disk features within the Galactic Center.

LiteBIRD is a planned JAXA-led CMB B-mode satellite experiment aiming for launch in the late 2020s, with a primary goal of detecting the imprint of primordial inflationary gravitational waves. Its current baseline focal-plane configuration includes 15 frequency bands between 40 and 402 GHz, fulfilling the mission requirements to detect the amplitude of gravitational waves with the total uncertainty on the tensor-to-scalar ratio, $\delta r$, down to $\delta r<0.001$. A key aspect of this performance is accurate astrophysical component separation, and the ability to remove polarized thermal dust emission is particularly important. In this paper we note that the CMB frequency spectrum falls off nearly exponentially above 300 GHz relative to the thermal dust SED, and a relatively minor high frequency extension can therefore result in even lower uncertainties and better model reconstructions. Specifically, we compare the baseline design with five extended configurations, while varying the underlying dust modeling, in each of which the HFT (High-Frequency Telescope) frequency range is shifted logarithmically towards higher frequencies, with an upper cutoff ranging between 400 and 600 GHz. In each case, we measure the tensor-to-scalar ratio $r$ uncertainty and bias using both parametric and minimum-variance component-separation algorithms. When the thermal dust sky model includes a spatially varying spectral index and temperature, we find that the statistical uncertainty on $r$ after foreground cleaning may be reduced by as much as 30--50 % by extending the upper limit of the frequency range from 400 to 600 GHz, with most of the improvement already gained at 500 GHz. We also note that a broader frequency range leads to better ability to discriminate between models through higher $\chi^2$ sensitivity. (abridged)

Understanding the sequence-dependent DNA damage formation requires probing a complete pool of sequences over a wide dose range of the damage-causing exposure. We used high throughput sequencing to simultaneously obtain the dose dependence and quantum yields for oligonucleotide damages for all possible 4096 DNA sequences with hexamer length. We exposed the DNA to ultraviolet radiation at 266 nm and doses of up to 500 absorbed photons per base. At the dimer level, our results confirm existing literature values of photodamage, whereas we now quantified the susceptibility of sequence motifs to UV irradiation up to previously inaccessible polymer lengths. This revealed the protective effect of the sequence context in preventing the formation of UV-lesions. For example, the rate to form dipyrimidine lesions is strongly reduced by nearby guanine bases. Our results provide a complete picture of the sensitivity of oligonucleotides to UV irradiation and allow us to predict their abundance in high-UV environments.

We analyse $Z^\prime$ contributions to FCNC processes at the one-loop level. In analogy to the CKM matrix we introduce two $3\times3$ unitary matrices $\hat\Delta_d(Z^\prime)$ and $\hat\Delta_u(Z^\prime)$ which are also hermitian. They govern the flavour interactions mediated by $Z^\prime$ between down-quarks and up-quarks, respectively, with $\hat\Delta_d(Z^\prime)=\hat\Delta_u(Z^\prime)\equiv \hat\Delta_L(Z^\prime)$ for left-handed currents due to the unbroken $\text{SU(2)}_L$ gauge symmetry. This assures the suppression of these contributions to all $Z^\prime$ mediated FCNC processes at the one-loop level. As, in contrast to the GIM mechanism, one-loop $Z^\prime$ contributions to flavour observables in $K$ and $B_{s,d}$ systems are governed by down-quark masses, they are ${\cal O}(m^2_b/M^2_{Z^\prime})$ and negligible. With the ${\cal O}(m^2_t/M^2_{Z^\prime})$ suppression they are likely negligible also in the $D$ system. We present an explicit parametrization of $\hat\Delta_L(Z^\prime)$ in terms of two mixing angles and two complex phases that distinguishes it profoundly from the CKM matrix. This framework can be generalized to purely leptonic decays with matrices analogous to the PMNS matrix but profoundly different from it. Interestingly, the breakdown of flavour universality between the first two generations and the third one, both for quark and lepton couplings to $Z^\prime$, is identified as a consequence of $\hat\Delta_L(Z^\prime)$ being hermitian. The importance of the unitarity for both $\hat\Delta_L(Z^\prime)$ and the CKM matrix in the light of the Cabibbo anomaly is emphasized.

Non-thermal emission from relativistic cosmic ray (CR) electrons gives insight into the strength and morphology of intra-cluster magnetic fields, as well as providing powerful tracers of structure formation shocks. Emission caused by CR protons on the other hand still challenges current observations and is therefore testing models of proton acceleration at intra-cluster shocks. Large-scale simulations including the effects of CRs have been difficult to achieve and have been mainly reduced to simulating an overall energy budget, or tracing CR populations in post-processing of simulation output and has often been done for either protons or electrons. We introduce CRESCENDO: Cosmic Ray Evolution with SpeCtral Electrons aND prOtons, an efficient on-the-fly Fokker-Planck solver to evolve distributions of CR protons and electrons within every resolution element of our simulation. The solver accounts for CR (re-)acceleration at intra-cluster shocks, based on results of recent particle-in-cell simulations, adiabatic changes, and radiative losses of electrons. We show its performance in test cases as well as idealized galaxy cluster (GC) simulations. We apply the model to an idealized GC merger following best-fitting parameters for CIZA J2242.4 + 5301-1 and study CR injection, radio relic morphology, spectral steepening, and synchrotron emission.

The Euclid space telescope will survey a large dataset of cosmic voids traced by dense samples of galaxies. In this work we estimate its expected performance when exploiting angular photometric void clustering, galaxy weak lensing, and their cross-correlation. To this aim, we implemented a Fisher matrix approach tailored for voids from the Euclid photometric dataset and we present the first forecasts on cosmological parameters that include the void-lensing correlation. We examined two different probe settings, pessimistic and optimistic, both for void clustering and galaxy lensing. We carried out forecast analyses in four model cosmologies, accounting for a varying total neutrino mass, M_ν, and a dynamical dark energy (DE) equation of state, w(z), described by the popular Chevallier-Polarski-Linder parametrization. We find that void clustering constraints on h and Ω_b are competitive with galaxy lensing alone, while errors on n_s decrease thanks to the orthogonality of the two probes in the 2D-projected parameter space. We also note that, as a whole, with respect to assuming the two probes as independent, the inclusion of the void-lensing cross-correlation signal improves parameter constraints by 10 − 15%, and enhances the joint void clustering and galaxy lensing figure of merit (FoM) by 10% and 25%, in the pessimistic and optimistic scenarios, respectively. Finally, when further combining with the spectroscopic galaxy clustering, assumed as an independent probe, we find that, in the most competitive case, the FoM increases by a factor of 4 with respect to the combination of weak lensing and spectroscopic galaxy clustering taken as independent probes. The forecasts presented in this work show that photometric void clustering and its cross-correlation with galaxy lensing deserve to be exploited in the data analysis of the Euclid galaxy survey and promise to improve its constraining power, especially on h, Ω_b, the neutrino mass, and the DE evolution.

This paper is published on behalf of the Euclid Consortium.

We investigate the accuracy and precision of triaxial dynamical orbit models by fitting 2D mock observations of a realistic N-body merger simulation resembling a massive early-type galaxy with a supermassive black hole. We show that we can reproduce the triaxial N-body merger remnant's correct black hole mass, stellar mass-to-light ratio and total enclosed mass (inside the half-light radius) for several different tested orientations with an unprecedented accuracy of 5-10 per cent. Our dynamical models use the entire non-parametric line-of-sight velocity distribution (LOSVD) rather than parametric LOSVDs or velocity moments as constraints. Our results strongly suggest that state-of-the-art integral-field projected kinematic data contain only minor degeneracies with respect to the mass and anisotropy recovery. Moreover, this also demonstrates the strength of the Schwarzschild method in general. We achieve the proven high recovery accuracy and precision with our newly developed modelling machinery by combining several advancements: (i) our new semiparametric deprojection code probes degeneracies and allows us to constrain the viewing angles of a triaxial galaxy; (ii) our new orbit modelling code SMART uses a 5-dim orbital starting space to representatively sample in particular near-Keplerian orbits in galaxy centres; (iii) we use a generalized information criterion AIC_p to optimize the smoothing and to compare different mass models to avoid biases that occur in χ²-based models with varying model flexibilities.

In the electroweak sector of the Standard Model, CP violation arises through a very particular interplay between the three quark generations, as described by the Cabibbo--Kobayashi--Maskawa (CKM) mechanism and the single Jarlskog invariant $J_4$. Once generalized to the Standard Model Effective Field Theory (SMEFT), this peculiar pattern gets modified by higher-dimensional operators, whose associated Wilson coefficients are usually split into CP-even and odd parts. However, CP violation at dimension four, i.e., at the lowest order in the EFT expansion, blurs this distinction: any Wilson coefficient can interfere with $J_4$ and mediate CP violation. In this paper, we study such interferences at first order in the SMEFT expansion, $\order{1/\Lambda^2}$, and we capture their associated parameter space via a set of 1551 linear CP-odd flavor invariants. This construction describes both new, genuinely CP-violating quantities as well as the interference between $J_4$ and CP-conserving ones. We call this latter possibility \textit{opportunistic CP violation}. Relying on an appropriate extension of the matrix rank to Taylor expansions, which we dub \emph{Taylor rank}, we define a procedure to organize the invariants in terms of their magnitude, so as to retain only the relevant ones at a given precision. We explore how this characterization changes when different assumptions are made on the flavor structure of the SMEFT coefficients. Interestingly, some of the CP-odd invariants turn out to be less suppressed than $J_4$, even when they capture opportunistic CPV, demonstrating that CP-violation in the SM, at dimension 4, is \textit{accidentally small.}

In the past few years, the Event Horizon Telescope (EHT) has provided the first-ever event horizon-scale images of the supermassive black holes (BHs) M87* and Sagittarius A* (Sgr A*). The next-generation EHT project is an extension of the EHT array that promises larger angular resolution and higher sensitivity to the dim, extended flux around the central ring-like structure, possibly connecting the accretion flow and the jet. The ngEHT Analysis Challenges aim to understand the science extractability from synthetic images and movies to inform the ngEHT array design and analysis algorithm development. In this work, we compare the accretion flow structure and dynamics in numerical fluid simulations that specifically target M87* and Sgr A*, and were used to construct the source models in the challenge set. We consider (1) a steady-state axisymmetric radiatively inefficient accretion flow model with a time-dependent shearing hotspot, (2) two time-dependent single fluid general relativistic magnetohydrodynamic (GRMHD) simulations from the H-AMR code, (3) a two-temperature GRMHD simulation from the BHAC code, and (4) a two-temperature radiative GRMHD simulation from the KORAL code. We find that the different models exhibit remarkably similar temporal and spatial properties, except for the electron temperature, since radiative losses substantially cool down electrons near the BH and the jet sheath, signaling the importance of radiative cooling even for slowly accreting BHs such as M87*. We restrict ourselves to standard torus accretion flows, and leave larger explorations of alternate accretion models to future work.

We study stellar population and structural properties of massive log (M_⋆/M_⊙) > 11 galaxies at z ≈ 2.7 in the Magneticum and IllustrisTNG hydrodynamical simulations and GAEA semi-analytic model. We find stellar mass functions broadly consistent with observations, with no scarcity of massive, quiescent galaxies at z ≈ 2.7, but with a higher quiescent galaxy fraction at high masses in IllustrisTNG. Average ages of simulated quiescent galaxies are between ≈0.8 and ${1.0\, \textrm {Gyr}}$, older by a factor ≈2 than observed in spectroscopically confirmed quiescent galaxies at similar redshift. Besides being potentially indicative of limitations of simulations in reproducing observed star formation histories, this discrepancy may also reflect limitations in the estimation of observed ages. We investigate the purity of simulated UVJ rest-frame colour-selected massive quiescent samples with photometric uncertainties typical of deep surveys (e.g. COSMOS). We find evidence for significant contamination (up to ${60\, \rm {per\, cent}}$) by dusty star-forming galaxies in the UVJ region that is typically populated by older quiescent sources. Furthermore, the completeness of UVJ-selected quiescent samples at this redshift may be reduced by $\approx {30\, \rm {per\, cent}}$ due to a high fraction of young quiescent galaxies not entering the UVJ quiescent region. Massive, quiescent galaxies in simulations have on average lower angular momenta and higher projected axis ratios and concentrations than star-forming counterparts. Average sizes of simulated quiescent galaxies are broadly consistent with observations within the uncertainties. The average size ratio of quiescent and star-forming galaxies in the probed mass range is formally consistent with observations, although this result is partly affected by poor statistics.

We discuss the gravitational wave spectrum produced by first-order phase transitions seeded by domain wall networks. This setup is important for many two-step phase transitions as seen for example in the singlet extension of the standard model. Whenever the correlation length of the domain wall network is larger than the typical bubble size, this setup leads to a gravitational wave signal that is shifted to lower frequencies and with an enhanced amplitude compared to homogeneous phase transitions without domain walls. We discuss our results in light of the recent PTA hints for gravitational waves.

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, for arbitrary fiducial background cosmologies, and for a large parameter space that covers the shape parameters ω_c, ω_b, and n_s, as well as the evolution parameters h, A_s, Ω_K, w₀, 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 two orders of magnitude, validation tests reveal an accuracy of $0.1\, {{\ \rm per\ cent}}$ for the monopole and quadrupole ($0.3\, {{\ \rm per\ cent}}$ 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 and Gaussian covariance matrices.

$Z^\prime$ models belong to the ones that can most easily explain the anomalies in $b\to s \mu^+\mu^-$ transitions. However, such an explanation by a single $Z^\prime$ gauge boson, as done in the literature, is severly constrained by the $B^0_s-\bar B_s^0$ mixing. Also the recent finding, hat the mass differences $\Delta M_s$, $\Delta M_d$, the CP-violating parameter $\varepsilon_K$, and the mixing induced CP-asymmetries $S_{\psi K_S}$ and $S_{\psi \phi}$ can be simultaneously well described within the SM without new physics (NP) contributions, is a challenge for $Z^\prime$ models with a single $Z^\prime$ contributing at tree-level to quark mixing. We point out that including a second $Z^\prime$ in the model allows to eliminate simultaneously tree-level contributions to the five $\Delta F=2$ observables used in the determination of the CKM parameters while leaving the room for NP in $\Delta M_K$ and $\Delta M_D$. The latter one can be removed at the price of infecting $\Delta M_s$ or $\Delta M_d$ by NP which is presently disfavoured. This pattern is transparently seen using the new mixing matrix for $Z^\prime$ interactions with quarks. This strategy allows significant tree-level contributions to $K$, $B_s$ and $B_d$ decays thereby allowing to explain the existing anomalies in $b\to s\mu^+\mu^-$ transitions and the anticipated anomaly in the ratio $\varepsilon'/\varepsilon$ much easier than in $Z^\prime$-Single scenarios. The proposed $Z^\prime$-Tandem mechanism bears some similarities to the GIM mechanism for the suppression of the FCNCs in the SM with the role of the charm quark played here by the second $Z^\prime$. However, it differs from the latter profoundly in that only NP contributions to quark mixing are eliminated at tree-level. We discuss briefly the implied flavour patterns in $K$ and $B$ decay observables in this NP scenario.

We compute bottom mass ($m_b$) corrections to the transverse momentum ($q_T$) spectrum of Higgs bosons produced by gluon fusion in the regime $q_T \sim m_b \ll m_H$ at leading power in $m_b/m_H$ and $q_T/m_H$, where the gluons couple to the Higgs via a top loop. To this end we calculate the quark mass dependence of the transverse momentum dependent gluon beam functions (aka gluon TMDPDFs) at two loops in the framework of SCET. These functions represent the collinear matrix elements in the factorized gluon-fusion cross section for small $q_T$. We discuss in detail technical subtleties regarding rapidity regulators and zero-bin subtractions in the calculation of the virtual corrections present for massive quarks. Combined with the known soft function for $m_b \neq 0$ our results allow to determine the resummed Higgs $q_T$ distribution at NNLL$^\prime$ (and eventually N$^3$LL) with full dependence on $m_b/q_T$. At fixed order the new corrections to the massless approximation lead to percent-level effects in the peak region of the Higgs $q_T$ spectrum. Upon resummation they may thus be relevant for state-of-the-art precision predictions for the LHC.