The AMBER experiment at CERN will measure the proton’s charge radius via muon-proton elastic scattering at high projectile energies and small momentum transfers to help to resolve the so-called proton radius puzzle, i.e., the discrepancy between charge radii measured with different experimental techniques. The core setup at AMBER consists of a hydrogen-filled time projection chamber (TPC). Tracking detectors upstream and downstream of the TPC measure the trajectories of the incoming and outgoing muons to determine their scattering angles. To resolve pile-up hits in the tracking detectors, we are constructing four high-granularity hodoscopes from 500-<math altimg="si1.svg" display="inline" id="d1e204"><mrow><mi mathvariant="normal">μ</mi><mi mathvariant="normal">m</mi></mrow></math> scintillating-plastic fibers and arrays of silicon photomultipliers. In this contribution, we present the design of the scintillating-fiber hodoscopes and first results of test-beam measurements with scaled-down prototypes. We will particularly emphasize how we managed to design detectors with a low material budget.

The combined use of the inverse kinematics technique and the advanced detection setup R³B (Reactions with Relativistic Radioactive Beams) at GSI/FAIR provides unique opportunities to study the fission process. This approach provides access to the complete isotopic identification of the two fission fragments, the precise determination of their velocities and the measurement of the neutrons and gammas emitted in coincidence, for a wide range of unstable fissile nuclei. In addition, quasi-free NN scattering represents a surrogate reaction to induce fission, allowing the complete identification of the fissioning system in terms of isotopic composition and excitation energy. The manuscript describes the technical realisation of these experiments as well as the physics programme and some preliminary results.

Neutron supermirrors are a crucial part of many scattering and particle physics experiments. So far, Ni(Mo)/Ti supermirrors have been used in experiments that require to transport a polarized neutron beam due to their lower saturation magnetization compared to Ni/Ti supermirrors. However, next generation <mml:math><mml:mi>β</mml:mi></mml:math> decay experiments require supermirrors that depolarize below <mml:math><mml:mrow><mml:mn>1</mml:mn><mml:msup><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mo>‑</mml:mo><mml:mn>4</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math> per reflection to reach their targeted precision. The depolarization of a polarized neutron beam due to reflection off Ni(Mo)/Ti supermirrors has not yet been measured to that precision. Recently, Cu/Ti supermirrors with a lower saturation magnetization compared to Ni(Mo)/Ti have been developed, and 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 curves of Ni(Mo) and Cu monolayers and <mml:math><mml:mrow><mml:mi>m</mml:mi><mml:mo>=</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:math> Ni(Mo)/Ti and Cu/Ti supermirrors measured at the neutron reflectometer SuperADAM and perform a full polarization analysis, with the aim to extract information about their magnetic moment. The results found, however, were inconclusive, since it seems a detection limit of this method for all measured samples was reached. At a second stage, we measured the depolarization (<mml:math><mml:mi>D</mml:mi></mml:math>) that a polarized neutron beam suffers after reflection off the same Ni(Mo)/Ti and Cu/Ti supermirrors by using the Opaque Test Bench setup. We find upper limits for the depolarization of <mml:math><mml:mrow><mml:msub><mml:mrow><mml:mi>D</mml:mi></mml:mrow><mml:mrow><mml:mtext>Cu/Ti(4N5)</mml:mtext></mml:mrow></mml:msub><mml:mo><</mml:mo><mml:mn>7</mml:mn><mml:mo>.</mml:mo><mml:mn>6</mml:mn><mml:mo>×</mml:mo><mml:mn>1</mml:mn><mml:msup><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mo>‑</mml:mo><mml:mn>5</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math>, <mml:math><mml:mrow><mml:msub><mml:mrow><mml:mi>D</mml:mi></mml:mrow><mml:mrow><mml:mtext>Ni(Mo)/Ti</mml:mtext></mml:mrow></mml:msub><mml:mo><</mml:mo><mml:mn>8</mml:mn><mml:mo>.</mml:mo><mml:mn>5</mml:mn><mml:mo>×</mml:mo><mml:mn>1</mml:mn><mml:msup><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mo>‑</mml:mo><mml:mn>5</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math>, and <mml:math><mml:mrow><mml:msub><mml:mrow><mml:mi>D</mml:mi></mml:mrow><mml:mrow><mml:mtext>Cu/Ti(2N6)</mml:mtext></mml:mrow></mml:msub><mml:mo><</mml:mo><mml:mn>6</mml:mn><mml:mo>.</mml:mo><mml:mn>0</mml:mn><mml:mo>×</mml:mo><mml:mn>1</mml:mn><mml:msup><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mo>‑</mml:mo><mml:mn>5</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math> at the <mml:math><mml:mrow><mml:mn>1</mml:mn><mml:mi>σ</mml:mi></mml:mrow></mml:math> confidence level, where (4N5) corresponds to a Ti purity of <mml:math><mml:mrow><mml:mn>99</mml:mn><mml:mo>.</mml:mo><mml:mn>995</mml:mn><mml:mspace></mml:mspace><mml:mstyle><mml:mi>%</mml:mi></mml:mstyle></mml:mrow></mml:math> and (2N6) to <mml:math><mml:mrow><mml:mn>99</mml:mn><mml:mo>.</mml:mo><mml:mn>6</mml:mn><mml:mspace></mml:mspace><mml:mstyle><mml:mi>%</mml:mi></mml:mstyle></mml:mrow></mml:math>. These results show that all three supermirrors are suitable for being used in next generation <mml:math><mml:mi>β</mml:mi></mml:math> decay experiments. We found no noticeable dependence of the depolarization on the <mml:math><mml:mi>q</mml:mi></mml:math> value or the magnetizing field, in which the samples were placed.

The diffuse supernova neutrino background (DSNB) created by stellar core-collapses throughout cosmic history is on the verge of discovery, with SK-Gd showing early deviations from the background expectation and JUNO starting to take data. However, the interpretation of early DSNB data will face significant challenges due to degeneracies between astrophysical parameters and uncertainties in supernova neutrino modeling. We explore how complementary astronomical observations can break these degeneracies and, in this context, we investigate whether early DSNB observations can constrain invisible supernovae, which have no optical emission but are powerful neutrino sources before being swallowed by a forming black hole. Leveraging the differences in the spectra between invisible and visible supernovae, we estimate the sensitivity of 1) detecting the existence of invisible supernovae, and 2) determining the fraction of invisible supernovae. Finally, we discuss how these conclusions depend on the spectral parameters of the black hole-forming component.

Dusty, submillimeter-selected galaxies without optical counterparts contribute a non-negligible fraction of the star formation in the early universe. However, such a population is difficult to detect through classical optical/UV-based surveys. We report the serendipitous discovery of such an optically dark galaxy, behind the quadruply-lensed $z=2.56$ quasar, H1413+117, offset to the north by 6\arcsec. From $^{12}$CO $J=4$--3, $J=6$--5, and part of the $J=13$--12 transitions, which all spatially coincide with a compact submillimeter continuum emission, we determine an unambiguous spectroscopic redshift, $z=3.386\pm 0.005$. This galaxy has a molecular mass $M_{\rm mol} \sim 10^{11}$ M$_\odot$ and a black hole mass $M_{\rm BH} \sim 10^{8}$ M$_\odot$, estimated from $^{12}$CO $J=4$--3 and archival {\it Chandra} X-ray data ($L_{\rm 2-10,keV} \sim 4 \times 10^{44}$\,erg\,s$^{-1}$), respectively. We also estimate a total infrared luminosity of $L_{\rm FIR} = (2.8\pm{2.3}) \times 10^{12}$ L$_\odot$ and a stellar mass of $M_* \lesssim 10^{11}$ M$_{\odot}$, from spectral energy distribution fitting. According to these simple mass estimations, this gas-rich and X-ray bright galaxy might be in a transition phase from starburst to quasar offering a unique case for studying galaxy-black hole co-evolution under extremely dusty conditions.

Aims. We investigate the James Webb Space Telescope (JWST) MIRI MRS gas molecular content of an externally irradiated Herbig disk, the F-type XUE 10 source, in the context of the eXtreme UV Environments (XUE) program. XUE 10 belongs to the massive star cluster NGC 6357 (1.69 kpc), where it is exposed to an external far-ultraviolet (FUV) radiation ≈10³ times stronger than in the solar neighborhood. Methods. We modeled the molecular features in the mid-infrared spectrum with local thermodynamic equilibrium (LTE) 0D slab models. We derived basic parameters of the stellar host from a VLT FORS2 optical spectrum using PHOENIX stellar templates. Results. We detected bright CO₂ gas with the first simultaneous detection (>5σ) of four isotopologues (¹²CO₂, ¹³CO₂, ¹⁶O¹²C¹⁸O, ¹⁶O¹²C¹⁷O) in a protoplanetary disk. We also detected faint CO emission (2σ) and the HI Pf α line (8σ). We placed strict upper limits on the water content, finding a total column density of ≲10¹⁸ cm^‑2. The CO₂ species trace low gas temperatures (300–370 K) with a range of column densities of 7.4 × 10¹⁷ cm^‑2 (¹⁶O¹²C¹⁷O)‑1.3 × 10²⁰ cm^‑2 (¹²CO₂) in an equivalent emitting radius of 1.15 au. The emission of ¹³CO₂ is likely affected by line optical depth effects. The ¹⁶O¹²C¹⁸O and ¹⁶O¹²C¹⁷O abundances may be isotopically anomalous compared to the ¹⁶O/¹⁸O and ¹⁶O/¹⁷O ratios measured in the interstellar medium and the Solar System. Conclusions. We propose that the mid-infrared spectrum of XUE 10 is explained by H₂O removal either via advection or strong photo-dissociation by stellar UV irradiation and enhanced local CO₂ gas phase production. Outer disk truncation supports the observed CO₂‑H₂O dichotomy. A CO₂ vapor enrichment in ¹⁸O and ¹⁷O can be explained by means of external UV irradiation and early (10^4–5 yr) delivery of isotopically anomalous water ice to the inner disk.

The obliquity between the stellar spin axis and the planetary orbit, detected via the Rossiter-McLaughlin (RM) effect, is a tracer of the formation history of planetary systems. While obliquity measurements have been extensively applied to hot Jupiters and short-period planets, they remain rare for cold and long-period planets due to observational challenges, particularly their long transit durations. We report the detection of the RM effect for the 19-hour transit of HIP 41378 f, a temperate giant planet on a 542-day orbit, observed through a worldwide spectroscopic campaign. We measured a slight projected obliquity of 21 ± 8° and a significant 3D spin-orbit angle of 52 ± 6°, based on the measurement of the stellar rotation period. HIP 41378 f is part a transiting system of five planets with planets close to mean motion resonances. The observed misalignment likely reflects a primordial tilt of the stellar spin axis relative to the protoplanetary disk, rather than dynamical interactions. HIP 41378 f is the first non-eccentric long-period planet (P>100 days) observed with the RM effect, opening new constraints on planetary formation theories. This observation should motivate the exploration of planetary obliquities across a longer range of orbital distances through international collaboration.

Cross-correlations techniques offer an alternative method to search for molecular species in JWST observations of exoplanet atmospheres. In a previous article, we applied cross-correlation functions for the first time to JWST NIRSpec/G395H observations of exoplanet atmospheres, resulting in a detection of CO in the transmission spectrum of WASP-39b and a tentative detection of CO isotopologues. Here we present an improved version of our cross-correlation technique and an investigation into how efficient the technique is when searching for other molecules in JWST NIRSpec/G395H data. Our search results in the detection of more molecules via cross-correlations in the atmosphere of WASP-39b, including $\rm H_{2}O$ and $\rm CO_{2}$, and confirms the CO detection. This result proves that cross-correlations are a robust and computationally cheap alternative method to search for molecular species in transmission spectra observed with JWST. We also searched for other molecules ($\rm CH_{4}$, $\rm NH_{3}$, $\rm SO_{2}$, $\rm N_{2}O$, $\rm H_{2}S$, $\rm PH_{3}$, $\rm O_{3}$ and $\rm C_{2}H_{2}$) that were not detected, for which we provide the definition of their cross-correlation baselines for future searches of those molecules in other targets. We find that that the cross-correlation search of each molecule is more efficient over limited wavelength regions of the spectrum, where the signal for that molecule dominates over other molecules, than over broad wavelength ranges. In general we also find that Gaussian normalization is the most efficient normalization mode for the generation of the molecular templates.

Generally, merger likelihood increases in denser environments; however, the large relative velocities at the centres of dense clusters are expected to reduce the likelihood of mergers for satellite galaxies. Tidal features probe the recent merger histories of galaxies. The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will produce an unprecedented sample of tidal features around millions of galaxies. We use LSST-like mock observations of galaxies at $z\sim0$ from the EAGLE, IllustrisTNG and Magneticum Pathfinder cosmological-hydrodynamical simulations to predict the occurrence rates of tidal features around satellite galaxies across group and cluster environments in the velocity-radius projected phase-space diagram to investigate the impact of these environments on tidal feature occurrence. We find that ancient infallers in the projected phase-space exhibit a decreasing tidal feature fraction with increasing halo mass, whereas recent infallers in the projected phase-space show unchanging tidal feature fractions with halo mass. Our results show, for the first time in cosmological simulations, a suppression of tidal feature fractions in the central regions of galaxy clusters, indicating a reduced merger rate due to higher cluster-centric velocities and lower galaxy total masses in the cluster centres. Using a toy model, we show that the presence of more tidal features in the recent infaller zone and cluster outskirts suggests that tidal features occur in interactions within infalling groups and dissipate by the time they are ancient infallers, indicating a $\lesssim3\pm2$ Gyr survival time of tidal features within clusters.

The substantial data volumes encountered in modern particle physics and other domains of fundamental physics research allow (and require) the use of increasingly complex data analysis tools and workflows. While the use of machine learning (ML) tools for data analysis has recently proliferated, these tools are typically special-purpose algorithms that rely, for example, on encoded physics knowledge to reach optimal performance. In this work, we investigate a new and orthogonal direction: Using recent progress in large language models (LLMs) to create a team of agents -- instances of LLMs with specific subtasks -- that jointly solve data analysis-based research problems in a way similar to how a human researcher might: by creating code to operate standard tools and libraries (including ML systems) and by building on results of previous iterations. If successful, such agent-based systems could be deployed to automate routine analysis components to counteract the increasing complexity of modern tool chains. To investigate the capabilities of current-generation commercial LLMs, we consider the task of anomaly detection via the publicly available and highly-studied LHC Olympics dataset. Several current models by OpenAI (GPT-4o, o4-mini, GPT-4.1, and GPT-5) are investigated and their stability tested. Overall, we observe the capacity of the agent-based system to solve this data analysis problem. The best agent-created solutions mirror the performance of human state-of-the-art results.

We study ultralight scalar fields with quadratic couplings to Standard-Model fermions and derive strong constraints from white-dwarf mass-radius data. Such couplings source scalar profiles inside compact stars, shift fermion masses, and can produce a new ground state of matter. We analyze couplings to electrons and to nucleons, incorporating composition and finite-temperature effects in white dwarf structure and equations of state. We identify two robust observables: (i) forbidden gaps - ranges of radii with no stable configurations - and (ii) characteristic shape distortions that drive white dwarf masses toward the Chandrasekhar limit (electron couplings) or shift the maximum mass (nucleon couplings). Confronting these predictions with precise measurements for Sirius B and Procyon B, together with the global white dwarf population, excludes large regions of unexplored parameter space and extends earlier QCD-axion-specific bounds to a broader class of scalar theories. Our stellar constraints rely only on sourcing and do not assume the scalar constitutes dark matter; where mass reductions are small, precision laboratory searches remain competitive. White-dwarf astrophysics thus provides a powerful, largely assumption-minimal probe of ultralight, quadratically coupled scalars.

