daily-arxiv

The quest for the origin of cosmic ray (CRs) is a fundamental issue in astrophysics. Shocks of supernova remnants (SNRs) have been considered as the dominant contributors to Galactic CRs below the spectral knee near $\sim 3$ petaelectronvolt (PeV). Whether SNRs are efficient accelerators of particles beyond PeV energies has long been debated. Here we report observations of very-high-energy $γ$-ray emission up to hundreds of TeV from two middle age shell-type SNRs, G150.3$+$4.5 and $γ$-Cygni, with the Large High Altitude Air Shower Observatory (LHAASO). Two (or three) distinct morphological/spectral components with convex spectral shapes are observed in both sources, with the low-energy one being more extended than the high-energy one. %Although it is possible that these high-energy components may be driven by powerful pulsars, The likely association of the high-energy component with molecular clouds at similar distances, and the weakness/absence of pulsar wind nebulae (PWNe) inside these SNRs clearly indicate for the first time that the highest energy emission is produced by collision of hadronic CRs up to PeV energies with the clouds. These results are compatible with the classic model prediction that PeV particles accelerated near the end of the free expansion phase of SNR evolution can illuminate nearby molecular clouds (MCs) to produce strong $γ$-ray emission.

Ringdown signals from perturbed black holes (BHs) offer a clean window into BH spacetime, strong-field gravity, and fundamental physics. Presently the quasi-normal modes of stellar-mass BH ringdowns have been successfully extracted in the ground-based gravitational wave (GW) observations. Looking ahead, the future space-borne observatories will listen to the ringdowns from massive BH binary coalescences more loudly and resolve multiple modes to unprecedented precision, which calls for efficient approaches to mitigate the sharply increasing computational burden. We develop a practical ringdown analysis pipeline for space-borne detectors by implementing FIREFLY, a novel acceleration algorithm validated in ground-based detectors, and for the first time demonstrate its compatibility and effectiveness with the time-delay interferometry (TDI) observables. With high fidelity, we achieve a $\sim 200$-fold speedup for a simulated ringdown signal including six modes, providing a viable and scalable route for multi-mode ringdown analysis in the space context. This new approach has sound statistical interpretation and is extensible to other GW sources in band.

We investigate the thermodynamic properties of static nonlinearly scalarized black holes in Einstein-scalar-Gauss-Bonnet theory with polynomial coupling functions. Based on the scalarized solutions constructed previously, we compute thermodynamical quantities of these scalarized black holes. Moreover, we examine the first law of black hole thermodynamics and consider the phase transitions between Schwarzschild and scalarized black holes. It shows that a phase transition from Schwarzschild black hole to scalarized black hole is a first-order with non-zero latent heat.

We present perturbative quantum chromodynamics (pQCD) predictions for high-momentum particle yield modification in very light ion collisions - ${}^{10}\mathrm{B}+{}^{10}\mathrm{B}$, ${}^{6}\mathrm{Li}+{}^{6}\mathrm{Li}$, ${}^{4}\mathrm{He}+{}^{4}\mathrm{He}$, and ${}^{3}\mathrm{He}+{}^{3}\mathrm{He}$ - with and without medium-induced energy loss. We find non-trivial suppression in symmetric systems from ${}^{208}\mathrm{Pb}+{}^{208}\mathrm{Pb}$ to ${}^{3}\mathrm{He}+{}^{3}\mathrm{He}$ and in asymmetric $A+B$ systems, with the suppression scaling approximately as $R_{AB} \simeq (\sqrt{AB})^{1/3}$. Further, we find that ${}^{3}\mathrm{He}$ and ${}^{6}\mathrm{Li}$ offer particularly clean environments for observing final-state partonic energy loss from quark-gluon plasma (QGP) formation in extremely small systems. Finally, we show that energy loss models generically predict $v_2\{\mathrm{SP}\} \approx 0$ in small systems, indicating that the large measured $v_2 > 0$ in $p+{}^{208}\mathrm{Pb}$ is not due to energy loss.

We investigate the properties of dense quark matter and strange quark stars within a nonperturbative, Poincaré-covariant framework. Employing a symmetry-preserving vector$\,\otimes\,$vector contact interaction model, we extend the quark gap equation to the regime of zero temperature and finite quark chemical potential. From the resulting momentum-independent quark propagator, we construct the equation of state (EOS) and solve the Tolman-Oppenheimer-Volkoff (TOV) equations to evaluate the mass-radius relations and tidal deformabilities of strange quark stars. We systematically analyze the sensitivity of the EOS and the macroscopic stellar properties to the model parameters, specifically the effective interaction strength and the ultraviolet cutoff. We demonstrate that reducing the coupling constant stiffens the EOS, whereas increasing the ultraviolet cutoff softens it. By confronting our predictions with multi-messenger astrophysical constraints-including pulsar mass measurements and gravitational-wave data-we identify parameter regimes that successfully describe current observations. Specifically, we find that parameter sets with $α_{ir}=0.735π$, $Λ_{uv}=0.905\,\mathrm{GeV}$ and $α_{ir}=0.588π$, $Λ_{uv}=0.9955\,\mathrm{GeV}$, alongside a vacuum bag pressure of $B \approx (0.106\,\mathrm{GeV})^4$, yield stellar properties in excellent agreement with empirical constraints.

The IceCube Neutrino Observatory has opened a new window into the high-energy Universe, providing measurements of neutrinos over a broad energy range. This contribution presents recent results, including a follow-up on the first identification of a steady neutrino source NGC 1068, measurements of the flavor composition of the diffuse astrophysical flux, limits on prompt atmospheric neutrinos, and searches for neutrinos from dark matter annihilation in the Sun. These measurements probe neutrino production mechanisms, fundamental particle interactions, and physics beyond the Standard Model. Looking forward, the recently deployed IceCube Upgrade will enhance sensitivity to lower-energy neutrinos and reduce systematic uncertainties, while the planned IceCube-Gen2 will expand the detector volume, increase the neutrino detection rate, and extend energy reach, enabling more detailed studies of cosmic sources and high-energy particle physics.

This work describes a few-body dynamics method based on the Faddeev integral equations in momentum space for determining the total cross sections of fusion and breakup reactions with two- and three-body final channels in the continuum, employing a cluster representation of the colliding nuclei. Total cross sections were obtained for the reactions $^3\text{He}(T,D)^4\text{He}$, $^3\text{He}(T,np)^4\text{He}$, $^3\text{He}(T,nD)^3\text{He}$, $^3\text{He}(^3\text{He},2p)^4\text{He}$, $^3\text{He}(^3\text{He},pD)^3\text{He}$, $^7\text{Li}(^3\text{He},\phantom{0}^4\text{He})^6\text{Li}$, $^7\text{Li}(^3\text{He},D^4\text{He})^4\text{He}$, and $^7\text{Li}(^3\text{He},T^3\text{He})^4\text{He}$, in which both the projectile and the target nucleus were treated in a cluster representation. The work also implements a two-potential method to determine the Coulomb $t-$matrix and to account for Coulomb effects in short-range dynamics in momentum space. Calculations of the initial-state Coulomb interaction were performed; furthermore, an estimate was obtained for the magnitude of the off-shell effect of the Coulomb $t-$matrix, as well as the magnitude of the atomic electron anti-screening effect on the Coulomb interaction of the colliding nuclei. The calculated contribution of the cluster mechanism to the total cross section of the considered fusion reactions in the kinetic energy range $T\in[1~\text{keV},20~\text{MeV}]$ is in good agreement with known experimental data.

