Papers for Tuesday, Jun 15 2021

We present the largest systematic search to date for luminous $z\gtrsim8$ galaxy candidates using ~1267 arcmin$^{2}$ of (pure-)parallel HST observations from the SuperBoRG data set, a compilation of 316 random sightlines with ACS and WFC3 observations, which together represent a factor ~1.4x larger than existing data sets. Using NIR color cuts and careful photo-$z$ analyses, we find 49 $z\sim8-12$ galaxy candidates over 44 unique sightlines, and derive global galaxy properties such as UV magnitudes and continuum slopes, sizes, and rest-frame optical properties (e.g., SFRs, stellar masses, $A_{\rm v}$). Taking advantage of the (pure-)parallel nature of our data set - making it one of the most representative thus far - and derived SFRs, we evaluate the cosmic star formation rate density for the bright end of the luminosity at $z\sim8-10$ and test the validity of luminosity function-derived results using a conversion factor. We find our method yields comparable results to those derived with luminosity functions. Furthermore, we present follow up observations of 4 (Super)BoRG targets with Keck/MOSFIRE, finding no evidence of Ly$\alpha$ in >3 hrs of $Y-$band observations in either, consistent with a largely neutral medium at $z\sim8$. Our results offer a definitive HST legacy on the bright end of the luminosity function and provide a valuable benchmark as well as targets for follow up with JWST.

Context.Coherence in the characteristics of neighboring sources in 2D and 3D space may suggest the existence of large-scale cosmic structures, which are useful for cosmological studies. Numerous works have been conducted to detect such features in global scalesas well as in confined areas of the sky. However, results are often contradictory and their interpretation remains controversial. Aims.We investigate the potential alignment of parsec-scale radio jets in localized regions of the coordinates-redshift space. Methods.We use data from the Astrogeo VLBI FITS image database to deduce jet directions of radio sources. We perform the search for statistical alignments between nearby sources and explore the impact of instrumental biases. Results.We unveil four regions for which the alignment between jet directions deviates from randomness at a significance level of more than 5 sigma and is unlikely due to instrumental systematics. Intriguingly, their locations coincide with other known large-scale cosmic structures and/or regions of alignments. Conclusions.If the alignments found are the result of physical processes, the discovered regions may designate some of the largest structures known to date.

White dwarf planetary science is a rapidly growing field of research featuring a diverse set of observations and theoretical explorations. Giant planets, minor planets, and debris discs have all been detected orbiting white dwarfs. The innards of broken-up minor planets are measured on an element-by-element basis, providing a unique probe of exoplanetary chemistry. Numerical simulations and analytical investigations trace the violent physical and dynamical history of these systems from au-scale distances to the immediate vicinity of the white dwarf, where minor planets are broken down into dust and gas and are accreted onto the white dwarf photosphere. Current and upcoming ground-based and space-based instruments are likely to further accelerate the pace of discoveries.

We present an analysis of \chandra/LETGS observations of the ultracompact X-ray binary (UCXB) 4U 1626$-$67, continuing our project to analyze the existing Chandra gratings data of this interesting source. The extremely low mass, hydrogen-depleted donor star provides a unique opportunity to study the properties and structure of the metal-rich accreted plasma. There are strong, double-peaked emission features of OVII-VIII and Ne IX-X, but no other identified emission lines are detected. Our spectral fit simultaneously models the emission line profiles and the plasma parameters, using a two-temperature collisionally-ionized plasma. Based on our line profile fitting, we constrain the inclination of the system to 25--60$^{\circ}$ and the inner disk radius to $\sim$1500 gravitational radii, in turn constraining the donor mass to $\lesssim$0.026 M_sun, while our plasma modeling confirms previous reports of high neon abundance in the source, establishing a Ne/O ratio in the system of $0.47 \pm 0.04$, while simultaneously estimating a very low Fe/O ratio of $0.0042 \pm 0.0008$ and limiting the Mg/O ratio to less than 1% by number. We discuss these results in light of previous work.

Paleo-detectors are a proposed experimental technique to search for dark matter (DM). In lieu of the conventional approach of operating a tonne-scale real-time detector to search for DM-induced nuclear recoils, paleo-detectors take advantage of small samples of naturally occurring rocks on Earth that have been deep underground ($\gtrsim 5$ km), accumulating nuclear damage tracks from recoiling nuclei for $\mathcal{O}(1)$ Gyr. Modern microscopy techniques promise the capability to read out nuclear damage tracks with nanometer resolution in macroscopic samples. Thanks to their $\mathcal{O}(1)$ Gyr integration times, paleo-detectors could constitute nuclear recoil detectors with keV recoil energy thresholds and 100 kilotonne-yr exposures. This combination would allow paleo-detectors to probe DM-nucleon cross sections orders of magnitude below existing upper limits from conventional direct detection experiments. In this article, we use improved background modeling and a new spectral analysis technique to update the sensitivity forecast for paleo-detectors. We demonstrate the robustness of the sensitivity forecast to the (lack of) ancillary measurements of the age of the samples and the parameters controlling the backgrounds, systematic mismodeling of the spectral shape of the backgrounds, and the radiopurity of the mineral samples. Specifically, we demonstrate that even if the uranium concentration in paleo-detector samples is $10^{-8}$ (per weight), many orders of magnitude larger than what we expect in the most radiopure samples obtained from ultra basic rock or marine evaporite deposits, paleo-detectors could still probe DM-nucleon cross sections below current limits. For DM masses $\lesssim 10$ GeV/$c^2$, the sensitivity of paleo-detectors could still reach down all the way to the conventional neutrino floor in a Xe-based direct detection experiment.

In the context of scaling relations between Supermassive Black Holes and host-galaxy properties, we aim to enhance the comparison between $M_{\bullet} - M_{G}\sigma^2$ and $M_{\bullet} - \sigma$ relations from a statistical point of view. First, it is suggested to take into account the predictive accuracy of the scaling relation, in addition to the classical measures of goodness of fit. Here, prediction accuracy is fairly evaluated according to a leave-one-out cross-validation strategy. Then, we spread more light on the analysis of residuals from the fitted scaling relation, in order to provide more useful information on the role played by the different variables in their correlation with the black hole mass. The findings from six samples are discussed.

We study the connection between the X-ray and UV properties of the broad absorption line (BAL) wind in the highly X-ray variable quasar PG 2112+059 by comparing Chandra-ACIS data with contemporaneous UV HST/STIS spectra in three different epochs. We observe a correlation whereby an increase in the equivalent-widths (EWs) of the BALs is accompanied by a redder UV spectrum. The growth in the BALs EWs is also accompanied by a significant dimming in soft X-ray emission (<2 keV), consistent with increased absorption. Variations in the hard X-ray emission (>2 keV) are only accompanied by minor spectral variations of the UV-BALs and do not show significant changes in the EW of BALs. These trends suggest a wind-shield scenario where the outflow inclination with respect to the line of sight is decreasing and/or the wind mass is increasing. These changes elevate the covering fraction and/or column densities of the BALs and are likely accompanied by a nearly contemporaneous increase in the column density of the shield.

In the standard paradigm of galaxy formation and evolution, the baryonic component of galaxies forms from the collapse and condensation of gas within dark matter haloes, and later grows from continuous accretion of gaseous mass, both in diffuse form and in mergers with other systems. After a first period of rapid and violent halo growth, the gas settles into a rotationally-supported structure, eventually giving rise to the formation of a stellar disc. Stars evolve and return chemically-processed gas and energy to the interstellar medium, mainly through Type II supernova explosions. In the disc region, the cosmological accretion of gas combines with the outflows resulting from supernovae, affecting the hydrodynamical and structural properties of the disc and producing gas flows in the vertical and radial directions. In this work, we use a simulation of the Auriga Project, a suite of magneto-hydrodynamical, zoom-in cosmological simulations of Milky Way-like galaxies, to study the temporal and radial dependencies of gas accretion onto the disc. We also investigate the disc evolution, focusing on the inside-out disc formation scenario, which is one of the fundamental hypotheses of chemical evolution models of the Galaxy.

The properties of the matter density field present in the initial conditions of a cosmological simulation have an impact on the features of the structures formed after running the simulation. Based on this fact, in this paper we use a random-forest classification algorithm to infer whether or not dark matter particles, traced back to the initial conditions, would end up in dark matter halos whose mass is above some threshold. This problem might be posed as a binary classification task, where the initial conditions of the matter density field are mapped to classification labels provided by a halo finder program. Our results show that random forests are useful tools to predict the output of cosmological simulations without running the full process. These techniques might be used in the future to save computational costs and to explore more efficiently the effect of different dark matter/dark energy candidates on the formation of cosmological structures.

