Luminosity Dependence Research Articles

The diverse star formation histories of early massive, quenched galaxies in modern galaxy formation simulations

Abstract We present a comprehensive study of the star formation histories of massive-quenched galaxies at z = 3 in three semi-analytic models (Shark, gaea, Galform) and three cosmological hydrodynamical simulations (Eagle, IllustrisTNG, Simba). We study the predicted number density and stellar mass function of massive-quenched galaxies, their formation and quenching timescales and star-formation properties of their progenitors. Predictions are disparate in all these diagnostics, for instance: (i) some simulations reproduce the observed number density of very massive-quenched galaxies ($>10^{11}\, \rm {\rm M}_{\odot }$) but underpredict the high density of intermediate-mass ones, while others fit well the lower masses but underpredict the higher ones; (ii) In most simulations, except for gaea and eagle, most massive-quenched galaxies had starburst periods, with the most intense ones happening at 4 < z < 5; however, only in Shark and IllustrisTNG we do find a large number of progenitors with star formation rates $>300\rm \, {\rm M}_{\odot }\, yr^{-1}$; (iii) quenching timescales are in the range ≈20 − 150 Myr depending on the simulation; among other differences. These disparate predictions can be tied to the adopted Active Galactic Nuclei (AGN) feedback model. For instance, the explicit black-hole (BH) mass dependence to trigger the ‘radio mode’ in IllustrisTNG and Simba makes it difficult to produce quenched galaxies with intermediate stellar masses, also leading to higher baryon collapse efficiencies (≈15 − 30 per cent); while the strong bolometric luminosity dependence of the AGN outflow rate in gaea leads to BHs of modest mass quenching galaxies. Current observations are unable to distinguish between these different predictions due to the small sample sizes. However, these predictions are testable with current facilities and upcoming observations, allowing a ‘true physics experiment’ to be carried out.

VLTI/GRAVITY interferometric measurements of the innermost dust structure sizes around active galactic nuclei

We present new Very Large Telescope Interferometer (VLTI)/GRAVITY near-infrared interferometric measurements of the angular size of the innermost hot dust continuum for 14 type 1 active galactic nuclei (AGNs). The angular sizes are resolved on scales of ∼0.7 mas and the inferred ring radii range from 0.028 to 1.33 pc, comparable to those reported previously and a factor of 10−20 smaller than the mid-infrared sizes in the literature. Combining our new data with previously published values, we compiled a sample of 25 AGNs with bolometric luminosity ranging from 1042 to 1047 erg s−1, with which we studied the radius-luminosity (R − L) relation for the hot dust structure. Our interferometric measurements of radius are offset by a factor of 2 from the equivalent relation derived through reverberation mapping. Using a simple model to explore the dust structure’s geometry, we conclude that this offset can be explained if the 2 μm emitting surface has a concave shape. Our data show that the slope of the relation is in line with the canonical R ∝ L0.5 when using an appropriately non-linear correction for bolometric luminosity. In contrast, using optical luminosity or applying a constant bolometric correction to it results in a significant deviation in the slope, suggesting a potential luminosity dependence on the spectral energy distribution. Over four orders of magnitude in luminosity, the intrinsic scatter around the R − L relation is 0.2 dex, suggesting a tight correlation between the innermost hot dust structure size and the AGN luminosity.

Luminosity dependence of the multiple cyclotron lines in 4U 0115+63

