Papers for Friday, Feb 12 2021

Recent observational constraints indicate that primordial black holes (PBHs) with the mass scale $\sim 10^{-12}M_{\odot}$ can explain most of dark matter in the Universe. To produce this kind of PBHs, we need an enhance in the primordial scalar curvature perturbations to the order of ${\mathcal{O}(10^{-2})}$ at the scale $ k \sim 10^{12}~\rm Mpc^{-1}$. Here, we investigate the production of PBHs and induced gravitational waves (GWs) in the framework of Galileon inflation. To this aim, we consider the Galileon term $G(X)=X/M^3$ as well as the $\alpha$-attractor base potential modified by a small local Gaussian bump. We solve numerically the Mukhanov-Sasaki equation to obtain the primordial scalar power spectrum. In addition, we estimate the PBHs abundance $f_{\text{PBH}}^{\text{peak}}$ as well as the energy density parameter $\Omega_{\rm GW,0}$ of induced GWs. Interestingly enough is that for a special set of model parameters, we estimate the mass scale and the abundance of PBHs as $\sim{\cal O}(10^{-12})M_{\odot}$ and $f_{\text{PBH}}^{\text{peak}}=0.92$, respectively. This confirms that the mechanism of PBHs production in Galileon inflation can justify most of dark matter. Furthermore, we evaluate the GWs energy density parameter and conclude that it behaves like a power-law function $\Omega_{\rm GW}\sim (f/f_c)^n$ where in the infrared limit $f\ll f_{c}$, the power index reads $n=3-2/\ln(f_c/f)$.

Luminous Red Novae (LRNe) are astrophysical transients associated with the partial ejection of a binary system's common envelope (CE) shortly before its merger. Here we present the results of our photometric and spectroscopic follow-up campaign of AT2018bwo (DLT18x), a LRN discovered in NGC45, and investigate its progenitor system using binary stellar-evolution models. The transient reached a peak magnitude of $M_r=-10.97\pm0.11$ and maintained this brightness during its optical plateau of $t_p = 41\pm5$days. During this phase, it showed a rather stable photospheric temperature of ~3300K and a luminosity of ~$10^{40}$erg/s. The photosphere of AT2018bwo at early times appeared larger and cooler than other similar LRNe, likely due to an extended mass-loss episode before the merger. Towards the end of the plateau, optical spectra showed a reddened continuum with strong molecular absorption bands. The reprocessed emission by the cooling dust was also detected in the mid-infrared bands ~1.5 years after the outburst. Archival Spitzer and Hubble Space Telescope data taken 10-14 years before the transient event suggest a progenitor star with $T_{prog}\sim 6500$K, $R_{prog}\sim 100R_{\odot}$ and $L_{prog}\sim 2\times10^4L_{\odot}$, and an upper limit for optically thin warm (1000 K) dust mass of $M_d<10^{-6}M_{\odot}$. Using stellar binary-evolution models, we determined the properties of binary systems consistent with the progenitor parameter space. For AT2018bwo, we infer a primary mass of 12-16 $M_{\odot}$, which is 9-45% larger than the ~11$M_{\odot}$ obtained using single-star evolution models. The system, consistent with a yellow-supergiant primary, was likely in a stable mass-transfer regime with -2.4<log ($\dot{M}/M_{\odot}$/yr)<-1.2 a decade before the main instability occurred. During the dynamical merger, the system would have ejected 0.15-0.5$M_{\odot}$ with a velocity of ~500 km/s.

Measuring the star-forming properties of AGN hosts is key to our understanding of galaxy formation and evolution. However, this topic remains debated, partly due to the difficulties in separating the infrared (i.e. 1--1000 $\mu$m) emission into AGN and star-forming components. Taking advantage of archival far-infrared data from Herschel, we present a new set of AGN and galaxy infrared templates, and introduce the spectral energy distribution fitting code IRAGNSEP. Both can be used to measure infrared host galaxy properties, free of AGN contamination. To build these, we used a sample of 100 local ($z$ < 0.3), low-to-high luminosity AGNs (i.e. $L_{\rm bol}~\sim~10^{42--46}~\rm erg~s^{-1}$), selected from the 105-month Swift - BAT X-ray survey, which have archival Spitzer - IRS spectra and Herschel photometry. We first built a set of seven galaxy templates using a sample of 55 star-forming galaxies selected via infrared diagnostics. Using these templates, combined with a flexible model for the AGN contribution, we extracted the intrinsic infrared emission of our AGN sample. We further demonstrate that we can reduce the diversity in the intrinsic shapes of AGN spectral energy distributions down to a set of three AGN templates, of which two represent AGN continuum, and one represents silicate emission. Our results indicate that, on average, the contribution of AGNs to the far-infrared ($\lambda~\gtrsim$ 50 $\mu$m) is not as high as suggested by some recent work. We further show that the need for two infrared AGN continuum templates could be related to nuclear obscuration, where one of our templates appears dominated by the emission of the extended polar dust.

Cosmological phase transitions in the primordial universe can produce anisotropic stochastic gravitational wave backgrounds (GWB), similar to the cosmic microwave background (CMB). For adiabatic perturbations, the fluctuations in GWB follow those in the CMB, but if primordial fluctuations carry an isocurvature component, this need no longer be true. It is shown that in non-minimal inflationary and reheating settings, primordial isocurvature can survive in GWB and exhibit significant non-Gaussianity (NG) in contrast to the CMB, while obeying current observational bounds. While probing such NG GWB is at best a marginal possibility at LISA, there is much greater scope at future proposed detectors such as DECIGO and BBO. It is even possible that the first observations of inflation-era NG could be made with gravitational wave detectors as opposed to the CMB or Large-Scale Structure surveys.

1ES 1927+654 is a nearby active galactic nucleus (AGN) which underwent a changing-look event in early 2018, developing prominent broad Balmer lines which were absent in previous observations. We have followed up this object in the X-rays with an ongoing campaign that started in May 2018, and that includes 265 NICER (for a total of 678ks) and 14 Swift/XRT (26ks) observations, as well as three simultaneous XMM-Newton/NuSTAR (158/169 ks) exposures. In the X-rays, 1ES 1927+654 shows a behaviour unlike any previously known AGN. The source is extremely variable both in spectral shape and flux, and does not show any correlation between X-ray and UV flux on timescales of hours or weeks/months. After the outburst the power-law component almost completely disappeared, and the source showed an extremely soft continuum dominated by a blackbody component. The temperature of the blackbody increases with the luminosity, going from $kT\sim 80$eV (for a 0.3--2keV luminosity of $L_{0.3-2}\sim 10^{41.5}\rm\,erg\,s^{-1}$) to $\sim 200$eV (for $L_{0.3-2}\sim 10^{44}\rm\,erg\,s^{-1}$). The spectra show evidence of ionized outflows, and of a prominent feature at $\sim 1$keV, which can be reproduced by a broad emission line. The unique characteristics of 1ES 1927+654 in the X-ray band suggest that it belongs to a new type of changing-look AGN. Future X-ray surveys might detect several more objects with similar properties.

Although tidal dissipation in binary stars has been studied for over a century, theoretical predictions have yet to match the observed properties of binary populations. This work quantitatively examines the recent proposal of tidal circularization by resonance locking, where tidal dissipation arises from resonances between the star's natural oscillation frequencies and harmonics of the orbital frequency, and where resonances are `locked' for an extended period of time due to concurrent stellar evolution. We focus on tidal resonances with axi-symmetric gravity-modes, and examine binaries with primary masses from one to two solar masses. We find that orbital evolution via resonance locking occurs primarily during the star's pre-main-sequence phase, with the main-sequence phase contributing negligibly. Resonance locking, ignoring non-linearity, can circularize binaries with peri-centre distances out to $\sim 10$ stellar radii, corresponding to circular periods of $\sim 4-6$ days. However, we find resonantly excited gravity-modes will become nonlinear in stellar cores, which prevents them from reaching their full, linear amplitudes. We estimate that such a `saturated resonance lock' reduces the circularization period by about a third, but resonance locking remains much more effective than the cumulative actions of equilibrium tides. In a companion paper, we examine recent binary data to compare against theory.

