Papers for Tuesday, Apr 06 2021

We model the effects of stellar collisions and close encounters on the stellar populations observed in the Milky Way nuclear stellar cluster. Our analysis is based on $N$-body simulations in which the nuclear stellar cluster forms by accretion of massive stellar clusters in the presence of a supermassive black hole. We attach stellar populations to our $N$-body particles and follow the evolution of their stars and the rate of collisions and close encounters between different types of stars. We find that the most common encounters are collisions between pairs of main-sequence stars, which lead to merger; destructive collisions of stars with compact objects are relatively rare. We find that the effects of collisions on the stellar populations are small for three reasons. First, the core-like density profile in our $N$-body models limits the stellar density. Secondly, the velocity dispersion in the nuclear stellar cluster is similar to the surface escape velocities of the stars, which minimises the collision rate. Finally, whilst collisions between main-sequence stars can destroy bright giants by accelerating the evolution of their progenitors, they can also create them by accelerating the evolution of lower-mass stars. To a good approximation these two effects cancel out. We also investigate whether the G2 cloud that made a close pass by Sgr A* in 2014 could be a fuzzball: a compact stellar core which has accreted a tenuous envelope in a close encounter with a red giant. We conclude that fuzzballs with cores below about $2\,\msun$ have thermal times-scales that are too short to be compatible with observations of G2. A fuzzball with a black-hole core could reproduce the surface properties of G2 but the production rate of such objects in our model is low.

A binary neutron star (BNS) merger can lead to various outcomes, from indefinitely stable neutron stars, through supramassive (SMNS) or hypermassive (HMNS) neutron stars supported only temporarily against gravity, to black holes formed promptly after the merger. Up-to-date constraints on the BNS total mass and the neutron star equation of state suggest that a long-lived SMNS may form in $\sim 0.45-0.9$ of BNS mergers. We find that a SMNS typically needs to lose $\sim 3-6\times 10^{52}$ erg of it's rotational energy before it collapses, on a fraction of the spin-down timescale. A SMNS formation imprints on the electromagnetic counterparts to the BNS merger. However, a comparison with observations reveals tensions. First, the distribution of collapse times is too wide and that of released energies too narrow (and the energy itself too large) to explain the observed distributions of internal X-ray plateaus, invoked as evidence for SMNS-powered energy injection. Secondly, the immense energy injection into the blastwave should lead to extremely bright radio transients which previous studies found to be inconsistent with deep radio observations of short gamma-ray bursts. Furthermore, we show that upcoming all-sky radio surveys will enable to constrain the distribution of extracted energies, independently of a GRB jet formation. Our results can be self-consistently understood, provided that BNS merger remnants collapse shortly after their formation (even if their masses are low enough to allow for SMNS formation). We briefly outline how this collapse may be achieved. Future simulations will be needed to test this hypothesis.

Magnetic fields are amplified as a consequence of galaxy formation and turbulence-driven dynamos. Galaxy mergers can potentially amplify the magnetic fields from their progenitors, making the magnetic fields dynamically important. However, the effect of mergers on magnetic fields is still poorly understood. We use thermal polarized emission observations to trace the magnetic fields in the molecular disk of the nearest radio active galaxy, Centaurus A, which is thought to be the remnant of a merger. Here, we detect that the magnetic field orientations in the plane of the sky are tightly following the $\sim3.0$ kpc-scale molecular warped disk. Our simple regular large-scale axisymmetric spiral magnetic field model can explain, to some extent, the averaged magnetic field orientations across the disk projected on the sky. Our observations also suggest the presence of small-scale turbulent fields, whose relative strength increases with velocity dispersion and column density. These results have strong implications for understanding the generation and role of magnetic fields in the formation of galaxies across cosmic time.

Quantitative estimates of the basic properties of the circumgalactic medium (CGM), such as its density and metallicity, depend on the spectrum of incident UV background radiation. Models of UV background are known to have large variations, mainly because they are synthesized using poorly constrained parameters, which introduce uncertainty in the inferred properties of the CGM. Here, we quantify this uncertainty using a large set of new UV background models with physically motivated toy models of metal-enriched CGM. We find that, the inferred density and metallicity of low-density ($10^{-5}$ cm$^{-3}$) gas is uncertain by factors of 6.3 and 3.7, whereas high density ($10^{-3}$ cm$^{-3}$) gas by factors of 3.3 and 2.2, respectively. The variation in the shape of the UV background models is entirely responsible for such a variation in the metallicity while variation in the density arises from both normalization and shape of the UV background. Moreover, we find a harder (softer) UV background infers higher (lower) density and metallicity. We also study warm-hot gas at $T= 10^{5.5}$ K and find that metallicity is robustly estimated but the inferred density is uncertain by a factor of 3 to 5.4 for low to high-density gas. Such large uncertainties in density and metallicity may severely limit the studies of the CGM and demand better observational constraints on the input parameters used in synthesizing UV background.

Stellar evolution and numerical hydrodynamics simulations depend critically on access to fast, accurate, thermodynamically consistent equations of state. We present Skye, a new equation of state for fully-ionized matter. Skye includes the effects of positrons, relativity, electron degeneracy, Coulomb interactions, non-linear mixing effects, and quantum corrections. Skye determines the point of Coulomb crystallization in a self-consistent manner, accounting for mixing and composition effects automatically. A defining feature of this equation of state is that it uses analytic free energy terms and provides thermodynamic quantities using automatic differentiation machinery. Because of this, Skye is easily extended to include new effects by simply writing new terms in the free energy. We also introduce a novel thermodynamic extrapolation scheme for extending analytic fits to the free energy beyond the range of the fitting data while preserving desirable properties like positive entropy and sound speed. We demonstrate Skye in action in the MESA stellar evolution software instrument by computing white dwarf cooling curves.

We aim to determine the mass, velocity anisotropy, and pseudo phase-space density profiles (M(r), beta(r), and Q(r), respectively) of clusters of galaxies at the highest redshifts investigated in detail so far. We combine the GOGREEN and GCLASS spectroscopic data-sets for 14 clusters with mass M200 > 10^14 Msolar at redshifts 0.9 < z < 1.4. We stack these 14 clusters into an ensemble cluster of 581 member galaxies with stellar mass > 10^9.5 M_solar. We use the MAMPOSSt method and the inversion of the Jeans equation technique to determine M(r) and beta(r). We then combine the results of the M(r) and beta(r) analysis to determine Q(r) for the ensemble cluster. The concentration c200 of the ensemble cluster M(r) is in excellent agreement with predictions from LambdaCDM cosmological numerical simulations, and with previous determinations for clusters of similar mass and at similar redshifts, obtained from gravitational lensing and X-ray data. We see no significant difference between the total mass density and either the galaxy number density distributions or the stellar mass distribution. Star-forming galaxies are spatially significantly less concentrated than quiescent galaxies. The orbits of cluster galaxies are isotropic near the center and more radial outside. Star-forming galaxies and galaxies of low stellar mass tend to move on more radially elongated orbits than quiescent galaxies and galaxies of high stellar mass. Q(r), determined either using the total mass or the number density profile, is very close to the power-law behavior predicted by numerical simulations. The internal dynamics of clusters at the highest redshift probed in detail so far are very similar to those of lower-redshift clusters, and in excellent agreement with predictions of numerical simulations. The clusters in our sample have already reached a high degree of dynamical relaxation. (Abridged)

Boltzmann solvers are an important tool for the computation of cosmological observables in the linear regime. They involve solving the Boltzmann equation, followed by an integration in momentum space, to arrive at the desired fluid properties. This is a cumbersome, computationally expensive procedure. In this work we introduce the so-called generalized Boltzmann hierarchy (GBH) for massive neutrinos in cosmology, a simpler alternative to the usual Boltzmann hierarchy, where the momentum dependence is integrated out leaving us with a two-parameter infinite set of ordinary differential equations. Along with the usual expansion in multipoles, there is now also an expansion in higher velocity weight integrals of the distribution function. We show that the GBH produces the density contrast neutrino transfer function to a per mille level accuracy at both large and intermediate scales compared to the neutrino free-streaming scale. Furthermore, by introducing a switch to a viscous fluid approximation after horizon crossing, we show that the GBH can achieve over all scales the same accuracy as the standard CLASS approach in its default precision settings. The GBH is then a powerful tool to include neutrino anisotropies in the computation of cosmological observables in linear theory, with integration being simpler and potentially faster than standard methods.

We provide a systematic study of the position-dependent correlation function in weak lensing convergence maps and its relation to the squeezed limit of the three-point correlation function (3PCF) using state-of-the-art numerical simulations. We relate the position-dependent correlation function to its harmonic counterpart, i.e., the position-dependent power spectrum or equivalently the integrated bispectrum. We use a recently proposed improved fitting function, BiHalofit, for the bispectrum to compute the theoretical predictions as a function of source redshifts. In addition to low redshift results ($z_s=1.0-2.0$) we also provide results for maps inferred from lensing of the cosmic microwave background, i.e., $z_s=1100$. We include a {\em Euclid}-type realistic survey mask and noise. In agreement with the recent studies on the position-dependent power spectrum, we find that the results from simulations are consistent with the theoretical expectations when appropriate corrections are included.

We present stellar evolution calculations from the Asymptotic Giant Branch (AGB) to the Planetary Nebula (PN) phase for models of initial mass 1.2 M\odot and 2.0 M\odot that experience a Late Thermal Pulse (LTP), a helium shell flash that occurs following the AGB and causes a rapid looping evolution between the AGB and PN phase. We use these models to make comparisons to the central star of the Stingray Nebula, V839 Ara (SAO 244567). The central star has been observed to be rapidly evolving (heating) over the last 50 to 60 years and rapidly dimming over the past 20 - 30 years. It has been reported to belong to the youngest known planetary nebula, now rapidly fading in brightness. In this paper we show that the observed timescales, sudden dimming, and increasing Log(g), can all be explained by LTP models of a specific variety. We provide a possible explanation for the nebular ionization, the 1980s sudden mass loss episode, the sudden decline in mass loss, and the nebular recombination and fading.

We review the needs of the supernova community for improvements in survey coordination and data sharing that would significantly boost the constraints on dark energy using samples of Type Ia supernovae from the Vera C. Rubin Observatories, the \textit{Nancy Grace Roman Space Telescope}, and the \textit{Euclid} Mission. We discuss improvements to both statistical and systematic precision that the combination of observations from these experiments will enable. For example, coordination will result in improved photometric calibration, redshift measurements, as well as supernova distances. We also discuss what teams and plans should be put in place now to start preparing for these combined data sets. Specifically, we request coordinated efforts in field selection and survey operations, photometric calibration, spectroscopic follow-up, pixel-level processing, and computing. These efforts will benefit not only experiments with Type Ia supernovae, but all time-domain studies, and cosmology with multi-messenger astrophysics.

We investigate the accretion induced spin up of the black hole via numerical simulations. Our method is based on general-relativistic magneto-hydrodynamics of the slowly-rotating flows in the Kerr metric, where possibly transonic shock fronts may form. We account for the changing black hole mass and spin during accretion which enforces dynamical evolution of the space-time metric. We first study non-magnetized flows with shocks, and we also include magnetic field endowed in the gas. The aim of this study is to verify whether the high mass black holes may be produced with large spins, even though at birth the collapsars might have contained slowly, or moderately spinning cores. In this way, we put constraints on the content of angular momentum in the collapsing massive stars. Our studies are also showing that shock fronts and magnetic fields may halt accretion and limit the black hole spin-up in the exploding supernovae.

