Papers for Monday, Dec 23 2024

Stellar flares are powerful bursts of electromagnetic radiation triggered by magnetic reconnection in the chromosphere of stars, occurring frequently and intensely on active M dwarfs. While missions like TESS and Kepler have studied regular and super-flares, their detection of flares with energies below $10^{30}$ erg remains incomplete. Extending flare studies to include these low-energy events could enhance flare formation models and provide insight into their impacts on exoplanetary atmospheres. This study investigates CHEOPS's capacity to detect low-energy flares in M dwarf light curves. Using CHEOPS's high photometric precision and observing cadence, along with a tailored wavelet-based denoising algorithm, we aim to improve detection completeness and refine flare statistics for low-energy events. We conducted a flare injection and recovery process to optimise denoising parameters, applied it to CHEOPS light curves to maximise detection rates, and used a flare breakdown algorithm to analyse complex structures. Our analysis recovered 349 flares with energies ranging from $2.2\times10^{26}$ to $8.1\times10^{30}$ erg across 63 M dwarfs, with $\sim$40% exhibiting complex, multi-peaked structures. The denoising algorithm improved flare recovery by $\sim$34%, though it marginally extended the lower boundary of detectable energies. For the full sample, the power-law index $\alpha$ was $1.92\pm0.07$, but a log-normal distribution fit better, suggesting multiple flare formation scenarios. While CHEOPS's observing mode is not ideal for large-scale surveys, it captures weaker flares than TESS or Kepler, expanding the observed energy range. Wavelet-based denoising enhances low-energy event recovery, enabling exploration of the micro-flaring regime. Expanding low-energy flare observations could refine flare generation models and improve the understanding of their role in star-planet interactions.

Gas cooling and heating rates are vital components of hydrodynamic simulations. However, they are computationally expensive to evaluate exactly with chemical networks or photoionization codes. We compare two different approximation schemes for gas cooling and heating in an idealized simulation of an isolated galaxy. One approximation is based on a polynomial interpolation of a table of Cloudy calculations, as is commonly done in galaxy formation simulations. The other approximation scheme uses machine learning for the interpolation instead on an analytic function, with improved accuracy. We compare the temperature-density phase diagrams of gas from each simulation run to assess how much the two simulation runs differ. Gas in the simulation using the machine learning approximation is systematically hotter for low-density gas with $-3 \lesssim \log{(n_b/\mathrm{cm}^{-3})} \lesssim -1$. We find a critical curve in the phase diagram where the two simulations have equal amounts of gas. The phase diagrams differ most strongly at temperatures just above and below this critical curve. We compare CII emission rates for collisions with various particles (integrated over the gas distribution function), and find slight differences between the two simulations. Future comparisons with simulations including radiative transfer will be necessary to compare observable quantities like the total CII luminosity.

We present observations of ASASSN-22ci (AT2022dbl), a nearby tidal disruption event (TDE) discovered by the All-Sky Automated Survey for Supernovae (ASAS-SN) at a distance of d$_L \simeq 125$ Mpc. Roughly two years after the initial ASAS-SN discovery, a second flare was detected coincident with ASASSN-22ci. UV/optical photometry and optical spectroscopy indicate that both flares are likely powered by TDEs. The striking similarity in flare properties suggests that these flares result from subsequent disruptions of the same star. Each flare rises on a timescale of $\sim$30 days, has a temperature of $\approx$30,000 K, a peak bolometric luminosity of $L_{UV/Opt} = 10^{43.6 - 43.9} \textrm{ erg} \textrm{ s}^{-1}$, and exhibits a blue optical spectrum with broad H, He, and N lines. No X-ray emission is detected during either flare, but X-ray emission with an unabsorbed luminosity of $L_{X} = 3\times10^{41} \textrm{ erg} \textrm{ s}^{-1}$ and $kT = 0.042$ eV is observed between the flares. Pre-discovery survey observations rule out the existence of earlier flares within the past $\approx$6000 days, indicating that the discovery of ASASSN-22ci likely coincides with the first flare. If the observed flare separation of $720 \pm 4.7$ days is the orbital period, the next flare of ASASSN-22ci should occur near MJD 61075 (2026 February 04). Finally, we find that the existing sample of repeating TDE candidates is consistent with Hills capture of a star initially in a binary with a total mass between $\sim$$1 - 4$ M$_{\odot}$ and a separation of $\sim$$0.01 - 0.1$ AU.

Detection and study of the intracluster light in rich clusters of galaxies has been a problem of long standing challenge and interest. Using the lowest surface brightness images of the Coma cluster of galaxies in the g and r bands, from the Halos and Environment of Nearby Galaxies (HERON) Coma Cluster Project, we obtained the most extensive image of intracluster light (ICL) in a single cluster to date, spreading over 1.5 Mpc from the cluster core. The unprecedented wealth of spectroscopic data made publicly available by the Dark Energy Spectroscopic Instrument (DESI) Early Data Release, complemented with a compilation from the NASA/IPAC Extragalactic Database and the literature, enabled the identification of 2,157 galaxy members within Coma, from which 42 distinct groups were identified. The synergy between these high-quality data allowed us to: 1) calculate ICL fractions of $19.9\pm0.5$\% and $19.6\pm0.6$\% in the g and r bands, respectively, consistent with a dynamically active cluster, 2) unveil Coma's faintest tidal features, and 3) provide a comprehensive picture of the dynamics and interactions within this complex system. Our findings indicate that the ICL connects several of these groups in a filamentous network, from which we infer the ongoing dynamical processes. In particular, we identified a faint stellar bridge linking the core of Coma with the galaxy NGC 4839, providing compelling evidence that this galaxy has already traversed the central region of the cluster.

Supermassive black holes are prevalent at the centers of massive galaxies, and their masses scale with galaxy properties, increasing evidence suggesting that these trends continue to low stellar masses. Seeds are needed for supermassive black holes, especially at the highest redshifts explored by the James Webb Space Telescope. We study the hierarchical merging of galaxies via cosmological merger trees and argue that the seeds of supermassive black holes formed in nuclear star clusters via stellar black hole mergers at early epochs. Observable tracers include intermediate-mass black holes, nuclear star clusters, and early gas accretion in host dwarf galaxies, along with a potentially detectable stochastic gravitational wave background, ejection of intermediate and supermassive black holes, and consequences of a significant population of tidal disruption events and extreme-mass ratio inspirals.

In this paper, we study how absorption-line systems affect the spectra and redshifts of quasars (QSOs), using catalogs of Mg II absorbers from the early data release (EDR) and first data release (DR1) of the Dark Energy Spectroscopic Instrument (DESI). We determine the reddening effect of an absorption system by fitting an un-reddened template spectrum to a sample of 50,674 QSO spectra that contain Mg II absorbers. We find that reddening caused by intervening absorbers (voff > 3500 km/s) has an average color excess of E(B-V) = 0.04 magnitudes. We find that the E(B-V) tends to be greater for absorbers at low redshifts, or those having Mg II absorption lines with higher equivalent widths, but shows no clear trend with voff for intervening systems. However, the E(B-V) of associated absorbers, those at voff < 3500 km/s, shows a strong trend with voff , increasing rapidly with decreasing voff and peaking (approximately 0.15 magnitudes) around voff = 0 km/s. We demonstrate that Mg II absorbers impact redshift estimation for QSOs by investigating the distributions of voff for associated absorbers. We find that at z > 1.5 these distributions broaden and bifurcate in a nonphysical manner. In an effort to mitigate this effect, we mask pixels associated with the Mg II absorption lines and recalculate the QSO redshifts. We find that we can recover voff populations in better agreement with those for z < 1.5 absorbers and in doing so typically shift background QSO redshifts by delta_z approximately equal to plus or minus 0.005.

