Papers for Wednesday, Aug 31 2022

We analyse the large-scale structure out to 100 Mpc around a sample of 16 confirmed fossil systems using spectroscopic information from the Sloan Digital Sky Survey Data Release 16. We compute the distance between our FGs and the centres of filaments and nodes presented in \citet{Chen2016}. We also study the density of bright galaxies, since they are thought to be good mass tracers, and the projected over densities of galaxies. Finally, we apply a FoF algorithm to detect virialised structures around our FGs, in order to have an estimate of the mass available in their surroundings. FGs are mainly found close to filaments, with a mean distance of $3.7 \pm 1.1$ R$_{200}$ and a minimum distance of 0.05 $R_{200}$. On the other hand, none of our FGs is found close to intersections, with a mean and minimum distance of $19.3 \pm 3.6$ and 6.1 $R_{200}$, respectively. There is a correlation for which FGs at higher redshifts are found in denser regions, when we use bright galaxies as tracers of the mass. At the same time, FGs with the largest magnitude gaps ($\Delta m_{12}$ > 2.5) are found in less dense environments and hosting, on average, smaller central galaxies. Our results suggest that FGs formed in a peculiar position of the cosmic web, close to filaments and far from nodes, in which their interaction with the cosmic web itself can be limited. We deduce that FGs with faint BCGs, large $\Delta m_{12}$, and low redshifts could be systems at the very last stage of their evolution. Moreover, we confirm theoretical predictions that systems with the largest magnitude gap are not massive.

High-redshift quasars ionize HeII into HeIII around them, heating the IGM in the process and creating large regions with elevated temperature. In this work, we demonstrate a method based on a convolutional neural network (CNN) to recover the spatial profile for $T_0$, the temperature at the mean cosmic density, in quasar proximity zones. We train the neural network with synthetic spectra drawn from a Cosmic Reionization on Computers simulation. We discover that the simple CNN is able to recover the temperature profile with an accuracy of $\approx 1400$ K in an idealized case of negligible observational uncertainties. We test the robustness of the CNN and discover that it is robust against the uncertainties in quasar host halo mass, quasar continuum and ionizing flux. We also find that the CNN has good generality with regard to the hardness of quasar spectra. Saturated pixels pose a bigger problem for accuracy and may downgrade the accuracy to $1700$ K in the outer parts of the proximity zones. Using our method, one could distinguish whether gas is inside or outside the HeIII region created by the quasar. Because the size of the HeIII region is closely related to the total quasar lifetime, this method has great potential in constraining the quasar lifetime on $\sim $Myr timescales.

SDSS J2222+2745, at z = 0.489, is one of the few currently known lens clusters with multiple images of a background (z = 2.801) quasar with measured time delays. We combine imaging from the Hubble Space Telescope (HST) with recent Multi Unit Spectroscopic Explorer (MUSE) spectroscopic data to securely identify 34 cluster members and 12 multiple images from 3 background sources. We measure the stellar velocity dispersions of 13 cluster galaxies, enabling an independent estimate of the contribution of the sub-halo mass component to the lens total mass. The projected total mass distribution of the lens cluster is best modelled with a single large-scale mass component, a galaxy-scale component, anchored by the MUSE kinematic information, and an external shear. The best-fit strong lensing model yields a root mean square separation between the model-predicted and observed positions of the multiple images of 0''.29. When analysing the impact of systematic uncertainties, stemming from modelling assumptions and used observables, we find that the projected total mass profile, relative weight of the sub-halo mass component, and critical lines are consistent, within the statistical uncertainties. The predicted magnification and time delay values are, instead, more sensitive to the local details of the lens total mass distribution, and vary significantly among lens models that are similarly good at reproducing the observed multiple image positions. Due to its complex morphology, the low number of point-like multiple images, and current model degeneracies, it becomes clear that additional information (from the observed surface brightness distribution of lensed sources and the measured time delays) needs to be included in the modelling of SDSS J2222+2745 for accurate and precise cosmological measurements. The full MUSE secure spectroscopic catalogue presented in this work is made publicly available.

[Abridged] Observations of protoplanetary disks suggest that they are depleted in gas-phase CO. It has been posed that gas-phase CO is chemically consumed and converted into less volatile species through gas-grain processes. Observations of interstellar ices reveal a CO$_2$ component within H$_2$O ice suggesting co-formation. The aim of this work is to experimentally verify the interaction of gas-phase CO with solid-state OH radicals above the sublimation temperature of CO. Amorphous solid water (ASW) is deposited at 15 K and followed by vacuum-UV (VUV) irradiation to dissociate H$_2$O and create OH radicals. Gas-phase CO is simultaneously admitted and only adsorbs with a short residence time on the ASW. Products in the solid state are studied with infrared spectroscopy and once released into the gas phase with mass spectrometry. Results show that gas-phase CO is converted into CO$_2$, with an efficiency of 7-27%, when interacting with VUV irradiated ASW. Between 40 and 90 K, CO$_2$ production is constant, above 90 K, O$_2$ production takes over. In the temperature range of 40-60 K, the CO$_2$ remains in the solid state, while at temperatures $\geq$ 70 K the formed CO$_2$ is released into the gas phase. We conclude that gas-phase CO reacts with solid-state OH radicals above its sublimation temperature. This gas-phase CO and solid-state OH radical interaction could explain the observed CO$_2$ embedded in water-rich ices. It may also contribute to the observed lack of gas-phase CO in planet-forming disks, as previously suggested. Our experiments indicate a lower water ice dissociation efficiency than originally adopted in model descriptions of planet-forming disks and molecular clouds. Incorporation of the reduced water ice dissociation and increased binding energy of CO on a water ice surfaces in these models would allow investigation of this gas-grain interaction to its full extend.

The investigation of the metal-poor tail in the Galactic bulge provides unique information on the early Milky Way assembly and evolution. A chemo-dynamical analysis of 17 very metal-poor stars (VMP, [Fe/H] $<-2.0$) selected from the Pristine Inner Galaxy Survey was carried out based on Gemini/GRACES spectra. The chemical abundances of $\alpha-$elements (Mg, Ca, and Ti), odd-Z (Na, K, and Sc), Fe-peak (Cr and Ni), and neutron-capture (Ba) elements are determined from model atmosphere analyses. The chemistry suggests that the majority of our stars are very similar to metal-poor stars in the Galactic halo. Orbits calculated from {\it Gaia} EDR3 imply these stars are brought into the bulge during the earliest Galactic assembly. Most of our stars have large [Na,Ca/Mg] abundances, and thus show little evidence of enrichment by pair-instability supernovae. Two of our stars (P171457, P184700) have chemical abundances compatible with second-generation globular cluster stars, suggestive of the presence of ancient and now dissolved globular clusters in the inner Galaxy. One of them (P171457) is extremely metal-poor ([Fe/H] $<-3.0$) and well below the metallicity floor of globular clusters, which supports the growing evidence for the existence of lower-metallicity globular clusters in the early Universe. A third star (P180956, [Fe/H] $\sim-2$) has low [Na,Ca/Mg] and very low [Ba/Fe] for its metallicity, which are consistent with formation in a system polluted by only one or a few low-mass supernovae. Interestingly, its orbit is confined to the Galactic plane, like other very metal-poor stars found in the literature, which have been associated with the earliest building blocks of the Milky Way.

Neutrino emission from gamma-ray bursts (GRBs) has been sought for a long time, and stringent limits on the most accredited GRB emission models have been obtained from IceCube. Multi-wavelength GRB observations of the last decades improved our knowledge of the GRB emission parameters, such as the Lorentz factor and the luminosity, which can vary from one GRB to another by several orders of magnitude. Empirical correlations among such parameters have been identified during the prompt phase, with direct implications on GRB models. In this work, we use the PSLab open-access code, developed for IceCube data analyses, to search for individual neutrino emission from the prompt and afterglow phases of selected GRBs, and for stacking emission from the ensemble of such GRBs. For the afterglow phase, we focus in particular on GRBs with X-ray flares and plateaus. While past stacking searches assumed the same GRB fluence at Earth, we present a stacking scheme based on physically motivated GRB weights. Moreover, we conceive a new methodology for the prompt phase that uses the empirical correlations to infer the GRB luminosity and Lorentz factor, when redshift measurements are not available. We do not observe any significant neutrino excess. Hence, we set constraints on the GRB neutrino fluxes and on relevant GRB parameters, including the magnetic field in the jet. Notably, the baryon loading is found to be <10 for typical GRB prompts, thus disfavoring a baryonic-dominated origin of the GRB ejecta.

The flux of high-energy astrophysical $\gamma$ rays is attenuated by the production of electron-positron pairs from scattering off of extragalactic background light (EBL). We use the most up-to-date information on galaxy populations to compute their contributions to the pair-production optical depth. We find that the optical depth inferred from $\gamma$-ray measurements exceeds that expected from galaxies at the $\sim2\sigma$ level. If the excess is modeled as a frequency-independent re-scaling of the standard contribution to the EBL from galaxies, then it is detected at the $2.7\sigma$ level (an overall $14-30\%$ increase of the EBL). If the frequency dependence of the excess is instead modeled as a two-photon decay of a dark-matter axion, then the excess is favored over the null hypothesis of no excess at the $2.1\sigma$ confidence level. While we find no evidence for a dark-matter signal, the analysis sets the strongest current bounds on the photon-axion coupling over the $8-25$ eV mass range.

Strong gravitational lensing has emerged as a promising approach for probing dark matter models on sub-galactic scales. Recent work has proposed the subhalo effective density slope as a more reliable observable than the commonly used subhalo mass function. The subhalo effective density slope is a measurement independent of assumptions about the underlying density profile and can be inferred for individual subhalos through traditional sampling methods. To go beyond individual subhalo measurements, we leverage recent advances in machine learning and introduce a neural likelihood-ratio estimator to infer an effective density slope for populations of subhalos. We demonstrate that our method is capable of harnessing the statistical power of multiple subhalos (within and across multiple images) to distinguish between characteristics of different subhalo populations. The computational efficiency warranted by the neural likelihood-ratio estimator over traditional sampling enables statistical studies of dark matter perturbers and is particularly useful as we expect an influx of strong lensing systems from upcoming surveys.

With the advent of the Square Kilometre Array Observatory (SKAO), scientists will be able to directly observe the Epoch of Reionization by mapping the distribution of neutral hydrogen at different redshifts. While physically motivated results can be simulated with radiative transfer codes, these simulations are computationally expensive and can not readily produce the required scale and resolution simultaneously. Here we introduce the Physics-Informed neural Network for reIONization (PINION), which can accurately and swiftly predict the complete 4-D hydrogen fraction evolution from the smoothed gas and mass density fields from pre-computed N-body simulation. We trained PINION on the C$^2$-Ray simulation outputs and a physics constraint on the reionization chemistry equation is enforced. With only five redshift snapshots and a propagation mask as a simplistic approximation of the ionizing photon mean free path, PINION can accurately predict the entire reionization history between $z=6$ and $12$. We evaluate the accuracy of our predictions by analysing the dimensionless power spectra and morphology statistics estimations against C$^2$-Ray results. We show that while the network's predictions are in good agreement with simulation to redshift $z>7$, the network's accuracy suffers for $z<7$ primarily due to the oversimplified propagation mask. We motivate how PINION performance can be drastically improved and potentially generalized to large-scale simulations.

