Papers for Wednesday, Mar 19 2025

A critical issue in studying the evolution of galaxy clusters is to find ways that enable meaningful comparisons of clusters observed at different redshifts, as well as in various stages of their growth. Studies in the past have typically suffered from uncertainties in cluster mass estimates due to the scatter between cluster observables and mass. Here we propose a novel and general approach that uses the probability distribution function of an observable-cluster mass relation, together with dark matter halo merger trees extracted from numerical simulations, such that one can trace the evolution in a self-contained fashion, for clusters chosen to lie in a specified range in mass and redshift. This method, when applied to clusters of different mass ranges, further allows one to examine the evolution of various observable-cluster mass scaling relations. We illustrate the potential of this method by studying the stellar mass content of red cluster member galaxies, as well as the growth of brightest cluster galaxies, from z=1.0 to z=0.3, using a large optically-detected cluster sample from the Subaru Hyper Suprime-Cam Survey, finding good agreement with previous studies.

We investigate the fundamental plane (FP) of selected early-type (ETG) member galaxies of the galaxy cluster PLCK G287.0+32.9 ($ z_c = 0.3833 $), exploring also four-dimensional hyperplane extensions. We measure ETGs structural parameters and photometry from Hubble Space Telescope (HST) observations. We use high-quality spectroscopic data from the Multi Unit Spectroscopic Explorer (MUSE) to measure the galaxy central stellar velocity dispersions and stellar population properties. With this data, we construct the FP through a robust fitting procedure and analyze its tilt and scatter. We then introduce two hyperplane extensions, one including the stellar mass ($M^\star$-HP) and another including the stellar over total mass fraction ($f_{\mathrm{e}}^\star$-HP), and compare their coefficients and scatter to those of the FP. The FP of PLCK G287.0+32.9 is found to have best-fit parameter values consistent with those in the literature ($\alpha = 1.2 \pm 0.1$ and $\beta = -0.75 \pm 0.05$), with a scatter of $0.09$ dex. The ($f_{\mathrm{e}}^\star$-HP) shows no tilt compared to the theoretical plane ($\alpha = 2.1 \pm 0.2$ and $\beta = -1.12 \pm 0.07$), with a scatter of $0.042$ dex, and the ($M^\star$-HP) reveals an even tighter relation, with a scatter of only $0.023$. Our findings support the idea that the FP is a lower-dimensional projection of a more complex hyperplane and confirm that the variations in the dark matter content contribute significantly to the tilt of the FP. Future studies incorporating larger samples of galaxies and additional physical parameters may further refine our understanding of the FP and its higher-dimensional extensions.

Observations indicate that central galaxies show a significant alignment of their main shape axes with other galaxies in their group, as well as with the large-scale structure of the universe. Simulations have corroborated this finding, providing further insights into how the shape of the stellar component aligns with the surrounding dark matter halo. Recent studies have also investigated the evolution of this alignment in bright central galaxies, revealing that the shapes of the dark matter halo and the stellar component can differ. In this work, we aim at gaining a deeper understanding of galaxy alignments by quantifying how this property is related to the mass of the haloes hosting central galaxies and to the large-scale environment measured at different scales. By studying different angles, we describe how the alignments of central galaxies depend on the mass of the haloes they inhabit. We explore how the main axes of central galaxies align across different scales, both in three-dimensional and two-dimensional projections. We examine how halo mass influences these alignments and how they vary in the surrounding large-scale environment. To conduct this study, we employ the TNG300 hydrodynamical simulations and compare our results with the spectroscopic data from the SDSS DR18. Three types of alignment were analysed: between stellar and dark matter components, between satellite galaxies and the central galaxy, and between the central galaxy and its host halo. The results show that the alignment increases with halo mass and varies with the environment (clusters, filaments, cluster periphery, and others). However, after controlling for local density, we found that most of the observed trends disappear. The SDSS observations confirm a mass dependence similar to the simulations, although observational biases limit the detection of differences between the different environments.

Protoplanetary disks naturally emerge during protostellar core-collapse. In their early evolutionary stages, infalling material dominates their dynamical evolution. In the context of planet formation, this means that the conditions in young disks are different from the typically considered disks where infall has subsided. High inward velocities are caused by the advection of accreted material which is deficient in angular momentum, rather than being set by viscous spreading, and accretion gives rise to strong velocity fluctuations. Therefore, we aim to investigate when it is possible for the first planetesimals to form and subsequent planet formation to commence. We analyze the disks obtained in numerical 3D nonideal magnetohydrodynamical simulations, which serve as a basis for 1D models representing the conditions during the Class 0/I evolutionary stages. We integrate the 1D models with an adapted version of the TwoPopPy code to investigate the formation of the first planetesimals via the streaming instability. In disks with temperatures such that the snow line is located at ~10 AU and where it is assumed that velocity fluctuations felt by the dust are reduced by a factor of 10 compared to the gas, ${\sim}10^{-3}M_\odot$ of planetesimals may be formed already during the first 100 kyr after disk formation, implying the possible early formation of giant planet cores. The cold-finger effect at the snow line is the dominant driver of planetesimal formation, which occurs in episodes and utilizes solids supplied directly from the envelope, leaving the disk solid reservoir intact. However, if the cold-finger effect is suppressed, early planetesimal formation is limited to cold disks with efficient dust settling whose dust-to-gas ratio is initially enriched to $\epsilon_0\geq 0.03$.

High-contrast observations with JWST can reveal key composition and vertical mixing dependent absorption features in the spectra of directly imaged planets across the 3-5 $\mu$m wavelength range. We present novel coronagraphic images of the HR 8799 and 51 Eri planetary systems using the NIRCam Long Wavelength Bar in an offset "narrow" position. These observations have revealed the four known gas giant planets encircling HR 8799, even at spatial separations challenging for a 6.5 m telescope in the mid-infrared, including the first ever detection of HR 8799 e at 4.6 $\mu$m. The chosen filters constrain the strength of CO, CH4, and CO2 absorption in each planet's photosphere. The planets display a diversity of 3-5 $\mu$m colors that could be due to differences in composition and ultimately be used to trace their formation history. They also show stronger CO2 absorption than expected from solar metallicity models, indicating that they are metal enriched. We detected 51 Eri b at 4.1 $\mu$m and not at longer wavelengths, which, given the planet's temperature, is indicative of out-of-equilibrium carbon chemistry and an enhanced metallicity. Updated orbits fit to the new measurement of 51 Eri b validate previous studies that find a preference for high eccentricities ($e{=}0.57_{-0.09}^{+0.03}$), which likely indicates some dynamical processing in the system's past. These results present an exciting opportunity to model the atmospheres and formation histories of these planets in more detail in the near future, and are complementary to future higher-resolution, continuum-subtracted JWST spectroscopy.

We search for evidence of rotational support by analyzing the thermodynamic profiles of the intracluster medium (ICM) in a sample of nearby, massive galaxy clusters. For each object of the XMM-Newton Cluster Outskirts Project (X-COP) sample, we present axisymmetric models of rotating ICM with composite polytropic distributions, in equilibrium in spherically symmetric dark halos, exploring cases both with and without turbulent support in the ICM. The profile of rotation velocity and the distribution of turbulent velocity dispersion are described with flexible functional forms, consistent with the properties of synthetic clusters formed in cosmological simulations. The models are tuned via a Markov Chain Monte Carlo algorithm to reproduce the radial profiles of the thermodynamic variables as resolved in the XMM-Newton and Planck maps, and to be consistent with the mass distributions estimated either from weak lensing observations (when available) or under the assumption of a "universal" value of the baryon fraction. Our models indicate that there is room for non-negligible rotation in the ICM of massive clusters, with typical peak rotation speed 300 km/s and peak rotation-velocity-to-velocity-dispersion ratio of 0.3. According to our models, the ICM in Abell 2255 can have a rotation speed as high as 500 km/s, corresponding to a rotation-velocity-to-velocity-dispersion ratio of 0.3, at a distance of 100 kpc from the center, where the X-ray emissivity is still high. This makes Abell 2255 a very promising candidate for the presence of rotation in the ICM that could be detected with the currently operating XRISM observatory, as we demonstrate computing and analyzing a mock X-ray spectrum.

In this work, we investigate the properties of hadronic and quark matter that would allow for a first order phase transition between them within neutron stars. To this end, we use a parameterizable Relativistic Mean-Field (RMF) description for the hadronic phase and a Renormalization Group-consistent Nambu-Jona-Lasino (RG-NJL) model for the quark phase. This also enables us to consider sequential phase transitions involving a two-flavor color-superconducting (2SC) and a color-flavor-locked (CFL) phase. We find large ranges for all parameters that facilitate a phase transition, even when constrained by current astrophysical data. We further attempt to filter out stars with a high chance of detectability by mass-radius measurement, i.e., stars with identical mass but different radii, so-called twin stars. However, we find that such configurations are outside the constrained parameter spaces. Instead, most of the mass-radius relations that feature a phase transition appear to be indistinguishable from a purely hadronic description.

Self-interacting dark matter (SIDM) has been proposed to address small-scale challenges faced by the cold dark matter (CDM) paradigm, such as the diverse density profiles observed in dwarf galaxies. In this study, we analyze the kinematics of dwarf galaxies by incorporating the effects of gravothermal core collapse into SIDM models using a semi-analytical subhalo framework. Our analysis covers the stellar kinematics of both classical and ultrafaint dwarf galaxies. The results indicate a bimodal preference for small and large self-interaction cross sections in ultrafaint dwarf galaxies, while in classical dwarfs, larger cross sections progressively decrease the model's statistical support. The combined analysis decisively prefers CDM to SIDM when the self-interaction cross section per unit mass, $\sigma/m$, exceeds $\sim$0.2 cm$^2$/g, if a velocity-independent cross section is assumed. Our study significantly enhances our understanding of dark matter dynamics on small scales.

