Papers for Wednesday, Apr 28 2021

This paper is an extension of the paper by Del Popolo, Chan, and Mota (2020) to take account the effect of dynamical friction. We show how dynamical friction changes the threshold of collapse, $\delta_c$, and the turn-around radius, $R_t$. We find numerically the relationship between the turnaround radius, $R_{\rm t}$, and mass, $M_{\rm t}$, in $\Lambda$CDM, in dark energy scenarios, and in a $f(R)$ modified gravity model. Dynamical friction gives rise to a $R_{\rm t}-M_{\rm t}$ relation differing from that of the standard spherical collapse. In particular, dynamical friction amplifies the effect of shear, and vorticity already studied in Del Popolo, Chan, and Mota (2020). A comparison of the $R_{\rm t}-M_{\rm t}$ relationship for the $\Lambda$CDM, and those for the dark energy, and modified gravity models shows, that the $R_{\rm t}-M_{\rm t}$ relationship of the $\Lambda$CDM is similar to that of the dark energy models, and small differences are seen when comparing with the $f(R)$ models. The effect of shear, rotation, and dynamical friction is particularly evident at galactic scales, giving rise to a difference between the $R_{\rm t}-M_{\rm t}$ relation of the standard spherical collapse of the order of $\simeq 60\%$. Finally, we show how the new values of the $R_{\rm t}-M_{\rm t}$ influence the constraints to the $w$ parameter of the equation of state.

Using our measurements of the H$\alpha$ emission line flux originating in the cool (T $\sim10^4$ K) gas that populates the halos of galaxies, we build a joint model to describe mass of the cool circumgalactic medium (CGM) as a function of galactic stellar mass ($10^{9.5} < ({\rm M_*/M}_\odot) < 10^{11}$) and environment. Because the H$\alpha$ emission correlates with the main cooling channel for this gas, we are able to estimate the rate at which the CGM cools and becomes fuel for star formation in the central galaxy. We describe this calculation, which uses our observations, previous measurements of some critical CGM properties, and modeling of the cooling mechanism using the \cloudy modeling suite. We find that the mass cooling rate is larger than the star formation rates of the central galaxies by a factor of $\sim 4 - 90$, empirically confirming that there is sufficient fuel to resolve the gas consumption problem and that feedback is needed to avoid collecting too much cold gas in galaxies. We find excellent agreement between our estimates of both the mass cooling rates and mass loading factors and the predictions of independent theoretical studies. The convergence in results that we find from several completely different treatments of the problem, particularly at the lower end of the galactic mass range, is a strong indication that we have a relatively robust understanding of the quantitative effects of feedback across this mass range.

Most of massive stars form in binary or higher-order systems in clumpy, sub-structured clusters. In the very first phases of their life, these stars are expected to interact with the surrounding environment, before being released to the field when the cluster is tidally disrupted by the host galaxy. We present a set of N-body simulations to describe the evolution of young stellar clusters and their binary content in the first phases of their life. To do this, we have developed a method that generates realistic initial conditions for binary stars in star clusters from hydrodynamical simulations. We considered different evolutionary cases to quantify the impact of binary and stellar evolution. Also, we compared their evolution to that of King and fractal models with different length scales. Our results indicate that the global expansion of the cluster from hydrodynamical simulations is initially balanced by the sub-clump motion and accelerates when a monolithic shape is reached, as in a post-core collapse evolution. Compared to the spherical initial conditions, the ratio of the 50% to 10% Lagrangian radius shows a very distinctive trend, explained by the formation of a hot core of massive stars triggered by the high initial degree of mass segregation. As for its binary population, each cluster shows a self-regulating behaviour by creating interacting binaries with binding energies of the order of its energy scales. Also, in absence of original binaries, the dynamically formed binaries present a mass dependent binary fraction, that mimics the trend of the observed one.

Satellites are not randomly distributed around their central galaxies but show polar and planar structures. In this paper, we investigate the axis-asymmetry or lopsidedness of satellite galaxy distributions around isolated galaxies in a hydrodynamic cosmological simulation. We find a statistically significant lopsided signal by studying the angular distribution of the satellite galaxies' projected positions around isolated central galaxies in a two-dimensional plane. The signal is dependent on galaxy mass, color and large-scale environment. Satellites that inhabit low-mass blue hosts, or located further from the hosts show the most lopsided signal. Galaxy systems with massive neighbors exhibit stronger lopsidedness. This satellite axis-asymmetry signal also decreases as the universe evolves. Our findings are in agreement with recent observational results, and they provide a useful perspective for studying galaxy evolution, especially on the satellite accretion, internal evolution and interaction with the cosmic large-scale structure.

We post-process galaxies in the IllustrisTNG simulations with SKIRT radiative transfer calculations to make predictions for the rest-frame near-infrared (NIR) and far-infrared (FIR) properties of galaxies at $z\geq 4$. The rest-frame $K$- and $z$-band galaxy luminosity functions from TNG are overall consistent with observations, despite a $\sim 0.4\,\mathrm{dex}$ underprediction at $z=4$ for $M_{\rm z}\lesssim -24$. Predictions for the JWST MIRI observed galaxy luminosity functions and number counts are given. We show that the next-generation survey conducted by JWST can detect 500 (30) galaxies in F1000W in a survey area of $500\,{\rm arcmin}^{2}$ at $z=6$ ($z=8$). As opposed to the consistency in the UV, optical and NIR, we find that TNG, combined with our dust modelling choices, significantly underpredicts the abundance of most dust-obscured and thus most luminous FIR galaxies. As a result, the obscured cosmic star formation rate density (SFRD) and the SFRD contributed by optical/NIR dark objects are underpredicted. The discrepancies discovered here could provide new constraints on the sub-grid feedback models, or the dust contents, of simulations. Meanwhile, although the TNG predicted dust temperature and its relations with IR luminosity and redshift are qualitatively consistent with observations, the peak dust temperature of $z\geq 6$ galaxies are overestimated by about $20\,{\rm K}$. This could be related to the limited mass resolution of our simulations to fully resolve the porosity of the interstellar medium (or specifically its dust content) at these redshifts.

We develop a new machine learning algorithm, Via Machinae, to identify cold stellar streams in data from the Gaia telescope. Via Machinae is based on ANODE, a general method that uses conditional density estimation and sideband interpolation to detect local overdensities in the data in a model agnostic way. By applying ANODE to the positions, proper motions, and photometry of stars observed by Gaia, Via Machinae obtains a collection of those stars deemed most likely to belong to a stellar stream. We further apply an automated line-finding method based on the Hough transform to search for line-like features in patches of the sky. In this paper, we describe the Via Machinae algorithm in detail and demonstrate our approach on the prominent stream GD-1. A companion paper contains our identification of other known stellar streams as well as new stellar stream candidates from Via Machinae. Though some parts of the algorithm are tuned to increase sensitivity to cold streams, the Via Machinae technique itself does not rely on astrophysical assumptions, such as the potential of the Milky Way or stellar isochrones. This flexibility suggests that it may have further applications in identifying other anomalous structures within the Gaia dataset, for example debris flow and globular clusters.

We present a spectroscopic determination of the redshift of the second source in the Jackpot gravitational lens system J0946+1006, for which only a photometric estimate of $z_{\rm phot}$ = 2.41$^{+0.04}_{-0.21}$ has previously been available. By visually inspecting an archival VLT X-Shooter observation, we located a single emission line from the source in the H-band. Among the possible options we find that this line is most likely to be [OIII] 5007 Ang at $z_{\rm spec} $ = 2.035. Guided by this proposal, we were able to detect the faint CIII] 1907,1909 Ang emission doublet in a deep VLT MUSE datacube. The CIII] emission is spatially coincident with the brightest parts of the second Einstein ring, and strongly supports the redshift identification. The spectroscopic redshift is only marginally consistent with the photometric estimate. Re-examining the cosmological constraints from J0946+1006, the revised measurement favours less negative values of the dark energy equation-of-state parameter $w$; when combined with a cosmic microwave background prior, we infer $w$ = $-1.04\pm0.20$. The revised redshift does not significantly help to reconcile the small discrepancy in the image positions for the even more distant third source in J0946+1006.

