Papers for Tuesday, Jul 01 2025

We present a comprehensive exploration of hierarchical triple black hole (BH) systems to address the "initial separation" problem in gravitational wave astrophysics. This problem arises because isolated BH binaries must have extremely small initial separations to merge within a Hubble time via gravitational wave (GW) emission alone, separations at which their stellar progenitors would have merged prematurely. Using a modified JADE secular code incorporating GW energy loss, we systematically investigate a seven-dimensional parameter space: masses of three BHs (5-100 $M_\odot$ inner binary, 1-200 $M_\odot$ tertiary), inner/outer semimajor axes (1-200 AU and 100-10,000 AU), outer orbit eccentricity (0-0.9), and mutual inclination (40°-80°). We employed an adaptive MCMC approach sampling the merger/nonmerger transition boundary across nearly 15 million simulations. Results reveal merger-conducive regions correspond to asymmetric inner binary masses, large inner separations where von Zeipel-Lidov-Kozai (ZLK) mechanism operates effectively without relativistic precession suppression, small outer separations providing stronger perturbations, and large outer eccentricities bringing the tertiary closer at pericenter. Merger probability correlates positively with mutual inclination. We developed a classification scheme for nonmerging systems based on GW emission and ZLK oscillations. A trained neural network predicts merger outcomes with 99% ROC score and 95% accuracy (99.7% for high-confidence predictions), enabling rapid population synthesis. Validation against N-body integrations showed 87% qualitative agreement, confirming our methodology captures essential triple BH dynamics while enabling unprecedented-scale exploration of configurations resolving the initial separation problem.

The increasing number of detailed exoplanet observations offers an opportunity to refine our analyses and interpretations, and uncover previously-hidden trends. One such trend is unexpectedly high water abundances ($\sim$1$-$10\%) for several hot Jupiters. Here, we show that such atmospheres may instead be helium-rich, and the perceived high water abundances simply biases in analyses compensating for ``ad hocly'' fixed helium-to-hydrogen ratio (He/H$_2$ = 0.157). As transmission spectra constrain the atmospheric scale height ($H$), a helium-enriched atmosphere can be misinterpreted as H$_2$-dominated water-rich to bring the mean molecular weight ($\mu$) to intermediate values ($\sim$4$-$15) when He/H$_2$ is fixed. We present a proof-of-concept reanalysis of HD~209458~b's JWST transmission spectrum treating He/H$_2$ as a free parameter, resulting in sets of solutions with helium enrichment and non-extreme water abundance. The He-enhancement found requires an order of magnitude less enrichment/evolution than does a water abundance of 10\%, suggesting that the He-rich hypothesis is a viable alternative. We argue that He enhancement must be considered to reliably constrain atmospheric composition, be sensitive to a more diverse planetary population, and ultimately yield robust metallicity trends to inform formation/evolution pathways. Looking ahead, we suggest leveraging insights from differences in pressure-broadening effects, outflow measurements, and atmospheric chemistry to help disentangle reliably between He- and volatile-rich atmospheres -- while recognizing that existing opacity/outflow/chemistry models may need targeted upgrades to reach the required fidelity level.

JWST has revealed a population of "dormant" galaxies at z>5, which have recently halted their star formation and are characterized by weak emission lines and significant Balmer breaks. So far, only three such galaxies have been published at z>5 with low stellar masses, $M_*<10^9M_\odot$ (so-called mini-quenched galaxies), and one massive quiescent galaxy at z>5 with $M_*=10^{10.2}M_\odot$, without any such galaxy reported at intermediate masses. Here, we present a systematic search for dormant galaxies at 5<z<7.4 which halted star formation at least 10 Myr before the time of observation. To do this, we make use of all the publicly available NIRSpec prism data in the DAWN JWST Archive and select galaxies showing low H{\alpha} equivalent widths ($EW_{0}<50$Å) and strong Balmer breaks ($F_{\nu,4200}/F_{\nu,3500}>1.4$). We find 14 dormant galaxies with stellar masses ranging from $10^{7.6}-10^{10.5}$, revealing an intermediate-mass population of such galaxies. By construction, the 14 sources reported are located about 1 dex below the star-forming main sequence. Their star formation histories suggest that they halted star formation between 10 and 25 Myr before the time of observation which, according to models, is comparable with the timescales of internally regulated bursts driving a "breathing" mode of star formation. Our results show that ~1% of the galaxies in the DJA are in a dormant phase of their star formation histories and span a wide stellar mass range. These galaxies can be empirically selected using only their spectral features in NIRSpec prism data.

The LISA datastream will be populated by large instrumental and astrophysical noises, both potentially exhibiting long-term non-stationarities. Modelling and inferring on them is a challenging task, central for accurate signal reconstruction. In this paper, we introduce $\texttt{bahamas}$, a codebase designed to characterize noises and stochastic gravitational wave backgrounds (SGWBs) in LISA. $\texttt{bahamas}$ adopts a time-frequency data representation, based on the Short Time Fourier Transform, to accurately describe the signal temporal evolution and accommodate for the presence of data gaps. In addition, $\texttt{bahamas}$ supports a variety of SGWB spectral models proposed in literature, enabling joint inference on them. Posterior sampling leverages No-U-Turn sampling an efficient variant of Hamiltonian Monte Carlo, inheriting the cross-hardware capabilities provided by NumPyro (CPU/GPU/TPU). We benchmark $\texttt{bahamas}$ performances on a simple test case, and present ongoing developments to appear in future releases.

The Laser Interferometer Space Antenna (LISA) is expected to detect thousands of individually resolved gravitational wave sources, overlapping in time and frequency, on top of unresolved astrophysical and/or primordial backgrounds. Disentangling resolved sources from backgrounds and extracting their parameters in a computationally intensive "global fit" is normally regarded as a necessary step toward reconstructing the properties of the underlying astrophysical populations. Here, we show that it is possible to infer the properties of the most numerous population of LISA sources - Galactic double white dwarfs - directly from the frequency (or, equivalently, time) strain series, by using a simulation-based approach that bypasses the global fit entirely. By training a normalizing flow on a custom-designed compression of simulated LISA frequency series from the Galactic double white dwarf population, we demonstrate how to infer the posterior distribution of population parameters (e.g., mass function, frequency, and spatial distributions). This allows for extracting information on the population parameters from both resolved and unresolved sources simultaneously and in a computationally efficient manner. Our approach to target population properties directly can be readily extended to other source classes (e.g., massive and stellar-mass black holes, extreme mass ratio inspirals), provided fast simulations are available, and to scenarios involving non-Gaussian or non-stationary noise (e.g., data gaps).

We present spatially-resolved rest-frame optical emission line maps of four galaxies at $z \sim 2$ observed with Keck/OSIRIS to study the physical conditions of the ISM at Cosmic Noon. Our analysis of strong emission line ratios in these galaxies reveals an offset from the local star-forming locus on the BPT diagram, but agrees with other star-forming galaxies at similar redshifts. Despite the offset towards higher [O III]$\lambda5008$/H$\beta$ and [N II]$\lambda6585$/H$\alpha$, these strong-line ratios remain consistent with or below the maximum starburst threshold even in the inner $\sim 1$ kpc region of the galaxies, providing no compelling evidence for central AGN activity. The galaxies also exhibit flat radial gas-phase metallicity gradients, consistent with previous studies of $z \sim 2$ galaxies and suggesting efficient radial mixing possibly driven by strong outflows from intense star formation. Overall, our results reveal the highly star-forming nature of these galaxies, with the potential to launch outflows that flatten metallicity gradients through significant radial gas mixing. Future observations with JWST/NIRSpec are crucial to detect fainter emission lines at higher spatial resolution to further constrain the physical processes and ionization mechanisms that shape the ISM during Cosmic Noon.

Microquasars, compact binary systems with an accreting stellar-mass black hole or neutron star, are promising candidates for high-energy particle acceleration. Recently, the LHAASO collaboration reported on the detection of $>100$ TeV $\gamma$-ray emission from five microquasars, suggesting that these sources are efficient particle accelerators. In microquasars, high-energy $\gamma$-rays can be produced in large-scale jets or winds. In this work, we explore the X-ray, $\gamma$-ray and neutrino emission from SS 433, V4641 Sgr and GRS 1905+105. We consider leptonic and hadronic scenarios to explain the spectra observed by LHAASO and other high-energy $\gamma$-ray detectors. We estimate the neutrino flux associated with the hadronic component and investigate the detectability of neutrinos from these sources in current and future neutrino telescopes. We find that among the three sources, V4641 Sgr has the best prospects of observation with a combined next-generation neutrino telescopes.

The distance duality relation (DDR) between luminosity and angular diameter distances holds if gravity is described by a metric theory and the universe is transparent. Recent cosmological inferences using Type Ia supernovae (SNe~Ia), baryon acoustic oscillation (BAO) and the cosmic microwave background (CMB) observations have suggested that dark energy may evolve in time. We test how the assumption of distance duality impacts dark energy inference. Marginalizing over the absolute SNe~Ia luminosity, we find no deviation from the DDR, independent of the SN~Ia compilation used, or the assumed dark energy model. This corresponds a maximum deviation in the SN~Ia luminosity of $\Delta m \sim 0.05$ mag at the highest redshift. Allowing for deviations in the DDR increases the errors in the dark energy equation of state parameters (EoS) by 30-50$\%$. For the Pantheon+ compilation the constraints on dark energy are within $2\sigma$ of $\Lambda$CDM when applying a more realistic minimum redshift cut $z_{\rm min} >0.023$. We constrain possible physical scenarios that can impact cosmic transparency, specifically photon-axion mixing and the presence of (gray) intergalactic (IG) dust, in the latter case limiting the dust density to $\Omega_{\rm dust}<2 \times 10^{-4}$ at 95\% C.L. When using the local Cepheid calibration of the SNe~Ia absolute luminosity, a significant deviation from the DDR relation (for which $\epsilon=0$) is preferred. However, the best fit parameter value $\epsilon = -0.091 \pm 0.024$ requires the SNe~Ia to be brighter than the case with DDR valid, therefore, neither of the physical scenarios which dim the SNe can explain the high $H_0$. Constraints on the photon-axion interaction length scale suggest a limit on the coupling constant of $g_{a\gamma} < 10^{-12}\,{\rm GeV}^{-1}$.

The Large Scale Polarization Explorer (LSPE) project, funded by the Italian Space Agency (ASI), includes the development of LSPE-Strip, a ground-based radio telescope for observing Cosmic Microwave Background (CMB) anisotropies. LSPE-Strip, nearing its construction phase, will operate from the Teide Observatory in Tenerife, employing 49 coherent polarimeters at 43 GHz to deliver critical data on CMB anisotropies and 6 channels at 95 GHz as atmospheric monitor. On-site characterization of such advanced instruments is crucial to detect possible systematic effects, such as gain fluctuations, beam distortions, and pointing errors, that can compromise performance by introducing spurious polarizations or radiation collection from unintended directions. To address these challenges, a drone-mounted Q-band test source for on-site characterization of LSPE-Strip's polarimeter array was developed. Modern Unmanned Aerial Vehicles (UAVs) offer a flexible approach for antenna pattern measurements, yet their use in high-frequency radio astronomy is not consolidated practice. In October 2022, a UAV-based measurement campaign was conducted with the TFGI instrument on the second QUIJOTE telescope in Tenerife, in collaboration with the Instituto de Astrofisica de Canarias. This pioneering effort aimed to validate UAV-based beam characterization methods and assess QUIJOTE's performance under operational conditions. Preliminary results demonstrated high measurement accuracy, leveraging QUIJOTE's dual-receiver configuration for beam validation. These findings provide valuable insights for optimizing UAV systems in preparation for LSPE-Strip's future characterization.

We predict the spectrum and composition of ultrahigh energy cosmic rays (UCRs for short) generated in a binary neutron star merger, exploiting the highly constrained initial conditions to make quantitative predictions for the cutoff energy for various nuclei. The key mechanism producing UCRs heavier than helium is acceleration in the magnetized turbulent outflow outside the jets. We find the rigidity cutoff of this component is $\mathcal{R}_{\rm cut} \equiv E_{\rm cut}/eZ \approx 6-9$ EV, consistent with the measured value $\mathcal{R}_{\rm cut} = 6.8^{+5.8}_{-2.8}\,$EV from fitting data. This agreement strengthens the case that BNS mergers are the main site of UCR production. The jets may also be capable of accelerating particles to ultrahigh energy: these UCRs would be protons and/or helium, in a ratio that depends on the neutron star equation of state. We estimate the cutoff energies of such a component to be $E_{\rm cut}\approx\! 10\,$EeV for p and $\approx\! 30\,$EeV for He. Such a jet component may explain the hint of a secondary light population at higher energy found in the analysis of Muzio et al. (2019). With more precise data, the locus of production of helium and intermediate mass nuclei can be inferred from their rigidity structure, and the predicted small deviations from a pure Peters Cycle can be checked. Predicting the relative abundances of different elements and the total energy in UCRs per merger event is feasible, but beyond the scope of this paper. This scenario predicts that each neutrino above 1 PEV is co-directional with a gravitational wave arriving ~1 day earlier, and that the highest energy UCRs have masses heavier than iron.

Gravitational wave (GW) neutron star mergers with an associated electromagnetic counterpart constitute powerful probes of binary evolution, the production sites of heavy elements, general relativity, and the expansion of the universe. Only a handful of candidate GW binary mergers during the fourth LIGO/Virgo/KAGRA observing run (O4) so far are believed to include a neutron star. We present optical-near infrared follow-up observations of the candidate neutron-star black hole GW merger S250206dm. This is the first high-significance mass gap neutron star-black hole candidate observed by multiple GW detectors (thus having a significantly smaller sky localization than one-detector events), offering the first opportunity to effectively follow up a GW event of this kind. Our GW MultiMessenger Astronomy DECam Survey (GW-MMADS) campaign consisted of a wide-field search using the Dark Energy Camera (DECam) and T80-South (T80S), as well as galaxy-targeted observations using the Southern Astrophysical Research (SOAR) imager and the Wendelstein 2.1m 3-channel camera. No viable kilonova counterpart was found in our observations. We use our observation depths to place competitive constraints on kilonova models similar to or brighter than the GW170817 kilonova AT 2017gfo within our observed fields, ruling out 100\% of such models with SOAR galaxy-targeted observations and $\sim43$\% (48\%) with DECam (DECam and T80S).

Classical Cepheids are key astrophysical laboratories for studying stellar structure and evolution, and for calibrating the cosmic distance scale. Despite major progress, uncertainties remain regarding their masses, luminosities, distances, and the role of processes like core overshooting, rotation, and mass loss. The high-precision data from ESA Gaia's third data release offer an opportunity to address these issues. This study aims to estimate the structural parameters and distances of a sample of Classical Cepheids using non-linear convective pulsation models. It also tests the consistency of Gaia parallaxes, constrains the Mass Luminosity relation, and investigates the dependence of the projection factor (p-factor) on pulsation period. We analysed 46 Cepheids with precise Gaia DR3 photometry and radial velocities. Model fitting was performed by comparing predicted and observed light curves in Gaia G, GBP, and GRP bands, and radial velocity curves. Distances from the models were compared with Gaia parallaxes, including published corrections. Predicted masses and luminosities were used to constrain the ML relation; radial velocities allowed independent pfactor estimates. Distances show statistical agreement with Gaia parallaxes, requiring no additional global offset. Inferred masses and luminosities are consistent with evolution models including mild core overshooting, mass loss, or rotation. No significant period dependence of the pfactor is found, with a constant value of p = 1.22+/-0.05, consistent with recent studies. Results also agree with the Period-Wesenheit-Metallicity relation from Gaia DR3.