There is a large consensus that gas in high-$z$ galaxies is highly turbulent, because of a combination of stellar feedback processes and gravitational instabilities driven by mergers and gas accretion. In this paper, we present the analysis of a sample of five Dusty Star Forming Galaxies (DSFGs) at $4 \lesssim z\lesssim 5$. Taking advantage of the magnifying power of strong gravitational lensing, we quantified their kinematic and dynamical properties from ALMA observations of their [CII] emission line. We combined the dynamical measurements obtained for these galaxies with those obtained from previous studies to build the largest sample of $z \sim 4.5$ galaxies with high-quality data and sub-kpc spatial resolutions, so far. We found that all galaxies in the sample are dynamically cold, with rotation-to-random motion ratios, $V/\sigma$, between 7 to 15. The relation between their velocity dispersions and their star-formation rates indicates that stellar feedback is sufficient to sustain the turbulence within these galaxies and no further mechanisms are needed. In addition, we performed a rotation curve decomposition to infer the relative contribution of the baryonic (gas, stars) and dark matter components to the total gravitational potentials. This analysis allowed us to compare the structural properties of the studied DSFGs with those of their descendants, the local early type galaxies. In particular, we found that five out of six galaxies of the sample show the dynamical signature of a bulge, indicating that the spheroidal component is already in place at $z \sim 4.5$.

Fast and automated inference of binary-lens, single-source (2L1S) microlensing events with sampling-based Bayesian algorithms (e.g., Markov Chain Monte Carlo; MCMC) is challenged on two fronts: high computational cost of likelihood evaluations with microlensing simulation codes, and a pathological parameter space where the negative-log-likelihood surface can contain a multitude of local minima that are narrow and deep. Analysis of 2L1S events usually involves grid searches over some parameters to locate approximate solutions as a prerequisite to posterior sampling, an expensive process that often requires human-in-the-loop and domain expertise. As the next-generation, space-based microlensing survey with the Roman Space Telescope is expected to yield thousands of binary microlensing events, a new fast and automated method is desirable. Here, we present a likelihood-free inference (LFI) approach named amortized neural posterior estimation, where a neural density estimator (NDE) learns a surrogate posterior $\hat{p}(\theta|x)$ as an observation-parametrized conditional probability distribution, from pre-computed simulations over the full prior space. Trained on 291,012 simulated Roman-like 2L1S simulations, the NDE produces accurate and precise posteriors within seconds for any observation within the prior support without requiring a domain expert in the loop, thus allowing for real-time and automated inference. We show that the NDE also captures expected posterior degeneracies. The NDE posterior could then be refined into the exact posterior with a downstream MCMC sampler with minimal burn-in steps.

We apply BayeSN, our new hierarchical Bayesian model for the SEDs of Type Ia supernovae (SNe Ia), to analyse the $griz$ light curves of 157 nearby SNe Ia ($0.015<z<0.08$) from the public Foundation DR1 dataset. We train a new version of BayeSN, continuous from 0.35--0.95 $\mu$m, which we use to model the properties of SNe Ia in the rest-frame $z$-band, study the properties of dust in their host galaxies, and construct a Hubble diagram of SN Ia distances determined from full $griz$ light curves. Our $griz$ Hubble diagram has a low total RMS of 0.13 mag using BayeSN, compared to 0.16 mag using SALT2. Additionally, we test the consistency of the dust law $R_V$ between low- and high-mass host galaxies by using our model to fit the full time- and wavelength-dependent SEDs of SNe Ia up to moderate reddening (peak apparent $B-V \lesssim 0.3$). Splitting the population at the median host mass, we find $R_V=2.84\pm0.31$ in low-mass hosts, and $R_V=2.58\pm0.23$ in high-mass hosts, both consistent with the global value of $R_V=2.61\pm0.21$ that we estimate for the full sample. For all choices of mass split we consider, $R_V$ is consistent across the step within $\lesssim1.2\sigma$. Modelling population distributions of dust laws in low- and high-mass hosts, we find that both subsamples are highly consistent with the full sample's population mean $\mu(R_V) = 2.70\pm0.25$ with a 95% upper bound on the population $\sigma(R_V) < 0.61$. The $R_V$ population means are consistent within $\lesssim1.2\sigma$. We find that simultaneous fitting of host-mass-dependent dust properties within our hierarchical model does not account for the conventional mass step.

We explore the kinematics of 27 z~6 quasar host galaxies observed in [CII]-158 micron ([CII]) emission with the Atacama Large Millimeter/sub-millimeter Array at a resolution of ~0.25''. We find that nine of the galaxies show disturbed [CII] emission, either due to a close companion galaxy or recent merger. Ten galaxies have smooth velocity gradients consistent with the emission arising from a gaseous disk. The remaining eight quasar host galaxies show no velocity gradient, suggesting that the gas in these systems is dispersion-dominated. All galaxies show high velocity dispersions with a mean of 129+-10 km/s. To provide an estimate of the dynamical mass within twice the half-light radius of the quasar host galaxy, we model the kinematics of the [CII] emission line using our publicly available kinematic fitting code, qubefit. This results in a mean dynamical mass of 5.0+-0.8(+-3.5) x 10^10 Msun. Comparison between the dynamical mass and the mass of the supermassive black hole reveals that the sample falls above the locally derived bulge mass--black hole mass relation at 2.4sigma significance. This result is robust even if we account for the large systematic uncertainties. Using several different estimators for the molecular mass, we estimate a gas mass fraction of >10%, indicating gas makes up a large fraction of the baryonic mass of z~6 quasar host galaxies. Finally, we speculate that the large variety in [CII] kinematics is an indication that gas accretion onto z~6 super massive black holes is not caused by a single precipitating factor.

Numerical experiments are the primary method of studying the evolution of circumbinary disks due to the strong nonlinearities involved. Many circumbinary simulations also require the use of numerical mass sinks: source terms which prevent gas from unphysically accumulating around the simulated point masses by removing gas at a given rate. However, special care must be taken when drawing physical conclusions from such simulations to ensure that results are not biased by numerical artifacts. We demonstrate how the use of improved sink methods reduces some of these potential biases in vertically-integrated simulations of aspect ratio 0.1 accretion disks around binaries with mass ratios between 0.1 and 1. Specifically, we show that sink terms that do not reduce the angular momentum of gas relative to the accreting object: 1) reduce the dependence on the sink rate of physical quantities such as the torque on the binary, distribution of accretion between binary components, and evolution of the binary semi-major axis; 2) reduce the degree to which the sink rate affects the structure of the accretion disks around each binary component; 3) alter the inferred variability of accretion onto the binary, making it more regular in time. We also investigate other potential sources of systematic error, such as the precise from of gravitational softening and previously employed simplifications to the viscous stress tensor. Because of the strong dependence of the orbital evolution of the binary on both the torque and distribution of mass between binary components, the sink methods employed can have a significant effect on the inferred orbital evolution of the binary.

We present high-resolution spectroscopy of two nearby white dwarfs with inconsistent spectroscopic and parallax distances. The first one, PG 1632+177, is a 13th magnitude white dwarf only 25.6 pc away. Previous spectroscopic observations failed to detect any radial velocity changes in this star. Here, we show that PG 1632+177 is a 2.05 d period double-lined spectroscopic binary (SB2) containing a low-mass He-core white dwarf with a more-massive, likely CO-core white dwarf companion. After L 870-2, PG 1632+177 becomes the second closest SB2 white dwarf currently known. Our second target, WD 1534+503, is also an SB2 system with an orbital period of 0.71 d. For each system, we constrain the atmospheric parameters of both components through a composite model-atmosphere analysis. We also present a new set of NLTE synthetic spectra appropriate for modeling high-resolution observations of cool white dwarfs, and show that NLTE effects in the core of the H$\alpha$ line increase with decreasing effective temperature. We discuss the orbital period and mass distribution of SB2 and eclipsing double white dwarfs with orbital constraints, and demonstrate that the observed population is consistent with the predicted period distribution from the binary population synthesis models. The latter predict more massive CO + CO white dwarf binaries at short ($<1$ d) periods, as well as binaries with several day orbital periods; such systems are still waiting to be discovered in large numbers.

