Papers for Monday, Feb 09 2026

In the framework of Milgromian dynamics (MOND), galaxy clusters are known to exhibit a residual missing mass problem, with the baryonic mass falling short of the dynamical mass by about a factor of two. The baryon content of clusters is dominated by the intracluster medium (ICM), while the stellar contribution depends sensitively on the assumed stellar initial mass function (IMF). We re-evaluate the stellar and remnant masses in galaxy clusters by adopting the integrated galaxy-wide initial mass function (IGIMF) theory, which accounts for the dependence of the IMF on galaxy properties and star formation histories. Massive elliptical galaxies, characterized by high metallicities and short formation timescales, are inferred to form with top-heavy IMFs, leading to a substantial population of stellar remnants. Using observational data from WINGS and 2MASS for 46 nearby (z < 0.1) galaxy clusters, we compute stellar, remnant, and intracluster light masses and combine them with previously derived ICM masses. The resulting total baryonic masses are compared to MOND dynamical masses inferred from hydrostatic equilibrium. We find that the baryonic mass in stars, remnants and the ICM accounts for at least $88^{+5+2}_{-4-1}\%$ of the MOND dynamical mass. This constrains the kick velocities of the remnants and substantially alleviates the missing mass problem for galaxy clusters in MOND.

Given the persistence of various tensions in the "Cosmic Concordance" -- such as the "Hubble Tension", and the possible departure from LambdaCDM time evolution -- seen when combining complementary data sets (CMB studies, Baryon Acoustic Oscillations, Type Ia Supernovae, etc.), it remains an ongoing possibility for these to have a real cosmological origin. If one assumes such deviations to be real, a model-independent formalism (cosmography) is useful for locating the source of the problem with concordance cosmology. The extraordinarily good fit of LambdaCDM to the CMB data shows that it was a successful model of the universe at high redshift. Yet at lower redshift -- when the dark energy density becomes significant, and its precise physical nature becomes important -- the universe may have gone off the track of simple LambdaCDM. Here we use linear cosmography fits to binned Supernova data to reconstruct the detailed temporal history of the Hubble parameter, thus probing for interesting time-dependent behaviors of the expansion rate during and after the onset of cosmic acceleration. Using combined Type Ia supernovae from the Dark Energy Survey 5-Year data release and the Union2.1 compilation, we find intriguing hints of an oscillatory pattern in the Hubble parameter during the acceleration era. While these hints are low-significance, and not robust under different redshift binnings, we present this work as a proof-of-concept demonstration of this method for reconstructing the Hubble parameter evolution, which may be useful for the voluminous Supernova data sets anticipated to become available during the next few years.

Calibrating out per-antenna signal chain effects is an essential step in analyzing radio interferometric data. For drift-scanning arrays, robustly calibrating the data is especially challenging due to the lack of the ability to track a calibration source. Consequently, calibration strategies for drift-scanning arrays are limited by our knowledge of the radio sky at large, as well as the direction-dependent instrument response. In the context of 21 cm cosmology, where small calibration errors can conspire to overwhelm the cosmological signal, it is therefore crucially important to develop calibration strategies that are capable of accurately calibrating the data in the presence of sky or instrument modeling errors. In this paper we present CorrCal, a covariance-based calibration strategy for redundant radio interferometric arrays. CorrCal is a hybrid calibration strategy that leverages the strengths of traditional sky-based calibration and redundant calibration in a computationally efficient framework that is fairly insensitive to modeling errors. We find that the calibration errors from CorrCal are unbiased and far below typical thermal noise thresholds across a wide range of modeling error scenarios. We show that CorrCal is computationally efficient: our implementation is capable of evaluating the likelihood and its gradient in less than a second for 1,000-element class arrays using just a single laptop core. Given CorrCal's computational efficiency and robustness to modeling errors, we anticipate that it will serve as a useful tool in the analysis of radio interferometric data from current and next-generation experiments targeting the cosmological 21 cm signal.

We present a unified Bayesian assessment of model comparison and data-set consistency for LCDM (cold dark matter plus a cosmological constant) and minimal extensions (neutrino mass, spatial curvature, constant or evolving dark energy) using cosmic microwave background (CMB), baryon acoustic oscillation (BAO), and type Ia supernova (SN) data. The major results are summarized in the first three figures. We quantify model preference with Bayesian evidence and assess consistency with complementary evidence- and likelihood-based diagnostics applied uniformly across data-set combinations. For the models considered, updated Planck processing systematically improves internal CMB consistency (low-$\ell$ versus high-$\ell$, and primary CMB versus CMB lensing). The preference for a closed geometry and an associated ``curvature tension'' with BAO and/or CMB lensing are largely confined to earlier Planck likelihood implementations and weaken substantially when using updated CMB processing and more recent BAO measurements. Apparent evidence for evolving dark energy in CMB+BAO+SN combinations depends sensitively on the specific pairing of CMB and SN likelihoods: plausible alternatives shift inferred tensions by more than $1\,\sigma$ and can completely reverse the preferred model. Allowing a free neutrino mass tends to absorb residual shifts without introducing new inconsistencies, and we do not find robust evidence for a standalone $\tau$-driven discrepancy once the full likelihood context is accounted for. We conclude that claims of a required update of our standard cosmological model from LCDM to $w_0w_a$CDM are premature.

Reverberation mapping (RM) determines the mass of black holes (BH) in active galactic nuclei (AGNs) by resolving the BH gravitational sphere of influence in the time domain. Recent RM campaigns yielded direct BH masses through dynamical modeling for a sample of 32 objects, spanning a wide range of AGN luminosities and BH masses. In addition, accurate BH masses have been determined by spatially resolving the broad-line region with GRAVITY for a handful of AGNs. Here, we present a detailed analysis of Hubble Space Telescope images using surface-brightness profile fitting with state-of-the-art programs. We derive AGN luminosity and host-galaxy properties, such as radii and luminosities for spheroid, disk, and bar (if present). The spheroid effective radii were used to measure stellar velocity dispersion from integral-field spectroscopy. Since the BH masses of our sample do not depend on any assumption of the virial factor needed in single-epoch spectroscopic mass estimates, we can show that the resulting scaling relations between the mass of the supermassive BHs and their host galaxies match those of quiescent galaxies, naturally extending to lower masses in these (predominantly) spiral galaxies. We find that the inner AGN orientation, as traced by the broad-line region inclination angle, is uncorrelated with the host-galaxy disk. Our sample has the most direct and accurate MBH measurements of any AGN sample and provides a fundamental local benchmark for studies of the evolution of massive black holes and their host galaxies across cosmic time.

Observations indicate that the disc population in galaxy clusters has undergone rapid evolution, transitioning from a dominance of blue spirals to red S0s over the past $\sim7$ Gyr. We build a simplified cluster evolutionary model in the $\Lambda$CDM framework to constrain the characteristic timescales of this transformation. In our model, field spirals joining the cluster are subject to ram-pressure stripping (RPS), which removes their gas reservoir leading to the quenching of their star formation on a timescale $t_{\rm s}$, and to an (initially) unspecified mechanism that transforms them into S0s on a timescale $t_{\rm m}$. We assume that $t_{\rm s}$ and $t_{\rm m}$ are independent and both power-law functions of $M_\star/M_{\rm cl}$, the galaxy-to-cluster mass ratio. We constrain our model using the observed distribution of spirals and S0s in a color-mass plane from the OmegaWINGS and EDisCS cluster surveys at $z\simeq0.055$ and $z\simeq0.7$. Our best-fit model reproduces the data remarkably well and predicts evolutionary trends for the main morphological fractions in agreement with previous studies. We find typical $t_{\rm s}$ between $0.1$ and $1$ Gyr, compatible with previous estimates. A surprisingly strong anti-correlation between $t_{\rm s}$ and $M_\star/M_{\rm cl}$ is required in order to suppress the formation of red, low-mass spirals at low redshift, which we interpret as driven by orbit anisotropy. Conversely, $t_{\rm m}$ depends very weakly on $M_\star/M_{\rm cl}$ and has typical values of a few Gyr. The inferred morphological evolution is compatible with that resulting from the ageing of the stellar populations in galaxies abruptly quenched by ram pressure stripping: we confirm spectrophotometric ageing as a key channel for the spiral-to-S0 transition in galaxy clusters, with secular evolution playing a secondary role.

High-energy neutrinos provide a potentially powerful and distinctive probe for dark matter (DM) - neutrino interactions, particularly in environments with enhanced DM densities, such as the DM spikes predicted to form around supermassive black holes (SMBHs) at the center of active galactic nuclei (AGN). Recent observations by the IceCube Neutrino Observatory, which identified four AGN, namely TXS 0506+056, NGC 1068, PKS 1424+240, and NGC 4151 as neutrino sources, provide a unique opportunity to search for signatures of these interactions. In this study, we use IceCube data to derive the most stringent constraints to date on both the energy-dependent and energy-independent DM-neutrino scattering cross-sections. We perform a statistical analysis using data from individual sources as well as a combined (stacked) analysis of all four sources. Our strongest limits arise from the stacking analysis, yielding an upper bound of $\sigma_{0} \le 8\times 10^{-39}$ cm$^2$ for an energy-independent cross-section and $\sigma_{0} \le 10^{-39}$ cm$^2$ for a linearly energy-dependent cross-section, both at 90\% confidence level, particularly in scenarios involving the adiabatic growth of black holes.

