Papers for Thursday, Feb 03 2022

The recent results of the BICEP and Keck collaborations have put stringent bounds on many inflationary models, including some well-motivated ones. This is certainly the case when gravity remains described by Einstein's theory up to the inflationary scale, but can be avoided by introducing quadratic-in-curvature terms that are effective at that scale. Recently it has also been shown that these terms can UV complete gravity respecting stability and unitarity. Here the predictions of such quadratic gravity are computed and compared with the BICEP/Keck constraints by focusing on some of the inflationary scenarios that are best-motivated from the particle physics point of view and are already ruled out in Einstein gravity: (critical) Higgs inflation and natural inflation. The first scenario can be considered as the most economical option as the inflaton is identified with the only known elementary scalar field in the Standard Model and the near criticality of the Standard Model is used to remain in the perturbative regime. In the second one a pseudo-Nambu-Goldstone boson contributes to the inflationary dynamics and its potential is naturally flat. It is shown that in both scenarios one can restore the agreement with the observational constraints in quadratic gravity.

We present a detailed spectral and imaging analysis of the central $15''$ radius ($\sim 7.5 \text{ kpc}$) region of the merger galaxy NGC 6240 that makes use of all the available \textit{Chandra}-ACIS data ($0.3 - 3 \text{ keV}$ effective exposure of $\sim 190 \text{ ks}$). This region shows extended X-ray structures with lower energy counterparts imaged in CO, [O III] and H$\alpha$ line emission. We find both photo-ionized phases of possible nuclear excitation and thermal shock-excited emission in the different large-scale components: the north-west "loop" detected in H$\alpha$, the region surrounding the two nuclei, the large outflow region to the north-east detected in [O III], and the southern X-ray extensions. The latter could be the ionization cone of the northern nucleus, with the N counterpart being obscured by the galaxy disk. The radial distribution of the X-ray surface brightness suggests a confined hot interstellar medium at $r < 2.5 \text{ kpc}$, with a free-flowing wind at larger radii; if the confinement is magnetic, we estimate B-field values of $\sim 100\,\mu \text{G}$ , similar to those measured in the halo of M82. The thermal gas of the extended halo at $kT \sim 1 \text{ keV}$ absorbs soft X-rays from the AGN, but not the extreme ultraviolet radiation leading to a rapid increase in $F_{\text{[O III]}}/F_{\text{X}}$ beyond $\sim 3 \text{ kpc}$. The $\alpha$ element to Fe abundance ratios of the thermal components in the different regions of the extended X-ray emission are generally compatible with SNe II yields, confirming the importance of the active star formation in NGC 6240.

The cold interstellar medium (ISM) plays a central role in the galaxy evolution process. It is the reservoir that fuels galaxy growth via star formation, the repository of material formed by these stars, and a sensitive tracer of internal and external processes that affect entire galaxies. Consequently, significant efforts have gone into systematic surveys of the cold ISM of the galaxies in the local Universe. This review discusses the resulting network of scaling relations connecting the atomic and molecular gas masses of galaxies with their other global properties (stellar masses, morphologies, metallicities, star formation activity...), and their implications for our understanding of galaxy evolution. Key take-home messages are as follows: (1) From a gas perspective, there are three main factors that determine the star formation rate of a galaxy: the total mass of its cold ISM, how much of that gas is molecular, and the rate at which any molecular gas is converted into stars. All three of these factors vary systematically across the local galaxy population. (2) The shape and scatter of both the star formation main sequence and the mass-metallicity relation are deeply linked to the availability of atomic and molecular gas. (3) Future progress will come from expanding our exploration of scaling relations into new parameter space (in particular the regime of dwarf galaxies), better connecting the cold ISM of large samples of galaxies with the environment that feeds them (the circumgalactic medium in particular), and understanding the impact of these large scales on the efficiency of the star formation process on molecular cloud scales.

Measuring the abundances of neutron-capture elements in Galactic disk stars is an important part of understanding key stellar and galactic processes. In the optical wavelength regime a number of different neutron-capture elements have been measured, however from the infrared H-band only the s-process dominated element cerium has been accurately measured for a large sample of disk stars. The more r-process dominated element ytterbium has only been measured in a small subset of stars so far. In this study we aim to measure the ytterbium (Yb) abundance of local disk giants using the Yb II line at $\lambda_\text{air}$=16498\AA. We also compare the resulting abundance trend with Ce and Eu abundances for the same stars to analyse the s- and r-process contributions. We analyse 30 K-giants with high-resolution H-band spectra using spectral synthesis. The very same stars have already been analysed using high-resolution optical spectra using the same method, but the abundance of Yb was not possible to determine from those spectra due to blending issues for stars with [Fe/H]>-1. In this present analysis, we utilise the stellar parameters determined from the optical analysis. We determined the Yb abundances with an estimated uncertainty for [Yb/Fe] of 0.1 dex. From comparison, the trend of [Yb/Fe] follows closely the [Eu/Fe] trend and has clear s-process enrichment in identified s-rich stars. From the comparison, both the validity of the Yb abundances are ensured, and the theoretical prediction of a roughly 40/60 s-/r-process contribution to Yb's origin is supported. These results show that with a careful and detailed analysis of infrared spectra, reliable Yb abundances can be derived for a wider sample of cooler giants in the range -1.1<[Fe/H]<0.3. This is promising for further studies of the production of Yb and for the r-process channel, key for Galactochemical evolution, in the infrared.

Image plane beam combination in optical interferometers multiplexes the interference fringes from multiple baselines onto a single detector. The beams of starlight are arranged in a non-redundant pattern at the entrance of the combiner so that the signal from each baseline can be separated from one another in the frequency domain. If the signals from different baselines overlap in the frequency domain, this can give rise to a systematic error in the fringe measurements known as baseline crosstalk. In this paper we quantify crosstalk arising from the combination of atmospheric seeing and beam propagation over distances of order hundreds of metres. We find that in idealised conditions atmospheric wavefront errors and beam propagation do not contribute to crosstalk. However, when aperture stops are included in the optical beam train we observe that wavefront errors can result in squared visibility errors arising from crosstalk as high as $\Delta V^{2} = 6.6\times10^{-3}$ under realistic observing conditions.

Distortions of the observed cosmic microwave background imprinted by the Sunyaev-Zel'dovich effect toward massive galaxy clusters due to inverse Compton scattering of microwave photons by high-energy electrons provide a direct measurement of the microwave background temperature at redshifts from 0 to 1. Some additional background temperature estimates exist at redshifts from 1.8 to 3.3 based on molecular and atomic line excitation temperatures in quasar absorption line systems, but are model dependent. To date, no deviations from the expected (1+z) scaling behavior of the microwave background temperature have been seen, but the measurements have not extended deeply into the matter-dominated era of the universe at redshifts z>3.3. Here we report the detection of sub-millimeter line absorption from the water molecule against the cosmic microwave background at z=6.34 in a massive starburst galaxy, corresponding to a lookback time of 12.8 Gyr. Radiative pumping of the upper level of the ground-state ortho-H2O(110-101) line due to starburst activity in the dusty galaxy HFLS3 results in a cooling to below the redshifted microwave background temperature, after the transition is initially excited by the microwave background. The strength of this effect implies a microwave background temperature of 16.4-30.2 K (1-sigma range) at z=6.34, which is consistent with a background temperature increase with redshift as expected from the standard CDM cosmology.

