Vote on papers for Monday, Oct 20 2025

The intracluster medium (ICM) records the history of galaxy clusters through its complex dynamical properties. To effectively interpret these properties, robust methods are needed to compare observational data with theoretical models. We present a novel machine learning framework for comparing ICM line-of-sight velocity maps derived from X-ray observations. Our approach uses convolutional and Siamese neural networks to identify similarities between different kinematic fields. We outline the architecture of this framework and perform a series of sanity checks to validate its performance. These checks demonstrate the model's ability to correctly identify and quantify kinematic features, establishing a powerful new tool for future comparative studies of the ICM.

Recent JWST/NIRSpec observations have revealed high-$z$ star-forming galaxies depart from the Fundamental Metallicity Relation (FMR), yet the $z = 0$ FMR has not been well-characterized in the low-mass regime ($\rm log(M_{\star}/M_{\odot}) \lesssim 9$) for an appropriate comparison of low- and high-$z$ systems. We attempt to rectify this limitation through a meta-analysis, providing a local, observational comparison for future high-$z$ FMR studies. We analyzed common FMR fitting methods for $\sim 700$ [OIII]$\lambda 4363$ emitters with $\rm log(M_{\star}/M_{\odot}) \lesssim 9$ at $z \sim 0$. We find no evidence of the FMR below $\rm log(M_{\star}/M_{\odot}) \lesssim 9$ through any method, suggesting that slowly-evolving, quasi-steady state gas reservoirs are not yet established. We simultaneously find a weak positive correlation between metallicity and star formation, and that these systems are gas-rich with substantial diversity in effective yields ($y_{\rm eff}$) spanning $\rm \sim 3~dex$. We demonstrate increasing $y_{\rm eff}$ correlates with decreasing FMR offsets, which in the context of the analytical and non-equilibrium gas models of Dalcanton et al. (2007), indicates a scenario where star formation bursts rapidly return and eject metals from the ISM before subsequent gas-balancing. Pristine infall diluting the ISM metal-content cannot lead to the $y_{\rm eff}$ diversity we measure, and thus is not the primary process behind FMR deviations. Our results suggest low-$\rm M_{\star}$ systems, regardless of redshift, depart from a steady-state gas reservoir shaping the canonical FMR, in which metallicity variations are primarily driven by star formation and enriched outflows. With this characterization, we demonstrate $z \gtrsim 3$ [OIII]$\lambda 4363$ systems are indeed more metal-poor than $z \sim 0$ counterparts ($\rm \Delta 12+log(O/H) = 0.3~dex$) at fixed $\rm M_{\star}$.

Black holes (BHs) play a major role in the structural and dynamical evolutions of Globular Clusters (GCs). Several recent works searched for BHs in Galactic GCs using scaling relations derived from numerical simulations. However, the conclusions drawn by such approaches are strongly dependent on the specific prescriptions adopted in numerical simulations. Therefore, we analyzed a survey of 101 Monte Carlo simulations to identify some observable parameters that allow us to probe the present-day BH population in GCs reliably. We thoroughly show that a single observable is not suited to infer the BH mass fraction in real GCs: similar values could be attained by systems with different BH mass fractions, depending on the specific dynamical evolution of the system. Finally, we present a combination of observable parameters that efficiently breaks this degeneracy. We also compare values from numerical simulations with a sample of Galactic GCs.

We present numerical relativity simulations of binary neutron star mergers incorporating neutrino flavor transformations triggered by fast flavor instability, quantum many-body effects, or potential beyond standard model physics. In both long-lived and short-lived remnant scenarios, neutrino flavor conversions modify species-dependent neutrino luminosities and mean energies, and drive the matter towards more neutron rich conditions. They produce up to $300\%$ more neutron rich ejecta and significantly boost the r-process yields, especially in low-density, near-equatorial outflows. We identify regions unstable to fast flavor instabilities and find that these instabilities persist despite flavor conversions. We further test the sensitivity to the equilibration timescale of the flavor conversions, finding that slower flavor conversions can interact with thermodynamic equilibration, and increase the neutron richness of the ejecta. Flavor conversions may also contribute to stronger gravitational wave and neutrino emissions, pointing to a correlation between neutrino transport and merger dynamics. These results highlight the potential impact of flavor conversions while motivating future work to improve on theoretical understanding of flavor instabilities in global simulations.

We investigate the orbital eccentricity evolution of supermassive black hole binaries within galactic environments. We analyze the dynamics in triaxial merger remnants and subsequent interactions with geometrically thick nuclear discs. We confirm that gravitational torques in triaxial potentials efficiently extract angular momentum, resulting in binary formation with high initial eccentricities. We then analyze the binary-disc interaction using a 3D analytical framework incorporating the Airy formalism and potential softening. We present a self-consistent derivation demonstrating that the 3D suppression of high-order torques leads to distinct scalings with disc thickness ($h$): migration rates $\tau_a^{-1} \propto h^{-3}$ and eccentricity damping rates $\tau_e^{-1} \propto h^{-5}$. This establishes a timescale hierarchy, $\tau_e/\tau_a \propto h^2$. For typical parameters ($h\approx 0.2$), eccentricity damping is significantly faster than orbital decay ($\tau_e \approx 0.04 \, \tau_a$). We further develop a wavelet-based formalism to quantify the impact of disc inhomogeneities arising from accretion feedback and turbulence. We derive the stochastic torque variance in the wavelet domain and employ a Fokker-Planck analysis to determine the equilibrium eccentricity distribution. While stochastic fluctuations counteract deterministic damping, the strong damping imposed by the thick disc geometry ensures the equilibrium eccentricity remains small unless the fluctuations are highly non-linear. Hence, even if born highly eccentric, SMBHBs are rapidly circularized. This circularization trap forces binaries to approach the gravitational wave-dominated regime on nearly circular orbits, prolonging the total merger timescale. This introduces a substantial cosmological delay governed by stellar relaxation, which impacts detection rates and the modeling of SMBH assembly in cosmological frameworks.

Globular clusters (GCs) host multiple stellar populations differing in their chemical and dynamical properties. A number of models for the formation of multiple populations predict that the subsystem of second generation (SG) stars is characterized by a more centrally concentrated spatial distribution and a more rapid rotation than the system of first generation (FG) stars. We present the results of N-body simulations exploring the long-term dynamical evolution of rotating multiple-population GCs. We study the evolution of systems starting with four different orientations of the GC's total internal angular momentum vector relative to the orbital angular momentum. We explore the evolution driven by two-body relaxation and the effects of the GC's interaction with the galactic tidal field. We focus on the kinematic differences between the two generations and we quantify them by exploring the FG and SG rotation velocity and angular momenta. We find that kinematic differences between the generations persist for most of the GCs' lifetimes, although the strength of these differences decreases after a few relaxation times. The differences can be seen most clearly in the lowest-mass stars. We find that the GCs' internal angular momentum gradually aligns with the orbital angular momentum, although there is little difference in this alignment between the FG and SG systems. We also find that stars in the GC's outer regions align with the orbital angular momentum vector more rapidly than those in the inner regions leading to a variation of the orientation of the internal angular momentum with the clustercentric distance. The alignment between internal angular momentum and orbital angular momentum occurs more rapidly for low-mass stars. We study the evolution of the anisotropy in the velocity distribution and find the SG to be characterized by a stronger radial anisotropy than the FG.(abridged)

We present a general analytic framework for the evolution of cosmic structure in multi-species dark matter models that simultaneously incorporates finite velocity dispersion and Poisson fluctuations. Our approach accommodates arbitrary numbers of dark matter components with distinct mass fractions, velocity distributions, and number densities -- ranging from cold particles to warm species and sparse populations such as primordial black holes or solitons. The framework is based on solving a truncated BBGKY hierarchy, whose solution is obtained by solving Volterra integral equations. We provide an efficient algorithm to solve for the total, as well as inter- and intra-species power spectra. Worked examples with two-component mixtures illustrate how isocurvature (initially Poisson) and adiabatic spectra evolve differently depending on the properties of the warm or sparse fraction. This evolution is controlled by the free-streaming and Jeans scales, and the results match analytic estimates and $N$-body simulations.

We introduce this http URL, an open-source and publicly available general relativistic ray-tracing toolkit for spectral modelling in arbitrary spacetimes. Our software is written in the Julia programming language, making use of forward-mode automatic differentiation for computing the Christoffel symbols during geodesic integration, and for propagating derivatives through the entire ray-tracer. Relevant numerical methods are detailed, and our models are validated using a number of tests and comparisons to other codes. The differentiability is used to optimally calculate Cunningham transfer functions -- used to efficiently pre-compute relativistic effects in spectral models. A method is described for calculating such transfer functions for disc with non-zero vertical height, including the treatment of self-obscuration. An extension of the transfer function formalism that includes timing information is described, and used to calculate high-resolution reverberation lag spectra for a lamppost corona. The lag-frequency and lag-energy spectra for a Shakura-Sunyaev accretion disc with various lamppost heights and Eddington ratios are calculated, and the general impact of disc thickness in reflection models is discussed.