Context. The Be X-ray binary 4U 0115+63 underwent a giant outburst in 2023 with the X-ray luminosity of the source reaching 1038 erg s−1. During the outburst, two target of opportunity observations were made with NuSTAR. Aims. The main goal of this work is to model the multiple cyclotron scattering features (CRSFs) present in 4U 0115+63 and study their dependence on the luminosity of the source. Methods. The 3–79 keV broadband X-ray coverage of NuSTAR allowed us to properly model the continuum and investigate the nature of the multiple cyclotron resonance scattering features present in the source spectrum. We used the epoch-folding technique to find the pulsation from the source and also studied the changes in the cyclotron line energy with an order of magnitude variation in the source luminosity. Results. We detected five cyclotron lines during the 2023 outburst near 12, 16, 24, 34, and 47 keV. The ∼16 keV cyclotron line cannot be harmonically related to the other cyclotron lines at ∼12 keV and ∼24 keV. This indicates the presence of two fundamental lines in the spectrum of 4U 0115+63 at 12 keV and 16 keV. Conclusions. With the inclusion of the two latest NuSTAR observations, we have expanded the data set of the CRSF line center to encompass a broad range of luminosity. This enables us to comprehensively investigate the relationship between the centroid energy of the cyclotron lines and luminosity. The CRSF line center shows no anticorrelation with luminosity, unlike previously reported. Instead, a weak positive correlation is found in four out of the five detected cyclotron lines of 4U 0115+63. The luminosity variation of the two fundamental CRSFs could be well explained with the prediction from the collisionless shock model. A tentative negative correlation was observed in the fundamental CRSF at 16 keV and its harmonics beyond a “critical luminosity” of 1038 erg/s. This behavior was not present for the 12 keV fundamental CRSF and its harmonic at 24 keV.

The star-forming and ionizing properties of dwarf z ~ 6–9 galaxies in JADES: insights on bursty star formation and ionized bubble growth

ABSTRACT Reionization is thought to be driven by faint star-forming galaxies, but characterizing this population has long remained very challenging. Here, we utilize deep nine-band JADES (JWST Advanced Deep Extragalactic Survey)/NIRCam (Near-Infrared Camera) imaging to study the star-forming and ionizing properties of 756 $z\sim 6-9$ galaxies, including hundreds of very ultraviolet (UV)-faint objects ($M_\mathrm{UV}\gt -18$). The faintest ($m\sim 30$) galaxies in our sample typically have stellar masses of $M_\ast \sim (1-3)\times 10^7\ \mathrm{ M}_\odot$ and young light-weighted ages ($\sim$50 Myr), though some show strong Balmer breaks implying much older ages ($\sim$500 Myr). We find no evidence for extremely massive galaxies ($\gt 3\times 10^{10}\ \mathrm{ M}_\odot$) in our sample. We infer a strong (factor $\gt $2) decline in the typical [O iii]$+$H $\beta$ equivalent widths (EWs) towards very faint $z\sim 6-9$ galaxies, yet a weak UV luminosity dependence on the H $\alpha$ EWs at $z\sim 6$. We demonstrate that these EW trends can be explained if fainter galaxies have systematically lower metallicities as well as more recently declining star formation histories relative to the most UV-luminous galaxies. Our data provide evidence that the brightest galaxies are frequently experiencing a recent strong upturn in star formation rate. We also discuss how the EW trends may be influenced by a strong correlation between $M_\mathrm{UV}$ and Lyman continuum escape fraction. This alternative explanation has dramatically different implications for the contribution of galaxies along the luminosity function to cosmic reionization. Finally, we quantify the photometric overdensities around two $z\,\gt\,7$ strong Ly $\alpha$ emitters. One Ly $\alpha$ emitter lies close to a strong photometric overdensity, while the other shows no significant nearby overdensity, perhaps implying that not all strong $z\,\gt\, 7$ Ly $\alpha$ emitters reside in large ionized bubbles.

Search for radio halos in starburst galaxies

Context. Starburst galaxies are undergoing intense episodes of star formation. In these galaxies, gas is ejected into the surrounding environment through winds created by the effect of hot stars and supernova explosions. When interacting with the intergalactic medium, these winds can produce strong shocks capable of accelerating cosmic rays. The radiation from these cosmic rays mainly occurs in radio and gamma rays. The radio halo can be characterized using the scale height, which is an important parameter for understanding cosmic ray acceleration and transport. Aims. We searched for the presence of radio halos in a sample of edge-on starburst galaxies gathered from the MeerKAT 1.28 GHz Atlas of Southern Sources in the IRAS Revised Bright Galaxy Sample. The investigation of how the radio halos relate to the global properties of the galaxies can shed light on the understanding of the halo origin and the underlying cosmic ray population. Methods. We selected a sample of 25 galaxies with inclinations i > 80° from the original sample and modeled their disk and halo contributions. We determined the scale heights and the radio luminosity of the halos when detected. Results. We have detected and characterized 11 radio halos from a sample of 25 edge-on galaxies. Seven of them are reported here for the first time. The average radio scale height is ∼1 kpc. We found that the halo scale heights increase linearly with the radio diameters and this relation does not depend on the star formation rate. All galaxies in our sample follow the radio-infrared relation with a q parameter value of 2.5 ± 0.1. The halo luminosity linearly increases with the infrared luminosity and star formation rate. Conclusions. The dependence of the halo luminosity on the star formation rate and the infrared luminosity supports the hypothesis that the radio halos are the result of synchrotron radiation produced by relativistic electrons and points toward the fact that the star formation activity plays a crucial role in halo creation. The average scale height of 1 kpc implies a dynamical range of 4 Myr, several orders of magnitude greater than the synchrotron losses for electrons of 10 TeV. This suggests that some process must exist to reaccelerate cosmic rays in the halo if gamma-ray emission of a leptonic origin is detected from the halo. According to the relation between the radio and gamma-ray luminosities, we found that NGC 4666 is a potential gamma-ray source for future observations.

