Papers for Monday, Feb 16 2026

Over the past few years, Astrophysics has experienced an unprecedented increase in research output, as is evident from the year-over-year increase in the number of research papers put onto the arXiv. As a result, keeping up with progress happening outside our respective sub-fields can be exhausting. While it is impossible to be informed on every single aspect of every sub-field, this paper aims to be the next best thing. We present a summary of statistics for every paper uploaded onto the Astrophysics arXiv over the past year - 2025. We analyse a host of metadata ranging from simple metrics like the number of pages and the most used keywords, as well as deeper, more interesting statistics like the distribution of journals to which papers are submitted, the most used telescopes, the most studied astrophysical objects including GW, GRB, FRB events, exoplanets and much more. We also indexed the authors' affiliations to put into context the global distribution of research and collaboration. Combining this data with the citation information of each paper allows us to understand how influential different papers have been on the progress of the field this year. Overall, these statistics highlight the general current state of the field, the hot topics people are working on and the different research communities across the globe and how they function. We also delve into the costs involved in publications and what it means for the community. We hope that this is helpful for both students and professionals alike to adapt their current trajectories to better benefit the field.

We present GAMERA-OP (Orthogonal-Plus), a three-dimensional finite-volume magnetohydrodynamics (MHD) solver for orthogonal curvilinear geometries. The solver advances magnetic fields using constrained transport to preserve $\nabla\!\cdot\!\mathbf{B}=0$ to machine precision and employs geometry-consistent high-order reconstruction with an enhanced Partial Donor Cell method (e-PDM) that accounts for geometry curvature. Flexible numerics include various numerical fluxes and time integrators. In axial symmetric coordinates, angular momentum are preserved to round-off, and a ring-averaging treatment near the axis relaxes CFL constraints while maintaining divergence-free magnetic fields. Optional capabilities include the semi-relativistic (Boris) correction, background-field splitting, and an anisotropic MHD formulation. Rewritten in C, the code adopts a modular design that simplifies case setup and facilitates the addition of physics modules and coupling to other first-principles codes. Standard benchmarks across multiple geometries verify the code's high accuracy, low numerical diffusion, and robust handling of coordinate singularities and rotating flows. GAMERA-OP provides a practical, high-order framework for space and astrophysical plasma applications where orthogonal curvilinear coordinates and exact angular-momentum conservation are advantageous.

We investigate the physical origin of star formation suppression in gas-rich early-type galaxies using five high-resolution hydrodynamical idealized galaxy simulations, performed with the moving-mesh code AREPO. These simulations include one Milky Way-like galaxy and four early-type galaxies, of which one early-type galaxy is found to have significantly less star formation despite a substantial molecular gas reservoir. We apply a modified virial theorem to the overdensities in each galaxy to quantify the forces regulating their stability and thus star formation. We find evidence that, in the suppressed galaxy, strong Coriolis forces driven by elevated galactic shear may inhibit gravitational collapse. This is caused by the galaxy's high central compactness, providing a physical mechanism for the suppression of star formation that does not require the removal of molecular gas. In contrast, less compact early-type galaxies host more gravity-dominated clouds and therefore exhibit higher star formation rates. However, we find that this gravitational stability occurs without significantly increasing the classical Toomre-Q parameter, and therefore a new criterion for suppressed star formation may be needed. We also discuss the impact of our choice of overdensity scale and connections to observations of molecular clouds.

We present a panchromatic optical-mm characterization of GNz7q, a recently identified X-ray weak, rapidly growing red quasar embedded within a dusty starburst galaxy at $z=7.1899$, using the full suite of JWST/NIRCam, NIRSpec, MIRI, and archival NOEMA observations. Our deep NIRSpec/G395M spectroscopy reveals unambiguous broad Balmer emission (FWHM $=2221\pm20$kms$^{-1}$), confirming a super-Eddington accreting black hole ($\lambda_{\rm Edd}=2.7\pm0.4$) with a mass of $\log(M_{\rm BH}/M_{\odot})=7.55\pm0.34$, using accretion-rate corrected BH mass estimators. After subtracting the point source, we robustly detect stellar emission from the host galaxy across multiple NIRCam and MIRI filters. Out joint morphological-spectral analysis yields a stellar mass of $\log (M_*/M_\odot)=10.5\pm0.4$ and an intense star formation rate of ${\rm SFR}=330\pm97\,M_\odot\,\rm yr^{-1}$, confirming the host as a massive, dusty starburst galaxy. We find that GNz7q lies on the local $M_{\rm BH}$-$M_*$ relation ($M_{\rm BH}/M_*\simeq 0.001$) and is well positioned to evolve into the locus of massive SDSS quasars with $\log (M_{\rm BH}/M_\odot)\approx 9$ and $M_*\approx 10^{11}\,M_\odot$ at $z\sim 6$, owing to its remarkably rapid growth in both the black hole and its host galaxy. This stands in stark contrast to many recently reported JWST AGN populations at similar redshifts, including the little red dots (LRDs), whose weak or undetected star formation makes it difficult for them to grow into the massive galaxies hosting SDSS-like quasars. These results suggest that GNz7q marks as a rare, pivotal phase of early BH-galaxy co-eolution, plausibly providing a crucial direct pathway to the supermassive black hole systems within the first billion years of the Universe.

The catalog of Gravitational Wave (GW) events is rapidly growing, providing key insights into the evolution of massive binaries and compact object formation. However, a key challenge is to explain the origin of exceptional events such as GW190814, among the most asymmetric mass-ratio mergers to date ($q\approx 0.1$). We show that it shares an evolutionary pathway with the most unequal mass Galactic High Mass X-ray Binary (HMXB) 4U 1700-37/ HD 153919. We demonstrate this unique connection by utilizing a rich set of existing observational constraints for the HMXB and compute detailed binary evolution models to explain its formation history. We find that conservative mass transfer, along with a directed natal kick are essential to explain its current state. We show that this system is unlikely to form a GW source due to a failed Common Envelope (CE) phase in the future, in agreement with previous work. With additional models, we show that a similar pathway naturally forms GW190814-like events, provided the first phase of mass transfer remains conservative, and the first-born (lower mass) compact object receives a large natal kick ($\gtrsim 100\,\mathrm{km/s}$) for the subsequent CE phase to be successful and form a asymmetric mass-ratio GW source. Anchored by the number of analogous Galactic HMXBs, we estimate rates for such GW events, which broadly agree with their observed rate. Our work demonstrates a unified formation pathway for highly asymmetric mass-ratio HMXBs and GW events. Moreover, it highlights the critical role of finding and characterizing local analogs in different evolutionary phases, and using them as a bridge to understand the origin of GW sources, especially the outliers like GW190814.

Combining deep Hubble Space Telescope (HST) images and extensive data from the Multi-Unit Spectroscopic Explorer, we present new mass models of the cluster MACS J1149.5+2223, strongly lensing the supernova (SN) Refsdal, fully exploiting the source surface-brightness distribution of the SN host for the first time. In detail, we incorporated 77,000 HST pixels, in addition to the known 106 point-like multiple images, in our modeling. We considered four different models to explore the effect of the relative weighting of the point-like multiple image positions and flux distribution of the SN host on the model optimization. When the SN host's extended image is included, we find that the statistical uncertainties of all 34 free model parameters are reduced by factors ranging from one to two orders of magnitude compared to the statistical uncertainty of the point-like only model, irrespective of the adopted different image weights. We quantified the remarkably increased level of precision with which the cluster's total mass and the predicted time delays of the SN Refsdal multiple image positions can be reconstructed. We also show the delensed image of the SN host, a spiral galaxy at zSN = 1.49, in multiple HST bands. In all those applications, we obtain a significant reduction of the statistical uncertainty, which is now below the level of even the small systematic uncertainty on the mass model that could be assessed by the different approaches. These results demonstrate that with extended image models of lensing clusters it is possible to measure the cluster's total mass distribution, the values of the cosmological parameters, and the physical properties of high-redshift sources with an unparalleled precision, making the typically not-quantified systematic uncertainties now crucial.

Polycyclic aromatic hydrocarbons (PAHs) are responsible for a variety of near- and mid-infrared spectral features in Galactic and extragalactic sources. A feature at 1.05 $\mu$m arising from electronic transitions in PAH cations is predicted by laboratory experiments but has never been observationally confirmed. We conduct a dedicated search for this feature in absorption on a highly-extinguished sight line toward BD+40 4223, a blue supergiant in Cyg OB2, using the TripleSpec spectrograph at Palomar Observatory. We place a $5\sigma$ upper limit on the feature strength of $\Delta\tau_{1.05}/A_{V} < 5.6 \times10^{-3}$, ruling out theoretical estimates with $> 10\sigma$ significance. We constrain the effective temperature of BD+40 4223 to be $\log_{10}\left(T_{\rm eff}\right)=4.41\pm0.03$ and infer that it is veiled by $6.39\pm0.05$ magnitudes of visual extinction, consistent with but more constraining than previous determinations. As dust on the sight line toward BD+40 4223 appears typical of the diffuse interstellar medium, this non-detection challenges existing models of PAH material properties and/or charge distribution.

Based on deep high resolution JWST images of AS1063, and after a careful masking of artifacts, extended features in the cluster, and background galaxies (including known lensed ones), we have identified tens of thousands of unresolved point sources in the central region of the galaxy cluster. We extended the identification of these point sources up to 1.18 Mpc from the center of the cluster using data in the second module. Most of these sources are expected to be globular clusters orbiting in the deep potential well of the cluster, but also the surviving compact cores of satellite galaxies. We study the distribution of the globular clusters and compared it with the distribution of mass from a lens model derived from the same JWST data. We find a very tight correlation between the two distributions, but also some differences, including a more concentrated distribution for the globular clusters than for dark matter. We explored the possibility of using the distribution of globular clusters as a proxy for the lensing mass. We find that a simple smoothing kernel can transform the discrete distribution of point sources into a continuous two-dimensional distribution that matches well the lensing convergence. This suggests that globular clusters can be used as tracers of the dark matter distribution in other massive clusters where gravitational lensing constraints are scarce but globular clusters can be detected more easily, for instance in low redshift galaxy clusters.

We present a simple analytic approximation for the covariance between pre-reconstruction galaxy power spectrum measurements and post-reconstruction two-point correlation functions. This cross-covariance is essential for joint analyses that combine full-shape clustering information with baryon acoustic oscillation (BAO) measurements, as commonly performed in modern spectroscopic surveys. Our model builds on the disconnected contribution to the covariance and accounts for the damping of correlations due to the BAO reconstruction process. We validate our analytic prescription against numerical simulations from the Dark Energy Spectroscopic Instrument (DESI), testing both idealized cubic geometries and realistic survey configurations including complex footprints and fiber assignment effects. Despite neglecting survey window functions in the analytic calculation, we find excellent agreement with simulation-based covariances and demonstrate that cosmological parameter constraints are virtually unchanged when using our approximation. Our results show that the pre-post cross-covariance is sufficiently small that even approximate treatments are adequate for cosmological inference, opening a pathway toward fully analytic covariance matrices for next-generation galaxy surveys.

