Papers for Wednesday, May 20 2026

The moving lens effect is a secondary CMB anisotropy induced by the transverse motion of gravitational potentials. We develop a Fourier-space cross-spectrum estimator that retains the scale dependence of the signal, and apply it to the Atacama Cosmology Telescope (ACT) DR6 CMB temperature maps and luminous red galaxies from the DESI Legacy Imaging Surveys. Using the foreground-reduced ACT NILC map, we find strong evidence for a non-zero amplitude of the cross-correlation $b_{\rm ML} = 1.24 \pm 0.26$ ($4.8\sigma$) for the extended sample and $0.93 \pm 0.25$ ($3.7\sigma$) for the main sample, both consistent with the halo-model prediction for the moving lens signal. Our Fourier-based pipeline enforces separation of scales between the reconstructed velocities and the cross-correlation, which we show is essential for foreground mitigation. The residual foreground contamination is expected to be significantly smaller than the signal from both simulations and the multi-frequency analysis presented in this paper. No curl-mode test exceeds $2\sigma$, and the results are robust across analysis variants. They constitute the first detection of the moving lens effect and unlock access to transverse velocities, a new cosmological probe. When combined with the kinematic Sunyaev-Zel'dovich effect, this provides a path toward mapping the three-dimensional velocity field of the Universe, opening a new avenue for probing the growth of structure and gravity on large scales.

Patchy reionization couples the ionized-bubble morphology to the underlying density field, making the CMB Thomson optical depth sensitive to both the global ionization history and anisotropic fluctuations on the sky. Using the large-volume radiation-hydrodynamical Lumina simulation, we compute $\tau_{\rm CMB}$ in two ways: (i) from global volume- and mass-weighted ionization histories, and (ii) from explicit line-of-sight integrations through on-the-fly light cones. We find that the sightline-averaged optical depth in the light cone, $\langle \tau_{\rm LOS} \rangle = 0.0550$, exceeds the value inferred from a global volume-weighted history, $\tau_{{\rm CMB},V} = 0.0515$, by $\approx 7\%$. This enhancement is largely captured by the global mass-weighted prediction, $\tau_{{\rm CMB},m} = 0.0544$, indicating that precision comparisons to CMB optical-depth constraints should use mass-weighted electron fractions or explicit light-cone integration rather than volume-weighted ionized fractions alone. The excess optical depth accumulates primarily near $z_{\rm LOS} = 8.0^{+1.9}_{-1.3}$, where the combination of high physical density and strong ionization-field patchiness is greatest. The resulting $\tau_{\rm LOS}$ field is non-Gaussian and exhibits $\gtrsim 5\%$ sightline-to-sightline scatter, with fluctuations tracing rare early-ionized overdensities and large-scale structure. Coarse-graining experiments show that smoothing the ionization field on $\gtrsim 3 {\rm cMpc}$ scales suppresses the density-ionization correlation and biases $\tau_{\rm CMB}$ low relative to the resolved calculation. Finally, angular power spectra and real-space correlation functions decomposed into HII, HeII, and HeIII auto- and cross-contributions reveal scale-dependent departures from simple hydrogen-helium co-tracing and evolving characteristic scales with redshift.

Dark matter halos destined to fall into a more massive host differ from typical field halos of the same mass even before infall. In cosmological simulations, we find that the progenitor mass functions of these "future subhalos" are systematically shifted toward higher masses, with the shift growing as infall approaches. The bias takes a compact form within extended Press-Schechter theory: the collapse barrier is multiplied by a function $\beta(D/D_\mathrm{infall},a)$, where $D$ is the linear growth factor at scale factor $a$ and $D_\mathrm{infall}$ is the growth factor at infall. We find $\beta(x,a)=(1-x)^{1.20+0.14a}$ for the $M_{200\mathrm{c}}$ mass definition and $(1-x)^{1.20+0.05a}$ for $M_{200\mathrm{m}}$; the explicit scale-factor dependence captures the late-time influence of dark energy. One consequence is that halos shortly before infall are 10-15% more centrally concentrated than typical field halos of the same mass.

We study the formation and evolution of bulgeless galaxies within the Milky Way-Andromeda analogue sample of the TNG50 simulation. Through kinematic decomposition with Mordor, we identified bulgeless galaxies with a bulge-to-disc mass ratio of B/D<0.08, in line with the Bulgeless Evolution And the Rise of Discs (BEARD) survey and Milky Way constraints. We compared them to bulge-dominated galaxies (B/D>1). We find that 74% of bulgeless galaxies experience at least one major merger (stellar mass ratio 1:4) over their lifetime. Bulgeless galaxies form later ($z_{50}\sim 0.7$) than bulge-dominated counterparts ($z_{50}\sim1.2$). Bulgeless galaxies have lower-mass haloes and higher specific stellar angular momentum, compatible with Milky Way observations. However, specific star formation rates and hydrogen gas fractions are slightly higher than Milky Way observations. Our analysis of the redshift evolution of stellar components reveals that bulgeless galaxies have gradual disc growth with high thin disc-to-total mass ratios (D/T>0.5) since $z\sim 1$ and minimal bulge growth (B/T<0.1) since $z\sim1.5$. In contrast, bulge-dominated galaxies have earlier disc formation, which is disrupted, resulting in higher morphology evolution. Bulgeless galaxies are more likely to undergo gas-rich, coplanar, and corotating mergers, promoting disc survival, compared to bulge-dominated galaxies that encounter a broader spectrum of mergers. We also observed differences in galaxy structure between bulgeless and bulge-dominated galaxies without major mergers, suggesting the relevance of early gas accretion and alignment. Bulgeless galaxies have younger stellar populations and more extended star formation histories than bulge-dominated galaxies, which rapidly quench and have older stellar populations. These findings elucidate the distinct merger-driven and secular pathways that give rise to Milky Way galaxies.

Mars is the Solar System's canonical small, rocky planet that transitioned from early geologic activity and surface liquid water to a cold and arid planet with a thin, cold, CO$_2$-dominated atmosphere. The evolution of Mars, in the context of such planetary parameters as size, mass, atmosphere, insolation flux, magnetosphere, and impact history, harbor important diagnostics regarding the development and sustainability of habitable surface conditions. In this work, we synthesize how the study of Mars contributes to our understanding of exoplanet processes, such as volatile delivery and loss, photochemistry, climate evolution (including CO$_2$ condensation and atmospheric loss), obliquity forcing, planetary architecture, and the role of intrinsic magnetism. We also evaluate optimal methods and prospects for detecting and characterizing potential Mars analogs beyond the Solar System. We focus on relevant results from planetary missions (e.g., Mars Reconnaissance Orbiter, MAVEN, Mars Science Laboratory, Mars2020) and observational studies of exoplanet atmospheres with the James Webb Space telescope (JWST) and future facilities. Through the convergence of these parallel pathways of inquiry, we describe the primary science questions and suggested avenues for characterizing small rocky planets that lie at the edge of potentially habitable conditions.

The Hertzsprung-Russell diagram (HRD) is central to stellar astrophysics but has rarely been used to interpret planet formation. We extend the HRD concept to forming planets and study how solid and gas accretion, cooling/contraction, and migration shape luminosity-temperature tracks in different formation scenarios. We compute planetary interior structures throughout formation and evolution with the Bern model and, for the first time, couple it to radiation-hydrodynamical simulations to obtain a time-dependent accretion-shock heating efficiency, helping to address the cold-/hot-start ambiguity. Planetary HRDs exhibit three branches corresponding to successive phases: (i) an ascending branch during solid-dominated growth, strongly set by the size of accreted bodies (and thus the solid accretion rate) and by migration; for in-situ planetesimal accretion we find analytically $L \propto T^8$. (ii) A near-horizontal branch beginning at detachment when gas accretion becomes disk-limited and contraction accelerates; hot accretion, higher masses, and pebble accretion bend tracks upward. Increasing electron degeneracy after detachment lowers interior temperatures and stabilises radii. (iii) A descending branch where accretion ends and planets join constant-mass cooling tracks with weak radius evolution and $L \sim T^4$. Our tracks agree well with synthetic populations and are broadly consistent with directly imaged planets. Populating the short-lived early branches observationally will be difficult, and embedded accreting planets require models including accretion-shock emission and circumplanetary-disk reprocessing.

We present analysis and modeling of optical and infrared (IR) spectroscopy of the Type II supernova (SN II) 2024ggi obtained with ground-based instruments and the James Webb Space Telescope (JWST) at phases of ~265 - 400 days. The near- and mid-IR spectra reveal diverse iron-group emission-line morphologies, including double-peaked profiles in [Ni I] 3.119 and 11.998 $\mu$m, [Fe II] 1.644 and 17.931 $\mu$m, and [Co I] 12.255 $\mu$m, alongside Gaussian profiles in [Ni II] 1.939 $\mu$m, [Co II] 10.520 $\mu$m, and [Ni I] 7.505 and 11.304 $\mu$m. These differences imply both chemical inhomogeneity and aspherical ionization of inner ejecta, consistent with expectations from the $^{56}$Ni bubble effect. Modeling of double-peaked profiles supports an ejecta distribution with polar enhancements as large as ~7 for Ni/Co/Fe-rich material and ~2 for intermediate-mass elements. LTE estimates imply a stable Ni mass of $M_{\rm Ni}\approx1.3\times10^{-3}$ M$_{\odot}$, but electron densities near critical values indicate departures from LTE. Comparisons to non-LTE radiative transfer models favor a progenitor mass of ~12 - 15.2 M$_{\odot}$. We show that a simple mapping between elemental mass distribution and projected velocity reproduces line profiles produced in a CMFGEN radiative transfer calculation. We apply this property to 3-D neutrino-driven explosion simulations and predict Ni emission profiles for varying viewing angles. We find that only energetic 3-D explosion models of high-mass progenitors reproduce the observed extent of Ni mixing in SN 2024ggi, conflicting with progenitor masses inferred from radiative transfer models. These results demonstrate the utility of resolved nebular IR lines as direct probes of the 3-D distribution of explosively synthesized material in core-collapse SNe.

Astrometric monitoring of stars provides a promising method for discovery of low-mass planets around nearby Sun-like stars. The astronomical community has proposed several telescopes designed to perform high-precision astrometric observations. One limiting factor intrinsic to stars is the astrometric noise - or "jitter" - induced by surface stellar activity such as starspots and faculae. Despite previous estimates, the relative size of this signal has not been empirically measured from direct photometric observations. We analyse high-resolution images of the Sun to quantify the photometric centroid jitter across three narrow wavelength regions over nearly a decade, spanning high and low activity periods of the Solar cycle. We compare our findings to previous theoretical estimates. We scale this jitter to simulate how a Solar-twin would appear at various distances, establishing an astrometric noise floor below which detection is significantly complicated by stellar activity. We also introduce starspot simulations that augment our data. We find the typical astrometric jitter of the Sun at \(\lambda = 607.2 \pm 0.25\text{nm}\) to be \(0.342\mu \text{as pc}\), ranging between \(0.058\mu \text{as pc}\) and \(1.294\mu \text{as pc}\) for low and high activity periods, respectively. This is lower than the expected \(\approx 3\mu \text{as}\) astrometric signal that an Earth-like planet would produce around a Sun-like star, at 1 pc. Therefore, the astrometric noise floor imposed by intrinsic stellar activity sets a detection limit below one Earth but greater than Mars around Solar-analog stars, making instrument precision the limiting factor for Earth-like exoplanet searches.

We present an automated pipeline for the detection of EVPA rotations in blazars, integrating correction of the 180$^\circ$ ambiguity, Bayesian Blocks segmentation, and statistical validation. Applied to RoboPol monitoring data, the method identified 48 rotations across 25 sources, including multiple events in RBPLJ2232+1143, RBPLJ1751+0939, RBPLJ1800+7828, and RBPLJ2253+1608. The rotations span amplitudes from 90.8$^\circ$ to 359.7$^\circ$, durations between 7.0 and 111.3 days, and rotation rates averaging 5.0$^\circ$/day. Comparison with previous catalogs reveals systematic differences: Bayesian Blocks rotations are on average $\sim$10\% larger in amplitude, about twice as long in duration, and roughly two-thirds slower in angular velocity, reflecting systematic biases between adaptive binning and manual segmentation. In addition, we report 27 previously unreported rotations, including 11 from the final 2016--2017 season. A correlation analysis with contemporaneous Fermi--LAT $\gamma$-ray light curves shows that longer rotations tend to coincide with enhanced $\gamma$-ray activity, while rotation amplitude alone is not predictive of $\gamma$-ray brightness. Our pipeline minimizes subjective biases, expands the list of known EVPA rotations, and provides a reproducible framework for future multiwavelength studies of blazar jet dynamics and particle acceleration.

Accretion onto supermassive black holes powers the most luminous persistent sources in the Universe, the so-called active galactic nuclei, whose emission is characterized by two distinct spectral components: thermal optical/ultraviolet radiation from an optically thick accretion disk and a power-law X-ray tail from a corona located in the innermost regions of the accretion flow. Yet, how radiatively efficient accretion disks develop and couple to the hot corona remains poorly understood. Using six years of simultaneous ultraviolet and X-ray monitoring of the nearby active galaxy ESO 511-G030, we witness a dramatic evolution of the broadband spectral energy distribution, driven by an increase of the ultraviolet flux from the disk by more than an order of magnitude over a time scale of less than three years. The overall behavior is unlikely to track an uncovering event, and is instead compatible with a progressive recovery of the optically thick component of the accretion flow. At accretion rates higher than approximately one per cent of the Eddington limit, ultraviolet and X-ray data are tightly coupled and follow the well-defined, non-linear correlation between disk and corona found in the more luminous quasars. Below this threshold, the relation apparently breaks down, as expected in case of evaporation of the inner accretion disk into a geometrically thick, optically thin hot flow. This is a strong hint of an accretion-state transition analogous to those observed in stellar-mass black holes, and confirms the need for a paradigm change in the models of radiatively efficient accretion flows around supermassive black holes.

