Papers for Wednesday, Aug 11 2021

We present long-duration numerical simulations of the tidal disruption of stars modelled with accurate stellar structures and spanning a range of pericentre distances, corresponding to cases where the stars are partially and completely disrupted. We substantiate the prediction that the late-time power-law index of the fallback rate $n_{\infty} \simeq -5/3$ for full disruptions, while for partial disruptions---in which the central part of the star survives the tidal encounter intact---we show that $n_{\infty} \simeq -9/4$. For the subset of simulations where the pericenter distance is close to that which delineates full from partial disruption, we find that a stellar core can reform after the star has been completely destroyed; for these events the energy of the zombie core is slightly positive, which results in late-time evolution from $n \simeq -9/4$ to $n \simeq -5/3$. We find that self-gravity can generate an $n(t)$ that deviates from $n_{\infty}$ by a small but significant amount for several years post-disruption. In one specific case with the stellar pericenter near the critical value, we find self-gravity also drives the re-collapse of the central regions of the debris stream into a collection of several cores while the rest of the stream remains relatively smooth. We also show that it is possible for the surviving stellar core in a partial disruption to acquire a circumstellar disc that is shed from the rapidly rotating core. Finally, we provide a novel analytical fitting function for the fallback rates that may also be useful in a range of contexts beyond TDEs.

Warm ionized and cold neutral outflows with velocities exceeding $100\,{\rm km\,s}^{-1}$ are commonly observed in galaxies and clusters. Theoretical studies however indicate that ram pressure from a hot wind, driven either by the central active galactic nucleus (AGN) or a starburst, cannot accelerate existing cold gas to such high speeds without destroying it. In this work we explore a different scenario, where cold gas forms in a fast, radiatively cooling outflow with temperature $T\lesssim 10^7\,{\rm K}$. Using 3D hydrodynamic simulations, we demonstrate that cold gas continuously fragments out of the cooling outflow, forming elongated filamentary structures extending tens of kiloparsecs. For a range of physically relevant temperature and velocity configurations, a ring of cold gas perpendicular to the direction of motion forms in the outflow. This naturally explains the formation of transverse cold gas filaments such as the blue loop and the horseshoe filament in the Perseus cluster. Based on our results, we estimate that the AGN outburst responsible for the formation of these two features drove bipolar outflows with velocity $>2,000\,{\rm km\,s}^{-1}$ and total kinetic energy $>8\times10^{57}\,{\rm erg}$ about $\sim10$ Myr ago. We also examine the continuous cooling in the mixing layer between hot and cold gas, and find that radiative cooling only accounts for $\sim10\%$ of the total mass cooling rate, indicating that observations of soft X-ray and FUV emission may significantly underestimate the growth of cold gas in the cooling flow of galaxy clusters.

I investigate the roles of cluster dynamics and massive binary evolution in producing stellar-remnant binary black hole (BBH) mergers over the cosmic time. To that end, dynamical BBH mergers are obtained from long-term direct N-body evolutionary models of $\sim10^4M_\odot$, pc-scale young massive clusters (YMC) evolving into moderate-mass open clusters (OC). Fast evolutionary models of massive isolated binaries (IB) yield BBHs from binary evolution. Population synthesis in a Model Universe is then performed, taking into account observed cosmic star-formation and enrichment histories, to obtain BBH-merger yields from these two channels observable at the present day and over cosmic time. The merging BBH populations from the two channels are combined by applying a proof-of-concept Bayesian regression chain, taking into account observed differential intrinsic BBH merger rate densities from the second gravitational-wave transient catalogue (GWTC-2). The analysis estimates an OB-star binary fraction of $f_{\rm Obin}\gtrsim60$% and a YMC formation efficiency of $f_{\rm YMC}\sim10^{-2}$, being consistent with recent optical observations and large scale structure formation simulations. The corresponding combined Model Universe present-day, differential intrinsic BBH merger rate density and the cosmic evolution of BBH merger rate density both agree well with those from GWTC-2. The analysis also suggests that despite significant 'dynamical mixing' at low redshifts, BBH mergers at high redshifts ($z_{\rm event}\gtrsim1$) could still be predominantly determined by binary-evolution physics. Caveats in the present approach and future improvements are discussed.

The interactions between radio jets and the interstellar medium play a defining role for the co-evolution of central supermassive black holes and their host galaxies, but observational constraints on these feedback processes are still very limited at redshifts $z > 2$. We investigate the radio-loud quasar PSO J352.4034-15.3373 at $z \sim 6$ at the edge of the Epoch of Reionization. This quasar is among the most powerful radio emitters and the first one with direct evidence of extended radio jets ($\sim$1.6 kpc) at these high redshifts. We analyze NOEMA and ALMA millimeter data targeting the CO (6-5) and [CII] far-infrared emission lines, respectively, and the underlying continuum. The broad $440\pm 80$ km s$^{-1}$ and marginally resolved [CII] emission line yields a systemic redshift of $z\!=\!5.832 \pm 0.001$. Additionally, we report a strong 215 MHz radio continuum detection, $88\pm 7$ mJy, using the GMRT. This measurement significantly improves the constraints at the low-frequency end of the spectral energy distribution of this quasar. In contrast to what is typically observed in high-redshift radio-quiet quasars, we show that cold dust emission alone cannot reproduce the millimeter continuum measurements. This is evidence that the strong synchrotron emission from the quasar contributes substantially to the emission even at millimeter (far-infrared in the rest-frame) wavelengths. This quasar is an ideal system to probe the effects of radio jets during the formation of a massive galaxy within the first Gyr of the Universe.

We present the color-magnitude diagrams and star formation histories (SFHs) of seven ultra-faint dwarf galaxies: Horologium 1, Hydra 2, Phoenix 2, Reticulum 2, Sagittarius 2, Triangulum 2, and Tucana 2, derived from high-precision Hubble Space Telescope photometry. We find that the SFH of each galaxy is consistent with them having created at least 80% of the stellar mass by $z\sim6$. For all galaxies, we find quenching times older than 11.5 Gyr ago, compatible with the scenario in which reionization suppresses the star formation of small dark matter halos. However, our analysis also reveals some differences in the SFHs of candidate Magellanic Cloud satellites, i.e., galaxies that are likely satellites of the Large Magellanic Cloud and that entered the Milky Way potential only recently. Indeed, Magellanic satellites show quenching times about 600 Myr more recent with respect to those of other Milky Way satellites, on average, even though the respective timings are still compatible within the errors. This finding is consistent with theoretical models that suggest that satellites' SFHs may depend on their host environment at early times, although we caution that within the error bars all galaxies in our sample are consistent with being quenched at a single epoch.

Hierarchical triple stars are ideal laboratories for studying the interplay between orbital dynamics and stellar evolution. Both stellar wind mass loss and three-body dynamics cooperate to destabilise triples, which can lead to a variety of astrophysical exotica. So far our understanding of their evolution was mainly built upon results from extensive binary-single scattering experiments. Starting from generic initial conditions, we evolve an extensive set of hierarchical triples using a combination of the triple evolution code TRES and an N-body code. We find that the majority of triples preserve their hierarchy throughout their evolution, which is in contradiction with the commonly adopted picture that unstable triples always experience a chaotic, democratic resonant interaction. The duration of the unstable phase is much longer than expected, so that stellar evolution cannot be neglected. Typically an unstable triple dissolve into a single star and a binary; sometimes democratically (the initial hierarchy is lost and the lightest body usually escapes), but also in a hierarchical way (the tertiary is ejected in a slingshot, independent of its mass). Collisions are common, and mostly involve the two original inner binary components still on the main-sequence. This contradicts the idea that collisions with a giant during democratic encounters dominate. Together with collisions in stable triples, we find that triple evolution is the dominant mechanism for stellar collisions in the Milky Way. Furthermore, our simulations produce runaway and walk-away stars with speeds up to several tens km/s, with a maximum of a few 100km/s. We suggest that destabilised triples can alleviate the tension behind the origin of the observed run-away stars. Lastly, we present a promising indicator to make general predictions for the fate of a specific triple, based on the initial inclination of the system.

NGC1275 is the Brightest Cluster Galaxy (BCG) in the Perseus cluster and hosts the active galactic nucleus (AGN) that is heating the central 100\,kpc of the intracluster medium (ICM) atmosphere via a regulated feedback loop. Here we use a deep 490ks Cycle-19 Chandra High-Energy Transmission Grating (HETG) observation of NGC1275 to study the anatomy of this AGN. The X-ray continuum is adequately described by an unabsorbed power-law with photon index $\Gamma\approx 1.9$, creating strong tension with the detected column of molecular gas seen via HCN and HCO$^+$ line absorption against the parsec-scale core/jet. This tension is resolved if we permit a composite X-ray source; allowing a column of $N_H\sim 8\times 10^{22}\,{\rm cm}^{-2}$ to cover $\sim 15$% of the X-ray emitter does produce a significant improvement in the statistical quality of the spectral fit. We suggest that the dominant unabsorbed component corresponds to the accretion disk corona, and the sub-dominant X-ray component is the jet working surface and/or jet cocoon that is expanding into clumpy molecular gas. We suggest that this may be a common occurence in BCG-AGN. We conduct a search for photoionized absorbers/winds and fail to detect such a component, ruling out columns and ionization parameters often seen in many other Seyfert galaxies. We detect the 6.4keV iron-K$\alpha$ fluorescence line seen previously by XMM-Newton and Hitomi. We describe an analysis methodology which combines dispersive HETG spectra, non-dispersive microcalorimeter spectra, and sensitive XMM-Newton/EPIC spectra in order to constrain (sub)arcsec-scale extensions of the iron-K$\alpha$ emission region.

