2025-12-12 | | Total: 16
Current pulsar timing array (PTA) analyses do not take full advantage of pulsar distance information, thereby missing out on improved angular resolution and on a potential factor-of-two gain in detection sensitivity for individual gravitational-wave (GW) sources. In this work, we investigate the impact of precise pulsar distance measurements on angular resolution as an extension to previous work measuring the angular resolution of a dense, isotropic PTA [Jow et al., 2025]. We present a coherent map-making technique that utilizes precise pulsar distance measurements to reach the diffraction-limited resolution of an individual source: $δθ_{\mathrm{diff}} \sim (1/\mathrm{SNR})(λ_{\mathrm{GW}}/r) \approx 2~\mathrm{arcmin}$, where the SNR refers to the detection strength of the source. With this level of angular resolution, identifying an EM counterpart may become feasible, enabling multi-messenger follow-up. We show that for $\rm SNR=10$, which may be the current sensitivity level using a coherent analysis, the diffraction limit is reached with roughly 9 pulsars. Moreover, angular resolution scales sharply with the number of known pulsar distances as $\sim (1/\mathrm{SNR})^{N_{\mathrm{dist}}/2}$. Thus, each additional pulsar with high signal-to-noise timing and precise distance measurement can improve PTA resolution by an order of magnitude. The distance to the best-timed millisecond pulsar (PSR J0437$-$4715) is already constrained to sub-parsec levels. We argue, therefore, that a coherent analysis of PTA data, fully incorporating pulsar distance information, is timely.
We present FAST-MEPSA, an optimised version of the MEPSA algorithm developed to detect peaks in uniformly sampled time series affected by uncorrelated Gaussian noise. Although originally conceived for the analysis of gamma-ray burst (GRB) light curves (LCs), MEPSA can be readily applied to other transient phenomena. The algorithm scans the input data by applying a set of 39 predefined patterns across multiple timescales. While robust and effective, its computational cost becomes significant at large re-binning factors. To address this, FAST-MEPSA introduces a sparser offset-scanning strategy. In parallel, building on MEPSA's flexibility, we introduce a 40th pattern specifically designed to recover a class of elusive peaks that are typically sub-threshold and lie on the rising edge of broader structures - often missed by the original pattern set. Both versions of FAST-MEPSA - with 39 and 40 patterns - were validated on simulated GRB LCs. Compared to MEPSA, the new implementation achieves a speed-up of nearly a factor 400 at high re-binning factors, with only a minor (~4%) reduction in the number of detected peaks. It retains the same detection efficiency while significantly lowering the false positive rate of low significance. The inclusion of the new pattern increases the recovery of previously undetected and sub-threshold peaks. These improvements make FAST-MEPSA an effective tool for large-scale analyses where a robust trade-off between speed, efficiency, and reliability is essential. The adoption of 40 patterns instead of the classical 39 is advisable when an enhanced efficiency in detecting faint events is desired. The code is made publicly available.
The Near Ultraviolet Transient Explorer (NUTEx) is a CubeSat-based near-ultraviolet (NUV) imaging payload designed for transient sky surveys and is currently under development. CubeSats are compact and cost-effective satellite platforms that have emerged as versatile tools for scientific exploration and technology demonstrations in space. NUTEx is an imaging telescope operating in the 200-300 nm wavelength range, intended for deployment on a micro-satellite bus. The optical system is based on a Ritchey Chretien (RC) telescope configuration, featuring a 146 mm primary mirror. The detector is a photon-counting microchannel plate (MCP) device with a solar-blind photocathode, paired with an in-house developed readout unit. The instrument has a wide field of view (FoV) of 4 deg, a peak effective area of approximately 18 sq cm at 260 nm, and can reach a sensitivity of 21 AB magnitude (SNR = 5) in a 200 second exposure. The primary scientific objective of NUTEx is to monitor the night sky for transient phenomena, such as supernova remnants, flaring M-dwarf stars, and other short-timescale events. The payload is currently scheduled for launch in Q2 2026. This paper presents the NUTEx instrument design, outlines its scientific goals and capabilities, and provides an overview of the electronics and mechanical subsystems, including structural analysis.
Stellar atmospheric parameters and elemental abundances are traditionally determined using template matching techniques based on high-resolution spectra. However, these methods are sensitive to noise and unsuitable for ultra-low-resolution data. Given that the Chinese Space Station Telescope (CSST) will acquire large volumes of ultra-low-resolution spectra, developing effective methods for ultra-low-resolution spectral analysis is crucial. In this work, we investigated the Fully Connected Residual Network (FCResNet) for simultaneously estimating atmospheric parameters ($T_\text{eff}$, $\log g$, [Fe/H]) and elemental abundances ([C/Fe], [N/Fe], [Mg/Fe]). We trained and evaluated FCResNet using CSST-like spectra (\textit{R} $\sim$ 200) generated by degrading LAMOST spectra (\textit{R} $\sim$ 1,800), with reference labels from APOGEE. FCResNet significantly outperforms traditional machine learning methods (KNN, XGBoost, SVR) and CNN in prediction precision. For spectra with g-band signal-to-noise ratio greater than 20, FCResNet achieves precisions of 78 K, 0.15 dex, 0.08 dex, 0.05 dex, 0.10 dex, and 0.05 dex for $T_\text{eff}$, $\log g$, [Fe/H], [C/Fe], [N/Fe] and [Mg/Fe], respectively, on the test set. FCResNet processes one million spectra in only 42 seconds while maintaining a simple architecture with just 348 KB model size. These results suggest that FCResNet is a practical and promising tool for processing the large volume of ultra-low-resolution spectra that will be obtained by CSST in the future.
We introduce the Morph approximation, a class of product approximations of probability densities that selects low-order disjoint parameter blocks by maximizing the sum of their total correlations. We use the posterior approximation via Morph as the importance distribution in optimal bridge sampling. We denote this procedure by MorphZ, which serves as a post-processing estimator of the marginal likelihood. The MorphZ estimator requires only posterior samples together with the prior and likelihood, and is fully agnostic to the choice of sampler. We evaluate MorphZ's performance across statistical benchmarks, pulsar timing array (PTA) models, compact binary coalescence (CBC) gravitational-wave (GW) simulations and the GW150914 event. Across these applications, spanning low to high dimensionalities, MorphZ yields accurate evidence at substantially reduced computational cost relative to standard approaches, and can improve these estimates even when posterior coverage is incomplete. Its bridge sampling relative error diagnostic provides conservative uncertainty estimates. Because MorphZ operates directly on posterior draws, it complements exploration-oriented samplers by enabling fast and reliable evidence estimation, while it can be seamlessly integrated into existing inference workflows.
A plethora of astronomical science cases can only be achieved with high angular resolution observations, and we can expect the number of these to grow as astronomers are constrained by the size limitations of single-aperture space telescopes, making space-based interferometry inevitable. However, the enabling technologies do not have flight heritage at the system level, and the concept remains immature to much of the community, meaning no direct-detection synthetic-aperture space-based interferometer has yet flown and an opportunity exists for the UK to take a world leading role. Here we propose the SPace-based InterFerometry Feasibility (SPIFF) Project as a program to address both issues simultaneously by: 1) completing a thorough survey of the science cases across the EM spectrum that would directly benefit from, or be impossible without, space-based interferometry; 2) down selecting key requirements via a Science Traceability Matrix mapping exercise; 3) host a workshop for the UK astronomical community to consolidate these findings; 4) build a technology demonstration mission to raise TRL and achieve flight heritage of critical technologies. Such a program positions the UK as the partner of choice for any future NASA or ESA space-based interferometry mission, allowing the UK to lead groundbreaking scientific discoveries, while also directly benefiting the UK industrial base by advancing domestic exportable technologies and building direct synergy with other UK space priorities. Indeed, the UK is uniquely positioned to lead in space-based interferometry, possessing a rare trifecta of domestic strengths: world-class expertise in ground-based interferometry and space-based instrumentation; commercial entities developing mission-critical technologies; and scientists whose research spans the full range of science cases that would benefit directly from space-based interferometry.
We propose Roman Coronagraph project HLC/575 nm observations of a newly-discovered brown dwarf (HIP 71618 B) from the Subaru/OASIS survey of young accelerating stars, which is supported by NASA headquarters with the directive to identify targets for the Roman Coronagraph that could fulfill TTR5 requirements and be observed during the technology demonstration phase. The target and multiple bright PSF references are within/close to the Roman Continuous Viewing Zone. A high SNR detection of this companion would singlehandedly fulfill TTR5 and would be the first optical detection of a companion at $<$10$^{-6}$ contrast. Roman CPP reference star vetting prioritizing stars that can be paired with HIP 71618 would aid the execution of a successful technology demonstration. Additional similar targets may be discovered from OASIS over the next few years that could increase CGI scheduling flexibility and enhance its scientific and technical return. A close collaborative partnership with the CPP team could ensure that they are schedulable.
Recent baryon acoustic oscillation measurements from the DESI, when combined with CMB data and Type Ia supernovae observations, indicate a preference for dynamical dark energy when considering the Chevallier-Polarski-Linder (CPL) model, over the standard ΛCDM or the wCDM model. However, the Bayes factor, a key metric for model comparison, remains inconclusive on which model is preferred. This paper applies the BF method, that integrates both Bayesian and frequentist approaches to DESI data to address the limitations of purely frequentist or Bayesian methods. It consists in considering the Bayes factor as a random variable and calculates its distribution, that results from values computed in a frequentist approach after perturbing the data following the model considered. We apply this hybrid method to DESI data, comparing the CPL and w models under various prior conditions, including weak and strong priors, and theory-informed priors. We find that, when the traditional bayes factor is considered, that weak priors favor the w model over CPL, while strong priors favor CPL. Additionally, theory-informed priors further enhance the preference for the w model. While when we apply the BF method, the preference for CPL over w is seen in all cases albeit with similar but reduced impact on the p-value from the different prior considerations. We also tried to generalize further, by perturbing as well the covariance matrix following the model considered, and found that, in general, the current data in that case is not stringent enough to disentangle between the two models. Our results demonstrate that varying the Bayes factor as a random variable, providing that the covariance matrix is kept as model independent, provides a robust model comparison, reducing the impact of prior dependence as well as offering quantitative assessment of the preferences of the competing models.(abridged)
The current fleet of X-ray telescopes produces a wealth of multi-dimensional data, allowing us to study sources in time, photon energy and polarization. At the same time, it has become increasingly clear that progress in our physical understanding will only come from studying these sources in multiple dimensions simultaneously. Enabling multi-dimensional studies of X-ray sources requires new theoretical models predicting these data sets, new methods to analyse them and a software framework to combine data, models and methods efficiently. In this paper, we introduce the alpha release of nDspec, a new python-based library designed to allow users to model one- and multi-dimensional datasets common to X-ray astronomy. In the alpha release, we focus on modelling time-averaged data as well as Fourier spectral-timing mode, but highlight how additional dimensions can be added. We discuss design philosophy and current features, and showcase an example use case by characterizing a NICER observation of a black hole X-ray binary. We also highlight current plans for extensions to other dimensions and new features.
Observing gravitational waves from sub-solar-mass, inspiraling compact binaries would provide almost smoking-gun evidence for primordial black holes. Here, we develop a method to search for ultra-compact binaries with chirp masses ranging from $[10^{-2},10^{-1}]M_\odot$. This mass range represents a previously unexplored gap in gravitational-wave searches for compact binaries: it was thought that the signals would too long for matched-filtering analyses but too short for time-frequency pattern-recognition techniques. Despite this, we show that a pattern-recognition technique, the Generalized frequency-Hough (GFH), can be employed with particular modifications that allow us to handle rapidly spinning-up binaries and to increase the statistical robustness of our method, and call this improved method BinaryGFH-v2. We then design a hypothetical search for binaries in this mass regime, compare the empirical and theoretical sensitivities of this method, and project constraints on formation rate densities and the fraction of dark matter that primordial black holes could compose in both current- and future-generation gravitational-wave detectors. Our results show that our method can be used to search for sub-solar-mass, ultra-compact objects in a mass regime that remains to-date unconstrained with gravitational waves.
The Babcock-Leighton dynamo, which relies on the generation of a poloidal field through the decay and dispersal of tilted bipolar magnetic regions (BMRs), is a promising paradigm for explaining the features of the solar magnetic cycle. In rapidly rotating stars, BMRs are expected to emerge at high latitudes, which are less efficient in generating the poloidal field due to poor cross-equatorial cancellation. The operation of the Babcock-Leighton dynamo in rapidly rotating stars is therefore questionable. We, for the first time, using a 3D kinematic dynamo model, STABLE, explore this question. By taking large-scale flows from mean-field hydrodynamics models for stars rotating at different speeds, We conduct a series of dynamo simulations in rapidly rotating stars, exploring the following four cases of spot deposition, each based on a different assumption about toroidal flux tube rise: (i) radial rise, (ii) parallel rise to the rotation axis, (iii) parallel rise combined with an increase in Joy's law slope with the stellar rotation rate, and (iv) increasing time delay and spot size. We find cyclic magnetic fields in all cases except case IV of the 1-day rotating star, for which the magnetic field is irregular. For the parallel-rise cases, the magnetic field becomes quadrupolar, and for all other cases, it is dipolar. Our work demonstrates that the Babcock-Leighton dynamo may operate even in rapidly rotating stars with starspots appearing at higher latitudes.
Precise stellar parameters for M-type stars, the Galaxy's most common stellar type, are crucial for numerous studies. In this work, we refined the LAMOST DR10 M-type star catalog through a two-stage process. First, we purified the catalog using techniques including deep learning and color-magnitude diagrams to remove 22,496 non-M spectra, correct 2,078 dwarf/giant classifications, and update 12,900 radial velocities. This resulted in a cleaner catalog containing 870,518 M-type spectra (820,493 dwarfs, 50,025 giants). Second, applying a label transfer strategy using values from APOGEE DR16 for parameter prediction with a ten-fold cross-validated CNN ensemble architecture, we predicted $T_\text{eff}$, $\log g$, [M/H], and [$α$/M] separately for M dwarfs and giants. The average internal errors for M dwarfs/giants are respectively: $T_\text{eff}$ 30/17 K, log $g$ 0.07/0.07 dex, [M/H] 0.07/0.05 dex, and [$α$/M] 0.02/0.02 dex. Comparison with APOGEE demonstrates external precisions of 34/14 K, 0.12/0.07 dex, 0.09/0.04 dex, and 0.03/0.02 dex for M dwarfs/giants, which represents precision improvements of over 20\% for M dwarfs and over 50\% for M giants compared to previous literature results. The catalog is available at https://nadc.china-vo.org/res/r101668/.
Gravity modes in main-sequence stars have traditionally been studied using a prewhitening approach, which iteratively identifies modes in the Fourier domain and subsequently tunes their frequencies, amplitudes, and phases through time-domain regression. While effective, this method becomes inefficient when analysing large volumes of long time-series data and often relies on subjective stopping criteria to determine the number of iterations. We aim to perform frequency extraction of gravity modes in slowly pulsating B (SPB) stars using a statistically robust, data-driven approach based on advanced power spectrum and harmonic analysis techniques. Our approach employs the multitaper non-uniform fast Fourier transform, mtNUFFT, a power spectrum estimator that addresses several statistical limitations of traditional methods such as the Lomb-Scargle periodogram. We apply its extension, the multitaper F-test, to extract coherent gravity modes from 4-year Kepler light curves of SPB stars and to search for period spacing patterns among the extracted modes. The multitaper F-test enables fast and accurate extraction of the properties of gravity modes with quasi-infinite lifetimes, preferentially selecting modes that exhibit purely periodic behaviour. Although the method typically extracts fewer frequencies than conventional prewhitening, it recovers most known modes and, in some cases, reveals new ones. We also find evidence for gravity modes with long but finite lifetimes, and detect more than one period spacing pattern in some of the studied SPB stars. Overall, the multitaper F-test offers a more objective and statistically sound alternative to prewhitening. It scales efficiently to large datasets containing thousands of pulsators, and has the potential to facilitate mode identification and to distinguish between the different excitation mechanisms operating in SPB stars.
Space-based time-domain telescopes such as CoRoT, Kepler/K2 and TESS have profoundly impacted astrophysics over the past two decades. Continuous light curves with high cadence and high photometric precision are now available for millions of sources within our galaxy and beyond. In addition to revolutionizing exoplanet science, the data have enabled breakthroughs ranging from the solar system to stellar interiors, the transient universe, and active galaxies. The key summary points of this review are: (1) Stellar astrophysics has been transformed by the ability to probe the internal structures of stars, test the physics of stellar convection, connect stellar rotation and magnetic activity, and reveal complex variability in young stars. (2) Ages of stellar populations probe the formation history of our Milky Way, and binary star variability enables the detection of "dark" galactic populations such as solar-mass black holes and neutron stars. (3) Early-time observations of explosive transients provide new insights into the progenitors of supernovae, while the quasi-periodic variability of galaxies probes the physics of accretion processes onto supermassive black holes and the tidal disruption of stars. (4) Observations of solar system objects reveal asteroid compositions through their rotation periods and amplitudes, constrain the cloud structure of ice giants, and allow the discovery of new objects in the outer solar system. (5) Open data policies and software have contributed to remarkable scientific productivity and enabled discoveries by citizen scientists, including new exoplanets and exotic variability in mature Sun-like stars.
Estimating the auto power spectrum of cosmological tracers from line-intensity mapping (LIM) data is often limited by instrumental noise, residual foregrounds, and systematics. Cross-power spectra between multiple lines offer a robust alternative, mitigating noise bias and systematics. However, inferring the auto spectrum from cross-correlations relies on two key assumptions: that all tracers are linearly biased with respect to the matter density field, and that they are strongly mutually correlated. In this work, we introduce a new diagnostic statistic, \(\mathcal{Q}\), which serves as a data-driven null test of these assumptions. Constructed from combinations of cross-spectra between four distinct spectral lines, \(\mathcal{Q}\) identifies regimes where cross-spectrum-based auto-spectrum reconstruction is unbiased. We validate its behavior using both analytic toy models and simulations of LIM observables, including star formation lines ([CII], [NII], [CI],[OIII]) and the 21-cm signal. We explore a range of redshifts and instrumental configurations, incorporating noise from representative surveys. Our results demonstrate that the criterion \( \mathcal{Q} \approx 1 \) reliably selects the modes where cross-spectrum estimators are valid, while significant deviations are an indicator that the key assumptions have been violated. The \( \mathcal{Q} \) diagnostic thus provides a simple yet powerful data-driven consistency check for multi-tracer LIM analyses.
The history of astronomical discovery shows that many of the most detectable phenomena, especially detection firsts, are not typical members of their broader class, but rather rare, extreme cases with disproportionately large observational signatures. Motivated by this, we propose the Eschatian Hypothesis: that the first confirmed detection of an extraterrestrial technological civilization is most likely to be an atypical example, one that is unusually "loud" (i.e., producing an anomalously strong technosignature), and plausibly in a transitory, unstable, or even terminal phase. Using a toy model, we derive conditions under which such loud civilizations dominate detections, finding for example that if a society is loud for only $10^{-6}$ of its lifetime, it must emit ${\gtrsim}1$% of its total observable energy budget during that phase to outrun quieter populations. The hypothesis naturally motivates agnostic anomaly searches in wide-field, multi-channel, continuous surveys as a practical strategy for a first detection of extraterrestrial technology.