2025-05-16 | | Total: 32
Recent X-ray observations have discovered a class of periodic X-ray flares in galactic nuclei known as quasi-periodic eruptions (QPEs). A promising explanation of QPEs is an emission produced when a stellar-mass object orbiting a central supermassive black hole crosses an accretion disk. If the companion is a compact object, such systems would be a prospective multi-messenger target for the space-based observatory LISA and its successors. Here we quantify the prospects for joint X-ray and GW detection of QPEs with orbital frequency in the mHz band using a minimal flare-emission model. Our analysis shows that X-ray observations are most effective at orbital frequencies up to roughly 1 mHz, whereas LISA is sensitive chiefly above about 1 mHz. Because the optimal sensitivity windows overlap only marginally, we predict at most one joint detection during LISA's nominal mission lifetime. Extending GW sensitivity into the sub-millihertz regime (< 0.1 mHz) would raise the possibility of the joint detection by an order of magnitude, enabling QPEs as an interesting multi-messenger target.
Compact binaries with accreting white dwarfs, neutron stars, or black holes are crucial for understanding accretion physics. In this study, we identify accreting compact binary candidates in the eROSITA DR1 by combining Gaia astrometry and ZTF time-domain photometry. We select accreting compact binary candidates whose X-ray emission is likely accretion-driven, based on their X-ray-to-optical flux ratios and optical colors, and further identify systems with short-period ZTF light curves, thereby ruling out AGNs. The resulting 22 candidates, including two known compact binaries, represent promising candidates for spectroscopic follow-up to confirm their accreting nature and explore their physical properties. Our results demonstrate the effectiveness of combining X-ray-to-optical flux ratio and optical color, and time-domain photometry in selecting accreting compact binaries. The method is scalable and adaptable to future multi-wavelength sky surveys, offering a promising path toward a more complete census of compact binaries in the Galaxy.
Fast radio bursts (FRBs) are mysterious millisecond-duration radio transients from the distant universe. Some of them repeat, while others do not. In order to explore their origin, periodic examinations have been conducted on repeating FRBs. Most of them show irregular properties, including burst rate, dispersion measure (DM), and rotation measure (RM). A notable exception is FRB~20180916B, which shows a significant 16-day periodic burst rate. Possible periodic activities have also been reported in FRB~20121102A and FRB~20240209A. However, periodic studies of other properties are sparse. FRB~20220529 was monitored by the Five-hundred-meter Aperture Spherical radio Telescope (FAST) for nearly three years, enabling periodic examinations of its properties. Here we report a possible period of ∼200 days in the RM evolution, with a significance of 4.2 \sigma estimated via the Lomb-Scargle algorithm and 3.5 \sigma with a phase-folding method. The burst rate was also examined for periodicity. This is consistent with the binary origin indicated by the significant RM increase and its prompt recovery.
We investigate the generation of magnetic fields above black hole accretion discs due to the non-zero curl of the disc radiation field. By self consistently computing the components of the radiation flux and their curl, we show that the rotational nature of the radiation field induces charge separation, leading to magnetic field generation in the plasma above the disc. Solving the magnetohydrodynamic equations, we derive the time evolution of these fields and demonstrate that they grow over astrophysically relevant timescales. For a standard Keplerian accretion disc, the produced magnetic fields remain weak, on the order of a few Gauss, consistent with previous predictions. However, when a luminous corona is present in the inner disc region (r_d < 3-10 r_g), the generated fields reach dynamically significant strengths of up to 10^5 Gauss where the magnetic energy density approaching few percentage of equipartition with the gas pressure. These fields develop within realistic growth timescales (such as viscous timescale) and can be dynamically significant in governing disc and jet evolution. Our findings suggest that radiation-driven magnetic fields play a crucial role in accretion flow magnetization, influencing both disc dynamics and observational signatures. The predicted field strengths could affect the thermal emission, synchrotron radiation, and polarization properties of black hole accretion systems, with implications for X-ray binaries, AGN, and jet formation. Future numerical simulations and high-resolution polarimetric observations, such as those from IXPE, eXTP, and EHT, may provide observational confirmation of our findings.
Interpreting the data from radio detectors for extensive air showers typically relies on Monte-Carlo based simulation codes, which, despite their accuracy are computationally expensive and present bottlenecks for analyses. To address this issue we developed a novel method called template synthesis, which synthesises the radio emission from cosmic ray air showers in seconds. This hybrid approach uses a microscopically simulated, sliced shower (the origin) as an input. It rescales the emission from each slice individually to synthesise the emission from a shower with different properties (the target). In order to be able to change the arrival direction during synthesis, we adjust the phases based on the expected geometrical delays. We benchmark the method by comparing synthesised traces to CoREAS simulations over a wide frequency range of [30, 500] MHz . The synthesis quality is primarily influenced by the difference in X_{max} between the origin and target shower. When \Delta X_{max} \leq 100 g/cm^2 , the scatter on the maximum amplitudes of the geomagnetic traces is at most 4%. For the traces from the charge-excess component this scatter is smaller than 6%. We also observe a bias with \Delta X_{max} up to 5% for both components, which appears to depend on the antenna position. Since the bias is symmetrical around \Delta X_{max} = 0 g/cm^2, we can use an interpolation approach to correct for it. We have implemented the template synthesis algorithm in a Python package called \texttt{SMIET}, which includes all the necessary parameters. This package has been successfully tested with air showers with zenith angles up to 50^{\circ} and can be used with any atmosphere, observation level and magnetic field. We demonstrate that the synthesis quality remains comparable to our main benchmarks across various scenarios and discuss use cases, including machine-learning-based analyses.
Recently, the gravitational wave (GW) event for black hole (BH)-BH mergers, S241125n, has been reported to be associated with a short gamma-ray burst (GRB) and X-ray afterglow emission. Such an association could potentially unveil the environments of mergers and provide attractive targets for multi-messenger observations. We summarize these observations and model it as a binary BH merger occurring within an active galactic nucleus (AGN), where the merger remnant accretes disk material at hyper-Eddington rates. The resulting jet could lead to the GRB associated with the GW event. The GRB, detected by Swift-BAT, exhibits a Comptonized spectrum with an unusually soft photon index, which is consistent with emission from the jet interacting the emission associated with the dense AGN disk environment. Its X-ray afterglow shows an unusually hard spectrum, indicating possible strong absorption by a high column density typical of AGN disks. We highlight the importance of identifying the host galaxy and conducting more sensitive infrared observations to test the model.
We investigate the origin of the elliptical ring structure observed in the images of the supermassive black hole M87*, aiming to disentangle contributions from gravitational, astrophysical, and imaging effects. Leveraging the enhanced capabilities of the Event Horizon Telescope (EHT) 2018 array, including improved (u,v)-coverage from the Greenland Telescope, we measure the ring's ellipticity using five independent imaging methods, obtaining a consistent average value of \tau = 0.08_{-0.02}^{+0.03} with a position angle \xi = 50.1_{-7.6}^{+6.2} degrees. To interpret this measurement, we compare against General Relativistic Magnetohydrodynamic (GRMHD) simulations spanning a wide range of physical parameters including thermal or non-thermal electron distribution function, spins, and ion-to-electron temperature ratios in both low and high-density regions. We find no statistically significant correlation between spin and ellipticity in GRMHD images. Instead, we identify a correlation between ellipticity and the fraction of non-ring emission, particularly in non-thermal models and models with higher jet emission. These results indicate that the ellipticity measured from the \m87 emission structure is consistent with that expected from simulations of turbulent accretion flows around black holes, where it is dominated by astrophysical effects rather than gravitational ones. Future high-resolution imaging, including space very long baseline interferometry and long-term monitoring, will be essential to isolate gravitational signatures from astrophysical effects.
Recent analyses on the properties of the central compact object in the HESS J1731-347 remnant and the PSR J1231-1411 pulsar indicated that these two compact objects are characterized by similar (low) masses and possibly different radii. This paper aims at reconciling the aforementioned measurements by utilizing the widely employed color-flavor locked (CFL) MIT bag model. The main objective is related to the examination of the acceptable values for the color superconducting gap \Delta and the bag parameter B. Furthermore, our analysis involves two distinct hypotheses for the nature of compact stars. Firstly, we considered the case of absolute stability for strange quark matter and we found that it is possible to explain both measurements, while also respecting the latest astronomical constraints on the masses and radii of compact stars. Secondly, we studied the case of hybrid stellar matter (transition from hadrons to quarks), and concluded that, when early phase transitions are considered, the simultaneous reconciliation of both measurements leads to results that are inconsistent to the existence of massive compact stars. However, we showed that all current constraints may be satisfied under the consideration that the HESS J1731-347 remnant contains a slow stable hybrid star.
The fundamental plane (FP) relation connects gamma-ray luminosity to intrinsic pulsar properties, offering the potential to estimate distances for radio-quiet (RQ) gamma-ray pulsars, where direct measurements are often unavailable. The Fermi Third Pulsar Catalog presents spectral data for 294 gamma-ray pulsars, including 72 RQ pulsars, of which only 9 have known distances. This study investigates the FP relation's potential to predict distances for RQ gamma-ray pulsars using machine learning (ML) approaches. Ordinary least-squares regression was employed alongside ML methods, including random forest and support vector regression, to predict RQ gamma-ray pulsar distances. To ensure robustness, the analysis considered spectral cutoff significance and outlier data. Results confirm that the FP relation is valid only for pulsars exhibiting significant gamma-ray cutoffs, with FP exponents closely matching theoretical expectations. It was shown that gamma-ray emission origins significantly favor the curvature radiation regime over the synchrotron regime at the 4.2 sigma level. The distances for 62 RQ pulsars, for which no traditional measurements exist, are predicted for the first time. The luminosity and spatial distributions of the considered RQ pulsar population match known pulsar distributions. Predictions for RQ pulsars with known distances are consistent with measured data. This study provides a proof of principle that the FP relation, combined with ML methods, is a promising tool for distance estimation in RQ gamma-ray pulsars, using gamma-ray emission directly for the first time in pulsar distance estimates. Future improvements can be achieved with larger data sets, precise spectral cutoff determinations, and corrections for key factors like the beaming factor.
The GeV-TeV {\gamma}-ray emission is observed from the direction of the source Sgr A^*, which is identified with the SMBH in the centre of our Galaxy. According to some models this {\gamma}-ray emission might originate in the very compact, central region identified with the direct surrounding of the SMBH. Sgr A^* is surrounded by a massive nuclear star cluster. Occasionally these stars might pass close to the line of sight of the observer, resulting in partial absorption of the sub-TeV {\gamma}-ray emission. We investigate the conditions at which an absorption feature appears in the {\gamma}-ray light curves from the Galactic Centre or nuclei of other galaxies containing SMBHs. The detection of such features would allow to obtain constraints on the emission site of {\gamma} rays in active galaxies. We calculate the optical depths for {\gamma} rays in the radiation of individual massive stars, or from the whole population of stars for different parameters of the star cluster. We show that the observer with a line of sight close to the orbital plane of the star can register a {\gamma}-ray absorption dip lasting from a fraction of a day up to a few tens of days. The combined effect of the bulk absorption on the whole population of stars instead can produce a flickering of the observed emission of a red-noise type in the power spectrum of the emission. Predicted absorption features in the sub-TeV {\gamma}-ray light curves from galaxies with SMBHs should be easily detectable by the Large-Sized Telescopes of the Cherenkov Telescope Array Observatory. The discovery of such absorption features provides a unique indication that the {\gamma}-ray production is occurring in a compact region, close to the horizon of the SMBH.
The recent detections of bright optical/infrared kilonova signals following two long-duration gamma-ray bursts (LGRBs), GRB 211211A and GRB 230307A, have significantly challenged the traditional classification of GRBs. These merger-driven LGRBs may represent a distinct GRB population. Since traditional GRB classification methods often struggle to distinguish merger-driven LGRBs from traditional merger-driven short-duration GRBs resulting from compact object mergers and collapse-driven LGRBs produced by massive stars, this work aims to explore the shared properties in terms of hardness, energy, and duration among currently observed merger-driven LGRB events, thereby identifying their observed differences from the traditional GRB population. We collect a sample of merger-driven LGRBs with known redshifts, including observed information on their main emission (ME) and whole emission (WE) phases. Treating ME and WE properties as two independent sets of information, we apply several GRB classification methodologies to explore their potential shared properties. Using the phenomenologically defined Energy-Hardness (EH) parameter, we identify a probable universal linear correlation across merger-driven LGRBs, regardless of whether their ME or WE phases are considered. We propose that such shared properties of merger-driven LGRBs are unlikely to arise from the low-redshift selection effect and they become particularly intriguing when compared with the relatively weak correlations or lack of correlation observed in traditional merger-driven short-duration GRBs (with or without extended emissions) and collapse-driven LGRBs. Our newly proposed correlation highlights the necessity for further investigation into the observations of merger-driven LGRBs and the physical mechanisms underlying the empirical correlation.
The wind-fed accretion of gas with an intrinsic magnetic field onto neutron stars in high-mass X-ray binaries (HMXBs) is discussed. The criteria for the implementation of quasi-spherical accretion in the system, the formation of a Keplerian accretion disk and a non-Keplerian magnetically arrested disk (MAD-accretion scenario) are formulated. It is shown that the period of the equilibrium rotation of a neutron star in the MAD-accretion scenario, in addition to the parameters of the binary system, is also determined by the velocity, temperature, and degree of magnetization of the stellar wind of its massive companion. Within the hypothesis of the equilibrium rotation, we have estimated the parameters of the stellar wind for the most fully-studied quasi-equilibrium accreting pulsars in the HMXBs corresponding to the criteria for the MAD-accretion scenario implementation. We show that the period of the equilibrium rotation of a neutron star in the system, all other things being equal, increases with the growth of the magnetization degree of the stellar wind component of its massive companion flowing out in the plane of the system's orbit. Therefore, the periods of accreting pulsars in systems with similar parameters can differ from each other by several orders of magnitude and, under favorable conditions, reach values of several thousand and even tens of thousands of seconds.
We investigate nonradial f-mode oscillations of protoneutron stars in full general relativity, employing equations of state described by the Brueckner-Hartree-Fock theory or the relativistic mean field model, while assuming isentropy and fixed lepton fractions for the internal structure. The validity of various universal relations for cold neutron stars involving f-mode characteristics and macroscopic properties of the star is confirmed for those isentropic protoneutron stars. Prospects of observations are also discussed. According to simulation results, we then model details of the thermal and trapping profiles in a PNS with the canonical mass. The corresponding f-mode frequencies and gravitational-wave strain amplitudes are presented. The validity of the universal relations during the evolution to the formation of a cold neutron star is confirmed.
Dark matter (DM) remains one of the most significant open questions in modern physics, with its nature and interactions largely unexplored. In this study, we investigate the behavior of massive fermionic DM particles in the context of neutron stars (NSs), extending prior studies which focused on the bosonic DM. By incorporating the motion of NSs in the Galaxy and considering scenarios with and without DM self-annihilation, we demonstrate their impact on the DM capture rate and the accumulation process inside NSs. Observational data from pulsars in the Milky Way are used to place constraints on DM properties, including the mass and the DM-nucleon scattering cross-section, offering a more comprehensive picture in probing DM interactions in astrophysical environments.
To understand physical processes such as mass transfer and binary evolution in X-ray binaries, the orbital parameters of the system are fundamental and crucial information. Cygnus X-3 is a high-mass X-ray binary composed of a compact object of unknown nature and a Wolf-Rayet star, which is of great interest in the context of wind-fed mass accretion and binary evolution. Here we present XRISM/Resolve high-resolution spectroscopy focusing on the Fe Ly\alpha lines in its hypersoft state. We perform an orbital phase-resolved spectral analysis of the lines to study the orbital modulation of the emission and absorption lines. It is found that the emission lines reflect the orbital motion of the compact object whose estimated velocity amplitude is 430^{~~+150}_{~~-140}~~\mathrm{km\,s^{~-1}}, while the absorption lines show a variation that can be interpreted as originating from the stellar wind. We discuss possible mass ranges for the binary components using the mass function with the estimated value of the velocity amplitude in this work, combined with the relation between the mass loss rate and the orbital period derivative and the empirical mass and mass loss rate relation for Galactic Wolf-Rayet stars. They are constrained to be (1.3\text{-}5.1)\,M_\odot and (9.3\text{-}12)\,M_\odot for the assumed inclination angle of i = 25 deg, which becomes more relaxed to (1.3\text{-}24)\,M_\odot and (9.3\text{-}16)\,M_\odot for i = 35 deg, respectively. Thus, it remains unclear whether the system harbors a black hole or a neutron star.
We report the detection of a spatial variation of X-ray polarization in the southwestern shell of SN 1006 (SN 1006 SW) using IXPE. The shell has an average X-ray polarization degree (PD) of 21.6\%\pm 4.5\% and polarization angle (PA) of -48^\circ \pm 5^\circ in the 2--4 keV energy band, similar to those in the northeastern shell. The PD varies along SN 1006 SW, with a peak PD= 40\%\pm 8\% in the south and a significantly lower PD \lesssim 27\% (99\% upper limit) in the west where the shell has been proposed to be interacting with an interstellar cloud. The correlation between the PD, which reflects the magnetic orderliness, and the preshock density provides observational evidence that magnetic turbulence and amplification are environment-dependent. The high PD detected in the southern region of the shell constrains the magnetic turbulence scale of \lesssim 0.1~pc. Moreover, by comparing the IXPE X-ray and MeerKAT radio polarization measurements for SN 1006 SW, we found that magnetic fields traced by X-ray polarization are nearly radially distributed, whereas those traced by radio polarization tend to follow a direction parallel to the Galactic plane. This suggests that the X-ray polarization probes freshly amplified magnetic fields from small-scale structures in the immediate postshock region, while the radio traces more extended regions influenced by the pre-existing ambient magnetic fields.
Although the X-ray spectra of Seyfert 1 galaxies exhibit absorption lines of He-like iron and H-like iron at blue-shifted velocities of approximately 500 \ \mathrm{km} \ \mathrm{s}^{-1}, the physical origin of these absorption lines remains uncertain. In this study, we performed X-ray radiative transfer based on the sub-parsec-scale thermally driven outflows. The initial step involved calculating the photoionization equilibrium using the Cloudy code, which is based on three-dimensional radiative hydrodynamic simulations. Subsequently, X-ray radiative transfer was performed using the Monte Carlo simulation for astrophysics and cosmology code. Our findings indicate that when the angle of inclination ranges from 55 \ \mathrm{degrees} to 65 \ \mathrm{degrees}, the transmitted component of the X-ray spectrum displays absorption lines of He-like and H-like iron, exhibiting a blue shift of approximately 500 \ \mathrm{km} \ \mathrm{s}^{-1}. The results suggest that the absorption lines are generated by a photoionized gas within 0.005 \ \mathrm{pc}. Additionally, the results indicate that the scattered component of the X-ray spectrum exhibits emission lines originating from neutral iron fluorescence, He-like iron, and H-like iron. The emission lines are broadened by approximately 7000 \ \mathrm{km} \ \mathrm{s}^{-1} due to the Keplerian rotation. Furthermore, the model reproduced the H-like iron and H-like iron absorption lines in NGC 3783 observed by the Chandra High Energy Transmission Grating.
We investigate the X-ray emission from the Galactic Centre (GC) region, focusing on the 6.4 keV fluorescent line of neutral or weakly ionised iron, which is commonly attributed to X-ray reflection from dense molecular clouds. Our goal is to separate the reflection signal from other physical X-ray components. We aim to produce a clean map of the 6.4 keV emission, thus providing a better understanding of the X-ray reflection processes in the GC. We utilised a deep mosaic of all available XMM-Newton observations, encompassing the central 40 square degrees of the Galaxy. The mosaics of two narrow bands centred at 6.7 keV and 6.4 keV, and a broader continuum band at lower energies (5-6.1 keV), provided valuable spatial and spectral information on the X-ray emission. These combined with the stellar mass distribution of our Galaxy enabled us to decompose the observed signal into physically meaningful components. Our analysis shows that the cleaned 6.4 keV band map, free from the contribution of bright and unresolved point sources, is predominantly shaped by X-ray reflection from dense molecular clouds. The spatial distribution of this emission, which strongly correlates with the molecular gas distribution in the Central Molecular Zone (CMZ), supports the interpretation that this map provides the best estimate of the X-ray reflection signal averaged over the last two decades. The cleaned reflection map produced could serve as a tool for future studies to quantify upper limits on the reflection contribution from low-energy cosmic rays in unilluminated regions. Moreover, we estimate that, on average within the CMZ, approximately 65% of the ridge emission contributes to the observed 6.4 keV emission, a factor that should be incorporated into upcoming investigations of the GC, such as polarisation studies of the reflected X-ray continuum from molecular clouds.
We present the results of four three-dimensional radiation magnetohydrodynamic simulations of accretion disks around a 10^8 solar mass black hole, which produce the far ultraviolet spectrum peak and demonstrate a robust physical mechanism to produce the extreme ultraviolet to soft X-ray power-law continuum component. The disks are fed from rotating tori and reach accretion rates ranging from 0.03 to 4 times the Eddington value. The disks become radiation pressure or magnetic pressure dominated depending on the relative timescales of radiative cooling and gas inflow. Magnetic pressure supported disks can form with or without net poloidal magnetic fields as long as the inflowing gas can cool quickly enough, which can typically happen when the accretion rate is low. We calculate the emerging spectra from these disks using multi-group radiation transport with realistic opacities and find that they typically peak around 10 eV. At accretion rates close to or above the Eddington limit, a power-law component can appear for photon energies between 10 eV and 1 keV with a spectral slope varying between L_\nu\propto\nu^{-1} and \nu^{-2}, comparable to what is observed in radio quiet quasars. The disk with 3\% Eddington accretion rate does not exhibit this component. These high energy photons are produced in an optically thick region \approx 30^{\circ}-45^{\circ} from the disk midplane by compressible bulk Comptonization within the converging accretion flow. Strongly magnetized disks that have a very small surface density will produce a spectrum that is very different from what is observed.
Correlated time-of-arrival measurements by pulsar timing arrays (PTAs) have provided a new means of constraining astrophysical or cosmological models that produce a gravitational wave (GW) background. For this work, we discuss the implications of PTA observations for Gauss-Bonnet (GB) inflationary models through the production and propagation of inflationary GWs. We show that our GB inflationary scenario is consistent with present PTA and cosmological data. A blue-tilted tensor power spectrum supported by PTAs can be naturally accommodated in GB inflation. Using observational constraints, we derive general conditions for the inflaton potential and the GB coupling function, suggesting that in GB inflation, the inflaton must climb up the potential before rolling downhill and reaching the end of inflation. We provide two concrete GB inflationary models to demonstrate the viability of this mechanism.
Very massive stars (VMS) play a fundamental role in astrophysics due to their winds and supernovae (SN), and their role as massive black hole (BH) progenitors. However, their origin and evolution remain a significant challenge. Recent theoretical work and observations suggest that VMS approaching the Eddington limit may experience mass loss above the standard wind predictions. This study investigates how enhanced winds influence single and binary VMS evolution, observable properties, and resulting BH populations. New stellar wind prescriptions, sensitive to the Eddington parameter (\Gamma_e) and the luminosity-to-mass ratio, were implemented into the stellar evolution code PARSEC v2.0. These updated single-star tracks (100 - 600 M_{\odot} at Z=0.006) were used to model the VMS population in the Tarantula Nebula and integrated into the SEVN binary evolution code. The \Gamma_e-enhanced single-star tracks match observed VMS properties better than standard models. Explaining the most massive star, R136a1, through a single-star origin suggests a zero-age main sequence (ZAMS) mass limit of < 400 M_{\odot} regardless of the wind recipe used. However, binary stellar mergers also offer a suitable origin for R136a1 and other observed VMS, potentially lowering the upper ZAMS mass limit by ~100 M_{\odot}. In binaries, enhanced winds inhibit main-sequence stellar mergers and limit BH production above the pair-instability mass gap's lower edge (~50 M_{\odot}). Binary BHs merging in a Hubble time with enhanced winds yield more primary BHs above 30 M_{\odot} and enable secondary BHs between 30-40 solar masses, a range not found with standard stellar winds at LMC metallicity. This study highlights the crucial role that stellar winds and binary interactions play in VMS evolution and offers predictions relevant for interpreting VMS observations and gravitational wave source origins.
We present full-polarization observations at \lambda = 0.87 mm (345 GHz) conducted with the Atacama Large Millimeter/submillimeter Array (ALMA) toward Messier 87 (M87) and seven other radio-loud active galactic nuclei (AGN). We determine polarization and Faraday properties of the targets, with linear polarization (LP) fractions spanning 1\% to 17\%, and Faraday rotation measures (RMs) exceeding 10^{5} rad m^{-2}. These RM values, 1-2 orders of magnitude higher than typically observed at \lambda > 3 mm, suggest denser Faraday screens or stronger magnetic fields in the (sub)millimeter emission regions. Additionally, we present the first submillimeter polarized images of the M87 jet and the observed AGN. In the M87 jet, we identify RM gradients and sign reversals, suggestive of a helical magnetic field structure extending over kiloparsec scales. These observations were part of the ALMA Band~7 Very Long Baseline Interferometry (VLBI) commissioning during the April 2021 Event Horizon Telescope (EHT) campaign. The results provide essential constraints for calibrating, analyzing, and interpreting VLBI data from the EHT at 345 GHz, representing a critical step toward submillimeter VLBI imaging.
Vela.jl is a package for performing Bayesian pulsar timing & noise analysis written in Julia and Python. In the wideband paradigm of pulsar timing, simultaneous time of arrival and dispersion measure measurements are derived from a radio observation using frequency-resolved integrated pulse profiles and templates without splitting the observation into multiple frequency sub-bands. We describe the implementation of the wideband timing paradigm in Vela.jl, and demonstrate its usage using the NANOGrav 12.5-year wideband data of PSR J1923+2515. Vela.jl is the first software package to provide this functionality.
Strontium (Sr) emissions have been observed across a wide range of astrophysical phenomena, from kilonovae (KNe) events to white dwarf (WD) stars. Precise and extensive atomic data for low ionisation stages of Sr is required for accurate theoretical modelling and to improve our understanding of evolutionary pathways. We calculated energy levels, Einstein A coefficients and electron-impact excitation collision strengths for Sr I, Sr III, Sr IV and Sr V at the temperature and density ranges of interest in KNe and WD research. We developed new target structures using the GRASP0 and AUTOSTRUCTURE packages. The energies and A-values arising from the new structures were found to be in good agreement with experimental and theoretical equivalents reported in the literature. Maxwellian averaged electron impact collision strengths were calculated using the R-matrix approach, as applied through the DARC and RMBP coding packages. These are presented in adf04 file format. The new data sets allowed us to construct synthetic spectra for the first five ionisation stages of Sr and probe possible density and temperature diagnostic lines. The synthetic spectra within the KNe regime revealed possible Sr IV and Sr V candidate lines at 1027.69nm and 1203.35nm respectively. These may provide useful benchmarks for determining the extent of Sr ionisation that can be reached in an evolving KNe event. Additional diagnostic lines were found to be poor across the Sr ion stages for both KNe and WD regimes due to most levels being in either coronal or Local Thermodynamic Equilibrium (LTE) conditions.
Quasars, powered by accretion onto supermassive black holes (SMBHs), exhibit significant variability, offering insights into the physics of accretion and the properties of the central engines. In this study, we analyze photometric variability and its correlation with key quasar properties, including black hole mass (M_{\mathrm{BH}}) and nuclear luminosities, using 915 quasars with 0\leq z<3.0 from the AQMES sample monitored within SDSS-V. Variability metrics were derived from approximately 6-year light curves provided by the Zwicky Transient Facility -- ZTF, while SMBH masses and luminosities were obtained from the SDSS DR16 quasar catalog of \citet{wu2022catalog}. We identify a strong anti-correlation between variability amplitude and luminosity, which strengthens with redshift, and a redshift-dependent trend for M_{\mathrm{BH}}: a positive correlation at low redshifts, no significant correlation at intermediate redshifts, and an anti-correlation the highest redshifts. Our main finding is a robust anti-correlation between photometric variability amplitude and Eddington ratio, consistent across redshift bins. We present a general equation encapsulating this relationship, that appears to be almost free of redshift dependence, enabling predictions of quasar variability based on accretion parameters {\bf or vice-versa}. The derived relation with the Eddington ratio provides a unified framework for interpreting variability in active galactic nuclei (AGN) and facilitates future studies of quasar variability using high-cadence surveys, such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time -- LSST.