Papers for Thursday, Oct 23 2025

Astrometry, one of the oldest branches of astronomy, has been revolutionized by missions like Hipparcos and especially Gaia, which mapped billions of stars with extraordinary precision. However, challenges such as detecting Earth-like exoplanets in nearby habitable zones and probing the influence of dark matter in galactic environments require sub-microarcsecond accuracy. With a 6--8 meter large-aperture telescope operating across at visible wavelengths, the Habitable Worlds Observatory by NASA can combine astrometry and direct imaging to detect rocky exoplanets within 10 parsecs and study their atmospheres. We consider here the scientific requirements and present a concept for a dedicated astrometric instrument on HWO. It is capable to produce diffraction-limited images of large fields, achieving a point-spread function (PSF) precision of 20 milliarcseconds. Equipped with a detector calibration system, HWO can perform high precision astrometry, and, detect and measure the orbit of Earth-mass planets in the habitable zone of Nearby Solar-type stars. HWO can dramatically improve current constraints on the self- interaction cross-section of heavy dark matter particles (WIMPs) and on the masses of ultra-high dark matter particles, through the study of stellar motions in galactic environments. The visible channel of the instrument features a large CMOS-based focal plane with stitched pixel arrays, enabling a large field of view. The ``Detector Calibration Unit'' system uses interferometric laser fringes to calibrate pixel positions. Using differential astrometry and pointed observations with a stable telescope design enables extended integration times, enhancing sensitivity to sub-microarcsecond precision for detecting exoplanets and studying dark matter through stellar motion.

We present a comprehensive investigation into the influence of stellar bars on star formation (SF) in galaxy pairs, using a large sample of low-redshift galaxies ($0.02$\,$<$\,z\,$<$\,$0.08$) from the DESI Legacy Imaging Surveys DR8. Our analysis examines whether bars enhance or suppress SF during pair interactions, and how these outcomes depend on the star-forming properties of companion galaxies. We find that bars either catalyze or inhibit SF in their host galaxies, depending on the companion's SF activity. In particular, barred galaxies paired with actively star-forming companions experience more pronounced central starbursts (with sSFR up to $\sim$\,2.5 dex higher) than unbarred counterparts, whereas those with passive companions often have suppressed SF (sometimes below isolated galaxy levels). The notion of the dual role of bars can reconcile conventional conflicting reports of bar-driven enhancement versus quenching of SF activity. Bars, well known to regulate kpc-scale dynamics, may also link to the impact of external environments: when a star-forming companion provides sufficient gas, bars drive central starbursts, whereas in gas-poor interactions, bars hasten gas depletion and contribute to SF suppression. This work highlights the necessity of accounting for both internal structure and companion properties to fully understand SF regulation in interacting galaxies.

Galactic warps are common in disk galaxies. While often attributed to galaxy--galaxy tides, a non-spherical dark matter (DM) halo has also been proposed as a driver of disk warping. We investigate links among warp morphology, satellite distribution, and large-scale structure using the Sloan Digital Sky Survey catalog of warped disks compiled by Zee et al.\ (2022). Warps are classified into 244 S and 127 U types, hosting 1,373 and 740 satellites, respectively, and are compared to an unwarped control matched in stellar mass, redshift, and local density. As an indirect, population-level proxy for the host halo's shape and orientation, we analyze the stacked spatial distribution of satellites. Warped hosts show a significant anisotropy: an excess at $45^{\circ}<\phi<90^{\circ}$ (measured from the host major axis), peaking at $P(\phi)\simeq 0.003$, versus nearly isotropic controls. Satellites of S-type warps preferentially align with the nearest cosmic filament, whereas those of U-type warps are more often perpendicular. The incidence of warps increases toward filaments ($r_{\rm fil}<4,{\rm Mpc},h^{-1}$), while the number of satellites around warped hosts remains approximately constant with filament distance, indicating a direct influence of the large-scale environment. We discuss possible links between galactic warps and the cosmic web, including anisotropic tidal fields and differences in evolutionary stage.

We use the THESAN radiation-hydrodynamics simulations to investigate how Lyman-$\alpha$ emitters (LAEs) trace ionized bubble sizes during the Epoch of Reionization. We generate realistic LAE catalogs by combining accurate intrinsic Ly$\alpha$ production and intergalactic transmission with an empirical model for dust absorption and gas outflows. By calibrating to observationally-constrained Ly$\alpha$ luminosity functions, we reproduce the rapid decline in Ly$\alpha$ visibility toward higher redshifts while revealing mild tensions in LAE fractions near the end of reionization. Before the midpoint of reionization, galaxies within larger line-of-sight bubbles ($\gtrsim 10$ cMpc) have higher observed Ly$\alpha$ luminosity and equivalent width (EW), demonstrating that the evolving LAE fraction provides a practical statistical tracer for bubble size. These correlations weaken as percolation progresses and the IGM becomes increasingly ionized. In LAE selected samples with $L_{\text{Ly}\alpha} > 10^{41.5}\ \text{erg s}^{-1}$, Ly$\alpha$ properties correlate with bubble size more strongly than UV magnitude, especially at $z \gtrsim 7$. This simulation-based framework maps LAE selections to bubble-size statistics, clarifies biases in more idealized models, and will supply public catalogs to interpret current and forthcoming JWST and narrow-band LAE surveys in terms of the evolving topology of reionization.

The distribution of molecular gas on small scales regulates star formation and the growth of supermassive black holes in galaxy centers, yet the role of active galactic nuclei (AGN) feedback in shaping this distribution remains poorly constrained. We investigate how AGN influence the small-scale structure of molecular gas in galaxy centers, by measuring the clumpiness of CO(3 - 2) emission observed with the Atacama Large Millimeter/submillimeter Array (ALMA) in the nuclear regions (50 - 200 pc from the AGN) of 16 nearby Seyfert galaxies from the Galaxy Activity, Torus, and Outflow Survey (GATOS). To quantify clumpiness, we apply three different methods: (1) the median of the pixel-by-pixel contrast between the original and smoothed maps; (2) the ratio of the total excess flux to the total flux, after substracting the background smoothed emission; and (3) the fraction of total flux coming from clumpy regions, interpreted as the mass fraction in clumps. We find a negative correlation between molecular gas clumpiness and AGN X-ray luminosity (L_X), suggesting that higher AGN activity is associated with smoother gas distributions. All methods reveal a turnover in this relation around L_X = 10^{42} erg/s, possibly indicating a threshold above which AGN feedback becomes efficient at dispersing dense molecular structures and suppressing future star formation. Our findings provide new observational evidence that AGN feedback can smooth out dense gas structures in galaxy centers.

We present a near-infrared census of stellar large-amplitude variables (LAVs) observed by the Palomar Gattini-IR (PGIR) surveyor from 2019-2021. Over the three-year time period, PGIR performed a brightness-limited survey of the Northern sky (18,000 sq. deg) to J-band AB magnitudes of $\sim 13$ within and $\sim 15$ outside the Galactic plane. From 70 million stars detected in PGIR reference images, we provide a spectral and photometric library of the 128 largest amplitude stellar variables detected to median SNR > 10 for more than 50 epochs with more than 5 high-amplitude detections, peak-to-peak magnitudes $\geq$ 2, and von Neumann ratios $\leq$ 0.2. We obtained medium-resolution near-infrared spectra with TripleSpec on the 200-inch Hale Telescope at Palomar Observatory and SpeX at NASA's Infrared Telescope Facility. The spectral census consists of 82 evolved and dust-obscured Asymptotic Giant Branch stars, 16 R Coronae Borealis stars, 13 young-stellar or pre-main-sequence objects, 8 symbiotic binaries, 7 erratic carbon- and oxygen-rich giants, and 2 RV Tauri supergiants. The spectral and photometric dataset serves as an atlas of near-infrared LAVs and a repository of evolved stars, eruptive variables, and binary systems for future deeper infrared surveys.

Machine learning has the potential to improve the reconstruction of the dark matter profile of galaxies with respect to traditional methods, like rotation curves. We demonstrate on the simulation suite Illustris-TNG that a steerable equivariant convolutional neural network (CNN) is able to infer the dark matter profiles within and around individual galaxies from photometric and interferometric data, improving on a standard CNN. Within the in silico environment of the simulations, our architecture is able to capture the dark matter distribution within galaxies without a parametrization of the profile. We perform an interpretability analysis to understand the internal mechanisms of the trained model and the most important data features used to estimate the dark matter profiles. The equivariant CNN recovers the dark matter profile of galaxies within the stellar mass range $[10^{10} - 10^{12} ]$ $M_{\odot}$ with excellent precision and accuracy: the mean squared error is reduced by a factor of ~ 3 from its value under the training distribution, demonstrating that the network has learnt from the data features. While this holds within the controlled 'in silico' environment of the simulation, we argue that few additional steps are needed before this method can be reliably applied to galaxies in the real field observations.

Stripped stars are hot, helium-rich stars formed when binary interactions remove a star's hydrogen envelope. While low-mass ($\lesssim 1\,M_\odot$) and high-mass ($\gtrsim 8\,M_\odot$) stripped stars are well studied as hot subdwarfs and Wolf-Rayet stars, their intermediate-mass counterparts ($1-8 \,M_\odot$) have only recently been discovered. The Stripped-Star Ultraviolet Magellanic Cloud Survey (SUMS) identified UV-excess sources in the Magellanic Clouds using Swift-UVOT photometry and selected 820 photometric stripped-star candidates. However, the selection function, completeness, and purity of this sample remain poorly understood. We forward model the population of stripped stars in the Magellanic Clouds using a binary population synthesis model combined with spatially resolved star formation histories and simulated UV photometry. To assess survey sensitivity, we inject simulated sources into real Swift-UVOT images and reproduce the SUMS selection process, including crowding, extinction, and photometric quality cuts. We recover $\sim 250$ simulated stripped stars with masses $\gtrsim 1\,M_{\odot}$, which corresponds to a recover rate of $9-13\,\%$. The rest are missed due to dilution by luminous companions, crowding, and high extinction. The observed population is biased toward systems with low-mass companions formed by common envelope evolution and toward systems with compact object companions. Of the stripped stars which show UV excess, $25-45\,\%$ are identifiable by SUMS; higher-resolution data or improved reddening corrections are needed to detect the rest. We predict contamination of the observed stripped star candidates by main-sequence stars with spurious UV excess due to crowding and provide guidelines for selecting higher-purity subsamples. These results enable tests of binary evolution models and realistic comparison of observed and simulated stripped-star populations.