Stellar high-energy radiation is a key driver of atmospheric erosion and evolution in exoplanets, directly affecting their long-term habitability. We present a comprehensive study on stellar high-energy radiation and its impact on exoplanetary atmospheres, leveraging data from the \textit{SRG/eROSITA} all-sky survey. Our sample consists of 3750 main-sequence stars identified by cross-matching with \textit{Gaia} DR3. Utilizing X-ray spectral fits from the \textit{eROSITA} catalog, we computed X-ray ($L_X$) and combined extreme-ultraviolet (EUV) luminosities ($L_{\mathrm{EUV}}$), which we used to derive XUV fluxes at the habitable zone ($F_{\mathrm{XUV,HZ}}$). We find that the majority of stars in our sample are significantly more XUV-active than the Sun, with habitable zone fluxes ranging from $10^0$ to $10^5$ erg~cm$^{-2}$~s$^{-1}$. The ratio of $L_{\mathrm{XUV}}/L_{\mathrm{bol}}$ is found to be higher for cooler, magnetically active stars, highlighting their potentially hazardous nature for planetary atmospheres. Applying the energy-limited escape model, we computed atmospheric mass-loss rates for hypothetical earth-like planets located at the habitable zone of each star. We also present local maps for distances up to $500$~pc of the average XUV flux, revealing ``hazard zones'' where stellar radiation could significantly influence planetary atmospheric evolution. This work demonstrates the power of X-ray surveys in constraining the high-energy environments of exoplanets and underscores the critical role of stellar activity in planetary habitability.

Our understanding of the formation and evolution of planetary systems has made major advances in the past decade. This progress has been driven in large part by the Atacama Large Millimeter/submillimeter Array (ALMA), which has given us an unprecedented view of Solar System bodies themselves, and of the structure and chemistry of forming exoplanetary systems. Within our own Solar System, ALMA has enabled the detection of new molecules and isotopologues across moons and comets, as well as placing new constraints on the compositions and histories of small bodies through thermal emission observations. In this article, we highlight some key areas where ALMA has contributed to a deeper understanding of our Solar System's formation and evolution, and place these discoveries in the context of our evolving understanding of protoplanetary disks.

Some cosmological hydrodynamical simulations predict that outflows driven by active galactic nuclei (AGN) play a key role in powering the Ly$\alpha$ nebulae observed around high-redshift quasars. In such simulations, AGN feedback seeded as powerful outflows leads to extended and luminous nebulae whose morphology and surface-brightness profiles accurately reproduce the observations, while suppressing AGN feedback leads to compact and faint nebulae. This link might arise from outflows opening up a channel for Ly$\alpha$ photons to escape from the galactic nucleus to the circumgalactic medium (CGM). The main aim of this paper is to test this theoretical prediction using observations, by comparing the physical properties of outflows and Ly$\alpha$ nebulae. We analyze integral-field unit data obtained with VLT/ERIS and GEMINI/GNIRS to trace the ionized gas in the interstellar medium (ISM) of a sample of six quasars at $z\sim2-3$, using the [O III] emission line. We detect powerful outflows in all the quasars of our sample, with velocities $>1500~\mathrm{km~s^{-1}}$ and kinetic energies $ \gtrsim 2\times10^{43}~\mathrm{erg~s^{-1}}$. Four of our quasars are spatially resolved and show signs of extended [O III] emission out to distances $>2$ kpc from the central supermassive black hole. When excluding one outlier, we find a positive monotonic correlation between the outflow power and the Ly$\alpha$ nebulae size ($\rho=0.89$, $p=0.03$) and luminosity ($\rho=0.6$, $p=0.28$). Additionally, we find evidence of spatial alignment between the ionization cone and the inner and brightest regions of the Ly$\alpha$ nebula. Our results provide tentative evidence in support of the theoretical prediction that AGN-driven outflows at ISM scales open a low-optical-depth path for central Ly$\alpha$ photons to reach the CGM and create extended nebulae.

The abundance of elements in the interstellar medium (ISM) is a key facet for many fields of astrophysical study. In the soft X-ray spectra, absorption by interstellar gas can result in deep absorption features that affect continuum measurements. In this paper, we focus on measuring the abundance of interstellar iron and neon from the column densities observed in soft spectra from XMM-Newton and Chandra for various low mass X-ray binaries (LMXBs), which allows for a direct probe of elemental abundances. As a noble gas, neon will not deplete into solid form, thus providing a benchmark with abundances determined via UV spectroscopy. We find that, when assuming Fe is 90\% depleted into grains, [Fe/Ne]$ = -0.523\pm0.025$, [Fe/H]$ + 12 = 7.482\pm0.016$, and [Ne/H]$ + 12 = 8.012\pm0.022$, which are the tightest observational constraints on these abundances to date, while being consistent with literature which uses protosolar abundances. We also test how depletion into solid grains and scattering affect the results. The choice of depletion fraction can affect the abundance measurement by roughly $5\%$, and that the inclusion of a scattering component can affect abundance measurements by $\sim1-7\%$.

(Abridged) A thorough interpretation of the current and upcoming generation of cosmological observations requires unprecedented large-scale, high-resolution simulations spanning multiple cosmological models and parameters. The realization of these computationally demanding simulations poses a crucial technical challenge. We present beyond - $\Lambda$CDM implementations of the Dynamic Zoom Simulations (DZS) method, a performance-enhancing technique tailored for large-scale simulations that produce lightcone-like outputs. This approach dynamically decreases the resolution of a simulation in the regions that are not in causal connection with the observer, saving computational resources without directly affecting the physical properties within the lightcone. We implemented the DZS algorithm in two state-of-the-art codes supporting non-standard cosmologies, namely modified $f(R)$ gravity in Arepo and dark sector interactions in Gadget4. We analyzed result accuracy and performance gains across resolution, simulation volume and model by comparing runs performed with and without the DZS algorithm. Our DZS reproduce the lightcone halo mass function, sky-projected massmaps, and matter and weak lensing convergence power spectra with an accuracy of $\simeq$ 0.1% or higher in most cases. In terms of performance, DZS runs in our test simulations can save up to $\sim$ 50% runtime compared to the non-DZS counterparts. A scaling to larger simulated volumes suggests that performance gains could improve by an additional $\sim$ 20% at the resolution levels of current state-of-the-art simulations. The validation of the DZS algorithm in non-standard models demonstrates that this technique can enable cost effective, large-scale ($\gtrsim$ 1 cGpc/h) simulations with state-of-the-art resolution, providing the computational framework needed to constrain and help the interpretation of forthcoming data.

Multi-messenger astronomical observations of neutron stars, together with more precise calculations and constraints coming from dense matter microphysics, are generating tension with regard to equations of state models used to describe neutron star cores. Assuming an abrupt first-order phase transition with a slow conversion speed between phases, we propose different slow stable hybrid star configurations aiming to reconcile all current constraints simultaneously; within this framework, we also introduce a novel non-CSS parametrization to the quark matter equation of state and discuss its strengths and limitations. We analyze our model results in conjunction with a review of other relevant theoretical possibilities existing in the literature. We found that modern neutron star observations seem to favor the existence of some type of exotic matter in the neutron star cores; in particular, our slow stable hybrid star scenario remains a proposal capable of satisfying these constraints. However, due both to the existing skepticism regarding some of the adopted hypotheses in most extreme neutron star measurements and to the precise adjustment needed for the equation-of-state parameters, significant tension and open questions remain.