We investigate the dynamics of transonic advective accretion flows around spinning black holes in the presence of viscosity. The spacetime of a Kerr black hole is approximated using a pseudo-potential. We study viscously driven shock oscillations over a range of black hole spin parameters. Our results show that the frequency range of quasi-periodic oscillations (QPOs) obtained from the power density spectra depends strongly on the black hole spin. Low-spin systems predominantly exhibit low-frequency QPOs, whereas rapidly rotating black holes (0.9) produce QPOs spanning a broad range from low to high frequencies, comparable to those observed in black hole X-ray binaries. We further obtain a correlation between the QPO frequency and the power-law photon index by computing the spectrum for a 10 solar mass black hole.

We use the two-flavor Linear Sigma Model with quarks as an effective description of QCD to investigate the nature of the chiral phase transition at finite baryon chemical potential and zero temperature. We work at one-loop order to set up and solve the system of self-consistent coupled equations for the particle pole masses. The chemical potential-dependent value of the chiral order parameter is obtained by minimizing the one-loop effective potential. This treatment goes beyond the conventional ring-diagram approximation and provides a description valid for arbitrary values of the chemical potential. We find that the phase transition is of first order, and occurs when the quark chemical potential reaches the value of the vacuum quark mass for the chosen set of parameters. The first order nature of the transition is signaled by the discontinuous behavior of the chiral condensate, the masses and the couplings. The thermodynamics of the system is readily implemented and in particular, we find that the square of the speed of sound exhibits a discontinuity at the phase transition and then smoothly approaches the conformal limit from below.

There has been a recent revival in understanding the spontaneous scalarization phenomenon in scalar-tensor gravity as a phase transition. Using the tools of the Landau theory, we now know that first-order transitions where scalarization occurs in a discontinuous manner is more prominent than what had been considered in the literature, and this might lead to novel observation channels. However, the examples so far have been restricted to specific quadratic scalar coupling terms and spherically symmetric stars. Here we explore the phase transition structure of scalarization for more general couplings, considering linear as well as quadratic terms in the conformal scaling factor of the theory. Moreover, we also investigate the effect of rotation on the scalarization phase transition. Both of these considerations are natural choices since the coupling in a scalar-tensor theory can appear at all orders, and astrophysical neutron stars commonly have angular momentum. The introduction of linear coupling leads to a complex solution space which is harder to explore. However, we demonstrate that the Landau model of scalarization enables us to systematically find the branches of scalarized solutions that are commonly overlooked in numerical searches, providing a novel tool. On the other hand, the main effect of stellar rotation is shifting the stellar masses at which the phase transition occurs to higher values, but the qualitative picture remains similar to what happens under spherical symmetry.

Compact stars serve as natural systems where matter exists at densities far beyond those achievable in laboratory experiments. Among them, magnetars are expected to possess interior magnetic fields that may reach values of the order of $10^{17}-10^{18}$ G. These extreme conditions are expected to alter the microscopic and macroscopic properties of dense matter. In this review, we examine how strong magnetic fields affect fermionic matter through mechanisms such as Landau quantization and anomalous magnetic moment interactions. We further discuss the behaviour of magnetized hadronic matter within relativistic mean-field approaches and consider the possible emergence of additional degrees of freedom, including hyperons, $Δ$ resonances, meson condensates and quark matter. The consequences of these effects for neutron-star structure and observational constraints are also briefly outlined.

The paper develops a model to understand the effective quantum geometry of a black hole horizon and the emission of Hawking spectrum in 2+1 dimensions. We argue that one may view the black hole horizon as formed out of quantised lengths of elementary quanta of value $8π\ell_{P}\, n$, where $n\in \mathbb{N}$, and $\ell_{P}$ is the Planck length. To an observer near the black hole horizon, the entropy (or length of horizon cross-section) is related to the black hole energy. Hence, one may develop a formulation of length ensemble (similar to the area canonical ensemble of Krasnov) from which the black body spectrum may be obtained directly. To this local observer, the temperature of the Hawking spectrum is modified due to the Tolman factor.

X-rays are emitted from the corona above the orbiting matter of the accretion disk and travel either directly to us or illuminate the disk. This illumination of the inner disk is enhanced by gravitational light bending, which focuses the rays towards the black hole and therefore towards the inner radii of the disk. These rays that hit the inner radii are reflected back to us, and we observe them in the X-ray reflection spectrum. In this work, we create novel general relativistic ray tracing simulations to investigate the effects of altering the geometry of the accretion disks of black holes on the most dominant part of the reflection spectrum, the iron K$α$ line. Work demonstrating the effect of disk geometry on the iron line has been performed, though many previous analyses have assumed a simplistic system, consisting of a point-source corona with a flat and infinitesimally thin accretion disk. We extend these models to more realistic accretion disk approximations. These include a constant aspect ratio disk, a radiation-pressure-dominated Shakura-Sunyaev disk, an expanded inner disk that has a non-negligible scale height in its inner regions due to radiation pressure, as well as various warped disks. Using measurement uncertainties from XRISM, we find that non-negligible thickness in accretion disks underestimates the black hole spin, corona height, and inclination angle if fitted with a flat disk model. The warped disk model could not be fit with the flat disk approximation.

The ground-state properties of superheavy $Z = 122$ isotopes are investigated using the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc). Bulk properties, including binding energies, Fermi energies, nucleon separation energies, quadrupole deformations, and root-mean-square radii, are calculated. The results are compared with those obtained from the relativistic continuum Hartree-Bogoliubov (RCHB) theory. By examining the dependence on the angular-momentum cutoff and the effects of triaxial and octupole deformations, a strategy for determining the ground states is suggested. Furthermore, based on an analysis of the Fermi and nucleon separation energies, the proton and neutron drip lines for $Z = 122$ isotopes are determined within both the DRHBc and RCHB frameworks. The possible magic numbers $N=184$, 258, and 350 are also suggested. Finally, the evolution of single-particle levels, deformation, charge and neutron radii as well as average pairing gaps with increasing neutron number, is discussed.

Supernovae are one of the most promising gravitational wave sources. But, since the system of the supernovae is nearly spherically symmetric, the expected gravitational waves from them are relatively weak, compared to the case of the compact binary mergers. Thus, at least using the current gravitational wave detectors, only the gravitational waves from a supernova that occurred in our galaxy could be detected. To reliably extract information from gravitational waves originating from such a low event rate, thorough preparation is essential. However, because supernova gravitational waves strongly depend on model parameters, such as progenitor mass and the equation of state for dense matter, it may be difficult to extract physical properties even if the gravitational waves are detected. The universal relations between gravitational-wave signals and physical properties, independent of model parameters, are important for solving this difficulty. To discuss such a universal relation, in this article, we systematically examine the protoneutron-star oscillation frequencies with the linear analysis, the so-called asteroseismology, and compare them with the gravitational wave signals in the simulations.

