Papers for Friday, Aug 12 2022

Optimal error estimation is key to achieve accurate photometry and astrometry. Stellar fluxes and positions in high angular resolution images are typically measured with PSF fitting routines, such as StarFinder. However, the formal uncertainties computed by these software packages tend to seriously underestimate the relevant uncertainties. We present a new approach to deal with this problem using a resampling method to obtain robust and reliable uncertainties without loss of sensitivity.

We constrain the orbital period distribution of low-mass detached main-sequence eclipsing binaries (EBs) with light curves from the Zwicky Transient Facility (ZTF), which provides a well-understood selection function and sensitivity to faint stars. At short periods ($P_{\rm orb}\lesssim 2$ days), binaries are predicted to evolve significantly due to magnetic braking (MB), which shrinks orbits and ultimately brings detached binaries into contact. The period distribution is thus a sensitive probe of MB. We find that the intrinsic period distribution of low-mass ($0.1\lesssim M_1/M_{\odot} < 0.9$) binaries is basically flat (${\rm d}N/{\rm d}P_{\rm orb} \propto P_{\rm orb}^0$), from $P_{\rm orb}=10$ days down to the contact limit. This is strongly inconsistent with predictions of classical MB models based on the Skumanich relation, which are widely used in binary evolution calculations and predict ${\rm d}N/{\rm d}P_{\rm orb} \propto P_{\rm orb}^{7/3}$ at short periods. The observed distributions are best reproduced by models in which the magnetic field saturates at short periods, with a MB torque that scales roughly as $\dot{J}\propto P_{\rm orb}^{-1}$, as opposed to $\dot{J} \propto P_{\rm orb}^{-3}$ in the standard Skumanich law. We also find no significant difference between the period distributions of binaries containing fully and partially convective stars. Our results confirm that a saturated MB law, which was previously found to describe the spin-down of rapidly rotating isolated M dwarfs, also operates in tidally locked binaries. We advocate using saturated MB models in binary evolution calculations. Our work supports previous suggestions that MB in cataclysmic variables (CVs) is much weaker than assumed in the standard evolutionary model, unless mass transfer in CVs significantly strengthens MB.

The correlation between the broad line region radius and continuum luminosity ($R-L$ relation) of active galactic nuclei (AGN) is critical for single-epoch mass estimates of supermassive black holes (SMBHs). At $z \sim 1-2$, where AGN activity peaks, the $R-L$ relation is constrained by the reverberation mapping (RM) lags of the Mg II line. We present 25 Mg II lags from the Australian Dark Energy Survey (OzDES) RM project based on six years of monitoring. We define quantitative criteria to select good lag measurements and verify their reliability with simulations based on both the damped random walk stochastic model and the re-scaled, re-sampled versions of the observed lightcurves of local, well-measured AGN. Our sample significantly increases the number of Mg II lags and extends the $R-L$ relation to higher redshifts and luminosities. The relative iron line strength $\mathcal{R}_{\rm Fe}$ has little impact on the $R-L$ relation. The best-fit Mg II $R-L$ relation has a slope $\alpha = 0.39 \pm 0.08$ with an intrinsic scatter $\sigma_{\rm rl} = 0.15^{+0.03}_{-0.02}$. The slope is consistent with previous measurements and shallower than the H$\beta$ $R-L$ relation. The intrinsic scatter of the new $R-L$ relation is substantially smaller than previous studies and comparable to the intrinsic scatter of the H$\beta$ $R-L$ relation. Our new $R-L$ relation will enable more precise single-epoch mass estimates and SMBH demographic studies at cosmic noon.

Planes of satellites are observed around many galaxies. However, these observations are still considered a point of tension for the $\Lambda$CDM paradigm. We use the fully hydrodynamical cosmological $\Lambda$CDM state-of-the-art simulation Magneticum Pathfinder to investigate the existence of such planes over a large range of haloes, from Milky Way to galaxy cluster masses. To this end, we develop the Momentum in Thinnest Plane (MTP) method to identify planes and quantify the properties of their constituent satellites, considering both position and momentum. We find that thin planes ($20\%\cdot R_\mathrm{halo}$) containing at least $50\%$ of the total number of satellites can be found in almost all systems. In Milky Way mass-like systems, around 86\% of such planes are even aligned in momentum ($90\%$ of the total satellite momentum), where the fraction is smaller if more satellites are required to be inside the plane. We further find a mass dependency, with more massive systems exhibiting systematically thicker planes. This may point towards the change from continuous accretion of small objects along filaments and sheets for less massive haloes to the accretion of large objects (e.g., major mergers) dominating the growth of more massive haloes. There is no correlation between the existence of a plane and the main galaxy's morphology. Finally, we find a clear preference for the minor axes of the satellite planes and the host galaxy to be aligned, in agreement with recent observations.

We use the unprecedented resolution and depth of the JWST NIRCam Early Release Observations (ERO) at 1-5$\mu m$ to study the stellar mass, age, and metallicity of compact star clusters in the neighborhood of the host galaxies in the SMACS J0723.3-7327 galaxy cluster field at z = 0.39. The measured colors of these star clusters show a similar distribution as quiescent galaxies at the same redshift, but are >3 magnitudes fainter than the current depths of wide-field galaxy survey. The star clusters are unresolved in the NIRCam/F150W data suggesting sizes smaller than 50pc. This is significantly smaller than star forming clumps or dwarf galaxies in local galaxies. From fitting their photometry with simple stellar population (SSP) models, we find stellar metallicities consistent with 0.2-0.3Z$_{\odot}$ and ages of $5.0^{+0.5}_{-1.1}$ Gyrs. We rule out ages of <3 Gyrs at a 2$\sigma$ confidence, and metallicities <0.2 Z$_{\odot}$ and solar/super-solar at 4$\sigma$ significance. Assuming the mass-to-light ratio (1.22 M$_{\odot}$/L$_{\odot}$) obtained from the best-fit SSP, we estimate stellar masses of $4.1^{+3.5}_{-1.8}\times10^6\,{\rm M_{\odot}}$. These are between average masses of local globular clusters and dwarf galaxies. Our analysis suggests middle-aged globulars with relatively recent formation times at z=1-2, which could be subsequently stripped away from their host galaxies due to interactions in the cluster environment. However, we cannot rule out these objects being compact cores of stripped dwarf galaxies.