Electron-capture supernovae (ECSNe) may arise from ONeMg-core collapse in super-asymptotic giant branch (sAGB) stars near the low-mass core-collapse limit ($\approx\!8$--$10$\,\Msun). At early times, models predict that ECSNe resemble low-mass red supergiant iron-core-collapse SNe (FeCCSNe), making the two channels difficult to distinguish. Nebular spectroscopy, however, can reveal differences in ejecta composition. We present a systematic sample of nebular spectra of 19 low-luminosity Type IIP (LLIIP) SNe from the ZTF CLU survey, obtained 115$-$450\,d after explosion. Their low velocities expose narrow lines blended in brighter SNe, which we identify and model to constrain progenitor properties. We find a strong correlation between the FWHM of H\,\textsc{i}\,$\lambda$6563 and peak luminosity, showing that LLIIP SNe occupy the low-energy end of the core-collapse population, but no correlation with plateau duration, suggesting that envelope and core properties are not tightly linked. Only one SN reaches the extremely low H\,\textsc{i}\,$\lambda$6563 widths predicted for the weakest $\sim$9\,M$_\odot$ explosion models, implying that such low-energy events are intrinsically rare. Combining our sample with 118 literature nebular spectra of Type II SNe, we infer an IMF slope of $2.1\pm1.2$. We also introduce an `ECSN score'' based on the absence of He- and O-shell emission lines, and identify two plausible ECSN candidates, SN~2023bvj and SN~2024btj. However, neither shows the extremely narrow nebular lines predicted by current ECSN models. If ECSNe arise predominantly through the LLIIP channel, we infer an upper limit on the ECSN rate of $\lesssim (5$--$8)\times10^{2}\,\mathrm{Gpc^{-3}\,yr^{-1}}$, corresponding to a narrow sAGB progenitor mass window of $\Delta M_{\rm sAGB} \lesssim 0.02$--$0.06\,\mathrm{M_\odot}$.

The extent to which galaxy structure is shaped by environment beyond the local universe, once stellar mass is controlled, remains an open question in galaxy evolution. We address this challenge using an unprecedentedly large sample of $\sim$3 million galaxies from the Hyper Suprime-Cam Subaru Strategic Program spanning $0.3 \leq z < 0.7$ with $\log(M/M_{\odot}) \geq 8.9$. We correlate a mass-independent bulge-to-total ratio statistic with large-scale overdensity maps and cluster catalogs, propagating structural parameter posteriors through a Monte Carlo framework to robustly assess significance. We confirm with $>5\sigma$ confidence that galaxy structure depends on environment at fixed stellar mass, but this dependence is secondary to stellar mass and varies with redshift, mass, and environmental scale. At $z < 0.5$, we detect no significant structural correlation with large-scale overdensity, but cluster galaxies show statistically significant bulge enhancement compared to mass-matched field galaxies, indicating cluster-specific processes such as ram-pressure stripping and cumulative tidal interactions dominate structural transformation at these epochs. At $z \geq 0.5$, massive galaxies exhibit bulge-enhancement across both cluster- and large-scale environments, while lower-mass systems show enhancement only in cluster environments. This indicates that environmental mechanisms operate across broader spatial scales at earlier cosmic epochs, and enhanced merger rates, group preprocessing, and cosmic web stripping augment cluster-specific processes. By separating into star-forming and quiescent subsamples, we find nearly flat trends within each, demonstrating that the observed environmental effects arise from coupled morphological and star formation transformations. These results collectively reveal the multi-scale, epoch-dependent nature of environmental effects on galaxy structure.

Flavor conversion can affect the neutrino-driven delayed explosion mechanism of collapsing massive stars, altering the efficiency of shock revival. We perform core-collapse supernova simulations in spherical symmetry for a set of progenitors with masses of $9.75\, M_\odot$, $11\, M_\odot$, $16.5\, M_\odot$, $28\, M_\odot$, $40\, M_\odot$, and $60\, M_\odot$, accounting for a mixing-length treatment for convection. Flavor conversion is modeled assuming instantaneous flavor equipartition below a critical baryon density, while conserving the lepton number. Regardless of the progenitor compactness, its mass, or the nuclear equation of state, we find that flavor conversion can increase heating (cooling) and enhance (hinder) the supernova explosion, if triggered near the gain (neutrino decoupling) region. Our findings suggest that the interplay among the region of the supernova core where flavor conversion occurs, the progenitor properties, and the nuclear equation of state is crucial in determining the fate of explosion and the properties of the compact remnant.

The association of broad-lined Type Ic supernovae with long-duration gamma-ray bursts (GRBs) has been known for 28 years. However, only about seventy gamma-ray burst supernovae (GRB-SNe) have been identified, of which only half have spectroscopic classifications. At $z=0.153$, GRB 260310A is the 12th spectroscopically confirmed GRB-SN discovered within 1 Gpc, offering a critical opportunity to follow one of these rare supernovae in detail. We present optical to near-infrared imaging and spectroscopy of GRB 260310A and SN 2026fgk out to 65 d after discovery. The optical afterglow is among the brightest ever observed from a GRB. Spectra obtained more than two weeks after the explosion reveal broad absorption features that securely identify SN 2026fgk as a Type Ic-BL supernova. Modeling of the multi-wavelength ($grizJK_s$) lightcurve shows that the supernova is approximately half the luminosity ($k_\textrm{98bw}=0.4-0.6$) of the canonical GRB-SN 1998bw. We derive a nickel mass of $M_\textrm{Ni}=0.4-0.5$ $M_\odot$ with a total ejected mass of $M_\textrm{ej}\approx4-6 $ $M_\odot$ and kinetic energy $E_\textrm{K}=(3-8)\times10^{51}$ erg. The GRB exploded at an extremely large offset of 15 kpc from its host galaxy. Long-slit spectra reveal a ``bridge'' of nebular emission extending along the galaxy's disk to the GRB location, which has a sub-solar metallicity ($\sim$\,$0.4Z_\odot$), compared to a near solar metallicity for the host galaxy. This indicates that the large offset arises from the galaxy's extended light profile rather than an isolated environment.

$Aims$. The existence of massive disc galaxies with little or no bulge challenges conventional $\Lambda$ cold dark matter model, which typically favours dynamically hot central structures due to early collapse and mergers. The study of these bulgeless disc galaxies is the aim of the Bulgeless Evolution And the Rise of Discs (BEARD) survey, as they offer a unique opportunity to investigate the link between galaxy morphology and the properties of their satellite systems. $Methods$. Using the high-resolution cosmological hydrodynamical simulation TNG50-1, we studied the satellite populations of 135 bulgeless galaxies. We compared their satellite properties to those of a bulge-dominated control sample with matched stellar masses. Our analysis focuses on satellite abundance, luminosity functions, spatial distribution, orbital alignment, and infall histories. $Results$. We find that satellite abundance is largely independent of host galaxy morphology. However, satellites around bulgeless galaxies exhibit luminosity functions with a steeper faint-end slope, are more centrally concentrated, and show stronger orbital alignment with the host disc plane. The orbital alignment originates from coherent post-infall dynamical evolution that depends on host galaxy morphology. The infall of more massive satellites can additionally perturb this process, contributing to a weakening or temporary stalling of the secular alignment. $Conclusions$. Due to the co-evolution of the host galaxy and the satellite system, the morphology of the central galaxy leaves a clear imprint on its satellite system. Bulgeless galaxies tend to have dynamically colder, more aligned, and more centrally concentrated satellite populations. These trends reflect a more quiet merger history and support the use of satellite properties as tracers of host galaxy formation pathways.

In the standard $\Lambda$ cold dark matter ($\Lambda$CDM) cosmology, the existence of massive pure-disc galaxies remains challenging within the hierarchical framework and is key to understanding the evolutionary history of Milky Way-like systems. In this work, we investigate the physical origin of the scatter in the stellar mass-size relation of massive spiral galaxies, with a particular focus on bulgeless systems. We analyse 22 nearby bulgeless galaxies from the Bulgeless Evolution And the Rise of Discs (BEARD) survey using deep $g$- and $r$-band imaging obtained with the 2.5 m Isaac Newton Telescope Wide Field Camera. We derive surface-brightness, colour, and stellar-mass-density radial profiles to measure $R_1$, the radius where $\Sigma_* = 1\,\mathrm{M}_\odot\,\mathrm{pc}^{-2}$, adopted here as a physically motivated size proxy. Point spread function (PSF) effects are corrected through star subtraction and wavelet deconvolution. BEARD bulgeless galaxies follow the tight stellar mass-$R_1$ relation defined in previous studies, with a similar scatter of $\sim 0.1$ dex. Using galaxies from the IllustrisTNG50 simulation, we find that the scatter is linked to morphology, with bulgeless, BEARD-like analogues and bulge-dominated galaxies tracing the upper and lower envelopes of the relation, respectively. This trend correlates with the specific central stellar mass density, $\Sigma^{\mathrm{spec}}_{1,\mathrm{kpc}}$, suggesting that differences between observations and simulations reflect the broader central-density range spanned by BEARD bulgeless galaxies. A deeper analysis of the physical driver of this morphological segregation reveals that the scatter in the mass-size relation is also related to the spatial configuration of merger events, rather than their frequency, with bulgeless systems tending to inhabit halos with a slightly higher spin. (abridged)

We investigate the role of ejective and preventive feedback in $\mathrm{\sim10^{10}-10^{11}\,M_\odot}$ dwarf halos using cosmological zoom-in simulations. These simulations use adaptive mesh refinement to capture high-specific-energy outflows, together with an implementation of discrete supernovae (SNe). We show that episodic, SNe-driven shock heating sustains the circumgalactic medium (CGM) at $\mathrm{\sim T_{vir}}$. This process also increases the ratio $\mathrm{t_{cool}/t_{ff} > 10}$ in the outer CGM and intergalactic medium (IGM), placing the gas in a radiatively stable regime. Hot outflows ($\mathrm{\gtrsim10^5\, K}$) dominate the energy budget, and their high specific energy allows them to traverse the CGM, escape the halo, and heat the IGM. In contrast, warm outflows ($\mathrm{\lesssim10^5\, K}$) dominate the mass budget and are largely recycled back into the interstellar medium (ISM), where they fuel future star formation. We identify a gradual transition at $\mathrm{\sim 5\, Gyr}$ that marks a shift in the balance between ejective and preventive feedback. At early times ($\mathrm{< 5\, Gyr}$), although the CGM cooling rate dominates for a larger fraction of time, the infrequent yet powerful SNe energy injection into the CGM is able to quickly dominate the cumulative energy balance. These outflows and their high specific energy are able to 'sweep' up mass in the CGM and IGM. At late times ($\mathrm{> 5\, Gyr}$), the CGM baryon fraction is only $\mathrm{\sim0.1}$, leading to a transition toward a preventive feedback mode in which SNe maintain $\mathrm{t_{cool}/t_{ff} > 10}$ and prevent $\mathrm{\sim75\%}$ of the expected baryon accretion rate.

The chemical composition of the intracluster medium (ICM) provides key insights into the enrichment history of galaxy clusters. However, high-resolution abundance measurements with X-ray microcalorimeters remain available for only a few systems. While most cool-core clusters exhibit near-Solar elemental abundance ratios relative to Fe, previous studies of the Virgo cluster suggested super-Solar ratios in its core. We investigate the chemical properties of the Virgo cluster core using XRISM/Resolve observations, focusing on precise measurements of Si, S, Ar, Ca, Cr, Fe, and Ni abundances. We aim to determine whether Virgo displays abundance patterns distinct from other nearby cool-core clusters and to explore the origin of any differences. We analysed XRISM/Resolve spectra in four regions of the Virgo core (center, east, northwest, and southwest) in the 1.7-11 keV band. Single-temperature, multi-temperature, and multi-abundance models were applied to characterize the thermal structure and derive elemental abundances. The resulting abundance ratios were compared between pointings, with previous studies of Virgo/M87, and with recent XRISM measurements of other clusters. All four pointings exhibit systematically super-Solar X/Fe ratios, although the northwest region shows values closer to Solar. Multi-abundance modeling of the eastern and southwestern regions reveals that cool, metal-rich gas uplifted by the AGN coexists with a hotter, more chemically homogeneous ambient ICM. The super-Solar ratios are robust against variations in bandpass and temperature structure. We interpret these enhanced ratios as reflecting the enrichment history of the old stellar population in M87 combined with a limited cold gas reservoir. The Virgo core exhibits spatially resolved chemical enrichment that differs from the abundance patterns observed in other cool-core clusters.

We measure the projected rotational velocities ($v \sin i$) of the solar-like blue straggler stars (BSSs) in the old ($\geq4$ Gyr) open clusters M67, NGC 188, and NGC 6791. We find that the BSS rotation distribution shows a Kraft break similar to that found in the field. The main-sequence progenitors of these BSSs were cooler than the Kraft break and have spun down by their age. The binary interactions that create BSSs are expected to spin up these progenitors, so current BSS rotation rates are due to transformation and any subsequent spin-down. We observe that BSSs hotter than the Kraft break are rapidly rotating, showing that binary evolution spins up these, and likely all, BSSs to initial rotational periods below two days, still below critical velocity. BSSs below the Kraft break currently have slow rotation rates, and those within the Kraft break have a mixture of rotation rates suggesting rotational transition. This dependence of rotation on effective temperature indicates that BSS envelopes behave like those of single stars, becoming convective and generating magnetic fields at the same temperatures. For globular cluster BSSs with [Fe/H]$\sim-1.5$, we find evidence of a BSS rotation transition region that is 100-250 K hotter than at solar metallicity. We find the $v \sin i$ distributions of BSSs in open clusters have similar characteristics to both high- and low-density globular clusters, indicating the density of environment is not the only factor that can determine rotational distributions. We suggest that velocity dispersion plays an important role.

The Galactic halo preserves a record of the Milky Way's earliest assembly and contains both in-situ stars and stars accreted from dwarf galaxies. Possible planets around these stars, therefore, probe formation in ancient, metal-poor environments, including systems of extragalactic origin. We present a search for short-period transiting planets around kinematically selected halo dwarfs using Gaia DR3 and TESS, focusing on planets with periods of $1 < P < 10$ days. We identify two hot-Jupiter (HJ) candidates, one in the in-situ and one in the accreted halo, although the latter is highly grazing and excluded from the occurrence analysis. The accreted candidate, if confirmed, would orbit the most metal-poor HJ host known ([Fe/H] $\approx -1$). Using injection--recovery tests and automated vetting, we constrain occurrence in the full halo, in-situ, and accreted samples. In the HJ regime ($8\,R_\oplus < R_{\rm p} < 22\,R_\oplus$, $1\,\text{day}\ < P < 10$ days), the non-grazing candidate implies an overall halo occurrence rate of $0.13^{+0.12}_{-0.07}\%$ if planetary, while the absence of confirmed detections gives a corresponding $1\sigma$ upper limit of $<0.14\%$. For the in-situ halo, we infer $0.17^{+0.17}_{-0.10}\%$ (or $<0.19\%$ assuming no detections), while for the accreted halo we derive an upper limit of $<0.56\%$. These rates lie well below the corresponding short-period giant-planet occurrence measured in the Galactic disc. A forward model assuming Kepler-like occurrence also predicts $10 \pm 3$ detections compared with at most one observed. We find no significant occurrence difference between the in-situ and accreted halo populations, strengthening the evidence that close-in giant planets are rare across the old, metal-poor halo.

We report the first detection of stratospheric propadiene (CH$_2$CCH$_2$) and propene (C$_3$H$_6$) at Jupiter's mid-to-high northern latitudes using IRTF-TEXES measurements recorded on March 5-6, 2025. Using radiative transfer software to quantitatively test for the presence of propadiene and propene, we report a $>$12-$\sigma$ detection of propadiene and a $>$17-$\sigma$ detection of propene inside Jupiter's northern auroral region (henceforth 'NAR'), where the species are most concentrated. For example, at 62$^\circ$N inside Jupiter's NAR, we derive a 1-mbar propadiene abundance of 2.0 $\pm$ 0.2 ppbv, which is 40 $\pm$ 3 higher than abundances predicted by the Moses & Poppe (2017) photochemical model (henceforth 'MP17'), and significantly higher than the 1.2-ppbv upper limit abundance derived at 42$^\circ$N (the lowest latitude sampled by the observations). Similarly, we derive a 1-mbar propene abundance of 8.1 $\pm$ 0.5 ppbv at 62$^\circ$N inside Jupiter's NAR, which is 28 $\pm$ 2 higher than the MP17 predicted abundance and higher than the 6-ppbv 1-mbar upper limit abundance derived at 42$^\circ$N. The fact that propadiene and propene are most enriched inside Jupiter's NAR strongly suggests that perturbations to the chemistry by auroral-related heating and exogenous ions/electrons are responsible for their significant enrichment. Spectral features of propane (C$_3$H$_8$) were not detected at any of the locations sampled by the data: 3-$\sigma$ upper limits of 10 ppbv were derived at the 10-mbar level at 62$^\circ$N inside Jupiter's NAR. The non-detection of propane could, in part, be explained by the vertical sensitivity of its spectral features to deeper pressures, where there is negligible auroral-related heating. The results of this work advocate for development of ion-neutral chemistry models of Jupiter's polar stratosphere.

Random telegraph noise (RTN) is a major contributor to read noise in many CMOS image sensors considered for astronomical use. While scientific CMOS image sensors deliver lower read noise than traditional charge-coupled devices, mitigating RTN would widen this gap and enable more precise photometry when using the fast readout rates achievable by CMOS image sensors. We report the levels of RTN in three CMOS image sensors used in astronomical instruments: the Sony IMX455, Gpixel GSENSE400, and Fairchild Imaging HWK4123. For the IMX455 in a high gain mode, RTN is the dominant source of pixels with high read noise and increases the overall read noise floor by >20%. RTN is present in the GSENSE400 and HWK4123 but to smaller effects. We compare two strategies for RTN mitigation: masking pixels exhibiting RTN or using a new algorithm for correcting RTN jumps. For faint (< 3 e-/pix/frame) observations of a stellar field with the IMX455, both masking and our algorithm improved the signal-to-noise ratio (SNR) of light curves by >5% on average. Larger improvements were achieved for sources falling on multiple RTN pixels. Our algorithm outperforms masking, especially when the point spread function is undersampled, masked pixels are near the source center, or read noise and shot noise are comparable. In such cases, masking may even deteriorate photometric precision. In other cases, masking remains an effective RTN mitigation technique. We have made available our software for identifying RTN pixels, parametrizing their bias level distributions, and applying our correction algorithm.

We present IceCube's latest release of muon track data for neutrino point-source searches, extending the previously published 10-year dataset to cover 14 years of observations (April 6, 2008 - May 23, 2022). This release features an updated event selection and improved detector calibration for data recorded after June 1, 2010. The release also includes binned instrument response functions and effective areas, enabling the community to perform sensitive searches for steady and transient neutrino sources. We report on key science results obtained with this dataset using internal IceCube analysis tools and compare them to those derived from analyses based on the binned response functions included in this public release. To facilitate reproducible research, we provide benchmark results obtained using this data release and publicly available software. This release represents IceCube's most sensitive and comprehensive publicly available all-sky muon track dataset to date and should be preferred over previous releases.

Magnetohydrodynamics (MHD) couples the Navier--Stokes and Maxwell equations into a nonlinear system of partial differential equations governing stellar interiors, astrophysical jets, fusion plasmas, and space weather. Numerical advances, including finite-volume Godunov schemes, constrained-transport algorithms, high-order spectral-element and discontinuous-Galerkin discretisations, and adaptive mesh refinement, have made MHD a predictive tool for solar eruptions, tokamak confinement, and magnetised turbulence. A fundamental barrier nevertheless remains. In three-dimensional MHD turbulence, the degrees of freedom required to resolve all active scales grow as $\mathcal{O}(\mathrm{Re}^{9/4})$ or faster, where $\mathrm{Re}$ is the Reynolds number. Direct numerical simulation is therefore intractable at astrophysical and fusion-relevant parameters, particularly when the Lundquist number $S$ exceeds $10^{10}$ and both viscous and resistive dissipation ranges must be resolved. Kinetic closures, radiation transport, and uncertainty quantification further increase the cost. This chapter examines how AI may help bridge this gap. We review physics-informed neural networks, Fourier neural operators and physics-informed neural operators, which learn solution operators across families of MHD problems; and hybrid operator-diffusion frameworks that combine deterministic surrogates with score-based generative models to recover broadband turbulent spectra. These developments are set within the wider landscape of exascale high-order solvers, GPU acceleration, task-based parallelism, data-driven sub-grid closures, and prospective quantum algorithms for implicit linear systems in resistive MHD. The central claim is that physics-informed AI, integrated with conventional solvers and trained on leadership-scale simulations, offers a credible route to regimes beyond the reach of classical discretisation alone.

The Near-Earth Object Wide-field Infrared Explorer (NEOWISE) mission provides a decade of all-sky time-series data at 3.4 and 4.6um and an unprecedented opportunity for the discovery and characterization of variable objects. This paper presents VarWISE, a catalog of infrared-variable objects discovered within the NEOWISE single-exposure data. We employ unique methodologies, including the spatial clustering of apparitions and the adoption of novel machine learning-based variable detection (VARnet) and classification (XGBoost) to identify and characterize significant variability. The catalog includes a prediction of variable object type and best-fit period values for each object, if its variations are cyclical, along with other calculated parameters to characterize the nature of the variability. The VarWISE Pure Catalog, containing only variables of highest confidence, has 457,080 objects, 49.81% of which are new discoveries; the VarWISE Extended Catalog, containing all sources, has 1,918,082 objects, 82.02% of which are new. We discuss caveats for each variable type and highlight a few new objects found during a quick perusal of the catalogs' contents.

The $98^\circ$ obliquity of Uranus is widely attributed to a giant impact that ejected material and formed a debris disk, which subsequently coalesced into its regular satellites. Previous Smoothed Particle Hydrodynamics (SPH) studies have yielded inconsistent disk compositions, a discrepancy often linked to the variety of numerical and physical modeling assumptions. We address this by presenting SPH simulations that systematically test three distinct EOS models alongside two SPH schemes (standard SPH, and the enhanced density-independent SPH). We utilized a $3M_{\oplus}$ impactor and explored a range of impact parameters which are capable of reproducing Uranus's current spin state. Our primary finding is that for impacts capable of reproducing Uranus's current rotation, the choice of EOS or SPH scheme barely affects macroscopic features such as the post-impact rotation period, disk mass, or disk size; these properties are primarily controlled by the impact's angular momentum. In contrast, the disk's rock fraction is highly EOS-dependent. Our results clarify that while disk mass and size are robust outcomes, the final disk composition is highly model-dependent. Therefore, accurate EOS modeling, integrated with detailed disk evolution studies, is essential to definitively validate the giant impact scenario for Uranus.

Two-dimensional relativistic particle-in-cell (PIC) simulations of radiative magnetic reconnection in pair plasmas with multiple interacting current sheets are carried out to mimic the dynamics in high-energy astrophysical environments, such as particle acceleration regions in pulsar wind nebulae and relativistic outflows, where the magnetic field is expected to reverse polarity multiple times. Initially, due to reconnection within each isolated sheet, particles are accelerated and synchrotron emission beyond the burn-off limit is confirmed, even if the particle distribution function shows steep slopes. After this phase, plasmoids lead to cross-sheet interactions and merging, with new current sheets formed. In this regime a Kolmogorov-like spectrum for the magnetic energy develops over a couple of decades, followed by a dissipation range starting around 5~$d_e$ (electron inertial lengths), showing that multi-sheet reconnection evolves nonlinearly into well-developed turbulence. This phase provides secondary acceleration and further cooling by synchrotron emission, with intermittent radiative bursts. We show that high energy accelerated particles by the primary current sheets are further energized during the turbulent phase, while the distribution of the most energetic particles remains steep.

We present a cloud-based tool that uses drones and machine learning to help recover instrumentally observed meteorite falls. We showcase a collection of improvements made upon previous iterations of our system, as well as detail the successes and limitations of this technique when applied to observed meteorite falls in South and Western Australia. This tool is available to the meteoritics research community upon request at this https URL.

We perform model-independent distance measurements on four Type Ia supernovae (SNe Ia) compilations (Pantheon, Pantheon+, DES-Dovekie, Union3) and compress each dataset into the values of $\log r_p(z)$ at eleven redshift knots, where $r_p(z)$ is a rescaled comoving distance. These Gaussian distributed compressed values, together with their full covariance, completely capture the distance-redshift relation information from each dataset. We demonstrate this by using these to perform an Markov Chain Monte Carlo (MCMC) likelihood analysis to infer cosmological parameters in flat $\Lambda$CDM, flat $w_0 w_a$CDM, and a non-parametric reconstruction of the dark-energy density $X(z) \equiv \rho_{\rm DE}(z)/\rho_{\rm DE}(0)$. Across all datasets and flux-averaging configurations and all three cosmological models, the resulting parameter contours and figures of merit reproduce the corresponding full distance-modulus analyses using the original SNe Ia data sets within the statistical sampling noise of the chains, demonstrating that the eleven $\log r_p$ data points are an operationally lossless compression of the cosmological information in the dataset. Our SN Ia data compression enables an analytic analysis that completes in $O(10^{-2})$ s per dataset and reduces the downstream cosmological MCMC to the fast evaluation of an $11$-dimensional Gaussian likelihood, with a per-step cost set by the number of knots and independent of the SNe Ia sample size. Our methodology will benefit the data analysis of future surveys from Euclid, Roman, and LSST, which will deliver SNe Ia samples one to three orders of magnitude larger than current ones.

The transient X-ray pulsar GRO J1008-57 was previously found to exhibit Type I outbursts occurring at stable orbital phases before its first observed Type II outburst in 2012. In this work, we extend the study to investigate the phase evolution after several Type II outbursts using long-term Swift/BAT and MAXI/GSC observations. Our results reveal that the orbital phases of Type I outbursts follow a step-like evolution: they remain largely stable over many orbital periods but undergo abrupt, small-amplitude jumps coincident with each Type II outburst. Such a step-like behavior is difficult to explain with the commonly proposed mechanisms involving a highly eccentric or precessing disk around the Be star. The energetics of Type I X-ray outbursts show a systematic increase before Type II outbursts, followed by a rapid decline and a subsequent gradual recovery. This behavior suggests cycles of disk depletion and reconstruction driven by Type II outbursts. Considering the small amplitude of each phase jump, we propose that this step-like phase evolution may be related to the long orbital period of GRO J1008-57, implying infrequent neutron star-disk interactions. After disk depletion by Type II outbursts, the disk around the Be star has enough time to rebuild its density and restore a geometric structure similar to its pre-Type II outburst state. Consequently, the orbital phases of subsequent Type I outbursts not only change very slightly but can also remain stable over many orbital periods until the next Type II-driven disk reconfiguration, yielding the observed step-like evolution.

While dark matter (DM) makes up roughly 80% of the total matter in the Universe, its microscopic properties remain one of the biggest questions in Cosmology today. Fortunately, those properties dictate the distribution and form of macro-scale gravitational structures in the universe, allowing for indirect studies which can distinguish between competing particle models. One such avenue for this research is via strong gravitational lensing systems, where dark halos in the lens substructure and along the line of sight perturb image positions and flux. However, the current population of sources suitable for this analysis is limited, especially at radio wavelengths where astrometric perturbations are observable. I will first discuss which properties of lens systems make them especially useful for DM constraints and examine the minimum amount of information necessary for such an experiment. Then, I present the results of two successful searches for new radio lenses in existing radio and optical surveys, utilizing a new method to expand the potential follow-up population for dark matter studies in the future. I conclude with a discussion of the completeness of this population.

We present a broadband spectro-polarimetric analysis of 1,565 polarized radio sources from the LOFAR Two-Metre Sky Survey Data Release 2 (LoTSS-DR2) RM Grid catalogue. This study uses frequency-dependent Stokes Q and U spectra across the 120-168 MHz LOFAR HBA band to investigate their polarization properties. The polarization behaviour of each source is modelled with multi-component Faraday depolarization models to investigate the magneto-ionic environments responsible for low-frequency depolarization. Significant Faraday complexity is observed throughout the sample, with 43.2% of sources requiring two or three Faraday components. External Faraday dispersion dominates the depolarization behaviour, with 54.1% of sources classified as external-screen dominated and approximately 60% showing statistically significant evidence for turbulent external Faraday-active media, while only 10.3% are consistent with pure internal differential Faraday rotation. The intrinsic polarization angle and RM separations between fitted components are generally small, suggesting that many RM components trace physically related emission regions embedded within common magneto-ionic environments. A weak but statistically significant anti-correlation is detected between the polarization spectral index, $\beta$, and weighted Faraday dispersion, $\sigma_{\rm RM,wtd}$, for two-component systems, whereas one- and three-component populations show no significant trend. The rest-frame Faraday dispersion, $\sigma_{\rm RM,rest}$, exhibits significant positive correlations with redshift for the external-screen dominated and mixed depolarization populations, even after controlling for radio luminosity, indicating increasingly turbulent or strongly magnetized environments surrounding radio AGN at earlier cosmic epochs.

The Laniakea Supercluster is the closest large scale structure of galaxies. Is such a structure expected in the standard cold dark matter model of cosmology? This would be a relatively simple question to answer, were it not for the fact that the Zone of Avoidance (ZOA) runs right through it. Recent improvements to this paucity of data in the innermost ZOA can be made from systematic 21 cm surveys using the MeerKAT telescope (e.g. Kraan-Korteweg et al. 2024), and implementing these HI-redshifts as an extension to the CosmicFlows4 database for reconstruction (Hollinger et al. 2026). In this paper we test the assumption that for the purpose of reconstruction, additional HI detected galaxies without peculiar velocity determinations could be placed at their Hubble distances. We present infrared photometry of 163 of these in HI detected MeerKAT ZOA galaxies, in addition to 2MASS Extended Sources in the ZOA to determine their peculiar velocities. Averaging these peculiar velocities into redshift bins, we find that peculiar velocity corrections in the Laniakea Supercluster ZoA region are not prohibitively large, and that one can proceed with its reconstruction using the copious redshift data now available.

The outer regions of AGN disks have temperatures similar to those of circumstellar disks, permitting dust condensation. Therefore, planet formation and growth could be active in these dust tori through similar mechanisms. We aim at quantifying the parameter space for the occurrence of streaming instability, and its outcomes in terms of the masses of the objects formed, their total number, and their continued growth via pebble accretion. We use a a recently proposed disk model with strong magnetization to keep the disk gravitationally stable. We find that the dust grain sizes required for streaming instability are easily attained through coagulation; the dust filaments it produces can contain solar masses, collapsing into tens of millions of planetesimals ranging from Earth to super-Jupiter masses. These planets are usually born in the 3D Bondi regime of pebble accretion, and have mass-doubling times from 10^3 to 10^7 yrs, though 3D Hill and geometric accretion are also realized. Gas accretion occurs concurrently, and crossover mass can be attained while still in the planetary mass range. As a result, vigorous accretion can occur, leading to objects with stellar masses -- defining a core accretion channel for star formation. The pebble isolation mass is beyond the hydrogen burning limit, so accretion is limited by stellar feedback instead of gap carving. Our model also predicts a population of exotic objects directly formed above the hydrogen burning limit, yet of pure dust. Our approximate model suggests that AGN dust tori host the largest populations of planets in the universe.

Galaxy assembly bias (GAB) is the dependence of galaxy clustering on secondary properties beyond halo mass. In this work, we study the connections between GAB and baryonic processes using the Galacticus semi-analytic model (SAM) for galaxy formation and evolution applied to the UNIT simulation. By generating hundreds of galaxy mocks with varying parameters governing gas cooling, star formation, stellar feedback, and AGN feedback, we employ a shuffling method to quantify the GAB signal and compare the contributions of halo concentration and local environment to GAB. Using the Random Forest algorithm, we evaluate the importance of different baryonic processes for GAB. We find that for stellar-mass-selected galaxies, the dominant baryonic processes are gas cooling and stellar feedback, and the result does not change significantly with the number density; for SFR-selected galaxies, the primary process shifts from star formation to gas cooling as the number density increases. These results establish a direct and quantitative link between baryonic physics and GAB, which can provide guidance for empirical GAB parameterizations in upcoming and future galaxy surveys.

The crystallinity of water ice not only records the thermal history experienced by an astronomical body, but also affects the composition of forming planets by controlling the trapping of volatile materials in amorphous ice and their subsequent transport. An additional structure within the 3~$\rm \mu m$ water-ice absorption band, known as the Fresnel feature, may serve as a diagnostic of ice crystallinity. Recent observations with the James Webb Space Telescope have detected a Fresnel peak in a debris disk and in Trans-Neptunian Objects (TNOs). Here, we propose a portable expression that translates the observed Fresnel peak strength into the degree of crystallinity of icy grains in debris disks. Our formula targets scattered light at around 90$^{\circ}$ angles, which are easily accessible for spatially resolved debris disks regardless of the inclination angle. Applying this expression, we derive the degree of crystallinity of a debris disk around HD 181327 to be 10-20%. We also study the Fresnel feature in protoplanetary disks and find that it is generally weaker than in debris disks even for the same crystallinity. We then analyzed a scattered light spectrum of the protoplanetary disk around d216-0939, which shows a weak crystalline feature, and inferred a crystallinity of $\sim$50%. We conclude that the Fresnel feature is a reliable observational tracer for ice crystallinity, and future near-IR spectroscopic observations will be crucial to elucidate the crystalline ice evolution.

Gravitational microlensing is a powerful technique for constraining the abundance of dark matter in asteroid mass to supermassive primordial black holes at masses of $-11 \lesssim \log M/\mathrm{M}_\odot \lesssim 5$. In this work, we introduce a new high-cadence stellar microlensing survey in the Large Magellanic Cloud, AMPM. The primary goal of AMPM is to place constraints in the asteroid-to-planetary-mass regime of primordial black hole dark matter. We present the five nights of survey data, the microlensing detection pipeline, and the microlensing efficiency of AMPM. We explore the impact of the stellar distribution in the Large Magellanic Cloud on the microlensing detection efficiency and conduct a detailed analysis of second-order microlensing effects and the impact on the primordial black hole dark matter constraints. Our findings indicate that these second-order effects shift the maximum sensitivity of AMPM toward the lunar-mass black hole regime at $10^{-8} - 10^{-6} \, M_{\odot}$. From the five nights of data, we detect a single microlensing candidate and find that AMPM can constrain at the 95\% C.L up to 30\% of the Galactic primordial black hole dark matter distribution.

Metal-poor stars are crucially important for understanding the early Galaxy, first stars, and the Universe. In this series of papers, we present a homogeneous non-local thermodynamic equilibrium (NLTE) abundances analysis of 12 elements for 103 very/extremely metal-poor (VMP/EMP) stars with metallicity down to $-4.3$ dex. The sample was selected from the LAMOST survey and observed by the high-resolution spectroscopy of Subaru. In this paper, we present the NLTE abundances and evolution of lithium in these stars. We report different lithium behaviors corresponding to different evolutionary stages and their signatures: 1) The Spite Plateau shows a slightly positive slope, indicating increasing lithium abundance with increasing metallicity. Most significantly, it appears to extend to lower metallicities as previously suggested, calling into question the reality of the so-called 'meltdown' at low metallicity; 2) We confirm a lithium plateau for lower red giant branch (LRGB) stars with A(Li) $= 1.13$ dex in our sample, while lithium abundance drops rapidly to A(Li)$<0.5$ as stars continue to evolve to higher stage. 3) We identify four Li-rich stars in our sample across different evolutionary stages, showing complex and multiple lithium production mechanisms in VMP/EMP stars. These findings suggest that early Galactic lithium enrichment results from a complex interplay between depletion and production processes.

Standard CMB analysis assumes a direct deterministic mapping between the multipole probed by the CMB $\ell$ and the primordial wavenumber $k$. Since the recombination era has a finite duration, this mapping is probabilistic by construction. We elevate the power spectrum of the primordial perturbations to a probability distribution caused by the finite duration of the recombination era. We show that a finite recombination width introduces a Gaussian smoothing scale in $\ln k$ with $\sigma_{\ln k} \sim \sigma_\eta / D_*$, leading to a probabilistic mapping from multipoles to inflationary e-folds. This effect is zero in standard power-law inflationary scenarios, but it may become relevant for scenarios with exotic oscillating features of the primordial power spectrum, which will be probed by the future CMB experiments. The observed effective power spectrum is the true primordial spectrum blurred by the uncertainty in scale reconstruction, which is mathematically identical to a Bayesian marginalization over a latent variable, and thus there is a propagation of the measurement error in the independent variable, which is another more formal way to view the smoothing effect. Our results indicate that the smoothing has quantifiable effects on the spectral index and its running, but more importantly the difference between the TT and EE inferred spectral indices, $n_s^{TT}-n_s^{EE}$, is non-trivial, in contrast to standard inflation without smoothing, and might become observable by future cosmic microwave background experiments. Any tension in $n_s^{TT}-n_s^{EE}$ could indicate oscillations in the primordial spectrum and the effects of the power spectrum smoothing. Finally, a minimal Fisher matrix analysis is performed to investigate the observability prospects of the smoothing effect.

Primordial Black Holes (PBH) are hypothesised to form during inflation and have long been considered a candidate for compact dark matter. Gravitational microlensing is known as a productive method for exoplanet discovery and characterisation, but also provides an experimental avenue to constrain the PBH abundance in the mass regime from $\sim 10^{-11}\ M_{\odot}$ to $\sim 10^5\ M_{\odot}$. We performed a high-cadence, optical microlensing survey with DECam over five nights towards the Large Magellanic Cloud, sensitive to microlensing timescales from minutes to days. Here, we report the discovery of an hour-long microlensing event. An optical depth probabilistic analysis indicates that the lensing object, which we refer to as Phoebe, is 5 orders of magnitude more likely to be part of the Milky Way's dark matter halo than part of the stellar content of the Milky Way and Large Magellanic Cloud. No matter the location of Phoebe, it is among the fastest and lowest mass microlensing signals ever detected, with an Einstein timescale of approximately 60 minutes. Using Bayesian modelling, we interpret Phoebe as a PBH with mass $0.032^{+0.227}_{-0.027} M_{\oplus}$, or approximately 3 lunar masses. Phoebe suggests a population of compact, lunar-mass objects associated with the dark matter distribution of the Milky Way, and potentially opens a new window to the physics of inflation.

Gaia21bja is a Gaia alerted young stellar object (YSO) that exhibits at least seven quasi-peridoic brightenings over a 20 year-long light curve with durations of 1.5-2 years and amplitudes up to $\sim$1.7 mag in the Gaia $G$-band. We analyze its optical and near-infrared photometry and spectra taken using the IRTF and VLT in its faint and bright states in order to characterize its physical properties. A Lomb-Scargle periodogram analysis results in a most significant period of $916\pm77$ days. We derived the stellar parameters as $R_\star= 0.78 \pm 0.04~R_\odot$, $L_\star=(4.5\pm0.3) \times 10^{-2}~L_\odot$, and $M_\star= 0.16 \pm 0.03~M_\odot$. The spectra taken during the burst are dominated by emission lines and are similar to those of EX Lupi-type eruptive young stars (EXors). We found that the accretion luminosity and mass accretion rate increased by a factor of $5.5-6$ during the burst. Based on this, and the quasi-periodic bursts, we suggest that Gaia21bja is an eruptive YSO, and is most consistent with the `Periodic' category of the Outbursting YSOs Catalogue.