Extreme Mass Ratio Inspirals (EMRIs) can be classified as dry EMRIs and wet EMRIs based on their formation mechanisms. Dry (or the "loss-cone") EMRIs, previsouly considered as the main EMRI sources for the Laser Interferometer Space Antenna, are primarily produced by multi-body scattering in the nuclear star cluster and gravitational capture. In this Letter, we highlight an alternative EMRI formation channel: (wet) EMRI formation assisted by the accretion flow around accreting galactic-center massive black holes (MBHs). In this channel, the accretion disk captures stellar-mass black holes that are intially moving on inclined orbits, and subsequently drives them to migrate towards the MBH - this process boosts the formation rate of EMRIs in such galaxies by orders of magnitude. Taking into account the fraction ($\mathcal O(10^{-2}-10^{-1})$) of active galactic nuclei where the MBHs are expected to be rapidly accreting, we forecast that wet EMRIs will contribute an important fraction of all EMRIs observed by spaceborne gravitational wave detectors and likely dominate for MBH hosts heavier than a few $10^5 M_\odot$.

The exquisite measurements of the cosmic microwave background (CMB) fluctuations by Planck allows us to tightly constrain the amplitude of matter fluctuations at redshift $\sim 1100$ in the $\Lambda$-cold dark matter ($\Lambda$CDM) model. This amplitude can be extrapolated to the present epoch, yielding constraints on the value of the $\sigma_8$ parameter. On the other hand the abundance of Sunyaev-Zeldovich (SZ) clusters detected by Planck, with masses inferred by using hydrostatic equilibrium estimates, leads to a significantly lower value of the same parameter. This discrepancy is often dubbed the "$\sigma_8$ tension" in the literature and is sometimes regarded as a possible sign of new physics. Here, we examine a direct determination of $\sigma_8$ at the present epoch in $\Lambda$CDM, and thereby the cluster mass calibrations using cosmological data at low redshift, namely the measurements of $f\sigma_8$ from the analysis of the completed Sloan Digital Sky Survey (SDSS): we combine redshift-space distortions measurements with Planck CMB constraints, X-ray, and SZ cluster counts within the $\Lambda$CDM framework, but leaving the present day amplitude of matter fluctuations as an independent parameter (i.e. no extrapolation is made from high-redshift CMB constraints). The calibration of X-ray and SZ masses are therefore left as free parameters throughout the whole analysis. Our study yields tight constraints on the aforementioned calibrations, with values entirely consistent with results obtained from the full combination of CMB and cluster data only. Such agreement suggests an absence of tension in the $\Lambda$CDM model between CMB-based estimates of $\sigma_8$ and constraints from low-redshift on $f\sigma_8$ but indicates a tension with the standard calibration of clusters masses.

The nearly logarithmic radiative impact of CO$_2$ means that planets near the outer edge of the liquid water habitable zone (HZ) require $\sim$10$^6$x more CO$_2$ to maintain temperatures conducive to standing liquid water on the planetary surface than their counterparts near the inner edge. This logarithmic radiative response also means that atmospheric CO$_2$ changes of a given mass will have smaller temperature effects on higher pCO$_2$ planets. Ocean pH is linked to atmospheric pCO$_2$ through seawater carbonate speciation and calcium carbonate dissolution/precipitation, and the response of pH to changes in pCO$_2$ also decreases at higher initial pCO$_2$. Here, we use idealized climate and ocean chemistry models to demonstrate that CO$_2$ perturbations large enough to cause catastrophic changes to surface temperature and ocean pH on low-pCO$_2$ planets in the innermost region of the HZ are likely to have much smaller effects on planets with higher pCO$_2$. Major bouts of extraterrestrial fossil fuel combustion or volcanic CO$_2$ outgassing on high-pCO$_2$ planets in the mid-to-outer HZ should have mild or negligible impacts on surface temperature and ocean pH. Owing to low pCO$_2$, Phanerozoic Earth's surface environment may be unusually volatile compared to similar planets receiving lower instellation.

Searches for continuous gravitational waves from $\textit{unknown}$ Galactic neutron stars provide limits on the shapes of neutron stars. A rotating neutron star will produce gravitational waves if asymmetric deformations exist in its structure that are characterized by the star's ellipticity. In this study, we use a simple model of the spatial and spin distribution of Galactic neutron stars to estimate the total number of neutron stars probed, using gravitational waves, to a given upper limit on the ellipticity. This may help optimize future searches with improved sensitivity. The improved sensitivity of third-generation gravitational wave detectors may increase the number of neutron stars probed, to a given ellipticity, by factors of 100 to 1000.

Liger is a next-generation near-infrared imager and integral field spectrograph (IFS) for the W.M. Keck Observatory designed to take advantage of the Keck All-Sky Precision Adaptive Optics (KAPA) upgrade. Liger will operate at spectral resolving powers between R$\sim$4,000 - 10,000 over a wavelength range of 0.8-2.4$\mu$m. Liger takes advantage of a sequential imager and spectrograph design that allows for simultaneous observations between the two channels using the same filter wheel and cold pupil stop. We present the design for the filter wheels and pupil mask and their location and tolerances in the optical design. The filter mechanism is a multi-wheel design drawing from the heritage of the current Keck/OSIRIS imager single wheel design. The Liger multi-wheel configuration is designed to allow future upgrades to the number and range of filters throughout the life of the instrument. The pupil mechanism is designed to be similarly upgradeable with the option to add multiple pupil mask options. A smaller wheel mechanism allows the user to select the desired pupil mask with open slots being designed in for future upgrade capabilities. An ideal pupil would match the shape of the image formed of the primary and would track its rotation. For different pupil shapes without tracking we model the additional exposure time needed to achieve the same signal to noise of an ideal pupil and determine that a set of fixed masks of different shapes provides a mechanically simpler system with little compromise in performance.

Liger is a next generation adaptive optics (AO) fed integral field spectrograph (IFS) and imager for the W. M. Keck Observatory. This new instrument is being designed to take advantage of the upgraded AO system provided by Keck All-Sky Precision Adaptive-optics (KAPA). Liger will provide higher spectral resolving power (R$\sim$4,000-10,000), wider wavelength coverage ($\sim$0.8-2.4 $\mu$m), and larger fields of view than any current IFS. We present the design and analysis for a custom-made dewar chamber for characterizing the Liger opto-mechanical system. This dewar chamber is designed to test and assemble the Liger imaging camera and slicer IFS components while being adaptable for future experiments. The vacuum chamber will operate below $10^{-5}$ Torr with a cold shield that will be kept below 90 K. The dewar test chamber will be mounted to an optical vibration isolation platform and further isolated from the cryogenic and vacuum systems with bellows. The cold head and vacuums will be mounted to a custom cart that will also house the electronics and computer that interface with the experiment. This test chamber will provide an efficient means of calibrating and characterizing the Liger instrument and performing future experiments.

Non-potential magnetic energy promptly released in solar flares is converted to other forms of energy. This may include nonthermal energy of flare-accelerated particles, thermal energy of heated flaring plasma, and kinetic energy of eruptions, jets, up/down flows, and stochastic (turbulent) plasma motions. The processes or parameters governing partitioning of the released energy between these components is an open question. How these components are distributed between distinct flaring loops and what controls these spatial distributions is also unclear. Here, based on multi-wavelength data and 3D modeling, we quantify the energy partitioning and spatial distribution in the well observed SOL2014-02-16T064620 solar flare of class C1.5. Nonthermal emissions of this flare displayed a simple impulsive single-spike light curves lasting about 20\,s. In contrast, the thermal emission demonstrated at least three distinct heating episodes, only one of which was associated with the nonthermal component. The flare was accompanied by up and down flows and substantial turbulent velocities. The results of our analysis suggest that (i) the flare occurs in a multi-loop system that included at least three distinct flux tubes; (ii) the released magnetic energy is divided unevenly between the thermal and nonthermal components in these loops; (iii) only one of these three flaring loops contains an energetically important amount of nonthermal electrons, while two other loops remain thermal; (iv) the amounts of direct plasma heating and that due to nonthermal electron loss are comparable; (v) the kinetic energy in the flare footpoints constitute only a minor fraction compared with the thermal and nonthermal energies.

In this work, we calibrate the relationship between Halpha emission and M dwarf ages. We compile a sample of 892 M dwarfs with Halpha equivalent width (HaEW) measurements from the literature that are either co-moving with a white dwarf of known age (21 stars) or in a known young association (871 stars). In this sample we identify 7 M dwarfs that are new candidate members of known associations. By dividing the stars into active and inactive categories according to their HaEW and spectral type (SpT), we find that the fraction of active dwarfs decreases with increasing age, and the form of the decline depends on SpT. Using the compiled sample of age-calibrators we find that HaEW and fractional Halpha luminosity (LHaLbol) decrease with increasing age. HaEW for SpT<M7 decreases gradually up until ~1Gyr. For older ages, we found only two early M dwarfs which are both inactive and seem to continue the gradual decrease. We also found 14 mid-type out of which 11 are inactive and present a significant decrease of HaEW, suggesting that the magnetic activity decreases rapidly after ~1Gyr. We fit LHaLbol versus age with a broken power-law and find an index of -0.11+0.02-0.01 for ages <~776Myr. The index becomes much steeper at older ages however a lack of field age-calibrators leaves this part of the relation far less constrained. Finally, from repeated independent measurements for the same stars we find that 94% of these has a level of HaEW variability <=5A at young ages (<1Gyr).

We report results from continued timing observations of PSR J0740+6620, a high-mass, 2.8-ms radio pulsar in orbit with a likely ultra-cool white dwarf companion. Our data set consists of combined pulse arrival-time measurements made with the 100-m Green Bank Telescope and the Canadian Hydrogen Intensity Mapping Experiment telescope. We explore the significance of timing-based phenomena arising from general-relativistic dynamics and variations in pulse dispersion. When using various statistical methods, we find that combining $\sim 1.5$ years of additional, high-cadence timing data with previous measurements confirms and improves upon previous estimates of relativistic effects within the PSR J0740+6620 system, with the pulsar mass $m_{\rm p} = 2.08^{+0.07}_{-0.07}$ M$_\odot$ (68.3\% credibility) determined by the relativistic Shapiro time delay. For the first time, we measure secular variation in the orbital period and argue that this effect arises from apparent acceleration due to significant transverse motion. After incorporating contributions from Galactic differential rotation and off-plane acceleration in the Galactic potential, we obtain a model-dependent distance of $d = 1.14^{+0.17}_{-0.15}$ kpc (68.3\% credibility). This improved distance confirms the ultra-cool nature of the white dwarf companion determined from recent optical observations. We discuss the prospects for future observations with next-generation facilities, which will likely improve the precision on $m_{\rm p}$ for J0740+6620 by an order of magnitude within the next few years.

In nominal mission scenarios, geostationary satellites perform end-of-life orbit maneuvers to reach suitable disposal orbits, where they do not interfere with operational satellites. This research investigates the long-term orbit evolution of decommissioned geostationary satellite under the assumption that the disposal maneuver does not occur and the orbit evolves with no control. The dynamical model accounts for all the relevant harmonics of the gravity field at the altitude of geostationary orbits, as well as solar radiation pressure and third-body perturbations caused by the Moon and the Sun. Orbit propagations are performed using two algorithms based on different equations of motion and numerical integration methods: (i) Gauss planetary equations for modified equinoctial elements with a Runge-Kutta numerical integration scheme based on 8-7th-order Dorman and Prince formulas; (ii) Cartesian state equations of motion in an Earth-fixed frame with a Runge-Kutta Fehlberg 7/8 integration scheme. The numerical results exhibit excellent agreement over integration times of decades. Some well-known phenomena emerge, such as the longitudinal drift due to the resonance between the orbital motion and Earth's rotation, attributable to the J22 term of the geopotential. In addition, the third-body perturbation due to Sun and Moon causes two major effects: (a) a precession of the orbital plane, and (b) complex longitudinal dynamics. This study proposes an analytical approach for the prediction of the precessional motion and shows its agreement with the orbit evolution obtained numerically. Moreover, long-term orbit propagations show that the above mentioned complex longitudinal dynamics persists over time scales of several decades. Frequent and unpredictable migrations toward different longitude regions occur, in contrast with the known effects due only to the J22 perturbation.