The measurement of gravitational waves produced by binary black-hole mergers at the Advanced LIGO has encouraged extensive studies on the stochastic gravitational wave background. Recent studies have focused on gravitational wave sources made of the same species, such as mergers from binary primordial black holes or those from binary astrophysical black holes. In this paper, we study a new possibility --- the stochastic gravitational wave background produced by mergers of one primordial black hole and one astrophysical black hole. Such systems are necessarily present if primordial black holes exist. We study the isotropic gravitational wave background produced through the history of the Universe. We find it is very challenging to detect such a signal. We also demonstrate that it is improper to treat the gravitational waves produced by such binaries in the Milky Way as a directional stochastic background, due to a very low binary formation rate.

We report flux-resolved spectroscopic analysis of the active galactic nucleus (AGN) ESO 103--035 using \textit{NuSTAR} observations. Following an earlier work, we fit the spectra using a thermal Comptonization model with a relativistic reflection component to obtain estimates of the coronal temperature for two flux levels. The coronal temperature was found to increase from 24.0$^{+6.8}_{-3.4}$ to 55.3$^{+54.6}_{-7.2}$ keV (errors at 1-$\sigma$ confidence level) as the flux increased from $9.8$ to $11.9 \times 10^{-11}$ erg cm$^{-2}$ s$^{-1}$ in the 3--78 keV band. A marginal variation in the high energy photon index allows for both, a non-varying optical depth and for the optical depth to have varied by a factor of $\sim$2. This is in contrast to a previous work on \textit{NuSTAR} flux resolved spectroscopy of the AGN, Ark 564, where the temperature was found to decrease with flux along with a $10$\% variation in the optical depth. The results maybe understood in a framework where AGN variability is either dominated by coronal heating variation leading to correlated increase of temperature with flux and the opposite effect being seen when the variability is dominated by changes in the seed photon flux.

We present spatially resolved measurements of cool gas flowing into and out of the Milky Way (MW), using archival ultraviolet spectra of background quasars from the Hubble Space Telescope/Cosmic Origins Spectrograph. We co-add spectra of different background sources at close projected angular separation on the sky. This novel stacking technique dramatically increases the signal-to-noise ratio of the spectra, allowing detection of low column density gas (down to $EW$ > 2 mA). We identify absorption as inflowing or outflowing, by using blue/redshifted high velocity cloud (HVC) absorption components in the Galactocentric rest frame, respectively. The mass surface densities of inflowing and outflowing gas both vary by more than an order of magnitude across the sky, with mean values of $\langle \Sigma_{in}\rangle \gtrsim 10^{4.6\pm0.1}$ $M_{\odot}\,\mathrm{kpc}^{-2}$ for inflowing gas and $\langle \Sigma_{out}\rangle \gtrsim 10^{3.5\pm 0.1}$ $M_{\odot}\,\mathrm{kpc}^{-2}$ for outflowing gas, respectively. The mass flow rate surface densities (mass flow rates per unit area) also show large variation across the sky with $\langle \dot{\Sigma}(d)_{in}\rangle \gtrsim (10^{-3.6\pm0.1})(d/12 \mathrm{kpc})^{-1} M_{\odot}\,\mathrm{kpc}^{-2} \mathrm{yr}^{-1}$ for inflowing and $\langle \dot{\Sigma}(d)_{out}\rangle \gtrsim (10^{-4.8\pm0.1})(d/12\,\mathrm{kpc})^{-1} M_{\odot}\,\mathrm{kpc}^{-2}\,\mathrm{yr}^{-1}$ for outflowing gas, respectively. The regions with highest surface mass density of inflowing gas are clustered at smaller angular scales ($\theta < 40^\circ$). This indicates that most of the mass in inflowing gas is confined to small, well-defined structures, whereas the distribution of outflowing gas is spread more uniformly throughout the sky. Our study confirms that the MW is predominantly accreting gas, but is also losing a non-negligible mass of gas via outflow.

We explore the observational appearance of the merger of a low-mass star with a white dwarf (WD) binary companion. We are motivated by Schreiber et al. (2016), who found that multiple tensions between the observed properties of cataclysmic variables (CVs) and standard evolution models are resolved if a large fraction of CV binaries merge as a result of unstable mass transfer. Tidal disruption of the secondary forms a geometrically thick disk around the WD, which subsequently accretes at highly super-Eddington rates. Analytic estimates and numerical hydrodynamical simulations reveal that outflows from the accretion flow unbind a large fraction >~ 90% of the secondary at velocities ~500-1000 km/s within days of the merger. Hydrogen recombination in the expanding ejecta powers optical transient emission lasting about a month with a luminosity > 1e38 erg/s, similar to slow classical novae and luminous red novae from ordinary stellar mergers. Over longer timescales the mass accreted by the WD undergoes hydrogen shell burning, inflating the remnant into a giant of luminosity ~300-5000 L_sun, effective temperature T_eff ~ 3000 K and lifetime ~1e4-1e5 yr. We predict that ~1e3-1e4 Milky Way giants are CV merger products, potentially distinguishable by atypical surface abundances. We explore whether any Galactic historical slow classical novae are masquerading CV mergers by identifying four such post-nova systems with potential giant counterparts for which a CV merger origin cannot be ruled out. We address whether the historical transient CK Vul and its gaseous/dusty nebula resulted from a CV merger.

The next generation of the IceCube Neutrino Observatory, IceCube-Gen2, will constitute a much larger detector, increasing the rate of high-energy neutrinos. IceCube-Gen2 will address the long-standing questions about astrophysical accelerators. The experiment will also include a surface air-shower detector which will allow for measurements of cosmic rays in the energy region where a transition between Galactic and extragalactic accelerators is expected. As a baseline design for the surface detector, we consider a surface array above the optical in-ice array consisting of the same type of stations used for the IceTop enhancement, i.e., scintillation detectors and radio antennas. In order to better understand the capabilities of such an array, we performed simulations of its response to air showers, including both detector types. We will show the results of this simulation study and discuss the prospects for the surface array of IceCube-Gen2.

The recent detection of a very high energy (VHE) emission from Gamma-Ray Bursts (GRBs) above 100 GeV performed by the MAGIC and H.E.S.S. collaborations, has represented a significant, long-awaited result for the VHE astrophysics community. Although these results' scientific impact has not yet been fully exploited, the possibility to detect VHE gamma-ray signals from GRBs has always been considered crucial for clarifying the poorly known physics of these objects. Furthermore, the discovery of high-energy neutrinos and gravitational waves associated with astrophysical sources have definitively opened the era of multi-messenger astrophysics, providing unique insights into the physics of extreme cosmic accelerators. In the near future, the Cherenkov Telescope Array (CTA) will play a major role in these observations. Within this framework, the Large Size Telescopes (LSTs) will be the instruments best suited to significantly impact on short time-scale transients follow-up thanks to their fast slewing and large effective area. The observations of the early emission phase of a wide range of transient events with good sensitivity below 100 GeV will allow us to open new opportunities for time-domain astrophysics in an energy range not affected by selective absorption processes typical of other wavelengths. In this contribution, we will report about the observational program and first transients follow-up observations performed by the LST-1 telescope currently in its commissioning phase on La Palma, Canary Islands, the CTA northern hemisphere site.

The Cherenkov Telescope Array (CTA) is a next generation ground-based very-high-energy gamma-ray observatory that will allow for observations in the >10 GeV range with unprecedented photon statistics and sensitivity. This will enable the investigation of the yet-marginally explored physics of short-time-scale transient events. CTA will thus become an invaluable instrument for the study of the physics of the most extreme and violent objects and their interactions with the surrounding environment. The CTA Transient program includes follow-up observations of a wide range of multi-wavelength and multi-messenger alerts, ranging from compact galactic binary systems to extragalactic events such as gamma-ray bursts (GRBs), core-collapse supernovae and bright AGN flares. In recent years, the first firm detection of GRBs by current Cherenkov telescope collaborations, the proven connection between gravitational waves and short GRBs, as well as the possible neutrino-blazar association with TXS~0506+056 have shown the importance of coordinated follow-up observations triggered by these different cosmic signals in the framework of the birth of multi-messenger astrophysics. In the next years, CTA will play a major role in these types of observations by taking advantage of its fast slewing (especially for the CTA Large Size Telescopes), large effective area and good sensitivity, opening new opportunities for time-domain astrophysics in an energy range not affected by selective absorption processes typical of other wavelengths. In this contribution we highlight the common approach adopted by the CTA Transients physics working group to perform the study of transient sources in the very-high-energy regime.

