Papers for Tuesday, Mar 10 2026

Current solar flare predictions often lack precise quantification of their reliability, resulting in frequent false alarms, particularly when dealing with datasets skewed towards extreme events. To improve the trustworthiness of space weather forecasting, it is crucial to establish confidence intervals for model predictions. Conformal prediction, a machine learning framework, presents a promising avenue for this purpose by constructing prediction intervals that ensure valid coverage in finite samples without making assumptions about the underlying data distribution. In this study, we explore the application of conformal prediction to regression tasks in space weather forecasting. Specifically, we implement full-disk solar flare prediction using images created from magnetic field maps and adapt four pre-trained deep learning models to incorporate three distinct methods for constructing confidence intervals: conformal prediction, quantile regression, and conformalized quantile regression. Our experiments demonstrate that conformalized quantile regression achieves higher coverage rates and more favorable average interval lengths compared to alternative methods, underscoring its effectiveness in enhancing the reliability of solar weather forecasting models.

With rapid improvements in the assembly of large samples of galaxy clusters, we are approaching the ability to study clusters at $z\gtrsim2$. Evolutionary studies comparing these distant clusters to the clusters in our local universe depend heavily on the reliability of low-redshift cluster samples, most of which are subject to X-ray selection effects, biasing them to relaxed, cool core clusters. Here, we introduce the Cluster Evolutionary Reference Ensemble At Low-$z$ (CEREAL) sample, composed of Chandra X-ray observations of 169 galaxy clusters that have been selected from the Planck Sunyaev-Zel'dovich catalog. CEREAL has a simple and well-understood selection function, spans an order of magnitude in mass at $z\sim0.15$, and has uniform, high-resolution X-ray follow-up. We present the full sample and provide results based on X-ray surface brightness properties, finding significantly more non-cool core systems than in X-ray-selected samples. We use surface brightness concentration (c$_\mathrm{SB}$) as a proxy for cool core strength and centroid shift ($w$) to measure dynamical state. Over the full sample, we find a cool core (c$_\mathrm{SB} > 0.075$) fraction of $0.39_{-0.04}^{+0.04}$, a strong cool core (c$_\mathrm{SB} > 0.155$) fraction of $0.13_{-0.03}^{+0.03}$, and a dynamically relaxed ($w<0.01$) fraction of $0.42_{-0.04}^{+0.04}$. We find no mass dependence in the fraction of clusters that appear relaxed or have cool cores. We quantify the rarity of X-ray-bright central point sources (L$_\mathrm{nuc,~2-10~keV} > 10^{43}$ erg s$^{-1}$), finding them to be intrinsically rare ($0.7_{-0.5}^{+1.2}$\% of massive, low-z clusters) with a notable increase in occurrence rate at the centers of cool cores.

We study the neutrino distributions that arise in a simulation of a neutron star merger that uses a Monte Carlo (MC) neutrino transport scheme. In a snapshot taken 1 ms after merger, we calculate relevant observables to test when neutrinos behave like a thermalized gas, and when a free-streaming picture is more appropriate. We find that in hot, dense regions where neutrino-matter interactions are frequent, MC neutrino and antineutrino distributions are consistent with thermalized neutrinos. In moderately warm regions, where neither approximation is expected to hold, we find significant departures from the predictions of the thermalized-neutrino approximation, particularly for the (anti)neutrino average opacity and net rate of absorption per baryon, even when average energies appear approximately thermal. At lower temperatures, MC results approach the free-streaming limit. Our results demonstrate that energy-averaged agreement with thermalized-neutrino assumptions does not guarantee accurate weak interaction rates. Non-equilibrium aspects of the neutrino distribution are therefore crucial for neutrino-mediated microphysics such as composition evolution in the early post-merger phase.

Although restoration of chiral symmetry is predicted by quantum chromodynamics to take place at high baryon density, most modeling of neutron star interiors disregards a chiral phase transition. We model neutron star cores with a parity doublet model, which allows for dynamical chiral symmetry restoration and predicts the appearance of the parity partners of nucleons and hyperons at large densities, as well as deconfined quark matter. We study the thermal evolution of neutron stars, focusing for the first time on the impact of Urca processes involving the parity partners in neutron star cooling simulations. We find that Urca processes for the parity partners of the nucleons significantly affect the thermal evolution of massive stars and allow for improved agreement with observed surface temperature and ages.

Molecular spectral cubes of prestellar cores encode the information on the physical and chemical properties of these objects along the line of sight. To retrieve this information, we need an interpretable model that reproduces the observed spectra. We designed a differentiable 3D geometrical model that produces synthetic observations from the parameterized density and velocity fields, and that can be efficiently optimized to reproduce the real data cubes. The model has been applied to p-NH2D and N2D+ spectral cubes in the prestellar core L1544. The optimized model suggests that to reproduce the observed velocity difference between p-NH2D and N2D+ in L1544, an asymmetric structure in density and velocity is necessary.

We study structure formation in alternative cosmological models constrained by background observations, including $\Lambda$CDM, wCDM, the Chevallier-Polarski-Linder parametrisation and a flexible Chebyshev expansion of the dark energy equation of state. The models are constrained using baryon acoustic oscillations, cosmic microwave background, cosmic chronometers and strong lensing measurements. Using the best-fitting parameters, we generate cosmology-dependent initial conditions and perform N-body simulations to analyse the matter power spectrum, halo mass function and halo density profiles. Although all models remain broadly consistent with $\Lambda$CDM at the background level, differences in the physical matter density $\Omega_{0m}h^2$ and in the expansion history $H(z)$ lead to distinct growth histories that are amplified by non-linear evolution. We find a clear hierarchy in the power spectrum amplitude and in $\sigma_8$, with the Chebyshev and CPL models exhibiting enhanced small-scale power, earlier halo formation at $z\gtrsim2$ and a migration of excess toward higher masses at late times. The wCDM model displays milder and partially compensating effects driven by its different expansion history. When expressed in terms of the scaled radius $r/R_{200c}$, halo density profiles show a high degree of universality across cosmologies, indicating that internal halo structure is largely governed by the same gravitational dynamics. These results demonstrate that even modest background-level variations in $w(z)$ can translate into coherent non-linear signatures, highlighting the constraining power of large-scale structure observables in extended dark energy models.

We present the first constraints on multipole-dependent cosmic birefringence using CMB polarization data from the BK18 dataset, which combines observations from BICEP2, Keck Array, and BICEP3 at frequencies of 95, 150, and 220 GHz. Photon coupling to an axion-like field leads to the rotation of CMB polarization, inducing non-zero EB cross-correlations. We show that a multipole-dependent rotation beta(l) imprints a distinct signature in the polarization spectra that can be constrained. Specifically, we consider an Early Dark Energy (EDE) scenario in which a pseudoscalar field couples to photons through a Chern-Simons interaction, generating a polarization rotation with multipole dependence. We introduce a phenomenological beta(l) as a step function, obtaining constraints on the step function size consistent with zero, with uncertainties less than 0.15 degrees (68% CL). In addition, using multi-frequency EE, BB, and EB cross-spectra, along with robust BICEP/Keck foreground treatment and likelihood framework, we derive constraints on the axion-photon coupling amplitude g for several choices of EDE parameters. For the baseline best-fit value f_EDE = 0.087 from the Planck 2018 analysis, we obtain g = 0.11 +/- 0.37 (68% CL), consistent with previous limits.

Galaxy clusters grow through the accretion of galaxies from groups, filaments, and other clusters. During this process, galaxies may undergo pre-processing in lower-density environments, where galaxy-galaxy mergers and other interactions can significantly alter their properties prior to cluster infall. We investigate the role of galaxy mergers in the pre-processing of galaxies prior to cluster infall by studying the spatial distribution of mergers across the cosmic web. We use a sample of 43,922 galaxies targeted by the 4MOST CHANCES survey in and around 33 low-redshift clusters (z < 0.07). Using Zoobot, a deep-learning framework trained on Galaxy Zoo data, we identify 698 galaxy mergers. We measure their distances to cosmic web filaments and compare them with those of non-merging galaxies. We find that galaxy mergers are significantly closer to filaments than the non-merging galaxy population, with this trend being strongest beyond the cluster virial radius. This suggests that filaments provide conditions conducive to mergers, possibly moderating relative velocities and enhancing gas availability. Our findings support a scenario in which filaments play a key role in transforming galaxies through pre-processing by promoting mergers before they enter cluster cores where star formation quenches.

The relationship between magnetic field strength and gas density is essential to understand the interstellar medium and star formation. Zeeman measurements in dense atomic and molecular gas phases have traditionally been used to directly probe magnetic field strengths in the Milky Way. This allowed derivation of a relationship between magnetic field strength $B$ and gas number density $n$. We recently generalized this relation as a two-part power-law with non-zero slopes and a transition density given as $B/B_0 \propto (n/n_0)^{\alpha_1}$ for $n \le n_0$ and $(n/n_0)^{\alpha_2}$ for $n > n_0$. Here, we extend our previous hierarchical Bayesian framework by incorporating a large body of pulsar observations that probe the diffuse interstellar medium and explicitly modelling density uncertainties through a global log-density correction parameter $R$ applied to all densities. We also account for magnetic field geometry and measurement uncertainties through a magnetic hyperparameter to estimate $B$. This results in a stronger constraint on the diffuse gas part of the $B$--$n$ relation. Our results confirm a non-zero exponent in the diffuse gas and a broad transition density with our best model and data set yielding maximum a posteriori results of $\alpha_1 = 0.18^{+0.02}_{-0.02}$, $\alpha_2 = 0.63^{+0.08}_{-0.05}$, $n_0 = 1630^{+2560}_{-1430}\,\text{cm}^{-3}$, and $B_0 = 7.60^{+2.00}_{-3.47}\,\mu\text{G}$.

We present the results from a 3-year-long Ultraviolet monitoring campaign of the broad line Seyfert 1 galaxy ESO 141-G55 using International Ultraviolet Explorer (IUE). By modelling all individual, extinction-corrected UV spectra in 1150-1978 A and 1850-3348 A wavelength range, we have observed a significant variability in both UV continuum and line fluxes. Variabilities due to ionised UV lines like SiIV, CIV and HeII are delayed with respect to the UV continuum by 2.92$^{+0.54}_{-0.61}$, 4.41$^{+0.44}_{-0.54}$, 4.11$^{+0.35}_{-0.81}$ days, respectively. At a distance of $\sim$0.004c, an outer accretion disc can be a possible site for the origin of UV lines.

Aims: We aim to characterise the multiscale statistical properties of the reconstructed velocity and density fields of the nearby universe, identify possible scaling regimes, quantify intermittency, and assess indications for the transition toward large-scale homogeneity within the range probed by current data. Methods: We analyse the CosmicFlows4 three-dimensional velocity and density-contrast cubes using absolute structure functions of arbitrary order, $q$. The analysis is performed within a volume extending to $z \lesssim 0.08$ ($\simeq 350~\mathrm{Mpc}$ $h^{-1}$). Structure function scaling exponents $\zeta(q)$ are estimated from configuration-space statistics. Intermittency is characterised using the Universal Multifractal formalism, and probability density functions of increments are examined. Results: Two regimes are detected. Small separations are dominated by reconstruction smoothing and show nearly linear $\zeta(q)$ behaviour. At larger separations, a scaling regime appears with $\zeta_\rho(1)\simeq0.3$ ($D_\rho\approx3.7$) and $\zeta_v(1)\simeq0.4$. The correlation function follows $\xi(r)\sim r^{-1.4}$ over $[45,250]~\mathrm{Mpc}\,h^{-1}$, implying $D_2\simeq1.6$. Non-linear $\zeta(q)$ and Lévy-stable increment PDFs indicate intermittency and strong non-Gaussianity. Velocity increments show a systematic negative skewness suggestive of a cascade-like asymmetry associated to amplification of negative, compressive gradients.

We present a spectroscopic analysis of XRISM and NuSTAR observations of the black hole X-ray binary GS~1354--64 during its 2026 outburst. A total number of 3.5 million photons are collected by the microcalorimeter Resolve on board XRISM, providing an unprecedented high-resolution view of the iron line profile. A clean broad iron line is found in the data, without significant narrow features. Modeling the broad iron line with relativistic reflection from the inner accretion disk suggests a rapidly spinning black hole (a>0.98) in the system. Measurements of the disk inclination angle from the reflection method are model-dependent. This work demonstrates the power of X-ray microcalorimeters in studying the inner accretion flow and constraining black hole parameters.

QSOs are essential for investigating the structure and evolution of the Universe. Historically, their identification has been concentrated in the northern hemisphere, primarily due to the sky coverage of major astronomical surveys. The QUBRICS survey, started in 2019 to address this asymmetry, has identified more than 1300 new bright (i<19.5) high-redshift (2.5<z<6) QSOs in the southern sky. We aim to quantify, using an independent QSO sample, the completeness and recall of the QUBRICS QSO selection methods, based on XGB (eXtreme Gradient Boosting) and PRF (Probabilistic Random Forest), since completeness is a fundamental metric for ensuring the statistical robustness of QSO-based cosmological investigations. A subset of Gaia DR3 sources with low-resolution spectra was analyzed, obtaining a sample of 3501 QSOs. To determine how many QSOs were correctly identified as candidates, we crossmatched this independent sample with the datasets used for selection: 894 QSOs with z>2.5 fell within the XGB dataset footprint, of which 152 were unclassified and thus eligible for completeness testing. Similarly, 675 QSOs with z>2.5 were within the PRF dataset footprint, including 69 unclassified objects. The XGB correctly identified as candidates 136 (89%) of the 152 QSOs with z>2.5 present in its dataset as unclassified objects. The PRF correctly identified as candidates 46 (66%) of the 69 QSOs with z>2.5 present in its dataset as unclassified objects. These findings confirm the high efficiency of the QUBRICS selection methods (recall=89%) and provide the completeness estimate for spectroscopically confirmed QSOs (82%), necessary for cosmological studies using QUBRICS data. This work also provides reliable redshifts for 1223 new QSOs (median redshift z=2.1 and magnitude G=17.8), that will help improve the performance of future selections.

