Papers for Friday, Jan 30 2026

The Habitable Worlds Observatory (HWO), planned for launch in the 2040s, represents the next major step in exoplanet characterisation. HWO will, for the first time, enable detailed studies of the atmospheres and surfaces of Earth-like exoplanets through high-contrast reflection spectroscopy across the UV, optical, and near-infrared. These wavelength ranges provide access to key molecular absorption features, including O2, O3, H2O, CO2, and CH4, as well as potential surface biosignatures such as the vegetation red edge or ocean glint, making HWO a cornerstone mission for assessing planetary habitability. Clouds are a dominant factor in determining planetary climate and observability, yet their properties remain highly degenerate when constrained using reflected flux alone. Spectropolarimetry, a measure of the polarisation state of reflected light as a function of wavelength and orbital phase, provides a powerful complementary diagnostic. Polarisation is highly sensitive to cloud particle size, composition, shape, vertical distribution, and surface type, enabling degeneracies between atmospheric and surface models to be broken. Numerous studies have demonstrated the value of polarimetry for characterising a wide range of exoplanets, from hot Jupiters to cooler potentially habitable worlds. HWO's proposed instrument suite includes a coronagraph, a high-resolution imager, and a candidate high-resolution spectropolarimeter, offering multiple pathways to exploit polarimetry across diverse planetary regimes. This white paper argues that incorporating polarimetric capability into HWO instruments would significantly enhance the mission's scientific return. We highlight the unique opportunity for UK leadership in both instrument development and theoretical modelling, and advocate for a strong UK role in shaping HWO's polarimetric capabilities to maximise its impact on exoplanet science.

The Python package teareduce has been developed to support teaching activities related to the reduction of astronomical data. Specifically, it serves as instructional material for students participating in practical classes on the processing of astronomical images acquired with various instruments and telescopes. These classes are part of the course Experimental Techniques in Astrophysics, which belongs to the Master's Degree in Astrophysics at the Complutense University of Madrid. The code is publicly available on GitHub, accompanied by a documentation page that includes Jupyter notebooks demonstrating the use of its various classes and functions.

Stars form as molecular clouds fragment into networks of dense cores, filaments, and subclusters. The characteristic spacing of these cores is a key observable imprint of fragmentation physics and is commonly measured using nearest-neighbour (NN) statistics. However, NN separations are derived from projected two-dimensional (2D) positions, while fragmentation occurs in three dimensions (3D). Using spherical and fractal toy models, we show that the standard geometric deprojection factor of $4/\pi\simeq1.27$ is inadequate because projection not only foreshortens separations but also rewires the NN network, while finite angular resolution merges close neighbours and inflates apparent spacings. We quantify these competing biases with Monte Carlo experiments spanning a wide range of morphologies, sample sizes, and effective resolutions. From these we derive an empirical correction factor that depends on both sample size and resolution: for small ($N\lesssim10$) or poorly resolved samples ($\lesssim$10 resolution elements across the field), intrinsic NN spacings exceed projected values by only 20 to 40%, whereas for well-sampled ($N\gtrsim100$), well-resolved data ($\gtrsim$30-50 resolution elements), true 3D separations are typically larger by a factor of $\sim$2. This calibration enables observers to convert measured 2D NN spacings into corresponding 3D estimates, with typical morphology-driven uncertainties of order 30 to 40%, and we demonstrate how it alters inferred fragmentation scales in observed and simulated core populations. [abridged]

A significant fraction of galaxy clusters show central cooling times of less than 1 Gyr and associated central cluster entropies below $30\,\mathrm{keV}\,\mathrm{cm}^2$. We provide a straight forward explanation for these low central entropies in cool core systems and how this is related to accretion onto supermassive black holes (SMBHs). Assuming a time-averaged equilibrium between active galactic nucleus (AGN) jet heating of the radiatively cooling intracluster medium (ICM) as well as Bondi accretion, we derive an equilibrium entropy that scales with the SMBH and cluster mass as $K\propto M_\bullet^{4/3}M_{500\mathrm{c}}^{-1}$. At fixed cluster mass, overly massive SMBHs would raise the central entropy above the cool core threshold, thus implying a novel way of limiting SMBH masses in cool core clusters. We find a limiting mass of $1.4\times10^{10}\,\mathrm{M}_\odot$ in a cool core cluster of mass $10^{15}\,\mathrm{M}_\odot$. We carry out three-dimensional hydrodynamical simulations of an idealized Perseus-like cluster with AGN jets and find that they reproduce the predictions of our analytic model, once corrections for elevated jet entropies are applied in calculating X-ray emissivity-weighted cluster entropies. Our findings have significant implications for modelling galaxy clusters in cosmological simulations: a combination of overmassive SMBHs and high heating efficiencies preclude the formation of cool core clusters.

Tidal disruption events (TDEs) have been proposed as candidate sources of high-energy neutrinos. Successful and choked jets, as well as the accretion disk, corona, wind, and outflow regions in a TDE have been examined and shown to produce TeV - PeV neutrinos. In this work, we use the IceCube 12.5 year high energy starting events (HESE) dataset and perform a maximum likelihood analysis to investigate the spatial and temporal correlations between HESE dataset and a selected sample of 89 TDEs. Our results indicate that the currently observed data do not show any significant correlation and hence is consistent with the background only hypothesis. Using this result, we place constraints on the fraction of TDEs harboring intrinsic jets ($f_{\rm jet}$) and the corresponding isotropic-equivalent cosmic ray (CR) energy ($\mathcal{E}_{\rm CR}$). We note that even with limited statistics, we can constrain the parameter space as $\mathcal{E}_{\rm CR} \lesssim 3 \times 10^{53}$ erg for $f_{\rm jet} \gtrsim 0.6$ at more than 90% C.L. Finally, we discuss the theoretical implications of our results and the limits on the all-sky diffuse neutrino flux from TDEs. With more observational data in the electromagnetic band for TDEs and neutrino observations from IceCube and KM3NeT, our analysis can be used to place stringent constraints on physical parameters associated with TDEs.

One of the main properties of galactic bars is their rotation (or pattern) speed, which is driven by both internal galactic properties, as well as external interactions. To assess the influence of these internal and external drivers on bar rotation in a cosmological setting, we use the Auriga suite of cosmological hydrodynamical zoom-in simulations. We calculate the bar pattern speed and the bar rotation rate - the ratio of corotation radius to bar length - at the time of bar formation and at z=0, and compare these to bar age, bar strength, baryon dominance, galaxy stellar mass, and the history of external galaxy interactions. We find that galaxies which are more baryon dominated at z=0 - and which lie above the observed stellar mass-halo mass abundance matching relation - host faster bars, while more dark matter dominated galaxies host slower bars. Baryon-dominated galaxies also form their bars earlier and their rotation rates stay constant or even decrease over time; this leads to older bars being faster than their younger counterparts - in contrast to the expectation of bar slow-down from dynamical friction imparted by the dark matter halo. We also find a trend in stellar mass, with 'faster' bars being hosted in more massive galaxies, which could be driven by the underlying higher baryon-dominance of more massive galaxies. Furthermore, we find that external interactions, such as mergers and flybys, correlate with lower bar rotation rates, particularly for strong interactions that occur around bar formation time. This correlation is relatively weak, leaving internal baryon-dominance as the main driver of fast bar rotation rates.

The Haro 11 galaxy merger is the closest known Lyman Continuum (LyC) leaker and a strong Lyman-$\alpha$ (Ly$\alpha$) emitter, making it an important analogue of the high-$z$ galaxies that reionised the early Universe. To investigate how Haro 11's properties arise, we perform a radiation hydrodynamics simulation of the merger, and create mock observations of LyC, Ly$\alpha$, and H$\alpha$, from which we compute their luminosities ($L$) and escape fractions ($f_{\rm esc}$). We track these quantities along multiple sightlines as the two progenitor galaxies merge, from the first interaction until the system resembles present-day Haro 11. We find that $L$ and $f_{\rm esc}$ vary by 1-2 orders of magnitude for LyC due to sightline variations. At the two pericentre passages, the total $f_{\rm esc}^{\rm LyC}$ increases by roughly an order of magnitude. Conversely, $f_{\rm esc}^{\rm Ly\alpha}$ shows a moderate increase at the pericentre passages, which affects the inference of LyC properties from Ly$\alpha$. We attribute this to a displacement of the LyC-emitting stars relative to the \Lya-emitting gas, combined with an increased density from gas compression. Furthermore, $f_{\rm esc}^{\rm LyC}$ is boosted during star formation bursts, likely due to stellar feedback. As direct comparison with Haro 11, the simulation qualitatively matches its morphology and luminosities. We find that among the dense stellar knots, knot C is the main contributor to both intrinsic and escaping LyC emission. Additionally, the Ly$\alpha$ spectra displays distinct features found in observations, implying similar gas conditions are present.

The Euclid Early Release Observations (ERO) programme targeted the Perseus galaxy cluster in its central region over 0.7deg$^2$. We combined the exceptional image quality and depth of the ERO-Perseus with FUV and NUV observations from GALEX and AstroSat/UVIT, as well as $ugrizH\alpha$ data from MegaCam at the CFHT, to deliver FUV-to-NIR magnitudes of the 87 brightest galaxies within the Perseus cluster. We reconstructed the star-formation history (SFH) of 59 early-type galaxies (ETGs) within the sample, through the spectral energy distribution (SED) fitting code CIGALE and state-of-the-art stellar population (SP) models to reproduce the galactic UV emission from hot, old, low-mass stars (i.e. the UV upturn). In addition, for the six most massive ETGs in Perseus [stellar masses $\log_{10}(M_{\ast}/M_{\odot}) \geq 10.3$], we analysed their spatially resolved SP through a radial SED fitting. In agreement with our previous work on Virgo ETGs, we found that (i) the majority of ETGs needs the presence of an UV upturn to explain their FUV emission, with temperatures $\langle T_{\rm UV}\rangle$~33800 K; (ii) ETGs have grown their stellar masses quickly, with SF timescales $\tau\lesssim 1500$ Myr. We found that all ETGs in the sample have formed more than about 30% of their stellar masses at z~5, up to ~100%. At z~5, the stellar masses of the most massive nearby ETGs, which have present-day stellar masses $\log_{10}(M_{\ast}/M_{\odot})\gtrsim 10.8$, are then found to be comparable to those of the red quiescent galaxies observed by JWST at similar redshifts (z>4.6). This study can be extended to ETGs in the 14000 deg$^2$ extragalactic sky that will soon be observed by Euclid, in combination with those from other major upcoming surveys (e.g. Rubin/LSST), and UV observations, to ultimately assess whether the nearby massive ETGs represent the progeny of the massive high-z JWST red quiescent galaxies.

We present an analysis of 22,656 narrow-line Seyfert 1 galaxies (NLSy1s) from Sloan Digital Sky Survey (SDSS) DR17 ($0.1\leq z\leq 0.9$), identifying a sample of spectroscopically anomalous sources. These anomalies were detected via the spectroscopic quasar anomaly detection (SQuAD) algorithm, which employed principal component analysis and hierarchical k-means clustering. Various physical diagnostic analysis were performed such as the color excess ($E_{(B-V)}$) calculations, Wide-field Infrared Survey Explorer (WISE) color analysis, probing O III, equivalent width as an inclination indicator, the BPT diagram and eigenvector 1 diagram. We detected 620 anomalous NLSy1s classified into two groups i.e. 246, Red NLSy1s, exhibiting host galaxy dominated spectra with a low luminosity active galactic nuclei (AGN) core revealed by the emission line widths. Another set of 374 Blue NLSy1s, strongly luminous galaxies with enhanced AGN activity, bluer continuum as compared to a typical NLSy1 and stronger Fe II emission. Finally, the third group of 257 outliers, identified as Intermediate Seyferts, a class of Seyfert galaxies identified by composite emission profiles, and extremely strong emission lines paired with virtually no continuum. These sources also exhibit rare and high ionization emission lines unseen in any other NLSy1 spectra (e.g. [Ne V]$\lambda3345$, Ne V $\lambda3426$, Ne III $\lambda3869$ etc). We conclude that the differentiating factor between red and blue NLSy1s is not dust obscuration or orientation effect, but intrinsic distinction in AGN activity. The resulting sample is presented as a value-added catalog.

Non-inflationary sources of gravitational waves in the early Universe generically predict causality-limited tensor power spectra at low frequencies. We report the first-ever constraints on such sources based on cosmic microwave background (CMB) $B$-mode polarization measurements. Using data from BICEP/Keck, SPTpol, SPT-3G, Planck, and WMAP, we constrain the amplitude of an early causal tensor (ECT) power spectrum parameterized by $r_{ect}$, the ratio of causal tensor power to total scalar power at $k~=~0.01$ Mpc$^{-1}$, and obtain a 95% CL upper limit of $r_{ect}<$ 0.0077. Since $r_{ect}$ can easily be related to the parameters of a given theory, our bound robustly constrains a broad class of well-motivated gravitational wave sources in the early universe, including first-order cosmological phase transitions, enhanced small-scale density perturbations, and various topological defects. Finally, we translate our limit into a bound on the present-day energy density in gravitational waves at ultra-low frequencies otherwise inaccessible to traditional gravitational wave detection strategies, including pulsar timing arrays, interferometers, and resonant cavities.

