Source Luminosity Research Articles

The potential role of binary neutron star merger afterglows in multimessenger cosmology

Binary neutron star mergers offer a new and independent means of measuring the Hubble constant H0 by combining the gravitational-wave inferred source luminosity distance with its redshift obtained from electromagnetic follow-up. This method is limited by the intrinsic degeneracy between the system distance and orbital inclination in the gravitational-wave signal. Observing the afterglow counterpart to a merger can further constrain the inclination angle, allowing this degeneracy to be partially lifted and improving the measurement of H0. In the case of the binary neutron star merger GW170817, afterglow light-curve and imaging modeling thus allowed the H0 measurement to be improved by a factor of three. However, systematic access to afterglow data is far from guaranteed. In fact, though each one allows a leap in H0 precision, these afterglow counterparts should prove rare in forthcoming multimessenger campaigns. We combine models for emission and detection of gravitational-wave and electromagnetic radiation from binary neutron star mergers with realistic population models and estimates for afterglow inclination angle constraints. Using these models, we quantify how fast H0 will be narrowed down by successive multimessenger events with and without the afterglow. We find that because of its rareness and though it greatly refines angle estimates, the afterglow counterpart should not significantly contribute to the measurement of H0 in the long run.

Discovery of 22 GHz Water Masers in the Serpens South Region

Abstract Using the Karl G. Jansky Very Large Array (VLA), we have conducted a survey for 22 GHz, 61,6–52,3 H2O masers toward the Serpens South region. The masers were also observed with the Very Long Baseline Array following the VLA detections. We detect for the first time H2O masers in the Serpens South region that are found to be associated to three Class 0–Class I objects, including the two brightest protostars in the Serpens South cluster, known as CARMA-6 and CARMA-7. We also detect H2O masers associated to a source with no outflow or jet features. We suggest that this source is most probably a background asymptotic giant branch star projected in the direction of Serpens South. The spatial distribution of the emission spots suggest that the masers in the three Class 0–Class I objects emerge very close to the protostars and are likely excited in shocks driven by the interaction between a protostellar jet and the circumstellar material. Based on the comparison of the distributions of bolometric luminosity of sources hosting 22 GHz H2O masers and 162 young stellar objects covered by our observations, we identify a limit of L Bol ≈ 10L ⊙ for a source to host water masers. However, the maser emission shows strong variability in both intensity and velocity spread, and therefore masers associated to lower-luminosity sources may have been missed by our observations. We also report 11 new sources with radio continuum emission at 22 GHz.

Observation of Gravitational Waves from Two Neutron Star–Black Hole Coalescences

Abstract We report the observation of gravitational waves from two compact binary coalescences in LIGO’s and Virgo’s third observing run with properties consistent with neutron star–black hole (NSBH) binaries. The two events are named GW200105_162426 and GW200115_042309, abbreviated as GW200105 and GW200115; the first was observed by LIGO Livingston and Virgo and the second by all three LIGO–Virgo detectors. The source of GW200105 has component masses 8 . 9 − 1 . 5 + 1. 2 and 1. 9 − 0. 2 + 0. 3 M ⊙ , whereas the source of GW200115 has component masses 5. 7 − 2 . 1 + 1. 8 and 1. 5 − 0. 3 + 0. 7 M ⊙ (all measurements quoted at the 90% credible level). The probability that the secondary’s mass is below the maximal mass of a neutron star is 89%–96% and 87%–98%, respectively, for GW200105 and GW200115, with the ranges arising from different astrophysical assumptions. The source luminosity distances are 280 − 110 + 110 and 300 − 100 + 150 Mpc , respectively. The magnitude of the primary spin of GW200105 is less than 0.23 at the 90% credible level, and its orientation is unconstrained. For GW200115, the primary spin has a negative spin projection onto the orbital angular momentum at 88% probability. We are unable to constrain the spin or tidal deformation of the secondary component for either event. We infer an NSBH merger rate density of 45 − 33 + 75 Gpc − 3 yr − 1 when assuming that GW200105 and GW200115 are representative of the NSBH population or 130 − 69 + 112 Gpc − 3 yr − 1 under the assumption of a broader distribution of component masses.