Fast radio bursts (FRBs) are millisecond-duration extragalactic radio transients of unknown origin. Rotation measures (RMs) probe their local magneto-ionic environments and provide important clues to their nature. While RM variability has been observed in several repeating FRBs, it is typically gradual or stochastic. Recently, observations of FRB~20220529 revealed an abrupt RM excursion followed by rapid recovery on week-long timescales, termed an ``RM flare'', suggesting a potentially distinct form of RM variability associated with localized magnetized plasma. In this work, we perform a systematic search for RM flare candidates in repeating FRBs with multi-epoch RM measurements. Using a $3σ$ significance threshold, we identify two candidates with multiple observational epochs (FRB~20121102A and FRB~20201124A) and two additional single-epoch candidates (FRB~20180916B), in addition to the event in FRB~20220529A. Our results suggest that RM flares, if confirmed, may not be rare among repeating FRBs and point to highly dynamic magnetized environments local to the sources. Future high-cadence polarimetric observations, particularly following the discovery of RM excursions, will be essential for confirming these candidates and constraining their physical origin.

The gravitational wave event GW231123, with component black hole masses lying within or above the pair-instability mass gap, poses a significant challenge to current stellar evolution models. In this Letter, we investigate its origin by coupling the galaxy formation model GAMESH with the cluster population synthesis code RAPSTER, and with two distinct binary population synthesis codes (SEVN and BSEEMP). This framework allows us, for the first time, to reconstruct the life cycle of GW231123-like candidates within the same cosmological simulation, enabling a self-consistent comparison between different formation channels. We find that, although both population synthesis codes can in principle produce black holes compatible with GW231123, isolated binary evolution fails to reproduce the inferred merger redshift. In SEVN, massive black hole binaries form with semi-major axes > 10^3 Rsun , preventing coalescences within a Hubble time. In BSEEMP, candidates arise only at extremely low metallicities (Z = 10^{-10}), which contribute negligibly to the star formation rate density in our overdense simulated volume. Our results instead strongly support a dynamical, hierarchical origin. The observed black hole masses are naturally reproduced through successive mergers in dense globular clusters. The high dimensionless spins reported by the LIGO-Virgo-KAGRA Collaboration are consistent with this hierarchical population. We find a local merger rate density of 0.78 Gpc^{-3} yr^{-1}, with a peak at z = 4 - 6, tracing the maximum formation rate of globular clusters in metal-poor environments (Z = 0.006). Overall, GW231123 may represent a benchmark event for a robust population of hierarchical black holes formed in the early Universe.

The supermassive black hole Sagittarius A* (Sgr A*) exhibits temporal and spectral variability across the electromagnetic spectrum. However, variability at radio frequencies below ~ 5 GHz for timescales shorter than a day remains largely unexplored. We investigate the variability of Sgr A* at 2.79 GHz on short timescales (1 min), to probe an under-explored regime of its emission process. Through point-source model fitting in the uv-domain, we analyse the flux density variation of Sgr A* over an 8 h observation. We detect flux variation on a few tens of minute timescale with a modulation index of 6.11 %, a mean flux density of ($827 \pm 0.1_{\mathrm{stat}} \pm 33_{\mathrm{sys}}) \, \mathrm{mJy}$, and a mean spectral slope of $0.08\pm0.03$. Furthermore, we measure the slope of the structure function of the observed light curve as $0.81 \pm 0.05$ with a characteristic timescale of about 120 min. Our study at low radio frequencies is a critical step toward constraining the physical mechanisms that drive Sgr A*'s variable emission and its spectral energy distribution. Our study suggests that variability at centimetre and millimetre wavelengths is likely more closely related than previously thought.

Gamma-ray bursts (GRBs) are generally classified as Type~I GRBs, which originate from compact binary mergers, and Type~II GRBs, which originate from massive collapsars. The traditional correspondence between short--Type~I GRBs and long--Type~II GRBs, separated by a duration of 2 seconds, has been challenged by recent observations of long GRBs associated with kilonovae (i.e., Type~I-L GRBs) and a short GRB associated with a supernova. In this paper, we focus on GRBs with precursor emission (PE) and compile 366 GRBs detected by Fermi/GBM. Applying the unsupervised machine learning methods t-SNE and UMAP, we are able to distinguish Type~I (including subclass Type~I-L) and Type~II GRBs for the first time and identify PE as a key feature for distinguishing GRBs of different origins. Inspired by results of machine learning, we propose a diagnostic parameter, the $E_{\rm p,ME}$-precursor index ($EPI$), defined as ${\rm log_{10}}(E_{\rm p,ME}^{2}/(T_{\rm 100,PE}T_{\rm 100,QE1}^{1/2}T_{\rm MVT,PE}))$, where most Type~I GRBs have $EPI > 6.2$ and most Type~II GRBs have $EPI < 6.2$. This parameter can help the community to diagnose the origin of any GRB with PE based solely on its prompt emission and rapidly plan for follow-up observations. The validation using Swift GRBs provides illustrative evidence that our method may also be applicable to GRBs observed by instruments other than Fermi.

Neutrino telescopes have detected astrophysical neutrinos with energies up to ${O}(100)$ PeV. Several current and proposed experiments aim to observe neutrinos at even higher energies, with the goal of detecting cosmogenic neutrinos. This increase in neutrino energy makes tests of Lorentz invariance violation (LIV) particularly appealing, since the effects of higher-dimension LIV operators on neutrino propagation grow rapidly with energy. In this work, we investigate the potential of the upcoming experiments GRAND and POEMMA to probe LIV in the neutrino sector through the detection of ultra-high-energy tau neutrinos. We generate the cosmogenic neutrino flux using SimProp and interface it with a calculation of neutrino flavor transition probabilities in the presence of LIV effects. Deviations from standard flavor transition probabilities manifest as changes in the expected tau neutrino event rates at GRAND and POEMMA. We first consider the case with a single nonzero LIV operator of various dimensions, and find that the projected sensitivities exceed existing limits from lower-energy probes by orders of magnitude. We then explore scenarios with multiple nonzero LIV parameters and show that their interplay can significantly modify the sensitivities compared to the single-parameter case. Overall, we find that upcoming observations of ultra-high-energy tau neutrinos will place some of the most stringent constraints on LIV.