In the next decade, cosmological surveys will have the statistical power to detect the absolute neutrino mass scale. N-body simulations of large-scale structure formation play a central role in interpreting data from such surveys. Yet these simulations are Newtonian in nature. We provide a quantitative study of the limitations to treating neutrinos, implemented as N-body particles, in N-body codes, focusing on the error introduced by neglecting special relativistic effects. Special relativistic effects are potentially important due to the large thermal velocities of neutrino particles in the simulation box. We derive a self-consistent theory of linear perturbations in Newtonian and non-relativistic neutrinos and use this to demonstrate that N-body simulations overestimate the neutrino free-streaming scale, and cause errors in the matter power spectrum that depend on the initial redshift of the simulations. For $z_{i} \lesssim 100$, and neutrino masses within the currently allowed range, this error is $\lesssim 0.5\%$, though represents an up to $\sim 10\%$ correction to the shape of the neutrino-induced suppression to the cold dark matter power spectrum. We argue that the simulations accurately model non-linear clustering of neutrinos so that the error is confined to linear scales.

To further shed light on whether pre-recombination models can resolve the Hubble tension, we explore constraints on the cosmic background evolution that are insensitive to early-universe physics. The analysis of the cosmic microwave background (CMB) anisotropy has been thought to highly rely on early-universe physics. However, we show that the difference between the sound horizon at recombination and that at the end of the drag epoch is small and insensitive to early-universe physics. This allows us to link the absolute sizes of the two horizons and treat them as free parameters. Jointly, the CMB peak angular size, Baryon Acoustic Oscillations (BAO), and Type Ia supernovae can be used as "early-universe-physics independent and uncalibrated cosmic standards", which measure the cosmic history from recombination to today. They can set strong and robust constraints on the post-recombination cosmic background, especially the matter density parameter with $\Omega_{\rm{m}}=0.302\pm0.008$ ($68\%$ C.L.). When combined with other non-local observations that are independent of or insensitive to early-universe physics, we obtain several constraints on $H_0$ that are currently more consistent with the Planck 2018 result than the local measurement results such as those based on Cepheids. These results suggest a tension between the post-recombination, but non-local, observations and the local measurements. This tension cannot be resolved by modifying the pre-recombination early universe physics.

The ease of interstellar rocket travel is an issue with implications for the long term fate of our own and other civilizations and for the much-debated number of technological civilizations in the Galaxy. We show that the physical barrier to interstellar travel can be greatly reduced if voyagers are patient, and wait for the close passage of another star. For a representative time of $\sim$1 Gyr, characteristic of the remaining time that Earth will remain habitable, one anticipates a passage of another star within $\sim 1500$~AU. This lowers the travel time for interstellar migration by $\sim$ two orders of magnitude compared with calculated travel times based on distances comparable to average interstellar separations (i.e., $\sim$1 pc) in the solar vicinity. We consider the implications for how long-lived civilizations may respond to stellar evolution, including the case of stars in wide binaries, and the difficulties of identifying systems currently undergoing a relevant close encounter. Assuming that life originates only around G-type stars, but migrates primarily to lower mass hosts when the original system becomes uninhabitable, the fraction of extant technological civilizations that exist as diaspora can be comparable to the fraction that still orbit their original host stars.

We investigate the impact of photochemical hazes and disequilibrium gases on the thermal structure of hot-Jupiters, using a detailed 1-D radiative-convective model. We find that the inclusion of photochemical hazes results in major heating of the upper and cooling of the lower atmosphere. Sulphur containing species, such as SH, S$_2$ and S$_3$ provide significant opacity in the middle atmosphere and lead to local heating near 1 mbar, while OH, CH, NH, and CN radicals produced by the photochemistry affect the thermal structure near 1 $\mu$bar. Furthermore we show that the modifications on the thermal structure from photochemical gases and hazes can have important ramifications for the interpretation of transit observations. Specifically, our study for the hazy HD 189733 b shows that the hotter upper atmosphere resulting from the inclusion of photochemical haze opacity imposes an expansion of the atmosphere, thus a steeper transit signature in the UV-Visible part of the spectrum. In addition, the temperature changes in the photosphere also affect the secondary eclipse spectrum. For HD 209458 b we find that a small haze opacity could be present in this atmosphere, at pressures below 1 mbar, which could be a result of both photochemical hazes and condensates. Our results motivate the inclusion of radiative feedback from photochemical hazes in general circulation models for a proper evaluation of atmospheric dynamics.

We search for unresolved X-ray emission from lensed sources in the FOV of 11 CLASH clusters with Chandra data. We consider the solid angle in the lens plane corresponding to a magnification $\mu>1.5$, that amounts to a total of ~100 arcmin$^2$. Our main goal is to assess the efficiency of massive clusters as cosmic telescopes to explore the faint end of X-ray extragalactic source population. We search for X-ray emission from strongly lensed sources identified in the optical, and perform an untargeted detection of lensed X-ray sources. We detect X-ray emission only in 9 out of 849 lensed/background optical sources. The stacked emission of the sources without detection does not reveal any signal in any band. Based on the untargeted detection, we find 66 additional X-ray sources that are consistent with being lensed sources. After accounting for completeness and sky coverage, we measure for the first time the soft- and hard-band number counts of lensed X-ray sources. The results are consistent with current modelization of the AGN population distribution. The distribution of de-lensed fluxes of the sources identified in moderately deep CLASH fields reaches a flux limit of ~$10^{-16}$ and ~$10^{-15}$ erg/s/cm$^{2}$ in the soft and hard bands, respectively. We conclude that, in order to match the depth of the CDFS exploiting massive clusters as cosmic telescopes, the required number of cluster fields is about two orders of magnitude larger than that offered by the 20 years Chandra archive. A significant step forward will be made when future X-ray facilities, with ~1' angular resolution and large effective area, will allow the serendipitous discovery of rare, strongly lensed high-$z$ X-ray sources, enabling the study of faint AGN activity in early Universe and the measurement of gravitational time delays in the X-ray variability of multiply imaged AGN.

We present the discovery of the fading radio transient FIRST J153350.8+272729. The source had a maximum observed 5-GHz radio luminosity of $8\times10^{39}$ erg s$^{-1}$ in 1986, but by 2019 had faded by a factor of nearly 400. It is located 0.15 arcsec from the center of a galaxy (SDSS J153350.89+272729) at 147 Mpc, which shows weak Type II Seyfert activity. We show that a tidal disruption event (TDE) is the preferred scenario for FIRST J153350.8+272729, although it could plausibly be interpreted as the afterglow of a long-duration gamma-ray burst. This is only the second TDE candidate to be first discovered at radio wavelengths. Its luminosity fills a gap between the radio afterglows of sub-relativistic TDEs in the local universe, and relativistic TDEs at high redshifts. The unusual properties of FIRST J153350.8+272729 (ongoing nuclear activity in the host galaxy, high radio luminosity) motivate more extensive TDE searches in untargeted radio surveys.