Multiple photometric studies have reported the presence of seemingly older accreting pre-main sequence stars (PMS) in optical colour-magnitude diagrams (CMDs). We investigate this phenomenon in the Orion Nebula, which harbors a subset of stars that show infrared excess detected by Spitzer and Halpha excess emission, yet display significantly older isochronal ages (>10 Myr) compared to the bulk population (~1-3 Myr) in the r, (r-i) CMD. We perform a detailed spectroscopic analysis of 40 Orion Nebula stars using VLT/X-Shooter, covering CMD-based isochronal ages from 1 to over 30 Myr. We derive extinction values, stellar properties, and accretion parameters by modeling the ultraviolet excess emission through a multicomponent fitting procedure. The sample spans spectral types from M4.5 up to K6, and masses in the range ~0.1-0.8 Msun. We demonstrate that, when extinction and, most importantly, accretion effects are accurately constrained, the stellar luminosity and effective temperature of the majority of the seemingly old stars become consistent with a younger population (~1-5 Myr). This is supported by strong lithium absorption, which corroborates their youth, and by the accretion-to-stellar luminosity ratios typical for young, accreting stars. Three of these sources, however, remain old even after our analysis, despite showing signatures consistent with ongoing accretion from a protoplanetary disc. More generally, our analysis indicates that excess continuum emission from accretion shocks affects the placement of PMS stars in the CMD, displacing sources towards bluer optical colours. This study highlights the critical role of accretion in shaping the stellar properties estimates (including age) derived from optical CMDs and emphasizes the need to carefully account for accretion effects when interpreting age distributions in star-forming regions.

Under the core-accretion model, gas giants form via runaway accretion. This process starts when the mass of the accreted envelope becomes equal to the mass of the core. Here, we model a population of warm sub-Saturns to search for imprints of their formation history in their internal structure. Using the GAS gianT modeL for Interiors (GASTLI), we calculate a grid of interior structure models on which we perform retrievals for our sample of 28 sub-Saturns to derive their envelope mass fractions ($f_{env}$). For each planet, we run three different retrievals assuming low (-2.0 < log(Fe/H) < 0.5), medium ( 0.5 < log(Fe/H) < 1.4), and high (1.4 < log(Fe/H) < 1.7) atmospheric metallicity. The distribution of $f_{env}$ in our sample is then compared to predictions of planet formation models. When compared to the outcomes of a planetesimal accretion model, we find that we require medium to high atmospheric metallicities to reproduce the simulated planet population. Additionally, we find a bimodal distribution of $f_{env}$ in our sample with a gap that is located at different values of $f_{env}$ for different atmospheric metallicities. For the high atmospheric metallicity case, the gap in the $f_{env}$ distribution is located between 0.5 and 0.7, which is consistent with assumptions by the core-accretion model where runaway accretion starts when $M_{env} \approx M_{core}$ ($f_{env} \sim 0.5$). We also find a bimodal distribution of the hydrogen and helium mass fraction ($f_{H/He}$) with a gap at $f_{H/He} = 0.3$. The location of this gap is independent of the assumed atmospheric metallicity. Lastly, we compare the distributions of our sub-Saturns in the Neptunian savanna to a population of sub-Saturns in the Neptune desert and ridge. We find that the observed $f_{env}$ distribution of savanna and ridge sub-Saturns is consistent with the planets coming from the same underlying population.

Planetesimal formation likely lasted for millions of years in the Solar nebula, and the cold classicals in the Kuiper belt are suggested to be the direct products of streaming instability. The presence of minor planetary bodies in the outer Solar System and the exo-Kuiper belts provide key constraints to planet formation models. In this work, we connected dust drift and coagulation, planetesimal formation, N-body gravity, pebble accretion, planet migration, planetary core accretion, gap opening, and internal photoevaporation in one modeling framework. We demonstrate that multiple classes of minor planets, or planetesimals, can form during disk dissipation and remain afterwards, including a scattered group, a resonant group and a dynamically cold group. Significant growth by pebble accretion was prevented by both dynamical heating due to the giant planet in the system and rapid dispersal of the disk towards the end of its lifetime. We also conducted a parameter study which showed that this is not a universal case, where the outcome is determined by the competition for dust between planetesimal formation and pebble accretion. Combining this scenario with sequential planet formation, this model provides a promising pathway towards an outer Solar System formation model.

We study magnetic conversion of ultra-relativistic axion-like particles (ALPs) into photons in compact-star environments, focusing on the hot, transient conditions of core-collapse supernova (SN) remnants and neutron-star mergers (NSMs). We address previously overlooked uncertainties, particularly the suppression caused by ejected matter near the stellar surface, a region crucial to the conversion process. We derive analytical expressions for the transition rate; they reveal the influence of key parameters and their uncertainties. We update constraints using historical gamma-ray data from SN~1987A and find $g_{aγ}<5\times10^{-12}~{\rm GeV}^{-1}$ for $m_a\lesssim10^{-9}$ meV. We also forecast sensitivities for a future Galactic SN and for NSMs, assuming observations with Fermi-LAT or similar gamma-ray instruments. We distinguish ALPs -- defined as coupling only to photons and produced via Primakoff scattering -- from axions, which also couple to nucleons and emerge through nuclear bremsstrahlung. We omit pionic axion production due to its large uncertainties and inconsistencies, though it could contribute comparably to bremsstrahlung under optimistic assumptions. For the compact sources, we adopt time-averaged one-zone models, guided by numerical simulations, to enable clear and reproducible parametric studies.

Dense young star clusters (YSCs) are ideal environments for dynamical interactions between stars and stellar mass black holes (BHs). In such dense environments, stars can undergo close encounters with BHs and fall within their tidal radius, resulting in micro-tidal disruption events (micro-TDEs), transient phenomena with potential multi-messenger signatures. We performed a suite of direct N-body simulations using the PETAR code, to which we implemented new prescriptions for modeling micro-TDEs. We constructed a set of realistic YSC models including primordial binaries, based on the observed Milky Way population. Our simulations incorporate stellar and binary evolution, supernova kicks, and stellar winds using the BSE code, and account for the Galactic tidal field via the GALPY library. We identify three dynamical channels for micro-TDE production: single star-single BH encounters, binary-mediated interactions (including supernova-kick triggers), and interactions involving higher-order multiple systems such as hierarchical triples and quadruples, as well as chaotic few-body interactions. Multiple encounters are the most efficient production channel, which dominates the total rate: 350-450 Gpc$^{-3}$ yr$^{-1}$. Micro-TDEs from YSCs are expected to be detectable by upcoming surveys, particularly the Legacy Survey of Space and Time, with detection rates potentially up to hundreds per year. The gravitational wave (GW) signals expected from micro-TDEs peak in the deci-Hertz band, making them accessible to future instruments such as the Lunar Gravitational Wave Antenna and the Deci-Hertz Interferometer Gravitational-wave Observatory. Micro-TDEs emerge as promising multi-messenger sources, potentially offering unique insights into star cluster dynamics, stellar collisions, and the population of dormant stellar-mass BHs, through both electromagnetic and GW observations.

As the mathematical capabilities of large language models (LLMs) improve, it becomes increasingly important to evaluate their performance on research-level tasks at the frontier of mathematical knowledge. However, existing benchmarks are limited, as they focus solely on final-answer questions or high-school competition problems. To address this gap, we introduce IMProofBench, a private benchmark consisting of 39 peer-reviewed problems developed by expert mathematicians. Each problem requires a detailed proof and is paired with subproblems that have final answers, supporting both an evaluation of mathematical reasoning capabilities by human experts and a large-scale quantitative analysis through automated grading. Furthermore, unlike prior benchmarks, the evaluation setup simulates a realistic research environment: models operate in an agentic framework with tools like web search for literature review and mathematical software such as SageMath. Our results show that current LLMs can succeed at the more accessible research-level questions, but still encounter significant difficulties on more challenging problems. Quantitatively, Grok-4 achieves the highest accuracy of 52% on final-answer subproblems, while GPT-5 obtains the best performance for proof generation, achieving a fully correct solution for 22% of problems. IMProofBench will continue to evolve as a dynamic benchmark in collaboration with the mathematical community, ensuring its relevance for evaluating the next generation of LLMs.

One of the most puzzling properties of the high-redshift AGN population recently discovered by JWST, including both broad-line and narrow-line sources, is their X-ray weakness. With very few exceptions, and regardless of the optical classification, they are undetected at the limits of the deepest Chandra fields, even when stacking signals from tens of sources in standard observed-frame energy intervals (soft, hard, and full bands). It has been proposed that their elusive nature in the X-ray band is due to heavy absorption by dust-free gas or intrinsic weakness, possibly due to high, super-Eddington accretion. In this work, we perform X-ray stacking in three customized rest-frame energy ranges (1-4, 4-7.25, and 10-30 keV) of a sample of 50 Type 1 and 38 Type 2 AGN identified by JWST in the CDFS and CDFN fields. For the Type 2 sub-sample, we reach a total of about 210 Ms exposure, and we report a significant ($\sim 3σ$) detection in the hardest (10-30 keV rest frame) band, along with relatively tight upper limits in the rest frame softer energy bands. The most straightforward interpretation is in terms of heavy obscuration due to gas column densities well within the Compton thick regime ($> 2 \times 10^{24} $cm$^{-2}$) with a large covering factor, approaching 4$π$. The same procedure applied to the Type 1 sub-sample returns no evidence for a significant signal in about 140 Ms stacked data in any of the adopted bands, confirming their surprisingly elusive nature in the X-ray band obtained with previous stacking experiments. A brief comparison with the current observations and the implications for the evolution of AGN are discussed.

We further scrutinize the evidence for a recently suggested pseudo-scalar particle, the electroweak $η_{\rm w}$-meson. Its existence is demanded by matching the removal of the weak vacuum angle $θ_{\rm w}$ by the anomalous $B+L$ - symmetry with a massive pole in the topological susceptibility of the vacuum. We specifically focus on the possibility of the emergence of $η_{\rm w}$ as a collective excitation of the phase of the condensate of the 't Hooft fermion determinant, generated by the electroweak instantons, which breaks the $B+L$ - symmetry spontaneously. We argue that the generation of the 't Hooft vertex is in one-to-one correspondence with its non-zero vacuum expectation value which is cutoff insensitive. We outline certain puzzles about the nature of the emergent $η_{\rm w}$ which require further investigations.

We embed the perturbative Fock state of the Yang-Mills BV-multiplet in the vertex operator algebra of the path-integral for the $\mathcal{N}=2$ supersymmetric world line and evaluate the pull-back of the latter to an integral form on supermoduli space. Choosing a suitable Poincaré dual on the latter, we show that this integral form describes an extension of Yang-Mills theory. Upon projection back to the Fock space, we recover the Yang-Mills action from the world line. This furthermore gives an a priori justification for the construction of Yang-Mills equations of motion as emerging from deformations of the BRST differential.

Luminosities of pre-main sequence stars evolve during the protoplanetary disc lifetime. This has a significant impact on the heating of their surrounding protoplanetary disks, the natal environments of planets. Moreover, stars of different masses evolve differently. However, this is rarely accounted for in planet formation models. We carry out pebble-driven core accretion planet formation modelling with focus on the impact of pre-main sequence stellar luminosity evolution on giant planet formation around host stars in the range of <inline-formula><tex-math id="TM0001" notation="LaTeX">$1{-}2.4\ \rm M_{\odot }$</tex-math></inline-formula>. We find that giant planet formation is sensitive to the evolution of stellar luminosity, specifically the locations and times at which giant planet formation can occur depend on it. High stellar luminosity causes an increase in the scale height of the gas and pebbles, which may decrease the efficiency of pebble accretion, making it more challenging to form giant planets. This has important consequences for the composition of these giant planets, stressing the need to incorporate such aspects into planet formation models.

Context. Our knowledge of the initial conditions of terrestrial planet formation is mainly based on the study of protoplanetary disks around nearby isolated low-mass stars. However, most young stars and therefore planetary systems form in high-mass star-forming regions and are exposed to ultraviolet radiation, affecting the protoplanetary disk. These regions are located at large distances and only now with JWST has it become accessible to study the inner disks surrounding young stars. Aims. We present the eXtreme UV Environments (XUE) program, which provides the first detailed characterization of the physical and chemical properties of the inner disks around young intermediate-mass (1–4 M_⊙) stars exposed to external irradiation from nearby massive stars. We present high-signal-to-noise MIRI-MRS spectroscopy of 12 disks located in three subclusters of the high-mass star-forming region NGC 6357 (d ~ 1690 pc). Methods. Based on their mid-infrared spectral energy distribution, we classified the XUE sources into Group I and II based on the Meeus scheme. We analyzed their molecular emission features, and compared their spectral indices and 10 μm silicate emission profiles to the ones of nearby Herbig and intermediate T Tauri (IMTT) disks. Results. The XUE program provides the first detailed characterization of the rich molecular inventory in IMTT disks, including water, CO, CO₂, HCN, and C₂H₂. In the XUE sample, the detected emission likely originates from within 10 au, although this inner disk origin may not be typical for all externally irradiated disks. Despite being more massive, the XUE stars host disks with a molecular richness comparable to isolated T Tauri systems. The spectral indices are also consistent with similar-mass stars in nearby regions. The 10 μm silicate features in the XUE sample exhibit lower F_11.3/F_9.8 ratios at a given F_peak, suggesting that the disk surfaces may be dominated by smaller grains compared to nearby disks. However, uncertainties in extinction prevent us from drawing firm conclusions about their inner disk properties. The majority of disks display water emission from the inner disk, suggesting that even in these extreme environments rocky planets can form in the presence of water. Only one object shows PAH emission, contrasting with the higher PAH detection rates in IMTT surveys from lower-UV environments. Conclusions. The absence of strong line fluxes and other irradiation signatures suggests that the XUE disks have been truncated by external UV photons. However, this truncation does not appear to significantly impact the chemical richness of their inner regions. These findings indicate that even in extreme environments, IMTT disks can retain the ingredients necessary for rocky planet formation, comparable to the ones of lower-mass T Tauri disks in low-mass star-forming regions.

We investigate the influence of supernova (SN) feedback on the satellites of elliptical host galaxies using hydrodynamic simulations. Utilizing a modified version of the GADGET-3 code, we perform cosmological zoom-in simulations of 11 elliptical galaxies with stellar masses in the range 10¹¹M_⊙ < M_* < 2 × 10¹¹M_⊙. We conduct two sets of simulations with identical initial conditions: a fiducial model, which includes a three-phase SN mechanical wind, and a weak SN feedback model, where nearly all SN energy is released as thermal energy with a reduced SN wind velocity. Our comparison shows minimal differences in the elliptical host galaxies, but significant variations in the physical properties of satellite galaxies. The weak SN feedback model produces a larger number of satellite galaxies compared to the fiducial model, and significantly more than observed. For satellite galaxies with stellar masses above 10⁸M_⊙, the weak SN feedback model generates approximately 5 times more satellites than observed in the Extending the Satellites Around Galactic Analogs Survey (or xSAGA) survey. Most of these overproduced satellites have small stellar masses, below 10¹⁰M_⊙. Additionally, satellites in the weak SN feedback model are about 3.5 times more compact than those observed in the SAGA survey and the fiducial model, with metallicities nearly 1 dex higher than observed values. In conclusion, the satellite galaxies in the fiducial model, which includes mechanical SN feedback, exhibit properties that more closely align with observations. This underscores the necessity of incorporating both mechanical active galactic nuclei and SN feedback to reproduce the observed properties of elliptical galaxies and their satellites in simulations.

(Sub-)millimetre single-dish telescopes feature faster mapping speeds and access larger spatial scales than their interferometric counterparts. However, atmospheric fluctuations tend to dominate their signals and complicate recovery of the astronomical sky. Here we develop a framework for Gaussian process-based sky reconstruction and separation of the atmospheric emission from the astronomical signal based on Numerical Information Field Theory (\texttt{NIFTy}). To validate this novel approach, we use the \textit{maria} software to generate synthetic time-ordered observational data mimicking the MUSTANG-2 bolometric array. This approach leads to significantly improved sky reconstructions versus traditional methods.