Sign-switching dark energy provides a novel mechanism for modifying the late-time expansion history of the Universe without invoking additional fields or finely tuned initial conditions. In this work, we investigate a class of background--level cosmological models in which the dark energy contribution changes sign at a transition redshift $z_\dagger$, producing a sharp deviation from standard $\Lambda$CDM dynamics. We confront these models with a comprehensive set of cosmological observations, including \textit{Planck 18} cosmic microwave background (CMB) measurements, \textit{\textit{DESI DR2}} Baryonic Acoustic Oscillation (BAO) data and the \textit{Pantheon+} $\&$ \textit{SH0ES} Type Ia supernova sample (SN). Using a full Markov Chain Monte Carlo (MCMC) analysis, we find that the sign-switching scenario significantly alleviates the Hubble tension while obtaining better results when statistically comparing with $\Lambda$CDM, as quantified by the Akaike and Bayesian information Criteria. Although the model is explored only at the background level, the improvement in the inferred Hubble constant demonstrates that sign-switching dark energy offers a promising and physically economical pathway toward resolving late-universe discrepancies.

The excitation and damping mechanisms for oscillation modes of gas giant planets are undetermined. We show that differential rotation may greatly enhance convective viscosity in giant planets, resulting in damping times of $t_{\rm damp} \sim 10^5-10^6 \, {\rm years}$ for f modes and low-order p modes. Radiative diffusion damps p modes on time scales of $t_{\rm damp} \sim 10^3-10^7 \, {\rm years}$. While the lethargic convective motions cannot effectively excite f mode or p modes, storms driven by condensation of water and/or silicates may play a role. High-order p modes are most effectively excited by cometary/asteroid impacts. Applying these calculations to solar system planets, water storms, rock storms, and impacts may all contribute to exciting the observed f modes amplitudes of Saturn via ring seismology. Similar f mode amplitudes with fractional gravitational perturbations of $\delta \Phi/\Phi \sim 10^{-10}-10^{-9}$ are expected for Jupiter and Uranus, apart from their lowest $\ell$ f modes which could have larger gravitational perturbations of $\delta \Phi/\Phi \sim 10^{-7}$. Rock storms may contribute to mode driving in Jupiter, while water storms are more important for Uranus. The highest-amplitude p modes are predicted to have periods of $\sim$10-30 minutes, with surface velocities of $\sim$10 cm/s for Jupiter and Saturn, and $\sim$1 cm/s for Uranus. These oscillation modes may be detectable with radial velocity measurements, ring seismology, or spacecraft Doppler tracking. However, both the damping and excitation physics are uncertain by orders of magnitude, so more careful examination of the relevant physics is required for robust estimates.

We present our analysis of how the detectability of carbon dioxide (CO2) on an Earth-like planet varies with respect to signal-to-noise ratio (SNR), wavelength, and molecular abundance. Using the Bayesian Analysis for Remote Biosignature Identification on exoEarths (BARBIE) methodology, we can inform the optimal long-wavelength cut-off for the future Habitable Worlds Observatory (HWO) coronagraph. We test 25 evenly-spaced 20% bandpasses between 0.8-2.0{\mu}m, and simulate data spanning a range of SNRs and molecular abundance to analyze the relationship between wavelength and detectability for different planetary archetypes. We examine abundance levels from varying Earth epochs and a Venus-like archetype to investigate how detectability would change throughout the evolution of a rocky planet. Here, we present our results on the planetary conditions and technological requirements to strongly detect CO2. In addition, we analyze the degeneracy of CO2 with carbon monoxide (CO), methane (CH4), and water (H2O). We determine that any abundance of CO does not achieve strong detections and that CH4 and H2O play a pivotal role in the ability to detect CO2. We conclude that the optimal long-wavelength cut-off for the Habitable Worlds Observatory coronagraph should be 1.68{\mu}m.

Context. Tidal disruption events (TDEs) are rare transients that provide important insights into the physics of galactic nuclei. A recently identified feature in their optical light curves is the presence of early bump-like structures (precursors) that appear before the onset of the main flare or during its rise. Aims. We aim to build and study the first sample of precursor TDEs in order to improve our understanding of these features, which could be key to revealing the origin of the optical emission in TDEs. Methods. We compiled all known precursor TDEs from the literature, searched for additional candidates, and analyzed them as a sample. Results. We find that precursor TDEs predominantly fall within the repeating TDE, fast TDE, and TDE in active galactic nucleus (AGN) subclasses. We reveal a positive correlation between the occurrence time of the precursors relative to the main peak and the central black hole mass. Conclusions. We suggest that the precursors appear due to interactions between the incoming stellar debris and the disk or leftover material from an earlier disruption (repeating and fast TDEs) or a stable pre-existing disk (TDEs in AGNs). Precursors are therefore potentially key signatures of repeating partial TDEs in previously quiescent galaxies.

The analysis of the large-scale clustering of quasars (QSO) observed by the Dark Energy Spectroscopic Instrument (DESI) represents a promising avenue for constraining local Primordial non-Gaussianity (PNG), parameterized by $f_{\rm NL}$. The signal to be constrained is the scale-dependent bias induced in the 2-point clustering of the considered tracer sample. The resulting constraints on $f_{\rm NL}$, however, are fully degenerate with the local PNG bias parameter $b_{\phi}$, dependent on the assembly bias parameter $p$. Using IllustrisTNG hydrodynamical simulations, we select a QSO sample reflecting the selection criteria and properties of DESI QSOs, and provide a robust prior for $p$, and thus for $b_{\phi}$, building on the findings of Fondi et al. 2024. We find a distribution with mean $\bar{p}\simeq1.4$ with weak redshift dependence, stable to selection noise and consistent with the expected recent merger history typical of quasar-hosting halos. By comparing with the CAMELS simulations we demonstrate that this prior is robust to astrophysical assumptions and cosmic variance. Finally, applying this prior to the DESI DR1 dataset, we derive updated constraints on local PNG, obtaining $f_{\rm NL}=-3.3\pm9.2$.

[Abridged] We investigate how the identification of Hickson-like CGs depends on the observer's reference frame, quantifying how frequently the same system would be recognised from different vantage points. Using a mock lightcone built from the Millennium I Simulation plus a semi-analytic model of galaxy formation, we identified 7709 CGs when applying the standard Hickson-like criteria. For each CG, we placed 1000 random observers on a surrounding sphere and reapplied the velocity and compactness requirements to test recoverability. We also examined the variation of population and local isolation. The velocity concordance criterion shows modest sensitivity to the observer's location: 10% of CGs fail for some observers, typically groups with members with high peculiar velocities (>1000 km/s). The compactness requirement is far more fragile, as 44% of CGs are missed by most observers, and these systems are very elongated or are chance alignments in real space. Tightening selection limits reduces this dependence. Lowering the surface brightness threshold to $\mu \leq 23 \ mag/arcsec^2$ reduces the compactness dependence to 16%, while reducing the velocity limit to $\Delta V\leq 250 \ km/s$ lowers velocity-driven failures to less than 4%. Applying both cuts simultaneously yields up to 84% observer-independent groups, although with a substantially smaller sample. Population and isolation are affected by bright interlopers seen from different directions. While such interlopers are common, they have only a minor effect on the compactness and velocity concordance criteria; however, the local isolation is commonly broken. Observer frame effects, dominated by the compactness criterion, can significantly bias Hickson-like CG samples. However, adjusting surface brightness and velocity difference thresholds allows users to balance the physical reliability according to their specific scientific goals.

3I/ATLAS is the third known interstellar object to pass through our Solar System. NASA's Transiting Exoplanet Survey Satellite (TESS) made dedicated observations of 3I/ATLAS between 15 -- 22 January 2026 (Sector 1751), capturing high-cadence observations at 200s and 20s cadence. We present two High Level Science Products (HLSPs): (1) comet-centered image time series, corrected for background scattered light and stars; and (2) aperture light curves extracted from the corrected images. We created these data products using the official TESS products and they are publicly available at the Mikulski Archive for Space Telescopes (MAST). TESS's high-precision, near-continuous photometry will provide unique insights into the comet's activity following its closest approach to the Sun. The TESS Science Support Center (TSSC) has created these data products to facilitate scientific analyses by the TESS and Solar System communities.

We present new time-dependent inverse sensitivities for the WFC3/IR channel. These were calculated using the sensitivity change slopes measured by \citet{2024wfc..rept....6M} and photometry of five CALSPEC standards (the white dwarfs GRW+70~5824, GD~153, GD~71, G191B2B, and the G-type star P330E) collected from 2009 to 2023. The new inverse sensitivities account for losses of 1-2\% over 15 years, depending on wavelength, and provide an internal photometric precision better than 0.5\% for all wide--, medium--, and narrow-band filters. An updated version of \texttt{calwf3} (v3.7.2) has been developed for use with a new time-dependent image photometry table (IMPHTTAB) and will be used to update the image header photometric keywords following MAST reprocessing, expected in late-2024. Alternatively, the new inverse sensitivities may be computed by the user for a specific observation date by running \texttt{stsynphot}.

X-ray pulse-profile modeling of millisecond pulsars offers a direct route to measuring neutron star masses and radii, thereby constraining the dense-matter equation of state. However, standard analyses typically rely on \emph{ad hoc} hotspot parameterizations rather than self-consistent physical models. While connecting surface heating directly to the magnetospheric geometry provides a more natural physical pathway, computing global magnetospheric solutions is too computationally expensive to perform on-the-fly during parameter inference. In this work, we bridge this gap by deriving fully analytic, first-principles expressions for surface return currents in mixed dipole--quadrupole magnetospheres. Working within force-free electrodynamics, we generalize the field-aligned current invariant $\Lambda$, the crucial scalar that maps the far-zone magnetic structure to the near-zone heating rate, from the standard dipole approximation to arbitrary quadrupolar configurations. We demonstrate that even when the quadrupole component is sub-dominant in the far zone (the mixing regime), using a dipole-based heating prescription fails to capture the significant enhancement or suppression of the return-current density on the polar cap. Our consistent quadrupole-aware framework reveals that these multipolar currents redistribute the surface heating, leading to systematic discrepancies in predicted pulse profiles that are amplified by atmosphere beaming and can reach $\sim 30\%$ near pulse peaks. These results provide a rigorous analytic foundation for mapping global magnetic geometry to surface heating in multipolar magnetospheres, enabling physically consistent inference beyond the idealized dipole approximation.

