Papers for Wednesday, May 03 2023

We investigate the ability of simple diagnostics based on Lagrangian descriptor (LD) computations of initially nearby orbits to detect chaos in conservative dynamical systems with phase space dimensionality higher than two. In particular, we consider the recently introduced methods of the difference ($D_L^n$) and the ratio ($R_L^n$) of the LDs of neighboring orbits, as well as a quantity ($S_L^n$) related to the finite-difference second spatial derivative of the LDs, and use them to determine the chaotic or regular nature of ensembles of orbits of a prototypical area-preserving map model, the 4-dimensional (4D) symplectic standard map. Using the distributions of the indices' values we determine appropriate thresholds to discriminate between regular and chaotic orbits, and compare the obtained characterization against that achieved by the Smaller Alignment Index (SALI) method of chaos detection, by recording the percentage agreement $P_A$ between the two classifications. We study the influence of various factors on the performance of these indices, and show that the increase of the final number of orbit iterations T and the order n of the indices (i.e. the dimensionality of the space where the considered nearby orbits lie), as well as the decrease of the distance $\sigma$ of neighboring orbits, increase the $P_A$ values along with the required computational effort. Balancing between these two factors we find appropriate T, n and $\sigma$ values, which allow the efficient use of the $D_L^n$, $R_L^n$ and $S_L^n$ indices as short time and computationally cheap chaos diagnostics achieving $P_A \gtrsim 90 \%$, with $D_L^n$ and $S_L^n$ having larger $P_A$ values than $R_L^n$. Our results show that the three LDs-based indices perform better for systems with large percentages of chaotic orbits.

Developments in the stability of modern spectrographs have led to extremely precise instrumental radial velocity (RV) measurements. For most stars, the detection limit of planetary companions with these instruments is expected to be dominated by astrophysical noise sources such as starspots. Correlated signals caused by rotationally-modulated starspots can obscure or mimic the Doppler shifts induced by even the closest, most massive planets. This is especially true for young, magnetically active stars where stellar activity can cause fluctuation amplitudes of $\gtrsim$0.1 mag in brightness and $\gtrsim$100 m s$^{-1}$ in RV semi-amplitudes. Techniques that can mitigate these effects and increase our sensitivity to young planets are critical to improving our understanding of the evolution of planetary systems. Gaussian processes (GPs) have been successfully employed to model and constrain activity signals in individual cases. However, a principled approach of this technique, specifically for the joint modeling of photometry and RVs, has not yet been developed. In this work, we present a GP framework to simultaneously model stellar activity signals in photometry and RVs that can be used to investigate the relationship between both time series. Our method, inspired by the $\textit{FF}^\prime$ framework of (Aigrain et al. 2012), models spot-driven activity signals as the linear combinations of two independent latent GPs and their time derivatives. We also simulate time series affected by starspots by extending the $\texttt{starry}$ software (Luger et al. 2019) to incorporate time evolution of stellar features. Using these synthetic datasets, we show that our method can predict spot-driven RV variations with greater accuracy than other GP approaches.

We report a detailed analysis of the XMM-Newton spectra of six Narrow-Line Seyfert 1 (NLS1) galaxies at redshift z = 0.35--0.92. Compared with the NLS1s at lower redshift in the previously most-studied sample, these NLS1s have larger black hole (BH) masses ($\log\,M_\text{BH}>7.5$) with similar or even lower Eddington ratios. Our extended XMM-Newton sample of NLS1s shows strong soft X-ray excess emission below 2 keV. The quantified soft excess strength does not show an obvious discrepancy from previous studies of the lower-redshift NLS1s. The systematic effect in the measurement of the Eddington ratio mainly lies in the bolometric correction factor. We also tentatively fit the spectra assuming two more physical models for the soft excess: warm Comptonization and relativistic reflection from the inner accretion disk. In the first scenario, we confirm the ubiquity of a warm and optically thick corona. The behavior of a single source can be better explained by relativistic reflection, although we cannot distinguish which model is a more favorable explanation for the soft excess based on the best-fit statistics.

In the coming years, Sunyaev-Zel'dovich (SZ) measurements can dramatically improve our understanding of the Intergalactic Medium (IGM) and the role of feedback processes on galaxy formation, allowing us to calibrate important astrophysical systematics in cosmological constraints from weak lensing galaxy clustering surveys. However, the signal is only measured in a two-dimensional projection, and its correct interpretation relies on understanding the connection between observable quantities and the underlying intrinsic properties of the gas, in addition to the relation between the gas and the underlying matter distribution. One way to address these challenges is through the use of hydrodynamical simulations such as the high-resolution, large-volume MillenniumTNG suite. We find that measurements of the optical depth, $\tau$, and the Compton-y parameter, $Y$, receive large line-of-sight contributions which can be removed effectively by applying a Compensated Aperture Photometry (CAP) filter. In contrast with other $\tau$ probes (e.g., X-rays and Fast Radio Bursts), the kSZ-inferred $\tau$ receives most of its signal from a confined cylindrical region around the halo due to the velocity decorrelation along the line-of-sight. Additionally, we perform fits to the $Y-M$ and $\tau-M$ scaling relations and report best-fit parameters adopting the smoothly broken power law (SBPL) formalism. We note that subgrid physics modeling can broaden the error bar on these by 30\% for intermediate-mass halos ($\sim$$10^{13} \, {\rm M}_{\odot}$). The scatter of the scaling relations can be captured by an intrinsic dependence on concentration, and an extrinsic dependence on tidal shear. Finally, we comment on the effect of using galaxies rather than halos in real observations, which can bias the inferred SZ profiles by $\sim$20\% for $L_\ast$-galaxies.

Dark matter subhaloes are key for the predictions of simulations of structure formation, but their existence frequently ends prematurely due to two technical issues, namely numerical disruption in N-body simulations and halo finders failing to identify them. Here we focus on the second issue, using the phase-space friends-of-friends halo finder ROCKSTAR as a benchmark (though we expect our results to translate to comparable codes). We confirm that the most prominent cause for losing track of subhaloes is tidal distortion rather than a low number of particles. As a solution, we present a flexible post-processing algorithm that tracks all subhalo particles over time, computes subhalo positions and masses based on those particles, and progressively removes stripped matter. If a subhalo is lost by the halo finder, this algorithm keeps tracking its so-called ghost until it has almost no particles left or has truly merged with its host. We apply this technique to a large suite of N-body simulations and restore lost subhaloes to the halo catalogues, which has a dramatic effect on key summary statistics of large-scale structure. Specifically, the subhalo mass function increases by about 50% and the halo correlation function increases by a factor of two at small scales. While these quantitative results are somewhat specific to our algorithm, they demonstrate that particle tracking is a promising way to reliably follow haloes and reduce the need for orphan models. Our algorithm and augmented halo catalogues are publicly available.

Proximity zones of high-redshift quasars are unique probes of their central supermassive black holes as well as the intergalactic medium in the last stages of reionization. We present 22 new measurements of proximity zones of quasars with redshifts between 5.8 and 6.6, using the enlarged XQR-30 sample of high-resolution, high-SNR quasar spectra. The quasars in our sample have UV magnitudes of $M_{1450}\sim -27$ and black hole masses of $10^9$$\unicode{x2013}$$10^{10}$ M$_\odot$. Our inferred proximity zone sizes are 2$\unicode{x2013}$7 physical Mpc, with a typical uncertainty of less than 0.5 physical Mpc, which, for the first time, also includes uncertainty in the quasar continuum. We find that the correlation between proximity zone sizes and the quasar redshift, luminosity, or black hole mass, indicates a large diversity of quasar lifetimes. Two of our proximity zone sizes are exceptionally small. The spectrum of one of these quasars, with $z=6.02$, displays, unusually for this redshift, damping wing absorption without any detectable metal lines, which could potentially originate from the IGM. The other quasar has a high-ionization absorber $\sim$0.5 pMpc from the edge of the proximity zone. This work increases the number of proximity zone measurements available in the last stages of cosmic reionization to 87. This data will lead to better constraints on quasar lifetimes and obscuration fractions at high redshift, which in turn will help probe the seed mass and formation redshift of supermassive black holes.

JWST has opened up a new observational probe of our Universe. The early data release by JWST have revealed several high redshift massive galaxy candidates by photometry, and some of them have been confirmed spectroscopically. We use these observations to study their implications on the primordial power spectrum. In the first part of this work, we use the data from the CEERS photometric survey, along with respective spectroscopic updates, to compute the cumulative comoving stellar mass density. We find that a very high star formation efficiency (unlikely in various theoretical scenarios) is required to explain these observations within $\Lambda$CDM cosmology. We show that the tension can be eased if the primordial power spectrum has a blue tilt on small length scales. The required blue tilt depends on the currently unknown star formation efficiency in these galaxy candidates. In the second part of this work, we study the spectroscopically confirmed galaxies reported in the JADES survey at redshift $z \gtrsim 10$, which have been shown to be consistent with $\Lambda$CDM cosmology. We investigate the implications of these measurements on a red-tilted primordial power spectrum. For these galaxies, we estimate the star formation efficiency from an earlier observation of galaxies (with similar redshifts) by the Spitzer telescope. We find that the star formation efficiency is an order of magnitude smaller than that required to explain the CEERS photometric observations mentioned earlier. Using the estimated star formation efficiency, we find the strongest constraints on the red tilt of the power spectrum on certain length scales. Our study shows that JWST observations will be an excellent probe of the power spectrum and can lead to novel discoveries.

