ALBERT

Description: The physical properties of the accretion flow and of the X-ray emitting plasma, in supermassive black holes accreting at extreme Eddington rates, are still very unclear. Here we present the analysis of simultaneous XMM-Newton and NuSTAR observations of the hyper-Eddington Seyfert 1 galaxy IRAS 04416+1215, carried out in 2020. The main goal of these observations is to investigate the properties of the X-ray corona, as well as the structure of the accretion flow and of the circumnuclear environment, in this regime of extreme accretion. IRAS 04416+1215 has one of the highest Eddington ratio (λEdd ≃ 472) in the local Universe. It shows an interesting spectral shape, very similar to the standard Narrow Line Seyfert 1 galaxy’s spectra, with the presence of multi-phase absorption structure composed of three phases, whose estimate of the minimum and maximum distances suggests two different interpretations, one consistent with the three X-ray winds being co-spatial, and possibly driven by magnetohydrodynamical processes, the other consistent with the multi-phase winds being also multi-scale. The X-ray spectrum of IRAS 04416+1215 also has a prominent soft excess component and a hard X-ray emission dominated by a reflection component. Moreover, our detailed spectral analysis shows that IRAS 04416+1215 has the lowest coronal temperature measured so far by NuSTAR (kTe = 3 − 22 keV, depending on the model). This is consistent with a hybrid coronal plasma, in which the primary continuum emission is driven by pair production due to high-energy tail of the energy distribution of non-thermal electrons.

Description: We investigate the properties of the dark matter (DM) halo surrounding the nearby galaxy group NGC 1600. Through the use of deep (252 ks) Chandra observations and 64.3 ks of XMM–Newton observations, we construct surface brightness profiles in multiple energy bands in order to perform hydrostatic equilibrium analysis of the hot plasma within NGC 1600. Regardless of the DM model profile assumed, we measure a halo concentration (c200) that is an extreme, positive outlier of the ΛCDM c200-M200 relation. For a typical NFW DM profile, we measure c200 = 26.7 ± 1.4 and M200 = (2.0 ± 0.2) × 1013 M⊙; assuming a similar halo mass, the average concentration expected is c200 = 6 − 7 for the theoretical ΛCDM c-M relation. Such a high concentration is similar to that of well-known fossil groups MRK 1216 and NGC 6482. While NGC 1600 exhibits some properties of a fossil group, it fails to meet the X-ray luminosity threshold of LX 〉 5 × 1041 erg s−1. Whether or not it is considered a fossil group, the high concentration value makes it part of a select group of galaxy groups.

Description: The continuum emission from radio galaxies can be generally classified into different morphological classes such as FRI, FRII, Bent, or Compact. In this paper, we explore the task of radio galaxy classification based on morphology using deep learning methods with a focus on using a small scale dataset (∼2000 samples). We apply few-shot learning techniques based on Twin Networks and transfer learning techniques using a pre-trained DenseNet model with advanced techniques like cyclical learning rate and discriminative learning to train the model rapidly. We achieve a classification accuracy of over 92% using our best performing model with the biggest source of confusion being between Bent and FRII type galaxies. Our results show that focusing on a small but curated dataset along with the use of best practices to train the neural network can lead to good results. Automated classification techniques will be crucial for upcoming surveys with next generation radio telescopes which are expected to detect hundreds of thousands of new radio galaxies in the near future.

Description: Over the past decade, studies of dust in the Andromeda galaxy (M31) have shown radial variations in the dust emissivity index (β). Understanding the astrophysical reasons behind these radial variations may give clues about the chemical composition of dust grains, their physical structure, and the evolution of dust. We use 12CO(J=1-0) observations taken by the Combined Array for Research in Millimeter Astronomy (CARMA) and dust maps derived from Herschel images, both with an angular resolution of 8″ and spatial resolution of 30 pc, to study variations in β across an area of ≈ 18.6 kpc2 in M31. We extract sources, which we identify as molecular clouds, by applying the astrodendro algorithm to the 12CO and dust maps, which as a byproduct allows us to compare continuum emission from dust and CO emission as alternative ways of finding molecular clouds. We then use these catalogues to investigate whether there is evidence that β is different inside and outside molecular clouds. Our results confirm the radial variations of β seen in previous studies. However, we find little difference between the average β inside molecular clouds compared to outside molecular clouds, in disagreement with models which predict an increase of β in dense environments. Finally, we find some clouds traced by dust with very little CO which may be either clouds dominated by atomic gas or clouds of molecular gas that contain little CO.

Description: Baryonic matter can be accreted on to primordial back holes (PBHs) formed in the early Universe. The radiation from accreting PBHs is capable of altering the evolution of the intergalactic medium (IGM), leaving marks on the global 21 cm signal in the dark ages. For accreting PBHs with mass MPBH = 103(104) M⊙ and mass fraction fPBH = 10−1(10−3), the brightness temperature deviation ΔδTb reaches $sim 18~(26)~ m mK$ at redshift z ∼ 90 ($ u sim 16~ m MHz$), and the gradient of the brightness temperature dδTb/dν reaches $sim 0.8~(0.5)~ m mK~MHz^{-1}$ at frequency $ u sim 28~ m MHz$ (z ∼ 50). For larger PBHs with higher mass fraction, the brightness temperature deviation is larger in the redshift range z ∼ 30–300 ($ u sim 5!-!46~ m MHz$), and the gradient is lower at the frequency range $ u sim 20!-!60~ m MHz$ (z ∼ 23–70). It is impossible to detect these low-frequency radio signals from the Earth due to the influence of the Earth’s ionosphere. However, after taking care of the essential factors properly, e.g. the foreground and interference, a future radio telescope in lunar orbit or on the farside surface of the Moon has a chance of detecting the global 21 cm signals impacted by accreting PBHs and distinguishing them from the standard model.