We investigate the impact of anomalous microwave emission (AME) on the radio-millimeter spectral energy distribution for three typical interstellar medium (ISM) conditions surrounding star-forming regions -- cold neutral medium, warm neutral medium, and photodissociation region -- by comparing the emissivities of three major contributors: free-free, thermal dust emission, and AME. In particular, for spinning nanoparticles (i.e., potential carriers of AME), we consider a known grain destruction mechanism due to a centrifugal force from spin-up processes caused by collisions between dust grains and supersonic neutral streams in a magnetized shock (C-shock). We demonstrate that, if the ISM in a magnetic field is impacted by a C-shock developed by a supernova explosion in the early phase of massive star-formation ($\lesssim 10$ Myr), AME can be significantly or almost entirely suppressed relative to free-free and thermal dust continuum emission if the grain tensile strength is small enough. This study may shed light on explaining the rare observations of AME from extragalactic star-forming regions preferentially observed from massive star clusters and suggest a scenario of "the rise and fall of AME" in accordance with the temporal evolution of star-forming regions.

arped disk galaxies are classified into two morphologies: S- and U-types. Conventional theories routinely attribute both types to galactic tidal interaction and/or gas accretion, but reproducing of U-types in simulations is extremely challenging. Here we investigate whether both types are governed by the same mechanisms using the most extensive sample of $\sim$8000 nearby (0.02\,$<$\,z\,$<$\,0.06) massive ($M_{*}/M_{\odot}$\,$>$\,$10^9$) edge-on disks from SDSS. We find that U-types show on average bluer optical colors and higher specific star formation rate (sSFR) than S-types, with more strongly warped U-types having higher sSFR. We also find that while the S-type warp properties correlate with the tidal force by the nearest neighbor regardless of the environment, there is no such correlation for U-types in groups/clusters, suggesting a non-tidal environmental could be at play for U-types, such as ram pressure stripping (RPS). Indeed, U-types are more common in groups/clusters than in fields and they have stellar mass, gas fraction, sSFR enhancement and phase-space distribution closely analogous to RPS-induced jellyfish galaxies in clusters. We furthermore show that the stellar disks of most RPS galaxies in the IllustirsTNG simulation are warped in U-shape and bent in opposite direction of stripped gas tails, satisfying theoretical expectations for stellar warps embeded in jellyfishes. We therefore suggest that despite the majority of U-types that live in fields being still less explained, RPS can be an alternative origin for those in groups/clusters.

The detected high-energy pulsars' population is growing in number, and thus, having agile and physically relevant codes to analyze it consistently is important. Here, we update our existing synchro-curvature radiation model by including a better treatment of the particle injection, particularly where the large pitch angle particles dominate the spectra, and by implementing a fast and accurate minimization technique. The latter allows a large improvement in computational cost, needed to test model enhancements and to apply the model to a larger pulsar population. We successfully fit the sample of pulsars with X-ray and $\gamma$-ray data. Our results indicate that, for every emitting particle, the spatial extent of their trajectory where the pitch angle is large and most of the detected X-ray radiation is produced is a small fraction of the light cylinder. We also confirm with this new approach that synchrotron radiation is not negligible for most of the gamma-ray pulsars detected. In addition, with the results obtained, we argue that J0357+3205 and J2055+2539 are MeV-pulsar candidates and are suggested for exhaustive observations in this energy band.

With observational efforts moving from the discovery into the characterisation mode, systematic campaigns that cover large ranges of global stellar and planetary parameters will be needed. We aim to uncover cloud formation trends and globally changing chemical regimes due to the host star's effect on the thermodynamic structure of their atmospheres. We aim to provide input for exoplanet missions like JWST, PLATO, and Ariel, as well as potential UV missions ARAGO, PolStar or POLLUX. Pre-calculated 3D GCMs for M, K, G, F host stars are the input for our kinetic cloud model. Gaseous exoplanets fall broadly into three classes: i) cool planets with homogeneous cloud coverage, ii) intermediate temperature planets with asymmetric dayside cloud coverage, and iii) ultra-hot planets without clouds on the dayside. In class ii),} the dayside cloud patterns are shaped by the wind flow and irradiation. Surface gravity and planetary rotation have little effect. Extended atmosphere profiles suggest the formation of mineral haze in form of metal-oxide clusters (e.g. (TiO2)_N). The dayside cloud coverage is the tell-tale sign for the different planetary regimes and their resulting weather and climate appearance. Class (i) is representative of planets with a very homogeneous cloud particle size and material compositions across the globe (e.g., HATS-6b, NGTS-1b), classes (ii, e.g., WASP-43b, HD\,209458b) and (iii, e.g., WASP-121b, WP0137b) have a large day/night divergence of the cloud properties. The C/O ratio is, hence, homogeneously affected in class (i), but asymmetrically in class (ii) and (iii). The atmospheres of class (i) and (ii) planets are little affected by thermal ionisation, but class (iii) planets exhibit a deep ionosphere on the dayside. Magnetic coupling will therefore affect different planets differently and will be more efficient on the more extended, cloud-free dayside.

Finding evidence of extraterrestrial life would be one of the most profound scientific discoveries ever made, advancing humanity into a new epoch of cosmic awareness. The Venus Life Finder (VLF) missions feature a series of three direct atmospheric probes designed to assess the habitability of the Venusian clouds and search for signs of life and life itself. The VLF missions are an astrobiology-focused set of missions, and the first two out of three can be launched quickly and at a relatively low cost. The mission concepts come out of an 18-month study by an MIT-led worldwide consortium.

Mounting evidence of chemical disequilibria in the Venusian atmosphere has heightened interest in the search for life within the planet's cloud decks. Balloon systems are currently considered to be the superior class of aerial platform for extended atmospheric sampling within the clouds, providing the highest ratio of science return to risk. Balloon-based aerial platform designs depend heavily on payload mass and target altitudes. We present options for constant- and variable-altitude balloon systems designed to carry out science operations inside the Venusian cloud decks. The Venus Life Finder (VLF) mission study proposes a series of missions that require extended in situ analysis of Venus cloud material. We provide an overview of a representative mission architecture, as well as gondola designs to accommodate a VLF instrument suite. Current architecture asserts a launch date of 30 July 2026, which would place an orbiter and entry vehicle at Venus as early as November 29 of that same year.