The IceCube experiment discovered PeV-energy neutrinos originating beyond our Galaxy with an energy flux that is comparable to that of TeV-energy gamma rays and EeV-energy cosmic rays. Neutrinos provide the only unobstructed view of the cosmic accelerators that power the highest energy radiation reaching us from the universe. We will review the rationale for building kilometer-scale neutrino detectors that led to the IceCube project, which transformed a cubic kilometer of deep transparent natural Antarctic ice into a neutrino telescope of such a scale. We will summarize the results from the first decade of operations: the status of the observations of cosmic neutrinos and of their first identified source, the supermassive black hole TXS 0506+056. Subsequently, we will introduce the phenomenology associated with cosmic accelerators in some detail. Besides the search for the sources of Galactic and extragalactic cosmic rays, the scientific missions of IceCube and similar instruments under construction in the Mediterranean Sea and Lake Baikal include the observation of Galactic supernova explosions, the search for dark matter, and the study of neutrinos themselves. This review resulted from notes created for summer school lectures and should be accessible to nonexperts.

We present a detailed analysis of the X-ray, infrared, and carbon monoxide (CO) emission for a sample of 35 local ($z \leq 0.15$), active ($L_X \geq 10^{42}$ erg s$^{-1}$) galaxies. Our goal is to infer the contribution of far-ultraviolet (FUV) radiation from star formation (SF), and X-ray radiation from the active galactic nuclei (AGN), respectively producing photodissociation regions (PDRs) and X-ray dominated regions (XDRs), to the molecular gas heating. To this aim, we exploit the CO spectral line energy distribution (CO SLED) as traced by Herschel, complemented with data from single-dish telescopes for the low-J lines, and high-resolution ALMA images of the mid-J CO emitting region. By comparing our results to the Schmidt-Kennicutt relation, we find no evidence for AGN influence on the cold and low-density gas on kpc-scales. On nuclear (r = 250 pc) scales, we find weak correlations between the CO line ratios and either the FUV or X-ray fluxes: this may indicate that neither SF nor AGN radiation dominates the gas excitation, at least at r = 250 pc. From a comparison of the CO line ratios with PDR and XDR models, we find that PDRs can reproduce observations only in presence of extremely high gas densities ($n > 10^5$ cm$^{-3}$). In the XDR case, instead, the models suggest moderate densities ($n \approx 10^{2-4}$ cm$^{-3}$). We conclude that a mix of the two mechanisms (PDR for the mid-J, XDR or possibly shocks for the high-J) is necessary to explain the observed CO excitation in active galaxies.

Features such as photon rings, jets, or hot. spots leave particular topological signatures in a black hole image. As such, topological data analysis can be used to characterize images resulting from high resolution observations (synthetic or otherwise) of black holes in the electromagnetic sector. We demonstrate that persistent homology allows for this characterization to be made automatically by counting the number of connected components and one-dimensional holes. Further, persistent homology also allows for the distance between connected components or diameter of holes to be extracted from the image. In order to apply persistent homology on synthetic black hole images, we also introduce metronization, a new algorithm to prepare black hole images into a form that is suitable for topological analysis.

The induced gravitational wave (GW) background from enhanced primordial scalar perturbations is one of the most promising observational consequences of primordial black hole (PBH) formation from inflation. We investigate the induced GW spectrum $\Omega_{\textrm{IGW}}$ from single-field inflation in the general ultra-slow-roll (USR) framework, restricting the peak frequency band to be inside $10^{-3}$-$1$ Hz and saturating PBH abundance to comprise all dark matter (DM) in the ultralight asteroid-mass window. By invoking successful baryogenesis driven by USR inflation, we verify the viable parameter space for the specific density ratio between baryons and PBH DM observed today, the so-called "cosmic coincidence." We show that the cosmic coincidence requirement bounds the spectral index $n_{\rm UV}$ in the high frequency limit, $\Omega_{\textrm{IGW}}(f\gg 1)\propto f^{-2n_{\rm UV}}$, into $0 < n_{\rm UV} < 1$, which implies that PBHs with masses in the range of $10^{-16}$- $3.6 \times 10^{-15} M_\odot$ could have been excluded by the non-detection of a stochastic GW background at LIGO and Virgo. The induced GW background sourced from the remaining mass window, $10^{-14}$-$10^{-12} M_\odot$, can be tested by upcoming Advanced LIGO and Virgo data and next generation experiments such as LISA, Einstein Telescope and DECIGO.

With the growing number of spectroscopic observations and observational platforms capable of exoplanet atmospheric characterization, there is a growing need for analysis techniques that can distill information about a large population of exoplanets into a coherent picture of atmospheric trends expressed within the statistical sample. In this work, we develop a Hierarchical Bayesian Atmospheric Retrieval (HBAR) model to infer population-level trends in exoplanet atmospheric characteristics. We demonstrate HBAR on the case of inferring a trend in atmospheric CO2 with incident stellar flux, predicted by the presence of a functioning carbonate-silicate weathering negative feedback cycle, an assumption upon which all calculations of the habitable zone (HZ) rest. Using simulated transmission spectra and JWST-quality observations of rocky planets with H2O, CO2, and N2 bearing atmospheres, we find that the predicted trend in CO2 causes subtle differences in the spectra of order 10 ppm in the 1-5 um range, underscoring the challenge inherent to testing this hypothesis. In the limit of highly precise data (100 stacked transits per planet), we show that our HBAR model is capable of inferring the population-level parameters that characterize the trend in CO2, and we demonstrate that the null hypothesis and other simpler trends can be rejected at high confidence. Although we find that this specific empirical test of the HZ may be prohibitively challenging in the JWST era, the HBAR framework developed in this work may find a more immediate usage for the analysis of gas giant spectra observed with JWST, Ariel, and other upcoming missions.

Gas mass is a fundamental quantity of protoplanetary disks that directly relates to their ability to form planets. Because we are unable to observe the bulk H$_2$ content of disks directly, we rely on indirect tracers to provide quantitative mass estimates. Current estimates for the gas masses of the observed disk population in the Lupus star-forming region are based on measurements of isotopologues of CO. However, without additional constraints, the degeneracy between H$_2$ mass and the elemental composition of the gas leads to large uncertainties in such estimates. Here we explore the gas compositions of seven disks from the Lupus sample representing a range of CO-to-dust ratios. With Band 6 and 7 ALMA observations, we measure line emission for HCO$^+$, HCN, and N$_2$H$^+$. We find a tentative correlation among the line fluxes for these three molecular species across the sample, but no correlation with $^{13}$CO or sub-mm continuum fluxes. For the three disks where N$_2$H$^+$ is detected, we find that a combination of high disk gas masses and sub-interstellar C/H and O/H are needed to reproduce the observed values. We find increases of $\sim$10-100$\times$ previous mass estimates are required to match the observed line fluxes. This study highlights how multi-molecular studies are essential for constraining the physical and chemical properties of the gas in populations of protoplanetary disks and that CO isotopologues alone are not sufficient for determining the mass of many observed disks.

Gravitational waves provide a novel and powerful way to test astrophysical models of compact objects, early universe processes, beyond the Standard Model particle physics, dark matter candidates, Einstein's theory of General Relativity and extended gravity models, and even quantum gravity candidate theories. A short introduction to the gravitational-wave background and the method we are using to detect it will be presented. Constraints on various astrophysical/cosmological models from the non-detectability of the gravitational-wave background will be discussed. Gravitational waves from transients will be highlighted and their physical implications will be summarised.