We propose a new method for probing inflationary models of primordial black hole (PBH) production, using only CMB physics at relatively large scales. In these scenarios, the primordial power spectrum profile for curvature perturbations is characterized by a pronounced dip, followed by a rapid growth towards small scales, leading to a peak responsible for PBH formation. We focus on scales around the dip that are well separated from the peak to analytically compute expressions for the curvature power spectrum and bispectrum. The size of the squeezed bispectrum is enhanced at the position of the dip, and it acquires a characteristic scale dependence that can be probed by cross-correlating CMB $\mu$-distortions and temperature fluctuations. We quantitatively study the properties of such cross-correlations and how they depend on the underlying model, discussing how they can be tested by the next generation of CMB $\mu$-distortion experiments. This method allows one to experimentally probe inflationary PBH scenarios using well-understood CMB physics, without considering non-linearities associated with PBH formation and evolution.

In Dou et al. (2021), we introduced the Fundamental Formation Relation (FFR), a tight relation between specific SFR (sSFR), H$_2$ star formation efficiency (SFE$_{\rm H_2}$), and the ratio of H$_2$ to stellar mass. Here we show that atomic gas HI does not follow a similar FFR as H$_2$. The relation between SFE$_{\rm HI}$ and sSFR shows significant scatter and strong systematic dependence on all of the key galaxy properties that we have explored. The dramatic difference between HI and H$_2$ indicates that different processes (e.g., quenching by different mechanisms) may have very different effects on the HI in different galaxies and hence produce different SFE$_{\rm HI}$-sSFR relations, while the SFE$_{\rm H_2}$-sSFR relation remains unaffected. The facts that SFE$_{\rm H_2}$-sSFR relation is independent of other key galaxy properties, and that sSFR is directly related to the cosmic time and acts as the cosmic clock, make it natural and very simple to study how different galaxy populations (with different properties and undergoing different processes) evolve on the same SFE$_{\rm H_2}$-sSFR $\sim t$ relation. In the gas regulator model (GRM), the evolution of a galaxy on the SFE$_{\rm H_2}$-sSFR($t$) relation is uniquely set by a single mass-loading parameter $\lambda_{\rm net,H_2}$. This simplicity allows us to accurately derive the H$_2$ supply and removal rates of the local galaxy populations with different stellar masses, from star-forming galaxies to the galaxies in the process of being quenched. This combination of FFR and GRM, together with the stellar metallicity requirement, provide a new powerful tool to study galaxy formation and evolution.

We develop a new phenomenological model that addresses current tensions between observations of the early and late Universe. Our scenario features: (i) a decaying dark energy fluid, which undergoes a transition at $z \sim 5,000$, to raise today's value of the Hubble parameter -- addressing the $H_0$ tension, and (ii) an ultra-light axion, which starts oscillating at $z\sim 16,000$, to suppress the matter power spectrum -- addressing the $S_8$ tension. Our Markov Chain Monte Carlo analyses show that such a Dark Sector model fits a combination of early time datasets slightly better than the $\Lambda$CDM model, while reducing both the $H_0$ and $S_8$ tensions to $\lesssim 3\sigma$ level. Combined with measurements from cosmic shear surveys, we find that the discrepancy on $S_8$ is reduced to the $1.4\sigma$ level, and the value of $H_0$ is further raised. Adding local supernovae measurements, we find that the $H_0$ and $S_8$ tensions are reduced to the $1.5\sigma$ and $1.1\sigma$ level respectively, with a significant improvement $\Delta\chi^2\simeq -17$ compared to the $\Lambda$CDM model. We discuss a possible particle physics realization of this model, with a dark confining gauge sector and its associated axion, although embedding the full details within microphysics remains an urgent open question. Our scenario will be decisively probed with future CMB surveys.

We present results on global very long baseline interferometry (VLBI) observations at 327 MHz of eighteen compact steep-spectrum (CSS) and GHz-peaked spectrum (GPS) radio sources from the 3C and the Peacock & Wall catalogues. About 80 per cent of the sources have a 'double/triple' structure. The radio emission at 327 MHz is dominated by steep-spectrum extended structures, while compact regions become predominant at higher frequencies. As a consequence, we could unambiguously detect the core region only in three sources, likely due to self-absorption affecting its emission at this low frequency. Despite their low surface brightness, lobes store the majority of the source energy budget, whose correct estimate is a key ingredient in tackling the radio source evolution. Low-frequency VLBI observations able to disentangle the lobe emission from that of other regions are therefore the best way to infer the energetics of these objects. Dynamical ages estimated from energy budget arguments provide values between 2x10^3 and 5x10^4 yr, in agreement with the radiative ages estimated from the fit of the integrated synchrotron spectrum, further supporting the youth of these objects. A discrepancy between radiative and dynamical ages is observed in a few sources where the integrated spectrum is dominated by hotspots. In this case the radiative age likely represents the time spent by the particles in these regions, rather than the source age.

We have followed up two ultra-diffuse galaxies (UDGs), detected adjacent to stellar streams, with Hubble Space Telescope (HST) imaging and HI mapping with the Jansky Very Large Array (VLA) in order to investigate the possibility that they might have a tidal origin. With the HST F814W and F555W images we measure the globular cluster (GC) counts for NGC 2708-Dw1 and NGC 5631-Dw1 as $2^{+1}_{-1}$ and $5^{+1}_{-2}$, respectively. NGC 2708-Dw1 is undetected in HI down to a 3$\sigma$ limit of $\log (M_\mathrm{HI}/\mathrm{M_\odot}) = 7.3$, and there is no apparent HI associated with the nearby stellar stream. There is a 2$\sigma$ HI feature coincident with NGC 5631-Dw1. However, this emission is blended with a large gaseous tail emanating from NGC 5631 and is not necessarily associated with the UDG. The presence of any GCs and the lack of clear HI connections between the UDGs and their parent galaxies strongly disfavor a tidal dwarf galaxy origin, but cannot entirely rule it out. The GC counts are consistent with those of normal dwarf galaxies, and the most probable formation mechanism is one where these UDGs were born as normal dwarfs and were later tidally stripped and heated. We also identify an over-luminous ($M_\mathrm{V} = -11.1$) GC candidate in NGC 2708-Dw1, which may be a nuclear star cluster transitioning to an ultra-compact dwarf as the surrounding dwarf galaxy gets stripped of stars.