Venus is known for its extreme surface temperature and its sulfuric acid clouds. But the cloud layers on Venus have similar temperature and pressure conditions to those on the surface of Earth and are conjectured to be a possible habitat for microscopic life forms. We propose a mission concept to explore the clouds of Venus for up to 30 days to evaluate habitability and search for signs of life. The baseline mission targets a 2026 launch opportunity. A super-pressure variable float altitude balloon aerobot cycles between the altitudes of 48 and 60 km, i.e., primarily traversing the lower, middle, and part of the upper cloud layers. The instrument suite is carried by a gondola design derived from the Pioneer Venus Large Probe pressure vessel. The aerobot transmits data via an orbiter relay combined with a direct-to-Earth link. The orbiter is captured into a 6-h retrograde orbit with a low, roughly 170-degree, inclination. The total mass of the orbiter and entry probe is estimated to be 640 kg. An alternate concept for a constant float altitude balloon is also discussed as a lower complexity option compared to the variable float altitude version. The proposed mission would complement other planned missions and could help elucidate the limits of habitability and the role of unknown chemistry or possibly life itself in the Venus atmosphere.

We study the effects of multiple transitions in the vacuum dark energy density on the $H_0$ tension problem. We consider a phenomenological model in which the vacuum energy density undergoes multiple transitions in the early as well as the late universe and compare the model's predictions using the three sets of data from CMB+BAO+SN. The transient dark energy can be either positive (dS-like) or negative (AdS-like). We conclude that a transient late-time AdS-type vacuum energy typically yields the higher value of $H_0$ which can alleviate the $H_0$ tension. In addition, to obtain a value of $H_0$ comparable to the value obtained from the local cosmological measurements the spectral index $n_s$ moves towards its Harrison-Zel'dovich scale invariant value

The solar convection zone rotates differentially, with its equatorial region rotating more rapidly than the polar regions. This form of differential rotation, also observed in many other low-mass stars, is understood to arise when Coriolis effects are stronger than those associated with buoyant driving of the convection. When buoyancy dominates, a so-called antisolar state of differential rotation results, characterized by rapidly-rotating poles and a slow equator. The transition between these two states has been shown to occur when the intensity of these two forces is roughly equal or, equivalently, when the convective Rossby number of the system is unity. Here we consider an alternative view of the transition that relates this phenomenon to convective structure and convective-zone depth. Using a series of 3-D rotating convection-zone simulations, we demonstrate that the solar/antisolar transition occurs when the columnar convective structures characteristic of rotating convection attain a diameter roughly equivalent to the shell depth. When the characteristic convective wavelength exceeds twice the shell depth, we find that the coherent convective structures necessary to sustain an equatorward Reynolds stress are lost, and an antisolar state results. We conclude by presenting a force-balance analysis that relates this geometric interpretation of the transition to the convective-Rossby-number criteria identified in previous studies.

The atmospheres of gas giant planets are thought to be inhomogeneous due to weather and patchy clouds. We present two full nights of coronagraphic observations of the HR 8799 planets using the CHARIS integral field spectrograph behind the SCExAO adaptive optics system on the Subaru Telescope to search for spectrophomometric variability. We did not detect significant variability signals, but placed the lowest variability upper limits for HR 8799 c and d. Based on injection-recovery tests, we expected to have a 50% chance to detect signals down to 10% H-band photometric variability for HR 8799 c and down to 30% H-band variability for HR 8799 d. We also investigated spectral variability and expected a 50% chance to recovery 20% variability in the H/K flux ratio for HR 8799 c. We combined all the data from the two nights to obtain some of the most precise spectra obtained for HR 8799 c, d, and e. Using a grid of cloudy radiative-convective-thermochemical equilibrium models, we found all three planets prefer supersolar metallicity with effective temperatures of ~1100 K. However, our high signal-to-noise spectra show that HR 8799 d has a distinct spectrum from HR 8799 c, possibly preferring more vertically extended and uniform clouds and indicating that the planets are not identical.

We study the spectral properties and accretion flow behavior of an ultraluminous X-ray source M82\,X-1 using {\it NuSTAR} observations. We use the physical two component advective flow (TCAF) model to fit the data and to derive the accretion flow properties of the source. From the model fitted parameters, we found that M82\,X-1 is harboring an intermediate mass black hole at its centre, where the mass varies from $156.04^{+13.51}_{-15.30}$ to $380.96^{28.38}_{-29.76}$ M$_\odot$. The error weighted average mass of the black hole is $273\pm43$ M$_\odot$, which accreted in nearly super-Eddington rate. The Compton cloud was compact with a size of $\sim13 r_g$ and the shock compression ratio had \textcolor{black}{nearly intermediate values except for the epoch four}. These indicate a possible significant mass outflow from the inner region of the disk. The quasi periodic oscillation (QPO) frequencies estimated from the model fitted parameters can reproduce the observed QPOs. The robustness of the model parameters is verified by drawing the confidence contours among them.

Damped Lyman-$\alpha$ Absorber(DLA), or HI 21cm absorber, is an important probe to directly measure the acceleration of spectroscopic velocity $v_\mathrm{S}$ via the Sandage-Loeb(SL) effect. Confined by the shortage of actual DLAs samples and the coarse background radio sources assignment, the detectable amount of Damped Lyman-$\alpha$ Absorption System(DLAS) is ambiguous in most cases. After differing the unmeasurable, global and physical $\ddot{a}$ from the observed and local $\dot{v}_\mathrm{S}$, we make a statistical investigation of the components of DLASs. We use Kernel Density Estimation(KDE) to depict a general redshift distribution of background radio sources via three radio deep survey datasets, CENSORS, LBDS-Hercules and CoNFIG-4, and provide a multi-Gaussian expression. Testing the generation process of DLA redshift number density in literature, we try to make a modified power law fitting of low-redshift($z\lesssim1.65$) DLA preselected by MgII absorption and analysis its defects. Finally, we present a simple DLASs number estimation of FAST, ASKAP and SKA-Mid when considering a blind HI absorption survey with our derived radio number density and the previous DLA one in literature. For comparability, our FAST prediction gives a practical amount of 100, and an optimistic amount of 470, while our latter amount and previous predictions are within an order of magnitude.

Ultra-light primordial black holes (PBHs) in the mass range of $10^{16} - 10^{22}$ g are allowed by current observations to constitute a significant fraction, if not all, of the dark matter in the Universe. In this work, we present limits on ultra-light, non-rotating PBHs which arise from the non-detection of the Hawking radiation signals from such objects in the keV-MeV energy band. Namely, we consider observations from the current-generation missions XMM-Newton and INTEGRAL/SPI and discuss the observational perspectives of the future missions Athena, eXTP, and THESEUS for PBH searches. Based on 3.4 Msec total exposure time XMM-Newton observations of Draco dwarf spheroidal galaxy, we conclude that PBH with masses $\lesssim 10^{16}$ g can not make all dark matter at 95% confidence level. Our ON-OFF-type analysis of >100 Msec of INTEGRAL/SPI data on the Milky Way halo puts significantly stronger constraints. Only $\lesssim 10$% dark matter can be presented by PBHs with masses $\lesssim 3 \cdot 10^{16}$ g while the majority of dark matter can not be presented by PBHs lighter than $7\cdot 10^{16}$ g at 95% confidence level. We discuss the strong impact of systematic uncertainty related to the variations of instrumental and astrophysical INTEGRAL/SPI background on the derived results and estimate its level. We also show that future large-field-of-view missions such as THESEUS/X-GIS will be able to improve the constraints by a factor of $10 - 100$ depending on the level of control under the systematics of these instruments.

