Papers for Tuesday, Aug 24 2021

The massive black hole (BH) population in dwarf galaxies ($M_{\rm BH} \lesssim 10^5~M_\odot$) can provide strong constraints on the origin of BH seeds. However, traditional optical searches for active galactic nuclei (AGNs) only reliably detect high-accretion, relatively high-mass BHs in dwarf galaxies with low amounts of star formation, leaving a large portion of the overall BH population in dwarf galaxies relatively unexplored. Here, we present a sample of 81 dwarf galaxies ($M_\star \le 3 \times 10^9~M_\odot$) with detectable [Fe X]$\lambda$6374 coronal line emission indicative of accretion onto massive BHs, only two of which were previously identified as optical AGNs. We analyze optical spectroscopy from the Sloan Digital Sky Survey and find [Fe X]$\lambda$6374 luminosities in the range $L_{\rm [Fe\,X]}\approx10^{36}$-$10^{39}$ erg s$^{-1}$, with a median value of $1.6 \times 10^{38}$ erg s$^{-1}$. The [Fe X]$\lambda$6374 luminosities are generally much too high to be produced by stellar sources, including luminous Type IIn supernovae (SNe). Moreover, based on known SNe rates, we expect at most 8 Type IIn SNe in our sample. On the other hand, the [Fe X]$\lambda$6374 luminosities are consistent with accretion onto massive BHs from AGNs or tidal disruption events (TDEs). We find additional indicators of BH accretion in some cases using other emission line diagnostics, optical variability, X-ray and radio emission (or some combination of these). However, many of the galaxies in our sample only have evidence for a massive BH based on their [Fe X]$\lambda$6374 luminosities. This work highlights the power of coronal line emission to find BHs in dwarf galaxies missed by other selection techniques and to probe the BH population in bluer, lower mass dwarf galaxies.

We present deep Hubble Space Telescope imaging of five faint dwarf galaxies associated with the nearby spiral NGC 253 (D$\approx$3.5 Mpc). Three of these are newly discovered ultra-faint dwarf galaxies, while all five were found in the Panoramic Imaging Survey of Centaurus and Sculptor (PISCeS), a Magellan$+$Megacam survey to identify faint dwarfs and other substructures in resolved stellar light around massive galaxies outside of the Local Group. Our HST data reach $\gtrsim$3 magnitudes below the tip of the red giant branch for each dwarf, allowing us to derive their distances, structural parameters, and luminosities. All five systems contain predominantly old, metal-poor stellar populations (age$\sim$12 Gyr, [M/H]$\lesssim$$-$1.5) and have sizes ($r_{h}$$\sim$110-3000 pc) and luminosities ($M_V$$\sim$$-7$ to $-12$ mag) largely consistent with Local Group dwarfs. The three new NGC 253 satellites are among the faintest systems discovered beyond the Local Group. We also use archival HI data to place limits on the gas content of our discoveries. Deep imaging surveys such as our program around NGC 253 promise to elucidate the faint end of the satellite luminosity function and its scatter across a range of galaxy masses, morphologies, and environments in the decade to come.

High eccentricity tidal migration (HEM) is a promising channel for the origins of hot Jupiters and hot Neptunes. In the typical HEM scenario, a planet forms beyond the ice line, but alternatively a planet can disk migrate or form warm and undergo a short final stretch of HEM. At the warm origin point, general relavistic precession can reduce the amplitude of Kozai-Lidov oscillations driven by an outer companion. We show that warm planets that achieve HEM under these conditions -- and with common types of planetary and stellar companions -- tend to end up with near-polar spin-orbit alignments (psi = 50-130 degrees) instead of concentrated at 40 and 140 degrees. Thus short distance, GR-reduced HEM is a possible explanation for the observed population of perpendicular planets.

Anomalies among the daughter nuclei of the extinct short-lived radionuclides (SLRs) in the calcium-aluminum-rich inclusions (CAIs) indicate that the Solar System must have been born near a source of the SLRs so that they could be incorporated before they decayed away. $\gamma$-rays from one such living SLR, $^{26}$Al, are detected in only a few nearby star-forming regions. Here we employ multi-wavelength observations to demonstrate that one such region, Ophiuchus, containing many pre-stellar cores that may serve as analogs for the emerging Solar System, is inundated with $^{26}$Al from the neighboring Upper-Scorpius association, and so may provide concrete guidance for how SLR enrichment proceeded in the Solar System complementary to the meteoritics. We demonstrate via Bayesian forward modeling drawing on a wide range of observational and theoretical results that this $^{26}$Al likely 1) arises from supernova explosions, 2) arises from multiple stars, 3) has enriched the gas prior to the formation of the cores, and 4) gives rise to a broad distribution of core enrichment spanning about two orders of magnitude. This means that if the spread in CAI ages is small, as it is in the Solar System, protoplanetary disks must suffer a global heating event.

The Cygnus Cocoon is the first gamma-ray superbubble powered by a massive stellar association, the OB2 association. It was postulated that the combined effects of the stellar winds of all the massive O-type stars of the OB2 association can accelerate the cosmic rays to PeV energy in the Cocoon. The conclusive proof of acceleration to PeV energy in the Cocoon will identify the stellar association as a PeV cosmic-ray accelerator, known as PeVatron. However, the Cocoon has been previously studied only up to 10 TeV. In this contribution, using 1343 days of High Altitude Water Cherenkov (HAWC) observatory data, we present the morphological and spectral study of the Cocoon above 1 TeV to beyond 100 TeV. The analysis at higher TeV energies reveals a softer spectrum compared to the GeV gamma-ray observation. This result suggests that the accelerator's efficiency decreases around hundreds of TeV, or after being accelerated, the highest-energy protons escape the region. The study above 10 TeV presented here demonstrates how CR accelerators operate in these extreme energies and how particle transport impacts high-energy emission.

1A 0535+262 is a Be X-ray binary pulsar and one of the only galactic pulsar systems to show radio jet emission. Characterizing the very high energy emission (VHE, >100 GeV) in these extreme microquasars is critical to understanding their contribution to the origin of galactic cosmic rays. The 2020 giant outburst of this system, where X-ray fluxes exceeded 12 Crab, marked a rare opportunity to investigate the gamma-ray and rapid optical variability of these transient systems while in such an extreme state. This month of activity marked the brightest flare measured in this system. VERITAS's developing optical capabilities in tandem with the ability to measure TeV gamma rays allowed for a unique campaign to be undertaken. VERITAS's observations of this system during the outburst will be presented in the context of observations at lower energies and previous observations of this system by imaging atmospheric Cherenkov telescopes.

Fast radio burst (FRBs) are an exciting class of bright, extragalactic, millisecond radio transients. The recent development of large field-of-view (FOV) radio telescopes has caused a rapid rise in the number of identified single burst and repeating FRBs. This has allowed for the extensive multi-wavelength follow-up to search for the potential counterparts predicted by theoretical models. New observations of similar radio transients in Galactic magnetars like SGR 1935+2154 have continued to motivate the search for rapid optical and very-high-energy (VHE, >100 GeV) counterparts. Since 2016 VERITAS has engaged in an FRB observing campaign to search for the prompt optical, and VHE emission from multiple repeating FRBs. We present these new results from VERITAS observations of five repeating sources including data taken simultaneously with bursts observed by the CHIME radio telescope.

We present the results of an investigation of the relation between space-weather parameters and cosmic ray (CR) intensity modulation using algorithm-selected Forbush decreases (FDs) from Moscow (MOSC) and Apatity (APTY) neutron monitor (NM) stations during solar cycle 23. Our FD location program detected 408 and 383 FDs from MOSC and APTY NM stations respectively. A coincident computer code employed in this work, detected 229 FDs that were observed at the same universal Time (UT) at the two stations. Out of the 229 simultaneous FDs, we formed a subset of 139 large FDs(\%) $\leq-4$ at Moscow station. We peformed a two dimensional regression analysis between the FD magnitudes and the space-weather data on the two samples. We find that there were significant space-weather disturbances at the time of the CR flux depressions. The correlation between the space-weather parameters and galactic cosmic ray (GCR) intensity decreases at the two NM stations are statistically significant. The implications of the present space-weather data on cosmic ray (CR) intensity depressions are highlighted.

It is widely accepted that coronal magnetic flux ropes are the core structures of large-scale solar eruptive activities, which inflict dramatic impacts on the solar-terrestrial system. Previous studies have demonstrated that varying magnetic properties of a coronal flux rope system could result in a catastrophe of the rope, which may trigger solar eruptive activities. Since the total mass of a flux rope also plays an important role in stabilizing the rope, we use 2.5-dimensional magnetohydrodynamic (MHD) numerical simulations in this letter to investigate how a flux rope evolves as its total mass varies. It is found that an unloading process that decreases the total mass of the rope could result in an upward (eruptive) catastrophe in the flux rope system, during which the rope jumps upward and the magnetic energy is released. This indicates that mass unloading processes could initiate the eruption of the flux rope. Moreover, when the system is not too diffusive, there is also a downward (confined) catastrophe that could be caused by mass loading processes, via which the total mass accumulates. The magnetic energy, however, is increased during the downward catastrophe, indicating that mass loading processes could cause confined activities that may contribute to the storage of energy before the onset of coronal eruptions.