The Search for Extraterrestrial intelligence (SETI) is a scientific and cultural effort seeking evidence of intelligent life beyond earth. Radio SETI observes the radio spectrum for ''technosignatures" that could be produced by an advanced ET society. This work models radio SETI as an end-to-end system, and focuses on narrow-band intentional transmissions. We look at strategies to maximize the expected number of detections per year (DPY) of search. Assuming that ET civilizations will be associated with star systems, we want to maximize the number of stars that may be observed at one time. Assuming a representative star density, this requires maximizing the search volume in a cone defined by the detection range and field of view (FOV). The parameter trades are modified from the case where one simply maximizes signal-to-noise ratio. Instead, a joint optimization between FOV and sensitivity is needed. Some implications: 1) Instead of focusing on the terrestrial microwave window of 1-10 GHz, frequencies below 1 GHz may be optimal for detection rate due to the larger field of view; 2) Arrays of smaller dishes should be favored compared to a single dish of equivalent area; 3) Aperture arrays are desirable due to their large potential FOV. Many radio telescopes under development will provide both high sensitivity and large FOV, and should offer much improved SETI detection rates. Still higher DPY is needed, however, to achieve results in reasonable time horizons, which should be possible by greatly expanding computation capability to the next-generation wide-FOV antenna arrays.

Secondary halo properties beyond mass, such as the mass accretion rate (MAR), concentration, and the half mass scale, are essential in understanding the formation of large-scale structure and dark matter halos. In this paper, we study the impact of secondary halo properties on the galaxy-galaxy lensing observable, $\Delta\Sigma$. We build an emulator trained on N-body simulations to model $\Delta\Sigma$ and quantify the impact of different secondary parameters on the $\Delta\Sigma$ profile. We focus on the impact of MAR on $\Delta\Sigma$. We show that a 3$\sigma$ detection of variations in MAR at fixed halo mass could be achieved with the Hyper Suprime Cam survey in combination with a proxy for MAR with scatter $\sigma_{\Gamma_\mathrm{dyn}|\mathrm{obs}}<1.5$. We show that the full radial profile of $\Delta\Sigma$ depends on secondary properties at fixed halo mass. Consequently, an emulator that can perform full shape fitting yields better than 2 times improvement upon the constraints on MAR than only using the outer part of the halo. Finally, we highlight that miscentering and MAR impact the radial profile of $\Delta\Sigma$ in a similar fashion, implying that miscentering and MAR need to be modeled jointly for unbiased estimates of both effects. We show that present-day lensing data sets have the statistical capability to place constraints on halo MAR. Our analysis opens up new possibilities for observationally measuring the assembly history of the dark matter halos that host galaxies and clusters.

We apply the automated AnomalyFinder algorithm of Paper I (Zang et al. 2021b) to 2018-2019 light curves from the $\simeq 13\,{\rm deg}^2$ covered by the six KMTNet prime fields, with cadences $\Gamma \geq 2\,{\rm hr}^{-1}$. We find a total of 10 planets with mass ratios $q<2\times 10^{-4}$, including five newly discovered planets, one planet that was reported in Paper I, and recovery of four previously discovered planets. One of the new planets, OGLE-2018-BLG-0977Lb, is in a planetary-caustic event, while the other four (OGLE-2018-BLG-0506Lb, OGLE-2018-BLG-0516Lb, OGLE-2019-BLG-1492Lb, and KMT-2019-BLG-0253) are revealed by a ``dip'' in the light curve as the source crosses the host-planet axis on the opposite side of the planet. These subtle signals were missed in previous by-eye searches. The planet-host separations (scaled to the Einstein radius), $s$, and planet-host mass ratios, $q$, are, respectively, $(s,q\times 10^5) = (0.88, 4.1)$, $(0.96\pm 0.10, 8.3)$, $(0.94\pm 0.07, 13)$, $(0.97\pm 0.07, 18)$, and $(0.97\pm0.04,4.1)$, where the ``$\pm$'' indicates a discrete degeneracy. The ten planets are spread out over the range $-5<\log q < -3.7$. Together with the two planets previously reported with $q\sim 10^{-5}$ from the 2018-2019 non-prime KMT fields, this result suggests that planets toward the bottom of this mass-ratio range may be more common than previously believed.

The presence of non-zero helicity in intergalactic magnetic fields is a smoking gun for their primordial origin since they have to be generated by processes that break CP invariance. As an experimental signature for the presence of helical magnetic fields, an estimator $Q$ based on the triple scalar product of the wave-vectors of photons generated in electromagnetic cascades from, e.g., TeV blazars, has been suggested previously. We propose to apply deep learning to helicity classification employing Convolutional Neural Networks and show that this method outperforms the $Q$ estimator.

We implement a test of the variability of the per-cycle annual modulation amplitude in the different phases of the DAMA/LIBRA experiment using Bayesian model comparison. Using frequentist methods, a previous study (Kelso et al 2018) had demonstrated that the DAMA amplitudes spanning over the DAMA/NaI and the first phase of the DAMA/LIBRA phases, show a mild preference for time-dependence in multiple energy bins. With that motivation, we first show using Bayesian techniques that the aforementioned data analyzed in Kelso et al, show a moderate preference for exponentially varying amplitudes in the 2-5 and 2-6 keV energy intervals. We then carry out a similar analysis on the latest modulation amplitudes released by the DAMA collaboration from the first two phases of the upgraded DAMA/LIBRA experiment. We also analyze the single-hit residual rates released by the DAMA collaboration to further look for any possible time-dependency. However, we do not find any evidence for variability of either of the two datasets by using Bayesian model selection. All our analysis codes and datasets have been made publicly available.

Current and future optical and near-infrared wide-field surveys have the potential of finding kilonovae, the optical and infrared counterparts to neutron star mergers, independently of gravitational-wave or high-energy gamma-ray burst triggers. The ability to discover fast and faint transients such as kilonovae largely depends on the area observed, the depth of those observations, the number of re-visits per field in a given time frame, and the filters adopted by the survey; it also depends on the ability to perform rapid follow-up observations to confirm the nature of the transients. In this work, we assess kilonova detectability in existing simulations of the LSST strategy for the Vera C. Rubin Wide Fast Deep survey, with focus on comparing rolling to baseline cadences. Although currently available cadences can enable the detection of more than 300 kilonovae out to 1400 Mpc over the ten-year survey, we can expect only 3-32 kilonovae similar to GW170817 to be recognizable as fast-evolving transients. We also explore the detectability of kilonovae over the plausible parameter space, focusing on viewing angle and ejecta masses. We find that observations in redder izy bands are crucial for identification of nearby (within 300 Mpc) kilonovae that could be spectroscopically classified more easily than more distant sources. Rubin's potential for serendipitous kilonova discovery could be increased by gain of efficiency with the employment of individual 30s exposures (as opposed to 2x15s snap pairs), with the addition of red-band observations coupled with same-night observations in g- or r-bands, and possibly with further development of a new rolling-cadence strategy.

Differential enrichment between $\alpha$- and Fe-peak elements is known to be strongly connected with the shape of the star formation history (SFH), the star formation efficiency (SFE), the inflow and outflow of material, and even the shape of the Initial Mass Function (IMF). However, beyond the Local Group detailed explorations are mostly limited to early-type galaxies due to the lack of a good proxy for [$\alpha$/Fe] in late-type ones, limiting our understanding of the chemical enrichment process. We intent to extend the explorations of [$\alpha$/Fe] to late-type galaxies, in order to understand the details of the differential enrichment process. We compare the gas phase oxygen abundance with the luminosity weighted stellar metallicity in an extensive catalog of $\sim$25,000 H ii regions extracted from the Calar Alto Legacy Integral Field Area (CALIFA) survey, an exploration using integral field spectroscopy of $\sim$900 galaxies, covering a wide range of masses and morphologies. This way we define [O/Fe] as the ratio between both parameters, proposing it as an indirect proxy of the [$\alpha$/Fe] ratio. Results. We illustrate how the [O/Fe] parameter describes the chemical enrichment process in spiral galaxies, finding that: (i) it follows the decreasing pattern with [Fe/H] reported for the [$\alpha$/Fe] ratio and (ii) its absolute scale depends of the stellar mass and the morphology. We reproduce both patterns using two different chemical evolution models (ChEM), considering that galaxies with different stellar mass and morphology present (i) different SFHs, SFEs and different inflow/outflow rates, or (ii) a different maximum stellar mass cut for the IMF. We will explore the differential chemical enrichment using this new proxy galaxy by galaxy and region by region in further studies.

