We present analysis of the plateau and late-time phase properties of a sample of 39 Type II supernovae (SNe II) that show narrow, transient, high-ionization emission lines (i.e., "IIn-like") in their early-time spectra from interaction with confined, dense circumstellar material (CSM). Originally presented by Jacobson-Galán et al 2024a, this sample also includes multicolor light curves and spectra extending to late-time phases of 35 SNe with no evidence for IIn-like features at <2 days after first light. We measure photospheric phase light-curve properties for the distance-corrected sample and find that SNe II with IIn-like features have significantly higher luminosities and decline rates at +50 days than the comparison sample, which could be connected to inflated progenitor radii, lower ejecta mass, and/or persistent CSM interaction. However, we find no statistical evidence that the measured plateau durations and $^{56}$Ni masses of SNe II with and without IIn-like features arise from different distributions. We estimate progenitor zero-age main sequence (ZAMS) masses for all SNe with nebular spectroscopy through spectral model comparisons and find that most objects, both with and without IIn-like features, are consistent with progenitor masses <12.5 M$_{\odot}$. Combining progenitor ZAMS masses with CSM densities inferred from early-time spectra suggests multiple channels for enhanced mass loss in the final years before core collapse such as a convection-driven chromosphere or binary interaction. Finally, we find spectroscopic evidence for ongoing ejecta-CSM interaction at radii $>10^{16}$ cm, consistent with substantial progenitor mass-loss rates of $\sim 10^{-4}$--$10^{-5}$ M$_{\odot}$ yr$^{-1}$ ($v_w < 50$ km/s) in the final centuries to millennia before explosion.

Effective Field Theory (EFT) modeling is expected to be a useful tool in the era of future higher-redshift galaxy surveys such as DESI-II and Spec-S5 due to its robust description of various large-scale structure tracers. However, large values of EFT bias parameters of higher-redshift galaxies could jeopardize the convergence of the perturbative expansion. In this paper we measure the bias parameters and other EFT coefficients from samples of two types of star-forming galaxies in the state-of-the-art MilleniumTNG and Astrid hydrodynamical simulations. Our measurements are based on the field-level EFT forward model that allows for precision EFT parameter measurements by virtue of cosmic variance cancellation. Specifically, we consider approximately representative samples of Lyman-break galaxies (LBGs) and Lyman-alpha emitters (LAEs) that are consistent with the observed (angular) clustering and number density of these galaxies at $z=3$. Reproducing the linear biases and number densities observed from existing LAE and LBG data, we find quadratic bias parameters that are roughly consistent with those predicted from the halo model coupled with a simple halo occupation distribution model. We also find non-perturbative velocity contributions (Fingers of God) of a similar size for LBGs to the familiar case of Luminous Red Galaxies. However, these contributions are quite small for LAEs despite their large satellite fraction values of up to $\sim 30\%$. Our results indicate that the effective momentum reach $k_{\rm{Max}}$ at $z=3$ for LAEs (LBGs) will be in the range $0.3-0.6 ~h\rm{Mpc}^{-1}$ ($0.2-0.8~h\rm{Mpc}^{-1}$), suggesting that EFT will perform well for high redshift galaxy clustering. This work provides the first step toward obtaining realistic simulation-based priors on EFT parameters for LAEs and LBGs.

Nuclear fission leads to the splitting of a nucleus into two fragments^1,2. Studying the distribution of the masses and charges of the fragments is essential for establishing the fission mechanisms and refining the theoretical models^3,4. It has value for our understanding of r-process nucleosynthesis^5,6, in which the fission of nuclei with extreme neutron-to-proton ratios is pivotal for determining astrophysical abundances and understanding the origin of the elements⁷ and for energy applications^8,9. Although the asymmetric distribution of fragments is well understood for actinides (elements in the periodic table with atomic numbers from 89 to 103) based on shell effects¹⁰, symmetric fission governs the scission process for lighter elements. However, unexpected asymmetric splits have been observed in neutron-deficient exotic nuclei¹¹, prompting extensive further investigations. Here we present measurements of the charge distributions of fission fragments for 100 exotic fissioning systems, 75 of which have never been measured, and establish a connection between the neutron-deficient sub-lead region and the well-understood actinide region. These new data comprehensively map the asymmetric fission island and provide clear evidence for the role played by the deformed Z = 36 proton shell of the light fragment in the fission of sub-lead nuclei. Our dataset will help constrain the fission models used to estimate the fission properties of nuclei with extreme neutron-to-proton ratios for which experimental data are unavailable.

In some scenarios, the dark matter relic abundance is set by the semi-annihilation of two dark matter particles into one dark matter particle and one Standard Model particle. These semi-annihilations might still be occurring today in the galactic center at a significant rate, generating a flux of boosted dark matter particles. We investigate the possible signals of this flux component in direct detection and neutrino experiments for sub-GeV dark matter masses. We show that for typical values of the semi-annihilation cross-section, the sensitivity of current experiments to the spin-independent dark matter-proton scattering cross-section can be several orders of magnitude larger than current constraints from cosmic-ray boosted dark matter. We also argue that the upcoming DARWIN and DUNE experiments may probe scattering cross-sections as low as <mml:math><mml:msup><mml:mrow><mml:mn>10</mml:mn></mml:mrow><mml:mrow><mml:mo>‑</mml:mo><mml:mn>37</mml:mn></mml:mrow></mml:msup><mml:mspace></mml:mspace><mml:msup><mml:mrow><mml:mi>cm</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:math> for masses between 30 MeV and 1 GeV.

Turbulent flows are known to produce enhanced effective magnetic and passive scalar diffusivities, which can fairly accurately be determined with numerical methods. It is now known that, if the flow is also helical, the effective magnetic diffusivity is reduced relative to the nonhelical value. Neither the usual second-order correlation approximation nor the various τ approaches have been able to capture this. Here we show that the helicity effect on the turbulent passive scalar diffusivity works in the opposite sense and leads to an enhancement. We have also demonstrated that the correlation time of the turbulent velocity field increases with the kinetic helicity. This is a key point in the theoretical interpretation of the obtained numerical results. Simulations in which helicity is being produced self-consistently by stratified rotating turbulence resulted in a turbulent passive scalar diffusivity that was found to be decreasing with increasing rotation rate.

We investigated the ability of the Euclid telescope to detect galaxy-scale gravitational lenses. To do so, we performed a systematic visual inspection of the 0.7 deg² Euclid Early Release Observations data towards the Perseus cluster using both the high-resolution I_E band and the lower-resolution Y_E , J_E, and H_E bands. Each extended source brighter than magnitude 23 in I_E was inspected by 41 expert human classifiers. This amounts to 12086 stamps of 10″ × 10″. We found 3 grade A and 13 grade B candidates. We assessed the validity of these 16 candidates by modelling them and checking that they are consistent with a single source lensed by a plausible mass distribution. Five of the candidates pass this check, five others are rejected by the modelling, and six are inconclusive. Extrapolating from the five successfully modelled candidates, we infer that the full 14 000 deg² of the Euclid Wide Survey should contain 100 000_{‑30 000}^{+ 70 000} galaxy-galaxy lenses that are both discoverable through visual inspection and have valid lens models. This is consistent with theoretical forecasts of 170 000 discoverable galaxy-galaxy lenses in Euclid. Our five modelled lenses have Einstein radii in the range 0'.'68 < θ_E < 1″.24, but their Einstein radius distribution is on the higher side when compared to theoretical forecasts. This suggests that our methodology is likely missing small-Einstein-radius systems. Whilst it is implausible to visually inspect the full Euclid dataset, our results corroborate the promise that Euclid will ultimately deliver a sample of around 10⁵ galaxy-scale lenses. ★This paper is published on behalf of the Euclid Consortium.

Mass transfer in binary systems is the key process in the formation of various classes of objects, including merging binary black holes (BBHs) and neutron stars. Orbital evolution during mass transfer depends on how much mass is accreted and how much angular momentum is lost $-$ two of the main uncertainties in binary evolution. Here, we demonstrate that there is a fundamental limit to how close binary systems can get via stable mass transfer (SMT), that is robust against uncertainties in orbital evolution. Based on detailed evolutionary models of interacting systems with a BH accretor and a massive star companion, we show that the post-interaction orbit is always wider than $\sim10R_{\odot}$, even when extreme shrinkage due to L2 outflows is assumed. Systems evolving towards tighter orbits become dynamically unstable and result in stellar mergers. This separation limit has direct implications for the properties of BBH mergers: long delay times ($\gtrsim1 \rm Gyr$), and no high BH spins from the tidal spin-up of helium stars. At high metallicity, the SMT channel may be severely quenched due to Wolf-Rayet winds. The reason for the separation limit lies in the stellar structure, not in binary physics. If the orbit gets too narrow during mass transfer, a dynamical instability is triggered by a rapid expansion of the remaining donor envelope due to its near-flat entropy profile. The closest separations can be achieved from core-He burning ($\sim8-15R_{\odot}$) and Main Sequence donors ($\sim15-30R_{\odot}$), while Hertzsprung Gap donors lead to wider orbits ($\gtrsim30-50R_{\odot}$) and non-merging BBHs. These outcomes and mass transfer stability are determined by the entropy structures, which are governed by internal composition profiles. The formation of compact binaries is thus sensitive to chemical mixing in stars and may relate to the blue-supergiant problem.

RNA and proteins are the foundation of life and a natural starting point to explore its origins. However, the prebiotic relationship between the two is asymmetric. While RNA evolved to assemble proteins from amino acids, a significant mirror-symmetric effect of amino acids to trigger the synthesis of RNA was missing. We describe ambient alkaline conditions where amino acids, without additional chemical activators, promote RNA copolymerisation more than 100-fold, starting from prebiotically plausible ribonucleoside-2′,3′-cyclic phosphates. The observed effect is explained by acid-base catalysis, with optimal efficiency at pH values near the amine pKaH. The fold-change in oligomerisation yield is nucleobase-selective, resulting in increased compositional diversity necessary for subsequent molecular evolution and favouring the formation of natural 3′−5′ linkages. The elevated pH offers recycling of oligonucleotide sequences back to 2′,3′-cyclic phosphates, providing conditions for high-fidelity replication by templated ligation. The findings reveal a clear functional role of amino acids in the evolution of RNA earlier than previously assumed.

The R³B (Reactions with Relativistic Radioactive Beams) experiment as a major instrument of the NUSTAR collaboration for the research facility FAIR in Darmstadt is designed for kinematically complete studies of reactions with high-energy radioactive beams. Part of the broad physics program of R³B is to constrain the asymmetry term in the nuclear equation-of-state and hence improve the description of highly asymmetric nuclear matter (e.g., in neutron stars). For a precise determination of the neutron-skin thickness – an observable which is directly correlated with the symmetry energy in theoretical calculations – by measuring absolute fragmentation cross sections, it is essential to quantify the uncertainty and challenge the reaction model under stable conditions. During the successful FAIR Phase-0 campaign of R³B, we precisely measured the energy dependence of total interaction cross sections in ¹²C+¹²C collisions, for a direct comparison with calculations based on the eikonal reaction theory.

We compute the partition function for the $N=1$ spinning particle, including pictures and the large Hilbert space, and show that it counts the dimension of the BRST cohomology in two- and four-dimensional target space. We also construct a quadratic action in the target space. Furthermore, we find a consistent interaction as a derived bracket based on the associative product of world line fields, leading to an interacting theory of multiforms in space-time. Finally, we comment on the equivalence of the multiform theory with a Dirac fermion. We also identify the chiral anomaly of the latter with a Hodge anomaly for the multiform theory, which manifests itself as a deformation of the gauge fixing.

So far, even the highest resolution galaxy formation simulations with gravitational softening have failed to reproduce realistic life cycles of star clusters. We present the first star-by-star galaxy models of star cluster formation to account for hydrodynamics, star formation, stellar evolution, and collisional gravitational interactions between stars and compact remnants using the updated SPHGAL + KETJU code, part of the GRIFFIN project. Gravitational dynamics in the vicinity of <inline-formula><tex-math id="TM0001" notation="LaTeX">$>3$</tex-math></inline-formula> M<inline-formula><tex-math id="TM0002" notation="LaTeX">$_\odot$</tex-math></inline-formula> stars and their remnants are solved with a regularized integrator (KETJU) without gravitational softening. Comparisons of idealized star cluster evolution with SPHGAL + KETJU and direct N-body show broad agreement and the failure of simulations that use gravitational softening. In the hydrodynamical simulations of idealized dwarf galaxies run with SPHGAL + KETJU, clusters up to <inline-formula><tex-math id="TM0003" notation="LaTeX">$\sim 900$</tex-math></inline-formula> M<inline-formula><tex-math id="TM0004" notation="LaTeX">$_\odot$</tex-math></inline-formula> form compact (effective radii 0.1-1 pc) and their sizes increase by up to a factor of 10 in agreement with previous N-body simulations and the observed sizes of exposed star clusters. The sizes increase rapidly once the clusters become exposed due to photoionizing radiation. On average 63 per cent of the gravitationally bound clusters disrupt during the first 100 Myr of evolution in the galactic tidal field. The addition of collisional dynamics reduces the fraction of supernovae in bound clusters by a factor of <inline-formula><tex-math id="TM0005" notation="LaTeX">$\sim 1.7$</tex-math></inline-formula>; however, the global star formation and outflow histories change by less than 30 per cent. We demonstrate that the accurate treatment of gravitational encounters with massive stars enables more realistic star cluster life cycles from the earliest stages of cluster formation until disruption in simulated low-mass galaxies.

Galaxy and halo scaling relations, connecting a broad range of parameters, are well established from observations. The origin of many of these relations and their scatter is still a matter of debate. It remains a sizable challenge for models to simultaneously and self-consistently reproduce as many scaling relations as possible. We introduce the Magneticum Pathfinder hydrodynamical cosmological simulation suite, to date the suite that self-consistently covers the largest range in box volumes and resolutions. It is the only cosmological simulation suite that is tuned on the hot gas content of galaxy clusters instead of the stellar mass function. By assessing the successes and shortcomings of tuning to the hot gas component of galaxy clusters, we aim to further our understanding of the physical processes shaping the Universe. We analyze the importance of the hot and cold gas components for galaxy and structure evolution. We analyze 28 scaling relations, covering large-scale global parameters as well as internal properties for halos ranging from massive galaxy clusters down to galaxies, and show their predicted evolution from z=4 to z=0 in comparison with observations. These include the halo-to-stellar-mass and Kennicutt--Schmidt relations, the cosmic star formation rate density as well as the Fundamental Plane. Magneticum Pathfinder matches a remarkable number of the observed scaling relations from z=4 to z=0, including challenging relations like the number density of quiescent galaxies at cosmic dawn, the mass--size evolution, the mass--metallicity relation, the Magorrian relation, and the temperature--mass relation. We compile our data to allow for straightforward future comparisons. Galaxy properties and scaling relations arise naturally and the large scatter in observables at high redshift is crucial to distinguish the various galaxy formation models reproducing the z=0 relations.

The decay of metastable 'false vacuum' states via bubble nucleation plays a crucial role in many cosmological scenarios. Cold-atom analog experiments will soon provide the first empirical probes of this process, with potentially far-reaching implications for early-Universe cosmology and high-energy physics. However, an inevitable difference between these analog systems and the early Universe is that the former have a boundary. We show, using a combination of Euclidean calculations and real-time lattice simulations, that these boundaries generically cause rapid bubble nucleation on the edge of the experiment, obscuring the bulk nucleation that is relevant for cosmology. We demonstrate that implementing a high-density 'trench' region at the boundary completely eliminates this problem, and recovers the desired cosmological behavior. Our findings are relevant for ongoing efforts to probe vacuum decay in the laboratory, providing a practical solution to a key experimental obstacle.

