Papers for Wednesday, Apr 30 2025

Icy satellites host topography at many length scales, from rifts and craters on the small end to equatorial-pole shell thickness differences that are comparable to these bodies' circumference. The rate of topographic evolution depends on the rheology of the ice shell, the difference in density between the shell and potential global subsurface oceans, the ice shell thickness distribution, as well as the rate of tidal heat dissipation in the shell. Here we analyze relaxation of topography starting from the Stokes equations for viscous fluid flow. We write out the shallow-shell models (with and without a rigid crust) and compare the resulting time scales to previous work. Using three-dimensional linearized flow for small-amplitude topographic perturbations, we recover the time scales of relaxation from the shallow models as well as shells that are not shallow. For a shell of constant viscosity, we find that the topographic relaxation time scale is constant for wavelengths larger than the ice thickness, an idea that runs contrary to the existing paradigm. For a shell with a viscosity that decreases exponentially with depth, we show numerically that there is a regime where the larger viscosity outer crust acts as a nearly rigid boundary. In this case, the relaxation time scale depends on the wavelength. For the largest spatial scales, the time scale becomes independent of wavelength again and the value is set by the average shell viscosity. However, the spatial scale that this transition occurs at becomes larger as the viscosity contrast increases, limiting the applicability of the scale-independent relaxation rate. These results for the relaxation of topography have implications for interpreting relaxed crater profiles, inferences of ice shell thickness from topography, and upcoming observations from missions to the outer solar system.

Field-level inference has emerged as a promising framework to fully harness the cosmological information encoded in next-generation galaxy surveys. It involves performing Bayesian inference to jointly estimate the cosmological parameters and the initial conditions of the cosmic field, directly from the observed galaxy density field. Yet, the scalability and efficiency of sampling algorithms for field-level inference of large-scale surveys remain unclear. To address this, we introduce a standardized benchmark using a fast and differentiable simulator for the galaxy density field based on \texttt{JaxPM}. We evaluate a range of sampling methods, including standard Hamiltonian Monte Carlo (HMC), No-U-Turn Sampler (NUTS) without and within a Gibbs scheme, and both adjusted and unadjusted microcanonical samplers (MAMS and MCLMC). These methods are compared based on their efficiency, in particular the number of model evaluations required per effective posterior sample. Our findings emphasize the importance of carefully preconditioning latent variables and demonstrate the significant advantage of (unadjusted) MCLMC for scaling to $\geq 10^6$-dimensional problems. We find that MCLMC outperforms adjusted samplers by over an order-of-magnitude, with a mild scaling with the dimension of our inference problem. This benchmark, along with the associated publicly available code, is intended to serve as a basis for the evaluation of future field-level sampling strategies. The code is readily open-sourced at this https URL

We explore the potential of optical intensity interferometry to extract angularly resolved information from supernova explosions, introducing the "expanding ejecta method" (EEM) as a robust alternative to the classical expanding photosphere method (EPM). Foreseeing future improvements to intensity interferometers of large light collection area ($25\pi\,\rm{m}^2$ per telescope) equipped with spectral multiplexing ($10^4$ spectral resolution) and fast photodetectors ($10\,\mathrm{ps}$ timing resolution, $50\%$ overall efficiency), we demonstrate that high signal-to-noise measurements of the visibility modulus are achievable for Type IIP (Type Ia) supernovae out to $3~(12)\,\mathrm{Mpc}$. By focusing on generic line emission and absorption in ballistic ejecta, the EEM can relax assumptions about spherical symmetry, blackbody radiation, and extinction. The EEM enables angular diameter distances to be determined with $\sim2\%$ precision for supernovae of apparent magnitude $m = 12$ from a 60-hour observation by an intensity interferometer with those same instrumental specifications. We argue that the EEM is significantly more robust to modeling uncertainties and systematic effects than (variants of) the EPM. In a companion paper, we show how the EEM can be used to provide geometric anchors for cosmic distance ladder calibration, or to construct a wholly independent Hubble diagram based on angular diameter distances.

The large-scale morphology of Milky Way (MW)-mass dark matter (DM) halos is shaped by two key processes: filamentary accretion from the cosmic web and interactions with massive satellites. Disentangling their contributions is essential for understanding galaxy evolution and constructing accurate mass models of the MW. We analyze the time-dependent structure of MW-mass halos from zoomed cosmological-hydrodynamical simulations by decomposing their mass distribution into spherical harmonic expansions. We find that the dipole ($\ell=1$) and quadrupole ($\ell=2$) moments dominate the gravitational power spectrum, encoding key information about the halo's shape and its interaction with the cosmic environment. While the dipole reflects transient perturbations from infalling satellites and damps on dynamical timescales, the quadrupole -- linked to the halo's triaxiality -- is a persistent feature. We show that the quadrupole's orientation aligns with the largest filaments, imprinting a long-lived memory on the halo's morphology even in its inner regions ($\sim30$ kpc). At the virial radius, the quadrupole distortion can reach 1-2 times the spherical density, highlighting the importance of environment in shaping MW-mass halos. Using multivariate Singular Spectrum Analysis, we successfully disentangle the effects of satellite mergers and filamentary accretion on quadrupole. We find that the quadrupolar response induced by LMC-mass satellites has an order of magnitude larger gravitational power than in spherically symmetric MW models. This highlights the need for models that incorporate the MW's asymmetry and time-evolution, with direct consequences for observable structures such as disk warps, the LMC-induced wake, and stellar tracers -- particularly in the era of precision astrometry.

Aims. Outflows from active galactic nuclei are invoked as the principal feedback process regulating the co-evolution of supermassive black holes and their host galaxies. Because of their multi-phase and multi-scale nature, an exhaustive description of these winds should exploit multiple tracers. However, connecting various outflow features remains a challenge. The aim of this work is to provide a complete characterisation of outflows in a sample of z$\sim$5 quasars, by exploiting the combination of different emission and absorption tracers. Methods. We analysed the UV/optical and FIR continuum, line emission, and absorption in a sample of 39 z$\sim$5 quasars observed with VLT/X-Shooter and ALMA (available for six objects). We identified broad and narrow absorption lines associated with the quasar and emission lines to determine black hole masses and bolometric luminosities. Results. Our sample encompasses massive (log($M_{\text{BH,MgII}}/M_\odot$) = 8.5-10) and luminous (log($L_{\text{bol}}$/(erg/s)) = 46.9-48) quasars at redshift 5-5.7. They display powerful ionised outflows detected in both emission and absorption, with velocities exceeding 48,000 km/s in some cases, and lie above the local black hole - host galaxy mass relation, exhibiting a behaviour similar to that of z$\gtrsim$6 quasars. These findings suggest a phase of efficient black hole feedback occurring at redshift z$\gtrsim$6 and likely persisting down to z$\sim$5, characterised by rapid black hole growth exceeding that of the host galaxy. The fraction of quasars with outflow detections in absorption is higher for larger CIV-MgII velocity shifts, suggesting that while the physical mechanisms powering the two outflow phenomena detected in emission and absorption may differ, a correlation exists between them.

The evolution of galaxies is significantly influenced by the environments they inhabit. While high-density regions, such as clusters of galaxies have been widely studied, the dynamics and quenching processes in intermediate environments remain less explored. These systems provide a valuable context for understanding the transition of galaxies from active star formation to quiescence. This study aims to characterise the astrophysical properties of galaxies in intermediate-mass galaxy groups $(13.5 \leq \log(M_{200}/M_{\odot}) \leq 13.7)$, with a focus on their evolutionary pathways and the key processes driving their transition through the green valley (GV) and green zone (GZ). Using the Illustris TNG300-1 hydrodynamical cosmological simulations, we classified galaxies based on their trajectories and environment into five categories: group galaxies (GRs), backsplash galaxies (BSs), recent infallers (RINs), infall galaxies (INs), and field galaxies (FGs). We examined their optical colours in the $(u-r)$-stellar mass diagram, specific star formation rates (sSFRs), gas fractions, and stellar mass evolution from $z=0.5$ to $z=0$. At $z=0$, FGs dominate the blue cloud, while GRs show progressive reddening. BSs exhibit the highest fraction of green galaxies, highlighting their transitional nature, whereas RINs show a rapid quenching upon entering $R_{200}$. Our results reveal that the timing of group entry and environmental effects, such as gas depletion, play a critical role in galaxy quenching. Green BS and RINs follow distinct evolutionary tracks, with the latter undergoing more rapid changes due to later infall into more massive systems.

Understanding the post-merger evolution of binary neutron star merger (BNSM) requires accurate modeling of neutrino transport and microphysics including neutrino flavor conversions. Many previous studies have suggested that fast flavor instability (FFI) and collisional flavor instability (CFI) pervade inner regions of BNSM remnant, and they could impact on fluid dynamics and r-process nucleosynthesis. In this work, we re-examine prospects of occurrences of FFI and CFI using Boltzmann neutrino transport, assuming a frozen fluid background obtained from a numerical relativity simulation of BNSM. We pay special attention to a case involving a long-lived ($>1\,$ s) hypermassive neutron star (HMNS). Apart from confirming the claim that these flavor instabilities can occur in BNSM remnants, some new insights are revealed. We identify multiple mechanisms responsible for generating electron neutrino lepton number angular crossings, corresponding to a key indicator of FFI onset, which differ notably from those in black hole (BH) accretion disk systems. We also find that the appearance of positive chemical potential of electron-type neutrinos is a necessary, though not sufficient, condition for FFI, which offers a valuable diagnostic tool for predicting FFI based solely on fluid properties. For CFI, their growth rates are generally lower than FFI, but they can persistently occur in most of the accretion disk up to $\sim 1\,$ s. We also find that neglecting contributions of heavy-leptonic neutrinos results in overestimating growth rate and area of unstable regions of CFI. Our result suggests that FFI (CFI) tends to occur transiently (persistently) and locally (widespread in the disk), and FFI is more sensitive to the central compact object (HMNS or BH) than CFI, though more self-consistent simulations with incorporating effects of flavor conversions are needed to confirm these claims.

At the Galactic Centre, we can expect a population of a few tens of early extreme-mass ratio inspirals (E-EMRIs) and extremely large mass ratio inspirals (XMRIs). Depending on their evolutionary stage, they can be highly eccentric, with moderate signal-to-noise ratios (SNRs) of tens or hundreds, or nearly circular, with SNRs as large as a few thousand. Their individual signals combine into a common signal, which can complicate the resolution of other types of sources. We have calculated the foreground signal of continuous E-EMRIs and XMRIs using a catalog based on the expected number of sources and a realistic phase-space distribution. The forest of E-EMRIs will cover a large portion of the LISA sensitivity curve, obscuring the signals of some massive black hole binaries, verification binaries, and harmonics of EMRIs in their polychromatic phase. The combined signal from XMRIs will be much weaker but still affect intermediate-mass black hole binaries. Due to the large SNR, this forest can be also found in other galactic nuclei, such as that of the Andromeda galaxy. Even under conservative assumptions, the forest created by E-EMRIs and XMRIs in our Galactic Centre will likely pose a challenge for resolving other types of sources, as their contribution is non-coherent and exhibits large SNRs.

We present the implementation of a fully time-dependent relativistic disk model - based on the light curve fitting package FitTeD - into the X-ray spectral fitting environment, pyXspec. This implementation enables simultaneous fitting of multi-epoch and multi-wavelength spectral data, where the only free parameters are those describing the black hole and the initial conditions, while the subsequent evolution is fully governed by the dynamical equations of an evolving accretion flow. We apply the model to fit seven epochs of X-ray spectra and two epochs of UV spectra of the 'long-lived' tidal disruption event (TDE) and quasi-periodic eruption (QPE) source GSN 069, spanning from its discovery in 2010 through late-2019. Our results show that such 'long-lived', X-ray-bright TDEs - of which GSN 069 is a prime, but not unique, example-can be naturally explained within the same framework as events with shorter-lived X-ray emission, such as ASASSN-14li and AT2019dsg. Their key distinction lies in the 'viscous' timescale parameter - tied to the disk's stress and angular momentum transport efficiency - which should be treated as a free parameter when modeling the disk evolution of transient events. We examine the implications for QPE models by tracking the time evolution of disk properties such as mass surface density and accretion rate. We show that existing QPE models cannot reproduce the observed connection between the presence (e.g., 2018) or absence (e.g., 2014) of eruptions and disk properties. In the context of orbiter-disk collision models, the change in mass surface density appears insufficient to explain the needed variation in the eruption's temperature. The absence of eruptions in GSN 069 in 2014 remains a major challenge for QPE models.

