Papers for Friday, Aug 04 2023

We present a comparison of the Milky Way's star formation rate (SFR) surface density ($\Sigma_{\rm SFR}$) obtained with two independent state-of-the-art observational methods. The first method infers $\Sigma_{\rm SFR}$ from the observations of the dust thermal emission from interstellar dust grains in far infrared wavelengths registered in the Herschel Infrared Galactic Plane Survey (Hi-GAL), as presented in Elia et al. (2022). The second method obtains $\Sigma_{\rm SFR}$ by modeling the current population of O-, B-, and A-type stars in a 6 kpc $\times$ 6 kpc area around the Sun, as presented in Zari et al. (2023). We found an agreement between the two methods within a factor of two for the mean SFRs and the SFR surface density profiles. Given the broad differences between the observational techniques and the independent assumptions in the methods to compute the SFRs, this agreement constitutes a significant advance in our understanding of the star formation of our Galaxy and implies that the local SFR has been roughly constant over the past 10 Myr.

The dominant formation mechanism of brown dwarfs and planetary mass objects in star-forming regions is presently uncertain. Do they form like stars, via the collapse and fragmentation of cores in Giant Molecular clouds, or do they form like planets in the discs around stars and are ejected via dynamical interactions? In this paper, we quantify the spatial distribution of substellar objects in NGC1333, in particular focusing on planetary-mass objects that have been the target of recent deep imaging observations. We find that these objects have a spatial distribution that is indistinguishable from the stars, and more massive brown dwarfs. We also analyse N-body simulations and find that a population of ejected planets would have a significantly different spatial and kinematic distribution to stars, and brown dwarfs that also formed through gravitational collapse and fragmentation. We therefore conclude that the low-mass substellar objects in NGC1333 formed more like stars than planets, although we predict that a population of hitherto undetected ejected planetary mass objects may be lurking in this, and other star-forming regions.

Nitrogen hydrides such as NH$_3$ and N$_2$H$^+$ are widely used by Galactic observers to trace the cold dense regions of the interstellar medium. In external galaxies, because of limited sensitivity, HCN has become the most common tracer of dense gas over large parts of galaxies. We provide the first systematic measurements of N$_2$H$^+$(1-0) across different environments of an external spiral galaxy, NGC6946. We find a strong correlation ($r>0.98,p<0.01$) between the HCN(1-0) and N$_2$H$^+$(1-0) intensities across the inner $\sim8\mathrm{kpc}$ of the galaxy, at kiloparsec scales. This correlation is equally strong between the ratios N$_2$H$^+$(1-0)/CO(1-0) and HCN(1-0)/CO(1-0), tracers of dense gas fractions ($f_\mathrm{dense}$). We measure an average intensity ratio of N$_2$H$^+$(1-0)/HCN(1-0)$=0.15\pm0.02$ over our set of five IRAM-30m pointings. These trends are further supported by existing measurements for Galactic and extragalactic sources. This narrow distribution in the average ratio suggests that the observed systematic trends found in kiloparsec-scale extragalactic studies of $f_\mathrm{dense}$ and the efficiency of dense gas (SFE$_\mathrm{dense}$) would not change if we employed N$_2$H$^+$(1-0) as a more direct tracer of dense gas. At kiloparsec scales our results indicate that the HCN(1-0) emission can be used to predict the expected N$_2$H$^+$(1-0) over those regions. Our results suggest that, even if HCN(1-0) and N$_2$H$^+$(1-0) trace different density regimes within molecular clouds, subcloud differences average out at kiloparsec scales, yielding the two tracers proportional to each other.

Many post-processing algorithms have been developed in order to better separate the signal of a companion from the bright light of the host star, but the effect of such algorithms on the shape of exoplanet spectra extracted from integral field spectrograph data is poorly understood. The resulting spectra are affected by noise that is correlated in wavelength space due to both optical and data processing effects. Within the framework of Bayesian atmospheric retrievals, we aim to understand how these correlations and other systematic effects impact the inferred physical parameters. We consider three algorithms (KLIP, PynPoint and ANDROMEDA), optimizing the choice of algorithmic parameters using a series of injection tests into archival SPHERE and GPI data of the HR 8799 system. The wavelength-dependent covariance matrix is calculated to provide a measure of instrumental and algorithmic systematics. We perform atmospheric retrievals using petitRADTRANS on optimally extracted spectra to measure how these data processing systematics influence the retrieved parameter distributions. The choice of data processing algorithm and parameters significantly impact the accuracy of retrieval results, with the mean posterior parameter bias ranging from 1 to 3 $\sigma$ from the true input parameters. Including the full covariance matrix in the likelihood improves the accuracy of inferred parameters, and cannot be accounted for using ad hoc scaling parameters in the retrieval framework. Using the Bayesian information criterion and other statistical measures as a heuristic goodness-of-fit metrics, the retrievals including the full covariance matrix are favoured when compared to using only the diagonal elements.

We use the spectro-photometric information of ~219 million stars from Gaia's DR3 to calculate synthetic, narrow-band, metallicity-sensitive CaHK magnitudes that mimic the observations of the Pristine survey, a survey of photometric metallicities of Milky Way stars that has been mapping more than 6,500 deg^2 of the northern sky with the CFHT since 2015. These synthetic magnitudes are used for an absolute re-calibration of the deeper Pristine photometry and, combined with broadband Gaia information, synthetic and Pristine CaHK magnitudes are used to estimate photometric metallicities over the whole sky. The resulting metallicity catalogue is accurate down to [Fe/H]~-3.5 and is particularly suited for the exploration of the metal-poor Milky Way ([Fe/H]<-1.0). We make available here the catalogue of synthetic CaHK_syn magnitudes for all stars with BP/RP information in Gaia DR3, as well as an associated catalogue of more than ~30 million photometric metallicities for high S/N FGK stars. This paper further provides the first public DR of the Pristine catalogue in the form of higher quality recalibrated Pristine CaHK magnitudes and photometric metallicities for all stars in common with the BP/RP information in Gaia DR3. We demonstrate that, when available, the much deeper Pristine data greatly enhances the quality of the derived metallicities, in particular at the faint end of the catalogue (G_BP>16). Combined, both catalogues include more than 2 million metal-poor star candidates as well as more than 200,000 and ~8,000 very and extremely metal-poor candidates. Finally, we show that these metallicity catalogues can be used efficiently, among other applications, for Galactic archaeology, to hunt for the most metal-poor stars, and to study how the structure of the Milky Way varies with metallicity, from the flat distribution of disk stars to the spheroid-shaped metal-poor halo. (Shortened)

The connection between X-ray weakness and powerful X-ray outflows is both expected in a scenario where outflows are connected with radiation-driven winds, and observed in several sources, both in the local Universe and at high redshift. Here I present the first results of a new study of this possible connection based on a search for SDSS quasars with weak X-ray emission in serendipitous XMM-Newton observations. The selected objects have a "normal" optical/UV blue continuum, but a flat and extraordinarily weak X-ray spectrum. The availability of rest-frame optical/UV spectra allows to check for the signature of outflows in the absorption lines and/or in the profiles of the emission lines. This method could reveal the presence of a population of so-far overlooked outflowing quasars and confirm the connection between winds and X-ray weakness in quasars.

Minimising the impact of stellar variability in Radial Velocity (RV) measurements is a critical challenge in achieving the 10 cm s$^{-1}$ precision needed to hunt for Earth twins. Since 2012, a dedicated programme has been underway with HARPS-N, to conduct a blind RV Rocky Planets Search (RPS) around bright stars in the Northern Hemisphere. Here we describe the results of a comprehensive search for planetary systems in two RPS targets, HD 166620 and HD 144579. Using wavelength-domain line-profile decorrelation vectors to mitigate the stellar activity and performing a deep search for planetary reflex motions using a trans-dimensional nested sampler, we found no significant planetary signals in the data sets of either of the stars. We validated the results via data-splitting and injection recovery tests. Additionally, we obtained the 95th percentile detection limits on the HARPS-N RVs. We found that the likelihood of finding a low-mass planet increases noticeably across a wide period range when the inherent stellar variability is corrected for using scalpels U-vectors. We are able to detect planet signals with $M\sin i \leq 1$ M$_\oplus$ for orbital periods shorter than 10 days. We demonstrate that with our decorrelation technique, we are able to detect signals as low as 54 cm s$^{-1}$, which brings us closer to the calibration limit of 50 cm s$^{-1}$ demonstrated by HARPS-N. Therefore, we show that we can push down towards the RV precision required to find Earth analogues using high-precision radial velocity data with novel data-analysis techniques.

This tutorial is an introduction to High-Contrast Imaging, a technique that enables astronomers to isolate light from faint planets and/or circumstellar disks that would otherwise be lost amidst the light of their host stars. Although technically challenging, high-contrast imaging allows for direct characterization of the properties of detected circumstellar sources. The intent of the article is to provide newcomers to the field a general overview of the terminology, observational considerations, data reduction strategies, and analysis techniques high-contrast imagers employ to identify, vet, and characterize planet and disk candidates.