We investigate the interplay between morphology, specific star formation rate (sSFR), and local environment using a sample of 7,408 galaxies from the SDSS-IV MaNGA survey. Our analysis spans stellar masses from dwarf to massive galaxies, enabling a unified view of how stellar mass and environment regulate galaxy evolution. Galaxies are classified by morphology (ellipticals (E), lenticulars (S0s), early-type spirals (ETS), and late-type spirals (LTS)) and local environmental density, with star formation activity traced using sSFR. Low-mass galaxies ($\log (M_{\star}/M_{\odot}) < 10$) are predominantly star-forming and dominated by LTS, whereas high-mass galaxies ($\log (M_{\star}/M_{\odot}) \geq 10$) are dominated by ETS and are largely quenched. By separating dwarf ($\log (M_{\star}/M_{\odot}) \leq 9.5$) and intermediate-mass galaxies ($9.5 < \log (M_{\star}/M_{\odot}) < 10$), we find that dwarf galaxies remain predominantly star-forming with only weak environmental dependence, whereas intermediate-mass galaxies exhibit clearer environmental trends toward quenching. Using the D4000 index as a tracer of long-term stellar population aging, we further show that dwarf E and S0s host systematically younger stellar populations than their intermediate-mass counterparts, implying reduced quenching efficiency and more gradual environmental processing in the dwarf regime. This distinction is not evident among spiral galaxies, whose stellar population properties are comparatively insensitive to the dwarf versus non-dwarf classification. Overall, these results indicate that the commonly defined low-mass galaxy population is not homogeneous and that dwarf and intermediate-mass galaxies show systematically different evolutionary trends. Treating them separately is therefore essential for interpreting galaxy evolution in the low-mass regime.