We compared high-contrast near-infrared images of the core of R136 taken by VLT/SPHERE, in two epochs separated by 3.06 years. For the first time we monitored the dynamics of the detected sources in the core of R136 from a ground-based telescope with adaptive optics. The aim of these observations was to search for High prOper Motion cAndidates (HOMAs) in the central region of R136 (r<6") where it has been challenging for other instruments. Two bright sources (K<15mag and V<16mag) are located near R136a1 and R136c (massive WR stars) and have been identified as potential HOMAs. These sources have significantly shifted in the images with respect to the mean shift of all reliable detected sources and their neighbours, and six times their own astrometric errors. We calculate their proper motions to be 1.36\pm0.22 mas/yr (321\pm52 km/s) and 1.15\pm0.11 mas/yr (273\pm26 km/s). We discuss different possible scenarios to explain the magnitude of such extreme proper motions, and argue for the necessity to conduct future observations to conclude on the nature of HOMAs in the core of R136.

We compute successfully the launching of two magnetic winds from two circumbinary disks formed after a common envelope event. The launching is produced by the increase of magnetic pressure due to the collapse of the disks. The collapse is due to internal torques produced by a weak poloidal magnetic field. The first wind can be described as a wide jet, with an average mass-loss rate of $\sim 1.3 \times 10^{-7}$ \Moy\ and a maximum radial velocity of $\sim 230$ \kms. The outflow has a half-opening angle of $\sim 20^{\circ}$. Narrow jets are also formed intermittently with velocities up to 3,000 \kms, with mass-loss rates of $\sim 6 \times 10^{-12} $ \Moy\ during short periods of time. The second wind can be described as a wide X-wind, with an average mass-loss rate of $\sim 1.68 \times 10^{-7}$ \Moy\ and a velocity of $\sim 30$ \kms. A narrow jet is also formed with a velocity of 250 \kms, and a mass-loss rates of $\sim 10^{-12} $ \Moy. The computed jets are used to provide inflow boundary conditions for simulations of proto-planetary nebulae. The wide jet evolves into a molecular collimated outflow within a few astronomical units, producing proto-planetary nebulae with bipolar, elongated shapes, whose kinetic energies reach $\sim 4 \times 10^{45}$ erg at 1,000 years. Similarities with observed features in W43A, OH231.8+4.2, and Hen 3-1475 are discussed. The computed wide X-wind produces proto-planetary nebulae with slower expansion velocities, with bipolar and elliptical shapes, and possible starfish type and quadrupolar morphology.

Observatory publications comprise the work of local astronomers from observatories around the world and are traditionally exchanged between observatories through libraries. However, large collections of observatory publications seem to be rare; or at the least rarely digitally described or accessible on the Internet. Notable examples to the contrary are the Woodman Astronomical Library at Wisconsin-Madison and the Dudley Observatory in Loudonville, New York both in the US. Due to the irregularities in receiving material, the collections are generally often incomplete both with respect to the observatories included as well as volumes. In order to assess the unique properties of the collections, we summarize and compare observatories present in our own as well as the collections from the Woodman Library and the Dudley Observatory.

In the current paradigm of planet formation research, it is believed that the first step to forming massive bodies (such as asteroids and planets) requires that small interstellar dust grains floating through space collide with each other and grow to larger sizes. The initial formation of these pebbles is governed by an integro-differential equation known as the Smoluchowski coagulation equation, to which analytical solutions are intractable for all but the simplest possible scenarios. While brute-force methods of approximation have been developed, they are computationally costly, currently making it infeasible to simulate this process including other physical processes relevant to planet formation, and across the very large range of scales on which it occurs. In this paper, we take a machine learning approach to designing a system for a much faster approximation. We develop a multi-output random forest regression model trained on brute-force simulation data to approximate distributions of dust particle sizes in protoplanetary disks at different points in time. The performance of our random forest model is measured against the existing brute-force models, which are the standard for realistic simulations. Results indicate that the random forest model can generate highly accurate predictions relative to the brute-force simulation results, with an $R^{2}$ of 0.97, and do so significantly faster than brute-force methods.

Coronal Mass Ejections (CMEs) are highly dynamic events originating in the solar atmosphere, that show a wide range of kinematic properties and are the major drivers of the space weather. The angular width of the CMEs is a crucial parameter in the study of their kinematics. The fact that whether slow and fast CMEs (as based on their relative speed to the average solar wind speed) are associated with different processes at the location of their ejection is still debatable. Thus, in this study, we investigate their angular width to understand the differences between the slow and fast CMEs. We study the width distribution of slow and fast CMEs and find that they follow different power law distributions, with a power law indices ($\alpha$) of -1.1 and -3.7 for fast and slow CMEs respectively. To reduce the projection effects, we further restrict our analysis to only limb events as derived from manual catalog and we find similar results. We then associate the slow and fast CMEs to their source regions, and classified the sources as Active Regions (ARs) and Prominence Eruptions (PEs). We find that slow and fast CMEs coming from ARs and PEs, also follow different power laws in their width distributions. This clearly hints towards a possibility that different mechanisms might be involved in the width expansion of slow and fast CMEs coming from different sources.These results are also crucial from the space weather perspective since the width of the CME is an important factor in that aspect.

We present results from \textit{Chandra} X-ray observations and 325 MHz Giant Metrewave Radio Telescope (GMRT) observations of the massive and X-ray luminous cluster of galaxies Abell S1063. We report the detection of large-scale \lq\lq excess brightness\rq\rq\ in the residual \textit{Chandra} X-ray surface brightness map, which extends at least 2.7 Mpc towards the north-east from the center of the cluster. We also present a high fidelity X-ray flux and temperature map using \textit{Chandra} archival data of 122 ksec, which shows the disturbed morphology in the cluster. The residual flux map shows the first observational confirmation of the merging axis proposed by earlier simulation by \citet{Gomez2012AJ....144...79G}. The average temperature within $R_{500}$ is $11.7 \pm 0.56$ keV, which makes AS1063 one of the hottest clusters in the nearby Universe. The integrated radio flux density at 325 MHz is found to be $62.0\pm6.3$ mJy. The integrated spectrum of the radio halo follows a power-law with a spectral index $\alpha=-1.43\pm 0.13$. The radio halo is found to be significantly under-luminous, which favored for both the hadronic as well as the turbulent re-acceleration mechanism for its origin.

We investigate the origin of $\Lambda$CDM parameter constraints in weak lensing, with a focus on the Hubble constant. We explain why current cosmic shear data are sensitive to the parameter combination $S_8 \propto \sigma_8 \Omega_m^{0.5}$, improving upon previous studies through use of the halo model. Motivated by the ongoing discrepancy in measurements of the Hubble constant from high and low redshift, we explain why cosmic shear provides almost no constraint on $H_0$ by showing how the lensing angular power spectrum depends on physical length scales in the dark matter distribution. We derive parameter constraints from galaxy lensing in KiDS and cosmic microwave background weak lensing from Planck and SPTpol, separately and jointly, showing how degeneracies between $\sigma_8$ and $\Omega_m$ can be broken. Using lensing to calibrate the sound horizon measured in projection by baryon acoustic oscillations gives $H_0 = 67.4 \pm 0.9 \; \mathrm{km} \, \mathrm{s}^{-1} \, \mathrm{Mpc}^{-1}$, consistent with previous results from Planck and the Dark Energy Survey. We find that a toy Euclid-like lensing survey provides only weak constraints on the Hubble constant due to degeneracies with other parameters that affect the shape of the lensing correlation functions. If external priors on $n_s$, the baryon density, and the amplitude of baryon feedback are available then sub-percent $H_0$ constraints are achievable with forthcoming lensing surveys.

Comoving pairs, even at the separations of $\mathcal{O}(10^6)\,$AU, are a predicted reservoir of conatal stars. We present detailed chemical abundances of 62 stars in 31 comoving pairs with separations of $10^2 - 10^7\,$AU and 3D velocity differences $< 2 \mathrm{\ km \ s^{-1}}$. This sample includes both bound comoving pairs/wide binaries and unbound comoving pairs. Observations were taken using the MIKE spectrograph on the Magellan/Clay Telescope at high resolution ($\mathrm{R} \sim 45,000$) with a typical signal-to-noise ratio of 150 per pixel. With these spectra, we measure surface abundances for 24 elements, including Li, C, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Sr, Y, Zr, Ba, La, Nd, Eu. Taking iron as the representative element, our sample of wide binaries is chemically homogeneous at the level of $0.05$ dex, which agrees with prior studies on wide binaries. Importantly, even systems at separations $2\times10^5-10^7\,$AU are homogeneous to $0.09$ dex, as opposed to the random pairs which have a dispersion of $0.23\,$dex. Assuming a mixture model of the wide binaries and random pairs, we find that $73 \pm 22\%$ of the comoving pairs at separations $2\times10^5-10^7\,$AU are conatal. Our results imply that a much larger parameter space of phase space may be used to find conatal stars, to study M-dwarfs, star cluster evolution, exoplanets, chemical tagging, and beyond.

We performed timing and spectral analysis of multi-epoch Suzaku, XMM-Newton and NuSTAR observations of the ultraluminous X-ray source (ULX) Circinus ULX5, to put constraints on the mass of the central object and the accretion mode operating in this source. We aim to answer whether the source contains a stellar mass black hole with a super-Eddington accretion flow or an intermediate mass black hole accreting matter in a sub-Eddington mode. Moreover, we search for major observed changes in spectra and timing occur, and determine if they are associated with major structural changes in the disk, similar to those in black hole X-ray binaries. We performed timing and spectral analysis to study the relation between luminosity and inner disk temperature. We constructed the hardness ratio versus intensity diagram to reveal spectral state transitions in Circinus ULX5. Our spectral analysis revealed at least three distinctive spectral states of Circinus ULX5, in analogy to state transitions in Galactic black hole X-ray binaries. Disk-dominated spectra are found in high flux states and the power-law dominated spectra are found in lower flux states. The source was also observed in an intermediate state, where the flux was low, but spectrum is dominated by a disk component. Over eighteen years of collected data, ULX5 appeared two times in the high, three times in the low, and two times in the intermediate state. The fastest observed transition was $\sim$7 months. Our analysis suggests that the central object in Circinus ULX5 is a stellar mass BH ($<10\ \rm M_{\odot}$), or possibly a NS even though we do not detect pulsations in the lightcurves. Fractional variability amplitudes are consistent with state transitions in Circinus ULX5 wherein higher variability from the power law-like Comptonized emission gets suppressed in the thermal disk-dominated state.

Supermassive binary black holes (SMBBHs) are laboratories par excellence for relativistic effects, including precession effects in the Kerr metric and the emission of gravitational waves. Binaries form in the course of galaxy mergers, and are a key component in our understanding of galaxy evolution. Dedicated searches for SMBBHs in all stages of their evolution are therefore ongoing and many systems have been discovered in recent years. Here we provide a review of the status of observations with a focus on the multiwavelength detection methods and the underlying physics. Finally, we highlight our ongoing, dedicated multiwavelength program MOMO (for Multiwavelength Observations and Modelling of OJ 287). OJ 287 is one of the best candidates to date for hosting a sub-parsec SMBBH. The MOMO program carries out a dense monitoring at >13 frequencies from radio to X-rays and especially with Swift since 2015. Results so far included: (1) The detection of two major UV-X-ray outbursts with Swift in 2016/17 and 2020; exhibiting softer-when-brighter behaviour. The non-thermal nature of the outbursts was clearly established and shown to be synchrotron radiation. (2) Swift multi-band dense coverage and XMM-Newton spectroscopy during EHT campaigns caught OJ 287 at an intermediate flux level with synchrotron and IC spectral components. (3) Discovery of a remarkable, giant soft X-ray excess with XMM and NuSTAR during the 2020 outburst. (4) Spectral evidence (at 2sigma) for a relativistically shifted iron absorption line in 2020. (5) The non-thermal 2020 outburst is consistent with an after-flare predicted by the SMBBH model of OJ 287.

We present the results obtained from detailed timing and spectral studies of a black hole candidate MAXI~J1813-095 using {\it Swift}, {\it NICER}, and {\it NuSTAR} observations during its 2018 outburst. The timing behaviour of the source is mainly studied by using {\it NICER} light curves in a $0.5-10$ keV range. We did not find any signature of quasi-periodic oscillations in the power density spectra of the source. We carry out spectral analysis with a combined disk blackbody \& power-law model, and a physical two-component advective flow (TCAF) model. From the combined {\tt disk blackbody} \& {\tt power-law} model, we extracted thermal and non-thermal fluxes, photon index, and inner disk temperature. We also find evidence for weak reflection in the spectra. We have tested the physical TCAF model on a broadband spectrum from {\it NuSTAR} and {\it Swift}/XRT. The parameters like mass accretion rates, the size of the Compton clouds and the shock strength are extracted. Our result shows that the source remained in the hard state during the entire outburst which indicates a `failed' outburst. We estimate the mass of the black hole as $7.4 \pm 1.5$ $M_{\odot}$ from the spectral study with the TCAF model. We use {\tt LAOR} model for the Fe K$\alpha$ line emission. From this, the spin parameter of the black hole is estimated as $a^* > 0.76$. The inclination angle of the system is estimated to be in the range of $28^{\circ} - 45^{\circ}$ from the reflection model. We estimate the source distance to be $\sim 6$ kpc.

The degree of porosity in interstellar dust-grain material is poorly defined, although recent work has suggested that the grains could be highly porous. Aside from influencing the optical properties of the dust, porosity has the potential to affect the chemistry occurring on dust-grain surfaces, via increased surface area, enhanced local binding energies, and the possibility of trapping of molecules within the pores as ice mantles build up on the grains. Through computational kinetics simulations, we investigate how interstellar grain-surface chemistry and ice composition are affected by the porosity of the underlying dust-grain material. Using a simple routine, idealized three-dimensional dust-grains are constructed, atom by atom, with varying degrees of porosity. Diffusive chemistry is then simulated on these surfaces using the off-lattice microscopic Monte Carlo chemical kinetics model, MIMICK, assuming physical conditions appropriate to dark interstellar clouds. On the porous grain surface, the build-up of ice mantles, mostly composed of water, leads to the covering over of the pores, leaving empty pockets. Once the pores are completely covered, the chemical and structural behavior is similar to non-porous grains of the same size. The most prominent chemical effect of the presence of grain porosity is the trapping of molecular hydrogen, formed on the grain surfaces, within the ices and voids inside the grain pores. Trapping of H2 in this way may indicate that other volatiles, such as inert gases not included in these models, could be trapped within dust-grain porous structures when ices begin to form.

We present SNIascore, a deep-learning based method for spectroscopic classification of thermonuclear supernovae (SNe Ia) based on very low-resolution (R $\sim100$) data. The goal of SNIascore is fully automated classification of SNe Ia with a very low false-positive rate (FPR) so that human intervention can be greatly reduced in large-scale SN classification efforts, such as that undertaken by the public Zwicky Transient Facility (ZTF) Bright Transient Survey (BTS). We utilize a recurrent neural network (RNN) architecture with a combination of bidirectional long short-term memory and gated recurrent unit layers. SNIascore achieves a $<0.6\%$ FPR while classifying up to $90\%$ of the low-resolution SN Ia spectra obtained by the BTS. SNIascore simultaneously performs binary classification and predicts the redshifts of secure SNe Ia via regression (with a typical uncertainty of $<0.005$ in the range from $z = 0.01$ to $z = 0.12$). For the magnitude-limited ZTF BTS survey ($\approx70\%$ SNe Ia), deploying SNIascore reduces the amount of spectra in need of human classification or confirmation by $\approx60\%$. Furthermore, SNIascore allows SN Ia classifications to be automatically announced in real-time to the public immediately following a finished observation during the night.

Analyzes of next-generation galaxy data require accurate treatment of systematic effects such as the bias between observed galaxies and the underlying matter density field. However, proposed models of the phenomenon are either numerically expensive or too inaccurate to achieve unbiased inferences of cosmological parameters even at mildly-nonlinear scales of the data. As an alternative to constructing accurate galaxy bias models, requiring understanding galaxy formation, we propose to construct likelihood distributions for Bayesian forward modeling approaches that are insensitive to linear, scale-dependent bias and provide robustness against model misspecification. We use maximum entropy arguments to construct likelihood distributions designed to account only for correlations between data and inferred quantities. By design these correlations are insensitive to linear galaxy biasing relations, providing the desired robustness. The method is implemented and tested within a Markov Chain Monte Carlo approach. The method is assessed using a halo mock catalog based on standard full, cosmological, N-body simulations. We obtain unbiased and tight constraints on cosmological parameters exploiting only linear cross-correlation rates for $k\le 0.10$ Mpc/h. Tests for halos of masses ~10$^{12}$ M$_\odot$ to ~10$^{13}$ M$_\odot$ indicate that it is possible to ignore all details of the linear, scale dependent, bias function while obtaining robust constraints on cosmology. Our results provide a promising path forward to analyzes of galaxy surveys without the requirement of having to accurately model the details of galaxy biasing but by designing robust likelihoods for the inference.

We present detailed timing and spectral studies of the black hole candidate MAXI J0637$-$430 during its 2019-2020 outburst using observations with the {\em Neutron Star Interior Composition Explorer (NICER)} and the {\em Neil Gehrels Swift Observatory}. We find that the source evolves through the soft-intermediate, high-soft, hard-intermediate and low-hard states during the outburst. No evidence of quasi-periodic oscillations is found in the power density spectra of the source. Weak variability with fractional rms amplitude $<5\%$ is found in the softer spectral states. In the hard-intermediate and hard states, high variability with the fractional rms amplitude of $>20\%$ is observed. The $0.7-10$ keV spectra with {\em NICER} are studied with a combined disk-blackbody and nthcomp model along with the interstellar absorption. The temperature of the disc is estimated to be $0.6$ keV in the rising phase and decreased slowly to $0.1$ keV in the declining phase. The disc component was not detectable or absent during the low hard state. From the state-transition luminosity and the inner edge of the accretion flow, we estimate the mass of the black hole to be in the range of 5$-$12 $M_{\odot}$, assuming the source distance of $d<10$ kpc.

Simulating the irradiation of planetary atmospheres by cosmic ray particles requires, among others, the ability to understand and to quantify the interactions of charged particles with planetary magnetic fields. Here we present a process that is very often ignored in such studies; the dispersion and focusing of cosmic ray trajectories in magnetospheres. The calculations were performed using our new code CosmicTransmutation, which has been developed to study cosmogenic nuclide production in meteoroids and planetary atmospheres and which includes the computation of the irradiation spectrum on top of the atmosphere. Here we discuss effects caused by dispersion and focusing of cosmic ray particle trajectories.

We present five simultaneous UV/X-ray observations of IC4329A by AstroSat, performed over {a five-month} period. We utilize the excellent spatial resolution of the Ultra-Violet Imaging Telescope (UVIT) onboard AstroSat to reliably separate the intrinsic AGN flux from the host galaxy emission and to correct for the Galactic and internal reddening, as well as the contribution from the narrow and broad-line regions. We detect large-amplitude UV variability, which is unusual for a large black hole mass AGN, like IC4329A, over such a small period. In fact, the fractional variability amplitude is larger in the UV band than in the X-ray band. This demonstrates that the observed UV variability is intrinsic to the disk, and is not due to X-ray illumination. The joint X-ray spectral analyses of five SXT and LAXPC spectral data reveal a soft-X-ray excess component, a narrow iron-line (with no indication of a significant Compton hump), and a steepening power-law ($\Delta\Gamma\sim 0.21$) with increasing X-ray flux. The soft excess component could arise due to thermal Comptonization of the inner disk photons in a warm corona with $kT_e\sim 0.26$ keV. The UV emission we detect acts as the primary seed photons for the hot corona, which produces the broadband X-ray continuum. The X-ray spectral variability is well described by the cooling of this corona from $kT_e\sim42$ keV to $\sim 32$ keV with increasing UV flux, while the optical depth remains constant at $\tau\sim 2.3$.

The direct imaging of rocky exoplanets is one of the major science goals for upcoming large telescopes. The contrast requirement for imaging such planets is challenging. However, the mid-IR (InfraRed) regime provides the optimum contrast to directly detect the thermal signatures of exoplanets in our solar neighbourhood. We aim to exploit novel fast chopping techniques newly developed for astronomy with the aid of adaptive optics to look for thermal signatures of exoplanets around bright stars in the solar neighbourhood. We use the upgraded VISIR (Very Large Telescope Imager and Spectrometer for the mid-InfraRed) instrument with high contrast imaging (HCI) capability optimized for observations at 10~$\mu$m to look for exoplanets around five nearby ($d$ < 4 pc) stars. The instrument provides an improved signal-to-noise (S/N) by a factor of $\sim$4 in the N-band compared to standard VISIR for a given S/N and time. In this work we achieve a detection sensitivity of sub-mJy, which is sufficient to detect few Jupiter mass planets in nearby systems. Although no detections are made we achieve most sensitive limits within $<2''$ for all the observed targets compared to previous campaigns. For $\epsilon$ Indi A and $\epsilon$ Eri we achieve detection limits very close to the giant planets discovered by RV, with the limits on $\epsilon$ Indi A being the most sensitive to date. Our non-detection therefore supports an older age for $\epsilon$ Indi A. The results presented here show the promise for high contrast imaging and exoplanet detections in the mid-IR regime.

We trained deep neural networks (DNNs) as a function of the neutrino energy density, flux, and the fluid velocity to reproduce the Eddington tensor for neutrinos obtained in our first-principles core-collapse supernova (CCSN) simulations. Although the moment method, which is one of the most popular approximations for neutrino transport, requires a closure relation, none of the analytical closure relations commonly employed in the literature captures all aspects of the neutrino angular distribution in momentum space. In this paper, we developed a closure relation by using the DNN that takes the neutrino energy density, flux, and the fluid velocity as the input and the Eddington tensor as the output. We consider two kinds of DNNs: a conventional DNN named a component-wise neural network (CWNN) and a tensor-basis neural network (TBNN). We found that the diagonal component of the Eddington tensor is reproduced better by the DNNs than the M1-closure relation especially for low to intermediate energies. For the off-diagonal component, the DNNs agree better with the Boltzmann solver than the M1 closure at large radii. In the comparison between the two DNNs, the TBNN has slightly better performance than the CWNN. With the new closure relations at hand based on the DNNs that well reproduce the Eddington tensor with much smaller costs, we opened up a new possibility for the moment method.

This study uses HI image data from the WALLABY pilot survey with the ASKAP telescope, covering the Hydra cluster out to 2.5$r_{200}$. We present the projected phase-space distribution of HI-detected galaxies in Hydra, and identify that nearly two thirds of the galaxies within $1.25r_{200}$ may be in the early stages of ram pressure stripping. More than half of these may be only weakly stripped, with the ratio of strippable HI (i.e., where the galactic restoring force is lower than the ram pressure in the disk) mass fraction (over total HI mass) distributed uniformly below 90%. Consequently, the HI mass is expected to decrease by only a few 0.1 dex after the currently strippable portion of HI in these systems has been stripped. A more detailed look at the subset of galaxies that are spatially resolved by WALLABY observations shows that, while it typically takes less than 200 Myr for ram pressure stripping to remove the currently strippable portion of HI, it may take more than 600 Myr to significantly change the total HI mass. Our results provide new clues to understanding the different rates of HI depletion and star formation quenching in cluster galaxies.

We investigate the possibility of the growth of magnetorotational instability (MRI) in disks around Class 0 protostars. We construct a disk model and calculate the chemical reactions of neutral and charged atoms, molecules and dust grains to derive the abundance of each species and the ionization degree of the disk. Then, we estimate the diffusion coefficients of non-ideal magnetohydrodynamics effects such as ohmic dissipation, ambipolar diffusion and the Hall effect. Finally, we evaluate the linear growth rate of MRI in each area of the disk. We investigate the effect of changes in the strength and direction of the magnetic field in our disk model and we adopt four different dust models to investigate the effect of dust size distribution on the diffusion coefficients. Our results indicate that an MRI active region possibly exists with a weak magnetic field in a region far from the protostar where the Hall effect plays a role in the growth of MRI. On the other hand, in all models the disk is stable against MRI in the region within $<20$ au from the protostar on the equatorial plane. Since the size of the disks in the early stage of star formation is limited to $\lesssim 10-$$20$ au, it is difficult to develop MRI-driven turbulence in such disks.

We present the discovery of another Odd Radio Circle (ORC) with the Australian Square Kilometre Array Pathfinder (ASKAP) at 944 MHz. The observed radio ring, ORC J0102-2450, has a diameter of ~70 arcsec or 300 kpc, if associated with the central elliptical galaxy DES J010224.33-245039.5 (z ~ 0.27). Considering the overall radio morphology (circular ring and core) and lack of ring emission at non-radio wavelengths, we investigate if ORC J0102-2450 could be the relic lobe of a giant radio galaxy seen end-on or the result of a giant blast wave. We also explore possible interaction scenarios, for example, with the companion galaxy, DES J010226.15-245104.9, located in or projected onto the south-eastern part of the ring. We encourage the search for further ORCs in radio surveys to study their properties and origin.

We present the low-resolution transmission spectra of the puffy hot Jupiter HAT-P-65b (0.53 M$_\mathrm{Jup}$, 1.89 R$_\mathrm{Jup}$, $T_\mathrm{eq}=1930$ K), based on two transits observed using the OSIRIS spectrograph on the 10.4 m Gran Telescopio CANARIAS (GTC). The transmission spectra of the two nights are consistent, covering the wavelength range 517--938 nm and consisting of mostly 5 nm spectral bins. We perform equilibrium-chemistry spectral retrieval analyses on the jointly fitted transmission spectrum and obtain an equilibrium temperature of $1645^{+255}_{-244}$ K and a cloud coverage of $36^{+23}_{-17}$%, revealing a relatively clear planetary atmosphere. Based on free-chemistry retrieval, we report strong evidence for TiO. Additional individual analyses in each night reveal weak-to-moderate evidence for TiO in both nights, but moderate evidence for Na or VO only in one of the nights. Future high-resolution Doppler spectroscopy as well as emission observations will help confirm the presence of TiO and constrain its role in shaping the vertical thermal structure of HAT-P-65b's atmosphere.

Context: In the context of the solar coronal heating problem, one possible explanation for the high coronal temperature is the release of energy by magnetohydrodynamic (MHD) waves. The energy transfer is believed to be possible, among others, by the development of the Kelvin-Helmholtz instability (KHI) in coronal loops. Aims: Our aim is to determine if standing slow waves in solar atmospheric structures such as coronal loops, and also prominence threads, sunspots, and pores, can trigger the KHI due to the oscillating shear flow at the structure's boundary. Methods: We used linearized nonstationary MHD to work out an analytical model in a cartesian reference frame. The model describes a compressible plasma near a discontinuous interface separating two regions of homogeneous plasma, each harboring an oscillating velocity field with a constant amplitude which is parallel to the background magnetic field and aligned with the interface. The obtained analytical results were then used to determine the stability of said interface, both in coronal and photospheric conditions. Results: We find that the stability of the interface is determined by a Mathieu equation. In function of the parameters of this equation, the interface can either be stable or unstable. For coronal as well as photospheric conditions, we find that the interface is stable with respect to the KHI. Theoretically, it can, however, be unstable with respect to a parametric resonance instability, although it seems physically unlikely. We conclude that, in this simplified setup, a standing slow wave does not trigger the KHI without the involvement of additional physical processes.

We present the first comprehensive spectroscopic study of the Andromeda galaxy's Eastern Extent. This ~4 degree long filamentary structure. located 70-90 kpc from the centre of M31, lies perpendicular to Andromeda's minor axis and the Giant Stellar Stream and overlaps Stream C. In this work, we explore the properties of the Eastern Extent to look for possible connections between it, the Giant Stellar Stream and Stream C. We present the kinematics and photometry for ~50 red giant branch stars in 7 fields along the Eastern Extent. We measure the systemic velocities for these fields and find them to be -368 km/s < v < -331 km/s with a slight velocity gradient of -0.51 +/- 0.21 km/s/kpc towards the Giant Stellar Stream. We derive the photometric metallicities for stars in the Eastern Extent finding them to be metal-poor with values of -1.0 < [Fe/H]phot < -0.7 with a <[Fe/H]phot> ~-0.9. We find consistent properties for the Eastern Extent, Stream B and one of the substructures in Stream C, Stream Cr, plausibly linking these features. Stream Cp and its associated globular cluster, EC4, have distinctly different properties indicative of a separate structure. When we compare the properties of the Eastern Extent to those of the Giant Stellar Stream, we find them to be consistent, albeit slightly more metal-poor, such that the Eastern Extent could plausibly comprise stars stripped from the progenitor of the Giant Stellar Stream.

Cosmic-ray transport in astrophysical environments is often dominated by the diffusion of particles in a magnetic field composed of both a turbulent and a mean component. This process needs to be understood in order to properly model cosmic-ray signatures. One of the most important aspects in the modeling of cosmic-ray diffusion is that fully resonant scattering, the most effective such process, is only possible if the wave spectrum covers the entire range of propagation angles. By taking the wave spectrum boundaries into account, we quantify cosmic-ray diffusion parallel and perpendicular to the guide field direction at turbulence levels above 5% of the total magnetic field. We apply our results of the parallel and perpendicular diffusion coefficient to the Milky Way. We show that simple purely diffusive transport is in conflict with observations of the inner Galaxy, but that just by taking a Galactic wind into account, data can be matched in the central 5 kpc zone. Further comparison shows that the outer Galaxy at $>5$ kpc, on the other hand, should be dominated by perpendicular diffusion, likely changing to parallel diffusion at the outermost radii of the Milky Way.

NASA's Spitzer Infrared Spectrometer (IRS) acquired mid-infrared (5-37 microns) disc-averaged spectra of Uranus very near to its equinox in December 2007. A mean spectrum was constructed from observations of multiple central meridian longitudes, spaced equally around the planet, which has provided the opportunity for the most comprehensive globally-averaged characterisation of Uranus' temperature and composition ever obtained (Orton et al., 2014 a [arXiv:1407.2120], b [arXiv:1407.2118]). In this work we analyse the disc-averaged spectra at four separate central meridian longitudes to reveal significant longitudinal variability in thermal emission occurring in Uranus' stratosphere during the 2007 equinox. We detect a variability of up to 15% at wavelengths sensitive to stratospheric methane, ethane and acetylene at the ~0.1-mbar level. The tropospheric hydrogen-helium continuum and deuterated methane absorption exhibit a negligible variation (less than 2%), constraining the phenomenon to the stratosphere. Building on the forward-modelling analysis of the global average study, we present full optimal estimation inversions (using the NEMESIS retrieval algorithm, Irwin et al., 2008 [10.1016/j.jqsrt.2007.11.006]) of the Uranus-2007 spectra at each longitude to distinguish between thermal and compositional variability. We found that the variations can be explained by a temperature change of less than 3 K in the stratosphere. Near-infrared observations from Keck II NIRC2 in December 2007 (Sromovsky et al., 2009 [arXiv:1503.01957], de Pater et al., 2011 [10.1016/j.icarus.2011.06.022]), and mid-infrared observations from VLT/VISIR in 2009 (Roman et al., 2020 [arXiv:1911.12830]), help to localise the potential sources to either large scale uplift or stratospheric wave phenomena.

We investigate magnetic activity properties of 21 stars via medium resolution optical spectra and long-term photometry. Applying synthetic spectrum fitting method, we find that all targets are cool giant or sub-giant stars possessing overall [M/H] abundances between $0$ and $-0.5$. We find that six of these targets exhibit only linear trend in mean brightness while eight of them clearly shows cyclic mean brightness variation. Remaining seven target appear to exhibit cyclic mean brightness variation but this can not be confirmed due to the long time scales of the predicted cycle compared to the current time range of the photometric data. We further determine seasonal photometric periods and compute average photometric period of each target. Analysed sample in this study provides a quantitative representation of a positive linear correlation between the inverse of the rotation period and the cycle period normalized to the rotation period, on the log-log scale. We also observe no correlation between the activity cycle length and the relative surface shear, indicating that the activity cycle must be driven by a parameter rather than the differential rotation. Our analyses show that the relative surface shear is positively correlated with the rotation period and there is a noticeable separation between main sequence stars and our sample. Compared to our sample, the relative surface shear of a main sequence star is larger for a given rotation period. However, dependence of the relative surface shear on the rotation period appears stronger for our sample. Analysis of the current photometric data indicates that the photometric properties of the observed activity cycles in 8 targets seem dissimilar to the sunspot cycle.

Context: Vela X-1 is one of the best studied X-ray binaries. Frequently though, specific values for its parameters have been used in subsequent studies without considering alternatives. Aims: We aim to provide a robust compilation and synthesis of the accumulated knowledge about Vela X-1 as a solid baseline for future studies and identify specific avenues of possible future research. Methods: We explore the literature for Vela X-1 and on modelling efforts, describing the evolution of the system knowledge. We also add information derived from public data, especially the Gaia EDR3 release. Results: We update the distance to Vela X-1, the spectral classification for HD 77518 and find that the supergiant may be very close to filling its Roche lobe. Constraints on the clumpiness of the stellar wind have improved. The orbit is very well determined, but the uncertain inclination limits information on the neutron star mass. Estimates for the stellar wind have evolved towards lower velocities, supporting the idea of transient wind-captured disks around the neutron star. Hydrodynamic models and observations are consistent with an accretion wake trailing the neutron star. Conclusions: Vela X-1 is an excellent laboratory, but a lot of room remains to improve. Well-coordinated multi-wavelength observations and campaigns addressing the intrinsic variability are required. New opportunities will arise through new instrumentation. Models of the stellar wind should account for the orbital eccentricity and the non-spherical shape of HD 77581. Realistic multi-dimensional models of radiative transfer in the UV and X-rays are needed, but remain very challenging. Improved MHD models covering a wide range of scales would be required to improve understanding of the plasma-magnetosphere coupling. A full characterization of the accretion column remains another open challenge. (Abbreviated for arXiv)

We present a full analysis of a broadband spectral line survey of Sagittarius B2 (Main), one of the most chemically rich regions in the Galaxy located within the giant molecular cloud complex Sgr B2 in the Central Molecular Zone. Our goal is to derive the molecular abundances and temperatures of the high-mass star-forming region Sgr B2(M) and thus its physical and astrochemical conditions. Sgr B2(M) was observed using the Heterodyne Instrument for the Far-Infrared (HIFI) on board the Herschel Space Observatory in a spectral line survey from 480 to 1907 GHz at a spectral resolution of 1.1 MHz, which provides one of the largest spectral coverages ever obtained toward this high-mass star-forming region in the submillimeter with high spectral resolution and includes frequencies > 1 THz unobservable from the ground. We model the molecular emission from the submillimeter to the far-IR using the XCLASS program. For each molecule, a quantitative description was determined taking all emission and absorption features of that species across the entire spectral range into account. Additionally, we derive velocity resolved ortho / para ratios for those molecules for which ortho and para resolved molecular parameters are available. Finally, the temperature and velocity distributions are analyzed and the derived abundances are compared with those obtained for Sgr B2(N) from a similar HIFI survey. A total of 92 isotopologues were identified, arising from 49 different molecules, ranging from free ions to complex organic compounds and originating from a variety of environments from the cold envelope to hot and dense gas within the cores. Sulfur dioxide, methanol, and water are the dominant contributors. For the ortho / para ratios we find deviations from the high temperature values between 13 and 27 %. In total 14 % of all lines remain unidentified.

In 3 previous papers I showed that the series of the midpoints of the times of all the X-ray flares of Sgr A* that were detected so far harbors a statistical trend termed pacemaker regularity. It means that X-ray flares are detected more frequently around time points that are subset of a periodic grid on the time axis of period P_X=0.1032 day=149 min. The series of the times of detection of the peaks of NIR flares of the object are also regulated by a pacemaker with a period P_IR=0.028 day=41 min. Here I show that the series of the midpoints of the times of recorded NIR flare is also regulated by a pacemaker of the period P_IRM=0.039 day=56min. The 2 pacemakers found in the previous papers were interpreted as signals of a star that revolves around the blackhole of Sgr A* in orbit with a mean radius of ~3.2 Schwarzschild radii of the blackhole, here corrected to ~3.13. The finding of the period of the third pacemaker is consistent with the suggested revolving star model. Here I present the specific orbit of the star as well as a plausible description of its sidereal rotation. The model also implies that the star has an unusual internal structure. I show that the discovery of the GRAVITY collaboration of motion of hot spots at distances from the blackhole that are of the order of very few Schwarzschild radii of it may well be understood within the context of the revolving star model.

As a rule, both lobes of Fanaroff-Riley (FR) type-II radio sources are terminated with hotspots, but the 3C328 radio galaxy is a specimen of an FR II-like object with a hotspot in only one lobe. A conceivable reason for such asymmetry is that the nucleus of 3C328 was temporarily inactive. There was no energy transfer from it to the lobes during the period of quiescence, and so they began to fade out. However, under the assumption that the axis connecting the two lobes makes an appreciable angle with the sky plane, and hence one is considerably farther from the observer than the other, the lobes are observed at two distinct stages of evolution due to the light-travel lag. While the far-side lobe is still perceived as being of the FR II type with a hotspot, decay of the near-side lobe is already apparent. No jets are visible in the VLA images, but the VLBA observations of the inverted-spectrum core component of 3C328 have revealed that it has a jet of a sub-arcsecond length pointing towards the lobe that shows evidence of decay. Since the jet always points to the near side, its observed orientation is in line with the scenario proposed here. The presence of the jet supports the inference that the nucleus of 3C328 is currently active; however, given the fact that the jet is short (approx. 200 pc in projection), the activity must have restarted very recently. The lower and upper limits of the quiescent period length have been calculated.

This is the third of a series of papers that presents an algorithm to search for close binaries with massive, possibly compact, unseen secondaries. The detection of such a binary is based on identifying a star that displays a large ellipsoidal periodic modulation, induced by tidal interaction with its companion. In the second paper of the series we presented a simple approach to derive a robust modified minimum mass ratio (mMMR), based on the observed ellipsoidal amplitude, without knowing the primary mass and radius, assuming the primary fills its Roche lobe. The newly defined mMMR is always smaller than the actual mass ratio. Therefore, a binary with an mMMR larger than unity is a good candidate for having a massive secondary, which might be a black hole or a neutron star. This paper considers 10,956 OGLE short-period ellipsoidals observed towards the Galactic Bulge. We re-analyse their modulation and identify 136 main-sequence systems with mMMR significantly larger than unity as candidates for having compact-object secondaries, assuming their observed periodic modulations reflect indeed the ellipsoidal effect. Obviously, one needs follow-up observations to find out the true nature of these companions.

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

Recent progress in constraining the massive accretions (>1:10) experienced by the Milky Way (MW) and the Andromeda galaxy (M31) offers an opportunity to understand the dwarf galaxy population of the Local Group. Using zoom-in dark matter-only simulations of MW-mass haloes and concentrating on subhaloes that are thought to be capable of hosting dwarf galaxies, we demonstrate that the infall of a massive progenitor is accompanied with the accretion and destruction of a large number of subhaloes. Massive accreted progenitors do not increase the total number of infalling subhaloes onto a MW-mass host, but instead focus surrounding subhaloes onto the host causing a clustering in the infall time of subhaloes. This leads to a temporary elevation in the number of subhaloes as well as changes in their cumulative radial profile within the virial radius of the host. Surviving associated subhaloes with a massive progenitor have a large diversity in their orbits. We find that the star formation quenching times of Local Group dwarf spheroidal galaxies ($10^{5} \mathrm{M_{\odot}} \lesssim \mathrm{M}_{*} \lesssim 10^{7} \mathrm{M_{\odot}}$) are clustered around the times of the most massive accretions suffered by the MW and M31. Our results imply that a) the quenching time of dwarf spheroidals is a good proxy of their infall time and b) the absence of recently quenched satellites around M31 suggests that M33 is not on its first infall and was accreted much earlier.

The underground muon detector of the Pierre Auger Observatory is aimed at attaining direct measurements of the muonic component of extensive air showers produced by cosmic rays with energy from $10^{16.5}$ eV up to the region of the ankle (around $10^{18.7}$ eV). It consists of two nested triangular grids of underground scintillators with 433 m, and 750 m spacings and a total of 71 positions, each with 192 scintillator strips (30 m$^2$) deployed 2.3 m underground. The light produced by impinging muons in the scintillators is propagated with optical fibers towards an array of silicon photomultipliers. In this work, we present the development, validation, and performance of an end-to-end tool for simulating the response of the underground muon detector to single-muon signals, which constitutes the basis for further simulations of the whole array. Laboratory data and simulation outcomes are found consistent, showing that with the underground muon detector we can measure single muons, with an efficiency of 99 %, up to about 1050 particles arriving at exactly the same time in 30 m$^2$ of scintillator.

The SRG observatory, equipped with the X-ray telescopes Mikhail Pavlinsky ART-XC and eROSITA, was launched by Roscosmos to the L2 point on July 13, 2019. The launch was carried out from Baikonur by a Proton-M rocket with a DM-03 upper stage. The German telescope eROSITA was installed on SRG under agreement between Roskosmos and DLR. In December 2019, SRG started to scan the celestial sphere in order to obtain X-ray maps of the entire sky in several energy bands (from 0.3 to 8 keV, eROSITA, and from 4 to 30 keV, ART-XC). By mid-December 2020, the second full-sky scan had been completed. Over 4 years, 8 independent maps of the sky will be obtained. Their sum will reveal more than three million quasars and over one hundred thousand galaxy clusters and groups. The availability of 8 sky maps will enable monitoring of long-term variability (every six months) of a huge number of extragalactic and Galactic X-ray sources, including hundreds of thousands of stars. Rotation of the satellite around the axis directed toward the Sun with a period of 4 hours makes it possible to track faster variability of bright X-ray sources. The chosen scanning strategy leads to the formation of deep survey zones near both ecliptic poles. We present sky maps obtained by the telescopes aboard SRG during the first scan of the sky and a number of results of deep observations performed during the flight to L2, demonstrating the capabilities of the Observatory in imaging, spectroscopy and timing. In December 2023 the Observatory will switch for at least two years to observations of the most interesting sources in the sky in triaxial orientation mode and deep scanning of selected fields with an area of up to 150 sq. deg. These modes of operation were tested during the Performance Verification phase. Every day, SRG data are dumped onto the largest antennae of the Russian Deep Space Network in Bear Lakes and near Ussuriysk.

As a follow-up of the optical spectroscopic campaign aimed at achieving completeness in the Third Catalog of Hard Fermi-LAT Sources (3FHL), we present here the results of a sample of 28 blazars of uncertain type observed using the 4m telescope at Cerro Tololo Inter-American Observatory (CTIO) in Chile. Out of these 28 sources, we find that 25 are BL Lacertae objects (BL Lacs) and 3 are Flat Spectrum Radio Quasars (FSRQs). We measure redshifts or lower limits for 16 of these blazars, whereas it is observed that the 12 remaining blazars have featureless optical spectra. These results are part of a more extended campaign of optical spectroscopy follow-up of 3FHL blazars, where until now 51 blazars of uncertain type have been classified into BL Lac and FSRQ categories. Further, this campaign has resulted in redshift measurements and lower limits for 15 of these sources. Our results contribute towards attaining a complete sample of blazars above 10 GeV, which then will be crucial in extending our knowledge on blazar emission mechanisms and the extragalactic background light.

Low-frequency radio selection finds radio-bright galaxies regardless of the amount of obscuration by gas and dust. We report \chandra\ observations of a complete 178~MHz-selected, and so orientation unbiased, sample of 44 $0.5<z<1$ 3CRR sources. The sample is comprised of quasars and narrow-line radio galaxies (NLRGs) with similar radio luminosities, and the radio structure serves as both an age and an orientation indicator. Consistent with Unification, intrinsic obscuration (measured by \nh, X-ray hardness ratio, and X-ray luminosity) generally increases with inclination. However, the sample includes a population not seen in high-$z$ 3CRR sources: NLRGs viewed at intermediate inclination angles with \nh~$<10^{22}$~cm$^{-2}$. Multiwavelength analysis suggests these objects have lower $L/L_{\rm Edd}$ than typical NLRGs at similar orientation. Thus both orientation and $L/L_{\rm Edd}$ are important, and a "radiation-regulated Unification" provides a better explanation of the sample's observed properties. In comparison with the 3CRR sample at $1<z<2$, our lower-redshift sample shows a higher fraction of Compton-thin NLRGs (45\% vs.\ 29\%) but similar Compton-thick fraction (20\%), implying a larger covering factor of Compton-thin material at intermediate viewing angles and so a more "puffed-up" torus atmosphere. We posit that this is due to a range of $L/L_{\rm Edd}$ extending to lower values in this sample. In contrast, at high redshifts the narrower range and high $L/L_{\rm Edd}$ values allowed orientation (and so simple Unification) to dominate the sample's observed properties.

Fuzzy dark matter (FDM) has been a promising alternative to standard cold dark matter. The model consists of ultralight bosons with mass $m_b \sim 10^{-22}$ eV and features a quantum-pressure-supported solitonic core that oscillates. In this work, we show that the soliton density oscillations persist even after significant tidal stripping of the outer halo. We report two intrinsic yet distinct timescales associated, respectively, with the ground-state soliton wavefunction $\tau_{00}$ and the soliton density oscillations $\tau_\text{soliton}$, obeying $\tau_\text{soliton} /\tau_{00} \simeq 2.3$. The central star cluster (SC) in Eridanus II has a characteristic timescale $\tau_\text{soliton} / \tau_\text{SC} \sim 2$ to $3$ that deviates substantially from unity. As a result, we demonstrate, both analytically and numerically with three-dimensional self-consistent FDM simulations, that the gravitational heating of the SC owing to soliton density oscillations is negligible irrespective of $m_b$. We also show that the subhalo mass function to form Eridanus II does not place a strong constraint on $m_b$. These results are contrary to the previous findings by Marsh & Niemeyer (2019).

Ultralight axions with axion-photon couplings $g_{a\gamma\gamma} \sim {\rm few} \times 10^{-11}$ GeV$^{-1}$ may resolve a number of astrophysical anomalies, such as unexpected ~TeV transparency, anomalous stellar cooling, and X-ray excesses from nearby neutron stars. We show, however, that such axions are severely constrained by the non-observation of X-rays from the magnetic white dwarf (MWD) RE J0317-853 using ~40 ks of data acquired from a dedicated observation with the Chandra X-ray Observatory. Axions may be produced in the core of the MWD through electron bremsstrahlung and then convert to X-rays in the magnetosphere. The non-observation of X-rays constrains the axion-photon coupling to $g_{a\gamma\gamma} \lesssim 5.5 \times 10^{-13} \sqrt{C_{a\gamma\gamma}/C_{aee}}$ GeV$^{-1}$ at 95% confidence for axion masses $m_a \lesssim 5 \times 10^{-6}$ eV, with $C_{aee}$ and $C_{a\gamma\gamma}$ the dimensionless coupling constants to electrons and photons. Considering that $C_{aee}$ is generated from the renormalization group, our results robustly disfavor $g_{a\gamma\gamma} \gtrsim 4.4 \times 10^{-11}$ GeV$^{-1}$ even for models with no ultraviolet contribution to $C_{aee}$.

The relic gravitational wave background due to tensor linear perturbations generated during Higgs inflation is computed. Both the Standard Model and a well-motivated phenomenological completion (that accounts for all the experimentally confirmed evidence of new physics) are considered. We focus on critical Higgs inflation, which improves on the non-critical version and features an amplification of the tensor fluctuations. The latter property allows us to establish that future space-borne interferometers, such as DECIGO, BBO and ALIA, may detect the corresponding primordial gravitational waves.

Using the holographic correspondence as a tool, we determine the steady-state velocity of expanding vacuum bubbles nucleated within chiral finite temperature first-order phase transitions occurring in strongly-coupled large $N$ QCD-like models. We provide general formulae for the friction force exerted by the plasma on the bubbles and for the steady-state velocity. In the top-down holographic description, the phase transitions are related to changes in the embedding of $Dq$-${\bar Dq}$ flavor branes probing the black hole background sourced by a stack of $N$ $Dp$-branes. We first consider the Witten-Sakai-Sugimoto $D4$-$D8$-$\bar D8$ setup, compute the friction force and deduce the equilibrium velocity. Then we extend our analysis to more general setups and to different dimensions. Finally, we briefly compare our results, obtained within a fully non-perturbative framework, to other estimates of the bubble velocity in the literature.

We have studied in the mechanical and chemical instabilities as well as the liquid-gas phase transition in isospin asymmetric quark matter based on the NJL and the pNJL model. Areas of the mechanical instability region and the liquid-gas coexistence region are seen to be enlarged with a larger quark matter symmetry energy or in the presence of strange quarks. Our study shows that the light cluster yield ratio observed in relativistic heavy-ion collisions may not be affected much by the uncertainties of the isospin effect, while the hadron-quark phase transition in compact stars as well as their mergers is likely to be a smooth one.

We construct Barrow holographic dark energy in the case of non-flat universe. In particular, considering closed and open spatial geometry we extract the differential equations that determine the evolution of the dark-energy density parameter, and we provide the analytical expression for the corresponding dark energy equation-of-state parameter. We show that the scenario can describe the thermal history of the universe, with the sequence of matter and dark energy epochs. Comparing to the flat case, where the phantom regime is obtained for relative large Barrow exponents, the incorporation of positive curvature leads the universe into the phantom regime for significantly smaller values. Additionally, in the case of negative curvature we find a reversed behavior, namely for increased Barrow exponent we acquire algebraically higher dark-energy equation-of-state parameters. Hence, the incorporation of slightly non-flat spatial geometry to Barrow holographic dark energy improves the phenomenology while keeping the new Barrow exponent to smaller values.