Papers for Thursday, May 26 2022

Gaia Data Release 3 will provide the astronomical community with the largest stellar spectroscopic survey to date ($>$ 100 million sources). The low resolution (R$\sim$50) blue photometer (BP) and red photometer (RP) spectra will allow for the estimation of stellar parameters such as effective temperature, surface gravity and metallicity. We create mock Gaia BP/RP spectra and use Fisher information matrices to probe the resolution limit of stellar parameter measurements using BP/RP spectra. The best-case scenario uncertainties that this analysis provides are then used to produce a mock observed stellar population in order to probe the false positive rate (FPR) of identifying extremely metal-poor (EMP) stars. We find that we can confidently identify metal-poor stars at magnitudes brighter than $G = 16$. At fainter magnitudes true detections start to be overwhelmed by false positives. When adopting the commonly-used $G < 14$ magnitude limit for metal-poor star searches, we find a FPR for the low-metallicity regimes [Fe/H] < -2, -2.5 and -3 of just 14$\%$, 33$\%$ and 56$\%$ respectively, offering the potential for a significant improvement on previous targeting campaigns. Additionally, we explore the chemical sensitivity obtainable directly from BP/RP spectra for Carbon and $\alpha$-elements. We find an absolute Carbon abundance uncertainty of $\sigma_{A(C)} < 1$ dex for Carbon-enriched metal-poor (CEMP) stars, indicating the potential to identify a CEMP stellar population for follow-up confirmation with higher resolution spectroscopy. Finally, we find that large uncertainties in $\alpha$-element abundance measurements using Gaia BP/RP spectra means that efficiently obtaining these abundances will be challenging.

The center of the Milky Way hosts the closest supermassive black hole, SgrA$^*$. Decades of near-infrared observations of our Galactic Center have shown the presence of a small population of stars (the so called S-star cluster) orbiting SgrA$^*$, which were recently reported to be arranged in two orthogonal disks. In this case, the timescale for Lense-Thirring precession of S-stars should be longer than their age, implying a low spin for SgrA$^*$. In contrast, the recent results by the Event Horizon Telescope favor a highly-spinning SgrA$^*$, which seem to suggest that the S-stars could not be arranged in disks. Alternatively, the spin of SgrA$^*$ must be small, suggesting that the models for its observed image are incomplete.

We present our sixth work in a series dedicated to variability studies of active galactic nuclei (AGN) based on the survey of the COSMOS field by the VLT Survey Telescope (VST). Its 54 r-band visits over 3.3 yr and single-visit depth of 24.6 r-band mag make this dataset a valuable scaled-down version that can help forecast the performance of the Rubin Observatory Legacy Survey of Space and Time (LSST). This work is centered on the analysis of the structure function (SF) of VST-COSMOS AGN, investigating possible differences in its shape and slope related to how the AGN were selected, and explores possible connections between the ensemble variability of AGN and black-hole mass, accretion rate, bolometric luminosity, redshift, and obscuration of the source. Given its features, our dataset opens up the exploration of samples ~2 mag fainter than most of the literature to date. We identify several samples of AGN - 677 in total - obtained by a variety of selection techniques which partly overlap. Our analysis compares results for the various samples. We split each sample in two based on the median of the physical property of interest, and analyze differences in the shape and slope of the SF, and possible causes. While the shape of the SF does not change with depth, it is highly affected by the type of AGN (unobscured/obscured) included in the sample. Where a linear region can be identified, we find that the variability amplitude anticorrelates with accretion rate and bolometric luminosity, consistent with previous literature on the topic, while no dependence on black-hole mass emerges from this study. With its longer baseline and denser and more regular sampling, the LSST will allow an improved characterization of the SF and its dependencies on the mentioned physical properties over much larger AGN samples.

We study the origin of the gamma rays from the supernova remnant (SNR) RX J1713.7-3946. Using an analytical model, we calculate the distribution of cosmic rays (CRs) around the SNRs. Motivated by the results of previous studies, we assume that the SNR is interacting with two-phase interstellar medium (ISM), where dense clumps are surrounded by tenuous interclump medium. We also assume that only higher-energy protons ~>TeV) can penetrate the dense clumps. We find that pi^0-decay gamma rays produced by protons reproduce the observed gamma-ray spectrum peaked at ~TeV. On the other hand, it has recently been indicated that the observed ISM column density (N_p), the X-ray surface brightness (I_X), and the gamma-ray surface brightness (I_g) at grid points across the SNR form a plane in the three-dimensional (3D) space of (N_p, I_X, I_g). We find that the planar configuration is naturally reproduced if the ISM or the CR electron-to-proton ratio is not spherically uniform. We show that the shift of the observed data in the 3D space could be used to identify which of the quantities, the ISM density, the CR electron-to-proton ratio, or the magnetic field, varies in the azimuthal direction of the SNR.

In the coming decade, thousands of stellar streams will be observed in the halos of external galaxies. What fundamental discoveries will we make about dark matter from these streams? As a first attempt to look at these questions, we model Magellan/Megacam imaging of the Centaurus A's (Cen A) disrupting dwarf companion Dwarf 3 (Dw3) and its associated stellar stream, to find out what can be learned about the Cen A dark-matter halo. We develop a novel external galaxy stream-fitting technique and generate model stellar streams that reproduce the stream morphology visible in the imaging. We find that there are many viable stream models that fit the data well, with reasonable parameters, provided that Cen A has a halo mass larger than M$_{200}$ $>4.70\times 10^{12}$ M$_{\odot}$. There is a second stream in Cen A's halo that is also reproduced within the context of this same dynamical model. However, stream morphology in the imaging alone does not uniquely determine the mass or mass distribution for the Cen A halo. In particular, the stream models with high likelihood show covariances between the inferred Cen A mass distribution, the inferred Dw3 progenitor mass, the Dw3 velocity, and the Dw3 line-of-sight position. We show that these degeneracies can be broken with radial-velocity measurements along the stream, and that a single radial velocity measurement puts a substantial lower limit on the halo mass. These results suggest that targeted radial-velocity measurements will be critical if we want to learn about dark matter from extragalactic stellar streams.

Machine learning can play a powerful role in inferring missing line-of-sight velocities from astrometry in surveys such as Gaia. In this paper, we apply a neural network to Gaia Early Data Release 3 (EDR3) and obtain line-of-sight velocities and associated uncertainties for ~92 million stars. The network, which takes as input a star's parallax, angular coordinates, and proper motions, is trained and validated on ~6.4 million stars in Gaia with complete phase-space information. The network's uncertainty on its velocity prediction is a key aspect of its design; by properly convolving these uncertainties with the inferred velocities, we obtain accurate stellar kinematic distributions. As a first science application, we use the new network-completed catalog to identify candidate stars that belong to the Milky Way's most recent major merger, Gaia-Sausage-Enceladus (GSE). We present the kinematic, energy, angular momentum, and spatial distributions of the ~450,000 GSE candidates in this sample, and also study the chemical abundances of those with cross matches to GALAH and APOGEE. The network's predictive power will only continue to improve with future Gaia data releases as the training set of stars with complete phase-space information grows. This work provides a first demonstration of how to use machine learning to exploit high-dimensional correlations on data to infer line-of-sight velocities, and offers a template for how to train, validate and apply such a neural network when complete observational data is not available.

We present photometric and spectroscopic observations of the nearby ($D\approx28$ Mpc) interacting supernova (SN) 2019esa, discovered within hours of explosion and serendipitously observed by the Transiting Exoplanet Survey Satellite (TESS). Early, high cadence light curves from both TESS and the DLT40 survey tightly constrain the time of explosion, and show a 30 day rise to maximum light followed by a near constant linear decline in luminosity. Optical spectroscopy over the first 40 days revealed a highly reddened object with narrow Balmer emission lines seen in Type IIn supernovae. The slow rise to maximum in the optical lightcurve combined with the lack of broad H$\alpha$ emission suggest the presence of very optically thick and close circumstellar material (CSM) that quickly decelerated the supernova ejecta. This CSM was likely created from a massive star progenitor with an $\dot{M}$ $\sim$ 0.3 M$_{\odot}$ yr$^{-1}$ lost in a previous eruptive episode 3--4 years before eruption, similar to giant eruptions of luminous blue variable stars. At late times, strong intermediate-width Ca II, Fe I, and Fe II lines are seen in the optical spectra, identical to those seen in the superluminous interacting SN 2006gy. The strong CSM interaction masks the underlying explosion mechanism in SN 2019esa, but the combination of the luminosity, strength of the H$\alpha$ lines, and mass loss rate of the progenitor all point to a core collapse origin.

