Papers for Wednesday, Nov 03 2021

With a single training image and using wavelet phase harmonic augmentation, we present polarized Cosmic Microwave Background (CMB) foreground marginalization in a high-dimensional likelihood-free (Bayesian) framework. We demonstrate robust foreground removal using only a single frequency of simulated data for a BICEP-like sky patch. Using Moment Networks we estimate the pixel-level posterior probability for the underlying {E,B} signal and validate the statistical model with a quantile-type test using the estimated marginal posterior moments. The Moment Networks use a hierarchy of U-Net convolutional neural networks. This work validates such an approach in the most difficult limiting case: pixel-level, noise-free, highly non-Gaussian dust foregrounds with a single training image at a single frequency. For a real CMB experiment, a small number of representative sky patches would provide the training data required for full cosmological inference. These results enable robust likelihood-free, simulation-based parameter and model inference for primordial B-mode detection using observed CMB polarization data.

Faint extended elliptically-shaped ultra-diffuse galaxies and slightly brighter and more compact dwarf elliptical and lenticular stellar systems are common in galaxy clusters. Their poorly constrained evolutionary paths can be studied by identifying young ultra-diffuse galaxy and dwarf elliptical analogs populated with bright, massive stars. Using data mining we identified 11 such low-mass ($2\times10^8<M_*<2\times10^9 M_{\odot}$) galaxies with large half-light radii ($2.0<r_e<5$~kpc) and recently quenched star formation in the Coma and Abell~2147 galaxy clusters. All galaxies happen to have ram-pressure-stripped tails with signs of current or recent star formation. Deep spectroscopic observations revealed rotating stellar discs containing 70--95\% dark matter by mass. A large fraction of the disc stars (10--60\%) formed in intense star bursts 180--970~Myr ago, probably triggered by ram pressure. Observed global gradients of stellar age corroborate this scenario. Passive evolution in the next 10~Gyr will transform 9 of the 11 galaxies into ultra-diffuse galaxies. If we assume a constant rate of galaxy infall, 44$\pm$16~\%\ of the most luminous present-day ultra-diffuse galaxies in Coma must have formed via ram pressure stripping of disky progenitors.

Users of the JKTEBOP code to solve the light curves of eclipsing binaries often confuse the gravity darkening coefficients, $y(\lambda)$, with the bolometric gravity darkening exponents, $\beta$. JKTEBOP requires the wavelength-dependent coefficients. I show that the numerical values of $y(\lambda)$ and $\beta$ can be rather different, leading to potential biases in the solution if the wrong quantities are used.

Observations have shown that spatially extended "TeV halos" are a common (and potentially generic) feature surrounding young and middle-aged pulsars. However, their morphology is not understood. They are larger than the "compact" region where the stellar remnant dominates the properties of the interstellar medium, but smaller than expected in models of cosmic-ray diffusion through the standard interstellar medium. Several explanations have been proposed, but all have shortcomings. Here, we revisit a class of models where the cosmic-ray gradient produced by the central source induces a streaming stability that "self-confines" the cosmic-ray population. We find that previous studies significantly underpredicted the degree of cosmic-ray confinement and show that corrected models can significantly inhibit cosmic-ray diffusion throughout the TeV halo, especially when similar contributions from the coincident supernova remnant are included.

We analyze new multi-filter Hubble Space Telescope (HST) photometry of the normal Type Ia supernova (SN Ia) 2011fe out to $\approx 2400~$days after maximum light, the latest observations to-date of a SN Ia. We model the pseudo-bolometric light curve with a simple radioactive decay model and find energy input from both $^{57}$Co and $^{55}$Fe are needed to power the late-time luminosity. This is the first detection of $^{55}$Fe in a SN Ia. We consider potential sources of contamination such as a surviving companion star or delaying the deposition timescale for $^{56}$Co positrons but these scenarios are ultimately disfavored. The relative isotopic abundances place direct constraints on the burning conditions experienced by the white dwarf and require a central density of $\rho_c \sim 3\times10^8~\rm{g}~\rm{cm}^{-3}$. Only 2 classes of explosion models are currently consistent with all observations of SN2011fe: 1) the delayed detonation of a low-$\rho_c$, near-$\rm{M}_{\rm{Ch}}$ ($1.2-1.3~M_\odot$) WD, or 2) a sub-$\rm{M}_{\rm{Ch}}$ ($1.0-1.1~M_\odot$) WD experiencing a thin-shell double detonation.

We present the first detailed 3D kinematic analysis of a sample of 3133 white dwarfs that use Gaia astrometry plus radial velocities, which were measured either by Gaia or by ground based spectroscopic observations. The studied sample includes either isolated white dwarfs, which have direct radial velocity measurements, or white dwarfs that belong to common proper motion pairs that contain non-degenerate companions with available radial velocities. A subset of common proper motion pairs also have metal abundances that have been measured by large scale spectroscopic surveys or by our own follow-up observations. We used the white dwarfs as astrophysical clocks, by determining their masses and total ages via interpolation with state-of-the-art evolutionary models, and we used the non-degenerate companions in common proper motions for chemical tagging of the population. Combining accurate radial velocities with Gaia astrometry and proper motions, we derived the velocity components of our sample in the Galactic rest frame and their Galactic orbital parameters. The studied sample is mostly located within ~300 pc from the Sun. It contains predominantly (90-95 %) thin disk stars with close-to-circular Galactic orbits, while the remaining 5-10 % of stars have more eccentric trajectories and belong to the thick disk. We identified seven isolated white dwarfs and two common proper motion pairs as halo members. We determined the age velocity-dispersion relation for the thin disk members, which is in agreement with previous results that were achieved from different white dwarf samples without published radial velocities and shows signatures of dynamical heating and saturation after 4-6 Gyr. [Abridged]