At relatively high frequencies, highly sensitive grating sidelobes occur in the primary beam patterns of low frequency aperture arrays (LFAA) such as the Murchison Widefield Array (MWA). This occurs when the observing wavelength becomes comparable to the dipole separation for LFAA tiles, which for the MWA occurs at approximately 300 MHz. The presence of these grating sidelobes has made calibration and image processing for 300 MHz MWA observations difficult. This work presents a new calibration and imaging strategy which employs existing techniques to process two example 300 MHz MWA observations. Observations are initially calibrated using a new 300 MHz sky-model which has been interpolated from low frequency and high frequency all-sky surveys. Using this 300 MHz model in conjunction with the accurate MWA tile primary beam model, we perform sky-model calibration for the two example observations. After initial calibration a self-calibration loop is performed by all-sky imaging each observation with WSCLEAN. Using the output all-sky image we mask the main lobe of the image. Using this masked image we perform a sky-subtraction by estimating the masked image visibilities using WSCLEAN. We then image the main lobe of the observations with WSCLEAN. This results in high dynamic range images of the two example observation main lobes. These images have a resolution of 2.4 arcminutes, with a maximum sensitivity of 31 mJy/beam. The calibration and imaging strategy demonstrated in this work, opens the door to performing science at 300 MHz with the MWA, which was previously an inaccessible domain. With this paper we release the code described below and the cross-matched catalogue along with the code to produce a sky-model in the range 70-1400 MHz.

We study waiting times of Gamma-ray flares of Flat Spectrum Radio Quasars. We adopt here a Scan-Statistic driven clustering method iSRS, (Pacciani 2018) to recognize flaring states within the FERMI-LAT archival data. We found that flares waiting times can be described with overlapping bursts of flares, with an average burst duration of ~0.6 year, and with an average burst rate of ~1.3/y. We also found evidence for a second population of waiting times for times shorter than several days (blue component). We found a few tens of cases of short waiting times. In our analysis, CTA 102 showed the large majority of short waiting times. Interestingly, the period of conspicuous detection of the blue component of waiting times for CTA 102 coincides with the crossing time of the K1 feature with the C1 stationary feature in radio reported in (Jorstad 2017,Casadio 2019). We interpret the obtained distribution of waiting times of Gamma-ray flares as originating from cold plasma moving along the jet for a deprojected length of ~30-60 pc (assuming a bulk Gamma=10), that sporadically produce Gamma-ray flares. Duration and Burst rate is in agreement with distribution of fading time and ejection rate of traveling structures observed with VLBA at 43 GHz.

A novel application of machine-learning (ML) based image processing algorithms is proposed to analyze an all-sky map (ASM) obtained using the Fermi Gamma-ray Space Telescope. An attempt was made to simulate a one-year ASM from a short-exposure ASM generated from one-week observation by applying three ML based image processing algorithms: dictionary learning, U-net, and Noise2Noise. Although the inference based on ML is less clear compared to standard likelihood analysis, the quality of the ASM was generally improved. In particular, the complicated diffuse emission associated with the galactic plane was successfully reproduced only from one-week observation data to mimic a ground truth (GT) generated from a one-year observation. Such ML algorithms can be implemented relatively easily to provide sharper images without various assumptions of emission models. In contrast, large deviations between simulated ML maps and GT map were found, which are attributed to the significant temporal variability of blazar-type active galactic nuclei (AGNs) over a year. Thus, the proposed ML methods are viable not only to improve the image quality of an ASM, but also to detect variable sources, such as AGNs, algorithmically, i.e., without human bias. Moreover, we argue that this approach is widely applicable to ASMs obtained by various other missions; thus, it has the potential to examine giant structures and transient events, both of which are rarely found in pointing observations.

We present the development of a machine learning based pipeline to fully automate the calibration of the frequency comb used to read out optical/IR Microwave Kinetic Inductance Detector (MKID) arrays. This process involves determining the resonant frequency and optimal drive power of every pixel (i.e. resonator) in the array, which is typically done manually. Modern optical/IR MKID arrays, such as DARKNESS (DARK-speckle Near-infrared Energy-resolving Superconducting Spectrophotometer) and MEC (MKID Exoplanet Camera), contain 10-20,000 pixels, making the calibration process extremely time consuming; each 2000 pixel feedline requires 4-6 hours of manual tuning. Here we present a pipeline which uses a single convolutional neural network (CNN) to perform both resonator identification and tuning simultaneously. We find that our pipeline has performance equal to that of the manual tuning process, and requires just twelve minutes of computational time per feedline.

We here report a monitor of the BL Lac object 1ES 1218+304 in both B- and R-bands by the GWAC-F60A telescope in eight nights, when it was triggerd to be at its highest X-ray flux in history by the VERITAS Observatory and Swift follow-ups. Both ANOVA and $\chi^2$-test enable us to clearly reveal an intra-day variability in optical wavelengths in seven out of the eight nights. A bluer-when-brighter chromatic relationship has been clearly identified in five out of the eight nights, which can be well explained by the shock-in-jet model. In addtion, a quasi-periodic oscilation phenomenon in both bands could be tentatively identified in the first night. A positive delay between the two bands has been revealed in three out of the eight nights, and a negative one in the other nights. The identfied minimum time delay enables us to estimate the $M_{\mathrm{BH}}=2.8\times10^7 \rm M_{\odot}$that is invalid.

Narrowband spikes are observed in solar flares for several decades. However, their exact origin is still discussed. To contribute to understanding of these spikes, we analyze the narrowband spikes observed in the 800-2000 MHz range during the impulsive phase of the November 7, 2013 flare. In the radio spectrum, the spikes started with typical broadband clouds of spikes, and then their distribution in frequencies changed into unique, very narrow bands having non-integer frequency ratios. We successfully fitted frequencies of these narrow spike bands by those, calculating dispersion branches and growth rates of the Bernstein modes. For comparison, we also analyzed the model, where the narrow bands of spikes are generated at the upper-hybrid frequencies. Using both models, we estimated the plasma density and magnetic field in spike sources. Then the models are discussed, and arguments in favor of the model with the Bernstein modes are presented. Analyzing frequency profiles of this spike event by the Fourier method, we found the power-law spectra with the power-law indices varying in the -0.8 -- -2.75 interval. Because at some times this power-law index was close to the Kolmogorov spectral index (-5/3), we propose that the spikes are generated through the Bernstein modes in turbulent plasma reconnection outflows or directly in the turbulent magnetic reconnection of solar flares.

Gravitational microlensing has become a mature technique for discovering small gravitational lenses in the Universe which are otherwise beyond our detection limits. Similarly, plasma microlensing can help us explore cosmic plasma lenses. Both pulsar scintillation and extreme scattering events of compact radio sources suggest the existence of ~ AU scale plasma lenses in the ionized intersteller medium (IISM), whose astrophysical correspondence remains a mystery. We propose that plasma microlensing events by these plasma lenses recorded in the form of wide-band dynamic spectra are a powerful probe of their nature. Using the recently developed Picard-Lefschetz integrator for the Kirchhoff-Fresnel integral, we simulate such dynamic spectra for a well-motivated family of single-variable plasma lenses. We demonstrate that the size, strength, and shape of the plasma lens can be measured from the location of the cusp point and the shape of spectral caustics respectively, with a combination of distances and the effective velocity known a priori or measured from the widths of the interference pattern. Future wide-band observations of pulsars, whose plasma microlensing events may be predictable from parabolic arc monitoring, are the most promising ground to apply our results for a deeper insight into the fine structures in the IISM.

We offer an exact analytical equation for the void central region. We show that the central density is solely determined by the amplitude of the initial perturbation. Our results suggest that N-body simulations somewhat overestimate the emptiness of voids: the majority of them should have the central underdensity $\delta_c > -73\%$, and there should be almost no voids with $\delta_c < -88\%$. The central region of a void is a part of an open Friedmann's 'universe', and its evolution differs drastically from the Universe evolution: there is a long stage when the curvature term dominates, which prevents the formation of galaxy clusters and massive galaxies inside voids. The density profile in the void center should be very flat. We discuss some void models obtained by N-body simulations and offer some ways to improve them. We also show that the dark energy makes the voids less underdense.

We carry out a test of the cosmic distance duality relation using a sample of 58 SPT-SZ clusters, along with X-ray measurements from XMM-Newton. To carry out this test, we need an estimate of the luminosity distance ($D_L$) at the redshift of the cluster. We use three independent methods for this purpose: directly using $D_L$ from the closest Type Ia Supernovae from the Union 2.1 sample, non-parametric reconstruction of $D_L$ using the same Union 2.1 sample, and finally using $H(z)$ measurements from cosmic chronometers and reconstructing $D_L$ using Gaussian Process regression. We use four different functions to characterize the deviations from CDDR. All our results for these ($4 \times 3$) analyses are consistent with CDDR to within $1-3\sigma$.