Magnetic reconnection is one of the prime candidate mechanisms that may energize the plasma emitting the strongly polarized X-ray emission from black hole X-ray binaries (BHXRBs) in their hard states. The mechanism requires strong magnetic fields in the upstream plasma entering the reconnection layer, and weaker, but still substantial, magnetic fields in the downstream regions. In this Letter, we estimate the coronal magnetic fields for three different magnetic energy dissipation mechanisms: plasmoid-dominated magnetic reconnection, fast collisionless reconnection, and magnetic field relaxation. We show that the lack of strong Faraday depolarization constrains viable models and can be used to benchmark numerical accretion flow models. We conclude by discussing the difficulties of disentangling the various effects that can depolarize the signals from BHXRBs at low energies. We furthermore emphasize that Faraday rotation is unlikely to play a role in the polarization of the coronal X-ray emission of active galactic nuclei.

The inner Galaxy is the most complex region of the Milky Way, comprising the bulge, inner thin and thick discs, and inner halo; the formation of the bar transferred gas and stars from the disc inward. Accretion of dwarf galaxies also occurred over the Galaxy's lifetime, merging with the original bulge. In this work, we constrain the metal-rich stars of the earliest spheroidal bulge. To study the oldest bulge stars, distributed in a spheroid, we applied kinematical and dynamical criteria in the metal-rich range [Fe/H] > -0.8. This complements our previous analysis of a symmetric sample with [Fe/H] < -0.8. We derived individual abundances through spectral synthesis for C, N, O, Al, P, S, K, Mn, and Ce using stellar parameters from APOGEE DR17, and compared the results with literature data and chemical-evolution models. The alpha elements Mg, Si, and Ca, and iron-peak elements V, Cr, Co, and Ni follow the expected trends relative to the models. Mn shows secondary behaviour. S and K display significant star-to-star scatter but remain broadly compatible with predictions. Phosphorus and cerium show an excess around [Fe/H] $\sim$ -0.7, more pronounced than in the metal-poor sample, suggesting a distinctive signature of the earliest bulge population. Diagrams of [Mg/Mn] versus [Al/Fe] and [Ni/Fe] versus [(C+N)/O] indicate an in situ origin for most stars. At super-solar metallicities, a subset shows enhanced K and Mn (possibly S) with low [Ce/Fe], hinting at enrichment linked to the nuclear disc and bar, and tracing a chemically distinct population shaped by the innermost Galaxy.

We report on our investigation of the NuSTAR and AstroSat observations along with simultaneous NICER observations of the accreting millisecond X-ray pulsar SAX J1808.4$-$3658, obtained during its tenth outburst from 2022. The NuSTAR observation captured the source near the outburst peak, while AstroSat observed it during the decay phase. Coherent pulsations at $\sim$401 Hz were detected throughout the outburst, with the fundamental amplitude in the 3--30 keV range increasing from $\sim$4% near the peak to $\sim$6% during the decay. The pulsations display strong energy dependence and negative time lags of $\sim$0.2--0.3 ms, with harder photons leading softer ones. The broadband spectra in both epochs are well described by a soft thermal component and Comptonized continuum, together with a prominent relativistic reflection component. As the outburst evolved, the continuum softened ($\Gamma$ increasing from $\sim$1.88 to $\sim$1.99) and the coronal electron temperature decreased ($kT_{\rm e}$ from $\sim$31 to $\sim$18 keV), consistent with enhanced Compton cooling at lower accretion rates. The ionization parameter declined ($\log \xi$ from $\sim$3.4 to $\sim$1.8) while the reflection fraction increased, suggesting a changing accretion geometry with a more compact corona and a larger disk covering fraction during the decay phase. The X-ray luminosity decreased by a factor of $\sim$3 between the two epochs. Our results suggest the coupled evolution of the corona, disk, and magnetosphere as the mass accretion rate declines.

Reliable predictions of weak rapid neutron capture ($\textit{r}$-process) abundances require a systematic treatment of nuclear physics uncertainties, especially neutron capture rates far from stability. We employ new neutron capture rates from cross sections calculated with Yet Another Hauser-Feshbach Code ($\texttt{YAHFC}$) using an uncertainty-quantified Koning-Delaroche potential modified for use on neutron-rich systems. Using these rates as a baseline, we perform Monte-Carlo studies with independent rate variations (uncorrelated Monte Carlo) and find correlations between specific neutron capture rates and the resulting elemental abundances for the three weak r-process scenarios: two separate simulations of neutron star merger remnant accretion disks and a simulation of a magnetorotational supernova. We discuss the underlying nuclear dynamics that give rise to these correlations and the role of astrophysical conditions in them. We demonstrate how reducing the uncertainty in these rates would improve the prospects for conducting precision $\textit{r}$-process studies in the future. We additionally present a correlated Monte Carlo study, which incorporates the full covariance matrix that describes the relationships between individual neutron capture rates that arise from an uncertainty-quantified optical potential. We find that the magnitude of the uncertainty in the abundance pattern is similar to that produced by an uncorrelated Monte Carlo that employs only the on-diagonal components of the covariance matrix. We show how correlations restructure how the abundances co-vary, but do not necessarily decrease the overall uncertainty envelope.

We present a new Pop III + Pop II subgrid framework implemented in the moving-mesh code {\sc arepo}, designed to study the impact of Pop III feedback on star formation in the early universe. The framework combines primordial non-equilibrium chemistry, metal-line cooling, IMF-sampled stellar evolution with SN feedback, and approximate Lyman-Werner (LW) and ionizing radiation transport. We run a suite of $1c{\rm Mpc}/h$ box simulations with different initial conditions and resolutions from $z=127$ to $z=10$. The highest gas mass and spatial resolution in the fiducial simulation reach $\sim10\,{\rm M_{\odot}}$ and $\sim4\,{\rm pc}$, respectively. The model successfully reproduces the UV-inferred Pop II star formation rate density (SFRD) from recent JWST observations across all initial conditions, with only minor variation driven by local halo interactions and LW irradiation. We find that the volume filling factor of metal-enriched gas converges to $\sim1\%$ at $z=10$. Convergence is achieved once subhalos with $M_{\rm subhalo}\gtrsim 10^{6.5}\,{\rm M_{\odot}}$ are resolved, and the total stellar mass at $z=10$ is largely insensitive to initial conditions or the resolution considered in this work. A fiducial simulation requires $\sim 10^4$ CPU hours, making the framework computationally tractable for larger box simulations and enabling future large parameter studies of stellar physics or environment effects such as Pop III IMF variations, X-ray radiation, or the streaming velocity at high redshift.

The physical nature and origins of the newly discovered class of Odd Radio Circles (ORCs) remain unclear. We investigate a model whereby ORCs are synchrotron-emitting vortex rings formed by the Richtmyer-Meshkov instability (RMI) when a shock interacts with a low-density fossil radio lobe in the intergalactic medium using 3D magnetohydrodynamic simulations. These rings initially exhibit oscillatory behavior that damps over time. We implement a new method to model Inverse-Compton cooling and synchrotron cooling at high frequencies in a scale-free manner, enabling us to test a wide range of model parameters against the observational constraints. We find that shock strengths of Mach 2-4 are consistent with the data, as expected in accretion, merger-driven, or active galactic nuclei-driven shocks. We find that the initial size of the bubbles required to explain the rings ranges from 140 to 250 kpc, with initial energy in the bubble of order $10^{57}-10^{59}$ erg, consistent with fossil lobes inflated by moderately powerful radio galaxies. Derived ambient pressures and densities place ORCs in low density environments, such as the outskirts of galaxy groups with ages of order 70-200 Myr. Our synthetic radio maps match the polarization properties of ORC1 and predict a dependency of the tangential magnetic field angle on the aspect ratio of ORCs. A key distinguishing trait of the RMI-driven vortex ring model is that it does not require the ORC to be centered on its host galaxy and is therefore redshift agnostic.

Ground Level Enhancement (GLE) in large solar energetic particle (SEP) events is indicative of protons accelerated to GeV energies. Almost all GLE events are associated with sustained gamma-ray emission (SGRE) from the Sun because the latter require >300 MeV protons that are readily present during GLEs. Here we report on the 2024 June 8 GLE event, which has the distinction of not being associated with an SGRE event. All the associated phenomena typical of SGRE events were present: a fast and wide CME, a major solar flare, and an intense type II radio bursts that extend from the metric to kilometric wavelength domains. There was a data gap of ~51 min, but the SGRE is expected to last for hours. We suggest the east-west asymmetry in the flow of energetic particles from the shock is likely to be the reason for the lack of SGRE emission.

Traditional interstellar polycyclic aromatic hydrocarbon (PAH) diagnostics rely on empirical band ratios, which often suffer from information loss and sample-selection bias. We introduce a machine learning framework that bypasses these limitations by treating the complete 2.75-20 micron emission spectrum as a high-dimensional fingerprint. Using a Random Forest classifier trained on a dataset of 23,653 spectra, we achieve a robust classification F1-score of about 0.96 across 12 size and charge categories. Our model maintains high performance on synthetic mixtures of unseen molecules. Feature importance analysis reveals that PAH size diagnostics are not universal but highly charge-dependent; while neutral size is traced mainly by C-H stretching modes, sizing ionized species also relies on the morphology of 6-8 micron C-C complexes, with the 12.5 micron feature emerging as a robust cross-charge tracer. This approach provides a robust, data-driven pathway for decoding the physical conditions of the interstellar medium.

We present Guangqi, a new two-dimensional, finite-volume radiation hydrodynamics code designed for high-performance astrophysical simulations. The code simultaneously resolves the hydrodynamic equations for complex equations of state (EoS) and implicit radiation transport under the flux-limited diffusion approximation. Written in Fortran and parallelized via the Message Passing Interface. Guangqi supports analytic hydrogen and helium EoS under the assumption of local thermal and chemical equilibrium. The framework is compatible with both Cartesian and spherical-polar geometries -- utilizing non-uniform grid spacing -- and incorporates static (SMR) and adaptive mesh refinement to optimize computational efficiency. To address the inherent challenges of angular momentum conservation in spherical-polar coordinates, we implement a robust and consistent "passive scalar angular momentum algorithm" (PSAMA). Domain decomposition is managed through both Z-order and Hilbert space-filling curves to ensure scalability. The code has been rigorously verified against a suite of standard benchmarks and newly designed test cases specifically intended to diagnose the non-linear coupling between gas dynamics, intricate EoS, radiation transport, and angular momentum conservation.

