Papers for Friday, Mar 27 2026

Superluminous supernovae (SLSNe) are rare transients that are $\sim 10 - 100$ times more luminous than ordinary stellar explosions, reaching peak optical luminosities $\sim 10^{44} - 10^{45}$ erg s$^{-1}$. The energy source powering SLSNe remains uncertain. In this work, we explore the multi-wavelength and multi-messenger signatures of the scenario in which SLSNe are powered by a newly born millisecond magnetar. We model the dynamical evolution and emission from the coupled system comprised of the magnetar, wind, nebula, and supernova ejecta, consistently evaluating the pair multiplicity of the wind and nebula regions, and the bulk wind Lorentz factor governing the $e^+ - e^-$ injection spectra in the nebula. We compute the thermal and non-thermal electromagnetic signatures, neutrino signatures, and investigate their detection prospects. For SN 2017egm, the nearest observed SLSNe, our prediction for high-energy gamma rays matches the recent detection by Fermi LAT. For neutrinos, using SN 2017egm a canonical SLSNe, we find that in the era of the Vera C. Rubin Observatory, a stacking analysis with upcoming neutrino observatories can lead to $3\sigma$ detection significance of neutrino events from a population of SLSNe within a decade of operation.

We present BLAZ4R, the first living catalog of confirmed $z>4$ blazars, with a focus on the contribution of eROSITA, on board of the Spectrum Roentgen Gamma (SRG) spacecraft. Blazars at $z>4$ are rare but powerful probes of how active supermassive black holes evolve in connection to relativistic jets, in the first 2 billion years of cosmic history. At these redshifts, X-ray observations are essential for constraining blazars jet power and orientation, enabling effective trace of their parent population. The all-sky surveys provided by eROSITA ensure X-ray detection for BLAZ4R sources and, in some cases, allow spectral and temporal studies of their jetted emission. BLAZ4R includes 54 confirmed blazars, characterized through their X-ray properties, radio spectra and morphology, and multiwavelength profiles. We confirm that jetted sources are significantly more numerous relative to non-jetted counterparts at high-$z$, and that blazars (and therefore the overall jetted population) do not exhibit significantly different features compared to the total active galactic nuclei population in the early Universe. Fast accretion processes that involve relativistic jets are clearly required to justify the existence of fully formed jetted AGN at $z>4$. However, the diverse multiwavelength properties characterizing BLAZ4R do not yet allow us to identify the specific signatures of these processes. We will continue updating BLAZ4R to search for such signatures and ultimately understand the early formation of jetted AGN.

Variability monitoring provides an unparalleled insight into the atmospheric processes of brown dwarfs and directly imaged exo-planets. Inhomogeneous clouds, aurorae and magnetic spots have all been postulated as potential drivers of variability. While objects at the L/T transition have had their variability studied extensively, the variability of early L-dwarfs remains an understudied region of the parameter space. We use observations from the Hubble Space Telescope in the near-infrared, using WFC3/G141 to disentangle the drivers of variability in three known variable early L-dwarfs: 2MASS J1721039+334415, 2MASS J00361617+1821104 and 2MASS J19064801+4011089. We find that all three objects exhibit significant variability at all wavelengths, with white-light amplitudes of 0.53-1.41 %. We find that their colour variations are brighter and bluer compared to later spectral types, except for 2MASSJ19064801+4011089 which exhibits largely grey variations. We report a new period for 2MASS J1721039+334415, of 4.9^{+0.4}_{-0.2} hours. We find evidence of long term light curve stability in each object, which may indicate the presence of long lived features on their surfaces. We create a flexible modelling framework to model three potential drivers of variability: clouds, aurorae and magnetic spots. We fit our models to the spectral variability amplitude from 1.1-1.67 {\mu}m of each object. We find that changing cloud properties or magnetic spots are the most likely drivers of variability in each object. Auroral models do not reproduce the variability within the HST wavelengths, however future observations at longer wavelengths that probe higher in the atmosphere may be more sensitive to auroral effects. This work provides a foundation for future variability studies of early L-dwarfs and directly imaged exoplanets to disentangle auroral, cloud and magnetic spot driven variability.

We present a probabilistic autoencoder (PAE) framework for galaxy spectral energy distribution (SED) modeling and redshift estimation, applied to synthetic SPHEREx 102-band spectrophotometry. Our PAE learns a compact latent representation of rest-frame galaxy SEDs transformed to a simple Gaussian base density using a normalizing flow, combined with an explicit forward model enabling joint Bayesian inference over intrinsic SED parameters and redshift with well-defined priors. In controlled tests on simulated SPHEREx spectra, our PAE improves on template fitting (TF) in source recovery, outlier rate, and posterior calibration, with trade-offs in redshift performance that depend on the assumed priors. A simple cut on the ratio of PAE and TF uncertainties identifies sources that are overwhelmingly TF outliers, which can be used to clean existing TF samples while retaining the vast majority of well-recovered sources. By directly profiling over PAE latent variables, we show these cases correspond to shallow likelihood surfaces where the PAE's continuous SED manifold produces broader likelihoods that more faithfully reflect the lack of constraining power in the data, whereas the TF discrete model grid yields artificially confident but incorrect redshift estimates. Lastly, we present an alternative, simulation-based inference approach using a Transformer encoder and conditional normalizing flow, which provides similar redshift performance to the PAE but with $\sim200\times$ faster inference throughput. Our implementation, \texttt{PAESpec}, is publicly available and provides a foundation for principled redshift estimation in modern photometric surveys.

We present the results of a multi-cycle Chandra program to systematically monitor the X-ray variability of 10 weak-line quasars (WLQs) that previously had limited multi-epoch X-ray observations. Three new Chandra 2.8 to 8.2 ks observations were obtained for each WLQ with C$\,$IV rest-frame equivalent widths (REWs) $\lesssim 10$ Å, substantially improving the monitoring data quality of WLQs and our ability to characterize their long-term X-ray variability behavior. We observe recurrent extreme X-ray variability in the historically variable WLQ SDSS J1539+3954, with an X-ray flux rise of a factor of $\gtrsim 6$ between 2023 and 2024 ($\gtrsim 21$ relative to 2013). Another previously X-ray weak WLQ in the sample, SDSS J0825+1155, underwent a significant X-ray flux variation by a factor of $\gtrsim 14$ between 2019 and 2023. We find the fraction of WLQs exhibiting evidence of extreme X-ray variability to be $0.20^{+0.17}_{-0.07}$. In the context of the thick disk and outflow (TDO) model, the substantial fraction of WLQs displaying extreme X-ray variability may suggest that the variability is driven by the intrinsic motion of the TDO wind rather than changes in the height of the TDO disk. We performed a statistical comparison between the distribution of variability amplitudes of WLQs and general radio-quiet quasars. We find that these underlying distributions are statistically different, with WLQs having $\approx 6.8$ times higher odds of exhibiting an extreme X-ray variability event than the general radio-quiet quasar population.

Galaxies with intense star formation often host multiphase, galaxy-scale winds powered by supernovae and fast stellar winds. These are strong enough to disrupt the star-forming interstellar medium, and they chemically enrich the surrounding circumgalactic medium. However, their launching mechanism remains unknown. Here we show that thermal gas pressure is sufficient to drive the multiphase wind in the prototypical starburst galaxy M82. Using a high energy-resolution ($\Delta E = 4.5$ eV) XRISM Resolve spectrum, including detections of FeXXV 6.7 keV, ArXVII 3.1 keV, and SXVI 2.6 keV, we measure the temperature ($T = 2.3^{+0.5}_{-0.2} \times 10^7$ K) and mass ($M \approx 6 \pm 2 \times 10^5$ M$_\odot$) of the hot gas in the starburst and provide the first direct measurement of its line-of-sight velocity dispersion ($\sigma = 595^{+464}_{-128}$ km s$^{-1}$). These values are consistent with a freely-expanding wind exceeding the galactic escape velocity. The size of the FeXXV-emitting region suggests a hot gas outflow rate of $\dot{M} \approx 4$ M$_\odot$ yr$^{-1}$, carrying a total energy of $\dot{E} \approx 2 \times 10^{42}$ erg s$^{-1}$. This is sufficient to drive the molecular, atomic, and ionized outflows while transporting up to $\approx 2$ M$_\odot$ yr$^{-1}$ of hot gas to the intergalactic medium. The estimated supernova rate implies that $\approx$ 60% of the supernova energy must be thermalized in hot gas. Our results suggest that additional driving mechanisms, such as cosmic-ray pressure, are not required to launch the wind.

The largest Lyapunov exponent is widely used to diagnose chaos in gravitational dynamics, but in mixed phase spaces and finite-N systems it does not always provide a complete description of orbital complexity and phase-space transport. Entropy-based diagnostics may offer a complementary perspective. We investigate whether trajectory-based information entropy can provide a useful diagnostic of chaos in gravitational systems and how it relates to the largest Lyapunov exponent as a function of orbital energy and of the number of degrees of freedom. We computed the largest Lyapunov exponent and a coarse-grained Shannon entropy for ensembles of trajectories in the Henon-Heiles potential and for test-particle orbits in live N-body realizations of a Plummer model. We then compared the dependence of both quantities on orbital energy and, for the N-body case, on particle number. In the Henon-Heiles system, the Shannon entropy follows the transition from weak to widespread chaos and exhibits an energy dependence that closely mirrors that of the largest Lyapunov exponent. For test-particle orbits in live N-body potentials, both diagnostics indicate stronger chaos for more tightly bound trajectories. However, their dependence on N differs: the largest Lyapunov exponent remains nearly constant over the explored range of particle numbers, whereas the Shannon entropy decreases monotonically as N increases. These results show that the information entropy can complement the largest Lyapunov exponent and may better capture changes in global phase-space mixing, especially in systems where the leading Lyapunov exponent alone is not sufficiently informative. It therefore provides a promising alternative for diagnosing chaos when tangent-space dynamics is unavailable or computationally expensive, and it is naturally suited to systems with densely sampled trajectories, such as minor bodies in the Solar System.

We establish the multi-messenger mechanics of episodic mass transfer in supermassive black hole binaries stalled within circumbinary discs. Utilizing continuous wavelet transforms, we isolate localized gas clumps at the cavity edge and track their evolution. By regularizing the forced fluid equations at Lindblad resonances via the inhomogeneous Airy differential equation, we bypass linear singularities to extract the finite wave amplitudes that trigger non-linear shock formation. These shocks produce bounded accretion bursts. We model the time-domain thermal luminosity, deriving an analytical power spectral density that forms a harmonic cascade. The superposition of the accretion streams generates a spectral beat frequency, providing an exact mathematical extraction of the binary mass ratio. The radiative cooling of shock-accelerated electrons produces a multi-wavelength spectral energy distribution from a synchrotron radio continuum to an inverse-Compton gamma-ray tail. We identify a relativistic signature: a discontinuous, high-frequency gravitational wave sideband termed the ``background gas humming''. This emission arises from the highly asymmetric, transient fluid geometry of the accretion shocks. Evaluating the asymptotic properties of the Airy regularization, we show that this humming manifests as a sequence of discrete high-frequency bursts with temporal quiescence gaps that systematically compress as the cavity shrinks. We show that the instantaneous mass of the gas actively trapped within the cavity violently amplifies prior to decoupling, culminating in a terminal burst near 4.0 mHz that serves as a multi-messenger precursor to the final vacuum inspiral.

We use WIYN/NEID radial velocity measurements to confirm the planetary nature and measure the mass of the TESS transiting exoplanet candidate around the mid-K dwarf GJ 523 ($V=9.23$, $K=6.525$). We find that GJ 523b is on a 17.75 day orbit and has a radius of $2.55\pm0.15\,R_\oplus$, a mass of $23.5\pm3.3\,M_\oplus$, and a zero-albedo equilibrium temperature of 538 K. GJ 523b's high bulk density of $7.8\pm1.8$ g cm$^{-3}$ and position on a mass-radius diagram implies a surprising low atmospheric mass fraction despite its relatively large mass. Additionally, we determine that the system has an age of $169^{+100}_{-48}$ Myr through a gyrochronological analysis of GJ 523 and its comoving companions. We also use the SED-derived stellar radius, the photometric rotation period, and the spectroscopic $v\sin i_\star$ to derive a stellar inclination of $17.6\pm5.0$ degrees, implying that GJ 523b has a minimum orbital obliquity of $71.4_{-5.0}^{+4.7}$ degrees. GJ 523b's high mass, apparent lack of a gas envelope, young age, and high orbital obliquity present a challenge to typical planet formation pathways, and at the moment there is not enough data on the system to definitively determine how GJ 523b formed. Finally, we present a new observational classification for ultra-dense, sub-Neptune-sized exoplanets similar to GJ 523b: the mega-Earths, planets with $R_p \geq2.1\,R_\oplus$ and $\rho_p \geq 5.5$ g cm$^{-3}$.

Black hole X-ray binaries (BHXRBs) have traditionally been discovered by X-ray surveys with cadences of hours to days. However, large optical time-domain surveys now provide novel avenues for early detection and insights into their elusive outburst triggering mechanisms. We present early-time light curves of the BHXRB AT 2019wey serendipitously observed by the Transiting Exoplanet Survey Satellite (TESS). The TESS images are sampled at 30 minute cadence from $\approx2$ d prior to $\approx25$ d after outburst, providing the highest time resolution optical rising phase observations of any known BHXRB. We fit a piece-wise power law to the rising light curve, finding an outburst onset time of MJD $58817.86\pm0.09$ and power-law rise index $n=0.74\pm0.04$. The onset time precedes all ground-based optical detections, and suggests that the optical rise began after the start of the faint X-ray brightening in MAXI data. We search for periodic high frequency modulation and detect none exceeding amplitude $\approx0.48$ mJy at periods of $\gtrsim1$ h at 90% confidence.

