2026-05-06 | | Total: 75
We present a data-driven fuzzy set framework for classifying galaxies into the red sequence, blue cloud, and green-valley populations using multiple observables from the Sloan Digital Sky Survey (SDSS DR18). Unlike traditional methods based on hard boundaries in colour or stellar mass, our approach assigns continuous membership degrees using sigmoidal functions derived from bimodal galaxy properties, including $(u-r)$ colour, specific star formation rate (sSFR), and $D4000$. Membership functions are constructed via Gaussian mixture modeling and combined using a conservative fuzzy minimum operator. Applying this method to a volume-limited sample of 88,579 galaxies, we compare with the empirical classification of \citet{schawinski14}. The fuzzy approach reduces contamination in the red and green-valley populations and yields more physically consistent distributions of star formation and morphology. Red galaxies show a unimodal low-sSFR distribution, while green-valley galaxies exhibit clearer signatures of morphological evolution. We also examine the dependence of active galactic nucleus (AGN) fraction on stellar mass and find no significant differences between methods, indicating robust global AGN trends. However, clustering analysis reveals subtle differences: fuzzy-classified red galaxies show enhanced large-scale clustering, suggesting a stronger association with highly biased dark matter halos. These results demonstrate that fuzzy classification provides a flexible, physically motivated alternative to hard-cut methods, enabling a more accurate and interpretable view of galaxy populations and their evolution.
The largest known molecule in space, C60 , has been detected in its neutral and cationic form through its vibrational, UV-driven fluorescence emission spectrum and its electronic absorption spectrum, respectively. The detection of several polycyclic aromatic hydrocarbon molecules through their pure rotation spectrum in cold, dense, molecular cloud cores suggests that C60 might be present in these environments as well. The low flux of UV pumping photons in molecular cloud cores and the absence of suitably bright background stars, make detection of C60 and its cation through the commonly used methods impractical. As C60 has no permanent dipole moment, its pure rotational transitions are forbidden and its presence must be inferred from the rotational transitions of C60 derivatives with permanent dipole moments. Here, we present a study of the predicted rotational spectrum of protonated C60 that has a sizeable permanent dipole moment. Protonation of C60 reduces the icosahedral symmetry to Cs and results in a dipole moment of about 3.8 Debye. The resulting C60H+ is a closed shell system
Stellar abundance analysis relies on flexible, high-performance spectral synthesis. To meet these needs, we present PySME v1.0, an updated Python implementation of Spectroscopy Made Easy (SME) designed for precise and survey-scale modelling of stellar spectra.A central challenge in SME based synthesis is the efficient treatment of very large line lists, including both the preselection of negligible lines and the subsequent formal synthesis. PySME v1.0 introduces a revised line-selection framework based on opacity ratio and line depth, together with dynamic line list construction and control of the effective wavelength span over which each line contributes to the synthetic spectrum. These workflows support parallel preprocessing of weak-line selection and reduce the line list passed to the synthesis core, thereby improving scalability while preserving synthetic accuracy. PySME v1.0 also incorporates an updated equation-of-state treatment that improves the modelling of hydrogen lines, particularly Balmer features, while maintaining close agreement with previous SME results for metal lines. The Python interface has further been extended to support parameter-dependent derived quantities updated during optimisation, and PySME provides non-local thermodynamic equilibrium (NLTE) departure-coefficient grids for 17 elements. Together, these developments establish PySME v1.0 as a robust and efficient framework for high-precision stellar abundance analyses in large spectroscopic surveys.
Past work has shown that ocean salinity and planetary obliquity both influence the climates of Earth-like exoplanets throughout the habitable zone of Sun-like stars. The effects of salinity and obliquity can be profound, with low vs. high salinity or obliquity resulting in distinct climate states in some scenarios. However, past work has considered salinity or obliquity in isolation and has not explored how each may modulate the effects of the other. We investigate how ocean salinity and planetary obliquity jointly impact climate and habitability using the ROCKE-3D coupled ocean-atmosphere general circulation model. We find that salinity and obliquity have a greater combined impact on planetary climate than the sum of their effects in isolation. This synergy between salinity and obliquity arises due to the ice-albedo feedback, producing distinct climate states that range from ice-free to globally glaciated while having same initial atmospheric conditions and receiving the same instellation. Consequently, ocean salinity and planetary obliquity can together lead to divergent habitability outcomes for otherwise identical planetary scenarios and initial conditions. Salinity and obliquity can jointly increase the planetary fractional habitability across oceans and continents, especially for cold exoplanets. Although neither ocean salinity nor planetary obliquity can be reliably predicted or observationally constrained, their synergistic effects must be considered in future studies of planetary climate and exoplanet observations, especially when characterizing planetary habitability.
We compare images of core-collapse supernova (CCSN) remnants (CCSNRs) and jet-shaped planetary nebulae (PNe) that have a narrow, faint zone extending from side to side, termed a pipe, with a hydrodynamical numerical simulation exploding a massive star with three pairs of jets in the framework of the jittering jets explosion mechanism (JJEM), and conclude that jets shaped the pipes in these CCSNRs and PNe. We present two jet-shaped PNe with a pipe and three PNe with two opposite narrow jet-shaped lobes, and argue that in some cases the two opposite narrow lobes might merge to form one long, faint zone extending from side to side of the PN, namely, a pipe. From the qualitative similarity of the pipe morphology of the two CCSNRs we analyze with the pipe of the PNe, we suggest that jets also shaped the pipe of these CCSNRs. We strengthen this conclusion with a three-dimensional hydrodynamic simulation that reproduces two opposite narrow lobes, similar to those observed in PNe with lobes. These lobes can merge later to form a pipe. This paper is another in a series that strengthen the case for the JJEM as the primary explosion mechanism of CCSNe by comparing CCSNR morphologies with those of jet-shaped PNe.
Neutron-star mergers can launch mildly relativistic to moderately relativistic outflows whose interaction with the ejecta can reshape kilonova emission. We parametrically study magnetically powered outbursts from long-lived merger remnants, such as flare-like eruptions and collapse-driven shocks, and quantify their impact on ejecta dynamics, composition, and observables. Using two-dimensional special-relativistic magnetohydrodynamic simulations, we follow magnetized blast waves injected into expanding merger ejecta for early- and late-launch scenarios across a range of shock strengths. We then post-process Lagrangian tracer histories with the nuclear reaction network WinNet and the radiative-transfer code SuperNu with realistic opacities, to connect shock heating directly to nucleosynthesis and kilonova light curves. We find that sufficiently strong shocks can reheat portions of the ejecta to nuclear statistical equilibrium, increase the electron fraction in the shocked material, and deposit entropy, leading to systematic changes in $r$-process yields. These thermodynamic and compositional changes can leave observable imprints on kilonova emission -- especially in color evolution and late-time light-curve behavior -- indicating that magnetically driven remnant variability can potentially contribute to kilonova diversity.
We investigate the plausible history of the XUV luminosity evolution of the planet-hosting M4 star GJ 1132 (~0.2 solar masses) to infer the cumulative incident XUV flux intercepted by the short-period (~1.6 d) Earth-sized transiting planet GJ 1132 b. We include the dominant observational uncertainties, compare two quiescent XUV luminosity evolution models, and simulate the XUV luminosity evolution from flares based on TESS data and a re-analysis of Kepler stars. We find only 4 flares in GJ 1132's TESS 123 day lightcurve, which is relatively few for M dwarfs and, in conjunction with the ~125 day period, suggests that this star is many Gyr old. We find that all model permutations predict that the planet has at least a 95% chance of receiving more than 50 times as much XUV flux as modern Earth, confirming that this planet is a good candidate for permanent atmospheric loss. We also find that an empirical XUV model for M dwarfs predicts 2-3 times more total XUV flux than a commonly used solar twin model and that the empirical model's distribution is 2-3 times narrower. Flares contribute about 20% of the cumulative XUV flux on planet b, which, while modest, ensures the planet lies firmly on the atmosphere-free side of the "cosmic shoreline."
The power spectrum (PS) of the density field in supersonic turbulence is a fundamental quantity that characterizes the statistical properties of the structures formed in compressible flows. It is also widely used to estimate the Mach number in the interstellar medium from simulation-derived relations. We provide here a first quantitative explanation for the evolution of the slope of the PS of the density field with the Mach number in homogeneous isotropic isothermal turbulence using a time-invariant quantity derived by Chandrasekhar (1951). For simulated turbulent flows, the model reproduces very well the measured slopes for different widths of the inertial range and density variances. Our model also provides a comprehensive interpretation of the characteristic slopes of the PS of the density field measured in the interstellar medium. Based on these results, we stress that the Mach number cannot be reliably deduced from the slope of the PS of the density field. We finally discuss a resolution criterion that must be fulfilled to correctly simulate a turbulent flow with a given density PS slope.
Quasar microlensing is both a very useful tool in cosmology and astrophysics, and a source of uncertainty in some studies like the determination of the Hubble constant from lensed quasars. Microlensing probability and time-scales have been statistically studied using as a reference scale the Einstein ring crossing time of an isolated mass. Our goal is to extend the statistical analysis of microlensing to all currently known lensed quasars with available data, considering realistic optical depths and the gravitational effect of the lens galaxy. We take into account new observational results about quasar sizes and peculiar velocities of lens galaxies. We apply automatic lens modeling to the 204 systems available. For each image, we compute microlensing magnification maps and histograms. Using thin disk source sizes scaled to take into account recent measurements of accretion disk sizes, we find a mean source crossing time of $2.59\pm 0.07$ years. The mean Einstein radius crossing time is $ 11.29 \pm 0.05$ years. When a fraction of mass in microlenses $α=0.2$ is adopted, we find a good matching between the modeled histogram of mean microlensing magnifications for the images in our sample and the experimental histogram of microlensing magnifications. From the modeling of microlensing magnification histograms, we estimate the average half-light radius of the quasar source, $R_{1/2}=5.4\pm 2.7$ light-days, and a lower limit to the mass fraction in microlenses, $α\ge 0.15$. From the microlensing magnification maps, we find that a lensed quasar image has a mean probability of approximately 9% of being involved in a high-magnification event ($Δm \le -0.32$). We select a group of images with the largest probabilities and the smallest crossing times.
Disk-mediated accretion is central to theories of massive star formation, setting the initial conditions for their evolution. Yet observations of Keplerian disks around early O-type protostars remain scarce, as they are often blended into complex surrounding structures. We report ALMA Band 6 observations (300 au resolution) of an accretion disk surrounding a high-mass protostar in the Sagittarius C (Sgr C) cloud in the Central Molecular Zone (CMZ) around the Galactic Center. We identify spectral lines and analyze the spatial distribution of the emission of the complex organic molecules. We use a dynamical model with an inner Keplerian disk and an outer free-fall envelope to fit the three-dimensional position-position-velocity data of the stacked CH$_3$OCHO molecular lines and constrain the mass of the central protostar to be $\sim40^{+2}_{-3} M_{\odot}$. The fitting results additionally show that the disk has a centrifugal radius at about 1300 au. Considering the infall velocity, radius, and mass of the envelope, we estimate the accretion rate from the envelope onto the disk to be $\sim7\times 10^{-3}\ M_{\odot}\,\mathrm{yr^{-1}}$. We also identify spiral-like structures in the disk that can be described by free-falling streamers. Our results highlight the critical role of accretion disks and streamers in the mass accumulation of early O-type stars in the CMZ.
New basis functions for 2 dimensional (2D) image analysis with a circular boundary (referred to as multipole analysis) are derived which are equal-area weighted. We present open access Python code hosted by GitHub, with which users can apply the multipole analysis to images. The new multipole analysis is applied to a set of 28 supernova remnants (SNRs) which are selected to have both radio and X-ray images, and have been identified as Type Ia or Type CC. Each pair of SNR images (radio and X-ray) was convolved to the same spatial resolution prior to analysis. The resulting multipole radial powers and angular powers, from order 0 to 5, for a given SNR are different for different multipoles and for a given multipole are different between X-ray and radio images. The X-ray radial powers (for orders >0) are larger on average than the radio radial powers (more radial structure in X-rays than radio). The angular powers are smaller than the radial powers on average (more radial structure than angular structure). Comparing Type Ia and Type CC populations, the radial powers (for orders >0) are on average larger for Type CC than Type Ia for X-ray and radio images, with larger difference for X-ray images. The angular powers (for orders >0) are similar between Type Ia and Type CC for both radio and X-ray images.
Context. Ellerman bombs (EBs) are observational signatures of small-scale magnetic reconnection, key to understanding the lower solar atmosphere. While their role in active regions has been widely studied using the H$α$ line, near-ultraviolet (NUV) spectra routinely observed by the Interface Region Imaging Spectrograph (IRIS) offer a promising alternative for EB identification, enabling large-scale studies. Aims. We aim to identify the most important spectral signatures of EBs in the IRIS NUV spectra. With this, we seek to develop a robust criterion for their detection solely using the IRIS NUV spectra. In parallel, we determine the typical atmospheric stratification associated with EBs. Methods. We used four coordinated observations between the Swedish 1-m Solar Telescope (SST) and IRIS. Using the H$α$ line as a reference, we detected 18 different EBs and studied their associated IRIS NUV spectra. In addition, we used the IRIS$^{2+}$ inversion tool to infer the temperature, line-of-sight velocity, and non-thermal broadening from the EB spectra. Results. The defining feature of EBs in the IRIS NUV is the enhancement of the wings of the subordinated Mg II triplet in between the Mg II h&k lines. Inversions reveal that these signatures are produced by localized temperature increase of $Δ$ T~1650 K around log$τ$=-3.8. Using only the Mg II triplet signatures, we found a detection criterion that successfully recovered 14 of 18 H$α$-detected EBs. In addition, the shape of the Mg II h&k lines in relation to the Mg II triplet can serve as a proxy for the EB formation height. Conclusions. The NUV spectrum observed by IRIS is a good candidate for detecting EBs, opening the doors to large-scale studies across the extensive IRIS database, removing the dependence on H$α$ observations.
Thermohaline convection is a standard chemical mixing process in stellar interiors, yet its mixing efficiency is not fully settled. Competing theories predict turbulent diffusion coefficients, $D_μ$, that can differ by orders of magnitude, leading to uncertainties in stellar models and interpretations of observations. This paper explores a potential resolution to existing discrepancies. We first complete the linear stability theory and identify two types of unstable modes: slow growing modes at large length scales and fast growing modes at small length scales. We then reevaluate $D_μ$ considering the full spectrum of unstable modes and find that it can self-consistently interpolate between previously proposed theoretical scalings across the instability parameter space. The question of thermohaline mixing efficiency in stars may be settled by future simulations that quantify the scale-dependent contributions of fast and slow modes to $D_μ$ and determine how the modes dominating the transport change across parameter space.
We present new constraints on the local-type primordial non-Gaussianity parameter, $f_\mathrm{NL}^\mathrm{local}$, through analysis of the scale-dependent bias effect on the cosmic infrared background (CIB). To avoid biases from galactic dust contamination on large scales, we use cross-correlations between the CIB and Planck cosmic microwave background (CMB) lensing maps to constrain non-Gaussianity. Our measurement employs new dust-cleaned CIB maps that have been designed to be unbiased on large scales, which allows us to improve our constraining power on $f_\mathrm{NL}^\mathrm{local}$ by a factor of $\sim 2$ over previous CIB analyses. We derive a constraint of $f_\mathrm{NL}^\mathrm{local}=43 \pm 23$, matching the precision of the tightest existing constraints from cross-correlation methods. Consistency- and null-tests demonstrate that our results are robust to modeling assumptions and residual dust contamination.
We present VERSUS, a publicly available, fast void-finding algorithm designed to identify spherical underdensities in the density field that can be accurately described by excursion set predictions of the void size function. We validate the algorithm against both a synthetic distribution of particles designed to trace a known input void population, and mock galaxy sample built from a $(2\ h^{-1}\text{Gpc})^3$ AbacusSummit simulation populated with a realistic galaxy-halo connection, including systematic effects designed to mimic real survey data. In all cases, VERSUS demonstrates excellent performance, achieving strong agreement with theoretical predictions for the void size function across the range $25 < R \,[\ h^{-1}\text{Mpc}] < 61$ without requiring any post-processing of the void catalogue. The code is user-friendly, modular, and readily applicable to observational survey data. Its computational efficiency further enables the use of simulation-based modelling approaches, facilitating robust and consistent cosmic void analyses with Stage-IV surveys.
We compute the thermal activation rate of metastable self-gravitating Bose-Einstein condensates with attractive self-interaction (e.g., dilute axion stars) by using the instanton theory. Explicit analytical results are given close to the maximum mass $M_{\rm max}$ [P.H. Chavanis, Phys. Rev. D 84, 043531 (2011)] by using the normal form of the saddle-node bifurcation close to that point. We show that the lifetime of metastable states is extremely long, scaling as $t_{\rm life}\sim e^N\, t_D$, where $N$ is the number of bosons in the system and $t_D$ is the dynamical time ($N\sim 10^{57}$ and $t_D\sim 10\, {\rm hrs}$ for typical QCD axion stars; $N\sim 10^{96}$ and $t_D\sim 100\, {\rm Myrs}$ for the quantum core of a dark matter halo made of ultralight axions). Therefore, metastable equilibrium states can be considered as stable equilibrium states in practice. We compare our results with similar results obtained for Bose-Einstein condensates in laboratory, globular clusters and self-gravitating Brownian particles in astrophysics, the Brownian mean field model (BMF) in statistical mechanics, and bacterial populations in biology. Our presentation parallels the calculation of the quantum tunneling rate of dilute axion stars given in a previous paper [P.H. Chavanis, Phys. Rev. D 102, 083531 (2020)]. These calculations can find application in various domains of physics and astrophysics.
The Strouhal number, ${\rm{St}}=t_{\rm cor}/t_{\rm out}$, measures the temporal coherence of turbulent driving relative to the outer-scale eddy turnover time. In turbulence-box models one commonly sets ${\rm{St}}=1$, although recent work by \citet{Grete2025_density_distribution} and \citet{Scannapieco2025_density_distribution} has shown that turbulence statistics, especially the mass-density distribution in compressively driven turbulence, are sensitive to this choice. In this Letter, we compute ${\rm{St}}$ directly from the measured two-time correlation tensor and outer-scale eddy time in stratified multiphase ISM simulations of Milky Way-like and starburst disks. We find isotropic median values ${\rm{St}}=0.26^{+0.30}_{-0.16}$ for the Milky Way-like model and ${\rm{St}}=0.25^{+0.11}_{-0.12}$ for the starburst model. These values are consistent with the picture that supernova remnants (SNRs) drive turbulence locally near $R_{\rm cool}$, where the unstable contact discontinuity in the expanding SNR sets comparable forcing and eddy times, ${\rm{St}}(R_{\rm cool})\approx 1$. The reconstructed scale-dependent curves reach ${\rm{St}}=1$ at a nearly universal outer-scale fraction, $\ell_\ast/\ell_{\rm out}\approx0.12\text{--}0.13$ ($\ell_\ast\approx25\text{--}32\,\rm{pc}$), so the standard ${\rm{St}}=1$ prescription is not an outer-scale model of SN-driven ISM turbulence, but a local-scale approximation tied to injection near the cooling radius of the SNR.
Modeling the interior of a rocky or water-rich exoplanet is a thermodynamic closure problem: every layer's density, temperature gradient, and phase must follow from an equation of state (EoS) that remains self-consistent across the pressure-temperature range from surface to core. Existing EoS span disciplines, use different formalisms, and rarely supply the full thermodynamic quantities needed by evolutionary models of interior phase transitions. We present PALEOS (Planetary Assemblage Layers: Equations of State), an open-source toolkit consolidating EoS for iron, magnesium silicate (MgSiO$_3$), and water (H$_2$O) into a unified, phase-aware, thermally responsive framework spanning 17 phases. PALEOS derives density, energy, entropy, heat capacities, thermal expansion, and the adiabatic gradient analytically via Maxwell relations, and is released as lookup tables on regular P-T grids. We validate it against the Preliminary Reference Earth Model, recovering Earth's radius to 0.3% and lower-mantle densities to 3%, and compute 17,900 mass-radius relations from 0.1 to 100 $M_\oplus$ for rocky (Fe + MgSiO$_3$) and water-rich (Earth-like core + H$_2$O envelope) compositions at 300-4000 K. Continuous solid-to-melt EoS let thermal expansion span the fully-solid to magma-ocean regime: the radius offset exceeds 1% above 1500 K and reaches 16% at 4000 K for low-mass silicate planets, comparable to composition degeneracy and transit-radius uncertainties. We demonstrate this on two ultrashort-period super-Earths, WASP-47 e and TOI-1807 b: each admits two purely rocky solutions indistinguishable in mass and radius but in radically different states, one fully solid with no dynamo, the other hosting a deep magma ocean and a liquid iron core capable of sustaining a magnetic field. Phase-aware, thermally resolved EoS are essential for translating astronomical observations into exoplanetary geophysics.
Observations indicate that the total abundance of S-bearing species in dense clouds is orders of magnitude lower than the cosmic sulfur abundance. Addressing this "missing sulfur problem" requires a combination of astronomical observations, laboratory experiments, and theoretical models. In this work, we use the pyRate astrochemical model to simulate the VUV photon irradiation of a CO$_2$:CS$_2$ ice mixture at 10 K in the laboratory, with the goal of supporting the interpretation of the experimental results and testing our current understanding of the sulfur evolution in interstellar ices. For this purpose, the astrochemical model was adapted to the experimental conditions, and the chemical network was compiled from several sources to ensure that all known reactions involving sulfur species were included. The results indicate that nondiffusive chemistry is necessary to reproduce the formation of S-bearing species observed in the experiment. However, some discrepancies were found in the major S-bearing ice chemistry products predicted by the model and the experiment. The compounds OCS, CS, and SO are overpredicted by the model, while it falls short in accounting for $\rm SO_2$ and sulfur allotropes. These discrepancies are likely due to a combination of an incomplete knowledge of the chemical reactions at play (either because of missing reactions and/or because of unconstrained reaction barriers), and uncertainties in the experimental analysis. This work represents the first effort to model the chemistry of a multicomponent ice analog with a rate-equation based code, and highlights the complementary nature of theoretical and experimental astrochemistry to disentangle the chemical evolution of sulfur in the interstellar medium.
The ESA Gaia mission is a 10+ year astrometric whole-sky scan, demanding consistent data quality over the whole timespan of operations Aims. The Gaia First Look (FL) is a system whose aim is monitoring the data quality to identify problems, which includes early warning capabilities for potential upcoming issues. Methods. In order to achieve its goals, the Gaia FL implemented its own limited astrometric solution, and used the daily calibrations from other segments of the Data Processing and Analysis Consortium (DPAC), as well as the diagnostic data from the satellite itself, in order to obtain a complete picture of the situation of the Gaia satellite on a daily basis. This led to a short-term health and data quality check, but also to a broader overview of the longer-term trends and evolutions within the payload. Potential issues that were encountered were reported to other groups within DPAC for further analysis purposes. When required, ways to mitigate the problems were discussed, and implemented. Results. We show a number of findings by the Gaia FL concerning longer-term evolution, individual but common effects, as well as detrimental impacts, all of which occurred over the operational phase of the Gaia mission
Detecting and characterising exoplanet atmospheres remains challenging because atmospheric signals can be comparable to residual noise and instrumental/astrophysical systematics. Spectral features span from a few ppm for small planets up to $\sim 10^3$ ppm for warm/hot giants, while high-quality JWST time-series spectroscopy typically reaches $\sim 10$--$50$ ppm (occasionally $\sim 100$--$200$ ppm in the presence of stellar variability or stronger systematics), making correlated noise across temporal and spectral dimensions a key limitation. With JWST delivering an increasing volume of high-precision transmission spectra, and Ariel set to extend this to a homogeneous survey of $\sim 10^3$ exoplanet atmospheres, robust benchmarking resources with known ground truth are essential to develop and validate data-driven (including ML-based) detrending approaches. As a major step towards this goal, we use ExoSim2 and TauREx to generate one of the most comprehensive public datasets based on the current payload design of the ESA Ariel mission, specifically intended to benchmark detrending algorithms. We also provide a deep neural network baseline for time-series reduction, and use it to highlight the limitations of ML based detrendng methods, i.e. the risks posed by dataset shift when observed distributions diverge from those of the training set, a scenario likely to arise in real observations. This dataset is featured in the Ariel Data Challenge 2024 on Kaggle and has been field-tested for robustness and simulation fidelity. By making these resources publicly available, we aim to support the community in developing, comparing, and stress-testing scalable and reliable methods for exoplanet transmission spectroscopy.
Context. Green Pea galaxies (GPs) are local starburst galaxies serving as analogues for high-redshift star-forming galaxies, particularly Lyman continuum leakers. It remains debated whether their starbursts are driven by internal secular processes or external triggers. Aims. We aim to constrain the role of environment in this triggering, testing whether external influence comes from close interactions or diffuse processes like gas accretion. Methods. We analyse VLT/MUSE observations of 24 GPs at $z \sim 0.2$ to identify companions via spectral line features. We derive key physical properties (extinction, SFR, stellar mass, age, metallicity) for GPs and companions, and estimate group dynamical masses. Results. We identify 22 emission-line galaxies, 11 being companions ($|Δv| \leq 500$ km s$^{-1}$). We find a high companion fraction ($33^{+11}_{-8}$%) and a $\sim$1 dex number density excess compared to the field, confirming GPs reside in overdense environments. Companions typically lie at projected separations of $\sim$100 kpc with no evidence of ongoing interactions. Physically, GPs form a homogeneous class of young (mass-weighted age $\sim$230 Myr), metal-poor, high-sSFR starbursts with elevated velocity dispersions. In contrast, companions are more evolved ($\sim$1.6 Gyr) and heterogeneous in stellar mass, metallicity, and dust attenuation. Inferred group dynamical masses are $\sim$3 dex higher than total stellar masses, suggesting significant dark matter and neutral gas. Conclusions. GPs do not appear triggered by ongoing major mergers with close (10-30 kpc) companions. Results favor a scenario where GPs are transient starbursts in overdense regions, plausibly sustained by gas accretion. Limited spatial resolution prevents ruling out very close mergers ($\lesssim 10$ kpc). High dynamical-to-stellar mass ratios imply substantial non-stellar mass in these systems.
JWST's discovery of unexpectedly bright $z>10$ galaxies has triggered claims that standard $Λ$CDM cannot reproduce their abundances, while estimates of the ionizing escape fraction $f_{\rm esc}$ at $z>6$ have spanned a factor of four for over a decade. Here we show that both tensions arise from a structural degeneracy in reionization equations: global observables constrain only the product $f_{\rm esc}\times f_{\star,0}$ (peak star formation efficiency), not individual parameters. We demonstrate that this degeneracy, previously considered a limitation, provides a precise diagnostic framework. By leveraging JWST UV luminosity function shapes to independently constrain $f_{\star,0}$, we derive robust bounds on $f_{\rm esc}$. Joint profile-likelihood analysis across Gaussian, log-normal, and duty-cycle burst scatter models excludes the proposed crisis threshold ($\varepsilon > 3.5\%$) at $4.5σ$ confidence, with stochastic star formation histories strengthening rather than weakening the result. Combining these constraints with constant and evolving $f_{\star,0}$ measurements yields the first empirical reconstruction of $f_{\rm esc}(z)$ across $z=7$--$12$. A constant-efficiency scenario ($f_{\rm esc}\approx 10$--$16\%$) connects smoothly to low-redshift direct detections, whereas an evolving scenario ($f_{\rm esc} \approx 6\%$ at $z=12$) conflicts with low-metallicity ISM porosity expectations. JWST Cycle 3--4 will distinguish these pathways at $>2σ$, transforming a long-standing fundamental inference barrier into a powerful quantitative probe of early-universe physics.
Among the few ways to probe the early Universe, neutrinos offer a particular window on high energy phenomena occurring before recombination. We discuss the opportunities of observing primordial high energy neutrinos (Phenus): neutrinos produced before or around recombination from the decay or annihilation of long-lived relics, arriving at detectors today with energies in the GeV-PeV range. We summarise the results of a general study of this scenario, covering the sharp spectral features such fluxes would display, the theoretical (BBN and CMB) and experimental constraints on the source particle parameter space, and the regions that could realistically be probed by current and future neutrino telescopes. We also present a dedicated Monte Carlo code for computing the distortion of the Phenu spectrum by final state radiation and interactions with the cosmic neutrino background during propagation, and apply it to assess the primordial origin hypothesis for the KM3-230213A ultrahigh energy neutrino event.
Robust reddening maps of the Large and Small Magellanic Clouds (LMC/SMC) are crucial for a wide range of astrophysical studies, including the calibration of the cosmic distance ladder, investigations of stellar populations in low-metallicity environments, and the characterization of interstellar dust properties. We aim to construct reddening maps of the Magellanic Clouds using spectral energy distribution (SED) fitting, and to investigate the impact of different stellar atmosphere models on the resulting maps. We combined optical ($ugriz$) photometry from the SMASH survey with near-infrared ($YJK_{\rm s}$) photometry from the VMC survey for red giant branch (RGB) stars. Observed SEDs were matched to synthetic photometry derived from three atmosphere model grids. Our maps cover 34.5 deg$^2$ of the LMC and 24.5 deg$^2$ of the SMC at 4 arcmin resolution. We find mean reddening values of $E(B-V)=0.076 \pm 0.022$ mag for the LMC and $0.058 \pm 0.024$ mag for the SMC. We found that employing different atmospheric models results in differences up to 0.03 mag in the mean reddening. Canonical $R_V$ values for the Magellanic Clouds (3.41 for LMC and 2.74 for SMC, Gordon et al. 2003) provide results consistent with previous studies. We confirm higher and more structured reddening in the LMC compared to the SMC, with 30 Doradus standing out as the dominant high-reddening region. Our results show that the absolute reddening scale depends on the choice of stellar atmosphere models, while the relative spatial structure of the reddening maps remains stable.