Interstellar objects provide the only directly observable samples of icy planetesimals formed around other stars, and can therefore provide insight into the diversity of physical and chemical conditions occurring during exoplanet formation. Here we report isotopic measurements of the interstellar comet 3I/ATLAS, which reveal an elemental composition unlike any Solar System body. The water in 3I/ATLAS is enriched in deuterium, at a level of D/H = (0.95 +- 0.06)%, which is more than an order of magnitude higher than in known comets, while its range of 12C/13C ratios (141-191 for CO2 and 123-172 for CO) exceeds typical values found in the Solar System, as well as nearby interstellar clouds and protoplanetary disks. Such extreme isotopic signatures indicate formation at temperatures $\lesssim30$ K in a relatively metal-poor environment, early in the history of the Galaxy. When interpreted with respect to models for Galactic chemical evolution, the carbon isotopic composition implies that 3I/ATLAS accreted roughly 10-12 billion years ago, following an early period of intense star formation. 3I/ATLAS thus represents a preserved fragment of an ancient planetary system, and provides direct evidence for active ice chemistry and volatile-rich planetesimal formation in the young Milky Way.

Halos surrounding spiral galaxies act as the bridges connecting the galactic disk and the intergalactic medium (IGM). They host a significant fraction of the baryonic mass in the Universe, and feedback from star formation (SF) or active galactic nuclei (AGN) likely plays an important role in regulating this vertical baryonic component. Despite its importance, the contribution of extraplanar molecular gas remains poorly understood. We aim to characterize the vertical extent and the kinematics of molecular gas traced by CO(2-1) emission in the nearby (D = 9.5 Mpc) spiral galaxy NGC 891, one of the best studied edge-on galaxies. We also compare our results with HI, H$\alpha$-traced DIG and dust maps from the literature. Our analysis is based on new CO(2-1) observations of NGC 891 obtained with the IRAM 30m telescope. We mapped two 6 kpc $\times$ 6 kpc regions on the northeastern side and the area surrounding the galactic center. We apply a careful method to estimate and remove the residual contribution of the error beam to the CO cube. The vertical extent of the molecular gas is best described by a two-component Gaussian fit, consisting of a bright thin disk component (deconvolved FWHM $\simeq$ 360 pc) and a fainter thick disk component (deconvolved FWHM $\simeq$ 1.1 kpc). Statistically significant CO(2-1) emission is detected up to 1.3-1.4 kpc above the disk midplane. We estimate that the thick molecular disk component contains up to 27% of the total molecular gas mass of the galaxy. Our results demonstrate that SF-driven feedback in a non-starburst galaxy can lift significant amounts of molecular gas to large vertical distances. We interpret the presence of extraplanar molecular gas in NGC 891 in the framework of a galactic fountain scenario, in which material is expelled from star-forming regions and transported toward the outer halo.

3C 286, a well-known radio calibrator, exhibits the stability in both of total flux density (FD) and polarization parameters. However, its stable and luminous interstellar radio signal may encounter interplanetary scintillation (IPS) due to density irregularities in the solar wind within the heliosphere. In this work, we analyze high-time-resolution observations of 3C 286 obtained with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) from 2019 to 2023. Our analysis reveals that IPS affects the polarized flux densities of the Stokes I, Q, and U parameters, whereas Stokes V shows no detectable IPS-induced variations. The IPS variations detected in Stokes I are synchronous with those in Stokes U, while those in Stokes Q exhibit greater randomness. The crosscorrelation function (CCF) results indicate no time delay between Stokes I and U, but a delay of approximately 2.8 seconds between Stokes I and Q. This suggests that the different polarized radio emissions of 3C 286 originate from distinct emission regions, specifically the core and the southwestern jet. Furthermore, the projections of the radio core and jet component onto the scintillation screen at 1 AU yield a solar wind plasma speed of $\sim 637$ km/s.

Different models of filament formation predict distinct patterns of angular momentum redistribution toward embedded cores, set by the underlying velocity-field structure, which can set the initial conditions for a preferential orientation between protostellar outflows and filaments. However, the absence of a dominant alignment in observations keeps this connection open to debate. We investigate whether gravity-driven longitudinal flows along filaments can redistribute angular momentum (AM) toward collapse centers and influence outflow-filament alignment. To this end, we analyze the distributions of 3D and 2D-projected angles between sink angular momentum vectors and host filament orientations in an SPH simulation of giant molecular cloud and filament formation. We also characterize the filament velocity field by measuring the angles between SPH particle velocity vectors and filament axes, and the degree of convergent flow toward filament density peaks. No preferred alignment between the sinks' AM and the filament direction is found at early evolutionary stages, neither in 3D nor in 2D. Later, however, a predominantly perpendicular configuration emerges in 3D. Tracking individual sinks indicates that this alignment is not primordial but develops as gravity strengthens. In individual filaments, the onset of perpendicular alignment coincides with the development of convergent longitudinal flows. Finally, we estimate the minimum fraction of perpendicular 3D angles required to reveal a perpendicular 2D alignment for a given sample size. While longitudinal flows develop over extended timescales, once established, they can rapidly reorient the angular momentum vector of the sinks, enabling perpendicular alignments to arise within typical outflow lifetimes.

Fast radio bursts (FRBs) provide a sensitive probe of ionized baryons through their dispersion measure (DM). In addition to slowly evolving cosmological terms, at least two repeaters now show clear secular DM-decrease episodes: FRB~20190520B and FRB~20121102 , supporting a dense, dynamically evolving local environment. We adopt a \emph{forward-modeling} approach and use time-dependent 1D SNR simulations for a young magnetar embedded in SN ejecta, combining single-star and binary-stripped progenitors with HD+NEI calculations to follow shock structure, ionization, and electron density. The shocked region contributes only limited DM ($\lesssim10\,{\rm pc\,cm^{-3}}$), while the dominant time-varying component is the unshocked ejecta, whose early behavior follows ${\rm DM}\propto t^{-\alpha}$ with $\alpha\simeq1.8$--$1.9$. Although shocked-region DM is small, shock-amplified magnetic fields can still generate substantial RM; in our shock-only RM framework, only the $11\,M_\odot$ SS model reproduces the FRB~20121102 RM evolution. Binary-stripped progenitors generally yield smaller DM than single-star models at fixed $M_{\rm ZAMS}$, with composition-dependent mean molecular weights introducing non-monotonic mass trends. Matching the observed ${\rm dDM}/{\rm d}t$ of FRB~20190520B (and the late-stage slope of FRB~20121102), we infer local SNR DM contributions of tens to hundreds ${\rm pc\,cm^{-3}}$. We also find GHz escape is allowed in most models, with $\tau_{\rm ff}=1$ typically reached by $t_{\rm esc}\lesssim70$ yr; for weakly ionized ejecta, the source can be nearly transparent from very early times. These results support a young CCSN/SNR origin for a substantial fraction of ${\rm DM}_{\rm source}$ and highlight that physically consistent local-environment modeling is essential for robust FRB cosmological DM inferences.

All detected water reservoirs in the solar system exhibit a deuterium enrichment that links back to the physical environment at the time of stellar birth. Gas-phase and ice-grain D-enrichments occur through chemical processes that operate at low temperatures ($<$~30~K) pointing towards an origin in the pre-stellar molecular cloud or in the outer parts of the protoplanetary disk. However, not all stars are born in environments similar to our Sun, nor do their subsequent evolutionary histories follow the same path. These environmental differences can be traced by the water D/H ratio. Here we use ALMA observations of the interstellar comet 3I/ATLAS to constrain the water D/H ratio in extrasolar cometary material. With a water D/H value of [D/H]$_{\mathrm{H_2O}} > 6.6\times10^{-3}$, 3I/ATLAS shows a deuterium enrichment exceeding Earth's ocean value by more than a factor of $\gtrsim40$ and typical Solar System cometary values by more than a factor of $\gtrsim30$. The elevated deuterium enrichment points to water that formed under colder, less irradiated conditions and from less thermally processed material, consistent with an origin in a planetary system that formed under different physical and chemical conditions than our own.

We report the first NuSTAR hard X-ray observations of the low-luminosity active galactic nucleus NGC 4278. The source is clearly detected beyond 10 keV with a hard X-ray spectrum consistent with a power law of photon index between 2.2 and 2.5 without evidence for a high-energy cutoff. The X-ray flux is low compared to the active state in 2021, but exhibits variability by a factor of ~2 on a timescale of a month. We discuss the origin of the hard X-ray emission and explore its connection to gamma rays and high-energy neutrinos. We explain the X-ray data, including both quiescent and active states, using a radiatively inefficient accretion flow (RIAF) model with a variable accretion rate. We also show that TeV gamma rays cannot escape from the RIAF disk, and very high-energy gamma rays observed in LHAASO are likely to originate from outer regions such as jets and winds, which is consistent with our results favoring a magnetically arrested disk. We also discuss hidden neutrino emission from RIAFs together with possible connections to coronae of active galactic nuclei with standard, radiatively efficient disks.

We present a new calibration of the broad-line region (BLR) size-luminosity-velocity relation using a sample of 157 AGNs with reliable Hbeta time-delay (\lag) measurements from Wang & Woo 2024. By incorporating the Eddington ratio as a third parameter, we effectively correct the systematic offset of high-Eddington AGNs in the traditional BLR size-luminosity relation. The resulting three-parameter fit defines a fundamental plane in the 3-D space of the \lag, optical luminosity, and Hbeta velocity, with an intrinsic scatter of 0.21 dex. This tight correlation reflects the coupled effects of gas kinematics, photoionization, and BLR geometry. In turn, we develop a new method to infer \lag\ from the combination of optical luminosity and Hbeta velocity, and derive single-epoch black hole mass estimators by adopting either the full-width-at-half-maximum (FWHM) or line dispersion ($\sigma$) of the Hbeta line profile as the velocity indicator. The derived \lag shows a ~0.1 dex scatter, depending on the choice of calibrations. We show that the previous mass estimates based on the two-parameter size-luminosity relation with a 0.5 slope can be overestimated by up to 0.5 dex, demonstrating that the new mass estimator substantially changes the cosmic black hole mass density and the growth of black hole seeds in the early universe.

Pulsar halos are produced by electrons and positrons diffusing in the interstellar medium around their parent pulsar wind nebulae. Recent observations by HAWC and LHAASO have revealed asymmetric morphologies in the halos surrounding Geminga and Monogem. The anisotropic diffusion model provides a natural explanation for such asymmetries, where the morphology is determined by the viewing angle of the mean magnetic field, the Alfvénic Mach number ($M_{\rm A}$), and the pulsar distance. In this work, we model the measured morphologies based on this framework and constrain the properties of interstellar magnetic turbulence. We find that the mean magnetic field orientations within the two halos are different, implying that they reside in different magnetic coherence regions, whereas the Alfvénic Mach numbers are relatively close ($M_{\rm A}\sim 0.2$). The results suggest a local magnetic field coherence length of approximately 100pc. Our study demonstrates that the morphology of pulsar halos serves as a powerful diagnostic tool for the properties of interstellar magnetic fields, highlighting the need for more accurate morphological measurements and sophisticated diffusion modeling in future studies.

We present that supernovae interacting with a dense Si and S-rich circumstellar matter like SN 2021yfj can originate from mergers of two white dwarfs. A C+O white dwarf accreting He from its non-degenerate He companion star can initiate a C burning frame at its surface propagating inward under certain conditions. Such a burning frame synthesizes intermediate mass elements such as Si and S, forming a hybrid WD with an outer Si+S-rich layer. After the He star companion becomes a white dwarf, the two white dwarfs can eventually merge. During the merger, the outer layers of the hybrid white dwarf can be tidally stripped, forming a dense Si and S-rich circumstellar matter. If a thermonuclear explosion is triggered in the merging white dwarfs, an explosion within a dense Si and S-rich circumstellar matter can be realized, resulting in SN 2021yfj-like events. We argue that the properties of SN 2021yfj can be reproduced by a dense Si and S-rich circumstellar matter having ~ 0.3 Msun within which an explosion having kinetic energy of ~ 4e50 erg and ejecta mass of ~ 0.3 Msun occurred. These properties are consistent with the double white dwarf merger scenario. This scenario can naturally explain the existence of He observed in SN 2021yfj. Because white dwarf mergers can also lead to the formation of He and C+O dense circumstellar matter, some Type Ibn and Icn supernovae may also originate from a similar evolutionary path.

The existence of dark matter is supported by multiple astrophysical observations, yet its particle nature remains unknown. The development of gravitational wave astronomy, especially with future space-based detectors such as LISA, provides new opportunities to study the interactions between dark matter and compact-object systems. This review summarizes the main dark matter candidates and their macroscopic distributions, and highlights three mechanisms through which dark matter can affect gravitational wave observations: (1) modifications to compact-object orbits and the dynamics of systems such as extreme mass-ratio inspirals, including dark matter spikes, dynamical friction, and potential perturbations; (2) gravitational lensing effects induced by the spatial distribution of dark matter, altering waveform amplitudes and phases; and (3) direct couplings between ultralight dark matter fields and detectors. As low-frequency gravitational wave detection techniques are proposed and continue to develop, these effects may offer a novel avenue for probing the properties of dark matter, and combining precise waveform modeling with multi-messenger observations could reveal insights into its microscopic structure.