On large ground-based telescopes, the combination of extreme adaptive optics (ExAO) and coronagraphy with high-dispersion spectroscopy (HDS), sometimes referred to as high-dispersion coronagraphy (HDC), is starting to emerge as a powerful technique for the direct characterisation of giant exoplanets. The high spectral resolution not only brings a major gain in terms of accessible spectral features but also enables a better separation of the stellar and planetary signals. Ongoing projects such as Keck/KPIC, Subaru/REACH, and VLT/HiRISE base their observing strategy on the use of a few science fibres, one of which is dedicated to sampling the planet's signal, while the others sample the residual starlight in the speckle field. The main challenge in this approach is to blindly centre the planet's point spread function (PSF) accurately on the science fibre, with an accuracy of less than 0.1 $\lambda/D$ to maximise the coupling efficiency. In the context of the HiRISE project, three possible centring strategies are foreseen, either based on retro-injecting calibration fibres to localise the position of the science fibre or based on a dedicated centring fibre. We implemented these three approaches, and we compared their centring accuracy using an upgraded setup of the MITHiC high-contrast imaging testbed, which is similar to the setup that will be adopted in HiRISE. Our results demonstrate that reaching a specification accuracy of 0.1 $\lambda/D$ is extremely challenging regardless of the chosen centring strategy. It requires a high level of accuracy at every step of the centring procedure, which can be reached with very stable instruments. We studied the contributors to the centring error in the case of MITHiC and we propose a quantification for some of the most impacting terms.

We study diffuse radio emission in the galaxy cluster A1550, with the aim of constraining particle re-acceleration in the intra-cluster medium. We exploit observations at four different frequencies: 54, 144, 400 and 1400 MHz. To complement our analysis, we make use of archival Chandra X-ray data. At all frequencies we detect an ultra-steep spectrum radio halo ($S_\nu \propto \nu^{-1.6}$) with an extent of 1.2 Mpc at 54 MHz. Its morphology follows the distribution of the thermal intra-cluster medium inferred from the Chandra observation. West of the centrally located head-tail radio galaxy, we detect a radio relic with projected extent of 500 kpc. From the relic, a 600 kpc long bridge departs and connect it to the halo. Between the relic and the radio galaxy, we observe what is most likely a radio phoenix, given its curved spectrum. The phoenix is connected to the tail of the radio galaxy through two arms, which show a nearly constant spectral index for 300 kpc. The halo could be produced by turbulence induced by a major merger, with its axis lying in the NE-SW direction. This is supported by the position of the relic, whose origin could be attributed to a shock propagating along the merger axis. It is possible that the same shock has also produced the phoenix through adiabatic compression, while the bridge could be generated by electrons which were pre-accelerated by the shock, and then re-accelerated by turbulence. Finally, we detect hints of gentle re-energisation in the two arms which depart from the tail of the radio galaxy.

Mounting evidence suggests that Luminous Fast Blue Optical Transients (LFBOTs) are powered by a compact object, launching an asymmetric and fast outflow responsible for the radiation observed in the ultraviolet, optical, infrared, radio, and X-ray bands. Proposed scenarios aiming to explain the electromagnetic emission include an inflated cocoon, surrounding a jet choked in the extended stellar envelope. In alternative, the observed radiation may arise from the disk formed by the delayed merger of a black hole with a Wolf-Rayet star. We explore the neutrino production in these scenarios, i.e. internal shocks in a choked jet and interaction between the outflow and the circumstellar medium (CSM). The choked jet provides the dominant contribution to the neutrino fluence. Intriguingly, the IceCube upper limit on the neutrino emission inferred from the closest LFBOT, AT2018cow, excludes a region of the parameter space otherwise allowed by electromagnetic observations. After correcting for the Eddington bias on the observation of cosmic neutrinos, we conclude that the emission from a choked jet and CSM interaction is compatible with the detection of two track-like neutrino events observed by the IceCube Neutrino Observatory in coincidence with AT2018cow, and otherwise considered to be of atmospheric origin. While the neutrino emission from LFBOTs does not constitute the bulk of the diffuse background of neutrinos observed by IceCube, detection prospects of nearby LFBOTs with IceCube and the upcoming IceCube-Gen2 are encouraging; neutrinos could be observed up to $300$ Mpc and $10^{4}$ Mpc from the CSM interaction and choked jet, respectively. Follow-up neutrino searches will be crucial for unravelling the mechanism powering this emergent transient class.

The successive addition of H atoms to CO in the solid phase has been hitherto regarded as the primary route to form methanol in dark molecular clouds. However, recent Monte Carlo simulations of interstellar ices alternatively suggested the radical-molecule H-atom abstraction reaction CH3O + H2CO -> CH3OH + HCO, in addition to CH3O + H -> CH3OH, as a very promising and possibly dominating (70 - 90 %) final step to form CH3OH in those environments. Here, we compare the contributions of these two steps leading to methanol by experimentally investigating hydrogenation reactions on H2CO and D2CO ices, which ensures comparable starting points between the two scenarios. The experiments are performed under ultrahigh vacuum conditions and astronomically relevant temperatures, with H:H2CO (or D2CO) flux ratios of 10:1 and 30:1. The radical-molecule route in the partially deuterated scenario, CHD2O + D2CO -> CHD2OD + DCO, is significantly hampered by the isotope effect in the D-abstraction process, and can thus be used as an artifice to probe the efficiency of this step. We observe a significantly smaller yield of D2CO + H products in comparison to H2CO + H, implying that the CH3O-induced abstraction route must play an important role in the formation of methanol in interstellar ices. Reflection-Absorption InfraRed Spectroscopy (RAIRS) and Temperature Programmed Desorption-Quadrupole Mass Spectrometry (TPD-QMS) analyses are used to quantify the species in the ice. Both analytical techniques indicate constant contributions of ~80 % for the abstraction route in the 10 - 16 K interval, which agrees well with the Monte Carlo conclusions. Additional H2CO + D experiments confirm these conclusions.

We present a time variability analysis of broad absorption lines (BAL; spread over the velocity range of 5800-29000 km/s) seen in the spectrum of J132216.25+052446.3 (z(em)= 2.04806) at ten different epochs spanning over 19 years. The strongest absorption component (BAL-A; spread over 5800-9900 km/s) is made up of several narrow components having velocity separations close to C IV doublet splitting. The C IV, N V, and Si IV absorption from BAL-A show correlated optical depth variability without significant changes in the velocity structure. A very broad and shallow absorption (BAL-C; spread over the velocity range 15000-29000 km/s) emerged during our monitoring period coinciding with a dimming episode of J1322+0524. All the identified absorption lines show correlated variability with the equivalent widths increasing with decreasing flux. This together with the C IV emission-line variability is consistent with ionization being the main driver of the correlated variability. The observed UV-continuum variations are weaker than what is required by the photo-ionization models. This together with a scatter in the C iv equivalent width at a given continuum flux can be understood if variations of the C IV ionizing photons are much larger than that of the UV continuum, the variations in the ionizing photon and UV fluxes are not correlated and/or the covering factor of the flow varies continuously. We suggest BAL-A is produced by a stable clumpy outflow located beyond the broad emission line region and BAL-C is a newly formed wind component located near the accretion disk and both respond to changes in the ionizing continuum.

The 2017 detection of a kilonova coincident with gravitational-wave emission has identified neutron star mergers as the major source of the heaviest elements, and dramatically constrained alternative theories of gravity. Observing a population of such sources has the potential to transform cosmology, nuclear physics, and astrophysics. However, with only one confident detection currently available, modelling the diversity of signals expected from such a population requires improved theoretical understanding. In particular, models which are quick to evaluate, and are calibrated with more detailed multi-physics simulations, are needed to design observational strategies for kilonovae detection, and to obtain rapid-response interpretations of new observations. We use grey-opacity models to construct populations of kilonovae, spanning ejecta parameters predicted by numerical simulations. Our modelling focuses on wavelengths relevant for upcoming optical surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST). In these simulations, we implement heating rates that are based on nuclear reaction network calculations. We create a Gaussian-process emulator for kilonova grey opacities, calibrated with detailed radiative transfer simulations. Using recent fits to numerical relativity simulations, we predict how the ejecta parameters from BNS mergers shape the population of kilonovae, accounting for the viewing-angle dependence. Our simulated population of binary neutron star (BNS) mergers produce peak i-band absolute magnitudes $-17 \leq M_i \leq -11$. A comparison with detailed radiative transfer calculations indicates that further improvements are needed to accurately reproduce spectral shapes over the full light curve evolution.