Multiply-imaged supernovae (SNe) provide a novel means of constraining the Hubble constant ($H_0$). Such measurements require a combination of precise models of the lensing mass distribution and an accurate estimate of the relative time delays between arrival of the multiple images. Only two multiply-imaged SNe, Refsdal and H0pe, have enabled measurements of $H_0$ thus far. Here we detail the third such measurement for SN Encore, a $z=1.95$ SNIa discovered in JWST/NIRCam imaging. We measure the time delay, perform simulations of additional microlensing and millilensing systematics, and combine with the mass models of Suyu et al. in a double-blind analysis to obtain our $H_0$ constraint. Our final time-delay measurement is $Δt_{1b,1a}=-39.8_{-3.3}^{+3.9}$ days, which is combined with seven lens models weighted by the likelihood of the observed multiple image positions for a result of $H_0=66.9_{-8.1}^{+11.2} \rm{km} \rm{s}^{-1}\rm{Mpc}^{-1}$. The uncertainty on this measurement could be improved significantly if template imaging is obtained. Remarkably, a sibling to SN Encore (SN "Requiem") was discovered in the same host galaxy, making the MACS J0138.0-2155 cluster the first system known to produce more than one observed multiply-imaged SN. SN Requiem has a fourth image that is expected to appear within a few years, providing an unprecedented decade-long baseline for time-delay cosmography and an opportunity for a high-precision joint estimate of $H_0$.

GRB 250702B is an exceptional transient that produced multiple episodes of luminous gamma-ray radiation lasting for $>25$ ks, placing it among the class of ultra-long gamma-ray bursts (GRBs). However, unlike any known GRB, a soft X-ray precursor was discovered by the Einstein Probe up to 24 hours before the gamma-ray triggers. We present comprehensive X-ray observations of the transient's afterglow obtained with the Neil Gehrels Swift Observatory, the Nuclear Spectroscopic Telescope Array, and the Chandra X-ray Observatory between 0.5 to 65 days (observer frame) after the initial high-energy trigger. The X-ray emission decays steeply as $\sim t^{-1.9}$, and shows short timescale X-ray variability ($ΔT/T < 0.03$) in both Swift and NuSTAR, consistent with flares superposed on an external shock continuum. Serendipitous detections by the Swift Burst Alert Telescope (BAT) out to $\sim$0.3 days and continued NuSTAR variability to $\sim$2 days imply sustained central engine activity; including the precursor, the required engine duration is $\gtrsim 3$ days. Afterglow modeling favors the combination of forward and reverse shock emission in a wind-like ($k \approx 2$) environment. These properties, especially the long-lived engine and soft X-ray precursor, are difficult to reconcile with a collapsar origin, and GRB 250702B does not fit neatly with canonical ultra-long GRBs or relativistic tidal disruption events (TDEs). A hybrid scenario in which a star is disrupted by a stellar-mass black hole (a micro-TDE) provides a plausible explanation, although a relativistic TDE from an intermediate-mass black hole remains viable. Decisive discrimination between progenitors will require sensitive late-time X-ray observations.

The polarized light of the cosmic microwave background is sensitive to new physics that violates parity symmetry. For example, the interaction of photons with the fields of elusive dark matter and dark energy could cause a uniform rotation of the plane of linear polarization across the sky, an effect known as cosmic birefringence. We extract the cosmological rotation angle, $β$, using Bayesian analysis of parity-violating correlations, $EB$ and $TB$, of polarization data from the Atacama Cosmology Telescope (ACT) Data Release 6. We use prior probabilities for instrumental miscalibration angles derived from the optics model for the ACT telescope and instruments, and marginalize over a residual intensity-to-polarization leakage. We measure $β= 0.215^\circ\pm 0.074^\circ$ (68\% confidence level), which excludes $β=0$ with a statistical significance of $2.9σ$. Although there remain systematics in the ACT data that are not understood and do not allow us to draw strong cosmological conclusions, this result is consistent with previous independent results from the \wmap\ and \planck\ missions. It is suggestive that independent data sets and analyses using different methodologies have yielded the same sign and comparable magnitudes for $β$.

Earthshine observations offer a unique opportunity to study Earth as an exoplanet seen from the Moon. As the Sun-Earth-Moon geometry changes, Earth can be observed as a spatially unresolved exoplanet at different phase angles, providing important context for future observations of Earth-like exoplanets. Here, we present a catalog of Earthshine polarization spectra obtained with FORS2 on the VLT, covering diverse scenes, surface conditions, cloud properties, and weather patterns for over a decade. For the first time, we model this extensive dataset in detail using a homogeneous modeling framework. Previous efforts to model some of these spectra struggled to reproduce the observed polarization continuum, even with advanced 3D radiative transfer models incorporating satellite-derived surface and atmospheric data. We improve upon this with a 3D model that includes subgrid cloud variability, wavelength-dependent surface albedo maps, and an accurate treatment of ocean glint. Our simulations successfully reproduce most observed spectra to a much higher precision than previously possible. Our statistical analysis reveals that the spectral slope in the visible can distinguish between ocean and mixed surfaces in both reflected and polarized light, which is not possible using broadband filters alone. Polarized light at large phase angles, beyond the Rayleigh scattering regime, is particularly effective in differentiating oceans from land, unlike reflected light. We also identify correlations between cloud optical thickness and the polarized spectral slope, and between cloud cover and broadband B-R differences in reflected light, demonstrating the diagnostic power of these observations. This work highlight the potential of polarization for characterizing Earth-like exoplanets. From polarization alone, we can infer oceans, vegetation, and an active water cycle, key indicators of a habitable planet.

GRB 250702B was the longest gamma-ray burst ever observed, with a duration that challenges standard collapsar models and suggests an exotic progenitor. We collected a rich set of optical and infrared follow-up observations of its rapidly fading afterglow using a suite of telescopes including the W. M. Keck Observatory, the Gemini telescopes, the Magellan Baade Telescope, the Victor M. Blanco 4-meter telescope, and the Fraunhofer Telescope at Wendelstein Observatory. Our analysis reveals that the afterglow emission is well described by forward shock emission from a highly obscured relativistic jet. Deep photometric observations of the host galaxy reveal a massive 10^10.66 solar masses, dusty, and extremely asymmetric system that is consistent with two galaxies undergoing a major merger. The galactocentric offset, host galaxy properties, and jet characteristics do not definitively distinguish between competing progenitor scenarios. We find that the afterglow and host are consistent with a range of progenitors including a collapsar, a merger between a helium star and a stellar mass black hole, the disruption of a star by a stellar mass compact object, and the tidal disruption of a star by an off-nuclear intermediate mass black hole.

Aims. We investigate the role of cosmic ray (CR) halos in shaping the physical properties of starburst-driven galactic outflows. Methods. We constructed a model for galactic outflows driven by a continuous central injection of energy, gas, and CRs, where the treatment of CRs accounts for the effect of CR pressure gradients on the flow dynamics. The model parameters were set by the effective properties of a starburst. By analyzing the asymptotic behavior of our model, we derived the launching criteria for starburst-driven galactic outflows and determined their corresponding outflow velocities. Results. We find that in the absence of CRs, stellar feedback can only launch galactic outflows 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. Outflows driven by CRs dominate in systems with SFR surface densities below the critical threshold, but their influence diminishes in highly star-forming systems. However, in older systems with established CR halos, the CR contribution to outflows weakens once the outflow reaches the galactic scale height, making CRs ineffective in sustaining outflows in such environments. 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 yet to build significant CR halos. In contrast, fast outflows in starburst galaxies, where the SFR surface density exceeds the critical threshold, are primarily driven by thermal energy and remain largely unaffected by CR halos.

We utilize the Magneticum suite of hydrodynamical simulations to investigate the formation and evolution of cosmic voids from $z = 5.04$ to present day, using cold dark matter and (sub-) halo tracers in high-density samples. This includes the evolution of their global properties, such as size, shape, inner density, and average density, as well as their radial density profiles. Our results provide several key conclusions for void analyses in modern surveys. We demonstrate that a relative size framework is required, mitigating methodological selection effects and revealing the true physical evolution of densities around halo-defined voids. This necessity arises from our findings that a void's properties are more fundamentally tied to its rank within its contemporary population than to its absolute size. Using this framework, we show that the evolution of halo voids stabilizes at redshifts below $z \simeq 1$, driven primarily by cosmic expansion rather than ongoing halo formation. We further find that the matter evolution around these stable voids is remarkably well-described by linear growth theory, with deviations appearing as non-linear growth on small scales and suppressed growth in the largest voids, potentially driven by the influence of dark energy. This late-time stability and the predictable evolution confirm voids as pristine laboratories for probing the nature of dark energy with upcoming surveys.

The core of a massive star (M > 8 Msun) eventually collapses. This implosion usually triggers a supernova (SN) explosion that ejects most of the stellar envelope and leaves behind a neutron star (NS) with a mass of up to about 2 Msun. Sometimes the explosion fails and a black hole forms instead. The NS radiates its immense binding energy (some 10% of its rest mass or $2-4\times10^{53}$ erg) almost entirely as neutrinos and antineutrinos of all flavors with typical energies of some 10 MeV. This makes core-collapse SNe the most powerful neutrino factories in the Universe. Such a signal was observed once - with limited statistics - from SN 1987A in the Large Magellanic Cloud. Today, however, many large neutrino detectors act as SN observatories and would register a high-statistics signal. A future Galactic SN, though rare (1-3 per century), would produce a wealth of astrophysical and particle-physics information, including possible signatures for new particles. Neutrinos are key to SN dynamics in the framework of the Bethe-Wilson delayed explosion paradigm. After collapse, they are trapped in the core for a few seconds, forming a dense neutrino plasma that can exhibit collective flavor evolution caused by the weak interaction, a subject of intense theoretical research.

Supermassive black hole feedback is the currently favoured mechanism to regulate the star formation rate of galaxies and prevent the formation of ultra-massive galaxies (M_⋆ > 10¹² M_⊙). However, the mechanism through which the outflowing energy is transferred to the surrounding medium strongly varies from one galaxy evolution model to another, such that a unified model for active galactic nucleus (AGN) feedback does not currently exist. The hot atmospheres of galaxy groups are highly sensitive laboratories of the feedback process, as the injected black hole energy is comparable to the binding energy of halo gas particles. Here we report multi-wavelength observations of the fossil galaxy group SDSSTG 4436. The hot atmosphere of this system exhibits a highly relaxed morphology centred on the giant elliptical galaxy NGC 3298. The X-ray emission from the system features a compact core (< 10 kpc) and a steep increase in the entropy and cooling time of the gas, with the cooling time reaching the age of the Universe ∼15 kpc from the centre of the galaxy. The observed entropy profile implies a total injected energy of ∼1.5 × 10⁶¹ ergs, which given the high level of relaxation could not have been injected by a recent merging event. Star formation in the central galaxy NGC 3298 is strongly quenched and its stellar population is very old (∼10.6 Gyr). The currently detected radio jets have low power and are confined within the central compact core. All the available evidence implies that this system was affected by giant AGN outbursts that raised the entropy of the neighbouring gas to the point that the gas no longer efficiently cools. Our findings imply that AGN outbursts can be energetic enough to unbind gas particles and lead to the disruption of cool cores.

We present new observations that densely sample the microwave (4-360 GHz) continuum spectra from eight young systems in the Taurus region. Multi-component, empirical model prescriptions were used to disentangle the contributions from their dust disks and other emission mechanisms. We found partially optically thick, free-free emission in all these systems, with positive spectral indices (median at 10 GHz) and contributing 5-50% of the 43 GHz fluxes. There is no evidence for synchrotron or spinning dust grain emission contributions for these targets. The inferred dust disk spectra all show substantial curvature: their spectral indices decrease with frequency, from -4.0 around 43 GHz to 1.7-2.1 around 340 GHz. This curvature suggests that a substantial fraction of the (sub)millimeter ( ≳ 200 GHz) dust emission may be optically thick, and therefore the traditional metrics for estimating dust masses are flawed. Assuming the emission at lower frequencies (43 GHz) is optically thin, the local spectral indices and fluxes were used to constrain the disk-averaged dust properties and estimate corresponding dust masses. These masses are roughly an order of magnitude higher ( ≈1000M⊕) than those found from the traditional approach based on (sub)millimeter fluxes. These findings emphasize the value of broad spectral coverage - particularly extending to lower frequencies ( ∼cm-band) - for accurately interpreting dust disk emission; such observations may help reshape our perspective on the available mass budgets for planet formation.

Broad absorption line (BAL) quasars are often considered X-ray weak relative to their optical/UV luminosity, whether intrinsically (i.e. the coronal emission is fainter) or due to large column densities of absorbing material. The SDSS-V is providing optical spectroscopy for samples of quasar candidates identified by eROSITA as well as Chandra, XMM, or Swift, making the resulting data sets ideal for characterizing the BAL quasar population within an X-ray selected sample. We use the Balnicity Index (BI) to identify the BAL quasars based on absorption of the C IV<inline-formula><tex-math>$\lambda \, 1549$</tex-math></inline-formula> emission line in the optical spectra, finding 143 BAL quasars in our sample of 2317 X-ray selected quasars within <inline-formula><tex-math>$1.5\le z \le 3.5$</tex-math></inline-formula>. This observed BAL fraction of <inline-formula><tex-math>$\approx$</tex-math></inline-formula> 6 per cent is comparable to that found in optically selected samples. We also identify absorption systems via the Absorption Index (AI) which includes mini-BALs and NALs, finding 954 quasars with AI <inline-formula><tex-math>$>0$</tex-math></inline-formula>. We consider the C IV emission space (equivalent width versus blueshift) to study the BAL outflows within the context of the radiatively driven accretion disc–wind model. X-ray selection excludes the highest outflow velocities in emission but includes the full range of absorption velocities which we suggest is consistent with the BAL gas being located further from the X-ray corona than the emitting gas. We observe both X-ray weak and X-ray strong BALs (via the optical-to-X-ray spectral slope, <inline-formula><tex-math>$\alpha _\text{ox}$</tex-math></inline-formula>) and detect little evidence for differing column densities between the BAL and non-BAL quasars, suggesting the BALs and non-BALs have the same shielding gas and intrinsic X-ray emission.

The goal of this work is to estimate the Hubble constant <inline-formula><tex-math>$H_0$</tex-math></inline-formula> through the time-delay cosmographic study of the quadruply lensed, variable quasi-stellar objects (QSO) SDSSJ1433+6007. We combine multifilter, archival Hubble Space Telescope data for lens modelling with a dedicated 3-yr long time-delay monitoring campaign using the 2.1 m Fraunhofer telescope at the Wendelstein Observatory. The lens modelling is performed with the public LENSTRONOMY PYTHON package individually for the infrared data, utilizing the higher resolution of the optical data to constrain image positions a priori. This approach revealed two luminous contaminants in one of the near-infrared exposures, which would bias the lensing potentials and cosmological inference if left unaccounted. After masking these contaminants, we repeated the modelling and combined the lens posteriors, obtaining a constraint on the Fermat potential with a statistical uncertainty of <inline-formula><tex-math>$2.6\, {{\rm per\ cent}}$</tex-math></inline-formula>. The g'-band Wendelstein light-curve data are reduced and then analysed using a free-knot spline fitting method implemented in the public PYTHON PYCS3 tools, accounting for microlensing correction. We obtain a precision of <inline-formula><tex-math>$6.5\, {{\rm per\ cent}}$</tex-math></inline-formula> for the time delays between the QSO images. By combining the posteriors for the Fermat potential differences and time delays, and assuming a flat Lambda-cold dark matter cosmology, we infer a Hubble constant of <inline-formula><tex-math>$H_0=71.7^{+3.9}_{-3.6}\, {\rm{km}}\,{\mathrm{Mpc}^{-1}~\mathrm{s}^{-1}}$</tex-math></inline-formula>, achieving <inline-formula><tex-math>$5.3\, {{\rm per\ cent}}$</tex-math></inline-formula> purely statistical uncertainty for this single system. Complementary observations and further study are required to address the systematic errors fully.