The detection of gravitational waves from the merger of binary neutron stars events (GW170817, GW190425) and subsequent estimations of tidal deformability play a key role in constraining the behaviour of dense matter. In addition, massive neutron star candidates ($\sim 2 M_{\odot}$) along with NICER mass-radius measurements also, set sturdy constraints on the dense matter equation of state. Strict bounds from gravitational waves and massive neutron stars observations constrain the theoretical models of nuclear matter comportment at large density regimes. On the other hand, model parameters providing the highly dense matter response are bounded by nuclear saturation properties. This work analyses coupling parametrizations from two classes based on covariant density functional models: Non-Linear and Density-Dependent schemes. Considering these constraints together, we study possible models and parametrization schemes with the feasibility of exotic degrees of freedom in dense matter which go well with the astrophysical observations as well as the terrestrial laboratory experiments. We show that most parametrizations with non-linear schemes do not support the observations and experiments while density-dependent scheme goes well with both. Astrophysical observations are well explained if the inclusion of heavier non-strange baryons is considered as one fraction of the dense matter particle spectrum.

We examine the spatial clustering of blue horizontal branch (BHB) stars from the $\textit{u}$-band of the Canada-France Imaging Survey (CFIS, a component of the Ultraviolet Near-Infrared Optical Northern Survey, or UNIONS). All major groupings of stars are associated with previously known satellites, and among these is NGC 5466, a distant (16 kpc) globular cluster. NGC 5466 reportedly possesses a long stellar stream, although no individual members of the stream have previously been identified. Using both BHBs and more numerous red giant branch stars cross-matched to $\textit{Gaia}$ Data Release 2, we identify extended tidal tails from NGC 5466 that are both spatially and kinematically coherent. Interestingly, we find that this stream does not follow the same path as the previous detection at large distances from the cluster. We trace the stream across 31$^{\circ}$ of sky and show that it exhibits a very strong distance gradient ranging from 10 $<$ R$_{helio}$ $<$ 30 kpc. We compare our observations to simple dynamical models of the stream and find that they are able to broadly reproduce the overall path and kinematics. The fact that NGC 5466 is so distant, traces a wide range of Galactic distances, has an identified progenitor, and appears to have recently had an interaction with the Galaxy's disk, makes it a unique test-case for dynamical modelling of the Milky Way.

We study the radial and azimuthal mass distribution of the lensing galaxy in WFI2033-4723. Mindful of the fact that modeling results depend on modeling assumptions, we examine two very different recent models: simply parametrized (SP) models from the H0LiCOW collaboration, and pixelated free-form (FF) GLASS models. In addition, we fit our own models which are a compromise between the astrophysical grounding of SP, and the flexibility of FF approaches. Our models consist of two offset parametric mass components, and generate many solutions, all fitting the quasar point image data. Among other results, we show that to reproduce point image properties the lensing mass must be lopsided, but the origin of this asymmetry can reside in the main lens plane or along the line of sight. We also show that there is a degeneracy between the slope of the density profile and the magnitude of external shear, and that the models from various modeling approaches are connected not by the mass sheet degeneracy, but by a more generalized transformation. Finally, we discuss interpretation degeneracy which afflicts all mass modeling: inability to correctly assign mass to the main lensing galaxy vs. nearby galaxies or line of sight structures. While this may not be a problem for the determination of $H_0$, interpretation degeneracy may become a major issue for the detailed study of galaxy structure.

Images of the solar corona by extreme-ultraviolet (EUV) telescopes reveal elegant arches of glowing plasma that trace the corona's magnetic field. Typically, these loops are preferentially illuminated segments of an arcade of vaulted field lines and such loops are often observed to sway in response to nearby solar flares. A flurry of observational and theoretical effort has been devoted to the exploitation of these oscillations with the grand hope that seismic techniques might be used as probes of the strength and structure of the corona's magnetic field. The commonly accepted viewpoint is that each visible loop oscillates as an independent entity and acts as a one-dimensional (1D) wave cavity for magnetohydrodynamic (MHD) kink waves. We argue that for many events, this generally accepted model for the wave cavity is fundamentally flawed. In particular, the 3D magnetic arcade in which the bright loop resides participates in the oscillation. Thus, the true wave cavity is larger than the individual loop and inherently multidimensional. We derive the skin depth of the near-field response for an oscillating loop and demonstrate that most loops are too close to other magnetic structures to oscillate in isolation. Further, we present a simple model of a loop embedded within an arcade and explore how the eigenmodes of the arcade and the eigenmodes of the loop become coupled. In particular, we discuss how distinguishing between these two types of modes can be difficult when the motions within the arcade are often invisible.

We use the angular Two Point Correlation Function (TPCF) to investigate the hierarchical distribution of young star clusters in 12 local (3--18 Mpc) star-forming galaxies using star cluster catalogues obtained with the \textit{Hubble Space Telescope} (\textit{HST}) as part of the Treasury Program LEGUS (Legacy ExtraGalactic UV Survey). The sample spans a range of different morphological types, allowing us to infer how the physical properties of the galaxy affect the spatial distribution of the clusters. We also prepare a range of physically motivated toy models to compare with and interpret the observed features in the TPCFs. We find that, conforming to earlier studies, young clusters ($T \la 10\, \mathrm{Myr}$) have power-law TPCFs that are characteristic of fractal distributions with a fractal dimension $D_2$, and this scale-free nature extends out to a maximum scale $l_{\mathrm{corr}}$ beyond which the distribution becomes Poissonian. However, $l_{\mathrm{corr}}$, and $D_2$ vary significantly across the sample, and are correlated with a number of host galaxy physical properties, suggesting that there are physical differences in the underlying star cluster distributions. We also find that hierarchical structuring weakens with age, evidenced by flatter TPCFs for older clusters ($T \ga 10\, \mathrm{Myr}$), that eventually converges to the residual correlation expected from a completely random large-scale radial distribution of clusters in the galaxy in $\sim 100 \, \mathrm{Myr}$. Our study demonstrates that the hierarchical distribution of star clusters evolves with age, and is strongly dependent on the properties of the host galaxy environment.

Fireball networks are used to recover meteorites, with the context of orbits. Observations from these networks cover the bright flight, where the meteoroid is luminescent, but to recover a fallen meteorite, these observations must often be predicted forward in time to the ground to estimate an impact position. This darkflight modelling is deceptively simple, but there is hidden complexity covering the precise interactions between the meteorite and the (usually active) atmosphere. We describe the method and approach used by the Desert Fireball Network, detailing the issues we have addressed, and the impact that factors such as shape, mass and density have on the predicted fall position. We illustrate this with a case study of Murrili meteorite fall that occurred into Lake Eyre-Kati Thanda in 2015. The fall was very well observed from multiple viewpoints, and the trajectory was steep, with a low altitude endpoint, such that the darkflight was relatively short. Murrili is 1.68 kg with a typical ordinary chondrite density, but with a somewhat flattened shape compared to a sphere, such that there are discrepancies between sphere-based predictions and the actual recovery location. It is notable that even in this relatively idealised darkflight scenario, modelling using spherical shaped projectiles resulted in a significant distance between predicted fall position and recovered meteorite.

We study the parameter space of the effective (with two scalars) models of cosmological inflation and primordial black hole (PBH) formation in the modified $(R+R^2)$ supergravity. Our models describe double inflation, whose first stage is driven by Starobinsky's scalaron coming from the $R^2$ gravity, and whose second stage is driven by another scalar belonging to the supergravity multiplet. The ultra-slow-roll regime between the two stages leads a large peak (enhancement) in the power spectrum of scalar perturbations, which results in efficient PBH formation. Both inflation and PBH formation are generic in our models, while those PBH can account for a significant part or the whole of dark matter. Some of the earlier proposed models in the same class are in tension (over $3\sigma$) with the observed value of the scalar tilt $n_s$, so that we study more general models with more parameters, and investigate the dependence of the cosmological tilts $(n_s,r)$ and the scalar power spectrum enhancement upon the parameters. The PBH masses and their density fraction (as part of dark matter) are also calculated. A good agreement (between $2\sigma$ and $3\sigma$) with the observed value of $n_s$ requires fine tuning of the parameters, and it is only realized in the so-called $\delta$-models. Our models offer the (super)gravitational origin of inflation, PBH and dark matter together, and may be confirmed or falsified by future precision measurements of the cosmic microwave background radiation and PBH-induced gravitational waves.

We present observations of the Eridanus supergroup obtained with the Australian Square Kilometre Array Pathfinder (ASKAP) as part of the pre-pilot survey for the Widefield ASKAP L-band Legacy All-sky Blind Survey (WALLABY). The total number of detected HI sources is 55, of which 12 are background galaxies not associated with the Eridanus supergroup. Two massive HI clouds are identified and large HI debris fields are seen in the NGC 1359 interacting galaxy pair, and the face-on spiral galaxy NGC 1385. We describe the data products from the source finding algorithm and present the basic parameters. The presence of distorted HI morphology in all detected galaxies suggests ongoing tidal interactions within the subgroups. The Eridanus group has a large fraction of HI deficient galaxies as compared to previously studied galaxy groups. These HI deficient galaxies are not found at the centre of the group. We find that galaxies in the Eridanus supergroup do not follow the general trend of the atomic gas fraction versus stellar mass scaling relation, which indicates that the scaling relation changes with environmental density. In general, the majority of these galaxies are actively forming stars.

