Papers for Friday, Jun 04 2021

Ultra-compact X-ray binaries (UCXBs) are low-mass X-ray binaries with hydrogen-deficient mass-donors and ultra-short orbital periods. They have been suggested to be the potential Laser Interferometer Space Antenna (LISA) sources in the low-frequency region. Several channels for the formation of UCXBs have been proposed so far. In this article, we carried out a systematic study on the He star donor channel, in which a neutron star (NS) accretes matter from a He main-sequence star through Roche-lobe overflow, where the mass-transfer is driven by the gravitational wave radiation. Firstly, we followed the long-term evolution of the NS+He main-sequence star binaries by employing the stellar evolution code Modules for Experiments in Stellar Astrophysics, and thereby obtained the initial parameter spaces for the production of UCXBs. We then used these results to perform a detailed binary population synthesis approach to obtain the Galactic rates of UCXBs through this channel. We estimate the Galactic rates of UCXBs appearing as LISA sources to be $\sim3.1-11.9\times10^{-6} \rm yr^{-1}$ through this channel, and the number of such UCXB-LISA sources in the Galaxy can reach about $100-390$. The present work indicates that the He star donor channel cannot be ignored in forming UCXB-LISA sources. We found that the evolutionary tracks of UCXBs through this channel can account for the location of the five transient sources with relatively long orbital periods quite well. We also found that such UCXBs can be identified by their locations in the period-mass transfer rate diagram.

Tidal disruption events (TDEs) occur when a star gets torn apart by the strong tidal forces of a supermassive black hole, which results in the formation of a debris stream that partly falls back towards the compact object. This gas moves along inclined orbital planes that intersect near pericenter, resulting in a so-called "nozzle shock". We perform the first dedicated study of this interaction, making use of a two-dimensional simulation that follows the transverse gas evolution inside a given section of stream. This numerical approach circumvents the lack of resolution encountered near pericenter passage in global three-dimensional simulations using particle-based methods. As it moves inward, we find that the gas motion is purely ballistic, which near pericenter causes strong vertical compression that squeezes the stream into a thin sheet. Dissipation takes place at the resulting nozzle shock, inducing a rise in pressure that causes the collapsing gas to bounce back, although without imparting significant net expansion. As it recedes to larger distances, this matter continues to expand while remaining thin despite the influence of pressure forces. This gas evolution specifies the strength of the subsequent self-crossing shock, which we find to be more affected by black hole spin than previously estimated. We also evaluate the impact of general-relativistic effects, viscous dissipation, magnetic fields and radiative processes on the nozzle shock. This study represents an important step forward in the theoretical understanding of TDEs, bridging the gap between our robust knowledge of the fallback rate and the more complex following stages, during which most of the emission occurs.

Observed rotation curves in star-forming galaxies indicate a puzzling dearth of dark matter in extended flat cores within haloes of mass $\geq\! 10^{12}M_\odot$ at $z\!\sim\! 2$. This is not reproduced by current cosmological simulations, and supernova-driven outflows are not effective in such massive haloes. We address a hybrid scenario where post-compaction merging satellites heat up the dark-matter cusps by dynamical friction, allowing AGN-driven outflows to generate cores. Using analytic and semi-analytic models (SatGen), we estimate the dynamical-friction heating as a function of satellite compactness for a cosmological sequence of mergers. Cosmological simulations (VELA) demonstrate that satellites of initial virial masses $>\!10^{11.3}M_\odot$, that undergo wet compactions, become sufficiently compact for significant heating. Constituting a major fraction of the accretion onto haloes $\geq\!10^{12}M_\odot$, these satellites heat-up the cusps in half a virial time at $z\!\sim\! 2$. Using a model for outflow-driven core formation (CuspCore), we demonstrate that the heated dark-matter cusps develop extended cores in response to removal of half the gas mass, while the more compact stellar systems remain intact. The mergers keep the dark matter hot, while the gas supply, fresh and recycled, is sufficient for the AGN outflows. AGN indeed become effective in haloes $\geq\!10^{12}M_\odot$, where the black-hole growth is no longer suppressed by supernovae and its compaction-driven rapid growth is maintained by a hot CGM. For simulations to reproduce the dynamical-friction effects, they should resolve the compaction of the massive satellites and avoid artificial tidal disruption. AGN feedback could be boosted by clumpy black-hole accretion and clumpy response to AGN.

The flux ratios of high-ionization lines are commonly assumed to indicate the metallicity of the broad emission line region in luminous quasars. When accounting for the variation in their kinematic profiles, we show that the NV/CIV, (SiIV+OIV])/CIV and NV/Lya line ratios do not vary as a function of the quasar continuum luminosity, black hole mass, or accretion rate. Using photoionization models from CLOUDY , we further show that the observed changes in these line ratios can be explained by emission from gas with solar abundances, if the physical conditions of the emitting gas are allowed to vary over a broad range of densities and ionizing fluxes. The diversity of broad line emission in quasar spectra can be explained by a model with emission from two kinematically distinct regions, where the line ratios suggest that these regions have either very different metallicity or density. Both simplicity and current galaxy evolution models suggest that near-solar abundances, with parts of the spectrum forming in high-density clouds, are more likely. Within this paradigm, objects with stronger outflow signatures show stronger emission from gas which is denser and located closer to the ionizing source, at radii consistent with simulations of line-driven disc-winds. Studies using broad-line ratios to infer chemical enrichment histories should consider changes in density and ionizing flux before estimating metallicities.

Context. We report the discovery of VVV-CL160, a new nearby globular cluster (GC) with extreme kinematics, located in the Galactic plane at $l = 10.1477$ deg, $b = 0.2999$ deg. Aims. We aim to characterize the physical properties of this new GC and place it in the context of the Milky Way, exploring its possible connection with the known GC NGC 6544 and with the Hrid halo stream. Methods. VVV-CL160 was originally detected in the VISTA Variables in the V\'ia L\'actea (VVV) survey. We use the proper motions (PMs) from the updated VVV Infrared Astrometric Catalog (VIRAC2) to select GC members and make deep near-infrared color-magnitude diagrams (CMDs) to study the cluster properties. We also fit King models to the decontaminated sample to determine the GC structural parameters. Results. VVV-CL160 has an unusually large PM for a Galactic GC as measured with VIRAC2 and Gaia EDR3: $\mu_{\alpha}\cos(\delta)$ = $-2.3 \pm 0.1 $ mas yr$^{-1}$ and $\mu_{\delta}$ = $-16.8 \pm 0.1 $ mas yr$^{-1}$. The kinematics are similar to those of the known GC NGC 6544 and the Hrid halo stream. We estimate a reddening of $E(J-K) = 1.95$ mag and an extinction of $A_{k}= 1.40$ mag for VVV-CL160. We also measure a distance modulus of $(m-M) = 13.01$ mag and a distance of $D_{\odot} = 4.0 \pm 0.5$ kpc. This places the GC at $z=29$ pc above the Galactic plane and at a galactocentric distance of $R_G=4.2$ kpc. We also measure a metallicity of $[Fe/H] = -1.4 \pm 0.2$ dex for an adopted age of $t=12$ Gyr; King model fits of the PM-decontaminated sample reveal a concentrated GC, with core radius $r_{c}= 22.8"$ and tidal radius $r_{t}= 50'$. .... We also explore the possible association of this new GC with other GCs and halo streams. Conclusions. Based on the locations and kinematics, we suggest that VVV-CL160, along with NGC 6544, may be associated with the extension of the Hrid halo stream.