A defining prediction of the cold dark matter (CDM) cosmological model is the existence of a very large population of low mass haloes, down to planet-size masses. However, their fate as they are accreted onto haloes many orders of magnitude more massive remains fundamentally uncertain. A number of numerical explorations have found subhaloes to be very resilient to tides, but resolution limits make it difficult to explore tidal evolution at arbitrarily low masses. What are the structural properties of heavily stripped subhaloes? Do tidal effects destroy low-mass CDM subhaloes? Here we focus on cosmologically motivated subhaloes orbiting CDM hosts, and show that subhaloes of any initial mass can be stripped to arbitrarily small mass fractions. We show that previous numerical results can be reproduced by a simple model that describes tidal evolution as a progressive `peeling' in energy space and subsequent re-virialization. Tidal heating can effectively be neglected at all masses, because its importance i) does not increase for subhaloes of decreasing mass at accretion, and ii) it decreases as stripping proceeds. This allows us to predict analytically the structural properties of subhaloes with any initial mass and arbitrary degrees of mass loss. Under the hypotheses that CDM haloes have centrally divergent density profiles and approximately isotropic phase space distributions, our results prove that subhaloes of very low masses are a robust prediction of the CDM model: as for haloes, CDM subhalo populations extend and are abundant down to very low masses.

Core-collapse supernovae (CCSNe) produce large ($\gtrsim 0.1 \, {\rm M}_\odot$) masses of dust, and are potentially the primary source of dust in the Universe, but much of this dust may be destroyed before reaching the interstellar medium. Cassiopeia A (Cas A) is the only supernova remnant where an observational measurement of the dust destruction efficiency in the reverse shock is possible at present. We determine the pre- and post-shock dust masses in Cas A using a substantially improved dust emission model. In our preferred models, the unshocked ejecta contains $0.6-0.8 \, {\rm M}_\odot$ of $0.1 \, {\rm \mu m}$ silicate grains, while the post-shock ejecta has $0.02-0.09 \, {\rm M}_\odot$ of $5-10 {\, {\rm nm}}$ grains in dense clumps, and $2 \times 10^{-3} \, {\rm M}_\odot$ of $0.1 \, {\rm \mu m}$ grains in the diffuse X-ray emitting shocked ejecta. The implied dust destruction efficiency is $74-94 \%$ in the clumps and $92-98 \%$ overall, giving Cas A a final dust yield of $0.05-0.30 \, {\rm M}_\odot$. If the unshocked ejecta grains are larger than $0.1 \, {\rm \mu m}$, the dust masses are higher, the destruction efficiencies are lower, and the final yield may exceed $0.5 \, {\rm M}_\odot$. As Cas A has a dense circumstellar environment and thus a much stronger reverse shock than is typical, the average dust destruction efficiency across all CCSNe is likely to be lower, and the average dust yield higher. This supports a mostly-stellar origin for the cosmic dust budget.

Radio free-free emission is considered to be one of the most reliable tracers of star formation in galaxies. However, as it constitutes the faintest part of the radio spectrum -- being roughly an order of magnitude less luminous than radio synchrotron emission at the GHz frequencies typically targeted in radio surveys -- the usage of free-free emission as a star formation rate tracer has mostly remained limited to the local Universe. Here we perform a multi-frequency radio stacking analysis using deep Karl G. Jansky Very Large Array observations at 1.4, 3, 5, 10 and 34 GHz in the COSMOS and GOODS-North fields to probe free-free emission in typical galaxies at the peak of cosmic star formation. We find that $z \sim 0.5 - 3$ star-forming galaxies exhibit radio emission at rest-frame frequencies of $\sim 65 - 90$ GHz that is $\sim 1.5 - 2\times$ fainter than would be expected from a simple combination of free-free and synchrotron emission, as in the prototypical starburst galaxy M82. We interpret this as a deficit in high-frequency synchrotron emission, while the level of free-free emission is as expected from M82. We additionally provide the first constraints on the cosmic star formation history using free-free emission at $0.5 \lesssim z \lesssim 3$, which are in good agreement with more established tracers at high redshift. In the future, deep multi-frequency radio surveys will be crucial in order to accurately determine the shape of the radio spectrum of faint star-forming galaxies, and to further establish radio free-free emission as a tracer of high-redshift star formation.

Multi-band images of galaxies reveal a huge amount of information about their morphology and structure. However, inferring properties of the underlying stellar populations such as age, metallicity or kinematics from those images is notoriously difficult. Traditionally such information is best extracted from expensive spectroscopic observations. Here we present the $Painting\, IntrinsiC\, Attributes\, onto\, SDSS\, Objects$ (PICASSSO) project and test the information content of photometric multi-band images of galaxies. We train a convolutional neural network on 27,558 galaxy image pairs to establish a connection between broad-band images and the underlying physical stellar and gaseous galaxy property maps. We test our machine learning (ML) algorithm with SDSS $ugriz$ mock images for which uncertainties and systematics are exactly known. We show that multi-band galaxy images contain enough information to reconstruct 2d maps of stellar mass, metallicity, age and gas mass, metallicity as well as star formation rate. We recover the true stellar properties on a pixel by pixel basis with only little scatter, $\lesssim20\%$ compared to $\sim50\%$ statistical uncertainty from traditional mass-to-light-ratio based methods. We further test for any systematics of our algorithm with image resolution, training sample size or wavelength coverage. We find that galaxy morphology alone constrains stellar properties to better than $\sim20\%$ thus highlighting the benefits of including morphology into the parameter estimation. The machine learning approach can predict maps of high resolution, only limited by the resolution of the input bands, thus achieving higher resolution than IFU observations. The network architecture and all code is publicly available on GitHub.

We provide a highly concise overview of what X-ray surveys and their multiwavelength follow-up have revealed about the nature of the cosmic X-ray background (CXRB) and its constituent sources. We first describe early global studies of the CXRB, the development of imaging CXRB surveys, and the resolved CXRB fraction. Second, we detail the sources detected in CXRB surveys describing their identification, classification, and basic nature. Third, since active galactic nuclei (AGNs) are the main contributors to the CXRB, we discuss some key insights about their demographics, physics, and ecology that have come from CXRB surveys. Finally, we highlight future prospects for the field.

The abundance of the faintest galaxies provides insight into the nature of dark matter and the process of dwarf galaxy formation. In the LCDM scenario, low mass halos are so numerous that the efficiency of dwarf formation must decline sharply with decreasing halo mass in order to accommodate the relative scarcity of observed dwarfs and satellites in the Local Group. The nature of this decline contains important clues to the mechanisms regulating the onset of galaxy formation in the faintest systems. We explore here two possible models for the stellar mass ($M_*$)-halo mass ($M_{200}$) relation at the faint end, motivated by some of the latest LCDM cosmological hydrodynamical simulations. One model includes a sharp mass threshold below which no luminous galaxies form, as expected if galaxy formation proceeds only in systems above the Hydrogen-cooling limit. In the second model, $M_*$ scales as a steep power-law of $M_{200}$ with no explicit cutoff, as suggested by recent semianalytic work. Although both models predict satellite numbers around Milky Way-like galaxies consistent with current observations, they predict vastly different numbers of ultra-faint dwarfs and of satellites around isolated dwarf galaxies. Our results illustrate how the satellite mass function around dwarfs may be used to probe the $M_*$-$M_{200}$ relation at the faint end and to elucidate the mechanisms that determine which low-mass halos "light up" or remain dark in the LCDM scenario.

The Orion Nebula Cluster (ONC) is the most massive region of active star formation within a kpc of the Sun. Using Gaia EDR3 parallaxes and proper motions, we examine the bulk motions of stars radially and tangentially relative to the cluster center. We find an age gradient with distance to the stars in the ONC, from 385 pc for the oldest stars, to 395 pc for the younger stars, indicating that the star forming front is propagating into the cloud. We find an organized signature of rotation of the central cluster, but it is present only in stars younger than 2 Myr. We also observe a net infall of young stars into the center of the ONC's deep gravitational potential well. The infalling sources lie preferentially along the filament, on the other hand, outflowing sources are distributed spherically around the cluster, and they have larger velocity dispersion. We further propose a solution to a long-standing question of why the ONC shows a weak signature of expansion even though the cluster is likely bound: much of this expansion is likely driven by unstable N-body interactions among stars, resulting in low-velocity ejections. Finally we observe a significant infall of stars in various low-mass star-forming regions towards the Orion Complex at distances as far away as 200 pc, presumably due to a strong gravitational potential of Orion. Through analyzing signatures imprinted on stellar dynamics across different spatial scales, these observation shed new light on the signatures of formation and evolution of young clusters.

We introduce the ASTRID simulation, a large-scale cosmological hydrodynamic simulation in a $250$ Mpc/h box with $2\times 5500^3$ particles. ASTRID contains a large number of high redshift galaxies, which can be compared to future survey data, and resolves galaxies in halos more massive than $2\times 10^9 M_\odot$. ASTRID has been run from $z=99$ to $z=3$. As a particular focus is modelling the high redshift Universe, it contains models for inhomogeneous hydrogen and helium reionization, baryon relative velocities and massive neutrinos, as well as supernova and AGN feedback. The black hole model includes mergers driven by dynamical friction rather than repositioning. We briefly summarise the implemented models, and the technical choices we took when developing the simulation code. We validate the model, showing good agreement with observed UV luminosity functions, galaxy stellar mass functions and specific star formation rates. We show that the redshift at which a given galaxy underwent hydrogen reionization has a large effect on the halo gas fraction. Finally, at $z=6$, halos with $M \sim 2\times 10^9 M_\odot$ which have been reionized have a star formation rate $1.5$ times greater than those which have not yet been reionized.

To recognize gravitational wave lensing events and being able to differentiate between similar lens models will be of crucial importance once one will be observing several lensing events of gravitational waves per year. In this work, we study the lensing of gravitational waves in the context of LISA sources and wave-optics regime. While different papers before ours studied microlensing effects enhanced by simultaneous strong lensing, we focus on frequency (time) dependent phase effects produced by one lens that will be visible with only one lensed image. We will show how, in the interference regime (i.e. when interference patterns are present in the lensed image), we are able to i) distinguish a lensed waveform from an unlensed one, and ii) differentiate between different lens models. In pure wave-optics, on the other hand, the feasibility of the study depends on the SNR of the signal and/or the magnitude of the lensing effect. To achieve these goals we study the phase of the amplification factor of the different lens models and its effect on the unlensed waveform, and we exploit the signal-to-noise calculation for a qualitative analysis.

Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within $\lesssim 3\sigma$ of the mean of all simulations on a given grid ($128^3$ and $256^3$ grid cells), where $\sigma$ describes the statistical variation due to the flow's time dependence. They also agree well with a $512^3$ reference run. The $128^3$ and $256^3$ simulations agree within $9\%$ and $4\%$, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study.

We demonstrate the use of an eigenbasis that is derived from principal component analysis (PCA) applied on an ensemble of random-noise images that have a "red" power spectrum; i.e., a spectrum that decreases smoothly from large to small spatial scales. The pattern of the resulting eigenbasis allows for the reconstruction of images with a broad range of image morphologies. In particular, we show that this general eigen basis can be used to efficiently reconstruct images that resemble possible astronomical sources for interferometric observations; even though the images in the original ensemble used to generate the PCA basis are significantly different from the astronomical images. We further show that the efficiency and fidelity of the image reconstructions depends only weakly on the particular parameters of the red-noise power spectrum used to generate the ensemble of images.

The precision measurements of the spectra of cosmic ray nuclei and leptons in recent years have revealed the existence of multiple features, such as the spectral break at $\sim 300$ GV rigidity seen by PAMELA and AMS-02 and more recently confirmed by DAMPE and CALET, the softening in the spectra of H and He nuclei at $\sim 10$ TV reported by DAMPE, confirming previous hints by NUCLEON and CREAM, a tiny change of slope at $\sim 40$ GeV in the electron spectrum, revealed by AMS-02, and the large spectral break at $\sim$ TeV reported by indirect (HESS, MAGIC and VERITAS) and direct (DAMPE, CALET) measurements of the total (electrons+positrons) lepton spectrum. In all these cases, the possibility has been suggested that these features might reflect the occasional presence of a local cosmic ray source, inducing a noticeable reshaping of the average expected spectra. All these proposals have to face the question of how likely it is for such a source to exist, a question that we address here in a quantitative way. We study the statistical properties of random distribution of sources in space and time, and the effect of the spiral structure of our Galaxy for both the spectra of light nuclei (p and He) and leptons (electrons and positrons) in different energy regions.

We report spectroscopic confirmations of 15 Lyman-alpha galaxies at $z\sim7$, implying a spectroscopic confirmation rate of $\sim$80% on candidates selected from LAGER (Lyman-Alpha Galaxies in the Epoch of Reionization), which is the largest (24 deg$^2$) survey aimed at finding Lyman-alpha emitters (LAEs) at $z\sim7$ using deep narrow-band imaging from DECam at CTIO. LAEs at high-redshifts are sensitive probes of cosmic reionization and narrow-band imaging is a robust and effective method for selecting a large number of LAEs. In this work, we present results from the spectroscopic follow-up of LAE candidates in two LAGER fields, COSMOS and WIDE-12, using observations from Keck/LRIS. We report the successful detection of Ly$\alpha$ emission in 15 candidates (11 in COSMOS and 4 in WIDE-12 fields). Three of these in COSMOS have matching confirmations from a previous LAGER spectroscopic follow-up and are part of the overdense region, LAGER-$z7$OD1. Additionally, two candidates that were not detected in the LRIS observations have prior spectroscopic confirmations from Magellan. Including these, we obtain a spectroscopic confirmation success rate of $\sim$$80$% for LAGER LAE candidates. Apart from Ly$\alpha$, we do not detect any other UV nebular lines in our LRIS spectra; however, we estimate a 2$\sigma$ upper limit for the ratio of NV/Ly$\alpha$, $f_{NV}/f_{Ly\alpha} \lesssim 0.27$, which implies that ionizing emission from these sources is mostly dominated by star formation. Including confirmations from this work, a total of 33 LAE sources from LAGER are now spectroscopically confirmed. LAGER has more than doubled the sample of spectroscopically confirmed LAE sources at $z\sim7$.

Bubble chambers using liquid xenon (and liquid argon) have been operated (resp. planned) by the Scintillating Bubble Chamber (SBC) collaboration for GeV-scale dark matter searches and CE$\nu$NS from reactors. This will require a robust calibration program of the nucleation efficiency of low-energy nuclear recoils in these target media. Such a program has been carried out by the PICO collaboration, which aims to directly detect dark matter using $\mathrm{C_3 F_8}$ bubble chambers. Neutron calibration data from mono-energetic neutron beam and Am-Be source has been collected and analyzed, leading to a global fit of a generic nucleation efficiency model for carbon and fluorine recoils, at thermodynamic thresholds of $2.45$ and $3.29\,\mathrm{keV}$. Fitting the many-dimensional model to the data ($34$ free parameters) is a non-trivial computational challenge, addressed with a custom Markov Chain Monte Carlo approach, which will be presented. Parametric MC studies undertaken to validate this methodology are also discussed. This fit paradigm demonstrated for the PICO calibration will be applied to existing and future scintillating bubble chamber calibration data.

Precise measurements of black hole (BH) masses are essential to understanding the coevolution of these sources and their host galaxies. In this work, we develop a novel approach to compute BH virial masses using measurements of continuum luminosities and emission line widths from partially-overlapping, narrow-band observations of quasars; we refer to this technique as single-epoch photometry. This novel method relies on forward-modelling quasar observations to estimate the previous properties, which enables accurate measurements of emission line widths even for lines poorly resolved by narrow-band data. We assess the performance of this technique using quasars from the Sloan Digital Sky Survey (SDSS) observed by the miniJPAS survey, a proof-of-concept project of the J-PAS collaboration covering $\simeq1\,\mathrm{deg}^2$ of the northern sky using the 56 J-PAS narrow-band filters. We find remarkable agreement between BH masses from single-epoch SDSS spectra and single-epoch miniJPAS photometry, with no systematic difference between these and a scatter ranging from 0.4 to 0.07 dex for masses from $\log(M_\mathrm{BH}/\mathrm{M}_\odot)\simeq8$ to 9.75, respectively. Reverberation mapping studies show that single-epoch masses approximately present 0.4 dex precision, letting us conclude that our novel technique delivers BH masses with only mildly worse precision than single-epoch spectroscopy. The J-PAS survey will soon start observing thousands of square degrees without any source preselection other than the photometric depth in the detection band, and thus single-epoch photometry has the potential to provide details on the physical properties of quasar populations not satisfying the preselection criteria of previous spectroscopic surveys.

Next-generation surveys will provide photometric and spectroscopic data of millions to billions of galaxies with unprecedented precision. This offers a unique chance to improve our understanding of the galaxy evolution and the unresolved nature of dark matter (DM). At galaxy scales, the density distribution of DM is strongly affected by the astrophysical feedback processes, which are difficult to fully account for in classical techniques to derive mass models. In this work, we explore the capability of supervised learning algorithms to predict the DM content of galaxies from luminous observational-like parameters, using the public catalog of the TNG100 simulation. In particular, we use Photometric, Structural and Kinematic parameters to predict the total DM mass, DM half-mass radius, DM mass inside one and two stellar half-mass radii. We adopt the coefficient of determination, $R^2$, as a reference metric to evaluate the accuracy of these predictions. We find that the Photometric features alone are able to predict the total DM mass with fair accuracy, while Structural and Photometric features together are more effective to determine the DM inside the stellar half mass radius, and the DM within twice the stellar half mass radius. However, using all observational quantities together (Photometry, Structural and Kinematics) incredibly improves the overall accuracy for all DM quantities. This first test shows that Machine Learning tools are promising approaches to derive predictions of the DM in real galaxies. The next steps will be to improve observational realism of the training sets, by closely select samples which accurately reproduce the typical observed luminous scaling relations. The trained pipelines will be suitable for real galaxy data collected from the next-generation surveys like Rubin/LSST, Euclid, CSST, 4MOST, DESI, to derive, e.g., the properties of their central DM fractions.

The high-precision Juno gravitational measurements allow us to infer the structure of Jupiter's deep atmospheric zonal flow. Since this inference is nonunique, it is important to explore the space of possible solutions. In this paper, we consider a model in which Jupiter's deep atmospheric zonal flow is barotropic, or invariant along the direction of the rotation axis, until it is truncated at depth by some dynamical process (e.g., Reynolds stress, Lorentz or viscous force). We calculate the density perturbation produced by the $z$-invariant part of the flow using the thermal wind equation and compare the associated odd zonal gravitational harmonics ($J_{3}$, $J_{5}$, $J_{7}$, $J_{9}$) to the Juno-derived values. Most of the antisymmetric gravitational signal measured by Juno can be explained by extending observed winds between $20.9^{\circ}\rm{S}-26.4^{\circ}\rm{N}$ to depths of $\sim 1000$ km. Because the small-scale features of the mid/high latitude zonal flow may not persist to depth, we allow the zonal flow in this region to differ from the observed surface winds. We find that the Juno odd zonal gravitational harmonics can be fully explained by $\sim 1000$ km deep barotropic zonal flows involving the observed winds between $20.9^{\circ}\rm{S}-26.4^{\circ}$N and a few broad mid/high latitude jets.

Contribution submitted for inclusion in `Advances in Very High Energy Astrophysics`, Mukherjee & Zanin, World Scientific (2022) Cosmic gamma rays are a perfect probe to search for fundamental physics. Some of these exotic and exciting scenarios are the subject of this contribution to the book. All current Imaging Atmospheric Cherenkov Telescopes (IACTs) have invested a great deal of time and resources in scouting these out. In this chapter, we focus on the case of indirect search for weakly interacting massive particles (Sec. 8.1) and that of search for axion-like particles (Sec. 8.2). We continue with less debated, yet interesting studies on the search for primordial black holes (Sec. 8.3), tau-neutrinos (8.4), and magnetic monopoles (8.5). Other scenarios of fundamental physics like the search for Lorentz Invariance violations are treated elsewhere in the book. We structured the contribution providing a self-contained (yet minimal) theoretical framework, as well as a complete report of all IACTs' published contribution to the topics under consideration. Our aim is to provide a reference review as well as take a photograph of the effort of the current generation of IACTs and the challenges taken, in preparation of the next-generation instrument to come, the Cherenkov Telescope Array.

Gravitational weak lensing by dark matter halos leads to a measurable imprint in the shear correlation function of galaxies. Fuzzy dark matter (FDM), composed of ultralight axion-like particles of mass $m\sim 10^{-22}\text{ eV}$, suppresses the matter power spectrum and shear correlation with respect to standard cold dark matter. We model the effect of FDM on cosmic shear using the optimised halo model HMCode, accounting for additional suppression of the mass function and halo concentration in FDM as observed in $N$-body simulations. We combine Dark Energy Survey year 1 (DES-Y1) data with the Planck cosmic microwave background anisotropies to search for shear correlation suppression caused by FDM. We find no evidence of suppression compared to the preferred cold DM model, and thus set a new lower limit to the FDM particle mass. Using a log-flat prior and marginalising over uncertainties related to the non-linear model of FDM, we find a new, independent 95\% C.L. lower limit $\log_{10}m>-23$ combining Planck and DES-Y1 shear, an improvement of almost two orders of magnitude on the mass bound relative to CMB-only constraints. Our analysis is largely independent of baryonic modelling, and of previous limits to FDM covering this mass range. Our analysis highlights the most important aspects of the FDM non-linear model for future investigation. The limit to FDM from weak lensing could be improved by up to three orders of magnitude with $\mathcal{O}(0.1)$ arcmin cosmic shear angular resolution, if FDM and baryonic feedback can be simultaneously modelled to high precision in the halo model.

In data from the Kepler mission, the normal F3V star KIC 8462852 (Boyajian's star) was observed to exhibit infrequent dips in brightness that have not been satisfactorily explained. A previous paper reported the first results of a search for other similar stars in a limited region of the sky around the Kepler field. This paper expands on that search to cover the entire sky between declinations of +22 degrees and +68 degrees. Fifteen new candidates with low rates of dipping, referred to as "slow dippers" in Paper I, have been identified. The dippers occupy a limited region of the HR diagram and an apparent clustering in space is found. This latter feature suggests that these stars are attractive targets for SETI searches.

We use 3mm continuum NOEMA and NH$_3$ VLA observations towards the First Hydrostatic Core (FHSC) candidate CB 17 MMS to reveal the dust structure and gas properties down to 600-1,100 au scales and constrain its evolutionary stage. We do not detect any compact source at the previously identified 1.3 mm point source, despite expecting a minimum signal-to-noise of 9. The gas traced by NH$_3$ exhibits subsonic motions, with an average temperature of 10.4 K. A fit of the radial column density profile derived from the ammonia emission finds a flat inner region of radius $\sim$1,800 au and a central density of $\sim$6$\times10^5$ cm$^{-3}$. Virial and density structure analysis reveals the core is marginally bound ($\alpha_{vir}$= 0.73). The region is entirely consistent with a young starless core, hence ruling out CB 17 MMS as a FHSC candidate. Additionally, the core exhibits a velocity gradient aligned with the major axis, showing an arc-like structure in the p-v diagram and an off-center region with high-velocity dispersion caused by two distinct velocity peaks. These features could be due to interaction with the nearby outflow, which appears to deflect due to the dense gas near the NH$_3$ column density peak. We investigate the specific angular momentum profile of the starless core, finding that it aligns closely with previous studies of such radial profiles in Class 0 sources. This similarity to more evolved objects suggests that motions at 1,000 au scales are determined by large-scale dense cloud motions and may be preserved through the early stages of star formation.

We present optical UBVRI photometry and low-to-medium resolution spectroscopic observations of type Iax SN 2020sck spanning -5.5 d to +67 d from maximum light in the B-band. From the photometric analysis we find $\Delta m_{\rm{B}}$(15) = 2.03$\pm$0.05 mag and $M_{\rm{B}}$=-17.81$\pm$0.22 mag. Radiation diffusion model fit to the quasi-bolometric light curve indicates 0.13$\pm$0.02 $M_{\odot}$ of $^{56}$Ni and 0.34 $M_{\rm \odot}$ of ejecta are synthesized in the explosion. Comparing the observed quasi-bolometric light curve with angle-averaged bolometric light curve of three-dimensional pure deflagration explosion of $M_{\rm{ch}}$ carbon-oxygen white dwarf, we find agreement with a model in which 0.16 $M_{\odot}$ of $^{56}$Ni and 0.37 $M_{\odot}$ of ejecta is formed. By comparing the +1.4 day spectrum of SN 2020sck with synthetic spectrum generated using SYN++, we find absorption features due to C II, C III and O I. These are unburned materials in the explosion and indicate a C-O white dwarf. One dimensional radiative transfer modeling of the spectra with TARDIS shows higher density in the ejecta near the photosphere and a steep decrease in the outer layers with an ejecta composition dominated mostly by C, O, Si, Fe, and Ni. The star formation rate of the host galaxy computed from the luminosity of the H$\alpha$ ($\lambda$6563) line is 0.09 $M_{\odot}$ yr$^{-1}$ indicating a relatively young stellar environment.

We present the description of the project \texttt{SCORPIO}, a Python package for retrieving images and associated data of galaxy pairs based on their position, facilitating visual analysis and data collation of multiple archetypal systems. The code ingests information from SDSS, 2MASS, and WISE surveys based on the available bands and is designed for studies of galaxy pairs as natural laboratories of multiple astrophysical phenomena such as tidal force deformation of galaxies, pressure gradient induced star formation regions, morphological transformation, to name a few.

Odd Radio Circles (ORCs) are unexpected faint circles of diffuse radio emission discovered in recent wide deep radio surveys. They are typically about one arcmin in diameter, and may be spherical shells of synchrotron emission about a million light years in diameter, surrounding galaxies at a redshift of ~0.2-0.6. Here we study the properties and environment of the known ORCs. All three known single ORCs either lie in a significant overdensity or have a close companion. If the ORC is caused by an event in the host galaxy, then the fact that they tend to be in an overdensity, or have a close companion, may indicate that the environment is important in creating the ORC phenomenon, possibly because of an increased ambient density or magnetic field.

High velocity clouds moving toward the disk will reach the Galactic plane and will inevitably collide with the disk. In these collisions a system of two shocks is produced, one propagating through the disk and the other develops within the cloud. The shocks produced within the clouds in these interactions have velocities of hundreds of kilometers per second. When these shocks are radiative they may be inefficient in accelerating fresh particles, however they can reaccelerate and compress Galactic cosmic rays from the background. In this work we investigate the interactions of Galactic cosmic rays within a shocked high velocity cloud, when the shock is induced by the collision with the disk. This study is focused in the case of radiative shocks. We aim to establish under which conditions these interactions lead to significant nonthermal emission, especially gamma rays. We model the interaction of cosmic ray protons and electrons reaccelerated and further energized by compression in shocks within the clouds, under very general assumptions. We also consider secondary electron-positron pairs produced by the cosmic ray protons when colliding with the material of the cloud. We conclude that nearby clouds reaccelerating Galactic cosmic rays in local shocks can produce high-energy radiation that might be detectable with existing and future gamma-ray detectors. The emission produced by electrons and secondary pairs is important at radio wavelengths, and in some cases it may be relevant at hard X-rays. Concerning higher energies, the leptonic contribution to the spectral energy distribution is significant at soft gamma rays.

We investigate the dynamics of interstellar dust particles in moderately high resolution ($512^3$ grid points) simulations of forced compressible transonic turbulence including self-gravity of the gas. Turbulence is induced by stochastic compressive forcing which is delta-correlated in time. By considering the nearly Jeans-unstable case, where the scaling of the simulation is such that a statistical steady state without any irreversible collapses is obtained, we obtain a randomly varying potential, acting as a second stochastic forcing. We show that, in this setting, low-inertia grains follow the gas flow and cluster in much the same way as in a case of statistical steady-state turbulence without self-gravity. Large, high-inertia grains, however, are accelerated to much higher mean velocities in the presence of self-gravity. Grains of intermediate size also show an increased degree of clustering. We conclude that self-gravity effects can play an important role for aggregation/coagulation of dust even in a turbulent system which is not Jeans-unstable. In particular, the collision rate of large grains in the interstellar medium can be much higher than predicted by previous work.

We observed ASASSN-19ax during the long outburst in 2021 September-October. The object has been confirmed to be an SU UMa-type dwarf nova with a superhump period of 0.1000-0.1001 d. This object showed two post-superoutburst rebrightenings both in the 2019 and 2021 superoutbursts. These observations have established that ASASSN-19ax belongs to a group of long-period SU UMa-type dwarf novae which show multiple rebrightenings. This phenomenon probably arises from premature quenching of the superoutburst due to the weak 3:1 resonance near the stability border of the resonance, resulting in a considerable amount of disk mass after the superoutburst. We noted that ASASSN-19ax is very similar to QZ Ser, an SU UMa-type dwarf nova with an orbital period of 0.08316 d and an anomalously hot, bright secondary star, in that both objects showed multiple post-superoutburst rebrightenings at least once and that they are bright in quiescence. We expect that the core of the secondary in ASASSN-19ax may be evolved as in QZ Ser.

We investigate differences in Spitzer/IRAC 3.6 and 4.5micron photometry that depend on observing strategy. Using archival calibration data we perform an in-depth examination of the measured flux densities ("fluxes") of ten calibration stars, observed with all the possible observing strategies. We then quantify differences in the measured fluxes as a function of 1) array mode (full or subarray), 2) exposure time, and 3) dithering versus staring observations. We find that the median fluxes measured for sources observed using the full array are 1.6% and 1% lower than those observed with the subarray at [3.6] and [4.5], respectively. Additionally, we found a dependence on the exposure time such that for [3.6] observations the long frame times are measured to be lower than the short frame times by a median value of 3.4% in full array and 2.9% in subarray. For [4.5] observations the longer frame times are 0.6% and 1.5% in full and subarray respectively. These very small variations will likely only affect science users who require high-precision photometry from multiple different observing modes. We find no statistically significant difference for fluxes obtained with dithered and staring-modes. When considering all stars in the sample, the fractional well depth of the pixel is correlated with the different observed fluxes. We speculate the cause to be a small non-linearity in the pixels at the lowest well depths where deviations from linearity were previously assumed to be negligible.

We present new stellar evolutionary sequences of very metal-rich stars evolved with the Monash Stellar Structure code and with MESA. The Monash models include masses of $1-8M_{\odot}$ with metallicities $Z=0.04$ to $Z=0.1$ and are evolved from the main sequence to the thermally-pulsing asymptotic giant branch (AGB). These are the first $Z=0.1$ AGB models in the literature. The MESA models include intermediate-mass models with $Z=0.06$ to $Z=0.09$ evolved to the onset of the thermally-pulsing phase. Third dredge-up only occurs in intermediate-mass models $Z \le 0.08$. Hot bottom burning (HBB) shows a weaker dependence on metallicity, with the minimum mass increasing from 4.5$M_{\odot}$ for $Z=0.014$ to $\approx 5.5 M_{\odot}$ for Z = 0.04, $6M_{\odot}$ for $ 0.05 \le Z \le 0.07$ and above 6.5$ M_{\odot}$ for $Z\ge 0.08$. The behaviour of the $Z=0.1$ models is unusual; most do not experience He-shell instabilities owing to rapid mass-loss on the early part of the AGB. Turning off mass-loss produces He-shell instabilities, however thermal pulses are weak and result in no third dredge-up. The minimum mass for carbon ignition is reduced from 8$M_{\odot}$ for $Z=0.04$ to 7$M_{\odot}$ for $Z=0.1$, which implies a reduction in the minimum mass for core-collapse supernovae. MESA models of similarly high metallicity ($Z=0.06 - 0.09$) show the same lowering of the minimum mass for carbon ignition: carbon burning is detected in a $6 M_{\odot}$ model at the highest metallicity ($Z=0.09$) and in all $7 M_{\odot}$ models with $Z \ge 0.06$. This demonstrates robustness of the lowered carbon burning threshold across codes.

The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure phase measurements on all six linearly independent non-trivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of $\sim3-5^\circ$. The only triangles that exhibit substantially higher variability ($\sim90-180^\circ$) are the ones with baselines that cross visibility amplitude minima on the $u-v$ plane, as expected from theoretical modeling. We used two sets of General Relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black-hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black-hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas.

Non-canonical cosmology with an uplifted Higgs vacuum expectation value (Higgs-VEV) in the early universe is believed to provide the solution for existing tensions within the $\Lambda$CDM regime. We recently proposed a theoretical model called axi-Higgs to explore this framework. The axi-Higgs model features an ultralight axion with mass $m_a \sim 10^{-29}$ eV, which couples to the Higgs field such that the Higgs-VEV is driven by the axion background evolution. In this paper, we perform Markov Chain Monte Carlo (MCMC) analyses with our modified Boltzmann solver to investigate the parameter space of axi-Higgs beside other models including $\Lambda$CDM, $\Lambda$CDM+$m_e$, $\Lambda$CDM+$\omega_a$. Combining cosmological data from Cosmic Microwave Background (CMB), Baryon Acoustic Oscillations (BAO), weak-lensing (WL) cosmic shear survey, we found $H_0 = 69.3^{+1.2}_{-1.4}$ km/s/Mpc and $S_8 = 0.797 \pm 0.012$, which reduces the Hubble tension to approximately $2.5\sigma$ and slightly alleviates the $S_8$ tension. The presence and behavior of this Higgs-VEV-driving axion may be tested by atomic clock measurements in laboratory and/or quasar spectral measurements in the near future.

We present line and continuum observations (resolution ~0.3"-3.5") made with the Atacama Large Millimeter/submillimeter Array (ALMA), Submillimeter Array, and Very Large Array of a young O-type protostar W42-MME (mass: 19-4 Msun). The ALMA 1.35 mm continuum map (resolution ~1") shows that W42-MME is embedded in one of the cores (i.e., MM1) located within a thermally supercritical filament-like feature (extent ~0.15 pc) containing three cores (mass ~1-4.4 Msun). Several dense/hot gas tracers are detected toward MM1, suggesting the presence of a hot molecular core with the gas temperature of ~38-220~K. The ALMA 865 micron continuum map (resolution ~0.3") reveals at least five continuum sources/peaks ("A-E") within a dusty envelope (extent ~9000 AU) toward MM1, where shocks are traced in the SiO(8-7) emission. The source "A" associated with W42-MME is seen almost at the center of the dusty envelope, and is surrounded by other continuum peaks. The ALMA CO(3-2) and SiO(8-7) line observations show the bipolar outflow extended below 10000 AU, which is driven by the source "A". The ALMA data hint the episodic ejections from W42-MME. A disk-like feature (extent ~2000 AU; mass ~1 Msun) with velocity gradients is investigated in the source "A" (dynamical mass ~9 Msun) using the ALMA H13CO+ emission, and is perpendicular to the CO outflow. A small-scale feature (below 3000 AU) probably heated by UV radiation from the O-type star is also investigated toward the source "A". Overall, W42-MME appears to gain mass from its disk and the dusty envelope.

The MAMMOTH-1 nebula at $z=2.317$ is an enormous Ly$\alpha$ nebula (ELAN) extending to a $\sim$440 kpc scale at the center of the extreme galaxy overdensity BOSS 1441. In this paper, we present observations of the $\rm CO(3-2)$ and 250 GHz dust-continuum emission from the MAMMOTH-1 using the IRAM NOrthern Extended Millimeter Array. Our observations show that $\rm CO(3-2)$ emission in this ELAN has not extended widespread emission into the circum- and inter-galactic media. We also find a remarkable concentration of six massive galaxies in $\rm CO(3-2)$ emission in the central $\sim$100 kpc region of the ELAN. Their velocity dispersions suggest a total halo mass of $M_{200c} \sim 10^{13.1} M_{\odot}$, marking a possible protocluster core associated with the ELAN. The peak position of the $\rm CO(3-2)$ line emission from the obscured AGN is consistent with the location of the intensity peak of MAMMOTH-1 in the rest-frame UV band. Its luminosity line ratio between the $\rm CO(3-2)$ and $\rm CO(1-0)$ $r_{3,1}$ is 0.61$\pm$0.17. The other five galaxies have $\rm CO(3-2)$ luminosities in the range of (2.1-7.1)$\times 10^9$ K $\rm km\,s^{-1}$ pc$^2$, with the star-formation rates derived from the 250GHz continuum of ($<$36)-224 $M_{\odot}$yr$^{-1}$. Follow-up spectroscopic observations will further confirm more member galaxies and improve the accuracy of the halo mass estimation.

We show that the star formation history, the reionization history and the luminosity function of galaxies are successfully predicted in a simple gravitational collapse model within the $\Lambda$CDM regime to almost a quantitative accuracy, when the physical conditions, the Jeans criterion and the cooling process, are taken into account: it predicts that the reionisation takes place sharply at around redshift $1+z\simeq 7.5$, and the resulting luminisity function turns off at $L_B\simeq 10^{10.7}L_\odot$, showing the consistency between the star formation history and reionisation of the Universe. The total amount of stars is predicted to be $\Omega_\mathrm{star}=0.004$ in units of the critical density compared to the observation $0.0044$ with the recycling factor $1.6$ included. The additional assumption needed for successful prediction is that the star formation efficiency is not halo mass independent but becomes maximum at the halo mass $\simeq 10^{12}M_\odot$ and is suppressed for both smaller and larger mass haloes. It is underlied that there are no adjustable parameters other than the suppression of the star formation rate in the model.

We investigate the multi-wavelength emission from hadronic and leptonic cosmic rays (CRs) in bubbles around galaxies, analogous to the Fermi bubbles of the Milky Way. The bubbles are modeled using 3D magnetohydrodynamical (MHD) simulations, and are driven by a 0.3 Myr intense explosive outburst from the nucleus of Milky Way-like galaxies. We compute their non-thermal emission properties at different stages throughout their evolution, up to 7 Myr, by post-processing the simulations. We compare the spectral and spatial signatures of bubbles with hadronic, leptonic and hybrid hadro-leptonic CR compositions. These each show broadly similar emission spectra, comprised of radio synchrotron, inverse Compton and bremsstrahlung components. However, hadronic and hybrid bubbles were found to be slightly dimmer than their leptonic counterparts, with a large part of their emission being driven by secondary electrons formed in hadronic interactions. The hadronic and hybrid bubbles also showed an additional $\pi^0$ decay emission component, which dominated their $\gamma$-ray emission from GeV energies to 10s of TeV. Conversely, we found that high-energy $\gamma$-ray emission from leptonic bubbles was inverse Compton-dominated, which fades quickly as the bubble expands and ages. This difference could be used to distinguish hadronic and leptonic bubbles around other galaxies, with older hadronic or hybrid bubbles being detectable in $\gamma$-rays to greater distances.

Thermoluminescence (TL) research provides a powerful tool for characterizing radiation-induced processes in extraterrestrial matter. One of the challenges in studying the spectral features of the natural TL of stony meteorites is its weak intensity. The present work showcases the capabilities of a high-sensitive original module for measuring the spectrally resolved TL characteristics of the Chelyabinsk and Tsarev chondrites. We have analyzed the emission spectra and glow curves of natural and induced TL over the 300 - 650 nm and RT - 873 K ranges. A quasi-continuous distribution of traps active within the 350 - 650 K range was found in the silicate substructure of both meteorites under study. Based on the general order kinetic formalism and using the natural TL data, we also estimated the activation energies of 0.86 and 1.08 eV for the Chelyabinsk and Tsarev chondrites, respectively.

Several methods are used to evaluate, from observational data, the dynamical state of galaxy clusters. Among them, the morphological analysis of cluster images is well suited for this purpose. We report a new approach to the morphology, which consists in analytically modelling the images with a set of orthogonal functions, the Zernike polynomials (ZPs). We validated the method on mock high-resolution Compton parameter maps of synthetic galaxy clusters from THE THREE HUNDRED project. To classify the maps for their morphology we defined a single parameter, $\mathcal{C}$, by combining the contribution of some ZPs in the modelling. We verify that $\mathcal{C}$ is linearly correlated with a combination of common morphological parameters and also with a proper 3D dynamical-state indicator available for the synthetic clusters we used. We also show the early results of the Zernike modelling applied on Compton parameter maps of local clusters ($z < 0.1$) observed by the $\textit{Planck}$ satellite. At last, we report the preliminary results of this kind of morphological analysis on mock X-ray maps of THE THREE HUNDRED clusters.

There is a growing interest of utilizing intrinsic alignment (IA) of galaxy shapes as a geometric and dynamical probe of cosmology. In this paper we present the first measurements of IA in a modified gravity model using the gravitational shear-intrinsic ellipticity correlation (GI) and intrinsic ellipticity-ellipticity correlation (II) functions of dark-matter halos from $f(R)$ gravity simulations. By comparing them with the same statistics measured in $\Lambda$CDM simulations, we find that the IA statistics in different gravity models show distinguishable features, with a trend similar to the case of conventional galaxy clustering statistics. Thus, the GI and II correlations are found to be useful in distinguishing between the $\Lambda$CDM and $f(R)$ gravity models. More quantitatively, IA statistics enhance detectability of the imprint of f(R) gravity on large scale structures by $\sim20\%$ when combined with the conventional halo clustering in redshift space. Our results demonstrate the usefulness of IA statistics as a probe of gravity beyond a consistency test of $\Lambda$CDM and general relativity.

The mm-to-cm range of the Spectral Energy Distribution of spiral galaxies remains largely unexplored. Its coverage is required to disentangle the contribution of dust emission, free-free and synchrotron radiation and can provide constraints on dust models, star-formation rates and ISM properties. We present the case for a synergy between NIKA2 observations of nearby spirals and those from planned and current instrumentation at the Sardinia Radio Telescope, and report on a pilot K-band program to search for Anomalous Microwave Emission, an elusive emission component which is presumably related to dust.

We present preliminary results from an on-going program that aims at mapping the intracluster medium (ICM) temperature of high redshift galaxy clusters from the MaDCoWS sample using a joint analysis of shallow X-ray data obtained by $Chandra$ and high angular resolution Sunyaev-Zel'dovich (SZ) observations realized with the NIKA2 and MUSTANG-2 cameras. We also present preliminary results from an on-going Open Time program within the NIKA2 collaboration that aims at mapping the ICM temperature of a galaxy cluster at $z=0.45$ from the resolved detection of the relativistic corrections to the SZ spectrum. These studies demonstrate how high angular resolution SZ observations will play a major role in the coming decade to push the investigation of ICM dynamics and non-gravitational processes to high redshift before the next generation X-ray observatories come into play.

We search for H$\alpha$ emitters at $z\sim7.8$ in four gravitationally lensed fields observed in the Hubble Frontier Fields program. We use the Lyman break method to select galaxies at the target redshift, and make the photometry in {\it Spitzer}/IRAC 5.8 $\mu$m band to detect the H$\alpha$ emission from the candidate galaxies. We find no significant detection of counterparts in the IRAC 5.8 $\mu$m band, and this gives a constraint on the H$\alpha$ luminosity function (LF) at $z\sim7.8$. We compare the constraint with previous studies on rest-frame UV and FIR observation using the correlation between the H$\alpha$ luminosity and the star formation rate. Additionally, we convert the constraint on the H$\alpha$ LF into an upper limit for the star formation rate density (SFRD) at this epoch assuming the shape of the LF. We examine two types of parameterization of the LF, and obtain an upper limit for the SFRD of $\log_{10}(\rho_{\rm SFR}\ [M_\odot\ \mathrm{yr^{-1}\ Mpc^{-3}}])\lesssim-1.1$ at $z\sim7.8$. With this constraint on the SFRD, we give an independent probe into the total star formation activity including the dust-obscured and unobscured star formation at the Epoch of Reionization.

We describe the preliminary on-sky results of the KIDs Interferometer Spectrum Survey (KISS), a spectral imager with a 1 deg field of view (FoV). The instrument operates in the range 120-180 GHz from the 2.25 m Q-U-I JOint TEnerife telescope in Teide Observatory (Tenerife, Canary Islands), at 2 395 m altitude above sea level. Spectra at low resolution, up to 1.45 GHz, are obtained using a fast (3.72 Hz mechanical frequency) Fourier transform spectrometer, coupled to a continuous dilution cryostat with a stabilized temperature of 170 mK that hosts two 316-pixel arrays of lumped-element kinetic inductance detectors. KISS generates more than 3 000 spectra per second during observations and represents a pathfinder to demonstrate the potential for spectral mapping with large FoV. We give an overall description of the spectral mapping paradigm and we present recent results from observations, in this paper.

In this talk I focus on how the modelling of the mass-observable relation and the halo mass function can impact the accuracy and precision of cosmological constraints inferred from galaxy clusters. I present a new analysis of clusters detected in mm wavelengths by the Planck satellite, highlighting the need of an improved description and calibration for the mass-observable relation.

We report the search for 7Be isotope in the outbursts of the classical nova V6595 Sgr by means of high resolution UVES observations taken at the ESO VLT in April 2021, about two weeks after discovery and under difficult circumstances due to the pandemic. Narrow absorption components with velocities at about -2620 and -2820 km/s, superposed on broader and shallow absorption, are observed in the outburst spectra for the 7BeII 313.0583, 313.1228 nm doublet resonance lines, as well as in several other elements such as CaII, FeI, MgI, NaI, HI but LiI. Using CaII K line as a reference element, we infer N(7Be)/N(H) ~ 7.4 x 10^{-6}, or ~ 9.8 x 10^{-6} when the 7Be decay is taken into account. The 7Be abundance is about half of the value most frequently measured in novae. The possible presence of over-ionization in the layers where 7Be is detected is also discussed. Observations taken at the Telescopio Nazionale Galileo (TNG) in La Palma 91 days after discovery showed prominent emission lines of Oxygen and Neon which allow to classify the nova as ONe type. Therefore, although 7Be is expected to be higher in CO novae, it is found at comparable levels in both nova types.

While third-generation CMB experiments have allowed to release the first maps of Compton-$y$ distortion due to thermal Sunyaev-Zeldovich (SZ) effect, next-generation CMB experiments should allow us to map also the electron gas temperature, $T_{\rm e}$, across the sky through the detection of relativistic corrections to the thermal SZ effect. We discuss about experimental requirements to break the $y$-$T_{\rm e}$ degeneracy of the observed SZ intensity, and propose a new component separation approach based on moment expansion to disentangle the $y$ and $T_{\rm e}$ observables of the relativistic SZ effect while mitigating foregrounds. We show how our approach offers a new spectroscopic view of the clusters not only across frequencies but now also across temperatures. We also show how the relativistic electron temperature power spectrum provides a new cosmological observable which may complement the Compton-$y$ map power spectrum to break some of the parameter degeneracies in future cosmological SZ analyses.

Icy pebbles may play an important role in planet formation close to the water ice line of protoplanetary discs. There, dust coagulation is more efficient and re-condensation of vapor on pebbles may enhance their growth outside the ice line. Previous theoretical studies showed that disruption of icy pebbles due to sublimation increases the growth rate of pebbles inside and outside the ice line, by freeing small silicate particles back in the dust reservoir of the disc. However, since planet accretion is dependent on the Stokes number of the accreting pebbles, the growth of planetesimals could be enhanced downstream of the ice line if pebbles are not disrupting upon sublimation. We developed two experimental models of icy pebbles using different silicate dusts, and we exposed them to low-temperature and low-pressure conditions in a vacuum chamber. Increasing the temperature inside the chamber, we studied the conditions for which pebbles are preserved through sublimation without disrupting. We find that small silicate particles ($<50$um) and a small quantity of ice (around $15$% pebble mass) are optimal conditions for preserving pebbles through sublimation. Furthermore, pebbles with coarse dust distribution ($100-300$um) do not disrupt if a small percentage ($10-20$% mass) of dust grains are smaller than $50$um. Our findings highlight how sublimation is not necessarily causing disruption, and that pebbles seem to survive fast sublimation processes effectively.

Primordial magnetic fields are generated during inflation by considering actions that break the conformal invariance of the electromagnetic field. To break the conformal invariance, the electromagnetic fields are coupled either to the inflaton or to the scalar curvature. Also, a parity violating term is often added to the action in order to enhance the amplitudes of the primordial electromagnetic fields. In this work, we examine the effects of deviations from slow roll inflation on the spectra of non-helical as well as helical electromagnetic fields. We find that, in the case of the coupling to the scalar curvature, there arise certain challenges in generating electromagnetic fields of the desired shapes and strengths even in slow roll inflation. When the field is coupled to the inflaton, it is possible to construct model-dependent coupling functions which lead to nearly scale invariant magnetic fields in slow roll inflation. However, we show that sharp features in the scalar power spectrum generated due to departures from slow roll inflation inevitably lead to strong features in the power spectra of the electromagnetic fields. Moreover, we find that such effects can also considerably suppress the strengths of the generated electromagnetic fields over the scales of cosmological interest. We illustrate these aspects with the aid of specific inflationary models that have been considered to produce specific features in the scalar power spectrum. Further, we find that, in such situations, if the strong features in the electromagnetic power spectra are to be undone, the choice of the coupling function requires considerable fine tuning. We discuss wider implications of the results we obtain.

Context. Serendipitous X-ray surveys have proven to be an efficient way to find rare objects, for example tidal disruption events, changing-look active galactic nuclei (AGN), binary quasars, ultraluminous X-ray sources (ULXs), and intermediate mass black holes. With the advent of very large X-ray surveys, an automated classification of X-ray sources becomes increasingly valuable.Aims. This work proposes a revisited naive Bayes classification of the X-ray sources in the Swift-XRT and XMM-Newton catalogs into four classes -- AGN, stars, X-ray binaries (XRBs) and cataclysmic variables (CVs) -- based on their spatial, spectral and timing properties and their multiwavelength counterparts. An outlier measure is used to identify objects of other natures. The classifier is optimized to maximize the classification performance of a chosen class (here XRBs) and it is adapted to data mining purposes.Methods. We augmented the X-ray catalogs with multiwavelength data, source class, and variability properties. We then built a reference sample of about 25000 X-ray sources of known nature. From this sample the distribution of each property was carefully estimated and taken as reference to assign probabilities of belonging to each class. The classification was then performed on the whole catalog, combining the information from each property.Results. Using the algorithm on the Swift reference sample we retrieved 99%, 98%, 92% and 34% of AGN, stars, XRBs, and CVs, respectively, and the false positive rates are 3%, 1%, 9% and 15%. Similar results are obtained on XMM sources. When applied to a carefully selected test sample, representing 55% of the X-ray catalog, the classification gives consistent results in terms of distributions of source properties. A substantial fraction of sources not belonging to any class is efficiently retrieved using the outlier measure, as well as AGN and stars with properties deviating from the bulk of their class. Our algorithm is then compared to a random forest method; the two showed similar performances but the algorithm presented in this paper improved insight into the grounds of each classification.Conclusions. This robust classification method can be tailored to include additional or different source classes and can be applied to other X-ray catalogs. The transparency of the classification compared to other methods makes it a useful tool in the search for homogeneous populations or rare source types, including multi-messenger events. Such a tool will be increasingly valuable with the development of surveys of unprecedented size, such as LSST, SKA and Athena, and the search for counterparts of multi-messenger events.

The gas mass fraction in galaxy clusters is a convenient tool to use in the context of cosmological studies. Indeed this quantity allows to constrain the universal baryon fraction $\Omega_b/\Omega_m$, as well as other parameters like the matter density $\Omega_m$, the Hubble parameter $h$ or the Equation of State of Dark Energy $w$. This gas mass fraction is also sensitive to baryonic effects that need to be taken into account, and that translate into nuisance parameters. Two of them are the depletion factor $\Upsilon$, and the hydrostatic mass bias $B = (1 - b)$. The first one describes how baryons are depleted in clusters with respect to the universal baryon fraction, while the other encodes the bias coming from the fact that the mass is deduced from X-ray observations under the hypothesis of hydrostatic equilibrium. We will show preliminary results, obtained using the {\it Planck}-ESZ clusters observed by XMM-{\it Newton}, on both cosmological and cluster parameters. We will notably discuss the investigation on a possible redshift dependence of the mass bias, which is considered to be non-existent in hydrodynamic simulations based on $\Lambda$-CDM, and compare our results with other studies.

The precise astrometric observation of small near-Earth objects (NEOs) is an important observational research topic in the astrometric discipline, which greatly promotes multidisciplinary research, such as the origin and evolution of the solar system, the detection and early warning of small NEOs, and deep-space navigation. The characteristics of small NEOs, such as faintness and fast moving speed, restrict the accuracy and precision of their astrometric observations. In the paper, we present a method to improve the accurate and precise astrometric positions of NEOs based on image fusion technique. The noise analysis and astrometric test from the observed images of the open cluster M23 are given. Using the image fusion technique, we obtain the sets of superimposed images and original images containing reference stars and moving targets respectively. The final fused image set includes background stars with high signal-to-noise ratios and ideal NEO images simultaneously and avoids the saturation of background stars. Using the fused images, we can reduce the influence of telescope tracking and NEO ephemeris errors on astrometric observations, and our results indicate that the accuracy and precision of NEO Eros astrometry are improved obviously after we choose suitable image fuse mode.

The CaFe Project involves the study of the properties of the low ionization emission lines (LILs) pertaining to the broad-line region (BLR) in active galaxies. These emission lines, especially the singly-ionized iron (Fe II) in the optical and the corresponding singly-ionized calcium (Ca II) in the near-infrared (NIR) are found to show a strong correlation in their emission strengths, i.e. with respect to the broad H$\beta$ emission line, the latter also belonging to the same category of LILs. The origin of this correlation is attributed to the similarity in the physical conditions necessary to emit these lines - especially in terms of the strength of the ionization from the central continuum source and the local number density of available matter in these regions. In this paper, we focus on the issue of the spectral energy distribution (SED) characteristic to a prototypical Type-1 Narrow-line Seyfert galaxy (NLS1) - I Zw 1. We extract the continuum from quasi-simultaneous spectroscopic measurements ranging from the near-UV ($\sim$1200A) to the near-infrared ($\sim$24000A) to construct the SED and supplement it with archival X-ray measurements available for this source. Using photoionization code CLOUDY, we assess and compare the contribution of the prominent "Big Blue Bump" seen in our SED versus the SED used in our previous work, wherein the latter was constructed from archival, multi-epoch photometric measurements. Following the prescription from our previous work, we constrain the physical parameter space to optimize the emission from these LILs and discuss the implication of the use of a "better" SED.

Warps are observed in a large fraction of disc galaxies, and can be due to a large number of different processes. Some of these processes might also cause vertical heating and flaring. Using a sample of galaxies simulated in their cosmological context, we study the connection between warping and disc heating. We analyse the vertical stellar density structure within warped stellar discs, and monitor the evolution of the scale-heights of the mono-age populations and the geometrical thin and thick disc during the warp's lifetime. We also compare the overall thickness and the vertical velocity dispersion in the disc before and after the warp. We find that for warps made of pre-existing stellar particles shifted off-plane, the scale-heights do not change within the disc's warped region: discs tilt rigidly. For warps made of off-plane new stellar material (either born in-situ or accreted), the warped region of the disc is not well described by a double $\mathrm{sech^2}$ density profile. Yet, once the warp is gone, the thin and thick disc structure is recovered, with their scale-heights following the same trends as in the region that was never warped. Finally, we find that the overall thickness and vertical velocity dispersion do not increase during a warp, regardless of the warp's origin. This holds even for warps triggered by interactions with satellites, which cause disc heating but before the warp forms. Our findings suggest that the vertical structure of galaxies does not hold any memory of past warps.

There are several strong K I lines found in the spectra of M dwarfs, among them the doublet near 7700 AA and another doublet near 12 500 AA. We study these optical and near-infrared doublets in a sample of 324 M dwarfs, observed with CARMENES, the high-resolution optical and near-infrared spectrograph at Calar Alto, and investigate how well the lines can be used as photospheric and chromospheric diagnostics. Both doublets have a dominant photospheric component in inactive stars and can be used as tracers of effective temperature and gravity. For variability studies using the optical doublet, we concentrate on the red line component because this is less prone to artefacts from telluric correction in individual spectra. The optical doublet lines are sensitive to activity, especially for M dwarfs later than M5.0 V where the lines develop an emission core. For earlier type M dwarfs, the red component of the optical doublet lines is also correlated with H$\alpha$ activity. We usually find positive correlation for stars with H$\alpha$ in emission, while early-type M stars with H$\alpha$ in absorption show anti-correlation. During flares, the optical doublet lines can exhibit strong fill-in or emission cores for our latest spectral types. On the other hand, the near-infrared doublet lines very rarely show correlation or anti-correlation to H$\alpha$ and do not change line shape significantly even during the strongest observed flares. Nevertheless, the near-infrared doublet lines show notable resolved Zeeman splitting for about 20 active stars which allows to estimate the magnetic fields B.

We use deep Hubble Space Telescope imaging to derive a distance to the Virgo Cluster ultradiffuse galaxy (UDG) VCC 615 using the tip of the red giant branch (TRGB) distance estimator. We detect 5,023 stars within the galaxy, down to a 50% completeness limit of F814W = 28.0, using counts in the surrounding field to correct for contamination due to background sources and Virgo intracluster stars. We derive an extinction-corrected F814W tip magnitude of $m_{\rm tip,0} = 27.19^{+0.07}_{-0.05}$, yielding a distance of $d=17.7^{+0.6}_{-0.4}$ Mpc. This places VCC 615 on the far side of the Virgo Cluster ($d_{\rm Virgo} = 16.5 Mpc$), at a Virgocentric distance of 1.3 Mpc and near the virial radius of the main body of Virgo. Coupling this distance with the galaxy's observed radial velocity, we find that VCC 615 is on an outbound trajectory, having survived a recent passage through the inner parts of the cluster. Indeed, our orbit modeling gives a 50% chance the galaxy passed inside the Virgo core (r<620 kpc) within the past Gyr, although very close passages directly through the cluster center (r<200 kpc) are unlikely. Given VCC 615's undisturbed morphology, we argue that the galaxy has experienced no recent and sudden transformation into a UDG due to the cluster potential, but rather is a long-lived UDG whose relatively wide orbit and large dynamical mass protect it from stripping and destruction by Virgo cluster tides. Finally, we also describe the serendipitous discovery of a nearby Virgo dwarf galaxy projected 90 arcseconds (7.2 kpc) away from VCC 615.

The hierarchy of motions that we are participating is well known, from the Earth's motion around the Sun and Sun's motion in the Milky Way, up to the Local Group's motion within the Virgo Supercluster of galaxies. The dipole anisotropy of the Cosmic Microwave Background (CMB) enables to define Sun's motion with respect to the CMB "absolute" frame. We now present evidence for CMB Galactic dipole signal due to the Sun's motion in the Galaxy, by means of Planck 2018 data and Gaia Early Data Release 3. The signal is weak and frequency depended, the strongest is at 30 GHz, up to 7.6\sigma confidence. The amplitude of the signal interpreted as Doppler caused, corresponds to velocity $v \approx 225.5 \pm 16.2 \, km\, sec^{-1}$, in agreement with Solar system's velocity with respect to the Galactic center. While the revealing of precise coordinates of the apex will need further refined analysis at various bands, the detected weak signal can indicate the appearance of a new cosmic scaling in CMB, thus opening a link to a bunch of physical effects.

Very long baseline interferometric (VLBI) localisations of repeating fast radio bursts (FRBs) have demonstrated a diversity of local environments: from nearby star-forming regions to globular clusters. Here we report the VLBI localisation of FRB 20201124A using an ad-hoc array of dishes that also participate in the European VLBI Network (EVN). In our campaign, we detected 18 total bursts from FRB 20201124A at two separate epochs. By combining the visibilities from both observing epochs, we were able to localise FRB 20201124A with a 1-$\sigma$ error of 4.5 milliarcseconds (mas). We use the relatively large burst sample to investigate astrometric accuracy, and find that for $\gtrsim20$ baselines ($\gtrsim7$ dishes) that we can robustly reach milliarcsecond precision even using single-burst data sets. Sub-arcsecond precision is still possible for single bursts, even when only $\sim$ six baselines (four dishes) are available. We explore two methods for determining the individual burst positions: the peaks of the dirty maps and a Gaussian fit to the cross fringe pattern on the dirty maps. We found the latter to be more reliable due to the lower mean and standard deviation in the offsets from the FRB position. Our VLBI work places FRB 20201124A 705$\pm$26 mas (1-$\sigma$ errors) from the optical centre of the host galaxy, and consistent with originating from within the recently-discovered extended radio structure associated with star-formation in the host galaxy. Future high-resolution optical observations, e.g. with Hubble Space Telescope, can determine the proximity of our FRB 20201124A VLBI position to nearby knots of star formation.

In this report we summarize the activities of the International Astronomical Consortium for High Energy Calibration (IACHEC) and the work done since the last in-person meeting in Japan (Shonan Village Center), May 2019, through two virtual meetings that were held in November 2020 and May 2021. The on-line only meetings divided the contents of the usual in-person workshop between mission updates and working group updates. The November meeting was dedicated to mission calibration updates and the current status of the cross-calibration between NuSTAR, Swift, and NICER, which frequently join together in observations of bright transients, and a review of the XMM-Newton and Chandra cross-calibration. Results between \nustar\ and \swift\ overall show good agreement, but issues persist in the overlap region 3--5 keV for bright source with large dust scattering halos. The NICER cross-calibration is still progressing and evolving, while for the XMM-Newton and Chandra cross-calibration systematic differences both in the absolute flux and spectral shape determination still exists on different classes of sources. The meeting in May was focused on the Working Group progress and reports summarized here.

We present observations obtained with the VLT/MUSE optical integral field spectrograph of the radio source 3C277.3, located at a redshift of 0.085 and associated with the galaxy Coma A. An emission line region fully enshrouds the double-lobed radio source, which is ~60 kpc x 90 kpc in size. Based on the emission line ratios, we identified five compact knots in which the gas ionization is powered by young stars located as far as ~60 kpc from the host. The emission line filaments surrounding the radio emission are compatible with ionization from fast shocks (with a velocity of 350-500 km/s), but a contribution from star formation occurring at the edges of the radio source is likely. Coma A might be a unique example in the local Universe in which the expanding outflow triggers star formation throughout the whole radio source.

Supermassive black hole masses (MBH) can dynamically be estimated with various methods and using different kinematic tracers. Different methods have only been cross-checked for a small number of galaxies and often show discrepancies. To understand these discrepancies, detailed cross-comparisons of additional galaxies are needed. We present the first part of our cross-comparison between stellar- and gas-based MBH estimates in the nearby fast-rotating early-type galaxy NGC 6958. The measurements presented here are based on ground-layer adaptive optics-assisted Multi-Unit Spectroscopic Explorer (MUSE) science verification data at around 0.6 arcsec spatial resolution. The spatial resolution is a key ingredient for the measurement and we provide a Gaussian parametrisation of the adaptive optics-assisted point spread function (PSF) for various wavelengths. From the MUSE data, we extracted the stellar kinematics and constructed dynamical models. Using an axisymmetric Schwarzschild technique, we measured an MBH of (3.6+2.7-2.4)\times 10^8 Msun at 3\sigma significance taking kinematical and dynamical systematics (e.g.,radially-varying mass-to-light ratio) into account. We also added a dark halo, but our data does not allow to constrain the dark matter fraction. Adding dark matter with an abundance matching prior results in a 25 per cent more massive black hole. Jeans anisotropic models return MBH of (4.6+2.5-2.7) \times 10^8 Msun and (8.6+0.8-0.8) \times 10^8 Msun at 3\sigma confidence for spherical and cylindrical alignment of the velocity ellipsoid, respectively. In a follow-up study, we will compare the stellar-based MBH with those from cold and warm gas tracers, which will provide additional constraints for the MBH for NGC 6958, and insights into assumptions that lead to potential systematic uncertainty.

The most massive halos of matter in the Universe grow via accretion and merger events throughout cosmic times. These violent processes generate shocks at many scales and induce large-scale bulk and turbulent motions. These processes inject kinetic energy at large scales, which is transported to the viscous dissipation scales, contributing to the overall heating and virialisation of the halo, and acting as a source of non-thermal pressure in the intra-cluster medium. Characterizing the physical properties of these gas motions will help us to better understand the assembly of massive halos, hence the formation and the evolution of these large-scale structures. We base this characterization on the study of the X-ray and Sunyaev-Zel'dovich effect brightness fluctuations. Our work relies on three complementary samples covering a wide range of redshifts, masses and dynamical states of clusters. We present the results of our X-ray analysis for the low redshift sample, X-COP, and a subsample of higher redshift clusters. We investigate the derived properties according to the dynamical state of our clusters, and the possibility of a self-similar behaviour based on the reconstructed gas motions power-spectra and the correlation with various morphological indicators.

Small scale CMB angular power spectrum contains not only primordial CMB information but also many contaminants coming from secondary anisotropies. Most of the latter depend on the cosmological model but are often marginalised over. We propose a new analysis of the SPT data focusing on the cosmological contribution of the thermal Sunyaev Zel'dovich (tSZ) effect. We model the tSZ angular spectrum with the halo model and train a random forest algorithm to speed up its computation. We show that using the cosmological information of the tSZ on top of the primordial CMB one contained in SPT data bring more constraints on cosmological parameters. We also combine for the first time Planck tSZ angular power spectrum with SPT ones to put further constraints. This proof of concept study shows how much a proper modelling of the foregrounds in the cosmological analyses is needed.

We present the results of two years of radio and X-ray monitoring of the magnetar XTE J1810$-$197 since the radio re-activation in late 2018. Single pulse analysis of radio observations from the Lovell and MkII telescopes at 1564 MHz and the Effelsberg telescope at 6 GHz has resulted in the detection of a total of 91 giant pulses (GPs) between MJDs 58858 and 59117. These GPs appear to be confined to two specific phase ranges (0.473 <= \phi <= 0.502$ and 0.541 <= \phi <= 0.567). We also observe that the first detection of GP emission corresponds to a minimum in the spin-down rate. Simultaneous radio and X-ray observations were performed on MJDs 59009 and 59096. The 0.5-10 keV X-ray spectrum from NICER is well characterised by a two component blackbody model that can be interpreted as two hot spots on the polar cap of the neutron star. The blackbody temperature decreases with time, consistent with the previous outburst, while the change in the pulsed fraction does not follow the same trend as was seen in the previous outburst. The radio and X-ray flux of XTE J1810-197 are correlated during the initial phase of the outburst (MJD 58450 - MJD 58550) and an increase in the radio flux is observed later that may be correlated to the onset of GPs. We argue that the disparity in the evolution of the current outburst compared to the previous one can be attributed to a change in geometry of the neutron star.

I use recent observations of circumstellar matter around type Ia supernovae (SNe Ia) to estimate the fraction of SNe Ia that explode into a planetary nebula (PN) and to suggest a new delay time distribution from the common envelope evolution (CEE) to the SN Ia explosion for SNe Ia that occur shortly after the CEE. I crudely estimate that about 50 per cent of all SNe Ia are SNe Ia inside PNe (SNIPs), and that the explosions of most SNIPs occur within a CEE to explosion delay (CEED) time of less than about ten thousand years. I also estimate that the explosion rate of SNIPs, i.e., the CEED time distribution (CEEDTD), is roughly constant within this time scale of ten thousand years. I argue that the short CEED time suggests than the majority, and even all, of SNIPs come from the core-degenerate (CD) scenario where the merger of the core with the white dwarf takes place at the end of the CEE. I list some further observations that might support or reject my claims, and the challenge to theoretical studies to find a process to explain a merger to explosion delay (MED) time of up to ten thousand years or so. A long MED will apply also to the double degenerate scenario.

One of the key elements needed to perform the cosmological exploitation of a cluster survey is the relation between the survey observable and the cluster masses. Among these observables, the integrated Compton parameter $Y$ is a measurable quantity in Sunyaev-Zeldovich (SZ) surveys, which tightly correlates with cluster mass. The calibration of the relation between the Compton parameter $Y_{500}$ and the mass $M_{500}$ enclosed within radius $R_{500}$ is one of the scientific goals of the NIKA2 SZ Large Program (LPSZ). We present an ongoing study to forecast the constraining power of this program, using mock simulated datasets that mimic the large program sample, selection function, and typical uncertainties on $Y_{500}$ and $M_{500}$. We use a Bayesian hierarchical modeling that enables taking into account a large panel of systematic effects. Our results show that the LPSZ can yield unbiased estimates of the scaling relation parameters for realistic input parameter values. The relative uncertainties on these parameters is $\sim 10\%$ for the intercept and slope of the scaling relation, and $\sim 34\%$ for its intrinsic scatter, foreshadowing precise estimates to be delivered by the LPSZ.

High-mass stars ($m_* \gtrsim 8 \, M_\odot$) play a crucial role in the evolution of galaxies, and so it is imperative that we understand how they are formed. We have used the New IRAM KIDs Array 2 (NIKA2) camera on the Institut de Radio Astronomie Millim\'{e}trique (IRAM) 30-m telescope to conduct high-sensitivity continuum mapping of $\sim2$ deg$^2$ of the Galactic plane (GP) as part of the Galactic Star Formation with NIKA2 (GASTON) large program. We have identified a total of 1467 clumps within our deep 1.15 mm continuum maps and, by using overlapping continuum, molecular line, and maser parallax data, we have determined their distances and physical properties. By placing them upon an approximate evolutionary sequence based upon 8 $\mu$m $\textit{Spitzer}$ imaging, we find evidence that the most massive dense clumps accrete material from their surrounding environment during their early evolution, before dispersing as star formation advances, supporting clump-fed models of high-mass star formation.

The Milky Way was shaped by the mergers with several galaxies in the past. We search for remnant stars that were born in these foreign galaxies and assess their ages in an effort to put upper limits on the merger times and thereby better understand the evolutionary history of our Galaxy. Using 6D-phase space information from Gaia eDR3 and chemical information from APOGEE DR16, we kinematically and chemically select $23$ red giant stars belonging to former dwarf galaxies that merged with the Milky Way. With added asteroseismology from Kepler and K2, we determine the ages of the $23$ ex-situ stars and $55$ in-situ stars with great precision. We find that all the ex-situ stars are consistent with being older than $8$ Gyr. While it is not possible to associate all the stars with a specific dwarf galaxy we classify eight of them as Gaia-Enceladus/Sausage stars, which is one of the most massive mergers in our Galaxy's history. We determine their mean age to be $9.5^{+1.2}_{-1.3}$ Gyr consistent with a merger time of $8$-$10$ Gyr ago. The rest of the stars are possibly associated with Kraken, Thamnos, Sequoia, or another extragalactic progenitor. The age determination of ex-situ stars paves the way to more accurately pinning down when the merger events occurred and hence provide tight constraints useful for simulating how these events unfolded.

Gravitationally lensed supernovae (glSNe) are of interest for time delay cosmology and SN physics. However, glSNe detections are rare, owing to the intrinsic rarity of SN explosions, the necessity of alignment with a foreground lens, and the relatively short window of detectability. We present the Las Cumbres Observatory Lensed Supernova Search, LCOLSS, a targeted survey designed for detecting glSNe in known strong-lensing systems. Using cadenced $r^\prime$-band imaging, LCOLSS targeted 112 galaxy-galaxy lensing systems with high expected SN rates, based on estimated star formation rates. No plausible glSN was detected by LCOLSS over two years of observing. The analysis performed here measures a detection efficiency for these observations and runs a Monte Carlo simulation using the predicted supernova rates to determine the expected number of glSN detections. The results of the simulation suggest an expected number of detections and $68\%$ Poisson confidence intervals, $N_{SN} = 0.20, [0,2.1] $, $N_{Ia} = 0.08, [0,2.0]$, $N_{CC} = 0.12, [0,2.0]$, for all SN, Type Ia, and core-collapse (CC) SNe respectively. These results are broadly consistent with the absence of a detection in our survey. Analysis of the survey strategy can provide insights for future efforts to develop targeted glSN discovery programs. We thereby forecast expected detection rates for the Rubin observatory for such a targeted survey, finding that a single visit depth of 24.7 mag with the Rubin observatory will detect $0.63 \pm 0.38$ SNe per year, with $0.47 \pm 0.28$ core collapse SNe per year and $0.16 \pm 0.10$ Type Ia SNe per year.

We present a multi-probe analysis of the well-known galaxy cluster CL J1226.9+3332 as a proof of concept for multi-wavelength studies within the framework of the NIKA2 Sunyaev-Zeldovich Large Program (LPSZ). CL J1226.9+3332 is a massive and high redshift (z = 0.888) cluster that has already been observed at several wavelengths. A joint analysis of the thermal SZ (tSZ) effect at millimeter wavelength with the NIKA2 camera and in X-ray with the XMM-Newton satellite permits the reconstruction of the cluster thermodynamical properties and mass assuming hydrostatic equilibrium. We test the robustness of our mass estimates against different definitions of the data analysis transfer function. Using convergence maps reconstructed from the data of the CLASH program we obtain estimates of the lensing mass, which we compare to the estimated hydrostatic mass. This allows us to measure the hydrostatic-to-lensing mass bias and the associated systematic effects related to the NIKA2 measurement. We obtain M500HSE = (7.65 +- 1.03) 1014 Msun and M500lens = (7.35 +- 0.65) 1014 Msun, which implies a HSE-to-lensing bias consistent with 0 within 20 percent.

Starting from the clusters included in the NIKA sample and in the NIKA2 Sunyaev-Zel'dovich Large Program (LPSZ) we have selected a sample of six common objects with the Cluster Lensing And Supernova survey with Hubble (CLASH) lensing data. For the LPSZ clusters we have at our disposal both high-angular resolution observations of the thermal SZ with NIKA and NIKA2 and X-ray observations with XMM-Newton from which hydrostatic mass estimates can be derived. In addition, the CLASH dataset includes lensing convergence maps that can be converted into lensing estimates of the total mass of the cluster. One-dimensional mass profiles are used to derive integrated mass estimates accounting for systematic effects (data processing, modeling, etc.). Two-dimensional analysis of the maps can reveal substructures in the cluster and, therefore, inform us about the dynamical state of each system. Moreover, we are able to study the hydrostatic mass to lensing mass bias, across different morphology and a range of redshift clusters to give more insight on the hydrostatic mass bias. The analysis presented in this proceeding follows the study discussed in Ferragamo et al. 2021.

The irrefutable successes of MOND are predicated upon the idea that a critical gravitational acceleration scale, $a_0$, exists. But, beyond its role in MOND, the question: 'Why should a critical gravitational acceleration scale exist at all?' remains unanswered. There is no deep understanding about what is going on. Over roughly the same period that MOND has been a topic of controversy, Baryshev, Sylos Labini, Pietronero and others have been arguing, with equal controversy in earlier years, that, on medium scales at least, material in the universe is distributed in a quasi-fractal $D\approx 2$ fashion. There is a link: if the idea of a quasi-fractal $D \approx 2$ universe on medium scales is taken seriously then there is an associated characteristic mass surface density scale, $\Sigma_F$ say, and an associated characteristic gravitational acceleration scale, $a_F = 4 \pi G \Sigma_F$. If, furthermore, the quasi-fractal structure is taken to include the inter-galactic medium, then it is an obvious step to consider the possibility that $a_0$ and $a_F$ are the same thing. Through the lens of very old ideas rooted in a Leibniz-Mach worldview we obtain a detailed understanding of the critical acceleration scale which, applied to the SPARC sample of galaxies with a stellar MLR, $\Upsilon_* \in (0.5, 1.0)$, and using standard photometric mass models, provides a finite algorithm to recover the information that $a_F \approx 1.2 \times 10^{-10}\, mtrs/sec^2$. This, combined with the fact that the Baryonic Tully-Fisher Relationship (BTFR) arises directly from the same source, but with $a_0$ replaced by $a_F$, leads to the unambiguous conclusion that $a_0$ and $a_F$ are, in fact, one and the same thing.

Context. Protoplanetary disks around young stars often contain substructures like rings, gaps, and spirals that could be caused by interactions between the disk and forming planets. Aims. We aim to study the young (1-3 Myr) star DR Tau in the near-infrared and characterize its disk, which was previously resolved through sub-millimeter interferometry with ALMA, and to search for possible sub-stellar companions embedded into it. Methods. We observed DR Tau with VLT/SPHERE both in polarized light (H broad band) and total intensity (in Y, J, H, and K spectral bands). We also performed L' band observations with LBTI/LMIRCam on the Large Binocular Telescope (LBT). Results. We found two previously undetected spirals extending north-east and south of the star, respectively. We further detected an arc-like structure north of the star. Finally a bright, compact and elongated structure was detected at separation of 303 +/- 10 mas and position angle 21.2 +/- 3.7 degrees, just at the root of the north-east spiral arm. Since this feature is visible both in polarized light and in total intensity and has a flat spectrum it is likely caused by stellar light scattered by dust. Conclusions. The two spiral arms are at different separation from the star, have very different pitch angles, and are separated by an apparent discontinuity, suggesting they might have a different origin. The very open southern spiral arm might be caused by infalling material from late encounters with cloudlets into the formation environment of the star itself. The compact feature could be caused by interaction with a planet in formation still embedded in its dust envelope and it could be responsible for launching the north-east spiral. We estimate a mass of the putative embedded object of the order of few M_Jup .

We show that a Denoising Diffusion Probabalistic Model (DDPM), a class of score-based generative model, can be used to produce realistic yet fake images that mimic observations of galaxies. Our method is tested with Dark Energy Spectroscopic Instrument grz imaging of galaxies from the Photometry and Rotation curve OBservations from Extragalactic Surveys (PROBES) sample and galaxies selected from the Sloan Digital Sky Survey. Subjectively, the generated galaxies are highly realistic when compared with samples from the real dataset. We quantify the similarity by borrowing from the deep generative learning literature, using the `Fr\'echet Inception Distance' to test for subjective and morphological similarity. We also introduce the `Synthetic Galaxy Distance' metric to compare the emergent physical properties (such as total magnitude, colour and half light radius) of a ground truth parent and synthesised child dataset. We argue that the DDPM approach produces sharper and more realistic images than other generative methods such as Adversarial Networks (with the downside of more costly inference), and could be used to produce large samples of synthetic observations tailored to a specific imaging survey. We demonstrate two potential uses of the DDPM: (1) accurate in-painting of occluded data, such as satellite trails, and (2) domain transfer, where new input images can be processed to mimic the properties of the DDPM training set. Here we `DESI-fy' cartoon images as a proof of concept for domain transfer. Finally, we suggest potential applications for score-based approaches that could motivate further research on this topic within the astronomical community.