We compare here the colours of simulated satellite galaxies from the Auriga project to observed satellite galaxies from the Exploration of Local VolumE Satellites (ELVES) Survey of satellite galaxies around local Milky Way analogues. Our goal is to understand the origin of the observed colour distribution and star-forming properties of ELVES satellites. We find that the satellite populations in the Auriga simulations, which was originally designed to model Milky Way-like host galaxies, resemble the populations in ELVES in their luminosity function, quenched fraction, and colour-magnitude distribution. We find that satellites transition from blue colours to red colours at the luminosity range $-15 \lesssim M_g \lesssim -12$ in both the simulations and observations and we show that this shift is driven by environmental effects in the simulations. We also demonstrate that the colour distribution in both simulations and observations can be decomposed into two populations based on their morphological type or star-forming status that are statistically distinct. In the simulations, these two populations also have statistically distinct infall time distributions. The transition in star-forming state appears to be consistent between Auriga and ELVES but in contrast with the Satellites Around Galactic Analogs (SAGA) survey. The comparison presented here seems to indicate that this tension is resolved by the improved target selection of ELVES, but there are still tensions in understanding the colours of faint galaxies, of which ELVES appears to have a significant population of faint blue satellites not recovered in Auriga.

Empirical constraints of fundamental properties of protoplanetary disks are essential for understanding planet formation and planetary properties (1,2). Carbon monoxide (CO) gas is often used to constrain disk properties (3). However, estimates show that the CO gas abundance in disks is depleted relative to expected values (4,5,6,7), and models of various disk processes impacting the CO abundance could not explain this depletion on observed 1Myr timescales (8,9,10,11,12,13,14). Here we demonstrate that surface energy effects on particles in disks, such as the Kelvin effect, that arise when ice heterogeneously nucleates onto an existing particle can efficiently trap CO in its ice phase. In previous ice formation models, CO gas was released when small ice-coated particles were lofted to warmed disk layers. Our model can reproduce the observed abundance, distribution and time evolution of gaseous CO in the four most studied protoplanetary disks (7). We constrain the solid and gaseous CO inventory at the midplane and disk diffusivities and resolve inconsistencies in estimates of the disk mass -- three crucial parameters that control planetary formation.