We present an analysis of medium/high-resolution JWST/NIRSpec spectra for 11 LRDs at $z \sim 4$. By decomposing the broad and narrow components of the Balmer emission lines, we investigate the connection between line emission and UV/optical continua for the LRD population. We find that the broad H$\alpha$ luminosity strongly correlates with the optical continuum (but not with the UV), indicating a common AGN origin for both. In contrast, the [O III] line strength is correlated with the UV continuum rather than the optical. Using the width and luminosity of the broad H$\alpha$ line, we estimate central black hole masses of $10^6-10^8 M_{\odot}$ accreting at high Eddington ratios, consistent with an early ($\lambda_{\rm Edd} \sim 0.6$), rapid-growth phase of AGN evolution. Assuming a constant mass accretion rate in the framework of slim-disk models, we infer growth timescales of $\sim 10^5-10^7\rm yr$, and suggest LRDs may evolve into narrow-line Seyfert 1 galaxies. Upper limits from our spectra indicate that LRDs exhibit intrinsically weak optical Fe II emission compared to typical AGN. To simultaneously account for the inferred broad-line region size and observed luminosity, we propose a "Clumpy Envelope" model in which the optical emission arises from an extended, clumpy gas with a characteristic radius of tens of light-days. The diversity in observed optical continuum shapes can be explained by radial temperature gradients and self-absorption effects within this structure. Our results demonstrate the power of JWST high-resolution spectroscopy in probing the central engines and physical nature of the LRD population.

JAXA's Hayabusa2 mission successfully returned samples from the asteroid Ryugu in December 2020. It executed two touchdowns to collect the surface and subsurface materials, one close to the crater created by an artificial impactor. The onboard camera system, Optical Navigation Camera (ONC), with two wide-angle cameras and one narrow-angle camera with seven color filters, was crucial for mapping geomorphology and composition such as hydrated minerals during navigation and scientific observation. More than 8,300 images revealed Ryugu's spinning-top shape and boulder-covered surface. However, most high-resolution images captured during descent/touchdown operations lacked precise location data and camera position/orientation information. Image geometry was refined using photogrammetric bundle adjustment. This method enabled the refinement of all high-resolution images captured during descent/touchdown operations. Furthermore, map-projected GeoTIFF images in GIS format containing geographic metadata were created for all ONC images, and these were integrated to construct global and regional mosaic maps. To facilitate scientific research on Ryugu, these refined image geometry information, maps, and mosaics are publicly available via this https URL

Polycyclic aromatic hydrocarbons (PAHs) are recognized as the primary contributors to the aromatic infrared bands (AIBs) widely observed in space. However, analyzing these AIBs remains challenging because of the immense structural diversity within the PAH family, which makes the computation of reliable reference spectra difficult. To address this, we developed an efficient graph neural network (GNN) framework that can predict PAH absorption spectra up to 10,000 times faster than traditional quantum chemical methods. We evaluated four representative GNN architectures, including graph convolutional network (GCN), graph attention network (GAT), message passing neural network (MPNN), and attentive fingerprint (AFP). The AFP model is found to deliver the best overall performance and is further trained using five different spectral distance metrics as loss functions, among which the Jensen-Shannon divergence yields the most accurate and stable results. The model performs best for PAHs containing 20-40 carbon atoms, while accuracy decreases for larger molecules, reflecting the limited availability of training data. Overall, this framework offers a fast method to generate approximate reference spectra for small- to medium-sized PAHs, supporting future AIB analysis.

We report the discovery and characterization of a planetary companion in the microlensing event KMT-2016-BLG-1337, which was produced by a binary system of low-mass stars. The light curve of the event exhibits a short-term anomaly superposed on the profile of a binary-lens single-source (2L1S) model. To investigate the nature of this anomaly, we performed detailed modeling under both the binary-lens binary-source (2L2S) and triple-lens single-source (3L1S) interpretations. The 3L1S model provides a substantially better fit to the data, strongly favoring the presence of a planetary companion in the lens system. Two viable $3L1S$ solutions describe the event nearly equally well. In one solution, the planet has a mass of $M_3 \sim 0.3~M_{\mathrm{J}}$ and lies at a projected separation of $a_{\perp,3} \sim 4~{\rm au}$ from the heavier member of the host binary. In the alternative solution, the planet has a mass of $M_3 \sim 7~M_{\mathrm{J}}$ and a projected separation of $a_{\perp,3} \sim 1.5~{\rm au}$. The host binary consists of early M-type dwarfs with masses of $M_1 \sim 0.54~M_\odot$ and $M_2 \sim 0.40~M_\odot$, separated in projection by $a_{\perp,2} \sim 3.5~{\rm au}$. The system is located at a distance of $D_{\rm L} \sim 7~{\rm kpc}$ toward the Galactic bulge. This event demonstrates the sensitivity of microlensing to planets in dynamically complex stellar environments, including systems beyond the reach of other detection techniques. It thereby contributes to a more comprehensive understanding of planet formation in multiple-star systems.

Radio relics in merging galaxy clusters are widely interpreted as synchrotron emission from relativistic electrons accelerated at large-scale shocks. However, the efficiency of diffusive shock acceleration (DSA) is expected to be reduced in the low-Mach-number, weakly turbulent environments characteristic of cluster merger shocks, and recent results suggest that DSA itself may not constitute a viable physical mechanism. In this work, we investigate ballistic surfing acceleration (BSA) as an electrodynamically grounded mechanism for electron energization that does not rely on prescribed diffusion coefficients. We formulate BSA under typical cluster shock conditions and derive the balance between coherent acceleration by the shock convection electric field and radiative losses due to synchrotron and inverse-Compton cooling. This balance determines both the maximum electron energy and the resulting steady-state spectrum. By forward-modeling the associated synchrotron emission and comparing it with integrated radio observations of the Sausage and Toothbrush relics, we find that the observed spectral curvature and high-frequency steepening can be reproduced when only a very small fraction ($\sim 10^{-9} - 10^{-8}$) of the available BSA acceleration capacity contributes to systematic electron energization. Despite this extremely small efficiency, it is sufficient to accelerate electrons to Lorentz factors $\gamma \sim 10^4 - 10^5$ under cluster conditions. These results suggest that radio relics provide a promising astrophysical laboratory for probing coherent acceleration, and that the BSA framework may account for the production of relativistic electrons in cluster shocks.

Infall and outflows, coupled with magnetic fields, rapidly structure the gas around newborn protostars. Shocks from interacting components encode the temperature and density distribution, offering a direct probe of the earliest evolution history. However, interferometric observations characterizing warm envelopes using high-excitation lines remain scarce. We present ALMA Band 9 observations of the Taurus dense core MC 27/L1521F, which hosts a Class 0 protostar, targeting the CO($J$=6-5) line at an angular resolution of $\sim$2\arcsec\ ($\approx$300 au). We detect an off-centered ring-like structure with a diameter of $\sim$1000 au that was not identifiable in previous low-$J$ CO data, where emission close to the systemic velocity is strongly affected by optical depth. The ring shows a typical peak brightness temperature of $\sim$3 K at our resolution. Excitation considerations indicate that the detected CO($J$=6-5) emission likely arises from relatively warm ($T \gtrsim 20$ K) and dense ($n({\rm H_2}) \gtrsim 10^{5}$ cm$^{-3}$) gas embedded within the surrounding cold, dense core. The morphology and kinematics suggest an energetic and localized shock-heating event, potentially linked to dynamical gas--magnetic-field interactions in the earliest protostellar phase. Our results demonstrate that high-$J$ CO observations provide a powerful new window on warm and dense gas components, enabling a more direct view of the physical processes operating at the onset of star formation.

Spectral-line results from a new cryogenic phased array feed (cryoPAF) on the Murriyang telescope at Parkes are presented. This array offers a significant improvement in field of view, aperture efficiency, bandwidth, chromaticity and survey speed compared with conventional horn-fed receivers. We demonstrate this with measurements of sky calibrators and observations of 21-cm neutral hydrogen (HI) in the LMC and the nearby galaxy NGC 6744. Within 0.3 deg of the optical axis, the ratio of system temperature to dish aperture efficiency is 25 K and the ratio with beam efficiency is 21 K (at 1.4 GHz). For the previously measured $T_{sys} = 17$ K, respective efficiency values 0.7 and 0.8 are derived. Our HI observational results are in good agreement with previous results, although detailed comparison with multibeam observations of the LMC suggests that the earlier observations may have missed an extended component of low-column-density gas ($8\times 10^{18}$ cm$^{-2}$). We use the cryoPAF zoom-band and wideband data to make a preliminary investigation of whether the large number of simultaneous beams (72) permits the use of novel data reduction methods to reduce the effects of foreground/background continuum contamination and RFI. We also investigate if these methods can better protect against signal loss for the detection of faint, extended cosmological signals such as HI intensity maps. Using robust higher-order singular value decomposition (SVD) techniques, we find encouraging results for the detection of both compact and extended sources, including challenging conditions with high RFI occupancy and significant sky continuum structure. Examples are shown that demonstrate that 3D SVD techniques offer a significant improvement in noise reduction and signal capture compared with more traditional layered 2D techniques.

In this whitepaper, we outline how recent technological advances and ongoing developments open qualitatively new science opportunities in cosmology, fundamental physics, and quantum astrophysics. First, intensity interferometry can contribute to one of the most foundational observables in cosmology: the expansion rate of the Universe. Its angular resolution allows it to resolve the angular extent of extragalactic objects such as supernovae or quasars; combined with a physical scale local to the source, this yields an angular diameter distance and hence a 'Hubble diagram'. Second, the nature of dark matter can be probed via the astrometric lensing signatures of tiny dark matter halos. Third, intensity interferometry gives direct access to second-order coherence properties of astrophysical emission, opening a window onto genuinely quantum aspects of astrophysical light.

Asymptotic Giant Branch (AGB) stars significantly contribute to the chemical composition of the universe. In their outflows, complex chemistry takes place, which critically depends on the local temperature. Therefore, if we want to accurately model the AGB environment, we need accurate cooling rates. The CO molecule is abundant in AGB outflows, and has a dipole moment, which enables it to cool through emission from its rotational transitions. We therefore expect it to significantly contribute to cooling in this environment, even at low temperatures ($10$ K $\leqslant T\leqslant 3000$ K). Currently, CO cooling rates are available for ISM-like conditions, which encompasses a different parameter regime, with generally lower densities and velocity gradients, compared to AGB winds. Therefore, these ISM cooling rates might not be applicable to the AGB regime. In this paper, we compute CO cooling rates for hydrodynamics simulations of AGB outflows. To evaluate the net cooling rate, we calculate the energy level distribution of CO self-consistently, using the non Local Thermodynamical Equilibrium (NLTE) line radiative transfer code Magritte. We verify whether already existing CO cooling rate prescriptions for the interstellar medium (ISM) are applicable for this regime. We noticed minor differences between these prescriptions and our calculated cooling rates in general. However, when used far outside their originally intended parameter regimes, significant differences occur. Therefore, we propose a new CO cooling rate prescription for the AGB environment and we study how the computed cooling rate varies depending on input parameters.