Simulating deep solar convection and its coupled mean-field motions is a formidable challenge where few observational results constrain models that suffer from the non-physical influence of the grid resolution. We present hydrodynamic global Implicit Large-Eddy simulations (ILES) of deep solar convection performed with the EULAG-MHD code, and explore the effects of grid resolution on the properties of rotating and non-rotating convection. The results, based on low-order moments and turbulent spectra reveal that convergence could be achieved in non-rotating simulations provided sufficient resolution in the radial direction. The flow is highly anisotropic, with the energy contained in horizontal divergent motions exceeding by more than three orders of magnitude their radial counterpart. By contrast, in rotating simulations the largest energy is in the toroidal part of the horizontal motions. As the grid resolution increases, the turbulent correlations change in such a way that a solar-like differential rotation, obtained in the simulation with the coarser grid, transitions to the anti-solar differential rotation. The reason for this change is the contribution of the effective viscosity to the balance of the forces driving large-scale flows. As the effective viscosity decreases, the angular momentum balance improves, yet the force balance in the meridional direction lessens, favoring a strong meridional flow that advects angular momentum towards the poles. The results suggest that obtaining the correct distribution of angular momentum may not be a mere issue of numerical resolution. Accounting for additional physics, such as magnetism or the near-surface shear layer, may be necessary in simulating the solar interior.

My thesis work aims to study the inter-relation between various physical and chemical conditions in a wide range of astrophysical environments. Our studied regions range from the super-hot regions (i.e., nebular, photon-dominated, or photodissociation regions, diffuse area, through which the lights of the background stars can reach us) to the super-cold regions (i.e., dense molecular clouds, proto-planetary disks, etc. where interstellar dust particles absorb all background visible and ultra-violet lights). The chemical complexity of the interstellar cloud gradually evolves due to the evolution in physical conditions. The dense molecular clouds are the birth sites of star-formation, where a wide variety of complex organic molecules are observed. Dust particles play an essential role in the formation of these complex organic molecules. During the warm-up stage of a star-forming region, the molecules formed during the cold phase start to return to the gas phase by various thermal and non-thermal evaporation processes. These complex molecules again freeze out to the outer part of the proto-planetary disk during the further evolved stage according to their condensation temperature and form the so-called snow-lines. The binding energies of these molecules with the prevailing dust particles play a crucial role in determining the structural information of this disk. Thus the binding energy of the molecules is vital to understand several critical aspects of the star and planet formation processes. In this thesis, I will discuss the chemical complexity obtained in a wide range of the star-forming region and whether this chemical complexity lead to biomolecules in space.

Highly irradiated planets in the hot Neptune desert are usually either small (R < 2 Rearth) and rocky or they are gas giants with radii of >1 Rjup. Here, we report on the intermediate-sized planet TOI-2196 on a 1.2 day orbit around a G-type star discovered by TESS in sector 27. We collected 42 radial velocity measurements with the HARPS spectrograph to determine the mass. The radius of TOI-2196 b is 3.51 +/- 0.15 Rearth, which, combined with the mass of 26.0 +/- 1.3 Mearth, results in a bulk density of 3.31+0.51-0.43 g/cm3. Hence, the radius implies that this planet is a sub-Neptune, although the density is twice than that of Neptune. A significant trend in the HARPS radial velocities points to the presence of a distant companion with a lower limit on the period and mass of 220 days and 0.65 Mjup, respectively, assuming zero eccentricity. The short period of planet b implies a high equilibrium temperature of 1860 +/- 20 K, for zero albedo and isotropic emission. This places the planet in the hot Neptune desert, joining a group of very few planets in this parameter space discovered in recent years. These planets suggest that the hot Neptune desert may be divided in two parts for planets with equilibrium temperatures of > 1800 K: a hot sub-Neptune desert devoid of planets with radii of 1.8-3 Rearth and a sub-Jovian desert for radii of 5-12 Rearth. More planets in this parameter space are needed to further investigate this finding. Planetary interior structure models of TOI-2196 b are consistent with a H/He atmosphere mass fraction between 0.4 % and 3 %, with a mean value of 0.7 % on top of a rocky interior. We estimated the amount of mass this planet might have lost at a young age, and we find that while the mass loss could have been significant, the planet had not changed in terms of character: it was born as a small volatile-rich planet, and it remains one at present.

The detection of coherent X-ray pulsations at ~314 Hz (3.2 ms) classifies MAXI J1957+032 as a fast-rotating, accreting neutron star. We present the temporal and spectral analysis performed using NICER observations collected during the latest outburst of the source. Doppler modulation of the X-ray pulsation revealed the ultra-compact nature of the binary system characterised by an orbital period of ~1 hour and a projected semi-major axis of 14 lt-ms. The neutron star binary mass function suggests a minimum donor mass of 1.7e-2 Msun, assuming a neutron star mass of 1.4 Msun and a binary inclination angle lower than 60 degrees. This assumption is supported by the lack of eclipses or dips in the X-ray light curve of the source. We characterised the 0.5-10 keV energy spectrum of the source in outburst as the superposition of a relatively cold black-body-like thermal emission compatible with the emission from the neutron star surface and a Comptonisation component with photon index consistent with a typical hard state. We did not find evidence for iron K-alpha lines or reflection components.

The formation of molecules in the interstellar medium (ISM) remains a complex and unresolved question in astrochemistry. A group of molecules of particular interest involves the linkage between a -carboxyl and -amine group, similar to that of a peptide bond. The detection of molecules containing these peptide-like bonds in the ISM can help elucidate possible formation mechanisms, as well as indicate the level of molecular complexity available within certain regions of the ISM. Two of the simplest molecules containing a peptide-like bond, formamide (NH2CHO) and acetamide (CH3CONH2), have previously been detected toward the star forming region Sagittarius B2 (Sgr B2). Recently, the interstellar detection of propionamide (C2H5CONH2) was reported toward Sgr B2(N) with ALMA observations at millimeter wavelengths. Yet, this detection has been questioned by others from the same set of ALMA observations as no statistically significant line emission was identified from any uncontaminated transitions. Using the PRrbiotic Interstellar MOlecule Survey (PRIMOS) observations, we report an additional search for C2H5CONH2 at centimeter wavelengths conducted with the Green Bank Telescope. No spectral signatures of C2H5CONH2 were detected. An upper limit for C2H5CONH2 at centimeter wavelengths was determined to be less than 1.8e14 cm-2 and an upper limit to the C2H5CONH2/CH3CONH2 ratio is found to be less than 2.34. This work again questions the initial detection of C2H5CONH2 and indicates that more complex peptide-like structures may have difficulty forming in the ISM or are below the detection limits of current astronomical facilities. Additional structurally related species are provided to aid in future laboratory and astronomical searches.