We investigated the dependence of the parameters of the segments of spiral arms of the Galaxy on the age of classical Cepheids. The catalog of Cepheids (Mel'nik et al. 2015) was divided into two samples$-$relatively young ($P>9^\text{d}$) and relatively old ($P<9^\text{d}$) objects. The parameters of the spiral structure were determined both for two samples separately and jointly for the combination of two systems of segments traced by young and old objects. For most of the segments, their parameters for young and old objects differ significantly. Taking into account the difference between the two segment systems, we obtained the estimate $R_0$ equal to $7.23^{+0.19}_{-0.18}$ kpc, which in the modern LMC calibration corresponds to the value of $R_0={8.08^{+0.21}_{-0.20}}|_{\text{stat.}} {}^{+0.38}_{-0.36}|_{\text{cal.}}$ kpc. It is shown that the displacement between the segments is not reduced to the effect of differential rotation only. To interpret this displacement for objects of Perseus and Sagittarius-2 segments we carried out a dynamic modeling of the change in the position of the segment points when moving in the smooth gravitational field of the Galaxy. At the angular velocity of rotation of the spiral pattern $\Omega_{\text{p}} = 25.2 \pm 0.5$ km s$^{-1}$ kpc$^{-1}$ (Dambis et al. 2015) the observed displacement between segments on young and old objects can be explained by the amplitude of spiral perturbations of the radial velocity of $u = 10 \pm1.5$ km s$^{-1}$. For the constructed double system of spiral segments, it is demonstrated that the assumption of constancy of the pitch angles within each segment and the assumption that the pole of the spiral pattern is in the direction of the nominal center of the Galaxy do not contradict the data within the range of uncertainty.

We use linear perturbation theory to study perturbations in dynamical dark energy models. We compare quintessence and tachyonic dark energy models with identical background evolution. We write the corresponding equations for different models in a form that makes it easier to see that the two models are very hard to distinguish in the linear regime, especially for models with $(1 + w) \ll 1$. We use Cosmic Microwave Background data and parametric representations for the two models to illustrate that they cannot be distinguished for the same background evolution with existing observations. Further, we constrain tachyonic models with the Planck data. We do this analysis for exponential and inverse square potentials and find that the intrinsic parameters of the potentials remain very weakly constrained. In particular, this is true in the regime allowed by low redshift observations.

One of the primary sources of uncertainties in modeling the cosmic-shear power spectrum on small scales is the effect of baryonic physics. Accurate cosmology for Stage-IV surveys requires knowledge of the matter power spectrum deep in the nonlinear regime at the percent level. Therefore, it is important to develop reliable mitigation techniques to take into account baryonic uncertainties if information from small scales is to be considered in the cosmological analysis. In this work, we develop a new mitigation method for dealing with baryonic physics for the case of the shear angular power spectrum. The method is based on an extended covariance matrix that incorporates baryonic uncertainties informed by hydrodynamical simulations. We use the results from 13 hydrodynamical simulations and the residual errors arising from a fit to a $\Lambda$CDM model using the extended halo model code {\tt HMCode} to account for baryonic physics. These residual errors are used to model a so-called theoretical error covariance matrix that is added to the original covariance matrix. In order to assess the performance of the method, we use the 2D tomographic shear from four hydrodynamical simulations that have different extremes of baryonic parameters as mock data and run a likelihood analysis comparing the residual bias on $\Omega_m$ and $\sigma_8$ of our method and the HMCode for an LSST-like survey. We use different modelling of the theoretical error covariance matrix to test the robustness of the method. We show that it is possible to reduce the bias in the determination of the tested cosmological parameters at the price of a modest decrease in the precision.

Over the last two decades, a large population of close-in planets has been detected around a wide variety of host stars. Such exoplanets are likely to undergo planetary migration through magnetic and tidal interactions. We aim to follow the orbital evolution of a planet along the structural and rotational evolution of its host star by taking into account simultaneously tidal and magnetic torques, in order to explain some properties of the distribution of observed close-in planets. We rely on a numerical model of a coplanar circular star-planet system taking into account stellar structural changes, wind braking and star-planet interactions, called ESPEM. We browse the parameter space of star-planet systems' configurations and assess the relative influence of magnetic and tidal torques on secular evolution. We then synthesize star-planet populations and confront their distribution in orbital and stellar rotation periods to Kepler satellite data. First, we find that after the dissipation of the protoplanetary disk, both types of interactions can dominate secular evolution depending on the initial configuration of the system and the evolutionary phase considered. Moreover, different populations of star-planet systems emerge from the combined action of both kinds of interactions, according to the evolutionary phase during which the planet migrates significantly. This may affect significantly the detectability of star-planet systems as well as the validity of gyrochonology. All in all, star-planet interactions significantly impact the global distribution in orbital periods during the main sequence, while the global distribution in stellar rotation periods is marginally affected. More precisely, star-planet magnetic interactions significantly affect the distribution of super-Earths around slowly rotating stars, while tidal effects are found to shape the distribution of giant planets.

Oxygenic photosynthesis is the most important biochemical process in Earth biosphere and likely very common on other habitable terrestrial planets, given the general availability of its input chemical ingredients and of light as source of energy. It is therefore important to evaluate the effective possibility of oxygenic photosynthesis on planets around stars as a function of their spectral type and the planet-star separation. We aim at estimating the photon flux, the exergy and the exergetic efficiency of the radiation in the wavelength range useful for the oxygenic photosynthesis as a function of the host star effective temperature and planet-star separation. We compute analytically these quantities and compare our results with the estimates for the small sample of known Earth-like planets. We find that exergy is an increasing function of the star effective temperature, within the range 2600-7200 K. It depends both on the star-planet separation and the star effective temperature. Biospheres on exoplanets around cool stars might be generally light-limited. So far, we have not observed terrestrial planets comparable to Earth in terms of useful photon flux, exergy and exergetic efficiency.

We consider a machine learning algorithm to detect and identify strong gravitational lenses on sky images. First, we simulate different artificial but very close to reality images of galaxies, stars and strong lenses, using six different methods, i.e. two for each class. Then we deploy a convolutional neural network architecture to classify these simulated images. We show that after neural network training process one achieves about 93 percent accuracy. As a simple test for the efficiency of the convolutional neural network, we apply it on an real Einstein cross image. Deployed neural network classifies it as gravitational lens, thus opening a way for variety of lens search applications of the deployed machine learning scheme.

Fragmentation often occurs in disk-like structures, both in the early Universe and in the context of present-day star formation. Supermassive black holes (SMBHs) are astrophysical objects whose origin is not well understood; they weigh millions of solar masses and reside in the centers of galaxies. An important formation scenario for SMBHs is based on collisions and mergers of stars in a massive cluster, in which the most massive star moves to the center of the cluster due to dynamical friction. This increases the rate of collisions and mergers since massive stars have larger collisional cross sections. This can lead to runaway growth of a very massive star which may collapse to become an intermediate-mass black hole. Here we investigate the dynamical evolution of Miyamoto-Nagai models that allow us to describe dense stellar clusters, including flattening and different degrees of rotation. We find that the collisions in these clusters depend mostly on the number of stars and the initial stellar radii for a given radial size of the cluster. By comparison, rotation seems to affect the collision rate by at most $20\%$. For flatness, we compared spherical models with systems that have a scale height of about $10\%$ of their radial extent, in this case finding a change in the collision rate of less than $25\%$. Overall, we conclude that the parameters only have a minor effect on the number of collisions. Our results also suggest that rotation helps to retain more stars in the system, reducing the number of escapers by a factor of $2-3$ depending on the model and the specific realization. After two million years, a typical lifetime of a very massive star, we find that about $630$ collisions occur in typical models with $N=10^4$, $R=100$ $\rm~R_\odot$ and a half-mass radius of $0.1$ $\rm~pc$, leading to a mass of about $6.3\times10^3$ $\rm~M_\odot$ for the most massive object.

For centuries, the rich nocturnal environment of the starry sky could be modelled only by analogue tools such as paper planispheres, atlases, globes and numerical tables. The immersive sky simulator of the twentieth century, the optomechanical planetarium, provided new ways for representing and teaching about the sky, but the high construction and running costs meant that they have not become common. However, in recent decades, "desktop planetarium programs" running on personal computers have gained wide attention. Modern incarnations are immensely versatile tools, mostly targeted towards the community of amateur astronomers and for knowledge transfer in transdisciplinary research. Cultural astronomers also value the possibilities they give of simulating the skies of past times or other cultures. With this paper, we provide an extended presentation of the open-source project Stellarium, which in the last few years has been enriched with capabilities for cultural astronomy research not found in similar, commercial alternatives.

In the majority of population synthesis calculations of close binary stars, the common envelope (CE) phase is modeled using a standard prescription based upon the conservation of energy, known as the alpha prescription. In this prescription, the orbital separation of the secondary and giant core at the end of the CE phase is taken to be the orbital separation when the envelope becomes unbound. However, recent observations of planetary nebulae with binary cores (BPNe), believed to be the immediate products of CE evolution, indicate orbital periods that are significantly shorter than predicted by population synthesis models using the alpha prescription. We argue that unbinding the envelope provides a necessary, but not sufficient, condition to escape a merger during CE evolution. The spiral-in of the secondary must also be halted. This requires the additional dynamical constraint that the frictional torque on the secondary be reduced to approximately zero. In this paper, we undertake a preliminary examination of the effect of adding this dynamical constraint in population synthesis calculations of BPNe. We assume that the frictional torque will be sufficiently reduced when the secondary enters a region within the giant where the mass-radius profile is flat. We crudely estimate the location of this region as a function of the core mass based upon existing stellar models of AGB stars between 1 and 7 solar masses. We calculate a theoretical orbital period distribution of BPNe using a population synthesis code that incorporates this dynamical constraint along with the alpha prescription.

Gas is now discovered commonly in exoplanetary systems with planetesimal belts. In this paper, we verify whether gas is also expected in the Kuiper belt (KB) in our Solar System. We predict that gas is still produced in the KB right now at a rate of $2 \times 10^{-8}$ M$_\oplus$/Myr for CO. Once released, the gas is quickly pushed out by the Solar wind. Therefore, we predict a gas wind in our Solar System starting at the KB location and extending far beyond with a current total CO mass of $\sim 2 \times 10^{-12}$ M$_\oplus$ (i.e. 20 times the CO quantity that was lost by the Hale-Bopp comet during its sole 1997 passage) and CO density in the belt of $3 \times 10^{-7}$ cm$^{-3}$. We also predict the existence of a slightly more massive atomic gas wind made of carbon and oxygen (neutral and ionized) with a mass of $\sim 10^{-11}$ M$_\oplus$. We predict that this gas could be observed with future in-situ missions and may play an important role for, e.g., planetary atmosphere formation in the Solar System history.

We simultaneously and successfully fit the multi-epoch X-ray spectra of the tidal disruption event (TDE) 3XMM J215022.4-055108 using a modified version of our relativistic slim disk model that now accounts for angular momentum losses from radiation. We explore the effects of different disk properties and of uncertainties in the spectral hardening factor fc and redshift z on the estimation of the black hole mass M and spin a. Across all choices of theoretical priors, we constrain M to less than 2.2e4 Ms at 1 sigma confidence. Assuming that the TDE host is a star cluster associated with the adjacent, giant, barred lenticular galaxy at z=0.055, we constrain M and a to be (1.75+0.45-0.05)e4 Ms and 0.8+0.12-0.02, respectively, at 1 sigma confidence. The high, but sub-extremal, spin suggests that, if this intermediate mass black hole (IMBH) has grown significantly since formation, it has acquired its last e-fold in mass in a way incompatible with both the standard and chaotic limits of gas accretion. Ours is the first clear IMBH with a spin measurement. As such, this object represents a novel laboratory for astro-particle physics; its M and a place tight limits on the existence of ultralight bosons, ruling out those with masses 1.0e-15 to 1.0e-16 eV.

Stellar population studies provide unique clues to constrain galaxy formation models. So far, detailed studies based on absorption line strengths have mainly focused on the optical spectral range although many diagnostic features are present in other spectral windows. In particular, the near-infrared (NIR) can provide a wealth of information about stars, such as evolved giants, that have less evident optical signatures. Due to significant advances in NIR instrumentation and extension of spectral libraries and stellar population synthesis (SPS) models to this domain, it is now possible to perform in-depth studies of spectral features in the NIR to a high level of precision. In the present work, taking advantage of state-of-the-art SPS models covering the NIR spectral range, we introduce a new set of NIR indices constructed to be maximally sensitive to the main stellar population parameters, namely age, metallicity and initial mass function (IMF). We fully characterize the new indices against these parameters as well as their sensitivity to individual elemental abundance variations, velocity dispersion broadening, wavelength shifts, signal-to-noise ratio and flux calibration. We also present, for the first time, a method to ensure that the analysis of spectral indices is not affected by sky contamination, which is a major challenge when dealing with NIR spectroscopy. Moreover, we discuss two main applications: (i) the ability of some NIR spectral indices to constrain the shape of the low-mass IMF and (ii) current issues in the analysis of NIR spectral indices for future developments of SPS modelling.

In this paper, the equation of state (EOS) of deconfined quark stars is studied in the framework of the two-flavor NJL model, and the self-consistent mean field approximation is employed by introducing a parameter $\alpha$ combining the original Lagrangian and the Fierz-transformed Lagrangian, $\mathcal{L}_R= (1-\alpha)\mathcal{L}+\alpha\mathcal{L}_F$, to measure the weights of different interaction channels. It is believed that the deconfinement of phase transition happens along with the chiral phase transition. Thus, due to the lack of description of confinement in the NJL model, the vacuum pressure is set to confine quarks at low densities, which is the pressure corresponding to the critical point of chiral phase transition. We find that deconfined quark stars can reach over two-solar-mass, and the bag constant therefore shifts from $(130 ~\mathrm{MeV})^4$ to $(150 ~\mathrm{MeV})^4$ as $\alpha$ grows. In addition, the tidal deformability $\Lambda$ is yielded ranging from 253 to 482 along with the decrease of $~\alpha$, which satisfies the astronomical constraint of $\Lambda<800$ for 1.4-solar-mass neutron stars.

The membership determination for open clusters in noisy environments of the Milky Way is still an open problem. In this paper, our main aim is provide the membership probability of stars using proper motions and parallax values of stars using Gaia EDR3 astrometry. Apart from the Gaia astrometry, we have also used other photometric data sets like UKIDSS, WISE, APASS and Pan-STARRS1 in order to understand cluster properties from optical to mid-infrared regions. We selected 438 likely members with membership probability higher than $50\%$ and G$\le$20 mag. We obtained the mean value of proper motion as $\mu_{x}=1.27\pm0.001$ and $\mu_{y}=-0.73\pm0.002$ mas yr$^{-1}$. The cluster's radius is determined as 7.5 arcmin (5.67 pc) using radial density profile. Our analysis suggests that NGC 1348 is located at a distance of $2.6\pm0.05$ kpc. The mass function slope is found to be $1.30\pm0.18$ in the mass range 1.0$-$4.1 $M_\odot$, which is in fair agreement with Salpeter's value within the 1$\sigma$ uncertainty. The present study validates that NGC 1348 is a dynamically relaxed cluster. We computed the apex coordinates $(A, D)$ for NGC 1348 as $(A_\circ, D_\circ)$ = $(-23^{\textrm{o}}.815\pm 0^{\textrm{o}}.135$, $-22^{\textrm{o}}.228\pm 0^{\textrm{o}}.105)$. In addition, calculations of the velocity ellipsoid parameters (VEPs), matrix elements $\mu_{ij}$, direction cosines ($l_j$, $m_j$, $n_j$) and the Galactic longitude of the vertex have been also conducted in this analysis.

AstroSat is India's first dedicated multi-wavelength space observatory launched by the Indian Space Research Organisation (ISRO) on 28 September 2015. After launch, the AstroSat Science Support Cell (ASSC) was set up as a joint venture of ISRO and the Inter-University Centre for Astronomy and Astrophysics (IUCAA) with the primary purpose of facilitating the use of AstroSat, both for making observing proposals and for utilising archival data. The ASSC organises meetings, workshops and webinars to train users in these activities, runs a help desk to address user queries, provides utility tools and disseminates analysis software through a consolidated web portal. It also maintains the AstroSat Proposal Processing System (APPS) which is deployed at ISSDC, a software platform central to the workflow management of AstroSat operations. This paper illustrates the various aspects of ASSC functionality.

Infrared quasi-stellar objects (IR QSOs) are a rare subpopulation selected from ultraluminous infrared galaxies (ULIRGs) and have been regarded as promising candidates of ULIRG-to-optical QSO transition objects. Here we present NOEMA observations of the CO(1-0) line and 3 mm continuum emission in an IR QSO IRAS F07599+6508 at $z=0.1486$, which has many properties in common with Mrk 231. The CO emission is found to be resolved with a major axis of $\sim$6.1 kpc that is larger than the size of $\sim$4.0 kpc derived for 3 mm continuum. We identify two faint CO features located at a projected distance of $\sim$11.4 and 19.1 kpc from the galaxy nucleus, respectively, both of which are found to have counterparts in the optical and radio bands and may have a merger origin. A systematic velocity gradient is found in the CO main component, suggesting that the bulk of molecular gas is likely rotationally supported. Based on the radio-to-millimeter spectral energy distribution and IR data, we estimate that about 30$\%$ of the flux at 3 mm arises from free-free emission and infer a free-free-derived star formation rate of 77 $M_\odot\ {\rm yr^{-1}}$, close to the IR estimate corrected for the AGN contribution. We find a high-velocity CO emission feature at the velocity range of about -1300 to -2000 km s$^{-1}$. Additional deep CO observations are needed to confirm the presence of a possible very high-velocity CO extension of the OH outflow in this IR QSO.

In this paper, we use large scale structure observations to test the cosmic distance duality relation (CDDR), $D_{\rm L}(1+z)^{-2}/D_{\rm A}=\eta=1 $, with $D_{\rm L}$ and $D_{\rm A}$, being the luminosity and angular diameter distances, respectively. In order to perform the test, the following data set are considered: strong lensing systems and galaxy cluster measurements (gas mass fractions). No specific cosmological model is adopted, only a flat universe is assumed. By considering two $\eta(z)$ parametrizations, we obtained the validity of the CDDR within $1.5\sigma$ which is in full agreement with other recent tests involving cosmological data. It is worth to comment that our results are independent of the baryon budget of galaxy clusters.

Using Planck polarization data, we search for and constrain spatial variations of the polarized dust foreground for cosmic microwave background (CMB) observations, specifically in its spectral index, $\beta_d$. Failure to account for such variations will cause errors in the foreground cleaning that propagate into errors on cosmological parameter recovery from the cleaned CMB map. It is unclear how robust prior studies of the Planck data which constrained $\beta_d$ variations are due to challenges with noise modeling, residual systematics, and priors. To clarify constraints on $\beta_d$ and its variation, we employ two pixel space analyses of the polarized dust foreground at $>3.7^{\circ}$ scales on $\approx 60\%$ of the sky at high Galactic latitudes. A template fitting method, which measures $\beta_d$ over three regions of $\approx 20\%$ of the sky, does not find significant deviations from an uniform $\beta_d = 1.55$, consistent with prior Planck determinations. An additional analysis in these regions, based on multifrequency fits to a dust and CMB model per pixel, puts limits on $\sigma_{\beta_d}$, the Gaussian spatial variation in $\beta_d$. At the highest latitudes, the data support $\sigma_{\beta_d}$ up to $0.45$, $0.30$ at mid-latitudes, and $0.15$ at low-latitudes. We also demonstrate that care must be taken when interpreting the current Planck constraints, $\beta_d$ maps, and noise simulations. Due to residual systematics and low dust signal to noise at high latitudes, forecasts for ongoing and future missions should include the possibility of large values of $\sigma_{\beta_d}$ as estimated in this paper, based on current polarization data.

Context. Centaurs go around the Sun between the orbits of Jupiter and Neptune. Only a fraction of the known centaurs has been found to display comet-like features. Comet 29P/Schwassmann-Wachmann 1 is the most remarkable active centaur. It orbits the Sun just beyond Jupiter in a nearly circular path. Only a handful of known objects follow similar trajectories. Aims. We present photometric observations of 2020 MK4, a recently found centaur with an orbit not too different from that of 29P, and perform a preliminary exploration of its dynamical evolution. Methods. We analyzed broadband Cousins R and Sloan g', r', and i' images of 2020 MK4 acquired with the Jacobus Kapteyn Telescope and the IAC80 telescope to search for cometary-like activity, and to derive its surface colors and size. Its orbital evolution was studied using direct N-body simulations. Results. Centaur 2020 MK4 is neutral-gray in color and has a faint, compact cometary-like coma. The values of its color indexes, (g'-r')=0.42+/-0.04 and (r'-i')=0.17+/-0.04, are similar to the solar ones. A lower limit for the absolute magnitude of the nucleus is Hg=11.30+/-0.03 mag that, for an albedo in the range 0.1-0.04, gives an upper limit for its size in the interval (23, 37) km. Its orbital evolution is very chaotic and 2020 MK4 may be ejected from the solar system during the next 200 kyr. Comet 29P experienced relatively close flybys with 2020 MK4 in the past, sometimes when they were temporary Jovian satellites. Conclusions. We confirm the presence of a coma of material around a central nucleus. Its surface colors place this centaur among the most extreme members of the gray group. Although its past, present, and future dynamical evolution resembles that of 29P, more data are required to confirm or reject a possible connection between the two objects and perhaps others.

In this work, two different models, one with cosmological constant $\Lambda$, and baryonic and dark matter (with $\omega_{dm} \neq 0$), and the other with an $X$ dark energy (with $\omega_{de} \neq -1$), and baryonic and dark matter (with $\omega_{dm} \neq 0$), are investigated and compared. Using Bayesian machine learning analysis, constraints on the free parameters of both models are obtained for the three redshift ranges: $z\in [0,2]$, $z\in [0,2.5]$, and $z\in [0,5]$, respectively. For the first two redshift ranges, high-quality observations of the expansion rate $H(z)$ exist already, and they are used for validating the fitting results. Additionally, the extended range $z\in [0,5]$ provides predictions of the model parameters, verified when reliable higher-redshift $H(z)$ data are available. This learning procedure, based on the expansion rate data generated from the background dynamics of each model, shows that, at cosmological scales, there is a deviation from the cold dark matter paradigm, $\omega_{dm} \neq 0$, for all three redshift ranges. The results show that this approach may qualify as a solution to the $H_{0}$ tension problem. Indeed, it hints at how this issue could be effectively solved (or at least alleviated) in cosmological models with interacting dark energy.

It is not currently possible to create a living organism ab initio due to the overwhelming complexity of biological systems. In fact, the origin of life mechanism, this being how biological organisms form from non-living matter, is unknown. In an attempt to better understand how abiogenesis can occur, some researchers have taken water out of their models and instead opted for more exotic approaches. These assumptions will have strong implications for astronomical observations and potential future space exploration. By breaking down water's properties to the physical, chemical and biological level, herewith it is demonstrated to be the most adequate medium for the formation of life.

We have measured the scattering timescale, $\tau$, and the scattering spectral index, $\alpha$, for 84 single-component pulsars. Observations were carried out with the MeerKAT telescope as part of the Thousand-Pulsar-Array programme in the MeerTime project at frequencies between 0.895 and 1.670 GHz. Our results give a distribution of values for $\alpha$ (defined in terms of $\tau$ and frequency $\nu$ as $\tau\propto\nu^{-\alpha}$) for which, upon fitting a Gaussian, we obtain a mean and standard deviation of $\langle\alpha\rangle = 4.0 \pm 0.6$. This is due to our identification of possible causes of inaccurate measurement of $\tau$, which, if not filtered out of modelling results, tend to lead to underestimation of $\alpha$. The pulsars in our sample have large dispersion measures and are therefore likely to be distant. We find that a model using an isotropic scatter broadening function is consistent with the data, likely due to the averaging effect of multiple scattering screens along the line of sight. Our sample of scattering parameters provides a strong data set upon which we can build to test more complex and time-dependent scattering phenomena, such as extreme scattering events.

SARAS is an ongoing experiment aiming to detect the redshifted global 21-cm signal expected from Cosmic Dawn (CD) and the Epoch of Reionization (EoR). Standard cosmological models predict the signal to be present in the redshift range $z \sim $6--35, corresponding to a frequency range 40--200~MHz, as a spectral distortion of amplitude 20--200~mK in the 3~K cosmic microwave background. Since the signal might span multiple octaves in frequency, and this frequency range is dominated by strong terrestrial Radio Frequency Interference (RFI) and astrophysical foregrounds of Galactic and Extragalactic origin that are several orders of magnitude greater in brightness temperature, design of a radiometer for measurement of this faint signal is a challenging task. It is critical that the instrumental systematics do not result in additive or multiplicative confusing spectral structures in the measured sky spectrum and thus preclude detection of the weak 21-cm signal. Here we present the system design of the SARAS~3 version of the receiver. New features in the evolved design include Dicke switching, double differencing and optical isolation for improved accuracy in calibration and rejection of additive and multiplicative systematics. We derive and present the measurement equations for the SARAS~3 receiver configuration and calibration scheme, and provide results of laboratory tests performed using various precision terminations that qualify the performance of the radiometer receiver for the science goal.

The era of atmospheric characterization of terrestrial exoplanets is just around the corner. Modeling prior to observations is crucial in order to predict the observational challenges and to prepare for the data interpretation. This paper presents the report of the TRAPPIST Habitable Atmosphere Intercomparison (THAI) workshop (14-16 September 2020). A review of the climate models and parameterizations of the atmospheric processes on terrestrial exoplanets, model advancements and limitations, as well as direction for future model development was discussed. We hope that this report will be used as a roadmap for future numerical simulations of exoplanet atmospheres and maintaining strong connections to the astronomical community.

We present the analytical framework for converting projected light distributions with a S\'ersic profile into three-dimensional light distributions for stellar systems of arbitrary triaxial shape. The main practical result is the definition of a simple yet robust measure of intrinsic galaxy size: the median radius $r_\mathrm{med}$, defined as the radius of a sphere that contains 50% of the total luminosity or mass, that is, the median distance of a star to the galaxy center. We examine how $r_\mathrm{med}$ depends on projected size measurements as a function of S\'ersic index and intrinsic axis ratios, and demonstrate its relative independence of these parameters. As an application we show that the projected semi-major axis length of the ellipse enclosing 50% of the light is an unbiased proxy for $r_\mathrm{med}$, with small galaxy-to-galaxy scatter of $\sim$10% (1$\sigma$), under the condition that the variation in triaxiality within the population is small. For galaxy populations with unknown or a large range in triaxiality an unbiased proxy for $r_\mathrm{med}$ is $1.3\times R_{e}$, where $R_{e}$ is the circularized half-light radius, with galaxy-to-galaxy scatter of 20-30% (1$\sigma$). We also describe how inclinations can be estimated for individual galaxies based on the measured projected shape and prior knowledge of the intrinsic shape distribution of the corresponding galaxy population. We make the numerical implementation of our calculations available.

There is a special interest to understand the dynamical properties of the accretion disk created around the newly formed black hole due to the supermassive black hole binaries which merge inside the gaseous disk. The newly formed black hole would have a kick velocity that drives a perturbation on a newly accreted torus around the black hole. In this paper, the effects of the kicked black holes onto the accreted torus are studied by using the general relativistic hydrodynamical code, focusing on changing the dynamics of the accretion disk during the accretion disk-black hole interaction. We have found the non-axisymmetric global mode m=1 Papaloizou-Pringle Instability (PPI) produced on the torus due to kicked black hole. The higher the perturbation velocity produced by the kicked black hole, the longer time the torus takes to reach the saturation point. The created funnel walls are also observed from the numerical simulations. These funnels are responsible for accreting matter toward the black hole. At the later time of evolution, the accretion through the funnel is stopped and PPI is developed on the torus around the black hole.

We study the properties of the binary neutron star (BNS) systems in the inspiral phase. To calculate the equation of state (EOS) of the neutron star (NS), we take the relativistic mean-field (RMF) model. The RMF model, namely NL3 (stiff) and two extended RMF model IOPB-I (less stiff) and G3 (soft) are taken to explore the properties of the NS. We assume that the dark matter (DM) particles are accreted inside the NS due to its enormous gravitational field. Different macroscopic properties of the NS such as mass $M$, radius $R$, tidal deformability $\lambda$ and dimensionless tidal deformability $\Lambda$ are calculated at different DM fractions. With the addition of DM inside the NS, the value of the quantities like $M$, $R$, $\lambda$ and $\Lambda$ decreases. To explore the BNS properties in the inspiral phase, the post-Newtonian (PN) formalism is considered because it is suitable up to the last orbits in the inspiral phase. We calculated the strain amplitude of the polarization waveforms $h_+$ and $h_\times$, (2,2) mode waveform $h_{22}$, orbital phase $\Phi$, frequency of the gravitational wave $f$ and PN parameter $x$ with DM as an extra candidate inside the NS. We find that the BNS with soft EOS sustains more time in their inspiral phase as compare to stiff EOS. In the case of DM admixed NS, the BNS with high DM fractions survives more time in the inspiral phase than lesser fraction of DM. The magnitude of $f$, $\Phi$ and $x$ are almost the same for all the assumed parameter sets, but their inspiral time in the last orbit is different. We find a significant change in the BNS systems properties in the inspiral phase with DM inside the NS.

In this contribution I present results achieved recently in the field of the OT forecast that push further the limit of the accuracy of the OT forecasts and open to new perspectives in this field.

This paper separately evaluates the effects of inclination and asymmetry of solar coronal loops on the resonant absorption of kink magnetohydrodynamic (MHD) oscillations. We modelled a typical coronal loop by a straight and axisymmetric cylindrical magnetic flux tube filled with cold plasma. We solved the dispersion relation numerically for different values of the longitudinal mass density stratification. We show that, in inclined and asymmetric loops, the frequencies and their corresponding damping rates of the fundamental and first-overtone modes of kink oscillations are smaller in comparison with semi-circular uninclined loops with the same lengths. The results also indicate that, the period ratio $P_1/P_2$, increases with increasing the inclination of the loop, but it decreases less than $2\%$ while imposing the asymmetry to each loop side, up to $9.66\%$ of the loop length. The ratio of each mode frequency to its corresponding damping rate remain unchanged approximately while the inclination or the asymmetry imposed. Hence, we conclude that these ratios are reliable for inferring the physical parameters of coronal loops and coronal medium, regardless of the loop shape or the state of its inclination. In addition, in contrast with the effect of asymmetry which is not significant on the period ratio $P_1/P_2$, when an observed oscillating loop has a smaller apex height, the state of its inclination is an important factor that should be considered, especially when the period ratio $P_1/P_2$, is taken into consideration for coronal seismology.

Precise measurements of neutron star (NS) velocities provide critical clues to the supernova physics and evolution of binary systems. Based on Gaia Data Release 2 (DR2), we selected a sample of 24 young (<3 Myr) pulsars with precise parallax measurements and measured the velocity of their local standard of rest (LSR) and the velocity dispersion among their respective local stellar groups. The median velocity difference between thus calculated LSRs and the Galactic rotation model is ~7.6km/s, small compared to the typical velocity dispersion of ~27.5km/s. For pulsars off the Galactic plane, such differences grow significantly to as large as ~40 km/s. More importantly, the velocity dispersion of stars in the local group of low-velocity pulsars can be comparable to their transverse velocities, suggesting that the intrinsic velocities of NS progenitors should be taken into account when we consider their natal kicks and binary evolution. We also examined the double NS systems J0737- 3039A/B, and measured its transverse velocity to be 26(+18, -13) km/s assuming nearby Gaia sources being representative of its birth environment. This work demonstrated the feasibility and importance of using Gaia data to study the velocity of individual systems and velocity distribution of NSs.

We identify 228 periodic variables in the region of young open cluster NGC 281 using time series photometry carried out from 1 m class ARIES telescopes, Nainital. The cluster membership of these identified variables is determined on the basis colour-magnitude, two colour diagrams and kinematic data. Eighty one variable stars are found to be members of the cluster NGC 281. Of 81 variables, 30 and 51 are probable main-sequence and pre-main-sequence members, respectively. Present study classifies main-sequence variable stars into different variability types according to their periods of variability, shape of light curves and location in the Hertzsprung-Russell diagram. These identified main-sequence variables could be $\beta$ Cep, $\delta$ Scuti, slowly pulsating B type and new class variables. Among 51 pre-main-sequence variables, majority of them are weak line T Tauri stars. The remaining 147 variables could belong to the field population. The variability characteristics of the field population indicate that these variables could be RR Lyrae, $\delta$ Scuti and binaries type variables.

We present the XFaster analysis package. XFaster is a fast, iterative angular power spectrum estimator based on a diagonal approximation to the quadratic Fisher matrix estimator. XFaster uses Monte Carlo simulations to compute noise biases and filter transfer functions and is thus a hybrid of both Monte Carlo and quadratic estimator methods. In contrast to conventional pseudo-$C_\ell$ based methods, the algorithm described here requires a minimal number of simulations, and does not require them to be precisely representative of the data to estimate accurate covariance matrices for the bandpowers. The formalism works with polarization-sensitive observations and also data sets with identical, partially overlapping, or independent survey regions. The method was first implemented for the analysis of BOOMERanG data, and also used as part of the Planck analysis. Here, we describe the full, publicly available analysis package, written in Python, as developed for the analysis of data from the 2015 flight of the SPIDER instrument. The package includes extensions for self-consistently estimating null spectra and for estimating fits for Galactic foreground contributions. We show results from the extensive validation of XFaster using simulations, and its application to the SPIDER data set.

We present a study on the spatially scanned spectroscopic observations of the transit of GJ 1132 b, a warm ($\sim$500 K) Super-Earth (1.13 R$_\oplus$) that was obtained with the G141 grism (1.125 - 1.650 $\mu$m) of the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope. We used the publicly available Iraclis pipeline to extract the planetary transmission spectra from the five visits and produce a precise transmission spectrum. We analysed the spectrum using the TauREx3 atmospheric retrieval code with which we show that the measurements do not contain molecular signatures in the investigated wavelength range and are best-fit with a flat-line model. Our results suggest that the planet does not have a clear primordial, hydrogen-dominated atmosphere. Instead, GJ 1132 b could have a cloudy hydrogen-dominated envelope, a very enriched secondary atmosphere, be airless, or have a tenuous atmosphere that has not been detected. Due to the narrow wavelength coverage of WFC3, these scenarios cannot be distinguished yet but the James Webb Space Telescope may be capable of detecting atmospheric features, although several observations may be required to provide useful constraints.

A detailed X-ray Soft Excess of the narrow-line Seyfert-1 galaxy Tonantzintla S180 (Ton s180) present. The study used four XMM-Newton observations for period from the year 2000 to the year 2016 taken from the XMM-Newton archive. Two different models were used to treat the X-ray Soft Excess for each data separately, the first model is power law component and the second one is two black body components. We found that, for all observations, the spectra fitted by using a power law component are poor and the output parameters much different from the previous studies, this model unsuccessfully reproduces the continuum shape of the spectrum, including the soft excess, of all observations (Nh = 3.62 +-1020 cm^{-2} and photon index for the hard and soft bands are 3.31 and 1.57). The study proved that it can rely upon the two black body components in treating X-ray Soft Excess of the narrow-line Seyfert-1 galaxy (Nh= 1.395 +- 10^20 cm^{-2} and The temperature of the two black-body components (KT) are 0.075 and 0.17 keV).

Obtaining accurate photometric redshift estimations is an important aspect of cosmology, remaining a prerequisite of many analyses. In creating novel methods to produce redshift estimations, there has been a shift towards using machine learning techniques. However, there has not been as much of a focus on how well different machine learning methods scale or perform with the ever-increasing amounts of data being produced. Here, we introduce a benchmark designed to analyse the performance and scalability of different supervised machine learning methods for photometric redshift estimation. Making use of the Sloan Digital Sky Survey (SDSS - DR12) dataset, we analysed a variety of the most used machine learning algorithms. By scaling the number of galaxies used to train and test the algorithms up to one million, we obtained several metrics demonstrating the algorithms' performance and scalability for this task. Furthermore, by introducing a new optimisation method, time-considered optimisation, we were able to demonstrate how a small concession of error can allow for a great improvement in efficiency. From the algorithms tested we found that the Random Forest performed best in terms of error with a mean squared error, MSE = 0.0042; however, as other algorithms such as Boosted Decision Trees and k-Nearest Neighbours performed incredibly similarly, we used our benchmarks to demonstrate how different algorithms could be superior in different scenarios. We believe benchmarks such as this will become even more vital with upcoming surveys, such as LSST, which will capture billions of galaxies requiring photometric redshifts.

Fast and precise propagation of satellite orbits is required for mission design, orbit determination and payload data analysis. We present a method to improve the computational performance of numerical propagators and simultaneously maintain the accuracy level required by any particular application. This is achieved by determining the positional accuracy needed and the corresponding acceptable error in acceleration on the basis of the mission requirements, removing those perturbation forces whose effect is negligible compared to the accuracy requirement, implementing an efficient and precise algorithm for the harmonic synthesis of the geopotential gradient (i.e., the gravitational acceleration) and adjusting the tolerance of the numerical propagator to achieve the prescribed accuracy level with minimum cost. In particular, to achieve the optimum balance between accuracy and computational performance, the number of geopotential spherical harmonics to retain is adjusted during the integration on the basis of the accuracy requirement. The contribution of high-order harmonics decays rapidly with altitude, so the minimum expansion degree meeting the target accuracy decreases with height. The optimum degree for each altitude is determined by making the truncation error of the harmonic synthesis equal to the admissible acceleration error. This paper presents a detailed description of the technique and test cases highlighting its accuracy and efficiency.

The papers studies the highly misknown works of french astronomers on ephemerides and lunar distances on the period 1742-1785. These works will have a strong influence on the development of the Nautical Almanac by Nevil Maskelyne in 1766-1767.

The recent discovery of a Galactic fast radio burst (FRB) occurring simultaneously with an X-ray burst (XRB) from the Galactic magnetar SGR J1935+2154 implies that at least some FRBs arise from magnetar activities. We propose that FRBs are triggered by crust fracturing of magnetars, with the burst event rate depending on the magnetic field strength in the crust. Crust fracturing produces Alfv\'en waves, forming a charge starved region in the magnetosphere and leading to non-stationary pair plasma discharges. An FRB is produced by coherent plasma emission due to nonuniform pair production across magnetic field lines. Meanwhile, the FRB-associated XRB is produced by the rapid relaxation of the external magnetic field lines. In this picture, the sharp-peak hard X-ray component in association with FRB 200428 is from a region between adjacent trapped fireballs, and its spectrum with a high cutoff energy is attributed to resonant Compton scattering. The persistent X-ray emission is from a hot spot heated by the magnetospheric activities, and its temperature evolution is dominated by magnetar surface cooling. Within this picture, magnetars with stronger fields tend to produce brighter and more frequent repeated bursts.

We highlight shortcomings of the dynamical dark energy (DDE) paradigm. For parametric models with equation of state (EOS), $w(z) = w_0 + w_a f(z)$ for a given function of redshift $f(z)$, we show that the errors in $w_a$ are sensitive to $f(z)$: if $f(z)$ increases quickly with redshift $z$, then errors in $w_a$ are smaller, and vice versa. As a result, parametric DDE models suffer from a degree of arbitrariness and focusing too much on one model runs that risk that DDE may be overlooked. In particular, we show the ubiquitous Chevallier-Polarski-Linder model is one of the least sensitive to DDE. We also comment on ``wiggles" in $w(z)$ uncovered in non-parametric reconstructions. Concretely, we isolate the most relevant Fourier modes in the wiggles, model them and fit them back to the original data to confirm the wiggles at $\lesssim2\sigma$. We delve into the assumptions going into the reconstruction and argue that the assumed correlations, which clearly influence the wiggles, place strong constraints on field theory models of DDE.

This paper analytically and numerically investigates misalignment and mode-mismatch induced power coupling coefficients and losses as a function of Hermite Gauss (HG) mode order. We show that higher-order HG modes are more susceptible to beam perturbations: the misalignment and mode-mismatch induced power coupling losses scale linearly and quadratically with respect to the mode indices respectively. As a result, the mode-mismatch-limited tolerance for $\mathrm{HG}_{3,3}$ mode for example is reduced to a factor of 0.28 compared against the currently-used fundamental Gaussian laser mode. This is a potential hurdle for replacing the fundamental mode with higher-order modes in future gravitational-wave detectors.

We observationally examine cosmological models based on primordial power spectra with quantized wavevectors. Introducing a linearly quantized power spectrum with $k_0=3.225\times10^{-4}\mathrm{Mpc}^{-1}$ and spacing $\Delta k = 2.257 \times 10^{-4} \mathrm{Mpc}^{-1}$ provides a better fit to the Planck 2018 observations than the concordance baseline, with $\Delta \chi^2 = -8.55$. Extending the results of Lasenby et al [1], we show that the requirement for perturbations to remain finite beyond the future conformal boundary in a universe containing dark matter and a cosmological constant results in a linearly quantized primordial power spectrum. It is found that the infrared cutoffs for this future conformal boundary quantized cosmology do not provide cosmic microwave background power spectra compatible with observations, but future theories may predict more observationally consistent quantized spectra.

Warm Jupiters -- defined here as planets larger than 6 Earth radii with orbital periods of 8--200 days -- are a key missing piece in our understanding of how planetary systems form and evolve. It is currently debated whether Warm Jupiters form in situ, undergo disk or high eccentricity tidal migration, or have a mixture of origin channels. These different classes of origin channels lead to different expectations for Warm Jupiters' properties, which are currently difficult to evaluate due to the small sample size. We take advantage of the \TESS survey and systematically search for Warm Jupiter candidates around main-sequence host stars brighter than the \TESS-band magnitude of 12 in the Full-Frame Images in Year 1 of the \TESS Prime Mission data. We introduce a catalog of 55 Warm Jupiter candidates, including 19 candidates that were not originally released as \TESS Objects of Interest (TOIs) by the \TESS team. We fit their \TESS light curves, characterize their eccentricities and transit-timing variations (TTVs), and prioritize a list for ground-based follow-up and \TESS Extended Mission observations. Using hierarchical Bayesian modeling, we find the preliminary eccentricity distributions of our Warm-Jupiter-candidate catalog using a Beta distribution, a Rayleigh distribution, and a two-component Gaussian distribution as the functional forms of the eccentricity distribution. Additional follow-up observations will be required to clean the sample of false positives for a full statistical study, derive the orbital solutions to break the eccentricity degeneracy, and provide mass measurements.

We characterise the conditional distributions of the HI gas-to-stellar mass ratio, $R_{HI}\equiv M_{HI}/M_{\ast}$, given the stellar mass, $M_{\ast}$, of local galaxies from $M_{\ast}\sim 10^7$ to $10^{12}$ $M_{\odot}$ separated into centrals and satellites as well as into late- and early-type galaxies (LTGs and ETGs, respectively). To do so, we use 1) the homogeneous "eXtended GALEX Arecibo SDSS Survey", xGASS (Catinella et al. 2018), by re-estimating their upper limits and taking into account them in our statistical analysis; and 2) the results from a large compilation of HI data reported in Calette et al. (2018). We use the $R_{HI}$ conditional distributions combined with the Galaxy Stellar Mass Function to infer the bivariate $M_{HI}$ and $M_{\ast}$ distribution of all galaxies as well of the late/early-type and central/satellite subsamples and their combinations. Satellites are on average less HI gas-rich than centrals at low and intermediate masses, with differences being larger for ETGs than LTGs; at $M_{\ast}>3-5\times 10^{10}$ $M_{\odot}$ the differences are negligible. The differences in the HI gas content are much larger between LTGs and ETGs than between centrals and satellites. Our empirical HI Mass Function is strongly dominated by central galaxies at all masses. The empirically constrained bivariate $M_{HI}$ and $M_{\ast}$ distributions presented here can be used to compare and constrain theoretical predictions as well as to generate galaxy mock catalogs.