We report on a new method for incorporating interactions between dark matter (DM) and baryons in cosmological simulations, capable of handling the challenging regime in which the dark matter particle mass is comparable to or lighter than the baryon mass. The method hybridizes two distinct approaches, as gas particles receive momentum and energy transfer according to a mean-field calculation while DM particles undergo Monte Carlo scatterings, which are derived from the Boltzmann equation and proved to be statistically equivalent. We present an open-source implementation of this method in the simulation code GIZMO. As a first application, we investigate the effects of DM--baryon interactions on an isolated Milky Way-like disk galaxy for dark matter having twice the proton mass, which roughly maximizes the average energy transfer per collision. For cross sections of order 1 barn ($10^{-24}\, \mathrm{cm}^2$), these interactions cause strong changes to the mass distribution in the center of the galaxy in less than 1 Gyr, even when bar formation is suppressed by hand. For cross sections typical of hadronic interactions, $\lesssim 30 \, {\rm mb}$, high-fidelity galaxy formation simulations will be needed to assess the effects on observable features of galaxies.

The relation between the maximum stellar mass in a very young cluster (mmax) and the total stellar mass of the cluster (Mecl), known as the mmax-Mecl relation, remains debated in the literature. To test the validity of this relation, we modelled young star clusters with masses between 102.5 and 105.0 M_sun and ages of 1-4 Myr using the galIMF code, in which stellar masses are optimally sampled from a varying initial stellar mass function. We compared the results with literature observations of extragalactic young star clusters. We incorporated stellar evolution via PARSEC and COLIBRI tracks and computed Halpha luminosities using the Pegase code. To account for dynamical ejections, we stochastically removed stars based on their spectral type, following previous N-body simulations. Additional sources of scatter, including uncertainties in age determination and contamination by field stars, were considered. Our results indicate that, under the assumptions explored here, optimal sampling is consistent with the extragalactic star cluster observations considered, whereas purely random sampling produces distributions that are not in agreement. These findings support a highly self-regulated interpretation of cluster formation in which stellar masses align optimally with the initial mass function rather than being drawn independently at random.

More than half of Solar-type stars are found in binary systems. The numbers of exoplanets within binary systems in s-type orbits now numbers over 700. However, whilst the numbers have increased, there still does not exist a global model of planet formation for wide binary systems, where there does for single stars and circumbinary systems. As a precursor to such a model, that includes the necessary physical planet formation processes, it is important to understand how an outer binary companion affects the evolution of circumstellar discs, and the formation of planetesimals and planetary embryos. The main mechanism for which these processes are affected, is through truncation of the protoplanetary disc outer edges. In this paper, we determine these effects, whilst also comparing them to the effects of external photoevaporation that competes to truncate protoplanetary discs. We find that disc truncation by both a binary companion and external photoevaporation significantly reduces the efficiency to which planetary embryos are able to accrete pebbles and grow into terrestrial mass planets. This is due to the pebble supply being cut off as the pebble production front reaches the disc outer edge before planets are able to significantly increase in mass. This hindrance to planet formation occurs when the truncation radius due to the binary companion is below $\sim 30$ au, corresponding to binary separations of $\sim90$ au for equal mass, circular binary stars. For separations greater than 300 au, planet formation operates similar to that around single stars. Our results highlight the detrimental effects of a binary companion for intermediate binaries, that can provide possible explanations for the dearth of multiple planets within binary systems of separations $<100$ au

Estimating galaxy redshifts is crucial for constraining key physical quantities like those in the equation of state of dark energy. Modern telescopes such as the James Webb Space Telescope, the Euclid Space Telescope, and the NASA Nancy Grace Roman Space Telescope are producing massive amounts of spectroscopic data that enable precise redshift estimation. However, a galaxy's redshift can be estimated only when emission lines are present in the observed spectrum, which is unknown a priori. A novel Bayesian approach to estimating redshift and simultaneously testing for the presence of emission lines is developed. Although modern spectroscopic surveys involve millions of spectra and give rise to highly multimodal posterior distributions, the proposed framework remains computationally efficient, admitting a parallelizable implementation suitable for large-scale inference.

The central supermassive black hole of the galaxy M87 is currently a target for precision spin measurement using high-resolution, horizon-scale imaging. Such observations aim to resolve the first lensed (${n}~{=}~{1}$) sub-image of the photon ring from the broader direct image. In this work, we identify a concrete observable -- the displacement between the centers of the ${n}~{=}~{1}$ photon-ring sub-image and the direct image -- and propose its use in a simple spin-measurement technique. Leveraging the assumption that the observed large-scale jet of M87 is aligned with the black-hole spin axis, we separate the relative position of the photon ring into components parallel and transverse to the projected spin axis, normalizing both components with respect to the measured diameter of the ${n}~{=}~{1}$ sub-image. We show that the parallel shift is primarily determined by inclination and emission radius, while the transverse shift is tightly correlated with inclination and spin. We demonstrate these effects both in a simple geometric model (to explain the underlying physics) and in GRMHD simulations with magnetically arrested disks (to provide realistic instantiations of the effect). We find that a relative astrometric resolution of ${\lesssim}~{0.1\;\mu\rm{as}}$ is sufficient to constrain the spin to better than 9% if the accretion flow is prograde or 22% if the flow is retrograde. If the direction of the accretion flow is undetermined, the spin can be constrained to within 26%. More generally, this identifies relative photon ring astrometry as a promising method to constrain the underlying spacetime geometry and introduces a spin-constraint technique that does not rely on geometric modeling of the observed emission.

Red quasars exhibit a higher incidence of compact (galaxy-scale or smaller) radio emission than blue quasars, arising from systems near the radio-loud/radio-quiet threshold. In this paper we select quasars from SDSS ($0.2 <z <2.4$), and use archival radio data (FIRST, VLASS, LoTSS) to visually determine the radio morphologies of 573 red quasars compared to a control sample of 1278 typical blue quasars. We find an excess of steep-slope radio emission ($\alpha_{1.4-3\text{ GHz}}\sim-1$, where $S_\nu \propto \nu^\alpha$) from red quasars with compact ($<6''$) radio morphologies over 144 MHz, 1.4 GHz, and 3 GHz. This excess steep radio emission signature is not seen in normal blue quasars (radio compact or extended) or red quasars with extended low-frequency radio emission, which instead show a broad range of radio spectral slopes consistent with a range of different physical processes. We show that the strength of the excess steep spectral slope component increases with dust extinction, along with an overall increase in the radio-detection fraction. We argue that this excess steep-slope radio emission is due to shocks between quasar-driven winds/jets and the dusty nuclear-host galaxy environment. The majority ($86^{+5}_{-21}\%$) of the dustiest quasars ($E(B-V)>0.4$ mag) with steep slopes have radio luminosities consistent with the prediction from a wind-shock model with wind efficiencies of up to 7%. This agrees with the scenario where these compact red quasars are undergoing a "dusty blow-out" phase, where compact jets and/or AGN-driven winds interact with a dusty ISM, causing shocks, leading to steep spectral slopes and enhanced radio detection rates.

For a decade there has been a factor of 2.5 gap in luminosity between the 275K WISE J085510.83-071442.5 (Luhman 2014) and all other Y dwarfs, with Teff >= 350K. Recently three objects were found which may fall in this gap. Two are companions to Y dwarfs: WISE J033605.05-014350.4B (Calissendorff et al. 2023) and CWISEP J193518.58-154620.3B (De Furio et al. 2025); the third is MEAD 62B, a candidate companion to a white dwarf (Albert et al. 2025). Evolutionary models calculate a tight relationship between luminosity and Teff for Y dwarfs. I determine luminosities and hence Teff for three Y dwarfs (WISE J085510.83-071442.5, WISE J173835.53+273259.0, WISE J182831.08+265037.7). I derive relationships between Teff and mid-infrared colors using these together with 22 T and Y dwarfs from Beiler et al. (2024) with luminosity-based Teff values. These relationships are used to explore the Teff distribution for Y dwarfs. A sample of 31 Y dwarfs within ~20 pc is presented with 275 < Teff K < 425. The JWST colors for WISE J053516.80-750024.9 and WISE J182831.08+265037.7 support previous suggestions that they are unresolved binaries, the former a 480K and 340K dwarf pair and the latter a pair of 387K dwarfs. Five other dwarfs have unusual colors; two are likely high gravity and/or metal-poor (WISE J024714.52+372523.5, WISEA J215949.54-480855.2), two low gravity and/or metal-rich (CWISEP J104756.81+545741.6, WISE J150115.92-400418.4), and the fifth cannot be interpreted (WISE J043052.92+463331.6). An Appendix provides colors which can be used as a reference for searches for brown dwarfs in JWST data.

Diffuse gamma-ray emission is a key tracer of cosmic rays (CRs) in galaxies, encoding information about their transport, energetics, and interactions with the interstellar medium. Interpreting the Milky Way gamma-ray sky is challenging because the observed emission depends on the three-dimensional distributions of CRs and gas, as well as the observer location within the Galaxy. Using the Rhea suite of CR-MHD simulations of a Milky Way analogue, we study how pion-decay gamma-ray emission varies with galactic environment, local conditions, and CR transport physics. Emission is computed in post-processing under steady-state assumptions, enabling analysis of luminosities, spectra, full-sky maps, and angular power spectra (APS) for multiple observer positions, including those inside Local-Bubble-like cavities. The simulated galaxy reproduces Milky Way-like gamma-ray luminosities and spectral slopes without parameter tuning. While total luminosity is relatively stable, the morphology of the gamma-ray sky varies strongly with observer location due to the complex local gas distribution, consistent with observations. For all observers, the APS traces gas column density rather than the smoother CR energy density, in agreement with previous studies. Comparisons with Fermi-LAT data show good agreement in both the all-sky spectrum and APS, with a diffusion coefficient scaling proportional to E^0.5 providing the best match. These results show that key features of Galactic gamma-ray emission arise naturally in self-consistent CR-MHD simulations. Gas density fluctuations primarily shape emission morphology, while CR transport governs spectral and structural properties. The Rhea simulations provide a physically grounded framework for interpreting diffuse gamma-ray observations and highlight the importance of local environment in tracing Galactic CR physics.

The merger between the carbon (C) and oxygen (O) shells hours to days before the collapse of a massive star significantly changes its nucleosynthesis, which is reflected in the elemental ratios observed in supernova remnants (SNRs). We present a nucleosynthesis study of $^{44}$Ti production in core-collapse supernovae (CCSNe), highlighting large silicon (Si), sulfur (S), calcium (Ca), and, most importantly, argon (Ar) to neon (Ne) ratios as diagnostics for carbon-oxygen (C--O) shell mergers. We compare yields from eight different sets of CCSNe models to observations of Cassiopeia A (Cas A), and show that C--O shell mergers are consistently the models that best match X-ray and infrared observations. These models produce high Ar/Ne ratios ($\gtrsim 0.1$), due to $^{20}$Ne depletion and production of $^{36}$Ar and $^{38}$Ar, while lower ratios are obtained from non-merger cases. Based on the Ar/Ne diagnostic, we compare the range of expected $^{44}$Ti produced by C--O shell mergers, which is up to $\sim 20 - 30 \%$ of the overall $^{44}$Ti, but expected to be located outside the reverse shock. Based on the sets of models considered, the photon flux expected from the $^{44}$Ti synthesized in the C--O shell merger in Cas A is below the $NuSTAR$ and $COSI$ detection limits, compatible with current limits locating most of the $^{44}$Ti interior to the reverse shock, but might be detectable from proposed missions like $ASCENT$. Finally, for the SNR of 1987A, a dominant C--O merger origin of the observed $^{44}$Ti is unlikely based on the observed redshift in its $^{44}$Ti line.

{Supernova (SN) siblings are powerful tools used to calibrate and improve distance measurement methods, and to make the systematic uncertainty to distances to their host galaxies considerably lower compared to other techniques.} {In this paper we present distance estimates to NGC6951, a galaxy that hosted the Type IIP SN~2020dpw, the Type Ib SN~2021sjt, and three other SNe.} {Photometric observations of the two objects were carried out using two 80cm RC telescopes located in Hungary, while spectra were obtained from the LCO and the WiseRep database. For the distance estimates, we applied the expanding photosphere method (EPM), which connects the observed angular radius ($\theta$) of a SN to its physical radius and is related to the velocity of the photosphere ($v_{\rm ph}$). Although the EPM is mostly applied to derive the distance of Type IIP SNe, in the literature there are several examples of this technique being used for Type IIn and stripped-envelope SNe as well. Therefore, we made another attempt to infer the distance of the Type Ib SN~2021sjt by applying the EPM together with its Type IIP sibling SN~2020dpw. } {Our analysis resulted in a distance of $25.76 \pm 0.34 (\rm random) \pm 5.51$ (systematic) Mpc and $24.57 \pm 1.27 (\rm random) \pm 4.64$ (systematic) Mpc for SN~2020dpw and SN~2021sjt, respectively. Systematic errors were estimated with respect to the used dilution factor, the interstellar reddening, and the date of the explosion (which was fixed to a value between the last non-detection and the first detection for each object).} {The obtained distance values agree with each other and with the literature, which shows the validity of the methods used. In this way, new and perhaps improved distance estimates to NGC 6951 were obtained, and the applicability of the EPM for Type Ib SNe was tested.}

The thermal evolution and interior structure of giant exoplanets are sensitive to the treatment of radiative opacity. At temperatures of ~2000 K, depletion of alkali metals can create a window of reduced opacity, potentially giving rise to deep radiative zones. While such zones have been discussed for Jupiter, their role in the evolution and characterisation of warm giant exoplanets has not been systematically investigated. We investigate how opacity windows and the resulting deep radiative zones affect the cooling, radius evolution, and the characterisation of interiors and atmospheres of giant exoplanets. We computed thermal evolution models for warm Jupiters spanning masses of 0.3 to 1.0 Jupiter masses with equilibrium temperatures of 200 to 800 K, with a parametrised reduction of the radiative opacity near ~2000 K. Deep radiative zones develop in moderately irradiated Jupiters older than ~4 Gyr even with unmodified opacities, and earlier and more extensively when the opacity is reduced. A deep opacity window accelerates the planetary cooling, reducing predicted radii by up to 5% and interior temperatures on the order of a few 10%. We show that this translates to a ~10 percentage point difference in inferred bulk metallicity. Deep radiative zones are likely common in warm giant exoplanets and could decouple atmospheric composition from bulk interior composition, complicating the interpretation of atmospheric observations. We suggest that the opacity treatment introduces significant uncertainties in atmospheric and interior characterisation.