The anecdotal connection between an interest in science fiction and career aspirations in astrophysics is well established. However strong statistical evidence for such a connection, and a quantitative assessment of its prevalence, has been missing. Here I report the results of two surveys examining the connection between science fiction enthusiasm and astronomical careers - first a case study of the University of Warwick Astronomy and Astrophysics group, carried out in February 2021, and second a larger survey of attendees at the UK National Astronomy Meeting in July 2022. In both surveys, a significant majority of respondents expressed an interest in science fiction. In the larger survey, 93% of UK astronomers (223 of 239 respondents) expressed an interest in science fiction, while 69% (164) stated that it had influenced their life or career choices. This study provides strong statistical evidence for the role of science fiction in influencing the adoption of astronomical careers.

We present follow-up results from the first Fundamental Reference AGN Monitoring Experiment (FRAMEx) X-ray/radio snapshot program of a volume-complete sample of local hard X-ray-selected active galactic nuclei (AGNs). Here, we added 9 new sources to our previous volume-complete snapshot campaign, two of which are detected in the 6 cm Very Long Baseline Array (VLBA) observations. We also obtained deeper VLBA observations for a sample of 9 AGNs not detected by our previous snapshot campaign. We recovered 3 sources with approximately twice the observing sensitivity. In contrast with lower angular resolution Very Large Array (VLA) studies, the majority of our sources continue to be undetected with the VLBA. The sub-parsec radio (6 cm) and X-ray (2-10 keV) emission show no significant correlation, with L_R/L_X ranging from 10^-8 to 10^-4, and the majority of our sample lies well below the fiducial 10^-5 relationship for coronal synchrotron emission. Additionally, our sources are not aligned with any of the proposed "fundamental" planes of black hole activity, which purport to unify black hole accretion in the M_BH-L_X-L_R parameter space. The new detections in our deeper observations suggest that the radio emission may be produced by the synchrotron radiation of particles accelerated in low luminosity outflows. Non-detections may be a result of synchrotron self-absorption at 6 cm in the radio core, similar to what has been observed in X-ray binaries (XRBs) transitioning from the radiatively inefficient state to a radiatively efficient state.

The Pierre Auger and the Telescope Array observatories have measured independent and statistical significant anisotropy in the arrival direction of ultra-high-energy cosmic rays (UHECR). Three hotspot regions with relative excess of events and a dipole signal have been identified in different regions of the sky and energy ranges. In this paper, we investigate the conditions under which these anisotropy signal could be generated by nearby (<23 Mpc) active galactic nuclei (AGN) and/or starburst galaxies (SBG). We studied a wide range of possibilities including injected nuclei (p, He, N, Si, and Fe), three UHECR luminosity proxies and three extragalactic magnetic field models. The results shows that both local AGN and SBG are needed to describe all the anisotropy signal. The contribution of AGN to hotspots and to the generation of the dipole is dominant in most cases. SBG is required only to explain the hotspot measured by the Telescope Array Observatory.

We continue the search for luminous blue variables (LBVs) in Local Volume galaxies in order to study their fundamental parameters. In this paper, we report the discovery of two new LBVs in the dwarf irregular galaxy NGC 1156. Both stars exhibit spectral variability simultaneously with strong brightness variations: $\Delta \text{R}_c = 0.84 \pm 0.23^m$ for J025941.21+251412.2 and $\Delta \text{R}_c = 2.59 \pm 0.10^m$ for J025941.54+251421.8. The bolometric luminosities of the stars are in the range of L$_\text{Bol} \approx (0.8 - 1.6) \times 10^6$ L$_\odot$. These values are corrected for reddening A$_\text{V} \approx 0.9$ and are given for the distance to the galaxy D=7.0$\pm$0.4 Mpc, which we have measured by the TRGB method. Both stars are above the Humphreys-Davidson limit in the region of relatively low temperatures, T$_\text{eff} \lesssim 10$ kK on the temperature-luminosity diagram. J025941.54+251421.8 had a temperature below the hydrogen ionisation threshold at maximum brightness, exhibiting behaviour very similar to that of the known LBV R71 during its 2012 outburst. We have estimated the masses of the detected LBVs and studied the properties of their stellar environment. We discuss our results within the framework of both a single star and a binary system evolution scenario for LBVs.