Context. Supergiant B[e] (sgB[e]) stars are exceptionally rare objects, with only a select number of confirmed examples in the Milky Way. The evolutionary pathways leading to the sgB[e] phase remain largely debated, highlighting the need for additional observations. The sgB[e] star Wd1-9, located in the massive cluster Westerlund 1 (Wd1), is enshrouded in a dusty cocoon – likely the result of past eruptive activity – leaving its true nature enigmatic. Aims. We present the most detailed X-ray study of Wd1-9 to date, using X-rays that pierce through its cocoon with the aim of uncovering its nature and evolutionary state. Methods. We utilised 36 Chandra observations of Wd1 from the 'Extended Westerlund 1 and 2 Open Clusters Survey' (EWOCS), plus eight archival datasets, totalling 1.1 Ms. We used this dataset to investigate long-term variability and periodicity in Wd1-9, and analysed X-ray colours and spectra over time to uncover patterns that shed light on its nature. Results. Wd1-9 exhibits significant long-term X-ray variability, within which we identify a strong ∼14-day periodic signal. We interpret this as the orbital period, marking the first period determination for the system. The X-ray spectrum of Wd1-9 is thermal and hard (kT ∼ 3.0 keV), resembling the spectra of bright Wolf-Rayet (WR) binaries in Wd1, while a strong Fe emission line at 6.7 keV indicates hot plasma from a colliding-wind X-ray binary. Conclusions. Wd1-9, with evidence of past mass loss, circumbinary material, a hard X-ray spectrum, and a newly detected 14-day period, displays all the hallmarks of a binary – likely a WR+OB – that recently underwent early Case B mass transfer. Its sgB[e] classification is likely phenomenological, reflecting emission from the dense circumbinary material. This places Wd1-9 in a rarely observed phase, possibly revealing a newly formed WN star, bridging the gap between immediate precursors and later evolutionary stages in Wd1.

Based on the scale-free nature of gravity, the structure in the Universe is expected to be self-similar on large scales. However, this self-similarity eventually breaks down due to small-scale gas physics such as star formation, active galactic nucleus (AGN) and stellar feedback, and non-linear effects gaining importance relative to linear structure formation. In this work, we investigate the large-scale matter flows that connect collapsed structures to their cosmic environments. Specifically, we focus on their agreement with self-similarity in various properties. For this purpose we used the full power of the hydrodynamical cosmological simulation suite Magneticum Pathfinder to precisely calculate the instantaneous inflow and outflow rates of structures on a large range of masses and redshifts. We find a striking self-similarity across the whole mass range and through time that only breaks down in the outflowing regime due to the different outflow driving mechanisms for galaxies versus galaxy clusters. We additionally performed a geometrical analysis of the patterns of inflow versus outflow to demonstrate how the inflows organise into anisotropic filaments driven by the tidal distortions of the environment, while the outflows are fairly isotropic due to their thermal nature. This also manifests in the differences in the thermal and chemical properties of the gas in the inflowing and outflowing component: While the inflowing gas is pristine and colder, encountering the accretion shock surfaces and entering the influence region of AGN and stellar feedback heats the gas up into a diffuse metal-enriched hot atmosphere. Overall the differences between outflowing and infalling gas are enhanced at the galaxy cluster scale compared to the galaxy scale due to the strong accretion shocks that reach out to large radii for these objects. An individual study of the gas motions in the outskirts of one of the most massive clusters in the simulations we carried out demonstrates these results to greater detail: Gas found in the outer (r > 1.2r_vir) hot atmosphere at z = 0 falls in and is completely enriched early in the assembly process before being shock heated and expanding.

Early stages of stellar birth comprise of a two-step process involving the formation of two hydrostatic cores. The second step of gravitational collapse sets the radiative efficiency and accretion rate of the young protostar. These two parameters, of prime importance for protostellar evolution, dictate the luminosities and thus play a key role in deciphering the current discrepancy between observational surveys and theoretical models. In this letter, we provide quantitative estimates on the evolution of the radiative efficiency and accretion rate obtained from self-consistent, high-resolution, radiative hydrodynamic simulations performed using the codes PLUTO and RAMSES. The main highlight of our result is that the radiative efficiency reaches unity, that is, supercriticality, relatively quickly after protostellar birth. Supercriticality at the accretion shock is a necessary condition for cold accretion. Our results thus support a rapid transition to the cold accretion scenario, which is one of the assumptions used in Pre-Main Sequence (PMS) models working towards solutions to explain observational data. We briefly discuss the implications of the time evolution of the radiative efficiency factor in the context of the luminosity problem, the Protostellar Luminosity Function (PLF), PMS evolution, accurate sink properties, and the stellar Initial Mass Function (IMF).

Active galactic nuclei (AGN) emit radiation via accretion across the entire energy spectrum. While the standard disk and corona model can somewhat describe this emission, it fails to predict specific features such as the soft X-ray excess, the short-term optical/UV variability, and the observed UV/X-ray correlation in AGN. In this context, the fraction of AGN emission in different bands (i.e., bolometric corrections) can be useful to better understand the accretion physics of AGN. Past studies have shown that the X-ray bolometric corrections are strongly dependent on the physical properties of AGN, such as their luminosities and Eddington ratios. However, since these two parameters depend on each other, it has been unclear which is the main driver of the X-ray bolometric corrections. We present here results from a large study of hard-X-ray-selected (14–195 keV) nearby (z < 0.1) AGN. Based on our systematic analysis of the simultaneous optical-to-X-ray spectral energy distributions of 236 unobscured AGN, we found that the primary parameter controlling the X-ray bolometric corrections is the Eddington ratio. Our results show that, while the X-ray bolometric correction increases with the bolometric luminosity for sources with intermediate Eddington ratios (0.01–1), this dependence vanishes for sources with lower Eddington ratios (<0.01). This could be used as evidence for a change in the accretion physics of AGN at low Eddington ratios.

Recent advances in quantum technologies have enabled quantum simulation of gauge theories -- some of the most fundamental frameworks of nature -- in regimes far from equilibrium, where classical computation is severely limited. These simulators, primarily based on neutral atoms, trapped ions, and superconducting circuits, hold the potential to address long-standing questions in nuclear, high-energy, and condensed-matter physics, and may ultimately allow first-principles studies of matter evolution in settings ranging from the early universe to high-energy collisions. Research in this rapidly growing field is also driving the convergence of concepts across disciplines and uncovering new phenomena. In this Review, we highlight recent experimental and theoretical developments, focusing on phenomena accessible in current and near-term quantum simulators, including particle production and string breaking, collision dynamics, thermalization, ergodicity breaking, and dynamical quantum phase transitions. We conclude by outlining promising directions for future research and opportunities enabled by available quantum hardware.

Dark matter decays into invisible particles can leave an imprint in large-scale structure surveys due to a characteristic redshift-dependent suppression of the power spectrum. We present a model with two quasi-degenerate singlet fermions, $χ_1$ and $χ_2$, in which the heavier state decays as $χ_2 \to \bar χ_1 νν$ on cosmological time-scales, and that also accommodates non-zero neutrino masses. Remarkably, for parameters that yield the correct dark matter abundance via freeze-in and reproduce the observed neutrino masses, dark matter decay can produce detectable signals in forthcoming large-scale structure surveys, a diffuse anti-neutrino flux accessible to JUNO, and a gamma-ray line within the energy range probed by COSI. Both the cosmological lifetime of $χ_2$ as well as the small (radiatively induced) mass splitting among $χ_{1,2}$ are a natural consequence of the mechanism of neutrino mass generation within this model. This highlights the potential role of large-scale structure surveys in probing some classes of neutrino mass models.

The de novo synthesis of life from non-living matter represents a bold scientific challenge, advancing our understanding of life's minimal requirements and offering revolutionary applications in biotechnology. We explore fuel-dependent synthetic cells based on complex coacervate droplets, which lack membranes and readily take up reactants. Given their fuel-dependent nature, these droplets emerge and grow when fuel is abundant but dissolve under starvation conditions, mimicking the non-equilibrium nature of life. However, their ability to produce offspring—a key requirement for life—has remained elusive. Moreover, their rescue in repetitive fueling-starvation experiments has not been demonstrated. Our work elucidates a mechanism of producing offspring by synthetic cells driven by solid-like speckles in droplets liberated as offspring. By fine-tuning parameters, we control offspring number and survival. Finally, refueling sustains second-generation synthetic cells. This system provides a platform for coupling offspring production with self-replicating molecules, paving the way for synthetic cells capable of Darwinian evolution.

In the standard formation models of terrestrial planets in the solar system and close-in super-Earths in non-resonant orbits recently discovered by exoplanet observations, planets are formed by giant impacts of protoplanets or planetary embryos after the dispersal of protoplanetary disk gas in the final stage. This study aims to theoretically clarify a fundamental scaling law for the orbital architecture of planetary systems formed by giant impacts. In the giant impact stage, protoplanets gravitationally scatter and collide with one another to form planets. Using {\em N}-body simulations, we investigate the orbital architecture of planetary systems formed from protoplanet systems by giant impacts. As the orbital architecture parameters, we focus on the mean orbital separation between two adjacent planets and the mean orbital eccentricity of planets in a planetary system. We find that the orbital architecture is determined by the ratio of the two-body surface escape velocity of planets $v_\mathrm{esc}$ to the Keplerian circular velocity $v_\mathrm{K}$, $k$ = The mean orbital separation and eccentricity are about $2 ka$ and $0.3 k$, respectively, where $a$ is the system semimajor axis. With this scaling, the orbital architecture parameters of planetary systems are nearly independent of their total mass and semimajor axis.

A five-dimensional minimal supergravity theory coupled to vector and hypermultiplets is specified by a set of trilinear couplings, given by an intersection form $C_{IJK}$, and gravitational couplings specified by an integer-valued vector $a_I$ and is consistent when these data define an integral cubic form. For every Calabi-Yau threefold reduction of M-theory, this condition is satisfied automatically. Via suitable redefinitions of the basis of 5D vectors, this is also shown to be the case for the circle reductions of six-dimensional anomaly-free (1,0) theories. When the 6D theory has a $\mathbb{Z}_k$ gauge symmetry, we point out that the consistency of the circle reduction with nontrivial $\mathbb{Z}_k$ holonomy is closely related to 6D constraints derived by Monnier and Moore. These constraints are extended to semidirect products with continuous gauge groups $\mathbb{Z}_k \ltimes G$ and CHL-like circle compactifications. When $\mathbb{Z}_k$ acts on anti-self-dual tensor fields of 6D supergravity, there should be a nontrivial action of holonomy on the topological Green-Schwarz terms.

Stellar shells and streams are remnants of satellite galaxies visible around galaxies. Advances in low-surface-brightness observations and increasing resolution of cosmological simulations now allow investigating the properties and origin of these features. The metallicity, age, and velocity dispersion of shells and streams are investigated to infer their progenitor galaxies properties. We employed the hydrodynamical cosmological simulations Magneticum Pathfinder to extract these properties and identify the progenitors of the shells and streams. We compared to observational results from surveys and individual galaxies, matching and testing the methodology used in observations. Mock observations of shells and streams agree well with observational data regarding their morphology and spatial distribution. We find that both types of features are associated with localized depressions in stellar velocity dispersion compared to the surrounding regions. They are not as clearly distinct in metallicity and ages, though overall shells and more metal rich and streams are younger. We confirm results from idealized models that shells form commonly from radial major mergers but also through minor mergers, while streams usually form from minor mergers on circular orbits. We do not find the widths of streams to correlate with the half-mass radii of their progenitors, but the progenitors follow the mass-metallicity relation. On average, the masses measured for shells and streams approximately corresponds to 20% of the progenitor mass. We introduce a class of star-forming streams, which originate from in-situ star formation rather than the disruption of a satellite galaxy. Measuring stellar population properties of shells and streams provides the means to reconstruct the progenitor properties, and especially distinguish those streams that are not made through the disruption of a galaxy but formed in-situ.

In this paper we argue that the information load carried by a black hole affects its classical perturbations. We refer to this phenomenon as the ``swift memory burden effect" and show that it is universal for objects of high efficiency of information storage. The effect is expected to have observable manifestations, for example, in mergers of astrophysical black holes in Einstein gravity. The black holes with different information loads, although degenerate in the ground state, respond very differently to perturbations. The strength of the imprint is controlled by the memory burden parameter which measures the fraction of the black hole's memory space occupied by the information load. This represents a new macroscopic quantum characteristics of a black hole. We develop a calculable theoretical framework and derive some master formulas which we then test on explicit models of black holes as well as on solitons of high capacity of information storage. We show that the effect must be significant for the spectroscopy of both astrophysical and primordial black holes and can be potentially probed in gravitational wave experiments. We also provide a proposal for the test of the memory burden phenomenon in a table-top laboratory setting with cold bosons.

The dispersion measure (DM) of fast radio bursts (FRBs) is sensitive to the electron distribution in the Universe, making it a promising probe of cosmology and astrophysical processes such as baryonic feedback. However, cosmological analyses of FRBs require knowledge of the contribution to the observed DM coming from the FRB host. The size and distribution of this contribution is still uncertain, thus significantly limiting current cosmological FRB analyses. In this study, we extend the baryonification (BCM) approach to derive a physically-motivated, analytic model for predicting the host contribution to FRB DMs. By focusing on the statistical properties of FRB host DMs, we find that our simple model is able to reproduce the probability distribution function (PDF) of host halo DMs measured from the CAMELS suite of hydrodynamic simulations, as well as their mass- and redshift dependence. Furthermore, we demonstrate that our model allows for self-consistent predictions of the host DM PDF and the matter power spectrum suppression due to baryonic effects, as observed in these simulations, making it promising for modelling host-DM-related systematics in FRB analyses. In general, we find that the shape of the host DM PDF is determined by the interplay between the FRB and gas distributions in halos. Our findings indicate that more compact FRB profiles require shallower gas profiles (and vice versa) in order to match the observed DM distributions in hydrodynamic simulations. Furthermore, the analytic model presented here shows that the shape of the host DM PDF is highly sensitive to the parameters of the BCM. This suggests that this observable could be used as an interesting test bed for baryonic processes, complementing other probes due to its sensitivity to feedback on galactic scales. We further discuss the main limitations of our analysis, and point out potential avenues for future work.

It has been suggested some time ago that many black hole properties are not specific to gravity, but rather are shared by a large class of objects, the so-called saturons, that saturate the quantum field theoretic upper bound on microstate degeneracy. By now, various aspects of this universality have been understood and demonstrated in a number of explicit examples. In the present paper, we add one more brick to the building by showing that the decay of a simple two-dimensional saturated soliton copies some key aspects of the black hole decay as well as of the information retrieval. In particular, we study the evaporation process of a classically-stable vacuum bubble of a spontaneously broken $SU(N)$-symmetry, coupled to massless fermions. We show that the decay rate as well as the characteristic energy of the emitted quanta are given by the inverse size of the object, in striking similarity with the Hawking evaporation of a black hole. The time-scale of information retrieval also matches the one previously suggested for a black hole by Page. We give the semiclassical derivation of the phenomenon as well as its fully quantum resolution as a decaying coherent state of Goldstone bosons. The universal nature of the effect and its microscopic understanding support the analogous quantum portrait of a black hole as a saturated coherent state of gravitons.

Context. The dust-to-gas ratio in the protoplanetary disc, which is likely imprinted into the host star metallicity, is a property that plays a crucial role during planet formation. On the observational side, statistical studies based on large exoplanet datasets have determined various correlations between planetary characteristics and host star metallicity. Aims. We aim to constrain planet formation and evolution processes by statistically analysing planetary systems produced at different metallicities by a theoretical model, and we compare them with the correlations derived from observational samples. Methods. We used the Generation III Bern model of planet formation and evolution to generate synthetic planetary systems at different metallicities. This global model incorporates the accretion of planetesimals and gas, planetary migration, N-body interactions between embryos, giant impacts, and protoplanetary disc evolution, as well as the planets' long-term contraction and atmospheric loss of gaseous envelopes. Using synthetic planets biased to observational completeness, we analysed the impact of stellar metallicity on planet occurrence rates, orbital periods, eccentricities, and the morphology of the radius valley. Results. Based on our nominal model, we find that (1) the occurrence rates of large giant planets and Neptune-sized planets are positively correlated with [Fe/H], while small sub-Earths exhibit an anti-correlation. In between, at radii of 1 to 3.5 R_⊕, the occurrence rate first increases and then decreases with increasing [Fe/H], with an inflection point at ~0.1 dex. (2) Planets with orbital periods shorter than ten days are more likely to be found around stars with a higher metallicity, and this tendency weakens with increasing planet radius. (3) Both giant planets and small planets exhibit a positive correlation between the eccentricity and [Fe/H], which could be explained by the self-excitation and perturbation of outer giant planets. (4) The radius valley deepens and becomes more prominent with increasing [Fe/H], accompanied by a lower super-Earth-to-sub-Neptune ratio. Furthermore, the average radius of the planets above the valley (2.1–6 R_⊕) increases with [Fe/H]. Conclusions. Our nominal model successfully reproduces many observed correlations with stellar metallicity either quantitatively or qualitatively, and supports the description of physical processes and parameters included in the Bern model. Quantitatively, the dependence of orbital eccentricity and period on [Fe/H] predicted by the synthetic population, however, is significantly weaker than observed. This discrepancy likely arises because the model only accounts for planetary interactions for the first 100 Myr and neglects the effects of the stellar environment (e.g. clusters, binaries). This suggests that long-term dynamical interactions between planets, along with the impact of binaries and/or companions, can drive the system towards a dynamically hotter state.