Future telescopes will characterize rocky exoplanets in reflected light, revealing their albedo, which depends on surface, cloud, and atmospheric properties. Identifying these features is crucial for assessing habitability. We present reference spectra and phase curves for an unresolved Earth-like exoplanet in reflected and polarized light, showing how phase- and wavelength-dependent reflectance reveals key planetary properties. Using the 3D radiative transfer code MYSTIC, we enhance surface and cloud modeling with validated, wavelength-dependent albedo maps of Earth's seasonal and spectral features, alongside a novel treatment of sub-grid cloud variability using ERA5 reanalysis data. Our models, incorporating high-resolution 3D cloud structures, show that sub-grid cloud variability reduces total reflectance and increases phase curve variability, especially at large phase angles where ocean glint dominates. We also find that neglecting wavelength-dependent albedo maps overestimates the vegetation red edge in spectra. Comparing an Ocean planet to an Earth-like planet with seasonal cloud variability, we show that polarization is more sensitive than intensity alone in distinguishing both cases. Moreover, polarization captures richer surface details, making it a crucial tool for resolving retrieval degeneracies. Our simulations serve as a reference for observing Earth as an exoplanet and provide benchmarks for optimizing observational strategies and retrieval frameworks for future telescopes targeting small, rocky exoplanets.

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.

We show experimentally that an effect of motion of ions, observed in a plasma-based accelerator, depends inversely on the plasma ion mass. The effect appears within a single wakefield event and manifests itself as a bunch tail, occurring only when sufficient motion of ions suppresses wakefields. Wakefields are driven resonantly by multiple bunches, and simulation results indicate that the ponderomotive force causes the motion of ions. In this case, the effect is also expected to depend on the amplitude of the wakefields, experimentally confirmed through variations in the drive bunch charge.

We calculate two-loop renormalization group equations (RGEs) in the Standard Model Effective Field Theory (SMEFT) with right-handed neutrinos, i.e., the so-called $ν$SMEFT, up to dimension five. Besides the two-loop RGEs of dimension-five (dim-5) operators, we also present those of the renormalizable couplings, including contributions from dim-5 operators. We check consistency relations among the first and second poles of $\varepsilon \equiv (4-d)/2$ with $d$ being the space-time dimension for all renormalization constants and find that those for lepton doublet and right-handed neutrino wave-function renormalization constants, as well as for renormalization constants of charged-lepton and neutrino Yukawa coupling matrices, do not hold. This leads to divergent RG functions for these fields and Yuwawa coupling matrices. We figure out that such infinite RG functions arise from the non-invariance of fields and Yukawa coupling matrices under field redefinitions, considering that flavor transformations are a kind of linear field redefinitions. Those infinite RG functions will disappear once one restores contributions from the derivative of renormalization constants with respect to the Wilson coefficients of redundant operators or, alternatively, considers the RGEs of flavor invariants, which are physical quantities and remain invariant under field redefinitions.

Context. Recent observations with the James Webb Space Telescope (JWST) have revealed unprecedented details of an intricate filamentary structure of unshocked ejecta within the young supernova remnant (SNR) Cassiopeia A (Cas A), offering new insights into the mechanisms governing supernova (SN) explosions and the subsequent evolution of ejecta. Aims. We aim to investigate the origin and evolution of the newly discovered web-like network of ejecta filaments in Cas A. Our specific objectives are: (i) to characterize the three-dimensional (3D) structure and kinematics of the filamentary network and (ii) to identify the physical mechanisms responsible for its formation. Methods. We performed high-resolution, 3D hydrodynamic (HD) and magneto-hydrodynamic (MHD) simulations to model the evolution of a neutrino-driven SN from the explosion to its remnant with the age of 1000 years. The initial conditions, set shortly after the shock breakout at the stellar surface, are based on a 3D neutrino-driven SN model that closely matches the basic properties of Cas A. Results. We found that the magnetic field has little impact on the evolution of unshocked ejecta, so we focused most of the analysis on the HD simulations. A web-like network of ejecta filaments, with structures compatible with those observed by JWST (down to scales ≈0.01 pc), naturally forms during the SN explosion. The filaments result from the combined effects of processes occurring soon after the core collapse, including the expansion of neutrino-heated bubbles formed within the first second after the explosion, hydrodynamic instabilities triggered during the blast propagation through the stellar interior, and the Ni-bubble effect following the shock breakout. The interaction of the reverse shock with the ejecta progressively disrupts the filaments through the growth of hydrodynamic instabilities. By around 700 years, the filamentary network becomes unobservable. Conclusions. According to our models, the filaments observed by JWST in Cas A most likely preserve a "memory" of the early explosion conditions, reflecting the processes active during and immediately after the SN event. Notably, a filamentary network closely resembling that observed in Cas A is naturally produced by a neutrino-driven SN explosion.

AT 2019aalc is a peculiar sequence of highly variable emission events observed towards the nucleus of the broad-line AGN SDSS J152416.66+045119.0. The system exhibited two distinct UV/optical flares (the first detected in 2019, the second one in 2023). Spectra obtained following the detection of the second flare revealed prominent Bowen fluorescence and high-ionization coronal emission lines, which were much weaker, if at all detectable, in a spectrum taken following the first flare. We present and analyze a large set of multi-wavelength, multi-epoch data for this source, with particular emphasis on optical spectroscopic monitoring conducted with the Las Cumbres Observatory network. During the relatively slow dimming that followed the second optical flare, the UV/optical light-curve shows a sequence of minor rebrightening events, while the Bowen fluorescence and the coronal lines vary (roughly) in tandem with these "bumps" in the broad-band light-curve. Most of the observed behavior of AT 2019aalc links it to the growing class of Bowen fluorescence flares (BFFs) while setting it apart from canonical tidal disruption events. However, AT 2019aalc has some outstanding peculiarities, including two short flares seen in its soft X-ray light-curve during the dimming phase of the second optical flare, and which do not seem to be linked to the emission line variations. We discuss the optical and X-ray properties of the source and possible scenarios of the origin of the flare, in particular radiation pressure instabilities in the (pre-existing) AGN accretion disk.

Through gravitational lensing, galaxy clusters can magnify supernovae (SNe) and create multiple images of the same SN. This enables measurements of cosmological parameters, which will be increasingly important in light of upcoming telescopic surveys. We study the prospects of detecting strongly lensed SNe in cluster fields with the Nancy Grace Roman Space Telescope (Roman)'s High Latitude Time Domain Survey (HLTDS) and the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST). We employed two approaches: one focusing on known multiply imaged galaxies behind clusters, along with the SN rates specific to those galaxies, and another based on the expected number of lensed SNe exploding in a given volume behind a galaxy cluster. We collected all the clusters in the literature that feature a well-constrained lens model and multiply imaged galaxies behind clusters with high-quality data for the lensed galaxies. This allowed us to determine the SN rate for each galaxy. We provide predictions for 46 clusters visible to the Vera C. Rubin Observatory, as well as for 9 observable by Roman's HLTDS, depending on whether the clusters fall within the survey's observing field. We predict that the number of multiply imaged SNe discovered by LSST in its first three years is $3.95 \pm 0.89$ from the first approach or $4.94 \pm 1.02$ from the second. For the HLTDS, the expected number of multiply imaged SNe ranges from $0.38 \pm 0.15$ to $5.2 \pm 2.2$, depending on the specific cluster observed, however, the fields to be targeted remain a matter of discussion. We conclude that LSST offers great prospects for detecting multiply imaged SNe. Our predictions are effectively lower limits, as we only considered the most massive and well-studied clusters. We provide a recommendation for HLTDS observing field selection, namely: either MACS J0553.4-3342 or Abell 1758a should be observed by the survey.

During the accretion phase of a core-collapse supernova (SN), dark-photon (DP) cooling can be largest in the gain layer below the stalled shock wave. In this way, it could counteract the usual shock rejuvenation by neutrino energy deposition and thus prevent the explosion. This peculiar energy-loss profile derives from the resonant nature of DP production. The largest cooling and thus strongest constraints obtain for DP masses of 0.1–0.4 MeV, a range corresponding to the photon plasma mass in the gain region. Electron-capture supernovae, once observationally unambiguously identified, could provide strong bounds even down to nearly 0.01 MeV. For a coupling strength so small that neutrino-driven explosions are expected to survive, the DP cooling of the core is too small to modify the neutrino signal, i.e., our new argument supersedes the traditional SN1987A cooling bound.

The Euclid Wide Survey (EWS) is predicted to find approximately 170 000 galaxy-galaxy strong lenses from its lifetime observation of 14 000 deg² of the sky. Detecting this many lenses by visual inspection with professional astronomers and citizen scientists alone is infeasible. As a result, machine learning algorithms, particularly convolutional neural networks (CNNs), have been used as an automated method of detecting strong lenses, and have proven fruitful in finding galaxy-galaxy strong lens candidates, such that the usage of CNNs in lens identification has increased. We identify the major challenge to be the automatic detection of galaxy-galaxy strong lenses while simultaneously maintaining a low false positive rate, thus producing a pure and complete sample of strong lens candidates from Euclid with a limited need for visual inspection. One aim of this research is to have a quantified starting point on the achieved purity and completeness with our current version of CNN-based detection pipelines for the VIS images of EWS. This work is vital in preparing our CNN-based detection pipelines to be able to produce a pure sample of the >100 000 strong gravitational lensing systems widely predicted for Euclid. We select all sources with VIS I_E < 23 mag from the Euclid Early Release Observation imaging of the Perseus field. We apply a range of CNN architectures to detect strong lenses in these cutouts. All our networks perform extremely well on simulated data sets and their respective validation sets. However, when applied to real Euclid imaging, the highest lens purity is just ∼11%. Among all our networks, the false positives are typically identifiable by human volunteers as, for example, spiral galaxies, multiple sources, and artifacts, implying that improvements are still possible, perhaps via a second, more interpretable lens selection filtering stage. There is currently no alternative to human classification of CNN-selected lens candidates. Given the expected ∼10⁵ lensing systems in Euclid, this implies 10⁶ objects for human classification, which while very large is not in principle intractable and not without precedent. ★This paper is published on behalf of the Euclid Consortium.

In the presence of a weak gravitational wave (GW) background, astrophysical binary systems act as high-quality resonators, with efficient transfer of energy and momentum between the orbit and a harmonic GW leading to potentially detectable orbital perturbations. In this work, we develop and apply a novel modeling and analysis framework that describes the imprints of GWs on binary systems in a fully time-resolved manner to study the sensitivity of lunar laser ranging, satellite laser ranging, and pulsar timing to both resonant and nonresonant GW backgrounds. We demonstrate that optimal data collection, modeling, and analysis lead to projected sensitivities which are orders of magnitude better than previously appreciated possible, opening up a new possibility for probing the physics-rich but notoriously challenging to access $μ\mathrm{Hz}$ frequency GWs. We also discuss improved prospects for the detection of the stochastic fluctuations of ultra-light dark matter, which may analogously perturb the binary orbits.

The M-theoretic emergence proposal claims that in an isotropic decompactification limit to M-theory the full effective action is generated via quantum effects by integrating out only the light towers of states of the theory. In the BPS particle sector, these include transversally wrapped $M2$- and $M5$-branes possibly carrying Kaluza-Klein momentum. This implies that a longitudinally wrapped $M5$-brane, i.e. a wrapped $D4$-brane, is not to be included in emergence computations. In this work we collect explicit evidence supporting this point by examining an $F^4$ gauge coupling in six dimensions, making use of the duality between heterotic string theory on $T^4$ and strongly coupled type IIA on K3. In this instance, the M-theoretic emergence proposal can be viewed as a tool for making predictions for the microscopic behavior of string theoretic amplitudes.

The authors present a new method to analytically prove global stability in ghost‑ridden dynamical systems. The proposal encompasses all prior results and consequentially extends them. In particular, it is shown that stability can follow from a conserved quantity that is unbounded from below, contrary to expectation. Novel examples illustrate all of the results. The findings take root on a careful examination of the literature, here comprehensively reviewed for the first time. This work lays the mathematical basis for ulterior extensions to field theory and quantization, and it constitutes a gateway for inter‑disciplinary research in dynamics and integrability.

Quasi-periodic eruptions (QPEs) are rapid, recurring X-ray bursts from supermassive black holes, believed to result from interactions between accretion disks and surrounding matter. The galaxy SDSS1335+0728, previously stable for two decades, exhibited an increase in optical brightness in December 2019, followed by persistent active galactic nucleus (AGN)-like variability for 5 yr, suggesting the activation of a ~10⁶-M_⊙ black hole. Since February 2024, X-ray emission has been detected, revealing extreme ~4.5-d QPEs with high fluxes and amplitudes, long timescales, large integrated energies and a ~25-d superperiod. Low-significance UV variations are reported, probably related to the long timescales and large radii from which the emission originates. This discovery broadens the possible formation channels for QPEs, suggesting that they are linked not solely to tidal disruption events but more generally to newly formed accretion flows, which we are witnessing in real time in a turn-on AGN candidate.

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

Finding the best parametrization for cosmological models in the absence of first-principle theories is an open question. We propose a data-driven parametrization of cosmological models given by the disentangled "latent" representation of a variational autoencoder (VAE) trained to compress cosmic microwave background (CMB) temperature power spectra. We consider a broad range of <inline-formula><mml:math><mml:mrow><mml:mi>Λ</mml:mi></mml:mrow></mml:math></inline-formula>-cold-dark-matter (<inline-formula><mml:math><mml:mi>λ</mml:mi><mml:mi>CDM</mml:mi></mml:math></inline-formula>) and beyond-<inline-formula><mml:math><mml:mi>Λ</mml:mi><mml:mi>CDM</mml:mi></mml:math></inline-formula> cosmologies with an additional early dark energy (EDE) component. We show that these spectra can be compressed into five (<inline-formula><mml:math><mml:mi>Λ</mml:mi><mml:mi>CDM</mml:mi></mml:math></inline-formula>) or eight (EDE) independent latent parameters, as expected when using temperature power spectra alone, and which reconstruct spectra at an accuracy well within the Planck errors. These latent parameters have a physical interpretation in terms of well-known features of the CMB temperature spectrum: these include the position, height and even-odd modulation of the acoustic peaks, as well as the gravitational lensing effect. The VAE also discovers one latent parameter which entirely isolates the EDE effects from those related to <inline-formula><mml:math><mml:mi>Λ</mml:mi><mml:mi>CDM</mml:mi></mml:math></inline-formula> parameters, thus revealing a previously unknown degree of freedom in the CMB temperature power spectrum. We further showcase how to place constraints on the latent parameters using Planck data as typically done for cosmological parameters, obtaining latent values consistent with previous <inline-formula><mml:math><mml:mi>Λ</mml:mi><mml:mi>CDM</mml:mi></mml:math></inline-formula> and EDE cosmological constraints. Our work demonstrates the potential of a data-driven reformulation of current beyond-<inline-formula><mml:math><mml:mi>Λ</mml:mi><mml:mi>CDM</mml:mi></mml:math></inline-formula> phenomenological models into the independent degrees of freedom to which the data observables are sensitive.

In some scenarios for the early universe, non-relativistic thermal dark matter chemically decouples from the thermal environment once the temperature drops well below the dark matter mass. The value at which the energy density freezes out depends on the underlying model. In a simple setting, we provide a comprehensive study of heavy fermionic dark matter interacting with the light degrees of freedom of a dark thermal sector whose temperature T decreases from an initial value close to the freeze-out temperature. Different temperatures imply different hierarchies of energy scales. By exploiting the methods of non-relativistic effective field theories at finite T, we systematically determine the thermal and in-vacuum interaction rates. In particular, we address the impact of the Debye mass on the bound-state formation cross section and the bound-state dissociation and transition widths, and ultimately on the dark matter relic abundance. We numerically compare the corrections to the present energy density originating from the resummation of Debye mass effects with the corrections coming from a next-to-leading order treatment of the bath-particle interactions. We observe that the fixed-order calculation of the inelastic heavy-light scattering at high temperatures provides a larger dark matter depletion, and hence an undersized yield for given benchmark points in the parameter space, with respect to the calculation where Debye mass effects are resummed.

