Papers for Tuesday, Dec 14 2021

Making the most of the rapidly increasing population of gravitational-wave detections of black hole (BH) and neutron star (NS) mergers requires comparing observations with population synthesis predictions. In this work we investigate the combined impact from the key uncertainties in population synthesis modelling of the isolated binary evolution channel: the physical processes in massive binary-star evolution and the star formation history as a function of metallicity, $Z$, and redshift $z, \mathcal{S}(Z,z)$. Considering these uncertainties we create 560 different publicly available model realizations and calculate the rate and distribution characteristics of detectable BHBH, BHNS, and NSNS mergers. We find that our stellar evolution and $\mathcal{S}(Z,z)$ variations can impact the predicted intrinsic and detectable merger rates by factors $10^2$-$10^4$. We find that BHBH rates are dominantly impacted by $\mathcal{S}(Z,z)$ variations, NSNS rates by stellar evolution variations and BHNS rates by both. We then consider the combined impact from all uncertainties considered in this work on the detectable mass distribution shapes (chirp mass, individual masses and mass ratio). We find that the BHNS mass distributions are predominantly impacted by massive binary-star evolution changes. For BHBH and NSNS we find that both uncertainties are important. We also find that the shape of the delay time and birth metallicity distributions are typically dominated by the choice of $\mathcal{S}(Z,z)$ for BHBH, BHNS and NSNS. We identify several examples of robust features in the mass distributions predicted by all 560 models, such that we expect more than 95% of BHBH detections to contain a BH $\gtrsim 8\,\rm{M}_{\odot}$ and have mass ratios $\lesssim 4$. Our work demonstrates that it is essential to consider a wide range of allowed models to study double compact object merger rates and properties.

Precise reconstruction of the cosmic microwave background lensing potential can be achieved with deep polarization surveys by iteratively removing lensing-induced $B$ modes. We introduce a lensing spectrum estimator and its likelihood for such optimal iterative reconstruction. Our modelling share similarities to the state-of-the-art likelihoods for quadratic estimator-based (QE) lensing reconstruction. In particular, we generalize the $N^{(0)}$ and $N^{(1)}$ lensing biases, and design a realization-dependent spectrum debiaser, making this estimator robust to uncertainties in the data modelling. We demonstrate unbiased recovery of the cosmology using map-based reconstructions. For a CMB-S4 survey, this spectrum estimator and likelihood can double the constraints on the lensing amplitude compared to the QE on a wide range of scales, while keeping numerical cost under control and being robust to errors.

We compare observations of HI from the Very Large Array (VLA) and the Arecibo Observatory and observations of HCO$^+$ from the Atacama Large Millimeter/submillimeter Array (ALMA) and the Northern Extended Millimeter Array (NOEMA) in the diffuse ($A_V\lesssim1$) interstellar medium (ISM) to predictions from a photodissociation region (PDR) chemical model and multi-phase ISM simulations. Using a coarse grid of PDR models, we estimate the density, FUV radiation field, and cosmic ray ionization rate (CRIR) for each structure identified in HCO$^+$ and HI absorption. These structures fall into two categories. Structures with $T_s<40~\mathrm{K}$, mostly with $N(\mathrm{HCO^+})\lesssim10^{12}~\mathrm{cm^{-2}}$, are consistent with modest density, FUV radiation field, and CRIR models, typical of the diffuse molecular ISM. Structures with spin temperature $T_s>40~\mathrm{K}$, mostly with $N(\mathrm{HCO^+})\gtrsim10^{12}~\mathrm{cm^{-2}}$, are consistent with high density, FUV radiation field, and CRIR models, characteristic of environments close to massive star formation. The latter are also found in directions with a significant fraction of thermally unstable HI. In at least one case, we rule out the PDR model parameters, suggesting that alternative mechanisms (e.g., non-equilibrium processes like turbulent dissipation and/or shocks) are required to explain the observed HCO$^+$ in this direction. Similarly, while our observations and simulations of the turbulent, multi-phase ISM agree that HCO$^+$ formation occurs along sightlines with $N(\mathrm{HI})\gtrsim10^{21}~\mathrm{cm^{-2}}$, the simulated data fail to explain HCO$^+$ column densities $\gtrsim\rm{few}\times10^{12}~\mathrm{cm^{-2}}$. Since a majority of our sightlines with HCO$^+$ had such high column densities, this likely indicates that non-equilibrium chemistry is important for these lines of sight.

We present a deep machine learning (ML) approach to constraining cosmological parameters with multi-wavelength observations of galaxy clusters. The ML approach has two components: an encoder that builds a compressed representation of each galaxy cluster and a flexible CNN to estimate the cosmological model from a cluster sample. It is trained and tested on simulated cluster catalogs built from the Magneticum simulations. From the simulated catalogs, the ML method estimates the amplitude of matter fluctuations, sigma_8, at approximately the expected theoretical limit. More importantly, the deep ML approach can be interpreted. We lay out three schemes for interpreting the ML technique: a leave-one-out method for assessing cluster importance, an average saliency for evaluating feature importance, and correlations in the terse layer for understanding whether an ML technique can be safely applied to observational data. These interpretation schemes led to the discovery of a previously unknown self-calibration mode for flux- and volume-limited cluster surveys. We describe this new mode, which uses the amplitude and peak of the cluster mass PDF as anchors for mass calibration. We introduce the term "overspecialized" to describe a common pitfall in astronomical applications of machine learning in which the ML method learns simulation-specific details, and we show how a carefully constructed architecture can be used to check for this source of systematic error.

Planetesimals inevitably bear the signatures of their natal environment, preserving in their composition a record of the metallicity of their system's original gas and dust, albeit one altered by the formation process. When planetesimals are dispersed from their system of origin, this record is carried with them. As each star is likely to contribute at least $10^{12}$ interstellar objects, the Galaxy's drifting population of interstellar objects (ISOs) provides an overview of the properties of its stellar population through time. Using the EAGLE cosmological simulation and models of protoplanetary formation, our modelling predicts an ISO population with a bimodal distribution in their water mass fraction. Objects formed in low-metallicity, typically older, systems have a higher water fraction than their counterparts formed in high-metallicity protoplanetary disks, and these water-rich objects comprise the majority of the population. Both detected ISOs seem to belong to the lower water fraction population; these results suggest they come from recently formed systems. We show that the population of ISOs in galaxies with different star formation histories will have different proportions of objects with high and low water fractions. This work suggests that it is possible that the upcoming Vera C. Rubin Observatory Legacy Survey of Space and Time will detect a large enough population of ISOs to place useful constraints on models of protoplanetary disks, as well as galactic structure and evolution.

The strong intervening absorption system at redshift 1.15 towards the very bright quasar HE 0515$-$4414 is the most studied absorber for measuring possible cosmological variations in the fine-structure constant, $\alpha$. We observed HE 0515$-$4414 for 16.1$\,$h with the Very Large Telescope and present here the first constraint on relative variations in $\alpha$ with parts-per-million (ppm) precision from the new ESPRESSO spectrograph: $\Delta\alpha/\alpha = 1.3 \pm 1.3_{\rm stat} \pm 0.4_{\rm sys}\,{\rm ppm}$. The statistical uncertainty (1$\sigma$) is similar to the ensemble precision of previous large samples of absorbers, and derives from the high S/N achieved ($\approx$105 per 0.4$\,$km$\,$s$^{-1}$ pixel). ESPRESSO's design, and calibration of our observations with its laser frequency comb, effectively removed wavelength calibration errors from our measurement. The high resolving power of our ESPRESSO spectrum ($R=145000$) enabled the identification of very narrow components within the absorption profile, allowing a more robust analysis of $\Delta\alpha/\alpha$. The evidence for the narrow components is corroborated by their correspondence with previously detected molecular hydrogen and neutral carbon. The main remaining systematic errors arise from ambiguities in the absorption profile modelling, effects from redispersing the individual quasar exposures, and convergence of the parameter estimation algorithm. All analyses of the spectrum, including systematic error estimates, were initially blinded to avoid human biases. We make our reduced ESPRESSO spectrum of HE 0515$-$4414 publicly available for further analysis. Combining our ESPRESSO result with 28 measurements, from other spectrographs, in which wavelength calibration errors have been mitigated, yields a weighted mean $\Delta\alpha/\alpha = -0.5 \pm 0.5_{\rm stat} \pm 0.4_{\rm sys}\,$ppm at redshifts 0.6-2.4.

The self-lensing of a massive black hole binary (MBHB), which occurs when the two BHs are aligned close to the line of sight, is expected to produce periodic, short-duration flares. Here we study the shapes of self-lensing flares (SLFs) via general-relativistic ray tracing in a superimposed binary BH metric, in which the emission is generated by geometrically thin accretion flows around each component. The suite of models covers eccentric binary orbits, black hole spins, unequal mass binaries, and different emission model geometries. We explore the above parameter space, and report how the light curves change as a function of, e.g., binary separation, inclination, and eccentricity. We also compare our light curves to those in the microlensing approximation, and show how strong deflections, as well as time-delay effects, change the size and shape of the SLF. If gravitational waves (GWs) from the inspiraling MBHB are observed by LISA, SLFs can help securely identify the source and localizing it on the sky, and to constrain the graviton mass by comparing the phasing of the SLFs and the GWs. Additionally, when these systems are viewed edge-on the SLF shows a distinct dip that can be directly correlated with the BH shadow size. This opens a new way to measure BH shadow sizes in systems that are unresolvable by current VLBI facilities.

Supermassive black hole (BH) binaries are thought to produce self-lensing flares (SLF) when the two BHs are aligned with the line-of-sight. If the binary orbit is observed nearly edge-on, we find a distinct feature in the light curve imprinted by the BH shadow from the lensed BH. We study this feature by ray-tracing in a binary model and predict that 1\% of the current binary candidates could show this feature. Our BH tomography method proposed here could make it possible to extract BH shadows that are spatially unresolvable by high-resolution VLBI.

We quantitatively address the conjecture that magnetic helicity must be shed from the Sun by eruptions launching coronal mass ejections in order to limit its accumulation in each hemisphere. By varying the ratio of guide and strapping field and the flux rope twist in a parametric simulation study of flux rope ejection from approximately marginally stable force-free equilibria, different ratios of self- and mutual helicity are set and the onset of the torus or helical kink instability is obtained. The helicity shed is found to vary over a broad range from a minor to a major part of the initial helicity, with self helicity being largely or completely shed and mutual helicity, which makes up the larger part of the initial helicity, being shed only partly. Torus-unstable configurations with subcritical twist and without a guide field shed up to about two-thirds of the initial helicity, while a highly twisted, kink-unstable configuration sheds only about one-quarter. The parametric study also yields stable force-free flux rope equilibria up to a total flux-normalized helicity of 0.25, with a ratio of self- to total helicity of 0.32 and a ratio of flux rope to external poloidal flux of 0.94. These results numerically demonstrate the conjecture of helicity shedding by coronal mass ejections and provide a first account of its parametric dependence. Both self- and mutual helicity are shed significantly; this reduces the total initial helicity by a fraction of $\sim\!0.4\mbox{--}0.65$ for typical source region parameters.

We present the results of our detailed light curve analysis of the ZZ Ceti star HS 1625+1231. We collected photometric time series data at Konkoly Observatory on 14 nights, and performed Fourier analysis of these data sets. We detected 11 significant frequencies, where six of them are found to be independent pulsation modes in the 514 - 881 s period range. By utilising these frequencies, we performed preliminary asteroseismic investigations to give constraints on the main physical parameters, and to derive seismic distances for the star. Finally, we compared the astrometric distance provided by the Gaia EDR3 data with those seismic distances. Our selected model, considering both the spectroscopic measurements and the distance value provided by Gaia, has $T_{eff}$ = 11 000 K and $M_*$ = 0.60 $M_{\odot}$.

Using 250 neutron star merger simulations with microphysics, we explore for the first time the role of nuclear incompressibility in the prompt collapse threshold for binaries with different mass ratios. We demonstrate that observations of prompt collapse thresholds, either from binaries with two different mass ratios or with one mass ratio but combined with the knowledge of the maximum neutron star mass or compactness, will constrain the incompressibility at the maximum neutron star density, $K_{\rm max}$, to within tens of percent. This, otherwise inaccessible, measure of $K_{\rm max}$ can potentially reveal the presence of hyperons or quarks inside neutron stars.

Spherically symmetric and static dark matter halos in hydrostatic equilibrium demand that dark matter should have an effective pressure that compensates the gravitational force of the mass of the halo. An effective equation of state can be obtained for each rotational velocity profile of the stars in galaxies. In this work, we study one of this dark matter equation of state obtained for the Universal Velocity Profile and analyze the properties of the self-gravitating structures that emerges from this equation of state. The resulting configurations explaining the observed rotational speeds are found to be unstable. We conclude that either the halo is not in hydrostatic equilibrium, or it is non spherically symmetric, or it is not static if the Universal Velocity profile should be valid to fit the rotational velocity curve of the galaxies.

We re-analyze the pristine eclipsing binary data from the $\textit{Kepler}$ and TESS missions, focusing on eccentricity measurements at short orbital periods to emperically constrain tidal circularization. We find an average circularization period of $~$6 days, as well as a short circularization period of $\sim$3 days for the $\textit{Kepler}$/TESS field binaries. We argue previous spectroscopic binary surveys reported longer circularization periods due to small sample sizes, which were contaminated by an abundance of binaries with circular orbits out to $\sim$10 days, but we re-affirm their data shows a difference between the eccentricity distributions of young ($<$1 Gyr) and old ($>$3 Gyr) binaries. Our work calls into question the long circularization periods quoted often in the literature.