The surface abundances of the light elements lithium (Li) and beryllium (Be) reveal information about the physical processes taking place in stellar interiors. The investigation of the amount of these two elements in stars in open clusters shows the effect of age on those mechanisms. We have obtained spectra of both Li and Be in main-sequence stars in NGC 752 at high spectral resolution and high signal-to-noise ratios with HIRES on the Keck I telescope. In order to make meaningful comparisons with other clusters, we have determined the stellar parameters on a common scale. We have found abundances of Li and Be by spectral synthesis techniques. NGC 752 is twice the age of the well-studied Hyades cluster. We find that 1) The Li dip centered near 6500 K is wider in NGC 752, having expanded toward cooler temperatures; 2) The Be dip is deeper in the older NGC 752; 3) The Li "peak" near 6200 K is lower by about 0.3 dex; 4) Although there is little Be depletion in the cooler stars, it is possible that Be may be lower in NGC 752 than in the Hyades; 5) The Li content in both clusters declines with decreasing temperature, but there is less Li in NGC 752 at a given temperature by $\sim$0.4 dex. These differences are consistent with the transport of the light-element nuclei below the surface convection zone as predicted by theory. That connection to rotational spin-down is indicated by the pattern of rotation with temperature in the two clusters.

Artificially increasing the luminosity and the thermal diffusivity of a model is a common tactic adopted in hydrodynamical simulations of stellar convection. In this work, we analyse the impact of these artificial modifications on the physical properties of stellar interiors and specifically on internal gravity waves. We perform two-dimensional simulations of solar-like stars with the MUSIC code. We compare three models with different luminosity enhancement factors to a reference model. The results confirm that properties of the waves are impacted by the artificial enhancement of the luminosity and thermal diffusivity. We find that an increase in the stellar luminosity yields a decrease in the bulk convective turnover timescale and an increase in the characteristic frequency of excitation of the internal waves. We also show that a higher energy input in a model, corresponding to a larger luminosity, results in higher energy in high frequency waves. Across our tests with the luminosity and thermal diffusivity enhanced together by up to a factor of 104, our results are consistent with theoretical predictions of radiative damping. Increasing the luminosity also has an impact on the amplitude of oscillatory motions across the convective boundary. One must use caution when interpreting studies of internal gravity waves based on hydrodynamical simulations with artificially enhanced luminosity.

Our ability to calibrate current kilometer-scale interferometers can potentially confound the inference of astrophysical signals. Current calibration uncertainties are well described by a Gaussian process. I exploit this description to analytically examine the impact of calibration uncertainty. I derive closed-form expressions for the conditioned likelihood of the calibration error given the observed data and an astrophysical signal (astrophysical calibration) as well as for the marginal likelihood for the data given a signal (integrated over the calibration uncertainty). I show that calibration uncertainty always reduces search sensitivity and the amount of information available about astrophysical signals. Additionally, calibration uncertainty will fundamentally limit the precision to which loud signals can be constrained, a crucial factor when considering the scientific potential of proposed third-generation interferometers. For example, I estimate that with $1\%$ uncertainty in the detector response's amplitude and phase, one will only be able to measure the leading-order tidal parameter ($\tilde\Lambda$) for a 1.4+1.4$\,M_\odot$ system to better than $\pm 1$ ($\sim 0.2\%$ relative uncertainty) for signals with signal-to-noise ratios $\gtrsim 10^4$. At this signal-to-noise ratio, calibration uncertainty increases $\sigma_{\tilde\Lambda}$ by a factor of $2$ compared to stationary Gaussian noise alone. Furthermore, 1\% calibration uncertainty limits the precision to always be $\sigma_{\tilde\Lambda} \gtrsim 0.5$. I also show how to best select the frequencies at which calibration should be precisely constrained in order to minimize the information lost about astrophysical parameters. It is not necessary to constrain the calibration errors to be small at all frequencies to perform precise astrophysical inference for individual signals.

Galaxies lose mass as a result of their luminosity or gaseous outflows. I calculate the resulting radial migration of stars outwards and show that it could potentially be measured with high resolution spectrographs on the next generation of large telescopes. Substantial accretion of matter in dense cosmic environments could trigger inward stellar migration that would be even more easily measurable.

Using spectra obtained with the VLT/FORS2 and Gemini-S/GMOS instruments we have investigated carbon, nitrogen and sodium abundances in a sample of red giants in the Small Magellanic Cloud cluster Kron 3. The metallicity and luminosity of this cluster are comparable to those of Galactic globular clusters although with a notably younger age of $\sim$ 6.5 Gyr. Specifically we have investigated the strengths of the CH ($\lambda$ 4300 A) and CN ($\lambda$ 3800, $\lambda$ 4215) molecular bands finding a bimodality of CN band-strengths and a CH/CN anti-correlation. Application of spectrum synthesis techniques reveals a large spread ($\sim$1.2 dex) in nitrogen abundance and a spread in [C/Fe] of $\sim$0.3 dex after applying corrections for evolutionary mixing. We have also estimated sodium abundances from the strengths of the Na D lines finding a range of $\sim$0.8 dex in [Na/Fe] that correlates positively with the N abundances. This is the first star-by-star spectroscopic demonstration of correlated Na, N abundance variations in an intermediate-age star cluster, adding to existing photometric and spectroscopic indications of the presence of multiple populations in such clusters with masses in excess of $\sim 10^5$ solar masses. Our results confirm that the mechanism(s) responsible for the multiple populations observed in globular clusters cannot be an early cosmological effect applying only in old clusters, and provide a key additional factor in the quest to understand the origin of the abundance anomalies.

The dynamical friction timescale of massive globular clusters (GCs) in the inner regions of cuspy dark haloes in dwarf spheroidal (dSph) galaxies can be much shorter than the Hubble time. This implies that a small fraction of the GCs is expected to be caught close to the centre of these galaxies. We compare the radial distribution of GCs predicted in simple Monte Carlo models with that of a sample of $38$ spectroscopically confirmed GCs plus 17 GC candidates, associated mainly to low-luminosity dSph galaxies. If dark matter haloes follow an NFW profile, the observed number of off-center GCs at projected distances less than one half the galaxy effective radius is significantly higher than models predict. This timing problem can be viewed as a fine-tuning of the starting GC distances. As a result of the short sinking timescale for GCs in the central regions, the radial distribution of GCs is expected to evolve significantly during the next 1-2 Gyr. However, dark matter haloes with cores of size comparable to the galaxy effective radii can lead to a slow orbital in-spiral of GCs in the central regions of these galaxies, providing a simple solution to the timing problem. We also examine any indication of mass segregation in the summed distribution of our sample of GCs.