That neutrinos carry a nonvanishing rest mass is evidence of physics beyond the Standard Model of elementary particles. Their absolute mass holds relevance in fields from particle physics to cosmology. We report on the search for the effective electron antineutrino mass with the KATRIN experiment. KATRIN performs precision spectroscopy of the tritium β-decay close to the kinematic endpoint. On the basis of the first five measurement campaigns, we derived a best-fit value of <inline-formula><mml:math><mml:mrow><mml:msubsup><mml:mi>m</mml:mi><mml:mi>ν</mml:mi><mml:mn>2</mml:mn></mml:msubsup><mml:mo>=</mml:mo><mml:mo>‑</mml:mo><mml:msubsup><mml:mrow><mml:mn>0.14</mml:mn></mml:mrow><mml:mrow><mml:mo>‑</mml:mo><mml:mn>0.15</mml:mn></mml:mrow><mml:mrow><mml:mo>+</mml:mo><mml:mn>0.13</mml:mn></mml:mrow></mml:msubsup></mml:mrow></mml:math></inline-formula> eV², resulting in an upper limit of m_ν < 0.45 eV at 90% confidence level. Stemming from 36 million electrons collected in 259 measurement days, a substantial reduction of the background level, and improved systematic uncertainties, this result tightens KATRIN's previous bound by a factor of almost two.

We study stellar core growth in simulations of merging massive (<inline-formula><tex-math id="TM0001" notation="LaTeX">$M_\star \gt 10^{11}\, \mathrm{M}_{\odot }$</tex-math></inline-formula>) elliptical galaxies by a supermassive black hole (SMBH) displaced by gravitational wave induced recoil velocity. With controlled, dense sampling of the SMBH recoil velocity, we find the core radius originally formed by SMBH binary scouring can grow by a factor of 2-3 when the recoil velocity exceeds <inline-formula><tex-math id="TM0002" notation="LaTeX">$\sim 50$</tex-math></inline-formula> per cent of the central escape velocity, and the mass deficit grows by up to a factor of <inline-formula><tex-math id="TM0003" notation="LaTeX">$\sim 4$</tex-math></inline-formula>. Using Bayesian inference we predict the distribution of stellar core sizes formed through this process to peak at <inline-formula><tex-math id="TM0004" notation="LaTeX">$\sim 1\, \mathrm{kpc}$</tex-math></inline-formula>. An orbital decomposition of stellar particles within the core reveals that radial orbits dominate over tube orbits when the recoil velocity exceeds the velocity dispersion of the core, whereas tube orbits dominate for the lowest recoil kicks. A change in orbital structure is reflected in the anisotropy parameter, with a central tangential bias present only for recoil velocities less than the local stellar velocity dispersion. Emulating current integral field unit observations of the stellar line-of-sight velocity distribution, we uncover a distinct signature in the Gauss-Hermite symmetric deviation coefficient <inline-formula><tex-math id="TM0005" notation="LaTeX">$h_4$</tex-math></inline-formula> that uniquely constrains the core size due to binary scouring. This signature is insensitive to the later evolution of the stellar mass distribution due to SMBH recoil. Our results provide a novel method to estimate the SMBH recoil magnitude from observations of local elliptical galaxies, and implies these galaxies primarily experienced recoil velocities less than the stellar velocity dispersion of the core.

The synthesis of life from non-living matter has captivated and divided scientists for centuries. This bold goal aims at unraveling the fundamental principles of life and leveraging its unique features, such as its resilience, sustainability, and ability to evolve. Synthetic life represents more than an academic milestone—it has the potential to revolutionize biotechnology, medicine, and materials science. Although the fields of synthetic biology, systems chemistry, and biophysics have made great strides toward synthetic life, progress has been hindered by social, philosophical, and technical challenges, such as vague goals, misaligned interdisciplinary efforts, and incompletely addressing public and ethical concerns. Our perspective offers a roadmap toward the synthesis of life based on discussions during a 2-week workshop with scientists from around the globe.

We present a formulation of coherent states as of consistent quantum description of classical configurations in the Becchi-Rouet-Stora-Tyutin (BRST)-invariant quantization of electrodynamics. The quantization with proper gauge-fixing is performed on the vacuum of the theory, whereas other backgrounds are obtained as BRST-invariant coherent states. One of the key insights is the possibility of constructing the coherent states of pure-gauge configurations. This provides a coherent state understanding of topologically nontrivial configurations in gauge theories and makes a number of features, such as the suppression of transitions between topologically distinct sectors, very transparent at full quantum level. As an example, we construct the Nielsen-Olesen string as a BRST-invariant coherent state. The Abelian pure-gauge configurations can also be viewed as useful analogs for a set of space-times related by coordinate reparametrizations in general relativity.

Historically, various methods have been employed to understand the origin of the elements, including observations of elemental abundances which have been compared to Galactic Chemical Evolution (GCE) models. It is also well known that 1D Local Thermodynamic Equilibrium (LTE) measurements fail to accurately capture elemental abundances. Non-LTE (NLTE) effects may play a significant role, and neglecting them leads to erroneous implications in galaxy modelling. In this paper, we calculate 3D NLTE abundances of seven key iron-peak and neutron-capture elements (Mn, Co, Ni, Sr, Y, Ba, Eu) based on carefully assembled 1D LTE literature measurements, and investigate their impact within the context of the OMEGA+ GCE model. Our findings reveal that 3D NLTE abundances are significantly higher for iron-peak elements at [Fe/H] <-3, with (for the first time ever) [Ni/Fe] and (confirming previous studies) [Co/Fe] on average reaching 0.6-0.8 dex, and [Mn/Fe] reaching -0.1 dex, which current 1D core-collapse supernova (CCSN) models cannot explain. We also observe a slightly higher production of neutron-capture elements at low metallicities, with 3D NLTE abundances of Eu being higher by +0.2 dex at [Fe/H] =-3. 3D effects are most significant for iron-peak elements in the very metal-poor regime, with average differences between 3D NLTE and 1D NLTE reaching up to 0.15 dex. Thus, ignoring 3D NLTE effects introduces significant biases, so including them should be considered whenever possible.

We introduce a novel orbit superposition method designed to reconstruct the stellar density structure, kinematics, and chemical abundance distribution of the entire Milky Way by leveraging 6D phase-space information from its resolved stellar populations, limited by the spatial coverage of APOGEE DR17.

Central regions of nearby disc galaxies display a large variety of structures in their stellar and gas disc that illustrates the outcome of a complex and dynamic evolution. The most visible central structure inside the bar region is the inner molecular disc (the so-called Central Molecular Zone (CMZ) in the Milky Way). Recent observational campaigns have shown that those inner molecular discs have a typical size ranging from a few hundred parsecs to one kiloparsec and tend to appear at the centre of barred discs. The physical phenomena involved in the building, consumption (e.g., star formation) and long-term evolution of those inner gas structures are still strongly debated. It is commonly accepted that the bar plays a role in the fuelling of gas from the large few kiloparsec scale to the inner molecular disc region. However, the exact physical processes involved in the loss of gas angular momentum (e.g., gravitational torques, shear, feedback) and its transport to the centre are not fully understood. Moreover, inner gas discs are intermediate-scale structures which connect the large kiloparsec scale and the subparsec scale physics (e.g., stellar-driven feedback, magnetic torques) involved in the fuelling of the central supermassive black hole (SMBH). Therefore, those inner discs can be considered as ‘gas reservoirs’ and may be efficient suppliers of material for the flickering of the central black hole, the so-called Active Galactic Nuclei (AGN). [...]

This study presents a novel method using spin quantum sensors to explore temporal variations of fundamental constants, significantly expanding the frequency range and providing constraints on scalar dark matter.

Understanding the ages of stars is crucial for unraveling the formation history and evolution of our Galaxy. Traditional methods for estimating stellar ages from spectroscopic data often struggle with providing appropriate uncertainty estimations and are severely constrained by the parameter space. In this work, we introduce a new approach using normalizing flows, a type of deep generative model, to estimate stellar ages for evolved stars with improved accuracy and robust uncertainty characterization. The model is trained on stellar masses for evolved stars derived from asteroseismology and predicts the relationship between the carbon and nitrogen abundances of a given star and its age. Unlike standard neural network techniques, normalizing flows enable the recovery of full likelihood distributions for individual stellar ages, offering a richer and more informative perspective on uncertainties. Our method yields age estimations for 378,720 evolved stars and achieves a typical absolute age uncertainty of approximately 2 Gyr. By intrinsically accounting for the coverage and density of the training data, our model ensures that the resulting uncertainties reflect both the inherent noise in the data and the completeness of the sampled parameter space. Applying this method to data from the SDSS-V Milky Way Mapper, we have produced the largest stellar age catalog for evolved stars to date.

We have carried out a systematic search for galaxy-scale lenses exploiting multi-band imaging data from the third public data release of the Hyper Suprime-Cam (HSC) survey with the focus on false-positive removal, after applying deep learning classifiers to all 110 million sources with i-Kron radius above 0.8". To improve the performance, we tested the combination of multiple networks from our previous lens search projects and found the best performance by averaging the scores from five of our networks. Although this ensemble network leads already to a false-positive rate (FPR) of 0.01% at a true-positive rate (TPR) of 75% on known real lenses, we have elaborated techniques to further clean the network candidate list before visual inspection. In detail, we tested the rejection using SExtractor and the modeling network from HOLISMOKES IX, which resulted together in a candidate rejection of 29% without lowering the TPR. We carried out a comprehensive multi-stage visual inspection involving eight individuals and identified 95 grade A (average grade G >2.5) and 503 grade B (2.5 >G >1.5) lens candidates, including 92 discoveries reported for the first time. This inspection also incorporated a novel environmental characterization using histograms of photometric redshifts. We publicly release the average grades, mass model predictions, and environment characterization of all visually inspected candidates, while including references for previously discovered systems, which makes this catalog one of the largest compilation of known lenses. The results demonstrate that (1) the combination of multiple networks enhances the selection performance and (2) both automated masking tools as well as modeling networks, which can be easily applied to hundreds of thousands of network candidates, help reduce the number of false positives that is the main limitation in lens search to date.

Context. Theoretical models of structure formation predict the presence of a hot gaseous atmosphere around galaxies. While this hot circumgalactic medium (CGM) has been observationally confirmed through UV absorption lines, the detection of its direct X-ray emission remains scarce. Recent results from the eROSITA collaboration have claimed the detection of the CGM out to the virial radius for a stacked sample of Milky Way-mass galaxies. Aims. 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 the gas temperature, hot gas fraction, and stellar mass. Methods. We selected a sample of central galaxies from the ultra-high-resolution cosmological volume (48 cMpc h^‑1) of the Magneticum Pathfinder set of hydrodynamical cosmological simulations. We classified them as star-forming (SF) or quiescent (QU) based on their specific star formation rate (SFR). For each galaxy, we generated X-ray mock data using the X-ray photon simulator PHOX, from which we obtained SB profiles out to the virial radius for different X-ray emitting components; namely, gas, active galactic nuclei (AGNs), and X-ray binaries (XRBs). We fit a β-profile to the gas component of each galaxy and observed trends between its slope and global quantities of the simulated galaxy. Results. We found marginal differences among the average total SB profile in SF and QU galaxies beyond r > 0.05 R_vir. The relative contribution from hot gas exceeds 70% and is non-zero (≲10%) for XRBs in both galaxy types. At small radii (r < 0.05 R_vir), XRBs dominate the SB profile over the hot gas for QU galaxies. We found positive correlations between the galaxies' global properties and the normalization of their SB profiles. The fitted β-profile slope is correlated with the total gas luminosity, which, in turn, shows strong connections to the current accretion rate of the central supermassive black hole (SMBH). We found the halo scaling relations to be consistent with the literature.

We introduce an extension of the evolution mapping framework to cosmological models that include massive neutrinos. The original evolution mapping framework exploits a degeneracy in the linear matter power spectrum when expressed in ${\rm Mpc}$ units, which compresses its dependence on cosmological parameters into those that affect its shape and a single extra parameter $\sigma_{12}$, defined as the RMS linear variance in spheres of radius $12 {\rm Mpc}$. We show that by promoting the scalar amplitude of fluctuations, $A_{\rm s}$, to a shape parameter, we can additionally describe the suppression due to massive neutrinos at any redshift to sub-0.01\% accuracy across a wide range of masses and for different numbers of mass eigenstates. This methodology has been integrated into the public COMET package, enhancing its ability to emulate predictions of state-of-the-art perturbative models for galaxy clustering, such as the effective field theory (EFT) model. Additionally, the updated software now accommodates a broader cosmological parameter space for the emulator, enables the simultaneous generation of multiple predictions to reduce computation time, and incorporates analytic marginalisation over nuisance parameters to expedite posterior estimation. Finally, we explore the impact of different infrared resummation techniques on galaxy power spectrum multipoles, demonstrating that any discrepancies can be mitigated by EFT counterterms without impacting the cosmological parameters.

Solid-state phonon and charge detectors probe the scattering of weakly interacting particles, such as dark matter and neutrinos, through their low recoil thresholds. Recent advancements have pushed sensitivity to eV-scale energy depositions, uncovering previously-unseen low-energy excess backgrounds. While some arise from known processes such as thermal radiation, luminescence, and stress, others remain unexplained. This review examines these backgrounds, their possible origins, and parallels to low-energy effects in solids. Their understanding is essential for interpreting particle interactions at and below the eV-scale.

Dark matter (DM) environments around black holes (BHs) can influence their mergers through dynamical friction, causing gravitational wave (GW) dephasing during the inspiral phase. While this effect is well studied for collisionless dark matter (CDM), it remains unexplored for self-interacting dark matter (SIDM) due to the typically low DM density in SIDM halo cores. In this work, we show that SIDM models with a massive force mediator can support dense enough DM spikes, significantly affecting BH mergers and producing a distinct GW dephasing. Using ${N}$-body simulations, we analyze GW dephasing in binary BH inspirals within CDM and SIDM spikes. By tracking the binary's motion in different SIDM environments, we show that the Laser Interferometer Space Antenna (LISA) can distinguish DM profiles shaped by varying DM interaction strengths, revealing detailed properties of SIDM.