We make use of neural networks to accelerate the calculation of power spectra required for the analysis of galaxy clustering and weak gravitational lensing data. For modern perturbation theory codes, evaluation time for a single cosmology and redshift can take on the order of two seconds. In combination with the comparable time required to compute linear predictions using a Boltzmann solver, these calculations are the bottleneck for many contemporary large-scale structure analyses. In this work, we construct neural network-based surrogate models for Lagrangian perturbation theory (LPT) predictions of matter power spectra, real and redshift space galaxy power spectra, and galaxy--matter cross power spectra that attain $\sim 0.1\%$ (at one sigma) accuracy over a broad range of scales in a $w$CDM parameter space. The neural network surrogates can be evaluated in approximately one millisecond, a factor of 1000 times faster than the full Boltzmann code and LPT computations. In a simulated full-shape redshift space galaxy power spectrum analysis, we demonstrate that the posteriors obtained using our surrogates are accurate compared to those obtained using the full LPT model. We make our surrogate models public at https://github.com/sfschen/EmulateLSS, so that others may take advantage of the speed gains they provide to enable rapid iteration on analysis settings, something that is essential in complex contemporary large-scale structure analyses.

Astronomical surveys continue to provide unprecedented insights into the time-variable Universe and will remain the source of groundbreaking discoveries for years to come. However, their data throughput has overwhelmed the ability to manually synthesize alerts for devising and coordinating necessary follow-up with limited resources. The advent of Rubin Observatory, with alert volumes an order of magnitude higher at otherwise sparse cadence, presents an urgent need to overhaul existing human-centered protocols in favor of machine-directed infrastructure for conducting science inference and optimally planning expensive follow-up observations. We present the first implementation of autonomous real-time science-driven follow-up using value iteration to perform sequential experiment design. We demonstrate it for strategizing photometric augmentation of Zwicky Transient Facility Type Ia supernova light-curves given the goal of minimizing SALT2 parameter uncertainties. We find a median improvement of 2-6% for SALT2 parameters and 3-11% for photometric redshift with 2-7 additional data points in g, r and/or i compared to random augmentation. The augmentations are automatically strategized to complete gaps and for resolving phases with high constraining power (e.g. around peaks). We suggest that such a technique can deliver higher impact during the era of Rubin Observatory for precision cosmology at high redshift and can serve as the foundation for the development of general-purpose resource allocation systems.

NASA's Juno mission recently reported Jupiter's high-degree (degree $\ell$, azimuthal order $m$ $=4,2$) Love number $k_{42}=1.289\pm0.063$ ($1\sigma$), an order of magnitude above the hydrostatic $k_{42}$ obtained in a nonrotating Jupiter model. After numerically modeling rotation, the hydrostatic $k_{42}=1.743\pm0.002$ is still $7\sigma$ away from the observation, raising doubts about our understanding of Jupiter's tidal response. Here, we use first-order perturbation theory to explain the hydrostatic $k_{42}$ result analytically. We use a simple Jupiter equation of state ($n=1$ polytrope) to obtain the fractional change in $k_{42}$ when comparing a rotating model with a nonrotating model. Our analytical result shows that the hydrostatic $k_{42}$ is dominated by the tidal response at $\ell=m=2$ coupled into the spherical harmonic $\ell,m=4,2$ by the planet's oblate figure. The $\ell=4$ normalization in $k_{42}$ introduces an orbital factor $(a/s)^2$ into $k_{42}$, where $a$ is the satellite semimajor axis and $s$ is Jupiter's average radius. As a result, different Galilean satellites produce a different $k_{42}$. We conclude that high-degree tesseral Love numbers ($\ell> m$, $m\geq2$) are dominated by lower-degree Love numbers and thus provide little additional information about interior structure, at least when they are primarily hydrostatic. Our results entail important implications for a future interpretation of the currently observed Juno $k_{42}$. After including the coupling from the well-understood $\ell=2$ dynamical tides ($\Delta k_2 \approx -4\%$), Jupiter's hydrostatic $k_{42}$ requires an unknown dynamical effect to produce a fractional correction $\Delta k_{42}\approx-11\%$ in order to fit Juno's observation within $3\sigma$. Future work is required to explain the required $\Delta k_{42}$.

The distinctive set of Saturnian small satellites, Aegaeon, Methone, Anthe, and Pallene, constitutes an excellent laboratory to understand the evolution of systems immersed in co-orbital dusty rings/arcs, subjected to perturbations from larger satellites and non-gravitational forces. In this work, we carried out a comprehensive numerical exploration of the long-term evolution of Pallene and its ring. Through frequency map analysis, we characterised the current dynamical state around Pallene. A simple tidal evolution model serves to set a time frame for the current orbital configuration of the system. With detailed short and long-term N-body simulations we determine whether Pallene is currently in resonance with one or more of six of Saturn's major moons. We analysed a myriad of resonant arguments extracted from the direct and indirect parts of the disturbing function, finding that Pallene is not in mean motion resonance from the present up to 5~Myr into the future; nonetheless, some resonant arguments exhibit intervals of libration and circulation at different timescales and moon pairings. We studied the dynamical evolution of micrometric particles forming the ring, considering gravitational and non-gravitational forces. Non-gravitational forces are responsible for particles vertical excursions and outward migration. By estimating the satellite's mass production rate, we find that Pallene could be responsible for keeping its ring in steady-state only if it is mainly composed of large micrometre-sized particles. If mainly composed of particles with a few micrometres for which Pallene is the only source, the ring will spread out, both radially and vertically, until it finally disappears.

Observational pre-cursors of large solar flares provide a basis for future operational systems for forecasting. Here, we study the evolution of the normalized emergence (EM), shearing (SH) and total (T) magnetic helicity flux components for 14 flaring with at least one X-class flare) and 14 non-flaring ($<$ M5-class flares) active regions (ARs) using the Spaceweather Helioseismic Magnetic Imager Active Region Patches vector magnetic field data. Each of the selected ARs contain a $\delta$-type spot. The three helicity components of these ARs were analyzed using wavelet analysis. Localised peaks of the wavelet power spectrum (WPS) were identified and statistically investigated. We find that: i) the probability density function of the identified WPS peaks for all the EM, SH and T profiles can be fitted with a set of Gaussian functions centered at distinct periods between $\sim$ 3 to 20 hours. ii) There is a noticeable difference in the distribution of periods found in the EM profiles between the flaring and non-flaring ARs, while no significant difference is found in the SH and T profiles. iii) In flaring ARs, the distributions of the shorter EM/SH/T periods ($<$ 10 hrs) split up into two groups after flares, while the longer periods ($>$ 10 hrs) do not change. iv) When the EM periodicity does not contain harmonics, the ARs do not host a large energetic flare. Finally, v) significant power at long periods ($\sim$ 20 hour) in the T and EM components may serve as pre-cursor for large energetic flares.

In this chapter, we discuss the contributions of gamma-ray astronomy at TeV energies to our understanding of the visible content and structure of the universe. We start from the present epoch with the second most intense electromagnetic background field after the CMB: the extragalactic background light (EBL). The EBL is composed of all the light emitted by stars and galaxies since the beginning of reionization, including light absorbed and re-emitted by dust. As such, the EBL traces the history of radiating matter in the universe. We then further dive into the large voids of the universe to study the large-scale magnetic fields that should permeate them. These fields could originate from the onset of structure formation or early phase transitions, bringing us back to the infancy of the universe. We conclude by looking back to the elusive Planck time scale, where the standard models of cosmology and particle physics are no longer applicable. Observations with current-generation gamma-ray astronomy experiments have now started to scratch the surface of cosmology, as we will show in this chapter.

We introduce a novel compact 4-channel beam splitter which is based on a combination of dichroic coatings and internal total reflection, similar in concept to the interference double-prism invented by K\"osters 90 years ago. Used with a rapidly-slewing 50 cm telescope in space, this would allow to double the presently known gamma-ray bursts at high (>5) redshift within 2 years.

Luminosity bursts in young FU Orionis-type stars warm up the surrounding disks of gas and dust, thus inflicting changes on their morphological and chemical composition. In this work, we aim at studying the effects that such bursts may have on the spatial distribution of dust grain sizes and the corresponding spectral index in protoplanetary disks. We use the numerical hydrodynamics code FEOSAD, which simulates the co-evolution of gas, dust, and volatiles in a protoplanetary disk, taking dust growth and back reaction on gas into account. The dependence of the maximum dust size on the water ice mantles is explicitly considered. The burst is initialized by increasing the luminosity of the central star to 100-300 L_sun for a time period of 100 yr. The water snowline shifts during the burst to a larger distance, resulting in the drop of the maximum dust size interior to the snowline position because of more efficient fragmentation of bare grains. After the burst, the water snowline shifts quickly back to its preburst location followed by renewed dust growth. The timescale of dust regrowth after the burst depends on the radial distance so that the dust grains at smaller distances reach the preburst values faster than the dust grains at larger distances. As a result, a broad peak in the radial distribution of the spectral index in the millimeter dust emission develops at \approx 10 au, which shifts further out as the disk evolves and dust grains regrow to preburst values at progressively larger distances. This feature is most pronounced in evolved axisymmetric disks rather than in young gravitationally unstable counterparts, although young disks may still be good candidates if gravitational instability is suppressed. Abridged.

We present a RR Lyrae (RRL) catalogue based on the combination of the third data release of the Zwicky Transient Facility (ZTF DR3) and \textit{Gaia} EDR3. We use a multi-step classification pipeline relying on the Fourier decomposition fitting to the multi-band ZTF light curves and random forest classification. The resulting catalogue contains 71,755 RRLs with period and light curve parameter measurements and has completeness of 0.92 and purity of 0.92 with respect to the SOS \textit{Gaia} DR2 RRLs. The catalogue covers the Northern sky with declination $\geq -28^\circ$, its completeness is $\gtrsim 0.8$ for heliocentric distance $\leq 80$~kpc, and the most distant RRL at 132~kpc. Compared with several other RRL catalogues covering the Northern sky, our catalogue has more RRLs around the Galactic halo and is more complete at low Galactic latitude areas. Analysing the spatial distribution of RRL in the catalogue reveals the previously known major over-densities of the Galactic halo, such as the Virgo over-density and the Hercules-Aquila Cloud, with some evidence of an association between the two. We also analyse the Oosterhoff fraction differences throughout the halo, comparing it with the density distribution, finding increasing Oosterhoff I fraction at the elliptical radii between 16 and 32 kpc and some evidence of different Oosterhoff fractions across various halo substructures.

The collision history of asteroids is an important archive of inner Solar System evolution. Evidence for these collisions is brought to Earth by meteorites, which can preserve impact-reset radioisotope mineral ages. However, as meteorites often preserve numerous mineral ages, their interpretation is controversial. Here, we combine analysis of phosphate U-Pb ages and mineral microtextures to construct a collision history for the highly shocked Chelyabinsk meteorite. We show that phosphate U-Pb ages in the meteorite are independent of thermal history at macro-to-microscales, correlating instead with phosphate microtexture. Isotopic data from pristine phosphate domains is largely concordant, whereas fracture-damaged domains universally display discordance. Combining both populations best constrains upper (4,473 +/- 11 Ma) and lower intercept (-9 +/- 55 Ma, i.e., within error of the present day) U-Pb ages for Chelyabinsk phosphates. We conclude that all phosphate U-Pb ages were completely reset during an ancient high energy collision. Fracture-damaged phosphate domains experienced further Pb-loss during mild collisional heating in the geologically recent past, and must be targeted to properly constrain a lower intercept age. Targeting textural sub-populations of phosphate grains can significantly improve the calculation and interpretation of U-Pb ages, permitting more robust reconstruction of both ancient and recent asteroidal collision histories.

We provide stability estimates, obtained by implementing the Nekhoroshev theorem, in reference to the orbital motion of a small body (satellite or space debris) around the Earth. We consider a Hamiltonian model, averaged over fast angles, including the $J_2$ geopotential term as well as third-body perturbations due to Sun and Moon. We discuss how to bring the Hamiltonian into a form suitable for the implementation of the Nekhoroshev theorem in the version given by P\"oschel(1993) for the `non-resonant' regime. The manipulation of the Hamiltonian includes i) averaging over fast angles, ii) a suitable expansion around reference values for the orbit's eccentricity and inclination, and iii) a preliminary normalization allowing to eliminate particular terms whose existence is due to the non-zero inclination of the invariant plane of secular motions known as the `Laplace plane'. After bringing the Hamiltonian to a suitable form, we examine the domain of applicability of the theorem in the action space, translating the result in the space of physical elements. We find that the necessary conditions for the theorem to hold are fulfilled in some non-zero measure domains in the eccentricity and inclination plane (e, i) for a body's orbital altitude (semi-major axis) up to about 20000 km. For altitudes around 11000 km we obtain stability times of the order of several thousands of years in domains covering nearly all eccentricities and inclinations of interest in applications of the satellite problem, except for narrow zones around some so-called `inclination-dependent' resonances. On the other hand, the domains of Nekhoroshev stability recovered by the present method shrink in size as the semi-major axis a increases (and the corresponding Nekhoroshev times reduce to hundreds of years), while the stability domains practically all vanish for a > 20000 km.

Data acquired by the Virgo interferometer during the second part of the O3 scientific run, referred to as O3b, were analysed with the aim of characterising the onset and time evolution of scattered light noise in connection with the variability of microseismic noise in the environment surrounding the detector. The adaptive algorithm used, called pytvfemd, is suitable for the analysis of time series which are both nonlinear and nonstationary. It allowed to obtain the first oscillatory mode of the differential arm motion degree of freedom of the detector during days affected by scattered light noise. The mode's envelope i.e., its instantaneous amplitude, is then correlated with the motion of the West end bench, a known source of scattered light during O3. The relative velocity between the West end test mass and the West end optical bench is used as a predictor of scattered light noise. Higher values of correlation are obtained in periods of higher seismic noise in the microseismic frequency band. This is also confirmed by the signal-to-noise ratio (SNR) of scattered light glitches from GravitySpy for the January-March 2020 period. Obtained results suggest that the adopted methodology is suited for scattered light noise characterisation and monitoring in gravitational wave interferometers.