Detection of X-ray pulsations at the lowest observed luminosity of Be/X-ray binary pulsar EXO 2030+375 with AstroSat

We present the results obtained from timing and spectral studies of Be/X-ray binary pulsar EXO 2030+375 using observations with the Large Area Xenon Proportional Counters and Soft X-ray Telescope of AstroSat, at various phases of its Type-I outbursts in 2016, 2018, and 2020. The pulsar was faint during these observations as compared to earlier observations with other observatories. At the lowest luminosity of 2.5\(\times \)10\(^{35}\) erg s\(^{-1}\) in 0.5–30 keV energy range, \(\approx \)41.3 s pulsations were clearly detected in the X-ray light curves. This finding establishes the first firm detection of pulsations in EXO 2030+375 at an extremely low mass accretion rate to date. The shape of the pulse profiles is complex due to the presence of several narrow dips. Though pulsations were detected up to \(\sim \)80 keV when the source was brighter, pulsations were limited up to \(\sim \)25 keV during the third AstroSat observation at lowest source luminosity. A search for quasi-periodic oscillations in \(2\times 10^{-4}\) Hz to 10 Hz yielded a negative result. Spectral analysis of the AstroSat data showed that the spectrum of the pulsar was steep with a power-law index of \(\sim \)2. The values of photon-indices at observed low luminosities follow the known pattern in sub-critical regime of the pulsar.

Study of the X-ray Pulsar XTE J1946+274 with NuSTAR