Context. Computational techniques are essential for mining large databases produced in modern surveys with value-added products. Aims. This paper presents a machine learning procedure to carry out simultaneously galaxy morphological classification and photometric redshift estimates. Currently, only spectral energy distribution (SED) fitting has been used to obtain these results all at once. Methods. We used the ancillary data gathered in the OTELO catalog and designed a non-sequential neural network that accepts optical and near-infrared photometry as input. The network transfers the results of the morphological classification task to the redshift fitting process to ensure consistency between both procedures. Results. The results successfully recover the morphological classification and the redshifts of the test sample, reducing catastrophic redshift outliers produced by SED fitting and avoiding possible discrepancies between independent classification and redshift estimates. Our technique may be adapted to include galaxy images to improve the classification.

Generalised Faraday rotation can induce frequency-dependent conversion between the linear and circular polarisation spectra of compact radio sources such as pulsars, fast radio bursts and active galactic nuclei. I devise a simple phenomenological model that can be used to measure the effects of generalised Faraday rotation on the linearly and circularly polarised spectra of these sources. The model is theory-agnostic, with an arbitrary wavelength dependence, and hence can accommodate for a variety of potential generalised Faraday rotation inducing media. It can also be combined with additional observables to infer the physical properties of the intervening medium.

The quantification of all possible waves and instabilities in a given system is of paramount importance, and knowledge of the full magnetohydrodynamic (MHD) spectrum allows one to predict the (in)stability of a given equilibrium state. This is highly relevant in many (astro)physical disciplines, and when applied to the solar atmosphere it may yield various new insights in processes like prominence formation and coronal loop oscillations. In this work we present a detailed, high-resolution spectroscopic study of the solar atmosphere, where we use our newly developed Legolas code to calculate the full spectrum with corresponding eigenfunctions of equilibrium configurations that are based on fully realistic solar atmospheric models, including gravity, optically thin radiative losses and thermal conduction. Special attention is given to thermal instabilities, known to be responsible for the formation of prominences, together with a new outlook on the thermal and slow continua and how they behave in different chromospheric and coronal regions. We show that thermal instabilities are unavoidable in our solar atmospheric models and that there exist certain regions where both the thermal, slow and fast modes all have unstable wave mode solutions. We also encounter regions where the slow and thermal continua become purely imaginary and merge on the imaginary axis. The spectra discussed in this work illustrate clearly that thermal instabilities (both discrete and continuum modes) and magneto-thermal overstable propagating modes are ubiquitous throughout the solar atmosphere, and may well be responsible for much of the observed fine-structuring and multi-thermal dynamics.

In the present paper we study the possibility of a simultaneous generation of radio waves and soft $X$-rays by means of the quasi-linear diffusion (QLD) in the anomalous pulsar AXP 4U 0142+61. Considering the magnetosphere composed of the so-called beam component and the plasma component respectively, we argue that the frozen-in condition will inevitably lead to the generation of the unstable cyclotron waves. These waves, via the QLD, will in turn influence the particle distribution function, leading to certain values of the pitch angles, thus to an efficient synchrotron mechanism, producing soft $X$-ray photons. We show that for physically reasonable parameters of magnetospheric plasmas, the QLD can provide generation of radio waves in the following interval $40$ MHz-$111$ MHz connected to soft $X$-rays for the domain $0.3$keV-$1.4$keV.

The transition from supernovae (SNe) to supernova remnants (SNRs) remains poorly understood, given the age gap between well-studied examples of the two. In order to bridge this gap, we analysed archival Chandra data for some of the oldest supernovae detected in X-rays, in order to extend their light curves out to late times. We fitted the spectra with thermal models. All the SNe with multiple X-ray data points were found to have similar X-ray luminosity, which was decreasing with time. The X-ray luminosity will likely continue to decrease while the SNe are evolving in a wind medium, but is anticipated to increase in the Sedov phase when the SNe are interacting with a constant density interstellar medium, bringing it in line with observed SNRs.

This paper explores the problem of analytically approximating the orbital state for a subset of orbits in a rotating potential with oblateness and ellipticity perturbations. This is done by isolating approximate differential equations for the orbit radius and other elements. The conservation of the Jacobi integral is used to make the problem solvable to first-order in the perturbations. The solutions are characterized as small deviations from an unperturbed circular orbit. The approximations are developed for near-circular orbits with initial mean motion $n_{0}$ around a body with rotation rate $c$. The approximations are shown to be valid for values of angular rate ratio $\Gamma = c/n_{0} > 1$, with accuracy decreasing as $\Gamma \rightarrow 1$, and singularities at and near $\Gamma = 1$. Extensions of the methodology to eccentric orbits are discussed, with commentary on the challenges of obtaining generally valid solutions for both near-circular and eccentric orbits.

We analyse a sample of twelve galaxy clusters, from the Kapteyn IAC WEAVE INT Cluster Survey (KIWICS) looking for dwarf galaxy candidates. By using photometric data in the $r$ and $g$ bands from the Wide Field Camera (WFC) at the 2.5-m Isaac Newton telescope (INT), we select a sample of bright dwarf galaxies (M$_r$ $\leq$ -15.5 mag) in each cluster and analyse their spatial distribution, stellar colour, and as well as their S\'ersic index and effective radius. We quantify the dwarf fraction inside the $R_{200}$ radius of each cluster, which ranges from $\sim$ 0.7 to $\sim$ 0.9. Additionally, when comparing the fraction in the inner region with the outermost region of the clusters, we find that the fraction of dwarfs tends to increase going to the outer regions. We also study the clustercentric distance distribution of dwarf and giant galaxies (M$_r$ $<$ -19.0 mag), and in half of the clusters of our sample, the dwarfs are distributed in a statistically different way as the giants, with the giant galaxies being closer to the cluster centre. We analyse the stellar colour of the dwarf candidates and quantify the fraction of blue dwarfs inside the $R_{200}$ radius, which is found to be less than $\sim$ 0.4, but increases with distance from the cluster centre. Regarding the structural parameters, the S\'ersic index for the dwarfs we visually classify as early type dwarfs tends to be higher in the inner region of the cluster. These results indicate the role that the cluster environment plays in shaping the observational properties of low-mass halos.

Abstract We investigate a short-period W UMa binary KIC 9026766 with an orbital period of 0.2721278 days in the Kepler field of view. By using an automated q-search for the folded light curve and producing a synthetic light curve for this object based on the PHOEBE code, we calculate the fundamental stellar parameters. We also analyze the O-C curve of the primary minima. The orbital period changes can be attributed to the combination of an upward quadratic function and light-travel time effect due to a possible third body with a minimum mass of 0.029 solar mass and an orbital period of 972.5866 days. The relative luminosity of the primary and secondary eclipses (Min I Min II) is calculated. The periodogram of the residuals of the LTT effect, and Min I Min II show peaks with the same period of 0.8566 days. The background effect of two nearby stars on our target is the possible reason for this signal. By considering the amplitudes and periods of the remaining signals in the OC curve of minima, spot motion is possible.

We report on the continuum and polarization observations of the Cygnus Loop supernova remnant (SNR) conducted by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). FAST observations provide high angular resolution and high sensitivity images of the SNR, which will help to disentangle its nature. We obtained Stokes I, Q and U maps over the frequency range of 1.03 - 1.46 GHz split into channels of 7.63 kHz. The original angular resolution is in the range of ~3 arcmin - ~3.8 arcmin, and we combined all the data at a common resolution of 4 arcmin. The temperature scale of the total intensity and the spectral index from the in-band temperature-temperature plot are consistent with previous observations, which validates the data calibration and map-making procedures. The rms sensitivity for the band-averaged total-intensity map is about 20 mK in brightness temperature, which is at the level of confusion limit. For the first time, we apply rotation measure (RM) synthesis to the Cygnus Loop to obtain the polarization intensity and RM maps. The rms sensitivity for polarization is about 5 mK, far below the total-intensity confusion limit. We also obtained RMs of eight extra-galactic sources, and demonstrate that the wide-band frequency coverage helps to overcome the ambiguity of RM determinations.

Cosmological shock waves are ubiquitous to cosmic structure formation and evolution. As a consequence, they play a major role in the energy distribution and thermalization of the intergalactic medium (IGM). We analyze the Mach number distribution in the Dianoga simulations of galaxy clusters performed with the SPH code GADGET-3. The simulations include the effects of radiative cooling, star formation, metal enrichment, supernova and active galactic nuclei feedback. A grid-based shock-finding algorithm is applied in post-processing to the outputs of the simulations. This procedure allows us to explore in detail the distribution of shocked cells and their strengths as a function of cluster mass, redshift and baryonic physics. We also pay special attention to the connection between shock waves and the cool-core/non-cool core (CC/NCC) state and the global dynamical status of the simulated clusters. In terms of general shock statistics, we obtain a broad agreement with previous works, with weak (low-Mach number) shocks filling most of the volume and processing most of the total thermal energy flux. As a function of cluster mass, we find that massive clusters seem more efficient in thermalising the IGM and tend to show larger external accretion shocks than less massive systems. We do not find any relevant difference between CC and NCC clusters. However, we find a mild dependence of the radial distribution of the shock Mach number on the cluster dynamical state, with disturbed systems showing stronger shocks than regular ones throughout the cluster volume.