Relativistic macroscopic plasma dynamics can be described by general-relativistic magnetohydrodynamics. In many high-energy astrophysical settings, such as the interior dynamics of magnetized stars, the ideal GRMHD approximation, in which we assume infinite conductivity, provides an excellent description. However, ideal GRMHD neglects resistive effects that are essential for processes such as magnetic reconnection, dissipation, and magnetospheric dynamics. Incorporating resistivity into astrophysical plasma models accounts for the fact that plasmas in such environments are not perfect conductors. We present a resistive version of the GPU-accelerated GRMHD code GRaM-X, which evolves the full resistive GRMHD equations using the Z4c formalism for Einstein's equations. We implement a second-order implicit-explicit Runge-Kutta scheme to handle stiff source terms, obtain the primitive quantities from the conserved quantities using a one-dimensional recovery method, and employ the HLLE Riemann solver in combination with TVD and WENO reconstruction schemes. We validate the module using a range of standard tests, including 1D shocktubes, current sheets, Alfvén waves, 2D cylindrical explosions, and 3D TOV stars. The results of these tests demonstrate accurate recovery of the ideal MHD limit, correct resistive behavior, and stable evolution in dynamical spacetimes. Leveraging the GPU-accelerated resistive version of GRaM-X enables efficient large-scale simulations, paving the way for realistic studies of binary mergers, accretion flows, and relativistic jets within the framework of multi-messenger astrophysics.

Recent studies on planet-dominated Type II migration demonstrated the presence of a correlation between the direction of planet migration and the parameter K describing the depth of the planetary gap. It was found that high (low) value for K correspond to outward (inward) migration. In this paper we aim at understanding the mechanism driving inward/outward migration and why it correlates with the gap depth. We performed a suite of 2D, live-planet, long-term simulations of massive planets migrating in discs with the hydro-code Fargo3D. We focus on a range of planet masses (1-13 m_J) and disc aspect ratios (0.03-0.1) and analyze the evolution of orbital elements and gap structure. We also study the torque contributions from outer Lindblad resonances to investigate their role in the migration outcome. We find that, while all planets initially migrate inwards, those with high enough K eventually enter a phase in which the torque reverses sign and migration becomes outwards, until eventually stalling. This behavior is associated with eccentricity growth in the outer disc and changes in the gap structure. We identify the surface density ratio at the 1:2 and 1:3 outer Lindblad resonances as a key output diagnostic that correlates with the migration direction. This ratio regulates the migration for all the cases where the massive planet remains in an almost circular orbit and the outer gap region exhibits moderate eccentricity. This characteristic sequence of inward-reversal-outwards-stalling occurs for a variety of K values and thus further work is required to identify the simulation input parameters that determine the onset of this sequence. Our results suggest that outward migration in the planet-dominated regime is primarily governed by the relative importance of the 1:2 and 1:3 resonances and, therefore, the gap profile plays a crucial role in determining the direction of migration.

We use the TNG50 cosmological hydrodynamic simulation to study the accreted stellar component and stellar haloes of isolated galaxies spanning a wide range of masses ($10^8<M_*/M_\odot<10^{11}$). We find that stars formed in the main progenitor (i.e., in-situ stars) typically dominate the inner regions as far as $\sim$10 half-light radii from the centre, implying that detecting uncontrovertible evidence for the presence of an accreted stellar halo requires probing the far outskirts of a galaxy. Stars from accreted, disrupted satellites (i.e., ex-situ stars) dominate beyond that radius (roughly $25\%$ of the virial radius, $r_{200}$), which we identify as the inner boundary of the outer stellar halo. The fraction of accreted stars decreases monotonically with decreasing galaxy mass, $M_*$, from $\sim$$20\%$ on average in $\sim$$2\times 10^{12}\, M_\odot$ haloes ($M_*\sim$$10^{11}\, M_\odot$) to $2$-$3\%$ in $\sim$$2\times 10^{10}\, M_\odot$ haloes ($M_*\sim$$10^{8}\, M_\odot$). The outer halo has a mass comparable to roughly $10\%$ of all accreted stars. Fewer than $\sim$$30\%$ of stars in the outer halo are in-situ stars, many of which originate from star-forming satellites during the late stages of disruption, especially in low-mass systems. Accreted stars are systematically more metal poor in less massive systems, which makes the outer haloes of dwarf galaxies a fertile hunting ground for extremely metal-poor stars. At given galaxy mass, the more massive stellar haloes are systematically more concentrated (smaller $R_{\rm eff}$) and have steeper density profiles (larger $n$). Our results provide a blueprint for interpreting observations of the outskirts of isolated galaxies in terms of their assembly histories.

We present a new formalism to separate large- and small-scale contributions to cosmic shear through $\textit{lensing counterterms}$ (LCT) inspired by effective field theory (EFT). Marginalizing over these LCTs isolates the large-scale cosmological signal in weak lensing power spectra while simultaneously constraining the impact of baryonic feedback or new physics (e.g. axion dark matter) at small scales. Our formalism removes the need for hard scale cuts in standard analyses, even when theoretical predictions are limited to below a physical cutoff $\Lambda$, resulting in significant improvements in constraining power -- up to $5\times$ smaller in the case of a LSST-Y10-like analysis without marginalizing over baryons when the analysis cutoff is set to $\Lambda = 1.0h$ Mpc$^{-1}$. We conduct a proof-of-principle analysis on the publicly available DES Y3 data, finding $S_8= 0.761\pm 0.045$ and $S_8 = 0.794\pm 0.040$ for analyses with cutoffs of $\Lambda = 0.5h$ Mpc$^{-1}$ and $1.0 h$ Mpc$^{-1}$, respectively, with no detection of modifications to small-scale clustering at $k > \Lambda$ beyond the predictions of collisionless dark matter in a $\Lambda$CDM universe. We make our $\texttt{JAX}$-based pipeline, $\texttt{gholax}$, integrated with intrinsic alignment predictions from the EFT of large-scale structure at 1-loop, publicly available.

Using HST/STIS observations, we present the highest-spatial-resolution spectroscopic study to date of four tidal disruption event (TDE) host galaxies, with the best observed being the post-starburst (PSB) host of ASASSN-14li. The stellar population of ASASSN-14li's host, within 44 pc of the nucleus, reveals a younger recent starburst ($\sim$340 Myr) compared to the population at an offset radius of 88 pc that excludes the nucleus ($\sim$550 Myr), a radial age gradient suggesting gas inflows from a minor merger. We estimate a stellar density of $\sim5900 \pm 800 \, M_\odot / \mathrm{pc}^3$ within 30 pc of the nucleus of ASASSN-14li's host, exceeding densities expected for nuclear star clusters. High-ionization ``coronal" emission lines, [Fe VI] $\lambda 5677$, [Fe VII] $\lambda 6087$, and [Fe X] $\lambda 6375$, are also detected within the nuclear spectra of the hosts of ASASSN-14li and PTF09ge, importantly alongside the non-detection of [O III] at the same scale. We similarly do not detect [O III] in the nuclear region of ASASSN-14ae's host despite its presence in the SDSS spectrum. The different ionization radiation levels detected at various radii from TDE host nuclei may indicate echoes of earlier accretion episodes, including, potentially, a prior TDE. We posit that a minor merger driving gas inflow to the nucleus could drive the enhanced TDE rates in post-starburst galaxies, inducing variation in nuclear gas properties and star formation history on $<$150 pc scales in TDE hosts.

Observations have shown that the optical colors of Galactic cirrus clouds differ significantly from those of extragalactic sources; thus, they can be used to distinguish Galactic cirrus from extragalactic low surface brightness (LSB) features. To understand these properties, we calculate radiative transfer models in dust clouds, where photons are incident from the ambient interstellar medium (ISM). Dust clouds are modeled to mimic a turbulent medium using a fractional Brownian motion algorithm, resulting in a lognormal density distribution and a power-law power spectral density that are appropriate for the ISM. The results are compared with optical observations of cirrus clouds in the Stripe 82 region. The observed color--color ($g-r$, $r-i$, and $i-z$) diagrams of the cirrus clouds can be reproduced by scattered light if the interstellar radiation field (ISRF) of Mathis et al. (as updated by Draine) is modified, either by reducing the intensities in the $i$ and $z$ bands or by enhancing those in the $g$ and $r$ bands. Similar results can also be obtained by adjusting the scattering albedos at the corresponding wavelengths. This demonstrates that the color--color diagrams are effective not only for identifying extragalactic LSB features but also for studying the ISRF and the properties of interstellar dust.

A photometric study in combination with existing stellar models has revealed details of this eclipsing post-mass-transfer binary. The shell star has an equatorial/ polar radius of ~2.60/1.90 Rsun at an equatorial rotational velocity of ~430 km s-1, an effective mean temperature Teff of ~13300 K and a mass of ~3.42 Msun. This former accretor star is surrounded by a large decretion disk of ~47 Rsun. The secondary star is a helium white dwarf precursor with a radius of 0.98 Rsun, a Teff of ~17100 K and a mass of 0.29 Msun. The parameters of this former donor star indicate an age of the binary system of only ~1.8 Myr after the end of mass transfer. The results fit to a sub-solar metallicity of Z = 0.007.