We scrutinize the light sail scenario of the first interstellar object (ISO) 1I/2017 U1 (`Oumuamua) by making comparisons between physical models and observational constraints. These analyses can be generalized for future surveys of `Oumuamua-like objects. The light sail goes through a drift in the interstellar space due to the magnetic field and gas atoms, which poses challenges to the navigation system. When the light sail enters the inner solar system, the sideways radiation pressure leads to a considerable non-radial displacement. The immensely high dimensional ratio and the tumbling motion could cause a light curve with an extremely large amplitude and could even make the light sail invisible from time to time. These observational features allow us to examine the light sail scenario of interstellar objects. Our results show that the drift of the freely rotating light sail in the interstellar medium is $\sim 100\,$au even if the travel distance is only 1 pc. The probability is < 1.5\% for the expected brightness modulation of the light sail to match with `Oumuamua's observed variation amplitude ($\sim$ 2.5 -- 3). In addition, the probability is 0.4% for the tumbling light-sail to be visible (brighter than V=27) in all 55 observations spread over two months after discovery. Radiation pressure could cause a larger displacement that is normal to the orbital plane for a lightsail than that for `Oumuamua. Also, the ratio of anti-solar to sideways acceleration of `Oumuamua deviates from that of the light sail by ~ 1.5{\sigma}. We suggest that `Oumuamua is unlikely a light sail. The dynamics of an intruding light sail, if exist, has distinct observational signatures, which can be quantitatively identified and analyzed with our methods in future surveys.

We report on the follow-up of the extremely bright gamma-ray burst GRB 210619B with optical polarimetry. We conducted optopolarimetric observations of the optical afterglow of GRB 210619B in the SDSS-r band in the time window ~ 5967 - 8245 seconds after the burst, using the RoboPol instrument at the Skinakas observatory. We find signs of variability of the polarization degree as well as the polarization angle during the time of observations. We also note a significant rise in polarization value and a significant change in the polarization angle towards the end of our observations. This is the first time such behavior is observed in this timescale.

This paper shows how a self-consistent dynamical model can be obtained by fitting the gravitational potential of the Milky Way to the stellar kinematics and densities from Gaia data. Using the Besancon Galaxy Model we derive a potential and the disc stellar distribution functions are computed based on three integrals of motion to model stationary stellar discs. The gravitational potential and the stellar distribution functions are built self-consistently, and then adjusted to be in agreement with the kinematics and the density distributions obtained from Gaia observations. A Markov chain Monte Carlo (MCMC) is used to fit the free parameters of the dynamical model to Gaia parallax and proper motion distributions. The fit is done on several sets of Gaia eDR3 data, widely spread in longitudes and latitudes. We are able to determine the velocity dispersion ellipsoid and its tilt for sub-components of different ages, both varying with R and z. The density laws and their radial scale lengths, for the thin and thick disc populations are also obtained self-consistently. This new model has some interesting characteristics, such as a flaring thin disc. The thick disc is found to present very distinctive characteristics from the old thin disc, both in density and kinematics. This well supports the idea that thin and thick discs were formed in distinct scenarios as the density and kinematics transition between them is found to be abrupt. The dark matter halo is shown to be nearly spherical. We also derive the Solar motion to be (10.79 $\pm$ 0.56, 11.06 $\pm$ 0.94, 7.66 $\pm$ 0.43) km/s, in good agreement with recent studies. The resulting fully self-consistent gravitational potential, still axisymmetric, is a good approximation of a smooth mass distribution in the Milky Way and can be used for further studies, including to compute orbits for real stars in our Galaxy (abridged).

Instrumental artifacts in gravitational-wave strain data can overlap with gravitational-wave detections and significantly impair the accuracy of the measured source localizations. These biases can prevent the detection of any electromagnetic counterparts to the detected gravitational wave. We present a method to mitigate the effect of instrumental artifacts on the measured source localization. This method uses inpainting techniques to remove data containing the instrumental artifact and then correcting for the data removal in the subsequent analysis of the data. We present a series of simulations using this method using a variety of signal classes and inpainting parameters which test the effectiveness of this method and identify potential limitations. We show that in the vast majority of scenarios, this method can robustly localize gravitational-wave signals even after removing portions of the data. We also demonstrate how an instrumental artifact can bias the measured source location and how this method can be used to mitigate this bias.

We investigate the evolution of cosmic voids in the Schrodinger Poisson formalism, finding wave mechanical solutions for the dynamics in a standard cosmological background with appropriate boundary conditions. We compare the results in this model to those obtained using the Zel'dovich approximation. We discuss the advantages of studying voids in general and the advantages of the Schrodinger Poisson description over other approaches. In particular, emphasizing the utility of the free particle approximation. We also discuss a dimensionless number, similar to the Reynolds number, for this system which allows our void solutions to be scaled to systems of different physical dimensions.

Direct imaging characterization of extrasolar planets is often done at low spectral resolution. We model the spectrograph for the Gemini Planet Imager upgrade (GPI 2.0) and assess the instrument's potential for allowing observers to constrain exoplanet properties through analysis of near-infrared spectra. We simulated noisy observations followed by calculations of posterior distributions from maximum likelihood comparison with the Sonora 2018 model grid. Preliminary results suggest that GPI 2.0 should allow observers to constrain temperature with sufficient accuracy, but gravity remains largely uncertain. We also explore the effects of incorporating convolution with the instrument line spread function into our simulation and compare the results with our preliminary findings.

The DESI Legacy Survey is a digital sky survey with a large footprint compared to other Earth-based surveys, covering both the Northern and Southern hemispheres. This paper shows the distribution of the spin directions of spiral galaxies imaged by DESI Legacy Survey. A simple analysis of dividing nearly 1.3$\cdot10^6$ spiral galaxies into two hemispheres shows a higher number of galaxies spinning counterclockwise in the Northern hemisphere, and a higher number of galaxies spinning clockwise in the Southern hemisphere. That distribution is consistent with previous observations, but uses a far larger number of galaxies and a larger footprint. The larger footprint allows a comprehensive analysis without the need to fit the distribution into an a priori model, making this study different from all previous analyses of this kind. Fitting the spin directions of the galaxies to cosine dependence shows a dipole axis alignment with probability of $P<10^{-5}$. The analysis is done with a trivial selection of the galaxies, as well as simple explainable annotation algorithm that does not make use of any form of machine learning, deep learning, or pattern recognition. While further work will be required, these results are aligned with previous studies suggesting the possibility of a large-scale alignment of galaxy angular momentum.

Recently, several eclipsing millisecond pulsars have been shown to experience strong and apparent weak lensing from the outflow of their ionized companions. Lensing can be a powerful probe of the ionized plasma, with the strongest lenses potentially resolving emission regions of pulsars, and understanding lensing in the `laboratory-like' conditions of an eclipsing pulsar may be analogously applied to fast radio bursts, which reside in dense, magnetized environments. We discover clear evidence of the two regimes of lensing, strong and apparent weak, and find variable dispersion measure (DM), absorption, scattering, and flux density in the original Black Widow pulsar PSR B1957+20 through an eclipse at the Arecibo Observatory at 327 MHz. We show that the flux density variations in the apparently weak lensing regime can be modeled directly from variations of DM, using geometric optics. The mean effective velocities in the ingress, $954\pm 99$ km/s, and egress $604\pm 47$ km/s cannot be explained by orbital motions alone, but are consistent with significant outflow velocity of material from the companion. We also show that geometric optics can predict when and where the lensing regime-change between weak and strong occurs, and argue that the apparent weak lensing to be due averaging of many images. Our framework can be applied to other eclipsing pulsars, which are likely lensed by intrabinary material near their eclipses to predict weak and strong lensing, and in principle, to independently constrain inclinations.

Aims. We present the implementation of the treatment of particle ejection and dust nucleation in the smoothed particle hydrodynamics (SPH) code phantom. These developments represent the first step toward a more complete modeling of dust-driven winds emanating from AGB stars. Methods. The AGB outflow is modeled by injecting the SPH particles from a spherical inner boundary. This boundary is a series of concentric shells, with the AGB star at its center, and the particles are positioned on these shells on the vertices of an isocahedron geodesic surface. The outermost shell is ejected with a predefined radial velocity, and subsequent lower shells replenish the ejected ones, all rotated randomly to improve the isotropy of the outflow. The physical properties of the particles on these shells are set by solving the 1D analytic steady wind equations. The formation of dust is calculated starting from a compact chemical network for carbon-rich material, which creates the building blocks of the solid-state particles. Subsequently, the theory of the moments is used to obtain dust growth rates, without requiring knowledge on the grain size distribution. Results. We tested our implementation against a series of 1D reference solutions. We demonstrate that our method is able to reproduce Parker-type wind solutions. For the trans-sonic solution, small oscillations are present in the vicinity of the sonic point, but these do not impact the trans-sonic passage or terminal wind velocity. Supersonic solutions always compare nicely with 1D analytic profiles. We also tested our implementation of dust using two formalisms: an analytic prescription for the opacity devised by Bowen and the full treatment of carbon-dust formation. Both simulations reproduce the 1D analytic solution displaying the expected additional acceleration when the gas temperature falls below the condensation temperature.

Particle acceleration physics at supernova remnant (SNR) shocks is one of the most intriguing problems in astrophysics. SNR RCW~86 provides a suitable environment for understanding the particle acceleration physics because one can extract the information of both accelerated particles and acceleration environment at the same regions through the bright X-ray emission. In this work, we study X-ray proper motions and spectral properties of the southwestern region of RCW~86. The proper motion velocities are found to be $\sim 300$--2000~km~s$^{-1}$ at a distance of 2.8~kpc. We find two inward-moving filaments, which are more likely reflected shocks rather than reverse shocks. Based on the X-ray spectroscopy, we evaluate thermal parameters such as the ambient density and temperature, and non-thermal parameters such as the power-law flux and index. From the flux decrease in time of several non-thermal filaments, we estimate the magnetic field amplitudes to be $\sim 30$--100~$\mu$G. Gathering the physical parameters, we then investigate parameter correlations. We find that the synchrotron emission from thermal-dominated filaments is correlated with the ambient density $n_{\rm e}$ as $\text{(power-law flux)} \propto n_{\rm e}^{1.0 \pm 0.2}$ and $\text{(power-law index)} \propto n_{\rm e}^{0.38 \pm 0.10}$, not or only weakly with the shock velocity and shock obliquity. As an interpretation, we propose a shock-cloud interaction scenario, where locally enhanced magnetic turbulence levels have a great influence on local acceleration conditions.

We present the first gri-band period-luminosity (PL) and period-Wesenheit (PW) relations for the fundamental mode anomalous Cepheids. These PL and PW relations were derived from a combined sample of five anomalous Cepheids in globular cluster M92 and the Large Magellanic Cloud, both of which have distance accurate to ~1% available from literature. Our g-band PL relation is similar to the B-band PL relation as reported in previous study. We applied our PL and PW relations to anomalous Cepheids discovered in dwarf galaxy Crater II, and found a larger but consistent distance modulus than the recent measurements based on RR Lyrae. Our calibrations of gri-band PL and PW relations, even though less precise due to small number of anomalous Cepheids, will be useful for distance measurements to dwarf galaxies.

Tidal disruption events (TDEs) occur when a star passes too close to a supermassive black hole and is destroyed by tidal gravitational forces. Radio observations of TDEs trace synchrotron emission from outflowing material that may be ejected from the inner regions of the accretion flow around the SMBH or by the tidal debris stream. Radio detections of tidal disruption events are rare, but provide crucial information about the launching of jets and outflows from supermassive black holes and the circumnuclear environment in galaxies. Here we present the radio detection of the TDE AT2020opy, including three epochs of radio observations taken with the Karl G. Jansky's Very Large Array (VLA), MeerKAT, and upgraded Giant Metrewave Radio telescope. AT2020opy is the most distant thermal TDE with radio emission reported to date, and from modelling the evolving synchrotron spectra we deduce that the host galaxy has a more dense circumnuclear medium than other thermal TDEs detected in the radio band. Based on an equipartition analysis of the synchrotron spectral properties of the event, we conclude that the radio-emitting outflow was likely launched approximately at the time of, or just after, the initial optical flare. We find no evidence for relativistic motion of the outflow. The high luminosity of this event supports that a dense circumnuclear medium of the host galaxy produces brighter radio emission that rises to a peak more quickly than in galaxies with lower central densities.

The total available sample of fast radio bursts (FRBs) has been growing steadily in recent years, facilitating the study of FRBs from a statistical point of view. At the same time, the classification of FRBs is currently an imperative issue. We propose that the brightness temperature of bursts can serve as an ideal criterion for classification. In this work, we gather the available data for all localized FRBs and we find a positive relation between the intrinsic pulse width and burst energy, $T_{\rm i}\propto E_\nu^{0.25}$, for three repeating FRBs that is similar to that of our previous work using FRB 20121102A data alone. The critical line $T_{\rm B,cri}$ is found to vary for different FRBs, which may reflect the differences in source properties. This relation can put strong constraints on mainstream radiation mechanisms. It is evident that neither the coherent curvature radiation or synchrotron maser radiation have the capability to reach the high brightness temperature required to reproduce this relation.

Centers of starburst galaxies may be characterized by a specific gas and ice chemistry due to their gas dynamics and the presence of various ice desorption mechanisms. This may result in a peculiar observable composition. We analyze abundances of $CO_2$, a reliable tracer of ice chemistry, from data collected as part of the ALMA large program ALCHEMI, a wide-frequency spectral scan toward the starburst galaxy NGC~253 with an angular resolution of 1.6$''$. We constrain the $CO_2$ abundances in the gas phase using its protonated form $HOCO^+$. The distribution of $HOCO^+$ is similar to that of methanol, which suggests that $HOCO^+$ is indeed produced from the protonation of $CO_2$ sublimated from ice. The $HOCO^+$ fractional abundances are found to be $(1-2)\times10^{-9}$ at the outer part of the central molecular zone (CMZ), while they are lower ($\sim10^{-10}$) near the kinematic center. This peak fractional abundance at the outer CMZ is comparable to that in the Milky Way CMZ, and orders of magnitude higher than that in Galactic disk star-forming regions. From the range of $HOCO^+/CO_2$ ratios suggested from chemical models, the gas-phase $CO_2$ fractional abundance is estimated to be $(1-20)\times10^{-7}$ at the outer CMZ, and orders of magnitude lower near the center. We estimate the $CO_2$ ice fractional abundances at the outer CMZ to be $(2-5)\times10^{-6}$ from the literature. A comparison between the ice and gas $CO_2$ abundances suggests an efficient sublimation mechanism. This sublimation is attributed to large-scale shocks at the orbital intersections of the bar and CMZ.

The magnetic field governs the corona; hence it is a crucial parameter to measure. Unfortunately, existing techniques for estimating its strength are limited by strong assumptions and limitations. These techniques include photospheric or chromospheric field extrapolation using potential or non-linear-force-free methods, estimates based on coronal seismology, or by direct observations via, e.g., the Cryo-NIRSP instrument on DKIST which will measure the coronal magnetic field, but only off the limb. Alternately, in this work we investigate a recently developed approach based on the magnetic-field-induced (MIT) transition of the \fex~257.261~\AA. In order to examine this approach, we have synthesized several \fex\ lines from two 3D magnetohydrodynamic simulations, one modeling an emerging flux region and the second an established mature active region. In addition, we take bound-free absorption from neutral hydrogen and helium and singly ionised helium into account. The absorption from cool plasma that occurs at coronal heights has a significant impact on determining the magnetic field. We investigate in detail the challenges of using these \fex\ lines to measure the field, considering their density and temperature dependence. We present a novel approach to deriving the magnetic field from the MIT using inversions of the differential emission measure as a function of the temperature, density, and magnetic field. This approach successfully estimates the magnetic field strength (up to \%18 relative error) in regions that do not suffer from significant absorption and that have relatively strong coronal magnetic fields ($>250$~G). This method allows the masking of regions where absorption is significant.

Gamma-ray bursts (GRBs) are intense bursts of high-energy photons (prompt emissions) caused by relativistic jets. After the emissions, multi-wavelength afterglows, from radio to very-high-energy (VHE) gamma-rays, last for more than a few days. In the past three years, the VHE gamma-ray photons from four GRBs (GRBs 180720B, 190114C, 190829A and 201216C) were detected by ground-based Imaging Atmospheric Cherenkov Telescopes, such as the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) telescopes and the High Energy Stereoscopic System (H.E.S.S.). One of them, GRB 190829A, had some peculiar features of showing achromatic peaks in X-ray and optical bands at $1.4\times10^3$ s and being classified as low-luminosity GRBs. Previously, we proposed a two-component jet model, which has `narrow jet' with a small initial jet opening half-angle $\theta_0(=0.015$ rad) and large bulk Lorentz factor $\Gamma_0=350$ and `wide jet' with $\theta_0=0.1$ rad and $\Gamma_0=20$. The narrow jet explained the early X-ray and optical emissions and apparently small isotropic gamma-ray energy and peak energy in the off-axis viewing case. Furthermore, the late X-ray and radio (1.3 and 15.5 GHz) afterglows were emitted from the wide jet. Here, we calculate the VHE gamma-ray flux by the synchrotron self-Compton (SSC) emission. The multi-wavelength afterglows of GRB 190829A including the VHE gamma-ray emission are well explained by our two-component jet model. The afterglow emissions from our two-component jet are also consistent with the observational results of GRBs 180720B, 190114C and 201216C. Furthermore, we discuss the detectability of off-axis orphan afterglows by the Cherenkov Telescope Array (CTA).

Applying functional differentiation to the density field with Newtonian gravity, we obtain the static, nonlinear equation of the three-point correlation function $\zeta$ of galaxies, to the third order density perturbations. We make the equation closed and perform renormalization of the mass and the Jeans wavenumber. Using the boundary condition inferred from observations, we obtain the third order solution $\zeta(r, u, \theta)$ at fixed $u=2$, which is positive, exhibits a $U$-shape along the angle $\theta$, and decreases monotonously along the radial $r$ up to the range $r \leq 30\, h^{-1}$Mpc in our computation. The corresponding reduced $Q(r, u, \theta)$ deviates from 1 of the Gaussian case, has a deeper $U$-shape along $\theta$, and varies non-monotonously along $r$. The third order solution agrees with the SDSS data of galaxies, quite close to the previous second order solution, especially at large scales. This indicates that the equations of correlation functions with increasing orders of density perturbation provide a stable description of the nonlinear galaxy system.

We present a set of peculiar radio sources detected using an unsupervised machine learning method. We use data from the Australian Square Kilometre Array Pathfinder (ASKAP) telescope to train a self-organizing map (SOM). The radio maps from three ASKAP surveys, Evolutionary Map of Universe pilot survey (EMU-PS), Deep Investigation of Neutral Gas Origins pilot survey (DINGO) and Survey With ASKAP of GAMA-09 + X-ray (SWAG-X), are used to search for the rarest or unknown radio morphologies. We use an extension of the SOM algorithm that implements rotation and flipping invariance on astronomical sources. The SOM is trained using the images of all "complex" radio sources in the EMU-PS which we define as all sources catalogued as "multi-component". The trained SOM is then used to estimate a similarity score for complex sources in all surveys. We select 0.5\% of the sources that are most complex according to the similarity metric, and visually examine them to find the rarest radio morphologies. Among these, we find two new odd radio circle (ORC) candidates and five other peculiar morphologies. We discuss multiwavelength properties and the optical/infrared counterparts of selected peculiar sources. In addition, we present examples of conventional radio morphologies including: diffuse emission from galaxy clusters, and resolved, bent-tailed, and FR-I and FR-II type radio galaxies. We discuss the overdense environment that may be the reason behind the circular shape of ORC candidates.

We probe four cosmological models which, potentially, can solve the Hubble tension according to the dark energy equation of state. In this context, we demonstrate that the Einstein Telescope is capable of achieving a relative accuracy below $1\%$ on the Hubble constant independently of the specific dark energy model. We firstly build mock catalogs containing gravitational wave events for one, five and ten years of observations, and above Signal-to-Noise Ratio equal to nine. From these catalogs, we extract the events which are most likely associated with possible electromagnetic counterpart detected by THESEUS. Finally, we select four dark energy models, namely a non-flat $\omega$CDM, an interacting dark energy, an emergent dark energy, and a time varying gravitational constant model, to forecast the precision down to which the Einstein Telescope can bound the corresponding cosmological parameters. We foresee that the Hubble constant is always constrained with less than $1\%$ uncertainty, thereby offering a potential solution to the Hubble tension. The accuracy on the other cosmological parameters is at most comparable with the one currently obtained using multiple probes, except for the emergent dark energy model for which the Einstein Telescope alone will be able to improve the current limits by more than one order of magnitude.

Environment plays a critical role in the star formation history of galaxies. Tidal and hydrodynamical stripping, prominent in cluster environment, can remove the peripheral gas of galaxies and star formation may thus be environmentally suppressed from the outside-in. We revisit the environmental dependence of the radial gradient of specific star formation rate (sSFR) profile. We probe the radial gradient by using the archival spectral indices D4000n and HdA measured from SDSS fiber spectra, to indicate central sSFR, and the total sSFR from fitting the spectral energy distribution. Despite the low spatial resolution, the wealth of SDSS data allows to disentangle the dependences on stellar mass, sSFR, and environment. We find that low-mass satellite galaxies in the mass range 9 < log M/M_solar < 9.8 on average quench in more inside-out pattern compared to isolated galaxies matched in mass, sSFR, and fiber coverage. This environmental effect is particularly strong for galaxies below the star formation main sequence, and peaks for those in the core of massive clusters where the phase-space diagram reveals clear links between the inside-out quenching and orbital properties. Our results suggest that both tidal and hydrodynamical interactions in cluster environment suppress the star formation of satellites mainly from the inside-out. As accreted gas of low angular momentum from hot gas halos is an important source for replenishing central gas reservoir, we discuss how gas stripping in clusters may lead to starvation and cause inside-out quenching when the outer star-forming discs are not significantly affected.

The three-body problem is famously chaotic, with no closed-form analytical solutions. However, hierarchical systems of three or more bodies can be stable over indefinite timescales. A system is considered hierarchical if the bodies can be divided into separate two-body orbits with distinct time- and length-scales, such that one orbit is only mildly affected by the gravitation of the other bodies. Previous work has mapped the stability of such systems at varying resolutions over a limited range of parameters, and attempts have been made to derive analytic and semi-analytic stability boundary fits to explain the observed phenomena. Certain regimes are understood relatively well. However, there are large regions of the parameter space which remain un-mapped, and for which the stability boundary is poorly understood. We present a comprehensive numerical study of the stability boundary of hierarchical triples over a range of initial parameters. Specifically, we consider the mass ratio of the inner binary to the outer third body ($q_{\rm out}$), mutual inclination ($i$), initial mean anomaly and eccentricity of both the inner and outer binaries ($e_{\rm in}$ and $e_{\rm out}$ respectively). We fit the dependence of the stability boundary on $q_{\rm out}$ as a threshold on the ratio of the inner binary's semi-major axis to the outer binary's pericentre separation $a_{\rm in}/R_{\rm p, out} \leq 10^{-0.6 + 0.04q_{\rm out}}q_{\rm out}^{0.32+0.1q_{\rm out}}$ for coplanar prograde systems. We develop an additional factor to account for mutual inclination. The resulting fit predicts the stability of $10^4$ orbits randomly initialised close to the stability boundary with $87.7\%$ accuracy.

Previous observations suggested that Ap and Bp stars exhibit a bimodal distribution of surface magnetic field strengths and that actually only few or no stars exist with magnetic dipole field strengths below 300 G down to a few Gauss. As the number of Ap and Bp stars currently known to possess weak magnetic fields is not large, it is necessary to carry out additional spectropolarimetric studies of Ap and Bp stars to prove whether the assumption of the existence of a critical value for the stability of magnetic fields is realistic. In this study, we present high-resolution HARPSpol magnetic field measurements for a sample of Ap stars with sharp spectral lines with a view to characterize the strengths of their magnetic fields. Out of the studied seven sharp-lined stars, two stars, HD 174779 and HD 203932, exhibit a rather weak longitudinal magnetic field with $\left< B_{\rm z}\right>=-45\pm3$ G and $\left< B_{\rm z}\right>=21\pm4$ G, respectively. Additionally, TESS observations were used to test previous conclusions on the differentiation of rotation periods of Ap and Bp stars. Apart from HD 189832 and HD 203932, all other studied sharp-lined stars have long rotation periods. Since an explanation for the slow rotation of Ap stars is currently missing, additional studies of slowly rotating Ap and Bp stars are necessary to improve our understanding of the formation and evolution of Ap and Bp stars.

The combination of multi-band imaging from HST with MUSE integral field spectroscopy, obtained at the VLT, has recently driven remarkable progress in strong lensing (SL) modeling of galaxy clusters. From a few tens of multiple images with photometric redshifts per cluster, a new generation of high-precision SL models have recently been developed, by exploiting in some cases over a hundred of spectroscopically confirmed multiple images and cluster member galaxies. A further step forward is expected with JWST observations of SL clusters (from hundreds to possibly a thousand of multiple images). In this context, we present a new, state-of-the-art SL model of the galaxy cluster MACS J0416.1-2403, utilizing 237 spectroscopically confirmed multiple images, which is the largest sample of secure multiply lensed sources utilized to date. This model incorporates stellar kinematics information of 64 cluster galaxies and the hot-gas mass distribution of the cluster determined from Chandra X-ray observations. The observed positions of the many multiple images are reproduced with a remarkable accuracy of 0.43 arcsec. To further assess the reliability of this lens model and to highlight the improvement over previously published models, we show the extended surface brightness reconstruction of several lensed galaxies through a newly developed forward modeling software. The comparison with other SL models of the same cluster demonstrates that this new model is better suited to accurately reproduce the positions, shapes and fluxes of the observed multiple images. Besides a robust characterization of the total mass distribution of the cluster, our model can provide accurate and precise magnification maps that are key to studying the intrinsic physical properties of faint, high-redshift lensed sources. The model is made publicly available through our newly developed Strong Lensing Online Tool (or SLOT).

Compact binary mergers involving neutron stars can eject a fraction of their mass to space. Being extremely neutron rich, this material undergoes rapid neutron capture nucleosynthesis, and the resulting radioactivity powers fast, short-lived electromagnetic transients known as kilonova or macronova. Such transients are exciting probes of the most extreme physical conditions and their observation signals the enrichment of the Universe with heavy elements. Here we review our current understanding of the mass ejection mechanisms, the properties of the ejecta and the resulting radioactive transients. The first well-observed event in the aftermath of GW170817 delivered a wealth of insights, but much of today's picture of such events is still based on a patchwork of theoretical studies. Apart from summarizing the current understanding, we also point out questions where no consensus has been reached yet, and we sketch possible directions for the future research. In an appendix, we describe a publicly available heating rate library based on the WinNet nuclear reaction network, and we provide a simple fit formula to alleviate the implementation in hydrodynamic simulations.

We have studied the kinematics of Galactic masers and radio stars with measured VLBI trigonometric parallaxes and proper motions. We have considered masers with relative trigonometric parallax errors less than 10\% and determined the Galactic rotation parameters from them. In particular, the linear rotation velocity of the Galaxy at the solar distance $R_0$ has been found to be $244.4\pm4.3$ km s$^{-1}$ (for the adopted $R_0=8.1\pm0.1$ kpc). We have performed a joint and separate spectral analysis of the radial, residual tangential, and vertical maser velocities. For example, from the vertical maser velocities we have estimated the velocity perturbation amplitude $f_W=5.2\pm1.5$ km s$^{-1}$ with the wavelength $\lambda_W=2.6\pm0.7$ kpc, arguing for the influence of the spiral density wave on the vertical stellar velocities. Based on 104 masers within 3 kpc of the Sun, as a result of the joint solution, we have estimated the radial, $f_R=6.7\pm1.1$ km s$^{-1}$, and tangential, $f_\theta=2.6\pm1.2$ km s$^{-1}$, velocity perturbations, the perturbation wavelength $\lambda=2.1\pm0.3$ kpc, and the Sun's phase in the Galactic spiral density wave $\chi_\odot=-148\pm15^\circ$. We have confirmed the presence of the Radcliffe wave in the spatial distribution of masers and radio stars belonging to the Local Arm.

A well-known precursor of an imminent solar eruption is the appearance of a hot S-shaped loop, also known as sigmoid, in an active region (AR). Classically, the formation of such an S-shaped loop is envisaged to be implemented by magnetic reconnection of two oppositely oriented J-shaped loops. However, the details of reconnection are elusive due to weak emission and subtle evolution during the pre-eruptive phase. In this paper, we investigate how a single J-shaped loop transforms into an S-shaped one through the slippage of one of its footpoints in NOAA AR 11719 on 2013 April 11. During an interval of about 16 min, the J-shaped loop slips through a low-corona region of strong electric current density in a bursty fashion, reaching a peak apparent speed as fast as over 1000 km s$^{-1}$, at the slipping footpoint. The enhancement of electric current density, as suggested by non-linear force-free field modeling, indicates that the "non-idealness" of coronal plasma becomes locally important, which may facilitate magnetic reconnection. The loop segment undergoing slipping motions is heated; meanwhile, above the fixed footpoint coronal emission dims due to a combination effect of the lengthening and heating of the loop, the latter of which is manifested in the temporal variation of dimming slope and of emission measure. These features together support an asymmetric scenario of sigmoid formation through slipping reconnection of a single J-shaped loop, which differs from the standard tether-cutting scenario involving a double J-shaped loop system.

Spectral synthesis is a powerful tool for interpreting the physical properties of galaxies by decomposing their spectral energy distributions into the main luminosity contributors (e.g. stellar populations or ionised gas). However, the impact nebular emission has on the inferred properties of star-forming (SF) galaxies has been largely overlooked over the years. The objective of this work is to estimate the relations between stellar properties of SF galaxies from SDSS DR7 by simultaneously fitting the stellar and nebular continua with FADO and comparing them to the results derived using STARLIGHT, a representative of purely stellar population synthesis codes. Differences between codes regarding average mass, mean age and mean metallicity values can go as high as $\sim$0.06 dex for the overall population of galaxies and $\sim$0.12 dex for SF galaxies (galaxies with EW(H$\alpha$)>3 \AA), with the most prominent difference between both codes in the light-weighted mean stellar age. A closer look into the average light- and mass-weighted star formation histories of intensively SF galaxies (EW(H$\alpha$)>75 \AA) suggests that STARLIGHT is underestimating the average light-weighted age of intensively SF galaxies by up to $\sim$0.17 dex and overestimating the light-weighted metallicity by up to $\sim$0.13 dex compared to FADO (or vice versa). The comparison between the average stellar properties of passive, SF and intensively SF galaxy samples also reveals that differences between codes increase with increasing EW(H$\alpha$) and decreasing total stellar mass. This work finds indirect evidence that a purely stellar population synthesis approach negatively impacts the inferred stellar properties of galaxies with relatively high star formation rates. In turn, this can bias interpretations of fundamental relations such as the mass-age or mass-metallicity.

The bubble size distribution is a summary statistics that can be computed from the observed 21-cm signal from the Epoch of Reionization. As it depends only on the ionization field and is not limited to gaussian information, it is an interesting probe, complementary to the power spectrum of the full 21-cm signal. Devising a flexible and reliable theoretical model for the bubble size distribution paves the way for using it for astrophysical parameters inference. The proposed model is built from the excursion set theory and a functional relation between the bubble volume and the collapsed mass in the bubble. Unlike previous models it accommodates any functional relation or distributions. Using parameterized relations allows us to test the predictive power of the model by performing a minimization best-fit to the bubble size distribution obtained from a high resolution, fully coupled radiative hydrodynamics simulations, HIRRAH-21. Our model is able to provide a better fit to the numerical bubble size distribution at ionization fraction of $x_{\text{H}_{\text{II}}} \sim 1\%$ and $3\%$ than other existing models. Moreover, the bubble volume to collapsed mass relation corresponding to the best-fit parameters, which is not an observable, is compared to numerical simulation data. A good match is obtained, confirming the possibility to infer this relation from an observed bubble size distribution using our model. Finally we present a simple algorithm that empirically implements the process of percolation. We show that it extends the usability of our bubble size distribution model up to $x_{\text{H}_{\text{II}}} \sim 30\%$.

We investigated the effect of magnetic fields on the collision process between dense molecular cores. We performed three-dimensional magnetohydrodynamic simulations of collisions between two self-gravitating cores using the Enzo adaptive mesh refinement code. The core was modeled as a stable isothermal Bonnor-Ebert (BE) sphere immersed in uniform magnetic fields. Collisions were characterized by the offset parameter $b$, Mach number of the initial core $\mathcal{M}$, magnetic field strength $B_{0}$, and angle $\theta$ between the initial magnetic field and collision axis. For head-on ($b = 0$) collisions, one protostar was formed in the compressed layer. The higher the magnetic field strength, the lower the accretion rate. For models with $b = 0$ and $\theta = 90^{\circ}$, the accretion rate was more dependent on the initial magnetic field strength compared with $b = 0$ and $\theta = 0^{\circ}$ models. For off-center ($b = 1$) collisions, the higher specific angular momentum increased; therefore, the gas motion was complicated. In models with $b = 1$ and $\mathcal{M} = 1$, the number of protostars and gas motion highly depended on $B_{0}$ and $\theta$. For models with $b = 1$ and $\mathcal{M} = 3$, no significant shock-compressed layer was formed and star formation was not triggered.

(Abridged) We study the distribution of galaxies across the ageing diagram (AD): the connection between the fraction of stars formed during the last few Myr -- traced by EW(H$\alpha$) -- and the amount of stellar mass formed on scales of ~Gyr -- using the dust-corrected optical colour $(g-r)_0$ or the Balmer break $D_n(4000)$ as observational proxies. By means of this diagram, it is possible to track recent changes on the star formation history (SFH) of galaxies and discriminate between galaxies governed by secular evolution (ageing) and systems whose star formation was suddenly interrupted (quenching). In order to provide a physical characterization of the location of galaxies across the AD, we use Pipe3D estimates of the SFH of two galaxy samples selected from the CALIFA and MaNGA surveys, in combination with the predictions from IllustrisTNG. We find that the distribution of galaxies in the AD can be described by four regimes, that strongly correlate with the stellar mass fractions formed in the last 300 Myr and 3 Gyr. Ageing systems, whose star formation rate changes on scales comparable to the age of the Universe, account for $70-80\%$ of the galaxy population in all our samples. Objects whose SFH was abruptly truncated arrange along a detached Quenched sequence that represents $\sim 5-10\%$ by (volume-corrected) number. Undetermined systems represent an intermediate population between the Ageing and Quenched regimes. Finally, Retired galaxies are located at the region in the AD where the Ageing and Quenched sequences converge. We define different star formation activity levels in terms of the birth rate parameter $b\equiv \frac{SFR}{\langle SFR \rangle}$ and find that galaxies transit from the Ageing to the Quenched sequence on scales ~500 Myr. We conclude that the ageing diagram provides a useful tool to discern recently Quenched galaxies from the dominant Ageing population.

We investigate the curvature perturbations induced by the randomness of the quantum tunneling process during cosmological first-order phase transitions (PTs) and for the first time ultilize curvature perturbations to constrain the PT parameters. We find that the observations of the cosmic microwave background spectrum distortion and the ultracompact minihalo abundance can give strict constraints on the PTs below 100GeV, especially for the low-scale PTs and the weak PTs. The current constraint on the PT parameters is largely extended by the results in this work.

Recent debates around the testability of the inflationary paradigm raise the question of how to model-independently discriminate it from competing scenarios. We argue that a detection of the Cosmic Graviton Background (CGB), the relic radiation from gravitons decoupling around the Planck time, would rule out the inflationary paradigm, as realistic inflationary models would dilute the CGB to an unobservable level. The CGB contribution to the effective number of relativistic species, $\Delta N_{{\rm eff},g} \approx 0.054$, is well within the reach of next-generation cosmological probes. We argue that detecting the high-frequency stochastic gravitational wave background associated to the CGB will be challenging but potentially feasible. We briefly discuss expectations within alternatives to inflation, focusing on bouncing cosmologies and emergent scenarios.

We analyze the reasons for the correlation between the temperature, direction, and speed of the interstellar neutral gas inflow into the heliosphere, obtained in analyzes of observations performed by the IBEX-Lo instrument onboard Interstellar Boundary Explorer (IBEX). We point out that this correlation is the combined result of the inability to measure the speed of the atoms that enter the instrument and the restriction of the observations to a short orbital arc around the Sun performed by the instrument during observation. We demonstrate that without the capability to measure the speed, but with the ability to perform observations along longer orbital arcs, or from at least two distant locations on the orbit around the Sun, it is possible to break the parameter correlation. This, however, requires a capability to adjust the boresight of the instrument relative to the spacecraft rotation axis, such as that of the planned IMAP-Lo camera onboard the Interstellar Mapping and Acceleration Probe (IMAP).

While the successful launch and operation in space of the Gas Pixel Detectors onboard the PolarLight cubesat and the Imaging X-ray Polarimetry Explorer demonstrate the viability and the technical soundness of this class of detectors for astronomical X-ray polarimetry, it is clear that the current state of the art is not ready to meet the challenges of the next generation of experiments, such as the enhanced X-ray Timing and Polarimetry mission, designed to allow for a significantly larger data throughput. In this paper we describe the design and test of a new custom, self-triggering readout ASIC, dubbed XPOL-III, specifically conceived to address and overcome these limitations. While building upon the overall architecture of the previous generations, the new chip improves over its predecessors in several, different key areas: the sensitivity of the trigger electronics, the flexibility in the definition of the readout window, as well as the maximum speed for the serial event readout. These design improvements, when combined, allow for almost an order of magnitude smaller dead time per event with no measurable degradation of the polarimetric, spectral, imaging or timing capability of the detector, providing a good match for the next generation of X-ray missions.

We show that measuring the proper motion of ${{\sim 2000}}$ stars within a dwarf galaxy, with an uncertainty of 1 km/s at most, can establish whether the Dark Matter (DM) density profile of the dwarf has a central core or cusp. We derive these limits by building mock star catalogues similar to those expected from future astrometric {\it Theia}-like missions and including celestial coordinates, radial velocity and proper motion of the stars. The density field of the DM halo of the dwarf is sampled from an extended Navarro-Frank-White (eNWF) spherical model, whereas the number density distribution of the stars is a Plummer sphere. The velocity field of the stars is set according to the Jeans equations. A Monte Carlo Markov Chain algorithm applied to a sample of $N\gtrsim 2000$ stars returns unbiased estimates of the eNFW DM parameters within $10\%$ of the true values and with $1\sigma$ relative uncertainties $\lesssim 20$\%. The proper motions of the stars lift the degeneracy among the eNFW parameters which appears when the line-of-sight velocities alone are available. {Our analysis demonstrates that, by estimating the log-slope of the mass density profile estimated at the half-light radius, a sample of $N=2000$ stars can distinguish between a core and a cusp at more than $8\sigma$.} Proper motions also return unbiased estimates of the dwarf mass profile with $1\sigma$ uncertainties that decrease, on average, from 2.65 dex to 0.15 dex when the size of the star sample increases from $N=100$ to $N=6000$ stars. The measure of the proper motions can thus strongly constrain the distribution of DM in nearby dwarfs and provides a fundamental contribution to understanding the nature and the properties of DM.

We measure optical colors for the bulges of 312 disk galaxies from the Carnegie-Irvine Galaxy Survey and convert their previously available $R$-band structural parameters to stellar mass parameters. We also measure their average stellar mass surface density in the central 1 kpc ($\Sigma_{1}$). Comparing the mass-based Kormendy relation with the original one based on flux, we find that the majority of the classifications into classical and pseudo bulges, as well as their overall statistical properties, remain essentially unchanged. While the bulge type classifications of the Kormendy relation are robust against stellar population effects, the mass-based classification criteria do produce better agreement between bulge structural properties and their stellar populations. Moreover, the mass-based Kormendy relation reveals a population of ultra-dense bulges akin to high-$z$ compact early-type galaxies, which are otherwise hidden in the original Kormendy relation. These bulges are probably relics of spheroids assembled in the early Universe, although for some we cannot rule out some contribution from secular growth. We confirm previous studies that $\Sigma_1$ correlates well with bulge surface densities.

With the goal of a thorough investigation of the ionised gas and its origin in early-type group-dominant galaxies, we present archival MUSE data for 18 galaxies from the Complete Local-Volume Groups Sample (CLoGS). This data allowed us to study the spatially-resolved warm gas properties, including the morphology of the ionised gas, EW(H$\alpha$) and kinematics as well as the gas-phase metallicity (12 + log(O/H)) of these systems. In order to distinguish between different ionisation mechanisms, we used the emission-line ratios [O III]/H$\beta$ and [N II]/H$\alpha$ in the BPT diagrams and EW(H$\alpha$). We find that the ionisation sources in our sample have variable impacts at different radii, central regions are more influenced by low-luminosity AGN, while extended regions of LINER-like emission are ionised by other mechanisms with pAGBs photoionisation likely contributing significantly. We classified our sample into three H$\alpha$+[N II] emission morphology types. We calculate the gas-phase metallicity assuming several methods and ionisation sources. In general, 12 + log(O/H) decreases with radius from the centre for all galaxies, independently of nebular morphology type, indicating a metallicity gradient in the abundance profiles. Interestingly, the more extended filamentary structures and all extranuclear star-forming regions present shallow metallicity gradients. Within the uncertainties these extended structures can be considered chemically homogeneous. We suggest that group-dominant galaxies in our sample likely acquired their cold gas in the past as a consequence of one or more mechanisms, e.g. gas-clouds or satellite mergers/accretion and/or cooling flows that contribute to the growth of the ionised gas structures.

We study off-limb emission of the lower solar atmosphere using high-resolution imaging spectroscopy in the H$\beta$ and Ca II 8542 \r{A} lines obtained with the CHROMospheric Imaging Spectrometer (CHROMIS) and the CRisp Imaging SpectroPolarimeter (CRISP) on the Swedish 1-m Solar Telescope. The H$\beta$ line wing images show the dark intensity gap between the photospheric limb and chromosphere which is absent in the Ca II images. We calculate synthetic spectra of the off-limb emissions with the RH code in the one-dimension spherical geometry and find good agreement with the observations. The analysis of synthetic line profiles shows that the gap in the H$\beta$ line wing images maps the temperature minimum region between the photosphere and chromosphere due to the well known opacity and emissivity gap of Balmer lines in this layer. However, observed gap is detected farther from the line core in the outer line wing positions than in the synthetic profiles. We found that an increased microturbulence in the model chromosphere is needed to reproduce the dark gap in the outer line wing, suggesting that observed H$\beta$ gap is the manifestation of the temperature minimum and the dynamic nature of the solar chromosphere. The temperature minimum produces a small enhancement in synthetic Ca II line-wing intensities. Observed off-limb Ca II line-wing emissions show similar enhancement below temperature minimum layer near the edge of the photospheric limb.

A significant number of double white dwarfs (DWDs) are believed to merge within the Hubble time due to the gravitational wave (GW) emission during their inspiraling phase. The outcome of a DWD system is either a type Ia Supernova as the double-degenerate model, or a massive, long-lasting merger remnant. Expected multi-messenger signals of these events will help us to distinguish detailed merging physical processes. In this work, we aim to provide a generic scenario of DWD merging, investigate the emission of all major messengers, with a focus on GWs and neutrinos. Our goal is to provide some guidance for current and future (collaborative) efforts of multi-messenger observations. Throughout the merging evolution of a DWD system, different messengers (GW, neutrino and electromagnetic wave) will dominate at different times. In this work, we show that DWD merger events located at the distance of 1 kpc can indeed produce detectable signals of GWs and neutrinos. The GW frequency are in 0.8-1.8 Hz band around 10 days before tidal disruption begin. We estimate that in optimistic situations, the neutrino number detected by upcoming detectors such as JUNO and Hyper-Kamiokande can reach O(1) for a DWD merging event at $\sim$ 1 kpc.

We present the first set of constant rest-mass ultra-massive oxygen/neon white dwarf cooling tracks with masses larger than 1.29 Msun which fully take into account the effects of general relativity on their structural and evolutionary properties. We have computed the full evolution sequences of 1.29, 1.31, 1.33, 1.35, and 1.369 Msun white dwarfs with the La Plata stellar evolution code, LPCODE. For this work, the standard equations of stellar structure and evolution have been modified to include the full effects of general relativity. For comparison purposes, the same sequences have been computed but for the Newtonian case. According to our calculations, the evolutionary properties of the most massive white dwarfs are strongly modified by general relativity effects. In particular, the resulting stellar radius is markedly smaller in the general relativistic case, being up to 25% smaller than predicted by the Newtonian treatment for the more massive ones. We find that oxygen/neon white dwarfs more massive than 1.369 Msun become gravitationally unstable with respect to general relativity effects. When core chemical distribution due to phase separation on crystallization is considered, such instability occurs at somewhat lower stellar masses, greater than 1.36 Msun. In addition, cooling times for the most massive white dwarf sequences result in about a factor of two smaller than in the Newtonian case at advanced stages of evolution. Finally, a sample of white dwarfs has been identified as ideal candidates to test these general relativistic effects. We conclude that the general relativity effects should be taken into account for an accurate assessment of the structural and evolutionary properties of the most massive white dwarfs.

In the current multi-messenger astronomy era, it is important that information about joint gravitational wave (GW) and electromagnetic (EM) observations through short gamma-ray burst (sGRBs) remains easily accessible. The possibility for non-experts to execute quick computations of joint GW-sGRB detections should be facilitated. We construct a sGRB model and add this to the framework of the previously-built Gravitational Wave Universe Toolbox. We provide expected joint GW-sGRB detection rates for different combinations of GW detectors and high-energy (HE) instruments. We employ and adapt a generic GRB model to create a top-hat jet model suitable for the Toolbox. We simulate a population of binary neutron stars (BNSs) observed by a user-specified GW detector. Our model predicts the properties of a resulting sGRB, as well as its detectability. We report predicted joint detection rates for combinations of GW detectors with HE instruments. Our findings stress the significance of the impact of the Einstein Telescope (ET); ET will observe BNSs at such a rate that the vast majority of detected sGRBs will have an observed GW counterpart. Additionally, given the limited LIGO horizon, a search for sub-threshold GW signals at higher redshifts using sGRB information from HE detectors has the potential to be very successful. Equivalently, during the ET era, GW data can assist in finding sub-threshold sGRBs, potentially increasing e.g. the number of joint ET-Fermi/GBM observations by $\sim$270%. Lastly, we find that our top-hat jet model underestimates the number of joint detections that include an off-axis sGRB. We correct for this with a second, wider and weaker jet component. We find that the majority of joint detections during the current era will include an off-axis sGRB, making GRB170817A not as unlikely as one would think. In the ET era, most joint detections will contain an on-axis sGRB.

Superconducting Transition-Edge Sensor (TES) bolometers are used for cosmic microwave background (CMB) observations. We used a testbed to evaluate the thermal performance of TES bolometers in regard to the saturation power Psat and intrinsic thermal time constant tau0. We developed an evaluation method that is complementary to methods with electrothermal feedback. In our method, the antenna termination resistor of the bolometer is directly biased with DC or AC electric power to simulate optical power, and the TES is biased with small power, which allows Psat and tau0 to be determined without contribution from the negative electrothermal feedback. We describe the method and results of the measurement using it. We evaluated Psat of five samples by applying DC power and confirmed the overall trend between Psat and the inverse leg length. We evaluated tau0 of the samples by applying DC plus AC power, and the measured value was reasonable in consideration of the expected values of other TES parameters. This evaluation method enables us to verify whether a TES has been fabricated with the designed values and to provide feedback for fabrication for future CMB observations.

We present a catalog of X-ray Detected Be Stars (XDBS) with 161 Be stars from the Be Star Spectra (BeSS) database having X-ray counterparts in the Chandra Source Catalog v2.0, XMM-Newton 4XMM-DR11 Catalog, or Swift 2SXPS Catalog. The multi-wavelength catalog includes accurate optical positions, X-ray properties (fluxes, photon indices and hardness ratios), optical, near-infrared and infrared photometry, source classifications (when available), and other properties including proper motions, effective temperatures, X-ray to optical flux ratios as well. We also provide a convenient graphical user interface which allows for easy visualization of the catalog content.

In analogy with GRB 190114C, we here analyze the ultrarelativistic prompt emission (UPE) of GRB 180720B observed in the rest-frame time interval $t_{\rm rf}=4.84$--$10.89$~s by Fermi-GBM. We reveal the UPE hierarchical structure from the time-resolved spectral analysis performed in time sub-intervals: the spectrum in each shorter time interval is always fitted by a composite blackbody plus cutoff power-law model. We explain this structure with the \textit{inner engine} of binary-driven hypernova (BdHN) model operating in a quantum electrodynamics (QED) regime. In this regime, the electric field induced by the gravitomagnetic interaction of the newborn Kerr BH with the surrounding magnetic field is overcritical, i.e., $|{\bf E}|\geq E_c$, where $E_c=m_e^2 c^3/(e\hbar)$. The overcritical field polarizes the vacuum leading to an $e^+~e^-$ pair plasma that loads baryons from the surroundings during its expansion. We calculate the dynamics of the self-acceleration of the pair-electromagnetic-baryon (PEMB) pulses to their point of transparency. We characterize the quantum vacuum polarization process in the sequences of decreasing time bins of the UPE by determining the radiation timescale, Lorentz factors, and transparency radius of the PEMB pulses. We also estimate the strength of the surrounding magnetic field $\sim 10^{14}$ G, and obtain a lower limit to the BH mass, $M=2.4~M_\odot$, and correspondingly an upper limit to the spin, $\alpha = 0.6$, from the conditions that the UPE is powered by the Kerr BH extractable energy and its mass is bound from below by the NS critical mass.

We analyse the observational signatures of galactic magnetic fields that are self-consistently generated in magnetohydrodynamic simulations of the interstellar medium through turbulence driven by supernova (SN) explosions and differential rotation. In particular, we study the time evolution of the Faraday rotation measure (RM), synchrotron radiation, and Stokes parameters by characterising the typical structures formed in the plane of observation. We do this by defining two distinct models for both thermal and cosmic ray (CR) electron distributions. Our results indicate that the maps of RM have structures which are sheared and rendered anisotropically by differential rotation and that they depend on the choice of thermal electrons model as well as the SN rate. Synchrotron maps are qualitatively similar to the maps of the mean magnetic field along the line of sight and structures are only marginally affected by the CR model. Stokes parameters and related quantities, such as the degree of linear polarisation, are highly dependent on both frequency and resolution of the observation.

Nowadays we know that the origin of the Cosmic X-ray Background (CXB) is due to the integrated emission of nearby active galactic nuclei. Thus, to obtain a precise estimate of the contribution of different source classes to the CXB it is crucial to fully characterize the hard X-ray sky. We present a multifrequency analysis of all sources listed in the 3d release of the Palermo Swift-BAT hard X-ray catalog (3PBC) with the goal of (i) identifying and classifying the largest number of sources adopting multifrequency criteria, with particular emphasis on extragalactic populations and (ii) extracting Seyfert galaxies to present here the 2nd release of the Turin-SyCAT catalog. We outline a classification scheme based on radio, infrared and optical criteria that allows us to distinguish between unidentified and unclassified hard X-ray sources, and classify and classify the remaining ones discriminating between Galactic and extragalactic classes. Our revised version of the 3PBC lists 1176 classified, 820 extragalactic, and 356 Galactic, 199 unclassified, and 218 unidentified sources. According to our analysis, the hard X-ray sky is mainly populated by Seyfert galaxies and blazars. For the Seyfert galaxies, we present the 2nd release of the Turin-SyCAT including a total of 633 Seyfert galaxies, with 282 new sources, which is an 80% increment from the previous release. We confirm, that there is no clear difference between the flux distribution of the infrared-to-hard X-ray flux ratio of Seyfert galaxies Type 1 and Type 2. However, we confirm a significant trend between the mid-IR flux and hard X-ray flux.

One of the models explaining the high luminosity of pulsing ultra-luminous X-ray sources (pULXs) was suggested by Mushtukov et al. (2015). They showed that the accretion columns on the surfaces of highly magnetized neutron stars can be very luminous due to opacity reduction in the high magnetic field. However, a strong magnetic field leads also to amplification of the electron-positron pairs creation. Therefore, increasing of the electron and positron number densities compensates the cross-section reduction, and the electron scattering opacity does not decrease with the magnetic field magnification. As a result, the maximum possible luminosity of the accretion column does not increase with the magnetic field. It ranges between 10$^{40} - 10^{41}$ erg s$^{-1}$ depending only slightly on the magnetic field strength.

For more than 20 years, the community has heavily relied on CORSIKA for the simulation of extensive air showers, their Cherenkov light emission and their radio signals. While tremendously successful, the Fortran-based monolithic design of CORSIKA up to version 7 limits adaptation to new experimental needs, for example in complex scenarios where showers transition from air into dense media, and to new computing paradigms such as the use of multi-core and GPU parallelization. With CORSIKA 8, we have reimplemented the core functionality of CORSIKA in a modern, modular, C++-based simulation framework, and successfully validated it against CORSIKA 7. Here, we discuss the philosophy of CORSIKA 8, showcase some example applications, and present the current state of implementation as well as the plans for the future.

The HAWC observatory is an air-shower detector, which is designed to study both astrophysical gamma-rays in the TeV region and galactic cosmic rays in the energy interval from 1 TeV to 1 PeV. This energy regime is quite interesting for cosmic ray research, since indirect observations overlap with direct measurements, which offers the opportunity for cross calibration and studies of experimental systematic errors in both techniques. One quantity that could help for this purpose is the all-particle energy spectrum of cosmic rays. In this work, we present an update of HAWC measurements on the total cosmic-ray energy spectrum between 10 TeV and 1 PeV. The spectrum was obtained from an unfolding analysis of almost two years of HAWC's data, which was collected from January, 2018 to December, 2019. For the energy estimation, we employed the high-energy hadronic interaction model QGSJET-II-04. The results show the presence of a knee-like structure around tens of TeVs, which was previously reported by the HAWC collaboration in 2017.

We present a reconstruction algorithm for extensive air showers with zenith angles between 65$^\circ$ and 85$^\circ$ measured with radio antennas in the 30-80 MHz band. Our algorithm is based on a signal model derived from CoREAS simulations which explicitly takes into account the asymmetries introduced by the superposition of charge-excess and geomagnetic radiation as well as by early-late effects. We exploit correlations among fit parameters to reduce the dimensionality and thus ensure stability of the fit procedure. Our approach reaches a reconstruction efficiency near 100% with an intrinsic resolution for the reconstruction of the electromagnetic energy of well below 5\%. It can be employed in upcoming large-scale radio detection arrays using the 30-80 MHz band, in particular the AugerPrime Radio detector of the Pierre Auger Observatory, and can likely be adapted to experiments such as GRAND operating at higher frequencies.

In a previous work, we investigated the evolution of the flow field around sunspots during sunspot decay and compared it with the flow field of supergranular cells. The decay of a sunspot proceeds as it interacts with its surroundings. This is manifested by the changes observed in the flow field surrounding the decaying spot. We now investigate in detail the evolution of the flow field in the direct periphery of the sunspots of the same sample and aim to provide a complete picture of the role of large-scale flows present in sunspot cells. We analyse the horizontal velocity profiles of sunspots obtained from observations by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO). We follow their evolution across the solar disc from their stable phase to their decay and their final disappearance. We find two different scenarios for the evolution of the flow region surrounding a spot in the final stage of its decay: (i) either the flow cell implodes and disappears under the action of the surrounding supergranules or (ii) it outlives the spot. In the later case, an inwards flow towards the remaining naked spot develops in the vicinity closest to the spot followed by an outflow further out. These findings provide observational evidence to theoretical predictions by realistic magnetohydrodynamic (MHD) sunspot and moat region simulations. The Evershed flow and the moat flow, both connected to the presence of fully fledged sunspots in a spot cell, vanish when penumbrae decay. Moat flows decline into supergranular flows. The final fate of a spot cell depends on its interaction with the surrounding supergranular cells. In the case of non-imploding spot cells, the remaining naked spot develops a converging inflow driven by radiative cooling and a geometrical alignment of granules in its periphery which is similar to that observed in pores.