A discrete tomography method has been developed that is able to reconstruct three-dimensional coronal mass ejection (CME) density structure. We test the method by producing synthetic coronagraph imagery for three CME events using the CORona--HELiosphere (CORHEL) model. We combine images from different numbers of observing spacecraft, ranging from three to seven, and we perform the method separately using polarimetric and non-polarimetric reconstructions, as a means to test their relative effectiveness. In all cases, we show that increasing the number of observing spacecraft reduces the mean relative absolute error (MRAE) between the simulated and reconstructed density. Furthermore, the MRAE is generally lower when using polarimetric reconstructions compared to non-polarimetric reconstructions. Methods applied to localise the CME front work well for all spacecraft configurations, and are improved when using polarimetric, over non-polarimetric, reconstructions. The presence of a CME front in the simulations can be identified with an accuracy of $(72\pm9)\%$, $(70\pm8)\%$ and $(52\pm12)\%$ for CME1, CME2 and CME3 via polarimetric reconstructions using only three spacecraft at L1, L4 and L5. The radial position of the CME front can be constrained to a high level of precision when using polarimetric reconstructions using the same three spacecraft; $0.003\pm0.004$\,au, $0.004\pm0.005$\,au and $0.005\pm0.004$\,au for CME1, CME2 and CME3, respectively. We expect that at least four spacecraft are required in order to derive accurate information about 3D CME structure. We find no strong evidence of improvement when including out-of-ecliptic observers, but that their inclusion increases the volume of space within which the inversion can be performed.

We investigate how diversity in the core-halo mass relation affects constraints on the fuzzy dark matter particle mass $m_\psi$ inferred from the internal kinematics of dwarf galaxies. Using stellar line-of-sight velocities and projected positions for eight Milky Way dwarf spheroidal galaxies, we model their dark matter halos as solitonic cores embedded within outer NFW envelopes. We apply both second- and fourth-order Jeans analyses to derive the posterior distribution of $m_\psi$. Our results show that there are two ranges of $m_\psi$ consistent with the observed kinematics: $-20.3 < \log_{10}(m_\psi/\mathrm{eV}) < -19.2$, and a narrower small-mass window $-22.1 < \log_{10}(m_\psi/\mathrm{eV}) < -21.5$, both within the 68\% credible intervals. The latter becomes prominent only if core-halo diversity is taken into account. These constraints pose a challenge to fuzzy dark matter, as the small-$m_\psi$ window is in conflict with Milky Way satellite abundances, and our upper bound largely excludes the parameter space permitted by Lyman-$\alpha$ forest constraints.

Interstellar objects provide a unique opportunity to further our understanding of the planetary formation process by studying in detail material formed around another star. Their ices contain precious clues about the environment and conditions prevailing in their home system. As fractionation processes can be sensitive to the temperature and radiation environment, isotopic ratios are powerful tracers of the origin and evolution of different species. While isotopic ratios have been measured in solar system comets, previously detected interstellar objects have been too faint to measure isotopic ratios. Here we report the measurement of two ratios in 3I/ATLAS from observations of the CN molecule: $^{12}$C/$^{13}$C and $^{14}$N/$^{15}$N. We report $^{12}$C/$^{13}$C=$147^{+87}_{-40}$ and $^{14}$N/$^{15}$N=$343^{+454}_{-124}$. The $^{14}$N/$^{15}$N is higher than the value of $\sim$~150 usually measured for solar system comets, close to the values measured in the interstellar medium, pre-stellar phases or the outside of protoplanetary discs. The $^{12}$C/$^{13}$C is marginally higher than the values usually measured for solar system comets and in the interstellar medium. These measurements could indicate an origin of 3I in the outer disc around an older low-metallicity star.

We investigate the magnetic origin of the coronal mass ejection (CME) that occurred on November 28, 2023, at 19:50UT from active region (AR) 13500 located near the solar disk-center. The eruption was associated with an S-shaped sigmoidal structure formed by the inner AR polarities along a sheared polarity inversion line, while the outer polarities evolved through proper motions. During November 26-28, the AR exhibited a decrease in net magnetic flux while progressively injecting magnetic helicity and energy into the corona toward the eruption onset, highlighting the key role of helicity-injection in triggering eruptions. To simulate this magnetic evolution, we employed a data-driven magnetofrictional (MF) simulation starting 2.8 days prior to the eruption. The energy input for the model was constrained using the observed energy injection through an ad-hoc parameter. The initial potential-field configuration gradually evolved into a sheared-arcade and eventually developed into a twisted flux rope (FR) over the observed time-scale. Proxy emission maps based on electric currents show remarkable morphological agreement between the simulated and observed sigmoidal structure. The average FR-core twist increasingly builds-up leading the FR to initiate slow-rise motion of FR top from 50Mm until its eruption onset at 80Mm. Importantly, the ratio of current-carrying to total relative-helicity increased from 0.13 at FR formation to 0.30 at eruption, when the FR core entered the torus-unstable regime, suggesting an association between torus-instability and a threshold helicity ratio. These results demonstrate that data-driven MF simulations can successfully reproduce the evolving coronal magnetic configuration and may provide a robust tool for assessing the eruptive potential of ARs, particularly the helicity ratio.

The environments around star clusters evolve as stellar feedback reshapes the interstellar medium and dynamical processes reorganize the structure of the surrounding stellar field. As approximately single-age populations, star clusters can serve as clocks to trace these environmental changes. In this exploratory study, we test whether convolutional neural networks (CNNs) can identify age-dependent changes in cluster environments. We take cluster ages as given from basic SED fitting of five-band UV-optical aperture photometry from the PHANGS (Physics at High Angular resolution in Nearby GalaxieS) HST survey. We first show that CNNs can be trained on image cutouts centered on clusters to recover ages directly from imaging. This demonstration provides the foundation for this study, which examines whether the information used by CNNs to predict age is coherent and physically meaningful. We perform controlled image occlusion experiments as an explainable AI method. These show that the CNNs extract age-predictive environmental cues in the absence of cluster light and when information on SED shape is removed by combining the five filters into one image. We find that reliance on environmental information increases at the youngest (<10 Myr) and oldest (>1 Gyr) ages, where clusters can exhibit similarly red colors. Our results are consistent with the long-recognized picture that cluster environments evolve systematically with age. We demonstrate that this information is encoded at a level detectable by machine-learning and recoverable from broadband imaging. This establishes a path for using new techniques to connect image-based age inference to the physical evolution of cluster environments.

We show that the flip-flop transitions in X-ray binaries (rapid cycling between different spectral states which are sometimes seen near the global state transition) show a series of analogies to the changing state phenomena (rapid changes in the emission line properties that seem to be driven by changes in the central engine) in active galactic nuclei (AGN). Specifically, (1) the timescales for the transitions scale approximately linearly with mass and (2) both phenomena occur at a few percent of the Eddington luminosity. Because most accretion physics is expected to be scale-free, it is likely that these represent two manifestations of the same phenomena. Demonstrating this would allow the use of a much wider range of observational techniques, on a much wider range of characteristic timescales, and provide a clearer pathway toward understanding these rapid transitions than is currently available. We discuss potential means to establish the connection more firmly, and to use the combination of the observational advantages of both classes of systems to develop a better understanding of the phenomenon.

One of the most interesting questions that astronomy can hope to answer is: are we alone in our Milky Way galaxy? A detection of an electromagnetic (EM) signal generated by an extraterrestrial technological intelligence (ETI), or the presence in our solar system of an alien probe, would answer this question in the affirmative. Purposeful interstellar communication is a 2-way street - the transmitting and receiving technological intelligence (TI) both need to do its part. As the receiving TI, our EM search programs should incorporate a model of what a transmitting TI is likely to be doing. Published searches for extraterrestrial technological intelligence (SETI) have generally not done so and, thus, have often been sub-optimally designed. We propose an improved search technique that more closely corresponds to astronomical surveys that have been undertaken for reasons that have nothing to do with SETI. Published non-SETI radio and optical surveys are sufficiently extensive that they already supply meaningful constraints on the prevalence of nearby purposely communicative alien civilizations. Purposeful communication can also include the sending of spaceships (probes). The absence of evidence for alien probes in the solar system suggests that no alien civilization has passed within 100 light-years of Earth during the past few billion years.

Black hole-neutron star mergers, together with binary neutron star mergers, are key laboratories for neutron star physics. They enable us to probe merger dynamics imprinted in gravitational waves and potential electromagnetic counterparts. These systems link microphysics and macrophysics by placing constraints on the dense matter equations of state, potentially revealing the imprint of hadron-quark phase transitions, clarifying the role of neutrino irradiation in shaping the ejecta, its r-process nucleosynthesis, and kilonova emission, as well as assessing how magnetically driven instabilities affect mass ejection and possible electromagnetic signatures. Despite their importance, black hole-neutron star mergers remain relatively less studied and therefore not yet well understood, largely due to the lack of publicly available numerical relativity setups suitable for such investigations. In this work, we present a fully reproducible black hole-neutron star merger simulation performed exclusively using Einstein Toolkit thorns, targeting the detected event \texttt{GW230529}. The simulations are carried out at three resolutions with finest grid spacings of $162$, $222$ and $310$ meters to assess numerical robustness. The entire setup, from initial data to a parameter file with some of the analysis scripts, is publicly released as a new Einstein Toolkit gallery example and will be distributed as part of the Hypatia release, establishing a reference black hole-neutron star merger configuration within the Einstein Toolkit.

Catastrophic planetesimal disruptions offer a unique opportunity to study and characterize large planetesimal populations in exoplanetary systems that are not currently detectable by modern observatories. The unexpected discovery of a second collision event in the Fomalhaut system raises important questions about the planetesimal population and dynamical state inside the Fomalhaut main belt that led to two collisions in 20 years. We present a statistical model developed and applied to the archetypal Fomalhaut system to provide new constraints on the bulk properties of the planetesimals in Fomalhaut's main belt. Utilizing the constraints provided by the spatially resolved Fomalhaut cs1 and cs2 collision events, we retrieve the belt parameters that best reproduce the observed collision rate while remaining consistent with the system's age and dust mass. Our best-fit model suggests a total main belt mass of 200-360 $M_{\oplus}$, with the transition from a collisionally evolved to a primordial planetesimal population occurring at a radius of $115_{-10}^{+30}$ km and a maximum planetesimal radius of $380_{-202}^{+643}$ km. We estimate a catastrophic collision rate of $0.086_{-0.048}^{+0.067}$ collision events per year for planetesimals with radii $\ge$ 100 km in the region interior to the main belt. Our findings show that further observable collisions are likely, motivating continued monitoring of Fomalhaut and other nearby debris disks.

We present an analysis that investigates the $H_0$ and $S_8$ tensions by considering a dark energy model. The latter is a late-time model characterized by a future abrupt event known as the Little Rip (LR) model and characterised by one extra parameter, $\beta$, compared to the standard model, $\Lambda$CDM. To test this approach, we perform a statistical analysis by the MCMC method using the most recent observational data. We obtain a positive correlation in ($H_0$, $\beta$) plane. We also note that the Hubble tension is less than $3\sigma$ when using early measurements, i.e., Cosmic Microwave Background (CMB) data, and when combining it with Baryon Acoustic Oscillation (BAO) data, but it is no longer so when we combine early and late measurements (i.e. PantheonPlus (PP)). In addition, we test the model with DESI-DR2 combined with CMB and recent SNIa measurements. We notice that our model shifts toward the quintessence field. For a complete statistical analysis, we use the Akaike Information Criteria and Bayesian analysis of the evidence. According to Bayes factors, we find that the LR model provides an improved fit only to CMB data.

Gamma-Ray Bursts (GRBs) are the most energetic electromagnetic explosions in the universe. Recently, the Large High Altitude Air Shower Observatory (LHAASO) reported the breakthrough observation of GRB 221009A with gamma-ray energies beyond 13 TeV. This discovery, together with the previous GRB detection well above 100 GeV, confirms the production of very-high-energy (VHE, $\gtrsim 100$ GeV) radiation which might be a common component of all bright GRBs. It is reasonable to expect that bright GRBs are important targets for ground-based gamma-ray experiments. In this work, we estimate the detection rate for current and upcoming ground-based gamma-ray observatories including LHAASO, Large Array of Imaging Atmospheric Cherenkov Telescopes (LACT) and the Southern Wide-field Gamma-ray Observatory (SWGO) under two emission models with GRB~221009A as the template: first, that they all share the same intrinsic VHE spectral shape; second, they have the same environmental parameter and electron spectral index, governing their synchrotron self-Compton (SSC) emission. Using the long GRB luminosity and redshift distribution function obtained from the Fermi-GBM GRB samples, and accounting for the cosmological effects and extra-galactic background light (EBL) absorption, we derive the expected VHE flux at Earth. The sensitivity analysis for LHAASO, the upcoming LACT, and SWGO to evaluate their detection potential across specific redshift and luminosity ranges has been performed. The corresponding 5$\sigma$ detection rates of 221009A-like GRBs for the two emission models are: LHAASO, 0.04-0.05 yr$^{-1}$; LACT, 0.03-0.06 yr$^{-1}$; SWGO, 0.2-0.4 yr$^{-1}$. These rates can vary by up to $\approx 24\%$ due to different EBL models.

As the oldest building blocks of our Galaxy, globular clusters retain the archaeological footprint of the early stellar environments. Accurate absolute ages of globular clusters are required to interpret this ancient record. Existing dating techniques often produce precise but discordant ages, suggestive of systematic errors in excess of 1 Gyr. The James Webb Space Telescope (JWST) has unlocked a new dating method that leverages the cooling behavior of previously unobservable brown dwarf members. With a largely independent set of systematic errors, this new method provides a new consistency test for more established methodologies. I present a likelihood-based histogram-free method to derive globular cluster ages from multi-band JWST photometry of cluster members near and below the hydrogen-burning limit. By applying the method to a large set of simulated observations, I establish that formal age errors (i.e. errors based on measurement uncertainties alone) under 0.2 Gyr are attainable for nearby globular clusters. I also evaluate the significance of associated systematic effects, including the chemical heterogeneity of globular clusters (multiple populations), unresolved binary systems and uncertainties in brown dwarf cooling rates. As with other methods of age determination, systematic effects dominate the error budget (in selected cases, by over an order of magnitude), but may be reduced with more sophisticated analysis. Finally, I provide a lookup table for determining the number of observations, exposure times and temporal baselines required to estimate the age of a given cluster to a prescribed precision.

