Papers for Wednesday, Dec 15 2021

Black hole (BH) shadows can be used to probe new physics in the form of ultra-light particles via the phenomenon of superradiant instability. By directly affecting the BH mass and spin, superradiance can lead to a time evolution of the BH shadow, which nonetheless has been argued to be unobservable through Very Long Baseline Interferometry (VLBI) over realistic observation timescales. We revisit the superradiance-induced BH shadow evolution including the competing effects of gas accretion and gravitational wave (GW) emission and, as a first step towards modelling realistic new physics scenarios which predict the existence of multiple ultra-light species, we study the system in the presence of two ultra-light bosons, whose combined effect could help reducing the shadow evolution timescale. We find that accretion and GW emission play a negligible role in our results (justifying previous simplified analyses), and that contrary to our intuition the inclusion of an additional ultra-light boson does not shorten the BH shadow evolution timescale and hence improve detection prospects. However, we point out an important subtlety concerning the observationally meaningful definition of the superradiance-induced BH shadow evolution timescale, which reduces the latter by about an order of magnitude, opening up the possibility of observing the superradiance-induced BH shadow evolution with upcoming VLBI arrays, provided angular resolutions just below the $\mu{\rm as}$ level can be reached. As a concrete example, we show that the angular size of the shadow of SgrA$^*$ can change by up to $0.6\,\mu{\rm as}$ over a period as short as $16$ years, which further strengthens the scientific case for targeting the shadow of SgrA$^*$ with next-generation VLBI arrays.

Recent measurements from the CMB and from high-redshift galaxy observations have placed rough constraints on the midpoint and duration of the Epoch of Reionization. Detailed measurements of the ionization history remain elusive, although two proposed probes show great promise for this purpose: the 21cm global signal and the kinetic Sunyaev-Zel'dovich (kSZ) effect. We formally confirm the common assumption that these two probes are highly complementary, with the kSZ being more sensitive to extended ionization histories and the global signal to rapidly evolving ones. We do so by performing a Karhunen-Lo\`{e}ve (KL) transformation, which casts the data in a basis designed to emphasize the information content of each probe. We find that reconstructing the ionization history using both probes gives significantly more precise results than individual constraints, although carefully chosen, physically motivated priors play a crucial part in obtaining a bias-free reconstruction. Additionally, in the KL basis, measurements from one probe can be used to detect the presence of residual systematics in the other, providing a safeguard against systematics that would go undetected when data from each probe is analyzed in isolation. Once detected, the modes contaminated by systematics can be discarded from the data analysis to avoid biases in reconstruction.

Over the past decades, several studies have discovered a population of galaxies undergoing very strong star formation events, called extreme emission line galaxies (EELGs). In this work, we exploit the capabilities of the Javalambre Photometric Local Universe Survey (J-PLUS), a wide field multifilter survey, with 2000 square degrees observed. We use it to identify EELGs at low redshift by their [OIII]5007 emission line. We intend to provide with a more complete, deep, and less biased sample of local EELGs. We select objects with an excess of flux in the J-PLUS mediumband $J0515$ filter, which covers the [OIII] line at z$<$0.06. We remove contaminants (stars and higher redshift systems) using J-PLUS and WISE infrared data, with SDSS spectra as a benchmark. We perform spectral energy distribution fitting to estimate the properties of the galaxies: line fluxes, equivalent widths (EWs), masses, etc. We identify 466 EELGs at ${\rm z} < 0.06$ with [OIII] EW over 300 \text{\AA} and $r$-band mag. below 20, of which 411 were previously unknown. Most show compact morphologies, low stellar masses ($\log (M_{\star}/M_{\odot}) \sim {8.13}^{+0.61}_{-0.58}$), low dust extinction ($E(B-V)\sim{0.1}^{+0.2}_{-0.1}$), and very young bursts of star formation (${3.0}^{+2.7}_{-2.0}$ Myr). Our method is up to $\sim$ 20 times more efficient detecting EELGs per Mpc$^3$ than broadband surveys, and as complete as magnitude-limited spectroscopic surveys (and reaching fainter objects). The sample is not directly biased against strong H$\alpha$ emitters, in contrast with broadband surveys. We demonstrate the capability of J-PLUS to identify, following a clear selection process, a large sample of previously unknown EELGs showing unique properties. A fraction of them are likely similar to the first galaxies in the Universe, but at a much lower redshift, which makes them ideal targets for follow-up studies.

We study the streaming instability of GeV-100 GeV cosmic rays (CRs) and its damping in the turbulent interstellar medium (ISM). We find that the damping of streaming instability is dominated by ion-neutral collisional damping in weakly ionized molecular clouds, turbulent damping in the highly ionized warm medium, and nonlinear Landau damping in the Galactic halo. Only in the Galactic halo, is the streaming speed of CRs close to the Alfv\'{e}n speed. Alfv\'{e}nic turbulence plays an important role in both suppressing the streaming instability and regulating the diffusion of streaming CRs via magnetic field line wandering, with the effective mean free path of streaming CRs in the observer frame determined by the Alfv\'{e}nic scale in super-Alfv\'{e}nic turbulence. The resulting diffusion coefficient is sensitive to Alfv\'{e}n Mach number, which has a large range of values in the multi-phase ISM. Super-Alfv\'{e}nic turbulence contributes to additional confinement of streaming CRs, irrespective of the dominant damping mechanism.

Ultra-light axions (ULAs) are a promising and intriguing set of dark-matter candidates. We study the prospects to use forthcoming measurements of 21-cm fluctuations from cosmic dawn to probe ULAs. We focus in particular on the velocity acoustic oscillations (VAOs) in the large-scale 21-cm power spectrum, features imprinted by the long-wavelength ($k\sim0.1\,{\rm Mpc}^{-1}$) modulation, by dark-matter--baryon relative velocities, of the small-scale ($k\sim 10-10^3\, {\rm Mpc}^{-1}$) power required to produce the stars that heat the neutral hydrogen. Damping of small-scale power by ULAs reduces the star-formation rate at cosmic dawn which then leads to a reduced VAO amplitude. Accounting for different assumptions for feedback and foregrounds, experiments like HERA may be sensitive to ULAs with masses up to $m_{\alpha}\approx 10^{-18}\text{eV}$, two decades of mass higher than current constraints.