We present PrestaLine: Gorynych, a comprehensive numerical tool designed to model the dynamical, chemical, and spectral evolution of collapsing molecular cores from the prestellar phase to protostellar accretion. The code integrates three key components: (1) Kamelung, a 1D hydrodynamics module simulating gravitational collapse up to first hydrostatic core (FHSC) formation followed by an accretion of envelope onto a young star; (2) Presta, a chemical evolution module post-processing density and temperature profiles to compute time-dependent molecular abundances; and (3) Uran(IA), a radiative transfer module generating synthetic molecular line spectra. We apply Gorynych to compare low-mass (5 Msun) and high-mass (50 Msun) cores, finding that their dynamical evolution is remarkably similar, with differences primarily in spatial scaling. Chemical evolution reveals sharp abundance changes during the FHSC transition, particularly for CO, H2O, and HCO+, though pre-collapse chemical initialization has minimal impact on most species. Spectral maps of 13CO(2-1), HCO+(3-2), and H2O(110--101) lines show distinct kinematic signatures of infall and depletion, with high-mass cores exhibiting spatially extended emission. Our results highlight Gorynych's utility to couple theoretical collapse models and observations, providing a framework to diagnose core evolution and initial conditions from molecular line data.

We present an overview of CHIME/Slow, a real-time transient search backend under development to search for second-duration radio transients using the CHIME telescope, and results obtained from a pilot survey carried out using the prototype version of the search pipeline. The prototype CHIME/Slow pipeline was tested on archival data obtained in December 2022, January 2023 and February 2023 with a total on-sky time of 17 days with an instantaneous Field of View (FoV) of $\sim$13 deg$^2$ . In this pilot survey, we detected nine bursts, one from a new non-repeating source and eight from the known hyperactive repeating source FRB 20220912A. Out of these nine bursts, two bursts from the repeater were not detected by CHIME/FRB, while the non-repeater was detected in the side-lobe of a beam in the CHIME/FRB exhibiting shorter pulse width and narrower bandwidth compared to the CHIME/Slow detection. Here we report properties of the bursts, discuss the sensitivity and completeness of the current version of the CHIME/Slow pipeline, and outline future development to improve its performance. Finally, based on these results, we report the all-sky rate (95% credible region) of radio transients with pulse widths between 16 ms to 5 s, fluence above 5 Jy ms and observing frequency of 600 MHz to be between 184 and 4556 bursts sky$^{-1}$ day$^{-1}$.

We report the multi-temperature structure of the intracluster medium (ICM) in the Centaurus cluster core observed with XRISM/Resolve. Thanks to its high energy resolution, Resolve enables us to measure fine structures of highly ionized emission lines from Si to Fe and to directly determine the excitation temperature and the ionization temperature from the emission line ratio diagnostics. The observed spectrum in the Centaurus core is well-represented by a double-temperature thermal plasma at collisional ionization equilibrium state rather than an isothermal one. The line ratio diagnostics also support this biphasic temperature structure. Particularly, the observed line ratios show a trend of increasing ionization temperature with atomic mass, while the ionization and excitation temperatures of Fe show nearly the same temperature. The resultant line ratios, which are well-represented by the two temperatures ICM, ~ 1.6 and ~ 3 keV, are also fairly consistent with the expected numbers when assuming the radial single-temperature ICM was projected in the cluster core along the line of sight. Due to the limited low-energy sensitivity of the Resolve with the gate valve closed, we investigated the effect of the cool component using the XMM-Newton/RGS spectrum, but it ultimately did not affect our results. The observed flux ratio between the Fe XXV He alpha resonance and forbidden lines shows an about 20% reduction, suggesting the presence of resonant scattering.