We report the analysis of four unambiguous planets and one possible planet from the subprime fields ($\Gamma \leq 1~{\rm hr}^{-1}$) of the 2017 Korea Microlensing Telescope Network (KMTNet) microlensing survey, to complete the KMTNet AnomalyFinder planetary sample for the 2017 subprime fields. They are KMT-2017-BLG-0849, KMT-2017-BLG-1057, OGLE-2017-BLG-0364, and KMT-2017-BLG-2331 (unambiguous), as well as KMT-2017-BLG-0958 (possible). For the four unambiguous planets, the mean planet-host mass ratios, $q$, are $(1.0, 1.2, 4.6, 13) \times 10^{-4}$, the median planetary masses are $(6.4, 24, 76, 171)~M_{\oplus}$ and the median host masses are $(0.19, 0.57, 0.49, 0.40)~M_{\odot}$ from a Bayesian analysis. We have completed the AnomalyFinder planetary sample from the first 4-year KMTNet data (2016--2019), with 112 unambiguous planets in total, which nearly tripled the microlensing planetary sample. The ``sub-Saturn desert'' ($\log q = \left[-3.6, -3.0\right]$) found in the 2018 and 2019 KMTNet samples is confirmed by the 2016 and 2017 KMTNet samples.

Exoplanets classified as super-Earths are commonly observed on short period orbits, close to their host stars, but their abundance on wider orbits is poorly constrained. Gravitational microlensing is sensitive to exoplanets on wide orbits. We observed the microlensing event OGLE-2016-BLG-0007, which indicates an exoplanet with a planet-to-star mass ratio roughly double the Earth-Sun mass-ratio, on an orbit longer than Saturn's. We combine this event with a larger sample from a microlensing survey to determine the distribution of mass ratios for planets on wide orbits. We infer there are $\sim 0.35$ super-Earth planets per star on Jupiter-like orbits. The observations are most consistent with a bimodal distribution, with separate peaks for super-Earths and gas giants. We suggest that this reflects differences in their formation processes.

We compare the self-similar baryonic mass fraction scaling relations between galaxy clusters from the South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) survey and the IllustrisTNG state-of-the-art magnetohydrodynamical cosmological simulations. Using samples of 218 (TNG100) and 1605 (TNG300) friends-of-friends (FoF) haloes within $0.0 \leq z \leq 1.5$ and $M_{200c} \geq 7 \times 10^{13} M_{\odot}$, we fit the scaling relations using Simple Power Law (SPL), Broken Power Law (BPL), and General Double Power Law (GDPL) models through non-linear least squares regression. The SPL model reveals null slopes for the baryonic fraction as a function of redshift, consistent with self-similarity. Observations and simulations agree within $1{-}2\sigma$, suggesting comparable baryonic scaling slopes. We identify $\sim$13.8$-$14.1 per cent of baryons as "missing", primarily in the form of intracluster light (ICL) across all halo masses and warm gas in low-mass haloes. High-mass haloes ($\log_{10}(M_{500c}/M_{\odot}) \geq 14$) adhere to self-similarity, while low-mass haloes exhibit deviations, with the breakpoint occurring at $\log_{10}(M_{500c}/M_{\odot}) \sim 14$, where baryons are redistributed to the outskirts. Our findings suggest that the undetected warm-hot intergalactic medium (WHIM) and baryon redistribution by feedback mechanisms are complementary solutions to the "missing baryon" problem.

The variation of metal production over time and its dilution in the interstellar medium depend on the star formation and gas accretion rates. Measuring age-chemistry relations across the Milky Way disk provides key constraints on the gas accretion and star formation histories, and offers insight into the birth locations of stars. We present a study based on a sample of nearly 30000 dwarf stars from the APOGEE DR17 survey within 2 kpc of the Sun, for which we measure accurate ages. Various parameter combinations are tested to optimize stellar age determination from isochrones. The resulting age-chemistry relations for a selected subsample of 12000 stars are interpreted with the aid of a chemical evolution model. The data reveal a well-defined, tight thick disk sequence, characterized by high [alpha/Fe], subsolar metallicities, and ages older than 8 Gyr. The thin disk, with lower [alpha/Fe] and younger ages, exhibits a wide spread in metallicity at all ages, with apparent structures. Dividing the sample by guiding radius into inner, intermediate, and outer disks shows distinct chemical evolution patterns. The inner disk displays a monotonic, homogeneous chemical evolution with little dispersion, while the outer disk shows little metallicity increase over the past 8 Gyr. The solar neighborhood appears as a mixture, not only due to stellar migration but also because the chemical evolution of the ISM in this intermediate region results from the mixing of gas from the inner and outer disks. In particular, we demonstrate that the solar vicinity experienced a decrease in the mean ISM metallicity 7-9 Gyr ago. A plausible explanation involves a radial inflow of lower-metallicity gas from the outer disk at that time, which diluted the gas leftover by the thick disk formation, contributing to the observed metallicity gradient in the intermediate region.

The COSMOS HI Large Extragalactic Survey (CHILES) Continuum & Polarization (CHILES Con Pol) survey is an ultra-deep continuum imaging study of the COSMOS field conducted using the Karl G. Jansky Very Large Array. We obtained 1000 hours of L-band ($\lambda = 20$ cm) observations across four spectral windows (1.063-1.831 GHz) on a single pointing and produced a confusion limited image with an apparent RMS noise of 1.67 $\mu$Jy beam$^{-1}$ with a synthesized beam of 5$.\!\!^{\prime\prime}$5$\times$5$.\!\!^{\prime\prime}$0. This paper reports a 1.4 GHz radio continuum source catalog containing 1678 sources detected above 7$\sigma$ (flux densities greater than 11.7 $\mu$Jy), identified using two independent source extraction programs applied to the Stokes $I$ image. Resolved sources dominate at flux density S$_{1.4GHz} \ge 42 $\mu$Jy. Radio spectral index for each source was derived using a power-law fit across the four spectral windows, and we found that a robust spectral index measurement requires a total signal-to-noise ratio of at least 20. Comparisons with previous 1.4 GHz radio continuum surveys show good overall consistency, but evidence for a high degree of catalog incompleteness and the effects of source confusion are evident for some of the earlier studies.

We present an X-ray/UV reverberation analysis of NGC 7469 across 210 days, using daily NICER observations with contemporaneous monitoring by Swift UVOT+XRT. We model the X-ray spectrum with a power law continuum and a soft excess during each NICER epoch. These emission sources demonstrate correlated flux variability with a lag consistent with zero days. We find that the power law emission is consistent with a compact X-ray corona, and that the soft excess can be explainedby reflected coronal emission from the inner accretion disk. We test the relationship between changes in the flux of the X-ray corona and the UVW2 continuum, finding strong correlation and a negative X-ray lag of less than one day. This is consistent with a scenario in which the X-ray corona drives the UVW2 light curve through thermal reprocessing

We introduce and describe the CHILES Continuum \& Polarization (CHILES Con Pol) Survey, a 1000 hour 1.4 GHz wideband full polarization radio continuum deepfield with the Very Large Array (VLA), commensurate with the CHILES HI deepfield. We describe the observational configuration, outline the calibration of the data, and discuss the effect of Radio Frequency Interference on different observing epochs. In addition, we present a novel radio continuum imaging strategy, using well known baseline subtraction techniques in radio spectral data, and discuss the applications to the removal of artifacts from sources far from the field center. Additionally, we discuss the nature of a low-level image-wide offset, the so-called ``negative bowl" and simulate our observations in order that we may both properly understand and correct for this artifact. Using these techniques, we present the first total intensity image of the data which achieves an r.m.s. noise of 1.3 $\mu$Jy beam$^{-1}$ with a synthesized beam of 4.5\arcsec x 4.0\arcsec, the most sensitive L-band image ever taken at this resolution. We then place this image into the broader context of 1.4 GHz radio continuum surveys in the literature, in terms of image sensitivity and fidelity, and $\mu$Jy level source counts and P(D) analysis.

Metal-poor M subdwarfs are among the oldest stellar populations and carry valuable information about the chemical enrichment history of the Milky Way. The measurements of chemical abundances of these stars therefore provide essential insights into the nucleosynthesis in the early stages of the Galaxy's formation. We present the detailed spectroscopic analysis of a nearby metal-poor M subdwarf, LHS 174 from its high-resolution optical spectrum, and apply our previously developed spectral fitting code, \texttt{AutoSpecFit}, to measure the abundances of five elements:[O/H]=$-$0.519$\pm$0.081, [Ca/H]=$-$0.753$\pm$0.177, [Ti/H]=$-$0.711$\pm$0.144, [V/H]=$-$1.026$\pm$0.077, and [Fe/H]=$-$1.170$\pm$0.135. We compare the abundances of O, Ti, and Fe derived from this work and those from previous studies and demonstrate the observed data is clearly better matched with the synthetic model generated based on our abundances than those from the other analyses. The accuracy of inferred stellar abundances strongly depends on the accuracy of physical parameters, which motivates us to develop a reliable technique to determine the parameters of low-mass M dwarfs more accurately ever than before and infer abundances with smaller uncertainties.

Two classes of Gamma-ray bursts (GRBs), corresponding to the short/hard and the long/soft events, with a putative intermediate class, are typically considered in the observer frame. However, when considering GRB characteristics in the cosmological rest frame, the boundary between the classes becomes blurred. The goal of this research is to check for the evidences of a third "intermediate" class of GRBs and investigate how the transformation from the observer to the rest frame affects the hardness-duration-based classification. We applied fits with skewed and non-skewed (symmetric) Gaussian and Student distributions to the sample of 409 GRBs with reliably measured redshifts to cluster the bursts on the hardness (Ep) - duration (T90) plane. We found that, based on AIC/BIC criteria, the statistically preferred number of clusters on the GRB rest-frame hardnesses-duration plane does not exceed two. We also assess the robustness of the clustering technique. We did not find any solid evidence of an intermediate GRB class on the rest-frame hardness-duration plane.

Time-domain surveys have advanced astronomical research by revealing diverse variable phenomena, from stellar flares to transient events. The scale and complexity of survey data, along with the demand for rapid classification, present significant challenges for analysis. While machine learning offers solutions, most existing models are tailored to single tasks, struggle to generalize, and depend heavily on large, accurately labeled datasets. We introduce FALCO, a foundation model for astronomical light curve analysis in time-domain astronomy. This work presents the initial version of FALCO trained via self-supervised learning on unlabeled Kepler light curves using a Transformer-based architecture. The model has been evaluated on three distinct tasks and demonstrates strong generalization: achieving 95 percent accuracy in stellar variability classification across eight classes, an overall RMSE of 0.1305 dex in surface gravity estimation (notably improved to below 0.08 dex when log g is less than 1, and approximately 0.02 dex near log g equals 3), and 87 percent precision in flare identification. These results highlight the model's versatility and ability to learn generalizable representations from light curves, enabling straightforward adaptation to diverse tasks. We further analyze the impact of model scaling and sequence length, finding performance improves with larger models and longer input sequences. We also apply FALCO to derive surface gravity (log g) measurements for 179,732 Kepler stars from their light curves.