This paper aims to derive a map of relative planet occurrence rates that can provide constraints on the overall distribution of terrestrial planets around FGK stars. Based on the planet candidates in the Kepler DR25 data release, I first generate a continuous density map of planet distribution using a Gaussian kernel model and correct the geometric factor that the discovery space of a transit event decreases along with the increase of planetary orbital distance. Then I fit two exponential decay functions of detection efficiency along with the increase of planetary orbital distance and the decrease of planetary radius. Finally, the density map of planet distribution is compensated for the fitted exponential decay functions of detection efficiency to obtain a relative occurrence rate distribution of terrestrial planets. The result shows two regions with planet abundance: one corresponds to planets with radii between 0.5 and 1.5 R_Earth within 0.2 AU, the other corresponds to planets with radii between 1.5 and 3 R_Earth beyond 0.5 AU. It also confirms the features that may be caused by atmospheric evaporation: there is a vacancy of planets of sizes between 2.0 and 4.0 R_Earth inside of ~ 0.5 AU, and a valley with relatively low occurrence rates between 0.2 and 0.5 AU for planets with radii between 1.5 and 3.0 R_Earth.

Most active supermassive black holes (SMBH) in present-day galaxies are underfed and consist of low-luminosity active galactic nuclei (LLAGN). They have multiwavelength broadband spectral energy distributions (SED) dominated by non-thermal processes which are quite different from those of the brighter, more distant quasars. Modelling the observed SEDs of LLAGNs is currently challenging, given the large computational expenses required. In this work, we used machine learning (ML) methods to generate model SEDs and fit sparse observations of LLAGNs. Our ML model consisted of a neural network and reproduced with excellent precision the radio-to-X-rays emission from a radiatively inefficient accretion flow around a SMBH and a relativistic jet, at a small fraction of the computational cost. The ML method performs the fit $4 \times 10^5$ times faster than previous semianalytic models. As a proof-of-concept, we used the ML model to reproduce the SEDs of the LLAGNs M87, NGC 315 and NGC 4261.

We report detection of a line-like feature in the $\gamma$-ray spectrum of the blazar B0516$-$621, for which the data obtained with the Large Area Telescope onboard {\it Fermi Gamma-Ray Space Telescope (Fermi)} areanalyzed. The feature is at $\sim$7\,GeV and cannot be explained with radiation processes generally considered for jet emission of blazars. Instead, it could be a signal due to the oscillations between photons and axion-like particles (ALPs) in the source's jet. We investigate this possibility by fitting the spectrum with the photon-ALP oscillation model, and find thatthe parameter space of ALP mass $m_a\leq 10^{-8}$\,eV and the coupling constant (between photons and ALPs) $g_{a\gamma}$=1.17--1.49$\times10^{-10}$\,GeV$^{-1}$ can provide a fit to the line-like feature, while the magnetic field at the emission site of $\gamma$-rays is fixed at 0.7\,G. The ranges for $m_a$ and $g_{a\gamma}$ are in tension with those previously obtained from several experiments or methods, but on the other hand in line with some of the others. This spectral-feature case and its possible indication for ALP existence could be checked from similar studies of other blazar systems and also suggest a direction of effort for building future high-energy facilities that would have high sensitivities and spectral resolutions for searching for similar features.

We study the evolution of rapid neutron-capture process (r-process) isotopes in the Galaxy. We analyze relative contributions from core collapse supernovae (CCSNe), neutron star mergers (NSMs) and collapsars under a range of astrophysical conditions and nuclear input data. Although the r-process in each of these sites can lead to similar (or differing) isotopic abundances, our simulations reveal that the early contribution of r-process material to the Galaxy was dominated by CCSNe and collapsar r-process nucleosynthesis, while the NSM contribution is unavoidably delayed even under the assumption of the shortest possible minimum merger time.

We report an analysis of the dust disk around DM~Tau, newly observed with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.3 mm. The ALMA observations with high sensitivity (8.4~$\mu$Jy/beam) and high angular resolution (35~mas, 5.1~au) detect two asymmetries on the ring at $r\sim$20~au. They could be two vortices in early evolution, the destruction of a large scale vortex, or double continuum emission peaks with different dust sizes. We also found millimeter emissions with $\sim$50~$\mu$Jy (a lower limit dust mass of 0.3~$M_{\rm Moon}$) inside the 3-au ring. To characterize these emissions, we modeled the spectral energy distribution (SED) of DM~Tau using a Monte Carlo radiative transfer code. We found that an additional ring at $r=$ 1~au could explain both the DM~Tau SED and the central point source. The disk midplane temperature at the 1-au ring calculated in our modeling is less than the typical water sublimation temperature of 150~K, prompting the possibility of forming small icy planets there.

We make extensive numerical studies of masses and radii of proto-neutron stars during the first second after their birth in core-collapse supernova events. We use a quasi-static approach for the computation of proto-neutron star structure, built on parameterized entropy and electron fraction profiles, that are then evolved with neutrino cooling processes. We vary the equation of state of nuclear matter, the proto-neutron star mass and the parameters of the initial profiles, to take into account our ignorance of the supernova progenitor properties. We show that if masses and radii of a proto-neutron star can be determined in the first second after the birth, e.g. from gravitational wave emission, no information could be obtained on the corresponding cold neutron star and therefore on the cold nuclear equation of state. Similarly, it seems unlikely that any property of the proto-neutron star equation of state (hot and not beta-equilibrated) could be determined either, mostly due to the lack of information on the entropy, or equivalently temperature, distribution in such objects.

This work includes a comprehensive analysis of the Kepler detached eclipsing binary system KIC 6629588 that aims to the detailed study of the oscillation properties of its pulsating component. Ground-based spectroscopic observations were obtained and used to classify the components of the system. The spectroscopic results were used as constrains for the modelling of the short-cadence Kepler light curves and for the estimation of the absolute parameters of the components. Furthermore, the light curves residuals are analyzed using Fourier transformation techniques in order to search for pulsation frequencies. The primary component of the system is identified as a $\delta$ Scuti star that pulsates in seven eigenfrequencies in the range 13-22 d$^{-1}$, while more than 270 combination frequencies were also detected. The absolute and the oscillation parameters of this pulsating star are compared with those of other $\delta$ Scuti stars-members of detached binary systems using evolutionary and correlation diagrams. Finally, the distance of the system is also estimated.

We present the results of high signal-to-noise ratio VLT spectropolarimetry of a representative sample of 25 bright type 1 AGN at z<0.37, of which nine are radio-loud. The sample covers uniformly the 5100 A optical luminosity at $L_{5100}\sim 10^{44}-10^{46}$ erg s$^{-1}$, and H$\alpha$ width at FWHM$\sim 1000-10,000$~ km/s. We derive the continuum and the H$\alpha$ polarization amplitude, polarization angle, and angle swing across the line, together with the radio properties. We find the following: 1. The broad line region (BLR) and continuum polarization are both produced by a single scattering medium. 2. The scattering medium is equatorial, and at right angle to the system axis. 3. The scattering medium is located at or just outside the BLR. The continuum polarization and the H$\alpha$ polarization angle swing, can both serve as an inclination indicator. The observed line width is found to be affected by inclination, which can lead to an underestimate of the black hole mass by a factor of $\sim 5$ for a close-to face-on view. The line width measured in the polarized flux overcomes the inclination bias, and provides a close-to equatorial view of the BLR in all AGN, which allows to reduce the inclination bias in the BLR based black hole mass estimates.

We inspect the evolution of SIRs from Earth to Mars (distance range 1-1.5 AU) over the declining phase of solar cycle 24 (2014-2018). So far, studies only analyzed SIRs measured at Earth and Mars at different times. We compare existing catalogs for both heliospheric distances and arrive at a clean dataset for the identical time range. This allows a well-sampled statistical analysis and for the opposition phases of the planets an in-depth analysis of SIRs as they evolve with distance. We use in-situ solar wind data from OMNI and the MAVEN spacecraft as well as remote sensing data from SDO. A superposed epoch analysis is performed for bulk speed, proton density, temperature, magnetic field magnitude and total perpendicular pressure. Additionally, a study of events during the two opposition phases of Earth and Mars in the years 2016 and 2018 is conducted. SIR related coronal holes with their area as well as their latitudinal and longitudinal extent are extracted and correlated to the maximum bulk speed and duration of the corresponding high speed solar wind streams following the stream interaction regions. We find that while the entire solar wind HSS shows no expansion as it evolves from Earth to Mars, the crest of the HSS profile broadens by about 17%, and the magnetic field and total pressure by about 45% around the stream interface. The difference between the maximum and minimum values in the normalized superposed profiles increases slightly or stagnates from 1-1.5 AU for all parameters, except for the temperature. A sharp drop at zero epoch time is observed in the superposed profiles for the magnetic field strength at both heliospheric distances. Maximum solar wind speed has a stronger dependence on the latitudinal extent of the respective coronal hole than on its longitudinal extent. We arrive at an occurrence rate of fast forward shocks three times as high at Mars than at Earth.