Weak gravitational lensing is a powerful probe for constraining cosmological parameters, but its success relies on accurate shear measurements. In this paper, we use image simulations to investigate how a joint analysis of high-resolution space-based and deep ground-based imaging can improve shear estimation. We simulate two scenarios: a grid-based setup, where galaxies are placed on a regular grid to mimic an idealised, blending-free scenario, and a random setup, where galaxies are randomly distributed to capture the impact of blending. Comparing these cases, we find that blending introduces significant biases, particularly in LSST-like data due to its larger point spread function. This highlights the importance of including realistic blending effects when evaluating the performance of joint analyses. Using simulations that account for blending, we find that the most effective catalogue-level synergy is achieved by combining all galaxies detected in either survey. This approach yields an effective galaxy number density of $44.08~\rm arcmin^{-2}$ over the magnitude range of 20.0 to 27.5, compared to $39.17~\rm arcmin^{-2}$ for LSST-like data alone and $30.31~\rm arcmin^{-2}$ for \textit{Euclid}-like data alone. Restricting the analysis to only the overlapping sources detected in both surveys results in a gain of ${\sim}12\%$ compared to using either survey alone. While this joint-object approach is suboptimal at the catalogue level, it may become more effective in pixel-level analyses, where a joint fit to individual galaxy shapes can better leverage the complementary strengths of both data sets. Studying such pixel-level combinations, with realistic blending effects properly accounted for, remains a promising direction for future work.

Studying compact-object binary mergers in star clusters is crucial for understanding stellar evolution and dynamical interactions in galaxies. Open clusters in particular are more abundant over cosmic time than globular clusters, however, previous research on low-mass clusters with $\lesssim 10^3~\textrm{M}_{\odot}$ has focused on binary black holes (BBHs) or black hole-neutron star (BH-NS) binaries. Binary mergers of other compact objects, such as white dwarfs (WDs), are also crucial as progenitors of transient phenomena such as Type Ia supernovae and Fast Radio Bursts. We present simulations of three types of open clusters with masses of $10^2$, $10^3$, and $10^4~\mathrm{M}_{\odot}$. In massive clusters with $\gtrsim 10^4~\textrm{M}_{\odot}$, BBHs are dynamically formed, however, less massive compact binaries such as WD-WD and WD-NS are perturbed inside the star clusters, causing them to evolve into other objects. We further find BH-NS mergers only in $10^3~\textrm{M}_{\odot}$ clusters. Considering star clusters with a typical open cluster mass, we observe that WD-WD merger rates slightly increase for $10^3~\textrm{M}_{\odot}$ clusters but decrease for $10^2~\textrm{M}_{\odot}$ clusters. Since the host clusters are tidally disrupted, most of them merge outside of the clusters. Our WD-WD merger results have further implications for two classes of transients. Super-Chandrasekhar WD-WD mergers are present in our simulations, demonstrating potential sources of Fast Radio Bursts at a rate of 70-780 Gpc$^{-3}$yr$^{-1}$, higher than the rate estimated for globular clusters. Additionally, we find that Carbon-Oxygen WD-WD mergers in our open clusters (34-640 $\textrm{Gpc}^{-3}$yr$^{-1}$) only account for 0.14-2.6% of the observed Type Ia supernova rate in our local Universe.

The UTMOST-NS pulsar timing programme operated at the Molonglo Observatory Synthesis Telescope from April 2021 to June 2023, observing 173 pulsars with an average cadence of 50 pulsars per day. An overview of the programme is presented, detailing the hardware, software, and observing strategy. Pulsar timing results are discussed, focusing on timing noise and glitches. It is shown that the scaling of residuals due to timing noise with pulsar parameters and observing timespan is consistent with earlier studies, and that the recovered timing noise parameters remain consistent as the observing timespan is increased. Second frequency derivatives are investigated, and it is shown that the uncertainty on $\ddot{\nu}$ is sensitive to the frequency cutoff in the timing noise model, varying by three-fold approximately depending on whether Fourier modes with frequency lower than the reciprocal of the observing timespan are included. We measure 39 non-zero values of $\ddot{\nu}$ when considering both models with and without low-frequency modes. An analytic scaling relating anomalous braking indices to timing noise amplitude is also validated. Glitches in the sample are discussed, including three detected by an ``online'' glitch detection pipeline using a hidden Markov model (HMM). In total 17 glitches are discussed, one of which, in PSR J1902+0615, has not been reported elsewhere. An ``offline'' glitch search pipeline using the HMM framework is used to search for previously undetected glitches. Systematic upper limits are set on the size of undetected glitches. The mean upper limit is $\Delta\nu^{90\%}/\nu = 6.3 \times 10^{-9}$ at 90\% confidence.

We present a new, state-of-the-art computation of the Diffuse Supernova Neutrino Background (DSNB), where we use neutrino spectra from multi-dimensional, multi-second core collapse supernova simulations - including both neutron-star and black-hole forming collapses - and binary evolution effects from modern population synthesis codes. Large sets of numerical results are processed and connected in a consistent manner, using two key quantities: the mass of the star's Carbon-Oxygen (CO) core at an advanced pre-collapse stage - which depends on binary evolution effects - and the compactness parameter, which is the main descriptor of the post-collapse neutrino emission. The method enables us to model the neutrino emission of a very diverse, binary-affected population of stars, which cannot unambiguously be mapped in detail by existing core collapse simulations. We find that including black hole-forming collapses enhances the DSNB by up to 50% at energies greater than 30-40 MeV. Binary evolution effects can change the total rate of collapses and generate a sub-population of high core mass stars that are stronger neutrino emitters. However, the net effect on the DSNB is moderate - up to a 15% increase in flux - due to the rarity of these super-massive cores and to the relatively modest dependence of the neutrino emission on the CO core mass. The methodology presented here is suitable for extensions and generalizations, and therefore it lays the foundation for modern treatments of the DSNB.

We apply a halo velocity bias model, $\gamma_{f}$, within the Aemulus simulation suite for General Relativity (GR) to investigate its efficacy in identifying the signature of assembly bias and Modified Gravity (MG). In the investigation of assembly bias, utilizing galaxy clustering data ranging from scales of $0.1 - 60 \text{Mpc}\,h^{-1}$, we discover that our emulator model accurately recreates the cosmological parameters, $\Omega_m$ and $\sigma_8$, along with the velocity bias $\gamma_{f}$, staying well within the 1-$\sigma$ error margins, provided that assembly bias is considered. Ignoring assembly bias can considerably alter our model constraints on parameters $\Omega_m$ and $\sigma_8$ if the test sample includes assembly bias. Using our emulator for MG simulations, which encompasses two Dvali-Gabadadze-Porrati models (DGP; N1, N5) and two $f(R)$ models (F4, F6), we can effectively identify a robust signature of modified gravity, for models such as DGP-N1 and $f(R)$-F4, as indicated by a noticeable deviation of $\gamma_{f}$ from unity. Using the velocity dispersion of dark matter halos to effectively represent the underlying strength of the velocity field of these MG simulations, we find that the simple $\gamma_{f}$ model can recover the truth with minimal bias. These evaluations indicate that our simple halo-velocity bias model is capable of detecting significant MG characteristics, although additional methodologies should be pursued to improve model constraints.

Inspired by recent observations suggesting that the retrograde precession of the brown dwarf binary 2M1510 AB is consistent with induction by a polar circumbinary planet, we investigate the formation of such planets by studying the evolution of a primordial misaligned circumbinary disc around a brown dwarf binary. Analytical calculations show that a critical tilt angle of $i_{\rm crit} \gtrsim 50^\circ$ for moderately eccentric binaries is needed for polar alignment of circumbinary discs in systems with low disc-to-binary angular momentum ratios. For higher ratios, this angle converges to the Kozai-Lidov instability threshold of $\sim 39^\circ$. We identify disc parameters, such as viscosity ($\alpha = 10^{-4}$) and aspect ratio ($H/r = 0.05$), that enable polar alignment within typical disc lifetimes. Notably, a circumbinary disc around a low-mass binary, such as a brown dwarf binary, will require more time to achieve polar alignment compared to higher-mass systems. A hydrodynamical simulation confirms that an initially inclined disc around a brown dwarf evolves towards a polar state, creating favorable conditions for polar planet formation. Using these results, we finish by placing 2M1510 AB into a wider context and speculate why such a polar circumbinary configuration has not been identified before.

We present a statistically grounded methodology for recovering physically consistent initial conditions from historical two-line element sets (TLEs), enabling accurate long-arc trajectory reconstruction for distant and highly eccentric Earth satellites. The approach combines numerical averaging of osculating orbital elements with Gaussian-mixture-model (GMM) filtering with robust fallback strategies to isolate a dynamically coherent subset of mean elements at a common reference epoch. From this filtered ensemble, a representative osculating element is reconstructed, yielding a recovered Cartesian state vector with predictive capability far exceeding that of raw TLE-based propagations and existing approaches. We apply this method to the case of OGO-1 (1964-054A), a spacecraft launched into a cislunar orbit and tracked intermittently over five decades. Despite large observational gaps and significant secular evolution, our recovered state, used within an unscented Kalman filter (UKF), accurately reproduces the full trajectory, including its atmospheric reentry in August 2020, to within one day of the true decay time. This result demonstrates the viability of filtered TLE statistics as a proxy for precise orbit determination, particularly in dynamical regimes where resonances and long-period perturbations dominate. The techniques presented here provide a framework for trajectory reconstruction and prediction using only publicly available data and are broadly applicable to the study of high-altitude debris objects and legacy space missions whose original tracking and covariance data are unavailable.

The analysis of the meridional displacement velocity of individual solar pores and sunspots has been performed. In the period May 2010\,-\,March 2025 of observations in the continuum of the \textit{Solar Dynamics Observatory/Helioseismic and Magnetic Imager} (SDO/HMI), we identified more than $3.6\cdot 10^{5}$ sunspots and pores for analysis and tracked their displacement. The velocity of the meridional displacement of spots $v_{\rm m}$ depends on their magnetic polarity, latitude, and stage of development. For sunspots and pores of trailing polarity, the velocity of movement is on average directed toward the poles. For such spots, the dependence of the velocity on latitude can be represented by linear regressions for pores: $v^{\rm pr}_{\rm tr} \approx 2.0+0.62\cdot \theta^{\rm o}$ m s$^{-1}$; for sunspots: $v^{\rm sp}_{\rm tr}\approx 0.02+0.94\cdot \theta ^{\rm o}$ m s$^{-1}$. For sunspots and pores of leading polarity, the dependence is non-monotonic in nature on latitude. For pores: $v^{\rm pr}_{\rm ld}\approx 0.35-11.7\cdot {\rm sin}(\theta)+16.5\cdot {\rm sin}^{\rm 2} (\theta) +76.5\cdot {\rm sin}^3 (\theta)-32.7\cdot {\rm sin}^{4}(\theta )$ m s$^{-1}$; for sunspots: $v^{\rm sp}_{\rm ld}\approx -0.35-18.3\cdot {\rm sin}(\theta)+32.2\cdot {\rm sin}^{\rm 2}(\theta) +71.4\cdot {\rm sin}^{\rm 3} (\theta)-6.7\cdot {\rm sin}^{\rm 4}(\theta)$ m s$^{-1}$. The highest speed of meridional movement to the poles is observed for sunspots of trailing polarity during the phase of growth of the sunspot area. The velocity of the meridional movement depends on their area, reaching a maximum for an area of $S\approx 80 - 100$ $\mu$sh.

We conducted a timing analysis of PSR J0922+0638 (B0919+06) using data from the Nanshan 26 m radio telescope and the MeerKAT telescope, spanning from January 2001 to March 2023. During this 22-year period, we discovered a previously unreported small glitch (glitch 1) before the well-known large glitch (glitch 2), occurring at ${\rm MJD} \sim 53325(3)$, with a frequency jump amplitude of $\Delta \nu/\nu \sim 0.79(6) \times 10^{-9}$. We also identified ten slow glitch events, half of which were newly detected. These slow glitches occurred quasi-periodically, with an average interval of approximately 553(21) days, fractional frequency changes ranging from $\Delta \nu/\nu \sim 1.13(1) \times 10^{-9}$ to $4.08(5) \times 10^{-9}$, and a maximum fractional change in the first derivative of the frequency of $\Delta \dot{\nu}/\dot{\nu} \sim -4.6 \times 10^{-3}$. Additionally, our timing noise analysis reveals a change in the spectral index for noise power before and after glitch 2, with values of $-6.0$ and $-5.3$, respectively, likely due to this large glitch. Throughout the entire observation period, the first derivative of the spin frequency ($\dot{\nu}$) showed a periodic structure. The possible modulation period was estimated to be 537(24) days before the 700-day data gap at MJD 56716 and 600(58) days afterward. We discuss the periodic oscillations in pulsar rotation as a possible manifestation of spin-down noise and quasi-periodic slow glitches.

We constrain the quark matter equation of state (EOS) with color superconductivity (2SC, 2SC+s, CFL phases) and QCD corrections using Bayesian inference informed by multi-messenger observations from the NICER pulsars, HESS J1731-347 and GW170817. Our results reveal that the bag constant $B$ = 47$\text{--}$55 MeV/fm$^{3}$ are tightly constrained, the QCD correction strength $a_4$ is phase-dependent ($a_4 \approx$ 0.9 for 2SC, $\sim$0.6 for 2SC+s). A large surface tension is required for 2SC phase ($\sigma \gtrsim 95$ MeV/fm$^2$) according to the stability constraints of ${^{266}\mathrm{Hs}}$. The mass-radius relations support quark stars as explanations for HESS J1731-347. The squared speed of sound converges to $c^2/3$, ruling out extreme stiffening. The inferred maximum quark star mass $M_{\mathrm{max}}< 2.28M_\odot$ (90\% confidence level) rules out the secondary component of GW190814 (2.59$\pm$0.08$M_\odot$) as a quark star. In particular, the quark stars (2SC phase) is strongly favored when including HESS J1731-347 data with the bayes factor $B_{\text{2SC, NS}}$ = 10.6, but neutron stars are moderately preferred with $B_{\text{NS, 2SC}}$ = 3.22 if HESS J1731-347 data is excluded. This demonstrates low-mass compact objects critically discriminate between quark matter phases and between quark and neutron-star models.

Understanding the origin of galactic angular momentum and its connection to the cosmic web remains a pivotal issue in galaxy formation. Using kinematic data from the MaNGA survey, we investigate the alignment between the spin directions of spiral galaxies and their host cosmic filaments. By incorporating filament spin measurements derived from redshift asymmetry across filament spines, we reveal a mass-dependent anti-parallel correlation: low-mass spiral galaxies ($\log_{10}(M_*/M_\odot) \lesssim 10$) exhibit a statistically significant anti-parallel alignment between their stellar/gas spins and filament spins, while high-mass spirals show no such trend. Spatial analysis further indicates that high-mass spirals preferentially reside near filament spines, whereas low-mass spirals occupy filament outskirts. These findings extend previous alignment studies that neglected directional spin correlations and provide new insights into how cosmic environments shape galactic angular momentum. The observed anti-parallel trend suggests a critical role for filament spin in regulating the angular momentum acquisition of low-mass spirals. This anti-parallel alignment is significantly enhanced for low-mass spirals residing in dynamically cold filaments, highlighting the importance of filament properties in shaping galaxy spin.