DAMPE (Dark Matter Particle Explorer) is a satellite-born experiment launched in 2015 in a sun-synchronous orbit at 500 km altitude, and it has been taking data in stable conditions ever since. Its main goals include the spectral measurements of cosmic electrons and positrons, protons, nuclei and gamma rays, up to very high energies. The detector's main features include the 32 radiation lengths deep calorimeter and large geometric acceptance, making DAMPE one of the most powerful space instruments in operation, covering with high statistics and small systematics the high energy frontier up to several hundreds TeV. The results of different species spectral measurements will be shown and discussed.

Evidence of disequilibrium chemistry due to vertical mixing in the atmospheres of many T and Y-dwarfs has been inferred due to enhanced mixing ratios of CO and reduced NH$_3$. Atmospheric models of planets and brown dwarfs typically parameterize this vertical mixing phenomenon with the vertical eddy diffusion coefficient, $K_{\rm zz}$. While $K_{\rm zz}$ can perhaps be approximated in the convective regions in the atmosphere with mixing length theory, in radiative regions the strength of vertical mixing is uncertain by many orders of magnitude. With a new grid of self-consistent 1D model atmospheres from $T_{\rm eff}$ of 400 - 1000 K, computed with a new radiative-convective equilibrium python code PICASO 3.0, we aim to assess how molecular abundances and corresponding spectra can be used as a probe of depth-dependent $K_{\rm zz}$. At a given surface gravity, we find non-monotonic behavior in the CO abundance as a function of $T_{\rm eff}$, as chemical abundances are sometimes quenched in either of two potential atmospheric convective zones, or quenched in either of two possible radiative zones. The temperature structure and chemical quenching behavior also changes with gravity. We compare our models with available near-infrared and M-band spectroscopy of several T and Y-dwarfs and assess their atmospheric vertical mixing profiles. We also compare to color-magnitude diagrams and make predictions for JWST spectra. This work yields new constraints, and points the way to significant future gains, in determining $K_{\rm zz}$, a fundamental atmospheric parameter in substellar atmospheres, with significant implications for chemistry and cloud modeling.

