Papers for Wednesday, Nov 01 2023

We describe several updates to Korg, a package for 1D LTE spectral synthesis of FGKM stars. Built-in functions to fit observed spectra via synthesis or equivalent widths make it easy to take advantage of Korg's automatic differentiation. Comparison to a past analysis of 18 Sco shows that we obtain significantly reduced line-to-line abundance scatter with Korg. Fitting and synthesis are facilitated by a rigorously-tested model atmosphere interpolation method, which introduces negligible error to synthesized spectra for stars with $T_\mathrm{eff} \gtrsim 4000\,\mathrm{K}$. For cooler stars, atmosphere interpolation is complicated by the presence of molecules, though we demonstrate an adequate method for cool dwarfs. The chemical equilibrium solver has been extended to include polyatomic and charged molecules, extending Korg's regime of applicability to M stars. We also discuss a common oversight regarding the synthesis of hydrogen lines in the infrared, and show that Korg's Brackett line profiles are a much closer match to observations than others available. Documentation, installation instructions, and tutorials are available at https://github.com/ajwheeler/Korg.jl.

Reviewing the empirical and theoretical parameter relationships between various parameters is a good way to understand more about contact binary systems. In this investigation, two-dimensional (2D) relationships for P-M_V(system), P-L_1,2, M_1,2-L_1,2, and q-L_ratio were revisited. The sample used is related to 118 contact binary systems with an orbital period shorter than 0.6 days whose absolute parameters were estimated based on the Gaia Data Release 3 (DR3) parallax. We reviewed previous studies on 2D relationships and updated six parameter relationships. Therefore, Markov chain Monte Carlo (MCMC) and Machine Learning (ML) methods were used, and the outcomes were compared. We selected 22 contact binary systems from eight previous studies for comparison, which had light curve solutions using spectroscopic data. The results show that the systems are in good agreement with the results of this study.

ESA's Hera mission aims to visit binary asteroid Didymos in late 2026, investigating its physical characteristics and the result of NASA's impact by the DART spacecraft in more detail. Two CubeSats on-board Hera plan to perform a ballistic landing on the secondary of the system, called Dimorphos. For these types of landings the translational state during descent is not controlled, reducing the spacecrafts complexity but also increasing its sensitivity to deployment maneuver errors and dynamical uncertainties. This paper introduces a novel methodology to analyse the effect of these uncertainties on the dynamics of the lander and design a trajectory that is robust against them. This methodology consists of propagating the uncertain state of the lander using the non-intrusive Chebyshev interpolation (NCI) technique, which approximates the uncertain dynamics using a polynomial expansion, and analysing the results using the pseudo-diffusion indicator, derived from the coefficients of the polynomial expansion, which quantifies the rate of growth of the set of possible states of the spacecraft over time. This indicator is used here to constrain the impact velocity and angle to values which allow for successful settling on the surface. This information is then used to optimize the landing trajectory by applying the NCI technique inside the transcription of the problem. The resulting trajectory increases the robustness of the trajectory compared to a conventional method, improving the landing success by 20 percent and significantly reducing the landing footprint.

In the recent years, primordial black holes (PBHs) have emerged as one of the most interesting and hotly debated topics in cosmology. Among other possibilities, PBHs could explain both some of the signals from binary black hole mergers observed in gravitational wave detectors and an important component of the dark matter in the Universe. Significant progress has been achieved both on the theory side and from the point of view of observations, including new models and more accurate calculations of PBH formation, evolution, clustering, merger rates, as well as new astrophysical and cosmological probes. In this work, we review, analyse and combine the latest developments in order to perform end-to-end calculations of the various gravitational wave signatures of PBHs. Different ways to distinguish PBHs from stellar black holes are emphasized. Finally, we discuss their detectability with LISA, the first planned gravitational-wave observatory in space.

We investigate the formation of dense stellar clumps in a suite of high-resolution cosmological zoom-in simulations of a massive, star forming galaxy at $z \sim 2$ under the presence of strong quasar winds. Our simulations include multi-phase ISM physics from the Feedback In Realistic Environments (FIRE) project and a novel implementation of hyper-refined accretion disk winds. We show that powerful quasar winds can have a global negative impact on galaxy growth while in the strongest cases triggering the formation of an off-center clump with stellar mass ${\rm M}_{\star}\sim 10^{7}\,{\rm M}_{\odot}$, effective radius ${\rm R}_{\rm 1/2\,\rm Clump}\sim 20\,{\rm pc}$, and surface density $\Sigma_{\star} \sim 10^{4}\,{\rm M}_{\odot}\,{\rm pc}^{-2}$. The clump progenitor gas cloud is originally not star-forming, but strong ram pressure gradients driven by the quasar winds (orders of magnitude stronger than experienced in the absence of winds) lead to rapid compression and subsequent conversion of gas into stars at densities much higher than the average density of star-forming gas. The AGN-triggered star-forming clump reaches ${\rm SFR} \sim 50\,{\rm M}_{\odot}\,{\rm yr}^{-1}$ and $\Sigma_{\rm SFR} \sim 10^{4}\,{\rm M}_{\odot}\,{\rm yr}^{-1}\,{\rm kpc}^{-2}$, converting most of the progenitor gas cloud into stars in $\sim$2\,Myr, significantly faster than its initial free-fall time and with stellar feedback unable to stop star formation. In contrast, the same gas cloud in the absence of quasar winds forms stars over a much longer period of time ($\sim$35\,Myr), at lower densities, and losing spatial coherency. The presence of young, ultra-dense, gravitationally bound stellar clumps in recently quenched galaxies could thus indicate local positive feedback acting alongside the strong negative impact of powerful quasar winds, providing a plausible formation scenario for globular clusters.

We present the DECam Ecliptic Exploration Project (DEEP) survey strategy including observing cadence for orbit determination, exposure times, field pointings and filter choices. The overall goal of the survey is to discover and characterize the orbits of a few thousand Trans-Neptunian Objects (TNOs) using the Dark Energy Camera (DECam) on the Cerro Tololo Inter-American Observatory (CTIO) Blanco 4 meter telescope. The experiment is designed to collect a very deep series of exposures totaling a few hours on sky for each of several 2.7 square degree DECam fields-of-view to achieve a magnitude of about 26.2 using a wide VR filter which encompasses both the V and R bandpasses. In the first year, several nights were combined to achieve a sky area of about 34 square degrees. In subsequent years, the fields have been re-visited to allow TNOs to be tracked for orbit determination. When complete, DEEP will be the largest survey of the outer solar system ever undertaken in terms of newly discovered object numbers, and the most prolific at producing multi-year orbital information for the population of minor planets beyond Neptune at 30 au.

Measuring the mass-radius relation of individual white dwarfs is an empirically challenging task that has been performed for only a few dozen stars. We measure the white dwarf mass-radius relation using gravitational redshifts and radii of 137 white dwarfs in wide binaries with main sequence companions. We obtain the space velocities to these systems using the main sequence companion, and subtract these Doppler redshifts from the white dwarfs' apparent motions, isolating their gravitational redshifts. We use Gaia data to calculate the surface temperatures and radii of these white dwarfs, thereby deriving an empirical gravitational redshift-radius relation. This work demonstrates the utility of low-resolution Galactic surveys to measure the white dwarf equation of state. Our results are consistent with theoretical models, and represent the largest sample of individual white dwarf gravitational redshift measurements to date.

Common-envelope evolution - where a star is engulfed by a companion - is a critical but poorly understood step in, e.g., the formation pathways for gravitational-wave sources. However, it has been extremely challenging to identify observable signatures of such systems. We show that for systems involving a neutron star, the hypothesized super-Eddington accretion onto the neutron star produces MeV-range, months-long neutrino signals within reach of present and planned detectors. While there are substantial uncertainties on the rate of such events (0.01-1/century in the Milky Way) and the neutrino luminosity (which may be less than the accretion power), this signal can only be found if dedicated new searches are developed. If detected, the neutrino signal would lead to significant new insights into the astrophysics of common-envelope evolution.

4-Dimensional Einstein-Gauss-Bonnet (4DEGB) gravity has garnered significant attention in the last few years as a phenomenological competitor to general relativity. We consider the theoretical and observational implications of this theory in both the early and late universe, (re-)deriving background and perturbation equations and constraining its characteristic parameters with data from cosmological probes. Our investigation surpasses the scope of previous studies by incorporating non-flat spatial sections. We explore consequences of 4DEGB on the sound and particle horizons in the very early universe, and demonstrate that 4DEGB can provide an independent solution to the horizon problem for some values of its characteristic parameter $\alpha$. Finally, we constrain an unexplored regime of this theory in the limit of small coupling $\alpha$ (empirically supported in the post-Big Bang Nucleosynthesis era by prior constraints). This version of 4DEGB includes a geometric term that resembles dark radiation at the background level, but whose influence on the perturbed equations is qualitatively distinct from that of standard forms of dark radiation. In this limit, only one beyond-$\Lambda$CDM degree of freedom persists, which we denote as $\tilde{\alpha}_C$. Our analysis yields the estimate $\tilde{\alpha}_C = (-9 \pm 6) \times 10^{-6}$ thereby providing a new constraint of a previously untested sector of 4DEGB.

The success of galactic archaeology and the reconstruction of the formation history of our galaxy critically relies on precise ages for large populations of stars. For evolved stars in the red clump and red giant branch, the carbon to nitrogen ratio ([C/N]) has recently been identified as a powerful diagnostic of mass and age that can be applied to stellar samples from spectroscopic surveys such as SDSS/APOGEE. Here, we show that at least 10\% of red clump stars and %$\approx 10\%$ of red giant branch stars deviate from the standard relationship between [C/N] and mass. {We use the APOGEE-\kepler\ (APOKASC) overlap sample to show that binary interactions are %the majority contributors to these responsible for the majority of these outliers and that stars with %any indicators of current or previous binarity should be excluded from galactic archaeology analyses that rely on [C/N] abundances to infer stellar masses. We also show that the %standard DR14 APOGEE analysis overestimates the surface gravities for even moderately rotating giants (vsini$>2$ km/s)}

Strong magnetic fields are expected to significantly modify the pulsation frequencies of waves propagating in the cores of red giants or in the radiative envelopes of intermediate- and high-mass main-sequence stars. We calculate the g-mode frequencies of stars with magnetic dipole fields which are aligned with their rotational axes, treating both the Lorentz and Coriolis forces non-perturbatively. We provide a compact asymptotic formula for the g-mode period spacing, and universally find that strong magnetism decreases this period spacing substantially more than is predicted by perturbation theory. These results are validated with explicit numerical mode calculations for realistic stellar models. The approach we present is highly versatile: once the eigenvalues $\lambda$ of a certain differential operator are precomputed as a function of the magnetogravity and rotational frequencies (in units of the mode frequency), the non-perturbative impact of the Coriolis and Lorentz forces is understood under a broad domain of validity, and is readily incorporated into asteroseismic modeling.

In this paper, we present updated estimates of the velocity of the neutron star (NS) in the supernova remnant (SNR) Cassiopeia A using over two decades of Chandra observations. We use two methods: 1.) recording NS positions from dozens of Chandra observations, including the astrometric uncertainty estimates on the data points but not correcting the astrometry of the observations, and 2.) correcting the astrometry of the 13 Chandra observations that have a sufficient number of point sources with identified Gaia counterparts. For method #1, we find velocity of 280 $\pm$ 123 km s$^{-1}$, with an angle of 87 $\pm$ 22 degrees south of east. For method #2, we find a velocity of 445 $\pm$ 90 km s$^{-1}$ at an angle of 68 $\pm$ 12 degrees south of east. Both of these results match with the explosion-center-estimated velocity of $\sim$350 km s$^{-1}$ and the previous 10 year baseline proper motion measurement of 570 $\pm$ 260 km s$^{-1}$, but our use of additional data over a longer baseline has led to a smaller uncertainty by a factor of 2$\unicode{x2013}$3. Our estimates rule out velocities $\gtrsim$600 km s$^{-1}$ and better match with simulations of Cassiopeia A that include NS kick mechanisms.

SDSS J1004+4112 is a well studied gravitational lens with a recently measured time delay between its first and fourth arriving quasar images. Using this new constraint, we present updated free-form lens reconstructions using the lens inversion method {\tt GRALE}, which only uses multiple image and time delay data as inputs. In addition, we obtain hybrid lens reconstructions by including a model of the brightest cluster galaxy (BCG) as a Sersic lens. For both reconstructions, we use two sets of images as input: one with all identified images, and the other a revised set leaving out images that have been potentially misidentified. We also develop a source position optimization MCMC routine, performed on completed {\tt GRALE} runs, that allows each model to better match observed image positions and time delays. All the reconstructions produce similar mass distributions, with the hybrid models finding a steeper profile in the center. Similarly, all the mass distributions are fit by the Navarro-Frenk-White (NFW) profile, finding results consistent with previous parametric reconstructions and those derived from Chandra X-ray observations. We identify a $\sim 5 \times 10^{11} M_{\odot}$ substructure apparently unaffiliated with any cluster member galaxy and present in all our models, and study its reality. Using our free-form and hybrid models we predict a central quasar image time delay of $\sim 2980 \pm 270$ and $\sim 3280 \pm 215$ days, respectively. A potential future measurement of this time delay will, while being an observational challenge, further constrain the steepness of the central density profile.

FR0 galaxies constitute the most abundant jet population in the local Universe. With their compact jet structure, they are broadband photon emitters and have been proposed as multi-messenger sources. Recently, these sources have been detected for the first time in $\gamma$ rays. Using a revised FR0 catalog, we confirm that the FR0 population as a whole are $\gamma$-ray emitters, and we also identify two significant sources. For the first time, we find a correlation between the 5 GHz core radio luminosity and $\gamma$-ray luminosity in the 1 - 800 GeV band, having a 4.5$\sigma$ statistical significance. This is clear evidence that the jet emission mechanism is similar in nature for FR0s and the well-studied canonical FR (FRI and FRII) radio galaxies. Furthermore, we perform broadband SED modeling for the significantly detected sources as well as the subthreshold source population using a one-zone SSC model. Within the maximum jet power budget, our modeling shows that the detected gamma rays from the jet can be explained as inverse Compton photons. To explain the multi-wavelength observations for these galaxies, the modeling results stipulate a low bulk Lorentz factor and a jet composition far from equipartition, with the particle energy density dominating over the magnetic field energy density.