The interiors of giant planets are commonly assumed to be convective and therefore adiabatic, making the adiabatic temperature gradient a key ingredient in interior and evolution models. While there are multiple numerically distinct ways to compute this gradient, their impact on inferred planetary structure and radius has not been systematically assessed. We investigate how the numerical treatment of adiabatic temperature profiles affects inferred planetary radii and internal structure, quantifying the impact of different methods for calculating the adiabatic gradient and different forms of the temperature differential equation on static interior models. We computed static interior models of a one Jupiter mass planet using a state-of-the-art hydrogen-helium equation of state, comparing five methods for evaluating the adiabatic gradient against a ground-truth isentropic baseline, for both the logarithmic and non-logarithmic forms of the temperature equation. The choice of numerical method significantly impacts the inferred interior structure and the radius. With the logarithmic temperature equation, central temperatures deviate by several thousand Kelvin and surface radii differ by up to 3.4%, exceeding the ~1% precision of current giant exoplanet radius measurements threefold. The non-logarithmic form reduces deviations to below ~1% for most methods. We therefore recommend spline derivatives to evaluate the adiabatic gradient via the triple-product rule, combined with the non-logarithmic temperature equation. Finite differencing and direct use of tabulated gradients or derivatives should be avoided.

We investigated if signatures of coronal mass ejections (CMEs) can be retrieved in the differential emission measure (DEM) from Sun-as-a-star extreme-ultraviolet (EUV) observations. We analyzed 16 major, eruptive (i.e., associated with CMEs) flares. For each flare we constructed light curves of the average intensity of full-disk images taken by the SDO/AIA mission in EUV channels centered at 94, 131, 171, 193, 211, and 335 Angstrom. We also corrected the light curves for the flare gradual phase. From the analysis of the light curves we find that all the studied flares exhibit dimmings, where the intensity decreases with respect to the pre-flare phase, mainly in the 171, 193, and 211 Angstrom channels. The dimmings in these channels become more pronounced upon applying the gradual-phase correction. Calculation of the DEM from the six AIA EUV channels shows that during all the observed dimmings, the DEM decreases with respect to its value in the pre-flare phase in the temperature range $10^{5.7}-10^{6.3}$ K. The signature of the dimming is more pronounced in the range $10^{5.7}-10^{6.0}$ K for the DEMs calculated with the original light curves, and in the range $10^{6.0}-10^{6.3}$ K for the DEMs calculated by taking into account the gradual-phase correction. For a sample event, we also calculated DEMs from EVE and spatially resolved AIA observations of the source region to assess the impact of spectral resolution and full-disk averaging. For both these cases the temperature range where the dimming in the DEM is more pronounced is similar to that resulting from the analysis of the spatially averaged AIA data; the magnitude of the dimming is similar for the EVE and larger for the spatially resolved AIA observations. Coronal dimmings associated with CMEs can be detected in Sun-as-a-star DEMs. The flare gradual phase can lead to an underestimation of the magnitude of the dimming.

We aim to search for satellites and circumplanetary or circumsubstellar disks around directly imaged substellar companions, exploring their immediate environment to constrain the conditions for satellites and disk formation. We conducted a dedicated survey of twelve planets and brown dwarfs with VLT/SPHERE using a novel application of the star hopping technique. By building libraries of contemporaneous point spread function (PSF) references from nearby stars, we applied a frame-by-frame subtraction of each companion's flux using the Negative Fake Companion method (NEGFC). This approach mitigates temporal PSF variability and enhances sensitivity to faint circumplanetary features. We derived contrast curves, translated them into mass detection limits using evolutionary models, and constrained the dynamically stable regions through estimates of Hill radii from orbital fits. Our analysis yields stringent limits on the presence of massive satellites, generally excluding companions more massive than a few Jupiter masses at separations beyond 1-5 au, depending on each system's Hill radius. In most cases, no convincing point-like or extended residuals were found. However, we identify promising signals for three systems: extended residuals consistent with a circumplanetary disk around CT Cha b, tentative repeated residuals near TYC 8047-232-1 B that may trace a bound satellite companion of 3-6 MJup, and marginal residual signals at the location of the previously reported candidate around DH Tau b, whose interpretation, however, remains uncertain due to possible contamination by instrumental effects. These results confirm the power of star hopping in reducing PSF-related artifacts and provide some of the most stringent constraints to date on the mass and location of potential satellites and disks around directly imaged substellar companions.

Magnetospheric accretion is the paradigm for accretion in Classical T-Tauri Stars (CTTS). However, the standard, one-flow magnetospheric accretion model fails to replicate important characteristics such as the observed Balmer decrements. We address this limitation by adopting a model with two axisymmetric magnetospheric accretion flows of different accretion rates and geometries. We calculate the fluxes of the hydrogen $H_\alpha$, $H_\beta$, and $H_\gamma$ lines of each flow with the magnetospheric accretion model and use Bayesian statistics to fit the Balmer line fluxes of 139 CTTS in the Orion OB1b subassociation, and in the Upper Scorpius, Lupus and Chamaeleon I star-forming regions. We find that the Balmer decrements and line fluxes can be fitted by two distinct but coexisting flows: a compact, high accretion rate flow, close to the star and narrow (mean inner radius $R_i \sim 2.9 R_*$ and mean width $\Delta R \sim 0.7 R_*$), covering a few percent of the emitting area, and a more spread out flow, thicker ($\Delta R \sim 1.2 R_*$), and larger ($R_i \sim 3.7 R_*$), with lower accretion rate, encompassing the rest of the emitting area. The two-flow model can also reproduce the empirical correlation between the luminosity in $H_\alpha$ and the accretion luminosity. Overall, our findings suggest that a multicolumn approach provides a more accurate representation of the observed Balmer line emission, in agreement with results of numerical simulations.

With the advanced capabilities of ground- and space-based telescopes, exoplanet science is beginning to characterize the physics and chemistry of exoplanet atmospheres. However, interpreting exoplanet observations requires sophisticated modeling tools that rely on laboratory data to provide critical inputs and constraints. In preparation for the influx of observational data that the coming decades will bring, laboratory experiments that simulate the diverse conditions expected in exoplanet atmospheres, surfaces and interiors are vital to advancing models and thereby our understanding of these worlds. Here we discuss the key areas where laboratory experiments are providing fundamental data for exoplanet models, particularly for low-mass planets from rocky worlds to sub-Neptunes. First, we present a series of experiments to measure outgassing and volatile solubilities that are essential for establishing the connection between low-mass exoplanet interiors and their observable atmospheres. We then discuss additional laboratory techniques that can be used to understand the interior-atmosphere connection and simulate the high pressure-high temperature conditions of exoplanet interiors. Next, we summarize the experimental methods used to constrain the spectroscopic properties of atmospheric gases and aerosols along with their formation and reaction mechanisms. We also discuss how similar methods can be used to constrain exoplanet surface compositions, which is important for interpreting observations of atmosphere-less worlds. Finally, we conclude by presenting several examples of astrobiology experiments that constrain how life can modify the atmosphere and surface of rocky exoplanets. Together, these laboratory efforts are crucial to maximizing the scientific yield of exoplanet observations in the coming decades.

Several theoretical and observational studies have shown that new waves of triggered star-formation can be induced by the feedback from newly formed massive protostars, due to the expansion of H II regions. We used the millimeter dust continuum data of the ALMAGAL survey and the Anderson et al. 2014 catalog of H II regions and selected one ALMAGAL source for ALMA follow-up observations. In fact, in source AG286.0716$-$1.8229 six cores were detected at a resolution of $\sim7600$ au, but only two at a higher resolution. The 4 cores not detected at higher resolution are prestellar core candidates. We used archival data from the SMGPS and RACS to confirm whether an H II region is present in the field. We observed the source with with ALMA in Band 4, covering the emission of DCO$^+$ (2$-$1), N$_2$D$^+$ (2$-$1), DCN (2$-$1), and CH$_3$CCH (9$-$8), to estimates whether these cores are in an early phase of the star-formation process. The new Band 4 continuum image revealed three cores outside of the ALMAGAL field of view, for a total of 9 cores in the region, 8 of which are located along an arch of radius $\sim0.75$ pc. We have derived a spectral index between -0.14 and -0.4, in the frequency range of 0.8-1.6 GHz for the candidate H II region, which is consistent with optically thin free-free emission. Using plausible temperature ranges, based on the information from chemical tracers and the dust continuum, we derived mass ranges for the cores ($\sim2-16\,$M$_{\odot}$) and ranges for the virial parameter ($\sim0.3-5$). All the cores along the arch have virial parameters $\lesssim$2, with only one exception. Comparing the typical separation and mass of the cores with those expected in the case of the collect and collapse scenario and with the thermal Jean length and mass, the best agreement is found with the characteristic scales in the case of triggered star formation.

In this work we investigate the weak lensing convergence using an end-to-end nonlinear general relativistic framework. Combining numerical relativity simulations of large-scale structure formation with general relativistic ray-tracing, we compare our nonlinear calculation to the expectation based on perturbation theory for a set of 20 synthetic observers. We focus on large angular scales $\ell < 100$ across a broad range of redshifts with $0.05<z<3$. We confirm the importance of Doppler lensing for redshifts below $z\sim$0.6, as predicted by previous works. On average across our observers, linear perturbation theory predicts the nonlinear convergence to within 3-30% across all redshifts and angular scales we study. In general, we find smaller angular scales are better matched by linear theory than larger angular scales. While we cannot definitively identify the source of the discrepancy, for our particular study of redshift slices on observers' light cones the differences are mostly below the level of cosmic variance.

We consider the application of a ubercalibration-like relative flux calibration to the grism observations of the Roman High Latitude Wide Area Survey (HLWAS). We propose a simplified model of the calibration with an independent flat field for each detector in each exposure of the focal plane. In addition, we include two wavelength dependent components: a single wavelength throughput curve, modulated by a simple parabolic model for the throughput as a function of a source's focal plane position. We consider the impact of the dither scale, as well as the calibrator magnitude cuts. We show that the width of the calibration residuals can be reduced to less than 1.5 mmag, or 0.15% in flux, within the optimal dither range 50-240". This wide range allows for significant flexibility in optimising other parts of the observing program without diminishing the effectiveness of the relative flux calibration. We also discuss some improvements to the methodology that must be strongly considered before the calibration can be applied to real data. Finally, although we focused on spectroscopic component of the HLWAS here, our formalism and results should carry over to the imaging surveys as well.

The increasing detection of new molecules in the interstellar medium (ISM) shows that stereoisomerism is a fundamental contributor to interstellar molecular complexity. This work presents the first comprehensive overview of interstellar stereoisomerism. A total of 16 stereoisomeric pairs have been identified (13 conformational and 3 geometric), spanning molecules with 5-12 atoms and energy separations from 10 K to 2667 K. They were observed across diverse astrophysical environments with kinetic temperatures ranging from low to high values (7.5-300 K). The observed stereoisomeric ratios (OSR) - defined as the column density ratio of the higher-energy isomer divided by that of the lower-energy isomer - vary widely (0.009-4). While systems with small energy differences (1.2 kcal/mol) in hot environments (> 100 K) generally follow thermodynamic expectations (often assisted by tunneling-driven interconversion), many stereoisomers - particularly those in cold clouds or with larger energy separations - exhibit abundances far exceeding equilibrium values. This demonstrates that thermodynamics alone cannot explain interstellar stereoisomerism. Instead, stereoselective formation/destruction pathways (in the gas phase and/or in the surface of dust grains), photoisomerization, and chemical rearrangement during desorption must play a dominant role. Stereoisomeric ratios thus provide powerful constraints on interstellar chemical pathways, and about the physico/chemical conditions of the ISM. This review highlights the need for stereochemistry-sensitive astrochemical models. Progress in this field requires expanded laboratory spectroscopy of higher-energy stereoisomers, dedicated quantum chemical studies of isomerization processes, and the explicit inclusion of stereoselective chemistry in chemical networks.

While primordial black holes (PBHs) have long been a benchmark target for microlensing searches, the modern landscape of dark matter models suggests other, distinct, formation channels for compact objects made of dark matter. In the large class of self-interacting, dissipative models, dark matter has cooling channels that can enable fragmentation and gravitational collapse of some dark matter into compact objects, including black holes. The resulting populations have mass distributions, bias parameters, and abundance, spatial profile and velocity dispersion within the Milky Way that all differ from those of PBHs. We demonstrate that these population-level differences can leave imprints in the space of microlensing observables, with the differences in how the populations trace the dark matter giving the primary distinguishing lever. We discuss the possible overlap of microlensing signals from dark and baryonic lenses, and the complementarity of microlensing detection or constraints with other gravitational probes of novel populations of dark matter origin.

We reassess the direct collapse black hole (DCBH) interpretation of UHZ1 (UNCOVER-26185), a gravitationally lensed galaxy at $z_\mathrm{spec}=10.054$. That interpretation rests on a hard ($2-7$ keV) X-ray excess detected with Chandra, attributed to a Compton-thick AGN with an inferred $2-10$ keV luminosity of $L_\mathrm{X,int}\sim10^{46}~\mathrm{erg~s^{-1}}$ (Bogdan et al. 2024). The resulting extreme X-ray to rest-frame optical-IR ratio was taken as the hallmark signature of an "outsize black hole galaxy" at cosmic dawn. We analyse the full 2.2 Ms Chandra imaging dataset -- including 0.95 Ms of unpublished observations -- and present new JWST/MIRI photometry at $\lambda_\mathrm{obs}>5~\mu\mathrm{m}$. Across the full range of plausible Chandra data reductions, the $2-7$ keV excess at the position of UHZ1 reaches a significance of only $2.3-2.9\sigma$; the originally reported $4.2-4.4\sigma$ detection is sensitive to the specific astrometric alignment adopted and is not robustly reproducible. Moreover, the hard X-ray signal does not grow with the additional exposure, contrary to expectations for a steady source, indicating that any excess is not persistent. UHZ1 is also undetected in all nine MIRI imaging bands. Fitting red/obscured AGN SED templates to the tightest MIRI upper limit, we constrain the bolometric luminosity of any buried AGN to $L_\mathrm{bol}<1.3\times10^{45}~\mathrm{erg~s^{-1}}$. These conclusions are further supported by independent JWST spectroscopy (Alvarez-Marquez et al. 2026), which reveals no AGN signatures in the rest-frame UV or optical. Taken together, the multiwavelength data paint a consistent picture of UHZ1 as a low-mass, metal-poor, star-forming galaxy in the early Universe, with no compelling evidence for a luminous obscured AGN, regardless of its proposed formation channel.