The corona is an integral component of Active Galactic Nuclei (AGN) which produces the bulk of the X-ray emission above 1--2 keV. However, many of its physical properties and the mechanisms powering this emission remain a mystery. In particular, the temperature of the coronal plasma has been difficult to constrain for large samples of AGN, as constraints require high quality broadband X-ray spectral coverage extending above 10 keV in order to measure the high energy cutoff, which provides constraints on the combination of coronal optical depth and temperature. We present constraints on the coronal temperature for a large sample of Seyfert 1 AGN selected from the Swift/BAT survey using high quality hard X-ray data from the NuSTAR observatory combined with simultaneous soft X-ray data from Swift/XRT or XMM-Newton. When applying a physically-motivated, non-relativistic disk reflection model to the X-ray spectra, we find a mean coronal temperature kT $=$ 84$\pm$9 keV. We find no significant correlation between the coronal cutoff energy and accretion parameters such as the Eddington ratio and black hole mass. We also do not find a statistically significant correlation between the X-ray photon index, $\Gamma$, and Eddington ratio. This calls into question the use of such relations to infer properties of supermassive black hole systems.

N-body simulation is an invaluable tool to understand cosmic velocity field. However, simulation samples velocity only at particle locations that leads to information loss when projecting velocity fields onto structured grid. Here we extract two-point statistics without field projection. These statistics, i.e. correlation/structure/dispersion functions in real space and spectrum functions in Fourier space, are modeled on both small and large scales. Kinematic relations between statistical measures are fully developed for incompressible, constant divergence, and irrotational flow. The nature of flow can be identified by these relations. Much more complex than incompressible flow, peculiar velocity of dark matter flow is of constant divergence on small scale and irrotational on large scale. Incompressible and constant divergence flow share same kinematic relations for even order correlations. The limiting correlation of velocity $\rho_L=1/2$ on the smallest scale ($r=0$) is a unique feature of collisionless flow ($\rho_L=1$ for incompressible flow). Assuming gravity is the only interaction and no radiation produced, this feature leads to an increase in particle mass converted from kinetic energy upon "annihilation". On large scale, transverse velocity correlation has an exponential form with a comoving scale $r_2$=21.3Mpc/h that maybe related to the size of sound horizon. All other correlation/structure/dispersion/spectrum functions for velocity/density/potential are derived analytically from kinematic relations for irrotational flow. On small scale, longitudinal structure function follows one-fourth law of $S^l_2\propto r^{1/4}$. All other statistical measures are obtained from kinematic relations for constant divergence flow. Vorticity is negatively correlated for $r$ between 1 and 7Mpc/h. Divergence is negatively correlated for $r$>30Mpc/h that leads to negative density correlation.

High-contrast imaging studies of debris disks have revealed a significant diversity in their morphologies, including large-scale asymmetries. Theories involving stellar flybys, an external source of gravitational disturbance, have offered a plausible explanation for the origin of these morphological variations. Our study is an experiment to gain empirical evidence that has been lacking from such theories. We explore this paradigm by using astrometric and radial velocity measurements from the Gaia DR2 and ground-based observations to trace the trajectories of 625 stars in the Sco-Cen OB Association from 5 Myr in the past to 2 Myr in the future. We identified 119 stars that had at least one past flyby event occurring within one Hill radius, and 23 of these experienced flybys within 0.5 Hill radius. We found no evidence of a significant correlation between the presence of flyby events and infrared excess detections, although the sample is not uniformly sensitive to infrared excess emission. Ten stars that had past flyby events host resolved circumstellar disks that appear relatively symmetric in the existing data except for the circumbinary disk surrounding HD 106906. We determined the trajectory and relative velocity of each of these flyby events, and compared these to the geometry of the spatially-resolved disks. Future work is needed to measure the kinematics of lower mass stars and to improve sensitivity to circumstellar disks for the entire sample.

We present a spectral study of the single magnetically active K giant OP And in the period 1979 -- 2018, monitoring the variability of the activity indicator line H${\alpha}$. Original data obtained in the period 2015 -- 2018 with the echelle spectrograph ESpeRo at the 2m telescope of the National Astronomical Observatory Rozhen in Bulgaria, previously unpublished original data obtained in the period 1997 -- 2007 and on one night in 2013 with the Coude spectrograph at the same telescope, as well as data from the literature are presented in this study. The variability of the H${\alpha}$ line reveals that the activity level of OP And is higher in the period 1993 -- 2000, while during the period 2008 -- 2010 it is lower, possibly close to a minimum. Also, our data for the period 2015 -- 2018 show that the activity level is increasing again. Spectral observations of the activity indicators CaII H&K lines and CaII IR triplet are sparse during the studied period. We use such ones when possible to confirm the detection of some flare events. The structure of H${\alpha}$ changes with the activity level: when the activity is higher, we observe a blue-shifted component of this line, interpreted as an expanding area above the photosphere, but during a lower activity period it is almost absent. Our results are in a good agreement with the idea that the magnetic field controls the mass outflow in this giant. More years of observations are necessary to determine the eventual activity cycle of OP And.

The radio-loud quasar CTD 135 (2234+282, J2236+2828) has been proposed as a candidate compact symmetric object (CSO), based on its symmetric radio structure revealed by multi-frequency very long baseline interferometry (VLBI) imaging observations on milliarcsec angular scales. CSOs are known as young jetted active galactic nuclei (AGN) whose relativistic plasma jets are misaligned with respect to the line of sight. The peculiarity of CTD 135 as a CSO candidate was its detection in gamma-rays, while the vast majority of known gamma-ray emitting AGN are blazars with jets pointing close to our viewing direction. Since only a handful of CSOs are known as gamma-ray sources, the unambiguous identification of a single candidate is important for studying this rare class of objects. By collecting and interpreting observational data from the recent literature, we revisit the classification of CTD 135. We present evidence that the object, based on its flat-spectrum radio core with high brightness temperature, variability at multiple wavebands, and infrared colours should be classified as a blazar rather than a CSO

We study a parametrization of the deceleration parameter in a tilted universe, namely a cosmological model equipped with two families of observers. The first family follows the smooth Hubble flow, while the second are the real observers residing in a typical galaxy inside a bulk flow and moving relative to the smooth Hubble expansion with finite peculiar velocity. We use the compilation of Type Ia Supernovae (SnIa) data, as described in the Pantheon dataset, to find the quality of fit to the data and study the redshift evolution of the deceleration parameter. In so doing, we consider two alternative scenarios, assuming that the bulk-flow observers live in the $\Lambda$CDM and in the Einstein-de Sitter universe. We show that a tilted Einstein-de Sitter model can reproduce the recent acceleration history of the universe, without the need of a cosmological constant or dark energy, by simply taking into account linear effects of peculiar motions. By means of a Markov Chain Monte Carlo (MCMC) method, we also constrain the magnitude and the uncertainties of the parameters of the two models. From our statistical analysis, we find that the tilted Einstein-de Sitter model, equipped with one or two additional parameters that describe the assumed large-scale velocity flows, performs similar to the standard $\Lambda$CDM paradigm in the context of model selection criteria (Akaike Information Criterion and Bayesian Information Criterion).

An original method of estimating the projected rotational velocity, $vsini$, and the macroturbulent velocity, $v_{\rm mac}$, of evolved M giant stars is presented. It is based on the use of spectrum synthesis and multi-line analysis tools. The goal is to fit the mean line profile of observations with that of synthetic spectra. The method is applied to the red giant star RZ Ari and the results $v\sin i$ = 6.0 $\pm$ 0.5 km/s and $v_{\rm mac}$ = 2.0 $\pm$ 1.0 km/s are obtained.