Cosmological models predict that galaxies forming in the early Universe experience a chaotic phase of gas accretion and star formation, followed by gas ejection due to feedback processes. Galaxy bulges may assemble later via mergers or internal evolution. Here we present submillimeter observations (with spatial resolution of 700 parsecs) of ALESS 073.1, a starburst galaxy at redshift z~5, when the Universe was 1.2 billion years old. This galaxy's cold gas forms a regularly rotating disk with negligible noncircular motions. The galaxy rotation curve requires the presence of a central bulge in addition to a star-forming disk. We conclude that massive bulges and regularly rotating disks can form more rapidly in the early Universe than predicted by models of galaxy formation.

AIMS: A few well studied cataclysmic variables (CVs) have shown discrete characteristic frequencies of fast variability; the most prominent ones are around log(f/Hz) $\simeq$ -3. Because we still have only small number statistics, we obtained a new observation to test whether this is a general characteristic of CVs, especially if mass transfer occurs at a high rate typical for dwarf nova in outbursts, in the so called "high state". METHODS: We analyzed optical Kepler data of the quiescent nova and intermediate polar V4743 Sgr. This system hosts a white dwarf accreting through a disk in the high state. We calculated the power density spectra, and searched for break or characteristic frequencies. Our goal is to assess whether the mHz frequency of the flickering is a general characteristic. RESULTS: V4743 Sgr has a clear break frequency at log(f/Hz) $\simeq$ -3. This detection increases the probability that the mHz characteristic frequency is a general feature of CVs in the high state, from 69% to 91%. Furthermore, we propose the possibility that the variability is generated by similar mechanism as in the nova-like system MV Lyr, which would make V4743 Sgr unique.

We present the sharpest and deepest near infrared photometric analysis of the core of R136, a newly formed massive star cluster at the centre of the 30 Doradus star forming region in the Large Magellanic Cloud. We used the extreme adaptive optics of the SPHERE focal instrument implemented on the ESO Very Large Telescope and operated in its IRDIS imaging mode, for the second time with longer exposure time in the H- and K filters. Our aim was to (i) increase the number of resolved sources in the core of R136, and (ii) to compare with the first epoch to classify the properties of the detected common sources between the two epochs. Within the field of view (FOV) of 10.8"x12.1" (2.7pc x3.0pc), we detected 1499 sources in both H and K filters, for which 76% of these sources have visual companions closer than 0.2". The larger number of detected sources, enabled us to better sample the mass function (MF). The MF slopes are estimated at ages of 1, 1.5 and 2 Myr, at different radii, and for different mass ranges. The MF slopes for the mass range of 10-300 solar-mass are about 0.3 dex steeper than the mass range of 3-300 solar-mass, for the whole FOV and different radii. Comparing the JHK colours of 790 sources common in between the two epochs, 67% of detected sources in the outer region (r >3") are not consistent with evolutionary models at 1-2 Myr and with extinctions similar to the average cluster value, suggesting an origin from ongoing star formation within 30 Doradus, unrelated to R136.

It is not yet entirely clear whether Mars began as a warm and wet planet that evolved towards the present-day cold and dry body or if it always was cold and dry with just some sporadic episodes of liquid water on its surface. An important clue into this question can be gained by studying the earliest evolution of the Martian atmosphere and whether it was dense and stable to maintain a warm and wet climate or tenuous and susceptible to strong atmospheric escape. We discuss relevant aspects for the evolution and stability of a potential early Martian atmosphere. This contains the solar EUV flux evolution, the formation timescale and volatile inventory of the planet including volcanic degassing, impact delivery and removal, the loss of a catastrophically outgassed steam atmosphere, atmosphere-surface interactions, and thermal and non-thermal escape processes affecting any secondary atmosphere. While early non-thermal escape at Mars before 4 billion years ago (Ga) is poorly understood, particularly in view of its ancient intrinsic magnetic field, research on thermal escape processes indicate that volatile delivery and volcanic degassing cannot counterbalance the strong thermal escape. Therefore, a catastrophically outgassed steam atmosphere of several bars of CO2 and H2O, or CO and H2 for reduced conditions, could have been lost within just a few million years (Myr). Thereafter, Mars likely could not build up a dense secondary atmosphere during its first ~400 Myr but might only have possessed an atmosphere sporadically during events of strong volcanic degassing, potentially also including SO2. This indicates that before ~4.1 Ga Mars indeed might have been cold and dry. A denser CO2- or CO-dominated atmosphere, however, might have built up afterwards but must have been lost later-on due to non-thermal escape processes and sequestration into the ground.

The KM3NeT research infrastructure is under construction in the Mediterranean Sea. It consists of two water Cherenkov neutrino detectors, ARCA and ORCA, aimed at neutrino astrophysics and oscillation research, respectively. Instrumenting a large volume of sea water with $\sim$ 6,200 optical modules comprising a total of $\sim$ 200,000 photomultiplier tubes, KM3NeT will achieve sensitivity to $\sim$ 10 MeV neutrinos from Galactic and near-Galactic core-collapse supernovae through the observation of coincident hits in photomultipliers above the background. In this paper, the sensitivity of KM3NeT to a supernova explosion is estimated from detailed analyses of background data from the first KM3NeT detection units and simulations of the neutrino signal. The KM3NeT observational horizon (for a $5\,\sigma$ discovery) covers essentially the Milky-Way and for the most optimistic model, extends to the Small Magellanic Cloud ($\sim$ 60 kpc). Detailed studies of the time profile of the neutrino signal allow assessment of the KM3NeT capability to determine the arrival time of the neutrino burst with a few milliseconds precision for sources up to 5$-$8 kpc away, and detecting the peculiar signature of the \textit{standing accretion shock instability} if the core-collapse supernova explosion happens closer than 3$-$5 pc, depending on the progenitor mass. KM3NeT's capability to measure the neutrino flux spectral parameters is also presented.

We report the detection of a large sample of high-$\alpha$-metal-rich stars on the low giant branch with $2.6<logg<3.3$ dex in the LAMOST-MRS survey. This special group corresponds to an intermediate-age population of $5-9$ Gyr based on the $[Fe/H]$-$[C/N]$ diagram and age-$[C/N]$ calibration. A comparison group is selected to have solar $\alpha$ ratio at super metallicity, which is young and has a narrow age range around 3 Gyr. Both groups have thin-disk like kinematics but the former shows slightly large velocity dispersions. The special group shows a larger extension in vertical distance toward 1.2 kpc, a second peak at smaller Galactic radius and a larger fraction of super metal rich stars with $[Fe/H]>0.2$ than the comparison group. These properties strongly indicate its connection with the outer bar/bulge region at $R=3-5$ kpc. A tentative interpretation of this special group is that its stars were formed in the X-shaped bar/bulge region, close to its corotation radius, where radial migration is the most intense, and brings them to present locations at 9 kpc and beyond. Low eccentricities and slightly outward radial excursions of its stars are consistent with this scenario. Its kinematics (cold) and chemistry ($[\alpha/Fe]$ $\sim 0.1$) further support the formation of the instability-driven X-shaped bar/bulge from the thin disk.