Strong gravitational lensing, which can make a background source galaxy appears multiple times due to its light rays being deflected by the mass of one or more foreground lens galaxies, provides astronomers with a powerful tool to study dark matter, cosmology and the most distant Universe. PyAutoLens is an open-source Python 3.6+ package for strong gravitational lensing, with core features including fully automated strong lens modeling of galaxies and galaxy clusters, support for direct imaging and interferometer datasets and comprehensive tools for simulating samples of strong lenses. The API allows users to perform ray-tracing by using analytic light and mass profiles to build strong lens systems. Accompanying PyAutoLens is the autolens workspace (see https://github.com/Jammy2211/autolens_workspace), which includes example scripts, lens datasets and the HowToLens lectures in Jupyter notebook format which introduce non experts to strong lensing using PyAutoLens. Readers can try PyAutoLens right now by going to the introduction Jupyter notebook on Binder (see https://mybinder.org/v2/gh/Jammy2211/autolens_workspace/master) or checkout the readthedocs (see https://pyautolens.readthedocs.io/en/latest/) for a complete overview of PyAutoLens's features.

The recent discovery of electromagnetic signals in coincidence with neutron-star mergers has solidified the importance of multimessenger campaigns in studying the most energetic astrophysical events. Pioneering multimessenger observatories, such as LIGO/Virgo and IceCube, record many candidate signals below the detection significance threshold. These sub-threshold event candidates are promising targets for multimessenger studies, as the information provided by them may, when combined with contemporaneous gamma-ray observations, lead to significant detections. Here we describe a new method that uses such candidates to search for transient events using archival very-high-energy gamma-ray data from imaging atmospheric Cherenkov telescopes (IACTs). We demonstrate the application of this method to sub-threshold binary neutron star (BNS) merger candidates identified in Advanced LIGO's first observing run. We identify eight hours of archival VERITAS observations coincident with seven BNS merger candidates and search them for TeV emission. No gamma-ray emission is detected; we calculate upper limits on the integral flux and compare them to a short gamma-ray burst model. We anticipate this search method to serve as a starting point for IACT searches with future LIGO/Virgo data releases as well as in other sub-threshold studies for multimessenger transients, such as IceCube neutrinos. Furthermore, it can be deployed immediately with other current-generation IACTs, and has the potential for real-time use that places minimal burden on experimental operations. Lastly, this method may serve as a pilot for studies with the Cherenkov Telescope Array, which has the potential to observe even larger fields of view in its divergent pointing mode.

Schneider et al. (2020) presented the discovery of WISEA J041451.67-585456.7 and WISEA J181006.18-101000.5, which appear to be the first examples of extreme T-type subdwarfs (esdTs; metallicity <= -1 dex, T_eff <= 1400 K). Here we present new discoveries and follow-up of three T-type subdwarf candidates, with an eye toward expanding the sample of such objects with very low metallicity and extraordinarily high kinematics, properties that suggest membership in the Galactic halo. Keck/NIRES near-infrared spectroscopy of WISEA J155349.96+693355.2, a fast-moving object discovered by the Backyard Worlds: Planet 9 citizen science project, confirms that it is a mid-T subdwarf. With H_W2 = 22.3 mag, WISEA J155349.96+693355.2 has the largest W2 reduced proper motion among all spectroscopically confirmed L and T subdwarfs, suggesting that it may be kinematically extreme. Nevertheless, our modeling of the WISEA J155349.96+693355.2 near-infrared spectrum indicates that its metallicity is only mildly subsolar. In analyzing the J155349.96+693355.2 spectrum, we present a new grid of low-temperature, low-metallicity model atmosphere spectra. We also present the discoveries of two new esdT candidates, CWISE J073844.52-664334.6 and CWISE J221706.28-145437.6, based on their large motions and colors similar to those of the two known esdT objects. Finding more esdT examples is a critical step toward mapping out the spectral sequence and observational properties of this newly identified population.

We present three-dimensional simulations of core-collapse supernovae using the FLASH code that follow the progression of the explosion to the stellar surface, starting from neutrino-radiation hydrodynamic simulations of the first seconds performed with the CHIMERA code. We consider a 9.6-$M_{\odot}$ zero-metallicity progenitor starting from both 2D and 3D CHIMERA models, and a 10-$M_{\odot}$ solar-metallicity progenitor starting from a 2D CHIMERA model, all simulated until shock breakout in 3D while tracking 160 nuclear species. The relative velocity difference between the supernova shock and the metal-rich Rayleigh-Taylor (R-T) "bullets" determines how the ejecta evolves as it propagates through the density profile of the progenitor and dictates the final morphology of the explosion. We find maximum $^{56}\rm{Ni}$ velocities of ${\sim} 1950~\rm{km~s}^{-1}$ and ${\sim} 1750~\rm{km~s}^{-1}$ at shock breakout from 2D and 3D 9.6-$M_{\odot}$ CHIMERA models, respectively, due to the bullets' ability to penetrate the He/H shell. When mapping from 2D, we find that the development of higher velocity structures is suppressed when the 2D CHIMERA model and 3D FLASH model meshes are aligned. The development of faster growing spherical-bubble structures, as opposed to the slower growing toroidal structure imposed by axisymmetry, allows for interaction of the bullets with the shock and seeds further R-T instabilities at the He/H interface. We see similar effects in the 10-$M_{\odot}$ model, which achieves maximum $^{56}\rm{Ni}$ velocities of ${\sim} 2500~\rm{km~s}^{-1}$ at shock breakout.

Since July 2014, the Gaia mission has been engaged in a high-spatial-resolution, time-resolved, precise, accurate astrometric, and photometric survey of the entire sky. Aims: We present the Gaia Science Alerts project, which has been in operation since 1 June 2016. We describe the system which has been developed to enable the discovery and publication of transient photometric events as seen by Gaia. Methods: We outline the data handling, timings, and performances, and we describe the transient detection algorithms and filtering procedures needed to manage the high false alarm rate. We identify two classes of events: (1) sources which are new to Gaia and (2) Gaia sources which have undergone a significant brightening or fading. Validation of the Gaia transit astrometry and photometry was performed, followed by testing of the source environment to minimise contamination from Solar System objects, bright stars, and fainter near-neighbours. Results: We show that the Gaia Science Alerts project suffers from very low contamination, that is there are very few false-positives. We find that the external completeness for supernovae, $C_E=0.46$, is dominated by the Gaia scanning law and the requirement of detections from both fields-of-view. Where we have two or more scans the internal completeness is $C_I=0.79$ at 3 arcsec or larger from the centres of galaxies, but it drops closer in, especially within 1 arcsec. Conclusions: The per-transit photometry for Gaia transients is precise to 1 per cent at $G=13$, and 3 per cent at $G=19$. The per-transit astrometry is accurate to 55 milliarcseconds when compared to Gaia DR2. The Gaia Science Alerts project is one of the most homogeneous and productive transient surveys in operation, and it is the only survey which covers the whole sky at high spatial resolution (subarcsecond), including the Galactic plane and bulge.

In the Milky Way, $\sim$18 Wolf-Rayet+O (WR+O) binaries are known with estimates of their stellar and orbital parameters. Whereas black hole+O (BH+O) binaries are thought to evolve from the former, only one such system is known in the Milky Way. To resolve this disparity, it was suggested that upon core collapse, the WR stars receive large kicks such that most of the binaries are disrupted. We reassess this issue, with emphasis on the uncertainty in the formation of an accretion disk around wind-accreting BHs in BH+O binaries, which is key to identifying such systems. We follow the methodology of previous work and apply an improved analytic criterion for the formation of an accretion disk around wind accreting BHs. We then use stellar models to predict the properties of the BH+O binaries which are expected to descend from the observed WR+O binaries, if the WR stars would form BHs without a natal kick. We find that disk formation depends sensitively on the O stars' wind velocity, the specific angular momentum carried by the wind, the efficiency of angular momentum accretion by the BH, and the spin of the BH. We show that the assumption of a low wind velocity may lead to predicting that most of the BH+O star binaries will have an extended X-ray bright period. However, this is not the case when typical wind velocities of O stars are considered. We find that a high spin of the BH can boost the duration of the X-ray active phase as well as the X-ray brightness during this phase, producing a strong bias for detecting high mass BH binaries in X-rays with high BH spin parameters. We conclude that large BH formation kicks are not required to understand the sparsity of X-ray bright BH+O stars in the Milky Way. Probing for a population of X-ray silent BH+O systems with alternative methods can inform us about BH kicks and the conditions for high energy emission from high mass BH binaries. (Abridged)

We performed a systematic study of 12 active regions (ARs) with a broad range of areas, magnetic flux and associated solar activity in order to determine whether there are upflows present at the AR boundaries and if these upflows exist, whether there is a high speed asymmetric blue wing component present in the upflows. To identify the presence and locations of the AR upflows we derive relative Doppler velocity maps by fitting a Gaussian function to {\it Hinode}/EIS Fe XII 192.394\,\AA\ line profiles. To determine whether there is a high speed asymmetric component present in the AR upflows we fit a double Gaussian function to the Fe XII 192.394\,\AA\ mean spectrum that is computed in a region of interest situated in the AR upflows. Upflows are observed at both the east and west boundaries of all ARs in our sample with average upflow velocities ranging between -5 to -26~km s$^{-1}$. A blue wing asymmetry is present in every line profile. The intensity ratio between the minor high speed asymmetric Gaussian component compared to the main component is relatively small for the majority of regions however, in a minority of cases (8/30) the ratios are large and range between 20 to 56~\%. These results suggest that upflows and the high speed asymmetric blue wing component are a common feature of all ARs.

We present the flux density measurements of the pulsars observed in the Meterwavelength single-pulse polarimetric emission survey. The average flux densities were estimated in 113 pulsars at two frequencies of 325 and 610 MHz using interferometric imaging. The average profile and single pulse emission in each pulsar were calibrated using the estimated flux density. We have used the flux calibrated average profile to study the variation of the spectral index across the emission beam in 21 pulsars where the core, inner cone and the outer conal components could be clearly identified. The central core component showed a steeper increase in emission at the lower frequency compared with conal emission, with an average difference in spectral index $\delta\alpha_{core-cone}\sim-0.7$ between the core and the conal components in this frequency range. In contrast the inner conal components had positive difference in their spectral index compared to the outer cones with average difference $\delta\alpha_{in-out}\sim+0.3$. The variation in the spectral index across the pulse window should provide valuable inputs for constraining the radio emission processes. The single pulse emission showed the presence of emission mode changing in 12 pulsars with 3 cases where the phenomenon is being reported for the first time. In addition we have also detected enhanced emission for short durations or flaring, in parts or across the entire emission window in 14 pulsars. The sudden changes in the emission during mode changing as well as these bursting states are unrelated to the emission mechanism and suggest the presence of rapid and repetitive changes during the plasma generation process.

Among the most extreme objects in the Universe, active galactic nuclei (AGN) are luminous centers of galaxies where a black hole feeds on surrounding matter. The variability patterns of the light emitted by an AGN contain information about the physical properties of the underlying black hole. Upcoming telescopes will observe over 100 million AGN in multiple broadband wavelengths, yielding a large sample of multivariate time series with long gaps and irregular sampling. We present a method that reconstructs the AGN time series and simultaneously infers the posterior probability density distribution (PDF) over the physical quantities of the black hole, including its mass and luminosity. We apply this method to a simulated dataset of 11,000 AGN and report precision and accuracy of 0.4 dex and 0.3 dex in the inferred black hole mass. This work is the first to address probabilistic time series reconstruction and parameter inference for AGN in an end-to-end fashion.

We investigate the orientation of the photospheric magnetic fields in the solar polar region using observations from the Helioseismic and Magnetic Imager (HMI). Inside small patches of significant polarization, the inferred magnetic field vectors at $1''$ scale appear to systematically deviate from the radial direction. Most tilt towards the pole; all are more inclined toward the plane of sky compared to the radial vector. These results, however, depend on the "filling factor" $f$ that characterizes the unresolved magnetic structures. The default, uninformative $f\equiv1$ for HMI will incur larger inclination and less radial fields than $f<1$. The observed trend may be a systematic bias inherent to the limited resolution.

Mid-infrared (MIR) observations are typically accomplished from the ground through oscillating the secondary mirror a few times a second. This chopping serves to remove the fast time variable components of (a) sky variation and (b) array background. However, there is a significant price to pay for this, including reduced on-object photon collection time, stringent demands on the secondary mirror, nodding the telescope to remove the radiative offset imprinted by the chopping, and an often-fixed chop-frequency regardless of the sky conditions in the actual observations. Worse, in the era of 30m telescopes it is wholly impracticable to chop the secondary mirror. If the array is stable enough, drift scanning holds the promise to remove the necessity of chopping. In this paper we report our experiments using the CanariCam MIR instrument on the 10.4m GranTeCan and the implications to future instruments and experiments.

Although galactic winds play a critical role in regulating galaxy formation, hydrodynamic cosmological simulations do not resolve the scales that govern the interaction between winds and the ambient circumgalactic medium (CGM). We implement the Physically Evolved Wind (PhEW) model of Huang et al. (2020) in the GIZMO hydrodynamics code and perform test cosmological simulations with different choices of model parameters and numerical resolution. PhEW adopts an explicit subgrid model that treats each wind particle as a collection of clouds that exchange mass, metals, and momentum with their surroundings and evaporate by conduction and hydrodynamic instabilities as calibrated on much higher resolution cloud scale simulations. In contrast to a conventional wind algorithm, we find that PhEW results are robust to numerical resolution and implementation details because the small scale interactions are defined by the model itself. Compared to conventional wind simulations with the same resolution, our PhEW simulations produce similar galaxy stellar mass functions at $z\geq 1$ but are in better agreement with low-redshift observations at $M_* < 10^{11}M_\odot$ because PhEW particles shed mass to the CGM before escaping low mass halos. PhEW radically alters the CGM metal distribution because PhEW particles disperse metals to the ambient medium as their clouds dissipate, producing a CGM metallicity distribution that is skewed but unimodal and is similar between cold and hot gas. While the temperature distributions and radial profiles of gaseous halos are similar in simulations with PhEW and conventional winds, these changes in metal distribution will affect their predicted UV/X-ray properties in absorption and emission.

We present observations from the First Light Infrared TEst CAMera (FLITECAM) on board the Stratospheric Observatory for Infrared Astronomy (SOFIA), the Spitzer Infrared Array Camera (IRAC) and the Spitzer Infrared Spectrograph (IRS) SH mode in three well-known Photodissocation Regions (PDRs), the reflection nebulae (RNe) NGC 7023 and NGC 2023 and to the southeast of the Orion Bar, which are well suited to probe emission from Polycyclic Aromatic Hydrocarbon molecules (PAHs). We investigate the spatial behaviour of the FLITECAM 3.3 um filter as a proxy for the 3.3 um PAH band, the integrated 11.2 um PAH band, and the IRAC 8.0 um filter as a proxy for the sum of the 7.7 and 8.6 um PAH bands. The resulting ratios of 11.2/3.3 and IRAC 8.0/11.2 provide an approximate measure of the average PAH size and PAH ionization respectively. In both RNe, we find that the relative PAH ionization and the average PAH size increases with decreasing distance to the illuminating source. The average PAH sizes derived for NGC 2023 are greater than those found for NGC 7023 at all points. Both results indicate that PAH size is dependent on the radiation field intensity. These results provide additional evidence of a rich carbon-based chemistry driven by the photo-chemical evolution of the omnipresent PAH molecules within the interstellar medium. In contrast, we did not detect a significant variation in the average PAH size found in the region southeast of the Orion Bar and report a peculiar PAH ionization radial profile.

The discovery of amino acids in meteorites has presented two clues to the origin of their processing subsequent to their formation: a slight preference for left-handedness in some of them, and isotopic anomalies in some of their constituent atoms. In this article we present theoretical results from the Supernova Neutrino Amino Acid Processing (SNAAP) model, which uses electron anti-neutrinos and the magnetic fields from source objects such as supernovae or colliding neutron stars to selectively destroy one amino acid chirality and to create isotopic abundance shifts. For plausible magnetic fields and electron anti-neutrino fluxes, non-zero, positive enantiomeric excesses, $ee$s, defined to be the relative left/right asymmetry in an amino acid population, are reviewed for two amino acids, and conditions are suggested that would produce $ee>0$ for all of the $\alpha$-amino acids. The relatively high energy anti-neutrinos that produce the $ee$s would inevitably also produce isotopic anomalies. A nuclear reaction network was developed to describe the reactions resulting from them and the nuclides in the meteorites. At similar anti-neutrino fluxes, assumed recombination of the detritus from the anti-neutrino interactions is shown to produce appreciable isotopic anomalies in qualitative agreement with those observed for D/$^1$H and $^{15}$N/$^{14}$N. The isotopic anomalies for $^{13}$C/$^{12}$C are predicted to be small, as are also observed. Autocatalysis may be necessary for any model to produce the largest $ee$s observed in meteorites. This allows the constraints of the original SNAAP model to be relaxed, increasing the probability of meteoroid survival in sites where amino acid processing is possible. These results have obvious implications for the origin of life on Earth.

The recent publication by Siraj & Loeb (2021; Nature Scientific Reports 11, 3803) attempts to revive the debate over whether the Chicxulub impactor was a comet or an asteroid. They calculate that ~20% of long-period comets impacting Earth will have first been disrupted by passage inside the Sun's Roche limit, generating thousands of fragments, each the needed size of the Chicxulub impactor. This would increase the impact rate of comets by a factor ~15, making them as likely to hit the Earth as an asteroid. They also argue that a comet would be a factor of 10 more likely to match the geochemical constraints, which indicate the Chicxulub impactor was carbonaceous chondrite-like. These conclusions are based on misinterpretations of the literature. Siraj & Loeb [1] overestimate the number of fragments produced during tidal disruption of a comet: tens of fragments are produced, not thousands. They also conflate 'carbonaceous chondrite' with specific types of carbonaceous chondrite, and ignore the evidence of iridium, making comets seem more likely than asteroids to match the Chicxulub impactor, when in fact they likely can be ruled out. Rather than a comet, an asteroidal impactor similar to CM or CR carbonaceous chondrites is strongly favored.

In this paper, a fast and parallelizable method based on Gaussian Processes (GPs) is introduced to emulate computer models that simulate the formation of binary black holes (BBHs) through the evolution of pairs of massive stars. Two obstacles that arise in this application are the a priori unknown conditions of BBH formation and the large scale of the simulation data. We address them by proposing a local emulator which combines a GP classifier and a GP regression model. The resulting emulator can also be utilized in planning future computer simulations through a proposed criterion for sequential design. By propagating uncertainties of simulation input through the emulator, we are able to obtain the distribution of BBH properties under the distribution of physical parameters.

The Zwicky Transient Facility recently announced the detection of an optical transient AT2020blt at redshift $z=2.9$, consistent with the afterglow of a gamma-ray burst. No prompt emission was observed. We analyse AT2020blt with detailed models, showing the data are best explained as the afterglow of an on-axis long gamma-ray burst, ruling out other hypotheses such as a cocoon and a low-Lorentz factor jet. We search \textit{Fermi} data for prompt emission, setting deeper upper limits on the prompt emission than in the original detection paper. Together with \konus{} observations, we show that the gamma-ray efficiency of AT2020blt is $\lesssim 2.8\%$, lower than $98.4\%$ of observed gamma-ray bursts. We speculate that AT2020blt and AT2021any belong to the low-efficiency tail of long gamma-ray burst distributions that are beginning to be readily observed due to the capabilities of new observatories like the Zwicky Transient Facility.

Planets grow via the collisional accretion of small bodies in a protoplanetary disk. Such small bodies feel strong gas drag and their orbits are significantly affected by the gas flow and atmospheric structure around the planet. We investigate the gas flow in the protoplanetary disk perturbed by the gravity of the planet by three-dimensional hydrodynamic simulation. We then calculate the orbital evolutions of particles in the gas structure obtained from the hydrodynamic simulation. Based on the orbital calculations, we obtain the collision rate between the planet and centimeter to kilometer sized particles. Our results show that meter-sized or larger particles effectively collide with the planet due to the atmospheric gas drag, which significantly enhances the collision rate. On the other hand, the gas flow plays an important role for smaller particles. Finally, considering the effects of the atmosphere and gas flow, we derive the new analytic formula for the collision rate, which is in good agreement with our simulations. We estimate the growth timescale and accretion efficiency of drifting bodies for the formation of a gas-giant solid core using the formula. We find the accretion of sub-kilometer sized bodies achieve a short growth timescale (~0.05 Myr) and a high accretion efficiency (~1) for the core formation at 5 au in the minimum mass solar nebula model.

Classifying the morphologies of galaxies is an important step in understanding their physical properties and evolutionary histories. The advent of large-scale surveys has hastened the need to develop techniques for automated morphological classification. We train and test several convolutional neural network architectures to classify the morphologies of galaxies in both a 3-class (elliptical, lenticular, spiral) and 4-class (+irregular/miscellaneous) schema with a dataset of 14034 visually-classified SDSS images. We develop a new CNN architecture that outperforms existing models in both 3 and 4-way classification, with overall classification accuracies of 83% and 81% respectively. We also compare the accuracies of 2-way / binary classifications between all four classes, showing that ellipticals and spirals are most easily distinguished (>98% accuracy), while spirals and irregulars are hardest to differentiate (78% accuracy). Through an analysis of all classified samples, we find tentative evidence that misclassifications are physically meaningful, with lenticulars misclassified as ellipticals tending to be more massive, among other trends. We further combine our binary CNN classifiers to perform a hierarchical classification of samples, obtaining comparable accuracies (81%) to the direct 3-class CNN, but considerably worse accuracies in the 4-way case (65%). As an additional verification, we apply our networks to a small sample of Galaxy Zoo images, obtaining accuracies of 92%, 82% and 77% for the binary, 3-way and 4-way classifications respectively.

XRISM is an X-ray astronomical mission by the JAXA, NASA, ESA and other international participants, that is planned for launch in 2022 (Japanese fiscal year), to quickly restore high-resolution X-ray spectroscopy of astrophysical objects. To enhance the scientific outputs of the mission, the Science Operations Team (SOT) is structured independently from the instrument teams and the Mission Operations Team. The responsibilities of the SOT are divided into four categories: 1) guest observer program and data distributions, 2) distribution of analysis software and the calibration database, 3) guest observer support activities, and 4) performance verification and optimization activities. As the first step, lessons on the science operations learned from past Japanese X-ray missions are reviewed, and 15 kinds of lessons are identified. Among them, a) the importance of early preparation of the operations from the ground stage, b) construction of an independent team for science operations separate from the instrument development, and c) operations with well-defined duties by appointed members are recognized as key lessons. Then, the team structure and the task division between the mission and science operations are defined; the tasks are shared among Japan, US, and Europe and are performed by three centers, the SOC, SDC, and ESAC, respectively. The SOC is designed to perform tasks close to the spacecraft operations, such as spacecraft planning, quick-look health checks, pre-pipeline processing, etc., and the SDC covers tasks regarding data calibration processing, maintenance of analysis tools, etc. The data-archive and user-support activities are covered both by the SOC and SDC. Finally, the science-operations tasks and tools are defined and prepared before launch.