We demonstrate that chiral symmetry breaking occurs in the confining regime of QCD-like theories with <mml:math altimg="si1.svg"><mml:msub><mml:mrow><mml:mi>N</mml:mi></mml:mrow><mml:mrow><mml:mi>c</mml:mi></mml:mrow></mml:msub></mml:math> colors and <mml:math altimg="si2.svg"><mml:msub><mml:mrow><mml:mi>N</mml:mi></mml:mrow><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub></mml:math> flavors. Our proof is based on a novel strategy, called 'downlifting', by which solutions of the 't Hooft anomaly matching and persistent mass conditions for a theory with <mml:math altimg="si3.svg"><mml:msub><mml:mrow><mml:mi>N</mml:mi></mml:mrow><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub><mml:mo linebreak="goodbreak" linebreakstyle="after">‑</mml:mo><mml:mn>1</mml:mn></mml:math> flavors are constructed from those of a theory with <mml:math altimg="si2.svg"><mml:msub><mml:mrow><mml:mi>N</mml:mi></mml:mrow><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub></mml:math> flavors, while <mml:math altimg="si1.svg"><mml:msub><mml:mrow><mml:mi>N</mml:mi></mml:mrow><mml:mrow><mml:mi>c</mml:mi></mml:mrow></mml:msub></mml:math> is fixed. By induction, chiral symmetry breaking is proven for any <mml:math altimg="si22.svg"><mml:msub><mml:mrow><mml:mi>N</mml:mi></mml:mrow><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub><mml:mo>≥</mml:mo><mml:msub><mml:mrow><mml:mi>p</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:math> in the confining regime, where <mml:math altimg="si5.svg"><mml:msub><mml:mrow><mml:mi>p</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:math> is the smallest prime factor of <mml:math altimg="si1.svg"><mml:msub><mml:mrow><mml:mi>N</mml:mi></mml:mrow><mml:mrow><mml:mi>c</mml:mi></mml:mrow></mml:msub></mml:math>. The proof can be extended to <mml:math altimg="si6.svg"><mml:msub><mml:mrow><mml:mi>N</mml:mi></mml:mrow><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub><mml:mo linebreak="goodbreak" linebreakstyle="after"><</mml:mo><mml:msub><mml:mrow><mml:mi>p</mml:mi></mml:mrow><mml:mrow><mml:mi>m</mml:mi><mml:mi>i</mml:mi><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:math> under the additional assumption on the absence of phase transitions when quark masses are sent to infinity. Our results do not rely on assumptions on the spectrum of massless bound states other than the fact that they are color-singlet hadrons.

We investigate the energy release in the interacting magnetospheres of binary neutron stars (BNSs) with global 3D force-free electrodynamics simulations. The system dynamics depend on the inclinations χ₁ and χ₂ of the stars' magnetic dipole moments relative to their orbital angular momentum. The simplest aligned configuration (χ₁ = χ₂ = 0^∘) has no magnetic field lines connecting the two stars. Remarkably, it still develops separatrix current sheets warping around each star and a dissipative region at the interface of the two magnetospheres. A Kelvin–Helmholtz (KH)–type instability drives significant dissipation at the magnetospheric interface, generating local Alfvénic turbulence and escaping fast magnetosonic waves. Binaries with inclined magnetospheres release energy in two ways: via KH instability at the interface and via magnetic reconnection flares in the twisted flux bundles connecting the companions. Outgoing compressive waves occur in a broad range of BNS parameters, possibly developing shocks and sourcing fast radio bursts. We discuss implications for X-ray and radio precursors of BNS mergers.

Robust modeling of non-linear scales is critical for accurate cosmological inference in Stage IV surveys. For weak lensing analyses in particular, a key challenge arises from the incomplete understanding of how non-gravitational processes, such as supernovae and active galactic nuclei — collectively known as baryonic feedback — affect the matter distribution. Several existing methods for modeling baryonic feedback treat it independently from the underlying cosmology, an assumption which has been found to be inaccurate by hydrodynamical simulations. In this work, we examine the impact of this coupling between baryonic feedback and cosmology on parameter inference at LSST Y1 precision. We build mock 3×2pt data vectors using the Magneticum suite of hydrodynamical simulations, which span a wide range of cosmologies while keeping subgrid parameters fixed. We perform simulated likelihood analyses for two baryon mitigation techniques: (i) the Principal Component Analysis (PCA) method which identifies eigenmodes for capturing the effect baryonic feedback on the data vector and (ii) HMCODE2020 [1] which analytically models the modification in the matter distribution using a halo model approach. Our results show that the PCA method is more robust than HMCODE2020 with biases in Ω_m-S ₈ up to 0.3σ and 0.6σ, respectively, for large deviations from the baseline cosmology. For HMCODE2020, the bias correlates with the input cosmology while for PCA we find no such correlation.

Cyanopolyynes are among the largest and most commonly observed interstellar complex organic molecules in star-forming regions. They are believed to form primarily in the gas phase, but their formation routes are not well understood. We present a comprehensive study of the gas-phase formation network of cyanobutadiyne, HC<inline-formula><tex-math id="TM0001" notation="LaTeX">$_5$</tex-math></inline-formula>N, based on new theoretical calculations, kinetics experiments, astronomical observations, and astrochemical modelling. We performed new quantum mechanics calculations for six neutral-neutral reactions in order to derive reliable rate coefficients and product branching fractions. We also present new CRESU data on the rate coefficients of three of these reactions (C<inline-formula><tex-math id="TM0002" notation="LaTeX">$_3$</tex-math></inline-formula>N + C<inline-formula><tex-math id="TM0003" notation="LaTeX">$_2$</tex-math></inline-formula>H<inline-formula><tex-math id="TM0004" notation="LaTeX">$_2$</tex-math></inline-formula>, C<inline-formula><tex-math id="TM0005" notation="LaTeX">$_2$</tex-math></inline-formula>H + HC<inline-formula><tex-math id="TM0006" notation="LaTeX">$_3$</tex-math></inline-formula>N, CN + C<inline-formula><tex-math id="TM0007" notation="LaTeX">$_4$</tex-math></inline-formula>H<inline-formula><tex-math id="TM0008" notation="LaTeX">$_2$</tex-math></inline-formula>) obtained at temperatures as low as 24 K. In practice, six out of nine reactions currently used in astrochemical models have been updated in our reviewed network. We also report the tentative detection of the <inline-formula><tex-math id="TM0009" notation="LaTeX">$^{13}$</tex-math></inline-formula>C isotopologues of HC<inline-formula><tex-math id="TM0010" notation="LaTeX">$_5$</tex-math></inline-formula>N in the L1544 prestellar core. We derived a lower limit of <inline-formula><tex-math id="TM0011" notation="LaTeX">$^{12}$</tex-math></inline-formula>C/<inline-formula><tex-math id="TM0012" notation="LaTeX">$^{13}$</tex-math></inline-formula>C > 75 for the HC<inline-formula><tex-math id="TM0013" notation="LaTeX">$_5$</tex-math></inline-formula>N isotopologues, which does not allow to bring new constraints to the HC<inline-formula><tex-math id="TM0014" notation="LaTeX">$_5$</tex-math></inline-formula>N chemistry. Finally, we verified the impact of the revised reactions by running the GRETOBAPE astrochemical model. We found good agreement between the HC<inline-formula><tex-math id="TM0015" notation="LaTeX">$_5$</tex-math></inline-formula>N predicted and observed abundances in cold (<inline-formula><tex-math id="TM0016" notation="LaTeX">$\sim$</tex-math></inline-formula>10 K) objects, demonstrating that HC<inline-formula><tex-math id="TM0017" notation="LaTeX">$_5$</tex-math></inline-formula>N is mainly formed by neutral-neutral reactions in these environments. In warm molecular shocks, instead, the predicted abundances are a factor of ten lower with respect to observed ones. In this environment possessing an higher gas ionization fraction, we speculate that the contribution of ion-neutral reactions could be significant.

A possible way of constructing polylogarithms on Riemann surfaces of higher genera facilitates integration kernels, which can be derived from generating functions incorporating the geometry of the surface. Functional relations among polylogarithms rely on identities for those integration kernels. In this article, we derive identities for Enriquez' meromorphic generating function and investigate the implications for the associated integration kernels. The resulting identities are shown to be exhaustive and therefore reproduce all identities for Enriquez' kernels conjectured in arXiv:2407.11476 recently.

The inner Solar System is depleted in refractory carbon in comparison to the interstellar medium and the depletion likely took place in the protoplanetary disk phase of the Solar System. We study the effect of photolysis of refractory carbon in the upper layers of the protosolar disk and its interplay with dust collisional growth and vertical mixing. We make use of a 1D Monte Carlo model to simulate dust coagulation and vertical mixing. To model the FUV flux of the disk, we use a simple analytical prescription and benchmark it with data from a radiative transfer simulation. We study the effects of fragmentation and bouncing on dust distribution and the propagation of carbon depletion. We find that when bouncing is included, the size distribution is truncated at smaller sizes than fragmentation-limited size distributions but there is a loss of small grains as well. The population of small grains is reduced due to fewer fragmentation events and this reduces the effectiveness of photolysis. We find that dust collisional growth and vertical mixing increase the effectiveness of carbon depletion by efficiently replenishing carbon to the upper regions of the disk with higher FUV flux. It takes around 100-300 kyr to reach the measured carbon abundances at 1 au, depending on the strength of the turbulence in the disk. These timescales are faster than reported by previous studies. Collisional redistribution and turbulent mixing are important aspects of dust evolution that should be included when modeling dust chemistry as they can influence the efficiency of chemical processes. Photolysis, along with another process such as sublimation, most likely played a key role in refractory carbon depletion that we see around us in the inner Solar System.

We present a systematic study of one-loop quantum corrections in scalar effective field theories from a geometric viewpoint, emphasizing the role of field-space curvature and its renormalisation. By treating the scalar fields as coordinates on a Riemannian manifold, we exploit field redefinition invariance to maintain manifest coordinate independence of physical observables. Focusing on the non-linear sigma model (NLSM) and \(\phi^4\) theory, we demonstrate how loop corrections induce momentum- and scale-dependent shifts in the curvature of the field-space manifold. These corrections can be elegantly captured through the recently proposed geometry-kinematics duality, which generalizes the colour-kinematics duality in gauge theories to curved field-space backgrounds. Our results highlight a universal structure emerging in the contractions of Riemann tensors that contribute to renormalisation of the field-space curvature. In particular, we find explicit expressions and a universal structure for the running curvature and Ricci scalar in simple models, illustrating how quantum effects reshape the underlying geometry. This geometric formulation unifies a broad class of scalar EFTs, providing insight into the interplay of curvature, scattering amplitudes, and renormalisation.

This paper presents a quantitative analysis of the stellar content in the Local Group dwarf irregular galaxy NGC 6822 by comparing stellar evolution models and observations in color-magnitude diagrams (CMDs) and color-color diagrams (CC-Ds). Our analysis is based on optical ground-based g,r,i photometry, and deep archive HST photometry of two fields in the galaxy disk. We compared young, intermediate-age, and old stellar populations with isochrones from the BaSTI-IAC library and found that NGC 6822 hosts a quite metal-rich ([Fe/H] = -0.7 to -0.4) young component with an age ranging from 20 to 100 Myr. The intermediate-age population experienced a modest chemical enrichment between 4 and 8 Gyr ago while stars older than 11 Gyr have a low metal abundance ([Fe/H] ~ -1.70). We also identified the AGB clump population with a luminosity peak at i ~ 23.35 mag. Our analysis of both the CMD and the optical-NIR-MIR CC-Ds of AGB oxygen- and carbon-rich stars, using the PARSEC+COLIBRI isochrones with and without circumstellar dust, reveal that this stellar component exhibits a spread in age from 1 to 2 Gyr and in metallicity between [Fe/H]=-1.30 and -1.70. The stellar models we used reproduce very well the two distinct color sequences defined by AGB O- and C-rich stars in the various optical-NIR-MIR CC-Ds, suggesting that they are reliable diagnostics to identify and characterise intermediate-age stellar populations. However, we also find that evolutionary prescriptions in the optical i-(r-i) CMDs predict, at fixed color, systematically lower luminosities than observed AGB stars.

Star clusters can interact and merge in galactic discs, halos, or centers. We present direct N-body simulations of binary mergers of star clusters with M_⋆ = 2.7 × 10⁴ M_⊙ each, using the N-body code BIFROSTwith subsystem regularisation and post-Newtonian dynamics. We include 500 M_⊙ massive black holes (MBHs) in the progenitors to investigate their impact on remnant evolution. The MBHs form hard binaries interacting with stars and stellar black holes (BHs). A few Myr after the cluster merger, this produces sizable populations of runaway stars (~800 with v_ej ≳ 50kms^-1) and stellar BHs (~30) escaping within 100 Myr. The remnants lose ~30% of their BH population and ~3% of their stars, with ~30 stars accelerated to high velocities ≳ 300kms^-1. Comparison simulations of isolated clusters with central hard MBH binaries and cluster mergers without MBHs show that the process is driven by MBH binaries, while those with a single 1000 M_⊙ MBH in isolated or merging clusters produce fewer runaway stars at lower velocities. Low-eccentricity merger orbits yield rotating remnants (v_rot ~ 3kms^-1) , but probing the presence of MBHs via kinematics alone remains challenging. We expect the binary MBHs to merge within a Hubble time, producing observable gravitational-wave (GW) events detectable by future GW detectors such as the Einstein Telescope and LISA. The results suggest that interactions with low-mass MBH binaries formed in merging star clusters are an important additional channel for producing runaway and high-velocity stars, free-floating stellar BHs and compact objects.

Ongoing and upcoming wide-field surveys at different wavelengths will measure the distribution of galaxy clusters with unprecedented precision, demanding accurate models for the two-point correlation function (2PCF) covariance. In this work, we assess a semi-analytical framework for the cluster 2PCF covariance that employs three nuisance parameters to account for non-Poissonian shot noise, residual uncertainties in the halo bias model, and subleading noise terms. We calibrate these parameters on a suite of fast approximate simulations generated by PINOCCHIO as well as full $N$-body simulations from OpenGADGET3. We demonstrate that PINOCCHIO can reproduce the 2PCF covariance measured in OpenGADGET3 at the few percent level, provided the mass functions are carefully rescaled. Resolution tests confirm that high particle counts are necessary to capture shot-noise corrections, especially at high redshifts. We perform the parameter calibration across multiple cosmological models, showing that one of the nuisance parameters, the non-Poissonian shot-noise correction $\alpha$, depends mildly on the amplitude of matter fluctuations $\sigma_8$. In contrast, the remaining two parameters, $\beta$ controlling the bias correction and $\gamma$ controlling the secondary shot-noise correction, exhibit more significant variation with redshift and halo mass. Overall, our results underscore the importance of calibrating covariance models on realistic mock catalogs that replicate the selection function of forthcoming surveys and highlight that approximate methods, when properly tuned, can effectively complement full $N$-body simulations for precision cluster cosmology.