A signification fraction of Galactic massive stars (> 8Mo) are ejected from their parent cluster and supersonically sail away through the interstellar medium (ISM). The winds of these fast-moving stars blow asymmetric bubbles thus creating a circumstellar environment in which stars eventually die with a supernova explosion. The morphology of the resulting remnant is largely governed by the circumstellar medium of the defunct progenitor star. In this paper, we present 2D magneto-hydrodynamical simulations investigating the effect of the ISM magnetic field on the shape of the supernova remnants of a 35Mo star evolving through a Wolf-Rayet phase and running with velocity 20 and 40 km/s, respectively. A 7 microG ambient magnetic field is sufficient to modify the properties of the expanding supernova shock front and in particular to prevent the formation of filamentary structures. Prior to the supernova explosion, the compressed magnetic field in the circumstellar medium stabilises the wind/ISM contact discontinuity in the tail of the wind bubble. A consequence is a reduced mixing efficiency of ejecta and wind materials in the inner region of the remnant, where the supernova shock wave propagates. Radiative transfer calculations for synchrotron emission reveal that the non-thermal radio emission has characteristic features reflecting the asymmetry of exiled core-collapse supernova remnants from Wolf-Rayet progenitors. Our models are qualitatively consistent with the radio appearance of several remnants of high-mass progenitors, namely the bilateral G296.5+10.0 and the shell-type remnants CTB109 and Kes 17, respectively.

The prediction of solar activity is important for advanced technologies and space activities. The peak sunspot number (SSN), which can represent the solar activity, has declined continuously in the past four solar cycles (21$-$24), and the Sun would experience a Dalton-like minimum, or even the Maunder-like minimum, if the declining trend continues in the following several cycles, so that the predictions of solar activity for cycles 25 and 26 are crucial. In Qin & Wu, 2018, ApJ, we established an SSN prediction model denoted as two-parameter modified logistic prediction (TMLP) model, which can predict the variation of SSNs in a solar cycle if the start time of the cycle has been determined. In this work, we obtain a new model denoted as TMLP-extension (TMLP-E). If the start time of a cycle $n$ is already known, TMLP-E can predict the variation of SSNs in the cycle $n+1$. Cycle 25 is believed to start in December 2019, so that the predictions of cycles 25 and 26 can be made with our models. It is found that the predicted solar maximum, ascent time, and cycle length are 115.1, 4.84 yr, and 11.06 yr, respectively, for cycle 25, and 107.3, 4.80 yr, and 10.97 yr, respectively, for cycle 26. The solar activities of cycles 25 and 26 are predicted to be at the same level as that of cycle 24, but will not decrease further. We therefore suggest that the cycles 24$-$26 are at a minimum of Gleissberg cycle.

Aims. We investigate the configuration of a complex flux rope above a {\delta} sunspot region in NOAA AR 11515, and its eruptive expansion during a confined M5.3-class flare. Methods. We study the formation of the {\delta} sunspot using continuum intensity images and photospheric vector magnetograms provided by SDO/HMI. We use EUV and UV images provided by SDO/AIA, and hard X-ray emission recorded by RHESSI to investigate the eruptive details. The coronal magnetic field is extrapolated with a non-linear force free field (NLFFF) method, based on which the flux rope is identified by calculating the twist number Tw and squashing factor Q. We search the null point via a modified Powell hybrid method. Results. The collision between two emerging spot groups form the {\delta} sunspot. A bald patch (BP) forms at the collision location, above which a complex flux rope is identified. The flux rope has multiple layers, with one compact end and one bifurcated end, having Tw decreasing from the core to the boundary. A null point is located above the flux rope. The eruptive process consists of precursor flaring at a 'v'-shaped coronal structure, rise of the filament, and flaring below the filament, corresponding well with the NLFFF topological structures, including the null point and the flux rope with BP and hyperbolic flux tube (HFT). Two sets of post-flare loops and three flare ribbons support the bifurcation configuration of the flux rope. Conclusions. The precursor reconnection, which occurs at the null point, weakens the overlying confinement to allow the flux rope to rise, fitting the breakout model. The main phase reconnection, which may occur at the BP or HFT, facilitates the flux rope rising. The results suggest that the {\delta} spot configuration presents an environment prone to the formation of complex magnetic configurations which will work together to produce activities.

The mismatch in the value of the Hubble constant from low- and high-redshift observations may be recast as a discrepancy between the low- and high-redshift determinations of the luminosity of Type Ia supernovae, the latter featuring an absolute magnitude which is $\approx 0.2$~mag lower. Here, we propose that a rapid transition in the value of the relative effective gravitational constant $\mu_G\equiv\frac{G_{\rm eff}}{G_N}$ at $z_t\simeq 0.01$ could explain the lower luminosity (higher magnitude) of local supernovae, thus solving the $H_0$ crisis. A model that features $\mu_G = 1$ for $z \lesssim 0.01$ but $\mu_G \simeq 0.9$ for $z \gtrsim 0.01$ is trivially consistent with local gravitational constraints but would raise the Chandrasekhar mass and so decrease the absolute magnitude of Type Ia supernovae at $z \gtrsim 0.01$ by the required value of $\approx 0.2$~mag. Such a rapid transition of the effective gravitational constant would not only resolve the Hubble tension but it would also help resolve the growth tension as it would reduce the growth of density perturbations without affecting the Planck/$\Lambda$CDM background expansion.

The recently-discovered giant exoplanet HR5183b exists on a wide, highly-eccentric orbit ($a=18$\,au, $e=0.84$). Its host star possesses a common proper-motion companion which is likely on a bound orbit. In this paper, we explore scenarios for the excitation of the eccentricity of the planet in binary systems such as this, considering planet-planet scattering, Lidov-Kozai cycles from the binary acting on a single-planet system, or Lidov-Kozai cycles acting on a two-planet system that also undergoes scattering. Planet-planet scattering, in the absence of a binary companion, has a $2.8-7.2\%$ probability of pumping eccentricities to the observed values in our simulations, depending on the relative masses of the two planets. Lidov-Kozai cycles from the binary acting on an initially circular orbit can excite eccentricities to the observed value, but require very specific orbital configurations for the binary and overall there is a low probability of catching the orbit at the high observed high eccentricity ($0.6\%$). The best case is provided by planet-planet scattering in the presence of a binary companion: here, the scattering provides the surviving planet with an initial eccentricity boost that is subsequently further increased by Kozai cycles from the binary. We find a success rate of $10.4\%$ for currently observing $e\ge0.84$ in this set-up. The single-planet plus binary and two-planet plus binary cases are potentially distinguishable if the mutual inclination of the binary and the planet can be measured, as the latter permits a broader range of mutual inclinations. The combination of scattering and Lidov-Kozai forcing may also be at work in other wide-orbit eccentric giant planets, which have a high rate of stellar binary companions.

Young exoplanets can offer insight into the evolution of planetary atmospheres, compositions, and architectures. We present the discovery of the young planetary system TOI 451 (TIC 257605131, Gaia DR2 4844691297067063424). TOI 451 is a member of the 120-Myr-old Pisces--Eridanus stream (Psc--Eri). We confirm membership in the stream with its kinematics, its lithium abundance, and the rotation and UV excesses of both TOI 451 and its wide binary companion, TOI 451 B (itself likely an M dwarf binary). We identified three candidate planets transiting in the TESS data and followed up the signals with photometry from Spitzer and ground-based telescopes. The system comprises three validated planets at periods of 1.9, 9.2 and 16 days, with radii of 1.9, 3.1, and 4.1 Earth radii, respectively. The host star is near-solar mass with V=11.0 and H=9.3 and displays an infrared excess indicative of a debris disk. The planets offer excellent prospects for transmission spectroscopy with HST and JWST, providing the opportunity to study planetary atmospheres that may still be in the process of evolving.