Low multipoles ($l$) in cosmic microwave background (CMB) temperature anisotropies have shown some `peculiarities' when examined since the release of the full sky CMB maps, using a variety of tests. In this paper, we concern ourselves with the very first peculiarities seen in CMB data viz., a breakdown of statistical isotropy in the form of axiality and planarity of these low-$l$ modes, and preferred alignments among them. We scrutinize the latest CMB data from ESA's Planck mission, PR4, to evaluate the current status of these deviations. We employ the Power tensor method which allows an invariant characterization of the distribution of power in a given multipole, and apply it to probe the first sixty multipoles i.e., $l=2$ to 61. We find that there are significant number of modes that are intrinsically anisotropic with a cumulative probability of $0.3\%$. However since the planarity study reveals that those modes that are unusually planar are subset of these anisotropic modes, we conclude that they may not be intrinsically planar. The quadrupole is still well aligned with the octopole. Besides, $l=3$, higher multipoles aligned with quadrupole are found to be insignificant. Interestingly, the collective alignment axis of the first sixty multipoles is found to be broadly closer to the axis of dipole, quadrupole, octopole and other modes aligned with $l=2$.

Maps of cosmic microwave background (CMB) are extracted from multi-frequency observations using a variety of cleaning procedures. However, in regions of strong microwave emission, particularly in the galactic plane from our own galaxy Milky Way and some extended or point sources, the recovered CMB signal is not reliable. Thus, a galactic mask is provided along with the \emph{cleaned} CMB sky for use with that CMB map which excises sky regions that may still be potentially contaminated even after cleaning. So, to avoid bias in our inferences, we impose such a foreground mask. In this paper, we analyze a cleaned CMB map from Planck PR4 to probe for any foreground residuals that may still be present \emph{outside} the galactic mask where the derived CMB sky is considered clean. To that end, we employ a local cross-correlation coefficient statistic where we cross-correlate widely used foreground templates that trace galactic synchrotron, free-free, and thermal dust emission from our galaxy with the cleaned CMB sky. Using simulations, we find that few regions of the derived CMB sky are still contaminated and have to be omitted. Based on this study, we derived a mask that could be used in conjugation with the standard mask to further improve the purpose of galactic masks.

We present a study aimed at quantifying the detectability of radio counterparts of binary neutron star (BNS) mergers with total masses $\lesssim 3$\,M$_{\odot}$, which may form neutron star (NS) remnants. We focus on mergers localized by gravitational-wave (GW) observations to sky areas $\lesssim 10$\,deg$^2$, a precision that greatly facilitates optical counterpart identification and enables radio discovery even without detections at other wavelengths. Widely separated GW detectors are essential for building samples of well-localized BNS mergers accessible to US-based radio telescopes, with minimum yearly detection rates (assuming the smallest values of the BNS local merger rate) ranging from a few with current GW detectors to hundreds with next-generation GW instruments. Current GW networks limit well-localized detections to $z\lesssim 0.2$, while next-generation GW detectors extend the reach to $z\lesssim 0.8$, encompassing the median redshift of short gamma-ray bursts (GRBs). With next-generation radio arrays operating at a several tens of GHz and providing an order of magnitude improvement in sensitivity compared to the most sensitive ones available today, short GRB-like jet afterglows can be detected for a large fraction of the considered BNS mergers. At lower radio frequencies, detections with current radio interferometric arrays are feasible, though subject to synchrotron self-absorption effects. The enhanced sensitivity and survey speed of future radio interferometers operating at a few GHz, combined with the higher detection rate of well-localized BNSs enabled by next-generation GW observatories, are key to probing disk-wind and dynamical ejecta afterglows, as well as remnant diversity.

The primary scientific objective of this Habitable Worlds Observatory (HWO) Science Case Development Document (SCDD) is to measure planetary rotation rates of transiting exoplanets to determine the structure, composition, circulation, and aerosol properties of their planetary atmospheres. For this analysis, HWO would obtain spectroscopic phase curves for planets with orbital periods of 5 - 20+ days, to assess tidal locking radius assumptions. Extending phase curve studies out to longer orbital periods than accessible with current and near-future telescopes will enable detailed investigation of atmospheric structure, composition, and circulation for planets that are much cooler than the more highly irradiated planets accessible with JWST phase curve observations (i.e., Teq < 500 K for HWO versus 1400 K <= Teq <= 2600 K for JWST). Broad wavelength coverage extending from the UV to the NIR would capture both reflected light and thermal emission, enabling HWO to conduct comprehensive characterization of planetary atmospheres. UV observations would probe high altitudes, thereby providing valuable insights into atmospheric (dis)equilibrium, aerosol properties, and the effects of photochemical processes on atmospheric composition. We also discuss the role of polarimetry in the classification of aerosols and the associated role they play in the atmospheric energy budget that directly ties them to the chemistry and circulation structure of the atmosphere.

Based on a unified description of various accretion flows, we find a long-ignored solution - the effectively optically thin accretion flow, occurring at accretion rates around Eddington value. As a consequence of radiation-pressure dominance, the density in a standard thin disc (SSD) decreases with the increase of accretion rates, making the innermost region effectively optically thin. Further increase in accretion rate leads to a rise of the temperature so that the Compton cooling is able to balance the accretion released energy. We demonstrate that the effectively optically thin flow is characterized by moderate temperature and large scattering optical depth, producing a multi-color Wien spectrum. For an appropriate accretion rate, the accretion flow transforms from an outer SSD to an inner effectively optically thin flow. Thus, the spectra of the whole accretion flow exhibit two components, i.e., a multi-color Wien spectrum at higher frequency and a multi-color blackbody, the former could provide an alternative origin of soft X-ray excess or formation of warm corona in active galactic nuclei (AGNs). Our stability analysis proves it is thermally stable and viscously unstable, indicating its existence in accreting systems. We show that effectively optically thin accretion flow exists in supermassive black holes for accretion rates around 0.1 to 10 times Eddington value, bridging the SSD at low accretion rates and slim disc at high rates. By comparing the predictions and average spectra of AGN, we constrain the viscosity parameter to be $\alpha \sim 0.03$, in good agreement with that derived from observed variability.

The association between long gamma-ray bursts (LGRBs) and core-collapse supernovae (SNe) has been well established since the discovery of SN 1998bw, which was linked to the low-luminosity LGRB 980425. However, long-term monitoring of several well-localized, low-redshift LGRBs has yielded compelling evidence for the absence of accompanying supernovae. Notably, two long bursts, GRB 211211A and GRB 230307A, show signatures consistent with kilonova emission from compact binary mergers, indicating that at least some long events may originate from progenitors other than core-collapse supernovae. In this study, we conduct a comparative analysis of two samples of LGRBs, i.e., LGRBs with and without supernova associations, to investigate the differences that may reveal intrinsic distinctions in their progenitors. A detailed examination of their prompt emission properties, host galaxy environments, and event rates is performed. While the two samples exhibit considerable overlap in most observed properties, a significant discrepancy in their event rate is revealed. LGRBs without supernova association have an event rate that aligns well with the star formation rate, whereas that of SN-associated LGRBs differs significantly. It indicates that LGRBs without a supernova association may constitute a distinct subclass with intrinsically different progenitors.

PSR J0537$-$6910 is the most frequently glitching pulsar, exhibiting large glitches at a rate of roughly thrice per year. It displays a negative long-term effective braking index $n' = -1.234(9)$ and roughly constant positive frequency second derivatives extending all the way from one glitch to the next. We use the published RXTE and NICER timing data (i) to explore whether `persistent shifts' (non-relaxing parts of the glitch $\Delta\dot{\nu}$ in the spin-down rate, like those observed in the Crab pulsar) can explain the negative effective braking index; (ii) to ascribe part of the inter-glitch relaxation at constant frequency second derivative $\ddot{\nu}$ to internal torques due to the coupling between the observed surface and the interior of the neutron star like deduced for the Vela pulsar; (iii) to demonstrate, as a proof of concept, that after taking account of the persistent shifts and internal torques, the true braking index $n$ associated with the pulsar braking torque is $n\sim 3 $, and (iv) to explore various empirical or model dependent correlations between glitch parameters and the time to the next glitch. The average persistent shift derived for PSR J0537$-$6910 is in agreement with the persistent shift values for the Crab pulsar. By applying the crustquake model, we infer a broken plate size of $D=80$ m, similar to the value obtained for the minimum glitch size observed in the Crab pulsar. The inter-glitch $\ddot{\nu}$ values inferred for the internal torques are commensurate with those in the Vela pulsar. Delineating the effects of persistent shifts and inter-glitch relaxation due to internal torques as observed in the Crab, Vela and other pulsars, we obtained a true braking index value $n=2.75(47)$ for PSR J0537$-$6910, similar to braking indices observed from other young pulsars. PSR J0537$-$6910 is not a likely source of gravitational waves.

In this paper, we employ two model-independent approaches, including redshift binning method and polynomial interpolation method, to reconstruct dark energy (DE) equation of state (EoS) $w(z)$ and DE density function $f(z)$. Our analysis incorporates data from the Dark Energy Spectroscopic Instrument (DESI) data release 2, Cosmic Microwave Background (CMB) distance priors from Planck 2018 and Atacama Cosmology Telescope data release 6, and three Type Ia supernovae (SN) compilations (PantheonPlus, Union3, and DESY5). To ensure model independence, we adopt three redshift binning schemes (n=3, 4, 5) and three polynomial interpolation schemes with the same number of nodes (n=3, 4, 5). Our main conclusions are as follows: 1) After taking into account DESI data, there is a clear trend that DE should evolve with redshift (with deviations from the cosmological constant reaching at least a $2.13\sigma$ confidence level), indicating that current observations favor a dynamical DE. 2) In the redshift range $0.5 < z < 1.5$, the DE EoS w(z) exhibits a decreasing trend and crosses the phantom divide $w=-1$, suggesting quintom-like behavior. 3) The DE density f(z) first increases at low redshift, reaching a hump around $z\approx 0.5$, and then decreases at $0.5 < z < 1.5$, with a rapid decrease at $z>1.0$. 4) For $z > 1.5$, current data are insufficient to place strong constraints on the evolution of DE, resulting in large uncertainties in the DE reconstruction. It must be emphasized that, these four main conclusions are independent of specific reconstruction models, and are insensitive to the choice of SN compilations.

The transient Galactic source MAXI J1744-294 went into an outburst in 2025 for the very first time. We study the spectral properties of this source during this outburst using archival data from the XRISM satellite for both of its Resolve and Xtend instruments. We have analyzed the source during one epoch, on March 03, 2025, or MJD 60737, on which XRISM data were available. Using both phenomenological and physical model fitting approaches for continuum emissions, along with line emission and interstellar absorption models, we analyzed the spectral data in the broad 2-10 keV energy band. From our spectral analysis, we have found the existence of multiple iron lines, which are different components of the Fe XXV emission. These line complexes arise from two highly ionized plasmas with ionization rate ~ 10000 erg cm/s with distinct turbulent velocities: one broad (vturb ~ 2513 km/s) from hot gas at the inner accretion disk, and one narrow (vturb ~ 153 km/s ) scattered by nearby photoionized gas. The source is a moderately spinning black hole with a spin of 0.63-0.70, a mass of 5.7-10.1 Solar masses, and a disk inclination angle of 19-24 degrees. The spectral model fitted parameters suggest that the source is in the soft spectral state. The source is situated in a crowded field near the Galactic center, resulting in a very high hydrogen column density.

The recent measurement of baryonic acoustic oscillations (BAO) by the Dark Energy Spectroscopic Instrument (DESI) reveals a mild tension with observations of the cosmic microwave background (CMB) within the standard $\Lambda$CDM cosmological model. This discrepancy leads to a preference for a total neutrino mass that is lower than the minimum value inferred from neutrino oscillation experiments. Alternatively, this tension can be eased within $\Lambda$CDM by assuming a higher optical depth ($\tau \simeq 0.09$), but such a value conflicts with large-scale CMB polarization data. We point out that cosmic birefringence, as suggested by recent Planck reanalyses, resolves this discrepancy if the birefringence angle varies significantly during reionization. Specifically, we consider the fact that the measured cosmic birefringence angle $\beta_0=0.34\pm0.09\,(1\,\sigma)\,$deg has the phase ambiguity, i.e., the measured rotation angle is described by $\beta=\beta_0+180n\,$deg ($n\in \mathbb{Z}$). We show that cosmic birefringence induced by axion-like particles with nonzero $n$ suppresses the reionization bump, allowing a higher $\tau$ consistent with data. We provide a viable parameter space where the birefringence effect simultaneously accounts for the low-$\ell$ polarization spectra, the Planck $EB$ correlations, and the elevated value of $\tau$, suggesting a key role for cosmic birefringence in current cosmological tensions.

During the Cosmic Dawn (CD), the HI 21-cm optical depth ($\tau$ ) in the intergalactic medium can become significantly large. Consequently, the second and higher-order terms of $\tau$ appearing in the Taylor expansion of the HI 21-cm differential brightness temperature ($\delta T_{\rm b}$ ) become important. This introduces additional non-Gaussianity into the signal. We study the impact of large $\tau$ on statistical quantities of HI 21-cm signal using a suite of standard numerical simulations that vary X-ray heating efficiency and the minimum halo mass required to host radiation sources. We find that the higher order terms suppress statistical quantities such as skewness, power-spectrum and bispectrum. However, the effect is found to be particularly strong on the non-Gaussian signal. We find that the change in skewness can reach several hundred percent in low X-ray heating scenarios. For moderate and high X-ray heating models changes are around $\sim40\%$ and $60\%$, respectively, for $M_{\rm h,min}=10^{9}\, {\rm M}_{\odot}$. This change is around $\sim 75\%$, $25\%$ and $20\%$ for low, moderate and high X-ray heating models, respectively, for $M_{\rm h,min}=10^{10}\, {\rm M}_{\odot}$. The change in bispectrum in both the halo cutoff mass scenarios ranges from $\sim 10\%$ to $\sim 300\%$ for low X-ray heating model. However, for moderate and high X-ray heating models the change remains between $\sim 10\%$ to $\sim 200\%$ for both equilateral and squeezed limit triangle configuration. Finally, we find that up to third orders of $\tau$ need to be retained to accurately model $\delta T_{\rm b}$, especially for capturing the non-Gaussian features in the HI 21-cm signal.

Dark matter is fundamental to the composition, structure, and evolution of the universe. Combining the ACT's cosmic microwave background, DESI's baryon acoustic oscillations with DESY5 type Ia supernova observations, we find a $3.4\,\sigma$ evidence for dynamical dark matter (DDM) with an equation of state, $\omega_{dm}(a)=\omega_{dm0}+\omega_{dma}(1-a)$. Independent of the Planck measurements, the ACT data confirms the linear relation $\omega_{dma}=-\omega_{dm0}$, inducing that the equation of state of dark matter is directly proportional to the scale factor $a$. Furthermore, the effects of DDM on the large-scale structure observables are thoroughly studied. Our findings are of great significance for understanding cosmic acceleration, structure growth, and the fate of the universe.

The combination of the high-resolution ALMA, JWST and HST observations provides unprecedented insights into the connection between individual molecular clouds and their internal stellar populations in nearby galaxies. The molecular clouds in five nearby galaxies were identified based on the integrated intensity maps of CO (2$-$1) emission from ALMA observations. We used the JWST 21 $\mu$m data to estimate the star formation rate (SFR) surface density of the clouds and calculate the masses of the embedded stellar populations in the clouds. After matching the star cluster and stellar association catalogs derived from the HST observations with the identified molecular clouds, we found clear correlations between the physical parameters of molecular clouds and their internal stellar populations. Based on the masses of the total stellar populations and their corresponding clouds, we obtained a typical value of the cloud-scale star formation efficiency (SFE), $\approx$1.4\%. The mass of the most massive cluster ($M_{\rm cluster, max}$) in a cloud is positively proportional to the mass ($M_{\rm cloud}$), the column density, the SFR sand the SFR surface density of the cloud. The observed $M_{\rm cluster, max}$-$M_{\rm cloud}$ relation can be interpreted theoretically on the basis of the integrated cloud-wide IMF (ICIMF) theory, which provides a quantitative framework for understanding the correlations between molecular clouds and their internal stellar populations.