Since its launch, the Alpha Magnetic Spectrometer-02 (AMS-02) has delivered outstanding quality measurements of the spectra of cosmic-ray (CR) species, $\bar{p}$, $e^{\pm}$, and nuclei (H-O, Ne, Mg, Si, Fe), which resulted in a number of breakthroughs. The most recent AMS-02 result is the measurement of the spectrum of CR fluorine up to $\sim$2 TV. Given its very low solar system abundance, fluorine in CRs is thought to be mostly secondary, produced in fragmentations of heavier species, predominantly Ne, Mg, and Si. Similar to the best-measured secondary-to-primary boron to carbon nuclei ratio that is widely used to study the origin and propagation of CR species, the precise fluorine data would allow the origin of Si-group nuclei to be studied independently. Meanwhile, the secondary origin of CR fluorine has never been tested in a wide energy range due to the lack of accurate CR data. In this paper, we use the first ever precise measurements of the fluorine spectrum by AMS-02 together with ACE-CRIS and Voyager 1 data to actually test this paradigm. Our detailed modeling shows an excess below 10 GV in the fluorine spectrum that is most likely due to the primary fluorine component. We also provide an updated local interstellar spectrum (LIS) of fluorine in the rigidity range from few MV to $\sim$2 TV. Our calculations employ the self-consistent GalProp-HelMod framework that has proved to be a reliable tool in deriving the LIS of CR $\bar{p}$, $e^{-}$, and nuclei $Z\le28$.