Stars in the Galactic disk have mean radial velocities $\overline{v}_R$ that oscillate as a function of angular momentum $J_\varphi$. This `$J_\varphi$-${\overline{v}}_R$ wave' signal also exhibits a systematic phase shift when stars are binned by their dynamical temperatures. However, the origin of the wave is unknown. Here we use linear perturbation theory to derive a simple analytic formula for the $J_\varphi$-$\overline{v}_R$ signal that depends on the equilibrium properties of the Galaxy and the history of recent perturbations to it. The formula naturally explains the phase shift, but also predicts that different classes of perturbation should drive $J_\varphi$-$\overline{v}_R$ signals with very different morphologies. Ignoring the self-gravity of disk fluctuations, it suggests that neither a distant tidal kick (e.g., from the Sgr dwarf) nor a rigidly-rotating Galactic bar can produce a qualitatively correct $J_\varphi$-$\overline{v}_R$ wave signal. However, short-lived spiral arms can, and by performing an MCMC fit we identify a spiral perturbation that drives a $J_\varphi$-${\overline{v}}_R$ signal in reasonable agreement with the data. We verify the analytic formula with test particle simulations, finding it to be highly accurate when applied to dynamically cold stellar populations. More work is needed to deal with hotter orbits, and to incorporate the fluctuations' self-gravity and the role of interstellar gas.

X-ray binaries, powered by black holes, neutron stars, or white dwarfs accreting matter from a companion star, are among the brightest beacons in galaxies, outshining the Sun by a factor of millions. Most emit primarily above 0.3 keV in X-rays, but cooler thermal sources peaking in the extreme ultraviolet (EUV) would be much more difficult to detect due to astronomy's critical blind spot in EUV. Here, we report the discovery of a remarkable new class of luminous, point-like, non-nuclear X-ray objects in galaxies-hypersoft X-ray sources -- that have been missed by all previous surveys to date. Detected primarily or exclusively below 0.3 keV, with 0.15-0.3 keV to 0.3-1.0 keV photon ratios >8, the most luminous examples radiate >1E38 erg/s in the narrow X-ray band, with spectral models indicating even greater bolometric luminosities, largely emitted in the EUV. They rank among the most energetic sources in galaxies, yet their EUV-peaking spectra evaded earlier detections. We propose that hypersoft sources are X-ray binaries spanning multiple physical classes, including accreting white dwarfs or post-nova systems-potential Type Ia supernova progenitors-and systems hosting accreting black holes. Beyond their elusive nature, they may play crucial role in ionizing gas within galaxies.

We revisit the luminosity function (LF) of compact HII regions in the context of the inertial--inflow model, in which massive stars assemble over extended, mass-dependent timescales. The comparison of the compact-HII-region LF with that of OB stars is traditionally used to estimate the compact-HII-phase lifetime and is often cited as evidence for the classical ``lifetime problem" of HII regions. We show that once stellar growth during the ionizing phase is included, the LF comparison instead constrains massive-star formation timescales, so the lifetime problem turns into evidence for prolonged growth. We illustrate the principle with a simple analytic model and then forward-model the two LFs with Monte Carlo realizations. We also derive revised Galactic LFs for compact HII regions and OB stars from the Red MSX Source survey and the Alma Luminous Star catalogue. The joint LF constraints imply a growth law where the formation time is about 2 Myr for a $60\,M_\odot$ star, with a square-root dependence on mass. The revised OB-star LF exhibits a statistically significant knee at $\log_{10}(L_{\rm k}/L_\odot)=5.0$, while the HII-region LF knee occurs at lower luminosity, as expected in the interpretation that HII regions are powered by stars that are still growing in mass. We conclude that massive stars in the Milky Way form over Myr timescales that increase with their final mass.

The modern generation of wide-field galaxy surveys, such as LSST, Euclid and Roman, will enable studies of dwarf galaxies $(10^6 \leq M_\ast / M_\odot \leq 10^9)$ beyond the Local Group (LG) in unprecedented detail. Improved theoretical understanding of this population is necessary to guide these observations, since predictions in this regime are generally limited to specific environments like the LG. We present predictions for the population of dwarf galaxies from the TNG50 run of the IllustrisTNG suite of cosmological hydrodynamical simulations, focusing on the environments within $1 < D / \mathrm{Mpc} < 25$ of LG analogues at $z = 0$. In the simulated sample, there are $\sim 1,000$ and $\sim 12,000$ dwarf galaxies within $10$ and $25$ Mpc, respectively. We compare our results with the 50 Mpc Galaxy Catalog and estimate that current observations are highly incomplete at low masses: for $10^6 \leq M_\ast / M_\odot \leq 10^7$ $(-13 \lesssim M_r \lesssim -10)$, we find completeness fractions of $\sim 23 \%$ within $10$ Mpc and $\sim 4 \%$ within $25$ Mpc. The simulated galaxies below the completeness limits of the observations exist in a range of environments, with notable populations of field dwarfs at all distances and satellites around centrals with masses $10^8 \lesssim M_\ast / M_\odot \lesssim 10^{11}$ within $10-25$ Mpc. We find that there are $\sim 8$ times more quiescent dwarf galaxies in the TNG50 sample than are currently cataloged. Our results suggest that upcoming observations should uncover a substantial population of dwarf galaxies, and that $\gtrsim 15 \%$ of these will be red, currently quenched galaxies in the field.

Mass transfer between non-degenerate stars is a fundamental but still poorly understood process in binary evolution. The commonly used rotationally limited accretion prescription in detailed binary evolution simulations that account for stellar rotation generally yields low accretion efficiencies that are difficult to reconcile with several observational constraints. We present a physically-motivated mass-accretion prescription in which accretion or decretion disks regulate the angular momentum transported to the accretor, thereby allowing for continued accretion at near-critical rotation. The accretion efficiency can be calculated from the conservation of the mass and the angular momentum of the disk. Analytical estimates show that the accretion efficiency depends on stellar rotation and mass ratio for direct impact accretion, and additionally on stellar radius and orbital separation in the disk accretion regime. The overall mass-weighted accretion efficiencies are close to the values expected near the threshold rotation rate, where the accreted specific angular momentum declines sharply. Applying this model to binary evolution simulations, we find that rotationally limited accretion systematically underestimates Be-star masses in Be+subdwarf O/B-type star (sdOB) systems, whereas the disk-star coupling model can produce more massive Be stars that are consistent with observations. The final binary component masses depend not only on accretion efficiency but also core-envelope mass ratio, which itself depends sensitively on the assumed overshooting. We find that our new disk-star coupling model with reduced overshooting yields component masses for Be+sdOB systems that are in closer agreement with observations.

We explore the notion that cosmological models that modify the late-time expansion history cannot simultaneously fit the SH0ES collaboration's measurements of the Hubble constant, DESI baryon acoustic oscillations data, and Type Ia supernova distances. Adopting a few simple phenomenological models, we quantitatively demonstrate that a satisfactory fit with a model with late-time expansion history can only be achieved if one of the following is true: 1) there is a sharp step in the absolute magnitude of Type Ia supernovae at very low redshift, $z\sim 0.01$, or 2) the distance duality relation, $d_L(z)=(1+z)^2d_A(z)$, is broken. Both solutions are trivial in that they effectively decouple the calibrated SNIa measurements from other data, and this qualitatively agrees with previous work built on studying specific dark-energy models. We also identify a less effective class of late-time solutions with a transition at $z\simeq 0.15$ that lead to a more modest improvement in fit to the data than models with a very low-z transition. Our conclusions are largely unchanged when we include surface brightness fluctuation distance measurements, with their current systematic uncertainties, to our analysis. We finally illustrate our findings by studying a physical model which, when equipped with the ability to smoothly change the absolute magnitude of Type Ia supernovae, partially resolves the Hubble tension.

I describe an `oct-tree' N-body code which randomly shifts, reorients, and resizes the root cell at each time step. Averaging over a plurality of root cell positions and orientations statistically restores translational and rotational invariance. The potentials and forces which result can be much more accurate than those obtained from a single force calculation. In this paper, the principle of averaging is tested on static configurations. The next paper will show how this technique can substantially improve global energy, momentum, and angular momentum conservation at a negligible computational cost.

Using real metre-sized asteroid Earth impactors from the last decade, we ask the question: ``If the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) had been operating over the last 10 years, how many imminent impactors would it have observed and discovered pre-impact?'' We feed 216 large impacts detected by global coverage orbital sensors through the LSST Solar System Survey Simulator Sorcha, and find that it would have made 99 observations of 28 unique objects, and discovered one object four days pre-impact. Recently proposed pipeline modifications would increase this discovery number to seven. Scaling our results to take into account the biases of our reference population, we estimate that LSST will discover 12 +/- 3 imminent impactors over its nominal 10 year survey, with an average warning time of 3.5 days. While this is at the low end of previous estimates of 1 - 10 discoveries per year, the significant increase in warning time compared to the current average (9 hours across 11 impactors) will bring significant opportunities for follow-up telescopic observations, deployment of specialised equipment for fireball observations, and planetary defence operations. We also show that the LSST will provide substantial precovery data for impactors discovered by other surveys, instantly lengthening observation arcs and thereby reducing the orbital and impact location uncertainties. In some cases, these observations may also enable the linkage of telescopic observations with observed fireballs post-impact, providing valuable pre-impact astrometric and photometric data. This has significant implications for both asteroid research and planetary defence.

The rotations of cometary nuclei are known to change in response to outgassing torques. The nucleus of comet 41P/Tuttle-Giacobini-Kresak exhibited particularly dramatic rotational changes when near perihelion in 2017 April. Here, we use archival Hubble Space Telescope observations from 2017 December to study the post-perihelion lightcurve of the nucleus and to assess the nucleus size. From both Hubble photometry and non-gravitational acceleration measurements we find a diminutive nucleus with effective radius r = 500+/-100 m. Systematic optical variations are consistent with a two-peaked (i.e., rotationally symmetric) lightcurve with period 0.60+/-0.01 days, substantially different from periods measured earlier in 2017. The spin of the nucleus likely reversed between perihelion in 2017 April and December as a result of the strong outgassing torque. We infer a dimensionless moment arm k = 0.013, about twice the median value in short-period comets. The lightcurve range of 0.4 magnitudes indicates a projected nucleus axis ratio greater than 1.4:1, while the active fraction of the nucleus decreased from 2.4 in 2001 (suggesting augmentation of the gas production by sublimating coma ice grains) to 0.14 in 2017, a result of long-term modification of the surface. We find that the physical lifetime of this small nucleus to spin-up is short compared to the reported 1500 year dynamical time spent in the current orbit. Two limiting reconciliations of this inequality are suggested. The nucleus could be in a state of unusually strong activity, leading us to over-estimate the average mass loss rate and outgassing torque and so to under-estimate the physical lifetime. Alternatively, the nucleus could be the surviving remnant of a once larger body for which outgassing torques were less effective in changing the spin.

Femtolensing of gamma ray bursts (GRBs) are vastly studied to constrain the primordial black hole (PBHs) lighter than $10^{-13}$ solar mass and may close the window for PBH dark matter. In this case, wave optics formalism is required and carefully implemented in our analysis. Incorporating the GRB observational data from Swift XRT, we perform the statistic analysis of PBH lensing, comparing with null hypothesis where BAND model is used to parametrize the GRB spectrum. We found few GRB data manifest the spectral fringe which characterize the feature of femtolensing by PBHs, and the analysis shows moderate statistical preference in terms of goodness of fit. Conversely, since most of the fitting to GRB spectral data do not improved with PBH lensing, we utilize to obtain upper bound on the PBH fractional abundance with respect to dark matter. However, the robust constraint cannot be achieved, unless the size of GRBs are smaller than $5\times10^{7}$ m for PBH mass around $5\times10^{-15}$ solar mass.

There is increasing evidence that, in the very late phase of stellar evolution before core collapse, massive stars have winds with large mass loss rates that give rise to a dense circumstellar medium (CSM) surrounding the progenitor star. After core collapse, a shock wave forms when the supernova ejecta interacts with this CSM. In such an interaction, the nuclei in the CSM can undergo diffusive shock acceleration and reach very high energies. We consider such a model, which includes magnetic field amplification from the non-resonant streaming instability, enhancement to the abundance of heavy-ions, and composition-dependent acceleration. Applying this to several supernova subclasses, we find that IIn supernovae can supply a dominant fraction of the observed super-knee cosmic-ray (CR) flux from $\sim{\rm few}\times10^{15}\,{\rm eV}$ to $\sim{\rm few}\times10^{17}\,{\rm eV}$ and is consistent with recent LHAASO measurements above the CR knee. This systematic model also explains the increasingly heavy nuclear composition in this energy range.