We review and update constraints on the Early Dark Energy (EDE) model from cosmological data sets, in particular Planck PR3 and PR4 cosmic microwave background (CMB) data and large-scale structure (LSS) data sets including galaxy clustering and weak lensing data from the Dark Energy Survey, Subaru Hyper Suprime-Cam, and KiDS+VIKING-450, as well as BOSS/eBOSS galaxy clustering and Lyman-$\alpha$ forest data. We detail the fit to CMB data, and perform the first analyses of EDE using the CAMSPEC and Hillipop likelihoods for Planck CMB data, rather than Plik, both of which yield a tighter upper bound on the allowed EDE fraction than that found with Plik. We then supplement CMB data with large-scale structure data in a series of new analyses. All these analyses are concordant in their Bayesian preference for $\Lambda$CDM over EDE, as indicated by marginalized posterior distributions. We perform a series of tests of the impact of priors in these results, and compare with frequentist analyses based on the profile likelihood, finding qualitative agreement with the Bayesian results. All these tests suggest prior volume effects are not a determining factor in analyses of EDE. This work provides both a review of existing constraints and several new analyses.

We present highly sensitive measurements taken with MeerKAT at 1280 MHz as well as archival GBT, MWA and VLA images at 333, 88 and 74 MHz. We report the detection of synchrotron radio emission from the infrared dark cloud (IRDC) associated with the halo of the Sgr B complex on a scale of ~60 pc. A strong spatial correlation between low-frequency radio continuum emission and dense molecular gas, combined with spectral index measurements, indicates enhanced synchrotron emission by cosmic-ray electrons. Correlation of the FeI 6.4 keV Kalpha line and synchrotron emission provides compelling evidence that the low energy cosmic-ray electrons are responsible for producing the Kalpha line emission. The observed synchrotron emission within the halo of the Sgr B cloud complex has mean spectral index alpha -1+/-1 gives the magnetic field strength ~100 muG for cloud densities nH = 10^4-10^5 cm-3 and estimate cosmic-ray ionization rates between 10^-13 and 10^-14 s^-1. Furthermore, the energy spectrum of primary cosmic-ray electrons is constrained to be E^-3 +/-1 for typical energies of few hundred MeV. The extrapolation of this spectrum to higher energies is consistent with X-ray and gamma-ray emission detected from this cloud. These measurements have important implications on the role that high cosmic-ray electron fluxes at the Galactic center play in production of radio synchrotron emission, the FeI Kalpha line emission at 6.4 keV and ~GeV gamma-ray emission throughout the central molecular zone (CMZ).

Modern cosmological hydrodynamical galaxy simulations provide tens of thousands of reasonably realistic synthetic galaxies across cosmic time. However, quantitatively assessing the level of realism of simulated universes in comparison to the real one is difficult. In this paper of the ERGO-ML series (Extracting Reality from Galaxy Observables with Machine Learning), we utilize contrastive learning to directly compare a large sample of simulated and observed galaxies based on their stellar-light images. This eliminates the need to specify summary statistics and allows to exploit the whole information content of the observations. We produce survey-realistic galaxy mock datasets resembling real Hyper Suprime-Cam (HSC) observations using the cosmological simulations TNG50 and TNG100. Our focus is on galaxies with stellar masses between $10^9$ and $10^{12} M_\odot$ at $z=0.1-0.4$. This allows us to evaluate the realism of the simulated TNG galaxies in comparison to actual HSC observations. We apply the self-supervised contrastive learning method NNCLR to the images from both simulated and observed datasets (g, r, i - bands). This results in a 256-dimensional representation space, encoding all relevant observable galaxy properties. Firstly, this allows us to identify simulated galaxies that closely resemble real ones by seeking similar images in this multi-dimensional space. Even more powerful, we quantify the alignment between the representations of these two image sets, finding that the majority ($\gtrsim 70$ per cent) of the TNG galaxies align well with observed HSC images. However, a subset of simulated galaxies with larger sizes, steeper Sersic profiles, smaller Sersic ellipticities, and larger asymmetries appears unrealistic. We also demonstrate the utility of our derived image representations by inferring properties of real HSC galaxies using simulated TNG galaxies as the ground truth.

Reconstructing astrophysical and cosmological fields from observations is challenging. It requires accounting for non-linear transformations, mixing of spatial structure, and noise. In contrast, forward simulators that map fields to observations are readily available for many applications. We present a versatile Bayesian field reconstruction algorithm rooted in simulation-based inference and enhanced by autoregressive modeling. The proposed technique is applicable to generic (non-differentiable) forward simulators and allows sampling from the posterior for the underlying field. We show first promising results on a proof-of-concept application: the recovery of cosmological initial conditions from late-time density fields.

A supermassive black hole (SMBH) surrounded by a dense, nuclear star cluster resides at the center of many galaxies. In this dense environment, high-velocity collisions frequently occur between stars. About $10 \%$ of the stars within the Milky Way's nuclear star cluster collide with other stars before evolving off the main-sequence. Collisions preferentially affect tightly-bound stars, which orbit most quickly and pass through regions of the highest stellar density. Over time, collisions therefore shape the bulk properties of the nuclear star cluster. We examine the effect of collisions on the cluster's stellar density profile. We show that collisions produce a turning point in the density profile which can be determined analytically. Varying the initial density profile and collision model, we characterize the evolution of the stellar density profile over $10$ Gyr. We find that old, initially cuspy populations exhibit a break around $0.1$ pc in their density profile, while shallow density profiles retain their initial shape outside of $0.01$ pc. The initial density profile is always preserved outside of a few tenths of parsec irrespective of initial conditions. Lastly, we comment on the implications of collisions for the luminosity and color of stars in the collisionly-shaped inner cluster.

The standard paradigm for the formation of the Universe suggests that large structures are formed from hierarchical clustering by the continuous accretion of less massive galaxy systems through filaments. In this context, filamentary structures play an important role in the properties and evolution of galaxies by connecting high-density regions, such as nodes, and being surrounded by low-density regions, such as cosmic voids. The availability of the filament and point critic catalogues extracted by \textsc{DisPerSE} from the \textsc{Illustris} TNG300-1 hydrodynamic simulation allows a detailed analysis of these structures. The halo occupation distribution (HOD) is a powerful tool for linking galaxies and dark matter halos, allowing constrained models of galaxy formation and evolution. In this work we combine the advantage of halo occupancy with information from the filament network to analyse the HOD in filaments and nodes. In our study, we distinguish the inner regions of cosmic filaments and nodes from their surroundings. The results show that the filamentary structures have a similar trend to the total galaxy sample indicating that, although the filaments span a wide range of densities, they may represent regions of average density. In the case of the nodes sample, an excess of faint and blue galaxies is found for the low-mass nodes suggesting that these structures are not virialised and that galaxies may be continuously falling through the filaments. Instead, the higher-mass halos could be in a more advanced stage of evolution showing features of virialised structures.

We collect near-infrared spectra ($\sim0.75-2.55\ \mu m$) of four Jovian irregular satellites and visible spectra ($\sim0.32-1.00\ \mu m$) of two Jovian irregular satellites, two Uranian irregular satellites, and four Neptune Trojans. We find close similarities between observed Jovian irregular satellites and previously characterized Jovian Trojans. However, irregular satellites' unique collisional histories complicate comparisons to other groups. Laboratory study of CM and CI chondrites show that grain size and regolith packing conditions strongly affect spectra of dark, carbonaceous materials. We hypothesize that different activity histories of these objects, which may have originally contained volatile ices that subsequently sublimated, could cause differences in regolith grain-size or packing properties and therefore drive spectral variation. The Uranian satellites Sycorax and Caliban appear similar to TNOs. However, we detect a feature near 0.7 $\mu m$ on Sycorax, suggesting the presence of hydrated materials. While the sample of Neptune Trojans have more neutral spectra than the Uranian satellites we observe, they remain consistent with the broad color distribution of the Kuiper belt. We detect a possible feature near 0.65-0.70 $\mu m$ on Neptune Trojan 2006 RJ103, suggesting that hydrated material may also be present in this population. Characterizing hydrated materials in the outer solar system may provide critical context regarding the origins of hydrated CI and CM chondrite meteorites. We discuss how the hydration state(s) of the irregular satellites constrains the thermal histories of the interiors of their parent bodies, which may have formed among the primordial Kuiper belt.

Core-collapse supernovae (SNs) are one of the most powerful cosmic sources of neutrinos, with energies of several MeV. The emission of neutrinos and antineutrinos of all flavors carries away the gravitational binding energy of the compact remnant and drives its evolution from the hot initial to the cold final states. Detecting these neutrinos from Earth and analyzing the emitted signals present a unique opportunity to explore the neutrino mass ordering problem. This research outlines the detection of neutrinos from SNs and their relevance in understanding the neutrino mass ordering. The focus is on developing a model-independent analysis strategy, achieved by comparing distinct detection channels in large underground detectors. The objective is to identify potential indicators of mass ordering within the neutrino sector. Additionally, a thorough statistical analysis is performed on the anticipated neutrino signals for both mass orderings. Despite uncertainties in supernova explosion parameters, an exploration of the parameter space reveals an extensive array of models with significant sensitivity to differentiate between mass orderings. The assessment of various observables and their combinations underscores the potential of forthcoming supernova observations in addressing the neutrino mass ordering problem.

The recent serendipitous discovery of a new population of short duration X-ray transients, thought to be associated with collisions of compact objects or stellar explosions in distant galaxies, has motivated efforts to build up statistical samples by mining X-ray telescope archives. Most searches to date however, do not fully exploit recent developments in the signal and imaging processing research domains to optimise searches for short X-ray flashes. This paper addresses this issue by presenting a new source detection pipeline, STATiX (Space and Time Algorithm for Transients in X-rays), which directly operates on 3-dimensional X-ray data cubes consisting of two spatial and one temporal dimension. The algorithm leverages wavelet transforms and the principles of sparsity to denoise X-ray observations and then detect source candidates on the denoised data cubes. The light curves of the detected sources are then characterised using the Bayesian blocks algorithm to identify flaring periods. We describe the implementation of STATiX in the case of XMM-Newton data, present extensive validation and performance verification tests based on simulations and also apply the pipeline to a small subset of seven XMM-Newton observations, which are known to contain transients sources. In addition to known flares in the selected fields we report a previously unknown short duration transient found by our algorithm that is likely associated with a flaring Galactic star. This discovery demonstrates the potential of applying STATiX to the full XMM-Newton archive.

The existence of a plateau in the short-duration tail of the observed distribution of cosmological Long-soft Gamma Ray Bursts (LGRBs) has been argued as the first direct evidence of Collapsars. A similar plateau in the short-duration tail of the observed duration distribution of Short-hard Gamma Ray Bursts (SGRBs) has been suggested as evidence of compact binary mergers. We present an equally plausible alternative interpretation for this evidence, which is purely statistical. Specifically, we show that the observed plateau in the short-duration tail of the duration distribution of LGRBs can naturally occur in the statistical distributions of strictly-positive physical quantities, exacerbated by the effects of mixing with the duration distribution of SGRBs, observational selection effects and data aggregation (e.g., binning) methodologies. The observed plateau in the short-duration tail of the observed distributions of SGRBs can similarly result from a combination of sample incompleteness and inhomogeneous binning of data.

The near infrared background (NIRB) is the collective light from unresolved sources observed in the band 1-10 $\mu$m. The measured NIRB angular power spectrum on angular scales $\theta \gtrsim 1$ arcmin exceeds by roughly two order of magnitudes predictions from known galaxy populations. The nature of the sources producing these fluctuations is still unknown. Here we test primordial black holes (PBHs) as sources of the NIRB excess. Considering PBHs as a cold dark matter (DM) component, we model the emission of gas accreting onto PBHs in a cosmological framework. We account for both accretion in the intergalactic medium (IGM) and in DM haloes. We self consistently derive the IGM temperature evolution, considering ionization and heating due to X-ray emission from PBHs. Besides $\Lambda$CDM, we consider a model that accounts for the modification of the linear matter power spectrum due to the presence of PBHs; we also explore two PBH mass distributions, i.e. a $\delta$-function and a lognormal distribution. For each model, we compute the mean intensity and the angular power spectrum of the NIRB produced by PBHs with mass 1-$10^3~\mathrm{M}_{\odot}$. In the limiting case in which the entirety of DM is made of PBHs, the PBH emission contributes $<1$ per cent to the observed NIRB fluctuations. This value decreases to $<0.1$ per cent if current constraints on the abundance of PBHs are taken into account. We conclude that PBHs are ruled out as substantial contributors to the NIRB.

Screening mechanisms allow light scalar fields to dynamically avoid the constraints that come from our lack of observation of a long-range fifth force. Galactic scale tests are of particular interest when the light scalar is introduced to explain the dark matter or dark energy that dominates our cosmology. To date, much of the literature that has studied screening in galaxies has described screening using simplifying approximations. In this work, we calculate numerical solutions for scalar fields with screening mechanisms in galactic contexts, and use these to derive new, precise conditions governing where fifth forces are screened. We show that the commonly used binary screened/unscreened threshold can predict a fifth force signal in situations where a fuller treatment does not, leading us to conclude that existing constraints might be significantly overestimated. We show that various other approximations of the screening radius provide a more accurate proxy to screening, although they fail to exactly reproduce the true screening surface in certain regions of parameter space. As a demonstration of our scheme, we apply it to an idealised Milky Way and thus identify the region of parameter space in which the solar system is screened.