We analytically study the dynamical and thermal properties of the optically-thin gases at the parsec-scale when they are spherically accreted onto low luminous active galactic nuclei (LLAGNs). The falling gases are irradiated by the central X-ray radiation with the Compton temperature of 5--15$\times10^7$ K. The radiative heating/cooling and the bulge stellar potential in galaxies are taken into account. We analyze the effect of accretion rate, luminosity, gas temperature, and Compton temperature on steady solutions of dynamical and thermal properties. The steady solutions are obviously different from Bondi solution. Compared to our models, the Bondi model underestimates the accretion rate. We give the boundary between thermal stability and instability. The boundary is significantly affected by Compton temperature. When Compton temperature is higher, the falling gases easily thermally unstable. When thermal instability takes place in the irradiated gases, the gases become two phases (i.e. hot gases and cool gases) and the hot gases may become outflows. This effect may reduce the accretion rates.

Since the gravitational waves were detected by LIGO and VIRGO, it has been promising that lots of information about the primordial Universe could be learned by further observations on stochastic gravitational waves background. The studies on gravitational waves induced by primordial curvature perturbations are of great interest. The aim of this paper is to investigate the third order induced gravitational waves. Based on the theory of cosmological perturbations, the first order scalar induces the second order scalar, vector and tensor perturbations. At the next iteration, the first order scalar, the second order scalar, vector and tensor perturbations all induce the third order tensor perturbations. We present the energy density spectrum of the third order gravitational waves for a monochromatic primordial power spectrum. The shape of the energy density spectrum of the third order gravitational waves is different from that of the second order scalar induced gravitational waves. And it is found that the third order gravitational waves sourced by the second order scalar perturbations dominate the energy density spectrum.

NGC 253 is a starburst galaxy of SAB(s)c type with increasing interest because of its high activity at unrivaled closeness. Its energetic event is manifested as the vertical gas features in its central molecular zone, for which stellar feedback was proposed as the driving engine. In order to pursue details of the activity, we have undertaken a kinematic analysis of the ALMA archive data of CO($J$=3--2) emission at the highest resolution $\sim$3 pc. We revealed that one of the non-rotating gas components in the central molecular zone shows a loop-like distribution of $\sim$200 pc radius. The loop is associated with a star cluster, whereas the cluster is not inside the loop and is not likely as the driver of the loop formation. Further, we find that the bar potential of NGC 253 seems to be too weak to drive the gas motion by the eccentric orbit. As an alternative we frame a scenario that magnetic acceleration by the Parker instability is responsible for the creation of the loop. We show that the observed loop properties are similar to those in the Milky Way, and argue that recent magnetro-hydrodynamics simulations lend support for the picture having the magnetic field strength of $\gtrsim$100 $\mu$G. We suggest that cluster formation was triggered by the falling gas to the footpoint of the loop, which is consistent with a typical dynamical timescale of the loop $\sim$1 Myr.

Jets and outflows are thought to play important roles in regulating star formation and disk evolution. HD 163296 is a well-studied Herbig Ae star that hosts proto-planet candidates, a protoplanetary disk, a protostellar jet, and a molecular outflow, which makes it an excellent laboratory for studying jets. We aim to characterize the jet at the inner regions and check if there are large differences with the features at large separations. A secondary objective is to demonstrate the performance of Multi Unit Spectroscopic Explorer (MUSE) in high-contrast imaging of extended line emission. MUSE in the narrow field mode (NFM) can provide observations at optical wavelengths with high spatial ($\sim$75 mas) and medium spectral ($R\sim$2500) resolution. With the high-resolution spectral differential imaging (HRSDI) technique, we can characterize the kinematic structures and physical conditions of jets down to 100 mas. We detect multiple atomic lines in two new knots, B3 and A4, at distances of <4" from the host star with MUSE. The derived $\dot{M}_{\rm jet} / \dot{M}_{\rm acc}$ is about 0.08 and 0.06 for knots B3 and A4, respectively. The observed [Ca II]/[S II] ratios indicate that there is no sign of dust grains at distances of <4". Assuming the knot A4 traces the streamline, we set an upper limit of 2.2 au on the size of the launching region. Although MUSE has the ability to detect the velocity shifts caused by high- and low-velocity components, we found no significant evidence of velocity decrease transverse to the jet direction. Our work demonstrates the capability of using MUSE NFM observations for the detailed study of stellar jets in the optical down to 100~mas. The derived $\dot{M}_{\rm jet} / \dot{M}_{\rm acc}$, no dust grain, and jet radius at the star support the magneto-centrifugal models as a launching mechanism for the jet.

Starquakes probably occur in rapidly spinning or ultra high field neutron stars. In this short article, we argue that highly compressed gas containing electron-positron pairs could evaporate and erupt from inside the neutron star when a crack forms and then heals during a starquake. Under the influence of the existing oscillation modes of the star, the crack may have sufficiently large size and long lifetime. An appropriate amount of gas can erupt into the magnetosphere with relativistic and nearly uniform velocity, producing various transient and bursting phenomena.

We calculate Galactic Chemical Evolution (GCE) of Mo and Ru by taking into account the contribution from $\nu p$-process nucleosynthesis. We estimate yields of $p$-nuclei such as $^{92,94}\mathrm{Mo}$ and $^{96,98}\mathrm{Ru}$ through the $\nu p$-process in various supernova (SN) progenitors based upon recent models. In particular, the $\nu p$-process in energetic hypernovae produces a large amount of $p$-nuclei compared to the yield in ordinary core-collapse SNe. Because of this the abundances of $^{92,94}\mathrm{Mo}$ and $^{96,98}\mathrm{Ru}$ in the Galaxy are significantly enhanced at [Fe/H]=0 by the $\nu p$-process. We find that the $\nu p$-process in hypernovae is the main contributor to the elemental abundance of $^{92}$Mo at low metallicity [Fe/H$]<-2$. Our theoretical prediction of the elemental abundances in metal-poor stars becomes more consistent with observational data when the $\nu p$-process in hypernovae is taken into account.

This paper presents the construction of the RapidXMM database that is available through the XMM-Newton Science Archive and offers access to upper limits and aperture photometry across the field of view of the XMM-Newton Pointed and Slew Survey observations. The feature of RapidXMM is speed. It enables the fast retrieval of X-ray upper limits and photometry products in three energy bands (0.2-2, 2-12, 0.2-12 keV) for large numbers of input sky positions. This is accomplished using the Hierarchical Equal Area Iso Latitude pixelation of the sphere (HEALPix). The pre-calculated upper-limits and associated X-ray photometry products are reprojected into the HEALPix grid of cells before being ingested into the RapidXMM database. This results in tables of upper limits and aperture photometry within HEALPix cells of size ~3 arcsec (Pointed Observations) and ~6 arcsec (Slew Survey). The database tables are indexed by the unique integer number of the HEALPix cells. This reduces spatial nearest-neighbor queries by sky position to an integer-matching exercise and significantly accelerates the retrieval of results. We describe in detail the processing steps that lead from the science products available in the XMM-Newton archive to a database optimised for sky queries. We also present two simple show-case applications of RapidXMM for scientific studies: searching for variable X-ray sources, and stacking analysis of X-ray faint populations