We propose a unified model for jet formation applicable to active galactic nuclei, young stellar objects, and X-ray binaries. In this model, the binding energy released from the accretion disk is primarily stored as turbulence rather than being radiated away, leading to the formation of advection-dominated accretion flows. Near the central object, a thick accretion disk with funnel-like structures develops. Within the turbulent flows, the smallest stable blobs can be accelerated beyond the escape velocity through a mechanism involving the combined effects of inward pressure forces and angular momentum conservation. These rapidly moving blobs may exit through the funnels, collectively forming two opposing jets. This model predicts that jets originate from the innermost region of the thick disk surrounding the central object. It can be extended to account for jet formation in active galactic nuclei, young stellar objects, X-ray binaries, and other analogous astronomical systems.

The rapid expansion of low-Earth-orbit (LEO) megaconstellations introduces new risks to radio astronomy from unintended electromagnetic radiation (UEMR). In this work, we present an attempt to search for UEMR from Starlink satellites using the 21 Centimeter Array (21CMA). Because the sensitivity of a single pod observation is limited, we focus on developing a robust observing and detection pipeline. Using Two-Line Element (TLE) data, we predict satellite transit times to guide the observations, and we define entry into the field of view (FoV) as an apparent declination greater than $85^{\circ}$ with respect to the 21CMA. We analyze the system equivalent flux density (SEFD) and the resulting single-pod sensitivity limits, which explain the detection of emission originating from the ORBCOMM satellites, rather than any detectable broadband UEMR in our dynamic spectra. To validate the methodology, we developed a Python package, orbdemod, to demodulate ORBCOMM downlink signals in our data. The recovered satellite ID agrees with the satellite predicted by our maximum-declination analysis, thereby validating the accuracy of our transit prediction and identification framework. Furthermore, via modulation power spectrum analysis, we show that the impulsive broadband bursts are produced by power line arcing near the array rather than by satellite UEMR.

A growing number of white dwarfs (WDs) exhibit one or more signs of remnant planetary systems, including transits, infrared excesses, and atmospheric metal pollution. WD 1054-226 stands out for its unique, highly structured, and persistent photometric variability. We aim to investigate the long-term stability and nature of the periodic signals observed in WD 1054-226 to better understand the origin and evolution of its transiting material. We analyse all available TESS light curves from Sectors 9, 36, 63, and 90 using Lomb-Scargle (LS), Box-Least-Squares (BLS), and Gaussian process (GP) periodogram analyses. We complement these with multiband, high-cadence ground-based photometry from LCOGT, MuSCAT2, ALFOSC, and ProEM to test for colour dependence and confirm the periodicities. We confirm the persistence of the previously-reported 25.01 h and 23.1 min periodicities over a six-year baseline. The 25.01 h signal shows some temporal evolution, while the 23.1 min dips are highly coherent on long timescales. A transient 11.4 h feature, previously reported, is detected only in early TESS sectors and is absent in recent data. No significant colour dependence is found in the ground-based observations. The stability of both the 25.01 h and 23.1 min signals indicates a long-lived, dynamically sculpted debris structure around WD 1054-226. The lack of colour dependence implies high optical depth, consistent with an opaque, edge-on debris ring rather than an optically thin dust population. This makes WD 1054-226 a key laboratory for testing models of remnant planetary systems around white dwarfs.

Realistic modeling of stellar spectra requires accurate radiative opacity coefficients. Owing to the fragmentary nature of existing data from rigorous quantum-mechanical calculations, photoionization coefficients based on the rigid-wavefunction approximation remain the only practical option for studies of magnetic white dwarfs. Although variants of this approach have been widely used in spectral analyses for decades, a complete and explicit treatment of degeneracy-level breaking has not previously been presented. In this work, we provide a comprehensive description of this procedure, including explicit expressions for the photoionization probability of individual bound-free transitions as functions of magnetic field strength and radiation polarization. We also evaluate the occupation numbers of bound states in a magnetized gas under ionization equilibrium, enabling the calculation of absolute photoionization opacities. Because high-lying atomic states are strongly perturbed by the magnetic field and ultimately dissolved, substantial modifications of the monochromatic absorption are found even for fields below 10 MG--a regime where fully rigorous quantum calculations are numerically demanding and have not yet been applied. Over a wide range of magnetic field strengths, pronounced dichroic features appear in the hydrogen continuum absorption.

Theoretical formation models and exoplanet detection surveys indicate that systems with multiple giant planets are common. We investigate how multiple super-thermal mass planets embedded in a circumstellar disk shape the dust distribution and examine the consequences for interpreting disk substructures and inferring planetary properties. We perform two-dimensional hydrodynamical simulations with a modified PLUTO code, treating dust as Lagrangian particles in a wide range of sizes. We analyze systems with two planets of different masses and orbital separations, comparing them to the single-planet scenario. We generate synthetic ALMA continuum maps using RADMC-3D and compute the relative impact velocities of dust particles to assess potential limitations to grain growth. Dust morphologies in multi-planet systems cannot be described as a simple superposition of single-planet gaps. Secular planetary perturbations can generate multiple dust traps and asymmetric structures, while also exciting significant eccentricities in dust particle orbits. As a consequence, the locations and widths of dust rings and gaps depend on the size of the particles, the masses of the planet, and the orbital configurations. Synthetic continuum images may hide gaps carved by multiple planets, thereby complicating the interpretation of observed substructures. In addition, eccentricities induced in dust orbits lead to stronger gas drag, reducing the Stokes number for a given particle size, and the enhanced relative velocities associated with eccentric orbits can further suppress grain growth, promoting fragmentation and replenishment of small dust grains.

The chiral plasma instability (CPI) has been invoked as a possible mechanism for generating primordial magnetic fields in the universe and ultrastrong fields in neutron stars. We investigate chiral dynamos where the chirality imbalance is pumped by a source on a timescale $t_0$ and show that the CPI rate $\gamma$ is limited to $\gamma_0/(1+{\cal Q}^2)$, where ${\cal Q}= (\gamma_0 t_0)^{1/3}$ and $\gamma_0$ corresponds to models with instantaneously created chirality imbalance $(t_0=0)$. We then find that chiral flipping with rate $\Gamma_{\mathrm f}$ hinders the chiral dynamo if $\Gamma_{\mathrm f} >\gamma_0/(1+{\cal Q}^2)$ and completely suppresses it if $\Gamma_{\mathrm f} >\gamma_0/(1+{\cal Q}^{3/2})$. Realistic $t_0$ typically give ${\cal Q}\gg 1$, which makes the dynamo greatly vulnerable to the suppression by chiral flipping. The suppression is strong in protoneutron stars and may be (barely) avoided near the electroweak transition in the early universe.

We present multi-epoch optical spectroscopy of the interstellar comet 3I/ATLAS obtained from December 2025 to January 2026 (heliocentric distances 1.8 to 3.3 au). From these spectra, we derive post-perihelion production rates and/or mixing ratios for daughter molecules (CN, C$_3$, C$_2$, and CH) and gaseous metals (Fe I and Ni I). We also estimate a lower limit on the CO abundance based on the [O I] $\lambda6300$ line. The resulting outgassing profiles reveal a pronounced perihelion asymmetry, with volatile production declining more gradually outbound than during the inbound phase. In addition, the coma becomes less depleted in C$_2$ following perihelion, and shows a substantial increase, relative to H$_2$O, in the production rates of metals and potentially also CO. These trends may indicate the activation of subsurface material or compositional heterogeneity revealed by seasonal effects. The inferred high CO abundance further suggests possible decoupling of CO from H$_2$O and CO$_2$, similar to the behavior observed in solar system comets. The potential post-perihelion enhancement of both CO and metal production, if confirmed, would also be consistent with metal carbonyls contributing to the release of gaseous metals in cometary comae.

While space-based telescopes offer unparalleled precision for asteroseismology, ground-based observations remain crucial for identifying compact pulsator candidates and enabling their pulsational study through multi-colour photometry. The BlackGEM telescope array significantly enhances the detection and characterisation possibilities for compact pulsators that tend to be much fainter than dwarf or giant pulsators. Using BlackGEM multi-colour photometry of the hot pre-white dwarf PG 1159-035, we demonstrate its capability to detect short, multi-periodic pulsations with amplitudes down to a few milli-magnitudes. The primary aim of this study is to establish the feasibility of pulsation mode identification in hot subdwarfs and white dwarfs via mode amplitude-ratio analysis derived from BlackGEM multi-colour observations. Pulsation frequencies were extracted from our target using iterative pre-whitening analysis. To validate our data-driven mode identification concept using multi-colour photometry, we used the well-studied hot pre-white dwarf PG 1159-035, with previously identified pulsation modes, as a prototypical object that served for validation. The pre-whitening analysis using BlackGEM's standard q-, u-, and i-band light curves of PG 1159-035 revealed pulsation frequencies of l=1 and l=2 modes, consistent with values obtained from the literature. Subsequent amplitude ratio calculations resulted in discernible distributions for the l=1 and l=2 modes. The future assembly of more BlackGEM amplitude ratios for well-known white dwarfs with already identified modes will lead to density estimators suitable for identifying newly detected modes in known or as-yet-undiscovered pulsators. Our proof-of-concept study paves the way for large-scale asteroseismic analyses of optically faint compact pulsating stars using ground-based facilities, such as BlackGEM.

We present a sample of large double radio sources hosted by spiral galaxies (Spiral Double Radio Active Galactic Nuclei, SDRAGNs). Candidates were selected during Radio Galaxy Zoo, and refined using Sloan Digital Sky Survey images. The most promising were targeted in the Zoo Gems Hubble Space Telescope program, yielding images for 36 candidates. We assess the likelihood of each spiral galaxy being the genuine host of the radio emission finding 15 new high-probability SDRAGNs. SDRAGN hosts are seen preferentially close to edge-on. SDRAGNs predominantly show FR II radio structures and optical pseudobulges. Accounting for sample selection effects, the radio-jet axes lie preferentially near the poles of the galaxy disks; we find a constant probability distribution for intrinsic pole-jet angles < 30 degrees, ramping to zero at 60 degrees. We have obtained optical spectra for all these new SDRAGNs. Among both previous and new SDRAGN samples, 8/25 show Seyfert 2 signatures, 6/25 show central star formation, and 5/25 show LINER emission strong enough to indicate AGN or shock ionization, broadly similar to radio galaxies in elliptical hosts with the addition of star formation (diluting or masking weak AGN signatures). SDRAGNs include FR II sources seen at unusually low radio power, and preferentially occur in significant galaxy overdensities on 1-Mpc scales. Our "false alarms" - systems where HST data show the spiral to not be the actual host galaxy - include radio sources seen through large parts of foreground spiral disks, potentially useful for Faraday-rotation studies of disk magnetic fields.

We report the discovery of a new long-period radio transient, ASKAP J142431.2-612611, with a 36 minute period, identified in the Australian SKA Pathfinder Evolutionary Map of the Universe survey. We detected pulsed emission from ASKAP J142431.2-612611 over a period of eight days during follow-up observations with the Australia Telescope Compact Array, after which the source appears to have switched off. No optical or near-infrared counterpart is detected in archival surveys or in targeted Gemini South FLAMINGOS-2 observations. During its active state, the source exhibits a stable pulse profile with fractional polarisation consistent with 100%, evolving from elliptically to linearly polarised and tracing a well-defined great-circle trajectory on the Poincaré sphere. We show that this behaviour is consistent with fully linearly polarised intrinsic emission modified by propagation through a linearly polarised birefringent medium. This discovery expands the known population of long-period transients and highlights the intermittent nature of their activity. We discuss the implications for proposed models of long-period transients and outline future observations needed to constrain the origin of their intermittency and polarisation properties.

We present a particle physics motivated realization of interacting dark energy in which a radiatively stable dark energy sector couples to weakly interacting massive particle dark matter through pure momentum exchange. The dark energy field arises as a pseudo-Nambu-Goldstone Boson from a complex scalar singlet charged under a softly broken global $U(1)_S$, while dark matter is identified with an inert scalar doublet stabilized by a discrete $Z_4$ symmetry. This symmetry structure allows renormalizable dark matter-dark energy portal operators; however, requiring the dark energy field to emerge as a radiatively stable pseudo-Nambu-Goldstone Boson necessitates their absence, leaving derivative interactions as the leading coupling. As a result, energy transfer between the dark sectors is absent at the background level, while momentum exchange modifies the evolution of cosmological perturbations. We implement the resulting interacting dark energy model self-consistently in the Boltzmann code CLASS and study its impact on the growth of structure. We find that, despite sizeable momentum exchange, the suppression of the clustering amplitude $\sigma_8$ saturates above the level required to fully resolve current low-redshift tensions. Our results demonstrate that symmetry-protected, momentum-exchange-only dark sector interactions possess an intrinsic limit on structure suppression, providing a theoretically controlled benchmark for interacting dark energy scenarios.