We examine the origin of dynamical friction using a non-perturbative, orbit-based approach. Unlike the standard perturbative approach, in which dynamical friction arises from the LBK torque due to pure resonances, this alternative, complementary view nicely illustrates how a massive perturber significantly changes the energies and angular momenta of field particles on near-resonant orbits, with friction arising from an imbalance between particles that gain energy and those that lose energy. We treat dynamical friction in a spherical host system as a restricted three-body problem. This treatment is applicable in the `slow' regime, in which the perturber sinks slowly and the standard perturbative framework fails due to the onset of non-linearities. Hence it is especially suited to investigate the origin of core-stalling: the cessation of dynamical friction in central constant-density cores. We identify three different families of near-co-rotation-resonant orbits that dominate the contribution to dynamical friction. Their relative contribution is governed by the Lagrange points (fixed points in the co-rotating frame). In particular, one of the three families, which we call Pac-Man orbits because of their appearance in the co-rotating frame, is unique to cored density distributions. When the perturber reaches a central core, a bifurcation of the Lagrange points drastically changes the orbital make-up, with Pac-Man orbits becoming dominant. In addition, due to relatively small gradients in the distribution function inside a core, the net torque from these Pac-Man orbits becomes positive (enhancing), thereby effectuating a dynamical buoyancy. We argue that core stalling occurs where this buoyancy is balanced by friction.

Observations of TeV emission from early gamma-ray burst (GRB) afterglows revealed the long sought for inverse Compton (IC) upscattering of the lower energy synchrotron. However, it turned out that the long hoped for ability to easily interpret the synchrotron-self-Compton (SSC) spectra didn't materialize. The TeV emission is in the Klein-Nishina (KN) regime and the simple Thomson regime SSC spectrum is modified, complicating the scene. We describe here a methodology, based on Nakar et al. 2009, to obtain an analytic approximation to an observed spectrum and infer the conditions at the emitting region. The methodology is general and can be used in any such source. As a test case, we apply it to the observations of GRB 190114C. We find that the procedure of fitting the model parameters using the analytic SSC spectrum suffers from some generic problems. However, at the same time, it conveniently gives a useful insight into the conditions that shape the spectrum. Once we introduce a correction to the standard KN approximation, the best fit solution is consistent with the one found in detailed numerical simulations. As in the numerical analysis, we find a family of solutions that provide a good approximation to the data and satisfy roughly $B\propto \Gamma^{-3}$ between the magnetic field and the bulk Lorentz factor, and we provide a tentative explanation why such a family arises.

Efficient searches for electromagnetic counterparts to gravitational wave, high-energy neutrino, and gamma-ray burst events, demand rapid processing of image arithmetic and geometry set operations in a database to cross-match galaxy catalogs, observation footprints, and all-sky images. Here we introduce HEALPix Alchemy, an open-source, pure Python implementation of a set of methods that enables rapid all-sky geometry calculations. HEALPix Alchemy is built upon HEALPix, a spatial indexing strategy that is in wide use in astronomical databases and also the native format of LIGO-Virgo-KAGRA gravitational-wave sky localization maps. Our approach leverages new multirange types built into the PostgreSQL 14 database engine. This enables fast all-sky queries against probabilistic multimessenger event localizations and telescope survey footprints. Questions such as "What are the galaxies contained within the 90% credible region of an event?" and "What is the rank-ordered list of the fields within an observing footprint with the highest probability of containing the event?" can be performed in less than a few seconds on commodity hardware using off-the-shelf cloud-managed database implementations without server-side database extensions. Common queries scale roughly linearly with the number of telescope pointings. As the number of fields grows into the hundreds or thousands, HEALPix Alchemy is orders of magnitude faster other implementations. HEALPix Alchemy is now used as the spatial geometry engine within SkyPortal, which forms the basis of the Zwicky Transient Facility transient marshal called Fritz.

Massive stars that become stripped of their hydrogen envelope through binary interaction or winds can be observed either as Wolf-Rayet stars, if they have optically thick winds, or as transparent-wind stripped-envelope stars. We approximate their evolution through evolutionary models of single helium stars, and compute detailed model grids in the initial mass range 1.5 to 70 M$_{\odot}$ for metallicities between 0.01 and 0.04, from core helium ignition until core collapse. Throughout their lifetime, some stellar models expose the ashes of helium burning. We propose that models that have nitrogen-rich envelopes are candidate WN stars, while models with a carbon-rich surface are candidate WC stars during core helium burning, and WO stars afterwards. We measure metallicity dependance of the total lifetime of our models and the duration of their evolutionary phases. We propose an analytic estimate of the wind optical depth to distinguish models of Wolf-Rayet stars from transparent-wind stripped-envelope stars, and find that the luminosity ranges at which WN, WC and WO type stars can exist is a strong function of metallicity. We find that all carbon-rich models produced in our grids have optically thick winds and match the luminosity distribution of observed populations. We construct population models and predict the numbers of transparent-wind stripped-envelope stars and Wolf-Rayet stars, and derive their number ratios at different metallicities. We find that as metallicity increases, the number of transparent-wind stripped-envelope stars decreases and the number of Wolf-Rayet stars increases. At high metallicities WC and WO type stars become more common. We apply our population models to nearby galaxies, and find that populations are more sensitive to the transition luminosity between Wolf-Rayet stars and transparent-wind helium stars than to the metallicity dependent mass loss rates.

The precise nature of Type Ia supernova (SN Ia) progenitors remains a mystery, but the detonation of a sub-Chandrasekhar-mass white dwarf (WD) has become a promising candidate. There is a growing body of work suggesting that the carbon core detonation of a sub-Chandrasekhar-mass WD can be triggered by the detonation of a helium shell accreted from a companion WD, through either inward shock convergence near the center or direct edge-lit detonation. This "double-detonation" SN Ia can be triggered by a small helium shell and is therefore well approximated by the detonation of a bare carbon-oxygen white dwarf (C/O WD). The impacts of uncertainties in experimentally and theoretically determined nuclear reaction rates on nucleosynthesis in the detonations of sub-Chandrasekhar-mass WDs have not yet been fully explored. We investigate the sensitivity of this model to nuclear reaction rate uncertainties to better constrain the nucleosynthetic yields resulting from these phenomena and identify the nuclear reaction rates whose uncertainties have the most significant impacts on nucleosynthesis. We find that the chemical abundances, and in particular those of the iron-group elements, are relatively insensitive to these nuclear reaction rate uncertainties.

We examine the contribution of high-redshift (z>6) active galactic nuclei (AGNs) to cosmic hydrogen reionization, by tracing the growth and ionizing output of the first generation of supermassive black holes (SMBHs). Our calculations are anchored to the observed population of z~6 quasars, and trace back the evolving spectral energy distributions (SEDs) of the accretion flows that power these early AGNs and consider a variety of growth histories, including super-Eddington accretion. Compared to a fixed-shape SED, the evolving thin disks can produce ionizing radiation that is higher by up to ~80%. Across a variety of SMBH growth scenarios, the contribution of AGNs to reionization is limited to late epochs (z<7), and remains sub-dominant compared to star-forming galaxies. This conclusion holds irrespective of the (still unknown) space density of low-luminosity z=6 AGNs, and for growth scenarios that allow super-Eddington accretion. The contribution of AGNs to reionization can extend to earlier epochs (z>8) in scenarios with relative slow SMBH mass growth, i.e., for low accretion rates and/or high spins. We finally demonstrate that our framework can reproduce the observed quasar proximity zone sizes, and that compact proximity zones around z=6 quasars can be explained by the late onset of super-Eddington accretion.

We present the result of a detailed analysis of HST UV and optical deep images of the massive and young (~1.5 Gyr) stellar cluster NGC 1783 in the Large Magellanic Cloud. This system does not show evidences of multiple populations (MPs) along the red giant branch (RGB) stars. However, we find that the cluster main-sequence (MS) shows evidence of a significant broadening (50% larger than what is expected from photometric errors) along with hints of possible bimodality in the MP sensitive (F343N - F438W, F438W) color-magnitude-diagram (CMD). Such an effect is observed in all color combinations including the F343N filter, while it is not found in the optical CMDs. This observational evidence suggests we might have found light-element chemical abundance variations along the MS of NGC 1783, which represents the first detection of MPs in a system younger than 2 Gyr. A comparison with isochrones including MP-like abundances shows that the observed broadening is compatible with a N abundance enhancement of Delta([N/Fe]) ~0.3. Our analysis also confirms previous results about the lack of MPs along the cluster RGB. However, we find that the apparent disagreement between the results found on the MS and the RGB is compatible with the mixing effects linked to the first dredge-up. This study provides new key information about the MP phenomenon and suggests that star clusters form in a similar way at any cosmic age.

We present an absolute calibration of the J-region Asymptotic Giant Branch (JAGB) method using published photometry of resolved stars in 20 nearby galaxies observed with HST using the WFC3-IR camera and the F110W (Broad J-Band) filter. True distance moduli for each of the galaxies are based on the Tip of the Red Giant Branch (TRGB) method as uniformly determined by Dalcanton et al. (2012). From a composite color-magnitude diagram composed of over 6 million stars, leading to a sample of 453 JAGB stars in these galaxies, we find M_{F110W}{JAGB} = -5.77 +/- 0.02 mag(statistical error on the mean). The external scatter seen in a comparison of the individual TRGB and the JAGB moduli is +/-0.081 mag (or 4% in distance). Some of this scatter can be attributed to small-number statistics arising from the sparse JAGB populations found in the generally low-luminosity galaxies that comprise the particular sample studied here. However, if this inter-method scatter is shared equitably between the JAGB and TRGB methods that implies that each are good to +/-0.06 mag, or better than 3% in distance.

We present unprecedented high fidelity radio images of the M87 jet. We analyzed Jansky Very Large Array (VLA) broadband, full polarization, radio data from 4 to 18 GHz. The observations were taken with the most extended configuration (A configuration), which allow the study of the emission of the jet up to kpc scales with a linear resolution $\sim$10 pc. The high sensitivity and resolution of our data allow to resolve the jet width. We confirm a double-helix morphology of the jet material between $\sim$300 pc and $\sim$1 kpc. We found a gradient of the polarization degree with a minimum at the projected axis and maxima at the jet edges, and a gradient in the Faraday depth with opposite signs at the jet edges. We also found that the behavior of the polarization properties along the wide range of frequencies is consistent with internal Faraday depolarization. All these characteristics strongly support the presence of a helical magnetic field in the M87 jet up to 1 kpc from the central black hole although the jet is most likely particle dominated at these large scales. Therefore, we propose a plausible scenario in which the helical configuration of the magnetic field has been maintained to large scales thanks to the presence of Kelvin-Helmholtz instabilities.

A new time--distance far-side imaging technique was recently developed by utilizing multiple multi-skip acoustic waves. The measurement procedure is applied to 11 years of Doppler observations from the Solar Dynamics Observatory / Helioseismic and Magnetic Imager, and over 8000 far-side images of the Sun have been obtained with a 12-hour temporal cadence. The mean travel-time shifts in these images unsurprisingly vary with the solar cycle. However, the temporal variation does not show good correlations with the magnetic activity in their respective northern or southern hemisphere, but show very good anti-correlation with the global-scale magnetic activity. We investigate four possible causes of this travel-time variation. Our analysis demonstrates that the acoustic waves that are used for mapping the Sun's far side experience surface reflections around the globe, where they may interact with surface or near-surface magnetic field, and carry travel-time deficits with them. The mean far-side travel-time shifts from these acoustic waves therefore vary in phase with the Sun's magnetic activity.

Empirically, the estimated lifetime of a typical protoplanetary disk is $<5-10$ Myr. However, the disk lifetimes required to produce a variety of observed exoplanetary systems may exceed this timescale. Some hypothesize that this inconsistency is due to estimating disk fractions at the cores of clusters, where radiation fields external to a star-disk system can photoevaporate the disk. To test this, we have observed a field on the western outskirts of the IC 1396 star-forming region with \textit{XMM-Newton} to identify new Class III YSO cluster members. Our X-ray sample is complete for YSOs down to $1.8\,M_{\odot}$. We use a subset of these X-ray sources that have near- and mid-infrared counterparts to determine the disk fraction for this field. We find that the fraction of X-ray-detected cluster members that host disks in the field we observe is $17_{-7}^{+10}\%$ (1$\sigma$), comparable with the $29_{-3}^{+4}\%$ found in an adjacent field centered on the cometary globule IC 1396A. We re-evaluate YSO identifications in the IC 1396A field using \textit{Gaia} parallaxes compared to previous color-cut-only identifications, finding that incorporating independent distance measurements provides key additional constraints. Given the existence of at least one massive star producing an external radiation field in the cluster core, the lack of statistically significant difference in disk fraction in each observed field suggests that disk lifetimes remain consistent as a function of distance from the cluster core.

By means of three-dimensional hydrodynamical simulations, we investigate the formation of second generation (SG) stars in young globular clusters of different masses. We consider clusters with a first generation of asymptotic giant branch (AGB) stars with mass 10^5 and 10^6 Msun moving at constant velocity through a uniform gas with density 10^(-24) and 10^(-23) g cm^(-3). Our setup is designed to reproduce the encounter of a young cluster with a reservoir of dense gas, e. g. during its orbital motion in the host galaxy. In the low-density models, as a result of the cooling AGB ejecta which collect in the centre, weakly perturbed by the external ram pressure, a compact central He-rich SG stellar component is formed on a timescale which decreases with increasing initial cluster mass. Our high-density models are subject to stronger ram pressure, which prevents the accumulation of the most He-rich AGB ejecta in the cluster centre. As a result, the SG is more extended and less He-enhanced than in the low-density models. By combining our results with previous simulations, we are able to study relevant, cluster-related scaling relations across a dynamical range of two orders of magnitude in mass (from 10^5 Msun to 10^7 Msun). In agreement with current observationally-based estimates, we find positive correlations between the SG-to-total number ratio and maximum He enhancement in SG stars as a function of the initial cluster mass.

We present the package Gravelamps which is designed to analyse lensed gravitational wave signals in order to constrain the mass density profile of the lensing object. Gravelamps does this via parameter estimation using the framework of bilby, which enables estimation of both the lens and the source parameters. The package can be used to study both microlensing and macrolensing cases -- where the lensing mass distribution is described by a point mass and extended mass density profile respectively -- and allows the user to easily and freely switch between the full wave optics and approximate geometric optics descriptions, which are applied to each of these cases. The performance of Gravelamps is demonstrated via simulated analysis of both a microlensing and macrolensing event, illustrating its capability for both parameter estimation and model selection. To further demonstrate the utility of the package, the real gravitational-wave event GW170809 was analysed using Gravelamps; this event was found to yield no strong evidence supporting the lensing hypothesis, consistent with previously published results.

Magnetars are the most promising progenitors of Fast Radio Bursts (FRBs). Strong FRB radio waves propagating through the magnetar wind are subject to non-linear effects, including modulation/filamentation instabilities. In pair plasmas, spatial modulations of the wave intensity grow exponentially because the ponderomotive force pushes particles outside regions of enhanced radiation intensity. Then in these regions the plasma refractive index increases, which creates a converging lens that further enhances the radiation intensity. Since in magnetically-dominated plasmas the ponderomotive force pushes particles primarily along the field lines, the FRB radiation intensity develops sheets perpendicular to the direction of the wind magnetic field. The radiation sheets are eventually scattered due to diffraction. The FRB scattering timescale depends on the properties of the magnetar wind. In a cold wind, the typical scattering timescale is $\tau_{\rm sc}\sim{\rm\; \mu s-ms}$ at the frequency $\nu\sim 1{\rm\; GHz}$. The scattering timescale increases at low frequencies, with the scaling $\tau_{\rm sc}\propto\nu^{-2}$. The frequency-dependent broadening of the brightest pulse of FRB 181112 is consistent with this scaling. From the scattering timescale of the pulse, one can estimate that the wind Lorentz factor is larger than a few tens. In a warm wind, the scattering timescale can approach $\tau_{\rm sc}\sim{\rm\; ns}$. Then scattering produces a frequency modulation of the observed intensity with a large bandwidth, $\Delta\nu\sim 1/\tau_{\rm sc}\gtrsim 100{\rm\; MHz}$. Broadband frequency modulations observed in FRBs could be due to scattering in a warm magnetar wind.

We present deep Hubble Space Telescope near-infrared (NIR) observations of the magnetar SGR 1935+2154 from June 2021, approximately 6 years after the first HST observations, a year after the discovery of fast radio burst like emission from the source, and in a period of exceptional high frequency activity. Although not directly taken during a bursting period the counterpart is a factor of ~1.5 to 2.5 brighter than seen at previous epochs with F140W(AB) = $24.65\pm0.02$ mag. We do not detect significant variations of the NIR counterpart within the course of any one orbit (i.e. on minutes-hour timescales), and contemporaneous X-ray observations show SGR 1935+2154 to be at the quiescent level. With a time baseline of 6 years from the first identification of the counter-part we place stringent limits on the proper motion of the source, with a measured proper motion of ${\mu} = 3.1\pm1.5$ mas/yr. The direction of proper motion indicates an origin of SGR 1935+2154 very close to the geometric centre of SNR G57.2+08, further strengthening their association. At an adopted distance of $6.6\pm0.7$ kpc, the corresponding tangential space velocity is ${\nu_T} = 97\pm48$ km/s (corrected for differential Galactic rotation and peculiar Solar motion), although its formal statistical determination may be compromised owing to few epochs of observation. The current velocity estimate places it at the low end of the kick distribution for pulsars, and makes it among the lowest known magnetar kicks. When collating the few-magnetar kick constraints available, we find full consistency between the magnetar kick distribution and the much larger pulsar kick sample

We study the striking case of a blue narrow stream with a possible globular cluster-like progenitor around the Milky Way-size galaxy NGC 7241 and its foreground dwarf companion. We present a follow-up spectroscopic study of this stream based on data taken with the MEGARA instrument at the 10.4-m Gran Telescopio Canarias using the integral field spectroscopy mode. Although our data suggest that this compact object in the stream is actually a foreground Milky Way halo star, we detect emission lines overlapping a less compact, bluer and fainter blob of the stream that is clearly visible in both ultra-violet and optical deep images. From its heliocentric systemic radial velocity derived from the [OIII] 5007A lines (V_syst= 1548.58+/-1.80 km\s^-1) and new UV and optical broad-band photometry, we conclude that this over-density could be the actual core of the stream, with an absolute magnitude of Mg~ -10 and a g-r = 0.08+/- 0.11, consistent with a remnant of a low-mass dwarf satellite undergoing a current episode of star formation. From the width of the stream, we calculate that the progenitor mass is between 6.4 x 10^6 Mo -2.7 x 10^7 Mo, which is typical of a dwarf galaxy. These estimates suggest that this is one of the lowest mass streams detected so far beyond the Local Group. We find that blue stellar streams containing star formation regions are commonly predicted by high-resolution cosmological simulations of galaxies lighter than the Milky Way. This scenario is consistent with the processes explaining the bursty star formation history of some dwarf satellites, which are followed by a gas depletion and a fast quenching once they enter within the virial radius of their host galaxies. Thus, it is likely that the stream's progenitor is suffering a star-formation burst comparable to those that have shaped the star-formation history of several Local Group dwarfs in the last few Gigayears.

Observed chemical species in the Venusian mesosphere show local-time variabilities. SO2 at the cloud top exhibits two local maxima over local time, H2O at the cloud top is uniformly distributed, and CO in the upper atmosphere shows a statistical difference between the two terminators. In this study, we investigated these local-time variabilities using a three-dimensional (3D) general circulation model (GCM) in combination with a two-dimensional (2D) chemical transport model (CTM). Our simulation results agree with the observed local-time patterns of SO2, H2O, and CO. The two-maximum pattern of SO2 at the cloud top is caused by the superposition of the semidiurnal thermal tide and the retrograde superrotating zonal (RSZ) flow. SO2 above 85 km shows a large day-night difference resulting from both photochemistry and the sub-solar to anti-solar (SS-AS) circulation. The transition from the RSZ flows to SS-AS circulation can explain the CO difference between two terminators and the displacement of the CO local-time maximum with respect to the anti-solar point. H2O is long-lived and exhibits very uniform distribution over space. We also present the local-time variations of HCl, ClO, OCS and SO simulated by our model and compare to the sparse observations of these species. This study highlights the importance of multidimensional CTMs for understanding the interaction between chemistry and dynamics in the Venusian mesosphere.

Tidal dissipation is thought to be responsible for the observed high heat loss on Enceladus. Forced librations can enhance tidal dissipation in the ice shell, but how such librations affect the thermal state of Enceladus has not been investigated. Here we investigate the heating effect of forced librations using the model of Van Hoolst et al. (2013), which includes the elasticity of the ice shell. We find that libration heating in the ice shell is insufficient to match the inferred conductive heat loss of Enceladus. This suggests that either Enceladus is not in a thermal steady state, or additional heating mechanisms beneath the ice shell are contributing the bulk of the power. In the presence of such an additional heat source, Enceladus resides in a stable thermal equilibrium, resisting small perturbations to the shell thickness. Our results do not support the occurrence of a runaway melting process proposed by Luan and Goldreich (2017). In our study, the strong dependence of conductive loss on shell thickness stabilizes the thermal state of Enceladus's ice shell. Our study implies that thermal runaway (if it occurred) or episodic heating on Enceladus is unlikely to originate from librations of the ice shell.

We present a narrowband [O III] imaging survey of 111 AGN hosts and 17 merging-galaxy systems, in search of distant extended emission-line regions (EELRs) around AGN (either extant or faded). Our data reach deeper than detection from the broadband SDSS data, and cover a wider field than some early emission-line surveys used to study extended structure around AGN. Spectroscopic followup confirms two new distant AGN-ionized clouds, in the merging systems NGC 235 and NGC 5514, projected at 26 and 75 kpc from the nuclei (respectively). We also recover the previously-known region in NGC 7252. These results strengthen the connection between EELRs and tidal features; kinematically quiescent distant EELRs are virtually always photoionized tidal debris. We see them in ~10% of the galaxies in our sample with tidal tails. Energy budgets suggest that the AGN in NGC 5514 has faded by >3 times during the extra light-travel time ~250,000 years from the nucleus to the cloud and then to the observer; strong shock emission in outflows masks the optical signature of the AGN. For NGC 235 our data are consistent with but do not unequivocally require variation over ~85,000 years. In addition to these very distant ionized clouds, luminous and extensive line emission within four galaxies - IC 1481, ESO 362-G08, NGC 5514, and NGC 7679. IC 1481 shows apparent ionization cones, a rare combination with its LINER AGN spectrum. In NGC 5514, we measure a 7-kpc shell expanding at ~370 km/s west of the nucleus.

We review the Seyfert 1.5 Galaxy ESO 362-G18 for exploring the origin of the soft X-ray excess. The Warm Corona and Relativistic Reflection models are two main scenarios to interpret the soft X-ray excess in AGNs at present. We use the simultaneous X-ray observation data of XMM-Newton and NuSTAR on Sep. 24th, 2016 to perform spectral analysis in two steps. First, we analyze the time-average spectra by using Warm Corona and Relativistic Reflection models. Moreover, we also explore the Hybrid model, Double Reflection model and Double Warm Corona model. We find that both of Warm Corona and Relativistic Reflection models can interpret the time-average spectra well but cannot be distinguished easily based on the time-averaged spectra fit statistics. Second, we add the RMS and covariance spectra to perform the spectral analysis with time-average spectra. The result shows that the warm corona could reproduce all of these spectra well. The the hot, optical thin corona and neutral distant reflection will increase their contribution with the temporal frequency, meaning that the corona responsible for X-ray continuum comes from the inner compact X-ray region and the neutral distant reflection is made of some moderate scale neutral clumps.

In this study, we propose a three-stage training approach of neural networks for both photometric redshift estimation of galaxies and detection of out-of-distribution (OOD) objects. Our approach comprises supervised and unsupervised learning, which enables using unlabeled (UL) data for OOD detection in training the networks. Employing the UL data, which is the dataset most similar to the real-world data, ensures a reliable usage of the trained model in practice. We quantitatively assess the model performance of photometric redshift estimation and OOD detection using in-distribution (ID) galaxies and labeled OOD (LOOD) samples such as stars and quasars. Our model successfully produces photometric redshifts matched with spectroscopic redshifts for the ID samples and identifies well the LOOD objects with more than 98% accuracy. Although quantitative assessment with the UL samples is impracticable due to the lack of labels and spectroscopic redshifts, we also find that our model successfully estimates reasonable photometric redshifts for ID-like UL samples and filter OOD-like UL objects. The code for the model implementation is available at https://github.com/GooLee0123/MBRNN_OOD.

Clusters of galaxies are the largest virialized objects in the Universe. Because mergers of these objects are the most energetic events in the Universe--driving shocks and turbulence that heat the gas and accelerate non-thermal phenomena in ways that are yet to be understood--they are inherently interesting. The galaxy cluster Abell 3395 is at an early stage merger containing two subclusters and is also connected to Abell 3391 via an intercluster filament. In this paper, we study the connection between Abell 3395 and the intercluster filament with NuSTAR, XMM-Newton, and Suzaku data. Since the NuSTAR observation is moderately contaminated by scattered light, we present a novel technique developed for disentangling this background from the intracluster medium emission. We find that the interface of the cluster and the intercluster filament does not show any signs of heated plasma, as was previously thought. This interface has low temperature, high density and low entropy, thus we suggest that the gas is cooling, being enhanced by the turbulent or tidal 'weather' driven during the early stage of the merger. Furthermore, our temperature results from the NuSTAR data are in agreement with that of XMM-Newton, and joint NuSTAR and XMM-Newton analysis for a region with ~25% scattered light contamination within 1 sigma. We show that the temperature constraint of the intracluster medium is valid even when the data is contaminated up to ~25% for ~5 keV cluster emission.

System equivalent flux density (SEFD) is an important figure of merit of a radio telescope. This paper aims to derive a general expression for SEFD of a polarimetric tripole interferometer. The derivation makes only two basic and reasonable assumptions. First, the noise under consideration is zero mean and when expressed in complex phasor domain, has independent and identically distributed (iid) real and imaginary components. Correlated and non-identically distributed noise sources are allowed as long as the real and imaginary components remain iid. Second, the system noise is uncorrelated between the elements separated by a baseline distance. The SEFD expression is derived from first principles, that is the standard deviation of the noisy flux estimate in a target direction due to system noise. The resulting SEFD expression is expressed as a simple matrix operation that involves a mixture of the system temperatures of each antenna and the Jones matrix elements. It is not limited to tripoles, but rather, fully extensible to multipole antennas; it is not limited to mutually orthogonal antennas. To illustrate the usefulness of the expression and how the formula is applied, we discuss an example calculation based on a tripole interferometer on lunar orbit for ultra-long wavelengths observation. We compared the SEFD results based on a formula assuming short dipoles and the general expression. As expected, the SEFDs converge at the ultra-long wavelengths where the dipoles are well-approximated as short dipoles. The general SEFD expression can be applied to any multipole antenna systems with arbitrary shapes.

While stellar rotation periods $P_\mathrm{rot}$ may be measured from broadband photometry, the photometric modulation becomes harder to detect for slower rotators, which could bias measurements of the long-period tail of the $P_\mathrm{rot}$ distribution. Alternatively, the $P_\mathrm{rot}$ distribution of stars can be inferred from their projected rotation velocities $v\sin i$ and radii $R$, without being biased against photometrically quiet stars. We solve this inference problem using a hierarchical Bayesian framework, which (i) is applicable to heteroscedastic measurements of $v\sin i$ and $R$ with non-Gaussian uncertainties and (ii) does not require a simple parametric form for the true $P_\mathrm{rot}$ distribution. We test the method on simulated data sets and show that the true $P_\mathrm{rot}$ distribution can be recovered from $\gtrsim 100$ sets of $v\sin i$ and $R$ measured with precisions of $1\,\mathrm{km/s}$ and $4\%$, respectively, unless the true distribution includes sharp discontinuities. We apply the method to a sample of 144 late-F/early-G dwarfs in the Kepler field with $v\sin i$ measured from Keck/HIRES spectra, and find that the typical rotation periods of these stars are similar to the photometric periods measured from Kepler light curves: we do not find a large population of slow rotators that are missed in the photometric sample, although we find evidence that the photometric sample is biased for young, rapidly-rotating stars. Our results also agree with asteroseismic measurements of $P_\mathrm{rot}$ for Kepler stars with similar ages and effective temperatures, and show that $\approx 1.1\,M_\odot$ stars beyond the middle of their main-sequence lifetimes rotate faster than predicted by standard magnetic braking laws.

We report the detection of radial velocity variations in nine evolved G- and K-type giant stars. The observations were conducted at Okayama Astrophysical Observatory. Planets or planet candidates can best explain these regular variations. However, a coincidence of near 280-day variability among five of them prevents us from fully ruling out stellar origins for some of the variations, since all nine stars behave similarly in stellar properties. In the planet hypotheses to the RV variations, the planets (including one candidate) may survive close to the boundary of the so-called "planet desert" around evolved stars, having orbital periods between 255 and 555 days. Besides, they are the least-massive giant planets detected around G- and K-type giant stars, with minimum masses between 0.45$M_{\rm{J}}$ and 1.34$M_{\rm{J}}$. We further investigated other hypotheses for our detection, yet none of them can better explain regular RV variation. With our detection, it is convinced that year-long regular variation with amplitude down to 15 $\rm{m\ s^{-1}}$ for G- and K-type giant stars is detectable. Moreover, we performed simulations to further confirm the detectability of planets around these stars. Finally, we explored giant planets around intermediate-mass stars, and likewise found a 4 Jupiter mass gap (e.g. Santo et al. 2017), which is probably a boundary of the giant planet population.

The sky at MeV energies is currently poorly explored. Here we present an innovative mission concept and we outline the scientific motivation for combining a coded mask and a Compton telescope. The Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO) is a novel concept for a next-generation telescope covering hard X-ray and soft gamma-ray energies. The potential and importance of this approach that will bridge the observational gap in the MeV energy range are presented. With the unprecedented angular resolution of the coded mask telescope combined with the sensitive Compton telescope, a mission such as GECCO will finally disentangle the discrete sources from the truly diffuse emission, unveiling the origin of the gamma-ray Galactic center excess and the Fermi Bubbles, and uncovering properties of low-energy cosmic rays, and their propagation in the Galaxy. Individual Galactic and extragalactic sources will be detected, which will also allow studies of source populations. Nuclear and annihilation lines will be spatially and spectrally resolved from the continuum emission and from sources, addressing the role of low-energy cosmic rays in star formation and galaxy evolution, the origin of the 511 keV positron line, fundamental physics, and the chemical enrichment in the Galaxy. It will also detect explosive transient gamma-ray sources, which will enable identifying and studying the astrophysical objects that produce gravitational waves and neutrinos in a multi-messenger context. A GECCO mission will provide essential contributions to the areas "New Messengers and New Physics" and "Unveiling the Drivers of Galaxy Growth" emphasized in the Decadal Report on Astronomy and Astrophysics 2020.

The occurrence of dual active galactic nuclei (AGN) on scales of a few tens of kpc can be used to study merger-induced accretion on massive black holes (MBHs) and to derive clues on MBH mergers, using dual AGN as a parent population of precursors. We investigate the properties of dual AGN in the cosmological simulation Horizon-AGN. We create catalogs of dual AGN selected with distance and luminosity criteria, plus sub-catalogs where further mass cuts are applied. We divide the sample into dual AGN hosted in different galaxies, on the way to a merger, and into those hosed in one galaxy, after the galaxy merger has happened. We find that the most luminous AGN in a pair has higher Eddington ratio and mass than the general AGN population, but that the relation between MBH and galaxy mass is similar to that of all AGN. The typical mass ratio of galaxy mergers associated to dual AGN is 0.2, with mass loss in the smaller galaxy decreasing the mass ratio as the merger progresses. The connection between dual AGN and MBH mergers is weaker. Between 30 and 80 per cent of dual AGN with separations between 4 and 30 kpc can be matched to an ensuing MBH merger. The dual AGN fraction increases with redshift and with separation threshold, although above 50 kpc the increase of multiple AGN limits that of duals. Multiple AGN are generally associated with massive halos, and mass loss of satellites shapes the galaxy-halo relation.

The Quark Nugget (QN) model of dark matter suggests that the dark matter may consist of compact composite objects of quark matter. Although such composite particles can strongly interact with visible matter, they may remain undetected because of a small cross section to mass ratio. We focus on anti-QNs made of antiquarks since they are heated by annihilation with visible matter and radiate. We study the radiation spectrum and power from anti-QNs in our galaxy and compare them with satellite observations. Thermal radiation from anti-QN is produced by fluctuations of the positron density. We calculate the thermal radiation of anti-QNs with the use of the Mie theory and found its ratio to the black-body radiation. This allows us to find the equilibrium temperature of anti-QNs in the interstellar medium and determine their contribution to the observed diffuse background radiation in our galaxy in different frequency intervals, from radio to UV. We also consider non-thermal radiations from anti-QNs which are produced by products of annihilations of particles of the interstellar gas with anti-QNs. Such radiations include photons from decays of $\pi^0$ mesons, synchrotron, bremsstrahlung and transition radiations from $\pi^\pm$ mesons, electrons and positrons. Synchrotron radiation in MHz frequency range and flux of photons from $\pi^0$ decays may be above the detection threshold in such detectors as Fermi-LAT.

We propose a random forest (RF) machine learning approach to determine the accreted stellar mass fractions ($f_\mathrm{acc}$) of central galaxies, based on various dark matter halo and galaxy features. The RF is trained and tested using 2,710 galaxies with stellar mass $\log_{10}M_\ast/M_\odot>10.16$ from the TNG100 simulation. For galaxies with $\log_{10}M_\ast/M_\odot>10.6$, global features such as halo mass, size and stellar mass are more important in determining $f_\mathrm{acc}$, whereas for galaxies with $\log_{10}M_\ast/M_\odot \leqslant 10.6$, features related to merger histories have higher predictive power. Galaxy size is the most important when calculated in 3-dimensions, which becomes less important after accounting for observational effects. In contrast, the stellar age, galaxy colour and star formation rate carry very limited information about $f_\mathrm{acc}$. When an entire set of halo and galaxy features are used, the prediction is almost unbiased, with root-mean-square error (RMSE) of $\sim$0.068. If only using observable features, the RMSE increases to $\sim$0.104. Nevertheless, compared with the case when only stellar mass is used, the inclusion of other observable features does help to decrease the RMSE by $\sim$20%. Lastly, when using galaxy density, velocity and velocity dispersion profiles as features, which represent approximately the maximum amount of information one can extract from galaxy images and velocity maps, the prediction is only slightly improved. Hence, with observable features, the limiting precision of predicting $f_\mathrm{acc}$ is $\sim$0.1, and the multi-component decomposition of galaxy images should have similar or larger uncertainties. If the central black hole mass and the spin parameter of galaxies can be accurately measured in future observations, the RMSE is promising to be further decreased by $\sim$20%.

For the subclass of Horndeski theory of gravity, we investigate the effects of reheating on the predictions of natural inflation. In the presence of derivative self-interaction of a scalar field and its kinetic coupling to the Einstein tensor, the gravitational friction to inflaton dynamics is enhanced. As a result, the tensor-to-scalar ratio $r$ is suppressed. We place the observational constraints on a natural inflation model and show that the model is now consistent with the observational data for some plausible range of the model parameter $\Delta$, mainly due to the suppressed tensor-to-scalar ratio. To be consistent with the data at the $1\sigma$ ($68\%$ confidence) level, a slightly longer $N_k\gtrsim60$ duration of inflation than usually assumed is preferred. Since the duration of inflation, for any specific inflaton potential, is related to reheating parameters, including the duration $N_{re}$, temperature $T_{re}$, and equation-of-state $\omega_{re}$ parameter during reheating, we imposed the effects of reheating to the inflationary predictions to put further constraints. The results show that the duration of inflation $N_k$ is affected by considerations of reheating, mainly by the $\omega_{re}$ and $T_{re}$ parameters. If reheating occurs instantaneously for which $N_{re}=0$ and $\omega_{re}=1/3$, the duration of inflation is estimated to be $N_k\simeq57$, where the exact value is less sensitive to the model parameter $\Delta$ compatible with the CMB data. The duration of inflation is longer (or shorter) than $N_k\simeq57$ for the equation of state larger (or smaller) than 1/3 hence $N_{re}\neq0$. The maximum temperature at the end of reheating is $T_{re}^\text{max}\simeq3\times 10^{15}$ GeV, which corresponds to the instantaneous reheating. The low reheating temperature, as low as a few MeV, is also possible when $\omega_{re}$ is closer to $1/3$.

We study the effect of Bethe-Heitler (BeHe) pair production on a proton synchrotron model for the prompt emission in gamma-ray bursts (GRBs). The possible parameter space of the model is constrained by consideration of the synchrotron radiation from the secondary BeHe pairs. We find two regimes of interest. 1) At high bulk Lorentz factor, large radius and low luminosity, proton synchrotron emission dominates and produces a spectrum in agreement with observations. For part of this parameter space, a subdominant (in the MeV band) power-law is created by the synchrotron emission of the BeHe pairs. This power-law extends up to few tens or hundreds of MeV. Such a signature is a natural expectation in a proton synchrotron model, and it is seen in some GRBs, including GRB 190114C recently observed by the MAGIC observatory. 2) At low bulk Lorentz factor, small radius and high luminosity, BeHe cooling dominates. The spectrum achieves the shape of a single power-law with spectral index $\alpha = -3/2$ extending across the entire GBM/Swift energy window, incompatible with observations. Our theoretical results can be used to further constrain the spectral analysis of GRBs in the guise of proton synchrotron models.

"Changing-look" active galactic nuclei (CL-AGNs) are a newly-discovered class of AGNs that show the appearance (or disappearance) of broad emission lines within a short time scale (months to years) and are often associated with the dramatic change of their continuum emissions. They provide us an unprecedented chance to directly investigate the host galaxy properties with minimal contamination from the luminous central engine during the "turn-off" state, which is difficult for normal luminous AGNs. In this work, for the first time, we systematically characterize the stellar populations and star formation histories (SFHs) of host galaxies for 26 turn-off CL-AGNs using the stellar population synthesis code STARLIGHT. We find that the stellar populations of CL-AGNs are similar to that of normal AGNs, excepts that the intermediate stellar populations contribute more fraction. We estimate their stellar velocity dispersions ($\rm \sigma_{\star}$) and black hole masses ($\rm M_{BH,vir}$) and find that CL-AGNs also follow the overall $\rm M_{BH}-\sigma_{\star}$ relationship. We also confirm the previous claim that CL-AGNs tend to be biased towards lower Eddington ratio, and their extreme variabilities are more likely due to the intrinsic changes of accretion rates. In addition, CL-AGNs with recent star formations (SF) tend to have higher Eddington ratio. Compared with previous studies, our analysis suggests that there may be a correlation between the CL-AGN host galaxy properties and their CL phenomena.

Millilensing of Gamma-Ray Bursts (GRBs) is expected to manifest as multiple emission episodes in a single triggered GRB with similar light-curve patterns and similar spectrum properties. Identifying such lensed GRBs could help improve constraints on the abundance of compact dark matter. Here we present a systemic search for millilensing among 3000 GRBs observed by the \textit{Fermi} GBM up to 2021 April. Eventually we find 4 interesting candidates by performing auto-correlation test, hardness test, and time-integrated/resolved spectrum test to the whole sample. GRB 081126A and GRB 090717A are ranked as the first class candidate based on their excellent performance both in temporal and spectrum analysis. GRB 081122A and GRB 110517B are ranked as the second class candidates (suspected candidates), mainly because their two emission episodes show clear deviations in part of the time-resolved spectrum or in the time-integrated spectrum. Considering a point mass model for the gravitational lens, our results suggest that the density parameter of lens objects with mass $M_{\rm L}\sim10^{6} M_{\odot}$ is larger than $1.5\times10^{-3}$.

The motion of Mercury using numerical methods in the framework of a model including only the non-relativistic Newtonian gravitational interactions of the solar system, 9 planets in translation (including Pluto) around the sun has been studied. Since the true trajectory of Mercury is an open, non-planar curve, we have paid special attention to the exact definition of the advance of Mercury's perihelion. For this purpose, we have introduced the notions of an extended and a geometrical perihelion. In addition, for each orbital period, a mean ellipse was fitted to the trajectory of Mercury. I have shown that the perihelion advance of Mercury deduced from the behavior of the Laplace-Runge-Lenz vector, as well as the extended and geometrical perihelion advance depend on the fitting time interval and for intervals of the order of 1000 years converge to a value of 532.1 arcseconds per century. The behavior of the perihelia, either extended or geometrical, is strongly impacted by the influence of Jupiter. The advance of the extended perihelion depends strongly on the time step used in the calculations, while the advance of the geometrical perihelion and that deduced by the rotation of the Laplace-Runge-Lenz vector depends only slightly on it.

Radio mini-haloes are a poorly-understood class of moderately-extended diffuse radio source that trace the presence of magnetic fields and relativistic electrons on scales of hundreds of kiloparsecs, predominantly in relaxed clusters. With relatively few confirmed detections to-date, many questions remain unanswered. This paper presents new radio observations of the galaxy cluster MS1455.0$+$2232 performed with MeerKAT at 1283 MHz and LOFAR at 145 MHz, the first results from a homogeneously selected mini-halo census. We find that the mini-halo in MS1455.0$+$2232 extends for around 590 kpc at 1283 MHz, significantly larger than previously believed, and has a flatter spectral index ($\alpha = -0.97 \pm 0.05$) than is typically expected. Our X-ray analysis clearly reveals a large-scale ($\sim254$ kpc) sloshing spiral in the intracluster medium. We perform a point-to-point analysis, finding a tight single correlation between radio and X-ray surface brightness, indicating a strong link between the thermal and non-thermal components of the intracluster medium. Conversely, in the spectral index/X-ray surface brightness plane, we find that regions inside and outside the sloshing spiral follow different correlations. We find compelling evidence for multiple sub-components in this mini-halo for the first time. While both the turbulent (re-)acceleration and hadronic scenarios are able to explain some observed properties of the mini-halo in MS1455.0$+$2232, neither scenario is able to account for all the evidence presented by our analysis.

The European Virtual Observatory (VO) initiative organises regular VO schools since 2008. The goals are twofold: i) to expose early-career European astronomers to the variety of currently available VO tools and services so that they can use them efficiently for their own research and; ii) to gather their feedback on the VO tools and services and the school itself. During the schools, VO experts guide participants on the usage of the tools through a series of predefined real science cases, an activity that took most of the allocated time. Participants also have the opportunity to develop their own science cases under the guidance of VO tutors. These schools have demonstrated to be very useful for students, since they declare to regularly use the VO tools in their research afterwards, and for us, since we have first hand information about the user needs. Here, we introduce our VO schools, the approach we follow, and present the training materials that we have developed along the years.

Numerous X-ray spectral models have been developed to model emission reprocessed by the torus of an active galactic nucleus (AGN), e.g., UXCLUMPY, CTORUS, and MYTORUS. They span a range of assumed torus geometries and morphologies-some posit smooth gas distributions, and others posit distributions of clouds. It is suspected that given the quality of currently available data, certain model parameters, such as coronal power law photon index and parameters determining the morphology of the AGN torus, may be poorly constrained due to model degeneracies. In this work, we test the reliability of these models in terms of recovery of parameters and the ability to discern the morphology of the torus using XMM-Newton and NuSTAR spectral data. We perform extensive simulations of X-ray spectra of Compton-thick AGNs under six X-ray spectral models of the torus. We use Bayesian methods to investigate degeneracy between model parameters, distinguish models and determine the dependence of the parameter constraints on the instruments used. For typical exposure times and fluxes for nearby Compton-thick AGN, we find that several parameters across the models used here cannot be well constrained, e.g., the distribution of clouds, the number of clouds in the radial direction, even when the applied model is correct. We also find that Bayesian evidence values can robustly distinguish between a correct and a wrong model only if there is sufficient energy coverage and only if the intrinsic flux of the object is above a particular value determined by the instrument combination and the model considered.

This paper presents the chemical abundance analysis of 217 stars in the metal-poor globular cluster NGC 6752, distributed from the turn-off to the lower red giant branch. Aluminium and Lithium abundances were derived through spectral synthesis applied to spectra collected with FLAMES, both in GIRAFFE and UVES modes. The work aims to gain insight into the nature of the polluter(s) responsible for the abundance variations and C-N, Na-O, Al-Mg anti-correlations associated with the multiple population phenomenon. We found a plateau at A(Li)=2.33$\pm$0.06 dex in unevolved stars, with the average Li content decreasing continuously down to $\sim$1.25 dex at the bottom of the red giant branch. As expected in the classic anti-correlation scenario, we found stars low in Al and high Li abundance and stars high in Al and low in Li. However, in addition, we also found evidence of Al-rich, second-generation stars with high Li content. This finding suggests the need for Li production, known to happen in intermediate-mass ($\sim$4-8M$_{sun}$) AGB stars through the Cameron-Fowler mechanism. It is worth noticing that the Li abundance observed in Al-rich stars never exceeds that in Al-poor stars.

Previous studies using Kepler data suggest that planets orbiting the same star tend to have similar sizes. However, due to the faintness of the stars, only a few of the planets were also detected with radial velocity follow-ups, and therefore the planetary masses were mostly unknown. It therefore yet to be determined whether planetary systems indeed behave as "peas in a pod". Follow-up programs of TESS targets significantly increased the number of confirmed planets with mass measurements, allowing for a more detailed statistical analysis of multi-planet systems. In this work we explore the similarity in radii, masses, densities, and period ratios of planets within planetary systems. We show that planets in the same system that are similar in radii could be rather different in mass and vice versa and that typically the planetary radii of a given planetary system are more similar than the masses. We also find that a transition in the "peas in the pod" pattern for planets more massive than ~100 Me and larger than ~10 Re. Planets below these limits are found to be significantly more uniform. We conclude that other quantities like the density may be crucial to fully understand the nature of planetary systems and that, due to the diversity of planets within a planetary system, increasing the number of detected systems is crucial for understanding the exoplanetary demographics.

We present a catalogue of 551,214 main-sequence stars in the local ($d < 2$ kpc) Galactic thick disk and halo, based on a search of stars with large proper motions ($\mu_{\mathrm{tot}} > 40.0$ mas/yr) in the Gaia Early Data Release 3. We derive photometric metallicity calibrated from the colour-luminosity-metallicity distribution of 20,047 stars with spectroscopic metallicities, collected from various spectroscopic surveys, including SDSS SEGUE/APOGEE, GALAH DR3, and LAMOST DR6. We combine these results to construct an empirical colour-magnitude-metallicity grid, which can be used to estimate photometric metallicities for low-mass metal-poor stars of K and M subtypes from their absolute $G$ magnitude and colour values. We find that low-mass, high-velocity stars in our catalogue share similar kinematics as reported in recent studies of more luminous Galactic halo stars. The pseudo-kinematic analysis of our sample recovers the main local halo structures, including the Gaia-Enceladus Stream and the Helmi stream; aside from these the local halo stars appear to show a remarkably smooth distribution in velocity space. Since the future Gaia data release will provide radial velocity measurements for only a small number of our sample, our catalogue provides targets of high interest for the future spectroscopic observation programs like SDSS-V, DESI MW survey, WEAVE, and/or 4MOST.

We introduce an ALMA band 3 spectroscopic survey, targeting the brightest submillimeter galaxies (SMGs) in the COSMOS field. Here we present the first results based on the 18 primary SMGs that have 870 $\mu$m flux densities of $S_{870}=12.4-19.3$ mJy and are drawn from a parent sample of 260 ALMA-detected SMGs from the AS2COSMOS survey. We detect emission lines in 17 and determine their redshifts to be in the range of $z=2-5$ with a median of $3.2\pm0.3$. We confirm that SMGs with brighter $S_{870}$ are located at higher redshifts. The data additionally cover five fainter companion SMGs, and we obtain line detection in one. Together with previous studies, our results indicate that for SMGs that satisfy our selection, their brightest companion SMGs are physically associated with their corresponding primary SMGs in $\ge40$% of the time, suggesting that mergers play a role in the triggering of star formation. By modeling the foreground gravitational fields, $<10$% of the primary SMGs can be strongly lensed with a magnification $\mu>2$. We determine that about 90% of the primary SMGs have lines that are better described by double Gaussian profiles, and the median separation of the two Gaussian peaks is 430$\pm$40 km s$^{-1}$. This allows estimates of an average baryon mass, which together with the line dispersion measurements puts our primary SMGs on the similar mass-$\sigma$ correlation found on local early-type galaxies. Finally, the number density of our $z>4$ primary SMGs is found to be $1\times10^6$ cMpc$^{-3}$, suggesting that they can be the progenitors of $z\sim3-4$ massive quiescent galaxies.

Silicon and Strontium are key elements to explore the nucleosynthesis and chemical evolution of the Galaxy by measurements of very metal-poor stars. There are, however, only a few useful spectral lines of these elements in the optical range that are measurable for such low-metallicity stars. Here we report on abundances of these two elements determined from near-infrared high-resolution spectra obtained with the Subaru Telescope Infrared Doppler instrument (IRD). Si abundances are determined for as many as 26 Si lines for six very and extremely metal-poor stars (-4.0<[Fe/H]<-1.5), which significantly improves the reliability of the abundance measurements. All six stars, including three carbon-enhanced objects, show over-abundances of Si ([Si/Fe]~+0.5). Two stars with [Fe/H]~-1.5 have relatively small over-abundances. The [Mg/Si] ratios agree with the solar value, except for one metal-poor star with carbon excess. Strontium abundances are determined from the triplet lines for four stars, including two for the first time. The consistency of the Sr abundances determined from near-infrared and optical spectra require further examination from additional observations.

Context. Solar Orbiter and PSP jointly observed the solar wind for the first time in June 2020, capturing data from very different solar wind streams, calm and Alfv\'enic wind as well as many dynamic structures. Aims. The aim here is to understand the origin and characteristics of the highly dynamic solar wind observed by the two probes, in particular in the vicinity of the heliospheric current sheet (HCS). Methods. We analyse the plasma data obtained by PSP and Solar Orbiter in situ during the month of June 2020. We use the Alfv\'en-wave turbulence MHD solar wind model WindPredict-AW, and perform two 3D simulations based on ADAPT solar magnetograms for this period. Results. We show that the dynamic regions measured by both spacecraft are pervaded with flux ropes close to the HCS. These flux ropes are also present in the simulations, forming at the tip of helmet streamers, i.e. at the base of the heliospheric current sheet. The formation mechanism involves a pressure driven instability followed by a fast tearing reconnection process, consistent with the picture of R\'eville et al. (2020a). We further characterize the 3D spatial structure of helmet streamer born flux ropes, which seems, in the simulations, to be related to the network of quasi-separatrices.

The properties of the interstellar medium (ISM) of the highest-redshift galaxies and quasars provide important indications of the complex interplay between the accretion of baryons onto galaxies, the physics that drives the build-up of stars out of this gas, the subsequent chemical evolution and feedback processes and the reionisation of the Universe. The Atacama Large Millimeter/submillimeter Array (ALMA) continues to play a pivotal role in the characterisation of the ISM of high-redshift galaxies. Observations of the dust continuum emission, atomic fine-structure and molecular lines arising from high-redshift galaxies are now carried out routinely, providing ever more constraints on the theoretical models of galaxy formation and evolution in the early Universe. The European Astronomical Society's EAS 2021 symposium dedicated to the exploration of the high-redshift Universe with ALMA provided a forum for the observational and theoretical high-redshift ALMA communities to exchange their views and recent results in this rapidly evolving field. The article summarises the exciting results that were presented at the meeting.

We study the capture of interstellar objects (ISOs) by a planet-star binary with mass ratio $q\ll1$, semi-major axis $a_p$, orbital speed $v_c$, and eccentricity $e_p$. Very close (slingshot) and wide encounters with the planet are amenable to analytic treatment, while numerically obtained capture cross-sections $\sigma$ closely follow the analytical results even in the intermediate regime. Wide interactions can only generate energy changes $\Delta E\lesssim qv_c^2$, when $\sigma\propto v_\infty^{-2} |\ln\Delta E|^{2/3}$ (with $v_\infty$ the ISO's incoming speed far away from the binary), which is slightly enhanced for $e_p>0$. Energy changes $\Delta E\gtrsim qv_c^2$, on the other hand, require close interactions when $\sigma\propto (v_\infty \Delta E)^{-2}$ hardly depending on $e_p$. Finally, at $\Delta E\gtrsim v_c^2$, the cross-section drops to zero, depending on the planet's radius $R_p$ through the Safronov number $\Theta=qa_p/R_p$. We also derive the cross-sections for collisions of ISOs with planets or moons.

This study builds on previous works, producing the most extensive prediction of the properties of such a hypothetical population to date, taking into account both Chandrasekhar and non-Chandrasekhar mass events. These results are then used to define criteria for membership of this population and characterise putative subpopulations This study contains 6x10^6 individual ejection trajectories out of the Galactic plane calculated with the stellar kinematics framework SHyRT, which are analysed with regard to their bulk observational properties. These are then put into context with the only previously identified population member US\,708 and applied to a number of other possible candidate objects. We find that two additional previously observed objects possess properties to warrant a designation as candidate objects. Characterisation of these object with respect to the predicted population finds all of them to be extreme in at least one astrometric observable. We find that current observations support a Galactic SN rate on the order of ~3x10^-7/yr to ~2x10^-6/yr, three orders of magnitude below the inferred Galactic SN Ia rate and two orders of magnitude below the formation rate of predicted He-donor progenitors. The number of currently observed population members suggests that the He-donor scenariois not a dominant contributor to the number of observed SNe Ia. However, we find that, even at the low event rate suggested, the majority of possibly detectable population members is still undetected. The extreme nature of current population members suggests that a still larger number of objects has simply evaded detection up to this point, hinting at a higher contribution than is currently supported by observation. - abridged -

Stellar orbits at the Galactic Center provide a very clean probe of the gravitational potential of the supermassive black hole. They can be studied with unique precision, beyond the confusion limit of a single telescope, with the near-infrared interferometer GRAVITY. Imaging is essential to search the field for faint, unknown stars on short orbits which potentially could constrain the black hole spin. Furthermore, it provides the starting point for astrometric fitting to derive highly accurate stellar positions. Here, we present $\mathrm{G^R}$, a new imaging tool specifically designed for Galactic Center observations with GRAVITY. The algorithm is based on a Bayesian interpretation of the imaging problem, formulated in the framework of information field theory and building upon existing works in radio-interferometric imaging. Its application to GRAVITY observations from 2021 yields the deepest images to date of the Galactic Center on scales of a few milliarcseconds. The images reveal the complicated source structure within the central $100\,\mathrm{mas}$ around Sgr A*, where we detected the stars S29 and S55 and confirm S62 on its trajectory, slowly approaching Sgr A*. Furthermore, we were able to detect S38, S42, S60, and S63 in a series of exposures for which we offset the fiber from Sgr A*. We provide an update on the orbits of all aforementioned stars. In addition to these known sources, the images also reveal a faint star moving to the west at a high angular velocity. We cannot find any coincidence with any known source and, thus, we refer to the new star as S300. From the flux ratio with S29, we estimate its K-band magnitude as $m_\mathrm{K}\left(\mathrm{S300}\right)\simeq 19.0 - 19.3$. Images obtained with CLEAN confirm the detection.

The stars orbiting the compact radio source Sgr A* in the Galactic Centre are precision probes of the gravitational field around the closest massive black hole. In addition to adaptive optics assisted astrometry (with NACO / VLT) and spectroscopy (with SINFONI / VLT, NIRC2 / Keck and GNIRS / Gemini) over three decades, since 2016/2017 we have obtained 30-100 mu-as astrometry with the four-telescope interferometric beam combiner GRAVITY / VLTI reaching a sensitivity of mK = 20 when combining data from one night. We present the simultaneous detection of several stars within the diffraction limit of a single telescope, illustrating the power of interferometry. The new data for the stars S2, S29, S38 and S55 yield significant accelerations between March and July 2021, as these stars pass the pericenters of their orbits between 2018 and 2023. This allows for a high-precision determination of the gravitational potential around Sgr A*. Our data are in excellent agreement with general relativity orbits around a single central point mass, M = 4.30 x 10^6 M_sun with a precision of about +-0.25%. We improve the significance of our detection of the Schwarzschild precession in the S2 orbit to 7 sigma. Assuming plausible density profiles, an extended mass component inside S2's apocentre (= 0.23" or 2.4 x 10^4 R_S) must be 3000 M_sun (1 sigma), or 0.1% of M. Adding the enclosed mass determinations from 13 stars orbiting Sgr A* at larger radii, the innermost radius at which the excess mass beyond Sgr A* tentatively is seen is r = 2.5" >= 10x the apocentre of S2. This is in full harmony with the stellar mass distribution (including stellar-mass black holes) obtained from the spatially resolved luminosity function.

The radioactively powered transient following a binary neutron star (BNS) merger, known as a kilonova (KN), is expected to enter the steady-state nebular phase a few days after merger. Steady-state holds until thermal reprocessing time-scales become long, at which point the temperature and ionisation states need to be evolved time-dependently. We study the onset and significance of time-dependent effects using the non-local thermodynamic equilibrium (NLTE) spectral synthesis code SUMO. We employ a simple single-zone model with an elemental composition of Te, Ce, Pt and Th, scaled to their respective solar abundances. The atomic data are generated using the Flexible Atomic Code (FAC), and consist of energy levels and radiative transitions, including highly forbidden lines. We explore the KN evolution from 5 to 100 days after merger, varying ejecta mass and velocity. We also consider variations in the degree of electron magnetic field trapping, as well as radioactive power generation for alpha and beta decay (but omitting fission products). We find that the transition time, and magnitude of steady-state deviations are highly sensitive to these parameters. For typical KN ejecta, the deviations are minor within the time-frame studied. However, low density ejecta with low energy deposition show significant differences from $\sim 10$ days. Important deviation of the ionisation structure solution impacts the temperature by altering the overall line cooling. Adiabatic cooling becomes important at $t \geq 60$ days which, in addition to the temperature and ionisation effects, lead to the bolometric light curve deviating from the instantaneous radioactive power deposited.

Capture of interstellar objects (ISOs) into the Solar system is dominated by ISOs with asymptotic incoming speeds $v_\infty<4\,$km\,s$^{-1}$. The capture rate is proportional to the ISO phase-space density in the Solar vicinity and does not vary along the Sun's Galactic orbit, i.e.\ is not enhanced during a passage through a cloud of ISOs (in contrast to previous suggestions). Most bound orbits crossing those of Jupiter and Saturn are fully mixed with unbound phase space, implying that they hold the same ISO phase-space density. Assuming an interstellar number density $n_{iso}\sim0.1\,$au$^{-3}$, we estimate that in 1000 years the planets capture $\sim2$ ISOs (while $\sim17$ fall into the Sun), resulting in a population of $\sim8$ captured ISOs within 5\,au of the Sun at any time, less than the number of visiting ISOs passing through the same volume on hyperbolic orbits. In terms of phase-space volume, capture onto and ejection from the Solar system are equal, such that on average ISOs will not remain captive at $a\lesssim2000\,$au for extensive periods.

Proper photometric data are challenging to obtain in the K2 mission of the Kepler space telescope due to strong systematics caused by the two-wheel-mode operation. It is especially true for variable stars wherein physical phenomena occur on timescales similar to the instrumental signals. We originally developed a method with the aim to extend the photometric aperture to be able to compensate the motion of the telescope which we named Extended Aperture Photometry (EAP). Here we present the outline of the automatized version of the EAP method, an open-source pipeline called autoEAP. We compare the light curve products to other photometric solutions for examples chosen from high-amplitude variable stars. Besides the photometry, we developed a new detrending method, which is based on phase dispersion minimization and is able to eliminate long-term instrumental signals for periodic variable stars.

In this article we perform a simulation-based study of the shape of the electromagnetic wavefront of air showers induced by cosmic rays. We show that for showers with zenith angles larger than 60{\deg}, a point-source like description of the wave-front is sufficient and that the reconstructed position of this point-source is a potential proxy for the determination of the nature of the cosmic rays initiating the showers, with performances similar to those obtained using the depth of shower maximum.

The number of stars with a measured value of [Fe/H] is considerably increasing thanks to spectroscopic surveys. However different methodologies, inputs and assumptions used in spectral analyses lead to different precisions in [Fe/H] and possibly to systematic differences that need to be evaluated. It is essential to understand the characteristics of each survey to fully exploit their potential, in particular if the surveys are combined. The purpose of this study is to compare [Fe/H] determinations from the largest spectroscopic surveys (APOGEE, GALAH, the Gaia ESO survey, RAVE, LAMOST, SEGUE ) to other catalogues taken as reference. Offsets and dispersions of the residuals are examined as well as their trends with other parameters. We use reference samples providing independent determinations of [Fe/H] which are compared to those from the surveys for common stars. The distribution of the residuals is assessed through simple statistics that measures the offset between two catalogues and the dispersion representative of the precision of both catalogues. When relevant, linear fits are performed. A large sample of FGK-type stars with [Fe/H] based on high-resolution, high signal to noise spectroscopy was built from the PASTEL catalogue to provide a reference sample. We also use FGK members in open and globular clusters to assess the internal consistency of [Fe/H] of each survey. The agreement of median [Fe/H] values for clusters observed by different surveys is discussed. All the surveys overestimate the low metallicities, and some of them also underestimate the high metallicities. They perform well in the most populated intermediate metallicity range, whatever the resolution. In most cases the typical precision that we deduce from the comparisons is in good agreement with the uncertainties quoted in the catalogues. Some exceptions to this general behaviour are discussed.

We investigate the primordial power spectra for general kinetic inflation models that support a period of kinetic dominance in the case of curved universes. We present derivations of the Mukhanov-Sasaki equations with a non-standard scalar kinetic Lagrangian which manifests itself through the inflationary sound speed $c_s^2$. We extend the analytical approximations exploited in Contaldi et al [1] and Thavanesan et al [2] to general kinetic Lagrangians and show the effect of k-inflation on the primordial power spectra for models with curvature. In particular, the interplay between sound speed and curvature results in a natural low wavenumber cutoff for the power spectra in the case of closed universes. Using the analytical approximation, we further show that a change in the inflationary sound speed between different epochs of inflation results in non-decaying oscillations in the resultant power spectra for the comoving curvature perturbation.

We present a modelization of the non-thermal pressure, $P_{NT}$, and we apply it to the X-ray (and Sunayev-Zel'dovich) derived radial profiles of the X-COP galaxy clusters. We relate the amount of non-thermal pressure support to the hydrostatic bias, $b$, and speculate on how we can interpret this $P_{NT}$ in terms of the expected levels of turbulent velocity and magnetic fields. Current upper limits on the turbulent velocity in the intracluster plasma are used to build a distribution $\mathcal{N}(<b) - b$, from which we infer that 50 per cent of local galaxy clusters should have $b < 0.2$ ($b<0.33$ in 80 per cent of the population). The measured bias in the X-COP sample that includes relaxed massive nearby systems is 0.03 in 50% of the objects and 0.17 in 80% of them. All these values are below the amount of bias required to reconcile the observed cluster number count in the cosmological framework set from Planck.

The Extreme Universe Space Observatory Supper Pressure Balloon 2 (EUSO-SPB2) is under development, and will prototype instrumentation for future satellite-based missions, including the Probe of Extreme Multi-Messenger Astrophysics (POEMMA). EUSO-SPB2 will consist of two telescopes. The first is a Cherenkov telescope (CT) being developed to identify and estimate the background sources for future below-the-limb very high energy (E>10 PeV) astrophysical neutrino observations, as well as above-the-limb cosmic ray induced signals (E>1 PeV). The second is a fluorescence telescope (FT) being developed for detection of Ultra High Energy Cosmic Rays (UHECRs). In preparation for the expected launch in 2023, extensive simulations tuned by preliminary laboratory measurements have been preformed to understand the FT capabilities. The energy threshold has been estimated at $10^{18.2}$ eV, and results in a maximum detection rate at $10^{18.6}$ eV when taking into account the shape of the UHECR spectrum. In addition, onboard software has been developed based on the simulations as well as experience with previous EUSO missions. This includes a level 1 trigger to be run on the computationally limited flight hardware, as well as a deep learning based prioritization algorithm in order to accommodate the balloon's telemetry budget. These techniques could also be used later for future, space-based missions.

A major source of uncertainty in the age determination of old ($\sim10$ Gyr) integrated stellar populations is the presence of hot horizontal branch (HB) stars. Here, we describe a simple approach to tackle this problem, and show the performance of this technique that simultaneously models the age, abundances and HB properties of integrated stellar populations. For this we compare the results found during the fits of the integrated spectra of a sample of stellar population benchmarks, against the values obtained from the analysis of their resolved CMDs. We find that the ages derived from our spectral fits for most (26/32) of our targets are within 0.1 dex to their CMDs values. Similarly, for the majority of the targets in our sample we are able to recover successfully the flux contribution from hot HB stars (within $\sim0.15 {\rm ~dex}$ for 18/24 targets) and their mean temperature (14/24 targets within $\sim30 \%$). Finally, we present a diagnostic that can be used to detect spurious solutions in age, that will help identify the few cases when this method fails. These results open a new window for the detailed study of globular clusters beyond the Local Group.

The objects that emit extragalatic fast radio bursts (FRBs) remain unidentified. Studies of the host galaxies and environments of accurately localised ($\lesssim1$ arcsec) FRBs promise to deliver critical insights into the nature of their progenitors. Here we demonstrate the effects of observational selection biases on analyses of the distributions of FRB host-galaxy properties (including star-formation rate, SFR, and stellar mass, $M_{*}$), and on the distributions of FRB offsets from the centres of their hosts. We consider the effects of "radio selection", wherein FRBs with larger dispersion measures and scattering timescales are less likely to be detected, and the effects of "optical selection", wherein FRBs with fainter host galaxies are more likely to have unidentified or mis-identified hosts. We develop a plausible, illustrative model for these effects in observations of FRBs and their host galaxies by combining the output catalogues of a semi-analytic galaxy formation model with a recently developed algorithm to associate FRBs with host galaxies (PATH). We find that optical selection biases are most important for the host-galaxy $M_{*}$ and SFR distributions, and that radio selection biases are most important for the distribution of FRB projected physical offsets. For our fiducial simulation of FRBs at $z<0.5$, the selection biases cause the median host-galaxy SFR to be increased by $\sim0.3$ dex, and the median $M_{*}$ by $\sim0.5$ dex. The median projected physical offset is increased by $\sim2$ kpc ($\sim0.25$ dex). These effects are sufficiently large so as to merit careful consideration in studies of localised FRBs, and our simulations provide a guide towards their mitigation.

In this paper we first review the results obtained by the INTEGRAL mission in the domain of Gamma-Ray Bursts (GRBs), thanks to the INTEGRAL Burst Alert System, which is able to deliver near real-time alerts for GRBs detected within the IBIS field of view. More than 120 GRBs have been detected to date and we summarize their properties here. In the second part of this review we focus on the polarimetric results obtained by IBIS and SPI on GRBs and Galactic compact objects.

In modified gravity models that allow for additional non-compact spacetime dimensions, energy from gravitational waves can leak into these extra spacetime dimensions, leading to a reduction in the amplitude of the observed gravitational waves, and thus a source of potential systematics in the inferred luminosity distances to gravitational wave sources. Since binary black hole (BBH) mergers are standard sirens, we use the pair-instability supernova (PISNe) mass gap and its predicted features to determine a mass scale and thus be able to break the mass-redshift degeneracy. We simultaneously fit for the BBH population and the extra spacetime dimensions parameters from gravitational leakage models using BBH observations from the recently released GWTC-3 catalog. We set constraints on the number of spacetime dimensions and find that $D= 3.95^{+0.09}_{-0.07}$ at $68\%$ C.L. for models that are independent of a screening scale, finding that the GWTC-3 constraint is as competitive as that set from GW170817 and its electromagnetic counterpart. For models where gravity leaks below a certain screening scale $R_c$, we find $D=4.23^{+1.50}_{-0.57}$ and $\log_{10} R_c/{\rm Mpc}= 4.14^{+0.55}_{-0.86}$ with a transition steepness $\log_{10} n = 0.86^{+0.73}_{-0.84}$ for the leakage, which for the first are constrained jointly with the BBH population at cosmological distances. These constraints are consistent with General Relativity (GR) where gravitational waves propagate in $D=3+1$ spacetime dimensions. Using the BBH population to probe modifications to standard cosmological models provides an independent test of GR that does not rely on any electromagnetic information but purely on gravitational wave observations.

The joint detection of gravitational waves and the gamma-ray counterpart of a binary neutron star merger event, GW170817, unambiguously validates the connection between short gamma-ray bursts and compact binary object (CBO) mergers. We focus on a special scenario where short gamma-ray bursts produced by CBO mergers are embedded in disks of active galactic nuclei (AGN), and we investigate the $\gamma$-ray emission produced in the internal dissipation region via synchrotron, synchrotron self-Compton and external inverse-Compton (EIC) processes. In this scenario, isotropic thermal photons from the AGN disks contribute to the EIC component. We show that a low-density cavity can be formed in the migration traps, leading to the embedded mergers producing successful GRB jets. We find that the EIC component would dominate the GeV emission for typical CBO mergers with an isotropic-equivalent luminosity of $L_{j,\rm iso}=10^{48.5}~\rm erg~s^{-1}$ which are located close to the central supermassive black hole. Considering a long-lasting jet of duration $T_{\rm dur}\sim10^2-10^3$ s, we find that the future CTA will be able to detect its $25-100$ GeV emission out to a redshift $z=1.0$. In the optimistic case, it is possible to detect the on-axis extended emission simultaneously with GWs within one decade using MAGIC, H.E.S.S., VERITAS, CTA, and LHAASO-WCDA. Early diagnosis of prompt emissions with $Fermi$-GBM and HAWC can provide valuable directional information for the follow-up observations.

We explore the significance of noise from thermal Sunyaev-Zel{'}dovich (tSZ) signals for cluster detection using cosmic microwave background (CMB) surveys. The noise arises both from neighboring objects and also from haloes below the detection limit. A wide range of surveys are considered: SPT-SZ, SPTpol, and SPT-3G from the South Pole Telescope; SO-Baseline and SO-Goal configurations for Simons Observatory; CMB-S4's wide area (S4-Wide) and deep (S4-Ultra deep) surveys; and the futuristic CMB-HD experiment. We find that the noise from tSZ signals has a significant impact on CMB-HD and to some extent on S4-Ultra deep. For other experiments, the effect is negligible as the noise in the tSZ map is dominated by residual foregrounds or experimental noise. In the limit when the noise from tSZ signals is important, we find that removing the detected clusters and rerunning the cluster finder allows us to find a new set of less massive and distant clusters. Since the detected clusters are the dominant source of the tSZ power, removing them reduces the power at $\ell = 3000$ by: $\times5$ for CMB-HD; $\times 3.1$ of S4-Ultra deep; $\times2.4$ for S4-Wide and SPT-3G; $\times1.5$ for SO-Goal and SPTpol; $\times1.35$ for SO-Baseline; and $\times1.08$ for SPT-SZ. We forecast the expected number of clusters and also derive parameter constraints by combining cluster counts with primary CMB and tSZ power spectra finding that the future surveys can reduce the error on the dark energy equation of state parameter to sub-percent levels and can also enable $\ge3\sigma$ detection of the sum of neutrino masses. The simulation products and results can be downloaded from https://github.com/sriniraghunathan/tSZ_cluster_forecasts .

Gravitational acceleration fields can be deduced from the collisionless Boltzmann equation, once the distribution function is known. This can be constructed via the method of normalizing flows from datasets of the positions and velocities of stars. Here, we consider application of this technique to the solar neighbourhood. We construct mock data from a linear superposition of multiple `quasi-isothermal' distribution functions, representing stellar populations in the equilibrium Milky Way disc. We show that given a mock dataset comprising a million stars within 1 kpc of the Sun, the underlying acceleration field can be measured with excellent, sub-percent level accuracy, even in the face of realistic errors and missing line-of-sight velocities. The effects of disequilibrium can lead to bias in the inferred acceleration field. This can be diagnosed by the presence of a phase space spiral, which can be extracted simply and cleanly from the learned distribution function. We carry out a comparison with two other popular methods of finding the local acceleration field (Jeans analysis and 1D distribution function fitting). We show our method most accurately measures accelerations from a given mock dataset, particularly in the presence of disequilibria.