Galaxies exhibit a tight correlation between their star-formation rate and stellar mass over a wide redshift range known as the star-forming main sequence (SFMS). With JWST, we can now investigate the SFMS at high redshifts down to masses of $\sim10^6$ M$_{\odot}$, using sensitive star-formation rate tracers such as H$\alpha$ emission -- which allow us to probe the variability in star formation histories. We present inferences of the SFMS based on 316 H$\alpha$-selected galaxies at $z\sim4$-$5$ with $\log(\rm M_\star/M_\odot) = 6.4$ -$10.6$. These galaxies were identified behind the Abell 2744 lensing cluster with NIRCam grism spectroscopy from the ``All the Little Things'' (ALT) survey. At face value, our data suggest a shallow slope of the SFMS (SFR $\propto \mathrm{M}_\star^\alpha$, with $\alpha=0.45$). After correcting for the H$\alpha$-flux limited nature of our survey using a Bayesian framework, the slope steepens to $\alpha = 0.59^{+0.10}_{-0.09}$, whereas current data on their own are inconclusive on the mass dependence of the scatter. These slopes differ significantly from the slope of $\approx1$ expected from the observed evolution of the galaxy stellar mass function and from simulations. When fixing the slope to $\alpha=1$, we find evidence for a decreasing intrinsic scatter with stellar mass (from $\approx 0.5$ dex at M$_\star=10^8$ M$_\odot$ to $0.4$ dex at M$_\star=10^{10}$ M$_\odot$). This tension might be explained by a (combination of) luminosity-dependent SFR(H$\alpha$) calibration, a population of (mini)-quenched low-mass galaxies, or underestimated dust attenuation in high-mass galaxies. Future deep observations across facilities can quantify these processes, enabling better insights into the variability of star formation histories.

We explore the Hubble tension within an anisotropic cosmological framework by revisiting the Bianchi type-I model introduced in Le Delliou et al. 2020. Motivated by ongoing debates surrounding back-reaction effects and observed anomalies in the cosmic microwave background (CMB), we investigate whether a departure from isotropy in the late Universe could reconcile the observed discrepancies in Hubble constant measurements. Using a Bayesian inference framework, we constrain the model parameters employing multiple nested sampling algorithms: bilby, PyMultiNest, and nessai. We perform the analysis under both uniform and Gaussian priors, allowing us to systematically assess the sensitivity of the inferred cosmological parameters to different prior assumptions. This dual-prior strategy balances agnostic parameter exploration with constraints informed by theory and observation. Our findings demonstrate the reliability of our inference pipeline across different samplers and emphasize the crucial role of prior selection in non-standard cosmological model testing. The results suggest that anisotropic models remain viable contenders in addressing current cosmological tensions: even though the present model does not show alleviation of the Hubble tension, the data points towards anisotropies. Future work may extend this methodology to more complex anisotropic scenarios and incorporate additional cosmological probes such as CMB polarization and gravitational wave standard sirens.

Using very long baseline interferometry (VLBI) observations at (sub)millimeter wavelengths, the Event Horizon Telescope (EHT) currently achieves the finest angular resolution of any astronomical facility, necessary for imaging the horizon-scale structure around supermassive black holes. A significant calibration challenge for high-frequency VLBI stems from rapid variations in the atmospheric water vapor content above each telescope in the array, which induce corresponding fluctuations in the phase of the correlated signal that limit the coherent integration time and thus the achievable sensitivity. In this paper, we introduce a model that describes station-based phase corruptions jointly with a parameterization for the source structure. We adopt a Gaussian Process (GP) prescription for the time evolution of these phase corruptions, which provides sufficient flexibility to capture even highly erratic phase behavior. The use of GPs permits the application of a Kalman filtering algorithm for numerical marginalization of these phase corruptions, which permits efficient exploration of the remaining parameter space. Our model also removes the need to specify an arbitrary ``reference station'' during calibration, instead establishing a global phase zeropoint by enforcing the GPs at all stations to have fixed mean and finite variance. We validate our method using a real EHT observation of the blazar 3C 279, demonstrating that our approach yields calibration solutions that are consistent with those determined by the EHT Collaboration. The model presented here can be straightforwardly extended to incorporate frequency-dependent phase behavior, such as is relevant for the frequency phase transfer calibration technique.

Surrounding Sgr A*, a cluster of young and massive stars coexist with a population of dust-enshrouded objects, whose astrometric data can be used to scrutinize the nature of Sgr A*. An alternative to the black hole (BH) scenario has been recently proposed in terms of a supermassive compact object composed of self-gravitating fermionic dark matter (DM). Such horizon-less configurations can reproduce the relativistic effects measured for S2 orbit, while being part of a single continuous configuration whose extended halo reproduces the latest GAIA-DR3 rotation curve. In this work, we statistically compare different fermionic DM configurations aimed to fit the astrometric data of S2, and five G-sources, and compare with the BH potential when appropriate. We sample the parameter spaces via Markov Chain Monte Carlo statistics and perform a quantitative comparison estimating Bayes factors for models that share the same likelihood function. We extend previous results of the S2 and G2 orbital fits for 56 keV fermions (low core-compactness) and show the results for 300 keV fermions (high core-compactness). For the selected S2 dataset, the former model is slightly favoured over the latter. However, more precise S2 datasets, as obtained by the GRAVITY instrument, remain to be analysed in light of the fermionic models. For the G-objects, no conclusive preference emerges between models. For all stellar objects tested, the BH and fermionic models predict orbital parameters that differ by less than 1%. More accurate data, particularly from stars closer to Sgr A*, is necessary to statistically distinguish between the models considered.

Nova Scorpii 2023 was first detected as a luminous supersoft X-ray source (SSS) 93 days after outburst and continued emitting soft X-rays for over two months, until it was too close to the Sun to observe. The nova was monitored with the Swift X-ray Telescope (XRT) and the Neutron Star Interior Composition Explorer (NICER) on the International Space Station, and in long exposures with the Chandra High Resolution Camera (HRC) and Low Energy Transmission Grating (LETG) on days 128, 129, and 183-185 after optical maximum. Swift detected a rapidly decaying SSS when observations resumed, constraining the constant bolometric luminosity phase to 9 months. The SSS flux was irregularly variable. A nearly three-fold increase in flux was observed between August and October 2023 in the 15 to 35 Angstrom range, from 3.5 x 10^(-11) to 9.4 x 10^(-11) erg cm^(-2) s^(-1). The SSS duration and effective temperature derived from the October LETG spectra indicate a massive white dwarf with temperature fitting nova evolutionary tracks for a 1.2 solar mass WD; emission lines superimposed on the WD continuum are attributed to surrounding shocked ejecta. We present a timing study based on Chandra and archival NICER data. The irregular variability timescale was days, but a 77.9 second periodic modulation in the SSS flux with varying amplitude was measured in many observations. Our analysis shows that this period was stable; short drifts derived with NICER, but not in long, uninterrupted Chandra exposures, are artifacts of measuring variable amplitude modulation. We suggest the modulations are associated with the WD rotation.

Understanding how active-region properties influence coronal mass ejection (CME) dynamics is essential for constraining eruption models and improving space-weather prediction. Magnetic diagnostics derived above polarity inversion lines (PILs), including the critical height ($h_{\rm crit}$) of torus instability onset, the overlying field strength ($B_{\rm t}$), and ribbon flux ($R_{\rm f}$), provide physically motivated measures of eruption onset. The two main aims of this work are to (i) show that $h_{\rm crit}$ and $B_{\rm t}$ can equally well predict CME speeds when evaluated over the region of interest (ROI) not directly above the PIL, and (ii) assess the value of $h_{\rm crit}$, $B_{\rm t}$ and $R_{\rm f}$ in predicting CME speed. Photospheric magnetograms are modeled with potential-field extrapolations to obtain decay index profiles. Critical heights above PILs correlate strongly with 3D CME speed ($r = 0.71$). Using ROIs of $\approx$ 1.8, 3.7, and 7.3 Mm), centered on the PIL, weighted $h_{\rm crit}$ from the 7.3x7.3 ROI provides the strongest correlation ($r = 0.73$), while mean $B_{\rm t}$ at 150 Mm is weaker ($r = 0.33$). Combining both offers little improvement ($r = 0.74$), confirming $h_{\rm crit}$ as the dominant predictor. CME speed correlates moderately with $B_{\rm t} \times R_{\rm f}$ ($r = 0.44$), and highest when combined with $h_{\rm crit}$ ($r = 0.76$). Thus, in potential field models, ROI-based critical heights are as predictive as those above the PIL, indicating that the broader active-region field structure is equally valid as a diagnostic. When all parameters are considered together, $h_{\rm crit}$ alone consistently shows the highest predictive power for CME speed.

We present the photometric and spectroscopic analysis of the short-duration GRB 250221A ($T_{90}=1.80\pm0.32$ s), using a data set from the optical facilities COLIBRÍ, the Harlingten 50 cm Telescope, and the Very Large Telescope. We complement these observations with data from the \textit{Neil Gehrels Swift Observatory} and the \textit{Einstein Probe}, as well as radio observations from the Very Large Array. GRB 250221A is among the few short GRBs with direct afterglow spectroscopy, which gives a secure redshift determination of $z=0.768$ and allows the unambiguous identification of the host as a galaxy with a star-formation rate of $\sim3\,M_\odot\,{\rm yr}^{-1}$. The X-ray and optical light curves up to $T_0+10$ ks (where $T_0$ refers to the GRB trigger time) are well described by forward-shock synchrotron emission in the slow-cooling regime within the standard fireball framework. However, at $T_0+0.6$ days, both the X-ray and optical bands exhibit an excess over the same interval, which we interpret as evidence of energy injection into a jet with a half-opening angle of $\theta_j=11.5^{\circ}$ through a refreshed shock powered by late central engine activity or a radially stratified ejecta. The burst properties (duration, spectral hardness, peak energy, and location in the Amati plane) all favour a compact binary merger origin. However, our modelling of the afterglow suggests a dense circumburst medium ($n\sim80$ cm$^{-3}$), which is more typical of a Collapsar environment. This tension over the classification of this burst (short-hard vs. long-soft) as inferred from the prompt and afterglow emissions makes GRB~250221A an unusual event and underscores the limitations of duration-based classifications and the importance of multi-wavelength, time-resolved follow-up observations.