The Epoch of Reionization marks the last major phase transition in the early Universe, during which the majority of neutral hydrogen once filling the intergalactic medium was ionized by the first galaxies. The James Webb Space Telescope is now identifying promising galaxy candidates capable of producing sufficient ionizing photons to drive this transformation. However, the fraction of these photons that escape into intergalactic space--the escape fraction--remains highly uncertain. Stellar feedback is thought to play a critical role in carving low-density channels that allow ionizing radiation to escape, but the dominant mechanisms, their operation, and their connection to observable signatures are not well understood. Local analogs of high-redshift galaxies offer a powerful alternative for studying these processes, since ionizing radiation is unobservable at high redshift due to intergalactic absorption. However, current UV space-based instrumentation lacks the spatial resolution and sensitivity required to fully address this problem. The core challenge lies in the multiscale nature of LyC escape: ionizing photons are generated on scales of 1-100 pc in super star clusters but must traverse the circumgalactic medium which can extend beyond 100 kpc. A UV integral field unit (IFU) spectrograph capable of resolving galaxies across these scales is necessary--and uniquely achievable with the proposed Habitable Worlds Observatory. In this article, we outline the scientific motivation, observables, and observational capabilities needed to make progress on these fundamental questions.

Future near-infrared spectroscopic galaxy surveys will target high-redshift emission-line galaxies (ELGs) to test cosmological models. Deriving optimal constraints from emission-line galaxy clustering hinges on a robust understanding of their environmental dependence. Using the TNG300-1 simulation, we explore the correlation between properties of H$\alpha$ emitters and their environment anisotropy rather than traditional density-based measures. Our galactic H$\alpha$ emission model includes contributions from the warm interstellar medium. The environment anisotropy and type are assigned using a halo mass-dependent smoothing scale. We find that most luminous ELGs ($L_{\rm{H}\alpha}>10^{42}\ \rm{erg\ s^{-1}}$) reside in filaments and knots. More generally, ELGs are more biased in strongly anisotropic environments. While correlations with galactic properties are found to be weak, they are statistically significant for host halo masses $M\lesssim 10^{12}\ M_\odot/h$. Our analysis motivates further investigation into how environmental anisotropy influences galaxy evolution, and highlights the potential for leveraging these effects in the analyses of upcoming cosmological surveys.

The origin of extragalactic, almost dark HI clouds with extreme gas-to-stellar mass ratios remains poorly understood. We investigate the nature and fate of the "almost dark" cloud AGC 226178, projected within the Virgo cluster, with an HI-to-stellar mass ratio of ~1000. We present deep single-dish HI mapping from the Five-hundred-meter Aperture Spherical Telescope (FAST), complemented by high-resolution interferometric data from the Very Large Array (VLA), as part of the Atomic gas in Virgo Interacting Dwarf galaxies (AVID) project. These observations provide the highest-quality HI analysis to date of such a cloud, combining resolution and sensitivity. FAST data reveal a short, low-velocity tail toward the dwarf galaxy VCC 2034, previously proposed as a possible origin for AGC 226178. However, VCC 2034 shows a line-of-sight asymmetric HI feature and cometary morphology indicating a stripping event unrelated to AGC 226178. VLA data reveal a velocity gradient across AGC 226178 and a clumpy internal structure. The velocity dispersion exceeds the thermal linewidth, implying turbulence or unresolved motions. The cloud cannot be gravitationally bound by atomic gas alone. The resolved HI clumps follow standard HI mass-star formation rate and mass-size relations, with those forming stars reaching surface densities above the threshold for self-shielding. We conclude that AGC 226178 is a free-floating HI cloud of unknown origin. The system appears to be in the process of disintegration. It is likely located well outside the Virgo cluster, as the preservation of its extended HI morphology within the cluster environment would otherwise require a substantial reservoir of unseen molecular gas with a mass exceeding that of the observed HI content. While confinement pressure from the hot intracluster medium may aid its stability, it is unlikely to be the dominant factor preventing its disruption.

The circumstellar environments of PNe are wonderful chemically rich astrophysical laboratories in which the ionization of atoms and the formation of simple and complex molecules can be studied. The new high-sensitivity receivers open the possibility to carry out deep observations, essential to unveil weak atomic and molecular spectra in the mm range. The main goal of this work is to study the emission lines detected in the spectra of the bright C-rich PNe IC 418 and NGC 7027 and to identify all those emission features associated with radio recombination lines (RRLs) of light elements. We aim to analyze the RRLs detected on each source, and to model the sources and derive their physical parameters. This work allows us to provide the most complete and updated catalog of RRLs in space, carried out at 2, 3 and 7 mm with the IRAM 30m and the Yebes 40m radio telescopes. We compare these observational data sets with synthetic models produced with the radiation transfer code Co3RaL. Our observations reveal the presence of several H and He I RRLs at mm wavelengths in the spectra of IC 418 and NGC 7027 and also of HeIIi RRLs in the spectrum of NGC 7027. Many of these lines had remained undetected until now due to their weakness and the lack of high-sensitivity observations at these frequencies. The data also confirm the absence of molecular emission towards IC 418, above a detection level of ~3 mK [Tmb]. These mm observations represent the most extended RRL line survey of two C-rich PNe carried out so far, with most of the lines never reported before. These extremely complete catalogs evidence the importance of high-sensitivity observations and are expected to be very helpful in the line identification process in mm observations, in particular for still unknown or poorly characterized molecular species existing in the vicinity of ionized environments.

In this paper we identify and explore a previously unidentified, multi meter-per-second, systematic correlation between time and RVs inferred through TM and LBL methods. We evaluate the influence of the data-driven stellar template in the RV bias and hypothesize on the possible sources of this effect. We first use the s-BART pipeline to extract RVs from three different datasets gathered over four nights of ESPRESSO and HARPS observations. Then, we demonstrate that the effect can be recovered on a larger sample of 19 targets, totaling 4124 ESPRESSO observations spread throughout 38 nights. We also showcase the presence of the bias in RVs extracted with the SERVAL and ARVE pipelines. Lastly, we explore the construction of the stellar template through the 5 years of ESPRESSO observations of HD10700, totalling more than 2000 observations. We find that a systematic quasi-linear bias affects the RV extraction with slopes that vary from -0.3 m/s-1/h-1 to -52 m/s-1/h-1 in our sample. This trend is not observed in CCF RVs and appears when all observations of a given star are collected within a short time-period (timescales of hours). We show that this systematic contamination exists in the RV time-series of two different template-matching pipelines, one line-by-line pipeline, and that it is agnostic to the spectrograph. We also find that this effect is connected with the construction of the stellar template, as we are able to mitigate it through a careful selection of the observations used to construct it. Our results suggest that a contamination of micro-telluric features, coupled other sources of correlated noise, could be the driving factor of this effect. We also show that this effect does not impact the usual usage of template-matching for the detection and characterization of exoplanets. Other short-timescale science cases can however be severely affected.

Mapping the spatial distributions and abundances of complex organic molecules in hot cores and hot corinos is crucial for understanding the astrochemical pathways and the inheritance of prebiotic material by nascent planetary systems. However, the line-rich spectra from these sources pose significant challenges for robustly fitting molecular parameters due to severe line blending and unidentified lines. We present an efficient framework, Bayesian Active Spectral-cube Inference and Learning (BASIL), for estimating molecular parameter maps for hundreds of molecules based on the local thermodynamic equilibrium (LTE) model, applied to wideband spectral datacubes of line-rich sources. We adopted stochastic variational inference to infer molecular parameters from spectra at individual positions, balancing between fitting accuracy and computational speed. For obtaining parameter maps, instead of querying every location or pixel, we introduced an active learning framework based on Bayesian quadrature and its parallelization. Specifically, we assessed and selected the locations or pixels of spectrum that are most informative for estimating the entire set of parameter maps by training a Gaussian processes model. By greedily selecting locations with maximum information gain, we achieve sublinear convergence. We benchmarked BASIL on a large synthetic datacube and demonstrated that it produces accurate 468 molecular parameter maps from 117 molecules within ~180 hours, orders of magnitude faster than traditional pixel-by-pixel fitting using Markov chain Monte Carlo methods, with visually reliable results emerging in just ~20 hours. Additional training iterations provide progressively more accurate results. This quick visualization meets the demands of big data in modern astronomical surveys.

We evaluated the jet power and the density of ambient matter in 3C 84 by using the momentum balance along the jet axis and the transonic condition for the cocoons observed at two different scales (approximately 1 and 6 parsec scales). For the inner cocoon, we precisely determined the ratio of jet power to ambient density $L_{\rm j}/n_{\rm a}$ to be $(0.3-0.7)\times 10^{43}\,{\rm erg}\,{\rm s}^{-1}\,{\rm cm}^3$. Similarly, for the outer cocoon, we found that this value is more than an order of magnitude larger at $(0.9-3.7)\times 10^{44}\,{\rm erg}\,{\rm s}^{-1}\,{\rm cm}^3$. This indicates that the outer cocoon is formed by a powerful jet that propagates through an ambient density of $20-300\,{\rm cm}^{-3}$ with a jet power of $10^{45-46.5}\,{\rm erg}\,{\rm s}^{-1}$. On the other hand, the inner cocoon is formed by a weaker jet with a power of $10^{43-44}\,{\rm erg}\,{\rm s}^{-1}$, propagating through a relatively low-density environment of $6-20\, {\rm cm}^{-3}$. These results suggest that: 1) with respect to the difference in $n_{\rm a}$, it appears to support the hypothesis that the inner cocoon, recently formed about 10 years ago, is expanding in the low-density cocoon created by the jet emitted about 25-50 years ago. 2) to achieve the short-lived and high $L_{\rm j}$ that generated the outer cocoon, a large mass accretion rate must be required over a short period to activate the jet. These may imply the extreme accretion event driven by the tidal disruption events (TDEs) of massive stars and/or the disk instability.

Gravitational systems in astrophysics often comprise a body -- the primary -- that far outweights the others, and which is taken as the centre of the reference frame. A fictitious acceleration, also known as the indirect term, must therefore be added to all other bodies in the system to compensate for the absence of motion of the primary. In this paper, we first stress that there is not one indirect term but as many indirect terms as there are bodies in the system that exert a gravitational pull on the primary. For instance, in the case of a protoplanetary disc with two planets, there are three indirect terms: one arising from the whole disc, and one per planet. We also highlight that the direct and indirect gravitational accelerations should be treated in a balanced way: the indirect term from one body should be applied to the other bodies in the system that feel its direct gravitational acceleration, and only to them. We point to situations where one of those terms is usually neglected however, which may lead to spurious results. These ideas are developed here for star-disc-planets interactions, for which we propose a recipe for the force to be applied onto a migrating planet, but they can easily be generalized to other astrophysical systems.

The recent detection of the ultra-high-energy neutrino event KM3-230213A ($\sim$220 PeV) by KM3NeT telescope poses a challenge to conventional astrophysical models, particularly in light of the absence of similar $\gtrsim$100 PeV events in IceCube data, despite its larger exposure. We propose a novel mechanism in which binary black hole mergers act as transient neutrino sources via gravitationally induced electroweak vacuum instability. In this scenario, the extreme spacetime curvature near the horizons during the final inspiral phase destabilizes the Higgs vacuum, triggering nucleation of true-vacuum bubbles. Collisions between these bubbles produce microscopic black holes that rapidly evaporate via Hawking radiation, emitting intense, short-lived bursts of neutrinos with energies exceeding 100 PeV. The resulting neutrino fluence follows a heavy-tailed distribution, allowing rare but highly luminous sources to account for events like KM3-230213A while remaining consistent with IceCube's non-detections. This framework links gravitational wave sources to ultra-high-energy neutrino production and suggests that future multi-messenger observations may detect electromagnetic signatures from microscopic black hole evaporation.

Impulsive radio signals such as fast radio bursts (FRBs) are imprinted with the signatures of multi-path propagation through ionised media in the form of frequency-dependent temporal broadening of the pulse profile (scattering). The dominant source of scattering for most FRBs is expected to be within their host galaxies, an assumption which can be tested by examining potential correlations between properties of the FRBs and global properties of their hosts. Using results from the CRAFT survey, we investigate correlations across a range of host galaxy properties against attributes of the FRB that encode propagation effects: scattering timescale tau, polarisation fractions, and absolute Faraday rotation measure. From 21 host galaxy properties considered, we find three correlated with tau, including the stellar surface density (or compactness; Pearson p-value p = 0.002 and Spearman p = 0.010), mass-weighted age (Spearman p-value p = 0.009), and a weaker correlation with gas-phase metallicity (Spearman p = 0.017). Weakly significant correlations are also found with Halpha equivalent widths and gravitational potential. From 10,000 trials of reshuffled datasets, we expect 2 strong Spearman correlations only 2% of the time, and three weaker correlations in 6.6% of cases. Compact host galaxies may have more ionised content which scatters the FRB further. No correlation is seen with host galaxy inclination, which weakens the case for an inclination bias, as previously suggested for samples of localised FRBs. A strong (p = 0.002) correlation is found for absolute rotation measure with optical disc axis ratio b/a; greater rotation measures are seen for edge-on host galaxies. Further high-time resolution FRB detections, coupled with localisation and detailed follow-up on their host galaxies, are necessary to corroborate these initial findings and shed further light into the FRB mechanism.

We introduce a novel dynamical model, named empirical triaxial orbit-superposition model, for the Milky Way halo. This model relies on minimal physical assumptions that the system is stationary, meaning the distribution function in 6D phase-space does not change when the stars orbiting in the correct gravitational potential. We validate our method by applying it to mock datasets that mimic the observations of the Milky Way halo from LAMOST + Gaia with stars' 3D position and 3D velocity observed. By removing the stellar disk and substructures, correcting the selection function, we obtain a sample of smooth halo stars considered as stationary and complete. We construct a gravitational potential including a highly flexible triaxial dark matter halo with adaptable parameters. Within each specified gravitational potential, we integrate orbits of these halo stars, and build a model by superposing the orbits together taking the weights of stars derived from the selection function correction. The goodness of the models are evaluated by comparing the density distributions as well as 3D velocity distributions numerically represented in the model to that in the data. The shape and radial density distribution of the underlying dark matter halo can be constrained well simultaneously. We apply it to three mock galaxies with different intrinsic shapes of their dark matter halos and achieved accurate recovery of the 3D dark matter density distributions for all.

Context: A large number of isolated black holes (IBHs) are expected to populate the Galaxy. However, only one has been confirmed by the analysis of a microlensing event, and no confirmed emission detection from an IBH has been reported so far. Aims: We analysed the detectability of electromagnetic signatures from IBHs moving in the Galaxy. Methods: We considered accretion from the interstellar medium onto an IBH and assumed the formation of an outflow. We modelled the accretion process and the interaction of the outflow with the surrounding medium on large scales, including mechanical feedback on the accretion process. Furthermore, we also calculated the emission from three different regions: the accretion region, the radiation from the outflow medium interaction structure, and the emission of relativistic particles that diffuse in the surrounding medium. Results: Multiwavelength emission associated with Galactic IBHs can be detected in systems moving through a very dense medium. Thermal emission from accretion could be observed in the mid infrared and in hard X rays with current and forthcoming observatories. Thermal and non thermal emission from the outflow medium shock could also be detected in the radio and millimetre ranges. Moreover, detection of the emission from particles diffusing in a dense medium could be feasible in gamma rays. Applying our model to the IBH associated with the gravitational microlensing event MOA2011BLG191 OGLE2011BLG0462, we inferred that radio and infrared detection of the IBH is plausible. Also, we derived that IBHs could be modest Galactic cosmic ray contributors, potentially reaching a 1% contribution at 1 PeV. Finally, by extending our model to primordial black holes, we conclude that efficient leptonic acceleration in their outflow medium interactions would rule them out as a major dark matter component.