Cosmic shear, galaxy clustering, and the abundance of massive halos each probe the large-scale structure of the Universe in complementary ways. We present cosmological constraints from the joint analysis of the three probes, building on the latest analyses of the lensing-informed abundance of clusters identified by the South Pole Telescope (SPT) and of the auto- and cross-correlation of galaxy position and weak lensing measurements (<inline-formula><mml:math display="inline"><mml:mn>3</mml:mn><mml:mo>×</mml:mo><mml:mn>2</mml:mn><mml:mi>pt</mml:mi></mml:math></inline-formula>) in the Dark Energy Survey (DES). We consider the cosmological correlation between the different tracers and we account for the systematic uncertainties that are shared between the large-scale lensing correlation functions and the small-scale lensing-based cluster mass calibration. Marginalized over the remaining <inline-formula><mml:math display="inline"><mml:mi mathvariant="normal">Λ</mml:mi></mml:math></inline-formula> cold dark matter (<inline-formula><mml:math display="inline"><mml:mi mathvariant="normal">Λ</mml:mi><mml:mi>CDM</mml:mi></mml:math></inline-formula>) parameters (including the sum of neutrino masses) and 52 astrophysical modeling parameters, we measure <inline-formula><mml:math display="inline"><mml:msub><mml:mi mathvariant="normal">Ω</mml:mi><mml:mi mathvariant="normal">m</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn>0.300</mml:mn><mml:mo>±</mml:mo><mml:mn>0.017</mml:mn></mml:math></inline-formula> and <inline-formula><mml:math display="inline"><mml:msub><mml:mi>σ</mml:mi><mml:mn>8</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mn>0.797</mml:mn><mml:mo>±</mml:mo><mml:mn>0.026</mml:mn></mml:math></inline-formula>. Compared to constraints from Planck primary cosmic microwave background (CMB) anisotropies, our constraints are only 15% wider with a probability to exceed of 0.22 (<inline-formula><mml:math display="inline"><mml:mn>1.2</mml:mn><mml:mi>σ</mml:mi></mml:math></inline-formula>) for the two-parameter difference. We further obtain <inline-formula><mml:math display="inline"><mml:msub><mml:mi>S</mml:mi><mml:mn>8</mml:mn></mml:msub><mml:mo>≡</mml:mo><mml:msub><mml:mi>σ</mml:mi><mml:mn>8</mml:mn></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi mathvariant="normal">Ω</mml:mi><mml:mi mathvariant="normal">m</mml:mi></mml:msub><mml:mo>/</mml:mo><mml:mn>0.3</mml:mn><mml:msup><mml:mo stretchy="false">)</mml:mo><mml:mn>0.5</mml:mn></mml:msup><mml:mo>=</mml:mo><mml:mn>0.796</mml:mn><mml:mo>±</mml:mo><mml:mn>0.013</mml:mn></mml:math></inline-formula> which is lower than the Planck measurement at the <inline-formula><mml:math display="inline"><mml:mn>1.6</mml:mn><mml:mi>σ</mml:mi></mml:math></inline-formula> level. The combined SPT cluster, DES <inline-formula><mml:math display="inline"><mml:mn>3</mml:mn><mml:mo>×</mml:mo><mml:mn>2</mml:mn><mml:mi>pt</mml:mi></mml:math></inline-formula>, and Planck datasets mildly prefer a nonzero positive neutrino mass, with a 95% upper limit <inline-formula><mml:math display="inline"><mml:mo>∑</mml:mo><mml:msub><mml:mi>m</mml:mi><mml:mi>ν</mml:mi></mml:msub><mml:mo><</mml:mo><mml:mn>0.25</mml:mn><mml:mtext> </mml:mtext><mml:mtext> </mml:mtext><mml:mi>eV</mml:mi></mml:math></inline-formula> on the sum of neutrino masses. Assuming a <inline-formula><mml:math display="inline"><mml:mi>w</mml:mi><mml:mi>CDM</mml:mi></mml:math></inline-formula> model, we constrain the dark energy equation of state parameter <inline-formula><mml:math display="inline"><mml:mi>w</mml:mi><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>1.1</mml:mn><mml:msubsup><mml:mn>5</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn>0.17</mml:mn></mml:mrow><mml:mrow><mml:mo>+</mml:mo><mml:mn>0.23</mml:mn></mml:mrow></mml:msubsup></mml:math></inline-formula> and when combining with Planck primary CMB anisotropies, we recover <inline-formula><mml:math display="inline"><mml:mi>w</mml:mi><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>1.2</mml:mn><mml:msubsup><mml:mn>0</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn>0.09</mml:mn></mml:mrow><mml:mrow><mml:mo>+</mml:mo><mml:mn>0.15</mml:mn></mml:mrow></mml:msubsup></mml:math></inline-formula>, a <inline-formula><mml:math display="inline"><mml:mn>1.7</mml:mn><mml:mi>σ</mml:mi></mml:math></inline-formula> difference with a cosmological constant. The precision of our results highlights the benefits of multiwavelength multiprobe cosmology and our analysis paves the way for upcoming joint analyses of next-generation datasets.

We investigate cosmological correlators for conformally coupled ϕ⁴ theory in four-dimensional de Sitter space. These in-in correlators differ from scattering amplitudes for massless particles in flat space due to the spacelike structure of future infinity in de Sitter. They also require a regularization which preserves de Sitter-invariance, which makes the flat space limit subtle to define at loop-level. Nevertheless we find that up to two loops, the in-in correlators are structurally simpler than the wave function and have the same transcendentality as flat space amplitudes. Moreover, we show that their loop integrands can be recast in terms of flat space integrands and can be derived from a novel recursion relation.

Neutron supermirrors (SMs) are a crucial part of many scattering and particle physics experiments. So far, Ni(Mo)/Ti SMs have been used in experiments that require to transport a polarized neutron beam due to their lower saturation magnetization compared to Ni/Ti SMs. However, next generation $\beta$ decay experiments require SMs that depolarize below $10^{-4}$ per reflection to reach their targeted precision. The depolarization of a polarized neutron beam due to reflection from Ni(Mo)/Ti SMs has not yet been measured to that precision. Recently developed Cu/Ti SMs with a very low saturation magnetization compared to Ni(Mo)/Ti may serve as an alternative. In this paper, we test the performance of both mirrors. At a first stage, we present four-states polarized neutron reflectivity (PNR) curves of Ni(Mo) and Cu monolayers measured at the neutron reflectometer SuperADAM and perform a full polarization analysis, showing a difference between the magnetic scattering length density (mSLD) of both materials, with Cu having a lower mSLD than Ni(Mo). These results are corroborated with the full polarization analysis of four-states PNR curves of $m=2$ Ni(Mo)/Ti and Cu/Ti SMs. In a second stage, we measured the depolarization ($D$) that a polarized neutron beam suffers after reflection from the same Ni(Mo)/Ti and Cu/Ti SMs by using the Opaque Test Bench setup. We find upper limits for the depolarization of $D_\text{Cu/Ti(4N5)}<7.6\times 10^{-5}$, $D_\text{Ni(Mo)/Ti}<8.5\times 10^{-5}$, and $D_\text{Cu/Ti(2N6)}<6.0\times 10^{-5}$ at the $1\sigma$ confidence level, where (4N5) corresponds to a Ti purity of $99.995\%$ and (2N6) to $99.6\%$. The uncertainties are statistical. These results show that all three SMs are suitable for being used in next generation $\beta$ decay experiments. We found no noticeable dependence of $D$ on the $q$ value or the magnetizing field, in which the samples were placed.

Motivated by JWST observations of dense, clumpy and clustered high redshift star formation, we simulate the hierarchical assembly of nine $M_{\mathrm{cl}}=10^6 M_\odot$ star clusters using the BIFROST N-body code. Our low metallicity models ($0.01Z_\odot$) with post-Newtonian equations of motion for black holes include evolving populations of single, binary and triple stars. Massive stars grow by stellar collisions and collapse into intermediate mass black holes (IMBHs) up to $M_\mathrm{\bullet}\sim6200 M_\odot$, stellar multiplicity boosting the IMBH masses by a factor of $2$--$3$. The IMBHs tidally disrupt (TDE) $\sim50$ stars in $10$ Myr with peak TDE rates of $\Gamma\sim10^{-5}$ yr$^{-1}$ per cluster. These IMBHs are natural seeds for supermassive black holes (SMBHs) and the hierarchical assembly frequently leads to $>2$ SMBH seeds per cluster and their rapid mergers ($t<10$ Myr). We propose that a gravitational wave (GW) driven merger of IMBHs with $1000 M_\odot \lesssim M_\bullet \lesssim 10000 M_\odot$ with comparable masses is a characteristic GW fingerprint of SMBH seed formation at redshifts $z>10$, and IMBH formation in similar environments at lower redshifts. Massive star clusters provide a unique environment for the early Universe GW studies for the next-generation GW observatories including the Einstein Telescope and the Laser Interferometer Space Antenna.

[Abridged] Cassiopeia A (Cas A) provides a unique opportunity to study supernova (SN) dynamics and interactions with the circumstellar medium (CSM). Recent JWST observations revealed the "Green Monster" (GM), a structure with a likely CSM origin. We investigate its pockmarked morphology, characterized by circular holes and rings, by examining the role of small-scale ejecta structures interacting with a dense circumstellar shell. We adopted a neutrino-driven SN model to trace the evolution of its explosion from core collapse to the age of the Cas A remnant using high-resolution 3D magnetohydrodynamic simulations. Besides other processes, the simulations include self-consistent calculations of radiative losses, accounting for deviations from electron-proton temperature equilibration and ionization equilibrium, as well as the ejecta composition derived from the SN. The GM's morphology is reproduced by dense ejecta clumps and fingers interacting with an asymmetric, forward-shocked circumstellar shell. The clumps and fingers form by hydrodynamic instabilities growing at the interface between SN ejecta and shocked CSM. Radiative cooling accounting for effects of non-equilibrium of ionization enhances the ejecta fragmentation, forming dense knots and thin filamentary structures that penetrate the shell, producing a network of holes and rings with properties similar to those observed. The origin of the holes and rings in the GM can be attributed to the interaction of ejecta with a shocked circumstellar shell. By constraining the timing of this interaction and analyzing the properties of these structures, we provide a distinction of this scenario from an alternative hypothesis, which attributes these features to fast-moving ejecta knots penetrating the shell ahead of the forward shock.

Aims. We use Gaia DR3 astrometry and photometry to analyze the spatial distribution of the young stellar populations and stellar clusters and to search for new OB star candidates in the Carina Nebula complex and the full extent (∼5°, corresponding to ∼200 pc) of the Car OB1 association. Methods. We first performed a new census of high-mass stars in Car OB1 and compiled a comprehensive catalog of 517 stars with known spectral types (128 O-type, WR, and supergiant stars, and 389 B-type stars) that have Gaia DR3 parallaxes consistent with membership in the association. We applied the clustering algorithm DBSCAN on the Gaia DR3 data of the region to find stellar clusters, determine their distances and kinematics, and estimate ages. We also used Gaia astrometry and the additional astrophysical_parameters table to perform a spatially unbiased search for further high-mass members of Car OB1 over the full area of the association. Results. Our DBSCAN analysis finds 15 stellar clusters and groups in Car OB1, four of which were not known before. Most clusters (80%) show signs of expansion or contraction, four of them with a ≥2σ significance. We find a global expansion of the Car OB1 association with a velocity of v_out = 5.25 ± 0.02 km s^‑1. A kinematic traceback of the high-mass stars shows that the spatial extent of the association was at a minimum 3–4 Myr ago. Using astrophysical parameters by Gaia DR3, we identified 15 new O-type and 589 new B-type star candidates in Car OB1. The majority (≳54%) of the high-mass stars constitute a non-clustered distributed stellar population. Based on our sample of high-mass stars, we estimate a total stellar population of at least ∼8 × 10⁴ stars in Car OB1. Conclusions. Our study is the first systematic astrometric analysis that covers the full spatial extent of the Car OB1 association, and it therefore substantially increases the knowledge of the distributed stellar population and spatial evolution of the entire association. Our results suggest suggests Car OB1 to be the most massive known star-forming complex in our Galaxy.

Radiation is crucial not only for observing astrophysical objects, but also for transporting energy and momentum. However, accurate on-the-fly radiation transport in astrophysical simulations is challenging and computationally expensive. Here we introduce AREPO-IDORT (Implicit Discrete Ordinates Radiation Transport), a scheme coupled to the explicit magnetohydrodynamic (MHD) solver in the 3D moving-mesh code AREPO. The discrete ordinates scheme means we directly solve for the specific intensities along discrete directions. We solve the time-dependent relativistic radiation transport equation via an implicit Jacobi-like iterative finite-volume solver, which overcomes the small radiation time-steps needed by explicit methods. Compared to commonly-used moment-based methods, e.g. flux-limited diffusion or M1 closure, this scheme has the advantage of correctly capturing the direction of radiation in both optically-thick and thin regions. It is based on the scheme by Jiang 2021 for the adaptive mesh refinement code ATHENA++, but we generalize the scheme to support (1) an unstructured moving-mesh, (2) local time-stepping, and (3) general equations of state. We show various test problems that commonly-used moment-based methods fail to reproduce accurately. To apply the scheme to a real astrophysics problem, we show the first global 3D radiation hydrodynamic simulation of the entire convective envelope of a red supergiant star. (abridged) For this problem, the radiation module only takes less than half of the total computational cost. Our current scheme assumes grey radiation, is first-order accurate in both time and space (abridged). We expect our scheme will enable more accurate multi-scale radiation MHD simulations involving supersonic bulk motions, ranging from planet formation in protoplanetary disks, stars and associated transients, to accretion flows near black holes.

Context. Brown dwarfs are the bridge between low-mass stars and giant planets. One way of shedding light on their dominant formation mechanism is to study them at the earliest stages of their evolution, when they are deeply embedded in their parental clouds. Several works have identified pre- and proto-brown dwarf candidates using different observational approaches. Aims. The aim of this work is to create a database of all the objects classified as very young substellar candidates in the literature in order to study them homogeneously. Methods. We gathered all the information about very young substellar candidates available in the literature until 2020. We retrieved their published photometry from the optical to the centimetre regime, and we wrote our own codes to derive their bolometric temperatures and luminosities, and their internal luminosities. We also populated the database with other parameters extracted from the literature, such as the envelope masses, their detection in some molecular species, and the presence of outflows. Results. The result of our search is the SUbstellar CANdidates at the Earliest Stages (SUCANES) database, containing 174 objects classified as potential very young substellar candidates in the literature. We present an analysis of the main properties of the retrieved objects. Since we updated the distances to several star forming regions, we were able to reject some candidates based on their internal luminosities. We also discuss the derived physical parameters and envelope masses for the best substellar candidates isolated in SUCANES. As an example of a scientific exploitation of this database, we present a feasibility study for the detection of radio jets with upcoming facilities: the next generation Very Large Array and the Square Kilometer Array interferometers. The SUCANES database is accessible through a graphical user interface, and it is open to any potential user.

We show how a method to construct canonical differential equations for multi-loop Feynman integrals recently introduced by some of the authors can be extended to cases where the associated geometry is of Calabi-Yau type and even beyond. This can be achieved by supplementing the method with information from the mixed Hodge structure of the underlying geometry. We apply these ideas to specific classes of integrals whose associated geometry is a one-parameter family of Calabi-Yau varieties, and we argue that the method can always be successfully applied to those cases. Moreover, we perform an in-depth study of the properties of the resulting canonical differential equations. In particular, we show that the resulting canonical basis is equivalent to the one obtained by an alternative method recently introduced in the literature. We apply our method to non-trivial and cutting-edge examples of Feynman integrals necessary for gravitational wave scattering, further showcasing its power and flexibility.

Cosmic-ray physics in the GeV-to-TeV energy range has entered a precision era thanks to recent data from space-based experiments. However, the poor knowledge of nuclear reactions, in particular for the production of antimatter and secondary nuclei, limits the information that can be extracted from these data, such as source properties, transport in the Galaxy and indirect searches for particle dark matter. The Cross-Section for Cosmic Rays at CERN workshop series has addressed the challenges encountered in the interpretation of high-precision cosmic-ray data, with the goal of strengthening emergent synergies and taking advantage of the complementarity and know-how in different communities, from theoretical and experimental astroparticle physics to high-energy and nuclear physics. In this paper, we present the outcomes of the third edition of the workshop that took place in 2024. We present the current state of cosmic-ray experiments and their perspectives, and provide a detailed road map to close the most urgent gaps in cross-section data, in order to efficiently progress on many open physics cases, which are motivated in the paper. Finally, with the aim of being as exhaustive as possible, this report touches several other fields -- such as cosmogenic studies, space radiation protection and hadrontherapy -- where overlapping and specific new cross-section measurements, as well as nuclear code improvement and benchmarking efforts, are also needed. We also briefly highlight further synergies between astroparticle and high-energy physics on the question of cross-sections.