We present the largest catalog to \modify{date} of triply ionized carbon (CIV) absorbers detected in quasar spectra from the Dark Energy Spectroscopic Instrument. Using an automated matched-kernel convolution method with adaptive signal-to-noise thresholds, we identify $101,487$ CIV systems in the redshift range $ 1.4 < z < 4.8 $ from $300,637$ quasar spectra. Completeness is estimated via Monte Carlo simulations and catalog is $50\%$ complete at $\mathrm{EW}_{\mathrm{CIV}} \geq 0.4$ Angstrom. The differential equivalent width frequency distribution declines exponentially and shows weak redshift evolution. The absorber incidence per unit comoving path increases by a factor of $2-5$ from $ z \approx 4.5 $ to $ z \approx 1.4 $, with stronger redshift evolution for strong systems. Using column densities derived from the apparent optical depth method, we constrain the cosmic mass density of CIV, $\Omega_{\mathrm{CIV}}$, which increases by a factor of $\sim 3.6$ from $ (0.46 \pm 0.01) \times 10^{-8} $ at $ z \approx 4.5 $ to $ (1.65 \pm 0.04) \times 10^{-8} $ at $ z \approx 1.4 $. From $ \Omega_{\rm CIV} $, we estimate a lower limit on intergalactic medium metallicity $ \log(Z_{\rm IGM}/Z_{\odot}) \gtrsim -3.5 $ at $ z \sim 2.3 $, with a smooth decline at higher redshifts. These trends trace the cosmic star formation history and HeII photoheating rate, suggesting a link between CIV enrichment, star formation, and UV background over $\sim 3$ Gyr. The catalog also provides a critical resource for future studies connecting circumgalactic metals to galaxy evolution, especially near cosmic noon.