Here we present the results of our search through the Early Data Release (EDR) of the Dark Energy Spectroscopic Instrument (DESI) for extreme coronal line emitters (ECLEs) - a rare classification of galaxies displaying strong, high-ionization iron coronal emission lines within their spectra. With the requirement of a strong X-ray continuum to generate the coronal emission, ECLEs have been linked to both active galactic nuclei (AGNs) and tidal disruption events (TDEs). We focus our search on identifying TDE-linked ECLEs. We identify three such objects within the EDR sample, highlighting DESI's effectiveness for discovering new nuclear transients, and determine a galaxy-normalized TDE-linked ECLE rate of $R_\mathrm{G}=5~^{+5}_{-3}\times10^{-6}~\mathrm{galaxy}^{-1}~\mathrm{yr}^{-1}$ at a median redshift of z = 0.2 - broadly consistent with previous works. Additionally, we also identify more than 200 AGNs displaying coronal emission lines, which serve as the primary astrophysical contaminants in searches for TDE-related events. We also include an outline of the custom python code developed for this search.

We examine the possibility that evidence for a non-zero value of the CP violating $ \frac{e^2 }{32\pi^2}\theta_{EM} \int d^4 x {\vec E}\cdot {\vec B}$ coupling might be extracted from Event Horizon Telescope observations of the black holes SgA* and M87*. The Fischler-Kundu\cite{FK} effect predicts a universal Hall current in the relaxation of charge falling onto the black hole horizon. We argue that this leads to a non-zero value of a certain CP-violating observable ${\cal C}$, defined below. The effect can be masked by parity violating plasma currents. In particular, evidence for polarization flips \cite{flip} in the signals from M87* indicate strong plasma effects in the data. We suggest that time averaging the data over periods including the flip might leave over a residual that would be an indicator of the FK signal. In addition, similarities in the polarization patterns between the two very different black holes, and a part of the signal that is uniform in frequency, might enable us to distinguish the universal topological signal from source and frequency dependent plasma effects. Current data does not appear to be sufficient to perform such a test.

Detecting and characterizing the atmospheres of rocky exoplanets has proven to be challenging for JWST. Transit spectroscopy of the TRAPPIST-1 planets has been impacted by the effects of spots and faculae on the host star. Secondary eclipses have detected hot rocks, but evidence for atmospheres has been difficult to obtain. However, there is a third option that we call color phase curves. This method will apply to synchronously rotating non-transiting planets as well as transiting planets. A color phase curve uses photometry at a long-IR wavelength near the peak of the planetary thermal emission (e.g., 21 microns) divided by photometry at a shorter wavelength where the star dominates more strongly (e.g., 12 microns). We avoid wavelengths having potentially strong molecular absorption (e.g., 15 microns) to minimize degeneracies in the color phase curve, and we aim to detect and characterize the planetary atmosphere via its longitudinal heat transfer. The ratio of two wavelengths observed nearly simultaneously is designed to isolate thermal emission from the planet, discriminate against the star, and largely cancel instrumental systematic effects. Moreover, we show that invoking mass-radius relations, and using self-consistent physical models, will permit the longitudinal heat transfer to be measured independent of the orbital inclination. Radial velocity surveys are detecting many new exoplanets, including temperate rocky worlds with Earth-like masses. Most of those planets will not transit, but color phase curves have the potential to detect and characterize their atmospheres.

Many early universe scenarios predict post-inflationary tensor perturbations from causality-limited, sub-horizon sources. While the microphysical details may differ, as long as these sources are bounded in duration and correlation length, their tensor power spectra exhibit a universal scaling behavior at small wavenumber: $P_h(k) \propto k^3$, corresponding to white noise on super-horizon scales at the time of production. If these early causal tensor sources (ECTs) exclusively produce gravitational waves before redshift $z \sim 10^5$, this scaling is realized on all of the scales observed in the cosmic microwave background (CMB), and thus yields a universal multipole distribution for the $B$-mode angular power spectrum. Unlike the scale-invariant distributions of inflationary $B$ modes, ECTs generically predict enhanced power on small scales and suppressed power on large scales, which allows these source classes to be distinguished given measurements over a sufficient range of angular scales. In this paper, we introduce a unified framework for characterizing ECTs and demonstrate how their universal infrared scaling manifests in low-frequency observables, including CMB $B$ modes and stochastic gravitational wave spectral densities. We illustrate this mapping with representative case studies of this universality class involving first-order phase transitions, topological defects, and enhanced scalar perturbations, which source tensor modes at second order in perturbation theory.

We investigate the hypothesis that the evolution of the Universe can be described by a single dark fluid whose effective equation of state (EoS), $\omega_{\rm{eff}}$, is a linear combination of a logarithmic term and a power law term, both involving the Lambert $W$ function. This particular form of EoS was first proposed by S. Saha and K. Bamba in 2019 and has two parameters, $\theta_1$ and $\theta_2$, which must be determined from observations. To this end, we place limits on these parameters by combining recent baryon acoustic oscillation (BAO) data -- including measurements from the Dark Energy Spectroscopic Instrument (DESI) -- with Type Ia supernova observations from the Pantheon+ compilation, along with direct determinations of the Hubble parameter. From this combined analysis, we obtain a best-fit value for the Hubble parameter, $H_0 = 67.4 \pm 1.2~\text{km\,s}^{-1}\text{Mpc}^{-1}$, while current measurements of the sound horizon at the baryon drag epoch yield $r_d = 146\pm 2.5$~Mpc. Furthermore, we study the evolution of the deceleration parameter, the effective EoS, and the jerk parameter, and support our findings using the $Om(z)$ diagnostic. The model exhibits noticeable deviations from the predictions of the concordance $\Lambda$CDM model. Despite these differences, our results indicate that the model provides a coherent description of late-time cosmic evolution and the observed accelerated expansion of the Universe. Finally, we assess the observational viability of the model using information criteria, particularly the Akaike Information Criterion (AIC) and the Bayesian Information Criterion (BIC), and compare these results with those obtained for the $\Lambda$CDM model, which serves as our reference.

We measure the radial $g-i$ colour profiles of $\sim$32,000 galaxies drawn from the Hyper Suprime-Cam Subaru Strategic Program optical imaging survey, including 1415 exhibiting tidal features. We compare the colour profiles of galaxies with and without tidal features to extract information about the properties of the mergers that created these features. We find negative colour gradients for both galaxies with and without tidal features and find that tidal feature-hosting red sequence galaxies have redder outskirts than their non-tidal feature hosting counterparts, consistent with the outskirts of these galaxies being dominated by stars accreted from gas-poor minor mergers. We find decreasing mass ratios of tidal features-to-host galaxy with increasing galaxy stellar mass, suggesting that less massive galaxies undergo mergers with companions closer in mass than more massive galaxies. Galaxies exhibiting streams have bluer outskirts than those hosting shells, and shells around red sequence galaxies tend to be more massive and have higher mass ratios to their hosts than streams, consistent with streams being formed from mergers with satellites less massive than those responsible for shells. The agreement between our findings and those of other observational and simulation-based works confirms the validity of our methodology and highlights the value of tidal features colours as a probe into the process through which galaxies evolve.

We present a study of the linear polarization properties of radio sources within the 10 deg$^2$ Wide Chandra Deep Field South (W-CDFS) in S-band (2-4 GHz). Our W-CDFS image has an angular resolution of 15 arcsec and a 1$\sigma$ RMS in Stokes $I$ of $\approx$50 $\mu$Jy/beam. We detect 1920 distinct source components in Stokes $I$ and 175 in linear polarization. We examine the polarized source counts, Faraday Rotation measures, and fractional polarization of the sources in the survey. We show that sources with a total intensity above $\approx$10mJy have a mean fractional polarization value of $\approx$3% from modeling the polarized counts. We also calculate an estimate for the limit on the fractional polarization level of sources with a total intensity below 1mJy (mostly star-forming galaxies) of $\stackrel{<}{_{\sim}}$3% using stacking. The mean Faraday Rotation we measure is consistent with that due to the Milky Way. We also show that fractional polarization is correlated with in-band spectral index, consistent with a lower mean fractional polarization for the flat-spectrum population. In addition to characterizing the S-band polarization properties of sources in the W-CDFS, this study will be used to validate the shallower, but higher angular resolution S-band polarimetric information that the VLA Sky Survey will provide for the whole sky above Declination -40 degrees over the next few years.

We report the results of a pilot study that searched for dual active galactic nuclei (AGN) in local ($z<$0.25) galaxies hosting double-peaked narrow emission lines in their optical spectra. We present high-resolution $L-$band (1.5 GHz or 18 cm) continuum images from the Very Long Baseline Array (VLBA) as well as WFC3/IR F160W images from the \textit{Hubble Space Telescope} of two candidate dual AGN systems: J0948+6848 and J1223+5409. In both targets, we detected compact non-thermal radio emission that is approximately co-spatial with the near-infrared AGN. Both systems host two high brightness temperature ($>10^{8}$ K) radio sources that indicate the presence of either a parsec-scale-separation dual AGN ($d_{\text{sep}} \sim 90$ pc and $\sim 56$ pc, respectively) or a radio jet. Matched-resolution multi-band radio observations are necessary to further characterize the AGN activity in these systems.

The stellar initial mass function (sIMF) is often treated as a stochastic probability distribution, yet such an interpretation implies Poisson noise that is inconsistent with growing observational evidence. In particular, the observed relation between the mass of the most massive star formed in an embedded cluster and the cluster's total stellar mass supports a deterministic sampling process, known as optimal sampling. However, the physical origin of optimal sampling has not been formally established in the literature. In this work, we show that the stellar mass distribution implied by optimal sampling emerges from applying the Maximum Entropy principle to the fragmentation of star-forming clumps, whose structure is set by density-dependent cooling in the optically thin regime. Here, the maximum entropy leads to unbiased distributions. By applying calculus of variations to minimize the entropy functional obtained assuming fragmentation, we recover the power-law form of the sIMF, and we show that any distribution deviating from the sIMF violates the Maximum Entropy principle. This work provides a first-principles foundation for the deterministic nature of star formation. Thus, the sIMF is the distribution resulting from a maximally unbiased system.

Infrared observations can probe photometric variability across the full evolutionary range of young stellar objects (YSOs), from deeply embedded protostars to pre-main-sequence stars with dusty disks. We present 3-8 micron light curves extending 27 years from 1997 to 2024 obtained with three space-based IR telescopes: ISO, Spitzer and WISE. Although unevenly sampled with large gaps in coverage, these light curves show variability on time scales ranging from days to decades. We focus on the Spitzer-identified YSOs with disks and envelopes that exhibit variations of a factor of two or more in this wavelength range. We identified seven YSOs where the light curves are dominated by bursts of sustained (> 5 yr) high flux, including four that show a steep decay ending the burst and three that are ongoing as of the final observation. We find six YSOs that are undergoing declines, which may be the end of bursts that began before 1997. The most common form of variability, exhibited by 26 YSOs in our sample, show variations over time intervals of years to months but do not exhibit sustained bursts or fades. The Spitzer [3.6]-[4.5] and WISE [3.5]-[4.6] colors either increase or remain constant with increasing brightness, inconsistent with dust extinction as being the primary source of the large-amplitude variability.

Targeted astrophysical surveys are limited by the amount of telescope time available, which makes it impossible to observe every single object of interest. In order to maximize the scientific return, we need a well thought strategy for selecting the observational targets, in our case exoplanets. This study evaluates various strategies for selecting exoplanet targets within limited observation windows, focusing specifically on the selection of exoplanet targets for Tier 2 transit spectroscopy with ESA's upcoming Ariel mission. We define three distinct selection criteria -- sample size, variance, and leverage -- and translate them into objective functions compatible with modern optimization algorithms. Specifically, we test five heuristics for maximizing sample leverage: leverage greedy, simulated annealing, K-means clustering, regular classes, and quantile classes. The performance of these methods is demonstrated through three practical exercises across one, two, and three parameters of diversity. Each criterion represents a unique trade-off between sample size, diversity, and total observation time. While a time-greedy approach maximizes the quantity of planets, it fails to capture diversity. Conversely, variance-greedy selection prioritizes diversity but introduces significant drawbacks: it oversamples rare cases and undersamples typical planets, ultimately reducing the total number of targets observed. Leverage-based selections emerge as the most effective middle ground, successfully balancing sample diversity with a robust sample size. This work supports the broader community effort to ensure that Ariel delivers the most diverse and scientifically valuable sample of exoplanet atmospheres within mission limits.

The Nuclear Stellar Disc (NSD) of the Milky Way is a dense, rotating stellar system in the central 200 pc. The NSD is thought to be primarily fuelled by bar-driven gas inflows from the inner Galactic disc. As part of the LEGARE project, we construct the first chemical evolution models for the NSD using a Bayesian approach tailored to reproduce the observed metallicity distribution functions (MDFs) and compared with the available abundance ratios for Mg, Si, Ca relative to Fe. We adopt a state-of-the-art chemical evolution model in which the gas responsible for the formation of the NSD is assumed to be driven by the Galactic bar-induced inflows. The chemical composition of the accreted material is assumed to reflect that of the Galactic disc at a radius of 4 kpc. A Bayesian MCMC framework is used to fit the MDFs of different samples of NSD stars. If we take the NSD data at face value, without considering a possible contamination from bulge stars, we find that a formation scenario based on the inner disc flowing gas is inconsistent with the low metallicity tail of the observed MDF. This is because the inner disc metallicity, at the epoch of bar formation, was already near solar. On the other hand, models invoking dilution from additional metal-poor inflows successfully reproduce the observations. The best-fit model requires inflow metallicity 5 times lower than the inner disc and a moderate star formation efficiency. The same model successfully reproduces the observed [$\alpha$/Fe] vs. [Fe/H] ratios and predicts a star formation history consistent with the most recent estimates. However, if we assume that the MDF is contaminated by metal poor bulge stars and restricted to [Fe/H] > -0.3 dex, gas dilution is no longer required. In this case, the best-fit model has a very low star formation efficiency and a mild galactic wind.