From the thousands of known exoplanets, those that transit bright host stars provide the greatest accessibility toward detailed system characterization. The first known such planets were generally discovered using the radial velocity technique, then later found to transit. HD 17156b is particularly notable among these initial discoveries because it diverged from the typical hot Jupiter population, occupying a 21.2 day eccentric ($e = 0.68$) orbit, offering preliminary insights into the evolution of planets in extreme orbits. Here we present new data for this system, including ground and space-based photometry, radial velocities, and speckle imaging, that further constrain the system properties and stellar/planetary multiplicity. These data include photometry from the Transiting Exoplanet Survey Satellite (TESS) that cover five transits of the known planet. We show that the system does not harbor any additional giant planets interior to 10 AU. The lack of stellar companions and the age of the system indicate that the eccentricity of the known planet may have resulted from a previous planet-planet scattering event. We provide the results from dynamical simulations that suggest possible properties of an additional planet that culminated in ejection from the system, leaving a legacy of the observed high eccentricity for HD 17156b.

We analytically solve the problem of Bose star growth in the bath of gravitationally interacting particles. We find that after nucleation of this object, the bath is described by a self-similar solution of the kinetic equation, which is an attractor. Together with the conservation laws, this fixes mass evolution of the Bose star. Our results explain slowdown of the star growth at a certain "core-halo" mass, but also predict formation of the heavier and lighter objects in magistral dark matter models.

A long-standing observed curiosity of globular clusters (GCs) has been that both the number and total mass of GCs in a galaxy are linearly correlated with the galaxy's virial mass, whereas its stellar component shows no such linear correlation. This work expands on an empirical model for the numbers and ages of GCs in galaxies presented by Valenzuela et al. (2021) that is consistent with recent observational data from massive elliptical galaxies down to the dwarf galaxy regime. Applying the model to simulations, GC numbers are shown to be excellent tracers for the dark matter (DM) virial mass, even when distinct formation mechanisms are employed for blue and red GCs. Furthermore, the amount of DM smooth accretion is encoded in the GC abundances, therefore providing a measure for an otherwise nearly untraceable component of the formation history of galaxies.

Solar and heliospheric parameters can depict notable differences between the northern and southern hemisphere. Although hemispheric asymmetries of some heliospheric parameters vary systematically with Hale cycle, this is not common for solar parameters. Also, no physical mechanism exists which can explain systematic hemispheric asymmetries. We use a novel method of high heliolatitudes to increase the fraction of one hemisphere in solar 10.7cm radio fluxes and sunspot numbers. We calculate sets of hemispheric high-latitude radio fluxes and sunspot numbers with increasing heliographic latitude during the last 75 years. We also normalise these fluxes by yearly means in order to study their continuous variation. We find that cycle maximum radio fluxes and sunspot numbers in each odd cycle (19, 21, 23) are larger at northern high latitudes, while in all even cycles (18, 20, 22 24) they are larger at southern latitudes. This alternation indicates a new form of Hale cycle variation in solar activity. Hemispheric differences at cycle maxima are 15% for radio flux and 23% for sunspot numbers. The difference is largest during cycle 19 and smallest in cycle 24. Continuous fluxes depict a Hale cycle in both hemispheres, with an opposite phase and amplitude of 5% in north and 4% in south. Hemispheric Hale cycle can be explained if there is a northward directed relic magnetic field, which is shifted northward. In odd cycles, the northern hemisphere is enhanced more than the southern hemisphere and, in even cycles, the northern hemisphere is reduced more than the southern hemisphere. The decrease of asymmetry during the 7 cycles can be explained if the relic shift oscillates at the 210-year Suess/deVries period. Gleissberg cycle consists of one off-equator excursion of the relic. Relic field in the Sun also offers a possibility for century-scale forecasting of solar activity.

Estimates suggest that while FRII jets appear to have lifetimes constrained to hundreds of millions of years, radio galaxies with FRI jets appear to be longer lived. We illustrate the nature of this time constraint from model perspectives, showing how compatibility between theory and data match in a way suggesting a key difference between active galaxies whose engines are characterized by accretion onto co-rotating versus counter-rotating black holes. We calculate a range of timescales for counter-rotating black holes for a range of accretion rates compatible with theory which we then compare to data. The validity of these timescales constitutes the most powerful recent piece of evidence for considering counter-rotation between black holes and accretion disks in high energy astrophysics.

The Cosmology Large Angular Scale Surveyor (CLASS) is a telescope array that observes the cosmic microwave background over 75\% of the sky from the Atacama Desert, Chile, at frequency bands centered near 40, 90, 150, and 220~GHz. This paper describes the CLASS data pipeline and maps for 40~GHz observations conducted from August 2016 to May 2022. We demonstrate how well the CLASS survey strategy, with rapid ($\sim10\,\mathrm{Hz}$) front-end modulation, recovers the large-scale Galactic polarization signal from the ground: the mapping transfer function recovers $\sim75$\% of $EE$, $BB$, and $VV$ power at $\ell=20$ and $\sim45$\% at $\ell=10$. We present linear and circular polarization maps over 75\% of the sky. Simulations based on the data imply the maps have a white noise level of $110\,\mathrm{\mu K\, arcmin}$ and correlated noise component rising at low-$\ell$ as $\ell^{-2.2}$. The transfer-function-corrected low-$\ell$ component is comparable to the white noise at the angular knee frequencies of $\ell\approx16$ (linear polarization) and $\ell\approx12$ (circular polarization). Finally, we present simulations of the level at which expected sources of systematic error bias the measurements, finding sub-percent bias for the $\Lambda\mathrm{CDM}$ $EE$ power spectra. Bias from $E$-to-$B$ leakage due to the data reduction pipeline and polarization angle uncertainty approaches the expected level for an $r=0.01$ $BB$ power spectrum. Improvements to the instrument calibration and the data pipeline will decrease this bias.

We present Australia Compact Telescope Array (ATCA) 21-cm observations of the nearby low-excitation radio galaxy (LERG) NGC 3100. This is the brightest galaxy of a loose group and hosts a young ($\sim 2$ Myr) radio source. The ATCA observations reveal for the first time the presence of neutral hydrogen (HI) gas in absorption in the centre of this radio galaxy, and in emission in two low-mass galaxies of the group and in a diffuse dark cloud in the proximity of NGC 3100. The sensitivity to low-column density gas ($N_{\rm HI}\sim 10^{19}$ cm$^{-2}$) allows us to reveal asymmetries in the periphery of most the HI-detected galaxies, suggesting that tidal interactions may be on-going. The diffuse cloud does not show a stellar counterpart down to $27$ mag/arcsec$^2$ and could be the remnant of these interactions. The analysis of the HI absorption complex in NGC 3100 indicates that the atomic phase of the hydrogen is distributed as its molecular phase (observed at arcsecond resolution through several carbon monoxide emission lines). We suggest that the interactions occurring within the group are causing turbulent cold gas clouds in the intra-group medium to be slowly accreted towards the centre of NGC 3100. This caused the recent formation of the cold circum-nuclear disk which is likely sustaining the young nuclear activity.

The most extensive survey of carbon monoxide (CO) gas in the Taurus molecular cloud relied on $^{12}$CO and $^{13}$CO $J=1 \rightarrow 0$ emission only, distinguishing the region where $^{12}$CO is detected without $^{13}$CO (named mask 1 region) from the one where both are detected (mask 2 region). We have taken advantage of recent $^{12}$CO $J=3\rightarrow2$ JCMT observations where they include mask 1 regions to estimate density, temperature, and $N$(CO) with a LVG model. This represents 1395 pixels out of $\sim$1.2 million in the mark 1 region. Compared to Pineda et al. (2010) results, and assuming a $T_\textrm{kin}$ of 30 K, we find a higher volume density of molecular hydrogen of 3.3$\rm \times\ 10^3$ $\textrm{cm}^{-3}$, compared to their 250--700 $\textrm{cm}^{-3}$ and a CO column density of 5.7$\rm \times\ 10^{15}\ \textrm{cm}^{-2}$, about a quarter of their value. The differences are important and show the necessity to observe several CO transitions to better describe the intermediate region between the dense cloud and the diffuse atomic medium. Future observations to extend the $^{12}$CO $J=3\rightarrow2$ mapping further away from the $^{13}$CO--detected region comprising mask 1 are needed to revisit our understanding of the diffuse portions of dark clouds.