We reconstruct the morphology and kinematics of a series of small transients that erupt from the Sun on 2021 April 24 using observations primarily from Parker Solar Probe (PSP). These sequential small coronal mass ejections (CMEs) may be the product of continuous reconnection at a current sheet, a macroscopic example of the more microscopic reconnection activity that has been proposed to accelerate the solar wind more generally. These particular CMEs are of interest because they are the first CMEs to hit PSP and be simultaneously imaged by it, using the Wide-field Imager for Solar Probe (WISPR) instrument. Based on imaging from WISPR and STEREO-A, we identify and model six discrete transients, and determine that it is the second of them (CME2) that first hits PSP, although PSP later more obliquely encounters the third transient as well. Signatures of these encounters are seen in the PSP in situ data. Within these data, we identify six candidate magnetic flux ropes (MFRs), all but one of which are associated with the second transient. The five CME2 MFRs have orientations roughly consistent with PSP encountering the right sides of roughly E-W oriented MFRs, which are sloping back towards the Sun.

Ultraviolet spectra were taken of 25 Detached Eclipsing Binaries (DEBs) with spectral types O, B, and early A with the International Ultraviolet Explorer (IUE) satellite in the 1150 to 1900 ${\AA}$ region. The spectra were compared with BOSZ model atmospheres (Bohlin, et al. 2017). The composite spectra of the DEBs were modeled by a combination of models representing the hot and cool components, and the temperatures of the hottest components of the systems were determined. From these temperatures a direct Mass-Temperature relation was obtained for stars close to the main sequence with solar metallicity for B and early A stars: log M/Msun = -5.90 $\pm$ 0.27 + (1.56 $\pm$ 0.07) x log T This relation allows a mass to be inferred for comparable stars from an ultraviolet spectrum. The five chemically peculiar Am stars in the sample have larger radii than normal A stars of the same mass.

We return to interpreting the historical SN~1987A neutrino data from a modern perspective. To this end, we construct a suite of spherically symmetric supernova models with the Prometheus-Vertex code, using four different equations of state and five choices of final baryonic neutron-star (NS) mass in the 1.36-1.93 M$_\odot$ range. Our models include muons and proto-neutron star (PNS) convection by a mixing-length approximation. The time-integrated signals of our 1.44 M$_\odot$ models agree reasonably well with the combined data of the four relevant experiments, IMB, Kam-II, BUST, and LSD, but the high-threshold IMB detector alone favors a NS mass of 1.7-1.8 M$_\odot$, whereas Kam-II alone prefers a mass around 1.4 M$_\odot$. The cumulative energy distributions in these two detectors are well matched by models for such NS masses, and the previous tension between predicted mean neutrino energies and the combined measurements is gone, with and without flavor swap. Generally, our predicted signals do not strongly depend on assumptions about flavor mixing, because the PNS flux spectra depend only weakly on antineutrino flavor. While our models show compatibility with the events detected during the first seconds, PNS convection and nucleon correlations in the neutrino opacities lead to short PNS cooling times of 5-9 s, in conflict with the late event bunches in Kam-II and BUST after 8-9 s, which are also difficult to explain by background. Speculative interpretations include the onset of fallback of transiently ejected material onto the NS, a late phase transition in the nuclear medium, e.g., from hadronic to quark matter, or other effects that add to the standard PNS cooling emission and either stretch the signal or provide a late source of energy. More research, including systematic 3D simulations, is needed to assess these open issues.

(Abridged) Molecular lines are commonly detected towards protostellar sources. However, to get a better understanding of the chemistry of these sources we need unbiased molecular surveys over a wide frequency range for as many sources as possible to shed light on the origin of this chemistry, particularly any influence from the external environment. We present results from the PILS-Cygnus survey of ten intermediate- to high-mass protostellar sources in the nearby Cygnus-X complex, through high angular resolution interferometric observations over a wide frequency range. Using the Submillimeter Array (SMA), a spectral line survey of ten sources was performed in the frequency range 329-361 GHz, with an angular resolution of $\sim$1\farcs5, ($\sim$2000 AU, source distance of 1.3 kpc). Spectral modelling was performed to identify molecular emission and determine column densities and excitation temperatures for each source. We detect CH$_3$OH towards nine of the ten sources, CH$_3$OCH$_3$ and CH$_3$OCHO towards three sources, and CH$_3$CN towards four sources. Towards five sources the chemistry is spatially differentiated (different species peak at different positions and are offset from the peak continuum emission). The chemical properties of each source do not correlate with their position in the Cygnus-X complex, nor do the distance or direction to the nearest OB associations. However, the five sources located in the DR21 filament do appear to show less line emission compared to the five sources outside the filament. This work shows how important wide frequency coverage observations are combined with high angular resolution observations for studying the protostellar environment. Based on the ten sources observed here, the external environment appears to only play a minor role in setting the chemical environment on these small scales ($<$ 2000 AU).

We present an analysis of the first connection mosaic made by the SPICE instrument on board of the ESA / NASA Solar Orbiter mission on March 2$^{nd}$, 2022. The data will be used to map coronal composition that will be compared with in-situ measurements taken by SWA/HIS to establish the coronal origin of the solar wind plasma observed at Solar Orbiter. The SPICE spectral lines were chosen to have varying sensitivity to the First Ionization Potential (FIP) effect, and therefore the radiances of the spectral lines will vary significantly depending on whether the elemental composition is coronal or photospheric. We perform temperature diagnostics using line ratios and Emission Measure (EM) loci, and compute relative FIP biases using three different approaches (two line ratio (2LR), ratios of linear combinations of spectral lines (LCR), and differential emission measure (DEM) inversion) in order to perform composition diagnostics in the corona. We then compare the SPICE composition analysis and EUI data of the potential solar wind sources regions to the SWA / HIS data products. Radiance maps are extracted from SPICE spectral data cubes, with values matching previous observations. We find isothermal plasma of around LogT=5.8 for the active region loops targeted, and that higher FIP-bias values are present at the footpoints of the coronal loops associated with two active regions.

The bar structure in disk galaxies models is formed by different families of orbits; however, it is not clear how these families of orbits support the bar throughout its secular evolution. Here, we analyze the orbital structure on three stellar disk N-body models embedded in a live dark matter halo. During the evolution of the models, disks naturally form a bar that buckles out of the galactic plane at different ages of the galaxy evolution generating boxy, X, peanut, and/or elongated shapes. To understand how the orbit families hold the bar structure, we evaluate the orbital evolution using the frequency analysis on phase space coordinates for all disk particles at different time intervals. We analyze the density maps morphology of the 2:1 family as the bar potential evolves. We showed that the families of orbits providing bar support exhibit variations during different stages of its evolutionary process, specifically prior to and subsequent to the buckling phase, likewise in the secular evolution of the bar. The disk-dominated model develops an internal boxy structure after the first Gyr. Afterwards, the outer part of the disk evolves into a peanut-shape, which lasts till the end of the simulation. The intermediary model develops the boxy structure only after 2 Gyr of evolution. The peanut shape appears 2 Gyr later and evolves slowly. The halo-dominated model develops the boxy structure much later, around 3 Gyr, and the peanut morphology is just incipient at the end of the simulation.

We investigate optical characteristics of flashes caused by impacting meter- to decameter-sized outer solar system objects on Jupiter and contributions of reflected light from surface clouds at visible wavelengths to estimate more accurate bulk parameters such as the luminous energy of the flash, the kinetic energy, the mass, and the size of the impact object. Based on the results of recent reflectivity studies of the Jovian surface, we develop a cloud reflection model that calculates the contribution of the reflected light relative to that directly from the flash. We compare the apparent luminous energy of the previously reported flashes with the expected cloud reflection contributions to obtain their revised bulk parameters. We found that the cloud reflection contributions can be up to 200% of the flux directly from the flash and thus can be the most significant uncertainty in the measurement of the bulk parameters. The reflection contributions strongly depend on wavelength. With our cloud reflection correction, the revised bulk parameters of the previously reported flashes are obtained. Our cloud reflection correction provides a better understanding of the properties of impacting objects on Jupiter and is crucial for ongoing detailed investigations using high-sensitivity and multi-wavelength observation systems such as PONCOTS. It will also be useful for understanding other optical transients in Jupiter's upper atmosphere, such as the recently discovered sprite-like events.