Very little work has been done on star formation in dwarf lenticular galaxies (S0s). We present 2D-spectroscopic and millimetre observations made by Centro Astronomico Hispano Aleman (CAHA) 3.5 m optical and the IRAM-30 m millimetre telescopes, respectively, for a sample of four dwarf S0 galaxies with multiple star formation regions in the field environment. We find that although most of the sources deviate from the star forming main sequence relation, they all follow the Kennicutt-Schmidt law. After comparing the stellar and Halpha kinematics, we find that the velocity fields of both stars and ionized gas do not show regular motion and the velocity dispersions of stars and ionized gas are low in the regions with high star formation, suggesting these star-forming S0 galaxies still have significant rotation. This view can be supported by the result that most of these dwarf S0 galaxies are classified as fast rotators. The ratio of average atomic gas mass to stellar mass (~ 47%) is much greater than that of molecular gas mass to stellar mass (~ 1%). In addition, the gas-phase metallicities in the star-forming regions are lower than that of the non-star-forming regions. These results indicate that the extended star formation may originate from the combination of abundant atomic hydrogen, long dynamic time scale and low-density environment.

Two groups of astronomers used large telescopes Keck and VLT for decades to observe trajectories of bright stars near the Galactic Centre. Based on results of their observations astronomers concluded that trajectories of the stars are roughly elliptical and foci of the orbits are approximately coincide with the Galactic Centre position. It gives an opportunity to claim that the Newtonian potential of point like mass around $4.3\times 10^6 M_\odot$ is a good initial approximation for the gravitational potential near the Galactic Centre. In the last years, the astronomers found that gravitational redshift of S2 star near pericenter passage in May 2018 is in accordance with general relativity predictions. In 2020 the GRAVITY team found that the observed relativistic precession of S2 star orbit is also consistent with theoretical estimates calculated for a weak gravitational field approximation in a Schwarzschild black hole. In last years a a self-gravitating dark matter core--halo distribution suggested by Ruffini, Arguelles and Rueda (MNRAS, 2015) (RAR model) was proposed and recently Becerra-Vergara et al. (MNRAS, 2021) claimed that this model provides a better fit of trajectories of bright stars in comparison with the conventional model with the supermassive black hole. We confirm that in the case of this dark matter distribution model for a dense core trajectories of test bodies are elliptical but in this case centers (not foci) of these ellipses should coincide with the Galactic Centre and orbital periods do not depend on semi-major axis and it contradicts observational data and therefore, we concluded supermassive black hole is a preferable model in comparison with the a dense core--diluted halo density profile for the Galactic Centre.

In this manuscript, an interesting blue Active Galactic Nuclei (AGN) SDSS J154751.94+025550 (=SDSS J1547) is reported with very different line profiles of broad Balmer emission lines: double-peaked broad H$\beta$ but single-peaked broad H$\alpha$. SDSS J1547 is the first AGN with detailed discussions on very different line profiles of the broad Balmer emission lines, besides the simply mentioned different broad lines in the candidate for a binary black hole (BBH) system in SDSS J0159+0105. The very different line profiles of the broad Balmer emission lines can be well explained by different physical conditions to two central BLRs in a central BBH system in SDSS J1547. Furthermore, the long-term light curve from CSS can be well described by a sinusoidal function with a periodicity about 2159days, providing further evidence to support the expected central BBH system in SDSS J1547. Therefore, it is interesting to treat different line profiles of broad Balmer emission lines as intrinsic indicators of central BBH systems in broad line AGN. Under assumptions of BBH systems, 0.125\% of broad line AGN can be expected to have very different line profiles of broad Balmer emission lines. Future study on more broad line AGN with very different line profiles of broad Balmer emission lines could provide further clues on the different line profiles of broad Balmer emission lines as indicator of BBH systems.

The VERITAS Imaging Air Cherenkov Telescope (IACT) array was augmented in 2019 with high-speed focal plane electronics to allow its use for Stellar Intensity Interferometry (SII) observations. Since January 2019, the VERITAS Stellar Interferometer (VSII) recorded more than 250 hours of moonlit observations on 39 different bright stars and binary systems ($m_V < 3.74$) at an effective optical wavelength of 416 nm. These observations resulted in the measurement of the diameters of several stars with better than 5% resolution. This talk will describe the status of the VSII survey and analysis.

A recent measurement of the proper motion of PSR J0908-4913 shows that it is a fast moving object at a distance of some 3 kpc. Here we present evidence that the pulsar is located at the edge of a previously unknown, filled-centre supernova remnant, G270.4-1.0. The velocity vector of the pulsar points directly away from the centre of the remnant. The putative association of the pulsar with SNR G270.4-1.0 implies the pulsar is ~12kyr old, significantly less than its characteristic age of 110kyr. We show that the rotation axis of the pulsar points in the direction of the proper motion. Rotation measure and dispersion measure variations are seen over time, likely indicating the pulsar is passing behind a filament of the remnant.

We present optical photometry and spectroscopy of the superluminous SN 2002gh from maximum light to $+202$ days, obtained as part of the Carnegie Type II Supernova (CATS) project. SN 2002gh is among the most luminous discovered supernovae ever, yet it remained unnoticed for nearly two decades. Using Dark Energy Camera archival images we identify the potential SN host galaxy as a faint dwarf galaxy, presumably having low metallicity, and in an apparent merging process with other nearby dwarf galaxies. We show that SN 2002gh is among the brightest hydrogen-poor SLSNe with $M_{V} = -22.40 \pm 0.02$, with an estimated peak bolometric luminosity of $2.6 \pm 0.2 \times 10^{44}$ erg s$^{-1}$. We discount the decay of radioactive nickel as the main SN power mechanism, and assuming that the SN is powered by the spin down of a magnetar we obtain two alternative solutions. The first case, is characterized by significant magnetar power leakage, and $M_{\mathrm{ej}}$ between 0.8 and 1.6 $M_{\odot}$, $P_{\mathrm{spin}} = 3.4$ ms, and $B = 5 \times 10^{13}$ G. The second case does not require power leakage, resulting in a huge ejecta mass of about 30 $M_{\odot}$, a fast spin period of $P_{\mathrm{spin}} \sim 1$ ms, and $B\sim 1.6 \times 10^{14}$ G. We estimate a zero-age main-sequence mass between 16 and 19 $M_{\odot}$ for the first case and of about 135 $M_{\odot}$ for the second case. The latter case would place the SN progenitor among the most massive stars observed to explode as a SN.

Since the collisional mean free path of charged particles in hot accretion flows can be significantly larger than the typical length-scale of the accretion flows, the gas pressure is anisotropic to magnetic field lines. For such a large collisional mean free path, the resistive dissipation can also play a key role in hot accretion flows. In this paper, we study the dynamics of resistive hot accretion flows in the presence of anisotropic pressure. We present a set of self-similar solutions where the flow variables are assumed to be a function only of radius. Our results show that the radial and rotational velocities and the sound speed increase considerably with the strength of anisotropic pressure. The increase in infall velocity and in sound speed are more significant if the resistive dissipation is taken into account. We find that such changes depend on the field strength. Our results indicate that the resistive heating is $10\%$ of the heating by the work done by anisotropic pressure when the strength of anisotropic pressure is 0.1. This value becomes higher when the strength of anisotropic pressure reduces. The increase in disk temperature can lead to heating and acceleration of the electrons in such flows. It helps us to explain the origin of phenomena such as the flares in Galactic Center Sgr A*.

The existence of a giant planet beyond Neptune -- referred to as Planet Nine (P9) -- has been inferred from the clustering of longitude of perihelion and pole position of distant eccentric Kuiper belt objects (KBOs). After updating calculations of observational biases, we find that the clustering remains significant at the 99.6\% confidence level. We thus use these observations to determine orbital elements of P9. A suite of numerical simulations shows that the orbital distribution of the distant KBOs is strongly influenced by the mass and orbital elements of P9 and thus can be used to infer these parameters. Combining the biases with these numerical simulations, we calculate likelihood values for discrete set of P9 parameters, which we then use as input into a Gaussian Process emulator that allows a likelihood computation for arbitrary values of all parameters. We use this emulator in a Markov Chain Monte Carlo analysis to estimate parameters of P9. We find a P9 mass of $6.2^{+2.2}_{-1.3}$ Earth masses, semimajor axis of $380^{+140}_{-80}$ AU, inclination of $16\pm5^\circ$ and perihelion of $300^{+85}_{-60}$ AU. Using samples of the orbital elements and estimates of the radius and albedo of such a planet, we calculate the probability distribution function of the on-sky position of Planet Nine and of its brightness. For many reasonable assumptions, Planet Nine is closer and brighter than initially expected, though the probability distribution includes a long tail to larger distances, and uncertainties in the radius and albedo of Planet Nine could yield fainter objects.

Fast Radio Burst (FRB) dispersion measures (DMs) record the presence of ionized baryons that are otherwise invisible to other techniques enabling resolution of the matter distribution in the cosmic web. In this work, we aim to estimate the contribution to FRB 180924 DM from foreground galactic halos. Localized by ASKAP to a massive galaxy, this sightline is notable for an estimated cosmic web contribution to the DM ($\rm DM_{cosmic} = 220~pc~cm^{-3}$), which is less than the average value at the host redshift ($\rm z = 0.3216$) estimated from the Macquart relation ($280~\rm pc~cm^{-3}$). In the favored models of the cosmic web, this suggests few intersections with foreground halos at small impact parameters ($\lesssim 100$ kpc). To test this hypothesis, we carried out spectroscopic observations of the field galaxies within $\sim$1' of the sightline with VLT/MUSE and Keck/LRIS. Furthermore, we developed a probabilistic methodology that leverages photometric redshifts derived from wide-field DES and WISE imaging. We conclude that there is no galactic halo that closely intersects the sightline and also that the net DM contribution from halos, $\rm DM_{halos}< 45~pc~cm^{-3}$ (95 % c.l.). This value is lower than the $\rm DM_{halos}$ estimated from an "average" sightline ($121~\rm pc~cm^{-3}$) using the Planck $\Lambda CDM$ model and the Aemulus halo mass function and reasonably explains its low $\rm DM_{cosmic}$ value. We conclude that FRB 180924 represents the predicted majority of sightlines in the universe with no proximate foreground galactic halos. Our framework lays the foundation for a comprehensive analysis of FRB fields in the near future.