We present the Australian Square Kilometre Array Pathfinder (ASKAP) WALLABY pre-pilot observations of two `dark' HI sources (with HI masses of a few times 10^8 Msol and no known stellar counterpart) that reside within 363 kpc of NGC 1395, the most massive early-type galaxy in the Eridanus group of galaxies. We investigate whether these `dark' HI sources have resulted from past tidal interactions or whether they are an extreme class of low surface brightness galaxies. Our results suggest that both scenarios are possible, and not mutually exclusive. The two `dark' HI sources are compact, reside in relative isolation and are more than 159 kpc away from their nearest HI-rich galaxy neighbour. Regardless of origin, the HI sizes and masses of both `dark' HI sources are consistent with the HI size-mass relationship that is found in nearby low-mass galaxies, supporting the possibility that these HI sources are an extreme class of low surface brightness galaxies. We identified three analogues of candidate primordial `dark' HI galaxies within the TNG100 cosmological, hydrodynamic simulation. All three model analogues are dark matter-dominated, have assembled most of their mass 12-13 Gyr ago, and have not experienced much evolution until cluster infall 1-2 Gyr ago. Our WALLABY pre-pilot science results suggest that the upcoming large area HI surveys will have a significant impact on our understanding of low surface brightness galaxies and the physical processes that shape them.

We use high-resolution ASKAP observations of galaxies in the Eridanus supergroup to study their HI, angular momentum and star formation properties, as part of the WALLABY pre-pilot survey efforts. The Eridanus supergroup is composed of three sub-groups in the process of merging to form a cluster. The main focus of this study is the Eridanus (or NGC 1395) sub-group. The baryonic specific angular momentum - baryonic mass ($j_{\mathrm{b}} - M_{\mathrm{b}}$) relation for the Eridanus galaxies is observed to be an unbroken power law of the form $j_{\mathrm{b}} \propto M_{\mathrm{b}}^{0.57 \pm 0.05}$, with a scatter of $\sim 0.10 \pm 0.01$ dex, consistent with previous works. We examine the relation between the atomic gas fraction, $f_{\mathrm{atm}}$, and the integrated atomic disc stability parameter $q$ (the $f_{\mathrm{atm}} - q$ relation), and find that the Eridanus galaxies deviate significantly from the relation owing to environmental processes such as tidal interactions and ram-pressure affecting their HI gas. We find that a majority of the Eridanus galaxies are HI deficient compared to normal star-forming galaxies in the field. We also find that the star formation among the Eridanus galaxies may be suppressed owing to their environment, thus hinting at significant levels of pre-processing within the Eridanus sub-group, even before the galaxies have entered a cluster-like environment.

Exoplanets that receive stellar irradiance of approximately Earth's or less have been discovered and many are suitable for spectral characterization. Here we focus on the temperate planets that have massive H2-dominated atmospheres, and trace the chemical reactions and transport following the photodissociation of H2O, CH4, NH3, and H2S, with K2-18 b, PH2 b, and Kepler-167 e representing temperate/cold planets around M and G/K stars. We find that NH3 is likely depleted by photodissociation to the cloud deck on planets around G/K stars but remains intact in the middle atmosphere of planets around M stars. A common phenomenon on temperate planets is that the photodissociation of NH3 in presence of CH4 results in HCN as the main photochemical product. The photodissociation of CH4 together with H2O leads to CO and CO2, and the synthesis of hydrocarbon is suppressed. Temperate planets with super-solar atmospheric metallicity and appreciable internal heat may have additional CO and CO2 from the interior and less NH3 and thus less HCN. Our models of K2-18 b can explain the transmission spectrum measured by Hubble, and indicate that future observations in 0.5-5.0 um would provide the sensitivity to detect the equilibrium gases CH4, H2O, and NH3, the photochemical gas HCN, as well as CO2 in some cases. Temperate and H2-rich exoplanets are thus laboratories of atmospheric chemistry that operate in regimes not found in the Solar System, and spectral characterization of these planets in transit or reflected starlight promises to greatly expand the types of molecules detected in exoplanet atmospheres.

We have studied a sample of 101 bright 2MASS galaxies from the Large Galaxy Atlas (LGA), whose morphologies span from early to late-types. We have generated estimates for structural parameters through a two-dimensional (2D) surface brightness photometric decomposition in the three 2MASS bands (J, H, Ks). This work represents a detailed multi-component photometric study of nearby galaxies. We report total magnitudes, effective radii, concentration indices, among other parameters, in the three 2MASS bands. We found that the integrated total magnitudes of early-type galaxies (ETGs) measured on 2MASS LGA mosaics are ~0.35 mag dimmer, when compared with images generated from IRSA image tiles service; nevertheless, when comparing late-type galaxies (LTGs) we did not find any difference. Therefore, for ETGs we present the results derived on IRSA image tiles, while for LTGs we used data from the LGA mosaics. Additionally, by combining these structural parameters with scaling relations and kinematic data, we separated classical bulges from pseudobulges. We found that ~40 % of the objects in our sample are classified as pseudobulges, which are found preferentially in LTGs. Also, our findings confirm trends reported earlier in the distributions for some physical parameters, such as S\'ersic index, B/T and q ratios. In general, our results are in agreement with previous one-dimensional studies. In a companion paper, we revise some of the scaling relations among global galaxy properties, as well as their interrelation with Supermassive Black Holes.

The Cherenkov Telescope Array (CTA) Observatory, with dozens of telescopes located in both the Northern and Southern Hemispheres, will be the largest ground-based gamma-ray observatory and will provide broad energy coverage from 20 GeV to 300 TeV. The large effective area and field-of-view, coupled with the fast slewing capability and unprecedented sensitivity, make CTA a crucial instrument for the future of ground-based gamma-ray astronomy. To maximise the scientific return, the array will send alerts on transients and variable phenomena (e.g. gamma-ray burst, active galactic nuclei, gamma-ray binaries, serendipitous sources). Rapid and effective communication to the community requires a reliable and automated system to detect and issue candidate science alerts. This automation will be accomplished by the Science Alert Generation (SAG) pipeline, a key system of the CTA Observatory. SAG is part of the Array Control and Data Acquisition (ACADA) working group. The SAG working group develops the pipelines to perform data reconstruction, data quality monitoring, science monitoring and real-time alert issuing during observations to the Transients Handler functionality of ACADA. SAG is the system that performs the first real-time scientific analysis after the data acquisition. The system performs analysis on multiple time scales (from seconds to hours). \abrb{SAG must issue candidate science alerts within} 20 seconds from the data taking and with sensitivity at least half of the CTA nominal sensitivity. These challenging requirements must be fulfilled by managing trigger rates of tens of kHz from the arrays. Dedicated and highly optimised software and hardware architecture must thus be designed and tested. In this work, we present the general architecture of the ACADA-SAG system.

We study the evolution of extremely metal-poor AGB stars, with metallicities down to [Fe/H]=-5, to understand the main evolutionary properties, the efficiency of the processes able to alter their surface chemical composition and to determine the gas and dust yields. We calculate two sets of evolutionary sequences of stars in the 1-7.5Msun mass range, evolved from the pre-main sequence to the end of the AGB phase. To explore the extremely metal-poor chemistries we adopted the metallicities Z=3x10^{-5} and Z=3x10^{-7} which correspond, respectively to [Fe/H]=-3 and [Fe/H]=-5. The results from stellar evolution modelling are used to calculate the yields of the individual chemical species. We also modelled dust formation in the wind, to determine the dust produced by these objects. The evolution of AGB stars in the extremely metal-poor domain explored here proves tremendously sensitive to the initial mass of the star. M<2Msun stars experience several third dredge-up events, which favour the gradual surface enrichment of C12 and the formation of significant quantities of carbonaceous dust, of the order of 0.01Msun. The C13 and nitrogen yiel are found to be significantly smaller than in previous explorations of low-mass, metal-poor AGB stars, owing to the weaker proton ingestion episodes experienced during the initial AGB phases. M>5Msun stars experience hot bottom burning and their surface chemistry reflects the equilibria of a very advanced proton-capture nucleosynthesis; little dust production takes place in their wind. Intermediate mass stars experience both third dredge-up and hot bottom burning: they prove efficient producers of nitrogen, which is formed by proton captures on C12 nuclei of primary origin dredged-up from the internal regions.

Recent observations have shown that repeating Fast Radio Bursts (FRBs) exhibit band-limited emission, whose frequency-dependent amplitude can be modeled using a Gaussian function. In this analysis, we show that banded emission of FRBs can lead to incompleteness across the observing band. This biases the detected sample of bursts and can explain the various shapes of cumulative energy distributions seen for repeating FRBs. We assume a Gaussian shape of the burst spectra and used simulations to demonstrate the above bias using an FRB 121102-like example. We recovered energy distributions that showed a break in power-law and flattening of power-law at low energies, based on the fluence threshold of the observations. We provide recommendations for single-pulse searches and analysis of repeating FRBs to account for this incompleteness. Primarily, we recommend that burst spectra should be modeled to estimate the intrinsic fluence and bandwidth of the burst robustly. Also, bursts that lie mainly within the observing band should be used for analyses of energy distributions. We show that the bimodality reported in the distribution of energies of FRB 121102 by Li et al. (2021) disappears when burst bandwidth, instead of the center frequency of the observation, is used to estimate energy. Sub-banded searches will also aid in detecting band-limited bursts. All the analysis scripts used in this work are available in a Github repository.