We investigate dense-matter equations of state (EOSs) within a Bayesian framework, with particular emphasis on whether recent small-radius compact-star candidates can be accommodated in a twin-star scenario. For the hadronic sector, we adopt a meta-modeling EOS constrained by the NICER mass--radius measurements of PSR J0030$+$0451, PSR J0437$-$4715, PSR J0614$-$3329, and the massive pulsar PSR J0740$+$6620. The hadronic inference indicates that PSR J0614$-$3329 favors a somewhat softer EOS than the other two \(\sim1.4\,M_\odot\) pulsars, while the \(\sim2\,M_\odot\) constraint prevents the EOS from becoming too soft. We then introduce a strong first-order phase transition through a constant-speed-of-sound quark-matter segment. Using HESS J1731$-$347 and XTE J1814$-$338 to constrain the phase-transition parameters, we find a preferred transition density of \(n_\mathrm{t}\sim2.7\text{--}2.8\,n_0\), a sizable energy-density jump of \(600\text{--}700\) MeV, and a relatively large post-transition sound speed of \(c_s^2/c^2\sim0.85\). Such a phase transition generates a disconnected hybrid branch with radii of about \(6\text{--}7\) km at masses around \(1.2\text{--}1.4\,M_\odot\), and strongly suppresses the dimensionless tidal deformability relative to the purely hadronic branch. This pronounced change in tidal deformability is a characteristic signature of the twin-star mechanism and may provide an important observational tool for identifying phase transitions in neutron-star matter in future multimessenger measurements. These results show that small-radius compact stars can provide direct constraints on both the strength of a first-order phase transition and the stiffness of the post-transition phase in dense matter.

We present the first high-resolution spectroscopy of the black hole high-mass X-ray binary Cygnus X-1 with XRISM, including orbital-phase-resolved analyses and tentative evidence of short-term variability in the Fe-K band on second timescales. Using data from the Performance Verification phase in April 2024, we analyzed spectral variability across orbital phases with the Resolve microcalorimeter and the Xtend CCD imager. The unprecedented resolution of Resolve reveals variability in highly ionized Fe absorption lines. The absorption features show orbital-phase-dependent variability in column density, ionization state, and blueshifted velocity, suggesting structural variations in the focused stellar wind along the line of sight. We also find indications of subtle broadening of the neutral Fe emission profile. In addition, intensity-sorted spectroscopy during dip phases suggests possible variability on timescales of a few seconds in the absorption features, consistent with cooler, denser and lower-ionized gas clumps. Although the statistical significance is limited, these results hint that the stellar wind and the X-rays from the accretion disk around the black hole may interact on timescales as short as a few seconds. These XRISM results constrain wind-fed accretion in Cyg X-1 and highlight Resolve's capability to probe plasma environments in high-mass X-ray binaries.

We present 3D kilonova radiative transfer simulations for a series of binary neutron star merger models. The masses of the neutron stars are varied as well as the total mass of the system and two different equations of state were used (SFHO and DD2), producing a range in dynamical ejecta masses and elemental abundance patterns. In this paper, we focus on the bolometric light curves and spectra in the polar direction for comparison with observations of the kilonova AT2017gfo. We calculate line-by-line opacities and include new calibrated lanthanide atomic data. All of the simulated spectra show strong features from Sr II, La III, Gd III and Ce III, which appear to correspond to features identified in AT2017gfo, although the simulated features are generally more blueshifted. The models with the lowest lanthanide fraction in the polar direction also show a Y II feature. Ce III, Ce II, Nd III and Nd II play an important role in shaping the spectral continuum. While the bolometric luminosities in the polar direction vary with the ejecta mass of each model, we find only little sensitivity of the spectral properties to the merger configuration. Our study demonstrates that dynamical ejecta alone can reproduce (although at earlier times) many spectral properties of AT2017gfo, suggesting dynamical ejecta may have a strong impact on the early spectral evolution. However, future simulations are needed to also elucidate the role of other ejecta components for shaping the kilonova spectrum.

By monitoring the times of arrival of radio pulses from millisecond pulsars, Pulsar Timing Arrays (PTAs) serve as unique gravitational wave (GW) laboratories in the nanohertz band. To date, the primary astrophysical sources of GWs targeted in this frequency range have been inspiraling supermassive black hole binaries (SMBHBs) on circular and eccentric orbits. In this work, we demonstrate that, thanks to the so-called pulsar term in the timing residual waveform of GW signals, PTAs can probe individual SMBHBs that merged before timing observations began. We refer to the latter as \emph{zombie binaries}. Using SMBHB population models consistent with current PTA constraints, we find that while the probability of detecting such systems in existing PTA datasets remains low, the Square Kilometer Array observatory is expected to achieve sufficient sensitivity to have a few zombie binaries with a signal-to-noise ratio exceeding 3 in its data. Although their confident identification might be challenging, this new class of PTA sources opens a novel window for studying the most massive SMBHBs in our local universe.

We present multi-wavelength study of the $γ$/X-ray transient EP250416a (also designated GRB 250416C), triggered by the Einstein Probe (EP) Wide-field X-ray Telescope and also by SVOM and Konus-Wind. Observations spanning the gamma-ray, X-ray, and optical bands facilitated detailed analysis of the burst's prompt emission, afterglow evolution, and physical origin. EP250416a exhibits a burst duration of 30 s in X-ray and 17.7 s in gamma-rays, with joint spectral fitting of 0.5-5000 keV data gives $E\rm_{peak}=342_{-232}^{+90}$ keV. Optical spectroscopy of the afterglow, acquired with the Gemini Multi-Object Spectrograph (GMOS) on Gemini South, yielded a redshift of $z=0.963$. Accounting for the measured redshift, the isotropic energies are $E\rm_{X,iso}=2.7_{-0.5}^{+0.9}\times10^{50}$ erg and $E\rm_{γ,iso}=7.34_{-2.1}^{+5.1}\times10^{51}$ erg, aligning with the Amati relation for long GRBs. The fluence ratio $\rm S(25-50~keV)/S(50-100~keV)=0.78_{-0.15}^{+0.1}$ classifies EP250416a as an X-ray rich (XRR) GRB. The X-ray afterglow shows an initial shallow decay ($α\approx -0.5$) transitioning to a canonical decay phase ($α\approx -1$), with a very late jet break at $t\sim 1.5\times 10^6$ s, corresponding to a jet half-opening angle of $θ_j=10.6_{-1.8}^{+1.9}$ degrees. EP250416a is optically dark, as it shows only a faint $r$-band detection ($r=24.16$ mag) from Gemini South-GMOS and a low optical-to-X-ray spectral index $β_{\rm OX} = 0.3$. This may be attributed to significant host-galaxy extinction, with a required $A_V^{\text{host}}=5.5\ \text{mag}$ derived from the extinction curve model.

We use GRAVITY near-infrared (NIR) flare astrometry to test whether Sagittarius A* could be a solitonic boson star. We consider five spherically symmetric solitonic boson-star models with different effective radii, together with the Schwarzschild black hole. Treating the flares as hot spots on circular equatorial orbits, we analyze their centroid motions and images in these spacetimes and use them for parameter fitting. We perform the fitting using both $χ^2$ analysis and Markov Chain Monte Carlo (MCMC) methods, which yield consistent results: the inferred masses of boson-star models are systematically larger than the established value of $4.3\times10^6M_\odot$. Notably, more diffusive boson stars exhibit imaging properties closer to those of a black hole, leading to mass estimates that are correspondingly closer to the established value. Overall, our results place stringent constraints on solitonic boson star interpretations of Sagittarius A*, although do not completely rule them out.