We show, for the first time, radio measurements of the depth of shower maximum ($X_\text{max}$) of air showers induced by cosmic rays that are compared to measurements of the established fluorescence method at the same location. Using measurements at the Pierre Auger Observatory we show full compatibility between our radio and the previously published fluorescence data set, and between a subset of air showers observed simultaneously with both radio and fluorescence techniques, a measurement setup unique to the Pierre Auger Observatory. Furthermore, we show radio $X_\text{max}$ resolution as a function of energy and demonstrate the ability to make competitive high-resolution $X_\text{max}$ measurements with even a sparse radio array. With this, we show that the radio technique is capable of cosmic-ray mass composition studies, both at Auger and at other experiments.

The Auger Engineering Radio Array (AERA), part of the Pierre Auger Observatory, is currently the largest array of radio antenna stations deployed for the detection of cosmic rays, spanning an area of $17$ km$^2$ with 153 radio stations. It detects the radio emission of extensive air showers produced by cosmic rays in the $30-80$ MHz band. Here, we report the AERA measurements of the depth of the shower maximum ($X_\text{max}$), a probe for mass composition, at cosmic-ray energies between $10^{17.5}$ to $10^{18.8}$ eV, which show agreement with earlier measurements with the fluorescence technique at the Pierre Auger Observatory. We show advancements in the method for radio $X_\text{max}$ reconstruction by comparison to dedicated sets of CORSIKA/CoREAS air-shower simulations, including steps of reconstruction-bias identification and correction, which is of particular importance for irregular or sparse radio arrays. Using the largest set of radio air-shower measurements to date, we show the radio $X_\text{max}$ resolution as a function of energy, reaching a resolution better than $15$ g cm$^{-2}$ at the highest energies, demonstrating that radio $X_\text{max}$ measurements are competitive with the established high-precision fluorescence technique. In addition, we developed a procedure for performing an extensive data-driven study of systematic uncertainties, including the effects of acceptance bias, reconstruction bias, and the investigation of possible residual biases. These results have been cross-checked with air showers measured independently with both the radio and fluorescence techniques, a setup unique to the Pierre Auger Observatory.

Big Bang Nucleosynthesis (BBN) is a strong probe for constraining new physics including gravitation. $f(R)$ gravity theory is an interesting alternative to general relativity which introduces additional degrees of freedom known as scalarons. In this work we demonstrate the existence of black hole solutions in $f(R)$ gravity and develop a relation between scalaron mass and black hole mass. We have used observed bound on the freezeout temperature to constrain scalaron mass range by modifying the cosmic expansion rate at the BBN epoch. The mass range of primordial black holes (PBHs) which are astrophysical dark matter candidates is deduced. The range of scalaron mass which does not spoil the BBN era is found to be $10^{-16}-10^4 \text{ eV}$. The scalaron mass window $10^{-16}-10^{-14}\text{ eV}$ is consistent with the $f(R)$ gravity PPN parameter derived from solar system experiments. The PBH mass range is obtained as $10^6-10^{-14}\text{ }M_{\odot}$. Scalarons constrained by BBN are also eligible to accommodate axion like dark matter particles. Estimation of deuterium (D) fraction and relative D+$^3$He abundance in the $f(R)$ gravity scenario shows that the BBN history mimics that of general relativity. While the PBH mass range is eligible for non-baryonic dark matter, the BBN bounded scalarons provide with an independent strong field test of $f(R)$ gravity.

The era of precision cosmology allows us to test the composition of the dark matter. Mixed ultralight or fuzzy dark matter (FDM) is a cosmological model with dark matter composed of a combination of particles of mass $m\leq 10^{-20}$ eV, with an astrophysical de Broglie wavelength, and particles with a negligible wavelength sharing the properties of cold dark matter (CDM). In this work, we simulate cosmological volumes with a dark matter wave function for the ultralight component coupled gravitationally to CDM particles. We investigate the impact of a mixture of CDM and FDM in various proportions $(0\%,\;1\%,\;10\%,\;50\%,\;100\%)$ and for ultralight particle masses ranging over five orders of magnitude $(2.5\times 10^{-25}\;\mathrm{eV}-2.5\times 10^{-21}\;\mathrm{eV})$. To track the evolution of density perturbations in the non-linear regime, we adapt the simulation code AxioNyx to solve the CDM dynamics coupled to a FDM wave function obeying the Schr\"odinger-Poisson equations. We obtain the non-linear power spectrum and study the impact of the wave effects on the growth of structure on different scales. We confirm that the steady-state solution of the Schr\"odinger-Poisson system holds at the center of halos in the presence of a CDM component when it composes $50\%$ or less of the dark matter but find no stable density core when the FDM accounts for $10\%$ or less of the dark matter. We implement a modified friends-of-friends halo finder and find good agreement between the observed halo abundance and the predictions from the adapted halo model axionHMCode.

For diffraction-limited optical systems an accurate physical optics model is necessary to properly evaluate instrument performance. Astronomical observatories outfitted with coronagraphs for direct exoplanet imaging require physical optics models to simulate the effects of misalignment and diffraction. Accurate knowledge of the observatory's PSF is integral for the design of high-contrast imaging instruments and simulation of astrophysical observations. The state of the art is to model the misalignment, ray aberration, and diffraction across multiple software packages, which complicates the design process. Gaussian Beamlet Decomposition (GBD) is a ray-based method of diffraction calculation that has been widely implemented in commercial optical design software. By performing the coherent calculation with data from the ray model of the observatory, the ray aberration errors can be fed directly into the physical optics model of the coronagraph, enabling a more integrated model of the observatory. We develop a formal algorithm for the transfer-matrix method of GBD, and evaluate it against analytical results and a traditional physical optics model to assess the suitability of GBD for high-contrast imaging simulations. Our GBD simulations of the observatory PSF, when compared to the analytical Airy function, have a sum-normalized RMS difference of ~10^-6. These fields are then propagated through a Fraunhofer model of a exoplanet imaging coronagraph where the mean residual numerical contrast is 4x10^-11, with a maximum near the inner working angle at 5x10^-9. These results show considerable promise for the future development of GBD as a viable propagation technique in high-contrast imaging. We developed this algorithm in an open-source software package and outlined a path for its continued development to increase the fidelity and flexibility of diffraction simulations using GBD.

We study the evolution of heavy stars ($M\ge40{\rm M}_\odot$) undergoing pair-instability in the presence of annihilating dark matter. Focusing on the scenario where the dark matter is in capture-annihilation equilibrium, we model the profile of energy injections in the local thermal equilibrium approximation. We find that significant changes to masses of astrophysical black holes formed by (pulsational) pair-instability supernovae can occur when the ambient dark matter density $ \rho_{\rm DM} \gtrsim10^9 \rm \, GeV \, cm^{-3}$. There are two distinct outcomes, depending on the dark matter mass. For masses $m_{\rm DM}\gtrsim1$ GeV the DM is primarily confined to the core. The annihilation increases the lifetime of core helium burning, resulting in more oxygen being formed, fueling a more violent explosion during the pair-instability-induced contraction. This drives stronger pulsations, leading to lighter black holes being formed than predicted by the standard model. For masses $m_{\rm DM}\lesssim0.5$ GeV there is significant dark matter in the envelope, leading to a phase where the star is supported by the energy from the annihilation. This reduces the core temperature and density, allowing the star to evade the pair-instability allowing heavier black holes to be formed. We find a mass gap for all models studied.

Solar-like oscillations have been detected in thousands of stars thanks to modern space missions. These oscillations have been used to measure stellar masses and ages, which have been widely applied in Galactic archaeology. One of the pillars of such applications is the $\nu_{\max}$ scaling relation: the frequency of maximum power $\nu_{\max}$, assumed to be proportional to the acoustic cut-off frequency, $\nu_{\rm ac}$, scales with effective temperature and surface gravity. However, the theoretical basis of the $\nu_{\max}$ scaling relation is uncertain, and there is an ongoing debate about whether it can be applied to metal-poor stars. We investigate the metallicity dependence of the $\nu_{\max}$ scaling relation by carrying out 3D near-surface convection simulations for solar-type stars with [Fe/H] between -3 and 0.5 dex. Firstly, we found a negative correlation between $\nu_{\rm ac}$ and metallicity from the 3D models. This is in tension with the positive correlation identified by studies using 1D models. Secondly, we estimated theoretical $\nu_{\max}$ values using velocity amplitudes determined from first principles, by quantifying the mode excitation and damping rates with methods validated in our previous works. We found that at solar effective temperature and surface gravity, $\nu_{\max}$ does not show correlation with metallicity. This study opens an exciting prospect of testing the asteroseismic scaling relations against realistic 3D hydrodynamical stellar models.

Spiral structures are important drivers of the secular evolution of disc galaxies, however, the origin of spiral arms and their effects on the development of galaxies remain mysterious. In this work, we present two three-armed spiral galaxies at z~0.3 in the Middle Age Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey. Taking advantage of the high spatial resolution (~0.6'') of the Multi-Unit Spectroscopic Unit (MUSE), we investigate the two-dimensional distributions of different spectral parameters: Halpha, gas-phase metallicity, and D4000. We notice significant offsets in Halpha (~0.2 dex) as well as gas-phase metallicities (~0.05 dex) among the spiral arms, downstream and upstream of MAGPI1202197197 (SG1202). This observational signature suggests the spiral structure in SG1202 is consistent with arising from density wave theory. No azimuthal variation in Halpha or gas-phase metallicities is observed in MAGPI1204198199 (SG1204), which can be attributed to the tighter spiral arms in SG1204 than SG1202, coming with stronger mixing effects in the disc. The absence of azimuthal D4000 variation in both galaxies suggests the stars at different ages are well-mixed between the spiral arms and distributed around the disc regions. The different azimuthal distributions in Halpha and D4000 highlight the importance of time scales traced by various spectral parameters when studying 2D distributions in spiral galaxies. This work demonstrates the feasibility of constraining spiral structures by tracing interstellar medium (ISM) and stellar population at z~0.3, with a plan to expand the study to the full MAGPI survey.

Time dependent photoionization modeling of warm absorber outflows in active galactic nuclei can play an important role in understanding the interaction between warm absorbers and the central black hole. The warm absorber may be out of the equilibrium state because of the variable nature of the central continuum. In this paper, with the help of time dependent photoionization modeling, we study how the warm absorber gas changes with time and how it reacts to changing radiation fields. Incorporating a flaring incident light curve, we investigate the behavior of warm absorbers using a photoionization code that simultaneously and consistently solves the time dependent equations of level population, heating and cooling, and radiative transfer. We simulate the physical processes in the gas clouds, such as ionization, recombination, heating, cooling, and the transfer of ionizing radiation through the cloud. We show that time dependent radiative transfer is important and that calculations which omit this effect quantitatively and systematically underestimate the absorption. Such models provide crucial insights into the characteristics of warm absorbers and can constrain their density and spatial distribution.

We present a hierarchical Bayesian pipeline, BP3M, that measures positions, parallaxes, and proper motions (PMs) for cross-matched sources between Hubble~Space~Telescope (HST) images and Gaia -- even for sparse fields ($N_*<10$ per image) -- expanding from the recent GaiaHub tool. This technique uses Gaia-measured astrometry as priors to predict the locations of sources in HST images, and is therefore able to put the HST images onto a global reference frame without the use of background galaxies/QSOs. Testing our publicly-available code in the Fornax and Draco dSphs, we measure accurate PMs that are a median of 8-13 times more precise than Gaia DR3 alone for $20.5<G<21~\mathrm{mag}$. We are able to explore the effect of observation strategies on BP3M astrometry using synthetic data, finding an optimal strategy to improve parallax and position precision at no cost to the PM uncertainty. Using 1619 HST images in the sparse COSMOS field (median 9 Gaia sources per HST image), we measure BP3M PMs for 2640 unique sources in the $16<G<21.5~\mathrm{mag}$ range, 25% of which have no Gaia PMs; the median BP3M PM uncertainty for $20.25<G<20.75~\mathrm{mag}$ sources is $0.44~$mas/yr compared to $1.03~$mas/yr from Gaia, while the median BP3M PM uncertainty for sources without Gaia-measured PMs ($20.75<G<21.5~\mathrm{mag}$) is $1.16~$mas/yr. The statistics that underpin the BP3M pipeline are a generalized way of combining position measurements from different images, epochs, and telescopes, which allows information to be shared between surveys and archives to achieve higher astrometric precision than that from each catalog alone.

One of the fundamental questions in astronomy is how planetary systems form and evolve. Measuring the planetary occurrence and architecture as a function of time directly addresses this question. In the fourth paper of the Planets Across Space and Time (PAST) series, we investigate the occurrence and architecture of Kepler planetary systems as a function of kinematic age by using the LAMOST-Gaia-Kepler sample. To isolate the age effect, other stellar properties (e.g., metallicity) have been controlled. We find the following results. (1) The fraction of stars with Kepler-like planets ($F_{\text{Kep}}$) is about 50% for all stars; no significant trend is found between $F_{\text{Kep}}$ and age. (2) The average planet multiplicity ($\bar{N}_p$) exhibits a decreasing trend (~2$\sigma$ significance) with age. It decreases from $\bar{N}_p$~3 for stars younger than 1 Gyr to $\bar{N}_p$~1.8 for stars about 8 Gyr. (3) The number of planets per star ($\eta=F_{\text{Kep}}\times\bar{N}_p$) also shows a decreasing trend (~2-3$\sigma$ significance). It decreases from $\eta$~1.6-1.7 for young stars to $\eta$~1.0 for old stars. (4) The mutual orbital inclination of the planets ($\sigma_{i,k}$) increases from $1.2^{+1.4}_{-0.5}$ to $3.5^{+8.1}_{-2.3}$ as stars aging from 0.5 to 8 Gyr with a best fit of $\log{\sigma_{i,k}}=0.2+0.4\times\log{\frac{\text{Age}}{\text{1Gyr}}}$. Interestingly, the Solar System also fits such a trend. The nearly independence of $F_{\text{Kep}}$~50% on age implies that planet formation is robust and stable across the Galaxy history. The age dependence of $\bar{N}_p$ and $\sigma_{i,k}$ demonstrates planetary architecture is evolving, and planetary systems generally become dynamically hotter with fewer planets as they age.