Observations of $\rm SgrA^*$ have provided a lot of insight on low-luminosity accretion, with a handful of bright flares accompanied with orbital motion close to the horizon. It has been proposed that gas supply comes from stellar winds in the neighborhood of the supermassive black hole. We here argue that the flow at the vicinity of the black hole has a low magnetization and a structure of alternating polarity totally dictated by the well studied and long-ago proposed MRI turbulent process. This can be the case, provided that in larger distances from the black hole magnetic diffusivity is dominant and thus the magnetic field will never reach equipartition values. For $\rm SgrA^*$, we show the immediate consequences of this specific magnetic field geometry, which are: (i) an intermittent flow that passes from quiescent states to flaring activity, (ii) no quasi-steady-state jet, (iii) no possibility of a magnetically arrested configuration. Moreover a further distinctive feature of this geometry is the intense magnetic reconnection events, occurring as layers of opposite magnetic polarity are accreted, in the vicinity of the black hole. Finally, we argue that the absence of a jet structure in such case will be a smoking gun in 43 \& 86 GHz observations.

Cosmic rays (CRs) play an important role in many astrophysical systems. Acting on plasma scales to galactic environments, CRs are usually modeled as a fluid, using the CR energy density as the evolving quantity. This method comes with the flaw that the corresponding CR evolution equation is not in conservative form as it contains an adiabatic source term that couples CRs to the thermal gas. In the absence of non-adiabatic changes, instead evolving the CR entropy density is a physically equivalent option that avoids this potential numerical inconsistency. In this work, we study both approaches for evolving CRs in the context of magneto-hydrodynamic (MHD) simulations using the massively parallel moving-mesh code AREPO. We investigate the performance of both methods in a sequence of shock-tube tests with various resolutions and shock Mach numbers. We find that the entropy-conserving scheme performs best for the idealized case of purely adiabatic CRs across the shock while both approaches yield similar results at lower resolution. In this setup, both schemes operate well and almost independently of the shock Mach number. Taking active CR acceleration at the shock into account, the energy-based method proves to be numerically much more stable and significantly more accurate in determining the shock velocity, in particular at low resolution, which is more typical for astrophysical large-scale simulations. For a more realistic application, we simulate the formation of several isolated galaxies at different halo masses and find that both numerical methods yield almost identical results with differences far below common astrophysical uncertainties.

Massive Black Hole (MBH) binaries are considered to be one of the most important sources of Gravitational Waves (GW) that can be detected by GW detectors like LISA. However, there are a lot of uncertainties in the dynamics of MBH binaries in the stages leading up to the GW-emission phase. It has been recently suggested that Nuclear Star Clusters (NSCs) could provide a viable route to overcome the final parsec problem for MBH binaries at the center of galaxies. NSCs are collisional systems where the dynamics would be altered by the presence of a mass spectrum. In this study, we use a suite of $N$-body simulations to understand how collisional relaxation under the presence of a mass spectrum of NSC particles affects the dynamics of the MBH binary under the merger of two NSCs. To understand the differences, we consider MBH binaries with different mass ratios and additional non-relaxed models. We find that mass-segregation driven by collisional relaxation can lead to accelerated hardening in lower mass ratio binaries but has the opposite effect in higher mass ratio binaries.The relaxed models demonstrate a lower mass density compared to the non-relaxed models within the influence radius and higher density only very close to the MBH which leads to a less efficient orbital decay in higher mass-ratio models. The results are robust and highlight the importance of subtle differences in changing the dynamics of the binary. Our models are fully collisional and use high enough particle numbers to model NSCs realistically.

We discuss the production of $\gamma$-rays from cosmic rays (CR) in the circumgalactic medium (CGM) of Andromeda (M31) in light of the recent detection of $\gamma$-rays from an annular region of $\sim 5.5-120$ kpc away from the M31 disc. We consider the CRs accelerated as a result of the star-formation in the M31 disk, which are lifted to the CGM by advection due to outflow and CR diffusion. The advection time scale due to bulk flow of gas triggered by star formation activity in the M31 disc is comparable ($\sim$ Gyr) to the diffusion time scale with diffusion coefficient $\ge10^{29}$ cm$^2$ s$^{-1}$ for the propagation of CR protons with energy $\sim 412$ GeV that are responsible for the highest energy photons observed. We show that a leptonic origin of the $\gamma$-rays from cosmic ray (CR) electrons has difficulties, as the inverse Compton time scale ($\sim$Myr) is much lower than advection time scale ($\sim$Gyr) to reach $120$ kpc. Invoking CR electrons accelerated by accretion shocks in the CGM at $\sim100-120$ kpc does not help since it would lead to diffuse X-ray features that are not observed. We, therefore, study the production of $\gamma$-rays via hadronic interaction between CR protons and CGM gas with the help of numerical two-fluid (thermal + CR) hydrodynamical simulation. We find that a combination of these mechanisms, that are related to the star formation processes in M31 in the last $\sim $ Gyr, along with diffusion and hadronic interaction, can explain the observed flux from the CGM of M31.

Aims. We study the dynamical evolution of a system consisting of the giant planets and a massive planetesimal disk over the age of the Solar System. The main question addressed in this study is whether distant trans-Neptunian objects could have come about as a result of the combined action of planetary perturbations and the self-gravity of the disk. Methods. We carried out a series of full N-body numerical simulations of gravitational interactions between the giant planets and a massive outer disk of planetesimals. Results. Our simulations show that the collective gravity of the giant planets and massive planetesimals produces distant trans-Neptunian objects across a wide range of the initial disk mass. The majority of objects that survive up through the age of the Solar System have perihelion distances of q > 40 au. In this region, there is a tendency toward a slow decrease in eccentricities and an increase in perihelion distances for objects with semimajor axes a > 150 au. Secular resonances between distant planetesimals play a major role in increasing their perihelion distances. This explains the origin of Sedna-type objects. In our integrations for the age of the Solar System, we registered times with both high and low clustering of longitudes of perihelion and arguments of perihelion for objects with q > 40 au, a > 150 au. The resulting distribution of inclinations in our model and the observed distribution of inclinations for distant trans-Neptunian objects have similar average values of around 20 degrees. Conclusions. Distant trans-Neptunian objects are a natural consequence in the models that include migrating giant planets and a self-gravitating planetesimal disk.

Population inference of gravitational-wave catalogs is a useful tool to translate observations of black-hole mergers into constraints on compact-binary formation. Different formation channels predict identifiable signatures in the astrophysical distributions of source parameters, such as masses and spins. One example within the scenario of isolated binary evolution is mass-ratio reversal: even assuming efficient core-envelope coupling in massive stars and tidal spin-up of the stellar companion by the first-born black hole, a compact binary in which the second- (first-) born black hole is more (less) massive and (non-) spinning can still form through mass transfer from the initially more to less massive progenitor. Using current LIGO/Virgo observations, we measure the fraction of sources in the underlying population with this mass-spin combination and interpret it as a constraint on the occurrence of mass-ratio reversal in massive binary stars. We modify commonly-used population models by including negligible-spin subpopulations and, most crucially, nonidentical component spin distributions. We do not find evidence for subpopulations of black holes with negligible spins and measure the fraction of massive binary stars undergoing mass-ratio reversal to be consistent with zero and $<32\%$ ($99\%$ confidence). The dimensionless spin peaks around $0.2\unicode{x2013}0.3$ appear robust, however, and are yet to be explained by progenitor formation scenarios.