Upcoming deep imaging surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time will be confronted with challenges that come with increased depth. One of the leading systematic errors in deep surveys is the blending of objects due to higher surface density in the more crowded images; a considerable fraction of the galaxies which we hope to use for cosmology analyses will overlap each other on the observed sky. In order to investigate these challenges, we emulate blending in a mock catalogue consisting of galaxies at a depth equivalent to 1.3 years of the full 10-year Rubin Observatory that includes effects due to weak lensing, ground-based seeing, and the uncertainties due to extraction of catalogues from imaging data. The emulated catalogue indicates that approximately 12% of the observed galaxies are "unrecognized" blends that contain two or more objects but are detected as one. Using the positions and shears of half a billion distant galaxies, we compute shear-shear correlation functions after selecting tomographic samples in terms of both spectroscopic and photometric redshift bins. We examine the sensitivity of the cosmological parameter estimation to unrecognized blending employing both jackknife and analytical Gaussian covariance estimators. A $\sim0.02$ decrease in the derived structure growth parameter $S_8 = \sigma_8 (\Omega_{\rm m}/0.3)^{0.5}$ is seen due to unrecognized blending in both tomographies with a slight additional bias for the photo-$z$-based tomography. This bias is about 2$\sigma$ statistical error in measuring $S_8$.

