Papers for Wednesday, Sep 28 2022

In galaxy clusters, the intracluster medium (ICM) is expected to host a diffuse, long-lived, and invisible population of "fossil" cosmic-ray electrons (CRe) with 1-100 MeV energies. These CRe, if "re-accelerated" by 100x in energy, can contribute synchrotron luminosity to cluster radio halos, relics, and phoenices. Re-acceleration may be aided by CRe scattering upon the ion-Larmor-scale waves that spawn when ICM is compressed, dilated, or sheared. We study CRe scattering and energy gain due to ion cyclotron (IC) waves generated by continuously-driven compression in 1D fully kinetic particle-in-cell simulations. We find that pitch-angle scattering of CRe by IC waves induces energy gain via magnetic pumping. In an optimal range of IC-resonant momenta, CRe may gain up to ~10-30% of their initial energy in one compress/dilate cycle with magnetic field amplification ~3-6x, assuming adiabatic decompression without further scattering and averaging over initial pitch angle.

Distinctive large-scale structures have been identified in the spatial distribution of optical galaxies up to redshift z ~ 1. In the more distant universe, the relationship between the dust-obscured population of star-forming galaxies observed at millimetre wavelengths and the network of cosmic filaments of dark matter apparent in all cosmological hydrodynamical simulations is still under study. Using the NIKA2 dual-band millimetre camera, we mapped a field of ~ 90 arcminutes^2 in the direction of the star GJ526 simultaneously in its 1.15-mm and 2.0-mm continuum wavebands to investigate the nature of the quasi-alignment of five sources found ten years earlier with the MAMBO camera at 1.2 mm. We find that these sources are not clumps of a circumstellar debris disc around this star as initially hypothesized. Rather, they must be dust-obscured star-forming galaxies, or sub-millimetre galaxies (SMGs), in the distant background. The new NIKA2 map at 1.15 mm reveals a total of seven SMGs distributed in projection on the sky along a filament-like structure crossing the whole observed field. Furthermore, we show that the NIKA2 and supplemental Herschel photometric data are compatible with a model of the spectral energy distributions (SEDs) of these sources when a common redshift of 2.5 and typical values of the dust parameters for SMGs are adopted. Hence, we speculate that these SMGs might be located in a filament of the distant `cosmic web'. The length of this candidate cosmic filament crossing the whole map is at least 4 cMpc (comoving), and the separations between sources are between 0.25 cMpc and 1.25 cMpc at this redshift, in line with expectations from cosmological simulations. Nonetheless, further observations to determine the precise spectroscopic redshifts of these sources are required to definitively support this hypothesis of SMGs embedded in a cosmic filament of dark matter.

We analyze the structure of 15 protocluster forming regions in the Milky Way using their $1.3$ mm continuum emission maps from the ALMA-IMF large program. The analysis of the clouds structure is performed using the delta-variance spectrum technique. The calculated spectra display a self-similar regime on small scales as well as the presence of a prominent bump on larger scales and whose physical size, $L_{\rm hub}$, falls in the range of $\approx 7000$ au to $60000$ au. These scales correspond to the sizes of the most compact clumps within the protocluster forming clouds. A significant correlation is found between $L_{\rm hub}$ and the surface density of the free-free emission estimated from the integrated flux of the H41$\alpha$ recombination line $\left(\Sigma_{\rm H41\alpha}^{\rm free-free}\right)$ as well as a significant anti-correlation between $L_{\rm hub}$ and the ratio of the 1.3 mm to 3 mm continuum emission fluxes $\left(S_{\rm 1.3 mm}^{\rm cloud}/S_{\rm 3 mm}^{\rm cloud}\right)$. Smaller values of $\left(S_{\rm 1.3 mm}^{\rm cloud}/S_{\rm 3 mm}^{\rm cloud}\right)$ and larger values of $\Sigma_{\rm H41\alpha}^{\rm free-free}$ correspond to more advanced evolutionary stages of the protocluster forming clumps. Hence, our results suggest that the sizes of the densest regions in the clouds are directly linked to their evolutionary stage and to their star formation activity with more evolved clouds having larger protocluster forming clumps. This is a an indication that gravity pays a vital role in regulating the size and mass growth of these clumps with ongoing gas accretion.

The ESA Euclid mission will provide high-quality imaging for about 1.5 billion galaxies. A software pipeline to automatically process and analyse such a huge amount of data in real time is being developed by the Science Ground Segment of the Euclid Consortium; this pipeline will include a model-fitting algorithm, which will provide photometric and morphological estimates of paramount importance for the core science goals of the mission and for legacy science. The Euclid Morphology Challenge is a comparative investigation of the performance of five model-fitting software packages on simulated Euclid data, aimed at providing the baseline to identify the best suited algorithm to be implemented in the pipeline. In this paper we describe the simulated data set, and we discuss the photometry results. A companion paper (Euclid Collaboration: Bretonni\`ere et al. 2022) is focused on the structural and morphological estimates. We created mock Euclid images simulating five fields of view of 0.48 deg2 each in the $I_E$ band of the VIS instrument, each with three realisations of galaxy profiles (single and double S\'ersic, and 'realistic' profiles obtained with a neural network); for one of the fields in the double S\'ersic realisation, we also simulated images for the three near-infrared $Y_E$, $J_E$ and $H_E$ bands of the NISP-P instrument, and five Rubin/LSST optical complementary bands ($u$, $g$, $r$, $i$, and $z$). To analyse the results we created diagnostic plots and defined ad-hoc metrics. Five model-fitting software packages (DeepLeGATo, Galapagos-2, Morfometryka, ProFit, and SourceXtractor++) were compared, all typically providing good results. (cut)

$\kappa$CrB is a ${\sim}2.5\,$Gyr old K1 sub-giant star, with an eccentric exo-Jupiter at ${\sim}2.8\,$au and a debris disc at tens of au. We present ALMA Band 6 ($1.3\,$mm) and HST scattered light ($0.6\,\mu$m) images, demonstrating $\kappa$CrB's broad debris disc, covering an extent $50{-}180\,$au in the millimetre (peaking at $110\,$au), and $51{-}280\,$au in scattered light (peaking at $73\,$au). By modelling the millimetre emission, we estimate the dust mass as ${\sim}0.016\,M{\oplus}$, and constrain lower-limit planetesimal sizes as $D_{\rm{max}}{>}1\,$km and the planetesimal belt mass as $M_{\rm{disc}}{>}1\,M_{\oplus}$. We constrain the properties of an outer body causing a linear trend in 17 years of radial velocity data to have a semi-major axis $8{-}66\,$au and a mass $(0.4{-}120)\,M_{\rm{Jup}}$. There is a large inner cavity seen in the millimetre emission, which we show is consistent with carving by such an outer massive companion with a string of lower mass planets. Our scattered light modelling shows that the dust must have a high anisotropic scattering factor ($g{\sim}0.8{-}0.9$) but an inclination ($i{\sim}30{-}40\,$degrees) that is inferred to be significantly lower than the $i{\sim}61\,$degrees millimetre inclination. The origin of such a discrepancy is unclear, but could be caused by a misalignment in the micron and millimetre sized dust. We place an upper limit on the CO gas mass of $M_{\rm{CO}}{<}(4.2{-}13){\times}10^{-7}\,M_{\oplus}$, and show this to be consistent with levels expected from planetesimal collisions, or from CO-ice sublimation as $\kappa$CrB begins its giant branch ascent.

As a candidate bound for the Binary Black Hole (BBH) merger events detected by LIGO/Virgo, Primordial Black Holes (PBHs) provide a useful tool to investigate the primordial curvature perturbations on small scales. Using the GWTC-1 to GWTC-3 catalogs, under the scenario that PBHs originate from large primordial curvature perturbations on small scales during inflationary epoch, we for the first time reconstruct the power spectrum of primordial curvature perturbations on small scales. It is found that the value of the amplitude of the primordial power spectrum is enhanced to $\mathcal{O}(10^{-2})$ on scales $\mathcal{O}(1)$ pc. This may imply the validity of PBH as a possible BBH merger candidate.

We present deep imaging of the ultra-diffuse Andromeda XIX dwarf galaxy from the Advance Camera for Surveys on the Hubble Space Telescope which resolves its stellar populations to below the oldest main sequence turn-off. We derive a full star formation history for the galaxy using MATCH, and find no evidence of star formation in the past 8 Gyr. We calculate a quenching time of $\tau_{90}=9.7\pm0.2$~Gyr, suggesting Andromeda~XIX ceased forming stars very early on. This early quenching, combined with its extremely large half-light radius, low density dark matter halo and lower than expected metallicity make it a unique galaxy within the Local Group and raises questions about how it formed. The early quenching time allows us to rule out feedback from bursty star formation as a means to explain its diffuse stellar population and low density dark matter halo. We find that the extended stellar population, low density halo and star formation could be explained by either tidal interactions (such as tidal shocking) or by late dry mergers, with the latter also explaining its low metallicity. Proper motions and detailed abundances would allow us to distinguish between these two scenarios.

We present an analysis of the optical and UV properties of AT 2020wey, a faint and fast tidal disruption event (TDE) at 124.3 Mpc. The light curve of the object peaked at an absolute magnitude of $M_{g} = -17.45$ mag and a maximum bolometric luminosity of $L_{\rm peak}=(8.74\pm0.69)\times10^{42}$ erg s$^{-1}$, making it comparably faint with iPTF16fnl, the faintest TDE to date. The time from the last non-detection to the $g$-band peak is 22.94 $\pm$ 2.03 days and the rise is well described by $L\propto t^{1.8}$. The decline of the bolometric light curve is described by a sharp exponential decay steeper than the canonical $t^{-5/3}$ power law, making AT 2020wey the fastest declining TDE to date. Multi-wavelength fits to the light curve indicate a complete disruption of a star of $M_*=0.11M_{\odot}$ by a black hole of $M_{\rm BH}=10^{6.46}M_{\odot}$. Our spectroscopic dataset reveals broad ($\sim10^{4}$ km s$^{-1}$) Balmer and He II $\lambda$4686 lines, with H$\alpha$ reaching its peak with a lag of $\sim8.2$ days compared to the continuum. In contrast to previous faint and fast TDEs, there are no obvious Bowen fluorescence lines in the spectra of AT 2020wey. There is a strong correlation between the MOSFIT-derived black hole masses of TDEs and their decline rate. However, AT 2020wey is an outlier in this correlation, which could indicate that its fast early decline may be dictated by a different physical mechanism than fallback. After performing a volumetric correction to a sample of 30 TDEs observed between 2018 and 2020, we conclude that faint TDEs are not rare by nature and that they should constitute up to $\sim$ 50 - 60 % of the entire population and their numbers could alleviate some of the tension between the observed and theoretical TDE rate estimates. We calculate the optical TDE luminosity function and we find a steep power-law relation $dN/dL_{g} \propto {L_{g}}^{-2.36}$.

We demonstrate that the James Webb Space Telescope (JWST) can detect infrared (IR) excess from the blended light spectral energy distribution of spatially unresolved terrestrial exoplanets orbiting nearby white dwarfs. We find that JWST is capable of detecting warm (habitable-zone; T$_{\rm eq}$=287 K) Earths or super-Earths and hot (400-1000 K) Mercury analogs in the blended light spectrum around the nearest 15 isolated white dwarfs with 10 hrs of integration per target using MIRI's Medium Resolution Spectrograph (MRS). Further, these observations constrain the presence of a CO$_2$-dominated atmosphere on these planets. The technique is nearly insensitive to system inclination, and thus observation of even a small sample of white dwarfs could place strong limits on the occurrence rates of warm terrestrial exoplanets around white dwarfs in the solar neighborhood. We find that JWST can also detect exceptionally cold (100-150 K) Jupiter-sized exoplanets via MIRI broadband imaging at $\lambda = 21\,\mathrm{\mu m}$ for the 34 nearest ($<13$ pc) solitary white dwarfs with 2 hrs of integration time per target. Using IR excess to detect thermal variations with orbital phase or spectral absorption features within the atmosphere, both of which are possible with long-baseline MRS observations, would confirm candidates as actual exoplanets. Assuming an Earth-like atmospheric composition, we find that the detection of the biosignature pair O$_3$+CH$_4$ is possible for all habitable-zone Earths (within 6.5 pc; six white dwarf systems) or super-Earths (within 10 pc; 17 systems) orbiting white dwarfs with only 5-36 hrs of integration using MIRI's Low Resolution Spectrometer (LRS).

We present the public data release of the Uchuu-UM galaxy catalogues by applying the \UM algorithm to assign galaxies to the dark matter halos in the Uchuu $N$-body cosmological simulation. It includes a variety of baryonic properties for all galaxies down to $\sim 5\times10^8 M_{\odot}$ with halos in a mass range of $10^{10}<M_{\rm halo}/M_{\odot}<5\times10^{15}$ up to redshift $z=10$. Uchuu-UM includes more than $10^{4}$ cluster-size halos in a volume of $\sim 25.7\,{\rm Gpc}^3$, and it reproduces observed stellar mass functions across the redshift range of $z=0-7$, galaxy quenched fractions, and clustering statistics at low redshifts. Compared to the previous largest UM catalogue, the Uchuu-UM catalogue includes significantly more massive galaxies hosted by large-mass dark matter halos. Overall, the number density profile of galaxies follows the dark matter profile, with the profile becoming steeper around the splashback radius and flattening at larger radii. The number density profile of galaxies tends to be steeper for larger stellar masses and depends on the color of galaxies, with red galaxies having steeper slopes at all radii compared with blue galaxies. The quenched fraction exhibits a strong dependence on the stellar mass and increases toward the inner regions of clusters. The publicly available Uchuu-UM galaxy catalogue presented here represents a step forward in modeling ongoing and upcoming large galaxy surveys.

The long-term evolution of the internal, strong magnetic fields of neutron stars needs a specific numerical modelling. The diversity of the observed phenomenology of neutron stars indicates that their magnetic topology is rather complex and three-dimensional simulations are required, for example, to explain the observed bursting mechanisms and the creation of surface hotspots. We present MATINS, a new three dimensions numerical code for magneto-thermal evolution in neutron stars, based on a finite-volume scheme that employs the cubed-sphere system of coordinates. In this first work, we focus on the crustal magnetic evolution, with the inclusion of realistic calculations for the neutron star structure, composition and electrical conductivity assuming a simple temperature evolution profile. MATINS follows the evolution of strong fields (1e14-1e15 Gauss) with complex non-axisymmetric topologies and dominant Hall-drift terms, and it is suitable for handling sharp current sheets. After introducing the technical description of our approach and some tests, we present long-term simulations of the non-linear field evolution in realistic neutron star crusts. The results show how the non-axisymmetric Hall cascade redistributes the energy over different spatial scales. Following the exploration of different initial topologies, we conclude that during a few tens of kyr, an equipartition of energy between the poloidal and toroidal components happens at small-scales. However, the magnetic field keeps a strong memory of the initial large-scales, which are much harder to be restructured or created. This indicates that large-scale configuration attained during the neutron star formation is crucial to determine the field topology at any evolution stage.

The ingot WFS is a new kind of wavefront sensor specifically designed to deal with the elongation of LGS reference sources on ELT-class telescopes. Like the pyramid, it belongs to the family of pupil plane wavefront sensors and can be considered as a generalization of the pyramid WFS for extended, three-dimensional elongated sources. The current design uses a simple, reflective roof-shaped prism to split the light into three pupils that are used to retrieve the wavefront shape. A test-bench has been realized at the INAF-Padova laboratories to test the alignment and functioning of the ingot. The bench is equipped with a deformable lens, conjugated to the pupil plane, able to apply low-order aberrations and with a hexapod for the precise alignment of the ingot prism. In this work we present a robust and fully automated Python-code alignment procedure, which is able, by using the optical feedback from the I-WFS, to adjust its 6-degrees of freedom. Moreover, we report on the tests conducted with the deformable lens to characterize the ingot WFS response to low-order aberrations in terms of sensitivity and linearity. The results are used as a comparison for simulations to validate the ray-tracing modeling approach with the future goal of optimizing the procedure adopted for signal calculation and phase retrieval.

The Ingot WFS represents an innovative and indispensable class of sensors conceived to overcome some limitations due to the LGSs geometry, which is significantly different from the point-object originated by a NGS. Here we overview the project, aiming at investigating the performance of an ELT-like telescope equipped with the Ingot WFS, facing different aspects of the program: the needs for numerical simulations and laboratory experiments, the prototype and, finally the future plan for the verification on sky.

The three-dimensional distribution of the Ly$\alpha$ forest has been extensively used to constrain cosmology through measurements of the baryon acoustic oscillations (BAO) scale. However, more cosmological information could be extracted from the full shapes of the Ly$\alpha$ forest correlations through the Alcock-Paczy\'nski (AP) effect. In this work, we prepare for a cosmological analysis of the full shape of the Ly$\alpha$ forest correlations by studying synthetic data of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We use a set of one hundred eBOSS synthetic data sets in order to validate such an analysis. These mocks undergo the same analysis process as the real data. We perform a full-shape analysis on the mean of the correlation functions measured from the one hundred eBOSS realizations, and find that our model of the Ly$\alpha$ correlations performs well on current data sets. We show that we are able to obtain an unbiased full-shape measurement of $D_M/D_H(z_\mathrm{eff})$, where $D_M$ is the transverse comoving distance, $D_H$ is the Hubble distance, and $z_\mathrm{eff}$ is the effective redshift of the measurement. We test the fit over a range of scales, and decide to use a minimum separation of $r_\mathrm{min}=25\ h^{-1}\text{Mpc}$. We also study and discuss the impact of the main contaminants affecting Ly$\alpha$ forest correlations, and give recommendations on how to perform such analysis with real data. While the final eBOSS Ly$\alpha$ BAO analysis measured $D_M/D_H(z_\mathrm{eff}=2.33)$ with $4\%$ statistical precision, a full-shape fit of the same correlations could provide a $\sim2\%$ measurement.

We aim to study spatial and kinematic substructures of the region encompassed by Upper Scorpius and Ophiuchus star-forming regions, determining dynamical traceback ages and studying the star formation history of the complex. We identified seven different groups in this region. Four groups (nu Sco, beta Sco, sigma Sco and delta Sco) are part of Upper Scorpius, two groups (rho Oph and alpha Sco) are in Ophiuchus, and another group (pi Sco) is a nearby young population. We found an age gradient from the rho Oph group (the youngest) to the delta Sco group (<~5 Myr), showing that star formation was a sequential process for the past 5 Myr. Our traceback analysis shows that Upper Scorpius and rho Oph groups share a common origin. The closer group of pi Sco is probably older, and the traceback analysis suggests that this group and the alpha Sco group have a different origin, likely related to other associations in the Sco-Cen complex. Our study shows that this region has a complex star formation history that goes beyond the current formation scenario, likely a result of stellar feedback from massive stars, supernovae explosions, and dynamic interactions between stellar groups and the molecular gas. In particular, we speculate that photo-ionisation from the massive delta Sco star could have triggered star formation first in the beta Sco group and then in the nu Sco group. The perturbations of stellar orbits due to stellar feedback and dynamical interactions could also be responsible for the 1-3 Myr difference that we found between dynamical traceback ages and isochronal ages.

Recent studies have shown that galaxies at cosmic noon are redder in the center and bluer in the outskirts, mirroring results in the local universe. These color gradients could be caused by either gradients in the stellar age or dust opacity; however, distinguishing between these two causes is impossible with rest-frame optical photometry alone. Here we investigate the underlying causes of the gradients from spatially-resolved rest-frame $U-V$ vs. $V-J$ color-color diagrams, measured from early observations with the James Webb Space Telescope. We use $1\, \mu m - 4\, \mu m$ NIRCam photometry from the CEERS survey of a sample of 54 galaxies with $M_* / M_\odot>10$ at redshifts $1.7<z<2.3$ selected from the 3D-HST catalog. We model the light profiles in the F115W, F200W and F356W NIRCam bands using \texttt{imcascade}, a Bayesian implementation of the Multi-Gaussian expansion (MGE) technique which flexibly represents galaxy profiles using a series of Gaussians. We construct resolved rest-frame $U-V$ and $V-J$ color profiles. The majority of star-forming galaxies have negative gradients (i.e. redder in the center, bluer in the outskirts) in both $U-V$ and $V-J$ colors consistent with radially decreasing dust attenuation. A smaller population (roughly 15\%) of star-forming galaxies have positive $U-V$ but negative $V-J$ gradients implying centrally concentrated star-formation. For quiescent galaxies we find a diversity of UVJ color profiles, with roughly one-third showing star-formation in their center. This study showcases the potential of JWST to study the resolved stellar populations of galaxies at cosmic noon.

The Gemini Planet Imager (GPI) is an integral field spectrograph (IFS) and coronagraph that is one of the few current generation instruments optimized for high-contrast direct imaging of substellar companions. The instrument is in the process of being upgraded and moved from its current mount on the Gemini South Observatory in Cerro Pachon, Chile, to its twin observatory, Gemini North, on Mauna Kea (a process colloquially dubbed 'GPI 2.0'). We present the designs that have been developed for the part of GPI 2.0 that pertains to upgrading various optical components of the GPI coronagraphic system. More specifically, we present new designs for the apodizer and Lyot stop (LS) that achieve better raw contrast at the inner working angle of the dark zone as well as improved core throughput while retaining a similar level of robustness to LS misalignment. To generate these upgraded designs, we use our own publicly available software package called APLC-Optimization that combines a commercial linear solver (Gurobi) with a high contrast imaging simulation package (HCIPy) in order to iteratively propagate light through a simulated model of an apodized phase lyot coronagraph (APLC), optimizing for the best coronagraph performance metrics. The designs have recently finished being lithographically printed by a commercial manufacturer and will be ready for use when GPI 2.0 goes on-sky in 2023.

Aims. With an observing time span of more than 20 years, the CORALIE radial-velocity survey is able to detect long-term trends in data corresponding to companions with masses and separations accessible to direct imaging. Combining exoplanet detection techniques such as radial-velocities from the CORALIE survey, astrometric accelerations from Hipparcos and Gaia eDR3, and relative astrometry from direct imaging, removes the degeneracy of unknown orbital parameters. This allows precise model-independent masses of detected companions to be derived, which provides a powerful tool to test models of stellar and substellar mass-luminosity relations. Methods. Long-term precise Doppler measurements with the CORALIE spectrograph reveal radial-velocity signatures of companions on long-period orbits. The long baseline of radial-velocity data allows the detectability of such companion candidates to be assessed with direct imaging. We combine long-period radial-velocity data with absolute astrometry from Hipparcos and Gaia eDR3 and relative astrometry derived from new direct imaging detections with VLT/SPHERE to fit orbital parameters and derive precise dynamical masses of these companions. Results. In this paper we report the discovery of new companions orbiting HD~142234, HD~143616, and HIP~22059, as well as the first direct detection of HD~92987~B, and update the dynamical masses of two previously directly imaged companions; HD~157338~B and HD~195010~B. The companions span a period range of 32 to 273 years and are all very low mass stellar companions, ranging from 225 to 477~$M_{\rm{Jup}}$. We compare the derived dynamical masses to mass-luminosity relations of very low mass stars (<0.5~$M_{\odot}$), and discuss the importance of using precursor radial-velocity and astrometric information to inform the future of high-contrast imaging of exoplanets and brown dwarfs

Young stellar populations provide a record of past star formation, and by establishing their members' dynamics and ages, it is possible to reconstruct the full history of star formation events. Gaia has greatly expanded the number of accessible stellar populations, with one of the most notable recently-discovered associations being Cepheus Far North (CFN), a population containing hundreds of members spanning over 100 pc. With its proximity (d $\lesssim$ 200 pc), apparent substructure, and relatively small population, CFN represents a manageable population to study in depth, with enough evidence of internal complexity to produce a compelling star formation story. Using Gaia astrometry and photometry combined with additional spectroscopic observations, we identify over 500 candidate CFN members spread across 7 subgroups. Combining ages from isochrones, asteroseismology, dynamics, and lithium depletion, we produce well-constrained ages for all seven subgroups, revealing a largely continuous 10 Myr star formation history in the association. By tracing back the present-day populations to the time of their formation, we identify two spatially and dynamically distinct nodes in which stars form, one associated with $\beta$ Cephei which shows mostly co-spatial formation, and one associated with EE Draconis with a more dispersed star formation history. This detailed view of star formation demonstrates the complexity of the star formation process, even in the smallest of regions.

We characterize the intrinsic properties of any FRB using its redshift $z$, spectral index $\alpha$ and energy $E_{33}$ in units of $10^{33} \, {\rm J}$ emitted across $2128 - 2848\; {\rm MHz}$ in the FRB's rest frame. Provided that $z$ is inferred from the measured extra-galactic dispersion measure $DM_{\rm Ex}$, the fluence $F$ of the observed event defines a track in $(\alpha, E_{33})$ space which we refer to as the "energy track". Here we consider the energy tracks for a sample of $254$ non-repeating low dispersion measure FRBs from the CHIME/FRB Catalog-1, and use these to determine $n(E_{33} \mid \alpha)$ the conditional energy distribution i.e. the number of FRBs in the interval $\Delta E_{33}$ given a value of $\alpha$. Considering $-10 \le \alpha \le 10$, we find that the entire energy scale shifts to higher energy values as $\alpha$ is increased. For all values of $\alpha$, we can identify two distinct energy ranges indicating that there are possibly two distinct FRB populations. At high energies, the distribution is well fitted by a modified Schechter function whose slope and characteristic energy both increase with $\alpha$. At low energies, the number of FRBs are in excess of the predictions of the modified Schechter function indicating that we may have a distinctly different population of low-energy FRBs. We have checked that our main findings are reasonably robust to the assumptions regarding the Galactic Halo and Host galaxy contributions to the dispersion measure.

We present the first comprehensive study of NaCl and KCl gases in Io's atmosphere in order to investigate their characteristics, and to infer properties of Io's volcanoes and subsurface magma chambers. In this work, we compile all past spectral line observations of NaCl and KCl in Io's atmosphere from the Atacama Large Millimeter/submillimeter Array (ALMA) and use atmospheric models to constrain the physical properties of the gases on several dates between 2012 and 2018. NaCl and KCl appear to be largely spatially confined and for observations with high spectral resolution, the temperatures are high (~500-1000 K), implying a volcanic origin. The ratio of NaCl:KCl was found to be ~5-6 in June 2015 and ~3.5-10 in June 2016, which is consistent with predictions based on observations of Io's extended atmosphere, and less than half the Na:K ratio in chondrites. Assuming these gases are volcanic in origin, these ratios imply a magma temperature of ~1300 K, such that the magma will preferentially outgas KCl over NaCl.

After 16 years of on-sky operation, Subaru Telescope's facility adaptive optics AO188 is getting several major upgrades to become the extreme-AO AO3000 (3000 actuators in the pupil compared to 188 previously). AO3000 will provide high-Strehl images for several instruments from visible to mid-infrared, notably the Infrared Camera and Spectrograph (IRCS), and the Subaru Coronagraphic Extreme Adaptive Optics (SCExAO). For this upgrade, the original 188-element deformable mirror (DM) will be replaced with ALPAO's $64\times64$ DM. The visible wavefront sensor will also be upgraded at a later date, but in the meantime we are adding a near-infrared Wavefront Sensor (NIR WFS), using either a double roof prism pyramid mode or a focal plane WFS mode. This new wavefront sensor will use for the first time First Light's C-RED ONE camera, allowing for a full control of the $64\times64$ DM at up to 1.6 kHz. One of the challenges is the use of non-destructive reads and a rolling shutter with the modulated pyramid. This upgrade will be particularly exciting for SCExAO, since the extreme-AO loop will focus more on creating high-contrast dark zones instead of correcting large atmospheric residuals. It will be the first time two extreme-AO loops will be combined on the same telescope. Finally, the setup AO3000+SCExAO+IRCS will serve as a perfect demonstrator for the Thirty Meter Telescope's Planetary Systems Imager (TMT-PSI). We will present here the design, integration and testing of AO3000, and show the first on-sky results.

We explore a re-parameterization of the lensing amplitude tension between weak lensing (WL) and cosmic microwave background (CMB) data and its implications for a joint resolution with the Hubble tension. Specifically, we focus on the lensing amplitude over a scale of 12 Mpc in absolute distance units using a derived parameter $S_{12}$ and show its constraints from recent surveys in comparison with Planck 2018. In WL alone, we find that the absolute distance convention correlates $S_{12}$ with $H_0$. Accounting for this correlation in the 3D space $S_{12}\times \omega_m \times h$ reproduces the usual levels of $2\sim 3\sigma$ tension inferred from $S_8\times\Omega_m$. Additionally, we derive scaling relations in the $S_8\times h$ and $S_{12}\times h$ planes that are allowed by $\Lambda$CDM and extrapolate target scalings needed to solve the $H_0$ and lensing-amplitude tensions jointly in a hypothetical beyond-$\Lambda$CDM model. As a test example, we quantify how the early dark energy scenario compares with these target scalings. Useful fitting formulae for $S_8$ and $S_{12}$ as a function of other cosmological parameters in $\Lambda$CDM are provided, with 1% precision.

The accretion disks of active galactic nuclei (AGN) are promising locations for the merger of compact objects detected by gravitational wave (GW) observatories. Embedded within a baryon-rich, high density environment, mergers within AGN are the only GW channel where an electromagnetic (EM) counterpart must occur (whether detectable or not). Considering AGN with unusual flaring activity observed by the Zwicky Transient Facility (ZTF), we describe a search for candidate EM counterparts to binary black hole (BBH) mergers detected by LIGO/Virgo in O3. After removing probable false positives, we find nine candidate counterparts to BBH mergers mergers during O3 (seven in O3a, two in O3b) with a $p$-value of 0.019. Based on ZTF sky coverage, AGN geometry, and merger geometry, we expect $\approx 3(N_{\rm BBH}/83)(f_{\rm AGN}/0.5)$ potentially detectable EM counterparts from O3, where $N_{\rm BBH}$ is the total number of observed BBH mergers and $f_{\rm AGN}$ is the fraction originating in AGN. Further modeling of breakout and flaring phenomena in AGN disks is required to reduce our false positive rate. Two of the events are also associated with mergers with total masses $> 100M_\odot$, which is the expected rate for O3 if hierarchical (large mass) mergers occur in the AGN channel. Candidate EM counterparts in future GW observing runs can be better constrained by coverage of the Southern sky as well as spectral monitoring of unusual AGN flaring events in LIGO/Virgo alert volumes. A future set of reliable AGN EM counterparts to BBH mergers will yield an independent means of measuring cosmic expansion ($H_0$) as a function of redshift.

We consider the potential for line intensity mapping (LIM) of the rotational CO(1-0), CO(2-1) and CO(3-2) transitions to detect deviations from General Relativity from $0 < z < 3$ within the framework of a very general class of modified gravity models, called Horndeski theories. Our forecast assumes a multi-tracer analysis separately obtaining information from the matter power spectrum and the first two multipoles of the redshift space distortion power spectrum. To achieve $\pm 0.1$ level constraints on the slope of the kinetic gravity braiding and Planck mass evolution parameters, a mm-wave LIM experiment would need to accumulate $\approx 10^8-10^9$ spectrometer hours, feasible with instruments that could be deployed in the 2030s. Such a measurement would constrain large portions of the remaining parameter space available to Scalar-Tensor modified gravity theories. Our modeling code is publicly available.

Although coronal mass ejections (CMEs) resembling flux ropes generally expand self-similarly, deformations along their fronts have been reported in observations and simulations. We present evidence of one CME becoming deformed after a period of self-similarly expansion in the corona. The event was observed by multiple white-light imagers on January 20-22, 2021. The change in shape is evident in observations from the heliospheric imagers from the Wide-Field Imager for Solar Probe Plus (WISPR), which observe this CME for $\sim$ 44 hours. We reconstruct the CME using forward-fitting models. In the first hours, observations are consistent with a self-similar expansion but later on the front flattens forming a dimple. Our interpretation is that the CME becomes deformed at $\sim0.1\ au$ due to differences in the background solar wind speeds. The CME expands more at higher latitudes, where the background solar wind is faster. We consider other possible causes for deformations, such as loss of coherence and slow-mode shocks. The CME deformation seems to cause a time-of-arrival error of 16 hours at $\sim 0.5\ au$. The deformation is clear only in the WISPR observations and, it thus, would have been missed by 1~AU coronagraphs. Such deformations may help explain the time-of-arrival errors in events where only coronagraph observations are available.

We present here for the first time a study of a peculiar giant radio galaxy (GRG) - GRG-J223301+131502 using deep multi-frequency radio observations from GMRT (323, 612, and 1300 MHz) and LOFAR (144 MHz) along with optical spectroscopic observations with the WHT 4.2m optical telescope. Our observations have firmly established its redshift of 0.09956 and unveiled its exceptional jet structure extending more than $\sim$200 kpc leading to a peculiar `kink' structure of $\sim$100 kpc. We measure the overall size of this GRG to be 1.83 Mpc, which exhibits lobes without any prominent hotspots and closely resembles a `Barbell'. Our deep low-frequency radio maps clearly reveal the steep-spectrum diffuse emission from the lobes of the GRG. The magnetic field strength of $\sim$ 5 $\mu$G and spectral ages between about 110 to 200 mega years for the radio lobes were estimated using radio data from LOFAR 144 MHz observations and GMRT 323 and 612 MHz observations. We discuss the possible causes leading to the formation of the observed `kink' feature for the GRG, which includes precession of the jet axis, development of instabilities and magnetic reconnection. Despite its enormous size, the Barbell GRG is found to be residing in a low mass (M$_{200} \sim 10^{14}$ $\rm M_{\odot}$) galaxy cluster. This GRG with two-sided, large-scale jets with a kink and diffuse outer lobes residing in a cluster environment, provides an opportunity of exploring the structure and growth of GRGs in different environments.

We systematically investigate the mergers of neutron star-white dwarf binaries from beginning to end, with focus on the properties of the inflows and outflows in accretion disks and their electromagnetic emissions. Using population synthesis models, we determine a subset of these binaries in which the white dwarf companion undergoes unstable mass transfer and complete tidal disruption, forming a large accretion disk around the neutron star. The material evolves according to an advection-dominated accretion disk model with nuclear burning, neutrino-emissions, and disk-surface wind ejection. The extreme dynamics of the entire process has proven difficult for analytic analysis, and thus currently the properties are poorly understood. The outflows from the mergers are iron- and nickel-rich, giving rise to optical and infrared emissions powered from the decay of the radioactive iron-type isotopes, calculated via the SuperNu light-curve code. We find these systems capable of powering bright, yet short-lived, optical transients with the potential to power gamma-ray bursts.

1H 0323+342 is optically one of the nearest and brightest very radio loud narrow line Seyfert 1 galaxies (vRL NLSy1). It is also one of the first vRLNLSy1s detected at gamma-ray energies by the Fermi-LAT. We report the results of monitoring the optical flux of 1H 0323+342 during more than six and a half years. In some cases, we, for the first time, simultaneously use two telescopes to monitor the optical flux of 1H 0323+342 on timescales ranging from minutes to hours, demonstrating the reality of low amplitude microvariability whole events with durations of a few hours for this object. Based on the present results, as well as those of earlier studies, we suggest that this represents a fundamental timescale associated with the underlying source region. We also present an enhancement of Howell's comparison star method for detecting microvariability.

We use multi-object near-infrared (NIR) spectroscopy with VLT/KMOS to investigate the role of the environment in the evolution of the ionized gas properties of narrow-band selected H$\alpha$ emitters (HAEs) in the Spiderweb protocluster at $z=2.16$. Based on rest-frame optical emission lines, H$\alpha$ and [NII]$\lambda$6584, we confirm the cluster membership of 39 of our targets (i.e. 93% success rate), and measure their star-formation rates (SFR), gas-phase oxygen abundances and effective radius. We parametrize the environment where our targets reside by using local and global density indicators based on previous samples of spectroscopic and narrow-band cluster members. We find that star-forming galaxies embedded in the Spiderweb protocluster display SFRs compatible with those of the main sequence and morphologies comparable to those of late-type galaxies at $z=2.2$ in the field. We also report a mild gas-phase metallicity enhancement ($0.6\pm0.3$ dex) at intermediate stellar masses. Furthermore, we identify two UVJ-selected quiescent galaxies with residual H$\alpha$-based star formation and find signs of extreme dust obscuration in a small sample of SMGs based on their FIR and H$\alpha$ emission. Interestingly, the spatial distribution of these objects differs from the rest of HAEs, avoiding the protocluster core. Finally, we explore the gas fraction-gas metallicity diagram for 7 galaxies with molecular gas masses measured by ATCA using CO(1-0). In the context of the gas-regulator model, our objects are consistent with relatively low mass-loading factors, suggesting lower outflow activity than field samples at the cosmic noon and thus, hinting at the onset of environmental effects in this massive protocluster.

Reliable localization of high-energy transients, such as Gamma-ray Bursts (GRBs) and Soft Gamma-ray Repeaters (SGRs), is the prerequisite for characterizing the burst properties (e.g. spectrum) and implementing the follow-up observations in the multi-wavelength and multi-messenger. Localization based on the relative counts of different detectors has been widely used for all-sky gamma-ray monitors, such as {\it CGRO}/BATSE, {\it Fermi}/GBM, POLAR, and GECAM. There are two major statistical frameworks for counts distribution localization methods: $\chi^{2}$ and Bayesian. Here, we studied and compared several localization methods based on these two statistical frameworks, by directly checking the localization probability map and credible region with comprehensive simulations. We find that the Bayesian method is generally more applicable for various bursts than $\chi^{2}$ method. We also proposed a location-spectrum iteration approach based on the Bayesian inference, which could not only take advantage of straightforward calculation but also alleviate the problems caused by the spectral difference between the burst and location templates.

We present photometric observations and numerical simulations of 2016 SD$_{106}$, a low inclination ($i=4.8^{\circ}$) extreme trans-Neptunian Object with a large semi-major axis ($a=350$ au) and perihelion ($q= 42.6$ au). This object possesses a peculiar neutral color of $g-r = 0.45\pm0.05$ and $g-i=0.72\pm0.06$, in comparison with other distant trans-Neptunian objects, all of which have moderate-red to ultra-red colors. A numerical integration based on orbital fitting on astrometric data covering eight years of arc confirms that 2016 SD$_{106}$ is a metastable object without significant scattering evolution. Each of the clones survived at the end of the 1 Gyr simulation. However, very few neutral objects with inclinations $<5^{\circ}$ have been found in the outer solar system, even in the main Kuiper belt. Furthermore, most mechanisms which lift perihelion distances are expected to produce a very low number of extreme objects with inclinations $<5^{\circ}$. We thus explored the possibility that a hypothetical distant planet could increase the production of such objects. Our simulations show that no 2016 SD$_{106}$-like orbits can be produced from three Kuiper belt populations tested (i.e. plutinos, twotinos, and Haumea Family) without the presence of an hypothetical planet, while a few similar orbits can be obtained with it; however, the presence of the additional planet produces a wide range of large semimajor-axis / large perihelion objects, in apparent contradiction with the observed scarcity of objects in those regions of phase space. Future studies may determine if there is a connection between the existence of a perihelion gap and a particular orbital configuration of an hypothetical distant planet.

We study the relationship between Galactic location ($R, Z$) and photometric activity for 3.6 million M dwarf stars within 1 kpc of the Sun. For this purpose, we identify 906 unique flare events as a proxy for magnetic activity from the SkyMapper Southern Survey DR3. We adopt vertical distance $|Z|$ from the Galactic disc as a proxy for age and confirm a strong trend of flaring fraction decreasing with growing stellar age. Among M dwarfs within 50 pc of the Sun, we find a flaring fraction of 1-in-1,500, independent of spectral type from M2 to M7, suggesting that these stars are all in a flare-saturated young evolutionary stage. We find a hint of a kink in the slope of the overall flare fraction near 100 pc from the plane, where a steep decline begins; this slope change is visible for mid-type M dwarfs (M3--M5), suggesting it is not an artefact of mixing spectral type. Together with SDSS H$\alpha$ emission, this trend is additional evidence that the activity fraction of M dwarfs depends on Galactic height and activity lifetime. While there is a hint of flattening of the overall activity fraction above $|Z|\approx$ 500 pc, our data do not constrain this further. Within $\sim$500 pc distance from the Sun, we find no sign of radial disk gradients in flare activity, which may only be revealed by samples covering a larger radial range.

We present mid-infrared spectroscopic observations of the nucleus of the nearby Seyfert galaxy NGC 7469 taken with the MIRI instrument on the James Webb Space Telescope (JWST) as part of Directors Discretionary Time Early Release Science (ERS) program 1328. The high resolution nuclear spectrum contains 19 emission lines covering a wide range of ionization. The high ionization lines show broad, blueshifted emission reaching velocities up to 1700 km s$^{-1}$ and FWHM ranging from $\sim500 - 1100$ km s$^{-1}$. The width of the broad emission and the broad to narrow line flux ratios correlate with ionization potential. The results suggest a decelerating, stratified, AGN driven outflow emerging from the nucleus. The estimated mass outflow rate is one to two orders of magnitude larger than the current black hole accretion rate needed to power the AGN. Eight pure rotational H$_{2}$ emission lines are detected with intrinsic widths ranging from FWHM $\sim 125-330$ km s$^{-1}$. We estimate a total mass of warm H$_{2}$ gas of $\sim1.2\times10^{7}$M$_{\odot}$ in the central 100 pc. The PAH features are extremely weak in the nuclear spectrum, but a $6.2\mu$m PAH feature with an equivalent width $\sim0.07\mu$m and a flux of $2.7\times10^{-17}$ W m$^{-2}$ is detected. The spectrum is steeply rising in the mid-infrared, with a silicate strength $\sim0.02$, significantly smaller than seen in most PG QSOs, but comparable to other Seyfert 1's. These early MIRI mid-infrared IFU data highlight the power of JWST to probe the multi-phase interstellar media surrounding actively accreting supermassive black holes.

We present the first quantitative spectral analysis of blue supergiant stars in the nearby galaxy NGC 2403. Out of a sample of 47 targets observed with the LRIS spectrograph at the Keck I telescope we have extracted 16 B- and A-type supergiants for which we have data of sufficient quality to carry out a comparison with model spectra of evolved massive stars and infer the stellar parameters. The radial metallicity gradient of NGC 2403 that we derive has a slope of -0.14 (+/- 0.05) dex/r_e, and is in accordance with the analysis of H II region oxygen abundances. We present evidence that the stellar metallicities that we obtain in extragalactic systems in general agree with the nebular abundances based on the analysis of the auroral lines, over more than one order of magnitude in metallicity. Adopting the known relation between stellar parameters and intrinsic luminosity we find a distance modulus m-M = 27.38 +/- 0.08 mag. While this can be brought into agreement with Cepheid-based determinations, it is 0.14 mag short of the value measured from the tip of the red giant branch. We update the mass-metallicity relation secured from chemical abundance studies of stars in resolved star-forming galaxies.

We examine the temporary evolution of axisymmetric magnetospheres around rapidly rotating black holes (BHs), by applying our two-dimensional particle-in-cell simulation code. Assuming a stellar-mass BH, we find that the created pairs fail to screen the electric field along the magnetic field, provided that the mass accretion rate is much small compared to the Eddington limit. Magnetic islands are created by reconnection near the equator and migrate toward the event horizon, expelling magnetic flux tubes from the BH vicinity during a large fraction of time. When the magnetic islands stick to the horizon due to redshift and virtually vanish, a strong magnetic field penetrates the horizon, enabling efficient extraction of energy from the BH. During this flaring phase, a BH gap appears around the inner light surface with a strong meridional return current toward the equator within the ergosphere. If the mass accretion rate is 0.025 percent of the Eddington limit, the BH's spin-down luminosity becomes 16-19 times greater than its analytical estimate during the flares, although its long-term average is only 6 percent of it. We demonstrate that the extracted energy flux concentrates along the magnetic field lines threading the horizon in the middle latitudes. It is implied that this meridional concentration of the Poynting flux may result in the formation of limb-brightened jets from low-accreting BH systems.

In this work, I present a qualitative discussion on the prospect of production of ultra-high photons in blazars. The sources are a subclass of active galactic nuclei which host supermassive black holes and fire relativistic jets into the intergalactic medium. The kpc scale jets are believed to be dominated by Poynting flux and constitute one of the most efficient cosmic particle accelerators, that potentially are capable of accelerating the particles up to EeV energies. Recent IceCube detection of astrophysical neutrino emission in coincidence with the enhanced gamma-ray from Tev blazar TXS 0506 + 056 further supports hadronic models of blazar emissions in which particle acceleration processes such as relativistic shocks, magnetic re-connection, and relativistic turbulence could energize hardrons, e. g. protons, up to energies equivalent to billions of Lorentz factors. The ensuing photo-pionic processes may then result in gamma-rays accompanied by neutrino flux. Furthermore, the fact that blazars are the dominant source of observed TeV emission encourages us to search for signatures of acceleration scenarios that would lead to the creation of ultra-high energy photons.

In recent years, scientific CMOS (sCMOS) sensors have found increasing applications to X-ray detection, including X-ray astronomical observations. In order to examine the performance of sCMOS sensors, we have developed X-ray cameras based on sCMOS sensors. Two cameras, CNX22 and CNX 66, have been developed using sCMOS sensors with a photosensitive area of 2 cm * 2 cm and 6 cm * 6 cm, respectively. The designs of the cameras are presented in this paper. The CNX22 camera has a frame rate of 48 fps, whereas CNX66 has a frame rate of currently 20 fps, that can be boosted to 100 fps in the future. The operating temperature of the sCMOS sensor can reach to -20C for CNX22 and -30C for CNX66 with a peltier cooler device. In addition to the commonly used mode of saving original images, the cameras provide a mode of real-time extraction of X-ray events and storage their information, which significantly reduces the requirement for data storage and offline analysis work. For both cameras, the energy resolutions can reach less than 200 eV at 5.9 keV using single-pixel events. These cameras are suitable for X-ray spectroscopy applications in laboratories and calibration for the space X-ray telescopes.

One of the exciting results of the last years concerning the study of the large-scale structure is the discovery of long (3-5 Mpc) bridges of radio emission connecting pairs of interacting clusters. This is the first direct evidence of the existence of particle acceleration and magnetic field amplification mechanisms outside galaxy clusters. Non-thermal components spread over such vast extents probe the dynamics of large-scale structures and the mechanisms of energy dissipation therein. In my talk, I will discuss recent results on radio bridges, showing how their observation represents a step forward in the search of the radio signature of the magnetized cosmic web.

Gaia DR3 and further releases have the potential to identify and categorise new emission-line stars in the Galaxy. We perform a comprehensive validation of astrophysical parameters from Gaia DR3 with the spectroscopically estimated emission-line star parameters from LAMOST OBA emission catalogue. We compare different astrophysical parameters provided by Gaia DR3 with those estimated using LAMOST spectra. By using a larger sample of emission-line stars, we perform a global polynomial and piece-wise linear fit to update the empirical relation to convert Gaia DR3 pseudo-equivalent width to observed equivalent width, after removing the weak emitters from the analysis. We find that the emission-line source classifications given by DR3 is in reasonable agreement with the classification from LAMOST OBA emission catalogue. The astrophysical parameters estimated by esphs module from Gaia DR3 provides a better estimate when compared to gspphot and gspspec. A second-degree polynomial relation is provided along with piece-wise linear fit parameters for the equivalent width conversion. We notice that the LAMOST stars with weak H{\alpha} emission are not identified to be in emission from BP/RP spectra. This suggests that emission-line sources identified by Gaia DR3 is incomplete. In addition, Gaia DR3 provides valuable information about the binary and variable nature of a sample of emission-line stars.

Variability is a commonly observed property of classical Be stars (CBe) stars. In extreme cases, complete disappearance of the H{\alpha} emission line occurs, indicating a disc-less state in CBe stars. The disc-loss and reappearing phases can be identified by studying the H{\alpha} line profiles of CBe stars on a regular basis. In this paper, we present the study of a set of selected 9 bright CBe stars, in the wavelength range of 6200 - 6700 {\AA}, to better understand their disc transient nature through continuous monitoring of their H{\alpha} line profile variations for 5 consecutive years (2015 -- 2019). Based on our observations, we suggest that 4 of the program stars (HD 4180, HD 142926, HD 164447 and HD 171780) are possibly undergoing disc-loss episodes, whereas one other star (HD 23302) might be passing through disc formation phase. The remaining 4 stars (HD 237056, HD 33357, HD 38708 and HD 60855) have shown signs of hosting a stable disc in recent epochs. Through visual inspection of the overall variation observed in the H{\alpha} EW for these stars, we classified them into groups of growing, stable and dissipating discs, respectively. Moreover, our comparative analysis using the BeSS database points out that the star HD 60855 has passed through a disc-less episode in 2008, with its disc formation happening probably over a timescale of only 2 months, between January and March 2008.

Context: We use phase curves of small bodies to measure absolute magnitudes and, together with complementary theoretical and laboratory results, to understand their surfaces' micro and macroscopic properties. Although we can observe asteroids up to phase angles of about 30 deg, the range of phase angles covered by outer solar system objects usually does not go further than 7 to 10 deg for centaurs and 2 deg for trans-Neptunian objects, and a linear relation between magnitude and phase angle may be assumed. Aims: We aim at directly comparing data taken for objects in the inner solar system (inside the orbit of Jupiter) with data of centaurs and trans-Neptunian objects. Methods: We use the SLOAN Moving Objects Catalog data to construct phase curves restricted to phase angles less than or equal to 7.5 deg, compatible with the angles observed for the trans-Neptunian/Centaur population. We assume a linear model for the photometric behavior to obtain absolute magnitudes and phase coefficients in the ugirz, V, and R filters. Results: We obtained absolute magnitudes in seven filters for $>4000$ objects. Our comparison with outer solar system objects points to a common property of the surfaces: intrinsically redder objects become blue with increasing phase angle, while the opposite happens for intrinsically bluer objects.

Over its lifespan, the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will monitor millions of supernovae (SNe) from explosion to oblivion, yielding an unprecedented ugrizy photometric dataset on their late-time evolution. Here, we show that the photometric evolution of Type Ibc SNe can be used to constrain numerous properties of their ejecta, without the need for expensive spectroscopic observations. Using radiative-transfer simulations for explosions of He-star progenitors of different initial masses, we show that the g-band filter follows primarily the strength of the FeII emission, the r-band [OI]6300-6364A and [NII]6548-6583A, the i-band [CaII]7291,7323A, and the z-band the CaII NIR triplet, and hence provides information on nucleosynthetic yields. Information on weaker lines, which may be used, for example, to constrain clumping, is absent. However, this deficiency may eventually be cured by improving the physical realism of radiative-transfer simulations through a closer connection to physically consistent 3D explosion models, and by the judicial selection of a much smaller set of spectral observations. Degeneracies inherent to the SN radiation will affect the interpretation of photometric measures, but line fluxes from nebular-phase spectra are similarly compromised. Importantly, our ``family'' of Type Ibc SN models follows a distinct trajectory in color-color magnitude diagrams as the ejecta evolve from 100 to 450d, allowing one to disentangle different progenitors or explosions. This photometric procedure provides a promising approach to study statistical samples of SNe Ibc and to confront them to ever improving progenitor and explosion models, to capture the onset of late-time interaction with circumstellar material, or to identify events currently unknown.

In this paper, we use three different kinds of observational data, including 130 strong gravitational lensing (SGL) systems, type Ia supernovae (SNeIa: Pantheon and Union2.1) and 31 Hubble parameter data points ($H(z)$) from cosmic chronometers to constrain the phenomenological model ($\rho_x\varpropto\rho_m a^{\xi}$). By combining these three kinds of data (Union2.1+SGL+$H(z)$), we get the parameter value at the confidence interval of $2\sigma$, $\Omega_{X,0} = 0.69\pm0.34$, $\omega_x = -1.24\pm0.61$, $\xi = 3.8\pm3.9$ and $H_0 = 70.22\pm0.86$ kms$^{-1}$Mpc$^{-1}$. According to our results, we find that the $\Lambda$CDM model is still the model which is in best agreement with the observational data at present, and the coincidence problem is not alleviated. In addition, the $\Omega_X$ and $\Omega_m$ have the same order of magnitude in $0<z<1.26$. At last, we obtain the transition redshift $z_T=0.645$. If the transition occurs in $z>0.645$, it is necessary to introduce the dark energy interacting with dark matter.

Pulsar halos are extended gamma-ray structures generated by electrons and positrons escaping from pulsar wind nebulae (PWNe), considered a new class of gamma-ray sources. They are ideal indicators for cosmic-ray propagation in localized regions of the Galaxy and particle escape process from PWNe. The cosmic-ray diffusion coefficient inferred from pulsar halos is more than two orders of magnitude smaller than the average value in the Galaxy, which has been arousing extensive discussion. We review the recent advances in the study of pulsar halos, including the characteristics of this class of sources, the known pulsar halos, the possible mechanisms of the extremely slow diffusion, the critical roles of pulsar halos in the studies of cosmic-ray propagation and electron injection from PWNe, and the implications on the problems of the cosmic positron excess and the diffuse TeV gamma-ray excess. Finally, we give prospects for the study in this direction based on the expectation of a larger sample of pulsar halos and deeper observations for bright sources.

Using VLBI to measure a so-called core shift effect is a common way of obtaining estimates of the jet magnetic field strength. The VLBI core is typically identified as the bright feature at the jet's base, and the position of the core changes with the observed frequency, $r_\mathrm{core} \propto \nu^{-1/k_r}$. In this work, we investigated the time variability of the core-shift effect in the blazar 3C 454.3. We employed self-referencing analysis of multi-frequency (5, 8, 15, 22-24, and 43 GHz) VLBA data covering 19 epochs from 2005 until 2010. We found significant core shift variability ranging from 0.27 to 0.86 mas between 5 and 43 GHz, confirming the core-shift variability phenomenon observed before. Time variability of the core-shift index ($k_r$) was found typically below one, with an average value of $0.85 \pm 0.08$ and a standard deviation of $0.30$. $k_r<1$ values were found during flaring and quiescent states and our results indicate that commonly assumed conical jet shape and equipartition conditions do not always hold simultaneously. Still, these conditions are often assumed when deriving magnetic field strengths from core shift measurements, leading to unreliable results if $k_r$ significantly deviates from unity. Therefore, it is important to verify that $k_r = 1$ holds before using core shift values and the equipartition assumption to derive physical parameters in the jets. When $k_r = 1$ epochs are selected in the case of 3C 454.3, the magnetic field estimates are indeed quite consistent, even though the core shift varies with time. Additionally, our estimations of the jet's magnetic flux in 3C 454.3 show that the source is indeed in the magnetically arrested disk state. Finally, we found a good correlation of the core position with the core flux density, $r_\mathrm{core}\propto S_\mathrm{core}^{0.7}$, which is consistent with increased particle density during the flares.

(Abridged) In this work, we study the 2.8x2.6 deg2 region in the emblematic Rosette Nebula, centred at the young cluster NGC 2244, with the aim of constructing the most reliable candidate member list to date, determining various structural and kinematic parameters, and learning about the past and the future of the region. Starting from a catalogue containing optical to mid-infrared photometry, as well as positions and proper motions from Gaia EDR3, we apply the Probabilistic Random Forest algorithm and derive membership probability for each source. Based on the list of almost 3000 probable members, of which about a third are concentrated within the radius of 20' from the centre of NGC 2244, we identify various clustered sources and stellar concentrations, and estimate the average distance of 1489+-37 pc (entire region), 1440+-32 pc (NGC 2244) and 1525+-36 pc (NGC 2237). The masses, extinction, and ages are derived by SED fitting, and the internal dynamic is assessed via proper motions relative to the mean proper motion of NGC 2244. NGC 2244 is showing a clear expansion pattern, with an expansion velocity that increases with radius. Its IMF is well represented by two power laws (dN/dM\propto M^{-\alpha}), with slopes \alpha = 1.05+-0.02 for the mass range 0.2 - 1.5 MSun, and \alpha = 2.3+-0.3 for the mass range 1.5 - 20 MSun, in agreement with other star forming regions. The mean age of the region is ~2 Myr. We find evidence for the difference in ages between NGC 2244 and the region associated with the molecular cloud, which appears slightly younger. The velocity dispersion of NGC 2244 is well above the virial velocity dispersion derived from the total mass (1000+-70 MSun) and half-mass radius (3.4+-0.2 pc). From the comparison to other clusters and to numerical simulations, we conclude that NGC 2244 may be unbound, and possibly even formed in a super-virial state.

Yukawa interactions can mediate relatively long-range attractive forces between fermions in the early universe. Such a globally attractive interaction creates an instability that can result in the growth of structure in the affected species even during the radiation dominated era. The formation and collapse of fermionic microhalos can create hot fireballs at the sites of the collapsing halos which inject energy into the cosmic plasma. We show that the injected energy can be partially converted into primordial magnetic fields and we estimate the correlation scale and the power spectrum of these fields. We show that they may be the seeds of the observed astrophysical magnetic fields.

To understand the Active Galactic Nuclei (AGN) phenomenon and their impact on the evolution of galaxies, a complete AGN census is required; however, finding heavily obscured AGNs is observationally challenging. Here we use the deep and extensive multi-wavelength data in the COSMOS field to select a complete sample of 578 infrared (IR) quasars ($L_{\rm AGN,IR}>10^{45}\rm \: erg\: s^{-1}$) at $z<3$, with minimal obscuration bias, using detailed UV-to-far IR spectral energy distribution (SED) fitting. We complement our SED constraints with X-ray and radio observations to further investigate the properties of the sample. Overall, 322 of the IR quasars are detected by Chandra and have individual X-ray spectral constraints. From a combination of X-ray stacking and $L_{\rm 2-10\rm keV}$ - $L_{\rm 6\: \mu m}$ analyses, we show that the majority of the X-ray faint and undetected quasars are heavily obscured (many are likely Compton thick), highlighting the effectiveness of the mid-IR band to find obscured AGNs. We find that 355 ($\approx$61%) IR quasars are obscured ($N_{\rm H}>10^{22}\rm \: cm^{-2}$) and identify differences in the average properties between the obscured and unobscured quasars: (1) obscured quasars have star-formation rates $\approx 3$ times higher than unobscured systems for no significant difference in stellar mass and (2) obscured quasars have stronger radio emission than unobscured systems, with a radio-loudness parameter $\approx 0.2 \rm \: dex$ higher. These results are inconsistent with a simple orientation model but in general agreement with either extreme host-galaxy obscuration towards the obscured quasars or a scenario where obscured quasars are an early phase in the evolution of quasars.

We present the pulsar_spectra software repository, an open-source pulsar flux density catalogue and automated spectral fitting software that finds the best spectral model and produces publication-quality plots. The Python-based software includes features that enable users in the astronomical community to add newly published spectral measurements to the catalogue as they become available. The spectral fitting software is an implementation of the method described in Jankowski et al. (2018) which uses robust statistical methods to decide on the best-fitting model for individual pulsar spectra. pulsar_spectra is motivated by the need for a centralised repository for pulsar flux density measurements to make published measurements more accessible to the astronomical community and provide a suite of tools for measuring spectra.

SVOM is a Sino-French space mission targeting high-energy transient astrophysical objects such as gamma-ray bursts. The soft X-ray part of the spectrum is covered by the Micro-channel X-ray Telescope (MXT) which is a narrow field telescope designed to precisely localize X-ray sources. This paper presents the method implemented on-board to characterize and localize X-ray sources with the MXT. A specific localization method was developed to accommodate the optical system of the MXT, which is based on "Lobster-Eye" grazing incidence micro-pore optics. For the first time, the algorithm takes advantage of cross-correlation techniques to achieve a localization accuracy down to 2 arcmin with less than 200 photons, which guarantees a rapid follow-up for most of the gamma-ray bursts that SVOM will observe. In this paper, we also study the limitations of the algorithm and characterize its performance.

The Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect is a net torque caused by solar radiation directly reflected and thermally re-emitted from the surface of small asteroids and is considered to be crucial in their dynamical evolution. By long-term photometric observations of selected near-Earth asteroids, it's hoped to enlarge asteroid samples with a detected YORP effect to facilitate the development of a theoretical framework. Archived light-curve data are collected and photometric observations are made for (1685) Toro and (85989) 1999 JD6, which enables measurement of their YORP effect by inverting the light curve to fit observations from a convex shape model. For (1685) Toro, a YORP acceleration $\upsilon=(3.2\pm0.3)\times10^{-9}\ \rm{rad\cdot d^{-2}}(1\sigma\ error)$ is updated, which is consistent with previous YORP detection based on different light-curve data; for (85989) 1999 JD6, it is determined that the sidereal period is $7.667749\pm 0.000009$ h, the rotation pole direction locates is at $\lambda=232\pm 2^{\circ},\ \beta = -59\pm 1^{\circ}$, the acceleration is detected to be $\upsilon = (2.4\pm0.3)\times10^{-8}\ \rm{rad\cdot d^{-2}}(1\sigma\ error)$ and in addition to obtaining an excellent agreement between the observations and model. YORP should produce both spin-up and spin-down cases. However, including (85989) 1999 JD6, the $\rm{d}\omega/\rm{d}t$ values of eleven near-Earth asteroids are positive totally, which suggests that there is either a bias in the sample of YORP detections or a real feature needs to be explained.

Aims. With the large amount of molecular emission data from (sub)millimeter observatories and incoming James Webb Space Telescope infrared spectroscopy, access to fast forward models of the chemical composition of protoplanetary disks is of paramount importance. Methods. We used a thermo-chemical modeling code to generate a diverse population of protoplanetary disk models. We trained a K-nearest neighbors (KNN) regressor to instantly predict the chemistry of other disk models. Results. We show that it is possible to accurately reproduce chemistry using just a small subset of physical conditions, thanks to correlations between the local physical conditions in adopted protoplanetary disk models. We discuss the uncertainties and limitations of this method. Conclusions. The proposed method can be used for Bayesian fitting of the line emission data to retrieve disk properties from observations. We present a pipeline for reproducing the same approach on other disk chemical model sets.

We report the characterisation of a multi-planetary system composed of five exoplanets orbiting the K-dwarf HD~23472 (TOI-174). In addition to the two super-Earths that were previously confirmed, we confirm and characterise three Earth-size planets in the system using ESPRESSO radial velocity observations. The planets of this compact system have periods of $P_d \sim 3.98\,$, $P_e \sim 7.90\,$, $P_f \sim 12.16\,$, $P_b \sim 17.67,\,$ and $P_c \sim 29.80\,$days and radii of $R_d \sim 0.75\,$ , $R_e \sim 0.82,$, $R_f \sim 1.13\,$, $R_b \sim 2.01,\,$ and, $R_c \sim 1.85\,$ \REarth. Because of its small size, its proximity to planet d's transit, and close resonance with planet d, planet e was only recently found. The planetary masses were estimated to be $M_d =0.54\pm0.22$, $M_e =0.76\pm0.30$, $M_f =0.64_{-0.39}^{+0.46}$, $M_b = 8.42_{-0.84}^{+0.83}$, and $M_c = 3.37_{-0.87}^{+0.92}$ \MEarth. These planets are among the lightest planets, with masses measured using the radial velocity method, demonstrating the very high precision of the ESPRESSO spectrograph. We estimated the composition of the system's five planets and found that their gas and water mass fractions increase with stellar distance, suggesting that the system was shaped by irradiation. The high density of the two inner planets ($\rho_d = 7.5_{-3.1}^{+3.9}$ and $\rho_e = 7.5_{-3.0}^{+3.9}\, \mathrm{g.cm^{-3}}$) indicates that they are likely to be super-Mercuries. This is supported by the modelling of the internal structures of the planets, which also suggests that the three outermost planets have significant water or gas content If the existence of two super-Mercuries in the system is confirmed, this system will be the only one known to feature two super-Mercuries, making it an excellent testing bed for theories of super-Mercuries formation. (Abridged)

The impact of feedback from Active Galactic Nuclei (AGN) on the cosmological evolution of the large scale structure is a long studied problem. However, it is still not well understood how the feedback energy couples to the ambient medium to influence the properties of AGN host galaxies and dark matter halos. In this work we investigate different modes of AGN feedback and their effect on the surrounding medium by probing the diffuse X-ray emission from the hot gas inside galaxy groups and clusters. For this purpose, we use the cosmological hydrodynamic simulation SIMBA to theoretically calculate the X-ray emission from simulated galaxy clusters/groups with the help of the Astrophysical Plasma Emission Code (APEC). We also perform synthetic observations of these systems with the Chandra X-ray telescope using the ray-tracing simulator Model of AXAF Response to X-rays (MARX). Our results show that in addition to the radiative wind mode of feedback from the AGN, jet and X-ray mode of feedback play significant roles in suppressing the X-ray emission from the diffuse gas in the vicinity of the black hole. Our mock observational maps suggest that the signatures of AGN feedback from the high redshift objects may not be detected with the instrumental resolution of current X-ray telescopes like Chandra, but provide promising prospect for detection of these features with potential X-ray missions such as Lynx.

We compute the radiative ro-vibrational emission spectrum of H2 involving quasibound states via a simple numerical method of resolution of the Schr\"odinger equation by introducing a modifed effective molecular potential. The comparison of the eigenvalues obtained with our approximation and other theoretical methods based on scattering resonance properties is excellent. Electric quadrupole and magnetic dipole contributions are calculated and we confirm the previous computations of Forrey of the electric quadrupole transition Einstein coeffcients. The astrophysical relevance of such quasibound levels is emphasized

The Pierre Auger Observatory is the world's largest detector for observation of ultra-high-energy cosmic rays (UHECRs) (above the energy of $10^{17}$ eV). It consists of a Fluorescence Detector (FD) and an array of particle detectors known as the Surface Detector (SD). Observations of extensive air showers by the Observatory can be used to probe hadronic interactions at high energy, in a kinematic and energy region inaccessible to experiments at man-made accelerators and to measure the muon component of the shower. Air showers induced by different primaries have different muon contents. With increasing mass of the primary cosmic ray particle, it is expected that the muon content in the corresponding air showers should also increase. Recent results obtained from the Pierre Auger Observatory and other experiments indicate that all the shower simulations underestimate the number of muons in the showers compared to the data. This is the so-called muon deficit. In this paper we briefly review the muon measurements, and present in more detail recent results on fluctuations in the muon number. These results provide new insights into the origin of the muon deficit in air shower simulations and constrain the models of hadronic interactions at ultrahigh energies. With the current design of the surface detectors it is also difficult to reliably separate the contributions of muons to the SD signal from the contributions of photons, electrons, and positrons. Therefore, we also present a new method to extract the muon component of the signal time traces recorded by each SD station using recurrent neural networks. The combination of such algorithms, with the future data collected by the upgraded Pierre Auger Observatory, will be a major step forward, as we are likely to achieve an unprecedented resolution in mass estimation on an event-by-event basis.

The nature of dark matter (DM) is still under intense debate. Sub-galactic scales are particularly critical, as different, currently viable DM models make diverse predictions on the expected abundance and density profile of DM haloes on these scales. We investigate the ability of sub-galactic DM haloes to act as strong lenses on background compact sources, producing gravitational lensing events on milli-arcsecond scales (milli-lenses), for different DM models. For each DM scenario, we explore whether a sample of $\sim$ 5000 distant sources is sufficient to detect at least one milli-lens. We develop a semi-analytical model to estimate the milli-lensing optical depth as a function of the source's redshift for various DM models. We employ the Press-Schechter formalism, as well as results from recent N-body simulations to compute the halo mass function, taking into account the appropriate spherically averaged density profile of haloes for each DM model. We treat the lensing system as a point-mass lens and invoke the effective surface mass density threshold to calculate the fraction of a halo that acts as a gravitational lens. We study three classes of dark matter models: cold DM, warm DM, and self-interacting DM. We find that haloes consisting of warm DM turn out to be optically thin for strong gravitational milli-lensing (zero expected lensing events). CDM haloes may produce lensing events depending on the steepness of the concentration-mass relation. Self-interacting DM haloes can efficiently act as gravitational milli-lenses only if haloes experience gravothermal collapse, resulting in highly dense central cores.

Simultaneous analysis of the C_2 and CN molecular bands in the 5100-5200 and 7930-8100A spectral regions is a promising alternative for the accurate determination of the carbon (C) and nitrogen (N) abundance in the atmospheres of the solar-like stars. Practical implementation of this new method became possible after recent improvements of the molecular constants for both molecules. The new molecular data predicted the correct line strength and line positions; therefore, they were included in the Vienna Atomic Line Database (VALD), which is widely used by astronomers and spectroscopists. In this paper, we demonstrate that the molecular data analysis provides C and, in particular, N abundances consistent with those derived from the atomic lines. We illustrate this by performing the analysis for three stars. Our results provide strong arguments for using the combination of C_2 and CN molecular lines for accurate nitrogen abundance determination keeping in mind the difficulties of using the N1 lines in the observed spectra of the solar-like stars.

We present an updated measurement of orbital period evolution of LMXB XTE J1710-281 by using eclipse timing technique. Using data obtained with XMM-Newton, Suzaku, RXTE, Chandra and AstroSat observatories, we report 21 new measurements of X-ray mid-eclipse times. We have discovered a third orbital period glitch in XTE J1710-281 with an F-test false alarm probability of ~0.7% for occurrence of the third glitch and report detection of four distinct epochs of orbital period in this system. This work presents a more robust estimation of occurrence of the second orbital period glitch. However, the epoch of occurrence of the third glitch is poorly constrained, between MJD 55726 to 56402. We have put lower limits of 1.48 ms, 0.97 ms and 0.45 ms, on sudden changes in orbital period between the successive epochs. We discuss the implications of our findings in context of magnetic nature of the companion star and possible scattering events with circum-binary objects around this binary system.

Observations of the South Polar Residual Cap suggest a possible erosion of the cap, leading to an increase of the global mass of the atmosphere. We test this assumption by making the first comparison between Viking 1 and InSight surface pressure data, which were recorded 40 years apart. Such a comparison also allows us to determine changes in the dynamics of the seasonal ice caps between these two periods. To do so, we first had to recalibrate the InSight pressure data because of their unexpected sensitivity to the sensor temperature. Then, we had to design a procedure to compare distant pressure measurements. We propose two surface pressure interpolation methods at the local and global scale to do the comparison. The comparison of Viking and InSight seasonal surface pressure variations does not show changes larger than +-8 Pa in the CO2 cycle. Such conclusions are supported by an analysis of Mars Science Laboratory (MSL) pressure data. Further comparisons with images of the south seasonal cap taken by the Viking 2 orbiter and MARCI camera do not display significant changes in the dynamics of this cap over a 40 year period. Only a possible larger extension of the North Cap after the global storm of MY 34 is observed, but the physical mechanisms behind this anomaly are not well determined. Finally, the first comparison of MSL and InSight pressure data suggests a pressure deficit at Gale crater during southern summer, possibly resulting from a large presence of dust suspended within the crater.

The general consensus is that LTGs undergo intense star-formation activity, while ETGs are mostly inactive. We question this general rule and investigate the existence of star-forming ETGs and quiescent LTGs in the local Universe. By computing the physical properties of 2,209 such galaxies in the GAMA survey being morphologically classified and using information on their structural properties as well as the density of their local environment, we seek for understanding the differences from their 'typical' counterparts. We separate galaxies into subsets based on their dominant ionising process, making use of criteria based on the WH$_{\alpha}$ width and the [NII/H$_{\alpha}$] ratio. Taking advantage of the SED fitting code CIGALE we derive galaxy properties, such as the $M_\text{star}$, $M_\text{dust}$, and SFR and also estimate the unattenuated and the dust-absorbed stellar emission, for both the young and old stellar populations. Ongoing star-formation activity is found in 47% of ETGs and 8% of LTGs are quiescent. The star-forming E galaxies, together with the LBSs, constitute a population that follows very well the SFMS of spiral galaxies. The fraction of the luminosity originating from young stars in the star-forming ETGs is quite substantial ($\sim$ 25%) and similar to that of the star-forming LTGs. Investigating possible differences between star-forming and quiescent galaxies we find that the intrinsic shape of the SED of the star-forming galaxies is, on average, very similar for all morphological types. Concerning their structural parameters, quiescent galaxies tend to show larger values of the S\'ersic index and larger $R_\text{eff}$ (compared to star-forming galaxies). Finally, we find that star-forming galaxies preferably reside in lower-density environments compared to the quiescent ones, which exhibit a higher percentage of sources being members of groups.

We present a complete halo model for mixed dark matter composed of cold dark matter (CDM) and ultralight axion-like particles (ULAs). Our model treats ULAs as a biased tracer of CDM, in analogy to treatments of massive neutrinos and neutral hydrogen. The model accounts for clustering of ULAs around CDM host halos, and fully models the cross correlations of both components. The model inputs include the ULA Jeans scale, and soliton density profile. Our model can be used to predict the matter power spectrum, $P(k)$, on non-linear scales for sub-populations of ULAs across the mass range $10^{-33}\text{ eV}\leq m\leq 10^{-21}\text{ eV}$, and can be calibrated against future mixed DM simulations to improve its accuracy. The mixed DM halo model also allows us to assess the importance of various approximations.

Spectroscopic measurements of iron abundances are prone to systematic modelling errors. We present 3D non-LTE calculations across 32 STAGGER-grid models with effective temperatures from 5000 K to 6500 K, surface gravities of 4.0 dex and 4.5 dex, and metallicities from $-$3 dex to 0 dex, and study the effects on 171 Fe I and 12 Fe II optical lines. In warm metal-poor stars, the 3D non-LTE abundances are up to 0.5 dex larger than 1D LTE abundances inferred from Fe I lines of intermediate excitation potential. In contrast, the 3D non-LTE abundances can be 0.2 dex smaller in cool metal-poor stars when using Fe I lines of low excitation potential. The corresponding abundance differences between 3D non-LTE and 1D non-LTE are generally less severe but can still reach $\pm$0.2 dex. For Fe II lines the 3D abundances range from up to 0.15 dex larger, to 0.10 dex smaller, than 1D abundances, with negligible departures from 3D LTE except for the warmest stars at the lowest metallicities. The results were used to correct 1D LTE abundances of the Sun and Procyon (HD 61421), and of the metal-poor stars HD 84937 and HD 140283, using an interpolation routine based on neural networks. The 3D non-LTE models achieve an improved ionisation balance in all four stars. In the two metal-poor stars, they remove excitation imbalances that amount to 250 K to 300 K errors in effective temperature. For Procyon, the 3D non-LTE models suggest [Fe/H] = 0.11 $\pm$ 0.03, which is significantly larger than literature values based on simpler models. We make the 3D non-LTE interpolation routine for FG-type dwarfs publicly available, in addition to 1D non-LTE departure coefficients for standard MARCS models of FGKM-type dwarfs and giants. These tools, together with an extended 3D LTE grid for Fe II from 2019, can help improve the accuracy of stellar parameter and iron abundance determinations for late-type stars.

In the quest for high-energy neutrino sources, the Astrophysical Multimessenger Observatory Net- work (AMON) has implemented a new search by combining data from the High Altitude Water Cherenkov (HAWC) observatory and the Astronomy with a Neutrino Telescope and Abyss environ- mental RESearch (ANTARES) neutrino telescope. Using the same analysis strategy as in a previous detector combination of HAWC and IceCube data, we perform a search for coincidences in HAWC and ANTARES events that are below the threshold for sending public alerts in each individual detector. Data were collected between July 2015 and February 2020 with a livetime of 4.39 years. Over this time period, 3 coincident events with an estimated false-alarm rate of $< 1$ coincidence per year were found. This number is consistent with background expectations.

Flux variability is a remarkable property of black hole (BH) accreting systems, and a powerful tool to investigate the multi-scale structure of the accretion flow. The X-ray band is where some of the most rapid variations occur, pointing to an origin in the innermost regions close to the BH. The study of fast time variability provides us with means to explore the accretion flow around compact objects in ways which are inaccessible via spectral analysis alone, and to peek at regions which cannot be imaged with the currently available instrumentation. In this chapter we will discuss fast X-ray variability in stellar-mass BH systems, namely binary systems containing a star and a BH, occasionally contrasting it with observations of supermassive BHs in active galactic nuclei. We will explore how rapid variations of the X-ray flux have been used in multiple studies as a diagnostic of the innermost regions of the accretion flow in these systems. To this aim we will provide an overview of the currently most used analysis approaches for the study of X-ray variability, describe observations of both aperiodic and quasi-periodic phenomena, and discuss some of the proposed models.

Next-generation wide-field spectroscopic surveys will observe the infall regions around large numbers of galaxy clusters with high sampling rates for the first time. Here we assess the feasibility of extracting the large-scale cosmic web around clusters using forthcoming observations, given realistic observational constraints. We use a sample of 324 hydrodynamic zoom-in simulations of massive galaxy clusters from TheThreeHundred project to create a mock-observational catalogue spanning $5R_{200}$ around 160 analogue clusters. These analogues are matched in mass to the 16 clusters targetted by the forthcoming WEAVE Wide-Field Cluster Survey (WWFCS). We consider the effects of the fibre allocation algorithm on our sampling completeness and find that we successfully allocate targets to 81.7 $\% \pm$ 1.3 of the members in the cluster outskirts. We next test the robustness of the filament extraction algorithm by using a metric, $D_{\text{skel}}$, which quantifies the distance to the filament spine. We find that the median positional offset between reference and recovered filament networks is $D_{\text{skel}} = 0.13 \pm 0.02$ Mpc, much smaller than the typical filament radius of $\sim$ 1 Mpc. Cluster connectivity of the recovered network is not substantially affected. Our findings give confidence that the WWFCS will be able to reliably trace cosmic web filaments in the vicinity around massive clusters, forming the basis of environmental studies into the effects of pre-processing on galaxy evolution.

Standard solar models (SSMs) constructed in accordance with old solar abundances are in reasonable agreement with seismically inferred results, but SSMs with new low-metal abundances disagree with the seismically inferred results. The constraints of neutrino fluxes on solar models exist in parallel with those of helioseismic results. The solar neutrino fluxes were updated by Borexino Collaboration. We constructed rotating solar models with new low-metal abundances where the effects of enhanced diffusion and convection overshoot were included. A rotating model using OPAL opacities and the Caffau abundance scale has better sound-speed and density profiles than the SSM with the old solar abundances and reproduces the observed $p$-mode frequency ratios $r_{02}$ and $r_{13}$. The depth and helium abundance of the convection zone of the model agree with the seismically inferred ones at the level of $1\sigma$. The updated neutrino fluxes are also reproduced by the model at the level of $1\sigma$. The effects of rotation and enhanced diffusion not only improve the model's sound-speed and density profiles but bring the neutrino fluxes predicted by the model into agreement with the detected ones. Moreover, the calculations show that OP may underestimate opacities for the regions of the Sun with $T\gtrsim5\times10^{6}$ K by around $1.5\%$, while OPAL may underestimate opacities for the regions of the Sun with $2\times10^{6}$ K $\lesssim T \lesssim 5\times10^{6}$ K by about $1-2\%$.

We report a detection of a 0.153 days period in the classical nova V1674 Her using the NICER observations taken within a month since the outburst (i.e., about MJD 59405). The X-ray period is consistent with the orbital period previously found in the optical band, strongly suggesting the NICER signal as the X-ray orbital periodicity of the system. A seemingly double-humped profile was obtained by folding the detrended X-ray light curve with the period after removing the rotational X-ray pulsations of the nova. The profile may be caused by occultation by the companion or the accretion disk, possibly indicating a high inclination of the system. The gamma-ray emission of V1674 Her with a significance level > 5 sigma was detected by Fermi-LAT close to its optical peak and the emission faded away within 1 day, which is the shortest duration known for a gamma-ray nova. Folded on 0.153days, a marginal gamma-ray variability can be also seen in the LAT light curve, but without the double-hump feature observed in X-rays. If the gamma-ray modulation is real, its origin is probably different from that observed in the X-ray and optical bands.

The orbital decay of massive black holes in galaxies in the aftermath of mergers is at the heart of whether massive black holes successfully pair and merge, leading to emission of low-frequency gravitational waves. The role of dynamical friction sourced from the gas distribution has been uncertain because many analytical and numerical studies have either focused on a homogeneous medium or have not reached resolutions below the scales relevant to the problem, namely the Bondi-Hoyle-Lyttleton radius. We perform numerical simulations of a massive black hole moving in a turbulent medium in order to study dynamical friction from turbulent gas. We find that the black hole slows down to the sound speed, rather than the turbulent speed, and that the orbital decay is well captured if the Bondi-Hoyle-Lyttleton radius is resolved with at least five resolution elements. We find that the larger the turbulent eddies, the larger the scatter in dynamical friction magnitude, because of the stochastic nature of the problem, and also of the larger over- and under-densities encountered by the black hole along its trajectory. Compared to the classic solution in a homogeneous medium, the magnitude of the force depends more weakly on the Mach number, and dynamical friction is overall more efficient for high Mach numbers, but less efficient towards and at the transonic regime.

A gravitational wave background from a first order phase transition in the early universe may be observable at millihertz gravitational wave (GW) detectors such as the Laser Interferometer Space Antenna (LISA). In this paper we introduce and test a method for investigating LISA's sensitivity to gravitational waves from a first order phase transition using parametrised templates as an approximation to a more complete physical model. The motivation for developing the method is to provide a less computationally intensive way to perform Markov Chain Monte Carlo (MCMC) inference on the thermodynamic parameters of a first order phase transition, or on generally computationally intensive models. Starting from a map between the physical parameters and the parameters of an empirical template, we first construct a prior on the empirical parameters that contains the necessary information about the physical parameters; we then use the inverse mapping to reconstruct approximate posteriors on the physical parameters from a fast MCMC on the empirical template. We test the method on a double broken power law approximation to spectra in the sound shell model. The reconstruction method substantially reduces the proposal evaluation time, and despite requiring some precomputing of the mapping, this method is still cost-effective overall. In two test cases, with signal-to-noise $\sim 40$, the method recovers the physical parameters and the spectrum of the injected gravitational wave power spectrum to $95\%$ confidence. In previous Fisher matrix analysis we found the phase boundary speed $v_{\rm w}$ was expected to be the best constrained of the thermodynamic parameters. In this work, for an injected phase transition GW power spectrum with $v_{\rm w} = 0.55$, with a direct sample on the thermodynamic parameters we recover $0.630^{+0.17}_{-0.059}$ and for our reconstructed sample $0.646^{+0.098}_{-0.075}$.

We conducted three-dimensional hydrodynamical simulations of common envelope evolution (CEE) of a neutron star (NS) or a black hole (BH) inside a red supergiant (RSG) envelope and find that the jets that we expect the NS/BH to launch during the CEE spin-up the common envelope. We find that when the NS/BH launches jets that are exactly perpendicular to the orbital plane (the jets are aligned with the orbital angular momentum) the jets deposit angular momentum to the envelope that is aligned with the orbital angular momentum. When the jets' axis is inclined to the orbital angular momentum axis so is the angular momentum that the jets deposit to the envelope. Such tilted jets might be launched when the NS/BH has a close companion when it enters the RSG envelope. We did not allow for spiralling-in and could follow the evolution for only three orbits. The first orbit mimics the plunge-in phase of the CEE, when the NS/BH rapidly dives in, while the third orbit mimics the self-regulated phase when spiralling-in is very slow. We find that the jets deposit significant amount of angular momentum only during the plunge-in phase. A post-CEE core collapse supernova explosion will leave two NS/BH, bound or unbound, whose spin might be misaligned to the orbital angular momentum. Our results strengthen an earlier claim that inclined-triple-star CEE might lead to spin-orbit misalignment of NS/BH-NS/BH binary systems.

The $\gamma-$ray emission from flat-spectrum radio quasars (FSRQs) is thought to be dominated by the inverse Compton scattering of the external sources of photon fields, e.g., accretion disk, broad-line region (BLR), and torus. FSRQs show strong optical emission lines and hence can be a useful probe of the variability in BLR output, which is the reprocessed disk emission. We study the connection between the optical continuum, H$\gamma$ line, and $\gamma-$ray emissions from the FSRQ PKS~1222+216, using long-term ($\sim$2011-2018) optical spectroscopic data from Steward Observatory and $\gamma-$ray observations from $Fermi$-LAT. We measured the continuum ($F_{C,opt}$) and H$\gamma$ ($F_{H\gamma}$) fluxes by performing a systematic analysis of the 6029-6452 \r{A} optical spectra. We observed stronger variability in $F_{C,opt}$ than $F_{H\gamma}$, an inverse correlation between H$\gamma$ equivalent width and $F_{C,opt}$, and a redder-when-brighter trend. Using discrete cross-correlation analysis, we found a positive correlation (DCF$\sim$0.5) between $F_{\gamma-ray>100MeV}$ and $F_{C,opt}$ (6024-6092 \r{A}) light curves with time-lag consistent with zero at 2$\sigma$ level. We found no correlation between $F_{\gamma-ray>100MeV}$ and $F_{H\gamma}$ light curves, probably dismissing the disk contribution to the optical and $\gamma$-ray variability. The observed strong variability in the $Fermi$-LAT flux and $F_{\gamma-ray>100MeV}-F_{C,opt}$ correlation could be due to the changes in the particle acceleration at various epochs. We derived the optical-to-$\gamma$-ray spectral energy distributions (SEDs) during the $\gamma$-ray flaring and quiescent epochs that show a dominant disk component with no variability. Our study suggests that the $\gamma$-ray emission zone is likely located at the edge of the BLR or in the radiation field of the torus.

The Event Horizon Telescope (EHT) observation unveiled the first image of supermassive black hole Sgr A* showing a shadow of diameter $\theta_{sh}= 48.7 \pm 7\,\mu$as with fractional deviation from the Schwarzschild black hole shadow diameter $\delta = -0.08^{+0.09}_{-0.09}~\text{(VLTI)},-0.04^{+0.09}_{-0.10}~\text{(Keck)}$. The Sgr A* shadow size is within $~10\%$ of the Kerr predictions, providing us with another tool to investigate the nature of strong-field gravity. We use the Sgr A* shadow observables to constrain metrics of four independent and well-motivated, parametrically different from Kerr spacetime, rotating regular spacetimes, and the corresponding no-horizon spacetimes. We present constraints on the deviation parameter $g$ of rotating regular black holes. The shadow angular diameter $\theta_{sh}$ within $1 \sigma$ region, places bounds on the parameters $a$ and $g$. Together with EHT bounds on $\theta_{sh}$ and $\delta$ of Sgr A*, our analysis concludes that the three rotating regular black holes, viz., Bardeen Hayward, and Simpson-Visser black holes, and corresponding no-horizon spacetimes agree with the EHT results of Sgr A*. Thus, these three rotating regular spacetimes and Kerr black holes are indiscernible in some parameter space, and one can not rule out the possibility of the former being strong candidates for astrophysical black holes.

The local detection rate of interstellar objects can allow for estimations of the total number of similar objects bound by the Milky Way thin disk. If interstellar objects of artificial origin are discovered, the estimated total number of objects can be lower by a factor of about $10^{16}$ if they target the habitable zone around the Sun. We propose a model for calculating the quantity of natural or artificial interstellar objects of interest based on the object's velocity and observed density. We then apply the model to the case of chemically propelled rockets from extraterrestrial civilizations. Finally, we apply the model to three previously discovered interstellar objects -- the object 'Oumuamua of unknown origin and the first interstellar meteors CNEOS 2014-01-08 and CNEOS 2017-03-09.

The Event Horizon Telescope (EHT) collaboration's image of the compact object at the galactic centre is the first direct evidence of the supermassive black hole Sgr A$^*$. The shadow of Sgr A$^*$ has an angular diameter $d_{sh}= 48.7 \pm 7\,\mu$as with fractional deviation from the Schwarzschild black hole shadow diameter $\delta= -0.08^{+0.09}_{-0.09}\,,-0.04^{+0.09}_{-0.10}$ (for the VLTI and Keck mass-to-distance ratios). Sgr A$^*$'s shadow size is within $~10\%$ of the Kerr predictions, equipping us with yet another tool to analyse the gravity in the strong-field regime, including testing loop quantum gravity (LQG). We use Sgr A$^*$'s shadow to constrain the metrics of two well-motivated LQG-inspired rotating black holes (LIRBHs) models characterized by an additional deviation parameter $L_q$, which recover the Kerr spacetime in the absence of quantum effects ($L_q \to 0$). We use the astrophysical observables shadow area $A$ and oblateness $D$ to estimate the black hole parameters. When increasing the size of the quantum effects through $L_q$, the black hole shadow size increases monotonically, while the shape gets more distorted, allowing us to constrain the fundamental parameter $L_q$. While the EHT observational results completely rule out the wormhole region in the second LIRBH, a substantial parameter region of the generic black holes in both models agree with the EHT results. We find upper bounds on $L_q$ from the shadow of Sgr A$^*$: $L_q \lesssim 0.0423$ and $L_q \lesssim 0.0821$ for the two LIRBHs respectively, both more stringent than those obtained with the EHT image of M87$^*$.

In the presence of a background magnetic field, axion dark matter induces an electric field and can thus excite phonon-polaritons in suitable materials. We revisit the calculation of the axion-photon conversion power output from such materials, accounting for finite volume effects, and material losses. Our calculation shows how phonon-polaritons can be converted to propagating photons at the material boundary, offering a route to detecting the signal. Using the dielectric functions of GaAs, Al$_2$O$_3$, and SiO$_2$, a fit to our loss model leads to a signal of lower magnitude than previous calculations. We demonstrate how knowledge of resonances in the dielectric function can directly be used to calculate the sensitivity of any material to axion dark matter. We argue that a combination of low losses encountered at $\mathcal{O}(1)$ K temperatures and near future improvements in detector dark count allow one to probe the QCD axion in the mass range $m_a\approx 100$ meV. This provides further impetus to examine novel materials and further develop detectors in the THz regime. We also discuss possible tuning methods to scan the axion mass.

We study the dynamics of domain wall solitons in $(2+1)d$ field theories. These objects are extended along one of the spatial directions, so they also behave as strings; hence the name of domain wall strings. We show analytically and numerically that the amount of radiation from the propagation of wiggles on these objects is negligible except for regions of high curvature. Therefore, at low curvatures, the domain wall strings behave exactly as the Nambu-Goto action predicts. We show this explicitly with the use of several different numerical experiments of the evolution of these objects in a lattice. We then explore their dynamics in the presence of internal mode excitations. We do this again by performing field theory simulations and identify an effective action that captures the relevant interactions between the different degrees of freedom living on the string. We uncover a new parametric resonance instability that transfers energy from the internal mode to the position of the domain wall. We show that this instability accelerates the radiation of the internal mode energy. We also explore the possibility of exciting the internal mode of the soliton with the collision of wiggles on the domain wall. Our numerical experiments indicate that this does not happen unless the wiggles have already a wavelength of the order of the string thickness. Finally, we comment on the possible relevance of our findings to cosmological networks of defects. We argue that our results cast some doubts on the significance of the internal modes in cosmological applications beyond a brief transient period right after their formation. This, however, should be further investigated using cosmological simulations of our model.

Realistic estimations on the elastic properties of neutron star matter are carried out with a large strain ($\varepsilon \lesssim 0.5$) in the framework of relativistic-mean-field model with Thomas-Fermi approximation, where various crystalline configurations are considered in a fully three-dimensional geometry with reflection symmetry. Our calculation confirms the validity of assuming Coulomb crystals for the droplet phase above neutron drip density, which nonetheless does not work at large densities since the elastic constants are found to be decreasing after reaching their peaks. Similarly, the analytic formulae derived in the incompressible liquid-drop model gives excellent description for the rod phase at small densities, which overestimates the elastic constants at larger densities. For slabs, due to the negligence on the variations of their thicknesses, the analytic formulae from liquid-drop model agree qualitatively but not quantitatively with our numerical estimations. By fitting to the numerical results, these analytic formulae are improved by introducing dampening factors. The impacts of nuclear symmetry energy are examined adopting two parameter sets, corresponding to the slope of symmetry energy $L = 41.34$ and 89.39 MeV. Even with the uncertainties caused by the anisotropy in polycrystallines, the elastic properties of neutron star matter obtained with $L = 41.34$ and 89.39 MeV are distinctively different, results in detectable differences in various neutron star activities.

Next-generation terrestrial gravitational-wave observatories will detect $\mathcal{O}(10^{5})$ signals from compact binary coalescences every year. These signals can last for several hours in the detectors' sensitivity band and they will be affected by multiple unresolved sources contributing to a confusion-noise background in the data. Using an information-matrix formalism, we estimate the impact of the confusion noise power spectral density in broadening the parameter estimates of a GW170817-like event. If our estimate of the confusion noise power spectral density is neglected, we find that masses, spins, and distance are biased in about half of our simulations under realistic circumstances. The sky localization, while still precise, can be biased in up to $80\%$ of our simulations, potentially posing a problem in follow-up searches for electromagnetic counterparts.

This paper studies the oscillation properties of relativistic, non-self-gravitating tori in the background of a distorted deformed compact object. This work concentrates on the static and axially symmetric metric containing two quadrupole parameters; relating to the central object and the external fields. This metric may associate the observable effects to these parameters as dynamical degrees of freedom. The astrophysical motivation for choosing such a field is the possibility of constituting a reasonable model for an actual scenario occurring in the vicinity of compact objects. This paper aims to investigate the radial epicyclic frequency in a perfect fluid disk and not a test particle scenario via a local analysis. To achieve this goal, we employ the vertically integrated technique to able to treat the equation analytically. The tori are also modelled with Keplerian and non-Keplerian distributions of specific angular momentum, and we discuss the dependence of oscillation properties on the variable of the model related to angular momentum distribution and quadrupoles. In the present contribution, we further explore these properties with the possibility of relating oscillatory frequencies to some high-frequency quasi-periodic oscillations models and observed data.

Equatorial Plasma Bubbles (EPBs) are plumes of low density plasma that rise up from the bottomside of the F layer towards the exosphere. EPBs are known causes of radio wave scintillations which can degrade communications with spacecraft. We build a random forest regressor to predict and forecast the probability of an EPB [0-1] detected by the IBI processor on-board the SWARM spacecraft. We use 8-years of Swarm data from 2014 to 2021 and transform the data from a time series into a 5 dimensional space consisting of latitude, longitude, mlt, year, and day-of-the-year. We also add Kp, F10.7cm and solar wind speed. The observations of EPBs with respect to geolocation, local time, season and solar activity mostly agrees with existing work, whilst the link geomagnetic activity is less clear. The prediction has an accuracy of 88% and performs well across the EPB specific spatiotemporal scales. This proves that the XGBoost method is able to successfully capture the climatological and daily variability of SWARM EPBs. Capturing the daily variance has long evaded researchers because of local and stochastic features within the ionosphere. We take advantage of Shapley Values to explain the model and to gain insight into the physics of EPBs. We find that as the solar wind speed increases the probability of an EPB decreases. We also identify a spike in EPB probability around the Earth-Sun perihelion. Both of these insights were derived directly from the XGBoost and Shapley technique.

The impact of axion-like particles on the light polarization around the horizon of suppermassive black hole (SMBH) is discussed in the light of the latest polarization measurement of the Event Horizon Telescope (EHT). We investigate different sources of the polarization due to axion interaction with photons and the magnetic field of SMBH. These can modify the linear and circular polarization parameters of the emitted light. We have shown that a significant circular polarization can be produced via the photon scattering from the background magnetic field with axions as off-shell particles. This can further constrain the parameter space of ultralight axion-like particles and their couplings with photons. The future precise measurements of circular polarization can probe the features of ultralight axions in the near vicinity of SMBH.