The paper discusses a model of the bombardment of the Earth and the Moon by small bodies when these planets were formed. It is shown that the total ice mass delivered with the bodies to the Earth from the feeding zone of the giant planets and the outer asteroid belt could have been comparable to the total mass of the Earth's oceans. Objects that initially crossed Jupiter's orbit could become Earth-crossers mainly within the first one million years. Most collisions of bodies originally located at a distance of 4 to 5 AU (astronomical units) from the Sun with the Earth occurred during the first ten million years. Some bodies from the Uranus and Neptune zones could fall onto the Earth in more than 20 million years. From their initial distances from the Sun of about 3 to 3.5 AU, some bodies could fall onto the Earth and Moon in a few billion years for the model that takes into account only the gravitational influence of the planets. The ratio of the number of bodies that collided with the Earth to the number of bodies that collided with the Moon varied mainly from 20 to 40 for planetesimals from the feeding zone of the terrestrial planets. For bodies originally located at a distance of more than 3 AU from the Sun, this ratio was mainly in the range between 16.4 and 17.4. The characteristic velocities of collisions of planetesimals from the feeding zones of the terrestrial planets with the Moon varied from 8 to 16 km/s, depending on the initial values of the semi-major axes and eccentricities of orbits of the planetesimals. The collision velocities of bodies that came from the feeding zones of Jupiter and Saturn with the Moon were mainly from 20 to 23 km/s.

Betelgeuse is the nearest red supergiant, one of the brightest stars in our sky, and statistically speaking it would be expected to be "typical". Yet it exhibits many features that seem "curious", to say the least. For instance it has a high proper motion. It rotates fast. It has little dust. It dimmed unexpectedly. Is any of these, and other, phenomena atypical, and taken together does it make Betelgeuse atypical? This is important to know, because we need to know whether Betelgeuse might be a prototype of red supergiants in general, or certain subclasses of red supergiants, since we can study it in such great detail. It is also important to know as it may be a link to understanding other, apparently atypical cases such as supernova 1987A, and maybe even such exotica as Thorne-\.Zytkov objects. Studying this question in itself helps us understand how we deal with rarity and coincidence in understanding the Universe we live in.

Continuous wavelet analysis has been increasingly employed in various fields of science and engineering due to its remarkable ability to maintain optimal resolution in both space and scale. Here, we extend wavelet-based statistics, including the wavelet power spectrum, wavelet cross-correlation and wavelet bicoherence, to the new designed continuous wavelet function -- {\em Gaussian-derived wavelet}. In this paper, these statistics are introduced to analyze the large-scale clustering of matter. For this purpose, we perform wavelet transforms on the density distribution obtained from the one-dimensional (1D) Zel'dovich approximation and then measure the wavelet power spectra and wavelet bicoherences of this density distribution. Our results suggest that the wavelet power spectrum and wavelet bicoherence can identify the effects of local environments on the clustering at different scales. Moreover, to reveal the usefulness of the wavelet cross-correlation, we apply it to the 1D projected density fields of the IllustrisTNG simulation at $z=0$ for different matter components. We find that wavelet cross-correlations between different matter components converge to one on large scales, while biases between them become significant on small scales. In addition, measurements of the wavelet power spectra show that clustering of the total matter is suppressed on scales $k\gtrsim 1 h\mathrm{Mpc}^{-1}$ relative to that of the corresponding dark matter-only simulation. The wavelet bicoherence of the total matter is enhanced on wide scales due to baryonic physics. These results are qualitatively consistent with those from three-dimensional Fourier analyses.

We explored simultaneous observations of chromospheric evaporation and condensation during the impulsive phase of a C6.7 flare on 9 May 2019. The solar flare was simultaneously observed by multiple instruments, i.e., the New Vacuum Solar Telescope (NVST), the Interface Region Imaging Spectrograph, the Atmospheric Imaging Assembly (AIA), the Fermi, the Mingantu Spectral Radioheliograph, and the Nobeyama Radio Polarimeters. Using the single Gaussian fitting and the moment analysis technique, redshifted velocities at slow speeds of 15-19 km/s are found in the cool lines of C II and Si IV at one flare footpoint location. Red shifts are also seen in the H-alpha line-of-sight (LOS) velocity image measured by the NVST at double footpoints. Those red shifts with slow speeds can be regarded as the low-velocity downflows driven by the chromospheric condensation. Meanwhile, the converging motions from double footpoints to the loop top are found in the high-temperature EUV images, such as AIA 131 A, 94 A, and 335 A. Their apparent speeds are estimated to be roughly 126-210 km/s, which could be regarded as the high-velocity upflows caused by the chromospheric evaporation. The nonthermal energy flux is estimated to be about 5.7x10^10 erg/s/cm^2. The characteristic timescale is roughly equal to 1 minute. All these observational results suggest an explosive chromospheric evaporation during the flare impulsive phase. While a HXR/microwave pulse and a type III radio burst are found simultaneously, indicating that the explosive chromospheric evaporation is driven by the nonthermal electron.

Recent observations of active galactic nuclei (AGNs) have shown a high Fe~II/Mg~II line-flux ratio in their broad-line regions, nearly independent of redshift up to $z \gtrsim 6$. The high flux ratio requires rapid production of iron in galactic nuclei to reach an abundance ratio of ${\rm [Fe/Mg]} \gtrsim 0.2$ as high as those observed in matured galaxies in the local universe. We propose a possible explanation of rapid iron enrichment in AGNs by massive star formation that follows a top-heavy initial mass function (IMF) with a power-law index of $\Gamma$ larger than the canonical value of $\Gamma=-2.35$ for a Salpeter IMF. Taking into account metal production channels from different types of SNe, we find that the high value of ${\rm [Fe/Mg]} \gtrsim 0.2$ requires the IMF to be characterized with $\Gamma \gtrsim -1$ ($\Gamma \gtrsim 0$) and a high-mass cutoff at $M_{\rm max} \simeq 100$--$150~{\rm M_\odot}$ $(M_{\rm max} \gtrsim 250~{\rm M_\odot})$. Given the conditions, core-collapse SNe with $M_\ast \gtrsim 70~{\rm M_\odot}$ and pair-instability SNe give a major contribution for iron enrichment. Such top-heavy stellar IMFs would be a natural consequence from mass growth of stars formed in dense AGN disks under Bondi-like gas accretion that is regulated by feedback at $M_\ast \gtrsim 10~{\rm M_\odot}$. The massive stellar population formed in AGN disks also leave stellar-mass black hole remnants, whose mergers associated with gravitational-wave emission account for at most 10 \% of the merger rate inferred from LIGO/Virgo observations to simultaneously explain the high ${\rm [Fe/Mg]}$ ratio with metal ejection.

We investigate the evolution of dust and gas in the vicinity of local pressure enhancements (pressure bumps) in a protoplanetary disc (PPD) with turbulence due to the Vertical Shear Instability (VSI). We perform global 2D axisymmetric and 3D simulations of dust and gas for a range of values for Z (ratio of dust-to-gas (d/g) surface mass densities or metallicity), particle Stokes numbers tau, and pressure bump amplitude A. Dust feedback onto the gas is included. For the first time we demonstrate in global 3D simulations the collection of dust in long-lived vortices induced by the VSI. Without a pressure bump and for Z~0.01 and tau~0.01 we find that such vortices reach d/g density ratios slightly below unity in the discs mid-plane, while for Z>0.05 long-lived vortices are largely absent. In presence of a pressure bump, for Z~0.01 & tau~0.01 a dusty vortex forms reaching d/g ratios of a few times unity, such that the SI is expected to develop, before it eventually shears out into a turbulent dust ring. For Z~0.03 this occurs for tau~0.005, with a weaker, more short-lived vortex, while for larger tau only a turbulent dust ring forms. For Z>0.03 we find that the dust ring becomes increasingly axisymmetric for increasing tau and d/g ratios reach ~1 for tau>0.005. Furthermore, the disc's vertical mass flow profile is strongly affected by dust for Z>0.03, such that gas is transported inward near the mid-plane and outward at larger heights, i.e. the reversed situation compared to simulations with zero or small amounts of dust. Viscous $\alpha$-values drop moderately as 0.001-0.0001 for increasing Z=0-0.05. Our results suggest that the VSI can play an active role in planetesimal formation through the formation of vortices for plausible values of Z and tau. Also it may provide a natural explanation for the presence/absence of asymmetries of observed dust rings in PPDs, depending on the value of Z.

Relatively little is understood about the atmospheric composition of temperate to warm exoplanets (equilibrium temperature $T_{\rm eq}<$ 1000 K), as many of them are found to have uncharacteristically flat transmission spectra. Their flattened spectra are likely due to atmospheric opacity sources such as planet-wide photochemical hazes and condensation clouds. We compile the transmission spectra of 23 warm exoplanets previously observed by the \textit{Hubble Space Telescope} and quantify the haziness of each exoplanet using a normalized amplitude of the water absorption feature ($A_{\rm H}$). By examining the relationships between $A_{\rm H}$ and various planetary and stellar forcing parameters, we endeavor to find correlations of haziness associated with planetary properties. Our analysis shows that the previously identified linear trends between $A_{\rm H}$ and $T_{\rm{eq}}$ or hydrogen-helium envelope mass fraction (f$_{\rm{HHe}}$) break down with the addition of new exoplanet data. Among all the parameters we investigated, atmospheric scale height ($H$), planet gravity ($g_{\rm p}$), and planet density ($\rho_{\rm p}$) hold the most statistically significant linear or linear logarithmic correlations with $A_{\rm H}$ ($p\leq0.02$). We also tentatively identified positive correlations for eccentricity ($e$) and stellar age ($t_{\rm age}$) with $A_{\rm H}$. Specifically, lower $H$, higher $g_{\rm p}$, $\rho_{\rm p}$, $e$, or $t_{\rm age}$ lead to clearer atmospheres. However, none of the parameters show very strong linear correlations with $A_{\rm H}$, suggesting that haziness in warm exoplanets is not simply controlled by any single planetary/stellar parameter. Additional observations and laboratory experiments are needed to fully understand the complex physical and chemical processes that lead to the hazy/cloudy atmospheres in warm exoplanets.

We review and discuss recent results on the search for correlations between astrophysical neutrinos and gamma-ray-detected sources, with many extra-galactic studies reporting potential associations with different types of blazars. We investigate possible dependencies on blazar sub-classes by using the largest catalogues and all the multi-frequency data available. Through the study of similarities and differences in these sources we conclude that blazars come in two distinct flavors: LBLs and IHBLs (low-energy-peaked and intermediate-high-energy-peaked objects). These are distinguished by widely different properties such as the overall spectral energy distribution shape, jet speed, cosmological evolution, broad-band spectral variability, and optical polarization properties. Although blazars of all types have been proposed as neutrino sources, evidence is accumulating in favor of IHBLs being the counterparts of astrophysical neutrinos. If this is indeed the case, we argue that the peculiar observational properties of IHBLs may be indirectly related to proton acceleration to very high energies.

Using X-ray observations from the NuSTAR and Swift satellites, we present temporal and spectral properties of an intermediate polar (IP) IGR J16547-1916. A persistent X-ray period at ~ 546 s confirming the optical spin period obtained from previous observations is detected. The detection of a strong X-ray spin pulse reinforces the classification of this system as an intermediate polar. The lack of orbital or side-band periodicities in the X-rays implies that the system is accreting predominantly via a disk. A variable covering absorber appears to be responsible for the spin pulsations in the low energy range. In the high energy band, the pulsations are likely due to the self occultation of tall shocks above the white dwarf surface. The observed double-humped X-ray spin pulse profile indicates two-pole accretion geometry with tall accretion regions in short rotating IP IGR J16547-1916. We present the variation of the spin pulse profile over an orbital phase to account for the effects of orbital motion on the spin pulsation. X-ray spectra obtained from the contemporaneous observations of Swift and NuSTAR in the 0.5-78.0 keV energy band are modeled with a maximum temperature of 31 keV and a blackbody temperature of 64 eV, along with a common column density of $1.8\times10^{23} cm^{-2}$ and a power-law index of -0.22 for the covering fraction. An additional Gaussian component and a reflection component are needed to account for a fluorescent emission line at 6.4 keV and the occurrence of X-ray reflection in the system. We also present the spin phase-resolved spectral variations of IGR J16547-1916 in the 0.5-78.0 keV energy band and find dependencies in the X-ray spectral parameters during the rotation of the white dwarf.

We present the evolution of the star-formation dispersion - stellar mass relation ($\sigma_{SFR}$-M$_{\star}$) in the DEVILS D10 region using new measurements derived using the ProSpect spectral energy distribution fitting code. We find that $\sigma_{SFR}$-M$_{\star}$ shows the characteristic 'U-shape' at intermediate stellar masses from 0.1<z<0.7 for a number of metrics, including using the deconvolved intrinsic dispersion. A physical interpretation of this relation is the combination of stochastic star-formation and stellar feedback causing large scatter at low stellar masses and AGN feedback causing asymmetric scatter at high stellar masses. As such, the shape of this distribution and its evolution encodes detailed information about the astrophysical processes affecting star-formation, feedback and the lifecycle of galaxies. We find that the stellar mass that the minimum ${\sigma}_{SFR}$ occurs evolves linearly with redshift, moving to higher stellar masses with increasing lookback time and traces the turnover in the star-forming sequence. This minimum ${\sigma}_{SFR}$ point is also found to occur at a fixed specific star-formation rate (sSFR) at all epochs (sSFR~10$^{-9.6}$yr$^{-1}$). The physical interpretation of this is that there exists a maximum sSFR at which galaxies can internally self-regulate on the tight sequence of star-formation. At higher sSFRs, stochastic stellar processes begin to cause galaxies to be pushed both above and below the star-forming sequence leading to increased SFR dispersion. As the Universe evolves, a higher fraction of galaxies will drop below this sSFR threshold, causing the dispersion of the low-stellar mass end of the star-forming sequence to decrease with time.