The topic of this review paper is on the influence of solar wind turbulence on shock propagation and its consequence on the acceleration and transport of energetic particles at shocks. As the interplanetary shocks sweep through the turbulent solar wind, the shock surfaces fluctuate and ripple in a range of different scales. We discuss particle acceleration at rippled shocks in the presence of ambient solar-wind turbulence. This strongly affects particle acceleration and transport of energetic particles (both ions and electrons) at shock fronts. In particular, we point out that the effects of upstream turbulence is critical for understanding the variability of energetic particles at shocks. Moreover, the presence of pre-existing upstream turbulence significantly enhances the trapping near the shock of low-energy charged particles, including near the thermal energy of the incident plasma, even when the shock propagates normal to the average magnetic field. Pre-existing turbulence, always present in space plasmas, provides a means for the efficient acceleration of low-energy particles and overcoming the well known injection problem at shocks.

In the present work we study Hinode/EIS observations of an active region taken before, during and after a small C2.0 flare in order to monitor the evolution of the magnetic field evolution and its relation to the flare event. We find that while the flare left the active region itself unaltered, the event included a large Magnetic Field Enhancement (MFE), which consisted of a large increase of the magnetic field to strengths just short of 500~G in a rather small region where no magnetic field was measured before the flare. This MFE is observed during the impulsive phase of the flare at the footpoints of flare loops, its magnetic energy is sufficient to power the radiative losses of the entire flare, and has completely dissipated after the flare. We argue that the MFE might occur at the location of the reconnection event triggering the flare, and note that it formed within 22 minutes of the flare start (as given by the EIS raster return time). These results open the door to a new line of studies aimed at determining whether MFEs 1) can be flare precursor events, 2) can be used for Space Weather forecasts; and 3) what advance warning time they could allow; as well as to explore which physical processes lead to their formation and dissipation, whether such processes are the same in both long-duration and impulsive flares, and whether they can be predicted by theoretical models.

In the frequency-domain multiplexing (FDM) scheme, transition-edge sensors (TES) are individually coupled to superconducting LC filters and AC biased at MHz frequencies through a common readout line. To make efficient use of the available readout bandwidth and to minimize the effect of non-linearities, the LC resonators are usually designed to be on a regular grid. The lithographic processes however pose a limit on the accuracy of the effective filter resonance frequencies. Off-resonance bias carriers could be used to suppress the impact of intermodulation distortions, which nonetheless would significantly affect the effective bias circuit and the detector spectral performance. In this paper we present a frequency shift algorithm (FSA) to allow off-resonance readout of TES's while preserving the on-resonance bias circuit and spectral performance, demonstrating its application to the FDM readout of a X-ray TES microcalorimeter array. We discuss the benefits in terms of mitigation of the impact of intermodulation distortions at the cost of increased bias voltage and the scalability of the algorithm to multi-pixel FDM readout. We show that with FSA, in multi-pixel and frequencies shifted on-grid, the line noises due to intermodulation distortion are placed away from the sensitive region in the TES response and the X-ray performance is consistent with the single-pixel, on-resonance level.

The shell type supernova remnant (SNR) Cas A exhibits structures at nearly all angular scales. Previous studies show the angular power spectrum $(C_{\ell})$ of the radio emission to be a broken power law, consistent with MHD turbulence. The break has been identified with the transition from 2D to 3D turbulence at the angular scale corresponding to the shell thickness. Alternatively, this can also be explained as 2D inverse cascade driven by energy injection from knot-shock interactions. Here we present $C_{\ell}$ measured from archival VLA $5$GHz (C band) data, and Chandra X-ray data in the energy ranges ${\rm A}=0.6-1.0 \, \, {\rm keV}$ and ${\rm B} =4.2-6.0 \, \, {\rm keV}$, both of which are continuum dominated. The different emissions all trace fluctuations in the underlying plasma and possibly also the magnetic field, and we expect them to be correlated. We quantify this using the cross $C_{\ell}$ between the different emissions. We find that X-ray B is strongly correlated with both radio and X-ray A, however X-ray A is only very weakly correlated with radio. This supports a picture where X-ray A is predominantly thermal bremsstrahlung whereas X-ray B is a composite of thermal bremsstrahlung and non-thermal synchrotron emission. The various $C_{\ell}$ measured here, all show a broken power law behaviour. However, the slopes are typically shallower than those in radio and the position of the break also corresponds to smaller angular scales. These findings provide observational inputs regarding the nature of turbulence and the emission mechanisms in Cas A.

Post-reionisation 21cm intensity mapping experiments target the spectral line of neutral hydrogen (HI) resident in dark matter haloes. According to the halo model, these discrete haloes trace the continuous dark matter density field down to a certain scale, which is dependent on the halo physical size. The halo physical size defines an exclusion region which leaves imprints on the statistical properties of HI. We show how the effect of exclusion due to the finite halo size impacts the HI power spectrum, with the physical boundary of the host halo given by the splashback radius. Most importantly, we show that the white noise-like feature that appears in the zero-momentum limit of the power spectrum is exactly cancelled when the finite halo size is taken into consideration. This cancellation in fact applies to all tracers of dark matter density field, including galaxies. Furthermore, we show that the exclusion due to finite halo size leads to a sub-Poissonian noise signature on large scales, consistent with the results from N-body simulations

The discovery of periodicity in the arrival times of the fast radio bursts (FRB) from two repeating sources poses a potential challenge to oft studied magnetar scenarios. However, models which postulate that FRB emission results from relativistic magnetized shocks, or magnetic reconnection in a striped outflow, are not necessarily specific to magnetar engines, instead requiring only the impulsive injection of relativistic energy into a dense magnetized medium. Motivated thus, we outline a new scenario in which FRBs are powered by short-lived relativistic outflows ("flares") from accreting black hole or neutron star systems, which propagate into the cavity of the pre-existing ("quiescent") jet. In order to reproduce FRB luminosities and rates, we are driven to consider binaries of stellar-mass compact objects undergoing super-Eddington mass transfer, similar to those which characterize some ultra-luminous X-ray (ULX) sources. Indeed, the host galaxies of FRBs, and their spatial offsets within their hosts, show broad similarities to those of ULX. Periodicity on timescales of days to years could be attributed to precession (e.g., Lens-Thirring) of the polar accretion funnel, along which the FRB emission is geometrically and relativistically beamed, across the observer line of sight. Accounting for the most luminous FRBs via accretion power may require a population of binaries undergoing brief-lived phases of unstable (dynamical timescale) mass transfer. This could lead to secular evolution in the burst properties of some repeating FRB sources on timescales as short as months to years, followed by a transient optical/IR counterpart akin to a luminous red nova or dusty common envelope transient. We encourage targeted FRB searches of known ULX sources.

Interacting dark matter (DM) - dark energy (DE) models have been intensively investigated in the literature for their ability to fit various data sets as well as to explain some observational tensions persisting within the $\Lambda$CDM cosmology. In this work, we employ Gaussian processes (GP) algorithm to perform a joint analysis by using the geometrical cosmological probes such as Cosmic chronometers, Supernova Type Ia, Baryon Acoustic Oscillations and the H0LiCOW lenses sample to infer a reconstruction of the coupling function between the dark components in a general framework, where the DE can assume a dynamical character via its equation of state. In addition to the joint analysis with these data, we simulate a catalog with standard siren events from binary neutron star mergers, within the sensitivity predicted by the Einstein Telescope, to reconstruct the dark sector coupling with more accuracy in a robust way. We find that the particular case, where $w = -1$ is fixed on the DE nature, has a statistical preference for an interaction in the dark sector at late times. In the general case, where $w(z)$ is analyzed, we find no evidence for such dark coupling, and the predictions are compatible with the $\Lambda$CDM paradigm. When the mock events of the standard sirens are considered to improve the kernel in GP predictions, we find preference for an interaction in the dark sector at late times.

This work evaluates the viability of polyimide and parylene-C for passivation of lithium-drifted silicon (Si(Li)) detectors. The passivated Si(Li) detectors will form the particle tracker and X-ray detector of the General Antiparticle Spectrometer (GAPS) experiment, a balloon-borne experiment optimized to detect cosmic antideuterons produced in dark matter annihilations or decays. Successful passivation coatings were achieved by thermally curing polyimides, and the optimized coatings form an excellent barrier against humidity and organic contamination. The passivated Si(Li) detectors deliver $\lesssim\,4$ keV energy resolution (FWHM) for 20$-$100 keV X-rays while operating at temperatures of $-$35 to $-45\,^{\circ}$C. This is the first reported successful passivation of Si(Li)-based X-ray detectors operated above cryogenic temperatures.