Space-based gravitational wave (GW) detectors, such as LISA, are expected to detect thousands of Galactic close white dwarf binaries emitting nearly monochromatic GWs. In this study, we demonstrate that LISA is reasonably likely to detect higher harmonic GW signals, particularly the $(l, |m|) = (3, 3)$ mode, from a limited sample of nearby close white dwarf binaries, even with small orbital velocities $v/c$ of order $10^{-3}$. The amplitudes of these post-Newtonian modes provide robust probes of mass asymmetry in such systems, making them valuable observational targets, especially in mass-transferring binaries. Long-term, coordinated detector operations will further improve the prospects for detecting these informative signals.

Understanding Titan's planetary boundary layer (PBL) -- the lowest region of the atmosphere influenced by surface conditions -- remains challenging due to Titan's thick atmosphere and limited observations. Previous modeling studies have also produced inconsistent estimates of surface temperature, a critical determinant of PBL behavior, often without clear explanations grounded in surface energy balance. In this study, we develop a theoretical framework and apply a three-dimensional dry general circulation model (GCM) to investigate how surface thermal inertia influences surface energy balance and temperature variability across diurnal and seasonal timescales. At diurnal timescales, lower thermal inertia surfaces exhibit larger temperature swings and enhanced sensible heat fluxes due to inefficient subsurface heat conduction. In contrast, at seasonal timescales, surface temperature variations show weak sensitivity to thermal inertia, as atmospheric damping tends to dominate over subsurface conduction. The PBL depth ranges from a few hundred meters to 1,000 m on diurnal timescales, while seasonal maxima reach 2,000--3,000 m, supporting the interpretation from a previous study that the Huygens probe captured the two PBL structures. Simulated seasonal winds at the Huygens landing site successfully reproduce key observed features, including near-surface retrograde winds and meridional wind reversals within the lowest few kilometers, consistent with Titan's cross-equatorial Hadley circulation. Simulations at the planned Dragonfly landing site predict shallower thermal PBLs with broadly similar wind patterns. This work establishes a physically grounded framework for understanding Titan's surface temperature and boundary layer variability, and offers a unified explanation of Titan's PBL behavior that provides improved guidance for future missions.

We demonstrate that Type Ia supernovae (SNe Ia) observations impose the critical constraint for resolving the Hubble tension through late-time expansion modifications. Applying the Fisher-bias optimization framework to cosmic chronometers (CC), baryon acoustic oscillations (BAO) from DESI DR2, Planck CMB, and Pantheon+ data, we find that: (i) deformations in $H(z \lesssim 3)$ (via $w(z)$ reconstruction) can reconcile tensions between CC, Planck, DESI BAO, and SH0ES measurements while maintaining or improving fit quality ($\Delta\chi^2 < 0$ relative to $\Lambda$CDM); (ii) In the neighborhood of Planck best-fit $\Lambda$CDM model, no cosmologically viable solutions targeting $H_0 \gtrsim 69$ satisfy SNe Ia constraints. MCMC validation confirms the maximum achievable $H_0 = 69.09\pm0.30$ ($\chi^2_{\rm BF} \approx \chi^2_{\Lambda\rm CDM}$) across all data combinations, indicating that the conflict between late-time $w(z)$ modifications and SNe Ia observations prevents complete resolution of the Hubble tension.

The unstable isotope $^{60}$Fe, with a half-life of 2.6 million years, is produced primarily in supernova explosions. The observed presence of $^{60}$Fe in cosmic rays and its detection in deep-sea crusts and sediments suggest two possible scenarios: either the direct acceleration of $^{60}$Fe from supernova ejecta or its enrichment in the circumstellar material surrounding supernova progenitors, which indicates cosmic ray production in clusters of supernovae. Focusing on the latter scenario, we consider an environment shaped by successive supernova explosions, reminiscent of the Local Bubble around the time of the most recent supernova explosion. We independently tracked the evolution of the $^{60}$Fe mass ratio within the Local Bubble using passive scalars. To investigate the spectra of protons and $^{60}$Fe, we explicitly modelled cosmic-ray acceleration and transport at the remnant of the last supernova by simultaneously solving the hydrodynamical equations for the supernova outflow and the transport equations for cosmic rays, scattering turbulence, and large-scale magnetic field, using the time-dependent acceleration code RATPaC. The main uncertainty in our prediction of the local $^{60}$Fe flux at about $pc=1$ GeV/nuc is the magnetic-field structure in the Local Bubble and the cosmic-ray diffusion beyond the approximately $100$ kyr of evolution covered by our study. We found that if the standard galactic propagation applies, the local $^{60}$Fe flux would be around 3% of that measured. If there is a sustained reduction in the diffusion coefficient at and near the Local Bubble, then the expected $^{60}$Fe flux could be up to 30% of that measured.

We investigate the breakdown of the separate-universe approximation during transitions in transient ultra-slow-roll inflation by analyzing the evolution of the comoving curvature perturbation ${\cal R}$ and its conjugate momentum $\Pi$. It is well known that spatial gradient terms lead to a failure of this approximation, particularly at the transition from slow-roll to ultra-slow-roll phase. We show that a similar breakdown also occurs during the second transition back to slow-roll when considering the evolution of $\Pi$. Interestingly, while the homogeneous solution for $\Pi$ accurately captures the dynamics across the first transition, it is the homogeneous solution for ${\cal R}$ that becomes valid across the second. Furthermore, we demonstrate that the spatial curvature term introduced in the extended $\delta N$ formalism of \cite{Artigas:2024ajh} can be interpreted as arising from the contribution of $\Pi$ to the energy density perturbation. Importantly, this modification of the local Hubble parameter is valid only when the first slow-roll parameter is both small and strictly constant.

We present the first radio--continuum detection of the circumstellar shell around the well-known WN8 type Wolf-Rayet star WR16 at 943.5\,MHz using the \ac{ASKAP} \ac{EMU} survey. At this frequency, the shell has a measured flux density of 72.2$\pm$7.2\,mJy. Using previous \ac{ATCA} measurements at 2.4, 4.8, and 8.64~GHz, as well as the \ac{EMU} observations of the star itself, we determine a spectral index of $\alpha\,=\,+0.74\pm0.02$, indicating thermal emission. We propose that the shell and star both exhibit thermal emission, supported by the its appearance in near-infrared and H$\alpha$ observations. The latest \textit{Gaia} parallax is used to determine a distance of 2.28$\pm$0.09\,kpc. This star is well-known for its surrounding circular nebulosity, and using the distance and an angular diameter of 8\farcm42, we determine the shell size to be 5.57$\pm$0.22~pc. We use the \textit{Gaia} \acp{PM} of WR16 to determine peculiar velocities of the star as $V_{\alpha}(pec) =$ --45.3$\pm$5.4\,\kms\ and $V_{\delta}(pec) =$ 22.8$\pm$4.7\,\kms, which indicates that the star is moving in a north-west direction, and translates to a peculiar tangential velocity to be 50.7$\pm$6.9\,\kms. We also use these \acp{PM} to determine the shell's origin, estimate an age of $\sim 9500\pm 1300$\,yr, and determine its average expansion velocity to be $280\pm40$\,\kms. This average expansion velocity suggests that the previous transitional phase is a \ac{LBV} phase, rather than a \ac{RSG} phase. We also use the measured flux at 943.5~MHz to determine a mass-loss rate of $1.753\times 10^{-5}~M_\odot~$yr$^{-1}$, and use this to determine a lower-limit on ionising photons of $N_{UV} > 1.406\times 10^{47}~$s$^{-1}$.

Using both analytical and numerical means, we demonstrate that linear stability analysis of a hydrodynamic stratified atmosphere or a 1D coronal loop model in non-adiabatic settings features a thermal continuum corresponding to highly localized eigenfunctions. This thermal continuum can be precomputed, involving the net heat-loss function and its partial derivatives, and is the generalization of the thermal instability introduced by~\citet{Parker1953}. We account for a thermal imbalance, directly affecting thermal instability growthrates. We present completely general equations that govern all eigenmodes, including non-adiabatically affected p- and g-modes of the stratified settings. We intend to clarify how linear thermal instability is relevant for solar loops that show spontaneous in-situ condensations, and eliminate recent confusion on specific isochoric routes to linear instability alongside other thermal instability channels. The thermal continuum, previously identified as a crucial ingredient in magnetohydrodynamic eigenmode spectra for coronal loops and atmospheres, drives multithermal aspects across our universe, such as forming solar coronal rain and prominences, or cold cloud creation in intracluster to interstellar medium environments.

In this proceedings, we review the importance and complexity of detecting and characterizing supermassive binary black holes using conventional techniques (spectroscopy, timing, imaging). We show how spectropolarimetric data can strengthen or discard binary black hole candidates on solid grounds, and present what kind of instrument could perform this task. Namely, a high spectral resolution polarimeter operating from the far-ultraviolet (FUV) to the near-infrared (NIR), such as the POLLUX prototype, would perfectly fit within the payload of the Habitable Worlds Observatory and help revolutionize the study of extragalactic objects.

Numerical simulations can follow the evolution of fluid motions through the intricacies of developed turbulence. However, they are rather costly to run, especially in 3D. In the past two decades, generative models have emerged which produce synthetic random flows at a computational cost equivalent to no more than a few time-steps of a simulation. These simplified models qualitatively bear some characteristics of turbulent flows in specific contexts (incompressible 3D hydrodynamics or magnetohydrodynamics), but generally struggle with the synthesis of coherent structures. We aim at generating random fields (e.g. velocity, density, magnetic fields, etc.) with realistic physical properties for a large variety of governing partial differential equations and at a small cost relative to time-resolved simulations. We propose a set of approximations applied to given sets of partial differential equations, and test the validity of our method in the simplest framework: 2D decaying incompressible hydrodynamical turbulence. We compare results of 2D decaying simulations with snapshots of our synthetic turbulence. We assess quantitatively the difference first with standard statistical tools: power spectra, increments and structure functions. These indicators can be reproduced by our method during up to about a third of the turnover time scale. We also consider recently developed scattering transforms statistics, able to efficiently characterise non-Gaussian structures. This reveals more significant discrepancy, which can however be bridged by bootstrapping. Finally, the number of Fourier transforms necessary for one synthesis scales logarithmically in the resolution, compared to linearly for time-resolved simulations. We have designed a multiscale turbulence synthesis (MuScaTS) method to efficiently short-circuit costly numerical simulations to produce realistic instantaneous fields.

JWST/NIRCam provides rest-frame near-IR photometry of galaxies up to $z=2.5$ with exquisite depth and accuracy. This affords an unprecedented view of the evolution of the UV-optical-near-IR color distribution and its interpretation in terms of the evolving dust attenuation, $A_{\mathrm{V}}$. We use the value-added data products (photometric redshift, stellar mass, rest-frame $U-V$ and $V-J$ colors, and $A_{\rm V}$) provided by the public DAWN JWST Archive. This data product derives from fitting the spectral energy distributions obtained from multiple NIRCam imaging surveys, augmented with pre-existing HST imaging data. Our sample consists of a stellar mass complete sample of $\approx 28,000$ $M_\star> 10^{9}~M_\odot$ galaxies in the redshift range $0.5<z<2.5$. The $V-J$ color distribution of star-forming galaxies evolves strongly, in particular for high-mass galaxies ($M_\star>3\times 10^{10}~M_\odot$), which have a pronounced tail of very red galaxies reaching $V-J> 2.5$ at $z>1.5$ that does not exist at $z<1$. Such red $V-J$ can only be explained by dust attenuation, with typical values for $M_\star \approx 10^{11}~M_\odot$ galaxies in the range $A_{\mathrm{V}}\approx 1.5-3.5$ at $z\approx 2$. This redshift evolution went largely unnoticed before because the photometric redshift estimates for the reddest ($V-J>2.5$), most attenuated galaxies has markedly improved thanks to the new, precise photometry. Despite the increased attenuation, $U-V$ colors across the entire mass range are slightly bluer at higher $z$. In conclusion, whereas the rest-frame UV-optical color distribution evolves remarkably little from $z=0.5$ to $z=2.5$, the rest-frame optical-near-IR color distribution evolves strongly, primarily due to a very substantial increase with redshift in dust attenuation for massive galaxies. (Abbr.)

Mitigating data gaps in Gamma-ray bursts (GRBs) light curves (LCs) holds immense value for its application in cosmological research because it provides more precise measurements of the parameter of interest of the two-dimensional Dainotti relation which is a relation among the end time of the plateau emission, Ta, its respective luminosity, La which is calculated from the fluxes at the end of the plateau, Fa. This study extends the work done by arXiv:2305.12126; arXiv:2412.20091v4 on the 545 GRB sample by introducing six different models: Deep Gaussian Process (DGP), Temporal Convolutional Network (TCN), Hybrid model of Convolutional Neural Network with Long Short-Term Memory (CNN-LSTM), Bayesian Neural Network (BNN), Polynomial Curve Fitting and Isotonic Regression. Our findings demonstrate that Isotonic Regression achieves the highest uncertainty reduction for all three parameters (36.3% for log Ta, 36.1% for log Fa, and 43.6% for {\alpha}) outperforming all the other models. The CNN- LSTM model shows consistent improvements across all GRB parameters with the lowest outlier rate for {\alpha} (0.550%), surpassing the performance of the LSTM model in arXiv:2412.20091v4. The DGP model offers reliable uncertainty reduction across all parameters and improves upon the single-layer GP baseline. These advancements are essential for using GRBs as theoretical model discriminators via the parameters of their LCs and standard candles in cosmology, investigating theoretical models, and predicting GRB redshifts through recent cutting-edge machine-learning analysis (arXiv:2411.10736,arXiv:2405.02263; arXiv:2410.13985).

In solar-like oscillators, acoustic waves are excited by turbulent motion in the convective envelope and propagate inward, generating a variety of standing pressure modes. When combining together the power of several solar acoustic modes, an excess not compatible with pure stochastic excitation was found in some studies. This could be the signature of a second mode excitation source. With over 27 years of helioseismic data from the Sun as a star observations by the Solar and Heliospheric Observatory (SoHO), we aim to study the variation in mode energy over this period, covering solar Cycles 23, 24, and the beginning of Cycle 25. We focus on the possible sources of high peaks in the mode-energy time series, i.e. instrumental problems or other exciting mechanisms, such as flares or Coronal Mass Ejections. We reconstruct the energy time series for each mode with a sampling time of 1.45 days. By combining the small-time-scale variations in energy for several low-degree modes in the 2090-3710uHz range, we study the correlation between the modes and their compatibility with the hypothesis that modes are only stochastically excited by convection. The observed excitation rate significantly deviates from what would be expected in the case of a purely stochastic excitation. Our results indicate that this energy excess cannot be only attributed to instrumental effects and does not exhibit a cyclic variation. Although high-energy excesses are occasionally associated with observations of flares or CMEs, no consistent pattern could be identified. The excitation is slightly more frequent for modes probing the upper layer of the convective zone. Furthermore, the energy supply rate seems to vary over time with the mean value following a modulation that can match the Quasi-Biennial Oscillation (QBO) observed in other solar indicators, and the variance being anti-correlated with the cycle.

In recent years supergranulation has emerged as one of the biggest challenges for the detection of Earth-twins in radial velocity planet searches. We used eight years of Sun-as-a-star radial velocity observations from HARPS-N to measure the quiet-Sun's granulation and supergranulation properties of most of its 11-year activity cycle, after correcting for the effects of magnetically active regions using two independent methods. In both cases, we observe a clear, order of magnitude variation in the time-scale of the supergranulation component, which is largest at activity minimum and is strongly anti-correlated with the relative Sunspot number. We also explored a range of observational strategies which could be employed to characterise supergranulation in stars other than the Sun, showing that a comparatively long observing campaign of at least 23 nights is required, but that up to 10 stars can be monitored simultaneously in the process. We conclude by discussing plausible explanations for the "supergranulation" cycle.