Separation ($\rho$) and Position Angle (PA) measurements are reported of 10 pairs which measures where last reported in the WDS +20 years from epoch of observation 2021.066. Measurements were obtained by direct imaging and are presented with associated measurement uncertainties, as well as, comparisons to measurements determined from Gaia DR2 & EDR3 and historic data extrapolation at epoch of J2000.0.

We present optical follow-up imaging obtained with the Katzman Automatic Imaging Telescope, Las Cumbres Observatory Global Telescope Network, Nickel Telescope, Swope Telescope, and Thacher Telescope of the LIGO/Virgo gravitational wave (GW) signal from the neutron star-black hole (NSBH) merger GW190814. We searched the GW190814 localization region (19 deg$^{2}$ for the 90th percentile best localization), covering a total of 51 deg$^{2}$ and 94.6% of the two-dimensional localization region. Analyzing the properties of 189 transients that we consider as candidate counterparts to the NSBH merger, including their localizations, discovery times from merger, optical spectra, likely host-galaxy redshifts, and photometric evolution, we conclude that none of these objects are likely to be associated with GW190814. Based on this finding, we consider the likely optical properties of an electromagnetic counterpart to GW190814, including possible kilonovae and short gamma-ray burst afterglows. Using the joint limits from our follow-up imaging, we conclude that a counterpart with an $r$-band decline rate of 0.68 mag day$^{-1}$, similar to the kilonova AT 2017gfo, could peak at an absolute magnitude of at most $-17.8$ mag (50% confidence). Our data are not constraining for ''red'' kilonovae and rule out ''blue'' kilonovae with $M>0.5 M_{\odot}$ (30% confidence). We strongly rule out all known types of short gamma-ray burst afterglows with viewing angles $<$17$^{\circ}$ assuming an initial jet opening angle of $\sim$$5.2^{\circ}$ and explosion energies and circumburst densities similar to afterglows explored in the literature. Finally, we explore the possibility that GW190814 merged in the disk of an active galactic nucleus, of which we find four in the localization region, but we do not find any candidate counterparts among these sources.

Significant clustering around the rarest luminous quasars is a feature predicted by dark matter theory combined with number density matching arguments. However, this expectation is not reflected by observations of quasars residing in a diverse range of environments. Here, we assess the tension in the diverse clustering of visible $i$-band dropout galaxies around luminous $z\sim6$ quasars. Our approach uses a simple empirical method to derive the median luminosity to halo mass relation, $L_{c}(M_{h})$ for both quasars and galaxies under the assumption of log-normal luminosity scatter, $\Sigma_{Q}$ and $\Sigma_{G}$. We show that higher $\Sigma_{Q}$ reduces the average halo mass hosting a quasar of a given luminosity, thus introducing at least a partial reversion to the mean in the number count distribution of nearby Lyman-Break galaxies. We generate a large sample of mock Hubble Space Telescope fields-of-view centred across rare $z\sim6$ quasars by resampling pencil beams traced through the dark matter component of the BlueTides cosmological simulation. We find that diverse quasar environments are expected for $\Sigma_{Q}>0.4$, consistent with numerous observations and theoretical studies. However, we note that the average number of galaxies around the central quasar is primarily driven by galaxy evolutionary processes in neighbouring halos, as embodied by our parameter $\Sigma_{G}$, instead of a difference in the large scale structure around the central quasar host, embodied by $\Sigma_{Q}$. We conclude that models with $\Sigma_{G}>0.3$ are consistent with current observational constraints on high-z quasars, and that such a value is comparable to the scatter estimated from hydrodynamical simulations of galaxy formation.

Exoplanets observed by the {\it Kepler} telescope exhibit a bi-modal, radius distribution, which is known as the radius gap. We explore an origin of the radius gap, focusing on multi-planet systems. Our simple theoretical argument predicts that type I planetary migration produces different configurations of protoplanets with different masses and such different configurations can result in two distinguishable populations of small-sized multi-planet systems. We then perform an observational analysis to verify this prediction. In the analysis, multiple Kolmogorov-Smirnov tests are applied to the observed systems, using the statistical measures that are devised to systematically characterize the properties of multi-planet systems. We find with 99.5\% confidence that the observed, small-sized multi-planet systems are divided into two distinct populations. The distinction likely originates from different spatial distributions of protoplanets, which are determined by type I migration and subsequently trigger giant impact. We also show that these distinct populations are separated around the radius gap when the gas surface density of protoplanetary disks is $\sim 10^2$ g cm$^{-2}$ in the vicinity of the host stars. This work therefore emphasizes the importance of planetary migration and the inner disk properties.

Given their prominent role in galaxy evolution, it is of paramount importance to unveil galaxy interactions and merger events and to investigate the underlying mechanisms. The use of high-resolution data makes it easier to identify merging systems, but it can still be challenging when the morphology does not show any clear galaxy-pair or gas bridge. Characterising the origin of puzzling kinematic features can help to reveal complicated systems. Here, we present a merging galaxy, MaNGA 1-114955, in which we highlighted the superimposition of two distinct rotating discs along the line of sight. These counter-rotating objects both lie on the star-forming main sequence but display perturbed stellar velocity dispersions. The main galaxy presents off-centred star formation as well as off-centred high-metallicity regions supporting the scenario of recent starbursts, while the secondary galaxy hosts a central starburst which coincides with an extended radio emission, in excess with respect to star formation expectations. Stellar mass as well as dynamical mass estimates agree towards a mass ratio within the visible radius of 9:1 for these interacting galaxies. We suggest we are observing a pre-coalescence stage of a merger. The primary galaxy has accreted gas through a past first pericentre passage about 1 Gyr ago, and more recently from the secondary gas-rich galaxy, which exhibits an underlying active galactic nucleus (AGN). Our results demonstrate how a galaxy can hide another one and the relevance of a multi-component approach to study ambiguous systems. We anticipate our method to be efficient at unveiling the mechanisms taking place in a sub-sample of galaxies observed by the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, all exhibiting kinematic features of puzzling origin in their gas emission lines.

Matter distribution models of the Milky Way galaxy are usually stationary, although there are known to be wave-like perturbations in the disc at $\sim10\%$ level of the total density. Modelling of the overall acceleration field by allowing non-equilibrium is a complicated task. We must learn to distinguish whether density enhancements are persistent or not by their nature. In the present paper, we elaborate our orbital arc method to include the effects of massless perturbations and non-stationarities in the modelling. The method is tested by modelling of simulation data and shown to be valid. We apply the method to the Gaia DR 2 data within a region of $\sim 0.5$ kpc from the Sun and confirm that acceleration field in the Solar Neighbourhood has a perturbed nature -- the phase space density along the orbits of stars grow in the order of $h\lesssim 5\%$ per Myr due to non-stationarity. This result is a temporally local value and can be used only within the timeframe of a few Myrs. An attempt to pinpoint the origin of the perturbation shows that the stars having larger absolute angular momentum are the main carriers of the local perturbation. As they are faster than the average thin disc star, they are either originating further away and are close in their pericentre or they are perturbed locally by a fast co-moving perturber, such as gas disc inhomogenities.

Recent observations using several different telescopes and sky surveys showed patterns of asymmetry in the distribution of galaxies by their spin directions as observed from Earth. These studies were done with data imaged from the Northern hemisphere, showing excellent agreement between different telescopes and different analysis methods. Here, data from the DESI Legacy Survey was used. The initial dataset contains $\sim2.2\cdot10^7$ galaxy images, reduced to $\sim8.1\cdot10^5$ galaxies annotated by their spin direction using a symmetric algorithm. That makes it not just the first analysis of its kind in which the majority of the galaxies are in the Southern hemisphere, but also by far the largest dataset used for this purpose to date. The results show strong agreement between opposite parts of the sky, such that the asymmetry in one part of the sky is similar to the inverse asymmetry in the corresponding part of the sky in the opposite hemisphere. Fitting the distribution of galaxy spin directions to cosine dependence shows a dipole axis with probability of 4.66$\sigma$. Interestingly, the location of the most likely axis is within close proximity to the CMB Cold Spot. The profile of the distribution is nearly identical to the asymmetry profile of the distribution identified in Pan-STARRS, and it is within 1$\sigma$ difference from the distribution profile in SDSS and HST. All four telescopes show similar large-scale profile of asymmetry.