Astrophysical neutrinos from Active Galactic Nuclei (AGN) offer a unique window into high-energy particle acceleration in obscured environments. The nearby Type II Seyfert galaxy NGC 1068 is a compelling example, exhibiting evidence of a high-energy neutrino excess without an associated TeV $\gamma$-ray counterpart. This suggests that hadronic processes may occur within an inner, magnetically dominated region, where the TeV emission is suppressed by $\gamma \gamma$ absorption and reprocessed via electromagnetic cascades in the dense, obscured environment. Building on our framework, which establishes turbulence-driven magnetic reconnection as the main driver for particle acceleration in this source, we present a refined lepto-hadronic model based on de Gouveia Dal Pino & Lazarian (2005) and Kadowaki et al. (2015). In these proceedings, we adopt a conservative inner disk radius compared to our previous results, moving the acceleration region further from the innermost stable circular orbit. We estimate the high-energy neutrino emission from hadronic and photo-hadronic processes, constrained by the acceleration timescale for first-order Fermi acceleration within the turbulent current sheet. The estimated model reproduces the IceCube neutrino flux excess, providing an essential technical complement and validation for our forthcoming comprehensive publication.

The Galactic Bulge, as the center of the Galaxy, is the closest laboratory for studying galaxy formation and evolution. However, its study faces significant challenges due to heavy dust extinction. This paper is devoted to deriving the average age of the Galactic Bulge and investigating its spatial distribution. We utilize a high-precision PSF-fitting photometric catalogue in the $J$ and $K_{\mathrm{s}}$ bands observed by VISTA to study the average stellar ages within the Bulge. Red giant stars are employed as tracers, with their average distances determined using red clump stars as references. The average ages are fitted with stellar models. Our analysis reveals a systematic age gradient across the Galactic Bulge ($2^{\circ} < |b| < 8^{\circ}$). The mean stellar age increases significantly with galactic latitude, shifting from a younger population ($\sim 4.69^{+0.97}_{-0.81}$ Gyr) prevalent near the plane to a predominantly older population ($\sim 10.48^{+0.93}_{-0.85}$ Gyr) at higher latitudes. We hypothesize that the young stellar population at low latitudes is predominantly composed of a pseudo-bulge formed via disk/bar processes (incorporating contributions from recent star-forming activity in the Galactic center), whereas the older stellar population is associated with spheroidal bulges generated through early-stage collapse or accretion of debris from merged dwarf galaxies.

The black hole X-ray binary (BHXRB) Swift J151857.0-572147 was discovered during its first outburst in March 2024. We review the archive of NICER observations from this outburst, focusing on the soft states. We select spectra for which the disk to total flux ratio exceeds 0.8 and the coronal scattering fraction fsc is less than 25%, conditions under which the accretion disk is expected to extend to the innermost stable circular orbit (ISCO) and remains geometrically thin. Through a continuum fitting analysis, we explore the dependence of spin on mass, inclination, and distance. We constrain the spin within the parameter space typical of stellar-mass black holes (sBHs), assuming a mass around 10 solar masses, inclination angles between 20 and 60 degree, and distances between 4 and 16 kpc. For fiducial parameters: a mass of 10 solar masses, a distance of 10 kpc, and an inclination angle of 40 degree, a moderate spin of approximately 0.7 is obtained. However, precise determination of the spin will require accurate measurements of these parameters. Our analysis provides a framework to infer the spin and estimate its uncertainties once more precise measurements of mass, distance, and inclination become available. As we demonstrate, lower inclination angles, greater distances, or larger black hole masses result in higher spin values.

Bipolar Light Bridges (BLBs) are bright regions located between sunspot umbrae of opposite magnetic polarity. They are typically characterized by strong magnetic fields and intense flows, which are believed to be closely associated with major solar flares. Despite their importance, their fine structure, formation and evolution remain poorly understood. In this work, we analyze the observations of a well-defined BLB obtained by the Goode Solar Telescope at the Big Bear Solar Observatory and the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory. The high-resolution GST observations reveal that the BLB is composed of fine, penumbral filament-like structures with widths of approximately 100-150 km. The corresponding Doppler velocity maps present a stable pattern of spatially adjacent red- and blueshifted patches within the BLB throughout the 5.5-hour GST observation. HMI observations show that the BLB arises from the converging and shearing motions of sunspots with opposite polarities. Penumbral regions originating from different polarities gradually evolve and interact, eventually forming the BLB. The observed Doppler velocity pattern, characterized by red- and blueshifted patches, can be interpreted as a projection effect of the Evershed flow within the penumbrae. Therefore, we argue that the BLB is formed through the compression and stretching of penumbral structures from oppositely polarized sunspots.

The spaceborne laser emission telescope is a core and critical component of the space gravitational wave detection this http URL with ground-based telescopes, the on-orbit space environment is more complex and harsh, presenting higher technical challenges for the design of the optical system and structure - both optical design and structural design face considerable difficulties. To meet the requirements of space gravitational wave detection, this paper designs a laser emission telescope based on an off-axis four-mirror configuration, with a capture field of view of 300{\mu}rad, an optical transmission efficiency of 86.3%, and an optical path stability index of TTL<0.025 nm/{\mu}rad. During the design process, based on existing theories and engineering experience, the primary mirror thickness optimization and lightweight structural design were completed, and a flexible support scheme was adopted to achieve a primary mirror surface figure accuracy of 9.42 nm; the total mass of the entire telescope (excluding mirrors) is only 3.845 kg. Multi-dimensional finite element analysis was conducted on the telescope under actual working conditions: the strength of the telescope's support materials was verified under self-weight and 10G gravity loads; after removing the rigid body displacement of the mirrors using Zernike polynomials, the surface deformation of the primary mirror was controlled within 1/30 wavelength. In the thermal stability analysis, the structural deformation of the telescope under a temperature change of 100 degree celsius was simulated, and key indicators such as eccentricity and tilt between the mirrors all meet the optical design requirements. In the modal analysis, the first-order natural frequency of the telescope reaches 200 Hz under both self-weight and weightless conditions, demonstrating excellent dynamic stability.

Low-surface-brightness (LSB) structures provide critical insights into the hierarchical formation of galaxies and galaxy clusters. The KASI Deep Rolling Imaging Fast Telescope (K-DRIFT) is designed to detect such diffuse features through deep, wide-field optical imaging with a surface brightness reaching $\sim$$30~\rm{mag}~\rm{arcsec}^{-2}$. To interpret the observation data expected from K-DRIFT, we have developed the Galaxy Replacement Technique (GRT), an $N$-body simulation framework optimized for tracing the gravitational evolution of stellar components. The GRT works by inserting high-resolution galaxy models, including a dark matter (DM) halo and stellar disk, in place of multiple low-resolution DM halos in the base $N$-body cosmological simulation. It allows us to achieve very high mass ($m_{star}=5.4\times10^4\msun\ h^{-1}$) and spatial resolution (10~$\rm{pc}~h^{-1}$) with shorter computation time compared to full hydrodynamic cosmological simulations. Therefore, this technique is particularly well-suited for studying LSB structures, with a surface brightness reaching $\sim$$31~\rm{mag}~\rm{arcsec}^{-2}$. In this paper, we present the motivation and methodology of the GRT, summarize key results from previous studies, and highlight its synergy with K-DRIFT observations. We further discuss planned science cases using the GRT, aiming to build a theoretical basis for interpreting LSB features in various environments.

We examine to what extent the galaxy spin parameter proxy ($\lambda_R$) is affected by bar formation and how it is related to the strong and weak classifications of stellar bars. By creating mock observations of a simulated galaxy, we show that the emergence of a stellar bar can cause mass-weighted $\lambda_R$ to decrease by up to 16%, depending on the bar's orientation. This decrease can be exaggerated if there is a burst of star formation due to the bar driving gas to the center of the galaxy. We use the SAMI galaxy survey to show that weakly barred galaxies have statistically significant younger average stellar populations, higher galaxy spin proxy and higher specific star formation rates compared to strongly barred galaxies within one effective radius. If we consider galaxies with average light-weighted stellar population age less than 3 Gyr within one effective radius, we still find weakly barred galaxies to have a higher galaxy spin proxy than strongly barred galaxies. Based on these trends found from the SAMI galaxy survey, we suggest weakly barred galaxies are rapidly forming, similar to the bar formation process seen in simulations, while strongly barred galaxies are undergoing slower (secular) evolution.

Red supergiants (RSGs) are crucial for studying the properties and evolution of massive stars. It is representative to conduct a census of RSGs across the Local Group, which spans a broad metallicity range. However, identifying RSGs in distant and metal-poor galaxies remains challenging mainly due to contamination of foreground dwarfs and observational limitations. In this work, we perform PSF photometry on publicly released JWST/NIRCam images of five Local Group galaxies: NGC 6822, Sextans~A, NGC 300, WLM, and IC 1613 using the DOLPHOT NIRCam module. We find an optimal color-color diagram (CCD) for metal-poor environments, that is F115W $-$ F200W versus F356W $-$ F444W, which clearly separates RSGs from foreground dwarfs. By using the CCD, we identify 208, 135, and 22 RSG candidates in NGC 6822, Sextans A, and NGC 300, respectively, free from contamination by foreground dwarfs and oxygen-rich asymptotic giant branch stars (O-AGBs). In addition, 40 and 14 RSG candidates are directly selected on the CMD in WLM and IC 1613, respectively. Compared with previous works, the number of RSG candidates within the same luminosity range and sky region increases significantly, demonstrating the advantages of JWST in constructing a more complete RSG sample in the Local Group thanks to its high spatial resolution and photometric quality. In addition, catalogs of O-AGBs and carbon-rich AGBs (C-AGBs) are provided as by-products.

We report the discovery of CGG-z7, the most compact galaxy group at $z\gtrsim7$ identified to the north of the GOODS-North field, observed by the JWST POPPIES program. The system consists of at least six members within a projected size of $7.8\times5.7$ kpc$^2$, four of which are spectroscopically confirmed via [O III] and H$\beta$ emission. The group exhibits a low line-of-sight velocity dispersion ($\approx93.7$ km s$^{-1}$) relative to its substantial stellar mass ($M_* \approx 10^{9.8} M_{\odot}$), yielding a stellar-to-dynamical mass ratio of $M_*/M_{\mathrm{vir}} \approx 0.15$. This ratio, exceeding typical values for virialized halos by a factor of $3$, indicates that the system is highly likely not in dynamical equilibrium. Instead, we interpret CGG-z7 as a pre-virialized structure, likely a major merger caught near apocenter -- capturing the rapid, chaotic formation of a massive "Red Nugget". Spectroscopic analysis reveals extreme ionization conditions and low metallicity across the group. In particular, the central galaxy reaches an extraordinarily high [O III]/H$\beta$ ratio of $\sim18$, which is likely indicative of an obscured AGN. CGG-z7 thus serves as a unique laboratory for the physics of pre-virialized galaxy assembly, bridging the gap between turbulent high-$z$ assembly and the quiescent galaxies seen at cosmic noon.

Recent obeservations have revealed significant soft X-ray polarizations from several quiescent magnetars, including the intriguing $90^°$ polarization angle (PA) swing as a function of photon energy for some sources. We present a general semi-analytical framework for calculating energy-dependent soft X-ray polarization signatures from magnetars, consistently incorporating both QED vacuum resonance in the atmosphere and resonant Compton scattering (RCS) in the magnetosphere. Starting from the polarized radiative transfer equation for RCS and treating vacuum-resonance-induced mode conversion as an input, we employ a first-order approximation in RCS optical depth to evaluate the effect of different magnetospheric plasma density (which depends on magnetic twist), drift velocity and temperature, and viewing geometry on the observed radiation. Our analysis reveals that magnetic twist and plasma drift velocity are the critical parameters controlling the impact of RCS on both the absolute polarization degree and its variation across the soft X-ray spectrum. We find that sufficiently strong RCS can wash out the PA swing caused by vacuum resonance. Furthermore, in addition to the QED vacuum resonance effect, significant relativistic signatures arising from plasma drift velocity ($\beta_0 \gtrsim 0.5$) may introduce an extra $90^\circ$ PA swing in the spectrum. Our calculation framework, based on single-scattering approximation, bypasses the need for complex, multi-dimensional Monte Carlo simulations, providing an analytical pathway for modeling full-surface emission and rotational-phase-resolved radiation from magnetic neutron stars, in support of current and future X-ray polarization missions.

In this review we outline our knowledge on slow neutron captures, concentrating on its main part occurring during the final stages of stellar evolution for low or intermediate-mass stars when they evolve during the Asymptotic Giant Branch, or AGB, stars. We focus our attention on how, in this field, studies passed from a first era of inquiries based on nuclear systematics, to numerical nucleosynthesis computations performed in stellar codes. We then discuss how these last were forced, by observational constraints, to almost abandon, for the synthesis of nuclei between Sr and Pb, the rather naturally activated Ne22 neutron source (operating efficiently at T > 30 keV, and producing a neutron density N_n > 5 10^8 cm^-3). This implied considering the alternative reaction 13C(alpha,n)16O, that can be activated locally after each of the TDU. The mentioned crucial reaction occurs at T< 8 keV, in the time intervals separating two subsequent thermal pulses (TP). The layers where 13C(alpha,n)16O operates are characterized by a radiative equilibrium and their low temperature also yields low values for the neutron density (N_n < 10^7 cm^-3).

We report new high-precision observations of the polarization of light scattered from the atmosphere of Venus, made 100 years after the pioneering studies by Bernard Lyot. The new observations include disk-integrated observations in a range of filters as well as imaging polarimetry. We compare the new results with past observations and models. We have reproduced the 1974 modelling of the Venus polarization by Hansen and Hovenier using modern radiative transfer codes. We show that the new models are in good agreement with the originals, and enable us to calculate the polarization for wavelengths not covered by the original study and to model the polarization distribution across the disk. The new observations are in good agreement with past determinations of the size distribution of the predominant particle mode. They agree with past studies in showing variability of the phase curve between synodic cycles and also polarization variability on short timescales, particularly at higher phase angles (crescent phases). Imaging polarimetry observations show good agreement with models for the redder wavelengths. However, observations in the ultraviolet show very different polarization behavior in the polar regions (within about 30 degrees of the north and south poles). The simplest explanation of this result is that there is a larger Rayleigh scattering component in the polar regions than in the equatorial and mid-latitudes and this could be explained by a lower cloud-top height in agreement with previous spacecraft observations. These ultraviolet polarization observations are inconsistent with horizontally homogeneous atmospheric models.