A recent determination of the growth index indicates a value significantly higher than the $\Lambda$CDM prediction, suggesting that alternative scenarios to $\Lambda$CDM may be required. In this work, we investigate whether a time-varying Newton's constant, $G_N$, can account for such a high growth index, $\gamma=0.063\pm0.025$. Adopting a phenomenological approach, we study two parameterizations of the effective gravitational coupling, $G_{\rm eff}$, one based on a Taylor expansion and another linked to the energy density parameter of Dark Energy. We constrain the models with Cosmic Chronometers (CC), Dark Energy Spectroscopic Instrument baryon acoustic oscillations (data release 2), CMB priors, and a gaussian likelihood for the growth index. We show that the constant $\gamma$ approximation is accurate for the parametrization linked to the energy density parameter of dark energy, but presents a non-negligible error for the other case, which we treat as a systematic error in the analysis. We find a $2.4\sigma-3.4\sigma$ tension level with constant $G_{\rm eff}$, depending on the parametrization. The results indicate that $G_{\rm eff}<G_N$ around the period of accelerated expansion, corresponding to a weaker effective gravitational interaction on cosmological scales, which leads to a suppression of the growth of cosmological structures.

Disc accretion onto astrophysical objects with a material surface proceeds through the boundary layer (BL) -- a radially narrow region in the inner disc where the incoming gas must slow down its rotation before settling onto the surface of the accretor. Here we numerically study a BL in which the angular momentum transport in the layer is accomplished via the excitation of global acoustic waves. While the earlier studies of such wave-mediated BLs typically modeled the internal structure of the central object as a globally isothermal sphere with sharply rising density profile, here we explore the effect of other internal density and temperature profiles on the mode operation. We model the inner structure of an accretor as a polytropic sphere, allowing a shallower increase of density and a non-trivial temperature profile inside the object. While the mix of acoustic modes observed in our long-duration (1000 inner orbits long) 2D hydrodynamic simulations is a weak function of the polytropic index $n$ of the accretor's structure, the mass accretion rate and the angular momentum flux across the BL show a clear dependence on $n$, both decreasing in amplitude as $n$ is lowered. Interestingly, in 2D these transport metrics are better correlated not with $n$ but with a total mass inside the central object contained within the simulation domain. These results improve our understanding of the wave-mediated BL accretion by quantifying the effect of the inner structure of the accretor on the excitation and propagation of acoustic modes mediating the BL transport.

Super-Earths and mini-Neptunes are missing from our Solar System, yet they appear to be the most abundant planetary types in our Galaxy. A detailed characterization of key planets within this population is important for understanding the formation mechanisms of rocky and gas giant planets and the diversity of planetary interior structures. In 2019, NASA's TESS satellite found a transiting planet candidate in a 17.6-day orbit around the star TOI-283. We started radial velocity (RV) follow-up observations with ESPRESSO to obtain a mass measurement. Mass and radius are measurements critical for planetary classification and internal composition modeling. We used ESPRESSO spectra to derive the stellar parameters of the planet candidate host star TOI-283. We then performed a joint analysis of the photometric and RV data of this star, using Gaussian processes to model the systematic noise present in both datasets. We find that the host is a bright K-type star ($d = 82.4$ pc, $\mathrm{T}_\mathrm{eff} = 5213 \pm 70$ K, $V = 10.4$ mag) with a mass and radius of $\mathrm{M}_\star = 0.80 \pm 0.01\; \mathrm{M}_\odot$ and $\mathrm{R}_\star = 0.85 \pm 0.03\; \mathrm{R}_\odot$. The planet has an orbital period of $P = 17.617$ days, a size of $\mathrm{R}_\mathrm{p} = 2.34 \pm 0.09\; \mathrm{R}_\oplus$, and a mass of $\mathrm{M}_\mathrm{p} = 6.54 \pm 2.04\; \mathrm{M}_\oplus$. With an equilibrium temperature of $\sim$600 K and a bulk density of $\rho_\mathrm{p} = 2.81 \pm 0.93$ g cm$^{-3}$, this planet is positioned in the mass-radius diagram where planetary models predict H$_2$O- and H/He-rich envelopes. The ESPRESSO RV data also reveal a long-term trend that is probably related to the star's activity cycle. Further RV observations are required to confirm whether this signal originates from stellar activity or another planetary body in the system.

We present an atlas of full-scan spectra from the Short-Wavelength Spectrometer (SWS) aboard the Infrared Space Observatory (ISO) after reprocessing and improving an earlier version published 22 years ago. The SWS spectra cover the wavelength range from 2.35 to 45.3 {\mu}m. They include scans in 12 separate bands, and we have updated the methods used to combine those bands into a single continuous spectrum. The main improvement comes from applying multiple constraints, including new photometry and spectra from the Infrared Spectrograph (IRS) on the Spitzer Space Telescope that have become available since the release of the original products, and individualized attention to each spectrum, to renormalize the separate bands into a more consistent single spectrum. In particular this removed unphysical negative fluxes that were common in the original data products. The new database, with 1035 reprocessed spectra, will be available to the community at IRSA, which also hosts the original processing.

Gravitational waves from core-collapse supernovae provide a unique probe of the equation of state (EOS) of high density matter. In this work, we focus on the bounce signal from numerical simulations of rotating supernovae and explore its potential for EOS inference. We employ a support vector machine, previously shown to perform best among tested methods, to classify GW signals simulated for 18 EOS models. For optimally oriented sources, we estimate that the Advanced LIGO A+ detector can probe the EOS for Galactic events, while third-generation observatories such as the Einstein Telescope and Cosmic Explorer can reach substantially farther. For randomly oriented sources, only these next-generation detectors are expected to have sufficient sensitivity. These results represent the potential observational range for probing the nuclear EOS, although, due to the simplifying assumptions adopted, they should be regarded as approximate upper limits.

Direction of arrival (DOA) estimation is mostly performed using specialized arrays that have carefully designed receiver spacing and layouts to match the operating frequency range. In contrast, radio interferometric arrays are designed to optimally sample the Fourier space data for making high quality images of the sky. Therefore, using existing radio interferometric arrays (with arbitrary geometry and wide frequency variation) for DOA estimation is practically infeasible except by using images made by such interferometers. In this paper, we focus on low cost DOA estimation without imaging, using a subset of a radio interferometric array, using a fraction of the data collected by the full array, and, enabling early determination of DOAs. The proposed method is suitable for transient and low duty cycle source detection. Moreover, the proposed method is an ideal follow-up step to online radio frequency interference (RFI) mitigation, enabling the early estimation of the DOA of the detected RFI.

While massive quiescent galaxies are known to be poor in atomic hydrogen (HI), their true HI content remains poorly constrained due to the limited sensitivity and morphological biases of existing surveys. We present deep HI observations using the Five-hundred-meter Aperture Spherical radio Telescope (FAST) for a representative sample of 78 low-redshift massive quiescent galaxies, selected by stellar mass ($M_\ast > 10^{10} M_\odot$) and color (NUV$-r > 5$). Our observations reach an exceptional detection limit of $\log(M_{HI}/M_\ast) = -3$. We find that one-third of the sample follows the predicted HI mass fraction distribution based on previous surveys. However, the majority ($\sim$2/3) exhibit significantly lower HI content, with a $3\sigma$ upper limit of $\log(M_{HI}/M_\ast) < -3.46$ derived from stacking the non-detections. As a consequence of this enormous dynamic range and the high fraction of non-detections, the HI mass fraction shows no clear correlation with parameters tracing star formation, structure, or morphology. Our FAST sample shows remarkable similarity to the $ATLAS^{3D}$ sample which only includes early-type galaxies, both in its high fraction of HI-poor galaxies and its high satellite fraction among HI-poor galaxies, as indicated by three different environmental characteristics: projected cross-correlation functions, background-subtracted neighbour counts, and central/satellite classification in the SDSS galaxy group catalog. These results suggest that while early-type morphology and environment play crucial roles in suppressing HI reservoirs, neither factor alone fully explains the observed gas depletion, indicating that additional physical mechanisms must be responsible for the extreme HI deficiency prevalent in massive quiescent galaxies.

