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

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

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

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

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

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

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

The essence of the \textit{memory burden} effect is that a load of information carried by a system stabilizes it. This universal effect is especially prominent in systems with a high capacity of information storage, such as black holes and other objects with maximal microstate degeneracy, the entities universally referred to as \textit{saturons}. The phenomenon has several implications. The memory burden effect suppresses a further decay of a black hole, the latest, after it has emitted about half of its initial mass. As a consequence, the light primordial black holes (PBHs), that previously were assumed to be fully evaporated, are expected to be present as viable dark matter candidates. In the present paper, we deepen the understanding of the memory burden effect. We first identify various memory burden regimes in generic Hamiltonian systems and then establish a precise correspondence in solitons and in black holes. We make transparent, at a microscopic level, the fundamental differences between the stabilization by a quantum memory burden versus the stabilization by a long-range classical hair due to a spin or an electric charge. We identify certain new features of potential observational interest, such as the model-independent spread of the stabilized masses of initially degenerate PBHs.

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

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

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

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

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

We present new VLT/MUSE observations of the Hubble Frontier Field (HFF) galaxy cluster MACS J1149.5+2223, lensing the well-known supernova "Refsdal" into multiple images, which enabled the first cosmological applications with a strongly lensed supernova. Thanks to these data, targeting a northern region of the cluster and thus complementing our previous MUSE program on the cluster core, we release a new catalog containing 162 secure spectroscopic redshifts. We confirm 22 cluster members, which were previously only photometrically selected, and detect ten additional ones, resulting in a total of 308 secure members, of which 63% are spectroscopically confirmed. We further identify 17 new spectroscopic multiple images belonging to 6 different background sources. By exploiting MUSE data, in combination with the deep HFF images, we develop an improved total mass model of MACS J1149.5+2223. This model includes 308 total mass components for the member galaxies and requires four additional mass profiles, one of which is associated with a cluster galaxy overdensity identified in the North, representing the DM mass distribution on larger scales. The values of the resulting 34 free parameters are optimized based on the observed positions of 106 multiple images from 34 different families, that cover the redshift range between 1.240 and 5.983. Our final model has a multiple image position rms value of 0.39", which is well in agreement with that of other cluster lens models. With this refined mass model, we pave the way towards even better strong-lensing analyses that will exploit the deep and high resolution observations with HST and JWST on a pixel level in the region of the supernova Refsdal host. This will increase the number of observables by around two orders of magnitudes, thus offering us the opportunity of carrying out more precise and accurate cosmographic measurements.

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

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

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

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

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

Context. Merging compact objects such as binary black holes provide a promising probe for the physics of dark matter (DM). The gravitational waves emitted during inspiral potentially allow to detect DM spikes around black holes. This is because the dynamical friction force experienced by the inspiraling black hole alters the orbital period and thus the gravitational wave signal. Aims. The dynamical friction arising from DM can potentially differ from the collisionless case when DM is subject to self-interactions. This paper aims to understand how self-interactions impact dynamical friction. Methods. To study the dynamical friction force, we use idealized N-body simulations, where we include self-interacting dark matter. Results. We find that the dynamical friction force for inspiraling black holes would be typically enhanced by DM self-interactions compared to a collisionless medium (ignoring differences in the DM density). At lower velocities below the sound speed, we find that the dynamical friction force can be reduced by the presence of self-interactions. Conclusions. DM self-interactions have a significant effect on the dynamical friction for black hole mergers. Assuming the Chandrasekhar formula may underpredict the deceleration due to dynamical friction.

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

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

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

The analysis of several spectroscopic surveys indicates the presence of a bimodality between the disc stars in the abundance ratio space of [${\alpha}$/Fe] versus [Fe/H]. The two stellar groups are commonly referred to as the high-${\alpha}$ and low-${\alpha}$ sequences. Some models capable of reproducing such a bimodality, invoke the presence of a hiatus in the star formation history in our Galaxy, whereas other models explain the two sequences by means of stellar migration. Our aim is to show that the existence of the gap in the star formation rate between high-$\alpha$ and low-$\alpha$ is evident in the stars of APOGEE DR17, if one plots [Fe/$\alpha$] versus [$\alpha$/H], thus confirming previous suggestions by Gratton et al. (1996) and Fuhrmann (1998). Then we try to interpret the data by means of detailed chemical models. We compare the APOGEE DR17 red giant stars with the predictions of a detailed chemical evolution model based on the two-infall paradigm, taking also into account possible accretion of dwarf satellites. The APOGEE DR17 abundance ratios [Fe/$\alpha$] versus [$\alpha$/H] exhibit a sharp increase of [Fe/$\alpha$] at a nearly constant [$\alpha$/H] (where $\alpha$ elements considered are Mg, Si, O) during the transition between the two disc phases. This observation strongly supports the hypothesis that a hiatus in star formation occurred during this evolutionary phase. Notably, the most pronounced growth in the [Fe/$\alpha$] versus [$\alpha$/H] relation is observed for oxygen, as this element is exclusively synthesised in core-collapse supernovae. A chemical model predicting a stop in the star formation of a duration of roughly 3.5 Gyr, and where the high-$\alpha$ disc starts forming from pre-enriched gas by a previous encounter with a dwarf galaxy can well explain the observations.

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

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

The metallicity distribution function (MDF) of the Galactic bulge features a multi-peak shape, with a metal-poor peak at [Fe/H]=-0.3 dex and a metal-rich peak at [Fe/H]=+0.3 dex. This bimodality is also seen in [alpha/Fe] versus [Fe/H] ratios, indicating different stellar populations in the bulge. We aim to replicate the observed MDF by proposing a scenario where the metal-poor bulge stars formed in situ during an intense star formation burst, while the metal-rich stars formed during a second burst and/or were accreted from the inner Galactic disk due to a growing bar. We used a chemical evolution model that tracks various chemical species with detailed nucleosynthesis, focusing on Fe production from both Type Ia supernovae and massive stars, including rotating massive stars with varying velocities. Our model also accounts for gas infall, outflow, and the effect of stellar migration. Results are compared to 13,000 stars from the SDSS/APOGEE survey within 3.5 kpc of the Galactic center. Our model successfully reproduces the double-peak shape of the bulge MDF and the alpha-element abundance trends relative to Fe by assuming (i) a multi-burst star formation history with a 250 Myr quenching of the first burst and (ii) stellar migration from the inner disk due to a growing bar. We estimate that about 40% of the bulge-bar's stellar mass originates from the inner disk. Nucleosynthesis models that assume either no rotation for massive stars or a rotational velocity distribution favoring slow rotation at high metallicities best match the observed MDF and [alpha/Fe] and [Ce/Fe] versus [Fe/H] abundance patterns.