The Large Magellanic Cloud (LMC), being a nearby and actively star-forming satellite galaxy of the Milky Way, is an ideal site to observe the multi-phase interstellar medium (ISM) of a galaxy across the electromagnetic spectrum. We aim to exploit the available SRG/eROSITA all-sky survey data to study the distribution, composition and properties of the diffuse X-ray emitting hot gas in the LMC. We construct multi-band X-ray images of the LMC, reflecting the morphology and temperatures of the diffuse hot gas. By performing spatially resolved X-ray spectroscopy of 175 regions, we constrain the distribution, physical state, and composition of the hot ISM phase throughout the LMC. We combine our constraints with multiwavelength data to obtain a comprehensive view of the different ISM phases. We measure a total X-ray luminosity of the hot ISM phase of $1.9\times10^{38}\,\mathrm{erg\,s^{-1}}$ ($0.2-5.0\,\mathrm{keV}$ band), and constrain its thermal energy to around $5\times10^{54}\,\mathrm{erg}$. The typical density and temperature of the X-ray emitting plasma are around $5\times10^{-3}\,\mathrm{cm^{-3}}$ and $0.25\,\mathrm{keV}$, respectively, with both exhibiting broad peaks in the southeast of the LMC. The observed degree of X-ray absorption correlates strongly with the distribution of foreground HI gas, whereas a spatial anticorrelation between the hot and cold ISM phases is visible on sub-kpc scales within the disk. The abundances of light metals show a strong gradient throughout the LMC, with the north and east exhibiting a strong $\alpha$-enhancement, as expected from observed massive stellar populations there. In contrast, the enigmatic ``X-ray spur'' exhibits a local deficit in $\alpha$-elements, and a peak in hot-gas pressure at $P/k\sim10^5\,\mathrm{K\,cm^{-3}}$, consistent with a dominant energy input through tidally driven gas collisions.

The NEID Earth Twin Survey (NETS) has been delivering a rich set of precise radial velocity (RV) measurements for 41 bright, nearby main sequence stars. Here, we describe the status of the survey after three years on sky and we present the full set of RV measurements and accompanying stellar activity indicators. We discuss intermediate survey diagnostics, including calibration of the known RV zero point offset introduced following the Contreras fire in 2022 and the identification of an undiagnosed and previously unknown zero point offset in 2021. An analysis of our data set using exttt{RVSearch} demonstrates that for these target stars, NEID is independently sensitive to nearly all known planets with periods shorter than the NETS observing baseline. We also highlight a number of newly detected RV signals, which present exciting opportunities for future investigations.

The phase diagrams of ultralight dark matter (DM), modeled as a self-gravitating Bose-Einstein condensate with axionlike interaction, are studied. We classify stable, metastable, and unstable DM states over a wide range of condensate wave function amplitudes. It is shown that the axionlike interaction causes instability and an imaginary speed of sound at low amplitudes, whereas, in a specific high-amplitude band, DM attains a stable state capable of forming a dense solitonic core and suppressing quantum fluctuations in the surrounding galactic DM halo. These findings are corroborated by evaluating thermodynamic functions for DM in the dwarf galaxy NGC2366 and its hypothetical analogs with different core-to-halo mass ratios. Distinct DM phase compositions respond differently to fluctuation-induced partial pressure, resulting in a first-order phase transition in a certain range of an interaction parameter. While the DM properties in NGC2366 lie within the supercritical regime, the phase transition nonetheless provides a thermodynamic marker separating stable from unstable DM configurations. Once a dense core forms - reaching a threshold of about 12% of the total mass - the enhanced gravitation stabilizes the DM halo against fluctuations, while the internal pressure ensures core stability. In particular, we find that NGC2366's dense DM comprises roughly 19% of the DM mass while occupying only 4.7% of its total volume.

Observational manifestations of disturbances in a protoplanetary disk caused by a collision with a massive planet are studied. It is assumed that the planet moves along a parabolic trajectory that intersects the disk plane near the star. Gas-dynamic simulation is performed using the finite volume method on a long time scale. On its basis, images of the disk observed from the pole and edge-on are constructed in the infrared and submillimeter ranges. A wide range of planet orbit parameters is considered. The approach of the planet was considered both prograde and retrograde with the respect to the disk rotation. Calculations have shown that in the images of the disk seen pole-on, two spiral arms can be observed in case of the prograde fall, and one with retrograde case. In the case of observations of a disk whose plane is inclined at a small angle to the line of sight, distortions of the disk plane can be noticeable. In addition, a gas tail is extended from the disk in the direction of the planet's motion, which can also be identified in observations.

We investigate the impact of statistical inhomogeneity and anisotropy in primordial scalar perturbations on the scalar-induced gravitational waves (SIGW). Assuming inflationary quantum fluctuations originate from a coherent state, the resulting primordial scalar perturbations acquire a non-zero space-dependent mean, violating statistical homogeneity, statistical isotropy, and parity. As a consequence of statistical inhomogeneities, SIGW acquires distinct scale-dependent features in its correlation function. Statistical anisotropies further lead to possible parity violation and correlation between different polarization modes in the tensor perturbations. Therefore, detection of these signatures in the stochastic gravitational wave background would offer probes to the statistical nature of primordial scalar perturbations beyond the scales accessible to CMB observations.

We present a comprehensive thermodynamic derivation and observational validation of the third-order galaxy correlation function, extending the seminal work of Iqbal et al [ 1, 2]. Using the statistical mechanics framework for gravitating systems, we derive the governing partial differential equation and provide explicit analytical solutions. Our enhanced framework includes: (1) comparison with N-body simulations and SDSS/BOSS galaxy survey data showing excellent agreement (\c{hi}2/dof = 1.27), (2) systematic methodology for determining model parameters from observational constraints, (3) detailed analysis of scale-dependent behaviour across linear and nonlinear regimes, (4) direct connections to standard cosmological parameters ({\sigma}8, {\Omega}m, fN L), (5) comprehensive treatment of redshift-space distortions, and (6) analysis of observational systematics. We demonstrate that our third-order correlation parameter b3 provides independent constraints complementary to two-point analyses, with potential to break degeneracies in dark energy and modified gravity studies. Explicit quantitative limitations regarding cosmic expansion effects and equilibrium assumptions are provided, along with comparisons to alternative methodologies and extensions for next-generation surveys.

The persistent radio source (PRS) associated with FRB 20121102A, the first precisely localized repeating fast radio burst (FRB), provides key constraints on both its local environment and the nature of the underlying FRB engine. We present a seven-year (2016-2023) temporal analysis of the PRS, combining new uGMRT observations with archival data across L-band frequencies. We find no statistically significant long-term trend in its L-band flux density. The observed variability is consistent with refractive interstellar scintillation, and the data do not require the PRS to be a source exhibiting strong intrinsic variability. This stability challenges models predicting rapid fading from evolving magnetized outflows, such as those powered by young magnetars or relativistic shocks. Our low-frequency observations show no evidence for spectral evolution between 1.4 GHz and 745 MHz, with a measured spectral index of $\alpha = -0.15 \pm 0.08$, in agreement with values reported from earlier observations in 2016-2017. The PRS remains compact, exhibits a flat radio spectrum, and-if powered by an intermediate-mass black hole accreting at a low Eddington ratio-its radio and X-ray properties are broadly consistent with the fundamental plane of radio-loud AGNs. While not conclusive, this scenario represents a viable alternative to magnetar wind nebula models and warrants further investigation. Furthermore, we find no statistically significant correlation between FRB burst activity and the luminosity of associated PRSs among repeating sources. This apparent decoupling challenges simple progenitor models that directly link bursts and persistent emission. Together, these results suggest that the FRB engine and PRS may in some systems originate from physically distinct sources, underscoring the need for flexible models to explain the diverse environments of repeating FRBs.

The inference of the magnetic field vector from spectropolarimetric observations is crucial for understanding the physical processes governing the solar corona. We investigate which information on the magnetic fields of coronal bright points (CBP) can be gained from the intensity and polarization of the Fe XIII 10747 A, Fe XIV 5303 A, Si X 14301 A, and Si IX 39343 A forbidden lines. We apply the P-CORONA synthesis code to a CBP model in the very low corona, obtained with the Bifrost code, and to a larger global model to study the impact of the outer coronal material along the line of sight (LoS). The enhanced density within the CBP produces an intensity brightening, but suppresses the linear polarization. The circular polarization from such regions often approaches 0.1% of the intensity. The contribution from the coronal material along the LoS depends strongly on its temperature, and is weaker for lines with a peak response at higher temperatures (Fe10747 at 1.7 MK; Fe5303 at 2 MK). The weak field approximation (WFA) provides information on the longitudinal magnetic fields in the strongest-emitting spatial intervals along the LoS, and is more reliable in the regions of the CBP where the field does not change sign. This tends to coincide with the regions where there is a strong correlation between the circular polarization and the wavelength derivative of the intensity. Considering roughly 30 minutes of time evolution, the CBP signals are somewhat attenuated but are still identifiable, while the area where the WFA can be suitably applied remains substantial. The circular polarization of the Fe5303 and especially Fe10747 lines are valuable diagnostics for the magnetic fields in the higher-temperature regions of the CBP, which could be exploited with future coronagraphs with similar capabilities to Cryo-NIRSP/DKIST, but designed to observe below 1.05 solar radii.

Modeling the broadband emission of blazars has become increasingly challenging with the advent of multimessenger observations. Building upon previous successes in applying convolutional neural networks (CNNs) to leptonic emission scenarios, we present an efficient CNN-based approach for modeling blazar emission under proton synchrotron and hybrid lepto-hadronic frameworks. Our CNN is trained on extensive numerical simulations generated by SOPRANO, which span a comprehensive parameter space accounting for the injection and all significant cooling processes of electrons and protons. The trained CNN captures complex interactions involving both primary and secondary particles, effectively reproducing electromagnetic and neutrino emissions. This allows for rapid and thorough exploration of the parameter space characteristic of hadronic and hybrid emission scenarios. The effectiveness of the trained CNN is demonstrated through fitting the spectral energy distributions of two prominent blazars, TXS 0506+059 and PKS 0735+178, both associated with IceCube neutrino detections. The modeling is conducted under assumptions of constant neutrino flux across distinct energy ranges, as well as by adopting a fitting that incorporates the expected neutrino event count through a Poisson likelihood method. The trained CNN is integrated into the Markarian Multiwavelength Data Center (MMDC; this https URL), offering a robust tool for the astrophysical community to explore blazar jet physics within a hadronic framework.

Avalanche models running in a self-organized critical regime have proven powerful in reproducing the power-law distributions and scale invariance that characterize the statistical properties of solar flares. They are often interpreted as representing an individual active region of the Sun. As a result, this class of models has rarely been applied to describe sympathetic flares $\unicode{x2014}$ solar eruptions that occur in close spatial and temporal proximity, seemingly driven by their mutual interaction. In this study, we investigate the phenomenon of sympathetic flaring using avalanche models and compare their statistical properties with observations of sympathetic flares on the Sun. We developed a novel avalanche model featuring two connected lattices, each representing a distinct active region. This connectivity allows the transfer of nodal variable between the lattices, simulating the non-local effects expected to occur during sympathetic flares. Our results show that under strong connectivity, the lattices exhibit temporal synchronization, with correlations between their avalanche energies. Furthermore, increasing the connectivity between the lattices results in an excess of avalanches at short waiting times. A quantitative comparison with observational data suggests that only a weak connectivity allows our model to replicate the observed solar waiting time distributions. Consequently, we propose that if magnetic connectivity between distinct active regions drives sympathetic flaring on the Sun, it must remain relatively weak.

The low equilibrium temperatures of temperate sub-Neptunes lead to extremely long chemical timescales in their upper atmospheres, causing the abundances of chemical species to be strongly shaped by atmospheric transport. Here, we used a three-dimensional (3D) general circulation model involving a passive tracer to investigate the atmospheric circulation and 3D transport of temperate gas-rich sub-Neptunes, using K2-18b as an example. We model K2-18b as a synchronous or asynchronous rotator, exploring spin-orbit resonances (SOR) of 2:1, 6:1, and 10:1. We find that the strong absorption of CO$_2$ and CH$_4$ induces a detached convective zone between 1 and 5 bar, resulting in strong vertical mixing at these levels. The upper atmosphere is dominated by eastward winds (an equatorial superrotating jet present in all simulations), leading to warmer evening terminators and approximately 20% higher passive tracer mass mixing ratios compared to the morning terminators. Rotation rates have minimal impact on the strength of global mean vertical mixing, but significantly influence the latitudinal distribution of passive tracers. For synchronous, 2:1 SOR, and 6:1 SOR simulations, passive tracers are more abundant in the upwelling branches at latitudes within 60$^\circ$, while for the 10:1 SOR simulation, strong transient eddies at high latitudes (>70$^\circ$) between 0.1 to 1 bar can transport passive tracers upward from the deep atmosphere, making them more abundant there, despite their alignment with the downwelling branch of the large-scale circulation. This study focuses on the atmospheric dynamics and its influence on passive tracer transport, while a follow-up paper will incorporate active chemical species.

Active Galactic Nuclei (AGN) are believed to be powered by accretion of matter onto a supermassive black hole. A fundamental ingredient in shaping our understanding of AGN is their variability across the entire electromagnetic spectrum. Variability studies have the potential to help us understand the geometry of the emitting regions (in various energy bands), their causal relations, and the physics of the accretion processes. This review focuses on the observational properties of AGN variability in the optical/UV/X-ray bands (where most of the AGN luminosity is emitted) and their dependence on the AGN physical parameters (i.e. mass, luminosity, accretion rate). We also discuss possible interpretations in the context of accreting compact systems, and we review the use of variability as a tool to discover AGN and trace their properties across cosmic time, using both ground and space facilities. Finally, we discuss the opportunities and challenges provided by current and next-generation optical/X-ray surveys, to use variability as an effective tool to probe the growth of super massive black holes in the Universe.

[abridged] 1H 0419-577 is a highly-accreting, luminous BLS1 AGN. This study aims to characterise its disc-corona system using, for the first time, simultaneous XMM-Newton and NuSTAR observations, performed in May and November 2018. We conducted high-resolution grating spectroscopy to identify potential soft X-ray absorption and emission features. To measure the hot corona temperatures from the spectral analysis above 3 keV, we also included data from a previous NuSTAR observation from June 2015. We characterised the disc-corona system properties by analysing the broadband spectra and the SED from UV to hard X-rays. 1H 0419-577 was observed in a bare-like high-flux state at both epochs, with negligible neutral and ionised absorption along its line of sight at both Galactic and AGN rest-frames. However, several soft X-ray emission lines were detected, notably a broad and intense OVII line indicating an accretion disc origin at only a few tens of gravitational radii. The broadband X-ray spectra revealed a prominent, absorption-free smooth soft X-ray excess, a weak Fe Kalpha complex, and a lack of a Compton hump. Fitting data above 3 keV yielded apparent moderate hot corona temperatures of ~20-30 keV for the 2018 and 2015 observations, depending on the model applied. The 2018 X-ray broadband spectra were well reproduced by either a relativistic reflection model with a high-density accretion disc (~10^18 cm^-2), or a hybrid model combining warm and hot coronae with relativistic reflection. We performed the SED analysis for the latter scenario, which indicated that both the hot and warm coronae would have a small spatial extent. Both scenarios can successfully reproduce the two 2018 observations of 1H 0419-577, but they imply very different physical conditions, for example, in terms of disc density, temperature and accretion power released in the hot corona and the origin of the UV emission.