Building upon [2308.02636], this article investigates the potential constraining power of persistent homology for cosmological parameters and primordial non-Gaussianity amplitudes in a likelihood-free inference pipeline. We evaluate the ability of persistence images (PIs) to infer parameters, compared to the combined Power Spectrum and Bispectrum (PS/BS), and we compare two types of models: neural-based, and tree-based. PIs consistently lead to better predictions compared to the combined PS/BS when the parameters can be constrained (i.e., for $\{\Omega_{\rm m}, \sigma_8, n_{\rm s}, f_{\rm NL}^{\rm loc}\}$). PIs perform particularly well for $f_{\rm NL}^{\rm loc}$, showing the promise of persistent homology in constraining primordial non-Gaussianity. Our results show that combining PIs with PS/BS provides only marginal gains, indicating that the PS/BS contains little extra or complementary information to the PIs. Finally, we provide a visualization of the most important topological features for $f_{\rm NL}^{\rm loc}$ and for $\Omega_{\rm m}$. This reveals that clusters and voids (0-cycles and 2-cycles) are most informative for $\Omega_{\rm m}$, while $f_{\rm NL}^{\rm loc}$ uses the filaments (1-cycles) in addition to the other two types of topological features.

The long-duration balloon platform for radio detection of energetic neutrinos, pioneered by ANITA, affords large instantaneous effective areas but has limited livetime. Conversely, tethered balloons, traditionally used for radio detection of non-science objectives (such as electronic surveillance), allow for much longer livetimes, albeit at significantly lower altitudes. In this contribution, a tethered balloon platform for neutrino detection is considered, including estimates of the neutrino sensitivity and a discussion the feasibility of such a platform. Both the Askaryan and tau-neutrino induced Extensive Air Shower channels are explored as target science data for the platform. Tethered balloons locations on ice sheets, land, and in steep valleys or fjords will be considered. Ground-based array calibration and cosmic-ray air shower use cases will also be briefly commented on.

In models of dark matter composed of feebly interacting ultralight bosons in the state of Bose-Einstein condensate, the dynamical friction force acting on circularly moving globular clusters modelled as Plummer spheres is determined. Analytic expressions for both radial and tangential components of the dynamical friction force are given. We reveal that the dynamical friction force for the Plummer sphere deviates from that for a point probe of the same mass for the significantly large ratio of the Plummer sphere radius to its orbital radius as well as for large values of the Mach number.

We report the detection of the [O III] auroral line in 42 galaxies within the redshift range of $3 < z < 10$. These galaxies were selected from publicly available JWST data releases, including the JADES and PRIMALsurveys, and observed using both the low-resolution PRISM/CLEAR configuration and medium-resolution gratings. The measured electron temperatures in the high-ionization regions of these galaxies range from $T_e$([O III]) = 12,000 to 24,000 K, consistent with temperatures observed in local metal-poor galaxies and previous JWST studies. In 10 galaxies, we also detect the [O II] auroral line, allowing us to determine electron temperatures in the low-ionization regions, which range between $T_e$([O II]) = 10,830 and 20,000 K. The direct-$T_e$-based metallicities of our sample span from 12 + log(O/H) = 7.2 to 8.4, indicating these high-redshift galaxies are relatively metal-poor. By combining our sample with 25 galaxies from the literature, we expand the dataset to a total of 67 galaxies within $3 < z < 10$, effectively more than doubling the previous sample size for direct-$T_e$ based metallicity studies. This larger dataset allow us to derive empirical metallicity calibration relations based exclusively on high-redshift galaxies, using six key line ratios: R3, R2, R23, Ne3O2, O32, and O3N2. Notably, we derive a novel metallicity calibration relation for the first time using high-redshift $T_e$-based metallicities: $\hat{R}$ = 0.18log $R2$ + 0.98log $R3$. This new calibration significantly reduces the scatter in high-redshift galaxies compared to the $\hat{R}$ relation previously calibrated for low-redshift galaxies.

Binary neutron stars (BNSs) detected in the Milky Way have the total masses distributing narrowly around $\sim2.6-2.7M_\odot$, while the BNS merger GW190425 detected via gravitational wave has a significantly larger mass ($\sim3.4M_\odot$). This difference is not well understood, yet. In this paper, we investigate the BNS spin evolution via an improved binary star evolution model and its effects on the BNS observability, with implementation of various relevant astrophysical processes. We find that the first-born neutron star component in low-mass BNSs can be spun up to millisecond pulsars by the accretion of Roche-lobe overflow from its companion and its radio lifetime can be comparable to the Hubble time. However, most high-mass BNSs have substantially shorter radio lifetime than the low-mass BNSs, and thus smaller probability being detected via radio emission. Adopting the star formation and metal enrichment history of the Milky Way given by observations, we obtain the survived Galactic BNSs with pulsar components from our population synthesis model and find that their distributions on the diagrams of spin period versus spin-period-time-derivative ($P-\dot{P}$) and orbital period versus eccentricity ($P_{\rm orb}-e$) can well match those of the observed Galactic BNSs. The total mass distribution of the observed Galactic BNSs can also be matched by the model. A significant fraction ($\sim19\%-22\%$) of merging BNSs at redshift $z\sim0$ have masses $\gtrsim3M_\odot$, which seems compatible with the GW observations. Future radio observations may detect many more Galactic BNSs, which will put strong constraint on the spin evolution of BNSs during their formation processes.

We simulate the formation of Fuzzy Dark Matter (FDM) cores in the presence of a Black Hole (BH) to explore whether BHs can serve as seeds for FDM core condensation. Our analysis is based on the core-condensation via the kinetic relaxation process for random initial conditions of the FDM. We show that in general the BH merges with pre-collapsed mini-clusters and once they share location the BH oscillates within the core. The condensation takes place around the black hole and the FDM acquires a density profile consistent with the density of the stationary solution of the FDM+BH eigenvalue problem in average. The central density of the resulting core depends on the mass of the BH, which due to its motion relative to the FDM cloud produces a smaller time-averaged densities for bigger BH masses, which lead to a new diversity of central FDM core densities. Our results indicate that BHs can indeed act as focal points for FDM core condensation. As a collateral result, for our analysis we revised the construction of stationary solutions of FDM+BH and found a phenomenological formula for the FDM density that can be used to fit FDM cores around BHs.

In this manuscript, optical quasi-periodic oscillations (QPOs) with 550 day periodicity related to a candidate of sub-pc binary black hole (BBH) system are reported in the reverberation mapped broad line quasar PG 1411+442 but with different line profile of broad H$\alpha$ from that of broad H$\beta$ in its rms spectrum. First, considering sine function to describe the 18.8years-long light curves from the CSS, ASAS-SN and ZTF, 550days periodicity can be confirmed with confidence level higher than 5$\sigma$. Second, the stable 550days optical QPOs can be re-confirmed with confidence levels higher than 5$\sigma$ by the Generalized Lomb-Scargle periodogram, the sine-like phase-folded light curves and the WWZ technique determined power maps. Third, based on simulated light curves by CAR process, confidence level higher than $3.5\sigma$ can be confirmed for the optical QPOs not related to intrinsic AGN variability. Moreover, considering spatial separation of central two BH accreting systems smaller than expected sizes of broad emission line regions (BLRs), central total BH mass higher than $10^6{\rm M_\odot}$ could lead to few effects of supposed BBH systems on estimated virial BH masses. Meanwhile, disk precession is not preferred due to the similar estimated sizes of optical and NUV emission regions, and jet precession can be ruled out due to PG 1411+442 as a radio quiet quasar. The results strongly indicate it would be practicable by applying very different line profiles of broad Balmer emission lines to detect candidates of BBH systems in normal broad line AGN in the near future.