[abridged] Accurately accounting for the AGN phase in galaxy evolution requires a large, clean AGN sample. This is now possible with SRG/eROSITA. The public Data Release 1 (DR1, Jan 31, 2024) includes 930,203 sources from the Western Galactic Hemisphere. The data enable the selection of a large AGN sample and the discovery of rare sources. However, scientific return depends on accurate characterisation of the X-ray emitters, requiring high-quality multiwavelength data. This paper presents the identification and classification of optical and infrared counterparts to eRASS1 sources using Gaia DR3, CatWISE2020, and Legacy Survey DR10 (LS10) with the Bayesian NWAY algorithm and trained priors. Sources were classified as Galactic or extragalactic via a Machine Learning model combining optical/IR and X-ray properties, trained on a reference sample. For extragalactic LS10 sources, photometric redshifts were computed using Circlez. Within the LS10 footprint, all 656,614 eROSITA/DR1 sources have at least one possible optical counterpart; about 570,000 are extragalactic and likely AGN. Half are new detections compared to AllWISE, Gaia, and Quaia AGN catalogues. Gaia and CatWISE2020 counterparts are less reliable, due to the surveys shallowness and the limited amount of features available to assess the probability of being an X-ray emitter. In the Galactic Plane, where the overdensity of stellar sources also increases the chance of associations, using conservative reliability cuts, we identify approximately 18,000 Gaia and 55,000 CatWISE2020 extragalactic sources. We release three high-quality counterpart catalogues, plus the training and validation sets, as a benchmark for the field. These datasets have many applications, but in particular empower researchers to build AGN samples tailored for completeness and purity, accelerating the hunt for the Universe most energetic engines.

Characterizing exoplanets' spectra is a crucial step in understanding the chemical and physical processes shaping their atmospheres and constraining their formation and evolutionary history. The X-SHYNE library is a homogeneous sample of 43 medium-resolution (R_λ ~ 8000) infrared (0.3–2.5 μm) spectra of young (<500 Myr), low-mass (<20 M_Jup), and cold (T_eff ~600–2000 K) isolated brown dwarfs and wide-separation companions observed with the VLT/X-Shooter instrument. To characterize our targets, we performed a global comparative analysis. We first applied a semiempirical approach. By refining their age and bolometric luminosity, we derived key atmospheric and physical properties, such as T_eff, mass, surface gravity (g), and radius, using the evolutionary model COND03. These results were then compared with the results from a synthetic analysis based on three self-consistent atmospheric models: the cloudy models Exo-REM and Sonora Diamondback, and the cloudless model ATMO. To compare our spectra with these grids we used the Bayesian inference code ForMoSA. We found similar L_bol estimates between both approaches, but an underestimated T_eff from the cloudy models, likely due to a lack of absorbers that could dominate the J and H bands of early L. We also observed a discrepancy in the log(g) estimates, which are dispersed between 3.5 and 5.5 dex for mid-L objects. We interpret this as a bias caused by a range of rotational velocities leading to cloud migration toward equatorial latitudes, combined with a variety of viewing angles that result in different observed atmospheric properties (cloud column densities, atmospheric pressures, etc.). This interpretation is supported by the correlation of the color anomaly Δ(J–K) of each object with log(g) and the parameter f_sed that drives the sedimentation of the clouds. Finally, while providing robust estimates of [M/H] and C/O for individual objects remains challenging, the X-SHYNE library globally suggests solar values that are consistent with a formation via stellar formation mechanisms. This study highlights the strength of homogeneous datasets in performing comparative analyses, reducing the impact of systematics, and ensuring robust conclusions while avoiding overinterpretation.

We study the implications of non-invertible chiral symmetry in a four-dimensional U(1) gauge theory coupled to massless fermions with electromagnetic $SL(2,\mathbb{Z})$ duality. This is done by deriving the Adler-Bell-Jackiw anomaly of massless QED in the dual frame that is used to explicitly construct the symmetry defect operator as well as the conserved two-form symmetry current. As expected, the non-invertible chiral symmetry is covariant under the duality transformation. This has implications for understanding the nature of kinetic and topological terms in the dual frame and for axion electrodynamics. In particular, we show that to generate an axion potential from a dyon loop, the one-form magnetic symmetry must be explicitly broken by a mutually non-local charged state with nonzero pairwise helicity.

Mean-field dynamo theory, describing the evolution of large-scale magnetic fields, has been the mainstay of theoretical interpretation of magnetism in astrophysical objects such as the Sun for several decades. More recently, three-dimensional magnetohydrodynamic simulations have reached a level of fidelity where they capture dynamo action self-consistently on local and global scales without resorting to parametrization of unresolved scales. Recent global simulations also capture many of the observed characteristics of solar and stellar large-scale magnetic fields and cycles. Successful explanation of the results of such simulations with corresponding mean-field models is a crucial validation step for mean-field dynamo theory. Here the connections between mean-field theory and current dynamo simulations are reviewed. These connections range from the numerical computation of turbulent transport coefficients to mean-field models of simulations, and their relevance to the solar dynamo. Finally, the most notable successes and current challenges in mean-field theoretical interpretations of simulations are summarized.

Context. Fast X-ray transients (FXTs) are a rare and poorly understood phenomenon with a variety of possible progenitors. The launch of the Einstein Probe (EP) mission has facilitated a rapid increase in the real-time discovery and follow-up of FXTs. Aims. We focus on the recent EP discovered transient EP241021a, which shows a peculiar panchromatic behavior, with the aim of understanding its origin. Methods. We obtained optical and near-infrared multiband imaging and spectroscopy with the Fraunhofer Telescope at Wendelstein Observatory, the Hobby-Eberly Telescope, and the Very Large Telescope of the newly discovered EP transient EP241021a over the first 100 days of its evolution. Results. EP241021a was discovered by EP as a soft X-ray trigger, but was not detected at gamma-ray frequencies. The observed soft X-ray prompt emission spectrum is consistent with nonthermal radiation, which requires at least a mildly relativistic outflow with a bulk Lorentz factor Γ ≳ 4. The optical and near-infrared light curve displays a two-component behavior, where an initially fading component, ∼ t^‑1, transitions to a rise steeper than ∼ t³ after a few days, before peaking at an absolute magnitude of M_r ≈ ‑21.8 mag and quickly returning to the initial decay. Standard supernova models are unable to reproduce either the absolute magnitude or the rapid timescale (< 2 d) of the rebrightening. The X-ray, optical and near-infrared spectral energy distributions display a red color, r ‑ J ≈ 0.8 mag, and point to a nonthermal origin (∼ ν^‑1) for the broadband emission. Considering a gamma-ray burst as a plausible scenario, we favor a refreshed shock as the cause of the rebrightening. This is consistent with the inference of an at least mildly relativistic outflow based on the prompt trigger. Conclusions. Our results suggest a link between EP-discovered FXTs and gamma-ray bursts, despite the lack of gamma-ray detections for the majority of EP transients.

Quasi-periodic eruptions (QPEs) are recurring X-ray bursts originating from the vicinity of supermassive black holes, but their driving mechanisms remain under debate. This study analyzes new NICER observations of QPEs in Ansky (a transient event in the nucleus of the galaxy SDSS J1335+0728), taken between January and June 2025. By examining flare durations, peak-to-peak recurrence times, and profiles, we compare the 2025 data with those from 2024 to investigate changes in energy, timescales, and flare shapes. The 2025 QPEs are found to be four times more energetic, with recurrence times of approximately 10 days and flare durations ranging from 2.5 to 4 days, making them both about twice as long as in 2024. Additionally, the flare profiles have become more asymmetric, showing longer decays. We explore different theoretical scenarios to explain the observed properties of the QPEs in Ansky, including evolving stream-disk interactions in an extreme mass-ratio inspiral (EMRI) system as a potential mechanism behind the observed changes in recurrence time and energetics, while also considering alternative models based on mass transfer and accretion disk instabilities. Continued observational efforts will be crucial for unveiling the nature of Ansky.

The origin of the high-<inline-formula><tex-math>$\alpha$</tex-math></inline-formula> component of the Galactic bulge remains debated, unlike the bar-driven origin of the low-<inline-formula><tex-math>$\alpha$</tex-math></inline-formula> bulge. We examine the metallicity-dependent dynamical properties of high-[Mg/Fe] stars in the bar region, using samples of low- and high-[Mg/Fe] stars from APOGEE DR17, complemented by the PIGS catalogue of <inline-formula><tex-math>${\rm [Fe/H]}< -1$</tex-math></inline-formula> stars. The mean Galactocentric rotational velocity <inline-formula><tex-math>$\overline{V}_{\phi }(R)$</tex-math></inline-formula> is nearly cylindrical for both low- and high-[Mg/Fe] stars across the bulge and outer bar. <inline-formula><tex-math>$\overline{V}_{\phi }(R)$</tex-math></inline-formula> of high-[Mg/Fe] stars with <inline-formula><tex-math>${\rm [Fe/H]}\ge -0.6$</tex-math></inline-formula> is similar within errors to low-[Mg/Fe] stars in the bulge, and 10–20 per cent lower in the outer bar. The mean radial velocity field of these stars exhibits a quadrupole pattern similar to low-[Mg/Fe] stars. Integrating orbits in realistic barred Galactic potentials, these model-independent properties correspond to a peanut bulge in the orbital density distributions for high-[Mg/Fe] stars with <inline-formula><tex-math>${\rm [Fe/H]}\ge -0.6$</tex-math></inline-formula>, transitioning toward a more spheroidal structure at lower metallicities. Additionally, <inline-formula><tex-math>$\overline{V}_{\phi }({\rm [Fe/H]})$</tex-math></inline-formula> for stars increases steeply as metallicity increases from about [Fe/H] <inline-formula><tex-math>$\sim -1.3$</tex-math></inline-formula>, resembling the spin-up observed at larger Galactic radii. This is accompanied by a transition in the dominant orbit families, from co- and counter-rotating <inline-formula><tex-math>${\rm cloud\, A}$</tex-math></inline-formula> and <inline-formula><tex-math>${\rm x_4}$</tex-math></inline-formula> orbits at low metallicities to co-rotating bar-supporting <inline-formula><tex-math>${\rm x_1}$</tex-math></inline-formula> family tree, <inline-formula><tex-math>${\rm box}$</tex-math></inline-formula>, and <inline-formula><tex-math>${\rm cloud\, A}$</tex-math></inline-formula> orbits at solar metallicity. Our results strengthen the case that the bulk of the high-[Mg/Fe] component in the bar region evolved from an <inline-formula><tex-math>$\alpha$</tex-math></inline-formula>-enhanced disc, while metal-poor stars with <inline-formula><tex-math>${\rm [Fe/H]}< -1$</tex-math></inline-formula> trace a more turbulent origin.

We propose a new model to explain the KM3NeT neutrino event through a low reheating scenario with a suppression in the GW spectrum originating from cosmic string networks. To achieve this, we extend the SM gauge sector by an abelian gauge symmetry and a singlet scalar. Once the abelian gauge symmetry spontaneously breaks, the extra gauge boson acquires mass and becomes a suitable Dark Matter (DM) candidate. Due to the kinetic mixing with the hypercharge gauge group, DM can decay into SM particles. To explain the KM3NeT signal, we need $\mathcal{O}(100)$ PeV DM, which can be produced in the correct order of DM density in a low reheating scenario. In this scenario, the overabundance issue of heavy DM can be tackled by diluting its abundance through the continuous injection of entropy when the matter-like inflaton decays into the SM bath. Using the low reheating scenario, we can obtain the correct value of DM density both for freeze-out and freeze-in mechanisms for super-heavy DM. Moreover, we have studied the Gravitational Waves (GWs) produced from cosmic strings, which fall within the detectable range of future proposed GW experiments. Additionally, the dominance of a quadratic inflaton potential before the reheating temperature changes the temperature-scale factor relation, which suppresses the GW spectrum at higher frequencies. Choosing an arbitrarily low reheating temperature provides only a tiny fraction of the DM density due to dilution from entropy injection. This fraction of the vector DM suggests that only the extragalactic contribution is relevant in the KM3NeT event because DM lifetime is shorter than the age of the Universe.

We study the physics potential of heavy QCD axions at high-energy muon colliders. Unlike typical axion-like particles, heavy QCD axions solve the strong CP problem with phenomenology driven by the anomalous gluon ($aG\widetilde G$) couplings. Several ultraviolet scenarios are presented in which QCD axions with TeV-scale masses and decay constants arise consistently with a solution to both the strong CP problem and the axion quality problem. We perform a detailed collider analysis for both a 3 and 10~TeV muon collider, focusing on hadronic axion decays that gives rise to a dijet-resonance signature. Our projections for the axion discovery reach in the multi-TeV mass range demonstrate that a muon collider can significantly extend sensitivity to heavy QCD axions compared to existing experiments.

A first-order phase transition could occur in the late universe when vacuum energy begins dominating the energy density ($z \lesssim 0.3$) and convert some latent heat into other forms such as invisible radiation. This generic possibility also has concrete motivation in particle physics models which invoke a multitude of vacua to address theoretical puzzles. The naïve constraint on such an event comes from measurements of the Hubble expansion rate, but this can only probe transitions involving $\mathcal{O}(10)\%$ of the dark energy. In this work, we show that significantly tighter constraints appear when accounting for phase transition fluctuations affecting CMB photon propagation anisotropically, akin to the integrated Sachs-Wolfe effect. For instance, if a completed phase transition has $β/H_\star\lesssim 25$, current CMB data limits the associated vacuum energy released to less than $1\%$ of the dark energy. A transition to negative vacuum energy (quasi-anti-de Sitter) is allowed only for $β/H_\star \gtrsim 300$. For $β/H_\star \lesssim 500$, the universe will not crunch for at least $14$ Gyr.

A crucial ingredient affecting fast neutrino flavor conversion in core-collapse supernovae (SNe) is the shape of the angular distribution of the electron-neutrino lepton number (ELN). The presence of an ELN crossing signals favorable conditions for flavor conversion. However, the dependence of ELN crossings on the SN properties is only partially understood. We investigate a suite of 12 spherically symmetric neutrino-hydrodynamics simulations of the core collapse of a SN with a mass of <inline-formula><mml:math><mml:mrow><mml:mn>18.6</mml:mn><mml:msub><mml:mrow><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mo>⊙</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>; each model employs different microphysics (i.e., three different nuclear equations of state, with and without muon creation) and includes or not a mixing-length treatment for proto-neutron star convection. We solve the Boltzmann equations to compute the neutrino angular distributions relying on static fluid properties extracted from each of the SN simulations in our suite for six selected postbounce times. We explore the dependence of the ELN distributions on the SN microphysics and proto-neutron star convection. We find that the latter shifts the proto-neutron star radius outward, favoring the appearance of ELN crossings at larger radii. On the other hand, muon creation causes proto-neutron star contraction, facilitating the occurrence of ELN crossings at smaller radii. These effects mildly depend on the nuclear equation of state. Our findings highlight the subtle impact of the SN microphysics, proto-neutron star convection, and neutrino transport on the ELN angular distributions.

We present a view of the stellar halo in the inner-central regions of the Milky Way (R ≲ 10 kpc) mapped by RR Lyrae stars. The combined BRAVA-RR/APOGEE RR Lyrae catalog is used to obtain a sample of 281 RR Lyrae stars located in the bulge region of the Galaxy, but with orbits indicating they belong to the inner-central halo. The RR Lyrae stars in the halo are more metal-poor than the bulge RR Lyrae stars and have pulsation properties more consistent with an accreted population. We use the Milky Way-like zoom-in cosmological simulation Auriga to compare the properties of the RR Lyrae stars to those expected from the "Gaia-Enceladus-Sausage" (GES) merger. The integrals of motions and eccentricities of the RR Lyrae stars are consistent with a small fraction of 6–9% ± 2% of the inner-central halo RR Lyrae population having originated from GES. This fraction, lower than what is seen in the solar neighborhood, is consistent with trends seen in the Auriga simulation, where a GES-like merger would have a decreasing fraction of GES stars at small Galactocentric radii compared to other accreted populations. Very few of the Auriga inner Galaxy GES-18 particles have properties consistent with belonging to a bulge population with (z_max < 1.1 kpc), indicating that no (or very few) RR Lyrae stars with bulge orbits should have originated from GES.