We present an analysis aimed at combining cosmological constraints from number counts of galaxy clusters identified through the Sunyaev-Zeldovich effect, obtained with the South Pole Telescope (SPT), and from Lyman-$\alpha$ spectra obtained with the MIKE/HIRES and X-shooter spectrographs. The SPT cluster analysis relies on mass calibration based on weak lensing measurements, while the Lyman-$\alpha$ analysis is built over a suite of hydrodynamical simulations for the extraction of mock spectra. The resulting constraints exhibit a tension ($\sim 3.3\sigma$) between the low $\sigma_8$ values preferred by the low-redshift cluster data, $\sigma_8=0.74 ^{+0.03}_{-0.04}$, and the higher one preferred by the high-redshift Lyman-$\alpha$ data, $\sigma_8=0.91 ^{+0.03}_{-0.03}$. We present a detailed analysis in order to understand the origin of this tension and, in particular, to establish whether it arises from systematic uncertainties related to the assumptions underlying the analyses of cluster counts and/or Lyman-$\alpha$ forest. We found this tension to be robust with respect to the choice of modeling of the IGM, even when including possible systematics from unaccounted sub-Damped Lyman-$\alpha$ (DLA) and Lyman-limit systems (LLS) in the Lyman-$\alpha$ data. We conclude that to solve this tension from the SPT side would require a large bias on the cluster mass estimate, or from the Lyman-$\alpha$ side large unaccounted errors on the Lyman-$\alpha$ mean fluxes, respectively. Our results have important implications for future analyses based on cluster number counts from future large photometric surveys (e.g. Euclid and LSST) and on larger samples of high-redshift quasar spectra (e.g. DESI and WEAVE surveys). If confirmed at the much higher statistical significance reachable by such surveys, this tension could represent a significant challenge for the standard $\Lambda$CDM paradigm.

We investigated the physical and chemical properties of the gas and dust components in a carbon-rich planetary nebula (PN) IC2165 using two-dimensional (2-D) emission-line maps with superior resolution. The extinction map is generated in a self-consistent and assumption-free manner. The circumstellar gas-to-dust mass ratio (GDR) map ranges radially from 1210 in the central nebula filled with hot gas plasma to 120 near the ionisation front. The determined GDR is comparable to ~400, which is commonly adopted for carbon-rich asymptotic giant branch (AGB) stars, and ~100 for ISM. Except for the inner regions, the GDR in IC2165 is nearly the same as in such AGB stars, indicating that most dust grains withstand the harsh radiation field without being destroyed. The gas and dust mass distributions concentrated in the equatorial plane may be related to the nonisotropic mass loss during the AGB phase and nebula shaping. The spatial distributions of electron densities/temperatures and ionic/elemental abundances were investigated herein. We determined 13 elemental abundances using PSF-matched spatially integrated multiwavelength spectra extracted from the same aperture. Their values are consistent with values predicted by a theoretical model for stars of initially 1.75 Msun and Z = 0.003. Finally, we constructed the photoionisation model using our distance measurement to be consistent with all derived quantities, including the GDR and gas and dust masses and post-AGB evolution. Thus, we demonstrate the capability of Seimei/KOOLS-IFU and how the spatial variation of the gas and dust components in PNe derived from IFU observations can help understand the evolution of the circumstellar/interstellar medium.

This paper investigates the impact of radiative and mechanical feedback from O-type stars on their parent molecular clouds and the triggering of formation of future generation of stars. We study the infrared bubble S111 created by the embedded massive stellar cluster G316.80-0.05. A significant fraction of gas in shells created due to the compression of the ambient medium by expanding bubbles is photodissociated by the stellar radiation. The kinematics of the shells are thus best studied using spectroscopic observations of singly ionized carbon, the most dominant species. We have used the velocity-resolved maps of the $^2{\rm P}_{3/2}\rightarrow ^2{\rm P}_{1/2}$ transition of [C II] at 158 micron, the J=2-1 transition of 13CO and C18O, and the J=1-0 transition of HCO^+ to study the rim of the bubble S111 that partly coincides with the southern part of the infrared dark ridge G316.75. The [C II] spectra conclusively show evidence of a shell expanding with a moderate velocity of ~7 km/s, which amounts to a kinetic energy that is ~0.5-40 times the thermal energy of the H II region. The pressure causing the expansion of the H II region arises mainly from the hydrogen ionization and the dust-processed radiation. Among the far-infrared sources located in the compressed shells, we find the core G316.7799-0.0942 to show broad spectral features consistent with outflow activity and conclude that it is a site of active star formation. Based on the age of the H II region we conclude that this expanding H II region is responsible for the triggering of the current star formation activity in the region.

Natural Inflation with non-minimal coupling NMC to gravity $(\xi \phi^2)$ is investigated in the context of an extended gravity of the form $R+ \alpha R^2$. The treatment is performed in the Palatini formalism.We discuss various limitations of the model "$\alpha \gg 1$ and $\alpha \ll 1$" in light of two scenarios of inflation: a "Slow roll" and a "Constant roll".By analyzing the observational consequences of the model, our results show a significant improvement regarding compatibility between the theoretical results of this model and the experimental ones from Planck 2018, as exemplified by the tensor-to-scalar ratio and spectral index. Furthermore, a broader range of parameter space for natural inflation is now compatible with the confidence contours of Planck results. The notable effect of both NMC and $\alpha R^2$'s contributions makes it a significant improvement. $\alpha R^2$ gravity influences scalar-tensor ratio values. In contrast, NMC has a more significant impact on the spectral index values. Contributions to this area allow more excluded intervals to be included and are compatible with observational data.

Stellar rotation is a fundamental observable that drives different aspects of stellar and planetary evolution. In this work, we present an unprecedented manifold analysis of 160 B-type stars with light curves collected by the TESS space mission using three different procedures (Fast Fourier Transform, Lomb-Scargle, and wavelet techniques), accompanied by rigorous visual inspection in the search for rotation periodicities. This effort provides rotational periodicities for 6 new TESS B-type stars and confirmed periodicities for 22 targets with rotation periods previously listed in the literature. For other 61 stars, already classified as possible rotational variables, we identify noisy, pulsational, binarity, or ambiguous variability behavior rather than rotation signatures. The total sample of 28 potential rotators shows an overlap of different classes of rotational variables, composed of $\alpha^2$ Canum Venaticorum, rotating ellipsoidal and SX Arietis stars. The combination of the three techniques applied in our analysis offers a solid path to overcome the challenges in the discrimination of rotation from other variabilities in stellar light curves, such as pulsation, binarity or other effects that have no physical meaning. Finally, the rotational periodicities reported in the present study may represent important constraints for improving stellar evolution models with rotation, as well as asteroseismic studies of hot stars.

In this work, we have performed numerical simulations of primordial black hole (PBH) formation in the Friedman-Lema\^itre-Robertson-Walker universe filled by radiation fluid, introducing the local-type non-Gaussianity to the primordial curvature fluctuation. We have compared the numerical results from simulations with previous analytical estimations on the threshold value for PBH formation done in the previous paper arXiv:2109.00791, particularly for negative values of the non-linearity parameter $f_{\rm NL}$. Our numerical results show the existence of PBH formation of (the so-called) type I also in the case $f_{\rm NL} \lesssim -0.336$, which was not found in the previous analytical expectations using the critical averaged compaction function. In particular, although the universal value for the averaged critical compaction function $\bar{\mathcal{C}}_{c}=2/5$ found previously in the literature is not satisfied for all the profiles considered in this work, an alternative direct analytical estimate has been found to be roughly accurate to estimate the thresholds, which gives the value of the critical averaged density with a few $\%$ deviation from the numerical one for $f_{\rm NL}\gtrsim -1$.