Different Dependence of Narrow Hα Line Luminosity on Optical Continuum Luminosity between Star-forming Galaxies and Type 2 Active Galactic Nuclei: Globally Negative AGN Feedback in Local Type 2 AGN?

In this manuscript, clues are provided to support globally negative active galactic nuclei (AGN) feedback on star formation in the host galaxies of the local low-redshift Sloan Digital Sky Survey (SDSS) Type 2 AGN, based on the different dependence of narrow Hα line luminosity L Hα on optical continuum luminosity λ L cont between star-forming galaxies and Type 2 AGN. Through the measured L Hα and λ L cont in SDSS star-forming galaxies, there is a strong linear correlation between λ L cont and L Hα , accepted as a standard correlation without the effects of AGN activity. Meanwhile, considering the apparent contribution of AGN activity to the narrow Hα line emissions in Type 2 AGN, the correlation between λ L cont and L Hα in the SDSS Type 2 AGN leads to a statistically lower L Hα in Type 2 AGN than in star-forming galaxies, with a significance level higher than 5σ, even after considering necessary effects (including effects of host galaxy properties), leading to the accepted conclusion on the globally negative AGN feedback in the local Type 2 AGN. Meanwhile, the properties of Dn(4000) and Hδ A can provide indirect clues to support the globally negative AGN feedback in local Type 2 AGN, due to older stellar ages in Type 2 AGN. Moreover, it is interesting to expect more than 50% narrow Hα emissions globally suppressed in the host galaxies of Type 2 AGN relative to the star-forming galaxies. The results not only support globally negative AGN feedback in local Type 2 AGN, but also show further clues on the quantification of suppressions of star formation by the globally negative AGN feedback.

Simulating emission line galaxies for the next generation of large-scale structure surveys

ABSTRACT We investigate emission line galaxies across cosmic time by combining the modified L-Galaxies semi-analytical galaxy formation model with the JiuTian cosmological simulation. We improve the tidal disruption model of satellite galaxies in L-Galaxies to address the time dependence problem. We utilize the public code cloudy to compute emission line ratios for a grid of H ii region models. The emission line models assume the same initial mass function as that used to generate the spectral energy distribution of semi-analytical galaxies, ensuring a coherent treatment for modelling the full galaxy spectrum. By incorporating these emission line ratios with galaxy properties, we reproduce observed luminosity functions for H α, H β, [O ii], and [O iii] in the local Universe and at high redshifts. We also find good agreement between model predictions and observations for autocorrelation and cross-correlation functions of [O ii]-selected galaxies, as well as their luminosity dependence. The bias of emission line galaxies depends on both luminosity and redshift. At lower redshifts, it remains constant with increasing luminosity up to around $\sim 10^{42.5}\, {\rm erg\, s^{-1}}$ and then rises steeply for higher luminosities. The transition luminosity increases with redshift and becomes insignificant above z = 1.5. Generally, galaxy bias shows an increasing trend with redshift. However, for luminous galaxies, the bias is higher at low redshifts, as the strong luminosity dependence observed at low redshifts diminishes at higher redshifts. We provide a fitting formula for the bias of emission line galaxies as a function of luminosity and redshift, which can be utilized for large-scale structure studies with future galaxy surveys.

Cyclotron line evolution revealed with pulse-to-pulse analysis in the 2020 outburst of 1A 0535+262