Gaseous circumbinary disks (CBDs) that are highly inclined to the binary orbit are commonly observed in nature. These disks harbor particles that can reach large mutual inclinations as a result of nodal precession once the gas disk has dissipated. With n-body simulations that include fragmentation we demonstrate that misaligned disks of particles can be efficient progenitors of interstellar asteroids (ISAs). Collisions that take place between particles with large mutual inclinations have large impact velocities which can result in mass ejection, with a wide range of fragment sizes and ejection velocities. We explore the binary parameters for which the majority of the terrestrial planet forming material is ejected rather than accreted into planets. The misalignment required to eject significant material decreases with binary eccentricity. If the distribution of binary eccentricity is uniform and the initial particle CBD orientation relative to the binary orbit is isotropic, about 59% of binaries are more likely to eject the majority of their CBD terrestrial planet disk mass through high velocity body-body collisions rather than retain this material and build terrestrial planets. However, binary--disk interactions during the gas disk phase with non-zero disk viscosity will reduce this fraction. The composition, small size, highly elongated shape, and tumbling motion of `Oumuamua is consistent with ISAs generated by misaligned CBDs.

Context. Blob detection is a common problem in astronomy. One example is in stellar population modelling, where the distribution of stellar ages and metallicities in a galaxy is inferred from observations. In this context, blobs may correspond to stars born in-situ versus those accreted from satellites, and the task of blob detection is to disentangle these components. A difficulty arises when the distributions come with significant uncertainties, as is the case for stellar population recoveries inferred from modelling spectra of unresolved stellar systems. There is currently no satisfactory method for blob detection with uncertainties. Aims. We introduce a method for uncertainty-aware blob detection developed in the context of stellar population modelling of integrated-light spectra of stellar systems. Methods. We develop theory and computational tools for an uncertainty-aware version of the classic Laplacian-of-Gaussians method for blob detection, which we call ULoG. This identifies significant blobs considering a variety of scales. As a prerequisite to apply ULoG to stellar population modelling, we introduce a method for efficient computation of uncertainties for spectral modelling. This method is based on the truncated Singular Value Decomposition and Markov Chain Monte Carlo sampling (SVD-MCMC). Results. We apply the methods to data of the star cluster M54. We show that the SVD-MCMC inferences match those from standard MCMC, but are a factor 5-10 faster to compute. We apply ULoG to the inferred M54 age/metallicity distributions, identifying between 2 or 3 significant, distinct populations amongst its stars.

We continue to search for LBV stars in galaxies outside the Local Group. In this work, we have investigated four luminous stars in NGC 4449. Multiple spectral observations carried out for J122810.94+440540.6, J122811.70+440550.9, and J122809.72+440514.8 revealed the emission features in their spectra that are characteristic of LBVs. Photometry showed noticeable brightness changes of J122809.72+440514.8 ($\Delta I=0.69\pm0.13^m$) and J122817.83+440630.8 ($\Delta R=2.15\pm0.13^m$), while the variability of J122810.94+440540.6 and J122811.70+440550.9 does not exceed $0.3^m$ regardless of the filter. We have obtained estimates of the interstellar reddening, photosphere temperatures, and bolometric luminosities $\log(\text{L}_\text{Bol}/\text{L}_{\odot}) \approx 5.24-6.42$. Using the CMFGEN code, we have modelled the spectrum of the cold state of J122809.72+440514.8 ($T_{\text{eff}}=9300\,$K) and have obtained possible value of the mass loss rate $\dot{M} = 5.2\times10^{-3}\,M_{\odot}\,yr^{-1}$. Based on the observational properties, J122809.72+440514.8 and J122817.83+440630.8 were classified as LBVs, while the other two stars were classified as LBV candidates or B[e]-supergiants candidates

We present Chandra X-ray Observatory observations and Space Telescope Imaging Spectrograph spectra of NGC 5972, one of the 19 "Voorwerpjes" galaxies. This galaxy contains an Extended Emission Line Region (EELR) and an arc-second scale nuclear bubble. NGC 5972 is a faded AGN, with EELR luminosity suggesting a 2.1 dex decrease in L$_{\textrm{bol}}$ in the last $\sim5\times10^{4}$ yr. We investigate the role of AGN feedback in exciting the EELR and bubble given the long-term variability and potential accretion state changes. We detect broadband (0.3-8 keV) nuclear X-ray emission coincident with the [OIII] bubble, as well as diffuse soft X-ray emission coincident with the EELR. The soft nuclear (0.5-1.5 keV) emission is spatially extended and the spectra are consistent with two APEC thermal populations ($\sim$0.80,$\sim$0.10 keV). We find a bubble age >2.2 Myr, suggesting formation before the current variability. We find evidence for efficient feedback with L$_{\textrm{kin}}/L_{\textrm{bol}}\sim0.8\%$, which may be overestimated given the recent L$_{\textrm{bol}}$ variation. Kinematics suggest an out-flowing 300 km s$^{-1}$ high-ionization [OIII]-emitting gas which may be the line of sight component of a $\sim$780 km s$^{-1}$ thermal X-ray outflow capable of driving strong shocks that could photoionize the precursor material. We explore possibilities to explain the overall jet, radio lobe and EELR misalignment including evidence for a double SMBH which could support a complex misaligned system.

The X-ray variability of active galactic nuclei (AGN) carries crucial information about the X-ray radiation mechanism. We performed a systematic study of the X-ray short-term (1-100 ks timescale) variability for a large sample of 78 Seyferts with 426 deep XMM-Newton observations. In this paper, we present the time-averaged spectra and rms spectra for the entire sample, which show a variety of properties. Based on the spectral shape, we divide the rms spectra into five subtypes and the time-averaged spectra into four subtypes. The most common shape of the rms spectra is concave-down where the rms peaks at $\sim$ 1 keV. We find that different sources can show similar time-averaged spectra and rms spectra. However, there is no one-to-one mapping between the subtypes of the time-averaged spectra and rms spectra, as similar time-averaged spectra can be accompanied by different rms spectra, and vice versa. This is likely because different physical mechanisms can produce similar rms spectra. For every subtype of the time-averaged spectra, we report its preferred subtypes of the rms spectra in both low- and high-frequency bands. We also compare the statistical properties for different subtypes, such as the black hole mass and Eddington ratio. Finally, we investigate the rms in the Fe K$\alpha$ line regime and find that those with a broad and extended red-wing profile tend to show stronger variability than those showing a narrow or relatively symmetric profile. Our results demonstrate the necessity of performing joint spectral and variability modeling in order to understand the mechanism of the X-ray emission in AGN. All of the rms spectra have been made publicly available.

The data of cosmic ray NEVOD-DECOR experiment on the investigation of inclined muon bundles for a long time period (May 2012 - March 2021) are presented. The analysis showed that the observed intensity of muon bundles at primary cosmic ray energies of about 1 EeV and higher can be compatible with the expectation in frame of widely used hadron interaction models only under the assumption of an extremely heavy mass composition. This conclusion is consistent with data of several experiments on investigations of muon content in air showers, but contradicts the available fluorescence data on Xmax which favor a light mass composition at these energies. In order to clarify the nature of the "muon puzzle", investigations of the muon bundle energy deposit in the detector material were carried out. For the first time, experimental estimates of the average energy of muons in the bundles of inclined air showers initiated by primary particles with energies from 10 to 1000 PeV have been obtained.

We assess the robustness of $\Lambda$CDM results from the full-shape analysis of BOSS power spectrum using the one-loop prediction of the Effective Field Theory of Large-Scale Structure (EFTofLSS). The public likelihoods PyBird and CLASS-PT lead to results in agreement only at the $1\sigma$ level, despite the fact that they are derived from the same BOSS dataset and theory model. We perform a thorough comparison of the various analyses choices made between the two pipelines, and identify that the differences come from the choice of prior on the EFT parameters, dubbed "West-coast" (WC) and "East-coast" (EC) prior, respectively associated to PyBird and CLASS-PT. In particular, because posteriors are non-Gaussian, projection effects from the marginalization over the EFT parameters shift the posterior mean of the cosmological parameters with respect to the best-fit up to $1\sigma$ in the WC prior and up to $2\sigma$ in the EC prior. We quantify that best-fit cosmological parameters extracted from BOSS given the two prior choices are consistent at $\sim 1\sigma$. The consistency improves to $\sim 0.5\sigma$ when doubling the prior widths. While this reveals that current EFT analyses are subject to prior effects, we show that cosmological results obtained in combination with CMB, or from forthcoming large-volume data, are less sensitive to those effects. In addition, we investigate differences between BOSS measurements. We find broad agreements across all pre-reconstructed measurements considered ($<0.6\sigma$), but the two available BOSS post-reconstructed measurements in Fourier space, once combined with the EFT full-shape analysis, lead to discrepant Hubble parameter $H_0$ at $\sim 0.9\sigma$. Given the various effects we discuss, we argue that the clustering amplitude $\sigma_8$ measured with BOSS is not in statistical tension with that inferred from Planck under $\Lambda$CDM.

Analyses of the full shape of BOSS DR12 power spectrum using the one-loop prediction from the Effective Field Theory of Large-Scale Structure (EFTBOSS) have led to new constraints on extensions to the $\Lambda$CDM model, such as Early Dark Energy (EDE) which has been suggested as a resolution to the "Hubble tension". In this paper, we re-assess the constraining power of the EFTBOSS on EDE in light of a correction to the normalization of BOSS window functions. Overall we find that constraints from EFTBOSS on EDE are weakened, and represent a small change compared to constraints from Planck and the conventional BAO/$f\sigma_8$ measurements. The combination of Planck data with EFTBOSS provides a bound on the maximal fractional contribution of EDE $f_{\rm EDE}<0.083$ at 95% C.L. (compared to $<0.054$ with the incorrect normalization, and $<0.088$ without full-shape data) and the Hubble tension is reduced to $2.1\sigma$. However, the more extreme model favored by an analysis with just data from the Atacama Cosmology Telescope is disfavored by the EFTBOSS data. We also show that the updated Pantheon+ Type Ia supernova analysis can slightly increase the constraints on EDE. Yet, the inclusion of the SN1a magnitude calibration by SH0ES strongly increases the preference for EDE to above $5\sigma$, yielding $f_{\rm EDE}\sim 0.12^{+0.03}_{-0.02}$ around the redshift $z_c=4365^{+3000}_{-1100}$. Our results demonstrate that EFTBOSS data (alone or combined with Planck data) do not exclude the EDE resolution of the Hubble tension.

Galaxy formation and evolution involves a variety of effectively stochastic processes that operate over different timescales. The Extended Regulator model provides an analytic framework for the resulting variability (or `burstiness') in galaxy-wide star formation due to these processes. It does this by relating the variability in Fourier space to the effective timescales of stochastic gas inflow, equilibrium, and dynamical processes influencing GMC creation and destruction using the power spectral density (PSD) formalism. We use the connection between the PSD and auto-covariance function (ACF) for general stochastic processes to reformulate this model as an auto-covariance function, which we use to model variability in galaxy star formation histories (SFHs) using physically-motivated Gaussian Processes in log SFR space. Using stellar population synthesis models, we then explore how changes in model stochasticity can affect spectral signatures across galaxy populations with properties similar to the Milky Way and present-day dwarfs as well as at higher redshifts. We find that, even at fixed scatter, perturbations to the stochasticity model (changing timescales vs overall variability) leave unique spectral signatures across both idealized and more realistic galaxy populations. Distributions of spectral features including H$\alpha$ and UV-based SFR indicators, H$\delta$ and Ca-H,K absorption line strengths, D$_n$(4000) and broadband colors provide testable predictions for galaxy populations from present and upcoming surveys with Hubble, Webb \& Roman. The Gaussian process SFH framework provides a fast, flexible implementation of physical covariance models for the next generation of SED modeling tools. Code to reproduce our results can be found at https://github.com/kartheikiyer/GP-SFH