IXPE is a Small Explorer mission that was launched at the end of 2021 to measure the polarization of X-ray emission from tens of astronomical sources. Its focal plane detectors are based on the Gas Pixel Detector, which measures the polarization by imaging photoelectron tracks in a gas mixture and reconstructing their initial directions. The quality of the single track, and then the capability of correctly determining the original direction of the photoelectron, depends on many factors, e.g., whether the photoelectron is emitted at low or high inclination with respect to the collection plane or the occurrence of a large Coulomb scattering close to the generation point. The reconstruction algorithm used by IXPE to obtain the photoelectron emission direction, also calculates several properties of the shape of the tracks which characterize the process. In this paper we compare several such properties and identify the best one to weight each track on the basis of the reconstruction accuracy. We demonstrate that significant improvement in sensitivity can be achieved with this approach and for this reason it will be the baseline for IXPE data analysis.

We describe the Gems of the Galaxy Zoos (Zoo Gems) project, a gap-filler project using short windows in the Hubble Space Telescope's schedule. As with previous snapshot programs, targets are taken from a pool based on position; we combine objects selected by volunteers in both the Galaxy Zoo and Radio Galaxy Zoo citizen-science projects. Zoo Gems uses exposures with the Advanced Camera for Surveys (ACS) to address a broad range of topics in galaxy morphology, interstellar-medium content, host galaxies of active galactic nuclei, and galaxy evolution. Science cases include studying galaxy interactions, backlit dust in galaxies, post-starburst systems, rings and peculiar spiral patterns, outliers from the usual color-morphology relation, Green Pea compact starburst systems, double radio sources with spiral host galaxies, and extended emission-line regions around active galactic nuclei. For many of these science categories, final selection of targets from a larger list used public input via a voting process. Highlights to date include the prevalence of tightly-wound spiral structure in blue, apparently early-type galaxies, a nearly complete Einstein ring from a group lens, redder components at lower surface brightness surrounding compact Green Pea starbursts, and high-probability examples of spiral galaxies hosting large double radio sources.

Any dark matter spikes surrounding black holes in our Galaxy are sites of significant dark matter annihilation, leading to a potentially detectable neutrino signal. In this paper we examine $10-10^5 M_\odot$ black holes associated with dark matter spikes that formed in early minihalos and still exist in our Milky Way Galaxy today, in light of neutrino data from the ANTARES and IceCube detectors. In various regions of the sky, we determine the minimum distance away from the solar system that a dark matter spike must be in order to have not been detected as a neutrino point source for a variety of representative dark matter annihilation channels. Given these constraints on the distribution of dark matter spikes in the Galaxy, we place significant limits on the formation of the first generation of stars in early minihalos -- stronger than previous limits from gamma-ray searches in Fermi Gamma-Ray Space Telescope data. The larger black holes considered in this paper may arise as the remnants of Dark Stars after the dark matter fuel is exhausted; thus neutrino observations may be used to constrain the properties of Dark Stars. The limits are particularly strong for heavier WIMPs. For WIMP masses $\sim 5 \,$TeV, we show that $\lesssim 10 \%$ of minihalos can host first stars that collapse into BHs larger than $10^3 M_\odot$.

Releasing the 12.5-year pulsar timing array data, the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has recently reported the evidence for a stochastic common-spectrum which would herald the detection of a stochastic gravitational wave background (SGWB) for the first time. We investigate if the signal could be generated from the end of a few MeV but still phenomenologically viable double-field inflation when the field configuration settles to its true vacuum. During the double-field inflation at such scales, bubbles of true vacuum that can collapse to LIGO mass and heavier primordial black holes form. We show that the produced gravitational wave spectrum matches the NANOGrav SGWB signal only when this process happens through a first-order phase transition. Using LATTICEEASY, we also examine the previous observation in the literature that by lowering the scale of preheating, despite the shift of the peak frequency of the gravitational wave profile to smaller values, the amplitude of the SGWB could be kept almost constant. We notice that this observation breaks down at the preheating scale, $M\lesssim 10^{-14}~m_{{}_{\rm Pl}}$.

We investigate at what abundances various hydrocarbon molecules (e.g. acetylene (C$_2$H$_2$), ethylene (C$_2$H$_4$), and methane (CH$_4$)) become detectable when observing the atmospheres of various planets using JWST. We focused on atmospheric models based on the parameters of a small sample of planets: HD 189733b, HD 209458b, HD 97658b, and Kepler-30c. We computed model transmission spectra using the Bayesian retrieval package ARCiS. We simulated observed spectra using the PandExo package. We subsequently ran retrievals on these spectra to determine whether the molecular abundances can be accurately retrieved from these simulated spectra. We find that generally we can detect and retrieve abundances of the hydrocarbon species as long as they have a VMR above approximately 1x10$^{-7}$-1x10$^{-6}$, at least for the brighter targets. There are variations based on planet type and instrument(s) used, and these limits will likely change depending on the abundance of other strong absorbers. We also find scenarios where the presence of one hydrocarbon is confused with another; this is often improved when two instruments are combined. C$_2$H$_2$, CH$_4$, and C$_2$H$_4$ will all be detectable with JWST, provided they are present in high enough abundances, and that the optimal instruments are chosen. A combination of two instruments, either NIRSpec G395M and MIRI LRS, or NIRCam F322W2 and MIRI LRS, is best for observing these species in bright exoplanet systems, with NIRSpec G395M and MIRI LRS the best option for the HD 189733b-like atmosphere with clouds included. The use of NIRSpec Prism is tentatively found to be best for fainter targets, potentially in combination with the MIRI LRS slit mode, although the target we test is too faint to draw any strong conclusions. Instrument sensitivity, noise, and wavelength range are all thought to play a role in being able to distinguish spectral features.

In the context of the recently proposed contact-binary model (Sekanina 2021), I investigate the circumstances of the first perihelion passage of the Kreutz sungrazers in orbits with barycentric periods near 735 yr, following the initial near-aphelion splitting of the presumed progenitor, Aristotle's comet of 372 BC. Given favorable conditions at this breakup and at episodes of secondary fragmentation in its aftermath, the fragments should have arrived at their first perihelion nearly simultaneously, reminiscent of the anticipated outcome for the two-superfragment model's perihelion return of AD 356 (Sekanina & Chodas 2004). The relevant case of a swarm of Kreutz sungrazers is examined to appraise possible scientific ramifications of the brief remark by Ammianus Marcellinus, a Roman historian, that "in broad daylight comets were seen" in late AD 363, only seven years later. The tested scenario, which does not contradict Ammianus' narrative and is consistent with the contact-binary model, involves a set of ten sungrazers visible in the daytime, all reaching perihelion over a period of 4.6 days. As part of this work, I comment on the role of the rapidly developing, brilliant post-perihelion tail; revise the apparent magnitude typical for the first and last naked-eye sightings; compare the visibility conditions in full daylight, in twilight, and at night; and, for the first time, present circumstantial evidence that favors comet X/1106 C1 as the parent to C/1843 D1 rather than to C/1882 R1 and C/1965 S1.