We present an ensemble X-ray analysis of systematic perturbations in the central hot gas properties for a sample of 28 nearby strong cool-core systems selected from the HIghest X-ray FLUx Galaxy Cluster Sample (HIFLUGCS). We analyze their cool-core features observed with the Chandra X-ray Observatory. All individual systems in our sample exhibit at least a pair of positive and negative excess perturbations in the X-ray residual image after subtracting the global brightness profile. We extract and analyze X-ray spectra of the intracluster medium (ICM) in the detected perturbed regions. To investigate possible origins of the gas perturbations, we characterize thermodynamic properties of the ICM in the perturbed regions and characterize their correlations between positive and negative excess regions. The best-fit relations for temperature and entropy show a clear offset from the one-to-one relation, $T_\mathrm{neg}/T_\mathrm{pos}=1.20^{+0.04}_{-0.03}$ and $K_\mathrm{neg}/K_\mathrm{pos}=1.43\pm 0.07$, whereas the best-fit relation for pressure is found to be remarkably consistent with the one-to-one relation $P_\mathrm{neg}=P_\mathrm{pos}$, indicating that the ICM in the perturbed regions is in pressure equilibrium. These observed features in the HIFLUGCS sample are in agreement with the hypothesis that the gas perturbations in cool cores are generated by gas sloshing. We also analyze synthetic observations of perturbed cluster cores created from binary merger simulations, finding that the observed temperature ratio agrees with the simulations, $T_\mathrm{neg}/T_\mathrm{pos}\sim 1.3$. We conclude that gas sloshing induced by infalling substructures plays a major role in producing the characteristic gas perturbations in cool cores. The ubiquitous presence of gas perturbations in cool cores may suggest a significant contribution of gas sloshing to suppressing runaway cooling of the ICM.

Gamma-ray bursts - the most luminous explosions in the Universe - are produced as a result of cataclysmic events such as the collapse of a massive star or the merger of two neutron stars. We monitored the position of the close-by (about 370 Megaparsecs) gamma-ray burst GRB~190829A, which originated from a massive star collapse, through very long baseline interferometry (VLBI) observations with the European VLBI Network and the Very Long Baseline Array, involving a total of 30 telescopes across four continents. We carried out a total of 9 observations between 9 and 117 days after the gamma-ray burst at 5 and 15 GHz, reaching an overall excellent resolution. From a state-of-the art analysis of these data, we obtained valuable information on the source size and expansion rate. The measurements are in remarkable agreement with the size evolution entailed by a detailed modelling of the multi-wavelength light curves with a forward plus reverse shock model, which agrees with the observations across almost 18 orders of magnitude in frequency (including the High Energy Stereoscopic System data at teraelectronvolt photon energies) and more than 4 orders of magnitude in time. Thanks to the multi-wavelength, high-cadence coverage of the afterglow, inherent degeneracies in the afterglow model are broken to a large extent, allowing us to capture some unique physical insights: we find a low prompt emission efficiency $\lesssim 10^{-3}$; we constrain the fraction of electrons that are accelerated to relativistic speeds in the forward shock downstream to be $\chi_e<13\%$ at the 90\% confidence level; we find that the magnetic field energy density in the reverse shock downstream must decay rapidly after the shock crossing.

When binary black holes merge in dense star clusters, their remnants can pair up with other black holes in the cluster, forming heavier and heavier black holes in a process called hierarchical merger. The most important condition for hierarchical merger to occur is that remnants formed by mergers are retained by the host star cluster. Using the publicly available gravitational-wave event database, we infer the magnitudes of kick velocities imparted to the remnant black holes due to anisotropic emission of gravitational waves and use that to quantify the retention probability of each event as a function of the escape speed of the star cluster. Among the GWTC-2 events, GW190814 provides the tightest constraint on the kick magnitude with ${\rm V_{kick}}=73_{-8}^{+11}$ km/s at 90\% credible level. We find that star clusters with escape speeds of 500 km/s can retain about 50\% of the events in the second gravitational-wave transient catalog (GWTC-2). Using the escape speed distributions of nuclear star clusters and globular clusters, we find that $\sim 8$ (1) remnants of GWTC-2 may be retained by the host star cluster if all GWTC-2 events occurred in nuclear (globular) clusters. Our study demonstrates the importance of folding in kick velocity inferences in future studies of hierarchical mergers.

Stellar prolate rotation in dwarf galaxies is rather uncommon, with only two known galaxies in the Local Group showing such feature (Phoenix and And II). Cosmological simulations show that in massive early-type galaxies prolate rotation likely arises from major mergers. However, the origin of such kinematics in the dwarf galaxies regime has only been explored using idealized simulations. Here we made use of hydrodynamical cosmological simulations of dwarfs galaxies with stellar mass between $3\times10^5$ and $5\times10^8$ M$_{\odot}$ to explore the formation of prolate rotators. Out of $27$ dwarfs, only one system showed clear rotation around the major axis, whose culprit is a major merger at $z=1.64$, which caused the transition from an oblate to a prolate configuration. Interestingly, this galaxy displays a steep metallicity gradient, reminiscent of the one measured in Phoenix and And II: this is the outcome of the merger event that dynamically heats old, metal-poor stars, and of the centrally concentrated residual star formation. Major mergers in dwarf galaxies offer a viable explanation for the formation of such peculiar systems, characterized by steep metallicity gradients and prolate rotation.

In this dissertation, several components of large-scale solar flows are studied observationally: solar equatorial Rossby waves (waves of radial vorticity), large-scale convection, and surface flows around active regions. Maps of horizontal flows are derived from photospheric observations by the Helioseismic and Magnetic Imager (HMI) aboard the Solar Dynamics Observatory (SDO) using two different techniques: granulation tracking and local helioseismology. First, the eigenfunctions of solar Rossby waves are measured from helioseismic ring-diagram flow maps with a correlation method and a spectral analysis. Down to $9$ Mm below the surface, the dependence of the radial vorticity with radius $r$ is consistent with $r^{m-1}$, for a given longitudinal wavenumber $m$. At the surface, the eigenfunctions are complex-valued. The real part decreases away from the equator and switches sign around $\pm 20-30^\circ$. The imaginary part is small, but nonzero, and may be due to wave attenuation. This may have implications for the transport of angular momentum in the latitudinal direction. Second, we revisit previous measurements of power spectra of longitudinal velocities near the solar surface, obtained from time-distance and ring-diagram helioseismology. Several issues in these past helioseismic analyses are identified and corrected. The corrections are not sufficient to remove the discrepancy between the measurements. I thus present new velocity power spectra from granulation tracking and ring-diagram helioseismology. The two new measurements are close to each other near the solar surface, and the corresponding kinetic energy decreases with increasing spatial scale.

The accuracy of theoretical mass, radius and effective temperature values for M-dwarf stars is an active topic of debate. Differences between observed and theoretical values have raised the possibility that current theoretical stellar structure and evolution models are inaccurate towards the low-mass end of the main sequence. To explore this issue we use the CHEOPS satellite to obtain high-precision light curves of eclipsing binaries with low mass stellar companions. We use these light curves combined with the spectroscopic orbit for the solar-type companion to measure the mass, radius and effective temperature of the M-dwarf star. Here we present the analysis of three eclipsing binaries. We use the pycheops data analysis software to fit the observed transit and eclipse events of each system. Two of our systems were also observed by the TESS satellite -- we similarly analyse these light curves for comparison. We find consistent results between CHEOPS and TESS, presenting three stellar radii and two stellar effective temperature values of low-mass stellar objects. These initial results from our on-going observing programme with CHEOPS show that we can expect to have ~24 new mass, radius and effective temperature measurements for very low mass stars within the next few years.

The dark matter halo surface density, given by the product of the dark matter core radius ($r_c$) and core density ($\rho_c$) has been shown to be a constant for a wide range of isolated galaxy systems. Here, we carry out a test of this {\em ansatz} using a sample of 17 relaxed galaxy groups observed using Chandra and XMM-Newton, as an extension of our previous analysis with galaxy clusters. We find that $\rho_c \propto r_c^{-1.35^{+0.16}_{-0.17}}$, with an intrinsic scatter of about 27%, which is about 1.5 times larger than that seen for galaxy clusters. Our results therefore indicate that the surface density is discrepant with respect to a constant value by about 2$\sigma$. Furthermore, we also implement a test of the radial acceleration relation for this group sample. We find that the residual scatter in the radial acceleration relation is about 0.32 dex and a factor of three larger than that obtained using galaxy clusters. The acceleration scale which we obtain is in-between that seen for single galaxies and clusters.