We perform MHD simulations of AGN-blown bubbles in the Intercluster Medium (ICM) containing large-scale coherent magnetic fields. We assume that bubbles, created by the intermittent jets from Active Galactic Nuclei, quickly relax to the Woltjer-Taylor spheromak-like state, with internal plasma beta-parameter $\sim 1$. We demonstrate that such bubbles rising through hydrostatically-stratified atmosphere are magnetically stabilized against fluid interface instabilities, remaining coherent for a long time. Typical velocity is $ v /c_s \sim \sqrt{R/H} \leq 1 $ ($c_s$ is sound speed, $R$ is the bubble size, $H$ is the scale height). Current-driven instabilities (internal kinks) lead to bubble's tilting, but develop on long time scales, and remain unimportant, leading to minor modifications of the internal structure. Our results explain apparent long-term stability of ICM cavities. Subsonically rising stable bubbles dissipate in their wake approximately the energy initially injected by the jet, and may efficiently reheat the clusters cores in a ``gentle'' way.

The traditional definition of the circumstellar habitable zone (HZ) focuses on liquid water, but neglects the crucial role of ultraviolet (UV) radiation in prebiotic chemistry. Low-mass stars typically emit insufficient UV radiation for photochemistry throughout the liquid water HZs during quiescent states. However, frequent flares can provide substantial UV fluxes, potentially fostering habitable conditions. We refine the concept of a UV habitable zone (UV-HZ) by incorporating a temperature-dependent model for RNA precursor synthesis. Furthermore, we explore a parameterized spectral energy distribution model and adopt an empirical flare frequency distribution for flares on different stars to quantify their UV contribution. Applying this framework to different flaring stars, we find the UV-HZ around low-mass stars can extend to inner regions, and overlap with the traditional HZ in wide ranges. Apply the analysis to 9 planets around Kepler flaring stars, three planets are located within both the refined UV-HZ and liquid water habitable zone (LW-HZ) without causing ozone depletion. Our findings highlight the significant role of flares in expanding the potential for life around low-mass stars, offering a revised perspective on exoplanet habitability criteria.

This study investigates the X-ray characteristics of jet-associated radio-quiet AGNs across distinct optical/UV to X-ray correlation phases. Quasi-simultaneous optical/UV/X-ray observations of NGC 4051 from May-June 2009, obtained through Swift and XMM-Newton, reveal a temporal dichotomy: a strong optical/UV to X-ray correlation dominates the initial observation phase (before May 27), followed by an optical/UV flare event concurrent with X-ray flux suppression in the latter period. Our multi-method analysis of XMM-Newton data, incorporating short-term X-ray variability assessment, spectral decomposition, and RGS spectral analysis, identifies significant inter-phase X-ray emission disparities. During optical/UV flaring episodes, compared to the correlated phase, we observe: attenuated short-term X-ray variability amplitudes, enhanced soft X-ray absorption, suppressed intrinsic hard X-ray flux, and more prominent RGS emission-line features. Notably, these X-ray characteristics during optical/UV flaring intervals show no statistically significant deviations from pre-flare low-state X-ray emission patterns. These non-synchronous optical/UV-X-ray variations contradict predictions from both reprocessing models, starburst-driven emission scenarios, and the simplistic absorption models. While potential jet-related mechanisms remain ambiguous, our findings demonstrate strong consistency with predictions from the inhomogeneous accretion disk perturbation framework.

A bar-like central feature is commonly observed in both nearby and distant spiral-type galaxies, including the Milky Way. While many methods exist to categorise this morphology, no one method has emerged as the field-wide standard. To develop a rigorous and consistent method for identifying these bars, we investigate a classification scheme based on dynamical actions. In the Gaia era, actions can be estimated for individual stars in both observations and simulations, making this a natural and unifying diagnostic, assuming the intrinsic errors and selection functions are understood. Our approach is straightforward: stars that participate in the bar are subject to a strongly non-axisymmetric potential and, therefore, do not completely conserve their actions. We use this property to define a star-by-star criterion, formulated as an inequality and evaluated within measurement uncertainties, to identify bar members based on the degree to which their total action fails to be conserved. From tests on simulated galaxies, we find that the bar region is indeed characterised by a lower fraction of stars with conserved actions and that stars on bar orbits are represented by larger percentage changes in their actions. We are able to classify the spatial extent of barred region via the standard parameters of bar length and orientation, while also individually separating bar-located from bar-member stars on bar orbits. As proof of concept, our automated method based on dynamical actions robustly identifies bar parameters that closely match the eye's performance (average bar length variation ~9%) in barred snapshots of the test galaxy.

We present a search for neutrino counterparts to fast radio bursts (FRBs) using temporal and spatial cross-matching between the Second CHIME/FRB catalogue and the IceCube high-energy alert-track catalogue ICECAT1. Because current FRB--neutrino models do not provide a unique consensus on emission ordering, our primary significance test adopts a two-sided, order-agnostic temporal hypothesis. The analysis accounts for declination-dependent CHIME/FRB exposure and the look-elsewhere effect across multiple trials. No statistically significant FRB--neutrino association is found. The most significant pair is FRB\,20190630C--IC\,190629A, with a post-trial probability of $p=0.076$ ($1.43\sigma$), consistent with a chance coincidence. Within our statistical framework, a detectable physical association would require a time offset shorter than $\sim256$~s at $3\sigma$ or $\sim63$~ms at $5\sigma$. Using a population-level stacking analysis, we derive 90\% upper limits on the neutrino-to-radio luminosity ratio of FRBs, $\xi \lesssim 10^{8}-10^{11}$ for neutrino power-law spectral indices $\gamma=1.0-3.0$. These limits improve upon previous constraints by approximately two orders of magnitude and represent the most stringent bounds from FRB--neutrino coincidence searches to date. Although the current limits remain above the predictions of most magnetar-based models, they begin to constrain scenarios involving exceptionally efficient hadronic energy dissipation.

We present high-resolution optical spectroscopy of interstellar comet 3I/ATLAS (C/2025 N1) obtained with the High Dispersion Spectrograph mounted on the Subaru Telescope on UT 2026 January 7, when the comet was on its outbound trajectory at a heliocentric distance of $r_{\mathrm{h}} = 2.87$ au. The spectra cover the forbidden atomic oxygen lines, [O~I], at 557.7, 630.0, and 636.4 nm. The [O~I] red-doublet intensity ratio $I_{630.0}/I_{636.4} = 2.91 \pm 0.21$ matches the optically thin branching ratio ($\sim$3; \citealt{StoreyZeippen2000}), indicating that optical-depth effects are small and that our relative flux calibration is reliable. We measure a green-to-red [O~I] intensity ratio of $G/R = I_{557.7}/(I_{630.0} + I_{636.4}) = 0.339 \pm 0.027$. This value is higher than those of most Solar System comets at similar heliocentric distances, but comparable to that of the interstellar comet 2I/Borisov. From the measured $G/R$ ratio in 3I/ATLAS, we estimate the CO$_2$/H$_2$O abundance ratio under the assumption that H$_2$O and CO$_2$ are the dominant parents of O($^1$S) and O($^1$D), with other oxygen-bearing species expected to have a smaller influence under typical conditions (e.g., \citealt{FestouFeldman1981}). The derived ratio is significantly lower than the extremely CO$_2$-rich composition reported from infrared observations on the inbound trajectory at $r_{\mathrm{h}} \sim 3.3$ au, yet higher than typical values measured for comets in the Solar System. Together with published pre- and post-perihelion measurements, our result indicates that the CO$_2$/H$_2$O ratio decreased substantially across perihelion.

We present results from a mini-survey of dust spectral features arising in galaxies at redshifts $0.5 < z < 1.2$ in our James Webb Space Telescope (JWST) mid-infrared spectra of physically-unrelated background quasars. We analyze the JWST Mid-infrared Instrument (MIRI) Medium-Resolution Spectrometer (MRS) spectra of five quasars presented in Klimenko, Kulkarni, \& Aller 2025a (Paper I) to determine the best-fit silicate mineralogies. Template profile fits to the 10 $\mu$m feature suggest the possible presence of crystalline silicates in three of the galaxies. This contrasts with the predominately amorphous silicate grains in the Milky Way diffuse interstellar medium (ISM). We also measure the extinction curves using existing data from UV to mid-IR. Combining our results with past Spitzer IRS studies, we find that (i) the 10~$\mu$m silicate peak optical depth ($\tau_{10}$) is about three times stronger than expected for the local diffuse ISM over the range $A_V =0.1-2.0$, with $\tau_{10}/A_V$=$0.17\pm0.09$. (ii) The relative strength of the UV bump is similar to that in the local ISM. However, the ratio $\tau_{10}/A_{2175}$ is larger ($\sim0.1-1$), and appears to decrease with $A_V$, approaching the Galactic ISM value ($\sim 0.1$) at $A_V\sim1.5-2$. (iii) No significant correlation of $\tau_{10}/A_V$ with $R_V$. (iv) $\tau_{10}$ is strongly correlated with the gas-phase Mg~II absorption strength for the quasar sightlines. Possible interpretations include that some quasar sightlines probe dust in the circumgalactic medium (CGM), and that dust grains may have been significantly reprocessed in the ISM and CGM under conditions that may differ from those in the local ISM.

AstroSat/LAXPC and SXT observed the persistent neutron star low-mass X-ray binary 4U 1820-30 between 2016 and 2022. During these observations, the hardness-intensity diagram (HID) and color-color diagram (CCD) indicated that the source was in the banana state. We divided the CCD into 11 segments for spectral and timing analyses. For each segment in the CCD, we modeled the spectral data using two distinct approaches over the 0.7-20.0 keV band. A combination of a multi-color-disk component with an inner disk temperature of around 0.6 keV and Comptonized emission from the boundary layer (BL)/ hot corona provided the best description of the X-ray spectral data of this source. The truncation radius was found to be in the range of $\sim$ 19-40 km. The Comptonized component has an optical depth in the range of $\sim 7 - 13$ with electron temperature in the range of $\sim 2.5 - 3.8$ keV. The optical depth of the corona varies significantly along the position on the CCD, while $\sim$ 80\% of the X-ray flux comes from the Comptonized component. We discuss possible physical scenarios to explain the relationship between the spectral evolution and motion of the source along the CCD. The timing analysis revealed kHz QPOs peaks at $\sim 710$ Hz and $\sim 740$ Hz in the lower left banana branch. An energy-dependent study indicates that these QPOs are stronger in the high-energy band.

Cosmological filaments, galaxy clusters, and galaxies are magnetized reservoirs of cosmic rays (CRs). The exchange of CRs across these structures is usually modeled assuming that they remain charged and magnetically confined. At high energies, hadronic interactions can convert CR protons to neutrons. This physics is routinely included in air-shower and ultra-high-energy (UHE) CR propagation Monte Carlo simulations used for composition studies but is rarely treated explicitly in propagation models of CR transport and exchange between magnetized reservoirs. CR neutrons are not affected by magnetic fields and can propagate ballistically over kpc-Mpc distances before decaying back into protons, with relativistic time dilation extending their effective decay length. We show how such charged-neutral switching modifies CR confinement and escape in four representative environments: a Milky Way-like galaxy, a starburst galaxy, a galaxy cluster, and a cosmological filament. By solving the transport of a confined CR proton population in each structure using a diffusion/streaming propagation approach with hadronic pp and p$\gamma$ interactions, and treating neutron production and decay as a stochastic Poisson ''jump'' process, we find that neutron-mediated steps can allow additional CR escape from large-scale cosmological structures at energies where charged-particle transport alone would predict strong CR confinement and attenuation in ambient radiation fields. These effects imply a qualitative shift in how ultra-high-energy CRs are transferred from embedded sources into filaments and voids once intermediate neutron propagation is considered, with consequences for the partitioning of CRs across the large-scale structure of the Universe.

Tidal disruption events (TDEs) are powerful probes of supermassive black hole (SMBH) properties and accretion physics. The existing light curve fitting tools assume that the disrupted stars are on parabolic orbits, which may introduce systematic biases in derived parameters. In this work, we develop a non-parabolic TDE model that incorporates orbital energy of the disrupted star as a free parameter ($\tilde{\epsilon}_{\rm orb}$) to modify the debris mass distribution and mass fallback rate. We apply this model to 30 TDEs from the ZTF-I survey and compare the results with those from a standard parabolic model. We find that neglecting orbital energy leads to biased black hole mass estimates: for eccentric (hyperbolic) orbits, parabolic models systematically underestimate (overestimate) the black hole mass. Additionally, we measure orbital eccentricities ($e$) and penetration factors ($\beta$) of the disrupted stars in this sample, enabling an investigation of their origins via the $e$-$\beta$ parameter space. Most events (24/30) are consistent with production via two-body relaxation in spherical nuclear star clusters, but six outliers with high $\beta$ and $e<1$ suggest alternative mechanisms. Our results highlight the importance of accounting for orbital energy in TDE modeling to improve the accuracy of SMBH mass measurements and to better understand the dynamical origin of the disrupted stars.