We report the discovery and characterisation of a giant transiting planet orbiting around a nearby M3.5V dwarf (d = 80.4 pc, G = 15.1 mag, K=11.2 mag, R$_\star$ = 0.354 $\pm$ 0.011 R$_\odot$, M$_\star$ = 0.3400 $\pm$ 0.0086 M$_\odot$). Using the photometric time series from the Transiting Exoplanet Survey Satellite (TESS) sectors 10, 36, 46, and 63, and near-infrared spectrophotometry from ExTrA, we measured a planetary radius of 0.766 $\pm$ 0.026 R$_J$ and an orbital period of 1.52 days. With high-resolution spectroscopy taken by the CFHT/SPIRou and ESO/ESPRESSO spectrographs, we refined the host star parameters ([Fe/H] = 0.27 $\pm$ 0.12) and measured the mass of the planet (0.2729 $\pm$ 0.0058 M$_J$). Based on these measurements, TOI- 4860 b joins the small set of massive planets found around mid-to-late M dwarfs (< 0.4 R$_\odot$), providing both an interesting challenge to planet formation theory and a favourable target for further atmospheric studies with transmission spectroscopy. We identify an additional signal in the radial velocity data that we attribute to an eccentric (e = 0.657 $\pm$ 0.089) planet candidate with an orbital period of 426.9 $\pm$ 7.4 days and a minimum mass of 1.66 $\pm$ 0.26 M$_J$.

Many relaxed cool-core clusters host diffuse radio emission on scales of hundreds of kiloparsecs: mini-haloes. However, the mechanism responsible for generating them, as well as their connection with central active galactic nuclei, is elusive and many questions related to their physical properties and origins remain unanswered. This paper presents new radio observations of the galaxy cluster Abell 1413 performed with MeerKAT (L-band; 872 to 1712 MHz) and LOFAR HBA (120 to 168 MHz) as part of a statistical and homogeneous census of mini-haloes. Abell 1413 is unique among mini-halo clusters as it is a moderately-disturbed non-cool-core cluster. Our study reveals an asymmetric mini-halo up to 584 kpc in size at 1283 MHz, twice as large as first reported at similar frequencies. The spectral index is flatter than previously reported, with an integrated value of $\alpha = -1.01 \pm 0.06$, shows significant spatial variation, and a tentative radial steepening. We studied the point-to-point X-ray/radio surface brightness correlation to investigate the thermal/non-thermal connection: our results show a strong connection between these components, with a super-linear slope of $b = 1.63 \pm 0.10$ at 1283 MHz and $b = 1.20 \pm 0.12$ at 145 MHz. We also explore the X-ray surface brightness/radio spectral index correlation, finding a slope of $b = 0.59 \pm 0.11$. Both investigations support the evidence of spectral steepening. Finally, in the context of understanding the particle acceleration mechanism, we present a simple theoretical model which demonstrates that hybrid scenarios - secondary electrons (re-)accelerated by turbulence - reproduce a super-linear correlation slope.

We present an analysis of the Lyman continuum (LyC) escape fraction based on a physically motivated model, and applied to galaxies of the Illustris TNG50 simulation in the redshift range $z=5.2-20$. Our study reveals a bimodal nature of LyC escape, which is associated either to (a) high metallicity ($10^{-3.5}<Z<10^{-2}$), low mass ($M_\star<10^7\mathrm{M}_\odot$) galaxies exhibiting extended star formation, with photons escaping primarily from the outer regions of the galactic disk ({\it ext}-mode), or (b) low metallicity ($Z<10^{-3}$), moderately massive galaxies ($M_\star<10^8\mathrm{M}_\odot$) with localized LyC escape originating from small central regions of high star formation ({\it loc}-mode). While the {\it loc}-mode is present at all redshifts under investigation, the {\it ext}-mode becomes prominent in small galaxies at later cosmic times, once sufficient metal enrichment has occurred. Building on these findings, we develop an analytical fitting formula to determine the escape fraction of galaxies based on their stellar and gas mass, as well as redshift, providing a valuable subgrid modelling tool for future studies.

A mid-infrared nulling-space interferometer is a promising way to characterize thermal light from habitable planet candidates around Sun-like stars. However, one of the main challenges for achieving this ambitious goal is a high-precision stability of the optical path difference (OPD) and amplitude over a few days for planet detection and up to a few weeks for in-depth characterization. Here we propose a new method called phase-space synthesis decomposition (PSSD) to shorten the stability requirement to minutes, significantly relaxing the technological challenges of the mission. Focusing on what exactly modulates the planet signal in the presence of the stellar leak and systematic error, PSSD prioritizes the modulation of the signals along the wavelength domain rather than baseline rotation. Modulation along the wavelength domain allows us to extract source positions in parallel to the baseline vector for each exposure. The sum of the one-dimensional data converts into two-dimensional information. Based on the reconstructed image, we construct a continuous equation and extract the spectra through the singular value decomposition (SVD) while efficiently separating them from a long-term systematic stellar leak. We performed numerical simulations to investigate the feasibility of PSSD for the LIFE mission concept. We confirm that multiple terrestrial planets in the habitable zone around a Sun-like star at 10 pc can be detected and characterized despite high levels and long durations of systematic noise. We also find that PSSD is more robust against a sparse sampling of the array rotation compared to purely rotation-based signal extraction. Using PSSD as signal extraction method significantly relaxes the technical requirements on signal stability and further increases the feasibility of the LIFE mission.

Hot magnetic stars often exhibit incoherent circularly polarized radio emission thought to arise from gyro-synchrotron emission by energetic electrons trapped in the circumstellar magnetosphere. Theoretical scalings for electron acceleration by magnetic reconnection driven by centrifugal breakout match well the empirical scalings for observed radio luminosity with both the magnetic field strength and the stellar rotation rate. This paper now examines how energetic electrons introduced near the top of closed magnetic loops are subsequently cooled by the energy loss associated with their gyro-synchrotron radio emission. For sample assumed distributions for energetic electron deposition about the loop apex, we derive the spatial distribution of the radiated energy from such "gyro-cooling". For sub-relativistic electrons, we show explicitly that this is independent of the input energy, but also find that even extensions to the relativistic regime still yield a quite similar spatial distribution. However, cooling by coulomb collisions with even a modest ambient density of thermal electrons can effectively quench the emission from sub-relativistic electrons, indicating that the observed radio emission likely stems from relativistic electrons that are less affected by such collisional cooling. The overall results form an initial basis for computing radio emission spectra in future models that account for such cooling and multimode excitation about the fundamental gyro-frequency. Though motivated in the context of hot-stars, the basic results here could also be applied to gyro-emission in any dipole magnetospheres, including those of ultra-cool dwarfs and even (exo)-planets.

Announcing the final release, v8, of the Milliquas (Million Quasars) quasar catalogue which presents all published quasars to 30 June 2023, including quasars from the first releases of the Dark Energy Spectroscopic Instrument (DESI) and the SDSS-DR18 Black Hole Mapper. Its totals are 907,144 type-I QSOs/AGN and 66,026 high-confidence (~99% likelihood) radio/X-ray associated quasar candidates. Type-II and Bl Lac type objects are also included, bringing the total count to 1,021,800. Gaia-EDR3 astrometry is given for most objects. The catalogue is available on NASA HEASARC and CDS and on its home page.

Announcing the release v2 of the MORX (Millions of Optical-Radio/X-ray Associations) catalogue which presents probable (40%-100% likelihood) radio/X-ray associations, including double radio lobes, to optical objects. All the large radio/X-ray surveys to June 2023 are included, being VLASS, LoTSS, RACS, FIRST, NVSS, and SUMSS radio surveys and Chandra, XMM-Newton, Swift, and ROSAT X-ray surveys. The totals are 3,115,575 optical objects of all classifications (or unclassified) so associated. The catalogue is available on multiple sites.

The number of satellites launched into low earth orbit has almost tripled (to over 4000) in the last three years due to the increasing commercialisation of space. Satellite constellations with a total of over 400,000 satellites are proposed to be launched in the near future. Many of these satellites are highly reflective, resulting in a high optical brightness that affects ground-based astronomical observations across the electromagnetic spectrum. Despite this, the potential effect of these satellites on Imaging Atmospheric Cherenkov Telescopes (IACTs) has so far been assumed to be negligible due to their nanosecond integration times. This has, however, never been verified. We aim to identify satellite trails in data taken by the High Energy Stereoscopic System (H.E.S.S.) IACT array in Namibia, using Night Sky Background (NSB) data from the CT5 camera installed in 2019. We determine which observation times and pointing directions are affected the most, and evaluate the impact on Hillas parameters used for classification and reconstruction of high-energy Extensive Air Shower events. Finally, we predict how future planned satellite launches will affect gamma-ray observations with IACTs.

The globular cluster Messier 80 was monitored by the Kepler space telescope for 80 days during the K2 mission. Continuous, high-precision photometry of such an old, compact cluster allows us to study its variable star population in unprecedented detail. We extract light curves for 27 variable stars using differential-image photometry. A search for new variables in the images led to the discovery of two new variable stars: an RR Lyrae and a variable red giant star, respectively. Analysis of the RR Lyrae population reveals multiple RRc stars with additional modes and/or peculiar modulation cycles. We newly classify star V28 as a spotted extreme horizontal branch variable. Despite their faintness, we clearly detect the three SX Phe stars but we did not find new pulsation modes beyond the known ones in them. Spectra taken with the VLT and Magellan Clay telescopes, as well as absolute color-magnitude diagrams of the cluster based on Gaia and Pan-STARRS observations confirm the classification of the peculiar modulated variables as bona-fide RRc stars. We propose that they highlight a subgroup of overtone stars that may have been overlooked before. We fit MESA isochrones to the CMDs to estimate the age and metallicity of the cluster. We confirm that M80 is old and metal-poor, but show that isochrone fitting to old populations comes with numerous uncertainties.