ABSTRACT We present a detailed analysis of the X-ray luminosity (LX) dependence of the cyclotron absorption line energy (Ecyc) for the X-ray binary pulsar 1A 0535+262 during its 2020 giant outburst based on pulse-to-pulse analysis. By applying this technique to high cadence observations of Insight-Hard X-ray Modulation Telescope, we reveal the most comprehensive Ecyc–LX correlation across a broad luminosity range of ∼(0.03–1.3) × 1038 erg s−1. Apart from the positive and negative correlations between cyclotron line energy and luminosity at LX ∼ (1–3) × 1037 erg s−1 and ∼(7–13) × 1037 erg s−1, which are expected from the typical subcritical and supercritical accretion regimes, respectively, a plateau in the correlation is also detected at ∼(3–7) × 1037 erg s−1. Moreover, at the lowest luminosity level (LX ≲ 1037 erg s−1), the positive Ecyc–LX correlation seems to be broken, and the pulse profile also occurs a significant transition. These discoveries provide the first complete view on the correlation between luminosity and the centriod energy of the cyclotron line, and therefore are relevant for understanding how accretion on to magnetized neutron stars depends on luminosity.

Multiwavelength Spectral Energy Distribution Analysis of X-Ray Selected Active Galactic Nuclei at z = 0.2–0.8 in the Stripe 82 Region

We perform a systematic, multiwavelength spectral energy distribution (SED) analysis of X-ray detected active galactic nuclei (AGNs) at z = 0.2–0.8 with Sloan Digital Sky Survey (SDSS) counterparts in the Stripe 82 region, consisting of 60 type 1 and 137 type 2 AGNs covering a 2–10 keV luminosity range of 41.6<logLx<44.7 . The latest CIGALE code, where dusty polar components are included, is employed. To obtain reliable host and AGN parameters in type 1 AGNs, we utilize the image-decomposed optical SEDs of host galaxies by Li et al. based on the Subaru Hyper-Suprime Cam images. The mean ratio of black hole masses (M BH) and stellar masses (M stellar) of our X-ray detected type 1 AGN sample, log(MBH/Mstellar)=−2.7±0.5 , is close to the local relation between BH and stellar masses as reported by Li et al. for SDSS quasars. This ratio is slightly lower than that found for more luminous ( logLbol>45 ) type 1 AGNs at z ∼ 1.5. This can be explained by the AGN luminosity dependence of log(MBH/Mstellar) , which little evolves with redshift. We confirm the trend that the UV-to-X-ray slope (α OX) or X-ray-to-bolometric correction factor (κ 2–10) increases with AGN luminosity or Eddington ratio. We find that type 1 and type 2 AGNs with the same luminosity ranges share similar host stellar mass distributions, whereas type 2s tend to show smaller AGN luminosities than type 1s. This supports the luminosity-dependent (or Eddington-ratio-dependent) unified scheme.

Bunch lifetime analysis based on the injection interval at SuperKEKB

We study the bunch lifetime dependence on the bunch position at the SuperKEKB main ring. The bunch lifetime ratio is evaluated from the injection interval. The lifetime imbalance is observed between the forward and backward of the bunch train at the positron ring. It becomes larger when the beam current becomes larger. The luminosity dependence on the RF-bucket spacing can be explained with the lifetime dependency which we determined in this analysis. The result indicates the existence of the unknown phenomenon which varies the Touschek lifetime with the RF-bucket spacing in the electron ring.

Progenitor-age dependence of type Ia supernova luminosity and its cosmological implication

Progenitor-age dependence of type Ia supernova luminosity and its cosmological implication

Reconsidering photometric estimation of local star formation environment and its correlation with Type Ia supernova luminosity