The Massive and Distant Clusters of WISE Survey 2 (MaDCoWS2) is a WISE-selected catalog of galaxy clusters at $0.1<z<2$ covering an effective area of $>6000$ deg$^2$. In this paper, we derive splashback radii for this cluster ensemble from galaxy density profiles and constrain the mass threshold of the survey as a function of redshift. We use MaDCoWS2 cluster candidates at $0.4\leq z \leq 1.65$ divided into subsamples with different signal-to-noise (S/N$_{\rm P}$) and redshifts, cross-correlated with galaxies from the CatWISE2020 catalog, to obtain average surface density profiles. We perform a Markov Chain Monte Carlo analysis to derive parameter estimates for theoretical models consisting of orbiting and infalling terms. A distinct splashback feature is detected in all subsamples. The measured splashback radii span from $0.89^{+0.02}_{-0.02}h^{-1}$ comoving Mpc/cMpc ($0.61^{+0.02}_{-0.02}h^{-1}$ proper Mpc/pMpc) at $\overline{z}=0.45$ to $1.27^{+0.05}_{-0.05}h^{-1}$ cMpc ($0.53^{+0.04}_{-0.04}h^{-1}$ pMpc) at $\overline{z}=1.54$. We also find that splashback radii increase with $S/N_{\rm P}$ at fixed redshift. The resultant splashback radii constrain the redshift dependence of the mass of MaDCoWS2 clusters at fixed $S/N_{\rm P}$. We calculate $M_{\rm 200m}$ from the radii using a relation based on a cosmological simulation. MaDCoWS2 $M_{\rm 200m}$ values derived from the simulation-based relation are lower than the expected values based on weak-lensing observations. More robust mass constraints will come from calibrating splashback radii derived from galaxy density profiles with weak lensing shear profiles from facilities such as $\textit{Euclid}$, Rubin, and $\textit{Roman}$.

The splashback radius, the radius of the apocenter of the first orbit of infalling material, is a measurable quantity marking the boundary between a galaxy cluster and its infalling region. We report detections of splashback radii in total light stacks, i.e. image stacks centered on the cores of galaxy clusters. Our analysis uses Wide-field Infrared Survey Explorer (WISE) W1 and W2 images of 83,345 candidate clusters at $0.5 \lesssim z \lesssim 1.9$ from the Massive and Distant Clusters of WISE Survey 2 (MaDCoWS2). The clusters are organized in stacks by redshift and signal-to-noise ($S\slash N$) ratios. We adopt a statistical approach, using 1000 bootstrap realizations to determine the median projected splashback radius and its confidence interval in a given bin. We compare our splashback radii with the measurements made by K. Thongkham et al. on a similar sample of MaDCoWS2 clusters using galaxy-cluster cross-correlation and find that they are consistent, although our method yields larger error bars. Our main systematic error is the accuracy of the background subtraction, but its impact remains small: the consistency of K. Thongkham et al. and our results suggests that neither method suffers from large systematics. The sensitivity of total light stacking to the contribution of faint galaxies can be advantageous to locate splashback radii when only the brightest galaxies are detected in individual images, such as at high redshifts. We present a potential application of this new technique to probe the evolution of the stellar mass in cluster infalling regions.

We present a highly scalable, MPI-parallelized framework for reconstructing the initial cosmic density field, designed to meet the computational demands of next-generation cosmological simulations, particularly the upcoming ELUCID-DESI simulation based on DESI BGS data. Building upon the Hamiltonian Monte Carlo approach and the FastPM solver, our code employs domain decomposition to efficiently distribute memory between nodes. Although communication overhead increases the per-step runtime of the MPI version by roughly a factor of eight relative to the shared-memory implementation, our scaling tests-spanning different particle numbers, core counts, and node layouts-show nearly linear scaling with respect to both the number of particles and the number of CPU cores. Furthermore, to significantly reduce computational costs during the initial burn-in phase, we introduce a novel ``guess'' module that rapidly generates a high-quality initial density field. The results of the simulation test confirm substantial efficiency gains: for $256^3$ particles, 53 steps ($\sim$54 CPU hours) are saved; for $1024^3$, 106 steps ($\sim$7500 CPU hours). The relative gain grows with the number of particles, rendering large-volume reconstructions computationally practical for upcoming surveys, including our planned ELUCID-DESI reconstruction simulation with $8192^3$ particles, with a rough estimation of 720 steps ($\sim$37,000,000 CPU hours).

HP~Tau/G2 is a luminous, short-period, fast-rotating G-type weak-line T Tauri star with a large radius, an oblate shape with gravity-darkening, little circumstellar material, and centered in a slowly expanding cloud cavity. It is an X-ray source and a variable nonthermal radio source. It forms, together with the late-type T Tauri star KPNO 15, a pair of oppositely directed walkaway stars launched when a multiple system broke apart ~5600 yr ago. Momentum conservation indicates a mass of G2 of only ~0.7 Msun, much lower than the ~1.9 Msun determined from evolutionary models. G2 is virtually a twin of FK Com, the prototype of a class of evolved stars resulting from coalescence of W UMa binaries. We suggest that G2 became a very close and highly eccentric binary during viscous evolution in the protostellar stage and with KPNO 15 formed a triple system, which again was part of a larger unstable group including the binary G3 and the single G1. Dynamical evolution led to multiple bound ejections of KPNO 15 before it finally escaped after ~2 Myr. As a result the G2 binary recoiled and contracted 5600 yr ago, became Darwin unstable and merged in a major outburst ~2000 yr ago. The nearby compact triple system G1+G3 was also disturbed, and broke up 4900 yr ago, forming another walkaway pair. The G5 star HD 283572 has similar unusual properties, indicating that G2 is not a pathological case. G2 is now fading towards a new stable configuration. YSO mergers may be rather common and could explain some FUor eruptions.

This paper demonstrates that the stellar masses of galaxies in the Galaxy and Mass Assembly (GAMA) survey, originally derived via stellar population synthesis modelling, can be accurately predicted using only their absolute magnitudes and colour indices. A central contribution of this work is the demonstration that this long-standing inference problem can be solved using an exceptionally simple machine-learning model: a fully connected, feed-forward artificial neural network with a single hidden layer. The network is trained exclusively on synthetic galaxies generated by the SHARK semi-analytic model and is shown to transfer effectively to real observations. Across nearly 3.5 dex in stellar mass, the predicted values closely track the GAMA SED-derived masses, with a typical scatter of ~0.131 dex. These results demonstrate that complex deep-learning architectures are not a prerequisite for robust stellar mass estimation, and that simulation-trained, lightweight machine-learning models can capture the dominant physical information encoded in broad-band photometry. The method is further applied to 17,006 GAMA galaxies lacking SED-derived masses, with photometric uncertainties propagated through the network to provide corresponding error estimates on the inferred stellar masses. Overall, this work establishes a computationally efficient and conceptually transparent pathway for simulation-to-observation transfer learning in galaxy evolution studies.

The unpredictability of solar filament eruptions presents major challenges for forecasting space weather, as such eruptions frequently drive coronal mass ejections (CMEs) that impact the heliosphere. While nearby flux emergence is often linked to their destabilisation, the specific characteristics of both the emerging flux and the filament that determine whether an eruption occurs remain unclear. We report observations of a quiescent filament that did not erupt following the nearby emergence of active region NOAA 13270 and a subsequent C-class flare in April 2023. Our analysis combines multi-viewpoint extreme ultraviolet (EUV) imaging and X-ray imaging with EUV spectroscopy, radio imaging and measurements of, and extrapolations from, the photospheric magnetic field. We identify the formation of a coronal null point and fan-spine topology at the interface between the active region and filament which exhibited persistent slow reconnection, indicated by chromospheric brightenings, persistent radio emission, and plasma upflows. Our results indicate that ongoing reconnection and jets can relieve magnetic stress and enable filament stability, even when under strong perturbation. We suggest that the orientation of emerging flux relative to the ambient field is a critical parameter in filament evolution, and provide observational constraints for models of filament stability and eruption.

Adaptive optics can be used to mitigate the effects of atmospheric turbulence on imaging systems, but the correction is only partial, and deconvolution is often required to improve the resolution. This results in entire optical/digital systems, which are traditionally designed sequentially, i.e. , the adaptive optics system is optimised first, and the restoration algorithms are designed a second time. Studies on optical/digital systems have shown that jointly optimizing the whole system is a better alternative. We propose to extend these co-design strategies to the design of an adaptive optics-assisted imaging system. We derive a simple criterion that takes into account the source properties and the entire optical/ digital system performance. To illustrate its interest, we use it to optimize the wavelength distribution between the wavefront sensor and the imaging camera. In addition, we explore the potential of using multiple imaging channels operating at different wavelengths as a means of making an imaging system robust to turbulence strength and source magnitude variations. Later, any parameter of the optical/digital system, if not the entire system itself, could be optimized this way.

Context. Protoplanetary disks are the birthplace of planets. As such, they set the initial chemical abundances available for planetary atmosphere formation. Thus, studying elemental abundances, molecular compositions, and abundance ratios in protoplanetary disks is key to linking planetary atmospheres to their formation sites. Aims. We aim to derive the sulfur abundance and the C/O ratio in the AB Aur disk using interferometric observations of CS and SO. Methods. New NOEMA observations of CS 3-2 towards AB Aur are presented. We used velocity-integrated intensity maps to determine the inclination and position angles. Keplerian masks were constructed for all observed species to assess the presence of non-Keplerian motions. We use the CS/SO ratio to study the C/O ratio. We compare our present and previous interferometric observations of AB Aur with a NAUTILUS disk model to gain insight into the S elemental abundance and C/O ratio. Results. We derive an observational CS/SO ratio ranging from 1.8 to 2.6. Only NAUTILUS models with C/O > 1 can reproduce such ratios. The comparison with models points to strong sulfur depletion, with [S/H]=8e-8, but we note that no single model can simultaneously fit all observed species.