The Seyfert-1 galaxy NGC 3516 has undergone major spectral changes in recent years. In 2017 we obtained Chandra, NuSTAR, and Swift observations during its new low-flux state. Using these observations we model the spectral energy distribution (SED) and the intrinsic X-ray absorption, and compare the results with those from historical observations taken in 2006. We thereby investigate the effects of the changing-look phenomenon on the accretion-powered radiation and the ionized outflows. Compared to its normal high-flux state in 2006, the intrinsic bolometric luminosity of NGC 3516 was lower by a factor of 4 to 8 during 2017. Our SED modeling shows a significant decline in the luminosity of all the continuum components from the accretion disk and the X-ray source. As a consequence, the reprocessed X-ray emission lines have also become fainter. The Swift monitoring of NGC 3516 shows remarkable X-ray spectral variability on short (weeks) and long (years) timescales. We investigate whether this variability is driven by obscuration or the intrinsic continuum. We find that the new low-flux spectrum of NGC 3516, and its variability, do not require any new or variable obscuration, and instead can be explained by changes in the ionizing SED that result in lowering of the ionization of the warm-absorber outflows. This in turn induces enhanced X-ray absorption by the warm-absorber outflows, mimicking the presence of new obscuring gas. Using the response of the ionized regions to the SED changes, we place constraints on their density and location.

Synchronous satellites of Venus have long been thought unstable, but we use Poincare's surface of section technique to show that synchronous quasi-satellites orbiting just outside Venus' Hill sphere are quite stable, at least for centuries. Such synchrosats always remain within a few degrees of Venus' equator, and drift very slowly in longitude. These synchrosats could be useful for continuous monitoring of points on Venus' surface, such as active landforms or long-lived landers.

We compiled available mass and redshift z data for supermassive black holes (SMBHs) at z >7.5, ranging in mass over 2 orders of magnitude and in age by nearly 300 million years. The data reveal that a large subset covering the entire age spectrum has markedly similar masses. The most likely implication is that SMBHs within the subset had seeds with similar masses and formed essentially concurrently. Based on this inference, the data of the subset are used derive quantitative empirical relations that provide insights into the origins of black hole seeds and constraints for models of seed formation. The relationships are tested and applied to thousands of SMBHs at nearly all redshifts. The results show that the masses of SMBHs > a million solar masses are accounted for with seeds formed at or near z=30 and ranging from Sun's mass to about 50-thousand solar masses. Apparently, the seeds grew at an exponentially increasing accretion rate that reached a maximum near z=7 and decreased thereafter. From z=30 to 15, the average accretion rate either exceeded the Eddington limit by a factor of 2 or less, or the radiative efficiency was less than its canonical value of 0.1 and increased thereafter. About half of the growth apparently occurs from z=30 to 3.5 at an average rate of about 614 solar masses per million years per unit solar seed mass, and the rest in the succeeding 12-billion years. The maximum mass that a black hole can accrete is about 2.35-million times the seed mass, and the largest observable black hole should not exceed about 100-billion solar masses. The seed of Sagittarius A* is inferred to have had a mass a few times Sun's mass and may have stopped accreting recently having achieved its maximum growth potential.

Active galactic nuclei (AGN) are typically identified through radio, mid-infrared, or X-ray emission or through the presence of broad and/or narrow emission lines. AGN can also leave an imprint on a galaxy's spectral energy distribution (SED) through the re-processing of photons by the dusty torus. Using the SED fitting code ProSpect with an incorporated AGN component, we fit the far ultraviolet to far-infrared SEDs of $\sim$494,00 galaxies in the D10-COSMOS field and $\sim$230,000 galaxies from the GAMA survey. By combining an AGN component with a flexible star formation and metallicity implementation, we obtain estimates for the AGN luminosities, stellar masses, star formation histories, and metallicity histories for each of our galaxies. We find that ProSpect can identify AGN components in 91 per cent of galaxies pre-selected as containing AGN through narrow-emission line ratios and the presence of broad lines. Our ProSpect-derived AGN luminosities show close agreement with luminosities derived for X-ray selected AGN using both the X-ray flux and previous SED fitting results. We show that incorporating the flexibility of an AGN component when fitting the SEDs of galaxies with no AGN has no significant impact on the derived galaxy properties. However, in order to obtain accurate estimates of the stellar properties of AGN host galaxies, it is crucial to include an AGN component in the SED fitting process. We use our derived AGN luminosities to map the evolution of the AGN luminosity function for $0<z<2$ and find good agreement with previous measurements and predictions from theoretical models.

Astrometric detection of exoplanets by stellar reflex motion will be made possible for giant planets by the recent Gaia mission and for Earth-like planets by proposed missions such as LUVOIR. At the theoretical limits, astrometry would allow for the detection of smaller planets than previously seen by current exoplanet search methods, but stellar activity may make these theoretical limits unreachable. Astrometric jitter of a Sun-like star due to magnetic activity in its photosphere induces apparent variability in the photocenter of order $0.5\ \textrm{mR}_\odot$. This jitter creates a fundamental astrophysical noise floor preventing detection of lower mass planets in a single spectral band. By injecting planet orbits into simulated solar data at five different passbands, we investigate mitigation of this fundamental astrometric noise using correlations across passbands. For a true solar analog and a planet at 1 au semi-major axis, the 5-sigma detection limit set by stellar activity for an ideal telescope at the best single passband is $0.01$ Earth masses. We found that pairs of passbands with highly correlated astrometric jitter due to stellar activity but with less motion in the redder band enable higher precision measurements of the common signal from the planet. Using this method improves detectable planet masses at 1 au by up to a factor of $8$, corresponding to at best $0.004$ Earth masses for a Sun-like star with a perfect telescope. Given these results, we recommend that future astrometry missions consider proceeding with two or more passbands to reduce noise due to stellar activity.

Gliese 86 is a nearby K dwarf hosting a giant planet on a $\approx$16-day orbit and an outer white dwarf companion on a $\approx$century-long orbit. In this study we combine radial velocity data (including new measurements spanning more than a decade) with high angular resolution imaging and absolute astrometry from Hipparcos and Gaia to measure the current orbits and masses of both companions. We then simulate the evolution of the Gl 86 system to constrain its primordial orbit when both stars were on the main sequence; the closest approach between the two stars was then about $9\,$AU. Such a close separation limited the size of the protoplanetary disk of Gl 86 A and dynamically hindered the formation of the giant planet around it. Our measurements of Gl 86 B and Gl 86 Ab's orbits reveal Gl 86 as a system in which giant planet formation took place in a disk truncated at $\approx$2$\,$AU. Such a disk would be just big enough to harbor the dust mass and total mass needed to assemble Gl 86 Ab's core and envelope, assuming a high disk accretion rate and a low viscosity. Inefficient accretion of the disk onto Gl 86 Ab, however, would require a disk massive enough to approach the Toomre stability limit at its outer truncation radius. The orbital architecture of the Gl 86 system shows that giant planets can form even in severely truncated disks and provides an important benchmark for planet formation theory.

We demonstrate the potential of line intensity mapping to place constraints on the initial mass function (IMF) of Population III (Pop III) stars via measurements of the mean He II 1640A/H$\alpha$ emission line ratio. We extend the 21cmFAST code with modern high-redshift galaxy formation and photoionization models, and estimate the line emission from Pop II and Pop III galaxies at redshifts $5 \le z \le 20$. In our models, mean ratio values of He II/H$\alpha \gtrsim 0.1$ indicate top-heavy Pop III IMFs with stars of several hundred solar masses, reached at $z \gtrsim 10$ when Pop III stars dominate star formation. A next-generation space mission with capabilities moderately superior to those of CDIM will be able to probe this scenario by measuring the He II and H$\alpha$ fluctuation power spectrum signals and their cross-correlation at high significance up to $z\sim 20$. Moreover, regardless of the IMF, a ratio value of He II/H$\alpha \lesssim 0.01$ indicates low Pop III star formation and, therefore, it signals the end of the period dominated by this stellar population. However, a detection of the corresponding He II power spectrum may be only possible for top-heavy Pop III IMFs or through cross-correlation with the stronger H$\alpha$ signal. Finally, ratio values of $0.01 \lesssim$ He II/H$\alpha$ $\lesssim 0.1$ are complex to interpret because they can be driven by several competing effects. We discuss how various measurements at different redshifts and the combination of the line ratio with other probes can assist in constraining the Pop III IMF in this case.

The Universe is currently undergoing accelerated expansion driven by dark energy. Dark energy's essential nature remains mysterious: one means of revealing it is by measuring the Universe's size at different redshifts. This may be done using the Baryon Acoustic Oscillation (BAO) feature, a standard ruler in the galaxy 2-Point Correlation Function (2PCF). In order to measure the distance scale, one dilates and contracts a template for the 2PCF in a fiducial cosmology, using a scaling factor $\alpha$. The standard method for finding the best-fit $\alpha$ is to compute the likelihood over a grid of roughly 100 values of it. This approach is slow; in this work, we propose a significantly faster way. Our method writes the 2PCF as a polynomial in $\alpha$ by Taylor-expanding it about $\alpha = 1$, exploiting that we know the fiducial cosmology sufficiently well that $\alpha$ is within a few percent of unity. The likelihood resulting from this expansion may then be analytically solved for the best-fit $\alpha$. Our method is 48-85$\times$ faster than a directly comparable approach in which we numerically minimize $\alpha$, and $\sim$$12,000 \times$ faster than the standard iterative method. Our work will be highly enabling for upcoming large-scale structure redshift surveys such as that by Dark Energy Spectroscopic Instrument (DESI).

The gas supply from the cosmic web is the key to sustain star formation in galaxies. It remains to be explored how the cosmic large-scale structure (LSS) effects on galaxy evolution at given local environments. We examine galaxy specific star formation rate as a function of local density in a LSS at $z=0.735$ in the Extended Chandra Deep Field South. The LSS is mapped by 732 galaxies with $R<24$\,mag and redshift at $0.72\le z \le 0.75$ collected from the literature and our spectroscopic observations with Magellan/IMACS, consisting of five galaxy clusters/groups and surrounding filaments over an area of $23.9 \times22.7$\,co-moving\,Mpc$^2$. The spread of spectroscopic redshifts corresponds a velocity dispersion of 494\,km\,s$^{-1}$, indicating the LSS likely to be a thin sheet with a galaxy density $\gtrsim 3.9$ times that of the general field. These clusters/groups in this LSS mostly exhibit elongated morphologies and multiple components connected with surrounding filaments. Strikingly, we find that star-forming galaxies in the LSS keep star formation at the same level as field, and show no dependence on local density but stellar mass. Meanwhile, an increasing fraction of quiescent galaxies is detected at increasing local density in both the LSS and the field, consistent with the expectation that galaxy mass and local dense environment hold the key to quench star formation. Combined together, we conclude that the cosmic environment of the LSS overtakes the local environment in remaining galaxy star formation to the level of the field.

Super-massive black holes (SMBHs) spend most of their lifetime accreting from ambient matter at a rate well below the Eddington limit, manifesting themselves as low-luminosity active galactic nuclei (LLAGNs). The prevalence of a hot wind from LLAGNs is a generic prediction by theories and numerical simulations of black hole accretion and is recently becoming a crucial ingredient of AGN kinetic feedback in cosmological simulations of galaxy evolution. However, direct observational evidence for this hot wind is still scarce. In this work, from high-resolution {\it Chandra} grating spectra of the LLAGN in NGC 7213, a nearby Sa galaxy hosting a $\sim10^8\rm~M_\odot$ SMBH, we identify significant Fe XXVI Ly$\alpha$ and Fe XXV K$\alpha$ emission lines with a blueshifted line-of-sight velocity of $\sim1100\rm~km~s^{-1}$. The measured flux ratio between Fe XXVI Ly$\alpha$ and Fe XXV K$\alpha$ suggests that these lines arise from a $\sim16$ keV hot plasma. By confronting these spectral features with synthetic X-ray spectra based on our custom magnetohydrodynamical simulations, we find that the high-velocity, hot plasma in this LLAGN is naturally explained by the putative hot wind driven by the hot accretion flow onto the SMBH. Alternative plausible origins of this hot plasma, including stellar activities, AGN photoionization and the hot accretion flow itself, can be quantitatively ruled out. A mass outflow rate $\sim0.08{\rm~M_{\odot}~yr^{-1}}$ is inferred for the hot wind. This is comparable to an independent estimate of the mass inflow rate, consistent with the prediction of the theory of hot wind. The wind carries a kinetic energy of $\sim3\times10^{42}\rm~erg~s^{-1}$, accounting for 15% of the LLAGN's bolometric luminosity, and a momentum flux of $\sim4\times10^{33}\rm~g~cm~s^{-1}$, about 6 times the photon momentum flux. (abridged)

Radio Astronomy is undergoing a renaissance, as the next-generation of instruments provides a massive leap forward in collecting area and therefore raw sensitivity. However, to achieve this theoretical level of sensitivity in the science data products we need to address the much more pernicious systematic effects, which are the true limitation. These become all the more significant when we consider that much of the time used by survey instruments, such as the SKA, will be dedicated to deep surveys. CHILES is a deep HI survey of the COSMOS field, with 1,000 hours of VLA time. We present our approach for creating the image cubes from the first Epoch, with discussions of the methods and quantification of the data quality from 946 to 1420MHz -- a redshift range of 0.5 to 0. We layout the problems we had to solve and describe how we tackled them. These are of importance as CHILES is the first deep wideband multi-epoch HI survey and it has relevance for ongoing and future surveys. We focus on the accumulated systematic errors in the imaging, as the goal is to deliver a high-fidelity image that is only limited by the random thermal errors. To understand and correct these systematic effects we ideally manage them in the domain in which they arise, and that is predominately the visibility domain. CHILES is a perfect test bed for many of the issues we can expect for deep imaging with the SKA or ngVLA and we discuss the lessons we have learned.

We studied the presence and spatiotemporal characteristics and evolution of the variations in the differential rotation rates and radial magnetic fields in the Schwabe and Quasi-biennial-oscillation (QBO) timescales. To achieve these objectives, we used rotation rate residuals and radial magnetic field data from the Michelson Doppler Imager on the Solar and Heliospheric Observatory and the Helioseismic and Magnetic Imager on the Solar Dynamics Observatory, extending from May 1996 to August 2020, covering solar cycles 23 and 24, respectively. Under the assumption that the radial surface magnetic field is non-local and the differential rotation is symmetric around the equator, our results suggest that the source region of the Schwabe cycle is confined between $\sim$30$^{\circ}$ N and S throughout the convection zone. As for the source region of the QBO, our results suggest that it is below 0.78R$_{\odot}$.