As recent advancements in physics and astronomy rapidly rewrite textbooks, there is a growing need in keeping abreast of the latest knowledge in these fields. Reading preprints is one of the effective ways to do this. By having journal clubs where people can read and discuss journals together, the benefits of reading journals become more prevalent. We present an investigative study of understanding the factors that affect the success of preprint journal clubs in astronomy, more commonly known as Astro-ph/Astro-Coffee (hereafter called AC). A survey was disseminated to understand how institutions from different countries implement AC. We interviewed 9 survey respondents and from their responses we identified four important factors that make AC successful: commitment (how the organizer and attendees participate in AC), environment (how conducive and comfortable AC is conducted), content (the discussed topics in AC and how they are presented), and objective (the main goal/s of conducting AC). We also present the format of our AC, an elective class which was evaluated during the Spring Semester 2020 (March 2020 - June 2020). Our evaluation with the attendees showed that enrollees (those who are enrolled and are required to present papers regularly) tend to be more committed in attending compared to audiences (those who are not enrolled and are not required to present papers regularly). In addition, participants tend to find papers outside their research field harder to read. Finally, we showed an improvement in the weekly number of papers read after attending AC of those who present papers regularly, and a high satisfaction rating of our AC. We summarize the areas of improvement in our AC implementation, and we encourage other institutions to evaluate their own AC in accordance with the four aforementioned factors to assess the effectiveness of their AC in reaching their goals.

We report a $\gamma$-ray enhancement event detected from Tycho's supernova remnant (SNR), the outcome of a type Ia supernova explosion that occurred in year 1572. The event lasted for 1.5 years and showed a factor of 3.6 flux increase mainly in the energy range of 4--100 GeV. While several young SNRs (including Tycho's SNR) were previously found to show peculiar X-ray structures with flux variations in one- or several-year timescales, such an event at $\gamma$-ray energies is for the first time seen. The hard $\gamma$-ray emission and year-long timescale of the event necessitate a synchrotron radiation process, although the required conditions are either ultra-high energies for the electrons in the process, upto $\sim$10 PeV (well above the cosmic-ray "knee" energy), or high inhomogeneity of the magnetic field in the SNR. This event in Tycho's SNR is likely analogous to the $\gamma$-ray flares observed in the Crab nebula, the comparably short timescales of them both requiring a synchrotron process, and similar magnetohydrodynamic processes such as magnetic reconnection would be at work as well in the SNR to accelerate particles to ultra-relativistic energies. The event and its implications thus reveal the more complicated side of the physical processes that can occur in young SNRs.

We reconstruct the history of reionization using Gaussian process regression. Using the UV luminosity data compilation from Hubble Frontiers Fields we reconstruct the redshift evolution of UV luminosity density and thereby the evolution of the source term in the ionization equation. This model-independent reconstruction rules out single power-law evolution of the luminosity density but supports the logarithmic double power-law parametrization. We obtain reionization history by integrating ionization equations with the reconstructed source term. Using optical depth constraint from Planck Cosmic Microwave Background observation, measurement of UV luminosity function integrated till truncation magnitude of -17 and -15, and derived ionization fraction from high redshift quasar, galaxies and gamma-ray burst observations, we constrain the history of reionization. In the conservative case we find the constraint on the optical depth as $\tau =0.052\pm0.001\pm0.002$ at 68% and 95% confidence intervals. We find the redshift duration between 10% and 90% ionization to be $2.05_{-0.21-0.30}^{+0.11+0.37}$. Longer duration of reionization is supported if UV luminosity density data with truncation magnitude of -15 is used in the joint analysis. Our results point out that even in a conservative reconstruction, a combination of cosmological and astrophysical observations can provide stringent constraints on the epoch of reionization.

Sunspots are known to be strong absorbers of solar oscillation modal power. The most convincing way to demonstrate this is done via Fourier-Hankel decomposition (FHD), where the local oscillation field is separated into in- and outgoing waves, showing the reduction in power. Due to the Helioseismic and Magnetic Imager's high-cadence Doppler measurements, power absorption can be investigated at frequencies beyond the acoustic cutoff frequency. We perform an FHD on five sunspot regions and two quiet-Sun control regions and study the resulting absorption spectra $\alpha_\ell(\nu)$, specifically at frequencies $\nu$ > 5.3 mHz. We observe an unreported high-frequency absorption feature, which only appears in the presence of a sunspot. This feature is confined to phase speeds of one-skip waves whose origins coincide with the sunspot's center, with $v_{ph}$ = 85.7 km/s in this case. By employing a fit to the absorption spectra at a constant phase speed, we find that the peak absorption strength $\alpha_{max}$ lies between 0.166 and 0.222 at a noise level of about 0.009 (5%). The well-known absorption along ridges at lower frequencies can reach up to $\alpha_{max}\approx$ 0.5. Thus our finding in the absorption spectrum is weaker, but nevertheless significant. From first considerations regarding the energy budget of high-frequency waves, this observation can likely be explained by the reduction of emissivity within the sunspot. We derive a simple relation between emissivity and absorption. We conclude that sunspots yield a wave power absorption signature (for certain phase speeds only), which may help in understanding the effect of strong magnetic fields on convection and source excitation and potentially in understanding the general sunspot subsurface structure.

The full third Gaia data release will provide the calibrated spectra obtained with the blue and red Gaia slit-less spectrophotometers. The main challenge when facing Gaia spectral calibration is that no lamp spectra or flat fields are available during the mission. Also, the significant size of the line spread function with respect to the dispersion of the prisms produces alien photons contaminating neighbouring positions of the spectra. This makes the calibration special and different from standard approaches. This work gives a detailed description of the internal calibration model to obtain the spectrophotometric data in the Gaia catalogue. The main purpose of the internal calibration is to bring all the epoch spectra onto a common flux and pixel (pseudo-wavelength) scale, taking into account variations over the focal plane and with time, producing a mean spectrum from all the observations of the same source. In order to describe all observations in a common mean flux and pseudo-wavelength scale, we construct a suitable representation of the internally calibrated mean spectra via basis functions and we describe the transformation between non calibrated epoch spectra and calibrated mean spectra via a discrete convolution, parametrising the convolution kernel to recover the relevant coefficients. The model proposed here is able to combine all observations into a mean instrument to allow the comparison of different sources and observations obtained with different instrumental conditions along the mission and the generation of mean spectra from a number of observations of the same source. The output of this model provides the internal mean spectra, not as a sampled function (flux and wavelength), but as a linear combination of basis functions, although sampled spectra can easily be derived from them.

Many exoplanets are discovered in binary star systems in internal or in circumbinary orbits. Whether the planet can be habitable or not depends on the possibility to maintain liquid water on its surface, and therefore on the luminosity of its host stars and on the dynamical properties of the planetary orbit. The trajectory of a planet in a double star system can be determined, approximating stars and planets with point masses, by solving numerically the equations of motion of the classical three-body system. In this study, we analyze a large data set of planetary orbits, made up with high precision long integration at varying: the mass of the planet, its distance from the primary star, the mass ratio for the two stars in the binary system, and the eccentricity of the star motion. To simulate the gravitational dynamics, we use a 15th order integration scheme (IAS15, available within the REBOUND framework), that provides an optimal solution for long-term integration. In our data analysis, we evaluate if an orbit is stable or not and also provide the statistics of different types of instability: collisions with the primary or secondary star and planets ejected away from the binary star system. Concerning the stability, we find a significant number of orbits that are only marginally stable, according to the classification introduced by Musielak et al. 2005. For planets of negligible mass, we estimate the critical semi-major axis $a_c$ as a function of the mass ratio and the eccentricity of the binary, in agreement with the results of Holman and Wiegert 1999. However, we find that for very massive planets (Super-Jupiters) the critical semi-major axis decrease in some cases by a few percent, compared to cases in which the mass of the planet is negligible.

As we enter the era of large-scale imaging surveys with the up-coming telescopes such as LSST and SKA, it is envisaged that the number of known strong gravitational lensing systems will increase dramatically. However, these events are still very rare and require the efficient processing of millions of images. In order to tackle this image processing problem, we present Machine Learning techniques and apply them to the Gravitational Lens Finding Challenge. The Convolutional Neural Networks (CNNs) presented have been re-implemented within a new modular, and extendable framework, LEXACTUM. We report an Area Under the Curve (AUC) of 0.9343 and 0.9870, and an execution time of 0.0061s and 0.0594s per image, for the Space and Ground datasets respectively, showing that the results obtained by CNNs are very competitive with conventional methods (such as visual inspection and arc finders) for detecting gravitational lenses.

Gravitational lensing has long been used to measure or constrain cosmology models. Although the lensing effect of gravitational waves has not been observed by LIGO/Virgo, it is expected that there can be a few to a few hundreds lensed events to be detected by the future Japanese space-borne interferometers DECIGO and B-DECIGO, if they are running for 4 years. Given the predicted lensed gravitational wave events, one can estimate the constraints on the cosmological parameters via the lensing statistics and the time delay methods. With the lensing statistics method, the knowledge of the lens redshifts, even with the moderate uncertainty, will set the tight bound on the energy density parameter $\Omega_M$ for matter, that is, $0.288\lesssim\Omega_M\lesssim0.314$ at best. The constraint on the Hubble constant $H_0$ can be determined using the time delay method. It is found out that at $5\sigma$, $|\delta H_0|/H_0$ ranges from $3\%$ to $11\%$ for DECIGO, and B-DECIGO will give less constrained results, $8\%-15\%$. In this work, the uncertainties on the luminosity distance and the time delay distance are set to be $10\%$ and $20\%$, respectively. The improvement on measuring these distances will tighten the bounds.