In the current proceedings, we summarise the results presented during the mm Universe@NIKA2 conference, taken from our main results in \cite{Haridasu:2021hzq}. We test the Degenerate higher-order scalar-tensor(DHOST) theory as a generalised platform for scalar-tensor theory at galaxy cluster scales to predict in such static systems small scale modification to the gravitational potential. {DHOST theory is not only a good alternative to $\Lambda$CDM for the background evolution but also predicts small-scale modification to the gravitational potential in static systems such as galaxy clusters.} With a sample of 12 clusters with accurate X-ray Intra Cluster Medium (ICM) data (X-COP project) and Sunyaev-Zel'dovich (SZ) ICM pressure (Planck satellite), we place preliminary constraints on the DHOST parameters defining the deviation from GR. Moreover, we also collect a few supplementary analyses we have performed during the course: i) Gaussian process reconstruction without parametric assumptions, ii) $\Psz$-only data analysis not aided by the X-ray data. Finally, we present possible extensions to the current work which may benefit from future high sensitivity and spatial resolution observations such as the ongoing NIKA2 camera.

The Planck Collaboration has shown that the number of clusters as a function of their mass and redshift is an extremely powerful tool for cosmological analyses. However, the true cluster mass is not directly measurable. Among the possible approaches, clusters mass could be related to different observables via self similar scaling law. These observables are related to the baryonic components of which a cluster is composed. However, the theoretical relations that allow the use of these proxies often are affected by observational and physical biases, which impacts on the determination of the cluster mass. Fortunately, cosmological simulations are an extremely powerful tool to assess these problems. We present our calibration of the scaling relation between mass and velocity dispersion of galaxy members from the study of the simulated clusters of \THP{} project with mass above $10^{13} M_\odot$. In order to investigate the presence of a redshift dependence, we analyzed 16 different redshifts between $z = 0$ and $z = 2$. Finally, we investigated the impact of different AGN feedback models.