Large mm surveys of star forming regions enable the study of entire populations of planet-forming disks and reveal correlations between their observable properties. Population studies of disks have shown that the correlation between disk size and millimeter flux could be explained either through disks with strong substructure, or alternatively by the effects of radial inward drift of growing dust particles. This study aims to constrain the parameters and initial conditions of planet-forming disks and address the question of the need for the presence of substructures in disks and, if needed, their predicted characteristics, based on the large samples of disk sizes, millimeter fluxes, and spectral indices available. We performed a population synthesis of the continuum emission of disks, exploiting a two-population model (two-pop-py), considering the influence of viscous evolution, dust growth, fragmentation, and transport varying the initial conditions of the disk and substructure to find the best match to the observed distributions. We show that the observed distributions of spectral indices, sizes, and luminosities together can be best reproduced by disks with significant substructure, namely a perturbation strong enough to be able to trap particles, and that is formed early in the evolution of the disk, that is within 0.4Myr. Agreement is reached by relatively high initial disk masses ($10^{-2.3}M_{\star}\leqslant M_{disk}\leqslant10^{-0.5}M_{\star}$) and moderate levels of turbulence ($10^{-3.5}\leqslant\alpha\leqslant 10^{-2.5}$). Other disk parameters play a weaker role. Only opacities with high absorption efficiency can reproduce the observed spectral indices. Our results extend to the whole population that substructure is likely ubiquitous, so far assessed only in individual disks and implies that most "smooth" disks hide unresolved substructure.

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

We introduce a block encoding method for mapping discrete subgroups to qubits on a quantum computer. This method is applicable to general discrete groups, including crystal-like subgroups such as $\mathbb{BI}$ of $SU(2)$ and $\mathbb{V}$ of $SU(3)$. We detail the construction of primitive gates -- the inversion gate, the group multiplication gate, the trace gate, and the group Fourier gate -- utilizing this encoding method for $\mathbb{BT}$ and for the first time $\mathbb{BI}$ group. We also provide resource estimations to extract the gluon viscosity. The inversion gates for $\mathbb{BT}$ and $\mathbb{BI}$ are benchmarked on the $\texttt{Baiwang}$ quantum computer with estimated fidelities of $40^{+5}_{-4}\%$ and $4^{+5}_{-3}\%$ respectively.

Life continuously transduces energy to perform critical functions using energy stored in reactive molecules like ATP or NADH. ATP dynamically phosphorylates active sites on proteins and thereby regulates their function. Inspired by such machinery, regulating supramolecular functions using energy stored in reactive molecules has gained traction. Enzyme-free, synthetic systems that use dynamic phosphorylation to regulate supramolecular processes have not yet been reported, to our knowledge. Here, we show an enzyme-free reaction cycle that consumes the phosphorylating agent monoamidophosphate by transiently phosphorylating histidine and histidine-containing peptides. The phosphorylated species are labile and deactivate through hydrolysis. The cycle exhibits versatility and tunability, allowing for the dynamic phosphorylation of multiple precursors with a tunable half-life. Notably, we show the resulting phosphorylated products can regulate the peptide's phase separation, leading to active droplets that require the continuous conversion of fuel to sustain. The reaction cycle will be valuable as a model for biological phosphorylation but can also offer insights into protocell formation.

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

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

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

The main goal of the CRESST-III experiment is the direct detection of dark matter particles via their scattering off target nuclei in cryogenic detectors. In this work we present the results of a Silicon-On-Sapphire (SOS) detector with a mass of 0.6$\,$g and an energy threshold of (6.7$\, \pm \,$0.2)$\,$eV with a baseline energy resolution of (1.0$\, \pm \,$0.2)$\,$eV. This allowed for a calibration via the detection of single luminescence photons in the eV-range, which could be observed in CRESST for the first time. We present new exclusion limits on the spin-independent and spin-dependent dark matter-nucleon cross section that extend to dark matter particle masses of less than 100$\,$MeV/c$^{2}$.

One of the most promising approaches for the next generation of neutrino experiments is the realization of large hybrid Cherenkov/scintillation detectors made possible by recent innovations in photodetection technology and liquid scintillator chemistry. The development of a potentially suitable future detector liquid with particularly slow light emission is discussed in the present publication. This cocktail is compared with respect to its fundamental characteristics (scintillation efficiency, transparency, and time profile of light emission) with liquid scintillators currently used in large-scale neutrino detectors. In addition, the optimization of the admixture of wavelength shifters for a scintillator with particularly high light emission is presented. Furthermore, the pulse-shape discrimination capabilities of the novel medium was studied using a pulsed particle accelerator driven neutron source. Beyond that, purification methods based on column chromatography and fractional vacuum distillation for the co-solvent DIN (Diisopropylnaphthalene) are discussed.

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

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

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

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

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

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

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

We investigate the possibility of disentangling different new physics contributions to the rare meson decays $K\rightarrow\pi+\displaystyle{\not}E$ and $B\rightarrow K(K^*)+\displaystyle{\not}E$ through kinematic distributions in the missing energy $\displaystyle{\not}E$. We employ dimension-$6$ operators within the Low-Energy Effective Field Theory (LEFT), identifying the invisible part of the final state as either active or sterile neutrinos. Special emphasis is given to lepton-number violating (LNV) operators with scalar and tensor currents. We show analytically that contributions from scalar, vector, and tensor quark currents can be uniquely determined from experimental data of kinematic distributions. In addition, we present new correlations of branching ratios for $K$ and $B$-decays involving scalar and tensor currents. As there could a priori also be new invisible particles in the final states, we include dark-sector operators giving rise to two dark scalars, fermions, or vectors in the final state. In this context, we present new calculations of the inclusive decay rate $B\rightarrow X_s+\displaystyle{\not}E$ for dark operators. We show that careful measurements of kinematic distributions make it theoretically possible to disentangle the contribution from LEFT operators from most of the dark-sector operators, even when multiple operators are contributing. We revisit sum rules for vector currents in LEFT and show that the latter are also satisfied in some new dark-physics scenarios that could mimic LEFT. Finally, we point out that an excess in rare meson decays consistent with a LNV hypothesis would point towards highly flavor non-democratic physics in the UV, and could put high-scale leptogenesis under tension.

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

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

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

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

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

The Gaia mission has delivered hundreds of thousands of variable star light curves in multiple wavelengths. Recent work demonstrates that these light curves can be used to identify (non-)radial pulsations in the OBAF-type stars, despite the irregular cadence and low light curve precision of order a few mmag. With the considerably more precise TESS photometry, we revisit these candidate pulsators to conclusively ascertain the nature of their variability. We seek to re-classify the Gaia light curves with the first two years of TESS photometry for a sample of 58,970 p- and g- mode pulsators, encompassing gamma Dor, delta Scuti, SPB, and beta Cep variables. We also supply four new catalogues containing the confirmed pulsators, along with their dominant and secondary pulsation frequencies, the number of independent mode frequencies, and a ranking according to their usefulness for future asteroseismic ensemble analysis. We find that the Gaia photometry is exceptionally accurate for detecting the dominant and secondary frequencies, reaching approximately 80% accuracy in frequency for p- and g-mode pulsators. The majority of Gaia classifications are consistent with the classifications from the TESS data, illustrating the power of the low-cadence Gaia photometry for pulsation studies. We find that the sample of g-mode pulsators forms a continuous group of variable stars along the main sequence across B, A, and F spectral types, implying that the mode excitation mechanisms for all these pulsators need to be updated with improved physics. Finally, we provide a rank-ordered table of pulsators according to their asteroseismic potential for follow-up studies. Our catalogue offers a major increase in the number of confirmed gravity-mode pulsators with an identified dominant mode suitable for follow-up TESS ensemble asteroseismology of such stars.