Nuclear star clusters (NSCs) are dense stellar environments located in the center of most galaxies. NSCs are thought to form through two primary methods; through the inspiral of globular clusters (GCs) to the galactic center due to dynamical friction, and through in-situ star formation. Recent observations of dwarf elliptical galaxy UGC 7346 suggest that it might be undergoing NSC formation due to the presence of multiple GCs near its photometric center. We perform direct N-body simulations of nine GCs belonging to UGC 7346's GC system to investigate whether their eventual infall to the galactic center would result in the formation of a NSC. Our simulations indicate that GC inspiral leads to the formation of a central stellar over-density relative to the background profile of the host galaxy within $\sim$1.5 Gyr, corresponding to a NSC with a typical mass of $(4.1\text{-}4.5)\times 10^5 M_{\odot}$. Several key structural parameters of the newly formed NSC, including the Sérsic index, effective radius, and central stellar density, lie well within the range observed for NSCs. We also test a hypothetical scenario in which some of the infalling GCs have larger masses ($M \sim 10^6 M_{\odot}$), resulting in the formation of a more massive NSC whose mass and size are more consistent with observations. Our results suggest that inspiral of GCs is a viable channel for assembling a significant mass in the shape of NSC in the center of dwarf galaxies and that UGC 7346 will host a NSC at its center in about 2-3 Gyr.

Spectro Capture is a Python-based software system developed to automate small- observatory fibre-fed spectroscopy. The system integrates target selection, telescope slewing, guide star acquisition, fibre position restoration, calibration and science exposure sequencing within a single observing workflow. The paper describes the design and operational behaviour of the system at Shenton Park Observatory, where it has been used for unattended multi-target spectroscopic observing. Log analysis from January to May 2026 shows that the system completed 339 of 345 attempted science target blocks in primary unattended batch runs, corresponding to a completion rate of 98.3%. The results demonstrate that reliable unattended spectroscopy is practical at a small observatory when target acquisition, guiding, scheduling and calibration control are treated as an integrated software problem.

Primordial black holes (PBHs) have attracted intensive research interest as a promising candidate of dark matter. However, because of the Hawking radiation, the PBHs lighter than $10^{15}~\rm{g}$ have already evaporated before today. To extend the PBH mass window to small-mass range, two possible ingredients are explored. The first is the consideration of regular PBHs with non-singular metrics, which can decrease the Hawking temperature, thereby lowering black hole evaporation. The second is the incorporation of the memory burden (MB) effect, which can further suppress the evaporation rate, after regular PBHs have lost a certain amount of their initial masses. In this work, we combine these two ingredients and study the MB effects of three types of regular PBHs (the Hayward, Bardeen and Simpson--Visser black holes). Assuming a phenomenological self-similar evaporation, we find that the MB effect significantly relaxes the evaporation constraints. For a benchmark of the MB strength parameter $k=1$, a new PBH mass window opens at around $10^6$--$10^8$ g, where regular PBHs can compose all dark matter without violating the Big Bang nucleosynthesis bounds.

Sunspots are transient, magnetically intense features that host oscillations linked to magnetohydrodynamic (MHD) waves. These waves may contribute to plasma heating and drive mass flows in the solar wind. Beyond their energetic role, they serve as diagnostic tools for probing sunspot structure. In this study, we investigated chromospheric wave propagation in a sunspot using high-resolution, multi-wavelength observations from the Goode Solar Telescope at Big Bear Solar Observatory. Spectral analysis shows that the intensity at H$\alpha$ line core and its wings exhibited oscillatory signal at about 3 min. We performed a cross-wavelet analysis to examine the phase relationship between the wing-integrated and line-core intensity oscillations of the H$\alpha$ line and the centroid-derived H$\alpha$ Doppler velocity. We also analyze the phase relationships between intensity pairs from different passband combinations of the H$\alpha$ line. The results indicate the presence of slow magnetoacoustic modes manifesting standing waves along with upward propagating waves. The observed phase patterns suggest that umbral waves are confined within a non-ideal acoustic resonator, providing measurable wave properties that could serve as input for sunspot seismology and refine models of sunspot atmospheric structure.