Motivated by JWST observations of dense, clumpy, and clustered high-redshift star formation, we simulate the hierarchical assembly of nine <inline-formula><tex-math>$M_{\mathrm{cl}}={10^6}\:\mathrm{M_\odot }$</tex-math></inline-formula> star clusters using the BIFROST N-body code. Our low-metallicity models (<inline-formula><tex-math>$0.01Z_\odot$</tex-math></inline-formula>) 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 <inline-formula><tex-math>$M_\mathrm{\bullet }\sim {6200}\:\mathrm{M_\odot }$</tex-math></inline-formula>, stellar multiplicity boosting the IMBH masses by a factor of 2–3. The IMBHs tidally disrupt (TDE) <inline-formula><tex-math>$\sim 50$</tex-math></inline-formula> stars in 10 Myr with peak TDE rates up to <inline-formula><tex-math>$\Gamma \sim 5\times 10^{-5}$</tex-math></inline-formula> <inline-formula><tex-math>$\rm{yr}^{-1}$</tex-math></inline-formula> per cluster. These IMBHs are natural seeds for supermassive black holes (SMBHs) and the hierarchical assembly frequently leads to <inline-formula><tex-math>$>2$</tex-math></inline-formula> SMBH seeds per cluster and their rapid mergers (<inline-formula><tex-math>$t< 10$</tex-math></inline-formula> Myr). We propose that a gravitational wave (GW)-driven merger of IMBHs with <inline-formula><tex-math>${1000}\:\mathrm{M_\odot } \lesssim M_\bullet \lesssim {10\,000}\:\mathrm{M_\odot }$</tex-math></inline-formula> with comparable masses is a characteristic GW fingerprint of SMBH seed formation at redshifts <inline-formula><tex-math>$z>10$</tex-math></inline-formula>, 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.

Supernova (SN) 1987A is a celebrated laboratory in searches for gamma-ray flashes produced by the radiative decay of sub-GeV particles such as axion-like particles (ALPs), sterile neutrinos, and novel gauge bosons. At large couplings, however, particles decay rapidly inside the stellar envelope, which results in a suppression of the signal. Focusing on the prototypical example of ALPs with a photon coupling, we show that core-collapse SNe of Type Ic are much less affected by this attenuation, thanks to the compactness of their progenitors ensuing from the loss of their envelope. While Fermi-LAT may miss the brief gamma-ray flash from a single Type Ic SN, their high rate allows for a statistical approach: by stacking many events, we can obtain constraints that significantly surpass those from SN 1987A at large couplings. Our approach can be extended to any feebly interacting particle featuring a decay channel into photons.

We explore the possibility that the underlying flavour structure of the Standard Model could be determined by mass chains on a fractal geometry. We consider, as an example, the theory space on a Sierpinski-like geometry. The fermion mass chains on a Sierpinski-like geometry with three decorations (iterations) lead to three zero modes, which can be identified with the three generations of the Standard Model. This framework also reproduces the measured charged and neutral lepton masses and mixing angles with very few parameters. We also briefly discuss the possible extension to the quark sector.

The age of the Local Bubble (LB) can be used to constrain the timescales, on which the interstellar medium in the solar neighborhood is evolving. Previous estimates have put the age of the LB at $\gtrsim 14\,\text{Myr}$, and suggested that its expansion was powered by $\sim 15-20$ SNe, yet in a companion paper we have seen hints that this age might be too high. Following up on these hints, we aim to place new constraints on the age of the LB. We reconstruct the geometry and momentum of the LB using publicly available 3D dust maps to compare its geometry to that of the high-quality sample of simulated supernova remnants in the SISSI project. We find that, in contrast to previous estimates, $\gtrsim 20$ SNe over $\sim 4\,\text{Myr}$ are required to explain both the momentum and the size of the LB. The julia source-code for our analysis is made available at doi.org/10.5281/zenodo.17054923. Previous estimates of the age of the LB have seemingly overestimated its age and underestimated the number of SNe powering its expansion. Our results are in tension with the assumption that the LB is powered solely by SNe associated with the nearby Scorpius-Centraurus OB association, which appears to have stopped forming stars at about the same time as the LB began to expand. In light of this new evidence, our results cast serious doubts on the claim that star formation in the solar neighborhood was driven by the expansion of the LB, and might have instead quenched it.

Gravitational-wave (GW) neutron star mergers with an associated electromagnetic counterpart constitute powerful probes of binary evolution, the production sites of heavy elements, general relativity, and the expansion of the Universe. Only a handful of candidate GW binary mergers during the fourth LIGO/Virgo/KAGRA observing run (O4) so far are believed to include a neutron star. We present optical–near-infrared follow-up observations of the candidate neutron star–black hole GW merger S250206dm. This is the first high-significance mass-gap neutron star–black hole candidate observed by multiple GW detectors (thus having a significantly smaller sky localization than one-detector events), offering the first opportunity to effectively follow up a GW event of this kind. Our GW MultiMessenger Astronomy DECam Survey (GW-MMADS) campaign consisted of a wide-field search using the Dark Energy Camera (DECam) and T80-South (T80S), as well as galaxy-targeted observations using the Southern Astrophysical Research (SOAR) imager and the Fraunhofer Telescope at Wendelstein Observatory. No viable kilonova counterpart was found in our observations. We use our observation depths to place competitive constraints on kilonova models similar to or brighter than the GW170817 kilonova AT 2017gfo within our observed fields, ruling out 100% of such models with SOAR galaxy-targeted observations and ∼43% (48%) with DECam (DECam and T80S).

We present a global neutrino oscillation analysis of models with a single large extra dimension in which right-handed neutrinos possess bulk Dirac masses. Two scenarios are considered: large extra dimensions with bulk masses and the dark dimension framework, both predicting a tower of sterile Kaluza-Klein states that mix with active neutrinos. Using data from MINOS/MINOS+, KamLAND, and Daya Bay, we perform a joint likelihood analysis. No signatures of these theories were found. Therefore, we constrain the compactification radius under different bulk mass and Yukawa coupling assumptions. Large positive bulk masses or sizable Yukawas lead to strong bounds, while small couplings or negative bulk masses remain less constrained.

Context. The cosmic time evolution of the radial structure is one of the key topics in the investigation of disc galaxies. In the build-up of galactic discs, gas infall is an important ingredient and it produces radial gas inflows as a physical consequence of angular momentum conservation since the infalling gas onto the disc at a specific radius has lower angular momentum than the circular motions of the gas at the point of impact. NGC 300 is a well-studied isolated, bulgeless, and low-mass disc galaxy ideally suited for an investigation of galaxy evolution with radial gas inflows. Aims. Our aim is to investigate the effects of radial gas inflows on the physical properties of NGC 300, for example the radial profiles of HI gas mass and star formation rate (SFR) surface densities, specific star formation rate (sSFR), and metallicity, and to study how the metallicity gradient evolves with cosmic time. Methods. A chemical evolution model for NGC 300 was constructed by assuming its disc builds up progressively by the infalling of metal-poor gas and the outflowing of metal-enriched gas. Radial gas inflows were also considered in the model. We used the model to build a bridge between the available data (e.g. gas content, SFR, and chemical abundances) observed today and the galactic key physical processes. Results. Our model including the radial gas inflows and an inside-out disc formation scenario can simultaneously reproduce the present-day observed radial profiles of HI gas mass surface density, SFR surface density, sSFR, gas-phase, and stellar metallicity. We find that, although the value of radial gas inflow velocity is as low as ‑0.1 km s^‑1, the radial gas inflows steepen the present-day radial profiles of HI gas mass surface density, SFR surface density, and metallicity, but flatten the radial sSFR profile. Incorporating radial gas inflows significantly improves the agreement between our model predicted present-day sSFR profile and the observations of NGC 300. Our model predictions are also in good agreement with the star-forming galaxy main sequence and the mass-metallicity relation of star-forming galaxies. It predicts a significant flattening of the metallicity gradient with cosmic time. We also find that the model predicted star formation has been more active recently, indicating that the radial gas inflows may help to sustain star formation in local spirals, at least in NGC 300.

Galaxies with high star-formation rate surface densities often host large-scale outflows that redistribute energy, momentum, and baryons between the interstellar medium and the halo, making them a key feedback channel regulating galaxy evolution. Despite their importance, the driving physics behind galactic outflows and their interaction with the surrounding halo is yet to be fully understood. In particular, the influence of a pre-existing reservoir of cosmic rays (CRs) in galaxy halos has not been clearly established. We determine the conditions required to launch outflows in the presence of halo CRs and investigate how CR pressure gradients modify outflow speeds. We find that CR halos suppress the development of large-scale, CR-driven winds and redirect CR feedback toward local recycling flows. Slow outflows are therefore more likely in young galaxies lacking extended CR halos, while fast winds in intense starbursts are dominated by momentum injection and largely unaffected by halo CRs.

We study a family of generalized hypergeometric integrals defined on punctured Riemann surfaces of genus g. These integrals are closely related to g-loop string amplitudes in chiral splitting, where one leaves the loop-momenta, moduli and all but one puncture un-integrated. We study the twisted homology groups associated to these integrals, and determine their intersection numbers. We make use of these homology intersection numbers to write a double-copy formula for the "complex" version of these integrals -- their closed-string analogues. To verify our findings, we develop numerical tools for the evaluation of the integrals in this work. This includes the recently introduced Enriquez kernels -- integration kernels for higher-genus polylogarithms.

Fast neutrino flavor conversion may impact the explosion mechanism and nucleosynthesis in core-collapse supernovae. A necessary condition for fast flavor conversion is the presence of crossings in the angular distribution of the electron-neutrino lepton number (ELN) crossing. Because of the computational costs, flavor-dependent angular distributions are not computed by the vast majority of state-of-the-art hydrodynamical simulations; instead, angular distributions are reconstructed employing approximate methods in postprocessing. In this work, we evaluate the performance of four methods adopted to diagnose the existence of ELN crossings. For selected postbounce times, we extract the fluid and thermodynamic properties from spherically symmetric supernova simulations for an <inline-formula><mml:math><mml:mrow><mml:mn>18.6</mml:mn><mml:msub><mml:mrow><mml:mi>M</mml:mi></mml:mrow><mml:mrow><mml:mo>⊙</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> progenitor, testing cases with and without muons as well as with and without mixing-length treatment of protoneutron star convection. We compare the occurrence of crossings in the angular distributions obtained by solving the Boltzmann equations with those in distributions reconstructed from angular moments of our Boltzmann solutions by using the maximum entropy and Minerbo schemes, and also with crossings identified via a polynomial weighting function applied to the angular moments. Our results show that the polynomial method and the Minerbo closure scheme have severe limitations. The maximum entropy approach captures most of the forward crossings, although it fails to reproduce or misidentifies crossings in a subset of our models. These findings highlight the need for robust modeling of the neutrino angular properties in order to assess the impact of flavor conversion on the supernova mechanism.

In this paper, we present analytical results for the two-loop QCD corrections to the production of two partons or a photon and a parton in hadronic collisions, mediated by loops of massive quarks. These amplitudes involve Feynman integrals defined on an elliptic curve. We compute them by generalizing our recent results for the production of two photons to include additional crossings of the corresponding master integrals, which we compute in terms of the same basis of independent iterated integrals. We discuss the analytical properties of the amplitudes, highlighting the cancellations of a large number of elliptic differential forms in their finite remainders. Finally, we elaborate on a strategy for their numerical evaluation based on generalized series expansions at singular points of the physical amplitude, through the introduction of suitable sets of variables that allow us to resolve all singularities.

In April 2013, the TeV blazar Markarian~421 underwent one of its most powerful emission outbursts to date. An extensive multi-instrument campaign featuring MAGIC, VERITAS, and \textit{NuSTAR} provided comprehensive very-high-energy (VHE; $E > 100$\,GeV) and X-ray coverage over nine consecutive days. In this work, we perform a detailed spectral analysis of the X-ray and VHE emissions on sub-hour timescales throughout the flare. We identify several clockwise spectral hysteresis loops in the X-rays, revealing a spectral evolution more complex than a simple harder-when-brighter trend. The VHE spectrum extends beyond 10\,TeV, and its temporal evolution closely mirrors the behavior in the X-rays. We report the first evidence of VHE spectral hysteresis occurring simultaneously with the X-ray loops. To interpret these findings, we apply a time-dependent leptonic model to 240 broadband spectral energy distributions (SEDs) binned on a 15-minute scale, allowing us to self-consistently track the particle distribution's history. Our modeling shows that the majority of the sub-hour flux and spectral variations are driven by changes in the luminosity and slope of the injected electron distribution. The required variations in the electron slope are difficult to reconcile with magnetic reconnection but are consistent with a shock-acceleration scenario where the shock compression ratio evolves by a factor of $\sim2$. The model also points to a relatively stable magnetic field and emitting region size, favoring a scenario where the emission originates from a stationary feature in the jet, such as a recollimation shock. However, this scenario requires a jet Lorentz factor that significantly exceeds values from VLBI measurements to account for the high minimum electron energy implied by the lack of variability in the optical band.

This is a pedagogical review of some recent progress in rigorously proving chiral symmetry breaking in a class of QCD-like theories that closely resemble the real-world QCD, namely the $SU(N_c)$ Yang-Mills theory coupled to $N_f$ flavors of massless quarks in the fundamental representation. Based on 't Hooft anomaly matching and persistent mass conditions, a general no-go theorem is formulated: assuming that the theory flows in the infrared to a fully color-screened, infrared-free phase described by color-singlet hadrons, symmetry and anomaly constraints necessarily imply spontaneous chiral symmetry breaking; conversely, any phase with unbroken chiral symmetry must retain unscreened color charges, thereby ruling out a fully color-singlet hadron description in the infrared. While these results have been widely assumed, the recent developments reviewed here establish them with a new level of rigor. The persistent mass condition, carefully formulated here, plays a central role, just as it does in the Vafa-Witten theorem on unbroken vectorlike symmetries.

MACS J0138-2155 is the only known cluster to strongly lens two supernovae (SNe), Requiem and Encore, from the same host galaxy at z=1.949. We present seven independent mass models of the galaxy cluster built using six software packages. By conducting a blind analysis (no exchanges of results between modeling teams), we quantified uncertainties due to modeling and software. Through HST, JWST and MUSE observations, we assembled high-quality data products, including eight "gold" lensed image systems consisting of 23 images with secure spectroscopic redshifts, and one "silver" system with a likely redshift value. Restricting to the gold images, we obtain overall consistent model predictions of the positions, magnifications and time delays of SN Encore and SN Requiem images, especially for models with $χ^2 \leq 25$. We predict the appearance of the next images of SNe Encore and Requiem with a time delay of >~3000 days and of ~3700 to 4000 days, respectively, based on a fiducial cosmological model of $H_0 = 70 {\rm\ km\ s^{-1}\ Mpc^{-1}}$ and $Ω_{\rm m} = 0.3$. We obtain relations between $H_0$ and the time delays of SNe Encore and Requiem. In particular, for $H_0 = 73 {\rm\ km\ s^{-1}\ Mpc^{-1}}$, the four lowest $χ^2$ models predict SN Requiem to reappear in ~Apr-Dec 2026; for $H_0 = 67 {\rm\ km\ s^{-1}\ Mpc^{-1}}$, in ~Mar-Nov 2027. Using the newly measured time delay between the two detected images of SN Encore by Pierel et al. (submitted) and our mass models, we jointly infer $H_0 = {\rm 66.9^{+11.2}_{-8.1}\ km\ s^{-1}\ Mpc^{-1}}$, where the uncertainty is dominated by that of the time delay. The long delays of the next-appearing SN Requiem and SN Encore images provide excellent opportunities to measure $H_0$ with an uncertainty of 2-3%. Our mass models form the basis for cosmological inference from this unique lens cluster with two strongly lensed SNe. (Abridged)