The Cherenkov Telescope Array (CTA) will be the next generation ground-based observatory for very-high-energy (VHE) gamma-ray astronomy, with the deployment of tens of highly sensitive and fast-reacting Cherenkov telescopes. It will cover a wide energy range (20 GeV - 300 TeV) with unprecedented sensitivity. To maximize the scientific return, the observatory will be provided with an online software system that will perform the first analysis of scientific data in real-time. This study investigates the precision and accuracy of available science tools and analysis techniques for the short-term detection of gamma-ray sources, in terms of sky localization, detection significance and, if significant detection is achieved, a first estimation of the integral photon flux. The scope is to evaluate the feasibility of the algorithms' implementation in the real-time analysis of CTA. In this contribution we present a general overview of the methods and some of the results for the test case of the short-term detection of a gamma-ray burst afterglow, as the VHE counterpart of a gravitational wave event.

In the multi-messenger era, space and ground-based observatories usually develop real-time analysis (RTA) pipelines to rapidly detect transient events and promptly share information with the scientific community to enable follow-up observations. These pipelines can also react to science alerts shared by other observatories through networks such as the Gamma-Ray Coordinates Network (GCN) and the Astronomer's Telegram (ATels). AGILE is a space mission launched in 2007 to study X-ray and gamma-ray phenomena. This contribution presents the technologies used to develop two types of AGILE pipelines using the RTApipe framework and an overview of the main scientific results. The first type performs automated analyses on new AGILE data to detect transient events and automatically sends AGILE notices to the GCN network. Since May 2019, this pipeline has sent more than 50 automated notices with a few minutes delay since data arrival. The second type of pipeline reacts to multi-messenger external alerts (neutrinos, gravitational waves, GRBs, and other transients) received through the GCN network and performs hundreds of analyses searching for counterparts in all AGILE instruments' data. The AGILE Team uses these pipelines to perform fast follow-up of science alerts reported by other facilities, which resulted in the publishing of several ATels and GCN circulars.

In this letter we calculate the gravitational waves (GWs) emitted from a small binary (SB) by solving the Teukolsky equation in the background of a massive exotic compact object (ECO) which is phenomenologically described by a Schwarzschild geometry with a reflective boundary condition at its ``would-be'' horizon. The ``continuous echo" waves propagating to infinity due to reflectivity of ECO at its ``would-be'' horizon provide an exquisite probe to the nature of the ECO's horizon.

The Cherenkov Telescope Array (CTA) will be the next generation very-high-energy gamma-ray observatory. CTA is expected to provide substantial improvement in accuracy and sensitivity with respect to existing instruments thanks to a tenfold increase in the number of telescopes and their state-of-the-art design. Detailed Monte Carlo simulations are used to further optimise the number of telescopes and the array layout, and to estimate the observatory performance using updated models of the selected telescope designs. These studies are presented in this contribution for the two CTA stations located on the island of La Palma (Spain) and near Paranal (Chile) and for different operation and observation conditions.

The ASTRI Mini-Array is an international collaboration led by the Italian National Institute for Astrophysics (INAF), aiming to construct and operate an array of nine Imaging Atmospheric Cherenkov Telescopes (IACTs) to study gamma-ray sources at very high energy (TeV) and to perform stellar intensity interferometry observations. This contribution describes the design and the technologies used by the ASTRI team to implement the Online Observation Quality System (OOQS). The main objective of the OOQS is to perform data quality analyses in real-time during Cherenkov and intensity interferometry observations to provide feedback to both the Central Control System and the Operator. The OOQS performs the analysis of key data quality parameters and can generate alarms to other sub-systems for a fast reaction to solve critical conditions. The results from the data quality analyses are saved into the Quality Archive for further investigations. The Operator can visualise the OOQS results through the Operator Human Machine Interface as soon as they are produced. The main challenge addressed by the OOQS design is to perform online data quality checks on the data streams produced by nine telescopes, acquired by the Array Data Acquisition System and forwarded to the OOQS. In the current OOQS design, the Redis in-memory database manages the data throughput generated by the telescopes, and the Slurm workload scheduler executes in parallel the high number of data quality analyses.

Describing how the properties of the interstellar medium combine across size-scales is crucial for understanding star formation scaling laws and connecting Galactic and extragalactic data of molecular clouds. We describe how the statistical structure of clouds, and its connection to star formation, changes from sub-parsec to kiloparsec scales in a complete region within the Milky Way disk. We build a census of molecular clouds within 2 kpc from the Sun using literature. We examine the dust-based column density probability distributions (N-PDFs) of the clouds and their relation to star formation traced by young stellar objects (YSOs). We then examine our survey region from the outside, within apertures of varying sizes, and describe how the N-PDFs and their relation to star formation changes with the size-scale. The N-PDFs of the clouds are not well described by any single simple model; use of any single model may bias the interpretation of the N-PDFs. The top-heaviness of the N-PDFs correlates with star formation activity, and the correlation changes with Galactic environment (spiral-/inter-arm regions). We find that the density contrast of clouds may be more intimately linked to star formation than the dense gas mass fraction. The aperture-averaged N-PDFs vary with the size-scale and are more top-heavy for larger apertures. The top-heaviness of the aperture N-PDFs correlates with star formation activity up to roughly 0.5 kpc size-scale, depending on the environment. Our results suggest that the relations between cloud structure and star formation are environment specific and best captured by relative quantities (e.g., the density contrast). Finally, we show how the density structures of individual clouds give rise to a kpc-scale Kennicutt-Schmidt relationship as a combination of sampling effects and blending of different galactic environments.

The cosmological luminosity-distance can be measured from gravitational wave (GW) standard sirens, free of astronomical distance ladders and the associated systematics. However, it may still contain systematics arising from various astrophysical, cosmological and experimental sources. With the large amount of dark standard sirens of upcoming third generation GW experiments, such potential systematic bias can be diagnosed and corrected by statistical tools of the large scale structure of the universe. We estimate that, by cross-correlating the dark siren luminosity-distance space distribution and galaxy redshift space distribution, multiplicative error $m$ in the luminosity distance measurement can be constrained with $1\sigma$ uncertainty $\sigma_m\sim 0.1$. This is already able to distinguish some binary black hole origin scenarios unambiguously. Significantly better constraints and therefore more applications may be achieved by more advanced GW experiments.

Ground-based gamma-ray astronomy aims at reconstructing the energy and direction of gamma rays from the extensive air showers they initiate in the atmosphere. Imaging Atmospheric Cherenkov Telescopes (IACT) collect the Cherenkov light induced by secondary charged particles in extensive air showers (EAS), creating an image of the shower in a camera positioned in the focal plane of optical systems. This image is used to evaluate the type, energy and arrival direction of the primary particle that initiated the shower. This contribution shows the results of a novel reconstruction method based on likelihood maximization. The novelty with respect to previous likelihood reconstruction methods lies in the definition of a likelihood per single camera pixel, accounting not only for the total measured charge, but also for its development over time. This leads to more precise reconstruction of shower images. The method is applied to observations of the Crab Nebula acquired with the Large Size Telescope prototype (LST-1) deployed at the northern site of the Cherenkov Telescope Array.

The occurrence of quark matter at the center of neutron stars is still in debate. This study defines some semi-empirical parameters that quantify the occurrence and the amount of quark matter at star interiors. These parameters show semi-universal relations across all the EoS. One parameter depends on the shifting of the keplerian mass-radius curve from the static one and shows it is a constant across all EoS. The Z-parameter shows how tidal deformability depends on the quark content of the star and the stiffness of the EoS. The quark content of the star also affects the compactness of the star, and its dependence is almost universal. The empirical parameter gives a bound on the quark content of the star and shows that if the amount of the quark content increases, the stars are likely to collapse into a black hole. It is seen that the change in the mass and radius after PT is linearly proportional to the mass of the parent NS. Given a hadronic EoS, bag constant, and quark coupling constant, one can have a critical mass of the neutron star and the maximum mass of the hybrid star for phase transition without any baryonic mass loss.

A number of cosmic-ray observatories have measured a change in both phase and amplitude of the dipole component in the distribution of cosmic-ray arrival directions above a primary energy of 100 TeV. We focus on probing the cosmic-ray dipole and multipole evolution in the energy region of mutli TeV to beyond PeV with a future large-area gamma-ray observatory, such as the Southern Wide-field Gamma-ray Observatory (SWGO). The ability to discriminate between different mass groups is essential to understand the origin of this evolution. Through a consideration of the energy and mass resolution for cosmic-ray detection by such an observatory, we estimate its separation power for decomposing the full-particle anisotropy into mass groups. In particular, we explore the feasibility of probing the dipole evolution with rigidity with SWGO. In this way, we demonstrate the great potential that this instrument offers for providing a deeper understanding of the origin of the cosmic-ray anisotropy.