With the advent of large-scale astronomical surveys such as the Zwicky Transient Facility (ZTF), the number of alerts generated by transient, variable and moving astronomical objects is growing rapidly, reaching millions per night. Concerning solar system minor planets, their identification requires linking the alerts of many observations over a potentially large time, leading to a very large combinatorial number. This work aims to identify new candidates for solar system objects from massive alert data streams produced by large-scale surveys, such as the ZTF and the Vera C. Rubin Observatory's Legacy Survey of Space and Time. Our analysis used the Fink alert broker capabilities to reduce the 111,275,131 processed alerts from ZTF between November 2019 and December 2022 to only 389,530 new solar system alert candidates over the same period. We then implemented a linking algorithm, Fink-FAT, to create real-time trajectory candidates from alert data and extract orbital parameters. The analysis was validated on ZTF alert packets linked to confirmed solar system objects from the Minor Planet Center database. Finally, the results were confronted against follow-up observations. Between November 2019 and December 2022, Fink-FAT extracted 327 new orbits from solar system object candidates at the time of the observations, over which 65 were still unreported in the MPC database as of March 2023. After two late follow-up observation campaigns of six orbit candidates, four were associated with known solar system minor planets, and two remain unknown. Fink-FAT is deployed in the Fink broker and successfully analyzes in real time the alert data from the ZTF survey by regularly extracting new candidates for solar system objects. Our scalability tests also show that Fink-FAT can handle the even larger volume of alert data that the Rubin Observatory will send.

Using Low Brightness Surface Galaxies (LBSG) rotational curves we inferred the free parameters of $\ell$-boson stars as a dark matter component. The $\ell$-boson stars are numerical solutions to the non-relativistic limit of the Einstein-Klein-Gordon system, the Schr\"odinger-Poisson (SP) system. These solutions are parametrized by an angular momentum number $\ell = (N-1)/2$ and an excitation number $n$. We perform a bayesian analysis by modifying the SimpleMC code to perform the parameter inference, for the cases with $\ell = 0$, $\ell = 1$ and multistates of $\ell$-boson stars. We used the Akaike information criterion (AIC), Bayesian information criterion and the Bayes factor to compare the excited state ($\ell$=1) and the multistate case with the ground state ($\ell$=0) as the base model due to its simplicity. We found that the data in most galaxies in the sample favours the multistates case and that the scalar field mass tends to be slightly bigger than the ground state case.

We present, for the first time, an observational test of the consistency relation for the large-scale structure (LSS) of the Universe through a joint analysis of the anisotropic two- and three-point correlation functions (2PCF and 3PCF) of galaxies. We parameterise the breakdown of the LSS consistency relation in the squeezed limit by $E_{\rm s}$, which represents the ratio of the coefficients of the shift terms in the second-order density and velocity fluctuations. $E_{\rm s}\neq1$ is a sufficient condition under which the LSS consistency relation is violated. A novel aspect of this work is that we constrain $E_{\rm s}$ by obtaining information about the nonlinear velocity field from the quadrupole component of the 3PCF without taking the squeezed limit. Using the galaxy catalogues in the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12, we obtain $E_{\rm s} = -0.92_{-3.26}^{+3.13}$, indicating that there is no violation of the LSS consistency relation in our analysis within the statistical errors. Our parameterisation is general enough that our constraint can be applied to a wide range of theories, such as multicomponent fluids, modified gravity theories, and their associated galaxy bias effects. Our analysis opens a new observational window to test the fundamental physics using the anisotropic higher-order correlation functions of galaxy clustering.

The matter power spectrum has been strongly constrained by astronomical measurements at large scales, but only weakly at small scales. Compared with the standard scenario, the deviation of the matter power spectrum at small scales has influence on the cosmological structure formation, e.g., the comoving number density of dark matter halos. The thermal history of the intergalactic medium (IGM) can be changed if dark matter is made of weakly interacting massive particles and can annihilate into standard model particles. The changes of the evolution of IGM could leave imprints on the relevant astronomical observations. Taking into account the dark matter annihilation, we investigate the impact of the deviation of matter power spectrum at small scales on the global 21-cm signal. In view of the measurements of the global 21-cm signal by the EDGES experiment, we explore the allowed parameter space of $m_s$, which describes the degree of deviation, by requiring the differential brightness temperature of the global 21-cm signal $\delta T_{21} \le -50~\rm mK$ at redshift $z=17$.

The James Webb Space Telescope (JWST) near-infrared camera (NIRCam) has been found to exhibit serious wisp-like structures in four of its eight short-wavelength detectors. The exact structure and strength of these wisps is highly variable with the position and orientation of JWST, so the use of static templates is non-optimal. Here we investigate a dynamic strategy to mitigate these wisps using long-wavelength reference images. Based on a suite of experiments where we embed a worst-case scenario median-stacked wisp into wisp-free images, we define suitable parameters for our wisp removal strategy. Using this setup we re-process wisp-affected public Prime Extragalactic Areas for Reionization and Lensing Science (PEARLS) data in the North Ecliptic Pole Time Domain Field (NEP-TDF) field, resulting in significant visual improvement in our detector frames and reduced noise in the final stacked images.

I observed the very fast nova V838 Her (Nova Her 1991, optical peak at 5-5.4 mag) during the fading phase of the nova eruption in 1991. I detected eclipses, for the first time in the world in any nova during eruption, and the epochs of the eclipses were reported to IAU Circular No. 5262. Although these epochs have been referenced in the literature, the light curves of these eclipses remained unpublished. Here, I present these light curves. The phase-averaged light curve around 1991 April 21 (mean V=13.1, 27 d after the optical peak) showed an 0.14 mag primary eclipse and an 0.03 mag secondary eclipse. Combined with the subsequent literature, the eclipses likely appeared after 1991 April 14 (V=12.5). It had been suggested that the accretion disk had already been re-established before this epoch and I found no strong argument against this. The early appearance of the secondary minimum appears to be a phenomenon common to very fast novae and it looks likely to be explained, at least partly, by a strongly heated secondary. This observation reinforces the possible interpretation of the early presence of a transient luminous donor for the fastest nova V1674 Her (Nova Her 2021). As a comparison and my motivation of the observation of V838 Her, I briefly review the early history of V1500 Cyg (Nova Cyg 1975).

We present the results from our deep optical photometric observations of Bochum 2 (Boc2) star cluster obtained using the $1.3$m Devasthal Fast Optical Telescope along with archival photometric data from Pan-STARRS2/2MASS/UKIDSS surveys. We also used high-quality parallax and proper motion data from the $Gaia$ Data Release 3. We found that the Boc2 cluster has a small size ($\sim$1.1 pc) and circular morphology. Using $Gaia$ parallax of member stars and isochrone fitting method, the distance of this cluster is estimated as $3.8\pm0.4$ kpc. We have found that this cluster holds young ($\sim5$ Myr) and massive (O$7-$O$9$) stars as well as an older population of low mass stars. We found that the massive stars have formed in the inner region of the Boc2 cluster in a recent epoch of star formation. We have derived mass function slope ($\Gamma$) in the cluster region as $-2.42\pm0.13$ in the mass range $\sim0.72<$M/M$_{\odot}<2.8$. The tidal radius of the Boc2 cluster ($\sim7-9$) is much more than its observed radius ($\sim1.1$ pc). This suggests that most of the low-mass stars in this cluster are the remains of an older population of stars formed via an earlier epoch of star formation.

The possibility of obtaining symbolic expressions for cosmic backreaction is explored through a case study of so-called 2-region models. By using the publicly available symbolic regression algorithm AI Feynman, it is shown that the kinematical backreaction from a single 2-region model can be well described as a function of the mean redshift (or, equivalently, the volume averaged scale factor). A single expression depending on the redshift/scale factor as well as a model parameter, $f$, that can accurately describe the backreaction for a significant range of models is naturally more complicated but is also achieved with percent-level accuracy. \newline\indent Data sets of redshift drift in the 2-region models are also considered. Again utilizing AI Feynman, expressions for the redshift drift are found. In particular, an expression for the difference between the mean redshift drift and the drift of the mean redshift in terms of the kinematical backreaction is easily obtained for a single 2-region model. An accurate symbolic expression that describes this difference for an array of 2-region models is achieved by using the redshift as a feature instead of the kinematical backreaction.

Symbolic expressions for cosmic backreaction and mean redshift drift in a range of 2-region models in terms of average quantities are presented. The demonstration that these expressions can be obtained constitutes the opening of a new avenue towards understanding the effects of cosmic backreaction in our universe: With a symbolic expression for the redshift drift at hand, the redshift drift can be used to constrain cosmological parameters including the large-scale expansion rate and backreaction. In addition, by introducing symbolic expressions for cosmic backreaction, this quantity can be constrained with observations such as redshift-distance measures.