The scientific potential of X-ray polarimetry has long been recognized, but the challenges in measuring polarization have left it largely unexplored, particularly in the hard X-ray regime. While tremendous advancement has been made in soft X-ray polarimetery, the lack of sensitive hard X-ray polarimeters and polarisation measurements continues to limit our understanding of high-energy astrophysical processes. With the development of hard X-ray mirrors, it is now possible to develop a sensitive focal plane hard X-ray polarimeter. One such effort is CXPOL, a prototype developed at PRL, India, which consists of a plastic scintillator as active scatterer readout by PMT surrounded by CsI(Tl) scintillators in cylindrical array with SiPM readout from one side. First results of the prototype have been demonstrated in 20 to 80 keV energy range. The sensitivity of the instrument can be significantly enhanced using faster and better light yield scintillator like NaI as absorbers. Further, the use of a position-sensitive scatterer and absorbers, can also provide spectroscopic information by measuring the interaction position along the length and from the known energy depositions in the detectors. Position sensitive detectors are also helpful in mitigating the systematic effects introduced by the off-axis events in the polarisation measurements. Here, we demonstrate the detection sensitivity in the 100x20x5 mm^3 NaI(Tl) scintillator absorber readout on both ends by SiPM arrays operating in co-incidence. In this work, we characterize the first prototype of this detector system and investigate the variation in energy and position resolution, and light output with irradiation position along the length of the detector. The two end readout in co-incidence also reduces the overall SiPM background per absorber by an order of magnitude, further enhancing the polarimetric sensitivity of the instrument.

While there are a number of proposed formation channels for subsolar mass compact objects, including black holes formed primordially, or neutron stars that form in collapsar disks, there have yet to be any conclusive observations of such objects. Motivated by the possibility that, if such objects exist, gravitational waves from binary mergers may reveal them, we study binary neutron star mergers where one star has a subsolar-mass in order to determine how well such systems are described by current models, and when they could be distinguished from a system with a subsolar-mass black hole. We perform fully general-relativistic simulations of a $1.7\ M_{\odot}$ star merging with a $0.8\ M_{\odot}$ star, leading to tidal deformabilities of up to $\mathcal{O}(10^4)$ for the latter, and quantify how this affects the merger dynamics and associated gravitation and electromagnetic signals. In this regime, we find mass transfer between the stars, as well as significantly lower disruption frequencies. Though this is not captured by current gravitational waveform models, we conclude that this does not significantly impact the sensitivity of current gravitational wave detectors to these sources. Assuming design sensitivity of the LIGO and Virgo detectors, we find no biases in the recovered intrinsic parameters for signal-to-noise ratios $\lesssim 100$. We also find that the large deformabilities lead to a significant increase in the amount of dynamically ejected matter compared to equal mass systems, exceeding the predictions of current phenomenological models.

The discovery of the complex stellar populations hosted in two massive stellar systems in the Galactic bulge, namely Terzan5 and Liller 1, posed intriguing questions about their origin. Despite their globular cluster appearance, they host sub-populations with significantly different ages (several Gyrs) and metallicities (about 1 dex) tracing a chemical abundance pattern that is consistent only with that observed in the bulge. These surprising properties can be naturally explained in the context of a self-enrichment scenario, opening the possibility that they could be the remnants of primordial massive structures that contributed to the bulge formation (the so-called Bulge Fossil Fragments, BFFs) capable of retaining supernova ejecta within their potential well. In this paper we present a first attempt to quantify the expected contribution of BFFs to the gravitational wave emission. In particular, by adopting Terzan5 as prototype of BFF, using its chemical evolutionary model, and following a scaling relation derived for globular clusters, we present a first-guess estimate of the number of binary black hole (BH) mergers expected in this stellar system. Within the adopted simplifying assumptions and the uncertainties about the initial conditions of the proto-Terzan 5 system, we find that several hundreds of binary BH mergers are expected, a number that is between 15 and 250 times larger than that produced by a typical globular cluster. Hence, this study identifies in the BFF family a new population of stellar systems potentially able to produce a significant number of gravitational wave emitters, that has not been considered in any previous investigation. Moreover we speculate that they could also be the natural place where BHs with masses above 60 Msun and even intermediate-mass BHs can form via repeated dynamical interactions.

The East Asia VLBI Network (EAVN) has recently enabled dual-polarization observations at $22$ and $43\,\mathrm{GHz}$. We present the first systematic verification of its polarimetric performance using EAVN observations of M87, 3C 279, 3C 273, and OJ 287, calibrated with the GPCAL pipeline and evaluated against near-contemporaneous VLBA images at comparable frequencies. Most stations show stable polarimetric leakages with amplitudes of $5$-$10\%$ over monthly timescales. While several VERA stations exhibit D-term phase variations between epochs, we attribute these to field-rotator (FR) offsets and demonstrate that phase stability is restored after applying the analytically derived FR corrections. The resulting linear-polarization morphologies and EVPAs broadly agree with the VLBA results within uncertainties; fractional polarization measured by the EAVN tends to be slightly higher near polarization peaks. Although exact one-to-one comparisons are limited by moderate frequency and epoch differences, the combined evidence indicates robust EAVN polarimetric calibration and imaging capabilities at $22$ and $43\,\mathrm{GHz}$. These results support the scientific capability of EAVN polarimetry and lay the groundwork for expanded, higher-fidelity polarimetric studies in East Asia.

Earth's long-lived geodynamo is difficult to reconcile with recent high estimates of the core thermal conductivity, a problem known as the new core paradox. At the same time, the long-term thermal evolution of the mantle remains uncertain, largely due to the poorly constrained onset of modern-style plate tectonics, which marks the transition to efficient cooling of the interior through mobile-lid convection. Because core cooling -- and thus magnetic field generation -- depends on the efficiency with which the mantle extracts heat from the core, these two problems are closely linked. Here, we investigate the coupled thermal evolution of mantle and core using a 1D model that incorporates a parametrized transition transition from stagnant- to mobile-lid convection, defined by its onset time and with a prescribed duration. This framework allows us to assess how different tectonic histories influence Earth's thermal and magnetic evolution. We perform a Bayesian inversion using constraints from the palaeomagnetic record, mantle cooling history, and present-day thermal state. Our results favour a transition from stagnant- to mobile-lid convection during the Archean, which promotes core cooling and enables a geodynamo throughout Earth's history, even for core thermal conductivities in excess of 100 W/m/K. A delayed onset of mobile-lid convection provides thus a viable solution to the new core paradox.

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will produce unprecedentedly deep and wide photometric catalogs, enabling transformative studies of faint stellar systems such as the research of ultra-faint dwarf galaxies (UFDs). A critical challenge for these studies is reliable star-galaxy separation at faint magnitudes, where compact background galaxies increasingly contaminate stellar samples. This work aims to assess the performance of supervised machine-learning techniques for star-galaxy separation in LSST-like data, quantify the relative importance of morphological and photometric information, and identify the most effective combinations of input features for minimizing galaxy contamination while preserving stellar completeness in the faint regime relevant for UFD searches. We apply a Random Forest classifier to observations of the Extended Chandra Deep Field South from LSST Data Preview 1 (DP1), the deepest field observed within the DP1. We construct a curated sample of bona fide stars and galaxies using spectroscopic data, Gaia DR3, and multi-band photometric catalogs. We train and validate the classifier using several configurations of LSST-based input features, including multi-band colors, the LSST morphological parameter refExtendedness, and photometric uncertainties. We find that LSST multi-band photometry alone delivers a good star-galaxy separation, significantly outperforming morphology-based classification at faint magnitudes. Colors involving the u-band are essential to provide a robust star galaxy separation. Furthermore, explicitly including photometric uncertainties as input features yields the best overall performance. Across all configurations that include all the six LSST filters, galaxy contamination remains negligible almost the whole magnitude range probed in this work (i.e. r < 27.5 mag). [abridged]

ProtoPol is a medium-resolution echelle spectro-polarimeter developed for Physical Research Laboratory (PRL) 1.2m and 2.5m telescopes, Mt. Abu, India. Though initially conceived to evaluate the development methodology of the echelle spectro-polarimeter, it was subsequently elevated to the level of a full-fledged back-end instrument for PRL telescopes. ProtoPol is developed on the traditional concept of using a half-wave plate with Wollaston prism to achieve the separation of two mutually orthogonal polarized beams. These separated beams are modulated and directed into an echelle spectrometer which is employs an echelle grating and two plane reflection gratings as the cross-dispersers. Therefore, the cross-dispersed spectra for two orthogonal polarized beams are recorded in multiple orders on a CCD detector. ProtoPol is designed to operate in the visible and near IR spectral range, 4000 - 9600 angstrom, with a spectral resolution ($\delta$$\lambda$) around 0.4-0.75 angstrom. The uniqueness of ProtoPol lies in its design which has entirely been developed with commercially available off-the-shelf optical and opto-mechanical components. This feature makes ProtoPol a noteworthy development as it offers a cost-effective way to develop spectro-polarimeters with such resolutions for small-aperture (2-3m) telescopes around the world, in a much shorter development period. ProtoPol has been successfully developed and commissioned on PRL 1.2m and 2.5m telescopes since December 2023, and a variety of observations have been carried out for instrument characterization, performance verification, and scientific purposes. This is the first of the two-part research articles series, wherein we present the design and development methodology of ProtoPol, along with its laboratory characterization and performance.

Between 2.4 and 0.6 Gy ago, our planet underwent several episodes of global glaciations, including the Snowball Earth case that ended 635 My ago. Causes of this last Snowball event presumably included a decreased greenhouse gas concentration and high continental albedo, both associated with the passage of the super-continent Rodinia at equatorial latitudes. When large continental masses are in equatorial regions, silicate weathering is enhanced, leading to decreased atmospheric CO2 concentration, while the bare continental masses, which at the time hosted no vegetation, enhanced reflection of solar radiation. Since then, no other Snowball episodes were recorded. Here we numerically explore the climatic dynamics of a rocky planet for different values of solar output, continental configuration (current and Rodinia-like), CO2 concentration and continental albedo, simulating the effects of land vegetation. We found that for the solar input typical of 600-700 My ago (95% of the current value), the presence of bare continents with albedo 0.35 (granite) in the position estimated for Rodinia was sufficient to trigger a Snowball state for CO2 concentrations up to at least 1000 ppm. When bare continents are located in modern positions, Snowball could be triggered only for values of CO2 concentration below 400 ppm. At current solar input values, Snowball states appear only at or below 100 ppm. Thus, we found: a lower solar output is an essential component of the transition to Snowball; the presence of land vegetation is crucial and reduces the probability of entering a Snowball state; a low CO2 concentration was not needed for triggering a Snowball in bare Rodinia-like conditions and reduced solar output; current solar luminosity does not allow Snowball states, even for equatorial continents, unless continental albedo is that of granite and CO2 concentration is 100 ppm or less. [Abridged]

The Taiji mission, a pioneering Chinese space-borne gravitational wave observatory, requires ultra-precise calibration of its gravitational reference sensors (GRSs) to achieve its targeted sensitivity of $3\times10^{-15} \mathrm{\ m\ s^{-2}\ Hz^{-1/2}}$. Maintaining this precision is challenged by time-varying scale factors drifts and dynamic center-of-mass (c.m.) offsets between the test masses (TMs) and spacecraft, driven by factors such as propellant consumption, thermal effects and aging electronics. This paper develops an advanced in-orbit calibration framework that simultaneously estimates the GRS scale factors and c.m. offsets between TMs and spacecraft through a combination of spacecraft maneuvers and Kalman filter. By applying periodic torque signals to induce controlled spacecraft angular accelerations, we leverage star tracker and GRS readouts to disentangle coupled disturbances and achieve dual-parameter calibration with unprecedented precision, with scale factors errors below 0.2\% and c.m. offsets residuals within 100 $\mathrm{\mu}$m, satisfies the Taiji-2 calibration requirements. This method is robust across different satellite configurations. The results not only ensure the feasibility of Taiji-2's scientific objectives but also establish a scalable calibration paradigm for future missions such as Taiji-3, where sub-micrometer c.m. stability and ultra-low noise gravitational reference will be essential.

Flat-Spectrum Radio Quasars (FSRQs) are among the most energetic and powerful active galactic nuclei, often exhibiting jet powers comparable to or exceeding the Eddington luminosity. In this work, we performed broadband spectral energy distribution (SED) modeling of two FSRQs PKS 1441+25 and Ton 599, using Swift-XRT/UVOT, NuSTAR, Fermi-LAT and VERITAS observations during 2015 and 2021, respectively. We considered four particle distribution models: a broken power law, a log-parabola, and two energy-dependent models in which either the diffusion or acceleration timescale depends on energy. Our results show that the jet power estimates derived from models with intrinsic curvature, such as the log-parabola and energy-dependent models, are of the same order as those obtained with a broken power-law distribution. This contrasts with the case of High Synchrotron Peaked Blazars (HBLs), where the power estimates can differ by nearly two orders of magnitude between models. We attribute this difference to the lower electron break energies typically observed in FSRQs. Consequently, our findings suggest that, unlike in HBLs, the estimated jet powers in FSRQs are relatively insensitive to the assumed particle energy distribution, reflecting the dominance of external Compton processes and weaker dependence on spectral curvature.

We investigate the impact of observational and theoretical uncertainties in cluster gravitational redshifts as a probe of modified gravity using an end-to-end forecasting pipeline. We use a generative model to build a halo catalogue with $M_{500}\ge 3\times 10^{13}\,M_\odot$, populate haloes with member galaxies via a five-parameter halo occupation distribution (HOD), assign projected positions from radial density profiles, apply survey-like selections, and infer a linear rescaling of the gravitational potential, $\alpha_\mathrm{MG}$, to parameterise modifications to general relativity (GR). We vary redshift uncertainties, radial and mass-redshift completeness, member abundance, minimum mass and maximum redshift, as well as mis-specify the clusters density and velocity profiles, centres, and mass function. We find that the intracluster velocity dispersion sets an effective floor: improving redshift precision beyond $\sigma_z\sim 10^{-4}(1+z)$ brings no improvement in the precision of $\alpha_\mathrm{MG}$. Realistic redshift and mass cuts primarily remove low-mass haloes and have minimal impact on the $\alpha_\mathrm{MG}$ precision. In this setting, we find that shallow, narrower spectroscopic surveys are preferable to deep, wide photometric ones for precise modified gravity constraints. We further find that mis-centring can mimic significant departures from GR. Baryonic deviations from a Navarro-Frenk-White profile and velocity anisotropies do not introduce appreciable biases. In the high-S/N regime of upcoming surveys, accurate determination of cluster centres will be essential to avoid interpreting systematic effects as new physics. The Spectroscopic Stage-5 Experiment and the Widefield Spectroscopic Telescope provide a clear route toward establishing gravitational redshifts as a competitive probe of modified gravity.