Compact radio sources such as pulsars and FRBs undergo scintillation in the interstellar medium (ISM) when scattered images interfere at the observer. ``Scintillometry'' refers to the range of techniques to extract astrometric information -- such as the angular positions of the images and distances to the scattering screen and source -- from scintillation observations. Pulsar scintillometry has proven to be a powerful technique, revealing rich and unexpected scattering phenomenology in the ISM and also shedding light on the emission physics of pulsars. FRB scintillometry stands to be a similarly powerful probe of FRB emission, as well as structure on tiny scales in ionized media beyond our galaxy, such as the circumgalactic medium (CGM). However, nascent FRB scintillation studies are far from the sophisticated lensing geometry reconstructions that have been performed for scintillating pulsars. In this paper, we introduce a novel theoretical framework for scintillometry, demonstrating that the full astrometric content of scintillation observations is contained within a single underlying observable: the instantaneous spatial wavefield. We relate the instantaneous spatial wavefield to more familiar concepts from the pulsar scintillometry literature, such as the dynamic spectrum. Using this framework, we discuss prospects and limitations for FRB scintillometry, towards the goal of full astrometric reconstructions of FRB lensing geometries. We show how key degeneracies in two-screen scattering measurements can be ameliorated. In addition, we discuss the possibility of inferring dispersion measure gradients across scintillation screens, which may shed light on the highly unconstrained physics of the cool CGM phase on tiny ($\sim 100\,{\rm au}$) scales.

Blue lurkers (BLs) are rejuvenated main-sequence stars hidden among normal main-sequence stars on color-magnitude diagrams of star clusters. In comparison, the blue straggler stars, formed via similar mass transfers or mergers, occupy a distinct space in the color-magnitude diagrams. We compile a list of BL candidates in open clusters using available rotation catalogs. BLs can be identified using either unusually faster rotation compared to similar mass stars, which is a signature of recent accretion, or the presence of a companion, which can only be formed by mass donation, e.g., an extremely low mass white dwarf. Here, we searched for fast-rotating stars on the main sequence of open clusters using Kepler, TESS, and spectroscopic rotation indicators, such as rotation periods and $v\sin i$ measurements. We identified 97 new BL candidates across 35 open clusters, almost tripling the previously known sample of 36. Based on the estimated completeness of $\approx$3\%, thousands of BLs are likely hidden within the cluster population. Detailed spectroscopic and time-series analyses will be essential to confirm their mass-transfer histories.

Observations of exoplanetary atmospheres provide critical insights into their chemical composition, formation and evolution history. Ultra-hot Jupiters serve as excellent targets for atmospheric characterization; studies of these planets may yield key understanding of gas giant's formation and evolution history. We present a thermal emission study of WASP-33 b's dayside atmosphere, based on two secondary eclipse observations with CFHT/WIRCam in two specific narrow band filters, namely the CO and CH4$_{\rm on}$ filters, and archival data with HST/WFC3 and Spitzer. Stellar pulsations of the host star induce some quasi-periodic photometric variations, particularly in the CH4$_{\rm on}$ band, which are modelled and corrected in the high-precision differential light curves. An eclipse depth of $1565.2^{+228.6}_{-237.5}$ ppm and $914.3^{+56.1}_{-57.0}$ ppm is determined for the CO and CH4$_{\rm on}$ bands, respectively. Combined with HST/WFC3 and Spitzer data, our joint retrieval of WASP-33 b's dayside atmosphere reveals a high metallicity ([Fe/H] $= 1.52^{+0.35}_{-0.52}$), high C/O ratio (C/O $= 0.78^{+0.03}_{-0.04}$), and a thermal inversion layer, suggesting a formation history involving metal-rich gas accretion. We confirm the presence of the molecules H$_{2}$O, H$^{-}$ and CO, and report a tentative detection of TiO in the dayside atmosphere of WASP-33 b. Future higher precision observations with JWST may provide better understand constraints on the chemical abundances of oxygen and refractory element abundances to better WASP-33 b's formation and evolutionary pathway.

The origin of the radio emission in active galactic nuclei (AGN) is still debated. Multiple physical mechanisms can contribute to the spectrum at these frequencies, including relativistic jets, the jet base, outflows, star formation, and synchrotron emission from the hot corona. Recently, new extreme radio variability has been observed in the class of low-mass/high-Eddington AGN known as narrow-line Seyfert 1 (NLS1) galaxies, suggesting that another, more exotic mechanism may also play a role, especially at frequencies above 10 GHz. To investigate this relatively unexplored area of the radio spectrum, we observed a sample of 50 NLS1s with the Karl G. Jansky Very Large Array (JVLA), and 20 of them were observed twice. In this sample, 24 sources were not detected, while the others are typically characterized by a steep spectrum that can be modeled with a power law. We also identified two new candidate jetted NLS1s, including a high-frequency peaker, which is an extremely young relativistic jet. We found no significant variability in the sources observed twice. We conclude that the radio spectrum of NLS1s is typically dominated by optically thin emission, likely from low-power outflows, or by circumnuclear star formation, with a limited contribution from relativistic jets. Further studies at different spatial scales and at other wavelengths are necessary to fully constrain the origin of the radio emission in this class of active galaxies.

Characterizing protostellar outflows is fundamental to understanding star formation feedback, yet traditional methods are often hindered by projection effects and complex morphologies. We present a multi-modal deep learning framework that jointly leverages spatial and spectral information from CO observations to infer outflow mass, inclination, and position angle ($PA$). Our model, trained on synthetic ALMA observations generated from 3D magnetohydrodynamic simulations, utilizes a cross-attention fusion mechanism to integrate morphological and kinematic features with probabilistic uncertainty estimation. Our results demonstrate that Vision Transformer architectures significantly outperform convolutional networks, showing remarkable robustness to reduced spatial resolution. Interpretability analysis reveals a physically consistent hierarchy: spatial features dominate across all parameters, whereas spectral profiles provide secondary constraints for mass and inclination. Applied to observational ALMA data, the framework delivers stable mass and $PA$ estimates with exceptionally tightly constrained inclination angles. This study establishes multi-modal deep learning as a powerful, interpretable tool for overcoming projection biases in high-mass star formation studies.

Comet C/2023 P1 (Nishimura) was observed by the Solar Orbiter Heliospheric Imager (SoloHI), onboard the Solar Orbiter spacecraft, from 2023 September 1 to 14. During this period, the ion tail of the comet exhibited continual fluctuations and four tail disconnection events (TDEs), each coinciding with the passage of a coronal mass ejection (CME). In this work, we report on the ion tail dynamics of the best observed TDE, which occurred on September 11. The SoloHI white-light images reveal an abrupt bending, subsequent kinks, and severing of a downstream portion of the pre-existing ion tail. The onset of disconnection occurred $\sim$6.5 hours after the projected passage of the CME leading edge in the images, consistent with a CME flank encounter. After the disconnection, the ion tail reformed within $\sim$24 hours, with a regrowth rate of $\sim$86$\pm7~\mathrm{km\,s^{-1}}$, indicating the rate at which newly ionized material forms along the magnetic field draped around the comet's coma. After the TDE, the detached tail drifted anti-sunward at an estimated speed of $\sim$295$\pm20~\mathrm{km\,s^{-1}}$, comparable to the local CME flank's speed, suggesting that the severed plasma was most likely carried away from the comet by the CME. This study provides the first direct, quantitative characterization of comet-CME interactions and the subsequent regrowth phase of a cometary TDE. These measurements were achievable by SoloHI's unique inner-heliospheric coverage, thanks to a combination of high photometric sensitivity, short exposure times, and a wide field of view that preserves the fine-scale tail dynamics.

We have discovered a compact hierarchical triple main-sequence star system, which is cataloged as Gaia DR3 1010268155897156864 or TIC 21502513. Hereafter, we call it ``G1010''. G1010 consists of a primary (the most massive) star and inner binary that orbit each other. The primary star is a $0.85_{-0.03}^{+0.03}\;{\rm M}_\odot$ main-sequence (MS) star, and the inner binary components are $0.63_{-0.02}^{+0.02}$ and $0.61_{-0.02}^{+0.02}\;{\rm M}_\odot$ MS stars. The outer and inner orbital periods are $277.2_{-1.3}^{+1.6}$ and $\sim 18.26$ days, respectively. G1010 is categorized as a single-lined spectroscopic binary, and its orbital solution indicates that G1010 possibly accompanies a massive compact object, such as a neutron star or massive white dwarf. In order to confirm the presence of a massive compact object, we have performed several-times low signal-to-ratio (SNR) and one-time high SNR spectroscopic observations, and determined the outer orbital parameters. Moreover, we have deeply analyzed the high SNR spectroscopic data, and found that G1010 accompanies not a massive compact object, but an inner binary. We have investigated G1010's light curve in Transiting Exoplanet Survey Satellite (TESS), and concluded that the inner binary is actually an eclipsing binary, not included in TESS Eclipsing Binary Stars. We have obtained the inner orbital parameters from the TESS light curve. G1010 is similar to compact hierarchical triple star systems previously discovered by eclipse timing variation analysis. Our discovery has shown that such triple star systems can be discovered by combination of low- and high-SNR spectroscopic observations with the help of Gaia DR3 and the upcoming Gaia DR4/DR5.

We present echelle spectra of TIC 399725538 obtained in Korea and Thailand to investigate the physical properties and evolutionary scenarios of EL CVn-type binaries. The time-series spectra yielded the radial velocities (RVs) of the primary component and its atmospheric parameters, $T_{\rm eff,A}=7194\pm70$ K and $v_{\rm A}\sin i=68\pm9$ km s$^{-1}$. Joint modeling of our RVs and the archival TESS data provided component masses of $M_{\rm A}=1.930\pm0.054$ $M_\odot$ and $M_{\rm B}=0.211\pm0.005$ $M_\odot$, radii of $R_{\rm A}=1.922\pm0.020$ $R_\odot$ and $R_{\rm B}=0.207\pm0.005$ $R_\odot$, and luminosities of $L_{\rm A}=8.87\pm0.39$ $L_\odot$ and $L_{\rm B}=0.546\pm0.034$ $L_\odot$. The surface gravity and distance derived from this modeling are consistent with the model-independent $\log g_{\rm B}$ obtained from the direct observables and with the Gaia distance. The third light of $l_3=0.136\pm0.003$ comes mainly from two neighboring stars, TIC 399725539 and TIC 399725544. Comparison with stellar models indicates that TIC 399725538 A lies within the instability region of $\delta$ Sct$-$$\gamma$ Dor hybrids, whereas its extremely low-mass companion is markedly underluminous compared to theoretical white dwarf (WD) counterparts. Multifrequency analysis of the binary-subtracted residual lights revealed three significant signals, two of which correspond to aliases at two and four times the dominant frequency $f_1=0.0848$ day$^{-1}$. The $f_1$ frequency is likely a $\gamma$ Dor-type pulsation arising from the brighter A component, though further confirmation is required. Our results demonstrate that TIC 399725538 is a short-period EL CVn system belonging to the thick-disk population, consisting of a main-sequence $\gamma$ Dor pulsator and a helium-core WD precursor formed through stable mass transfer.

This study reports detection of rare p-process isotopes within presolar grains. Presolar grains are relic dust grains from dying stars. These microscopic dust particles are found in primitive solar system materials. Their distinct isotopic compositions record the nucleosynthetic processes in their parent stars and the Galactic chemical environment in which these stars formed. We studied presolar graphite grains of high-density type from the Murchison meteorite and found five grains with subgrains that show enrichments in 84Sr compared to the solar abundance.84Sr is the neutron-deficient isotope of strontium that can be produced in the deep oxygen-rich interior of high-mass stars that end their lives as core-collapse supernovae. The observed 84Sr excesses cannot be produced in low-mass asymptotic giant branch stars, the source of most high-density presolar graphites found in meteorites. High-density graphites with embedded 84Sr excesses are, instead, compatible with a core-collapse supernovae origin. The graphite subgrains condensed from carbon-rich materials in the outer layers of core-collapse supernovae, where 84Sr was destroyed by neutron captures during hydrostatic evolution of the stars and their final explosion. Based on current theoretical stellar models, a few percent of contribution from the inner regions of core-collapse supernovae, which are enriched in p-process nuclides, to the outer carbon-rich regions is the most likely explanation for the observed enrichment of 84Sr in the subgrains of the high-density graphites. In this study, we present the first observational evidence that core-collapse supernovae produce and eject isotopes made by the p-process.

Interstellar object (ISO) astronomy has rapidly emerged over the past decade as a new frontier in planetary astrophysics, yet current observations remain limited by short visibility windows, inference degeneracies and fragmented follow-up capabilities. We argue that these constraints are structural rather than incidental and motivate a coordinated, end-to-end observational strategy for future ISO studies. We propose the Comprehensive ISO Network (CISON) which combines dual hemisphere wide-field discovery with rapid high resolution characterization and selective escalation to interceptor missions. By coupling this architecture to the differential formulation of the Loeb Scale, ISO classification and risk assessment become predictive rather than reactive. This framework transforms ISO astronomy into a mature, scalable discipline capable of maximizing scientific return and informing planetary defense in the coming decades.

Empirical solar wind speed models play an important role in enabling space weather forecasting with low computational cost. Among these, one model called WS model is based on the asymptotic expansion factor. However, it is known that it fails in the case of pseudostreamers. In this study, as a first step toward constructing a solar wind speed empirical model based on physical parameters, we investigated the effect of the radial profile of flux-tube shape on the solar wind speed using one-dimensional numerical simulations. In the simulations, ad hoc Alfvén waves are injected from the photosphere at $r=R_\odot$ as the energy source, and the MHD equations are solved out to the interplanetary space at $r=70R_\odot$ to reproduce solar wind acceleration. As a result, even when the coronal base magnetic field and the asymptotic expansion factor are fixed, the final solar wind speed varies by approximately 300 km s$^{-1}$ depending on changes in the expansion height or non-monotonic expansion. Additionally, across all simulations performed, a better correlation is found with the quantities that reflect the information about the radial profile of flux-tube shape than the asymptotic expansion factor. Our results suggest that, as a physical characteristic parameter of the solar wind speed, an operation that can account for the expansion factor throughout the corona is necessary.