The NIKA2 Guaranteed-Time SZ Large Program (LPSZ) is dedicated to the high-angular resolution SZ mapping of a representative sample of 45 SZ-selected galaxy clusters drawn from the catalogues of the Planck satellite, or of the Atacama Cosmology Telescope. The LPSZ sample spans a mass range from $3$ to $11 \times 10^{14} M_{\odot}$ and a redshift range from $0.5$ to $0.9$, extending to higher redshift and lower mass the previous samples dedicated to the cluster mass calibration and universal properties estimation. The main goals of the LPSZ are the measurement of the average radial profile of the ICM pressure up to $R_{500}$ by combining NIKA2 with Planck or ACT data, and the estimation of the scaling law between the SZ observable and the mass using NIKA2, XMM-Newton and Planck/ACT data. Furthermore, combining LPSZ data with existing or forthcoming public data in lensing, optical/NIR or radio domains, we will build a consistent picture of the cluster physics and further gain knowledge on the mass estimate as a function of the cluster morphology and dynamical state.

Cosmic dust grains are one of the fundamental ingredients of the interstellar medium (ISM). In spite of their small contribution to the total mass budget, dust grains play a significant role in the physical and chemical evolution of galaxies. Over the past decades, a plethora of multi-wavelength data, from UV to far-infrared, has increased substantially our knowledge on the dust properties of nearby galaxies. Nevertheless, one regime of the spectrum, the mm range, remains relatively unexplored. Thanks to the new, high-resolution data in the mm range observed with the NIKA2 instrument and our radiative transfer framework, we aim to firmly characterise the physical properties of the very cold dust (<15K), and to quantify the importance of different emission mechanisms in the mm. So far, we have developed a methodology to use dust radiative transfer modelling and applied it to a small group of face-on spiral galaxies. The combination of the new NIKA2 data with our radiative transfer techniques would provide the right conditions to generate an accurate model of the interplay between starlight and dust in a sizeable sample of spatially-resolved nearby galaxies.