Aims: Even the most-precise radial-velocity instruments gather high-resolution spectra that present systematic errors that a data reduction pipeline cannot identify and correct for efficiently. In this paper, we aim at improving the radial-velocity precision of HARPS measurements by cleaning individual extracted spectra using the wealth of information contained in spectra time-series. Methods: We developed YARARA, a post-processing pipeline designed to clean high-resolution spectra from instrumental systematics and atmospheric contamination. Spectra are corrected for: tellurics, interference pattern, detector stitching, ghosts and fiber B contaminations as well as more advanced spectral line-by-line corrections. YARARA uses Principal Component Analysis on spectra time-series with prior information to disentangle contaminations from real Doppler shifts. We applied YARARA on three systems: HD10700, HD215152 and HD10180 and compared our results to the HARPS standard Data Reduction Software and the SERVAL post-processing pipeline. Results: On HD10700, we obtain radial-velocity measurements that present a rms smaller than 1 m/s over the 13 years of the HARPS observations, which is 20 and 10 % better than the HARPS Data Reduction Software and the SERVAL post-processing pipeline, respectively. We also injected simulated planets on the data of HD10700 and demonstrated that YARARA does not alter pure Doppler shifted signals. On HD215152, we demonstrated that the 1-year signal visible in the periodogram becomes marginal after processing with YARARA and that the signals of the known planets become more significant. Finally, on HD10180, the known six exoplanets are well recovered although different orbitals parameters and planetary masses are provided by the new reduced spectra.

Galaxies experiencing intense star-formation episodes are expected to be rich in energetic cosmic rays (CRs). These CRs undergo hadronic interactions with the interstellar gases of their host to drive $\gamma$-ray emission, which has already been detected from several nearby starbursts. Unresolved $\gamma$-ray emission from more distant star-forming galaxies (SFGs) is expected to contribute to the extra-galactic $\gamma$-ray background (EGB). However, despite the wealth of high-quality all-sky data from the Fermi-LAT $\gamma$-ray space telescope collected over more than a decade of operation, the exact contribution of such SFGs to the EGB remains unsettled. We investigate the high-energy $\gamma$-ray emission from SFGs up to redshift $z=3$ above a GeV, and assess the contribution they can make to the EGB. We show the $\gamma$-ray emission spectrum from a SFG population can be determined from just a small number of key parameters, from which we model a range of possible EGB realisations. We demonstrate that populations of SFGs leave anisotropic signatures in the EGB, and that these can be accessed using the spatial power spectrum. Moreover, we show that such signatures will be accessible with ongoing operation of current $\gamma$-ray instruments, and detection prospects will be greatly improved by the next generation of $\gamma$-ray observatories, in particular the Cherenkov Telescope Array.

We analyze the possible effect of rings on orbital velocities in galaxies. The superposition of the central force with the gravitational forces induced by the rings opens up various possibilities of the course of orbital velocities. The orbital velocity depends on the position of the star in the ring. We illustrate this dependence on several models, where we show the course of potential curves and the curves of field strength.

The recent advanced LIGO/Virgo detections of gravitational waves (GWs) from stellar binary black hole (BBH) mergers, in particular GW190521, which is potentially associated with a quasar, have stimulated renewed interest in active galactic nuclei (AGNs) as factories of merging BBHs. Compact objects evolving from massive stars are unavoidably enshrouded by a massive envelope to form accretion-modified stars (AMSs) in the dense gaseous environment of a supermassive black hole (SMBH) accretion disk. We show that most AMSs form binaries due to gravitational interaction with each other during radial migration in the SMBH disk, forming BBHs inside the AMS. When a BBH is born, its orbit is initially governed by the tidal torque of the SMBH. Bondi accretion onto BBH at a hyper-Eddington rate naturally develops and then controls the evolution of its orbits. We find that Bondi accretion leads to efficient removal of orbital angular momentum of the binary, whose final merger produces a GW burst. Meanwhile, the Blandford-Znajek mechanism pumps the spin energy of the merged BH to produce an electromagnetic counterpart (EMC). Moreover, hyper-Eddington accretion onto the BBH develops powerful outflows and triggers a Bondi explosion, which manifests itself as a EMC of the GW burst, depending on the viscosity of the accretion flow. Thermal emission from Bondi sphere appears as one of EMCs. BBHs radiate GWs with frequencies $\sim 10^{2}\,$Hz, which are accessible to LIGO.

This article describes our approach to quantifying the characteristics of meteors such as temperature, chemical composition, and others. We are using a new approach based on colourimetry. We analyze an image of Leonid meteor-6230 obtained by Mike Hankey in 2012. Analysis of the temporal features of the meteoroid trail is performed. For determining the meteor characteristics we use the "tuning technique" in combination with a simulation model of intrusion. The progenitor of the meteor was found as an object weighing 900 kg at a speed of 36.5 km/s. The meteoroid reached a critical value of the pressure at an altitude of about 29 km in a time of about 4.6 sec with a residual mass of about 20 kg, and a residual speed of about 28 km/s. At this moment, a meteoroid exploded and destroyed. We use the meteor multicolour light curves revealed from a DSLR image in the RGB colour standard. We switch from the RGB colour system to Johnson's RVB colour system introducing colour corrections. This allows one to determine the colour characteristics of the meteor radiation. We are using a new approach based on colourimetry. Colourimetry of BGR three-beam light curves allows the identification of the brightest spectral lines. Our approach based on colourimetry allows direct measurements of temperature in the meteor trail. We find a part of the trajectory where the meteoroid radiates as an absolutely black body. The R/G and B/G light curves ratio allow one to identify the wavelengths of the emission lines using the transmission curves of the RGB filters. At the end of the trajectory, the meteoroid radiates in the lines Ca II H, K 393, 397 nm, Fe I 382, 405 nm, Mg I 517 nm, Na I 589 nm, as well as atmospheric O I 779 nm.

V838 Mon erupted in 2002 quickly becoming the prototype of a new type of stellar eruptions known today as (luminous) red novae. The red nova outbursts are thought to be caused by stellar mergers. The merger in V838 Mon took place in a triple or higher system involving two B-type stars. We mapped the merger site with ALMA at a resolution of 25 mas in continuum dust emission and in rotational lines of simple molecules, including CO, SiO, SO, SO$_2$, AlOH, and H$_2$S. We use radiative transfer calculations to reproduce the remnant's architecture at the epoch of the ALMA observations. For the first time, we identify the position of the B-type companion relative to the outbursting component of V838 Mon. The stellar remnant is surrounded by a clumpy wind with characteristics similar to winds of red supergiants. The merger product is also associated with an elongated structure, $17.6 \times 7.6$ mas, seen in continuum emission, and which we interpret as a disk seen at a moderate inclination. Maps of continuum and molecular emission show also a complex region of interaction between the B-type star (its gravity, radiation, and wind) and the flow of matter ejected in 2002. The remnant's molecular mass is about 0.1 M$_{\odot}$ and the dust mass is 8.3$\cdot$10$^{-3}$ M$_{\odot}$. The mass of the atomic component remains unconstrained. The most interesting region for understanding the merger of V838 Mon remains unresolved but appears elongated. To study it further in more detail will require even higher angular resolutions. ALMA maps show us an extreme form of interaction between the merger ejecta with a distant (250 au) companion. This interaction is similar to that known from the Antares AB system but at a much higher mass loss rate. The B-type star not only deflects the merger ejecta but also changes its chemical composition with an involvement of circumstellar shocks.

We report ALMA observations of the dust continuum and {\cii} emission of the host galaxy of J0439+1634, a gravitationally lensed quasar at $z=6.5$. Gravitational lensing boosts the source-plane resolution to $\sim0\farcs15$ $(\sim0.8\text{ kpc})$. The lensing model derived from the ALMA data is consistent with the fiducial model in \citet{fan19} based on {\it HST} imaging. The host galaxy of J0439+1634 can be well-fitted by a S\'ersic profile consistent with an exponential disk, both in the far-infrared (FIR) continuum and the {\cii} emission. The overall magnification is $4.53\pm0.05$ for the continuum and $3.44\pm0.05$ for the {\cii} line. The host galaxy of J0439+1634 is a compact ultra-luminous infrared galaxy, with a total star formation rate (SFR) of $1.56\times10^{3}M_\odot/\text{year}$ after correcting for lensing and an effective radius of $0.74$ kpc. The resolved regions in J0439+1634 follow the ``{\cii} deficit," where the {\cii}-to-FIR ratio decreases with FIR surface brightness. The reconstructed velocity field of J0439+1634 appears to be rotation-like. The maximum line-of-sight rotation velocity of 130 km/s at a radius of 2 kpc. However, our data cannot be fit by an axisymmetric thin rotating disk, and the inclination of the rotation axis, $i$, remains unconstrained. We estimate the dynamical mass of the host galaxy to be $7.9\sin^{-2}(i)\times10^{9}M_\odot$. J0439+1634 is likely to have a high gas-mass fraction and an oversized SMBH compared to local relations. The SFR of J0439+1634 reaches the maximum possible values, and the SFR surface density is close to the highest value seen in any star-forming galaxy currently known in the universe.