We present the development of ProtoPol - a medium resolution echelle spectro-polarimeter for the PRL 1.2m and 2.5m telescopes at Mt Abu observatory, India. In this second and final part of the paper series, we report on the development of a dedicated data reduction pipeline of ProtoPol along with several characterization, performance evaluation, and scientific observations to quantify the performance of the instrument. ProtoPol provides a spectral resolution in the range of $\sim$0.4 - 0.75 angstrom across various orders in the visible wavelength range of 4000-9600 angstrom. On PRL 2.5m telescope, an SNR of 10 is achieved for $m_V\sim13.2$ source in 1 hour of integration time, and its throughput is estimated to be $\sim$6\% including all the contributing factors such as atmospheric transmission, telescope reflectivity, instrument's optics, CCD efficiency etc. ProtoPol achieved a linear polarization accuracy $\delta P \approx 0.1-0.2\%$ in 2 hours of integration time for a source with $m_V\approx8$. The instrumental polarization is determined to be around $0.1\%$. We also present the first science results with ProtoPol to demonstrate the capabilities of the instrument. A sample of Herbig Ae/Be stars, classical Herbig stars, Symbiotic stars, and AGB/post-AGB stars were observed over the period of one and half years for their spectro-polarimetry measurements covering various physical mechanisms such as intrinsic line polarization in Herbig and classical Be stars, Raman scattered features in Symbiotic stars, as well as continuum polarization in AGB/post-AGB stars to verify the polarization performance of the instrument.

We present a detailed analysis of the neutral atomic gas (HI) in the dwarf galaxy UGCA 320, observed with the MeerKAT telescope as part of the MHONGOOSE (MeerKAT HI Observations of Nearby Galactic Objects: Observing Southern Emitters) programme. In a small group consisting of three dwarf galaxies, all of which contain HI, it is the most massive. Detailed kinematic modelling shows that UGCA 320 contains a substantial amount of (kinematically) anomalous gas (>=20%), at least ~30% of which is likely the result of a tidal interaction with its neighbour UGCA 319. It also reveals that UGCA 320 likely harbours a star-formation driven outflow, and that ~10% of its HI is extra-planar and has a filamentary structure. Although UGCA 320 aligns with established scaling relations from the literature, its neutral hydrogen content is notably complex - shaped by its immediate environment. This underscores the importance of deep, resolved observations and detailed kinematic analyses to capture the nuances of galaxy evolution.

We focus on the automated classification of eclipsing binary stars using deep learning methods to handle the vast data generated by large-scale photometric sky surveys. These surveys produce extensive datasets that are impractical for manual analysis. By using machine learning to classify eclipsing binary stars based on light curve morphology, this study aims to contribute to the efforts to efficiently process and accurately interpret massive data from projects Kepler, TESS and Gaia missions.

Angular broadening observations of compact radio sources provide a powerful method for probing the solar corona and solar wind. Such observations enable studies of the phase structure function, turbulence amplitude, intermediate-scale density fluctuations, solar-wind heating rates, and dissipation scales. When a compact radio source is observed through coronal or solar-wind plasma, several observable effects can arise: (1) the apparent source size increases because of scattering by turbulent plasma, producing angular broadening; (2) the source flux density decreases because of scattering and absorption; (3) the observed angular broadening can be anisotropic, reflecting anisotropic turbulence in the corona and solar wind; and (4) the position angle of the anisotropy, measured from north through east, can help infer the orientation of the coronal magnetic field. These effects provide insights into the physical processes governing the solar wind and its interaction with electromagnetic waves, and they offer constraints on coronal turbulence and magnetic-field structure. At present, angular broadening studies remain limited and have mostly focused on very bright radio sources such as Tau A. The unprecedented sensitivity and angular resolution of the Square Kilometre Array are expected to greatly expand the number of suitable background sources, opening a new window on the solar corona, solar wind, and heliosphere.

Using Parker Solar Probe measurements of the solar wind, we demonstrate that $\beta_{\parallel}$ is the main driver that determines which instabilities limit proton temperature anisotropy. At radial distances from 10 to 30 solar radii, $\beta_{\parallel}<1$ drives electromagnetic ion-cyclotron and parallel firehose instabilities, in contrast to the situation at 1 astronomical unit, where, due to most $\beta_{\parallel}>1$, mirror and oblique firehose modes are dominant instead. Furthermore, we show that the temperature anisotropy radially evolves following the semi-empirical anti-correlation $T_\perp/T_\parallel\sim\beta_\parallel^{-0.55}$, consistent with observations at larger distances from the Sun.

Strong solar activity is often accompanied by a variety of radio bursts. These bursts are valuable diagnostics of coronal and heliospheric processes and also have potential applications in space weather monitoring and forecasting. However, space weather applications require low-latency, high-sensitivity radio burst recording and reporting capabilities, which have remained limited. In this work, we present the development of a near-realtime radio burst recording and reporting system using the Owens Valley Radio Observatory Long Wavelength Array. The system directly clips data from a realtime buffer and streams them as a live radio dynamic spectrogram. These spectrograms are then processed by a deep-learning-based burst identification module for type III radio bursts. The identifier is based on a YOLO (You Only Look Once) architecture and is trained on synthetic type III radio bursts generated using a physics-based model to achieve accurate and robust detection. This system enables continuous realtime radio spectrum streaming and automatic reporting of type III radio bursts within approximately 10 seconds of their occurrence.

We present HEALFormer, a transformer-based neural network architecture for weak gravitational lensing mass mapping that reconstructs convergence maps from incomplete and noisy shear observations on the celestial sphere. The model operates directly on the Hierarchical Equal Area isoLatitude Pixelization and employs learnable mask tokens to handle arbitrary survey geometries without requiring preprocessing. Through a progressive training strategy, HEALFormer efficiently processes high-resolution maps up to Nside = 1024 and demonstrates excellent performance across diverse survey footprints including KiDS, DES, DECaLS, and Planck. The model generalizes robustly to cosmological parameters beyond its training set, producing nearly unbiased reconstructions with superior noise suppression compared to traditional Kaiser-Squires and Wiener filter methods. Remarkably, HEALFormer exceeds the theoretical phase recovery limits of linear reconstruction methods at small scales, achieving a fundamental breakthrough in weak lensing analysis. The combination of computational efficiency, reconstruction accuracy, and adaptability to varying survey configurations makes HEALFormer well-suited for current and next-generation cosmological surveys. Code is available at GitHub.

ALMA observations have shown that substructures are ubiquitous in protoplanetary discs. A sub-group, the transition discs, shows large cavities and rings in dust continuum. Among these, some present very high contrast asymmetries possibly due to the presence of vortices. HD 34700A is a binary system featuring a cavity, a ring, and multiple spiral arms detected in scattered light, a prominent crescent in the ALMA continuum and a complex gas morphology possibly connected with ongoing infall. We present new ALMA band 6 (1.3 mm) continuum images of the circumbinary disc around HD 34700A and compare them with two other systems showcasing high ($\gtrsim30$, measured as the peak-to-azimuthal-average ratio) contrast continuum asymmetries, IRS 48 and HD 142527. We aim to characterise the crescent morphology and discuss their possible origin. We perform visibility modelling of the new high resolution (0.''11x0.''09) ALMA band 6 continuum data of HD 34700A, together with improved visibility modelling of the other two targets. Our visibility model is in remarkable agreement with the HD 34700A data, featuring only localised residuals in the region of the disc corresponding to the tail of the asymmetry. We reproduce the double-peaked emission in HD 142527, and recover the crescent shape in IRS 48. We then run a hydrodynamic model of a vortex with different dust fluids, reproducing the general asymmetric crescent morphology observed in the HD 34700A and IRS 48 systems. With a combination of visibility, dust evolution and hydrodynamical models, we have constrained the morphology of the dust continuum emission of HD 34700A for the first time, and improved existing models for IRS 48 and HD 142527. The high azimuthal contrast of the asymmetries rules out the orbit clustering of eccentric cavities scenario, while the dust evolution models we consider suggest that the vortex scenario is a plausible option.

Missing channels in radio-interferometric visibility data can introduce systematic artifacts into the estimated 21-cm power spectrum. A common workaround is to first estimate the two-frequency correlation $C(\Delta\nu)$ and then Fourier-transform it to obtain the power spectrum $P(k_\parallel)$. This procedure yields an unbiased estimate when the signal is statistically homogeneous (ergodic) along the line-of-sight, but it fails in the presence of non-ergodic foregrounds. Smooth Component Filtering (SCF) has recently been proposed as a solution to this problem, in which the dominant non-ergodic (spectrally smooth) component is removed prior to estimating $C(\Delta\nu)$. In existing implementations, the smooth component is estimated by convolving the visibilities with a Hann window along the frequency axis. We demonstrate that this Hann-based SCF performs adequately only when foregrounds are extremely spectrally smooth, i.e., when they possess a long frequency-correlation length. In contrast, it breaks down when foregrounds exhibit short correlation lengths, as is frequently encountered in real observations. We introduce a Bayesian extension, Bayes-SCF, based on Gaussian Process (GP) regression, which overcomes this limitation. Bayes-SCF models the smooth component via a covariance function with a fixed correlation length, enabling a controlled and data-driven filtering. Using simulated data, we show that Bayes-SCF robustly recovers the input model 21-cm power spectrum even in the presence of spectrally unsmooth foregrounds. Bayes-SCF is also effective in a delay-spectrum approach. The primary trade-off introduced by the Bayesian framework is increased computational cost; future work will focus on optimizing the algorithm and applying it to real MWA data.

Fuzzy Dark Matter (FDM), composed of ultra-light axions ($m_a \sim 1 \times 10^{-22}$ eV), exhibits wave-like properties that can significantly impact early-universe star formation. Using the arepo code with the axirepo module, we simulate the assembly of haloes across a range of axion masses ($1 \times 10^{-22}$ eV $\le m_a \le 7 \times 10^{-22} $ eV) and halo masses ($3 \times 10^{8} M_\odot \le M_h \le 8 \times 10^{9} M_\odot$). We investigate how small-scale dynamics of the FDM density field affect the accumulation of cold, dense gas. We find that the delay in star formation scales inversely with both halo mass and axion mass. While the static, cored geometry of the soliton primarily sets the timing of the initial collapse, we identify a secondary dynamical barrier driven by stochastic fluctuations that is most potent at the low-mass end of our parameter space. These dynamics dictate the spatial scale of dense gas by injecting kinetic energy and inducing significant angular momentum, which can rotationally stabilize gas out to the soliton radius. This wave-driven stirring leads to the formation of extended $\text{H}_2$ plumes and promotes dynamical mixing, effectively starving the central regions and forcing gas to cool in a more fragmented, diffuse manner. Our results indicate a shift from the monolithic central star formation seen in CDM toward lower-mass, fragmented clusters. These internal inefficiencies provide a physical mechanism for delaying Cosmic Dawn beyond the effects of the power spectrum cut-off, which is essential for refining observational constraints on the axion mass.

Particle physics and gravitational waves provide complementary probes of the deep structure of the Universe. Gravitational waves from the mergers of neutron stars and black holes are sensitive to the structure of dense quark matter and to different dark matter scenarios. Measurements of stochastic gravitational waves backgrounds can teach us about possible first order phase transitions in the early Universe, including providing sensitivity to the TeV scale which is of key interest to future particle collider experiments. Gravitational waves measurements will also give new probes of the evolution and expansion of the Universe, complementary to measurements with electromagnetic radiation. This Perspectives article explores the physics synergies between the science opportunities provided by next generation gravitational waves measurements and particle physics experiments. Gravitational waves can also probe deep into the early Universe reaching physics much above possible collider energies if the signals can be detected.

The properties and spatial distribution of large-scale structures of the solar corona determine the observed solar wind structure at 1 au. Coronal holes are a major source of fast solar wind, an important geo-effective component, and appear as large dark patches in extreme ultraviolet images. Solar observatories provide images of the solar corona at different wavelengths, enabling identification of coronal hole morphology and other large-scale structures along a given line of sight. The problem is that although models exist, few work in real time, separate coronal holes from other dark features, or are fully automatic and suitable for comparison with models. The main goal of this work is to develop an automatic threshold-based coronal hole detection tool across solar cycles 23, 24, and 25 using artificial intelligence. The only user input is the date, enabling retrieval of the threshold value used to detect coronal hole contours in line-of-sight extreme ultraviolet images from SDO/AIA and SoHO/EIT. We retrieve data affecting the threshold due to contrast changes from the Heliophysics Events Knowledge database for large-scale features such as active regions, solar flares, coronal mass ejections, and filaments, and engineer them to train the neural network model (POP-CORN). The model input comprises categorical features of large-scale structures in the solar corona, including spatial distribution and properties such as solar flare class by intensity. The neural network model is trained to achieve higher accuracy and determines the threshold needed to detect coronal holes, allowing their boundaries to be identified automatically and consistently. We conclude that properties of large-scale structures affect the determination of coronal hole regions, and incorporating these properties into training improves detection.

We investigate a class of interacting dark energy (IDE) models arising from density-driven spontaneous symmetry breaking in a conformally coupled scalar-tensor framework. In this construction, the dark matter-scalar interaction is dynamically activated as the cosmological density evolves, and the redshift dependence of the coupling follows a logistic profile whose steepness is determined by the local curvature of the symmetry-breaking potential. Working in the controlled adiabatic tracking regime, we implement the resulting epoch-dependent interaction in a perturbative background close to $\Lambda$CDM and confront the model with late-time cosmological data, including Planck 2018 CMB lensing reconstruction, redshift-space distortions, and Pantheon+SH0ES supernova data. We analyze realizations in which the activation index is allowed to vary and compare them with a restricted realization in which it is fixed to the canonical quadratic minimum value, thereby probing the structural role of the activation profile. We find no statistically significant preference for interaction over $\Lambda$CDM; current observations constrain the model to a hierarchical regime in which the scalar remains heavier than the Hubble scale at activation and background deformations remain perturbatively small. Allowing the activation index to vary preserves an extended degeneracy direction in parameter space, whereas fixing it removes this freedom and leads to a contraction of the allowed posterior region once geometric and growth data are combined. Our results delineate the viable parameter regime of symmetry-breaking IDE and clarify the structural distinction between microphysically motivated scalar-tensor realizations and phenomenological interacting models.