We present optical observations of Halley type comet 12P/Pons-Brooks on its approach to perihelion. The comet was active even in the first observations at about 8 au. Starting at 4 au, 12P exhibited an extraordinary series of outbursts, in which the brightness changed by a factor up to 100 and the coma morphology transformed under the action of radiation pressure into a distinctive ``horned'' appearance. Individual outburst dust masses are several x 1e9 kg, with kinetic energies 1e14 J, release times 1e4 s and effective power 1e10 W. These properties are most consistent with, although do not definitively establish, an origin by the crystallization of amorphous water ice with the related release of trapped supervolatile gases. This interpretation is supported by the observation that the specific outburst energy and the specific crystallization energy are comparable (both near 1e5 J kg/s.

We report five nearby ($d_{\mathrm{helio}} < 5$ kpc) stellar substructures in the Galactic halo from a subset of 138,661 stars in the Dark Energy Spectroscopic Instrument (DESI) Milky Way Survey Year 1 Data Release. With an unsupervised clustering algorithm, HDBSCAN*, these substructures are independently identified in Integrals of Motion ($E_{\mathrm{tot}}$, $L_{\mathrm z}$, $\log{J_r}$, $\log{J_z}$) space and Galactocentric cylindrical velocity space ($V_{R}$, $V_{\phi}$, $V_{z}$). We associate all identified clusters with known nearby substructures (Helmi streams, M18-Cand10/MMH-1, Sequoia, Antaeus, and ED-2) previously reported in various studies. With metallicities precisely measured by DESI, we confirm that the Helmi streams, M18-Cand10, and ED-2 are chemically distinct from local halo stars. We have characterised the chemodynamic properties of each dynamic group, including their metallicity dispersions, to associate them with their progenitor types (globular cluster or dwarf galaxy). Our approach for searching substructures with HDBSCAN* reliably detects real substructures in the Galactic halo, suggesting that applying the same method can lead to the discovery of new substructures in future DESI data. With more stars from future DESI data releases and improved astrometry from the upcoming {\it Gaia} Data Release 4, we will have a more detailed blueprint of the Galactic halo, offering a significant improvement in our understanding of the formation and evolutionary history of the Milky Way Galaxy.

The Dark Energy Spectroscopic Instrument (DESI) is a massively parallel spectroscopic survey on the Mayall telescope at Kitt Peak, which has released measurements of baryon acoustic oscillations determined from over 14 million extragalactic targets. We combine DESI Data Release 2 with CMB datasets to search for evidence of matter conversion to dark energy (DE), focusing on a scenario mediated by stellar collapse to cosmologically-coupled black holes (CCBH). In this physical model, which has the same number of free parameters as $\Lambda$CDM, DE production is determined by the cosmic star formation rate density (SFRD), allowing for distinct early- and late-time cosmologies. Using two SFRDs to bracket current observations, we find that the CCBH model: accurately recovers the cosmological expansion history, agrees with early-time baryon abundance measured by BBN, reduces tension with the local distance ladder, and relaxes constraints on the summed neutrino mass $\sum m_\nu$. For these SFRDs, we find a peaked positive $\sum m_\nu < 0.149\,\rm eV$ (95\% confidence) and $\sum m_\nu = 0.106^{+0.050}_{-0.069}\,\rm eV$ respectively, in good agreement with lower limits from neutrino oscillation experiments. A peak in $\sum m_\nu > 0$ results from late-time baryon consumption in the CCBH scenario and is expected to be a general feature of any model that converts sufficient matter to dark energy during and after reionization.

The detection of the gravitational wave (GW) event GW190521 raises questions about the formation of black holes within the pair-instability mass gap (PIBHs). We propose that Population III (Pop III) star clusters significantly contribute to events similar to GW190521. We perform $N$-body simulations and find that PIBHs can form from stellar collisions or binary black hole (BBH) mergers, with the latter accounting for 90\% of the contributions. Due to GW recoil during BBH mergers, approximately 10-50% of PIBHs formed via BBH mergers escape from clusters, depending on black hole spins and cluster escape velocities. The remaining PIBHs can participate in secondary and multiple BBH formation events, contributing to GW events. Assuming Pop III stars form in massive clusters (initially 100,000 $M_\odot$) with a top-heavy initial mass function, the average merger rates for GW events involving PIBHs with 0% and 100% primordial binaries are $0.005$ and $0.017$ $\text{yr}^{-1} \text{Gpc}^{-3}$, respectively, with maximum values of $0.030$ and $0.106$ $\text{yr}^{-1} \text{Gpc}^{-3}$. If Pop III stars form in low-mass clusters (initial mass of $1000M_\odot$ and $10000 M_\odot$), the merger rate is comparable with a 100% primordial binary fraction but significantly lower without primordial binaries. We also calculate the characteristic strains of the GW events in our simulations and find that about 43.4% (LISA) 97.8% (Taiji) and 66.4% (Tianqin) of these events could potentially be detected by space-borne detectors, including LISA, Taiji, and TianQin. The next-generation GW detectors such as DECIGO, ET, and CE can nearly cover all these signals.

We report the discovery of a progenitor candidate for the Type IIb SN 2024abfo using multi-epoch pre-explosion images from the Hubble Space Telescope and the Dark Energy Camera Legacy Survey. The progenitor exhibited a ~0.7 mag decline in F814W from 2001 to 2013, followed by significant brightness and color fluctuations in the g, r and z bands. This is the first time that a variable progenitor has been found for a SN IIb. We suggest that the variability is caused by intrinsic changes in the progenitor star instead of varying obscuration by circumstellar dust. Our results show that the progenitor of SN 2024abfo was likely a yellow supergiant star with an initial mass of 12--14, 14--15 and 15--18 Msun for different assumed values of circumstellar reddening of E(B-V)_CSM = 0, 0.1 and 0.2 mag. Our study underscores the critical role of multi-epoch imaging surveys in revealing the final stages of core-collapse supernovae progenitors.

We present the construction and early results of ESPARTACO 2, the new stellar spectrograph built for research and education at the Astronomical Observatory of the Universidad de los Andes in Bogotá, Colombia. This instrument offers several resolutions from 20,000 in first order using a 50 $\mu$m fiber, to 100,000 in second order in the near infrared. Precise radial-velocity measurements are made possible by simultaneous wavelength calibration. Combined with the 40-cm Meade telescope located at our facilities, a limiting magnitude of 6 is reached. This instrument is a considerable improvement over its predecessor in throughput, reliability and ease.

Sub-Neptunes, the most common planet type, remain poorly understood. Their atmospheres are expected to be diverse, but their compositions are challenging to determine, even with JWST. Here, we present the first JWST spectroscopic study of the warm sub-Neptune GJ3090b (2.13R$_\oplus$, Teq~700 K) which orbits an M2V star, making it a favourable target for atmosphere characterization. We observed four transits of GJ3090b; two each using JWST NIRISS/SOSS and NIRSpec/G395H, yielding wavelength coverage from 0.6-5.2 $\mu$m. We detect the signature of the 10833 Å metastable Helium triplet at a statistical significance of 5.5$\sigma$ with an amplitude of 434$\pm$79 ppm, marking the first such detection in a sub-Neptune with JWST. This amplitude is significantly smaller than predicted by solar-metallicity forward models, suggesting a metal-enriched atmosphere which decreases the mass-loss rate and attenuates the Helium feature amplitude. Moreover, we find that stellar contamination, in the form of the transit light source effect, dominates the NIRISS transmission spectra, with unocculted spot and faculae properties varying across the two visits separated in time by approximately six months. Free retrieval analyses on the NIRSpec/G395H spectrum find tentative evidence for highly muted features and a lack of CH4. These findings are best explained by a high metallicity atmosphere (>100x solar at 3$\sigma$ confidence, for clouds at $\sim \mu$bar pressures) using chemically-consistent retrievals and self-consistent model grids. Further observations of GJ3090b are needed for tighter constraints on the atmospheric abundances, and to gain a deeper understanding of the processes that led to its potential metal enrichment.

Recently, the James Webb Space Telescope has enabled detailed spectroscopic characterization of sub-Neptune atmospheres. With detections of carbon- and oxygen-bearing species such as CO, CO$_2$, CH$_4$, and H$_2$O, a central question is whether the atmospheric C/O ratio, commonly used to trace formation location in giant planets, can serve a similar diagnostic role for sub-Neptunes. We use the global chemical equilibrium framework of Schlichting & Young (2022) to quantify how magma ocean-atmosphere interactions affect the atmospheric C/O ratio. We find that the resulting C/O ratios range from several orders of magnitude below solar to a few times solar. The atmospheric C/O ratio in sub-Neptunes is therefore not inherited from the protoplanetary disk, but instead emerges from chemical equilibrium between the atmosphere and the underlying magma ocean. Planetary mass, atmospheric mass fraction, and thermal state all strongly influence the atmospheric C/O ratio. In addition, carbon partitioning into the metal phase typically reduces the atmospheric C/O ratio substantially, particularly for atmospheric mass fractions less than a few percent. Finally, we couple the deep equilibrium compositions to 1D atmospheric models that self-consistently solve for the pressure-temperature structure and chemical composition, including photochemistry. We find that the C/O ratio varies with altitude under low vertical mixing conditions (K$_\text{zz}=10^4$ cm$^2$s$^{-1}$), but remains constant under strong mixing (K$_\text{zz}=10^7$ cm$^2$s$^{-1}$). Our results imply that observed C/O ratios of sub-Neptunes can be used to probe their interiors. Specifically, C/O ratios much lower than host star values would imply an underlying magma ocean with iron metal having sequestered significant amounts of carbon.

This paper describes the identification of Gaia counterparts of a sample of oxygen-rich AGB stars with OH maser emission as a first step towards the compilation of a general Gaia Catalogue of Galactic AGB stars. With this catalogue, tests of evolutionary models for the AGB star population in the solar neighbourhood will become feasible. We cross-matched AGB star candidates showing OH maser emission with Gaia DR3 using a cross-match with AllWISE and 2MASS as intermediate steps to avoid ambiguities. With the help of the Virtual Observatory, we retrieved photometric data from the near-ultraviolet to the far-infrared and built spectral energy distributions (SEDs) of the sources. The SEDs were fitted with theoretical models. The fit results, together with information from the literature, allowed us to clean the sample from non-AGB stars and to obtain bolometric fluxes for the AGB stars. Distances based on Gaia parallaxes were used to derive the stellar luminosities. We identified unique Gaia counterparts for 1487 OH masers. Of these, 1172 had an unambiguous classification as AGB stars; they make up the Gaia OH/IR star sample. Parallaxes with relative errors $<$20\% and astrometric excess noise $<$1.5 mas were available for 222 OH/IR stars. We conclude that the study of the AGB population in the solar neighbourhood is limited by the obscuration by circumstellar dust, as \G\ DR3 only provides parallaxes for a few of our candidates. The location of the OH/IR stars matches that of LPV discovered by \G\ in the (BP$-$RP; G$_{abs}$) diagram, but the OH/IR star sample is biased towards redder colours (BP$-$RP$>$4) mag and larger amplitudes ($>$1 mag in the G-band), which are typical for periodic large-amplitude Mira variables.

The tomographic Alcock-Paczynski (AP) method, developed over the past decade, exploits redshift evolution for cosmological determination, aiming to mitigate contamination from redshift distortions and capture nonlinear scale information. Marked Correlation Functions (MCFs) extend information beyond the two-point correlation. For the first time, this study integrated the tomographic AP test with MCFs to constrain the flat $w$CDM cosmology model. Our findings show that multiple density weights in MCFs outperform the traditional two-point correlation function, reducing the uncertainties of the matter density parameter $\Omega_m$ and dark energy equation of state $w$ by 48\% and 45\%, respectively. Furthermore, we introduce a novel principal component analysis (PCA) compression scheme that efficiently projects high-dimensional statistical measurements into a compact set of eigenmodes while preserving most of the cosmological information. This approach retains significantly more information than traditional coarse binning methods, which simply average adjacent bins in a lossy manner. Applying PCA compression also enables the effective use of marked correlation functions in 2D $(s,\mu)$ space, yielding an additional $\sim 50\%$ reduction in error margins. To assess robustness, we incorporate realistic redshift errors expected in future spectroscopic surveys. While these errors modestly degrade cosmological constraints, our combined framework, which utiizes MCFs and PCA compression within tomographic AP tests, is less affected and always yield to tight cosmological constraints. This scheme remains highly promising for upcoming slitless spectroscopic surveys, such as the Chinese Space Station Telescope (CSST).

The distribution of metals and temperature in the intracluster medium (ICM) provides key insights into galaxy cluster evolution, revealing information about chemical enrichment and heating and cooling processes, respectively. To access this information, it is crucial to understand the transport processes in the ICM. Here, we systematically study the transport mechanisms in the ICM with tracer particle resimulations of the Omega500 cosmological hydrosimulation, using a sample of four galaxy clusters of comparable masses but different mass assembly histories. Through the analysis of particle pair dispersion statistics, we find a time-dependent scaling index linked to the cluster's dynamical state. It reaches or exceeds Richardson scaling briefly during major mergers but remains much lower in relaxed clusters. We identify a coherent transport mode during major mergers that causes directional flow in the ICM. Although coherent transport can move particles to outer regions, the particles transported to the cluster outskirts compose only a small fraction of the density there; thus the anisotropy it creates in the overall density distribution is limited. Moreover, strong turbulence generated by mergers quickly disperses these particles, further limiting this effect. We also provide useful statistics on the radial evolution of the ICM and the fraction of particles that ever reached the inner regions as a function of radius. Our results show that major mergers primarily drive particle transport, linking ICM transport to merger-driven dynamics, and highlighting the interplay between coherent and turbulent transport.

The millimetre spectrum of Galactic regions and galaxies is rich in continuum and molecular lines. That diversity is mostly tackled by either broad-band photometry or high resolution heterodyne spectroscopy. We aim to map the millimetre continuum emission of Galactic regions with a spectral resolution intermediate between broad-band photometry and heterodyne spectroscopy, enabling us to quickly cover large sky areas with spectroscopy. We report observations of the Orion Nebula with the CONCERTO instrument, which was installed at the APEX telescope focal plane from April 2021 to December 2023. We find that the spectrum of Orion is dominated by dust emission with an emissivity index from 1.3 to 2.0 along with the strong CO(2-1) and H$_2$O lines which are naturally separated from the continuum, thanks to CONCERTO spectral capabilities. Many regions show a strong free-free emission as well, at lower frequencies. We demonstrate the spectral capabilities of CONCERTO at intermediate spectral resolution, with a frequency coverage from 130 to 310 GHz. A sensitivity of 200 mK is achieved in one second, for one beam and a 6 GHz frequency width, over an 18 arcmin diameter field-of-view, which is within a factor 3 of the expectations. We show that we can spectrally disentangle the continuum from the CO line emission. The slope of the millimetre continuum is line-free mapped for the first time in Orion.

Pulsar Wind nebulae (PWNe), are among the most efficient particle accelerators in the Universe, however understanding the physical conditions and the magnetic geometry in their inner region has always proved elusive. X-ray polarization provides now a unique opportunity to investigate the magnetic field structure and turbulence properties close to where high energy particles are accelerated. Here we report on the recent X-ray polarization measurement of the PWN 3C 58 by the International X-ray Polarimeter Explorer (IXPE). 3C 58 is a young system displaying a characteristic jet-torus structure which, unlike other PWNe, is seen almost edge on. This nebula shows a high level of integrated polarization ~ 22% at an angle ~ 97deg, with an implied magnetic field oriented parallel to the major axis of the inner torus, suggesting a toroidal magnetic geometry with little turbulence in the interior, and an intrinsic level of polarization possibly approaching the theoretical limit for synchrotron emission. No significant detection of a polarized signal from the associated pulsar was found. These results confirm that the internal structure of young PWNe is far less turbulent than previously predicted, and at odds with multidimensional numerical simulations.

Brightest cluster galaxies (BCG) are today passive and very massive galaxies at the center of their clusters, still accreting mass through swallowing companions, and flows of cold gas, regulated by radio-mode active galactic nucleus (AGN) feedback. However, their formation history is still a matter of debate. We report new results based on millimeter observations performed with the Northern Extended Millimeter Array (NOEMA) interferometer, mapping the cold molecular gas (CO) that feeds the star formation of distant BCGs. We selected three among the strongest cool-core BCGs at intermediate redshifts ($z\simeq0.4$), namely RX 1532, MACS 1447, and CHIPS 1911. Previous unresolved millimeter observations and multi-wavelength analysis showed that they are among the most star forming (${\rm SFR}\simeq100~ M_\odot/{\rm yr}$) and gas rich ($M_{H_2}\simeq10^{11}~M_\odot$) BCGs at intermediate redshifts. The selected BCGs are thus caught in a phase of rapid mass assembly, which makes them ideal targets for high-resolution observations of their molecular gas. By combining NOEMA intensity and velocity maps with archival images from the Hubble Space Telescope, we detect in-situ star formation, filaments of accreting cold gas likely regulated by AGN feedback, disturbed morphology associated with tidal tails of molecular gas, as well as gas compression and tails originated from stripping of gas. While effective condensation of the intra-cluster medium is required to explain the large molecular gas reservoirs, the BCGs exhibit a broad variety of environment-driven mechanisms responsible for the processing of their cold gas: flows of cooling gas (RX 1532), ram pressure or sloshing of the intra-cluster medium (MACS 1447), and galactic tides (CHIPS 1911). This study thus sheds new insights on the physical mechanisms responsible for the mass assembly of galaxies hosting AGN at the center of clusters.

We present a comprehensive investigation into the properties of 1H 0323+342, a prominent jetted narrow-line Seyfert 1 galaxy. The primary objective is to understand the interplay between the relativistic jet, the hot corona, and the accretion disk around the supermassive black hole. This study spans the years 2006 to 2023, incorporating a rich dataset with 172 Swift observations, including the optical, UV, and X-ray bands, integrated with Fermi Large Area Telescope (LAT) observations. Spectral analysis was conducted on the X-ray observations using the XSPEC software, and the results were compared with optical, UV, and gamma-ray flux measurements and photon index values. Our key findings include the identification of three distinct zones in the X-ray photon index-flux plot, characterized by high flux and a hard photon index (zone 1), high flux and a soft photon index (zone 2), and low flux and a soft photon index (zone 3). Before 2017, 1H 0323+342 moved back and forth between zones 1 and 2; after that epoch, it transitioned to zones 2 and 3. Correspondingly, we observed a decreasing jet activity in the Fermi/LAT data and a reduction in the accretion rate in optical/UV data from Swift/UVOT. We interpret these observations in the framework of an intermittent jet scenario, driven by radiation-pressure instability in the accretion disk.

Models of astrophysical dust are key to understand several physical processes, from the role of dust grains as cooling agents in the ISM to their evolution in dense circumstellar disks, explaining the occurrence of planetary systems around many stars. Currently, most models aim at providing optical properties for dust grains in the diffuse ISM, and many do not account properly for complexity in composition and structure when dust is expected to evolve in dense astrophysical environments. Our purpose is to investigate, with a pilot sample of micron-size dust grains, the influence of dust structure, porosity, and composition when computing the optical properties of grown dust grains in the infrared and millimetre domains, where observations are widely used to constrain the dust properties. Starting from the small dust grains developed in the THEMIS 2.0 model, we use the Discrete Dipole Approximation to compute the optical properties of 1 um grains, varying the hypotheses made on their composition and structure. We look at the dust scattering, emission and extinction to isolate potential simplifications and unavoidable differences between grain structures. We note that both dust structure and porosity influence the dust properties in infrared and millimetre ranges, demonstrating that dust aggregates cannot be correctly approximated by compact or porous spheres. In particular, we show that the dust emissivity index in the millimetre can vary with fixed grain size. We shed light on the importance of taking the dust structure and porosity into account when interpreting observations in environments where dust grains may have evolved significantly. Efforts in carrying out physical models of grain growth, for instance, are required to establish realistic constraints on the structure of grown dust grains, and will be used in the future to build realistic dust models for the dense ISM.

T CrB is a symbiotic recurrent nova that last erupted in 1946. Given its recurrence timescale of approximately 80 years, the next outburst is eagerly anticipated by the astronomical community. In this work, we analyse the optical light curves of T CrB, comparing recent photometric evolution with historical data to evaluate potential predictive indicators of nova eruptions. Although the "super-active" phases preceding both the 1946 and anticipated eruptions are strikingly similar, the subsequent photometric behaviour differs. We find that the decline in brightness observed in 2023, interpreted by some as a "pre-eruption dip", deviates from the deep minimum recorded prior to the 1946 event and does not reliably predict the eruption timing. Recent photometric and spectroscopic observations indicate that the system is returning to a high-accretion state. Given this, an eruption may be imminent, even without distinct precursors. While the next eruption of T CrB will be a major scientific event, its expected peak brightness of $V \sim 2$ mag highlights the importance of setting realistic public expectations for what will be a visually modest, yet astrophysically very significant, celestial event.

Context. It is important to determine which Galactic cosmic-ray sources can accelerate particles to the knee of the cosmic ray spectrum at a few PeV, and in particular whether supernova remnants may contribute. Current models for particle acceleration in very young remnants assume the circumstellar material consists of smooth, freely expanding winds. There is strong evidence that some supernovae expand into much denser circumstellar material including dense shells ejected by eruptions shortly before explosion. Aims. We investigate the effects of dense circumstellar shells on particle acceleration in supernova shocks during the first few years post-explosion, to quantify whether such interaction supernovae may act as PeVatrons. Methods. We used the pion code to model the circumstellar medium around Luminous Blue Variables after having a brief episode with a mass-loss rate of up to dM/dt = 2Msol/yr. Consequently, we performed spherically symmetric 1-D simulations using our time-dependent acceleration-code RATPaC in which we simultaneously solve the transport equations for cosmic-rays, magnetic turbulence, and the hydrodynamical flow of the thermal plasma in the test-particle limit. Results. We find that the interaction with the circumstellar shells can significantly boost the maximum energy by enhancing particle escape during the onset of the shock-shell interaction followed by the reacceleration of the shock propagating into a medium with a pre-amplified field. Early interactions boost the maximum energy to a greater degree and interactions within the first 5 months after explosion can increase Emax to more then 1 PeV.

We present narrowband-filtered integral field unit (IFU) observations of the Hot Dust-Obscured Galaxy (Hot DOG) WISE J202615.27$+$071624.0 (hereafter W2026$+$0716) at redshift $z=2.570$ using Keck/OSIRIS. Our analysis reveals a multi-component ionized gas outflow structure in this heavily obscured AGN host galaxy. Multi-component Gaussian decomposition of the [O III] and H$\alpha$ emission lines uncovers extremely broad and asymmetric profiles, characteristic of AGN-driven outflows. Kinematic mapping shows spatially distinct structures: the [O III] and H$\alpha$ dominated components (with radii of $1.20 \pm 0.56$ kpc) are separated by a projected offset of $\sim 1.1$ kpc and exhibit divergent velocity regimes. The [O III] outflow reaches a velocity of 3210 $\pm$ 50 km s$^{-1}$, while the H$\alpha$ outflow component attains 2310 $\pm$ 840 km s$^{-1}$. Dynamical modeling supports a biconical outflow structure, with [O III] and H$\alpha$ emissions dominating separate cones and significant dust obscuration of the redshifted outflow. Their comparable momentum outflow rates and energy outflow rates suggest a potential physical connection in their driving mechanisms. Spectral energy distribution (SED) analysis reveals anomalous optical/UV excess, attributed to AGN photon scattering by dust or outflowing material, classifying W2026+0716 as a "Blue Hot DOG". The derived outflow timescale ($\sim10^{5}$ yr) aligns with the evolutionary phase of Blue Hot DOGs, suggesting AGN feedback operates persistently during this transitional stage.

We present a comprehensive analysis of the time-resolved spectral energy distributions (SEDs) of the blazar OJ 287 over a 15-year period (2008-2023), using multi-wavelength data. In the $\gamma$-ray band, multiple flaring episodes were observed, with the strongest flare reaching a peak flux of $(5.60\pm1.11)\times10^{-7}\:{\rm photons\:cm^{-2}\:s^{-1}}$ on MJD 55869.03 (04 November 2011). In the optical/UV band, the source was in an active state between MJD 57360 (04 December 2015) and 57960 (26 July 2017), during which the highest flux of $(1.07\pm0.02)\times10^{-10}\:{\rm erg\:cm^{-2}\:s^{-1}}$ was observed on MJD 57681.23 (20 October 2016). In the X-ray band, both the flux and spectral index exhibit variability. To investigate the origin of the broadband emission from OJ 287, we systematically modeled 739 quasi-simultaneous SEDs using a leptonic model that self-consistently accounts for particle injection and cooling. This analysis is possible thanks to the recent development of a surrogate neural-network-based model, trained on kinetic simulations. This innovative, time-resolved, neural network-based approach overcomes the limitations of traditional single-epoch SED modeling, enabling to explore the temporal evolution of key model parameters, such as the magnetic field strength, Doppler factor, and electron injection distribution, across different states of the source. We identified distinct emission states characterized by unique combinations of magnetic field $ B $, electron index $ p $, and Doppler boost $ \delta $, associated to different underlying mechanisms such as varying acceleration processes (e.g., shocks, turbulence) and magnetic confinement. The analysis provides insights into the jet physics processes, including particle acceleration mechanisms and dynamic changes in the jet structure.

The eruption of the filament/flux rope generates the coronal perturbations, which further form EUV waves. There are two types of EUV waves, including fast-mode magnetosonic waves and slow waves. In this paper, we first report an event showing the Quasiperiodic Slow-Propagating (QSP) EUV "wave" trains during an M6.4-class flare (SOL2023-02-25T18:40), using multiple observations from SDO/AIA, CHASE/HIS, ASO-S/FMG, SUTRI, and LASCO/C2. The QSP "wave" trains occurred as the filament showed a rapid rise. The QSP "wave" trains have the projected speeds of 50-130 km s$^{-1}$ on the plane of the sky, which is slower than the fast-mode magnetosonic speed in the solar corona. And the calculated period of the QSP wave trains is 117.9 s, which is in good agreement with the associated flare Quasi-Periodic Pulsation (140.3 s). The QSP wave trains could be observed during the entire impulsive phase of the flare and lasted about 30 minutes in the field of view (FOV) of SDO/AIA. About 30 minutes later, they appeared in the FOV of LASCO/C2 and propagated to the northwest. We suggest that the QSP wave trains are probably apparent waves that are caused by the successive stretching of the inclined field lines overlying the eruptive filament. The periodic pattern of the QSP wave trains may be related to the intermittent energy release during the flare.

Solar nanojets are small-scale jets generated by component magnetic reconnection, characterized by collimated plasma motion perpendicular to the reconnecting magnetic field lines. As an indicator of nanoflare events, they are believed to play a significant role in coronal heating. Using high-resolution extreme-ultraviolet (EUV) imaging observations from the Extreme Ultraviolet Imager (EUI) onboard the Solar Orbiter mission, we identified 27 nanojets in an erupting filament on September 30, 2024. They are potentially associated with the untwisting of magnetic field lines of the filament. Most nanojets exhibit velocities around 450 km s$^{-1}$, with the fastest reaching approximately 800 km s$^{-1}$, significantly higher than previously reported but comparable to the typical coronal Alfvén speed. To our knowledge, these are the highest speeds ever reported for small-scale jets (less than ~1 Mm wide) in the solar atmosphere. Our findings suggest that these nanoflare-type phenomena can be more dynamic than previously recognized and may contribute to the energy release process of solar eruptions and the heating of coronal active regions.

Recently, there has been considerable debate regarding potential evidence for the dynamical nature of dark energy (DE), particularly in light of Baryon Acoustic Oscillations (BAO) measurements released by DESI survey. In this work, we propose an agnostic test that simultaneously constrains the dark energy (DE) equation of state (EoS) and probes the possibility of a transition between the quintessence and phantom regimes, or vice versa. Our initial approach is independent of physical priors, allowing the data to determine which behavior best fits the parameters. We then consider a minimally modified gravity theory known as VCDM, into which we can map our initial approximation, placing it within a theoretically stable framework. To this end, we incorporate the most up-to-date datasets available, including BAO measurements from DESI-DR2, Type Ia Supernovae from the PantheonPlus, DESY5, and Union3 samples, as well as Cosmic Microwave Background (CMB) data from Planck. Our analysis reveals strong and statistically significant evidence for a quintessence-phantom transition across various data combinations. \textit{The strongest evidence is found for Planck+DESI+DESY5, with a significance exceeding $\sim$5$\sigma$ in favor of a quintessence-phantom transition at $z_† = 0.493^{+0.063}_{-0.081}$}. Beyond this redshift, the EoS remains within the phantom regime, while for $z < z_†$, it favors the quintessence regime. Despite this strong indication, \textit{we find that such transitions do not resolve the $H_0$ tension}.

Aims: We aim to study the near-infrared variability of the T8 dwarfs 2MASS J09393548-2448279 and EQ J1959-3338 by analyzing their $J$-band photometric signal, which can provide new insights into the atmospheric dynamics of cold brown dwarfs. Methods: We used FLAMINGOS-2 on the Gemini South telescope to perform $J$-band differential photometry continuously over 4 h for each target. The resulting light curves have a cadence of 20 s and a photometric uncertainty of 2-4 mmag. Results: We detect periodic variability in both T8 dwarfs, with an amplitude of $16.6\pm0.9$ mmag and a period of $1.364\pm0.012$ h for EQ J1959-3338, which we attribute to rotational modulation. For 2MASS J09393548-2448279, we observe an amplitude of $4.6\pm0.4$ mmag and a period of $1.733\pm0.040$ h, though this periodicity could represent a fraction of a longer period. Conclusions: With the detection of variability in 2MASS J09393548-2448279 and EQ J1959-3338, the number of known variable T8 dwarfs has doubled, making them prime candidates for infrared space-based monitoring and radio observations to investigate atmospheric dynamics and the influence of the magnetic field in very cool atmospheres.

The near-infrared YJKs Visual and Infrared Survey Telescope for Astronomy (VISTA) survey of the Magellanic Clouds (VMC) is complete and there are also data from additional programmes enhancing its quality over the original footprint. This work presents the final data release of the VMC survey, which includes additional observations and provides an overview of the scientific results. The overall data quality has been revised and reprocessed standard data products, that already appeared in previous data releases, are made available together with new data products. These include individual stellar proper motions, reddening towards red clump stars and source classifications. Several data products, such as the parameters of some variable stars and of background galaxies, from the VMC publications are associated to a data release for the first time. The data are processed using the VISTA Data Flow System and additional products, e.g. catalogues with point-spread function photometry or tables with stellar proper motions, are obtained with software developed by the survey team. This release supersedes all previous data releases of the VMC survey for the combined (deepstacked) data products, whilst providing additional (complementary) images and catalogues of single observations per filter. Overall, it includes about 64 million detections, split nearly evenly between sources with stellar or galaxy profiles. The VMC survey provides a homogeneous data set resulting from deep and multi-epoch YJKs-band imaging observations of the Large and Small Clouds, the Bridge and two fields in the Stream. The VMC data represent a valuable counterpart for sources detected at other wavelengths for both stars and background galaxies.

We present a timing analysis of observations with the Hard X-ray Modulation Telescope of the black hole X-ray transient Swift J1727.8-1613 during its 2023 outburst. We detect, for the first time in a black hole X-ray binary, a prominent dip at ~ 3-15 Hz in the real part of the cross spectrum between high-energy (>25 keV) and low-energy (<10 keV) photons in the Low Hard and Hard Intermediate States, during which the QPO frequency rapidly increases and then stabilizes at ~ 1.0-1.5 Hz. Remarkably, the real part of the cross spectrum reaches negative values at the frequencies around the minimum of the dip, indicative of a phase lag ranging between ${\pi}/2$ and ${\pi}$ in this frequency range. We fit the power spectra and the real and imaginary parts of the cross spectra simultaneously using a multi-Lorentzian model. Among the lag models, the Gaussian phase-lag model provides a slightly better reduced ${\chi}^2$ than the constant phase-lag and constant time-lag models, while it also alleviates the degeneracy associated with those models. From the parameters of the Lorentzian that fits the dip, we estimate the size of the accretion flow, which consistently exceeds 10,000 km as the QPO frequency increases from 0.13 Hz to 2.0 Hz. Furthermore, both the energy-dependent phase-lag and fractional-rms spectra of the dip exhibit a change in trend around 15 keV, with the phase lag dropping and rms reaching a local minimum. These spectra closely resemble the shapes predicted by the time-dependent Comptonization model, vKompth, for a low feedback factor, offering a pathway to explain the radiative properties of the corona. Additionally, the coherence function suggests a diversity of variability components, potentially arising from different parts of the corona.

We cross-correlated galaxies from the LOw-Frequency ARray (LOFAR) Two-metre Sky Survey (LoTSS) second data release (DR2) radio source with the extended Baryon Oscillation Spectroscopic Survey (eBOSS) luminous red galaxy (LRG) sample to extract the baryon acoustic oscillation (BAO) signal and constrain the linear clustering bias of radio sources in LoTSS DR2. In the LoTSS DR2 catalogue, employing a flux density limit of $1.5$ mJy at the central LoTSS frequency of 144 MHz and a signal-to-noise ratio (S/N) of $7.5$, additionally considering eBOSS LRGs with redshifts between 0.6 and 1, we measured both the angular LoTSS-eBOSS cross-power spectrum and the angular eBOSS auto-power spectrum. These measurements were performed across various eBOSS redshift tomographic bins with a width of $\Delta z=0.06$. By marginalising over the broadband shape of the angular power spectra, we searched for a BAO signal in cross-correlation with radio galaxies, and determine the linear clustering bias of LoTSS radio sources for a constant-bias and an evolving-bias model. Using the cross-correlation, we measured the isotropic BAO dilation parameter as $\alpha=1.01\pm 0.11$ at $z_{\rm eff}=0.63$. By combining four redshift slices at $z_{\rm eff}=0.63, 0.69, 0.75$, and $0.81$, we determined a more constrained value of $\alpha = 0.968^{+0.060}_{-0.095}$. For the entire redshift range of $z_{\rm eff}=0.715$, we measured $b_C = 2.64 \pm 0.20$ for the constant-bias model, $b(z)=b_C$, and then $b_D = 1.80 \pm 0.13$ for the evolving-bias model, $b(z) = b_D / D(z)$, with $D(z)$ denoting the growth rate of linear structures. Additionally, we measured the clustering bias for individual redshift bins.

We investigate the statistical distribution of source counts-in-cells in the second data release of the LOFAR Two-Metre Sky Survey (LoTSS-DR2) and we test a computationally cheap method based on the counts-in-cells to estimate the two-point correlation function. We compare three stochastic models for the counts-in-cells which result in a Poisson distribution, a compound Poisson distribution, and a negative binomial distribution. By analysing the variance of counts-in-cells for various cell sizes, we fit the reduced normalised variance to a single power-law model representing the angular two-point correlation function. Our analysis confirms that radio sources are not Poisson distributed, which is most likely due to multiple physical components of radio sources. Employing instead a Cox process, we show that there is strong evidence in favour of the negative binomial distribution above a flux density threshold of 2 mJy. Additionally, the mean number of radio components derived from the negative binomial distribution is in good agreement with corresponding estimates based on the value-added catalogue of LoTSS-DR2. The scaling of the counts-in-cells normalised variance with cell size is in good agreement with a power-law model for the angular two-point correlation. At a flux density threshold of 2 mJy and a signal-to-noise ratio of 7.5 for individual radio sources, we find that for a range of angular scales large enough to not be affected by the multi-component nature of radio sources, the value of the exponent of the power law ranges from -0.8 to -1.05. This closely aligns with findings from previous optical, infrared, and radio surveys of the large scale structure. The scaling of the counts-in-cells statistics with cell size provides a computationally efficient method to estimate the two-point correlation properties, offering a valuable tool for future large-scale structure studies.

Correlations between various multiwavelength (MWL) bands are an intermittent feature in blazar light curves; that is, they are observed in some instances but not in others. With the CTAO we will obtain detailed very-high-energy (VHE) gamma-ray light curves for many sources also during their low states, enabling detailed MWL correlation studies. For two blazars, the HBL Mrk\,421 and the FSRQ PKS\,1510-089, the long-term X-ray and optical light curves are used to induce variations in input parameters of the lepto-hadronic one-zone code OneHaLe. We show light curves in the CTA energy range for three different energy thresholds. The results are: 1) the presence of relativistic protons has a significant effect on the correlation of the light curves as the emerging pair cascade prolongs flaring states at the highest energies; and 2) comparison of the theoretical light curves with existing VHE gamma-ray data shows that both leptonic and hadronic models can only partially reproduce the data.