The unusually large, ``feeble dwarf" galaxy Crater II, with its small velocity dispersion, $\simeq 3$ km/s, defies expectations that low mass galaxies should be small and dense. Here we examine its unusual properties in the context of ``Wave Dark Matter", combining the latest stellar and velocity dispersion profiles for Crater II, finding a prominent dark core of radius $\simeq 0.71^{+0.09}_{-0.08}$ kpc, surrounded by a low density halo, with a visible transition between the core and the halo. This observed behaviour is very similar to the distinctive core-halo profile structure of dark matter as a Bose-Einstein condensate, $\psi$DM, where the ground state forms a prominent soliton core, surrounded by a tenuous halo of interfering waves, with a marked density transition predicted between the soliton and the halo. Crater II conforms well this distinctive $\psi$DM prediction, with consistency found between its large core and low velocity dispersion for a boson mass of $m_\psi c^2\simeq (1.9 \pm 0.3) \times10^{-22}$ eV. Similar core-halo structure is also apparent in most dwarf spheroidal galaxies (dSph), but with typically smaller cores, $\simeq 0.25$ kpc and higher velocity dispersions, $\simeq 9$km/s. We argue that Crater II may have have been a more typical dSph dwarf that has lost most of its halo mass to tidal stripping in the context of $\psi$DM, resulting in a factor 3 reduction in velocity dispersion causing a threefold expansion of the soliton core, following the inverse scaling between velocity and de Broglie wavelength required by the Uncertainty Principle. This tidal origin for Crater II is supported by its small pericenter of $\simeq 20$ kpc, now established by GAIA, implying significant tidal stripping by the Milky Way.

Sokolsky and D_Avignon have recently reported an examination of a variety of measurements that relate to the question as to how the mass composition of the highest-energy cosmic-rays evolves with energy. They assert that cosmic rays arriving from the Northern Hemisphere have a different mass composition from those arriving from the Southern Hemisphere, implying a diversity of sources of high-energy cosmic-rays in the two hemispheres. Were this conclusion to be correct, it would have profound implications for theories of cosmic-ray origin and would influence planning of future projects. Their claim thus merits careful scrutiny. In this paper their analysis is examined in detail with the verdict being that evidence for a North/South difference is not proven, a conclusion supported by other data from the Northern and Southern Hemispheres. However, what is of major importance is that the study of Sokolsky and D_Avignon provides long-awaited confirmation of the claim that the mean mass of cosmic rays increases with energy above ~3 E eV made by the Pierre Auger Collaboration in 2014.

Globular clusters are known to host an unusually large population of millisecond pulsar when compared to the Galactic disk. This is thanks to the high rate of dynamical encounters occurring in the clusters that can create the conditions to efficiently recycle neutron stars into millisecond pulsars. The result is a rich population of pulsars with properties and companions difficult or impossible to replicate in the Galactic disk. For these reasons, globular clusters have been and still are a prime target of searches for new and exciting pulsars. Because of their large distances, the limiting factor inhibiting these discoveries is the telescope sensitivity. The MeerKAT radio telescope, a 64-dish interferometer in South Africa, guarantees unrivalled sensitivity for globular clusters in the southern sky. Observations of well-studied globular clusters with MeerKAT have already returned more than 35 new pulsars with many more expected. These exciting discoveries will help us to understand more about the neutron star equation of state, stellar evolution, accretion physics and to hunt for intermediate mass black holes. In this talk I will present the prospects and current discoveries of the globular cluster working group in the MeerTIME and TRAPUM programmes.

We analyse the evolution of the oxygen abundance gradient of star-forming galaxies with stellar mass Mstar > 10^9 Mo in the EAGK simulation over the redshift range z=[0, 2.5]. We find that the median metallicity gradient of the simulated galaxies is close to zero at all z, whereas the scatter around the median increases with z. The metallicity gradients of individual galaxies can evolve from strong to weak and vice-versa, since mostly low-metallicity gas accretes onto the galaxy, resulting in enhanced star formation and ejection of metal enriched gas by energy feedback. Such episodes of enhanced accretion, mainly dominated by major mergers, are more common at higher z, and hence contribute to increasing the diversity of gradients. For galaxies with negative metallicity gradients, we find a redshift evolution of ~ -0.03 dex/kpc/\delta z$. A positive mass dependence is found at z< 0.5, which becomes slightly stronger for higher redshifts and, mainly, for Mstar < 10^9.5 Mo. Only galaxies with negative metallicity gradients define a correlation with galaxy size, consistent with an inside-out formation scenario. Our findings suggest that major mergers and/or significant gas accretion can drive strong negative or positive metallicity gradients. The first ones are preferentially associated with disc-dominated galaxies, and the second ones with dispersion-dominated systems. The comparison with forthcoming observations at high redshift will allow a better understanding of the potential role of metallicity gradients as a chemical probe of galaxy formation.

We review problems raised by Fe II emission in AGNs and address the question of its relationship to other broad-line region (BLR) lines. The self-shielding, stratified, BLR model of Gaskell, Klimek & Nazarova (2007; GKN) predicts that Fe II emission comes from twice the radius of H$\beta$, in agreement with widths of lines of Fe II lines being only 70% of the widths of H$\beta$. This disagrees with some reverberation mapping results which have suggested that Fe II and H$\beta$ arise at similar radii. The highest quality reverberation mapping, however, supports the predictions that Fe II comes from twice the radius of H$\beta$. We suggest that lower quality reverberation mapping of Fe II is biased to give too small lags. We conclude that, in agreement with the GKN model, the region emitting Fe II is the outermost part of the BLR just inside the surrounding dust. The model naturally gives the Doppler broadening require by models of Fe II emission. Optical Fe II emission implies typical reddenings of E(B-V) ~ 0.20. This helps explain the ratio of UV to optical Fe II emission. Simulations show that the amplitude of Fe II variability is consistent with being the same as for H$\beta$ variability. The FeII/H$\beta$ ratio is a good proxy for the Eddington ratio. The ratio might be driven in part by the strong soft X-ray excess because the X-rays destroy grains and release iron into the gas phase. We propose that the correlation of Fe II strength with radio and host galaxy properties is a result of AGN downsizing.

We study the tidal interaction of galaxies in the Eridanus supergroup, using HI data from the pre-pilot survey of WALLABY (Widefield ASKAP L-band Legacy All-sky Blind surveY). We obtain optical photometric measurements and quantify the strength of tidal perturbation using a tidal parameter $S_{sum}$. For low-mass galaxies of $M_* \lesssim 10^9 M_\odot$, we find a dependence of decreasing HI-to-optical disk size ratio with increasing $S_{sum}$, but no dependence of HI spectral line asymmetry with $S_{sum}$. This is consistent with the behavior expected under tidal stripping. We confirm that the color profile shape and color gradient depend on the stellar mass, but there is additional correlation of low-mass galaxies having their color gradients within $2R_{50}$ increasing with higher $S_{sum}$. For these low-mass galaxies, the dependence of color gradients on $S_{sum}$ is driven by color becoming progressively redder in the inner disk when tidal perturbations are stronger. For high-mass galaxies, there is no dependence of color gradients on $S_{sum}$, and we find a marginal reddening throughout the disks with increasing $S_{sum}$. Our result highlights tidal interaction as an important environmental effect in producing the faint end of the star formation suppressed sequence in galaxy groups.

The radial velocity method is amongst the most robust and most established means of detecting exoplanets. Yet, it has so far failed to detect circumbinary planets despite their relatively high occurrence rates. Here, we report velocimetric measurements of Kepler-16A, obtained with the SOPHIE spectrograph, at the Observatoire de Haute-Provence's 193cm telescope, collected during the BEBOP survey for circumbinary planets. Our measurements mark the first radial velocity detection of a circumbinary planet, independently determining the mass of Kepler-16~(AB)~b to be $0.313 \pm 0.039\,{\rm M}_{\rm Jup}$, a value in agreement with eclipse timing variations. Our observations demonstrate the capability to achieve photon-noise precision and accuracy on single-lined binaries, with our final precision reaching $\rm 1.5~m\,s^{-1}$ on the binary and planetary signals. Our analysis paves the way for more circumbinary planet detections using radial velocities which will increase the relatively small sample of currently known systems to statistically relevant numbers, using a method that also provides weaker detection biases. Our data also contain a long-term radial velocity signal, which we associate with the magnetic cycle of the primary star.

We report on our second campaign to search for old nova shells around cataclysmic variables (CVs). Our aim was to test the theory that nova eruptions cause cycles in the mass transfer rates of CVs. These mass transfer cycles change the behaviour of CVs during their inter-eruption periods. We examined H-alpha images of 47 objects and found no new shells around any of the targets. Combining our latest results with our previous campaign (Sahman et al. 2015), and the searches by Schmidtobreick et al. (2015) and Pagnotta & Zurek (2016), we estimate that the nova-like phase of the mass transfer cycle lasts approximately 3,000 years.

Ultra-short orbital period contact binaries (Porb < 0.26 d) host some of the smallest and least massive stars. These systems are faint and rare, and it is believed that they have reached a contact configuration after several Gyrs of evolution via angular momentum loss, mass transfer and mass loss through stellar wind processes. This study is conducted in the frame of Contact Binaries Towards Merging (CoBiToM) Project and presents the results from light curve and orbital analysis of 30 ultra-short orbital period contact binaries, with the aim to investigate the possibility of them being red nova progenitors, eventually producing merger events. Approximately half of the systems exhibit orbital period modulations, as a result of mass transfer or mass loss processes. Although they are in contact, their fill-out factor is low (less than 30 per cent), while their mass ratio is larger than the one in longer period contact binaries. The present study investigates the orbital stability of these systems and examines their physical and orbital parameters in comparison to those of the entire sample of known and well-studied contact binaries, based on combined spectroscopic and photometric analysis. It is found that ultra-short orbital period contact binaries have very stable orbits, while very often additional components are gravitationally bound in wide orbits around the central binary system. We confirmed that the evolution of such systems is very slow, which explains why the components of ultra-short orbital period systems are still Main Sequence stars after several Gyrs of evolution.

Traditionally, ground-based spectrophotometric observations probing transiting exoplanet atmospheres have employed a linear map between comparison and target star light curves (e.g. via differential spectrophotometry) to correct for systematics contaminating the transit signal. As an alternative to this conventional method, we introduce a new Gaussian Processes (GP) regression-based method to analyse ground-based spectrophotometric data. Our new method allows for a generalised non-linear mapping between the target transit light curves and the time series used to detrend them. This represents an improvement compared to previous studies because the target and comparison star fluxes are affected by different telluric and instrumental systematics, which are complex and non-linear. We apply our method to six Gemini/GMOS transits of the warm (T$_{\rm eq}$ = 990 K) Neptune HAT-P-26b. We obtain on average $\sim$20 % better transit depth precision and residual scatter on the white light curve compared to the conventional method when using the comparison star light curve as a GP regressor and $\sim$20% worse when explicitly not using the comparison star. Ultimately, with only a cost of 30% precision on the transmission spectra, our method overcomes the necessity of using comparison stars in the instrument field of view, which has been one of the limiting factors for ground-based observations of the atmospheres of exoplanets transiting bright stars. We obtain a flat transmission spectrum for HAT-P-26b in the range of 490-900 nm that can be explained by the presence of a grey opacity cloud deck, and indications of transit timing variations, both of which are consistent with previous measurements.

Blazars are known to emit exceptionally variable non-thermal emission over the wide range (from radio to $\gamma$-rays) of electromagnetic spectrum. We present here the results of our $\gamma$-ray flux and spectral variability study of the blazar Ton 599, which has been recently observed in the $\gamma$-ray flaring state. Using 0.1$-$300 GeV $\gamma$-ray data from the $\it{Fermi}$ Gamma-ray Space Telescope (hereinafter $\it{Fermi}$), we generated one-day binned light curve of Ton 599 for a period of about one-year from MJD 59,093 to MJD 59,457. During this one year period, the maximum $\gamma$-ray flux detected was 2.24 $\pm$ 0.25 $\times$ $10^{-6}$ ph cm$^{-2}$ s$^{-1}$ at MJD 59,399.50. We identified three different flux states, namely, epoch A (quiescent), epoch B (pre-flare) and epoch C (main-flare). For each epoch, we calculated the $\gamma$-ray flux variability amplitude (F$_{var}$) and found that the source showed largest flux variations in epoch C with F$_{var} \sim$ 35%. We modelled the $\gamma$-ray spectra for each epoch and found that the Log-parabola model adequately describes the $\gamma$-ray spectra for all the three epochs. We estimated the size of the $\gamma$-ray emitting region as 1.03 $\times$ $10^{16}$ cm and determined that the origin of $\gamma$-ray radiation, during the main-flare, could be outside of the broad line region.

NASA's Kepler mission observed background regions across its field of view for more than three consecutive years using custom designed super apertures (EXBA masks). Since these apertures were designed to capture a region of the sky rather than single targets, the Kepler Science Data Processing pipeline produced Target Pixel Files, but did not produce light curves for the sources within these background regions. In this work we produce light curves for $9,327$ sources observed in the EXBA masks. These light curves are generated using aperture photometry estimated from the instrument's Pixel Response Function (PRF) profile computed from Kepler's full-frame images. The PRF models enable the creation of apertures that follow the characteristic shapes of the PSF in the image and the computation of flux completeness and contamination metrics. The light curves are available at MAST as a High Level Science Product (kbonus-apexba). Alongside this dataset, we present kepler-apertures, a Python library to compute PRF models and use them to perform aperture photometry on Kepler-like data. Using light curves from the EXBA masks we found an exoplanet candidate around Gaia EDR3 2077240046296834304 consistent with a large planet companion with a $0.81 R_J$ radius. Additionally, we report a catalog of 69 eclipsing binaries. We encourage the community to exploit this new dataset to perform in depth time domain analysis, such as eclipsing binaries demographic and others.