We present the results of our spectral and timing analysis of the emission from the transient X-ray pulsar XTE J1946+274 based on the simultaneous NuSTAR and Swift/XRT observations in the broad energy range 0.3-79 keV carried out in June 2018 during a bright outburst. Our spectral analysis has confirmed the presence of a cyclotron absorption line at an energy $\sim38$ keV in both averaged and phase-resolved spectra of the source. Phase-resolved spectroscopy has also allowed the variation in spectral parameters with neutron star rotation phase, whose period is $\simeq15.755$ s, to be studied. The energy of the cyclotron line is shown to change significantly (from $\simeq34$ to $\simeq39$ keV) on the scale of a pulse, with the line width and optical depth also exhibiting variability. The observed behavior of the cyclotron line parameters can be interpreted in terms of the model of the reflection of emission from a small accretion column (the source's luminosity at the time of its observations was $\sim 3 \times 10^{37}$ erg s$^{-1}$) off the neutron star surface. The equivalent width of the iron line has been found to also change significantly with pulse phase. The time delay between the pulse and equivalent width profiles can be explained by the reflection of neutron star emission from the outer accretion disk regions.

Comprehensive Gas Characterization of a z = 2.5 Protocluster: A Cluster Core Caught in the Beginning of Virialization?

Abstract In order to connect galaxy clusters to their progenitor protoclusters, we must constrain the star formation histories within their member galaxies and the timescale of virial collapse. In this paper we characterize the complex star-forming properties of a z = 2.5 protocluster in the COSMOS field using ALMA dust continuum and new Very Large Array CO (1–0) observations of two filaments associated with the structure, sometimes referred to as the “Hyperion” protocluster. We focus in particular on the protocluster “core,” which has previously been suggested as the highest-redshift bona fide galaxy cluster traced by extended X-ray emission in a stacked Chandra/XMM image. We reanalyze these data and refute these claims, finding that at least 40% ± 17% of extended X-ray sources of similar luminosity and size at this redshift arise instead from inverse Compton scattering off recently extinguished radio galaxies rather than intracluster medium. Using ancillary COSMOS data, we also constrain the spectral energy distributions of the two filaments’ eight constituent galaxies from the rest-frame UV to radio. We do not find evidence for enhanced star formation efficiency in the core and conclude that the constituent galaxies are already massive (M ⋆ ≈ 1011 M ⊙), with molecular gas reservoirs >1010 M ⊙ that will be depleted within 200–400 Myr. Finally, we calculate the halo mass of the nested core at z = 2.5 and conclude that it will collapse into a cluster of (2–9) × 1014 M ⊙, comparable to the size of the Coma Cluster at z = 0 and accounting for at least 50% of the total estimated halo mass of the extended “Hyperion” structure.

Influence of Irradiation-driven Winds on the Evolution of Intermediate-mass Black Hole X-ray Binaries

Abstract In young dense clusters, an intermediate-mass black hole (IMBH) may get a companion star via exchange encounters or tidal capture and then evolve toward the IMBH X-ray binary by the Roche lobe overflow. It is generally thought that IMBH X-ray binaries are potential ultraluminous X-ray sources (ULXs); hence, their evolution is very significant. However, the irradiation-driven winds by the strong X-ray flux from the accretion disks around the IMBHs play an important role in determining the evolution of IMBH X-ray binaries and should be considered in the detailed binary evolution simulation. Employing the models with the MESA code, we focus on the influence of irradiation-driven winds on the evolution of IMBH X-ray binaries. Our simulations indicate that a high wind-driving efficiency (f = 0.01 for Z = 0.02, and f = 0.002 for Z = 0.001) substantially shortens the duration in the ULX stage of IMBH X-ray binaries with an intermediate-mass (5 M ⊙) donor star. However, this effect can be ignored for high-mass (10 M ⊙) donor stars. The irradiation effect (f = 0.01 or 0.002) markedly shrinks the initial parameter space of IMBH binaries evolving toward ULXs with high luminosity (L X > 1040 erg s−1) and hyperluminous X-ray sources in the donor-star mass versus orbital period diagram. Furthermore, the irradiation effect results in an efficient angular momentum loss, yielding to IMBH X-ray binaries with relatively close orbits. In our simulated parameter space, about 1% of IMBH binaries would evolve toward compact X-ray sources owing to short initial orbital periods, some of which might be detected as low-frequency gravitational-wave sources.

Tight multimessenger constraints on the neutron star equation of state from GW170817 and a forward model for kilonova light-curve synthesis

ABSTRACT We present a rapid analytic framework for predicting kilonova light curves following neutron star (NS) mergers, where the main input parameters are binary-based properties measurable by gravitational wave detectors (chirp mass and mass ratio, orbital inclination) and properties dependent on the nuclear equation of state (tidal deformability, maximum NS mass). This enables synthesis of a kilonova sample for any NS source population, or determination of the observing depth needed to detect a live kilonova given gravitational wave source parameters in low latency. We validate this code, implemented in the public mosfit package, by fitting it to GW170817. A Bayes factor analysis overwhelmingly (B > 1010) favours the inclusion of an additional luminosity source in addition to lanthanide-poor dynamical ejecta during the first day. This is well fit by a shock-heated cocoon model, though differences in the ejecta structure, opacity or nuclear heating rate cannot be ruled out as alternatives. The emission thereafter is dominated by a lanthanide-rich viscous wind. We find the mass ratio of the binary is q = 0.92 ± 0.07 (90 per cent credible interval). We place tight constraints on the maximum stable NS mass, MTOV $=2.17^{+0.08}_{-0.11}$ M⊙. For a uniform prior in tidal deformability, the radius of a 1.4-M⊙ NS is R1.4 ∼ 10.7 km. Re-weighting with a prior based on equations of state that support our credible range in MTOV, we derive a final measurement R1.4 $=11.06^{+1.01}_{-0.98}$ km. Applying our code to the second gravitationally detected NS merger, GW190425, we estimate that an associated kilonova would have been fainter (by ∼0.7 mag at 1 d post-merger) and declined faster than GW170817, underlining the importance of tuning follow-up strategies individually for each GW-detected NS merger.

The Hi-GAL compact source catalogue – II. The 360° catalogue of clump physical properties

ABSTRACT We present the 360° catalogue of physical properties of Hi-GAL compact sources, detected between 70 and 500 $\mu$m. This release not only completes the analogous catalogue previously produced by the Hi-GAL collaboration for −71° ≲ ℓ ≲ 67°, but also meaningfully improves it because of a new set of heliocentric distances, 120 808 in total. About a third of the 150 223 entries are located in the newly added portion of the Galactic plane. A first classification based on detection at 70 $\mu$m as a signature of ongoing star-forming activity distinguishes between protostellar sources (23 per cent of the total) and starless sources, with the latter further classified as gravitationally bound (pre-stellar) or unbound. The integral of the spectral energy distribution, including ancillary photometry from λ = 21 to 1100 $\mu$m, gives the source luminosity and other bolometric quantities, while a modified blackbody fitted to data for $\lambda \ge 160~\mu$m yields mass and temperature. All tabulated clump properties are then derived using photometry and heliocentric distance, where possible. Statistics of these quantities are discussed with respect to both source Galactic location and evolutionary stage. No strong differences in the distributions of evolutionary indicators are found between the inner and outer Galaxy. However, masses and densities in the inner Galaxy are on average significantly larger, resulting in a higher number of clumps that are candidates to host massive star formation. Median behaviour of distance-independent parameters tracing source evolutionary status is examined as a function of the Galactocentric radius, showing no clear evidence of correlation with spiral arm positions.

Comptonization as an Origin of the Continuum in Intermediate Polars

Abstract In this paper we test if the ∼0.3–15 keV XMM-Newton EPIC pn spectral continuum of IPs can be described by the thermal Comptonization compTT model. We used publicly observations of 12 IPs (AE Aqr, EX Hya, V1025 Cen, V2731 Oph, RX J2133.7+5107, PQ Gem, NY Lup, V2400 Oph, IGR J00234+6141, IGR J17195-4100, V1223 Sgr, and XY Ari). We find that our modeling is capable of fitting well the average spectral continuum of these sources. In this framework, UV/soft X-ray seed photons (with 〈kT s〉 of 0.096 ± 0.013 keV) coming presumably from the star surface are scattered off by electrons present in an optically thick plasma (with 〈kT e〉 of 3.05 ± 0.16 keV and optical depth 〈τ〉 of 9.5 ± 0.6 for plane geometry) located nearby (on top) to the more central seed photon emission regions. A soft blackbody (bbody) component is observed in 5 out of the 13 observations analyzed, with a mean temperature 〈kT bb 〉 of 0.095 ± 0.004 keV. We observed that the spectra of IPs show in general two photon indices Γ, which are driven by the source luminosity and optical depth. Low luminosity IPs show 〈Γ〉 of 1.83 ± 0.19, whereas high luminosity IPs show lower 〈Γ〉 of 1.34 ± 0.02. Moreover, the good spectral fits of PQ Gem and V2400 Oph indicate that the polar subclass of CVs may be successfully described by the thermal Comptonization as well.

ALMA/NICER observations of GRS 1915+105 indicate a return to a hard state

Context. GRS 1915+105 is a transient black hole X-ray binary consistently emitting 10–100% of the Eddington luminosity in the X-ray band over the last three decades until mid-2018 when the source luminosity suddenly decreased by an order of magnitude. This phase was followed by a change to a state with even lower average X-ray fluxes never seen before during the outburst but presenting renewed flaring activity at different wavelengths, albeit with mean fluxes still in decline. Aims. GRS 1915+105 has the longest orbital period known among low-mass X-ray binaries, the largest accretion disk size, and therefore the largest mass supply for accretion. The high inclination of the disk allows the study of geometrical effects of the accretion flow such as changes in the height-to-radius ratio or the effect of accretion disk winds on the intrinsic emission that is expected during the outburst decay. In addition, the transient jet is expected to change to a compact, self-absorbed, steady jet. Methods. We conducted two full polarization Atacama Large Millimeter Array observations to study the jet properties during the outburst decay by analyzing the spectral, polarization, and intra-epoch variability for both observation epochs. In addition, we analyzed almost daily Neutron Star Interior Composition Explorer pointing observations, modeling X-ray power spectral densities, spectral energy distributions, and light curves with a physically motivated model to follow the changing accretion disk properties throughout the outburst decay and relating them to the jet emission. Results. We show that the X-ray and millimeter (mm) spectral, timing, and polarization properties are consistent with those of a typical decaying X-ray binary outburst and that GRS 1915+105 has descended into the low-luminosity hard X-ray state. The jet emission in the mm is consistent with a compact, steady jet with ∼1% linear polarization, and the magnetic field is likely aligned with the jet position angle. Relating the mm emission to the X-ray emission reveals that the source has changed from a higher radio/X-ray correlation index to a lower one; Lradio ∝ LX0.6.

Radio loudness and classification for radio sources

Abstract The extragalactic radio sources are divided into two subclasses (radio-loud and radio-quiet sources) in the literature using radio loudness (R), which is defined as the ratio of radio emission to optical emission, but the boundary R-value separating the two classes is different in various sources. In this work, a sample of 2419 objects from the 13th catalog of quasars and active nuclei is used to build a boundary for the two subclasses. To do so, we compiled the radio and optical data, calculated their radio and optical indexes, made K-correction, obtained the radio loudness, and adopted a Bayesian analysis method to the logarithm of radio loudness for classification. We also investigated the correlations of radio loudness with radio/optical luminosities. Our main conclusions are summarized as follows: (1) The distribution of the logarithm of radio loudness (log R) is bimodal, the sources with log R < 1.26 are classified as radio-quiet sources, and those with log R > 1.26 are classified as radio-loud ones from the Bayesian analysis method. (2) The average radio-optical effective spectral index of radio-quiet sources is 〈αRO〉 = 0.05, while that of radio-loud sources is 〈αRO〉 = 0.55. (3) There are positive correlations between radio luminosity and radio loudness for both radio-loud sources and radio-quiet sources. (4) A dividing line of separating the distribution of the clusters on the diagram of radio loudness against radio luminosity was obtained statistically to set the boundary between radio-loud sources and radio-quiet sources, with an accuracy of $99.73\%$ based on the classification result from the Bayesian analysis method.

Simulating dust grain-radiation coupling on a moving mesh

ABSTRACT We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code arepo-rt, which solves the moment-based radiative transfer equations on an unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in arepo-rt, which treats UV, optical photons as single scattered and IR photons as multi scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.

Intergalactic Media Densities of Extragalactic Radio Sources and Source Ages

In this paper, we use analytical methods with some plausible assumptions to develop a mathematical model which may be used to obtain estimate of density of the intergalactic medium (IGM). This relation suggests that the plot of the linear sizes against the bolometric luminosities of the extragalactic radio sources in our sample should be able to furnish us with an estimate of the density of IGM in which the radio sources are domiciled. Moreover, for the purpose of obtaining an empirical relation similar to the theoretical model, we carry out linear regression analysis of linear sizes and luminosities of the 31 radio galaxies in our sample. The regression result shows an empirical relationship between observed source linear sizes and observed bolometric luminosities with slight correlation whose coefficient is 4. Moreover, we compared the theoretical relation and the empirical relation to obtain an estimate of the particle number density of the IGM in which the radio sources are located. The estimate indicates a very low density when compared with the values obtained by authors for the compact steep spectrum (CSS) sources. However, while the linear sizes of the CSS sources are of sub-galactic dimensions – they are located within the interstellar media (ISM) – those of the large extended extragalactic radio sources extend into the intergalactic media. The staggering difference in these obtained densities of IGM and ISM simply supports the general notion that there is a sharp decrease in the density at the IGM–ISM interface. Moreover, using theoretical approach again, with some plausible assumptions, we find another relation which may be used to estimate the ages of extragalactic radio sources. This relation possibly indicates that the plot of the linear sizes against velocities of the radio jets of the radio sources in our sample should be able to supply us with an estimate of the sources’ ages. We find from the simple linear regression analysis of linear sizes and jet velocities of the radio galaxies in our sample, a relation with correlation coefficient given as 4, which is also marginal. Comparison of the two relations – theoretical and the empirical relations – gives an estimate of age of the radio sources. The value shows that these large extended extragalactic radio sources are old sources.

Evolution of the truncated disc and inner hot-flow of GX 339–4

Aims. We study the changes in geometry of the truncated disc and the inner hot-flow of GX 339–4 by analysing the power spectral density (PSD) extracted from six XMM-Newton observations taken at the very end of an outburst. Methods. We developed a theoretical model of the PSD of GX 339–4 in the 0.3–0.7 keV (thermal-reverberation-dominated) and 0.7–1.5 keV (disc-continuum-dominated) energy bands. The model assumes the standard accretion disc to be truncated at a specific radius, inside of which are two distinct hot-flow zones: one spectrally soft and the other spectrally hard. The effects of disc fluctuations and thermal reverberation are taken into account. Results. This model successfully produces the traditional bumpy PSD profiles and provides good fits to the GX 339–4 data. The truncation radius is found to increase from rtrc ∼ 10 to 55rg as the source luminosity decreases, confirming that the truncation radius can be characterized as a function of luminosity. Keeping in mind the large uncertainty in previous measurements of the truncation radius, our values are larger than some obtained from spectroscopic analysis, but smaller than those implied by reverberation lag analysis. Furthermore, the size of two inner hot-flow zones that are spectrally hard and spectrally soft also increases from ∼5 to 27rg and from ∼3 to 26rg, respectively, as the flux decreases. We find that the radial range of the inner hard zone is always larger than the range of the soft hot-flow zone, but by a comparatively small factor of ∼1.1–2.2

Probing episodic accretion with chemistry: CALYPSO observations of the very-low-luminosity Class 0 protostar IRAM 04191+1522

Context. The process of mass accretion in the earliest phases of star formation is still not fully understood: does the accretion rate smoothly decline with the age of the protostar or are there short, intermittent accretion bursts? The latter option would also yield the possibility for very low-luminosity objects (VeLLOs) to be precursors of solar-type stars, even though they do not seem to have sufficiently high accretion rates to reach stellar masses during their protostellar lifetime. Nevertheless, probing such intermittent events in the deeply embedded phase is not easy. Chemical signatures in the protostellar envelope can trace a past accretion burst. Aims. We aim to explore whether or not the observed C18O and N2H+ emission pattern towards the VeLLO IRAM 04191+1522 can be understood in the framework of a scenario where the emission is chemically tracing a past accretion burst. Methods. We used high-angular-resolution Plateau de Bure Interferometer (PdBI) observations of C18O and N2H+ towards IRAM 04191+1522 that were obtained as part of the CALYPSO IRAM Large Program (Continuum And Lines in Young ProtoStellar Objects). We model these observations using a chemical code with a time-dependent physical structure coupled with a radiative transfer module, where we allow for variations in the source luminosity. Results. We find that the N2H+ line emission shows a central hole, with the N2H+ emission peaking at a radius of about 10′′ (1400 au) from the source, while the C18O emission is compact (1.3′′ FWHM, corresponding to 182 au). The morphology of these two lines cannot be reproduced with a constant luminosity model based on the present-day internal luminosity (0.08 L⊙). However, the N2H+ peaks are consistent with a constant-luminosity model of 12 L⊙. Using a model with time-dependent temperature and density profiles, we show that the observed N2H+ peak emission could indeed be caused by a past accretion burst with a luminosity 150 times higher than the present-day luminosity. Such a burst should have occurred a couple of hundred years ago. Conclusions. We suggest that an accretion burst occurred in IRAM 04191+1522 in the recent past. If such bursts are common and sufficiently long in VeLLOs, they could lead to higher accretion onto the central object than their luminosity suggests. For IRAM 04191 in particular, our results yield an estimated final mass of 0.2–0.25 M⊙ by the end of the Class 0 phase, which would make this object a low-mass star rather than a brown dwarf. More generally, our analysis demonstrates that the combination of observations of N2H+ and C18O is a more reliable diagnostic of past outburst activity than C18O or N2H+ emission alone.

FIRST Winged Radio Galaxies with X and Z Symmetry

Abstract X-shaped radio galaxies (XRGs) are a subclass of radio sources that exhibit a pair of secondary low surface brightness radio lobes oriented at an angle to the primary high surface brightness lobes. Sometimes, the secondary low brightened lobes emerge from the edges of the primary high brightened lobes and form a Z-symmetric morphology. We present a systematical search result for XRGs and Z-shaped radio galaxies (ZRGs) from the VLA Faint Images of the Radio Sky at Twenty-Centimeters (VLA FIRST) Survey at 1.4 GHz. Our search yields a total of 296 number of radio sources, out of which 161 are XRGs and 135 are ZRGs. We have also made optical identification of these sources from the different available literature. J1124 + 4325 and J1319 + 0502 are the farthest known XRG and ZRG, respectively. We have estimated spectral index and radio luminosity of these radio sources and made a comparative study with previously detected XRGs and ZRGs. The average value of luminosities for XRGs is higher than that of ZRGs. With the help of a large sample size of the newly discovered XRGs and ZRGs, various statistical properties of these sources are studied. Out of 161 XRGs presented in the current paper, 70% (113) are FR II radio galaxies and 13% (20) are FR I radio galaxies. For 28 XRGs, the morphology is complex and could not be classified. For XRGs, the statistical studies are done on the angle between the major axis and minor axis and the relative size of the major and minor axes. For the ZRGs a statistical study is done on the angular size.

Photoionization feedback in turbulent molecular clouds

ABSTRACT We present a study of the impact of photoionization feedback from young massive stars on the turbulent statistics of star-forming molecular clouds. This feedback is expected to alter the density structure of molecular clouds and affect future star formation. Using the AMUN-Rad code, we first generate a converged isothermal forced turbulent density structure inside a periodic box. We then insert an ionizing source in this box and inject photoionization energy using a two-temperature pseudo-isothermal equation of state. We study the impact of sources at different locations in the box and of different source luminosities. We find that photoionization has a minor impact on the 2D and 3D statistics of turbulence when turbulence continues to be driven in the presence of a photoionizing source. Photoionization is only able to disrupt the cloud if the turbulence is allowed to decay. In the former scenario, the presence of an H ii region inside our model cloud does not lead to a significant impact on observable quantities, independent of the source parameters.

Constraining photohadronic scenarios for the unified origin of IceCube neutrinos and ultrahigh-energy cosmic rays

The diffuse neutrino flux measured in IceCube is comparable with the ultrahigh-energy cosmic ray (UHECR) flux, which has led to the concept of a unified origin of high-energy neutrino and UHECR backgrounds. We construct a generic unification model of sources to explain UHECR data at $\gtrsim 10^{19}$ eV, and high-energy neutrinos with energies that exceed $\sim100$ TeV in the framework of photo-meson production processes, and provide general constraints on the source properties. A source environment with moderately efficient in-situ production of $\gtrsim 100$ TeV neutrinos with an optical depth of $0.1 \lesssim \tau_{p\gamma}\lesssim 0.6$ must be realized to accelerate cosmic rays to ultrahigh energies. The measured fluxes of cosmic rays and neutrinos set a bound on the source luminosity and its rate density. Although the results are rather general and applicable to unknown source population, among the proposed source candidates, low-luminosity gamma-ray bursts (GRBs) and tidal disruption events (TDEs) could satisfy the requirements if the Lorentz bulk factor of plasma outflow and the equipartition parameters for cosmic rays and magnetic field are appropriately selected.

AstroSat observations of the first Galactic ULX pulsar Swift J0243.6+6124

ABSTRACT Swift J0243.6+6124, the first Galactic ultraluminous X-ray pulsar, was observed during its 2017–2018 outburst with AstroSat at both sub- and super-Eddington levels of accretion with X-ray luminosities of LX ∼ 7 × 1037 and 6 × 1038 erg s−1, respectively. Our broad-band timing and spectral observations show that X-ray pulsations at ${\sim}9.85~\rm {s}$ have been detected up to 150 keV when the source was accreting at the super-Eddington level. The pulse profiles are a strong function of both energy and source luminosity, showing a double-peaked profile with pulse fraction increasing from ∼$10{{{\ \rm per\ cent}}}$ at $1.65~\rm {keV}$ to 40–80 ${{\ \rm per\ cent}}$ at $70~\rm {keV}$. The continuum X-ray spectra are well modelled with a high-energy cut-off power law (Γ ∼ 0.6–0.7) and one or two blackbody components with temperatures of ∼0.35 and $1.2~\rm {keV}$, depending on the accretion level. No iron line emission is observed at sub-Eddington level, while a broad emission feature at around 6.9 keV is observed at the super-Eddington level, along with a blackbody radius ($121\!-\!142~\rm {km}$) that indicates the presence of optically thick outflows.