A new method for determining the nature of the orbital motion (chaotic or regular) of globular clusters in the central region of the Galaxy with a radius of 3.5 kpc, which are most affected by the bar, is proposed. The method is based on calculating the orbital power spectrum as a function of time and calculating the entropy of the power spectrum as a measure of orbital chaos. The sample includes 45 globular clusters. To form the 6D phase space required for integrating the orbits, the most accurate astrometric data to date from the Gaia satellite (Vasiliev \& Baumgardt, 2021) were used, as well as new refined average distances (Baumgardt \& Vasiliev, 2021). Orbits of globular clusters are obtained in a non-axisymmetric potential with a bar in the form of a triaxial ellipsoid embedded in an axisymmetric potential, traditionally used by us to construct orbits of globular clusters, described in detail in the paper by Bajkova et al. (2023a). The following, most realistic, bar parameters are adopted: mass $10^{10} M_\odot$, length of the major semi-axis 5 kpc, angle of rotation of the bar axis 25$^o$, angular velocity of rotation 40 km s$^{-1}$ kpc$^{-1}$. A list of 23 globular clusters with regular dynamics and 22 globular clusters with chaotic dynamics is determined. The correlation of the obtained classification of globular clusters with the classification obtained by us using other methods in the work of Bajkova et al. (2024a) was determined.