Description: We present VIVACE, the VIrac VAriable Classification Ensemble, a catalogue of variable stars extracted from an automated classification pipeline for the Vista Variables in the Vía Láctea (VVV) infrared survey of the Galactic bar/bulge and southern disc. Our procedure utilises a two-stage hierarchical classifier to first isolate likely variable sources using simple variability summary statistics and training sets of non-variable sources from the Gaia early third data release, and then classify candidate variables using more detailed light curve statistics and training labels primarily from OGLE and VSX. The methodology is applied to point-spread-function photometry for ∼490 million light curves from the VIRAC v2 astrometric and photometric catalogue resulting in a catalogue of ∼1.4 million likely variable stars, of which ∼39, 000 are high-confidence (classification probability 〉0.9) RR Lyrae ab stars, ∼8000 RR Lyrae c/d stars, ∼187, 000 detached/semi-detached eclipsing binaries, ∼18, 000 contact eclipsing binaries, ∼1400 classical Cepheid variables and ∼2200 Type II Cepheid variables. Comparison with OGLE-4 suggests a completeness of around 90  per cent for RRab and ≲ 60 per cent for RRc/d, and a misclassification rate for known RR Lyrae stars of around 1 per cent for the high confidence sample. We close with two science demonstrations of our new VIVACE catalogue: first, a brief investigation of the spatial and kinematic properties of the RR Lyrae stars within the disc/bulge, demonstrating the spatial elongation of bar-bulge RR Lyrae stars is in the same sense as the more metal-rich red giant population whilst having a slower rotation rate of ∼40 km s−1kpc−1; and secondly, an investigation of the Gaia EDR3 parallax zeropoint using contact eclipsing binaries across the Galactic disc plane and bulge.

Description: In this paper, we report white dwarfs identified in the 5th Data Release of the Large Area Multi-Object fibre Spectroscopic Telescope, including spectral types of DA, DB, DC, DZ, and so on. There are 2 625 DA spectra of 2 281 DA stars, 182 DB spectra of 166 DB stars, 62 DC spectra of 58 DC stars, 36 DZ spectra of 33 DZ stars and many other types identified, in addition to our previous paper (Data Release 2). Among those sources, 393 DA stars and 46 DB stars are new identifications after cross-matching with the literature. In order to select DA candidates, we use the classification result from the LAMOST pipeline, colour-colour cut method and a random forest machine learning method. For DBs, since there is no template for DB in the pipeline model, a random forest machine learning method is chosen to select candidates. All the WD candidates have been visually checked individually. The parameters of effective temperature, surface gravity, mass, and cooling age have been estimated for relatively high signal-to-noise ratio DAs and DBs. The peaks of the DA and DB mass distributions are found to be around 0.62 M⊙ and 0.65 M⊙, respectively. Finally, the data and method we used to select white dwarf candidates for the second phase of LAMOST survey are also addressed in this paper.

Description: We use the k-nearest neighbor probability distribution function (kNN-PDF, Banerjee & Abel 2021) to assess convergence in a scale-free N-body simulation. Compared to our previous two-point analysis, the kNN-PDF allows us to quantify our results in the language of haloes and numbers of particles, while also incorporating non-Gaussian information. We find good convergence for 32 particles and greater at densities typical of haloes, while 16 particles and fewer appears unconverged. Halving the softening length extends convergence to higher densities, but not to fewer particles. Our analysis is less sensitive to voids, but we analyse a limited range of underdensities and find evidence for convergence at 16 particles and greater even in sparse voids.

Description: The Enthalpy Based Thermal Evolution of Loops (EBTEL) approximate model for static and dynamic coronal loops is developed to include the effect of a loop cross-sectional area which increases from the base of the transition region (TR) to the corona. The TR is defined as the part of a loop between the top of the chromosphere and the location where thermal conduction changes from an energy loss to an energy gain. There are significant differences from constant area loops due to the manner in which the reduced volume of the TR responds to conductive and enthalpy fluxes from the corona. For static loops with modest area variation the standard picture of loop energy balance is retained, with the corona and TR being primarily a balance between heating and conductive losses in the corona, and downward conduction and radiation to space in the TR. As the area at the loop apex increases, the TR becomes thicker and the density in TR and corona larger. For large apex areas, the coronal energy balance changes to one primarily between heating and radiation, with conduction playing an increasingly unimportant role, and the TR thickness becoming a significant fraction of the loop length. Approximate scaling laws are derived that give agreement with full numerical solutions for the density, but not the temperature. For non-uniform areas, dynamic loops have a higher peak temperature and are denser in the radiative cooling phase by of order 50% than the constant area case for the examples considered. They also show a final rapid cooling and draining once the temperature approaches 1 MK. Although the magnitude of the emission measure will be enhanced in the radiative phase, there is little change in the important observational diagnostic of its temperature dependence.

Description: Magnetars, the likely sources of Fast Radio Bursts (FRBs), produce both steady highly relativistic magnetized winds, and occasional ejection events. We demonstrate that the requirement of conservation of the magnetic flux dominates the overall dynamics of magnetic explosions. This is missed in conventional hydrodynamic models of the ejections as expanding shell with parametrically added magnetic field, as well as one-dimensional models of magnetic disturbances. Magnetic explosions from magnetars come into force balance with the pre-flare wind close to the light cylinder. They are then advected quietly with the wind, or propagate as electromagnetic disturbances. No powerful shock waves are generated in the wind.