Eighteen years after their discovery, the astronomical sources and radiation mechanisms of fast radio bursts remain mysterious. Their radiation is as bright as that of pulsars, with brightness temperatures as high as $\sim 10^{36}$ K, implying coherent emission, but the plasma physics that forms the coherent charge bunches, with net charges of order a Coulomb, is not understood. Some FRB have been identified with galaxies at redshifts of a few tenths, but one originated within a globular cluster in the galaxy M81 at a distance of 3.6 Mpc. A minority of FRB have been observed to repeat, in some cases thousands of times. The vast majority of FRB have not been observed to repeat, but it is not known if they are truly ``one-offs'' or repeat at unobservably long intervals. Some FRB originate within dense, rapidly varying, plasma environments, while others appear to be surrounded by high vacuum. Hypotheses for their sources include magnetars and black hole accretion discs.

The sound horizon scale $r_s$ is a key source of information for early-time $H_0$ measurements, and is therefore a common target of new physics proposed to solve the Hubble tension. We present a sub-2% measurement of the Hubble constant that is independent of this scale, using data from the first data release of the Dark Energy Spectroscopic Instrument (DESI DR1). Building on previous work, we remove dependency on the sound horizon size using a heuristic rescaling procedure at the power spectrum level. A key innovation is the inclusion of \emph{uncalibrated} (agnostic to $r_s$) post-reconstruction BAO measurements from DESI DR1, as well as using the CMB acoustic scale $\theta_*$ as a high-redshift anchor. Uncalibrated type-Ia supernovae are often included as an independent source of $\Omega_m$ information; here we demonstrate the robustness of our results by additionally considering two supernova-independent alternative datasets. We find somewhat higher values of $H_0$ relative to our previous work: $69.2^{+1.3}_{-1.4}$, $70.3^{+1.4}_{-1.2}$, and $69.6^{+1.3}_{-1.8}\,{\rm km\,s^{-1}\,Mpc^{-1}}$ respectively when including measurements from i) Planck/ACT CMB lensing $\times$ unWISE galaxies, ii) the DES Year 3 6$\times$2pt analysis, and iii) Planck/ACT CMB lensing + the DES Year 5 supernova analysis. These remarkably consistent constraints achieve better than 2% precision; they are among the most stringent sound horizon-independent measurements from LSS to date, and provide a powerful avenue for probing the origin of the Hubble tension.

Dipping neutron star low-mass X-ray binaries (NS LMXBs) are systems that exhibit periodic drops in their X-ray light curves. These are believed to be caused by material at the impact point of the gas stream onto the accretion disk, the bulge. Dipping systems are observed at high inclination and provide exceptional opportunities to address important open questions about accretion disks, such as the physical properties of the bulge, and the connection between disk atmospheres and disk winds. We aimed to characterize the accretion disk plasmas present in the 21h-period NS LMXB 4U 1624-490, and perform a detailed spectral analysis of the material present at the impact region. We used four XMM EPIC pn observations that were specifically targeting dips, and allow us to probe dipping activity over different timescales (i.e. consecutive orbits and $\sim$6 months). We use both time- and flux-resolved spectroscopic analysis to probe the structural properties of the bulge moving along the line of sight and its homogeneity, respectively. During dipping, the primary spectrum is modulated by an ionized (log$\xi\sim$ 3.4) absorber with varying column density and covering factor, as well as a colder absorber. This suggests that the bulge is a multiphase and clumpy absorbing medium. From size scale arguments, we estimate the number of clumps in the bulge to be $>$7$\times10^{3}$. A highly ionized disk atmosphere becomes evident only when different phases of absorption are analyzed individually. This work demonstrates the feasibility of constructing a physical picture of the bulge, and highlights how future research could reveal how its properties depend on system parameters, and whether the bulge could influence the dynamics of the accretion disk.

Constructing a general-purpose framework for mapping between dark matter simulations and observable hydrodynamical simulation outputs is a long-standing problem in modern astrophysics. In this work, we present a new approach utilizing stochastic interpolants to map between cheap fast particle mesh simulations and baryonic quantities in three dimensions, requiring a total of 7 GPU minutes per 256^3 grid size simulation. Using the CAMELS multifield dataset, we are able to condition our mapping on both cosmological and astrophysical properties. We focus this work on hydrodynamical quantities suitable for Lya observables finding excellent agreement up to small spatial scales, k ~ 10.0 (h^(-1) Mpc) at z=2.0, for Lya flux statistics. Our approach is fully convolutional, allowing training on comparatively small volumes and application to larger volumes, which was tested on TNG50.

We propose a two-field model where the inflaton $\chi$ is non-minimally coupled to the instanton $\phi$. By choosing an appropriate coupling function, we realize the scenario where the difference of the values of potential between false vacuum (FV) and true vacuum (TV) is maximized during inflation. Most of the bubbles are created at this time. After inflation ends, the potential value of FV drops below that of TV so that these bubbles collapse to form primordial black holes (PBHs). By tuning the parameters of our model, we analyze the Coleman-de Luccia (CDL) and Hawking-Moss (HM) process, finding that the corresponding mass function of PBHs is sharply peaked, implying that we can realize either PBHs as cold dark matter, sub-solar PBHs, or supermassive PBHs in this scenario without enhancement of primordial curvature perturbations.

SMSS J0521-4351 is reportedly the most luminous quasar known to date, and assuming a mean radiative efficiency of 0.1, it is inferred to be the fastest-growing black hole, accreting approximately one solar mass per day. Assessing the implications of this assumption on the seed mass and inception time of J0529-4351, we show that the inferred accretion rate is unreasonably high and that its radiative efficiency must be much greater than 0.1. Then, we derive its accretion rate and seed mass, and for comparison of three other similar-size (~1-2E+10 solar masses) black holes at various redshifts, using well-tested empirical scaling relations. The results indicate that J0529-4351 grew from a heavy seed (~2-3E+04 solar masses), and that its accretion rate (~10-13 solar masses/year) is the lowest of the four black holes. However, its radiative efficiency inferred from its bolometric luminosity and the derived accretion rate is the highest, which explains why it is the most luminous despite having the lowest accretion rate. This study challenges the prevailing notion that a higher luminosity or a higher Eddington ratio implies a higher accretion rate, highlights the dependence of a black hole's luminosity on radiative efficiency, reveals the pitfalls of inferring black hole properties assuming a standard value for radiative efficiency, and suggests that the Eddington ratios of high-luminosity BHs may be significantly overestimated.

Accurate identification of unobservable regions in nighttime is essential for autonomous scheduling and data quality control in this http URL methods-such as infrared sensing or photometric extinction-provide only coarse,non-spatial estimates of sky clarity,making them insufficient for real-time this http URL not only wastes observing time but also introduces contamination when telescopes are directed toward cloud-covered or moonlight-affected this http URL address these limitations,we propose a deep learning-based segmentation framework that provides pixel-level masks of unobservable areas using all-sky this http URL by a manually annotated dataset of nighttime images,our method enables precise detection of cloud- and moonlight-affected this http URL segmentation results are further mapped to celestial coordinates through Zenithal Equal-Area projection,allowing seamless integration with observation control systems (OCS) for real-time cloud-aware this http URL developed for the Mephisto telescope,the framework is generalizable and applicable to other wide-field robotic observatories equipped with all-sky monitoring.

We report wideband scatter-broadening estimates of 14 pulsars towards the Gum nebula region using the Band-3 of the upgraded GMRT. This work increases the measurements of frequency scaling index of scatter-broadening ($\alpha$) across the nebula by more than 3 times. A strong correlation between the distance and the scattering strength is observed for pulsars behind the nebula. It is also observed that for distant pulsars ($> 2 kpc$), the effect of the Gum nebula in DM and scattering strength is not substantial. We also report a much flatter $\alpha$ for the Vela pulsar and argue that its scattering is not caused by the Gum nebula, but the Vela supernova remnant.

We present a Bayesian latent model to describe the scaling relation between globular cluster populations and their host galaxies, updating the framework proposed in de Souza 2015. GC counts are drawn from a negative-binomial (NB) process linked to host stellar mass, augmented with a newly introduced Gaussian observation layer that enables efficient propagation of measurement errors. The revised formulation preserves the underlying NB process while improving computational tractability. The code snippets, implemented in Nimble and PyMC are released under the MIT license at this https URL

In this thesis, we introduce novel methods for analyzing pulsar populations using a variety of mathematical techniques. These tools-particularly graph theory-have been thoroughly validated in advanced mathematics, enabling us to overcome some of the constraints (even dimensional) inherent in conventional visualization approaches. This exploration benefits from dimensionality reduction techniques, which not only lessen computational demands but also highlight potential for describing physical characteristics. The resulting structures encode information about pulsar similarities that extend beyond standard spin parameters, revealing relationships that are not readily apparent in traditional diagrams. With a physically motivated topological perspective, we leverage the strengths of these methods and present results that span from prospective source classification and the emergence of new classes to catalog comparison, among other applications. This new approach enables fresh interpretations of longstanding problems, laying a new foundation for visualizing the pulsar population and categorizing sources. Building on this, we identify several sources as likely members of specific binary subclasses and investigate the potential transitional nature of others. Furthermore, we extend the use of graph theory to the boundary of machine learning, demonstrating its capability for binary separation in an unsupervised context. Finally, we introduce and apply an innovative, flexible time-series alignment technique to the field of gamma-ray astrophysics. The method identifies notable similarities among the light curves of gamma-ray pulsars. The results presented here are promising, offering a refreshing direction for the field and new pathways for rigorous mathematical analysis, ultimately providing meaningful alternatives to traditional approaches in high-energy astrophysics.

