Papers for Tuesday, Nov 25 2025

Black hole environments often host plasmas that are fully collisionless or contain intrinsically collisionless regions, including relativistic jets and coronae where particle energization is ubiquitous. Capturing the physics of these systems requires numerical methods capable of modeling relativistic, magnetized, collisionless plasmas in strong gravitational fields. In this work, we introduce the general-relativistic module for the Entity -- the first open-source, coordinate-agnostic performance-portable particle-in-cell code. The code enables fast axisymmetric simulations of collisionless plasmas around black holes on any modern high-performance computing architecture (both GPUs and CPUs).

Helium reionization is the most recent phase change of the intergalactic medium, yet its timing and main drivers remain uncertain. Among the probes to trace its unfolding, the 3.46 cm hyperfine line of singly-ionized helium opens the study of helium reionization to upcoming radio surveys. We aim to evaluate the detectability of the 3.46,cm signal with radio surveys and the possible constraints it can place on helium reionization, in particular whether it can distinguish between early and late helium reionization scenarios. Moreover, we perform a comprehensive study of the advantages of single-dish vs. interferometric setup. Using hydrodynamical simulations post-processed with radiative transfer, we construct mock data cubes for two models of helium reionization. We compute the power spectrum of the signal and forecast the signal-to-noise ratio for SKA-1 MID, DSA-2000, and a PUMA-like survey, in both observational setups. The two scenarios produce distinct power spectra, but the faintness of the signal, largely caused by weak coupling between the spin temperature and the kinetic temperature in low-density regions of the IGM, combined with high instrumental noise, makes detection very difficult within realistic integration times for current surveys. A PUMA-like survey operating in single-dish mode could, however, detect the 3.46 cm signal with an integrated signal-to-noise ratio of a few in < 1000 h in both scenarios. Distinguishing helium reionization scenarios with 3.46 cm line-intensity mapping therefore remains challenging for current facilities. Our results, however, indicate that next-generation, high-sensitivity surveys with optimized observing strategies, especially when combined with complementary probes of the IGM, could begin to place meaningful constraints on the timing and morphology of helium reionization.

Primordial dark matter halos are well understood from cold dark matter-only simulations. Since they can contract significantly as baryons settle into their centers, direct comparisons with observed galaxies are complicated. We present an approach to reversing the halo contraction by numerically calculating the halo response to baryonic infall and iterating the initial condition. This allowed us to derive spherically averaged primordial dark matter halos for observed galaxies. We applied this approach to the Milky Way and found that the latest Gaia measurements for the rotation velocities imply an odd primordial Galactic halo: Its concentration and total mass differ by more than 3$\sigma$ from the predictions, and the density profile presents an inner core that is too shallow and an outer decline that is too steep to be compatible with the cold dark matter paradigm.

Very long baseline interferometry (VLBI) achieves the highest angular resolution in astronomy. VLBI measures corrupted Fourier components, known as visibilities. Reconstructing on-sky images from these visibilities is a challenging inverse problem, particularly for sparse arrays such as the Event Horizon Telescope (EHT) and the Very Long Baseline Array (VLBA), where incomplete sampling and severe calibration errors introduce significant uncertainty in the image. To help guide convergence and control the uncertainty in image reconstructions, regularization on the space of images is utilized, such as enforcing smoothness or similarity to a fiducial image. Coupled with this regularization is the introduction of a new set of parameters that modulate its strength. We present a hierarchical Bayesian imaging approach (Hierarchical Interferometric Bayesian Imaging, HIBI) that enables the quantification of uncertainty for al parameters. Incorporating instrumental effects within HIBI is straightforward, allowing for simultaneous imaging and calibration of data. To showcase HIBI's effectiveness and flexibility, we build a simple imaging model based on Markov random fields and demonstrate how different physical components can be included, e.g., black hole shadow size, and their uncertainties can be inferred. For example, while the original EHT publications were unable to constrain the ring width of M87*, HIBI measures a width of $9.3\pm 1.3\,\mu{\rm as}$. We apply HIBI to image and calibrate EHT synthetic data, real EHT observations of M87*, and multifrequency observations of \oj287. Across these tests, HIBI accurately recovers a wide variety of image structures and quantifies their uncertainties. HIBI is publicly available in the this http URL VLBI software repository.

It is relatively rare for gas giant planets to have resonant or near-resonant companions, but these systems are particularly useful for constraining planet formation and migration models. In this study, we examine Kepler-1624b, a sub-Saturn orbiting an M dwarf that was previously found to exhibit transit timing variations with an amplitude of approximately 2 minutes, suggesting the presence of a nearby non-transiting companion. We reanalyze the transits from archival Kepler data and extend the TTV baseline by 11 years by combining TESS data with three new ground-based transit observations from Palomar and Las Cumbres Observatories. We jointly fit these datasets and find that the TTV amplitude is significantly weaker in our updated analysis. We calculate the Bayes factor for a one-planet versus two-planet model and find that the one-planet model is preferred. Our results highlight the need for careful analysis of systems with relatively low amplitude TTV signals that are identified in large automated catalogs.

Entity is a new-generation, fully open-source particle-in-cell (PIC) code developed to overcome key limitations in astrophysical plasma modeling, particularly the extreme separation of scales and the performance challenges associated with evolving, GPU-centric computing infrastructures. It achieves hardware-agnostic performance portability across various GPU and CPU architectures using the Kokkos library. Crucially, Entity maintains a high standard for usability, clarity, and customizability, offering a robust and easy-to-use framework for developing new algorithms and grid geometries, which allows extensive control without requiring edits to the core source code. This paper details the core general-coordinate special-relativistic module. Entity is the first PIC code designed to solve the Vlasov-Maxwell system in general coordinates, enabling a coordinate-agnostic framework that provides the foundational structure for straightforward extension to arbitrary coordinate geometries. The core methodology achieves numerical stability by solving particle equations of motion in the global orthonormal Cartesian basis, despite using generalized coordinates like Cartesian, axisymmetric spherical, and quasi-spherical grids. Charge conservation is ensured via a specialized current deposition technique using conformal currents. The code exhibits robust scalability and performance portability on major GPU platforms (AMD MI250X, NVIDIA A100, and Intel Max Series), with the 3D particle pusher and the current deposition operating efficiently at about 2 nanoseconds per particle per timestep. Functionality is validated through a comprehensive suite of standard Cartesian plasma tests and the accurate modeling of relativistic magnetospheres in curvilinear axisymmetric geometries.

We analyze the H2 emission observed in the HH46 Class I system as part of PROJECT-J (Protostellar Jets Cradle Tested with JWST), to investigate the origin and excitation of the warm molecular outflow. We used NIRSpec and MIRI spectral maps (1.6-27.9 microns) to trace the structure and physical conditions of the outflow. By fitting the H2 rotational diagrams with a multi-temperature gas model, we derived key physical parameters including temperature, extinction, column densities, and the ortho-to-para ratio. This information is combined with a detailed kinematical analysis and comparison with irradiated shock models. We find no evidence of H2 temperature or velocity stratification from the axis to the edge of the outflow, as would be expected in MHD disk-wind models and as observed in other outflows. Instead, the observations suggest that the H2 emission arises from shock interactions between jet bow shocks and/or wide-angle winds with the ambient medium and cavity walls. NIRSpec emission and velocity maps reveal expanding molecular shells, likely driven by the less luminous source in the binary system. We infer an accretion rate of less than 10^-9 solar masses per year for the secondary source, approximately one order of magnitude lower than that of the primary. The H2 emission is consistent with excitation by low-velocity (approximately 10 km/s) J-type shocks, irradiated by an external UV field that may originate from strong dissociative shocks driven by the atomic jet. Future JWST observations will further constrain the evolution of the expanding shell and the mechanisms driving the outflow.

The relationship between the composition of rocky exoplanets and their host stars is fundamental to understanding planetary formation and evolution. However, previous studies have been limited by inconsistent datasets, observational biases and methodological differences. This study investigates the compositional relationship between rocky exoplanets and their host stars, utilizing a self-consistent and homogeneous dataset of 21 exoplanets and their 20 host stars. By applying sophisticated interior structure modeling and comprehensive chemical analysis, we identify a potential 1:1 best-fit line between the iron-mass fraction of planets and their host stars equivalent with a slope of $m = 0.94^{+1.02}_{-1.07}$ and intercept of $c = -0.02^{+0.31}_{-0.29}$. This results are consistent at the 1$\sigma$ level with other homogeneous studies, but not with heterogeneous samples that suggest much steeper best-fit lines. Although, our results remain tentative due to sample size and data uncertainties, the updated dataset significantly reduces the number of super-Mercuries from four to one, but it remains that several high-density planets are beyond what a primordial origin would suggest. The planets in our sample have a wider range of compositions compared to stellar equivalent values, that could indicate formation pathways away from primordial or be the result of random scattering owing to current mass-radius uncertainties as we recover the observed outliers in mock population analysis $\sim15\%$ of the time. To truly determine whether the origin is primordial with a 1:1 true relation, we find that sample of at least 150 planets is needed and that stars that are iron enrich or depleted are high value targets.

Strong gravitational lensing is a powerful probe for studying the fundamental properties of dark matter on sub-galactic scales. Detailed analyses of galaxy-scale lenses have revealed localized gravitational perturbations beyond the smooth mass distribution of the main lens galaxy, largely attributed to dark matter subhalos and intervening line-of-sight halos. Recent studies suggest that, in contrast to subhalos, line-of-sight halos imprint distinct anisotropic features on the two-point correlation function of the effective lensing deflection field. These anisotropies are particularly sensitive to the collisional nature of dark matter, offering a potential means to test alternatives to the cold dark matter paradigm. In this study, we explore whether a neural density estimator can directly identify such anisotropic signatures from galaxy-galaxy strong lens images. We model the multipoles of the two-point function using a power-law parameterization and train a neural density estimator to predict the corresponding posterior distribution of lensing parameters, alongside parameter distributions for dark matter substructure. Our results show that recovering the dark matter substructure mass functions and mass-concentration parameters remains challenging, owing to difficulties in generating uniform training data set while using physically motivated priors. We also unveil an important degeneracy between the line-of-sight halo mass-function amplitude and the subhalo mass-function normalization. Furthermore, the network exhibits limited accuracy in predicting the two-point function multipole parameters, suggesting that both the training data and the adopted power-law fitting function may inadequately represent the true underlying structure of the anisotropic signal.

Low surface brightness galaxies (LSBs) are an important class of galaxies that allow us to broaden our understanding of galaxy formation and test various cosmological models. We present a survey of low surface brightness galaxies at $0.4 < z_{\rm phot} < 0.8$ in the GOODS-S field using JADES data. We model LSB surface brightness profiles, identifying those with $\bar{\mu}_{\rm eff} > 24$ mag arcsec$^{-2}$ in the F200W JWST/NIRCam filter. We study the spatial distribution, number density, Sérsic profile parameters, and rest-frame colours of these LSBs. We compare the photometrically-derived star formation histories, mass-weighted ages, and dust attenuations of these galaxies with a high surface brightness (HSB) sample at similar redshift and a lower redshift ($z_{\rm phot} < 0.4$) LSB sample, all of which have stellar masses $\lesssim 10^8 M_{\odot}$. We find that both the high and the low redshift LSB samples have low star formation (SFR$_{100} \lesssim 0.01$ $M_{\odot}$ yr$^{-1}$) compared with the HSB sample (SFR$_{100} \gtrsim 0.01$ $M_{\odot}$ yr$^{-1}$). The star formation histories show that the LSBs and HSBs possibly come from the same progenitors at $z \gtrsim 2$, though the histories are not well constrained for the LSB samples. The LSBs appear to have minimal dust, with most of our LSB samples showing $A_V < 1$ mag. JWST has pushed our understanding of LSBs beyond the local Universe.

We report a series of images of the Tycho supernova remnant at eight epochs extending over thirty years: 1986-2016. In addition to our H{\alpha} images, we have obtained matched continuum images which we subtract to reveal faint emission, including a far more extensive network of optical knots and filaments than reported previously. The deepest images also show an extremely faint, fairly diffuse arc of emission surrounding much of the circumference of Tycho to the southeast and south, coinciding with the rim of the radio/X-ray shell. We have measured proper motions for 46 filaments, including many fainter ones near the Tycho outer rim. Our measurements are generally consistent with previous ones by Kamper and vandenBergh (1978), but ours have far greater precision. Most optical filaments at the shell rim have expansion indices reasonably consistent with the Sedov value (0.40), while the interior filaments have somewhat smaller values, as expected. From the combination of proper motions of filaments at the shell rim and shock velocity at the same positions, one should be able to calculate the distance to Tycho by simple geometry. Determination of the shock velocity from broad Balmer-line profiles is subject to model uncertainties, but the availability of dozens of such filaments with a range of conditions offers the possibility to substantially improve the distance determination for Tycho.