Stellar-mass primordial black holes (PBHs) from the early Universe can directly contribute to the gravitational wave (GW) events observed by LIGO, but can only comprise a subdominant component of the dark matter (DM). The primary DM constituent will generically form massive halos around seeding stellar-mass PBHs. We demonstrate that gravitational lensing of sources at cosmological ($\gtrsim$Gpc) distances can directly explore DM halo dresses engulfing PBHs, challenging for lensing of local sources in the vicinity of Milky Way. Strong lensing analysis of fast radio bursts detected by CHIME survey already starts to probe parameter space of dressed stellar-mass PBHs, and upcoming searches can efficiently explore dressed PBHs over $\sim 10-10^3 M_{\odot}$ mass-range and provide a stringent test of the PBH scenario for the GW events. Our findings establish a general test for a broad class of DM models with stellar-mass PBHs, including those where QCD axions or WIMPs comprise predominant DM. The results open a new route for exploring dressed PBHs with various types of lensing events at cosmological distances, such as supernovae and caustic crossings.

With planned space-based and 3rd generation ground-based gravitational wave detectors (LISA, Einstein Telescope, Cosmic Explorer), and proposed DeciHz detectors (DECIGO, Big Bang Observer), it is timely to explore statistical cosmological tests that can be employed with the forthcoming plethora of data, $10^4-10^6$ mergers a year. We forecast the combination of the standard siren measurement with the weak lensing of gravitational waves from binary mergers. For 10 years of 3rd generation detector runtime, this joint analysis will constrain the dark energy equation of state with marginalised $1\sigma$ uncertainties of $\sigma(w_0)$~0.005 and $\sigma(w_a)$~0.04. This is comparable to or better than forecasts for future galaxy/intensity mapping surveys, and better constraints are possible when combining these and other future probes with gravitational waves. We find that combining mergers with and without an electromagnetic counterpart helps break parameter degeneracies. In the post-LISA era, we demonstrate for the first time how merging binaries could achieve a precision on the sum of neutrino masses of $\sigma(\Sigma m_{\nu})$~0.04 eV, and ~percent or sub-percent precision also on curvature, dark energy, and other parameters, independently from other probes. Finally, we demonstrate how the cosmology dependence in the redshift distribution of mergers can be exploited to improve dark energy constraints if the cosmic merger rate is known, instead of relying on measured distributions as is standard in cosmology. In the coming decades gravitational waves will become a formidable probe of both geometry and large scale structure.