We present an analysis of the time delay between the radio emission and the X-ray Compton luminosity during the 2010-2011 outburst of GX 339-4. Using the interpolated cross-correlation function (ICCF), we measure the time delay between the Compton luminosity and the radio luminosity, and find that during the rising hard state, the radio emission precedes the Compton luminosity by approximately 3 days. In contrast, in the decaying hard state, the radio emission lags behind the Compton luminosity by about 8 days. By estimating the mass accretion rate and the disk truncation radius, the calculated inner magnetic field can account for both the radio delay in the decaying hard state and the radio precedence in the rising hard state. The time delays observed in different outbursts across multiple sources are compared further, and the underlying physical mechanisms account for this difference are discussed. These results provide insights into the evolving coupling between the inner accretion flow and the jet in black hole X-ray binaries.

The cosmic dipole observed in the cosmic microwave background (CMB) is traditionally interpreted as being caused by the observer's motion relative to the background. However, tensions with dipole measurements from radio galaxy counts motivate the need for independent probes. This work investigates the feasibility of using strongly lensed gravitational wave (GW) events to measure the cosmic dipole. Strongly lensed GWs produce multiple time-delayed images, which can be used to infer the distances to both the lens and the source. These distances, associated with the observed redshifts of the lens and the source from galaxy catalogues, encode information about the background cosmology and cosmic dipole effects. By reconstructing a statistical sample of doubly lensed GW events based on the singular isothermal sphere lens model, the cosmic dipole can be estimated jointly with background cosmological parameters. Using realistic simulations for Einstein Telescope and Cosmic Explorer, we forecast that a dipole magnitude $g$ consistent with both the CMB and number count measurement could be detected with 10 years of observation. Furthermore, constraints on $g$ are greatly improved by combining constraints from doubly lensed events with those from triply or quadruply lensed events. In the most optimistic scenario, where we measure the number count dipole magnitude with 10 years of observation, we obtain $g = (2.45^{+1.53}_{-1.28}) \times 10^{-3}$ from the combined constraint, provided that systematic uncertainties are mitigated. Although challenging, strongly lensed GWs offer a novel approach to measuring the cosmic dipole, providing an independent consistency test with different systematics from electromagnetic probes.

Radio technosignature searches and radio-based ultra-high energy (UHE) neutrino experiments address different scientific questions, but share a closely related data analysis problem: identifying rare signals of unknown morphology within large datasets dominated by thermal noise and anthropogenic radio-frequency interference (RFI). UHE neutrino radio experiments (including ARA, RNO-G, ANITA, and PUEO) have developed advanced methodologies for continuous-wave (CW) mitigation and background characterization. This invited contribution makes that connection concrete through three points. First, we demonstrate that catalog-based time-domain sine subtraction -- the CW mitigation technique used in ANITA and ARA -- can be adapted for technosignature pipelines by restricting subtraction to documented persistent contaminants, improving broadband transient visibility while preserving uncataloged narrowband candidates. Second, we identify a structural equivalence between spatiotemporal clustering used in UHE neutrino experiments and direction/cadence-based RFI rejection in radio SETI, proposing a joint feature space incorporating direction, time, frequency, bandwidth, duration, and polarization. Third, we argue that background-only anomaly ranking is the natural second stage of this workflow, providing morphology-agnostic candidate triage. Together, these ideas motivate a 'preserve-then-rank' workflow for commensal rare-event discovery, opening a near-term path toward cross-community collaboration.

Precision geodesy relies on the stability of the International Celestial Reference Frame (ICRF), yet its reference sources, Active Galactic Nuclei (AGN), exhibit changes in source structure that can manifest as apparent shifts in their astrometric positions. The high-precision radio measurements used to maintain the ICRF therefore provide a means to investigate the astrophysical mechanisms driving these changes. In particular, the observed astrometric variability offers a unique opportunity to link positional shifts in AGN to high-energy astrophysical processes. We investigated the relationship between the astrometric positions of ICRF AGN and their Gamma-ray emission. We measured the positional offsets of radio cores relative to the ICRF3 at both S/X and K bands and compared them to Fermi-LAT Gamma-ray fluxes within +/-30 days. Out of 92 radio sources, we identified 57 that had enough overlapping data. We find a high incidence of statistically significant (p<0.05) power-law correlations, with ~90% of sources exhibiting this behaviour. The nature of this correlation is complex: we observe both positive and negative correlations, and the sign of the correlation can differ between the two frequency bands for the same source. To explain the correlations, we tested variable Gamma-ray emission locations, changes in nuclear opacity, and variations in jet position angle. Our analysis reveals no single explanation and suggests a complex interplay of multiple physical mechanisms. A search for time lags between the radio position offsets and Gamma-ray fluxes revealed tentative - and highly caveated - evidence for a time-delay in only five sources. A statistical comparison with the OCARS catalogue shows that, although our sample is biased towards optically brighter sources with better-constrained astrometric solutions, it remains representative of the broader AGN population in terms of redshift.

Solar internal rotation at high latitudes is poorly constrained by acoustic-mode helioseismology. Global inertial modes observed on the Sun are highly sensitive to solar differential rotation and may provide new diagnostics of rotation in these regions. We aim to constrain solar rotation with the measured frequency of the $m=1$ high-latitude inertial mode, starting from the HMI/SDO reference rotation profile given by p-mode helioseismology for 2010-2024. Using a validated and accurate eigenvalue solver, we compute the perturbation to the mode frequency resulting from localised changes in the differential rotation rate throughout the solar interior. We find that the linear sensitivity kernel of the $m=1$ high-latitude mode peaks at latitude 75 deg and radius $0.8 R_\odot$, with full widths of 7 deg and $0.13 R_\odot$. From the observed mode frequency in the Carrington frame, $-87.9 \pm 1.9$ nHz (retrograde, averaged over 2010-2024), we infer that the solar rotation rate near this location is $365.3\pm 2.0$ nHz, which exceeds the reference p-mode estimate by $8.1$ nHz. Additionally, we propose a latitudinally smooth, radially independent modification to the rotation rate at high latitudes beyond the linear (small-perturbation) regime. This work demonstrates that individual inertial modes can provide direct constraints on rotation in the bulk of the solar convection zone, well below the surface, representing the first example of spatially resolved inertial-mode helioseismology.

The Telescope Array (TA) is the largest ultra-high energy cosmic ray (UHECR) observatory in the Northern Hemisphere. Together with its extensions, TA Low Energy (TALE), TALE infill, and the TA$\times$4 array, it measures extensive air showers (EAS) initiated by UHECRs across an energy range spanning from $10^{15}$ eV to beyond $10^{20}$ eV. All components of the experiment employ a hybrid detection approach, combining plastic scintillator arrays that sample the EAS footprint at ground level with telescopes that record fluorescence and Cherenkov light from shower development in the atmosphere. The ongoing construction of TA$\times$4 will significantly increase statistics at the highest energies by expanding the surface detector area by a factor of four. In addition, the recently deployed TALE infill array further lowers the hybrid energy threshold of TALE. This presentation summarizes the current status of the TA experiment and highlights recent findings on the energy spectrum, mass composition, and anisotropy.

We introduce a fast ($\sim 1$-$50$ ms) and generalised framework for modelling interaction-powered transients. The framework solves the thin-shell equations of motion for ejecta colliding with circumstellar material (CSM), and supports arbitrary CSM density and velocity profiles, including steady winds, eruptions, and complex time-variable mass-loss histories. For optical/UV lightcurves, we implement two luminosity treatments: a fast one-zone mode based on the thin-shell shock power, and a finite-shell transport mode that evolves trapped radiation, photon diffusion, shock emergence, and post-emergence cooling for finite, static CSM shells. In a benchmark comparison, the transport calculation and an optional time-dependent shock-efficiency prescription reproduce the main qualitative and quantitative features of a one-dimensional radiation-hydrodynamical simulation. We use the same shock solution to post-process radio synchrotron and thermal bremsstrahlung X-ray predictions, enabling self-consistent multi-wavelength diagnostics. We show that the assumed CSM velocity structure can significantly affect inferred parameters even when the density profile at explosion is identical, and that aspherical CSM can mimic multiple spherical shells in bolometric lightcurves. We demonstrate the framework through recovery of a synthetic time-variable mass-loss history and applications to six transients: the Type IIn SN~2010jl, the rapidly evolving stripped-envelope SN~2023xgo, the Type Ia-CSM SN~2020aeuh, the hydrogen-poor superluminous SN~2015bn, the eruptive LBV-like transient SN~2009ip, and the long-duration interacting event iPTF14hls. The inferred CSM structures span steady or enhanced winds, thermonuclear interaction, eruptive density enhancements, and highly structured pre-supernova mass loss, illustrating the framework's utility for inference on upcoming large samples of interacting transients.

The stellar orbital anisotropy parameter ($\beta_{\rm ani}$) is a persistent systematic uncertainty in galaxy-scale strong gravitational lensing (SGL) cosmology. Typically fixed to isotropy or a local prior, it frequently degenerates with the lens density profile. We demonstrate this apparent redundancy largely arises from incomplete photometric constraints. We cross-matched 130 SGL systems with the Pantheon+ SN~Ia compilation, constructing a strictly matched sample of 107 SGL-SN pairs using a 5\% comoving-distance tolerance. Assuming a flat universe ($\Omega_k = 0$), the distance ratio is derived from apparent magnitude differences between paired SNe~Ia, eliminating $H_0$ and absolute magnitude dependence without fitting explicit dark-energy models. To break the kinematic degeneracy, we incorporate lens-by-lens luminosity density slopes ($\delta_i$) from high-resolution imaging. Adopting the quasi-model-independent P2 redshift-evolutionary framework ($\gamma(z) = \gamma_0 + \gamma_z z$), we find very strong statistical evidence for a free $\beta_{\rm ani}$. Fixing $\beta_{\rm ani}$ to isotropy ($\beta=0$) or a local prior ($\beta=0.18$) is strongly disfavored ($\Delta\bic = 14.2$ and $48.9$) and artificially inflates intrinsic scatter. A complementary P3 framework ($\gamma(z,\tilde{\Sigma}) = \gamma_0 + \gamma_z z + \gamma_s \log_{10}\tilde{\Sigma}$) confirms these penalties ($\Delta\bic = 13.5$ and $49.1$). Across all P2 variants, we consistently detect a negative redshift evolution of the density slope ($\gamma_z \approx -0.42$ to $-0.46$; ${\sim}1.5{-}2.0\sigma$), indicating ETG density profiles become shallower at higher redshifts. We conclude that when individual photometric constraints are incorporated, $\beta_{\rm ani}$ is statistically required as a free parameter to prevent severe dynamical modeling biases.