New observations of the edge-on galaxy NGC 891, at 1.15 and 2 mm obtained with the IRAM 30-m telescope and the NIKA2 camera, within the framework of the IMEGIN (Interpreting the Millimetre Emission of Galaxies with IRAM and NIKA2) Large Program, are presented in this work. By using multiwavelength maps (from the mid-IR to the cm wavelengths) we perform SED fitting in order to extract the physical properties of the galaxy on both global and local ($\sim$kpc) scales. For the interpretation of the observations we make use of a state-of-the-art SED fitting code, HerBIE (HiERarchical Bayesian Inference for dust Emission). The observations indicate a galaxy morphology, at mm wavelengths, similar to that of the cold dust emission traced by sub-mm observations and to that of the molecular gas. The contribution of the radio emission at the NIKA2 bands is very small (negligible at 1.15 mm and $\sim10\%$ at 2 mm) while it dominates the total energy budget at longer wavelengths (beyond 5 mm). On local scales, the distribution of the free-free emission resembles that of the dust thermal emission while the distribution of the synchrotron emission shows a deficiency along the major axis of the disc of the galaxy.

In this article we treat the non-adiabatic resonant photon-to-axion conversion of curvature radiation, synchrotron emission and inverse Compton scattering dominating the spectral density function of pulsars. First, we introduce state-of-the-art emission models and reference observational data. Relying on recurring idealizations in the literature, we derive an analytical expression that allows one to estimate the photon-axion oscillation probability across the light cylinder of the neutron star. Then, we estimate the dark matter flux induced by photon-to-axion conversion. We find that pulsars might produce axion overdensities many orders of magnitude over the occupation number of dark matter in the Galactic halo within a broad parameter space. We discuss axion-astronomy prospects and heterogeneous future lines.