There is clear observational evidence that the main Class 0/I stages of the star formation process are associated with powerful collimated outflows (jets), which sometimes propagate up to distances as large as $10^{4-5}$ au scales in molecular clouds. Additionally, intermediate high-mass and low-mass protostars have often been observed to form in crowded clusters, where the typical separation distance between any two cluster members is of the same order or smaller than the scale of the outflow length. Therefore, there must be an interaction between the molecular outflows of different protostars within the protostellar association. A good example of this is the case of Cepheus E-mm, which is a protostellar outflow extending over a few dozen au. At its core is a binary system consisting of two protostars, Cep E-A and Cep E-B, separated by about 1000 au. Both protostars eject molecular jets at velocities of ~100 km/s. The interaction between these molecular outflows provides an opportunity to study the effects of jet collisions in a clustered star-forming environment, as they may leave detectable imprints on the morphology of the main envelope of the system. Our work aims to study the effects of the collision of molecular jets associated with the components of the binary system Cep-A and Cep-E, analyzing the disruption or reduction of molecular emission in the main envelope of the system, which the molecular outflow { launched} by Cep-A presumably pushes. If we characterize the collision in this system, we can provide insights into the expected morphology and molecular emissions in collisions of molecular outflows { associated to star forming process.

The High Time Resolution Universe (HTRU) survey is an all-sky survey looking for pulsars and other radio transients. A new single-pulse (SP) search pipeline is presented, tailored to the northern part of the HTRU survey collected with the 100m Effelsberg Radio Telescope. In a selection of this data, synthetic SPs are injected with frequency-time structures resembling those of the detected Fast Radio Burst (FRB) population and processed by the pipeline to characterize its performance. Therefore, several new software toolkits have been developed (FRBfaker and RFIbye) to enable the injection of SPs with complex frequency-time structures and cope with the Radio Frequency Interference (RFI) in the survey's data. The operation of these toolkits is described alongside the overall functionality of the SP pipeline. Qualification of the pipeline confirmed that it is ready to process all the HTRU-North data. Additionally, the survey's sensitivity to SPs, the impact of RFI thereon, the performance of the deep-learning classifier FETCH, and some insights that may be used to improve the pipeline's performance in the future are determined. Within the small data sample analysed, 21 known pulsars and a RRAT are detected. In addition, eight faint SP trains that might originate from yet undiscovered neutron stars and 141 isolated SP candidates were discovered.

We present X-ray observations of the periodic optical source ZTF J185139.81+171430.3 (hereafter ZTF J1851) by the XMM, NICER and NuSTAR telescopes. The source was initially speculated to be a white dwarf (WD) pulsar system due to its short period ($P\sim12$ min) and highly-modulated optical lightcurves. Our observations revealed a variable X-ray counterpart extending up to 40 keV with an X-ray luminosity of $L_X \sim 3\times10^{33}$ erg s$^{-1}$ (0.3--40 keV). Utilizing timing data from XMM and NICER, we detected a periodic signal at $P_{\rm spin}=12.2640(7)\pm0.0583$ min with $>6\sigma$ significance. The pulsed profile displays $\sim 25\%$ and $\sim10\%$ modulation in the 0.3--2 and 2--10 keV bands, respectively. Broadband X-ray spectra are best characterized by an absorbed optically-thin thermal plasma model with $kT \approx 25$ keV and a Fe K-$\alpha$ fluorescent line at 6.4 keV. The bright and hard X-ray emission rules out the possibility of a WD pulsar or ultra-compact X-ray binary. The high plasma temperature and Fe emission lines suggest that ZTF J1851 is an intermediate polar spinning at 12.264 min. We employed an X-ray spectral model composed of the accretion column emission and X-ray reflection to fit the broadband X-ray spectra. Assuming spin equilibrium between the WD and the inner accretion disk, we derived a WD mass range of $M_{\rm WD}=(1.07\rm{-}1.32)M_{\odot}$ exceeding the mean WD mass of IPs ($\langle M_{\rm WD} \rangle = 0.8 M_\odot)$. Our findings illustrate that follow-up broadband X-ray observations could provide unique diagnostics to elucidate the nature of periodic optical sources anticipated to be detected in the upcoming Rubin all-sky optical surveys.

Following the detection of a radio-loud Active Galactic Nucleus (AGN) in the z=4.3 protocluster SPT2349-56, we have obtained additional observations with MeerKAT in S-band (2.4 GHz) with the aim of further characterizing radio emission from amongst the ~30 submillimeter (submm) galaxies (SMGs) identified in the structure. We newly identify three of the protocluster SMGs individually at 2.4GHz as having a radio-excess, two of which are now known to be X-ray luminous AGN. Two additional members are also detected with radio emission consistent with their star formation rate (SFR). Archival MeerKAT UHF (816 MHz) observations further constrain luminosities and radio spectral indices of these five galaxies. The Australia Telescope Compact Array (ATCA) is used to detect and resolve the central two sources at 5.5 and 9.0 GHz finding elongated, jet-like morphologies. The excess radio luminosities range from L1.4,rest = (1-20)x10^25 W/Hz, ~10-100x higher than expected from the SFRs, assuming the usual far-infrared-radio correlation. Of the known cluster members, only the SMG `N1' shows signs of AGN in any other diagnostics, namely a large and compact excess in CO(11-10) line emission. We compare these results to field samples of radio sources and SMGs. The overdensity of radio-loud AGN in the compact core region of the cluster may be providing significant heating to the recently discovered nascent intra-cluster medium (ICM) in SPT2349-56.

TIC 322208686 is known to be a detached binary that exhibits two types of variability: pulsation and eclipse. We present the physical properties of the target star using the short-cadence TESS data from sectors 24, 57, and 58, and our echelle spectra that show the presence of a tertiary companion. The spectral analysis led to the triple-lined radial velocities and the atmospheric parameters of the eclipsing components. Joint modeling of these observations reveals that the eclipsing pair contains two F-type stars with masses $1.564\pm0.012$ $M_\odot$ and $1.483\pm0.012$ $M_\odot$, radii $1.588\pm0.011$ $R_\odot$ and $1.500\pm0.012$ $R_\odot$, effective temperatures $7028\pm100$ K and $7020\pm110$ K, and luminosities $5.51\pm0.32$ $L_\odot$ and $4.90\pm0.32$ $L_\odot$. The light contributions of the three stars obtained from this modeling match well with those calculated from the observed spectra. The binary star parameters are in satisfactory agreement with evolutionary model predictions for age $t$ = 0.4 Gyr and metallicity $Z$ = 0.03. We extracted 11 significant frequencies from the TESS light residuals with the binary effects removed. Of these, five signals between 0.65 day$^{-1}$ and 1.89 day$^{-1}$ can be considered as $\gamma$ Dor pulsations originating mainly from the primary component, while the other frequencies are likely instrumental artifacts or combination terms. These results suggest that TIC 322208686 is a hierarchical triple, containing a pulsating eclipsing pair and a tertiary companion.

It is usually assumed that mutual synchronisation of a tidal two-body system happens through tidal recession, assuming the reduced Hill sphere is not reached. However, synchronisation can be achieved also via tidal approach, provided the Roche limit is not crossed. For each of the two scenarios, hereafter referred to as Scenario 1 and Scenario 2, respectively, we derive the condition under which the evolving synchronicity radius catches up with the tidally evolving orbit. We consider these two scenarios for the Pluto-Charon system, examine the impact origin hypothesis of Charon's formation, and propose that capture is a likelier option. We investigate Scenario 2, both analytically and numerically, where the orbital evolution of Charon starts at a higher altitude than present and undergoes tidal descent. In Scenario 2, the greater initial orbital separation between the partners reduces tidally induced thermal processes and fracturing, as compared to Scenario 1. In several study cases, we also observe temporary locking of Charon into higher spin-orbit resonances (3:2 to 7:2) in the first 0.5Myr of the system's evolution.

We forecast the detectability of the Doppler magnification dipole with a joint analysis of galaxy spectroscopic redshifts and size measurements. The Doppler magnification arises from an apparent size variation caused by galaxies' peculiar velocities when mapping them from redshift space to real space. This phenomenon is the dominant contribution to the convergence at low redshifts ($\lesssim$ 0.5). A practical observational strategy is to cross-correlate a galaxy number count tracer, e.g. from the Dark Energy Spectroscopic Instrument (DESI) Bright Galaxy Survey, with the convergence field reconstructed from galaxy size measurements obtained by the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST). To assess the achievable precision of galaxy size measurements, we simulate LSST Y1-quality galaxy images with Galsim and measure them with the Galight profile fitting package. Our investigations, based on galaxy populations from LSST's synthetic galaxy catalogue cosmoDC2, show that the variance due to intrinsic galaxy size variation dominates over size measurement errors as expected, but may be lower than previous studies have suggested. Under our analysis assumptions, the Doppler magnification dipole would be detectable with a signal-to-noise ratio $\geq 10$ in multiple redshift bins between $0.1 \leq z \leq 0.5$ with DESI spectroscopic redshifts and LSST imaging.

Earth-like exoplanets can create resonant structures in exozodiacal dust through mean motion resonances (MMRs). These structures not only suggest the presence of such planets, but also act as potential noise sources in future mid-infrared (MIR) nulling interferometry observations. We aim to investigate how resonant structures in exozodiacal dust vary across stellar spectral types (F4--M4), and to evaluate how stellar wind drag affects their morphology and brightness in mature planetary systems. We conducted numerical simulations of dust dynamics, extending earlier studies by including spectral type variation in stellar wind drag in addition to Poynting-Robertson (PR) drag. Our models represented systems of a few Gyr hosting an Earth-like exoplanet in the habitable zone (HZ). We produced spatially resolved maps of optical depth and thermal emission for different stellar spectral types. Our simulations showed that resonant ring structures were formed for all stellar spectral types considered. In particular, we found that stellar wind drag played a critical role in shaping dust dynamics around old M-type stars, where it could dominate over PR drag by a factor of approximately 44. This reduced the contrast of resonant rings relative to the background disk, compared to cases without spectral type variation in stellar wind. Across different spectral types, the optical depth contrast of the resonant ring increased for lower-mass stars, assuming a fixed background level. Asymmetric thermal emission distributions were derived across all spectral types, which peaked for K-type stars. Our findings highlight the importance of incorporating both resonant dynamics and stellar wind effects when modeling exozodiacal dust around stars of different spectral types.

We announce a recently detected outburst that is currently only a few months old, and probably of FU Orionis type. The progenitor to the outburst was an emission-line, flat-spectrum SED young stellar object located in the W5 region, though somewhat outside the main star formation action. We present optical, near-infrared, and mid-infrared lightcurves that illustrate the quiescent state of [KAG2008] 13656 and its subsequent$\Delta r \approx -4$ mag and $\Delta J\approx -3$ mag outburst over $\sim$75 days in late-2025. Follow-up optical and near-infrared spectroscopy confirms the expected features from an FU Ori disk and outflow.

PSR J2039-5617 is a redback millisecond pulsar binary system consisting of a compact star with a mass of 1.1-1.6 $M_\odot$ and a low-mass companion of 0.15-0.22 $M_\odot$. For this binary, we performed a timing analysis using 16 years of data from the Fermi Large Area Telescope, covering the period from 2008 August to 2024 October. Our analysis detected an orbital modulation with a period of 0.2279781 days at a significance level of $\sim4\sigma$, which is in good agreement with previous findings. However, unlike previous reports, we identified a transition in the orbital modulation around 2021 August, after which the orbital signal disappeared. We speculate that the system may be undergoing a transition from a rotation-powered to an accretion-powered state at this epoch. Additionally, we conducted the phase-resolved and spectral analyses, and in the phase-resolved results, we observed an anti-correlation between its $\gamma$-ray and X-ray emissions, which consistent with the predictions of high-energy radiation models for such systems. We provide some predictive discussions based on the results of $\gamma$-ray data analysis, and future Fermi-LAT observations will determine whether these predictions hold true.

We report a constraint on the cosmological variation of the proton g-factor, $g_p$. By comparing the measured redshifts between \mbox{H\,{\sc i}} 21 cm and OH 18 cm lines observed with the newly commissioned Five-hundred-meter Aperture Spherical radio Telescope (FAST) toward PKS 1413+135 at $z$ = 0.24671, we obtain $\Delta g_{p}/g_{p} = (-4.3\pm2.5)\times10^{-5}$, which is more than two orders of magnitude more sensitive than previous constraints. In addition, we obtain sensitive constraints of $\Delta (\mu\alpha^{2})/(\mu\alpha^{2}) = (2.0\pm1.2)\times10^{-5}$ and $\Delta (\mu\alpha^{2}g_{p}^{0.64})/(\mu\alpha^{2}g_{p}^{0.64}) = (-4.7\pm1.9)\times10^{-6}$.

AM CVn binaries are the most compact of accreting binaries having orbital periods in the range ~5-70 min. They consist of a white dwarf accreting hydrogen deficient material from a degenerate or semi-degenerate star and are predicted to be amongst the verification sources for future gravitational wave observatories such as LISA. Using the recent catalogue of Green et al (2025) I focus attention on the orbital period range in which outbursts are seen from AM CVn's. I examine in more detail the outburst properties of KL Dra which has an outburst every few months and has many sectors of TESS data as an open resource. Using observational data on the outbursting systems in general, I compare the outburst recurrence time, duration and amplitude as a function of orbital period with the predictions of the disc instability model. The recurrence time is well described, although there is some evidence that the amount of material in the disc at the end of the quiescence phase is less than earlier model assumptions. The distribution of the outburst duration appears to be dependent on the cadence of the observations and how it is defined. Similarly the amplitude distribution is dependent on cadence and the filter, which causes an apparent spread in distribution. Both of these features need to be systematically studied using consistent benchmarks. AM CVn binaries remain an excellent sources to test models which aim to predict the properties of disc accreting systems.

The Galactic Center (GC) is an extreme region of the Milky Way that is host to a complex set of thermal and non-thermal structures. In particular, the GC contains high-density gas and dust that is collectively referred to as the Central Molecular Zone (CMZ). In this work, we study a subset of HNCO filaments identified in band 3 ALMA observations of the GC obtained by the ALMA CMZ Exploration Survey (ACES) that are comparable to high density filaments identified in the Galactic Disk. We compare the orientation of the magnetic field derived from 214 um SOFIA and 850 um JCMT observations with the filament orientation to determine which mechanisms dominate the formation of these filaments. We observe a large range of magnetic orientations in our observed filaments indicating the complex environments the filaments are located in. We also compare the observational results to synthetic data sets created using an MHD model of the GC. Our analysis reveals that the dominant mechanisms local to the HNCO filaments vary throughout the GC with some filaments being dominated by supersonic turbulence and others by subsonic turbulence. The comparison to synthetic observations indicates that the observed filaments are in magnetically dominated environments that could be supporting these filaments against collapse. Our results on the CMZ filaments are also compared to results obtained on similar filaments located in the Galactic Disk, and we find that the filaments studied here are possible CMZ analogs to the dense filamentary "bones" observed previously in the Galactic Disk.

The atlas of the spectrum of the B[e] star CI Cam obtained with the 6-meter BTA telescope in combination with the NES spectrograph in the wavelength interval of 395$÷$780 nm with a resolution of $\lambda/\Delta\lambda\ge60\,000$ is presented. The atlas contains about 400 spectral features and illustrates the diversity of spectral features of the unique star, forming in a complex circumstellar environment. The spectrum is dominated by forceful single-peaked HI, HeI emissions and numerous double-peaked permitted and forbidden emissions of ions of chemical elements starting from the CNO triad up to metals (Mg, Al, Ti, V, Cr, Fe) with practically ``rectangular'' profiles. The FeII and [FeII] emissions predominate in the spectrum. However, several other double-peaked forbidden ion emissions were also detected: [VII], [Cr II], [Ni II]. The atlas is presented graphically, with a separate figure corresponding to each echelle order. A list of identified lines, including a number of known interstellar features is presented in table form. CI Cam supergiant status has been confirmed.

Recent studies on the dark photon (DP) production in collapsing stars argue that the cooling effect induced by DPs can hinder supernova explosions and lead to a ``failing supernova" constraint on the photon-DP mixing parameter $\epsilon$. In order to verify the idea, we perform two-dimensional neutrino-radiation hydrodynamic simulations coupled with the DP production with the mass of 0.3 MeV. We find that the shock revival does not happen until the end of the simulations when $\epsilon\gtrsim3\times10^{-9}$. The photon-DP mixing parameter above this value can be excluded by the failing supernova argument. Interestingly, our constraint is close to the one reported by the previous studies which adopted the post-processing framework. This result motivates one to investigate a wider parameter range of DPs with self-consistent simulations and evaluate uncertainties in the constraint.

The calibration of tomographic redshift distributions is essential for cosmological analysis of weak lensing data. In this work, we calibrate all four tomographic bins of the Hyper Suprime Camera (HSC) weak lensing catalog with the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 and 2 using the clustering redshifts technique. We include $z>1.2$ redshift sources such as emission line galaxies (ELG) and quasars (QSO) sources in our calibration, which were not available in the previous HSC calibration (Rau et al. 2022, arXiv:2211.16516), allowing a complete calibration of all the redshift bins. We find the first tomographic bin exhibits a small shift towards low redshifts. The second bin is in good agreement with the photometric calibration, while third and fourth bin exhibit a shift towards higher redshifts. However, these shifts are considerably smaller than the shifts obtained in the HSC Year 3 cosmic shear analyses. We evaluate the impact of galaxy bias and magnification effects from all the samples on the measurements, finding them to be small, and we propose corrections to reduce them further. We model the redshift distributions with splines and compare our results to previous analyses as well as to other parameterizations found in literature. For the two high-redshift tomographic bins, we find the shifts to higher redshifts with respect to the measurements performed in Rau+2022 to be $\Delta z_3=-0.039^{+0.020}_{-0.021}$ and $\Delta z_4=-0.048^{+0.012}_{-0.012}$.

We reanalyze cosmological constraints from Hyper Suprime-Cam (HSC) Y3 shear-shear correlation function using new calibration of the tomographic redshift distribution via the clustering redshifts method with DESI spectroscopy presented in Choppin de Janvry et al. (2025a). We present both importance sampling of the original MCMC chains by HSC, applying the weights of our newly calibrated $\Delta z$ priors, as well as full MCMC analysis with new photometric redshift distributions, finding consistent results between the two. We obtain the growth of structure parameter $S_8\equiv\sigma_8\sqrt{\Omega_m/0.3}=0.805\pm{0.018}$, compared to previous HSC Y3 result of $S_8=0.769^{+0.031}_{-0.034}$, which is a 1.8 reduction of error due to the improved clustering redshift calibrations, with the central value shifting considerably higher towards Planck cosmology. With the new photometric redshift calibration, HSC Y3 has comparable constraining power to the recent KIDS Legacy and DES Y3. Combining all three gives $S_8=0.813^{+0.009}_{-0.010}$, which can be compared to $S_8=0.828\pm 0.012$ from CMB lensing. Overall there is no evidence for deviation from Planck on $S_8$ in any of the weak lensing analyses, and combining galaxy lensing with CMB lensing gives a sub-percent constraint $S_8=0.818\pm 0.007$, comparable in both precision and value to the most recent Planck+ACT+DESI constraint $S_8=0.812\pm 0.007$.

The highest-energy cosmic neutrino detected by the ARCA detector of KM3NeT has reignited the quest to pinpoint the sources of ultrahigh-energy cosmic rays (UHECRs; $E\gtrsim 0.1$ EeV). By uncovering the associated multimessenger signals, we investigate the origin of the 220 PeV $\nu_\mu$ event KM3-230213A from a transient source that accelerated cosmic rays to $\sim 10$ EeV. UHECR protons that escape the source interact with the cosmic background radiation, producing a PeV-EeV cosmogenic neutrino spectrum. The secondary $e^\pm$ and $\gamma$-rays initiate an electromagnetic cascade, resulting in a cosmogenic $\gamma$-ray spectrum that peaks at a delayed time due to deflection of charged particles in the extragalactic magnetic field (EGMF). Our results shed light on the nature of the UHECR source for the $\nu_\mu$ event and provide crucial insights into the detection of multi-TeV $\gamma$-rays of cosmogenic origin from similar past cosmological transients. Using the $\gamma$-ray sensitivity of currently operating and next-generation imaging atmospheric Cherenkov telescopes, the flux and time-delay distribution can constrain the source distance. We further show that the detection of such a $\gamma$-ray signal above the background depends on the EGMF strength. Together with the non-detection of coincident spatial or temporal photon counterparts at the current epoch, this detection is the first compelling candidate for a sub-EeV cosmogenic neutrino.

JWST has already observed near-infrared transmission spectra of over a dozen super-Earths and sub-Neptunes. While some observations have allowed astronomers to characterize sub-Neptunes in unprecedented detail, small feature amplitudes and poorly-understood systematics have led to ambiguous results for others. Using the first seven targets from the COMPASS program, which is currently surveying 12 small planet atmospheres using NIRSpec/G395H, we investigate these timeseries systematics. We implement a model that uses the principle components of the normalized pixel fluxes to account for variations in the shape and position of the spectral trace. We find that observations with a smaller number of groups-per-integration benefit most profoundly from the use of this model, and that systematics are particularly strong between 2.8 and 3.5 $\mu$m. Despite these systematics, \texttt{pandexo} is a relatively accurate predictor of the precision of the spectra, with real error bars on average 5\% larger in NRS1 and 12\% larger in NRS2 than predicted. We compute new limits on metallicity and opaque pressure level for each target and compare these to previous results from the COMPASS program. Next, we co-add spectra from multiple targets to reduce the effective noise in the combined spectra in hopes of detecting transmission features in common between the targets, but this exercise does not yield compelling evidence any signals. We find that a handful of additional transits are sufficient to break the degeneracy between metallicity and aerosols for the majority of our targets, pointing towards the possibility of unraveling the mysteries of these worlds with future allocations of JWST time.

The inner stellar mass surface density within 1 kpc, Sigma1, has emerged as a suitable proxy for bulge growth and galaxy quenching. However, the dependence of cold gas content on Sigma1 has not been thoroughly explored. In this paper, we examine the relationship between Sigma1, as well as the mass-relative parameter delSigma1, and the atomic (fHI) and molecular (fH2) cold gas fractions in massive, nearby galaxies. We utilize a sample of 341 galaxies with HI data and 201 galaxies with H2 data from the xGASS and xCOLDGASS surveys, spanning 0.02 <= z <= 0.05 and a stellar mass range of 10^10 <= M_*/M_odot <= 10^11.5. While we observe that a decline in both fHI and fH2 is associated with increasing Sigma1, we find that fH2 shows a sharper decline above a threshold value of delSigma1 = 0. In addition, the fraction of galaxies with AGN activity (Seyferts and LINERs) increases with delSigma1, with the greatest increase occurring between 0 <= delSigma1 <= 0.2 dex. We propose an evolutionary track in the plane of fH2-delSigma1, whereby molecular gas depletion at fixed mass coincides with a rise in AGN activity. Our results suggest that central bulge growth is more tightly coupled to the depletion of molecular gas rather than atomic gas, with AGN feedback possibly contributing to this process. Our work highlights the utility of Sigma1 and delSigma1 as tracers of quenching in massive galaxies.

We present a comprehensive analysis of 475 ks (438 ks unpublished & 37 ks archival) XMM-Newton/EPIC-pn observation of a nearby, highly inclined, star-forming, luminous infrared galaxy NGC 3221 through spatial, temporal, and spectral information. We confirm the presence of a low-luminosity (presumably Compton-thick) AGN. The 0.4$-$12 keV luminosity and the hardness ratio of the six ultra-luminous X-ray sources (ULX) previously identified in Chandra data exhibit diverse variability on day-scale. The collective emission from unresolved sources exhibits a different day-scale variability. We have also discovered two new predominantly soft ($<1$ keV) sources. One of these has an enigmatic spectral shape featuring a soft component, which we interpret as a superbubble in NGC 3221, and a variable hard component from a compact object, which is unresolved from the superbubble. We do not confidently detect any X-ray emission from SN 1961L. The hot gas in the ISM (out to $\pm$6 kpc from the disk plane) and the extraplanar region (6$-$12 kpc) both require two thermal phases at $\sim 0.15$ keV and $\sim 0.55$ keV. The $\sim 0.55$ keV component is fainter in the ISM than the $\sim 0.15$ keV component, but the emission from the latter falls off more steeply with disk height than the former. This makes the extraplanar region hotter and less dense than the ISM. The proximity of NGC 3221 and the occurrence of the underluminous AGN offer a unique observing opportunity to study the hot diffuse medium in conjunction with nuclear and disk-wide point sources.

In this work we investigate curvature perturbations and non-Gaussianity arising from Higgs modulated reheating in the early Universe. We employ three different methods -- the period-averaging (PA) method, the exact method, and the non-perturbative $\delta N$ formalism -- to compute the power spectrum and bispectrum of curvature perturbations. Our results show that the non-perturbative $\delta N$ method provides a reliable estimate across a wide range of reheating time and Higgs field values, including regimes where the Higgs field oscillates significantly after inflation. We find that a smaller Higgs self-coupling ($\lambda$) leads to a larger curvature perturbation, with the non-Gaussianity predominantly taking a local shape. This highlights the importance of considering non-perturbative effects in calculating the curvature perturbation during Higgs modulated reheating, especially for smaller values of $\lambda$. Our findings offer valuable insights into the dynamics of reheating and the generation of primordial perturbations in the early Universe.

We present a comprehensive re-evaluation of pulsar radio spectra using the largest curated dataset of calibrated flux densities to date, comprising 897 pulsars, and employing a robust Bayesian framework for model comparison alongside frequentist methods. Contrary to the established consensus that pulsar spectra are predominantly simple power laws, our analysis reveals that complex spectral shapes with curvature or breaks are in fact the norm. The broken power law emerges as the most common spectral shape, accounting for 60\% of pulsars, while the simple power law describes only 13.5\%, with 68.8\% of pulsars decisively favoring curved or broken models. We further identify 74 confident gigahertz-peaked spectrum pulsars, and demonstrate that millisecond pulsars frequently exhibit spectral curvature. A key finding is that the previously reported dominance of the simple power law was largely a statistical artifact of the frequentist method used in earlier work. These findings substantially revise the prevailing view of pulsar spectra and establish a critical, model-classified foundation for future theoretical work.

The Space-based multiband astronomical Variable Objects Monitor (SVOM), detected its first short gamma-ray burst, GRB 240715A, in-flight, which was jointly observed by Fermi. Based on observational data of SVOM/GRM and Fermi/GBM, we perform a comprehensive temporal and spectral analysis for individual pulse in the prompt emission of this burst, and novel characteristics are revealed. Firstly, opposite evolutions of spectral lag are found in the first and third pulse of this burst. Second, the large negative lag of the first pulse is an outlier in short GRB sample, especially when the pulse duration is considered. Spectral analysis shows that the negative lag of the first pulse is caused by the evolution of spectrum index, and is irrelevant to Epeak, which is inconsistent with the previous study. The intrinsic mechanism is probably attributed to electron cooling in the decaying magnetic field, which leads to the continuous hardening of the spectrum index and results in negative lag. Furthermore, spectral analysis also shows that the third pulse is more likely to be described by a quasi-thermal spectrum, indicating the existence of photospheric emission. It is difficult to explain how the synchrotron radiation appears before photospheric emission in a single GRB and some assumptions are discussed.

By investigating the time-variable 230 GHz images using ray-tracing general relativistic radiative transfer calculation, we propose a novel method for estimating the spin parameter of the supermassive black hole at the M87 center by utilizing the sudden and short-term increase in emissivity in the innermost region of the accretion disk. It is found that the flux of the photon ring increases simultaneously as the flux of the direct ring, which brightens first, decreases, and then gradually diminishes, when the increase in emissivity persists for $15 t_{\rm g}$ with $t_{\rm g}$ being the light crossing time of the gravitational radius. The direct ring is formed by photons emitted from the vicinity of the innermost region of the disk and traveling directly to the observer without orbiting around the black hole, while the photon ring is formed by photons passing near the spherical photon orbit. The time-averaged width of the dark region between the direct ring and the photon ring (dark crescent) becomes thinner when the spin parameter is higher and the increase in the emissivity of the accretion disk is greater. The time variation of two rings also causes the intensity-weighted center to oscillate both in the direction of the black hole's angular momentum vector projected onto the screen ($Y$-direction) and in the perpendicular direction ($X$-direction). The amplitude of oscillatory time variation in the $X$-direction becomes large when the spin parameter is higher, and that in the $Y$-direction becomes large when the increase in the emissivity of the disk is large. The spin parameter can be estimated by combining the time-averaged dark crescent width and the ratio of the amplitudes in the $X$- and $Y$-directions. This method is applicable when the duration of the increase in emissivity of the accretion disk ranges at least from approximately 10-20 $t_{\rm g}$.

The BAO characteristic scale is a useful tool for understanding the evolution of the universe, especially the influence of dark energy on this evolution. In this work, we study the projection effect in transverse BAO, namely the mixing of BAO signals from different epochs caused by $z_c$ uncertainty within a chosen bin. We focus our forecast on two surveys of interest: the Square Kilometre Array (SKA) HI galaxy redshift survey and the Dark Energy Spectroscopic Instrument (DESI)'s final Luminous Red Galaxy (LRG) sample. We test the sensitivity in finding the transverse BAO depending of three bin configurations: a Gaussian, a top-hat and an intermediate of them semi-Gaussian. We also analyse the precision the bin widths $\sigma_z$ from smaller to wider bins $0.01<\sigma_z<0.1$. In order to correct these deviations, we propose a correction based on adjacent redshift to $z_c$, this would provide a semi-statistical correction instead of only relying on fiducial cosmology. Finally, we conclude that despite the higher shot-noise than the top-hat bin separation, the semi-Gaussian bin is the most accurate case to find the BAO signal and to constrain parameters through the angular power spectrum. A Gaussian binning gives the least precise parameter constraints compared to the other two cases.

Many circumstellar dust scattering regions have been detected and investigated with polarimetric imaging. However, the quantitative determination of the intrinsic polarization and of dust properties is difficult because of complex observational effects. This work investigates instrumental convolution and polarimetric calibration effect for high contrast imaging polarimetry with the aim to define procedures for accurate measurements of the circumstellar polarization. For this we simulate the instrumental convolution and polarimetric cancellation effects for a Gaussian PSF and an extended PSF_{AO} typical for a modern adaptive optics system. Further, polarimetric zero-point corrections (zp-corrections) are simulated for different cases like coronagraphic observations or systems with barely resolved circumstellar scattering regions. We find that the PSF convolution reduces the integrated azimuthal polarization Q_phi for the scattering region while the net Stokes signals Q and U are not changed. For non-axisymmetric systems a spurious U_phi-signal is introduced. These effects are strong for compact systems but scattering regions can still be detected down to small separations while unresolved scattering regions can be constrained by the central Stokes Q,U signal. The smearing by PSF_{AO} produces an extended, low surface brightness polarization signal changing the angular distribution of the polarization, but the initial signal can be recovered partly from the Stokes Q and U quadrant pattern. A polarimetric zp-correction applied for the removal of offsets from instrumental or interstellar polarization depends on the selected reference region and can also introduce strong bias effects for the azimuthal distribution of the polarization signal. Strategies for the zp-correction are described for coronagraphic data or observations of partly unresolved systems.

An approximate analytic one-dimensional model is constructed, for the accretion disk boundary layer surrounding a neutron star whose low magnetic field does not affect the process of accretion. A high luminosity model is considered, with radiation pressure dominant in the interior part of the disk.

We present a catalog of 3D Galactic HI filaments over 1/3 of the sky using Galactic Arecibo L-band Feed Array HI (GALFA-HI) data. The 3D filaments are defined to be linear HI features that are continuous in position-position-velocity (PPV) and are found with fil3d, an algorithm that expands on the 2D FilFinder. The catalog contains 3333 HI filaments between +/- 50 km/s at a range of Galactic positions. 1542 of the PPV filaments are identified as local at the distance of the wall of the Local Bubble, and 209 are likely at the disk-halo interface of our Galaxy. The catalog and properties of the PPV filaments are obtained after an unsharp mask (USM) is applied to the data. The widths of the filaments are consistently ~12' (0.34 pc at 100 pc), and constrained by the 4' resolution. The local filaments have median properties of N_HI = $6 \times 10^{18}$ cm$^{-2}$, M_HI = 0.17 M_sun, FWHM = 3.2 km/s, and length of 6.4 pc. The disk-halo population has similar column densities, but the median FWHM = 7.7 km/s, consistent with them being higher z-height, warmer structures. The L $\propto$ M$^{0.5}$ relationship found for the HI filaments and their bundling on the sky are consistent with a hierarchical structure, and is likely related to turbulence playing a role in their formation.

Accurate estimation of stellar parameters -- stellar age, lifetime, and evolutionary stage -- remains a fundamental challenge in astrophysics. We introduce a hybrid deep learning architecture combining multimodal spectroscopic and photometric data from SDSS DR17. The model comprises a Multi-Layer Perceptron for numerical features and a CNN with a Vision Transformer for spectra, with three output heads for age, lifetime, and evolutionary stage prediction. Training labels are derived from MIST v1.2 isochrones, with evolutionary stage binned into five classes (Hot, Medium, Cool, Subgiant, Red Giant). We conduct multi-phase evaluation: Phase I explores model architectures and data balancing strategies, Phase II tunes architectural complexity, and Phase III optimizes the multi-task loss composition. The final model achieves a balance between precision (Age RMSE 0.093 in $\log(\mathrm{yrs})$) and physical realism. Monte Carlo Dropout confirms well-calibrated uncertainties, enabling meaningful astrophysical interpretation and establishing a new benchmark for multimodal stellar parameter estimation.

Understanding dark energy requires precision measurements of the expansion history of the universe and the growth of large-scale structure. The 21 cm intensity mapping (21 cm IM) technique enables rapid large-area surveys that can deliver these measurements. China is constructing three hundred-meter-class single-dish radio telescopes, including the QiTai 110 m Radio Telescope (QTT), the 120 m Jingdong Radio Telescope (JRT), and the 120 m Huadian Radio Telescope (HRT), whose designs are well suited for 21 cm IM cosmology. We use a Fisher-to-MCMC forecasting framework to evaluate the baryon acoustic oscillations / redshift space distortions (BAO/RSD) measurement capabilities of QTT, JRT, and HRT and propagate them to dark-energy constraints in the $w_0w_a$CDM model. Our results show that achieving a redshift coverage up to $z_{\mathrm{max}} = 1$ is crucial for fully realising the potential of hundred-meter-class single-dish telescopes for 21 cm cosmology. If all three telescopes carry out 21 cm IM surveys over the same redshift range up to $z_{\mathrm{max}}=1$ and combine their BAO/RSD measurements, QTT+JRT+HRT yield $\sigma(w_0)=0.094$ and $\sigma(w_a)=0.487$, providing tighter constraints than DESI DR2 results.

Proximity of lensing critical curves features highly magnified portions of lensed galaxies. Accurate knowledge of the location and shape of the critical curve will be useful for understanding the nature of highly magnified stellar sources near critical curves and for revealing sub-galactic dark matter structures within the lens. In galaxy-cluster lenses, however, prediction of critical curves can be uncertain due to complexity in global mass modeling. Bracketing the critical curve with sparse lensed image pairs of compact continuum sources may not be sufficient as additional information available throughout the extended arc is not used. We explore and validate a kinematics-based method for locating the critical curve. This method leverages the continuous line-of-sight velocity profile of the lensed galaxy mapped through integral field spectroscopy of nebular emission lines, and combines an agnostic local lens model and a disk rotation model. Applying our method to a highly magnified region of the Dragon Arc in the Abell 370 cluster lensing field using archival VLT/MUSE IFU mapping of the H$\beta$ line, we constrain the critical curve to an uncertainty band with a half-width of 0.23" ($1\sigma$). This reveals locations of recently detected extremely magnified stars biased toward the negative-parity side of the critical curve, as predicted for intracluster microlensing. With future JWST/NIRSpec IFU mapping of the H$\alpha$ line at SNR $\simeq$ 10 (20), uncertainty may improve to 0.12" (0.08"). A measurement of this type with sufficiently small uncertainty may reveal small-scale wiggles in the shape of the critical curve where it crosses the giant arc, which can arise from the lensing perturbation of sub-galactic dark matter substructure. Our approach is generally applicable to caustic-crossing giant arcs in extragalactic lensing fields and can be incorporated into global lens modeling.

Ground-based radio astronomical observation at frequencies below 30 MHz is hampered by the Ionosphere and radio frequency interference (RFI). The Discovering Sky at the Longest wavelength (DSL) mission, also known as the Hongmeng mission, employs a linear array of satellites on a circular orbit around the Moon to make interferometric observations in this band. Though vastly different from the usual ground-based arrays, the interferometric visibility data collected by such an array is linearly related to the sky map, and the reconstruction is in principle an inversion problem of linear mapping. In this paper, we investigate a number of issues in the algorithm of global map reconstruction, focusing on the impact of sub-pixel noise induced by the finite pixelization of the sky, and errors due to regularization. We find that in the reconstruction process, if one builds up the beam matrix, which relates the sky pixels to the visibilities, by naively evaluating its elements at each of the pixel centers, then the sub-pixel noise can give rise to a significant aliasing effect. However, this effect can be effectively mitigated by a simple pixel-averaging method. Based on evaluation of the image quality using the correlation coefficient between the input and reconstructed map, and the signal-to-noise ratio, we discuss the selection strategy of the regularization parameter, and show that the sky can be well reconstructed with a reasonable choice of the regularization parameter.

Observations of transiting hot Jupiters have revealed a mismatch between the values of the Bond versus geometric albedos. In the planetary science literature, the ratio of these quantities is known as the phase integral. It has been extensively measured for the Solar System planets and shown to generally be non-unity in value. We use existing Cassini data of Jupiter to derive bandpass-integrated geometric albedos and phase integrals in the CHEOPS, TESS and Ariel bandpasses, demonstrating that these quantities vary markedly across these different wavelength ranges. By performing a population study of geometric albedos and phase integrals, we demonstrate that atmospheres with partial cloud cover may be identified using measurements of the phase integral if its measured uncertainty is $\sim 0.1$, which corresponds to an uncertainty of $\sim 3\%$ on the optical/visible secondary eclipse depth. The upcoming Ariel space mission will conduct an unprecedented statistical survey of cloud cover on hot Jupiters via the simultaneous measurement of $\sim 100$ infrared phase curves and optical secondary eclipses. Whenever available, the shape of optical phase curves of reflected light will directly constrain the phase integral, spherical albedo, degree of cloud cover and scattering asymmetry factor.

The line-of-sight magnetic field of galaxy clusters can be probed using Faraday rotation measure (RM) data. However, our understanding of cluster magnetism is limited due to the scarcity of polarized background radio sources, with most previous studies being constrained to $\sim 10$ sources per cluster. Leveraging the increased source density of the POlarisation Sky Survey of the Universe's Magnetism (POSSUM), we probe the magnetic field properties of the galaxy cluster Abell 3581 {(A3581)} with 111 RMs. We find that the standard deviation in the RM declines monotonically with increasing radius up to 0.75 Mpc{, agreeing with a radially declining magnetic field and electron density profile modeled as Gaussian and lognormal random fields, respectively. {We compare our observations of the inner 0.75 Mpc of A3581 to various semi-analytic models of the magnetic field and electron density, and obtain several best-fit models.} For the first time, we compare the observed RMs in a cluster to full magnetohydrodynamic simulated clusters from TNG-Cluster and find that the non-monotonic trend in RM standard deviation past 0.75 Mpc in A3581 is likely caused by past or present merger activity. We identify a possible candidate for a merger to be the galaxy group [DZ2015b] 276, which would be the first group detected in RMs that is not strongly emitting in X-rays. We find a possible merger axis of A3581 with this group at a position angle of $\theta = 52\pm 4$ deg.

Core-collapse supernovae mark the end of life of massive stars. However, despite their importance in astrophysics, their underlying mechanisms remain unclear. Neutrinos that emerge from the dense core of the star offer a promising way to study supernova dynamics. A strategy is presented to improve the potential of the KM3NeT neutrino telescope to detect core-collapse supernovae in our Galaxy or the Large Magellanic Cloud by further exploiting the properties of its optical modules equipped with multiple photomultipliers. A supernova burst is expected to produce a sudden hit rate increase in the KM3NeT detectors. New observables have been defined for individual optical modules that exploit the geometry and time distribution of the detected hits, enabling a better discrimination between signal and background signatures. In addition, a thorough investigation of the related systematic uncertainties is presented for the first time. When implemented, this new methodology allowed KM3NeT to probe 46% more Galactic core-collapse supernova candidates than with the previous trigger strategy, reaching the dense Galactic bulge. It is now expected that, once completed, KM3NeT will achieve full Galactic sensitivity to core-collapse supernovae.

We investigate nuclear reactions and feedback in hyperaccreting neutron star environments, considering accretion rates in the range 0.3 - $3\times10^4$ $M_\odot$ yr$^{-1}$, typical of short-period compact object binaries in common envelopes. Our mode ls account for weak reactions, neutrino energy loss, nuclear energy release, pair production, degenerate equations of state, and general relativistic hydrodynamics. Depending on accretion rates, these systems can develop both proton and neutron-rich atmospheres with strong convective instabilities linking the neutrino sphere to the outgoing accretion shock inside the radia tion trapping zone. Convection drives nucleons through multiple heating and cooling cycles, with photodisintegration dominat ing during the heating phase and heavy element synthesis during the cooling phase, ejecting material with abundances that dep end on the accretion rate and depth of the final decompression trajectory. The turbulent nature of convective currents is con ducive to creating a wide range of nuclear products through a variety of effects, including NSE freezeout and $r$, $p$ and $\ gamma$ processes. We also observe a novel multi-step process in reheated trajectories, consisting of proton-capture and photo-dissociation reactions operating on $r$-process seeds, producing overall neutron-deficient isotopes. A significant amount of infalling gas experiences high entropy and short (millisecond) freezeout timescales capable of making $r$-process elements w ith high over-abundances through a disequilibrium effect between neutrons and $\alpha$-particles that does not require an excess of neutrons.

Large photometric surveys will image billions of galaxies, but we currently lack quick, reliable automated ways to infer their physical properties like morphology, stellar mass, and star formation rates. Simulations provide galaxy images with ground-truth physical labels, but domain shifts in PSF, noise, backgrounds, selection, and label priors degrade transfer to real surveys. We present a preliminary domain adaptation pipeline that trains on simulated TNG50 galaxies and evaluates on real SDSS galaxies with morphology labels (elliptical/spiral/irregular). We train three backbones (CNN, $E(2)$-steerable CNN, ResNet-18) with focal loss and effective-number class weighting, and a feature-level domain loss $L_D$ built from GeomLoss (entropic Sinkhorn OT, energy distance, Gaussian MMD, and related metrics). We show that a combination of these losses with an OT-based "top_$k$ soft matching" loss that focuses $L_D$ on the worst-matched source-target pairs can further enhance domain alignment. With Euclidean distance, scheduled alignment weights, and top-$k$ matching, target accuracy (macro F1) rises from $\sim$46% ($\sim$30%) at no adaptation to $\sim$87% ($\sim$62.6%), with a domain AUC near 0.5, indicating strong latent-space mixing.

The Magellanic Stream (MS) is a vast gaseous structure in the Milky Way halo, containing most of its mass in ionized form and tracing the interaction between the Large and Small Magellanic Clouds and the Galaxy. Using HST/COS G160M spectra from the Plane Quasar Survey, we detect C IV absorbers likely associated with the MS, extending to the northern side of the Galactic plane, approximately 60$^\circ$ beyond its previously known ionized extent. These absorbers exhibit position and kinematic alignment and show consistent ionization trends with previously studied MS sight lines. The non-detection of low ions such as Al II and Si II, and the detection of C IV (and Si IV in some sightlines), indicates a highly ionized gas phase. The observed Si IV/C IV column density ratios suggest a gas temperature of $T \sim 10^{5.3}$~K and favor collisional ionization over photoionization. We estimate the newly detected extension increases the previous ionized gas mass of the MS, and its coherent kinematics suggest that it was stripped within the past few hundred Myr and has not yet mixed with the Milky Way halo. The existence of highly-ionized MS gas at a location above the Galactic Plane may constrain the orbital direction of the Magellanic Clouds.

The environmental density of galaxies within the cosmic web constrains their 3D locations in filaments, voids, groups, and clusters. It traces the distribution of baryons and the influence of dark-matter halos on galaxy evolution, and helps diagnose external processes such as supernova and AGN feedback, tidal interactions, ram-pressure stripping, and large-scale flows. We focus on isolated barred galaxies, the parent population that includes Milky Way analogs. To measure local environmental density and verify isolation (|Delta v| <= 500 km s^-1), we built a Python pipeline that works in two redshift regimes: low (z0 < 0.02) and high (z0 >= 0.02). Densities were estimated with k-nearest neighbors and Voronoi tessellations and classified as void (Sigma < 0.05), filament (0.05 <= Sigma < 0.5), group (0.5 <= Sigma < 2.0), and cluster (Sigma >= 2.0). Our northern-sky sample contains 311 isolated barred galaxies from 2MIG plus Milky Way analog systems (z < 0.07). We find 157 void, 84 filament, 27 group, and 11 cluster galaxies; 30 have no detected neighbors. Sixty-seven lie in extremely low-density regions (Sigma_3D < 0.01 gal Mpc^-3), and 22 have the nearest companion beyond 5 Mpc, indicating residence in extended voids. The Milky Way (Sigma_5NN ~ 0.13 gal Mpc^-3, R ~ 2.1 Mpc) and its close analog NGC 3521 both lie in filamentary environments at the edge of a nearby void. For z > 0.02 we identify three additional Milky Way analog candidates from the density metrics: CGCG 208-043 (3D Voronoi) and NGC 5231 and CGCG 047-026 (5th-nearest neighbor). These results show that local environmental density is an effective, physically motivated criterion for selecting Milky Way analogs.

The analysis of the optical spectra of PMN J0948+0022 showed significant variations in the spectral lines that, when combined with the Fermi $\gamma$-ray light curve and radio observations reported by other authors, were interpreted as the result of interactions between the relativistic jet and the narrow-line region (NLR). In this work, we present order-of-magnitude calculations of the energetics associated with this proposed jet-NLR interaction. We demonstrate that the observed outflows are capable of absorbing a fraction of the jet energy and converting it into kinetic energy. This mechanism provides a natural explanation for the optical spectral variability recorded with the X-shooter and Multi-Unit Spectroscopic Explorer (MUSE) instruments. Our results support the scenario in which feedback from the relativistic jet can dynamically influence the circumnuclear gas, offering new insights into the coupling between jets and the NLR in $\gamma$-ray-emitting narrow-line Seyfert 1 galaxies.

Primordial Black Holes (PBHs) provide a powerful probe of early-universe physics, linking inflationary fluctuations to observable cosmological phenomena. In this work, we use a bottom-up approach to study how PBHs with masses in the range $10^{8} \leq M \leq 10^{13}\,\mathrm{g}$ modify Big Bang Nucleosynthesis (BBN) through Hawking radiation. We incorporate PBH evaporation into a reaction-network code to evaluate its impact on light-element abundances. Our analysis shows that PBH evaporation acts as an entropy injection mechanism, increasing the comoving entropy density. To reproduce the observed comoving entropy density per baryon $(s/n_{\mathrm{b}})$ from the CMB, BBN simulations must therefore begin with a smaller initial entropy than in the standard scenario without PBHs. The results also reveal a threshold near $M \approx 10^{10}\,\mathrm{g}$ that separates two distinct regimes of BBN behavior. As an example, for $M \geq 10^{10}\,\mathrm{g}$, the $^4{\mathrm{He}}$ mass fraction $Y_{\mathrm{P}}$ increases monotonically with $\beta$, driven by the enhanced Hubble expansion from PBH energy density. In contrast, for $M \leq 10^{10}\,\mathrm{g}$, $Y_{\mathrm{P}}$ exhibits non-monotonic behavior shaped by the timing of PBH evaporation and its influence on nuclear reaction rates. These findings highlight the sensitivity of BBN to PBH evaporation and establish a framework for understanding how PBH populations influence the thermal history of the early universe.

We present the first local simulations of disc breaking/tearing in a warped accretion disc. Warps can arise due to a misalignment between the disc and the rotation axis of the central object, or a misalignment with the orbital plane of a binary (or planetary) companion. Warped discs can break into rings, as found in observations of circumbinary protoplanetary discs and global simulations of tilted discs around spinning black holes. In this work we isolate the mechanism of disc breaking in high-resolution, quasi-2D, local (shearing box), hydrodynamic simulations of a Keplerian disc. We consider the evolution of a free (unforced) warp in the wavelike ($\alpha < H/r$) regime. At large warp amplitudes ($\psi_{\text{max}} \gtrsim 1$) the disc breaks into four rings on timescales of around 20 orbits which are separated by gaps of around $\sim 10H_0$. The warp exhibits a rich tapestry of small-scale dynamics, including horizontal sloshing motions, vertical oscillations or bouncing, warp steepening, and shocks. The shocks act as a source of enhanced dissipation which facilitates gap opening and thus disc breaking. At smaller warp amplitudes $\psi_{\text{max}} \lesssim 1$, for which we also develop a quasi-linear theory, the disc does not break, but instead exhibits hydrodynamic parametric instability. We also investigate the effect of viscosity: at small warp amplitudes the parametric instability is damped and the warp propagates as a pure bending wave, while at large warp amplitudes the emerging gaps are partially filled by viscous diffusion.

The nature of dark energy is one of the big puzzling issues in cosmology. While $\Lambda$CDM provides a good fit to the observational data, evolving dark energy scenarios, such as the CPL parametrization, offer a compelling alternative. In this paper, we present a Bayesian model comparison of various dark energy parametrizations using a joint analysis of Cosmic Microwave Background data, DESI Baryon Acoustic Oscillation measurements, and the PantheonPlus (or Union3) Supernovae type Ia sample. We find that while the $\Lambda$CDM model is initially favored over a constant $w$CDM model, the CPL parametrization is significantly preferred over $w$CDM, reinforcing recent evidence for an evolving dark energy component, consistent with DESI collaboration findings. Crucially, when testing higher-order CPL extensions, the so-called CPL$^+$ and CPL$^{++}$, our Bayesian analysis shows that the observational data do not favor these more complex scenarios compared to the standard CPL. This result indicates that adding excessive complexity to the CPL form is unwarranted by current observations. Interestingly, similar to the CPL parametrization, alternative two-parameter forms, specifically $w_{de}(a) = w_0 + w_b(1-a)^2$ and $w_{de}(a) = w_0 + w_c(1-a)^3$, yield a better fit to observational data than the standard $\Lambda$CDM cosmology. Our results challenge the necessity for overly complex CPL extensions and confirm that well-chosen two-parameter $w_0w_a$ parametrizations effectively capture DE evolution with current cosmological data, supporting the recent signals for dynamical dark energy by DESI collaboration.

The detection of fast radio bursts (FRBs) in both young and old stellar populations suggests multiple formation pathways, beyond just young magnetars from core-collapse supernovae. A promising delayed channel involves the formation of FRB-emitting neutron stars through merger- or accretion-induced collapse of a massive white dwarf (WD). By simulating a realistic stellar population with both binaries and triples, we identify pathways to WD collapse that could produce FRB candidates. We find that (i) triple dynamics open new merger channels inaccessible to isolated binaries, significantly enhancing the overall merger rate; (ii) triple-induced mergers broaden the delay-time distribution, producing long-delay ($\gtrsim1$-8~Gyr) events largely independent of metallicity, alongside a shorter-delay population ($\lesssim100$~Myr) of rapid mergers; (iii) these long delays naturally yield FRBs in older environments such as quiescent host galaxies and galactic halos; (iv) when convolved with the cosmic star-formation history, binary channels track the star-formation rate ($z_{\rm peak} \sim 2$), while triple channels peak later ($z_{\rm peak} \sim 1$), giving a combined local source rate of $R_0 \approx 2\times10^4~{\rm Gpc^{-3}~yr^{-1}}$, consistent with observations; and (v) applying the same framework to Type~Ia supernovae, we find that triples extend the delay-time tail and roughly double the Ia efficiency relative to binaries, yielding rates and redshift evolution in good agreement with observations. If FRBs originate from the collapse of WDs, our results establish triples, alongside binaries, as a crucial and previously overlooked formation pathway whose predicted rates, host demographics, and redshift evolution offer clear tests for upcoming surveys.

The trace anomaly $\Delta$ is an important quantity that measures the broken conformal symmetry in neutron star matter. In this work, we present quasi-universal relations that connect the stellar profile of $\Delta$ to the compactness, moment of inertia, and tidal deformability of neutron stars. We apply the quasi-universal relations to determine the trace anomaly profiles for PSR J0030+0451 and PSR J0740+6620 based on their mass-radius measurements. We also analyze PSR J0737-3039A according to its moment of inertia inferred from Bayesian modeling of nuclear equation of state. A recent multimessenger constraint on the tidal deformability is also studied, resulting in an estimate value of the trace anomaly $\Delta_c = 0.1770^{+0.0365}_{-0.0432}$ at the center of a $1.4M_\odot$ canonical neutron star. It is expected that more precise observations from both electromagnetic and gravitational-wave channels in the future will provide tighter constraints on the behavior of $\Delta$ inside neutron stars.

Abridged. The measured braking index, $n=\nu \ddot{\nu}/\dot{\nu}^2$, of a rotation-powered pulsar with spin frequency $\nu$ and braking torque $K \nu^{n_{\rm pl}}$, features secular and stochastic anomalies arising from $\dot{K} \neq 0$ and random torque noise respectively. Previous studies quantified the variance $\langle n^{2} \rangle = (n_{\rm pl}+\dot{K}_{\rm dim})^{2}+\sigma_{\rm dim}^{2}$, where the secular anomaly, $\dot{K}_{\rm dim}$, is inversely proportional to the characteristic time-scale $\tau_{K}$ over which $K$ varies; the stochastic anomaly, $\sigma_{\rm dim}^{2} = \sigma_{\ddot{\nu}}^{2}\nu^{2}\gamma_{\ddot{\nu}}^{-2}\dot{\nu}^{-4}T_{\rm obs}^{-1}$, is a function of the timing noise amplitude $\sigma_{\ddot{\nu}}$, a damping time-scale $\gamma_{\ddot{\nu}}^{-1}$ and the total observing time $T_{\rm obs}$; and the average is taken over an ensemble of random realizations of the noise process. Here, we use a hierarchical Bayesian scheme, based on the formula for $\langle n^{2} \rangle$, to infer the population-level distribution of $n_{\rm pl}+\dot{K}_{\rm dim}$ for a sample of $68$ young radio pulsars, observed for $\gtrsim 10~{\rm years}$ with Murriyang, the 64-m Parkes radio telescope. Upon assuming that the $n_{\rm pl}+\dot{K}_{\rm dim}$ values are drawn from a population-level Gaussian, $N(\mu_{\rm pl}, \sigma_{\rm pl})$, the Bayesian scheme returns the mean $\mu_{\rm pl} = 9.95^{+5.58}_{-5.26}$ and standard deviation $\sigma_{\rm pl}=10.89^{+5.14}_{-3.69}$. At a per-pulsar level it returns posterior medians satisfying $-13.86 \leq n_{\rm pl}+\dot{K}_{\rm dim} \leq 30.38$. The secular anomaly dominates the stochastic anomaly, with posterior medians satisfying $|n_{\rm pl} + \dot{K}_{\rm dim}| \geq \sigma_{\rm dim}$ in 10 out of 68 objects.

Measuring the occurrence rates of celestial objects is a valuable way to study their origins and evolution. Giant planets and brown dwarfs produce large Doppler signatures that are easily detectable by modern instrumentation, and legacy radial velocity (RV) surveys have now achieved full orbital coverage for periods $\lesssim$30 years. However, the Doppler method's sensitivity to companion minimum mass $M_c \sin i$ -- as opposed to true mass $M_c$ -- prevents unambiguous characterization using RVs alone, as purported giant planets may be brown dwarfs or stars on inclined orbits. Here we combined legacy RVs with absolute astrometry to re-fit the orbits of 195 companions from the California Legacy Survey. Nearly 50% (8/18) of the ``brown dwarfs" ($M_c \sin i$=13--80 $M_{Jup}$) we refit had true masses above 80 $M_{Jup}$. We incorporated our orbital posteriors and target sensitivity maps into a Poisson likelihood model to calculate occurrence as a function of true companion mass $M_c$ (0.8--80 $M_{Jup}$) and separation $a$ (0.3--30 AU). The semi-major axis distributions of objects in this range vary smoothly with mass, with Jupiter analogs favoring an abrupt increase in occurrence near 1 AU and brown dwarfs exhibiting a gradual enhancement at wider separations. Marginalized companion occurrence between 1--10 AU decreases smoothly with mass, with brown dwarfs having the lowest occurrence rate: $0.9^{+0.5}_{-0.4}$%. Jupiter analogs are 10 times as common as brown dwarfs per mass interval in this range, demonstrating that the brown dwarf desert extends to 10 AU. The smooth variation in these distributions disfavors a sharp transition mass between ``bottom-up" core accretion and ``top-down" gravitational instability formation mechanisms, and rather suggests that these processes may produce companions in overlapping mass ranges.

We investigate the impact of a low-mass planet on dust coagulation, and its consequent feedback on planetary migration, using a linear analysis of the coupled dust-gas hydrodynamic equations. Dust coagulation is incorporated via a single-size approximation. In the co-orbital region of the planet, we find that the growth of dust size is significantly suppressed by planet-induced coagulation modes (CMs). This effect are less pronounced with smaller stopping times, stronger gaseous turbulence or imperfect sticking. Regarding planetary migration, we find that CMs make outward migration require $\tau \gtrsim 0.3$ ($\tau$ is dimensionless stopping time) with typical turbulent strength and dust coagulation efficiency. We demonstrate that the torque variations are reasonable and arise from phase shifts between the density and stopping time perturbations in the coagulation modes.

We revisit the tidal evolution of the WASP-12 system using direct numerical calculations with the GYRE-tides code. WASP-12b is a hot Jupiter on a 1.1-day orbit around a slightly evolved F-type star. Its observed orbital decay rate, $|\dot{P}_{\rm orb}/P_{\rm orb}| \approx 3.2\,\mathrm{Myr}^{-1}$, provides a strong constraint on stellar tidal dissipation. We confirm that linear tides with radiative damping and convective damping, as currently implemented, are not sufficient to reproduce the observed inspiral timescale. Nevertheless, our calculations, based on fully non-adiabatic forced oscillations in MESA stellar models with convective envelopes, yield dissipation rates that are consistent with previous semi-analytic and adiabatic estimates, confirming the robustness of our numerical framework. As the only open-source, actively maintained tool capable of computing orbital evolution in exoplanet systems, GYRE-tides provides a benchmark calculation for WASP-12 and future applications. Our results validate GYRE-tides as a tool for analyzing combined radiative and convective damping, and indicate that the observed decay rate requires tidal dissipation operating in or near the fully damped regime, which may be achieved through nonlinear damping. These contributions could also be evaluated by computing the wave luminosity at the radiative-convective boundary using our tool. GYRE-tides offers an open-source framework for computing tidal dissipation in short-period exoplanet systems, including the many systems expected to show orbital decay in upcoming Roman surveys.

In the outer solar system beyond Jupiter, water ice is a dominant component of planetary bodies, and most solid objects in this region are classified as icy bodies. Icy bodies display a remarkable diversity of geological, geophysical, and atmospheric processes, which differ fundamentally from those of the rocky terrestrial planets. Evidence from past and ongoing spacecraft missions has revealed subsurface oceans, cryovolcanic activity, and tenuous but persistent atmospheres, showing that icy bodies are active and evolving worlds. At the same time, major questions remain unresolved, including the chemical properties of icy materials, the geological histories of their surfaces, and the coupling between internal evolution and orbital dynamics. Current knowledge of the surfaces, interiors, and atmospheres of the principal icy bodies is built on spacecraft measurements, telescopic observations, laboratory experiments, and theoretical modeling. Recent contributions from Juno, JWST, and stellar occultation studies have added valuable constraints on atmospheric composition, interior structure, and surface activity. Looking ahead, missions such as JUICE, Europa Clipper, Dragonfly, and the Uranus Orbiter and Probe are expected to deliver substantial progress in the study of icy bodies. Their findings, combined with continued Earth- and space-based observations and laboratory studies, will be critical for assessing the potential habitability of these environments and for placing them within a broader framework of planetary system formation and evolution.

With the exponential growth of time-domain surveys, the volume of light curves has increased rapidly. However, many survey projects, such as Gaia, still rely on offline batch-processing workflows in which data are calibrated, merged, and released only after an observing phase is completed. This latency delays scientific analysis and causes many high-value transient events to be buried in archival data, missing the window for timely follow-up. While existing alert brokers handle heterogeneous data streams, it remains difficult to deploy a unified framework that combines high-performance incremental storage with real-time classification on local infrastructure. To address this challenge, we propose TDLight, a scalable system that adapts the time-series database TDengine (a high-performance IoT database) for astronomical data using a one-table-per-source schema. This architecture supports high-throughput ingestion, achieving 954,000 rows s^-1 for archived data and 541,000 rows s^-1 for incremental streams, while Hierarchical Equal Area isoLatitude Pixelization (HEALPix) indexing enables efficient cone-search queries. Building on this storage layer, we integrate the pre-trained hierarchical Random Forest classifier from the LEAVES framework to construct an incremental classification pipeline. Using the LEAVES dataset, we simulate data accumulation and evaluate a trigger-based strategy that performs early classification at specific observational milestones. In addition, by monitoring the evolution of classification probabilities, the system identifies "high-value candidates" -- sources that show high early confidence but later undergo significant label shifts. TDLight is released as an open-source Dockerized environment, providing a deployable infrastructure for next-generation time-domain surveys.

It is assumed that heavy dark matter $\phi$ with O(TeV) mass captured by the Earth may decay to relativistic light milli-charged particles (MCPs). These MCPs could be measured by the IceCube neutrino telescope. The massless hidden photon model was taken for MCPs to interact with nuclei with a running electromagnetic coupling constant, so that the numbers and fluxes of expected MCPs may be evaluated at IceCube. Meanwhile, the numbers of expected neutrino background events were also evaluated at IceCube. Based on the assumption that no events are observed at IceCube in 10 years, the corresponding upper limits on MCP fluxes were calculated at 90\% C. L.. These results indicated that the MCPs from the Earth's core could be directly detected at O(1TeV) energies at IceCube when $5.65\times10^{-5}\lesssim\epsilon^2\lesssim1.295\times10^{-3}$. And a new region of 4 GeV < $m_{MCP}$ < 100 GeV and $5.51\times10^{-2}$ < $\epsilon$ < 0.612 is ruled out in the $m_{MCP}$-$\epsilon$ plane with 10 years of IceCube data.

The earliest phase of the Milky Way's evolution involved a transition from a dispersion-supported proto-galaxy to a rotationally supported disk. A key chemical signature of this transition is the moderate rise in [$\alpha$/Fe] near $\mathrm{[Fe/H]}\approx-1.3$, which we previously interpreted as evidence for $\alpha$-enhanced gas accretion fueling early disk formation. However, this trend alone does not uniquely constrain the trade-off between initial gas mass, inflow rate, and star formation efficiency (SFE), leaving the physical condition of the proto-Milky Way uncertain. To break this degeneracy, we analyze the metal-poor tail ($-3<\mathrm{[Fe/H]}<-2$) of the Milky Way's metallicity distribution function (MDF) using Gaia DR3 BP/RP (XP) metallicities from ten catalogs. After applying recommended quality cuts, all catalogs exhibit a single-slope exponential tail with slopes $k\simeq0.5$--2.0. Comparison with one-zone galactic chemical-evolution (GCE) models that replicated the [$\alpha$/Fe]-rise from Paper I shows that shallow tails ($k\simeq0.6$) require a massive initial cold gas reservoir ($\gtrsim10^9\, \mathrm{M_\odot}$), while steeper tails ($k\gtrsim1$) arise from small reservoirs that built up over time with weak inflow. MDFs with $k \simeq 1.0$ are best reproduced under our GCE framework, which favor a proto-Galaxy with a moderate gas reservoir ($10^{8}$--$10^{9}\, \mathrm{M_\odot}$) sustained through weak continuous inflow ($\sim 2 \ \mathrm{M_\odot\,yr^{-1}}$) and SFE comparable to today's value (a few $\times 10^{-10}\,\mathrm{yr^{-1}}$) during the first Gyr. This scenario is reinforced by MDFs of 30 Milky Way analogs in the Auriga simulations, which exhibit similar slopes ($k\approx1.2$). The metal-poor MDF tail thus provides a quantitative constraint on the Milky Way's early gas accretion and star formation history.

This study broadens our comprehensive investigation of total-eclipse W Ursae Majoris-type contact binaries by analyzing eight additional systems, continuing our previous research. Multiband $BVR_cI_c$ photometric data were obtained at an observatory in Mexico, from which new times of minima were determined. All target systems also had available space-based TESS time-series data. Orbital period variations were studied for eight target systems, showing either linear or parabolic trends. The target systems exhibiting parabolic trends demonstrated a sustained decrease in their orbital periods over time. We modeled the light curves utilizing the PHOEBE Python code in combination with the BSN application. We revisited the relationship between orbital period and the temperature of the hotter component in contact binary systems using an empirical approach. Our analysis identified a clear break at P=0.27 days, separating the systems into two distinct groups for orbital periods shorter than 0.6 days. Following the determination of stellar extinction, absolute parameters for seven systems were estimated employing parallax measurements from Gaia DR3. Based on the components' effective temperatures and masses, the systems were classified into A- and W-subtypes. Their evolutionary states were illustrated using mass-radius and mass-luminosity diagrams.

Some of the most interesting insights into solar physics and space weather come from studying radio emissions associated with solar activity, which remain inherently unpredictable. Hence, a real-time triggering system is needed for solar observations with the versatile new-generation radio telescopes to efficiently capture these episodes of solar activity with the precious and limited solar observing time. We have developed such a system, Solar Triggered Observations of Radio bursts using MWA and Yamagawa (STROMY) for the Murchison Widefield Array (MWA), the precursor for the low frequency telescope of upcoming Square Kilometre Array Observatory (SKAO). It is based on near-real-time data from the Yamagawa solar spectrograph, located at a similar longitude to the MWA. We have devised, implemented, and tested algorithms to perform an effective denoising of the data to identify signatures of solar activity in the Yamagawa data in near real-time. End-to-end tests of triggered observations have been successfully carried out at the MWA. STORMY is operational at the MWA for the routine solar observations, a timely development in the view of the ongoing solar maximum. We present this new observing framework and discuss how it can enable efficient capturing of event-rich solar data with existing instruments, like the LOw Frequency ARray (LOFAR), Owens Valley Radio Observatory - Long Wavelength Array (OVRO-LWA) etc., and pave the way for triggered observing with the SKAO, especially the SKA-Low.

Ultralight axions (ULAs) are hypothetical particles which can behave like dark matter (DM) or dark energy (DE) depending on masses generated at the symmetry-breaking scale. It remains a mystery whether the ULAs can make up a fraction of DM or DE. Although theoretical predictions indicate that the ULAs may leave distinct imprints on cosmological signals, these signatures may exist in a broad spatial and temporal scales, and may be degenerate with the known effects of the standard model. The ULA signatures are extremely subtle and the observational evidence of the ULAs remain elusive. In this work, we infer the ULA properties using both the early and late universe observations from the cosmic microwave background (CMB) and baryon acoustic oscillations (BAO). We validate modeling of the ULA effects using the CMB and BAO mock data and perform different tests to cross-check the results. By analyzing the Planck 2018 CMB measurements and the BAO measurements from the Data Release 2 of Dark Energy Spectroscopic Instrument (DESI), we constrain the energy density fraction ratio of the ULAs to total dark matter $\Omega_a/\Omega_d$ and obtain a new upper bound of $\Omega_a/\Omega_d$. Future CMB and BAO measurements will achieve unprecedented precision and will be crucial for understanding the nature of the ULAs.

Whether planetary motions influence the solar magnetic cycle has remained an open question due to the lack of a rigorous physical mechanism. Here we develop a Lagrangian framework based on the Virial theorem to show that planetary orbital angular momentum modulates the Sun's rotation through barycentric dynamics, consistent with angular momentum conservation. Using Singular Spectrum Analysis (SSA) applied to sunspot number (SSN) and terrestrial length-of-day (LOD) records spanning the past 250 years, we identify two non-linear trends plus 11 common pseudo-cycles whose periods match those of planetary resonances (LR$^{*}$). We then demonstrate that a reduced nonlinear $\alpha-\Omega$ dynamo model, forced solely by the summed planetary right ascension RA(t), reproduces both the Schwabe cycle and its multi-decadal envelope, including the Dalton Minimum and the Modern Maximum. These results provide strong evidence, based on a physically-grounded model, that the solar dynamo is weakly forced by planetary motion, suggesting that the solar dynamo, classically viewed as autonomous, may in fact behave as a synchronized system.

With their jet pointing towards us, blazars are ideal tools to study the physics and structure of extragalactic jets. Their powerful jets are cosmic particle accelerators and are alleged to be one of the production sites of the high-energy neutrinos detected by the IceCube Observatory. Doppler beaming of the jet nonthermal radiation increases blazar brightness, blue-shifts their emission, and shortens their variability time scales, which are observed to range from years down to minutes. This review will focus on blazar flux, spectral, and polarization variability across the electromagnetic spectrum. Interpretation of blazar variability calls into question both intrinsic and extrinsic mechanisms. Shock waves, magnetic reconnection, and turbulence can accelerate particles inside the jet, while jet precession, rotation, and twisting can produce variations in Doppler beaming. Microvariability puts strong constraints on the size of the emitting regions, suggesting a multizone emitting jet. Twisting jets have been proposed to explain the long-term multiwavelength variability. They are supported by radio observations of bent or helical jets, and by results of relativistic magnetohydrodynamics simulations of plasma jets. Detection of (quasi)periodic behaviour has been ascribed to orbital motion in black hole binary systems, jet precession, kink instabilities developing inside the jet, or perturbations in the accretion disc. Gravitational microlensing has been suggested to explain blazar behaviour in some cases. Polarization provides information on the structure and behaviour of the magnetic field in the emission zones. Both the degree and angle of polarization can show strong and fast variability, which is sometimes correlated with flux. The interpretation of flux, spectral, and polarization variability within a consistent picture challenges current models of blazar variability.

The inner regions of planet-forming disks hold invaluable insights for our understanding of planet formation. The disk around the Herbig star HD 100453 presents one such environment, with an inner disk that is significantly misaligned with respect to the outer disk. This paper expands the existing H-band (PIONIER) and K-band (GRAVITY) interferometric studies of the HD 100453 inner disk to the L-band with the MATISSE VLTI instrument. With snapshot data spanning approximately four years we aim for a more comprehensive understanding of the inner disk structures and their potential time evolution. Based on the MATISSE data obtained, we use a combination of analytical models and image reconstruction to constrain the disk structure. Additionally, we fit a temperature gradient model to the selected wavelength range of PIONIER, GRAVITY and MATISSE to derive physical properties of the inner regions. Our parametric model finds an inclination of $\sim 47.5^\circ$ and a position angle of $\sim 83.6^\circ$, which corroborates the case of strong inner-outer disk misalignment. From the symmetric temperature gradient we derive an inner disk radius around $\sim0.27$ au, with dust surface densities of $\Sigma_{\rm{subl}} \approx 10^{-3.2}$ g/cm$^2$ and vertical optical depth $\tau_{\rm{z, subl}} \approx 0.1-0.06$. Same-night MATISSE and GRAVITY observations indicate the necessity for higher-order asymmetries to explain the interferometric signals, which is further supported by a MATISSE snapshot image reconstruction. The chromatic interferometric data reveal the likely need for higher-order asymmetries to explain the inner disk of HD~100453, suggesting a possible origin in dynamic interactions or disk instabilities. Coordinated multi-wavelength infrared interferometric observations with GRAVITY and MATISSE will be crucial to confirm these findings and uncover its underlying nature.

Transmission spectroscopy of sub-Neptunes was expected to reveal their compositions and hence origins, yet many show flat near- to mid-infrared spectra. Such spectra can be explained either by metal dominated atmospheres or by high-altitude, grey aerosols. Observations of escaping hydrogen and helium from several of these planets rule out metal dominated atmospheres, while homogeneous distributions of small aerosols cannot produce flat spectra and large particles require unphysically high production rates. We investigate the role of heterogeneous, "clumpy" aerosol distributions in shaping transmission spectra. Modestly optically thick clumps at high altitudes can produce flat spectra even with small particles and physically realistic production rates. Clumping increases the effective photon mean-free path while reducing wavelength dependence, allowing the aerosol distribution to behave as an effective grey absorber. Applying this framework to the sub-Neptune TOI-776c, we show that clumpy aerosols can reconcile the observed flattening of its transmission spectrum with a primordial H/He-dominated atmosphere. We further discuss implications for emission spectra, where enhanced stellar radiation penetration and altered scattering in a clumpy medium could produce observable signatures. These results suggest that clumpy aerosol distributions naturally resolve the tension between flat spectra and low-metallicity atmospheres and may be a common feature of sub-Neptune exoplanets. More broadly, our results highlight the need to consider aerosol heterogeneity when interpreting high-altitude microphysics and the spectral appearance of exoplanet atmospheres with JWST, and motivate theoretical work to identify the physical mechanisms capable of generating clumpy aerosol distributions.

Traditionally, fast radio transient searches are conducted on dedispersed time series using thresholding techniques based on the statistical properties of the data. However, peaks in dedispersed time series do not directly provide information on the nature of the source. In the DM-time domain, the S/N variation of real, dispersed astrophysical signals forms a characteristic bow tie shape, whereas radio frequency interference (RFI) can take multiple different forms. We have developed a method that bypasses the thresholding step of traditional single-pulse searches in favour of a direct DM-time domain image analysis. The backbone of our pipeline is a Mask R-CNN, a deep learning model designed for object detection, enabling it to identify the bow tie signature and distinguish real sources from RFI. Previous deep learning models often include a snippet of the DM-time domain in their input. We have trained the model on simulated bursts injected on top of real MeerKAT noise observations. We tested the model on MeerKAT follow-up observations of the repeater FRB20240114A and we were able to recover all bursts with a signal-to-noise above the traditional threshold, while detecting two bursts that were fainter. Our new approach considerably reduces the number of candidates above a nominal threshold while being capable of running in real time for typical surveys. We also propose a modified version of the traditional dedispersion plan optimised for this method.

Tidal disruption events (TDEs) correspond to the destruction of a star by the tidal forces around a black hole, leading to outbursts which can last from months to years. These transients are rare, and increasing the current sample is paramount to understand them. As part of the Fink alert broker, we have developed an early detection system for TDEs for the Zwicky Transient Facility (ZTF) data. In this paper, we report on the optical transients we found either during the development of this tool, or when applying the classifier to the existing archive. We use this sample to anticipate what improvements to the TDE detection systems will need to be implemented for future surveys. For all the transients, we present optical and infrared archival photometry from ZTF, WISE, and Catalina, and assess the previous nuclear activity of the host. We fit the ZTF lightcurves with both a phenomenological and a physically-motivated model. We report on a total of 19 optical nuclear transients, out of which nine are in passive galaxies, eight in active galaxies, and two for which the activity of the host is uncertain. Two transients are newly discovered repeated TDE candidates, and we compare them to the current sample of repeated optical nuclear transients. One transient is exceptionally long-lived (over 5 years), in an until-now passive galaxy. Three of the TDE-like flares in active galaxies have absolute g-band magnitudes brighter than -24, making them new Extreme Nuclear Transient (ENT) candidates. One seemingly repeated object was revealed to be two independent supernovae in the same galaxy. This sample shows both the potential of our detection system for future discovery, and the relevance of archival searches to reveal overlooked transients. It also raises several points of concern and avenues of improvement for current and future classifiers.

Kodaikanal Solar Observatory (KoSO) is one of the oldest solar observatories, possessing an archive of multi-wavelength solar observations, including white light, Ca II K, and H-alpha images spanning over a century. In addition to these observations, KoSO has preserved hand-drawn suncharts (1904-2022), on which various solar features such as sunspots, plages, filaments, and prominences are marked on the Stonyhurst grid with distinct colour coding. In this study, we present the first comprehensive result that includes the entire data set from these suncharts using a supervised Machine Learning model called "Convolutional Neural Networks (CNNs)", firstly to identify the solar disks from the charts (1909-2007), secondly to identify the plages, spanning 9 solar cycles (1916-2007). We train the CNN with the manually identified solar disk and plage. We first detect the solar limb and the North-South line in the suncharts, which enables the extraction of disk centre coordinates, radius, and P-angle. Following that, we use a CNN similar architecture to achieve accurate image segmentation for the identification of plages. We compare plage areas derived from the suncharts with those obtained from Ca II K full-disk observations, and find good agreement that demonstrates the potential application of such an ML technique for historical data. The results of this study further demonstrate the potential application of sunchart data to fill the existing data gaps in the KoSO multi-wavelength observations and contribute toward constructing a composite series over the last century.

Chondrules are spherical or subspherical particles of crystallized or partially crystallized liquid silicates that constitute large-volume fractions of most chondritic meteorites. Chondrules typically range $0.1-2\,$mm in size and solidified with cooling rates of $10-1000\,{\rm K\,h^{-1}}$, yet these characteristics prove difficult to reconcile with proposed formation models. We numerically show that collisions among planetesimals containing volatile material naturally explain both the sizes and cooling rates of chondrules. We show that the high-velocity collisions with volatile-rich planetesimals first induced in the solar nebula by Jupiter's formation produced increasing amounts of silicate melt for increasing impact velocities above $2\,{\rm km\,s^{-1}}$. We propose that the expanding gas formed from volatile materials by collisional heating dispersed and cooled the silicate melt, resulting in droplet sizes and cooling rates consistent with the observed sizes and inferred cooling rates. We further show that the peak melt production is linked to the onset of Jupiter's runaway gas accretion, and argue that the peak age of chondrules points to Jupiter's birth dating 1.8 Myr after CAIs.

We tested whether simultaneous spectropolarimetric imaging in two magnetically sensitive optical spectral lines, which probe two different layers of the solar atmosphere (the photosphere and the temperature minimum), can help constrain the depth variation of horizontal flows. We first tested the feasibility of our method using Fourier local correlation tracking (FLCT) to track physical quantities at different optical depths ($\log\tau_{500}={-1,-2,-3,-4}$) in an atmosphere simulated with the MURaM code. We then inferred the horizontal distribution of the LOS magnetic field component from synthetic spectropolarimetric observations of Fe I 525.0 nm and Mg I b2 spectral lines, applied FLCT to the time sequence of these synthetic magnetograms, and compared our findings with the original height-dependent horizontal velocities. Tracking the LOS magnetic field component (which coincides with the vertical component at the disk center) yields horizontal velocities that, after appropriate temporal and spatial averaging, agree excellently with the horizontal component of the simulated velocities, both calculated at constant $\tau_{500}$ surfaces, up to the temperature minimum ($\log\tau_{500}=-3$). When tracking the temperature at constant $\tau_{500}$ surfaces, this agreement already breaks down completely at the mid photosphere ($\log\tau_{500}=-2$). Tracking the vertical component of the magnetic field inferred from synthetic observations of the Fe I 525.0 nm and the Mg I b2 spectral lines yields a satisfactory inference of the horizontal velocities in the mid-photosphere ($\log\tau_{500}\approx-1$) and the temperature minimum ($\log\tau_{500}\approx-3$), respectively. Our results indicate that high-spatial-resolution spectropolarimetric imaging in solar spectral lines can provide meaningful information about the horizontal plasma velocities over a range of heights.

Solar flares are explosive releases of magnetic energy stored in the solar corona, driven by magnetic reconnection. These events accelerate electrons, generating hard X-ray emissions and often display Quasi Periodic Pulsations (QPPs) across the energy spectra. However, the energy transfer process remains poorly constrained, with competing theories proposing different acceleration mechanisms. We investigate electron acceleration and transport in a flaring coronal loop by solving a time-dependent Fokker-Planck equation. Our model incorporates transient turbulent acceleration, simulating the effects of impulsive energy input to emulate the dynamics of time-dependent reconnection processes. We compute the density-weighted electron flux, a diagnostic directly comparable to observed X-ray emissions, across the energy and spatial domains from the corona to the chromosphere. We investigate different time-dependent functional forms of the turbulent acceleration, finding that the functional form of the acceleration source maintains its signature across energy bands (1 to 100 keV) with a response time that is energy dependent (with higher energy bands displaying a longer response time). In addition, we find that (a) for a square pulse the switch on and off response time is different; (b) for a sinusoidal input the periodicity is preserved; and (c) for a damped sinusoidal the decay rate increases with density and higher energy bands lose energy faster. This work presents a novel methodology for analyzing electron acceleration and transport in flares driven by time-dependent sources.

3I/ATLAS is only the second confirmed cometary object known to enter the Solar System from interstellar space. Cosmogonic considerations suggest that this body may possess relatively high tensile strength and a substantial metal fraction. We present photometric observations along its inbound trajectory toward perihelion, together with a spectroscopic comparison to pristine carbonaceous chondrites from the NASA Antarctic collection. The spectral similarities indicate that 3I/ATLAS may be a primitive carbonaceous object, likely enriched in native metal and undergoing significant aqueous alteration during its approach to the Sun, experiencing cryovolcanism as we could expect for a pristine Trans-Neptunian Object. We propose that the combination of elevated metal abundance and abundant water ice can account for the unusual coma morphology and chemical products reported to date. To do so, corrosion of fine-grained metal grains can originate energetic Fischer-Tropsch reactions, generating specific chemical products into the coma that are not so common in other comets because most of them formed in the outer solar system and didn't inherited so much metal. Interstellar objects such as 3I/ATLAS provide rare opportunities to investigate physical and chemical processes in distant minor bodies of our own Solar System, including trans-Neptunian objects and Oort Cloud comets.

Protoplanetary disks inherit their chemical composition from their natal molecular cloud, but the extent to which this material is preserved versus reset through chemical reprocessing remains an open question. Understanding this balance is a major topic in astrochemistry. Comparing the chemical composition of the envelope and the protoplanetary disk is key to solving the topic. The goal of this paper is to investigate the chemical differences between the disk and the surrounding envelope by comparing the column density ratios of a few selected species in each region. The source we focus on is HL Tau. We present new NOEMA observations of HL Tau targeting the following species: CS, H2CO, H2S, and SO2. We produced zeroth-, first-, and second-moment maps for the species where emission was detected and used them to analyze the spatial distribution and kinematic properties of the different molecules in the disk and the envelope. We derived the column densities and compared the values derived for the envelope and disk. We also computed the rotational diagram for the SO2 detected transitions. Assuming two different temperature regimes, 17 and 58 K, we derived column densities for the species surveyed in the disk and compared them with values derived for the envelope. We find large differences in the derived column density ratios of the surveyed molecules, especially for N(CS)/N(H2S), which is 40 to 50 times larger in the envelope. We attribute these variations to the different excitation and UV-irradiation regimes in the disk and envelope. We also note strong gradients in the ratios between different positions of the disk and tentatively attribute them to different levels of turbulence at different azimuths. The observed differences in molecular ratios in the envelope and the disk are suggestive of chemical reprocessing of the gas during the formation and evolution of the protoplanetary disk.

To efficiently probe primordial non-Gaussianity using Cosmic Microwave Background (CMB) data, we require theoretical predictions that are factorizable, \textit{i.e.}\ those whose kinematic dependence can be separated. This property does not hold for many models, hindering their application to data. In this work, we introduce a general framework for constructing separable approximations to primordial bispectra, enabling direct CMB constraints on arbitrary models including those computed using numerical tools. In contrast to other approaches such as modal decompositions, we learn the basis functions directly from the data, allowing high-fidelity representations with just a handful of terms. This is practically implemented using machine-learning techniques, utilizing neural network basis functions and a loss function designed to mimic the CMB cosine similarity. We validate our pipeline using a variety of input bispectra, demonstrating that the approximations are $>99.5\%$ correlated with the truth with just three terms. By incorporating the neural basis into the \textsc{PolySpec} code, we derive KSW-type CMB estimators, which reproduce local- and equilateral-type non-Gaussianity to within $0.1\sigma$. As a proof-of-concept, we constrain two inflationary bispectra from the `cosmological collider' scenario; these feature an additional strongly-mixed particle sector and cannot be computed analytically. By combining the numerical predictions from \textsc{CosmoFlow} with our factorizable approach (with just three terms), we place novel constraints on the collider models using \textit{Planck} PR4 data, finding no detection of non-Gaussianity. Our method facilitates detailed studies of the inflationary paradigm, connecting modern theoretical tools with high-resolution observational data.

We present the commissioning results and first scientific observations from the newly installed TWIST observatory - a 50 cm telescope equipped with an sCMOS camera providing a $36.1\times24.1$ arcmin$^2$ field of view - housed in the former SuperWASP-North enclosure. We conducted a 67-night, 199-day baseline white-light monitoring campaign centred on the globular cluster M3 aimed at characterizing stellar variability within the cluster while also assessing the photometric performance of the newly commissioned system. We report the discovery of four new SX Phoenicis variables (V301-304), confirm their cluster membership, and identify fundamental-mode pulsation in one, allowing an independent period-luminosity-based distance estimate to M3. We revisited 231 previously known RR Lyrae stars, providing updated period measurements for 203 and white-light amplitudes for 198. We detected Blazhko-like modulation in 53 stars and characterized the modulation parameters for 28. Notably, we measure periods and amplitudes for the unclassified variables V286 and V287 for the first time. We also identify three foreground flaring M dwarfs, and assess the feasibility of detecting microlensing events in M3, concluding that expected rates are negligible. Alongside the scientific results, we introduce a new correction technique for flat-field images affected by scattered light and present a full characterization of the observatory's photometric capabilities. These results demonstrate the scientific utility of TWIST for high-cadence time-domain surveys using modest-aperture instrumentation.

$E_G$ is a summary statistic that combines cosmological observables to achieve a test of gravity that is relatively model-independent. Here, we consider the power of a measurement of $E_G$ using galaxy-galaxy lensing and galaxy clustering with sources from the Rubin Observatory's Legacy Survey of Space and Time (LSST), and lenses from the Dark Energy Spectroscopic Instrument (DESI). We first update the theoretical framework for the covariance of $E_G$ to accommodate this Stage IV scenario. We then demonstrate that $E_G$ offers in principle a model-agnostic test of gravity using only linear-scale information, with the caveat that a careful treatment of galaxy bias is required. We finally address the persistent issue of $E_G$'s theoretical dependence on the measured value of $\Omega_{\rm M}^0$. We propose a framework that takes advantage of the posterior predictive test to consistently incorporate uncertainty on $\Omega_{\rm M}^0$ in tests of gravity with $E_G$, which should be of general use beyond the LSST+DESI scenario. Our forecasting study using this method shows that the prior information available for $\Omega_{\rm M}^0$ is instrumental in determining the power of $E_G$ in the LSST+DESI context. For the full survey dataset, with priors on $\Omega_{\rm M}^0$ from existing CMB data, we find that for some modified gravity scenarios considered, we are likely to be able to reject the GR null hypothesis.

We present the scientific case for exploiting the capabilities of the PLATO mission to study bright pulsating white dwarfs across a wide spectral range, including hydrogen-deficient types (GW Vir and DBV stars) and hydrogen-rich classes (classical DAVs, pulsating extremely low-mass DA white dwarfs, and ultra-massive DA white dwarfs). PLATOs exceptional photometric precision, long-duration continuous monitoring, and extensive sky coverage promise transformative advances in white dwarf asteroseismology. Our key objectives include probing the internal structure and chemical stratification of white dwarfs, detecting secular changes in pulsation modes over extended timescales, and discovering rare or previously unknown classes of pulsators. To assess feasibility, we constructed a sample of 650 white dwarf candidates identified within PLATOs Southern LOPS2 field using the PLATO complementary science catalogue combined with Gaia DR3, and derived atmospheric parameters through photometric modeling. This sample comprises 118 DA white dwarfs (including 23 ZZ Ceti candidates), and 41 non-DAs (including 35 DBV candidates). Simulated observations using PlatoSim demonstrate that PLATO will be capable of detecting white dwarf pulsation modes with amplitudes as low as 0.1 mma depending on stellar magnitude, observation duration, pixel location, and the number of contributing cameras. We provide detailed detection limits and visibility forecasts for known pulsators across a representative range of these parameters. Furthermore, we emphasize strong synergies with Gaia astrometry, TESS photometry, and targeted spectroscopic campaigns, which together will enable robust mode identification and detailed stellar modeling. Collectively, these efforts will unlock unprecedented insights into white dwarf origins, evolution and internal physics, and the fate of their planetary systems.

We present a data-driven method for determining the inclination and position angle (PA) of disc galaxies using a Fourier-Laguerre basis decomposition of imaging data. We define a dimensionless metric, $\eta$, that characterises the ratio of the quadrupole and monopole coefficients in the Fourier-Laguerre basis function expansion. This metric serves as a robust measure which is related to the inclination of a galaxy. We find an empirical relationship between $\eta$ and inclination which is agnostic to the galaxy morphology. The PA is derived directly from the phase of the quadrupolar Fourier-Laguerre functions. Across a benchmark sample of galaxies, the method reproduces published inclination and PA values to within a median of 10$^\circ$ and 5$^\circ$, respectively, while also demonstrating essentially zero catastrophic failures. Applying this pipeline to galaxies from the Siena Galaxy Atlas (SGA), we report measurements of $\eta$, scale length and PA for three different bands of 133,942 disc galaxies. Our computationally inexpensive technique automates parametrisation analysis and returns reproducible results for large surveys. We release a Python package ready for application to next generation surveys.

Sunspots depict large variability during the last 100 years, a period called the Modern Maximum (MM). However, other variables are needed to study the long-term evolution, e.g., of weaker fields and different radiative emissions. Recently, the relation between sunspots and the F10.7 and F30 radio fluxes and the MgII index (proxies of EUV irradiance) was found to vary during the last 70 years so that a relative sunspot dominance over EUV in the 1950s-1960s changed to EUV dominance in the 2000s (Mursula et al., 2024). Here we use data from eight long-operating observatories to calculate the yearly range of daily variation of the geomagnetic Y-component, the rY index, for the last 137 years. The rY index correlates very well with the MgII index and the solar F30 radio flux. These three indices have no trend relative to each other. On the other hand, the F10.7 flux has a significant trend with respect to the three co-varying EUV indices (MgII, F30, rY). Therefore, the rY index replaces F10.7 as the best long-term EUV proxy, and extends the MgII index by 90 years. We verify that all the four EUV proxies (rY, MgII, F30, F10.7) have an increasing trend with respect to sunspots during the last 50-70 years. This is valid both for sunspot numbers and group numbers. We find that the relation between rY index and sunspots has a quadratic evolution over the MM. The Sun has more sunspots relative to EUV irradiance during the growth and maximum of the MM, while the opposite is true during its decay. We estimate that the MgII index increases by 24\% of its solar cycle variation with respect to the sunspot number during the last 70 years. Our results indicate a systematic difference in the evolution between sunspots (photosphere) and plages (chromosphere) with long-term solar activity. The implied varying spot-facula ratio has consequences to the stellar evolution of the Sun and Sun-like stars.

Luminous Red Novae (LRNe) have been argued to be related to the ejection of common envelopes (CEs) in binary star systems. Ejection of CEs leads to tightened stellar orbits capable of forming compact binaries that merge in Hubble time. As these mergers are seen by gravitational-wave (GW) detectors such as LIGO, Virgo and KAGRA (LVK), we ask what the merger rates of compact binaries in LVK tell us about the fraction of LRNe that lead to the formation of compact binaries that merge in Hubble time. Using the observed volumetric rates of LRNe from the Zwicky Transient Facility (ZTF) and of compact binary mergers from LVK observations, we derive limits on the fraction of LRNe that produce compact binaries that merge in Hubble time. Assuming the LRNe rate closely follows the star formation rate at any redshift, we use the delay time distribution models for compact binaries to compute the compact binary merger rate. A comparison of this merger rate with the latest volumetric rates of compact binary mergers from the fourth GW transient catalog (GWTC-4) at the present epoch of LVK allows us to constrain the above fraction. We find that only a fraction as small as $\sim 10^{-3}$ (median) of the LRNe correspond to the GW-observed binary neutron star (BNS) and neutron star-black hole (NSBH) mergers. This potentially implies that the majority of the LRNe population will not lead to mergers of compact objects, but other end products, such as stellar mergers.

Cosmological N-body simulations have long suggested that the pseudo phase-space density (PPSD), $\rho/\sigma^3$, of cold dark matter halos follows the universal relation $\rho/\sigma^3 \propto r^{\chi}$, with $\chi \approx -1.875$, as predicted by spherical secondary-infall similarity solutions. This power law appears to hold despite the fact that neither the density $\rho(r)$ nor velocity dispersion $\sigma(r)$ follow universal power law relations individually, even at fixed mass. We analyze 246 host halos from the Symphony suite of high-resolution cosmological zoom-in simulations, to consistently measure PPSD profiles across host masses from $10^{11}$ to $10^{15} M_\odot$. We find that the PPSD systematically deviates from a power law, and that halos with larger deviations from Jeans equilibrium systematically develop steeper average PPSD slopes. This result suggests that the PPSD is not universal; instead, it is linked to a halo's degree of dynamical equilibrium, which is ultimately set by halo formation history. As a result, we show that secondary halo properties such as concentration and accretion rate inherit significant correlations with the PPSD slope. Moreover, our hosts' PPSD profiles are remarkably consistent with predictions from 1D self-similar fluid collapse models, indicating that three-dimensional structure, velocity anisotropy, and cosmological environment all play negligible roles in shaping the PPSD. These findings imply that the PPSD does not follow a universal power law, but is instead determined by halo mass assembly history alone.

Stellar variability and transient events provide critical insights into astrophysics, accelerated by missions like CoRoT, Kepler, and K2. NASA's Transiting Exoplanet Survey Satellite (TESS) adds a unique combination of long baseline and all-sky coverage, though extracting light curves from full-frame images (FFIs) is challenging due to scattered light and blending. We processed TESS FFIs to produce TEQUILA (TESS quick-look and light curve analysis), a comprehensive catalog of variable point sources from the prime mission, enabling diverse studies without requiring raw data processing. We used difference image analysis, constructing reference images from quality-filtered FFIs for each CCD across sectors 1-26. Iterative subtraction mitigated systematics, and light curves were created using aperture photometry for sources varying in residual images. The pipeline yields over six million light curves, including stellar variables, transients, systematics, and moving objects. Approximately $6 \times 10^5$ span multiple sectors, with roughly $10^3$ from continuous viewing zones. We achieve median differential variability noise of $10^{-3}$ to $10^0$ for sources between 5.0 and 16.0 Tmag, while typical photometric RMS variability ranges from $10^{-2}$ to $10^1$. A convolutional neural network identifies light curves caused by instrumental noise, assigning a confidence score to each classification. To avoid confusion with astrophysical variables, we also flag light curves prompted by known Solar System objects (SSOs). All light curves are accessible via MAST as a High-Level Science Product. This catalog serves as a discovery tool for new variables; future work will refine methods and extend coverage to the TESS extended mission.

Model-independent approaches have gained increasing attention as powerful tools to investigate persistent tensions between cosmological observations and the predictions of $\Lambda$CDM. Notably, recent DESY5 Type Ia Supernovae (SNIa) and DESI Baryon Acoustic Oscillation (BAO) data challenge the validity of the cosmological constant, and they remain in tension with SH0ES local distance ladder measurements under standard pre-recombination physics. Building on our previous work, MNRAS 523 (2023) 3, 3406-3422, we present a follow-up analysis of the model-independent calibration of the local and inverse distance ladders using cosmic chronometers (CCH) data and Gaussian Processes. We jointly constrain the SNIa absolute magnitude, $M$, the comoving sound horizon at the baryon-drag epoch, $r_d$, and the spatial curvature parameter, $\Omega_k$, using CCH with DESY5 and DESI DR1/DR2. We find this data combination compatible with a flat universe at $\sim1.7\sigma$, with $\Omega_k=-0.143\pm0.085$, showing weaker compatibility than with Pantheon+, while the ladder calibrators read $M=-19.324_{-0.095}^{+0.092}$ and $r_d=(144.00^{+5.38}_{-4.88}$) Mpc. Although current uncertainties limit the precision of our constraints and prevent us from arbitrating the Hubble tension, it is nevertheless instructive to explore the constraining power of our methodology with future SNIa, CCH, and BAO from surveys such as LSST, Euclid, and DESI. We present the first forecast analysis for the triad $(M,\Omega_k,r_d)$, finding that, in an optimistic scenario, upcoming data will improve agnostic constraints on $M$ by $\sim$54% and on $r_d$ by $\sim$66%, enabling a $\sim2$% determination of $H_0$. Precision on $\Omega_k$ will increase by $\sim50$%. Our analysis outlines which improvements in future data - whether in quality, quantity, or redshift coverage - are likely to most effectively tighten these constraints.[abridged]

Statistically-significant differences in the value of the Hubble parameter are found depending on the measurement method that is used, a result known as the Hubble tension. A variety of ways of comparing, grouping, and excluding measurements have been used to try to explain this, either in terms of physical effects or systematic errors. We present a systematic 'Bayesian jackknife' analysis of 16 independent measurements of the Hubble parameter in an attempt to identify whether the measurements fall into meaningful clusters that would help explain the origin of the tension. After evaluating evidence ratios for the commonly-used split into early- vs late-time measurements, we then study a range of simplified alternative physical scenarios that reflect different physical origins of an apparent bias or shift in the value of $H_0$, assigning phenomenological population parameters to each subset. These include scenarios where specific subsets are biased (e.g. due to unrecognised experimental systematics in the local distance ladder or cosmic microwave background measurements), as well as more cosmologically-motivated cases involving modifications to the expansion history. Many of these scenarios have similar marginal likelihood, but the model where no measurements are biased is strongly disfavoured. Finally, we marginalise over all these scenarios to estimate the 'model agnostic' posterior distribution of $H_0$. The resulting distribution is mildly multi-modal, but modestly favours values near $H_0=68$ km/s/Mpc, with a 95\% credible region of $66.7 < H_0 < 72.7$ km/s/Mpc.

Jets from stellar-mass and supermassive black holes provide the unique opportunity to study similar processes in two very different mass regimes. Historically, the apparent speeds of black hole x-ray binary (BHXRBs) jets have been observed to be lower than jet speeds from active galactic nuclei (AGN) and specifically blazars. In this work, we show that selection effects could be the primary cause of the observed population differences. For the first time, it is possible to perform a statistical analysis of the underlying BHXRB jet Lorentz factor distribution. We use both the Anderson-Darling test and apply nested sampling to this problem. With Bayes factors, we confirm that the Lorentz factor distribution of BHXRBs is best described with a power law, the same model that has been applied to AGN jets. For a Lorentz factor distribution following $\rm N(\Gamma) \propto \Gamma^b$ we find a value for the exponent of $b=-2.64_{-0.55}^{+0.46}$. This exponent is consistent with values found in AGN population studies, within $1\sigma$ for \textit{Swift}-BAT and \textit{Fermi}-LAT selected AGN. The best-fit exponent for the radio selected MOJAVE sample is just above our $2 \sigma$ limit. This is a remarkable agreement given the different scales at which the jets are observed. The observed slower apparent speeds in BHXRBs are largely due to the much larger inclinations in this sample. Furthermore, nested sampling confirms that $\Gamma_{\rm max}$ is completely unconstrained using this method. Therefore, based on kinematics alone, BHXRB jets are broadly consistent with being just as relativistic as those from supermassive black holes.

The historical reconstruction of the Sun's surface magnetic field remains a persistent challenge, limiting our ability to investigate the long-term global properties of the Sun, including the evolution of the large-scale magnetic field, solar cycle prediction, reconstruction of total solar irradiance (TSI), and secular solar variability. In this study, we employ the Advective Flux Transport (AFT) model in conjunction with our newly developed Synthetic Active Region Generator (SARG) to construct a catalog of synthetic active regions (ARs) spanning Solar Cycles 1-24 (1755-2020). We use the SIDC/SILSO sunspot number version 2.0 as the sole input governing the properties of the synthetic ARs in this catalog. This SARG catalog is then incorporated into the AFT model, which simulates the emergence of new ARs on the Sun, which are then transported under the influence of surface flows to produce maps of the full-Sun radial photospheric magnetic field over the entire 265-year period. We modulate the active region tilt for each cycle in order to ensure that the polar fields are consistent with the solar cycle amplitudes. We find that the polar fields derived from these simulations exhibit excellent correlation (r > 0.8) with observational proxies, including polar faculae counts and Ca ii K polar network indices. Daily synchronic maps from these simulations for the entire 265-year period are made publicly available to support a wide range of applications beyond those presented in this work.

Dark matter is approximately five times more abundant than baryonic matter in the universe, but its physical nature continues to elude physicists. One potential candidate for dark matter is a weakly-interacting massive particle (WIMP), which is predicted by various extensions to the Standard Model (SM) of particle physics. After becoming gravitationally bound in cosmic structures, WIMPs can self-annihilate and produce SM particles including neutrinos, which are observable by detectors like IceCube. We present a search for neutrinos from low-mass $(\leq 300 \, \mathrm{GeV})$ WIMP annihilation in dwarf spheroidal galaxies with over seven years of IceCube livetime. We find no statistically significant evidence of neutrinos produced by WIMP annihilation, and therefore set upper limits on the velocity-averaged annihilation cross section $\left<\sigma v\right>$. Our strongest upper limits at the 90\% confidence level are $\mathcal{O}\!\left(10^{-22} \, \mathrm{{cm}^{3} \, s^{-1}}\right)$ for WIMP annihilation directly into neutrino-antineutrino pairs. For our least sensitive channel, the corresponding limits are $\mathcal{O}\!\left(10^{-20} \, \mathrm{{cm}^{3} \, s^{-1}}\right)$, which is an improvement of over two orders of magnitude over previous IceCube limits from dwarf galaxies at the upper end of our mass range.

The nature of the compact objects within the supposed "lower mass gap" remains uncertain. Observations of GW190814 and GW230529 highlight the challenges gravitational waves face in distinguishing neutron stars from black holes. Interpreting these systems is especially difficult because classifications depend simultaneously on measurement noise, compact binary population models, and equation of state (EOS) constraints on the maximum neutron star mass. We analyze 66 confident events from GWTC-3 to quantify how the probability of a component being a neutron star, P(NS), varies across the population. The effects are substantial, the dominant drivers of classification are the pairing preferences of neutron stars with other compact objects, and the neutron star spin distributions. The data reveals that P(NS) varies between 1% - 67% for GW230529's primary and between 51% - 100% for GW190425's primary. By contrast, P(NS) for GW190814's secondary varies by <10%, demonstrating robustness from its high signal-to-noise ratio and small mass ratio. Analysis using EOS information tends to affect P(NS) through the inferred maximum neutron star mass rather than the maximum spin. As it stands, P(NS) remains sensitive to numerous population parameters, limiting its reliability and potentially leading to ambiguous classifications of future GW events.

The observation of multiple images from a strongly lensed gravitational wave (GW) source provides the observer with a stereoscopic view of the source. This allows for a measure of its relative proper motion by comparing the induced GW Doppler shifts between the different images. In addition, if the GW source is in a dynamical environment it will be subject to an acceleration, which will show up as a time dependent Doppler shift in each individual image. In this work we quantify for the first time how a joint detection of these effects can be used to constrain the underlying dynamics and environment of the lensed GW source. We consider a range of different astrophysical environments, from massive clusters to stellar triples, and find that binary black hole (BBH) mergers in Active Galactic Nuclei disks (AGN-disks) are particularly likely to have orbital parameters that can be constrained through our considered lensing setup. Applying these methods to the upcoming catalog of cosmologically strongly lensed GW sources will open up new possibilities for probing their origin and underlying formation mechanisms.

We present the characterization of macrolensing properties of the gravitationally lensed system PKS 1830-211, utilizing data from the Fermi Large Area Telescope. While at gamma-rays we can not spatially resolve the lensed images, a macrolensing-induced time pattern is expected in the blazar's lightcurve, resulting from the delay between variable gamma-ray components originating from its two brightest lensed images. Compared to our previous study, here we employ high-quality lightcurves coupled with prolonged outburst activity, and improved time-series techniques. Analyzing six independent data segments, we identified a delay of 20.26 +/- 0.62 days (statistical and stochastic uncertainty), with a chance detection probability at the 2.5 x $10^{-5}$ level (post-trial). We also present a novel approach to the magnification ratio estimate based on a comparison with simulated data. Our work suggests that the gamma-ray flux ratio between the two main lens components is $\mu_{\gamma} \lesssim$ 1.8. We do not observe convincing evidence of microlensing effects, as previously claimed. The measured gamma-ray time delay is in 2-sigma tension with radio-based estimates, suggesting either distinct emission sites, underestimated radio uncertainties, or gamma-ray production in a region opaque to radio. Our study highlights the potential of well sampled lightcurves and advanced time-series techniques to distinguish true lensing-induced delays from stochastic variability. When combined with improved radio constraints and refined lens models, PKS 1830-211 and similar sources constitute promising systems for time-delay cosmography, offering new insights into both jet structure and cosmological parameters.

We investigate the properties of neutron stars within the framework of $f(Q)$ gravity by incorporating rotational effects through a slowly rotating metric. We derive the modified TOV equations and calculate the angular velocity profiles and moments of inertia (MOI) for linear, quadratic, exponential, and logarithmic $f(Q)$ models. Our results show that deviations in the MOI are more pronounced than those in the stellar mass profiles, suggesting that rotational observables are highly sensitive to geometric corrections. We also calculate a quasi-universal relation between the dimensionless MOI and compactness ($\bar{I}$-$C$). The linear and quadratic models are generally consistent with observational data from PSR J0737-3039A, although the deviations are small and difficult to distinguish from General Relativity due to inherent EoS variability. On other hand, the logarithmic and exponential models show larger deviations (over 20 %), exceeding the EoS-induced uncertainty reported by Suleiman & Read (2024), highlighting the relation's sensitivity to the $f(Q)$ gravity model. These results indicate that $f(Q)$ gravity could potentially be tested in the strong-field regime and point to a direction for future studies, such as investigating EoS-insensitive quasi-universal relations, like the $\bar{I}(\Lambda)$ relations, within the $f(Q)$ framework. Such relations may provide a clearer pathway for exploring possible signatures in strong-field gravity when combined with more precise future observations.

This study presents a bidirectional Long Short-Term Memory (LSTM) neural network for classifying transient astronomical object light curves from the Photometric LSST Astronomical Time-series Classification Challenge (PLAsTiCC) dataset. The original fourteen object classes were reorganized into five generalized categories (S-Like, Fast, Long, Periodic, and Non-Periodic) to address class imbalance. After preprocessing with padding, temporal rescaling, and flux normalization, a bidirectional LSTM network with masking layers was trained and evaluated on a test set of 19,920 objects. The model achieved strong performance for S-Like and Periodic classes, with ROC area under the curve (AUC) values of 0.95 and 0.99, and Precision-Recall AUC values of 0.98 and 0.89, respectively. However, performance was significantly lower for Fast and Long classes (ROC AUC of 0.68 for Long class), and the model exhibited difficulty distinguishing between Periodic and Non-Periodic objects. Evaluation on partial light curve data (5, 10,and 20 days from detection) revealed substantial performance degradation, with increased misclassification toward the S-Like class. These findings indicate that class imbalance and limited temporal information are primary limitations, suggesting that class balancing strategies and preprocessing techniques focusing on detection moments could improve performance.

Private lunar missions are faced with the challenge of robust autonomous navigation while operating under stringent constraints on mass, power, and computational resources. This work proposes a motion-field inversion framework that uses optical flow and rangefinder-based depth estimation as a lightweight CPU-based solution for egomotion estimation during lunar descent. We extend classical optical flow formulations by integrating them with depth modeling strategies tailored to the geometry for lunar/planetary approach, descent, and landing, specifically, planar and spherical terrain approximations parameterized by a laser rangefinder. Motion field inversion is performed through a least-squares framework, using sparse optical flow features extracted via the pyramidal Lucas-Kanade algorithm. We verify our approach using synthetically generated lunar images over the challenging terrain of the lunar south pole, using CPU budgets compatible with small lunar landers. The results demonstrate accurate velocity estimation from approach to landing, with sub-10% error for complex terrain and on the order of 1% for more typical terrain, as well as performances suitable for real-time applications. This framework shows promise for enabling robust, lightweight on-board navigation for small lunar missions.

The IceCube Neutrino Observatory is a cubic-kilometer Cherenkov detector embedded in the Antarctic ice at the South Pole. Its densely instrumented sub-array and dedicated low-energy analyses provide sensitivity to neutrinos in the 5-100 GeV range, enabling precision studies of neutrino oscillations and searches for new physics. This work focuses specifically on this low-energy regime, where sparse hit patterns limit the performance of topology-based reconstruction and classification methods. We introduce Waveform-based Particle Identification (WavePID), a statistically rigorous and interpretable likelihood-ratio discriminator for track-cascade separation, built from Monte Carlo templates in timing-aware, physics-motivated observables and validated through dedicated simulations. Applied to both Monte Carlo and 11.1 years of IceCube data, WavePID suggests improved cascade purity by about 5 percentage points at a fixed 20% down-selection rate relative to the current leading cascade selection, while maintaining Data-MC agreement within detector systematics. The approach is compact and robust to sparse observations, demonstrating the value of waveform-level timing for low-energy reconstruction.

In this paper, we derive the $\mu$-deformed Einstein field equations from the generalized thermodynamic functions of the $\mu$-deformed analog of Bose gas model, applying the (adapted) Verlinde's approach. The basic role of deformation parameter is shown: it provides the possibility to vary the value of the cosmological constant. Due to this, we suggest an interesting treatment of the cosmological constant (CC) problem within the framework of $\mu$-deformation. Namely, viewing the derived $\mu$-deformed CC as an effective one and varying the parameter $\mu$ appropriately, we gain the possibility to drastically reduce the CC, so as to get for it the realistic value. Also, the relation with dark matter is discussed.

Precise alignment sensing and control are essential for maintaining the stability of laser interferometric gravitational-wave detectors. Conventional Wave Front Sensing technique (WFS), which relies on the beat between the carrier and phase-modulated (PM) sidebands, is dominated by arm-axis signals when the carrier resonates in the full interferometer. This dominance limits the detection of other optical axes, such as the Power Recycling Cavity (PRC) and incident beam axes. To address this problem, we propose a novel sensing technique, "Phase-Modulated-sideband $\times$ Phase-Modulated-sideband Wave Front Sensing" (PMPMWFS), which demodulates the beat signal at the difference frequency between two anti-resonant PM sidebands. We derived the theoretical response of PMPMWFS and experimentally demonstrated it using the Power-Recycled X-arm (PRXARM) configuration of KAGRA. The results show that PMPMWFS effectively decouples angular fluctuation signals of the PRC and incident beam from those of the arm cavity and provides orthogonal signal components for the end mirror of the arm cavity. Furthermore, feedback control using PMPMWFS achieved stable interferometer locking for over one hour. These results demonstrate that PMPMWFS offers an effective sensing method for decoupling multiple alignment degrees of freedom in future gravitational-wave detectors.

In this work we examine the internal structure of compact stars within an extended gravitational framework described by the function $f(\mathcal{R},\mathcal{G},\mathcal{T})$. Throughout this work, the quantity $\mathcal{R}$ refers to the curvature scalar formed from the Ricci tensor. The term $\mathcal{G}$ denotes the Gauss--Bonnet curvature invariant, while $\mathcal{T}$ corresponds to the trace obtained by contracting the matter energy-momentum tensor. Our analysis is directed toward massive radio pulsars with masses above $1.8\,M_{\odot}$, which provide an exceptional testing ground for gravity under conditions inaccessible to laboratory experiments. Adopting the linear form $f(\mathcal{R},\mathcal{G},\mathcal{T})=\mathcal{R}+\alpha\,\mathcal{G}+\beta\,\mathcal{T}$ where $\alpha$ and $\beta$ are parameters of suitable dimensionality,\footnote{$\alpha$ has dimensions of $[L^{2}]$ and $\beta$ carries units of $[N^{-1}]$.} we obtain an exact analytic solution for static anisotropic stellar matter in hydrostatic equilibrium. This solution allows all physical quantities to be expressed in terms of the dimensionless parameters $ \alpha_{1}=\alpha/R^{2},\qquad \beta_{1}=\beta/\kappa^{2}$ together with the compactness $C=2GM/(Rc^{2})$. We constraint the two parameters $\alpha$ and $\beta$ by matching the model with the mass and radius of pulsar \textit{U1724} requires restricting these parameters to $\alpha_{1}=\pm0.023$ and $\beta_{1}=\pm0.001$, where $\kappa^{2}=8\pi G/c^{4}$ is the standard Einstein coupling. The resulting stellar configuration satisfies the causal bound on the radial sound speed, $c_{s}^{2}<c^{2}/3$, distinguishing it from the corresponding behaviour in general relativity.

We propose and analyze an autocatalytic mechanism in which bursts of Hawking radiation from evaporating micro-primordial black holes (PBHs) trigger the collapse of near-critical plasma overdensities. In a primordial plasma seeded with such patches, this feedback self organizes into a traveling ignition front that successively forms new PBHs and then self-quenches as the Universe expands. A minimal reaction-diffusion model yields conservative criteria for ignition and freeze-out and predicts a stochastic gravitational-wave background with a sharp causal low-frequency edge set by the freeze-out correlation length and largely insensitive to Planck-scale PBH endpoint microphysics. The resulting sub-Hz-to-audio band and amplitudes satisfy cosmological energy-injection bounds, providing a clean, testable target for forthcoming gravitational-wave observatories.

In this paper we discuss how the Cosmic Neutrino Background can affect the measured phase difference in a matter interferometer. This phase is proportional to a difference in potential energies along the two interferometer paths. The relevant potentials here are the well-known neutrino matter potential and a potential related to the Stodolsky effect. We show how they can be rewritten in terms of scalar potentials, pseudo magnetic fields and spin-spin interactions. Unfortunately, current technology is unlikely to detect this effect and we discuss prospects for the future. We also briefly comment on fermionic Dark Matter which can give rise to very similar effects which can easily be larger than the neutrino case.

We study the influence of a massive vector (Proca) field on the energy fluxes from extreme-mass-ratio inspirals (EMRIs) around a slowly rotating Kerr black hole. The secondary compact object, carrying a Proca hair, emits additional dipolar radiation that alters total energy flux relative to general relativity (GR). These modifications induce a secular drift in the orbital evolution of circular geodesic orbits, leading to measurable dephasing in the resulting EMRIs waveforms. By evaluating waveform mismatches between the Einstein-Proca framework and its GR counterpart, we show that the Laser Interferometer Space Antenna (LISA) can distinguish the signatures of a light Proca field when black hole rotation is included. Furthermore, using a Fisher information matrix analysis, we forecast LISA's capability to place stringent constraints on the Proca mass with EMRIs signal from slowly rotating Kerr black holes. For representative EMRIs configurations, we find that LISA can detect or constrain Proca masses down to $\mu_v\sim 10^{-20}$eV, with typical fractional uncertainties at the level of tens percent, depending on the black-hole spin.

Geostationary hyperspectral satellites generate terabytes of data daily, creating critical challenges for storage, transmission, and distribution to the scientific community. We present a variational autoencoder (VAE) approach that achieves x514 compression of NASA's TEMPO satellite hyperspectral observations (1028 channels, 290-490nm) with reconstruction errors 1-2 orders of magnitude below the signal across all wavelengths. This dramatic data volume reduction enables efficient archival and sharing of satellite observations while preserving spectral fidelity. Beyond compression, we investigate to what extent atmospheric information is retained in the compressed latent space by training linear and nonlinear probes to extract Level-2 products (NO2, O3, HCHO, cloud fraction). Cloud fraction and total ozone achieve strong extraction performance (R^2 = 0.93 and 0.81 respectively), though these represent relatively straightforward retrievals given their distinct spectral signatures. In contrast, tropospheric trace gases pose genuine challenges for extraction (NO2 R^2 = 0.20, HCHO R^2 = 0.51) reflecting their weaker signals and complex atmospheric interactions. Critically, we find the VAE encodes atmospheric information in a semi-linear manner - nonlinear probes substantially outperform linear ones - and that explicit latent supervision during training provides minimal improvement, revealing fundamental encoding challenges for certain products. This work demonstrates that neural compression can dramatically reduce hyperspectral data volumes while preserving key atmospheric signals, addressing a critical bottleneck for next-generation Earth observation systems. Code - this https URL

Gravitational wave detection has transformed astrophysics, granting us direct access to black hole mergers, neutron star collisions, and the cataclysms of stellar death. Yet the great observatories of today, LIGO, Virgo, KAGRA, and the planned Einstein Telescope, rest on Michelson interferometers that, despite their triumphs, confront fundamental barriers of scale, cost, and environmental vulnerability. We envision a new path, a Sagnac-based fiber interferometer that leverages reciprocity and inherent robustness. Its meter-scale, modular design compact enough to fit within a small facility, offers dramatic gains in scalability and affordability over kilometer scale Michelson systems. Tunable to frequency bands where conventional detectors lose sensitivity, it opens the door to compact, versatile, and accessible GW observatories, empowering universities and research centers worldwide. Linked together in a global network, such facilities could transcend mere detection, they could localize cosmic sources and reconstruct them into images, much as black holes were first directly revealed, ushering in a new era of gravitational-wave astronomy and multi-messenger discovery.

We explore a novel cosmological model based on coupled fields in the framework of scalar tensor theories, considering the specific interplay between gravity and scalar fields. The model further extends a recent axion-dilaton system by introducing viable couplings with the space--time geometry encoded into the scalar curvature. After briefly introducing the action and the corresponding field equations, we employ linear stability theory to investigate the physical properties. The analysis showed the compatibility of the current theoretical model with the recent history of the Universe, obtaining viable constraints for the model's parameters in some specific cases. In the present setup, the axion--dilaton system is non--minimally coupled with gravity in an independent manner, leading to distinct physical features in the phase-space structure, possible alleviating the cosmic coincidence problem.

The current KM3NeT/ORCA neutrino telescope, still under construction, has not yet reached its full potential in neutrino reconstruction capability. When training any deep learning model, no explicit information about the physics or the detector is provided, thus they remain unknown to the model. This study leverages the strengths of transformers by incorporating attention masks inspired by the physics and detector design, making the model understand both the telescope design and the neutrino physics measured on it. The study also shows the efficacy of transformers on retaining valuable information between detectors when doing fine-tuning from one configurations to another.

In shift-symmetric Einstein-scalar-Gauss-Bonnet gravity, stationary black holes have a non-vanishing scalar charge. During the inspiral, the phase evolution is modified by several effects, primarily an additional scalar dipole radiation, which enters at -1PN order. This effect accelerates the inspiral when compared to general relativity, when including corrections up to 2PN. Using fully non-linear numerical simulations of quasi-circular, comparable mass binaries, we find that in the late stages the orbital dynamics are altered so that the overall effect is instead a decelerated merger phase for the modified gravity case. We attribute this to a change in the conservative dynamics, and show that at the late inspiral stage more energy must be emitted in scalar-Gauss-Bonnet gravity to induce a given change in frequency. In longer signals, this should lead to a distinctive switch between a faster and slower frequency evolution relative to general relativity as the binary approaches merger. This work suggests we should revisit existing constraints on the theory that are obtained assuming PN approximations apply up to merger, or based on order by order approximations that neglect backreaction effects on the metric, and shows the importance of including non-linear effects that modify the gravitational sector in the strong field regime.

We revisit the framework of Newer General Relativity, defined by all independent quadratic invariants of the non-metricity tensor, including the unique quadratic parity-violating term. We analyze linear perturbations around a flat FLRW background and find that the theory generically exhibits ghost instabilities, except for the special case of Symmetric Teleparallel Equivalent of General Relativity (STEGR) supplemented by the parity-violating operator. Using both the Dirac-Bergmann procedure and the Cartan-Kuranishi algorithm, we show that this remaining case propagates eight degrees of freedom at the fully non-linear level. This result implies that the parity-violating extension of STEGR is strongly coupled around FLRW and therefore lies beyond the regime of validity of standard linear perturbation theory.

We analyze the impact of breaking diffeomorphism invariance in the inflaton sector. In particular, we consider inflaton models which are invariant under the subgroup of transverse diffeomorphisms and address the possibility of implementing a slow-roll phase. We obtain the corresponding expressions for relevant quantities such as the slow-roll parameters and the number of $e$-folds, and derive the primordial power-spectrum of curvature perturbations. The scalar spectral index features modifications which are confronted with CMB data from Planck and ACT. We study in detail the quadratic potential model, combining asymptotic and numerical analysis. We show that the post-inflationary behavior can be drastically different from the diffeomorphism-invariant case, exhibiting novel dynamical regimes.

Physics simulators are essential in science and engineering, enabling the analysis, control, and design of complex systems. In experimental sciences, they are increasingly used to automate experimental design, often via combinatorial search and optimization. However, as the setups grow more complex, the computational cost of traditional, CPU-based simulators becomes a major limitation. Here, we show how neural surrogate models can significantly reduce reliance on such slow simulators while preserving accuracy. Taking the design of interferometric gravitational wave detectors as a representative example, we train a neural network to surrogate the gravitational wave physics simulator Finesse, which was developed by the LIGO community. Despite that small changes in physical parameters can change the output by orders of magnitudes, the model rapidly predicts the quality and feasibility of candidate designs, allowing an efficient exploration of large design spaces. Our algorithm loops between training the surrogate, inverse designing new experiments, and verifying their properties with the slow simulator for further training. Assisted by auto-differentiation and GPU parallelism, our method proposes high-quality experiments much faster than direct optimization. Solutions that our algorithm finds within hours outperform designs that take five days for the optimizer to reach. Though shown in the context of gravitational wave detectors, our framework is broadly applicable to other domains where simulator bottlenecks hinder optimization and discovery.

We investigate an one loop radiative neutrino mass model, where the loop particles, notably a singlet fermion ($\chi$), a doublet fermion ($\Psi$) and three generations of singlet scalars ($\phi_i, i=\{1,2,3\}$) are assumed to be odd under an additional $\mathcal{Z}_2$-symmetry. In this setup, the singlet fermion mixes with the neutral component of the doublet to give rise singlet-doublet Majorana dark matter. The addition of $\mathcal{Z}_2$ odd scalars in the model provides rich phenomenological implications. We find that the quartic interaction terms between the SM Higgs and $\phi_i$s play a significant role in modifying the scalar potential to have a first-order phase transition (FOPT) leading to observable gravitational waves (GWs) spectra. We also examine the non-trivial role played by the singlet-doublet fermion DM and the scalars in loop-induced neutrino mass, $(g-2)_\mu$, and lepton flavor violation. We find that the model is predictive due to the combined constraints and can be verified at different terrestrial experiments.

Conditional diffusion models provide a natural framework for probabilistic prediction of dynamical systems and have been successfully applied to fluid dynamics and weather prediction. However, in many settings, the available information at a given time represents only a small fraction of what is needed to predict future states, either due to measurement uncertainty or because only a small fraction of the state can be observed. This is true for example in solar physics, where we can observe the Sun's surface and atmosphere, but its evolution is driven by internal processes for which we lack direct measurements. In this paper, we tackle the probabilistic prediction of partially observable, long-memory dynamical systems, with applications to solar dynamics and the evolution of active regions. We show that standard inference schemes, such as autoregressive rollouts, fail to capture long-range dependencies in the data, largely because they do not integrate past information effectively. To overcome this, we propose a multiscale inference scheme for diffusion models, tailored to physical processes. Our method generates trajectories that are temporally fine-grained near the present and coarser as we move farther away, which enables capturing long-range temporal dependencies without increasing computational cost. When integrated into a diffusion model, we show that our inference scheme significantly reduces the bias of the predicted distributions and improves rollout stability.

Neutrinos can acquire "refractive masses" as a consequence of their interactions with ultralight dark matter (DM). We explore a model with two additional sterile neutrinos and an ultralight scalar field which acts as DM and interacts with all five neutrinos. We show that the effective $5 \times 5$ Hamiltonian for neutrino propagation can be diagonalized by a unitary matrix $\mathbb{P}$ parametrized by 6 mixing angles and 1 complex phase. When active-sterile mixing angles are small, we identify a parametrization for $\mathbb{P}$ that reduces neutrino propagation inside the Sun to a two-flavor problem for a uniform DM background. In the presence of a DM halo inside the Sun, however, the propagation shows additional features in the region of halo dominance. We derive approximate analytic expressions for the electron neutrino survival probability in the presence of the DM halo. We show that this probability has a strong dependence on the neutrino production region even for a fixed energy, and numerically calculate the effects of averaging over these production regions. Comparisons with the re-interpreted solar data, in the light of possible active-sterile neutrino conversions, would allow putting bounds on the halo parameters. Finally, we examine the possibility of reviving the dark-LMA solution in this context, where the survival probability spectrum can have attractive features aligned with the measurements at Super-Kamiokande.

It is known that Yang-Mills theories, especially in the large-$N$ limit, exhibit a $\theta$-vacuum structure with a multi-branched vacuum energy. In this work, we demonstrate that this multi-branch structure can play a crucial role in axion cosmology when the axion acquires its mass from the Yang-Mills sector, even when that sector is never reheated by the inflaton. The key observation is that the axion potential is directly tied to the tunneling rate between adjacent branches. We find qualitatively new phenomena, including a new class of first-order phase transitions, bouncing bubbles, and nested "bubbles-within-bubbles." When the axion has a decay constant around the Planck scale, as motivated by the string Axiverse, the axion can be driven near the hilltop by inflationary dynamics, allowing the phase transition to be triggered. The associated energy release can be large enough to generate a significant stochastic gravitational-wave background, produce primordial black holes, or populate the Yang-Mills sector with particles. These phenomena represent novel predictions of the Axiverse and should be taken into account when assessing the cosmological impact of axions or axion-like particles. To recover conventional axion cosmology, one must suppress or avoid the dynamics discussed in this work.