AGN bubbles could play an important role in accelerating high-energy CRs and galactic feedback. Only in nearby galaxies could we have high enough angular resolution in multi-wavelengths to study the sub-kpc environment of the AGN, where the bubbles are produced and strongly interact with the surrounding ISM. In this paper, we present the latest Chandra observations of the Virgo cluster galaxy NGC 4438, which hosts multi-scale bubbles detected in various bands. The galaxy also has low current star formation activity, so these bubbles are evidently produced by the AGN rather than a starburst. We present spatially resolved spectral analysis of the Chandra data of the $\sim3^{\prime\prime}\times5^{\prime\prime}$ ($\sim200{\rm~pc}\times350\rm~pc$) nuclear bubble of NGC 4438. The power law tail in the X-ray spectra can be most naturally explained as synchrotron emission from high-energy CR leptons. The hot gas temperature increases, while the overall contribution of the non-thermal X-ray emission decreases with the vertical distance from the galactic plane. We calculate the synchrotron cooling timescale of the CR leptons responsible for the non-thermal hard X-ray emission to be only a few tens to a few hundreds of years. The thermal pressure of the hot gas is about three times the magnetic pressure, but the current data cannot rule out the possibility that they are still in pressure balance. The spatially resolved spectroscopy presented in this paper may have important constraints on how the AGN accelerates CRs and drives outflows. We also discover a transient X-ray source only $\sim5^{\prime\prime}$ from the nucleus of NGC 4438. The source was not detected in 2002 and 2008, but became quite X-ray bright in March 2020, with an average 0.5-7 keV luminosity of $\sim10^{39}\rm~ergs~s^{-1}$.

Cosmological simulations predict more classical bulges than their observational counterpart in the local Universe. Here, we quantify evolution of the bulges since $z=0.1$ using photometric parameters of nearly 39,000 unbarred disc galaxies from SDSS DR7 which are well represented by two components. We adopted a combination of the S\'ersic index and Kormendy relation to separate classical bulges and disc-like pseudo-bulges. We found that the fraction of pseudo-bulges (classical bulges) smoothly increases (decreases) as the Universe gets older. In the history of the Universe, there comes a point ($z \approx 0.016$) when classical bulges and pseudo-bulges become equal in number. The fraction of pseudo-bulges rises with increasing bulge to disc half-light radius ratio until R$_{\rm e}$/R$_{\rm hlr} \approx 0.6$ suggesting concentrated disc is the most favourable place for pseudo-bulge formation. The mean ellipticity of pseudo-bulges is always greater than that of classical bulges and it decreases with decreasing redshift indicating that the bulges tend to be more axisymmetric with evolution. Also, the massive bulges are progressing towards axisymmetry at steeper rate than the low-mass bulges. There is no tight correlation of bulge S\'ersic index evolution with other photometric properties of the galaxy. Using the sample of multi-component fitting of $S^4G$ data and $N-$body galaxy models, we have verified that our results are consistent or even more pronounced with multi-component fitting and high-resolution photometry.

Fast Radio Bursts must be powered by uniquely energetic emission mechanisms. This requirement has eliminated a number of possible source types, but several remain. Identifying the physical nature of Fast Radio Burst (FRB) emitters arguably requires good localisation of more detections, and broadband studies enabled by real-time alerting. We here present the Apertif Radio Transient System (ARTS), a supercomputing radio-telescope instrument that performs real-time FRB detection and localisation on the Westerbork Synthesis Radio Telescope (WSRT) interferometer. It reaches coherent-addition sensitivity over the entire field of the view of the primary dish beam. After commissioning results verified the system performed as planned, we initiated the Apertif FRB survey (ALERT). Over the first 5 weeks we observed at design sensitivity in 2019, we detected 5 new FRBs, and interferometrically localised each of these to 0.4--10 sq. arcmin. All detections are broad band and very narrow, of order 1 ms duration, and unscattered. Dispersion measures are generally high. Only through the very high time and frequency resolution of ARTS are these hard-to-find FRBs detected, producing an unbiased view of the intrinsic population properties. Most localisation regions are small enough to rule out the presence of associated persistent radio sources. Three FRBs cut through the halos of M31 and M33. We demonstrate that Apertif can localise one-off FRBs with an accuracy that maps magneto-ionic material along well-defined lines of sight. The rate of 1 every ~7 days next ensures a considerable number of new sources are detected for such study. The combination of detection rate and localisation accuracy exemplified by the 5 first ARTS FRBs thus marks a new phase in which a growing number of bursts can be used to probe our Universe.

Thermonuclear supernovae are the result of the violent unbinding of a white dwarf, but the precise nature of the explosion mechanism(s) is a matter of active debate. To this end, several specific scenarios have been proposed to explain the observable traits of SNe Ia. A promising pathway is the double-detonation scenario, where a white dwarf accretes a shell of helium-rich material from a companion and a detonation in the resulting helium shell is the primary cause of the explosion. Through a set of two-dimensional grid-based simulations of this scenario we clearly distinguish three phases of evolution: external helium-rich detonation, core compressive heating, and a final core carbon burn. Though final disruption of the whole system is achieved at all resolutions, only models with minimum resolutions of 4~km and better exhibit all three phases. Particularly, core compression heating is only observed for higher resolutions, producing qualitatively different nucleosynthetic outcomes. We identify the effect of finer spatial resolution on the mixing of hot silicon at the interface between the detonating helium layer and the underlying C/O WD as a primary driver of these dynamic differences.

Spectral data cubes of the interacting pair of galaxies NGC 2535 and NGC 2536 (the Arp 82 system) targeting bright emission lines in the visible band, obtained with the imaging Fourier transform spectrometer (iFTS) SITELLE attached to the Canada-France-Hawaii Telescope (CFHT), are presented. Analysis of H$\upalpha$ velocity maps reveals a bar in $\rm NGC\,2536$. In $\rm NGC\,2535$, we find strong non-circular motions outside the ocular ring, in the elliptical arc and tidal tails of $\rm NGC\,2535$ and a misalignment between the kinematic and photometric position angles. We detect 155 HII region complexes in the interacting pair of galaxies and determine oxygen abundances for 66 of them using different calibrators. We find, regardless of the indicator used, that the oxygen abundance distribution in $\rm NGC\,2536$ is shallow whereas, in $\rm NGC\,2535$, it is best fitted by two slopes, the break occurring beyond the ocular ring. The inner slope is comparable to the one observed in isolated normal star-forming galaxies but the outer slope is shallow. We present a numerical simulation of the interaction that reproduces the observed tidal features, kinematics, and metallicity distribution, to investigate the effect of the interaction on the galaxies. The model indicates that the galaxies have undergone a close encounter, strongly prograde for the primary, and are half way in their course to a second close encounter.

The possible encounter of the meteoric material from 73P/Schmassmann--Wachmann~3 produced during the comet's 1995 outburst provides a rare and valuable opportunity to understand a fragmenting comet. Here we explore various ejection configurations and their impact on the possible meteor outburst in the early hours of UT 2022 May 31. Based on available evidence, it is likely that only sub-millimeter-class meteoroids will be delivered to the Earth, producing an outburst rich in extremely faint meteors difficult to be detected by conventional techniques. However, we confirm that the meteoroids only need to be ejected as little as 40\% faster in order for larger, millimeter-class meteoroids to reach the Earth and generate a more visible, albeit likely still small, outburst. Other effects such as an enhanced lunar sodium tail and a visible glow from the meteoroid trail may also occur during the encounter.

The proper characterization of the general statistical behavior of these fluctuations, from a limited sample of observations or simulations, is of prime importance to understand the process of star formation. In this article, we use the ergodic theory for any random field of fluctuations, as commonly used in statistical physics, to derive rigorous statistical results. We outline how to evaluate the autocovariance function (ACF) and the characteristic correlation length of these fluctuations. We then apply this statistical approach to astrophysical systems characterized by a field of density fluctuations, notably star-forming clouds. When it is difficult to determine the correlation length from the empirical ACF, we show alternative ways to estimate the correlation length. We show that the statistics of the column-density field is hampered by biases introduced by integration effects along the line of sight and we explain how to reduce these biases. The statistics of the probability density function (PDF) ergodic estimator also yields the derivation of the proper statistical error bars. We provide a method that can be used by observers and numerical simulation specialists to determine the latter. We show that they (i) cannot be derived from simple Poisson statistics and (ii) become increasingly large for increasing density contrasts, severely hampering the accuracy of the low and high end part of the PDF because of a sample size that is too small. As templates of various stages of star formation in MCs, we then examine the case of the Polaris and Orion B clouds in detail. We calculate, from the observations, the ACF and the correlation length in these clouds and show that the latter is on the order of $\sim$1\% of the size of the cloud.