This paper explores how time-varying increases in mass accretion onto rapidly spinning black holes influence their long-term spin evolution when affected by superradiance - a process where energy is extracted from the black hole by a surrounding axion field. Using simulations the study tracks how sudden accretion boosts affect a critical spin-down phase (the superradiance drop) during which the black hole's spin rapidly decreases while its mass remains nearly constant. The black hole spin evolution is controlled by the competition between two processes: how fast angular momentum is added through accretion, and how fast it is removed by the axion cloud. One major conclusion is that boosts to the accretion rate before the superradiance drop have the strongest effect, as they can delay or reshape the drop and significantly shrink the region of the mass-spin plane depopulated due to the superradiance. In particular, a super-Eddington accretion rate of 5 times Eddington accretion, lasting for 4 Myr and occurring 30 Myr before the superradiance drop can reduce the superradiance exclusion region in the mass-spin plane by 40 percent. In contrast, boosts to the accretion rate after the superradiance drop only cause temporary changes in the black hole spin. The study also shows that black holes with lighter axion clouds are more sensitive to these early boosts and can show observable spin changes lasting tens to hundreds of millions of years. Heavier axion clouds, however, require much stronger or longer-lasting boosts to produce similar effects, making them more stable under variable accretion.

Early interaction of supernova blast waves with CSM has the potential to accelerate particles to PeV energies, although this has not yet been detected. Current models for this interaction assume the shock expands into a smooth stellar wind, although observations of many SNe do not support this assumption. We extend previous work by considering shocks expanding into complex density profiles consisting of smooth winds with dense CSM shells at various distances from the progenitor star. We aim to predict the gamma-ray and multiwavelength signatures of CSM interaction. We used the PION code to model the CSM around LBV including a brief episode of enhanced mass-loss and to simulate the formation of photoionization-confined shells around RSGs. Consequently, we used the time-dependent acceleration-code RATPaC to study the acceleration of cosmic rays in SNe expanding into these media and to evaluate the emitted radiation across the whole electromagnetic spectrum. We find that the interaction with the CSM shells can significantly boost the gamma-ray emission, with the emission peaking weeks to years after the explosion. The peak luminosity for Type-IIP and Type-IIn remnants can exceed the luminosity expected for smooth winds by orders of magnitude. For Type-IIP explosions, the light-curve peak is only reached years after the explosion. We evaluate the multiwavelength signatures expected from the interaction of the blast wave with a dense CSM shell from radio, over optical, to thermal X-rays. We identify high-cadence optical surveys and continuous monitoring of nearby SN in radio and mm wavelengths as the best-suited strategies for identifying targets that should be followed-up by gamma-ray observatories. We predict that gamma-rays from interaction with dense CSM shells may be detectable out to a few Mpc for late interaction, and tens of Mpc for early interaction.

Using the Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) we performed an untargeted search for H I-bearing ultra-diffuse galaxies (UDGs). We identified a core sample of 10 UDGs defined by $\mu_{g,0}\ge24$ mag arcsec$^{-2}$ and $R_{e}\ge1.5$ kpc, and a broader sample including 12 additional faint diffuse galaxies ($\mu_{g,0}\ge23.7$ mag arcsec$^{-2}$ and $R_{e}\ge1.3$ kpc). Within the core sample, we highlight the first discovery of a UDG pair. Their projected separation is just 75 arcsec (22 kpc at 61.9 Mpc), with a central H I velocity difference of 34 km s$^{-1}$. The North-Western UDG (WALLABY J104513-262755-UDG-1) has a larger H I reservoir, $\log_{10}(M_{HI}/\rm M_{\odot}) = 8.95\pm0.03$, compared to the South-Eastern UDG (WALLABY J104513-262755-UDG-2), $\log_{10}(M_{HI}/\rm M_{\odot}) = 8.60\pm0.04$. UDG-1's stellar mass and star formation rate are also approximately an order of magnitude larger at $\log_{10}(M_*/\rm M_{\odot}) = 8.07\pm0.12$ and $\log_{10}(SFR/\rm M_{\odot}~yr^{-1}) = -1.26\pm0.12$ respectively. The pair has an isolated local environment, with no other galaxies or H I sources within 30 arcmin (525 kpc) and $\pm1000$ km s$^{-1}$. However, in the context of the larger-scale structure, the pair is located outside the virial radius of the Hydra cluster, with its position on the phase-space diagram indicating that it is infalling into the cluster. The identification of this H I-bearing UDG pair raises important questions around the formation of such a unique system and the evolution of UDGs in a transitional phase before ram pressure stripping and cluster infall.

Numerous astronomical and cosmological observations point to the existence of dark matter, which constitutes about 27% of the Universe. Despite extensive efforts, only the DAMA/LIBRA experiment, using NaI(Tl) detectors at Gran Sasso National Laboratory, has reported a positive dark matter signal. To independently verify this result, using the same NaI target is essential. This is the goal of ANAIS-112, operating with 112.5 kg of NaI(Tl) scintillators since 2017 at the Canfranc Underground Laboratory (LSC). This thesis presents work within ANAIS-112, focusing on data analysis and the development of Geant4-based simulations to reduce systematics and increase sensitivity. Efforts have focused on understanding the ANAIS-112 crystals' response to nuclear recoils. Systematic uncertainties related to the scintillation quenching factors (QFs) of sodium and iodine recoiling in NaI(Tl) must be addressed to ensure reliable comparisons with DAMA/LIBRA and other experiments. For that purpose, onsite neutron calibrations in ANAIS-112 have been performed since 2021 at LSC using 252Cf sources. This work has compared calibration data with Geant4-based neutron simulations, which has revealed Geant4 deficiencies in certain decay processes and cross sections in some versions. The simulations have proven highly sensitive to the QF used, favouring models in which the QF increases with energy and disfavouring DAMA values. In parallel, this thesis has contributed to improving the ANAIS-112 background model through a multiparametric fit of its background components. New physics searches have been conducted using ANAIS-112 data, including a reanalysis of the annual modulation signal and a search for solar axions using six years of exposure. Finally, this thesis also includes work within the COSINUS experiment, focused on background modelling and the impact of internal backgrounds on the experiment's sensitivity.

In the vicinity of neutron star mergers (NSMs), it is possible for the neutrino self-interaction potential to cancel with the matter potential leading to matter neutrino resonance (MNR). MNR is one of the most interesting mechanisms by which neutrino flavor evolution can occur in dense astrophysical environments. Previous studies have typically assumed that the neutrino flavor field evolves to a steady state -- a simplification also used in other self-interaction models such as the neutrino-bulb model. Here, we perform reproducible calculations of MNR using both time-independent and time-dependent formalisms and show that they yield qualitatively different flavor survival probabilities. The time-independent approach produces unstable steady-state solutions that differ fundamentally from the dynamical behavior captured in time-dependent simulations. These results demonstrate that the steady-state assumption is generally invalid, and physical interpretations based on time-independent calculations of dense neutrino systems require re-evaluation.

Jets from active galactic nuclei (AGNs) are expected to heat the surrounding intracluster medium (ICM). We investigate how the interaction between jets and the ICM appears in high-resolution X-ray observations using mock X-ray observations based on two-dimensional hydrodynamic simulations. We constructed a model of an active galactic nucleus (AGN) similar to Cygnus A (Cyg A), a powerful FR II radio galaxy. Our simulations model bipolar jets propagating into a stratified ICM, forming forward shocks and low-density cocoons. Based on these results, we generate synthetic spectra that incorporate both shocked and unshocked ICM components. Then, we perform mock observations using the XRISM/Resolve X-ray spectrometer. We focus particularly on viewing angle effects. Our mock observations revealed that the smallest line broadening, observed as velocity dispersion, associated with the cocoon's bulk expansion occurs when observing along the jet direction, where the expansion velocity is highest. Although this may appear counterintuitive, it occurs because the rapidly expanding jet head contributes little to X-ray emission due to its high temperature and low density. Our results highlight the importance of considering the temperature and density structure of AGN-driven shocks and cocoons when interpreting XRISM data. These findings lay the groundwork for XRISM's observations of AGN jets and will improve our understanding of AGN feedback processes in galaxy clusters.

Powerful ionized accretion disk winds are often observed during episodic outbursts in Galactic black hole transients. Among those X-ray absorbers, \fexxvi\ doublet structure (Ly$\alpha_1$+Ly$\alpha_2$ with $\sim 20$eV apart) has a unique potential to better probe the underlying physical nature of the wind; i.e. density and kinematics. We demonstrate, based on a physically-motivated magnetic disk wind scenario of a stratified structure in density and velocity, that the doublet line profile can be effectively utilized as a diagnostics to measure wind density and associated velocity dispersion (due to thermal turbulence and/or dynamical shear motion in winds). Our simulated doublet spectra with post-process radiative transfer calculations indicate that the profile can be (1) broad with a single peak for higher velocity dispersion ($\gsim 5,000$ km~s$^{-1}$), (2) a standard shape with 1:2 canonical flux ratio for moderate dispersion ($\sim 1,000-5,000$ km~s$^{-1}$) or (3) double-peaked with its flux ratio approaching 1:1 for lower velocity dispersion ($\lsim 1,000$ km~s$^{-1}$) in optically-thin regime, allowing various line shape. Such a diversity in doublet profile is indeed unambiguously seen in recent observations with XRISM/Resolve at microcalorimeter resolution. We show that some implications inferred from the model will help constrain the local wind physics where \fexxvi\ is predominantly produced in a large-scale, stratified wind.

Atmospheric tomography, the problem of reconstructing atmospheric turbulence profiles from wavefront sensor measurements, is an integral part of many adaptive optics systems. It is used to enhance the image quality of ground-based telescopes, such as for the Multiconjugate Adaptive Optics Relay For ELT Observations (MORFEO) instrument on the Extremely Large Telescope (ELT). To solve this problem, a singular-value decomposition (SVD) based approach has been proposed before. In this paper, we focus on the numerical implementation of the SVD-based Atmospheric Tomography with Fourier Domain Regularization Algorithm (SAFR) and its performance for Multi-Conjugate Adaptive Optics (MCAO) systems. The key features of the SAFR algorithm are the utilization of the FFT and the pre-computation of computationally demanding parts. Together, this yields a fast algorithm with less memory requirements than commonly used Matrix Vector Multiplication (MVM) approaches. We evaluate the performance of SAFR regarding reconstruction quality and computational expense in numerical experiments using the simulation environment COMPASS, in which we use an MCAO setup resembling the physical parameters of the MORFEO instrument of the ELT.