Redshift-independent distances underpin much of astrophysics, and there exists a plethora of methods to estimate them. However, the extent to which the distances they imply are consistent, while crucial for the integrity of the distance ladder, has been little explored. We construct a statistical framework to assess both internal (between measurements with the same method) and external (between-method) consistency by comparing differences between distances to their quoted statistical uncertainties in the NASA/IPAC Extragalactic Database of Distances (NED-D). 66 of the 76 indicators in NED-D are amenable to a consistency test by having at least two measurements to the same galaxy or at least one measurement to a galaxy also measured by another method. We find that only 12 of these methods produce self-consistent distances across literature determinations, of which 7 are also consistent with distances to the same galaxies measured by all other methods. The most consistent 6 methods (M-stars luminosity, Novae, Masers, Globular Cluster Fundamental Plane, O- and B-type Supergiants and BL Lac Luminosity) also give similar average distances to the mean of all indicators, while the 7th (Proper Motion) underestimates distances relative to the mean by 17.1%. We also investigate consistency of Cepheid distances in the SH0ES 2022 catalogue, finding no evidence for unaccounted-for systematics. Our NED-D results imply that considerable work remains to obtain reliable distances by a multitude of methods, a crucial endeavour for constructing a multiply cross checked and fully robust distance ladder.

We extend the description of equivalent-barotropic equations for exoplanets to the diabatic case -- that is, with explicit heating and/or cooling representation, rather than with a stationary deflection of the bottom bounding surface. In the diabatic case, the equation for potential temperature (or entropy) is directly forced and cannot be decoupled from the equations for momentum and nonlinear pressure, the mass-like variable; and, the isentropic surfaces do not remain coincident with material surfaces. Here the formalism is presented for an atmosphere with the Lamb vertical structure, as the formalism is substantially simplified under the structure. The equations presented set the stage for accurate global simulations which permit small-scale vortices, gravity waves, and fronts observed in current three-dimensional global simulations to be studied in detail.

The Imaging X-ray Polarimetry Explorer (IXPE) observations of accreting X-ray pulsars (XRPs) continue to provide novel insights into the physics and geometry of these sources. We present the first X-ray polarimetric study of the persistent wind-fed XRP 4U 1538-52, based on five IXPE observations totaling 360 ks, conducted in March and October 2024. We detect marginally significant polarization in the combined data set in the full 2--8 keV energy band, with a polarization degree (PD) of 3.0+-1.1% and polarization angle (PA) of -18 degrees. The energy-resolved analysis shows a clear energy dependence of the polarization properties, with a remarkable ~70 degrees switch in PA between low and high energies. Similarly, the pulse phase-resolved spectro-polarimetric analysis reveals different signatures at low and high energies. At low (2--3 keV) energies, the PD ranges between ~2% and ~18%, with the PA varying between -16 and 70 degrees. At higher (4--8 keV) energies, the PD varies between ~3% and ~12%, with a drastically different PA behavior. Fitting the rotating vector model to the pulse phase dependence of the PA at the lower energies, we constrain the geometric configuration of the pulsar. The analysis favors a high spin-axis inclination of >50 which agrees with both previous pulse-phase-dependent spectral fitting of the cyclotron line region and the known high orbital inclination of the binary system. The magnetic obliquity is estimated to be 30 degrees and the spin position angle to be 19 degrees. A sharp switch in PA around 3 keV presents a particular theoretical challenge, as it is not consistent with the right-angle switch that was only seen in one other pulsar Vela X-1.

The massive, hot galaxy cluster PSZ2 G286.98+32.90 (hereafter PLCKG287, z=0.383) hosts a giant radio halo and two prominent radio relics which are signs of a disturbed dynamical state. However, despite optical and radio observations indicate a clear multiple merger, the X-ray emission of the cluster, derived from XMM-Newton observations, shows only moderate disturbance. We present new 200 ks Chandra observations of PLCKG287. We detect a shock front to the NW direction at a distance of ~390 kpc from the X-ray peak, characterized by a Mach number M~1.3, as well as a cold front at a distance of ~300 kpc from the X-ray peak, nested in the same direction of the shock in a typical configuration expected by a merger. We also find evidence for X-ray depressions to the E and W, that could be the signature of feedback from the active galactic nucleus (AGN). The radial profile of the thermodynamic quantities show a temperature and abundance peak in the cluster center, where also the pressure and entropy have a rapid increase. Based on these properties, we argue that PLCKG287 is what remains of a cool core after a heating event. We estimate that both the shock energy and the AGN feedback energy, implied by the analysis of the X-ray cavities, are sufficient to heat the core to the observed temperature of ~17 keV in the central ~160 kpc. We discuss the possible origin of the detected shock by investigating alternative scenarios of merger and AGN outburst, finding that they are both energetically viable. However, no single model seems able to explain all the X-ray features detected in this system. This suggests that the combined action of merger and central AGN feedback is likely necessary to explain the reheated cool core, the large-scale shock and the cold front. The synergy of these two processes may act in shaping the distribution of cool core and non cool core clusters. [Abridged]

[Abridged] Compact groups (CGs) of galaxies have proven to be unique environments for studying galaxy interactions. We propose a detailed analysis of the galaxy evolution to disentangle the relationship between the two first-ranked galaxies in CGs throughout their history as a function of the assembly channels of their hosts. Our study was performed from a semi-analytical point of view, using more than 20000 CGs. We based our analysis on studying the first- (1R) and second-ranked (2R) galaxies in CGs, where the ranking is determined using the galaxy stellar mass. The 1R galaxies have significantly reduced their star-forming capacity over time, reaching a quenching stage and often becoming bulge-dominated or elliptical. Notably, this transformation occurred earlier for 1R galaxies in early-formed CGs (around 5 to 8 Gyrs ago), while those in recently-formed CGs experienced this change more recently (around 2 to 3 Gyrs ago). The analysis of the time evolution of a variant of the Tremaine & Richstone statistics showed that the 1R galaxy in early-formed CGs began to stand out for its dominant properties around 6 Gyrs ago, almost 5 Gyrs earlier than the 1R inhabiting recently formed CGs. We observed that a large majority of the 1R galaxies have experienced at least one major merger event during their life, while we observed this only for a third of 2R galaxies. The 1R galaxies can also display several of these events and most of their last major merger events can be described as the addition of a progenitor that is the second most massive galaxy in their surroundings at the time of the merger. We find that the semi-analytical framework explored in this work describes a scenario where galaxy mergers are the main driving force in shaping the properties of the 1R galaxies in CGs. We note that this scenario is especially intensive when those galaxies inhabit CGs that had formed early on.

The Ultraviolet Near-Infrared Optical Northern Survey (UNIONS) is a "collaboration of collaborations" that is using the Canada-France-Hawai'i Telescope, the Pan-STARRS telescopes, and the Subaru Observatory to obtain $ugriz$ images of a core survey region of 6250 deg$^2$ of the northern sky. The $10\sigma$ point source depth of the data, as measured within a 2-arcsecond diameter aperture, are $[u,g,r,i,z] = [23.7, 24.5, 24.2, 23.8, 23.3]$\ in AB magnitudes. UNIONS is addressing some of the most fundamental questions in astronomy, including the properties of dark matter, the growth of structure in the Universe from the very smallest galaxies to large-scale structure, and the assembly of the Milky Way. It is set to become the major ground-based legacy survey for the northern hemisphere for the next decade and provides an essential northern complement to the static-sky science of the Vera C. Rubin Observatory's Legacy Survey of Space and Time. UNIONS supports the core science mission of the {\it Euclid} space mission by providing the data necessary in the northern hemisphere for the calibration of the wavelength dependence of the {\it Euclid} point-spread function and derivation of photometric redshifts in the North Galactic Cap. This region contains the highest quality sky for {\it Euclid}, with low backgrounds from the zodiacal light, stellar density, extinction, and emission from Galactic cirrus. Here, we describe the UNIONS survey components, science goals, data products, and the current status of the overall program.

Filament eruptions are magnetically driven violent explosions commonly observed on the Sun and late-type stars, sometimes leading to monster coronal mass ejections that directly affect the nearby planets' environments. More than a century of research on solar filaments suggests that the slow evolution of photospheric magnetic fields plays a decisive role in initiating filament eruptions, but the underlying mechanism remains unclear. Using high-resolution observations from the \textit{Chinese H$\alpha$ Solar Explorer}, the \textit{Solar Upper Transition Region Imager}, and the \textit{Solar Dynamics Observatory}, we present direct evidence that a giant solar filament eruption is triggered by a series of minifilament eruptions occurring beneath it. These minifilaments, which are homologous to the giant filament but on a smaller tempo-spatial scale, sequently form and erupt due to extremely weak mutual flux disappearance of opposite-polarity photospheric magnetic fields. Through multi-fold magnetic interactions, these erupting minifilaments act as the last straw to break the force balance of the overlying giant filament and initiate its ultimate eruption. The results unveil a possible novel pathway for small-scale magnetic activities near the stellar surface to initiate spectacular filament eruptions, and provide new insight into the magnetic coupling of filament eruptions across different tempo-spatial scales.

This paper describes the development, design, ground verification, and in-orbit verification, performance measurement, and calibration of the timing system for the X-Ray Imaging and Spectroscopy Mission (XRISM). The scientific goals of the mission require an absolute timing accuracy of 1.0~ms. All components of the timing system were designed and verified to be within the timing error budgets, which were assigned by component to meet the requirements. After the launch of XRISM, the timing capability of the ground-tuned timing system was verified using the millisecond pulsar PSR~B1937+21 during the commissioning period, and the timing jitter of the bus and the ground component were found to be below $15~\mu$s compared to the NICER (Neutron star Interior Composition ExploreR) profile. During the performance verification and calibration period, simultaneous observations of the Crab pulsar by XRISM, NuSTAR (Nuclear Spectroscopic Telescope Array), and NICER were made to measure the absolute timing offset of the system, showing that the arrival time of the main pulse with XRISM was aligned with that of NICER and NuSTAR to within $200~\mu$s. In conclusion, the absolute timing accuracy of the bus and the ground component of the XRISM timing system meets the timing error budget of $500~\mu$s.