Dark matter direct detection experiments are designed to look for the scattering of dark matter particles that are assumed to move with virial velocities $\sim 10^{-3}$. At these velocities, the energy deposition in the detector is large enough to cause ionization/scintillation, forming the primary class of signatures looked for in such experiments. These experiments are blind to a large class of dark matter models where the dark matter has relatively large scattering cross-sections with the standard model, resulting in the dark matter undergoing multiple scattering with the atmosphere and the rock overburden, and thus slowing down considerably before arriving at underground detectors. We propose to search for these kinds of dark matter by looking for the anomalous heating of a well shielded and sensitive calorimeter. In this detector concept, the dark matter is thermalized with the rock overburden but is able to pierce through the thermal shields of the detector causing anomalous heating. Using the technologies under development for EDELWEISS and SuperCDMS, we estimate the sensitivity of such a calorimetric detector. In addition to models with large dark matter - standard model interactions, these detectors also have the ability to probe dark photon dark matter.

We examine the sensitivity of the Large Hadron Collider (LHC) to light lepton portal dark matter with its mass below 10 GeV. The model features an extra doublet scalar field and singlet Dirac dark matter, which have Yukawa interactions with left-handed leptons. To correctly produce the dark matter abundance via the thermal freeze-out, a large mass splitting among the extra scalars is required, thus providing a light neutral scalar below ${\cal O}(10)$GeV and heavy neutral and charged scalars at the electroweak scale. In this paper, we focus on the electroweak pair-production of the extra scalars with subsequent model-specific scalar decays and evaluate the current constraints with the LHC Run 2 data and the discovery potential at the High Luminosity LHC (HL-LHC). It turns out that a large part of the theoretically allowed parameter space can be tested at the HL-LHC by taking into account complementarity between slepton searches and mono-$Z$ plus missing transverse energy search. We also discuss same-sign charged scalar production as a unique prediction of the model, and the implication of the collider searches in the thermal dark matter scenario.

We investigate the impact of interactions between hidden sectors and the discovered Higgs boson $h_{125}$, allowing for additional invisible decay channels of $h_{125}$. We perform $\chi^2$-fits to the measurements of the Higgs-boson cross sections as a function of the invisible branching ratio and different combinations of coupling modifiers, where the latter quantify modifications of the couplings of $h_{125}$ compared to the predictions of the Standard Model. We present generic results in terms of exclusion limits on the coupling modifiers and the invisible branching ratio of $h_{125}$. Additionally, we apply our results to a variety of concrete model realizations containing a hidden sector: dark matter within Higgs- and singlet-portal scenarios, models featuring (pseudo) Nambu-Goldstone bosons and two Higgs doublet extensions. One of the main conclusions of our work is that in a wide class of models the indirect constraints resulting from the measurements of the cross sections of $h_{125}$ provide substantially stronger constraints on the invisible Higgs-boson branching ratio compared to the direct limits obtained from searches for the invisible decay of $h_{125}$. However, we demonstrate that the presence of an invisible decay mode of $h_{125}$ can also open up parameter space regions which otherwise would be excluded as a result of the indirect constraints. As a byproduct of our analysis, we show that in light of the new results from the LZ collaboration a fermionic DM candidate within the simplest Higgs-portal scenario is completely ruled out.

We present what we know on nucleosynthesis in the Universe and hypotheses that have been made in this regard. A brief description of the Universe's evolution during its different stages is offered, indicating which are the periods and mechanisms of element formation. A critical prospective on future research is formulated to validate, modify, or reject the hypotheses formulated. These will involve joint observations that encompass finer measurements of cosmic background radiation, galaxy clusters, and gravitational waves produced by neutron star collisions. The information thus obtained will be combined with restrictions given by theoretical models. Perhaps many current doubts will be clarified, but new questions will arise.

In this work we investigate the quantitative effects of the misalignment kinetic axion on $R^2$ inflation. Due to the fact that the kinetic axion possesses a large kinetic energy which dominates its potential energy, during inflation its energy density redshifts as stiff matter fluid and evolves in a constant-roll way, making the second slow-roll index to be non-trivial. At the equations of motion level, the $R^2$ term dominates the evolution, thus the next possible effect of the axion could be found at the cosmological perturbations level, via the second slow-roll index which is non-trivial. As we show, the latter elegantly cancels from the observational indices, however, the kinetic axion extends the duration of the inflationary era to an extent that it may cause a 15$\%$ decrease in the tensor-to-scalar ratio of the vacuum $R^2$ model. This occurs because as the $R^2$ model approaches its unstable quasi-de Sitter attractor in the phase space of $F(R)$ gravity due to the $\langle R^2 \rangle $ fluctuations, the kinetic axion dominates over the $R^2$ inflation and in effect the background equation of state is described by a stiff era, or equivalently a kination era, different from the ordinary radiation domination era. This in turn affects the duration of the inflationary era, increasing the $e$-foldings number up to $5$ $e$-foldings in some cases, depending on the reheating temperature, which in turn has a significant quantitative effect on the observational indices of inflation and especially on the tensor-to-scalar ratio.

An axion cloud surrounding a supermassive black hole can be naturally produced through the superradiance process. Its existence can be examined by the axion induced birefringence effect. It predicts an oscillation of the electric vector position angle of linearly polarized radiations. Stringent constraints of the existence of the axion in a particular mass window has been obtained based on the recent Event Horizon Telescope measurement on M87$^\star$. The future Very-Long-Baseline Interferometry (VLBI) observations will be able to measure the vicinity of many supermassive black holes, thus it opens the possibility to search for the existence of axions in a wide mass regime. In this paper, we study how different black hole properties and accretion flows influence the signatures of the axion induced birefringence. We include the impacts of black hole inclination angles, spins, magnetic fields, plasma velocity distributions, the thickness of the accretion flows. We pay special attention to characterize the washout effects induced by the finite thickness of the accretion flows and the lensed photons. Based on this study, we give prospects on how to optimize the axion search using future VLBI observations, such as the next-generation Event Horizon Telescope, to further increase the sensitivity.

Dark matter may only interact with the visible sector through mediators which are more massive than the inflaton, such as those at the Planck scale or the grand unification scale. In such a scenario, the dark matter is mainly produced out of equilibrium during the period of reheating, often referred to as UV freeze-in. We evaluate the abundance of the dark matter generated from bremsstrahlung off the inflaton decay products assuming no direct coupling between the inflaton and the dark matter. This process generally dominates the production of dark matter for low reheating temperatures where the production through the annihilations of particle in the thermal plasma becomes inefficient. We find that the bremsstrahlung process dominates for reheating temperatures $T_{\mathrm{RH}} \lesssim 10^{10}$ GeV, and produces the requisite density of dark matter for a UV scale $\simeq 10^{16}$ GeV. As examples, we calculate numerically the yield of the dark matter bremsstrahlung through gravitation and dimension-6 vector portal effective interactions.