2025-06-06 | | Total: 92
Multiple clumpy wind components (vout∼0.2−0.3c) in the luminous quasar PDS 456 have recently been resolved by XRISM in the Fe-K band for the first time. In this paper, we investigate the structure of ultra-fast outflows (UFOs) using coordinated observations from XRISM, XMM-Newton, and NuSTAR, along with the self-consistently calculated photoionization model \texttt{PION}. Our results reveal a stratified ionization structure likely driven by the radiation field, characterized by a relation between wind velocity and ionization parameter vout∝ξ(0.38±0.06). To evaluate the impact of the screening effect, we tested all possible order permutations of six \texttt{PION} components. We find that highly ionized UFOs (logξ>4.5) are insensitive to their relative positions, whereas the soft X-ray UFO (logξ∼3 and vout∼0.27c) and the lowest-ionized hard X-ray UFO (logξ∼4.1 and vout∼0.23c) are statistically favored -- based on the evidence from both the C-statistic and Bayesian analysis -- to occupy the middle and innermost layers, respectively. This suggests a possible trend where slower UFOs are launched from regions closer to the supermassive black hole (SMBH). The soft X-ray UFO is found to be thermally unstable, regardless of its relative position. However, its location remains unclear. Our sequence analysis and its similarity to hard X-ray UFOs suggest that they may be co-spatial, while variability constraints support its location within the broad-line region at sub-parsec scales. Simulations with the gate-valve opened XRISM show that high-resolution soft X-ray data can enhance the reliability of our results. Furthermore, simulations with the future X-ray mission NewAthena demonstrate its capability to resolve the absorber sequence and spatial distributions, enabling the determination of UFO structures and their roles in AGN feedback.
The likely JWST detection of vibrationally excited H3+ emission in the irradiated disk d203-506, reported by Schroetter et al., raises the question of whether cosmic-ray ionization is enhanced in disks within clustered star-forming regions, or if alternative mechanisms contribute to H3+ formation and excitation. We present a detailed model of the photodissociation region (PDR) component of a protoplanetary disk-comprising the outer disk surface and the photoevaporative wind-exposed to strong external far-ultraviolet (FUV) radiation. We investigate key gas-phase reactions involving excited H2 that lead to the formation of H3+ in the PDR, including detailed state-to-state dynamical calculations of reactions H2 + HOC+ -> H3+ + CO and H2 + H+ -> H2+ + H. We also consider the effects of photoionization of vibrationally excited H2(v>=4), a process that has not previously been included in disk models. We find that these FUV-driven reactions dominate the formation of H3+ in the PDR of strongly irradiated disks, largely independently of cosmic-ray ionization. The predicted H3+ abundance peaks at x(H3+)~10^-8, coinciding with regions where HOC+ is also abundant. The high abundances of H3+ and HOC+ are ultimately linked to the strength of the external FUV field (G0), the presence of C+, and an enhanced reservoir of hot H2O. The predicted H3+ column densities (~10^13 cm^-2) are consistent with the presence of H3+ in the PDR of irradiated disks. We also find that formation pumping, resulting from exoergic reactions between excited H2 and HOC+, drive the vibrational excitation of H3+. We expect this photochemistry to be highly active in disks where G0 > 10^3. The H3+ formation pathways studied here may also be relevant in the inner disk region (near the host star), in exoplanetary ionospheres, and in the early Universe.
Most protoplanetary disks experience a phase in which they are subjected to strong ultraviolet radiation from nearby massive stars. This UV radiation can substantially alter their chemistry by producing numerous radicals and molecular ions. In this Letter we present detailed analysis of the JWST-NIRSpec spectrum of the d203-506 obtained as part of the PDRs4All Early Release Science program. Using state-of-the-art spectroscopic data, we searched for species using a multi-molecule fitting tool, PAHTATmol, that we developed for this purpose. Based on this analysis, we report the clear detection of ro-vibrational emission of the CH radical and likely detection of the H+3 molecular ion, with estimated abundances of a few times 10−7 and approximately 10−8, respectively. The presence of CH is predicted by gas-phase models and well explained by hydrocarbon photochemistry. H+3 is usually formed through reactions of H2 with H+2 originating from cosmic ray ionization of H2. However, recent theoretical studies suggest that H+3 also forms through UV-driven chemistry in strongly irradiated (G0>103), dense (nH>106 cm−3) gas. The latter is favored as an explanation for the presence of ``hot'' H+3 (Tex≳1000 K) in the outer disk layers of d203-506, coinciding with the emission of FUV-pumped H_2 and other ``PDR species'', such as CH^+, CH_3^+, and OH. Our detection of infrared emission from vibrationally excited H_3^+ and CH raises questions about their excitation mechanisms and, underscore that UV radiation can have a profound impact on the chemistry of planet forming disks. They also demonstrate the power of JWST pushing the limit for the detection of elusive species in protoplanetary disks.
We present a simulation of the time-domain catalog for the Nancy Grace Roman Space Telescope's High-Latitude Time-Domain Core Community Survey. This simulation, called the Hourglass simulation, uses the most up-to-date spectral energy distribution models and rate measurements for ten extra-galactic time-domain sources. We simulate these models through the design reference Roman Space Telescope survey: four filters per tier, a five day cadence, over two years, a wide tier of 19 deg^2 and a deep tier of 4.2 deg^2, with \sim20% of those areas also covered with prism observations. We find that a science-independent Roman time-domain catalog, assuming a S/N at max of >5, would have approximately 21,000 Type Ia supernovae, 40,000 core-collapse supernovae, around 70 superluminous supernovae, \sim35 tidal disruption events, 3 kilonovae, and possibly pair-instability supernovae. In total, Hourglass has over 64,000 transient objects, 11 million photometric observations, and 500,000 spectra. Additionally, Hourglass is a useful data set to train machine learning classification algorithms. We show that SCONE is able to photometrically classify Type Ia supernovae with high precision (\sim95%) to a z > 2. Finally, we present the first realistic simulations of non-Type Ia supernovae spectral-time series data from Roman's prism.
Binary systems in the Asymptotic Giant Branch (AGB) phase are widely recognized as a leading theoretical framework underpinning the observed asymmetric morphologies of planetary nebulae. However, the detection of binary companions in AGB systems is severely hampered by the overwhelming brightness and variability of the evolved primary star, which dominate the photo-metric and spectroscopic signatures. Ultraviolet (UV) excess emission has been proposed as a candidate diagnostic for the presence of binary companions in AGB systems. This paper evaluates the Chinese Space Station Telescope's (CSST) ability to detect UV excess emission in AGB stars, leveraging its unprecedented UV sensitivity and wide-field survey capabilities. We employed synthetic spectral libraries of M0-M8 type giants for primary stars and the ATLAS 9 atmospheric model grid for companion stars spanning a temperature range of 6500 K to 12000 K. By convolving these model spectra with the CSST multi-band filter system, we computed color-color diagrams (g-y versus NUV-u) to construct a diagnostic grid. This grid incorporates interstellar extinction corrections and establishes a framework for identifying AGB binary candidates through direct comparison between observed photometry and theoretical predictions. Furthermore, we discuss the physical origins of UV excess in AGB stars. This study pioneers a diagnostic framework leveraging CSST's unique multi-band UV-visible synergy to construct color-color grids for binary candidate identification, overcoming limitations of non-simultaneous multi-instrument observations.
We present DELTA (Data-Empiric Learned Tidal Alignments), a deep learning model that isolates galaxy intrinsic alignments (IAs) from weak lensing distortions using only observational data. The model uses an Equivariant Graph Neural Network backbone suitable for capturing information from the local galaxy environment, in conjunction with a probabilistic orientation output. Unlike parametric models, DELTA flexibly learns the relationship between galaxy shapes and their local environments, without assuming an explicit IA form or relying on simulations. When applied to mock catalogs with realistic noisy IAs injected, it accurately reconstructs the noise-free, pure IA signal. Mapping these alignments provides a direct visualization of IA patterns in the mock catalogs. Combining DELTA with deep learning interpretation techniques provides further insights into the physics driving tidal relationships between galaxies. This new approach to understanding and controlling IAs is suitable for application to joint photometric and spectroscopic surveys such as the combination of upcoming Euclid, Rubin, and DESI datasets.
We discuss the requirements, concepts, simulations, implementation, and calibration of two dual Fabry-Perot based imaging spectropolarimeters, CRISP and CHROMIS, at the Swedish 1-meter Solar Telescope, and CRISP2 that is under construction. These instruments use a combination of a high-resolution and a low-resolution etalon together with an order-sorting prefilter to define the bandpass. The overall design is made robust and stable by tailoring the low-resolution etalon reflectivity to accommodate expected cavity errors from both etalons, and by using a compact optical design that eliminates the need for folding mirrors. By using a telecentric design based on lenses rather than mirrors, image degradation by the FPI system is negligible, as shown in a previous publication, and the throughput of the system is maximised. Initial alignment, and maintaining that alignment over time, is greatly simplified. The telecentric design allows full calibration and/or modelling of essential system parameters to be carried out without interfering with the optical setup. We also discuss briefly the polarimeters developed for CRISP and CHROMIS. The high performance of CRISP and CHROMIS has been demonstrated in an earlier publication through measurements of the granulation contrast and comparisons with similar measurements simultaneously made through broadband continuum filters. Here, we focus on the aspects of the design that are central to enabling high performance and robustness, but also discuss the calibration and processing of the data, and use a few examples of processed data to demonstrate the achievable image and data quality. We put forward a proposal for a similar conceptual design for the European Solar Telescope and conclude by discussing potential problems of the proposed approach to designs of this type. Some aspects of these FPI systems may be of interest also outside the solar community.
Context. Local young stellar associations (LYSAs <50 Myr and <150 pc) are important laboratories to test predictions from star-formation theories. Estimating their ages through various dating techniques with minimal biases is thus of paramount importance. Aims. We aim at determining the ages of LYSAs with the expansion rate dating technique. Methods. We estimate the ages of the LYSAs using literature membership lists, publicly available data (astrometry and radial velocities), and a recent open-source Bayesian code that implements the expansion rate method. This code in combination with simple statistical assumptions allow us to decontaminate, identify possible substructures or populations, and estimate expansion ages. Results. We derive the largest and most methodological homogeneous set of ages of LYSAs. We rediscover three and discover four associations hidden within the literature membership lists of the classical ones. Conclusions. The expansion ages we report here are compatible with literature age estimates. Moreover, our analysis shows that previous age tensions can be explained, in most cases, by the presence of unidentified populations or substructures.
Context. Determining the ages of young stellar systems is fundamental to test and validate current star-formation theories. Aims. We aim at developing a Bayesian version of the expansion rate method that incorporates the a priori knowledge on the stellar system's age and solves some of the caveats of the traditional frequentist approach. Methods. We upgrade an existing Bayesian hierarchical model with additional parameter hierarchies that include, amongst others, the system's age. For this later, we propose prior distributions inspired by literature works. Results. We validate our method on a set of extensive simulations mimicking the properties of real stellar systems. In stellar associations between 10 and 40 Myr and up to 150 pc the errors are <10%. In star forming regions up to 400 pc, the error can be as large as 80% at 3 Myr but it rapidly decreases with increasing age. Conclusions. The Bayesian expansion rate methodology that we present here offers several advantages over the traditional frequentist version. In particular, the Bayesian age estimator is more robust and credible than the commonly used the frequentist ones. This new Bayesian expansion rate method is made publicly available as a module of the free and open-source code Kalkayotl.
T~Coronae Borealis is the nearest symbiotic recurrent nova. Twice in the last two centuries, in 1866 and 1946, the accreted material ignited on the surface of the white dwarf via runaway thermonuclear fusion reactions and produced a nova eruption. Both eruptions occurred approximately midway through a transient state of high luminosity. A possible explanation of such a state is a dwarf-nova-like outburst, which may arise from a transient increase in the mass-transfer rate of the donor star. We simulate the response of an accretion disk to an event of enhanced mass-transfer that is ``interrupted'' by a pre-eruption dip associated to the convective phase leading to the thermonuclear runaway, and model the resulting optical light curve using the parameters of the T~CrB binary. Our model represents the first attempt to reproduce the transient high-accretion state. The observed brightening can be satisfactorily reproduced by models of an accretion disk with a viscosity parameter \alpha = 3, an event of enhanced mass-transfer with a duration of \Delta t = 15\,yr, and quiescent and high-state mass-transfer rates of 2.0 \times 10^{-9} \, M_\odot\,yr^{-1} and 1.9 \times 10^{-7} \, M_\odot\,yr^{-1}, respectively, while the pre-eruption dip can be reproduced by the small, accelerated expansion of the inner disk radius, at an average velocity of 0.02\,km\,s^{-1}. Our model is also capable of reproducing the observed changes in color of T~CrB throughout the transient event.
A major weakness in one-dimensional (1D) stellar structure and evolution modeling is the simplified treatment of convection, which leads to erroneous near-surface stratification and considerable uncertainties in predicted effective temperatures and luminosities of low-mass stars. In a series of preceding works, a novel method for coupling 1D stellar structural models with a grid of 3D surface convection simulations during stellar evolution was developed, at solar metallicity. This 1D-3D coupling method slightly shifts evolutionary tracks relative to standard calculations, meanwhile providing oscillation frequencies that agree more closely with asteroseismic observations. Here we extend this method to model metal-poor and metal-rich FGK-type stars, by implementing interpolations on-the-fly across metallicity (\rm -3 < [Fe/H] < 0.5) for mean 3D models during stellar evolution. We demonstrate quantitatively that the fundamental stellar parameters modeled within our framework are insensitive to the mixing-length parameter. A 20% change in the mixing-length parameter results in evolutionary tracks with a temperature shift of less than 30 K, compared to a difference of over 200 K in standard evolution calculations. Our extension is validated against eclipsing binary systems with extremely precise observational constraints as well as stars in binaries with asteroseismic data. Using a fixed mixing-length parameter that merely governs convective heat transport in the near-adiabatic layers, the 1D-3D coupling method successfully reproduces most observational constraints for all target stars. Coupling 1D stellar evolution models with 3D simulations greatly reduces uncertainties associated with the choice of atmosphere boundary conditions and mixing-length parameters, hence offering a powerful tool for characterizing stars with seismic measurements and determining ages for globular clusters.
We present a cosmology analysis of weak lensing convergence maps using the Neural Field Scattering Transform (NFST) to constrain cosmological parameters. The NFST extends the Wavelet Scattering Transform (WST) by incorporating trainable neural field filters while preserving rotational and translational symmetries. This setup balances flexibility with robustness, ideal for learning in limited training data regimes. We apply the NFST to 500 simulations from the CosmoGrid suite, each providing a total of 1000 square degrees of noiseless weak lensing convergence maps. We use the resulting learned field compression to model the posterior over \Omega_m, \sigma_8, and w in a wCDM cosmology. The NFST consistently outperforms the WST benchmark, achieving a 16% increase in the average posterior probability density assigned to test data. Further, the NFST improves direct parameter prediction precision on \sigma_8 by 6% and w by 11%. We also introduce a new visualization technique to interpret the learned filters in physical space and show that the NFST adapts its feature extraction to capture task-specific information. These results establish the NFST as a promising tool for extracting maximal cosmological information from the non-Gaussian information in upcoming large-scale structure surveys, without requiring large simulated training datasets.
Recent works have found evidence of significant intrasystem uniformity in planet properties such as radius, mass, and orbital spacing, collectively termed 'peas in a pod' trends. In particular, correlations in radius and mass have been interpreted as implying uniformity in planet bulk density and composition within a system. However, the samples used to assess trends in mass tend to be small and biased. In this paper, we re-evaluate correlations in planet properties in a large sample of systems with at least two planets for which mass and radius have been directly measured, and therefore bulk density can be calculated. Our sample was assembled using the most up-to-date exoplanet catalogue data, and we compute the relevant statistics while using a procedure to 'weight' the data points according to measurement precision. We find a moderate correlation in radius and a weak correlation in the densities of adjacent planets. However, masses of neighbouring planets show no overall correlation in our main sample and a weak correlation among pairs of planets similar in size or pairs restricted to Mp<100 M_Earth, Rp<10 R_Earth. Similarly, we show that the intrasystem dispersion in radius is typically less than that in mass and density. We identify ranges in stellar host properties that correlate with stronger uniformity in pairs of adjacent planets: low Teff for planet masses, and low metallicity and old age for planet densities. Furthermore, we explore whether peas in a pod trends extend into planet compositions or interior structures. For small neighbouring planets with similar radii, we show that their masses and interior structures are often disparate, indicating that even within the same system, similarity in radii is not necessarily a good proxy for similarity in composition or the physical nature of the planets.
Wide Binaries (WBs) are interesting systems to test Newton-Einstein gravity in low potentials. The basic concept is to verify whether the difference in velocity between the WB components is compatible with what is expected from the Newton law. Previous attempts, based solely on Gaia proper motion differences scaled to transverse velocity differences using mean parallax distances, do not provide conclusive results. Here we add to the Gaia transverse velocities precise measurements of the third velocity component, the radial velocity (RV), in order to identify multiple stars, and to improve the reliability of the test by using velocity differences and positions in three dimensions. We use the HARPS spectra to determine accurate RV difference between the WB components, correcting the observed velocities for gravitational redshift and convective shift. We exploit the Gaia distance distributions to determine the projected and intrinsic separations s and r and the 3-dimensional velocity differences of the binaries. Of the 44 pairs observed with HARPS, 27% show sign of multiplicity or are not suitable for the test, and 32 bona-fide WBs survive our selection. Their projected separation s is up to 14 kAU, or 0.06 parsec. We determine distances, eccentricities and position angles to reproduce the velocity differences according to Newton's law, finding reasonable solutions for all WBs but one, and with some systems possibly too near pericenter and/or at too high inclination. Our (limited) number of WBs does not show obvious trends with separation or acceleration and is consistent with Newtonian dynamics. We are collecting a larger sample of this kind to robustly assess these results.
We report seismic analyses of five pulsating subdwarf B (sdBV) stars observed during Kepler's K2 mission, each with a white dwarf companion. We find three of the five to be g-mode-dominated hybrid pulsators. For the other two, we only detect g modes. We determine rotation periods from frequency multiplets for four stars and each rotates subsynchronously to its binary period, including PG 0101+039 and PG 0902+124 both with binary periods near 0.57 days and spin periods near 9 days. We detect frequency multiplets in both p and g modes for PG 0101+039 and LT Cnc and determine that PG 0101+039 rotates like a solid body while LT Cnc rotates differentially radially with the envelope spinning faster than deeper layers. Mostly we find these five stars to be quite similar to one another, spectroscopically and seismically. We find the p modes of the three hybrid pulsators to have gaps between regions of power, which we interpret as overtones and apply a technique to assign modes. We examine their g mode period spacings and deviations thereof and again, find the stars to be similar with period spacings near the average of 250 s and deviations mostly under 25 s. We compare Kepler-observed sdBV stars of different binary types and likely-single pulsators.
Euclid mission is designed to understand the dark sector of the universe. Precise redshift measurements are provided by H2RG detectors. We propose an unbiased method of fitting the flux with Poisson distributed and correlated data, which has an analytic solution and provides a reliable quality factor - fundamental features to ensure the goals of the mission. We compare our method to other techniques of signal estimation and illustrate the anomaly detection on the flight like detectors. Although our discussion is focused on Euclid NISP instrument, much of what is discussed will be of interest to any mission using similar near-infrared sensors
We present the readout noise reduction methods and the 1/f noise response of an 2Kx2K HgCdTe detector similar to the detectors that will be used in the Near Infrared Spectrometer Photometer - one of the instruments of the future ESA mission named Euclid. Various algorithms of common modes subtraction are defined and compared. We show that the readout noise can be lowered by 60% using properly the references provided within the array. A predictive model of the 1/f noise with a given frequency power spectrum is defined and compared to data taken in a wide range of sampling frequencies. In view of this model the definition of ad-hoc readout noises for different sampling can be avoided
A new generation of optical intensity interferometers are emerging in recent years taking advantage of the existing infrastructure of Imaging Atmospheric Cherenkov Telescopes (IACTs). The MAGIC SII (Stellar Intensity Interferometer) in La Palma, Spain, has been operating since its first successful measurements in 2019 and its current design allows it to operate regularly. The current setup is ready to follow up on bright optical transients, as changing from regular gamma-ray observations to SII mode can be done in a matter of minutes. A paper studying the system performance, first measurements and future upgrades has been recently published. MAGIC SII's first scientific results are the measurement of the angular size of 22 stars, 13 of which with no previous measurements in the B band. More recently the Large Sized Telescope prototype from the Cherenkov Telescope Array Observatory (CTAOLST1) has been upgraded to operate together with MAGIC as a SII, leading to its first correlation measurements at the beginning of 2024. MAGIC+CTAO-LST1 SII will be further upgraded by adding the remaining CTAOLSTs at the north site to the system (which are foreseen to be built by the end of 2025). MAGIC+CTAO-LST1 SII shows a feasible technical solution to extend SII to the whole CTAO.
Fast radio bursts (FRBs) are enigmatic extragalactic radio transients with unknown origins. We performed comprehensive Monte Carlo simulations based on the first CHIME/FRB catalog to test whether the FRB population tracks the cosmic star formation history directly or requires a delay. By fully considering CHIME's complex selection effects and beam response, we find that the hypothesis that the FRB population tracks the SFH is not ruled out by the current data, although a small delay is preferred. This is consistent with the scenario in which young magnetars formed through core-collapse supernovae serve as the progenitors of FRBs. However, we estimate the local volumetric rate of FRB sources with energy above 10^{38} erg to be 2.3^{+2.4}_{-1.2} \times 10^5~\rm{Gpc}^{-3}~\rm{yr}^{-1}, which is consistent with previous results. This high volumetric rate means the core-collapse magnetar scenario alone cannot fully account for the observed population. Further theoretical efforts are required to explore alternative or additional progenitor channels for FRBs.
Albeit at a lower frequency than around hotter stars, short-period gas giants around low-mass stars (T_\mathrm{eff} < 4965 K) do exist, despite predictions from planetary population synthesis models that such systems should be exceedingly rare. By combining data from TESS and ground-based follow-up observations, we seek to confirm and characterize giant planets transiting K dwarfs, particularly mid/late K dwarfs. Photometric data were obtained from the TESS mission, supplemented by ground-based imaging- and photometric observations, as well as high-resolution spectroscopic data from the CORALIE spectrograph. Radial velocity (RV) measurements were analyzed to confirm the presence of companions. We report the confirmation and characterization of three giants transiting mid-K dwarfs. Within the TOI-2969 system, a giant planet of 1.16\pm 0.04\,M_\mathrm{Jup} and a radius of 1.10 \pm 0.08\,R_\mathrm{Jup} revolves around its K3V host in 1.82 days. The system of TOI-2989 contains a 3.0 \pm 0.2\,M_\mathrm{Jup} giant with a radius of 1.12 \pm 0.05\,R_\mathrm{Jup}, which orbits its K4V host in 3.12 days. The K4V TOI-5300 hosts a giant of 0.6 \pm 0.1\,M_\mathrm{Jup} with a radius of 0.88 \pm 0.08\,R_\mathrm{Jup} and an orbital period of 2.3 days. The equilibrium temperatures of the companions range from 1001 to 1186 K, classifying them as Hot Jupiters. However, they do not present radius inflation. The estimated heavy element masses in their interior, inferred from the mass, radius, and evolutionary models, are 90 \pm 30\,M_\oplus, 114 \pm 30\,M_\oplus, and 84 \pm 21\,M_\oplus, respectively. The heavy element masses are significantly higher than most reported heavy elements for K-dwarf Hot Jupiters. These mass characterizations contribute to the poorly explored population of massive companions around low-mass stars.
We present photometric time series data spanning 2018-2024 that show the effects of temporal dynamics in the binary system KH 15D, a member of the NGC 2264 star forming region. This source exhibits complex eclipsing behavior due to a precessing circumbinary disk or ring that is slightly inclined relative to the orbital plane of the binary. Using g-band and r-band observations from the Zwicky Transient Facility (ZTF) over seven observing seasons, we follow the evolution of the KH 15D lightcurve as it continues to emerge from its deepest observed photometric minimum about 15 years ago. Our observations are consistent with previous models that propose a precessing, warped circumbinary disk orbiting KH 15D. We verify the gradual precession of the disk by quantifying the times of eclipse ingresses and egresses. We also examine the central re-brightening within the minima of the phased lightcurve. This feature has increased in amplitude over our observing seasons, and we measure its evolution in both amplitude and phase from year to year. Finally, we assess color as a function of phase and brightness. Our findings support the assertion that line-of-sight variations in disk density and structure, possibly due to clumping, coupled with a precessing circumbinary disk are responsible for the central re-brightening event.
We derive the full covariance matrix formulae are derived for proper treatment of correlations in signal fitting procedures, extending the results from previous publications. The straight line fits performed with these matrices demonstrate that a significantly higher signal to noise is obtained when the fluence exceeds 1 e/sec/pix in particular in long (several hundreds of seconds) spectroscopic exposures. The improvement arising from the covariance matrix is particularly strong for the initial intercept of the fit at t=0, a quantity which provides a useful redundancy to cross check the signal quality. We demonstrate that the mode that maximizes the signal to noise ratio in all ranges of fluxes studied in this paper is the one that uses all the frames sampled during the exposure. While at low flux there is no restriction on the organization of frames within groups for fluxes lower than 1 e/sec/pix, for fluxes exceeding this value the coadding of frames shell be avoided.
Ultralight axion-like dark matter (ALDM) is a leading candidate in the dark matter realm, characterized by its prominent wave properties on astronomical scales. Pulsar Timing Arrays (PTAs) and Pulsar Polarization Arrays (PPAs) aim to detect this dark matter through timing and polarization measurements, respectively, of pulsars. The PTA relies on gravitational effects, as the ALDM halo perturbs the spacetime metric within the Milky Way, while the PPA detects non-gravitational effects, namely cosmological birefringence induced by the ALDM Chern-Simons coupling with photons. These two methods complement each other, synergistically enhancing the pulsar array's capability to identify the ALDM signals in the data. In this article, we provide a foundational development of this synergy. We begin by revisiting previously derived two-point correlation functions for both PTA and PPA, and expand our analysis to include correlations between timing and polarization signals. We then construct likelihood functions for PTA and combined PTA-PPA analyses within a Bayesian framework, aimed at detecting the characteristic correlations of ALDM signals. We emphasize the non-Gaussianity of the ALDM timing signals, which arises from their non-linear dependence on the field, in contrast to the Gaussian nature of its polarization signals. To address the complexities introduced, we approach this investigation in two ways: one involves a Gaussian approximation with proper justifications, while the other derives the formalism from the generic Gaussian characteristics of the ALDM field. We anticipate that these efforts will lead to further developments in PTA and PTA-PPA analysis methods.
The Chinese Pulsar Timing Array (CPTA) has collected observations from 57 millisecond pulsars using the Five-hundred-meter Aperture Spherical Radio Telescope (FAST) for close to three years, for the purpose of searching for gravitational waves (GWs). To robustly search for ultra-low-frequency GWs, pulsar timing arrays (PTAs) need to use models to describe the noise from the individual pulsars. We report on the results from the single pulsar noise analysis of the CPTA data release I (DR1). Conventionally, power laws in the frequency domain are used to describe pulsar red noise and dispersion measurement (DM) variations over time. Employing Bayesian methods, we found the choice of number and range of frequency bins with the highest evidence for each pulsar individually. A comparison between a dataset using DM piecewise measured (DMX) values and a power-law Gaussian process to describe the DM variations shows strong Bayesian evidence in favour of the power-law model. Furthermore, we demonstrate that the constraints obtained from four independent software packages are very consistent with each other. The short time span of the CPTA DR1, paired with the large sensitivity of FAST, has proved to be a challenge for the conventional noise model using a power law. This mainly shows in the difficulty to separate different noise terms due to their covariances with each other. Nineteen pulsars are found to display covariances between the short-term white noise and long-term red and DM noise. With future CPTA datasets, we expect that the degeneracy can be broken. Finally, we compared the CPTA DR1 results against the noise properties found by other PTA collaborations. While we can see broad agreement, there is some tension between different PTA datasets for some of the overlapping pulsars. This could be due to the differences in the methods and frequency range compared to the other PTAs.
Context. The high precision of recent asteroseismic observations of red-giant stars has revealed the presence of mixed dipole modes in their oscillation spectra. These modes allow for a look inside the stars. Among the parameters used to characterize mixed modes is the coupling strength q, which is sensitive to the stellar structure in the evanescent zone near the bottom of the convective envelope. Aims. The aim of this work is to probe the validity of the weak and strong coupling approximations, commonly used to calculate q, during stellar evolution along the red-giant branch (RGB). Methods. To test the approximations empirically, we calculate q-values in both, the weak and strong limit for stellar models on the RGB and compare them to the coupling derived from the mixed mode frequency pattern obtained from numerical solutions to the oscillation equations. Results. We find good agreement with the strong coupling approximation on the early RGB, when the evanescent zone lies in the radiative layer right above the hydrogen-burning shell; and with the weak coupling approximation once the evanescent zone is situated in the convective envelope. This is consistent with earlier studies. Additionally, we find that it is viable to use the weak coupling approximation as an estimate for q in the intermediate regime, in the mass range considered in this work (1.00 Msun <= M <= 2.00 Msun). Conclusions. The width of the evanescent zone serves as a good measure for which approximation to use. The serendipitous alignment of the weak coupling approximation with the observable q in the regime where neither approximation is expected to be valid simplifies the asymptotic calculation of mixed mode properties.