Contact binary star systems represent the long-lived penultimate phase of binary evolution. Population statistics of their physical parameters inform understanding of binary evolutionary pathways and end products. We use light curves and new optical spectroscopy to conduct a pilot study of ten (near-)contact systems in the long-period (P>0.5 d) tail of close binaries in the Kepler field. We couple PHOEBE light curve models with Markov-Chain Monte Carlo analyses to compute Bayesian probabilities on inclinations, fillout factors, mass ratios, third-light fractions, and component temperature ratios. Mass ratios and third-light contributions measured from spectra agree well with those inferred from the light curves. Most binaries in the pilot study have extreme mass ratios q<0.32. At least eight are probable triples. A Bayesian analysis of all 783 Kepler 0.15 d<P<2 d (near-)contact binaries forming an unbiased sample of unprecedented size and photometric precision results in 178 probable contact systems, 114 probable detached systems, and 491 ambiguous systems. We report best-fitting and 16th/50th/84th percentile parameters of modeled systems. Contact systems are rare at periods P>0.5 d, as are systems with mass ratios near unity (q>0.8). There exists an empirical mass ratio lower limit $q_{min}$(P)~0.05-0.15 below which contact systems are absent, in agreement with a new set of theoretical predictions we obtain by modeling the evolution of contact systems under the constraints of mass and angular momentum conservation. Pre-merger systems should be found at at long periods and near this mass ratio lower limit, which rises from q=0.044 for P=0.74 d to q=0.15 at P=2.0 d. These findings support a scenario whereby mass transfer drives systems toward extreme $q$ and larger $P$ until the onset of the Darwin instability at $q_{min}$ precipitates a merger.

Gamma-ray bursts are traditionally classified as short and long bursts based on their $T_{\rm 90}$ value (the time interval during which an instrument observes $5\%$ to $95\%$ of gamma-ray/hard X-ray fluence). However, $T_{\rm 90}$ is dependent on the detector sensitivity and the energy range in which the instrument operates. As a result, different instruments provide different values of $T_{\rm 90}$ for a burst. GRB 210217A is detected with different duration by {\it Swift} and {\it Fermi}. It is classified as a long/soft GRB by {\it Swift}-BAT with a $T_{\rm 90}$ value of 3.76 sec. On the other hand, the sub-threshold detection by {\it Fermi}-GBM classified GRB 210217A as a short/hard burst with a duration of 1.024 sec. We present the multi-wavelength analysis of GRB 210217A (lying in the overlapping regime of long and short GRBs) to identify its actual class using multi-wavelength data. We utilized the $T_{\rm 90}$-hardness ratio, $T_{\rm 90}$-\Ep, and $T_{\rm 90}$-$t_{\rm mvts}$ distributions of the GRBs to find the probability of GRB 210217A being a short GRB. Further, we estimated the photometric redshift of the burst by fitting the joint XRT/UVOT SED and place the burst in the Amati plane. We found that GRB 210217A is an ambiguous burst showing properties of both short and long class of GRBs.

The Sydney University Giant Air-shower Recorder (SUGAR) measured the muon component of extensive air showers with a unique array of muon detectors. The SUGAR data allows us to reconstruct the empirical dependence of muon density on the distance from the axis of the shower, the lateral distribution function (LDF). We compare the shape of this function with the predictions of hadronic-interaction models, QGSJET-II-04 and EPOS-LHC. We find a difference between the observed data and the simulation: the observed muon density falls faster with the increased core distance than it is predicted in simulations. This observation may be important for interpretation of the energy-dependent discrepancies in the simulated and observed numbers of muons in air showers, known as the "muon excess".

We analyze the properties that any late-time modification of the $\Lambda$CDM expansion history must have in order to consistently solve both the $H_0$ and the $\sigma_8$ tensions. Taking a model-independent approach, we obtain a set of necessary conditions that can be applied to generic late-time extensions. Our results are fully analytical and merely based on the assumptions that the deviations from the $\Lambda$CDM background remain small. For the concrete case of a dark energy fluid with equation of state $w(z)$, we derive the following general requirements: (i) Solving the $H_0$ tension demands $w(z)<-1$ at some $z$ (ii) Solving both the $H_0$ and $\sigma_8$ tensions requires $w(z)$ to cross the phantom divide. Finally, we also allow for small deviations on the effective gravitational constant. In this case, our method is still able to constrain the functional form of these deviations.

In this paper, we continue our study of the motion of spinning test bodies orbiting Kerr black holes. Non-spinning test bodies follow geodesics of the spacetime in which they move. A test body's spin couples to the curvature of that spacetime, introducing a ``spin-curvature force'' which pushes the body's worldline away from a geodesic trajectory. The spin-curvature force is an important example of a post-geodesic effect which must be modeled carefully in order to accurately characterize the motion of bodies orbiting black holes. One motivation for this work is to understand how to include such effects in models of gravitational waves produced from the inspiral of stellar mass bodies into massive black holes. In this paper's predecessor, we describe a technique for computing bound orbits of spinning bodies around black holes with a frequency-domain description which can be solved very precisely. In that paper, we present an overview of our methods, as well as present results for orbits which are eccentric and nearly equatorial (i.e., the orbit's motion is no more than $\mathcal{O}(S)$ out of the equatorial plane). In this paper, we apply this formulation to the fully generic case -- orbits which are inclined and eccentric, with the small body's spin arbitrarily oriented. We compute the trajectories which such orbits follow, and compute how the small body's spin affects important quantities such as the observable orbital frequencies $\Omega_r$, $\Omega_\theta$ and $\Omega_\phi$.

A photon sphere is known as the geometrical structure shaping a black hole shadow. The mechanism is well understood for static or stationary black hole spacetimes such as the Schwarzschild and the Kerr spacetimes. In this paper, we investigate and explicitly specify a photon sphere that shapes a black hole shadow in a dynamical spacetime while taking the global structure of the spacetime into account. We consider dynamical and eternal black hole cases of the Vaidya spacetime, which represents a spherically symmetric black hole with accreting null dust. First, we numerically show that there are the dynamical photon sphere and photon orbits corresponding to the shadow edge in a moderate accretion case. Second, the photon spheres are derived analytically in special cases. Finally, we discuss the relation between our photon sphere and the several notions defined as a photon sphere generalization.