Some compact stars may contain deconfined quark matter, forming hybrid stars or quark stars. If the quark matter forms an inhomogeneous condensate in the crystalline color superconducting phase, its rigidity may be high enough to noticeably alter the stellar properties. In this paper, we investigate whether these elastic stars follow the universal relations, i.e., relations insensitive to equations of state, that have been well established for fluid stars. We improve upon previous studies by allowing quark matter in the background, static, and spherically symmetric configuration to be sheared. Such background shear can be treated in the form of an effective pressure anisotropy. We then calculate the moment of inertia $I$, tidal deformability ${\lambda}_2$, and spin-induced quadrupole moment $Q$ of these models with pressure anisotropy. The $I$-${\lambda}_2$-$Q$ universal relations for the elastic hybrid (quark) star models are valid up to a variation of $\approx2\,(3)\%$, larger than that for typical fluid star models, when the maximal magnitude of quark matter shear modulus is considered in the crystalline color superconducting phase from realistic calculations. The uncertainty in universal relations related to the stellar compactness for these elastic star models, on the other hand, remain comparable to those for typical fluid star models. Our results demonstrate the validity of universal relations for hybrid stars and quark stars with a realistic degree of pressure anisotropy due to the crystalline color superconducting quark matter.

The next-generation radio-interferometric (RI) telescopes require imaging algorithms capable of forming high-resolution high-dynamic-range images from large data volumes spanning wide frequency bands. Recently, AIRI, a plug-and-play (PnP) approach taking the forward-backward algorithmic structure (FB), has demonstrated state-of-the-art performance in monochromatic RI imaging by alternating a data-fidelity step with a regularisation step via learned denoisers. In this work, we introduce HyperAIRI, its hyperspectral extension, underpinned by learned hyperspectral denoisers enforcing a power-law spectral model. For each spectral channel, the HyperAIRI denoiser takes as input its current image estimate, alongside estimates of its two immediate neighbouring channels and the spectral index map, and provides as output its associated denoised image. To ensure convergence of HyperAIRI, the denoisers are trained with a Jacobian regularisation enforcing non-expansiveness. To accommodate varying dynamic ranges, we assemble a shelf of pre-trained denoisers, each tailored to a specific dynamic range. At each HyperAIRI iteration, the spectral channels of the target image cube are updated in parallel using dynamic-range-matched denoisers from the pre-trained shelf. The denoisers are also endowed with a spatial image faceting functionality, enabling scalability to varied image sizes. Additionally, we formally introduce Hyper-uSARA, a variant of the optimisation-based algorithm HyperSARA, promoting joint sparsity across spectral channels via the l2,1-norm, also adopting FB. We evaluate HyperAIRI's performance on simulated and real observations. We showcase its superior performance compared to its optimisation-based counterpart Hyper-uSARA, CLEAN's hyperspectral variant in WSClean, and the monochromatic imaging algorithms AIRI and uSARA.

Precise frequency values have been determined for H$_2^{~16}$O radio lines appearing in protected line lists of the International Astronomical Union and the Panel on Frequency Allocations of the US National Academy of Sciences. The improved precision is attributable to a spectroscopic network built from a large set of near-infrared Lamb-dip lines augmented with a handful of ultrahigh-accuracy rotational transitions. The ultraprecise H$_2^{~16}$O network contains 376 Lamb-dip lines recorded previously via our frequency-comb locked cavity-enhanced spectrometer. During the present study, altogether 55 target lines have been (re)measured at high accuracy. Due to our network-assisted measurements, the accuracy has been significantly improved with respect to previous direct radio-frequency measurements for all the protected lines of H$_2^{~16}$O above 750 GHz. Furthermore, 43 of these protected transitions now benefit from the accuracy of the new near-infrared Lamb dips reported in this paper.

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will commence full-scale operations in 2026, yielding an unprecedented volume of astronomical images. Constructing an astronomical catalog, a table of imaged stars, galaxies, and their properties, is a fundamental step in most scientific workflows based on astronomical image data. Traditional deterministic cataloging methods lack statistical coherence as cataloging is an ill-posed problem, while existing probabilistic approaches suffer from computational inefficiency, inaccuracy, or the inability to perform inference with multiband coadded images, the primary output format for LSST images. In this article, we explore a recently developed Bayesian inference method called neural posterior estimation (NPE) as an approach to cataloging. NPE leverages deep learning to achieve both computational efficiency and high accuracy. When evaluated on the DC2 Simulated Sky Survey -- a highly realistic synthetic dataset designed to mimic LSST data -- NPE systematically outperforms the standard LSST pipeline in light source detection, flux measurement, star/galaxy classification, and galaxy shape measurement. Additionally, NPE provides well-calibrated posterior approximations. These promising results, obtained using simulated data, illustrate the potential of NPE in the absence of model misspecification. Although some degree of model misspecification is inevitable in the application of NPE to real LSST images, there are a variety of strategies to mitigate its effects.

The correlation between optical starlight polarization and polarized thermal dust emission can be used to infer intrinsic dust properties. This correlation is quantified by the ratio Rp/p, which has been measured to be 5.42 +/- 0.05 MJy sr^-1 at 353 GHz when averaged over large areas of the sky. We investigate this correlation using newly published stellar polarimetric data densely sampling a continuous sky region of ~4 square degrees at intermediate Galactic latitude. We combine RoboPol optical polarization measurements for 1,430 stars with submillimeter data from the Planck satellite at 353 GHz. We perform linear fits between the Planck (Qs, Us) and optical (qv, uv) Stokes parameters, accounting for the differences in resolution between the two datasets as well as the distribution of clouds along the line of sight. We find that in this region of the sky the Rp/p value is 3.67 +/- 0.05 MJy sr^-1, indicating a significantly shallower slope than that found previously using different stellar samples. We also find significant differences in the fitted slopes when fitting the Qs-qv and Us-uv data separately. We explore two explanations using mock data: miscalibration of polarization angle and variations in Rp/p along the line of sight due to multiple clouds. We show that the former can produce differences in the correlations of Qs-qv and Us-uv, but large miscalibration angles would be needed to reproduce the magnitude of the observed differences. Our simulations favor the interpretation that Rp/p differs between the two dominant clouds that overlap on the sky in this region. The difference in Rp/p suggests that the two clouds may have distinct dust polarimetric properties. With knowledge from the tomographic decomposition of the stellar polarization, we find that one cloud appears to dominate the correlation of Us-uv, while both clouds contribute to the correlation of the Qs-qv data.

Multiplicity is ubiquitous among massive stars. While the stellar components usually display similar masses, some binaries with extremely low mass ratios were also observed. Some of them are primordial, while others arise from binary interactions. The identification of systems with extreme mass ratios brings valuable information, notably on the origin of fast rotation in massive stars. We identify new short-period systems with extreme mass ratios through the detection of eclipses and reflection effects. The physical properties of a dozen newly identified cases were precisely evaluated through high-quality photometry and spectroscopy. In addition to characterizing these binaries, we found a clear signature of apsidal motion in one system, and three other systems display long-term shifts in eclipse times. All systems we reported here are composed of a massive star and a cool low-mass companion. They are therefore primordial cases. This doubles the known number of these systems in the Galaxy. In this context, it is important to note that most massive stars in these systems, as well as in previous systems reported in the literature, rotate fast (supersynchronous compared to the orbital motion). The high incidence of fast rotation in these nascent binaries provides strong constraints for star formation models.

Type I X-ray bursts (XRBs) are thermonuclear runaways on the surface of accreting neutron stars, powered by rapid proton-capture and alpha-capture processes on neutron-deficient nuclei. Uncertainties in the corresponding reaction rates remain a major limitation in modeling burst light curves and ashes. We present a systematic study of the sensitivity of XRB models to uncertainties in charged-particle-induced reaction rates across a broad parameter space of accretion rates and fuel compositions in low-mass X-ray binaries. The study proceeds in two stages: ignition conditions are first determined with a semi-analytic framework coupled to a full reaction network, followed by a sensitivity analysis using the ONEZONE model with individual rate variations. We identify 41 reactions that alter the burst light curve and 187 that significantly impact final abundances. Reactions on bottleneck isotopes in the alpha-p- and rp-process paths strongly affect both observables, while most (p, gamma) reactions on medium-mass (A > 32) and heavy-mass (A > 55) nuclei influence only the final composition. Medium-mass cases dominate in He-rich bursts, where the reaction flow terminates earlier, while heavy-mass cases appear in mixed H and He bursts with extended rp-process paths reaching A ~ 110. We identify a subset of reactions whose rate uncertainties exert influence on the final 12C yield in helium-rich bursts, which could have important consequences for the mechanism of ignition of carbon superbursts. Our results identify key targets for nuclear reaction experiments to reduce nuclear physics uncertainties in XRB models.