55 Cnc e is a transiting super-Earth (radius $1.88\rm\,R_\oplus$ and mass $8\rm\, M_\oplus$) orbiting a G8V host star on a 17-hour orbit. Spitzer observations of the planet's phase curve at 4.5 $\mu$m revealed a time-varying occultation depth, and MOST optical observations are consistent with a time-varying phase curve amplitude and phase offset of maximum light. Both broadband and high-resolution spectroscopic analyses are consistent with either a high mean molecular weight atmosphere or no atmosphere for planet e. A long term photometric monitoring campaign on an independent optical telescope is needed to probe the variability in this system. We seek to measure the phase variations of 55 Cnc e with a broadband optical filter with the 30 cm effective aperture space telescope CHEOPS and explore how the precision photometry narrows down the range of possible scenarios. We observed 55 Cnc for 1.6 orbital phases in March of 2020. We designed a phase curve detrending toolkit for CHEOPS photometry which allows us to study the underlying flux variations of the 55 Cnc system. We detected a phase variation with a full-amplitude of $72 \pm 7$ ppm but do not detect a significant secondary eclipse of the planet. The shape of the phase variation resembles that of a piecewise-Lambertian, however the non-detection of the planetary secondary eclipse, and the large amplitude of the variations exclude reflection from the planetary surface as a possible origin of the observed phase variations. They are also likely incompatible with magnetospheric interactions between the star and planet but may imply that circumplanetary or circumstellar material modulate the flux of the system. Further precision photometry of 55 Cnc from CHEOPS will measure variations in the phase curve amplitude and shape over time this year.

We study cosmic evolution based on the fixed points in the dynamical analysis of the Degenerate Higher-Order Scalar-Tensor (DHOST) theories. We consider the DHOST theories in which the propagation speed of gravitational waves is equal to the speed of light, the tensor perturbations do not decay to dark energy perturbations, and the scaling solutions exist. The scaling fixed point associated with late time acceleration of universe can be either stable or saddle depending on the parameters of the theory. For some ranges of the parameters, the accelerated scaling point and the field dominated point can be simultaneously stable. Cosmic evolution will reach the accelerated scaling point if the time derivative of the scalar field in the theory is positive during the matter domination. If the time derivative of the scalar field is negative during the matter domination, the background universe will evolve towards the field dominated point. The density parameter of the matter can be larger than unity before reaching the scaling fixed point if the deviation from the Einstein theory of gravity is too large and the initial conditions for the dynamical variables during the matter domination are significantly different from the accelerated scaling point. The stabilities of $\phi$MDE fixed point are similar to the coupled dark energy models. In our consideration, the universe can only evolve from the $\phi$MDE to the field dominated point.

Primordial black holes may have been produced in the early stages of the thermal history of the Universe after cosmic inflation. If so, dark matter in the form of elementary particles can be subsequently accreted around these objects, in particular when it gets non-relativistic and further streams freely in the primordial plasma. A dark matter mini-spike builds up gradually around each black hole, with density orders of magnitude larger than the cosmological one. We improve upon previous work by carefully inspecting the computation of the mini-spike radial profile as a function of black hole mass, dark matter particle mass and temperature of kinetic decoupling. We identify a phase-space contribution that has been overlooked and that leads to changes in the final results. We also derive complementary analytical formulae using convenient asymptotic regimes, which allows us to bring out peculiar power-law behaviors for which we provide tentative physical explanations.

We present a detailed photometric and kinematical analysis of poorly studied open cluster IC 1434 using CCD VRI, APASS, and Gaia DR2 database for the first time. BY determining the membership probability of stars, we identified the 238 most probable members with a probability higher than 60% by using proper motion and parallax data as taken from the Gaia DR2 catalog. The mean proper motion of the cluster is obtained as emu_x= - 3.89 +/- 0.19 and emu_y= - 3.34 +/- 0.19 mas/yr in both the directions of right ascension and declination. The radial distribution of member stars provides cluster extent as 7.6 arcmin. We have estimated the interstellar reddening (E(B-V)) as 0.34 mag using the transformation equations from literature. We obtained the values of cluster age and distance are 631 +/- 73 Myr and 3.2 +/- 0.1 Kpc.

We propose a common envelope evolution (CEE) scenario where a red giant branch (RGB) star engulfs a planet during its core helium flash to explain the puzzling system WD 1856+534 where a planet orbits a white dwarf (WD) of mass 0.52Mo with an orbital period of 1.4 day. At the heart of the scenario is the recently proposed assumption that the vigorous convection that core helium flash of RGB stars drive in the core excite waves that propagate and deposit their energy in the envelope. Using the binary-MESA stellar evolution code we show that this energy deposition substantially reduces the binding energy of the envelope and causes its expansion. We propose that in some cases RGB stars might engulf massive planets of ~0.01Mo during their core helium flash phase, and that the planet can unbind most of the mass of the bloated envelope. We show that there is a large range of initial orbital radii for which this scenario might take place under our assumptions. This scenario is relevant to other systems of close sub-stellar objects orbiting white dwarfs, like the brown dwarf-WD system ZTFJ003855.0+203025.5.

We propose an inflationary primordial feature model that can explain both the large and small-scale anomalies in the currently measured cosmic microwave background (CMB) anisotropy spectra, revealing a clip of adventurous history of the Universe during its primordial epoch. Although the model is currently statistically indistinguishable from the Standard Model, we show that future observations such as the Simons Observatory and LiteBIRD will complement each other in distinguishing the model differences due to their accurate E-mode polarization measurements, and the PICO mission, if funded, can put stringent constraints on all characteristic properties. The model predicts a signal of classical primordial standard clock, which can also be used to distinguish the inflation and alternative scenarios in a model-independent fashion.

Geometry of relativistic jets in active galaxies provides important information about mechanisms of launching, collimation and acceleration of plasma flow. We propose a new method to probe a boundary shape of a jet on parsec scales -- in the vicinity of its radio core. Apparent speed of an outflow is derived from variability time delays and core shifts measured at the same jet region, providing a self-consistent estimate of the Lorentz factor $\Gamma$. We link together the distance along the jet z with its transverse size assuming a constant flow acceleration. Our results indicate that jets have parabolic shape and sustain an effective acceleration in the core region, consistent with the Lorentz factor dependency $\Gamma\propto z^{0.5}$. The proposed method can be applied to the sources observed at small viewing angles as well as to the distant sources when direct measurements are impossible due to a limited angular resolution.

Eruptions of coronal mass ejections (CMEs) from the Sun are usually associated with a number of signatures that can be identified in solar disc imagery. However, there are cases in which a CME that is well observed in coronagraph data is missing a clear low-coronal counterpart. These events have received attention during recent years, mainly as a result of the increased availability of multi-point observations, and are now known as 'stealth CMEs'. In this work, we analyse examples of stealth CMEs featuring various levels of ambiguity. All the selected case studies produced a large-scale CME detected by coronagraphs and were observed from at least one secondary viewpoint, enabling a priori knowledge of their approximate source region. To each event, we apply several image processing and geometric techniques with the aim to evaluate whether such methods can provide additional information compared to the study of "normal" intensity images. We are able to identify at least weak eruptive signatures for all events upon careful investigation of remote-sensing data, noting that differently processed images may be needed to properly interpret and analyse elusive observations. We also find that the effectiveness of geometric techniques strongly depends on the CME propagation direction with respect to the observers and the relative spacecraft separation. Being able to observe and therefore forecast stealth CMEs is of great importance in the context of space weather, since such events are occasionally the solar counterparts of so-called 'problem geomagnetic storms'.

In the standard (classic) approach, galaxy clustering measurements from spectroscopic surveys are compressed into baryon acoustic oscillations and redshift space distortions measurements, which in turn can be compared to cosmological models. Recent works have shown that avoiding this intermediate step and fitting directly the full power spectrum signal (full modelling) leads to much tighter constraints on cosmological parameters. Here we show where this extra information is coming from and extend the classic approach with one additional effective parameter, such that it captures, effectively, the same amount of information as the full modelling approach, but in a model-independent way. We validate this new method (ShapeFit) on mock catalogs, and compare its performance to the full modelling approach finding both to deliver equivalent results. The ShapeFit extension of the classic approach promotes the standard analyses at the level of full modelling ones in terms of information content, with the advantages of i) being more model independent; ii) offering an understanding of the origin of the extra cosmological information; iii) allowing a robust control on the impact of observational systematics.