Mapping exoplanets across phases and during secondary eclipse is a powerful technique for characterizing Hot Jupiters in emission. Since these planets are expected to rotate about axes normal to their orbital planes, with rotation periods synchronized with their orbital periods, mapping provides a direct correspondence between orbital phase and planetary longitude. For planets with fewer constraints on their spin properties, the relationship between the shape of the eclipse light curve and the visible portion of the surface is more complex. We develop a framework to understand the information content of planets where the rotation rate and/or axis orientation are not well constrained, by constructing a basis of emission time series ("light curves") that are orthogonal in integrated flux across secondary eclipse. The most orthogonal bases in eclipse consist of periodic functions - akin to sinusoids - and at slow enough rotation rates these follow a monotonic series in frequency. If only data during eclipse are considered, we show that very similar light curves at a given frequency can be generated by maps of similar complexity at a wide range of spin axis orientations. Constraining spin axis orientations from eclipse data alone may therefore depend on strong prior knowledge of plausible emission map structures and/or rotation rates. When the planetary rotation period becomes shorter than the total eclipse duration, there will be an additional ambiguity between the complexity of the emission map and observable variations due to rotation, as both can in principle generate similar signals. By modeling example eclipse observations of the Warm Jupiter HAT-P-18 b, we demonstrate that the available signal-to-noise for $\sim 10$ orbits is just sufficient to derive map structure beyond the eclipse depth.