A fireball of radiation plasma created near the surface of a neutron star (NS) expands under its own pressure along magnetic field lines, and produces photon emission and relativistic matter outflow. We comprehensively classify the expanding fireball evolution into five cases and obtain the photospheric luminosity and the kinetic energy of the outflow, taking into account key processes; lateral diffusion of photons escaping from a magnetic flux tube, effects of strong magnetic field, baryon loading from the NS surface, and radiative acceleration via cyclotron resonant scattering, some of which have not been considered in the context of gamma-ray bursts. Applying our model to magnetar bursts with fast radio bursts (FRBs), in particular the X-ray short bursts from SGR 1935+2154 associated with the Galactic FRB 20200428A, we show that the burst radiation can accelerate the outflow to high Lorentz factor with sufficient energy to power FRBs.

We studied an X1.6 solar flare produced by NOAA AR 12602 on 2014 October 22. The entirety of this event was covered by RHESSI, IRIS, and Hinode/EIS, allowing analysis of the chromospheric response to a nonthermal electron driver. We derived the energy contained in nonthermal electrons via RHESSI spectral fitting, and linked the time-dependent parameters of this call to the response in Doppler velocity, density, and nonthermal width across a broad temperature range. The total energy injected was $4.8\times10^{30}$ erg, and lasted $352$ seconds. This energy drove explosive chromospheric evaporation, with a delineation in both Doppler and nonthermal velocities at the flow reversal temperature, between 1.35--1.82 MK. The time of peak electron injection (14:06 UT) corresponded to the time of highest velocities. At this time, we found 200 km s$^{-1}$ blueshifts in the core of Fe XXIV, which is typically assumed to be at rest. Shortly before this time, the nonthermal electron population had the shallowest spectral index ($\approx$ 6), corresponding to the peak nonthermal velocity in Si IV and Fe XXI. Nonthermal velocities in Fe XIV, formed near the flow reversal temperature were low, and not correlated with density or Doppler velocity. Nonthermal velocities in ions with similar temperatures were observed to increase and correlate with Doppler velocities, implying unresolved flows surrounding the flow reversal point. This study provides a comprehensive, time-resolved set of chromospheric diagnostics for a large X-class flare, along with a time-resolved energy injection profile, ideal for further modeling studies.