Solar flare is one of the most important solar activities which emit all electromagnetic waves in gigantic burst. The radio emission can be used to determine the physical properties of the solar flares. The e-CALLISTO worldwide network is designed to detect the radio emission of the solar flares and this study used the spectroscopic data from the e-CALLISTO system. Among the five types of solar radio bursts, this study was focused on type II radio bursts. The spectroscopic analysis estimated the shock speed of type II radio bursts using the uniform electron density model and the nonuniform electron density model of the sun. The shock speed is proportional to the electron density (Ne) and inversely proportional to the rate of change in electron density with altitude (dNe/dr). The determined shock speed at the altitude of one solar radius is 2131 km/s for uniform model and 766 km/s for non-uniform model. Although the uniform electron density model is widely used we attempted the non-uniform electron density since in the active region of the sun, the electron densities are non-uniform. The estimated shock speeds of uniform density model is relatively high so that it is reasonable to use non-uniform electron density model for shock speed estimation of type II radio bursts.

We present a comprehensive investigation of HI (super)clouds, molecular clouds (MCs), and star formation in the Carina spiral arm of the outer Galaxy. Utilizing HI4PI and CfA CO survey data, we identify HI clouds and MCs based on the ($l$, ${v_\mathrm{LSR}}$) locations of the Carina arm. We analyzed 26 HI clouds and 48 MCs. Most of the identified HI clouds are superclouds, with masses exceeding $10^6~{\mathrm{M_\odot}}$. We find that 15 of these superclouds have associated MC(s) with ${M_\mathrm{HI}} \gtrsim 10^6~{\mathrm{M_\odot}}$ and ${\Sigma_\mathrm{HI+H_2}} \gtrsim$ 50 ${\mathrm{M_\odot}} \rm pc^{-2}$. Our virial equilibrium analysis suggests that these CO-bright HI clouds are gravitationally bound or marginally bound. We report an anti-correlation between molecular mass fractions and Galactocentric distances, and a correlation with total gas surface densities. Nine CO-bright HI superclouds are associated with HII regions, indicating ongoing star formation. We confirm the regular spacing of HI superclouds along the spiral arm, which is likely due to some underlying physical process, such as gravitational instabilities. We observe a strong spatial correlation between HII regions and MCs, with some offsets between MCs and local HI column density peaks. Our study reveals that in the context of HI superclouds, the star formation rate surface density is independent of HI and total gas surface densities but positively correlates with molecular gas surface density. This finding is consistent with both extragalactic studies of the resolved Kennicutt-Schmidt relation and local giant molecular clouds study of Lada et al. (2013), emphasizing the crucial role of molecular gas in regulating star formation processes.

A CALLISTO system was set up at the Arthur C Clarke Institute and connected to the e-CALLISTO global network which observes the solar radio bursts in 24 hours. CALLISTO is the foremost observation facility to investigate celestial objects in radio region in Sri Lanka. The system consists of the CALLISTO spectrometer and controlling software,logarithmic periodic antenna and pre-amplifier. CALLISTO spectrometer is able to detect solar radio bursts in the frequency range of 45 MHz to 870 MHz with a channel resolution of 62.5 kHz.The log-periodic antenna was designed for 7 dBi gain and achieved the voltage standing wave ratio, less than 1.5 which is acquired by the overall impedance of the antenna, 49.3 ohms. The linear polarized antenna is pointing to zenith and the dipoles directed to north-south direction. The system detects solar radio emissions originated by solar flares and corona mass ejections. The radio bursts occurs as emission stripes in the radio spectra and classify from type I to V mainly on drift rate and band width. The system observed a type III solar radio burst on 5th July 2013 and a type II burst on 25th October 2013 which was originated by X1.7 solar flare. The type II bursts characterize with narrow bandwidth and drift slowly from higher to lower frequencies while the main features of type III bursts are high drift rate and broad bandwidth.

We use the Horizon Run 5 cosmological simulation to study the effect of galaxy intrinsic properties and the local environment on AGNs characterized by their threshold of the accretion rate. We select galaxies in the stellar mass range $10^{9.5} \le M^{}{*}/M^{}{\odot} \le 10^{10.5}$ in the snapshot at redshift $z$=0.625. Among various intrinsic properties, we find that the star formation rate of the host galaxy is most correlated to the AGN activity. To quantify the environment, we use background galaxy number density (large-scale environment) and distance and morphological type of the nearest neighbors (small-scale environment), and study their relative effects on the AGN properties. We find that, compared to the background density, the nearest neighbor environment is the dominant quantity determining the bolometric luminosity, star formation rate, and kinematic properties of AGNs and better dictates the gas mass of the host galaxy. We show that the cold gas content in the host galaxies is crucial in triggering AGN activity. However, when the nearest neighbor environment effects start to act at the neighbor distance of less than about half the virial radius of the neighbor, the neighbor environmental effects are the most dominant factor for quasar activity.

Gravitational waves from isolated sources have eluded detection so far. The upper limit of long-lasting continuous gravitational wave emission is now at the stage of probing physically-motivated models with the most optimistic being strongly constrained. One potential avenue to remedy this is to relax the assumption of the gravitational wave being quasi-infinite in duration, leading to the idea of transient continuous gravitational waves. In this paper, we outline how to get transient continuous waves from magnetars (or strongly-magnetised neutron stars) that exhibit glitches and/or antiglitches. We put forward a toy model whereby at a glitch or antiglitch, mass is ejected from the magnetar but becomes trapped on its outward journey through the magnetosphere. Depending on the specific values of the height of the trapped ejecta and the magnetic inclination angle, we are able to reproduce both glitches and antiglitches from simple angular momentum arguments. The trapped ejecta sets the magnetar into precession causing gravitational waves to be emitted at once and twice the magnetar's spin frequency, for a duration equal to however long the ejecta is trapped for. We find that the gravitational waves are more likely to be detectable when the magnetar is: closer, rotating faster, or has larger glitches/antiglitches. Specific to this model, we find that the detectability improves when the ejecta height and magnetic inclination angle have values near the boundary in the parameter space that separates glitches and antiglitches, though this requires more mass to be ejected to remain consistent with the observed glitch/antiglitch.

Accurate polarimetric calibration of the radio pulse profiles from pulsars is crucial for studying their radiation properties at these wavelengths. Inaccurate calibration can also distort recorded pulse profiles, introducing noise in time of arrival (TOA) data and thus degrading pulsar timing analyses. One method for determining the full polarimetric response of a given telescope is to conduct observations of bright polarized pulsars over wide ranges of parallactic angles, to sample different orientations of their polarization angle and determine the cross-couplings between polarization feeds. The Nan\c{c}ay decimetric Radio Telescope (NRT) is a 94m equivalent meridian telescope, capable of tracking a given pulsar for approximately one hour around transit. In November 2019, we began conducting regular observations of the bright and highly linearly polarized pulsar PSR~J0742$-$2822, in a special mode where the feed horn rotates by $\sim 180^\circ$ over the course of the one hour observation, mimicking wide parallactic angle variations and enabling us to determine the polarimetric response of the NRT at 1.4~GHz. The improved polarimetric response of the NRT as determined from these observations was applied to observations of a selection of MSPs with published polarimetric properties. We find that the new polarimetric profiles and polarization position angles are consistent with previous findings, unlike NRT polarimetric results obtained with the previously used method of calibration. The analysis of timing data on J1730$-$2304, J1744$-$1134, and J1857+0953 shows that the new calibration method improves the quality of the timing, and the Matrix Template Matching (MTM) method proves very effective at reducing noise from imperfect calibration. For pulsars with sufficient degrees of polarization, the MTM method appears to be the preferred method for extracting TOAs from NRT observations.

We present a method to measure the the oblateness parameter q of the dark matter halos of gas rich galaxies that have extended HI disks. We have applied our model to a sample of 20 nearby galaxies that are gas rich and close to face-on, of which 6 are large disk galaxies, 8 have moderate stellar masses and 6 are low surface brightness (LSB) dwarf galaxies. We have used the stacked HI velocity dispersion and HI surface densities to derive q in the outer disk regions. Our most important result is that gas dominated galaxies (such as LSB dwarfs) that have M(gas)/M(baryons)>0.5 have oblate halos (q<0.55), whereas stellar dominated galaxies have a range of q values from 0.2 to 1.3. We also find a significant positive correlation between q and stellar mass, which indicates that galaxies with massive stellar disks have a higher probability of having halos that are spherical or slightly prolate, whereas low mass galaxies preferably have oblate halos. We briefly also discuss how the halo shape affects the disks of galaxies, especially the oblate halos.