Binaries harboring a millisecond pulsar (MSP) and a black hole (BH) are a key observing target for current and upcoming pulsar surveys. We model the formation and evolution of such binaries in isolation at solar metallicity using the next-generation binary population synthesis code POSYDON. We examine neutron star (NS)-BH binaries where the NS forms first (labeled NSBH), as the NS must be able to spin-up to MSP rotation periods before the BH forms in these systems. We find that NSBHs are very rare and have a birth rate < 1 Myr$^{-1}$ for a Milky Way-like galaxy in our typical models. The NSBH birth rate is 2-3 orders of magnitude smaller than that for NS-BHs where the BH forms first (labeled BHNS). These rates are also sensitive to model assumptions about the supernova (SN) remnant masses, natal kicks, and common-envelope efficiency. We find that 100% of NSBHs undergo a mass ratio reversal before the first SN and up to 64% of NSBHs undergo a double common envelope phase after the mass ratio reversal occurs. Most importantly, no NSBH binaries in our populations undergo a mass transfer phase, either stable or unstable, after the first SN. This implies that there is no possibility of pulsar spin-up via accretion, and thus MSP-BH binaries cannot form. Thus, dynamical environments and processes may provide the only formation channels for such MSP-BH binaries.

The DESI Legacy Imaging Surveys (DESI-LIS) comprise three distinct surveys: the Dark Energy Camera Legacy Survey (DECaLS), the Beijing-Arizona Sky Survey (BASS), and the Mayall z-band Legacy Survey (MzLS).The citizen science project Galaxy Zoo DECaLS 5 (GZD-5) has provided extensive and detailed morphology labels for a sample of 253,287 galaxies within the DECaLS survey. This dataset has been foundational for numerous deep learning-based galaxy morphology classification studies. However, due to differences in signal-to-noise ratios and resolutions between the DECaLS images and those from BASS and MzLS (collectively referred to as BMz), a neural network trained on DECaLS images cannot be directly applied to BMz images due to distributional this http URL this study, we explore an unsupervised domain adaptation (UDA) method that fine-tunes a source domain model trained on DECaLS images with GZD-5 labels to BMz images, aiming to reduce bias in galaxy morphology classification within the BMz survey. Our source domain model, used as a starting point for UDA, achieves performance on the DECaLS galaxies' validation set comparable to the results of related works. For BMz galaxies, the fine-tuned target domain model significantly improves performance compared to the direct application of the source domain model, reaching a level comparable to that of the source domain. We also release a catalogue of detailed morphology classifications for 248,088 galaxies within the BMz survey, accompanied by usage recommendations.

In this brief note, we tentatively investigate the possibility that the radio filaments are produced when the Galactic center wind washes over magnetic field structures. The electrons and ions, with their disparate charge-to-mass ratios, are deflected differently by the magnetic field, and a current results. The current is subsequently Z-pinched into filaments, creating an electron-accelerating electric field along the way, because the magnetic field necessarily rearranges during the dynamic constriction process. An axial magnetic field also arises, possibly via the diocotron channel, to eventually quench the pinching and stabilize the filaments against a variety of instabilities.

Most previous efforts for hydrodynamic studies on detonation in the context of Type Ia supernovae did not take into account the scale of the cellular structure for a criterion in initiation, propagation, quenching, and the resolution requirement of detonation, whereas it is quite common to consider cell sizes in the discussion on terrestrial detonation in chemically reactive systems. In our recent study, the terrestrial cell-based theories, which incorporates the cell-size data acquired in 2D simulations of helium detonation in the double-detonation model, were demonstrated to be a powerful diagnostics in reproducing the thresholds in the initiation and quenching provided by previous studies. In the present study, 2D simulation results of the cellular detonation in the base of white-dwarf (WD) envelope are described in detail, in terms of the dynamic wave morphology and chemical abundance structure. The cellular structure is observed at a range of upstream density and envelope composition explored in the present work. C/O contamination by the WD core material reduces the cell width rapidly, as accelerated by the {\alpha}-capture reaction. It is also indicated that nickel production could be significantly delayed for the C/O-rich composition. The small cell width makes it extremely demanding to resolve the detonation structure in full-star simulations of SNe Ia; this could raise a concern on the robustness of the outcomes of some numerical simulations in terms of the success and failure of detonation. This issue may be overcome by sub-grid modeling that incorporates the cellular dynamics acquired in resolved simulations.

The Cosmological Gravitational Wave Background (CGWB) anisotropies contain valuable information about the physics of the early universe. Given that General Relativity is intrinsically nonlinear, it is important to look beyond first-order contributions in cosmological perturbations. In this work, we present a non-perturbative approach for the computation of CGWB anisotropies at large scales, providing the extension of the initial conditions and the Sachs-Wolfe effect for the CGWB, which encodes the full non-linearity of the scalar metric perturbations. We also derive the non-perturbative expression for three-point correlation of the gravitational wave energy density perturbation in the case of an inflationary CGWB with a scale-invariant power spectrum and negligible primordial non-Gaussianity. We show that, under such conditions, the gravitational wave energy density perturbations are lognormally distributed, leading to an interesting effect such as intermittency.

We have studied the accreting black hole binary GX 339-4 using two highly accurate broad-band X-ray data sets in very soft spectral states from simultaneous NICER and NuSTAR observations. Simultaneous fitting of both data sets with relativistic models of the disk, its Comptonization and reflection allows us to relatively accurately determine the black-hole mass and spin, and the distance and inclination. However, we find the measured values strongly depend on the used disk model. With the widely used thin-disk Kerr models kerrbb and kerrbb2 (which employ color corrections), we find relatively low masses and strongly negative spins. Then, the models utilizing detailed disk atmospheric spectra, bhspec and slimbh, predict moderately positive spins and high masses. When adding a warm corona above the disk (as proposed before for both AGNs and accreting binaries), we find the spin is weakly constrained, but consistent with zero. In all cases, the fitted inclination is low, $\approx$30-$34^\circ$. For the spin aligned with the binary orbit, the mass function for this binary implies large values of the mass, consistent only with those obtained with either slimbh or warm corona. We also test the disk models for an assumed set of mass, distance and inclination. We find that, e.g., kerrbb yields values of the spin parameter lower than bhspec or slimbh by $\sim$0.2-0.3. Our results confirm previously found strong disk-model dependencies of the measured black-hole spin, now for a low-mass X-ray binary.

We test and compare coasting cosmological models with curvature parameters ${k=\left\{ -1,0,+1 \right\}}$ in ${H_0^2 c^{-2}}$ units and the flat $\Lambda$CDM model by fitting them to cosmic chronometers (CC), the Pantheon+ sample of type Ia supernovae (SNe), and standardized quasars (QSOs). We used the \texttt{emcee} code for fitting CC data, a custom Markov Chain Monte Carlo implementation for SNe and QSOs, and Anderson-Darling tests for normality on normalized residuals for model comparison. Best-fit parameters are presented, constrained by data within redshift ranges $z\leq 2$ for CCs, $z\leq 2.3$ for SNe, and $z\leq 7.54$ for QSOs. Coasting models, particularly the flat coasting model, are generally favored over the flat $\Lambda$CDM model. The overfitting of the flat $\Lambda$CDM model to Pantheon+ SNe and the large intrinsic scatter in QSO data suggest a need to refine error estimates in these datasets. We also highlight the seemingly fine-tuned nature of either the CC data or $\Omega_{\mathrm{m},0}$ in the flat $\Lambda$CDM model to an ${H_1=H_0}$ coincidence when fitting ${H(z)=H_1z+H_0}$, a natural feature of coasting models.