Our peculiar motion in a homogeneous and isotropic universe imprints a dipole in the cosmic microwave background (CMB) temperature field and similarly imprints a dipole in the distribution of extragalactic radio sources on the sky. Each of these effects have been measured, however each of these measurements give different results for the velocity of our motion through the Universe: the radio dipole measurements finds the speed of our motion to be around three times larger than that of the CMB. Here we show the effects of the previously unconstrained lensing dipole, whereby necessarily local structures (required for large angular lensing scales) will distort the distribution of radio sources on the sky. We find that the inclusion of these effects does not reduce the tension between the CMB and radio source dipole measurements however without their inclusion future extragalactic number counts could lead to incorrect inferences of our peculiar motion. In addition we can constrain the size of the lensing dipole to be $\kappa < 3 \cdot 10^{-2}$ at the $2 \sigma$ level.

Accurate computations of spectral distortions of the cosmic microwave background (CMB) are required for constraining energy release scenarios at redshifts $z\gtrsim 10^3$. The existing literature focuses on distortions that are small perturbations to the background blackbody spectrum. At high redshifts ($z\gtrsim 10^6$), this assumption can be violated, and the CMB spectrum can be significantly distorted at least during part of its cosmic evolution. In this paper, we carry out accurate thermalization computations, evolving the distorted CMB spectrum in a general, fully non-linear way, consistently accounting for the time-dependence of the injection process, modifications to the Hubble expansion rate and relativistic Compton scattering. Specifically, we study single energy injection and decaying particle scenarios, discussing constraints on these cases. We solve the thermalization problem using two independent numerical approaches that are now available in {\tt CosmoTherm} as dedicated setups for computing CMB spectral distortions in the large distortion regime. New non-linear effects at low frequencies are furthermore highlighted, showing that these warrant a more rigorous study. This work eliminates one of the long-standing simplifications in CMB spectral distortion computations, which also opens the way to more rigorous treatments of distortions induced by high-energy particle cascade, soft photon injection and in the vicinity of primordial black holes.

Gravitational-wave population studies have become more important in gravitational-wave astronomy because of the rapid growth of the observed catalog. In recent studies, emulators based on different machine learning techniques are used to emulate the outcomes of the population synthesis simulation with fast speed. In this study, we benchmark the performance of two emulators that learn the truncated power-law phenomenological model by using Gaussian process regression and normalizing flows techniques to see which one is a more capable likelihood emulator in the population inference. We benchmark the characteristic of the emulators by comparing their performance in the population inference to the phenomenological model using mock and real observation data. Our results suggest that the normalizing flows emulator can recover the posterior distribution by using the phenomenological model in the population inference with up to 300 mock injections. The normalizing flows emulator also underestimates the uncertainty for some posterior distributions in the population inference on real observation data. On the other hand, the Gaussian process regression emulator has poor performance on the same task and can only be used effectively in low-dimension cases.

Recently, a low-$z$ measurement of the Hubble constant, $H_0 = 73.04 \pm 1.04 {\rm ~km/s/Mpc}$, was reported by the SH0ES Team. The long-standing Hubble tension, i.e. the difference between the Hubble constant from the local measurements and that inferred from the cosmic microwave background data based on the $\Lambda$CDM model, was further strengthened. There are many cosmological models modifying the cosmology after and around the recombination era to alleviate this tension. In fact, some of the models alter the small-scale fluctuation amplitude relative to larger scales, and thus require a significant modification of the primordial density perturbation, especially the scalar spectral index, $n_s$. In certain promising models, $n_s$ is favored to be larger than the $\Lambda$CDM prediction, and even the scale-invariant one, $n_s=1$, is allowed. In this Letter, we focus on the very early Universe models to study the implication of such unusual $n_s$. In particular, we find that an axiverse with an axion in the equilibrium distribution during inflation can be easily consistent with $n_s = 1$. This is because the axion behaves as a curvaton with mass much smaller than the inflationary Hubble parameter. We also discuss other explanations of $n_s$ different from that obtained based on the $\Lambda$CDM.

Fast radio bursts (FRBs) are millisecond-duration signals, to be highly dispersed at distant galaxies, the physical origin of which is still challenging. Coherent curvature emission by bunches powered, e.g., by starquakes, has already been proposed for repeating FRBs, with a friendly nature of understanding the narrowband radiation exhibiting time-frequency drifting. Recently, a highly active FRB source, FRB 20201124A, was reported to enter a newly active episode and emit at least some highly circle-polarized bursts. We revisit the polarized FRB emission here, investigating especially the production mechanisms of a high circular polarization by deriving intrinsic and propagation effects. The intrinsic mechanisms by invoking charged bunches are approached in two scenarios of coherence: curvature radiation (CR) and inverse Compton scattering (ICS), and consequently, a high circular polarization could naturally be explained by the coherent summation of outcome waves, generated or scattered by bunches, with different phases and electric vectors. Cyclotron resonance can result in an absorption of R-mode photons at lower altitude region of magnetosphere, and an FRB should then be emitted from a higher region if the waves are of strong linear polarization. Circularly polarized components could be produced from Faraday conversion exhibiting a $\lambda^3$-oscillation, but the average circular polarization fraction depends only on the income wave, indicating a possibility of highly circle-polarized income wave. The analysis could be welcome if extremely high (e.g., almost 100\%) circular polarization from repeating FRBs would be detected in the future. The production of a bulk of energetic bunches in pulsar-like magnetosphere is discussed finally, which is relevant to the nature of FRB central engine.

Remnant radio galaxies (RRGs), characterized by the cessation of AGN activity, represent a short-lived last phase of radio galaxy's life-cycle. Hitherto, searches for RRGs, mainly based on the morphological criteria, have identified large angular size sources resulting into a bias towards the remnants of powerful FR-II radio galaxies. In this study we make the first attempt to perform a systematic search for RRGs of small angular sizes ($<$30$^{\prime\prime}$) in the XMM-LSS field. By using spectral curvature criterion we discover 48 remnant candidates exhibiting strong spectral curvature ${\it i.e.}$; ${\alpha}_{\rm 150~MHz}^{\rm 325~MHz}$ - ${\alpha}_{\rm 325~MHz}^{\rm 1.4~GHz}$ $\geq$0.5. Spectral characteristics at higher frequency regime ($>$1.4 GHz) indicate that some of our remnant candidates can depict recurrent AGN activity with an active core. We place an upper limit on the remnant fraction ($f_{\rm rem}$) to be 3.9$\%$, which increases to 5.4$\%$ if flux cutoff limit of S$_{\rm 150~MHz}$ $\geq$10 mJy is considered. Our study unveils, hitherto unexplored, a new population of small-size ($<$200 kpc) remnant candidates that are often found to reside in less dense environments and at higher redshifts ($z$) $>$1.0. We speculate that a relatively shorter active phase and/or low jet power can be plausible reasons for the small size of remnant candidates.

(Abridged) Recent surveys of young star formation regions have shown that the average Class II object does not have enough dust mass to make the cores of giant planets. Younger Class 0/I objects have enough dust in their embedded disk, which begs the questions: can the first steps of planet formation occur in these younger systems? The first step is building the first planetesimals, generally believed to be the product of the streaming instability. Hence the question can be restated: are the physical conditions of embedded disks conducive to the growth of the streaming instability? Here we model the collapse of a `dusty' proto-stellar cloud to show that if there is sufficient drift between the falling gas and dust, regions of the embedded disk can become sufficiently enhanced in dust to drive the streaming instability. We include four models, three with different dust grain sizes and one with a different initial cloud angular momentum to test a variety of collapse trajectories. We find a `sweet spot' for planetesimal formation for grain sizes of a few 10s of micron since they fall sufficiently fast relative to the gas to build a high dust-to-gas ratio along the disk midplane, but have slow enough radial drift speeds in the embedded disk to maintain the high dust-to-gas ratio. Unlike the gas, which is held in hydrostatic equilibrium for a time due to gas pressure, the dust can begin collapsing from all radii at a much earlier time. The streaming instability can produce at least between 7-35 M$_\oplus$ of planetesimals in the Class 0/I phase of our smooth embedded disks, depending on the size of the falling dust grains. This first generation of planetesimals could represent the first step in planet formation, and occurs earlier in the lifetime of the young star than is traditionally thought.

Dynamically excited objects within the Kuiper belt show a bimodal distribution in their surface colors, and these differing surface colors may be a tracer of where these objects formed. In this work we explore radial color distributions in the primordial planetesimal disk and implications for the positions of ice line/color transitions within the Kuiper belt's progenitor populations. We combine a full dynamical model of the Kuiper belt's evolution due to Neptune's migration with precise surface colors measured by the Colours of the Outer Solar System Origins Survey in order to examine the true color ratios within the Kuiper belt and the ice lines within the primordial disk. We investigate the position of a dominant, surface color changing ice-line, with two possible surface color layouts within the initial disk; (1) inner neutral surfaces and outer red, and (2) inner red surfaces and outer neutral. We performed simulations with a primordial disk that truncates at 30 au. By radially stepping the color transition out through 0.5 au intervals we show that both disk configurations are consistent with the observed color fraction. For an inner neutral, outer red primordial disk we find that the color transition can be at $28^{+2}_{-3}$ au at a 95% confidence level. For an inner red, outer neutral primordial disk the color transition can be at $27^{+3}_{-3}$ au at a 95% confidence level.

Lorentz invariance violation (LIV) is often described by dispersion relations of the form $E_i^2=m_i^2+p_i^2+\delta_{i,n} E^{2+n}$ with delta different based on particle type $i$, with energy $E$, momentum $p$ and rest mass $m$. Kinematics and energy thresholds of interactions are modified once the LIV terms become comparable to the squared masses of the particles involved. Thus, the strongest constraints on the LIV coefficients $\delta_{i,n}$ tend to come from the highest energies. At sufficiently high energies, photons produced by cosmic ray interactions as they propagate through the Universe could be subluminal and unattenuated over cosmological distances. Cosmic ray interactions can also be modified and lead to detectable fingerprints in the energy spectrum and mass composition observed on Earth. The data collected at the Pierre Auger Observatory are therefore possibly sensitive to both the electromagnetic and hadronic sectors of LIV. In this article, we explore these two sectors by comparing the energy spectrum and the composition of cosmic rays and the upper limits on the photon flux from the Pierre Auger Observatory with simulations including LIV. Constraints on LIV parameters depend strongly on the mass composition of cosmic rays at the highest energies. For the electromagnetic sector, while no constraints can be obtained in the absence of protons beyond $10^{19}$ eV, we obtain $\delta_{\gamma,0} > -10^{-21}$, $\delta_{\gamma,1} > -10^{-40}$ eV$^{-1}$ and $\delta_{\gamma,2} > -10^{-58}$ eV$^{-2}$ in the case of a subdominant proton component up to $10^{20}$ eV. For the hadronic sector, we study the best description of the data as a function of LIV coefficients and we derive constraints in the hadronic sector such as $\delta_{\mathrm{had},0} < 10^{-19}$, $\delta_{\mathrm{had},1} < 10^{-38}$ eV$^{-1}$ and $\delta_{\mathrm{had},2}< 10^{-57}$ eV$^{-2}$ at 5$\sigma$ CL.

The astrophysics of cosmic dawn, when star formation commenced in the first collapsed objects, is predicted to be revealed as spectral and spatial signatures in the cosmic radio background at long wavelengths. The sky-averaged redshifted 21-cm absorption line of neutral hydrogen is a probe of cosmic dawn. The line profile is determined by the evolving thermal state of the gas, radiation background, Lyman-$\alpha$ radiation from stars scattering off cold primordial gas and the relative populations of the hyperfine spin levels in neutral hydrogen atoms. We report a radiometer measurement of the spectrum of the radio sky in the 55--85~MHz band, which shows that the profile found by Bowman et al. in data taken with the Experiment to Detect the Global Epoch of Reionization Signature (EDGES) low-band instrument is not of astrophysical origin; their best-fitting profile is rejected with 95.3\% confidence. The profile was interpreted to be a signature of cosmic dawn; however, its amplitude was substantially higher than that predicted by standard cosmological models. Explanations for the amplitude of the profile included non-standard cosmology, additional mechanisms for cooling the baryons, perhaps via interactions with millicharged dark matter and an excess radio background at redshifts beyond 17. Our non-detection bears out earlier concerns and suggests that the profile found by Bowman et al. is not evidence for new astrophysics or non-standard cosmology.

A spatio-temporal analysis of IRIS spectra of MgII, CII, and SiIV ions allows us to study the dynamics and the stratification of the flare atmosphere along the line of sight during the magnetic reconnection phase at the jet base. Strong asymmetric MgII and CII line profiles with extended blue wings observed at the reconnection site are interpreted by the presence of two chromospheric temperature clouds: one explosive cloud with blueshifts at 290 km/s and one cloud with smaller Doppler shift (around 36 km/s). Simultaneously at the same location a mini flare was observed with strong emission in multi temperatures (AIA), in several spectral IRIS lines (e.g. Oiv and Siiv, Mgii), absorption of identified chromospheric lines in Siiv line profile, enhancement of the Balmer continuum and X-ray emission by FERMI/GBM. With the standard thick-target flare model we calculate the energy of non thermal electrons observed by FERMI and compare it to the energy radiated by the Balmer continuum emission. We show that the low energy input by non thermal electrons above 20 keV was still sufficient to produce the excess of Balmer continuum.

Continued observations of the Double Pulsar, PSR J0737-3039A/B, consisting of two radio pulsars (A and B) that orbit each other with a period of 2.45hr in a mildly eccentric (e=0.088) binary system, have led to large improvements in the measurement of relativistic effects in this system. With a 16-yr data span, the results enable precision tests of theories of gravity for strongly self-gravitating bodies and also reveal new relativistic effects that have been expected but are now observed for the first time. These include effects of light propagation in strong gravitational fields which are currently not testable by any other method. We observe retardation and aberrational light-bending that allow determination of the pulsar's spin direction. In total, we have detected seven post-Keplerian (PK) parameters, more than for any other binary pulsar. For some of these effects, the measurement precision is so high that for the first time we have to take higher-order contributions into account. These include contributions of A's effective mass loss (due to spin-down) to the observed orbital period decay, a relativistic deformation of the orbit, and effects of the equation of state of super-dense matter on the observed PK parameters via relativistic spin-orbit coupling. We discuss the implications of our findings, including those for the moment of inertia of neutron stars. We present the currently most precise test of general relativity's (GR's) quadrupolar description of gravitational waves, validating GR's prediction at a level of $1.3 \times 10^{-4}$ (95% conf.). We demonstrate the utility of the Double Pulsar for tests of alternative theories by focusing on two specific examples and discuss some implications for studies of the interstellar medium and models for the formation of the Double Pulsar. Finally, we provide context to other types of related experiments and prospects for the future.