The `Inverse Distance Ladder' uses relative-distance measurements with type Ia supernovae (SNe Ia), anchored to an absolute distance scale from Baryon Acoustic Oscillations (BAO) and the cosmic microwave background (CMB), to provide an alternative measurement technique for the local expansion rate $H_0$. Using SNe Ia from the Dark Energy Survey and BAO measurements from the Dark Energy Spectroscopic Instrument, we show that the inverse distance ladder is unable to explain the Hubble Tension, anchored as it is to the CMB and its value of $H_0 = 67.4 \pm 0.5$ kms$^{-1}$ Mpc$^{-1}$. To do so, we first show that the suite of systematics considered in cosmology analyses with SNe Ia only move the inferred $H_0$ by $<0.1$kms$^{-1}$ Mpc$^{-1}$, and second, we investigate the scale of redshift-dependent magnitude changes necessary to change the inferred inverse distance ladder $H_0$ from $67$ to the local network of distance measurements value of $73$, and the impact that this would have on other cosmological inferences with SNe Ia. We find that a change of $d\mu/dz = 0.2$ mag would be necessary to infer an $H_0$ in concordance with local distance measurements, and that this $d\mu/dz$ value would result in a Flat $\Lambda$CDM $\Omega_M = 0.23$, $10\sigma$ discrepant with other cosmological probes, {indicating that the precision of SNe Ia measurements preclude the necessary redshift evolution for an $H_0$ of 74 kms$^{-1}$ Mpc$^{-1}$}. Therefore, we conclude that current SN Ia cosmology leaves little freedom for the inverse distance ladder to yield $H_0$ values significantly different from $67$ kms$^{-1}$ Mpc$^{-1}$.

The IceCube neutrino telescope has provided upper limits on neutrino emission from gamma ray bursts. These constraints provided by the IceCube detector have been instrumental in investigating the properties of the GRB jet and its emission models. During the prompt phase of gamma ray burst emission, intense radiation components are generated that interact with the shock-accelerated particles within the jet. We study various GRB emission models, such as the internal shock model, the photospheric models, and also include a model-independent case. Based on these models, we calculate the neutrino fluence using the photo-hadronic interaction process. We estimate the bulk Lorentz factor using the well-known correlations between prompt phase observables, which is then used to calculate the emission site for the model-dependent scenarios. For GRB 221009A, we find that a low baryon loading scenario is consistent with the IceCube upper limits; however, for GRB 240825A, a higher value of baryon loading is preferred. Also, the values of the microphysical parameters $\epsilon_e$ and $\epsilon_B$ for GRB 240825A are lower by factors of approximately 10 and 100, respectively, compared to those of GRB 221009A. Further, using neutrino upper limits for these two sources, we estimate the lower limits on the dissipation radius for our models. The current TeV PeV upper limits for GRB 221009A are already useful for constraining parameter space for the BPH and MPH models.

The main objective is to calibrate the OPT-IR luminosity relation for quasars, focusing on accurate estimations of dusty torus and accretion disk luminosities. We analyzed contaminations related to host galaxies, particularly from polar dust, the interstellar medium, and stellar emission that affect the optical and infrared. We used a sample of nearly 400 quasars with photometrical observations and spectroscopical redshift divided into four redshift bins (0.7-2.4). Full spectral energy distribution (SED) fitting was performed with the CIGALE code, and results were compared with simplified photometric luminosity estimates. The impact of non-active galactic nucleus components and the role of polar dust in the fitting process were assessed. We show that for sources with a disk luminosity above 10^45 [erg/s], the photometric estimates are consistent with SED-based values. While polar dust contributes marginally to luminosity, its presence significantly alters SED fitting, particularly the torus opening angle and cold dust properties. In the optical domain, stellar emission is the dominant contamination. In the infrared, disk emission and cold dust play major roles. We propose two empirical calibrations for the OPT-IR relation. We conclude that the optical band is dominated by the accretion disk component above 10^45 or 10^46 [erg/s] depending on redshift, while IR luminosity is dominated by the dusty torus emission above 1.6 $\times$ 10^45 or 2 $\times$ 10^46 [erg/s] depending on the redshift. In this high-luminosity regime, simplified photometric methods yield reliable disk and torus luminosity estimates. The aim of the analysis we present is to test the parameter space in order to build a well behaving OPT-IR nonlinear luminosity relation for quasars that could serve as a cosmological probe.

We present new constraints on ultralight axion-like particles (ALPs) through multi-epoch measurements of differential birefringence induced due to a coupling ($g_{a\gamma}$) between the ALP and electromagnetic fields. Broadband polarimetric observations in the 2-8 GHz range of the gravitationally lensed system CLASS B1152+199 were carried out over five epochs spanning three months with a cadence of roughly 20 days, and the differential birefringence angle ($\Delta\,\theta_{a,{\rm lens}}$) between the lensed images were estimated. We also combined an archival observation that effectively increases the span to 9.5 yr to probe the effect of an oscillating ALP field imprinted as oscillating ${\Delta\,\theta_{a,{\rm lens}}}$ over time. Here we present a new technique for combining multi-epoch measurements of ${\Delta\,\theta_{a,{\rm lens}}}$ by considering the coherence of the ALP field, such that, ${\Delta\,\theta_{a,{\rm lens}}}$ over these observations are related. The time scale of coherence depends on the mass of the ALP field ($m_a$). With these new observations, we constrain $g_{a\gamma} \leq 9.0\times 10^{-12} \,\left( {\rho_{a,\text{em}}}/{20 \text{ GeV cm}^{-3}} \right)^{-1/2}\;\mathrm{GeV}^{-1}$ to $\leq 3.5\times 10^{-8} \,\left( {\rho_{a,\text{em}}}/{20 \text{ GeV cm}^{-3}} \right)^{-1/2}\;\mathrm{GeV}^{-1}$ for $m_a$ between $1.6\times 10^{-22}\;\mathrm{eV}$ and $3.8\times 10^{-18}\;\mathrm{eV}$, where $\rho_{a,{\rm em}}$ is the density of the ALP field at emission. This improves over the constraint provided by the CERN Axion Solar Telescope by up to an order of magnitude in the $m_a$ range $1.6\times 10^{-22}\;\mathrm{eV}$ to $3\times 10^{-21}$ eV.

We presented the first high-precision, detailed photometric analysis of four W Ursae Majoris (W UMa)-type contact binaries, Linear 10772300, Linear 11150338, Linear 20372537 and DM Cir. In addition to ground-based multiband photometric observations, data from the Transiting Exoplanet Survey Satellite (TESS) were employed for the analysis of the DM Cir system. New ephemeris and linear fit to the O-C diagrams were derived using extracted times of minima and additional literature. The light curve modeling was performed using the PHysics Of Eclipsing BinariEs (PHOEBE) Python code and the BSN application, employing a Markov Chain Monte Carlo approach. In each systems, the two stellar components exhibited minimal temperature differences ($\Delta T<150$ K), confirming efficient energy exchange within their common convective envelopes. Absolute parameters were estimated using the Gaia Data Release 3 (Gaia DR3) parallax and astrophysical equations. Based on effective temperatures and component masses, two systems were classified as W-subtype systems, while others belonged to the A-subtype. We computed the initial masses of the primary ($M_{1i}$) and secondary ($M_{2i}$) components for four target systems using a method based on the observational properties of overluminous secondary components. We found initial primary masses in the range 0.6-1.0$M_\odot$ and initial secondary masses in the range 0.9-1.7$M_\odot$ with mass loss $<1.0M_{\odot}$. We investigated the relative energy transfer rates ($U_{1}$ and $U_{2}$) and nuclear luminosities ($L_{10}$ and $L_{20}$) based on the physical parameters of 411 W UMa-type contact binaries, including the four systems analyzed in this study, through wide range of mass ratios. The results for all systems provided a comprehensive view of energy transfer behavior throughout different evolutionary stages of contact binaries.

Genetic Algorithms (GA) are a powerful tool for stochastic optimisation and non-parametric symbolic regression, already widely used in cosmology. They are capable of reconstructing analytical functions directly from data points without introducing new physical models. A limitation of this approach is that while the reconstructed function is very efficient at reproducing the behaviour of the data points, non-observable quantities involving derivatives are particularly sensitive to stochasticity, hyperparameters, and to the choice of the best-fit function obtained by the GA, which implies the risk of the algorithm getting stuck in a local minimum. In this work we propose an update to the GA methodology for the reconstruction of analytical functions that involves computing a weighted average of an ensemble of GA configurations (\texttt{GAME}). We define the weights via a quantity that accounts for both the goodness-of-fit of the points and the smoothness of the resulting function. We also present a practical method to analytically estimate and correct the errors on the averaged function by combining a path-integral approach with an ensemble variance. We demonstrate the improvement offered by \texttt{GAME} methodology on a generic test function. We then apply the new methodology to a non-parametric reconstruction of the Hubble rate $H(z)$ using Cosmic Chronometers data and, assuming a flat Friedmann-Lemaître-Robertson-Walker background and General Relativity, we infer the corresponding dark energy equation of state $w(z)$. Through consistency tests, we show that current data produces results compatible with $\Lambda$CDM, and that Stage IV cosmology surveys will allow GA reinforced with \texttt{GAME} methodology to become an even more competitive tool for discriminating between different models.

Wolf-Rayet (WR) stars have recently attracted attention as possible drivers of early chemical enrichment, including the production of fluorine, whose nucleosynthetic origin remains debated. To test the contribution of massive stars to fluorine production in the early Universe, we conducted Atacama Large Millimeter/submillimeter Array Band~5 spectroscopy of the HF(1-0) absorption line toward a dusty star-forming galaxy at $z=6.024$. This galaxy has a known gas-phase metallicity and is too young for low-mass AGB stars to have contributed significantly, providing a clean environment to isolate massive-star yields. We do not detect significant HF absorption ($\sim2\sigma$) and derive a conservative 5$\sigma$ upper limit of $N_\mathrm{HF}/N_\mathrm{H_2} < 2.2\times10^{-9}$. This limit is about an order of magnitude below typical local measurements, indicating inefficient fluorine enrichment $\sim0.9$\,Gyr after the Big Bang. Comparison with chemical evolution models shows that our constraint is consistent with scenarios without WR yields at this epoch. Expanding the sample of HF absorption measurements in high-redshift galaxies with well-characterized metallicities will be crucial for tracing the onset of WR enrichment and fluorine production across cosmic time.

Substructures such as rings, gaps, and cavities are commonly observed in protoplanetary discs and are thought to play a key role in dust evolution and planet formation. However, a fraction of the extended discs (68% dust radii > 30 AU) in nearby star-forming regions remain unresolved, leaving their substructure content uncertain and thereby limiting our understanding of dust evolution and the initial conditions for planet formation across the full disc population. We aim to investigate the presence of substructures in previously unresolved, extended discs to assess whether all extended protoplanetary discs in the Solar neighbourhood exhibit substructures. We present new high-resolution ($\sim$0.12") ALMA Band 6 continuum observations at 1.33 mm of 26 previously unresolved, extended discs within 200 pc, completing the high-resolution sample of extended discs in Taurus, Ophiuchus, Chamaeleon, Lupus, Upper Scorpius, Upper Centaurus-Lupus and Lower Centaurus-Crux. We analyse radial intensity profiles using Frankenstein and Galario to detect substructures. Seventeen discs show clear substructures, while nine appear compact and structureless, smooth or ambiguous due to inclination or possible binarity/late-stage infall. We detect $^{12}$CO J=2-1 emission in 15 discs, with extended emission in four. Combined with literature data, our complete sample of 730 protoplanetary discs reveals that nearly all extended discs exhibit substructures, $\sim$91% detected in the full sample, and up to $\sim$98% when correcting for high-inclination systems where substructures may be hidden. Substructures are a near-universal feature of extended protoplanetary discs. They are more commonly detected in larger, massive discs and around higher-mass stars, and structured discs retain their dust mass over time. This supports the scenario in which dust traps, possibly induced by giant planets, shape disc morphologies.

We present a feasibility study exploring the implementation of optical interferometry and speckle techniques with the 1.3-m Devasthal Fast Optical Telescope (DFOT) at ARIES, which is currently dedicated to photometric observations. Using the sCMOS camera as the DFOT backend, we perform interferometric speckle observations of several binary stars. Standard Speckle Interferometry (SI) algorithms are applied to analyze the recorded data. While this study does not aim to achieve the diffraction limit of DFOT or address a full science-driven resolution case, it serves as a crucial testbed for instrumentation, data acquisition, and analysis of Speckles with DFOT. Notably, we successfully identify and correct tracking-related positional errors in the observed binary systems, demonstrating the viability of the approach. These results provide strong motivation for more systematic observations and future implementation of optical interferometry techniques at meter-class telescopes.

The Magellanic Clouds (MCs) offer a unique laboratory for studying galaxy interaction and the evolution of dwarf galaxies. By investigating when and how stars formed, the star formation history (SFH) is a powerful tool to provide constraints for dynamical modeling of the system's past interactions and understand the processes of stripping and triggered star formation in tidally influenced environments. We aim to reconstruct the SFH of the Magellanic Bridge, the gaseous and stellar stream connecting the two Clouds. We used data from the deep optical STEP survey, which covers 54 $\mathrm{deg\, {^{2}}}$ across the Small Magellanic Cloud (SMC) and the Bridge, reaching stars below the oldest main sequence turnoff at the distance of the MCs. We applied the synthetic color-magnitude diagram (CMD) technique to 14 deg$^2$ of STEP data. We constructed two libraries of synthetic stellar populations based on the PARSEC-COLIBRI and BaSTI stellar evolutionary models, with metallicities in the range $-2.0\leq[$Fe/H$]\leq0$ across the whole Hubble time. We find a clear peak of recent star formation $\sim100$ Myr ago in the Magellanic Bridge, which becomes increasingly pronounced toward the SMC. The low metallicity of this population suggests that it formed from gas stripped from the SMC during its most recent close encounter with the LMC. In the eastern part of the Bridge (LMC side), the star formation peaks at earlier times, around 10 Gyr and 2 Gyr ago. We estimate a total stellar mass in the Bridge of $ (5.1 \pm 0.2) \times 10^5 M_\odot$ and a present-day stellar metallicity of $[$Fe/H$]\sim0.6$ dex, close to SMC value.

We present time-series medium-band (R~20-40) observations of the third interstellar object 3I/ATLAS (C/2025 N1) obtained with the 7-Dimensional Telescope (7DT), enabling spatially resolved monitoring of its gas and dust activity from 2025 July to September. The m400-band image (lambda_c = 400 nm, Delta lambda approx 25 nm) reveals the emergence of pronounced and spatially extended CN emission at heliocentric distances r_h < 3 au. This onset is consistently identified across multiple diagnostics, including a break in the light-curve evolution, excess reflectance, inward expansion of annular excess beyond 10,000-20,000 km, growth of the coma half-light radius from ~11,000 to ~19,000 km, and a rapid rise in the CN production rate Q_CN relative to Af rho. We further separate the CN-emitting and dust-scattered components through two-dimensional surface-brightness fitting into inner (dust) and outer (gas) components. The outer component preserves a nearly constant profile shape, varying only in normalization, implying relatively fast expansion of CN-bearing molecules. Together, these results reveal a transition in the optical from dust-dominated scattering at large heliocentric distances to volatile-driven, gas-dominated activity as 3I/ATLAS enters the inner Solar System. The timing and characteristics of the CN activation resemble the volatile enhancement observed in 2I/Borisov, suggesting that both known active interstellar objects exhibit comparable activation behavior at heliocentric distances of ~2-3 au.

Aim. We present and interpret the data of the nearby hydrogen-deficient but helium-rich superluminous supernova SN~2021bnw which reached a magnitude of -20.7 at maximum luminosity in g band. Methods. We discuss the light curves and spectra of SN 2021bnw based on its spectro-photometric follow up exploiting different observational facilities. We reproduce the NIR spectrum of SN 2021bnw with TARDIS to inspect the chemical composition at late photospheric phases and identify helium features. We also use a STELLA model coupling hydrodynamics and radiation transport to constrain the physical parameters of the explosion assmunig a 56Ni+CSM scenario. Results. We suggest that SN 2021bnw was mainly powered by the interaction of the ejecta with a previously lost He-rich circumstellar material, coupled with a central power source. Conclusions. This work expands the data sample of He-rich superluminous supernovae rich (SLSNe Ib) and, assuming a single progenitor scenario, can constrain the masses and the physics of their progenitors.

Modern time-domain surveys produce alert streams at a scale that makes exhaustive manual inspection infeasible, requiring automated methods to identify unusual transients for follow-up. In this work, we present an unsupervised anomaly detection pipeline applied to the ZTF alert stream using the Lasair broker. We define normal objects as SN Ia, SN II, and SN Ib/c. Anomalous objects include (i) more exotic transients (AGN, TDEs, SLSNe, CVs, and nuclear transients) and (ii) supernova-labeled objects, either spectroscopically or by Lasair, with anomalous properties, such as incorrect or absent host associations, or non-supernova-like light curves. Our pipeline consists of three independently trained simple autoencoders operating on distinct alert stream data products: object features, triplet image cutouts, and light curves. Each model is trained on predominantly normal transients, and performance is assessed using the recall of exotic objects and the purity of all anomalous objects across both a spectroscopically classified held-out test set and the live alert stream. In the test set, performance is evaluated at a fixed rank corresponding to the top ten scoring candidates, while in the alert stream it is evaluated using an anomaly threshold defined from test set behavior. Across both settings, the algorithms consistently recover exotic transients and anomalous supernovae among their top-ranked candidates. Over 25 days of live alert stream application, we identify 87 unusual supernova candidates for follow-up. The overlap between anomalies flagged by different autoencoders in the test set is non-existent, and in the alert stream is small, with maximum overlap between any two algorithms being 11 objects. The framework is data-efficient, requiring only a few thousand training examples, making it well suited for early and ongoing application to the Rubin Observatory alert stream.

The recent discovery of three off-nuclear tidal disruption events (EP240222a, AT2024tvd, and AT2025abcr) - following the first such source, 3XMM J2150$-$05 - reveals a small but robust population of off-nuclear, or `wandering', black holes (WBHs) with masses $M_\bullet > 10^4 M_\odot$. Two demographic trends are already apparent: (i) all events occur in massive, early-type parent galaxies with stellar masses $10.8 \lesssim \log_{10}(M_\star/M_\odot) \lesssim 11.1$; and (ii) events at larger halo-centric radii ($R_{\rm TDE}/R_{200}$) are associated with dwarf satellites ($M_\star \sim 10^7 M_\odot$), while those closer to halo centers lack detected stellar counterparts. Using results from the \texttt{ROMULUS} cosmological simulation, we show that both trends naturally arise from hierarchical galaxy formation. By combining the simulation with empirical constraints on the local galaxy population, we compute the volumetric density of WBHs, $\phi_{\rm WBH}(M_\star)$, finding that it peaks at $\log_{10}(M_\star/M_\odot)=11.10^{+0.05}_{-0.10}$ and that more than half of all WBHs in the local Universe reside in galaxies with $10.7 \lesssim \log_{10}(M_\star/M_\odot) \lesssim 11.2$, explaining (i) and predicting its persistence as the sample grows. We further show that ii), i.e., the observed link between detection of stellar counterparts and $R_{\rm TDE}/R_{200}$, is also expected from tidal stripping. These results demonstrate that off-nuclear TDEs are powered by the population of WBHs long predicted by cosmological simulations.