During its annual conference in 2025, the French Society of Astronomy \& Astrophysics (SF2A) hosted, for the fifth time, a special session dedicated to discussing the environmental transition within the French A\&A research community. During the 2025 workshop, the goal was to review four contemporary topics within the context of environmental transition actions and discussions: (1) institutional actions, (2) the early-career researchers singularity, (3) research infrastructures and tools, and (4) the geopolitical conditions under which A\&A research remains possible. The workshop concluded with a round-table discussion that brought together the various speakers so that every participant could express their ideas.

Neutron monitors are a standard tool for high-precision continuous monitoring of galactic cosmic ray flux variations arising from variations in heliospheric conditions and solar activity for space weather applications. These measurements form the basis for solving the inverse problem of determining the cosmic ray anisotropy vector beyond the magnetosphere. To support such studies, periodic latitude measurements are necessary to determine the coupling functions of primary and secondary cosmic rays variations. The aim of this work is to develop and characterize a modernized standard neutron monitor based on a CHM-15 boron thermal neutron counter and a data acquisition system designed for marine expeditionary studies of cosmic ray variations. Modern nuclear physics experimental methods and the principles of microprocessor-based data acquisition systems were used to solve this problem. The results of test trials and of continuous monitoring showed that the characteristics of the upgraded and standard neutron monitor are similar, and the ease of use, compactness, and stability allow us to conclude that the mobile neutron detector can be used in expeditionary conditions with limited access for maintenance personnel.

The magnetar Swift J1834.9-0846 presents a significant challenge to neutron star spin-down models. It exhibits two key anomalies: an insufficient rotational energy loss rate to power its observed X-ray luminosity, and a braking index of $ = 1.08\pm 0.04$, which starkly contradicts the canonical magnetic dipole value of $n=3$. To explain these anomalies, we develop a unified spin-evolution model that self-consistently integrates magnetic dipole radiation, gravitational wave emission, and wind braking. Within this framework, we constrain the wind braking parameter to $\kappa \in [13, 37]$ from the nebular properties, finding it contributes substantially (17%-51%) to the current spin-down torque. Bayesian inference reveals that the birth period is poorly constrained by present data and is prior-dependent, indicating a millisecond birth is allowed but not required. Furthermore, we constrain the number of precession cycles to $\xi \sim 10^{4}$--$10^{5}$, and our analysis favors a toroidally-dominated internal magnetic field configuration as the most self-consistent explanation for the low braking index. Finally, we assess the continuous gravitational-wave detectability. The present-day signal is undetectable. However, the early-time signal might have reached the projected sensitivity of next-generation gravitational-wave observatories, such as the Advanced Laser Interferometer Gravitational-Wave Observatory (aLIGO) and the Einstein Telescope (ET), although a confident detection would require exceptionally stable rotation, an assumption considered highly optimistic for a young magnetar. This work establishes a unified framework that links magnetar spin-down with their interior physics and multi-messenger observables, providing a physically consistent interpretation for Swift J1834.9-0846 and a new tool for understanding similar extreme neutron stars.

We construct exotic compact objects composed entirely of self-interacting asymmetric fermionic dark matter governed by a repulsive Yukawa potential with massive dark interaction boson. By considering the structural, tidal, and rotational properties of solar mass self-gravitating dark matter systems, and contrasting them against purely baryonic neutron stars, described by the well understood SLy4 equation of state, we hope to shed some light on the place of dark compact systems in the context of gravitational wave astronomy, specifically due to the difficulty parsing mass and radius data from events with no electromagnetic counterpart. Here we consider systems composed of 1 GeV and 10 GeV dark matter. Relevant compact objects are then analysed and simulated as both static bodies, and rotating systems governed by the Hartle-Thorne formalism to second order. Here within we highlight the differences in key tidal and rotational properties encoded in gravitational wave signals, and analyse how dark objects may mimic or distinguish themselves to current and future gravitational wave observatories.

Numerical end-to-end simulation in Adaptive Optics (AO) is a key tool in the development of complex systems, from the initial design to the commissioning phase. Based on our previous experience with PASSATA, we decided to develop a new AO simulation framework in Python language, naming it SPECULA (short for: Scalable Parallel Execution of Computations Upscaling Large Adaptive optics simulations). Following an object-oriented approach, the physical entities are modeled as processing objects connected to each other to exchange data objects. A simulation is run by providing its description instead of writing and executing a specific script. The Python language and its library flexibility allowed us to write one single code that can be run on CPU and GPU platforms. We put a strong focus on computational efficiency, relying on CuPy and its interface to access the CUDA-stream mechanism. Moreover, SPECULA is capable of distributed computations over multiple processing nodes, making it suitable to run in an HPC environment, as tested on the Italian supercomputer Leonardo. SPECULA can also be used in laboratory environment to implement a hybrid simulation, allowing us to interface simulated and concrete objects: this feature was demonstrated in the Adaptive Optics laboratories at Arcetri Observatory. In this paper, we describe the main characteristics of SPECULA, show some relevant examples of its use, and finally draw our goals for the future.

The stellar mass-size relation is a sensitive probe of how environment shapes galaxy structure. We analyse this relation in the local Universe for galaxies in compact groups (CGs), low-mass groups ($M_{\rm vir} \leq 10^{13}~M_{\odot}$), and high-mass groups, comparing them to field galaxies using data from the Southern Photometric Local Universe Survey. Galaxies are classified as early types (ETGs; $n \geq 2.5$, $(u-r)_0 \geq 2.3$), late types (LTGs; $n < 2.5$, $(u-r)_0 < 2.3$), transition galaxies (TGs; $n < 2.5$, $(u-r)_0 \geq 2.3$), and others (OGs; $n \geq 2.5$, $(u-r)_0 < 2.3$). We find that ETGs and OGs show no significant environmental dependence: their mass-size slopes and intercepts are statistically consistent across CGs, groups, and the field. LTGs also follow similar relations in the field and in most groups, with only a modest tendency for LTGs in CGs to be smaller at fixed stellar mass. By contrast, TGs display a clear environmental signal: in groups the slope steepens to $\alpha \sim 0.4$ (versus $\alpha \sim 0.2$ in the field) and their sizes are smaller than in the field, with non-overlapping 95\% posterior intervals. These trends suggest that TGs in denser environments are more structurally evolved, likely owing to enhanced bulge prominence and fading of the outer disc, consistent with the Sérsic-index distributions, which show an excess of TGs with $n_r \gtrsim 1.5$ in groups and CGs. Our findings highlight TGs as an environmentally sensitive population, providing insight into the structural transformation of galaxies in group environments.

The JWST provides exoplanet transit observations with unprecedented spectral coverage, enabling detailed atmospheric characterisation. However, systematics introduced during data reduction can lead to small but significant uncertainties that propagate into atmospheric retrievals, making it essential to assess their impact on inferred properties. We aim to quantify the impact of different JWST/NIRSpec PRISM data-reduction processes as well as the relevance of saturation on the retrieved atmospheric properties of WASP-39b. We also assess whether or not these biases are comparable to those introduced by assumptions made in atmospheric modelling, particularly in the treatment of aerosol extinction. We perform nested-sampling Bayesian retrievals using MultiNest and forward models generated with the Planetary Spectrum Generator. Six independently reduced spectra are analysed, considering both the full wavelength range and versions excluding the saturated region. We further test the effect of including three different cloud-opacity parameterisations. Differences among JWST/NIRSpec data-reduction pipelines lead to substantial variations in retrieved atmospheric properties of WASP-39b, often exceeding one order of magnitude, comparable to uncertainties from modelling assumptions. Excluding the saturated region reduces inter-pipeline dispersion but increases parameter degeneracies. This highlights the need for robust and homogeneous calibration strategies. The results also confirm that JWST data possess the sensitivity required to probe aerosol spectral behaviour, although such constraints remain strongly dependent on the adopted data-reduction strategy.

A large proportion of observed white dwarfs (WDs) show evidence of debris disks, remnants of the former planetary systems, and/or signatures of heavy elements in their atmospheres, induced by the accretion of planetary matter onto their surfaces. The observed abundances are the result of the balance between the accretion flux and the dilution of this planetary material by internal transport processes. A recent study showed that more massive DA WDs are less polluted than smaller mass ones. It was suggested that the reason could be related to the formation of planetary systems when these stars were on the main sequence. The aim of this work is to test how internal dilution processes, including thermohaline convection, change with WD masses, and whether such an effect could account for variations in the observed pollution. We computed the efficiency of atomic diffusion and thermohaline convection after the accretion of heavy elements onto WDs using static DA models with various masses, effective temperatures, and hydrogen contents. We confirm that thermohaline convection is always more efficient in diluting accreted elements than atomic diffusion, as previously shown in the literature. However, we find that element dilution by thermohaline convection is less efficient in massive WDs than in smaller mass ones, due to their larger internal density. We showed that the differences in observed heavy element pollution in WDs according to their masses cannot be explained by the dilution induced by atomic diffusion and thermohaline mixing alone. Indeed, the pollution by planetary system accretion should be more easily detectable in massive WDs than in low-mass ones. We discuss other processes that should be taken into account before drawing any conclusion about the occurrences of planetary systems according to the mass of the star on the main sequence.