We put together the experimental results on muon component of extensive air showers (EAS) which were gained with various techniques at the detector complex of the Tien Shan mountain station. According to this comparison, the problem of the EAS muon content in the range of primary cosmic ray energies (1-100)PeV seems to be more complicated than it was usually supposed. Generally, from the models of nuclear interaction it follows that the EAS which have produced gamma-hadron families in the Tien Shan X-ray emulsion chamber did preferably originate from interaction of the light cosmic ray nuclei, such that their muon abundance must be ~1.5 times below an average calculated over all showers. In contrary, the experimental muon counts in the EAS with families demonstrate a (1.5-2)-fold excess above the average, and this difference starts to be observable in the showers with the energy above the 3PeV knee of the primary cosmic ray spectrum. Later on, the rise of muon production in EAS after the knee was confirmed at Tien Shan by another experiment on detection of the neutrons stemmed from interaction of cosmic ray muons. Thus, the results obtained by the two completely different methods do mutually agree with each other but contradict to the common models of hadron interaction.

The abundance of deuterium in giant planet atmospheres provides constraints on the reservoirs of ices incorporated into these worlds during their formation and evolution. Motivated by discrepancies in the measured deuterium-hydrogen ratio (D/H) on Jupiter and Saturn, we present a new measurement of the D/H ratio in methane for Saturn from ground-based measurements. We analysed a spectral cube (covering 1151-1160 cm$^{-1}$ from 6 February 2013) from the Texas Echelon Cross Echelle Spectrograph (TEXES) on NASA's Infrared Telescope Facility (IRTF) where emission lines from both methane and deuterated methane are well resolved. Our estimate of the D/H ratio in stratospheric methane, $1.65_{-0.21}^{+0.27} \times 10^{-5}$ is in agreement with results derived from Cassini CIRS and ISO/SWS observations, confirming the unexpectedly low CH$_{3}$D abundance. Assuming a fractionation factor of $1.34 \pm 0.19$ we derive a hydrogen D/H of $1.23_{-0.23}^{+0.27} \times 10^{-5}$. This value remains lower than previous tropospheric hydrogen D/H measurements of (i) Saturn $2.10 (\pm 0.13) \times 10^{-5}$, (ii) Jupiter $2.6 (\pm 0.7) \times 10^{-5}$ and (iii) the proto-solar hydrogen D/H of $2.1 (\pm 0.5) \times 10^{-5}$, suggesting that the fractionation factor may not be appropriate for stratospheric methane, or that the D/H ratio in Saturn's stratosphere is not representative of the bulk of the planet.

Dust clouds are ubiquitous in the atmospheres of hot Jupiters and affect their observable properties. The alignment of dust grains in the clouds and resulting dust polarization is a promising method to study magnetic fields of exoplanets. Moreover, the grain size distribution plays an important role in physical and chemical processes in the atmospheres, which is rather uncertain in atmospheres. In this paper, we first study grain alignment of dust grains in the atmospheres of hot Jupiters by RAdiative Torques (RATs). We find that silicate grains can be aligned by RATs with the magnetic fields (B-RAT) due to strong magnetic fields of hot Jupiters, but carbonaceous grains of diamagnetic material tend to be aligned with the radiation direction (k-RAT). At a low altitude of $r<2R_{\rm p}$ with $R_{\rm p}$ being the planet radius, only large grains can be aligned, but tiny grains of $a\sim 0.01\mu$m can be aligned at a high altitude of $r>3R_{\rm p}$. We then study rotational disruption of dust grains by the RAdiative Torque Disruption (RATD) mechanism. We find that large grains can be disrupted by RATD into smaller sizes. Grains of high tensile strength are disrupted at an altitude of $r>3R_{\rm p}$, but weak grains can be disrupted at a lower altitude. We suggest that the disruption of large grains into smaller ones can facilitate dust clouds to escape to high altitudes due to lower gravity and may explain the presence of high-altitude clouds in hot Jupiter as well as super-puff atmospheres.

As the precision frontier of astrophysics advances towards the one millimagnitude level, flux calibration of photometric instrumentation remains an ongoing challenge. We present the results of a lab-bench assessment of the viability of monocrystalline silicon solar cells to serve as large-aperture (up to 125mm diameter), high-precision photodetectors. We measure the electrical properties, spatial response uniformity, quantum efficiency (QE), and frequency response of 3$^{rd}$ generation C60 solar cells, manufactured by Sunpower. Our new results, combined with our previous study of these cells' linearity, dark current, and noise characteristics, suggest that these devices hold considerable promise, with QE and linearity that rival those of traditional, small-aperture photodiodes. We argue that any photocalibration project that relies on precise knowledge of the intensity of a large-diameter optical beam should consider using solar cells as calibrating photodetectors.

We analyze the photometric variability of over 6,000 active galactic nuclei (AGNs) from the Sloan Digital Sky Survey Stripe 82. We recover the spectral energy distribution (SED) of the variable flux as a function of wavelength. For rest wavelengths longer than $\sim 2200$ Angstroms we find that the SED of the variable component of the bluest AGNs is in agreement with the $F_{\nu} \propto \nu^{+1/3}$ spectrum predicted for an externally-illuminated accretion disc. We confirm there is some residual variable emission corresponding to the "small blue bump" and other broad-line region variability. We interpret steeper optical spectra of the variable component as being due to intrinsic reddening. This is supported by the correlation of the Balmer decrement with the colour excess of the variable component. We find that the median reddening of the SDSS AGNs in Stripe 82 is $E(B-V) \thickapprox 0.17$ in agreement with the reddening derived from the Balmer decrement.

The Fermi Bubbles were discovered about a decade ago in the {\it Fermi}-LAT data as a double-lobe structure extending up to 55 deg. in Galactic latitudes above and below the Galactic Center. At the moment their origin is still unknown. The H.E.S.S. collaboration is currently performing the first ever survey in TeV gamma rays of the Milky Way inner region: the Inner Galaxy Survey. This survey is intended to achieve the best sensitivity to faint and diffuse emissions in a region of several degrees around the Galactic Centre. It provides an unprecedented sensitivity to dark matter signals, new diffuse emissions, and TeV outflows from the Galactic Centre. Understanding the properties of the Fermi Bubbles at low Galactic latitudes will provide key insights into their origin. We search for TeV emission at the base of the Fermi Bubbles using low-latitude spatial templates. The first results obtained with the 2014-2020 H.E.S.S. observations will be reported.

We report the discovery of CXOU J191100-595621 and CXOU J191012-595619, two galaxy clusters serendipitously detected in the direction of globular cluster NGC 6752, based on archival {\it Chandra} observations with a total exposure time of $\sim 344$ ks. The deep {\it Chandra} X-ray data enabled us to measure properties of both systems, which result in a redshift of $z=0.239\pm0.013$ and $z=0.375\pm0.016$, a temperature of $kT=3.32^{+0.57}_{-0.46}$ keV and $kT=3.71^{+1.18}_{-0.86}$ keV, an iron abundance of $Z_{\rm Fe}=0.64^{+0.34}_{-0.29}Z_{\rm Fe\odot}$ and $Z_{\rm Fe}=1.29^{+0.97}_{-0.65}Z_{\rm Fe\odot}$, and a rest-frame full band (0.5-7 keV) luminosity of $L_{\rm X}=9.2^{+1.2}_{-1.1}\times 10^{43} {\rm \, erg\, s^{-1}}$ and $L_{\rm X}=9.9^{+2.7}_{-2.2}\times 10^{43} {\rm \, erg\, s^{-1}}$ for CXOU J191100-595621 and CXOU J191012-595619, respectively. The temperature profile of CXOU J191100-595621 is found to decreases with decreasing radius, indicating a cool core in this cluster. The hydrostatic equilibrium estimation suggests the clusters are moderately weighted, with $M_{500}=(1.3\pm0.4)\times 10^{14}\, M_{\odot}$ and $M_{500}=(2.0\pm1.5)\times 10^{14}\, M_{\odot}$, respectively. We search for optical and radio counterparts of X-ray point sources in the clusters. Three active galactic nuclei are found, among which one is identified with a narrow-angle-tail radio galaxy, and one is found to associated with the brightest central galaxy (BCG) of CXOU J191100-595621.