Stellar dust grains are predominantly composed of mineralic, anorganic material forming in the circumstellar envelopes of oxygen-rich AGB stars. However, the initial stage of the dust synthesis, or its nucleation, is not well understood. In particular, the chemical nature of the nucleating species, represented by molecular clusters, is uncertain. We investigated the vertical and adiabatic ionization energies of four different metal-oxide clusters by means of density functional theory. They included clusters of magnesia (MgO)$_n$, silicon monoxide (SiO)$_n$, alumina (Al$_2$O$_3$)$_n$, and titania (TiO$_2$)$_n$ with stoichiometric sizes of $n$=1$-$8. The magnesia, alumina, and titania clusters showed relatively little variation in their ionization energies with respect to the cluster size n: 7.1$-$8.2 eV for (MgO)$_n$, 8.9$-$10.0 eV for (Al$_2$O$_3$)$_n$, and 9.3$-$10.5 eV for (TiO$_2$)$_n$. In contrast, the (SiO)$_n$ ionization energies decrease with size $n$, starting from 11.5 eV for $n$=1, and decreasing to 6.6 eV for $n$=8. Therefore, we set constraints on the stability limit for neutral metal-oxide clusters to persist ionization through radiation or high temperatures and for the nucleation to proceed via neutral-neutral reactions.

Spiral structure (both flocculent and Grand Design types) is very rarely observed in dwarf galaxies because the formation of spiral arms requires special conditions. In this work we analyze the sample of about 40 dS-galaxies found by scanning by eye the images of late-type galaxies with $m_B<15^m$ and $M_B>-18^m$ and photometric diameter $D_{25}<12$~kpc. We found that apart from the lower average gas (HI) fraction the other properties of dS-galaxies including the presence of a bar and the isolation index do not differ much from those for dwarf Irr or Sm-types of similar luminosity and rotation velocity (or specific angular momentum).There are practically no dS-galaxies with rotation velocity below 50\,--\,60~km\,sec$^{-1}$. To check the conditions of formation of spiral structure in dwarf galaxies we carried out a series of N-body/hydrodynamic simulations of low-mass stellar-gaseous discy galaxies by varying the model kinematic parameters of discs, their initial thickness, relative masses and scale lengths of stellar and gaseous disc components, and stellar-to-dark halo masses. We came to conclusion that the gravitational mechanism of spiral structure formation is effective only for thin stellar discs, which are non-typical for dwarf galaxies, and for not too slowly rotating galaxies. Therefore, only a small fraction of dwarf galaxies with stellar/gaseous discs have spiral or ring structures. The thicker stellar disc, the more gas is required for the spiral structure to form. The reduced gas content in many dS-galaxies compared to non-spiral ones may be a result of more efficient star formation due to a higher volume gas density thank to the thinner stellar/gaseous discs.

Radio detection of air showers produced by ultra-high energy cosmic rays is a cost-effective technique for the next generation of sparse arrays. The performance of this technique strongly depends on the environmental background, which has different constituents, namely anthropogenic radio frequency interference, synchrotron galactic radiation and others. These components have recognizable features, which can help for background suppression. A powerful method for handling this is the application of convolution neural networks with a specific architecture called autoencoder. By suppressing unwanted signatures, the autoencoder keeps the signal-like ones. We have successfully developed and trained an autoencoder, which is now applied to the data from Tunka-Rex. We show the procedures of the training and optimization of the network including benchmarks of different architectures. Using the autoencoder, we improved the standard analysis of Tunka-Rex in order to lower the threshold of the detection. This enables the reconstructing of sub-threshold events with energies lower than 0.1 EeV with satisfactory angular and energy resolutions.

Air-shower radio arrays operate in low signal-to-noise ratio conditions, which complicates the autonomous measurement of air-shower signals without using an external trigger from optical or scintillator detectors. A simple threshold trigger for radio detector can be efficiently applied onlyin radio-quiet conditions, because for other cases this trigger detects a high fraction of noise pulses. In the present work, we study aspects of independent air-shower detection by dense antenna clusters with a complex real-time trigger system. For choosing the optimal procedures for the real-time analysis, we study the dependence between trigger efficiency, count rate, detector hardware and geometry. For this study, we develop a framework for testing various methods of signal detection and noise filtration for arrays with various specifications and the hardware implementation of these methods based on field programmable gate arrays. The framework provides flexible settings for the management of station-level and cluster-level steps of detecting the signal, optimized for the hardware implementation for real-time processing. It includes data-processing tools for the initialconfiguration and tests on pre-recorded data, tools for configuring the trigger architecture andtools for preliminary estimates of the trigger efficiency at given thresholds of cosmic-ray energyand air-shower pulse amplitude. We show examples of the trigger pipeline developed with this framework and discuss the results of tests on simulated data.

Tunka-Rex (Tunka Radio Extension) was a detector for ultra-high energy cosmic rays measuring radio emission for air showers in the frequency band of 30-80 MHz, operating in 2010s. It provided an experimental proof that sparse radio arrays can be a cost-effective technique to measure the depth of shower maximum with resolutions competitive to optical detectors. After the decommissioning of Tunka-Rex, as last phase of its lifecycle and following the FAIR (Findability - Accessibility - Interoperability - Reuse) principles, we publish the data and software under free licenses in the frame of the TRVO (Tunka-Rex Virtual Observatory), which is hosted at KIT under the partnership with the KCDC and GRADLCI projects. We present the main features of TRVO, its interface and give an overview of projects, which benefit from its open software and data.

Since 2015, the direct detection of Gravitational Waves (GWs) became possible with ground-based interferometers like LIGO and Virgo. GWs became the center of attention of the astronomical community and electromagnetic observatories took a particular interest in follow-up observations of such events. The main setback of these observations is the poor localization of GW events. In fact, GW localization uncertainties can span tens to hundreds of deg$^{2}$ the sky even with the advanced configurations of current GW interferometers. In this contribution, we present five follow-up strategies developed for the High Energy Stereoscopic System (H.E.S.S.) and assess their performances. We show how a 2D and 3D galaxy targeted search approach exploiting the integral probability inside the instruments field of view are best suited for medium field of view instruments like H.E.S.S. We also develop an automatic response scheme within the H.E.S.S. Transient Follow-up system that is optimized for fast response and is capable of responding promptly to all kind of GW alerts. GW events are filtered by the developed scheme and prompt and afterglow observations are automatically scheduled. The H.E.S.S. response latency to prompt alerts is measured to be less than 1 minute. With this continually optimized GW response scheme, H.E.S.S. scheduled several GW follow-up observations during the second and third LIGO/Virgo observation runs.

The recent progress in the radio detection technique for air showers paves the path to future cosmic-ray radio detectors. Digital radio arrays allow for a measurement of the air-shower energy and depth of its maximum with a resolution comparable to those of the leading optical detection methods. One of the remaining challenges regarding cosmic-ray radio instrumentation is an accurate estimation of their efficiency and aperture. We present a probabilistic model to address this challenge. We use the model to estimate the efficiency and aperture of the Tunka-Rex radio array. The basis of the model is a parametrization of the radio footprint and a probabilistic treatment of the detection process on both the antenna and array levels. In this way, we can estimate the detection efficiency for air showers as function of their arrival direction, energy, and impact point on the ground. In addition, the transparent internal relationships between the different stages of the air-shower detection process in our probabilistic approach enable to estimate the uncertainty of the efficiency and, consequently, of the aperture of radio arrays. The details of the model will be presented in the contribution.

We have searched quasi-periodic oscillations (QPOs) for BL Lac PKS J2134-0153 in the 15 GHz radio light curve announced by the Owens Valley Radio Observatory 40-m telescope during the period from 2008-01-05 to 2019-05-18, utilizing the Lomb-Scargle periodogram (LSP) and the weighted wavelet Z-transform (WWZ) techniques. This is the first time that to search for periodic radio signal in BL Lac PKS J2134-0153 by these two methods. These two methods consistently reveal a QPO of 4.69 $\pm$ 0.14 years (>5 $\sigma$ confidence level). We discuss possible causes for this QPO, and we expected that the binary black holes scenario, where the QPO is caused by the precession of the binary black holes, is the most likely explanation. BL Lac PKS J2134-0153 thus could be a good binary black hole candidate. In the binary black holes scenario, the distance between the primary black hole and the secondary black hole is 1.83$\times$10$^{16}$ cm.

We monitored BL Lacertae in the B, V, R and I bands for 14 nights during the period of 2016-2018. The source showed significant intraday variability on 12 nights. We performed colour-magnitude analysis and found that the source exhibited bluer-when-brighter chromatism. This bluer-when-brighter behavior is at least partly caused by the larger variation amplitude at shorter wavelength. The variations at different wavelengths are well correlated and show no inter-band time lag.