We present sensipy, an open-source Python toolkit for simulating observations of transient astrophysical sources, particularly in the high-energy (HE, keV-GeV) and very-high-energy (VHE, GeV-TeV) gamma-ray ranges. The most explosive events in our universe are often short-lived, emitting the bulk of their energy in a relatively narrow time window. Due to often rapidly fading emission profiles, understanding how and when to observe these sources is crucial both to test theoretical predictions and efficiently optimize available telescope time. The information extracted from the tools included in sensipy can be used to help astronomers investigate the detectability of sources considering various theoretical assumptions about their emission processes and mechanisms. This information can further help to justify the feasibility of proposed observations, estimate detection rates (events/year) for various classes of sources, and provide scheduling insight in realtime during gamma-ray and multi-messenger observational campaigns.

We present results of a pilot observational wide-field VLBI campaign on probing scattering properties of the partly ionized interstellar medium towards the Galactic plane sky region between $28^\circ<l<36^\circ$ and $|b|<1^\circ$. This covers the region where the Galactic bar connects to the spiral arms and where a lot of star formation is currently ongoing. The Very Long Baseline Array (VLBA) observations of the whole region were performed in a special mode with multiple phase centers at L-band (1.4 -- 1.8 GHz) during April-June 2022 and a year later complemented by sessions at S (2.2 -- 2.4 GHz) and C-band (4.6 -- 5.0 GHz) partially covering the pilot region. We found compelling evidence that target sources are subject to scattering. The total detection rate in L, S and C-bands is 1.5, 3.4 and 9.2 per cent, respectively, and approximately scales with the square of the observation frequency. The low rate values imply that scattering is strong. Its power is non-uniform across the Galactic plane and it can be approximated by a Gaussian with a width of about $2^\circ$ peaking at the Galactic mid-plane. One of the brightest sources of the field shows anisotropic scattering, with a $\lambda^2$ dependence of its observed angular size, along a position angle of $26^\circ$ aligned with the line of constant Galactic latitude. We estimate the turbulence dissipation scale $r_\text{in}\approx1500$ km toward the source J1833+0015.

We present a comprehensive framework for predicting the detection prospects of supermassive black hole binaries (SMBHBs) by future gravitational wave (GW) observatories, examining both space-borne detectors (LISA, Taiji, TianQin) and next-generation pulsar timing array (PTA) combined with the Square Kilometre Array (SKA-PTA). Leveraging dual active galactic nucleus (AGN) fractions and AGN X-ray luminosity functions, we systematically evaluate the detectable SMBHB populations with a detection threshold of signal-to-noise ratio $\geq 5$ for each GW observatory. Our analysis reveals that space-borne detectors are expected to identify approximately $\sim 1 \text{--} 2$ to $\sim 20$ events per year, depending on the SMBHB orbital evolution prescriptions. On the other hand, SKA-PTA demonstrates the potential to reach the first GW detection from individual SMBHBs within a few years of observation and achieve detectable GW source counts of $10^2 \text{--} 10^3$ after about 10 years, depending on PTA configurations. These facilities will significantly improve SMBHB detectability and enable characterization of their properties across different frequency bands.

Strong outflows from active galactic nuclei are frequently observed in objects with lower coronal X-ray luminosity. This intrinsic X-ray weakness is considered a requirement for the formation of radiatively driven winds. To obtain an unbiased view on the connection between X-ray emission and the presence of powerful winds in the most luminous quasar phase, we present an X-ray analysis of a sample of extremely luminous, radio-quiet quasars with signatures of strong outflows in their rest-frame ultraviolet (UV) emission spectra. We study the $Chandra$ X-ray spectral properties of 10 objects, selected from the Sloan Digital Sky Survey Data Release 16 quasar catalogue based on their UV luminosities and ${\rm C}_{\rm IV}$ emission line blueshifts, comparing them to typical optically blue quasars. Our analysis reveals that seven out of 10 quasars in our sample have photon indices $\Gamma>1.7$. Only two out of 10 objects exhibiting outflows with velocities exceeding 1400 km/s are X-ray 'weak', consistent with the fraction of X-ray 'weak' objects generally observed in quasar populations. Notably, one of the objects identified as X-ray 'weak' is likely an intrinsically X-ray 'normal' quasar that is heavily obscured. We observe a tentative indication at a $\sim$2$\sigma$ confidence level that the correlation between the excessively low X-ray flux level and the presence of ${\rm C}_{\rm IV}$ emission-line outflows might emerge at wind velocities greater than 3000 km/s. Our study provides additional evidence that the relationship between X-ray emission and the presence of winds is intricate. Our findings emphasise the need for X-ray observations of a larger sample of UV-selected quasars with confirmed strong emission-line outflows to unravel the nuanced interplay between winds and X-ray emission.

We derive robust constraints on the initial abundance of evaporating primordial black holes (PBHs) using the reionization history of the Universe as a cosmological probe. We focus on PBHs that inject electromagnetic (EM) energy into the intergalactic medium (IGM) after recombination, in the mass range $3.2\times 10^{13}\,\mathrm{g} \lesssim M_{\rm PBH} \lesssim 5\times 10^{14}\,\mathrm{g}$. For each PBH mass, we compute the redshift-dependent energy injection from Hawking evaporation using \texttt{BlackHawk}, fully accounting for the time evolution of the PBH mass and the complete spectrum of emitted Standard Model particles and gravitons. The resulting photons and electrons are propagated through the primordial plasma using \texttt{DarkHistory}, which self-consistently models EM cascades and determines the fraction of injected energy deposited into ionization, excitation, and heating of the IGM. These modifications to the ionization and thermal histories are incorporated into a Gaussian Process reconstruction of the free-electron fraction based on low-$\ell$ CMB polarization data from the \textit{Planck} mission. This non-parametric approach allows for a statistically well-defined separation between exotic high-redshift energy injection and late-time astrophysical reionization, allowing PBH evaporation to be constrained through its contribution to the high-redshift optical depth. Requiring consistency with current CMB measurements, we obtain upper limits on the initial PBH abundance that are robust against reionization modeling uncertainties and systematically more conservative than existing bounds, reflecting the fully numerical and time-dependent treatment of Hawking evaporation and energy deposition. Our results demonstrate the power of reionization observables as a precision probe of PBH evaporation and other scenarios involving late-time energy injection.

Host galaxies of tidal disruption events (TDEs) show enhanced central stellar concentration and are preferentially found in post-starburst and green valley populations. This connection has led to the proposal that TDE host galaxies likely have gone through recent mergers. We conduct a new morphological study of 14 TDE host galaxies, using the r-band images from Sloan Digital Sky Survey (SDSS), Dark Energy Camera Legacy Survey (DECaLS), and Ultraviolet Near-Infrared Optical Northern Survey (UNIONS), with the images from the latter two surveys having much higher depth and resolution than SDSS. We examine galaxy structures using conventional methods and also apply diagnostics of merger activity from machine learning models. Consistent with previous studies, our results show that TDE host galaxies are ~16% more centrally concentrated when compared to non-TDE-host controls. However, surprisingly, TDE hosts lack any indication of recent merger activity from both morphological analysis and our machine learning merger classifier. Instead, our results reveal that TDE host galaxies are approximately 1.5 to 2.5 times more likely to have bar-like or ring-like structures than their controls. This enhancement is even more prominent for TDEs in the green valley, with the factor reaching almost 3. Based on these results, we propose that bar-driven secular evolution, instead of merger, likely dominates the recent evolution of TDE hosts found in the green valley, which can simultaneously explain their distinctive nuclear properties and enhanced TDE rates.

The availability of asteroseismic constraints for tens of thousands of red giant (RG) stars has opened the door to robust age estimates, enabling time-resolved studies of different populations of stars in the Milky Way. This study leverages data from Kepler, K2, and TESS, in conjunction with astrometric data from Gaia DR3 and spectroscopic constraints from APOGEE DR17 and GALAH DR3, to infer parameters for over 17,000 RGs. We use the code PARAM to homogeneously infer stellar properties considering in detail the sensitivity of our results to different choices of observational constraints. We focus on age estimation, identifying potentially unreliable age determinations, and highlight stars with unreliable $\Delta\nu$ measurements based on comparisons using Gaia luminosities. These are particularly relevant in K2 data due to the short duration of the observations of each campaign, and therefore important to characterise for Galactic archaeology studies where the spatial range of K2 is a benefit. Thanks to the combination of data from different missions we explore trends in age, mass, and orbital parameters such as $R_\mathrm{g}$ and $Z_\mathrm{max}$, and examine time-resolved [$\alpha$/M]-[Fe/H] planes across different Galactic regions. Additionally, we compare age distributions in low- and high-$\alpha$ populations and chemically selected ex situ stars. The study also extends known mass-[C/N] ratio relationships to lower masses. The catalogues resulting from this work will be instrumental in addressing key questions in Galactic archaeology and stellar evolution, and to improve training sets for machine-learning-based age estimations.

The LISA mission will be the first observatory to detect gravitational waves from space within the millihertz frequency band. Magnetic forces have an important impact on the instrument's sensitivity below the millihertz. Hence, monitoring the magnetic environment within each of the LISA spacecrafts is of utmost importance. In this Letter we present the characterization of the coils that were used in LISA Pathfinder (LPF) when operating as magnetic sensors in the audio frequency band. The necessity of implementing this type of magnetometer is presented in order to monitor high frequency magnetic signals from the electronics on-board. We show that the LPF coils have a performance one order of magnitude better than the current requirements set by the LISA mission at the low end of the audio-band frequency. The LPF coils are able to measure a magnetic noise level of 1.45 $\rm pT/\sqrt{\rm{Hz}}$ at 50 Hz and 0.17 $\rm pT/\sqrt{\rm{Hz}}$ at 500 Hz. Additionally, the LPF coils can reach a magnetic noise floor of 0.1 $\rm pT/\sqrt{\rm{Hz}}$ at frequencies above 1 kHz.