We report the highest spatial resolution measurement of magnetic fields in M17 using thermal dust polarization taken by SOFIA/HAWC+ centered at 154 $\mu$m wavelength. Using the Davis-Chandrasekhar-Fermi method, we found the presence of strong magnetic fields of $980 \pm 230\;\mu$G and $1665 \pm 885\;\mu$G in lower-density (M17-N) and higher-density (M17-S) regions, respectively. The magnetic field morphology in M17-N possibly mimics the fields in gravitational collapse molecular cores while in M17-S the fields run perpendicular to the matter structure and display a pillar and an asymmetric hourglass shape. The mean values of the magnetic field strength are used to determine the Alfv\'enic Mach numbers ($\mathcal{M_A}$) of M17-N and M17-S which turn out to be sub-Alfv\'enic, or magnetic fields dominate turbulence. We calculate the mass-to-flux ratio, $\lambda$, and obtain $\lambda=0.07$ for M17-N and $0.28$ for M17-S. The sub-critical values of $\lambda$ are in agreement with the lack of massive stars formed in M17. To study dust physics, we analyze the relationship between the dust polarization fraction, $p$, and the thermal emission intensity, $I$, gas column density, $N({\rm H_2})$, and dust temperature, $T_{\rm d}$. The polarization fraction decreases with intensity as $I^{-\alpha}$ with $\alpha = 0.51$. The polarization fraction also decreases with increasing $N(\rm H_{2})$, which can be explained by the decrease of grain alignment by radiative torques (RATs) toward denser regions with a weaker radiation field and/or tangling of magnetic fields. The polarization fraction tends to increase with $T_{\rm d}$ first and then decreases when $T_ {\rm d} > 50$ K. The latter feature seen in the M17-N, where the gas density changes slowly with $T_{d}$, is consistent with the RAT disruption effect.

The polarisation Sagnac speedmeter interferometer has the potential to replace the Michelson interferometer as the instrumental basis for future generations of ground-based gravitational wave detectors. The quantum noise benefit of this speedmeter is dependent on high-quality polarisation optics, the polarisation beam-splitter (PBS) and quarter-waveplate (QWP) optics that are key to this detector configuration and careful consideration of the effect of birefringence in the arm cavities of the interferometer. A PBS with an extinction ratio of better than 4000 in transmission and 700 in reflection for a $41^{\circ}$ angle of incidence was characterised along with a QWP of birefringence of $\frac{\lambda}{4} + \frac{\lambda}{324}$. The cavity mirror optics of a 10m prototype polarisation Sagnac speedmeter were measured to have birefringence in the range $1\times10^{-3}$ to $2\times10^{-5}$ radians. This level of birefringence, along with the QWP imperfections, can be canceled out by careful adjustment of the QWP angle, to the extent that the extinction ratio of the PBS is the leading limitation for the polarisation Sagnac speedmeter in terms of polarisation effects.

Weak lensing by large-scale structure is a powerful probe of cosmology if the apparent alignments in the shapes of distant galaxies can be accurately measured. We study the performance of a fully data-driven approach, based on MetaDetection, focusing on the more realistic case of observations with an anisotropic PSF. Under the assumption that PSF anisotropy is the only source of additive shear bias, we show how unbiased shear estimates can be obtained from the observed data alone. To do so, we exploit the finding that the multiplicative shear bias obtained with MetaDetection is nearly insensitive to the PSF ellipticity. In practice, this assumption can be validated by comparing the empirical corrections obtained from observations to those from simulated data. We show that our data-driven approach meets the stringent requirements for upcoming space and ground-based surveys, although further optimisation is possible.

Traditionally, 1D models based on scaling laws have been used to parameterized convective heat transfer rocks in the interior of terrestrial planets like Earth, Mars, Mercury and Venus to tackle the computational bottleneck of high-fidelity forward runs in 2D or 3D. However, these are limited in the amount of physics they can model (e.g. depth dependent material properties) and predict only mean quantities such as the mean mantle temperature. We recently showed that feedforward neural networks (FNN) trained using a large number of 2D simulations can overcome this limitation and reliably predict the evolution of entire 1D laterally-averaged temperature profile in time for complex models [Agarwal et al. 2020]. We now extend that approach to predict the full 2D temperature field, which contains more information in the form of convection structures such as hot plumes and cold downwellings. Using a dataset of 10,525 two-dimensional simulations of the thermal evolution of the mantle of a Mars-like planet, we show that deep learning techniques can produce reliable parameterized surrogates (i.e. surrogates that predict state variables such as temperature based only on parameters) of the underlying partial differential equations. We first use convolutional autoencoders to compress the temperature fields by a factor of 142 and then use FNN and long-short term memory networks (LSTM) to predict the compressed fields. On average, the FNN predictions are 99.30% and the LSTM predictions are 99.22% accurate with respect to unseen simulations. Proper orthogonal decomposition (POD) of the LSTM and FNN predictions shows that despite a lower mean absolute relative accuracy, LSTMs capture the flow dynamics better than FNNs. When summed, the POD coefficients from FNN predictions and from LSTM predictions amount to 96.51% and 97.66% relative to the coefficients of the original simulations, respectively.

The residuals of the power spectra of WMAP and Planck's cosmic microwave background (CMB) anisotropies data are known to exhibit a few interesting anomalies at different scales with marginal statistical significance. Combining bottom-up and top-down model-building approaches and using a pipeline that efficiently compares model predictions with data, we construct a model of primordial standard clock that is able to link and address the anomalies at both the large and small scales. This model, and its variant, provide some of the best fits to the feature anomalies in CMB. According to Bayes evidences, these models are currently statistically indistinguishable from the Standard Model. We show that the difference between them will soon become statistically significant with various higher quality data on the CMB polarization. We demonstrate that such a model-building and data-analyses process may be used to uncover a portion of detailed evolutionary history of our universe during its primordial epoch.

Silicon photomultipliers (SiPMs) have become the baseline choice for cameras of the small-sized telescopes (SSTs) of the Cherenkov Telescope Array (CTA). On the other hand, SiPMs are relatively new to the field and covering large surfaces and operating at high data rates still are challenges to outperform photomultipliers (PMTs). The higher sensitivity in the near infra-red and longer signals compared to PMTs result in higher night sky background rate for SiPMs. However, the robustness of the SiPMs represents a unique opportunity to ensure long-term operation with low maintenance and better duty cycle than PMTs. The proposed camera for large size telescopes will feature $0.05 degree pixels, low power and fast front-end electronics and a fully digital readout. In this work, we present the status of dedicated simulations and data analysis for the performance estimation. The design features and the different strategies identified, so far, to tackle the demanding requirements and the improved performance are described.

The world is changing fast, and so is the space sector. Planning for large scientific experiments two decades ahead may no longer be the most sensible approach. I develop the argument that large science experiments are becoming comparable to terrestrial civil infrastructures in terms of cost. As a result, these should incorporate plans for a return on investment (or impact, not necessarily economic), require a different approach for inter-division coordination within the European Space Agency(ESA), and a broader participation of all society stakeholders (civil society representatives, and the broader public). Defining which experiments will be relevant two decades ahead adds rigidity and quenches creativity to the development of cutting edge science and technology. This is likely to discourage both senior and earlier career professionals into supporting such long-term (and often precarious) plans. A more sensible strategy would be increasing the rate of smaller well understood experiments, engage more society sectors in the development of a truly space-bound infrastructure, and formulate a strategy more in tune with the challenges faced by our society and planet. We argue that such strategy would lead to equally large -- even larger -- scale experiments in the same time-scale, while providing economic returns and a common sense of purpose. A basic but aggressive road map is outlined.

We use the Sloan-Digital Sky Survey quasar catalog to statistically infer the local abundance of black holes heavier than $10^8M_\odot$, which allows us to estimate the detection prospect of super-massive black holes by future observational campaigns. We find that the upcoming James Webb Space Telescope (JWST) and the Extremely Large Telescope (ELT) should be able to resolve, with integral field spectroscopy techniques, the gravitational influence of $\sim10^3$ black holes within a sphere of $\sim50\;\rm Mpc$. A Very-Long Baseline (VLB) observatory with one receiver placed in a geostationary orbit, is predicted to capture $\sim10$ images of the silhouette of a black hole, similar to the image of $\rm M87^*$ recently performed by the Event Horizon Telescope.

The Cherenkov Telescope Array (CTA) will be the next-generation observatory in the field of very-high-energy (20 GeV to 300 TeV) gamma-ray astroparticle physics. Classically, data analysis in the field maximizes sensitivity by applying quality cuts on the data acquired. These cuts, optimized using Monte Carlo simulations, select higher quality events from the initial dataset. Subsequent steps of the analysis typically use the surviving events to calculate one set of instrument response functions (IRFs). An alternative approach is the use of event types, as implemented in experiments such as the Fermi-LAT. In this approach, events are divided into sub-samples based on their reconstruction quality, and a set of IRFs is calculated for each sub-sample. The sub-samples are then combined in a joint analysis, treating them as independent observations. This leads to an improvement in performance parameters such as sensitivity, angular and energy resolution. Data loss is reduced since lower quality events are included in the analysis as well, rather than discarded. In this study, machine learning methods will be used to classify events according to their expected angular reconstruction quality. We will report the impact on CTA high-level performance when applying such an event-type classification, compared to the classical procedure.