A complete sample of red supergiant stars (RSGs) is important for studying their properties. Identifying RSGs in extragalatic field first requires removing the Galactic foreground dwarfs. The color-color diagram (CCD) method, specifically using $r-z/z-H$ and $J-H/H-K$, has proven successful in several studies. However, in metal-poor galaxies, faint RSGs will mix into the dwarf branch in the CCD and would be removed, leading to an incomplete RSG sample. This work attempts to improve the CCD method in combination with the Gaia astrometric measurement to remove foreground contamination in order to construct a complete RSG sample in metal-poor galaxies. The empirical regions of RSGs in both CCDs are defined and modified by fitting the locations of RSGs in galaxies with a range of metallicity. The metal-poor galaxy NGC 6822 is taken as a case study for its low metallicity ([Fe/H] $\approx$ -1.0) and moderate distance (about 500 kpc). In the complete sample, we identify 1,184 RSG, 1,559 oxygen-rich AGB (O-AGBs), 1,075 carbon-rich AGB (C-AGBs), and 140 extreme AGB (x-AGBs) candidates, with a contamination rate of approximately 20.5%, 9.7%, 6.8%, and 5.0%, respectively. We also present a pure sample, containing only the sources away from the dwarf branch, which includes 843 RSG, 1,519 O-AGB, 1,059 C-AGB, and 140 x-AGB candidates, with a contamination rate of approximately 6.5%, 8.8%, 6.1%, and 5.0%, respectively. About 600 and 450 RSG candidates are newly identified in the complete and pure sample, respectively, compared to the previous RSG sample in NGC 6822.

A gravitational wave (GW) passing through an astrometric observer causes periodic shifts of the apparent star positions measured by the observer. For a GW of sufficient amplitude and duration, and of suitable frequency, these shifts might be detected with a Gaia-like astrometric telescope. This paper aims to analyse in detail the effects of GWs on an astrometric solution based on Gaia-like observations, which are one-dimensional, strictly differential between two widely separated fields of view and following a prescribed scanning law. We present a simple geometric model for the astrometric effects of a plane GW in terms of the time-dependent positional shifts. Using this model, the general interaction between the GW and a Gaia-like observation is discussed. Numerous Gaia-like astrometric solutions are made, taking as input simulated observations that include the effects of a continuous plain GW with constant parameters and periods ranging from ~50 days to 100 years. The resulting solutions are analysed in terms of the systematic errors on astrometric and attitude parameters, as well as the observational residuals. It is found that a significant part of the GW signal is absorbed by the astrometric parameters, leading to astrometric errors of a magnitude (in radians) comparable to the strain parameters. These astrometric errors are in general not possible to detect, because the true (unperturbed) astrometric parameters are not known to corresponding accuracy. The astrometric errors are especially large for specific GW frequencies that are linear combinations of two characteristic frequencies of the scanning law. Nevertheless, for all GW periods smaller than the time span covered by the observations, significant parts of the GW signal also go into the astrometric residuals. This fosters the hope for a GW detection algorithm based on the residuals of standard astrometric solutions.

The availability of accurate and timely state predictions for objects in near-Earth orbits is becoming increasingly important due to the growing congestion in key orbital regimes. The Two-line Element Set (TLE) catalogue remains, to this day, one of the few publicly-available, comprehensive sources of near-Earth object ephemerides. At the same time, TLEs are affected by measurement noise and are limited by the low accuracy of the SGP4 theory, introducing significant uncertainty into state predictions. Previous literature has shown that filtering TLEs with batch least squares methods can yield significant improvements in long-term state prediction accuracy. However, this process can be highly sensitive to TLE quality which can vary throughout the year. In this study, it is shown that either extended-duration fit windows of the order of months, or the removal of systematic biases in along-track position prior to state estimation can produce significant reductions in post-fit position errors. Simple models for estimating these systematic biases are shown to be effective without introducing the need for high-complexity Machine Learning (ML) models. Furthermore, by establishing a TLE-based error metric, the need for high accuracy ephemerides is removed when creating these models. For selected satellites in the Medium Earth Orbit (MEO) regime, post-fit position errors are reduced by up to 80 %, from approximately 5 km to 1 km; meanwhile, for selected satellites in the Geostationary Earth Orbit (GEO)/Geosynchronous Earth Orbit (GSO) regime, large oscillations in post-fit position error can be suppressed.

The Imaging X-ray Polarimetry Explorer (IXPE) observed for the first time highly polarized X-ray emission from the magnetar 1E 1841-045, targeted after a burst-active phase in August 2024. To date, IXPE has observed four other magnetars during quiescent periods, highlighting substantially different polarization properties. 1E 1841-045 exhibits a high, energy-dependent polarization degree, which increases monotonically from ~15% at 2-3 keV up to ~55% at 5.5-8 keV, while the polarization angle, aligned with the celestial North, remains fairly constant. The broadband spectrum (2-79 keV) obtained by combining simultaneous IXPE and NuSTAR data is well modeled by a blackbody and two power-law components. The unabsorbed 2-8 keV flux (~2E-11 erg/cm2/s) is about 10% higher than that obtained from archival XMM-Newton and NuSTAR observations. The polarization of the soft, thermal component does not exceed ~25%, and may be produced by a condensed surface or a bombarded atmosphere. The intermediate power law is polarized at around 30%, consistent with predictions for resonant Compton scattering in the star magnetosphere; while, the hard power law exhibits a polarization degree exceeding 65%, pointing to a synchrotron/curvature origin.

The BL Lacertae object VER J0521+211 underwent a notable flaring episode in February 2020. A short-term monitoring campaign, led by the MAGIC (Major Atmospheric Gamma Imaging Cherenkov) collaboration, covering a wide energy range from radio to very-high-energy (VHE, 100 GeV < E < 100 TeV) gamma rays was organised to study its evolution. These observations resulted in a consistent detection of the source over six consecutive nights in the VHE gamma-ray domain. Combining these nightly observations with an extensive set of multiwavelength data made modelling of the blazar's spectral energy distribution (SED) possible during the flare. This modelling was performed with a focus on two plausible emission mechanisms: i) a leptonic two-zone synchrotron-self-Compton scenario, and ii) a lepto-hadronic one-zone scenario. Both models effectively replicated the observed SED from radio to the VHE gamma-ray band. Furthermore, by introducing a set of evolving parameters, both models were successful in reproducing the evolution of the fluxes measured in different bands throughout the observing campaign. Notably, the lepto-hadronic model predicts enhanced photon and neutrino fluxes at ultra-high energies (E > 100 TeV). While the photon component, generated via decay of neutral pions, is not directly observable as it is subject to intense pair production (and therefore extinction) through interactions with the cosmic microwave background photons, neutrino detectors (e.g. IceCube) can probe the predicted neutrino component. Finally, the analysis of the gamma-ray spectra, as observed by MAGIC and the Fermi-LAT telescopes, yielded a conservative 95\% confidence upper limit of z \leq 0.244 for the redshift of this blazar.