The growing interest in the interactions between dark matter particles and electrons has received a further boost by the observation of an excess in electron recoil events in the XENON1T experiment. Of particular interest are dark matter models in which the scattering process is inelastic, such that the ground state can upscatter into an excited state. The subsequent exothermic downscattering of such excited states on electrons can lead to observable signals in direct detection experiments and gives a good fit to the XENON1T excess. In this work, we study terrestrial upscattering, i.e. inelastic scattering of dark matter particles on nuclei in the Earth, as a plausible origin of such excited states. Using both analytical and Monte Carlo methods, we obtain detailed predictions of their density and velocity distribution. These results enable us to explore the time dependence of the flux of excited states resulting from the rotation of the Earth. For the case of XENON1T, we find the resulting daily modulation of the electron recoil signal to be at the level of 10% with a strong dependence on the dark matter mass.

We show that macroscopic dark matter (DM) impacts on the degenerate helium cores of red-giant branch (RGB) stars can ignite helium fusion via DM-baryon elastic scattering. The onset of helium burning leads to a characteristic drop in luminosity and rise in temperature that marks the transition to a horizontal branch star. We show that such impacts can alter the RGB luminosity function of globular clusters (GCs), focusing in particular on the GC M15. Using models of M15 stars constructed with the stellar simulation code MESA, we compute the expected DM-ignition event rates and the theoretical RGB luminosity functions under the null and signal hypotheses. We constrain DM with masses $10^{17}\ {\rm g} \lesssim m_{\chi} \lesssim 10^{20}\ \rm{g}$ and geometric cross sections $10^2\ {\rm cm}^2 \lesssim \sigma_{\chi n} \lesssim 10^{7}\ \rm{ cm}^2 $ assuming that the DM in M15 is sourced by the background Milky Way halo. We also place more stringent constraints assuming that M15 formed in a DM subhalo that survives today.