The physical sizes of TDE accretion discs are regularly inferred, from the modelling of the TDEs X-ray spectrum as a single temperature blackbody, to be smaller than the plausible event horizons of the black holes which they occur around - a clearly unphysical result. In this Letter we demonstrate that the use of single-temperature blackbody functions results in the systematic underestimation of TDE accretion disc sizes by as much as an order-of-magnitude. In fact, the radial `size' inferred from fitting a single temperature blackbody to an observed accretion disc X-ray spectrum does not even positively correlate with the physical size of that accretion disc. We further demonstrate that the disc-observer inclination angle and absorption of X-ray photons may both lead to additional underestimation of the radial sizes of TDE discs, but by smaller factors. To rectify these issues we present a new fitting function which accurately reproduces the size of an accretion disc from its 0.3-10 keV X-ray spectrum. Unlike traditional approaches, this new fitting function does not assume that the accretion disc has reached a steady state configuration, an assumption which is unlikely to be satisfied by most TDEs. An XSPEC implementation of this new fitting function is available at github.com/andymummeryastro/TDEdiscXraySpectrum.

We examine archival XMM-Newton data on the extremely variable narrow-line Seyfert 1 (NLS1) active galactic nucleus (AGN) 1H 0707-495. We construct fractional excess variance (Fvar) spectra for each epoch, including the recent 2019 observation taken simultaneously with eROSITA. We explore both intrinsic and environmental absorption origins for the variability in different epochs, and examine the effect of the photoionised emission lines from outflowing gas. In particular, we show that the unusual soft variability first detected by eROSITA in 2019 is due to a combination of an obscuration event and strong suppression of the variance at 1 keV by photoionised emission, which makes the variance below 1 keV appear more extreme. We also examine the variability on long timescales, between observations, and find that it is well described by a combination of intrinsic variability and absorption variability. We suggest that the typical extreme high frequency variability which 1H 0707-495 is known for is intrinsic to the source, but the large amplitude, low frequency variability that causes prolonged low-flux intervals is likely dominated by variable low-ionisation, low velocity absorption.

We present near-infrared K-band spectra for a sample of 7 Class 0 protostars in the Perseus and Orion star-forming regions. We detect Br gamma, CO overtone, and H2 emission, features that probe the near circumstellar environment of the protostar and reveal evidence of magnetospheric accretion, a hot inner disk atmosphere, and outflows, respectively. Comparing the properties of these features with those of Class I sources from the literature, we find that their Br gamma and CO emission are generally consistent in strength and velocity width. The Br gamma line profiles are broad and centrally peaked, with FWHMs of 200 km/s and wings extending to 300 km/s. The line ratios of our H2 emission features, which are spatially extended for some sources, are consistent with shock excitation and indicate the presence of strong jets or a disk wind. Within our small sample, the frequency of CO band emission (67%) is high relative to that of Class I samples (15%), indicating that Class 0s have high inner disk accretion rates, similar to those of the most actively accreting Class I sources. Collectively, our results suggest that Class 0 sources have similar accretion mechanisms to the more evolved classes, with strong organized stellar magnetic fields established at the earliest observable stage of evolution.

Anisotropies of the cosmic optical background (COB) and cosmic near-IR background (CNIRB) are capable of addressing some of the key questions in cosmology and astrophysics. In this work, we measure and analyze the angular power spectra of the simulated COB and CNIRB in the ultra-deep field of the China Space Station Telescope (CSST-UDF). The CSST-UDF covers about 9 square degrees, with magnitude limits ~28.3, 28.2, 27.6, 26.7 for point sources with 5-sigma detection in the r (0.620 um), i (0.760 um), z (0.915 um), and y (0.965 um) bands, respectively. According to the design parameters and scanning pattern of the CSST, we generate mock data, merge images and mask the bright sources in the four bands. We obtain four angular power spectra from l=200 to 2,000,000 (from arcsecond to degree), and fit them with a multi-component model including intrahalo light (IHL) using the Markov chain Monte Carlo (MCMC) method. We find that the signal-to-noise ratio (SNR) of the IHL is larger than 8 over the range of angular scales that are useful for astrophysical studies (l~10,000-400,000). Comparing to previous works, the constraints on the model parameters are improved by factors of 3~4 in this study, which indicates that the CSST-UDF survey can be a powerful probe on the cosmic optical and near-IR backgrounds.

The identification of the first confirmed neutron star - black hole (NS-BH) binary mergers by the LIGO, Virgo and KAGRA collaboration provides the opportunity to investigate the observed properties of the early sample of confirmed and candidate events. Here, we focus primarily on the tilt angle of the black hole's spin relative to the orbital angular momentum vector of the binary, and the implications for understanding the physical processes that determine this tilt. While the tilt angle of GW200105 is unconstrained, the posterior tilt distributions of both GW200115 and the candidate event GW190426_152155 peak at significantly anti-aligned orientations (though both display wide distributions). If taken close to the peaks of their posteriors, both would be significant outliers from the posterior predictive tilt distribution of binary black hole mergers. Producing these tilts would require stronger natal kicks than are typically considered (and preferentially-polar kicks would be ruled out), and/or an additional source of tilt such as stable mass transfer. The early sample of NS-BH events are less massive than expected for classical formation channels, and may provide evidence for efficient mass transfer that results in the merger of more massive NS-BH binaries before their evolution to the compact phase is complete. We predict that future gravitational-wave detections of NS-BH events will continue to display total binary masses of $\approx 7$M$_{\odot}$ and mass ratios of $q \sim 3$ if this is the correct interpretation. Large tilts in a significant fraction of merging NS-BH systems would weaken the prospects for electromagnetic detection. However, we show that EM observations, including non-detections, can significantly tighten the constraints on spin and mass ratio.

It was first observed in the 1970s that the dwarf galaxies surrounding our Milky Way, so-called satellites, appear to be arranged in a thin, vast plane. Similar discoveries have been made around additional galaxies in the local Universe such as Andromeda, Centaurus A, and potentially M83. In the specific cases with available kinematic data, the dwarf satellites also appear to preferentially co-orbit their massive host galaxy. Planes of satellites are rare in the lambda cold dark matter ($\Lambda$CDM) paradigm, although they may be a natural consequence of projection effects. Such a phase-space correlation, however, remains difficult to explain. In this work we analyzed the 2D spatial distribution of 2210 dwarf galaxies around early-type galaxies (ETGs) in the low-to-medium density fields of the "Mass Assembly of early-Type GaLAxies with their fine Structures" (MATLAS) survey. Under the assumption that the dwarfs are satellite members of the central massive ETG, we identified flattened structures using both a variation in the Hough transform and total least square (TLS) fitting. In 119 satellite systems, we find 31 statistically significant flattened dwarf structures using a combination of both methods with subsequent Monte Carlo (MC) simulations with random data. The vast majority of these dwarf structures lie within the estimated virial radii of the massive host. The major axes of these systems are aligned better than 30{\deg} with the estimated orientation of the large-scale structure in nine (50%) cases. Additional distance measurements and future kinematic studies will be required to confirm the planar nature of these structures and to determine if they are corotating systems.

Context. We present precise radial-velocity measurements of five solar-type stars observed with the HARPS Echelle spectrograph mounted on the 3.6-m telescope in La Silla (ESO, Chile). With a time span of more than 10 years and a fairly dense sampling, the survey is sensitive to low mass planets down to super-Earths on orbital periods up to 100 days. Aims. Our goal was to search for planetary companions around the stars HD39194, HD93385, HD96700, HD154088, and HD189567 and use Bayesian model comparison to make an informed choice on the number of planets present in the systems based on the radial velocity observations. These findings will contribute to the pool of known exoplanets and better constrain their orbital parameters. Methods. A first analysis was performed using the DACE (Data & Analysis Center for Exoplanets) online tools to assess the activity level of the star and the potential planetary content of each system. We then used Bayesian model comparison on all targets to get a robust estimate of the number of planets per star. We did this using the nested sampling algorithm PolyChord. For some targets, we also compared different noise models to disentangle planetary signatures from stellar activity. Lastly, we ran an efficient MCMC (Markov chain Monte Carlo) algorithm for each target to get reliable estimates for the planets' orbital parameters. Results. We identify 12 planets within several multiplanet systems. These planets are all in the super-Earth and sub-Neptune mass regime with minimum masses ranging between 4 and 13 M$_\oplus$ and orbital periods between 5 and 103 days. Three of these planets are new, namely HD 93385 b, HD 96700 c, and HD 189567 c.

Dark spots on the surface of active stars produce changes in the shapes of the spectral lines that mimic spurious Doppler shifts, compromising the detection of small planets by means of the radial velocity (RV) technique. Modelling the spot-driven RV variability (known as ``jitter'') and how it affects the RV data sets is therefore crucial to design efficient activity-filtering techniques and inform observing strategies. Here, we characterise starspots and simulate the radial velocity curves induced by them to determine typical jitter amplitudes for a representative sample of 15 known host stars spanning between F and M spectral type. We collect information on the $\log R'_{\mathrm{HK}}$ activity index from the literature for 205 stars and, due to a lack of data in the temperature range 4000-4500 K, we measure it for ten stars using archival data. Additional stellar parameters required for the simulations are collected from the literature or constrained by observational data, in order to derive realistic estimates. Our results can be used as reference to determine typical peak-to-peak spot-induced RV jitter in the visible domain that can be expected when targeting host stars with different properties.

The reionization process is expected to be prolonged by the small-scale absorbers (SSAs) of ionizing photons, which have been seen as Lyman-limit systems in quasar absorption line observations. We use a set of semi-numerical simulations to investigate the effects of absorption systems on the reionization process, especially their impacts on the neutral islands during the late epoch of reionization (EoR). Three model are studied, i.e. the extreme case of no-SSA model with a high level of ionizing background, the moderate-SSA model with a relatively high level of ionizing background, and the dense-SSA model with a low level of ionizing background. We find that while the characteristic scale of neutral regions decreases during the early and middle stages of reionization, it stays nearly unchanged at about 10 comoving Mpc during the late stage for the no-SSA and moderate-SSA models. However, in the case of weak ionizing background in the dense-SSA model, the characteristic island scale shows obvious evolution, as large islands break into many small ones that are slowly ionized. The evolutionary behavior of neutral islands during the late EoR thus provides a novel way to constrain the abundance of SSAs. We discuss the 21-cm observation with the upcoming Square Kilometre Array (SKA). The different models can be distinguished by the 21-cm power spectrum measurement, and it is also possible to extract the characteristic island scale from the imaging observation with a proper choice of the 21-cm brightness threshold.

In Huybrighs et al., 2020 we investigated energetic proton depletions along Galileo's Europa flyby E26. Based on a particle tracing analysis we proposed that depletions are caused by perturbed electrogmagnetic fields combined with atmospheric charge exchange and possible plumes. One depletion feature identified as a plume signature was shown to be an artefact Jia et al., 2021. Despite that, here we emphasize that Huybrighs et al., 2020 demonstrates that plumes can cause proton depletions and that these features should be sought after. Furthermore, the conclusions on the importance of perturbed electromagnetic fields and atmospheric charge exchange on the depletions are unaffected. We suggest that the artefact's cause is a mistagging of protons as heavier ions by EPD. The artefact prevents us from confirming or excluding that there is a plume associated depletion. We also address comments on the MHD simulations and demonstrate that 540-1040 keV losses are not necessarily inconsistent with 115-244 keV losses by plume associated charge exchange.

We present a comprehensive analysis of the compact planetary nebula M2-31 investigating its spectral properties, spatio-kinematical structure and chemical composition using GTC MEGARA integral field spectroscopic observations and NOT ALFOSC medium-resolution spectra and narrow-band images. The GTC MEGARA high-dispersion observations have remarkable tomographic capabilities, producing an unprecedented view of the morphology and kinematics of M2-31 that discloses a fast spectroscopic bipolar outflow along position angles 50$^\circ$ and 230$^\circ$, an extended shell and a toroidal structure or waist surrounding the central star perpendicularly aligned with the fast outflows. These observations also show that the C II emission is confined in the central region and enclosed by the [N II] emission. This is the first time that the spatial segregation revealed by a 2D map of the C II line implies the presence of multiple plasma components. The deep NOT ALFOSC observations allowed us to detect broad WR features from the central star of M2-31, including previously undetected broad O VI lines that suggest a reclassification as a [WO4]-type star.

Imaging atmospheric Cherenkov telescopes are continuously exposed to varying weather conditions that have short and long-term effects on their response to Cherenkov light from extensive air showers. This work presents the implementation of a throughput calibration method for the VERITAS telescopes taking into account changes in the optical response and detector performance over time. Different methods to measure the total throughput of the instrument, which depend on mirror reflectivites and PMT camera gain and efficiency, are discussed as well as the effect of its evolution on energy thresholds, effective collection areas, and energy reconstruction. The application of this calibration in the VERITAS data analysis chain is discussed, including the validation using Monte Carlo simulations and observations of the Crab Nebula.

In this present article, we study different accretion properties regarding viscous accretion of dark energy. Modified Chaplygin gas is chosen as the dark energy candidate. Viscosity is encountered with the help of Shakura-Sunyaev viscosity parameter. We study sonic speed vs radial distance curves. We compare between adiabatic and dark energy dominated cases and follow that sonic speed falls as we go nearer to the central gravitating object. As viscosity is imposed, a threshold drop in accretion sonic speed is followed. Average rate of fall in accretion sonic speed is increased with black hole's spin. This is signifying that this kind of accretion is weakening the overall matter/energy infall. Specific angular momentum to Keplerian angular momentum ratio is found to fall as we go far from the black hole. Accretion Mach number turns high as we go towards the inner region and high wind Mach number is not allowed as we are going out. Combining, we conclude that the system weakens the feeding process of accretion.

Numerous protoplanetary discs show distinct spiral arms features. While possibly caused by a range of processes, detailed pattern analysis points at close stellar flybys as cause for some of them. Surprisingly, these discs reside in young low-mass clusters, where close stellar flybys are expected to be rare. This fact motivated us to take a fresh look at the frequency of close flybys in low-mass clusters. In the solar neighbourhood, low-mass clusters have smaller half-mass radii than their more massive counterparts. We show that this observational fact results in the mean and central stellar density of low-mass clusters being approximately the same as in high-mass clusters, which is rarely reflected in theoretical studies. We perform N-body simulations of the stellar dynamics in young clusters obeying the observed mass-radius relation. Taking the mean disc truncation radius as a proxy for the degree of influence of the environment, we find that the influence of the environment on discs is more or less the same in low- and high-mass clusters. Even the fraction of small discs($<$ 10 au) is nearly identical. Our main conclusion is that the frequency of close flybys seems to have been severely underestimated for low-mass clusters. A testable prediction of this hypothesis is that low-mass clusters should contain 10%-15% of discs smaller than 30 au truncated by flybys. These truncated discs should be distinguishable from primordially small discs by their steep outer edge.

The presence of dark matter (DM) is suggested by a wealth of astrophysical and cosmological measurements. However, its underlying nature is yet unknown. Among the most promising candidates are weakly interacting massive particles (WIMPs): particles with mass and coupling strength at the electroweak scale and thermally produced in the early universe have a present relic density consistent with that observed today. WIMP self-annihilation would produce Standard Model particles including gamma-rays, which have been long-time recognized as a prime messenger to indirectly detect dark matter signals. The centre of the Milky Way is predicted as the brightest source of DM annihilations. The H.E.S.S. collaboration is currently performing a survey of the inner region of the Milky Way, the Inner Galaxy Survey (IGS), intended to achieve the best sensitivity to faint and diffuse emissions in a region of several degrees around the Galactic Centre. We analyzed 2014-2020 observations taken with the five-telescope array to search for a DM annihilation signal. With the current dataset of about 550 hours, we found no significant excess and therefore derived strong constraints on the velocity-weighted annihilation cross-section. TeV thermal WIMPs can be probed in different annihilation channels.

We present the discovery of a low-mass comoving system found by means of the NOIRLab Source Catalog (NSC) DR2. The system consists of the high proper-motion star LEHPM 5005 and an ultracool companion 2MASS J22410186-4500298 with an estimated spectral type of L2. The primary (LEHPM 5005) is likely a mid-M dwarf but over-luminous for its color, indicating a possible close equal mass binary. According to the Gaia EDR3 parallax of the primary, the system is located at a distance of $58\pm2$ pc. We calculated an angular separation of 7.2" between both components, resulting in a projected physical separation of 418 AU.

We consider the possibility that the gauge hierarchy is a byproduct of the metastability of the electroweak vacuum, i.e., that whatever mechanism is responsible for the latter also sets the running Higgs mass to a value smaller than its natural value by many orders of magnitude. This perspective is motivated by the early-time framework for eternal inflation put forth recently, which favors vacua that are relatively short-lived, but applies more generally to any theoretical approach predicting that our vacuum should be metastable. We find that the metastability of the electroweak vacuum, together with the requirement that such a non-trivial vacuum exists, requires the Higgs mass to be smaller than the instability scale by around one order of magnitude. While this bound is quite weak in the Standard Model (SM), as the instability scale is $\sim 10^{11}$ GeV, simple and well-motivated extensions of the SM - concretely, the $\nu$MSM with an approximate $B-\tilde{L}$ symmetry and the minimal SU(4)/Sp(4) composite Higgs model - can significantly tighten the bound by lowering the instability scale. We find that the bound can be brought down to $\simeq 10$ TeV where our perturbative treatment of the decay rate becomes unreliable. Our results imply that, assuming the SM symmetry breaking pattern, small running Higgs masses are a universal property of theories giving rise to metastability, suggesting a common origin of the two underlying fine-tunings and providing a strong constraint on any attempt to explain metastability.

The first-order phase transitions in the early universe are one of the well-known sources which release the stochastic background of gravitational waves (GWs). In this paper, we study the contribution of an external static and strong magnetic field on the stochastic background of gravitational waves (GWs) expected during QCD phase transition. In the light of the strongly magnetized hot QCD Equation of State which deviated from the ideal gas up to one-loop approximation, we estimate two phenomenologically important quantities: peak-frequency redshifted to today ($f_{\rm peak}$) and GW strain amplitude ($h^2 \Omega_{gw}$). The trace anomaly induced by the magnetized hot QCD matter around phase transition generates the stochastic background of GW with the peak-frequencies lower than the ideal gas-based signal (around nHz). Instead, the strain amplitudes corresponding to the peak frequencies are of the same order of magnitudes of the expected signal from ideal gas. This may be promising in the sense that although the strong magnetic field could mask the expected stochastic background of GWs but merely by upgrading the frequency sensitivity of detectors in the future, the magnetized GW is expected to be identified. Faced with the projected reach of detectors EPTA, IPTA, and SKA, we find that for the tail of the magnetized GW signals there remains a mild possibility of detection as it can reach the projected sensitivity of SKA.

We propose a novel, frequency-domain approach to the analysis of the gravitational-wave ringdown signal of binary black holes and the identification of quasinormal mode frequencies of the remnant. Our approach avoids the issues of spectral leakage that would normally be expected (associated with the abrupt start of the ringdown) by modeling the inspiral and merger parts of the signal using a flexible sum of sine-Gaussian wavelets truncated at the onset of the ringdown. Performing the analysis in the frequency domain allows us to use standard (and by now well-established) Bayesian inference pipelines for gravitational wave data as well as giving us the ability to readily search over the sky position and the ringdown start time, although we find that it is necessary to use an informative prior for the latter. We test our method by using it to analyze several simulated signals with varying signal-to-noise ratios injected into two- and three-detector networks. We find that our frequency-domain approach is generally able to place tighter constraints on the remnant black-hole mass and spin than a standard time-domain analysis.

We investigate the thermodynamics of FRW (Friedmann-Robertson-Walker) universe in the extended phase space. We generalize the unified first law with a cosmological constant $\Lambda$ by using the Misner-Sharp energy. We treat the cosmological constant as the thermodynamic pressure of the system, and derive thermodynamic equation of state $P = P(V, T)$ for the FRW universe. To clarify our general result, we present two applications of this thermodynamic equation of state, including Joule-Thomson expansions and efficiency of the Carnot heat engines. These investigations lead to physical insights of the evolution of the universe in view of thermodynamics.

Gravitational lensing is one of the most impressive celestial phenomena, which has interesting behaviors in its strong field limit. Near such limit, Bozza finds that the deflection angle of light is well-approximated by a logarithmic term and a constant term. In this way he explicitly derived the analytic expressions of deflection angles for a few types of black holes. In this paper, we study the explicit calculation to two new types of metrics in the strong field limit: (i) the Schwarzschild metric extended with an additional $r^{-n}(n\geq 3)$ term in the metric function; (ii) the Reissner-Nordstrom metric extended with an additional $r^{-6}$ term in the metric function. With such types of metrics, Bozza's original way of choosing integration variables may lead to technical difficulties in explicitly expressing the deflection angles, and we use a slightly modified version of Bozza's method to circumvent the problem.

In this paper, we study the thermodynamics especially the $P$-$V$ criticality of the Friedmann-Robertson-Walker (FRW) universe in the novel 4-dimensional Gauss-Bonnet gravity, where we define the thermodynamic pressure $P$ from the cosmological constant $\Lambda$ as $P=-\frac{\Lambda}{8\pi}$. We obtain the first law of thermodynamics and equation of state of the FRW universe. We find that, if the Gauss-Bonnet coupling constant $\alpha$ is positive, there is no $P$-$V$ phase transition. If $\alpha$ is negative, there are $P$-$V$ phase transitions and critical behaviors within $-1/3\leq\omega\leq1/3$. Particularly, there are two critical points of the $P$-$V$ criticality in the case $\alpha<0,~-1/3<\omega<1/3$. We investigate these $P$-$V$ criticality around the critical points, and calculate the critical exponents. We find that these critical exponents in the $-1/3<\omega\leq1/3$ case are consistent with those in the mean field theory, and hence satisfy the scaling laws.

A mechanism for the formation of primordial black holes is proposed. Here, heavy quarks of a confining gauge theory produced by de Sitter fluctuations are pushed apart by inflation and get confined after horizon re-entry. The large amount of energy stored in the colour flux tubes connecting the quark pair leads to black-hole formation. These are much lighter and can be of higher spin than those produced by standard collapse of horizon-size inflationary overdensities. Other difficulties exhibited by such mechanisms are also avoided. Phenomenological features of the new mechanism are discussed as well as accounting for both the entirety of the dark matter and the supermassive black holes in the galactic centres. Under proper conditions, the mechanism can be realised in a generic confinement theory, including ordinary QCD. We discuss a possible string-theoretic realisation via $D$-branes. Interestingly, for conservative values of the string scale, the produced gravity waves are within the range of recent NANOGrav events. Simple generalisations of the mechanism allow for the existence of a significant scalar component of gravity waves with distinct observational signatures.

A canonical description of a corotating solar wind high speed stream, in terms of velocity profile, would indicate three main regions:a stream interface or corotating interaction region characterized by a rapid flow speed increase and by compressive phenomena due to dynamical interaction between the fast wind flow and the slower ambient plasma;a fast wind plateau characterized by weak compressive phenomena and large amplitude fluctuations with a dominant Alfv\'enic character;a rarefaction region characterized by a decreasing trend of the flow speed and wind fluctuations dramatically reduced in amplitude and Alfv\'enic character, followed by the slow ambient wind. Interesting enough, in some cases the region where the severe reduction of these fluctuations takes place is remarkably short in time, of the order of minutes, and located at the flow velocity knee separating the fast wind plateau from the rarefaction region. The aim of this work is to investigate which are the physical mechanisms that might be at the origin of this phenomenon. We firstly looked for the presence of any tangential discontinuity which might inhibit the propagation of Alfv\'enic fluctuations from fast wind region to rarefaction region. The absence of a clear evidence for the presence of this discontinuity between these two regions led us to proceed with ion composition analysis for the corresponding solar wind, looking for any abrupt variation in minor ions parameters (as tracers of the source region) which might be linked to the phenomenon observed in the wind fluctuations. In the lack of a positive feedback from this analysis, we finally propose a mechanism based on interchange reconnection experienced by the field lines at the base of the corona, within the region separating the open field lines of the coronal hole, source of the fast wind, from the surrounding regions mainly characterized by closed field lines.

We review results from QCD axion string and domain wall simulations and propagate the associated uncertainties into the calculation of the axion relic density. This allows us to compare different results in the literature and, using cosmological constraints, to perform statistical inference on the axion mass window in the post-inflationary Peccei-Quinn symmetry breaking scenario. For dark matter axions, this leads to a median dark matter axion mass of 0.50 meV, while the 95% credible interval at highest posterior density is between 0.48 and 0.52 meV. For simulations including string-domain wall decays, these numbers are 0.22 meV and [0.16, 0.27] meV. Relaxing the condition that axions are all of the dark matter, the axion mass window is completed by an upper bound of around 80 meV, which comes from hot dark matter constraints. This demonstrates, at least from the statistical perspective, that the axion mass can be constrained rather precisely once it is possible to overcome the much larger systematic uncertainties.

The synthesis of heavy, proton rich isotopes is a poorly understood astrophysical process. Thermonuclear (type Ia) supernova explosions are among the suggested sites and the abundance of some isotopes present in the early solar system may be used to test the models. 92Nb is such an isotope and one of the reactions playing a role in its synthesis is 91Zr(p,gamma)92Nb. As no experimental cross sections were available for this reaction so far, nucleosynthesis models had to solely rely on theoretical calculations. In the present work the cross section of 91Zr(p,gamma)92mNb has been measured at astrophysical energies by activation. The results excellently confirm the predictions of cross sections and reaction rates for 91Zr(p,gamma)92Nb, as used in astrophysical simulations.

Tomographic imaging with radionuclides commonly used in nuclear medicine, such as $^{111}$In (171 and 245 keV) and $^{131}$I (364 keV), is in high demand for medical applications and small animal imaging. The Si/CdTe Compton camera with its high angular and high energy resolutions is an especially promising detector to extend the energy coverage for imaging to the range that covers gamma-ray emitted from these radionuclides. Here, we take the first steps towards short-distance imaging by conducting experiments using three-dimensional phantoms composed of multiple sphere-like solutions of $^{111}$In and $^{131}$I with a diameter of 2.7 mm, placed at a distance of 41 mm. Using simple back-projection methods, the positions of the sources are reproduced with a spatial resolution of 11.5 mm and 9.0 mm (FWHM) for $^{111}$In and $^{131}$I, respectively. We found that a LM-MLEM method gives a better resolution of 4.0 mm and 2.7 mm (FWHM). We resolve source positions of a tetrahedron structure with a source-to-source separation of 28 mm. These findings demonstrate that Compton Cameras have the potential of close-distance imaging of radioisotopes distributions in the energy range below 400 keV.

The spatial distribution and polarization of Saturn narrowband (NB) emissions have been studied by using Cassini Radio and Plasma Wave Sciences data and goniopolarimetric data obtained through an inversion algorithm with a preset source located at the center of Saturn. From 2004 January 1 to 2017 September 12, NB emissions were selected automatically by a computer program and rechecked manually. The spatial distribution shows a preference for high latitude and intensity peaks in the region within 6 Saturn radii for both 5 and 20 kHz NB emissions. 5 kHz NB emissions also show a local time preference roughly in the 18:00-22:00 sector. The Enceladus plasma torus makes it difficult for NB emissions to propagate to the low latitude regions outside the plasma torus. The extent of the low latitude regions where 5 and 20 kHz NB emissions were never observed is consistent with the corresponding plasma torus density contour in the meridional plane. 20 kHz NB emissions show a high circular polarization while 5 kHz NB emissions are less circularly polarized with |V|<0.6 for majority of the cases. And cases of 5kHz NB emissions with high circular polarization are more frequently observed at high latitude especially at the northern and southern edges of the Enceladus plasma torus.

We consider modified gravity cosmological models that can be transformed into two-field chiral cosmological models by the conformal metric transformation. For the $R^2$ gravity model with an additional scalar field and the corresponding two-field model with the cosmological constant and nonstandard kinetic part of the action, the general solutions have been obtained in the spatially flat FLRW metric. We analyze the correspondence of the cosmic time solutions obtained and different possible evolutions of the Hubble parameters in the Einstein and Jordan frames.

I selectively review the theoretical properties and observational limits pertaining to dark atoms, as well as composite dark matter candidates bound by a confining gauge interaction: dark glueballs, glueballinos, mesons and baryons. Emphasis is given to cosmological, direct and indirect detection constraints.