We present zephyr, a novel method that integrates cutting-edge normalizing flow techniques into a mixture density estimation framework, enabling the effective use of heterogeneous training data for photometric redshift inference. Compared to previous methods, zephyr demonstrates enhanced robustness for both point estimation and distribution reconstruction by leveraging normalizing flows for density estimation and incorporating careful uncertainty quantification. Moreover, zephyr offers unique interpretability by explicitly disentangling contributions from multi-source training data, which can facilitate future weak lensing analysis by providing an additional quality assessment. As probabilistic generative deep learning techniques gain increasing prominence in astronomy, zephyr should become an inspiration for handling heterogeneous training data while remaining interpretable and robustly accounting for observational uncertainties.

We numerically construct a series of axisymmetric rotating magnetic wind solutions, aiming at exploring the observation properties of massive white dwarf (WD) merger remnants with a strong magnetic field, a fast spin, and an intense mass loss, as inferred for WD J005311. We investigate the magnetospheric structure and the resultant spin-down torque exerted to the merger remnant with respect to the surface magnetic flux $\Phi_*$, spin angular frequency $\Omega_*$ and the mass loss rate $\dot M$. We confirm that the wind properties for $\sigma \equiv \Phi^2_* \Omega_*^2/\dot M v_\mathrm{esc}^3 \gtrsim 1$ significantly deviate from those of the spherical Parker wind, where $v_\mathrm{esc}$ is the escape velocity at stellar surface. For such a rotating magnetic wind sequence, we find: (i) quasi-periodic mass eruption triggered by magnetic reconnection along with the equatorial plane (ii) a scaling relation for the spin-down torque $T \approx (1/2) \times \dot{M} \Omega_* R^2_* \sigma^{1/4}$. We apply our results to discuss the spin-down evolution and wind anisotropy of massive WD merger remnants, the latter of which could be probed by a successive observation of WD J005311 using Chandra.

We extend the analysis presented in Contini et al. 2023 to higher redshifts, up to $z=2$, by focusing on the relation between the intracluster light (ICL) fraction and the halo mass, its dependence with redshift, role played by the halo concentration and formation time, in a large sample of simulated galaxy groups/clusters with $13\lesssim \log M_{halo} \lesssim 15$. Moreover, a key focus is to isolate the relative contributions provided by the main channels for the ICL formation to the total amount. The ICL fraction at higher redshift is weakly dependent on halo mass, and comparable with that at the present time, in agreement with recent observations. Stellar stripping, mergers and pre-processing are the major responsible channels of the ICL formation, with stellar stripping that accounts for $\sim 90\%$ of the total ICL, regardless of halo mass and redshift. Pre-processing is an important process for clusters to accrete already formed ICL. The diffuse component forms very early, $z\sim 0.6$, and its formation depends on both concentration and formation time of the halo, with more concentrated and earlier formed haloes that assemble their ICL earlier than later formed ones. The efficiency of this process is independent of halo mass, but increases with decreasing redshift, which implies that stellar stripping becomes more important with time as the concentration increases. This highlights the link between the ICL and the dynamical state of a halo: groups/clusters that have a higher fraction of diffuse light are more concentrated, relaxed and in an advanced stage of growth.

The nearby dwarf galaxy POX 52 at $z = 0.021$ hosts an active galactic nucleus (AGN) with a black-hole (BH) mass of $M_{\rm BH} \sim 10^{5-6} M_\odot$ and an Eddington ratio of $\sim$ 0.1-1. This object provides the rare opportunity to study both AGN and host-galaxy properties in a low-mass highly accreting system. To do so, we collected its multi-wavelength data from X-ray to radio. First, we construct a spectral energy distribution, and by fitting it with AGN and host-galaxy components, we constrain AGN-disk and dust-torus components. Then, while considering the AGN-disk emission, we decompose optical HST images. As a result, it is found that a classical bulge component is probably present, and its mass ($M_{\rm bulge}$) is consistent with an expected value from a local relation. Lastly, we analyze new quasi-simultaneous X-ray (0.2-30 keV) data obtained by NuSTAR and XMM-Newton. The X-ray spectrum can be reproduced by multi-color blackbody, warm and hot coronae, and disk and torus reflection components. Based on this, the spin is estimated to be $a_{\rm spin} = 0.998_{-0.814}$, which could suggest that most of the current BH mass was achieved by prolonged mass accretion. Given the presence of the bulge, POX 52 would have undergone a galaxy merger, while the $M_{\rm BH}$-$M_{\rm bulge}$ relation and the inferred prolonged accretion could suggest that AGN feedback occurred. Regarding the AGN structure, the spectral slope of the hot corona, its relative strength to the bolometric emission, and the torus structure are found to be consistent with Eddington-ratio dependencies found for nearby AGNs.

The known near-Earth object (NEO) population consists of over 32,000 objects, with a yearly discovery rate of over 3000 NEOs per year. An essential component of the next generation of NEO surveys is an understanding of the population of known objects, including an accounting of the discovery rate per year as a function of size. Using a near-Earth asteroid (NEA) reference model developed for NASA's NEO Surveyor (NEOS) mission and a model of the major current and historical ground-based surveys, an estimate of the current NEA survey completeness as a function of size and absolute magnitude has been determined (termed the Known Object Model; KOM). This allows for understanding of the intersection of the known catalog of NEAs and the objects expected to be observed by NEOS. The current NEA population is found to be $\sim38\%$ complete for objects larger than 140m, consistent with estimates by Harris & Chodas (2021). NEOS is expected to catalog more than two thirds of the NEAs larger than 140m, resulting in $\sim76\%$ of NEAs cataloged at the end of its 5 year nominal survey (Mainzer et al, 2023}, making significant progress towards the US Congressional mandate. The KOM estimates that $\sim77\%$ of the currently cataloged objects will be detected by NEOS, with those not detected contributing $\sim9\%$ to the final completeness at the end its 5 year mission. This model allows for placing the NEO Surveyor mission in the context of current surveys to more completely assess the progress toward the goal of cataloging the population of hazardous asteroids.

The isotopic ratios are good tools for probing the stellar nucleosynthesis and chemical evolution. We performed high-sensitivity mapping observations of the J=7-6 rotational transitions of OCS, OC34S, O13CS, and OC33S toward the Galactic Center giant molecular cloud, Sagittarius B2 (Sgr B2) with IRAM 30m telescope. Positions with optically thin and uncontaminated lines are chosen to determine the sulfur isotope ratios. A 32S/34S ratio of 17.1\pm0.9 was derived with OCS and OC34S lines, while 34S/33S ratio of 6.8\pm1.9 was derived directly from integrated intensity ratio of OC34S and OC33S. With independent and accurate measurements of 32S/34S ratio, our results confirm the termination of the decreasing trend of 32S/34S ratios toward the Galactic Center, suggesting a drop in the production of massive stars at the Galactic centre.

Data from the Michelson Doppler Imager (MDI) and Helioseismic and Magnetic Imager (HMI) are analyzed from 1996 to 2023 to investigate tilt angles ($\gamma$) of bipolar magnetic regions and Joy's Law for Cycles 23, 24, and a portion of 25. The HMI radial magnetic field ($B_{r}$) and MDI magnetogram ($B_{los}$) data are used to calculate ($\gamma$) using the flux-weighted centroids of the positive and negative polarities. Each AR is only sampled once. The analysis includes only Beta ($\beta$)-class active regions since computing $\gamma$ of complex active regions is less meaningful. During the emergence of the ARs, we find that the average tilt angle ($\bar{\gamma}$) increases from 3.30$^{\circ}\pm$0.75 when 20\% of the flux has emerged to 6.79$^{\circ}\pm$0.66 when the ARs are at their maximum flux. Cycle 24 had a larger average tilt $\bar{\gamma}_{24}$=6.67$\pm$0.66 than Cycle 23, $\bar{\gamma}_{23}$=5.11$\pm$0.61. There are persistent differences in $\bar{\gamma}$ in the hemispheres with the southern hemisphere having higher ${\bar{\gamma}}$ in Cycles 23 and 24 but the errors are such that these differences are not statistically significant.

Currently astrometric microlensing is the only tool that can directly measure the mass of a single star, it can also help us to detect compact objects like isolated neutron stars and black holes. The number of microlensing events that are being predicted and reported is increasing. In the paper, the potential lens stars are selected from three types of stars, high-proper-motion stars, nearby stars and high-mass stars. For each potential lens star, we select a larger search scope to find possible matching sources to avoid missing events as much as possible. Using Gaia DR3 data, we predict 4500 astrometric microlensing events with signal>0.1mas that occur between J2010.0 and J2070.0, where 1664 events are different from those found previously. There are 293 lens stars that can cause two or more events, where 5 lens stars can cause more than 50 events. We find that 116 events have the distance of background stars from the proper motion path of lens stars more than 8 arcsec in the reference epoch, where the maximum distance is 16.6 arcsec, so the cone search method of expanding the search range of sources for each potential lens star can reduce the possibility of missing events.

Numerical simulations and observations show that galaxies are not uniformly distributed. In cosmology, the largest known structures in the universe are galaxy filaments formed from the hierarchical clustering of galaxies due to gravitational forces. These consist of "walls" and "bridges" that connect clusters. Here, we use graph theory to model these structures as Euclidean networks in three-dimensional space. Using percolation theory, cosmological graphs are reduced based on the valency of nodes to reveal the inner, most robust structural formation. By constraining the network, we then find thresholds for physical features, such as length-scale and density, at which galaxy filament clusters are identified.

An event-based maximum likelihood method for handling X-ray polarimetry data is extended to include the effects of background and nonuniform sampling of the possible position angle space. While nonuniform sampling in position angle space generally introduces cross terms in the uncertainties of polarization parameters that could create degeneracies, there are interesting cases that engender no bias or parameter covariance. When including background in Poisson-based likelihood formulation, the formula for the minimum detectable polarization (MDP) has nearly the same form as for the case of Gaussian statistics derived by Elsner et al. (2012) in the limiting case of an unpolarized signal. A polarized background is also considered, which demonstrably increases uncertainties in source polarization measurements. In addition, a Kolmogorov-style test of the event position angle distribution is proposed that can provide an unbinned test of models where the polarization angle in Stokes space depends on event characteristics such as time or energy.

The dust in planet-forming disks evolve rapidly through growth and radial drift, and external photoevaporation also contributes to this evolution in massive star-forming regions. We test whether the presence of substructures can explain the survival of the dust component and observed millimeter continuum emission in protoplanetary disks located within massive star-forming regions. We also characterize the dust content removed by the photoevaporative winds. For this, we performed hydrodynamical simulations of protoplanetary disks subject to irradiation fields of $F_{UV} = 10^2$, $10^3$, and $10^4\, G_0$, with different dust trap locations. We used the FRIED grid to derive the mass loss rate for each irradiation field and disk properties, and then measure the evolution of the dust mass over time. For each simulation we estimate continuum emission at $\lambda = 1.3\, \textrm{mm}$ along with the radii encompassing $90\%$ of the continuum flux, and characterize the dust size distribution entrained in the photoevaporative winds, along with the resulting far-ultraviolet (FUV) cross section. Our simulations show that the presence of dust traps can extend the lifetime of the dust component of the disk to a few millionyears if the FUV irradiation is $F_{UV} \lesssim 10^3 G_0$, but only if the dust traps are located inside the photoevaporative truncation radius. The dust component of a disk quickly disperse if the FUV irradiation is strong ($10^4\, G_0$) or if the substructures are located outside the photoevaporation radius. We do find however, that the dust grains entrained with the photoevaporative winds may result in an absorption FUV cross section of $\sigma \approx 10^{-22}\, \textrm{cm}^2$ at early times of evolution ($<$0.1 Myr), which is enough to trigger a self-shielding effect that reduces the total mass loss rate, and slow down the disk dispersal in a negative feedback loop process.