We construct an inspiral-merger-ringdown eccentric gravitational-wave (GW) model for binary black holes with non-precessing spins within the effective-one-body formalism. This waveform model, SEOBNRv4EHM, extends the accurate quasi-circular SEOBNRv4HM model to eccentric binaries by including recently computed eccentric corrections up to 2PN order in the gravitational waveform modes, notably the $(l,|m|)=(2,2),(2,1),(3,3),(4,4),(5,5)$ multipoles. The waveform model reproduces the zero eccentricity limit with an accuracy comparable to the underlying quasi-circular model, with the unfaithfulness of $\lesssim1\%$ against quasi-circular numerical-relativity (NR) simulations. When compared against 28 public eccentric NR simulations from the Simulating eXtreme Spacetimes catalog with initial orbital eccentricities up to $e\simeq0.3$ and dimensionless spin magnitudes up to $+0.7$, the model provides unfaithfulness $<1\%$, showing that both the $(2,|2|)$-modes and the higher-order modes are reliably described without calibration to NR datasets in the eccentric sector. The waveform model SEOBNRv4EHM is able to qualitatively reproduce the phenomenology of dynamical captures, and can be extended to include spin-precession effects. It can be employed for upcoming observing runs with the LIGO-Virgo-KAGRA detectors and used to re-analyze existing GW catalogs to infer the eccentricity parameters for binaries with $e\lesssim0.3$ (at 20 Hz or lower) and spins up to $\lesssim 0.9-0.95$. The latter is a promising region of the parameter space where some astrophysical formation scenarios of binaries predict mild eccentricity in the ground-based detectors' bandwidth. Assessing the accuracy and robustness of the eccentric waveform model SEOBNRv4EHM for larger eccentricities and spins will require comparisons with, and, likely, calibration to eccentric NR waveforms in a larger region of the parameter space.