It has recently been revealed that spinning black holes of the photon-fluid model can support acoustic `clouds', stationary density fluctuations whose spatially regular radial eigenfunctions are determined by the $(2+1)$-dimensional Klein-Gordon equation of an effective massive scalar field. Motivated by this intriguing observation, we use {\it analytical} techniques in order to prove a no-short hair theorem for the composed acoustic-black-hole-scalar-clouds configurations. In particular, it is proved that the effective lengths of the stationary bound-state co-rotating acoustic scalar clouds are bounded from below by the series of inequalities $r_{\text{hair}}>{{1+\sqrt{5}}\over{2}}\cdot r_{\text{H}}>r_{\text{null}}$, where $r_{\text{H}}$ and $r_{\text{null}}$ are respectively the horizon radius of the supporting black hole and the radius of the co-rotating null circular geodesic that characterizes the acoustic spinning black-hole spacetime.

Gravitational waves (GWs) generate oscillating electromagnetic effects in the vicinity of external electric and magnetic fields. We discuss this phenomenon with a particular focus on reinterpreting the results of axion haloscopes based on lumped-element detectors, which probe GWs in the 100 kHz-100 MHz range. Measurements from ABRACADABRA and SHAFT already place bounds on GWs, although the present strain sensitivity is weak. However, we demonstrate that the sensitivity scaling with the volume of such instruments is significant - faster than for axions - and so rapid progress will be made in the future. With no modifications, DMRadio-m$^3$ will have a GW strain sensitivity of $h \sim 10^{-20}$ at 200 MHz. A simple modification of the pickup loop used to readout the induced magnetic flux can parametrically enhance the GW sensitivity, particularly at lower frequencies.

The origin of water on the Earth is a long-standing mystery, requiring a comprehensive search for hydrous compounds, stable at conditions of the deep Earth and made of Earth-abundant elements. Previous studies usually focused on the current range of pressure-temperature conditions in the Earth's mantle and ignored a possible difference in the past, such as the stage of the core-mantle separation. Here, using ab initio evolutionary structure prediction, we find that only two magnesium hydrosilicate phases are stable at megabar pressures, $\alpha$-Mg$_2$SiO$_5$H$_2$ and $\beta$-Mg$_2$SiO$_5$H$_2$, stable at 262-338 GPa and >338 GPa,respectively (all these pressures now lie within the Earth's iron core). Both are superionic conductors with quasi-one-dimensional proton diffusion at relevant conditions. In the first 30 million years of Earth's history, before the Earth's core was formed, these must have existed in the Earth, hosting much of Earth's water. As dense iron alloys segregated to form the Earth's core, Mg$_2$SiO$_5$H$_2$ phases decomposed and released water. Thus, now-extinct Mg$_2$SiO$_5$H$_2$ phases have likely contributed in a major way to the evolution of our planet.

The NASA Astrophysics Data System (ADS), a critical research service for the astrophysics community, strives to provide the most accessible and inclusive environment for the discovery and exploration of the astronomical literature. Part of this goal involves creating a digital platform that can accommodate everybody, including those with disabilities that would benefit from alternative ways to present the information provided by the website. NASA ADS follows the official Web Content Accessibility Guidelines (WCAG) standard for ensuring accessibility of all its applications, striving to exceed this standard where possible. Through the use of both internal audits and external expert review based on these guidelines, we have identified many areas for improving accessibility in our current web application, and have implemented a number of updates to the UI as a result of this. We present an overview of some current web accessibility trends, discuss our experience incorporating these trends in our web application, and discuss the lessons learned and recommendations for future projects.

We discuss various methods for acquiring optical links in space using a dedicated acquisition sensor. Statistical models are developed and simple analytical equations derived that compare the performance between a single and dual spiral scan approach as well as between sequential and parallel acquisition of link chains. Simple derived analytical equations allow relating essential search parameters such as track width, variance of the uncertainty distribution, capture radius and scan speed to the probabilities of acquiring the links within a specific time. We also assess the probability of failing to acquire a link due to beam jitter and derive a simple analytical model that allows determining the maximum tolerable jitter for a given beam overlap and required probability of success. All results are validated by Monte Carlo simulations and applied to the concrete example of the GRACE FO mission.

We discuss gravitational physics in the Jordan and Einstein frames of $f (R)$ gravity coupled to the Standard Model. We elucidate the way in which the observed gravitational coupling arises in the Einstein frame for generic $f (R)$. We point out that the effect of "running units" in the Einstein frame is related to the fact that the explicit and implicit quantum parameters of the Standard Model, such as the Higgs vacuum expectation value and the parameter $\Lambda_\text{QCD}$, are modified by the conformal transformation of the metric and matter fields and become scalaron-dependent. Considering the scalaron of $f (R)$ gravity describing dark matter, we show that the effect of running units in this case is extremely weak, making two frames practically equivalent.

Gravitational waves from binary black hole and neutron star mergers are being regularly detected. As of 2021, ninety confident gravitational wave detections have been made by the LIGO and Virgo detectors. Work is ongoing to further increase the sensitivity of the detectors for the fourth observing run, including installing some of the A+ upgrades designed to lower the fundamental noise that limits the sensitivity to gravitational waves. In this review, we will overview how the LIGO detectors work, including their optical configuration and lock acquisition procedure, discuss the detectors' fundamental and technical noise limits and compare to the current measured sensitivity, and review the A+ upgrades currently being installed in the detectors.

According to the standard model of cosmology, LambdaCDM, the mass-energy budget of the current stage of the universe is not dominated by the luminous matter that we are familiar with, but instead by some form of dark matter (and dark energy). It is thus tempting to adopt scientific realism about dark matter. However, there are barely any constraints on the myriad of possible properties of this entity -- it is not even certain that it is a form of matter. In light of this underdetermination I advocate caution: we should not (yet) be dark matter realists. The "not(-yet)-realism" that I have in mind is different from Hacking's (1989) anti-realism, in that it is semantic rather than epistemological. It also differs from the semantic anti-realism of logical empiricism, in that it is naturalistic, such that it may only be temporary and does not automatically apply to all other unobservables (or even just to all other astronomical unobservables, as with Hacking's anti-realism). The argument is illustrated with the analogy of the much longer history of the concept of a gene, as the current state of the concept of dark matter resembles in some relevant ways that of the early concept of genes.

Time-delay interferometry (TDI) is a data processing technique for LISA designed to suppress the otherwise overwhelming laser noise by several orders of magnitude. It is widely believed that TDI can only be applied once all phase or frequency measurements from each spacecraft have been synchronized to a common time frame. We demonstrate analytically, using as an example the commonly-used Michelson combination X, that TDI can be computed using the raw, unsynchronized data, thereby avoiding the need for the initial synchronization processing step and significantly simplifying the initial noise reduction pipeline for LISA. Furthermore, the raw data is free of any potential artifacts introduced by clock synchronization and reference frame transformation algorithms, which allows to operate directly on the MHz beatnotes. As a consequence, in-band clock noise is directly suppressed as part of TDI, in contrast to the approach previously proposed in the literature, in which large trends in the beatnotes are removed before the main laser-noise reduction step and clock noise is suppressed in an extra processing step. We validate our algorithm with fullscale numerical simulations that use LISA Instrument and PyTDI and show that we reach the same performance levels as the previously proposed methods, ultimately limited by the clock sideband stability.

We derive exact general solutions (as opposed to attractor particular solutions) and corresponding first integrals for the evolution of a scalar field $\phi$ in a universe dominated by a background fluid with equation of state parameter $w_B$. In addition to the previously-examined linear [$V(\phi) = V_0 \phi$] and quadratic [$V(\phi) = V_0 \phi^2$] potentials, we show that exact solutions exist for the power law potential $V(\phi) = V_0 \phi^n$ with $n = 4(1+w_B)/(1-w_B) + 2$ and $n = 2(1+w_B)/(1-w_B)$. These correspond to the potentials $V(\phi) = V_0 \phi^6$ and $V(\phi) = V_0 \phi^2$ for matter domination and $V(\phi) = V_0 \phi^{10}$ and $V(\phi) = V_0 \phi^4$ for radiation domination. The $\phi^6$ and $\phi^{10}$ potentials can yield either oscillatory or non-oscillatory evolution, and we use the first integrals to determine how the initial conditions map onto each form of evolution. The exponential potential yields an exact solution for a stiff/kination ($w_B = 1$) background. We use this exact solution to derive an analytic expression for the evolution of the equation of state parameter, $w_\phi$, for this case.