The $^{63}$Ga(p,$\gamma$)$^{64}$Ge and $^{64}$Ge(p,$\gamma$)$^{65}$As thermonuclear reactions connect the ZnGa and GeAs cycles by diverting the flow of the rapid proton capture process from $^{63}$Ga to $^{65}$As. Changes in these two reaction rates regulate the ZnGa and GeAs cycles and may affect the modeled properties matching with the observed counterparts of a type I X-ray burster. We implement the latest $^{63}$Ga(p,$\gamma$)$^{64}$Ge and $^{64}$Ge(p,$\gamma$)$^{65}$As reaction rates to the state-of-the-art self-consistent one-dimensional multi-zone thermo-hydrodynamic code, KEPLER, to study the influence of these new reaction rates on the models of the GS 1826$-$24 clocked burster and SAX J1808.4$-$3658 photospheric radius expansion burster. Both new reaction rates obtained by Lu et al. [Phys. Rev. C 110, 065804 (2024)] are determined from complementing the experimental input with the nuclear spectroscopic information deduced from the full pf-shell space configuration-interaction shell-model calculations. By constraining the models on reproducing the observed burst peak, light-curve profile, fluence, and recurrence time, we find that the impact of the newly measured proton thresholds and respective proton-capture reactions on the burst light-curve profile of the GS 1826$-$24 clocked burster is, in fact, not as significant as claimed by Zhou et al. [Nat. Phys. 19, 1091 (2023)]. With or without the inclusion of the newly determined reaction rate of the highly influential $^{22}$Mg($\alpha$,p)$^{25}$Al reaction, the impact of the new $^{63}$Ga(p,$\gamma$)$^{64}$Ge and $^{64}$Ge(p,$\gamma$)$^{65}$As reaction rates on SAX J1808.4$-$3658 photospheric radius expansion bursts is evident. Our finding indicates that the models reproducing the 2002 October epoch of SAX J1808.4$-$3658 photospheric radius expansion burster is more sensitive to the uncertainties of thermonuclear reaction rates.

Accretion in black hole X-ray binaries is commonly believed to be supplied by the Roche lobe overflow or the stellar wind. The former is thought to form a geometrically thin disc while the diffuse wind could form a geometrically thick hot accretion flow. In this paper, we instead consider a more generalised case, i.e., accretion with both cold and hot gas supplies, which feed a disc and a corona respectively. We investigate the interaction of disc and corona by analysing the energy coupling and matter exchange, i.e. corona condensation/disc evaporation, with a semi-analytical method. It is found that the accretion geometry in the radial direction and the resultant emission spectrum depend strongly on both the total gas supply rate and the ratio of cold and hot gases. For gas supply rates of a few percent of the Eddington value, diverse geometries and spectral shapes are possible, depending on the fraction of cold gas supply. This provides an interpretation for the various spectra observed in intermediate states. However, at higher accretion rates, regardless of the form of the feeding gas, the inner accretion flow is always disc-dominated, implying an inevitable transition to the soft state, while at very low gas supply rates, hard state spectrum dominated by the hot flow is expected. We also present the predicted hardness-intensity correlation of Cygnus X-1, and constrain the value of the viscosity parameter of the accretion flow to the range of 0.25--0.35 by comparing our results with MAXI observations.

Launching in 2027 and 2029, respectively, Twinkle and Ariel will conduct the first large-scale homogeneous spectroscopic surveys of the atmospheres of hundreds of diverse exoplanets. This will fundamentally transition the field to an era of population-level characterisation. In this pilot study, we aim to explore possible synergies between Twinkle and Ariel to determine for instance whether prior Twinkle observations can substantially inform the target selection and observing strategy of Ariel. This study primarily aims to encourage further investigation by both consortium communities by showing what a potential scientific synergy would look like on a promising scientific case that requires further exploration. For this aim, we select a small subset of "cool" planets that are also particularly well-suited to be observed by Twinkle and therefore Ariel. By using representative noise estimates for both missions, we compute the number of visits required for an observation. Then, we simulate and retrieve transmission spectra of each target, assuming gaseous, H2/He-dominated atmospheres and various atmospheric models. For all candidates, we find that atmospheric parameters are generally retrieved well within 1-sigma to input values, with Ariel typically achieving tighter constraints. We demonstrate that for a small subset of cool gaseous planets, exploitable synergies exist between Twinkle and Ariel observations and Twinkle may very well provide a vantage point to plan Ariel observations. The true extent of the potential synergies, far beyond our considered sample, will be determined by the final target lists. Once Twinkle is operational and its performance is known, it could reliably inform Ariel's target prioritization and Ariel's capabilities which are already well-established can help define optimal targets and observational approaches for Twinkle.

We show that there is a strong dependence of the radio LDF electric field amplitudes at ground level on the position of $X_{\rm max}$ in the atmosphere, even accounting for differences in the EM energy of the showers. Since an $X_{\rm max}$ dependence leads to a primary composition dependence, this implies that information on the mass composition is encoded not only in the LDF shape but also in its amplitude. This $X_{\rm max}$ dependence can be explained in terms of two competing scalings of the measured electric field: One goes with $(1/\rho)^J$, where $\rho$ is the air density at $X_{\rm max}$ and $J$ is a zenith dependent non-linearity factor describing coherence loss. This density scaling tends to decrease the geomagnetic emission of deeper showers. The other scaling goes with $(1/R)$, where $R$ is the distance from $X_{\rm max}$ to the core at ground, and instead increases the measured electric field of deeper showers. At low zenith angles, the $(1/R)$ scaling is stronger and leads to larger measured electric fields as $X_{\rm max}$ increases. The picture at higher zeniths, i.e., lower densities, is more nuanced. In this region, the deflections due to the Lorentz force are much larger and introduce extra time delays between the particle tracks, decreasing the coherence of the emission. This loss of coherence is highly dependent on the strength of the geomagnetic field and can slow down, or even reverse the increase of the radio emission with decreasing air density. This strong, yet historically overlooked LDF amplitude dependence on $X_{\rm max}$/composition could be used to directly infer, even bypassing any $X_{\rm max}$ reconstruction, the cosmic ray primary composition on an event-by-event basis. It could also have some repercussions on other radio reconstruction methods, such as a possible $X_{\rm max}$/composition bias on shower electromagnetic energy reconstruction methods.

Time series analysis is fundamental to characterizing the variability inherent in multi-wavelength emissions from blazars. However, a major observational challenge lies in the need for well-sampled, temporally uniform data, which is often hindered by irregular sampling and data gaps. These gaps can significantly affect the reliability and accuracy of methods used to probe source variability. This paper investigates the impact of such observational gaps on time series analysis of blazar emissions. To do so, we systematically evaluate how these gaps alter observed variability patterns, mask genuine periodic signals, and introduce false periodicity detections. This evaluation is conducted using both simulated and real observational data. We assess a range of widely used time series analysis methods, including the Lomb-Scargle periodogram, Phase Dispersion Minimization, and the recently proposed Singular Spectrum Analysis (SSA). Our results demonstrate a clear and significant degradation in period detection reliability when the percentage of gaps exceeds 50\%. In such cases, the period-significance relationship becomes increasingly distorted, often leading to misleading results. Among the tested methods, SSA stands out for its ability to yield consistent and robust detections despite high data incompleteness. Additionally, the analyzed methods tend to identify artificial periodicities of around one year, likely due to seasonal sampling effects, which can result in false positives if not carefully recognized. Finally, the periods detected with $\geq$3$\sigma$ confidence are unlikely to result from stochastic processes or from the presence of gaps in the analyzed time series.

The Standard Model of particle physics successfully describes all known fundamental particles and their interactions; however, it leaves several unanswered questions. Theories beyond the Standard Model typically introduce new particles and symmetries to address these issues. In the early universe, when such particles become non-relativistic, or the symmetries are broken, there are associated reductions in the equation of state of the primordial plasma. These reductions lead to an exponential enhancement in the formation rate of primordial black holes. In this paper, we calculate the equation of state for several supersymmetric and composite Higgs models, which naturally predict a large number of additional degrees of freedom. Using these equations of state, we compute some example primordial black hole abundances, which we find can be enhanced by up to 20 orders of magnitude.

We present high-resolution (~7.5 au) ALMA observations at 1.3 and 3 mm of 16 disks around Class 0/I protostars across multiple star-forming regions and a variety of multiplicities, showing a range of disk sizes (~2-100 au) and including circumbinary disks (CBDs) in binaries with separations <100 au. The disk properties show similarities to Class II disks, including (a) low spectral index (SI) values (alpha=2.1) that increase with disk radius, (b) 3 mm disk sizes only marginally smaller than at 1.3 mm (<10%), and (c) radial intensity profiles well described by modified self-similar profiles. We also find key differences: (i) SI values increasing with radius, but exceeding 2 only at the disk edge (ii) higher brightness temperatures Tb, in some cases higher than the predicted temperatures due to irradiation, and (iii) ~10x higher luminosity at a given size compared to the Class II disks. These results confirm significant optical depth in the observed Class 0/I disks, at both 1.3 and 3 mm, helping to explain their higher luminosities, but higher temperatures are also required for the most compact (< 40 au) disks, suggesting additional viscous heating. Considering optical depth, most disk dust masses are estimated in the range 30-900 Mearth (0.01-0.3 Msun in gas), resulting in some disks reaching marginal gravitational instability. The median location of the water iceline is ~3 au, but it can extend beyond 10-20 au for the hottest disks. CBDs exhibit lower optical depths at both wavelengths and hence higher SI values (alpha=3.0), dust masses of 100 Mearth, and beta~1.5 (2 Class 0 CBDs) and beta~1 (1 Class I CBD), suggesting substantial grain growth only in the more evolved CBD. The inferred high optical depths provide a compelling explanation for the apparent scarcity of dust substructures in the younger disks at ~ 1 mm, despite mounting evidence for early planet formation (ABRIDGED).