The shock fronts of supernova remnants (SNRs) are believed to be significant sites of acceleration of cosmic ray particles. Previous researchers have shown that a particle distribution similar to a log-parabola can be generated when particles have an energy-dependent escape. We explore the acceleration of electrons at SNR shock fronts, and show that modeling this energy-dependent particle escape model can produce spectral energy distributions consistent with observations of two lepton-radiation-dominated SNRs: RX J1713.7-3946 and SN 1006. The model includes the evolution of both the electron distribution and photon spectra as a result of the combined effects of the SNR evolution and electron energy loss. The electron-escape energy dependence is quite weak, but the electron distribution and photon spectra turn out to be very sensitive to changes in the electron escape. We also explore how sensitive the spectra and electron distributions are to the parameters used in the log-parabola model.

Pulsar timing arrays record gravitational waves from supermassive black hole binaries at two spacetime points: an Earth term, measured when the wave passes the Earth, and a pulsar term, measured when the wave passed each pulsar at an earlier epoch. We show that a future $μ$Hz-band detection of a nearby massive binary by a mission such as $μ$Ares would turn PTA pulsar terms into targeted probes of binary evolution. In analogy with supernova light echoes, each pulsar term acts as a gravity echo: a dated snapshot of the binary at an earlier stage of its inspiral. Together, the $μ$Hz Earth-term measurement and the nHz pulsar-term echoes provide a temporal baseline that neither detector could access alone. For a fiducial equal-mass binary with total mass $10^9\,M_\odot$ at 80~Mpc, we find a combined pulsar timing array echo signal-to-noise ratio of 33, with up to 24 pulsars individually resolving the signal among pulsars with 50-year baselines. The angular dependence of the single-pulsar echo sensitivity alone enables independent sky localization of the source to $\sim$10--100~deg$^2$, and the resolved pulsar-term frequencies directly measure the binary inspiral rate hundreds to thousands of years ago. With sufficient pulsar distance precision, a small set of anchor pulsars could additionally phase-connect the array and trace the post-Newtonian evolution coherently over kpc baselines. The source population required for gravity echoes is drawn from the same massive-end census responsible for the observed nanoHertz stochastic background.

Around 10 % of hydrogen-poor supernovae explode inside compact ($\sim 10^{15}$ cm), massive ($\sim 0.1 \ \mathrm{M_\odot}$) circumstellar material (CSM), signalling an episode of enhanced pre-explosion mass loss whose mechanism remains unclear. The extreme members of this population are considered to constitute some of the Fast Blue Optical Transients (FBOTs), which exhibit rapid rise times of $\sim$ few days and high peak luminosity $\sim 10^{44} \ \mathrm{erg}$. Recent binary evolution calculations show that the expansion of helium stars during their latest evolutionary stages can trigger a rapid but stable mass-transfer episode that can form a dense circumbinary disc (CBD) that may explain the observed dense CSM. However, a detailed, quantitative analysis of this process and the resulting CBD properties such as its mass, radius and density profile has not yet been undertaken. We present a set of models that solve the viscous evolution of such a CBD under time-dependent mass injection. We find that although the injected mass is initially sub-Keplerian, a lower ``accretion eigenvalue'' $χ$ prevents more mass from falling back onto the central binary. For our fiducial set of models, the CBD immediately prior to the explosion reaches a mass of $0.07-0.20 \ \mathrm{M}_\odot$, a half-mass radius of $640 - 4000 \ \mathrm{R}_\odot$, and an aspect ratio of $θ= H/R \sim 0.1$. We also show that the interaction between SN ejecta and the CBD can power some of the fastest-evolving interacting Type Ibc SNe that can be classified as FBOTs, such as SN 2018gep or SN 2019jc. Despite uncertainties in the model parameters, our results demonstrate that CBD formation triggered by rapid, stable mass transfer is a viable mechanism to explain the dense circumstellar environments observed around rapid, hydrogen-poor interacting SNe. (abridged)

We propose a quantum-enhanced sensing scheme for the detection of wave-like dark matter and high-frequency gravitational waves using two-dimensional ion crystals in a Penning trap. The protocol employs spin-motion squeezed states to improve the signal-to-noise ratio and enable a super-Heisenberg scaling with respect to the number of ions over a broad parameter range. We analyze the sensitivity of the protocol to representative wave-like dark matter candidates, including the axion-like particle and the dark photon, as well as to high-frequency gravitational waves, taking into account the decoherence and dephasing of the ion spins. Our results indicate that two-dimensional ion crystals and this new protocol provide a promising platform for probing previously unexplored parameter space in searches for light dark matter and high-frequency gravitational waves.

This work presents a non-equilibrium framework for thermodynamicized black holes, inspired by the entropy-functional interpretation of emergent gravity and by residue-based methods in black hole thermodynamics. The main idea is to unify three components: an entropy functional principle for selecting physical on-shell backgrounds, a Euclidean and contour-based description of the horizon temperature through simple pole singularities, and a topological residue classification of multi-horizon black hole configurations. On this basis, the paper introduces a quasi-stationary non-equilibrium partition functional in which irreversible entropy production appears as an additional contribution to the singular action. The formalism reproduces the standard equilibrium relations in the adiabatic limit, while also extending them to dynamically driven black-hole systems with matter, charge, and rotational fluxes. The framework is then applied to Kerr Newman type black holes in constant curvature f(R) gravity, where the equilibrium entropy remains weighted by the derivative of f at the background curvature, while non-equilibrium corrections arise from flux-induced deformations of the effective thermodynamic action. The analysis further shows that the outer and inner horizons carry opposite topological orientations, so the non-extremal Kerr Newman family stays in the topological class W = 0 unless a horizon bifurcation or merger changes the singularity structure. Finally, several function plots are provided to illustrate the behavior of equilibrium and non-equilibrium free energy, the Kerr Newman temperature curve, and the entropy production law.

Reducing orbital eccentricity in numerical relativity simulations of binary black holes is essential for producing astrophysically relevant gravitational wave models, as many of these systems are expected to be near-circular in nature. Standard eccentricity reduction procedures rely on iterative schemes, often requiring four or more trial simulations to achieve desired thresholds. This approach is computationally expensive because each trial simulation adds ~10% to the total simulation run time of multiple weeks to months. We introduce a data-driven approach that accelerates this process by learning the values of the initial orbital frequency, Omega_0, and radial velocity, adot_0, that yield an evolution with small eccentricity. This is done using a Gaussian Process Regression model trained on an archive of previously eccentricity-reduced numerical relativity simulations. For all configurations tested, using the trained model consistently reduces the number of required eccentricity reduction iterations to just zero or one, significantly lowering computational costs relative to post-Newtonian initial guesses. These results demonstrate the power of data-driven methods in accelerating expensive numerical relativity simulations.

Ultra high energy cosmic rays have been see coming from the direction of the local cosmic void. We use this fact to argue that at least some of these these cosmic rays are relatively light magnetic monopoles and that their relative fraction above 1020 eV can be found from full sky observations.