In the first article of this series, we have used the ergodic theory to assess the validity of a statistical approach to characterize various properties of star-forming molecular clouds (MCs) from a limited number of observations or simulations. This allows the proper determination of confidence intervals for various volumetric averages of statistical quantities obtained form observations or numerical simulations. In this joint paper, we apply the same formalism to a different kind of (observational or numerical) study of MCs. Indeed, as observations cannot fully unravel the complexity of the inner density structure of star forming clouds, it is important to know whether global observable estimates, such as the total mass and size of the cloud, can give an accurate estimation of various key physical quantities that characterize the dynamics of the cloud. Of prime importance is the correct determination of the total gravitational (binding) energy and virial parameter of a cloud. We show that, whereas for clouds that are not in a too advanced stage of star formation, such as Polaris or Orion B, the knowledge of only their mass and size is sufficient to yield an accurate determination of the aforementioned quantities from observations (i.e. in real space). In contrast, we show that this is no longer true for numerical simulations in a periodic box. We derive a relationship for the ratio of the virial parameter in these two respective cases.

Follow-up observations of high-magnification gravitational microlensing events can fully exploit their intrinsic sensitivity to detect extrasolar planets, especially those with small mass ratios. To make followup more uniform and efficient, we develop a system, HighMagFinder, based on the real-time data from the Korean Microlensing Telescope Network (KMTNet) to automatically alert possible ongoing high-magnification events. We started a new phase of follow-up observations with the help of HighMagFinder in 2021. Here we report the discovery of two planets in high-magnification microlensing events, KMT-2021-BLG-0171 and KMT-2021-BLG-1689, which were identified by the HighMagFinder. We find that both events suffer the ``central-resonant'' caustic degeneracy. The planet-host mass-ratio is $q\sim4.7\times10^{-5}$ or $q\sim 2.2\times10^{-5}$ for KMT-2021-BLG-0171, and $q\sim2.5\times10^{-4}$ or $q\sim 1.8\times10^{-4}$ for KMT-2021-BLG-1689. Together with two events reported by Ryu et al. (2022), four cases that suffer such degeneracy have been discovered in the 2021 season alone, indicating that the degenerate solutions may have been missed in some previous studies. We also propose a new factor for weighting the probability of each solution from the phase-space. The resonant interpretations for the two events are disfavored under this consideration. This factor can be included in future statistical studies to weight degenerate solutions.

The next generation of cosmological observations will be sensitive to small deviations from a pure power law in the primordial power spectrum of the curvature perturbations. In the context of slow-roll inflation, these deviations are expected and correspond to the so-called running of the spectral index. Their measurement would bring as much information as the discovery of deviations from scale invariance. However, robust parameter inference requires to marginalize over any possible higher order uncertainties, which have been, up to now, not fully determined. We tackle this issue by deriving the inflationary scalar and tensor slow-roll power spectra at next-to-next-to-next to leading order (N3LO), fully expanded around an observable pivot wavenumber, for all single field inflationary models having minimal and non-minimal kinetic terms. Our result therefore encompasses string-inspired inflationary models having a varying speed of sound.

We studied the cloud-cloud collision candidate G323.18+0.15 based on signatures of induced filaments, clumps, and star formation. We used archival molecular spectrum line data from the SEDIGISM $^{13}$CO($J$\,=\,2--1) survey, from the Mopra southern Galactic plane CO survey, and infrared to radio data from the GLIMPSE, MIPS, Hi-GAL, and SGPS surveys. Our new result shows that the G323.18+0.15 complex is 3.55kpc away from us and consists of three cloud components, G323.18a, G323.18b, and G323.18c. G323.18b shows a perfect U-shape structure, which can be fully complemented by G323.18a, suggesting a collision between G323.18a and the combined G323.18bc filamentary structure. One dense compressed layer (filament) is formed at the bottom of G323.18b, where we detect a greatly increased velocity dispersion. The bridge with an intermediate velocity in a position-velocity diagram appears between G323.18a and G323.18b, which corresponds to the compressed layer. G323.18a plus G323.18b as a whole are probably not gravitationally bound. This indicates that high-mass star formation in the compressed layer may have been caused by an accidental event. The column density in the compressed layer of about $1.36 \times 10^{22}$cm$^{-2}$ and most of the dense clumps and high-mass stars are located there. The average surface density of classI and classII young stellar objects (YSOs) inside the G323.18+0.15 complex is much higher than the density in the surroundings. The timescale of the collision between G323.18a and G323.18b is $1.59$Myr. This is longer than the typical lifetime of classI YSOs and is comparable to the lifetime of classII YSOs.

The cosmological principle asserts that the Universe looks spatially homogeneous and isotropic on sufficiently large scales. Given the fundamental implications of the cosmological principle, it is important to empirically test its validity on various scales. In this paper, we use the Type Ia supernova (SN~Ia) magnitude-redshift relation, from both the Pantheon and JLA compilations, to constrain theoretically motivated anisotropies in the Hubble flow. In particular, we constrain the quadrupole moment in the effective Hubble parameter and the dipole moment in the effective deceleration parameter. We find no significant quadrupole term regardless of the redshift frame we use. Our results are consistent with the theoretical expectation of a quadrupole moment of a few percent at scales of $\sim 100 h^{-1}$ Mpc. We place an upper limit of a $\sim 10\%$ quadrupole amplitude relative to the monopole, $H_0$, at these scales. We find that we can detect a $\sim 7\%$ quadrupole moment at the 5$\sigma$ level, for a forecast low-$z$ sample of 1055 SNe~Ia. We find an exponentially decaying dipole moment of the deceleration parameter varies in significance depending on the redshift frame we use. In the heliocentric frame, as expected, it is detected at $\sim 3 \sigma$ significance. In the rest-frame of the cosmic microwave background (CMB), we find a marginal $\sim 2 \sigma$ dipole, however, after applying peculiar velocity corrections, the dipole is insignificant. Finally, we find the best-fit frame of rest relative to the supernovae to differ from that of the CMB.

Analyses of quasi-periodic oscillations (QPOs) are important to understanding the dynamic behaviour in many astrophysical objects during transient events like gamma-ray bursts, solar flares, magnetar flares and fast radio bursts. Astrophysicists often search for QPOs with frequency-domain methods such as (Lomb-Scargle) periodograms, which generally assume power-law models plus some excess around the QPO frequency. Time-series data can alternatively be investigated directly in the time domain using Gaussian Process (GP) regression. While GP regression is computationally expensive in the general case, the properties of astrophysical data and models allow fast likelihood strategies. Heteroscedasticity and non-stationarity in data have been shown to cause bias in periodogram-based analyses. Gaussian processes can take account of these properties. Using GPs, we model QPOs as a stochastic process on top of a deterministic flare shape. Using Bayesian inference, we demonstrate how to infer GP hyperparameters and assign them physical meaning, such as the QPO frequency. We also perform model selection between QPOs and alternative models such as red noise and show that this can be used to reliably find QPOs. This method is easily applicable to a variety of different astrophysical data sets. We demonstrate the use of this method on a range of short transients: a gamma-ray burst, a magnetar flare, a magnetar giant flare, and simulated solar flare data.

GRB 201015A is a long-duration Gamma-Ray Burst (GRB) which was detected at very high energies (> 100 GeV) using the MAGIC telescopes. If confirmed, this would be the fifth and least luminous GRB ever detected at this energies. We performed a radio follow-up of GRB 201015A over twelve different epochs, from 1.4 to 117 days post-burst, with the Karl G. Jansky Very Large Array, e-MERLIN and the European VLBI Network. We included optical and X-rays observations, performed with the Multiple Mirror Telescope and the Chandra X-ray Observatory respectively, together with publicly available data. We detected a point-like transient, consistent with the position of GRB 201015A until 23 and 47 days post-burst at 1.5 and 5 GHz, respectively. The source was detected also in both optical (1.4 and 2.2 days post-burst) and X-ray (8.4 and 13.6 days post-burst) observations. The multi-wavelength afterglow light curves can be explained with the standard model for a GRB seen on-axis, which expands and decelerates into a medium with a homogeneous density, while a circumburst medium with a wind-like profile is disfavoured. Notwithstanding the high resolution provided by the VLBI, we could not pinpoint any expansion or centroid displacement of the outflow. If the GRB is seen at the viewing angle which maximises the apparent velocity, we estimate that the Lorentz factor for the possible proper motion is $\Gamma_{\alpha}$ < 40 in right ascension and $\Gamma_{\delta}$ < 61 in declination. On the other hand, if the GRB is seen on-axis, the size of the afterglow is <5 pc and <16 pc at 25 and 47 days. Finally, the early peak in the optical light curve suggests the presence of a reverse shock component before 0.01 days from the burst.