We explore the impact of the reionization history on examining the shape of the power spectrum of the primordial gravitational waves (PGWs) with the cosmic microwave background (CMB) polarization. The large-scale CMB generated from the reionization epoch is important in probing the PGWs from all-sky experiments, such as LiteBIRD. The reionization model has been constrained by several astrophysical observations. However, its uncertainty could impact constraining models of the PGWs if we use large-scale CMB polarization. Here, by expanding the analysis of Mortonson & Hu (2007), we estimate how reionization uncertainty impacts constraints on a generic primordial tensor power spectrum. We assume that CMB polarization is measured by a LiteBIRD-like experiment and the tanh model is adopted for a theoretical template when we fit data. We show that constraints are almost unchanged even if the true reionization history is described by an exponential model, where all parameters are within 68% Confidence Level (CL). We also show an example of the reionization history that the constraints on the PGWs are biased more than 68% CL. Even in that case, using E-mode power spectrum on large scales would exclude such a scenario and make the PGW constraints robust against the reionization uncertainties.

Context. The Galactic ASKAP collaboration (GASKAP) is undertaking an HI emission survey of the 21cm line to map the Magellanic system and the Galactic plane with the Australian Square Kilometre Array Pathfinder (ASKAP). One of the first areas observed in the Pilot Phase I of the survey was the Small Magellanic Cloud (SMC). Previous surveys of the SMC have uncovered new structures in the periphery of the SMC, along relatively low column density lines of sight. Aims. In this work we aimed to uncover the phase distribution of three distinct structures in the periphery of the SMC. This work will add to the constraints we have on the existence and survival of the cold neutral medium (CNM) in the SMC. Methods. We used ROHSA, a Gaussian decomposition algorithm, to model the emission across each cloud and classify the HI emission into their respective phases based on the linewidths of the fitted Gaussians. We created maps of velocity and column density of each phase of the HI across these three clouds. We measured the HI mass and CNM number density for each cloud. We also compared the HI results across the different phases with other gas tracers. Results. We find that in two clouds, the ends of each cloud are almost completely CNM dominated. Analysis of these two clouds indicates they are experiencing a compressive force from the direction of the SMC main body. In the third cloud we find a uniform CNM distribution along one wall of what is likely a supershell structure. Comparison with previous measurements of CO clumps in two of the clouds show the CO and HI are co-moving within a few km/s in regions of high HI column density, particularly when considering just the CNM.

We present this http URL, an efficient, modular, easy-to-use Bayesian pulsar timing and noise analysis package written in Julia. this http URL provides an independent, efficient, and parallelized implementation of the full non-linear pulsar timing and noise model along with a Python binding named pyvela. One-time operations such as data file input, clock corrections, and solar system ephemeris computations are performed by pyvela with the help of the PINT pulsar timing package. Its reliability is ensured via careful design utilizing Julia's type system, strict version control, and an exhaustive test suite. This paper describes the design and usage of this http URL focusing on the narrowband paradigm.

Deriving physical parameters from integrated galaxy spectra is paramount to interpret the cosmic evolution of star formation, chemical enrichment, and energetic sources. We develop modeling techniques to characterize the ionized gas properties in the subset of 2052 star-forming galaxies from the volume-limited, dwarf-dominated, z~0 ECO catalog. The MULTIGRIS statistical framework is used to evaluate the performance of various models using strong lines as constraints. The reference model involves physical parameters distributed as power-laws with free parameter boundaries. Specifically, we use combinations of 1D photoionization models (i.e., considering the propagation of radiation toward a single cloud) to match optical HII region lines, in order to provide probability density functions of the inferred parameters. The inference predicts non-uniform physical conditions within galaxies. The integrated spectra of most galaxies are dominated by relatively low-excitation gas with a metallicity around 0.3 solar. Using the average metallicity in galaxies, we provide a new fit to the mass-metallicity relationship which is in line with direct abundance method determinations from the calibrated range at low metallicity to stacks at high metallicity. The average metallicity shows a weakly bimodal distribution which may be due related to external (e.g., refueling of non-cluster early-type galaxies above ~10^9.5 solar masses) or internal processes (more efficient star-formation in metal-rich regions). The specific line set used for inference affects the results and we identify potential issues with the use of the [SII] line doublet. Complex modelling approaches are limited by the inherent 1D model database as well as caveats regarding the gas geometry. Our results highlight, however, the possibility to extract useful and significant information from integrated spectra.