Extended <inline-formula><tex-math>$\mathrm{Ly\,\,\alpha }$</tex-math></inline-formula> emission is commonly observed around star-forming galaxies, opening a window for probing the neutral hydrogen gas in the circumgalactic medium (CGM). In this paper, we develop a prescription of spherically symmetric CGM gas properties and build emulators to model circularly averaged surface brightness (SB) profiles of the extended <inline-formula><tex-math>$\mathrm{Ly\,\alpha }$</tex-math></inline-formula> emission. With CGM gas properties parametrized by the density, velocity and temperature profiles, a self-shielding calculation is carried out to obtain the neutral gas distribution with ionizing photons from the ultraviolet (UV) background and star formation in the galaxy. Our calculation reveals three types of systems with distinct neutral gas distribution: non-shielded systems with the CGM being highly ionized across all radii, shielded systems with a neutral gas shell shielding the UV background, and transitional systems in between. <inline-formula><tex-math>$\mathrm{Ly\,\alpha }$</tex-math></inline-formula> SB profiles are obtained through <inline-formula><tex-math>$\mathrm{Ly\,\alpha }$</tex-math></inline-formula> radiative transfer (RT) simulations, performed for the CGM models with three kinds of <inline-formula><tex-math>$\mathrm{Ly\,\alpha }$</tex-math></inline-formula> sources: the star formation from central and satellite galaxies, and the recombination in the CGM. We build emulators to efficiently predict <inline-formula><tex-math>$\mathrm{Ly\,\alpha }$</tex-math></inline-formula> SB profiles for given model parameters and <inline-formula><tex-math>$\mathrm{Ly\,\alpha }$</tex-math></inline-formula> sources, based on Gaussian process regression. After being trained with only 180 RT simulations for each <inline-formula><tex-math>$\mathrm{Ly\,\alpha }$</tex-math></inline-formula> source, the emulators reach an overall accuracy at the level of <inline-formula><tex-math>$\sim 20$</tex-math></inline-formula> per cent. By applying the emulators to fit mock <inline-formula><tex-math>$\mathrm{Ly\,\alpha }$</tex-math></inline-formula> SB profiles constructed from our model, we find a reasonable recovery of model parameters, indicating the potential of extracting physical information of the CGM and galaxies from the observed extended <inline-formula><tex-math>$\mathrm{Ly\,\alpha }$</tex-math></inline-formula> emission.

We present the first gamma-ray burst (GRB) host galaxy with a measured absorption line and electron temperature (T<inline-formula><tex-math>$_e$</tex-math></inline-formula>) based metallicity, using the temperature sensitive [O III]<inline-formula><tex-math>$\lambda$</tex-math></inline-formula>4363 auroral line detected in the JWST/NIRSpec spectrum of the host of GRB 050505 at redshift <inline-formula><tex-math>$z=4.28$</tex-math></inline-formula>. We find that the metallicity of the cold interstellar gas, derived from the absorption lines in the GRB afterglow, of 12 + log(O/H) <inline-formula><tex-math>$\sim 7.7$</tex-math></inline-formula> is in reasonable agreement with the temperature-based emission line metallicity in the warm gas of the GRB host galaxy, which has values of 12 + log(O/H) = 7.80<inline-formula><tex-math>$\pm$</tex-math></inline-formula>0.19 and 7.96<inline-formula><tex-math>$\pm$</tex-math></inline-formula>0.21 for two common indicators. When using strong emission line diagnostics appropriate for high-z galaxies and sensitive to ionization parameter, we find good agreement between the strong emission line metallicity and the other two methods. Our results imply that, for the host of GRB050505, mixing between the warm and the cold interstellar medium along the line of sight to the GRB is efficient, and that GRB afterglow absorption lines can be a reliable tracer of the metallicity of the galaxy. If confirmed with a large sample, this suggest that metallicities determined via GRB afterglow spectroscopy can be used to trace cosmic chemical evolution to the earliest cosmic epochs and in galaxies far too faint for emission line spectroscopy, even for JWST.

Blazars are promising high-energy neutrino source candidates. Yet, leptohadronic models face challenges in describing neutrino emission within a viable energy budget, and their predictive power is limited by the commonly used single-zone approximation and the reliance on phenomenological parameters. In this work, we present a new leptohadronic model where a sub-Eddington jet evolves from magnetically- to kinetically dominated. A small fraction of the electrons and protons picked up by the jet are continuously accelerated to a power-law spectrum, estimated based on the local magnetic field strength, turbulence, and ambient density, for which we assume power-law profiles. The model parameters are thus directly tied to the jet physics and are comparable in number to typical single-zone models. We then numerically calculate the emission along the jet. Applying the model to the IceCube candidate TXS 0506+056, we find that protons are accelerated to EeV energies in the inner jet, producing a neutrino flux up to order 100 PeV that is consistent with the public 10 year IceCube point-source data. Proton emission at 0.1 pc describes the X-ray and gamma-ray data, while electron emission at the parsec scale describes the optical data. Protons carry a power of about 1% of the Eddington luminosity, showing that the model is energetically viable. The particle spectra follow $E^{-1.8}$, with diffusion scaling as $E^{0.3}$, ruling out Bohm-like diffusion. Additional particle injection near the broad line region can reproduce the 2017 flare associated to a high-energy neutrino. We also apply the model to blazar PKS 0605-085, which may be associated with a recent neutrino detected by KM3NeT above 100 PeV. The results suggest that blazars are efficient neutrino emitters at ultra-high energies, making them prime candidates for future experiments targeting this challenging energy range.

Context. The distribution of chemical elements in star-forming regions can store information on the chemical enrichment history of galaxies and particularly of recent events. Negative metallicity gradients are expected in galaxies forming inside-out. Azimuthal-averaged profiles are usually fit to the projected chemical distributions to quantify them. However, observations show that the metallicity profiles can be broken. Aims. We aim to study the diversity of metallicity profiles that can arise in the current cosmological context and compare them with available observations. Additionally, we seek to identify the physical processes responsible for breaks in metallicity profiles by using two galaxies as case studies. Methods. We analyzed central galaxies from the cosmological simulations of CIELO project, with stellar masses within the range of 10^8.5 to 10^10.5 M_⊙ at z = 0. A new algorithm, DB-A, was developed to fit multiple power laws to the metallicity profiles, enabling a flexible assessment of metallicity gradients in various galactic regions. The simulations include detailed modeling of gas components, metal-dependent cooling, star formation, and supernova feedback. Results. At z = 0, we find a diversity of shapes, with inner and outer drops and rises, and there are a few galaxies with double breaks. Inner, outer, and middle gradients are in agreement with observations. We also find that using a single linear regression to fit gradients usually traces the middle gradient well. A detailed temporal analysis of the main galaxies of a Local Group analog revealed the occurrence of inner and outer breaks at all cosmic times, with the latter being the most common feature during the evolution of our case studies. Significant variability in the metallicity gradients was found at high redshift, transitioning to more gradual evolution at lower redshifts. Most of the inner breaks have enhanced oxygen abundances in the center, which are linked to gas accretion followed by efficient star formation. Inner breaks with diluted oxygen abundances are less common and are found in galaxies with disrupted gas distributions which are affected by feedback-driven ejection of enriched gas. Outer breaks with high abundances are linked to processes such as the re-accretion of enriched material, extended star formation, and enhanced gas mixing from the circumgalactic medium. Outer breaks with diluted metallicities in the outskirts are found mainly at high redshift and are associated with the accretion of metal-poor gas from cold flows. We also highlight and illustrate the complex interplay of these processes which act often together.

In order to compress and more easily interpret Lyman-$α$ forest (Ly$α$F) datasets, summary statistics, e.g. the power spectrum, are commonly used. However, such summaries unavoidably lose some information, weakening the constraining power on parameters of interest. Recently, machine learning (ML)-based summary approaches have been proposed as an alternative to human-defined statistical measures. This raises a question: can ML-based summaries contain the full information captured by traditional statistics, and vice versa? In this study, we apply three human-defined techniques and one ML-based approach to summarize mock Ly$α$F data from hydrodynamical simulations and infer two thermal parameters of the intergalactic medium, assuming a power-law temperature-density relation. We introduce a metric for measuring the improvement in the figure of merit when combining two summaries. Consequently, we demonstrate that the ML-based summary approach not only contains almost all of the information from the human-defined statistics, but also that it provides significantly stronger constraints by a ratio of better than 1:3 in terms of the posterior volume on the temperature-density relation parameters.

Aims. We investigated the gas reservoirs, star formation properties, and environment of the ultra-diffuse galaxy GAMA 526784 to understand its formation history, the efficiency of molecular gas conversion into stars, and the possible role of an interacting companion in shaping its current morphology. Methods. We analysed low- and high-resolution CO observations to place constraints on the molecular gas content of the galaxy, compared them with HI data, and examined the star formation efficiency of GAMA 526784. The potential influence of a newly identified nearby dwarf galaxy was assessed using photometric and spatial information. Results. GAMA 526784 exhibits a regular HI reservoir (M_HI/M_⋆ = 2.88), but we are only able to place upper limits on its molecular gas mass (M_H2(5σ)/M_⋆ < 0.23). The galaxy's HI reservoir and CO non-detection can be explained by several mechanisms: (1) the predominance of CO-dark H₂, which remains invisible to CO observations but contributes to star formation; (2) a time delay in HI-to-H₂ conversion following a recent interaction; or (3) elevated turbulence inhibiting gas collapse. An identified companion, optically found at a projected distance of ∼48 kpc, shows similar colours and lies in the direction of the young star clusters in GAMA 526784, indicating a possible association. We hypothesise that this companion may have triggered the formation of the star clusters in GAMA 526784 through a high-velocity encounter. Conclusions. Our findings suggest that GAMA 526784 may have undergone a dwarf–dwarf interaction that significantly influenced its gas reservoirs and star formation activity. The presence of a nearby companion galaxy is consistent with predictions of a high-speed encounter, potentially offering a rare observational example of such an interaction in progress. We hypothesise that this encounter may have played a key role in shaping the system's recent evolution. Future observations, particularly targeting molecular gas tracers beyond CO and resolved HI observations, will be crucial in determining the true extent of GAMA 526784's cold gas reservoir and the nature of its recent star formation activity.

We investigate the high-ionization, narrow [Ne V] λ3427 line emission in a sample of over 340 ultrahard X-ray (14–195 keV) selected active galactic nuclei (AGN) drawn from the BAT AGN Spectroscopic Survey project. The analysis includes measurements in individual and stacked spectra and considers several key AGN properties such as X-ray luminosity, supermassive black hole (SMBH) mass, Eddington ratios, and line-of-sight column density. The [Ne V] λ3427 line is robustly detected in ≈43% (146/341) of the AGN in our sample, with no significant trends between the detection rate and key AGN/SMBH properties. In particular, the detection rate remains high even at the highest levels of obscuration (>70% for <inline-formula> <mml:math><mml:mi>log</mml:mi><mml:mfenced><mml:mrow><mml:msub><mml:mi>N</mml:mi><mml:mi>H</mml:mi></mml:msub><mml:mo>/</mml:mo><mml:msup><mml:mi>cm</mml:mi><mml:mrow><mml:mo>‑</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:mfenced><mml:mo>)</mml:mo><mml:mo>≳</mml:mo><mml:mn>23</mml:mn></mml:math> </inline-formula>). On the other hand, even some of our highest signal-to-noise spectra (S/N > 50) lack a robust [Ne v] detection. The typical (median) scaling ratios between [Ne v] line emission and (ultra)hard X-ray emission in our sample are <inline-formula> <mml:math><mml:mi>log</mml:mi><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mi>Ne</mml:mi><mml:mspace></mml:mspace><mml:mi>V</mml:mi><mml:mo>]</mml:mo></mml:mrow></mml:msub><mml:mo>/</mml:mo><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mn>14</mml:mn><mml:mo>‑</mml:mo><mml:mn>150</mml:mn><mml:mspace></mml:mspace><mml:mi>keV</mml:mi></mml:mrow></mml:msub><mml:mo>≃</mml:mo><mml:mo>‑</mml:mo><mml:mn>3.75</mml:mn></mml:math> </inline-formula> and <inline-formula> <mml:math><mml:mi>log</mml:mi><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mi>Ne</mml:mi><mml:mspace></mml:mspace><mml:mi>V</mml:mi><mml:mo>]</mml:mo></mml:mrow></mml:msub><mml:mo>/</mml:mo><mml:msub><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mo>‑</mml:mo><mml:mn>10</mml:mn><mml:mspace></mml:mspace><mml:mi>keV</mml:mi></mml:mrow></mml:msub><mml:mo>≃</mml:mo><mml:mo>‑</mml:mo><mml:mn>3.36</mml:mn></mml:math> </inline-formula>. The scatter on these scaling ratios, ≲0.5 dex, is comparable to, and indeed smaller than, what is found for other commonly used tracers of AGN radiative outputs (e.g., [O III] λ5007). Otherwise, we find no significant relations between the (relative) strength of [Ne v] and the basic AGN/SMBH properties under study, in contrast with simple expectations from models of SMBH accretion flows. Our results reaffirm the usability of [Ne v] as an AGN tracer even in highly obscured systems, including dual AGN and high-redshift sources.

Recent hydrodynamical simulations of isolated barred disc galaxies have suggested a structural change in the distribution of the interstellar medium (ISM) around a stellar mass M_* of 10¹⁰ M_⊙. In the higher-mass regime (M_∗ ≥ 10¹⁰ M_⊙), we observe the formation of a central gas and stellar disc with a typical size of a few hundred parsecs connected through lanes to the ends of the stellar bar. In the lower-mass regime (M_∗ < 10¹⁰ M_⊙), such an inner disc is absent and the gas component exhibits a more chaotic distribution. Observations of nearby star-forming galaxies support the existence of such a change. These inner gas discs may represent an important intermediate scale connecting the large kiloparsec-scale structures with the nuclear (sub-parsec) region, transporting gas inwards to fuel the central supermassive black hole (SMBH). For this work we used an extended set of high-resolution hydrodynamical simulations of isolated disc galaxies with initial properties (i.e. stellar mass, gas fraction, stellar disc scale length, and the bulge mass fraction) with properties covering the range of galaxies in the PHANGS sample to investigate this change of regime. We studied the physical properties of the star-forming ISM in both stellar mass regimes and extracted a few physical tracers: the inner Lindblad resonance (ILR), the probability distribution function (PDF), the virial parameter, and the Mach number. In line with observations, we confirm a structure transition in the simulations that occurs between a stellar mass of 10^9.5 and 10¹⁰ M_⊙. We show that the physical origin of this change of regime is driven by stellar feedback and its contribution relative to the underlying gravitational potential. With their shallower potential and typically higher gas mass fraction, lower-mass disc PHANGS galaxies combine two ingredients that significantly delay or even prevent the formation of a central gas (and stellar) disc. These results shed some light on the observed properties of star-forming barred galaxies and have implications for the star formation regimes, the growth of central structures, and the overall secular evolution of disc galaxies.

Effective Field Theories (EFTs) have become indispensable tools for physics beyond the

Standard Model (BSM) to parameterize new physics in a model-independent way. To

obtain consistent predictions, the underlying assumptions have to be specified, and the

EFT has to be applied systematically. For BSM physics, two particular EFTs are widely

used: the Standard Model Effective Field Theory (SMEFT) and the electroweak chiral

Lagrangian, often referred to as the Higgs Effective Field Theory (HEFT). These two

approaches are based on distinct organizing principles. This thesis is dedicated to inves-

tigating the systematic application of both frameworks. For SMEFT, we argue that a power-counting scheme based solely on canonical dimensions is insufficient. It must be supplemented by a loop-order counting scheme, conveniently expressed through the assignment of chiral dimensions. By accounting for both canonical and chiral dimensions, a clear hierarchical structure emerges among the operators, enabling a more meaningful identification of potentially dominant contributions in high-energy processes. As a concrete example, we explore the interplay between these two

counting schemes by matching the Two-Higgs Doublet Model (2HDM) to SMEFT in the

decoupling limit. [...]