We use the spectral lag catalog of 46 short GRBs obtained by Xiao et al (2022), between two fixed energy intervals in the source frame, to carry out an independent search for Lorentz Invariance violation (LIV). For this purpose, we use a power-law model as a function of energy for the intrinsic astrophysical induced spectral lags. The expansion history of the universe needed for the evaluation of LIV was obtained in a non-parametric method using cosmic chronometers. We use Bayesian model comparison to determine if the aforementioned spectral lags show evidence for LIV as compared to only astrophysically induced lags. We do any find any evidence for LIV, and obtain 95\% c.l. lower limits for the energy scale of LIV to be $4 \times 10^{15}$ GeV and $6.8 \times 10^{9}$ GeV for the linear and quadratic LIV models respectively.

Aiming at the detection of cosmological gas being accreted onto galaxies of the local Universe, we examined the Halpha emission in the halo of 164 galaxies in the field of view of the Multi-Unit Spectroscopic Explorer Wide survey (\musew ) with observable Halpha (redshift < 0.42). An exhaustive screening of the corresponding Halpha images led us to select 118 reliable Halpha emitting gas clouds. The signals are faint, with a surface brightness of 10**(-17.3 pm 0.3) erg/s/cm2/arcsec2. Through statistical tests and other arguments, we ruled out that they are created by instrumental artifacts, telluric line residuals, or high redshift interlopers. Around 38% of the time, the Halpha line profile shows a double peak with the drop in intensity at the rest-frame of the central galaxy, and with a typical peak-to-peak separation of the order of pm 200 km/s. Most line emission clumps are spatially unresolved. The mass of emitting gas is estimated to be between one and 10**(-3) times the stellar mass of the central galaxy. The signals are not isotropically distributed; their azimuth tends to be aligned with the major axis of the corresponding galaxy. The distances to the central galaxies are not random either. The counts drop at a distance > 50 galaxy radii, which roughly corresponds to the virial radius of the central galaxy. We explore several physical scenarios to explain this Halpha emission, among which accretion disks around rogue intermediate mass black holes fit the observations best.

Recent studies have proposed that GRAVITY+ instrument is able to trace the circular orbit of the subparsec close-binary supermassive black holes (CB-SMBHs) by measuring the photocentre variation of the hot dust emission. However, the CB-SMBHs orbit may become highly eccentric throughout the evolution of these objects, and the orbital period may be far longer than the observational time baseline. We investigate the problem of detecting the CB-SMBH with hot dust emission and high eccentricity (eCBSMBH, e=0.5) when the observed time baselines of their astrometric data and radial velocities are considerably shorter than the orbital period. The parameter space of the Keplerian model of the eCBSMBH is large for exploratory purposes. We therefore applied the Bayesian method to fit orbital elements of the eCBSMBH to combined radial velocity and astrometric data covering a small fraction of the orbital period. We estimate that a number of potential eCBSMBH systems within reach of GRAVITY + will be similar to the number of the planned circular targets. We show that using observational time baselines that cover ~ 10% of the orbit increases the possibility of determining the period, eccentricity, and total mass of an eCBSMBH. When the observational time baseline becomes too short (~ 5%), the quality of the retrieved eCBSMBH parameters degrades. We also illustrate how interferometry may be used to estimate the photo-centre at the eCBSMBH emission line, which could be relevant for GRAVITY+ successors. Even if the astrometric signal for eCBSMBH systems is reduced by a factor of sqrt{1-e^{2}} compared to circular ones, we find that the hot dust emission of eCBSMBHs can be traced by GRAVITY+ at the elementary level.

Surface-bound exospheres facilitate volatile migration across the surfaces of nearly airless bodies. However, such transport requires that the body can both form and retain an exosphere. To form a sublimation exosphere requires the surface of a body to be sufficiently warm for surface volatiles to sublime; to retain an exosphere, the ballistic escape and photodestruction rates and other loss mechanisms must be sufficiently low. Here we construct a simple free molecular model of exospheres formed by volatile desorption/sublimation. We consider the conditions for forming and retaining exospheres for common volatile species across the Solar System, and explore how three processes (desorption/sublimation, ballistic loss, and photodestruction) shape exospheric dynamics on airless bodies. Our model finds that the CO2 exosphere of Callisto is too dense to be sustained by impact-delivered volatiles, but could be maintained by only ~7 hectares of exposed CO2 ice. We predict the peak surface locations of Callisto's CO2 exosphere along with other Galilean moons, which could be tested by JUICE observations. Our model finds that to maintain Iapetus' two-tone appearance, its dark Cassini Regio likely has unresolved exposures of water ice, perhaps in sub-resolution impact craters, that amount to up to ~0.06% of its surface. In the Uranian system, we find that the CO2 deposits on Ariel, Umbriel, Titania, and Oberon are unlikely to have been delivered via impacts, but are consistent with both a magnetospheric origin or sourced endogenously. We suggest that exosphere-mediated volatile transport could produce these moons' leading/trailing CO2 asymmetries, and may be a seasonal equinox feature that could be largely erased by volatile migration during the Uranian solstices. We calculate that ~2.4-6.4 mm thick layer of CO2 could migrate about the surface of Uranus' large moons during a seasonal cycle.

Bayesian workflows often require the introduction of nuisance parameters, yet for core science modelling one needs access to a marginal posterior density. In this work we use masked autoregressive flows and kernel density estimators to encapsulate the marginal posterior, allowing us to compute marginal Kullback-Leibler divergences and marginal Bayesian model dimensionalities in addition to generating samples and computing marginal log probabilities. We demonstrate this in application to topical cosmological examples of the Dark Energy Survey, and global 21cm signal experiments. In addition to the computation of marginal Bayesian statistics, this work is important for further applications in Bayesian experimental design, complex prior modelling and likelihood emulation. This technique is made publicly available in the pip-installable code margarine.

The radio/X-ray correlation diagram, for Black Hole X-ray Binaries (BHXBs) in the hard state, depicts the connection that might exists between the radio jets and X-ray emitting accretion discs. The current version of the radio/X-ray correlation diagram shows two populations of BHXBs. These two populations evolve along two different correlation tracks, namely the standard track and the outliers. Over the past years, the key question has been to explain the existence of these two tracks. In this paper, we investigate the impact of the black hole mass on the radio/X-ray correlation for a sample of 17 BHXBs. We are led to conclude at least one of the following with full consideration of large uncertainties in mass estimates:(i) Most of the reported mass estimates of black holes are incorrect or insufficiently accurate to infer their impact on radio/X-ray correlation diagram. (ii) The estimated radio luminosities and X-ray luminosites are still not reliable enough for two reasons. One reason is due to the lack of associated errors in observational data. Another is that some sources might still transit from one track to the other. (iii) The mass of BH has a significant influence on which track, the source belongs to, on a radio/X-ray correlation diagram.

We present and apply a method to identify the stellar content of the ROSAT all-sky survey (RASS). We performed a crossmatch between the RASS sources and stellar candidates selected from Gaia Early Data Release 3 (EDR3) and estimated stellar probabilities for every RASS source from the geometric properties of the match and additional properties, namely the X-ray to G-band flux ratio and the counterpart distances. A comparison with preliminary detections from the first eROSITA all-sky survey (eRASS1) show that the positional offsets of the RASS sources are larger than expected from the uncertainties given in the RASS catalog. From the RASS sources with reliable positional uncertainties, we identify 28630 (24.9 %) sources as stellar; this is the largest sample of stellar X-ray sources to date. Directly from the stellar probabilities, we estimate the completeness and reliability of the sample to be about 93 % and confirm this value by comparing it to the identification of randomly shifted RASS sources, preliminary stellar eRASS1 identifications, and results from a previous identification of RASS sources. Our stellar RASS sources contain sources of all spectral types and luminosity classes. According to their position in the color-magnitude diagram, many stellar RASS sources are young stars with ages of a few $10^7$ yr or binaries. When plotting the X-ray to bolometric flux ratio as a function of the color, the onset of convection and the saturation limit are clearly visible. We note that later-type stars reach continuously higher $F_X/F_{bol}$ values, which is probably due to more frequent flaring. The color distribution of the stellar RASS sources clearly differs from the unrelated background sources. We present the three-dimensional distribution of the stellar RASS sources that shows a clear increase in the source density near known stellar clusters.