The circumgalactic medium (CGM) regulates galaxy growth and retains the imprint of feedback from supernovae and supermassive black holes. However, the bulk of the hot CGM produces little X-ray emission and is challenging to study with X-ray telescopes. We propose a novel method for evaluating the CGM mass using resonant scattering of the helium-like oxygen (OVII) resonant line at $E=574$ eV. In a spherically symmetric and static CGM halo with a sharp central X-ray peak, the number of OVII ions within an outer radial shell can be calculated from the ratio of the two directly observable quantities: the OVII flux from the bright inner region and the scattered OVII flux from the shell (where the scattered flux can be much higher than the intrinsic emission). To evaluate the accuracy of this geometric estimate for realistic galaxies -- with satellites, asymmetries, and gas velocities -- we use a sample of galaxies from the TNG50 cosmological simulation. We find that, when the most irregular systems are excluded based on their X-ray observables, we accurately predict the OVII mass in the outer halo (e.g., in an $r=R_{\rm 500c}-R_{\rm 200c}$ shell) from the ratio of the fluxes in the corresponding annulus and the central peak region ($r<0.2R_{\rm 500c}$), with only a 10% bias and an rms scatter of $\sim 0.2$ dex. As OVII mass strongly correlates with the total oxygen and gas mass, this direct OVII-counting method enables indirect estimates of those quantities by future X-ray microcalorimeter missions, such as NewAthena and HUBS.

As the (JWST) pursues its observing journey, several thousands of icy-grain spectra are expected to be measured and analysed. The inventory of ices in particular, via the observations of background sources, is accessible for hundreds of lines of sight (LOSs) per molecular-cloud region, opening the possibility to add strong constraints on the solid phase chemistry in a vast domain of cloud densities. SynthIceSpec is a synthetic infrared (IR) spectrum generator that has been designed as a tool to support observing proposals and to test the outcome of chemical models. It is based on laboratory measurements of pure and mixed ices, where each vibrational component is fitted by a sum of Gaussian profiles. Given an initial ice chemical composition (either set by the user or the outputs of a chemical model), a full JWST spectrum is generated, to which the contribution of silicates; continuum, stellar photospheric absorption bands; and extinction law can be added. For the continuum, stellar photospheric models for a wide range of spectral types can be selected by the program, or, Spectral Energy Distribution (SEDs). We present a few use cases of SynthIceSpec: we probed the impact of dust temperature on CO_2 ice formation using IR data and gas-grain modelling. Next, we used SynthIceSpec to explore the detectability of the main feature of CH_3CN at 4.45 um in a cold core environment with the JWST, which was previously tentatively detected in YSOs. The detection thresholds we derive are reasonably low and observable, but identification is directly impacted by the photosphere absorptions that can greatly hinder identification. For some photostellar types, it could remain feasible. Coupled with the Estimated Time Calculator of the Space Telescope Science Institute, SynthIceSpec can be used to find the optimum observational setup for new observations.

Upcoming radio surveys will probe the sky with unprecedented depth and sky coverage, enabling a broad range of cosmological and astrophysical applications, as well as powerful synergies with experiments at other wavelengths. The preparation and scientific exploitation of these surveys require realistic mock catalogues that capture the complexity of the radio sky and the interplay of its emitting components. We present a modular and extensible algorithm for generating empirical simulations over the full radio sky, i.e. a solid angle of $4\pi$ steradians ($f_{\rm sky}=1$), down to redshift $z=5$, comprising both radio continuum and line emission. The framework combines a simulated dark-matter light-cone with empirically sampled galaxy populations and a probabilistic galaxy-halo assignment scheme, producing self-consistent mock catalogues including multiple radio populations on the same light-cone. We release two public catalogues: a shallow catalogue, fully constrained by existing observational data and limited to flux thresholds of $S_\text{1.4 GHz}^\text{lim} \sim 8\times10^{-5}\ \text{Jy}$ at $1.4\ \text{GHz}$ and $S_\text{21}^\text{lim} \sim 2\ \text{Jy}\cdot\text{Hz}$ for the HI 21 cm line; and a deep catalogue extending the calibrated empirical model to better sensitivities, broadly matching future SKAO surveys, with flux limits of $S_\text{1.4 GHz}^\text{lim} \sim 4\times10^{-5}\ \text{Jy}$ and $S_\text{21}^\text{lim} \sim 0.3\ \text{Jy}\cdot\text{Hz}$. The catalogues include radio continuum active galactic nuclei and star-forming galaxies, together with HI-emitting galaxies, for a total of more than 260 million sources in the shallow catalogue and more than 1 billion in the deep catalogue. We validate the simulations by analysing their statistical properties: the mocks reproduce the targeted clustering and population statistics while retaining minimal physical assumptions.

The adiabatic equation of state $P \propto n^{\Gamma}$ describes the pressure evolution of highly collisional, isotropic plasmas in terms of their density, providing a possible closure of the fluid moment hierarchy in the absence of heat fluxes and dissipation. An analogous closure exists for collisionless, magnetised plasmas, whose pressure tensor is anisotropic with respect to the magnetic field, and the closure is therefore double-adiabatic, prescribing the evolution of the parallel and perpendicular pressures in terms of the magnetic-field strength and density. Here, we present a general first-principle formalism to derive adiabatic laws using the symmetries of the system. With this theory we recover the adiabatic equation of state $P \propto n^{\Gamma}$ for isotropic plasmas and the double-adiabatic equations of state for collisionless, magnetised plasmas. We extend the latter to the relativistic regime, finding that their exact functional form depends on the pressure anisotropy and is not a simple power law. Our double-adiabatic equations of state describe simple geometries, like magnetic mirrors or compressed homogeneous plasmas, as well as complex high-energy astrophysical processes, such as the evolution of plasmoid structures formed during magnetic reconnection.

Observations with the Canadian Hydrogen Intensity Mapping Experiment (CHIME) have been used to measure the 21 cm intensity mapping auto power spectrum, at $z\sim 1$, over a frequency range from 608.2 MHz to 707.8 MHz at wavenumbers $0.4~h~{\rm Mpc}^{-1} \lesssim k \lesssim 1.5~h~{\rm Mpc}^{-1}$. In this paper, we present the results of two different approaches to interpreting this measurement. In the first approach, we use a parametric power spectrum model to constrain an amplitude parameter, defined as $\mathcal{A}^2_{\rm HI} \equiv 10^6 \Omega_{\rm HI}^2(b^2_{\rm HI}+\langle f \mu^2\rangle)^2$, where $\Omega_{\rm HI}$ is the cosmological density parameter for atomic hydrogen ($\rm HI$), $b_{\rm HI}$ is the linear bias for $\rm HI$, and $\langle f \mu^2\rangle$ incorporates the dominant large-scale impact of redshift-space distortions on the angle-averaged power spectrum. Imposing an additional prior on either $\Omega_{\rm HI}$ or $b_{\rm HI}$, based on values in the literature, allows us to break the pairwise degeneracy between those two parameters. In the second approach, we compare CHIME's measurement with predictions for the power spectrum of $\rm HI$ from the IllustrisTNG simulations, finding that the measurement disagrees with the TNG100 run at $3.1\sigma$ and the TNG300 run at $4.0\sigma$. This disagreement is most likely attributable to the strength of nonlinear redshift-space clustering of $\rm HI$ in the simulations, rather than the total abundance of $\rm HI$, and invites further investigation of the physical processes in the simulations that determine the behavior of $\rm HI$ at nonlinear scales. These results exemplify the ability of 21 cm intensity mapping to provide astrophysical information using measurements at nonlinear scales.

We present joint measurements of the pre- and post-reconstruction power spectra, $P_{\rm pre}$ and $P_{\rm post}$, together with their cross-power spectrum, $P_{\rm cross}$, for the Luminous Red Galaxies (LRGs) in the DESI Data Release 1 (DR1). We jointly analyse these observables with an emulator-based full-shape modeling framework, thereby, for the first time, we extract complementary nonlinear information from the galaxy density field before and after reconstruction in real survey data. Specifically, including $P_{\rm post}$ and $P_{\rm cross}$ in addition to $P_{\rm pre}$ (hereafter $P_{\rm all}$) yields an improvement of approximately $18$-$27\%$ in the $\sigma_8$ constraint in both $\Lambda$CDM and $w$CDM, depending on the redshift bin, relative to the $P_{\rm pre}$-only analysis with the cosmic microwave background distance priors (hereafter CMB). In $w$CDM, the joint CMB+$P_{\rm all}$ analysis can tighten the constraints on $w$ by approximately $5$-$15\%$ across the two LRG redshift bins, compared to the CMB+$P_{\rm pre}$ combination. Further incorporating the Type Ia supernova dataset and comparing the cosmological constraints in $w$CDM from each individual power-spectrum component with those from the full combination, we find that $P_{\rm all}$ consistently provides the tightest constraints. From the joint CMB+$P_{\rm all}$+DES-Dovekie dataset, we obtain $\Omega_m = 0.314 \pm 0.0048$ and $w = -0.988 \pm 0.023$ for the \texttt{LRG1} sample, and $\Omega_m = 0.318 \pm 0.0046$ and $w = -0.988 \pm 0.025$ for \texttt{LRG2}. These results demonstrate that combining pre- and post-reconstruction power spectra with their cross-correlation enables DESI to harvest additional nonlinear information, leading to tighter constraints on cosmological parameters.

Characterizing the surface and atmosphere of Earth-like planets in reflected light is a key goal for upcoming direct imaging surveys. NASA's next flagship-class astrophysics mission concept, the Habitable Worlds Observatory (HWO), is a space-based Ultraviolet/Optical/Near-Infrared observatory with a mission design requirement to reach the $10^{-10}$ contrast necessary to characterize Earth-like planets around Sun-like stars. While reflected light from planetary surfaces provides a unique opportunity to constrain the coverage of surface materials and biopigments, detailed predictions of HWO's ability to retrieve surface fractions are necessary but have not been conducted. Here, we model photon-counting noise from astrophysical, instrumental, and post-processing sources for the HWO Exploratory Analytic Case 5 design equipped with a charge-6 vector-vortex coronagraph. By combining our photon-counting noise with five distinct modern Earth models at quadrature, we simulate single-visit HWO observations and perform spectral retrievals using the open-source code $\texttt{POSEIDON}$ to assess our ability to constrain both the surface and atmospheric composition. We find that degeneracies between planetary radius, surface pressure, surface material, and cloud coverage in reflected-light retrievals can significantly complicate the classification of surface features. These degeneracies can complicate the detection of surface biopigments, such as the chlorophyll-induced red edge on modern Earth. Our work shows that developing concrete strategies for detecting surface features and breaking degeneracies in reflected-light observations of Earth-like planets is a critical priority for mission design and data analysis.

We present a joint Bayesian analysis to impose constraints on the generalized $\alpha$-Starobinsky inflationary model using the high-precision cosmological datasets: Planck, CMB lensing from ACT DR6, and Baryon Acoustic Oscillations (BAO) from DESI DR2. For the parameter inference, we introduce an alternative sampling approach. Rather than imposing priors on the cosmological parameters of the inflationary potential $(V_0, \, \alpha, \, N_*)$, we place priors directly on the primordial physical observables $(A_s,\, n_s,\, r)$ through using analytical slow-roll consistency relations, our pipeline internally maps these sampled observables to the corresponding $\alpha$-Starobinsky parameters. These values are then passed to a modified version of $\tt{CLASS}$, which solves the exact inflationary dynamics fully numerically. This pipeline ensures that the final reported posteriors for the observables are computed exactly, completely free from the slow-roll approximation. Applying this methodology, we explore the viability of the $\alpha$-Starobinsky model. We show that, when the full combined dataset is considered, the pure Starobinsky model (i.e., the canonical limit $\alpha = 1$) faces an apparent discrepancy: it requires a large number of $e$-folds of inflation after horizon crossing ($N_* > 60$) due to the shift in the scalar spectral index, $n_s$. In contrast, allowing the deformation parameter $\alpha$ as a free parameter yields a clear $1\sigma$ preference for $\log_{10} \alpha > 0$ present across all datasets. By favoring a broader inflationary plateau, the $\alpha$-Starobinsky model elegantly reconciles theoretically sound expansion histories with empirical data. Notably, we also show that the addition of ACT DR6 lensing data introduces no significant impact on these primordial constraints, confirming that our robust posteriors are primarily driven by Planck and DESI measurements.

We present profile-likelihood constraints on velocity-independent dark matter-proton scattering, including cases in which only a fraction of dark matter has such non-gravitational interactions. Frequentist profile-likelihood techniques provide prior-independent constraints, circumventing prior-volume effects that we show arise in Bayesian constraints on this model. In the limit where the scattering cross section or the fraction of interacting dark matter approaches zero, the other interacting dark matter model parameters become unconstrained, causing the posterior distribution to favor that region of parameter space. Using Planck 2018 cosmic microwave background anisotropy data, we find a clear impact of prior-volume effects on the posteriors used to place constraints on dark matter scattering. Compared to the frequentist analysis, the Bayesian method consistently overestimates the constraints on the cross section. Given the potentially biased upper limits on models subject to prior-volume effects, such as this one, we recommend supplementing Bayesian constraints with frequentist statistics to better assess the impact of priors.