We propose a simple fit function, $L_{\nu_i}(t) = C\, t^{-\alpha}\, e^{-(t/\tau)^{n}}$, to parametrize the luminosities of neutrinos and antineutrinos of all flavors during the protoneutron star (PNS) cooling phase at post-bounce times $t \gtrsim 1$ s. This fit is based on results from a set of neutrino-hydrodynamics simulations of core-collapse supernovae in spherical symmetry. The simulations were performed with an energy-dependent transport for six neutrino species and took into account the effects of convection and muons in the dense and hot PNS interior. We provide values of the fit parameters $C$, $\alpha$, $\tau$, and $n$ for different neutron star masses and equations of state as well as correlations between these fit parameters. Our functional description is useful for analytic supernova modeling, for characterizing the neutrino light curves in large underground neutrino detectors, and as a tool to extract information from measured signals on the mass and equation of state of the PNS and on secondary signal components on top of the PNS's neutrino emission.

Water-based Liquid Scintillator (WbLS) is a novel detector medium for particle physics experiments. Applications range from the use as hybrid Cherenkov/scintillation target in low-energy and accelerator neutrino experiments to large-volume neutron vetoes for dark matter detectors. Here, we present a novel WbLS featuring new components (the surfactant Triton-X and vitamin C for long-term stability), a new production recipe, and a thorough characterization of its properties. Moreover, based on neutron scattering data we are able to demonstrate that the pulse shape discrimination capabilities of this particular LS are comparable to fully-organic LAB based scintillators.

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

Cosmic rays are charged particles that are accelerated to relativistic speeds by astrophysical shocks. Numerical models have been successful in confirming the acceleration process for (quasi-)parallel shocks, which have the magnetic field aligned with the direction of the shock motion. However, the process is less clear when it comes to (quasi-)perpendicular shocks, where the field makes a large angle with the shock-normal. For such shocks, the angle between the magnetic field and flow ensures that only highly energetic particles can travel upstream at all, reducing the upstream current. This process is further inhibited for relativistic shocks, since the shock can become superluminal when the required particle velocity exceeds the speed of light, effectively inhibiting any upstream particle flow. In order to determine whether such shocks can accelerate particles, we use the particle-in-cell (PIC) method to determine what fraction of particles gets reflected initially at the shock. We then use this as input for a new simulation that combines the PIC method with grid-based magnetohydrodynamics to follow the acceleration (if any) of the particles over a larger time-period in a two-dimensional grid. We find that quasi-perpendicular, relativistic shocks are capable of accelerating particles through the DSA process, provided that the shock has a sufficiently high Alfvénic Mach number.

In this work we revisit the problem of the dynamical stability of hierarchical triple systems with applications to circumbinary planetary orbits. We carry out more than 3 10^8 numerical simulations of planets between the size of Mercury and the lower fusion boundary (13 Jupiter masses) which revolve around the center of mass of a stellar binary over long timescales. For the first time, three dimensional and eccentric planetary orbits are considered. We explore systems with a variety of binary and planetary mass ratios, binary and planetary eccentricities from 0 to 0.9 and orbital mutual inclinations ranging from 0 to 180 degrees. The simulation time is set to 10^6 planetary orbital periods. We classify the results of those long term numerical integrations into three categories: stable, unstable and mixed. We provide empirical expressions in the form of multidimensional, parameterized fits for the two borders that separate the three dynamical domains . In addition, we train a machine learning model on our data set in order to have an alternative tool of predicting the stability of circumbinary planets. Both the empirical fits and the machine learning model are tested against randomly generated circumbinary systems with very good results regarding the predictions of orbital stability. The empirical formulae are also applied to the Kepler and TESS circumbinary systems, confirming the stability of the planets in these systems. Finally, we present a REST API with a web based application for convenient access of our simulation data set.

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

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

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

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

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

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

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

In recent times, the sensitivity of low-mass direct dark matter searches has been limited by unknown low energy backgrounds close to the energy threshold of the experiments known as the low energy excess (LEE). The CRESST experiment utilises advanced cryogenic detectors constructed with different types of crystals equipped with Transition Edge Sensors (TESs) to measure signals of nuclear recoils induced by the scattering of dark matter particles in the detector. In CRESST, this low energy background manifests itself as a steeply rising population of events below 200 eV. A novel detector design named doubleTES using two identical TESs on the target crystal was studied to investigate the hypothesis that the events are sensor-related. We present the first results from two such modules, demonstrating their ability to differentiate between events originating from the crystal's bulk and those occurring in the sensor or in its close proximity.