The potential of prebiotic chemistry to propagate on an exoplanet fundamentally depends on whether the atmospheric conditions can facilitate the production of prebiotic feedstock molecules. Photochemical simulations of exoplanet atmospheres can be used to explore this potential atmospheric synthesis, but require a comprehensive chemical network. We present the implementation of the CRAHCN-O network, constructed to simulate the formation of feedstock molecules such as HCN, H$_2$CO, and simple hydrocarbons, into the VULCAN photochemical kinetics code. We investigate the production of feedstock molecules driven by M-star radiation and compare these to predictions by the N-C-H-O network in VULCAN, for N$_2$-dominated atmospheres with C/O ratios between 0.5-1.5. Predicted abundances are similar for C/O${=}$0.5. Once CH$_4$ is included (i.e., for C/O${>}$0.5), the abundance profiles diverge in the photochemical regions. By analysing the attenuation of UV radiation, we find that hydrocarbon photochemical shielding causes the diverging profiles. CRAHCN-O accumulates C$_2$H$_6$, while N-C-H-O accumulates C$_4$H$_3$ and C$_3$H$_4$. Importantly, C$_2$H$_6$ is photochemically active whereas C$_4$H$_3$ and C$_3$H$_4$ are assumed inactive. With mixing ratios up to a few percent in CRAHCN-O, C$_2$H$_6$ shields CH$_4$ and CO$_2$ from photodissociation and weakens the destruction of HCN and H$_2$CO. Maximum HCN mixing ratios reach 1000 ppm with CRAHCN-O compared to only 3 ppm with N-C-H-O. Other feedstock molecules like HC$_3$N and C$_2$H$_2$ form more efficiently in N-C-H-O. The shielding mechanism and its impact on feedstock molecules persist for radiation from distinct M-star types. These results demonstrate the crucial role of chemical kinetics in understanding prebiotic processes in exoplanet atmospheres, including important considerations for the construction and applicability of chemical networks.

We report optical time-resolved photometric observations of a newly-discovered outbursting system, ASASSN-25dc. Its 8-mag amplitude, 40-day duration, 1-mag dip in the outburst plateau, and positive superhumps are characteristic of a dwarf nova superoutburst in a non-magnetic cataclysmic variable. We establish its stage-A and stage-B superhump periods as 0.059387(5) d and 0.058864(3) d, respectively. The negative superhump period derivative ($-$1.4(2)$\times10^{-5}$ cycle$^{-1}$) during the stage-B superhumps and the empirical relation indicate the mass ratio is 0.054(7), below the period bounce range. The long outburst decline timescale (35.2(1) d mag$^{-1}$) and small superhump amplitude ($\simeq$0.08 mag) observed in ASASSN-25dc are also seen in some period bouncer systems, but not seen in systems well before the period bounce. Despite its short superhump period and indicated small mass ratio, we find no evidence of the excitement of the 2:1 tidal resonance. Moreover, its outburst rise timescale (1.62(9) d mag$^{-1}$) is significantly longer than those measured at less than 0.4 d mag$^{-1}$ in other dwarf nova outbursts around the period minimum. Overall, an inside-out dwarf nova outburst from a massive disc in a system with a mass ratio around or even below the period minimum, but lacking the 2:1 tidal resonance, may explain all these observations. However, this challenges the existing models of dwarf nova superoutbursts, which do not predict these outburst properties in low-mass-ratio systems.

We propose that the recently discovered luminous hypersoft X-ray sources can be explained by accretion onto a binary black hole in a hierarchical triple system. For black hole masses $\sim 15 M_\odot$, the orbital separation of the internal binary might be $\sim 0.01 $~AU. If the donor provides $\gtrsim 10^{-8} M_\odot$ yr$^{-1}$, then the circumbinary accretion disc can explain the observed properties of the most luminous supersoft sources.

Event horizon-scale observations of low-luminosity black hole accretion flows favor magnetically arrested disks, characterized by dynamically important magnetic fields ($\beta\lesssim1$, where $\beta$ is the ratio of plasma thermal pressure to magnetic pressure) and a two-temperature transrelativistic plasma. Motivated by plasma conditions in the synchrotron-emitting regions of these models, we perform 2D particle-in-cell simulations of electron-ion plasmas with a realistic mass ratio, subject to continuous compression perpendicular to the mean magnetic field $\boldsymbol{B}_0$. Conservation of particle magnetic moments drives pressure anisotropy $P_{\perp}>P_{\parallel}$, triggering anisotropy-driven instabilities. For ion plasma beta $\beta_{i0}=0.5$ and ion temperature $k_{\text{B}}T_{i0}/m_i c^2=0.05$, the ion pressure anisotropy is regulated by the ion cyclotron instability, while the mirror mode influences the late-time electron anisotropy. Both species develop nonthermal components at high energies, consistent with stochastic acceleration by cyclotron-scale fluctuations. We characterize how the onset and time evolution of the plasma instabilities, as well as the resulting ion and electron anisotropies and energy spectra, vary with $\beta_{i0}$, $k_{\text{B}}T_{i0}/m_i c^2$, electron-to-ion temperature ratio $T_{e0}/T_{i0}$, and the compression rate. Increasing the thermal energy toward relativistic values raises the anisotropy thresholds for all instabilities observed in our simulations, allowing larger anisotropies to develop. For $T_{e0}/T_{i0}<1$, as expected in collisionless two-temperature accretion flows, the growth of mirror and whistler instabilities is delayed or suppressed, leading to increasingly adiabatic evolution of the electrons. Our findings can be used to inform global fluid models of black hole accretion.

The embedded disks surrounding individual Class 0 protostars are structureless. Disks surrounding Class I stars may be continuous or have a ring-gap substructure, whereas all disks around Class II stars have a ring-gap substructure that gradually disappear as the disks evolve into debris disks. This common sequence in young lone stars requires an explanation. This study aims to show that the physical model Reaction-Diffusion Systems with Moving Reaction Front can be used to describe and classify protostellar disks according to their structure. A comprehensive review of observations made with the ALMA radio telescope shows: first, that the protostar-disk system presents a geometry analogous to that of an reaction-diffusion system with two separate compartments, namely, protostar and disk. Second, that in the protostar, matter is processed at high temperature, resulting in a chemical composition different from that of the disk. Third, that the equatorial outflow emitted by the protostar, rich in highly reactive trihydrogen cation, acts as a moving reaction front, MRF, that triggers the formation of molecules and nuclei in the disk. The time lag of nucleation with respect to the passage of the MRf would be the cause of the formation of the gaps between the rings of particles that form in the disk. The MRF is a transient phenomenon and its passage causes the transformation of a continuous disk, Class 0, into a disk with a ring-gap structure, Class II, whose temporal evolution begins at the interface of the star and moves outwards.

As one of the most widely recognized high-quality standard stars, the Stetson standards have been extensively used as a photometric reference for calibrating other surveys. In this work, we present an independent validation and re-calibration of the Stetson standard star photometry using the BEST database. Based on typically 30,000-70,000 calibration stars per band, we find that the original Stetson photometry achieves field-to-field zero-point precisions of approximately 10--40\,mmag in the $UBVRI$-band. In addition, significant spatially dependent magnitude offsets are detected within individual Stetson fields for all bands, with magnitudes exceeding 1\%, probably caused by the calibration errors in the Stetson photometry. After correcting those systematic errors, the agreement between the Stetson and BEST photometry is improved to $\sim$5\,mmag for individual fields for $BVRI$-band. The re-calibrated photometry is further validated using the SCR standards, yielding agreement better than 10\,mmag for individual stars in the $BVRI$ bands and confirming zero-point precisions of 2--4\,mmag in the $BVI$ band. The precisions is further confirmed by checks using Gaia DR3 broad band colors. These results highlight the power of the BEST database for improving photometric calibration and suggest that, if feasible, it be incorporated into the calibration process of future releases of the Stetson standard catalog.

There is an ongoing discussion in the literature on the nature of long-period transients (LPTs), radio-emitting sources with periods ranging from hundreds to tens of thousands of seconds. Although some of these objects have been identified as white dwarf (WD) + M-dwarf binaries, this description currently does not fit the entire class. An example is GLEAM-X J162759.5-523504.3 (hereafter GLEAM-X J1627-5235), with a period of 1091 s, for which the lack of an optical counterpart disfavors the presence of such a binary system. In this case, GLEAM-X J1627-5235 could be interpreted as an isolated, massive, fast-rotating, and highly magnetized (~ 1e+9 G) WD pulsar. Its properties are consistent with a carbon-oxygen WD of mass ~1.3 Msun and radius ~2500 km, possibly supported by small-scale multipolar magnetosphere structures that keep it above the death line for WD-pulsars. We assess a double WD merger origin, modeling the post-merger rotational evolution under accretion, propeller, and magnetic braking torques. We find rotational age of ~572 Myr for GLEAM-X J1627-5235, i.e., the post-merger time required to reach its observed period. This result is consistent with current optical upper limits for GLEAM-X J1627-5235 and support the WD pulsar interpretation for this source. We also discuss how the same model can apply to other LPTs.

Globular clusters (GCs) exhibit varying degrees of flattening (ellipticity), which may provide insight into their internal dynamics and evolution histories. Commonly used methods to measure ellipticity, such as ellipse fitting of density contours and principal component analysis, often produce biased results, especially for clusters that are nearly round or have few observable stars. Using a combination of ground-based and space-based photometry, we investigate the shapes of 29 Galactic GCs. To that end, we test two commonly used methods: an ellipse fit to a kernel density profile and a principal component analysis. We find that both methods suffer from bias arising when the number of stars is small or the cluster is close to round. To solve this issue, we develop a robust method to measure the ellipticity of GCs, test it extensively on mock data, and apply it to the 29 Milky Way GCs in our sample. Using the $V/\sigma$ diagram used in the isotropic oblate rotator framework, we examine potential causes for the flattening, including rotation and velocity anisotropy. For ten clusters: NGC~104, NGC~1261, NGC~2808, NGC 3201, NGC 5286, NGC 5904, NGC 5986, NGC 6205, NGC 6341, and NGC 7078 we identify a very good agreement between the rotation angle and semi-minor axis of the ellipse, further corroborating the findings that rotation is the main driver of the ellipticity. The $V/\sigma$ diagram reveals that velocity anisotropy or tides could also be important in shaping the GCs. The robust method developed provides reliable measurements of the ellipticity of GCs, emphasising the importance of taking into account the flattening in theoretical models and simulations. It also offers a promising way to investigate the shapes of multiple stellar populations within GCs, where only small samples are usually available.

In recent years, massive star cluster environments have proved to be bright sources of very-high energy gamma-rays, in particular young clusters which are powered by the winds interacting in their cores. In order to understand how these winds can accelerate particles up to very-high energies, it is necessary to model their interactions from small (sub-pc) to large (10s of pc) scales over several millions of years. A key open question concerns the structure and properties of the resulting wind termination shock. By performing 3D magnetohydrodynamic simulations of clustered winds embedded in a superbubble cavity, we demonstrate that the dynamics of stellar wind interactions and the resulting shock structure solely depends on the density and pressure of the cavity. This implies that the initial conditions of the simulation can be tuned in order to simulate star clusters of arbitrary age at a reduced computational cost. This novel method is validated using a toy cluster hosting 30 identical stars. We discuss the properties of the resulting cluster-wind termination shock under various assumptions. In particular, we are able for the first time to obtain a fully decoupled spherical wind termination shock for a 5 Myr old cluster. We further show that radiative cooling increases the sphericity of the shock. In general, the morphology of the outflow depends on the number of dominant stars, on the power of the stars sitting at the edge of the cluster core, and on the compactness of the cluster. We additionally show how a semi-analytical model can be used in order to estimate key morphological properties of the outflow without relying on large-scale simulations.

Eclipsing close double white dwarf (WD) systems provide a unique opportunity to directly constrain hydrogen-envelope retention and test common-envelope (CE) evolution in low-mass stars, since they allow precise determinations of stellar masses and radii. We analyze J2102-4145, an eclipsing binary composed of two low-mass helium-core white dwarfs in a 2.4-hour orbit. By comparing the observed radii and effective temperatures with updated evolutionary models for CE evolution and stable Roche-lobe overflow (SRLOF), we confirm that both stars are helium-core white dwarfs. The primary, with a mass of 0.375 solar masses, is consistent with SRLOF models that retain thick hydrogen envelopes and sustain residual nuclear burning, whereas the secondary, with a mass of 0.314 solar masses, can only be reproduced by CE models in which the hydrogen envelope is almost completely removed. The inferred cooling ages (approximately 220 Myr for the secondary and between about 260 and 510 Myr for the primary, depending on the contribution of residual nuclear burning) support a formation sequence in which the primary formed first through SRLOF, followed by a CE phase that produced the compact secondary. Reconstruction of the CE energy budget yields progenitor and orbital parameters consistent with this evolutionary picture. The unusually small radius of the secondary requires an extremely thin hydrogen envelope, with a mass below about 10e-7 solar masses, well below the values predicted by standard bifurcation criteria. J2102-4145 therefore provides one of the strongest observational constraints on the hydrogen-envelope mass of post-CE low-mass white dwarfs and represents a benchmark challenge for current prescriptions of envelope ejection.