Stellar flares are intense bursts of radiation caused by magnetic reconnection events on active stars. They are especially frequent on M dwarfs, where they can strongly influence planetary habitability. Flare frequency distributions (FFDs) are usually modeled as power laws, but recent studies have proposed alternatives such as lognormal distributions, implying different flare mechanisms and planetary impacts. This work investigates which statistical distribution best describes flare occurrences on M dwarfs, considering both equivalent duration (ED), the quantity directly measured from photometry, and bolometric energy, which is more relevant for habitability assesments. We analyzed 110 M dwarfs observed with TESS and CHEOPS, detecting 5,620 flares. Complex events were decomposed, detection biases corrected, and FFDs from both missions scaled to build a combined distribution spanning nearly 10 orders of magnitude in bolometric energy. ED-based FFDs follow a power law, reflecting intrinsic photometric flare occurrence. However, bolometric energy-based FFDs deviate from a pure power law, being better described by a lognormal distribution or, more accurately, by a truncated power law with a break near $10^{33}$ erg, the typical superflare threshold. This truncation suggests a change in flare-generation physics between regular flares and superflares, with implications for the cumulative impact on exoplanetary atmospheres. The apparent low-energy flattening previously attributed to lognormal behavior arises from observational biases, while the drop in flare frequency above $10^{35}$ erg remains unexplained, possibly reflecting an intrinsic cutoff or current observational limits. The upcoming PLATO mission will be well suited to probe both regimes.

Imaging atmospheric Cherenkov telescopes (IACTs) are the main technique for detecting gamma rays with energies between tens of GeV and hundreds of TeV. Amongst them, the High Energy Stereoscopic System (H.E.S.S.) has pioneered the use of different telescope types to achieve an energy range as broad as possible. A large, 28 m diameter telescope is used in monoscopic mode to access the lowest energies ($E \gtrsim 30$ GeV), while the four smaller, 12 m diameter telescopes are used in stereoscopic mode to study energies between 150 GeV and 100 TeV. Nevertheless, a combination of both telescope types and trigger strategies has proven to be challenging. In this work, we propose for the first time an analysis based on event types capable of exploiting both telescope types, trigger strategies, and the whole energy range of the experiment. Due to the large differences between monoscopic and stereoscopic reconstructions, the types are defined based on Hillas parameters of individual events, resulting in three types (Type M, Type B, and Type A), each dominating over a different energy range. The performances of the new analysis configurations are compared to the standard configurations in the H.E.S.S. Analysis Package (HAP), Mono and Stereo. The proposed analysis provides optimal sensitivity over the whole energy range, in contrast to Mono and Stereo, which focus on smaller energy ranges. On top of that, improvements in sensitivity of 25-45% are found for most of the energy range. The analysis is validated using real data from the Crab Nebula, showing the application to data of an IACT analysis capable of combining significantly different telescope types with significantly different energy ranges. Larger energy coverage, lower energy threshold, smaller statistical uncertainty, and more robustness are observed. The need for a run-by-run correction for the observation conditions is also highlighted.

We study inflation driven by a scalar potential arising from composite-sector dynamics, inspired by generalized composite Higgs models. The introduction of a non-minimal coupling, possessing the same functional form as the potential, induces a flattening at large field values that enables successful inflation. We analyze the conditions for ultra-slow-roll inflation, which leads to enhanced curvature perturbations, by combining analytical criteria near the inflection point with comprehensive numerical scans of the parameter space. The region consistent with Cosmic Microwave Background constraints and yielding approximately $N_e \approx 55\text{-}60$ e-folds also predicts primordial black holes with masses in the range $10^3\text{-}10^5\,\mathrm{g}$. Although such ultra-light primordial black holes are typically expected to have evaporated, recent proposals invoking evaporation suppression via memory-burden effects could allow their survival as viable dark matter candidates. Under this assumption, the induced gravitational wave signal peaks in a frequency gap between standard interferometers and resonant-cavity experiments, motivating exploration of this observational frontier.

Like light, gravitational waves are gravitationally lensed by intervening massive astrophysical objects, such as galaxies, clusters, black holes, and stars, resulting in a variety of potentially observable gravitational-wave lensing signatures. Searches for gravitational-wave lensing by the LIGO-Virgo-KAGRA (LVK) collaboration have begun. One common method focuses on strong gravitational-wave lensing, which produces multiple "images": repeated copies of the same gravitational wave that differ only in amplitude, arrival time, and overall "Morse phase." The literature identifies two separate approaches to identifying such repeated gravitational-wave events based on frequentist and Bayesian approaches. Several works have discussed selection effects and identified challenges similar to the well-known "birthday problem", namely, the rapidly increasing likelihood of false alarms in an ever-growing catalogue of event pairs. Here, we discuss these problems, unify the different approaches in Bayesian language, and derive the posterior odds for strong lensing. In particular, the Bayes factor and prior odds are sensitive to the number of gravitational-wave events in the data, but the posterior odds are insensitive to it once strong lensing time delays are accounted for. We confirm Lo et al.'s (2020) finding that selection effects enter the Bayes factor as an overall normalisation constant. However, this factor cancels out in the posterior odds and does not affect frequentist approaches to strong lensing detection.

X-ray observations are essential to achieve a deeper understanding of the broadband emission mechanism in blazars. Here, we present a long-term spectral and temporal analysis of X-ray and optical observations of 1E 0229+200 collected with the Neil Gehrels Swift Observatory from 2008 to 2024, complemented by hard X-ray observations from the Nuclear Spectroscopic Telescope Array NuSTAR. The blazar 1ES 0229+200 is a high-frequency, peaked BL Lac object, known for its exceptionally hard very high-energy (VHE) $\gamma$-ray spectrum extending up to 10 TeV. In August 2021, NuSTAR observed the source in a low X-ray state, revealing a concave spectral shape with a distinct upturn around 25 keV. This feature contrasts with previous observations performed with NuSTAR and Swift-BAT, which showed no such spectral upturn. Previous observations of 1ES 0229+200 and broadband SED (spectral energy distribution) modelling suggest that its X-ray emission extends beyond 100 keV without a significant cutoff. The newly detected spectral upturn may indicate a transition between the synchrotron and inverse Compton components or could be linked to photohadronic processes involving high-energy neutrinos. We discuss the implications of this finding in the context of blazar spectral energy distributions, particularly the potential existence of a third SED bump in the kiloelectronvolt to megaelectronvolt range. The observed spectral features support the hypothesis that 1ES 0229+200 could be a source of high-energy neutrino emission.

The knowledge of the Venus near-surface atmosphere is sparse. Few spacecrafts landed on the surface and measured winds with amplitudes below 1 m/s. The diurnal cycle of the wind amplitude and orientation is not known. Recent numerical simulations showed that slope winds along topographic structures could strongly impact the direction of winds. This study presents the first mesoscale modelling of such winds on Venus. A change of direction is occurring during the day in the main slopes, with upslope winds at noon due to solar heating and downslope winds at night. This is due to efficient IR cooling of the surface during the night, being colder than its surroundings slope atmospheric environment and leading to displacement of air. The temperature is impacted by the adiabatic cooling/warming induced by those winds. A strong heating effect is occurring for the downslope winds, leading to an anti-correlation between the surface temperature diurnal amplitude and the topography. This diurnal amplitude reaches 4 K in the plains and below 1 K in the mountains. The saltation of sediment by those winds was also quantified, with a higher probability at night along the slopes on the western flanks.

All future lunar missions require a definition of the lunar reference system and a realization in the form of the lunar reference frame to ensure consistent products for positioning, navigation, cartography, and timing. This paper defines the origin, orientation, and scale of the Lunar Reference System (LRS), as well as provides numerical solutions for the first realization of the International Lunar Reference Frame (ILRF). ILRF is defined as the Principal Axis (PA) system, attached to the surface and co-rotating with the Moon, with its origin in the lunar center of mass (lunocenter). The ILRF realization is based on variance component estimation of the three lunar ephemeris solutions: INPOP21a, DE430, and EPM2021 for the series of the position of the lunar center of mass and rotation Euler angles -- precession, nutation, and proper rotation. The solution is valid starting with the period covered by Lunar Laser Ranging (LLR) data in 1970 and ending with extrapolated ILRF realizations in 2052 for future lunar missions. Results. The combined ILRF is characterized by the mean error of 17.6 cm for 2010-2030, where 15.3 cm comes from the origin and 8.6 cm from the orientation realization. The error in the realization of the origin is mainly caused by a poor geometry of the retroreflector network, resulting in a high correlation between the scale and the X component of the lunocenter in PA. The LLR post-fit residuals in ILRF are at the level of 2-3 cm in terms of the standard deviations of one-way ranges for best-performing LLR stations. The mean errors of the transformation between ILRF and other reference frame realizations in PA are at the level of 3 cm, whereas the mean transformation error to the DE421 Mean Earth frame equals 5 cm.