Within the quiet Sun corona imaged at 1 MK, much of the field of view consists of diffuse emission that appears to lack the spatial structuring that is so evident in coronal loops or bright points. We seek to determine if these diffuse regions are categorically different in terms of their intensity fluctuations and spatial configuration from the more well-studied dynamic coronal features. We analyze a time series of observations from Solar Orbiter's High Resolution Imager in the Extreme Ultraviolet to quantify the characterization of the diffuse corona at high spatial and temporal resolutions. We then compare this to the dynamic features within the field of view, mainly a coronal bright point. We find that the diffuse corona lacks visible structuring, such as small embedded loops, and that this is persistent over the 25 min duration of the observation. The intensity fluctuations of the diffuse corona, which are within +/-5%, are significantly smaller in comparison to the coronal bright point. Yet, the total intensity observed in the diffuse corona is of the same order as the bright point. It seems inconsistent with our data that the diffuse corona is a composition of small loops or jets or that it is driven by discrete small heating events that follow a power-law-like distribution. We speculate that small-scale processes like MHD turbulence might be energizing the diffuse regions, but at this point we cannot offer a conclusive explanation for the nature of this feature.

We establish a sample of 370 Mira variables that are likely near the Galactic center (GC). The sources have been selected from the OGLE and BAaDE surveys based on their sky coordinates, OGLE classifications, and BAaDE maser-derived line-of-sight velocities. As the distance to the GC is known to a high accuracy, this sample is a test bed for reddening and extinction studies toward the GC and in Mira envelopes. We calculated separate interstellar- and circumstellar-extinction values for individual sources, showing that there is a wide range of circumstellar extinction values (up to four magnitudes in the K$_s$ band) in the sample, and that circumstellar reddening is statistically different from interstellar reddening laws. Further, the reddening laws in the circumstellar environments of our sample and the circumstellar environments of Large Magellanic Cloud (LMC) Miras are strikingly similar despite the different metallicities of the samples. Period-magnitude relations for the mid-infrared (MIR) WISE and MSX bands are also explored, and in the WISE bands we compare these to period-magnitude relationships derived from Miras in the LMC as it is important to compare these LMC relations to those in a higher metallicity environment. Emission from the envelope itself may contaminate MIR magnitudes altering the relations, especially for sources with thick envelopes.

We develop a new theoretical framework for studying the pairwise and cross-pairwise polarised kinetic Sunyaev Zeldovich (pkSZ) effect arising from the transverse peculiar velocity of galaxy clusters. The pkSZ effect is second order in peculiar velocities and has a spectrum that can be decomposed into y-type and blackbody components, whereas the unpolarised linear kSZ effect has only the blackbody component. Thus, the detectability of the pkSZ effect depends only on the sensitivity and the number of frequency channels of the survey and not on the other primary and secondary CMB anisotropies. We consider pairing of clusters with other clusters as well as cross-pairing of clusters with galaxies from spectroscopic galaxy surveys. The pairwise pkSZ signal is a function of intra-pair spatial separation. We develop and compare estimators of the pairwise pkSZ effect and study the detectability of the pairwise signal with cluster catalogs consisting of a few hundred thousand clusters expected from surveys such as eROSITA and CMB-S4. We find that cross-pairing clusters with galaxies from a large overlapping spectroscopic survey having a few $\mathrm{billion}$ galaxies will enable us to detect the pairwise pkSZ effect with CMB-S4. The pairwise pkSZ effect will thus open up a new window into the large-scale structure of the Universe in the coming decades.

The cores of massive neutron stars provide a unique environment for the dense nuclear matter in the universe. The global properties of a neutron star and gravitational waves emitted from the binary neutron star merger carry information about dense nuclear matter. We study in this paper the effect of the possible hadron-quark transition on the properties of the neutron star and the gravitational waves emitted from the binary neutron star merger by using the equations of state constructed from the Maxwell ansatz, Gibbs ansatz and, the crossover scenario. Our results show that the short period of the inspiral phase and the earlier collapse to a black hole indicate a soft equation of state. In combination with the future detection of the $10$kHz gravitational waves emitted from the binary neutron star merger and the signals from the electromagnetic counterparts, we expect the present study could reveal some characters of the dense nuclear matter.

Primordial black holes, which could have formed during the early Universe through overdensities in primordial density fluctuations during inflation, are potential candidates for dark matter. We explore the use of lensing parallax of Gamma ray bursts, which results in different fluxes being observed from two different vantage points, in order to probe the abundance of primordial black holes in the unexplored window within the mass range $[10^{-15}-10^{-11}]M_\odot$. We derive the optical depth for the detectability of lensing of GRBs with a distribution of source properties and realistic detector sensitivities. We comment on the ability of the proposed Indian twin satellite mission Daksha in its low earth orbit to conduct this experiment. If the two Daksha satellites observe 10000 GRBs simultaneously and the entirety of dark matter is made up of $[10^{-15}-10^{-12}]M_\odot$ black holes, Daksha will detect non-zero lensing events with a probability ranging from [80, 50] per cent. Non-detections will not conclusively rule out primordial black holes as dark matter in this mass range. However, we show that meaningful constraints can be obtained in such a case if the two satellites are separated by at least the Earth-Moon distance.

We present two independent measurements of stellar velocity dispersions ( $\sigma_\rm{\star}$ ) from the Ca\,H+K \& Mg\,\textsc{i} region (3880--5550~\AA) and the Calcium Triplet region (CaT, 8350--8750~\AA) for 173 hard X-ray-selected Type 1 AGNs ($z \leq$ 0.08) from the 105-month Swift-BAT catalog. We construct one of the largest samples of local Type 1 AGNs that have both single-epoch (SE) 'virial' black hole mass ($M_\rm{BH}$) estimates and $\sigma_\rm{\star}$ measurements obtained from high spectral resolution data, allowing us to test the usage of such methods for SMBH studies. We find that the two independent $\sigma_\rm{\star}$ measurements are highly consistent with each other, with an average offset of only $0.002\pm0.001$ dex. Comparing $M_\rm{BH}$ estimates based on broad emission lines and stellar velocity dispersion measurements, we find that the former is systematically lower by $\approx$0.12 dex. Consequently, Eddington ratios estimated through broad-line $M_\rm{BH}$ determinations are similarly biased (but in the opposite way). We argue that the discrepancy is driven by extinction in the broad-line region (BLR). We also find an anti-correlation between the offset from the $M_\rm{BH}$ - $\sigma_\rm{\star}$ relation and the Eddington ratio. Our sample of Type 1 AGNs shows a shallower $M_\rm{BH}$ - $\sigma_\rm{\star}$ relation (with a power law exponent of $\approx$3.5) compared with that of inactive galaxies (with a power-law exponent of $\approx$4.5), confirming earlier results obtained from smaller samples.

We study the correlation between the non-thermal velocity dispersion ($\sigma_{\rm nth}$) and the length-scale (L) in the neutral interstellar medium (ISM) using a large number of Hi gas components taken from various published Hi surveys and previous Hi studies. We notice that above the length-scale ($L$) of 0.40 pc, there is a power-law relationship between $\sigma_{\rm nth}$ and $L$. However, below 0.40 pc, there is a break in the power-law, where $\sigma_{\rm nth}$ is not significantly correlated with $L$. It has been observed from the Markov chain Monte Carlo (MCMC) method that for the dataset of $L > 0.40$ pc, the most probable values of intensity ($A$) and power-law index ($p$) are 1.14 and 0.55 respectively. Result of $p$ suggests that the power-law is steeper than the standard Kolmogorov law of turbulence. This is due to the dominance of clouds in the cold neutral medium. This is even more clear when we separate the clouds into two categories: one for $L$ is > 0.40 pc and the kinetic temperature ($T_k$ ) is < 250 K, which are in the cold neutral medium (CNM) and for other one where L is > 0.40 pc and T k is between 250 K and 5000 K, which are in the thermally unstable phase (UNM). Most probable values of $A$ and $p$ are 1.14 and 0.67 respectively in the CNM phase and 1.01 and 0.52 respectively in the UNM phase. A greater number of data points is effective for the UNM phase in constructing a more accurate estimate of $A$ and $p$, since most of the clouds in the UNM phase lie below 500 K. However, from the value of $p$ in the CNM phase, it appears that there is a significant difference from the Kolmogorov scaling, which can be attributed to a shock-dominated medium.