We investigate the evolution of initial fractal clusters at 3 kpc from the Galactic Center (GC) of the Milky Way and show how red supergiant clusters (RSGCs)-like objects, which are considered to be the result of active star formation in the Scutum complex, can form by 16 Myr. We find that initial tidal filling and tidal over-filling fractals are shredded by the tidal force, but some substructures can survive as individual subclusters, especially when the initial virial ratio is $\leq$this http URL surviving subclusters are weakly mass segregated and show a top-heavy mass function. This implies the possibility that a single substructured star cluster can evolve into multiple `star clusters'.

Astrometric measurements provide a unique avenue for constraining the stochastic gravitational wave background (SGWB). In this work, we investigate the application of two neural network architectures, a fully connected network and a graph neural network, for analyzing astrometric data to detect the SGWB. Specifically, we generate mock Gaia astrometric measurements of the proper motions of sources and train two networks to predict the energy density of the SGWB, $\Omega_\text{GW}$. We evaluate the performance of both models under varying input datasets to assess their robustness across different configurations. Our results demonstrate that neural networks can effectively measure the SGWB, showing promise as tools for addressing systematic uncertainties and modeling limitations that pose challenges for traditional likelihood-based methods.

Machine learning (ML) methods have become popular for parameter inference in cosmology, although their reliance on specific training data can cause difficulties when applied across different data sets. By reproducing and testing networks previously used in the field, and applied to 21cmFast and Simfast21 simulations, we show that convolutional neural networks (CNNs) often learn to identify features of individual simulation boxes rather than the underlying physics, limiting their applicability to real observations. We examine the prediction of the neutral fraction and astrophysical parameters from 21 cm maps and find that networks typically fail to generalise to unseen simulations. We explore a number of case studies to highlight factors that improve or degrade network performance. These results emphasise the responsibility on users to ensure ML models are applied correctly in 21 cm cosmology.

Parker Solar Probe's (PSP) discovery of the prevalence of switchbacks (SBs), localised magnetic deflections in the nascent solar wind, has sparked interest in uncovering their origins. A prominent theory suggests these SBs originate in the lower corona through magnetic reconnection processes, closely linked to solar jet phenomena. Jets are impulsive events, observed across scales and solar atmosphere layers, associated with the release of magnetic twist and helicity. This study examines whether self-consistent jets can form and propagate into the super-Alfvénic wind, assesses the impact of distinct Parker solar wind profiles on jet dynamics, and determines if jet-induced magnetic untwisting waves display signatures typical of SBs. We employed parametric 3D numerical MHD simulations using the ARMS code to model the self-consistent generation of solar jets. Our study focuses on the propagation of solar jets in distinct atmospheric plasma $\beta$ and Alfvén velocity profiles, including a Parker solar wind. Our findings show that self-consistent coronal jets can form and propagate into the super-Alfvénic wind. Notable structures such as the leading Alfvénic wave and trailing dense-jet region were consistently observed across diverse plasma $\beta$ atmospheres. The jet propagation dynamics are significantly influenced by atmospheric variations, with changes in Alfvén velocity profiles affecting the group velocity and propagation ratio of the leading and trailing structures. U-loops, prevalent at jet onset, do not persist in the low-$\beta$ corona, but magnetic untwisting waves associated with jets show SB-like signatures. However, full-reversal SBs were not observed. These findings may explain the absence of full reversal SBs in the sub-Alfvénic wind and illustrate the propagation of magnetic deflections through jet-like events, shedding light on possible SB formation processes.

We present a detailed time-resolved and phase-resolved polarimetric analysis of the transient X-ray pulsar RX J0440.9+4431/LS V +44 17, using data from the Imaging X-ray Polarimetry Explorer (IXPE) during the 2023 giant outburst. We conducted a time-resolved analysis by dividing the data into several intervals for each observation. This analysis reveals a continuous rotation of the phase-averaged polarization angle (PA) across the observations performed during the super-critical and sub-critical regimes. To investigate the origin of the PA rotation, we performed a pulse phase-resolved polarimetric analysis over four time intervals, each spanning approximately three days. Applying the rotating vector model (RVM), the geometric parameters of the system were determined for each interval. Despite the short time gap of just $\sim$ 20 days, we observed significant variation in the RVM parameters between the first interval and the subsequent three, indicating the presence of an additional polarized component alongside the RVM component. Using a two-polarized component model with the assumption that this additional component remains constant across pulse phases, we calculated the phase-averaged PA and polarized flux of both the variable and constant components. The phase-averaged PA of each component remained relatively stable over time, but the polarized flux of the constant component decreased, while that of the variable component increased. The observed rotation of the PA is attributed to the gradual shift in the polarized flux ratio between the two components and is not directly related to the different accretion regimes.

Two similar-looking, two-part interplanetary type II burst events from 2003 and 2012 are reported and analysed. The 2012 event was observed from three different viewing angles, enabling comparisons between the spacecraft data. In these two events, a diffuse wide-band type II radio burst was followed by a type II burst that showed emission at the fundamental and harmonic (F-H) plasma frequencies, and these emission bands were also slightly curved in their frequency-time evolution. Both events were associated with high-speed, halo-type coronal mass ejections (CMEs). In both events, the diffuse type II burst was most probably created by a bow shock at the leading front of the CME. However, for the later-appearing F-H type II burst there are at least two possible explanations. In the 2003 event there is evidence of CME interaction with a streamer, with a possible shift from a bow shock to a CME flank shock. In the 2012 event a separate white-light shock front was observed at lower heights, and it could have acted as the driver of the F-H type II burst. There is also some speculation on the existence of two separate CMEs, launched from the same active region, close in time. The reason for the diffuse type II burst being visible only from one viewing direction (STEREO-A), and the ending of the diffuse emission before the F-H type II burst appears, still need explanations.

Core-collapse supernovae constitute a unique laboratory for particle physics and astrophysics. They are powerful neutrino sources of all flavors, emitting essentially all the gravitational binding energy through neutrinos, at the end of their life. I will highlight how crucial is the observation of the next core-collapse supernova and of the diffuse supernova neutrino background, whose discovery might be imminent.

The Pierre Auger Observatory has detected downward terrestrial gamma-ray flashes (TGFs) with its Surface Detector. A key to understanding this high-energy radiation in thunderstorms is to combine such measurements with measurements of lightning processes in their earliest stages. With eleven modified Auger Engineering Radio Array (AERA) stations we can build an interferometric lightning detection array working in the bandwidth between 30 - 80 MHz inside the Surface Detector array to precisely measure lightning stepped leaders in 3D. These measurements allow us to decipher the cause of TGFs and clarify the reason for the observed high-energy particles in thunderstorms. We will present the current status of the detection plans including the configuration of the interferometric lightning detection array and the steps to take as well as the reconstruction characteristics obtained with AERA.

We have detected three new hydrogen-deficient (H < 0.001 mass fraction) pre-white dwarfs (WDs) with helium-dominated atmospheres. The first object is a relatively cool PG1159 star (effective temperature Teff = 72,000 K) that has the lowest surface gravity of any PG1159 star known (log g = 4.8). It is a PG1159 star in the earliest pre-WD phase. The second object is a hot subdwarf O (sdO) star (Teff = 50,000 K, log g = 5.3) with high carbon and oxygen abundances. It is only the third known member of the recently established CO-sdO spectral class, which comprises stars that are thought to be formed by a merger of a disrupted low-mass CO WD with a higher-mass He WD. The third object is one of the rare stars of spectral type O(He) (Teff = 90,000 K, log g = 5.5).

We consider the possibility that the stochastic gravitational wave (GW) background suggested by Pulsar Timing Array (PTA) datasets is sourced by Primordial Black Holes (PBHs). Specifically, we perform a Bayesian search in the International PTA Data Release 2 (IPTA DR2) for a combined GW background arising from scalar perturbations and unresolved PBH mergers, assuming a broad PBH mass distribution. In our analysis, we incorporate constraints on the curvature power spectrum from CMB $\mu$-distortions and the overproduction of PBHs, which significantly suppress the contribution of PBH mergers to the total GW background. We find that scalar-induced GWs dominate the nHz frequency range, while PBH mergers alone cannot account for the observed signal under the standard PBH formation scenario involving Gaussian perturbations, and including only Poissonian PBH clustering. However, specific PBH models, such as those with enhanced clustering, could yield a GW background dominated by PBH mergers. Overall, we find that the IPTA DR2 strongly favors an astrophysical origin for the reported common-spectrum process over the PBH models considered in this analysis.

We report on IXPE and NuSTAR observations that began forty days following the onset of the 2024 outburst of the magnetar 1E 1841-045, marking the first ever IXPE observation of a magnetar in an enhanced state. Our spectropolarimetric analysis indicates that a non-thermal double power-law (PL) spectral model can fit the phase-averaged intensity data well, with the soft and hard components dominating below and above around 5 keV, respectively. We find that the soft PL exhibits a polarization degree (PD) of about 20% while the hard X-ray PL displays a PD of about 50%; both components have a polarization angle (PA) compatible with 0 degree. These results are supported through model-independent polarization analysis which shows an increasing PD from about 15% to 70% in the 2-3 keV and 6-8 keV ranges, respectively, while the PA remains consistent with 0 degree. We find marginal evidence for variability in the polarization properties with pulse phase, namely a higher PD at spin phases coinciding with the peak in the hard X-ray pulse. We compare the hard X-ray PL to the expectation from direct resonant inverse Compton scattering (RICS) and secondary pair cascade synchrotron radiation from primary high-energy RICS photons, finding that both can provide reasonable spectropolarimetric agreement with the data, yet, the latter more naturally. Finally, we suggest that the soft power law X-ray component may be emission emanating from a Comptonized corona in the inner magnetosphere.

Recent measurements of the galaxy 4-Point Correlation Function (4PCF) have seemingly detected non-zero parity-odd modes at high significance. Since gravity, the primary driver of galaxy formation and evolution is parity-even, any parity violation, if genuine, is likely to have been produced by some new parity-violating mechanism in the early Universe. Here we investigate an inflationary model with a Chern-Simons interaction between an axion and a $U(1)$ gauge field, where the axion itself is the inflaton field. Evaluating the trispectrum (Fourier-space analog of the 4PCF) of the primordial curvature perturbations is an involved calculation with very high-dimensional loop integrals. We demonstrate how to simplify these integrals and perform all angular integrations analytically by reducing the integrals to convolutions and exploiting the Convolution Theorem. This leaves us with low-dimensional radial integrals that are much more amenable to efficient numerical evaluation. This paper is the first in a series in which we will use these results to compute the full late-time 4PCF for axion inflation, thence enabling constraints from upcoming 3D spectroscopic surveys such as Dark Energy Spectroscopic Instrument (DESI), Euclid, or Roman.

The fact that the standard dispersion relation for photons in vacuum could be modified because of their interaction with the quantum nature of spacetime has been proposed more than two decades ago. A quantitative model [Jacob \& Piran, JCAP 01, 031 (2008)], has been tested extensively using distant highly energetic astrophysical sources, searching for energy-dependent time delays in photon arrival times. Since no delay was firmly measured, lower limits were set on the energy scale $\Lambda$ related to these effects. In recent years, however, different but equally well-grounded expressions beyond the Jacob \& Piran model were obtained for the photon dispersion relation, leading to different expressions for the dependence of lag versus redshift. This article introduces a general parameterization of modified dispersion relations in cosmological symmetry, which directly leads to a general parameterized lag versus redshift dependence encompassing both existing and new models. This parameterization could be used in the future to compare the predicted time lags of the different models and test them against observations. To investigate this possibility, realistic data sets are simulated, mimicking different types of extragalactic sources as detected by current and future instruments. When no lag is injected in the simulated data, each lag-redshift model leads, as expected, to a different value for the limit on $\Lambda$, and the Jacob \& Piran model gives the most stringent bound. When a lag at $\Lambda \sim E_P$ in the Jacob \& Piran model is injected, it is detected for all the other lag-redshift relations considered, although leading to different values. Finally, the possibility to discriminate between several lag-redshift models is investigated, emphasizing the importance of an evenly distributed sample of sources across a wide range of redshifts.

The very high-energy afterglow in GRB 221009A, known as the `Brightest Of All Time' (B.O.A.T.), has been thoroughly analyzed in previous studies. In this paper, we conducted a statistical analysis of the waiting time behavior of 172 TeV photons from the B.O.A.T. observed by LHAASO-KM2A. The following results were obtained: (I) The waiting time distribution (WTD) of these photons deviates from the exponential distribution. (II) The behavior of these photons exhibits characteristics resembling those of a self-organized critical system, such as power-law distribution and scale-invariance features in the waiting time distribution. The power-law distribution of waiting times is consistent with the prediction of a non-stationary process. (III) The relationship between the power-law slopes of the WTD and the scale-invariant characteristics of the Tsallis q-Gaussian distribution deviates from existing theory. We suggest that this deviation is due to the photons not being completely independent of each other. In summary, the power-law and scale-free characteristics observed in these photons imply a self-organized critical process in the generation of TeV photons from GRB 221009A. Based on other relevant research, we propose that the involvement of a partially magnetically dominated component and the continuous energy injection from the central engine can lead to deviations in the generation of TeV afterglow from the simple external shock-dominated process, thereby exhibiting the self-organized critical characteristics mentioned above.