Volcanism on Venus has never been directly observed, but several measurements indicate present-day activity. Volcanism could potentially play a role in climatic processes on Venus, especially in the sulfur cycle like on Earth. Observation of volcanic activity is the primary objective of future Venus spacecraft. However, there are many unknowns regarding its Venusian characteristics, like the condition at the vent, the volatile content and composition. Past modelling efforts have only studied explosive volcanic plume propagation over a limited range of flow parameters at the vent and in an idealised Venus atmospheric configuration. We propose to use the 1D FPLUME volcanic plume model in a realistic Venusian environment. In similar Venusian conditions, the height of the plume is consistent with past modelling. The present study shows that explosive volcanism would preferably reach 15 km of altitude. Under certain conditions, plumes are able to reach the VenSpec-H tropospheric altitude range of observations and even the 45 km cloud floor. For the first time, the impact of wind was quantified, and the super-rotating winds have a substantial impact by plume-bending of reducing the height of plumes. Contrary to the Earth, the atmospheric heat capacity depends greatly on temperature, and will disadvantage lower plumes and allow larger plumes to propagate at higher altitudes. The high latitude atmospheric environment, due to the thermal profile and weaker winds, is favorable to plumes reaching higher altitudes.

We report our near-infrared (NIR) follow-up observations of the gravitational wave (GW) event S240422ed using the Subaru Telescope/MOIRCS. S240422ed was initially classified as a black hole-neutron star merger with $>$ 99% probability of electromagnetic wave emission. We started follow-up observations 7.8 hours after the event. Over two nights, we observed 206 nearby galaxies in $Y$ and $K_{\rm s}$ bands down to about 21.4 and 21.1 AB mag (3$\sigma$), respectively. The total completeness of our survey based on galaxy $B$-band luminosity is 22%. As a result of our observations, five candidate counterparts were identified. We show that properties of these five objects are not consistent with kilonova such as AT2017gfo. Four objects are consistent with known classes of transients such as supernovae or dwarf nova outbursts. On the other hand, the nature of the remaining one object, which shows a red color and rapid decline, remains unclear. Although later analyses of GW signal reclassified S240422ed as likely terrestrial noise, our NIR observations provide valuable lessons for future NIR surveys for GW sources. We demonstrate that deep NIR follow-up observations as presented in this work would effectively constrain the presence of red kilonova even at 200 Mpc distance. We also discuss the importance of deep and wide NIR reference images and of understanding the properties and frequency of Galactic transients.

Context. Cygnus OB2, in the Cygnus X complex -- one of the most active star-forming regions of the Galaxy -- hosts hundreds of O- and B-type stars at different evolutionary stages. This association provides a unique laboratory to study massive star evolution and dynamics. However, despite extensive studies, the absence of a fast-rotating group ($v\sin{i}>200\,\mathrm{km\,s^{-1}}$) among the O-type population of Cygnus OB2 challenges current models of massive star evolution. Aims. Stellar rotation strongly impacts spectral line shapes of O-type stars, and high rotation can potentially lead to misclassifications. We investigate whether some stars in Cygnus OB2, classified at low spectral resolution as B0, are actually rapidly rotating late-O types. Such cases could explain the observed lack of fast rotators in Cygnus OB2. Methods. Accounting for rotation, we reclassify the known B0 population in Cygnus OB2, using the MGB tool and both new and pre-existing optical spectroscopy. Finally, we compute projected rotational velocities using iacob-broad. Results. About $19\,\%$ of the initial B0 population in Cygnus OB2 are, in fact, late-O types. Only 6 stars in the entire dataset show $v\sin{i}>200\,\mathrm{km\,s^{-1}}$, with just 1 new O-type star exceeding this threshold. Conclusions. In our study of Cygnus OB2, we continue to find a notable lack of fast rotators among its O-type population. We propose a combination of three factors as the most likely explanation: (i) the young age of Cygnus OB2 may imply that fast rotators has not been produced yet due to binary interactions; (ii) fast rotators may have been dynamically ejected from the core as runaway stars; (iii) local star formation conditions may hinder binary formation (reducing spin-up interactions) or result in slower rotational velocities at birth.

Context. Planetary nebulae are shells ejected by low- and intermediate-mass stars. The slow wind ejected by the asymptotic giant branch star is compressed by a fast stellar wind to produce an expanding gaseous shell surrounding a hot bubble. The shell is a source of thermal radio emission which shows a spectral index between -0.1 and 2. Only two planetary nebulae are known to show non-thermal radio emission indicating magnetic fields and non-thermal electrons. Aims. The aim of this paper is verification of presence of magnetic fields of planetary nebulae. Magnetic fields can have a significant influence on shaping planetary nebulae. Methods. We observed a sample of northern planetary nebulae in radio continuum at 144 MHz with the Low Frequency Array. We combined our observations with archival observations at higher frequencies. Results. The spectral indices in 30 planetary nebulae were below -0.1, indicating non-thermal radio emission. The majority of this sample consists of bipolar planetary nebulae, which are known to originate from binary central stars. Most of the nebulae have sizes larger than 20 arcsec. Magnetic fields and nonthermal emission may be common in smaller planetary nebulae, but can be suppressed by thermal emission. Our results suggest that different mechanisms can be responsible for the origin of magnetic fields and non-thermal emission in planetary nebulae.

Aims. Studying protostellar objects in their earliest stages, particularly during the Class 0 phase, provides key insight into the beginnings of planet formation and dust evolution. Disentangling the various components, however, is particularly challenging. High spatial and spectral resolution observations of molecular line emission with the Atacama Large Millimeter/submillimeter Array (ALMA) are therefore crucial for probing their complex environments. Methods. In this work, we present high-resolution ($\sim$30 au) ALMA observations at 1.3 millimeters of the Class 0 protostellar system IRAS4A2. Results. We detected large, well-traced outflows in HCN (3-2), H$_{2}$CO$(2_{1~2}-1_{1~1})$, and HCO$^{+}$ (3-2), along with numerous complex organic molecules (COMs) tracing central, more compact regions. Using moment maps, we analyzed the kinematics and spatial distributions of the molecular emission, revealing a wide range of spatial scales, from compact structures within the IRAS4A2 core at $\sim$8 au in radius, to extended $\sim$5000 au outflow emission. Specifically, we find that CH$_{3}$CDO and CH$_{3}$OCHO could be both good tracers of the disk, possibly tracing its rotation. Lines of OCS (22-21), SO$_{2}~(13_{3~11}-13_{2~12})$, HCN, H$_{2}$CO, and HCO$^{+}$, show more extended structures around IRAS4A2, likely tracing the envelope, disk, accretion shocks, the base of an outflow, and the outflow itself. Conclusions. Most COMs appear to trace distinct inner regions near the central protostar, while other molecules trace more extended structures, such as the envelope or outflows. The kinematics, emission patterns, and position-velocity diagrams suggest that individual molecules trace multiple components simultaneously, making it challenging to disentangle their true origins. Altogether, these findings highlight the complex spatial distribution within the IRAS4A2 system. Abridged.

Context: Diffuse $\gamma$-ray emission from cosmic ray (CR) protons scattering off the gas in the intracluster and intergalactic medium (ICM and IGM) remains out of reach for current observations. Detecting this emission would provide constraints on intergalactic magnetic fields as well as dark matter interaction models. Aims: We aim to provide estimates for diffuse $\gamma$-ray emission in the Fermi-LAT band from galaxy clusters and the cosmic web in the local Universe. Methods: In this work we show results from the first cosmological MHD simulation with an on-the-fly spectral CR model. We model CR injection at shocks, account for adiabatic energy changes and advection of CR protons, and obtain their $\gamma$-ray emissivity directly from the simulated CR energy density and spectra. For this we use constrained initial conditions, which evolve in a field closely resembling that of the local Universe, allowing direct comparison to Fermi-LAT data on massive clusters. Results: We find CR proton acceleration at all structure formation and accretion shocks in galaxy clusters and cosmic web filaments. These protons provide the basis for diffuse $\gamma$-ray emission in these regimes. The absolute value of the diffuse $\gamma$-ray emission in our simulation lies a few orders of magnitude below the current upper limits found by Fermi-LAT. Under the assumption of our model, a sensitivity of $F_\gamma < 10^{-11} \: \gamma~ \text{s}^{-1}~\text{cm}^{-2}$ is required for a detection of diffuse emission in Coma.