We develop a machine learning model based on a structured variational autoencoder (VAE) framework to reconstruct and generate neutron star (NS) equations of state (EOS). The VAE consists of an encoder network that maps high-dimensional EOS data into a lower-dimensional latent space and a decoder network that reconstructs the full EOS from the latent representation. The latent space includes supervised NS observables derived from the training EOS data, as well as latent random variables corresponding to additional unspecified EOS features learned automatically. Sampling the latent space enables the generation of new, causal, and stable EOS models that satisfy astronomical constraints on the supervised NS observables, while allowing Bayesian inference of the EOS incorporating additional multimessenger data, including gravitational waves from LIGO/Virgo and mass and radius measurements of pulsars. Based on a VAE trained on a Skyrme EOS dataset, we find that a latent space with two supervised NS observables, the maximum mass $(M_{\max})$ and the canonical radius $(R_{1.4})$, together with one latent random variable controlling the EOS near the crust--core transition, can already reconstruct Skyrme EOSs with high fidelity, achieving mean absolute percentage errors of approximately $(0.15\%)$ for $(M_{\max})$ and $(R_{1.4})$ derived from the decoder-reconstructed EOS.

In generalized Proca theories, characterized by a vector field with broken $U(1)$ gauge invariance, late-time cosmic acceleration can be realized with a dark energy equation of state in the regime $w_{\rm DE} < -1$. In such scenarios, however, a phantom-divide crossing, as recently suggested by DESI observations, is not achieved without encountering theoretical inconsistencies. We incorporate a canonical scalar field with a potential, in addition to the vector field, and show that the phantom-divide crossing from $w_{\rm DE} < -1$ to $w_{\rm DE} > -1$ can occur at low redshifts. We propose a minimal model that admits such a transition and identify the region of parameter space in which all dynamical degrees of freedom in the scalar, vector, and tensor sectors are free from ghost and Laplacian instabilities. We further investigate the evolution of linear cosmological perturbations by applying the quasi-static approximation to modes well inside the Hubble radius. The dimensionless quantities $\mu$ and $\Sigma$, which characterize the growth of matter perturbations and the bending of light rays, respectively, depend on the sound speed $c_\psi$ of the longitudinal scalar perturbation associated with the vector field. Since $c_\psi$ is influenced by the transverse vector mode, the model exhibits sufficient flexibility to yield values of $\mu$ and $\Sigma$ close to 1. Consequently, unlike theories such as scalar Galileons, the present model can be consistent with observations of redshift-space distortions and integrated Sachs-Wolfe-galaxy cross-correlations.

Extremely metal-poor stars are intrinsically rare, but emerging methods exist to accurately classify them from all-sky Gaia XP low-resolution spectra. To assess their overall accuracy for targeting metal-poor stars, we present a high-resolution spectroscopic followup of 75 very metal-poor candidates selected from the catalog by R. Andrae, V. Chandra, and H. W. Rix. We discover 2 new extremely metal-poor ($\rm{[Fe / H]}<-3$) stars and 20 new very metal-poor ($\rm{[Fe/H]} < -2$) stars. Abundances of up to 22 elements are derived from 1D local thermodynamic equilibrium analysis and kinematic parameters are derived using Gaia astrometry and spectroscopic radial velocities. The chemodynamical properties are mostly consistent with expectations for halo stars, but we discover an Mg-enhanced CEMP star ($\mathrm{[Mg/Fe]} = 0.89$) and an Mg-poor star from an accreted ultra-faint dwarf galaxy. The Gaia XP metallicity estimates are consistent with our $\rm{[Fe/H]}$ measurements down to $\rm{[Fe/H]}\sim -3.0$, but estimates worsen in highly extincted regions. We find that 4 other XP-based metallicity catalogs succeed in mitigating contaminants and can also classify metal-poor stars robustly to $\rm{[Fe/H]}\sim -3.0$. Our results demonstrate the utility of Gaia XP spectra for identifying the most metal-poor stars across the Galaxy.

We determine the galaxy luminosity function of cluster galaxies in the nearby galaxy cluster Abell 2199 (A2199), focusing on the faint-end slope down to $M_r \sim -14.5$. To achieve this, we augment the existing dataset by adding redshift data from our deep MMT/Hectospec survey and from the Dark Energy Spectroscopic Instrument (DESI), significantly improving the spectroscopic completeness down to $r_{\mathrm{petro},0} = 20.8$ within the central $30^\prime$ region. The resulting luminosity function is well described by a Schechter function with a characteristic magnitude $M^* = -21.30 \pm 0.27$ and a faint-end slope $\alpha = -1.23 \pm 0.05$. This faint-end slope is consistent with those measured in the nearby Coma and Virgo clusters and in a cluster from the TNG50 cosmological simulation, and is slightly shallower than that of field galaxies. These findings indicate that the previously claimed steep faint-end upturn (with $\alpha \sim -2$) in nearby galaxy clusters is not supported. Instead, they indicate that environmental processes in dense cluster cores does not seem to trigger the formation or survival of low-mass galaxies, thereby preventing a steep faint-end upturn in the luminosity function.

We conduct a cross-correlation analysis between large-scale structures traced by the Two Micron All Sky Survey (2MASS) galaxy catalog and the cosmic dispersion measures of 133 localized fast radio bursts (FRBs). The cross-correlation signal is measured as a function of the comoving separation $R$ between 2MASS galaxies and background FRB sightlines, making full use of the available redshift information for both datasets. Our measurements are consistent with a null detection over the range $0.01 < R\, [h^{-1}\mathrm{Mpc}] < 1$. Using a halo-based model in which free-electron density profiles are drawn from the hydrodynamical simulation IllustrisTNG-300 (TNG300), we show that the null signal at $R \sim 0.01\, h^{-1}\mathrm{Mpc}$ is inconsistent with the TNG300 prediction. This discrepancy indicates that the hot-gas mass fraction in halos with masses of $10^{12-13}\, M_\odot$ hosting 2MASS galaxies must be lower than that predicted by TNG300. A simple phenomenological modification of the TNG300 model suggests that the hot-gas mass fraction in halos of $10^{12-13}\, M_\odot$ should be below $\sim 10\%$ of the global baryon fraction in the nearby universe, implying the need for stronger feedback in this mass range. Our constraints are consistent with those inferred from X-ray emission and Sunyaev-Zel'dovich measurements in galaxies, while providing a direct estimate of the hot-gas mass fraction that does not rely on electron-temperature measurements. These results demonstrate that galaxy-FRB cross correlations offer a powerful probe of feedback processes in galaxy formation.

The observed masses and radii of sub-Neptunes can be explained by a variety of bulk compositions, with the two leading scenarios being the gas dwarf and the water world scenario. The evolutionary history of sub-Neptunes on a population level has been proposed as a method to distinguish between the possible bulk compositions. Previous evolutionary models, however, neglected the crucial role of chemical interactions between the atmosphere and interior. We present a novel evolution framework for sub-Neptunes that not only considers the thermal evolution but also takes the chemical coupling of atmosphere and interior into account. Using this model, we examine how planets formed inside and outside the ice line can be observationally distinguished. Young sub-Neptunes store the majority of their volatile budget in the interior, independent of formation location and thus initial composition. Nevertheless, the atmospheric metallicity is a factor 4 higher for the planet formed outside the ice line. As the planet cools, hydrogen and oxygen exsolve from the interior, leading to an increase in atmosphere mass fraction for both planets, counteracting the contraction due to cooling. Consequently, radius evolution alone cannot distinguish sub-Neptunes formed inside the water ice line from water-rich planets formed outside of it. Instead, a key discriminator is the abundance of carbon-bearing species and the resulting atmospheric C/O ratio. For water-rich sub-Neptunes formed outside the \ice line, almost all carbon is in the gaseous phase. We find that high molar fractions of CH$_4$ ($\>10^{-2}$) and H$_2$O ($> 5\times10^{-2}$), and a high C/O ratio $(> 5\times10^{-1})$ are indicative of formation outside the ice line. In contrast, sub-Neptunes formed inside the ice line exhibit strongly suppressed CH$_4$ abundances, yielding C/O ratios ranging widely from $10^{-7}$ to $10^{-1}$. (Shortened version)

This work investigates the origin of the TeV emission detected by the Large High Altitude Air Shower Observatory (LHAASO) from NGC 4278, a galaxy hosting a low-luminosity active galactic nucleus (LLAGN). Considering two plausible scenarios, AGN jets and winds, we model the X-ray, GeV, and TeV emission during both TeV-low (quasi-quiet) and TeV-high (active) states. The spectral energy distributions can be explained either by single-zone leptonic emission from moderately relativistic jets or by lepto-hadronic emission from sub-relativistic winds. The best-fit parameters suggest that the transition from the quasi-quiet to the active state may be driven jointly by an enhanced accretion rate and the expansion of jets or winds. We further show that future MeV and very-high-energy $\gamma$-ray observations can discriminate between the {leptonic and lepto-hadronic scenarios}. Although the neutrino flux from NGC~1068 predicted by the wind model is too low to be detected with current neutrino observatories, a lepto-hadronic wind scenario can account for the PeV diffuse neutrino background when adopting a local LLAGN density corrected for the TeV duty cycle, $n_{\rm L,0}(\Delta T_{\rm TeV}/T) \sim 10^{-5}~\rm Mpc^{-3}$, as inferred from the LHAASO detection.

We present the cosmological analysis of the one-dimensional Lyman-$\alpha$ flux power spectrum from the first data release of the Dark Energy Spectroscopic Instrument (DESI). We capture the dependence of the signal on cosmology and intergalactic medium physics using an emulator trained on a cosmological suite of hydrodynamical simulations, and we correct its predictions for the impact of astrophysical contaminants and systematics, many of these not considered in previous analyses. We employ this framework to constrain the amplitude and logarithmic slope of the linear matter power spectrum at $k_\star=0.009\,\mathrm{km^{-1}s}$ and redshift $z=3$, obtaining $\Delta^2_\star=0.379\pm0.032$ and $n_\star=-2.309\pm0.019$. The robustness of these constraints is validated through the analysis of mocks and a large number of alternative data analysis variations, with cosmological parameters kept blinded throughout the validation process. We then combine our results with constraints from DESI BAO and temperature, polarization, and lensing measurements from Planck, ACT, and SPT-3G to set constraints on $\Lambda$CDM extensions. While our measurements do not significantly tighten the limits on the sum of neutrino masses from the combination of these probes, they sharpen the constraints on the effective number of relativistic species, $N_\mathrm{eff}=3.02\pm0.10$, the running of the spectral index, $\alpha_\mathrm{s}=0.0014\pm0.0041$, and the running of the running, $\beta_\mathrm{s}=-0.0006\pm0.0048$, by a factor of 1.18, 1.27, and 1.90, respectively. We conclude by outlining the improvements needed to fully reach the level of confidence implied by these uncertainties.

In this study, we present the chemical abundances of Fe, Mg, Al, Si, and K for a sample of 233 likely member stars of M 54. All the stars were observed with the FLAMES high-resolution multi-object spectrograph mounted at the VLT. Our analysis confirmed the presence of a large metallicity range in M 54, with the majority of the stars having -1.8 < [Fe/H] < -1.0 dex and few stars with [Fe/H] > -1.0 dex. The mean value of the total sample is [Fe/H] = -1.40 ({\sigma} = 0.22 dex). A Markov Chain Monte Carlo analysis revealed that the observed spread in [Fe/H] is compatible with a non-null intrinsic iron dispersion. We also found that the metallicity distribution function and the broadening of the red giant branch of M 54 are not compatible with a single age, but instead they suggest a wide age range from ~ 13 Gyr to ~ 1 - 2 Gyr or a smaller age range if a significant He enhancement (Y ~ 0.35/0.40) is present in the most metal-rich stars. We identified among the stars in M 54 the entire pattern of anticorrelations linked to the MgAl burning cycle. In particular, the metal-rich component displays a higher level of H-burning with the presence of more extended anticorrelations than the metal-poor component. No Mg-poor ([Mg/Fe]<0.0 dex) stars are identified in M 54. The evidence collected so far cannot be explained neither with a globular cluster-like scenario nor with a galactic chemical evolution. The chemical properties of M 54 can be explained within a scenario where this system formed through the merging of two globular clusters, the metal-poor one with standard characteristics and the more metal-rich one with more pronounced chemical anomalies, a possibly younger than the first one. M 54 is confirmed as a key stellar system for explaining the chemical evolution of a nuclear star cluster.

This work conducts a three-dimensional (3D), nonlinear magnetohydrodynamic (MHD) simulation to investigate wave generating, time-dependent reconnection around a magnetic null point. A non-periodic perturbation (in the $xz$-plane) triggers oscillatory reconnection (OR) at the 3D null, resulting in a self-sustained oscillation with a constant period $P$. We investigate the response of the system using three distinct wave proxies (compressible parallel, compressible transverse and incompressible parallel) as well as Spectral Proper Orthogonal Decomposition for decoupling and analyzing the resultant MHD wave behavior. We find that OR generates a slow magnetoacoustic wave of period $P$ that propagates outwards in all directions along the spine and fan plane of the 3D null point. We also find the generation of a propagating Alfvén wave of period $P$, exclusively along the $y$-axis in the fan plane, i.e. in the direction perpendicular to the spine motion. These findings provide new insights into waves generated from a 3D null point and their implications for coronal seismology.