Recent measurements from the Atacama Cosmology Telescope (ACT), combined with Planck and DESI data, suggest a higher value for the spectral index $n_s$. This places Starobinsky inflation at the edge of the $2\sigma$ constraints for a number of e-folds $N_\star$ around $60$ when using the usual analytical approximations. We present refined predictions for Starobinsky inflation that go beyond the commonly used analytical approximations. By evaluating the model with these improved expressions, we show that for $N_\star \gtrsim 60$ it remains consistent with current observational constraints at the $2\sigma$ level. Additionally, we examine the implications of the ACT results for post-inflationary reheating parameters. Specifically, we find a lower bound on the effective equation of state parameter during reheating of approximately $\omega \gtrsim 0.462$; this excludes purely perturbative reheating, which leads to $\omega \simeq 0$. We also show that the reheating temperature is constrained to be $T_{\text{rh}} \lesssim 2 \times 10^{12}~\text{GeV}$, assuming $\omega \leq 1$. Furthermore, we find that the predictions for the spectral index and tensor-to-scalar ratio can lie within $1\sigma$ of the recent ACT constraints if the reheating temperature satisfies $4~\text{MeV} \lesssim T_{\text{rh}} \lesssim 10~\text{GeV}$ for $0.8 \lesssim \omega \leq 1$.

Spectroscopy for faint dwarf galaxies outside of our own Local Group is challenging. Here, we present MUSE spectroscopy to study the properties of four known dwarf satellites and one stellar stream (KK208) surrounding the nearby grand spiral M83, which resides together with the lenticular galaxy Cen A in the Centaurus group. This data complete the phase-space information for all known dwarf galaxies around M83 down to a completeness of $-$10 mag in the $V$ band. All studied objects have an intermediate to old and metal-poor stellar population and follow the stellar luminosity-metallicity relation as defined by the Local Group dwarfs. For the stellar stream we serendipitously identify a previously unknown globular cluster, which is old and metal-poor. Two dwarf galaxies (NGC5264 and dw1341-29) may be a bound satellite of a satellite system due to their proximity and shared velocities. Having access to the positions and velocities of 13 dwarfs around M83, we estimate the mass of the group with different estimators. Ranging between 1.3 and $3.0 \times 10^{12}$ M$_\odot$ for the halo mass we find it to be larger than previously assumed. This may impact the previously reported tension for cold dark matter cosmology with the count of dwarf galaxies. In contrast to Cen A, we do not find a co-rotating plane-of-satellites around M83.

WZ Sge-type dwarf novae (DNe) form a subclass in cataclysmic variables, characterized by short-period variations called superhumps during an outburst. Here we present optical ground-based and TESS observations of ASASSN-24hd in its 2024-2025 outburst. ASASSN-24hd is the first reported WZ Sge-type DN outburst fully covered by TESS, providing a great opportunity to study the evolution of superhumps. Our observations establish its early and stage-A ordinary superhumps as 0.05711(4) and 0.05919(5) d, respectively, resulting in its mass ratio of 0.098(4). The TESS observations confirm that the evolution of its superhump period, amplitude, and profile after the appearance of ordinary superhumps is generally consistent with those of SU UMa-type DNe observed with Kepler and TESS. Furthermore, we find that ASASSN-24hd in outburst shares a great similarity to the 2010 superoutburst of an SU UMa-type DN V585 Lyr, observed by Kepler, particularly including the superhump evolution and the long waiting time ($\gtrsimeq$ 5 d) before the stage A--B transition of ordinary superhumps. The shorter superoutburst cycles and smaller outburst amplitude in V585 Lyr than those of ASASSN-24hd disfavor the interpretation that V585 Lyr is, in fact, a face-on WZ Sge-type DN where early superhumps are undetectable. Instead, one possibility of their critical differences is either low quiescence viscosity or inner disk truncation, which has been invoked to explain the extreme nature of WZ Sge-type DNe, but future observations in quiescence are vital to conclude. These findings emphasize the borderline between SU UMa-type and WZ Sge-type DNe.

Binary neutron star mergers provide insight into strong-field gravity and the properties of ultra-dense nuclear matter. These events offer the potential to search for signatures of physics beyond the standard model, including dark matter. We present the first numerical-relativity simulations of binary neutron star mergers admixed with dark matter, based on constraint-solved initial data. Modeling dark matter as a non-interacting fermionic gas, we investigate the impact of varying dark matter fractions and particle masses on the merger dynamics, ejecta mass, post-merger remnant properties, and the emitted gravitational waves. Our simulations suggest that the dark matter morphology - a dense core or a diluted halo - may alter the merger outcome. Scenarios with a dark matter core tend to exhibit a higher probability of prompt collapse, while those with a dark matter halo develop a common envelope, embedding the whole binary. Furthermore, gravitational wave signals from mergers with dark matter halo configurations exhibit significant deviations from standard models when the tidal deformability is calculated in a two-fluid framework. This highlights the need for refined models in calculating the tidal deformability when considering mergers with extended dark matter structures. These initial results provide a basis for further exploration of dark matter's role in binary neutron star mergers and their associated gravitational wave emission and can serve as a benchmark for future observations from advanced detectors and multi-messenger astrophysics.

Shocks associated with Interplanetary Coronal Mass Ejections are known to energize charged particles and give rise to Solar Energetic Particles. Many of these energetic particles move ahead of the shock to create a foreshock region. The foreshock region primarily consists of solar wind plasma, exhibiting turbulent velocity and magnetic fields. Such turbulent behavior results from inherent solar wind turbulence modified by energetic particles. We analyze magnetic field data from six such ICME shocks observed by the Wind spacecraft. The analysis of the shock upstream shows that the magnetic power spectral density (PSD) maintains a power-law slope of $-5/3$. We also identify clear intermittent peaks in the PSD. After characterizing these peaks, we investigate various possibilities for their generation. Our analysis indicates that these peaks in the PSD are due to the resonant interaction of Alfvén waves with the bulk solar wind protons and protons with energy up to $10$~keV. However, evidence of Alfvén wave interaction with highly energetic protons is not evident in our analysis, and we anticipate that such evidence is obscured by the prevailing solar wind turbulence in the shock upstream.

The radial acceleration relation (RAR) is a fundamental relation linking baryonic and dark matter in galaxies by relating the observed acceleration derived from dynamics to the one estimated from the baryonic mass. This relation exhibits small scatter, thus providing key constraints for models of galaxy formation and evolution -- allowing us to map the distribution of dark matter in galaxies -- as well as models of modified dynamics. However, it has only been extensively studied in the very local Universe with largely heterogeneous samples. We present a new measurement of the RAR, utilising a homogeneous, unbiased sample of H{\sc i}-selected galaxies. We introduce a novel approach of measuring resolved stellar masses using spectral energy distribution (SED) fitting across 10 photometric bands to determine the resolved mass-to-light ratio, which we show is essential for measuring the acceleration due to baryons in the low-acceleration regime. We find that the stellar mass-to-light ratio varies across galaxies and radially, favouring lower mass-to-light ratios compared to previous studies. Our results reveal a tight RAR with a low-acceleration power-law slope of $\sim 0.5$, consistent with previous studies. Adopting a spatially varying mass-to-light ratio yields the tightest RAR with an intrinsic scatter of only $0.045 \pm 0.022$ dex. We also find the first tentative evidence for redshift evolution in the acceleration scale, but more data will be required to confirm this. These results highlight the importance of resolved stellar mass measurements in accurately characterizing the gravitational contribution of the baryons in low-mass, gas-rich galaxies.

After decades of uncertainty about the origin of lithium, recent evidence suggests Galactic novae as its main astrophysical source. In this work, we present possible evidence for the first detection of the $^7$Be line at 478 keV, observed with the INTEGRAL satellite. The emission is temporally and spatially coincident with the outburst of the bright nova V1369 Cen, and line significance ranges from 2.5$\sigma$ to $\sim$1.9$\sigma$, depending on the detection methodology. A bootstrap analysis, assuming a fixed FWHM of 8 keV, provides a flux of $(4.9 \pm 2.0) \times 10^{-4}$ ph/cm$^2$/s centered at 479.0 $\pm$ 2.5 keV, with a 2.5$\sigma$ significant excess. This flux implies a total $^7$Be mass of $M_{^7Be} = (1.2^{+2.0}_{-0.6})$ $\times 10^{-8}$ M$_{\odot}$ at the distance determined using several indicators including the {\em Gaia} satellite. For a nova ejected mass estimated from radio observations, this result implies a $^7$Be=Li yield corresponding to $A(Li) = 7.1^{+0.7}_{-0.3}$. This value is comparable to those measured in a dozen novae through optical observations. Crucially, we confirm optically derived $^7$Li yields and demonstrate the groundbreaking potential of using gamma-ray data to measure Li abundances.