Studying the SDSS-DR16 quasar catalog, we detect a baryon acoustic oscillation (BAO) signal in the two-point angular correlation function with a statistical significance of $3\sigma$, at an effective redshift of $z_{\rm eff}=1.725$. Using a simple parameterization-comprising a polynomial plus a Gaussian function-we measure the transverse BAO scale as $\theta_{\rm BAO}=1.928^{\circ}\pm0.094^{\circ}$. This measurement is obtained from a narrow redshift shell, $z \in [1.72, 1.73]$ (i.e., $\Delta z=0.01$), thin enough that projection-effect corrections are negligible, making it only weakly dependent on the assumed fiducial cosmology. The only assumption adopted is isotropy in the computation of the correlation function, further ensuring that the result depends only weakly on specific cosmological-model hypotheses. We also investigate possible systematics that could affect the detection or significance of the BAO signal and find them to be subdominant or implausible. When combined with other transverse BAO measurements from the literature, our result shows good concordance-within the $1\sigma$ confidence level-with the cosmological parameter values reported by the Planck and DESI collaborations. This new measurement of the transverse BAO scale, obtained from the SDSS quasar sample with minimal cosmological-model assumptions, provides an additional independent constraint for updated statistical studies aimed at probing the nature of dark energy.

The question of whether genuinely young high-$\alpha$ stars exist has been discussed for over a decade since their discovery from asteroseismology of giant stars as it is challenging to break the degeneracy between the binary interaction and the genuinely young scenarios. Young high-$\alpha$ stars are hard to explain with traditional chemical evolution model as the high-$\alpha$ disk is typically associated with the early epoch of star formation in the Milky Way. Combined with recent advances of gyrochronology, and that $^7$Li can serve as an unambiguous indicator for identifying merger products in dwarfs thanks to its low burning temperature, we identified young high-$\alpha$ dwarf candidate stars through their fast rotation in a previous study. In this paper, we performed high-resolution spectroscopic follow-up of these candidates using Potsdam Echelle Polarimetric and Spectroscopic Instrument (PEPSI), and confirm 3 additional stars that are most likely genuinely young. Together with the star from the earlier paper, we find three out of four of them center around [Fe/H]=-0.5 dex, are ~5 Gyr old, and have a similar amount of elevated Li (~0.5 dex) and Al (~0.1 dex) compared to stars with matching $\log g$, $T_{\rm eff}$, Mg, and Fe within observational uncertainties, hinting at their common formation pathway.

Two main formation routes have been proposed for interstellar complex organic molecules (iCOMs): on dust grain surfaces and in the gas phase. Observing such molecules in protostellar outflow shock regions - provided that their ages are well-constrained - can help distinguish between these pathways by probing chemical evolution over time. This study focuses on the potential daughter-mother relationship of glycolaldehyde (CH$_2$OHCHO) and ethanol (C$_2$H$_5$OH), previously proposed in the literature. We test whether gas-phase reactions converting ethanol into glycolaldehyde derived in these works can explain the observed abundance of the latter in star-forming regions. We target the southern outflow of L1157, which hosts three shock regions, B0, B1 and B2, of increasing ages: about 900, 1500 and 2300 yr. We obtained high-resolution IRAM NOEMA maps of three lines of glycolaldehyde and one line of ethanol. We derived their abundances in the three shocks and used a pseudo time-dependent astrochemical model to simulate gas-phase and grain-surface formation scenarios for glycolaldehyde. Ethanol is assumed to form on grains and be released in the gas by shocks, where it is gradually converted into glycolaldehyde via the ethanol-tree reaction network. We present the first spatially resolved maps of glycolaldehyde and ethanol in the L1157 southern outflow, and more generally toward solar-like star forming regions. The abundance ratio [CH$_2$OHCHO]/[C$_2$H$_5$OH] increases from B1 to B2, consistent with model predictions. However, the model cannot reproduce all three shocked regions simultaneously, suggesting that one of the assumptions of our model, such as the same excitation temperature and grain composition in B0, B1 and B2, or gas temperature evolution, is wrong. Nonetheless, our modeling rules out the possibility that all the observed gaseous glycolaldehyde is a grain-surface product.

Millisecond pulsars (MSPs) are the oldest but fastest pulsars known to date. To explain how these pulsars could be formed, a new hypothesis was formulated: the recycling of pulsars, i.e the fact that a pulsar could accrete matter from a companion and been spun up. In this paper, we developed a population synthesis algorithm for pulsars which belong to a binary, in order to check whether most of the observed recycled pulsars were formed via an accretion mechanism and derive statistics about their properties, that are difficult to obtain through observations. We also make predictions for future surveys. Toward the presented objectives, we use the code Stellar EVolution for N-body (SEVN) to take into account all the binary processes and our own code to evolve each pulsar self-consistently by taking into account the secular evolution of a force-free magnetosphere, the magnetic field decay, gravitational braking and spatial evolution. Each pulsar is born in binary with a main sequence companion, and evolve to present time. The radio and $\gamma$-ray emission locations were modeled by the polar cap geometry and striped wind model, respectively. Our simulations seem to reproduce well the population of radio and $\gamma$-ray pulsars observed in the selected surveys. We also found that there should be less than $330$ unidentified pulsars in the Fourth Fermi-LAT catalogue of $\gamma$-ray sources (4FGL). High values of the viewing angle $\zeta$ seem to be needed to be able to observe the recycled pulsars, and it also seems difficult to observe recycled pulsars with an aligned rotation axis and magnetic axis (i.e., $\chi \leq 10$°). We find that only a small fraction, approximately $\sim 7.5\times10^{-3}$ %, of oxygen-neon white dwarfs (ONeWDs) in binary systems appear to contribute to the population of mildly recycled pulsars through accretion-induced collapse.

We report the serendipitous discovery of two RR Lyrae stars exhibiting significant s-process element enrichment, a rare class previously represented solely by TY Gruis. Our goal is to characterise these objects chemically and dynamically, exploring their origins and evolutionary histories. Using high-resolution spectroscopy from HERMES@AAT and UVES@VLT, we derived detailed chemical abundances of key s-process elements, carbon along with $\alpha$-elements. We also employed Gaia DR3 astrometric data to analyse their kinematics, orbital properties, and classify their Galactic population membership. We compared observational results with theoretical asymptotic giant branch nucleosynthesis models to interpret their enrichment patterns. Both stars exhibit clear signatures of s-process enrichment, with significant overabundances in second-peak elements such as Ba and La compared to first-peak Y and Zr. Comparison with AGB nucleosynthesis models suggests their progenitors experienced pollution of s-process-rich material, consistent with early binary interactions. However, notable discrepancies in dilution factors highlight the need for more refined low-metallicity asymptotic giant branch (AGB) models. We also explore and discuss alternative scenarios, including sub-luminous post-AGB-like evolution or double episodes of mass this http URL findings confirm the existence of s-process-enhanced RR Lyrae stars and demonstrate the importance of combining chemical and dynamical diagnostics to unveil their complex evolutionary pathways. Future detailed binary evolution modelling and long-term orbital monitoring are essential to resolve their formation scenarios and assess the role of binarity in the evolution of pulsating variables.

The alignments of galaxies across the large-scale structure of the Universe are known to be a source of contamination for gravitational lensing, but they can also probe cosmology and the physics of galaxy evolution in many ways. In this review, I cover developments in our understanding of intrinsic alignments over the past 25 years on: (1) different approaches to model intrinsic alignments across a range of scales, (2) existing observational constraints, (3) predictions from cosmological numerical $N$-body and hydrodynamical simulations, (4) mitigation strategies to account for their contamination to lensing observables and (5) cosmological and astrophysical applications. While the review focuses mostly on two-point statistics of intrinsic alignments, I also give a summary of other statistics beyond two-point. Finally, I point out some of the open problems hindering the understanding or application of intrinsic alignments and how they might be overcome in the future.

Dark-matter-dominated dwarf galaxies provide an excellent laboratory for testing dark matter models at small scale and, in particular, the ultralight dark matter (ULDM) class of models. Within the framework of self-interacting bosonic dark matter, we use the observed velocity-dispersion pro- files of seven dwarf spheroidal galaxies to constrain the parameters of ULDM. In our modeling, we account for the impact of the baryonic component on the velocity dispersion and ULDM halo structure. We find that in the regime of repulsively interacting ULDM, the self-interaction, which fits the observations, is almost negligible, consistent with non-interacting ULDM with a boson mass of approximately 1.6*10^(-22) eV. In contrast, for attractively interacting ULDM, the best fit corre- sponds to a smaller boson mass of about 1.3*10^(-22) eV, with self-interaction playing a significant role in shaping the dark-matter halo and thereby influencing the interpretation of observations.