The study of narrow-line Seyfert 1 galaxies (NLS1s) is now mostly limited to low redshift ($z<0.8$) because their definition requires the presence of the H$\beta$ emission line, which is redshifted out of the spectral coverage of major ground-based spectroscopic surveys at $z>0.8$. We studied the correlation between the properties of H$\beta$ and Mg II lines of a large sample of SDSS DR14 quasars to find high-$z$ NLS1 candidates. Based on the strong correlation of $\mathrm{FWHM(MgII)=(0.880\pm 0.005) \times FWHM(H\beta)+ (0.438\pm0.018)}$, we present a sample of high-$z$ NLS1 candidates having FWHM of Mg II $<$ 2000 km s$^{-1}$. The high-$z$ sample contains 2684 NLS1s with redshift $z=0.8-2.5$ with a median logarithmic bolometric luminosity of $46.16\pm0.42$ erg s$^{-1}$, logarithmic black hole mass of $8.01\pm0.35 M_{\odot}$, and logarithmic Eddington ratio of $0.02\pm0.27$. The fraction of radio-detected high-$z$ NLS1s is similar to that of the low-$z$ NLS1s and SDSS DR14 quasars at a similar redshift range, and their radio luminosity is found to be strongly correlated with their black hole mass.

We report the systematic analysis of knots, hotspots, and lobes in 57 active galactic nuclei (AGNs) to investigate the variation of the magnetic field along with the jet from the sub-pc base to the terminus in kpc-to-Mpc scales. Expanding the number of radio/X-ray samples in Kataoka & Stawarz (2005), we analyzed the data in 12 FR I and 30 FR II radio galaxies, 12 quasars, and 3 BL Lacs that contained 76 knots, 42 hotspots, and 29 radio lobes. We first derived the equipartition magnetic fields in the cores and then estimated those in various jet components by assuming $B_{\rm est}$ $\propto$ $d^{-1}$, where $d$ is the distance from the jet base. On the other hand, the magnetic field in large-scale jets (knots, hotspots, and lobes), $B_{\rm eq}$, can be estimated from the observed flux and spatial extent under the equipartition hypothesis. We show that the magnetic field decreases as the distance along the jet increases, but generally gentler than $\propto d^{-1}$. The increase in $B_{\rm eq}/B_{\rm est}$ at a larger $d$ may suggest the deceleration of the jet around the downstream, but there is no difference between FR I and FR II jets. Moreover, the magnetic fields in the hotspots are systematically larger than those of knots and lobes. Finally, we applied the same analysis to knots and lobes in Centaurus A to check whether the above discussion will hold even in a single jet source.

The temperatures of red supergiants (RSGs) are expected to depend on metallicity (Z) in such a way that lower-Z RSGs are warmer. In this work, we investigate the Z-dependence of the Hayashi limit by analysing RSGs in the low-Z galaxy Wolf-Lundmark-Mellote (WLM), and compare with the RSGs in the higher-Z environments of the Small Magellanic Cloud (SMC) and Large Magellanic Cloud (LMC). We determine the effective temperature ($T_{\textrm{eff}}$) of each star by fitting their spectral energy distributions, as observed by VLT+SHOOTER, with MARCS model atmospheres. We find average temperatures of $T_{\textrm{eff}_{\textrm{WLM}}}=4400\pm202$ K, $T_{\textrm{eff}_{\textrm{SMC}}}=4130\pm103$ K, and $T_{\textrm{eff}_{\textrm{LMC}}}=4140\pm148$ K. From population synthesis analysis, we find that although the Geneva evolutionary models reproduce this trend qualitatively, the RSGs in these models are systematically too cool. We speculate that our results can be explained by the inapplicability of the standard solar mixing length to RSGs.

We have carried out a new kinematical analysis of the molecular gas in the Sh2-233 region by using the CO $J$ = 2-1 data taken at $\sim$0.5 pc resolution. The molecular gas consists of a filamentary cloud of 5-pc length with 1.5-pc width where two dense cloud cores are embedded. The filament lies between two clouds, which have a velocity difference of 2.6 km s$^{-1}$ and are extended over $\sim$5 pc. We frame a scenario that the two clouds are colliding with each other and compressed the gas between them to form the filament in $\sim$0.5 Myr which is perpendicular to the collision. It is likely that the collision formed not only the filamentary cloud but also the two dense cores. One of the dense cores is associated with the high-mass protostellar candidate IRAS 05358+3543, a representative high-mass protostar. In the monolithic collapse scheme of high mass star formation, a compact dense core of 100 $M_\odot$ within a volume of 0.1 pc radius is assumed as the initial condition, whereas the formation of such a core remained unexplained in the previous works. We argue that the proposed collision is a step which efficiently collects the gas of 100 $M_\odot$ into 0.1 pc radius. This lends support for that the cloud-cloud collision is an essential process in forming the compact high-mass dense core, IRAS 05358+3543.

Similar to Earth, Saturn's largest moon, Titan, possesses a system of high-altitude zonal winds (or jets) that encircle the globe. Using the Atacama Large Millimeter/submillimeter Array (ALMA) in August 2016, Lellouch et al. (2019) discovered an equatorial jet at much higher altitudes than previously known, with a surprisingly fast speed of up to ~340 m/s, but the origin of such high velocities is not yet understood. We obtained spectrally and spatially resolved ALMA observations in May 2017 to map Titan's 3D global wind field and compare our results with a reanalysis of the August 2016 data. Doppler wind velocity maps were derived in the altitude range ~300-1000 km (from the upper stratosphere to the thermosphere). At the highest, thermospheric altitudes, a 47% reduction in the equatorial zonal wind speed was measured over the 9-month period (corresponding to L_s = 82-90 degrees on Titan). This is interpreted as due to a dramatic slowing and loss of confinement (broadening) of the recently-discovered thermospheric equatorial jet, as a result of dynamical instability. These unexpectedly-rapid changes in the upper-atmospheric dynamics are consistent with strong variability of the jet's primary driving mechanism.

Through the use of a High Band Antenna system, the Low Frequency Array Two-Meter Sky Survey (LoTSS) is an attempt to complete a high-resolution survey of the northern celestial sky. To date, thousands of radio sources have been classified by LOFAR with most of them consisting of active galactic nuclei (AGN). The strong AGN emissions detected by LoTSS, it is thought, are powered by supermassive black holes (SMBHs) at the center of galaxies. During an analysis of 1500 images of these AGNs, we identified 10 radio sources with radial spokes emitted from an unknown source. According to LOFAR, the radial spokes are artifacts due to calibration errors with no origin, and therefore they cannot be associated with an optical source. Our preliminary hypothesis for the artifacts was that they were produced by ionized jets emitted from quasi-stellar objects (QSOs). Specifically, that strong emissions from Type 1 broadline (BL) quasars were directed in the line of sight of the observer (i.e., LOFAR) and as a result produced the image artifacts. To test our hypotheses, we cross-referenced the ra and dec coordinates of the artifacts with the galactic coordinates indexed in the Sloan Digital Sky Survey (SDSS) and confirmed that the artifacts were associated with BL QSOs. Further analysis of the QSOs, moreover, demonstrated they exhibited prominent broad emission lines such as CIII and MgII, which are characteristic of Type 1 BL quasars. It is our interpretation, therefore, that the radial spokes characterized as artifacts by LOFAR were produced by the emission of Type 1 BL quasars in the line of sight of the radio telescope.

We present a new set of index-based measurements of [$\alpha$/Fe] for a sample of 2093 galaxies in the SAMI Galaxy Survey. Following earlier work, we fit a global relation between [$\alpha$/Fe] and the galaxy velocity dispersion $\sigma$ for red sequence galaxies, [$\alpha$/Fe]=(0.378$\pm$0.009)log($\sigma$/100)+(0.155$\pm$0.003). We observe a correlation between the residuals and the local environmental surface density, whereas no such relation exists for blue cloud galaxies. Returning to the full sample, we find that galaxies in high-density environments are $\alpha$-enhanced by up to 0.06 dex at galaxy velocity dispersions $\sigma$<100 km/s, compared to their counterparts in low-density environments. This $\alpha$-enhancement includes a dependence on morphology, with an offset of 0.057$\pm$0.014 dex for ellipticals, and decreasing along the Hubble sequence towards spirals, with an offset of 0.019$\pm$0.014 dex. Conversely, for galaxies with $\sigma$>100 km/s in low-density environments, the [$\alpha$/Fe]-$\sigma$ relation is consistent across all morphological types earlier than Sc. At low galaxy velocity dispersion and controlling for morphology, we therefore estimate that star formation in galaxies in high-density environments is truncated $\sim$1 Gyr earlier, compared to those in low-density environments. At the highest velocity dispersions, $\sigma$>200 km/s, we find no difference in the [$\alpha$/Fe] ratio of galaxies earlier than Sc. Hence, we infer that the integrated star-formation timescales cannot differ substantially between high-$\sigma$ galaxies across varied environments, supporting the relative dominance of mass-based quenching mechanisms at the highest mass scales.