Here we study the MonR2 star forming region, which has a rich network of filaments joining in a star cluster forming hub, aiming at understanding the hub structure and to examine the mass fraction residing in the hub and in the filaments, which is a key factor that influences massive star formation. We conducted a multi-scale, multi-component analysis of the Herschel column density maps (resolution of 18.2" or $\sim$0.07 pc at 830 pc) of the region using a newly developed algorithm "getsf" to identify the structural components, namely, extended cloud, filaments, and sources. We find that cascades of lower column density filaments coalesce to form higher density filaments eventually merging inside the hub (0.8 pc radius). As opposed to the previous view of the hub as a massive clump with $\sim$1 pc radius, we find it to be a network of short high-density filaments. The total mass reservoir in the MonR2 HFS (5 pc $\times$ 5 pc) is split between filaments (54%), extended cloud (37%) and sources (9%). The M/L of filaments increase from $\sim$ 10 Msun/pc at 1.5pc from the hub to $\sim$ 100 Msun/pc at its centre, while the number of filaments per annulus of 0.2pc width decreases from 20 to 2 in the same range. The observed radial column density structure of the HFS (filament component only) displays a power-law dependence of $N_{\mathrm{H}_2} \propto r^{-2.17}$ up to a radius of $\sim$2.5 pc from the central hub, resembling a global collapse of the HFS. We present a scenario where the HFS can be supported by magnetic fields which interact, merge and reorganize themselves as the filaments coalesce. In the new view of the hub as a network of high-density filaments, we suggest that only the stars located in the network can benefit from the longitudinal flows of gas to become massive, which may explain the reason for the formation of many low-mass stars in cluster centres.

Extragalactic radio sources appear under different morphologies, the most frequent ones are classified as Fanaroff-Riley type I (FR I), typically with lower luminosities, and Fanaroff-Riley type II, (FR II), typically more luminous. This simple classification, however, has many exceptions that we intend to investigate. Following previous analyses in the three-dimensional Hydrodynamic and Magneto-Hydrodynamic limits, we extend the numerical investigation to the Relativistic Magneto-Hydrodynamic regime, to include sources whose jet kinetic power sets in the range that separates FR Is from FR IIs. We consider weakly and mildly relativistic, underdense, supersonic jets that propagate in a stratified medium. In the model, the ambient temperature increases with distance from the jet origin maintaining constant pressure. We present three cases with low, high and intermediate kinetic luminosity that evolve into different morphologies. We find that the resulting morphology can be highly time dependent and that, apart from the jet power, the jet-to-ambient density ratio and the magnetization parameter play a crucial role in the jet evolution as well.

Models of stellar population synthesis (SPS) are the fundamental tool that relates the physical properties of a galaxy to its spectral energy distribution (SED). In this paper, we present DSPS: a python package for stellar population synthesis. All of the functionality in DSPS is implemented natively in the JAX library for automatic differentiation, and so our predictions for galaxy photometry are fully differentiable, and directly inherit the performance benefits of JAX, including portability onto GPUs. DSPS also implements several novel features, such as i) a flexible empirical model for stellar metallicity that incorporates correlations with stellar age, and ii) support for the diffstar model that provides a physically-motivated connection between the star formation history of a galaxy (SFH) and the mass assembly of its underlying dark matter halo. We detail a set of theoretical techniques for using autodiff to calculate gradients of predictions for galaxy SEDs with respect to SPS parameters that control a range of physical effects, including SFH, stellar metallicity, nebular emission, and dust attenuation. When forward modeling the colors of a synthetic galaxy population, we find that DSPS can provide a factor of 20 speedup over standard SPS codes on a CPU, and a factor of over 1000 on a modern GPU. When coupled with gradient-based techniques for optimization and inference, DSPS makes it practical to conduct expansive likelihood analyses of simulation-based models of the galaxy--halo connection that fully forward model galaxy spectra and photometry.

Most general-purpose classification methods, such as support-vector machine (SVM) and random forest (RF), fail to account for an unusual characteristic of astronomical data: known measurement error uncertainties. In astronomical data, this information is often given in the data but discarded because popular machine learning classifiers cannot incorporate it. We propose a simulation-based approach that incorporates heteroscedastic measurement error into any existing classification method to better quantify uncertainty in classification. The proposed method first simulates perturbed realizations of the data from a Bayesian posterior predictive distribution of a Gaussian measurement error model. Then, a chosen classifier is fit to each simulation. The variation across the simulations naturally reflects the uncertainty propagated from the measurement errors in both labeled and unlabeled data sets. We demonstrate the use of this approach via two numerical studies. The first is a thorough simulation study applying the proposed procedure to SVM and RF, which are well-known hard and soft classifiers, respectively. The second study is a realistic classification problem of identifying high-$z$ $(2.9 \leq z \leq 5.1)$ quasar candidates from photometric data. The data were obtained from merged catalogs of the Sloan Digital Sky Survey, the $Spitzer$ IRAC Equatorial Survey, and the $Spitzer$-HETDEX Exploratory Large-Area Survey. The proposed approach reveals that out of 11,847 high-$z$ quasar candidates identified by a random forest without incorporating measurement error, 3,146 are potential misclassifications. Additionally, out of ${\sim}1.85$ million objects not identified as high-$z$ quasars without measurement error, 936 can be considered candidates when measurement error is taken into account.

This paper presents an integrated modeling software to analyze the PSF of wide-field telescopes affected by misalignments. Even relatively small misalignments in the optical system of a telescope can significantly deteriorate the image quality by introducing large aberrations. In particular, wide-field telescopes are critically affected by these errors, insomuch that usually a closed-loop active optics system is adopted for a continuous correction, rather than for sporadic alignment procedures. Typically, a ray-tracing software such as Zemax OpticStudio is employed to accurately analyze the system during the optical design. However, an analytical model of the optical system is preferable when the PSF of the telescope must be reconstructed quickly for algorithmic purposes. Here the analytical model is derived through a hybrid approach and developed in a custom software package, designed to be general and flexible in order to be tailored to different optical configurations. First, leveraging on the Zemax OpticStudio API, the ray-tracing software is integrated into a Matlab pipeline. This allows to perform a statistical analysis by automatically simulating the system response in a variety of misaligned working conditions. Then, the resulting dataset is employed to populate a database of parameters describing the model.

Though cluster-summed luminosities have served as mass proxies, cluster-summed colors have received less attention. Since galaxy colors have given useful insights into dust content and specific star formation rates, this research investigates possible correlations between cluster-summed colors and various observable and intrinsic halo properties for clusters in subsamples of TNG, SDSS, and Buzzard. Cluster color--magnitude space shows a peak towards the red and bright corner, drawn there by bright red galaxies. Summing colors across a cluster reduces the scatter in color spaces, since magnitude summing acts somewhat like a weighted average. The correlation between these summed colors were $(73 \pm 24)\%$ across all three datasets. Summed colors and cluster properties typically had low correlations but ranged up to $\sim 40\%$. The correlation between color and mass didn't change significantly with richness threshold for TNG and Buzzard, but for SDSS the correlation decreased dramatically with increasing richness, passing from positive correlation to negative correlation near a richness threshold of ten. We also looked at mass proxy scaling relations with richness or magnitude and measured the reduction in mass scatter once we added cluster colors. The reduction was generally insignificant, but several large reductions in mass scatter occurred under certain circumstances: high-richness Buzzard mass--magnitude relation saw a reduction of $(19 \pm 28)\%$ while low-richness SDSS saw similar order reductions of $(16 \pm 8)\%$ and $(14 \pm 8)\%$ for the mass--richness and mass--magnitude relations respectively. This first look at summed cluster color shows potential in aiding mass proxies under certain circumstances, but more deliberate and thorough investigations are needed to better characterize and make use of cluster-summed colors.

We present the fourth Open Gravitational-wave Catalog (4-OGC) of binary neutron star (BNS), binary black hole (BBH) and neutron star-black hole (NSBH) mergers. The catalog includes observations from 2015-2020 covering the first through third observing runs (O1, O2, O3a, O3b) of Advanced LIGO and Advanced Virgo. The updated catalog includes 7 BBH mergers which were not previously reported with high significance during O3b for a total of 94 observations: 90 BBHs, 2 NSBHs, and 2 BNSs. The most confident new detection, GW200318_191337, has component masses $49.1^{+16.4}_{-12.0}~M_{\odot}$ and $31.6^{+12.0}_{-11.6}~M_{\odot}$; its redshift of $0.84^{+0.4}_{-0.35}$ ($90\%$ credible interval) may make it the most distant merger so far. We provide reference parameter estimates for each of these sources using an up-to-date model accounting for instrumental calibration uncertainty. The corresponding data release also includes our full set of sub-threshold candidates.

We study a Lagrangian with a cubic Galileon term and a standard scalar-field kinetic contribution with two exponential potentials. In this model the Galileon field generates scaling solutions in which the density of the scalar field $\phi$ scales in the same manner as the matter density at early-time. These solutions are of high interest because the scalar field can then be compatible with the energy scale of particle physics and can alleviate the coincidence problem. The phenomenology of linear perturbations is thoroughly discussed, including all the relevant effects on the observables. Additionally, we use cosmic microwave background temperature-temperature and lensing power spectra by Planck 2018, the baryon acoustic oscillations measurements from the 6dF galaxy survey and SDSS and supernovae type Ia data from Pantheon in order to place constraints on the parameters of the model. We find that despite its interesting phenomenology, the model we investigate does not produce a better fit to data with respect to $\Lambda$CDM, and it does not seem to be able to ease the tension between high and low redshift data.

The nucleosynthesis in classical novae, in particular that of radioactive isotopes, is directly measurable by its $\gamma$-ray signature. Despite decades of observations, MeV $\gamma$-rays from novae have never been detected -- neither individually at the time of the explosion, nor as a result of radioactive decay, nor the diffuse Galactic emission from the nova population. Thanks to recent developments in modeling of instrumental background for MeV telescopes such as INTEGRAL/SPI and Fermi/GBM, the prospects to finally detect these elusive transients are greatly enhanced. This demands for updated and refined models of $\gamma$-ray spectra and light curves of classical novae. In this work, we develop numerical models of nova explosions using sub- and near-Chandrasekhar CO white dwarfs as the progenitor. We study the parameter dependence of the explosions, their thermodynamics and energetics, as well as their chemical abundance patterns. We use a Monte-Carlo radiative transfer code to compute $\gamma$-ray light curves and spectra, with a focus on the early time evolution. We compare our results to previous studies and find that the expected 511-keV-line flash at the time of the explosion is heavily suppressed, showing a maximum flux of only $10^{-9}\,{\rm ph}\,$cm$^{-2}\,$s$^{-1}$ and thus making it at least one million times fainter than estimated before. This finding would render it impossible for current MeV instruments to detect novae within the first day after the outburst. Nevertheless, our time-resolved spectra can be used for retrospective analyses of archival data, thereby improving the sensitivity of the instruments.

Oscillations in the frequency profile of the stochastic gravitational wave background are a characteristic prediction of small-scale features during inflation. In this paper we present a first investigation of the detection prospects of such oscillations with the upcoming space-based gravitational wave observatory LISA. As a proof of principle, we show for a selection of feature signals that the oscillations can be reconstructed with LISA, employing a method based on principal component analysis. We then perform a Fisher forecast for the parameters describing the oscillatory signal. For a sharp feature we distinguish between the contributions to the stochastic gravitational wave background induced during inflation and in the post-inflationary period, which peak at different frequencies. We find that for the latter case the amplitude of the oscillation is expected to be measurable with $< 10\%$ accuracy if the corresponding peak satisfies $h^2 \Omega_\textrm{GW} \gtrsim 10^{-12}$-$10^{-11}$, while for inflationary-era gravitational waves a detection of the oscillations requires a higher peak amplitude of $h^2 \Omega_\textrm{GW}$, as the oscillations only appear on the UV tail of the spectrum. For a resonant feature the detection prospects with LISA are maximised if the frequency of the oscillation falls into the range $\omega_\textrm{log} = 4$ to $10$. Our results confirm that oscillations in the frequency profile of the stochastic gravitational wave background are a worthwhile target for future detection efforts and offer a key for experimentally testing inflation at small scales.

High-precision measurements of the pulsar dispersion measure (DM) are possible using telescopes with low-frequency wideband receivers. We present an initial study of the application of the wideband timing technique, which can simultaneously measure the pulsar times of arrival (ToAs) and DMs, for a set of five pulsars observed with the upgraded Giant Metrewave Radio Telescope (uGMRT) as part of the Indian Pulsar Timing Array (InPTA) campaign. We have used the observations with the 300-500 MHz band of the uGMRT for this purpose. We obtain high precision in DM measurements with precisions of the order 10^{-6}cm^{-3}pc. The ToAs obtained have sub-{\mu}s precision and the root-mean-square of the post-fit ToA residuals are in the sub-{\mu}s range. We find that the uncertainties in the DMs and ToAs obtained with this wideband technique, applied to low-frequency data, are consistent with the results obtained with traditional pulsar timing techniques and comparable to high-frequency results from other PTAs. This work opens up an interesting possibility of using low-frequency wideband observations for precision pulsar timing and gravitational wave detection with similar precision as high-frequency observations used conventionally.

We show for the first time how to conduct a direct search for dark matter using Gaia observations. Its public astrometric data may contain the signals of mirror stars, exotic compact objects made of atomic dark matter with a tiny kinetic mixing between the dark and SM photon. Mirror stars capture small amounts of interstellar material in their cores, leading to characteristic optical/IR and X-ray emissions. We develop the detailed pipeline for conducting a mirror star search using data from Gaia and other stellar catalogues, and demonstrate our methodology by conducting a search for toy mirror stars with a simplified calculation of their optical/IR emissions over a wide range of mirror star and hidden sector parameters. We also obtain projected exclusion bounds on the abundance and properties of mirror stars if no candidates are found, demonstrating that Gaia is a new and uniquely powerful probe of atomic dark matter. Our study provides the blueprint for a realistic mirror star search that includes a more complete treatment of the captured interstellar gas in the future.