Self-gravity is important in protoplanetary disks for planet formation through gravitational instability (GI). We study the cooling effect on GI in a thin two-dimensional protoplanetary disk. By solving the linear perturbation equations in global geometry, we obtain all the normal modes. Faster cooling leads to faster growth rate of GI with lower azimuthal wavenumber $m$. According to the spatial structure of normal modes at different mass-averaged Toomre number $\overline{Q}$ and dimensionless cooling timescale $\beta$, we identify three modes: local, transitional and global. The transitional modes are located in the outer disk while the other two modes in the inner disk. At $\beta\approx 1$ (the resonance of dynamical timescale and thermal timescale) the growth rate changes sharply and the transitional modes dominate. The disk $\alpha$ due to GI is much higher in the transitional modes than in the other two. Our result implies that the transitional modes at $\overline{Q}\approx 1$ and $\beta\approx 1$ can plausibly interpret the substructures and planet/brown dwarf formation in the outer disk.

Observations with the Atacama Large Millimeter/submillimeter Array (ALMA) and the Jansky Very Large Array (JVLA) have revealed many dust rings in protoplanetary disks, often interpreted as dust traps at gas pressure bumps. Previous studies have typically modeled these rings by assuming a single dust species in drift-diffusion equilibrium, neglecting dust size evolution resulting from coagulation and fragmentation. In this work, we perform numerical simulations that incorporate both dust-gas dynamics (drift and diffusion) and dust size evolution. Our results show that the radial distributions of different dust species are nearly identical, as dust growth dominates over drift and diffusion. Building on this finding, we develop a comprehensive, self-consistent analytical theory that describes the dust ring structure while explicitly accounting for size evolution effects. Our model provides a unified framework for interpreting multi-wavelength observations by linking the physical dust distribution to the observed ring properties, thus laying the foundation for future observational modeling.

We present polarization profiles of 23 pulsars exhibiting interpulse (IP) emissions using the Five-hundred-meter Aperture Spherical Radio Telescope (FAST). By applying the rotating vector model (RVM), we measured the inclination and impact angles for 16 pulsars, allowing us to investigate their beam geometries. Our analysis indicates that for 5 pulsars, the IP emissions likely originate from the same magnetic pole as the main pulse (MP), whereas for the remaining 11 pulsars, the IP and MP emissions originate from opposite magnetic poles. For the 7 pulsars that do not conform to the RVM, we are unable to determine whether the IP emissions originate from the same or opposite magnetic pole as the MPs. By analyzing the beam configurations of pulsars, we found that the emission within the beam is not fully active across both longitude and latitude. Filling factors ranging from 0.27 to 0.99 were obtained for pulsars with opposite pole IP emissions, suggesting an underestimation of emission height when applying the geometrical method. The emissions for MPs and IPs occur at different heights in the pulsar magnetosphere, with the difference in emission height ranging from tens to thousands of kilometers. We also found that some pulsars have wide emission beams, indicating that radio emissions may occur in regions of high altitude within the pulsar magnetosphere.

Solar flare above-the-loop-top (ALT) regions are vital for understanding solar eruptions and fundamental processes in plasma physics. Recent advances in 3D MHD simulations have revealed unprecedented details on turbulent flows and MHD instabilities in flare ALT regions. Here, for the first time, we examine the observable anisotropic properties of turbulent flows in ALT by applying a flow-tracking algorithm on narrow-band Extreme Ultraviolet (EUV) images that are observed from the face-on viewing perspective. First, the results quantitatively confirm the previous observation that vertical motions dominate and that the anisotropic flows are widely distributed in the entire ALT region with the contribution from both upflows and downflows. Second, the anisotropy shows height-dependent features, with the most substantial anisotropy appearing at a certain middle height in ALT, which agrees well with the MHD modeling results where turbulent flows are caused by Rayleigh-Taylor-type instabilities in the ALT region. Finally, our finding suggests that supra-arcade downflows (SADs), the most prominently visible dynamical structures in ALT regions, are only one aspect of turbulent flows. Among these turbulent flows, we also report the anti-sunward-moving underdense flows that might develop due to MHD instabilities, as suggested by previous three-dimensional flare models. Our results indicate that the entire flare fan displays group behavior of turbulent flows where the observational bright spikes and relatively dark SADs exhibit similar anisotropic characteristics.

Solar type-II radio bursts are coherent plasma emissions arising from magnetohydrodynamic shocks produced by either coronal mass ejections (CMEs) or flares. Type-II bursts sometimes show split-band emissions in the dynamic spectrum. When these split-band emissions come from regions just upstream and downstream of the shock, type-II band-splitting can be used as an important tool for estimating magnetic fields at the shock front. Earlier studies have shown that only $\sim$20\% of the type-IIs show morphologically similar split-bands. Imaging studies can unambiguously identify such instances, though they remain very rare. Here we suggest a useful approach to augment dynamic spectra-based studies by also examining the Gradient Dynamic Spectra (GraDS) of type-II emission. We also verified the conclusions of this approach against those from an imaging study.

We present results from a search for radio afterglows of compact object mergers conducted with the Australian SKA Pathfinder. We used data from four epochs of the Rapid ASKAP Continuum Survey to search compact binary merger localization regions observed during the LIGO/Virgo O2, and O3 observing runs. Our investigation focused on eleven events (published in the GWTC-1, GWTC-2, and GWTC-3 catalogues of gravitational-wave events) with 90\% posterior localisations smaller than $150\,°^2$ and $\ge$99\% probabilities of being of astrophysical origin, to identify potential radio afterglow-like transients up to $\lesssim$1500 days post-merger. We identified candidate afterglow-type variable sources in the 90\% localisation for events -- GW190503, GW200202 and GW200208, which were ruled out as unlikely to be related to the corresponding GW event on further analysis. Since we find no likely candidate counterparts, we constrain the inclination angle and the circum-merger density at isotropic equivalent energies ranging from $2\times10^{51} -1\times10^{54}\rm \:erg$. These constraints are based on the assumption that the electron energy distribution in the associated jets follows a power-law index of $ p = 2.2$, with 1% of the shock energy in the magnetic field ($ \epsilon_B = 0.01$) and 10% in the electrons ($\epsilon_e = 0.1$). We discuss the detectability of late-time afterglows as a function of merger distance and inclination angles with millijansky surveys.

The formation mechanisms of open cluster (OCs) groups remain unclear due to limited sample sizes and data precision. Recent advancements in Gaia astrometric data provide an unprecedented opportunity to study OC groups in greater detail. This study aims to extend the sample of OC groups and investigate their formation and evolution mechanisms, with a focus on the role of stellar feedback in triggering star formation. We identify four new OC groups based on Gaia data, whose member OCs are spatially proximate and kinematically coherent. Their age spreads are consistent with the timescale of continuous star formation, suggesting that their member OCs formed sequentially from the same molecular cloud. N-body simulation results further reveal that these groups will gradually disperse, evolving into independent OCs. By analyzing the correlation between OC ages and their separation from potential SN explosion sites, we predict SN explosion regions around the birthplaces of OC groups. The strong correlation between OC ages and predicted SN explosion sites supports a supernova-triggered star formation scenario. Additionally, we trace pulsar (PSR) orbits to examine their association with these regions. We detected three PSRs near Group 1 and 26 PSRs near Group 2, whose birthplaces align with the predicted SN explosions regions. The presence of PSRs associated with OC groups provides additional observational evidence for SN explosions in this region, further supporting a supernova-triggered star formation scenario for G1 and G2. We propose that multiple SN explosions in a short period triggered the formation of Group 1 and Group 2, reinforcing the hierarchical star formation model. These results highlight the multi-scale interactions driving star and OC formation and provide new insights into the role of stellar feedback in shaping OC groups.

Type Ia supernovae (SNe Ia) serve as crucial cosmological distance indicators due to their empirical consistency in peak luminosity and characteristic light-curve decline rates. These properties facilitate them to be standardized candles for the determination of the Hubble constant ($H_0$) within late-time universe cosmology. Nevertheless, a statistically significant difference persists between $H_0$ values derived from early and late-time measurements, a phenomenon known as the Hubble tension. Furthermore, recent observations have identified a subset of over-luminous SNe Ia, characterized by peak luminosities exceeding the nominal range and faster decline rates. These discoveries raise questions regarding the reliability of SNe Ia as standard candles in measuring cosmological distances. In this article, we present the Bayesian analysis of eight over-luminous SNe Ia and show that they yield a lower $H_0$ estimates, exhibiting closer concordance with $H_0$ estimates derived from early-universe data. This investigation potentially represent a step toward addressing the Hubble tension.

The variability of Young Stellar Objects (YSOs) is a crucial tool for understanding the mechanisms driving flux changes. In this study, we present an infrared variability analysis of a large sample of over 20,000 candidate YSOs, using data from the ALLWISE and NEOWISE surveys, which span around a decade with a 6-month cadence. We applied Lomb-Scargle Periodogram (LSP) analysis and linear fitting to the light curves, classifying them into distinct categories: {\it Secular} ({\it Linear}, {\it Curved}, and {\it Periodic}) and {\it Stochastic} ({\it Burst}, {\it Drop}, and {\it Irregular}). Our findings show that 5,467 (26.2$\pm$0.3\%) of the sources exhibit variability, with most (19.7$\pm$0.3\%) showing {\it Irregular} variations, followed by {\it Curved} and {\it Periodic} variations. In addition, 235 sources of {\it Bursts} and 122 {\it Drop} sources were identified. Variability is more pronounced in Class I sources with a higher fraction of variables (36.3$\pm$0.6\%) compared to Class II (22.1$\pm$0.4\%) and Class III (22.5$\pm$1.0\%) sources. The color (W1 $-$ W2) versus magnitude analysis (W2) using linear fitting shows that the trend ``redder-when-brighter" (RWB) is more prevalent (85.4$\pm$0.5\%) among YSOs. In contrast, the trend ``bluer-when-brighter" (BWB) is more common in younger sources compared to more evolved ones, having a BWB fraction of 29.0$\pm$1.1\% for Class I to 4.0$\pm$0.9\% for Class III.