This paper is the first in a series investigating magnetohydrodynamic instabilities that may contribute to angular-momentum transport and magnetic-field evolution in stellar radiative zones. We focus on shear-driven instabilities, specifically the Goldreich-Schubert-Fricke (GSF) instability and the magnetorotational instability (MRI), which are expected to play key roles in the internal dynamics of radiative regions. We carried out a local linear stability analysis using a numerical approach that extends beyond classical limiting cases and includes stabilizing effects such as stratification and magnetic tension, allowing the exploration of realistic flow regimes. These results were validated through a global mode analysis in a Taylor-Couette configuration. We recovered the standard MRI and azimuthal MRI stability criteria and quantified the effects of stratification, magnetic tension, and diffusion on their growth. In strongly sheared regimes, we derived a new criterion for the magnetised GSF (MGSF) instability and clarified the transition from SMRI to MGSF as stratification and magnetic effects narrow the unstable domain. We also provided approximate growth-time formulae that identify the dominant instability under given stellar conditions and can be implemented in 1D stellar evolution codes. Global Taylor-Couette calculations validate the local WKB analysis. Applied to subgiants and young red giants, our results show that shear-driven instabilities can grow rapidly for magnetic fields below 100 kG. Strong axial fields (100 kG) confined to the hydrogen-burning shell suppress instabilities unless the shear is sufficiently distant. These results support incorporating our criteria and growth estimates into stellar evolution models to assess the efficiency of shear-driven transport.

The discovery of 2024 YR4 presented the planetary defense community with the most significant impact threat in almost two decades, reaching level 3 on the Torino scale. The community, now mature and well-organized, responded with a global observational effort. Astrometric measurements, forming the basis for orbital refinement and impact prediction, were a central component of this response. In this paper, we present the astrometric data collected by the international community, from the time of discovery until the object became too faint for all existing observational assets, including JWST. We also discuss the coordination role played by the International Asteroid Warning Network, and the importance of publicly available image archives to enable precovery searches.

Galaxy clusters are tracers of the large scale structures of the Universe, making the time evolution of their mass function dependent on key cosmological parameters, such as the cosmic matter density or the amplitude of density fluctuations $\sigma_8$. Accurate measurements of cluster's total masses are therefore essential, yet they can be challenging, particularly for clusters with complex morphologies, as simple mass profiles are often adopted to fit the measurements. In this work, we focus on the Frontier Fields galaxy clusters: a sample of six extremely massive systems, that, in most cases, exhibit highly complex mass distributions. The BUFFALO survey extended the Hubble Space Telescope observations for the Frontier Fields galaxy clusters, providing high-resolution multi-band imaging within a few Mpc. Combining this high-quality imaging dataset with ancillary spectroscopy, we produce weak-lensing catalogues with very high source densities, about 50 sources/arcmin$^2$. This allows us to robustly estimate the individual weak-lensing cluster masses and quantify the sensitivity of these measurements on different factors, such as the cluster centring, the uncertainty on the redshift distribution or the foreground contamination and boost factor correction. This provides a data-driven analysis of the different sources of systematics that can impact such measurements. We find that the largest sources of systematic bias arise for the most disturbed clusters, such as the multi-modal, merging galaxy cluster Abell 2744. This analysis sets a comprehensive framework for assessing the impact of systematics on the weak-lensing estimates of cluster masses, and in particular, in the case of unrelaxed clusters. This can play a key role in forthcoming cosmological analyses based on wide-field surveys such as Euclid and the Legacy Survey of Space and Time of the Rubin Observatory.

The quasar OP 313 was discovered in December 2023 in very-high-energy $\gamma$ rays above 100 GeV, enabling for the first time a complete characterization of its emission. However, the lack of updated measurements of its accretion disk, broad line region and dusty torus hampers a detailed interpretation of the role of accretion in the observed $\gamma$-ray production. We intend to characterize, during high-activity states, the external photon fields contributing to the IR-to-UV emission$-$namely dusty torus, broad line region and accretion disk$-$and investigate potential variability and blurring effects on the broad emission lines. We present new spectroscopic observations of OP 313 with the NOT and TNG telescopes to characterize its optical spectrum and variability with respect to archival observations from SDSS. We measure the luminosity of different broad emission lines, characterizing the broad line region, accretion disk and dusty torus. We measure the Mg II emission line, with an average flux of $F_{\mathrm{Mg \ II}} = (0.85 \pm 0.11)\times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$. Its equivalent width and luminosity are consistent with a constant line with a variable non-thermal continuum. From the stable Mg II line we derive a constant luminosity of the thermal components, $\log(L_{\mathrm{BLR}} \ \mathrm{[erg \ s^{-1}]}) = 44.91 \pm 0.19$, $\log(L_{\mathrm{disk}} \ \mathrm{[erg \ s^{-1}]}) = 45.91 \pm 0.19$ and $\log(L_{\mathrm{torus}} \ \mathrm{[erg \ s^{-1}]}) = 44.70 \pm 0.16$, and estimated a BH mass of $\log(M_{BH}/M_{\odot})=8.36 \pm 0.18$, in line with with that derived from the C III] line. These characteristics and the indicator of the accretion rate from the disk/Eddington luminosity ratio $\lambda =L_{AD}/L_{Edd} = 0.23 \pm 0.10$, along with a high Compton dominance, favour a FSRQ-like nature, contrary to the argued changing-look nature.

We present an automated pipeline for operational short-term forecasting of coronal mass ejection (CME) magnetic field structure at L1, coupling arrival time prediction, in situ detection, and iterative flux rope reconstruction, following near-real-time remote-sensing CME identification. The system is triggered by new entries in the CCMC DONKI database and first applies the drag-based ELEvo model to determine whether an Earth impact is expected and estimate arrival time. This estimate defines a temporal window constraining the search for CME signatures in real-time L1 in situ solar wind data, where the magnetic obstacle (MO) is automatically detected using the deep learning model ARCANE. Upon MO onset, iterative reconstructions with the semi-empirical flux rope model 3DCORE are performed, using a Monte Carlo fitting scheme, producing continuously updated forecasts of the remaining magnetic field profile. We evaluate the pipeline using 3870 archived DONKI entries and archived NOAA real-time in situ data from 2013-2025, assessing forecast performance at different stages of MO observation. For 61 events with an associated ground-truth counterpart in the ICMECAT catalog, forecasts based on initial MO data already achieve performance comparable to full-event reconstructions. Typical errors are ~5 hours in timing of magnetic field extrema and ~10 nT in field strength metrics, with limited systematic improvement as more of the event is observed. Results show substantial event variability and systematic underestimation of extrema, indicating deviations from ideal flux rope assumptions. These findings demonstrate the potential of fully autonomous real-time forecasting while highlighting limitations imposed by event complexity and model representational capacity.

Mapping the physical conditions of the shocked plasma of young supernova remnants (SNR) is crucial for understanding their explosion mechanisms, ejecta structure, and large-scale asymmetries. Using $>350$ ks of XRISM/Resolve high spectral resolution observations of Cassiopeia A (Cas A), the youngest known Galactic core-collapse SNR, we present the first microcalorimeter-based plasma parameter maps of any SNR. We tessellate Cas A into $1'\times1'$ regions and fit the broadband spectra as thermal emission from two pure-metal ejecta components -- corresponding to intermediate-mass elements (IMEs) and iron-group elements (IGEs) -- plus nonthermal synchrotron radiation. For robust inference, we introduce UltraSPEX, a Bayesian framework that couples the SPEX plasma code with the UltraNest nested-sampling algorithm, yielding full posterior distributions and exploration of parameter degeneracies. Key findings include enhanced Ar/Si and Ca/Si abundance ratios near the base of the Si-rich jets, and a high Ni/Fe mass ratio ($0.08\pm0.015$) in the base of NE jet. IGEs ejecta exhibit systematically higher Doppler velocities and broadenings than IMEs ejecta in most regions, with maximum differences of $\sim800$ km/s and $\sim1200$ km/s, respectively; Ca shows distinct (faster) kinematics from other IMEs in several SE regions. The ionization timescale and electron temperature show a robust anti-correlation, particularly for IGEs. This relation and measured parameter values could be explained by semi-analytical models with significant ejecta clumping (overdensities of $\sim10$ for IGEs and up to $\sim100$ for IMEs) and reduced historical reverse-shock velocities. Nonthermal emission accounts for a substantial fraction, with at least 47% of the 4--6 keV continuum and dominates in the western regions, where the spectrum hardens.

Compact astrophysical objects, such as neutron stars and white dwarfs, can act as detectors of energetic particle fluxes originating from astrophysical accelerators. While most existing capture and heating calculations assume isotropic very low energetic incident fluxes from the halo dark matter, many realistic sources produce highly directional beams or jets, for which gravitational focusing, trajectory multiplicity, and local energy deposition must be treated consistently. In this work, we develop a general relativistic formalism to compute the local density, capture probability, and energy deposition of particles arriving as directed beams onto compact objects. The framework is based on the mapping of an asymptotic particle flux to local densities through geodesic congruences, allowing for gravitational focusing, multi-stream regions, and optical depth effects to be incorporated in a unified way. The formalism applies to arbitrary particle species and interaction models, and separates capture from through-going energy deposition in a frame-consistent manner. As an explicit application, we consider relativistic particle beams generated in astrophysical jets and evaluate their interaction with two compact objects samples: a white dwarf and a neutron star. In particular, we illustrate the framework using boosted dark matter produced in a list of 324 blazars as a representative case study, computing the resulting fluxes and the associated heating in the selected stars. Additional regimes such as the interaction roof and geometric limit are discussed, highlighting the conditions under which compact objects can efficiently convert incident beam energy into observable heating.