There are only few very-high-energy sources in our Galaxy which might accelerate particles up to the knee of the cosmic-ray spectrum. To understand the mechanisms of particle acceleration in these PeVatron candidates, \textit{Fermi}-LAT and H.E.S.S. observations are essential to characterize their $\gamma$-ray emission. HESS J1640$-$465 and the PeVatron candidate HESS J1641$-$463 are two neighboring (\ang[astroang]{0.25}) $\gamma$-ray sources, spatially coincident with the radio supernova remnants (SNRs) G338.3$-$0.0 and G338.5+0.1. Detected both by H.E.S.S. and \textit{Fermi}-LAT, we present here a morphological and spectral analysis of these two sources using 8 years of \textit{Fermi}-LAT data between 200 \si{\mega\electronvolt} and 1 \si{\tera\electronvolt} with multi-wavelength observations to assess their nature. The morphology of HESS J1640$-$465 is described by a 2D Gaussian ($\sigma=$ \ang[astroang]{0.053} $\pm$ \ang[astroang]{0.011}$_{stat}$ $ \pm$ \ang[astroang]{0.03}$_{syst}$) and its spectrum is modeled by a power-law with a spectral index $\Gamma = 1.8\pm0.1_{\rm stat}\pm0.2_{\rm syst}$. HESS J1641$-$463 is detected as a point-like source and its GeV emission is described by a logarithmic-parabola spectrum with $\alpha = 2.7 \pm 0.1_ {\rm stat} \pm 0.2_ {\rm syst} $ and significant curvature of $\beta = 0.11 \pm 0.03_ {\rm stat} \pm 0.05_ {\rm syst} $. Radio and X-ray flux upper limits were derived. We investigated scenarios to explain their emission, namely the emission from accelerated particles within the SNRs spatially coincident with each source, molecular clouds illuminated by cosmic rays from the close-by SNRs, and a pulsar/PWN origin. Our new \emph{Fermi}-LAT results and the radio and flux X-ray upper limits pose severe constraints on some of these models.