We have attempted to visualize transport coefficients like shear viscosity and electrical conductivity with respect to density of neutron star environment, whose core may expect a hadron-quark phase transition due to acquiring much higher density than nuclear saturation density. By making sandwich between MIT bag model for quark phase and two different effective hadronic models for hadronic phase, we have estimated the transport coefficients of two phases. During sketching of the transport coefficients, we have discussed their detailed density profile of phase-space part and relaxation time part. By calculating shear viscosity to density ratio, we have also explored the nearly perfect fluid domain along the density axis of hadron-quark phase diagram.

Liquid argon (LAr) is a common choice as detection medium in particle physics and rare-event searches. Challenges of LAr scintillation light detection include its short emission wavelength, long scintillation time and short attenuation length. The addition of small amounts of xenon to LAr is known to improve the scintillation and optical properties. We present a characterization campaign on xenon-doped liquid argon (XeDLAr) with target xenon concentrations ranging from 0 to 300 ppm by mass encompassing the measurement of the photoelectron yield $Y$ , effective triplet lifetime $\tau_3$ and effective attenuation length $\lambda_\mathrm{att}$. The measurements were conducted in the Subterranean Cryogenic ARgon Facility, SCARF, a 1 t (XeD)LAr test stand in the shallow underground laboratory (UGL) of TU-Munich. These three scintillation and optical parameters were observed simultaneously with a single setup, the Legend Liquid Argon Monitoring Apparatus, LLAMA. The actual xenon concentrations in the liquid and gaseous phases were determined with the Impurity DEtector For Investigation of Xenon, IDEFIX, a mass spectrometer setup, and successful doping was confirmed. At the highest dopant concentration we find a doubling of $Y$ , a tenfold reduction of $\tau_3$ to $\sim$ 90 ns and a tenfold increase of $\lambda_{att}$ to over 6 m.