Aims. We introduce the SISSI (Supernovae In a Stratified, Shearing Interstellar medium) simulation suite, which aims to enable a more comprehensive understanding of supernova remnants (SNRs) evolving in a complex interstellar medium (ISM) structured under the influence of galactic rotation, gravity and turbulence. Methods. We utilize zoom-in simulations of 30 SNRs expanding in the ISM of a simulated isolated disk galaxy. The ISM of the galaxy is resolved down to a maximum resolution of $\sim 12\,\text{pc}$, while we achieve a zoomed-in resolution of $\sim 0.18\, \text{pc}$ in the vicinity of the explosion sources. We compute the time-evolution of the SNRs' geometry and compare it to the observed geometry of the Local Bubble. Results. During the early stages of evolution, SNRs are well described by existing analytical models. On longer timescales, starting at about a percent of the orbital timescale, they depart from spherical symmetry and become increasingly prolate or oblate. The timescale for the departure from spherical symmetry is shorter than the expectation from a simple model for the deformation by galactic shear, suggesting that galactic shear alone cannot explain these differences. Yet, the alignment of the minor- and major axis of the SNRs is in line with expectations from said model, indicating that the deformation might have a shear-related origin. A comparison with the geometry of the Local Bubble reveals that it might be slightly younger than previously believed, but otherwise has a standard morphology for a SNR of its age and size. Conclusions. Studying the geometry of SNRs can reveal valuable insights about the complex interactions shaping their dynamical evolution. Future studies targeting the geometry of Galactic SNRs may use this insight to obtain a clearer picture of the processes shaping the Galactic ISM.

Neutrinos are the most abundant fundamental matter particles in the Universe and play a crucial role in particle physics and cosmology. Neutrino oscillation, discovered about 25 years ago, reveals that the three known species mix with each other. Anomalous results from reactor and radioactive-source experiments suggest a possible fourth neutrino state, the sterile neutrino, which does not interact via the weak force. The KATRIN experiment, primarily designed to measure the neutrino mass via tritium $\beta$-decay, also searches for sterile neutrinos suggested by these anomalies. A sterile-neutrino signal would appear as a distortion in the $\beta$-decay energy spectrum, characterized by a discontinuity in curvature (kink) related to the sterile-neutrino mass. This signature, which depends only on the shape of the spectrum rather than its absolute normalization, offers a robust, complementary approach to reactor experiments. KATRIN examined the energy spectrum of 36 million tritium $\beta$-decay electrons recorded in 259 measurement days within the last 40 electronvolt below the endpoint. The results exclude a substantial part of the parameter space suggested by the gallium anomaly and challenge the Neutrino-4 claim. Together with other neutrino-disappearance experiments, KATRIN probes sterile-to-active mass splittings from a fraction of an electron-volt squared to several hundred electron-volts squared, excluding light sterile neutrinos with mixing angles above a few percent.

The strongly lensed Supernova (SN) Encore at a redshift of $z = 1.949$, discovered behind the galaxy cluster MACS J0138$-$2155 at $z=0.336$, provides a rare opportunity for time-delay cosmography and studies of the SN host galaxy, where previously another SN, called SN Requiem, had appeared. To enable these studies, we combine new James Webb Space Telescope (JWST) imaging, archival Hubble Space Telescope (HST) imaging, and new Very Large Telescope (VLT) spectroscopic data to construct state-of-the-art lens mass models that are composed of cluster dark-matter (DM) halos and galaxies. We determine the photometric and structural parameters of the galaxies across six JWST and five HST filters. We use the color-magnitude and color-color relations of spectroscopically-confirmed cluster members to select additional cluster members, identifying a total of 84 galaxies belonging to the galaxy cluster. We construct seven different mass models using a variety of DM halo mass profiles, and explore both multi-plane and approximate single-plane lens models. As constraints, we use the observed positions of 23 multiple images from eight multiply lensed sources at four distinct spectroscopic redshifts. In addition, we use stellar velocity dispersion measurements to obtain priors on the galaxy mass distributions. We find that six of the seven models fit well to the observed image positions. Mass models with cored-isothermal DM profiles fit well to the observations, whereas the mass model with a Navarro-Frenk-White cluster DM profile has an image-position $\chi^2$ value that is four times higher. We build our ultimate model by combining four multi-lens-plane mass models and predict the image positions and magnifications of SN Encore and SN Requiem. Our work lays the foundation for building state-of-the-art mass models of the cluster for future cosmological analysis and SN host galaxy studies.

Precision spectroscopy of the electron spectrum of the tritium $\beta$-decay near the kinematic endpoint is a direct method to determine the effective electron antineutrino mass. The KArlsruhe TRItium Neutrino (KATRIN) experiment aims to determine this quantity with a sensitivity of better than 0.3$\,$eV (90$\,$% C.L.). An inhomogeneous electric potential in the tritium source of KATRIN can lead to distortions of the $\beta$-spectrum, which directly impact the neutrino-mass observable. This effect can be quantified through precision spectroscopy of the conversion-electrons of co-circulated metastable $^{83m}$Kr. Therefore, dedicated, several-weeks long measurement campaigns have been performed within the KATRIN data taking schedule. In this work, we infer the tritium source potential observables from these measurements, and present their implications for the neutrino-mass determination.

Understanding the coevolution of supermassive black holes and their host galaxies requires tracing their growth over time. Mass measurements of distant black holes have been limited to active nuclei and commonly rely on spatially unresolved observations, leading to large uncertainties. Accurate masses can be determined by resolving the kinematics of stars within the sphere of influence, which has heretofore been possible only in the local universe. Using JWST, we have measured the mass $M_{\bullet}=6.0^{+2.1}_{-1.7}\times10^9$ ${\rm M}_{\odot}$ of an inactive black hole in a gravitationally lensed quiescent galaxy at redshift $z=1.95$, along with detailed host properties. Comparisons to local galaxies suggest that the correlation between $M_{\bullet}$ and bulge mass has evolved substantially, whereas the correlation with stellar velocity dispersion may have been in place for 10 Gyr.

A puzzling population of extremely massive quiescent galaxies at redshifts beyond z = 3 has recently been revealed by JWST and the Atacama Large Millimeter/submillimeter Array, some of them with stellar ages that show their quenching times to be as high as z = 6, while their stellar masses are already above 5 × 10¹⁰ M _⊙. These extremely massive yet quenched galaxies challenge our understanding of galaxy formation at the earliest stages. Using the hydrodynamical cosmological simulation suite Magneticum Pathfinder, we show that such massive quenched galaxies at high redshifts can be successfully reproduced with similar number densities as observed. The stellar masses, sizes, formation redshifts, and star formation histories of the simulated quenched galaxies match those determined with JWST. Following these quenched galaxies at z = 3.4 forward in time, we find 20% to be accreted onto a more massive structure by z = 2, and from the remaining 80% about 30% rejuvenate up to z = 2, another 30% stay quenched, and the remaining 40% rejuvenate on a very low level of star formation. Stars formed through rejuvenation are mostly formed on the outer regions of the galaxies, not in the centers. Furthermore, we demonstrate that the massive quenched galaxies do not reside in the most massive nodes of the cosmic web, but rather live in side nodes of approximately Milky Way halo mass. Even at z = 0, only about 10% end up in small-mass galaxy clusters, while most of the quenched galaxies at z = 3.4 end up in group-mass halos, with about 20% actually not even reaching 10¹³ M _⊙ in halo mass.

Our picture of galaxy evolution currently assumes that galaxies spend their life on the star formation main sequence until they may eventually be quenched. However, recent observations show indications that the full picture might be more complicated. We reveal how the star formation rates of galaxies evolve, possible causes and imprints of different evolution scenarios on galactic features. We follow the evolution of central galaxies in the highest-resolution box of the Magneticum Pathfinder cosmological hydrodynamical simulations and classify their evolution scenarios with respect to the star formation main sequence. We find that a major fraction of the galaxies undergoes long-term cycles of quenching and rejuvenation on Gyr timescales. This expands the framework of galaxy evolution from a secular evolution to a sequence of multiple active and passive phases. Only 14% of field galaxies on the star formation main sequence at z~0 actually evolved along the scaling relation, while the bulk of star forming galaxies in the local universe have undergone cycles of quenching and rejuvenation. In this work we describe the statistics of these galaxy evolution modes and how this impacts their stellar masses, ages and metallicities today. Galaxies with rejuvenation cycles can be distinguished well from main-sequence-evolved galaxies in their features at z~0. We further explore possible explanations and find that the geometry of gas accretion at the halo outskirts shows a strong correlation with the star formation rate evolution, while the density parameter as a tracer of environment shows no significant correlation. A derivation of star formation rates from gas accretion with simple assumptions only works reasonably well in the high-redshift universe where accreted gas gets quickly converted into stars.

We present a systematic method for analytically computing time-dependent observables for a relativistic probe particle in Coulomb and Schwarzschild backgrounds. The method generates expressions valid both in the bound and unbound regimes, namely bound-unbound universal expressions. To demonstrate our method we compute the time-dependent radius and azimuthal angle for relativistic motion in a Coulomb background (relativistic Keplerian motion), as well as the electromagnetic field radiated by a relativistic Keplerian source. All of our calculations exhibit bound-unbound universality. Finally, we present an exact expression for the semi-classical wave function in Schwarzschild. The latter is crucial in applying our method to any time-dependent observable for probe-limit motion in Schwarzschild, to any desired order in velocity and the gravitational constant $G$.

By employing the potential non-relativistic quantum chromodynamics (pNRQCD) effective field theory within an open quantum system framework, we derive a Lindblad equation governing the evolution of the heavy-quarkonium reduced density matrix, accurate to next-to-leading order (NLO) in the ratio of the state's binding energy to the medium's temperature [1]. The derived NLO Lindblad equation provides a more reliable description of heavy-quarkonium evolution in the quark-gluon plasma at low temperatures compared to the leading-order truncation. For phenomenological applications, we numerically solve this equation using the quantum trajectories algorithm. By averaging over Monte Carlo-sampled quantum jumps, we obtain solutions without truncation in the angular momentum quantum number of the considered states. Our analysis highlights the importance of quantum jumps in the nonequilibrium evolution of bottomonium states within the quark-gluon plasma [2]. Additionally, we demonstrate that the quantum regeneration of singlet states from octet configurations is essential to explain experimental observations of bottomonium suppression. The heavy-quarkonium transport coefficients used in our study align with recent lattice QCD determinations.

Cosmic birefringence (CB) is the rotation of the photons' linear polarisation plane during propagation. Such an effect is a tracer of parity-violating extensions of standard electromagnetism and would probe the existence of a new cosmological field acting as dark matter or dark energy. It has become customary to employ cosmic microwave background (CMB) polarised data to probe such a phenomenon. Recent analyses on Planck and WMAP data provide a hint of detection of the isotropic CB angle with an amplitude of around $0.3^\circ$ at the level of $2.4$ to $3.6\sigma$. In this work, we explore the LiteBIRD capabilities in constraining such an effect, accounting for the impact of the more relevant systematic effects, namely foreground emission and instrumental polarisation angles. We build five semi-independent pipelines and test these against four different simulation sets with increasing complexity in terms of non-idealities. All the pipelines are shown to be robust and capable of returning the expected values of the CB angle within statistical fluctuations for all the cases considered. We find that the uncertainties in the CB estimates increase with more complex simulations. However, the trend is less pronounced for pipelines that account for the instrumental polarisation angles. For the most complex case analysed, we find that LiteBIRD will be able to detect a CB angle of $0.3^\circ$ with a statistical significance ranging from $5$ to $13 \, \sigma$, depending on the pipeline employed, where the latter uncertainty corresponds to a total error budget of the order of $0.02^\circ$.

Aims. We investigate the role of cosmic ray (CR) halos in shaping the properties of starburst-driven galactic outflows. Methods. We develop a microphysical model for galactic outflows driven by a continuous central feedback source, incorporating a simplified treatment of CRs. The model parameters are linked to the effective properties of a starburst. By analyzing its asymptotic behavior, we derive a criterion for launching starburst-driven galactic outflows and determine the corresponding outflow velocities. Results. We find that in the absence of CRs, galactic outflows can only be launched if the star-formation rate (SFR) surface density exceeds a critical threshold proportional to the dynamical equilibrium pressure. In contrast, CRs can always drive slow outflows. CRs dominate in systems with SFR surface densities below the critical threshold but become negligible in highly star-forming systems. However, in older systems with established CR halos, the CR contribution to outflows diminishes once the outflow reaches the galactic scale height, rendering CRs ineffective in sustaining outflows in such systems. Conclusions. Over cosmic time, galaxies accumulate relic CRs in their halos, providing additional non-thermal pressure support that suppresses low-velocity CR-driven outflows. We predict that such low-velocity outflows are expected only in young systems that have not yet built up significant CR halos. In contrast, fast outflows in starburst galaxies, where the SFR surface density exceeds the critical threshold, are primarily driven by momentum injection and remain largely unaffected by CR halos.

It has been shown that proton ingestion episodes can happen in the formation of hot-subdwarf stars, and that neutron-capture processes are possible in those cases. Moreover, some helium-rich hot subdwarfs display extraordinarily high abundances of heavy elements such as Zr, Yr and Pb on their surfaces. We explore under which conditions neutron-capture processes can occur in late helium core flashes, i.e. those occurring in the cores of stripped red-giant stars. We compute evolutionary models through the helium core flash and the subsequent hydrogen ingestion episode in stripped red giant stars. Stellar structure models are then used in post-processing to compute the detailed evolution of neutron-capture elements. We find that for metallicities of $10^{-3}$ and below, neutron densities can be as high as $10^{15}\,$cm$^{-3}$ and intermediate neutron capture processes occur in some of our models. The results depend very strongly on the H-envelope mass that survives after the stripping. Interestingly, we find that computed abundances in some of our models closely match the element abundances up to tin observed for EC 22536-5304, the only well-studied star for which the hot-flasher scenario assumed in our models is the most likely evolutionary path. Intermediate neutron capture processes can occur in the He-core flash experienced by the cores of some stripped red giants, and might be connected to the abundances of heavy elements observed in some helium-rich hot-subdwarf stars. The agreement between the observed abundances in EC 22536-5304 and those of our models offers support to our nucleosynthesis calculations. Moreover, if confirmed, the idea that heavy element abundances retain signatures of the different evolutionary channels opens the possibility that heavy element abundances in iHe-sdOB stars can be used to infer their evolutionary origin.

The latest generation of cosmic-ray direct detection experiments is providing a wealth of high-precision data, stimulating a very rich and active debate in the community on the related strong discovery and constraining potentials on many topics, namely dark matter nature, and the sources, acceleration, and transport of Galactic cosmic rays. However, interpretation of these data is strongly limited by the uncertainties on nuclear and hadronic cross-sections. This contribution is one of the outcomes of the \textit{Cross-Section for Cosmic Rays at CERN} workshop series, that built synergies between experimentalists and theoreticians from the astroparticle, particle physics, and nuclear physics communities. A few successful and illustrative examples of CERN experiments' efforts to provide missing measurements on cross-sections are presented. In the context of growing cross-section needs from ongoing, but also planned, cosmic-ray experiments, a road map for the future is highlighted, including overlapping or complementary cross-section needs from applied topics (e.g., space radiation protection and hadrontherapy).

Data from particle physics experiments are unique and are often the result of a very large investment of resources. Given the potential scientific impact of these data, which goes far beyond the immediate priorities of the experimental collaborations that obtain them, it is imperative that the collaborations and the wider particle physics community publish and preserve sufficient information to ensure that this impact can be realised, now and into the future. The information to be published and preserved includes the algorithms, statistical information, simulations and the recorded data. This publication and preservation requires significant resources, and should be a strategic priority with commensurate planning and resource allocation from the earliest stages of future facilities and experiments.