The importance of the chromosphere in the mass and energy transport within the solar atmosphere is now widely recognised. This review discusses the physics of magnetohydrodynamic (MHD) waves and instabilities in large-scale chromospheric structures as well as in magnetic flux tubes. We highlight a number of key observational aspects that have helped our understanding of the role of the solar chromosphere in various dynamic processes and wave phenomena, and the heating scenario of the solar chromosphere is also discussed. The review focuses on the physics of waves and invokes the basics of plasma instabilities in the context of this important layer of the solar atmosphere. Potential implications, future trends and outstanding questions are also delineated.

We study rotating Bose stars, i.e. gravitationally bound clumps of Bose-Einstein condensate composed of nonrelativistic particles with nonzero angular momentum $l$. We analytically prove that these objects are unstable at arbitrary $l\ne 0$ if particle self-interactions are attractive or negligibly small. In the latter case we calculate the profiles and dominant instability modes of the rotating stars: numerically at $1 \leq l \leq 15$ and analytically at $l\gg 1$. Notably, their lifetimes are always comparable to the inverse binding energies; hence, these objects cannot be considered long-living. Finally, we numerically show that in models with sufficiently strong repulsive self-interactions the Bose star with $l=1$ is stable.

Intense laser-plasma interactions are an essential tool for the laboratory study of ion acceleration at a collisionless shock. With two-dimensional particle-in-cell calculations of a multicomponent plasma we observe two electrostatic collisionless shocks at two distinct longitudinal positions when driven with a linearly-polarized laser at normalized laser vector potential a0 that exceeds 10. Moreover, these shocks, associated with protons and carbon ions, show a power-law dependence on a0 and accelerate ions to different velocities in an expanding upstream with higher flux than in a single-component hydrogen or carbon plasma. This results from an electrostatic ion two-stream instability caused by differences in the charge-to-mass ratio of different ions. Particle acceleration in collisionless shocks in multicomponent plasma are ubiquitous in space and astrophysics, and these calculations identify the possibility for studying these complex processes in the laboratory.