X-ray polarimetry is a new tool capable of probing the geometry of accretion onto weakly magnetized neutron stars. Here we present the first X-ray spectropolarimetric results from coordinated observations of the atoll source 4U 1735-44, conducted with the Imaging X-ray Polarimetry Explorer (IXPE), NICER, and NuSTAR. Over the 2-8 keV energy range, we obtained a marginal detection of polarization with the polarization degree of $1.4\%\pm0.7\%$ and polarization angle of $-29^\circ\pm14^\circ$, corresponding to a $3σ$ upper limit on the polarization degree of 3.5\%. The best-fit model to describe the spectrum comprises a thermal component associated with the accretion disk, a Comptonized blackbody component, and a relativistic reflection component. From the reflection model, we infer a disk inclination of $\sim 40^\circ$. The spectroscopic and polarimetric properties of 4U 1735-44 are consistent with those observed in other atoll sources studied by IXPE, with its low polarization likely due to its low inclination.

We propose a cosmological scenario in which, beyond matter and radiation, an additional barotropic fluid with positive equation of state $ω_s$ contributes to the cosmic energy budget, in contrast to Early Dark Energy (EDE). We investigate the theoretical implications of this framework, here dubbed the $Λ_{ω_s}$CDM model, at both the background and perturbative levels, exploring its impact on the expansion history and structure formation. We show that, while remaining subdominant at late times and therefore consistent with current observational bounds, the additional fluid modifies the early-time expansion rate, leading to a higher inferred value of the Hubble constant. Thus, we perform a full Bayesian analysis using a modified version of the \texttt{CLASS} Boltzmann code interfaced with \texttt{MontePython}, considering combinations of \textit{Planck} 2018 Cosmic Microwave Background (CMB) data, DESI DR2 Baryon Acoustic Oscillations (BAO) measurements, Pantheon Type Ia supernovae (SNe Ia), and SH0ES determinations of $H_0$. We find that the inclusion of the SH0ES prior, $H_0 = 73.04 \pm 1.04\,\mathrm{km/s/Mpc}$, leads to a preference for a nonvanishing barotropic fluid. In particular, we obtain $ω_s = 0.290^{+0.017(0.021)}_{-0.007(0.028)}$ and density $10^5Ω_s = 1.47^{+0.35(1.14)}_{-0.62(0.94)}$ for the dataset combination CMB + BAO + Pantheon + SH0ES, and $ω_s = 0.302^{+0.024(0.034)}_{-0.013(0.038)}$ and $10^5Ω_s = 1.21^{+0.31(1.10)}_{-0.65(0.86)}$ when BAO data are excluded. We further compare our scenario with the EDE framework and show that, statistically, no strong evidence is found against the $Λ_{ω_s}$CDM model. Finally, we provide a physical interpretation of our fluid in terms of matter with pressure, indicating that the standard cosmological model may be incomplete in its current minimal formulation.

General relativity is a background-independent theory of a dynamical classical spacetime geometry. Quantum theory is formulated in a classical spacetime, as an intrinsically probabilistic, contextual theory of non-classical, interfering probabilities, with a fixed Born rule for computing those probabilities. We argue that the quantum nature of spacetime, which includes a non-commutative dual companion to the (observed) classical spacetime, is the reason behind an intrinsically probabilistic and contextual nature of quantum theory, with the fixed Born rule. In quantum gravity, we claim, quantum theory is gravitized into a background-independent structure with dynamical and contextual quantum probabilities. This proposal implies intrinsic triple and higher-order interference in the context of massive quantum probes, which sheds light on string theory and the observed vacuum energy as well as the masses of elementary particles.

Little Red Dots (LRDs) have emerged as one of the central puzzles of the JWST era. Their spectra increasingly require dense gas close to the source, yet the physical origin of that cocoon-like structure remains unclear. We examine whether late pulsational mass loss from supermassive stars (SMS)leads to dense gas cocoons. We analyze five accreting GENEC models at different metallicities with characteristic masses of order $10^5\,M_\odot$, following them through post-accretion evolution with radial pulsation calculations and general relativistic (GR) stability diagnostics. Mass loss during the final stages of evolution occurs not as a steady wind, but through discrete strange-mode ejection episodes. In the $Z=10^{-2}\,Z_\odot$ model, which provides the clearest LRD analogue, four late episodes last $41$--$282$ yr and eject $10$--$348\,M_\odot$ each, for a total loss of $(4.8-10)\times10^2\,M_\odot$; the final episode alone contributes $\simeq 73\%$ of that budget. Since the last episode dominates the mass-loss, it is the only event sufficiently massive enough to leave behind a compact, optically thick shell extending out to 0.4 pc that reproduces the LRD dense gas cocoon. The final ejecta are H/He dominated but chemically distinctive, with a robust nitrogen-rich composition, $\log(\mathrm{N/O})\simeq0.13$ and $\log(\mathrm{C/O})\simeq-0.23$. The SMS reaches GR instability at an age of $\sim 1$ Myr and collapses in $\sim10^4$ s, retaining $\sim 99\%$ all of its mass. Across the full metallicity range from Pop III to $10^{-2}\,Z_\odot$, this shell-ejection channel persists. Pulsational mass-loss from SMSs therefore provides a physically motivated origin for the compact cocoon-like structure implied by LRDs, while remaining the natural progenitors of the massive black hole seeds invoked in direct collapse scenario.

G310.7-5.4 is a supernova remnant (SNR) candidate identified as a faint shell in the second epoch Molonglo Galactic Plane Survey (MGPS-2), but this has not been followed up with multi-wavelength observations until now. It is an example of an SNR at high Galactic latitude showing spatially coinciding $γ$-ray emission. Here, we make the first detailed investigation of the radio emission from the G310.7-5.4 region, aiming to characterise the radio structure, polarisation measurements and the coinciding GeV emission. We used recent radio continuum observations at 943.5 MHz from the EMU and the POSSUM surveys with ASKAP, as well as 16.5 years of Fermi-LAT observations. We furthermore considered the multiwavelength context of the object by investigating observations previously conducted with other instruments, such as infrared and X-ray surveys. We confirm the SNR candidate as a new supernova remnant, dubbed Abeona. We detect the presence of a faint, extended, bilateral radio shell of the size of around 30' diameter and ASKAP radio flux density of $1.5^{+1.5}_{-0.1}$ Jy with no obvious infrared counterparts. With a radio surface brightness of about $2.4^{+2.4}_{-0.1}\times10^{-22}$ W m$^{-2}$ Hz$^{-1}$ sr$^{-1}$, this SNR is one of the faintest radio SNRs known. The northern part of the shell shows linearly polarised radio emission, characteristic of synchrotron emission in SNRs. The physical size of the SNR is estimated to be around $42^{+42}_{-21}$ pc, which would give a distance of around $4.9^{+4.9}_{-2.5}$ kpc. Furthermore, the spatially coincident $γ$-ray source 4FGL J1413.9-6705 shows an energy flux of $1.26\pm0.35\times 10^{-6}$ MeV cm$^{-2}$ s$^{-1}$ with a significance of 5.7 $σ$ between 100 MeV and 100 GeV. The SNR is also put in context with known high-latitude SNRs with $γ$-ray counterparts and compared with their observational properties.