The innermost parts of the Milky Way (MW) are very difficult to observe due to the high extinction along the line of sight, especially close to the disc mid-plane. However, this region contains the most massive complex stellar component of the MW, the bulge, primarily composed of disc stars whose structure is (re-)shaped by the evolution of the bar. In this work, we extend the application of the orbit superposition method to explore the present-day 3D structure, orbital composition, chemical abundance trends and kinematics of the MW bulge. Thanks to our approach, we are able to transfer astrometry from Gaia and stellar parameters from APOGEE DR 17 to map the inner MW without obscuration by the survey footprint and selection function. We demonstrate that the MW bulge is made of two main populations originating from a metal-poor, high-α thick disc and a metal-rich, low-α thin disc, with a mass ratio of 4:3, seen as two major components in the metallicity distribution function (MDF). Finer MDF structures hint at multiple sub-populations associated with different orbital families of the bulge, which, however, have broad MDFs themselves. Decomposition using 2D Gaussian Mixture Models in the [Fe/H]-[Mg/Fe] plane identifies five components, including a population with ex-situ origin. Two dominant ones correspond to the thin and thick discs, and two in between trace the transition between them. We show that there is no universal metallicity gradient value that can characterise the MW bulge. The radial gradients closely trace the X-shaped bulge density structure, while the vertical gradient variations follow the boxy component. The MW bulge, while on average having subsolar metallicity, is more metal-rich compared to the surrounding disc populations, in agreement with extragalactic observations and state-of-the-art simulations, reinforcing its secular origin.

Aims. Pulsating variable stars are invaluable tracers for reconstructing the star formation history and chemical evolution of their host galaxies. In this work, we explore the variable star population of Leo II, a distant dwarf spheroidal satellite of the Milky Way. Methods. We analyse an extensive dataset of ground-based BVRI time-series photometry spanning over 35 years. By examining the properties of RR Lyrae stars, we constrain the early chemical enrichment and spatial variation within Leo II. Additionally, we investigate the anomalous Cepheids in order to connect their characteristics with the galaxy's prolonged star formation history, as revealed by deep HST/WFC3 colour–magnitude diagrams. Results. We identified and characterised 175 variable stars, with all but one associated with Leo II. Our work includes the discovery of 25 new RR Lyrae stars and two new anomalous Cepheids. Moreover, we reclassified V88 as a BL Her variable star. By employing multiple independent methods, including metallicity-luminosity relations for RR Lyrae stars and period–luminosity and period–Wesenheit relations for both RR Lyrae and anomalous Cepheids, we derived a true distance modulus of (m ‑ M)₀ = 21.60 ± 0.03 mag, corresponding to 209 ± 4 kpc. Furthermore, we discuss potential formation scenarios for anomalous Cepheids and suggest that, despite the extended star formation history, they are all compatible having originated from old binary stars, with no contribution from young, evolved single stars.

Nonlinear effects in chemical reactions, coupled with amplifying catalysis, can lead to remarkable phenomena like spontaneous symmetry breaking, central to the origin of biological homochirality. Soai's asymmetric autocatalysis is a prototypical reaction for this, where the enantiomeric excess of the product alcohol is amplified during alkylation of pyridyl and pyrimidyl carbaldehydes by diisopropylzinc. However, the complex equilibria and elusive intermediates make the mechanism difficult to clarify. Here we unravel the intricate dynamics of this reaction by in situ high-resolution mass spectrometry, kinetic analysis, and reaction profile simulations. We identify for both the pyrimidyl and the pyridyl systems transient hemiacetalate isopropyl zinc complexes, formed by the addition of the alcoholate product to the aldehyde, as key catalytic intermediates. These diastereomeric complexes enable dual stereocontrol, explaining the observed enantioselectivity. Our analysis confirms the structures of all intermediates and validates the autocatalytic cycle, offering insights into how substituent and structural variations influence reaction performance. This understanding guides the design of new, efficient asymmetric autocatalytic systems.

Red supergiants (RSGs), which are progenitors of hydrogen-rich Type II supernovae (SNe), have been known to pulsate from both observations and theory. The pulsations can be present at core collapse and affect the resulting SN. However, SN light curve models of such RSGs commonly use hydrostatic progenitor models and ignore pulsations. Here, we model the final stages of a 15 solar-mass RSG and self-consistently follow the hydrodynamical evolution. We find the growth of large amplitude radial pulsations in the envelope. After a transient phase where the envelope restructures, the pulsations settle to a steady and periodic oscillation with a period of 817 days. We show that they are driven by the $κγ$-mechanism, which is an interplay between changing opacities and the release of recombination energy of hydrogen and helium. This leads to complex and non-coherent expansion and contraction in different parts of the envelope, which greatly affect the SN progenitor properties, including its location in the Hertzsprung-Russell diagram. We simulate SN explosions of this model at different pulsations phases. Explosions in the compressed state result in a flat light curve (Type II-P). In contrast, the SN light curve in the expanded state declines rapidly, reminiscent of a Type II-L SN. For cases in between, we find light curves with various decline rates. Features in the SN light curves are directly connected to features in the density profiles. These are in turn linked to the envelope ionization structure, which is the driving mechanism of the pulsations. We predict that some of the observed diversity in Type II SN light curves can be explained by RSG pulsations. For more massive RSGs, we expect stronger pulsations that might even lead to dynamical mass ejections of the envelope and to an increased diversity in SN light curves.

Shocks are promising sites of particle acceleration in extragalactic jets. In electron-ion shocks, electrons can be heated up to large Lorentz factors, making them an attractive scenario to explain the high minimum electron Lorentz factors regularly needed to describe the emission of BL Lac objects. Still, the thermal electron component is commonly neglected when modelling the observations, although it holds key informations on the shock properties. We model the broadband emission of the HSP blazar Mrk421 employing particle distributions that include a thermal relativistic Maxwellian component at low energies followed by a nonthermal power-law, as motivated by PIC simulations. The observations in the optical/UV and MeV-GeV bands efficiently restrict the nonthermal emission from the Maxwellian electrons, which we use to derive constraints on the basic properties, such as the fraction $ε_e$ of the total shock energy stored in the nonthermal electrons. The best-fit model yields a nonthermal electron power-law with an index of ~2.4, close to predictions from shock acceleration. Successful fits are obtained when the ratio between the Lorentz factor at which the nonthermal distribution begins ($γ_{\rm nth}$) and the dimensionless electron temperature ($θ$) satisfies $γ_{\rm nth}/θ\lesssim 8$. Since $γ_{\rm nth}/θ$ controls $ε_e$, the latter limit implies that at least $ε_e \approx 10\%$ of the shock energy is transferred to the nonthermal electrons. These results are almost insensitive to the shock velocity $γ_{\rm sh}$, but radio observations indicate $γ_{\rm sh} \gtrsim 5$ since for lower shock velocities the radio fluxes are overproduced by the Maxwellian electrons. If shocks drive the particle energisation, our findings indicate that they operate in the mildly to fully relativistic regime with efficient electron acceleration.

The large $N$ analysis of QCD states that the potential for the $η'$ meson develops cusps at $η' = π/ N_f$, $3 π/N_f$, $\cdots$, with $N_f$ the number of flavors. Furthermore, the recent discussion of generalized anomalies tells us that even for finite $N$ there should be cusps if $N$ and $N_f$ are not coprime, as one can show that the domain wall configuration of $η'$ should support a Chern-Simons theory on it, i.e., domains are not smoothly connected. On the other hand, there is a supporting argument for instanton-like, smooth potentials of $η'$ from the analyses of softly-broken supersymmetric QCD for $N_f= N-1$, $N$, and $N+1$. We argue that the analysis of the $N_f = N$ case should be subject to the above anomaly argument, and thus there should be a cusp; while the $N_f = N \pm 1$ cases are consistent, as $N_f$ and $N$ are coprime. We discuss how this cuspy/smooth transition can be understood. For $N_f< N$, we find that the number of branches of the $η'$ potential is $\operatorname{gcd}(N,N_f)$, which is the minimum number allowed by the anomaly. We also discuss the condition for s-confinement in QCD-like theories, and find that in general the anomaly matching of the $θ$ periodicity indicates that s-confinement can only be possible when $N_f$ and $N$ are coprime. The s-confinement in supersymmetric QCD at $N_f = N+1$ is a famous example, and the argument generalizes for any number of fermions in the adjoint representation.

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

We present simulations of the supernova-driven turbulent interstellar medium (ISM) in a simulation domain of volume (256 pc)³ within which we resolve the formation of protostellar accretion disks and their stellar cores to spatial scales of ~10^-4 au, using the moving-mesh code AREPO. We perform simulations with no magnetic fields, ideal magnetohydrodynamics (MHD) and ambipolar diffusion, and compare the resulting first Larson cores and their associated structures, including the accretion disks, their location within the larger-scale structure and the streamers connecting these. We find that disks of sizes 10 - 100 au form early in the simulations without magnetic fields, while there are no disks larger than 10 au with ideal MHD. Ambipolar diffusion causes large disks to form in a subset of cases (two out of six cores), and generally reduces the strength of outflows, which are seen to play a central role. When they are able to carry away significant angular momentum, they prevent the formation of a rotationally supported disk. Magnetic fields strengths grow from 0.1 - 1 mG in the protostellar core to more than 10 G in the first Larson core in all simulations with ideal MHD. The rotationally supported disks which form can have rotation speeds >1 km s^-1 even out to further than 100 au from the centre, become gravitationally unstable and form complex spiral substructures with Toomre Q < 1. We conclude that the impact of magnetic fields and non-ideal MHD on the formation of protostellar disks is substantial in realistic formation scenarios from the turbulent ISM.

The stellar disc is the dominant luminous component of the Milky Way (MW). Although our understanding of its structure is rapidly expanding due to advances in large-scale surveys of stellar populations across the Galaxy, our picture of the disc remains substantially obscured by selection functions and an incomplete spatial coverage of observational data. In this work, we present the comprehensive chrono-chemo-kinematic structure of the MW disc, recovered using a novel orbit superposition approach combined with data from APOGEE DR 17. We detected periodic azimuthal metallicity variations within 6–8 kpc with an amplitude of 0.05–0.1 dex peaking along the bar major axis. The radial metallicity profile of the MW also varies with azimuth, displaying a pattern typical among other disc galaxies, namely: a decline outside the solar radius and an almost flat profile in the inner region, attributed to the presence of old, metal-poor high-α populations, comprising ≈40% of the total stellar mass. The geometrically defined thick disc and the high-α populations have comparable masses, but with differences in their stellar population content, which we quantified using the reconstructed 3D MW structure. The well-known [α/Fe]-bimodality in the MW disc, once it has been weighted by the stellar mass, is less pronounced at a given metallicity for the whole galaxy but distinctly visible in a narrow range of galactic radii (5–9 kpc), explaining its relative lack of prominence in external galaxies and galaxy formation simulations. Analysing a more evident double age–abundance sequence, we constructed a scenario for the MW disc formation, advocating for an inner and outer disc dichotomy genetically linked to the MW's evolutionary stages. In this picture, the extended solar vicinity is a transition zone that shares the chemical properties of both the inner (old age-metallicity sequence) and outer discs (young age-metallicity sequence).

AT 2019aalc is a peculiar sequence of highly variable emission events observed towards the nucleus of the broad-line active galactic nucleus (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 (BF) 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 BF and the coronal lines vary (roughly) in tandem with these "bumps" in the broadband light curve. Most of the observed behavior of AT 2019aalc links it to the growing class of BF flares 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 (preexisting) AGN accretion disk.

In this exploratory work we analyze the possible existence of a strangeness <inline-formula><mml:math><mml:mrow><mml:mi>S</mml:mi><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>3</mml:mn></mml:mrow></mml:math></inline-formula>, isospin <inline-formula><mml:math><mml:mrow><mml:mi>I</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:math></inline-formula> pentaquark state <inline-formula><mml:math><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>sss</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> generated dynamically from the <inline-formula><mml:math><mml:mrow><mml:mover><mml:mrow><mml:mi>K</mml:mi></mml:mrow><mml:mrow><mml:mo>¯</mml:mo></mml:mrow></mml:mover><mml:mi>Ξ</mml:mi></mml:mrow></mml:math></inline-formula> interaction. We employ a unitarized scheme in coupled channels based on the chiral Lagrangian expanded up to next-to-leading order (NLO), and show that the inclusion of the NLO terms can be crucial to provide the necessary attraction that favors the existence of such triply strange pentaquark. The <inline-formula><mml:math><mml:mrow><mml:mover><mml:mrow><mml:mi>K</mml:mi></mml:mrow><mml:mrow><mml:mo>¯</mml:mo></mml:mrow></mml:mover><mml:mi>Ξ</mml:mi></mml:mrow></mml:math></inline-formula> femtoscopic correlation functions are calculated as example of a possible experimental measurement in which a direct signal of the <inline-formula><mml:math><mml:msub><mml:mi>P</mml:mi><mml:mrow><mml:mi>sss</mml:mi></mml:mrow></mml:msub></mml:math></inline-formula> state could be observed.

Context. Dark matter (DM) particles can interact with particles characterised by the standard model. Although there are a number of constraints derived from direct and indirect detection experiments, the dynamical evolution of astrophysical objects could offer a promising probe for such interactions. Obtaining astrophysical predictions is challenging and primarily limited by our ability to simulate scattering between DM and baryonic particles within N-body and hydrodynamics simulations. Aims. We have developed the first scheme allowing for the simulation of these interacting dark matter (IDM) models, accurately accounting for their angular and velocity dependence, as well as the mass ratio between the DM and baryonic scattering partners. Methods. To describe DM-baryon interactions, we used an N-body code together with its implementation of smoothed-particle hydrodynamics (SPH) and meshless finite mass. The interaction itself was realised in a pairwise fashion by creating a virtual scattering partner from the baryonic particle and allowing it to interact with a DM particle using a scattering routine initially developed for self-interacting dark matter (SIDM). After the interaction, the virtual particle is rejoined with the baryonic particle, fulfilling the requirements of energy and momentum conservation. Results. Through several test problems, we demonstrated that we are able to reproduce the analytic solutions with our IDM scheme. This includes a test for scattering with a physical mass ratio of 1:1000, which is beyond the limits of current SIDM simulations. We comment on various numerical aspects and challenges, and we describe the limitations of our numerical scheme. Furthermore, we study the impact of IDM on halo formation with a collapsing over-density. Conclusions. We find that it is possible to accurately model IDM within N-body and hydrodynamics simulations commonly used in astrophysics. Finally, our scheme allows for novel predictions to be made and new constraints on DM-baryon scattering to be set.

We present a neural-network framework designed to reconstruct the properties of cosmic-ray nuclei traversing the scintillating-fiber tracking calorimeter of the RadMap Telescope. Employing the Geant4 simulation toolkit and a simplified model of the detector to generate training and test data, we achieve the spectroscopic capabilities required for an accurate determination of the biologically relevant dose that astronauts receive in space. We can reconstruct a particle's trajectory with an angular resolution of better than $1.4^\circ$ and achieve a charge separation of better than $95\%$ for nuclei with $Z\leq8$; specifically, we reach an accuracy of $99.8\%$ for hydrogen. The energy resolution is $<20\%$ for energies below 1 GeV/n and elements up to iron. We also discuss the limitations of our detector, the reconstruction framework, and this feasibility study, as well as possible improvements.

We present a microscopic model of the dark sector that resolves the Hubble tension within standard current datasets based on well-known fundamental principles, gauge symmetry and spontaneous symmetry breaking. It builds on the Hot New Early Dark Energy (Hot NEDE) setup, featuring a dark $SU(N)$ gauge symmetry broken to $SU(N-1)$ in a supercooled phase transition that creates a thermal bath of self-interacting dark radiation in the epoch between Big Bang Nucleosynthesis and recombination. Adding a fermion multiplet charged under the gauge symmetry provides a naturally stable component of dark matter that interacts with dark radiation. Spontaneous symmetry breaking predicts a decoupling of this interaction once the dark sector cools down, that we refer to as dark radiation matter decoupling (DRMD). We find agreement between the SH${}_0$ES determination of $H_0$ as well as combined Planck 2018, Pantheon+ and DESI baryon acoustic oscillation (BAO) data at 1.4$σ$ level, compared to a 5.7$σ$ tension in the $Λ$ Cold Dark Matter model. We also provide a simplified three-parameter DRMD model encoding the essential features, while the full model offers additional falsifiable predictions.