ABSTRACT Recent studies on the environmental dependence of Type Ia supernova (SN Ia) luminosity focus on the local environment where the SN exploded, considering that this is more directly linked to the SN progenitors. However, there is a debate about the local environmental, specifically local star formation rate (SFR), dependence of the SN Ia luminosity. A recent study claims that the dependence is insignificant (0.051 ± 0.020 mag; 2.6σ), based on the local SFR measurement by fitting local ugrizy photometry data. However, we find that this photometric local SFR measurement is inaccurate. We argue this based on the theoretical background of SFR measurement and the methodology used to make that claim with their local ugrizy photometry data, especially due to a limited range of extinction parameters used when fitting the data. Therefore, we reanalyse the same host galaxies with the same fitting code, but with more physically motivated extinction treatments and global ugriz photometry of host galaxies. We estimate global stellar mass and SFR. Then, local star formation environments are inferred by using the method which showed that SNe Ia in globally passive galaxies have locally passive environments, while those in globally star-forming low-mass galaxies have locally star-forming environments. We find that there is significant local environmental dependence of SN Ia luminosities: SNe Ia in locally star-forming environments are 0.072 ± 0.021 mag (3.4σ) fainter than those in locally passive environments, even though SN Ia luminosities have been further corrected by the BEAMS with Bias Corrections method that reduces the size of the dependence.

Quasar Factor Analysis—An Unsupervised and Probabilistic Quasar Continuum Prediction Algorithm with Latent Factor Analysis

Since their first discovery, quasars have been essential probes of the distant Universe. However, due to our limited knowledge of its nature, predicting the intrinsic quasar continua has bottlenecked their usage. Existing methods of quasar continuum recovery often rely on a limited number of high-quality quasar spectra, which might not capture the full diversity of the quasar population. In this study, we propose an unsupervised probabilistic model, quasar factor analysis (QFA), which combines factor analysis with physical priors of the intergalactic medium to overcome these limitations. QFA captures the posterior distribution of quasar continua through generatively modeling quasar spectra. We demonstrate that QFA can achieve the state-of-the-art performance, ∼2% relative error, for continuum prediction in the Lyα forest region compared to previous methods. We further fit 90,678 2 < z < 3.5, signal-to-noise ratio >2 quasar spectra from Sloan Digital Sky Survey Data Release 16 and found that for ∼30% quasar spectra where the continua were ill-determined with previous methods, QFA yields visually more plausible continua. QFA also attains ≲1% error in the 1D Lyα power spectrum measurements at z ∼ 3 and ∼4% in z ∼ 2.4. In addition, QFA determines latent factors representing more physical motivation than principal component analysis. We investigate the evolution of the latent factors and report no significant redshift or luminosity dependency except for the Baldwin effect. The generative nature of QFA also enables outlier detection robustly; we showed that QFA is effective in selecting outlying quasar spectra, including damped Lyα systems and potential Type II quasar spectra.

X-ray emission from Ae/Be Herbig stars due to disc–stellar magnetosphere interaction

ABSTRACT We reanalyse archival X-ray data of 16 Ae/Be Herbig stars obtained by the XMM–Newton and Chandra satellites. Stellar X-ray spectra in the energy range 0.2–8 keV were fitted with the use of APEC and MEKAL hot plasma emission models, and with models with an additional power-law component. We find that for Herbig stars, the dependence of the unabsorbed X-ray luminosity on stellar mass and radius, LX ∝ RαMβ with α ≈ 3 and β ≈ 2, is similar to that for T Tauri stars. The independently determined accretion rates, rotation periods, and the surface magnetic fields follow a tight correlation predicted by the standard magnetospheric accretion theory. We suggest that X-ray emission from Herbig stars is powered by magnetic reconnection events in the tenuous corona at the disc–magnetosphere boundary.

A Selection of Hα Emitters at z = 2.1–2.5 Using the K s -band Photometry of ZFOURGE

Large and less-biased samples of star-forming galaxies are essential to investigate galaxy evolution. Hα emission line is one of the most reliable tracers of star-forming galaxies because its strength is directly related to recent star formation. However, it is observationally expensive to construct large samples of Hα emitters by spectroscopic or narrowband imaging survey at high redshifts. In this work, we demonstrate a method to extract Hα fluxes of galaxies at z = 2.1–2.5 from K s broadband photometry of ZFOURGE catalog. Combined with 25–39 other filters, we estimate the emission-line fluxes by SED fitting with stellar population models that incorporate emission-line strengths. 2005 galaxies are selected as Hα emitters by our method and their fluxes show good agreement with previous measurements in the literature. On the other hand, there are more Hα luminous galaxies than previously reported. The discrepancy can be explained by extended Hα profiles of massive galaxies and a luminosity dependence of dust attenuation, which are not taken into account in the previous work. We also find that there are a large number of low-mass galaxies with a much higher specific star formation rate (sSFR) than expected from the extrapolated star formation main sequence. Such low-mass galaxies exhibit larger ratios between Hα and UV fluxes compared to more massive high sSFR galaxies. This result implies that a “starburst” mode may differ among galaxies: low-mass galaxies appear to assemble their stellar mass via short-duration bursts, while more massive galaxies tend to experience longer-duration (>10 Myr) bursts.