We estimate galaxy clustering under a modified gravitational potential. In particular, the modifications in gravitational potential energy occur due to a power-law and cosmological constant terms. We derive a canonical partition function for the system of galaxies interacting under such a modified gravitational potential. Moreover, we compute various thermodynamical equation of states for the system. We do comparative analysis in order to emphasize the effect of corrections on thermodynamics of the system. Interestingly, the modifications in thermodynamical quantities are embedded in clustering parameter only.

Recently, the LIGO-Virgo Collaboration (LVC) concluded that there is no evidence for lensed gravitational waves (GW) in the first half of the O3 run, claiming "We find the observation of lensed events to be unlikely, with the fractional rate at $\mu>2$ being $3.3\times 10^{-4}$". While we agree that the chance of an individual GW event being lensed at $\mu>2$ is smaller than $10^{-3}$, the number of observed events depends on the product of this small probability times the rate of mergers at high redshift. Observational constraints from the stochastic GW background indicate that the rate of conventional mass BBH mergers (8 < M (M$_{\odot}$) < 15) in the redshift range 1<z< 2 could be as high as O($10^7$) events per year, more than sufficient to compensate for the intrinsically low probability of lensing. To reach the LVC trigger threshold these events require high magnification, but would still produce up to 10 to 30 LVC observable events per year. Thus, all the LVC observed ordinary stellar mass BBH mergers from this epoch must be strongly lensed. By adopting low-rates at high redshift, LVC assumes that lensed events can not be taking place, thus incorrectly assigning them a closer distance and higher masses by a factor of a few (typically 2 to 5). The LVC adopted priors on time delay are in tension with the distribution of observed time delays in lensed quasars. Pairs of events like GW190421-GW190910 and GW190424-GW190910, which are directly assigned a probability of zero by LVC, should be instead considered as prime candidates to be strongly lensed GW pairs, since their separation in time is consistent with observations of time delays in lensed quasars. Correcting for the LVC wrong Bayesian priors, maximum merger rate of conventional mass BBH in 1<z<2, and gravitational lensing time-delay model, reverses the LVC conclusions and supports the strong gravitational lensing hypothesis.

We re-evaluate the status of supersonic electroweak baryogenesis using a generalized fluid Ansatz for the non-equilibrium distribution functions. Instead of truncating the expansion to first order in momentum, we allow for higher order terms as well, including up to 21 fluctuations. The collision terms are computed analytically at leading-log accuracy. We also point out inconsistencies in the standard treatments of transport in electroweak baryogenesis, arguing that one cannot do without specifying an Ansatz for the distribution function. We present the first analysis of baryogenesis using the fluid approximation to higher orders. Our results support the recent findings that baryogenesis may indeed be possible even in the presence of supersonic wall velocities.

Tidal perturbations play an important role in the study of the dynamics in the classical two-body system. Understanding tidal effects in strong-field regions may allow one to use gravitational-wave or electromagnetic observations to locate or constraint the location of possible companions. Here, we investigate how timelike and null geodesics of a Schwarzschild black hole are affected in the presence of a companion. There is a panoply of new effects. In some limiting cases, we find analytical solutions for closed null or timelike geodesics. Our results show that light ring period as measured by a far-away observer can be eiter shorter or longer, depending on the location of the companion. We also show that there are closed lightlike trajectories which are elliptic (for equatorial companions), and that timelike particles are affected in a similar manner. Finally, we attempt at estimating the ringdown from tidally perturbed geometries. Our results indicate that there are two stages in the relaxation of such geometries, one associated with a prompt decay of waves around the deformed photonsphere, and a later relaxation of the global geometry. These results are consistent with previous, full numerical studies.

Gravitational-wave cosmology began in 2017 with the observation of the gravitational waves emitted in the merger of two neutron stars, and the coincident observation of the electromagnetic emission that followed. Although only a $30\%$ measurement of the Hubble constant was achieved, future observations may yield more precise measurements either through other coincident events or through cross correlation of gravitational-wave events with galaxy catalogs. Here, we implement a new way to measure the Hubble constant without an electromagnetic counterpart and through the use of the binary Love relations. These relations govern the tidal deformabilities of neutron stars in an equation-of-state insensitive way. Importantly, the Love relations depend on the component masses of the binary in the source frame. Since the gravitational-wave phase and amplitude depend on the chirp mass in the observer (and hence redshifted) frame, one can in principle combine the binary Love relations with the gravitational-wave data to directly measure the redshift, and thereby infer the value of the Hubble constant. We implement this approach in both real and synthetic data through a Bayesian parameter estimation study in a range of observing scenarios. We find that for the LIGO/Virgo/KAGRA design sensitivity era, this method results in a similar measurement accuracy of the Hubble constant to those of current-day, dark-siren measurements. For third generation detectors, this accuracy improves to $\lesssim 10\%$ when combining measurements from binary neutron star events in the LIGO Voyager era, and to $\lesssim 2\%$ in the Cosmic Explorer era.

The dual-phase xenon time projection chamber (TPC) is a powerful tool for direct-detection experiments searching for WIMP dark matter, other dark matter models, and neutrinoless double-beta decay. Successful operation of such a TPC is critically dependent on the ability to hold high electric fields in the bulk liquid, across the liquid surface, and in the gas. Careful design and construction of the electrodes used to establish these fields is therefore required. We present the design and production of the LUX-ZEPLIN (LZ) experiment's high-voltage electrodes, a set of four woven mesh wire grids. Grid design drivers are discussed, with emphasis placed on design of the electron extraction region. We follow this with a description of the grid production process and a discussion of steps taken to validate the LZ grids prior to integration into the TPC.

Recently, the Thakurta metric has been adopted as a model of primordial black holes by several authors. We show that the spacetime described by this metric has neither black-hole event horizon nor black-hole trapping horizon and involves the violation of the null energy condition as a solution of the Einstein equation. Therefore, this metric does not describe a cosmological black hole in the early universe.

We discuss stiffening of matter in quark-hadron continuity. We introduce a model that relates quark wavefunctions in a baryon and the occupation probability of states for baryons and quarks in dense matter. In dilute regime, the confined quarks contribute to the energy density through the masses of baryons, but do not directly contribute to the pressure, hence the equations of state are very soft. This dilute regime continues until the low momentum states for quarks get saturated; this may happen even before baryons fully overlap, possibly at density slightly above the nuclear saturation density. After the saturation the pressure grows rapidly while changes in energy density are modest, producing a peak in the speed of sound. If we use baryonic descriptions for quark distributions near the Fermi surface, we reach a description similar to the quarkyonic matter model of McLerran and Reddy. With a simple adjustment of quark interactions to get the nucleon mass, our model becomes consistent with the constraints from 1.4-solar mass neutron stars, but the high density part is too soft to account for two-solar mass neutron stars. We delineate the relation between the saturation effects and short range interactions of quarks, suggesting interactions that leave low density equations of state unchanged but stiffen the high density part.

We continue our investigations of the optical properties of the solar gravitational lens (SGL). We treat the Sun as an extended axisymmetric body and model its gravitational field using zonal harmonics. We consider a point source that is positioned at a large but finite distance from the Sun and, using our new angular eikonal method, we established the electro-magnetic (EM) field on the image plane in the focal region behind the SGL and derive the SGL's impulse response in the form of its point-spread function (PSF). The expression that we derive describes the extended Sun in all regions of interest, including the regions of strong and weak interference and the region of geometric optics. The result is in the form of a single integral with respect to the azimuthal angle of the impact parameter, covering all lensing regimes of the SGL. The same expression can be used to describe gravitational lensing by a compact axisymmetric mass distribution, characterized by small deviations from spherical symmetry. It is valid in all lensing regimes. We also derive results that describe the intensity of light observed by an imaging telescope in the focal region. We present results of numerical simulations showing the view by a telescope that moves in the image plane toward the optical axis. We consider imaging of both point and extended sources. We show that while point sources yield a number of distinct images consistent with the caustics due to zonal harmonics of a particular order (e.g., Einstein cross), extended sources always result in forming an Einstein ring. These results represent the most comprehensive wave-theoretical treatment of gravitational lensing in the weak gravitational field of a compact axisymmetric gravitational lens.

A viable quantum theory does not allow curvature invariant terms of different higher orders to be accommodated in the gravitational action. We show that there is indeed a conflict between the curvature squared and Gauss-Bonnet squared term from the point of view of hermiticity. This means one should choose either, in addition to the Einstein-Hilbert term, but never the two together. The choice may be made from inflationary paradigm.