We present an overview of the main characteristics of several spectroscopic surveys designed to advance our understanding of the physical properties and evolution of massive stars. We also summarize key results obtained from the analysis of these datasets, highlighting how the interpretation of some observables in the framework of massive-star evolution is considerably more complex than previously anticipated.

Recent reverberation delay measurements have moved beyond the 10 keV X-ray range, providing evidence for the Compton hump (a.k.a. reflection hump) in the lag spectra. We report the relativistic reverberation of the reflection hump in the bright Seyfert\,1 galaxy IC\,4329A based on a long {\it NuSTAR} observation. We find a delayed response of the 20--30 keV X-ray band, with a lag time of $\sim1825$ s at frequencies $0.5-1.5 \times 10^{-4}$ Hz. The lag amplitude drops to $\sim195$ s as the frequencies increase to $(1.5-10)\times10^{-4}$ Hz. Including IC\,4329A, so far five sources have been explored for reflection hump reverberation. We perform reverberation modelling of the 3--50 keV lag-energy spectra using the general relativistic transfer function code, which provides independent timing-based measurements of the black hole mass $M_{\rm BH}=1.37_{-0.36}^{+0.33}\times10^8~M_{\odot}$ and the coronal height $h=2.45_{-2.36}^{+1.92}~R_{\rm g}$ (with uncertainties at 90\% confidence). Within the uncertainties, the measured mass is found to be consistent with the previous finding. Furthermore, we undertake reflection spectroscopy to account for the hump feature and the associated relativistic effect using the time-averaged flux spectrum. Further sampling of the {\it NuSTAR} data (with a bin width of 0.2/0.4 keV below and above 10 keV) that reshapes the spectral resolution allows us to constrain the coronal temperature at $50.26_{-4.03}^{+5.58}$ keV -- consistent with the previous result from the combined {\it Suzaku} and {\it NuSTAR} data.

Twenty co-orbital asteroids of Venus are currently known, several of which may evolve into potentially hazardous asteroids (PHAs) over timescales of thousands of years. We report the identification and first dynamical characterization of 2025 TV$_{10}$, a newly discovered Venus co-orbital asteroid, and assess its potential collisional hazard to Earth. We performed numerical simulations of a large number of asteroid clones, and we studied their close encounters with Venus and Earth. The asteroid may leave its co-orbital orbit on timescales of 1200 yr. The orbit of 2025 TV$_{10}$ is one of the closest of the absolute minimum of the "minimum orbital intersection distance" (MOID) with Earth for the Venus co-orbital asteroids known to date. Owing to its orbital parameters, 2025 TV$_{10}$ represents one of the most dynamically hazardous members of the Venus co-orbital population identified to date. Its faint magnitude and restricted observability windows make future observations challenging but essential for constraining its orbit.

DESHIMA (Deep Spectroscopic HIgh-redshift MApper) is a broadband integrated superconducting spectrometer (ISS) for millimeter (mm) / sub-millimeter (sub-mm) wave astronomy based on Kinetic Inductance Detectors (KIDs). This paper describes characterization of DESHIMA 2.0 in laboratory settings. The instrument features NbTiN superconducting microstrip (MS) filters with low-loss a-SiC:H dielectric and an ultra-wideband leaky-wave antenna. A laboratory setup was designed, incorporating the cryostat housing cryogenic optics and ISS chip comprising 339 KIDs connected to MS filters tuned for (sub-)mm wave frequencies. Room-temperature mirrors on a hexapod stage allowed precise positioning and alignment of optical elements. The sky-position chopper was positioned on a motor-controlled stage for fine-tuned control over its position and alignment. Thanks to the multiplexing capability of KIDs, we could simultaneously measure multiple performance metrics across the entire frequency range. We showed that DESHIMA 2.0 achieved significant improvements in performance compared to its predecessor (DESHIMA 1.0): measured instantaneous frequency coverage was 200$-$400 GHz with a mean filter $Q_{filter}$ of $340 \pm 50$; instrument efficiency reached $\sim8$ \%, indicating 4 times wider band coverage and 4 times higher sensitivity. The yield rate for MS filters exceeded 98 \%. The estimated aperture efficiency from measured beam patterns agreed well with the designed value of approximately 70 \%. The telescope far-field beam patterns calculated from measured beam patterns also exhibited good agreement with design specifications. We also demonstrated validity of a new method of absolute frequency calibration using the data from beam pattern measurement.

The Spectro-Photometer for the History of the Universe, Epoch of Reionization and Ices Explorer (SPHEREx) will provide all-sky spectral survey data covering optical to mid-infrared wavelengths with a spatial resolution of 6\farcs2, which can be widely used to study galaxy formation and evolution. We investigate the galaxy-galaxy blending in SPHEREx datasets using the mock galaxy catalogs generated from cosmological simulations and observational data. Only $\sim0.7\%$ of the galaxies will be blended with other galaxies in all-sky survey data with a limiting magnitude of 19 AB mag. However, the fraction of blended galaxies dramatically increases to $\sim7$--$9\%$ in the deep survey area around the ecliptic poles, where the depth reaches $\sim22$ AB mag. We examine the impact of the blending in the number count and luminosity function analyses using the SPHEREx data. We find that the number count can be overestimated by up to $10$--$20\%$ in the deep regions due to the flux boosting, suggesting that the impact of galaxy-galaxy blending on the number count is moderate. However, galaxy-galaxy blending can marginally change the luminosity function by up to 50\%\ over a wide range of redshifts. As we only employ the magnitude limit at $K_s$-band for the source detection, the blending fractions determined in this study should be regarded as lower limits.

Eta-Virginids is a less-known meteor shower active in March. We investigated the meteoroids of this shower using fireball data from the European Fireball Network supplemented by video data of faint meteors. We first derived the criteria for assigning meteors to this shower. A fragmentation model was then applied to selected shower fireballs with good deceleration data and light curves. Meteoroid fragmentation strengths and bulk densities were derived and compared with three other showers. We have confirmed the four year periodicity in the activity of eta-Virginids and their presence in the 3:1 mean motion resonance with Jupiter. The orbital period of four years was directly measured for the fireballs. Fainter meteors showed somewhat longer periods but the shower is poor in faint meteors. No member fainter than magnitude +1 was observed instrumentally. The physical properties of meteoroids are different from cometary streams and are similar to the Geminids. The limiting fragmentation strength of 0.5 MPa and typical bulk density of cm-sized meteoroids of 1500 kg/m3 suggest that the parent body is a carbonaceous asteroid. Besides Geminids, eta-Virginids is another stream of asteroidal origin. Some small meteoroids have densities around 2500 kg/m3. Three asteroids of the sizes between 40 - 120 meters have been found to have similar orbits but their relation to eta-Virginids remains uncertain.

We investigate how supermassive black hole (SMBH) seeds form in the early Universe at the centers of ultralight dark matter (ULDM) halos. Focusing on the ULDM Jeans scale, we identify the critical conditions under which high-redshift baryonic gas, strongly confined by central solitonic cores of the halos, undergoes direct and monolithic collapse. The solitonic potential naturally drives rapid inflow and shock heating, allowing the gas to exceed the critical atomic-cooling and fragmentation-suppression threshold of $\sim 3 \times 10^4 {\rm K}$ without invoking an external UV background. We derive semi-analytic relations for the halo mass, soliton mass, baryonic core radius, and thermodynamic state of the gas, including the effects of baryonic contraction. These relations simultaneously determine the minimum and maximum SMBH seed masses as functions of redshift. In this framework, pristine gas clouds that satisfy the temperature threshold collapse without fragmentation, forming SMBH seeds with characteristic masses of order $\sim 10^5M_\odot$, while systems below the threshold are expected to form compact star clusters instead. Our model also implies an upper limit on the attainable SMBH mass, predicting a maximum mass scale of order $\sim10^{10}M_\odot$, consistent with the most massive quasars observed to date. The ULDM particle mass required to reproduce the inferred seed mass scale, $m \simeq 10^{-22}{\rm eV}$, coincides with the value favored by galactic-scale observations, providing a unified explanation for the characteristic masses of both galactic cores and early SMBH seeds. Our model predicts efficient SMBH seed formation at redshifts $z \gtrsim 10$ and offers a natural interpretation of recently observed little red dots as SMBHs embedded in compact, hot, ionized gas clouds.

In this paper we introduce PERTURB-c, a correlation-aware framework for interpreting black box regression models with one-dimensional structured inputs. We demonstrate this framework on a simulated case study with machine learning based transit spectroscopy retrievals of exoplanet WASP-107b. Characterising many exoplanet atmospheres can answer important questions about planetary populations, but traditional retrievals are very resource intensive; machine learning based methods offer a fast alternative however (i) they require high volumes data (only obtainable through simulations) to train and (ii) their complexity renders them black-boxes. Better understanding how they reach predictions can allow us to inspect for biases, which is especially important with simulated data, and verify that predictions are made on the basis of physically plausible features. This ultimately improves the ease of adoption of machine learning techniques. The most used methods to explain machine learning model predictions (such as SHAP and other methods that rely on stochastic sample generation) suffer from high computational complexity and struggle to account for interactions between inputs. PERTURB-c addresses these issues by leveraging physical knowledge of the known spectral correlation. For visualisation of this analysis, we propose a heat-map-based representation which is better suited to large numbers of input features along a single dimension, and that is more intuitive to those who are already familiar with retrieval methods. Note that while we chose this exoplanet retrieval context to demonstrate our methodologies, the PERTURB-c framework is model agnostic and in a broader context has potential value across a plethora of adjacent regression problems.

Extended gamma-ray sources surrounding middle-aged pulsars, primarily observed at teraelectronvolt energies, have been interpreted as pulsar halos, where relativistic $e^\pm$ diffuse into the interstellar medium and produce inverse-Compton (IC) emission. HESS J1813-126, associated with the energetic, radio-quiet gamma-ray pulsar PSR J1813-1246, has been suggested as a candidate pulsar halo, though its nature remains uncertain. We interpreted the high-energy emission of PSR J1813-1246 using the synchro-curvature (SC) radiation model and tested whether the gamma-ray spectral energy distribution (SED) of HESS J1813-126 can be explained as a pulsar halo powered by PSR J1813-1246. We explain the X-ray and gamma-ray SEDs of the pulsar using the SC framework. We further computed the transport and losses of $e^\pm$ injected by the pulsar through time-dependent diffusion-loss equations, exploring various common pulsar halo transport models. The resulting IC emission was compared with \textit{Fermi}-LAT, H.E.S.S., HAWC, and LHAASO data. We present predictions for the surface brightness profiles (SBPs) and the aperture-dependent emission for the different transport models, providing key diagnostics for assessing the observability of HESS J1813-126 with current and future instruments. The SC framework successfully reproduces the emission of PSR J1813-1246. The SED of HESS J1813-126 can be consistently reproduced within different pulsar halo frameworks, albeit with distinct predictions across different transport models. The corresponding SBP predictions and aperture-dependent emission offer testable signatures for future imaging atmospheric Cherenkov telescopes, which will be crucial for discriminating between the transport models. We further examined the link between the pulsar central engine and its extended halo by comparing the pair multiplicities in the magnetospheric and halo regions.

Cool, dense material is frequently observed to permeate the hot, tenuous solar corona in the form of prominences, spicules and coronal rain. Both the cool material and surrounding corona exist at temperatures that are effectively thermally stable, in that their local radiative losses occur on relatively long timescales compared to the dynamics. However, as the solar atmosphere evolves, driving mixing between the condensations and surrounding hot material, intermediate temperatures form, which can become subject to highly efficient radiative losses. The thermal energy lost due to radiation can far exceed the turbulent heating thus the system undergoes mixing-induced cooling. Here, a 3D radiative MHD simulation is performed of the shear-driven Kelvin-Helmholtz Instability (KHI) occurring between a cool condensation and the hot solar corona. During the evolution, thermal instabilities form naturally within the mixing layer, and grow with time to produce long, narrow structures that extend perpendicular to the magnetic field. The thermal instabilities form self-consistently within the mixing layer as small isolated events, and are then stretched by the background flows to create long structures in relatively narrow planes. The turbulent flows agitate the condensations and cause them to fragment, creating smaller localised clumps of cool, dense (prominence-like) material that can merge and further fragment. In the presented simulation, the thermal instabilities act to replenish the cool, dense material lost due to mixing, with the total mass of cool material being approximately constant through time. By analysing the thermal energy loss due to optically-thin radiation, thermal instabilities are found to account for 15-20\% of all radiative losses in the turbulent plasma.

We present a systematic pulse-by-pulse analysis of 22 long-duration GRBs observed by \emph{Swift}, each exhibiting a well-separated precursor before the main burst. We compare duration, spectral hardness ratio, minimum variability timescale (MVT), and spectral lag between these components. Both precursors and main bursts have durations and hardness broadly consistent with Type II GRBs. However, precursors show longer MVTs (by factors of 3-10) and diverse lags with near-zero median values, while main bursts display variable MVTs and positive lags. These differences suggest precursors may originate from distinct dissipation conditions, possibly due to cocoon shock breakout or early magnetically dominated outflows. Despite temporal differences, both episodes are consistent with a single collapsar origin, providing no evidence for dual-progenitor events. Our findings support pulse-resolved classification and show that precursors offer critical insights into jet formation and pre-burst activity.