In our Bayesian and Statistical Analysis investigation of the open cluster NGC 6416, we utilized Gaia EDR3 astrometry data and ensemble-based unsupervised machine learning techniques to identify 406 cluster members. Using the MESA Isochrones and Stellar Tracks (MIST) with Gaia EDR3 data, we determined the following parameters for NGC 6416: a distance of approximately 1021pc, an age of about $12.58\pm 0.1$ Myr, a metallicity (z) of roughly $0.032\pm0.0015$, a binarity fraction near $0.419\pm0.021$, and an extinction ($A_V$) of approximately $0.995\pm0.058$ mag for an $R_V$ value of around $3.064\pm0.102$. We also fitted the radial surface density profile and conducted orbit analysis of the cluster using galpy. And finally found that the star formation scenarios are not observed in the open star cluster NGC 6416.

We present numerical simulations of discrete relativistic ejecta from an X-ray binary (XRB) with initial conditions directly informed by observations. XRBs have been observed to launch powerful discrete plasma ejecta during state transitions, which can propagate up to parsec distances. Understanding these ejection events unveils new understanding of jet-launching, jet power, and jet-ISM interaction among other implications. Multi-frequency quasi-simultaneous radio observations of ejecta from the black hole XRB MAXI J1820+070 produced both size and calorimetry constraints, which we use as initial conditions of a relativistic hydrodynamic simulation. We qualitatively reproduce the observed deceleration of the ejecta in a homogeneous interstellar medium (ISM). Our simulations demonstrate that the ejecta must be denser than the ISM, the ISM be significantly low-density, and the launch be extremely powerful, in order to propagate to the observed distances. The blob propagates and clears out a high-pressure low-density cavity in its wake, providing an explanation for this pre-existing low-density environment, as well as 'bubble-like' environments in the vicinity of XRBs inferred from other studies. As the blob decelerates, we observe the onset of instabilities and a long-lived reverse shock -- these mechanisms convert kinetic to internal energy in the blob, responsible for in-situ particle acceleration. We transform the outputs of our simulation into pseudo-radio images, incorporating the u,v coverage of the MeerKAT and e-MERLIN telescopes from the original observations with real-sky background. Through this, we maximize the interpretability of the results and provide direct comparison to current data, as well as provide prediction capabilities.

Un-doing the effect of gravitational lensing on the Cosmic Microwave Background ('de-lensing') is essential in shaping constraints on weak signals limited by lensing effects on the CMB, for example on a background of primordial gravitational waves. Removing these anisotropies induced by large-scale structures from the CMB maps also generally helps our view of the primordial Universe, by sharpening the acoustic peaks and the damping tail. However, delensing does transfer parts of these anisotropies to the noise maps. This will induce a new large scale 'mean-field' bias to any anisotropy estimator applied to the delensed CMB, and this bias directly traces large-scale structures. This paper analytically quantifies this delensed noise mean-field and its impact on quadratic (QE) and likelihood-based lensing estimators. We show that for Simons-Observatory-like surveys, this mean-field bias can reach 15% in cross-correlation with large-scale structures if unaccounted for. We further demonstrate that this delensed noise mean-field can be safely neglected in likelihood-based estimators without compromising the quality of lensing reconstruction or $B$-mode delensing, provided the resulting lensing map is properly renormalized.

The XRISM Resolve microcalorimeter array measured the velocities of hot intracluster gas at two positions in the Coma galaxy cluster: 3'x3' squares at the center and at 6' (170 kpc) to the south. We find the line-of-sight velocity dispersions in those regions to be sigma_z=208+-12 km/s and 202+-24 km/s, respectively. The central value corresponds to a 3D Mach number of M=0.24+-0.015 and the ratio of the kinetic pressure of small-scale motions to thermal pressure in the intracluster plasma of only 3.1+-0.4%, at the lower end of predictions from cosmological simulations for merging clusters like Coma, and similar to that observed in the cool core of the relaxed cluster A2029. Meanwhile, the gas in both regions exhibits high line-of-sight velocity differences from the mean velocity of the cluster galaxies, Delta v_z=450+-15 km/s and 730+-30 km/s, respectively. A small contribution from an additional gas velocity component, consistent with the cluster optical mean, is detected along a sightline near the cluster center. The combination of the observed velocity dispersions and bulk velocities is not described by a Kolmogorov velocity power spectrum of steady-state turbulence; instead, the data imply a much steeper effective slope (i.e., relatively more power at larger linear scales). This may indicate either a very large dissipation scale resulting in the suppression of small-scale motions, or a transient dynamic state of the cluster, where large-scale gas flows generated by an ongoing merger have not yet cascaded down to small scales.

Between the 4th and 6th of September 2024, the Astronomy & Astrophysics group at the University of Warwick held a meeting to celebrate 21 years of astronomy at Warwick and the scientific legacy of the late Prof. Tom Marsh, the group founder. More than a hundred people attended the meeting, with about half of the attendees being external delegates and coming from as far afield as the USA and South Africa. Tom Marsh moved to the University of Warwick from Southampton in 2003, after the Department of Physics decided to expand the scope of its research. From its humble beginnings with only two staff members, Tom himself and Boris Gänsicke, one postdoc and a couple of PhD students, the group has now grown to more than 95 members, including 25 staff. Tom pioneered the development of Doppler tomography, led key discoveries in the field of double-degenerate binary systems and made extensive contributions to instrumentation, primarily to developing the high-speed imaging photometers ULTRACAM, ULTRASPEC and HiPERCAM. This article provides a summary of Tom's legacy and Warwick's history as presented in the 21 years of Astronomy at Warwick meeting.