Super-Earths present compositions dominated by refractory materials. However, there is a degeneracy in their interior structure between a planet with no atmosphere and a small Fe content, and a planet with a thin atmosphere and a higher core mass fraction. To break this degeneracy, atmospheric characterization observations are required. We present a self-consistent interior-atmosphere model to constrain the volatile mass fraction, surface pressure, and temperature of rocky planets with water and CO2 atmospheres. These parameters obtained in our analysis can then be used to predict observations in emission spectroscopy and photometry with JWST, which can determine the presence of an atmosphere, and if present, its composition. To obtain the bolometric emission and Bond albedo for an atmosphere in radiative-convective equilibrium, we present the k-uncorrelated approximation for fast computations within our retrieval on planetary mass, radius and host stellar abundances. For the generation of emission spectra, we use our k-correlated atmospheric model. An adaptive MCMC is used for an efficient sampling of the parameter space at low volatile mass fractions. We show how to use our modelling approach to predict observations with JWST for TRAPPIST-1 c and 55 Cancri e. TRAPPIST-1 c's most likely scenario is a bare surface, although the presence of an atmosphere cannot be ruled out. If the emission in the MIRI F1500 filter is 731 ppm or higher, there would be a water-rich atmosphere. For fluxes between 730 and 400 ppm, no atmosphere is present, while low emission fluxes (300 ppm) indicate a CO2-dominated atmosphere. In the case of 55 Cancri e, a combined spectrum with NIRCam and MIRI LRS may present high uncertainties at wavelengths between 3 and 3.7 $\mu$m. However, this does not affect the identification of H2O and CO2 because they do not present spectral features in this wavelength range.

We present an innovative approach to constraining the non-cold dark matter model using a convolutional neural network (CNN). We perform a suite of hydrodynamic simulations with varying dark matter particle masses and generate mock 21cm radio intensity maps to trace the dark matter distribution. Our proposed method complements the traditional power spectrum analysis. We compare our CNN classification results with those from the power spectrum of the differential brightness temperature map of 21cm radiation, and find that the CNN outperforms the latter. Moreover, we investigate the impact of baryonic physics on the dark matter model constraint, including star formation, self-shielding of HI gas, and UV background model. We find that these effects may introduce some contamination in the dark matter constraint, but they are insignificant when compared to the realistic system noise of the SKA instruments.

We present constraints on models of the galaxy distribution in the cosmic web using a magnitude limited sample from the Galaxy and Mass Assembly (GAMA) survey. We model the redshift-space behaviour of the 2-point correlation function (2pcf) and the recently proposed Voronoi volume function (VVF) -- which includes information beyond 2-point statistics. We extend the standard halo occupation distribution model by introducing extra satellite degrees of freedom, and also by including two assembly bias parameters, $\alpha_{\rm cen}$ and $\alpha_{\rm sat}$, which respectively correlate the occupation numbers of central and satellite galaxies with their host halo's tidal environment. We measure $\alpha_{\rm sat}=1.44^{+0.25}_{-0.43}$ and $\alpha_{\rm cen}=-0.79^{+0.29}_{-0.11}$ using a combination of 2pcf and VVF measurements. These represent a detection of assembly bias at the 3.3$\sigma$ (2.4$\sigma$) significance level for satellite (central) galaxies, a result that remains robust to possible anisotropies in the halo-centric distribution of satellites as well as technicalities of estimating the data covariance. We show that the growth rate ($f\sigma_8$) deduced using models with assembly bias is about 7\% (i.e. $1.5\sigma$) lower than if assembly bias is ignored. Projecting $f\sigma_8$ onto the $\Omega_m$-$\sigma_8$ plane, we find that the model constraints without assembly bias overlap with Planck expectations, but that allowing assembly bias introduces significant tension with Planck, preferring either a lower $\Omega_m$ or a lower $\sigma_8$. We also study the effect of assembly bias on the weak lensing signal. While the all-galaxy lensing signal is unaffected, both central and satellite sub-populations individually show significantly different signals in the presence of assembly bias. [abridged]

Recent work has suggested that black holes (BHs) could be cosmologically coupled to the accelerated expansion of the universe, potentially becoming a candidate for dark energy. This would imply BH mass growth following the cosmological expansion, with the masses of individual BHs growing as $M_{\rm BH}\propto (1+z)^3$. In this letter, we discuss the binary systems Gaia BH1 and Gaia BH2, which contain $\sim 9\,M_{\odot}$ BHs orbited by $\sim 1\,M_{\odot}$ stars in widely-separated orbits. The ages of both systems can be constrained by the properties of the luminous stars. If BH masses are indeed growing as $(1+z)^3$, the masses of both BHs at formation would have been significantly smaller than today. We find a 77% probability that the mass of the BH in Gaia BH2 would have been below $2.2M_\odot$ at formation. This is below the classical Tolman-Oppenheimer-Volkov limit, though it is not yet clear if BHs subject to cosmological coupling should obey this limit. For Gaia BH1, the same probability is 70%. This analysis is consistent with results from two BHs in the globular cluster NGC3201, but unlike the NGC3201 BHs, the Gaia BHs have well-constrained inclinations and thus firm upper mass limits. The discovery of more BHs in binary systems with Gaia astrometry in the coming years will allow us to test the cosmological coupling hypothesis decisively.

The centrifugal force is often omitted in simulations of stellar convection. This force might be important in rapidly rotating stars such as solar analogues due to its $\Omega^2$ scaling, where $\Omega$ is the rotation rate of the star. We study the effects of the centrifugal force in a set of 21 semi-global stellar dynamo simulations with varying rotation rates. Among these, we include three control runs aimed at distinguishing the effects of the centrifugal force from the nonlinear evolution of the solutions. We solve the 3D MHD equations with the Pencil Code in a solar-like convective zone in a spherical wedge setup with a $2\pi$ azimuthal extent. We decompose the magnetic field in spherical harmonics and study the migration of azimuthal dynamo waves (ADWs), energy of different large-scale magnetic modes, and differential rotation. In the regime with the lowest rotation rates, $\Omega = 5-10\Omega_\odot$, where $\Omega_\odot$ is the rotation rate of the Sun, we see no marked changes in neither the differential rotation nor the magnetic field properties. For intermediate rotation with $\Omega = 20-25\Omega_\odot$ we identify an increase of the differential rotation as a function of centrifugal force. The axisymmetric magnetic energy tends to decrease with centrifugal force while the non-axisymmetric one increases. The ADWs are also affected, especially the propagation direction. In the most rapidly rotating set with $\Omega=30\Omega_\odot$, these changes are more pronounced and in one case the propagation direction of the ADW changes from prograde to retrograde. Control runs suggest that the results are a consequence of the centrifugal force and not due to the details of the initial conditions or the history of the run. We find that the differential rotation and properties of the ADWs change as a function of the centrifugal force only when rotation is rapid enough.

We present an intriguing, serendipitously-detected system consisting of an S0/a galaxy, which we refer to as the "Kite", and a highly-collimated tail of gas and stars that extends over 380 kpc and contains pockets of star formation. In its length, narrowness, and linearity the Kite's tail is an extreme example relative to known tails. The Kite (PGC 1000273) has a companion galaxy, Mrk 0926 (PGC 070409), which together comprise a binary galaxy system in which both galaxies host active galactic nuclei. Despite this systems being previously searched for signs of tidal interactions, the tail had not been discovered prior to our identification as part of the validation process of the SMUDGes survey for low surface brightness galaxies. We confirm the kinematic association between various H$\alpha$ knots along the tail, a small galaxy, and the Kite galaxy using optical spectroscopy obtained with the Magellan telescope and measure a velocity gradient along the tail. The Kite shares characteristics common to those formed via ram pressure stripping ("jellyfish" galaxies) and formed via tidal interactions. However, both scenarios face significant challenges that we discuss, leaving open the question of how such an extreme tail formed. We propose that the tail resulted from a three-body interaction from which the lowest-mass galaxy was ejected at high velocity.

The nature of the merger remnant of binary neutron star (BNS) remains an open question. From the theoretical point of view, one possible outcome is a supra-massive neutron star (SMNS), which is supported by rigid rotation and through its survival of hundreds of seconds before collapsing into a black hole (BH). If this is the case, the SMNS can emit continuous gravitational waves (GW) and electromagnetic (EM) radiation, particularly in the X-ray band. In this work, the ellipticity and initial frequency of SMNS are constrained with a Bayesian framework using simulated X-ray and GW signals, which could be detected by The Transient High Energy Sky and Early Universe Surveyor (THESEUS) and Einstein Telescope (ET), respectively. We found that only considering the X-ray emission can not completely constrain the initial frequency and ellipticity of the SMNS, but it can reduce the ranges of the parameters. Afterwards, we can use the posterior distribution of the X-ray parameter estimates as a prior for the GW parameter estimates. It was found that the 95$\%$ credible region of the joint X-ray-GW analysis was about $10^5$ times smaller than that of the X-ray analysis alone.