The Population III.1 theory for supermassive black hole (SMBH) formation predicts a very early ($z\sim20-25$), transient phase, ``The Flash'', of cosmic reionization powered by supermassive stars that are SMBH progenitors. The universe then quickly recombined to become mostly neutral, with this state persisting until galaxies begin to reionize intergalactic gas again at $z\sim 10$. The overall Thomson scattering optical depth, $\tau$, from The Flash has been shown to be $\tau_{\rm PopIII.1}\sim0.03$, leading to a total $\tau\sim0.08-0.09$. Such a value, while significantly larger than that previously inferred from {\it Planck} observations of the low-$l$ $EE$ polarization power spectrum of the CMB, can help relieve several ``tensions'' faced by the standard $\Lambda$CDM cosmological model, especially the preference for negative neutrino masses and dynamic dark energy. Here we compute $EE$ power spectra of example models of The Flash. We find that, because of its very high redshift, the contribution to $l\lesssim8$ modes is dramatically reduced compared to usual low-$z$ reionization models for the same value of $\tau$, while the power at $l\gtrsim8$ is boosted. Thus the Pop III.1 reionization scenario provides a natural way to increase $\tau$, while remaining closer to the latest CMB low-$l$ polarization observations.

The origin of TeV-PeV neutrinos detected by IceCube remains largely unknown. The most significant individual neutrino source is the close-by Seyfert galaxy NGC 1068 at 4.2$\sigma$ level with a soft spectral index. Another notable candidate is the Seyfert galaxy NGC 7469, which has been recently proposed as a potential neutrino emitter. The likelihood fit of the IceCube data for this source returned a very hard spectral index of ~ 1.9 and the excess is dominated by two high-energy events, issued as two neutrino alerts IC220424A and IC230416A. The energies of the two neutrinos are estimated to be 100-200 TeV, implying a maximum proton energy > 2 PeV, significantly higher than that in NGC 1068. The lack of lower-energy neutrinos from NGC 7469 also suggests a neutrino spectrum harder than that of NGC 1068. In this paper, we analyze the Fermi-LAT observations of NGC 7469, which yield non-detection. By requiring the cascade flux accompanying neutrino production not to exceed the upper limit of the GeV flux, the size of the neutrino-emitting region can be constrained when the neutrino flux takes a high value of the allowed range. We suggest that protons are accelerated to PeV energies via turbulence or magnetic reconnection in the corona of NGC 7469 and interact with OUV photons from the accretion disk and X-rays from the corona through the $p\gamma$ process, producing neutrinos with energy of 100-200 TeV. In the turbulence acceleration scenario, the required maximum proton energy can be achieved with a magnetization parameter close to unity ($\sigma\sim 1$), while in the reconnection scenario, a magnetization parameter with $\sigma\sim 10$ is needed. In both scenarios, a pair dominated composition for the corona is preferred. The difference in the neutrino spectrum between NGC 7469 and NGC 1068 could be due to a different magnetization despite that they belong to the same type of AGN.

We investigate whether combining gas and stellar kinematic maps provides measurable advantages in recovering galaxy mass profiles, compared to using single-component maps alone. While traditional methods struggle to integrate multi-tracer data effectively, we test whether deep learning models can leverage this joint information. We develop a probabilistic convolutional neural network (CNN) framework trained and tested on mock galaxy kinematic maps from multiple cosmological simulation suites. Our model is trained on gas-only, stars-only, and combined gas+stellar velocity maps, allowing direct comparison of performance across tracers. To assess robustness, we include simulations with differing feedback models and galaxy properties. Combining gas and stellar maps reduces the dispersion in the inferred mass profiles by up to a factor of $\sim$1.5 compared to models using either tracer independently. The CNN architecture effectively captures complementary information from the two components. However, we find limitations in generalizing between simulation suites, with reduced performance when applying models trained on one suite to galaxies from another.

This perspective offers a viewpoint on how the challenges of molecular scattering investigations of astrophysical interest have evolved in recent years. Computational progress has steadily expanded collisional databases and provided essential tools for modeling non-LTE astronomical regions. However, the observational leap enabled by the JWST and new observational facilities has revealed critical gaps in these databases. In this framework, two major frontiers emerge: the characterization of collisional processes involving heavy projectiles, and the treatment of ro-vibrational excitation. The significant computational effort of these investigations emphasizes the need to test and develop robust theoretical methods and approximations, capable of extending the census of collisional coefficients required for reliable astrophysical modeling. Recent developments in these directions are outlined, with particular attention to their application and their potential to broaden the coverage of molecular systems and physical environments.

The analysis of the composition of primitive C$-$complex asteroids is essential to understand the distribution of volatiles in the Solar System since its formation. Primitive low-albedo families within the inner main asteroid belt are of particular interest because they are a significant source of carbonaceous near-Earth asteroids, such as Ryugu and Bennu. This study, part of the JWST SAMBA3 project (Spectral Analysis of Main Belt Asteroids in the 3 $\mu$m region), report the first spectroscopic analysis of asteroid (84) Klio in the 3 $\mu$m region, in order to better constrain its composition. We analysed the infrared (0.97$-$5.10 $\mu$m) Spectrum of Klio measured by the NIRSpec instrument on board JWST. We used the NEATM thermal model to extract the reflectance spectrum of the asteroid. Several spectral features were then analysed in the 2.8, 3.4, and 3.9 $\mu$m regions by different Gaussian fitting. The Spectrum of Klio shows an absorption band at 2.776 $\pm$ 0.001 $\mu$m that we attributed to phyllosilicates. We compared the position and shape of the feature with that observed in primitive materials such as carbonaceous chondrites and returned samples from Ryugu and Bennu. The position and shape of the 2.8 $\mu$m band, as well as the presence of a 0.7 $\mu$m band in the visible, suggest that Klio's spectrum is similar to certain CM2 meteorites. We observed an absorption band around 3.9 $\mu$m, with a depth of $0.020 \pm 0.001$ that could be attributed to carbonates. We could not clearly detect any absorption associated with organics at 3.4 $\mu$m.

Aims: The polarization of the Ca II resonant doublet (H and K lines) and the subordinate infrared triplet lines are key observables for diagnosing solar chromospheric magnetism. It is thus necessary to understand the physical mechanisms that shape their Stokes profiles in magnetic environments. Methods: Using the spectral synthesis module of the HanleRT-TIC code, we study the effects of anisotropic radiation pumping with partial frequency redistribution (PRD) and J-state interference (JSI) in a plane-parallel semi-empirical static solar atmospheric model. We also analyze the sensitivity of these lines to magnetic fields of varying strengths and orientations, accounting for the combined action of the Hanle and Zeeman effects. Results: Including PRD is crucial to model the polarization in the core regions of the resonant lines, while JSI strongly affects their far wings. The metastable lower levels of the subordinate lines also influence the scattering polarization of the K line. With horizontal magnetic fields, the resonant lines respond to field strengths from sub-gauss to tens of gauss, whereas the infrared triplet scattering polarization is mainly sensitive to milligauss fields. At a near-limb line of sight (LOS) with $\mu = 0.1$, the Hanle effect modifies the scattering polarization via a depolarization and a rotation in the plane of linear polarization. At disk center, horizontal fields generate linear polarization in the 1D model: for the K line, the Hanle effect dominates from sub-gauss to a few tens of gauss, and the Zeeman effect dominates in stronger fields. For vertical fields, the Hanle effect vanishes, but magneto-optical effects affect the linear polarization wings. Finally, atomic level polarization impacts the outer circular polarization lobes of the resonant lines, and the weak-field approximation overestimates the LOS magnetic component in this frequency range.

We present an analysis of neutrino-driven magnetohydrodynamic (MHD) waves and instabilities in a rotating magnetoplasma with weak neutrino interactions. We show, for the first time, that neutrino-driven shear Alfv{é}n and oblique magnetosonic waves can be coupled by the Coriolis force, forming new wave modes affected by this force, as well as neutrino beam and two neutrino flavor oscillations. Our work extends previous theories by demonstrating that shear Alfv{é}n waves are influenced by neutrino effects and by identifying instabilities resulting from resonant interactions with both a streaming neutrino beam and flavor oscillations. We find that the Coriolis force, as well as plasma density and magnetic field strength, have a significant impact on the profiles of the instability growth rates. Our findings may shed new light on the physical mechanisms underlying core-collapse supernovae.

Radio observations provide a window into a planet's interior and play a crucial role in studying its atmosphere and surface, key factors to find potential habitability. The discovery of thousands of exoplanets, together with advances in radio astronomy through the Square Kilometre Array (SKA), motivates the search for planetary-scale radio emissions. Here, we employ the radiometric Bode's law (RBL) and machine learning techniques to analyze a dataset of 1330 confirmed exoplanets, aiming to estimate their potential radio emission. Permutation Importance (PI) and SHapley Additive exPlanations (SHAP) analyses indicate that a planet's mass, radius, orbital semi-major axis, and distance from Earth are sufficient to dependably forecast its radio flux and frequency. The random forest model accurately reproduces these radio characteristics, confirming its reliability for exoplanetary radio predictions. Considering observational constraints, we find that 64 exoplanets could generate signals detectable by the SKA, 52 of which remain observable in the intermediate AA* deployment. Among these, MASCARA-1 b stands out with a predicted flux of 7.209 mJy at 135.1 MHz, making it an excellent SKA-Low target. Meanwhile, WASP-18 b, with a flux of 18.638 mJy peaking at 812.9 MHz, is the most promising candidate for SKA-Mid. These results show that the SKA can detect gas giants, such as MASCARA-1 b (SNR>400) and WASP-18 b (SNR>4236), within feasible integration times. Additionally, we identify four candidates (HATS-18 b, WASP-12 b, WASP-103 b, and WASP-121 b) that are likely affected by radio quenching, highlighting the importance of considering this effect in target selection for observation campaigns.