We introduce the notion of a Bayesian analysis motivated `reliability' that gives a truer distinction of cusp-core and other halo-parameters (like mass-concentration) in an ensemble of observed galaxies. Our approach goes beyond the standard statistical techniques of parameter estimation and model fitting. We create hundreds of thousands of realistic mock SPARC RCs, with both cuspy and cored DM density profiles as model inputs. These RCs carefully incorporate the details of SPARC data such as the nature of observed uncertainties and different sources of scatters arising from observation, presence of baryons, DM mass-concentration, etc. Bayesian analysis of these mock RCs enables us to reconstruct and identify the parameter space in galaxy observable and theory where one can venture beyond best-fits to a preferred DM halo model or model selections between different density models. We find that it is imperative to choose low stellar surface density ($\Sigma_{\star}$) galaxies for reliable cusp-vs-core distinction; for example, RC data for galaxies with $\Sigma_{\star} \leq 2.5$ is needed for a 75\% confidence in distinguishing cusps from cores. Similarly, we also find that for correct estimations of the halo masses and concentrations, the RCs need to be measured to at least a radial distance $\geq 0.8r_s$ where $r_s$ is the scale radii of the corresponding DM halo density profiles. Out of the total $\sim$ 135 SPARC galaxies, using our reliability criteria, we find that only 21 RCs clear the bar to be used for any unbiased cusp-core distinction as well as DM halo mass-concentration estimates at $\geq$75\% reliability confidence level. With $\geq$66\% ( $\geq$50\%) reliability settings, the sample size increases to 44 (59). Interestingly, in the $\geq 75$\% reliable subsample, there are 5 times more galaxies that are reliably cored than cuspy. [Abridged]

In this study, we investigate how the baryonic effects vary with scale and local density environment mainly by utilizing a novel statistic, the environment-dependent wavelet power spectrum (env-WPS). With four state-of-the-art cosmological simulation suites, EAGLE, SIMBA, Illustris, and IllustrisTNG, we compare the env-WPS of the total matter density field between the hydrodynamic and dark matter-only (DMO) runs at $z=0$. We find that the clustering is most strongly suppressed in the emptiest environment of $\rho_\mathrm{m}/\bar\rho_\mathrm{m}<0.1$ with maximum amplitudes $\sim67-89$ per cent on scales $\sim1.86-10.96\ h\mathrm{Mpc}^{-1}$, and less suppressed in higher density environments on small scales (except Illustris). In the environments of $\rho_\mathrm{m}/\bar\rho_\mathrm{m}\geqslant0.316$ ($\geqslant10$ in EAGLE), the feedbacks also lead to enhancement features at intermediate and large scales, which is most pronounced in the densest environment of $\rho_\mathrm{m}/\bar\rho_\mathrm{m}\geqslant100$ and reaches a maximum $\sim 7-15$ per cent on scales $\sim0.87-2.62\ h\mathrm{Mpc}^{-1}$ (except Illustris). The baryon fraction of the local environment decreases with increasing density, denoting the feedback strength, and potentially explaining some differences between simulations. We also measure the volume and mass fractions of local environments, which are affected by $\gtrsim 1$ per cent due to baryon physics. In conclusion, our results reveal that the baryonic processes can change the overall cosmic structure greatly over the scales of $k>0.1\ h\mathrm{Mpc}^{-1}$. These findings warrant further investigation and confirmation by using much larger or more numerous simulations, comparing different cosmic structure classifications, or altering a specific process in the simulation, e.g. the star formation, black hole accretion, stellar feedback or active galactic nucleus feedback.

Stars and planets move supersonically in a gaseous medium during planetary engulfment, stellar interactions and within protoplanetary disks. For a nearly uniform medium, the relevant parameters are the Mach number and the size of the body, $R$, relative to its accretion radius, $R_A$. Over many decades, numerical and analytical work has characterized the flow, the drag on the body and the possible suite of instabilities. Only a limited amount of work has treated the stellar boundary as it is in many of these astrophysical settings, a hard sphere at $R$. Thus we present new 3-D Athena++ hydrodynamic calculations for a large range of parameters. For $R_A\ll R$, the results are as expected for pure hydrodynamics with minimal impact from gravity, which we verify by comparing to experimental wind tunnel data in air. When $R_A\approx R$, a hydrostatically-supported separation bubble forms behind the gravitating body, exerting significant pressure on the sphere and driving a recompression shock which intersects with the bow shock. For $R_A\gg R$, the bubble transitions into an isentropic, spherically-symmetric halo, as seen in earlier works. These two distinct regimes of flow morphology may be treated separately in terms of their shock stand-off distance and drag coefficients. Most importantly for astrophysical applications, we propose a new formula for the dynamical friction which depends on the ratio of the shock stand-off distance to $R_A$. That exploration also reveals the minimum size of the simulation domain needed to accurately capture the deflection of incoming streamlines due to gravity.

This paper investigates the spin-up of a mass-accreting star in a close binary system passing through the first stage of mass exchange in the Hertzsprung gap. Inside an accreting star, angular momentum is carried by meridional circulation and shear turbulence. The circulation carries part of the angular momentum entering the accretor to its surface. The greater the rate of arrival of angular momentum in the accretor is, the greater this part. It is assumed that this part of the angular momentum can be removed by the disk further from the accretor. If the angular momentum in the matter entering the accretor is more than half the Keplerian value, then the angular momentum obtained by the accretor during mass exchange stage does not depend on the rate of arrival of angular momentum. The accretor may have the characteristics of a Be-star immediately after the end of mass exchange.

Ncorpi$\mathcal{O}$N is a $N$-body software developed for the time-efficient integration of collisional and fragmenting systems of planetesimals or moonlets orbiting a central mass. It features a fragmentation model, based on crater scaling and ejecta models, able to realistically simulate a violent impact. The user of Ncorpi$\mathcal{O}$N can choose between four different built-in modules to compute self-gravity and detect collisions. One of these makes use of a mesh-based algorithm to treat mutual interactions in $\mathcal{O}(N)$ time. Another module, much more efficient than the standard Barnes-Hut tree code, is a $\mathcal{O}(N)$ tree-based algorithm called FalcON. It relies on fast multipole expansion for gravity computation and we adapted it to collision detection as well. Computation time is reduced by building the tree structure using a three-dimensional Hilbert curve. For the same precision in mutual gravity computation, Ncorpi$\mathcal{O}$N is found to be up to 25 times faster than the famous software REBOUND. Ncorpi$\mathcal{O}$N is written entirely in the C language and only needs a C compiler to run. A python add-on, that requires only basic python libraries, produces animations of the simulations from the output files. The name Ncorpi$\mathcal{O}$N, reminding of a scorpion, comes from the French $N$-corps, meaning $N$-body, and from the mathematical notation $\mathcal{O}(N)$, due to the running time of the software being almost linear in the total number $N$ of moonlets. Ncorpi$\mathcal{O}$N is designed for the study of accreting or fragmenting disks of planetesimal or moonlets. It detects collisions and computes mutual gravity faster than REBOUND, and unlike other $N$-body integrators, it can resolve a collision by fragmentation. The fast multipole expansions are implemented up to order six to allow for a high precision in mutual gravity computation.

Simulating the evolution of the gravitational N-body problem becomes extremely computationally expensive as N increases since the problem complexity scales quadratically with the number of bodies. We study the use of Artificial Neural Networks (ANNs) to replace expensive parts of the integration of planetary systems. Neural networks that include physical knowledge have grown in popularity in the last few years, although few attempts have been made to use them to speed up the simulation of the motion of celestial bodies. We study the advantages and limitations of using Hamiltonian Neural Networks to replace computationally expensive parts of the numerical simulation. We compare the results of the numerical integration of a planetary system with asteroids with those obtained by a Hamiltonian Neural Network and a conventional Deep Neural Network, with special attention to understanding the challenges of this problem. Due to the non-linear nature of the gravitational equations of motion, errors in the integration propagate. To increase the robustness of a method that uses neural networks, we propose a hybrid integrator that evaluates the prediction of the network and replaces it with the numerical solution if considered inaccurate. Hamiltonian Neural Networks can make predictions that resemble the behavior of symplectic integrators but are challenging to train and in our case fail when the inputs differ ~7 orders of magnitude. In contrast, Deep Neural Networks are easy to train but fail to conserve energy, leading to fast divergence from the reference solution. The hybrid integrator designed to include the neural networks increases the reliability of the method and prevents large energy errors without increasing the computing cost significantly. For this problem, the use of neural networks results in faster simulations when the number of asteroids is >70.

To probe star formation processes, we present a multi-scale and multi-wavelength investigation of the `Snake' nebula/infrared dark cloud G11.11$-$0.12 (hereafter, G11; length $\sim$27 pc). Spitzer images hint at the presence of sub-filaments (in absorption), and reveal four infrared-dark hub-filament system (HFS) candidates (extent $<$ 6 pc) toward G11, where massive clumps ($>$ 500 $M_{\odot}$) and protostars are identified. The $^{13}$CO(2-1), C$^{18}$O(2-1), and NH$_{3}$(1,1) line data reveal a noticeable velocity oscillation toward G11, as well as its left part (or part-A) around V$_{lsr}$ of 31.5 km s$^{-1}$, and its right part (or part-B) around V$_{lsr}$ of 29.5 km s$^{-1}$. The common zone of these cloud components is investigated toward the center's G11 housing one HFS. Each cloud component hosts two sub-filaments. In comparison to part-A, more ATLASGAL clumps are observed toward part-B. The JWST near-infrared images discover one infrared-dark HFS candidate (extent $\sim$0.55 pc) around the massive protostar G11P1 (i.e., G11P1-HFS). Hence, the infrared observations reveal multiple infrared-dark HFS candidates at multi-scale in G11. The ALMA 1.16 mm continuum map shows multiple finger-like features (extent $\sim$3500-10000 AU) surrounding a dusty envelope-like feature (extent $\sim$18000 AU) toward the central hub of G11P1-HFS. Signatures of forming massive stars are found toward the center of the envelope-like feature. The ALMA H$^{13}$CO$^{+}$ line data show two cloud components with a velocity separation of $\sim$2 km s$^{-1}$ toward G11P1. Overall, the collision process, the ``fray and fragment'' mechanism, and the ``global non-isotropic collapse'' scenario seem to be operational in G11.

This paper is a parametrization of the equation of state (EoS) parameter of dark energy (DE), which is parameterized using Square-Root (SR) form i.e. $\omega _{SR}=\text{$\omega _{0}$}+\text{$\omega _{1}$}\frac{z}{\sqrt{z^{2}+1}}$, where $\omega _{0}$ and $\omega _{1}$ are free constants. This parametrization will be examined in the context of the recently suggested $f(Q)$ gravity theory as an alternative to General Relativity (GR), in which gravitational effects are attributed to the non-metricity scalar $Q$ with the functional form $f(Q)=Q+\alpha Q^{n}$, where $\alpha$ and $n$ are arbitrary constants. We derived observational constraints on model parameters using the Hubble dataset with 31 data points and the Supernovae (SNe) dataset from the Pantheon samples compilation dataset with 1048 data points. For the current model, the evolution of the deceleration parameter, density parameter, EoS for DE, and $Om(z)$ diagnostic have all been investigated. It has been shown that the deceleration parameter favors the current accelerated expansion phase. It has also been shown that the EoS parameter for DE has a quintessence nature at this time.

We present a pan-chromatic study of AT2017bcc, a nuclear transient that was discovered in 2017 within the skymap of a reported burst-like gravitational wave candidate, G274296. It was initially classified as a superluminous supernova, and then reclassified as a candidate tidal disruption event. Its optical light curve has since shown ongoing variability with a structure function consistent with that of an active galactic nucleus, however earlier data shows no variability for at least 10 years prior to the outburst in 2017. The spectrum shows complex profiles in the broad Balmer lines: a central component with a broad blue wing, and a boxy component with time-variable blue and red shoulders. The H$\alpha$ emission profile is well modelled using a circular accretion disc component, and a blue-shifted double Gaussian which may indicate a partially obscured outflow. Weak narrow lines, together with the previously flat light curve, suggest that this object represents a dormant galactic nucleus which has recently been re-activated. Our time-series modelling of the Balmer lines suggests that this is connected to a disturbance in the disc morphology, and we speculate this could involve a sudden violent event such as a tidal disruption event involving the central supermassive black hole. Although we find that the redshifts of AT2017bcc ($z=0.13$) and G274296 ($z>0.42$) are inconsistent, this event adds to the growing diversity of both nuclear transients and multi-messenger contaminants.

We present self-consistent three-dimensional core-collapse supernova simulations of a rotating $20M_\odot$ progenitor model with various initial angular velocities from $0.0$ to $4.0$ rad s$^{-1}$ using a smoothed particle hydrodynamics code, SPHYNX, and a grid-based hydrodynamics code, FLASH. We identify two strong gravitational-wave features, with peak frequencies of $\sim300$ Hz and $\sim1.3$ kHz in the first $100$ ms postbounce. We demonstrate that these two features are associated with the $m=1$ deformation from the proto-neutron star (PNS) modulation induced by the low-$T/|W|$ instability, regardless of the simulation code. The $300$ Hz feature is present in models with an initial angular velocity between $1.0$ and $4.0$ rad s$^{-1}$, while the $1.3$ kHz feature is present only in a narrower range, from $1.5$ to $3.5$ rad s$^{-1}$. We show that the $1.3$ kHz signal originates from the high-density inner core of the PNS, and the $m=1$ deformation triggers a strong asymmetric distribution of electron anti-neutrinos. In addition to the $300$ Hz and $1.3$ kHz features, we also observe one weaker but noticeable gravitational-wave feature from higher-order modes in the range between $1.5$ and $3.5$ rad s$^{-1}$. Its peak frequency is around $800$ Hz initially and gradually increases to $900-1000$ Hz. Therefore, in addition to the gravitational bounce signal, the detection of the $300$ Hz, $1.3$ kHz, the higher-order mode, and even the related asymmetric emission of neutrinos, could provide additional diagnostics to estimate the initial angular velocity of a collapsing core.

Observations of clusters suffer from issues such as completeness, projection effects, resolving individual stars and extinction. As such, how accurate are measurements and conclusions are likely to be? Here, we take cluster simulations (Westerlund2- and Orion- type), synthetically observe them to obtain luminosities, accounting for extinction and the inherent limits of Gaia, then place them within the real Gaia DR3 catalogue. We then attempt to rediscover the clusters at distances of between 500pc and 4300pc. We show the spatial and kinematic criteria which are best able to pick out the simulated clusters, maximising completeness and minimising contamination. We then compare the properties of the 'observed' clusters with the original simulations. We looked at the degree of clustering, the identification of clusters and subclusters within the datasets, and whether the clusters are expanding or contracting. Even with a high level of incompleteness (e.g. $<2\%$ stellar members identified), similar qualitative conclusions tend to be reached compared to the original dataset, but most quantitative conclusions are likely to be inaccurate. Accurate determination of the number, stellar membership and kinematic properties of subclusters, are the most problematic to correctly determine, particularly at larger distances due to the disappearance of cluster substructure as the data become more incomplete, but also at smaller distances where the misidentification of asterisms as true structure can be problematic. Unsurprisingly, we tend to obtain better quantitative agreement of properties for our more massive Westerlund2-type cluster. We also make optical style images of the clusters over our range of distances.