Many supernovae (SNe) imply an interaction of the SN ejecta with matter (CSM) surrounding the progenitor star. This suggests that many massive stars may undergo various degrees of envelope stripping shortly before exploding, and produce a considerable diversity in their pre-explosion CSM properties. We explore a generic set of ~100 detailed massive binary evolution models to characterize the amount of envelope stripping and the expected CSM configurations. Our binary models were computed with the MESA stellar evolution code, considering an initial primary star mass of 12.6 Msun, and focus on initial orbital periods above 500 d. We compute these models up to the time of the primary's iron core collapse. We find that Roche lobe overflow often leads to incomplete stripping of the mass donor, resulting in a large variety of pre-SN envelope masses. Many of our models' red supergiant (RSG) donors undergo core collapse during Roche lobe overflow, with mass transfer and thus system mass loss rates of up to 0.01 Msun/yr at that time. The corresponding CSM densities are similar to those inferred for Type IIn SNe like 1998S. In other cases, the mass transfer turns unstable, leading to a common envelope phase at such late time that the mass donor explodes before the common envelope is fully ejected or the system has merged. We argue that this may cause significant pre-SN variability, as for example in SN 2020tlf. Other models suggest a common envelope ejection just centuries before core collapse, which may lead to the strongest interactions, as in superluminous Type IIn SNe like 1994W, or 2006gy. Wide massive binaries offer a natural framework to understand a broad range of hydrogen-rich interacting SNe. On the other hand, the flash features observed in many Type IIP SNe, like in SN 2013fs, may indicate that RSGs are more extended than currently assumed.

Through their radio loudness, lack of thermal UV emission from the accretion disk and power-law dominated spectra, Low Luminosity AGN (LLAGN) display similarity with the hard state of stellar-mass black hole X-Ray Binaries (BHBs). In this work we perform a systematic hard X-ray spectral study of a carefully selected sample of unobscured LLAGN using archival $NuSTAR$ data, to understand the central engine properties in the lower accretion regime. We analyze the $NuSTAR$ spectra of a sample of 16 LLAGN. We model the continuum emission with detailed Comptonization models. We find a strong anti-correlation between the optical depth and the electron temperature of the corona, previously also observed in the brighter Seyferts. This anti-correlation is present irrespective of the shape of the corona, and the slope of this anti-correlation in the log space for LLAGN (0.68-1.06) closely matches that of the higher accretion rate Seyferts (0.55-1.11) and hard state of BHBs ($\sim$0.87). This anti-correlation may indicate a departure from a fixed disk-corona configuration in radiative balance. Our result, therefore, demonstrates a possible universality in Comptonization processes of black hole X-ray sources across multiple orders of magnitude in mass and accretion rate.

The next Galactic core-collapse supernova (CCSN) presents a once-in-a-lifetime opportunity to make astrophysical measurements using neutrinos, gravitational waves, and electromagnetic radiation. CCSNe local to the Milky Way are extremely rare, so it is paramount that detectors are prepared to observe the signal when it arrives. The IceCube Neutrino Observatory, a gigaton water Cherenkov detector below the South Pole, is sensitive to the burst of neutrinos released by a Galactic CCSN at a level $>$10$\sigma$. This burst of neutrinos precedes optical emission by hours to days, enabling neutrinos to serve as an early warning for follow-up observation. IceCube's detection capabilities make it a cornerstone of the global network of neutrino detectors monitoring for Galactic CCSNe, the SuperNova Early Warning System (SNEWS 2.0). In this contribution, we describe IceCube's sensitivity to Galactic CCSNe and strategies for operational readiness, including "fire drill" data challenges. We also discuss coordination with SNEWS 2.0.

The overlap of galaxy surveys and CMB experiments presents an ideal opportunity for joint cosmological dataset analyses. In this paper we develop a halo-model-based method for the first joint analysis combining these two experiments using 10 correlated two-point functions (10x2pt) derived from galaxy position, galaxy shear, CMB lensing convergence, and Compton-y fields. We explore this method using the Vera Rubin Observatory Legacy Survey of Space and Time (LSST) and the Simons Observatory (SO) as examples. We find such LSSxCMB joint analyses lead to significant improvement in Figure-of-Merit of $\Omega_m$ and $S_8$ over the constraints from using LSS-only probes within $\Lambda$CDM. We identify that the shear-$y$ and $y$-$y$ correlations are the most valuable additions when tSZ is included. We further identify the dominant sources of halo model uncertainties in the small-scale modelling, and investigate the impact of halo self-calibration due to the inclusion of small-scale tSZ information.

Apsidal motion is the precession of the line of apsides in the orbit of a binary star due to perturbations from General Relativity (GR), tides, or third body interactions. The rate of precession due to tidal effects depends on the interior structures of the stars and, as a result, binaries where such precession occurs are of great interest. Apsidal motion is observed through the analysis of eclipse times, revealing small changes in the average interval between successive primary and secondary eclipses, taking all available observed times of eclipse and yielding an estimate of the apsidal rate. Given that this is a single observed quantity, various degeneracies are unavoidably present. Ideally, one would have a model that predicts eclipse times given the orbital and stellar parameters. These parameters for a given binary could then be computed using least squares, provided a suitably large number of eclipse times. Here we use the Eclipsing Light Curve (ELC) program as such a model. The Newtonian equations of motion with extra force terms accounting for GR contributions and tidal distortions are integrated, yielding precise sky positions as a function of time. Times of mid-eclipse and instantaneous orbital elements are computed as a function of time. In this paper we outline the method, and compare numerically computed apsidal rates with standard formulae using a set of 15 binaries based on real systems. For our simulated systems, the derived apsidal rates agree with the standard formula.

We present the construction of an image similarity retrieval engine for the morphological classification of galaxies using the Convolutional AutoEncoder (CAE). The CAE is trained on 90,370 preprocessed Sloan Digital Sky Survey galaxy images listed in the Galaxy Zoo 2 (GZ2) catalog. The visually similar output images returned by the trained CAE suggest that the encoder efficiently compresses input images into latent features, which are then used to calculate similarity parameters. Our Tool for Searching a similar Galaxy Image based on a Convolutional Autoencoder using Similarity (TSGICAS) leverages this similarity parameter to classify galaxies' morphological types, enabling the identification of a wider range of classes with high accuracy compared to traditional supervised ML techniques. This approach streamlines the researcher's work by allowing quick prioritization of the most relevant images from the latent feature database. We investigate the accuracy of our automatic morphological classifications using three galaxy catalogs: GZ2, Extraction de Formes Id\'ealis\'ees de Galaxies en Imagerie (EFIGI), and Nair $\&$ Abraham (NA10). The correlation coefficients between the morphological types of input and retrieved galaxy images were found to be 0.735, 0.811, and 0.815 for GZ2, EFIGI, and NA10 catalogs, respectively. Despite differences in morphology tags between input and retrieved galaxy images, visual inspection showed that the two galaxies were very similar, highlighting TSGICAS's superior performance in image similarity search. We propose that morphological classifications of galaxies using TSGICAS are fast and efficient, making it a valuable tool for detailed galaxy morphological classifications in other imaging surveys.

We present an analysis of Type Ia Supernovae (SNe~Ia) from both the Carnegie Supernova Project~I (CSP-I) and II (CSP-II), and extend the Hubble diagram from the optical to the near-infrared wavelengths ($uBgVriYJH$). We calculate the Hubble constant, $H_0$ using various distance calibrators: Cepheids, Tip of the Red Giant Branch (TRGB), and Surface Brightness Fluctuations (SBF). Combining all methods of calibrations, we derive $\rm H_0=71.43 \pm 0.62 \ (stat) \pm 2.43 \ (sys) \ km \ s^{-1} \ Mpc^{-1}$ from $B$-band, and $\rm H_0=72.65 \pm 0.63 \ (stat) \pm 2.88 \ (sys) \ km \ s^{-1} \ Mpc^{-1}$ from $H$-band. By assigning equal weight to the Cepheid, TRGB, and SBF calibrators, we derive the systematic errors required for consistency in the first rung of the distance ladder, resulting in an increased systematic error in $H_0$. As a result, the tension between the late-time $H_0$ we derive by combining the various distance calibrators and the early-time $H_0$ from the Cosmic Microwave Background is reduced. The highest precision in SN~Ia luminosity is found in the $Y$ band ($0.12\pm0.01$ mag), as defined by the intrinsic scatter ($\sigma_{int}$). We revisit SN~Ia Hubble residual-host mass correlations and recover previous results that these correlations do not change significantly between the optical and the near-infrared wavelengths. Finally, SNe~Ia that explode beyond 10 kpc from their host centers exhibit smaller dispersion in their luminosity, confirming our earlier findings. Reduced effect of dust in the outskirt of hosts may be responsible for this effect.

We carry out experiments of 104 m/s velocity oblique impacts into a granular medium (sand). Impact craters have nearly round rims even at a grazing angle of about $10^\circ$, however, the strength of seismic pulses excited by the impact is dependent upon impact angle, and the ratio between uprange and downrange velocity peaks can be as large as 5, particularly at shallow depths. Crater slope, an offset between crater center and impact site, crater volume, azimuthal variation in ejection angle, seismic pulse shapes and subsurface flow direction are also sensitive to impact angle, but to a much lower degree than subsurface pulse strength. Uprange and downrange pulse peak amplitudes can be estimated from the horizontal and vertical components of the momentum imparted to the medium from the projectile