Stars and other macroscopic objects may be surrounded by potentially large field configurations of very light scalars coupled to ordinary matter. If the star ends in a black hole, e.g. via a supernova or a neutron star merger, the source vanishes, and the field is released. In this paper, we improve on previous estimates for the field configurations arriving at large distances by including the effects of general relativity and an improved modelling of the initial field configurations. The total amount of energy released is typically of the same order of magnitude as suggested by simple flat space estimates. The spectrum receives noticeable corrections.

We study the impact of primordial black holes (PBHs) on the chiral phase transition and its associated stochastic gravitational-wave (GW) signals. Using the three-flavor Nambu-Jona-Lasinio model, we construct the chiral effective potential in a Schwarzschild spacetime background. We find that PBHs promote chiral symmetry restoration and induce a nontrivial local phase structure in the vicinity of the event horizon simultaneously. In particular, this structure exhibits a novel chiral symmetry breaking pattern involving both second- and first-order phase transitions, a feature absent in flat spacetime. We further demonstrate that PBHs act as genuine catalysts for the chiral phase transition by analyzing the bounce solution in curved spacetime. The presence of PBHs substantially enhances the inverse duration parameter $\beta/H$ while leaving the overall transition strength comparable to that in flat spacetime. As a consequence, even a small population of PBHs can induce $\mathcal{O}(1)$ shifts in both the peak frequency and the peak amplitude of the GW spectrum generated by the first-order dark chiral phase transition.

We construct a gravastar configuration within the Shtanov-Sahni (SS) braneworld scenario, characterized by a timelike extra dimension and negative brane tension. Unlike classical black holes, which inevitably culminate in central curvature singularities, our model demonstrates that the SS braneworld dynamics naturally regularize the interior geometry and prevent singularity formation. By solving the modified Einstein field equations induced on the brane, we obtain explicit interior, shell, and exterior solutions without invoking the idealized thin-shell approximation. The gravastar core is modeled as a Bose--Einstein condensate, while the intermediate shell consists of ultra-dense stiff matter. Bulk Weyl corrections induce anisotropic effective pressures on the brane, a feature that emerges intrinsically in this scenario and supports stability. We analyze the active gravitational mass, energy, entropy, and proper thickness of the shell, and establish the junction conditions at the interfaces. Our analysis reveals that the SS gravastar exhibits suppressed or even negative effective mass, reflecting the repulsive nature of the interior condensate, and admits stable equilibrium solutions consistent with energy conditions. This highlights the SS braneworld gravastar as a physically viable compact object and a compelling alternative to black holes. A key novelty of our construction is that the stabilizing pressure anisotropy and suppressed effective gravitational mass arise dynamically from higher-dimensional Weyl corrections, rather than being imposed through ad hoc matter sources or thin-shell idealizations. This provides the first fully analytic realization of a finite-thickness, stable gravastar in the Shtanov-Sahni braneworld, highlighting a genuinely geometric mechanism for singularity avoidance in compact objects.

We perform a Bayesian inference of the dense-matter equation of state (EOS) within a unified framework that incorporates hadronic matter, quark matter, and a smooth hadron--quark crossover. The EOS is constrained using physical consistency filters, gravitational-wave data from GW170817, NICER mass--radius measurements, and hypothetical future high-precision radius data. We find that current observations strongly constrain the density dependence of the nuclear symmetry energy, particularly its slope and curvature, while the highest-density hadronic parameters and quark-matter parameters remain only weakly constrained. The posterior distributions favor a crossover centered at an energy density $\overline{\varepsilon}\sim(4$--$6)\varepsilon_0$ with a width $\Gamma\sim(0.5$--$1.0)\varepsilon_0$, where $\varepsilon_0$ is the energy density of symmetric nuclear matter at saturation. The most probable radius of a canonical neutron star is $R_{1.4}\simeq 11.5-13.0$ km, and the maximum mass is $M_{\rm TOV}\simeq 2.2\,M_\odot$. Overall, present data primarily probe the low-to-intermediate density EOS and provide limited direct sensitivity to quark matter and genuinely high-density physics, highlighting the need for next-generation precision radius measurements.

The Spin Light of neutrino ($SL\nu$) is an electromagnetic radiation of the neutrino magnetic moment emitted when neutrino moves in external conditions (fields or matter). The effect can be of significance in the extremely dense matter of compact astrophysical objects such as neutron stars (NS). If detected, this radiation could provide a fair opportunity to study the properties of neutrinos and the medium through which they move, since the properties of the radiation depend on both. Motivated by the possibility of the nuclear matter spin-polarization, in this paper, we study the new properties to $SL\nu$ obtained under the influence of net matter polarization. We demonstrate that the polarization can enhance or completely suppress the radiation. Also, it introduces a characteristic asymmetry into the total radiation from the compact object, which could be an observable feature dependent on the matter polarization and the magnetic field inside the stellar (if the field is connected to the stellar matter polarization). The research may have implications for the physics of NS and magnetars, bringing us closer to the possibility of studying their internal structure.

Tilt-to-length (TTL) noise induced by angular jitter of spacecraft and test masses can affect the sensitivity of space-based gravitational-wave detectors such as LISA, Taiji, and TianQin. Such angular jitter can be measured using the differential wavefront sensing technique, enabling the modeling and subtraction of TTL noise from the data. However, owing to the multiple degrees of freedom of the detector constellation, a linear TTL model requires at least 24 parameters, while a higher-fidelity quadratic model involves up to 60 coefficients, rendering parameter estimation computationally expensive. To accelerate parameter determination, we propose a modified parameter set obtained via a linear transformation of the original angular coupling coefficients, which effectively reduces correlations among TTL noise components. In addition, we perform parameter fitting using an alternative second-generation time-delay interferometry configuration, PD4L, rather than the fiducial Michelson configuration. These two improvements enhance the convergence speed of the fitting procedure by a factor of approximately 10 for the linear model and approximately 18 for the quadratic model. The proposed approach can therefore substantially improve the efficiency of TTL noise calibration in space-based gravitational-wave detectors.

Visualizations have become an indispensable part of the scientific process. A vibrant ecosystem of visualization tools exists, catering to a wide variety of different needs. Real-time visualizations of numerical simulations offer scientists immediate feedback about the status of their simulations and can also be valuable educational and outreach tools. Developing a visualization tool with support for different operating systems, CPU/GPU architectures, and programming languages can be a challenge. It is common to use one or more graphics or UI libraries to act as abstraction layers and hide the underlying complexity. Whereas external libraries greatly simplify the initial programming effort, we argue that relying on them introduces new dependencies and problems, such as a higher entry barriers for new developers and users, and uncertainty regarding long-term support. In this paper we present a new approach for real time visualizations which we have implemented for the N-body package REBOUND. We propose to use a web browser to handle GPU accelerated rendering. This enables us to offer 3D, interactive visualizations on all major operating systems. What makes our new approach unique is that we achieve this without the need for any external libraries. We utilize WebAssembly to reuse existing OpenGL visualization code. Using communication via HTTP and a custom built-in web server, we are able to provide both local and remote real time visualizations. In addition to the browser based real time visualization, our approach offers other additional operating modes, including simulations running entirely within the browser, visualizations within jupyter notebooks, and traditional standalone visualizations using OpenGL. We focus on the implementation in REBOUND but the concepts and ideas discussed can be applied to many other areas in need of scientific and non-scientific real time visualizations.

In this work, we present the first extension of the Many-Body Forces (MBF) Model to finite temperature. The MBF Model describes nuclear matter in a relativistic quantum hadrodynamics formalism that takes many-body forces into account through a field dependence of the nuclear interaction coupling constants. Assuming nuclear matter to be charge neutral, beta-equilibrated, and populated by the baryon octet, electrons, and muons, we explore the parameters of the model, three different hyperon coupling schemes (also introduced here for the first time in MBF), and temperature effects to describe basic properties of nuclear matter, including the speed of sound, compressibility, and adiabatic index. We also investigate the mass-radius relation of compact stars by solving the Tolman-Oppenheimer-Volkoff equations at zero and finite temperature, including scenarios with fixed entropy per baryon. Our original results at finite temperature open the path to a new description of proto-neutron stars.

We demonstrate that dark matter axion production is enhanced through a natural and unavoidable mechanism: primordial temperature fluctuations periodically modulate the axion mass during the QCD phase transition, thereby triggering parametric resonance in axion field evolution. This interplay between parametric resonance and the misalignment mechanism shifts the predicted axion mass window for the observed dark matter abundance to $10^{-4}-10^{-3} \, \text{eV}$, displacing the canonical axion mass window to previously unexplored higher ranges.

We present a framework for systematic computations of scattering amplitudes for gravitational Raman scattering, -- the inelastic scattering of massless fields off compact relativistic objects. We focus on the small-frequency (post-Minkowskian, PM) regime relevant for the study of tidal effects, which can be mapped onto gravitational wave observables during the inspiraling phase of a merger. We demonstrate that this setup is ideal for systematic studies of tidal effects, in a way that is free from coordinate, gauge, and field redefinition ambiguities. We use a combination of worldline effective field theory, the background field method, and advanced scattering amplitude techniques to derive phase shifts for scattering of spin-$0,1,2$ fields off generic compact objects at third PM order. We demonstrate that the inclusion of the recoil of the object is crucial for consistency of this calculation. Focusing on a particular case of black holes, we extract the leading static and dynamical Love numbers of the spin-0 field and the static Love number of the spin-1 field in four dimensions by matching our EFT amplitudes and calculations in General Relativity. We show, fully on-shell, that the leading static Love numbers vanish identically, while the dynamical Love numbers are not zero and run logarithmically. The latter resolves the ambiguities of previous off-shell matching calculations. We also extend our results to seven dimensions, where spin-2 Love numbers undergo a renormalization group running at 2PM, which we compute explicitly. In addition, we extract the leading static Love numbers of spin-0 and spin-1 fields in five dimensions, which also run.