We perform a comprehensive analysis of lunar impact flash (LIF) light curve shapes and their dependence on the lunar terrain, using the large sample of LIFs detected by NELIOTA over the last 9 years. We classified 124 multi-frame light curves into mare, highland and `border' regions. Subsequently, we derived analytical expressions for single-size and dual-size ejecta cooling models, which were fitted to the observational data to estimate their physical properties. While impacts on both terrains yield similar peak magnitude distributions, their decay behaviour differs significantly; highland LIFs exhibit a shallower and longer-lasting decay compared to mare flashes, which are faster and steeper. The dual-size model suggests this extended duration is primarily driven by the fine droplets of the ejecta. The profile and duration of the LIF light curves represent the initial stages of the impact cratering process. The observed dichotomy between highland and mare LIFs demonstrates that the initial stages of the impact cratering process are fundamentally dependent on lunar lithology.

Intermediate-mass black holes (IMBHs) are the missing link in the cosmic hierarchy of black holes, bridging the gap between stellar-mass black holes and supermassive ones. They also serve as unique laboratories for testing strong-field gravity and are prime targets for future multi-messenger observations. However, IMBHs are a population that has remained notoriously difficult to detect. The microlensing effect of fast radio bursts (FRBs) can serve as a clean and powerful method to probe IMBHs. In this work, we develop a pipeline to search for microlensed FRBs based on their dynamic spectra and apply it to the CHIME/FRB Catalog 2. Two microlensing signatures have been identified in two separate sources, i.e. FRB~20190131D and FRB~20211115A. The inferred lens masses for these two signatures are $\sim[539-609]~M_{\odot}$ and $\sim[1544-2571]~M_{\odot}$, respectively. Here we interpret them as evidence for IMBHs. If there are no intervening structures-such as galaxies or clusters-along the line of sights for these two sources, the two identified IMBHs might be isolated and of primordial origins. In that case, we obtain primordial black holes (PBHs) within these two mass ranges would constitute $\sim4\%$ of dark matter. Moreover, if these two candidates are not genuine lensing signatures, the abundance of intermediate-mass PBHs with masses $>300,M_{\odot}$ is constrained to be $\sim13\%$ at $95\%$ confidence level. Therefore, more comprehensive observational information for FRBs, together with a deeper understanding of whether the intrinsic emission mechanisms of FRBs can produce lensing-like signals, will be crucial for establishing this effect as a powerful tool for probing (primordial) IMBHs.

Disentangling SF and AGN emission is essential for understanding galaxy evolution, yet remains challenging in merging systems where both processes are enhanced and spatially intertwined. Galaxy mergers drive gas inflows that simultaneously fuel nuclear SBs and BH accretion, shaping morphology from nuclear ($\lesssim 250$~pc) to large-scale ($\gtrsim 500$~pc) regions. Radio interferometry provides an unobscured view, but separating compact nuclear SBs, AGN, and diffuse SF requires multiscale, multi-frequency observations. We present a systematic method to characterise multiscale radio properties in 15 local ($z\lesssim 0.1$) U/LIRGs ($L_{\mathrm{IR}} > 10^{11}\mathrm{L}_{\odot}$). Using \emph{e}-MERLIN and VLA at 1.4, 6.0 and 33.0~GHz, we probe physical scales from $\sim 10$--$250$~pc to $\sim 0.5$--$3.0$~kpc. We decompose radio emission into nuclear (compact cores and nuclear extended) and large-scale (total and diffuse) components, comparing morphological properties (emission fractions, sizes, luminosities, surface densities) and investigating correlations with source classes, merger stages, and $L_{\mathrm{IR}}$. We find: i) nuclear emission contributes $\sim$50\% of total radio emission on average; ii) total multiscale diffuse emission (SF-related) contributes $\sim$80\% to total power; iii) nuclear emission components act together to correlate with total radio and infrared luminosities, which increase with merger stage, whilst diffuse emission at larger scales shows no clear dependence on nuclear processes; iv) sources with radio excess (lower $q_{\mathrm{IR}}$) show lower nuclear luminosity ratios $L_{\mathrm{R,33}}^{\mathrm{N}}/L_{\mathrm{R,6}}^{\mathrm{N}}$, indicating a deficit of high-$\nu$ radio emission; since 33.0~GHz traces recent SF, this suggests the radio excess is dominated by non-thermal emission at lower $\nu$, likely AGN-related, rather than enhanced SF.

We have carried out a stacking analysis with the COSMOS-Web catalogue on one of the deepest ever SCUBA-2 images at 850 microns, allowing us to estimate the mean submillimetre flux density for samples of galaxies split by stellar mass and morphological class over the redshift range 0<z<12. For all morphological classes, the mean star-formation rate estimated from the dust emission increases with redshift, reaching a value for the most massive galaxies (~10^11 soar masses) of >~80 solar masses per year at 2 < z < 4.5. In this redshift range, the mean star-formation rate for these galaxies falls along the Hubble sequence from ~280 solar masses per year for irregular galaxies at one end to ~80 solar masses per year for spheroids at the other end, which shows that quenching was already happening shortly after the emergence of the Hubble sequence. The decrease in the star-formation rate for the spheroidal galaxies can be reproduced with a `starvation' quenching model with a depletion time of ~10^{8.2} years. We also show that the transformation of `submillimetre galaxies' can reproduce the growth in number-density of massive bulge-dominated and spheroidal galaxies over the redshift range 1.5 <z < 4. As a side-project, we have used our stacking results to show that the ratio of dust mass to stellar mass in galaxies increases with redshift out to z~8 and to determine the relationship between the mean density of dust and redshift in the range 0 < z <12. We show that a chemical evolution model based on the `star-formation history' of the universe, with a gas outflow rate equal to the star-formation rate, can explain the monotonic rise in the dust-to-stellar mass ratio and reproduce the relationship between mean dust density and redshift remarkably accurately.

Elemental abundances provide a powerful diagnostic of the physical mechanisms and processes that heat the solar atmosphere and drive the solar wind. The First Ionisation Potential (FIP) effect and its inverse (IFIP) are observed both on the Sun and other stars however, the underlying fractionation mechanisms, their dependence on the magnetic field topology, and the role of wave dynamics and turbulence in the chromosphere are not entirely understood. To address these challenges, a focused team, including observers, theorists, modellers and instrument scientists, spanning a range of career stages and institutions, came together for the Royal Society Theo Murphy meeting ``Solar Atmospheric Abundances in Space and Time". As a result of this meeting, the team worked in collaboration to produce 16 publications for this Special Issue. These publications are introduced here, including a discussion of the open questions and future directions in the context of advances in numerical modelling and current and upcoming solar and stellar missions. This article is part of the Royal Society Theo Murphy Meeting Special Issue ``Solar Atmospheric Abundances in Space and Time''.

We present the Small Phased Array DEmonstrator (SPADE), a compact phased array spectrograph designed for the monitoring of solar activity in the decameter range, where a majority of bursts of interest for solar physics and space weather occur. Unlike other existing phased arrays, SPADE operations like Sun tracking and spectra productions are made entirely digitally, and rely on the use of open source library software and commercial generic software defined radio receivers. In this paper, we describe the instrument and its associated software and present the first observations of the Sun and of the jovian magnetosphere performed between September and December 2024. For the solar observations, in particular, we highlight the capacities of the instrument in terms of temporal and frequency resolution, making it capable of addressing science topics like turbulence in the coronal medium and propagation effects affecting radio waves, which are currently only accessible by larger but non-solar dedicated phased arrays.

A sub-class among Young Stellar Objects (YSOs), known as FU Ori type stars, undergo sudden rises in luminosity by several orders of magnitude on timescales of a few months to a few years, and decay back to quiescence on timescales of a few decades. Modelling the light curves of these objects is crucial to understanding how different components of these accretion disk systems evolve during outburst. For this purpose, we use a parametric model that couples the stellar photospheric emission, magnetospheric accretion shocks, an irradiated dust disk, and a viscously heated gas disk. We adopt time-dependent accretion rate profiles that mimic the observed morphologies of FU Ori outburst light curves, and we use the accretion model infrastructure to simulate multi-band light curves, as well as color curves. The model enables us to study how different components dominate the flux in each band over the course of an outburst, providing insight into star-magnetosphere-disk interactions throughout the outburst cycle. We find that throughout an accretion outburst, red optical and near-infrared lightcurves generally follow the same or very similar form as the input accretion profile, being sensitive to heating in the accretion shocks and inner gas disk, while mid-infrared lightcurves are more responsive to the location and heating of the innermost dust disk.

(1) Background: for the study of highly magnetized neutron stars observed as magnetars, and to quantify the effect of this intense magnetic field onto the star's structure and shape which can be particularly relevant for the study of emission of continuous gravitational waves, both numerical and perturbative approaches have been developed. (2) Methods: we compare these two approaches in General Relativity with the limitation to the case where the magnetic field has a purely poloidal structure. The perturbative one (Konno-99) assumes that the deformation induced by the magnetic field is small and that this field arises only from dipole currents. The full numerical one is based on the library LORENE. (3) Results: we have used both approaches to compute the magnetic field distribution and the deformation of the star, varying the value of the magnetic field at the pole, the compactness of the star and its equation of state. (4) Conclusions: whereas the perturbative approach breaks down for very high polar magnetic field values (typically above a few times $10^{16}$ G), it gives very good results for observed values, even in magnetars. On the contrary, the numerical code exhibits resolution problems for relatively low magnetic field values (typically $10^{10}$ G), which translates into imprecise computation of the star's deformation and mass quadrupole moment.

The lensing imprint of cosmic voids in the Cosmic Microwave Background (CMB) provides a powerful test of the $\Lambda$CDM model. However, recent studies report a "lensing-is-low" tension between observations and mock predictions. To investigate this, we measure the stacked CMB lensing signal of 3D cosmic voids identified in the DESI Legacy Survey DR9 Luminous Red Galaxy (LRG) sample, cross-correlated with the Planck 2018 lensing map. We compare our observations to $\Lambda$CDM templates derived from Buzzard mocks, critically calibrated using over one million DESI spectra to perfectly match the sparseness and photometric redshift error distributions of the observed data. By categorizing voids based on their gravitational potential, we disentangle the signals of void-in-voids and void-in-clouds, achieving two independent record detection significances of $\sim 17\sigma$. We find full agreement between observations and simulated $\Lambda$CDM templates across all void populations and redshift bins ($0.35\! <\! z\! <\! 0.95$), measuring an amplitude parameter $A_\kappa = 1.016 \pm 0.054$ for the full sample. This highlights the necessity of accurate systematic control, effectively resolving the reported lensing tension within this dataset. This proceeding summarizes the results presented in Sartori et al. (2025).

Primordial Black Holes (PBHs) forming prior to Big Bang Nucleosynthesis evolve during the reheating epoch, an environment dominated by an oscillating inflaton field decaying into a relativistic thermal bath. In this work, we track the complete lifecycle of PBHs within this coupled inflaton-radiation background. Utilizing $\alpha$-attractor E-models, we analytically anchor the reheating initial conditions directly to Cosmic Microwave Background observations. By matching exact scalar field solutions in a Schwarzschild spacetime to the cosmological far-zone, we derive the cycle-averaged mass accretion rate and couple it to the growing radiation bath. We find that this combined accretion induces a highly non-linear enhancement of the final PBH mass. Because the Hawking evaporation timescale scales cubically with mass, PBHs forming near their critical runaway limits experience a massive extension of their lifespans. Surviving deeper into the radiation-dominated era triggers a multi-order-of-magnitude amplification in their emitted Stochastic Gravitational Wave Background (SGWB).

We investigate a dark energy model driven by a planar domain-wall-like structure with a thickness comparable to, or larger than, the current Hubble radius, focusing on its intrinsic anisotropy and observational viability. Near the centre of the domain wall (DW), the spacetime is anisotropic, with distinct expansion rates parallel and perpendicular to the wall. This anisotropic structure induces direction-dependent cosmic expansion and modifies photon geodesics from cosmological sources, leaving characteristic signatures in cosmological observables. We confront the model with recent observational data. We first compute the anisotropic Cosmic Microwave Background (CMB) temperature multipoles generated by the DW and impose constraints from the Planck 2018 measurements. These constraints severely limit the allowed DW abundance, requiring the DW energy density to be less than $\mathcal{O}(10^{-5})$ of the current critical density in order to suppress the quadrupole contributions. We then perform a Markov Chain Monte Carlo (MCMC) analysis using Type Ia supernova (SNe Ia) data, including the Pantheon+ SH0ES and DESY5 samples, to compare the DW scenario with the standard $\Lambda$CDM model. We find that although the DW naturally realises anisotropic accelerated expansion, the combined constraints from the CMB and SNe Ia favour the $\Lambda$CDM limit, in which the DW contribution is negligible, and the universe is effectively isotropic. Our results demonstrate that a Hubble-scale domain wall is tightly constrained by current observations and can only play a subdominant role in the late-time cosmic acceleration.

We characterise position-error behaviour of Two-Line Element (TLE) propagation against operator-updated truth on Starlink, sweeping 24,641 next-TLE-truth pairs across 501 satellites stratified by altitude shell (540, 550, 560 km) and platform generation (v1.0, v1.5, v2-mini) over April 2026. Each pair is propagated with SGP4 and GMAT at high fidelity (EGM2008 $70\times70$, NRLMSISE-00 drag, Sun and Moon third-body gravity, conical-shadow SRP), then compared against the operator's next TLE as proxy truth. Three findings: First, position error follows a per-cell power law $\lVert\Delta\mathbf{r}(\Delta t)\rVert \approx A\,\Delta t^{k}$ with fitted exponents in $(1,2)$ on every v2-mini cell and on the high-fidelity v1.x cells at 540 and 560 km, while SGP4 v1.x and high-fidelity v1.x at 550 km are sub-linear ($k \lesssim 1$); the cohort-specific mix of mean-motion bias and unmodelled in-track acceleration sets the per-cell exponent. Pooled $L_{2}$ medians grow from $\sim 1$ km at 6 h to $\sim 38$ km (SGP4) / $\sim 76$ km (high-fid) at 7 d. Second, high-fidelity propagation from public-TLE inputs does not improve over SGP4 at any of the four staleness horizons; SGP4 wins on $\sim 65$--$75\%$ of pairs, with v2-mini at long $\Delta t$ the one regime where high-fidelity wins on a majority of pairs at both populated shells. The negative result reflects operator-OD residual dominance at epoch, SGP4-vs-SGP4 truth-construction kernel alignment, and spacecraft-property bias amplification on the high-fidelity arm. Third, the per-satellite SGP4 staleness coefficient regressed against F10.7 returns a positive slope clearing conventional significance at one shell (560 km) on the 30-day, $\sim 17$ sfu window -- direction-consistent with the LEO density-gradient expectation, not a calibrated F10.7-modulation measurement.

With the increasing precision of asteroseismic observations, it becomes possible to reliably measure oscillation properties of an increasing number of stars. Interpreting these measurements requires a good theoretical understanding of their link to fundamental stellar properties. In this study, we focus on the phase offset in gravity(g)-mode frequencies, which is imprinted in the asymptotic eigenfrequency pattern of mixed dipole modes observed in red--giant branch stars. We aim to unravel its physical origin and thus enable an informed interpretation of observations. Using stellar models, we empirically test the contribution of the g-mode offset $\varepsilon_\mathrm{g}$ (which is related to the wave reflection at cavity boundaries and commonly considered to be the dominant phase term) and glitches to the total observable phase. We find that, additionally to $\varepsilon_\mathrm{g}$, buoyancy glitches play an important role in the correct interpretation of the g--mode frequency phase. We further find that glitches in the evanescent zone also contribute to the phase, and we present a formalism to quantify this contribution. Finally, we propose a modification to the widely used formula for $\varepsilon_\mathrm{g}$. The g--mode frequency phase carries more information than previously considered. It has large analytic potential to study not only the reflection properties of the buoyancy cavity, but also the properties of glitches in the Brunt-Väisälä frequency.

Despite the low solar irradiation it receives, Neptune shows a very active atmosphere with some of the most intense dynamics observed in Solar System atmospheres. Characterizing the atmospheric temperature profiles of the planet is a key to understand these observed processes. In this work, we derived the Neptune pressure-latitude thermal field, using 2016 ALMA measurements of the CO(3-2) spectral line at 345.796 GHz, with a spatial resolution of about 0.37" on Neptune's 2.24" disk. To analyse the data, we developed MCMC retrieval methods to derive both the temperature profiles and the CO abundance profile. We find that our data probes the upper troposphere and the lower stratosphere of the planet, between 2 bar and 0.1 mbar. Although temperature and CO profile are strongly correlated, simultaneous retrievals of both parameters for disk-integrated observations reveal a factor of 2-3 larger CO abundance in the stratosphere than in the troposphere, reinforcing the hypothesis of CO delivery by a recent cometary impact. CO retrievals with fixed temperature profile do not fit the spatially-resolved observations, implying underlying temperature variations. By performing temperature retrievals with spatially constant CO, we find distinct trends of the thermal profile between the southern polar region, mid-latitudes and the equator. At 10--100 mbar, this structure is consistent with that observed by Voyager 2 in 1989, i.e. colder mid-latitudes, with a warmer equator and south polar regions. At 300-600 mbar, however, we find a cold layer of about 45 K at Southern polar regions (-80°) which disappears towards mid-latitudes and the equator. By probing pressure levels not easily accessible to other observing methods e.g. infrared sounding, submillimeter observations offer a new view of the complex thermal structure in the upper troposphere and lower stratosphere of the Icy Giants.

Accretion-disc corona sources are high-inclination low-mass X-ray binaries in which the innermost regions are hidden and the observed X-ray emission is dominated by radiation scattered above the disc. 4U 1822-37 is a key binary system of this class, but its geometry is still debated. X-ray polarimetry offers a direct probe of the scattering structure. We present the first X-ray spectro-polarimetric study of 4U 1822-37 and test whether its spectrum and polarisation can be explained within the accretion-disc-corona scenario. We analysed a coordinated campaign with IXPE, XMM-Newton, NuSTAR, and Swift. We performed broadband spectral modelling and model-independent, energy-resolved, and orbital-phase-resolved polarimetric analyses. The broadband spectrum requires a soft thermal component, a Comptonised continuum, a hard power-law tail, and relativistically blurred reflection. The observed 0.1-100 keV luminosity, $L_{\rm obs}\simeq6.1\times10^{36}\ {\rm erg\,s^{-1}}$, is far below the intrinsic luminosity expected from the orbital evolution, supporting a geometry in which only a small fraction of the intrinsic emission is scattered into the line of sight by an extended, optically thin corona. In the 2-8 keV band, IXPE measures ${\rm PD}=7.9\pm0.6\%$ and ${\rm PA}=-24^\circ\pm2^\circ$. The PD increases with energy, while the PA remains approximately constant. During eclipse, the PD decreases to ${\rm PD}=5.5\pm1.7\%$, with no significant PA variation. This behaviour is consistent with the companion occulting the most polarimetrically efficient part of the extended corona. The high PD, stable PA, energy-dependent polarisation, and eclipse behaviour support a picture in which 4U 1822-37 is observed in an extreme high-inclination, scattering-dominated regime. The extended corona is the main structure shaping both the observed X-ray emission and its polarisation.

Many barred galaxies, including the Large Magellanic Cloud (LMC), display strong lopsidedness and off-centre bars. The dynamical connection between bar-disc misalignments, internal mass asymmetries, and arm morphology is not yet fully characterised. We investigate how internal mass imbalances within the bar and global offsets between the bar and the centre of mass of the system modify the equilibrium-point structure and the invariant manifolds that organise arms. We construct a barred galaxy model which includes an off-centred and asymmetric in shape bar. Using numerical continuation, we track the position, stability, and bifurcations of the Lagrangian equilibrium points as functions of the displacement of the asymmetric mass component along the bar and of the offset between the bar and the system's centre of mass. For representative configurations we compute the invariant manifolds of planar Lyapunov orbits around unstable points and analyse the resulting arm structures. Internal bar lopsidedness and modest bar-disc offsets that keep the centre of mass inside the bar preserve the classical configuration with five equilibrium points, but strongly distort the associated invariant manifolds, producing two arms with different densities and shapes. The bar-disc offset reaches a threshold at the point at which the galactic centre of mass exits the bar ellipsoid, in which a pitchfork bifurcation removes the collinear unstable points and the system transitions to a three-equilibrium-point configuration in which a single unstable point and its associated manifold supports one arm. This framework is compatible with the observed correlation between off-centre bars and photometric lopsidedness, and it provides a dynamical explanation for the strongly asymmetric arm morphology of galaxies such as the LMC.

We present new multielement NanoSIMS isotopic measurements (C, N, Si, Mg-Al, Ti, and Ni) for four putative nova SiC grains and 79 AB SiC grains from the Murchison meteorite to reassess their stellar origins. High-resolution imaging and a revised Mg/Al relative sensitivity factor for SiC yield substantially improved 26Al/27Al ratios and the most reliable multielement characterization to date for 13C-rich presolar SiC grains. To interpret these data, we computed an expanded suite of hydrodynamic CO, ONe, and recurrent nova models spanning a range of white-dwarf masses and pre-enrichment parameters. When all isotopic systems are considered together (C, N, Mg-Al, Si, Ti, and Ni), the CO nova models provide the closest and most self-consistent match to both the putative nova grains and the subset of AB grains lacking s-process signatures. CO novae of low- to intermediate-mass naturally reproduce the observed 14N/15N-26Al/27Al trend, the Si isotope compositions of AB grains which dominantly reflect Galactic chemical evolution (GCE), and the mild Si isotope shifts in putative nova grains relative to the GCE trend defined by AB grains. In contrast, ONe and recurrent nova models fail multiple isotopic constraints simultaneously. These results demonstrate that low- to intermediate-mass CO novae (0.6-1.0 Msun) are the most plausible stellar sources of 13C-rich SiC dust lacking s-process signatures (1-2\% of all presolar SiC), and they establish a multielement, model-anchored framework for quantifying nova contributions to the dust reservoir in the interstellar medium.

TeV string-mass-scale strings have been excluded experimentally at colliders, as their effects have not been observed at the Large Hadron Collider (CERN). On the other hand, higher-scale string theory, with mass scales typically close to the Planck scale, is often regarded as experimentally inaccessible due to the enormous energies required for direct tests, and far beyond the reach of present or foreseeable particle accelerators. Nevertheless, the early Universe may provide an indirect observational window for high-string scale through imprints left on the Cosmic Microwave Background (CMB). In this work, building on previous studies, we reexamine temperature and polarization angular correlations as probes of the geometry and topology of the pre-inflationary Universe. We focus in particular on two-point correlation functions at large angular scales, where signatures of nontrivial spatial topology may survive as relics of primordial physics. We investigate the observational consequences of toroidal compactification and analyze their impact on the primordial power spectrum of the CMB provided by the Planck satellite. Within the current experimental and theoretical uncertainties, we identify a possible indication closely related to spatial-parity breaking, consistent with the presence of six spatial extra dimensions in the early Universe, compactified at the GUT epoch before the start of inflation. Finally, we extend our framework to B-mode polarization, highlighting its potential as a sensitive probe in forthcoming ground-based and space-borne experiments with unprecedented precision.

Type Ia supernovae are bright thermonuclear explosions of one or more white dwarf stars. The exact origin and explosion mechanism for these supernovae is still poorly understood. In the near-Chandrasekhar mass progenitor model, a simmering phase precedes the explosion. During this simmering phase, central carbon burning heats the core and drives convection. A poorly understood aspect of this phase is the convective Urca process, a linking of weak nuclear reactions and convective mixing. Convective Urca has the potential to alter characteristics of the convection zone and thus alter the evolution of the white dwarf. To study the convective Urca process, we use the low Mach number hydrodynamic code MAESTROeX to run 3D simulations of the convection zone. We build off previous work to implement a more comprehensive carbon burning network and include the A=21, A=23, and A=25 Urca pairs in the simulations. We compare simulations with and without the convective Urca process to isolate the direct effects the process has on the convection zone. We find the convective Urca process reduces the efficiency of convective mixing near the the convective boundary, but does not restrict the size of the convection zone. We additionally find the A=23 Urca pair to be the most important Urca pair to the convective Urca process in these simulations. All together, our results better inform our understanding of this complex phenomena as well as demonstrates the range of potential convective structures, particularly at the convective boundary, of a simmering white dwarf.

Machine-learning methods are increasingly applied to astronomical surveys, providing powerful tools for detecting and studying galaxy clusters. We investigate the mass-redshift properties and completeness of the ComPACT galaxy cluster catalogue, constructed using a convolutional neural network applied to publicly available combined ACT+Planck maps. The ComPACT catalogue contains 2,962 SZ-selected galaxy cluster candidates. We confirm clusters by estimating redshifts using literature information and photometric techniques based on DESI Legacy Imaging Surveys data. Cluster masses are derived from ACT+Planck and Planck Compton-y maps via SZ scaling relations. The completeness is assessed using simulated cluster injections into real microwave maps. We confirm approximately $\sim$60 % of the ComPACT candidates as galaxy clusters. The redshifts span the range $0.007 < z < 1.7$, including approximately 116 new measurements. Masses are obtained for 56 % of the sample, covering the range $(0.25 - 13.1) \times 10^{14} ~M_\odot$ and including 158 new mass determinations. We identify five previously unreported massive clusters ($M_{500c} > 6 \times 10^{14}~M_\odot$) at $z > 0.7$, increasing the known population of such systems by approximately 10 %. The ComPACT catalogue expands the SZ-selected Planck-like cluster population, especially at high redshift and high mass, demonstrating the effectiveness of deep-learning approaches for cluster detection in microwave data.

We present the discovery and characterization of TIC 295741342, a triply-eclipsing triple star system with a giant tertiary. The eclipsing binary consists of two similar main-sequence stars in a 4.75-day orbit. The binary is in a 412.8-day orbit with the giant tertiary. We found two degenerate solutions for the system: one where the tertiary is ascending the Red Giant Branch (RGB), and the other where the tertiary is on the Horizontal Branch (HB) and will eventually ascend the Asymptotic Giant Branch (AGB). In both solutions, the system is near-perfectly coplanar. In TESS Sector 33, the binary passes behind the giant tertiary, producing a distinctive "head-and-shoulders" eclipse that directly constrains the relative flux contributions and radii of all three stars. We modeled the system using a comprehensive spectro-photodynamical model that simultaneously fits the TESS lightcurve, eclipse times, spectral energy distribution, and radial velocities from 48 TRES spectra obtained over four years of observation resolving all three components. Evolutionary analysis using MIST tracks indicates that, in both solutions, the tertiary will overflow its Roche lobe, one in the RGB and the other in the AGB. The Roche lobe overflow will initiate either a stable mass transfer to the binary or a common envelope evolution that will likely result in ejections and/or mergers. Our models predict the midpoint of the next outer eclipse will occur on September 1, 2026 and we encourage follow-up observations with a $\pm$3 day window to observe the full event and further constrain the system parameters.

Spectral-siren cosmology constrains the Hubble constant $H_0$ using gravitational-wave observations of compact-binary coalescences. The method combines luminosity distances inferred from the waveform with redshift information statistically encoded in population features of the source-frame mass spectrum. Because the detector measures redshifted masses, structure in the intrinsic mass distribution acts as an internal ``ruler'', making the inference sensitive to assumptions about the population model. In particular, redshift evolution of the mass spectrum is widely discussed as a potential systematic for $H_0$ measurements. We revisit spectral-siren constraints with the GWTC-4.0 binary black hole catalog, explicitly allowing the main mass scales of a standard parametric mass model to evolve with redshift. We find no compelling evidence for evolution at current sensitivity. Allowing evolution produces a modest, non--statistically--significant shift of the $H_0$ posterior toward lower values, which we interpret with targeted posterior and event-level diagnostics. Importantly, the associated systematic uncertainty is subdominant to that induced by alternative redshift-independent descriptions of the mass spectrum, such as the number of spectral features and the functional form used to model them. Our results indicate that, at current sensitivity, spectral-siren constraints on $H_0$ are robust to redshift evolution of the mass spectrum within the flexibility explored here. Using injection studies, we show that this mild $H_0$ shift is reproduced when a non-evolving underlying population is analyzed with an evolving model, consistent with an over-flexible population description at the present signal-to-noise. The sign and magnitude of the shift can, however, depend on detector sensitivity and redshift reach as the population features become increasingly constrained directly by the data.

We use the latest ANTARES Galactic Ridge neutrino measurements to investigate their implications for indirect dark matter (DM) searches. We consider both annihilating and decaying DM scenarios, spanning a wide range of masses and final states, and systematically compare the resulting neutrino fluxes with the expected astrophysical Galactic diffuse emission. Furthermore, we compare the results for different DM density profiles allowed by the observations, from spike and cuspy to cored profiles. We do so for the WIMP model-independent scenario and explore two more specific models: branons and very heavy sterile neutrinos, where a cold DM candidate arises naturally from the theory. We show the potential neutrino measurements in the Galactic Ridge for DM and make predictions for future neutrino observatories.

Recent observations from the James Webb Space Telescope (JWST) have revealed proximate damped Lyman-$\alpha$ systems (PDLAs) in the foreground of high redshift galaxies ($z \gt 5$), which have been interpreted as neutral circumgalactic media (CGM). The ionization state of the CGM, potentially inferred from DLA strength, may serve as a probe to trace the progress of reionization, similarly to the ionization state of the intergalactic medium (IGM). To determine if this method has merit, we use the cosmological hydrodynamical simulation TECHNICOLOR DAWN to study simulated gas halos at redshifts $z = 10, 8, 6,$ and $5.5$. We investigate the reionization topology to determine whether the CGM and IGM have similar ionization histories, and we study the relation between column density of neutral hydrogen (observationally measured by DLA strength), neutral fraction, and gas mass fraction of the foreground gas to determine whether PDLAs can be used to trace the progress of reionization. We find an inside-out-middle reionization topology, where the CGM reionizes after the IGM and remains partially neutral at $ z= 5.5$. The foreground column density of neutral hydrogen depends mostly on halo mass, with a weak dependence on neutral fraction or redshift. Therefore, provided precise estimates of halo or stellar mass, PDLAs may be used to trace the progress of reionization particularly at high redshifts.

Brown dwarfs are compact objects that do not reach temperatures high enough to produce sustained hydrogen fusion. Consequently, they cool over time, gradually evolving through later spectral types. In fact, three new spectral types (L, T, and Y) were added to the Harvard sequence to accommodate the spectral features of brown dwarfs. During the cooling process, some brown dwarfs unexpectedly become bluer instead of redder (at optical and near-infrared wavelengths). This phenomenon, known as the bluing effect, is particularly noticeable at the L/T spectral transition. The aim of this work is to approximate the spectral type of brown dwarfs using only photometric data, in particular 2MASS and WISE magnitudes. We used two machine learning algorithms, Random Forest and Gaussian Processes, which were evaluated using a 70/30 train/test split. Both models were trained using 5-fold cross-validation and achieved F1-scores of 0.84 and 0.87, respectively, on the test set. After validating the reliability of the algorithms, we applied them to 21 isolated brown dwarfs without prior spectral type determinations. Our results indicate that 5 of these objects have a spectral type between L0 and L4, while the remaining 16 fall within the M6-M9 range. Machine learning algorithms, combined with multi-band photometry, are a powerful tool for estimating the spectral types of brown dwarfs.

We propose an exceptional-point (EP) framework for black-hole ringdown beyond the standard quasinormal-mode (QNM) paradigm. It provides a first-principles characterization of the resonance associated with avoided crossings near EPs, an effect that conventional QNM analysis cannot fully capture. Employing a phenomenological environmental black-hole model with the hyperboloidal framework, we identify near-coalescence of both QNM eigenvalues and eigenfunctions, and directly demonstrate that the resonance produces enhanced mode contributions in the time domain, resulting in characteristic departures from exponentially damped oscillations. Our formulation further reveals that the EP frequency, given by the average of the resonant modes, emerges as the physically relevant observable in the near-EP regime, and offers a robust foundation for modeling and extracting resonant ringdown signals.

The structure of the ground $(2^+)$ and excited $(1^+)$ bound states of the $^8$B and $^8$Li nuclei is studied within the framework of the $\alpha+^3$He($^3$H)+$p(n)$ three-body potential cluster model based on the hyperspherical Lagrange-mesh method. The two-body $\alpha-^3$He($^3$H), $\alpha$-N, and $^3$He($^3$H)-N realistic potentials have been applied from the literature. Convergent theoretical estimates for the three-body binding energy and matter radius have been obtained with the maximal hypermomentum $K_{max}=22$ and 28 for the ground and excited $1^+$ states, respectively. The ANC value of the virtual transition of the $^8$B nucleus is estimated self-consistently by matching the overlap integral of the $^8$B three-body and the $^7$Be two-body wave functions with it's asymptotics. The obtained values are $0.211$~fm$^{-1/2}$ and $0.739$~fm$^{-1/2}$ in the spin 1 and spin 2 channels, respectively. For the ANC values of the $^8$Li nucleus the estimates $0.220$~fm$^{-1/2}$ and $0.774$~fm$^{-1/2}$ are extracted. For the zero-energy astrophysical factor of the direct nuclear capture process $^7$Be(p,$\gamma)^8$B an estimate $22.492\pm0.014$ eV b was obtained based on the asymptotic theory developed by D. Baye [Phys. Rev. C {\bf 62}, 065803 (2000)]. The most important contribution comes from the spin 2 channel with $S^{(2)}_{17}(0)=20.838 \pm 0.014$ eV b, while the spin 1 channel yields $S^{(1)}_{17}(0)=1.654 \pm 0.003$ eV b. These results for $S_{17}(0)$ are in a good agreement with the estimate $20.8\pm0.7{\rm(th)}\pm1.4{\rm(exp)}$ eV b of the SF II, but larger than the recommended value $20.5\pm0.70$ eV b of the SF III. At the same time, our estimate is very close to the value 22.4 eV b used in the most successful Solar Model BAR2M [W.~Yang and Z.~Tian, AJ {\bf 970}, 38 (2024)].

The majoron is the Nambu-Goldstone boson associated with the spontaneous breaking of a global $B-L$ symmetry. Remarkably, the minimal majoron framework can simultaneously address three key empirical indications of physics beyond the Standard Model: neutrino masses, the matter-antimatter asymmetry, and dark matter. In this work, we identify the cosmologically viable region in which majoron dark matter and high-scale thermal leptogenesis can be realised simultaneously. We show that successful leptogenesis plays a central role in making this scenario predictive: by constraining the right-handed-neutrino mass scale, it determines the irreducible freeze-in contribution to the majoron abundance and fixes the size of the couplings relevant for visible dark matter decays. Combining the irreducible dark matter production mechanisms with warm dark matter limits and indirect searches for decaying dark matter, we map the resulting majoron cosmological window and show that future X- and gamma-ray telescopes can probe part of the surviving parameter space.

We investigate the dependence of the plasmoid-mediated magnetic reconnection rate on the magnetic Prandtl number using two-dimensional magnetohydrodynamic simulations of two coalescing magnetic islands. For Lundquist numbers below the onset of the plasmoid instability, the reconnection rate follows the expected Sweet-Parker scaling and decreases with increasing magnetic Prandtl number. However, once the current sheet becomes plasmoid unstable, the dependence on the magnetic Prandtl number weakens considerably. In the fully plasmoid-mediated regime, we find reconnection rates that remain nearly independent of the magnetic Prandtl number over the explored parameter range. We show that the largest reconnection rates are associated with strongly non-linear phases involving plasmoid interactions and mergers. We further compare our results with simulations of the boundary-driven Taylor problem, where previous studies reported a stronger magnetic Prandtl number dependence, and provide a possible explanation for the differing scalings obtained in the two setups. These results may have implications for reconnection-mediated decay in magnetically dominated turbulence and related astrophysical systems.

Proca and Maxwell fields in $d$-dimensional Schwarzschild black holes with anti-de Sitter (AdS) asymptotics are investigated through their linear perturbations and associated quasinormal modes (QNMs) with Dirichlet boundary conditions at infinity. The Proca field equations reduce to one decoupled and two coupled radial wave-like equations. We demonstrate how the Maxwell equations emerge from the zero-mass limit of the Proca system. Several analytical properties of the corresponding QNM spectrum are examined. To compute the QNM frequencies, we employ two complementary numerical methods particularly suited to asymptotically AdS spacetimes. Using these techniques, we determine the QNMs modes of Proca field perturbations in $4$, $5$, $6$, and $7$-dimensional Schwarzschild-AdS backgrounds. As a new result, we find numerically that scalar-type Maxwell perturbations in large $d\geq 5$ Schwarzschild-AdS black holes exhibit purely imaginary low-frequency modes, analogous to those found in vector-type gravitational perturbations. The presence of such modes is especially relevant within the AdS/CFT correspondence, as they correspond to the linearized hydrodynamic regime in the dual conformal field theory. We also analyze the influence of the Proca mass on the QNM spectrum, also emphasizing how Maxwell modes are recovered in the massless limit. The dependence of the spectrum on the black hole radius is explored. In addition, analytic expressions for the QNM frequencies of vector-type and monopole Proca perturbations, as well as Maxwell modes, are derived for small $d$-dimensional Schwarzschild-AdS black holes by matching asymptotic expansions using an intermediate region. These analytic results show good agreement with the numerical findings, confirming, in particular, the existence of purely imaginary low-frequency scalar-type Maxwell modes in large $d\geq 5$ Schwarzschild-AdS spacetimes.

We study the relation between geodesic completeness, the averaged null energy condition (ANEC), and spatial curvature in Friedmann--Robertson--Walker (FRW) cosmology within classical general relativity. Using the affinely parameterized ANEC along radial null geodesics, we prove that non-static flat or open FRW spacetimes in the regular classes considered here cannot be both null geodesically complete and ANEC-satisfying. Bounded oscillatory or cyclic flat/open models do not circumvent the obstruction: the negative affine-ANEC bulk term accumulates over infinitely many cycles, giving \(I_{\rm ANEC}=-\infty\) for non-static periodic cases. Equivalently, within these classes, non-static ANEC-satisfying flat or open models are null incomplete. The sign obstruction is absent in closed \((k=+1)\) FRW geometry, where the positive curvature term enters the affine ANEC identity with the opposite sign and can support nonsingular, geodesically complete cosmologies with ordinary NEC-respecting matter. We give explicit closed-FRW scalar-field constructions, including a fully analytic quadratic reconstruction and a cubic-root reconstruction in closed quadrature, and contrast them with their flat realizations, which require NEC-violating support. Furthermore, we quantify how positive curvature can bias flat-model reconstructions toward an effective phantom equation of state, finding only a percent-level effect under current curvature bounds. The result is a curvature classification of ANEC-compatible eternal FRW cosmology: flat and open branches are obstructed, while the closed branch admits explicit complete realizations, with global de Sitter appearing as the vacuum limiting representative.

Gravitational-wave observations of merging binary neutron stars and black holes are now routinely made by detectors in the Advanced LIGO-Virgo-KAGRA network. Neutron star binary systems may also produce detectable electromagnetic and particle emission over times scales ranging from seconds to years. Real-time gravitational-wave searches play a central role in enabling time-critical electromagnetic and/or neutrino follow-up observations. During the fourth observing run (O4) of the Advanced LIGO-Virgo-KAGRA network, multiple real-time searches operated continuously to identify candidate gravitational-wave events and publicly disseminate information about these discoveries. Here, the performance and results of the GstLAL real-time analysis are reported. The analysis is designed to identify candidates with low latency, high detection efficiency, and sustained operational uptime over long observing periods. Across O4, it produced initial candidate uploads with a median latency of 15.8 s while maintaining an effective uptime of 98% during the first two parts of the observing run. During the run, the analysis contributed to 250 candidates classified as astrophysically plausible, provided the first upload for 222 of these, and was the sole contributor for 75. Among Gravitational-Wave Transient Catalog events with a false-alarm rate below one per year, 88% were identified as significant in low latency and promoted for expert vetting and public dissemination. The low-latency astrophysical classifications agreed with the final catalog classifications for 93% of the events considered.

The vast chemical space of possible small molecules, estimated at 10^60 compounds for molecules composed of just C, N, O, and S, is only sparsely occupied by biology. We propose that where life selects molecules within this space constitutes a detectable ecological signature: a fingerprint not of specific compounds, but of the statistical structure of elemental composition across molecules sam-pled from ecological systems. Here we introduce a framework combining Van Krevelen diagrams and element scaling laws to characterize the elemental composition of regions of chemical space occupied by biological systems and contrast them with other chemical systems. Applying this framework to 11,834 microbial metagenomic samples, we show that microbial metabolisms occupy a region of chemical space, which is enriched in heteroatoms such as P, S, N, and O relative to C, shifted toward higher O:C and H:C ratios. We observe sublinear element scaling with system size, yielding insights into how elemental constraints dictate how biological systems occupy chemical space. These patterns are distinct from a sample of 18,000 compounds from the comprehensive Reaxys synthetic chemical database. Critically, datasets from molecules detected in planetary science mission data occupy statistically distinct regions from both terrestrial biological and Reaxys distributions, demonstrating that with standardized methods for data collection, the approach could be developed to discriminate biotic from abiotic chemical signatures in small molecule data from planetary science missions. Our work shows how a combination of Van Krevelen fingerprinting and elemental scaling laws can provide a new class of ecological biosignatures for life detection leveraging mass spectrometric data from planetary missions, which could generalize beyond Earth's specific biochemistry.

The gradient expansion and the separate universe approach provide an effective description of inflationary soft modes after coarse-graining shorter-wavelength degrees of freedom. We formulate a locality condition on the quantum state, requiring that the hard-mode state in each local universe depend on the soft modes only through the local soft-mode values in the same patch. When this condition is satisfied, the coarse-grained soft-mode dynamics remains local, and loop corrections from hard modes to superhorizon correlators of the adiabatic curvature perturbation are perturbatively suppressed. This provides a model-independent diagnosis of when enhanced corrections due to hard modes can invalidate the gradient expansion. We further show that the same locality condition implies a generalized soft theorem, from which the standard consistency relations follow under additional assumptions. This formulation clarifies the origin of possible deviations from the standard consistency relations in multi-field systems or in a non-attractor background. We also show that the locality condition guarantees the absence of infrared divergences for the correlators of operators invariant under a large gauge transformation. Thus, locality of the hard-mode state provides a unified criterion for the effective description of inflationary soft modes, generalized soft theorems, the suppression of hard-mode loop corrections, and the infrared regularity of observable correlators.

Owing to the emergence of attractive interactions between quarks, color superconductivity is expected to occur, with the color-flavor-locked (CFL) phase favored at high densities. This work investigates the absolute stability of beta-equilibrated CFL quark matter in bulk within the modified Nambu-Jona-Lasinio model, under color and electric charge neutrality conditions relevant to compact stars. Motivated by the possible existence of an ultra-low-mass central compact object in the supernova remnant HESS J1731-347 and the "mass-gap" secondary component in the GW190814 event, we systematically explore how vector repulsion, attractive diquark pairing, and nonperturbative vacuum effects influence the stiffness of CFL quark matter and its stability. Our findings suggest the existence of a physically viable region of parameter space in which the CFL phase is the true ground state of strongly interacting matter, thereby theoretically supporting the scenario of self-bound quark stars. This configuration is not only consistent with current astrophysical constraints from NICER and LIGO/Virgo observations, but also provides a possible explanation for both the $~2.6\ M_{\odot}$ secondary component in GW190814 and the ultra-low-mass compact object with $M = 0.77^{+0.20}_{-0.17}\ M_{\odot}$in HESS J1731-347.

Within general relativity, black holes are widely regarded as the ultimate benchmark for compactness in the Universe. Recently, however, neutron star models have been constructed in a higher-curvature theory -- quasi-topological gravity (QTG) -- whose compactness can exceed the black-hole limit~\cite{LD19666}. Here we present a detailed analysis of both the equilibrium structure and radial stability of such configurations in QTG. By examining several representative equations of state and different values of the gravitational coupling constant, we find that in the high-central-density regime the compactness exceeding the black-hole bound exhibits a universal behavior in QTG. We further show that QTG corrections grow increasingly significant at large central densities and can stabilize configurations that are radially unstable in general relativity over a broad parameter range. These results establish ultra-compact neutron stars in QTG as theoretically viable strong-field configurations and provide a foundation for further investigations of their dynamical and phenomenological implications.

We present a new accurate analysis of the $^3$He$(p,e^+\nu_e)$${}^4$He (''hep'') reaction at astrophysical energies. The S-factor is computed using a state-of-the-art method to calculate the four-nucleon scattering and bound-state wave functions (the hyperspherical harmonic expansion), and by using nuclear interactions and accompanying electroweak nuclear currents obtained within the chiral effective field theory framework. Our analysis includes a detailed examination of the theoretical uncertainties coming from two different sources: the truncation of the interaction and current chiral expansions, and the model dependence. Our recommended final theoretical value for the hep S-factor at zero energyis $S(0)=(8.7\pm 0.9)\times 10^{-20}$ keV b. We provide also the energy spectrum of the outgoing hep positrons which may be measured in future experiments. We include also an analysis of the ''sister'' reaction $^3$He$(n,\gamma)$${}^4$He (''hen'') at low energies, showing that the calculation well reproduce the total cross section from thermal energies to few MeV, validating our results on the hep reaction.

In this paper we propose a novel unified cosmological model that connects a pre-inflationary epoch, starting at the Planckian time, with the onset of inflation within a single scalar-field framework. The pre-inflationary phase is characterized by a decelerated expansion with an increasing comoving Hubble horizon, followed by a gradually transition to an accelerated inflationary regime. This early dynamics leads to a modified causal structure that naturally accounts for the suppression of large-angle $(\theta \gtrsim 60^\circ)$ correlations in the cosmic microwave background (CMB) reported by the satellite PLANCK. We study the quantum fluctuations of the scalar field using the Mukhanov-Sasaki formalism and a canonical quantization procedure based on energy minimization. We find that the vacuum state is well-defined only for sub-horizon modes at the onset of inflation, which induces a natural cutoff in the primordial power spectrum. The resulting spectrum exhibits a suppression at large scales while remaining nearly scale-invariant at small scales. In the appropriate limit, the model recovers the standard de Sitter result, in agreement with current observational constraints. These results highlight the relevance of pre-inflationary dynamics for addressing large-scale anomalies within a consistent inflationary framework.

Super-Kamiokande has reported a small excess of electron antineutrino events in the 20 MeV energy range, in the search for the diffuse supernova neutrino background. We interpret this signal as a possible indication of dark matter that annihilates dominantly into neutrinos, pointing to a thermal dark matter candidate with $s$-wave annihilation and with mass in the tens of MeV range. This mass scale naturally fits into rich dark sector extensions of the Standard Model. Neutrino experiments, including JUNO, will be able to test this hypothesis in the coming years.

In the standard thermal relic scenario, dark matter remains in chemical equilibrium with the Standard Model radiation bath until freeze-out occurs at $T \sim m_X/20$, where $m_X$ is the dark matter mass. In this familiar class of models, the observed relic density is obtained for annihilation cross sections of order $\sigma v \sim 10^{-26}$ cm$^3$/s. We show that comparable cross sections can naturally be realized in hidden-sector models in which the dark matter and Standard Model sectors decouple at a very high temperature, $T \gg m_X$, and subsequently evolve with separate thermal histories. Despite this decoupling, the two sectors have similar temperatures during freeze-out, leading to the usual thermal relic abundance. As a consequence, the coupling between the Standard Model and hidden sectors can be extremely small, potentially placing direct detection and collider signals far below foreseeable sensitivities.