Papers for Thursday, Oct 21 2021

The environment of NGC1052 has recently attracted much attention due to the presence of low surface brightness galaxies (LSBGs) with claimed exotic properties, making the detection of new objects in this region of high interest. We have used public deep photometric data from DECaLS to carry out a comprehensive search for LSBGs over a wide region of 6x6 degrees, equivalent to 2x2 Mpc at the distance of NGC1052. We detected 42 LSBGs with r_e > 5" and mu_g(0) > 24 mag arcsec-2, of which 20 are previously undetected objects. Among all the newly detected objects, RCP32 stands out with extreme properties: r_e = 23.0" and <mu_g>_e = 28.6 mag arcsec-2. This makes RCP32 one of the lowest LSBGs ever detected through integrated photometry, located at just 10 arcmin from the extensively studied NGC1052-DF2. We explored the presence of globular clusters (GCs) in the LSBGs. We marginally detected a GC system in RCP32, and argue the great interest of follow-up observations on RCP32 given its extremely low baryon density. After analyzing the distribution of galaxies with available spectroscopy, we identified a large-scale structure of approximately 1 Mpc, well isolated in redshift-space, centred on NGC1052. The spatial correlation analysis between the LSBGs and this large-scale structure suggests their association. However, when exploring the distribution of effective radius we found an overpopulation of large LSBGs (r_e > 15") located close to the line of sight of NGC1052. We argue that this is suggestive of a substructure with similar radial velocity in sight projection, but at a closer distance, to which some of these apparently larger LSBGs could be associated, however possible effects derived from tidal interactions are worth further study. Our work expands the catalogue of LSBGs with new interesting objects and provides a detailed environmental context for the study of LSBGs in this region.

We introduce the Galaxy Replacement Technique (GRT) that allows us to model tidal stripping of galaxies with very high-mass (m$_{\rm{star}}=5.4\times10^4$~M$_\odot$/h) and high-spatial resolution (10 pc/h), in a fully cosmological context, using an efficient and fast technique. The technique works by replacing multiple low-resolution DM halos in the base cosmological simulation with high-resolution models, including a DM halo and stellar disk. We apply the method to follow the hierarchical build-up of a cluster since redshift $\sim8$ to now, through the hierarchical accretion of galaxies, individually or in substructures such as galaxy groups. We find we can successfully reproduce the observed total stellar masses of observed clusters since redshift $\sim$1. The high resolution allows us to accurately resolve the tidal stripping process and well describe the formation of ultra-low surface brightness features in the cluster ($\mu_{V}<32$ mag arcsec$^{-2}$) such as the intracluster light (ICL), shells and tidal streams. We measure the evolution of the fraction of light in the ICL and brightest cluster galaxy (BCG) using several different methods. While their broad response to the cluster mass growth history is similar, the methods show systematic differences, meaning we must be careful when comparing studies that use distinct methods. The GRT represents a powerful new tool for studying tidal effects on galaxies and exploring the formation channels of the ICL in a fully cosmological context and with large samples of simulated groups and clusters.

We investigate the scaling relations between X-ray observables of the clusters detected in the eFEDS field using Spectrum-Roentgen-Gamma/eROSITA observations taking into account the selection effects and the distributions of observables with cosmic time. We extract X-ray observables (Lx, Lbol, T, Mgas, Yx) within R500 for the sample of 542 clusters in the eFEDS field. By applying detection and extent likelihoods, we construct a subsample of 265 clusters with a contamination level of <10% (including AGNs and spurious fluctuations) to be utilized in the scaling relation analysis. The selection function based on the state-of-the-art simulations of the eROSITA sky is fully accounted for in our work. We provide the X-ray observables in the core-included <R500 and core-excised 0.15*R500-R500 apertures for 542 galaxy clusters and groups detected in the eFEDS field. Additionally, we present our best-fit results for the normalization, slope, redshift evolution and intrinsic scatter parameters of the X-ray scaling relations between Lx-T, Lx-Mgas, Lx-Yx, Lbol-T, Lbol-Mgas, Lbol-Yx and Mgas-T. We find that the best-fit slopes significantly deviate from the self-similar model at a >3sigma confidence level however, our results are in good agreement with the simulations including non-gravitational physics and the recent results that take into account selection effects. Strong deviations we find from the self-similar scenario indicate that the non-gravitational effects play an important role in shaping the observed physical state of clusters. This work extends the scaling relations to low mass, low luminosity galaxy cluster and group regime using eFEDS observations, demonstrating eROSITA's ability to measure ICM emission out to R500 with survey-depth exposures and constrain the scaling relations in a wide mass-luminosity-redshift range.

We report the discovery of five new doubly-imaged lensed quasars from the first 2 500 square degrees of the ongoing Canada-France Imaging Survey (CFIS), which is a component of the Ultraviolet Near Infrared Optical Northern Survey (UNIONS), selected from initial catalogues of either $Gaia$ pairs or MILLIQUAS quasars. We take advantage of the deep, 0.6 arcsec median-seeing $r$-band imaging of CFIS to confirm the presence of multiple point sources with similar colour of $u-r$, via convolution of the Laplacian of the point spread function. Requiring similar-colour point sources with flux ratios less than 2.5 mag in $r$-band, reduces the number of candidates from 256 314 to 7 815. After visual inspection we obtain 30 high-grade candidates, and prioritise spectroscopic follow-up for those showing signs of a lensing galaxy upon subtraction of the point sources. We obtain long-slit spectra for 18 candidates with ALFOSC on the 2.56-m Nordic Optical Telescope (NOT), confirming five new doubly lensed quasars with $1.21<z<3.36$ and angular separations from 0.8 to 2.5 arcsec. One additional system is a probable lensed quasar based on the CFIS imaging and existing SDSS spectrum. We further classify six objects as nearly identical quasars -- still possible lenses but without the detection of a lensing galaxy. Given our recovery rate ($83\%$) of existing optically bright lenses within the CFIS footprint, we expect that a similar strategy, coupled with $u-r$ colour-selection from CFIS alone, will provide an efficient and complete discovery of small-separation lensed quasars of source redshifts below $z=2.7$ within the CFIS $r$-band magnitude limit of 24.1 mag.

The environment within dark matter haloes can quench the star formation of galaxies. However, environmental effects beyond the virial radius of haloes ($\gtrsim$ 1 Mpc) are less evident. An example is the debated correlation between colour or star formation in central galaxies and neighbour galaxies in adjacent haloes at large separations of several Mpc, referred to as two-halo galactic conformity. We use two galaxy catalogues generated from different versions of the semi-analytic model SAG applied to the MDPL2 cosmological simulation and the IllustrisTNG300 cosmological hydrodynamical simulation to study the two-halo conformity by measuring the quenched fraction of neighbouring galaxies as a function of the real-space distance from central galaxies. We find that low-mass central galaxies in the vicinity of massive systems ($M_{\rm 200c}$ $\geq$ 10$^{13}$ $h^{-1}~\rm M_{\odot}$) out to 5 $h^{-1}$ Mpc are preferentially quenched compared to other central galaxies at fixed stellar mass $M_{\star}$ or fixed host halo mass $M_{\rm 200c}$ at $z$ ~ 0. In all the galaxies catalogues is consistent that the low-mass ($M_{\star} < 10^{10}$ $h^{-1}~\rm M_{\odot}$ or $M_{\rm 200c} < 10^{11.8}$ $h^{-1}~\rm M_{\odot}$) central galaxies in the vicinity of groups and clusters of galaxies mostly produce the two-halo galactic conformity. On average, the quenched low-mass central galaxies are typically much closer to massive haloes than star-forming central galaxies of the same mass (by a factor of ~5). Our results support that the environmental influence of massive haloes can extend beyond the virial radius and affect nearby low-mass central galaxies.

Context. Strong lensing mass measurements require the knowledge of the redshift of both the lens and the source galaxy. Traditionally, spectroscopic redshifts are used for this purpose. Upcoming surveys, however, will lead to the discovery of $\sim10^5$ strong lenses, and it will be very difficult to obtain spectroscopic redshifts for most of them. Photometric redshift measurements will also be very challenging, due to the blending between lens and source light. Aims. The goal of this work is to demonstrate how to carry out an inference of the structural properties of the galaxy population from the analysis of a set of strong lenses with no individual source redshift measurements, and to assess the loss in precision compared to the case in which spectroscopic redshifts are available. Methods. Building on the formalism introduced in Paper III, I developed a method that allows to carry out a statistical strong lensing inference while marginalising over the source redshifts. This method, which relies on the knowledge of the properties of the unlensed background source population and of the selection function of the survey, generalises an approach known as photogeometric redshift, originally introduced by the Strong Lensing Legacy Survey collaboration. I tested the method on simulated data consisting of a subset of 137 strong lenses that is complete above a cut in observational space. Results. The method recovers the properties of the galaxy population, with a precision that is comparable to that attainable in the case in which individual source redshifts are known. Conclusions. The photogeometric redshift method is a viable approach for the analysis of large sets of strong lenses, provided that the background source population properties and lens selection function are well-known.

Galactic nuclei (GNs) are dense stellar environments abundant in gravitational-wave (GW) sources for LIGO, VIRGO, and KAGRA. The GWs may be generated by stellar-mass black hole (BH) or neutron star mergers following gravitational bremsstrahlung, dynamical scattering encounters, Kozai-Lidov type oscillations driven by the central supermassive black hole (SMBH), or gas-assisted mergers if present. In this paper, we examine a smoking gun signature to identify sources in GNs: the GWs scattered by the central SMBH. This produces a secondary signal, an astrophysical GW echo, which has a very similar time-frequency evolution as the primary signal but arrives after a time delay. We determine the amplitude and time-delay distribution of the GW echo as a function of source distance from the SMBH. Between $10\%-90\%$ of the detectable echoes arrive within $(1-100)M_6\,\mathrm{sec}$ after the primary GW for sources between $10-10^4\,r_{\rm S}$, where $r_{\rm S}=2GM/c^2$, $M$ is the observer-frame SMBH mass, and $M_6=M/(10^6\,M_{\odot})$. The echo arrival times are systematically longer for high signal-to-noise ratio (SNR) primary GWs, where the GW echo rays are scattered at large deflection angles. In particular, $10\%-90\%$ of the distribution is shifted to $(5-1800)M_6\,\mathrm{sec}$ for sources, where the lower limit of echo detection is $0.02$ of the primary signal amplitude. We find that $5\%-30\%$ ($1\%-7\%$) of GW sources have an echo amplitude larger than $0.2-0.05$ times the amplitude of primary signal if the source distance from the SMBH is $r=50\,r_{\rm S}$ ($200\,r_{\rm S}$). Non-detections can rule out that a GW source is near an SMBH.

The eROSITA Final Equatorial-Depth Survey (eFEDS) executed during the performance verification phase of the Spectrum- Roentgen-Gamma (SRG)/eROSITA telescope was completed in November 2019. One of the science goals of this survey is to demonstrate the ability of eROSITA to detect samples of clusters and groups at the final depth of the eROSITA All-Sky Survey. Due to the sizeable point-spread function of eROSITA, high redshift clusters of galaxies or compact nearby groups hosting bright active galactic nuclei (AGN) can be misclassified as point sources by the source detection algorithms. A total of 357 galaxy clusters and groups were identified based on their red sequence in the eFEDS point source catalog. In this work, we examine the multi-wavelength properties of these clusters and groups to understand the potential biases in our selection process and the completeness of the extent-selected sample. We find that the majority of the clusters and groups in the point source sample are indeed under-luminous and compact compared to the extent-selected sample. Their faint X-ray emission, well under the flux limit of the extent-selected eFEDS clusters, and their compact X-ray emission are likely to be the main reason behind this misclassification. In the sample, we confirm that 11% of the sources host AGN in their brightest cluster galaxies (BCGs). By studying their X-ray, optical, and infrared, and radio properties, we establish a method to identify the most extraordinary clusters and groups that host AGNs in their BCGs. We test this method through the SDSS optical spectroscopy successfully on the current point source catalog and find four low mass clusters and groups with active radio-loud AGNs. The method we developed will be used to identify high redshift clusters and AGN dominated groups and low mass clusters misclassified in the future eROSITA All-Sky Survey point source catalogs.

In recent years, many Type IIn supernovae have been found to share striking similarities with the peculiar SN 2009ip, whose true nature is still under debate. Here, we present 10 years of observations of SN 2011fh, an interacting transient with spectroscopic and photometric similarities to SN 2009ip. SN 2011fh had a M$_r \sim -16$ mag brightening event, followed by a brighter M$_r \sim -18$ mag explosion in August 2011. The spectra of SN 2011fh are dominated by narrow to intermediate Balmer emission lines throughout its evolution, with P Cygni profiles indicating fast-moving material at $\sim 6400 \ \textrm{km s}^{-1}$. HST/WFC3 observations from October of 2016 revealed a bright source with $\textrm{M}_{F814W} \approx -13.3$ mag, indicating that the star might be going through an eruptive phase five years after the bright explosion of 2011. Using HST photometry of the stellar cluster around SN 2011fh, we estimated an age of $\sim 4.5$ Myr for the progenitor, which implies a stellar mass of $\sim 60$ M$_\odot$. VLT/MUSE observations show that the source is located in a star-forming region with an oxygen abundance of $12 + \textrm{log}_{10}(\textrm{O/H}) \approx 8.39$. We also show that the progenitor of SN 2011fh exceeded the classical Eddington limit by a large factor in the months preceding the bright explosion of 2011, suggesting strong super-Eddington winds as a possible mechanism for the observed mass-loss. These findings favor an energetic explosion in a young and massive star, possibly a luminous blue variable.

The Universe is not just cold dark matter and dark energy, it also contains baryons, radiation and neutrinos. The presence of these components, beyond the pressure-less cold dark matter and the quasi-uniform dark energy ones, imply that the single clock assumption from inflation is no longer preserved. Here we quantify this effect and show that the single-clock symmetry is ensured only on scales where baryonic effects, neutrinos effects, or sound speed are zero. These scales depend on the cosmic epoch and the Universe composition. Hence for all use and purposes of interpreting state-of-the-art and possibly forthcoming surveys, in the accessible scales, single clock symmetry cannot be said to be satisfied. Breaking the single-clock symmetry has key consequences for the study of non-Gaussian features generated by pure single-field inflation which arise from non-linearities in the metric yielding non-Gaussianities of the local type: the $n_{s}-1$ and the relativistic $-5/3$ term.

We present a multi-wavelength characterisation of the planetary nebula (PN) NGC 6905 and its [Wolf-Rayet]-type ([WR]) central star (CSPN) HD 193949. Our Nordic Optical Telescope (NOT) Alhambra Faint Object Spectrograph and Camera (ALFOSC) spectra and images unveil in unprecedented detail the high-ionization structure of NGC 6905. The high-quality spectra of HD 193949 allowed us to detect more than 20 WR features including the characteristic O-bump, blue bump and red bump, which suggests a spectral type no later than a [WO2]-subtype. Moreover we detect the Ne VII and Ne VIII broad emission lines, rendering HD 193949 yet another CSPN with $T_\mathrm{eff}\lesssim150$ kK exhibiting such stellar emission lines. We studied the physical properties ($T_\mathrm{e}$ and $n_\mathrm{e}$) and chemical abundances of different regions within NGC 6905 including its low-ionization clumps; abundances are found to be homogeneous. We used the PoWR stellar atmosphere code to model the spectrum of HD 193949, which is afterwards used in a photoionization model performed with Cloudy that reproduces the nebular and dust properties for a total mass in the 0.31-0.47 M$_{\odot}$ range and a mass of C-rich dust of $\sim$2 $\times10^{-3}$ M$_{\odot}$. Adopting a current stellar mass of 0.6 M$_{\odot}$, our model suggests an initial mass $\sim$1 M$_\odot$ for HD 193949, consistent with the observations.

An only early or only late time alteration to $\Lambda$CDM has been inadequate at resolving both the $H_0$ and $S_8$ tensions simultaneously; however, a combination of early and late time alterations to $\Lambda$CDM can provide a solution to both tensions. As an illustration, we examine a combined Early Dark Energy - Decaying Dark Matter model. While early dark energy has the ability to resolve the $H_0$ tension, it leads to a discrepancy in $S_8$ measurements. We show that the addition of decaying dark matter helps resolve the $S_8$ discrepancy that would otherwise be enhanced in an early dark energy model, while the latter is able to relieve the $H_0$ disagreement to within the 95th percentile interval. Our results show a preference for the combined model over $\Lambda$CDM with $\Delta \rm{AIC} = -6.72$, hinting that both early and late universe modifications may be necessary to address the cosmological tensions.

TOI-216 is a pair of close-in planets with orbits deep in the 2:1 mean motion resonance. The inner, Neptune-class planet (TOI-216b) is near 0.12 au (orbital period $P_{\rm b} \simeq 17$ d) and has a substantial orbital eccentricity ($e_{\rm b} \simeq 0.16$), and large libration amplitude ($A_\psi \simeq 60^\circ$) in the resonance. The outer planet (TOI-216c) is a gas giant on a nearly circular orbit. We carry out $N$-body simulations of planet migration in a protoplanetary gas disk to explain the orbital configuration of TOI-216 planets. We find that TOI-216b's migration must have been halted near its current orbital radius to allow for a convergent migration of the two planets into the resonance. For the inferred damping-to-migration timescale ratio $\tau_e/\tau_a \simeq 0.02$, overstable librations in the resonance lead to a limit cycle with $A_\psi \simeq 80^\circ$ and $e_{\rm b}<0.1$. The system could have remained in this configuration for the greater part of the protoplanetary disk lifetime. If the gas disk was removed from inside out, this would have reduced the libration amplitude to $A_\psi \simeq 60^\circ$ and boosted $e_{\rm b}$ via the resonant interaction with TOI-216c. Our results suggest a relatively fast inner disk removal ($\sim 10^5$ yr). Another means of explaining the large libration amplitude is stochastic stirring from a (turbulent) gas disk. For that to work, overstable librations would need to be suppressed, $\tau_e/\tau_a \simeq 0.05$, and very strong turbulent stirring (or some other source of large stochastic forcing) would need to overcome the damping effects of gas. Hydrodynamical simulations can be performed to test these models.

We review all existing air-fluorescence measurements of the elongation rate of extensive air showers (slope of mean EAS shower maximum (Xmax) vs log of shower energy E) above 1017 eV. We find remarkable agreement for all current and historic experiments over a 30 year period for the energy range from 1017 to 3x1018 eV. The mean elongation rate in this energy interval is near 80 gm/cm2/decade Above this energy, experiments in the Northern hemisphere are in good agreement with an average elongation rate of 48 +/- 10 gm/cm2/decade while Southern hemisphere experiments have a flatter elongation rate of 26 +/- 2 gm/cm2/decade We point out that, given the agreement at lower energies, possible systematic reasons for this difference are unlikely. Given this, the world elongation rate data alone may indicate a composition difference of UHECR in the Northern and Southern hemisphere and thus a diversity of UHECR sources in the Northern and Southern sky.

The promise by the LIGO/Virgo/Kagra (LVK) collaboration to detect black hole-neutron star (BH-NS) mergers via gravitational wave (GW) emission has recently been fulfilled with the detection of GW200105 and GW200115. Mergers of BH-NS binaries are particularly exciting for their multi-messenger potential, since the GW detection can be followed by an electromagnetic (EM) counterpart (kilonova, gamma-ray burst, afterglow) that can reveal important information on the equation of state (EOS) of NSs and the nature of the BH spin. This can happen whenever the NS does not directly plunge into the BH, but rather is tidally disrupted leaving behind debris to accrete. We carry out a statistical study of the binary stars that evolve to form a BH-NS binary and compute the rate of merger events that can be followed by an EM counterpart. We find that $\gtrsim 50\%$ of the mergers can lead to an EM counterpart only in the case BHs are born highly spinning, otherwise this fraction never exceeds about $30\%$ for stiff NS EOSs and a few percent for soft NS EOSs. However, the possibilities that BHs are born with near-maximal spins and that NS internal structure is described by a stiff EOS are disfavored by current LVK constraints. Considering that these values only represent an upper limit to observe an EM counterpart due to current observational limitations, as in brightness sensitivity and sky localization, BH-NS mergers are unlikely multi-messenger sources.

We present the discovery of 2013 VZ70, the first known horseshoe coorbital companion of Saturn. Observed by the Outer Solar System Origins Survey (OSSOS) for 4.5 years, the orbit of 2013 VZ70 is determined to high precision, revealing that it currently is in `horseshoe' libration with the planet. This coorbital motion will last at least thousands of years but ends ~10 kyr from now; 2013 VZ70 is thus another example of the already-known `transient coorbital' populations of the giant planets, with this being the first known prograde example for Saturn (temporary retrograde coorbitals are known for Jupiter and Saturn). We present a theoretical steady state model of the scattering population of trans-Neptunian origin in the giant planet region (2--34 au), including the temporary coorbital populations of the four giant planets. We expose this model to observational biases using survey simulations in order to compare the model to the real detections made by a set of well-characterized outer Solar System surveys. While the observed number of coorbitals relative to the scattering population is higher than predicted, we show that the number of observed transient coorbitals of each giant planet relative to each other is consistent with a transneptunian source.

We explore the fascinating eclipses and dynamics of the compact hierarchical triple star system KOI-126 (KIC 5897826). This system is comprised of a pair of M-dwarf stars (KOI-126 B and C) in a 1.74 day orbit which revolve around an F-star (KOI-126 A) every 34 days. Complex eclipse shapes are created as the M stars transit the F star, due to two effects: (i) the duration of the eclipse is a significant fraction of the M-star orbital period, so the prograde or retrograde motion of the M stars in their orbit lead to unusually short or long duration eclipses; (ii) due to 3-body dynamics, the M-star orbit precesses with an astonishingly quick timescale of 1.74 years for the periastron (apsidal) precession, and 2.73 years for the inclination and nodal angle precession. Using the full Kepler data set, supplemented with ground-based photometry, plus 29 radial velocity measurements that span 6 years, our photodynamical modeling yields masses of $M_{A} = 1.2713 \pm 0.0047 M_{\odot}$ (0.37%), $M_{B} = 0.23529 \pm 0.00062 M_{\odot}$ (0.26%), and $M_{C} = 0.20739 \pm 0.00055 M_{\odot}$ (0.27%) and radii of $R_{A} = 1.9984 \pm 0.0027 R_{\odot}$ (0.14%), $R_{B}= 0.25504 \pm 0.00076 R_{\odot}$ (0.3%), and $R_{C} = 0.23196 \pm 0.00069 R_{\odot}$ (0.3%). We also estimate the apsidal motion constant of the M-dwarfs, a parameter that characterizes the internal mass distribution. While not particularly precise, we measure a mean apsidal motion constant, $\overline{k_{2}}$, of $ 0.046^{+0.046}_{-0.028}$, which is approximately 2-$\sigma$ lower than the theoretical model prediction of 0.150. We explore possible causes for this discrepancy.

We searched through the entire Gaia EDR3 candidate white dwarf catalogue for stars with proper motions and positions that are consistent with the stars having escaped from the Alpha Persei cluster within the past 81~Myr, the age of the cluster. In this search we found five candidate white dwarf escapees from Alpha Persei and obtained spectra for all five. We confirm that three are massive white dwarfs sufficiently young to have originated in the cluster. All three are more massive than any white dwarf previously associated with a cluster using Gaia astrometry, and possess some of the most massive progenitors. In particular, the white dwarf Gaia~EDR3~4395978097863572, which lies within 25~pc of the cluster centre, has a mass of about 1.20 solar masses and evolved from an 8.5 solar-mass star, pushing the upper limit for white dwarf formation from a single massive star, while still leaving a substantial gap between the resulting white dwarf mass and the Chandrasekhar mass.

We have performed a molecular line search toward the flaring 6.7-GHz masers G24.33+0.13 and G359.62-0.24 using the Australia Telescope Compact Array. We present spectra of the 6.7-GHz class~II methanol and 22.2-GHz water masers toward these sources and provide comparison with other recent flaring events these sources have experienced. We also detect the fourth example of a 23.4-GHz class~I methanol maser, and the eleventh example of a 4.8-GHz formaldehyde maser toward G24.33+0.13. Alongside these results, we observe the previously detected ammonia (3,3) emission and report upper limits on the presence of various other cm-wavelength methanol, ammonia and OH transitions. Our results are consistent with the flaring of G24.33+0.13 being driven by a variable accretion rate in the host high-mass young stellar object.

The halo-mediated inverse mass cascade is a key feature of the intermediate statistically steady state for self-gravitating collisionless flow (SG-CFD). A broad spectrum of halos and halo groups are necessary to form from inverse mass cascade for long-range interaction system to maximize its entropy. The limiting velocity ($\textbf X$), speed ($\textbf Z$), and energy ($\textbf E$) distributions of collisionless particles can be obtained analytically from a maximum entropy principle. Halo mass function, i.e. the spectrum of halo mass, is a fundamental quantity for structure formation and evolution. Instead of basing mass functions on simplified spherical/elliptical collapse models, it is possible to reformulate mass function as an intrinsic distribution to maximize system entropy during the everlasting statistically steady state. Starting from halo-based description of non-equilibrium collisionless flow, distributions of particle virial dispersion ($\textbf H$), square of particle velocity ($\textbf P$), and number of halos ($\textbf J$) are proposed. Their statistical properties and connections with the limiting velocity distribution ($\textbf X$) are well studied and established. With $\textbf H$ being essentially the halo mass function, two limiting cases of $\textbf H$ distribution are analyzed for large halos ($\textbf H_\infty$) and small halos ($\textbf H_s$), respectively. For large halos, $\textbf H_\infty$ is shown to also be a maximum entropy distribution. For small halos, $\textbf H_s$ approximates the $\textbf P$ distribution and recovers the Press-Schechter mass function. The full solution of $\textbf H$ distribution depends on the limiting distribution ($\textbf X$) that maximizes system entropy and the exact models of halo velocity dispersions.

Previous observations have shown that the $\lesssim$10 au, $\gtrsim$400 K hot inner disk of the archetypal accretion outburst young stellar object, FU Ori, is dominated by viscous heating. To constrain dust properties in this region, we have performed radio observations toward this disk using the Karl G. Jansky Very Large Array (JVLA) in 2020 June-July, September, and November. We also performed complementary optical photometric monitoring observations. We found that the dust thermal emission from the hot inner disk mid-plane of FU Ori has been approximately stationary and the maximum dust grain size is $\gtrsim$1.6 mm in this region. If the hot inner disk of FU Ori which is inward of the 150-170 K water snowline is turbulent (e.g., corresponding to a Sunyaev & Shakura viscous $\alpha_{t}\gtrsim$0.1), or if the actual maximum grain size is still larger than the lower limit we presently constrain, then as suggested by the recent analytical calculations and the laboratory measurements, water-ice free dust grains may be stickier than water-ice coated dust grains in protoplanetary disks. Additionally, we find that the free-free emission and the Johnson B and V bands magnitudes of these binary stars are brightening in 2016-2020. The optical and radio variability might be related to the dynamically evolving protostellar or disk accretion activities. Our results highlight that hot inner disks of outbursting objects are important laboratories for testing models of dust grain growth. Given the active nature of such systems, to robustly diagnose the maximum dust grain sizes, it is important to carry out coordinated multi-wavelength radio observations.

We have used the Australia Telescope Compact Array (ATCA) to make new observations of the 36.2-GHz ($4_{-1}\rightarrow3_0$E) methanol transition toward NGC 4945 and NGC 253. These observations have revealed the presence of new maser components toward these galaxies, and have provided the first clear evidence for variability in extragalactic class~I methanol masers. Alongside the new observations of NGC 4945 and NGC 253, we present the results of recent 36.2-GHz methanol maser searches toward 12 galaxies, placing upper limits on the emission from the 36.2-GHz class~I transition and the 37.7-GHz ($7_{2}\rightarrow8_1$E) class~II maser line toward these sources. Flux density values for the 7-mm continuum emission toward these sources are also reported where applicable. A re-analysis of the published 36.2-GHz methanol observations of Arp 220 undertaken as part of the search revealed some issues with previous imaging procedures. The re-analysis, combined with non-detections in independent follow-up observations suggest that there is no 36.2-GHz methanol emission toward Arp 220 stronger than 3.5 mJy in a 10 km s$^{-1}$ channel (5$\sigma$ upper limit).

In this study, we demonstrate the dependence of atmospheric dust size on the near-infrared spectra of ten L dwarfs, and constrain the sizes of dust grains in each L dwarf atmosphere. In previous studies, by comparing observed and modeled spectra, it was suggested that the deviations of their spectral shapes from theoretical prediction are general characteristics. Here, we focus on the dust size in brown dwarf atmospheres to understand the observed spectra. We confirm that changing the dust size changes the temperature-pressure structure of the atmosphere, with the shape of the spectrum changing accordingly. At the wavelength at which dust is the main absorber of radiation (the dust-dominated regime), a large dust opacity combined with a medium grain size, e.g., 0.1 $\mu$m, results in a low photospheric temperature, and thus a small flux. Conversely, for the wavelength at which gas absorption is dominant (the gas-dominated regime), a large dust opacity modifies the temperature-pressure structure, resulting in a high photospheric temperature, which corresponds to large flux emissions. Taking into account the size effect, we compare the model spectral fluxes in the wavelength range 1-5 $\mu$m with the observational ones to constrain the main dust size in the atmosphere of each of the ten L dwarfs observed with AKARI and SpeX or CGS4. Ultimately, we reveal that the observed data are reproduced with higher fidelity by models based on a medium dust size of 0.1-3.0 $\mu$m for six of these L dwarfs; therefore, we suggest that such atmospheric dust sizes apply to the majority of L dwarfs.

We reanalyze the mid-infrared (5-40 $\mu$m) Spitzer spectra of 86 low-redshift ($z < 0.5$) Palomar-Green quasars to investigate the nature of polycyclic aromatic hydrocarbon (PAH) emission and its utility as a star formation rate (SFR) indicator for the host galaxies of luminous active galactic nuclei (AGNs). We decompose the spectra with our recently developed template-fitting technique to measure PAH fluxes and upper limits, which we interpret using mock spectra that simulate the effects of AGN dilution. While luminous quasars can severely dilute and affect the detectability of emission lines, PAHs are intrinsically weak in some sources that are otherwise gas-rich and vigorously forming stars, conclusively demonstrating that powerful AGNs destroy PAH molecules. Comparing PAH-based SFRs with independent SFRs derived from the mid-infrared fine-structure neon lines and the total infrared luminosity reveals that PAHs can trace star formation activity in quasars with bolometric luminosities $\lesssim 10^{46}\, \rm erg\,s^{-1}$, but increasingly underestimate the SFR for more powerful quasars, typically by $\sim 0.5$ dex. Relative to star-forming galaxies and low-luminosity AGNs, quasars have a comparable PAH 11.3 $\mu$m/7.7 $\mu$m ratio but characteristically lower ratios of 6.2 $\mu$m/7.7 $\mu$m, 8.6 $\mu$m/7.7 $\mu$m, and 11.3 $\mu$m/17.0 $\mu$m. We suggest that these trends indicate that powerful AGNs preferentially destroy small grains and enhance the PAH ionization fraction.

PSR J1813--1749 has peculiarities that make it a very interesting object of study. It is one of the most energetic and the most scattered pulsar known. It is associated with HESS J1813--178, one of the brightest and most compact TeV sources in the sky. Recently, Ho et al. used archival X-ray Chandra observations separated by more than 10 years and determined that the total proper motion of PSR J1813--1749 is $\sim66$~mas~yr$^{-1}$, corresponding to a velocity of $\sim$1900 km s$^{-1}$ for a distance of 6.2 kpc. These results would imply that this pulsar is the fastest neutron star known in the Galaxy and, by estimating the angular separation with respect to the center of the associated supernova remnant, has an age of only $\sim300$ years, making it one of the youngest pulsars known. Using archival high-angular-resolution VLA observations taken over 12 years we have estimated the radio proper motions of PSR~J1813--1748 to be much smaller: ($\mu_{\alpha}\cdot\cos(\delta),\,\mu_\delta$)=($-5.0\pm3.7,\,-13.2\pm6.7$)~mas~yr$^{-1}$, or a total proper motion of $14.8\pm5.9$~mas~yr$^{-1}$. The positions referenced against quasars make our results reliable. We conclude that PSR J1813--1749 is not a very fast moving source. Its kinematic age using the new total proper motion is $\sim1350$~years. This age is consistent within a factor of a few with the characteristic age of the pulsar and with the age estimated from the broadband spectral energy distribution of HESS J1813--178, as well as the age of the associated SNR.

We study the quiescent luminosities of accreting neutron stars by a new mechanism as neutrino heating for additional deep crustal heating, where the neutrino heating is produced by charged pions decay from the nuclear collisions on the surface of neutron star during its active accretion. For low mass neutron star($\lesssim1.4~M_{\odot}$), as the neutrino heating is little($\lesssim1$ MeV per accreted nucleon) or there would be no neutrino heating, the quiescent luminsoity will be not affected or slightly affected. While for massive neutron star ($\gtrsim2~M_{\odot}$), the quiescent luminosity will be enhanced more obviously with neutrino heating in the range 2-6 MeV per accreted nucleon. The observations on cold neutron stars such as 1H 19605+00, SAX J1808.4-3658 can be explained with neutrino heating if a fast cooling and heavy elements surface are considered. The observations on a hot neutron star such as RX J0812.4-3114 can be explained with neutrino heating if the direct Urca process is forbidden for a massive star with light elements surface, which is different from the previous work that the hot observations should be explained with small mass neutron star and the effect of superfluidity.

We introduce and validate a delensing framework for the Simons Observatory (SO), which will be used to improve constraints on inflationary gravitational waves (IGWs) by reducing the lensing noise in measurements of the $B$-modes in CMB polarization. SO will initially observe CMB by using three small aperture telescopes and one large-aperture telescope. While polarization maps from small-aperture telescopes will be used to constrain IGWs, the internal CMB lensing maps used to delens will be reconstructed from data from the large-aperture telescope. Since lensing maps obtained from the SO data will be noise-dominated on sub-degree scales, the SO lensing framework constructs a template for lensing-induced $B$-modes by combining internal CMB lensing maps with maps of the cosmic infrared background from Planck as well as galaxy density maps from the LSST survey. We construct a likelihood for constraining the tensor-to-scalar ratio $r$ that contains auto- and cross-spectra between observed $B$-modes and lensing $B$-mode template. We test our delensing analysis pipeline on map-based simulations containing survey non-idealities, but that, for this initial exploration, do not include contamination from Galactic and extragalactic foregrounds. We find that the SO survey masking and inhomogeneous and atmospheric noise have very little impact on the delensing performance, and the $r$ constraint becomes $\sigma(r)\approx 0.0015$ which is close to that obtained from the idealized forecasts in the absence of the Galactic foreground and is nearly a factor of two tighter than without delensing. We also find that uncertainties in the external large-scale structure tracers used in our multi-tracer delensing pipeline lead to bias much smaller than the $1\,\sigma$ statistical uncertainties.

We present results from a search for the HI 21-cm line in absorption towards 16 bright radio sources with the 6-antenna commissioning array of the Australian Square Kilometre Array Pathfinder (ASKAP). Our targets were selected from the 2-Jy sample, a flux-limited survey of the southern radio sky with extensive multi-wavelength follow-up. Two sources were detected in HI absorption including a new detection towards the bright FRII radio galaxy PKS 0409-75 at a redshift of z=0.674. The HI absorption line is blueshifted by ~3300 km/s compared to the optical redshift of the host galaxy of PKS 0409-75 at z=0.693. Deep optical imaging and spectroscopic follow-up with the GMOS instrument on the Gemini-South telescope reveal that the HI absorption is associated with a galaxy in front of the southern radio lobe with a stellar mass of $3.2 - 6.8 \times 10^{11}M_\odot$, a star-formation rate of $\sim 1.24 M_\odot$ yr$^{-1}$ and an estimated HI column density of $2.16\times10^{21}$ cm$^{-2}$, assuming a spin temperature of $T_{\rm spin}=500$ K and source covering factor of $C_{\rm f}=0.3$. Using polarisation measurements of PKS 0409-75 from the literature we estimate the magnetic field of the absorbing galaxy to be ~14.5$\mu$G, consistent with field strengths observed in nearby spiral galaxies, but larger than expected for an elliptical galaxy. Results from this pilot study can inform future surveys as new wide-field telescopes allow us to search for 21-cm HI absorption towards all bright radio sources as opposed to smaller targeted samples.

A new, more comprehensive model of gas-grain chemistry in hot molecular cores is presented, in which nondiffusive reaction processes on dust-grain surfaces and in ice mantles are implemented alongside traditional diffusive surface/bulk-ice chemistry. We build on our nondiffusive treatments used for chemistry in cold sources, adopting a standard collapse/warm-up physical model for hot cores. A number of other new chemical model inputs and treatments are also explored in depth, culminating in a final model that demonstrates excellent agreement with gas-phase observational abundances for many molecules, including some (e.g.~methoxymethanol) that could not be reproduced by conventional diffusive mechanisms. Observed ratios of structural isomers methyl formate, glycolaldehyde and acetic acid are well reproduced by the models. The main temperature regimes are identified in which various complex organic molecules (COMs) are formed. Nondiffusive chemistry advances the production of many COMs to much earlier times and lower temperatures than in previous model implementations. Those species may form either as by-products of simple-ice production, or via early photochemistry within the ices while external UV photons can still penetrate. Cosmic ray-induced photochemistry is less important than in past models, although it affects some species strongly over long timescales. Another production regime occurs during the high-temperature desorption of solid water, whereby radicals trapped in the ice are released onto the grain/ice surface, where they rapidly react. Several recently-proposed gas-phase COM-production mechanisms are also introduced, but they rarely dominate. New surface/ice reactions involving CH and CH$_2$ are found to contribute substantially to the formation of certain COMs.

High mass proto-stars create Ultra-Compact Hii regions (UCHII) at the stage of evolution when most of the accretion is finished but the star is still heavily embedded in molecular material. The morphologies of UCHII regions reflect the interactions of stellar winds, stellar motions, and density structure in the molecular cloud; they are complex and it has been very difficult to interpret them. We here present data obtained with TEXES on the NASA IRTF of the [NeII] emission line in the proto-cluster of young OB stars in W33 Main. The data cube has a spatial resolution of ~ 1.4 arcsec and true velocity resolution ~ 5 km/s; with A ~ 0.02Av it is relatively unaffected by extinction. We have run 3D hydrodynamic and line profile simulations, using PLUTO and RADMC-3D, of the gas structures created by multiple windy stars moving relative to each other through the ambient cloud. By iterative comparison of the data cube and the simulations, we arrive at models that reproduce the different morphology and kinematic structure of each UCHII region in W33 Main. The results indicate that each sub-source probably holds multiple exciting stars, permitting an improved view of the stellar population, and finds the stellar trajectories, which may determine the dynamical development of the proto-cluster.

In this paper, we demonstrate the importance of embedding temporal information for an accurate prediction of solar irradiance. We have used two sets of models for forecasting solar irradiance. The first one uses only time series data of solar irradiance for predicting future values. The second one uses the historical solar irradiance values, together with the corresponding timestamps. We employ data from the weather station located at Nanyang Technological University (NTU) Singapore. The solar irradiance values are recorded with a temporal resolution of $1$ minute, for a period of $1$ year. We use Multilayer Perceptron Regression (MLP) technique for forecasting solar irradiance. We obtained significant better prediction accuracy when the time stamp information is embedded in the forecasting framework, as compared to solely using historical solar irradiance values.

One probe for systematic effects in gravitational lensing surveys is the presence of so-called B-modes in the cosmic shear two-point correlation functions \xi_\pm(\vt), since lensing is expected to produce only E-mode shear. Furthermore, there exist ambiguous modes which can not uniquely be assigned to either E- or B-mode shear. We derive explicit equations for the pure-mode shear correlation functions \xi_\pm^E/B(\vt) and their ambiguous components \xi_\pm^amb(\vt), that can be derived from the measured \xi_\pm(\vt) on a finite angular interval \tmin\le\vt\le\tmax, such that the latter can be decomposed uniquely into pure-mode functions as \xi_+=\xi_+^E+\xi_+^B+\xi_+^amb and \xi_-=\xi_-^E-\xi_-^B+\xi_-^amb. The derivation is obtained by defining a new set of COSEBIs, for which explicit relations are obtained, and which yields a smaller covariance between COSEBI modes. We derive the relation between \xi_\pm^E/B/amb and the underlying E-/B-mode power spectra. The pure-mode correlation functions can provide a diagnostics of systematics in configuration space. We then apply our results to SLICS simulations and the KiDS-1000 cosmic shear data, calculate the new COSEBIs and the pure-mode correlation functions, as well as the corresponding covariances, and show that the new statistics fit equally well to the best-fitting cosmological model as the previous KiDS-1000 analysis and recovers the same level of (insignificant) B-modes. We also consider in some detail the ambiguous modes at first- and second-order level, finding some surprising results; for example, the shear field of a point mass, when cut along a line through the center, can not be ascribed uniquely to an E-mode shear and is thus ambiguous, and the shear correlation functions resulting from a random ensemble of point masses, when measured over a finite angular range, correspond to an ambiguous mode.

We review recent research on Delta Scuti stars from an observer's viewpoint. First, some signposts helping to lead the way through the Delta Scuti jungle are placed. Then, some problems in studying individual pulsators in the framework of asteroseismology are given before a view on how the study of these variables has benefited (or not) from past and present high-precision asteroseismic space missions is presented. Some possible pitfalls in the analysis of data with a large dynamical range in pulsational amplitudes are pointed out, and a strategy to optimize the outcome of asteroseismic studies of Delta Scuti stars is suggested. We continue with some views on "hybrid" pulsators and interesting individual High Amplitude Delta Scuti stars, and then take a look on Delta Scuti stars in stellar systems of several different kinds. Recent results on pre-main sequence Delta Scuti stars are discussed as are questions related to the instability strip of these variables. Finally, some remarkable new theoretical results are highlighted before, instead of a set of classical conclusions, questions to be solved in the future, are raised.

3C 84 (NGC 1275) is the radio source at the center of the Perseus Cluster and exhibits a bright radio jet. We observed the source with the Global Millimeter VLBI Array (GMVA) between 2008 and 2015, with a typical angular resolution of $\sim$50 $\mu$as. The observations revealed a consistent double nuclear structure separated by $\sim$770 gravitational radii assuming a Black Hole mass of 3.2 $\times 10^{8}$ $M_{\odot}$. The region is likely too broad and bright to be the true jet base anchored in the accretion disk or Black Hole ergosphere. A cone and parabola were fit to the stacked (time averaged) image of the nuclear region. The data did not strongly prefer either fit, but combined with a jet/counter-jet ratio analysis, an upper limit on the viewing angle to the inner jet region of $\leq$35$^{\circ}$ was found. This provides evidence for a variation of the viewing angle along the jet (and therefore a bent jet) within $\sim$0.5 parsec of the jet launching region. In the case of a conical jet, the apex is located $\sim$2400 gravitational radii upstream of the bright nuclear region and up to $\sim$600 gravitational radii upstream in the parabolic case. We found a possible correlation between the brightness temperature and relative position angle of the double nuclear components, which may indicate rotation within the jet.

Far-infrared (far-IR)/sub-mm emission linked to AGN-heated dust has been a topic of contention for many years. Results have been diverse and various views have been presented. The empirical AGN SED derived by Symeonidis et al. (2016, hereafter S16) has more far-IR/sub-mm emission than other SEDs in the literature, and thus it is contested by other works which argue that its luminosity in that part of the spectrum is overestimated. Here, I investigate this topic and the concerns raised over the S16 AGN SED. I also examine the differences between the S16 AGN SED and other commonly-used empirical AGN SEDs. My findings show that the reasons proposed by other works as to why the S16 AGN SED is not a reasonable representation of AGN emission in the far-IR/sub-mm, do not hold.

A major boost in the understanding of the universe was given by the revelation of the first coalescence event of two neutron stars (GW170817) and the observation of the same event across the entire electromagnetic spectrum. With 3rd Generation gravitational wave detectors and the new astronomical facilities, we expect many multi messenger events of the same type. We anticipate the need to analyse the data provided to us by such events, to fulfill the requirements of real-time analysis, but also in order to decipher the event in its entirety through the information emitted in the different messengers using Machine Learning. We propose a change in the paradigm in the way we will do multi-messenger astronomy, using simultaneously the complete information generated by violent phenomena in the Universe. What we propose is the application of a multimodal machine learning approach to characterize these events.

The cataclysmic explosions of massive stars as supernovae are one of the key ingredients of galaxy formation. However, their evolution is not well understood in the presence of magnetic fields or cosmic rays (CRs). We study the expansion of individual supernova remnants (SNRs) using our suite of 3D hydrodynamical (HD), magnetohydrodynamical (MHD) and CRMHD simulations generated using RAMSES. We explore multiple ambient densities, magnetic fields and fractions of supernova energy deposited as CRs ($\chi_{\rm CR}$), accounting for cosmic ray anisotropic diffusion and streaming. All our runs have comparable evolutions until the end of the Sedov-Taylor phase. However, our CRMHD simulations experience an additional CR pressure-driven snowplough phase once the CR energy dominates inside the SNR. We present a model for the final momentum deposited by supernovae that captures this new phase: $p_{\rm SNR} = 2.87\times 10^{5} (\chi_{\text{CR}} + 1)^{4.82}\left(\frac{n}{\text{cm}^{-3}}\right)^{-0.196} M_{\odot}$ km s$^{-1}$. Assuming a 10% fraction of SN energy in CRs leads to a 50% boost of the final momentum, with our model predicting even higher impacts at lower ambient densities. The anisotropic diffusion of CRs assuming an initially uniform magnetic field leads to extended gas and cosmic ray outflows escaping from the supernova poles. We also study a tangled initial configuration of the magnetic field, resulting instead in a quasi-isotropic diffusion of CRs and earlier momentum deposition. Finally, synthetic synchrotron observations of our simulations using the POLARIS code show that the local magnetic field configuration in the interstellar medium modifies the overall radio emission morphology and polarisation.

We present the unique and challenging case of a radio galaxy in Abell 3266 observed as part of the MeerKAT Galaxy Cluster Legacy Survey. It has quasi-periodic bright patches along the tail which connect to never-before-seen thin transverse extensions, which we call `ribs', reaching up to approx. 50 kpc from the central axis of the tail. At a distance of approx. 400 kpc from the host (assuming the z=0.0594 redshift of Abell 3266) we find what appears to be a triple source with its own apparent host at a photometric redshift of 0.78. Mysteriously, the part of the tail far from the host and the triple are connected by a series of thin filaments, which we call "tethers." The far tail, tethers and triple also have similar spectra and Faraday rotation measures, suggesting that there is only one -- quite complicated -- source, with a serendipitous background AGN in the triple. We look at possible causes for the "rib" and "tether" structures, and the emerging phenomena of intracluster medium filaments associated with radio galaxies.

Context: Thermal non-equilibrium (TNE) produces several observables that can be used to constrain the spatial and temporal distribution of solar coronal heating. Its manifestations include prominence formation, coronal rain, and long-period intensity pulsations in coronal loops. The recent observation of abundant periodic coronal rain associated with intensity pulsations by Auch\`ere et al. allows to unify these two phenomena as the result of TNE condensation and evaporation cycles. On the other hand, many intensity pulsation events observed by Froment et al. show little to no coronal rain formation. Aims: Our goal is to understand why some TNE cycles produce such abundant coronal rain, while others produce little to no rain. Methods: We reconstruct the geometry of the event reported by Auch\`ere et al., using images from STEREO/SECCHI/EUVI and magnetograms from SDO/HMI. We then perform 1D hydrodynamic simulations of this event, for different heating parameters and variations of the loop geometry (9000 simulations in total). We compare the resulting behaviour to simulations of TNE cycles by Froment et al. that do not produce coronal rain. Results: Our simulations show that both prominences and TNE cycles (with and without coronal rain) can form within the same magnetic structure. We show that the formation of coronal rain during TNE cycles depends on the asymmetry of the loop and of the heating. Asymmetric loops are overall less likely to produce coronal rain, regardless of the heating. In symmetric loops, coronal rain forms when the heating is also symmetric. In asymmetric loops, rain forms only when the heating compensates the asymmetry.

A GOES M1.9 flare took place in active region AR 11153 on February 9,2011. With the resolution of 200 kHz and a time cadence of 80 ms, the reverse-drifting (RS) type III bursts, intermittent sequence of type U bursts, drifting pulsation structure (DPS), and fine structures were observed by the Yunnan Observatories Solar Radio Spectrometer(YNSRS). Combined information revealed by the multi-wavelength data indicated that after the DPS which observed by YNSRS, the generation rate of type U bursts suddenly increased 5 times than before. In this event, the generation rate of type U bursts may depend on the magnetic reconnection rate. Our observations are consistent with previous numerical simulations results. After the first plasmoid produced (plasma instability occurred), the magnetic reconnection rate increased suddenly 5-8 times than before. Furthermore, after the DPS, the frequency range of turnover frequency of type U bursts is obviously broadened 3 times than before, which indicates the fluctuation amplitude of the density in the loop-top. Our observations also support the numerical simulations during the flare impulsive phase. The turbulence occurs at the top of the flare loop, the plasmoids can trap the non-thermal particles and cause the density fluctuation at the loop-top. The observations are generally consistent with the results of numerical simulations, helping us to better understand the characteristics of the whole physical process of eruption.

We re-evaluate the outer edge of orbital stability for possible exomoons orbiting the radial velocity planet discovered in the HD 23079 system. In this system, a solar-type star hosts a Jupiter-mass planet in a nearly circular orbit in the outer stellar habitable zone. The outer stability limit of exomoons is deduced using $N$-body and tidal migration simulations considering a large range of initial conditions, encompassing both prograde and retrograde orbits. In particular, we extend previous works by evaluating many values in the satellite mean anomaly to identify and exclude regions of quasi-stability. Future observations of this system can make use of our results through a scale factor relative to the currently measured minimum mass. Using a constant time lag tidal model (Hut 1981), we find that plausible tidal interactions within the system are insufficient to induce significant outward migration toward the theoretical stability limit. While current technologies are incapable of detecting exomoons in this system, we comment on the detectability of putative moons through Doppler monitoring within direct imaging observations in view of future research capacities.

The coming generation of galaxy surveys will provide measurements of galaxy clustering with unprecedented accuracy and data size, which will allow us to test cosmological models at much higher precision than achievable previously. This means that we must have more accurate theoretical predictions to compare with future observational data. As a first step towards more accurate modelling of the redshift space distortions (RSD) of small-scale galaxy clustering in modified gravity (MG) cosmologies, we investigate the validity of the so-called Skew-T (ST) probability distribution function (PDF) of halo pairwise peculiar velocities in these models. We show that, combined with the streaming model of RSD, the ST PDF substantially improves the small-scale predictions by incorporating skewness and kurtosis, for both $\Lambda$CDM and two leading MG models: $f(R)$ gravity and the DGP braneworld model. The ST model reproduces the velocity PDF and redshift-space halo clustering measured from MG $N$-body simulations very well down to $\sim 5 \, h^{-1}\mathrm{Mpc}$. In particular, we investigate the enhancements of halo pairwise velocity moments with respect to $\Lambda$CDM for a larger range of MG variants than previous works, and present simple explanations to the behaviours observed. By performing a simple Fisher analysis, we find a significnat increase in constraining power to detect modifications of General Relativity by introducing small-scale information in the RSD analyses.

Mrk926 is known to be a highly variable AGN with very broad line profiles. We studied the continuum and line profile variations of this object with high temporal resolution in order to determine its broad-line region structure and to derive its BH mass. We carried out a high-cadence spectroscopic variability campaign of Mrk926 with the 10m HET telescope, aided by photometric V-band data taken at the Wise Observatory, over a period of about five months. We extracted spectroscopic continuum and line light curves, and computed CCFs as well as velocity-resolved CCFs with respect to the combined continuum V-band light curve. The continuum luminosity of Mrk926 showed a drastic decrease during our campaign: it dropped to less than 50% of its original value within only 2.5 months. The spectra show complex and very broad Balmer line profiles, including outer Balmer satellites ranging from +-5000 to +-13,000 kms-1. The Ha, Hb, and HeI5876 line light curves are delayed relative to the continuum light curve. The Ha and Hb lines show two velocity-delay structures in the central part of their line profile (within +- 5000 kms-1), at ~10 and ~57 ld and at ~5 and ~48 ld, respectively. These structures might be interpreted as the signature of a line-emitting ring, inclined by ~50 deg to the line of sight and orbiting the BH at radii, R, of 33.5 and 26.5 ld. We determined continuum luminosities, log(lambda\,L_(lambda)/erg s-1), of 43.68 to 44.13, which are in good agreement with the established R_(BLR)-L_(AGN) relation. We derive a black hole mass of 1.1 +- 0.2 * 10^8 M_solar; this indicates a low Eddington ratio, which decreased from 8 to 3 percent during our campaign. The Balmer satellite components show a response to the continuum variations on the order of only 3-5 days. We attribute this to the central line segment and the Balmer satellites having different, spatially distinct regions of origin.

Time series data mining is an important field of research in the era of "Big Data". Next generation astronomical surveys will generate data at unprecedented rates, creating the need for automated methods of data analysis. We propose a method of light curve characterisation that employs a pipeline consisting of a neural network with a Long-Short Term Memory Variational Autoencoder architecture and a Gaussian mixture model. The pipeline performs extraction and aggregation of features from light curve segments into feature vectors of fixed length which we refer to as light curve "fingerprints". This representation can be readily used as input of down-stream machine learning algorithms. We demonstrate the proposed method on a data set of Rossi X-ray Timing Explorer observations of the galactic black hole X-ray binary GRS 1915+105, which was chosen because of its observed complex X-ray variability. We find that the proposed method can generate a representation that characterises the observations and reflects the presence of distinct classes of GRS 1915+105 X-ray flux variability. We find that this representation can be used to perform efficient classification of light curves. We also present how the representation can be used to quantify the similarity of different light curves, highlighting the problem of the popular classification system of GRS 1915+105 observations, which does not account for intermediate class behaviour.

Metrewave solar type-III radio bursts offer a unique means to study the properties of turbulence across the coronal heights.Theoretical models have shown that the apparent intensity and size of the burst sources evolve at sub-second scales under the influence of local turbulence. The properties of the evolution varies with frequency. However, observational studies remained difficult due to the lack of high fidelity imaging capabilities at these fine temporal scales simultaneously across wide spectral bands. I present a spectroscopic snapshot imaging (0.5 s, 160 kHz resolution) study of a type-III burst event across 80 - 200 MHz band.By modelling the temporal variability of the source sizes and intensity at every observation frequency, the characteristics of coronal turbulence is studied across a heliocentric height range of ~1.54 - 1.75 $R_\odot$. To understand the morphological evolution of the type-III source, a 2D Gaussian fitting procedure is used. The observed trends in the source area and integrated flux density are analysed in the framework of theoretical and data driven models.Results.The strength of density fluctuations ($\delta N/N$) in the corona is derived as a function of height (R). Combined with the archival low frequency data, $\delta N/N$ values across ~1.5 - 2.2 $R_\odot$ agree within a few factors. The burst decay time ($\tau_{decay}$) and the FWHM of the source showed a power-law dependency with frequency, roughly consistent with the results from data driven models. However,the values of $\tau_{decay}$ across frequency are higher than expected. The intrinsic sizes of the burst source were derived correcting for scatter broadening. This roughly matched the expected size of flux tubes at the coronal heights explored. I also report the observation of an intrinsic anti-phased pulsation in area and flux density of the source.

Astroparticle physics is undergoing a profound transformation, due to a series of extraordinary new results, such as the discovery of high-energy cosmic neutrinos with IceCube, the direct detection of gravitational waves with LIGO and Virgo, and many others. This white paper is the result of a collaborative effort that involved hundreds of theoretical astroparticle physicists and cosmologists, under the coordination of the European Consortium for Astroparticle Theory (EuCAPT). Addressed to the whole astroparticle physics community, it explores upcoming theoretical opportunities and challenges for our field of research, with particular emphasis on the possible synergies among different subfields, and the prospects for solving the most fundamental open questions with multi-messenger observations.

In this work we discuss the application of the Hamilton-Jacobi formalism on the scalar field implementation of Generalized Chaplygin Gas models. This corresponds to a Generalised Born-Infeld action for the scalar field, which in an initial fast-rolling phase mimics a matter-like behavior and in the final slow-rolling phase mimics a cosmological constant. In order to enrich the phenomenology of the model, we add an extra functional freedom, specified through a scalar potential for the field. Interestingly, we find that, due to the lifting induced by the non-standard kinetic term, an asymptotic de Sitter-like configuration can be obtained even for negative potentials. We show that at the background level, this model can easily mimic the {\Lambda}CDM model both with and without independent baryonic and radiation components.

Jets which are powered by an AGN are a crucial element in the study of their central black holes (BH) and their immediate surroundings. The formation of such jets is the subject of intense research, mainly based on the dichotomy presented by the two main jet launching scenarios $-$ the one from Blandford & Payne (1982), and the one from Blandford & Znajek (1977). In this work we study the prominent and nearby radio galaxy 3C 84 (NGC 1275) with 15, 43, and 86 GHz quasi-simultaneous VLBI observations. From these we determine the jet apex to be located $83\pm7\,\mu$as ($0.028-0.11$pc) upstream of the 86 GHz VLBI core, applying a two dimensional cross-correlation analysis. A byproduct of this analysis are spectral index maps, in which we identify a robust spectral index gradient in the north-south direction, for the first time at such high resolution, for the 43-86 GHz pair. The magnetic field strength at distances from the VLBI core comparable to measurements from the literature ($\sim10$ Schwarzschild radii) for other prominent AGN, like NGC 1052 and M 87, is computed to be $70-100$G. Implications for the magnetic field topology are also discussed.

We constrain the matter density $\Omega_{\mathrm{m}}$ and the amplitude of density fluctuations $\sigma_8$ within the $\Lambda$CDM cosmological model with shear peak statistics and angular convergence power spectra using mass maps constructed from the first three years of data of the Dark Energy Survey (DES Y3). We use tomographic shear peak statistics, including cross-peaks: peak counts calculated on maps created by taking a harmonic space product of the convergence of two tomographic redshift bins. Our analysis follows a forward-modelling scheme to create a likelihood of these statistics using N-body simulations, using a Gaussian process emulator. We include the following lensing systematics: multiplicative shear bias, photometric redshift uncertainty, and galaxy intrinsic alignment. Stringent scale cuts are applied to avoid biases from unmodelled baryonic physics. We find that the additional non-Gaussian information leads to a tightening of the constraints on the structure growth parameter yielding $S_8~\equiv~\sigma_8\sqrt{\Omega_{\mathrm{m}}/0.3}~=~0.797_{-0.013}^{+0.015}$ (68\% confidence limits), with a precision of 1.8\%, an improvement of ~38\% compared to the angular power spectra only case. The results obtained with the angular power spectra and peak counts are found to be in agreement with each other and no significant difference in $S_8$ is recorded. We find a mild tension of $1.5 \thinspace \sigma$ between our study and the results from Planck 2018, with our analysis yielding a lower $S_8$. Furthermore, we observe that the combination of angular power spectra and tomographic peak counts breaks the degeneracy between galaxy intrinsic alignment $A_{\mathrm{IA}}$ and $S_8$, improving cosmological constraints. We run a suite of tests concluding that our results are robust and consistent with the results from other studies using DES Y3 data.

We present a cosmological analysis using the second and third moments of the weak lensing mass (convergence) maps from the first three years of data (Y3) data of the Dark Energy Survey (DES). The survey spans an effective area of 4139 square degrees and uses the images of over 100 million galaxies to reconstruct the convergence field. The second moment of the convergence as a function of smoothing scale contains information similar to standard shear 2-point statistics. The third moment, or the skewness, contains additional non-Gaussian information. The data is analysed in the context of the $\Lambda$CDM model, varying 5 cosmological parameters and 19 nuisance parameters modelling astrophysical and measurement systematics. Our modelling of the observables is completely analytical, and has been tested with simulations in our previous methodology study. We obtain a 1.7\% measurement of the amplitude of fluctuations parameter $S_8\equiv \sigma_8 (\Omega_m/0.3)^{0.5} = 0.784\pm 0.013$. The measurements are shown to be internally consistent across redshift bins, angular scales, and between second and third moments. In particular, the measured third moment is consistent with the expectation of gravitational clustering under the $\Lambda$CDM model. The addition of the third moment improves the constraints on $S_8$ and $\Omega_{\rm m}$ by $\sim$15\% and $\sim$25\% compared to an analysis that only uses second moments. We compare our results with {\it Planck} constraints from the Cosmic Microwave Background (CMB), finding a $2.2$ \textendash $2.8\sigma$ tension in the full parameter space, depending on the combination of moments considered. The third moment independently is in $2.8\sigma$ tension with {\it Planck}, and thus provides a cross-check on analyses of 2-point correlations.

The observed anomalous steady decrease in surface temperature of the supernova remnant Cassiopeia A (Cas A), which was reported about ten years ago, has generated much debate. Several exotic cooling scenarios have been proposed using non-standard assumptions about the physics and evolution of this neutron star (NS). At present, significant corrections have been made to the observational data, which make it possible to numerically simulate the Cas A NS cooling process in the framework of the scenario of minimal neutrino cooling. If there is an additional source of cooling, such as axion emission, the steepness of the Cas A NS surface temperature drop will increase with the growth of the axion-nucleon interaction strength. This makes it possible to limit the minimum value of the axion decay constant $f_a$ using the condition that the NS surface temperature should be within the 99% confidence interval obtained from the observational data. Two types of axion models are considered: the Kim-Shifman-Weinstein-Zakharov -- KSVZ model and the Dean-Fischler-Srednitsky-Zhitnitsky --DFSZ model. The above criterion gives a lower limit on the axion decay constant, $f_a>3\times 10^7$ GeV and $f_a>4.5\times 10^8$ GeV for KSVZ and DFSZ axions, respectively.

One important global topological property of a spacetime manifold is orientability. It is widely believed that spatial orientability can only be tested by global journeys around the Universe to check for orientation-reversing closed paths. Since such global journeys are not feasible, theoretical arguments that combine universality of physical experiments with local arrow of time, CP violation and CPT invariance are usually offered to support the choosing of time- and space-orientable spacetime manifolds. The nonexistence of globally defined spinor fields on a non-orientable spacetime is another theoretical argument for orientability. However, it is conceivable that orientability can be put to test by local physical effects. In this paper, we show that it is possible to locally access spatial orientability of a spatially flat Friedmann--Robertson-Walker spacetime through quantum vacuum electromagnestic fluctuations. We argue that a putative non-orientability of the spatial sections of spatially flat FRW spacetime can be ascertained by the study of the stochastic motions of a charged particle or a point electric dipole under quantum vacuum electromagnetic fluctuations. In particular, the stochastic motions of a dipole permit the recognition of a presumed non-orientability of $3-$space in itself.

The accumulated energy density of the excited entropy modes in multiple field inflationary scenarios can play the role of dark matter. In the usual case of a flat field space without any turning trajectory, only light superhorizon entropy modes can be excited through the gravitational instability. In the case of a negatively curved field space, we show that subhorizon entropy modes can be excited as well through the tachyonic instability induced by the negative curvature of the field space. The latter allows for production of entropy modes with masses larger than or at the order of the Hubble expansion rate during inflation, leading to a new dark matter scenario. Due to the contribution of subhorizon modes, the corresponding spectral density has a peak at a scale smaller than its counterpart in the models based on a flat field space. This difference makes our model observationally distinguishable.

We derive fourth post-Newtonian (4PN) contributions to the Keplerian-type parametric solution associated with the conservative dynamics of eccentric, non-spinning compact binaries. The solution has been computed while ignoring certain zero-average, oscillatory terms arising due to 4PN tail effects. We provide explicit expressions for the parametric solution and various orbital elements in terms of the conserved energy, angular momentum and symmetric mass ratio. Canonical perturbation theory (along with the technique of Pade approximant) is used to incorporate the 4PN nonlocal-in-time tail effects within the action-angles framework. We then employ the resulting solution to obtain an updated inspiral-merger-ringdown (IMR) waveform that models the coalescence of non-spinning, moderately eccentric black hole binaries, influenced by arXiv:1709.02007. Our updated waveform is expected to be valid over similar parameter range as the above reference. We also present a related waveform which makes use of only the post-Newtonian equations and thus is valid only for the inspiral stage. This waveform is expected to work for a much larger range of eccentricity ($e_t \lesssim 0.85$) than our full IMR waveform (which assumes circularization of the binaries close to merger). We finally pursue preliminary data analysis studies to probe the importance of including the 4PN contributions to the binary dynamics while constructing gravitational waveform templates for eccentric mergers.

If dark matter annihilates with a velocity-dependent cross section within a subhalo, then the magnitude and angular distribution of the resulting photon signal will change. These effects are encoded in the $J$-factor. In this work we compute the $J$-factor for a variety of choices for the cross section velocity-dependence, and for a variety of choices for the dark matter profile, including generalized Navarro-Frenk-White (NFW), Einasto, Burkert and Moore. We find that, for cuspy profiles, the angular distribution of a future signal can potentially be used to determine the velocity-dependence of the annihilation cross section, and that these results are robust to changes in the form of the profile. For a cored profile, determining the velocity-dependence of the cross section is more difficult, but potentially still possible. Interestingly, we find that for a density profile with an inner slope power law steeper than 4/3, Sommerfeld-enhanced annihilation in the Coulomb limit leads to a divergence at the center, requiring a more detailed treatment of departure from the Coulomb limit.

Motivated by various excesses observed by the Fermi-LAT and AMS collaborations, we perform a detailed analysis of QCD uncertainties on particle spectra from dark-matter annihilation (or decay) into jets. When annihilated to SM particles, the final-state annihilation products undergo various complicated processes such as QED and QCD bremsstrahlung, hadronisation, and hadron decays. These processes contain some intrinsic uncertainties which are usually difficult to model and which are neglected in physical analyses. First, we perform several re-tunings of the fragmentation function parameters. Then, we estimate two kinds of uncertainties: {\it (i)} perturbative from QCD showers and {\it (ii)} non-perturbative from hadronisation function. The results are tabulated for a wide range of dark matter masses, $m_\chi \in [10, 10^5]~{\rm GeV}$, and annihilation channels. They can be found on Zenodo: https://doi.org/10.5281/zenodo.3764809

The orbital period loss of the compact binary systems is the first indirect evidence of gravitational waves which agrees well with Einstein's general theory of relativity to a very good accuracy. However, there is less than one percent uncertainty in the measurement of orbital period loss from the general reltivistic prediction. Perihelion precession of planets, Gravitational light bending and Shapiro delay are three other successful tests of general relativity theory. Though there are uncertainties in the measurements of those observations from the general reltivistic predictions as well. To resolve these uncertainties, we assume radiation of ultralight axions and light gauge boson particles of $L_i-L_j$ type from those systems which can be a possible candidate of fuzzy dark matter. In this article, we obtain bounds on new physics parameters from those astrophysical observations.

In the presence of interactions between neutrinos and dark matter (DM), DM can potentially be produced via freeze-in from the neutrino sector. We investigate the implications of such a scenario for the evolution of both DM and neutrinos in the early Universe, and show that the future cosmic neutrino detection experiment PTOLEMY might be sensitive to neutrino signals that originate from DM annihilation in this model.

The non-thermal acceleration of electrons and ions at an oblique, non-relativistic shock is studied using large scale particle-in-cell (PIC) simulations in one spatial dimension. Physical parameters are selected to highlight the role of electron pre-heating in injection and subsequent acceleration of particles at high Mach number shocks. Simulation results show evidence for the early onset of the diffusive shock acceleration process for both electrons and ions at a highly oblique subluminal shock. Ion acceleration efficiencies of $\lesssim 5\%$ are measured at late times, though this is not the saturated value.

In this paper, we investigate the potential of dark sirens by the space-borne atom interferometric gravitational-wave detectors to probe the Hubble constant. In the mid-frequency band, the sources live a long time. The motion of a detector around the Sun as well as in Earth orbit would induce large Doppler and reorientation effects, providing a precise angular resolution. Such precise localization for the GW sources makes it possible to observe the dark sirens with only one potential host galaxy, which are dubbed "golden dark sirens". We construct the catalogs of golden dark sirens and estimate that there are around 79 and 35 golden dark sirens of binary neutron stars (BNS) and binary black holes (BBH) that would be pass the detection threshold of AEDGE in 5 years. Our results show that with 5, 10, and all 79 golden dark BNS tracked by AEDGE one can constrain $H_0$ at 5.5\%, 4.1\%, and 1.8\% precision levels. With 5, 10, and all 35 golden dark BBH one can constrain $H_0$ at 2.2\%, 1.8\%, and 1.5\% precision levels, respectively. It suggests that only 5-10 golden dark BBH by AEDGE are sufficient to arbitrate the current tension between local and high-$z$ measurements of $H_0$.

Effects of electroweak phase transition (EWPT) in balance between baryon excess and the excess of stable quarks of new generation is studied. With the conservation of SU(2) symmetry and other quantum numbers, it makes possible sphaleron transitions between baryons, leptons and new of leptons and quarks. A definite relationship between the excess relative to baryon asymmetry is established. In passing by we also show the small, yet negligible dilution in the pre-existing dark matter density due the sphaleron transition.

Non-standard interactions (NSI) between neutrinos and electrons can significantly modify the decoupling of neutrinos from the plasma. These interactions have two effects on the overall picture: (i) they alter neutrino oscillations though matter effects and (ii) they modify the scattering and annihilation processes involving neutrinos and electrons and positrons. We study the role of non-universal and flavour-changing NSI in the decoupling and how they impact the determination of the effective number of neutrinos, $N_{eff}$. We examine the degeneracies between NSI parameters and we compare the expected sensitivity from future cosmological surveys with the current limits from terrestrial experiments. We outline the complementarity between both approaches.

The scope of this study is to understand the strength behaviour and fragment size of rocks during indirect, quasi-static, and dynamic tensile tests. Four rocks with different lithological characteristics namely, basalt, granite, sandstone and marble are selected. The Brazilian disc experiments are performed over a range of strain rates from 10-5 /s to 2.7x101 /s using a hydraulic loading frame and a split-Hopkinson bar. Over the range of strain rates, our measurements of dynamic strength increase are in good agreement with the universal theoretical scaling relationship of Kimberley et al. (2013). Dynamic fragmentation during a split tension mode receives very little attention and not much information is available about the generated fragment size distributions. The fragments fall into two distinct groups based on the nature of failure, coarser primary fragments and finer secondary fragments. The degree of fragmentation is assessed in terms of characteristic strain rate and is compared with existing theoretical tensile fragmentation models. The primary fragment size are less sensitive to strain rate, particularly at lower strain rates. The size of secondary fragment has a strong strain rate dependency over the entire testing range. Marble and sandstone are found to generate more pulverized secondary debris when compared to basalt and granite. Further, it is shown that the mean fragment sizes of primary and secondary fragments are well described by a power law function of strain rate.

Searching for departures from general relativity (GR) in more than one post-Newtonian (PN) phasing coefficients, called a \emph{multi-parameter test}, is known to be ineffective given the sensitivity of the present generation of gravitational-wave (GW) detectors. Strong degeneracies in the parameter space make the outcome of the test uninformative. We argue that Principal Component Analysis (PCA) can remedy this problem by constructing certain linear combinations of the original PN parameters that are better constrained by gravitational-wave observations. By analyzing binary black hole events detected during the first and second observing runs (O1 and O2) of LIGO/Virgo, we show that the two dominant principal components can capture the essence of a multi-parameter test. Combining five binary black hole mergers during O1/O2, we find that the dominant linear combination of the PN coefficients obtained from PCA is consistent with GR within the 0.38 standard deviation of the posterior distribution. Furthermore, using a set of simulated \emph{non-GR} signals in the three-detector LIGO-Virgo network with designed sensitivities, we find that the method is capable of excluding GR with high confidence as well as recovering the injected values of the non-GR parameters with good precision.