We investigate the scalar sector of linear cosmological perturbations in quadratic gravity. Working in the Einstein frame, we derive the equations of motion in a gauge-independent manner and express them in terms of three sets of gauge-invariant variables. This approach allows us to distinguish genuine physical effects from gauge artefacts, which is particularly relevant for assessing the stability of perturbations in this theory. In the superhorizon limit, we obtain the leading-order behaviour of the relevant gauge-invariant variables and analyse the perturbations in several commonly used gauges. We find that the Newtonian gauge exhibits an apparent instability, characterised by the exponential growth of the metric perturbations. However, this growth is non-generic and gauge-dependent; in the other gauges analysed in this work, the perturbations remain well behaved within the perturbative regime. Physical observables can thus be consistently computed, and the apparent instability is identified as a gauge artefact rather than a pathology of the theory. Our analysis also demonstrates how the evolution behaviour of a gauge-invariant variable changes under the frame transformation and clarifies the relation between results obtained in the Jordan and Einstein frames.

We compute the abundance of cosmological gravitational waves produced during the evolution of an ultralight vector (spin-1) dark matter field. A homogeneous background vector field breaks spatial isotropy, requiring a Bianchi I geometry and inducing a mixing between the scalar, vector, and tensor perturbation sectors. We derive the perturbation equations in this background and show that, as a consequence of this mixing, scalar perturbations act as a source of tensor modes, generating a stochastic GW background. The production and cosmological evolution of these gravitational waves are implemented in a modified version of CLASS, from which we obtain the present-day spectrum.

In this work, we systematically investigate the quantum-information diagnostics of cosmological perturbations with a nontrivial sound speed, utilizing a normalized open two-mode squeezed-state framework. Rather than introducing new observables, our analysis focuses on how a modified sound speed dynamically reshapes the Schrödinger evolution of the squeezing parameters ($r_k$ and $φ_k$). We demonstrate how these dynamical changes are inherited by the reduced density matrix of the observable sector. By employing a sound-speed-resonance parametrization, we derive and evaluate the purity, von Neumann entropy, Rényi entropies, and logarithmic negativity. To overcome the intrinsic multiscale stiffness of the post-inflationary equations, we introduce a bounded variable $x = \tanh r_k$ as a partial regularization, which enables reliable numerical simulations exclusively within the inflationary regime. Our numerical results reveal that a nontrivial sound speed significantly suppresses the purity of the reduced state, indicating enhanced effective mixedness. Simultaneously, it strongly amplifies and modulates both the entropic and entanglement diagnostics. More precisely, a nontrivial sound speed postpones the onset of classicality by modulating the decoherence process. Ultimately, we show that a nontrivial sound speed leaves distinct and identifiable quantum-information signatures within the entanglement structure of the early universe.

Primordial black holes (PBHs) provide a compelling interpretation for the binary black holes (BBHs) observed by ground-based gravitational-wave (GW) detectors, especially for those BBHs in the theoretical mass gap. In the early Universe, the scalar perturbations required to produce such PBHs inevitably generate scalar-induced GWs (SIGWs). These SIGWs peak in the sub-nanohertz band, and manifest secularly as measurable jerk-like drifts in the second derivative of pulsar spin periods. In this Letter, we perform the first search for SIGWs using pulsar parameter drifts, and place a 95\% confidence-level upper limit on the PBH abundance of $f_{\mathrm{PBH}} < 10^{-10}$ over the mass range $[3 \times 10^{-1}, 4 \times 10^{4}] M_{\odot}$. Our results strongly disfavor a PBH origin for the BBHs currently detected by the LIGO-Virgo-KAGRA (LVK) Collaborations.

We model the multi-band lightcurves of 80 SNe Ic-BL, including 11 associated with lGRBs, using a magnetar engine model with $^{56}$Ni decay. We find that the data are all consistent with a magnetar central engine, and such a model yields high-quality fits across the sample. The medians with $1σ$ regions of the key parameters are $P_{\rm{i}}\sim2.04^{+1.84}_{-0.96}\,{\rm{ms}}$, $B_{\rm{p}}\sim3.96^{+3.28}_{-1.40}\times10^{15}\,{\rm{G}}$, $M_{\rm{ej}}\sim2.30^{+1.48}_{-1.02}\,M_\odot$, and $M_{\rm{Ni}}\sim0.18^{+0.14}_{-0.09}\,M_\odot$, with strong and statistically significant correlations observed for both $M_{\rm{ej}}-P_{\rm{i}}$ (anti-correlation) and $M_{\rm{Ni}}-M_{\rm{ej}}$ (correlation). Comparing the SN Ic-BL samples with and without lGRB association using fitting parameters, we find no significant difference between them, although the GRB-associated sample is slightly brighter, possibly due to an observational bias. Relative to ordinary SNe Ic, SNe Ic-BL have similar $^{56}$Ni and ejecta masses, suggesting comparable pre-SN progenitor properties, with differences possibly arising from the presence of a magnetar engine. In comparison with other possible magnetar-powered SESNe, including SLSNe Ic and FBOTs, we confirm a strong universal $M_{\rm{ej}}-P_{\rm{i}}$ correlation, indicating a common origin. SNe Ic-BL and SLSNe Ic have similar ejecta mass distributions, typically $M_{\rm ej}\gtrsim0.5\,M_\odot$, while FBOTs mostly lie below this value. Differences between SNe Ic-BL and SLSNe Ic may arise from magnetar properties, with SN Ic-BL magnetars rotating faster and having stronger fields. Moreover, the $P_{\rm{i}}-B_{\rm{p}}$ distribution of lGRB magnetars largely overlaps with that of SN Ic-BL magnetars. In connection with binary simulation results, we propose a unified physical classification and progenitor framework for magnetar-powered and ordinary SESNe.

We study short-timescale 340 GHz flux-density variability of Sgr A* using ALMA Cycle 3 observations. Careful self-calibration enabled 10 s snapshot imaging with very high effective image-domain SNR, allowing high-cadence monitoring of Galactic Center sources. To reduce atmospheric and instrumental effects, we measured Sgr A* relative to multiple non-variable sources in the same field and corrected apparent variability caused by time-dependent u-v coverage and PSF changes using simulations with a static input model. We then searched for characteristic timescales over 20 s < tau < Tobs/3 using structure functions, the Lomb--Scargle method, and state-space-model autoregressive spectral analysis. No dominant narrow periodicity is found. Instead, the data show a short-timescale flat, white-noise-like regime at tau below about 2.3--6.3 min, followed by red-noise-like behavior at longer timescales. This flat regime appears in both active and quiescent phases, suggesting statistically independent fluctuations on these timescales. We interpret its upper boundary as an empirical transition timescale between decorrelated short-timescale fluctuations and longer-timescale correlated variability. The physical origin of this flat component remains uncertain, since previous theoretical and numerical studies more commonly report red-noise-like or broken-power-law variability.