We present the discovery of TIC-65910228 b / NGTS-38 b, a giant exoplanet with a radius of $1.081\pm0.047$ R$_\text{J}$ and a mass of $4.78_{-0.37}^{+0.40}$ M$_\text{J}$ on a long-period ($180.52791\pm0.00038$ day), moderately eccentric ($e=0.308\pm0.011$) orbit transiting a bright (V=$10.230\pm0.020$ mag) metal rich ([Fe/H]=$0.33\pm0.09$, 'dex') F6V-F7V type host star. The planet was initially detected from a single transit in TESS Sector 33. A photometric monitoring campaign of 228 nights with NGTS detected a transit egress of the planet, which together with spectroscopic radial velocity monitoring with CORALIE and HARPS identified an orbital period of ~180.5,d. These radial velocity measurements also showed the mass of the companion to be planetary. Additional transit observations coordinated by the TESS follow-up observing program allowed further confirmation and refinement of this period. With its relatively cool equilibrium temperature of $458\pm11$,K, NGTS-38 b joins a small but growing population of well characterised transiting warm-Jupiters and has one of the longest periods of any discovered to date. The target is situated in the LOPS2 field of the upcoming PLATO mission which will allow for greater refinement of the system parameters and potential for the discovery of additional companions too small and/or too long-period to be seen by TESS or NGTS. NGTS-38 b's bright host star and wide orbital separation make it an attractive target for further study, including potential measurement of its spin-orbit alignment or targeted exomoon/ring searches.

In preparation for development of both key technologies and instrument concept studies to use those technologies, the Habitable Worlds Observatory Technology Maturation Project Office at the NASA Goddard Space Flight Center has compiled a series of statements of state of the art for those same key technologies. These statements are being provided to the public as exemplars and suggestions for possible future collaboration for those same instrument concept studies, but without mandate, to enable proposing teams to be able to find the technical solutions they need to field a compelling proposal. This information resides in the public domain and is presented without prejudice.

The unprecedented rest-frame UV and optical coverage provided by JWST enables simultaneous constraints on the electron density (n$_{\rm e}$) and temperature (T$_{\rm e}$) of ionized gas in galaxies at z>5. We present a self-consistent direct method based on multiple OIII]1661,66) and [OIII] ($\lambda$4363, and $\lambda$5007) transitions to characterize the physical conditions of the high-ionization zone. This new approach is insensitive to a wide range of n$_{\rm e}$ due to the high critical densities of the OIII] and [OIII] transitions. Applying this technique to six galaxies at z=5-9, we find electron densities up to n$_{\rm e}$$\sim 3\times 10^{5}$ cm$^{-3}$ and temperatures of T$_{\rm e}$ $\sim 20,000$ K in systems at $z>6$. Accounting for these self-consistent densities changes the derived T$_{\rm e}$ and modifies the inferred metallicities by up to 0.29 dex relative to previous estimates. We discuss the reported N/O overabundances in the high-$z$ galaxies from our sample, which arise entirely from the high N$^{3+}$/H$^{+}$ values inferred from NIV] lines. We point out that a T$_{\rm e}$-stratification, in which the N$^{3+}$ zone has a slightly higher T$_{\rm e}$ than T$_{\rm e}$([OIII]), could substantially reduce the inferred N/O. Quantitatively, if T$_{\rm e}$(N$^{3+}$) were 10\% higher than T$_{\rm e}$([OIII]), this could induce a systematic overestimation of N$^{3+}$/O$^{2+}$ of nearly 50\%. Classical N/O diagnostics such as N$^{+}$/O$^{+}$, due to their critical densities, can significantly impact the inferred N/O abundance in the presence of high-density gas, whereas N$^{2+}$/O$^{2+}$ place these galaxies closer to $z\sim0$ systems in the N/O-O/H plane. Future JWST programs with larger and more diverse samples will be essential to test the universality and robustness of these results.

Modern time-domain surveys, such as the Zwicky Transient Facility (ZTF), detect far more extragalactic transients than can be spectroscopically classified. Photometric classification offers a scalable alternative, enabling the identification of larger, fainter, and higher-redshift supernova samples suitable for applications such as Type Ia supernova (SN Ia) cosmology. We present a feature-based photometric classifier for SNe detected by ZTF, with the primary goal of constructing a photometric SN Ia sample for cosmological analyses. Our approach utilises the autoencoder architecture of ParSNIP (Boone 2021) to capture the intrinsic diversity of SN light curves. We trained the model on a spectroscopically classified ZTF SN sample, incorporating a realistic noise augmentation procedure that simulates the flux uncertainties of fainter sources. Light curve features were used to train a gradient-boosted decision tree classifier, implemented in both binary (SN Ia vs. non-Ia) and multi-class configurations. We validated our classifier on independent, fainter ZTF data with and without noise augmentation. To evaluate real-time performance, we also applied our classifier to live ZTF alerts and conducted a spectroscopic classification survey within the ePESSTO+ collaboration. We found that noise augmentation significantly improves classification performance, particularly for fainter sources. Our binary classifier achieves an SN Ia recall of (98.1 $\pm$ 0.4)%, averaged across five train-test splits. SN Ia recall exceeds 98% for events with a peak apparent magnitude up to 20 and more than 10 detections, and remains above 96% up to magnitude 20.5. Overall, 95% of sources were correctly classified in both binary and multi-class modes. Our classifier performs efficiently on real ZTF data and enables construction of a large photometric SN Ia sample for cosmology.

Observational appearance of debris disks is largely controlled by collisional grinding of their dust grains. However, the mechanical strength of dust at sizes in the micrometer to millimeter range is poorly known. Recent studies suggested that dust particles in the Solar system might have a higher critical fragmentation energy $Q_{D}^*$ value than previously anticipated. Another recent study considered the Fomalhaut debris disk and found lower $Q_{D}^*$ values to provide better fits to the data. In order to constrain the mechanical strength of dust, we investigate collisional evolution of debris disks with $Q_{D}^*$ prescriptions differing by $\sim 3$ orders of magnitude. We find that, above a certain threshold $Q_{D}^*$ value, the disk's collisional evolution is dominated by rebounding -- rather than disruptive or cratering -- collisions. Rebounding (a.k.a. bouncing) collisions are those in which both impactors survive, being only slightly eroded and producing fragments that only carry a minor fraction of their mass. We show that disks dominated by rebounding collisions would have brightness profiles increasing outward outside the parent belt. Since such profiles are not observed, this places an upper limit on how hard the debris dust is allowed to be in order not to violate the observations. We derive an approximate analytic expression for this limit: $Q_{D}^* \approx (1/8) v_{K}^2(r)$ for grains close to the radiation pressure blowout size, where $v_{K}$ in the Keplerian circular speed at a distance $r$ from the star. This implies $Q_{D}^* \lesssim 10^{9...10} \,\hbox{erg}\,\hbox{g}^{-1}$ for micrometer-sized grains in typical debris disks. Even though rebounding collisions are not expected to affect debris disk evolution significantly, we emphasize that these collisions are actually much more frequent than disruptive and cratering ones in all debris disks.

Post-starburst (PSB) galaxies, having recently experienced a starburst followed by rapid quenching, are excellent laboratories to probe physical mechanisms that drive starbursts and shutting down of star formation. Integral-field spectroscopy reveals the galaxies' spatially-resolved properties, where observed directional patterns can be linked to the galaxies' past evolution. We measure the resolved star-formation histories (SFHs), stellar metallicity evolution and dust properties of three local PSBs from the MaNGA survey, down to $0.5$" resolution ($\sim0.3\,$kpc) using a hierarchical Bayesian model. Local parameters were constrained simultaneously with parameters describing spatial trends. We found that all three galaxies first experienced an outer, weaker and slower quenching starburst, followed by a central, stronger and faster quenching starburst that peaked $\sim 1\,$Gyr after the first. The central starbursts induced a significantly stronger rise in stellar metallicity compared to the outer starbursts. These results are consistent with the effects of a recent gas-rich (wet) merger, where the first pericentre passage triggered starbursts in the outer regions, while the later coalescence triggers a stronger centralised starburst. We find non-axisymmetric features in the maps of burst mass fraction and dust attenuation in all galaxies, which could be caused by tidal effects during the recent merger. Comparisons with literature binary merger simulations suggests that the galaxies' rapid quenching was driven by gas consumption and the stabilisation against gas gravitational collapse by a growing spheroid, while AGN feedback was not necessarily a primary cause.

We investigate how the H$_2$ ortho-to-para ratio (OPR) and dueterium fractionation in star-forming regions are affected by nuclear spin conversion (NSC) on dust grains. Particular focus is placed on the rotational energy difference between ortho-H$_2$ (o-H$_2$) and para-H$_2$ (p-H$_2$) on grain surfaces. While the ground state of o-H$_2$ has a higher rotational energy than that of p-H$_2$ by 170.5 K in the gas phase, this energy difference is expected to become smaller on solid surfaces, where interactions between the surface and adsorbed H$_2$ molecules affect their rotational motion. A previous study by Furuya et al. (2019) developed a rigorous formulation of the rate for the temporal variation of the H$_2$ OPR via the NSC on grains, assuming that adsorbed o-H$_2$ has higher rotational energy than adsorbed p-H$_2$ by 170.5 K, as in the gas phase. In this work, we relax the assumption and re-evaluate the rate, varying the rotational energy difference between their ground states. The re-evaluated rate is incorporated into a gas-ice astrochemical model to study the evolution of the H$_2$ OPR and the deuterium fractionation in prestellar cores and the outer, cold regions of protostellar envelopes. The inclusion of the NSC on grains reduces the timescale of the H2 OPR evolution and thus the deuterium fractionation, at densities of >10$^4$ cm$^{-3}$ and temperatures of <14-16 K (depending on the rotational energy difference), when the ionization rate of H$_2$ is 10$^{-17}$ s$^{-1}$.