We discuss the current state of knowledge of terrestrial planet formation from the aspects of different planet formation models and isotopic data from 182Hf-182W, U-Pb, lithophile-siderophile elements, 48Ca/44Ca isotope samples from planetary building blocks, 36Ar/38Ar, 20Ne/22Ne, 36Ar/22Ne isotope ratios in Venus' and Earth's atmospheres, the expected solar 3He abundance in Earth's deep mantle and Earth's D/H sea water ratios that shed light on the accretion time of the early protoplanets. Accretion scenarios that can explain the different isotope ratios, including a Moon-forming event after ca. 50 Myr, support the theory that the bulk of Earth's mass (>80%) most likely accreted within 10-30 Myr. From a combined analysis of the before mentioned isotopes, one finds that proto-Earth accreted 0.5-0.6 MEarth within the first ~4-5 Myr, the approximate lifetime of the protoplanetary disk. For Venus, the available atmospheric noble gas data are too uncertain for constraining the planet's accretion scenario accurately. However, from the available Ar and Ne isotope measurements, one finds that proto-Venus could have grown to 0.85-1.0 MVenus before the disk dissipated. Classical terrestrial planet formation models have struggled to grow large planetary embryos quickly from the tiniest materials within the typical lifetime of protoplanetary disks. Pebble accretion could solve this long-standing time scale controversy. Pebble accretion and streaming instabilities produce large planetesimals that grow into Mars-sized and larger planetary embryos during this early accretion phase. The later stage of accretion can be explained well with the Grand-Tack, annulus or depleted disk models. The relative roles of pebble accretion and planetesimal accretion/giant impacts are poorly understood and should be investigated with N-body simulations that include pebbles and multiple protoplanets.

We extend the work of Ryan on mapping the spacetime of the central object of an extreme mass-ratio inspiral (EMRI) by using gravitational waves (GWs) emitted by the system, which may be observed in future missions such as LISA. Whether the central object is a black hole or not can be probed by observing the phasing of these waves, which carry information about its mass and spin multipole moments. We go beyond the phase terms found by Ryan, which were obtained in the quadrupolar approximation of the point-particle limit, and derive terms up to the fifth post-Newtonian (PN) order. Since corrections due to horizon absorption (i.e., if the central object is a black hole) and tidal heating appear by that order, at 2.5PN and 5PN, respectively, we include them here. Corrections due to the motion of the central object, which was addressed only partially by Ryan, are included as well. Additionally, we obtain the contribution of the higher order radiative multipole moments. For the tidal interaction, our results have been derived in the approximation of the Newtonian tidal field. Therefore, in the potential for tidal field only the contribution due to the mass of the central object has been included as well. Using these results we argue that it might be possible for LISA to probe if the central object in an EMRI has a horizon or not. We also discuss how our results can be used to test the No-hair theorem from the inspiral phase of such systems.

We study the quasi-normal modes (QNMs) of static, spherically symmetric black holes in $f(R)$ theories. We show how these modes in theories with non-trivial $f(R)$ are fundamentally different from those in General Relativity. In the special case of $f(R) = \alpha R^2$ theories, it has been recently argued that iso-spectrality between scalar and vector modes breaks down. Here, we show that such a break down is quite general across all $f(R)$ theories, as long as they satisfy $f''(0)/(1+f''(0)) \neq 0$, where a prime denotes derivative of the function with respect to its argument. We specifically discuss the origin of the breaking of isospectrality. We also show that along with this breaking the QNMs receive a correction that arises when $f''(0)/(1+f'(0)) \neq 0$ owing to the inhomogeneous term that it introduces in the mode equation. We discuss how these differences affect the "ringdown" phase of binary black hole mergers and the possibility of constraining $f(R)$ models with gravitational-wave observations. We also find that even though the iso-spectrality is broken in $f(R)$ theories, in general, nevertheless in the corresponding scalar-tensor theories in the Einstein frame it is unbroken.

Over 30 years ago, Barrow & Tipler proposed the principle according to which the action integrated over the entire 4-manifold describing the universe should be finite. Here we explore the cosmological consequences of a related criterion, namely that semi-classical transition amplitudes from the early universe up to current field values should be well defined. On a classical level, our criterion is weaker than the Barrow-Tipler principle, but it has the advantage of being sensitive to quantum effects. We find significant consequences for early universe models, in particular: eternal inflation and strictly cyclic universes are ruled out. Within general relativity, the first phase of evolution cannot be inflationary, and it can be ekpyrotic only if the scalar field potential is trustworthy over an infinite field range. Quadratic gravity eliminates all non-accelerating backgrounds near a putative big bang (thus imposing favourable initial conditions for inflation), while the expected infinite series of higher-curvature quantum corrections eliminates Lorentzian big bang spacetimes altogether. The scenarios that work best with the principle of finite amplitudes are the no-boundary proposal, which gives finite amplitudes in all dynamical theories that we have studied, and string-inspired loitering phases. We also comment on the relationship of our proposal to the swampland conjectures.

Dark matter may self-interact through a continuum of low-mass states. This happens if dark matter couples to a strongly-coupled nearly-conformal hidden sector. This type of theory is holographically described by brane-localized dark matter interacting with bulk fields in a slice of 5D anti-de Sitter space. The long-range potential in this scenario depends on a non-integer power of the spatial separation, in contrast to the Yukawa potential generated by the exchange of a single 4D mediator. The resulting self-interaction cross section scales like a non-integer power of velocity. We identify the Born, classical and resonant regimes and investigate them using state-of-the-art numerical methods. We demonstrate the viability of our continuum-mediated framework to address the astrophysical small-scale structure anomalies. Investigating the continuum-mediated Sommerfeld enhancement, we demonstrate that a pattern of resonances can occur depending on the non-integer power. We conclude that continuum mediators introduce novel power-law scalings which open new possibilities for dark matter self-interaction phenomenology.

Probing the QCD axion dark matter (DM) hypothesis is extremely challenging as the axion interacts very weakly with Standard Model particles. We propose a new avenue to test the QCD axion DM via transient radio signatures coming from encounters between neutron stars (NSs) and axion minihalos around primordial black holes (PBHs). We consider a general QCD axion scenario in which the PQ symmetry breaking occurs before (or during) inflation coexisting with a small fraction of DM in the form of PBHs. The PBHs will unavoidably acquire around them axion minihalos with the typical length scale of parsecs. The axion density in the minihalos may be much higher than the local DM density, and the presence of these compact objects in the Milky Way today provides a novel chance for testing the axion DM hypothesis. We study the evolution of the minihalo mass distribution in the Galaxy accounting for tidal forces and estimate the encounter rate between NSs and the dressed PBHs. We find that the encounters give rise to transient line-like emission of radio frequency photons produced by the resonant axion-photon conversion in the NS magnetosphere and the characteristic signal could be detectable with the sensitivity of current and prospective radio telescopes.

The excited-state structure of atomic nuclei can modify nuclear processes in stellar environments. In this work, we study the influence of nuclear excitations on Urca cooling (repeated back-and-forth beta decay and electron capture in a pair of nuclear isotopes) in the crust and ocean of neutron stars. We provide for the first time an expression for Urca process neutrino luminosity which accounts for excited states of both members of an Urca pair. We use our new formula with state-of-the-art nuclear structure inputs to compute neutrino luminosities of candidate Urca cooling pairs. Our nuclear inputs consist of the latest experimental data supplemented with calculations using the projected shell model. We show that, in contrast with previous results that only consider the ground states of both nuclei in the pair, our calculated neutrino luminosities for different Urca pairs vary sensitively with the environment temperature and can be radically different from those obtained in the one transition approximation.