We predicted, observed, and analyzed the multichord stellar occultation of the Second Gaia Data Release (Gaia DR2) source 3449076721168026624 (m$_v$ = 14.1 mag) by the plutino object 2003 VS$_2$ (hereafter, VS$_2$) on 2019 October 22. We also carried out photometric observations to derive the rotational light curve amplitude and rotational phase of VS$_2$ during the stellar occultation. Combining the results and assuming a triaxial shape, we derived the 3D shape of VS$_2$. Out of the 39 observatories involved in the observational campaign, 12 sites reported a positive detection; this makes it one of the best observed stellar occultations by a TNO so far. We obtained a rotational light curve amplitude of ${\Delta}$m = 0.264 $\pm$ 0.017 mag, a mean area-equivalent diameter of D$_{A_{eq}}$ = 545 $\pm$ 13 km, and a geometric albedo of 0.134 $\pm$ 0.010. The best triaxial shape obtained for VS$_2$ has semiaxes a = 339 $\pm$ 5 km, b = 235 $\pm$ 6 km, and c = 226 $\pm$ 8 km. The derived aspect angle is ${\theta}$ = 59${\deg} \pm$ 2${\deg}$ or its supplementary ${\theta}$ = 121${\deg} \pm$ 2${\deg}$, depending on the north-pole position. The spherical-volume equivalent diameter is D$_{V_{eq}}$ = 524 $\pm$ 7 km. If we consider large albedo patches on its surface, the semi-major axis of the ellipsoid could be ~10 km smaller. These results are compatible with the previous ones determined from the single-chord 2013 and four-chord 2014 stellar occultations and with the effective diameter and albedo derived from Herschel and Spitzer data. They provide evidence that VS$_2$'s 3D shape is not compatible with a homogeneous triaxial body in hydrostatic equilibrium, but it might be a differentiated body and/or might be sustaining some stress. No secondary features related to rings or material orbiting around VS$_2$ were detected.

Astronomers are increasingly faced with a deluge of information, and finding worthwhile targets of study in the sea of data can be difficult. Outlier identification studies are a method that can be used to focus investigations by presenting a smaller set of sources that could prove interesting because they do not follow the trends of the underlying population. We apply a principal component analysis (PCA) and an unsupervised random forest algorithm (uRF) to sources from the Chandra Source Catalog v.2 (CSC2). We present 119 high-significance sources that appear in all repeated applications of our outlier identification algorithm (OIA). We analyze the characteristics of our outlier sources and crossmatch them with the SIMBAD database. Our outliers contain several sources that were previously identified as unusual. This OIA leads to the identification of interesting targets that could motivate more detailed study.

Decaying or annihilating dark matter particles could be detected through gamma-ray emission from the species they decay or annihilate into. This is usually done by modelling the flux from specific dark matter-rich objects such as the Milky Way halo, Local Group dwarfs and nearby groups. However, these objects are expected to have significant emission from baryonic processes as well, and the analyses discard gamma-ray data over most of the sky. Here we construct full-sky templates for gamma-ray flux from the large-scale structure within $\sim$200 Mpc by means of a suite of constrained $N$-body simulations (CSiBORG) produced using the Bayesian Origin Reconstruction from Galaxies algorithm. Marginalising over uncertainties in this reconstruction, small-scale structure and parameters describing astrophysical contributions to the observed gamma ray sky, we compare to observations from the Fermi Large Area Telescope to constrain dark matter annihilation cross-sections and decay rates through a Markov Chain Monte Carlo analysis. We rule out the thermal relic cross-section for $s$-wave annihilation for all $m_\chi \lesssim 7 {\rm \, GeV}/c^2$ at 95% confidence if the annihilation produces $Z$ bosons, gluons or quarks less massive than the bottom quark. We infer a contribution to the gamma ray sky with the same spatial distribution as dark matter decay at $3.3\sigma$. Although this could be due to dark matter decay via these channels with a decay rate $\Gamma \approx 3 \times 10^{-28} {\rm \, s^{-1}}$, we find that a power-law spectrum of index $p=-2.75^{+0.71}_{-0.46}$, likely of baryonic origin, is preferred by the data.

M dwarfs are favorable targets for exoplanet detection with current instrumentation, but stellar companions can induce false positives and inhibit planet characterization. Knowledge of stellar companions is also critical to our understanding of how binary stars form and evolve. We have therefore conducted a survey of stellar companions around nearby M dwarfs, and here we present our new discoveries. Using the DSSI speckle imager at the 4.3-meter Lowell Discovery Telescope, and the similar NESSI instrument at the 3.5-meter WIYN telescope, we carried out a volume-limited survey of M-dwarf multiplicity to 15 parsecs, with a special emphasis on including the later M dwarfs that were overlooked in previous surveys. Additional brighter targets at larger distances were included for a total sample size of 1070 M dwarfs. Observations of these 1070 targets revealed 26 new companions; 22 of these systems were previously thought to be single. If all new discoveries are confirmed, then the number of known multiples in the sample will increase by 7.6%. Using our observed properties, as well as the parallaxes and 2MASS K magnitudes for these objects, we calculate the projected separation, and estimate the mass ratio and component spectral types, for these systems. We report the discovery of a new M-dwarf companion to the white dwarf Wolf 672 A, which hosts a known M-dwarf companion as well, making the system trinary. We also examine the possibility that the new companion to 2MASS J13092185-2330350 is a brown dwarf. Finally, we discuss initial insights from the POKEMON survey.

The thermal structure of planetary atmospheres is an essential input for predicting and retrieving the distribution of gases and aerosols, as well as the bulk chemical abundances. In the case of Jupiter, the temperature at a reference level {\hyphen} generally taken at 1 bar {\hyphen} serves as the anchor in models used to derive the planet{\quotesingle}s interior structure and composition. Most models assume the temperature measured by the Galileo probe (Seiff et al. 1998). However, those data correspond to a single location, an unusually clear, dry region, affected by local atmospheric dynamics. On the other hand, the Voyager radio occultation observations cover a wider range of latitudes, longitudes, and times (Lindal et al. 1981). The Voyager retrievals were based on atmospheric composition and radio refractivity data that require updating and were never properly tabulated{\colon} the few existing tabulations are incomplete and ambiguous. Here, we present a systematic electronic digitization of all available temperature profiles from Voyager, followed by their reanalysis, employing currently accepted values of the abundances and radio refractivities of atmospheric species. We find the corrected temperature at the 1 bar level to be up to 4 K greater than previously published values, i.e., 170.3{\pm}3.8 K at 12{\deg}S (Voyager 1 ingress) and 167.3{\pm}3.8 K at 0{\deg}N (Voyager 1 egress). This is to be compared with the Galileo probe value of 166.1{\pm}0.8 K at the edge of an unusual feature at 6.57{\deg}N. Altogether, this suggests that Jupiter{\quotesingle}s tropospheric temperatures may vary spatially by up to 7 K between 7{\deg}N and 12{\deg}S.

We report the results of analysis of the Hbeta emission-line region of a sample of thirty low-redshift (z<1) iron low-ionization broad absorption line quasars (FeLoBALQs). Eleven of these objects are newly classified as FeLoBALQs. A matched sample of 132 unabsorbed quasars was analyzed in parallel. The emission lines showed the well known anticorrelation between the [OIII] and FeII emission (Boroson & Green 1992). Using a summary statistic called E1 to quantify this anticorrelation, we found that while the distribution of E1 for the unabsorbed quasars has a single peak, the FeLoBALQs have a bimodal shape in this parameter. Previous studies have shown that the line emission properties of BAL and non-BAL quasars are consistent, and therefore the difference in the Hbeta region emission between FeLoBAL quasars and unabsorbed quasars is a new result. The two populations of FeLoBAL quasars are characterized by low and high bolometric luminosities and Eddington ratios. Some previous studies have suggested that BAL quasars are high accretion-rate objects, and therefore the discovery of the low accretion-rate branch of FeLoBAL quasars was unexpected. We also found that the Hbeta FWHM is systematically broader among the FeLoBALQs compared with the non-BAL quasars implying a higher inclination viewing angle or a dearth of low-velocity line-emitting gas.

