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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Full Tables 1 and 2 are available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (ftp://130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/687/A312

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

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

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

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

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

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

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

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

Context. The Kepler high-precision planetary sample has revealed a `radius valley' separating compact super-Earths from sub-Neptunes with lower densities. Super-Earths are generally assumed to be rocky planets that were probably born in situ, while the composition and formation of sub-Neptunes remains debated. Numerous statistical studies have explored planetary and stellar properties and their correlations to provide observational clues. However, no conclusive result on the origin of the radius valley or the composition of sub-Neptunes has been derived to date.
Aims: To provide more constraints, our aim is to investigate the distributions of the orbital spacing of sub-Neptunes and super-Earth planets in Kepler systems and compare their distributions with theoretical predictions of planet pairs of different formation pathways and compositions in synthetic planetary systems.
Methods: Based on the Kepler planetary sample, we derived the distributions of period ratios of sub-Neptune and super-Earth planet pairs. Using synthetic planetary systems generated by the Generation III Bern Model, we also obtained theoretical predictions of period ratio distributions of planet pairs of different compositions and origins.
Results: We find that Kepler sub-Neptune pairs show a significant preference to be near first-order mean motion resonances by a factor of 1.7<sub>−0.3</sub><sup>+0.3</sup>. This is smaller than the model predictions for `water-rich' pairs but larger than that of `water-poor' pairs by confidence levels of ~2σ. Kepler super-Earth pairs show no significant preference for mean motion resonances from a random distribution. The derived normalised fraction of near first-order resonances of actual Kepler super-Earth pairs is consistent with that of synthetic water-poor planet pairs but significantly (≳3σ) smaller than that of synthetic water-rich planet pairs.
Conclusions: The orbital migration has been more important for sub-Neptunes than for super-Earths, suggesting a partial ex situ formation of the former and an origin of the radius valley caused in part by distinct formation pathways. However, the model comparisons also show that sub-Neptunes in Kepler multiple systems are not likely (~2σ) to all be water-rich planets born ex situ but a mixture of the two (in situ and ex situ) pathways. Whereas, Kepler super-Earth planets are predominantly composed of water-poor planets that were born inside the ice line, likely through a series of giant impacts without large-scale migration.

Context. Gas in protoplanetary disks mostly cools via thermal accommodation with dust particles. Thermal relaxation is thus highly sensitive to the local dust size distributions and the spatial distribution of the grains. So far, the interplay between thermal relaxation and gas turbulence has not been dynamically modeled in hydrodynamic simulations of protoplanetary disks with dust.
Aims: We aim to study the effects of the vertical shear instability (VSI) on the thermal relaxation times, and vice versa. We are particularly interested in the influence of the initial dust grain size on the VSI and whether the emerging turbulence is sustained over long timescales.
Methods: We ran three axisymmetric hydrodynamic simulations of a protoplanetary disk including four dust fluids that initially resemble MRN size distributions of different initial grain sizes. From the local dust densities, we calculated the thermal accommodation timescale of dust and gas and used the result as the thermal relaxation time of the gas in our simulation. We included the effect of dust growth by applying the monodisperse dust growth rate and the typical growth limits.
Results: We find that the emergence of the VSI is strongly dependent on the initial dust grain size. Coagulation also counteracts the emergence of hydrodynamic turbulence in our simulations, as shown by others before. Starting a simulation with larger grains (100 μm) generally leads to a less turbulent outcome. While the inner disk regions (within ∼70 au) develop turbulence in all three simulations, we find that the simulations with larger particles do not develop VSI in the outer disk.
Conclusions: Our simulations with dynamically calculated thermal accommodation times based on the drifting and settling dust distribution show that the VSI, once developed in a disk, can be sustained over long timescales, even if grain growth is occurring. The VSI corrugates the dust layer and even diffuses the smaller grains into the upper atmosphere, where they can cool the gas. Whether the instability can emerge for a specific stratification depends on the initial dust grain sizes and the initial dust scale height. If the grains are initially ≳100 μm and if the level of turbulence is initially assumed to be low, we find no VSI turbulence in the outer disk regions.

We describe an attempt of string theoretic derivation of the Gibbons-Hawking entropy. Despite not admitting a de Sitter vacuum, the string theory, by the power of open-close correspondence, captures the Gibbons-Hawking entropy as the entropy of Chan-Paton species on a de Sitter-like state obtained via $D$-branes. Moreover, this derivation sheds a new light at the origin of the area-form, since the equality takes place for a critical 't Hooft coupling for which the species entropy of open strings saturates the area-law unitarity bound.

Observations and high-resolution hydrodynamical simulations indicate that massive star clusters assemble hierarchically from subclusters with a universal power-law cluster mass function. We study the consequences of such assembly for the formation of intermediate-mass black holes (IMBHs) at low metallicities (Z = 0.01 Z_⊙) with our updated N-body code BIFROST based on the hierarchical fourth-order forward integrator. BIFROST integrates few-body systems using secular and regularized techniques including post-Newtonian equations of motion up to order PN3.5 and gravitational-wave recoil kicks for BHs. Single stellar evolution is treated using the fast population synthesis code SEVN. We evolve three cluster assembly regions with N_tot = 1.70-2.35 × 10⁶ stars following a realistic IMF in ~1000 subclusters for t = 50 Myr. IMBHs with masses up to m_• ~ 2200 M_⊙ form rapidly mainly via the collapse of very massive stars (VMSs) assembled through repeated collisions of massive stars followed by growth through tidal disruption events and BH mergers. No IMBHs originate from the stars in the initially most massive clusters. We explain this by suppression of hard massive star binary formation at high velocity dispersions and the competition between core collapse and massive star lifetimes. Later the IMBHs form subsystems resulting in gravitational-wave BH-BH, IMBH-BH, and IMBH-IMBH mergers with a m_• ~ 1000 M_⊙ gravitational-wave detection being the observable prediction. Our simulations indicate that the hierarchical formation of massive star clusters in metal poor environments naturally results in formation of potential seeds for supermassive black holes.

We present a machine-learning-based model for the total density profiles of subhaloes with masses $M \gtrsim 7\times 10^8\, h^{-1}{\rm M}_\odot$ in the IllustrisTNG100 simulation. The model is based on an interpretable variational encoder (IVE) which returns the independent factors of variation in the density profiles within a low-dimensional representation, as well as the predictions for the density profiles themselves. The IVE returns accurate and unbiased predictions on all radial ranges, including the outer region profile where the subhaloes experience tidal stripping; here its fit accuracy exceeds that of the commonly used Einasto profile. The IVE discovers three independent degrees of freedom in the profiles, which can be interpreted in terms of the formation history of the subhaloes. In addition to the two parameters controlling the normalization and inner shape of the profile, the IVE discovers a third parameter that accounts for the impact of tidal stripping on to the subhalo outer profile; this parameter is sensitive to the mass loss experienced by the subhalo after its infall on to its parent halo. Baryonic physics in the IllustrisTNG galaxy formation model does not impact the number of degrees of freedom identified in the profile compared to the pure dark matter expectations, nor their physical interpretation. Our newly proposed profile fit can be used in strong lensing analyses or other observational studies which aim to constrain cosmology from small-scale structures.

Lambda- Cold Dark Matter (LambdaCDM) has been successful at explaining the large-scale structures in the universe but faces severe issues on smaller scales when compared to observations. Introducing self-interactions between dark matter particles claims to provide a solution to the small-scale issues in the LambdaCDM simulations while being consistent with the observations at large scales. The existence of the energy region in which these self-interactions between dark matter particles come close to saturating the S-wave unitarity bound can result in the formation of dark matter bound states called darkonium. In this scenario, all the low energy scattering properties are determined by a single parameter, the inverse scattering length gamma. In this work, we set bounds on gamma by studying the impact of darkonium on the observations at direct detection experiments using data from CRESST-III and XENON1T. The exclusion limits on gamma are then subsequently converted to exclusion limits on the self-interaction cross-section and compared with the constraints from astrophysics and N-body simulations.

RR Lyrae stars are excellent tracers of stellar populations for old, metal-poor components in the Milky Way Galaxy and the Local Group. Their luminosities have a metallicity-dependence, but determining spectroscopic [Fe/H] metallicities for RR Lyrae stars, especially at distances outside the solar neighbourhood, is challenging. Using 40 RRLs with metallicities derived from both Fe(II) and Fe(I) abundances, we verify the calibration between the [Fe/H] of RR Lyrae stars from the Calcium triplet. Our calibration is applied to all RR Lyrae stars with Gaia RVS spectra in Gaia DR3 as well as to 80 stars in the inner Galaxy from the BRAVA-RR survey. The co-added Gaia RVS RR Lyrae spectra provide RR Lyrae metallicities with an uncertainty of 0.25~dex, which is a factor of two improvement over the Gaia photometric RR Lyrae metallicities. Within our Galactic bulge RR Lyrae star sample, we find a dominant fraction with low energies without a prominent rotating component. Due to the large fraction of such stars, we interpret these stars as belonging to the $in-situ$ metal-poor Galactic bulge component, although we can not rule out that a fraction of these belong to an ancient accretion event such as Kraken/Heracles.

Life is a phenomenon far from equilibrium, which suggests that its emergence was also decisively shaped and driven by the non-equilibrium systems present around 4 billion years ago. These include geothermal systems with heat flows through thin, water-filled rock fractures, which drive the convection of water and the thermophoresis of dissolved substances. My thesis shows that their interplay drives a highly selective enrichment of components that offers unique opportunities for prebiotic systems on the path to the emergence of life.

Prebiotic chemistry usually requires specific compositions and high concentrations of compounds for the desired products to form at the end of sometimes complex synthesis pathways whilst suppressing undesired side reactions. A possible solution for this is provided by thermophoretic accumulation, which can robustly separate amino acids, nucleotides and RNA building blocks (Chapter 3) even for mass-identical compounds. In networks of connected rock cracks, as found in natural systems, this purification is further enhanced. The resulting diversity of compositions increases reaction yields by up to four orders of magnitude, for example for peptide polymerization, thus demonstrating how temperature gradients in rock cracks can foster prebiotic chemistry.[...]

. We study the excitation spectrum of light and strange mesons in diffractive scattering. We identify different hadron resonances through partial wave analysis, which inherently relies on analysis models. Besides statistical uncertainties, the model dependence of the analysis introduces dominant sys tematic uncertainties. We discuss several of their sources for the π⁻π⁻π⁺ and K⁰_S K^¯ final states and present methods to reduce them. We have developed a new approach exploiting a-priori knowledge of signal continuity over adjacent final-state-mass bins to stably fit a large pool of partial-waves to our data, al lowing a clean identification of very small signals in our large data sets. For two-body final states of scalar particles, such as K⁰_S K^¯, mathematical ambigui ties in the partial-wave decomposition lead to the same intensity distribution for different combinations of amplitude values. We will discuss these ambiguities and present solutions to resolve or at least reduce the number of possible so lutions. Resolving these issues will allow for a complementary analysis of the a_J-like resonance sector in these two final states.

We review observational, experimental and theoretical results related to Dark Matter.

Radio emission has been detected from tens of white dwarfs, in particular in accreting systems. Additionally, radio emission has been predicted as a possible outcome of a planetary system around a white dwarf. We searched for 3 GHz radio continuum emission in 846 000 candidate white dwarfs previously identified in Gaia using the Very Large Array Sky Survey (VLASS) Epoch 1 Quick Look Catalogue. We identified 13 candidate white dwarfs with a counterpart in VLASS within 2 arcsec. Five of those were found not to be white dwarfs in follow-up or archival spectroscopy, whereas seven others were found to be chance alignments with a background source in higher resolution optical or radio images. The remaining source, WDJ204259.71+152108.06, is found to be a white dwarf and M-dwarf binary with an orbital period of 4.1 d and long-term stochastic optical variability, as well as luminous radio and X-ray emission. For this binary, we find no direct evidence of a background contaminant, and a chance alignment probability of only ≈2 per cent. However, other evidence points to the possibility of an unfortunate chance alignment with a background radio and X-ray emitting quasar, including an unusually poor Gaia DR3 astrometric solution for this source. With at most one possible radio emitting white dwarf found, we conclude that strong (≳1-3 mJy) radio emission from white dwarfs in the 3 GHz band is virtually non-existent outside of interacting binaries.

Time-delay cosmography is a powerful technique to constrain cosmological parameters, particularly the Hubble constant ($H_{0}$). The TDCOSMO collaboration is performing an ongoing analysis of lensed quasars to constrain cosmology using this method. In this work, we obtain constraints from the lensed quasar~WGD$\,$2038$-$4008~using new time-delay measurements and previous mass models by TDCOSMO. This is the first TDCOSMO lens to incorporate multiple lens modeling codes and the full time-delay covariance matrix into the cosmological inference. The models are fixed before the time delay is measured, and the analysis is performed blinded with respect to the cosmological parameters to prevent unconscious experimenter bias. We obtain $D_{\Delta t} = 1.68^{+0.40}_{-0.38}$ Gpc using two families of mass models, a power-law describing the total mass distribution, and a composite model of baryons and dark matter, although the composite model is disfavored due to kinematics constraints. In a flat $\Lambda$CDM cosmology, we constrain the Hubble constant to be $H_{0} = 65^{+23}_{-14}\, \rm km\ s^{-1}\,Mpc^{-1}$. The dominant source of uncertainty comes from the time delays, due to the low variability of the quasar. Future long-term monitoring, especially in the era of the Vera C. Rubin Observatory's Legacy Survey of Space and Time, could catch stronger quasar variability and further reduce the uncertainties. This system will be incorporated into an upcoming hierarchical analysis of the entire TDCOSMO sample, and improved time delays and spatially-resolved stellar kinematics could strengthen the constraints from this system in the future.

Theoretical models of structure formation predict the presence of a hot gaseous atmosphere around galaxies. While this hot circum-galactic medium (CGM) has been observationally confirmed through UV absorption lines, the detection of its direct X-ray emission remains scarce. We investigate theoretical predictions of the intrinsic CGM X-ray surface brightness (SB) using simulated galaxies and connect them to their global properties such as gas temperature, hot gas fraction and stellar mass. We select a sample of galaxies from the ultra-high resolution ($48\ \rm{cMpc\, h^{-1}}$) cosmological volume of the Magneticum Pathfinder set of hydrodynamical cosmological simulations. We classify them as star-forming (SF) or quiescent (QU) based on their specific star-formation rate. For each galaxy we generate X-ray mock data using the X-ray photon simulator PHOX, from which we obtain SB profiles out to the virial radius for different X-ray emitting components, namely gas, active galactic nuclei and X-ray binaries (XRBs). We fit a $\beta$-profile to each galaxy and observe trends between its slope and global quantities of the simulated galaxy. We find marginal differences between the average total SB profile of the CGM in SF and QU galaxies, with the contribution from hot gas being the largest ($>50\%$) at radii $r>0.05\,R_{\rm{vir}}$. The contribution from X-ray binaries (XRBs) equals the gas contribution for small radii and is non-zero for large radii. The galaxy population shows positive correlations between global properties and normalization of the SB profile. The slope of fitted $\beta$-profiles correlates strongly with the total gas luminosity, which in turn shows strong connections to the current accretion rate of the central super-massive black hole (SMBH).

Context. Planets orbiting members of open or globular clusters offer a great opportunity to study exoplanet populations systematically, as stars within clusters provide a mostly homogeneous sample, at least in chemical composition and stellar age. However, even though there have been coordinated efforts to search for exoplanets in stellar clusters, only a small number of planets have been detected. One successful example is the seven-year radial velocity (RV) survey `Search for giant planets in M 67' of 88 stars in the open cluster M 67, which led to the discovery of five giant planets, including three close-in (P < 10 days) hot-Jupiters.
Aims: In this work, we continue and extend the observation of stars in M 67, with the aim being to search for additional planets.
Methods: We conducted spectroscopic observations with the Habitable Planet Finder (HPF), HARPS, HARPS-North, and SOPHIE spectrographs of 11 stars in M 67. Six of our targets showed a variation or long-term trends in their RV during the original survey, while the other five were not observed in the original sample, bringing the total number of stars to 93.
Results: An analysis of the RVs reveals one additional planet around the turn-off point star S1429 and provides solutions for the orbits of stellar companions around S2207 and YBP2018. S1429 b is a warm-Jupiter on a likely circular orbit with a period of [77.48_-0.19^+0.18] days and a minimum mass of M sin i = 1.80 ± 0.2 M_J. We update the hot-Jupiter occurrence rate in M 67 to include the five new stars, deriving [4.2_-2.3^+4.1 %] when considering all stars, and [5.4_-3.0^+5.1 %] if binary star systems are removed.

A growing number of supernovae (SNe) are now known to exhibit evidence for significant interaction with a dense, pre-existing, circumstellar medium (CSM). SNe Ibn comprise one such class that can be characterized by both rapidly evolving light curves and persistent narrow He I lines. The origin of such a dense CSM in these systems remains a pressing question, specifically concerning the progenitor system and mass-loss mechanism. In this paper, we present multiwavelength data of the Type Ibn SN 2020nxt, including HST/STIS ultraviolet spectra. We fit the data with recently updated CMFGEN models designed to handle configurations for SNe Ibn. The UV coverage yields strong constraints on the energetics and, when combined with the CMFGEN models, offer new insight on potential progenitor systems. We find the most successful model is a ≲4 M_⊙ helium star that lost its $\sim 1\, {\rm M}_\odot$ He-rich envelope in the years preceding core collapse. We also consider viable alternatives, such as a He white dwarf merger. Ultimately, we conclude at least some SNe Ibn do not arise from single, massive (>30 M_⊙) Wolf-Rayet-like stars.

Context. Disc winds and planet-disc interactions are two crucial mechanisms that define the structure, evolution, and dispersal of protoplanetary discs. While winds are capable of removing material from discs, eventually leading to their dispersal, massive planets can shape their disc by creating sub-structures such as gaps and spiral arms.
Aims: We studied the interplay between an X-ray photoevaporative disc wind and the sub-structures generated due to planet-disc interactions to determine how their mutual interactions affect the disc's and the planet's evolution.
Methods: We performed 3D hydrodynamic simulations of viscous discs (α = 6.9 × 10⁻⁴) that host a Jupiter-like planet and undergo X-ray photoevaporation. We traced the gas flows within the disc and wind and measured the rate of accretion onto the planet, as well as the gravitational torque that is acting on it.
Results: Our results show that the planetary gap removes the wind's pressure support, allowing wind material to fall back into the gap. This opens new pathways for material from the inner disc (and part of the outer disc) to be redistributed through the wind towards the gap. Consequently, the gap becomes shallower and the flow of mass across the gap in both directions is significantly increased, as is the planet's mass-accretion rate (by factors of ≈5 and ≈2, respectively). Moreover, the wind-driven redistribution results in a denser inner disc and a less dense outer disc, which, combined with the recycling of a significant portion of the inner wind, leads to longer lifetimes for the inner disc, contrary to the expectation in a planet-induced photoevaporation scenario that has been proposed in the past.

The SRG/eROSITA All-Sky Survey (eRASS) is expected to contain ~100 quasars that emitted their light when the universe was less than a billion years old, i.e. at z>5.6. By selection, these quasars populate the bright end of the AGN X-ray luminosity function and their count offers a powerful demographic diagnostic of the parent super-massive black hole population. Of the >~ 400 quasars that have been discovered at z>5.6 to date, less than 15 % have been X-ray detected. We present a pilot survey to uncover the elusive X-ray luminous end of the distant quasar population. We have designed a quasar selection pipeline based on optical, infrared and X-ray imaging data from DES DR2, VHS DR5, CatWISE2020 and the eRASS. The core selection method relies on SED template fitting. We performed optical follow-up spectroscopy with the Magellan/LDSS3 instrument for the redshift confirmation of a subset of candidates. We have further obtained a deeper X-ray image of one of our candidates with Chandra ACIS-S. We report the discovery of five new quasars in the redshift range 5.6 < z < 6.1. Two of these quasars are detected in eRASS and are by selection X-ray ultra-luminous. These quasars are also detected at radio frequencies. The first one is a broad absorption line quasar which shows significant X-ray dimming over 3.5 years, i.e. about 6 months in the quasar rest frame. The second radio-detected quasar is a jetted source with compact morphology. We show that a blazar configuration is likely for this source, making it the second most distant blazar known to date. With our pilot study, we demonstrate the power of eROSITA as a discovery machine for luminous quasars in the epoch of reionization. The X-ray emission of the two eROSITA detected quasars are likely to be driven by different high-energetic emission mechanisms a diversity which will be further explored in a future systematic full-hemisphere survey.

Context. Radio jets are present in a diverse sample of AGN. However, the mechanisms of jet powering are not fully understood, and it remains unclear to what extent they obey mass-invariant scaling relations similar to those found for the triggering and fuelling of X-ray-selected AGN.
Aims: We use the multi-wavelength data in the eFEDS field observed by eROSITA/Spectrum-Roentgen-Gamma (SRG) and LOFAR to study the incidence of X-ray and radio AGN as a function of several stellar mass (M_*)-normalised AGN power indicators.
Methods: From the LOFAR - eFEDS survey, we defined a new sample of radio AGN, with optical counterparts from Legacy Survey DR9, according to a radio-excess relative to their host star formation rate. We further divided the sample into compact and complex radio morphologies. In this work, we used the subset matching to the well-characterised, highly complete spectroscopic GAMA09 galaxies (0 < z < 0.4). We release this value-added LOFAR - eFEDS catalogue^*. We calculated the fraction of GAMA09 galaxies hosting radio, X-ray, and both radio and X-ray AGN as functions of the specific black hole kinetic (λ_Jet) and radiative (λ_Edd) power.
Results: Despite the soft-X-ray eROSITA-selected sample, the incidence of X-ray AGN as a function of λ_Edd shows the same mass-invariance and power law slope (−0.65) as that found in previous studies once corrected for completeness. Across the M_* range probed, the incidence of compact radio AGN as a function of λ_Jet is described by a power law with constant slope, showing that it is not only high mass galaxies hosting high power jets and vice versa. This slope is steeper than that of the X-ray incidence, which has a value of around −1.5. Furthermore, higher-mass galaxies are more likely to host radio AGN across the λ_Jet range, indicating some residual mass dependence of jet powering. Upon adding complex radio morphologies, including 34 FRIIs, three of which are giant radio galaxies, the incidence not only shows a larger mass dependence but also a jet power dependence, being clearly boosted at high λ_Jet values. Importantly, the latter effect cannot be explained by such radio AGN residing in more dense environments (or more massive dark matter haloes). The similarity in the incidence of quiescent and star-forming radio AGN reveals that radio AGN are not only found in "red and dead" galaxies. Overall, our incidence analysis reveals some fundamental statistical properties of radio AGN samples, but highlights open questions regarding the use of a single radio luminosity-jet power conversion. We explore how different mass and accretion rate dependencies of the incidence can explain the observed results for varying disk-jet coupling models.

The source catalogue is available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (ftp://130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/686/A43 or on the LOFAR Surveys DR website: https://lofar-surveys.org/efeds.html

The exploration of outflows in protobinary systems presents a challenging yet crucial endeavour, offering valuable insights into the dynamic interplay between protostars and their evolution. In this study, we examine the morphology and dynamics of jets and outflows within the IRAS 4A protobinary system. This analysis is based on ALMA observations of SiO(5-4), H₂CO(3_{0, 3}-2_{0, 3}), and HDCO(4_{1, 4}-3_{1, 3}) with a spatial resolution of ~150 au. Leveraging an astrochemical approach involving the use of diverse tracers beyond traditional ones has enabled the identification of novel features and a comprehensive understanding of the broader outflow dynamics. Our analysis reveals the presence of two jets in the redshifted emission, emanating from IRAS 4A1 and IRAS 4A2, respectively. Furthermore, we identify four distinct outflows in the region for the first time, with each protostar, 4A1 and 4A2, contributing to two of them. We characterize the morphology and orientation of each outflow, challenging previous suggestions of bends in their trajectories. The outflow cavities of IRAS 4A1 exhibit extensions of 10 and 13 arcsec with position angles (PA) of 0° and -12°, respectively, while those of IRAS 4A2 are more extended, spanning 18 and 25 arcsec with PAs of 29° and 26°. We propose that the misalignment of the cavities is due to a jet precession in each protostar, a notion supported by the observation that the more extended cavities of the same source exhibit lower velocities, indicating they may stem from older ejection events.

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

We estimate the efficiency of mitigating the lensing B-mode polarization, the so-called delensing, for the LiteBIRD experiment with multiple external data sets of lensing-mass tracers. The current best bound on the tensor-to-scalar ratio, r, is limited by lensing rather than Galactic foregrounds. Delensing will be a critical step to improve sensitivity to r as measurements of r become more and more limited by lensing. In this paper, we extend the analysis of the recent LiteBIRD forecast paper to include multiple mass tracers, i.e., the CMB lensing maps from LiteBIRD and CMB-S4-like experiment, cosmic infrared background, and galaxy number density from Euclid- and LSST-like survey. We find that multi-tracer delensing will further improve the constraint on r by about 20%. In LiteBIRD, the residual Galactic foregrounds also significantly contribute to uncertainties of the B-modes, and delensing becomes more important if the residual foregrounds are further reduced by an improved component separation method.

We present initial results from the Dark Energy Spectroscopic Instrument (DESI) complete calibration of the colour-redshift relation (DC3R2) secondary target survey. Our analysis uses 230 k galaxies that overlap with KiDS-VIKING ugriZYJHK_s photometry to calibrate the colour-redshift relation and to inform photometric redshift (photo-z) inference methods of future weak lensing surveys. Together with emission line galaxies (ELGs), luminous red galaxies (LRGs), and the Bright Galaxy Survey (BGS) that provide samples of complementary colour, the DC3R2 targets help DESI to span 56 per cent of the colour space visible to Euclid and LSST with high confidence spectroscopic redshifts. The effects of spectroscopic completeness and quality are explored, as well as systematic uncertainties introduced with the use of common Self-Organizing Maps trained on different photometry than the analysis sample. We further examine the dependence of redshift on magnitude at fixed colour, important for the use of bright galaxy spectra to calibrate redshifts in a fainter photometric galaxy sample. We find that noise in the KiDS-VIKING photometry introduces a dominant, apparent magnitude dependence of redshift at fixed colour, which indicates a need for carefully chosen deep drilling fields, and survey simulation to model this effect for future weak lensing surveys.

We study the possibility of using the LiteBIRD satellite B-mode survey to constrain models of inflation producing specific features in CMB angular power spectra. We explore a particular model example, i.e. spectator axion-SU(2) gauge field inflation. This model can source parity-violating gravitational waves from the amplification of gauge field fluctuations driven by a pseudoscalar "axionlike" field, rolling for a few e-folds during inflation. The sourced gravitational waves can exceed the vacuum contribution at reionization bump scales by about an order of magnitude and can be comparable to the vacuum contribution at recombination bump scales. We argue that a satellite mission with full sky coverage and access to the reionization bump scales is necessary to understand the origin of the primordial gravitational wave signal and distinguish among two production mechanisms: quantum vacuum fluctuations of spacetime and matter sources during inflation. We present the expected constraints on model parameters from LiteBIRD satellite simulations, which complement and expand previous studies in the literature. We find that LiteBIRD will be able to exclude with high significance standard single-field slow-roll models, such as the Starobinsky model, if the true model is the axion-SU(2) model with a feature at CMB scales. We further investigate the possibility of using the parity-violating signature of the model, such as the TB and EB angular power spectra, to disentangle it from the standard single-field slow-roll scenario. We find that most of the discriminating power of LiteBIRD will reside in BB angular power spectra rather than in TB and EB correlations.

We explore the capability of measuring lensing signals in LiteBIRD full-sky polarization maps. With a 30 arcmin beam width and an impressively low polarization noise of 2.16 μK-arcmin, LiteBIRD will be able to measure the full-sky polarization of the cosmic microwave background (CMB) very precisely. This unique sensitivity also enables the reconstruction of a nearly full-sky lensing map using only polarization data, even considering its limited capability to capture small-scale CMB anisotropies. In this paper, we investigate the ability to construct a full-sky lensing measurement in the presence of Galactic foregrounds, finding that several possible biases from Galactic foregrounds should be negligible after component separation by harmonic-space internal linear combination. We find that the signal-to-noise ratio of the lensing is approximately 40 using only polarization data measured over 80% of the sky. This achievement is comparable to Planck's recent lensing measurement with both temperature and polarization and represents a four-fold improvement over Planck's polarization-only lensing measurement. The LiteBIRD lensing map will complement the Planck lensing map and provide several opportunities for cross-correlation science, especially in the northern hemisphere.

In this work, we propose a novel approach to the homotopy transfer procedure starting from a set of homotopy data such that the first differential complex is a differential graded module over the second one. We show that the module structure may be used to induce an $A_\infty$-algebra on the second differential complex, constructed in a similar fashion to the homotopy transfer $A_\infty$-algebra. We prove that, under certain conditions, the $A_\infty$-algebras obtained with this procedure are quasi-isomorphic to the homotopy transfer one. On the other hand, when the side conditions do not hold, we find that there are cases where the existence of an $A_\infty$-quasi-isomorphism with the homotopy transfer $A_\infty$-algebra is obstructed. In other words, we obtain a new $A_\infty$-algebra on the second complex, inequivalent to the homotopy transfer one. Lastly, we prove that these $A_\infty$-algebras are not infinitesimal Hochschild deformations of the homotopy transfer $A_\infty$-algebra.

Simulating high-resolution detector responses is a computationally intensive process that has long been challenging in Particle Physics. Despite the ability of generative models to streamline it, full ultra-high-granularity detector simulation still proves to be difficult as it contains correlated and fine-grained information. To overcome these limitations, we propose Intra-Event Aware Generative Adversarial Network (IEA-GAN). IEA-GAN presents a Transformer-based Relational Reasoning Module that approximates an event in detector simulation, generating contextualized high-resolution full detector responses with a proper relational inductive bias. IEA-GAN also introduces a Self-Supervised intra-event aware loss and Uniformity loss, significantly enhancing sample fidelity and diversity. We demonstrate IEA-GAN's application in generating sensor-dependent images for the ultra-high-granularity Pixel Vertex Detector (PXD), with more than 7.5 M information channels at the Belle II Experiment. Applications of this work span from Foundation Models for high-granularity detector simulation, such as at the HL-LHC (High Luminosity LHC), to simulation-based inference and fine-grained density estimation.

Testing gravity and the concordance model of cosmology, $\Lambda$CDM, at large scales is a key goal of this decade's largest galaxy surveys. Here we present a comparative study of dark matter power spectrum predictions from different numerical codes in the context of three popular theories of gravity that induce scale-independent modifications to the linear growth of structure: nDGP, Cubic Galileon and K-mouflage. In particular, we compare the predictions from full $N$-body simulations, two $N$-body codes with approximate time integration schemes, a parametrised modified $N$-body implementation and the analytic halo model reaction approach. We find the modification to the $\Lambda$CDM spectrum is in $2\%$ agreement for $z\leq1$ and $k\leq 1~h/{\rm Mpc}$ over all gravitational models and codes, in accordance with many previous studies, indicating these modelling approaches are robust enough to be used in forthcoming survey analyses under appropriate scale cuts. We further make public the new code implementations presented, specifically the halo model reaction K-mouflage implementation and the relativistic Cubic Galileon implementation.

Axion quark nuggets (AQNs) are hypothetical objects with a mass greater than a few grams and sub-micrometer size, formed during the quark-hadron transition. Originating from the axion field, they offer a possible resolution of the similarity between visible and dark components of the Universe. These composite objects behave as cold dark matter, interacting with ordinary matter and resulting in pervasive electromagnetic radiation throughout the Universe. This work aims to predict the electromagnetic signature in large-scale structures from the AQN-baryon interaction, accounting for thermal and non-thermal radiations. We use Magneticum hydrodynamical simulations to describe the distribution and dynamics of gas and dark matter at cosmological scales. We calculate the electromagnetic signature from radio, starting at $\nu \sim$ 1 GHz, up to a few keV X-ray energies. We find that the AQNs signature is characterized by monopole and fluctuation signals. The amplitude of both signals strongly depends on the average AQN mass and the ionization level of the baryonic environment. We identify a most optimistic scenario with a signal often near the sensitivity limit of existing instruments, such as FIRAS and the South Pole Telescope for high-resolution. Fluctuations in the Extra-galactic Background Light caused by the AQN can be tested with space-based imagers Euclid and James Webb Space Telescope. We also identify a minimal configuration, still out of reach of existing instruments, but future experiments might be able to pose constraints on the AQN model. We conclude that this is a viable dark matter model, which does not violate the canons of cosmology, nor existing observations. The best chances for testing this model reside in 1) ultra-deep IR and optical surveys, 2) spectral distorsions of the CMB and 3) low-frequency (1 GHz < $\nu $ < 100 GHz) and high-resolution ($\ell > 10^4$) observations.

We use the reconstructed properties of the Amaterasu particle, the second-highest energy cosmic ray ever detected, to map out three-dimensional constraints on the location of its unknown source. We highlight possible astrophysical sources that are compatible with these regions and requirements. Among these, M82, a powerful starburst galaxy, stands out as a strong candidate due to its position and proximity. To derive our constraints, we use CRPropa 3 to model all relevant propagation effects, including deflections in the Galactic and extra-Galactic magnetic fields. We consider key input quantities such as source distance, position, energy, and the strength and coherence length of the extra-Galactic magnetic field as free parameters. We then infer constraints on these parameters by applying approximate Bayesian computation. We present our results, demonstrating the impact of different assumptions for the arrival mass of the Amaterasu particle and the systematic uncertainties on the energy scale.

Context. In this work, we test feedback propagation models on the test case of 2MASS 0918+2117 (2M0918), a z = 0.149 X-ray variable AGN that shows tentative evidence for nuclear ultra-fast outflows (UFOs) in a 2005 XMM-Newton observation. We also investigate whether UFOs can be related to the observed X-ray variability.
Aims: We observed 2M0918 with XMM-Newton and NuSTAR in 2020 to confirm the presence of and characterize the UFOs. We performed a kinematic analysis of the publicly available 2005 SDSS optical spectrum to reveal and measure the properties of galaxy-scale ionized outflows. Furthermore, we constructed 20-year-long light curves of observed flux, line-of-sight column density, and intrinsic accretion rate from the spectra of the first four SRG/eROSITA all-sky surveys and archival observations from Chandra and XMM-Newton.
Methods: We detect UFOs with v ∼ 0.16c and galaxy-scale ionized outflows with velocities of ∼700 km s⁻¹. We also find that the drastic X-ray variability (factors > 10) can be explained in terms of variable obscuration and variable intrinsic luminosity.
Results: Comparing the energetics of the two outflow phases, 2M0918 is consistent with momentum-driven wind propagation. 2M0918 expands the sample of AGN with both UFOs and ionized gas winds from 5 to 6 and brings the sample of AGN hosting multiscale outflows to 19, contributing to a clearer picture of feedback physics. From the variations in accretion rate, column density, and ionization level of the obscuring medium, we propose a scenario that connects obscurers, an accretion enhancement, and the emergence of UFOs.

Global magnetic fields of early-type stars are commonly characterised by the mean longitudinal magnetic field <B_z> and the mean field modulus , derived from the circular polarisation and intensity spectra, respectively. Observational studies often report a root mean square (rms) of <B_z> and an average value of . In this work, I used numerical simulations to establish statistical relationships between these cumulative magnetic observables and the polar strength, B_d, of a dipolar magnetic field. I show that in the limit of many measurements randomly distributed in rotational phase, <B_z>_rms = 0.179_−0.043^+0.031 B_d and _avg = 0.691_−0.023^+0.020 B_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.

Context. The rapidly improving quality and resolution of both low surface brightness observations and cosmological simulations of galaxies enable us to address the important question of how the formation history is imprinted in the outer unrelaxed regions of galaxies, and to inspect the correlations of these imprints with another tracer of galaxy formation, the internal kinematics.
Aims: Using the hydrodynamical cosmological simulation called Magneticum Pathfinder, we identified tidal tails, shells, streams, and satellite planes, and connected them to the amount of rotational support and the formation histories of the host galaxies. This presents the first combined statistical census considering all these four types of features in hydrodynamical cosmological simulations.
Methods: Tidal features were visually classified from a three-dimensional rendering of the simulated galaxies by several scientists independently. Only features that were identified by at least half of the participating individuals were considered to be existing features. The data on satellite planes and kinematic properties of the simulated galaxies were taken from previous work. The results were compared to observations, especially from the MATLAS survey.
Results: Generally, prominent features are much more common around elliptical than around disk galaxies. Shells are preferentially found around kinematically slowly rotating galaxies in both simulations and observations, while streams can be found around all types of galaxies, with a slightly higher probability to be present around less rotationally supported galaxies. Tails and satellite planes, however, appear independently of the internal kinematics of the central galaxy, indicating that they are formed through processes that have not (yet) affected the internal kinematics. Prolate rotators have the overall highest probability to exhibit tidal features, but the highest likelihood for a specific type of feature is found for galaxies with kinematically distinct cores (KDCs), nearly 20% of which exhibit streams.
Conclusions: As shells are formed through radial merger events while streams are remnants of circular merger infall, this suggests that the orbital angular momentum of the merger event plays a more crucial role in transforming the host galaxy than previously anticipated. The existence of a shell around a given slow rotator furthermore is a sign of a radial merger formation for this particular slow rotator because one-third of the galaxies with a shell were transformed into slow rotators by the merger event that also caused the shells. The appearance of a stream around a KDC is a direct indicator for the multiple merger formation pathway of that KDC as opposed to the major merger pathway.

Context. Recent observations of protostellar cores suggest that most of the material in the protostellar phase is accreted along streamers. Streamers in this context are defined as velocity coherent funnels of denser material potentially connecting the large-scale environment to the small scales of the forming accretion disk.
Aims: Using simulations that simultaneously resolve the driving of turbulence on the filament scale as well as the collapse of the core down to protostellar disk scales, we aim to understand the effect of the turbulent velocity field on the formation of overdensities in the accretion flow.
Methods: We performed a three-dimensional numerical study on a core collapse within a turbulent filament using the RAMSES code and analysed the properties of overdensities in the accretion flow.
Results: We find that overdensities are formed naturally by the initial turbulent velocity field inherited from the filament and subsequent gravitational collimation. This leads to streams that are not really filamentary but show a sheet-like morphology. Moreover, they have the same radial infall velocities as the low density material. As a main consequence of the turbulent initial condition, the mass accretion onto the disk does not follow the predictions for solid body rotation. Instead, most of the mass is funneled by the overdensities to intermediate disk radii.

We developed a grid of stellar rotation models for low-mass and solar-type classical T Tauri stars (CTTS; 0.3M _⊙ < M _* < 1.2M _⊙). These models incorporate the star–disk interaction and magnetospheric ejections to investigate the evolution of the stellar rotation rate as a function of the mass of the star M _*, the magnetic field (B _*), and stellar wind ( ). We compiled and determined stellar parameters for 208 CTTS, such as projected rotational velocity , mass accretion rate , stellar mass M _*, ages, and estimated rotational periods using Transiting Exoplanet Survey Satellite (TESS) data. We also estimated a representative value of the mass-loss rate for our sample using the [O I] λ 6300 spectral line. Our results confirm that measurements in CTTS agree with the rotation rates provided by our spin models in the accretion-powered stellar winds picture. In addition, we used the approximate Bayesian computation technique to explore the connection between the model parameters and the observational properties of CTTS. We find that the evolution of with age might be regulated by variations in (1) the intensity of B _* and (2) the fraction of the accretion flow ejected in magnetic winds, removing angular momentum from these systems. The youngest stars in our sample (∼1 Myr) show a median branching ratio and median B _* ∼ 2000 G, in contrast to ∼0.01 and 1000 G, respectively, for stars with ages ≳3 Myr.

We study two-loop corrections to the scattering amplitude of four massive leptons in quantum electrodynamics. These amplitudes involve previously unknown elliptic Feynman integrals, which we compute analytically using the differential equation method. In doing so, we uncover the details of the elliptic geometry underlying this scattering amplitude and show how to exploit its properties to obtain compact, easy-to-evaluate series expansions that describe the scattering of four massive leptons in QED in the kinematical regions relevant for Bhabha and Møller scattering processes.

In this work we demonstrate that significant gains in performance and data efficiency can be achieved in High Energy Physics (HEP) by moving beyond the standard paradigm of sequential optimization or reconstruction and analysis components. We conceptually connect HEP reconstruction and analysis to modern machine learning workflows such as pretraining, finetuning, domain adaptation and high-dimensional embedding spaces and quantify the gains in the example usecase of searches of heavy resonances decaying via an intermediate di-Higgs system to four b-jets. To our knowledge this is the first example of a low-level feature extraction network finetuned for a downstream HEP analysis objective.

Context. There are several symbiotic stars (e.g., BF Cyg, Z And, and FN Sgr) in which periodic signals of tens of minutes have been detected. These periods have been interpreted as the spin period of magnetic white dwarfs that accrete through a magnetic stream originating from a truncated accretion disc.
Aims: To shed light on the origin of magnetic symbiotic stars, we investigated the system FN Sgr in detail. We searched for a reasonable formation pathway to explain its stellar and binary parameters including the magnetic field of the accreting white dwarf.
Methods: We used the MESA code to carry out pre-CE and post-CE binary evolution and determined the outcome of CE evolution assuming the energy formalism. For the origin and evolution of the white dwarf magnetic field, we adopted the crystallization scenario.
Results: We found that FN Sgr can be explained as follows. First, a non-magnetic white dwarf is formed through CE evolution. Later, during post-CE evolution, the white dwarf starts to crystallize and a weak magnetic field is generated. After a few hundred million years, the magnetic field penetrates the white dwarf surface and becomes detectable. Meanwhile, its companion evolves and becomes an evolved red giant. Subsequently, the white dwarf accretes part of the angular momentum from the red giant stellar winds. As a result, the white dwarf spin period decreases and its magnetic field reaches super-equipartition, getting amplified due to a rotation- and crystallization-driven dynamo. The binary then evolves into a symbiotic star, with a magnetic white dwarf accreting from an evolved red giant through atmospheric Roche-lobe overflow.
Conclusions: We conclude that the rotation- and crystallization-driven dynamo scenario, or any age-dependent scenario, can explain the origin of magnetic symbiotic stars reasonably well. This adds another piece to the pile of evidence supporting this scenario. If our formation channel is correct, our findings suggest that white dwarfs in most symbiotic stars formed through CE evolution might be magnetic, provided that the red giant has spent ≳3 Gyr as a main-sequence star.

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.

We consider the strong coupling limit of lattice QCD with massless staggered quarks and study the resource requirements for quantum simulating the theory in its Hamiltonian formulation. The bosonic Hilbert space of the color-singlet degrees of freedom grows quickly with the number of quark flavors, making it a suitable testing ground for resource considerations across different platforms. In particular, in addition to the standard model of computation with qubits, we consider mapping the theory to qudits $(d>2)$ and qumodes, as used on trapped-ion systems and photonic devices, respectively.

Active and sterile neutrinos may acquire "dark moments" via one-loop diagrams with a massless dark photon and new light particles in the loop. Due to the kinetic mixing between the dark photon and the Standard Model photon, neutrinos would obtain effective electromagnetic moments. This mechanism allows for enhanced electromagnetic moments that can evade constraints on particles directly charged under electromagnetism. We show that in a wide region of parameter space, the model features testable predictions for the anapole and magnetic moment of active and sterile neutrinos with dark matter direct detection experiments sensitive to the solar neutrino flux.

The mass loss mechanism of red supergiant stars is not well understood, even though it has crucial consequences for their stellar evolution and the appearance of supernovae that occur upon core-collapse. We argue that outgoing shock waves launched near the photosphere can support a dense chromosphere between the star's surface and the dust formation radius at several stellar radii. We derive analytic expressions for the time-averaged density profile of the chromosphere, and we use these to estimate mass loss rates due to winds launched by radiation pressure at the dust formation radius. These mass loss rates are similar to recent observations, possibly explaining the upward kink in mass loss rates of luminous red supergiants. Our models predict that low-mass red supergiants lose less mass than commonly assumed, while high-mass red supergiants lose more. The chromospheric mass of our models is $$0.01 solar masses, most of which lies within a few stellar radii. This can help explain the early light curves and spectra of type-II P supernovae without requiring extreme pre-supernova mass loss. We discuss implications for stellar evolution, type II-P supernovae, SN 2023ixf, and Betelgeuse.

We compute the suppression of bottomonium in the quark-gluon plasma using the three-loop QCD static potential. The potential describes the spin-averaged bottomonium spectrum below threshold with a less than 1% error. Within potential nonrelativistic quantum chromodynamics and an open quantum systems framework, we compute the evolution of the bottomonium density matrix. The values of the quarkonium transport coefficients are obtained from lattice QCD measurements of the bottomonium in-medium width and thermal mass shift; we additionally include for the first time a vacuum contribution to the dispersive coefficient γ . Using the three-loop potential and the values of the heavy quarkonium transport coefficients, we find that the resulting bottomonium nuclear modification factor is consistent with experimental observations, while at the same time reproducing the lattice measurements of the in-medium width.

The static QCD force from the lattice can be used to extract Λ_{MS ¯}, which determines the running of the strong coupling. Usually, this is done with a numerical derivative of the static potential. However, this introduces additional systematic uncertainties; thus, we use another observable to measure the static force directly. This observable consists of a Wilson loop with a chromoelectric field insertion. We work in the pure SU(3) gauge theory. We use gradient flow to improve the signal-to-noise ratio and to address the field insertion. We extract Λ_MS^{¯ n_f=0} from the data by exploring different methods to perform the zero-flow-time limit. We obtain the value √{8 t₀ }Λ_MS^{¯ n_f=0}=0.62 9_-26⁺²² , where t₀ is a flow-time reference scale. We also obtain precise determinations of several scales: r₀/r₁, √{8 t₀ }/r₀, √{8 t₀ }/r₁, and we compare these to the literature. The gradient flow appears to be a promising method for calculations of Wilson loops with chromoelectric and chromomagnetic insertions in quenched and unquenched configurations.

Off-lightcone Wilson-line operators are constructed using local operators connected by time-like or space-like Wilson lines, which ensure gauge invariance. Off-lightcone Wilson-line operators have broad applications in various contexts. For instance, space-like Wilson-line operators play a crucial role in determining quasi-distribution functions (quasi-PDFs), while time-like Wilson-line operators are essential for understanding quarkonium decay and production within the potential non-relativistic QCD (pNRQCD) framework. In this work, we establish a systematic approach for calculating the matching from the gradient-flow scheme to the MS ¯ scheme in the limit of small flow time for off-lightcone Wilson-line operators. By employing the one-dimensional auxiliary-field formalism, we simplify the matching procedure, reducing it to the matching of local current operators. We provide one-loop level matching coefficients for these local current operators. For the case of hadronic matrix element related to the quark quasi-PDFs, we show at one-loop level that the finite flow time effect is very small as long as the flow radius is smaller than the physical distance z, which is usually satisfied in lattice gradient flow computations. Applications include lattice gradient flow computations of quark/gluon quasi-PDFs, gluonic correlators related to quarkonium decay and production in pNRQCD, and spin-dependent potentials in terms of chromoelectric and chromomagnetic field insertions into a Wilson loop.

We scrutinize the Sommerfeld enhancement in dark matter pair annihilation for p-wave and higher-ℓ partial waves. For the Yukawa potential these feature a super-resonant Breit-Wigner peak in their velocity-dependence close to Sommerfeld resonances as well as a universal scaling with velocity for all ℓ ≥ 1 that differs from the s-wave case. We provide a quantum mechanical explanation for these phenomena in terms of quasi-bound states sustained by the centrifugal barrier of the partial-wave potential, and give approximate WKB expressions capturing the main effects. The impact of quasi-bound states is exemplified for wino dark matter and models with light mediators, with a focus on indirect detection signals. We note that quasi-bound states can also explain similar peaks in the bound-state formation and self-scattering cross sections.

The reduction and distortion of quantum correlations in the presence of classical noise leads to varied levels of inefficiency in the availability of entanglement as a resource for quantum information processing protocols. While generically minimizing required entanglement for mixed quantum states remains challenging, a class of many-body Gaussian quantum states (denoted N IC ) is here identified that exhibits two-mode bipartite entanglement structure, resembling that of pure states, for which the logarithmic negativity entanglement measure remains invariant upon inclusion of the classical correlations and optimal entanglement resources can be clearly quantified. This subclass is found to be embedded within a broader class of many-body Gaussian states (denoted N -SOL) that retain two-mode entanglement structure for detection processes. These two entanglement classes are relevant in theoretical and experimental applications from the scalar field vacuum to the local axial motional modes of trapped ion chains. Utilizing the subspace that heralds inseparability in response to partial transposition, a minimum noise filtering process is designed to be necessary, sufficient, and computable for determining membership in these classes of entanglement structure. Application of this process to spacelike regions of the free scalar field vacuum is found to improve resource upper bounds, providing new understanding of the entanglement required for the quantum simulation of quantum fields as observed by arrays of local detectors.

We study the evolution of the two-pion correlation function parameters with collision energy in the context of relativistic heavy-ion collisions within the NICA energy range. To this end, we perform UrQMD simulations in the cascade mode to produce samples of pions from <inline-formula id="IEq1"><mml:math><mml:mrow><mml:mn>5</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn>10</mml:mn><mml:mn>6</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> Bi+Bi collisions for each of the studied energies. The effects of the quantum-statistical correlations are introduced using the correlation afterburner code CRAB. We fit the correlation function using Gaussian, exponential and symmetric Lévy shapes and show that for all collision energies the latter provides the best fit. We separate the sample into pions coming from primary processes and pions originating from the decay of long-lived resonances, and show that the source size for the latter is significantly larger than for the former. The source size for the secondaries, is similar but in general larger than the size for the whole pion sample. To further characterize the pion source, we also simulate the effects of a non-ideal detector introducing a momentum smearing parameter, representing the minimum pair momentum and thus a maximum source size that can be resolved. The values of the correlation function intercept parameter are therefore modified from the values they attain for the perfect detector case. Using the core-halo picture of the source, we show that the values of the intercept parameter are influenced by the presence of a significant fraction of core pions coming from the decay of long-lived but slow-moving resonances. These findings serve as a benchmark to compare with future Monte Carlo studies that consider an Equation of State and thus allow for a phase transition within the studied energy domain.

During the public Kaggle competition "IceCube – Neutrinos in Deep Ice", thousands of reconstruction algorithms were created and submitted, aiming to estimate the direction of neutrino events recorded by the IceCube detector. Here we describe in detail the three ultimate best, award-winning solutions. The data handling, architecture, and training process of each of these machine learning models is laid out, followed up by an in-depth comparison of the performance on the Kaggle datatset. We show that on cascade events in IceCube above 10 TeV, the best Kaggle solution is able to achieve an angular resolution of better than 5<inline-formula id="IEq1"><mml:math><mml:mmultiscripts><mml:mrow></mml:mrow><mml:mrow></mml:mrow><mml:mo>∘</mml:mo></mml:mmultiscripts></mml:math></inline-formula>, and for tracks correspondingly better than 0.5<inline-formula id="IEq2"><mml:math><mml:mmultiscripts><mml:mrow></mml:mrow><mml:mrow></mml:mrow><mml:mo>∘</mml:mo></mml:mmultiscripts></mml:math></inline-formula>. These results indicate that the Kaggle solutions perform at a level comparable to the current state-of-the-art in the field, and that they may even be able to outperform existing reconstruction resolutions for certain types of events.

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

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.

Context. It has been claimed for more than a decade that energies other than orbital and thermodynamic internal are required to explain post-common envelope (CE) binaries with sufficiently long orbital periods (≳1 d) hosting AFGK-type main-sequence stars (∼0.5 − 2.0 M_⊙) paired with oxygen-neon white dwarfs (≳1.1 M_⊙). This would imply a completely different energy budget during CE evolution for these post-CE binaries in comparison to the remaining systems hosting M dwarfs and/or less massive white dwarfs.
Aims: In this first in a series of papers related to long-period post-CE binaries, we investigated whether extra energy is required to explain the currently known post-CE binaries with sufficiently long orbital periods consisting of oxygen-neon white dwarfs with AFGK-type main-sequence star companions.
Methods: We carried out binary population simulations with the BSE code adopting empirically derived inter-correlated main-sequence binary distributions for the initial binary population and assuming that the only energy, in addition to orbital, that help to unbind the CE is thermal energy. We also searched for the formation pathways of the currently known systems from the zero-age main-sequence binary to their present-day observed properties.
Results: Unlike what has been claimed for a long time, we show that all such post-CE binaries can be explained by assuming inefficient CE evolution, which is consistent with results achieved for the remaining post-CE binaries. There is therefore no need for an extra energy source. We also found that for CE efficiency close to 100%, post-CE binaries hosting oxygen-neon white dwarfs with orbital periods as long as one thousand days can be explained. For all known systems we found formation pathways consisting of CE evolution triggered when a highly evolved (i.e. when the envelope mass is comparable to the core mass), thermally pulsing, asymptotic giant branch star fills its Roche lobe at an orbital period of several thousand days. Due to the sufficiently low envelope mass and sufficiently long orbital period, the resulting post-CE orbital period can easily be several tens of days.
Conclusions: We conclude that the known post-CE binaries with oxygen-neon white dwarfs and AFGK-type main-sequence stars can be explained without invoking any energy source other than orbital and thermal energy. Our results strengthen the idea that the most common formation pathway of the overall population of post-CE binaries hosting white dwarfs is through inefficient CE evolution.

We present the results of a search for Miras and long-period variables (LPVs) in M33 using griJHK_S archival observations from the Canada-France-Hawaii Telescope. We use multiband information and machine learning techniques to identify and characterize these variables. We recover ~1300 previously discovered Mira candidates and identify ~13 000 new Miras and LPVs. We detect for the first time a clear first-overtone pulsation sequence among Mira candidates in this galaxy. We use O-rich, fundamental-mode Miras in the LMC and M33 to derive a distance modulus for the latter of μ = 24.629 ± 0.046 mag.

The inverted pendulum is a mechanical system with a rapidly oscillating pivot point. Using techniques similar in spirit to the methodology of effective field theories, we derive an effective Lagrangian that allows for the systematic computation of corrections to the so-called Kapitza equation. The derivation of the effective potential of the system requires non-trivial matching conditions, which need to be determined order by order in the power-counting of the problem. The convergence behavior of the series is investigated on the basis of high-order results obtained by this method. The results from this analysis can be used to determine the regions of parameter space, in which the inverted position of the pendulum is stable or unstable to high precision.

Within the framework of non-relativistic QCD (NRQCD) effective field theory, we study the leptoproduction of $J/{\psi}$ at next-to-leading order in perturbative QCD for both unpolarized and polarized electron-ion collisions. We demonstrate that the $J/{\psi}$-tagged deep inelastic scattering in the future Electron-Ion Collider can be served as a golden channel for the reasons including constraining NRQCD long distance matrix elements, probing the nuclear gluon distribution functions, as well as investigating the gluon helicity distribution inside a longitudinal polarized proton.

We present forecasts for constraints on the Hu and Sawicki f (R ) modified gravity model using realistic mock data representative of future cluster and weak lensing surveys. We create mock thermal Sunyaev-Zel'dovich effect selected cluster samples for SPT-3G and CMB-S4 and the corresponding weak gravitational lensing data from next-generation weak-lensing (ngWL) surveys like Euclid and Rubin. We employ a state-of-the-art Bayesian likelihood approach that includes all observational effects and systematic uncertainties to obtain constraints on the f (R ) gravity parameter log₁₀|f_{R 0}|. In this analysis we vary the cosmological parameters [Ω_m,Ω_νh²,h ,A_s,n_s,log₁₀|f_{R 0}|], which allows us to account for possible degeneracies between cosmological parameters and f (R ) modified gravity. The analysis accounts for f (R ) gravity via its effect on the halo mass function which is enhanced on cluster mass scales compared to the expectations within general relativity (GR). Assuming a fiducial GR model, the upcoming cluster dataset SPT -3 G ×ngWL is expected to obtain an upper limit of log₁₀|f_{R 0}|<-5.95 at 95% credibility, which significantly improves upon the current best bounds. The CMB -S 4 ×ngWL dataset is expected to improve this even further to log₁₀|f_{R 0}|<-6.23 . Furthermore, f (R ) gravity models with log₁₀|f_{R 0}|≥-6 , which have larger numbers of clusters, would be distinguishable from GR with both datasets. We also report degeneracies between log₁₀|f_{R 0}| and Ω_m as well as σ₈ for log₁₀|f_{R 0}|>-6 and log₁₀|f_{R 0}|>-5 respectively. Our forecasts indicate that future cluster abundance studies of f (R ) gravity will enable substantially improved constraints that are competitive with other cosmological probes.

A relativistic microscopic optical model potential for nucleon-nucleus scattering is developed based on the ab initio relativistic Brueckner-Hartree-Fock (RBHF) theory with the improved local density approximation, which is abbreviated as the RBOM potential. Both real and imaginary parts of the single-particle potentials in symmetric and asymmetric nuclear matter at various densities are determined uniquely in the full Dirac space. The density distributions of the target nuclei are calculated by the covariant energy density functional theory with the density functional PC-PK1. The central and spin-orbit terms of the optical potentials are quantitatively consistent with the relativistic phenomenological optical potentials. The performance of the RBOM potential is evaluated by considering proton scattering with incident energy E ≤200 MeV on five target nuclei, ²⁰⁸Pb, ¹²⁰Sn, ⁹⁰Zr, ⁴⁸Ca, and ⁴⁰Ca. Scattering observables including the elastic scattering angular distributions, analyzing powers, spin rotation functions, and reaction cross sections are analyzed. Theoretical predictions show good agreements with the experimental data and the results derived from phenomenological optical potentials. We anticipate that the RBOM potential can provide reference for other phenomenological and microscopic optical model potentials, as well as reliable descriptions for nucleon scattering on exotic nuclei in the era of rare-isotope beams.

We report on the SRG/eROSITA detection of ultra-soft ($kT=47^{+5}_{-5}$ eV) X-ray emission (L_X =$2.5^{+0.6}_{-0.5} \times 10^{43}$ erg s^-1) from the tidal disruption event (TDE) candidate AT 2022dsb ~14 d before peak optical brightness. As the optical luminosity increases after the eROSITA detection, then the 0.2-2 keV observed flux decays, decreasing by a factor of ~39 over the 19 d after the initial X-ray detection. Multi-epoch optical spectroscopic follow-up observations reveal transient broad Balmer emission lines and a broad He II 4686 Å emission complex with respect to the pre-outburst spectrum. Despite the early drop in the observed X-ray flux, the He II 4686 Å complex is still detected for ~40 d after the optical peak, suggesting the persistence of an obscured hard ionizing source in the system. Three outflow signatures are also detected at early times: (i) blueshifted H α emission lines in a pre-peak optical spectrum, (ii) transient radio emission, and (iii) blueshifted Ly α absorption lines. The joint evolution of this early-time X-ray emission, the He II 4686 Å complex, and these outflow signatures suggests that the X-ray emitting disc (formed promptly in this TDE) is still present after optical peak, but may have been enshrouded by optically thick debris, leading to the X-ray faintness in the months after the disruption. If the observed early-time properties in this TDE are not unique to this system, then other TDEs may also be X-ray bright at early times and become X-ray faint upon being veiled by debris launched shortly after the onset of circularization.

Charged and quasi-neutral beams propagating through an unmagnetised plasma are subject to numerous collisionless instabilities on the small scale of the plasma skin depth. The electrostatic two-stream instability, driven by longitudinal and transverse wakefields, dominates for dilute beams. This leads to modulation of the beam along the propagation direction and, for wide beams, transverse filamentation. A three-dimensional spatiotemporal two-stream theory for warm beams with a finite extent is developed. Unlike the cold beam limit, diffusion due to a finite emittance gives rise to a dominant wavenumber, and a cut-off wavenumber above which filamentation is suppressed. Particle-in-cell simulations give excellent agreement with the theoretical model. This work provides deeper insights into the effect of diffusion on filamentation of finite beams, crucial for comprehending plasma-based accelerators in laboratory and cosmic settings.

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.

Polarization of the cosmic microwave background (CMB) can help probe the fundamental physics behind cosmic inflation via the measurement of primordial B modes. As this requires exquisite control over instrumental systematics, some next-generation CMB experiments plan to use a rotating half-wave plate (HWP) as polarization modulator. However, the HWP non-idealities, if not properly treated in the analysis, can result in additional systematics. In this paper, we present a simple, semi-analytical end-to-end model to propagate the HWP non-idealities through the macro-steps that make up any CMB experiment (observation of multi-frequency maps, foreground cleaning, and power spectra estimation) and compute the HWP-induced bias on the estimated tensor-to-scalar ratio, r. We find that the effective polarization efficiency of the HWP suppresses the polarization signal, leading to an underestimation of r. Laboratory measurements of the properties of the HWP can be used to calibrate this effect, but we show how gain calibration of the CMB temperature can also be used to partially mitigate it. On the basis of our findings, we present a set of recommendations for the HWP design that can help maximize the benefits of gain calibration.

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.

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.

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.

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.