We present photometric and spectroscopic observations of SN 2023ixf covering from day one to 442 days after explosion. SN 2023ixf reached a peak $V$-band absolute magnitude of $-18.2 \pm 0.07$, and light curves show that it is in the fast-decliner (IIL) subclass with a relatively short ``plateau'' phase (fewer than $\sim 70$ days). Early-time spectra of SN 2023ixf exhibit strong, very narrow emission lines from ionized circumstellar matter (CSM), possibly indicating a Type IIn classification. But these flash/shock-ionization emission features faded after the first week and the spectrum evolved in a manner similar to that of typical Type II SNe, unlike the case of most genuine SNe~IIn in which the ejecta interact with CSM for an extended period of time and develop intermediate-width emission lines. We compare observed spectra of SN 2023ixf with various model spectra to understand the physics behind SN 2023ixf. Our nebular spectra (between 200-400 d) match best with the model spectra from a 15 $\rm M_{\odot}$ progenitor which experienced enhanced mass loss a few years before explosion. A last-stage mass-loss rate of $\dot{M} = 0.01 \rm M_{\odot} yr^{-1}$ from the r1w6 model matches best with the early-time spectra, higher than $\dot{M} \approx 2.4 \times 10^{-3} \rm M_{\odot} yr^{-1}$ derived from the ionized H${\alpha}$ luminosity at 1.58 d. We also use SN 2023ixf as a distance indicator and fit the light curves to derive the Hubble constant by adding SN 2023ixf to the existing sample; we obtain H$_{0}=73.1^{+3.68}_{-3.50}$ km s$^{-1}$ Mpc$^{-1}$, consistent with the results from SNe~Ia and many other independent methods.

In this paper, we investigate the heliospheric modulation of cosmic rays in interplanetary space, focusing on their propagation times and energy losses over the solar cycle. To perform the calculations, we employed a data-driven model based on the stochastic method. Our model was calibrated using time-resolved and energy-resolved data from several missions including AMS-02, PAMELA, EPHIN/SOHO, BESS, and data from Voyager-1. This approach allows us to calculate probability density functions for the propagation time and energy losses of cosmic protons and antiprotons in the heliosphere. Furthermore, we explore the temporal evolution of these probabilities spanning from 1993 to 2018, covering a full 22-year cycle of magnetic polarity, which includes two solar minima and two magnetic reversals. Our calculations were carried out for cosmic protons and antiprotons, enabling us to investigate the role of charge-sign dependent effects in cosmic ray transport. These findings provide valuable insights into the physical processes of cosmic-ray propagation in the heliosphere and contribute to a deeper understanding of the solar modulation phenomenon.

Light sterile neutrinos, $\nu_s$, are often introduced to explain an anomalous deficit in the electron antineutrino flux from nuclear reactors. If they exist, sterile neutrinos would also be produced in collapsing massive stars through the active-sterile neutrino oscillation. In order to investigate the impacts of sterile neutrinos on supernova dynamics, we perform two-dimensional neutrino-radiation hydrodynamic simulations of stellar core-collapse coupled with the active-sterile oscillation. As the initial condition of our simulations, we adopt a blue supergiant model that is tuned to reproduce observational features of the SN 1987A progenitor to compare our models with observations of the event. It is found that the active-sterile oscillation reduces the $\nu_{e}$ and $\bar{\nu}_e$ fluxes and decreases the explosion energy. We also find that, if the mixing angle $\theta$ and the mass difference $\delta m_\mathrm{s}^2$ between $\nu_e$ and $\nu_s$ are large enough, the star fails to explode. This suggests that these mixing parameters relevant to sterile neutrinos could be constrained by supernova explodability, though other uncertainties in supernova theory need to be addressed to refine them. In addition, we predict neutrino signals from a nearby supernova event and find that the neutrino event number can significantly decrease because the $\nu_e$ and $\bar{\nu}_e$ fluxes are reduced. In particular, DUNE observations of $\nu_e$ will be useful to search for a signature of sterile neutrinos with a tiny mixing angle because a smaller mixing angle leads to a larger effect on the $\nu_e$ flux.

This is the second paper in a series that utilize IFS from MaNGA, NUV imaging from Swift/UVOT and NIR imaging from 2MASS to study dust attenuation properties on kpc scales in nearby galaxies. We apply the method developed in Paper I (Zhou et al. 2023) to the updated SWiM_v4.2 catalog, and measure the optical attenuation curve and the attenuation in three NUV bands for 2487 spaxels selected from 91 galaxies with S/N>20 and $A_V$>0.25. We classify all spaxels into two subsets: star-forming (SF) regions and non-SF regions. We explore the correlations of optical opacity ($A_V$) and the optical and NUV slopes of attenuation curves ($A_B/A_V$ and $A_{w2}/A_{w1}$) with a broad range of stellar and emission-line properties, including specific surface brightness of H$\alpha$ emission, stellar age, stellar and gas-phase metallicity, and diagnostics of recent star formation history. When comparing SF and non-SF regions, we find that $A_V$ and $A_B/A_V$ exhibit similar correlations with all the stellar population and emission-line properties considered, while the NUV slopes in SF regions tend to be flatter than those in non-SF regions. The NUV slope $A_{w2}/A_{w1}$ exhibits an anti-correlation with specific surface brightness of H$\alpha$ emission, a trend that is primarily driven by the positive correlation between $A_{w2}/A_{w1}$ and $\Sigma_\ast$. The NUV slope flattens in SF regions that contain young stellar populations and have experienced recent star formation, but it shows no obvious dependence on stellar or gas-phase metallicity. The spatially resolved dust attenuation properties exhibit no clear correlations with the inclination of host galaxies or the galactocentric distance of the regions. This finding reinforces the conclusion from Paper I that dust attenuation is primarily regulated by local processes on kpc scales or smaller, rather than by global processes at galactic scales.

The analysis of variations in the emission intensity of the pulsar B0950+08 from 2014 to 2022 with scales from minutes to years was carried out. The observations were obtained in a round-the-clock daily survey conducted on the Large Phased Array (LPA) radio telescope. The high variability of emission is shown not only from pulse to pulse, but also at scales greater than 3 min. The average value of the estimated amplitude of these variations in 3.2 minutes is 25~Jy, the modulation index is 1. The average relative amplitude of the interpulse (IP) is $2.00 \pm 0.28\%$ of the main pulse. In individual pulses, the amplitude of the interpulse may exceed the amplitude of the main pulse (MP), but this is a rare event. Emission is observed in almost the entire period of the pulsar. For the first time, the relative amplitude of emission between the main pulse and the interpulse (emission bridge) was measured. When averaging about 10 hours, it varies from $0.8\%$ to $1.31\%$ with an average value of $1.04 \pm 0.28\%$. A high correlation was found between MP and IP amplitude variations both when averaging profiles over 3.2 minutes and when averaging over years. This correlation is due to refractive interstellar scintillation. The frequency scale of IP diffraction interstellar scintillation was measured for the first time and it was shown that the spectral forms for IP and MP are well correlated and have the same frequency scale. There are strong variations in the frequency scale of scintillation $f_{dif}$ from session to session (time interval from one day) on scales of 200-800 kHz. The refractive scale of scintillation for 1-2 days has been determined. A modulation of emission with a characteristic scale of about 130 days was detected, which, apparently, is also associated with refractive scintillation.

We observed the new Long Period Comet C/2024 E1 (Wierzchos), inbound at 7 au from the Sun, using the NIRSpec integral field unit on JWST. The spectrum shows absorption features due to water ice in the coma and evidence for CO$_2$ driven activity, with a production rate of $Q(CO_2) = 2.546 \pm 0.019 \times 10^{25}$ molecules s$^{-1}$, and no emission features of water or CO. The latter is surprising, given that CO is more volatile than CO$_2$, and suggests that this comet may have lost its near-surface CO during its early evolution, before implantation in the Oort cloud.

The Gaia DR3 has provided a large sample of more than 6.6 million quasar candidates with high completeness but low purity. Previous work on the CatNorth quasar candidate catalog has shown that including external multiband data and applying machine-learning methods can efficiently purify the original Gaia DR3 quasar candidate catalog and improve the redshift estimates. In this paper, we extend the Gaia DR3 quasar candidate selection to the southern hemisphere using data from SkyMappper, CatWISE, and VISTA surveys. We train an XGBoost classifier on a unified set of high-confidence stars and spectroscopically confirmed quasars and galaxies. For sources with available Gaia BP/RP spectra, spectroscopic redshifts are derived using a pre-trained convolutional neural network (RegNet). We also train an ensemble photometric redshift estimation model based on XGBoost, TabNet, and FT-Transformer, achieving an RMSE of 0.2256 and a normalized median absolute deviation of 0.0187 on the validation set. By merging CatSouth with the previously published CatNorth catalog, we construct the unified all-sky CatGlobe catalog with nearly 1.9 million sources at $G<21$, providing a comprehensive and high-purity quasar candidate sample for future spectroscopic and cosmological investigations.