We present a summary of the results of the OJ 287 observational campaign, which was carried out during the 2021/2022 observational season. This season is special in the binary model because the major axis of the precessing binary happens to lie almost exactly in the plane of the accretion disc of the primary. This leads to pairs of almost identical impacts between the secondary black hole and the accretion disk in 2005 and 2022. In 2005, a special flare called "blue flash" was observed 35 days after the disk impact, which should have also been verifiable in 2022. We did observe a similar flash and were able to obtain more details of its properties. We describe this in the framework of expanding cloud models. In addition, we were able to identify the flare arising exactly at the time of the disc crossing from its photo-polarimetric and gamma-ray properties. This is an important identification, as it directly confirms the orbit model. Moreover, we saw a huge flare that lasted only one day. We may understand this as the lighting up of the jet of the secondary black hole when its Roche lobe is suddenly flooded by the gas from the primary disk. Therefore, this may be the first time we directly observed the secondary black hole in the OJ 287 binary system.

Many past studies have predicted the steady-state production and maintenance of abiotic O$_2$ and O$_3$ in the atmospheres of CO$_2$-rich terrestrial planets orbiting M dwarf stars. However, the time-dependent responses of these planetary atmospheres to flare events - and the possible temporary production or enhancement of false positive biosignatures therein - has been comparatively less well studied. Most past works that have modeled the photochemical response to flares have assumed abundant free oxygen like that of the modern or Proterozoic Earth. Here we examine in detail the photochemical impact of the UV emitted by a single flare on abiotic O$_2$/O$_3$ production in prebiotic, CO$_2$-dominated atmospheres of M dwarf planets with CO$_2$ levels ranging from 10% to 90% of 1 bar. We find that a single flare generally destroys O$_2$ while modestly enhancing O$_3$ column densities. We simulate the spectral observables of both the steady-state atmosphere and time-dependent spectral response over the flare window for both emitted and transmitted light spectra. Over the course of the flare, the O$_3$ UV Hartley band is modestly enhanced by a maximum of 6 ppm while the CO$_2$ molecular transit depths modestly decline by 7 ppm. In both emitted and transmitted light spectra, the 9.65 $\mu$m O$_3$ band is hidden by the overlapping 9.4 $\mu$m CO$_2$ band for all scenarios considered. Overall, we find that the possible enhancements of abiotic O$_3$ due to a single flare are small compared to O$_3$'s sensitivity to other parameters such as CO$_2$ and H$_2$O abundances or the availability of reducing gases such as H$_2$.

The relation between giant radio halos and mini-halos in galaxy clusters is not understood. The former are usually associated with merging clusters, the latter are found in relaxed systems. In the last years, the advent of low-frequency radio observations has challenged this dichotomy, finding intermediate objects with a hybrid radio morphology. We aim to investigate the presence of diffuse radio emission in the cluster Abell 1413 and determine its dynamical status. We used LOFAR HBA observations centred at 144 MHz to study the diffuse emission hosted by this cluster.To investigate the dynamical state of the system, we complete our study with newly analysed XMM-Newton archival data. A1413 shows features that are typically present in both relaxed (e.g., peaked x-ray surface brightness distribution and little large-scale inhomogeneities) and disturbed (e.g., flatter temperature and metallicity profiles) clusters.This evidence supports the scenario that A1413 is neither a disturbed nor fully relaxed object. We argue that it is an intermediate-phase cluster.Using radio observations at 144 MHz, we discover the presence of a wider diffuse component surrounding the previously reported mini-halo at the cluster centre. By fitting the radio surface brightness profile with a double-exponential model, we can disentangle the two components. We find an inner mini-halo with an e-folding radius r_e1=28 kpc and the extended component with r_e2 = 290 kpc. We also performed point-to-point correlations between radio and X-ray surface brightness, finding a sub-linear relation for the outer emission and a super-linear relation for the mini-halo.The mini-halo and the diffuse emission extend over different scales and show different features, confirming the double nature of the radio emission and suggesting that the mechanisms responsible for the re-acceleration of the radio-emitting particle might be different.

We performed a search for faint low surface brightness dwarf galaxies around the major spiral galaxy M\,101 and in the large rectangular area within SGL = [30 -- 80]$^{\circ}$, SGB =[10 -- 37]$^{\circ}$ spanning a chain of galaxies: M\,63, M\,51, M\,101, and NGC\,6503, based on the data from DESI Legacy Imaging Surveys. Six new supposed dwarf members of the complex were discovered. We present a list of 25 prospective members of the M\,101 group and estimate the total mass and the total-mass-to-$K$-band luminosity ratio of the group as $(1.02\pm0.42)\times10^{12}~M_{\odot}$ and $(16.0\pm6.5)~M_{\odot}/L_{\odot}$, respectively. We notice that the average dark mass-to-luminosity ratio in the groups around M\,63, M\,51, and M\,101 is $(12\pm4)M_{\odot}/L_{\odot}$ that almost an order of magnitude lower than the global cosmic ratio, $(102\pm5)M_{\odot}/L_{\odot}$.

This work aims at investigating the optical transmission system needed for such lightweight sail, taking into account the physical constraints of such unprecedented link and focusing on the optimal scheme for the optical signal emission. In particular, the optical signal is distributed to several emitters on the sail. The light diffraction resulting from the pattern of the emitters acting coherently determines the characteristics of the whole beam transmitted by the sail and of the received signal on the Earth. The performance of the digital communication system using pulse position modulation (PPM) can be assessed and channel coding schemes are proposed. We are using the paradigm for which the entire sail communication system is described as a Tree-of-light: the detectors, CPU, memory and laser transmitter are the central unit, representing the trunk of the tree. The branches of the tree are the waveguides, directed to the sail surface. By means of multimode splitters, the signal is further distributed via the petioles to the emitters, the leaves, realized by grating couplers (GCs), on which this work is more focused.

Despite their relatively high intrinsic brightness and the fact that they are more numerous than younger OB stars and kinematically colder than older red giants, A-type stars have rarely been used as Galactic tracers. They may, in fact, be used to fill the age gap between these two tracers, thereby allowing us to study the transition between them. We analyse Galactic disc structure and kinematic perturbations up to 6 kpc from the Sun based on observations of A-type stars. This work presents a catalogue of A-type stars selected using the IGAPS photometric survey. It covers the Galactic disc within $30^{o}\leq l\leq215^{o}$ and $|b|\leq5^{o}$ up to a magnitude of $r\leq19$ mag with about 3.5 million sources. We used Gaia Data Release 3 parallaxes and proper motions, as well as the line-of-sight velocities, to analyse the large-scale features of the Galactic disc. We carried out a study of the completeness of the detected density distributions, along with a comparison between the $b<0^{o}$ and $b>0^{o}$ regions. Possible biases caused by interstellar extinction or by the usage of some kinematic approximations were examined as well. We find stellar overdensities associated with the Local and the Perseus spiral arms, as well as with the Cygnus region. A-type stars also provide kinematic indications of the Galactic warp towards the anticentre, which displays a median vertical motion of ~6-7 km/s at a Galactocentric radius of R=14 kpc. It starts at R=12 kpc, which supports the scenario where the warp begins at larger radii for younger tracers when compared with other samples in the literature. We also detect a region with downward mean motion extending beyond 2 kpc from the Sun towards $60^{o}<l<75^{o}$ that may be associated with a compression breathing mode. Furthermore, A-type stars reveal very clumpy inhomogeneities and asymmetries in the $V_Z$-$V_{\phi}$ velocity space plane.

A new interpretation of the optical knot in the jet of M87, HST-1, is presented. High sensitivity 22 GHz Very Large Array images locate HST-1 to within 6 mas of the jet axis immediately upstream. 1.7 GHz Very Long Baseline Array images of a bright flare in 2005 indicates that the preponderance of emission in the early stages originates in an elongated region that is tilted $12.5^{\circ}$ from the jet axis. The superluminal motion, shape, location and the large jet-aligned optical/UV polarization suggest an identification with the putative relativistic spine of the jet. As such, energy flux estimates for HST-1, $\sim 870$ mas from the nucleus, published in 2006 indicate that the central engine injected $Q_{\rm{spine}}\approx 2.5 \times 10^{41}\rm{ergs/s}$ into the base of the spine $\sim 200$ years earler. Furthermore, previous studies reveal a tubular protonic jet on sub-mas scales that envelopes a low luminosity core, presumably the faint spine base. It was estimated that the central engine injected $Q_{\rm{tubular\,jet}}\approx 6.1\times 10^{41}\rm{ergs/s}$ $\sim 1.5$ years earlier. If one component of the jet is inherently more powerful, a firm constraint on total jet power in the recent past exists. If the emitted jet is inherently dominated by the spine (tubular jet) then the total bilaterally symmetric jet power emitted from the central engine was $<4Q_{\rm{spine}}\approx 1.0 \times 10^{42}\rm{ergs/s}$ ($< 4Q_{\rm{tubular\,jet}}\approx 2.4\times 10^{42}\rm{ergs/s}$) $\sim 200$ ($\sim 1.5$) years earlier. Assuming a nearly constant central engine injected jet power for $\sim 200$ years indicates a total jet power of $\lesssim 2\times 10^{42}$ ergs/s in epochs of modern observation or $\lesssim 3.5\%$ jet production efficiency for an accretion rate of 0.001$M_{\odot}$/yr.