We demonstrate that chiral symmetry breaking occurs in the confining phase of QCD-like theories with $N_c$ colors and $N_f$ flavors. Our proof is based on a novel strategy, called `downlifting', by which solutions of the 't Hooft anomaly matching and persistent mass conditions for a theory with $N_f-1$ flavors are constructed from those of a theory with $N_f$ flavors, while $N_c$ is fixed. By induction, chiral symmetry breaking is proven for any $N_f\geq p_{min}$, where $p_{min}$ is the smallest prime factor of $N_c$. The proof can be extended to $N_f <p_{min}$ under the additional assumption on the absence of phase transitions when quark masses are sent to infinity. Our results do not rely on ad-hoc assumptions on the spectrum of massless bound states.

This work aims at assessing the impact of DM self-interactions on the properties of galaxy clusters. In particular, the goal is to study the angular dependence of the cross section by testing rare (large angle scattering) and frequent (small angle scattering) SIDM models with velocity-dependent cross sections. We re-simulate six galaxy cluster zoom-in initial conditions with a dark matter only run and with a full-physics setup simulations that includes a self-consistent treatment of baryon physics. We test the dark matter only setup and the full physics setup with either collisionless cold dark matter, rare self-interacting dark matter, and frequent self-interacting dark matter models. We then study their matter density profiles as well as their subhalo population. Our dark matter only SIDM simlations agree with theoretical models, and when baryons are included in simulations, our SIDM models substantially increase the central density of galaxy cluster cores compared to full-physics simulations using collisionless dark matter. SIDM subhalo suppression in full-physics simulations is milder compared to the one found in dark matter only simulations, because of the cuspier baryionic potential that prevent subhalo disruption. Moreover SIDM with small-angle scattering significantly suppress a larger number of subhaloes compared to large angle scattering SIDM models. Additionally, SIDM models generate a broader range of subhalo concentration values, including a tail of more diffuse subhaloes in the outskirts of galaxy clusters and a population of more compact subhaloes in the cluster cores.

Understanding the explosion mechanism and hydrodynamic evolution of core-collapse supernovae is a long-standing quest in astronomy. The asymmetries caused by the explosion are encoded into the line profiles which appear in the nebular phase of the SN evolution -- with particularly clean imprints in He star explosions. Here, we carry out nine different supernova simulations of He-core progenitors, exploding them in 3D with parametrically varied neutrino luminosities using the $\texttt{Prometheus-HotB}$ code, hydrodynamically evolving the models to the homologeous phase. We then compute nebular phase spectra with the 3D NLTE spectral synthesis code $\texttt{ExTraSS}$ (EXplosive TRAnsient Spectral Simulator). We study how line widths and shifts depend on progenitor mass, explosion energy, and viewing angle. We compare the predicted line profile properties against a large set of Type Ib observations, and discuss the degree to which current neutrino-driven explosions can match observationally inferred asymmetries. With self-consistent 3D modelling -- circumventing the difficulties of representing $^{56}$Ni mixing and clumping accurately in 1D models -- we find that neither low-mass He cores exploding with high energies nor high-mass cores exploding with low energies contribute to the Type Ib SN population. Models which have line profile widths in agreement with this population give sufficiently large centroid shifts for calcium emission lines. Calcium is more strongly affected by explosion asymmetries connected to the neutron star kicks than oxygen and magnesium. Lastly, we turn to the NIR spectra from our models to investigate the potential of using this regime to look for the presence of He in the nebular phase.

Self-interacting dark matter (SIDM) has been proposed to solve small-scale problems in $\rm {\Lambda CDM}$ cosmology. In previous work, constraints on the self-interaction cross-section of dark matter have been derived assuming that the self-interaction cross-section is independent of velocity. However, a velocity-dependent cross-section is more natural in most theories of SIDM. Using idealized N-body simulations without baryons, we study merging clusters with velocity-dependent SIDM. In addition to the usual rare scattering in the isotropic limit, we also simulate these systems with anisotropic, small-angle (frequent) scatterings. We find that the collisionless brightest cluster galaxy (BCG) has an offset from the DM peak that grows at later stages. Finally, we also extend the existing upper bounds on the velocity-independent, isotropic self-interaction cross-section to the parameter space of rare and frequent velocity-dependent self-interactions by studying the central densities of dark matter-only isolated haloes. For these upper-bound parameters, the DM-BCG offsets just after the first pericentre in the dark matter-only simulations are found to be ≲10 kpc. On the other hand, because of BCG oscillations, we speculate that the distribution of BCG offsets in a relaxed cluster is a statistically viable probe. Therefore, this motivates further studies of BCG off-centring in hydrodynamic cosmological simulations.

Dark matter self-interactions may have the capability to solve or at least mitigate small-scale problems of the cosmological standard model, Lambda cold dark matter. There are a variety of self-interacting dark matter models that lead to distinguishable astrophysical predictions and hence varying success in explaining observations. Studies of dark matter (DM) density cores on various mass scales suggest a velocity-dependent scattering cross-section. In this work, we investigate how a velocity dependence alters the evolution of the DM distribution for frequent DM scatterings and compare to the velocity-independent case. We demonstrate that these cases are qualitatively different using a test problem. Moreover, we study the evolution of the density profile of idealized DM haloes and find that a velocity dependence can lead to larger core sizes and different time-scales of core formation and core collapse. In cosmological simulations, we investigate the effect of velocity-dependent self-interaction on haloes and satellites in the mass range of ≈10¹¹-$10^{14} \, \mathrm{M_\odot }$. We study the abundance of satellites, density, and shape profiles and try to infer qualitative differences between velocity-dependent and velocity-independent scatterings as well as between frequent and rare self-interactions. We find that a strongly velocity-dependent cross-section can significantly amplify the diversity of rotation curves, independent of the angular dependence of the differential cross-section. We further find that the abundance of satellites in general depends on both the velocity dependence and the scattering angle, although the latter is less important for strongly velocity-dependent cross-sections.

Despite being a fundamental property of galaxies that dictates the form of the potential, the 3D shape is intrinsically difficult to determine from observations. The improving quality of triaxial modeling methods in recent years has made it possible to measure these shapes more accurately. This study provides a comprehensive understanding of the stellar and dark matter (DM) shapes of galaxies and the connections between them as well as with other galaxy properties. Using the hydrodynamical cosmological simulation Magneticum Box4, we computed the stellar and DM shapes of galaxies at different radii. We determined their morphologies, their projected morphological and kinematic parameters, and their fractions of in-situ formed stars. The DM follows the stellar component in shape and orientation at $3R_{1/2}$, indicating that DM is heavily influenced by the baryonic potential in the inner parts of the halo. The outer DM halo is independent of the inner properties such as morphology, however, and is more closely related to the large-scale anisotropy of the gas inflow. The stellar shapes of galaxies are correlated with morphology: ellipticals feature more spherical and prolate shapes than disk galaxies. Galaxies with more rotational support are flatter, and the stellar shapes are connected to the mass distribution. In particular, more extended elliptical galaxies have larger triaxialities. Finally, the shapes can be used to constrain the in-situ fraction of stars when combined with the stellar mass. The found relations show that shapes depend on the details of the accretion history. The similarities between the inner DM and stellar shapes signal the importance of baryonic matter for DM in galaxies and will help improve dynamical models in the future. At large radii the DM shape is completely decoupled from the central galaxy and is coupled more to the large-scale inflow.

The goal of this white paper is to provide a snapshot of the data availability and data needs primarily for the Ariel space mission, but also for related atmospheric studies of exoplanets and brown dwarfs. It covers the following data-related topics: molecular and atomic line lists, line profiles, computed cross-sections and opacities, collision-induced absorption and other continuum data, optical properties of aerosols and surfaces, atmospheric chemistry, UV photodissociation and photoabsorption cross-sections, and standards in the description and format of such data. These data aspects are discussed by addressing the following questions for each topic, based on the experience of the "data-provider" and "data-user" communities: (1) what are the types and sources of currently available data, (2) what work is currently in progress, and (3) what are the current and anticipated data needs. We present a GitHub platform for Ariel-related data, with the goal to provide a go-to place for both data-users and data-providers, for the users to make requests for their data needs and for the data-providers to link to their available data. Our aim throughout the paper is to provide practical information on existing sources of data whether in databases, theoretical, or literature sources.

Draco dwarf Spheroidal galaxy (dSph) is one of the nearest and the most dark matter dominated satellites of the Milky Way. We obtained multi-epoch near-infrared (NIR, $JHK_s$) observations of the central region of Draco dSph covering a sky area of $\sim 21'\times21'$ using the WIRCam instrument at the 3.6-m Canada-France-Hawaii Telescope. Homogeneous $JHK_s$ time-series photometry for 212 RR Lyrae (173 fundamental-mode, 24 first-overtone, and 15 mixed-mode variables) and 5 Anomalous Cepheids in Draco dSph is presented and used to derive their period-luminosity relations at NIR wavelengths for the first-time. The small scatter of $\sim 0.05$~mag in these empirical relations for RR Lyrae stars is consistent with those in globular clusters and suggests a very small metallicity spread, up to $\sim0.2$~dex, among these centrally located variables. Based on empirically calibrated NIR period-luminosity-metallicity relations for RR Lyrae in globular clusters, we determined a distance modulus to Draco dSph of $\mu_\textrm{RRL} = 19.557 \pm 0.026$ mag. The calibrated $K_s$-band period-luminosity relations for Anomalous Cepheids in the Draco dSph and the Large Magellanic Cloud exhibit statistically consistent slopes but systematically different zero-points, hinting at possible metallicity dependence of $\sim-0.3$ mag~dex$^{-1}$. Finally, the apparent magnitudes of the tip of the red giant branch in $I$ and $J$ bands also agree well with their absolute calibrations with the adopted RR Lyrae distance to Draco. Our recommended $\sim1.5\%$ precise RR Lyrae distance, $D_\textrm{Draco} = 81.55 \pm 0.98 \textrm{(statistical)} \pm 1.17 \textrm{(systematic)}$~kpc, is the most accurate and precise distance to Draco dSph galaxy.

In fibers made from organic plastic scintillators, a combination of extrinsic and intrinsic effects results in the attenuation of light and thus in a position-dependent light yield. The trapping of photons can further be affected if fibers are coated with or wrapped in a light-absorbing or reflecting material to suppress optical cross-talk. These effects have frequently been studied for long (> 0.5 m) fibers, yet little data is available for shorter ones. We experimentally studied the position-dependent light yield of single-cladded Kuraray SCSF-78 fibers with lengths of < 10 cm and tested the effect of different cross-talk-preventing materials. Contrary to the often acceptable simplification that light is transmitted in the fiber core alone, we found that photons trapped by the protective cladding significantly contribute to the light transmission in short fibers. In this paper, we perform an in-depth characterization of the position-dependent light yield of fibers sputter-coated with aluminum and wrapped in aluminum foil using a double-exponential attenuation function. Finally, we compare our findings to a simple photon transport model.

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

Type Ia supernovae (SNe) remain poorly understood despite decades of investigation. Massive computationally intensive hydrodynamic simulations have been developed and run to model an ever-growing number of proposed progenitor channels. Further complicating the matter, a large number of subtypes of Type Ia SNe have been identified in recent decades. Due to the massive computational load required, inference of the internal structure of Type Ia SNe ejecta directly from observations using simulations has previously been computationally intractable. However, deep-learning emulators for radiation transport simulations have alleviated such barriers. We perform abundance tomography on 40 Type Ia SNe from optical spectra using the radiative transfer code TARDIS accelerated by the probabilistic DALEK deep-learning emulator. We apply a parametric model of potential outer ejecta structures to comparatively investigate abundance distributions and internal ionization fractions of intermediate-mass elements (IMEs) between normal and 1991T-like Type Ia SNe in the early phases. Our inference shows that the outer ejecta of 1991T-like Type Ia SNe is underabundant in the typical intermediate mass elements that heavily contribute to the spectral line formation seen in normal Type Ia SNe at early times. Additionally, we find that the IMEs present in 1991T-like Type Ia SNe are highly ionized compared to those in the normal Type Ia population. Finally, we conclude that the transition between normal and 1991T-like Type Ia SNe appears to be continuous observationally and that the observed differences come out of a combination of both abundance and ionization fractions in these SNe populations.

Context. The nuclear region of the Milky Way, within approximately −1° < l < +1° and −0.3° < b < +0.3° (i.e., |l|< 150 pc, |b|< 45 pc), is believed to host a nuclear stellar disk, co-spatial with the gaseous central molecular zone. Previous kinematical studies detected faster rotation for the stars belonging to the nuclear stellar disk, compared to the surrounding regions.
Aims: We analyze the rotation velocity of stars at the nuclear stellar disk, and compare them with its analog in a few control fields just outside this region. We limit our analysis to stars in the red clump of the color magnitude diagram, in order to be able to relate their mean de-reddened luminosity with distance along the line of sight.
Methods: We used a proper motion catalog, obtained from point spread function photometry on VISTA variables in the Vía Láctea images, to construct maps of the transverse velocity for these stars. We complemented our analysis with radial velocities from the 17th data release of the APOGEE survey.
Results: We find that the main difference between the nuclear stellar disk region and its surroundings is that at the former we see only stars moving eastward, which we believe are located in front of the Galactic center. On the contrary, in every other direction, we see the brightest red clump stars moving eastward, and the faintest ones moving westward, as expected for a rotating disk. We interpret these observations as being produced by the central molecular zone, hiding stars behind itself. What we observe is compatible with being produced by just the absence of the component at the back, without requiring the presence of a cold, fast rotating disk. This component is also not clearly detected in the newest release of the APOGEE catalog. In other words, we find no clear signature of the nuclear stellar disk as a distinct kinematical component.
Conclusions: This work highlights the need for nearby control fields when attempting to characterize the properties of the nuclear stellar disk, as the different systematics affecting this region, compared to nearby ones, might introduce spurious results. Deep, wide field and high resolution photometry of the inner 4 deg of the Milky Way is needed in order to understand the structure and kinematics of this very unique region of our Galaxy.

Full Table 1 is available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (ftp://130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/684/A214

Based on observations taken within the ESO VISTA Public Survey VVV, Program ID 179.B-2002.

We present initial results from a James Webb Space Telescope (JWST) survey of the youngest Galactic core-collapse supernova remnant, Cassiopeia A (Cas A), made up of NIRCam and MIRI imaging mosaics that map emission from the main shell, interior, and surrounding circumstellar/interstellar material (CSM/ISM). We also present four exploratory positions of MIRI Medium Resolution Spectrograph integral field unit spectroscopy that sample ejecta, CSM, and associated dust from representative shocked and unshocked regions. Surprising discoveries include (1) a weblike network of unshocked ejecta filaments resolved to ∼0.01 pc scales exhibiting an overall morphology consistent with turbulent mixing of cool, low-entropy matter from the progenitor's oxygen layer with hot, high-entropy matter heated by neutrino interactions and radioactivity; (2) a thick sheet of dust-dominated emission from shocked CSM seen in projection toward the remnant's interior pockmarked with small (∼1″) round holes formed by ≲0.″1 knots of high-velocity ejecta that have pierced through the CSM and driven expanding tangential shocks; and (3) dozens of light echoes with angular sizes between ∼0.″1 and 1' reflecting previously unseen fine-scale structure in the ISM. NIRCam observations place new upper limits on infrared emission (≲20 nJy at 3 μm) from the neutron star in Cas A's center and tightly constrain scenarios involving a possible fallback disk. These JWST survey data and initial findings help address unresolved questions about massive star explosions that have broad implications for the formation and evolution of stellar populations, the metal and dust enrichment of galaxies, and the origin of compact remnant objects.

We present the discovery of large radio shells around a massive pair of interacting galaxies and extended diffuse X-ray emission within the shells. The radio data were obtained with the Australian Square Kilometer Array Pathfinder (ASKAP) in two frequency bands centred at 944 MHz and 1.4 GHz, respectively, while the X-ray data are from the XMM-Newton observatory. The host galaxy pair, which consists of the early-type galaxies ESO 184-G042 and LEDA 418116, is part of a loose group at a distance of only 75 Mpc (redshift z = 0.017). The observed outer radio shells (diameter ~ 145 kpc) and ridge-like central emission of the system, ASKAP J1914-5433 (Physalis), are likely associated with merger shocks during the formation of the central galaxy (ESO 184-G042) and resemble the new class of odd radio circles (ORCs). This is supported by the brightest X-ray emission found offset from the centre of the Physalis system, instead centered at the less massive galaxy, LEDA 418116. The host galaxy pair is embedded in an irregular envelope of diffuse light, highlighting on-going interactions. We complement our combined radio and X-ray study with high-resolution simulations of the circumgalactic medium (CGM) around galaxy mergers from the Magneticum project to analyse the evolutionary state of the Physalis system. We argue that ORCs / radio shells could be produced by a combination of energy release from the central AGN and subsequent lightening up in radio emission by merger shocks traveling through the CGM of these systems.

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

Self-assembly is a fundamental concept in biology and of significant interest to nanotechnology. Significant progress has been made in characterizing and controlling the properties of the resulting structures, both experimentally and theoretically. However, much less is known about kinetic constraints and determinants of dynamical properties like time efficiency, although these constraints can become severe limiting factors of self-assembly processes. Here, we investigate how the time efficiency and other dynamical properties of reversible self-assembly depend on the morphology (shape) of the building blocks for systems in which the binding energy between the constituents is large. As paradigmatic examples, we stochastically simulate the self-assembly of constituents with triangular, square, and hexagonal morphology into two-dimensional structures of a specified size. We find that the constituents' morphology critically determines the assembly time and how it scales with the size of the target structure. Our analysis reveals three key structural parameters defined by the morphology: the nucleation size and attachment order, which describe the effective order of the chemical reactions by which clusters nucleate and grow, respectively, and the growth exponent, which determines how the growth rate of an emerging structure scales with its size. Using this characterization, we formulate an effective theory of the self-assembly kinetics, which we show exhibits an inherent scale invariance. This allows us to identify general scaling laws that describe the minimal assembly time as a function of the size of the target structure. We show how these insights on the kinetics of self-assembly processes can be used to design assembly schemes that could significantly increase the time efficiency and robustness of artificial self-assembly processes.

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

The phenomenon of multiple stellar populations is exacerbated in massive globular clusters, with the fraction of first-population (1P) stars a decreasing function of the cluster present-day mass. We decipher this relation in far greater detail than has been done so far. We assume (i) a fixed stellar mass threshold for the formation of second-population (2P) stars, (ii) a power-law scaling ${F}_{1{\rm{P}}}\propto {m}_{\mathrm{ecl}}^{-1}$ between the mass m _ecl of newly formed clusters and their 1P star fraction F _1P, and (iii) a constant F _1P over time. The F _1P(m _ecl) relation is then evolved up to an age of 12 Gyr for tidal field strengths representative of the entire Galactic halo. The 12 Gyr old model tracks cover the present-day distribution of Galactic globular clusters in the (mass, F _1P) space extremely well. The distribution is curtailed on its top right side by the scarcity of clusters at large Galactocentric distances and on its bottom left side by the initial scarcity of very high-mass clusters and dynamical friction. Given their distinct dissolution rates, "inner" and "outer" model clusters are offset from each other, as observed. The locus of Magellanic Clouds clusters in the (mass, F _1P) space is as expected for intermediate-age clusters evolving in a gentle tidal field. Given the assumed constancy of F _1P, we conclude that 2P stars do not necessarily form centrally concentrated. We infer a minimum mass of 4 · 10⁵ M _⊙ for multiple-population clusters at secular evolution onset. This high-mass threshold severely limits the number of 2P stars lost from evolving clusters, thereby fitting the low 2P star fraction of the Galactic halo field.

The Orion Nebula Cluster (ONC) is the closest site of very young ($\sim$ 1 Myrs) massive star formation. The ONC hosts more than 1600 young and X-ray bright stars with masses ranging from $\sim$ 0.1 to 35 $M_\odot$. The Chandra HETGS Orion Legacy Project observed the ONC with the Chandra high energy transmission grating spectrometer (HETGS) for $2.1\,$Ms. We describe the spectral extraction and cleaning processes necessary to separate overlapping spectra. We obtained 36 high resolution spectra which includes a high brilliance X-ray spectrum of $\theta^1$ Ori C with over 100 highly significant X-ray lines. The lines show Doppler broadening between 300 and $400\;\mathrm{km}\;\mathrm{s}^{-1}$. Higher spectral diffraction orders allow us to resolve line components of high Z He-like triplets in $\theta^1$ Ori C with unprecedented spectral resolution. Long term light curves spanning $\sim$20 years show all stars to be highly variable, including the massive stars. Spectral fitting with thermal coronal emission line models reveals that most sources show column densities of up to a few times $10^{22}\,$cm$^{-2}$ and high coronal temperatures of 10 to 90 MK. We observe a bifurcation of the high temperature component where some stars show a high component of 40 MK, while others show above 60 MK indicating heavy flaring activity. Some lines are resolved with Doppler broadening above our threshold of $\sim200\;\mathrm{km}\;\mathrm{s}^{-1}$, up to $500\;\mathrm{km}\;\mathrm{s}^{-1}$. This data set represents the largest collection of HETGS high resolution X-ray spectra from young pre-MS stars in a single star-forming region to date.

The calculation of loop corrections to the correlation functions of quantum fields during inflation or in the de Sitter background presents greater challenges than in flat space due to the more complicated form of the mode functions. While in flat space highly sophisticated approaches to Feynman integrals exist, similar tools still remain to be developed for cosmological correlators. However, usually only their late-time limit is of interest. We introduce the method-of-region expansion for cosmological correlators as a tool to extract the late-time limit, and illustrate it with several examples for the interacting, massless, minimally coupled scalar field in de Sitter space. In particular, we consider the in-in correlator «ϕ²(η, q)ϕ(η, k₁)ϕ(η, k₂)», whose region structure is relevant to anomalous dimensions and matching coefficients in Soft de Sitter effective theory.

Aims: We have performed the first broadband study of Mrk 421 from radio to TeV gamma rays with simultaneous measurements of the X-ray polarization from IXPE.
Methods: The data were collected as part of an extensive multiwavelength campaign carried out between May and June 2022 using MAGIC, Fermi-LAT, NuSTAR, XMM-Newton, Swift, and several optical and radio telescopes to complement IXPE data.
Results: During the IXPE exposures, the measured 0.2-1 TeV flux was close to the quiescent state and ranged from 25% to 50% of the Crab Nebula without intra-night variability. Throughout the campaign, the very high-energy (VHE) and X-ray emission are positively correlated at a 4σ significance level. The IXPE measurements reveal an X-ray polarization degree that is a factor of 2-5 higher than in the optical/radio bands; that implies an energy-stratified jet in which the VHE photons are emitted co-spatially with the X-rays, in the vicinity of a shock front. The June 2022 observations exhibit a rotation of the X-ray polarization angle. Despite no simultaneous VHE coverage being available during a large fraction of the swing, the Swift-XRT monitoring reveals an X-ray flux increase with a clear spectral hardening. This suggests that flares in high synchrotron peaked blazars can be accompanied by a polarization angle rotation, as observed in some flat spectrum radio quasars. Finally, during the polarization angle rotation, NuSTAR data reveal two contiguous spectral hysteresis loops in opposite directions (clockwise and counterclockwise), implying important changes in the particle acceleration efficiency on approximately hour timescales.

All data shown in Figs. 1, 2, 5, 7, and 8 are available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (ftp://130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/684/A127

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

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

Vorticity has recently been suggested to be a property of highly spinning black holes. The connection between vorticity and limiting spin represents a universal feature shared by objects of maximal microstate entropy, so-called saturons. Using Q -ball-like saturons as a laboratory for black holes, we study the collision of two such objects and find that vorticity can have a large impact on the emitted radiation as well as on the charge and angular momentum of the final configuration. As black holes belong to the class of saturons, we expect that the formation of vortices can cause similar effects in black hole mergers, leading to macroscopic deviations in gravitational radiation. This could leave unique signatures detectable with upcoming gravitational-wave searches, which can thereby serve as a portal to macroscopic quantum effects in black holes.

General circulation models of gas giant exoplanets predict equatorial jets that drive inhomogeneities across the planetary atmosphere. We studied the transmission spectrum of the hot Jupiter WASP-127b during one transit in the K band with CRIRES+. Telluric and stellar signals were removed from the data using SYSREM and the planetary signal was investigated using the cross-correlation (CCF) technique. After detecting a spectral signal indicative of atmospheric inhomogeneities, we employed a Bayesian retrieval framework with a 2D modelling approach tailored to address this scenario. We detected strong signals of H$_2$O and CO, which exhibited not one but two distinct CCF peaks. The double peaked signal can be explained by a supersonic equatorial jet and muted signals at the poles, with the two peaks representing the signals from the planet's morning and evening terminators, respectively. We calculated an equatorial jet velocity of $7.7\pm0.2$km/s from our retrieved overall equatorial velocity and the planet's tidally locked rotation, and derive distinct atmospheric properties for the two terminators as well as the polar region. The evening terminator is found to be hotter than the morning terminator by $175^{+116}_{-133}$K and the muted signals from the poles can be explained by significantly lower temperatures or a high cloud deck. Our retrieval yields a solar C/O ratio and metallicity and challenges previous studies of WASP-127b's atmosphere. The presence of a double peaked signal highlights the importance of accounting for planetary 3D structure during interpretation of atmospheric signals. The measured supersonic jet velocity and the lack of signal from the polar regions, representing a detection of latitudinal inhomogeneity in a spatially unresolved target, showcases the power of high-resolution transmission spectroscopy for the characterization of global circulation in exoplanets.

Due to rapidly improving quantum computing hardware, Hamiltonian simulations of relativistic lattice field theories have seen a resurgence of attention. This computational tool requires turning the formally infinite-dimensional Hilbert space of the full theory into a finite-dimensional one. For gauge theories, a widely used basis for the Hilbert space relies on the representations induced by the underlying gauge group, with a truncation that keeps only a set of the lowest dimensional representations. This works well at large bare gauge coupling, but becomes less efficient at small coupling, which is required for the continuum limit of the lattice theory. In this work, we develop a new basis suitable for the simulation of an SU(2) lattice gauge theory in the maximal tree gauge. In particular, we show how to perform a Hamiltonian truncation so that the eigenvalues of both the magnetic and electric gauge-fixed Hamiltonian are mostly preserved, which allows for this basis to be used at all values of the coupling. Little prior knowledge is assumed, so this may also be used as an introduction to the subject of Hamiltonian formulations of lattice gauge theories.

We review the phenomenology of classical Cepheids (CCs), Anomalous Cepheids (ACs) and type II Cepheids (TIICs) in the Milky Way (MW) and in the Magellanic Clouds (MCs). We also examine the Hertzsprung progression in different stellar systems by using the shape of I-band light curves (Fourier parameters) and observables based on the difference in magnitude and in phase between the bump and the minimum in luminosity. The distribution of Cepheids in optical and in optical-near infrared (NIR) color-magnitude diagrams is investigated to constrain the topology of the instability strip. The use of Cepheids as tracers of young (CCs), intermediate (ACs) and old (TIICs) stellar populations are brought forward by the comparison between observations (MCs) and cluster isochrones covering a broad range in stellar ages and in chemical compositions. The different diagnostics adopted to estimate individual distances (period-luminosity, period-Wesenheit, period-luminosity-color relations) are reviewed together with pros and cons in the use of fundamental and overtones, optical and NIR photometric bands, and reddening free pseudo magnitudes (Wesenheit). We also discuss the use of CCs as stellar tracers and the radial gradients among the different groups of elements (iron, α , neutron-capture) together with their age-dependence. Finally, we briefly outline the role that near-future space and ground-based facilities will play in the astrophysical and cosmological use of Cepheids.

The collective properties of star clusters are investigated using a simulation of the collision between two dwarf galaxies. The characteristic power law of the cluster mass function, N(M), with a slope dlog N/dlog M ~ -1, is present from cluster birth and remains throughout the simulation. The maximum mass of a young cluster scales with the star formation rate (SFR). The relative average minimum separation, R(M)= N(M)^{1/p}D_min(M)/D(M_low), for average minimum distance D_min(M) between clusters of mass M, and for lowest mass, M_low, measured in projection (p=2) or three dimensions (p=3), has a negative slope, dlog R/dlog M ~ -0.2, for all masses and ages. This agrees with observations of R(M) in low-mass galaxies studied previously. Like the slope of N(M), R}(M) is apparently a property of cluster birth for dwarf galaxies that does not depend on SFR or time. The negative slope for R(M) implies that more massive clusters are centrally concentrated relative to lower mass clusters throughout the entire mass range. Cluster growth through coalescence is also investigated. The ratio of the kinetic to potential energy of all near-neighbor clusters is generally large, but a tail of low values in the distribution of this ratio suggests that a fraction of the clusters merge, ~8% by number throughout the ~300 Myr of the simulation and up to 60% by mass for young clusters in their first 10 Myr, scaling with the SFR above a certain threshold.

Dust grains play a significant role in several astrophysical processes, including gas/dust dynamics, chemical reactions, and radiative transfer. Replenishment of small-grain populations is mainly governed by fragmentation during pair-wise collisions between grains. The wide spectrum of fragmentation outcomes, from complete disruption to erosion and/or mass transfer, can be modelled by the general non-linear fragmentation equation. Efficiently solving this equation is crucial for an accurate treatment of the dust fragmentation in numerical modelling. However, similar to dust coagulation, numerical errors in current fragmentation algorithms employed in astrophysics are dominated by the numerical over-diffusion problem -- particularly in 3D hydrodynamic simulations where the discrete resolution of the mass density distribution tends to be highly limited. With this in mind, we have derived the first conservative form of the general non-linear fragmentation with a mass flux highlighting the mass transfer phenomenon. Then, to address cases of limited mass density resolution, we applied a high-order discontinuous Galerkin scheme to efficiently solve the conservative fragmentation equation with a reduced number of dust bins. An accuracy of 0.1 -1% is reached with 20 dust bins spanning a mass range of 9 orders of magnitude.

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

It has long been known that the maximal cut of the equal-mass four-loop banana integral is a period of a family of Calabi-Yau threefolds that depends on the kinematic variable $z=m^2/p^2$. We show that it can also be interpreted as a period of a family of genus-two curves. We do this by introducing a general Calabi-Yau-to-curve correspondence, which in this case locally relates the original period of the family of Calabi-Yau threefolds to a period of a family of genus-two curves that varies holomorphically with the kinematic variable $z$. In addition to working out the concrete details of this correspondence for the equal-mass four-loop banana integral, we outline when we expect a correspondence of this type to hold.

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

While Bayesian inference techniques are standard in cosmological analyses, it is common to interpret resulting parameter constraints with a frequentist intuition. This intuition can fail, e.g. when marginalizing high-dimensional parameter spaces onto subsets of parameters, because of what has come to be known as projection effects or prior volume effects. We present the method of Informed Total-Error-Minimizing (ITEM) priors to address this. An ITEM prior is a prior distribution on a set of nuisance parameters, e.g. ones describing astrophysical or calibration systematics, intended to enforce the validity of a frequentist interpretation of the posterior constraints derived for a set of target parameters, e.g. cosmological parameters. Our method works as follows: For a set of plausible nuisance realizations, we generate target parameter posteriors using several different candidate priors for the nuisance parameters. We reject candidate priors that do not accomplish the minimum requirements of bias (of point estimates) and coverage (of confidence regions among a set of noisy realizations of the data) for the target parameters on one or more of the plausible nuisance realizations. Of the priors that survive this cut we select the ITEM prior as the one that minimizes the total error of the marginalized posteriors of the target parameters. As a proof of concept, we apply our method to the Density Split Statistics (DSS) measured in Dark Energy Survey Year 1 data. We demonstrate that the ITEM priors substantially reduce prior volume effects that otherwise arise and allow sharpened yet robust constraints on the parameters of interest.

The joint probability distribution of matter overdensity and galaxy counts in cells is a powerful probe of cosmology, and the extent to which variance in galaxy counts at fixed matter density deviates from Poisson shot noise is not fully understood. The lack of informed bounds on this stochasticity is currently the limiting factor in constraining cosmology with the galaxy-matter PDF. We investigate stochasticity in the conditional distribution of galaxy counts at fixed matter density and present a halo occupation distribution (HOD)-based approach for obtaining plausible ranges for stochasticity parameters. To probe the high-dimensional space of possible galaxy-matter connections, we derive HODs which conserve linear galaxy bias and number density to produce redMaGiC-like galaxy catalogs within the AbacusSummit suite of N-body simulations. We study the impact of individual HOD parameters and cosmology on stochasticity and perform a Monte Carlo search in HOD parameter space, subject to the constraints on bias and density. In mock catalogs generated by the selected HODs, shot noise in galaxy counts spans both sub-Poisson and super-Poisson values, ranging from 80% to 133% of Poisson variance at mean matter density. Nearly all derived HODs show a positive relationship between local matter density and stochasticity. For galaxy catalogs with higher stochasticity, quadratic galaxy bias is required for an accurate description of the conditional PDF of galaxy counts at fixed matter density. The presence of galaxy assembly bias also substantially extends the range of stochasticity in the super-Poisson direction. This HOD-based approach leverages degrees of freedom in the galaxy-halo connection to obtain informed bounds on model nuisance parameters and can be adapted to other parametrizations of stochasticity, in particular to motivate prior ranges for cosmological analyses.

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

Our Sun lies within 300 parsecs of the 2.7-kiloparsecs-long sinusoidal chain of dense gas clouds known as the Radcliffe Wave¹. The structure's wave-like shape was discovered using three-dimensional dust mapping, but initial kinematic searches for oscillatory motion were inconclusive^2-7. Here we present evidence that the Radcliffe Wave is oscillating through the Galactic plane while also drifting radially away from the Galactic Centre. We use measurements of line-of-sight velocity⁸ for ¹²CO and three-dimensional velocities of young stellar clusters to show that the most massive star-forming regions spatially associated with the Radcliffe Wave (including Orion, Cepheus, North America and Cygnus X) move as though they are part of an oscillating wave driven by the gravitational acceleration of the Galactic potential. By treating the Radcliffe Wave as a coherently oscillating structure, we can derive its motion independently of the local Galactic mass distribution, and directly measure local properties of the Galactic potential as well as the Sun's vertical oscillation period. In addition, the measured drift of the Radcliffe Wave radially outwards from the Galactic Centre suggests that the cluster whose supernovae ultimately created today's expanding Local Bubble⁹ may have been born in the Radcliffe Wave.

Neutrinos can undergo fast flavor conversions (FFCs) within extremely dense astrophysical environments, such as core-collapse supernovae (CCSNe) and neutron star mergers (NSMs). In this study, we explore FFCs in a multienergy neutrino gas, revealing that when the FFC growth rate significantly exceeds that of the vacuum Hamiltonian, all neutrinos (regardless of energy) share a common survival probability dictated by the energy-integrated neutrino spectrum. We then employ physics-informed neural networks (PINNs) to predict the asymptotic outcomes of FFCs within such a multienergy neutrino gas. These predictions are based on the first two moments of neutrino angular distributions for each energy bin, typically available in state-of-the-art CCSN and NSM simulations. Our PINNs achieve errors as low as ≲6 % and ≲18 % for predicting the number of neutrinos in the electron channel and the relative absolute error in the neutrino moments, respectively.

In a pedagogical manner, we review recent developments in the investigation of the Emergence Proposal. Although it is fair to say that this idea is still at an exploratory level and a fully coherent picture has yet to be developed, we put it into perspective to previous work on the swampland program and on emergence in QG. In view of the emergent string conjecture, we argue and provide evidence that it is not the emergent string but rather the decompactification limit which is a natural candidate for the potential realization of the Emergence Proposal. This resonates in a compelling way with old ideas of emergence in M(-atrix) theory and gives rise to a number of further speculations.

The deposition of energy and momentum by supernova explosions has been subject to numerous studies in the past few decades. However, while there has been some work that focused on the transition from the adiabatic to the radiative stage of a supernova remnant (SNR), the late radiative stage and merging with the interstellar medium (ISM) have received little attention. Here, we use three-dimensional, hydrodynamic simulations, focusing on the evolution of SNRs during the radiative phase, considering a wide range of physical explosion parameters ( ${n}_{{\rm{H}},\mathrm{ISM}}\in \left[0.1,100\right]{\mathrm{cm}}^{-3}$ and ${E}_{\mathrm{SN}}\in \left[1,14\right]\times {10}^{51}\,\mathrm{erg}$ ). We find that the radiative phase can be subdivided in four stages: A pressure-driven snowplow phase, during which the hot overpressurized bubble gas is evacuated and pushed into the cold shell; a momentum-conserving snowplow phase that is accompanied by a broadening of the shell; an implosion phase where cold material from the back of the shell is flooding the central vacuum; and a final cloud phase, during which the imploding gas is settling as a central, compact overdensity. The launching timescale for the implosion ranges from a few 100 kyr to a few Myr, while the cloud formation timescale ranges from a few to about 10 Myr. The highly chemically enriched clouds can become massive (M _cl ∼ 10³–10⁴ M _⊙) and self-gravitating within a few Myr after their formation, providing an attractive, novel pathway for supernova-induced star and planet formation in the ISM.

The apparent tension between the luminosity functions of red supergiant (RSG) stars and of RSG progenitors of Type II supernovae (SNe) is often referred to as the RSG problem and it motivated some to suggest that many RSGs end their life without an SN explosion. However, the luminosity functions of RSG SN progenitors presented so far were biased to high luminosities, because the sensitivity of the search was not considered. Here, we use limiting magnitudes to calculate a bias-corrected RSG progenitor luminosity function. We find that only (36 ± 11)% of all RSG progenitors are brighter than a bolometric magnitude of ‑7 mag, a significantly smaller fraction than (56 ± 5)% quoted by Davies & Beasor. The larger uncertainty is due to the relatively small progenitor sample, while uncertainties on measured quantities such as magnitudes, bolometric corrections, extinction, or SN distances, only have a minor impact, as long as they fluctuate randomly for different objects in the sample. The bias-corrected luminosity functions of RSG SN progenitors and Type M supergiants in the Large Magellanic Cloud are consistent with each other, as also found by Davies & Beasor for the uncorrected luminosity function. The RSG progenitor luminosity function, hence, does not imply the existence of failed SNe. The presented statistical method is not limited to progenitor searches, but applies to any situation in which a measurement is done for a sample of detected objects, but the probed quantity or property can only be determined for part of the sample.