We introduce a new, discretized grain size evolution model, incorporated into the GIZMO code and coupled with FIRE-3 stellar feedback and ISM physics, to investigate variations in dust abundance, chemical composition, and grain sizes observed in the Local Group. This model tracks the size evolution of specific dust species, and includes stellar production of dust, dust growth through gas-phase metal accretion, dust destruction by sputtering, SNe shocks, and astration, grain-grain collisional shattering and coagulation, and turbulent dust diffusion. Using idealized galaxy simulations, we test the dependence of MW dust properties on variations in each dust process and find that our model uniquely predicts a bimodal grain size distribution. This bimodality is due to our simulation's ability to resolve each dust process and where they occur in the ISM, unlike other works. We find that Local Group dust abundances are determined by dust growth and destruction, with little dependence on coagulation or shattering, explaining why models that do not include these processes can match abundance observations. We also find that variations in Local Group extinction curve slopes are determined by coagulation, with inefficient coagulation leading to steeper slopes. However, inefficient coagulation also results in stronger extinction curve bumps, which are not observed. We also do not predict a population of very small (${<}1$ nm) carbonaceous grains, required for MIR emission features, due to their rapid growth by accretion. These results highlight the possible necessity of ``top-down'' PAH formation from preexisting grains as a means to inhibit carbonaceous dust growth.

In ground-based high-contrast instruments, non-common path aberrations (NCPAs) limit detection performance, as they are unseen by the adaptive optics (AO) wavefront sensor but impact the astrophysical image, creating quasi-static speckles. SAXO+, the upgrade of the SAXO (SPHERE AO system) includes a second loop of AO downstream of the SAXO loop that is equipped with a near-infrared pyramid wavefront sensor whose nonlinearities, usually described with modal optical gains, might be challenging for removing quasi-static speckles. We investigated two methods of quasi-static speckle removal : NCPA compensation and a dark hole loop, behind a pyramid AO system, measuring the interest of compensating for the pyramid optical gains. We performed end-to-end numerical simulations under various astrophysical conditions. We offset the pyramid wavefront sensor operating point to apply both the speckle suppression methods, with or without optical gain calibration. We evaluated the performance by measuring the residual starlight in the coronagraph image. A by-product of our study is an on-sky calibration method of measuring the pyramid optical gains. NCPA compensation reduces the residual starlight in the coronagraph image by a factor of 20 for seeing between 0.7" and 1" for a bright star and a factor of 2 at 0.7" for a faint star. Optical gains compensation enhances the performance at poor seeing and small pyramid modulation radius with a bright star, but shows a useless or even negative impact due to estimation inaccuracies at faint targets. On the other hand, the dark hole loop reduces the residual starlight by a factor of 200. The optical gain calibration enhances the dark hole performance behind a single pyramid AO system but is useless behind the SAXO+ system. Our parametric study gives baseline values for the efficient control of the dark hole loop for the SAXO+ system.

Giant Molecular Clouds (GMCs) are dominated by supersonic turbulence, creating a complex network of shocks and filaments that regulate star formation. While the global inefficiency of star formation is well-observed, predicting exactly which gas parcels within a turbulent cloud will collapse to form stars remains a challenge. In this work, we present a supervised machine learning framework to forecast the Lagrangian history of prestellar cores in magnetohydrodynamic (MHD) turbulence. We utilize Extreme Gradient Boosting (XGBoost) to train a regression model on the trajectories of $\sim 2.1$ million tracer particles evolved within a self-gravitating, turbulent MHD simulation. By mapping the instantaneous phase-space state (position, velocity, and density) of gas parcels to their future coordinates, our model successfully predicts the 3D evolution of star-forming cores over a horizon of $\sim 0.45$ Myr ($0.25~t_{\rm ff}$). We achieve a global coefficient of determination of $R^2 > 0.99$ and demonstrate that the model captures the non-linear convergent flows characteristic of gravitational collapse. Crucially, we show that local phase-space information alone is sufficient to distinguish between transient density fluctuations and bound collapsing cores. This data-driven approach offers a computationally efficient alternative to traditional sink-particle algorithms and provides a pathway for developing high-fidelity subgrid models for galaxy-scale simulations.

On the Sun, the inertial mode with the largest observed amplitude (rms velocity exceeding $10$ m/s) is the high-latitude mode with longitudinal wavenumber $m=1$. In two dimensions, on the sphere, linear theory predicts that this mode is unstable due to a shear instability associated with latitudinal differential rotation (fast equator, slower polar regions). We investigate the evolution of this instability numerically and theoretically. The nonlinear vorticity equation is solved using direct numerical simulations in the time domain. The only control parameter is the Ekman number $E$. For $10^{-3}\lesssim E< E_c \approx 1.5\times10^{-3}$, only the high-latitude $m=1$ mode is unstable. We extract its saturation amplitude as a function of $E$ and compare the results with predictions from two perturbative approaches in nonlinear stability theory. The simulations reveal a supercritical Hopf bifurcation. Near onset, the mode amplitude is well described by the Landau equation $d|A|/dt=\sigma_I |A|+\beta_I |A|^3$, with a positive linear growth rate $\sigma_I$ and a negative nonlinear coefficient $\beta_I$. The coefficient $\beta_I$ depends weakly on $E$, implying that the saturated amplitude scales approximately as $|A|\propto\sigma_I^{1/2}$. The equilibrium mode contains the $m=1$ fundamental and harmonics $m=2$ and $m=3$, whose amplitudes scale as $\sigma_I^{m/2}$. Saturation results from Reynolds stresses that smooth the latitudinal differential rotation. For $E=4\times10^{-4}$, consistent with solar-like turbulent viscosity, the saturated velocity reaches $28$ m/s, comparable to solar observations. These results should be interpreted cautiously, since in three dimensions the instability is baroclinic and involves different physics.

The molecular composition inside the dust sublimation zones of protoplanetary disks is mostly unknown but important to understanding terrestrial planet formation. A few molecules have been observed from this region, specifically CO, H2O, OH and SiO. The small surface area makes observing this region difficult, hence modeling is required to disentangle the innermost disk from regions further out. We model a protoplanetary disk around a Herbig-type star including the dust depleted inner region (approx. 0.1-0.3 au) and aim to investigate the chemistry of this region and explain existing and future observations. Methods. We post-process the dust and gas distribution of a magnetohydrostatic model with the radiation thermochemical code ProDiMo to study the chemistry and to produce observables. We find that the dust free inner disk is a molecular rich environment, where besides CO we also find H2, H2O and SiO. The gas temperature profile is complex and fluctuates between 700 and 2000 K, which is warm enough to produce CO overtone line emission. Next to the CO overtone lines we also find strong high J-level fundamental CO lines between 4.3 and 4.6 micron. The elemental enrichment of Si due to dust sublimation leads to 2 orders of magnitude more SiO abundance. The SiO gas has average temperatures of approx. 1000 K resulting in strong SiO overtone emission in the spectral range between 4 and 4.3 micron. We predict that the gas density in the dust depleted inner disk is high enough to allow for H2 formation, resulting in an molecular rich environment. For our representative Herbig model, the dust-depleted inner disk is responsible for at least 90% of the line emission for CO and H2O between 1 and 28 micron. Next to CO overtone lines, SiO overtone lines are expected to be an important tracer of a dust free inner disk.

Most long gamma-ray bursts (GRBs) exhibit afterglows broadly consistent with external forward-shock emission, typically described by smooth broken power-law decays in the multiband light curve. However, a minority of well-sampled GRBs deviate from this behavior, including GRB 250129A. This object shows multiple late-time rebrightenings at X-ray and optical wavelengths. Rebrightenings are often attributed to energy injection from prolonged central engine activity, refreshed shocks from delayed shell collisions, density jumps in the ambient medium, or angular jet structure and viewing-angle effects. After analysing the prompt emission observed in gamma-rays and the near-infrared, we investigate the origin of X-ray and optical flaring episodes in GRB 250129A. Physical processes in the afterglow light curves were investigated using methods ranging from empirical fitting to Bayesian inference. The well-sampled flares and the connection between the prompt and afterglow emission allow us to test the consistency of the fireball model and alternative scenarios. Conducting the prompt and time-resolved analyses, we obtained an isotropic-equivalent energy of E_iso,gamma = (1.35 +/- 0.12) x 10^53 erg. By modeling the afterglow using an agnostic Bayesian framework (NMMA), we rule out both a single external-shock evolution and a one-time energy-injection scenario. Numerical calculations show that the rebrightening episodes are consistent with refreshed shocks from delayed collisions between relativistic shells. Based on the consistency between our analyses of the prompt and afterglow GRB 250129A data, we find that two statistically significant rebrightening episodes occur within 1.1 days post trigger and can be explained by a sequence of refreshed shocks. Temporally and spectrally rich GRB datasets such as the one presented in this work, provide a powerful means to test current modeling frameworks.

We present 5 GHz e-MERLIN radio images at 50 mas resolution of the nuclear regions of the Legacy e-MERLIN Multi-band Imaging of Nearby Galaxies survey (LeMMINGs), the deepest statistically complete radio-band survey of the local Universe (<120 Mpc), consisting of 280 galaxies spanning all morphological and nuclear types. We detect nuclear radio emission above a median 5 sigma threshold of 0.33 mJy beam^-1 in 68 of 280 sources (24 percent), with core luminosities in the range 10^35 to 10^41.9 erg s^-1. The radio emission is attributed to active galactic nuclei, circumnuclear star formation, or, in the case of NGC 3690, a tidal disruption event. The brightest radio nuclei, with brightness temperatures >=10^6 K, reside in optically active galaxies such as LINERs and Seyferts. The detection rate for inactive systems (H II and absorption-line galaxies), which may host low-luminosity active galactic nuclei, is 8 percent. Most detections (78 percent) are compact (<10 pc), while the remaining 22 percent show extended jet-like features up to 380 pc. Compared to the 1.5 GHz LeMMINGs data, the 5 GHz observations provide superior resolution and spatial filtering, resolving out large-scale structures and isolating genuine nuclear emission. Our results suggest that low-luminosity active galactic nuclei are the primary manifestation of black hole activity in the local Universe in the form of compact jets and cores, with a preference for early-type hosts. The two LeMMINGs campaigns indicate that up to 30 percent of the local galaxy population hosts a radio-active nucleus, highlighting the necessity of high-resolution, high-sensitivity imaging for uncovering nuclear emission at the lowest luminosities.

Nonequilibrium chemistry is central to many astrophysical environments but remains a major computational bottleneck in simulations because solving the associated stiff ODE systems is expensive. Neural surrogates promise large speedups, yet existing studies rarely provide systematic comparisons of architectures or rigorous optimization toward both accuracy and efficiency. We introduce CODES, a principled framework for optimizing and benchmarking astrochemical surrogate models. Using CODES, we compare four neural surrogate architectures across four KROME-generated datasets spanning primordial and molecular-cloud chemistry with up to 287 reactions across 37 species. Dual-objective optimization reveals pronounced accuracy-efficiency trade-offs across architectures. Fully connected models achieve the highest accuracy and most reliable uncertainty estimates, while latent-evolution models show improved robustness under iterative prediction. Our results highlight the importance of systematic optimization and architectural comparison. The datasets, metrics, and benchmarking procedure are publicly released within CODES to enable reproducible surrogate benchmarking.

We study the full evolution of a 1.313 Msun white dwarf star that descended from a 9 Msun main-sequence progenitor with an initial metallicity of Z=0.02. Using MESA r24.08.01, we calculate its entire evolution from pre-ZAMS to the WD cooling curve, including both the evolution through 139 thermal pulses and the post-AGB phase. The resulting remnant is an ultramassive H-deficient WD, for which the composition, in mass fraction, is 47.7% O16, 39.7% Ne20, 4.2% Mg24, 3.3% Na23 and 0.386% C12 -- corresponding to a total mass of 5 x 10^-3 Msun of C --, surrounded by a 1.5 x 10^-5 Msun He layer. We also investigate the effects of fully suppressing the TP-SAGB stage by adopting a high mass-loss rate only after the second dredge-up, and find only minor differences in the final mass and composition. In addition, we calculate models with and without phase separation during the WD stage, estimating a cooling delay of only 16 Myr. This is the first ultramassive white dwarf sequence for which both the TP-SAGB and post-AGB stages are calculated and, to our knowledge, the most massive WD model from complete evolution for which cooling times and detailed abundance profiles are published

Molecular hydrogen (H$_2$) and hydrogen deuteride (HD) are key coolants in primordial gas and regulate the formation of the first stars and proto-galaxies. Recent results from the James Webb Space Telescope provide striking insights into galaxies detected at high redshifts, which are found to be significantly more abundant and luminous than expected from galaxy formation models, thus suggesting a gap in our understanding of the early Universe. Standard pathways for H$_2$ formation in the early Universe proceed through the H$^-$ and H$_2^+$ intermediates, both of which are strongly suppressed at high redshift by the cosmic microwave background. We propose an additional pathway for H2 and HD formation that could be active as early as the end of the epoch of recombination and could enable the formation of the first stars earlier than the current prediction at redshift z ~ 30 - 20. The proposed pathway relies on the manifestation of Jahn-Teller dynamical coupling between electronic states of H$_3^+$. This coupling induces transient three-body recombination in H$^+$, H and H, and charge exchange within the charged atom-dimer complex that directly creates ground-state H$_2$ (and HD), bypassing the fragile intermediates that limit the standard primordial pathways. Our analysis shows that this mechanism could occur under the thermodynamic conditions of the post-recombination epoch, also suggesting that it might be playing a role in the active galactic nuclei feedback processes, regulating the formation rates of the first stars and the accretion rates of the first black holes. Though the global impact on galaxy formation and black-hole growth is not yet determined and will require quantitative assessment in future modeling, the mechanism offers an additional chemical route for H$_2$ and HD formation, with substantial cosmological relevance for primordial chemistry and early structure formation.