$Aims:$ We used the nearby Coma Cluster as a laboratory in order to probe the impact of ram pressure on star formation as well as to constrain the characteristic timescales and velocities for the stripping of the non-thermal ISM. $Methods:$ We used high-resolution ($6.5'' \approx 3\,\mathrm{kpc}$), multi-frequency ($144\,\mathrm{MHz} - 1.5\,\mathrm{GHz}$) radio continuum imaging of the Coma Cluster to resolve the low-frequency radio spectrum across the discs and tails of 25 ram pressure stripped galaxies. With resolved spectral index maps across these galaxy discs, we constrained the impact of ram pressure perturbations on galaxy star formation. We measured multi-frequency flux-density profiles along each of the ram pressure stripped tails in our sample. We then fit the resulting radio continuum spectra with a simple synchrotron aging model. $Results:$ We showed that ram pressure stripped tails in Coma have steep ($-2 \lesssim \alpha \lesssim -1$) spectral indices. The discs of galaxies undergoing ram pressure stripping have integrated spectral indices within the expected range for shock acceleration from supernovae ($-0.8 \lesssim \alpha \lesssim -0.5$), though there is a tail towards flatter values. In a resolved sense, there are gradients in spectral index across the discs of ram pressure stripped galaxies in Coma. These gradients are aligned with the direction of the observed radio tails, with the flattest spectral indices being found on the `leading half'. From best-fit break frequencies we estimated the projected plasma velocities along the tail to be on the order of hundreds of kilometers per second, with the precise magnitude depending on the assumed magnetic field strength.

As the most energetic explosions in the Universe, gamma-ray bursts (GRBs) are commonly believed to be generated by relativistic jets. Recent observational evidence suggests that the jets producing GRBs are likely to have a structured nature. Some studies have suggested that non-axisymmetric structured jets may be formed through internal non-uniform magnetic dissipation processes or the precession of the central engine. In this study, we analyze the potential characteristics of GRB afterglows within the framework of non-axisymmetric structured jets. We simplify the profile of the asymmetric jet as a step function of the azimuth angle, dividing the entire jet into individual elements. By considering specific cases, we demonstrate that the velocity, energy, and line-of-sight direction of each jet element can greatly affect the behaviour of the overall light curve. The radiative contributions from multiple elements may lead to the appearance of multiple distinct peaks or plateaus in the light curve. Furthermore, fluctuations in the rising and declining segments of each peak can be observed. These findings establish a theoretical foundation for future investigations into the structural characteristics of GRBs by leveraging GRB afterglow data.

Supernova explosions attributed to the unseen companion in several binary systems identified by the Third Gaia Data Release (Gaia DR3) may be responsible for a number of well-known and well-studied features in the radio sky, including the Low-Latitude-Intermediate-Velocity Arch and the North-Celestial-Pole Loop. Slices from the Longitude-Latitude-Velocity data cube of the $\lambda$-21-cm galactic neutral hydrogen HI4PI survey (HI4PI Collaboration et al. 2016) show multiple signatures of an expanding shell. The source of this expansion, which includes the Low-Latitude-Intermediate-Velocity Arch on the approaching side, may be the neutron star candidate in the Gaia DR3 1093757200530267520 binary. If we make the simplifying assumptions that the expansion of the cavity is uniform and spherically symmetric, then the explosion took place about 700,000 years ago. The momentum is in reasonable agreement with recent model estimates for a supernova this old. The HI on the receding side of this cavity is interacting with the gas approaching us on the near side of a second cavity. The North-Celestial-Pole Loop appears to be located at the intersection of these two expanding features. The neutron star candidate in the Gaia DR3 1144019690966028928 binary may be (in part) responsible for this cavity. Explosions from other candidates may account for the observed elongation along the line of sight of this second cavity. We can use the primary star in these binaries to anchor the distances to the Low-Latitude-Intermediate-Velocity Arch and North-Celestial-Pole Loop, which are about 167 and about 220 pc, respectively.

We describe the interaction of parallel-propagating Alfv\'en waves with ion-acoustic waves and other Alfv\'en waves, in magnetized, high-$\beta$ collisionless plasmas. This is accomplished through a combination of analytical theory and numerical fluid simulations of the Chew-Goldberger-Low (CGL) magnetohydrodynamic (MHD) equations closed by Landau-fluid heat fluxes. An asymptotic ordering is employed to simplify the CGL-MHD equations and derive solutions for the deformation of an Alfv\'en wave that results from its interaction with the pressure anisotropy generated either by an ion-acoustic wave or another, larger-amplitude Alfv\'en wave. The difference in timescales of acoustic and Alfv\'enic fluctuations at high-$\beta$ means that interactions that are local in wavenumber space yield little modification to either mode within the time it takes the acoustic wave to Landau damp away. Instead, order-unity changes in the amplitude of Alfv\'enic fluctuations can result after interacting with frequency-matched acoustic waves. Additionally, we show that the propagation speed of an Alfv\'en-wave packet in an otherwise homogeneous background is a function of its self-generated pressure anisotropy. This allows for the eventual interaction of separate co-propagating Alfv\'en-wave packets of differing amplitudes. The results of the CGL-MHD simulations agree well with these predictions, suggesting that theoretical models relying on the interaction of these modes should be reconsidered in certain astrophysical environments. Applications of these results to weak Alfv\'enic turbulence and to the interaction between the compressive and Alfv\'enic cascades in strong, collisionless turbulence are also discussed.

Anomalous Cepheids (ACEPs) are intermediate mass metal-poor pulsators mostly discovered in dwarf galaxies of the Local Group. However, recent Galactic surveys, including the Gaia DR3, found a few hundreds of ACEPs in the Milky Way. Their origin is not well understood. We aim to investigate the origin and evolution of Galactic ACEPs by studying for the first time the chemical composition of their atmospheres. We used UVES@VLT to obtain high-resolution spectra for a sample of 9 ACEPs belonging to the Galactic halo. We derived the abundances of 12 elements, including C, Na, Mg, Si, Ca, Sc, Ti, Cr, Fe, Ni, Y, and Ba. We complemented these data with literature abundances for an additional three ACEPs that were previously incorrectly classified as type II Cepheids, thus increasing the sample to a total of 12 stars. All the investigated ACEPs have an iron abundance [Fe/H]$<-1.5$ dex as expected from theoretical predictions for these pulsators. The abundance ratios of the different elements to iron show that the ACEP's chemical composition is generally consistent with that of the Galactic halo field stars, except the Sodium, which is found overabundant in 9 out of the 11 ACEPs where it was measured, in close similarity with second-generation stars in the Galactic Globular Clusters. The same comparison with dwarf and ultra-faint satellites of the Milky Way reveals more differences than similarities so it is unlikely that the bulk of Galactic ACEPs originated in such a kind of galaxies which subsequently dissolved in the Galactic halo. The principal finding of this work is the unexpected overabundance of Sodium in ACEPs. We explored several hypotheses to explain this feature, finding that the most promising scenario is the evolution of low-mass stars in a binary system with either mass transfer or merging. Detailed modelling is needed to confirm this hypothesis.

We have measured structural parameters and radial color profiles of 108 ultra-diffuse galaxies (UDGs), carefully selected from six distant massive galaxy clusters in the Hubble Frontier Fields (HFF) in redshift range from 0.308 to 0.545. Our best-fitting GALFIT models show that the HFF UDGs have a median S\'ersic index of 1.09, which is close to 0.86 for local UDGs in the Coma cluster. The median axis-ratio value is 0.68 for HFF UDGs and 0.74 for Coma UDGs, respectively. The structural similarity between HFF and Coma UDGs suggests that they are the same kind of galaxies seen at different times and the structures of UDGs do not change at least for several billion years. By checking the distribution of HFF UDGs in the rest-frame $UVJ$ and $UVI$ diagrams, we find a large fraction of them are star-forming. Furthermore, a majority of HFF UDGs show small $\rm U-V$ color gradients within \,1\,*\,$R_{e,SMA}$ region, the fluctuation of the median radial color profile of HFF UDGs is smaller than 0.1\,mag, which is compatible to Coma UDGs. Our results indicate that cluster UDGs may fade or quench in a self-similar way, irrespective of the radial distance, in less than $\sim$ 4 Gyrs.

Scientific data collected at ESO's observatories are freely and openly accessible online through the ESO Science Archive Facility. In addition to the raw data straight out of the instruments, the ESO Science Archive also contains four million processed science files available for use by scientists and astronomy enthusiasts worldwide. ESO subscribes to the FAIR (Findable, Accessible, Interoperable, Reusable) guiding principles for scientific data management and stewardship. All data in the ESO Science Archive are distributed according to the terms of the Creative Commons Attribution 4.0 International licence (CC BY 4.0).

Ultra-slow-roll~(USR) inflation predicts an exponential amplification of scalar perturbations at small scales, which leads to a stochastic gravitational wave background~(SGWB) through the coupling of the scalar and tensor modes at the second-order expansion of the Einstein equation. In this work, we search for such a scalar-induced SGWB from the NANOGrav 15-year (NG15) dataset, and find that the SGWB from USR inflation could explain the observed data. We place constraints on the amplitude of the scalar power spectrum to $P_{\mathrm{Rp}} > 10^{-1.80}$ at $95\%$ confidence level (C.L.) at the scale of $k\sim 20\, \mathrm{pc}^{-1}$. We find that $\log_{10} P_{\mathrm{Rp}}$ degenerates with the peak scale $\log_{10} k_{\mathrm{p}}$. We also obtain the parameter space allowed by the data in the USR inflationary scenario, where the $e$-folding numbers of the duration of the USR phase has a lower limit $\Delta N > 2.80$ ($95\%$ C.L.) when the USR phase ends at $N\approx 20$. Since the priors for the model parameters %in the USR model are uncertain, we do not calculate the Bayes factors. Instead, to quantify the goodness of fit, we calculate the maximum values of the log-likelihood for USR inflation, bubble collision of the cosmological phase transition, and inspiraling supermassive black hole binaries (SMBHBs), respectively. Our results imply that the SGWB from USR inflation can fit the data better than the one from SMBHBs.

PRIMA (The PRobe for-Infrared Mission for Astrophysics) is a concept for a far-infrared (IR) observatory. PRIMA features a cryogenically cooled 1.8 m diameter telescope and is designed to carry two science instruments enabling ultra-high sensitivity imaging and spectroscopic studies in the 24 to 235 microns wavelength range. The resulting observatory is a powerful survey and discovery machine, with mapping speeds better by 2 - 4 orders of magnitude with respect to its far-IR predecessors. The bulk of the observing time on PRIMA should be made available to the community through a General Observer (GO) program offering 75% of the mission time over 5 years. In March 2023, the international astronomy community was encouraged to prepare authored contributions articulating scientific cases that are enabled by the telescope massive sensitivity advance and broad spectral coverage, and that could be performed within the context of GO program. This document, the PRIMA General Observer Science Book, is the edited collection of the 76 received contributions.

Theoretical models of accretion discs and observational data indicate that the X-ray emission from the inner parts of an accretion disc can irradiate its outer regions and induce a thermal wind, which carries away the mass and angular momentum from the disc. Our aim is to investigate the influence of the thermal wind on the outburst light curves of black hole X-ray binary systems. We carry out numerical simulations of a non-stationary disc accretion with wind using upgraded open code freddi. We assume that the wind launches only from the ionised part of the disc and may turn off if the latter shrinks fast enough. Our estimates of the viscosity parameter $\alpha$ are shifted downward compared to a scenario without a wind. Generally, correction of $\alpha$ depends on the spectral hardness of central X-rays and the disc outer radius, but unlikely to exceed a factor of 10 in the case of a black hole low-mass X-ray binary (BH LMXB). We fit 2002 outburst of BH LMXB 4U 1543-47 taking into account the thermal wind. The mass loss in the thermal wind is of order of the accretion rate on the central object at the peak of the outburst. New estimate of the viscosity parameter $\alpha$ for the accretion disc in this system is about two times lower than the previous one. Additionally, we calculate evolution of the number of hydrogen atoms towards 4U 1543-47 due to the thermal wind from the hot disc.

Dust emission is an important tool in studies of star-forming clouds, as a tracer of column density and indirectly via the dust evolution that is connected to the history and physical conditions of the clouds. We examine radiative transfer (RT) modelling of dust emission over an extended cloud region, using a filament in the Taurus molecular cloud as an example. We examine how well far-infrared observations can be used to determine both the cloud and the dust properties. Using different assumptions of the cloud shape, radiation field, and dust properties, we fit RT models to Herschel observations of the Taurus filament. Further comparisons are made with measurements of the near-infrared extinction. The models are used to examine the degeneracies between the different cloud parameters and the dust properties. The results show significant dependence on the assumed cloud structure and the spectral shape of the external radiation field. If these are constrained to the most likely values, the observations can be explained only if the dust far-infrared (FIR) opacity has increased by a factor of 2-3 relative to the values in diffuse medium. However, a narrow range of FIR wavelengths provides only weak evidence of the spatial variations in dust, even in the models covering several square degrees of a molecular cloud. The analysis of FIR dust emission is affected by several sources of uncertainty. Further constraints are therefore needed from observations at shorter wavelengths, especially regarding the trends in dust evolution.

Magnetically arrested accretion disks (MADs) around a rapidly rotating black hole (BH) have been proposed as a model for jetted tidal disruption events (TDEs). However, the dynamics of strongly magnetized disks in a more realistic simulation which can mimic the chaotic dynamics during a TDE have previously been unexplored. Here we employ global GRMHD simulations of a pre-existing MAD disk interacting with an injected TDE stream with impact parameter $\beta\equiv R_t/R_p=4-7$ to investigate how strongly magnetized TDEs differ from the standard MAD picture. We demonstrate for the first time that a MAD or semi-MAD state can be sustained and jets powered by the BH spin are produced in a TDE. We also demonstrate that the strength of the self-intersection shock depends on how dense the disk is relative to the stream, or the density contrast $f_\rho=\rho_d/\rho_s$. The jet or funnel can become significantly tilted (by $10-30^\circ$) due to the self-intersection outflow when $f_\rho \leq 0.1$. In models with a powerful jet and $f_\rho\leq 0.01$, the tilted jet interacts with and ultimately tilts the disk by as much as 23 degrees from the incoming stream. We illustrate that as $f_\rho$ increases, the tilt of the jet and disk is expected to realign with the BH spin once $f_\rho \geq 0.1$. We illustrate how the tilt can rapidly realign if $f_\rho$ increases rapidly and apply this to TDEs which have shown X-ray evolution on timescales of days-weeks.

We examined the past history of the three most detached TransNeptunian Objects (TNOs) -- Sedna, 2012 VP113, and Leleakuhonua (2015 TG387) -- the three clearest members of the dynamical class known as sednoids, with high perihelia distances $q$. By integrating backward their nominal (and a set of cloned) orbits for the Solar System's age, we surprisingly find that the only time all their apsidal lines tightly cluster was 4.5 Gyr ago, at perihelion longitude $\varpi$ of 200{\deg}. This "primordial alignment" is independent of the observational biases that contribute to the current on-sky clustering in the large-semimajor axis Kuiper Belt. If future sednoid discoveries confirm these findings, this strongly argues for an initial event during the planet formation epoch which imprinted this particular apsidal orientation on the early detached TNO population and then subsequently modified only by the simple precession from the 4 giant planets. If other sednoids also cluster around the same primordial value, various models suggesting a still present planet in the outer Solar System would be incompatible with this alignment. We inspected two scenarios that could potentially explain the primordial alignment. First, a rogue planet model (where another massive planet raises perihelia near its own longitude until ejection) naturally produces this signature. Alternatively, a close stellar passage early in Solar System history raises perihelia, but it is poor at creating strong apsidal clustering. We show that all other known $35<q<55$ au TNOs are either too perturbed or orbits are still too uncertain to provide evidence for or against this paradigm.

NGC 1068 is a nearby widely studied Seyfert II galaxy presenting radio, infrared, X- and $\gamma$-ray emission as well as strong evidence for high-energy neutrino emission. Recently, the evidence for neutrino emission could be explained in a multimessenger model in which the neutrinos originate from the corona of the active galactic nucleus (AGN). In this environment $\gamma$-rays are strongly absorbed, so that an additional contribution from e.g. the circumnuclear starburst ring is necessary. In this work, we discuss whether the radio jet can be an alternative source of the $\gamma$-rays between about $0.1$ and $100$ GeV as observed by Fermi-LAT. In particular, we include both leptonic and hadronic processes, i.e. accounting for inverse Compton emission and signatures from $pp$ as well as $p\gamma$ interactions. In order to constrain our calculations, we use VLBA and ALMA observations of the radio knot structures, which are spatially resolved at different distances from the supermassive black hole. Our results show that the best leptonic scenario for the prediction of the Fermi-LAT data is provided by the radio knot closest to the central engine. For that a magnetic field strength $\sim 1\,\text{mG}$ is needed as well as a strong spectral softening of the relativistic electron distribution at $(1-10)\,\text{GeV}$. However, we show that neither such a weak magnetic field strength nor such a strong softening is expected for that knot. A possible explanation for the $\sim$ 10 GeV $\gamma$-rays can be provided by hadronic pion production in case of a gas density $\gtrsim 10^4\,\text{cm}^{-3}$. Nonetheless, this process cannot contribute significantly to the low energy end of the Fermi-LAT range. We conclude that the emission sites in the jet are not able to explain the $\gamma$-rays in the whole Fermi-LAT energy band.

We confirm the planetary nature of TOI-5344 b as a transiting giant exoplanet around an M0 dwarf star. TOI-5344 b was discovered with the Transiting Exoplanet Survey Satellite photometry and confirmed with ground-based photometry (the Red Buttes Observatory 0.6m telescope), radial velocity (the Habitable-zone Planet Finder), and speckle imaging (the NN-Explore Exoplanet Stellar Speckle Imager). TOI-5344 b is a Saturn-like giant planet ($\rho = 0.80^{+0.17}_{-0.15}\ \text{g cm}^{-3}$) with a planetary radius of $9.7 \pm \ 0.5 \ \text{R}_{\oplus}$ ($0.87 \pm \ 0.04 \ \text{R}_{\text{Jup}}$) and a planetary mass of $135^{+17}_{-18} \text{M}_{\oplus}$ ($0.42^{+0.05}_{-0.06} \ \text{M}_{\text{Jup}}$). It has an orbital period of $3.792622 \pm 0.000010$ days and an orbital eccentricity of $0.06^{+0.07}_{-0.04}$. We measure a high metallicity for TOI-5344 of [Fe/H] = $0.48 \pm 0.12$, where the high metallicity is consistent with expectations from formation through core accretion. We compare the metallicity of the M-dwarf hosts of giant exoplanets to that of M-dwarf hosts of non-giants ($\lesssim 8\ \text{R}_{\oplus}$). While the two populations appear to show different metallicity distributions, quantitative tests are prohibited by various sample caveats.

Accretion onto black holes is one of the most efficient energy source in the Universe. Black hole accretion powers some of the most luminous objects in the universe, including quasars, active galactic nuclei, tidal disruption events, gamma-ray bursts, and black hole X-ray transients. In the present review, we give an astrophysical overview of black hole accretion processes, with a particular focus on black hole X-ray binary systems. In Section 1, we briefly introduce the basic paradigms of black hole accretion. Physics related to accretion onto black holes are introduced in Section 2. Models proposed for black hole accretion are discussed in this section, from the Shakura-Sunyaev thin disk accretion to the advective-dominated accretion flow. Observational signatures that make contact to stellar-mass black hole accretion are introduced in Section 3, including the spectral and fast variability properties. A short conclusion is given in Section 4.

Constrained cosmological simulations play an important role in modelling the local Universe, enabling investigation of the dark matter content of local structures and their formation histories. We introduce a method for determining the extent to which individual haloes are reliably reconstructed between constrained simulations, and apply it to the Constrained Simulations in BORG (CSiBORG) suite of $101$ high-resolution realisations across the posterior probability distribution of initial conditions from the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm. The method is based on the overlap of the initial Lagrangian patch of a halo in one simulation with those in another, and therefore measures the degree to which the haloes' particles are initially coincident. By this metric we find consistent reconstructions of $M\gtrsim10^{14}~M_\odot / h$ haloes across the CSiBORG simulations, indicating that the constraints from the BORG algorithm are sufficient to pin down the masses, positions and peculiar velocities of clusters to high precision. The effect of the constraints tapers off towards lower mass however, and the halo spins and concentrations are largely unconstrained at all masses. We document the advantages of evaluating halo consistency in the initial conditions, describe how the method may be used to quantify our knowledge of the halo field given galaxy survey data analysed through the lens of probabilistic inference machines such as BORG, and describe applications to matched but unconstrained simulations.

The first year of JWST has revealed a surprisingly large number of luminous galaxy candidates beyond $z>10$. While some galaxies are already spectroscopically confirmed, there is mounting evidence that a subsample of the candidates with particularly red inferred UV colors are in fact lower redshift contaminants.These interlopers are often found to be `HST-dark' or `optically-faint' galaxies at $z\sim2-6$, a population key to understanding dust-obscured star formation throughout cosmic time. This paper demonstrates the complementarity of ground-based mm-interferometry and JWST infrared imaging to unveil the true nature of red 1.5-2.0 $\mu$m dropouts that have been selected as ultra-high-redshift galaxy candidates. We present NOEMA Polyfix follow-up observations of four JWST red 1.5-2.0 $\mu$m dropouts selected by Yan et al. 2023 as ultra-high-redshift candidates in the PEARLS field. The new NOEMA observations constrain the rest-frame far-infrared continuum emission and efficiently discriminate between intermediate- and high-redshift solutions. We report $>10\sigma$ NOEMA continuum detections of all our target galaxies at observed frequencies of $\nu$=236 and 252 GHz, with FIR slopes indicating a redshift $z<5$. We model their optical-to-FIR spectral energy distribution (SED) with multiple SED codes, and find that they are not $z>10$ galaxies but instead dust-obscured, massive star-forming galaxies at $z\sim 2-4$. The contribution to the cosmic star-formation rate density of such sources is not negligible at $z\simeq 3.5$ ($\phi\gtrsim(1.9-4.4)\times10^{-3}\ \rm{cMpc}^{-3}$), in line with previous studies of optically-faint/sub-millimeter galaxies. This work showcases a new way to select intermediate- to high-redshift dust-obscured galaxies in JWST fields with minimal wavelength coverage to open a new window on obscured star-formation at intermediate redshifts .[abridged]