Hot subdwarfs (hot sds) are compact, evolved stars near the Extreme Horizontal Branch (EHB) and are key to understanding stellar evolution and the ultraviolet excess in galaxies. We extend our previous analysis of Gaia XP spectra of hot subdwarf stars to a much larger sample, enabling a comprehensive study of their physical and binary properties. Our goal is to identify patterns in Gaia XP spectra, investigate binarity, and assess the influence of parameters such as temperature, helium abundance, and variability. We analyse approximately 20000 hot subdwarf candidates selected from the literature, combining Gaia XP data with published parameters. We apply Uniform Manifold Approximation and Projection (UMAP) to the XP coefficients, which represent the Gaia XP spectra in a compact, feature-based form, to construct a similarity map. We then use self-organizing maps (SOMs) and convolutional neural networks (CNNs) to classify spectra as binaries or singles, and as cool and helium-poor or hot and helium-rich. The spectra are normalised using asymmetric least squares baseline fitting to emphasise individual spectral features. The BP-RP colour dominates the similarity map, with additional influence from temperature, helium abundance, and variability. Most binaries, identified via the Virtual Observatory SED Analyser (VOSA), cluster in two filaments linked to main sequence companions. CNN classification suggests a strong correlation between variability and binarity, with binary fractions exceeding 60 percent for active hot subdwarfs. Gaia XP spectra combined with dimensionality reduction and machine learning effectively reveal patterns in hot subdwarf properties. Our findings indicate that binarity and environmental density strongly shape the evolutionary paths of hot subdwarfs, and we identify possible contamination by main sequence and cataclysmic variable stars in the base sample.

We set out to accurately measure gravity in the low-acceleration range $(10^{-11},10^{-9})$ m\,s$^{-2}$ from 3D motions of isolated wide binary stars. Gaia DR3 provides precise measurements of the four sky-plane components of the 3D relative displacement and velocity ($\mathbf{r}, \mathbf{v}$) for a wide binary, but not comparably precise line-of-sight (radial) separation and relative velocity $v_{r}$. Based on our new observations and the public databases/publications, we assemble a sample of 36 nearby (distance $<150$pc) wide binaries in the low-acceleration regime with accurate values of $v_{r}$ (uncertainty $< 100$ m\,s$^{-1}$). Kinematic contaminants such as undetected stellar companions are well under control using various observational diagnostics such as Gaia's ruwe parameter, the color-magnitude diagram, multi-epoch observations of radial velocities, Speckle interferometric follow-up observations, and requiring Hipparcos-Gaia proper motion consistency. For the parameter $\Gamma \equiv \log_{10}\sqrt{\gamma}$ with $\gamma \equiv G/G_{\rm N}$ (where $G$ is a parameter generalizing Newton's constant $G_{\rm N}$ in elliptical orbits), we find $\Gamma=0.102_{-0.021}^{+0.023}$, inconsistent with standard gravity at $4.9\sigma$, giving a gravity boost factor of $\gamma=1.600_{-0.141}^{+0.171}$. Four wide binaries have 3D relative velocities exceeding their estimated Newtonian escape velocities with $1<v_{\rm obs}/v_{\rm escN}\le1.2$. These systems are unlikely to be chance associations and are expected in a nonstandard paradigm such as Milgromian dynamics (MOND). The hypothesis that Newtonian gravity can be extrapolated to the low-acceleration limit is falsified by this independent study with accurate 3D velocities. Future radial velocity monitoring and Speckle interferometric imaging for larger samples will be useful to refine the present result.

Recent observational studies of red giant stars have estimated the visibility of their mixed oscillation modes, which is a proxy of the average energy of these modes. Among other things, they demonstrated that although the damping rate of the oscillations in the core of many red giants appears to be negligible, other red giants exhibit high core damping rates that are sometimes consistent with the infinite value limit. Up until now, it has not been possible to link the mixed mode visibilities to core damping rates in a quantitative way. In this study, we use the progressive wave picture to derive an analytical function expressing the approximate resonance pattern of red giants up to a proportionality factor. This function can model the influence of the damping on the oscillations, as well as take into account other effects such as mode asymmetries. In particular, this expression can be used to obtain a quantitative estimate for the visibility of mixed modes and to predict the detectability of mixed mode and multiplet signatures under different core damping rates. Here, we conduct a parameter study to investigate how the damping processes affect these aspects. We find that the visibility approaches the value expected for an infinite core damping rate already at finite values. Furthermore, we find that both the mixed mode and the multiplet signatures disappear at finite core damping rates. This implies that the observational characteristics of red giants with finite core damping rates can appear as if their core damping rate were infinite, providing an explanation for the observed populations. Moreover, we have used our method to quantitatively estimate the core damping rates of red giants with unusually low mixed mode amplitudes from their observed visibilities.

AGN variability offers a direct probe of accretion physics, disk structure, and black hole growth, but progress has been limited by sample size, cadence heterogeneity, and photometric systematics. The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will deliver multi-band light curves for millions of AGN, enabling variability studies at a true population scale. We synthesize recent results from the Zwicky Transient Facility (ZTF), which demonstrate that optical variability amplitudes and timescales are primarily regulated by accretion state, with secondary dependence on black hole mass and redshift, and establish the feasibility of survey-driven continuum reverberation mapping. ZTF measurements reveal optical continuum-emitting region sizes that often exceed standard thin disk predictions, implicating diffuse continuum emission from the broad line region as a significant contributor to observed inter-band lags. We evaluate the implications of LSST cadence and survey strategy, particularly the deep drilling fields, for continuum and emission line reverberation mapping, changing-look AGN, extreme variability quasars, and periodic variability searches. Key limitations of broadband photometric variability are identified, including variable emission line contamination, diffuse BLR continuum emission, and cadence-dependent lag recoverability. We argue that realizing LSST's full scientific potential requires community-scale, standardized variability metric pipelines, probabilistic classification integrated with alert brokers for follow-up triggering, and complementary medium-band photometric observations to isolate the accretion disk continuum. Together, these elements will enable LSST to convert photometric variability into quantitative constraints on accretion disks, BLR structure, and supermassive black hole growth across cosmic time.

The Hidden Gems project searches the TESS data for additional planets transiting low-mass stars in confirmed systems. Our goal is to identify planet candidates that are below the detection threshold set by the SPOC and QLP pipelines using sherlock, a specialized pipeline for robust detection and vetting of transit signals in TESS data. We present the discovery of two inner rocky planets in the TOI-237 and TOI-4336 A systems, confirmed with ground-based photometry from the TRAPPIST, SPECULOOS, ExTrA, and LCO facilities. TOI-237 c has a radius of 1.21 +/- 0.04Re, orbits its mid-M host star every 1.74 days, and is close to a 3:1 mean-motion resonance with TOI-237 b. TOI-4336 A c has a radius of 1.17 +/- 0.06Re, and orbits with a period of 7.58 days an M3.5 host star which is part of a hierarchical triple system. We performed model comparison to search for non-zero eccentricities, and found that the circular transit models are statistically favored. Dynamical simulations show that both systems are in stable configurations, and the TTVs expected for the TOI-237 system are of the order of seconds. TOI-237 c and TOI-4336 A c join the high-interest population of warm likely super-Earths below the so-called "radius valley". In particular, TOI-237 c is a good candidate for phase curve observations with JWST/MIRI thanks to the small radius of the host star and its short period.

The massive hot stars play crucial role in the dynamics of galaxies. These stars influence their surroundings through strong winds which are highly structured processes. The theoretical study of the non-symmetric phenomena of the stellar winds are becoming important these days, mainly because 1-D models are not sufficient enough. We present a new version of our Monte Carlo radiative transfer code Andy Antares with improved treatment of the velocity field for arbitrary geometries. Our aim is to develop a numerical scheme that can incorporate a general velocity field defined at discrete points. Our main objective is to calculate radiative transfer in a general input hydrodynamic model. The Andy Antares code currently calculates pure radiative transfer. The input model is pre-calculated by another hydrodynamical code. The whole radiative transfer calculation is then processed in a Cartesian grid. Radiative transfer is solved using the Monte Carlo approach in 3-D regardless of the input hydrodynamical model's dimension. The velocity field at any given point is interpolated using the trilinear interpolation. The optical depth is then integrated numerically along the photon's path. We verified the accuracy of the numerical velocity interpolation by comparison with results obtained for analytical velocity fields, achieving successful outcomes. We also tested the radiative transfer solution on a 3-D model generated from a 2-D hydrodynamic model, and obtained emergent radiation. The code is suitable for the numerical solution of radiative transfer in 3-D with arbitrary velocity fields.

The evolutionary origin of Wolf-Rayet (WR) stars at Solar metallicity is unclear. Single-star evolution from massive O stars, possibly via a Luminous Blue Variable phase, is challenged by binary period distributions of different WR subtypes. Wd1-72 is a WN7b+O binary embedded in the collective wind of the Galactic young massive cluster Westerlund 1 (Wd 1). It is surrounded by highly structured nebulosity, with cometary tails pointing away from Wd 1 and quasi-spherical droplets towards it. In this letter, we demonstrate that this morphology can be qualitatively reproduced by a hydrodynamic simulation of non-conservative Roche Lobe Overflow (RLOF) mass-loss into a cluster wind. Our model is based on a detailed binary evolution track consistent with key known properties of Wd1-72. Our work suggests Wd1-72 could be only ~10 kyr post-RLOF, and the hydrogen-free nature of Wd1-72 favours this being a second or subsequent RLOF episode. Follow-up observations could make Wd1-72 a valuable benchmark for probing mass-loss and mass-transfer in forming gravitational-wave binary-progenitor systems.

Gravitational lensing of gravitational wave (GW) will become the next frontier in studying cosmology and gravity. While time-delay cosmography using quadruply lensed GW events associated with optical images of the lens systems can provide precise measurement of the Hubble constant ($H_0$), they are considered to be much rarer than doubly lensed events. In this work, we analyze time-delay cosmography with doubly lensed GW events for the first time. We generate mock doubly lensed GW events with designed sensitivity of the LIGO-Virgo-KAGRA (LVK) O5 network, with LIGO post-O5 upgrade, and with Einstein Telescope (ET) + Cosmic Explorer (CE) respectively, and select the events that can be associated with future galaxy surveys. Over 1000 realizations, we find an average of 0.2(2.4) qualified events with the LVK O5(post-O5) network. Whereas with the ET+CE network, we find an average of 73.2 qualified events over 100 realizations. Using the Singular Isothermal Sphere (SIS) lens model, we jointly estimate waveform parameters and the impact parameter with doubly lensed GW signals, and then forecast the constraints on cosmological parameters and modified GW propagation by combining time-delay cosmography and the standard siren approach. The average posterior gives a constraint on $H_0$ with a relative uncertainty of $14\%$, $10\%$ and $0.42\%$ in the $\Lambda$CDM model for the LVK O5, LVK post-O5, and ET+CE network, respectively. While the LVK network gives uninformative constraints on the $(w_0,w_a)$ dynamical dark energy model, the ET+CE network yields a moderate constraint of $w_0=-1.02^{+0.31}_{-0.22}$ and $w_a=0.48^{+0.99}_{-1.54}$. In addition, our method can provide precise constraints on modified GW propagation effects jointly with $H_0$.

Particle-in-cell (PIC) numerical simulations are currently among the most advanced tools to investigate particle acceleration at relativistic shocks. Still, they come with limitations imposed by finite computing power, whose impact is not straightforward to evaluate a priori. Observational features are hence required as verification. energy electrons accelerated at external shocks, provides a testbed for such predictions. Current numerical studies suggest that in GRB afterglows the maximum synchrotron photon energy, which corresponds to the limit of electron acceleration, may fall within the $\sim$ 0.1--10 keV X-ray energy band at late times, $t\gtrsim 10^6 - 10^7$ s. To test this prediction, we analyzed the X-ray spectra of six GRBs with \emph{Swift}/XRT detections beyond $10^7$ s: our analysis reveals no clear evidence of a spectral cutoff. Using a model that accounts for the effect of the finite opening angle of the shock on the observed maximum synchrotron photon energy, we show that these observations are incompatible with PIC simulation predictions, unless one or more physical afterglow parameters attain values at odds with those typically inferred from afterglow modeling (small radiative efficiency, low ambient density, large equipartition fraction $\epsilon_{\rm B}$ of the magnetic field). These findings challenge existing numerical simulation results and imply a more efficient acceleration of electrons to high-energies than seen in PIC simulations, with important implications for our understanding of particle acceleration in relativistic shocks.

One of the major discoveries by the James Webb Space Telescope (JWST) is the identification of a large population of luminous galaxies at $z>10$, challenging theoretical models for early galaxy formation. The unexpectedly high number density of these systems has triggered intense debate about potential differences in the physical properties of galaxies at such extreme redshifts and those at lower redshift. However, progress has been limited by the lack of rest-frame optical diagnostics, which are critical for constraining the key properties. Here we present deep JWST/MIRI observations of a UV-luminous galaxy at $z=11.04$, CEERS2-588, only 400 Myr after the Big Bang. CEERS2-588 is detected in the MIRI F560W and F770W bands, while deep MIRI/MRS spectroscopy yields no detection of H$\alpha$ or [OIII]$\lambda5007$ line, revealing a prominent Balmer break detected for the first time at $z>10$. Spectral energy distribution (SED) fitting indicates an extended star formation history possibly reaching $z>15$, followed by rapid quenching within the recent $\sim10$ Myr, in stark contrast to other $z>10$ galaxies. The MIRI detections also significantly improve our stellar mass estimate to $\mathrm{log}(M_*/M_\odot)=9.1^{+0.1}_{-0.1}$, making CEERS2-588 the most massive galaxy securely confirmed at $z>10$. Remarkably, the inferred gas-phase metallicity is near solar, exceeding predictions from current theoretical models. These results suggest that efficient starbursts play a key role in producing the abundant luminous galaxy population in the early universe.