Several galaxy clusters host X-ray cavities, often filled with relativistic electrons emitting in the radio band. In the cluster MS 0735.6+7421 the cavities have been detected through the Sunyaev Zel'dovich (SZ) effect, but it has not been possible to determine if this effect is thermal (produced by a very high temperature gas filling the cavity) or non-thermal (produced by the relativistic electrons that produce the diffuse radio emission detected in the cavity). In this paper we discuss the role of the density of the high temperature gas inside the cavities in determining whether the dominant SZ effect is the thermal or the non-thermal one, and how it can be possible to distinguish between the two possibilities, discussing the role of observations at higher energy bands.

Thanks to the relative ease of finding and characterizing small planets around M dwarf stars, these objects have become cornerstones in the field of exoplanet studies. The current paucity of planets in long-period orbits around M dwarfs make such objects particularly compelling as they provide clues about the formation and evolution of these systems. In this study, we present the discovery of TOI-2257 b (TIC 198485881), a long-period (35 d) sub-Neptune orbiting an M3 star at 57.8pc. Its transit depth is about 0.4%, large enough to be detected with medium-size, ground-based telescopes. The long transit duration suggests the planet is in a highly eccentric orbit ($e \sim 0.5$), which would make it the most eccentric planet that is known to be transiting an M-dwarf star. We combined TESS and ground-based data obtained with the 1.0-m SAINT-EX, 0.60-m TRAPPIST-North and 1.2-m FLWO telescopes to find a planetary size of 2.2 $R_{\oplus}$ and an orbital period of 35.19 days. In addition, we make use of archival data, high-resolution imaging, and vetting packages to support our planetary interpretation. With its long period and high eccentricity, TOI-2257 b falls in a novel slice of parameter space. Despite the planet's low equilibrium temperature ($\sim$ 256 K), its host star's small size ($R_* = 0.311 \pm{0.015}$) and relative infrared brightness (K$_{mag}$ = 10.7) make it a suitable candidate for atmospheric exploration via transmission spectroscopy.