From a spectroscopic survey of candidate H II regions in the Andromeda galaxy (M31) with MMT/Hectospec, we have identified 294 H II regions using emission line ratios and calculated elemental abundances from strong-line diagnostics (values ranging from sub-solar to super-solar) producing both Oxygen and Nitrogen radial abundance gradients. The Oxygen gradient is relatively flat, while the Nitrogen gradient is significantly steeper, indicating a higher N/O ratio in M31's inner regions, consistent with recent simulations of galaxy chemical evolution. No strong evidence was found of systematic galaxy-scale trends beyond the radial gradient. After subtracting the radial gradient from abundance values, we find an apparently stochastic and statistically significant scatter of standard deviation 0.06 dex, which exceeds measurement uncertainties. One explanation includes a possible collision with M32 200 - 800 Myrs ago. Using the two-point correlation function of the Oxygen abundance, we find that, similar to other spiral galaxies, M31 is well-mixed on sub-kpc scales but less so on larger (kpc) scales, which could be a result of an exponential decrease in mixing speed with spatial scale, and the aforementioned recent merger. Finally, the MMT spectroscopy is complemented by a dust continuum and CO survey of individual Giant Molecular Clouds, conducted with the Submillimeter Array. By combining the MMT and SMA observations, we obtain a unique direct test of the Oxygen abundance dependence of the $\alpha^{\prime}(^{12}\mathrm{CO})$ factor which is crucial to convert CO emission to dust mass. Our results suggest that within our sample there is no trend of the $\alpha^{\prime}(^{12}\mathrm{CO})$ with Oxygen abundance.

Recent observations from the James Webb Space Telescope revealed a surprisingly large number of galaxies formed at high redshift. Along with strong lensing studies and nearby galaxy observations, these could challenge the standard Lambda Cold Dark Matter cosmology with a power-law primordial power spectrum. In this study, we conduct high-resolution cosmological zoom-in dark matter-only simulations of Milky Way host size halos with a blue, tilted primordial power spectrum ($P(k)\propto k^{m_s}$ with $m_s>1$ at small scales $>1~{\rm Mpc}^{-1}$). We find that the blue-tilted subhalo mass functions can be enhanced by more than a factor of two for subhalo masses $M_{\rm sub} \lesssim 10^{10}~ M_{\odot}$, whereas the subhalo $V_{\rm max}$ functions can be enhanced by a factor of four for maximum circular velocities $V_{\rm max}\lesssim 30 ~{\rm km/s}$. The blue-tilted scaled cumulative substructure fraction can be an order of magnitude higher at $\sim$10\% of the virial radius. The blue-tilted subhalos also have higher central densities, since the blue-tilted subhalos reach the same $V_{\rm max}$ at a smaller distance $R_{\rm max}$ from the center. We have also verified these findings with higher-resolution simulations.

In the context of Einstein-Maxwell-scalar theory with a nonminimal coupling between the electromagnetic and scalar field, we study linear (non)radial perturbations and nonlinear radial dynamics of spherically symmetric black holes. In a certain region of the parameter space, this theory admits hairy black holes with a stable photon sphere. This has a counterpart in the effective potential of linear perturbations, featuring multiple maxima and minima. The corresponding quasinormal mode spectrum contains long-lived modes trapped in the potential cavity and the time-domain linear response displays echoes, as previously observed for horizonless compact objects. Interestingly, the black-hole dynamics in this theory can be studied at the nonlinear level. By performing fully-fledged 1+1 simulations, we show that echoes are present even when the nonlinearities are significant. To our knowledge, this is the first example of echoes appearing in a consistent theory beyond a linearized analysis. In a follow-up work we will study whether this feature is also present in the post-merger signal from black hole collisions in this theory.

The cosmology of metric-affine gravity is studied for the general, parity preserving action quadratic in curvature, torsion and non-metricity. The model contains 27 a priori independent couplings in addition to the Einstein constant. Linear and higher order relations between the quadratic operators in a Friedmann--Lemaitre--Robertson--Walker spacetime are obtained, along with the modified Friedmann, torsion and non-metricity equations. Extra parameter constraints lead to two special branches of the model. Firstly, a branch is found in which the Riemannian spatial curvature (thought to be slightly closed or flat in the Lambda-CDM model of our Universe) is entirely screened from all the field equations, regardless of its true value. Secondly, an integrable branch is found which yields (anti) de Sitter expansion at late times. The particle spectra of these two branches are studied, and the need to eliminate higher-spin particles as well as ghosts and tachyons motivates further parameter constraints in each case. The most general model is also found which reproduces the exact Friedmann equations of general relativity. The full set of equations describing closed, open or flat cosmologies, for general parity-even quadratic metric-affine gravity, is made available for SymPy, Mathematica and Maple platforms.