We introduce a minimal two-parameter formulation of the dark energy (DE) density evolution normalized to its present-day value, $f_{\rm DE}(z) \equiv \rho_{\rm DE}(z)/\rho_{\rm DE,0}$, in terms of $f_p\equiv f_{\rm DE}(z_p)$ and the DE equation of state $w_p\equiv w(z_p)$, at a pivot redshift $z_p$. This provides an alternative framework for assessing the evidence for evolving DE, complementary to the established Chevallier-Polarski-Linder (CPL) parameterization. By parameterizing the DE density directly, the $(w_p,\,f_p)$ formulation avoids the approximate degeneracies intrinsic to the $(w_0,\,w_a)$ basis -- in particular the weak sensitivity of the expansion history to $w_a$ -- while reproducing the background evolution of representative quintessence models with equivalent accuracy. Confronting it with the latest baryon acoustic oscillation (BAO) measurements from DESI, a prior on early-universe parameters from Planck cosmic microwave background (CMB) observations, and Type Ia supernovae (SNe) data, we find that the $w_p$ and $f_p$ parameters are both tightly constrained and sensitive to distinct subsets of the data. Specifically, $w_p$ is measured to percent-level precision by BAO and CMB alone, while $f_p$ is pinned down by the independent matter density constraint that only SNe provide. Including the Pantheon+ SNe sample, we obtain $w_p = -1.04 \pm 0.04$ and $f_p = 1.07 \pm 0.04$, with similar results when using the DESY5 SNe sample. The preference for evolving DE over $\Lambda$CDM remains below $3\sigma$ across all dataset combinations, comparable to that obtained with CPL. Notably, the proximity of both $w_p$ and $f_p$ to their cosmological constant values of $(-1,1)$ -- precisely at the epoch where the data are most sensitive -- deepens the coincidence previously identified in the CPL framework, reinforcing the case for caution in interpreting the current evidence for dynamical DE.

We show that dynamical dark energy theories can imprint $O(1)$ modifications on the quasi-normal mode (QNM) spectrum characterising black hole ringdown. The time dependence of dynamical dark energy naturally gives rise to cosmological 'hair' around a black hole. Taking the cubic Galileon as a concrete example, which admits the only known stable solution of this kind, we parametrically connect the cosmological and black hole regimes, derive the induced QNM shifts and forecast the resulting dark energy constraints. We find that the dark energy field profile can be constrained with an accuracy of up to $10^{-2}$ for LVK and $10^{-4}$ for LISA.

Motivated by theoretical and observational developments of cosmological coupled black holes, we construct an exact analytical solution for a black hole immersed in an anisotropic dark sector background, adopting the framework established by [Cadoni et al., JCAP 03 (2024) 026]. By generalizing a static seed metric to a dynamical FLRW background, we derive a solution where the black hole mass co-evolves with the cosmic expansion. We then obtain the explicit form of the radius-dependent coupling exponent, revealing that the interaction is governed by the dark halo profile. Considering the ubiquity of the dark halos surrounding supermassive black holes, our model provides a potential realization of cosmological coupling, interpreting the mass growth as the dynamical response of the surrounding dark sector fluid to the Hubble flow, distinct from the method of modifying the black hole's internal equation of state.

We investigate the imprint of post-inflationary evolution on primordial cosmological correlators by relaxing the standard assumption of instantaneous reheating. We compute the power spectrum and bispectrum for a conformally coupled and non-minimally coupled ($\xi \neq \frac{1}{6},0$) scalar field with cubic self-interaction, across a finite reheating epoch characterized by an effective equation of state $w$, and reheating temperature $T_{\rm reh}$. We find that for a conformally coupled scalar field, the signatures of reheating, such as the modification to oscillatory features, are strictly confined to the subhorizon regime. On superhorizon scales, the correlators receive negligible corrections at leading order. In sharp contrast, for a non-minimally coupled field, we identify a distinct phenomenology where non-trivial values of $\xi$ and reheating equation of state $w$, can induce a tachyonic enhancement of the field modes on superhorizon scales. This enhancement modifies both the power spectrum and bispectrum, marking a sharp deviation from the instantaneous transition limit. Our results demonstrate that to the leading order while conformal coupling shields large-scale correlations from the expansion history, non-minimal coupling breaks this degeneracy and opens an observable window into the physics of reheating.

The Bondi--van der Burg--Metzner--Sachs (BMS) frame of gravitational waves produced by numerical relativity (NR) simulations is crucial for building accurate waveform models. A proper comparison of NR waveforms with other models requires fixing the arbitrary BMS frame. In this work we improve the center-of-mass (CoM) frame fixing for quasicircular, nonprecessing binary systems. Past work approximated the CoM motion with just a linear fit. We compute a post-Newtonian result of the boosted CoM charge to also capture its physical out-spiraling oscillations. We show that using the analytical results improves the robustness of the fit parameters -- translation and boost vectors -- to the choice of duration and time of the fitting window. Our analysis demonstrates a maximum improvement in robustness when the window is placed at the center of the inspiral. We quantified this improvement by computing the ratio of variances of fit parameters when the fit window size is varied. The largest improvement in robustness of parameters is by a factor of $\sim 25$ for the boost vector and $\sim 20$ for the translation vector. Finally, we incorporate this method into the BMS frame-fixing routine of the python package $\texttt{scri}$ for waveforms produced with Cauchy-characteristic evolution.

Classically scale-invariant (and perturbative) theories provide a way to understand large hierarchies, as scales are generated through dimensional transmutation. They always lead to first-order phase transitions, since symmetries are radiatively broken, and they generically feature quasi-flat potentials, which are suitable for inflation. We construct a simple but fully realistic model of this kind that accounts for all observational evidence of new physics and is remarkably compatible with the most recent constraints on inflationary observables from both the Planck/BICEP/Keck and the Atacama Cosmology Telescope (ACT) collaborations. This model illustrates how classical scale invariance generically leads to a non-standard cosmology in which inflation occurs in two stages: a slow-roll stage and a thermal stage, separated by a radiation-dominated era.

We consider a Universe with a three-torus topology which before inflation is devoid of any matter or radiation. Its pre-inflationary evolution is driven solely by Casimir energies of the existing fields, with a radiation-like equation of state. We show that, quite strikingly, with the appropriate number of fermions and bosons in the theory, such a Universe evolving from Planck size at Planck time has the correct critical energy density today. Moreover, assuming typical parameters for inflation and reheating, the three-torus Universe can be several Hubble radii across at present time, which is precisely the size hinted by the anomaly in the cosmic microwave background observed at low multipole moments. We derive a relation between the size of the Universe, the number of e-folds of inflation, and the energy density decrease during reheating, which we then use to determine the parameter values consistent with the anomaly and the lower bound on the size of the Universe from Planck.

We propose that Quintom dark energy, the simplest framework allowing crossing of the cosmological-constant boundary, admits a natural UV completion in a 5D anisotropic orbifold lattice: the Non-Perturbative Gauge-Higgs Unification (NPGHU) model. In this setup, a bulk 5D SU(2) gauge field projects on the 4D boundary to a complex scalar and a U(1) gauge field, identified with the dynamical dark-energy sector, while the Standard Model and dark matter remain localized in four dimensions. At late times, bulk-induced dimension-6 higher-derivative operators generate both physical and phantom scalar and gauge degrees of freedom. We show that the resulting 4D effective action is a modified Quintom model whose background equation of state can naturally realize Quintom-B behavior. A crucial contribution arises from the massive gauge ghost, allowing an excellent fit to DESI data with negligible fine-tuning, unlike standard Quintom scenarios. We further show that the inherited properties of the NPGHU construction e.g. absence of fundamental ghost instabilities, absence of potential terms and a finite low-energy cutoff $\Lambda$ associated with approximate Lorentz invariance, play a central role in the consistency of the effective theory under linear perturbations and vacuum decay. For the most natural regime, $\Lambda \approx {\cal O}(10)H_0$, the model remains robust despite the presence of IR phantom modes. Our results provide a natural and predictive framework in which Quintom dark energy can be consistently embedded in a fundamental theory.

We propose a novel cogenesis framework based on the type-I Dirac seesaw mechanism. The minimal type-I Dirac seesaw with three heavy vector like fermions $(N)$, one singlet scalar $(\eta)$ and the right-handed counterparts $(\nu_R)$ of the Standard Model (SM) neutrinos is extended to include a Dirac fermion dark matter (DM) $(\chi)$ and its heavier scalar companion ($\phi$). The out-of-equilibrium decays of the vector-like fermion generate asymmetries simultaneously in the visible sector, through decay channels involving $(\nu_R,\eta)$ or lepton, Higgs doublets in the SM, and in the dark sector via decaying into $(\chi,\phi)$. The resulting lepton asymmetry is partially converted into the observed baryon asymmetry by electroweak sphaleron processes, while the dark-sector asymmetry survives to constitute the present-day asymmetric DM relic. The generation of asymmetries in multiple sectors and their mutual washouts provide rich dynamics while also keeping the model testable at different observations involving DM, neutrinos, cosmic microwave background (CMB), as well as gravitational waves (GW). We find that successful cogenesis can be realized for DM masses in the range $100~\mathrm{MeV} \lesssim m_\chi \lesssim 39~\mathrm{TeV}$. The lower bound arises from the requirement that the symmetric component of DM annihilates efficiently before the big bang nucleosynthesis (BBN) epoch, while the upper bound is set by unitarity constraints on the asymmetric DM.

We study primordial black hole (PBH) formation in a minimally coupled $f(T)$ teleparallel cosmology that generates a transient departure from standard radiation domination. The model is constructed so that modified-gravity effects are negligible at early and late times, but become dynamically relevant over a finite epoch, during which an effective torsion component reduces the total equation-of-state parameter below 1/3.We show that this transient softening lowers the collapse threshold for overdensities at horizon re-entry, leading to an exponential enhancement of PBH formation. In addition, the modified background alters the relation between temperature and horizon mass, producing a localized feature in the PBH mass function. For representative parameters, PBHs with asteroid-scale masses can account for a significant fraction, or even the entirety, of dark matter for perturbation amplitudes $\sigma^2 \sim \mathcal{O}(10^{-3})$, while remaining consistent with current constraints. Our results demonstrate that modified gravity alone can efficiently generate PBHs, without requiring ad hoc modifications of the radiation sector.

In this work, we investigate the critical behavior of photon rings in the Kerr-Bertotti-Robinson spacetime, describing a rotating black hole immersed in a background magnetic field. We analyze the radial and angular motions of photons under the small magnetic field approximation. Focusing on unstable spherical orbits, we determine three key parameters, $\gamma$, $\delta$, and $\tau$, which characterize radial compression, azimuthal advancement, and time delay. We then examine how these parameters depend on the black hole spin, magnetic field strength, and observer inclination for both on-axis and off-axis observers, and we further analyze the properties of higher-order images through near-critical lens equations. The results show that the magnetic field modifies the geodesic structure, and leads to observable changes in the fine structure of photon rings, providing a useful framework for probing magnetized black hole environments.

A high-granularity telescope system with a large sensitive area and low material budget has been developed for high-energy heavy ion beam tests. The telescope consists of nine layers of silicon microstrip detectors (SSDs), whose performance was validated through a heavy ion beam test at the CERN SPS. A hybrid machine learning algorithm is proposed to address the challenges of nuclear charge measurement with SSDs. The system achieves a spatial resolution of $\mathcal{O}(1) \,$\SI{}{\micro\metre} and a charge resolution better than 0.16 charge units for nuclei from $Z = 1$ to $Z = 29$, with a sensitive area of $8 \times 8 \, \mathrm{cm}^2$. To the best of our knowledge, this represents the most precise charge and spatial resolution simultaneously achieved by a silicon telescope to date.

The ratio between the gravitational waves (GW) and electromagnetic waves (EMW) luminosity distance ratio is a key observable that allows to test the nature of gravity, using gravitational waves emitted from compact binary coalescences. We develop a new non parametric method for constraining the GW-EMW distance ratio, in order to perform model independent analysis of observational data, not based on any specific theoretical of phenomenological assumption. We apply the method to the analysis of binary black hole mergers data from the GWTC-3 catalogue, performing a joint analysis of cosmological and population parameters. The results are consistent with general relativity and with previous analyses based on parametric methods.

The astrophysical $S$-factor for the proton-proton fusion is calculated in the low-energy regime for a variety of nuclear interactions and consistent nuclear currents, derived within chiral effective field theory. We estimate, for the first time, the theoretical uncertainty on the $S$-factor due to the truncation of the chiral expansion of the currents using a Bayesian analysis. In order to reach an accuracy at the percent level in the calculation, the electromagnetic potential includes contributions beyond the leading Coulomb interaction, such as two-photon exchange and vacuum polarization. The initial proton-proton state is expanded in partial waves and only the ${}^1S_0$ contribution is included, as it is known that the other partial-waves effects are negligible. The low-energy constant entering the contact term in the weak axial current operator is calibrated to reproduce the Gamow-Teller matrix element in Tritium $\beta$-decay. The value $S(0)$ is found to be $S(0)=(4.068 \pm 0.025)\times 10^{-25} \: \text{MeV}\: \text{b}$.

The magnetospheres of magnetars, a class of highly magnetized neutron stars, host magnetic fields exceeding the Schwinger limit, where Quantum Electrodynamics (QED) becomes nonlinear. In such environments, QED scattering processes are strongly modified, which may affect plasma dynamics. In this work, we apply a formalism originally developed for the study of magnetic-field effects in hot quark-gluon plasma to strong-field QED. The method resums interactions between virtual electrons and the external magnetic field, consistently incorporating the finite decay widths of excited Landau levels derived from the fermion self-energy. Using this framework, we perform the first systematic analysis of tree-level QED scattering processes in strong magnetic fields, concentrating on the processes of highest relevance for the plasma dynamics of magnetars. All resulting cross sections are provided in an open-source Python package.

We revisit the double adiabatic evolution equations and extend them to the relativistic and ultrarelativistic regimes. We analytically solve the relativistic, time-dependent drift kinetic equation for a homogeneous, magnetized, collisionless plasma and obtain a solution explicitly dependent on the magnetic field and density variations. In the case of an initial relativistic Maxwellian distribution, a natural extension to an anisotropic Maxwell-Jüttner is obtained. We calculate the moments of this time-dependent solution and obtain analytical expressions for the evolution of the perpendicular and parallel pressures in the ultrarelativistic case. We numerically solve the moment equations in the relativistic case and obtain general expressions for the double-adiabatic equations in this regime. We confirm our results using fully kinetic particle-in-cell simulations of shearing and compressing boxes. Our findings can be readily applied to relativistic species including cosmic-rays and electron-positron pairs, present in astrophysical plasmas like pulsar wind nebulae, astrophysical jets, black hole accretion flows, and Van Allen radiation belts.