We present results from simulating slitless spectroscopic observations with the Nancy Grace Roman Space Telescope's (Roman) Wide-Field Instrument (WFI) P127 prism spanning 0.75 $\mu m$ to 1.8 $\mu m$. We quantify the efficiency of recovered Type Ia supernovae (SNe Ia) redshifts, as a function of P127 prism exposure time, to guide planning for future observing programs with the Roman prism. Generating the two-dimensional dispersed images and extracting one-dimensional spectra is done with the slitless spectroscopy package pyLINEAR along with custom-written software. From the analysis of 1698 simulated SN Ia P127 prism spectra, we show the efficiency of recovering SN redshifts to $z\lesssim3.0$, highlighting the exceptional sensitivity of the Roman P127 prism. Redshift recovery is assessed by setting a requirement of $\sigma_z = (\left|z - z_\mathrm{true} \right|)/(1+z) \leq 0.01$. We find that 3 hr exposures are sufficient for meeting this requirement, for $\gtrsim 50\%$ of the sample of mock SNe Ia at $z\approx2$ and within $\pm5$ days of rest-frame maximum light in the optical. We also show that a 1 hr integration of Roman can achieve the same precision in completeness to a depth of $24.4 \pm 0.06$ AB mag (or $z\lesssim 1$). Implications for cosmological studies with Roman P127 prism spectra of SNe Ia are also discussed.

In several approaches to evading the information paradox, the semiclassical black hole is replaced by an Exotic Compact Object (ECO). It has been conjectured that gravitational waves emitted by the merger of ECOs can reflect off the ECOs and produce a detectable `echo'. We argue that while a part of the wave can indeed reflect off the surface of an ECO, this reflected wave will get trapped by a new closed trapped surface produced by its own backreaction. Thus no detectable signal of the echo will emerge to infinity. The only assumption in this analysis is that causality is maintained to leading order in gently curved spacetime. Thus if echoes are actually detected, then we would face a profound change in our understanding of physics.

The current cosmic time evolution of the Universe is described by the General Relativity theory when a cosmological principle is considered under a flat space time landscape. The set of known as Friedmann equations, contain the principles that lead to the construction of the standard $\Lambda$CDM model. However, the current state-of-art regarding these equations, even if it is a fundamental method, is based in solving analytically the differential equations by considering several forms of matter/energy components or evaluating them in specific cosmic times where two or more components contribute at the same rate. This latter can be carry out through the approach of piecewice solutions, whose reduce the numerical integrals. In this paper we discuss new solutions through special analytical functions and constraint them with an updated compilation of observational Hubble observations in order to deal with the local $H_0$ tension reported.

We explore the possibility of gravitationally generated particle production in the scalar-tensor representation of $f(R,T)$ gravity. Due to the explicit nonminimal curvature-matter coupling in the theory, the divergence of the matter energy-momentum tensor does not vanish. We explore the physical and cosmological implications of this property by using the formalism of irreversible thermodynamics of open systems in the presence of matter creation/annihilation. The particle creation rates, pressure, temperature evolution and the expression of the comoving entropy are obtained in a covariant formulation and discussed in detail. Applied together with the gravitational field equations, the thermodynamics of open systems lead to a generalization of the standard $\Lambda$CDM cosmological paradigm, in which the particle creation rates and pressures are effectively considered as components of the cosmological fluid energy-momentum tensor. We also consider specific models, and compare the scalar-tensor $f(R,T)$ cosmology with the $\Lambda$CDM scenario, and the properties of particle creation rates, creation. pressures, and entropy generation through gravitational matter production in both low and high redshift limits.

Geomagnetically Induced Currents (GICs) arise from spatio-temporal changes to Earth's magnetic field which arise from the interaction of the solar wind with Earth's magnetosphere, and drive catastrophic destruction to our technologically dependent society. Hence, computational models to forecast GICs globally with large forecast horizon, high spatial resolution and temporal cadence are of increasing importance to perform prompt necessary mitigation. Since GIC data is proprietary, the time variability of horizontal component of the magnetic field perturbation (dB/dt) is used as a proxy for GICs. In this work, we develop a fast, global dB/dt forecasting model, which forecasts 30 minutes into the future using only solar wind measurements as input. The model summarizes 2 hours of solar wind measurement using a Gated Recurrent Unit, and generates forecasts of coefficients which are folded with a spherical harmonic basis to enable global forecasts. When deployed, our model produces results in under a second, and generates global forecasts for horizontal magnetic perturbation components at 1-minute cadence. We evaluate our model across models in literature for two specific storms of 5 August 2011 and 17 March 2015, while having a self-consistent benchmark model set. Our model outperforms, or has consistent performance with state-of-the-practice high time cadence local and low time cadence global models, while also outperforming/having comparable performance with the benchmark models. Such quick inferences at high temporal cadence and arbitrary spatial resolutions may ultimately enable accurate forewarning of dB/dt for any place on Earth, resulting in precautionary measures to be taken in an informed manner.

Cooling processes of brown dwarf stars and giant planets are studied in the framework of DHOST theories. We confirm the previous results in the field that the effect of modified gravity on substellar objects' age is pronounced the most.

Positivity bounds on scattering amplitudes provide a necessary condition for a low-energy effective field theory to have a consistent ultraviolet completion. Their extension to gravity theories has been studied in the past years aiming at application to the swampland program, showing that positivity bounds hold at least approximately even in the presence of gravity. An issue in this context is how much negativity is allowed for a given scattering process. In this paper we address importance of this rather technical issue by demonstrating that it is relevant to physics within the scope of ongoing experiments, especially in the context of dark sector physics. In particular, we provide detailed analysis of dark photon scenarios as an illustrative example. This motivates further studies on gravitational positivity bounds.

Quantum fluctuations can endow spacetime with a foamy structure. In this review article we discuss our various proposals to observationally constrain models of spacetime foam. One way is to examine if the light wave-front from a distant quasar or GRB can be noticeably distorted by spacetime-foam-induced phase incoherence. As the phase fluctuations are proportional to the distance to the source, but inversely proportional to the wavelength, ultra-high energy photons from distant sources are particularly useful. We elaborate on several proposals, including the possibility of detecting spacetime foam by observing "seeing disks" in the images of distant quasars and active galactic nuclei. We also discuss the appropriate distance measure for calculating the expected angular broadening. Then we discuss our more recent work in which we investigate whether wave-front distortions on small scales (due to spacetime foam) can cause distant objects become undetectable because the phase fluctuations have accumulated to the point at which image formation is impossible. Another possibility that has recently become accessible is to use interferometers to observe cosmologically distant sources, thereby giving a large baseline perpendicular to the local wave vector over which the wave front could become corrugated and thus distorted, reducing or eliminating its fringe visibility. We argue that all these methods ultimately depend on the availability of ways (if any) to carry out proper averaging of contributions from different light paths from the source to the telescope.

In this work, we present the results of searches for signatures of dark matter decay or annihilation into Standard Model particles, and secret neutrino interactions with dark matter. Neutrinos could be produced in the decay or annihilation of galactic or extragalactic dark matter. Additionally, if an interaction between dark matter and neutrinos exists then dark matter will interact with extragalactic neutrinos. In particular galactic dark matter will induce an anisotropy in the neutrino sky if this interaction is present. We use seven and a half years of the High-Energy Starting Event (HESE) sample data, which measures neutrinos in the energy range of approximately 60 TeV to 10 PeV, to study these phenomena. This all-sky event selection is dominated by extragalactic neutrinos. For dark matter of $\sim$ 1 PeV in mass, we constrain the velocity-averaged annihilation cross section to be smaller than $10^{-23}$cm$^3$/s for the exclusive $\mu^+\mu^-$ channel and $10^{-22}$ cm$^3$/s for the $b\bar b$ channel. For the same mass, we constrain the lifetime of dark matter to be larger than $10^{28}$ s for all channels studied, except for decaying exclusively to $b\bar b$ where it is bounded to be larger than $10^{27}$ s. Finally, we also search for evidence of astrophysical neutrinos scattering on galactic dark matter in two scenarios. For fermionic dark matter with a vector mediator, we constrain the dimensionless coupling associated with this interaction to be less than 0.1 for dark matter mass of 0.1 GeV and a mediator mass of $10^{-4}~$ GeV. In the case of scalar dark matter with a fermionic mediator, we constrain the coupling to be less than 0.1 for dark matter and mediator masses below 1 MeV.