The application of VLBI to scientific problems has undergone a relentless expansion since its conception, yet the potential for further expansion is still large. We are on the cusp of revolutionary progress given the arrival of a host of next-generation instruments. Over the last few years the community has been working hard to ensure the SKA design includes the capability to enable multiple simultaneous tied-array beams, which is a crucial technology to deliver ultra-precise astrometry and improve survey speed capabilities. However, to reach the full potential requires that the network of antennas is upgraded to match the SKA capabilities. We identify multiple-pixel technology, on large telescopes and connected arrays, as a crucial missing component and here will make recommendations for the upgrade path of the partner EVN (and other network) telescopes. Our feasibility studies on SKA-VLBI suggest an order of magnitude improvement in the precision and also in the frequency range at which astrometry can be performed today, if the full network has the required capabilities.

We explore an extension of the standard $\Lambda$CDM model by including an interaction between neutrinos and dark matter, and making use of the ground based telescope data of the Cosmic Microwave Background (CMB) from the Atacama Cosmology Telescope (ACT). An indication for a non-zero coupling between dark matter and neutrinos (both assuming a temperature independent and $T^2$ dependent cross-section) is obtained at the 1$\sigma$ level coming from the ACT CMB data alone and when combined with the Planck CMB and Baryon Acoustic Oscillations (BAO) measurements. This result is confirmed by both fixing the effective number of relativistic degrees of freedom in the early Universe to the Standard Model value of $N_{\rm eff}=3.044$, and allowing $N_{\rm eff}$ to be a free cosmological parameter. Furthermore, when performing a Bayesian model comparison, the interacting $\nu$DM (+$N_{\rm eff}$) scenario is mostly preferred over a baseline $\Lambda$CDM (+$N_{\rm eff}$) cosmology. The preferred value is then used as a benchmark and the potential implications of dark matter's interaction with a sterile neutrino are discussed.