The NASA Lucy mission is scheduled to fly-by the main belt asteroid (52246) Donaldjohanson on April 20, 2025. Donaldjohanson (DJ hereafter) is a member of the primitive (C-type class) Erigone collisional asteroid family located in the inner main belt in proximity of the source regions of asteroid (101955)~Bennu and (162173)~Ryugu, visited respectively by OSIRIS-REx and Hayabusa2 missions. In this paper we provide an updated model for the Erigone family age, and discuss DJ evolution resulting from non-gravitational forces (namely Yarkovsky and YORP), as well as its collisional evolution. We conclude the best-fit family age to be $\sim 155$~Myr, and that, on such timescales, both Yarkovsky and YORP effects may have affected the orbit and spin properties of DJ. Furthermore, we discuss how the NASA Lucy mission could provide independent insights on such processes, namely by constraining DJ shape, surface geology and cratering history.

Characterizing the temporal variability of astrophysical sources is key to understanding the underlying physical processes driving their emissions. This work introduces a gammapy_SyLC, a Python package that offers tools to simulate and fit time-domain data, with a focus on Active Galactic Nuclei (AGN) variability. The package was developed taking into account possible interactions with gammapy but does not directly depend on it. gammapy_SyLC incorporates optimized implementations of the Timmer & Koenig and Emmanoulopoulos algorithms for light curve simulation, capable of generating synthetic lightcurves from specified PSDs and amplitude distribution models. It also provides functionalities for PSD fitting, histogram-based PDF interpolation, and Monte Carlo-based parameter estimation, making it a full-stack tool for investigating variable phenomena and specifically the long-term behavior of AGNs. To showcase its capabilities, the package was applied to gamma-ray light curves from the Fermi Large Area Telescope repository, reconstructing PSDs and PDFs and constraining variability models for observed sources.

Stars in galactic centers are occasionally scattered so close to the central supermassive black hole that they are completely disrupted by tidal forces, initiating a transient accretion event. The aftermath of such a tidal disruption event (TDE) produces a bright-and-blue accretion flow which is known to persist for at least a decade (observationally) and can in principle produce ionizing radiation for hundreds of years. Tidal disruption events are known (observationally) to be overrepresented in galaxies which show extended emission line regions (EELRs), with no pre-TDE classical AGN activity, and to produce transient ``coronal lines'', such as [FeX] and [FeXIV]. Using coupled CLOUDY-TDE disk simulations we show that tidal disruption event disks produce a sufficient ionizing radiation flux over their lifetimes to power both EELR of radial extents of $r \sim 10^4$ light years, and coronal lines. EELRs are produced when the ionizing radiation interacts with low density $n_H \sim 10^1 - 10^3 \, {\rm cm}^{-3}$ clouds on galactic scales, while coronal lines are produced by high density $n_H \sim 10^6 - 10^8 \, {\rm cm}^{-3}$ clouds near the galactic center. High density gas in galactic centers will also result in the rapid switching on of narrow line features in post-TDE galaxies, and also various high-ionization lines which may be observed throughout the infrared with JWST. Galaxies with a higher intrinsic rate of tidal disruption events will be more likely to show macroscopic EELRs, which can be traced to originate from the previous tidal disruption event in that galaxy, which naturally explains why TDEs are more likely to be discovered in galaxies with EELRs. We further argue that a non-negligible fraction of so-called optically selected ``AGN'' are tidal disruption events.

We present the detection of a previously unknown giant molecular cloud (GMC) located at the midpoint of the Galactic Bar Dust Lanes (M4.7--0.8), using spectral line observations taken with the Green Bank Telescope (GBT). This $\sim$60 pc long GMC is associated with accreting material that is transitioning from the quieter Galactic disk environment to the more extreme central molecular zone (CMZ) environment. Our 24 GHz single-dish radio observations targeted the NH$_3$ (1,1)$-$(4,4) and HC$_5$N (9$-$8), known dense gas tracers. The observations reveal the main features of the GMC, which we have dubbed the `Nexus' and `Filament', covering a 0$.\!\!^\circ$5$\times$0$.\!\!^\circ$25 area at 31$''$ angular resolution. In this publication we investigate the gas kinematics within the observed region and compare the distribution of molecular emission to previous infrared surveys to better understand the dust component. The observed gas tracers show centrally condensed cores corresponding to the positions of high dust column densities and low dust temperatures. We report the detection of a previously unknown NH$_3$ (3,3) maser, along with a 70$\mu$m source association, which supports the identification of this region as being actively star-forming. Gas emission in this region shows broad linewidths, comparable to values seen in CMZ clouds. The overall description of this cloud that we present is that of a highly dynamic region comprising dense gas and dust. This encapsulates a wide range of features associated with star formation, in addition to material transport related to the CMZ.

Artificial intelligence (AI) is transforming scientific research, with deep learning methods playing a central role in data analysis, simulations, and signal detection across particle, nuclear, and astroparticle physics. Within the JENA communities-ECFA, NuPECC, and APPEC-and as part of the EuCAIF initiative, AI integration is advancing steadily. However, broader adoption remains constrained by challenges such as limited computational resources, a lack of expertise, and difficulties in transitioning from research and development (R&D) to production. This white paper provides a strategic roadmap, informed by a community survey, to address these barriers. It outlines critical infrastructure requirements, prioritizes training initiatives, and proposes funding strategies to scale AI capabilities across fundamental physics over the next five years.

The analysis of state-of-the-art cosmological surveys like the Dark Energy Spectroscopic Instrument (DESI) survey requires high-resolution, large-volume simulations. However, the computational cost of hydrodynamical simulations at these scales is prohibitive. Instead, dark matter (DM)-only simulations are used, with galaxies populated a posteriori, typically via halo occupation distribution (HOD) models. While effective, HOD models are statistical in nature and lack full physical motivation. In this work, we explore using neural networks (NNs) to learn the complex, physically motivated relationships between DM haloes and galaxy properties. Trained on small-volume, high-resolution hydrodynamical simulations, our NN predicts galaxy properties in a larger DM-only simulation and determines which galaxies should be classified as luminous red galaxies (LRGs). Comparing the original LRG sample to the one generated by our NN, we find that, while the subhalo mass distributions are similar, our NN selects fewer low-mass subhaloes as LRG hosts, possibly due to the absence of baryonic feedback effects in DM-only simulations. This feedback could brighten or redden galaxies, altering their classification. Finally, we generate a new LRG sample by fitting an HOD model to the NN-generated LRG sample. We verify that both the HOD- and NN-generated samples preserve a set of bias parameter relations, which assume that the higher-order parameters, $b_{s2}$ and $b_{3\rm{nl}}$, are determined by the linear bias parameter $b_{1}$. These relations are commonly used to simplify clustering analyses.

The XMM-Newton observatory has accumulated a vast archive of over 17,000 X-ray observations over the last 25 years. However, the standard data processing pipelines may fail to detect certain types of transient X-ray sources due to their short-lived or dim nature. Identifying these transient sources is important for understanding the full range of temporal X-ray behaviour, as well as understanding the types of sources that could be routinely detected by future missions such as Athena. This work aims to reprocess XMM-Newton archival observations using newly developed dedicated software in order to identify neglected and missed transient X-ray sources that were not detected by the existing pipeline. We use a new approach that builds upon previous methodologies, by transforming event lists into data cubes, which are then searched for transient variability in short time windows. Our method enhances the detection capabilities in the Poisson regime by accounting for the statistical properties of sparse count rates, and allowing for transient search in previously discarded periods of high background activity. Our reprocessing efforts identified 32,247 variable sources at the 3-sigma level and 4,083 sources at the 5-sigma level in 12,926 XMM archival observations. We highlight four noteworthy sources: A candidate quasi-periodic eruption (QPE), a new magnetar candidate, a previously undetected Galactic hard X-ray burst and a possible X-ray counterpart to a Galactic radio pulsar. Our method demonstrates a new, fast, and effective way to process event list data from XMM-Newton, which is efficient in finding rapid outburst-like or eclipsing behaviour. This technique can be adapted for use with future telescopes, such as Athena, and can be generalised to other photon counting instruments operating in the low-count Poisson regime.

In this work we present results of the first in-depth analysis of extra-tidal mock stars of Milky Way globular clusters recently generated by Grondin et al. (2024). Particularly, we selected a sample of globular clusters with a general consensus of being formed in the bulge or in the disk of the Milky Way. From the catalog we estimated the width and the dispersion in the z-component of the angular momentum and in the line-of-sight and tangential velocities of their tidal tails, and compared the results with those predicted by cosmological simulations of Malhan et al. (2021) and observations. We found that the resulting values of these four quantities are not in agreement with an in-situ formation of the associated globular clusters. On average, the resulting widths agree with an in-situ origin, while the dispersion in the z-component of the angular momentum, and the dispersion in the line-of-sight and in the tangential velocities fail in matching this formation scenario. The four quantities derived for globular clusters formed in the bulge or in the disk show similar correlations with the stream length, namely: the width and the dispersion in the z-component of the angular momentum increase with the stream length, while the bulk of dispersion values in the line-of-sight and in the tangential velocities is around 12 km s$^{\rm -1}$ along the mock stream.

The tilts of bipolar magnetic regions are believed to be caused by the action of Coriolis force on rising magnetic flux tubes. Here we analysed the combined Greenwich and Debrecen observatories sunspot-group data during the period 1874-2017 and the tilt angles of sunspot groups measured at Mt. Wilson Observatory during the period 1917-1986 and Debrecen Observatory during the period 1994-2013. We find that there exists about 8-solar cycle (Gleissberg cycle) trend in the long-term variation of the slope of Joy's law (increase of tilt angle with latitude). There exists a reasonably significant correlation between the slope/coefficient of Joy's law and the slope (namely, residual covariance) of the linear relationship between the rotation residuals and meridional motions of sunspot groups in the northern hemisphere and also in the southern hemisphere during Solar Cycles 16-21. We also find that there exists a good correlation between north--south difference (asymmetry) in the coefficient of Joy's law and that in the residual covariance. We consider the residual covariance represents tentatively the coefficient of angular momentum transport. These results suggest that there exists a relationship between the surface/subsurface poleward/equatorward angular momentum transport and the Joy's law. There is a suggestion of the strength of the Joy's law depends on the strength of the poleward angular momentum transport.