The center of our galaxy hosts Sagittarius~A*, a supermassive compact object of $\sim 4.3\times 10^6$ solar masses, usually associated with a black hole. Nevertheless, black holes possess a central singularity, considered unphysical, and an event horizon, which leads to loss of unitarity in a quantum description of the system. To address these theoretical inconsistencies, alternative models, collectively known as exotic compact objects, have been proposed. In this paper, we investigate the potential detectability of signatures associated with non-rotating exotic compact objects within the Sgr~A* polarized flares dataset, as observed through GRAVITY and future instruments. We examine a total of eight distinct metrics, originating from four different categories of static and spherically symmetric compact objects: Black Holes, Boson stars, Fluid spheres, and Gravastars. Our approach involves utilizing a toy model that orbits the compact object in the equatorial plane, at the Schwarzschild-Keplerian velocity. Using simulated astrometric and polarimetric data with present GRAVITY and future GRAVITY+ uncertainties, we fit the datasets across all metrics examined. We evaluated the detectability of the metric for each dataset based on the resulting $\chi^2_\mathrm{red}$ and BIC-based Bayes factors. Plunge-through images of ECOs affect polarization and astrometry. With GRAVITY's present uncertainties, only a compact boson-star model is discernible. GRAVITY+'s improved sensitivity allows detection of most exotic compact object models. However, enhancing the astrophysical complexity of the hot spot model diminishes these outcomes. Presently, GRAVITY's uncertainties limit us to detecting just one exotic compact object metric. With GRAVITY+'s enhanced sensitivity, we can expect to uncover additional exotic compact object models and use Sgr~A* as a laboratory for fundamental physics.

The Vera C. Rubin Observatory recently released Data Preview 1 (DP1) in advance of the upcoming Legacy Survey of Space and Time (LSST), which will enable boundless discoveries in time-domain astronomy over the next ten years. DP1 provides an ideal sandbox for validating innovative data analysis approaches for the LSST mission, whose scale challenges established software infrastructure paradigms. This note presents a pair of such pipelines for variability-finding using powerful software infrastructure suited to LSST data, namely the HATS (Hierarchical Adaptive Tiling Scheme) format and the LSDB framework, developed by the LSST Interdisciplinary Network for Collaboration and Computing (LINCC) Frameworks team. This note presents a pair of variability-finding pipelines built on LSDB, the HATS catalog of DP1 data, and preliminary results of detected variable objects, two of which are novel discoveries.

This work aims to demonstrate that two arrays of Lumped Element Kinetic Inductance Detectors (LEKIDs), when employed in filled array configuration and separated by an external linear polarizer oriented at 45 degrees, can achieve the precision required by next-generation cosmological experiments. The focus here is on validating their ability to meet stringent uncertainty requirements, in particular for polarization angle reconstruction. To achieve this, the uncertainties in the reconstruction of the polarization angle have been characterized in the laboratory using a dedicated closed-circuit 100 mK dilution cryostat. This is optically coupled to a Martin-Puplett interferometer and a custom-designed sky simulator equipped with both photometric and polarized sources, allowing one to reproduce realistic ground-based observation conditions. This experimental setup allows us to generate intensity and polarization maps with diffraction-limited resolution, allowing us to determine the polarization angles and their associated uncertainties. The results show performance in line with expectations for the next generation CMB experiments. The polarization angle was reconstructed with an uncertainty of 6.5 arcmin.

We forecast the impact of direct radial acceleration measurements, based on two redshift measurements of the same target one decade apart, on constraining the Yukawa correction to the Newtonian potential in the Milky Way. The Galaxy's matter distribution is modeled as the sum of a spherical bulge, a spherical dark matter halo, and two axially symmetric disks. Considering a sample of 165 Milky Way globular clusters, we find that the precision of next-generation spectrographs ($\sim$ 10 cm s$^{-1}$) is not sufficient to provide competitive constraints compared to rotation curve data using the same baryonic matter distribution. The latter sample only becomes competitive for a precision better than 0.6 cm s$^{-1}$. On the other hand, we find that adopting a population of $1.3 \times 10^5$ RR Lyrae stars as targets, a precision of $\sim$ 10 cm s$^{-1}$ can achieve constraints on the Yukawa parameters as strong as with the rotation curves.

The NASA Nancy Grace Roman Space Telescope (Roman) will be the first survey able to detect large numbers of both cold and hot exoplanets across Galactic distances: ~1,400 cold exoplanets via microlensing and ~200,000 hot, transiting planets. Differing sensitivities to planet bulk properties between the microlensing and transit methods require relations like a planet mass-radius relation (MRR) to mediate. We propose using instead a planet-host ratio relation (PHRR) to couple directly microlensing and transit observables. Unlike the MRR, a PHRR uses parameters that are always measured and so can leverage the full Roman exoplanet sample. Using 908 confirmed exoplanets from the NASA Exoplanet Archive, we show that transit depth, $\delta$, and planet-host mass ratio, $q$, obey a PHRR that is continuous over all planet scales. The PHRR is improved by including orbital period, $P$, and host effective temperature, $T_{\star}$. We compare several candidate PHRRs of the form $\delta (q,T_\star, P)$, with the Bayesian Information Criterion favouring power-law dependence on $T_\star$ and $P$, and broken power-law dependence on $q$. The break in $q$ itself depends on $T_\star$, as do the power-law slopes in $q$ either side of the break. The favoured PHRR achieves a fairly uniform $50\%$ relative precision in $\delta$ for all $q$. Approximately $5\%$ of the sample has a transit depth that is strongly under-predicted by the PHRR; around half of these are associated with large stars ($R_\star > 2.5 \, R_{\odot}$) potentially subject to Malmquist bias.

Determining stellar ages is challenging, as it depends on other stellar parameters in a non-linear way and often relies on stellar evolution models to infer the underlying relation between these parameters and age. This complexity increases when comparing different age-dating methods, as they rely on distinct indicators and are often applicable to non-overlapping regions of the color-magnitude diagram. Moreover, many empirical calibration methods rely on pre-determined ages, often from open clusters or asteroseismology, which only cover a limited parameter space. Fortunately, the age-velocity-dispersion relation (AVR), in which the velocity dispersion increases with age, is a universal feature among stars of all evolutionary stages. In this paper, we 1) explore the parameter space in which [C/N] and gyrochronology are applicable, extending beyond the domains probed by asteroseismology and open clusters, and 2) assess whether the traditionally assumed [C/N] and gyrochronology relations yield ages on a consistent physical scale, after calibrating both using the same AVR. We find gyrochronology can be applied to all partially convective stars after they have converged onto the slow rotating sequence and before they experience weakened magnetic braking; [C/N] can be used to infer ages for all giants with metallicity > -0.8 dex and [C/N] < -0.05 dex, and can be used as an age-indicator down to [Fe/H] of -1 dex if only selecting the low-$\alpha$ disk. Lastly, ages obtained from [C/N] and gyrochronology agree within uncertainty after accounting for systematic offsets.

Precise measurement of the Hubble parameter will enable stringent tests of the standard model for cosmology. Standard sirens, using the luminosity distances measured by gravitational-wave observations of compact binary mergers, are expected to provide such measurements independently in the next decade. With the ground- and space-based gravitational wave observatories, the LIGO-Virgo-KAGRA (LVK) network and the Laser Interferometer Space Antenna (LISA), different types of standard sirens altogether will place constraints across a wide redshift range. In this paper, we forecast the precisions of standard siren Hubble parameter measurements and compare various scenarios, accounting for the dominant sources of systematic uncertainty. Specifically, we find a $2\%$ constraint on $H_0$, a $1.5-3\%$ constraint on $H(z)$ at $z=1$, and a $3-5\%$ constraint on $H(z)$ at $z=7$ when combining LVK and LISA standard sirens with precise redshift measurements from electromagnetic counterpart observations. We do not find a significant improvement when including standard sirens with no EM counterpart, but which rely on features in the black hole mass distribution, and the potential systematics introduced by the possible redshift evolution of such features could further degrade the measurement accuracy if not properly accounted for.

We present a new constraint on the primordial helium abundance, $Y_\mathrm{P}$, based on Subaru observations. A major source of uncertainty in previous $Y_\mathrm{P}$ determinations is the lack of extremely metal-poor galaxies (EMPGs; $0.01-0.1Z_\odot$), which have metallicities a few to ten times lower than the metal-poor galaxies (MPGs; $0.1-0.4Z_\odot$) predominantly used in earlier studies, requiring substantial extrapolation to zero metallicity. Here, we perform Subaru near-infrared spectroscopy of 29 galaxies, including 14 EMPGs. By incorporating existing optical spectra, we derive He/H for each galaxy using photoionization modeling of helium and hydrogen emission lines, including the He I 10830 Å line to break the density-temperature degeneracy. After carefully selecting galaxies with robust He/H determinations, and adding 58 galaxies from previous studies, we obtain $Y_\mathrm{P} = 0.2387^{+0.0036}_{-0.0031}$. This $Y_\mathrm{P}$ value is $\sim1\sigma$ lower than most of the previous estimates, but agrees with recent determinations using EMPGs and the CMB constraint from the Atacama Cosmology Telescope (ACT) experiment. Our result indicates $N_\mathrm{eff} = 2.54^{+0.21}_{-0.20}$, showing a mild ($1-2\sigma$) tension with the Standard Model and Planck results. These tensions may suggest a nonzero lepton asymmetry, which would be alleviated with $\xi_\mathrm{e} = 0.05^{+0.02}_{-0.02}$. More observations of EMPGs and further assessments of systematic uncertainties are essential to test the potential tension more rigorously.

In this letter, we report a new phenomena of recurring regions when relating observables for hybrid neutron stars and hybrid neutron-star mergers. To describe dense matter within hybrid stars, we introduce a percolation to vary the size and characteristics of the deconfinement phase transition to quark matter. Before and after the percolation, we keep the hadronic and quark phases the same, described by different realistic models for the equation of state of beta-equilibrated, charge-neutral, zero-temperature matter. When solving spherical and deformed equations for neutron stars in general relativity, we find that: no matter the size or characteristics of the percolation region, or the order of the phase transition on either side (hadronic and quark), as long as we minimize the average sound speed from the beginning of the percolation region to the central density for a given star, we can produce equations of state that cross through the same, small recurring region within mass-radius and mass-tidal deformability diagrams. Our findings provide a new way to produce hybrid equations of state for dense matter that match a given observation of neutron stars or neutron star mergers.

Galaxies form due to the continuous accretion of gaseous material towards the centre of dark matter haloes, which gives rise to stellar discs once the gas reaches high density. Analysis of stellar age distribution shows that stars in the inner regions of the disc are born at earlier times, in contrast to those found in the outer regions, which constitute a younger population. This scheme, typical of galactic formation processes, is known as the "inside-out" scenario and is closely related to the global properties of a significant fraction of the population of spiral galaxies. In this study, we analyse model galaxies from the Auriga project, a set of high-resolution, magnetohydrodynamic cosmological simulations of Milky Way-mass galaxies, to investigate the distribution of metals at $z=0$. Since different galactic components were subjected to different formation histories, we focus on the distribution of metals in the halo, bulge, and disc using a dynamical decomposition method based on the degree of rotation of the stars and their gravitational potential.

Recent observations have revealed an intriguing abundance of polar-orbiting Neptune-sized planets, many of which exhibit unusually inflated radii. While such misaligned orbits point to a complex dynamical history, the connection between their orbital orientations and planetary structures remains poorly understood. In this study, we analyze a sample of 12 misaligned and 12 aligned planets using structure models that incorporate tidal heating. We use various statistical tests to demonstrate with at least $90\%$ confidence that misaligned planets experience more tidally-induced radius inflation compared to aligned planets. This inflation likely stems from their dynamically active histories, which often place them in close-in, eccentric, and highly inclined orbits. We further present a case study of WASP-107~b, an exceptionally inflated polar Neptune, and model its history using a simple coupled orbital and radius evolution approach. Our results place constraints on the planet's tidal quality factor that agree with recent JWST observations.

The neutrino telescopes have firmly established the existence of ultra-high-energy (UHE) neutrinos. Observations of these neutrinos offer a unique probe of neutrino self-interactions. This work investigates how the self-interacting neutrinos, mediated by scalar bosons, inject energy into the medium through radiative scattering with the cosmic neutrino background, leaving an imprint on the Cosmic Microwave Background (CMB) spectrum. The energy injection into plasma in redshift ranges, $5\times10^4\lesssim z\lesssim2\times10^6$ and $ z\lesssim5\times10^4$, leads to $\mu$-type and $y$-type of CMB spectral distortions, respectively. Using observational constraints from Cosmic Background Explorer/Far Infrared Absolute Spectrophotometer (COBE/FIRAS) and projected sensitivities from Primordial Inflation Explorer (PIXIE) experiments for $\mu$-type and $y$-type CMB distortions, we derive the stringent upper bounds on the self-interaction coupling strength as a function of mediator mass for neutrino interactions. We focus on flavor-specific self-interaction related to muon neutrinos and sub-GeV mass mediators ($m_{\phi}$). We find the upper bound on the self-interaction coupling strength to be $\sim 3.1\times 10^{-5}$ for the muon neutrino, considering UHE muon neutrino energy to be 1 PeV and PIXIE projected upper bounds on $\mu$-type CMB spectral distortion. The bound remains constant till the mediator mass reaches the center of mass energy, and after that, it gets relaxed and becomes proportional to the mediator mass. We have also compared our results with existing bounds in the literature. Our findings indicate that CMB spectral distortion could play a decisive role in exploring neutrino physics beyond the standard model of particle physics, and future missions like PIXIE can provide valuable insights.

We use [CII] observations of a large QSO sample to segregate sources by host galaxy mass, aiming to identify those in the most massive hosts. [CII] luminosity, a known tracer of molecular gas, is taken as a proxy for host mass and used to rank 190 QSOs at z>5.9, spanning a 6-mag UV luminosity range (-22<Muv<-28). Particularly valuable are ALMA data from a cycle-10 CISTERN program, providing [CII] coverage for 46 UV-faint (M_{UV,AB}>-24.5) and 25 especially UV-faint (Muv>-23.5) QSOs, improving statistics by 5x and 6x, respectively. Taking massive host galaxies to be those where L[CII]>1.8x10^9 Lsol (median L[CII] of UV-bright QSOs), we identify 61 QSOs, including 13 which are UV-faint and 7 especially UV-faint. Using these selections and recent QSO luminosity functions (LFs), we present the first characterization of UV luminosity distribution for QSOs in massive host galaxies and quantify [CII] LFs for both UV-bright and UV-faint QSOs. While ~3% of massive-host QSOs are UV-bright (Muv<-26), >~85% are UV-faint (Muv>-24.5). This wide dispersion in UV luminosities reflects variations in dust obscuration, accretion efficiency, and black hole mass. Though spectroscopy is needed for definitive conclusions, black hole mass appears to be the dominant factor driving variations in the UV luminosity, based on 34 [CII]-luminous (L[CII]>1.8x10^9 Lsol) QSOs distributed across a ~3-mag baseline in UV luminosity and with measured MBH. At Muv~-23, the median extrapolated log10 (MBH/Msol) is 8.1+/-0.4, consistent with the local relation. SMBHs in UV-bright QSOs thus appear to be ~15(-9)(+25)x more massive than typical for massive host galaxies at z~6.