Flat radio spectra of compact jets launched by both supermassive and stellar-mass black holes are explained by an interplay of self-absorbed synchrotron emission up to some distance along the jet and optically thin synchrotron at larger distances (Blandford & Konigl 1979). Their spatial structure is usually studied using core shifts, in which the position of the peak (core) of the emission depends on the frequency. Here, we propose a novel method to fit the spatial dependence of the flux density at a given frequency of the jet and counterjet (when observed) using the theoretical spatial dependencies, which we provide as simple analytical formulae. We apply our method to the spatial structure of the jets in the luminous hard spectral state of the black-hole X-ray binary Swift J1727.8--1613. It was the most resolved continuous jet from an X-ray binary ever observed. We find that the observed approaching jet is significantly intrinsically stronger than the receding one, which we attribute to an increase in the emission of both jets with time (observationally confirmed), together with the light travel effect, causing the receding jet to be observed at an earlier epoch than the approaching one. The jets are relatively slow, with the velocity $\sim(0.3$--$0.4)c$. Our findings imply that the magnetic field strength increased with time. Also, the magnetic flux is much lower than in jets launched by `Magnetically Arrested Disks'. Our method is general, and we propose that it be applied to jets launched by stellar-mass and supermassive black holes.

Ultra-short-period (USP) planets, with orbital periods shorter than one day, represent a unique class of exoplanets whose origin remains puzzling. Determining their age distribution and temporal evolution is vital for uncovering their formation and evolutionary pathways. Using a sample of over 1,000 short-period planets around Sun-like stars, we find that the host stars of USP planets are relatively older and have a higher prevalence in the Galactic thick disk compared to stars hosting other short-period planets. Furthermore, we find that the occurrence of USP planets increases with stellar age and uncover evidence indicating that USP planetary system architectures evolve on Gyr timescales. This includes a distinct dip-pileup in period distributions around ~1 day and an expansion of orbital spacings with time. In addition, younger USP planet systems are observed to have fewer multiple transiting planets, implying fewer nearby companions and/or larger mutual orbital inclinations. Our findings suggest that USP planets continuously form through inward migration driven by tidal dissipation over Gyr timescales, and that younger and older USP planets may have originated via different specific tidal migration pathways.

We analyse the large-scale angular clustering of quasars in the Gaia-unWISE quasar catalog, Quaia, and their cross-correlation with maps of the lensing convergence of the Cosmic Microwave Background (CMB), to constrain the level of primordial non-Gaussianity (PNG). Specifically, we target the scale-dependent bias that would be induced by PNG on biased tracers of the matter inhomogeneities on large scales. The Quaia sample is particularly well suited for this analysis, given the large effective volume covered, and our ability to map out the main potential sources of systematic contamination and mitigate their impact. Using the universality relation to characterise the response of the quasar overdensity to PNG ($p_\phi=1$), we report constraints on the local-type PNG parameter $f_{\rm NL}$ of $f_{\rm NL}=-20.5^{+19.0}_{-18.1}$ (68\% C.L.) by combining the quasar auto-correlation and its cross-correlation with CMB lensing in two tomographic redshift bins (or $f_{\rm NL}=-28.7^{+26.1}_{-24.6}$ if assuming a lower response for quasars, $p_\phi=1.6$). Using the CMB lensing cross-correlations alone, we find $f_{\rm NL}=-13.8^{+26.7}_{-25.0}$. These are the tightest constraints on $f_{\rm NL}$ to date from angular clustering statistics and cross-correlations with CMB lensing.

The Stochastic Gravitational Wave Background (SGWB) from cosmic superstrings offers one of the few known possibilities to test String Theory within current experimental reach. However, in order to be compatible with the existing constraints, the tension of a cosmic superstring network is required to lie several orders of magnitude below the Planck scale. This is naturally realized in string compactifications where the volume of the extra dimensions is parametrically large (in string units). We estimate the GW spectrum arising from cosmic superstrings in such scenarios, providing analytical formulae as well as numerical results. Crucially, we do so within a fully-fledged string cosmology, taking into account various modified cosmological epochs (such as kination or early matter domination) induced by the presence of moduli and a time-dependent string tension. We show that part of the spectrum generically lies within reach of LISA and ET, with a large class of models predicting a characteristic drop in the amplitude which may be robustly probed by LISA. The corresponding signal would encode information on the dynamics of moduli and reheating. On the other hand, the ultra-high frequency part of the spectrum can be significantly enhanced by a long, early phase of kination with time-varying tension, yielding a spectral index unique to this set-up.

The formation and evolution of ultra-short-period (USP) planets is poorly understood. However, it is widely thought that these planets migrated inwards through interactions with outer neighbours. We aimed to confirm and characterize the USP Earth-sized planet K2-157 b ($P_{\rm orb}$ = 8.8 h). To do so, We measured 49 radial velocities (RVs) with the ESPRESSO spectrograph and derived the properties of the system through a RV and transit model. We detect two additional super-Neptune-mass planets within the warm Neptunian savanna, K2-157 c ($P_{\rm orb, c}$ = $25.942^{+0.045}_{-0.044}$ d, $M_{\rm p, c} \, \textrm{sin} \, i$ = $30.8 \pm 1.9$ $\rm M_{\oplus}$), and K2-157 d ($P_{\rm orb, d}$ = $66.50^{+0.71}_{-0.59}$ d, $M_{\rm p,d} \, \textrm{sin}\,i$ = $23.3 \pm 2.5$ $\rm M_{\oplus}$). The mass of K2-157 b, $M_{\rm p,b}$ = $1.14^{+0.39}_{-0.41}$ $\rm M_{\oplus}$ ($<$ 2.3 $\rm M_{\oplus}$ at 3$\sigma$), together with its radius, $R_{\rm p}$ = 0.935 $\pm$ 0.090 $\rm R_{\oplus}$, make the planet compatible with a rocky composition. K2 data discard non-grazing transit configurations for K2-157 c ($i_{\rm c}$ $<$ 88.4$^{\circ}$ at 3$\sigma$), and ESPRESSO data constrain the eccentricities of K2-157 c and K2-157 d to $e_{\rm c}$ $<$ 0.2 and $e_{\rm d}$ $<$ 0.5 at 3$\sigma$. At a population level, we find that the trend that the closest USP planets tend to orbit late-type stars does not hold when scaling the orbital separation to the Roche limit, which suggests that the orbital distribution of the closest planets across spectral types is determined by tidal disruption. The orbital architecture of K2-157 is unusual, with only one similar case reported to date: 55 Cnc. The USP planets of these systems, being accompanied by massive, long-period, relatively spaced, and possibly misaligned neighbours, could have migrated inwards through eccentricity-based mechanisms triggered by secular interactions.

The formation of solitons (such as closed domain walls) in the super-Early Universe is predicted in a number of theories of the formation of primordial black holes. However, the interaction of particles of the surrounding medium with the solitons should affect their dynamics. In the paper, we consider the interaction between domain walls and scalar particles which can play a role of dark matter. It is shown that when the temperature of the scalar particle gas, caused by the expansion of the Universe, decreases below a certain threshold value, the wall abruptly becomes opaque and locks particles inside itself. We discuss the dynamics of a single domain wall taking into account pressure of scalar particles locked inside a closed wall. It is shown, this effect leads to a time delay of domain wall collapse and the deferred formation of primordial black holes.

We demonstrate that induced gravitational waves (IGWs) can naturally emerge within the frame- work of thermal leptogenesis models, thereby providing a robust probe for exploring this theory at remarkably high energy scales. To illustrate this principle, we put forth a basic leptogenesis model in which an early matter-dominated phase, tracing the leptogenesis scale, enhances the generation of gravitational waves induced by an early structure formation. Leveraging recent N-body and lat- tice simulation results for IGW computations in the non-linear regime, we show that it is possible to establish a direct link between the frequency and amplitude of these IGWs and the thermal leptogenesis scale.

We report the existence of two exotic compact objects in the leading relativistic model of modified Newtonian dynamics, namely aether-scalar-tensor theory. This model is consistent with precision cosmology and gravitational wave constraints on tensor speed. Black hole mimickers could subtly change observations: gravitational waves from their mergers might show unusual echoes or altered ringdown patterns, and images of their horizon-scale shadows might be slightly different from those of a true black hole. Shapiro-free lenses are massless objects that deflect light without any gravitational time delay, producing distinctive lensing events. These predictions connect to ongoing and future gravitational-wave searches, horizon-scale imaging, and time-domain lensing surveys.

Neutrino self-interactions ({\nu}SI) mediated by light bosonic particles can produce characteristic spectral dips in astrophysical neutrino fluxes, thereby modifying the expected energy spectrum. The high-energy astrophysical neutrino spectrum has been extensively used to probe {\nu}SI models through these distinctive features. The recent detection of the ultrahigh-energy event KM3-230213A offers a new opportunity to explore {\nu}SI phenomenology at extreme energies. In this work, we investigate two implications of this observation under the assumption that the event originates from a diffuse power-law spectrum. First, we find that {\nu}SI induces spectral distortions that can mildly alleviate tensions between the KM3-230213A detection and the previous non-observation of PeV-scale neutrinos in IceCube data. Second, we derive new constraints on the {\nu}SI coupling strength for mediator masses around 100 MeV. Our analysis shows that neutrino telescopes can surpass existing collider constraints in this mass range. In the near future, IceCube-Gen2 is expected to substantially enhance the sensitivity to {\nu}SI over a broader range of parameter space.

We consider a spinning primordial black hole (PBH) described by the Kerr metric in Kerr-Schild polar coordinates (KSP). We derive an analytical expression for the four-vector gravitational potential in the underlying Hermitian Dirac Hamiltonian using these coordinates. This gravitational potential introduces an axial vector term in the Dirac equation in curved spacetime. We find that the magnitudes of the temporal and spatial components of the four-vector gravitational potential are significantly affected by the angle of the position vector of the spinor with respect to the spin axis of the PBH, its radial distance from the PBH, and the strength of the specific angular momentum of the PBH. These potentials modify the effective mass matrix of the neutrino-antineutrino system and significantly affect the transition probabilities during the flavor oscillation of the neutrino-antineutrino system. We then use the transition probability to investigate the quantum speed limit time (QSLT) bound ratio for the two-flavor oscillation of the neutrino-antineutrino system in curved spacetime. This helps us estimate how quickly the initial neutrino flavor state evolves over time under the influence of the gravitational field. Finally, we discuss quantum correlations such as entanglement entropy during the two-flavor oscillation of the neutrino-antineutrino system near a spinning PBH.

We present a numerical relativity study of the accretion properties of a non-spinning black hole in a cosmology driven by a scalar field. The simulations are carried out with a modified moving-puncture gauge condition suitable for cosmological space-times. We considered a scalar field with potential $ V=\lambda \,\varphi^4/4$ and derived the black hole mass growth formula for this scenario using the dynamical horizon framework. As with perturbative studies, we find that the accretion rate $\dot M \propto M^2$ with $M$ the mass of the black hole, and that $\dot M \propto \dot\varphi^2$. We verify that the results of the simulations satisfy the mass growth formula. Unexpectedly, the dynamics of the scalar field in the neighborhood of the black hole is not significantly different from the behavior of the field far away from the hole. We found situations in which the black hole can growth $\sim 15\%$ of its initial mass before the scalar field reaches the bottom of its potential.

We consider the formation of metastable cosmic strings in a left-right unified theory. The produced monopoles are diluted by a stage of F-term hybrid inflation (FHI) which is realized consistently with the SUSY breaking and a global $U(1)$ $R$ symmetry in the context of a $U(1)_{\rm R}\times U(1)_{B-L}$ extension of MSSM. The hidden-sector Kaehler manifold enjoys an enhanced $SU(1,1)/U(1)$ symmetry with the scalar curvature determined by the achievement of a SUSY-breaking de Sitter vacuum without ugly tuning. FHI turns out to be compatible with data -- including the recent ACT results --, provided that the magnitude of the emergent soft tadpole term is confined in the range $(0.1-70)$ TeV, and it is accompanied with the production of cosmic strings. Their dimensionless tension $G\mu_{\rm cs}\simeq(1-11)\cdot10^{-8}$ interprets the present observations from PTA experiments on the stochastic background of gravitational waves. The $\mu$ parameter of MSSM arises by appropriately adapting the Giudice-Masiero mechanism and facilitates the out-of-equilibrium decay of the $R$ saxion at a reheat temperature lower than about $34$ GeV. The SUSY mass scale turns out to lie in the PeV region.

Free tropospheric (FT) nitrogen dioxide (NO2) plays a critical role in atmospheric oxidant chemistry as a source of tropospheric ozone and of the hydroxyl radical (OH). It also contributes significantly to satellite-observed tropospheric NO2 columns, and must be subtracted when using these columns to quantify surface emissions of nitrogen oxide radicals (NOx = NO + NO2). But large uncertainties remain in the sources and chemistry of FT NO2 because observations are sparse. Here, we construct a new cloud-sliced FT NO2 (700-300 hPa) product from the TEMPO geostationary satellite instrument over North America. This product provides higher data density and quality than previous products from low Earth orbit (LEO) instruments, with the first observation of the FT NO2 diurnal cycle across seasons. Combined with coincident observations from the Geostationary Lightning Mapper (GLM), the TEMPO data demonstrate the dominance of lightning as a source of FT NO2 in non-winter seasons. Comparison of TEMPO FT NO2 data with the GEOS-CF atmospheric chemistry model shows overall consistent magnitudes, seasonality, and diurnal variation, with a midday minimum in non-winter seasons from photochemical loss. However, there are major discrepancies that we attribute to GEOS-CF's use of a standard cloud-top-height (CTH)-based scheme for the lightning NOx source. We find this scheme greatly underestimates offshore lighting flash density and misrepresents the diurnal cycle of lightning over land. Our FT NO2 product provides a unique resource for improving the lightning NOx parameterization in atmospheric models and the ability to use NO2 observations from space to quantify surface NOx emissions.

We model the light HESS J1731-347 compact object (of known stellar mass and radius) within Einstein's General Relativity imposing the Karmarkar condition in gravity for anisotropic stars. The three free parameters of the analytic solution are determined imposing the matching conditions at the surface of the star for objects of known stellar mass and radius. Finally, using well established criteria it is shown that the solution is compatible with all requirements for well behaved and realistic solutions. Furthermore, we study the radial oscillation modes, and we compare to the ones corresponding to an isotropic star modeled by the Tolman IV exact analytic solution obtained long time ago. A comparison between the large frequency separations is made as well.

The most general bound binary black hole (BBH) system has an eccentric orbit and precessing spins. The detection of such a system with significant eccentricity close to the merger would be a clear signature of dynamical formation. In order to study such systems, it is important to be able to evolve their spins and eccentricity from the larger separations at which the binary formed to the smaller separations at which it is detected, or vice versa. Knowledge of the precessional evolution of the binary's orbital angular momentum can also be used to twist up aligned-spin eccentric waveform models to create a spin-precessing eccentric waveform model. In this paper, we present a new publicly available code to evolve eccentric, precessing BBHs using orbit-averaged post-Newtonian (PN) equations from the literature. The spin-precession dynamics is 2PN accurate, i.e., with the leading spin-orbit and spin-spin corrections. The evolution of orbital parameters (orbital frequency, eccentricity, and periastron precession), which follow the quasi-Keplerian parametrization, is 3PN accurate in the point particle terms and includes the leading order spin-orbit and spin-spin effects. All the spin-spin terms include the quadrupole-monopole interaction. The eccentricity enhancement functions in the fluxes use the high-accuracy hyperasymptotic expansions from Loutrel and Yunes [Classical Quantum Gravity {\bf 34} 044003 (2017)]. We discuss various features of the code and study the evolution of the orbital and spin-precession parameters of eccentric, precessing BBHs. In particular, we study the dependence of the spin morphologies on eccentricity, where we find that the transition point from one spin morphology to another can depend nonmonotonically on eccentricity, and the fraction of binaries in a given morphology at a given point in the evolution of a population depends on the instantaneous eccentricity.

Despite previous results that rule out the possibility of a static black hole horizon in cosmology we present a black hole metric that has a static horizon whilst mimicking the cosmological behavior at large scales away from the black hole horizon. By a suitable choice of coordinates, we show that it is possible to have a static black hole horizon in cosmology that does not suffer from issues like singularities in curvature invariants at the horizon. The resulting metric is consistent with the Schwarzschild-de Sitter metric for constant Hubble parameter. It should be noted that this metric does not lead to isotropic pressure, i.e. $p_\perp\neq p_r$ for radial distances comparable to Schwarzschild radius, but for large $r$ the stress tensor (assuming general relativity) tends to the cosmological value. Turning to a more realistic framework, we study the cases of black hole$+$cosmological matter and black hole$+$matter$+$cosmological constant where we assume the pressure component takes its cosmological value and show that although the resulting metric does not have a static horizon, the physical radius of the apparent horizon can tend to a constant value for large cosmological time. This metric does not suffer from issues like a singular horizon for asymptotically vanishing Hubble parameter.

Massless or massive particles in unstable orbits around a Kerr black hole exhibit chaotic motion when perturbed. They either plunge into the black hole or escape to infinity after making some oscillations around the equatorial plane. In both of these cases, chaotic motion causes information production. In the case of the photons that escape to infinity, it was recently suggested that this information can be used to resolve the subring structure of the shadow image and obtain more precise data about the black hole mass and spin. Here, we extend this method to obtain more precise results by including the non-equatorial contributions to the Lyapunov exponents. In the other case of massive particles that plunge into the Kerr black hole, we show that the associated Kolmogorov-Sinai entropy derived from the Lyapunov exponents can be interpreted in the context of black hole thermodynamics and that it obeys Bekenstein's bound on the entropy of a physical material system. Thus, the perturbed unstable orbits, either ending inside the black hole or at the observer's screen, have physical consequences.

The minimal warm inflation scenario proposed in Ref. [1] -- featuring an axion-like inflaton coupled to Standard Model (SM) gluons via the standard interaction $\phi G \tilde G$ -- offers a compelling bridge between inflationary dynamics and SM particle content. While the model retains only the inflaton as a beyond-SM field, its original analysis relied on some approximate treatments of warm inflation's (WI) dynamics. Here, we revisit this scenario using WI2easy, a precision computational tool for WI dynamics [2], to rigorously evaluate the model's viability and full range of model's parameters compatible with the observational parameters. Overall, we find that the results of Ref. [1] hold, but with significant differences in the weak and strong dissipative regimes of WI.