We develop a theoretical framework to provide observational constraints on the early Universe galaxy-halo connection by combining measurements of the UV luminosity function (UVLF) and galaxy clustering via the 2-point correlation function (2PCF). We implemented this framework in the FRESCO and CONGRESS JWST NIRCam/grism surveys by measuring the 2PCF of spectroscopically selected samples of H$\alpha$ and [OIII] emitters at $3.8<z<9$ in 124 arcmin$^2$ in GOODS-N and GOODS-S. By fitting the 2PCF and UVLF at $3.8<z<9$ we inferred that the H$\alpha$ and [OIII] samples at $\langle z \rangle \sim4.3, 5.4$ and $7.3$ reside in halos of masses of log$(M_{\rm h}/$M$_{\odot}) = 11.5$, $11.2$, $11.0$ respectively, while their galaxy bias increases with redshift with values of $b_{\rm g} = 4.0$, $5.0$, $7.6$. These halos do not represent extreme overdense environments at these epochs. We constrain the instantaneous star formation efficiency (SFE), defined as the ratio of the star formation rate over the baryonic accretion rate as a function of halo mass. The SFE rises with halo mass, peaks at $\sim20\%$ at $M_{\rm h} \sim 3 \times 10^{11}\, M_{\odot}$, and declines at higher halo masses. The SFE-$M_{\rm h}$ shows only a mild evolution with redshift with tentative indications that low mass halos decrease but the high mass halos increase in efficiency with redshift. The scatter in the $M_{\rm UV}-M_{\rm h}$ relation, quantified by $\sigma_{\rm UV}$, implies stochasticity in the UV luminosities of $\sim 0.7$ mag, relatively constant with z. Extrapolating our model to $z>9$ shows that a constant SFE-$M_{\rm h}$ fixed at $z=8$ cannot reproduce the observed UVLF and neither high maximum SFE nor high stochasticity alone can explain the high abundances of luminous galaxies seen by JWST. Extending the analysis of the UVLF and 2PCF to $z>9$ measured from wider surveys will be crucial in breaking degeneracies.

High-Energy Physics (HEP) and Gravitational Wave (GW) communities serve different scientific purposes. However, their methodologies might potentially offer mutual enrichment through common software developments. A suite of libraries is currently being prototyped and made available at this https URL, extending at no cost the CERN ROOT data analysis framework toward advanced signal processing. We will also present a performance benchmark comparing the FFTW and KFR library performances.

We report measurements of the linear polarisation degree (PD) and angle (PA) for hard X-ray emission from the Crab pulsar and wind nebula. Measurements were made with the XL-Calibur ($\sim$15-80 keV) balloon-borne Compton-scattering polarimeter in July 2024. The polarisation parameters are determined using a Bayesian analysis of Stokes parameters obtained from X-ray scattering angles. Well-constrained ($\sim$8.5$\sigma$) results are obtained for the polarisation of the $\sim$19-64 keV signal integrated over all pulsar phases: PD=(25.1$\pm$2.9)% and PA=(129.8$\pm$3.2)$^\circ$. In the off-pulse (nebula-dominated) phase range, the PD is constrained at $\sim$4.5$\sigma$ and is compatible with the phase-integrated result. The PA of the nebular hard X-ray emission aligns with that measured by IXPE in the 2-8 keV band for the toroidal inner region of the pulsar wind nebula, where the hard X-rays predominantly originate. For the main pulsar peak, PD=(32.8$^{+18.2}_{-28.5}$)% and PA=(156.0 $\pm$ 21.7)$^\circ$, while for the second peak (inter-pulse), PD=(0.0$^{+33.6}_{-0.0}$)% and PA=(154.5 $\pm$ 34.5)$^\circ$. A low level of polarisation in the pulsar peaks likely does not favour emission originating from the inner regions of the pulsar magnetosphere. Discriminating between Crab pulsar emission models will require deeper observations, e.g. with a satellite-borne hard X-ray polarimeter.

In this study, we conduct a comparative analysis of the properties of Blazhko and non-Blazhko RRab stars. We identified 1054 non-Blazhko and 785 Blazhko RRab stars in the photometric data observed by K2 mission, which, combined with those 37 stars observed in the original Kepler field, constituted our study sample. Using the Fourier Decomposition method, we calculated the pulsation parameters, including phase differences and amplitude ratios, for these RRab stars, revealing significant discrepancies in the pulsation parameters between Blazhko and non-Blazhko RRab stars. However, distinguishing between Blazhko and Non-Blazhko RRab stars based on Fourier parameters remains challenging due to the significant overlap in their distributions. By cross-matching our sample with the LRS of LAMOST DR12, we identified 147 Blazhko and 111 non-Blazhko RRab stars, which exhibit similar metallicity distributions. Furthermore, cross-matching with Gaia DR3 data yielded 766 Blazhko and 950 non-Blazhko RRab stars, showing differences in color indices but not in absolute magnitudes. Our findings suggested the Blazhko effect is linked to pulsation parameters and colors, rather than metallicities or absolute magnitude.

We combine parallax distances to nearby O stars with parsec-scale resolution three-dimensional dust maps of the local region of the Milky Way (within 1.25 kpc of the Sun) to simulate the transfer of Lyman continuum photons through the interstellar medium. Assuming a fixed gas-to-dust ratio, we determine the density of ionized gas, electron temperature, and H$\alpha$ emissivity throughout the local Milky Way. There is good morphological agreement between the predicted and observed H$\alpha$ all-sky map of the Wisconsin H$\alpha$ Mapper. We find that our simulation underproduces the observed H$\alpha$ emission while overestimating the sizes of HII regions, and we discuss ways in which agreement between simulations and observations may be improved. Of the total ionizing luminosity of $5.84 \times 10^{50}~{\rm photons~s^{-1}}$, 15% is absorbed by dust, 64% ionizes "classical'' HII regions, 11% ionizes the diffuse warm ionized medium, and 10% escapes the simulation volume. We find that 18% of the high altitude ($|b| > 30^{\circ}$) H$\alpha$ arises from dust-scattered rather than direct emission. These initial results provide an impressive validation of the three-dimensional dust maps and O-star parallaxes, opening a new frontier for studying the ionized ISM's structure and energetics in three dimensions.

The Magellanic Stream (MS) is a large tail of neutral and ionized gas originating from tidal and hydrodynamical interactions between the Magellanic Clouds as they orbit the Milky Way (MW). It carries a significant gas reservoir that could impact the future evolution of the MW. Despite its importance, no direct observational constraints on the Stream's distance have been previously published. In this study, we analyze Very Large Telescope/Ultraviolet and Visual Echelle Spectrograph (VLT/UVES) spectra of five blue horizontal branch (BHB) stars in the MW halo located at distances ranging from 13 to 56 kpc near two regions of the Stream, with the aim of detecting Ca II and Na I absorption. No Ca II or Na I absorption is detected at Stream velocities in any of the individual spectra, or in higher signal-to-noise stacks of the spectra. The resulting limits on the Ca II absorption are significantly lower than the Ca II columns measured in the Stream along extragalactic directions. These non-detections establish a firm lower distance limit of 20 kpc for the two regions of the Stream studied. Ca II non-detections in the most distant stars yield tentative lower distance limits of 42 kpc and 55 kpc for the two regions, but deeper spectra are needed to confirm this. Our results provide the first observational constraints on the gaseous Stream's distance.

With the state-of-the-art Gaia astrometry, the number of confirmed white dwarfs has reached a few hundred thousand. We have reached the era where small features in the white dwarf luminosity function (WDLF) of the solar neighbourhood can be resolved. We demonstrate how to apply Markov chain Monte Carlo sampling on a set of pre-computed partial-WDLFs to derive the star formation history of their progenitor stellar populations. We compare the results against many well-accepted and established works using various types of stars, including white dwarfs, main sequence stars, sub-giants and the entire stellar population. We find convincing agreements among most of the methods, particularly at the intermediate age of 0.1-9 Gyr.

We demonstrate by three-dimensional hydrodynamical simulations of energy deposition into the envelope of a red supergiant (RSG) model the inflation of a Rayleigh-Taylor unstable envelope that forms a compact clumpy circumstellar material (CSM). Our simulations mimic vigorous core activity years to months before a core-collapse supernova (CCSN) explosion that deposits energy to the outer envelope. The fierce core nuclear activity in the pre-CCSN explosion phase might excite waves that propagate to the envelope. The wave energy is dissipated where envelope convection cannot carry the energy. We deposit this energy into a shell in the outer envelope with a power of L(wave)=2.6e6Lo or L(wave)=5.2e5Lo for 0.32 year. The energy-deposition shell expands while its pressure is higher than its surroundings, but its density is lower. Therefore, this expansion is Rayleigh-Taylor unstable and develops instability fingers. Most of the inflated envelope does not reach the escape velocity in the year of simulation but forms a compact and clumpy CSM. The high density of the inflated envelope implies that if a companion is present in that zone, it will accrete mass at a very high rate and power a pre-explosion outburst.

Star formation is a fundamental, yet poorly understood, process of the Universe. It is important to study how star formation occurs in different galactic environments. Thus, here, in the first of a series of papers, we introduce the Low-Metallicity Star Formation (LZ-STAR) survey of the Sh2-284 (hereafter S284) region, which, at $Z\sim 0.3-0.5Z_\odot$, is one of the lowest-metallicity star-forming regions of our Galaxy. LZ-STAR is a multi-facility survey, including observations with {\it JWST}, {\it ALMA}, {\it HST}, {\it Chandra} and {\it Gemini}. As a starting point, we report {\it JWST} and {\it ALMA} observations of one of the most massive protostars in the region, S284p1. The observations of shock-excited molecular hydrogen reveal a symmetric, bipolar outflow originating from the protostar, spanning several parsecs, and fully covered by the {\it JWST} field of view and the {\it ALMA} observations of CO(2-1) emission. This allows us to infer that the protostar has maintained a relatively stable orientation of disk accretion over its formation history. The {\it JWST} near-IR continuum observations detect a centrally illuminated bipolar outflow cavity around the protostar, as well as a surrounding cluster of low-mass young stars. We develop new radiative transfer models of massive protostars designed for the low metallicity of S284. Fitting these models to the protostar's spectral energy distribution implies a current protostellar mass of $\sim11\:M_\odot$ has formed from an initially $\sim100\:M_\odot$ core over the last $\sim3\times10^5$ years. Overall, these results indicate that massive stars can form in an ordered manner in low-metallicity, protocluster environments.