We explore the observational prospects for detecting gravitational lensing induced by cosmological matter currents, a relativistic correction to the standard scalar lensing effect arising from the motion of matter. We propose to isolate this contribution by cross-correlating the weak-lensing convergence field with a reconstructed cosmic momentum field inferred from galaxy surveys. Using numerical simulations, we demonstrate that this reconstructed momentum field is uncorrelated with the scalar lensing signal, enabling a clean separation of the gravitomagnetic component. We then forecast the detectability of this signal for upcoming wide-field galaxy and weak-lensing surveys, showing that a statistically significant detection may be achievable under realistic observational conditions. Such a measurement would provide the first direct probe of the large-scale cosmic momentum field, offering a novel test of general relativity and Lorentz invariance on cosmological scales.

Aims: We carried out a detailed investigation of Lithium and CNO abundances, including carbon isotope ratios, in RS CVn stars to assess the role of magnetic activity in the mixing of stellar atmospheres. Methods: We obtained high-resolution spectra at the Moletai Astronomical Observatory. Lithium abundances were determined by spectral synthesis of the 6707 A line and the CNO abundances using the C2 band heads at 5135 and 5635.5 A CN bands at 6470- 6490 A and 7980 to 8005 A, and the [O I] line at 6300 A. By fitting the 13CN band at 8004.7 A, we determined the carbon isotope this http URL. We determined the main atmospheric parameters and investigated the chemical composition of 32 RS CVn stars. Lithium abundances were determined for 13 additional stars using archival spectra. We report that *iot Gem and HD 179094 have carbon isotope ratios already affected by extra-mixing, even though they are in the evolutionary stage below the red giant branch luminosity bump. About half of the low-mass giants, for which the lithium abundance was determined, follow the first dredge-up predictions; however, other stars show reduced Lithium abundances, as predicted by thermohaline-induced mixing. The intermediate-mass stars show reduced Lithium abundances reduced, as predicted by rotation-induced mixing. Conclusions. In low-mass, chromospherically active RS CVn stars, extra-mixing of lithium and carbon isotopes may begin earlier than in normal giants. The Li-rich RS CVn giant V*OP And has large C/N and carbon isotope ratios and raises questions about the origin of its lithium enhancement.

We study the substructure, connectivity, and galaxy content of galaxy clusters A1656 and 1367 in the Coma supercluster and of A1185 in the Leo supercluster with the aim of understanding the evolution of clusters from turnaround to virialisation. We used data from the SDSS DR10 MAIN galaxy sample and from DESI cluster catalogues. The projected phase space diagram and the distribution of mass were used to identify regions of various infall stages (early and late infall, and regions of ongoing infall, i.e. regions of influence), their characteristic radii, embedded mass, and density contrasts in order to study the evolution of clusters with the spherical collapse model. We determined the substructure of clusters using normal mixture modelling and their connectivity by counting filaments in the cluster's regions of influence, analysed galaxy content of clusters, and derived scaling relations between cluster masses. All three clusters have a substructure with two to five components and up to six filaments connected to them. The radii of regions of influence are $R_\mathrm{inf} \approx 4$ Mpc, and the density contrast at their borders is $\Delta\rho_{inf} \approx 50 - 60$. The scaling relations between the masses of clusters have a very small scatter. The galaxy content of the clusters and of their regions of influence vary from cluster to cluster. In superclusters the percentage of quiescent galaxies is higher than in low-density regions between superclusters. The collapse of the regions of influence of clusters started at redshifts $z \approx 0.4 - 0.5$. Clusters will be virialised in $\approx 3.3$ Gyrs. Clusters in superclusters will not merge, and their present-day turnaround regions will be virialised in $\approx 10$ Gyrs. The large variety of properties of clusters suggests that they have followed different paths during evolution.

The pre-main-sequence evolution of the chemically peculiar (CP) stars on the upper main sequence is still a vast mystery and not well understood. Our analysis of young associations and open clusters aims to find (very) young CP stars to try to put a lower boundary on the age of such objects. Using three catalogues of open clusters and associations, we determined membership probabilities using HDBSCAN. The hot stars from this selection were submitted to synthetic $\Delta a$ photometry, spectral, and light curve classification to determine which ones are CP stars and candidates. Subsequently, we used spectral energy distribution fitting and emission line analysis to check for possible PMS CP stars. The results were compared to the literature. We detected 971 CP stars and candidates {in 217 clusters and associations}. A relatively large fraction, $\sim$10\% of those, show characteristics of PMS CP stars. This significantly expands the known list of candidate PMS CP stars, bringing us closer to solving the mystery of their origin.

Light scalar particles arise naturally in many extensions of the Standard Model and are well-motivated dark matter candidates. Gravitational interactions near black holes can trigger the growth of dense scalar configurations that, if sustained during inspiral, alter binary dynamics and imprint signatures on gravitational-wave signals. Detecting such effects would provide a novel probe of fundamental physics and dark matter. Here we develop a semi-analytic waveform model for binaries in scalar environments, validated against numerical relativity simulations, and apply it in a Bayesian analysis of the LIGO-Virgo-KAGRA catalog. Our results set physically meaningful upper bounds on scalar environments around compact binaries. When superradiance priors are included, we find tentative evidence for such an environment in GW190728 with $\ln B_{\mathrm{vac}}^{\mathrm{env}} \approx 3.5$, which would correspond to the existence of a light scalar field with mass $\sim 10^{-12}\,\mathrm{eV}$.

Axions that couple to nuclear spins via the axial current interaction can be both produced and detected using nuclear magnetic resonance (NMR) techniques. In this scheme, nuclei driven by a real oscillating magnetic field in one device act as an axion source, which can drive NMR in a nearby spin-polarized sample interrogated with a sensitive magnetometer. We study the prospects for detecting axions through this method and identify two key characteristics that result in compelling detection sensitivity. First, the gradient of the generated axion field can be substantial, set by the inverse distance from the source. In the near zone, it reduces to the inverse of the source's geometric size. Second, because the generated axion field is produced at a known frequency, the detection medium can be tuned precisely to this frequency, enabling long interrogation times. We show that the experimental sensitivity of a pair of centimeter-scale NMR devices operating over a 15-day integration time can already surpass existing astrophysical bounds on the axion-nucleon coupling. A similar sensitivity can be achieved with 10 centimeter-scale NMR devices with only 1 hour of integration time. These dual NMR configurations are capable of probing a wide range of axion masses, up to values comparable to the inverse distance between the source and the sensor.

Dark matter may not be perfectly stable, and its decay could generate distinctive gravitational-wave signatures. In this work, we present model-independent predictions for the stochastic gravitational-wave background arising from the decay of ultralight dark matter into gravitons. Within this framework, we forecast the sensitivity reach of current and forthcoming gravitational-wave detectors to such signals.

We present a complete pipeline for detecting and characterising gravitational waves (GWs) produced by the inspiral of stellar-mass binary black holes in data from the Laser Interferometer Space Antenna (LISA). The analysis framework relies on an efficient time-frequency implementation of an adaptive semi-coherent detection statistic, which we show to be robust against non-stationary noise and the presence of gaps of varying duration and cadence. The search is able to detect signals down to a coherent signal-to-noise ratio $\approx$ 11 over the full parameter space of black holes with spins aligned to the orbital angular momentum and orbital eccentricity $\leq$ 0.01 when deployed on the 2-year-long LISA Data Challenge Yorsh. The search can be run within a day using $\approx$ 40 GPUs. The techniques presented here have wider applications in GW astronomy, in particular the search for extreme-mass-ratio inspirals in LISA data.

We numerically study the formation of the gauge-mediation type Q balls in the logarithmic square potential on three-dimensional lattices. We obtain the broad charge distribution of the Q ball of this type for the first time. The charge of the Q ball at the peak of the distribution is smaller than what we estimated as the average of the largest tens of the Q balls in the logarithmic potential for the same initial amplitude of the field at the onset of its oscillation. We also discuss some impacts of the broad distribution on cosmology and astrophysics. In the B ball (Q being the baryon number) case, the broad distribution would lead to the coexistence of both stable and unstable B balls. We find that stable B balls can account for the dark matter of the universe without affecting successful big bang nucleosynthesis by the decay of the unstable B balls, but the baryon number of the universe cannot be explained by them. On the other hand, the large L balls (Q being the lepton number) would be the dark matter as well while avoiding the constraints on the X and/or gamma rays from the decay of the smaller L balls.

We provide insight about the full form of the equations for matter density perturbations and the scalar Bardeen metric potentials in general $f(R)$ theories of gravity. When considering viable modifications to the standard $\Lambda$CDM background, the full scale-dependent equations for the metric perturbations are provided and are shown to match the ones obtained with the quasistatic approximation. We investigate the impact of the $n=2$ Hu-Sawicki model on the late-time growth of structures. We find that updated late-time growth of structure data imposes $|f_{R_0}|\lesssim10^{-6}-10^{-5}$ and thus conclude that the Hu-Sawicki $f(R)$ model contributes no significant phenomenology at both background and perturbative level beyond the effective cosmological constant encompassed in its definition. This conclusion points to the survival of the present tension between early and late measurements of $\sigma_8$, as the Hu-Sawicki model can only worsen this issue or at best reproduce the results from the current concordance cosmological model. The generalized perturbative method showcased in this work can be applied to more elaborate $f(R)$ models to isolate genuine higher-order signatures beyond the quasistatic approximation.

Dynamical captures of black holes are unique events that provide an exceptional opportunity to probe the strong-field regime of gravitational physics. In this article, we perform numerical relativity simulations to study the events of dynamical capture of two equal-mass non-spinning black holes. We consider a suite of scenarios within a range of initial linear momenta ($p/M=0.095-0.75$) and incidence angles ($\theta=6.36^\circ-2.83^\circ$), and study the emitted Weyl scalar ($\Psi_4$) of each case, as well as the spins and masses of the black holes before and after they merge. We provide a simple analytical model which accurately fits the gravitational-wave emission. We study the dependence of the time-interval between the capture and the merger emissions with respect to the incidence angle, which can be well parametrized by a first-order divergent behavior, allowing to find the angle that separates a scattering event from a dynamical capture. We also find that, in general, the parameters that model the first emission can be well described by linear or exponentially decaying functions in terms of the incidence angle, while others display more complex behaviors that offer valuable insights into the nature of these events.