The image of a supermassive black hole surrounded by an optically-thin, radiatively-inefficient accretion flow, like that observed with the Event Horizon Telescope, is characterized by a bright ring of emission surrounding the black-hole shadow. In the Kerr spacetime this bright ring, when narrow, closely traces the boundary of the shadow and can, with appropriate calibration, serve as its proxy. The present paper expands the validity of this statement by considering two particular spacetime geometries: a solution to the field equations of a modified gravity theory and another that parametrically deviates from Kerr but recovers the Kerr spacetime when its deviation parameters vanish. A covariant, axisymmetric analytic model of the accretion flow based on conservation laws and spanning a broad range of plasma conditions is utilized to calculate synthetic non-Kerr black-hole images, which are then analysed and characterized. We find that in all spacetimes: (i) it is the gravitationally-lensed unstable photon orbit that plays the critical role in establishing the diameter of the rings observed in black-hole images, not the event horizon or the innermost stable circular orbit, (ii) bright rings in these images scale in size with, and encompass, the boundaries of the black-hole shadows, even when deviating significantly from Kerr, and (iii) uncertainties in the physical properties of the accreting plasma introduce subdominant corrections to the relation between the diameter of the image and the diameter of the black-hole shadow. These results provide theoretical justification for using black-hole images to probe and test the spacetimes of supermassive black holes.

The detection of the gravitational wave counterpart GRB 170817A, underluminous compared to the cosmological GRB population by a factor of 10,000, motivates significant effort in detecting and localizing a dim, nearby, and slightly off-axis population of short GRBs. Swift/BAT is the most sensitive GRB detector in operation, and the only one that regularly localizes GRBs to arcminute precision, critical to rapid followup studies. However, the utility of BAT in targeted sub-threshold searches had been historically curtailed by the unavailability of the necessary raw data for analysis. The new availability of time-tagged event (TTE) data from the GUANO system (arXiv:2005.01751), motivates renewed focus on developing sensitive targeted search analysis techniques to maximally exploit these data. While computationally cheap, we show that the typical coded-mask deconvolution imaging is limited in its sensitivity due to several factors. We formalize a maximum likelihood framework for the analysis of BAT data wherein signals are forward modelled through the full instrument response, and -- coupled with the development of new response models -- demonstrate its superior sensitivity to typical imaging via archival comparisons, injection campaigns, and, after implementing as a targeted search, a large number of low-latency GRB discoveries and confirmed arcminute localizations to date. We also demonstrate independent localization of some out-of-FOV GRBs for the first time. NITRATES's increased sensitivity boosts the discovery rate of GRB 170817A-like events in BAT by a factor of at least $3-4$x, along with enabling joint analyses and searches with other GRB, GW, neutrino, and FRB instruments. We provide public access to the response functions and search pipeline code.

The results of a temporal analysis of observations for a sample of nine low-mass X-ray binaries (LMXBs) are presented. Of these sources, five host a neutron star (NS) primary (4U1608-52, Aql X-1, 4U1705-44, GX17+2, and Cyg X-2), and four host a black hole (BH) (GX339-4, XTE J1859+226, H1743-322, and MAXI J1659-152). The NS group includes three Atolls and two Z-type sources. We utilized archival Proportional Counter Array (PCA)/RXTE data to construct high-resolution lightcurves. A wavelet transform of the lightcurves is deployed to extract a minimal time-scale (MTS) associated with the spectral state of the sources. The MTS, together with the the fractional root-mean-square (RMS) and hardness ratios, is used to construct RMS-MTS and hardness-MTS diagrams that enable a direct comparison of the evolution of spectral transitions in the target sources. Observations with high fractional RMS and high hardness cluster in a broad region occupied jointly by BH and NS sources. For low fractional RMS observations, the MTS (on average) separates Atolls from Z-type and BH sources: small MTS corresponds to Z-types (and BHs), and large MTS to Atoll sources. Furthermore, in the hardness-MTS plane, BH sources are the sole occupiers of the low hardness, small MTS region thus potentially signaling a unique property for distinguishing BH and NS hosts in LMXBs.