Glitches are the observational manifestations of superfluidity inside neutron stars. The aim of this paper is to describe an automated glitch detection pipeline, which can alert the observers on possible real-time detection of rotational glitches in pulsars. Post alert, the pulsars can be monitored at a higher cadence to measure the post-glitch recovery phase. Two algorithms namely, Median Absolute Deviation (MAD) and polynomial regression have been explored to detect glitches in real time. The pipeline has been optimized with the help of simulated timing residuals for both the algorithms. Based on the simulations, we conclude that the polynomial regression algorithm is significantly more effective for real time glitch detection. The pipeline has been tested on a few published glitches. This pipeline is presently implemented at the Ooty Radio Telescope. In the era of upcoming large telescopes like SKA, several hundreds of pulsars will be observed regularly and such a tool will be useful for both real-time detection as well as optimal utilization of observation time for such glitching pulsars.

Different works have recently found an increase of the average X-ray-to-radio luminosity ratio with redshift in the blazar population. We evaluate here whether the inverse Compton interaction between the relativistic electrons within the jet and the photons of the Cosmic Microwave Background (IC/CMB) can explain this trend. Moreover, we test whether the IC/CMB model can also be at the origin of the different space density evolutions found in X-ray and radio selected blazar samples. By considering the best statistically complete samples of blazars selected in the radio or in the X-ray band and covering a large range of redshift (0.5$\lesssim$ z$\lesssim$ 5.5), we evaluate the expected impact of the CMB on the observed X-ray emission on each sample and then we compare these predictions with the observations. We find that this model can satisfactorily explain both the observed trend of the X-ray-to-radio luminosity ratios with redshift and the different cosmological evolutions derived from the radio and X-ray band. Finally, we discuss how currently on-going X-ray missions, like eROSITA, could help to further constrain the observed evolution at even higher redshifts (up to z$\sim$6-7).

The presence of stable, compact circumbinary discs of gas and dust around post-asymptotic giant branch (post-AGB) binary systems has been well established. We focus on one such system: IRAS 08544-4431. We present an interferometric multi-wavelength analysis of the circumstellar environment of IRAS 08544-4431. The aim is to constrain different contributions to the total flux in the H, K, L, and N-bands in the radial direction. The data from VLTI/PIONIER, VLTI/GRAVITY, and VLTI/MATISSE range from the near-infrared, where the post-AGB star dominates, to the mid-infrared, where the disc dominates. We fitted two geometric models to the visibility data to reproduce the circumbinary disc: a ring with a Gaussian width and a flat disc model with a temperature gradient. The flux contributions from the disc, the primary star (modelled as a point-source), and an over-resolved component are recovered along with the radial size of the emission, the temperature of the disc as a function of radius, and the spectral dependencies of the different components. The trends of all visibility data were well reproduced with the geometric models. The near-infrared data were best fitted with a Gaussian ring model while the mid-infrared data favoured a temperature gradient model. This implies that a vertical structure is present at the disc inner rim, which we attribute to a rounded puffed-up inner rim. The N-to-K size ratio is 2.8, referring to a continuous flat source, analogues to young stellar objects. By combining optical interferometric instruments operating at different wavelengths we can resolve the complex structure of circumstellar discs and study the wavelength-dependent opacity profile. A detailed radial, vertical, and azimuthal structural analysis awaits a radiative transfer treatment in 3D to capture all non-radial complexity.

We extend the local stellar galaxy-(sub)halo connection to the atomic hydrogen (HI) component by seeding semi-empirically galaxies into a large N-body dark matter (DM) simulation. The main input to construct the mock galaxy catalog are: our constrained stellar mass-to-(sub)halo circular velocity ($M_{\ast}$-$V_{\rm DM}$) relation, assuming a scatter independent of any galaxy property, and the empirical $M_{\rm HI}$ conditional probability distributions given $M_{\ast}$ for central and satellite galaxies. We find that the $\langle\log M_{\rm HI}\rangle-\log M_{\rm DM}$ relation is not a monotonic increasing function. It increases with mass up to $M_{\rm DM}\sim 10^{12}$ $M_{\odot}$, attaining a maximum of $\langle\log(M_{\rm HI}/M_{\odot})\rangle \sim 9.2$, and at higher (sub)halo masses, $\langle\log(M_{\rm HI})\rangle$ decreases slightly with $M_{\rm DM}$. The scatter around it is also large and mass dependent. The bivariate $M_{\rm HI}$ and $M_{\rm DM}$ distribution is broad and bimodal, specially at $M_{\rm DM}\gtrsim 10^{12}$ $M_\odot$, which is inherited from the input $M_{\rm HI}$ conditional distributions. We also report the total (central+satellites) HI gas mass within halos, $\langle M^{\rm tot}_{\rm HI}(M_{\rm DM})\rangle$, as a function of $M_{\rm DM}$. The mean $\langle\log M^{\rm tot}_{\rm HI}(M_{\rm DM})\rangle$ relation is an increasing monotonic function. The galaxy spatial clustering increases weakly as the $M_{\rm HI}$ threshold increases. Our HI mock galaxies cluster more in comparison to the blind HI ALFALFA survey but we show that it is mainly due to the selection effects. We discuss the implications of our results in the light of predictions from semi-analytical models and hydrodynamics simulations of galaxy evolution.

In this reply, we address the comment [arXiv:2105.14908] to our recent paper [arXiv:2105.09328], where we argued that the Thakurta metric does not describe cosmological black holes. We clarify that the mass growth of Thakurta black holes is due to an influx of energy (i.e. accretion), which, by definition, is not a feature of geometry. The conclusions of [arXiv:2105.09328] are independent of the interpretation of this energy flux. We show that the average energy density of primordial Thakurta black holes scales as $a^{-2}$ and requires an unrealistic and fine-tuned energy transfer from a smooth dark matter component to the primordial black hole sector.

The photometric analysis of sample Am stars is carried out to determine the stellar characteristics and to constrain the stellar dynamics. The spectroscopic analysis of the studied Am stars confirms their general characteristics of Am stars. The available data on elemental abundances for HD 113878 and HD 118660 have shown different characteristics during different epochs of observations. The basic stellar parameters (mass, luminosity, radius, life time, distance, proper-motion, etc.) are also determined to identify the stellar habitat zones for earth like exoplanet. Such information is important to identify suitable planets for human settlement in the near future. In this connection, the tidal radius and boundaries of the habitable zone of each star have been computed to support the search of an extra-terrestrial life around them. Asteroseismic mass scale test shows greater stellar masses comparable to the solar mass.

Inflationary model driven by a scalar field whose potential has a step in the second derivative with respect to the field is considered. For the best fit potential parameter values, the 3-point function and the non-Gaussianity associated with the featured model is calculated. We study the shape and scale dependence of the 3-point function. The distinctive feature of this model is its characteristic ringing behaviour of $f_{NL}$. We can see that the oscillations in $f_{NL}$ in this model last for a much longer range of k values, as compared to the previously studied models. In that sense, this model is potentially distinguishable from models with other features in the potential.

We show how the generation of right-handed neutrino masses in Majoron models may be associated with a first-order phase transition and accompanied by the production of a stochastic background of gravitational waves (GWs). We explore different energy scales with only renormalizable operators in the effective potential. If the phase transition occurs above the electroweak scale, the signal can be tested by future interferometers. We consider two possible energy scales for phase transitions below the electroweak scale. If the phase transition occurs at a GeV, the signal can be tested at LISA and provide a complementary cosmological probe to right-handed neutrino searches at the FASER detector. If the phase transition occurs below 100 keV, we find that the peak of the GW spectrum is two or more orders of magnitude below the putative NANOGrav GW signal at low frequencies, but well within reach of the SKA and THEIA experiments. We show how searches of very low frequency GWs are motivated by solutions to the Hubble tension in which ordinary neutrinos interact with the dark sector. We also present general calculations of the phase transition and Euclidean action that apply beyond Majoron models.

The possibility that rotating black holes could be natural particle accelerators has been subject of intense debate. While it appears that for extremal Kerr black holes arbitrarily high center of mass energies could be achieved, several works pointed out that both theoretical as well as astrophysical arguments would severely dampen the attainable energies. In this work we study particle collisions near Kerr--Newman black holes, by reviewing and extending previously proposed scenarios. Most importantly, we implement the hoop conjecture for all cases and we discuss the astrophysical relevance of these collisional Penrose processes. The outcome of this investigation is that scenarios involving near-horizon target particles are in principle able to attain, sub-Planckian, but still ultra high, center of mass energies of the order of $10^{21}-10^{23}$ eV. Thus, these target particle collisional Penrose processes could contribute to the observed spectrum of ultra high-energy cosmic rays, even if the hoop conjecture is taken into account, and as such deserve further scrutiny in realistic settings.

We perform a systematic study of the electric and magnetic dipole moments of dark matter (DM) that are induced at the one-loop level when DM experiences four-fermion interactions with Standard Model (SM) charged fermions. Related to their loop nature these moments can largely depend on the UV completion at the origin of the four-fermion operators. We illustrate this property by considering explicitly two simple ways to generate these operators, from $t$- or $s$-channel tree-level exchange. Fixing the strength of these interactions from the DM relic density constraint, we obtain in particular a magnetic moment that, depending on the interaction considered, lies typically between $10^{-20}$ to $10^{-23}$ ecm or identically vanishes. These non-vanishing values induce, via photon exchange, DM-nucleus scattering cross sections that could be probed by current or near future direct detection experiments.

In the framework of the Einstein-Maxwell-axion-aether theory we establish the model, the Lagrangian of which contains the sin-type generalization of the term describing the axion-photon coupling, and the axionically induced cosine-type modification of the term attributed to the dilaton-photon interactions. The extension of the axion-dilaton-aether electrodynamics is inspired by the Jackson's idea concerning the internal symmetry of the equations of electromagnetism. The application of the extended theory to the anisotropic homogeneous cosmological model of the Bianchi-I type is considered. The exact solutions to the model evolutionary equations are obtained for the case, when the axionic dark matter is in the state of equilibrium, which is characterized by vanishing potential of the pseudoscalar field and its first derivative. The state of the axion-photon system, which is of a new type and is indicated as a dynamic equilibrium, is studied in the framework of electrodynamics with axionic non-linearity. We show that the nonlinear axion-photon interactions can mimic the dilaton-photon coupling. We discuss the stability of the model with respect to homogeneous fluctuations of the axion field.