With the growing number of gamma-ray monitors in operation, several research teams have adopted a strategy of joint operation and scientific duty to improve efficiency. A successful example is the GECAM-HXMT-SVOM (GHS) constellation collaboration, which sets a precedent for other gamma-ray monitor constellations. However, joint duty also presents challenges to Burst Advocates (BAs), including the increased number of triggers and, more importantly, the frequent switching between various systems due to incompatibilities among different missions, which complicates the situation. To address the current requirements of multi-wavelength and multi-messenger astronomy, we developed a customized framework for unified trigger processing within the GHS joint duty, named "BA's Rapid Evaluation and Analysis Kit for Formulating Alerts and Summary Tools" (BREAKFAST). This framework incorporates a series of automated, semi-automated, and manual pipelines designed to rapidly process triggers of prompt emissions in the gamma-ray band from different instruments, while maintaining flexible compatibility for future missions. The pursuit of BREAKFAST goes beyond merely providing trigger processing for BAs. BREAKFAST also aims to filtering high-value targets and guiding follow-up telescopes through rapid analysis and reporting, thus serving as an important bridge between prompt emission observations and afterglow observations. To this end, a suite of comprehensive analysis modules is included in BREAKFAST, particularly the specially designed module that predicts X-ray afterglow brightness based on prompt emission properties. The framework's effectiveness has already been demonstrated in recent observational campaigns, and it is expected to play a significant role in the discovery and observation of peculiar transients in the future.

The anisotropies of the Cosmic Neutrino Background (C$\nu$B) offer an ideal tool to test non-standard neutrino interactions, since they directly trace the perturbations in the neutrino distribution function. Here, we study how invisible neutrino decays impact the C$\nu$B anisotropies, in a framework where neutrinos decay non-relativistically to dark radiation and lighter neutrinos in a manner consistent with the measured mass splittings. For this purpose, we perform the first implementation of such a late-time neutrino decay scenario within a linear Einstein-Boltzmann solver, and compute the C$\nu$B angular power spectra from the Boltzmann hierarchy solutions for a range of lifetimes and decay channels. We find that neutrino decays leave very strong signatures on the C$\nu$B angular spectra, about two orders of magnitude larger than on the CMB angular spectra, particularly for lifetimes comparable to the age of the Universe. We show that a future polarized tritium target run of the PTOLEMY experiment, with sufficient counting statistics to measure just the first $\sim 15$ multipoles of the neutrino sky map, could test neutrino decay models that remain undetectable with CMB data.

The lensed Sunburst Arc ($z = 2.369$) hosts a young ($\sim2$--$4\,\rm Myr$), massive ($M_\star \sim 10^7\,M_\odot$), compact ($R_{\rm eff} \sim 8\,\rm pc$) Lyman-continuum (LyC) leaking super star cluster, which powers a compact ($< 10\,\rm pc$), high-pressure nebula at sub-solar metallicity $\sim0.2\,Z_\odot$ and with an anomalously elevated nitrogen-to-oxygen ratio $\log({\rm N/O}) \sim -0.2$. We present semi-analytic models and 3D magnetohydrodynamic simulations with radiative feedback in an attempt to reproduce this system. The results indicate that the progenitor giant molecular cloud (GMC) may have $M_{\rm cloud} \gtrsim 3 \times 10^7\,M_\odot$ and $R_{\rm cloud} \sim 70\,\rm pc$, corresponding to a surface density $\sim10^3$--$10^4\,M_\odot\,{\rm pc}^{-2}$. Incorporating feedback from individual Very Massive Stars (VMSs; $\ge 100\,M_\odot$) sampled from the Kroupa initial mass function, we find that their winds rapidly enrich $\sim 10^4\,M_\odot$ of nearby gas with nitrogen ($\sim 1\,$dex) and helium ($\sim 0.1$--$0.2\,$dex). In the first $1$--$3\,$Myr, some cold gas falls to the system center where a central cluster builds up from sub-cluster mergers. There, the gas is photoionized, pressurized, and chemically enriched by the newly formed VMSs, before being radiatively expelled in the next $\sim1\,\rm Myr$. We find that both VMS feedback and a high-surface-density progenitor GMC are necessary to reproduce the observed nebular properties, such as high N/O, high pressure, and stellar proximity. Low metallicity ($Z \le 0.004$) may be essential to avoid overproduction of carbon from WC stars. Such enrichment processes localized to compact starburst events may have caused strong nitrogen emission from dense ionized gas as observed in high-redshift galaxies such as GN-z11 and GS_3073.

HD 100453 disk is a prototypical companion-disk interaction system hosting a pair of spirals and a substellar companion. We present new noncoronagraphic high-contrast imaging observations of HD 100453 with $V$ filter on SPHERE/ZIMPOL. We combined high-contrast imaging data of the reflected light from 0.55 to 2.2 $\mu m$ using the $V$, $I'$, $J$, and $Ks$ band data of ZIMPOL and IRDIS at VLT/SPHERE. For each observational epoch, we corrected for the smearing effect to derive the intrinsic disk-integrated polarized flux. We derived a steady increase with wavelengths from $\hat{Q}_{\varphi}/I_{\star}(V)=0.3\%$ to $\hat{Q}_{\varphi}/I_{\star}(K)=1.2\%$. We applied reference differential imaging to extract the disk intensity for the $V$ and $Ks$ bands using star hopping observations. We obtained the first $V$-band total intensity for HD 100453 with ZIMPOL star hopping. The integrated total flux $I_{\rm disk}/I_{\star}(V)=1.5\%$ increases to $I_{\rm disk}/I_{\star}(K)=4.8\%$. Both the total intensity and the polarization fraction show red colors, and the intrinsic maximum degree of polarization increases moderately from $40\%$ to $55\%$. We then used RADMC-3D radiative transfer modeling with a parametrized Henyey-Greenstein phase function to constrain the dust properties. From the $V$ to $Ks$ band, the dominating dust in the outer disk has an increasing scattering albedo and degree of polarization, while the asymmetry parameter slightly decreases. The outer disk of HD 100453 contains sub-micron-sized low porosity grains/aggregates. The cavity in scattered light is not empty and is replenished with optically thin dust with a maximum size of $\leq 0.1\mu m$. The linear polarization is higher in the spiral region than in the other regions, suggesting different dust properties in those regions.

Recent BAO observations from DESI DR2 either hint at a possible dynamical dark energy component, which would worsen the Hubble tension, or at a 95\% credible interval for the summed neutrino mass hardly compatible with neutrino oscillation experiments. In this context, it is interesting to investigate constraints on neutrino masses, dark energy and the Hubble parameter that are agnostic to some aspects of the cosmological model. Here we choose to be agnostic to the value of the sound horizon at recombination, while sticking to standard assumptions regarding the time of recombination and the growth of structures. To be consistent, we also disregard information on the full shape of the CMB temperature and polarization spectrum on sub-degree scale. With such agnostic and conservative assumptions, we find that: (i) the dark energy evolution is well constrained by uncalibrated data on angular and luminosity distances, with a mild preference for dynamical dark energy, independently of the value of the sound horizon; (ii) large values of the Hubble rate are favored, $H_0=76.0^{+4.5}_{-4.3}$ km/s/Mpc (68\%CL), together with low values of the sound horizon, $r_{\rm s}=129.2^{+6.7}_{-8.1}$ Mpc (68\%CL); the SH0ES value of $H_0$ is thus marginally preferred over the low value returned by the standard inverse distance ladder analysis; (iii) the cosmological neutrino mass bound relaxes to $\sum m_\nu <1.6$ eV (95\%CL) and becomes well compatible with the normal and inverted neutrino mass schemes.

The robust association of electromagnetic candidates discovered during follow-up of gravitational-wave alerts is challenging, not only due to the large sky areas and broad distance uncertainties, but also due to the tens to hundreds of unrelated optical transients that are observed per event. We present a Bayesian ranking method to identify electromagnetic counterparts to GW events using only location information. The framework combines three-dimensional gravitational wave skymaps with host-galaxy information, a morphology-aware host association, empirical offset priors, and peculiar velocity corrections. We apply the method to GW170817 where it ranks AT2017gfo as the top candidate and correctly selects NGC\,4993 as the host. The approach is directly applicable to transient candidates with only location information and enables more efficient follow-up with prioritized candidates and leads to more reliable counterpart identification in current and future observing runs.

We report the precise X-ray localization of GRB~230906A, a short duration ($T_{90}\sim$0.9 s) burst with no optical or radio counterpart. Deep imaging with the Hubble Space Telescope detects a faint galaxy (G$^\ast$; $F160W\simeq26$ AB mag) coincident with the sub-arcsecond X-ray position. Compared with standard GRB galaxies, its faintness, compact size and color would suggest a high redshift ($z\gtrsim$3) host. However, our observations also reveal the presence of a galaxy group at $z\!\sim$0.453, confirmed spectroscopically with VLT/MUSE, with clear signs of interactions and mergers among group members. The GRB and its putative host project onto an extended ($\approx$180 kpc) tidal tail emerging from the group's central galaxy. The probability of a chance alignment is small ($P_{cc}\!\lesssim\!4$\%), we thus argue that the GRB and its galaxy G$^*$ reside within the group. Their peculiar location along the tidal debris suggests that an enhanced burst of star formation, induced by the galaxy merger, might have formed the progenitor compact binary $\lesssim$700 Myr ago. The compact binary later evolved in a neutron star merger which produced GRB 230906A and injected $r$-process material into the surrounding circumgalactic medium.