The non-thermal secondary CMB anisotropies from a cosmic distribution of radio galaxy lobes

ABSTRACT Current and upcoming high angular resolution and multifrequency experiments are well poised to explore the rich landscape of secondary cosmic microwave background (CMB) anisotropies. In this context, we compute, for the first time, the power spectrum of CMB fluctuations from a cosmological distribution of evolving lobes of giant radio galaxies. We also explicitly take into account the non-thermal electron distribution, which has important implications for the inference of the CMB angular power spectrum. We calculate the mean global non-thermal y-distortion, 〈y〉NT. For observationally reasonable distribution of the jet luminosities in the range of 1045–1047 erg s−1, we find 〈y〉NT to be less than 10−5, and hence not violating the cosmic background explorer limit as previously claimed. Using the unique spectral dependence of the non-thermal Sunyaev–Zeldovich (SZ), we show that a detection of 〈y〉NT can be within reach at the level of ≳5σ from a future Primordial Inflation Explorer (PIXIE)-like experiment provided we understand the foregrounds precisely. The total non-thermal SZ power spectrum, $C^{\mathrm{ NT}}_\ell$, from the radio lobes peaks at ℓ ∼ 3000 with an amplitude $\sim 1{{\ \rm per\ cent}}$ of thermal SZ power spectrum from galaxy clusters. A detection of the $C^{\mathrm{ NT}}_\ell$, with a PIXIE-like sensitivity experiment, can lead to ∼5σ constraint on the mass dependence of the jet luminosity with the constraint becoming at least ten times better for the proposed more ambitious CMB-HD survey. This will further lead to the tightest constraint on the central black hole mass-to-host halo mass scaling relations.

First eROSITA study of nearby M dwarfs and the rotation-activity relation in combination with TESS

We present the first results with the ROentgen Survey with an Imaging Telescope Array (eROSITA) on board the Russian Spektrum-Roentgen-Gamma mission, and we combine the new X-ray data with observations with the Transiting Exoplanet Survey Satellite (TESS). We used theSUPERBLINKproper motion catalog of nearby M dwarfs as input sample to search for eROSITA and TESS data. We extractedGaiaDR2 data for the full M dwarf catalog, which comprises ~9000 stars, and we calculated the stellar parameters from empirical relations with optical/IR colors. Then we cross-matched this catalog with the eROSITA Final Equatorial Depth Survey (eFEDS) and the first eROSITA all-sky survey (eRASS1). After a meticulous source identification in which we associated the closestGaiasource with the eROSITA X-ray detections, our sample of M dwarfs is defined by 687 stars with SpT = K5..M7 (673 from eRASS1 and 14 from eFEDS). While for eRASSl we used the data from the source catalog provided by the eROSITA_DE consortium, for the much smaller eFEDS sample, we performed the data extraction, and we analyzed the X-ray spectra and light curves. This unprecedented data base for X-ray emitting M dwarfs allowed us to place a quantitative constraint on the mass dependence of the X-ray luminosity, and to determine the change in the activity level with respect to pre-main-sequence stars. TESS observations are available for 489 of 687 X-ray detected M dwarfs. By applying standard period search methods, we were able to determine the rotation period for 180 X-ray detected M dwarfs. This is about one-forth of the X-ray sample. With the joint eROSITA and TESS sample, and combining it with our compilation of historical X-ray and rotation data for M dwarfs, we examined the mass dependence of the saturated regime of the rotation-activity relation. A first comparison of eROSITA hardness ratios and spectra shows that 65% of the X-ray detected M dwarfs have coronal temperatures of ~0.5 keV. We performed a statistical investigation of the long-term X-ray variability of M dwarfs by comparing the eROSITA measurements to those obtained ~30 yr earlier during the ROSAT all-sky survey (RASS). Evidence for X-ray flares is found in various parts of our analysis: directly from an inspection of the eFEDS light curves, in the relation between RASS and eRASSl X-ray luminosities, and in a subset of stars that displays hotter X-ray emission than the bulk of the sample according to the hardness ratios. Finally, we point out the need to obtain X-ray spectroscopy for more M dwarfs to study the coronal temperature-luminosity relation, which is not well constrained by our eFEDS results.