The observation of high-energy neutrinos from the direction of the nearby active galaxy, NGC 1068, was a major step in identifying the origin of high-energy cosmic neutrinos. The multimessenger data imply that high-energy neutrinos originate from the hearts of active galaxies which are opaque to GeV-TeV $\gamma$-rays. This realization is reinforced by the excess of neutrinos in the direction of NGC 4151 and Circinus Galaxy, other nearby active galactic nuclei (AGNs). Understanding the vicinity of supermassive black holes with electromagnetic radiation is often challenging due to uncertainties associated with the absorption of emission in these dense environments, and neutrinos can be used as a powerful probe of the inner parts of the active galaxies. Considering the five brightest neutrino-active galaxies, NGC 1068, NGC 4151, CGCG 420-15, Circinus Galaxy, and NGC 7469, we employ the measured neutrino spectra together with the sub-GeV $\gamma$-ray emission measured by the {\em Fermi} satellite to break the degeneracy and narrow in on the parameter space of neutrino emission from turbulent coronae of AGNs. We also study contributions of jet-quiet AGNs, whose properties are similar to NGC 1068 and NGC 7469, to the isotropic neutrino background flux, through exploring possibilities that the neutrino luminosity function may deviate from the X-ray luminosity function. Our results will help estimate the prospects for identifying additional neutrino-active galaxies and guide future targeted analyses.

Context. Ultra fast outflows (UFOs) from active galactic nuclei (AGN) are expected to lead to the formation of sub-relativistic strong shocks expanding in a dense circumnuclear medium, and thus have the potential for being efficient particle accelerators, and to be proficient sources of gamma rays and neutrinos. Aims. We investigate the detectability of a sample of nearby identified UFOs in gamma rays and neutrinos with current and next- generation instruments. Methods. We model the acceleration of particles at the strong shocks of UFOs, and estimate the associated gamma-ray and neutrino signal. We adopt our model to investigate the prospects for detection with current and next-generation observatories. Results. We find that several UFOs could be detectable in the very-high-energy (VHE) domain - for example, by the Cherenkov Telescope Array Observatory (CTAO)- even if they remain undetected by Fermi-LAT in the high-energy range. Detectability is favored for hard proton spectra (spectral index {\alpha} \lessim 3.9), high acceleration efficiencies, and amplified magnetic fields. Our results suggest that next-generation VHE observatories could detect the first gamma-ray signatures of AGN UFOs, providing a new probe of particle acceleration in sub-relativistic shocks

Recent advances in quantum sensing -- optical clocks at $5.5\times 10^{-19}$ systematic uncertainty, frequency-dependent squeezing below the standard quantum limit, quantum magnetometers approaching fundamental sensitivity limits -- raise a natural question: can these technologies detect gravitational waves directly, or enhance existing detectors beyond current capabilities? We show that the answer is primarily determined by the \emph{coupling mechanism} between the gravitational wave and the sensor. Starting from the tidal Hamiltonian in Fermi normal coordinates, we identify three physically distinct coupling mechanisms and derive their transducer gains within linearized general relativity and non-relativistic quantum mechanics. Internal atomic coupling (tidal distortion of electronic wavefunctions) yields a transducer gain $G_A = 2.4\times 10^{-20}$, with vanishing first-order energy shifts for all $J=0$ clock states -- a $\sim\!10^{35}$ deficit relative to laser interferometry that exceeds any projected quantum enhancement. Center-of-mass coupling (Doppler shifts from geodesic motion) reaches strain sensitivities of $\sim\!10^{-18}$, still $10^4$ above LISA requirements. Light propagation coupling (phase accumulation over macroscopic baselines) provides the enormous transducer gain that makes laser interferometry -- and atom interferometry -- viable. For detectors exploiting this third mechanism, we quantify how much improvement quantum sensors can provide through the detector's noise architecture: LISA's noise budget is $\sim\!91\%$ classical, limiting combined quantum enhancement to $\mathcal{E} \approx 1.04$, while ground-based detectors in the shot-noise-dominated regime achieve $\mathcal{E} = 1.8$--$2.4$. Atom interferometers exploit the same light-propagation mechanism to uniquely access the 0.01--10~Hz band.

Spectral-sirens inference enables the extraction of cosmological parameters from gravitational-wave data alone, without electromagnetic counterparts or galaxy catalogs. We introduce new parametric mass functions for the binary black hole population that capture significant structure across the mass spectrum and are moderately favoured by Bayesian evidence over simpler models. Analysing the latest gravitational-wave transient catalog, GWTC-4.0, we show that powerlaws-only population models constrain the Hubble constant to $23\%$ precision, $H_0 = 53.3^{+14.0}_{-10.8} ~\rm km \,s^{-1} \,Mpc^{-1}$ at $68\%$ confidence level. This represents a $\sim 50\%$ improvement over the corresponding binary black hole-only analysis by the LIGO-Virgo-KAGRA collaboration, achieving precision comparable to their joint analyses including neutron stars and galaxy catalogs. We further test alternative cosmological models, establishing competitive constraints on modified gravitational-wave propagation, while bounds on the dark energy equation-of-state parameters remain uninformative. Projecting to future O5 observing run, we forecast substantial improvements in $H_0$ and modified propagation parameters with larger datasets at higher redshifts. Our results highlight the strong interplay between the black hole mass distribution and inferred cosmology, underscoring the need for suitable population models to fully exploit gravitational-wave data.

This study presents the trajectory design for a mission touring Saturn's Inner Large Moons (Rhea, Dione, Tethys, Enceladus, and Mimas) engineered to meet observational requirements, including full surface coverage, while ensuring low fuel consumption and compatibility with current power and propulsion technologies (radioisotope thermoelectric generators and Hall effect thrusters). The tour begins at Rhea and ends at Mimas, using a trajectory concept that alternates between extended observation phases around each moon and Saturn centered low-thrust spiral arcs to transition efficiently to the next target. The J2-perturbed Circular Restricted Three-Body Problem is adopted to design exploration paths, with halo orbits serving as staging points for heteroclinic and homoclinic loops that enable prolonged, repeated, and comprehensive surface reconnaissance (including critical regions such as Enceladus poles, where geological activity produces intense plumes). Stable and unstable hyperbolic invariant manifolds of the halo orbits act as departure and arrival gateways for propelled inter-moon transfers, modeled in an ephemeris-based framework including gravitational perturbations from the moons, the Sun, and Saturn's oblateness. The dynamical model setup is guided by a rigorous perturbation analysis to maximize computational efficiency while maintaining high fidelity trajectory design. A locally optimal guidance law minimizes propellant consumption. The proposed tour offers an alternative to traditional flyby missions, providing comparable total duration but greater observing time and reduced fuel requirements, and advances previous work by achieving both complete lunar surface coverage and high-fidelity modeling.

We adopt a general relativistic ray-tracing approach to study the shadows and polarization images of spherically symmetric Gauss-Bonnet (GB) black holes enveloped by geometrically thick accretion flows. Specifically, we adopt a phenomenological RIAF-like model and an analytical Hou disk model. In the RIAF-like model, increasing the GB coupling parameter $\lambda$ reduces both the size and brightness of the higher-order image, while increasing $\theta$ alters the shape of the higher-order image and obscures the horizon's outline. The main difference between isotropic and anisotropic emission is that the latter produce distortion of the high-order image in the vertical direction, leading to an elliptical morphology. For the Hou disk model, due to specific regions being geometrically thinner with the conical approximation, the high-order images are narrower with the increase in $\lambda$ than the RIAF model. While increasing $\theta$ enhances the brightness of the direct images outside the higher-order images, but hardly changes the size of the higher-order images, which is in sharp contrast to the RIAF model. Meanwhile, the Hou disk produces polarization patterns that trace the brightness configuration and are affected by $\lambda$ and $\theta$, reflecting the intrinsic structure of spacetime. These results illustrate that intensity and polarization in thick-disk models provide probes of GB black holes and near-horizon accretion dynamics.

Previous work on the gravitational-wave background generated in a two-scalar-field cosmological model, in which warm inflation and the dark sector are unified within a single framework, has shown that the resulting spectrum could be potentially detectable by planned next-generation gravitational-wave observatories. In this work, we extend this analysis to the weak dissipation regime of warm inflation, highlighting how the features of the inflationary scenario play a crucial role in the production of gravitational waves. The full gravitational-wave energy spectrum is calculated using the formalism of continuous Bogoliubov coefficients. By comparing our results with those obtained in the strong dissipation regime and with the sensitivity curves of future detectors, we find that the weak dissipation regime improves the prospects for observational detection.

Energy dissipation and particle acceleration in the collisionless magnetotail plasma remain incompletely understood. While Kinetic Alfvén Waves (KAWs) are widely hypothesized to mediate these processes, observational characterization of their spectral properties and dissipation signatures in magnetotail boundary layers remains limited. We report observations of KAW turbulence and parallel electric fields ($E_{\parallel}$) in the outer Plasma Sheet Boundary Layer (PSBL) using high-resolution burst-mode data from the Magnetospheric Multiscale (MMS) mission. For a crossing event on May 31, 2017, we identify broadband KAW turbulence characterized by a normalized electric-to-magnetic field ratio $\mathcal{R} = |\delta E_{\perp}|/(v_A|\delta B_{\perp}|) = 2.5 \pm 1.2$ exceeding the MHD limit, a spectral break near ion scales, a steep kinetic-range spectral slope ($\alpha = -3.48 \pm 0.13$), and low magnetic compressibility ($C_{\parallel} \approx 0.03$). We observe impulsive parallel electric field structures (up to 15~mV/m) and large-amplitude density fluctuations (up to 68\%) during intervals of enhanced wave activity. The steep spectral slope, steeper than theoretical predictions for undamped KAW cascades ($-7/3$ to $-8/3$), is consistent with substantial energy removal from the cascade at kinetic scales. The near-zero correlation between the $E_{\parallel}$ waveform and density fluctuations ($r \approx -0.03$) suggests that the observed $E_{\parallel}$ structures are not straightforwardly organized by compressive density variations, consistent with dissipation through direct wave--particle interaction. Attribution to a specific damping channel (e.g., Landau damping) is not uniquely constrained by the present diagnostics. These observations support collisionless damping of KAW turbulence at kinetic scales in the intermediate-beta, outer PSBL of the terrestrial magnetotail.

We show that the hidden symmetry and separability characteristic of Kerr geometry arise as an inevitable \emph{local} consequence of the Einstein equations for rotating spacetimes. Without assuming separability, algebraic speciality, Killing-Yano symmetry, or global boundary conditions, we analyze general stationary and axisymmetric geometries in a locally non-rotating orthonormal frame. Imposing a minimal physical requirement, the local equilibrium condition, we find that the mixed Einstein equations enforce a rigid projective alignment between the radial and angular sectors. This result does not rely on vacuum: the mixed equations are insensitive to the detailed form of a comoving stress-energy tensor. Consistency then requires equality of Schwarzian derivatives, leading to a universal classification of local solutions into Möbius, exponential, and trigonometric branches. Global regularity generically excludes the trigonometric branch, leaving precisely the Kerr-type sector and the emergence of Petrov type D structure. This provides a structural precursor to Kerr uniqueness, demonstrating that the kinematical core of Kerr geometry is encoded locally within the Einstein equations.

Previous studies have established that a weak-to-strong transition occurs in Alfvenic magnetohydrodynamic (MHD) turbulence as energy cascades from large to small scales. However, the spatiotemporal (frequency-wavenumber) properties of compressible MHD turbulence involving all eigenmodes, which encode the strength of nonlinear interactions, remain difficult to characterize observationally. Consequently, whether a similar weak-to-strong transition occurs in compressible turbulence remains elusive. Using a novel multi-spacecraft, polarization-based mode-decomposition technique with measurements from the Cluster spacecraft in Earth's magnetosheath, we obtain spatiotemporal power spectra of all MHD eigenmodes and present the first quantitative assessment of nonlinear frequency broadening. Our results show that slow modes exhibit a weak-to-strong transition, evolving from wave-like peaks to frequency-broadened spectra as nonlinearity increases, whereas fast modes remain weakly turbulent with narrow peaks near their eigenfrequencies. Both Alfvenic and compressible fluctuations contribute significantly to low-frequency, large-scale quasi-two-dimensional structures. These findings provide a comprehensive observational characterization of compressible turbulence across mode composition, spatiotemporal scales, and weak-strong turbulence regimes, with implications for energetic particle transport, turbulent dynamos, plasma heating, and solar wind-magnetosphere coupling.

One of the main targets for space-borne gravitational wave detectors is the detection of Extreme Mass Ratio Inspirals (EMRIs). The data analysis of EMRIs requires waveform models that are both accurate and fast. The major challenge for the fast generation of such waveforms is the generation of the Teukolsky amplitudes for generic (eccentric and inclined) Kerr orbits. The requirement for the modeling of $\sim10^5$ harmonic modes across a four-dimensional parameter space makes traditional approaches, including direct computation or dense interpolation, computationally prohibitive. To overcome this issue, we introduce a convolutional encoder-decoder architecture for a fast and end-to-end global fitting of the Teukolsky amplitudes. We also adopt a transfer learning strategy to reduce the size of the training dataset, and the model is trained gradually from the simplest Schwarzschild circular orbits to generic Kerr orbits step by step. Within this framework, we obtain a surrogate model based on a semi-analytical Post-Newtonian dataset, and the full harmonic amplitudes can be generated within milliseconds, while the median mode-distribution error for generic orbits is approximately $\sim10^{-3}$. This result indicates that the framework is viable for constructing efficient waveform models for EMRIs.

Black holes are the most efficient scramblers in nature. By mapping the instantaneous mass and angular momentum of two spinless black holes in a quasi-circular binary onto those of an effective Kerr black hole, we demonstrate that the final state of the merger remnant corresponds with remarkable accuracy to the configuration that renders null geodesics unstable at the highest possible rate. This suggests a deep connection between the properties of black holes resulting from binary mergers and their unstable null orbits.