We present a statistical analysis of the extreme brightness of the fast radio burst FRB 20250316A, a luminous, apparently non-repeating event detected by CHIME/FRB. Employing a model-agnostic framework based on the Generalized Extreme Value (GEV) distribution applied to the second CHIME/FRB catalog, we quantify its rarity within the current population. Bayesian fitting of GEV models to block maxima of peak flux and fluence data reveals FRB 20250316A to be a profound statistical outlier. For the peak flux, the analysis yields a return period of $\sim$ 600 years ($1\sigma$ credible level), with the underlying distribution being of the heavy-tailed, unbounded Fréchet type ($\xi > 0$). The fluence analysis indicates greater complexity: while the full sample suggests a Fréchet-type distribution with a $\sim50$-year return period in $1\sigma$ credible level, the removal of three other notable outliers points toward a light-tailed Weibull-type distribution ($\xi < 0$) with a finite upper bound far exceeded by FRB 20250316A. This dichotomy highlights the challenge in characterizing the tail of the FRB luminosity function with limited data. Although less extreme in recurrence time than the ``Brightest Of All Time'' gamma-ray burst GRB 221009A, FRB 20250316A constitutes a similarly exceptional event (a potential FRB ``BOAT'') within the short observational history of wide-field radio surveys. Our results underscore the existence of rare, highly luminous events that may probe the upper limits or distinct sub-populations of the FRB luminosity distribution.

Nuclear star clusters (NSCs) are among the densest stellar systems in the Universe and often coexist with supermassive black holes (SMBHs) at galaxy centres. While SMBH formation histories are essentially lost, NSCs preserve evolutionary imprints through their stellar populations and stellar kinematics, reflecting the cumulative effects of mergers, accretion, and internal dynamical evolution. We aim to investigate the orbital structure of the unusually large NSC in FCC 47 (NGC 1336) by decomposing its stellar orbits into dynamically distinct components. We extract stellar kinematics, and in particular the line-of-sight velocity distributions (LOSVDs), from VLT/MUSE integral-field spectroscopy using the non-parametric Bayes-LOSVD approach, and apply triaxial Schwarzschild orbit-superposition modelling with the DYNAMITE software. We decompose the orbit library into hot, warm, cold, and counter-rotating components. We detect triple-peaked LOSVDs in the nucleus, indicating a complex orbital structure. The NSC forms a counter-rotating, kinematically decoupled component. A hot pressure-supported component, a warm counter-rotating structure and a counter-rotating cold disk in the centre suggest hierarchical assembly via early star cluster accretion and later in situ star formation. Our orbital decomposition of FCC 47 supports a hybrid formation scenario for this NSC. Dynamically distinct substructures reflect the interplay of accretion and in situ star formation during galaxy evolution.

Recent studies by B. Villarroel and colleagues have assembled and analyzed datasets of unidentified features measured from digital scans of photographic plates captured by the first-epoch Palomar Observatory Sky Survey (POSS1) in the pre-Sputnik era. These studies have called attention to (i) a purported deficit of features within Earth's shadow; (ii) the sporadic presence of linear clusters; and (iii) a positive correlation between the timing of feature observations and nuclear tests as well as Unidentified Aerial Phenomena (UAP) sighting reports. These observations were cited as evidence that some fraction of the unidentified features represent glinting artificial objects near Earth. We have examined these claims using two related, previously published datasets. When analyzing the most vetted of these, we do not observe the reported deficit in the terrestrial shadow. We determine that a third of the features in the reported linear clusters were not confidently distinguished from catalog stars. We find that the reported correlation between the timing of feature observations and nuclear tests becomes insignificant after properly normalizing by the number of observation days, and is almost completely determined by the observation schedule of the Palomar telescope. We uncover important inconsistencies in the definitions of the datasets used in these studies, as well as the use of unvalidated datasets containing catalog stars, scan artifacts, and plate defects. It has not been shown that any of the features in these datasets represent optical transients. We examine the spatial distribution of the plate-derived features, finding an overall gradual increase in number density toward the corners and edges of plates, as well as examples of (i) empty north-south strips that span multiple plates; (ii) clusters and voids having geometric shapes; and (iii) amorphous clusters.

We present a catalog of broad H$\beta$ variability properties for all spectra of quasars with $z<0.8$ and at least two observations included in the Sloan Digital Sky Survey (SDSS) Data Release 16 quasar catalog. For each spectrum, we perform a spectral decomposition to isolate the broad H$\beta$ emission. We measure the luminosity, FWHM, equivalent width, centroid, and Pearson skewness coefficient of broad H$\beta$ and provide derived physical properties such as the single-epoch black hole mass and the bolometric luminosity. For each pair of spectra in the sample, we calculate the change in radial velocity of the centroid of broad H$\beta$ emission ($\Delta v_{rad}$) as well as other derived properties related to broad H$\beta$ shape variability. We use forward-modeling methods to estimate the uncertainty in our measurements and discuss an improved method for estimating the uncertainty in $\Delta v_{rad}$ in the case where a spectral decomposition is used to isolate the broad H$\beta$ emission. We find that $\Delta v_{rad}$ is not normally distributed and that the shape of the distribution depends on the interval between observations. We discuss the effect of the predominance of the Reverberation Mapping subsample in the sample of pairs of spectra in SDSS.

We use deep Hubble Space Telescope imaging to study stellar populations in the outer halo of the spiral galaxy M96, located in the dynamically active Leo I galaxy group. Our imaging targets two fields at a projected distance of 50 kpc from the galaxy's center, with a 50% photometric completeness limit of F814W = 28.0, nearly two magnitudes below the tip of the red giant branch. In both fields we find a clear detection of red giant stars in M96's halo, with a space density that corresponds to an equivalent broadband surface brightness of $\mu_V \approx $ 31.7 mag arcsec$^{-2}$. We find little evidence for any difference in the spatial density or color of the RGB stars in the two fields. Using isochrone matching we derive a median metallicity for the red giants of [M/H] = -1.36 with an interquartile spread of $\pm$0.75 dex. Adopting a power-law radial density profile, we also derive a total halo mass of $M_h = 7.8^{+17.4}_{-4.9}\times10^9$ M$_\odot$, implying a stellar halo mass fraction of $M_{*,halo}/M_{*,tot} = 15^{+33}_{-9}$%, on the high end for spiral galaxies, but with significant uncertainty. Finally, we find that M96 appears offset from the stellar halo mass-metallicity relationship for spirals, with a halo that is distinctly metal-poor for its halo mass. While a variety of systematic effects could have conspired to drive M96 off this relationship, if confirmed our results may argue for a markedly different accretion history for M96 compared to other spirals in the nearby universe.

We study the effects of emission geometry on the capability to constrain black hole parameters from measurements of the shadow and inner-shadow of a Reissner-Nordström black hole. We investigate the capability to constrain mass, charge, observer inclination, and emission co-latitude from images of black hole accretion flows that would arise from thick and thin accretion disks. We confirm previous studies that have shown that independent radii measurements of the shadow and inner-shadow can constrain black hole parameters if the viewing inclination is known, but find that it is only possible if the true emission geometry is also assumed. We study the constraining capabilities of the shadow and inner-shadow observations of M87* and Sgr A* like systems within the context of the BHEX and NgEHT future observatories.

Several methods for identifying Be star candidates are reviewed for observational bias with respect to system inclination, that is the angle between the stellar/disk rotation axis and the observer's line of sight, with focus on two photometric methods that leverage narrow-band filters centred on H$\alpha$ and a spectroscopic method using a H$\alpha$ peak-finding algorithm. Tests for bias were performed using a sample of 20,000 synthetic Be stars drawn from a Salpeter initial mass function and computed libraries of spectral energy distributions and H$\alpha$ profiles. The spectroscopic method showed substantial bias against high inclinations ($i > 80^\circ$). Both photometric methods were biased against low inclinations, with one also biased against inclinations above $80^\circ$, resulting in a surplus in the Be star candidate detection rate for moderate inclinations ($ 50^\circ < i < 80^\circ$). Inclination probability distributions, including the random $\sin i$ factor, are given for the three methods that can be applied to observational samples.

We present the first spectroscopic characterization of an interstellar object at mid-infrared wavelengths. Post-perihelion observations of 3I/ATLAS using the JWST/MIRI medium-resolution spectrometer were obtained on 2025 December 15--16 and 27 when the object was at heliocentric distances of 2.20 and 2.54 au, respectively. Our 5--28 micron spectra exhibit fluorescence features from several gaseous species, including the $\nu_2$ band of water at 5.8--7.0 microns. the primary $\nu_2$ and associated hot bands of carbon dioxide around 15 microns, and a forbidden transition of atomic nickel at 7.507 microns. We also report the first direct detection of methane in an interstellar object. The delayed onset of methane production relative to water suggests past depletion from the outermost layers, with the observed methane emerging from unprocessed subsurface material. Comparison of the volatile production rates measured during the two epochs indicate a significant reduction in overall outgassing over 12 days, with the measured water activity level dropping more steeply than other species. As shown through near-nucleus coma mapping, 3I continues to display an extended source of water production from icy grains entrained within the coma. Our production rate measurements confirm that 3I exhibits a strongly enhanced CO$_2$:H$_2$O mixing ratio relative to typical solar system comets, as well as a somewhat enriched CH$_4$:H$_2$O value.

The origin of the Hubble tension remains one of the central open problems in modern cosmology, with competing explanations invoking either early-Universe physics, late-time modifications of cosmic expansion, or unresolved observational systematics. In this Letter we propose a new, purely geometric null test of the late-time expansion history that is exactly independent of the Hubble constant. By combining strong-lensing time-delay distances with gravitational-wave standard-siren luminosity distances, we construct a dimensionless ratio that depends only on the redshift dependence of the expansion rate and can be both predicted from early-Universe data and measured directly at late times, without relying on the cosmic distance ladder or the sound horizon. We show that the comparison between the early- and late-time determinations of this ratio provides a transparent consistency test of the standard cosmological expansion. When combined with an independent standard-siren measurement of $H_{0}$, this framework allows one to unambiguously distinguish between early- and late-time origins of the Hubble tension. With the forthcoming detection of lensed gravitational-wave standard sirens, the proposed test provides a timely and robust framework for probing this long-standing cosmological puzzle.

We present a study of extended Ly$\alpha$ emission around four quasars hosted by dusty starbursts, which are composite systems thought to represent a transitional stage in quasar evolution. To extract faint CGM emission in the presence of bright point sources, we introduce {\it CubeCarve}, a dual-channel deconvolution algorithm that separates unresolved quasar emission from spatially extended structure. This approach enables reliable recovery of \Lya\ emission projected onto the quasar position without introducing subtraction artifacts. Using {\it CubeCarve}, we find that the \Lya\ surface brightness profiles of these systems are, on average, fainter and shallower than those of quasars of similar bolometric luminosities. We also find that the total integrated \Lya\ luminosities of the nebulae are lower in systems whose host galaxies exhibit brighter far-infrared emission. These results suggest that the CGM conditions in composite systems differ from those in the broader quasar population. Our study highlights both the physical diversity of quasar CGM environments and the effectiveness of {\it CubeCarve} for recovering diffuse emission in modern IFU datasets.

BE Lyncis (BE Lyn) is a well-studied high-amplitude $\delta$ Scuti variable star (HADS). Recently, Niu et al. (2026) analyzed a 39-year baseline of times of maximum light of BE Lyn, reporting that it is the most eccentric binary known ($e \approx 0.9989$) and hosts the nearest black hole (BH) known to date. We analyze Hipparcos and Gaia astrometry of BE Lyn, predicting what the observed proper motion anomaly (PMA) over the 25 year baseline between the two missions would be were the companion really a $\gtrsim 17.5\,M_{\odot}$ BH. We find that the predicted PMA is at least an order of magnitude larger than the observed value of $\approx 1.7 \pm 0.8$ mas yr$^{-1}$, regardless of the assumed orientation of the orbit. We predict the expected Gaia DR3 RUWE for different orientations of the putative BH binary, finding that it ranges from $\approx 2.5$-$4.0$, much larger than the reported value of $1.073$. The observed value is instead consistent with a low-mass secondary or a single star. We find that BE Lyncis would have received a 7-parameter acceleration solution if it were a BH binary, in contradiction with its absence from the Gaia DR3 non-single star catalogs. Finally, we show that the reported orbit is impossible because the luminous star would overflow its Roche lobe at periastron, irrespective of inclination. We recommend caution in interpreting light-travel time effect (LTTE) models that require very high eccentricities, face-on inclinations, or large companion masses. The observed pulsation timing variations are most likely simply a result of red noise or pulsation phase evolution.

The sound speed of dark energy encodes fundamental information about the microphysics underlying cosmic acceleration, yet remains essentially unconstrained by existing observations. We demonstrate that a lunar-based laser interferometer, such as the proposed Laser Interferometer Lunar Antenna (LILA), can directly probe the sound speed of dark energy by measuring the real-time evolution of horizon-scale gravitational potentials. Operating in the ultra-low-frequency gravitational band inaccessible from Earth, LILA is sensitive to scalar metric perturbations sourced by dark energy dynamics. Using both fluid and effective field theory descriptions, we develop a complete framework linking dark energy sound speed to observable strain signatures. We construct a likelihood pipeline and Fisher forecasts, showing that LILA can either detect clustering dark energy or exclude broad classes of models with unprecedented sensitivity. This establishes lunar interferometry as a novel and powerful probe of the physics driving cosmic acceleration.