Black-hole spacetimes that possess stationary equatorial matter rings are known to exist in general relativity. We here reveal the existence of black-hole spacetimes that support {\it non}-equatorial matter rings. In particular, it is proved that rapidly-rotating Kerr black holes in the dimensionless large-spin regime ${\bar a}>{\bar a}_{\text{crit}}= \sqrt{\big\{{{7+\sqrt{7}\cos\big[{1\over3}\arctan\big(3\sqrt{3}\big)\big]- \sqrt{21}\sin\big[{1\over3}\arctan\big(3\sqrt{3}\big)\big]}\big\}/12}}\simeq0.78$ can support a pair of non-equatorial massive scalar rings which are negatively coupled to the Gauss-Bonnet curvature invariant of the spinning spacetime (here ${\bar a}\equiv J/M^2$ is the dimensionless angular momentum of the central supporting black hole). We explicitly prove that these non-equatorial scalar rings are characterized by the dimensionless functional relation $-{{57+28\sqrt{21}\cos\big[{1\over3}\arctan\big({{1}\over{3\sqrt{3}}}\big)\big]} \over{8(1+\sqrt{1-{\bar a}^2})^6}} \cdot{{\bar\eta}\over{{\bar\mu}^2}}\to 1^{+}$ in the large-mass ${\bar\mu}\equiv M\mu\gg1$ regime (here $\{{\bar\eta}<0,\mu\}$ are respectively the non-trivial coupling parameter of the composed Einstein-Gauss-Bonnet-massive-scalar field theory and the proper mass of the supported non-minimally coupled scalar field).

The Aether Scalar Tensor (AeST) theory is an extension of general relativity(GR) successful at reproducing galactic rotational curves, gravitational lensing, linear large scale structure and cosmic microwave background power spectrum observations. We solve the most general static spherically symmetric vacuum equations in the strong-field regime of AeST and find two classes of stealth black hole solutions -- those with exact GR geometries -- containing non-trivial secondary hair. In particular, one of these can be continuously joined to the cosmological solution of AeST. We also derive a non-black hole solution with zero spatial component in the vector field. This result proves the existence of mathematically and observationally consistent candidates for black holes in AeST, and creates a basis for testing the theory in the strong-field regime.

Since the derivation of a well-defined $D\to4$ limit for 4D Einstein-Gauss-Bonnet (4DEGB) gravity coupled to a scalar field, there has been considerable interest in testing it as an alternative to Einstein's general theory of relativity. Past work has shown that this theory hosts interesting compact star solutions which are smaller in radius than a Schwarzschild black hole of the same mass in general relativity (GR), though the stability of such objects has been subject to question. In this paper we solve the equations for radial perturbations of neutron stars in the 4DEGB theory with SLy/BSk class EOSs, along with the MS2 EOS, and show that the coincidence of stability and maximum mass points in GR is still present in this modified theory, with the interesting additional feature of solutions re-approaching stability near the black hole solution on the mass-radius diagram. Besides this, as expected from past work, we find that larger values of the 4DEGB coupling $\alpha$ tend to increase the mass of neutron stars of the same radius (due to a larger $\alpha$ weakening gravity) and move the maximum mass points of the solution branches closer to the black hole horizon.

Composite asymmetric dark matter (ADM) is the framework that naturally explains the coincidence of the baryon density and the dark matter density of the Universe. Through a portal interaction sharing particle-antiparticle asymmetries in the Standard Model and dark sectors, dark matter particles, which are dark-sector counterparts of baryons, can decay into antineutrinos and dark-sector counterparts of mesons (dark mesons) or dark photon. Subsequent cascade decay of the dark mesons and the dark photon can also provide electromagnetic fluxes at late times of the Universe. We derive constraints on the lifetime of dark matter decay in the composite ADM scenario from the astrophysical observations of the $e^+$, $e^-$, and $\gamma$-ray fluxes. The constraints from cosmic-ray positron measurements by AMS-02 are the most stringent at $\gtrsim2$ GeV: a lifetime should be larger than the order of $10^{26}$ s, corresponding to the cutoff scale of the portal interaction of about $10^8 \text{--} 10^9 \, \mathrm{GeV}$. We also show the importance of neutrino observations with Super-Kamiokande and Hyper-Kamiokande, which give conservative bounds.

We demonstrate a photon proliferation effect from $N$-body dark matter (DM) annihilation in the early Universe, which can induce a drastic photon-temperature shift after neutrino decoupling. For pseudoscalar DM mass below the eV scale, we show that the photon proliferation effect becomes significant as the mass approaches the ultralight end, due to the huge enhancement from the background DM number density. This presents the leading constraints on the DM-photon coupling, DM self-interaction, and DM-electron coupling, which are stronger than the existing bounds up to several orders of magnitude. The present research can be extended to other interactions and DM candidates, and highlights the importance of multi-body processes in the early Universe.

Parity-violating electron scattering experiments on $\rm ^{48}Ca$ (CREX) and $\rm ^{208}Pb$ (PREX-II) offer valuable insight into the isovector properties of finite nuclei, providing constraints for the density dependence of the nuclear equation of state, which is crucial for understanding astrophysical phenomena. In this work, we establish functional dependencies between the properties of finite nuclei - such as weak charge form factors and neutron skin thickness - and the bulk properties of neutron stars, including tidal deformability from binary neutron star mergers and neutron star radii. The dependencies are formulated by introducing a family of $\beta$-equilibrated equations of state based on relativistic energy density functionals with point coupling interactions. The charge minus the weak form factors derived from CREX and PREX-II measurements, combined with the observational constraints on tidal deformability from the GW170817 event, are used to constrain the symmetry energy and neutron star radii. Notably, the energy density expanded up to the fourth order in symmetry energy yields larger radii compared to calculations limited to the second order term. However, the results reveal a discrepancy between the constraints provided by the CREX and PREX-II experiments. For a more quantitative assessment, higher precision parity-violating electron scattering data and neutron star observations are required.

We analyze the orbital dynamics of spherical test bodies in ``black hole surrounded by dark matter halo'' spherically symmetric spacetimes. When the test body pulsates periodically (such as a variable star), altering its quadrupole tensor, Melnikov's method shows that its orbital dynamics presents homoclinic chaos near the corresponding unstable circular orbits however small the oscillation amplitude is. Since for supermassive black holes the period of revolution of a star near the innermost stable circular orbit roughly spans time intervals from minutes to hours, the formalism can be applied in principle to the astrophysical scenario of a pulsating (variable) star inspiraling into a supermassive black hole, including the black hole SgrA* at the center of our Galaxy. The chaotic nature of its orbit, due to pulsation, is imprinted in the redshift time series of the emitted light and can, in principle, be observed in the corresponding light curves and even in gravitational-wave signals detected by future observatories such as the Laser Inteferometer Space Antenna. Also, although periodic with respect to the star's proper time, the chaotic orbital motion will produce an erratic light curve (and gravitational-wave signal) in terms of observed, coordinate time. Although our results were obtained for a specific exact solution, we argue that this phenomenon is generic for pulsating bodies immersed in black hole spacetimes surrounded by self-gravitating fluids.

The neutrino mass generation via conventional seesaw mechanism is realized at high scales around $O(10^{14})$GeV and probing new physics of the seesaw scale poses a great challenge. A striking fact is that the neutrino seesaw scale is typically around the cosmological inflation scale. In this work, we propose a framework incorporating inflation and neutrino seesaw in which the inflaton primarily decays into right-handed neutrinos after inflation. This decay process is governed by the inflaton interaction with the right-handed neutrinos that respects the shift symmetry. Under the neutrino seesaw mechanism, fluctuations of the Higgs field can modulate the inflaton decays, contributing to the curvature perturbation. We investigate the induced non-Gaussian signatures and demonstrate that such signatures provides an important means to probe the high-scale neutrino seesaw mechanism.