Due to non-zero neutrino rest masses we expect the energy density today in non-relativistic matter, $\omega_{\rm m}$, to be greater than the sum of baryon and cold dark matter densities, $\omega_{\rm cb}$. We also expect the amplitude of deflections of CMB photons due to gravitational lensing to be suppressed relative to expectations assuming massless neutrinos. The combination of CMB and BAO data, however, appear to be defying both of these expectations. Here we review how the neutrino rest mass is determined from cosmological observations, and emphasize the complementary roles played by BAO and lensing data in this process. We explain why, for current constraints on the sum of neutrino masses, the addition of BAO data to primary CMB data is much more informative than the addition of CMB lensing reconstruction data. We then use a phenomenological model to find that the preference from CMB and BAO data for a matter density that is below expectations from the CMB alone is at the $3\, \sigma$ level. We also show that if a fraction of the dark matter decays to dark radiation, the preference for $\omega_{\rm m} > \omega_{\rm cb}$ can be restored, but with a small increase to the CMB lensing excess.

Recently, the KM3NeT Collaboration released the detection of a 220 PeV neutrino from the celestial coordinates RA=94.3\degree~ and Dec.=-7.8\degree~ on 13 February 2023 at 01:16:47 UTC \cite{KM3NeT:2025npi}. The source for this extra-ordinary cosmic neutrino, designated KM3-230213A, is not identified yet but there has been a speculation that it might be associated with a gamma-ray burst GRB~090401B \cite{Amelino-Camelia:2025lqn}. The purpose of this report is to search the association of this 220 PeV neutrino with potential GRB sources from a more general consideration of Lorentz invariance violation (LV) without fixed LV scale. We try to associate this extra-ordinary neutrino with potential GRBs within angular separation of 1\degree, 3\degree~ and 5\degree~ respectively and the results are listed in Table 1. We find the constraints $E_{\rm{LV}}\leq 5.3\times 10^{18}$~GeV for subluminal LV violation and $E_{\rm{LV}}\leq 5.6\times 10^{19}$~GeV for superluminal LV violation if KM3-230213A is a GRB neutrino.

We consider tidal masses and ages of Milky Way open clusters, as well as a simple model of their distribution. Our aim is to investigate the space of model parameters and the correspondence between modelled and observed two-dimensional cluster age-mass distributions. The model for cluster evolution is comprised of a two-section cluster initial mass function, constant cluster formation rate, and a mass loss function. This mass loss function represents a supervirial phase after gas expulsion, mass loss due to stellar evolution, and gradual dissolution driven by internal dynamics and the Galactic tidal field. We construct different estimators of model fitness based on $\chi^2$-statistics, the Kullback-Leibler divergence (KLD) and a maximum-likelihood approach. Using these estimators and Markov Chain Monte Carlo sampling, we obtain best-fit values and posterior distributions for a selection of model parameters. The KLD returns a superior model compared to the other statistics. The cluster initial mass function is well constrained and we find a clear signature of an enhanced cluster mass loss in the first 50 Myr. In the KLD best model, clusters lose 72% of their initial mass in the violent relaxation phase, after which cluster mass loss slows down, allowing for a relatively low rate of cluster formation of $0.088\mathrm{M_\odot kpc^{-2} Gyr^{-1}}$. The observed upper limit of cluster ages at approx. 5 Gyr is reflected in the model by a shallow lifetime-mass relation for clusters with initial masses above $1000\mathrm{M_\odot}$. The application of the model to an independent cluster sample based on Gaia DR3 data yields similar results except for a systematic shift in age. The observed cluster age-mass distribution is compatible with a constant cluster formation rate. The enhanced number of young massive clusters observed requires an early violent relaxation phase of strong mass loss.