We present the characterization of two planetary systems orbiting the M dwarfs TOI-4336 A (M3.5V) and TOI-4342 (M0V), each hosting two transiting planets previously validated with TESS and ground-based observations. We refined the photometry of the TOI-4342 system using TESS and LCOGT data, and characterized the host stars with NIRPS and ESPRESSO spectroscopy. High-precision ESPRESSO radial velocities allowed us to constrain the planetary masses and investigate their potential compositions. The TOI-4336 A system is composed of a sub-Neptune with a period of 16.34 days, a radius of $2.14 \pm 0.08$ Re, and a mass of $3.33 \pm 0.36$ Me, along with an inner super-Earth on a 7.59-day orbit with a radius of $1.25 \pm 0.07$ Re and a mass of $1.55 \pm 0.13$ Me. The TOI-4342 system hosts two sub-Neptunes of similar sizes ($2.33 \pm 0.09$ Re and $2.35 \pm 0.09$ Re), with periods of 5.54 and 10.69 days. Their masses are measured to be $7.3 \pm 1.3$ Me and $4.8 \pm 1.4$ Me, respectively. The RVs also reveal a planet candidate around TOI-4342, likely non-transiting, with a period of 47.5 days and a minimum mass of $17.8 \pm 3.0$ Me. With precise radii and masses, we derived bulk densities and explored possible compositions. The TOI-4336 A sub-Neptune and super-Earth have densities of $1.87 \pm 0.30$ and $4.35 \pm 0.79$ g cm$^{-3}$, while the two similar-sized sub-Neptunes in TOI-4342 show distinct densities of $3.18 \pm 0.67$ and $2.01 \pm 0.63$ g cm$^{-3}$. All four planets are excellent targets for future atmospheric characterization with JWST, and their multi-planet nature makes them especially interesting for comparative planetology. Notably, TOI-4336 A b stands out as one of the best-known targets in its size and temperature regime, with a TSM of 138, comparable to benchmark planets such as K2-18 b and LHS 1140 b.

We compute deviations from observed gravitational wave signals, where the amplitude of the signal is unchanged. As an example, we consider the detectability of low lying dynamical tides in binary neutron star or neutron star black hole mergers. Tidal forces can excite oscillatory modes of one or both of the stars in the binary when the orbital frequency of the binary system sweeps through the resonant mode frequency dissipating energy into the vibrational mode. The orbital energy loss to the vibrational mode extracts energy from the orbital motion, advancing the time to merger. The inspiral then continues with an excess phase and a time advance. Both will cause a mismatch when fitting to a system that has not gone through the resonance. To resolve this effect, we compute the mismatch for current and planned detectors using both a quasi-analytical approach that relies on the computation of moment integrals and an optimized version of the standard numerical match function. We conclude that detectability can occur for time advances of the order of 1 ms with advanced LVK detectors for an excess energy-flux that is a few percent of the gravitational wave emission. Our results contrast with previous work, which model this effect solely as a phase shift of the waveform or by using the difference in the number of cycles induced by the resonant behavior. We show that tidal resonance effects primarily cause a time advance of the merger, rather than a phase difference, and that the single-frequency approximation commonly used in the literature significantly overestimates the detectability of this effect.

Exact, analytic, asymptotically flat rotating black-hole solutions are exceedingly rare, with only a handful of examples known. Using a Kerr-Schild ansatz, we derive a multitude of exact, analytic, asymptotically flat rotating black-hole solutions within a broad class of Generalized Proca theories. These black holes differ significantly from Kerr black holes, as they possess primary hair and are non-circular, thus breaking a symmetry that vacuum black holes exhibit in General Relativity.

General relativity has its successes at the local astrophysical level, however, it seems to be insufficient in describing the Universe at large scales. In this work we investigate how the most general field theories in the context of general relativity can accomodate a phantom-to-quintessence transition which may be essential element of realistic Dark Energy scenarios in the late Universe. As we demonstrate in a very detailed manner, this is impossible for a canonical and minimally coupled single scalar field theory, but it may be possible for ghost condensate theories like $k$-essence theories. We point out how the ghost instabilities may be eliminated, and we analyze the quantitative features of a $k$-essence theory that may realize a phantom-to-quintessence transition in the late Universe. We also qualitatively compare the difficulties and fine-tunings required for $k$-essence theories to realize a phantom-to-quintessence transition, and how such a transition is naturally realized in modified gravity, without unnecessary fine-tunings and ghost eliminations.

Phantom divide line transitions are not possible in the context of single scalar field scalar-tensor theories. In this article we study a combined framework of a tachyonic minimally coupled single scalar field theory in the presence of an $\mathcal{R}^2$ correction term and with a rescaled Einstein-Hilbert term of the form $\sim \lambda \frac{\mathcal{R}}{16\pi G}$. Such terms can be part of an $f(\mathcal{R})$ gravity which in the large curvature regime yields such correction terms effectively. Alternatively, such terms can simply be quantum corrections to the scalar field action. We aim to answer two questions, firstly if this framework can lead to phantom divide line transitions and secondly whether the resulting model can be compatible with the ACT data. The model we studied is an inverse square power-law model, well known from tachyon inflation models. As we show, the field equations can be cast in terms of the scalar field solely, however the resulting theory is distinct from a single scalar field theory, because the phantom divide line is crossed during inflation. Thus initially the tachyonic nature of the scalar field generates a phantom equation of state parameter, and during inflation the phantom divide line is crossed, with the effective equation of state parameter at the end of inflation being $w=-1/3$ which corresponds to the non-accelerating state of the Universe. The model is proved to be compatible with the ACT data, only when the gravity during inflation is stronger than Einstein-Hilbert gravity, with the effective gravitational constant during inflation being $\frac{G}{\lambda}$. The effective theory is valid only during inflation, thus Big-Bang nucleosynthesis is not affected by the rescaling of the Einstein-Hilbert gravity. The feature of a phantom crossing in $f(\mathcal{R},\phi)$ frameworks is new in the literature.

The origin of life is yet a compelling scientific mystery that has sometimes been attributed to high-pressure impacts by small solar system bodies such as comets, meteoroids, asteroids, and transitional objects. High-pressure torsion (HPT) is an innovative method with which to simulate the extreme conditions of astronomical impacts and offers insights relevant to prebiotic chemistry. In the present study, we investigated the polymerization and stability of adenosine monophosphate (AMP), a key precursor to ribonucleic acid (RNA), in dry and hydrated conditions (10 wt% water) under 6 GPa at ambient and boiling water temperatures. Comprehensive analyses with the use of X-ray diffraction, Raman spectroscopy, Fourier-transform infrared spectroscopy, nuclear magnetic resonance, scanning electron microscopy, and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry revealed no evidence of polymerization, while AMP partly transformed to other organic compounds such as nucleobase-derived fragments of adenine, phosphoribose fragments, dehydrated adenosine, protonated adenosine, and oxidized adenosine. The torque measurements during HPT further highlight the mechanical behavior of AMP under extreme conditions. These findings suggest that, while HPT under the conditions tested does not facilitate polymerization, the formation of various compounds from AMP confirms the significance of astronomical impacts on the prebiotic chemistry of RNA on early Earth. Keywords: Ribonucleic acid (RNA), Origin of life; Phase transformations; Chemical reactions, Small solar system bodies

The formation of primordial black holes (PBHs) generally requires large density perturbations, which is widely supported by researchers. This paper studies the local coupling properties of the Starobinsky potential and KKLT potential by introducing a linear Lorentzian-type coupling, which locally breaks the slow roll conditions. We found that both positive and negative coupling can form a considerable abundance of PBH. Additionally, we also studied the scalar-induced gravitational waves (SIGWs) generated by this model.

Gravitational wave astronomy typically relies on rigorous, computationally expensive Bayesian analyses. Several methods have been developed to perform rapid Bayesian inference, but they are not yet used to inform our full analyses. We present a novel approach for doing this whilst ensuring that the Bayesian prior remains independent of the data, providing a statistically rigorous way to leverage low-latency information to accelerate the final inference. By combining the fast constraints from the simple-pe algorithm with the nested sampling acceleration technique of posterior repartitioning, we demonstrate that our method can guide the nested sampler towards the most probable regions of parameter space more efficiently for signal-to-noise ratios (SNR) greater than 20, while mathematically guaranteeing that the final inference is identical to that of a standard, uninformed analysis. We validate the method through an injection study, demonstrating that it produces statistically robust and unbiased results, whilst providing speedups of up to $2.2\times$ for binaries with SNRs $< 150$. Importantly, we show that the performance gain provided by our method scales with SNR, establishing it as a powerful technique to mitigate the cost of analysing signals from current and future gravitational-wave observatories.

The unimodular version of the Kaluza-Klein theory is briefly recalled, and its projection on the $4$-dimensional spacetime is constructed. Imposing unimodularity condition on the $5$-dimensional Kaluza-Klein metric, det$g_{AB}=1$ is equivalent with introducing cosmological term in Einstein's equations in $4$ dimensions, and with scalar field of the Brans-Dicke type. Singularity-free cosmological solutions with scalar field and with matter sources are constructed, and their basic properties analyzed, along the results obtained in previous publications. In the present paper, attention is focussed on the perturbative analysis of cosmological solutions, providing a clue concerning their stability against small fluctuations.

An increasing number of surveys has been reporting large-scale peculiar motions with sizes and speeds in excess of those allowed by the concordance cosmological model. These are the so called bulk flows, the presence of which has come to be treated as a problem for the $\Lambda$CDM paradigm. However, the limits of the $\Lambda$CDM model are based on Newtonian studies, which predict the mediocre $v\propto t^{1/3}$ growth-rate for the peculiar-velocity field ($v$). Recently, a few fully relativistic treatments have appeared in the literature, arguing for a much stronger velocity growth that could explain the reported fast and deep bulk flows. What separates the Newtonian from the relativistic studies is the gravitational input of the peculiar flux, namely of the kinetic energy triggered by the moving matter. The latter has no direct gravitational contribution in Newtonian theory, but it does so in general relativity. This drastically changes the driving agent of the peculiar-velocity field and boosts its linear growth. The aim of this work is to directly compare the two treatments, as well as identify and discuss the reasons for their different results. In the process, we also demonstrate how one could recover the relativistic growth-rate from a Newtonian setup by selectively including certain (typically ignored) source-free terms into the Poisson equation. This way, we provide a unified covariant comparison of the Newtonian, the quasi-Newtonian and the fully relativistic studies.

Recent studies have shown that dark matter halos can support regular black holes or compact stars by assuming an anisotropic energy-momentum tensor. In this paper, we extend the analysis to the dark matter halo as an isotropic perfect fluid. By employing galactic dark matter profiles-specifically the Einasto and Dehnen models-as the mass-energy density source, we numerically solve the Einstein field equations and find a class of non-singular, horizonless compact star solutions. Moreover, these configurations remain stable against axial perturbations while satisfying the dominant energy condition.

We perform fully nonlinear, spherically symmetric numerical simulations of superhorizon false-vacuum-domain (FVD) collapse in a coupled gravity-scalar-fluid system to study primordial black hole (PBH) formation during delayed first-order phase transitions (FOPTs). Using adaptive mesh refinement to resolve the bubble wall, we identify three dynamical outcomes: type B (supercritical) PBHs with an interior baby universe and a bifurcating trapping horizon, type A (subcritical) PBHs with an apparent horizon formed by direct wall collapse, and dispersal with no PBH formation. To separate these three cases, we evaluate two commonly used PBH-formation criteria: the time scale ratio $t_\mathrm{H}/t_\mathrm{V}$ (horizon crossing time versus vacuum-energy domination time) and the local density contrast $\delta(t_\mathrm{H})$ at horizon crossing. For the parameter space explored, we find that $t_\mathrm{H}/t_\mathrm{V}$ is a more robust predictor of outcome: type B PBHs form when $t_\mathrm{H}/t_\mathrm{V} \gtrsim 1$ (critical range $\sim 1.1 - 1.6$ in our survey), type A PBHs arise when $t_\mathrm{H}/t_\mathrm{V}$ is below this threshold but remains above a lower bound (typical range $\sim 0.35 - 0.7$), and no-PBH dispersal occurs when $t_\mathrm{H}/t_\mathrm{V}$ falls below this lower bound. When a clear thin-wall FVD boundary exists, $\delta(t_\mathrm{H})$ can correspondingly distinguish different outcomes (roughly $\delta_c(t_\mathrm{H}) \sim 1 - 1.7$ for type B and $\delta_c(t_\mathrm{H}) \sim 0.35 - 0.5$ for type A), but is highly sensitive to wall structure and model details and thus less universal. These results offer new insights into the dynamics of FVD collapse, quantify practical PBH-formation thresholds, and pave the way for precise predictions of PBH abundance from delayed FOPTs.

We apply analytical models to study the property of neutron stars and dark stars. With the aim of exploring the global observable properties of those compact stars, we investigate the total masses and radii, the tidal deformabilities and especially the fundamental (f -) mode oscillations. While we choose two typical models in this work, this method applies to any analytical equations of state. By comparing with the multi-messenger observations, one can constrain the corresponding parameters in those models.

As the number of detected gravitational wave (GW) events increases with the improved sensitivity of the observatories, detecting strongly lensed pairs of events is becoming a real possibility. Identifying such lensed pairs, however, remains challenging due to the computational cost and/or the reliance on prior knowledge of source parameters in existing methods. This study investigates a novel approach, Optimal Cross-Correlation Analysis for Multiplets (OCCAM), applied to strain data from one or more detectors for Compact Binary Coalescence (CBC) events identified by GW searches, using an optimal, mildly model-dependent, low computation cost approach to identify strongly lensed candidates. This technique efficiently narrows the search space, allowing for more sensitive, but (much) higher latency, algorithms to refine the results further. We demonstrate that our method performs significantly better than other computationally inexpensive methods. In particular, we achieve 97 percent (80 percent) lensed event detection at a pairwise false positive probability of approximately 13 percent (7 percent) for a single detector with LIGO design sensitivity, assuming an SNR greater than or equal to 10 astrophysically motivated lensed and unlensed populations. Thus, this method, using a network of detectors and in conjunction with sky-localisation information, can enormously reduce the false positive probability, making it highly viable to efficiently and quickly search for lensing pairs among thousands of events, including the sub-threshold candidates.