State-of-the-art surveys reveal that most massive stars in the universe evolve in close binaries. Massive stars in such systems are expected to develop aspherical envelopes due to tidal interactions and/or rotational effects. Recently, it was shown that point explosions in oblate stars can produce relativistic equatorial ring-like outflows. Moreover, since stripped-envelope stars in binaries can expand enough to fill their Roche lobes anew, it is likely that these stars die with a greater degree of asphericity than the oblate spheroid geometry previously studied. We investigate the effects of this asymmetry by studying the gas dynamics of axisymmetric point explosions in stars in various stages of filling their Roche lobes. We find that point explosions in these pear-shaped stars produce trans-relativistic ejecta that coalesces into bullets pointed both toward and away from the binary companion. We present this result and comment on key morphological differences between core-collapse explosions in spherical versus distorted stars in binary systems, effects on gravitational wave sources, and observational signatures that could be used to glean these explosion geometries from current and future surveys.

Stellar feedback plays a crucial role in star formation and the life cycle of molecular clouds. The intense star formation region 30 Doradus, which is located in the Large Magellanic Cloud (LMC), is a unique target for detailed investigation of stellar feedback owing to the proximity of the hosting galaxy and modern observational capabilities that together allow us to resolve individual molecular clouds $-$ nurseries of star formation. We study the impact of large-scale feedback on the molecular gas using the new observational data in the $^{12}$CO(3$-$2) line obtained with the APEX telescope. Our data cover an unprecedented area of 13.9 sq. deg. of the LMC disc with a spatial resolution of 5 pc and provide an unbiased view on the molecular clouds in the galaxy. Using this data, we locate molecular clouds in the disc of the galaxy, estimate their properties such as the areal number density, relative velocity and separation, width of the line profile, CO-line luminosity, size, virial mass and compare these properties between the clouds of 30 Doradus and those in the rest of the LMC disc. We find that compared with the rest of the observed molecular clouds in the LMC disc, those in 30 Doradus show the highest areal number density; they are spatially more clustered, move faster with respect to each other and feature larger linewidths. In parallel, we do not find statistically significant differences in such properties as the CO-line luminosity, size, and virial mass between the clouds of 30 Doradus and the rest of the observed field. We interpret our results as signatures of gas dispersal and fragmentation due to high-energy large-scale feedback.

In this work, we reconstruct the Hubble diagram using various data sets, including correlated ones, in Artificial Neural Networks (ANN). Using ReFANN, that was built for data sets with independent uncertainties, we expand it to include non-Guassian data points, as well as data sets with covariance matrices among others. Furthermore, we compare our results with the existing ones derived from Gaussian processes and we also perform null tests in order to test the validity of the concordance model of cosmology.

Ultralight dark photons are compelling dark matter candidates, but their allowed kinetic mixing with the Standard Model photon is severely constrained by requiring that the dark photons do not collapse into a cosmic string network in the early Universe. Direct detection in minimal production scenarios for dark photon dark matter is strongly limited, if not entirely excluded; discovery of sub-meV dark photon dark matter would therefore point to a nonminimal dark sector. We describe a model that evades such constraints, capable of producing cold dark photons in any parameter space accessible to future direct detection experiments. The associated production dynamics yield additional signatures in cosmology and small-scale structure, allowing for possible positive identification of this particular class of production mechanisms.

A short phenomenological account of the genesis and evolution of the universe is presented with emphasis on the primordial phases as well as its physical composition, i.e. dark matter and dark energy. We discuss Einstein's theory of General Relativity and its consequences for the birth of modern relativistic astrophysics. We introduce the Big-Bang theory of Mons. Lemaitre as well as the competing theory of the Steady State Universe of Fred Hoyle. Since Big-Bang theory appeared quite in agreement with Christian doctrine of creation, Pope Pius XII delivered a message to the pontifical Academy of Sciences in 1951 claiming a certain agreement between the creation account in the book of Genesis and the Big-Bang theory (a concordist view), a position which he did not repeat later. On the other hand, Lemaitre always kept separate the scientific and theological planes as two parallel "lines" never intersecting, i.e., as two complementary "magisteria". Similar kind of tensions, between science and theology, emerge also today with the Hartle-Hawking solution to the Wheeler-DeWitt equation in quantum cosmology and its related speculations. To avoid some sort of confusion between theological and physics concepts, we, briefly, summarise the concept of creation in Christian theology.

The suppression of relic gravitational waves due to their conversion into electromagnetic radiation in a cosmological magnetic field is studied. It is shown that the subsequent elimination of photons from the beam due to their interaction with the primary plasma prevents from the inverse restoration of the gravitational waves by the photons. The coupled system of equations describing gravitational and electromagnetic wave propagation in an arbitrary curved space-time and in external magnetic field is derived. The system of equations is solved numerically in Friedmann- LeMaitre-Robertson-Walker metric for the upper limit of the intergalactic magnetic field strength of 1 nGs. We conclude that the gravitational wave conversion into photons in the intergalactic magnetic field can significantly change the amplitude of the relic gravitational wave and their frequency spectrum.