[abridged] AGN feedback is a crucial ingredient for understanding galaxy evolution. However, a complete quantitative time-dependent framework, including the dependence of such feedback on AGN, host galaxy, and host halo properties, is yet to be developed. Using the complete sample of 682 radio AGN from the LOFAR-eFEDS survey ($z0.4$), we derive the average jet power of massive galaxies and its variation as a function of stellar mass ($M_*$), halo mass ($M_h$) and radio morphology. We compare the incidence distributions of compact and complex radio AGN as a function of specific black hole kinetic power, $λ_{\rm Jet}$, and synthesise, for the first time, the radio luminosity function (RLF) by $M_*$ and radio morphology. Our RLF and derived total radio AGN kinetic luminosity density, $\log Ω_{\rm kin}/[\rm {W~Mpc^{-3}}]=32.15_{-0.34}^{+0.18}$, align with previous work. We find that kinetic feedback from radio AGN dominates over any plausible inventory of radiatively-driven feedback for galaxies with $\log M_*/M_\odot > 10.6$. More specifically, it is the compact radio AGN which dominate this global kinetic energy budget for all but the most massive galaxies ($10.6 < \log M_*/M_{\odot} < 11.5$). Subsequently, we compare the average injected jet energy against the galaxy and halo binding energy, and against the total thermal energy of the host gas within halos. We find that radio AGN cannot fully unbind their host galaxies nor host halos. However, they have enough energy to impact the global thermodynamical heating and cooling balance in small halos and significantly contribute to offsetting local cooling flows in even the most massive clusters cores. Overall, our findings provide important insights on jet powering, accretion processes and black hole-galaxy coevolution via AGN feedback, as well as a clear observational benchmark to calibrate AGN feedback simulations.

The measurement of the bound-state $β$ decay of $^{205}\mathrm{Tl}^{81+}$ at the Experimental Storage Ring at GSI, Darmstadt, has recently been reported with substantial impact on the use of $^{205}\mathrm{Pb}$ as an early Solar System chronometer and the low-energy measurement of the solar neutrino spectrum via the LOREX project. Due to the technical challenges in producing a high-purity $^{205}\mathrm{Tl}^{81+}$ secondary beam, a robust statistical method needed to be developed to estimate the variation in the contaminant $^{205}\mathrm{Pb}^{81+}$ produced during the fragmentation reaction. Here we show that Bayesian and Monte Carlo methods produced comparable estimates for the contaminant variation, each with unique advantages and challenges given the complex statistical problems for this experiment. We recommend the adoption of such methods in future experiments that exhibit unknown statistical fluctuations.

We study the deflection of light rays in a cold, nonmagnetized plasma using the worldline framework. Starting from Synge's Hamiltonian formalism, we construct a position-space action and use it perturbatively to calculate light bending angles. In the homogeneous case, the action reduces to that of a massive particle, allowing us to extract the bending angle of light in the presence of the medium using a well-known analogy. For the inhomogeneous case, we consider a power law model and construct Feynman rules in time to compute the purely plasma-induced corrections to the bending angle at next-to-leading-order.

Perturbative calculations for processes involving heavy flavours can be carried out using two approaches: the massive and the massless schemes. These schemes can also be combined to leverage their respective strengths. Additionally, both massive and massless frameworks can be supplemented by soft-gluon resummation. However, matching resummed calculations across the two schemes presents significant challenges, primarily due to the non-commutativity of the soft and small mass limits. The consistent resummation of mass and soft logarithms has been recently achieved at next-to-leading logarithmic (NLL) accuracy. In this paper, we consider heavy-quark fragmentation functions in electron-positron collisions and we extend this framework to achieve the so-called NLL<inline-formula><mml:math><mml:mmultiscripts><mml:mrow></mml:mrow><mml:mrow></mml:mrow><mml:mo>'</mml:mo></mml:mmultiscripts></mml:math></inline-formula> accuracy, which accounts for finite terms in the soft limit.

Cosmic Voids are a promising probe of cosmology for spectroscopic galaxy surveys due to their unique response to cosmological parameters. Their combination with other probes promises to break parameter degeneracies. Due to simplifying assumptions, analytical models for void statistics are only representative of a subset of the full void population. We present a set of neural-based emulators for void summary statistics of watershed voids, which retain more information about the full void population than simplified analytical models. We build emulators for the void size function and void density profiles traced by the halo number density using the Quijote suite of simulations for a broad range of the $\Lambda\mathrm{CDM}$ parameter space. The emulators replace the computation of these statistics from computationally expensive cosmological simulations. We demonstrate the cosmological constraining power of voids using our emulators, which offer orders-of-magnitude acceleration in parameter estimation, capture more cosmological information compared to analytic models, and produce more realistic posteriors compared to Fisher forecasts. We find that the parameters $\Omega_m$ and $\sigma_8$ in this Quijote setup can be recovered to $14.4\%$ and $8.4\%$ accuracy respectively using void density profiles; including the additional information in the void size function improves the accuracy on $\sigma_8$ to $6.8\%$. We demonstrate the robustness of our approach to two important variables in the underlying simulations, the resolution, and the inclusion of baryons. We find that our pipeline is robust to variations in resolution, and we show that the posteriors derived from the emulated void statistics are unaffected by the inclusion of baryons with the Magneticum hydrodynamic simulations. This opens up the possibility of a baryon-independent probe of the large-scale structure.

Collective oscillations in dense neutrino gases (flavor waves) are notable for their instabilities that cause fast flavor conversion. We develop a quantum theory of interacting neutrinos and flavor wave quanta, which are analogous to plasmons, but also carry flavor. The emission or absorption of such flavor plasmons $\psi$, or flavomons, changes the neutrino flavor. When an angular crossing occurs, the process $\nu_\mu\to\nu_e+\psi$ is more rapid than its inverse along the direction of the crossing, triggering stimulated $\psi$ emission and fast instability. Calculating the rate via Feynman diagrams matches the fast instability growth rate. Our novel $\nu$ and $\psi$ kinetic equations, corresponding to quasi-linear theory, describe instability evolution without resolving the small scales of the flavomon wavelength, potentially overcoming the main challenge of fast flavor evolution.

Ultra-hot Jupiters, an extreme class of planets not found in our solar system, provide a unique window into atmospheric processes. The extreme temperature contrasts between their day- and night-sides pose a fundamental climate puzzle: how is energy distributed? To address this, we must observe the 3D structure of these atmospheres, particularly their vertical circulation patterns, which can serve as a testbed for advanced Global Circulation Models (GCM) [e.g. 1]. Here, we show a dramatic shift in atmospheric circulation in an ultra-hot Jupiter: a unilateral flow from the hot star-facing side to the cooler space-facing side of the planet sits below an equatorial super-rotational jet stream. By resolving the vertical structure of atmospheric dynamics, we move beyond integrated global snapshots of the atmosphere, enabling more accurate identification of flow patterns and allowing for a more nuanced comparison to models. Global circulation models based on first principles struggle to replicate the observed circulation pattern [3], underscoring a critical gap between theoretical understanding of atmospheric flows and observational evidence. This work serves as a testbed to develop more comprehensive models applicable beyond our Solar System as we prepare for the next generation of giant telescopes.

We use the effective field theory approach to systematically study the dynamics of classical and quantum systems in an oscillating magnetic field. We find that the fast field oscillations give rise to an effective interaction which is able to confine charged particles as well as neutral particles with a spin magnetic moment. The effect is reminiscent of the renown dynamical stabilization of charges by the oscillating electric field and provides a foundation for a new class of magnetic traps. The properties characteristic to the dynamical magnetic confinement are reviewed.

Transit spectroscopy usually relies on the integration of one or several transits to achieve the S/N necessary to resolve spectral features. Consequently, high-S/N observations of exoplanet atmospheres are essential for disentangling the complex chemistry and dynamics beyond global trends. In this study, we combined two partial 4-UT transits of the ultrahot Jupiter WASP-121 b, observed with the ESPRESSO at the VLT in order to revisit its titanium chemistry. Through cross-correlation analysis, we achieved detections of H I, Li I, Na I, K I, Mg I, Ca I, Ti I, V I, Cr I, Mn I, Fe I, Fe II, Co I, Ni I, Ba II, Sr I, and Sr II. Additionally, narrow-band spectroscopy allowed us to resolve strong single lines, resulting in significant detections of H$\alpha$, H$\beta$, H$\gamma$, Li I, Na I, K I, Mg I, Ca II, Sr I, Sr II, and Mn I. Our most notable finding is the high-significance detection of Ti I ($\sim$ 5$\sigma$ per spectrum, and $\sim$ 19$\sigma$ stacked in the planetary rest frame). Comparison with atmospheric models reveals that Ti I is indeed depleted compared to V I. We also resolve the planetary velocity traces of both Ti I and V I, with Ti I exhibiting a significant blueshift toward the end of the transit. This suggests that Ti I primarily originates from low-latitude regions within the super-rotating jet observed in WASP-121 b. Our observations suggest limited mixing between the equatorial jet and the mid-latitudes, in contrast with model predictions from GCMs. We also report the non-detection of TiO, which we attribute to inaccuracies in the line list that could hinder its detection, even if present. Thus, the final determination of the presence of TiO must await space-based observations. We conclude that the 4-UT mode of ESPRESSO is an excellent testbed for achieving high S/N on relatively faint targets, paving the way for future observations with the ELT.

Observations of molecular lines are a key tool to determine the main physical properties of prestellar cores. However, not all the information is retained in the observational process or easily interpretable, especially when a larger number of physical properties and spectral features are involved. We present a methodology to link the information in the synthetic spectra with the actual information in the simulated models (i.e., their physical properties), in particular, to determine where the information resides in the spectra. We employ a 1D gravitational collapse model with advanced thermochemistry, from which we generate synthetic spectra. We then use neural network emulations and the SHapley Additive exPlanations (SHAP), a machine learning technique, to connect the models' properties to the specific spectral features. Thanks to interpretable machine learning, we find several correlations between synthetic lines and some of the key model parameters, such as the cosmic-ray ionization radial profile, the central density, or the abundance of various species, suggesting that most of the information is retained in the observational process. Our procedure can be generalized to similar scenarios to quantify the amount of information lost in the real observations. We also point out the limitations for future applicability.

We report the discovery and characterization of two sub-Saturns from the Transiting Exoplanet Survey Satellite (TESS) using high- resolution spectroscopic observations from the MaHPS spectrograph at the Wendelstein Observatory and the SOPHIE spectrograph at the Haute-Provence Observatory. Combining photometry from TESS, KeplerCam, LCOGT, and MuSCAT2, along with the radial velocity measurements from MaHPS and SOPHIE, we measured precise radii and masses for both planets. TOI-5108 b is a sub-Saturn, with a radius of 6.6 ± 0.1 R_⊕ and a mass of 32 ± 5 M_⊕. TOI-5786 b is similar to Saturn, with a radius of 8.54 ± 0.13 R_⊕ and a mass of 73 ± 9 M_⊕. The host star for TOI-5108 b is a moderately bright (Vmag 9.75) G-type star. TOI-5786 is a slightly dimmer (Vmag 10.2) F-type star. Both planets are close to their host stars, with periods of 6.75 days and 12.78 days, respectively. This puts TOI-5108 b just within the bounds of the Neptune desert, while TOI-5786 b is right above the upper edge. We estimated hydrogen-helium (H/He) envelope mass fractions of 38% for TOI-5108 b and 74% for TOI-5786 b. However, when using a model for the interior structure that includes tidal effects, the envelope fraction of TOI-5108 b could be much lower (~20%), depending on the obliquity. We estimated mass-loss rates between 1.0 x 10⁹ g/s and 9.8 x 10⁹ g/s for TOI-5108 b and between 3.6 x 10⁸ g/s and 3.5 x 10⁹ g/s for TOI-5786 b. Given their masses, both planets could be stable against photoevaporation. Furthermore, at these mass-loss rates, there is likely no detectable signal in the metastable helium triplet with the James Webb Space Telescope (JWST). We also detected a transit signal for a second planet candidate in the TESS data of TOI-5786, with a period of 6.998 days and a radius of 3.83 ± 0.16 R_⊕. Using our RV data and photodynamical modeling, we were able to provide a 3-σ upper limit of 26.5 M_⊕ for the mass of the potential inner companion to TOI-5786 b.

We transform the one-loop four-point type I open superstring gluon amplitude to correlation functions on the celestial sphere including both the (non-)orientable planar and non-planar sector. This requires a Mellin transform with respect to the energies of the scattered strings, as well as to integrate over the open-string worldsheet moduli space. After accomplishing the former we obtain celestial string integrands with remaining worldsheet integrals <mml:math altimg="si1.svg"><mml:mi mathvariant="normal">Ψ</mml:mi><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mi>β</mml:mi><mml:mo stretchy="true">)</mml:mo></mml:mrow></mml:math>, where β is related to the conformal scaling dimensions of the conformal primary operators under consideration. Employing an alternative approach of performing an <mml:math altimg="si2.svg"><mml:msup><mml:mrow><mml:mi>α</mml:mi></mml:mrow><mml:mrow><mml:mo>‧</mml:mo></mml:mrow></mml:msup></mml:math>-expansion of the open superstring amplitude first and Mellin transforming afterwards, we obtain a fully integrated expression, capturing the pole structure in the β-plane. The same analysis is performed at tree-level yielding similar results. We conclude by solving <mml:math altimg="si1.svg"><mml:mi mathvariant="normal">Ψ</mml:mi><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mi>β</mml:mi><mml:mo stretchy="true">)</mml:mo></mml:mrow></mml:math> for specific values of β, consistently reproducing the results of the <mml:math altimg="si2.svg"><mml:msup><mml:mrow><mml:mi>α</mml:mi></mml:mrow><mml:mrow><mml:mo>‧</mml:mo></mml:mrow></mml:msup></mml:math>-expansion ansatz. In all approaches we find that the dependence on <mml:math altimg="si2.svg"><mml:msup><mml:mrow><mml:mi>α</mml:mi></mml:mrow><mml:mrow><mml:mo>‧</mml:mo></mml:mrow></mml:msup></mml:math> reduces to that of a simple overall factor of <mml:math altimg="si3.svg"><mml:msup><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:msup><mml:mrow><mml:mi>α</mml:mi></mml:mrow><mml:mrow><mml:mo>‧</mml:mo></mml:mrow></mml:msup><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mrow><mml:mi>β</mml:mi><mml:mo linebreak="badbreak" linebreakstyle="after">‑</mml:mo><mml:mn>3</mml:mn></mml:mrow></mml:msup></mml:math> at loop and <mml:math altimg="si185.svg"><mml:msup><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:msup><mml:mrow><mml:mi>α</mml:mi></mml:mrow><mml:mrow><mml:mo>‧</mml:mo></mml:mrow></mml:msup><mml:mo stretchy="true">)</mml:mo></mml:mrow><mml:mrow><mml:mi>β</mml:mi></mml:mrow></mml:msup></mml:math> at tree level, consistent with previous literature.

We performed a detailed spectroscopic analysis of three extremely metal-poor RR Lyrae stars, exploring uncharted territories at these low metallicities for this class of stars. Using high-resolution spectra acquired with HARPS-N at TNG, UVES at VLT, and PEPSI at LBT, and employing Non-Local Thermodynamic Equilibrium (NLTE) spectral synthesis calculations, we provide abundance measurements for Fe, Al, Mg, Ca, Ti, Mn, and Sr. Our findings indicate that the stars have metallicities of [Fe/H] = ‑3.40 ± 0.05, ‑3.28 ± 0.02, and ‑2.77 ± 0.05 for HD 331986, DO Hya, and BPS CS 30317-056, respectively. Additionally, we derived their kinematic and dynamical properties to gain insights into their origins. Interestingly, the kinematics of one star (HD 331986) is consistent with the Galactic disc, while the others exhibit Galactic halo kinematics, albeit with distinct chemical signatures. We compared the [Al/Fe] and [Mg/Mn] ratios of the current targets with recent literature estimates to determine whether these stars were either accreted or formed in situ, finding that the adopted chemical diagnostics are ineffective at low metallicities ([Fe/H] ≲ ‑1.5). Finally, the established horizontal branch evolutionary models, indicating that these stars arrive at hotter temperatures on the Zero-Age Horizontal Branch (ZAHB) and then transition into RR Lyrae stars as they evolve, fully support the existence of such low-metallicity RR Lyrae stars. As a consequence, we can anticipate detecting more of them when larger samples of spectra become available from upcoming extensive observational campaigns. ⋆ Based on observations acquired at the Telescopio Nazionale Galileo under program A43DDT3, on a DDT program with PEPSI at LBT (2021-2022, PI Crestani) and on VLT ESO programs 69.C-0423(A) and 165.N-0276(A).