We study the signal of anisotropy in AGNs/quasars of CatWISE2020 catalogue using different observables. It has been reported earlier that this data shows a strong signal of dipole anisotropy in the source number counts. We test this claim using two independent data analysis procedures and find our number count dipole consistent with the earlier results. In addition to number counts, we test for the anisotropy signal in two other observables -- mean spectral index $\bar{\alpha}$ and mean flux density $\bar{B}$. We find a dipole signal of considerable strength both in the mean spectral index and the mean flux density. The dipole in mean flux density points towards the galactic center and becomes very weak after imposing a flux cut to remove sources with flux greater than 1 mJy. This can be attributed to the presence of some bright sources. The signal in mean spectral index, however, is relatively stable as a function of both flux and galactic cuts. The dipole in this observable points roughly opposite to the galactic center and hence most likely arises due to galactic bias. Hence, the signal in both the mean spectral index and mean flux density appears to be consistent with isotropy.

Extended sources of the stochastic gravitational backgrounds have been conventionally searched on the spherical harmonics bases. The analysis during the previous observing runs by the ground-based gravitational wave detectors, such LIGO and Virgo, have yielded the constraints on the angular power spectrum $C_\ell$, yet it lacks the capability of estimating model parameters. In this paper, we introduce an alternative Bayesian formalism to search for such stochastic signals with a particular distribution of anisotropies on the sky. This approach provides a Bayesian posterior of model parameters and also enables selection tests among different signal models. While the conventional analysis fixes the highest angular scale \textit{a priori}, here we show a more systematic and quantitative way to determine the cut-off scale based on a Bayes factor, which depends on the amplitude and the angular scale of observed signals. Also, we analyze the third observing runs of LIGO and Virgo for the population of milli-second pulsars and obtain the 95 % constrains of the signal amplitude, $\epsilon < 2.7\times 10^{-8}$.