Axion is a promising candidate for ultralight dark matter which may cause a polarization rotation of laser light. Recently, a new idea of probing the axion dark matter by optical linear cavities used in the arms of gravitational wave detectors has been proposed [Phys. Rev. Lett. 123, 111301 (2019)]. In this article, a realistic scheme of the axion dark matter search with the arm cavity transmission ports is revisited. Since photons detected by the transmission ports travel in the cavity for odd-number of times, the effect of axion dark matter on their phases is not cancelled out and the sensitivity at low-mass range is significantly improved compared to the search using reflection ports. We also take into account the stochastic nature of the axion field and the availability of the two detection ports in the gravitational wave detectors. The sensitivity to the axion-photon coupling, $g_{a\gamma}$, of the ground-based gravitational wave detector, such as Advanced LIGO, with 1-year observation is estimated to be $g_{a\gamma} \sim 3\times10^{-12}$ GeV$^{-1}$ below the axion mass of $10^{-15}$ eV, which improves upon the limit achieved by the CERN Axion Solar Telescope.

We present a systematic investigation of the possible locations for the special point (SP), a unique feature of hybrid neutron stars in the mass-radius. The study is performed within the two-phase approach where the high-density (quark matter) phase is described by the constant-sound-speed (CSS) equation of state (EoS) and the nuclear matter phase around saturation density is varied from very soft (APR) to stiff (DD2 with excluded nucleon volume. Different construction schemes for the deconfinement transition are applied: Maxwell construction, mixed phase construction and parabolic interpolation. We demonstrate for the first time that the SP is invariant not only against changing the nuclear matter EoS, but also against variation of the construction schemes for the phase transition. Since the SP serves as a proxy for the maximum mass and accessible radii of massive hybrid stars, we draw conclusions for the limiting masses and radii of hybrid neutron stars.

We discuss the energy loss due to gravitational radiation of binaries composed of exotic objects whose horizon boundary conditions are replaced with reflective ones. Our focus is on the extreme mass-ratio inspirals, in which the central heavier black hole is replaced with an exotic compact object. We show, in this case, a modulation of the energy loss rate depending on the evolving orbital frequency occurs and leads to two different types of modifications to the gravitational wave phase evolution; the oscillating part directly corresponding to the modulation in the energy flux, and the non-oscillating part coming from the quadratic order in the modulation. This modification can be sufficiently large to detect with future space-borne detectors.

Extreme mass-ratio inspirals (EMRIs) detectable by the Laser Interferometer Space Antenna (LISA) are unique probes of the nature of supermassive compact objects. We compute the gravitational-wave signal emitted by a stellar-mass compact object in a circular equatorial orbit around a Kerr-like horizonless supermassive object defined by an effective radius and a reflectivity coefficient. The Teukolsky equations are solved consistently with suitable (frequency-dependent) boundary conditions, and the modified energy and angular-momentum fluxes are used to evolve the orbital parameters adiabatically. The gravitational fluxes have resonances corresponding to the low-frequency quasinormal modes of the central object, which can contribute significantly to the gravitational-wave phase. Overall, the absence of a classical event horizon in the central object can affect the gravitational-wave signal dramatically, with deviations even larger than those previously estimated by a model-independent analysis of the tidal heating. We estimate that EMRIs could potentially place the most stringent constraint on the reflectivity of supermassive compact objects at the remarkable level of ${\cal O}(10^{-6})\%$ and would allow to constrain various models which are not ruled out by the ergoregion instability. In particular, an EMRI detection could allow to rule out (or provide evidence for) signatures of quantum black-hole horizons with Boltzmann reflectivity. Our results motivate performing rigorous parameter estimations to assess the detectability of these effects.

The COSINUS (Cryogenic Observatory for SIgnatures seen in Next-generation Underground Searches) experiment aims at the detection of dark matter-induced recoils in sodium iodide (NaI) crystals operated as scintillating cryogenic calorimeters. The detection of both scintillation light and phonons allows performing an event-by-event signal to background discrimination, thus enhancing the sensitivity of the experiment. The construction of the experimental facility is foreseen to start by 2021 at the INFN Gran Sasso National Laboratory (LNGS) in Italy. It consists of a cryostat housing the target crystals shielded from the external radioactivity by a water tank acting, at the same time, as an active veto against cosmic ray-induced events. Taking into account both environmental radioactivity and intrinsic contamination of materials used for cryostat, shielding and infrastructure, we performed a careful background budget estimation. The goal is to evaluate the number of events that could mimic or interfere with signal detection while optimising the geometry of the experimental setup. In this paper we present the results of the detailed Monte Carlo simulations we performed, together with the final design of the setup that minimises the residual amount of background particles reaching the detector volume.

Lovelock gravity in $D$-dimensional space-times is considered adopting Cartan's structure equations. In this context, we find out exact solutions in cosmological and spherically symmetric backgrounds. In the latter case, we also derive horizons and the corresponding Bekenstein-Hawking entropies. Moreover, we focus on the topological Chern-Simons theory, providing exact solutions in five dimensions. Specifically, it is possible to show that Anti-de Sitter invariant Chern-Simons gravity can be framed within Lovelock-Zumino gravity in $5$-dimensions, for particular choices of Lovelock parameters.

We investigate the gravitational wave spectrum resulted from the cosmological first-order phase transition. We compare two models; one is a scalar field model without gravitation, while the other is a scalar field model with gravitation. Based on the sensitivity curves of the eLISA space-based interferometer on the stochastic gravitational-wave background, we compare the difference between the gravitational wave spectra of the former and the latter cases resulted from the bubble collision process. Especially, we calculated the speed of the bubble wall before collision for the two models numerically. We show that the difference between the amplitudes of those spectra can clearly distinguish between the two models. We expect that the eLISA with C4 could observe the spectrum as the fast first-order phase transition or that with C1 as the slow first-order phase transition from the latter case.

For general relativistic, inviscid, axisymmetric flow around Kerr black hole one may choose different flow thickness. The stationary flow equations can be solved using methods of dynamical system to get transonic accretion flows , i.e, flow infalling in the blackhole that turns supersonic from subsonic with decreasing radial distance, or vice versa. This transonic flows are obtained by choosing the particular flow passing through critical points of phase portrait. For certain flow thickness like the one maintaining conical shape, the sonic point coincide with the critical point. But there are certain flows maintaining hydrostatic equilibrium, such as the one described by Novikov-Thorne, where the sonic point is not same as the critical point. We perturb the flow for both kind of flow and study the behaviour of linear perturbation which behaves like massless scalar field in some curved spacetime, known as, analogue space time. We draw the compactified causal structure, i.e, Penrose Carter diagram for both kind of analogue metric and prove that for both cases critical points are the acoustic horizons, whereas in the case where sonic points do not coincide with critical points, the sonic points are not the acoustic horizon, as one may expect from the definition of sound speed.

Analogue gravity models describe linear fluctuations of fluids as a massless scalar field propagating on stationary acoustic spacetimes constructed from the background flow. In this paper, we establish that this paradigm generalizes to arbitrary order nonlinear perturbations propagating on dynamical analogue spacetimes. Our results hold for all inviscid, spherically symmetric and barotropic non-relativistic flows in the presence of an external conservative force. We demonstrate that such fluids always admit a dynamical description governed by a coupled pair of wave and continuity equations. We provide an iterative approach to solve these equations about any known stationary solution to all orders in perturbation. In the process, we reveal that there exists a dynamical acoustic spacetime on which fluctuations of the mass accretion rate propagate. The dynamical acoustic spacetime is shown to have a well defined causal structure and curvature. In addition, we find a classical fluctuation relation for the acoustic horizon of the spacetime that admit scenarios wherein the horizon can grow as well as recede, with the latter being a result with no known analogue in black holes. As an example, we numerically investigate the Bondi flow accreting solution subject to exponentially damped time dependent perturbations. We find that second and higher order classical perturbations possess an acoustic horizon that oscillates and changes to a new stable size at late times. In particular, the case of a receding acoustic horizon is realized through `low frequency' perturbations. We discuss our results in the context of more general analogue models and its potential implications on astrophysical accretion flows.

We compare two competing relativistic approaches to the N-body simulation of the Universe large-scale structure. To this end, employing the corresponding alternative computer codes ("gevolution" and "screening"), we conduct a series of cosmological simulations in boxes of different sizes and calculate the power spectra of the scalar perturbation $\Phi$, the frame-dragging vector potential ${\bf B}$ and the difference between scalar modes $\chi=\Phi-\Psi$. We demonstrate that the corresponding power spectra are in very good agreement between the compared schemes. For example, the relative difference of the power spectra for $\Phi$ is 0.04% maximum. Since the perturbed Einstein equations have much simpler form in the screening approach, the simulation with this code consumes less computational time, saving almost 40% of CPU hours.