We aim to understand the onset of cloud formation {and study the formation of} TiO$_2$-CCNs. The formation of (TiO$_2$)$_{\rm N}$ clusters as precursors to extrasolar cloud formation is modelled by two different methods in order to understand their potential, identify underlying shortcomings, and to validate our methods. We propose potential spectral tracers for TiO$_2$-CCN formation. We applied three-dimensional Monte Carlo (3D MC) simulations to model the collision-induced growth of TiO$_2$-molecules to (TiO$_2$)$_{\rm N}$-clusters in the free molecular flow regime of an atmospheric gas. We derived individual, time-dependent (TiO$_2$)$_{\rm N}$ cluster number densities. For $T=1000$K, the results are compared to a kinetic approach that utilises thermodynamic data for individual (TiO$_2$)$_{\rm N}$ clusters. The {(TiO$_2$)$_{\rm N}$} cluster size distribution is temperature dependent and evolves in time until a steady state is reached. For $T=1000$K, the 3D MC and the kinetic approach agree well regarding the cluster number densities for $N=1\,\ldots\,10$, the vivid onset of cluster formation, and the long transition into a steady state. Collision-induced growth and evaporation simulated using a 3D MC approach enables a faster onset of cluster growth through nucleation bursts. Different size distributions develop for monomer-cluster and for cluster-cluster growth, with the largest clusters appearing for cluster-cluster growth. The (TiO$_2$)$_N$ cluster growth efficiency has a sweet-spot temperature at $\approx 1000$K at which CCN formation is triggered. The combination of local thermodynamic conditions and chemical processes therefore determines CCN formation efficiency.

The Baikal Gigaton Volume Detector (Baikal-GVD) is a km$^3$-scale neutrino detector currently under construction in Lake Baikal, Russia. The detector currently consists of 2304 optical modules arranged on 64 vertical strings. Further extension of the array is planned for March 2022. The data from the partially complete array have been analyzed using a $\chi^2$-based track reconstruction algorithm. After suppression of the downward-going atmospheric muon background, a flux of upward-going neutrino events is observed, dominated by the atmospheric neutrinos. The observed flux is in good agreement with Monte Carlo predictions.

Blazars are variable targets in the sky, whose variation mechanism remains an open question. In this work, we make a comprehensive study on the variation phenomena of the spectral index and polarization degree (PD) to deeply understand the variation mechanism of B2 1633+382 (4C 38.41). We use the local cross-correlation function (LCCF) to perform the correlation analysis between multi-wavelength light curves. We find that both $\gamma$-ray and optical $V$-band are correlated with the radio 15 GHz at the beyond 3$\sigma$ confidence level. Based on the lag analysis, the emitting regions of $\gamma$-ray and optical locate at $14.2_{-2.4}^{+0}$ pc and $14.2_{-8.3}^{+8.3}$ pc upstream of the core region of radio 15 GHz, and are far away from the broad-line region (BLR). The broad lines in the spectrum indicate the existence of the accretion disk component in the radiation. Thus, we consider the two-component (TC) model, which includes the relative constant background component and the varying jet component to study the variation behaviors. The Markov Chain Monte Carlo (MCMC) procedure is adopted to study the physical parameters of the jet and the background components. To some extent, the study of normalized residuals indicates that the TC model fits better than the linear fitting model. The jet with helical magnetic field is hopeful to explain the variation, and the shock in jet model is not completely ruled out.

We present a detailed analysis of the mid-infrared spectra obtained from the Spitzer Space Telescope of the dark globule, DC 314.8-5.1, which is at the onset of low-mass star formation. The cloud has a serendipitous association with a B-type field star, which illuminates a reflection nebula in the cloud, allowing us to investigate infrared characteristics not otherwise discernible in such systems until a later evolutionary stage. We focus specifically on the polycyclic aromatic hydrocarbon (PAH) emission features prevalent throughout the mid-infrared range. We find that the intensity profiles do not obey any single unique scaling, for example a monotonic decrease related to the decreasing starlight toward the cloud's central regions. We note that a diversity in trends over distance is also present in the intensity profiles of the molecular (H$_2$) and atomic (Ar, Ne, and S) emission lines, which are however much less prominent in the spectrum when compared with the PAH features. All in all, our analysis reveals that (i) there is a stratification in dust sizes within the reflection nebula, with larger grains dominating the PAH emission at the outskirts of the system, and (ii) the ionization level within the reflection nebula is fairly constant, and as such independent on the amount of the ionizing UV continuum from the neighbouring star, (iii) the intensity ratios of the prominent PAH features, do not follow correlations established for the reflection nebulae with active star formation.

The recent discovery and initial characterization of sub-Neptune-sized exoplanets that receive stellar irradiance of approximately Earth's raised the prospect of finding habitable planets in the coming decade, because some of these temperate planets may support liquid water oceans if they do not have massive H2/He envelopes and are thus not too hot at the bottom of the envelopes. For planets larger than Earth, and especially planets in the 1.7-3.5 R_Earth population, the mass of the H2/He envelope is typically not sufficiently constrained to assess the potential habitability. Here we show that the solubility equilibria vs. thermochemistry of carbon and nitrogen gases results in observable discriminators between small H2 atmospheres vs. massive ones, because the condition to form a liquid-water ocean and that to achieve the thermochemical equilibrium are mutually exclusive. The dominant carbon and nitrogen gases are typically CH4 and NH3 due to thermochemical recycling in a massive atmosphere of a temperate planet, and those in a small atmosphere overlying a liquid-water ocean are most likely CO2 and N2, followed by CO and CH4 produced photochemically. NH3 is depleted in the small atmosphere by dissolution into the liquid-water ocean. These gases lead to distinctive features in the planet's transmission spectrum, and a moderate number of repeated transit observations with the James Webb Space Telescope should tell apart a small atmosphere vs. a massive one on planets like K2-18 b. This method thus provides a way to use near-term facilities to constrain the atmospheric mass and habitability of temperate sub-Neptune exoplanets.

Thanks to recent technological development, a new generation of cosmic ray experiments have been developed with more sensitivity to study these particles in the primary energy interval from 10 TeV to 1 PeV, such as HAWC. Due to its design and high altitude, the HAWC gamma-ray and cosmic ray observatory can provide a bridge between the data from direct and indirect cosmic ray detectors. In 2017 the HAWC collaboration published its first result on the total energy spectrum of cosmic rays, which covers the range from 10 to 500 TeV. This work updates the previous result by extending the energy interval of the measured all-particle cosmic-ray energy spectrum up to 1 PeV. The energy spectrum was obtained from the analysis of two years of HAWC's data using an unfolding method. We employed the QGSJET-II-04 model for the energy calibration and the spectrum reconstruction. The results confirm the presence of a knee like feature at tens of TeV, as previously reported by the HAWC collaboration in 2017.

We present the first results of our ongoing project conducting simultaneous multiwavelength observations of flares on nearby active M dwarfs. We acquired data of the nearby dM3.5e star EV Lac using 5 different observatories: NASA's Transiting Exoplanet Survey Satellite (TESS), NASA's Neil Gehrels Swift Observatory (\textit{Swift}), NASA's Neutron Interior Composition Explorer (NICER), the University of Hawaii 2.2-m telescope (UH88) and the Las Cumbres Observatory Global Telescope (LCOGT) Network. During the $\sim$25 days of TESS observations, we acquired three simultaneous UV/X-ray observations using \textit{Swift} that total $\sim$18 ks, 21 simultaneous epochs totaling $\sim$98 ks of X-ray data using NICER, one observation ($\sim$ 3 hours) with UH88, and one observation ($\sim$ 3 hours) with LCOGT. We identified 56 flares in the TESS light curve with estimated energies in the range log $E_{\rm T}$ (erg) = (30.5 - 33.2), nine flares in the \textit{Swift} UVM2 light curve with estimated energies in the range log $E_{UV}$ (erg) = (29.3 - 31.1), 14 flares in the NICER light curve with estimated minimum energies in the range log $E_{N}$ (erg) = (30.5 - 32.3), and 1 flare in the LCOGT light curve with log $E_{L}$ (erg) = 31.6. We find that the flare frequency distributions (FFDs) of TESS and NICER flares have comparable slopes, $\beta_{T}$ = -0.67$\pm$0.09 and $\beta_{N}$ = -0.65$\pm$0.19, and the FFD of UVOT flares has a shallower slope ($\beta_{U}$ = -0.38$\pm$0.13). Furthermore, we do not find conclusive evidence for either the first ionization potential (FIP) or the inverse FIP effect during coronal flares on EV Lac.

We define various types of "phantom" stars that may appear in the TESS Input Catalog (TIC), and provide examples and lists of currently known cases. We present a methodology that can be used to check for phantoms around any object of interest in the TIC, and we present an approach for correcting the TIC-reported flux contamination factors accordingly. We checked all 2077 TESS Objects of Interest (TOIs) known as of July 21st 2020 (Sectors 1 to 24) and sent corrections for 291 stars to MAST where they are integrated into the publicly available TIC-8, updating it to TIC 8.1. We used the experience gained to construct an all-sky algorithm searching for "phantoms" which led to 34 million updates integrated into TIC 8.2.

In order to answer the open questions of modern cosmology and galaxy evolution theory, robust algorithms for calculating photometric redshifts (photo-z) for very large samples of galaxies are needed. Correct estimation of the various photo-z algorithms' performance requires attention to both the performance metrics and the data used for the estimation. In this work, we use the supervised machine learning algorithm MLPQNA to calculate photometric redshifts for the galaxies in the COSMOS2015 catalogue and the unsupervised Self-Organizing Maps (SOM) to determine the reliability of the resulting estimates. We find that for spec-z<1.2, photo-z predictions are on the same level of quality as SED fitting photo-z. We show that the SOM successfully detects unreliable spec-z that cause biases in the estimation of the photo-z algorithms' performance. Additionally, we use SOM to select the objects with reliable photo-z predictions. Our cleaning procedures allow to extract the subset of objects for which the quality of the final photo-z catalogs is improved by a factor of two, compared to the overall statistics.