We present constraints on the <inline-formula><mml:math display="inline"><mml:mi>f</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>R</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> gravity model using a sample of 1005 galaxy clusters in the redshift range 0.25–1.78 that have been selected through the thermal Sunyaev-Zel'dovich effect from South Pole Telescope data and subjected to optical and near-infrared confirmation with the multicomponent matched filter algorithm. We employ weak gravitational lensing mass calibration from the Dark Energy Survey Year 3 data for 688 clusters at <inline-formula><mml:math display="inline"><mml:mi>z</mml:mi><mml:mo><</mml:mo><mml:mn>0.95</mml:mn></mml:math></inline-formula> and from the Hubble Space Telescope for 39 clusters with <inline-formula><mml:math display="inline"><mml:mn>0.6</mml:mn><mml:mo><</mml:mo><mml:mi>z</mml:mi><mml:mo><</mml:mo><mml:mn>1.7</mml:mn></mml:math></inline-formula>. Our cluster sample is a powerful probe of <inline-formula><mml:math display="inline"><mml:mi>f</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>R</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> gravity, because this model predicts a scale-dependent enhancement in the growth of structure, which impacts the halo mass function (HMF) at cluster mass scales. To account for these modified gravity effects on the HMF, our analysis employs a semianalytical approach calibrated with numerical simulations. Combining calibrated cluster counts with primary cosmic microwave background temperature and polarization anisotropy measurements from the Planck 2018 release, we derive robust constraints on the <inline-formula><mml:math display="inline"><mml:mi>f</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>R</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> parameter <inline-formula><mml:math display="inline"><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mi>R</mml:mi><mml:mn>0</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula>. Our results, <inline-formula><mml:math display="inline"><mml:msub><mml:mi>log</mml:mi><mml:mn>10</mml:mn></mml:msub><mml:mo stretchy="false">|</mml:mo><mml:msub><mml:mi>f</mml:mi><mml:mrow><mml:mi>R</mml:mi><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo stretchy="false">|</mml:mo><mml:mo><</mml:mo><mml:mo>-</mml:mo><mml:mn>5.32</mml:mn></mml:math></inline-formula> at the 95% credible level, are the tightest current constraints on <inline-formula><mml:math display="inline"><mml:mi>f</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>R</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula> gravity from cosmological scales. This upper limit rules out <inline-formula><mml:math display="inline"><mml:mi>f</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>R</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:math></inline-formula>-like deviations from general relativity that result in more than a <inline-formula><mml:math display="inline"><mml:mo>∼</mml:mo><mml:mn>20</mml:mn><mml:mo>%</mml:mo></mml:math></inline-formula> enhancement of the cluster population on mass scales <inline-formula><mml:math display="inline"><mml:mrow><mml:msub><mml:mrow><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mn>200</mml:mn><mml:mi mathvariant="normal">c</mml:mi></mml:mrow></mml:msub><mml:mo>></mml:mo><mml:mn>3</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mrow><mml:mn>10</mml:mn></mml:mrow><mml:mrow><mml:mn>14</mml:mn></mml:mrow></mml:msup><mml:msub><mml:mrow><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false">⊙</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>.

In the present study, we employ three distinct, physically motivated speed of sound bounds to construct hybrid models, where the high-density phase is described by the maximally stiff equation of state. In particular, we consider the bounds related to special relativity, relativistic kinetic theory, and conformality. The low-density hadronic phase is described by a state-of-the-art microscopic relativistic Brueckner-Hartree-Fock theory. This work aims to access the effect of the different speed of sound constraints on the relevant parameter space of the key parameters of first-order phase transitions by utilizing recent astronomical data. This involves a systematic analysis that also includes two distinct schemes for the construction of hybrid models (abrupt and smooth). Finally, a relevant discussion is conducted on the possible occurrence of a thermodynamic inconsistency that is related to the stability of the high-density phase over hadronic matter at large densities.

Context. The intense X-ray and UV emission of some active M stars has raised questions about the habitability of planets around M-type stars. Aims. We aim to determine the unbiased distribution of X-ray luminosities in complete, volume-limited samples of nearby M dwarfs, and compare them to those of K and G dwarfs. Methods. We constructed volume-complete samples of 205 M stars with a spectral type ≤ M6 within 10 pc of the Sun, 129 K stars within 16 pc, and 107 G stars within 20 pc. We used X-ray data from Chandra, XMM-Newton, eROSITA, and ROSAT to obtain the X-ray luminosities of the stars. Results. Our samples reach an X-ray detection completeness of 85%, 86%, and 80% for M, K, and G stars, respectively. The fractional X-ray luminosities relative to the bolometric luminosities, log(L_X/L_bol), of the M stars show a bimodal distribution, with one peak at around ‑5, mostly contributed by early M stars (M0–M4), and another peak around ‑3.5, contributed mainly by M4–M6 stars. The comparison of the different spectral classes shows that 63% of all M stars in our sample (80% of the M stars with a spectral type < M4) have L_X/L_bol values that are within the central 80% quantile of the distribution function for G stars. In addition, 55% of all M stars in our sample (and 72% of the M stars with a spectral type < M4) have L_X/L_bol less than 10 times the solar value. Conclusions. The X-ray activity levels of the majority (≳60%) of nearby M dwarfs no later than M6 are actually not higher than the typical (80% quantile) levels for G-type stars. The X-ray irradiation of habitable-zone planets around these stars should therefore not present a specific problem for their habitability.

The formation details of globular clusters (GCs) are still poorly understood due to their old ages and the lack of detailed observations of their formation. A large variety of models for the formation and evolution of GCs have been created to improve our understanding of their origins, based on GC properties observed at <inline-formula><tex-math id="TM0001" notation="LaTeX">$z=0$</tex-math></inline-formula>. We present the first side-by-side comparison of six current GC formation models with respect to their predictions for the GC ages and formation redshifts in Milky Way (MW)-like galaxies. We find that all the models are capable of forming most of the surviving GCs at more than <inline-formula><tex-math id="TM0002" notation="LaTeX">$10 \,\mathrm{G}{\rm {yr}}$</tex-math></inline-formula> ago, in general agreement with the observation that most GCs are old. However, the measured MW GC ages are still systematically older than those predicted in the galaxies of four of the models. Investigating the variation of modelled GC age distributions for general MW-mass galaxies, we find that some of the models predict that a significant fraction of MW-mass galaxies would entirely lack a GC population older than <inline-formula><tex-math id="TM0003" notation="LaTeX">$10 \,\mathrm{G}{\rm {yr}}$</tex-math></inline-formula>, whereas others predict that all MW-mass galaxies have a significant fraction of old GCs. This will have to be further tested in upcoming surveys, as systems without old GCs in that mass range are currently not known. Finally, we show that the models predict different formation redshifts for the oldest surviving GCs, highlighting that models currently disagree about whether the recently observed young star clusters at high redshifts could be the progenitors of today's GCs.

Strongly lensed supernovae (SNe) are a rare class of transient that can offer tight cosmological constraints that are complementary to methods from other astronomical events. We present a follow-up study of one recently discovered strongly lensed SN, the quadruply imaged type Ia SN 2022qmx (aka "SN Zwicky"), at z = 0.3544. We measure updated, template-subtracted photometry for SN Zwicky and derive improved time delays and magnifications. This is possible because SNe are transient, fading away after reaching their peak brightness. Specifically, we measure point-spread-function photometry for all four images of SN Zwicky in three Hubble Space Telescope WFC3/UVIS passbands (F475W, F625W, and F814W) and one WFC3/IR passband (F160W), with template images taken ∼11 months after the epoch in which the SN images appear. We find consistency to within 2σ between lens-model-predicted time delays (≲1 day) and measured time delays with HST colors (≲2 days), including the uncertainty from chromatic microlensing that may arise from stars in the lensing galaxy. The standardizable nature of SNe Ia allows us to estimate absolute magnifications for the four images, with images A and C being elevated in magnification compared to lens model predictions by about 6σ and 3σ, respectively, confirming previous work. We show that millilensing or differential dust extinction is unable to explain these discrepancies, and we find evidence for the existence of microlensing in images A, C, and potentially D that may contribute to the anomalous magnification.

We identify a chain of galaxies along an almost straight line in the nearby Universe with a projected length of ~5 Mpc. The galaxies are distributed within projected distances of only 7–105 kpc from the axis of the identified filament. They have redshifts in a very small range of z = 0.0361‑0.0370 so that their radial velocities are consistent with galaxy proper motions. The filament galaxies are mainly star forming and have stellar masses in a range of 10^9.1‑10^10.7 M_⊙. We search for systems with similar geometrical properties in the full-sky mock galaxy catalog of the MillenniumTNG simulations and find that, although such straight filaments are unusual and rare, they are predicted by ΛCDM simulations (4% incidence). We study the cold H I gas in a 1.3 Mpc section of the filament through H I 21 cm emission line observations and detect 11 H I sources, many more than expected from the H I mass function in a similar volume. They have H I masses 10^8.5‑10^9.5 M_⊙ and are mostly within ~120 kpc projected distance from the filament axis. None of these H I sources has a confirmed optical counterpart. Their darkness together with their large H I 21 cm line widths indicates that they contain gas that might not yet be virialized. These clouds must be marking the peaks of the dark matter and H I distributions over large scales within the filament. The presence of such gas clouds around the filament spines is predicted by simulations, but this is the first time that the existence of such clouds in a filament is observationally confirmed.

Observed low-mass galaxies with nuclear star clusters (NSCs) can host accreting massive black holes (MBH). We present simulations of dwarf galaxies (<inline-formula><tex-math id="TM0001" notation="LaTeX">$M_{\mathrm{baryon}} \sim 0.6\!-\!2.4 \times 10^8 \rm \, M_\odot$</tex-math></inline-formula>) at solar mass resolution (<inline-formula><tex-math id="TM0002" notation="LaTeX">$0.5\rm \, M_\odot \lt \mathit{ m}_{\mathrm{gas}} \lt 4 \rm \, M_\odot$</tex-math></inline-formula>) with a multiphase interstellar medium (ISM) and investigate the impact of NSCs on MBH growth and nuclear star formation (SF). The GRIFFIN simulation model includes non-equilibrium low temperature cooling, chemistry and the effect of H II regions and supernovae (SNe) from massive stars. Individual stars are sampled down to 0.08 <inline-formula><tex-math id="TM0003" notation="LaTeX">$\rm M_\odot$</tex-math></inline-formula> and their non-softened gravitational interactions with MBHs are computed with the regularized KETJU integrator. MBHs with masses in the range of <inline-formula><tex-math id="TM0004" notation="LaTeX">$10^2 \!-\! 10^5 \, \rm M_\odot$</tex-math></inline-formula> are represented by accreting sink particles without feedback. We find that the presence of NSCs boost nuclear SF (i.e. NSC growth) and MBH accretion by funneling gas to the central few parsecs. Low-mass MBHs grow more rapidly on <inline-formula><tex-math id="TM0005" notation="LaTeX">$\sim 600$</tex-math></inline-formula> Myr time-scales, exceeding their Eddington rates at peak accretion. MBH accretion and nuclear SF is episodic (i.e. leads to multiple stellar generations), coeval and regulated by SN explosions. On 40-60 Myr time-scales the first SN of each episode terminates MBH accretion and nuclear SF. Without NSCs, low-mass MBHs do not grow and MBH accretion and reduced nuclear SF become irregular and uncorrelated. This study gives the first insights into the possible co-evolution of MBHs and NSCs in low-mass galaxies and highlights the importance of considering dense NSCs in galactic studies of MBH growth.

Using two-dimensional general relativistic resistive magnetohydrodynamic simulations, we investigate the properties of the sheath separating the black hole jet from the surrounding medium. We find that the electromagnetic power flowing through the jet sheath is comparable to the overall accretion power of the black hole. The sheath is an important site of energy dissipation as revealed by the copious appearance of reconnection layers and plasmoid chains. About 20% of the sheath power is dissipated between 2 and 10 gravitational radii. The plasma in the dissipative sheath moves along a nearly paraboloidal surface with transrelativistic bulk motions dominated by the radial component, whose dimensionless 4-velocity is ∼1.2 ± 0.5. In the frame moving with the mean (radially dependent) velocity, the distribution of stochastic bulk motions resembles a Maxwellian with an "effective bulk temperature" of ∼100 keV. Scaling the global simulation to Cygnus X-1 parameters gives a rough estimate of the Thomson optical depth across the jet sheath, ∼0.01–0.1, and it may increase in future magnetohydrodynamic simulations with self-consistent radiative losses. These properties suggest that the dissipative jet sheath may be a viable "coronal" region, capable of upscattering seed soft photons into a hard, nonthermal tail, as seen during the hard states of X-ray binaries and active galactic nuclei.

Cosmic-ray acceleration processes in astrophysical plasmas are often investigated with fully kinetic or hybrid kinetic numerical simulations, which enable us to describe a detailed microphysics of particle energization mechanisms. Tracing of individual particles in such simulations is especially useful in this regard. However, visually inspecting particle trajectories introduces a significant amount of bias and uncertainty, making it challenging to pinpoint specific acceleration mechanisms. Here, we present a novel approach utilizing neural networks to assist in the analysis of individual particle data. We demonstrate the effectiveness of this approach using the dataset from our recent particle-in-cell (PIC) simulations of non-relativistic perpendicular shocks, which consists of 252 000 electrons, each characterized by their position, momentum, and electromagnetic field at particle's position, recorded in a time series of 1200 time steps. These electrons cross a region affected by the electrostatic Buneman instability, and a small percentage of them attain high energies. We perform classification, regression, and anomaly detection algorithms on the dataset by using a convolutional neural network, a multi-layer perceptron, and an autoencoder. Despite the noisy and imbalanced dataset, all methods demonstrate the capability to differentiate between thermal and accelerated electrons with remarkable accuracy. The proposed methodology may considerably simplify particle classification in large-scale PIC and hybrid simulations.

The braking torque that dictates the timing properties of magnetars is closely tied to the large-scale dipolar magnetic field on their surface. The formation of this field has been a topic of ongoing debate. One proposed mechanism, based on macroscopic principles, involves an inverse cascade within the neutron star's crust. However, this phenomenon has not been observed in realistic simulations. In this study, we provide compelling evidence supporting the feasibility of the inverse cascading process in the presence of an initial helical magnetic field within realistic neutron star crusts and discuss its contribution to the amplification of the large-scale magnetic field. Our findings, derived from a systematic investigation that considers various coordinate systems, peak wavenumber positions, crustal thicknesses, magnetic boundary conditions, and magnetic Lundquist numbers, reveal that the specific geometry of the crustal domain–with its extreme aspect ratio–requires an initial peak wavenumber from small-scale structures for the inverse cascade to occur. However, this same aspect ratio confines the cascade to structures on the scale of the crust, making the formation of a large-scale dipolar surface field unlikely. Despite these limitations, the inverse cascade remains a significant factor in the magnetic field evolution within the crust and may help explain highly magnetized objects with weak surface dipolar fields, such as low-field magnetars and central compact objects.