Long period radio transients (LPTs) are the slowest radio-pulsing sources ever found, with the current population spanning periods of seven minutes to over six hours. Two of the thirteen published LPTs, ILT J1101+5521 and GLEAM-X J0704--37, have been associated with an M dwarf closely orbiting a white dwarf (WD) through optical spectroscopy. Here, we present new Keck I/LRIS optical spectroscopy of ILT J1101+5521, which reveals H$α$ emission from the M dwarf and confirms an orbital period nearly matching the radio period (2.092 hr). Radio pulses in both systems arrive just after maximum M dwarf redshift, assuming the radio period matches the orbital period. Based on Gaia proper motions and systemic velocities, we find that these systems are kinematically hotter and less concentrated in the Galactic plane than other LPTs. Both systems harbor unusually massive and cool WDs, with $M_\mathrm{WD} \approx 0.84-1.0 M_\odot$ and $T_\mathrm{eff} \approx 5200-7300$ K, implying that their carbon-oxygen cores are nearly entirely crystallized. Both systems are unusually close to being face-on binaries ($i=13^\circ-28^\circ$), signaling that the production of coherent radio pulses may be a strongly inclination-dependent phenomenon. We present MESA models that show that the M dwarf in each system will fill its Roche lobe within $\sim1$ Gyr, becoming a cataclysmic variable. Finally, we place lower limits on the space density of WD + M dwarf LPTs ($ρ\gtrsim 10^{-8}\;\mathrm{pc}^{-3}$); based on the broader population of WD + M dwarf binaries, we estimate that there are 100 (2000) WD + M dwarf LPTs within 2 kpc if current radio findings are 100% (10%) complete. Current and upcoming radio surveys will be sensitive to many such systems, and M dwarf optical counterparts out to $\sim$2 kpc will be detectable with the Rubin Observatory Legacy Survey of Space and Time (LSST).

Fermion condensate inflation, where inflation emerges from four-fermion interactions induced by spacetime torsion, removes the need for additional scalar fields beyond the Standard Model. In this framework, the fermion field can be decomposed into two distinguished sectors, each giving rise to bound states. After integrating out fermions, the bound fields play the roles of the inflaton and the auxiliary fields, resembling hybrid inflation with a waterfall mechanism. The inclusion of an axial chemical potential naturally introduces a mechanism to end inflation and trigger instant preheating. During the waterfall phase, the effective potential of the fermion condensate supports the formation of non-topological solitons such as Q-balls, which act as seeds of primordial black holes. This model is intrinsically connected to Chern-Simons gravity, which implies a parity-violating universe. Consequently, both the primordial black hole (PBH) dark-matter abundance and parity-violation signatures could provide observational tests of the model.

High energy neutrinos can be injected in the early Universe from the decay or annihilation of long lived primordial relics. We analyse the possibility that the ultrahigh energy neutrino event recently observed by the KM3NeT neutrino telescope could have such an origin. This possibility has the advantage of leading to a sharp spectral feature in a way that the neutrino flux can be small at all energies except at the KM3NeT event energy. Thus, along this scenario the tension with null results from other experiments is reduced with respect to the usual power law case analysed by the KM3NeT and IceCube experiments. At such energies and for an emission around the recombination time, interactions of these neutrinos with background neutrinos prove to be relevant and must be determined from the development of a dedicated code. These interactions, as well as final state radiation processes, modify the spectrum. Interestingly, it turns out that the scenario can also leave an imprint in the CMB that could be probed in the near future. Interestingly too, this scenario does not predict an associated $γ$-ray flux beyond observation. All in all we do find that the high energy neutrino could be a primordial high energy neutrino, provided it has been produced around the recombination time or later.

The multi-messenger exploration of dark matter and physics beyond the Standard Model has emerged as a central direction in modern astro-particle physics, particularly following the discovery of gravitational waves. In this work, we present a comprehensive review and forward-looking perspective on machine-learning-enhanced multi-messenger approaches, combining information from gravitational waves, cosmic rays, gamma rays, neutrinos, and collider experiments. We summarize the current state of the field, discuss recent methodological developments, and outline a coherent research program aimed at integrating heterogeneous datasets within a unified inference framework. Our collaboration proposes here a plan for forthcoming analyses aiming at extracting information on the properties and interactions of dark matter, and finally on its genesis, combining multi-messenger astronomy techniques and inputs from laboratory physics. The main objectives planned in this line of research comprise: i) the multi-messenger analysis of new physics in cosmology, including mainly, but not only, several different models of dark matter; ii) the phenomenology of new physics signatures in ground-based cosmic rays experiments, with cross-correlation to the corresponding physical, astrophysical and cosmological observations; iii) the development of machine learning methods for data analysis in ground-based cosmic rays experiments, in light of the new physics signatures. We note that several groups have explored the use of multi-messenger observations, including gravitational waves, to probe alternative dark matter candidates. The present work builds on these developments by focusing on the role of machine learning in integrating heterogeneous datasets. We foresee that such a cross-fertilizing approach will represent the right path to extract information about the main questions left in fundamental physics.

Cosmic rays (CRs) drive ionization and influence gas dynamics in molecular clouds (MCs), potentially impacting the resulting star formation outcomes. Although previous simulations of individual star formation have included methods for cosmic ray transport (CRT), none have been large enough to resolve the stellar initial mass function (IMF). We conduct numerical simulations following the collapse of a $20000 M_{\odot}$ MC and the subsequent star formation including CRT, both with and without CRs accelerated by winds from the young massive stars, and compare against a non-CRT simulation. We show that after the first massive stars form, the cavity produced by feedback is more pronounced in the CRT simulations because the external CRs are able to propagate inwards and compress the gas into higher density structures. This increases the subsequent star formation in the cloud; by the end of the simulation, the SFE in the CRT simulation including stellar wind CRs is 43 \% higher than the non-CRT simulation. The IMF is also top heavy in comparison, with a slope above 1 $M_{\odot}$ that is shallower by $\sim 20$ \%. These effects are also present in the simulation without wind-accelerated CRs, but they are not as pronounced; the SFE is only 16 \% higher than the non-CRT simulation, and the IMF high-mass slope is shallower by $\sim 10$ \%. These results may explain some of the observed top-heavy IMFs, which typically occur in high-CR environments such as the galactic center.

We investigate timelike geodesics in asymptotically flat regular black holes supported by a phantom scalar field characterized by a scalar charge $A$. This parameter removes the central singularity and continuously deforms the Schwarzschild geometry while preserving asymptotic flatness. We derive the equations of motion for massive test particles and classify bounded and unbounded trajectories in terms of the conserved energy and angular momentum. We determine circular and critical orbits, including the innermost stable circular orbit (ISCO), and analyze the transition between capture and scattering. We show that the scalar charge modifies the location of the unstable and stable circular orbits, the ISCO, and the threshold angular momentum for scattering, exhibiting a nontrivial dependence on the radial coordinate. Their physical scales are naturally described in terms of the invariant areal radius $R(r)=\sqrt{r^2+A^2}$. In the weak-field regime, we compute the perihelion precession and obtain corrections proportional to the scalar charge, allowing us to constrain the scalar charge from Solar System observations. We also analyze the motion with vanishing angular momentum and show that, while the qualitative structure of the trajectories remains connected to the Schwarzschild limit $A\to 0$, the quantitative deviations encode the geometric effects of the scalar hair.