$\gamma$-ray emission of blazars infer the presence of large-scale magnetic fields in the intergalactic medium, but their origin remains a mystery. Using recent data from MAGIC, H.E.S.S. and $\textit{Fermi}$-LAT, we investigate whether the large-scale magnetic fields in the intergalactic medium could have been generated by a first-order electroweak phase transition in the two-Higgs-doublet model (2HDM). We study two representative scenarios where we vary the initial conditions of the magnetic field and the plasma, assuming either a primordial magnetic field with maximal magnetic helicity or a primordial magnetic field with negligible magnetic helicity in a plasma with kinetic helicity. By considering a primordial magnetic field with maximal helicity and applying the conservative constraints derived from MAGIC and $\textit{Fermi}$-LAT data, we demonstrate that a first-order electroweak phase transition within the 2HDM may account for the observed intergalactic magnetic fields in the case of the strongest transitions. We show that this parameter space also predicts strong gravitational wave signals in the reach of space-based detectors such as LISA, providing a striking multi-messenger signal of the 2HDM.

We report on a new search for continuous gravitational waves from NS 1987A, the neutron star born in SN 1987A, using open data from Advanced LIGO and Virgo's third observing run (O3). The search covered frequencies from 35-1050 Hz, more than five times the band of the only previous gravitational wave search to constrain NS 1987A [B. J. Owen et al., ApJL 935, L7 (2022)]. It used an improved code and coherently integrated from 5.10 days to 14.85 days depending on frequency. No astrophysical signals were detected. By expanding the frequency range and using O3 data, this search improved on strain upper limits from the previous search and was sensitive at the highest frequencies to ellipticities of 1.6e-5 and r-mode amplitudes of 4.4e-4, both an order of magnitude improvement over the previous search and both well within the range of theoretical predictions.

Machine learning techniques can automatically identify outliers in massive datasets, much faster and more reproducible than human inspection ever could. But finding such outliers immediately leads to the question: which features render this input anomalous? We propose a new feature attribution method, Inverse Multiscale Occlusion, that is specifically designed for outliers, for which we have little knowledge of the type of features we want to identify and expect that the model performance is questionable because anomalous test data likely exceed the limits of the training data. We demonstrate our method on outliers detected in galaxy spectra from the Dark Energy Survey Instrument and find its results to be much more interpretable than alternative attribution approaches.

Andreev-Bashkin entrainment makes the hydrodynamics of the binary superfluid solution particularly interesting. We investigate stability and motion of quantum vortices in such system.

We perform the first numerical simulations modeling the inflow and outflow of magnetized plasma in the Kerr-Sen spacetime, which describes classical spinning black holes (BHs) in string theory. We find that the Blandford-Znajek (BZ) mechanism, which is believed to power astrophysical relativistic outflows or ``jets'', is valid even for BHs in an alternate theory of gravity, including near the extremal limit. The BZ mechanism releases outward Poynting-flux-dominated plasma as frame-dragging forces magnetic field lines to twist. However, for nonspinning BHs, where the frame-dragging is absent, we find an alternate powering mechanism through the release of gravitational potential energy during accretion. Outflows from non-spinning stringy BHs can be approximately $250\%$ more powerful as compared to Schwarzschild BHs, due to their relatively smaller event horizon sizes and, thus, higher curvatures. Finally, by constructing the first synthetic images of near-extremal non-Kerr BHs from time-dependent simulations, we find that these can be ruled out by horizon-scale interferometric images of accreting supermassive BHs.

The launching of astrophysical jets provides the most compelling observational evidence for direct extraction of black hole (BH) spin energy via the Blandford-Znajek (BZ) mechanism. Whilst it is know that spinning Kerr BHs within general relativity (GR) follow the BZ jet power relation, the nature of BH energy extraction in general theories of gravity has not been adequately addressed. This study performs the first comprehensive investigation of the BZ jet power relation by utilising a generalised BH spacetime geometry which describes parametric deviations from the Kerr metric of GR, yet recovers the Kerr metric in the limit that all deviation parameters vanish. Through performing and analysing an extensive suite of three-dimensional covariant magnetohydrodynamics (MHD) simulations of magnetised gas accretion onto these generalised BH spacetimes we find that the BZ jet power relation still holds, in some instances yielding jet powers far in excess of what can be produced by even extremal Kerr BHs. It is shown that the variation of the quadrupole moment of the BH can enhance or suppress the effects of BH spin, and by extension of frame-dragging. This variation greatly enhances or suppresses the observed jet power and underlying photon ring image asymmetry, introducing a previously unexplored yet important degeneracy in BH parameter inference.

Space-based gravitational wave detection is one of the most anticipated gravitational wave (GW) detection projects in the next decade, which will detect abundant compact binary systems. However, the precise prediction of space GW waveforms remains unexplored. To solve the data processing difficulty in the increasing waveform complexity caused by detectors' response and second-generation time-delay interferometry (TDI 2.0), an interpretable pre-trained large model named CBS-GPT (Compact Binary Systems Waveform Generation with Generative Pre-trained Transformer) is proposed. For compact binary system waveforms, three models were trained to predict the waveforms of massive black hole binary (MBHB), extreme mass-ratio inspirals (EMRIs), and galactic binary (GB), achieving prediction accuracies of 98%, 91%, and 99%, respectively. The CBS-GPT model exhibits notable interpretability, with its hidden parameters effectively capturing the intricate information of waveforms, even with complex instrument response and a wide parameter range. Our research demonstrates the potential of large pre-trained models in gravitational wave data processing, opening up new opportunities for future tasks such as gap completion, GW signal detection, and signal noise reduction.

The pursuit of understanding the mysteries surrounding dark energy has sparked significant interest within the field of cosmology. While conventional approaches, such as the cosmological constant, have been extensively explored, alternative theories incorporating scalar field-based models and modified gravity have emerged as intriguing avenues. Among these, teleparallel theories of gravity, specifically the $f(T,\phi)$ formulation, have gained prominence as a means to comprehend dark energy within the framework of teleparallelism. In this study, we investigate two well-studied models of teleparallel dark energy and examine the presence of cosmological singularities within these scenarios. Using the Goriely-Hyde procedure, we examine the dynamical systems governing the cosmological equations of these models. Our analysis reveals that both models exhibit Type IV singularities, but only for a limited range of initial conditions. These results could indicate a potential edge for teleparallel cosmological models over their other modified gravity counterparts, as the models we examine seem to be only allowing for weak singularities that too under non general conditions.

Background: The $^{12}{\rm C}(\alpha,\gamma)^{16}$O reaction, determining the survival of carbon in red giants, is of interest for nuclear reaction theory and nuclear astrophysics. Numerous attempts to obtain the astrophysical factor of the $^{12}{\rm C}(\alpha,\gamma)^{16}$O reaction, both experimental and theoretical, have been made for almost 50 years. The specifics of the $^{16}$O nuclear structure is the presence of two subthreshold bound states, (6.92 MeV, 2$^+$) and (7.12 MeV, 1$^-$), dominating the behavior of the low-energy $S$-factor. The strength of these subthreshold states is determined by their asymptotic normalization coefficients (ANCs) which need to be known with high accuracy. Recently, using the model-independent extrapolation method, Blokhintsev {\it et al.} [Eur. Phys. J. A {\bf 59}, 162 (2023)] determined the ANCs for the three subthreshold states in $^{16}$O. Purpose: In this paper, using these newly determined ANCs, we calculated the low-energy astrophysical $S$-factors for the $^{12}{\rm C}(\alpha,\gamma)^{16}$O radiative capture. Method: The $S$-factors are calculated within the framework of the $R$-matrix method using the AZURE2 code. Conclusion: Our total $S$-factor includes the resonance $E1$ and $E2$ transitions to the ground state of $^{16}$O interfering with the corresponding direct captures and cascade radiative captures to the ground state of $^{16}$O through four subthreshold states: $0_2^+,\,3^-,\, 2^+$ and $1^-$. Since our ANCs are higher than those used by deBoer {\it et al.} [Rev. Mod. Phys. {\bf 89}, 035007 (2017)], the present total $S$-factor at the most effective astrophysical energy of 300 keV is 174 keVb versus 137 keVb of that work. Accordingly, our calculated reaction rate at low temperatures ($T_{9} < 2$) is higher than the one given in the aforesaid paper.

We propose the possibility that compact extra dimensions can obtain large size by higher dimensional inflation, relating the weakness of the actual gravitational force to the size of the observable universe. Solution to the horizon problem implies that the fundamental scale of gravity is smaller than $10^{13}$ GeV which can be realised in a braneworld framework for any number of extra dimensions. However, requirement of (approximate) flat power spectrum of primordial density fluctuations consistent with present observations makes this simple proposal possible only for one extra dimension at around the micron scale. After the end of five-dimensional inflation, the radion modulus can be stabilised at a vacuum with positive energy of the order of the present dark energy scale. An attractive possibility is based on the contribution to the Casimir energy of right-handed neutrinos with a mass at a similar scale.

The cosmological evolution within the framework of exponential $F(R)$ gravity is analysed by assuming two forms for dark matter: (a) a standard dust-like fluid and (b) an axion scalar field. As shown in previous literature, an axion-like field oscillates during the cosmological evolution but can play the role of dark matter when approaching the minimum of its potential. Both scenarios are confronted with recent observational data including the Pantheon Type Ia supernovae, Hubble parameter estimations (Cosmic Chronometers), Baryon Acoustic Oscillations and Cosmic Microwave Background distances. The models show great possibilities in describing these observations when compared with the $\Lambda$CDM model, supporting the viability of exponential $F(R)$ gravity. The differences between both descriptions of dark matter is analysed.

We assess the impact of accurate, self-consistent modelling of thermal effects in neutron-star merger remnants in the context of third-generation gravitational-wave detectors. This is done through the usage, in Bayesian model selection experiments, of numerical-relativity simulations of binary neutron star (BNS) mergers modelled through: a) nuclear, finite-temperature (or ``tabulated'') equations of state (EoSs), and b) their simplifed piecewise (or ``hybrid'') representation. These cover four different EoSs, namely SLy4, DD2, HShen and LS220. Our analyses make direct use of the Newman-Penrose scalar $\psi_4$ outputted by numerical simulations. Considering a detector network formed by three Cosmic Explorers, we show that differences in the gravitational-wave emission predicted by the two models are detectable with a natural logarithmic Bayes Factor $\log{\cal{B}}\geq 5$ at average distances of $d_L \simeq 50$Mpc, reaching $d_L \simeq 100$Mpc for source inclinations $\iota \leq 0.8$, regardless of the EoS. This impact is most pronounced for the HShen EoS. For low inclinations, only the DD2 EoS prevents the detectability of such modelling differences at $d_L \simeq 150$Mpc. Our results suggest that the usage a self-consistent treatment of thermal effects is crucial for third-generation gravitational wave detectors.

The axion is a long-postulated boson that can simultaneously solve two fundamental problems of modern physics: the charge-parity symmetry problem in the strong interaction and the enigma of dark matter. In this work we estimate, by means of Monte Carlo simulations, the sensitivity of the Dark-photons$\&$Axion-Like particles Interferometer (DALI), a new-generation Fabry-P\'erot haloscope proposed to probe axion dark matter in the 25-250 $\mu$eV band.

The nature of the dark matter in the Universe is one of the hardest unsolved problems in modern physics. Indeed, on one hand, the overwhelming indirect evidence from astrophysics seems to leave no doubt about its existence; on the other hand, direct search experiments, especially those conducted with low background detectors in underground laboratories all over the world seem to deliver only null results, with a few debated exceptions. Furthermore, the lack of predicted candidates at the LHC energy scale has made this dichotomy even more puzzling. We will recall the most important phases of this novel branch of experimental astro-particle physics, analyzing the interconnections among the main projects involved in this challenging quest, and we will draw conclusions slightly different from how the problem is commonly understood.

We study the possibility that dark matter re-enters kinetic equilibrium with a radiation bath after kinetic decoupling, a scenario we dub kinetic recoupling. This naturally occurs, for instance, with certain types of resonantly-enhanced interactions, or as the result of a phase transition. While late kinetic decoupling damps structure on small scales below a cutoff, kinetic recoupling produces more complex changes in the power spectrum that depend on the nature and extent of the recoupling period. We explore the features that kinetic recoupling imprints upon the matter power spectrum, and discuss how such features can be traced to dark matter microphysics with future observations.

A recent research presented a solution for a Schwarzschild black hole with a force-free magnetic field showing how the canonical form of background metric changes, Found. Phys. 52 (2022) 4, 93. Therefore, it is logical that this research seeks modified solutions for the Kerr black hole as well. The goal of this paper is, thereby, to pinpoint an exact solution for the Kerr black holes perturbed by force-free magnetic fields. However, to do this we use a well-known tetrad formalism and obtain an explicit expression for the electromagnetic strength tensor in the background metric, that is Kerr, based on the tangent space calculations. Analyzing the stress-energy tensor reveals the perturbed factors in the metric that allow us to understand better the physics of the disks around massive and supermassive black holes. These results show that it is not enough to depend solely on perturbation theory set against the background metric, without considering the backreaction of the force-free magnetufluid on it. This research indicates the necessity of studying the impact of force-free magnetic fields on the structure and evolution of mysterious cosmic objects.

In this Reply we include the corrections suggested in the Comment [Phys. Rev. Lett. 131, 169001]. We show that their impact on our results is small, and that the overall conclusion of the Article [Phys. Rev. Lett. 129, 111102] are robust. As pointed out in the Article, it is crucial to account for the statistical uncertainty in the ringdown starting time, neglected in most previous studies. This uncertainty is ~40 times larger than the systematic shift induced by the software bug mentioned in the Comment. The remaining discrepancies between the Comment and the Article can be attributed to additional differences in the setup, notably the sampling rate and the noise estimation method (in the Article the latter was chosen to mimic the original methods of [Phys. Rev. Lett. 123, 111102]). Beyond data analysis considerations, the physics of the problem cannot be ignored. As shown in [arXiv:2302.03050], a model consisting of a sum of constant-amplitude overtones starting at the peak of the waveform introduces uncontrolled systematic uncertainties in the measurement due to dynamical and strong-field effects. These theoretical considerations imply that studies based on such models cannot be interpreted as black hole spectroscopy tests.