2026-09-25 | | Total: 101
Although the ultraviolet background from reionization is thought to conspire against the smallest halos, some gas-rich halos in the post-reionization Universe still form stars. We ask which properties, environmental or intrinsic, determine ignition. Using the FIREbox cosmological volume, we compare recently-ignited halos to a mass- and gas-matched sample of gas-resolved starless halos at $z=2$--6. We test three families of "third axis" properties beyond $M_{\rm vir}$, $M_{\rm gas}$, and $z$: geometric environment (the Perturbation Index and local galaxy density $ρ_5$), hydrodynamic environment (the $1$--$2,R_{\rm vir}$ gas-supply shell), and intrinsic halo properties (concentration and the neutral-gas reservoir). A nested random-forest classification finds that geometric environment adds little, while the hydrodynamic shell is modestly informative. Three complementary properties each discriminate ignited from starless halos beyond the mass--gas--redshift baseline (AUC $=0.635$, where $0.5$ is chance): total HI mass ($M_{\rm HI}$, AUC $=0.739$), the extent of the neutral gas ($r_{\rm HI,50}/R_{\rm vir}$, AUC $=0.738$), and the central HI density ($ρ_{\rm HI,0}$, AUC $=0.700$); jointly they reach AUC $=0.823$, the highest of any combination we test. At fixed halo mass, gas mass and redshift, recently-ignited halos harbor more HI over a larger radius than their matched controls, and these two enhancements are uncorrelated. The same reservoir sets when ignition can happen: from $z=6$ to $z=0$ the median starless halo keeps its gas ($+0.05$ dex) while its neutral fraction falls by $3.9$ dex, so the reservoir is ionized rather than exhausted and the ignition window closes below $z\simeq2$. In sum, a halo's own neutral-gas reservoir discriminates ignited from starless halos more sharply than its location in the cosmic web -- nature over nurture.
The distribution $f(e)=2e$ is commonly referred to as the `thermal eccentricity distribution' of Keplerian binaries. However, the result does not require thermal equilibrium, and its usual derivations obscure why the distribution is linear in eccentricity. We seek a simple geometrical interpretation of this result and its relation to more general eccentricity distributions. We consider bound Kepler orbits in a Euclidean space of dimension $D\geq2$, with a phase-space distribution depending only on energy, and derive the eccentricity distribution from the Liouville measure. We obtain $f_D(e)=(D-1) e (1-e^2)^{(D-3)/2}$, so that, within this family, $D=3$ is the only case for which the distribution function is linear in $e$. This follows from the $D-1$ transverse momentum dimensions and the Kepler period degeneracy. We further show that two known generalisations -- namely a power-law weighting of the normalised angular momentum and a constant velocity anisotropy -- are, in fact, two representations of the same one-parameter family of eccentricity distributions.
The thermodynamics of protoplanetary disks, governed by a complex interplay between gas, dust, and radiation, exerts a strong influence on their morphology and dynamics at both large and small scales. Historically, hydrodynamical simulations have treated thermodynamics with approximate prescriptions, such as local isothermality, parametrized $β$-cooling, and radiation hydrodynamics with gas-dust thermal equilibrium. However, a more comprehensive and self-consistent approach would more accurately reproduce the wealth of features found in observations, which in modern times have reached unprecedented spectral, angular, and vertical resolution. With this motivation, we have devised a radiation-hydrodynamics scheme for the PLUTO code with energy exchange (via absorption, emission, and collisions) for gas, radiation, and multiple species of dust. Dust-gas dynamics are handled by treating each dust species as a pressureless, diffusive fluid, ensuring that our scheme can represent changes in disk illumination caused by grain settling and trapping. We demonstrate the effectiveness of our scheme in several test problems, and conclude by discussing its relevance to open questions in protoplanetary disk studies.
We describe a probabilistic approach to gravitational clustering. Working purely in configuration space, we show that the two-point function can be formulated in terms of a transition probability describing the likely final position of a particle given its initial position. The transition probability is initially sharply peaked and then over time broadens out in a diffusion-like manner. The effect is more pronounced for widely-separated particles (e.g.\ at the baryon acoustic scale), leading to an effective suppression of the two-point function in what is normally referred to as infrared resummation.
Integral field spectroscopy (IFS) of the central kiloparsecs of active galactic nuclei (AGN) reveals a mixture of spatially coincident emission from star-formation and AGN feedback exciting the interstellar medium. Disentangling these components remains a challenge, as most spectral tracers are affected by both processes, limiting robust interpretation of kinematics and energetics. We present a new framework for decomposing IFS data into distinct components on a spaxel-by-spaxel basis using a multi-tracer, stacked kinematics approach. This method combines kinematic modelling with imposed flux decomposition per spaxel, quantifying the contribution of each component across the field of view. We apply this method to JWST IFS observations of the Seyfert galaxy NGC 5728 from the Galaxy Activity, Torus, and Outflow Survey (GATOS), analysing twelve mid-infrared fine-structure lines (4.49 < $λ$ < 25.89 $μ$m, 7.9 < IP < 126.2 eV). We find that the circumnuclear emission can be decomposed into two dominant components: a star-forming ring and an AGN-driven biconical outflow. Our method separates these structures and recovers their detailed spatial morphology. This framework provides a general and scalable method for physically motivated component separation in IFS data, applicable across many wavelength ranges and targets, enabling reliable interpretation of complex emission line structures (or morphologies) in active galaxies and beyond.
A large fraction of stars interact with a close companion during their lifetime, in which the transfer of mass and angular momentum shapes their evolution and final fate. Standard rotationally limited accretion models predict that accretors reach critical rotation after gaining only a small fraction of their mass, severely suppressing further accretion. This is in tension with observations of post-interaction systems that require substantial mass gain. We introduce a disc-mediated angular momentum transport prescription for mass-transfer onto stellar companions. This is based on a novel perturbative, analytic star-disc boundary model that allows for continued mass inflow as the accretor approaches critical rotation. Furthermore, the model reproduces previous numerical results in which perturbations to super-critical rotation result in negative torques exerted by the disc, extracting excess angular momentum while allowing continued mass inflow. We implement this mechanism in detailed binary evolution calculations with differential rotation using MESA, and compute grids spanning primary mass, orbital period, and mass ratio. Whereas rotationally limited models predict $β_{\rm eff}\lesssim 0.1$, the disc model yields sustained mass inflow near critical rotation, with effective mass-transfer efficiencies of $β_{\rm eff}\sim 0.4$-$1$ across much of the parameter space. The resulting accretor properties are broadly consistent with observed post-interaction sdOB+Be binaries. Disc-mediated angular momentum transport may therefore represent a key missing ingredient in standard binary evolution models, with important implications for rapidly rotating stars and compact-object progenitors.
High-performance optical systems, including space telescopes, require wavefront-sensing and correction capabilities to achieve resolution near the diffraction limit. Existing wavefront-sensing approaches span a broad range of architectures, but many require dedicated interferometric hardware, rely on computationally intensive nonlinear or iterative reconstruction algorithms or provide only limited wavefront information. Focal-plane wavefront-sensing methods are particularly attractive because they can operate without dedicated metrology hardware, and among these, Double Blind Fourier Holography (DBFH) stands out as a robust and computationally light alternative to nonlinear and iterative solvers. Here we introduce differential Double Blind Fourier Holography (diff-DBFH), a two-image variant of DBFH based on a small pupil-plane amplitude perturbation. Inspired by differential Optical Transfer Function (dOTF), diff-DBFH leverages the linear framework of DBFH to produce a more precise wavefront estimate. We derive the diff-DBFH formalism, compare its performance against dOTF using numerical simulations on a segmented aperture geometry inspired by the Habitable Worlds Observatory, analyze its sensitivity under noisy conditions, and demonstrate the method experimentally using an optical bench model of the same aperture.
We report spectroscopic observations during a prominence eruption on 9 October 2025 from the Sun's west limb, associated with an M2.0 flare that commenced at 12:11 UT. The event was captured in sit-and-stare mode by the Visible Emission Line Coronagraph (VELC) onboard Aditya-L1, providing continuous observations of the Fe XIV 5303 Åcoronal greenline. The spatially averaged green-line emission exhibits a clear correspondence with the EUV late-phase (ELP) emission observed in multiple SDO/AIA channels. In particular, an enhanced 5303 Åemission peak at 15:24 UT coincides with secondary peaks in the AIA 211, 171, 193, and 335 Åpassbands, occurring $\sim173$ minutes after the X-ray peak. The 5303 Åline-width information provides important constraints on the mechanism responsible for the ELP. The line-width decreases progressively by about 28\% toward the ELP peak, broadly consistent with the cooling of post-flare coronal loops. This behavior suggests that the ELP is predominantly associated with plasma cooling, with relatively weak additional energy input from magnetic reconnection, which does not produce a significant enhancement in line-widths. The estimated non-thermal velocity decreases from $\sim$36 km s$^{-1}$ during the impulsive phase to $\sim$21 km s$^{-1}$ during the gradual phase, indicating a reduction in turbulent plasma motions. And doppler velocities evolve from predominantly positive values during the eruptive rise phase to negative values during the later gradual phase, implying upward and subsequent downward plasma motions. To our knowledge, this study provides the first space-based spectroscopic evidence linking coronal green-line emission with the ELP and its origin during a solar flare.
JWST has increased both the number of high-redshift galaxies with high-quality spectral energy distributions (SEDs) and their information content. In parallel, wide-area surveys from Euclid, Rubin's LSST, and Roman will increase galaxy samples by orders of magnitude. Analysing these datasets requires stellar-population models that are both flexible and computationally efficient. We present CERIDWEN, a GPU-native SED fitting framework written in JAX, with an end-to-end differentiable forward model spanning stellar populations, nebular emission, dust attenuation and emission, and projection into the observer frame. Its vectorised, compiled architecture lets nested sampling replace a batch of live points in parallel on the GPU, which makes flexible stellar-population models tractable under full Bayesian inference. Automatic differentiation also provides exact gradients for the gradient-based samplers in the package. We jointly infer time-dependent chemical-enrichment histories instead of a single stellar metallicity, and demonstrate non-parametric star-formation histories (SFHs) with $\sim$120 age bins. Using $α$-enhanced stellar libraries from FSPS, CERIDWEN can sample stellar [$α$/Fe] jointly with [Fe/H], mass, and SFH, so that the joint posterior represents the [Fe/H]-[$α$/Fe] degeneracy explicitly. In controlled mocks, CERIDWEN recovers parameters with well-calibrated posterior uncertainties, while fits to real JWST observations reproduce posteriors from the established Prospector framework: on a single GPU, CERIDWEN completes a fit in a median sampling time of $\sim$4 min, $\sim$134$\times$ faster per fit than equivalent CPU-based Prospector runs. CERIDWEN therefore makes full Bayesian inference practical for larger galaxy samples and more flexible stellar-population models, reducing computational constraints on the physical complexity explored in SED fitting.
The Mg-bearing cations MgC$_4$H$^+$, MgC$_6$H$^+$, MgC$_3$N$^+$, and MgC$_5$N$^+$ have been recently observed in the carbon-rich envelope IRC\,+10216. These species are thought to form upon radiative association between Mg$^+$ and the corresponding neutral radical. This work aims to determine the radiative association rate coefficients for the cations Al$^+$, Mg$^+$, and Na$^+$ reacting with carbon chain radicals to estimate the relative importance of metal-bearing carbon cations for each of these metals. For this purpose we use statistical methods in combination with highly accurate ab initio calculations to obtain the radiative association rate coefficient. Moreover, anisotropic effects beyond mono-dimensional capture model are considered, providing a more reliable description of the simulated processes.} The radiative association rate coefficients to form Mg-, Al- and Na-C$_6$H$^+$, as well as MgC$_4$H$^+$, MgC$_3$N$^+$, and MgC$_5$N$^+$, are calculated with the new restrictions considered in this work. These new rate coefficients are included in a chemical model of the C-rich AGB envelope IRC\,+10216. The abundances of MgC$_4$H$^+$, MgC$_6$H$^+$, and MgC$_5$N$^+$ are consistent with observations, with differences no greater than one order of magnitude, corroborating radiative association as a plausible formation mechanism for these adducts. However, the calculated abundance of MgC$_3$N$^+$ is about two orders of magnitude lower than observed, which indicates that either the chemistry of this species is significantly different from the others or the calculated rate coefficient is too low. The abundances obtained for Na-C$_6$H$^+$ and Al-C$_6$H$^+$ are around two orders of magnitude lower than for Mg, making very difficult to detect the analogous metal-bearing cations in IRC\,+10216.
The $β$ Pictoris debris disk is a unique system where mid-infrared coronagraphic imaging with JWST/MIRI revealed a new spatially resolved substructure known as the Cat's Tail. We present deep HST/STIS coronagraphic imaging of $β$ Pic in scattered light, combining six epochs of observations between 2012 and 2025, in search of the scattered light component of the Cat's Tail. We detect the disk out to a projected separation of 500 au for the first time with HST/STIS, achieving SNR > 100 across much of the midplane between 50 and 200 au. We report the detection of the visible/near-infrared scattered light component of the Cat's Tail substructure. We find through injection-recovery analysis that the scattered light component of the Cat's Tail has a total flux of $0.1-0.4\%$ of the MIRI F1550C flux. Using the measured STIS-to-MIRI flux ratio, we model the Cat's Tail grain properties and find that the grains are highly porous and almost entirely composed of organic refractory material, which is in agreement with findings with JWST/MIRI. We use the organics content of Solar System dwarf planets and dust in the local ISM as a proxy to estimate the size of the colliding progenitors that can produce enough organic refractory material seen in the Cat's Tail. We find that each colliding progenitor must have a mass of at least $1-3\times 10^{21}$ kg, comparable to Charon and Makemake in our Kuiper Belt. Additionally, we find a prominent warped morphology and surface brightness asymmetry in the outer dust halo of $β$ Pic. We compare the observed halo morphology and asymmetry to predicted vertical structures, which may arise from planet--disk interactions.
Within the geometric optics approximation, the optical expansion scalar describes the local rate of change in the amplitude of a wave or the area of a light bundle propagating through the Universe. As such, it is fundamental for describing almost all the observations we make in cosmology. It was recently pointed out that this quantity can, in principle, be reconstructed from measurements of the angular diameter distance $D_A$ and the effective observed expansion rate $\mathfrak{H}$. Motivated by this possibility, we introduce the observable $\vartheta \equiv θ/ω_o$, where $θ$ is Sachs' optical expansion scalar, and $ω_o$ the observed photon frequency. We then derive the linear perturbative expression and angular power spectrum for $\vartheta$ in a flat FLRW universe with scalar perturbations. Comparing the power spectra with those of the fluctuations in $D_A$ and with power spectra from relativistic N-body simulations reveals no immediate striking advantage of studying the fluctuations of $\vartheta$ to those of $D_A$. However, the anisotropies of $\vartheta$ do exhibit several unique qualities such as an enhanced sensitivity to the peculiar velocity and relativistic potentials. We suggest that our results may hint at $\vartheta$ being useful for developing consistency tests of cosmological observations and possibly probing departures from FLRW cosmology and modified gravity theories. The formalism developed and presented here provides the foundation for future studies into these possibilities.
Relativistic jets in blazars exhibit strong variability across multiple timescales, but the connection between temporal variability and the evolution of jet structures remains poorly understood. In this work, we investigate the $γ$-ray variability of the BL Lac object S5 0716+714 from July 2014 to June 2015 and explore its possible connection with the parsec-scale jet evolution. We identify an escalating periodic oscillation (EPO) in the $γ$-ray light curve, characterized by a sequence of successive flux enhancements with progressively increasing temporal separations. Combining \textit{Fermi}-LAT data with contemporaneous 43 GHz Very Long Baseline Array (VLBA) observations, we find a temporal association between the EPO and the downstream propagation of a superluminal knot. A helical trajectory model simultaneously fitting the $γ$-ray variability and the two-dimensional VLBA component positions reproduces the observed flux modulation and component motion. These results suggest that geometric Doppler modulation caused by the changing viewing angle of a moving jet component provides a possible explanation for the observed EPO. This work provides a framework of connecting blazar variability to the knot evolution of relativistic jets.
Solar vortices are rotating plasma structures that play a key role in the transport of energy and magnetic helicity from the photosphere to the chromosphere and corona. While numerical simulations predict the existence of small-scale photospheric vortices with diameters of $50$--$100$ km, no direct observational evidence has been reported to date, primarily because of the limited spatial resolution of previous instrumentation. Here, we analyze high-resolution observations obtained with the 4-meter Daniel K. Inouye Solar Telescope (DKIST) at 416 nm, with a spatial resolution of 12.3 km. Using Fourier local correlation tracking and the swirling strength criterion ($λ_{\rm ci}$), we identify 201 vortex detections across a sequence of 143 reconstructed frames, spanning approximately 3.2 minutes, within a field of view of $4.92 \times 4.13$ Mm. The observations cover a photospheric plage region containing pore-like magnetic concentrations. The detected vortices have a median equivalent diameter of 38.6 km and a mean diameter of $39.5 \pm 3.3$ km. Virtually all detections are smaller than 50 km and therefore fall within the subgranular regime predicted by magnetohydrodynamic simulations but never previously observed. The median swirling period, $τ_{\rm ci}=37.3$ s, satisfies the $τ_{\rm ci}<100$ s criterion established by MURaM simulations. No statistically significant preference is found for prograde rotation (53.7\%) over retrograde rotation (46.3\%), consistent with the negligible influence of the Coriolis force at subgranular scales. These results demonstrate the unique capability of DKIST to resolve the population of photospheric convective vortices predicted by state-of-the-art magnetoconvection simulations, opening a new observational window into small-scale solar convective dynamics.
We aim to constrain the dynamical masses of the four confirmed transiting planets orbiting the infant ($20\pm10$ Myr) star V1298,Tau (orbital periods $P_b\sim$24.1 d, $P_c\sim$8.2 d, $P_d\sim$12.4 d, and $P_e\sim$48.7 d), through a new analysis of a homogeneous radial velocity dataset. This system is a key benchmark for testing formation and early-stage evolutionary models of planets. We analysed more than 380 high-resolution spectra collected with the HARPS-N spectrograph over a time span of 4.5,yr. To account for dominant stellar activity and transit-timing-derived planet-planet interactions, we analysed radial velocities extracted using two different approaches, and combined N-body simulations with multi-dimensional Gaussian process regression. We tested a set of spectroscopic and photometric activity diagnostics to mitigate the variability due to the stellar magnetic activity in the radial velocity time series. The mass determinations are sensitive to the choice of the radial velocity extraction method and activity diagnostics. For some specific models, we measured masses with a significance of 3.5--4$σ$ for planets b ($0.52^{+0.12}{-0.14}$, $0.65^{+0.21}{-0.18}$, $0.67\pm0.16$, and $0.78^{+0.19}{-0.20}$ $M_{\rm Jup}$) and e (0.40$^{+0.14}{-0.13}$, 0.66$^{+0.17}{-0.18}$, and $0.71\pm0.20$ $M_{\rm Jup}$). These values are compatible with masses of giant planets, and are notably higher than those recently derived from transit timing variations alone, which indicate that these are low-density ``super-puffy'' planets. Injection-and-recovery simulations moderately support this conclusion for planet b. For planets c and d, we primarily provide updated mass upper limits. Our analysis represents an exploratory study to address the complex issue of mass measurement in a multi-planetary system characterised by high levels of stellar activity.
We present a JWST mid-infrared (MIR) view of Pox 186, one of the best local analogs of reionization-era galaxies based on its extreme properties in the ultraviolet, optical, and far-infrared (FIR). We extracted a high signal-to-noise spectrum of Pox 186 from the spatially resolved MIR data taken with the Medium Resolution Spectrometer (MRS) on the MIR Instrument (MIRI), and compared its MIR properties with those of the other dwarf galaxies in the Herschel Dwarf Galaxy Survey. We probe the radiation hardness of Pox 186 using the so-called softness diagram, based on the MIR line ratios [NeIII]/[NeII] and [SIV]/[SIII], which suggest that Pox 186 hosts ionization sources with a hard radiation field. We probe the dominant ionization source in Pox 186 using diagnostic diagrams based on MIR line ratios. Our photoionization models provide evidence that stripped binary stars are the main cause of the extreme MIR line ratios. Similarly, a young starburst is needed to reproduce the MIR dust continuum. We probe the warm and cold molecular gas content using the MIR H$_2$ rotational lines (S(1)-S(7)) from JWST, and the archival sub-mm CO(2-1) line from ALMA, which indicate a high star-formation efficiency (SFE; $>$ 15 $\times$ 10$^{-9}$ yr$^{-1}$) and low gas depletion time comparable to that of the reionization-era galaxies. This analysis strongly suggests that star formation rather than active galactic nuclei drives the hard radiation field observed in Pox 186, and cautions against relying solely on the UV diagnostics to establish the nature of high-z galaxies.
We infer the Hubble constant $H_0$ from the ~400 Extragalactic Distance Database hosts with a tip of the red giant branch (TRGB) distance, calibrated on the geometric anchors of the Large Magellanic Cloud and NGC 4258. We forward-model the tip magnitudes and host redshifts including magnitude-limited selection, sky exposure, and density-dependent galaxy bias, conditioning on each of the 80 Manticore-Local realisations of the local density and peculiar-velocity fields in turn. As most hosts lie at $cz_{\rm CMB} \lesssim 1000$ km/s, where peculiar velocities are important, the inference is limited by the variance of these fields rather than by the TRGB calibration. For a fixed realisation, the redshift likelihood and density smoothing each move $H_0$ by less than 1 km/s/Mpc and a velocity monopole by $5.4 \pm 3.4$ km/s/Mpc, whereas the field-to-field scatter in $H_0$ is 5.2 km/s/Mpc. Weighting the posteriors from each realisation equally gives $H_0 = 68.4 \pm 5.9$ km/s/Mpc, while weighting each by its Bayesian evidence gives $72.3 \pm 2.6$ km/s/Mpc. The latter value is set by a single particularly high-evidence realisation, making it sensitive to the finite sampling of the reconstruction posterior. This ladder therefore cannot yet arbitrate the Hubble tension. Extending the TRGB data beyond the present ~1000 km/s, together with improved reconstructions of the local Universe, would reduce the field-to-field scatter well below the 5.2 km/s/Mpc we find here.
The physical conditions of the interstellar medium (ISM) are thought to differ greatly at high redshifts compared to local galaxies. Polycyclic aromatic hydrocarbons (PAHs) are key ISM constituents, mediating the heating and cooling of gas in photo-dissociation regions. However, PAHs are unexplored in main-sequence (MS) galaxies at z~1-2 (cosmic noon), limiting our view of the ISM's thermal and chemical balance. We present the PAH Spectroscopic Survey (PAHSPECS) JWST program, the first MIRI/MRS survey of the 3.3, 6.2, 7.7, 8.6 and 11.3$μ$m PAH features in MS galaxies at z~1. The five targets come from the deep ALMA CO and continuum flux-limited sample in the HUDF. We detect the 6.2, 7.7 and 11.3$μ$m PAHs in all galaxies, and the 3.3$μ$m feature in only two. Spectral energy distribution (SED) fitting to broad-band mid-IR photometry recovers the 6.2 and 7.7$μ$m PAH luminosities derived from spectral decomposition of the MRS spectra to within ~25% accuracy. However, the 3.3 and 11.3$μ$m features are overestimated by ~50% (up to 3-5x). Spatially resolved maps of ASPECS-6, the most extended source, show variations of the 7.7$μ$m PAH contribution to the MIRI F1500W flux, suggesting that while photometry can infer PAH luminosity accurately, IFU spectroscopy is necessary to study variations in PAH ratios and characterize galaxy dust properties. The PAH fraction of the total dust mass from photometric SED modeling roughly agrees with the $L_{\rm PAH}/L_{\rm IR}$ ratio measured from spectroscopy, but the correlation is not as tight as models predict. On average, PAHSPECS galaxies show an excess in 6.2$μ$m equivalent width relative to nearby dusty luminous systems, but lower 3.3$μ$m equivalent widths than local systems, similar to higher-mass starbursts at z~1-3. This points to differential changes in the size distributions of neutral and ionized PAHs in cosmic noon galaxies.
S4 0954+65 (4FGL J0958.7+6534) is a TeV-detected blazar at a redshift of $z = 0.3694 \pm 0.0011$, classified as an intermediate-synchrotron-peaked BL Lac object. In early 2025, it entered an exceptional $γ$-ray high state. We aim to investigate the origin and physical properties of the exceptional 2025 flare, and to constrain the emission processes responsible for this flaring activity. We performed a multi-wavelength analysis using $γ$-ray, X-ray, optical/UV, radio, and 43 GHz Very Long Baseline Array (VLBA) observations. We examined multi-band correlations with the z-transformed discrete correlation function (zDCF), analyzed the spectral evolution and parsec-scale jet kinematics, and modeled the broadband spectral energy distributions (SEDs). The $γ$-ray variations lead the optical and radio emission by 3.22 days and possibly 18.26 days, respectively, while no significant correlation is found between the $γ$-ray and X-ray emission. The $γ$-ray spectra show a harder-when-brighter behavior, consistent with enhanced particle acceleration during the active state. A similar spectral trend is also observed in the X-ray band, where the photon index is anti-correlated with flux. During the major $γ$-ray flare, VLBA images reveal the emergence of a new superluminal radio knot from the core, whose extrapolated ejection time is consistent with the peak of the $γ$-ray flare. The temporal and structural evolution during the 2025 flare is consistent with a shock-in-jet scenario, in which a newly emerging disturbance propagates downstream and drives the flare. Broadband SED modeling shows that a leptonic synchrotron self-Compton plus external Compton model with dusty torus seed photons reproduces the multi-wavelength observations with physically plausible parameters. A lepto-hadronic interpretation remains possible but requires substantially higher jet power.
Standard point spread function (PSF) matching algorithms consider one image as the reference and aim to find a convolution kernel such that, when applied to the reference image, the PSFs of the reference and non-reference images are matched. This method works well if one image is a degraded version of the other. However, this is not always the case: if the PSFs of the two images have different orientations or shapes, a simple convolution of one of the images does not lead to a satisfactory match. While this can be trivially solved by convolving the first PSF with the second and vice-versa, this results in unnecessarily large PSFs. In this paper, we present a new algorithm, Dual Kernel Optimization (DKO), that simultaneously solves for two convolution kernels, one for each image, such that the PSFs match and the sizes of the kernels are minimized. This results in a minimally-degraded final set of PSFs. Since this problem is highly degenerate, we use Stochastic Gradient Langevin Dynamics (SGLD) which rigorously explores the parameter space and converges on the globally-optimal solution.We discuss the algorithmic challenges for this problem and the modifications to standard SGLD used to constrain it. Finally, we include a discussion on the application of this methodology to real-world images, including cases where the PSF varies over the field of view.
We report the discovery of J1351+0039, a new isolated dark galaxy candidate identified through a blind HI survey with the Five-hundred-meter Aperture Spherical radio Telescope (FAST), combined with deep optical imaging from the DESI Legacy Survey. The object displays a compact HI structure with a systemic velocity of V_{LSR} = 429.0 \kms and an inferred HI mass of 4.8*10^5 Msun. No optical counterpart is detected in the DESI images, with a g-band surface brightness limit of 28.1 mag arcsec^{-2}, corresponding to a stellar mass limit of M_{star} of 1.5*10^{5} Msun. The dynamical mass is estimated to be M_{dyn} of 1.3*10^{8} Msun, yielding a dark-to-baryonic mass ratio of 200, implying that J1351+0039 is overwhelmingly dominated by dark matter. While week rotation is detected in J1351+0039, the kinematics are dominated by internal velocity dispersion, confirming pressure support and reminiscent of Reionization-Limited HI Clouds (RELHICs). J1351+0039 stands out as one of the most promising isolated RELHIC candidates, offering a unique laboratory for studying baryonic physics and dark matter in the absence of star formation.
The Vera C. Rubin Observatory will enable large-scale characterization of minor bodies through high-cadence, multi-band photometry. Early Science Validation (SV) observations from 2025-2026 provide one of the first extensive datasets for Solar System studies, but typically with sparser and more heterogeneous coverage than the Rubin First Look data, where hundreds of observations were available for well-characterized objects. In this work, we build on the methodology developed for the Rubin First Look data to extract asteroid rotation periods, amplitudes, colors, and taxonomic classifications from early Rubin SV data using multi-band (griz) photometry. Rotation periods are determined through a combination of high-order Fourier light-curve fitting and multi-band Lomb-Scargle analysis. A unified light-curve model is then used to derive rotation-corrected colors, which are converted into taxonomic classes using calibrated color indices and spectral slopes. We validate the method on objects with sufficient observational coverage and show that preliminary physical properties can be estimated from sparse and irregular datasets, while robust determinations require improved temporal coverage. By simulating sparsely observed data, we identify the observational conditions required for period determinations. The inclusion of z-band photometry extends traditional color-based taxonomy and improves compositional discrimination relative to gri-only analyses.
We present a multiwavelength analysis of the star $ρ$ Oph C, covering X-ray, optical, near-infrared, and radio observations to test for the presence of interferometrically-detected cool star companions. The X-ray observations from Chandra, XMM-Newton, and NuSTAR show flare-like events lasting minutes or days in time and spectral properties consistent with those of cooler stars instead of the primary magnetic B star. The near-infrared data is also consistent with the proposed trinary system properties due to the presence of CO-band heads. After removing the primary B star's rotation, the optical data from TESS shows a residual signal at the level of 0.8 mmags, matching the variations expected from cool stars in orbit around a bright B star. Finally, multi-epoch radio data reveals that $ρ$ Oph C exhibits significant large scale flux variations that are atypical for a magnetic massive star. Taken as a whole, the multi-epoch and -wavelength data presents a consistent picture of $ρ$\,Oph\,C as a rare trinary system composed of a hot star and two cool stars.
We analyze the possible observational signatures and evolution of global warps in the accretion flows produced in tidal disruption events. These warps are driven by Lense-Thirring torques resulting from the misalignment of the disk's rotation axis and the black hole's spin axis. In a previous paper we argued that these global warps should be generically present in TDE accretion flows as they enter a late-time thin disk phase of evolution. We isolate three possible observational signatures of warps: (i) anomalous (reddened) late time optical colors, resulting from irradiative heating of the warped outer disk by the inner flow, (ii) delayed X-ray rises, resulting from the inner disk being slowly revealed to a distant observer as the warped flow gradually aligns with the spin axis, and (iii) reverberation of turbulent X-ray flares into optical/UV frequencies as the inner disk heats the outer disk during a flaring state. We present preliminary observational evidence that all three effects may have already been observed in TDE systems. Observations of TDEs may well offer the cleanest astronomical signatures of the strong-gravity effects which cause disk warping, and provide unique constraints on the poorly understood evolution of warped accretion disks.
Context: Hot exozodiacal dust has been detected around approximately 20% of nearby main-sequence stars, indicating that it is a common feature of extrasolar planetary systems. However, its spatial distribution and grain properties remain poorly understood. Aims: We aim to show that the spatial distribution of hot dust around beta Pictoris, an iconic planetary system with two known planets and a population of exocomets, can be constrained through a detailed analysis of optical long-baseline interferometry data. Methods: We obtained L-band observations at approximately 3.5 microns of beta Pictoris with MATISSE at the VLTI. We developed a method to assess correlations in the visibility measurements by estimating the full data covariance matrix, and modeled the circumstellar emission with geometric brightness distributions. Results: We detect circumstellar emission around beta Pictoris consistent with hot exozodiacal dust. Modeling the emission with and without accounting for data correlations demonstrates their significant impact on the inferred dust properties. The best-fit model reveals a compact (approximately 0.05 au), highly inclined dust structure near the sublimation radius, misaligned by approximately 60 degrees on the sky with respect to the extended edge-on outer disk. It contributes approximately 4%-7.5% of the total flux at 3.4 microns and is composed of submicron-sized grains. The data are also compatible with an additional dust component near 1 au. Conclusions: For beta Pictoris, the two-decade-old paradigm placing hot exozodiacal dust near the sublimation radius and consisting of submicron-sized grains finds observational support in the detailed analysis of MATISSE data. The misalignment with the known debris disk further supports a link to exocometary activity.