We show that gravitational theories with a nonminimal coupling (NMC) to the matter fields lead to a violation of Etherington's distance-duality relation, which relates the luminosity and angular diameter distances. We derive constraints on power-law and exponential NMC models using existing measurements of type Ia supernovae and baryon acoustic oscillations throughout the redshift range $0<z<1.5$. These complement previous constrains derived from cosmic-microwave background radiation and big-bang nucleosynthesis data.

Radio astronomy commenced in earnest after World War II, with Australia keenly engaged through the Council for Scientific and Industrial Research. At this juncture, Australia's Commonwealth Solar Observatory expanded its portfolio from primarily studying solar phenomena to conducting stellar and extragalactic research. Subsequently, in the 1950s and 1960s, astronomy gradually became taught and researched in Australian universities. However, most scientific publications from this era of growth and discovery have no country of affiliation in their header information, making it hard to find the Australian astronomy articles from this period. In 2014, we used the then-new Astrophysics Data System (ADS) tool Bumblebee to overcome this challenge and track down the Australian-led astronomy papers published during the quarter of a century after World War II, from 1945 until the lunar landing in 1969. This required knowledge of the research centres and facilities operating at the time, which are briefly summarised herein. Based on citation counts -- an objective, universally-used measure of scientific impact -- we report on the Australian astronomy articles which had the biggest impact. We have identified the top-ten most-cited papers, and thus also their area of research, from five consecutive time-intervals across that blossoming quarter-century of astronomy. Moreover, we have invested a substantial amount of time researching and providing a small tribute to each of the 62 scientists involved, including several trail-blazing women. Furthermore, we provide an extensive list of references and point out many interesting historical connections and anecdotes.