A scientific instrument comprised of a global network of millions of independent, connected, remote devices presents unique data acquisition challenges. We describe the software design of a mobile application which collects data from smartphone cameras without overburdening the phone's CPU or battery. The deployed software automatically calibrates to heterogeneous hardware targets to improve the quality and manage the rate of data transfer, and connects to a cloud-based data acquisition system which can manage and refine the operation of the network.

High electron mobility transistors are widely used as microwave amplifiers owing to their low microwave noise figure. Electronic noise in these devices is typically modeled by noise sources at the gate and drain. While consensus exists regarding the origin of the gate noise, that of drain noise is a topic of debate. Here, we report a theory of drain noise as a type of partition noise arising from real-space transfer of hot electrons from the channel to the barrier. The theory accounts for the magnitude and dependencies of the drain temperature and suggests strategies to realize devices with lower noise figure.

We propose a novel scenario to explain the small cosmological constant (CC) by a finely tuned inflaton potential. The tuned shape is stable under radiative corrections, and our setup is technically natural. The peculiar po- tential approximately satisfies the following conditions: the inflation is eternal if CC is positive, and not eternal if CC is negative. By introducing a slowly varying CC from a positive value to a negative value, the dominant volume of the Universe after the inflation turns out to have a vanishingly small CC. The scenario does not require eternal inflation but the e-folding number is exponentially large and the inflation scale should be low enough. The scenario can have a consistent thermal history, but the present equation of state of the Universe is predicted to differ from the prediction of the {\Lambda}CDM model. A concrete model with a light scalar field is studied.

Electroweak baryogenesis (EWBG) is sourced by nonstandard $CP$-violating interactions of the Higgs boson with fermions, usually taken to be the top quark, enhanced by its large Yukawa coupling. Numerous papers have studied EWBG sourced by lighter fermions, including the tau lepton and off-diagonal quark mass terms. We critically reassess the viability of EWBG in these scenarios, comparing the predictions based on the semiclassical (WKB) formalism for the source term to those from the VEV insertion approximation (VIA), using updated values for the collision terms, and clarifying discrepancies in the definition of the weak sphaleron rate. The VIA systematically predicts a baryon asymmetry that is orders of magnitude larger than the WKB formalism. We trace this to the differing shapes of the $CP$-violating source terms in the two formalisms, showing that the additional spatial derivative in the WKB source term causes large cancellations when it is integrated over the bubble wall profile. An important exception is a source term from $c$-$t$ quark mixing, where the WKB prediction also allows for a realistically large baryon asymmetry. In contrast, the analogous $b$-$s$ mixing source is found to be orders of magnitude too small.

We study thermal axion production around the confinement scale. At higher temperatures, we extend current calculations to account for the masses of heavy quarks, whereas we quantify production via hadron scattering at lower temperatures. Matching our results between the two opposite regimes provides us with a continuous axion production rate across the QCD phase transition. We employ such a rate to quantify the axion contribution to the effective number of neutrino species.

Convection in the Sun occurs at Rayleigh numbers, $Ra$, as high as $10^{22}$, molecular Prandtl number, $Pr$, as low as $10^{-6}$, and occurs under conditions that are far from satisfying the Oberbeck-Boussinesq (OB) idealization. The effects of these extreme circumstances on turbulent heat transport are unknown, and no comparable conditions exist on Earth. Our goal is to understand how these effects scale (since we cannot yet replicate the Sun's conditions faithfully). We study thermal convection by using direct numerical simulations, and determine the variation with respect to $Pr$, up to $Pr$ as low as $10^{-4}$, of the turbulent Prandtl number, $Pr_t$, which is the ratio of turbulent viscosity to thermal diffusivity. The simulations are primarily two-dimensional but we draw upon some three-dimensional results as well. We focus on non-Oberbeck-Boussinesq (NOB) conditions of a certain type, but also study OB convection for comparison. The OB simulations are performed in a rectangular box of aspect ratio 2 by varying $Pr$ from $O(10)$ to $10^{-4}$ at fixed Grashof number $Gr \equiv Ra/Pr = 10^9$. The NOB simulations are done in the same box by letting only the thermal diffusivity depend on the temperature. Here, the Rayleigh number is fixed at the top boundary while the mean $Pr$ varies in the bulk from 0.07 to $5 \times 10^{-4}$. The three-dimensional simulations are performed in a box of aspect ratio 25 at a fixed Rayleigh number of $10^5$, and $0.005 < Pr < 7$. The principal finding is that $Pr_t$ increases with decreasing $Pr$ in both OB and NOB convection: $Pr_t \sim Pr^{-0.3}$ for OB convection and $Pr_t \sim Pr^{-1}$ for the NOB case. The $Pr_t$-dependence for the NOB case especially suggests that convective flows in the astrophysical settings behave effectively as in high-Prandtl-number turbulence.

We report the first (in)elastic scattering measurement of $^{25}\mathrm{Al}+p$ with the capability to select and measure in a broad energy range the proton resonances in $^{26}$Si contributing to the $^{22}$Mg$(\alpha,p)$ reaction at type I x-ray burst energies. We measured spin-parities of four resonances above the $\alpha$ threshold of $^{26}$Si that are found to strongly impact the $^{22}$Mg$(\alpha,p)$ rate. The new rate advances a state-of-the-art model to remarkably reproduce lightcurves of the GS 1826$-$24 clocked burster with mean deviation $<$9 % and permits us to discover a strong correlation between the He abundance in the accreting envelope of photospheric radius expansion burster and the dominance of $^{22}$Mg$(\alpha,p)$ branch.

In this article, we study the fundamental ($f$-)modes of non-rotating compact stars with realistic equations of state (EoS), extracted in the dynamical spacetime using numerical relativity simulations. We use a set of EoS with varying degree of stiffness and numerically evolve perturbed star models for several mass configurations (in the range of $1.2 - 2.0$ $M_{\odot}$) for each of these EoS. We find the $f$-mode frequency being lower for the stiffer EoS. We also see the increase in the frequency for higher mass systems as compared to the smaller mass cases. While the increase is linear for $M \leq 1.6 M_\odot$, it shows deviation from linearity with the sharper change and higher values of $f$-modes for some of the EoS for more massive systems. We notice that the frequencies are distinguishable for soft, intermediate and stiffer EoS, thus it might be possible to constraint the EoS based on the detected signal frequency from the binary neutron star merger. More specifically for the softer EoSs, the $f$-mode frequency is in the range of $1.8-2.2$ kHz for the masses between $1.2-1.8 M_{\odot}$. On the other hand, in the case of stiffer EoS, such as `BHB' and `DD2' frequency is shifted to lower values $1.55-1.8$ kHz for the same mass range.

We check whether General Relativity's field self-interaction alleviates the need for dark matter to explain the universe's large structure formation. We found that self-interaction accelerates sufficiently the growth of structures so that they can reach their presently observed density. No free parameters, dark components or modifications of the known laws of nature were required. This result adds to the other natural explanations provided by the same approach to the, $inter~alia$, flat rotation curves of galaxies, supernovae observations suggestive of dark energy, and dynamics of galaxy clusters, thereby reinforcing its credibility as an alternative to the dark universe model.

An exact Kerr-like solution has been obtained recently in Einstein-bumblebee gravity model where Lorentz symmetry is spontaneously broken. In this paper, we investigate the superradiance instability of the Kerr-like black hole under the perturbation of a massive scalar field. We find the Lorentz breaking parameter $L$ affects superradiance regime but not the regime of the bound states. We calculate the bound state spectrum via the continued-fraction method and show the influence of $L$ on the maximum binding energy and the damping rate. The superradiance instability could occur since the superradiance condition and the bound state condition could be both satisfied. Compared with Kerr black hole, the nature of the superradiance instability of this black hole depends non-monotonously not only on the rotation speed of the black hole $a$ and the product of the black hole mass $M$ and the field mass $\mu$, but also on the Lorentz breaking parameter $L$. Through the Monte Carlo method, we find that for $l=m=1$ state the most unstable mode occurs at $L=-0.79637$, $a/M=2.213$ and $M\mu=0.439$, with the maximum growth rate of the field $\omega_{I}M=1.676\times10^{-6}$, which is about 10 times of that in Kerr black hole.

A novel fast multi-impulse optimization method for long-duration perturbed orbit rendezvous is proposed. First, based on the analytically estimated impulses, the terminal rendezvous deviation with precise dynamics model can be predicted. Then, an analytical correction to the impulses using the deviations of orbit elements can be calculated based on the analytical J2 perturbed dynamics equation of a circular orbit. The iteration process repeating prediction and correction is then designed to quickly obtain a precise solution and trajectory. The simulation results proved that the iteration method adapts well to the analytical dynamics and high-precision dynamics. The deviation could always converge within five iterations.