Due to the different environments in the Milky Way's disk and halo, comparing wide binaries in the disk and halo is key to understanding wide binary formation and evolution. By using Gaia Early Data Release 3, we search for resolved wide binary companions in the H3 survey, a spectroscopic survey that has compiled ~150,000 spectra for thick-disk and halo stars to date. We identify 800 high-confidence (a contamination rate of 4%) wide binaries and two resolved triples, with binary separations mostly between $10^3$-$10^5$ AU and a lowest [Fe/H] of -2.7. Based on their Galactic kinematics, 33 of them are halo wide binaries, and most of those are associated with the accreted Gaia-Sausage-Enceladus galaxy. The wide binary fraction in the disk decreases toward the low metallicity end, consistent with Hwang et al. (2021). Our key finding is that the halo wide binary fraction is consistent with the thick-disk stars at a fixed [Fe/H]. There is no significant dependence of the wide binary fraction on the $\alpha$-capture abundance. Therefore, the wide binary fraction is mainly determined by the iron abundance, not their disk or halo origin nor the $\alpha$-captured abundance. Our results suggest that the formation environments play a major role for the wide binary fraction, instead of other processes like radial migration that only apply to disk stars.

Eccentricity of wide binaries is difficult to measure due to their long orbital periods. With Gaia's high-precision astrometric measurements, eccentricity of a wide binary can be constrained by the angle between the separation vector and the relative velocity vector (the $v$-$r$ angle). In this paper, by using the $v$-$r$ angles of wide binaries in Gaia Early Data Release 3, we develop a Bayesian approach to measure the eccentricity distribution as a function of binary separations. Furthermore, we infer the eccentricities of individual wide binaries and make them publicly available. Our results show that the eccentricity distribution of wide binaries at $10^2$ AU is close to uniform and becomes superthermal at $>10^{3}$ AU, suggesting two formation mechanisms dominating at different separation regimes. The close binary formation, most likely disk fragmentation, results in a uniform eccentricity distribution at $<10^{2}$ AU. The wide binary formation that leads to highly eccentric wide binaries at $>10^{3}$ AU may be turbulent fragmentation and/or the dynamical unfolding of compact triples. With Gaia, measuring eccentricities is now possible for a large number of wide binaries, opening a new window to understanding binary formation and evolution.

We explore superradiance for Alfv\'en waves (Alfv\'enic superradiance) in an axisymmetric rotating magnetosphere of a Kerr black hole within the force-free approximation. On the equatorial plane of the Kerr spacetime, the Alfv\'en wave equation is reduced to a one-dimensional Schr\"{o}dinger-type equation by separating variables of the wave function and introducing a tortoise coordinate mapping the inner and outer light surfaces to $-\infty$ and $+\infty$, respectively, and we investigate a wave scattering problem for Alfv\'en waves. An analysis of the asymptotic solutions of the wave equation and the conservation of the Wronskian give the superradiant condition for Alfv\'en waves, and it will be shown that the condition coincides with that for the Blandford-Znajek process. This means that when Alfv\'enic superradiance occurs, the Blandford-Znajek process also occurs in the force-free magnetosphere. Then, we evaluate the reflection rate of Alfv\'en waves numerically and confirm that Alfv\'enic superradiance is indeed possible in the Kerr spacetime. Moreover, we will discuss a resonant scattering of Alfv\'en waves, which is related to a "quasinormal mode" of the magnetosphere.

Kination denotes an era in the cosmological history corresponding to an equation of state $\omega=+1$ such that the total energy density of the universe redshifts as the sixth inverse power of the scale factor. This arises if the universe is dominated by the kinetic energy of a scalar field. It has often been motivated in the literature as an era following inflation, taking place before the radiation era. In this paper, we review instead the possibility that kination is disconnected from primordial inflation and occurs much later, inside the Standard Model radiation era. We study the implications on all main sources of primordial gravitational waves. We show how this leads to very distinctive peaked spectra in the stochastic background of long-lasting cosmological sources of gravitational waves, namely the irreducible gravitational waves from inflation, and gravitational waves from cosmic strings, both local and global, with promising observational prospects. We present model-independent signatures and detectability predictions at LIGO, LISA, ET, CE, BBO, as a function of the energy scale and duration of the kination era. We then argue that such intermediate kination era is in fact symptomatic in a large class of axion models. We analyse in details the scalar field dynamics, the working conditions and constraints in the underlying models. We present the gravitational-wave predictions as a function of particle physics parameters. We derive the general relation between the gravitational-wave signal and the axion dark matter abundance as well as the baryon asymmetry. We investigate the predictions for the special case of the QCD axion. The key message is that gravitational-waves of primordial origin represent an alternative experimental probe of axion models.

We consider imaging of fast moving small objects in space, such as low earth orbit satellites, which are also rotating around a fixed axis. The imaging system consists of ground based, asynchronous sources of radiation and several passive receivers above the dense atmosphere. We use the cross-correlation of the received signals to reduce distortions from ambient medium fluctuations. Imaging with correlations also has the advantage of not requiring any knowledge about the probing pulse and depends weakly on the emitter positions. We account for the target's orbital velocity by introducing the necessary Doppler compensation. To image a fast rotating object we also need to compensate for the rotation. We show that the rotation parameters can be extracted directly from the auto-correlation of the data before the formation of the image. We then investigate and analyze an imaging method that relies on backpropagating the cross-correlation data structure to two points rather than one, thus forming an interference matrix. The proposed imaging method consists of estimating the reflectivity as the top eigenvector of the migrated cross-correlation data interference matrix. We call this the rank-1 image and show that it provides superior image resolution compared to the usual single-point migration scheme for fast moving and rotating objects. Moreover, we observe a significant improvement in resolution due to the rotation leading to a diffraction limited resolution. We carry out a theoretical analysis that illustrates the role of the two point migration method as well as that of the inverse aperture and rotation in improving resolution. Extensive numerical simulations support the theoretical results.

We show how light scalar fields can account for the discrepancy between the theoretical and observed values of the anomalous magnetic moment of the (anti)muon. When coupled to both matter and photons, light scalar fields induce a change of the anomalous magnetic moment of charged particles. This arises from two concurrent effects. Classically, light scalars induce a change of the cyclotron frequency, complementing the electromagnetic effects coming from the magnetic and electric fields used experimentally. Light scalars also contribute to the anomalous magnetic moment quantum mechanically at the one-loop level. For unscreened scalar fields coupling with a Yukawa interaction to matter, these contributions are negligible after applying the Cassini bound on deviations from Newtonian gravity. On the other hand, screened scalars such as chameleons or symmetrons can couple strongly to matter in the laboratory and decouple in the Solar System. This allows us to probe branches of their parameter spaces where the recently measured anomalous magnetic moment of the (anti)muon can be accounted for in the chameleon and symmetron cases. This might be a hint that modified gravity is at play in the laboratory.

We revisit dark-matter production through freeze-in and freeze-out by solving the Boltzmann equations at the level of the phase-space distribution $f(p,t)$. Using the $2\to2$ annihilation and the $1\to2$ decay processes for illustration, we compare the resulting dark-matter relic abundance with that from the number-density approach. In the transition regime between freeze-in and freeze-out, we find the difference can be as large as $\sim$10% for $2\to2$ and $\sim$50% for $1\to2$, or even a factor of $\mathcal{O}(10)$ if the annihilation of dark-matter particles or the decaying mediator is neglected. The freeze-in production in the $2\to2$ and the $1\to2$ processes can also result in non-thermal phase-space distributions, or even multi-modal ones with out-of-equilibrium decay, which can potentially affect structure formation at late times. We also investigate how elastic scatterings can distort such non-thermal distributions.

We apply the formalism of dynamical system analysis to investigate the evolution of interacting dark energy scenarios at the background and perturbation levels in a unified way. Since the resulting dynamical system contains the extra perturbation variable related to the matter overdensity, the critical points of the background analysis split, corresponding to different behavior of matter perturbations, and hence to stability properties. From the combined analysis, we find critical points that describe the non-accelerating matter-dominated epoch with the correct growth of matter structure, and the fact that they are saddle provides the natural exist from this phase. Furthermore, we find stable attractors at late times corresponding to a dark energy-dominated accelerated solution with constant matter perturbations, as required by observations. Thus, interacting cosmology can describe the matter and dark energy epochs correctly, both at the background and perturbation levels, and since this is not possible in standard, i.e. non-interacting, quintessence, it reveals the crucial effect of the interaction.

By combining the noncanonical kinetic term with the nonminimal coupling between gravity and Higgs field, we propose a novel inflationary mechanism and show that the curvature perturbation at small scales can be enhanced by the fluctuations of the Higgs field while satisfying the cosmic microwave background observations at large scales. We find that the recent NANOGrav signal and the black holes in LIGO-Virgo events may originate from the vacuum fluctuations of Higgs field. The cosmological probe of Higgs field by space-based gravitational wave observatory in the future is proposed. We also consider the non-Gaussianity effect on the abundance of primordial black holes and scalar induced secondary gravitational waves. We find that primordial non-Gaussianity makes primordial black holes form more easily, but its effect on the energy density of scalar induced gravitational waves is negligible.

The Deep Space Network is NASA's international array of antennas that support interplanetary spacecraft missions. A track is a block of multi-dimensional time series from the beginning to end of DSN communication with the target spacecraft, containing thousands of monitor data items lasting several hours at a frequency of 0.2-1Hz. Monitor data on each track reports on the performance of specific spacecraft operations and the DSN itself. DSN is receiving signals from 32 spacecraft across the solar system. DSN has pressure to reduce costs while maintaining the quality of support for DSN mission users. DSN Link Control Operators need to simultaneously monitor multiple tracks and identify anomalies in real time. DSN has seen that as the number of missions increases, the data that needs to be processed increases over time. In this project, we look at the last 8 years of data for analysis. Any anomaly in the track indicates a problem with either the spacecraft, DSN equipment, or weather conditions. DSN operators typically write Discrepancy Reports for further analysis. It is recognized that it would be quite helpful to identify 10 similar historical tracks out of the huge database to quickly find and match anomalies. This tool has three functions: (1) identification of the top 10 similar historical tracks, (2) detection of anomalies compared to the reference normal track, and (3) comparison of statistical differences between two given tracks. The requirements for these features were confirmed by survey responses from 21 DSN operators and engineers. The preliminary machine learning model has shown promising performance (AUC=0.92). We plan to increase the number of data sets and perform additional testing to improve performance further before its planned integration into the track visualizer interface to assist DSN field operators and engineers.

In this study, we revisit the absorption and scattering process by which a massless bosonic field impinges on a charged dilatonic black hole. Using the partial wave method, we determine numerically the total absorption cross-section in terms of the decoupling parameter called $M\omega$, finding that the amplitude of the dilatonic black hole is lower than the Reissner-Nordstr\"om one for mild frequencies. In the limit of high-frequency, the absorption cross-section exhibits two different kinds of complex behaviors, one is referred to as the fine structure and the other as the hyperfine structure. To fully grasp the main properties of the charged dilatonic black hole, we consider a different framework where the compact object is impinged by a charged massive bosonic field. In the low-frequency limit, we argue that the absorption cross-section presents two different phases, which are indicated by the value of a critical velocity. Depending on the dark matter model and the black hole mass, we show that both phases are relevant. We verify that the superradiance scattering takes place, being enhanced by smaller scalar field mass but larger values of the scalar field charge. For intermediate frequency, the superradiant effect is lessened in relation to the Reissner-Nordstr\"om case. However, this effect does not necessarily imply the existence of any instability. In order to trigger the superradiant instability, two conditions must be met. There must be unstable modes that remain trapped outside the event horizon with a mechanism based on the reflecting-mirror boundary conditions. This mechanism allows for the system composed of a charged scalar field plus a charged black hole to produce a charged black hole bomb. We provide an analytic formula (lower bound) for the values of the charge field which can trigger this superradiant instability.

Experimental and theoretical studies are presented on the reactivity of methanimine and aminomethylene cations with ethene. Selective isomer generation is performed via dissociative photoionization of suitable neutral precursors and reactive cross sections and branching ratios are measured as a function of photon and collision energies.

We examine the dynamics of FLRW cosmologies in which the vacuum interacts with a perfect fluid through an energy exchange, exploring nonsingular configurations, including cyclic and bouncing models. We consider two specific choices for the energy transfer. In the first case, the energy transfer is proportional to a linear combination of the vacuum and fluid energy densities which makes the conservation equations exactly integrable. The resulting Friedmann equation can be interpreted as an energy constraint equation with an effective potential for the scale factor that may include an infinite barrier forcing a bounce at small values of the scale factor, as well as a potential well allowing for cycling solutions. In the second case, the energy transfer is a nonlinear combination of the vacuum and fluid energy densities. Nonetheless even in this case the dynamics can be partially integrated, leading to a first integral, reducing the number of degrees of freedom. We show that also in this nonlinear case bouncing and cycling cosmologies may arise. In both cases the structure of the resulting phase space allows for nonsingular orbits with an early accelerated phase around a single bounce, connected via a decelerated matter-dominated era to a late-time accelerated phase dominated by an effective cosmological constant.