MOBSTER – VI. The crucial influence of rotation on the radio magnetospheres of hot stars

ABSTRACT Numerous magnetic hot stars exhibit gyrosynchrotron radio emission. The source electrons were previously thought to be accelerated to relativistic velocities in the current sheet formed in the middle magnetosphere by the wind opening magnetic field lines. However, a lack of dependence of radio luminosity on the wind power, and a strong dependence on rotation, has recently challenged this paradigm. We have collected all radio measurements of magnetic early-type stars available in the literature. When constraints on the magnetic field and/or the rotational period are not available, we have determined these using previously unpublished spectropolarimetric and photometric data. The result is the largest sample of magnetic stars with radio observations that has yet been analysed: 131 stars with rotational and magnetic constraints, of which 50 are radio-bright. We confirm an obvious dependence of gyrosynchrotron radiation on rotation, and furthermore find that accounting for rotation neatly separates stars with and without detected radio emission. There is a close correlation between H α emission strength and radio luminosity. These factors suggest that radio emission may be explained by the same mechanism responsible for H α emission from centrifugal magnetospheres, i.e. centrifugal breakout (CBO), however, while the H α-emitting magnetosphere probes the cool plasma before breakout, radio emission is a consequence of electrons accelerated in centrifugally driven magnetic reconnection.

SOFIA/FORCAST Monitoring of the Dust Emission from R Aqr: Start of the Eclipse

Abstract We present mid-infrared (mid-IR) spectra from our continued monitoring of R Aquarii, the nearest symbiotic Mira, using the Stratospheric Observatory for Infrared Astronomy (SOFIA). New photometric and spectroscipic data were obtained with the Faint Object infraRed CAmera for the SOFIA Telescope in 2018 and 2019 after the system had started its “eclipse,” during which it became two magnitudes fainter in the visual. The mid-IR flux, in particular the 10 μm silicate feature, has strengthened compared with the previous cycles. Radiative transfer models for the circumstellar dust emission were calculated for the new spectra, and recalculated for those previously obtained using more appropriate values of the near-IR magnitudes to constrain the properties of the asymptotic giant branch spectra heating the dust. The modeling shows that the luminosity dependence on pulsation phase is not affected by the onset of the eclipse, and that the increase in the mid-IR flux is due to a higher dust density. The models also confirm our earlier results that micron-size grains are present, and that no changes in the grain composition are required to explain the variations in the spectra.

Study of timing and spectral properties of the X-ray pulsar 1A 0535+262 during the giant outburst in 2020 November–December

ABSTRACT We made a detailed study of the timing and spectral properties of the X-ray pulsar 1A 0535+262 during the recent giant outburst in 2020 November and December. The flux of the pulsar reached a record value of ∼12.5 Crab as observed by Swift/BAT (15–50 keV) and the corresponding mass accretion rate was ∼6.67 × 1017 g s−1 near the peak of the outburst. There was a transition from the subcritical to the supercritical accretion regime which allows exploring different properties of the source in the supercritical regime. A q-like feature was detected in the hardness–intensity diagram during the outburst. We observed high variability and strong energy dependence of pulse profiles during the outburst. Cyclotron Resonant Scattering Feature (CRSF) was detected at ∼44 keV from the NuSTAR energy spectrum in the subcritical regime and the corresponding magnetic field was B ≃ 4.9 × 1012 G. The energy of the CRSF was shifted towards lower energy in the supercritical regime. The luminosity dependence of the CRSF was studied and during the supercritical regime, a negative correlation was observed between the line energy and luminosity. The critical luminosity was ∼6 × 1037erg s−1 above which a state transition occurred. A reversal of correlation between the photon index and luminosity was observed near the critical luminosity. The NuSTAR spectra can be described by a composite model with two continuum components, a blackbody emission, cut-off power law, and a discrete component to account for the iron emission line at 6.4 keV. An additional cyclotron absorption feature was included in the model.