Double white dwarfs (DWDs) are by far the most common compact binaries in the Milky Way, are important low-frequency gravitational-wave sources, and in some cases merge to become Type Ia supernovae. So far, no DWD has been identified solely through relativistic Doppler beaming, even though the beaming amplitude directly relates to the radial velocity semi-amplitude. In this work, we initiate a comprehensive binary population synthesis using SeBa and incorporate the resulting binaries into a tripartite Galaxy model. Our proof-of-concept simulations demonstrate that the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) can reliably recover relatively bright ($r \lesssim20~$mag) unequal-mass binaries in compact orbits with P $\approx$ 10-600 minutes with moderate to high inclinations. We find that LSST can detect at least 287 short-period DWDs, of which 47 are LISA-detectable gravitational wave sources. LSST lightcurves allow us to readily determine the period and fully characterize the orbit, in contrast with the challenges of orbit determination for DWDs in spectroscopic searches. The formation of unequal mass, short-period DWDs strongly depends on the assumptions regarding the mass-transfer phases during binary population synthesis, and the total number and characteristics of Doppler-beamed DWD systems observed in LSST will provide new tests of models of stellar binary evolution. Here, we lay the foundation for the comprehensive integration of synthetic Galactic binary population into realistic LSST survey simulations, thereby enabling quantitative forecasts of the number and characteristics of any binary sub-population during the LSST era.

Fermi LAT observations of gamma-ray pulsars can be used to build a pulsar timing array (PTA) experiment to search for gravitational wave (GW) signals at nanohertz frequencies. At those frequencies, the dominant signal is expected to be a stochastic gravitational wave background (GWB) produced by the incoherent superposition of the quasi-monochromatic GW emissions from a population of supermassive black hole binaries. While the radio PTAs have recently announced compelling evidence for a GWB signal with a power law spectrum of strain amplitude $\approx2-3\times10^{-15}$ (at the frequency of $1 {\rm yr}^{-1}$), in 2022 an analysis of $12.5$ years of Fermi data for 35 pulsars led to an upper limit of $1\times10^{-14}$ for the GWB amplitude. The analysis was carried out on times-of-arrival (TOAs) obtained by folding from six months up to one year of photon observations. A photon-by-photon approach was also tested to infer constraints on the GWB amplitude from individual pulsars, but without accounting for the cross-pulsar correlations that a GWB would induce. Here, we reanalyse the same dataset using a regularized likelihood method that correctly models cross-pulsar correlations directly from the photons, while additionally marginalising over the uncertain pulse profile shape. While the two methods are not expected to have significant differences in sensitivity, we prove through simulations of gamma-ray PTA datasets that the photon-by-photon method for GWB recoveries is, statistically, more robust. The resulting upper limit obtained for the GWB strain amplitude is $1.2\times10^{-14}$, indicating that the improved method yields a consistent result with the previous analyses.

Using long-term observations from the Transiting Exoplanet Survey Satellite (TESS) along with spectroscopic observations from the 3.6 m Devasthal Optical Telescope (DOT), we present a comprehensive time-series and spectral analysis of the intermediate polar V515 And. Our analysis reveals that V515 And resides within the period gap, with the detection of its orbital period of 2.73116 h. Additionally, we confirm the earlier findings of the spin and beat periods to be 465.4721 s and 488.6067 s, respectively. The time-resolved timing analysis reveals that V515 And undergoes changes in its accretion geometry, not only between different TESS sectors but also within individual sector observations. The system exhibits a transition in the dominant accretion mode, switching between disc-fed and stream-fed accretion. In the TESS light curve, we identify two successive outburst-like episodes, each persisting for roughly a day and reaching peak luminosities of $2.7\times10^{33}$ and $1.9\times10^{33}$ erg s$^{-1}$. Our analysis suggests that these bursts belong to the recently proposed class of micronovae. The optical spectrum of V515 And is characterised by strong Balmer and He II emission lines and shows an inverse Balmer decrement indicating the magnetic nature of the source.

We introduce $\texttt{GPUmonty}$, a CUDA/C-based Monte Carlo radiative transfer code accelerated using graphics processing units (GPUs). $\texttt{GPUmonty}$ derives from the CPU-based code $\texttt{grmonty}$ and offloads the most computationally expensive stages of the calculation -- superphoton generation, sampling, tracking, and scattering -- to the GPU. Whereas $\texttt{grmonty}$ handles photons sequentially, $\texttt{GPUmonty}$ processes large numbers of superphotons concurrently, leveraging the single-instruction, multiple-thread (SIMT) execution model of modern GPUs. Benchmarks demonstrate a speedup of about $12\times$ relative to the original CPU implementation on a single GPU, with runtime limited primarily by register pressure rather than compute or memory bandwidth saturation. We validate the implementation through analytic tests for a optically thin synchrotron sphere, as well as comparisons with $\texttt{igrmonty}$ for scattering synchrotron sphere and GRMHD simulation data. Relative errors remain below a percent level and convergence is consistent with the expected $N_{\rm s}^{-1/2}$ Monte Carlo scaling. By significantly reducing computational costs, GPUmonty enables the extensive parameter space surveys and faster spectra modeling required to interpret horizon-scale observations of supermassive black holes. $\texttt{GPUmonty}$ is publicly available under the GNU General Public License.

We present a non-perturbative approach to correlation functions and polyspectra of locally non-Gaussian fields and develop a simple semi-perturbative framework that does not rely on the local expansion. As an example, we apply it to locally non-Gaussian fields possessing exponential tails and derive some exact analytic results in the strongly non-Gaussian limit.

Neutrino self-interactions (NSI) offer a potential pathway to address anomalies in standard cosmology and explain existing cosmological tensions. In this work, we present a novel framework to obtain the neutrino--neutrino collision term within the Boltzmann hierarchy, incorporating both neutrino and mediator masses as free parameters. Our calculations encompass both Dirac-like and Majorana neutrinos and distinguish between two neutrino mass eigenstates. This work provides a valuable tool for future analyses, should a NSI signal be detected. Remarkably, our results show a smooth transition from the light to the heavy mediator approximation as the Universe cools down for non-resonant cases. Thus, for the widely studied heavy mediator, our new scheme eliminates the need to approximate at high redshifts when the temperature increases above the mediator mass, and it provides the tools to test the threshold of validity of the heavy mediator paradigm. While this work focuses on NSI mediated by a scalar particle, the presented framework could be adapted to a broader range of neutrino NSI and possibly to warm dark matter self-interacting scenarios.

We study the effects of dark matter on the properties of hybrid neutron stars, in particular the influence on the mass-radius relation, the value of the maximum mass, and the hadron-quark phase transition. To single out the equilibrium configurations of dark-matter-admixed hybrid neutron stars (DHSs), we also study their radial oscillations. Both the stellar structure equations and the radial oscillation equations are solved for the two-fluid system, where the ordinary matter component and dark matter component couple only through gravity. For the ordinary matter components, we adopt the Brueckner-Hartree-Fock method for nuclear matter, and the Dyson-Schwinger or the field-correlator model for quark matter. For the dark matter component, we use a non-self-annihilating self-interacting fermionic model. We find that the presence of dark matter in DHSs leads to a decrease of the critical mass of the hadron-quark phase transition, a related possible onset of quark matter in dark-matter accreting stars, and a significant reduction of radial oscillation frequencies.

The ``hyperon puzzle'' remains a fundamental challenge in nuclear astrophysics. We investigate hyperon emergence in neutron star matter using the $SU(3)$ parity doublet model with chiral representation $(3,\bar{3}) + (\bar{3},3)$. This framework naturally incorporates chiral symmetry restoration and provides a systematic description of baryon masses in dense matter through the interplay between the chiral condensate and the chiral invariant mass $m_0$. We find that the hyperon onset density exhibits strong sensitivity to $m_0$: for $m_0 = 500$ MeV, hyperons first appear at $1.9n_0$ while for $m_0 \gtrsim 750$ MeV, hyperons emerge only above $5n_0$. This delayed onset arises from the weakened density dependence of baryon masses at larger $m_0$ values. When the hyperon onset density exceeds the expected quark-hadron transition range ($2$--$5n_0$), matter undergoes deconfinement before hyperons populate, avoiding the EoS softening while maintaining consistency with massive neutron star observations. Our results demonstrate that chiral dynamics provides a natural resolution to the hyperon puzzle without requiring ad hoc repulsive hyperon interactions.

We study the impact of light GeV-scale heavy neutral leptons (HNLs) on Big Bang nucleosynthesis (BBN) in the neutrino-extended Standard Model Effective Field Theory ($\nu$SMEFT). We show that, based on very general considerations, BBN constraints complement laboratory searches at colliders, beam dumps, and neutrinoless double beta decay, by providing an upper bound on the cut-off scale of the effective field theory for HNL masses above $\sim$100 MeV. We identify target regions for future laboratory probes of the $\nu$SMEFT parameter space that is bounded from above and below.

Time series analysis from gravitational-wave detectors often relies on the assumption that time chunks, or frequency bins, are uncorrelated. We discuss the validity of this approximation in the context of searches for stochastic gravitational-wave backgrounds. We examine the impact of averaging over time and frequency, a reduction technique commonly employed to minimize the computational expense of likelihood evaluations. We introduce an analytical tool based on Fisher information to quantify the error in parameter inference arising from ignoring these effects. Finally, we address the issue of locally stationary processes and optimal time chunking.

The fundamental properties of ice have always attracted a lot of interest due to omnipresence of ice in many different natural contexts. Since cubic ice recently become experimentally accessible from a low-density gas hydrate precursor [1, 2], it has been possible to measure its density as a function of temperature in the whole thermodynamic range of metastability. We found strong analogies with respect to the other ice I polytype, i.e., hexagonal ice Ih [3], including the presence of a negative thermal expansion behavior at low temperature. Based on these results, a new enthalpy calculation quantifies the metastable nature of the cubic form and, consequently its inaccessibility from a "normal" ice Ih precursor.

One of the promising dark matter (DM) candidates is a keV scale sterile neutrino. In the early universe the observed relic of the sterile neutrino DM is generated via the \textit{Dodelson-Widrow} mechanism. However, this production scenario is severely constraint by various astrophysical observations. Many non-standard interactions between active ($\nu_a$) and sterile ($\nu_s$) neutrino have been proposed to evade these astrophysical bounds. Here, we study sterile neutrino in the context of a mass-varying scenario by coupling both active and sterile neutrino to a scalar field. This novel mechanism opens up a new parameter space that generates the observed DM relic and alleviates the \textit{Hubble tension}. We find that the resulting parameter space can be fully probed by future X-ray missions.