Ultralight scalar fields such as axions can form clouds around rotating black holes (BHs) by the superradiant instability. It is important to consider the evolution of clouds associated with BH binaries for the detectability of the presence of clouds through gravitational wave signals and observations of the mass and spin parameters of BHs. We re-examine the impact on the axion cloud due to the tidal perturbation from the companion in a binary system taking into account the following points. First, we study the influence of higher-multipole moments. Second, we consider the backreaction due to the angular momentum transfer between the cloud and the orbital motion. This angular momentum transfer further causes the backreaction to the hyperfine split through the change in geometry. Finally, we calculate the particle number flux to the infinity induced by the tidal interaction. As a result, we found that the scalar field is not reabsorbed by the BH. Instead, the scalar particles are radiated away to evaporate during the inspiral, irrespectively of the direction of the orbital motion, for almost equal mass binaries.

We study the parametrized post-Newtonian (PPN) limit of higher-derivative-torsion Modified Teleparallel Gravity. We start from the covariant formulation of modified Teleparallel Gravity by restoring the spin connection of the theory. Then, we perform the post-Newtonian expansion of the tetrad field around the Minkowski background and find the perturbed field equations. We compute the PPN metric for the higher-order Teleparallel Gravity theories which allows us to show that at the post-Newtonian limit this more general class of theories are fully conservative and indistinguishable from General Relativity . In this way, we extend the results that were already found for $F(T)$ gravity in previous works. Furthermore, our calculations reveal the importance of considering a second post-Newtonian (2PN) order approximation or a parametrized post-Newtonian cosmology (PPNC) framework where additional perturbative modes coming from general modifications of Teleparallel Gravity could lead to new observable imprints.

A nonlocal gravity model (2.1) was introduced and considered recently [49], and two exact cosmological solutions in flat space were presented. The first solution is related to some radiation effects generated by nonlocal dynamics on dark energy background, while the second one is a nonsingular time symmetric bounce. In the present paper we investigate other possible exact cosmological solutions and find some the new ones in nonflat space. Used nonlocal gravity dynamics can change background topology. To solve the corresponding eqations of motion, we first look for a solution of the eigenvalue problem $\Box (R -4\Lambda) = q\ (R - 4\Lambda) .$ We also discuss possible extension of this model with nonlocal operator symmetric under $\Box \longleftrightarrow \Box^{-1}$ and its connection with another interesting nonlocal gravity model.

Our Universe is expected to finally approach a de Sitter universe whose horizon is considered to be in thermal equilibrium. In the present article, both the energy stored on the horizon and its thermodynamic fluctuations are examined through the holographic equipartition law. First, it is confirmed that a flat Friedmann--Robertson--Walker universe approaches a de Sitter universe, using a cosmological model close to lambda cold dark matter ($\Lambda$CDM) models. Then, based on the holographic equipartition law, the energy density of the Hubble volume is calculated from the energy on the Hubble horizon of a de Sitter universe. The energy density for a de Sitter universe is constant and the order of the energy density is consistent with the order of that for the observed cosmological constant. Second, thermodynamic fluctuations of energy on the horizon are examined, assuming stable fluctuations around thermal equilibrium states. A standard formulation of the fluctuations for a canonical ensemble is applied to the Hubble horizon of a de Sitter universe. The thermodynamic fluctuations of the energy are found to be a universal constant corresponding to the Planck energy, regardless of the Hubble parameter. In contrast, the relative fluctuations of the energy can be characterized by the ratio of the one-degree-of-freedom energy to the Planck energy. At the present time, the order of the relative fluctuations should be within the range of a discrepancy derived from a discussion of the cosmological constant problem, namely a range approximately from $10^{-60}$ to $10^{-123}$. The present results may imply that the energy stored on the Hubble horizon is related to a kind of effective dark energy, whereas the energy that can be `maximumly' stored on the horizon may behave as if it were a kind of effective vacuum-like energy in an extended holographic equipartition law.

The spherical modes of gravitational waves (GWs) have become a major focus of resent detection campaigns due to the additional information they can provide about different properties of the source. However, GW detection is restricted to only detecting one ray and hence it is not obvious how we can extract information about angular properties. In this paper we introduce a new gauge which makes visible GW detection does not only contain information on the second time derivative but also on the angular derivatives of the GW. In particular, we show that the angular derivatives are of the same order as the time derivatives of the wave thus allowing us to constrain the spherical modes.

High-order spatial derivatives are of crucial importance for constructing the low energy effective action of a Lorentz or parity violating theory of quantum gravity. One example is the Ho\v{r}ava-Lifshitz gravity, in one has to consider at least the sixth-order spatial derivatives in the gravitational action, in order to make the theory power-counting renormalizable. In this paper, we consider the Lorentz and parity violating effects on the propagation of GWs due to the fifth and sixth-order spatial derivatives respectively. For this purpose we calculate the corresponding Lorentz and parity violating waveforms of GWs produced by the coalescence of compact binaries. By using these modified waveforms, we perform the full Bayesian inference with the help of the open source software \texttt{Bilby} on the selected GW events of binary black hole (BBH) and binary neutron stars (BNS) merges in the LIGO-Virgo catalogs GWTC-1 and GWTC-2. Overall we do not find any significant evidence of Lorentz and parity violation due to the fifth and sixth-order spatial derivatives and thus place lower bounds on the energy scales $M_{\rm LV} > 2.4 \times 10^{-16} \; {\rm GeV}$ for Lorentz violation and $M_{\rm PV} > 1.0 \times 10^{-14} \; {\rm GeV}$ for parity violation at 90\% confidence level. Both constraints represent the first constraints on the fifth- and sixth-order spatial derivative terms respectively in the framework of spatial covariant gravity by using the observational data of GWs.

Motivated by the current strong constraints on the spin-independent dark matter (DM)-nucleus scattering, we investigate the spin-dependent (SD) interactions of the light Majorana DM with the nucleus mediated by an axial-vector boson. Due to the small nucleus recoil energy, the ionization signals have now been used to probe the light dark matter particles in direct detection experiments. With the existing ionization data, we derive the exclusion limits on the SD DM-nucleus scattering through Migdal effect in the MeV-GeV DM mass range. It is found that the lower limit of the DM mass can reach about several MeVs. Due to the momentum transfer correction induced by the light mediator, the bounds on the SD DM-nucleus scattering cross sections can be weakened in comparison with the heavy mediator.

Coherent WaveBurst is a generic, multi-detector gravitational wave burst search based on the excess power approach. The coherent WaveBurst algorithm currently employed in the all-sky short-duration gravitational wave burst search uses a conditional approach on a selected attributes in the multi-dimensional event attribute space to distinguish between noisy event from that of astrophysical origin. We have been developing a supervised machine learning approach based on the Gaussian mixture modeling to model the attribute space for signals as well as noise events to enhance the probability of burst detection. We further extend the GMM approach to the all-sky short-duration coherent WaveBurst search as a post-processing step on events from the first half of the third observing run (O3a). We show an improvement in sensitivity to generic gravitational wave burst signal morphologies as well as the astrophysical source such as core-collapse supernova models due to the application of our Gaussian mixture model approach to coherent WaveBurst triggers. The Gaussian mixture model method recovers the gravitational wave signals from massive compact binary coalescences identified by coherent WaveBurst targeted for binary black holes in GWTC-2, with better significance than the all-sky coherent WaveBurst search. No additional significant gravitational wave bursts are observed.

In this work we review the issue of imposing the conservation of the energy-momentum tensor as a necessary condition to recover the equivalence between the unimodular gravity and General Relativity (GR) equipped with a cosmological constant. This procedure is usually interpreted as an {\it ad hoc} imposition on the unimodular theory's structure. Whereas the consequences of avoiding the conservation of the total energy-momentum tensor has been already introduced in the literature, it has been not widely explored so far. We study an expanding universe sourced by a single effective perfect fluid such that the null divergence of its energy-momentum tensor is not imposed. As we shall show, in this scheme, the unimodular theory has its own conservation equation obtained from the Bianchi identities. We explore the evolution of the homogeneous and isotropic expanding background and show that a viable cosmological scenario exists. Also, we consider scalar perturbations with particular attention given to the gauge issue. We show that contrary to the traditional unimodular theory where the synchronous and longitudinal (newtonian) gauge for cosmological perturbations are not permitted, if the conservation of the energy-momentum is relaxed the scalar perturbations in the synchronous condition survive and present a growing mode behavior. We study therefore a new cosmological scenario in which the dynamics of the universe transits from the radiative phase directly to a accelerated one but allowing thus for structure formation.

Unimodular gravity is characterized by an extra condition with respect to General Relativity: the determinant of the metric is constant. This extra condition leads to a more restricted class of invariance by coordinate transformation. Even though, if the conservation of the energy-momentum tensor is imposed in unimodular gravity, the General Relativity theory is recovered with an additional integration constant which is associated to the cosmological term $\Lambda$. However, if the energy-momentum tensor does not conserve separately, a new geometric structure appears with potentially observational signatures. In this text, we consider the evolution of gravitational waves in the nonconservative unimodular gravity, showing how it differs from the usual signatures in the standard model.

Detectors based on liquid argon (LAr) often require surfaces that can shift vacuum ultraviolet (VUV) light and reflect the visible shifted light. For the LAr instrumentation of the LEGEND-200 neutrinoless double beta decay experiment, several square meters of wavelength-shifting reflectors (WLSR) were prepared: the reflector Tetratex (TTX) was in-situ evaporated with the wavelength shifter tetraphenyl butadiene (TPB). For even larger detectors, TPB evaporation will be more challenging and plastic films of polyethylene naphthalate (PEN) are considered as an option to ease scalability. In this work, we first characterized the absorption (and reflectivity) of PEN, TPB (and TTX) films in response to visible light. We then measured TPB and PEN coupled to TTX in a LAr setup equipped with a VUV sensitive photomultiplier tube. The effective light yield in the setup was first measured using an absorbing reference sample, and the VUV reflectivity of TTX quantified. The characterization and simulation of the setup along with the measurements and modelling of the optical parameters of TPB, PEN and TTX allowed to estimate the quantum efficiency (QE) of TPB and PEN in LAr (at 87K) for the first time: these were found to be above 67% and 49%, respectively (at 90% CL). These results provide relevant input for the optical simulations of experiments that use TPB in LAr, such as LEGEND-200, and for experiments that plan to use TPB or PEN to shift VUV scintillation light.

We formulate power-law holographic dark energy, which is a modified holographic dark energy model based on the extended entropy relation arising from the consideration of state mixing between the ground and the excited ones in the calculation of the entanglement entropy. We construct two cases of the scenario, imposing the usual future event horizon choice, as well as the Hubble one. Thus, the former model is a one-parameter extension of standard holographic dark energy, recovering it in the limit where power-law extended entropy recovers Bekenstein-Hawking one, while the latter belongs to the class of running vacuum models, a feature that may reveal the connection between holography and the renormalization group running. For both models we extract the differential equation that determines the evolution of the dark-energy density parameter and we provide the expression for the corresponding equation-of-state parameter. We find that the scenario can describe the sequence of epochs in the Universe evolution, namely the domination of matter followed by the domination of dark energy. Moreover, the dark-energy equation of state presents a rich behavior, lying in the quintessence regime or passing into the phantom one too, depending on the values of the two model parameters, a behavior that is richer than the one of standard holographic dark energy.

We introduce the Merger Entropy Index ($\mathcal{I}_\mathrm{BBH}$), a new parameter to measure the efficiency of entropy transfer for any generic binary black hole merger in General Relativity. We find that $\mathcal{I}_\mathrm{BBH}$ is bounded between an asymptotic maximum and minimum. For the observed population of mergers detected by LIGO and Virgo, we find that $\mathcal{I}_\mathrm{BBH}$ is $\lesssim30\%$ of its theoretical maximum. By imposing the thermodynamical consistency between the pre- and post-merger states through $\mathcal{I}_\mathrm{BBH}$, we showcase BRAHMA -- a novel framework to infer the properties and astrophysical implications of gravitational-wave detections. For GW190521 -- the heaviest confirmed binary black hole merger observed so far -- our framework rules out high mass-ratio, negative effective inspiral spin, and electromagnetic counterpart claims. Furthermore, our analysis provides an independent confirmation that GW190521 belongs to a separate population.

The ever-increasing number of detections of gravitational waves (GWs) from compact binaries by the Advanced LIGO and Advanced Virgo detectors allows us to perform ever-more sensitive tests of general relativity (GR) in the dynamical and strong-field regime of gravity. We perform a suite of tests of GR using the compact binary signals observed during the second half of the third observing run of those detectors. We restrict our analysis to the 15 confident signals that have false alarm rates $\leq 10^{-3}\, {\rm yr}^{-1}$. In addition to signals consistent with binary black hole (BH) mergers, the new events include GW200115_042309, a signal consistent with a neutron star--BH merger. We find the residual power, after subtracting the best fit waveform from the data for each event, to be consistent with the detector noise. Additionally, we find all the post-Newtonian deformation coefficients to be consistent with the predictions from GR, with an improvement by a factor of ~2 in the -1PN parameter. We also find that the spin-induced quadrupole moments of the binary BH constituents are consistent with those of Kerr BHs in GR. We find no evidence for dispersion of GWs, non-GR modes of polarization, or post-merger echoes in the events that were analyzed. We update the bound on the mass of the graviton, at 90% credibility, to $m_g \leq 1.27 \times 10^{-23} \mathrm{eV}/c^2$. The final mass and final spin as inferred from the pre-merger and post-merger parts of the waveform are consistent with each other. The studies of the properties of the remnant BHs, including deviations of the quasi-normal mode frequencies and damping times, show consistency with the predictions of GR. In addition to considering signals individually, we also combine results from the catalog of GW signals to calculate more precise population constraints. We find no evidence in support of physics beyond GR.