In this paper, we construct a bounce inflation cosmological scenario in the framework of the modified symmetric teleparallel gravity, namely f(Q) theory, and investigate the tensor perturbations therein. As is well-known, the tensor perturbations generated in the very early Universe (inflation and pre-inflation regions) can account for the primordial gravitational waves (PGWs) that are to be detected by the next generation of GW experiments. We discuss the stability condition of the tensor perturbations in the bounce inflation process and investigate in detail the evolution of the perturbation variable. The general form of the tensor power spectrum is obtained both for large as well as small scale modes. As a result, we show for both kinds of modes (short or long wavelength modes), the tensor spectrum may get a positive tilt in the parametric range where the tensor perturbation proves to be stable -- this interestingly hints an enhancement of gravitational waves' amplitude in the background of the f(Q) bounce-inflation scenario. Moreover, we study the LQC-like scenario as a specific case of our model, in which, the primordial tensor power spectrum turns out to be nearly scale-invariant on both small and large scales.

Various theoretical models predict the existence of exotic compact objects that can mimic the properties of black holes (BHs). Gravitational waves (GWs) from the mergers of compact objects have the potential to distinguish between exotic compact objects and BHs. The measurement of spin-induced multipole moments of compact objects in binaries provides a unique way to test the nature of compact objects. The observations of GWs by LIGO and Virgo have already put constraints on the spin-induced quadrupole moment, the leading order spin-induced moment. In this work, we develop a Bayesian framework to measure the spin-induced octupole moment, the next-to-leading order spin-induced moment. The precise measurement of the spin-induced octupole moment will allow us to test its consistency with that of Kerr BHs in general relativity and constrain the allowed parameter space for non-BH compact objects. For various simulated compact object binaries, we explore the ability of the LIGO and Virgo detector network to constrain spin-induced octupole moment of compact objects. We find that LIGO and Virgo at design sensitivity can constrain the symmetric combination of component spin-induced octupole moments of binary for dimensionless spin magnitudes $\sim 0.8$. Further, we study the possibility of simultaneously measuring the spin-induced quadrupole and octupole moments. Finally, we perform this test on selected GW events reported in the third GW catalog. These are the first constraints on spin-induced octupole moment using full Bayesian analysis.

There are observations indicating a possible anomalous transparency of intergalactic space (filled with infrared background light) for extragalactic gamma-rays of very high energy. The anomaly is usually associated with effects of some new physics. However, another explanation is possible -- as a manifestation relating to a cut-off of the zero-point vibration spectrum.

It is known that theories of phantom dark energy, considered as quantum fields, predict a continuous production of positive- plus negative-energy particles, from spontaneous decay of the vacuum. We show that this can be a new source of boosted dark matter or radiation, with consequences for direct detection. We set constraints on such models using data from the XENONnT experiment, and we show that recent excess events reported by the DAMIC experiment can be consistently described as coming from dark radiation, produced by vacuum decay, interacting with electrons.

Axion strings are horizon-size topological defects that may be produced in the early Universe. Ultra-light axion-like particles may form strings that persist to temperatures below that of big bang nucleosynthesis. Such strings have been considered previously as sources of gravitational waves and cosmic microwave background (CMB) polarization rotation. In this work we show, through analytic arguments and dedicated adaptive mesh refinement cosmological simulations, that axion strings deposit a sub-dominant fraction of their energy into high-energy Standard Model (SM) final states, for example, by the direct production of heavy radial modes that subsequently decay to SM particles. This high-energy SM radiation is absorbed by the primordial plasma, leading to novel signatures in precision big bang nucleosynthesis, the CMB power spectrum, and gamma-ray surveys. In particular, we show that CMB power spectrum data constrains axion strings with decay constants $f_a \lesssim 10^{12}$ GeV, up to model dependence on the ultraviolet completion, for axion masses $m_a \lesssim 10^{-29}$ eV; future CMB surveys could find striking evidence of axion strings with lower decay constants.

We explore the impact of highly excited bound states on the evolution of number densities of new physics particles, specifically dark matter, in the early Universe. Focusing on dipole transitions within perturbative, unbroken gauge theories, we develop an efficient method for including around a million bound state formation and bound-to-bound transition processes. This enables us to examine partial-wave unitarity and accurately describe the freeze-out dynamics down to very low temperatures. In the non-Abelian case, we find that highly excited states can prevent the particles from freezing out, supporting a continuous depletion in the regime consistent with perturbativity and unitarity. We apply our formalism to a simplified dark matter model featuring a colored and electrically charged $t$-channel mediator. Our focus is on the regime of superWIMP production which is commonly characterized by a mediator freeze-out followed by its late decay into dark matter. In contrast, we find that excited states render mediator depletion efficient all the way until its decay, introducing a dependence of the dark matter density on the mediator lifetime as a novel feature. The impact on the viable dark matter mass can amount to an order of magnitude, relaxing constraints from Lyman-$\alpha$ observations.

We revisit the formation of a thermal population of hadronic axions in nonstandard cosmologies, in light of the recent developments in obtaining continuous and smooth interaction rates for both the gluon and photon couplings. For certain cosmological histories, such as low-temperature reheating (LTR) and kination-like scenarios, the thermalization of the axion can be severely delayed to higher masses. In the case that thermal equilibrium is achieved, we improve the constraints on LTR for axion masses around the eV scale with respect to previous works and we constrain for the first time early matter-dominated (EMD) cosmologies. We also point out the possibility of having the co-existence of cold and warm dark matter populations of axions in kination-like scenarios in the eV mass range.

Dynamical cancellation frameworks present a potential means of mitigating the effect of a large vacuum energy, that would otherwise ruin the late-time, low energy dynamics of the Universe. Certain models in the literature, such as the Fab Four and Well Tempering, realize this idea by introducing some degeneracy in the dynamical equations. In this paper, we introduce a third potential route to self-tuning, and infer the existence of a new, exact Milne solution in the simplest tadpole plus cubic-Galileon scalar-tensor theory. We study the dynamics of the scalar field and metric in the vicinity of the Milne coordinate singularity, and find that the vacuum solution belongs to a more general family of Milne-like metrics. By numerically evolving the field equations for a range of initial conditions, we show that the Milne solution is not an attractor, and varying the initial scalar field data can lead to completely different asymptotic states; exponential growth of the scale factor, a static non-spatially flat metric or a severe finite-time instability in the scalar field and metric. We generalise the Milne solution to a class of FLRW spacetimes, finding that the tadpole-cubic Galileon model admits perfect-fluid-like solutions in the presence of matter. Finally, we present a second Horndeski model which also admits an exact Milne solution, hinting at the existence of a larger undiscovered model space containing vacuum-energy-screened solutions.

We propose a novel method for testing gravity models using seismic data from Earth. By imposing observational constraints on Earth's moment of inertia and mass, we rigorously limit the gravitational models' parameters within a $2\sigma$ accuracy. Our method constrains the parameters governing additional terms to the General Relativity Lagrangian to the following ranges: $-2\times10^9\lesssim\beta\lesssim 10^9 \text{m}^2$ for Palatini $f(R)$ gravity, $-8\times10^9\lesssim\epsilon\lesssim 4\times 10^9 \text{m}^2$ for Eddington-inspired Born-Infeld gravity, and $-10^{-3}\lesssim\Upsilon\lesssim10^{-3}$ for Degenerate Higher-Order Scalar-Tensor theories. We also discuss potential avenues to enhance the proposed method, aiming to impose even tighter constraints on gravity models.

We present the conservative effective two-body Hamiltonian at the third order in the post-Newtonian expansion with gravitoelectric quadrupolar dynamical tidal-interactions. Our derivation of the effective two-body Lagrangian is based on the diagrammatic effective field theory approach and it involves Feynman integrals up to three loops, which are evaluated within the dimensional regularization scheme. The elimination of the divergent terms occurring in the effective Lagrangian requires the addition of counterterms to ensure finite observables, thereby introducing a renormalization group flow to the post-adiabatic Love number. As a limiting case of the renormalized dynamical effective Hamiltonian, we also derive the effective Hamiltonian for adiabatic tides, and, in this regime, calculate the binding energy for a circular orbit, and the scattering angle in a hyperbolic scattering.