In this paper we explore the existing tensions in the local cosmological expansion rate, $H_0$, and amplitude of the clustering of large-scale structure at $8\, h^{-1}\mathrm{Mpc}$, $\sigma_8$, as well as models that claim to alleviate these tensions. We consider seven models: evolving dark energy ($w$CDM), extra radiation ($N_\mathrm{eff}$), massive neutrinos, curvature, primordial magnetic fields (PMF), self-interacting neutrino models, and early dark energy (EDE). We test these models against three data sets that span the full range of measurable cosmological epochs, have significant precision, and are well-tested against systematic effects: the Planck 2018 cosmic microwave background data, the Sloan Digital Sky Survey baryon acoustic oscillation scale measurements, and the Pantheon catalog of Type Ia supernovae. We use the recent SH0ES $H_0$ measurement and several measures of $\sigma_8$ (and its related parameter $S_8=\sigma_8\sqrt{\Omega_\mathrm{m}/0.3}$). We find that four models are above the "strong" threshold in Bayesian model selection, $w$CDM, $N_\mathrm{eff}$, PMF, and EDE. However, only EDE also relieves the $H_0$ tension in the full data sets to below 2$\sigma$. Contrarily, no model alleviates the $S_8/\sigma_8$ tension in the full data set, nor does better than $\Lambda$CDM in the combined case of both $H_0$ and $S_8/\sigma_8$ tensions.