The Navier-Stokes equations generate an infinite set of generalized Lyapunov exponents defined by different ways of measuring the distance between exponentially diverging perturbed and unperturbed solutions. This set is demonstrated to be similar, yet different, from the generalized Lyapunov exponent that provides moments of distance between two fluid particles below the Kolmogorov scale. We derive rigorous upper bounds on dimensionless Lyapunov exponent of the fluid particles that demonstrate the exponent's decay with Reynolds number $Re$ in accord with previous studies. In contrast, terms of cumulant series for exponents of the moments have power-law growth with $Re$. We demonstrate as an application that the growth of small fluctuations of magnetic field in ideal conducting turbulence is hyper-intermittent, being exponential in both time and Reynolds number. We resolve the existing contradiction between the theory, that predicts slow decrease of dimensionless Lyapunov exponent of turbulence with $Re$, and observations exhibiting quite fast growth. We demonstrate that it is highly plausible that a pointwise limit for the growth of small perturbations of the Navier-Stokes equations exists.

Using multipoint Magnetospheric Multiscale (MMS) observations in an unusual string-of-pearls configuration, we examine in detail observations of the reformation of a fast magnetosonic shock observed on the upstream edge of a foreshock transient structure upstream of Earth's bow shock. The four MMS spacecraft were separated by several hundred km, comparable to suprathermal ion gyro-radius scales or several ion inertial lengths. At least half of the shock reformation cycle was observed, with a new shock ramp rising up out of the "foot" region of the original shock ramp. Using the multipoint observations, we convert the observed time-series data into distance along the shock normal in the shock's rest frame. That conversion allows for a unique study of the relative spatial scales of the shock's various features, including the shock's growth rate, and how they evolve during the reformation cycle. Analysis indicates that: the growth rate increases during reformation, electron-scale physics play an important role in the shock reformation, and energy conversion processes also undergo the same cyclical periodicity as reformation. Strong, thin electron-kinetic-scale current sheets and large-amplitude electrostatic and electromagnetic waves are reported. Results highlight the critical cross-scale coupling between electron-kinetic- and ion-kinetic-scale processes and details of the nature of nonstationarity, shock-front reformation at collisionless, fast magnetosonic shocks.

Extending the Standard Model with three right-handed neutrinos and a simple QCD axion sector can account for neutrino oscillations, dark matter and baryon asymmetry; at the same time, it solves the strong CP problem, stabilize the electroweak vacuum and can implement critical Higgs inflation (satisfying all current observational bounds). We perform here a general analysis of dark matter (DM) in such a model, which we call the $a\nu$MSM. Although critical Higgs inflation features a (quasi) inflection point of the inflaton potential we show that DM cannot receive a contribution from primordial black holes in the $a\nu$MSM. This leads to a multicomponent axion-sterile-neutrino DM and allows us to relate the axion parameters, such as the axion decay constant, to the neutrino parameters. We include several DM production mechanisms: the axion production via misalignment and decay of topological defects as well as the sterile-neutrino production through the resonant and non-resonant mechanisms and in the recently proposed CPT-symmetric universe.

Macroscopic dark matter (or macro) provides a broad class of alternative candidates to particle dark matter. These candidates would transfer energy primarily through elastic scattering, and this linear energy deposition would produce observable signals if a macro were to traverse the atmosphere. We study the fluorescence emission produced by a macro passing through the atmosphere. We estimate the sensitivity of EUSO-SPB2 to constrain the two-dimensional parameter space ($\sigma$ vs. $M$), where $M$ is the macro mass and $\sigma$ its cross sectional area.

We present a study of neutron star models that contain dark matter (DM) in the core. The DM is assumed to have a particle nature and to be self-interacting. Using constraints on the mass and radius of neutron stars, we investigate the allowed properties of either bosonic or fermionic DM particles. For this purpose three different models of neutron stars are considered, the first involving nucleons only, the second including hyperons, and the last involving strange matter in the core. Different EoS are constructed for the various cases of fermionic and bosonic DM. These EoSs are solved for selected properties of the DM particles and the results are tested against mass, radius and tidal deformability constraints for neutron stars. The distribution of energy density of DM and ordinary matter inside the neutron stars is also presented. It is found that if the DM is fermionic in nature it does not just sit in the core but it is present everywhere in the star, from the centre to outside the surface and may even envelop.

We study the change in internal rotational energy in the transformation of protons to neutrons in the \b{eta}-decay reactions that take place in the collapse of the iron core of massive stars that precede type II supernova explosions. We consider an ensemble of electrons, protons, neutrons and neutrinos undergoing \b{eta}-decay reactions, treat the protons and neutrons as triatomic rotors, evaluate the equilibrium constant to obtain the change in rotational energy during the proton-to-neutron transformation. We estimate such change for a variety of conditions, and compare to the energy released in a supernova explosion.

A variety of kinetic waves develop in the solar wind. The relationship between these waves and larger-scale structures, such as current sheets and ongoing turbulence remain a topic of investigation. Similarly, the instabilities producing ion-acoustic waves in the solar wind remains an open question. The goals of this paper are to investigate kinetic electrostatic Langmuir and ion-acoustic waves in the solar wind at 0.5 AU and determine whether current sheets and associated streaming instabilities can produce the observed waves. The relationship between these waves and currents is investigated statistically. Solar Orbiter's Radio and Plasma Waves instrument suite provides high-resolution snapshots of the fluctuating electric field. The Low Frequency Receiver resolves the waveforms of ion-acoustic waves and the Time Domain Sampler resolves the waveforms of both ion-acoustic and Langmuir waves. Using these waveform data we determine when these waves are observed in relation to current structures in the solar wind, estimated from the background magnetic field. Langmuir and ion-acoustic waves are frequently observed in the solar wind. Ion-acoustic waves are observed about 1% of the time at 0.5 AU. The waves are more likely to be observed in regions of enhanced currents. However, the waves typically do not occur at current structures themselves. The observed currents in the solar wind are too small to drive instability by the relative drift between single ion and electron populations. When multi-component ion and/or electron distributions are present the observed currents may be sufficient for instability. Ion beams are the most plausible source of ion-acoustic waves. The spacecraft potential is confirmed to be a reliable probe of the background electron density by comparing the peak frequencies of Langmuir waves with the plasma frequency calculated from the spacecraft potential.

The no-hair conjecture in General Relativity (GR) states that a Kerr black hole (BH) is completely described by its mass and spin. As a consequence, the complex quasi-normal-mode (QNM) frequencies of a binary-black-hole (BBH) ringdown can be uniquely determined by the mass and spin of the remnant object. Conversely, measurement of the QNM frequencies could be an independent test of the no-hair conjecture. This paper extends to spinning BHs earlier work that proposed to test the no-hair conjecture by measuring the complex QNM frequencies of a BBH ringdown using parameterized inspiral-merger-ringdown waveforms in the effective-one-body formalism, thereby taking full advantage of the entire signal power and removing dependency on the predicted or estimated start time of the ringdown. Our method was used to analyze the properties of the merger remnants for BBHs observed by LIGO-Virgo in the first half of their third observing (O3a) run. After testing our method with GR and non-GR synthetic-signal injections in Gaussian noise, we analyze, for the first time, two BBHs from the first (O1) and second (O2) LIGO-Virgo observing runs, and two additional BBHs from the O3a run. We then provide joint constraints with published results from the O3a run. In the most agnostic and conservative scenario where we combine the information from different events using a hierarchical approach, we obtain, at $90\%$ credibility, that the fractional deviations in the frequency (damping time) of the dominant QNM are $\delta f_{220}=0.03^{+0.10}_{-0.09}$ ($\delta \tau_{220}=0.10^{+0.44}_{-0.39}$), respectively, an improvement of a factor of $\sim 4$ ($\sim 2$) over the results obtained with our model in the LIGO-Virgo publication. The single-event most-stringent constraint to date continues to be GW150914 for which we obtain $\delta f_{220}=0.05^{+0.11}_{-0.07}$ and $\delta \tau_{220}=0.07^{+0.26}_{-0.23}$.

On 14th September 2020, the Royal Astronomical Society made an official statement coupled with a webminar on the discovery of phosphine on Venus. Single-line millimetre-waveband spectral detections of phosphine (with a signal-to-noise ratio of $\approx$ 15$\sigma$) from the JCMT and ALMA telescopes indicated a phosphine abundance of 20 ppb (parts per billion), 1000 times more than that on the Earth. Phosphine is an important biomarker and immediate speculation in the media about indicators of life being found on Venus followed. This article presents an analysis of the study and the results on the observation of the spectral absorption feature of phosphine in the clouds of Venus, thus implying as a potential biosignature. If phosphine is produced through biotic, as opposed to abiotic pathways, the discovery could imply a significant biomass in the Venusian atmosphere. The discovery led to a major controversy with criticism of the analysis and results and responses to it. The issue remains unresolved, leading to a fresh interest in the study of Venus including ground-based observations as well as space-probes that can answer these questions conclusively.

One-loop correction to the power spectrum in generic single-field inflation is calculated extracting the most important operator, which also generates primordial non-Gaussianity, using soft effective field theory. Due to the enhancement inversely proportional to the observed red-tilt of the spectral index of curvature perturbation, the correction turns out to be much larger than previously anticipated. As a result, the primordial non-Gaussianity must be much smaller than the current observational bound in order to warrant the validity of cosmological perturbation theory.

The Higgs scalar which is the only experimentally verified fundamental cosmological scalar, is also known to produce a viable inflationary phenomenology. In this work we investigate the effects of the Higgs model on static neutron stars. Particularly we derive the Einstein frame Tolman-Oppenheimer-Volkoff equations, and by numerically integrating them for both the interior and the exterior of the neutron star, using a double shooting python 3 based numerical code, we extract the masses and radii of the neutron stars, along with the several other related physical quantities of interest. With regard to the equation of state for the neutron star, we use a piecewise polytropic equation of state with the central part being SLy, APR or the WFF1 equations of state. The resulting $M-R$ graphs are compatible with the observational bounds imposed by the GW170817 event which require the radius of a static $M\sim 1.6 M_{\odot}$ neutron star to be larger than $R=10.68^{+15}_{-0.04}$km and the radius of a static neutron star corresponding to the maximum mass of the star to be larger than $R=9.6^{+0.14}_{-0.03}$km. Moreover, the WFF1 EoS, which was excluded for static neutron stars in the context of general relativity, for the Higgs neutron star model provides realistic results compatible with the GW170817 event.

The NANOGrav collaboration has recently presented its pulsar timing array data which seem compatible with the presence of a stochastic gravity wave background emitted by cosmic strings with a dimensionless string tension $G\mu\simeq 2\times 10^{-11}-3\times 10^{-10}$ at $95\%$ confidence level ($G$ is Newton's constant and $\mu$ denotes the string tension). However, there is some tension between these results and the previous pulsar timing array bound $G\mu\lesssim 4\times 10^{-11}$ from the PPTA experiment. We propose a relaxation of this tension by invoking primordial inflation which partially inflates the string network. The latter re-enters the horizon at later times after the end of inflation, and thus the short string loops are not produced. This leads to a reduction of the gravity wave spectrum which is more pronounced at higher frequencies. The reconciliation of the NANOGrav results with the PPTA bound is possible provided that the strings re-enter the horizon at adequately late times. We consider an example of a realistic $SO(10)$ model incorporating successful inflation driven by a gauge singlet real scalar field with a Coleman-Weinberg potential. This model leads to the production of intermediate scale topologically stable cosmic strings that survive inflation. We show regions of the parameter space where the tension between NANOGrav and PPTA is alleviated. Finally, we present an example in which both monopoles and strings survive inflation with the above tension resolved.