Since the discovery of the first exoplanet orbiting a Sun-like star, the confirmation of nearly 6000 exoplanets to date - and their diversity - has revolutionized our knowledge of planetary systems in the past three decades. Nevertheless, the majority of these planets are around mature stars (${\gtrsim}1$ Gyr), where the planet birth environments have already dissipated. Indeed, we have only confirmed 2 forming planets (i.e., protoplanets; ${\lesssim}10$ Myr) residing in one single system. In comparison, we have imaged over 200 protoplanetary disks in the past decade, with many of them hosting substructures such as spirals and gaps which suggest the existence of protoplanets. To understand the early stages of planet formation, the Habitable Worlds Observatory (HWO) - with its high-contrast imaging and integral field spectroscopy capabilities - presents a unique opportunity to explore the demographics of the natal stages of planet formation and their birth environments. We propose to image protoplanets within substructured protoplanetary disks using HWO via direct imaging, and characterize them (i.e., protoplanets, protoplanetary disks, circumplanetary disks) using integral field spectroscopy and spectropolarimetry. This effort will dramatically extend current population of protoplanets, probing and characterizing over 200 protoplanets. By expanding the number of protoplanets by two orders of magnitude, these observations will test and refine planet formation theory and planet-disk interaction theory, and further motivate planet migration studies together with existing mature planets. The results will offer critical insight into planetary system formation and evolution, and help understand the origin of our own Solar System.

We investigate observational constraints on cubic curvature corrections to general relativity by analyzing quasi-periodic oscillations (QPOs) in accreting black hole systems. In particular, we study Kerr black hole solution corrected by cubic curvature terms parameterized by $\beta_5$ and $\beta_6$. While $\beta_6$ corresponds to a field-redefinition invariant structure, the $\beta_5$ term can in principle be removed via a field redefinition. Nonetheless, since we work in the frame where the accreting matter minimally couples to the metric, $\beta_5$ is in general present. Utilizing the corrected metric, we compute the QPO frequencies within the relativistic precession framework. Using observational data from GRO J1655$-$40 and a Bayesian analysis, we constrain the coupling parameters to $|\beta_5| < (27.73\,\mathrm{km})^4$ and $|\beta_6| < (15.34\,\mathrm{km})^4$. These bounds improve upon existing constraints from big-bang nucleosynthesis and the speed of gravitational waves.

We propose a multi-field fiber inflation scenario in type IIB perturbative large volume compactifications, showing that the multi-field dynamics suppresses trans-Planckian displacements of the canonical inflaton. Considering a concrete K3-fibred Calabi-Yau (CY) threefold with $h^{1,1}({\rm CY})=3$ and having certain underlying symmetries, we show that the presence of multi-fibre moduli creates an assisted inflation scenario where multiple moduli collectively help in producing the cosmological observables consistent with the experimental observations. We further argue that individual field range excursions $(\Delta\phi_n)$ corresponding to each of the inflaton fields can be estimated as $\Delta\phi_n = \Delta\phi/\sqrt{n}$, where $\Delta\phi$ denotes to the effective single-field inflaton range needed to generate the desired cosmological observables, and $n$ is the number of moduli assisting the multi-fibre inflation. We also present various numerical benchmark models consistently producing cosmological observables in the lights of the recent ACT experiments.

An independent indirect detection search for Dark Matter-Matter (DM-M) interactions is undertaken to set cross section limits based on neutron production data collected by the NMDS-II detector for 1440 hours at 1166 m.w.e. and 6504 hours at 583 m.w.e.. The detector system consists of a 30 cm cube Pb-target instrumented with 60 $^3$He neutron counters. The neutron detector system calibrated with a $^{252}$Cf source yields a single particle detection efficiency of 23.2\%$\pm$1.2\%. During data collection, the highest neutron multiplicity event observed 54 neutrons. The neutron multiplicity, n, distribution, fits well to a power law $k \times n^{-p}$, for both the data and cosmic ray muon induced neutron production in Geant4 simulations. Two DM-M interaction models were used to set limits. The first, a spallation model, assumes a single proton with kinetic energy equal to the DM-M interaction energy. The other, a fire-ball model assumes an annihilation between DM-M producing pions with a limiting Hagedorn temperature. The two extreme models produce similar upper DM-M cross section limits over the DM mass range between 300 MeV to 100 GeV. Limits assume all the DM energy is deposited in the Pb-target. Spin independent limits, proportional to A$^{-2}$, are at the level ~10$^{-45}~ cm^2$. Spin dependent limits, proportional to A$^{-1}$, are at the level, $2\times 10^{-42}~cm^2$.

This work introduces a Hamiltonian approach to regularization and linearization of central force particle dynamics through a new canonical extension of the so-called "projective decomposition". The regularization scheme is formulated within the framework of classic analytical Hamiltonian dynamics as a redundant-dimensional canonical/symplectic coordinate transformation, combined with an evolution parameter transformation, on extended phase space. By considering a generalized version of the standard projective decomposition, we obtain a family of such canonical transformations which differ at the momentum level. From this family of transformations, a preferred canonical coordinate set is chosen that possesses a simple and intuitive connection to the particle's local reference frame. Using this transformation, closed-form solutions are readily obtained for inverse square and inverse cubic radial forces (or any superposition thereof) on any finite-dimensional Euclidean space. From these solutions, a new set of orbit elements for Kepler-Coulomb dynamics is derived, along with their variational equations for arbitrary perturbations (singularity-free in all cases besides rectilinear motion). Governing equations are numerically validated for the classic two-body problem, incorporating the J_2 gravitational perturbation.

In this paper, we systematically investigate the optical appearance of a Schwarzschild black hole illuminated by three geometrically thin accretion disk models under varying observational inclination angles. Based on the geometric relationship between the black hole and observer, we first divide the accretion disk into co-side and counter-side semi-disks. Subsequently, we analyze their differences in photon orbit distributions, transfer functions, and received intensities. Our results reveal distinct inclination-dependent behaviors between the two semi-disks: when the inclination approaches $\pi/2$, the co-side semi-disk develops richer ring structures with the lensed emission dominating the total brightness. Furthermore, through explicit specification of the emission profiles of the three models, we present optical images for both optically thin and thick disk scenarios at different inclinations. The results demonstrate that: (i) the bright rings in all three models become progressively compressed and deviate from circularity as inclination increases; (ii) for thick disks, partial rings are obscured and the overall intensity is lower than thin disks. These results may advance our understanding of general black hole imaging processes and provide a new approach to test gravitational theories through optical morphology studies.

We show that a single strongly coupled gauge theory naturally gives rise to a unified description of the dark sector, simultaneously addressing the origin of dark energy at the meV scale and self-interacting dark matter at the TeV scale. Building on the density-responsive scalar field framework that resolves cosmological singularities and generates dark energy without fine-tuning, we show that the same underlying SU(3) gauge theory with $N_f = 10$ flavors predicts complete dark matter phenomenology. The theory generates not only the required anomalous dimension $\gamma \approx 0.50\pm0.05$ for the dark-energy mechanism, but also produces a composite particle spectrum featuring a TeV-scale dark matter candidate ($m_\chi \approx 600$~GeV) and a heavier scalar meson. Remarkably, the mass ratio $m_{\Phi_h}/m_\chi \approx 2$ naturally emerges from the dynamics of confinement, without any fine-tuning. This resonance condition, which initially appears as a coincidence, is revealed to be a robust prediction of the gauge dynamics, completely resolving the long-standing tension between achieving the correct relic abundance and solving small-scale structure problems through self-interactions. The framework makes sharp, testable predictions: a narrow scalar resonance at $\sim 1.2$~TeV accessible at HL-LHC, dark matter self-interactions of $\sigma_T/m \sim 1$~cm$^2$/g at dwarf galaxy velocities, and a stochastic gravitational wave background from the confinement phase transition. By deriving both dark energy and dark matter from a single microscopic theory, our work demonstrates how 95\% of the universe's energy content can emerge from one underlying strongly-coupled gauge theory, providing a complete and predictive framework for dark sector physics across 15 orders of magnitude in energy.

We introduce a novel class of algorithms, the ``spatially varying boost'', for generating dynamical field initial conditions with prescribed bulk velocities. Given (non-moving) initial field data, the algorithm generates new initial data with the given velocity profile by performing local Lorentz boosts. This algorithm is generic, with no restriction on the type of the field, the equation of motion, and can endow fields with ultra-relativistic velocities. This algorithm enables new simulations in different branches of physics, including cosmology and condensed matter physics. For demonstration, we used this algorithm to (1) boost two Sine-Gordon solitons to ultra-relativistic speeds for subsequent collision, (2) generate a relativistic transverse Proca field with random velocities, and (3) set up a spin-$1$ Schrödinger-Poisson field with velocity and density perturbations consistent with dark matter in matter dominated universe.

The Ukrainian satellite telescope of electrons and protons STEP-F conducted a scientific experiment in space on board the low Earth orbit spacecraft "CORONAS-Photon" in 2009. In the article it is described the background, stages of development, manufacturing and tests of weight-dimensional, breadboard, technological, flight and auxiliary models of the instrument, and of the control and test suites for the STEP-F. The principles of operation, design, and technical and scientific characteristics of the instrument have been described. Stages of adjustment, graduating and autonomous, docking, complex, ground pre-flight and flight tests have been shown. Scientific results obtained by the STEP-F instrument during the time of deep minimum of solar activity demonstrate the detection of new features in high-energy charged particle distributions inside radiation belts of the Earth, in the region of the Brazilian magnetic anomaly, and outside of the referred volumes. Keywords: Earth$'$s magnetosphere, satellite instrument, scintillation detector, silicon matrix, charged particles, radiation belt, solar activity, Brazilian magnetic anomaly.

The validity of inflation models is mainly evaluated according to the consistency of the predicted scalar spectral index $n_{\mathrm{s}}$, the tensor scalar ratio $r$, and the running scalar spectral index $\alpha_{\mathrm{s}}$ with cosmic microwave background observations. In warm inflation (WI) scenarios, one can find exact analytical solutions for $\alpha_{\mathrm{s}}$ in principle, but long expressions may be obtained. Previous studies for WI scenarios have only shown approximate analytical solutions or numerical results for $\alpha_{\mathrm{s}}$. In this study, we present a general analytical expression of $\alpha_{\mathrm{s}}$ without approximation in WI. The obtained analytical expression of $\alpha_{\mathrm{s}}$ is used in the study of warm natural inflation (WNI). Although $n_{\mathrm{s}}$ and $r$ have been previously investigated, $\alpha_{\mathrm{s}}$ is omitted in previous studies on WNI. Our study of $\alpha_{\mathrm{s}}$ completes previous phenomenological studies on WNI. In particular, the lower limit of the symmetry-breaking scale in WNI becomes more concrete in this study.

Pulsar timing array (PTA) collaborations recently presented evidence for a gravitational-wave background (GWB) signal at nanohertz frequencies. In this paper, we introduce new refitting techniques for PTA data analysis that elevate related techniques in the literature to a more rigorous level and thus provide the basis for fast and accurate Bayesian inference and physically intuitive model comparisons. The key idea behind our approach is to construct maps \Phi from GWB spectral models to a running-power-law (RPL) reference model, such that the pullback \Phi^* P_RPL of the RPL posterior density P_RPL induces a likelihood on the GWB model parameter space; in other words, we refit spectral models to the RPL posterior density. In order to construct \Phi, we introduce a matched-filtering approach in which \Phi follows from a \chi^2 minimization that accounts for the frequency dependence of PTA sensitivity curves. We validate and illustrate our techniques by three concrete examples: GWs from stable cosmic strings, GWs from metastable strings, and scalar-induced GWs.

We present a new method to search for gravitational waves from quasinormal modes in the ringdowns of the remnants of the mergers of the binary black hole systems. The method is based on maximum likelihood estimation. We derive a time-domain matched-filtering statistic that can be used to search for any number of modes in the data. The parameters of the modes can be estimated and the modes present in the data can be reconstructed. We perform Monte Carlo simulations of the method by injecting the quasinormal mode waveforms to noise. We analyze performance of the method for searches of quasinormal modes in the advanced detectors data like LIGO and Virgo, in the third generation of detectors like Einstein Telescope and Cosmic Explorer and in the space detector LISA data. We analyze ringdown of publicly available GW190521 event and we compare our results with analysis by other methods.

Ultra-supercooling phase transitions can generate large overdensities in the Universe, potentially leading to the formation of primordial black holes (PBHs), which can also be a dark matter candidate. In this work, we focus on the supercooling phase transition for the scale symmetry breaking based on the effective potential of the Coleman-Weinberg (CW) type. We investigate the effect on the PBH production in the presence of an additional mass term for the CW scalar field, what we call a soft-scale breaking term, which serves as the extra explicit-scale breaking term other than the quantum scale anomaly induced by the CW mechanism. We demonstrate that even a small size of the soft-scale breaking term can significantly affect the PBH production depending on its sign: a positive term slows down the phase transition, thereby enhancing the PBH abundance and improving the model's ability to account for dark matter; in contrast, a negative term suppresses the PBH formation. The inclusion of such soft-scale breaking terms broadens the viable parameter space and increases the flexibility of the framework. We further illustrate our results through two ultraviolet-complete realizations: i) a many-flavor QCD-inspired model as a reference model which can dynamically induce a positive-soft scale breaking; ii) a Higgs portal model with a $B-L$ scalar as the benchmark for the case where a negative-soft scale breaking is induced. Our study would provide a new testable link between PBH dark matter and gravitational wave signatures in the CW-type scenario.

We present a theoretical and experimental study of superconducting ring resonators as an initial step towards their application to superconducting electronics and quantum technologies. These devices have the potentially valuable property of supporting two orthogonal electromagnetic modes that couple to a common Cooper pair, quasiparticle, and phonon system. We present here a comprehensive theoretical and experimental analysis of the superconducting ring resonator system. We have developed superconducting ring resonator models that describe the key features of microwave behaviour to first order, providing insights into how transmission line inhomogeneities give rise to frequency splitting and mode rotation. Furthermore, we constructed signal flow graphs for a four-port ring resonator to numerically validate the behaviour predicted by our theoretical analysis. Superconducting ring resonators were fabricated in both coplanar waveguide and microstrip geometries using Al and Nb thin films. Microwave characterisation of these devices demonstrates close agreement with theoretical predictions. Our study reveals that frequency splitting and mode rotation are prevalent in ring systems with coupled degenerate modes, and these phenomena become distinctly resolved in high quality factor superconducting ring resonators.

The Dark Energy Survey (DES) and the Dark Energy Spectroscopic Instrument (DESI) measurements claim that the dark energy equation of state $w \ne -1$. This observation can be explained by the axion Dark Energy (aDE) model of an ultralight axion plus a cosmological constant $\Lambda$. Despite a relatively large degeneracy, there is a high probability that $\Lambda <0$. This negative $\Lambda$ leads the universe to end in a big crunch. Using the best-fit values of the model as a benchmark, we find the lifespan of our universe to be 33 billion years.

We propose a new axion dark matter detection strategy that employs optical readout of laser beam ellipticity modulations caused by axion-induced electric fields in a microwave cavity, using electro-optic (EO) crystals, enhanced by externally injected radio-frequency (RF) power. Building upon the variance-based probing method, PRD 107, 103005, (2023), we extend this concept to the optical domain: a weak probe laser interacts with an EO crystal placed inside a resonant microwave cavity at cryogenic temperatures, and the axion-induced electric field is revealed through induced ellipticity. The injected RF signal coherently interferes with that of the axion field, amplifying the optical response and significantly improving sensitivity. While our EO-based method employs a Fabry-Perot resonator, we do not require Michelson interferometers. Our method hence enables compact, high-frequency axion searches, across the 0.5-50\,GHz range. Operating at cryogenic temperatures not only suppresses thermal backgrounds but, critically, allows the probing method to mitigate the laser quantum noise. This approach offers a scalable path forward for axion detection over the $\sim (\text{few}-200$)~$\mu$eV mass range -- covering the preferred parameter space for the post-inflationary Peccei-Quinn axion dark matter -- using compact, tunable systems.