The black holes that have been detected via gravitational waves (GW) can have either astrophysical or primordial origin. Some GW events show significant spin for one of the components and have been assumed to be astrophysical, since primordial black holes are generated with very low spins. However, it is worth studying if they can increase their spin throughout the evolution of the universe. Possible mechanisms that have already been explored are multiple black hole mergers and gas accretion. We propose here a new mechanism that can occur in dense clusters of black holes: the spin-up of primordial black holes when they are involved in close hyperbolic encounters. We explore this effect numerically with the Einstein Toolkit for different initial conditions, including variable mass ratios. For equal masses, there is a maximum spin that can be induced on the black holes, $\chi = a/m \leq 0.2$. We find however that for large mass ratios one can attain spins up to $\chi \simeq 0.8$, where the highest spin is induced on the most massive black hole. For small induced spins we provide simple analytical expressions that depend on the relative velocity and impact parameter.

Motivated by the SU(2)$_{\rm CMB}$ modification of the cosmological model $\Lambda$CDM, we consider isolated fuzzy-dark-matter lumps, made of ultralight axion particles whose masses arise due to distinct SU(2) Yang-Mills scales and the Planck mass $M_P$. In contrast to SU(2)$_{\rm CMB}$, these Yang-Mills theories are in confining phases (zero temperature) throughout most of the Universe's history and associate with the three lepton flavours of the Standard Model of particle physics. As the Universe expands, axionic fuzzy dark matter comprises a three-component fluid which undergoes certain depercolation transitions when dark energy (a global axion condensate) is converted into dark matter. We extract the lightest axion mass $m_{a,e}= 0.675\times 10^{-23}\,$eV from well motivated model fits to observed rotation curves in low-surface-brightness galaxies (SPARC catalogue). Since the virial mass of an isolated lump solely depends on $M_P$ and the associated Yang-Mills scale the properties of an e-lump predict those of $\mu$- and $\tau$-lumps. As a result, a typical e-lump virial mass $\sim 6.3\times 10^{10}\,M_\odot$ suggests that massive compact objects in galactic centers such as Sagittarius A$^*$ in the Milky Way are (merged) $\mu$- and $\tau$-lumps. In addition, $\tau$-lumps may constitute globular clusters. SU(2)$_{\rm CMB}$ is always thermalised, and its axion condensate never has depercolated. If the axial anomaly indeed would link leptons with dark matter and the CMB with dark energy then this would demystify the dark Universe through a firmly established feature of particle physics.

Cosmic ray muons, that reach the earth's surface, provide a natural source of radiation that is used for radiography. In this paper, we show that radiography using cosmic radiation background provides a method that can be used to monitor bulk aspects of human anatomy. We describe a method that can be used to measure changes in patients as a function of time by radiographing them using cosmic-ray muons. This could provide hourly readouts of parameters such as lung density with sufficient sensitivity to detect time changes in inflammation of the lungs in, e.g., Covid patients.

We review the derivation of the pulsations equations for spherically symmetric boson stars and then make a thorough study of the radial oscillation frequencies for the fundamental and first excited modes. We do this for self-interacting boson stars and consider a range of values for the self-coupling constant. We also numerically evolve boson stars and Fourier transform the dynamic solutions. The Fourier transform gives an independent computation of the radial oscillation frequencies and allows us to verify our results obtained from the pulsation equations. We find excellent agreement between the two methods.

Aims: We present the first measurements of the solar-wind angular-momentum (AM) flux recorded by the Solar Orbiter spacecraft. Our aim is the validation of these measurements to support future studies of the Sun's AM loss. Methods: We combine 60-minute averages of the proton bulk moments and the magnetic field measured by the Solar Wind Analyser (SWA) and the magnetometer (MAG) onboard Solar Orbiter. We calculate the AM flux per solid-angle element using data from the first orbit of the mission's cruise phase during 2020. We separate the contributions from protons and from magnetic stresses to the total AM flux. Results: The AM flux varies significantly over time. The particle contribution typically dominates over the magnetic-field contribution during our measurement interval. The total AM flux shows the largest variation and is typically anti-correlated with the radial solar-wind speed. We identify a compression region, potentially associated with a co-rotating interaction region or a coronal mass ejection, that leads to a significant localised increase in the AM flux, yet without a significant increase in the AM per unit mass. We repeat our analysis using the density estimate from the Radio and Plasma Waves (RPW) instrument. Using this independent method, we find a decrease in the peaks of positive AM flux but otherwise consistent results. Conclusions: Our results largely agree with previous measurements of the solar-wind AM flux in terms of amplitude, variability, and dependence on radial solar-wind bulk speed. Our analysis highlights the potential for future, more detailed, studies of the solar wind's AM and its other large-scale properties with data from Solar Orbiter. We emphasise the need to study the radial evolution and latitudinal dependence of the AM flux in combination with data from Parker Solar Probe and assets at heliocentric distances of 1 au and beyond.

This chapter sets the stage for the rest of the book by exploring the role of intuition as a tool to deepen understanding in Einsteinian physics. Drawing on examples from the history of general relativity, we argue that the development of physical intuition is a crucial goal in physics education in parallel with any mathematical development of a physics subject. This chapter is for readers who wish to learn how expert physicists think conceptually about their subjects to understand them plus readers who wish to see how we can introduce Einsteinian physics to students by developing their intuition as well as teaching them the mathematics.

In this work we investigate neutron stars (NS) in $f(\mathtt{R,L_m})$ theory of gravity for the case $f(\mathtt{R,L_m}) = \mathtt{R} + \mathtt{L_m} + \sigma\mathtt{R}\mathtt{L_m}$, where $\mathtt{R}$ is the Ricci scalar and $\mathtt{L_m}$ the Lagrangian matter density. In the term $\sigma\mathtt{R}\mathtt{L_m}$, $\sigma$ represents the coupling between the gravitational and particles fields. For the first time the hydrostatic equilibrium equations in the theory are solved considering realistic equations of state and NS masses and radii obtained are subject to joint constrains from massive pulsars, the gravitational wave event GW170817 and from the PSR J0030+0451 mass-radius from NASA's Neutron Star Interior Composition Explorer (${\it NICER}$) data. We show that in this theory of gravity, the mass-radius results can accommodate massive pulsars, while the general theory of relativity can hardly do it. The theory also can explain the observed NS within the radius region constrained by the GW170817 and PSR J0030+0451 observations for masses around $1.4~M_{\odot}$.

We propose a new portal coupling to dark matter by taking advantage of the nonminimally coupled portal sector to the Ricci scalar. Such a portal sector conformally induces couplings to the trace of the energy-momentum tensor of matters including highly secluded dark matter particles. The portal coupling is so feeble that dark matter is produced by freeze-in processes of scatterings and/or the decay of the mediator. We consider two concrete realizations of the portal: conformally induced Higgs portal and conformally induced mediator portal. The former case is compatible with the Higgs inflation, while the latter case can be tested by dark matter direct detection experiments.

The core-envelope models presented in {Ref1}; {Ref2}, corresponding to the values of compactness parameter, $u \equiv M/a$ = 0.30 and 0.25 (mass to size ratio in geometrized units) have been studied under slow rotation. It is seen that these models are capable of explaining all the observational values of glitch healing parameter, $G_h = I_{\rm core}/I_{\rm total} < 0.55$ {Ref3} for the Vela pulsar. The models yield the maximum values of mass, $M$, surface redshift, $z_a$, and the moment of inertia, $I_{\rm Vela}$ for the Vela pulsar in the range $M = 3.079M_\odot - 2.263M_\odot$; $z_a = 0.581 - 0.414$ and $I_{\rm Vela,45} =6.9 - 3.5$ (where $I_{45}=I/10^{45}\rm g{cm}^2$) respectively for the values of $u = $ 0.30 and 0.25 and for an assigned value of the surface density, $E_a = 2\times 10^{14}\rm g{cm}^{-3}$ {Ref4}. The values of masses lower than the above mentioned values ( so called the realistic mass range, $M = 1.4\pm0.2 M_\odot$, in the literature) but significantly higher than that of the unrealistic mass range $M \leq 0.5M_\odot$ (obtained for the Vela pulsar in the literature on the basis of parametrized neutron star (NS) models based on equations of state (EOSs) of dense nuclear matter {Ref3}) and other parameters may be obtained likewise for the above mentioned range of the values of $G_h$ corresponding to the values of $u < 0.25$. The models are found to be causally consistent, gravitational bound and pulsationally stable. The upper bound on neutron star (NS) mass obtained in this study which is applicable for the Vela pulsar, in fact, corresponds to the mean value of the upper bound on NS mass obtained in the classical result {Ref5} and that obtained on the basis of modern EOSs for neutron star matter {Ref6} and is in a good agreement with the most recent theoretical estimate {Ref7}.

This work presents the construction of a novel spherical wavelet basis designed for incomplete spherical datasets, i.e. datasets which are missing in a particular region of the sphere. The eigenfunctions of the Slepian spatial-spectral concentration problem (the Slepian functions) are a set of orthogonal basis functions which exist within a defined region. Slepian functions allow one to compute a convolution on the incomplete sphere by leveraging the recently proposed sifting convolution and extending it to any set of basis functions. Through a tiling of the Slepian harmonic line one may construct scale-discretised wavelets. An illustration is presented based on an example region on the sphere defined by the topographic map of the Earth. The Slepian wavelets and corresponding wavelet coefficients are constructed from this region, and are used in a straightforward denoising example.