New, neutrinophilic mediators are one potential extension beyond the Standard Model of particle physics. Often, studies of neutrinophilic mediator consist of searching for direct evidence of its production and/or its tree-level virtual effect for generating strong neutrino self-interaction. In this work, we focus instead on the fact that such new mediators \textit{also} lead to deviations in neutrino-matter scattering via radiative corrections. With a mediator mass well below the electroweak scale, these effects are potentially observable in a variety of contexts, including coherent elastic neutrino-nucleus scattering (CEvNS), neutrino deeply-inelastic scattering ($\nu$DIS), and neutrino-electron scattering (e.g., at Borexino). Additionally, such effects lead to new contributions to the $Z$-boson decay width and to non-standard neutrino interactions relevant for long-baseline oscillation experiments. We explore all of these scenarios in some depth, building on the rich phenomenology associated with neutrinophilic mediators.

High-frequency gravitational waves represent an unexplored frontier of gravitational physics. While pulsar timing arrays reach nHz frequencies and ground-based interferometers probe the audio band, the regime above the MHz scale remains essentially untested. We propose a hybrid detection framework in which gravitons convert to photons via the Gertsenshtein effect in strong magnetic fields; the resulting electromagnetic signal is resonantly enhanced in a cavity and read out by atomic quantum sensors. Our approach combines multiple channels, including long-coherence Raman interferometry in alkali atoms, microwave Rydberg transitions, photoionization-based detection, and inner-shell electronic resonances. We derive the graviton-photon conversion probability, compute the induced photon flux, and estimate the shot-noise-limited spectral strain noise density. Depending on configuration, the projected shot-noise-limited sensitivity ranges from $\sqrt{S_{h, {\rm min}}} \sim 10^{-22}{\rm\,Hz^{-1/2}}$ in conservative microwave implementations to $\sqrt{S_{h, {\rm min}}} \sim 10^{-37}{\rm\,Hz^{-1/2}}$ in aggressive optical Raman schemes, potentially surpassing the cosmological bound from Big Bang Nucleosynthesis for stochastic backgrounds, while coherent bursts remain unconstrained and detectable. Such reach allows to test scenarios involving primordial black holes, topological defects, violent phase transitions, and (p)reheating, motivating advances in high-$Q$ cavities, strong-field magnets, and quantum-limited atomic sensors, with broad impact across quantum metrology and fundamental physics.

Our ability to infer the true source properties of colliding black holes from gravitational wave observations requires not only accurate waveform models but also their correct use. A key property when evaluating time-domain models is when to start the waveform: choosing a time that is too late can omit low-frequency power from higher order multipoles. By focusing on binary systems with total mass $\ge 200 \, M_{\odot}$, we show that current detectors are sensitive to this missing power and biased source properties can be obtained. We show that for systems with total mass $\lesssim 300 \, M_{\odot}$, mass ratio $\gtrsim 0.33$, and signal-to-noise ratio $\rho \gtrsim 20$, templates starting at $20 \, \mathrm{Hz}$ recover biased source properties. As the total mass increases, and the component masses become more asymmetric, templates starting from $13 \, \mathrm{Hz}$ recover biased properties. If the gravitational-wave signal is observed at signal-to-noise ratio $\rho < 20$, time-domain models can start from $20\, \mathrm{Hz}$ as statistical uncertainties dominate.

The effect of dissipation from the thermal bath of minimal warm inflation (MWI) on the temperature anisotropies of the cosmic microwave background (CMB) is investigated. A shift in the phase and peaks of the theoretical TT mode angular power spectrum of CMB for MWI with increasing dissipation strength hints at the dynamical nature of dark energy from warm inflation. Outcome of the study highlights the necessity of modifying the cosmological parameters, signalling a physics beyond the standard $\Lambda$CDM model of cosmology. The current findings may open a new direction for explaining the increased expansion rate in the early universe, thus alleviating some of the persistent tensions in cosmology, especially the Hubble tension.

We analyze radiative corrections to the Starobinsky model of inflation arising from self-interactions of the inflaton, and from its Yukawa couplings, $y$, to matter fermions, and dimensionful trilinear couplings, $\kappa$, to scalar fields, which could be responsible for reheating the Universe after inflation. The inflaton self-interactions are found to be of higher order in the Hubble expansion rate during inflation, and hence unimportant for CMB observations. In contrast, matter couplings to the Starobinsky inflaton can have significant effects on the spectral index of scalar CMB perturbations, $n_s$, and on the tensor-to-scalar ratio, $r$. Using a renormalization-group improved analysis of the effective inflationary potential, we find that the Planck measurement of $n_s$ constrains the inflaton coupling to light fermions in the Einstein frame: $y < 4.5 \times 10^{-4}$, corresponding to an upper limit on the reheating temperature $T_{\rm RH} < 2 \times 10^{11}~{\rm GeV}$, whereas the ACT DR6 measurement of $n_s$ corresponds to $3.8 \times 10^{-4} < y < 5.6 \times 10^{-4}$ and $1.7 \times 10^{11} ~{\rm GeV} < T_{\rm RH} < 2.8 \times 10^{11}~{\rm GeV}$, while the upper limits on $r$ provide weaker constraints. Planck data also imply a constraint on a trilinear inflaton coupling to light scalars in the Einstein frame: $\kappa \leq 4 \times 10^{12}~{\rm GeV}$, corresponding to $T_{\rm RH} \leq 4.2 \times 10^{13}~{\rm GeV}$. We further present constraints on inflaton couplings to massive fermions and scalars, and analyze constraints on couplings in the Jordan frame.

This study investigates the effects of strong magnetic fields on antikaon condensation in neutron star matter using the extended FSUGold model model. It is found that the presence of strong magnetic fields alters the threshold density of antikaon condensation significantly, which means the threshold density of antikaon condensation is shifted to higher density compared with the magnetic field-free case. In the presence of strong magnetic fields, the equation of state (EoS) becomes stiffer than that of the zero field case. The effects of the $\sigma$-cut scheme on the EoS are also researched when the appearance of antikaon condensation is occurred. Through careful choice of the parameter of the $\sigma$-cut scheme, we are able to produce a maximum mass neutron star heavier than 2$M_{sun}$.

In this work, we have constructed anisotropic bosonic dark-matter star (DMS) solutions in the context of a regularized four-dimensional Einstein$-$Gauss$-$Bonnet (4D EGB) gravity theory. Using dimensional regularization, we solve modified Tolman$-$Oppenheimer$-$Volkoff equations for a self-interacting complex scalar field in the dilute polytropic regime, $p_r = K \rho^2$, with anisotropy parameterized as $\sigma = \beta\, p_r \left( 1 - e^{-2\lambda} \right)$. We perform a comprehensive numerical analysis across the \((\alpha,\beta)\) parameter domain, where \(\alpha \in [0,8]~\mathrm{km}^2\) and \(\beta \in [-2,0]\), to examine mass$-$radius relations and evaluate multiple stability indicators including static equilibrium \(dM/dp_c\), sound-speed causality, the radial adiabatic index \(\Gamma_r\), and energy conditions. Positive Gauss$-$Bonnet coupling enhances both the maximum mass and compactness (e.g., \(M_{\rm max} \approx 1.62\, M_\odot\) at \(\alpha=0\) rising to \(\approx 2.09\, M_\odot\) at \(\alpha = 8~\mathrm{km}^2\)), while negative anisotropy reduces them (e.g., from \(\approx 2.21\, M_\odot\) at \(\beta=0\) to \(\approx 1.73\, M_\odot\) at \(\beta = -2\)). The resulting configurations remain statically stable up to the mass peak and satisfy physical criteria. This work extends previous isotropic boson-star analyses by systematically incorporating anisotropy within a regularized 4D EGB framework. These findings provide observationally relevant predictions for compact dark-matter objects under modified gravity.

Tidal forces acting on orbiting bodies arise from inhomogeneities in the gravitational field, generating stresses that can deform or even disrupt these objects. In this work, we analyze relativistic tidal forces associated with ultracompact objects described by static and spherically symmetric spacetimes, focusing on observers in circular geodesic motion. We show that, in contrast to the case of radial geodesics, tidal forces diverge as the orbit approaches null circular geodesics. As illustrative examples, we study two uniform-density stellar models: one isotropic and another supported purely by tangential stresses. We conjecture that the divergence of tidal forces near light rings may play a role in the nonlinear stability of ultracompact, horizonless objects.