We report a new imaging spectroscopic observation of Oort-cloud comet C/2017 K2 (hereafter K2) on its way to perihelion at 2.53 au, around a heliocentric distance where H2O ice begins to play a key role in comet activation. Normalized reflectances over 6 500--8 500 AA for its inner and outer comae are 9.7+/-0.5 and 7.2+/-0.3 % (10^3 AA)^-1, respectively, the latter being consistent with the slope observed when the comet was beyond the orbit of Saturn. The dust coma at the time of observation appears to contain three distinct populations: mm-sized chunks prevailing at <~10^3 km; a 10^5-km steady-state dust envelope; and fresh anti-sunward jet particles. the dust chunks dominate the continuum signal and are distributed over a similar radial distance scale as the coma region with redder dust than nearby. they also appear to be co-spatial with OI1D, suggesting that the chunks may accommodate H2O ice with a fraction (>~1 %) of refractory materials. The jet particles do not colocate with any gas species detected. The outer coma spectrum contains three significant emissions from C2(0,0) Swan band, OI1D, and CN(1,0 red band, with an overall deficiency in NH2. Assuming that all OI1D flux results from H2O dissociation, we compute an upper limit on the water production rate Q_H2O of ~7 x 10^28 molec s^-1 (with an uncertainty of a factor of two). the production ratio log[Q_C2/Q_CN] of K2 suggests that the comet has typical carbon-chain composition, with the value potentially changing with distance from the Sun. Our observations suggest that water ice-containing dust chunks (>0.1 mm) near K2's nucleus emitted beyond 4 au may be responsible for its very low gas rotational temperature and the discrepancy between its optical and infrared lights reported at similar heliocentric distances.

The atomic hydrogen (HI) properties of star-forming galaxies in the local Universe are known to correlate with other galaxy properties via the ``HI scaling relations''. The redshift evolution of these relations serves as an important constraint on models of galaxy evolution. However, until recently, there were no estimates of the HI scaling relations at cosmological distances. Using data from a deep Giant Metrewave Radio Telescope HI 21 cm survey of the Extended Groth Strip, and the technique of spectral line stacking, we determine the scaling relation between the HI mass and the stellar mass for star-forming galaxies at $z\approx0.35$. We use this measurement, along with the main-sequence relation in galaxies, to infer the dependence of the HI depletion timescale of these galaxies on their stellar mass. We find that massive star-forming galaxies at $z\approx0.35$, with stellar mass $\rm M_* \gtrsim10^{9.5}\:M_{\odot}$, are HI-poor compared to local star-forming galaxies of a similar stellar mass. However, their characteristic HI depletion time is lower by a factor of $\approx 5$ than that of their local analogues, indicating a higher star-formation efficiency at intermediate redshifts (similar to that at $z \approx 1$). While our results are based on a relatively small cosmic volume and could thus be affected by cosmic variance, the short characteristic HI depletion timescales ($\lesssim 3$ Gyr) of massive star-forming galaxies at $z \approx 0.35$ indicate that they must have acquired a significant amount of neutral gas through accretion from the circumgalactic medium over the past four Gyr, to avoid quenching of their star-formation activity.

On the Large Phased Array (LPA) of Lebedev Physics Institute (LPI), a search for pulsars outside the Galaxy plane was carried out in a 300 sq. deg area. The search with a sensitivity 5-10 times better than that of previously conducted surveys was at a frequency of 111 MHz. The search was carried out in the summed power spectra. With an accumulation equivalent to 100 hours of continuous observations for each point of the area, 5 known pulsars were detected with a signal-to-noise ratio (S/N) from 20 to 1300 in the first harmonic of the spectrum. Average profiles were obtained for the detected pulsars. Estimates of the peak and integral flux densities of the found pulsars are given for individual sessions and for the power spectra summarized over 5.5 years, obtained using the developed method based on measurements of the height of harmonics in the power spectrum. No new pulsars have been detected in the area. Apparently, when searching for pulsars in the area, we have approached the lower limit of the luminosity of the second pulsars. The completeness of the survey is at the level of 0.5 mJy.

We present the methodology and results of a simulation to determine the recoverability of LEO objects using a blind stacking technique. The method utilises sCMOS and GPU technology to inject and recover LEO objects in real observed data. We explore the target recovery fraction and pipeline run-time as a function of three optimisation parameters; number of frames per data-set, exposure time, and binning factor. Results are presented as a function of magnitude and velocity. We find that target recovery using blind stacking is significantly more complete, and can reach fainter magnitudes, than using individual frames alone. We present results showing that, depending on the combination of optimisation parameters, recovery fraction is up to 90% of detectable targets for magnitudes up to 13.5, and then falls off steadily up to a magnitude limit around 14.5. Run-time is shown to be a few multiples of the observing time for the best combinations of optimisation parameters, approaching real-time processing.

The Einstein Telescope (ET) is a third generation gravitational wave detector that includes a room-temperature high-frequency (ET-HF) and a cryogenic low-frequency laser interferometer (ET-LF). The cryogenic ET-LF is crucial for exploiting the full scientific potential of ET. We present a new baseline design for the cryogenic payload that is thermally and mechanically consistent and compatible with the design sensitivity curve of ET. The design includes two options for the heat extraction from the marionette, based on a monocrystalline high-conductivity marionette suspension fiber and a thin-wall titanium tube filled with static He-II, respectively. Following a detailed description of the design options and the suspension thermal noise (STN) modelling, we present the sensitivity curves of the two baseline designs, discuss the influence of various design parameters on the sensitivity of ET-LF and conclude with an outlook to future R&D activities.

The time delay between images of strongly gravitationally lensed quasars is an established cosmological probe. Its limitations, however, include uncertainties in the assumed mass distribution of the lens. We re-examine the methodology of a prior work presenting a geometric probe of cosmology independent of the lensing potential which considers differential time delays over images, originating from spatially-separated photometric signals within a strongly lensed quasar. We give an analytic description of the effect of the differential lensing on the emission line spectral flux for axisymmetric Broad Line Region geometries, with the inclined ring or disk, spherical shell, and double cone as examples. The proposed method is unable to recover cosmological information as the observed time delay and inferred line-of-sight velocity do not uniquely map to the three-dimensional position within the source.

Recently much progress has been made in probing the embedded stages of massive star formation, pointing to formation scenarios akin to a scaled up version of low-mass star formation. However, the latest stages of massive star formation have rarely been observed. Using 1st and 2nd overtone CO bandhead emission and near- to mid-infrared photometry we aim to characterize the remnant formation disks around 5 unique pre-main-sequence (PMS) stars with masses $6-12~\rm M_{\odot}$, that have constrained stellar parameters thanks to their detectable photospheres. We seek to understand this emission and the disks it originates from in the context of the evolutionary stage of the studied sources. We use an analytic LTE disk model to fit the CO bandhead and the dust emission, found to originate in different disk regions. For the first time we modeled the 2nd overtone emission. Furthermore, we fit continuum normalized bandheads and show the importance of this in constraining the emission region. We also include $^{13}\rm CO$ in our models as an additional probe of the young nature of the studied objects. We find that the CO emission originates in a narrow region close to the star (<1 AU) and under very similar disk conditions (temperatures and densities) for the different objects. This is consistent with previous modeling of this emission in a diverse range of young stellar objects. We discuss these results in the context of the positions of these PMS stars in the Hertzsprung-Russel diagram and the CO emission's association with early age and high accretion rates in (massive) young stellar objects. We conclude that, considering their mass range and for the fact that their photospheres are detected, the M17 PMS stars are observed in a relatively early formation stage. They are therefore excellent candidates for longer wavelength studies to further constrain the end stages of massive star formation.

Fundamental cosmology observations, such as the detection of the redshift drift and tests of the universality of physical laws, are key science and design drivers of the ArmazoNes high Dispersion Echelle Spectrograph (ANDES), an Extremely Large Telescope instrument. While separate forecasts for each of them have been reported, we have developed Fisher Matrix based forecast tools combining both of these observables. We demonstrate the synergies between the two ANDES datasets, quantifying the improvements in cosmology and fundamental physics parameter constraints for two separate theoretical paradigms. We publicly release this forecast code, which is one of the tools for the optimisation of the ANDES observing strategy.

Constraints on the interior structure of the Moon have been derived from its inductive response, principally as measured by the magnetic transfer function (TF) between the distantly orbiting Explorer 35 satellite and the Apollo 12 surface station. The most successful prior studies used a dataset 0.01-1 mHz, so the lunar response could be modeled as a simple dipole. However, earlier efforts also produced transfer functions up to 40 mHz. The smaller electromagnetic skin depth at higher frequency would better resolve the uppermost mantle - where key information about primitive lunar evolution may still be preserved - but requires a multipole treatment. I compute new profiles of electrical conductivity vs depth using both the low-frequency and the full-bandwidth ranges of published Apollo-Explorer TFs. I derive temperature profiles at depths >400 km (<1 mHz) consistent with conductive heat loss and expectations of the iron (and possibly water) content of the mantle. The near-constant iron fraction (Mg# 81 +/- 7) implies either efficient mixing, due to now-defunct convection or perhaps incomplete overturn of gravitationally unstable cumulates following crystallization of the magma ocean. In contrast, the full-bandwidth analysis produced a different conductivity profile that could not be realistically matched by conduction, convection, partial melting, or simple considerations of lateral heterogeneity. I conclude that the TF method at the Moon is unreliable >>1 mHz. Future EM sounding using the magnetotelluric method can operate up to 100s Hz and is largely insensitive to multipole effects, resolving structure to 100 km or less.

ALMA has spatially resolved over 200 annular structures in protoplanetary discs, many of which are suggestive of the presence of planets. Constraining the mass of these putative planets is quite degenerate for it depends on the disc physical properties, and for simplicity a steady-state is often assumed whereby the planet position is kept fixed and there is a constant source of dust at the outer edge of the disc. Here we argue against this approach by demonstrating how the planet and dust dynamics can lift degeneracies of such steady-state models. We take main disc parameters from the well-known protoplanetary disc HD 163296 with a suspected planet at $R\approx 86$~au as an example. By running gas and dust hydrodynamical simulations post-processed with dust radiative transfer calculations, we first find steady-state disc and planet parameters that reproduce ALMA continuum observations fairly well. For the same disc mass, but now allowing the planet to migrate in the simulation, we find that the planet undergoes runaway migration and reaches the inner disc in $\sim 0.2$ Myr. Further, decreasing the disc mass slows down planet migration, but it then also increases the dust's radial drift, thereby depleting the disc dust faster. We find that the opposing constraints of planet migration and dust drift require the disc mass to be at most $0.025~\msun$, must less massive than previously estimated, and for the dust to be porous rather than compact. We propose that similar analysis should be extended to other sources with suspected planetary companions.

Gravitational redshift is a fundamental parameter that allows us to determine the mass-to-radius ratio of compact stellar objects, such as black holes, neutron stars, and white dwarfs (WDs). In the X-ray spectra of the close binary system, RX J1712.6$-$2414, obtained from the Chandra High-Energy Transmission Grating observation, we detected significant redshifts for characteristic X-rays emitted from hydrogen-like magnesium, silicon ($\Delta E/E_{\rm rest} \sim 7 \times 10^{-4}$), and sulfur ($\Delta E/E_{\rm rest} \sim 15 \times 10^{-4}$) ions, which are over the instrumental absolute energy accuracy (${\Delta E/E_{\rm rest} \sim 3.3} \times 10^{-4}$). Considering some possible factors, such as Doppler shifts associated with the plasma flow, systemic velocity, and optical depth, we concluded that the major contributor to the observed redshift is the gravitational redshift of the WD harbored in the binary system, which is the first gravitational redshift detection from a magnetic WD. Moreover, the gravitational redshift provides us with a new method of the WD mass measurement by invoking the plasma-flow theory with strong magnetic fields in close binaries. Regardless of large uncertainty, our new method estimated the WD mass to be $M_{\rm WD}> 0.9\,M_{\odot}$.

Gamma-ray bursts (GRBs) have long been suspected as possible ultra-high energy cosmic ray (UHECR) accelerators. In this brief note, I report that GRB 980425/SN 1998bw falls within the region of interest (ROI) with the highest significance in an all-sky blind search for magnetically-induced effects in the arrival directions of UHECRs conducted by the Pierre Auger Collaboration with events detected up to to 2018 August 31. There is also report in the literature of delayed TeV emission in archival Fermi-LAT observations from the direction of GRB 980425/SN 1998bw. The combined probability that two distinct cosmic ray acceleration signatures in two different multimessenger experiments may appear at the same spatial location by chance is estimated to be between $1.62 \times 10^{-3}$ and 0.0157.

We study the one-loop corrections in power spectrum of long gravitational waves induced from small scale modes in the models of single field inflation undergoing a phase of ultra-slow-roll (USR). We show that the spectrum of long tensor perturbations are largely unaffected by the loop corrections from the short scalar modes. In particular, the spectrum of long tensor perturbations is insensitive to the sharpness of the transition from the USR phase to the final slow-roll phase. This is in contrast to the case of scalar power spectrum in which the loop corrections can be large for a sharp transition while it is slow-roll suppressed in a mild transition. We study the tensor-scalar-scalar bispectrum in the squeezed limit and demonstrate that the Maldacena consistency condition does hold.

We present an analytical and numerical study of a system composed of a stellar binary pair and a massless, locally isothermal viscous accretion disk that is coplanar to the binary orbital plane. Analytically, we study the effect of the binary's gravitational potential over short timescales through the study of stability for epicyclic orbits, and over long timescales by revisiting the concept of resonant torques. Numerically, we perform two-dimensional Newtonian numerical simulations of the disk-binary system over a range of binary mass ratios. We find that the results of our simulations are consistent with previous numerical studies. We additionally show, by comparison of the analytical and numerical results, that the circumbinary gap is maintained on the orbital timescale through the driving of epicyclic instabilities, and does not depend on resonant torquing, contrary to standard lore. While our results are applicable to any disk-binary system, we highlight the importance of this result in the search for electromagnetic and gravitational-wave signatures from supermassive black-hole binaries.

In an ideal germanium detector, fully-absorbed monoenergetic gamma-rays will appear in the measured spectrum as a narrow peak, broadened into a Gaussian of width determined only by the statistical properties of charge cloud generation and the electronic noise of the readout electronics. Multielectrode detectors complicate this picture. Broadening of the charge clouds as they drift through the detector will lead to charge sharing between neighboring electrodes and, inevitably, low-energy tails on the photopeak spectra. We simulate charge sharing in our germanium cross strip detectors in order to reproduce the low-energy tails due to charge sharing. Our goal is to utilize these simulated spectra to develop an analytical fit (shape function) for the spectral lines that provides a robust and high-quality fit to the spectral profile, reliably reproduces the interaction energy, noise width, and the number of counts in both the true photopeak and the low-energy tail, and minimizes the number of additional parameters. Accurate modeling of the detailed line profiles is crucial for both calibration of the detectors as well as scientific interpretation of measured spectra.

The slitless grism on the Nancy Grace Roman Space Telescope will enable deep near-infrared spectroscopy over a wide field of view. We demonstrate Roman's capability to detect Ly$\alpha$ galaxies at $z>7$ using a multi-position-angle (PA) observational strategy. We simulate Roman grism data using a realistic foreground scene from the COSMOS field. We also input fake Ly$\alpha$ galaxies spanning redshift z=7.5-10.5 and a line-flux range of interest. We show how a novel data cube search technique -- CUBGRISM -- originally developed for GALEX can be applied to Roman grism data to produce a Ly$\alpha$ flux-limited sample without the need for continuum detections. We investigate the impact of altering the number of independent PAs and exposure time. A deep Roman grism survey with 25 PAs and a total exposure time of $70$hrs can achieve Ly$\alpha$ line depths comparable to the deepest $z=7$ narrow-band surveys ($L_{{\rm{Ly}}\alpha}\gtrsim10^{43}$erg s$^{-1}$). Assuming a null result, where the opacity of the intergalactic medium (IGM) remains unchanged from $z\sim7$, this level of sensitivity will detect $\sim400$ deg$^{-2}$ Ly$\alpha$ emitters from $z=7.25-8.75$. A decline from this expected number density is the signature of an increasing neutral hydrogen fraction and the onset of reionization. Our simulations indicate that a deep Roman grism survey has the ability to measure the timing and magnitude of this decline, allowing us to infer the ionization state of the IGM and helping us to distinguish between models of reionization.

PySR is an open-source library for practical symbolic regression, a type of machine learning which aims to discover human-interpretable symbolic models. PySR was developed to democratize and popularize symbolic regression for the sciences, and is built on a high-performance distributed back-end, a flexible search algorithm, and interfaces with several deep learning packages. PySR's internal search algorithm is a multi-population evolutionary algorithm, which consists of a unique evolve-simplify-optimize loop, designed for optimization of unknown scalar constants in newly-discovered empirical expressions. PySR's backend is the extremely optimized Julia library SymbolicRegression.jl, which can be used directly from Julia. It is capable of fusing user-defined operators into SIMD kernels at runtime, performing automatic differentiation, and distributing populations of expressions to thousands of cores across a cluster. In describing this software, we also introduce a new benchmark, "EmpiricalBench," to quantify the applicability of symbolic regression algorithms in science. This benchmark measures recovery of historical empirical equations from original and synthetic datasets.

Constraining cosmological parameters from large-scale structure observations requires precise and accurate tools to compute its properties. While perturbation theory (PT) approaches can serve this purpose, exploration of large parameter space is challenging due to the potentially large computational cost of such calculations. In this study, we show that a response function approach applied to the regularized PT (RegPT) model at the 2-loop order, plus correction terms induced by redshift space distortion effects, can reduce the runtime by a factor of 50 compared to direct integration. We illustrate the performance of this approach by performing the parameter inference of five fundamental cosmological parameters from the redshift space power spectrum measured from $N$-body simulations as mock measurements, and inferred cosmological parameters are directly compared with parameters used to generate initial conditions of the simulations. From this \textit{PT challenge} analysis, the constraining power of cosmological parameters and parameter biases are quantified with the survey volume and galaxy number density expected for the \textit{Euclid} mission at the redshift $z=1$ as a function of the maximum wave-number of data points $k_\mathrm{max}$. We find that RegPT with correction terms reproduces the input cosmological parameters without bias up to maximum wave-number $k_\mathrm{max} = 0.18 \, h\,\mathrm{Mpc}^{-1}$. Moreover, RegPT+, which introduces one free parameter to RegPT to handle the damping feature on small scales, delivers the best performance among the examined models and achieves tighter constraints without significant parameter bias for higher maximum wave-number $k_\mathrm{max} = 0.21 \, h\,\mathrm{Mpc}^{-1}$.

This paper presents an analysis of XMM X--ray spectra of the quasar 1ES 1553+113, in search for absorption lines from the intervening warm--hot intergalactic medium. A search for OVII, OVIII and NeIX resonance absorption lines was performed at eight fixed redshifts that feature OVI or HI broad Lyman--$\alpha$ absorption lines that were previously detected from HST data. The search yielded one possible detection of OVII at a redshift z=0.1877 with an OVI prior, with a statistical significance that is equivalent to a 2.6-$\sigma$ confidence level. The spectra were also stacked at the wavelengths of the expected redshifted OVII and OVIII lines, but the analysis did not reveal evidence for the presence of additional X--ray absorbing WHIM. Moreover, the spectra were used to investigate two putative OVII absorption lines that were detected serendipitously in an earlier analysis of the same data by F. Nicastro and collaborators. The paper also presents a comprehensive statistical framework for cosmological inferences from the analysis of absorption lines, which makes use of cosmological simulations for the joint probability distributions of FUV and X--ray ions. Accordingly, we conclude that the new possible OVII absorption at z=0.1877 is consistent with a contribution from the hot WHIM to the baryon density in an amount of $\Omega_{WHIM,X}/\Omega_b = 44\pm22$\%. However, there are large systematic uncertainties associated with the temperature and abundances of the absorbers, and only a larger sample of X-ray sources can provide an accurate determination of the cosmological density of the WHIM.

I present a class of theories that generalize quasilinear MOND (QUMOND). Like QUMOND, these GQUMOND theories require solving only the linear Poisson equation (twice). Unlike QUMOND, their Lagrangian depends on higher derivatives of the Newtonian potential. They thus dictate different ``phantom'' densities as virtual sources in the Poisson equation for the MOND potential. These theories might open new avenues to more fundamental theories, and have much heuristic value. I use them to demonstrate that even within limited classes of modified-gravity formulations of MOND, theories can differ substantially on lower-tier MOND predictions. Such GQUMOND theories force, generically, the introduction of dimensioned constants other than the MOND acceleration, $a_0$, such as a length, a frequency, etc. As a result, some of these theories reduce to QUMOND itself only, e.g., on length scales (or, in other versions, dynamical times) larger than some critical value. But in smaller systems (or, alternatively, in ones with shorter dynamical times), MOND effects are screened, even if their internal accelerations are smaller than $a_0$. In such theories it is possible that MOND (expressed as QUMOND) applies on galactic scales, but its departures from Newtonian dynamics are substantially suppressed in some subgalactic systems -- such as binary stars, and open, or globular star clusters. The same holds for the effect of the galactic field on dynamics in the inner solar system, which can be greatly suppressed compared with what QUMOND predicts. Tidal effects of a galaxy on smaller subsystems are the same as in QUMOND, for the examples I consider. I also describe briefly versions that do not involve dimensioned constants other than $a_0$, and yet differ from QUMOND in important ways.

Earth-side observations of solar p modes can be used to image and monitor magnetic activity on the Sun's far side. Here we use magnetograms of the far side obtained by the Polarimetric and Helioseismic Imager (PHI) onboard Solar Orbiter (SO) to directly assess -- for the first time -- the validity of far-side helioseismic holography. We wish to co-locate the positions of active regions in helioseismic images and magnetograms, and to calibrate the helioseismic measurements in terms of magnetic field strength. We identify three magnetograms on 18 November 2020, 3 October 2021, and 3 February 2022 displaying a total of six active regions on the far side. The first two dates are from SO's cruise phase, the third from the beginning of the nominal operation phase. We compute contemporaneous seismic phase maps for these three dates using helioseismic holography applied to time series of Dopplergrams from the Helioseismic and Magnetic Imager (HMI) on the Solar Dynamics Observatory (SDO). Among the six active regions seen in SO/PHI magnetograms, five active regions are identified on the seismic maps at almost the same positions as on the magnetograms. One region is too weak to be detected above the seismic noise. To calibrate the seismic maps, we fit a linear relationship between the seismic phase shifts and the unsigned line-of-sight magnetic field averaged over the active region areas extracted from the SO/PHI magnetograms. SO/PHI provides the strongest evidence so far that helioseismic imaging provides reliable information about active regions on the far side, including their positions, areas, and mean unsigned magnetic field.

One of the most reliable means of studying the stellar interior is through the apsidal motion in double line eclipsing binary systems since these systems present errors in masses, radii, and effective temperatures of only a few per cent. On the other hand, the theoretical values of the apsidal motion to be compared with the observed values depend on the stellar masses of the components and more strongly on their radii (fifth power).The main objective of this work is to make available grids of evolutionary stellar models that, in addition to the traditional parameters (e.g. age, mass, log g, T$_{\rm eff}$), also contain the necessary parameters for the theoretical study of apsidal motion and tidal evolution. This information is useful for the study of the apsidal motion in eclipsing binaries and their tidal evolution, and can also be used for the same purpose in exoplanetary systems. All models were computed using the MESA package. We consider core overshooting for models with masses $\ge$ 1.2 M$_\odot$. For the amount of core overshooting we adopted a recent relationship for mass $\times$ core overshooting. We adopted for the mixing-length parameter $\alpha_{\rm MLT}$ the value 1.84 (the solar-calibrated value). Mass loss was taken into account in two evolutionary phases. The models were followed from the pre-main sequence phase to the white dwarf (WD) stage.The evolutionary models containing age,luminosity, log g, and Teff, as well as the first three harmonics of the internal stellar structure (k$_2$, k$_3$, and k$_4$), the radius of gyration $\beta$ y, and the dimensionless variable $\alpha$, related to gravitational potential energy, are presented in 69 tables covering three chemical compositions: [Fe/H] = -0.50, 0.00, and 0.50. Additional models with different input physics are available.

The second postulate of special relativity states that the speed of light in vacuum is independent of the emitter's motion. Though this claim has been verified in various experiments and observations involving electromagnetic radiation with very high accuracy, such a test for gravitational radiation still needs to be explored. We analyzed data from the LIGO and Virgo detectors to test this postulate for gravitational radiation within the ambit of \textit{emission models}, where the speed of gravitational waves emitted by a source moving with a velocity $v$ relative to a stationary observer is given by ${c' = c + k\,v}$, where $k$ is a constant. We have estimated the upper bound on the 90\% credible interval over $k$ that parameterizes the deviation from the second postulate to be ${k \leq 8.3 \times {10}^{-18}}$ which is several orders of magnitude more stringent compared to previous bounds obtained from electromagnetic observations. The Bayes' factor supports the second postulate, with very strong evidence that the data is consistent with the null hypothesis $k = 0$. This confirms that the speed of gravity is independent of the motion of the emitter, upholding the principle of relativity for gravitational interactions.

We use multimessenger observations of the neutron star merger event GW170817 to derive new constraints on axion-like particles (ALPs) coupling to photons. ALPs are produced via Primakoff and photon coalescence processes in the merger, escape the remnant and decay back into two photons, giving rise to a photon signal approximately along the line-of-sight to the merger. We analyze the spectral and temporal information of the ALP-induced photon signal, and use the Fermi-LAT observations of GW170817 to derive our new ALP constraints. We also show the improved prospects with future MeV gamma-ray missions, taking the spectral and temporal coverage of AMEGO-X as an example.

We show via both analytical calculation and numerical simulation that the optimal cross-correlation statistic (OS) for stochastic gravitational-wave-background (GWB) searches using data from pulsar timing arrays follows a generalized chi-squared (GX2) distribution-i.e., a linear combination of chi-squared distributions with coefficients given by the eigenvalues of the quadratic form defining the statistic. This observation is particularly important for calculating the frequentist statistical significance of a possible GWB detection, which depends on the exact form of the distribution of the OS signal-to-noise ratio (S/N) $\hat\rho \equiv \hat A_{\rm gw}^2/\sigma_0$ in the absence of GW-induced cross correlations (i.e., the null distribution). Previous discussions of the OS have incorrectly assumed that the analytic null distribution of $\hat\rho$ is well-approximated by a zero-mean unit-variance Gaussian distribution. Empirical calculations show that the null distribution of $\hat\rho$ has "tails" which differ significantly from those for a Gaussian distribution, but which follow (exactly) a GX2 distribution. So, a correct analytical assessment of the statistical significance of a potential detection requires the use of a GX2 distribution.

In extreme astrophysical environments such as core-collapse supernovae and binary neutron star mergers, neutrinos play a major role in driving various dynamical and microphysical phenomena, such as baryonic matter outflows, the synthesis of heavy elements, and the supernova explosion mechanism itself. The interactions of neutrinos with matter in these environments are flavor-specific, which makes it of paramount importance to understand the flavor evolution of neutrinos. Flavor evolution in these environments can be a highly nontrivial problem thanks to a multitude of collective effects in flavor space, arising due to neutrino-neutrino ($\nu$-$\nu$) interactions in regions with high neutrino densities. A neutrino ensemble undergoing flavor oscillations under the influence of significant $\nu$-$\nu$ interactions is somewhat analogous to a system of coupled spins with long-range interactions among themselves and with an external field ('long-range' in momentum-space in the case of neutrinos). As a result, it becomes pertinent to consider whether these interactions can give rise to significant quantum correlations among the interacting neutrinos, and whether these correlations have any consequences for the flavor evolution of the ensemble. In particular, one may seek to utilize concepts and tools from quantum information science and quantum computing to deepen our understanding of these phenomena. In this article, we attempt to summarize recent work in this field. Furthermore, we also present some new results in a three-flavor setting, considering complex initial states.

In this work, we calculate the average equation of state for the oscillating inflaton field of the simplest $\alpha$-attractor E-model. We show that the average equation of state can be solved analytically. We discover that when $\alpha$ is small, the average equation of state of the oscillating inflaton field approaches that of a cosmological constant. This is the phenomenon of oscillating inflation.

We consider cosmology with an inflaton scalar field with an additional quartic kinetic term. Such a theory can be motivated by Palatini $R+R^2$ modified gravity. Assuming a runaway inflaton potential, we take the Universe to become dominated by the kinetic energy density of the scalar field after inflation. Initially, the leading kinetic term is quartic and we call the corresponding period hyperkination. Subsequently, the usual quadratic kinetic term takes over and we have regular kination, until reheating. We study, both analytically and numerically, the spectrum of primordial gravitational waves generated during inflation and re-entering the horizon during the subsequent eras. We demonstrate that the spectrum is flat for modes re-entering during radiation domination and hyperkination and linear in frequency for modes re-entering during kination: kinetic domination boosts the spectrum, but hyperkination truncates its peak. As a result, the effects of the kinetic period can be extended to observable frequencies without generating excessive gravitational waves, which could otherwise destabilise the process of Big Bang Nucleosynthesis. We show that there is ample parameter space for the primordial gravitational waves to be observable in the near future. If observed, the amplitude and `knee' of the spectrum will provide valuable insights into the background theory.

We discuss the production of primordial black holes in an early matter dominated era, which typically takes place in string inspired early universe cosmological models. In particular, we consider a pre-big bang scenario (extending previous results regarding formation in the radiation dominated era) where the enhancement of curvature perturbations is induced by a variation of the sound-speed parameter c_s during the string phase of high-curvature inflation. After imposing all relevant observational constraints, we find that the considered class of models is compatible with the production of a large amount of primordial black holes, in the mass range relevant to dark matter, only for a small range of the parameters space. On the other hand, we find that a huge production of light primordial black holes may occur both in such matter dominated era and in the radiation dominated one.

Numerical simulations of merging compact objects and their remnants form the theoretical foundation for gravitational wave and multi-messenger astronomy. While Cartesian-coordinate-based adaptive mesh refinement is commonly used for simulations, spherical-like coordinates are more suitable for nearly spherical remnants and azimuthal flows due to lower numerical dissipation in the evolution of fluid angular momentum, as well as requiring fewer numbers of computational cells. However, the use of spherical coordinates to numerically solve hyperbolic partial differential equations can result in severe Courant-Friedrichs-Lewy (CFL) stability condition timestep limitations, which can make simulations prohibitively expensive. This paper addresses this issue for the numerical solution of coupled spacetime and general relativistic magnetohydrodynamics evolutions by introducing a double FFT filter and implementing it within the fully MPI-parallelized SphericalNR framework in the Einstein Toolkit. We demonstrate the effectiveness and robustness of the filtering algorithm by applying it to a number of challenging code tests, and show that it passes these tests effectively, demonstrating convergence while also increasing the timestep significantly compared to unfiltered simulations.

We outline a partial historical summary of the steps through which the nucleosynthesis phenomena induced by {\it slow} neutron captures (the {\it s-process}) were clarified, a scientific achievement in which Franz K\"appeler played a major role. We start by recalling the early phenomenological approach, which yielded a basic understanding of the subject even before models for the parent stellar evolutionary stages were developed. Through such a tool, rough limits for the neutron density and exposure were set, and the crucial fact was understood that more than one nucleosynthesis component is required to account for solar abundances of $s$-process nuclei up to the Pb-Bi region. We then summarize the gradual understanding of the stellar processes actually involved in the production of nuclei from Sr to Pb (the so-called {\it Main Component}, achieved in the last decade of the past century and occurring in red giants of low and intermediate mass, ($M \lesssim$ 8 $M_{\odot}$), populating, in the {\it HR} diagram, the {\it Asymptotic Giant Branch} or {\it AGB} region. We conclude by giving some details on more recent research concerning mixing mechanisms inducing the activation of the main neutron source, $^{13}$C($\alpha$,n)$^{16}$O.

The anisotropies of the Cosmological Gravitational Wave Background (CGWB) retain information about the primordial mechanisms that source the gravitational waves and about the geometry and the particle content of the universe at early times. In this work, we discuss in detail the computation of the angular power spectra of CGWB anisotropies and of their cross correlation with Cosmic Microwave Background (CMB) anisotropies, assuming different processes for the generation of these primordial signals. We present an efficient implementation of our results in a modified version of CLASS which will be publicly available. By combining our new code GW_CLASS with MontePython, we forecast the combined sensitivity of future gravitational wave interferometers and CMB experiments to the cosmological parameters that characterize the cosmological gravitational wave background.

In this paper we study the viability of an entropic cosmological model. The effects of entropic gravity are derived from a modified entropy-area relationship with a volumetric entropy term. This model describes a late time limit {cosmic acceleration}, whose origin is related to a volumetric term in the entropy. Moreover, we analyze the phenomenological implications of the entropic model using the Supernovae {\it Pantheon} compilation and the observational Hubble parameter data to find consistency with cosmological observations. Finally, we show the equivalence between the entropic model and a brane world cosmological model, by means of an effective geometrical construction.