We propose and apply a new parameterization of the modified chiral effective model to study rotating neutron stars with hyperon cores in the framework of the relativistic mean-field theory. The inclusion of mesonic cross couplings in the model has improved the density content of the symmetry energy slope parameters, which are in agreement with the findings from recent terrestrial experiments. The bulk viscosity of the hyperonic medium is analyzed to investigate its role in the suppression of gravitationally driven $r$-modes. The hyperonic bulk viscosity coefficient caused by non-leptonic weak interactions and the corresponding damping timescales are calculated and the $r$-mode instability windows are obtained. The present model predicts a significant reduction of the unstable region due to a more effective damping of oscillations. We find that from $\sim 10^8$ K to $\sim 10^{9}$ K, hyperonic bulk viscosity completely suppresses the $r$-modes leading to a stable region between the instability windows. Our analysis indicates that the instability can reduce the angular velocity of the star up to $\sim$0.3~$\Omega_K$, where $\Omega_K$ is the Kepler frequency of the star.

We examine effective field theories (EFTs) with a continuum sector in the presence of gravity. We first explain, via arguments based on central charge and species scale, that an EFT with a free continuum cannot consistently couple to standard (i.e. 4D Einstein) gravity. It follows that EFTs with a free, or nearly-free, continuum must either have a finite number of degrees of freedom or nonstandard gravity. We demonstrate the latter through holographically-defined continuum models, focusing on a class of 5D dilaton-graviton systems giving rise to a gapped continuum (i.e. the linear dilaton background). In the simplest version of the model we find an $R^{-2}$ deviation from the Newtonian potential. At finite temperature, we find an energy density with $a^{-5}$ scaling law (i.e. $w=\frac{2}{3}$) in the brane Friedmann equation, induced by the horizon in the bulk. We also present a slightly more evolved model for which these exotic deviations transition into those from pure AdS. Brane cosmology in dilaton-gravity backgrounds could be explored along these lines.

Massive spinning particles acquire helicity-dependent chemical potentials during the inflation from axion-type couplings. Such spinning fields can mediate sizable inflaton correlators which we call the helical inflation correlators. Helical inflaton correlators are approximately scale invariant, dS boost breaking, parity-violating, and are promising observables of cosmological collider physics. In this work, we present complete and analytical results for 4-point helical inflation correlators with tree-level exchanges of massive spinning particles, including both the smooth background and the oscillatory signals. We compute the bulk Schwinger-Keldysh integrals in two independent ways, including the partial Mellin-Barnes representation and solving bootstrap equations. We also present new closed-form analytical results for 3-point functions with massive scalar or helical spinning exchanges. The analytical results allow us to concretely and efficiently explore the phenomenological consequences of helicity-dependent chemical potentials. In particular, we show that the chemical potential can exponentially enhance oscillatory signals of both local and nonlocal types, but only affects the background in a rather mild way. Our results extend the de Sitter bootstrap program to include nonperturbative breaking of de Sitter boosts. Our results also explicitly verify the recently proposed cutting rule for cosmological collider signals.

Large classes of standard single-field slow-roll inflationary models consistent with the required number of e-folds, the current bounds on the spectral index of scalar perturbations, the tensor-to-scalar ratio, and the scale of inflation can be efficiently constructed using genetic algorithms. The setup is modular and can be easily adapted to include further phenomenological constraints. A semi-comprehensive search for sextic polynomial potentials results in roughly O(300,000) viable models for inflation. The analysis of this dataset reveals a preference for models with a tensor-to-scalar ratio in the range 0.0001 < r < 0.0004. We also consider potentials that involve cosine and exponential terms. In the last part we explore more complex methods of search relying on reinforcement learning and genetic programming. While reinforcement learning proves more difficult to use in this context, the genetic programming approach has the potential to uncover a multitude of viable inflationary models with new functional forms.

In the present work, we investigate the neutral-current neutrino-nucleon scattering in the nuclear medium using various energy-density functional (EDF) models such as the KIDS (Korea-IBS-Daegu-SKKU) and SLy4, together with the quark-meson coupling (QMC) model for the nucleon form factors at finite density. The differential cross section (DCS) and neutrino mean free path (NMFP) are computed numerically, considering the density-dependent nucleon form factors (DDFF) and neutrino structural properties such as the neutrino magnetic moment (NMM) and its electric charge radius (NCR). It turns out that the DDFF decreases the scattering cross-section, while the NCR increases it considerably. The effect of the NMM turns out to be almost negligible. We also observe that the value of the neutron effective mass is of importance in the neutron-star cooling process, indicating that for the neutron effective mass larger than the mass in free space, the neutrino can interact with matter at densities $\rho \gtrsim 1.5 \rho_0$ in the neutron star with radius 13 km.

We consider the angular distribution of the photon signal which could arise from velocity-dependent dark matter annihilation within the Galactic bulge. We find that, for the case of Sommerfeld-enhanced annihilation, dark matter annihilation within the bulge is dominated by slow speed particles which never leave the bulge, allowing one to find a simple analytic relationship between the dark matter profile within the Galactic bulge and the angular distribution. On the other hand, for the case $p$- or $d$-wave annihilation, we find that the small fraction of high-speed particles which can leave the bulge provide a significant, often dominant, contribution to dark matter annihilation within the bulge. For these scenarios, fully understanding dark matter annihilation deep within the Galactic bulge, and the angular distribution of the resulting photon signal, requires an understanding of the dark matter profile well outside the bulge. We consider the Galactic Center excess in light of these results, and find that an explanation of this excess in terms of $p$-wave annihilation would require the dark matter profile within the bulge to have a much steeper profile than usually considered, but with uncertainties related to the behavior of the profile outside the bulge.

Most all-sky searches for continuous gravitational waves assume the source to be isolated. In this paper, we allow for an unknown companion object in a long-period orbit and opportunistically use previous results from an all-sky search for isolated sources to constrain the continuous gravitational wave amplitude over a large and unexplored range of binary orbital parameters without explicitly performing a dedicated search for binary systems. The resulting limits are significantly more constraining than any existing upper limit for unknown binary systems, albeit the latter apply to different orbital parameter ranges.

This brief paper generalizes dark matter (DM) constraints recently derived in [A.E.Egorov, Phys. Rev. D 106, 023023 (2022)] by radio observations of M31 to all possible annihilation channels (except the cases of neutrinos and photons as primary annihilation products). All the methodology here exactly repeats that in the mentioned paper, where only two representative and popular annihilation channels - $\chi\chi \rightarrow b\overline{b}$ and $\chi\chi \rightarrow \tau^+\tau^-$ - were considered. It is confirmed here, that in the case of light primary annihilation products $\chi\chi \rightarrow \tau^+\tau^-, \mu^+\mu^-,gg,c\overline{c},u\overline{u},d\overline{d},s\overline{s},e^+e^-,b\overline{b}$ (and an arbitrary combination of them) the fiducial (averaged over uncertainties) lower mass limit for the thermal weakly interacting massive particle (WIMP) is confined in the range $\approx$(40-70) GeV, which is bounded by $\chi\chi \rightarrow b\overline{b},\tau^+\tau^-$ cases. Heavier WIMP, which can annihilate to $W^+W^-,Z^0Z^0,t\overline{t},hh$; can not be probed at the level of thermal cross section, unless one assumes the optimistic cases of DM density and magnetic field distributions in M31. In conclusion, $m_x \gtrsim 40$ GeV represents the fiducial channel-independent mass limit for the thermal WIMP with the full uncertainty range estimated to be $\approx$(20-90) GeV.

Primordial Black Holes (PBHs) in the mass range $\sim 10^{17}- 10^{22}$g are currently unconstrained, and can constitute the full Dark Matter (DM) density of the universe. Motivated by this, in the current work, we aim to relate the existence of PBHs in the said mass range to the production of observable Gravitational Waves (GWs) in the upcoming GW detectors. We follow a relatively model-independent approach assuming that the PBHs took birth in a radiation dominated era from enhanced primordial curvature perturbation at small scales produced by inflation. We show that the constraints from CMB and BAO data allow for the possibility of PBHs being the whole of DM density of the universe. Finally, we derive the GW spectrum induced by the enhanced curvature perturbations and show that they are detectable in the future GW detectors like eLISA, BBO and DECIGO.

Modified gravity models often contain modes that couple to normal matter and propagate with slightly less than the speed of light. High-energy cosmic rays then lose energy due to Cherenkov radiation, which constrains such models. This is also true for some MOND (Modified Newtonian Dynamics) models. However, these constraints are difficult to make precise because MOND is inherently non-linear and because the results may depend on the specific high-acceleration behavior of these models, i.e. the behavior outside the MOND regime. Recently, various hybrid MOND dark matter models were proposed, where cold dark matter (CDM) phenomenology on cosmological scales and MOND phenomenology on galactic scales share a common origin. Such models typically contain a mode that is directly coupled to matter (for MOND), but with non-relativistic sound speed (for CDM). Thus, even non-relativistic objects like stars can emit gravitational Cherenkov radiation. We calculate a lower bound on the associated energy loss. We use a controlled approximation that depends only on the MOND regime of these models. We apply our results to three concrete models: For the original superfluid dark matter model (SFDM), we rule out a part of the parameter space, including the most commonly used parameters. For two-field SFDM, we find no constraint since the matter coupling of the relevant mode is suppressed by mixing. For the recently-proposed model by Skordis and Z{\l}o\'snik, we find no constraint since the matter coupling is suppressed in non-static situations.