We investigate the velocity structure of nebular gas in the central galaxies of four clusters: Abell 1835, PKS 0745-191, Abell 262, and RXJ0820.9+0752, using data from the Keck Cosmic Web Imager (KCWI). Velocity structure functions (VSFs) of the [OII] emission line are compared to VSFs of molecular clouds observed with the Atacama Large Millimeter/submillimeter Array (ALMA). Apart from Abell 262 where the gas is located in a circumnuclear disk, the nebular gas in the remaining galaxies lies in off-nuclear filamentary structures with VSFs steeper than the Kolmogorov slope. This steepening may be plausibly attributed to gravity although other factors, such as magnetic stresses and bulk motion,} may be significant. The VSFs of CO and [OII] emission are similar in RXJ0820 and Abell 262, indicating close coupling of the nebular and molecular gases. In contrast, the nebular and molecular gases are differentiated on most scales in PKS 0745 and Abell 1835. This discrepancy is likely due to the radio-AGN churning the gas. We compare the scale-dependent velocity amplitudes of the hot atmospheres constrained by X-ray surface brightness fluctuation analysis using Chandra observations to the nebular VSFs. The large-scale consistency in Abell 1835 and RXJ0820 is consistent with condensation from the hot atmospheres. {We explore substantial systematic biases, including projection effects, windowing, and smoothing effects when comparing VSFs using different telescopes and instruments.

We discuss the predictions in the simplest theory for neutrino masses based on the spontaneous breaking of local lepton number. This theory provides a simple theoretical framework to understand the possible relation between the origin of neutrino masses and the nature of the dark matter. In this theory, one of the fields needed for anomaly cancellation is a dark matter candidate and the local lepton number is broken at the low scale. We discuss in great detail the dark matter properties showing the allowed parameter space by the relic density bounds and the predictions for direct detection. The predictions for gamma and neutrino lines from dark matter annihilation are investigated. In the case of Dirac neutrinos, the bound on the effective number of relativistic degrees of freedom plays an important role and the predictions for gamma lines could be tested in the near future. We discuss the predictions in the case of Majorana neutrinos where the dark matter candidate has extra annihilation channels and compare all the predictions to the case with Dirac neutrinos.

Quasi-Periodic Oscillations (QPOs) are an important phenomenon commonly observed in the X-ray radiation of black holes and neutron stars, closely related to the dynamics of accretion disks around compact objects and general relativistic effects. The objective of this study is to use the QPO phenomenon to distinguish between dark matter-black hole systems and naked singularities, as well as to investigate the effects of different dark matter models (Cold Dark Matter, CDM, and Scalar Field Dark Matter, SFDM) on the accretion disk dynamics. By introducing a dark matter correction model within the framework of general relativity, we systematically investigate the differences in dragging effects, characteristic frequency distribution, and the innermost stable circular orbit (ISCO) radius between dark matter-black hole systems and naked singularities, while analyzing the potential coupling between QPO frequencies and dark matter distribution. The main results of this study are as follows: $\nu_r$ and $\nu_\theta$ in dark matter-black hole systems can be identified as HFQPOs, while for lower spins ($a < 0.5$), $\nu_\text{nod}$ can be identified as LFQPOs, and for higher spins ($1 > a \geq 0.5$), $\nu_\text{nod}$ falls within the HFQPO observation range. Cold Dark Matter (CDM) and Scalar Field Dark Matter (SFDM) modulate the accretion disk dynamics at the order of $10^{-6}$.

Lorentz symmetry is a cornerstone of both the General relativity and Standard Model and its experimental verification deepens our understanding of nature. This paper focuses on the investigation of Lorentz violations with the context of clock comparison experiments in the framework of Standard Model Extension (SME). Considering matter-gravity coupling sector, we provide a generic frame to study the sensitivities of Lorentz-violating coefficients for three distinct types of clock redshift tests, including the traditional gravitational redshift test, null-redshift test I and null-redshift test II. Each of these tests is sensitivity to different combinations of Lorentz-violating coefficients. By using the current clock comparison results, we estimate the limits of SME coefficients at level of parts in $10^{4}$ down to parts in $10^{7}$. Better sensitivity may be achieved in the clock comparisons by using the state-of-the-art optical clocks. Additionally considering relativistic factors in null-redshift I, the frequency comparison result of E2 and E3 transitions of Yb$^{+}$ can set the limit $c^{e}_{00}=(7.4\pm9.3)\times10^{-9}$ in the electron sector. Our analysis demonstrates that clock-comparison redshift experiments may contribute to explore the vast parameters space on searching for the Lorentz violation.

Focusing on supercooled phase transitions in models with classical scale symmetry, we formulate a state-of-the art framework for computing the bubble-nucleation rate, accounting for the presence of various energy scales. In particular, we examine the limitations of derivative expansions in constructing a thermal effective field theory for bubble nucleation. We show that for gauge field fluctuations, derivative expansions diverge after the leading two orders due to the strong variation in gauge field masses between the high- and low-temperature phases. By directly computing these contributions using the fluctuation determinant, we capture these effects while also accounting for large explicit logarithms at two loops, utilising the exact renormalisation group structure of the EFT. Finally, we demonstrate how this approach significantly improves nucleation rate calculations compared to leading-order results, providing a more robust framework for predicting gravitational-wave signals from supercooled phase transitions in models such as the SU(2)cSM.

Heavy sterile neutrinos can be produced in core-collapse supernovae (CCSNe), which are superb particle generators because of their high densities and temperatures. If the sterile neutrinos are long-lived, these may be produced inside the supernova core and escape the stellar envelope, later decaying into SM particles like photons and neutrinos. In this work, we first improve the calculation of the $\gamma$-ray fluxes. We then revisit the bounds on the sterile neutrino parameter space from the non-observation of $\gamma$-rays from SN1987A by the Solar Maximum Mission (SMM) and constraints from the diffuse $\gamma$-ray background arising from sterile neutrino decays. We find that the constraints arising from both the SMM data and the diffuse $\gamma$-ray background are weaker than those that have previously appeared in the literature. Finally, we study the sensitivity of several present and near-future $\gamma$-ray telescopes such as e-ASTROGAM and Fermi-LAT, assuming a nearby future galactic CCSN. We show that future observations can probe mixing angles as low as $|U_{\tau/\mu4}|^2\sim 5\times10^{-17}$.

We report the measurement of the electric dipole moment of aluminum monochloride (AlCl) using a cryogenic buffer-gas beam source. Our measurements provide values for the dipole moments of the two lowest vibrational states of the $X^1\Sigma^+$ and the $A^1\Pi$ electronic states. We also show that spin-orbit coupling with an extended number of spin states is essential in the ab initio calculation to correctly describe both the dipole moment and the Te energy of AlCl. We further lay out the implications of these results for astrophysical models of stellar and planetary evolution that have used a substitute value for the dipole moment of AlCl until now.

Feebly interacting particles, such as sterile neutrinos, dark photons, and axions, can be abundantly produced in the proto-neutron star (PNS) formed in core-collapse supernovae (CCSNe). These particles can decay into photons or charged leptons, depositing energy outside the PNS. Strong bounds on new particles can thus be derived from the observed luminosity of CCSNe, with even tighter bounds obtained from low-energy SNe observations. For the first time we highlight that, at sufficiently large couplings, particle production \textit{outside} the PNS must also be considered. Using the prototypical case of axions coupling to two photons, we show that at large couplings the energy transfer from PNS to its surroundings is diffusive rather than ballistic, substantially reducing the deposited energy. Our findings have implications for the parameter space of particles probed in beam dump experiments and for dark matter models involving a sub-GeV mediator.

Finite mean free paths of light particles, like photons and neutrinos, lead to dissipative effects and damping of small-scale density fluctuations in the early universe. We study the impact of damping on the spectral density of gravitational waves induced by primordial fluctuations in the radiation-dominated universe. We show that the most important effects of damping are $(i)$ regularization of the resonant frequency and $(ii)$ a far low-frequency tail with no logarithmic running. The exact location of the break frequency below which the logarithmic running is lost depends on the damping rate. Both effects stem from the effective finite lifetime of the gravitational wave source caused by damping. Interestingly, we find that, for the standard model of particles, the effects of damping are most relevant at around or below the nHz frequencies. Our results showcase the importance of including the damping of primordial fluctuations in future analysis of induced gravitational waves. We provide detailed analytical formulas and approximations for the kernel of induced gravitational waves. Lastly, we discuss possible implications of damping in alleviating the gauge issue of induced gravitational waves and in suppressing the so-called poltergeist mechanism.

We investigate the dynamics of charged particles in the spacetime of a global monopole swallowed by a Reissner-Nordström (RN) black hole in the presence of a external weak asymptotically homogeneous magnetic field. We carefully analyze and deduce the conditions to have such a magnetic field around this black hole and show that this is indeed possible in the small but nontrivial charge and monopole term limit. We obtain general equations of motion and analyze them for special cases of circular orbits, focusing on the inner-most stable circular orbit (ISCO) of this configuration. The richness of the parameters and complicated forms of the resulting equations of motion necessitate a numerical approach. Hence, we have presented our results with numerous graphs, which help to understand the evolution of ISCO as a function of the external test magnetic field and the monopole term depending on the parameters of the black hole, such as its electrical charge as well as the properties of the test particle such as its specific charge, angular momentum, and energy. We have also analyzed the effective potential that these fields generate and deduced results for the aforementioned values of external and internal parameters of spacetime.

Infrasound sensing offers critical capabilities for detecting and geolocating bolide events globally. However, the observed back azimuths, directions from which infrasound signals arrive at stations, often differ from the theoretical expectations based on the bolide's peak brightness location. For objects with shallow entry angles, which traverse longer atmospheric paths, acoustic energy may be emitted from multiple points along the trajectory, leading to substantial variability in back azimuth residuals. This study investigates how the entry angle of energetic bolides affects the back azimuth deviations, independent of extrinsic factors such as atmospheric propagation, station noise, and signal processing methodologies. A theoretical framework, the Bolide Infrasound Back-Azimuth EXplorer Model (BIBEX-M), was developed to compute predicted back azimuths solely from geometric considerations. The model quantifies how these residuals vary as a function of source-to-receiver distance, revealing that bolides entering at shallow angles of 15 degrees can experience average residuals of 12.4 degrees, with deviations reaching up to 155 degrees at distances below 1500 km, and remaining significant even at 10,000 km. In contrast, bolides with steeper entry angles (greater than 60 degrees) show much smaller deviations, typically under 10 degrees at 1000 km and diminishing to less than two degrees beyond 5000 km. These findings attest to the need for careful interpretation when evaluating signal detections and estimating bolide locations. This work is not only pertinent to bolides but also to other high-energy, extended-duration atmospheric phenomena such as space debris and reentry events, where similar geometric considerations can influence infrasound arrival directions.

The thermal relaxation time of neutron stars, typically defined by a sudden drop in surface temperature, is usually on the order of 10 to 100 years. In this study, we investigate neutron star thermal relaxation by incorporating nucleon superfluidity and non-nucleonic particles, specifically considering hyperons as a representative case. We find that rapidly cooling neutron stars driven by neutron superfluidity and direct Urca processes demonstrate delayed thermal relaxation under specific physical conditions. The former acquires that the neutron $^3P_2$ critical temperature is small enough, whereas the latter depends on the presence of a small core that permits direct Urca processes. To explore these scenarios, we propose simple theoretical frameworks to describe these delayed thermal relaxation behaviors and discuss how an recently-established enhanced modified Urca rate influences the relaxation time. By confronting the theoretical results with the observation of Cassiopeia A, we can effectively constrain the maximum neutron $^3P_2$ critical temperature.

This work investigates a single-field inflationary model, a specific class of the K-essence models where a coupling term exists between canonical Lagrangian and the potential. This coupling term has many effects on key inflationary parameters consisting of the power spectral, the spectral index, the tensor-to-scalar ratio, the Hubble parameter, the equation of state parameter, and the slow-roll parameter. By solving the equations numerically and deriving analytical results, how this modification affects inflationary dynamics can be analyzed. Our results show that the coupling term, $\alpha$, decreases the inflationary parameters, such as the tensor-to-scalar ratio, $r$, and improves the consistency with observational constraints from Planck and BICEP/Keck at the $68 \%$ and $95 \%$ confidence. These findings indicate that the studied model provides a promising alternative to the early universe dynamics while aligning with recent cosmological observations.

Due to the progressive increase in size of the latest Cherenkov-type detectors, it is becoming increasingly important to design a suitable compensation system based on coils of the Earth's magnetic field to ensure the correct operation of the photomultipliers (PMTs). Until now, most studies have assessed the correct functioning of such a system by the proportion of PMTs experiencing more than 100 mG of magnetic field perpendicular to their axis. In the present study, we discuss whether this evaluation parameter is the most appropriate and propose the average residual perpendicular magnetic field $<B_{perp}>$ as an alternative that more closely reflects the loss of detection efficiency of PMTs. A compensation system design is also proposed that offers good results as well as being economical to optimise this parameter.

Super-strongly magnetized plasmas play a crucial role in extreme environments of magnetar and laboratory laser experiments, demanding comprehensive understanding of how quantum electrodynamic (QED) effects influence plasma behaviour. Earlier analytical and semi-analytical calculations have shown that QED effects can significantly modify the plasma polarization mode behaviour around magnetars using analytical and semi-analytical calculations. In this work, we present the first electromagnetic field solver that is valid beyond the Schwinger limit. QED vacuum polarization in super-strong magnetic fields are modeled with nonlinear Maxwell equations. We show that electromagnetic waves in simulations follow the analytical solutions well and reproduce the birefringence effects of electromagnetic wave modes between the $O$ and $X$ polarizations of perpendicular electromagnetic waves and those between $L$ and $R$ polarizations of parallel waves. This new framework can be applied to kinetic as well as in other types of computer simulations. The solver's key advantage lies in its versatility, allowing it to be used in gyro-motion, gyro-center, and gyro-kinetic simulations, which do not resolve the cyclotron motion, or in plasma studies with ground-level Landau quantization.

Current data on ultra-high-energy (UHE) cosmic rays suggest they are predominantly made of heavy nuclei. This indicates that the flux of neutrinos produced from proton collisions on the cosmic microwave background is small and hard to observe. Motivated by the recent extremely-high-energy muon event reported by KM3NeT, we explore the possibility of enhancing the energy-flux of cosmogenic neutrinos through nuclear photodisintegration in the presence of new physics. Specifically, we speculate that UHE neutrons may oscillate into a new state, dark (or mirror) neutron $n'$ that in turn decays injecting large amount of energy to neutrinos, $n\to n'\to \nu_\text{UHE}$. While this mechanism does not explain the tension between the KM3NeT event and null results from IceCube, it reconciles the experimental preference for a heavier cosmic ray composition with a large diffuse cosmogenic flux of UHE neutrinos.

We define a relativistic version of the global symmetries responsible for the restricted mobility of fracton quasiparticles. The theories have a symmetry current that is proportional to a vector field that spontaneously breaks Lorentz boost symmetry. We argue that the existence of a pressureless dust in the early universe could be a consequence of this symmetry. We provide an example of a fractonic scalar field with a quartic self-interaction evolving on a Friedmann-Robertson-Walker background and show that the interaction gives rise to a separately conserved fluid with equation of state $w=1$.