The LISA Pathfinder (LPF) mission succeeded outstandingly in demonstrating key technological aspects of future space-borne gravitational-wave detectors, such as the Laser Interferometer Space Antenna (LISA). Specifically, LPF demonstrated with unprecedented sensitivity the measurement of the relative acceleration of two free-falling cubic test masses. Although most disruptive non-gravitational forces have been identified and their effects mitigated through a series of calibration processes, some faint transient signals of yet unexplained origin remain in the measurements. If they appear in the LISA data, these perturbations (also called glitches) could skew the characterization of gravitational-wave sources or even be confused with gravitational-wave bursts. For the first time, we provide a comprehensive census of LPF transient events. Our analysis is based on a phenomenological shapelet model allowing us to derive simple statistics about the physical features of the glitch population. We then implement a generator of synthetic glitches designed to be used for subsequent LISA studies, and perform a preliminary evaluation of the effect of the glitches on future LISA data analyses.

This paper presents a new method to process the data recorded with muon telescopes. We have developed this processing method for the plastic scintillator-based hodoscopes located around the volcano La Soufri\`ere de Guadeloupe, in the French Lesser Antilles, in order to perform muon radiographies of the lava dome region, strongly impacted by the volcanic hydrothermal activity. Our method relies on particle trajectory reconstruction, performing a fit of the recorded hits in the impacted scintillator bars using a Random Sample Consensus algorithm. This algorithm is specifically built to discriminate outlier points, usually due to noise hits, in the data. Thus, it is expected to significantly improve the signal/noise separation in muon track hits and to obtain higher quality estimates of the particles' incident trajectories in our detectors. The first analysis of the RANSAC-reconstructed events offers promising results in terms of average density maps. To illustrate the performances of this algorithm, we provide angular resolution and reconstruction efficiency estimates using a GEANT4 simulation of a telescope equipped with four detection matrices. In addition, we also show preliminary results from open-sky data recorded with such telescope at La Soufri\`ere de Guadeloupe volcano.

We summarize our work on the generation of primordial black holes in a type of two-field inflationary models. The key ingredient is a sharp turn of the background trajectory in field space. We show that certain classes of solutions to the equations of motion exhibit precisely this kind of behavior. Among them we find solutions, which describe a transition between an ultra-slow roll and a slow roll phases of inflation.

Gravitational lensing describes the bending of the trajectories of light and gravitational waves due to the gravitational potential of a massive object. Strong lensing by galaxies can create multiple images with different overall amplifications, arrival times, and image types. If, furthermore, the gravitational wave encounters a star along its trajectory, microlensing will take place. Our previous research studied the effects of microlenses on strongly-lensed type I images. We extend our research to type II strongly-lensed images to complete the story. Our results are broadly consistent with prior work by other groups. As opposed to being magnified, the type II images are typically demagnified. Moreover, the type II images produce larger mismatches than type I images. Similar to our prior work, we find that the wave optics effects significantly suppress microlensing in the stellar-mass limit. We also discuss the implication of our results for a particularly promising method to use strong lensing to detect microlensing. In the future, it will be crucial to incorporate these microlensed waveforms in gravitational-wave lensing searches.

We study fermions derivatively coupled to axion-like or pseudoscalar fields, and show that the axial vector current of the fermions is not conserved in the limit where the fermion is massless. This violation of the classical chiral symmetry is due to the background axion field. We compute the contributions to this anomalous Ward identity due to the pseudoscalar field alone, which arise in Minkowski space, as well as the effects due to interaction with an external gravitational field. In all cases, these interactions induce terms in the axion effective action that can be removed by the addition of local counterterms. We perform our computations both perturbatively using Feynman graphs, as well as by studying the transformation properties of the path integral measure. Using the heat kernel method, we include the effects of gravity as well as gauge fields, and compute the anomaly. Finally, we verify our relation by considering derivatively coupled fermions during pseudoscalar-driven inflation and computing the divergence of the axial current in de Sitter spacetime.